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 °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)
more »
« less
Preferred orientations of garnet porphyroblasts reveal previously cryptic templating during nucleation
Abstract Electron back scattered diffraction data of garnet crystals from the Nelson Aureole, British Columbia and the Mosher’s Island formation, Nova Scotia, reveals that 22 garnet crystals are all oriented with one of three crystal directions parallel to the trace of the foliation plane in thin section. Structural models suggest that these relationships are due to preferential garnet nucleation onto muscovite, with the alignment of repeating rows of Al octahedra and Si tetrahedra in each leading to inheritance of garnet orientation from the muscovite. These results highlight that epitaxial nucleation may be a prevalent process by which porphyroblast minerals nucleate during metamorphism and carry implications for the role that non-classic nucleation pathways play in the crystallization of metamorphic minerals, the distribution of porphyroblasts in metamorphic rocks, and, in cases in which nucleation is the rate limiting step for crystallization, the energetics of metamorphic reactions.
more »
« less
- PAR ID:
- 10262184
- Date Published:
- Journal Name:
- Scientific Reports
- Volume:
- 11
- Issue:
- 1
- ISSN:
- 2045-2322
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Metamorphic rocks from the Connecticut Valley Trough (CVT), Vermont, and Massachusetts, have been examined using quartz‐in‐garnet (QuiG) and conventional thermobarometry, thermodynamic reaction modelling, diffusion modelling, and40Ar/39Ar thermochronology to constrain theirP–T–tpaths during Acadian metamorphism and subsequent exhumation. Numerous samples, collected in the vicinity of the Acadian domes, contain garnet porphyroblasts that display cloudy zones characterized by numerous fluid inclusions and modified garnet compositions associated with the replacement of the original garnet by biotite±muscovite±plagioclase±quartz±lowXgrs/enrichedXsps. QuiG and conventional thermobarometry constrain both the conditions of garnet nucleation and peakP–Tconditions to have occurred at ~0.85–1.05 GPa, ~550–600°C. Most notably, QuiG barometry was performed on inclusions adjacent to these reaction zones in conjunction with Gibbs method reaction modelling to reveal that these dissolution–reprecipitation reactions occurred during nearly isothermal decompression from the peakP–Tconditions to around ~0.3 GPa, 550°C. Diffusion modelling reveals that the Mn zoning profiles created during garnet resorption that accompanied decompression formed in less thanc. 3 Ma, which constrains the tectonic exhumation to have occurred at 8–10 mm/year. Subsequent cooling to 500°C occurred rapidly at a rate of 100°C/Ma, followed by slower cooling reaching 1.7°C /Ma by the mid Carboniferous. This is the first reported example of QuiG barometry revealing a multi‐stage metamorphic history and highlights the utility of this method for unravelling complex metamorphic terranes.more » « less
-
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 zones (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).more » « less
-
Sublithospheric diamonds and the inclusions they may carry crystallize in the asthenosphere, transition zone, or uppermost lower mantle (from 300 to ∼800 km), and are the deepest minerals so far recognized to form by plate tectonics. These diamonds are distinctive in their deformation features, low nitrogen content, and inclusions of these major mantle minerals: majoritic garnet, clinopyroxene, ringwoodite, CaSi perovskite, ferropericlase, and bridgmanite or their retrograde equivalents. The stable isotopic compositions of elements within these diamonds (δ11B, δ13C, δ15N) and their inclusions (δ18O, δ56Fe) are typically well outside normal mantle ranges, showing that these elements were either organic (C) or modified by seawater alteration (B, O, Fe) at relatively low temperatures. Metamorphic minerals in cold slabs are effective hosts that transport C as CO3 and H as H2O, OH, or CH4 below the island arc and mantle wedge. Warming of the slab generates carbonatitic melts, supercritical aqueous fluids, or metallic liquids, forming three types of sublithospheric diamonds. Diamond crystallization occurs by movement and reduction of mobile fluids as they pass through host mantle via fractures—a process that creates chemical heterogeneity and may promote deep focus earthquakes. Geobarometry of majoritic garnet inclusions and diamond ages suggest upward transport, perhaps to the base of mantle lithosphere. From there, diamonds are carried to Earth's surface by eruptions of kimberlite magma. Mineral assemblages in sublithospheric diamonds directly trace Earth's deep volatile cycle, demonstrating how the hydrosphere of a rocky planet can connect to its solid interior.▪Sublithospheric diamonds from the deep upper mantle, transition zone, and lower mantle host Earth's deepest obtainable mineral samples.▪Low-temperature seawater alteration of the ocean floor captures organic and inorganic carbon at the surface eventually to become some of the most precious gem diamonds.▪Subduction transports fluids in metamorphic minerals to great depth. Fluids released by slab heating migrate, react with host mantle to induce diamond crystallization, and may trigger earthquakes.▪Sublithospheric diamonds are powerful tracers of subduction—a plate tectonic process that deeply recycles part of Earth's planetary volatile budget.more » « less
-
The Appalachian Mountains expose one of the most complete deeply exhumed orogenic belts in the world. These rocks provide the opportunity to understand tectonic processes in the mid- to lower- crust that can be linked to upper crustal processes interpreted from less exhumed orogenic belts. However, 3 Paleozoic orogenies (Taconic, Neoacadian, Alleghanian) in the southern Appalachians produced a complicated thermal-metamorphic history that is poorly understood. Recently obtained monazite U-Pb ages in the western, central, and eastern Blue Ridge of Tennessee and the Carolinas range from 459 to 441 Ma, indicating that this part of the Blue Ridge preserves Taconic (Ordovician) metamorphic mineral assemblages and were not significantly reheated during Neoacadian (Devonian) or Alleghanian (Mississippian) orogenesis. Five published garnet Sm-Nd ages from the eastern Blue Ridge in Alabama and Georgia of 331 to 320 Ma indicate widespread Alleghanian metamorphism. The northwestern extent of these Alleghanian metamorphic rocks is constrained by a garnet Sm-Nd age of 357±3 Ma from NW of the transtensional Goodwater-Enitachopco fault. However, published metamorphic age constraints are lacking SE of and along strike to the NE of the Alleghanian rocks. We report new garnet Sm-Nd ages for northern Georgia that constrain the extent of the Alleghanian metamorphic rocks. Garnet-staurolite-hornblende gneiss in the Pumpkinvine Creek Formation yields an Alleghanian age of 323±3 Ma (MSWD=6.6, N=7). To the NE, garnet-muscovite-biotite gneiss from within the structural window at Brasstown Bald and migmatiticsillimanite- and spinel-bearing garnet-biotite neiss from outside the window at Blood Mountain have ages of 446±6 (MSWD=0.7, N=4) and 448±8 (MSWD=6, N=7) Ma, respectively. These 2 indistinguishable ages confirm the premetamorphic stacking of thrust sheets exposed in the structural window. Comparison of these new ages indicates post metamorphic displacement on the Allatoona fault between the Dahlonega terrane and the western Blue Ridge. Additional garnet ages spatially distributed across the Piedmont of east central Alabama and the Murphy belt of NE Georgia extent are currently in-progress. The full data set will be used to test tectonic models including possible out-of-sequence thrusting and crustal channel flow.more » « less
An official website of the United States government

