An official website of the United States government
Abstract Ultrahigh‐temperature (UHT; >900°C) metamorphism drives crustal differentiation and is widely recognized in the rock record, but its geodynamic causes are debated. Previous work on granulite‐facies metapelite xenoliths from San Luis Potosí, Mexico suggests the lower crust experienced a protracted UHT metamorphic event that coincided with the onset of regional extension. To determine the duration, conditions, and heat sources of UHT metamorphism recorded by these xenoliths, this study characterizes the major‐element, trace‐element, and U‐Pb isotopic systematics of quartz, rutile, feldspar, garnet, and zircon by in situ electron microprobe (EPMA) and laser‐ablation inductively coupled‐plasma mass spectrometry (LA‐ICP‐MS), and augments these data with detailed petrography, thermobarometry, phase equilibria modeling, and diffusion modeling. Thermobarometry and phase equilibria modeling suggest peak metamorphic conditions exceeded 0.7 GPa and 900°C. Zircon petrochronology confirms >15 Myr of UHT conditions since its onset at ∼30 Ma. A small population of zircon record elevated temperatures following transition from backarc compression to extension during the waning stages of orogenesis (60–37 Ma). Garnet preserves trace‐element zoning and mineral inclusions consistent with suprasolidus garnet growth and subsequent compositional modification by intracrystalline rare‐earth element diffusion during protracted heating, with diffusion chronometry timescales in agreement with zircon data, followed by fluid‐driven remobilization of trace elements along now‐healed fractures within ∼1 Myr of eruption. In sum, these data are most compatible with lithospheric mantle attenuation or removal as the dominant heat transport mechanism driving synextensional UHT metamorphism and crustal melting, which has bearing on models for crustal differentiation and formation of modern and ancient granulite terranes globally.
more »
« less
This content will become publicly available on October 1, 2026
Coupled Diffusion and Stress in Garnet from Holsnøy, Norway: Implications for Mineral Replacement Reactions and the Timescales of Episodic Heating and Fluid Flow during Subduction
Abstract Holsnøy, Norway, offers a world-class natural laboratory for studying the impact of fluid on subducting lower crust. Holsnøy is composed of dry, metastable lower crustal granulite that was infiltrated by fluids along shear zones and seismic fractures during subduction. The infiltration facilitated the localized growth of eclogite facies mineral assemblages along the fluid flow pathways. The duration of the eclogite facies metamorphism, however, remains uncertain. Previous garnet diffusion chronometry studies have estimated timescales ranging from hundreds of years to millions of years based on diffusional relaxation between metastable granulite facies garnet cores and eclogite facies garnet rims and fractures. The shorter timescales are inferred from extremely sharp Ca gradients across chemical contacts present in some garnets whereas the longer timescales are from wider Mg and Fe profiles present in all garnets. The different timescale estimates have led to divergent models for the region’s tectonometamorphic evolution. Here we show that the sharp Ca contacts can be explained by diffusion-induced compositional stress. As Ca is significantly larger than Mg and Fe, its movement strains the crystal lattice and generates stress that limits the relaxation of sharp chemical contacts. When compositional stress is accounted for, the sharp contacts yield timescales that are consistent with the wider Mg and Fe diffusion profiles. We determine that eclogite facies conditions (670–700 °C, 1.5–2.2 GPa) lasted a maximum of c. 300 kyr. The relatively short duration of eclogite facies conditions requires that multiple transient heating events were superimposed on a longer (>106 yr) overall timescale of metamorphism. Granulite facies garnet cores are surrounded by multiple generations of eclogite facies rims formed by interface-coupled dissolution–reprecipitation (ICDR) reactions. The garnet rims indicate two rapid, regional-scale fluid pulses and additional smaller, more localized pulses. The fluid pulses may be linked to episodes of seismic moment release as well as transient heating via exothermic hydration reactions and/or shear deformation. Our model results predict up to 400 MPa of differential stress at the garnet core–rim contacts, consistent with observed eclogite facies microfractures that extend into relic granulite facies garnet cores. The microfractures indicate that ICDR was aided by compositional stress: diffusion ahead of the reaction front generated stress and fracturing that created porosity for further ICDR. Thus, compositional stress can markedly impact both diffusion systematics and intracrystalline deformation. Together, these results show that despite their brevity, transient thermal, fluid flux, and/or baric episodes may exert the primary controls on the mineralogical and rheological development of subducted lithologies, in contrast to the long, slow burial and exhumation typically envisioned for regional metamorphism.
more »
« less
- Award ID(s):
- 2208229
- PAR ID:
- 10642887
- Publisher / Repository:
- Journal of Petrology, Oxford University Press
- Date Published:
- Journal Name:
- Journal of Petrology
- Volume:
- 66
- Issue:
- 10
- ISSN:
- 0022-3530
- Page Range / eLocation ID:
- 1-29
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Garnet ages for eclogite and granulite from the Western Fiordland Orthogneiss (WFO) provide a precise age for high-grade metamorphism and partial melting of the lower crust in a Cretaceous magmatic arc currently exposed in Fiordland, New Zealand. U/Pb zircon ages and pluton areas indicate that a high magmatic flux event between 118 and 115 Ma added >3,000 km2 of mid- to lower-crustal plutons. The high flux event was followed by high temperature metamorphism and partial melting which resulted in pervasive leucosomes, and trondhjemite layers and veins. At least 1,800 km2 of the newly added crust was metamorphosed to garnet granulite facies orthogneiss. Thermobarometry and phase diagram models indicate that garnet grew at 850 to 1,000°C and 12 to 14 kbar in this monzodiorite and diorite gneiss of the Misty, Malaspina, and Breaksea plutons. Sm-Nd garnet-rock isochrons for these three plutons of the WFO (>700 km2of lower crust) indicate that peak temperatures were reached at 111.7±1.0 Ma (N=16). The isotopic and chemical composition of zircon indicate that the Cretaceous arc flare-up was most likely triggered by partial melting and hybridization of subducted oceanic crust and enriched subcontinental lithospheric mantle directly prior to cessation of arc magmatism. The driving mechanism for the terminal magmatic surge is inferred to be propagation of a discontinuous slab tear beneath the arc, or a ridge-trench collision event between 136 and 128 Ma. The lack of ca. 112 Ma plutons in the western part of Fiordland negates a magmatic heat source for garnet granulite metamorphism. Therefore, we infer that high heat flow associated with mantle advection at the base of the arc after the magmatic surge continued for several m.y., heating the lower crust to granulite facies temperaturesmore » « 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
-
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 crust 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.more » « less
-
Zircon U-Pb, and garnet Sm-Nd and Lu-Hf dates provide important constraints on local and orogenic scale processes in lower-crustal rocks. However, in high-temperature metamorphic rocks these isotopic systems typically yield significant ranges reflecting both igneous and metamorphic processes. Therefore, linking dates to specific aspects of rock history can be problematic. In Fiordland, New Zealand, granulite-facies orthogneiss is cut by leucosomes that are bordered by garnet clinopyroxene reaction zones (garnet reaction zones). In both host orthogneiss and garnet reaction zones, zircon are typically anhedral with U-Pb dates ranging from 118.30 ± 0.13 to 115.70 ± 0.18 Ma (CA-ID-TIMS) and 121.4 ± 2.0 to 109.8 ± 1.8 Ma (SHRIMP-RG). Zircon dates in host and garnet reaction zone do not define distinct populations. In addition, the dates cannot be readily grouped based on external morphology or internal CL zoning. Zircon trace-element concentrations indicate two distinct crystallization trends, clearly seen in Th and U. Garnet occurs in selvages to the leucosome veins and in the adjacent garnet reaction zones. In selvages and host orthogneiss, garnet is generally 0.5 to 1 cm diameter and euhedral and is 0.1 to 0.5 cm diameter and subhedral in garnet reaction zones. Garnet Sm-Nd and Lu-Hf dates range from ca. 115 to 101 Ma (including uncertainties) and correlate with grain size. We interpret the CA-ID-TIMS zircon dates to record the age of magma emplacement and the SHRIMP-RG dates to record a range from igneous crystallization to metamorphic dissolution and reprecipitation and/or local Pb loss. Zircon compositional trends within the garnet reaction zone and host are compatible with locally isolated melt and/or separate intrusive magma batches for the two samples described here. Dates for the largest, ~1 cm, garnet of ~113 Ma record growth during metamorphism, while the smaller grains with younger dates reflect high-temperature intracrystalline diffusion and isotopic closure during cooling. The comprehensive geochronological data set for a single location in the Malaspina Pluton illustrates a complex and protracted geologic history common in granulite facies rocks, estimates lower crustal cooling rates of ~20 °C/m.y., and underlines the importance of multiple chronometers and careful textural characterization for assigning meaningful ages to lower-crustal rocks. Numerous data sets from single locations, like the one described here, are needed to evaluate the spatial extent and variation of cooling rates for Fiordland and other lower crustal exposures.more » « less
