Skip to main content
Log in

Thermodynamic modeling of high-grade metabasites: a case study using the Tso Morari UHP eclogite

  • Original Paper
  • Published:
Contributions to Mineralogy and Petrology Aims and scope Submit manuscript

Abstract

Thermodynamic modeling is an important technique to simulate the evolution of metamorphic rocks, particularly the poorly preserved prograde metamorphic reactions. The development of new thermodynamic modeling techniques and availability of updated thermodynamic databases and activity–composition (a–X) relations, call for an evaluation of best practices for modeling pressure–temperature (P–T) paths of metabasites. In this paper, eclogite from the Tso Morari UHP terrane, NW India, is used as a representative metabasite to directly compare the outputs (pseudosections and P–T paths) generated from recent versions of the widely used THERMOCALC and Theriak-Domino programs. We also evaluate the impact of using the most updated thermodynamic database (ds 62, Holland and Powell in J Metamorph Geol 29(3):333–383, 10.1111/j.1525-1314.2010.00923.x, 2011) relative to an older version (ds 55, Holland and Powell in J Metamorph Geol 16(3):309–343, 10.1111/j.1525-1314.1998.00140.x, 1998), and the effect of the user’s choice of mineral a–X relations while considering the effect of garnet fractionation on the rock’s effective bulk composition. The following modeling protocols were assessed: (1) TC33; THERMOCALC version 3.33 with database ds 55 and garnet a–X relations of White et al. (J Metamorph Geol 25(5):511–527, 10.1111/j.1525-1314.2007.00711.x, 2007); (2) TC47; THERMOCALC version 3.47 with database ds 62 and garnet a–X relations of White et al. (J Metamorph Geol 32(3):261–286, 10.1111/jmg.12071, 2014a); (3) TDG; Theriak-Domino with database ds 62 and garnet a–X relations of White et al. (2014a), and (4) TDW; Theriak-Domino with database ds 62 and garnet a–X relations of White et al. (2007). TC47 and TDG modeling yield a similar peak metamorphic P–T of 34 ± 1.5 kbar at 544 ± 15 °C and 551 ± 12 °C, respectively. The results are 5–8 kbar higher in pressure than that determined from TC33 modeling (26 ± 1 kbar at 565 ± 8 °C), and TDW modeling (28.5 ± 1.5 kbar at 563 ± 13 °C). Results indicate that all four modeling protocols generally provide consistent metamorphic phase relations and thermodynamic simulations regarding fractionation of the bulk composition and prograde metamorphism within uncertainty. In all model calculations, the initial bulk composition measured by XRF does not represent the effective bulk composition at the time of garnet nucleation. The choice of garnet a–X relations can affect predictions of peak pressure, regardless of program choice. This study illustrates the importance of careful consideration of which a–X relations one chooses, as well as the need for comparison between modeling predictions and evidence from the geochemistry and petrography of the rock(s) themselves.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Fig. 1
Fig. 2
Fig. 3
Fig. 4
Fig. 5
Fig. 6
Fig. 7
Fig. 8
Fig. 9
Fig. 10
Fig. 11
Fig. 12
Fig. 13
Fig. 14
Fig. 15

Similar content being viewed by others

References

Download references

Acknowledgements

We thank Emilee Darling for preparing the samples. C. Menold thanks Dennis Donaldson and Alex Webb for help in the field. The manuscript was greatly improved with the help of detailed reviews from two anonymous reviewers, to whom we are grateful. We also thank the CMP editorial team for their assistance in the review process. This research was funded by NSF EAR -1822524 to C. Macris and C. Menold.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Catherine A. Macris.

Additional information

Communicated by Mark S Ghiorso.

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Electronic supplementary material

Below is the link to the electronic supplementary material.

410_2020_1717_MOESM1_ESM.tif

Supplementary file1 Figure S1. Pseudosections and garnet phase boundaries constructed by (a) TC33, (b) TC47, (c) TDG, and (d) TDW, using incrementally fractionated bulk composition. 14 steps of garnet removal and EBC calculations have been performed. (TIF 20196 kb)

Supplementary file2 (DOCX 22 kb)

Supplementary file3 (DOCX 14 kb)

Supplementary file4 (DOCX 29 kb)

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Pan, R., Macris, C.A. & Menold, C.A. Thermodynamic modeling of high-grade metabasites: a case study using the Tso Morari UHP eclogite. Contrib Mineral Petrol 175, 78 (2020). https://doi.org/10.1007/s00410-020-01717-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1007/s00410-020-01717-w

Keywords

Navigation