skip to main content

Title: Thermodynamic modelling on the UHP metamorphism and fluid infiltration of the Tso Morari coesite-bearing eclogite in NW India
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-T 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) more » 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
; ;
Award ID(s):
Publication Date:
Journal Name:
Annual VM Goldschmidt Conference
Sponsoring Org:
National Science Foundation
More Like this
  1. 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 significantlymore »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 °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
  2. 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
  3. Abstract. Retrograde metamorphic rocks provide key insights into the pressure–temperature (P–T) evolution of exhumed material, and resultant P–T constraints have direct implications for the mechanical and thermal conditions of subduction interfaces. However, constraining P–T conditions of retrograde metamorphic rocks has historically been challenging and has resulted in debate about the conditions experienced by these rocks. In this work, we combine elastic thermobarometry with oxygen isotope thermometry to quantify the P–T evolution of retrograde metamorphic rocks of the Cycladic Blueschist Unit (CBU), an exhumed subduction complex exposed on Syros, Greece. We employ quartz-in-garnet and quartz-in-epidote barometry to constrain pressures of garnetmore »and epidote growth near peak subduction conditions and during exhumation, respectively. Oxygen isotope thermometry of quartz and calcite within boudin necks was used to estimate temperatures during exhumation and to refine pressure estimates. Three distinct pressure groups are related to different metamorphic events and fabrics: high-pressure garnet growth at ∼1.4–1.7 GPa between 500–550 ∘C, retrograde epidote growth at ∼1.3–1.5 GPa between 400–500 ∘C, and a second stage of retrograde epidote growth at ∼1.0 GPa and 400 ∘C. These results are consistent with different stages of deformation inferred from field and microstructural observations, recording prograde subduction to blueschist–eclogite facies and subsequent retrogression under blueschist–greenschist facies conditions. Our new results indicate that the CBU experienced cooling during decompression after reaching maximum high-pressure–low-temperature conditions. These P–T conditions and structural observations are consistent with exhumation and cooling within the subduction channel in proximity to the refrigerating subducting plate, prior to Miocene core-complex formation. This study also illustrates the potential of using elastic thermobarometry in combination with structural and microstructural constraints, to better understand the P–T-deformation conditions of retrograde mineral growth in high-pressure–low-temperature (HP/LT) metamorphic terranes.« less
  4. The large range in oxidation states of sulfur (-II to +VI) provides it with a large oxidation potential in rocks, even at relatively low concentrations. Most importantly, the transition from sulfide to sulfate species in rocks and silicate melts occurs in the same approximate fO2 region (for a given temperature) as the transition from ferrous to ferric iron, and reduced S species can coexist with oxidized Fe and vice versa. The result is a large potential for reactions involving sulfur to oxidize or reduce Fe in silicate minerals, since Fe only occurs in two oxidation states (+II and +III). Inmore »order for sulfur to be released during slab dehydration, sulfur in sulfide must be converted into an easily dissolved species, such as SO42− or H2S, through either oxidation or reduction. We propose that oxidation of sulfur in sulfide follows the generalized reaction: 8Fe3+SiaOb(OH)c +S2− = 8Fe2+SidOe +SO42− +(H2O)f (1) In this type of reaction, sulfur participates in the dehydration of greenschist- or blueschist-facies hydrous silicates during transition to the eclogite facies: ferric Fe in Fe-bearing silicates (chlorite, amphibole, epidote) is reduced to ferrous Fe in anhydrous ferromagnesian silicates (pyroxene, garnet). At the same time, the reaction consumes sulfide by oxidation of S2− to produce SO42−, which is readily dissolved in the fluid produced during dehydration. Additionally, a similar redox reaction could oxidize sulfur by reducing ferric Fe in oxides. It is important to note that one mole of S has the same redox potential as 8 moles of Fe. The molar ratio of 8 moles of Fe per 1 mole of S translates to a mass ratio of approximately 14; therefore, small concentrations of sulfur can have a large impact on reduction/oxidation of the silicate assemblage. Our observations show that sulfide minerals that can be identified as primary or related to the peak metamorphic stage are rare in eclogites and restricted to inclusions in garnet, consistent with reaction (1). Thermodynamic modeling is currently underway to assess the influence of sulfur on the phase equilibria of silicate phases during high pressure metamorphism.« 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