skip to main content

Title: Thermodynamic Modeling of the Tso Morari UHP Eclogite, NW Himalaya
Thermodynamic modeling is an important technique to interpret metamorphic phase relations and calculate model pressure-temperature (P-T) paths for metamorphic rocks. This study uses representative, coesite-bearing eclogites from the Tso Morari UHP terrane of the NW Himalaya to simulate its prograde metamorphism using multiple modeling programs and thermobarometry. Our modeling yields a peak metamorphism P-T of ~32-33 kbar and ~560-570 °C by the THERMOCALC345 and Theriak-Domino programs (Green et al., 2016), which is ~5 kbar higher in pressure and ~15 °C lower in temperature than that determined by using THERMOCALC333 (White et al., 2007) (~27.8 kbar and ~580 °C). The significantly higher pressure obtained using the THERMOCALC345 and Theriak-Domino is likely a result of the upgrade of thermodynamic parameters of minerals (i.e. garnet Wpy-gr and agr) in the newer a-x relations. The modeled effective bulk compositions and mineral stabilities along the calculated P-T path show different patterns under the two modeling techniques. Modeling by the Theriak-Domino programs is preferred in this case because the results are more consistent with the measured mineral compositions of our rocks. Multiple thermobarometers by garnet-omphacite-phengite, garnet-omphacite, garnet-phengite on the garnet rim, high-Si phengite and matrix omphacite yield a peak metamorphism of ~ 28.5-29.0 kbar and ~ 650-728 more » °C, which is generally consistent with the modeled P-T path. Based on our model calculations, the initial bulk composition measured by XRF does not represent the reactant bulk composition at the time of garnet nucleation, and this compositional discrepancy possibly is caused by the crystallization of pre-garnet minerals (i.e. hematite), reaction overstepping, or partial reequilibration. In summary, by implementing and evaluating multiple modeling strategies and considering the petrography and metamorphic mineralogy of the rocks, this study finds that the eclogite modeling using Theriak-Domino programs in the Tso Morari terrane provide more consistent metamorphic phase relations and more reasonable thermodynamic simulations regarding fractionation of the bulk composition and prograde metamorphism. References: Green et al. J Metamorph Geol 34, 845-869 (2016) White et al. J Metamorph Geol 25, 511-527 (2007) « less
; ;
Award ID(s):
Publication Date:
Journal Name:
AGU Fall Meeting 2020
Sponsoring Org:
National Science Foundation
More Like this
  1. The Tso Morari terrane within the Himalayan orogeny underwent ultrahigh-pressure (UHP) metamorphism due to northward subduction under the Eurasian continent during the early Eocene. The advancement of computational petrology and availability of relevant thermodynamic databases provide the mechanism to more precisely quantify metamorphic processes. In this study, we model the eclogite’s prograde pressure-temperature (P-T) path as well as multiple fluid infiltration events during exhumation using Theriak-Domino with dataset ds62 and garnet[1] and other metabasic mineral activity-composition relations. The effect of garnet fractionation on the rock’s effective bulk composition is considered in simulating prograde garnet growth. A “fishhook” shape clockwise P-Tmore »path is obtained with a peak pressure of ~28.5 kbar at ~563 °C, followed by a peak temperature of ~613 °C at ~24.5 kbar[2]. Thermodynamic modelling using P-M(H2O) pseudosections on Tso Morari eclogites indicates three distinct phases of fluid infiltration during exhumation. Fluid infiltration Ⅰ occurs at ~610 °C and ~23.5 kbar with ~3.1 mol % fluid expulsion due to the destabilization of lawsonite. The modelling results are consistent with petrographic observations in the eclogite: we found ~6.0 vol % epidote and ~21.0 vol % amphibole and the possible pre-existence of lawsonite evidenced by its pseudomorph (as epidote and paragonite aggregates) in a garnet core and rim[3], and CNASH modelling on the epidote and its inclusion paragonite. Fluid infiltration Ⅱ occurs at ~9.2 kbar and ~608 °C with >2.6 mol % fluid infiltration at amphibolite-facies. This phase of fluid infiltration is characterized by aggressive amphibolization from the boudin core to rim. Fluid infiltration Ⅲ occurs at ~610 °C and ~8.7 kbar, caused by breakdown of phengite as predicted through modelling the symplectitic association (plagioclase, biotite, and amphibole) surrounding omphacite. In summary, this study not only illustrates the application of thermodynamic modelling in quantifying metamorphic processes, but also the need of comparison between modeling predictions and petrographic observations. [1] White et al. (2007), J Metamorph Geol 25, 511–527. [2] Pan et al. (2020), Contrib Mineral Petrol 175, 1–28. [3] St-Onge et al. (2013), J Metamorph Geol 31, 469–504.« less
  2. We quantify the metamorphic pressure-temperature (P-T) conditions for a newly discovered silica-undersaturated high-pressure granulite (HPG) from the Central Maine Terrane (CMT) in northeastern Connecticut, U.S.A. The rocks lie within the Acadian-Neoacadian orogenic belt (Devonian) and form part of the Brimfield Schist. The Brimfield and the adjacent Bigelow Brook Formation contain silica-saturated rocks that have previously been shown to have undergone ~1000 °C metamorphism. The pressure was less well constrained at ≥ ~1 GPa. Silica-undersaturated rocks hold underutilized potential for pinpointing peak metamorphic conditions, particularly pressure, because of their resilience to melting and the variety of refractory minerals they contain. Themore »typical silica-undersaturated mineral assemblage is garnet + spinel + corundum + plagioclase + K-feldspar + biotite + ilmenite. Leucosomes are syenites consisting of two feldspars ± biotite. Plagioclase is commonly antiperthitic, particularly in feldspathic domains surrounding peritectic garnet; such garnet crystals reach ~10 cm in diameter. Alkali feldspars are perthitic. The rocks contain remarkable ellipsoidal spinels as much as 5.5 cm long comprising discrete crystallographic domains hosting crystallographically oriented lamellae of a Fe-Ti phase, most likely ilmenite. Corundum is usually colorless, but can also be found as sapphire in shades of pink, purple, and blue, particularly in antiperthite-rich domains surrounding large garnets. Some sapphires are concentrically color zoned. We carried out a P-T estimation using ternary feldspar reintegration thermometry of metamorphic antiperthites together with pseudosection modeling. Samples texturally and chemically record near-eclogite facies equilibration at minimum conditions of ~1040 °C and ~1.8 GPa, establishing the CMT in northeastern CT as the first known HPG locality in the U.S. These results are consistent with high P2O5 levels found in garnet (0.18 wt%), Ti-in-biotite thermometry, regional sillimanite pseudomorphs after kyanite, and preliminary experimental work on melt inclusions in garnet (Ferrero et al. 2017). The leucosomes provide strong evidence that partial melting of silica-undersaturated rocks at HPG conditions can produce syenitic magmata. Strongly melt-depleted silica-undersaturated rocks may also be protoliths for garnet + spinel + corundum xenoliths reported from kimberlites. The presence of HPG gneisses demonstrates that the large-scale thrusts of the CMT sample the deepest roots of the orogenic belt (60–70 km), and perhaps even deeper subduction zone lithologies as well.« less
  3. Examination of a global suite of eclogite-facies metabasites and metasediments suggests that eclogites tend to exhibit reduced mineral assemblages relative to their protoliths. High-pressure rocks tend to lack sulfides and Fe3+-bearing oxides in the eclogite facies. We suggest that eclogite-facies mineral assemblages are consistent with prograde reactions that balance the oxidation of S2- or S- to S6+ by reducing Fe3+in silicates or oxides: (1)8Fe3+Si O (OH) +S2-=8Fe2+Si O +SO 2-+(H O) abc de42f The oxidation of one mole of S2-or S-is balanced by the reduction of 7 to 8 moles of Fe3+, and typical S concentrations in the oceanic crustmore »are capable of fully reducing the entire Fe3+ budget of metabasites. As most eclogite facies rocks do not preserve peak metamorphic sulfides, petrographic evidence for prograde S oxidation reactions are cryptic; however, textures associated with sulfate reduction in response to influx of external fluids are common (reaction 1 in reverse). These reactions produce Fe3+-rich phases and are observed in both metasedimentary and metabasic rocks across a range of retrograde P-T paths (blueschist to granulite facies). For example, high-P calc- schists exhibit reaction textures that suggest the breakdown of garnet and white mica to produce pyrite + chalcopyrite + epidote + biotite + magnetite. Our thermodynamic models of aS2 and aO2 at subduction zone P-T conditions suggest assemblages of this type are indicative of aO2 0.7 to 4.5 log units above the quartz-fayalite-magnetite buffer. In rehydrated eclogites, pyrite is commonly associated with the breakdown of garnet + omphacite to amphibole + pyrite. Additionally, direct precipitation of sulfide from sulfate is observed in two samples: 1) The retrograde assemblage pyrite + ilmenite + gypsum occurs in one retrogressed metagabbroic eclogite, and 2) Coronas of secondary pyrite + barite + gypsum enclose early retrograde pyrite in a retrogressed garnet blueschist. In many eclogites, S- is reduced to S2- as pyrite is replaced by pyrrhotite, chalcopyrite, and mixed valence Co-Ni sulfides. These reactions are balanced by oxidation of divalent to trivalent Fe-Co-Ni. Reactions of this type are consistent with increasing aS2 during retrograde metamorphism. Thus, ample evidence exists for oxidized S-bearing fluids released from subducting slabs.« less
  4. Interpretation of geochronological and petrological data from partially-melted granulite is challenging. However, integration of multiple chronometers and mineral assemblage diagrams (MAD) can be used to estimate the nature and duration of processes. Excellent lower-crustal exposures of garnet granulite from the Malaspina Pluton, Fiordland New Zealand provide an ideal place to employ this kitchen sink approach. We use zircon U-Pb ages from LA-ICPMS, SHRIMP-RG, and CA-TIMS, garnet Lu-Hf and Sm-Nd ages, and MAD in order to evaluate local partial melting vs. melt injection, equilibrium volumes, P-T conditions, and the duration of lower crustal thermal events. Host diorite (H), garnet-clinopyroxene reaction zonesmore »(GRZ), coarse garnet selvages, and tonalite veins provide a record of intrusion and granulite facies partial melting. Zircon U-Pb ages range from 123 to 107 Ma (all); LA-ICPMS ages contain the entire range; CA-TIMS ages range from 118.30±0.13 to 115.7±0.18 Ma; and SHRIMP-RG ages range from 121.4±2 to 109.8±1.8 Ma. The latter two techniques are interpreted to indicate primary igneous crystallization from ~119 to ~116 Ma and the youngest ~110 Ma ages are interpreted as metamorphic zircon growth. Garnet ages for ~1 cm grains are ~113 Ma (Lu-Hf & Sm-Nd) and record metamorphic growth, and <0.3 mm grains with Sm-Nd ages from 113 to 104 Ma reflect high temperature intracrystalline diffusion and isotopic closure during cooling to amphibolite facies. Zircon trace-element compositions indicate 2 distinct crystallization trends reflecting evolution of primary magma batches. MAD indicate that garnet was not in equilibrium with sampled rock compositions. Instead, garnet shows apparent equilibrium with a modeled mixture of the GRZ and the H and grew in equilibrium with an effective bulk composition that shifted toward the leucosome. This would produce the observed increase in garnet grossular content. We conclude that: Malaspina rocks from Crooked Arm preserve evidence for 2 igneous layers which evolved as discrete magmas, igneous crystallization lasted 2 to 3 m.y., granulite metamorphism peaked ~ 3 m.y. after intrusion, metamorphism lasted ≥3 m.y., cooling occurred at ~20°C/m.y., and granulite minerals equilibrated with a mixture of solid phases and melt at ~14 kbar and 920°C (based on garnet compositions and MAD).« less
  5. Abstract Garnet is a common metamorphic and igneous mineral with extensive solid solution that can be stable to mantle depths ≥400 km. High-T and/or high-P garnet may contain oriented lamellae of other minerals, most commonly simple oxides (e.g., rutile, ilmenite), apatite, and, in ultrahigh-P cases, silicates including pyroxene and amphibole. Lamellae have classically been considered to be precipitation features preserving a record of former garnet chemistry richer in the lamellae nutrients (e.g., Ti4+). Such microtextural origins in precipitation systems (e.g., alloys) have long been studied via the crystallographic orientation relationships (COR) that form between a host and a separating phase,more »and by the shape-preferred orientation (SPO) of the lamellae. Recently, however, alternative hypotheses to precipitation have been suggested that require emplacement of lamellae in garnet by fluids, or co-growth, overgrowth, or inheritance mechanisms. These hypotheses posit that lamellae cannot be used to study former garnet chemistry. Moreover, they predict that lamellae phases, SPO, and COR should differ widely between localities, as lamellae formation will be controlled by various local rock-specific factors such as fluid presence, fluid chemistry, or mineral growth sequence. On the other hand, if lamellae characteristics are largely consistent between localities, it likely reflects control by precipitation energetics, rather than external factors. There have been few comparative COR studies in geologic systems, but the integrative assessment of COR, SPO, and lamellae assemblages should fingerprint lamellae growth process. To test the precipitation and alternative hypotheses, we collected large electron backscatter diffraction (EBSD) data sets for rutile, ilmenite, and apatite lamellae in garnet from the Brimfield Schist, Connecticut (≥1000 °C metamorphism; Central Maine Terrane, U.S.A.). We analyzed these data alongside published EBSD data for rutile, ilmenite, and corundum from metapegmatites metamorphosed in the eclogite facies from the Austrian Alps (Griffiths et al. 2016). The apatite data set is the first of its kind, and reveals that apatite preferentially aligns its close-packed direction parallel to that of garnet (c-axisapatite//<111>garnet). We also recognize a rutile-garnet COR related to those in meteorites with Widmanstätten patterns that are unequivocal products of exsolution. This is the first identification of direct similarities between silicate-oxide and metal-metal COR of which we are aware. Significantly, this rutile-garnet COR is found in diverse geologic settings including Connecticut and Idaho (U.S.A.), Austria, Germany, Greece, and China over a broad range of bulk-rock compositions. Results for all lamellae minerals show that COR are largely consistent between localities and, furthermore, are shared between apatite, ilmenite, and corundum. Moreover, between 70% and 95% of lamellae have COR and there is a dominant COR for each lamellae phase. Calculations show that d-spacing ratios of host-lamellae pairs can successfully predict the most commonly observed specific COR (those COR with two or more axial alignments with the host). These results, especially similarity of COR from markedly different geologic settings and a low diversity of lamellae minerals, are fully consistent with lamellae formation by precipitation (likely via exsolution). In contrast, the alternative hypotheses remain unsupported by COR results as well as by mineralogical and petrological evidence. Lamellae with similar traits as those in this work should thus be considered precipitates formed during unmixing of garnet compositions originally stable at elevated or extreme pressures and temperatures.« less