Coupled garnet Lu–Hf and monazite U–Pb geochronology constrain early convergent margin dynamics in the Ross orogen,Antarctica

The Ross orogen of Antarctica is an extensive (>3000 km‐long) belt of deformed and metamorphosed sedimentary rocks and granitoid batholiths, which formed during convergence and subduction of palaeo‐Pacific lithosphere beneath East Gondwana in the Neoproterozoic–early Palaeozoic. Despite its prominent role in Gondwanan convergent tectonics, and a well‐established magmatic record, relatively little is known about the metamorphic rocks in the Ross orogen. A combination of garnet Lu–Hf and monazite U–Pb (measured by laser‐ablation split‐stream ICP‐MS) geochronology reveals a protracted metamorphic history of metapelites and garnet amphibolites from a major segment of the orogen. Additionally, direct dating of a common rock‐forming mineral (garnet) and accessory mineral (monazite) allows us to test assumptions that are commonly used when linking accessory mineral geochronology to rock‐forming mineral reactions. Petrography, mineral zoning, thermobarometry and pseudosection modelling reveal a Barrovian‐style prograde path, reaching temperatures of ~610–680 °C. Despite near‐complete diffusional resetting of garnet major element zoning, the garnet retains strong rare earth element zoning and preserves Lu–Hf dates that range from c. 616–572 Ma. Conversely, monazite in the rocks was extensively recrystallized, with concordant dates that span from c. 610–500 Ma, and retain only vestigial cores. Monazite cores yield dates that overlap with the garnet Lu–Hf dates and typically have low‐Y and heavy rare earth element (HREE) concentrations, corroborating interpretations of low‐Y and low‐HREE monazite domains as records of synchronous garnet growth. However, ratios of REE concentrations in garnet and monazite do not consistently match previously reported partition coefficients for the REE between these two minerals. High‐Y monazite inclusions within pristine, crack‐free garnet yield U–Pb dates significantly younger than the Lu–Hf dates for the same samples, indicating recrystallization of monazite within garnet. The recrystallization of high‐Y and high‐HREE monazite domains over >50 Ma likely records either punctuated thermal pulses or prolonged residence at relatively high temperatures (up to ~610–680 °C) driving monazite recrystallization. One c. 616 Ma garnet Lu–Hf date and several c. 610–600 Ma monazite U–Pb dates are tentatively interpreted as records of the onset of tectonism metamorphism in the Ross orogeny, with a more robust constraint from the other Lu–Hf dates (c. 588–572 Ma) and numerous c. 590–570 Ma monazite U–Pb dates. The data are consistent with a tectonic model that involves shortening and thickening prior to widespread magmatism in the vicinity of the study area. The early tectonic history of the Ross orogen, recorded in metamorphic rocks, was broadly synchronous with Gondwana‐wide collisional Pan‐African orogenies.     

Microsampling Lu–Hf geochronology on mm‐sized garnet in eclogites constrains early garnet growth and timing of tectonometamorphism in the North Qilian orogenic belt

This study presents Lu–Hf geochronology of zoned garnet in high‐P eclogites from the North Qilian orogenic belt. Selected samples have ~mm‐sized garnet grains that have been sampled with a micro‐drill and analysed for dating. The Lu–Hf dates of bulk garnet separates, micro‐drilled garnet cores and the remnant, rim‐enriched garnet were determined by two‐point isochrons, with cores being consistently older than the bulk‐ and rim‐enriched garnet. The bulk garnet separates of each sample define identical garnet–whole rock isochron date of c. 457 Ma. Consistent U–Pb zircon dates of 455 ± 8 Ma were obtained from the eclogite. The Lu–Hf dates of the drilled cores and rim‐rich separates suggest a minimum garnet growth interval of 468.9 ± 2.4 and 452.1 ± 1.6 Ma. Major and Lu element profiles in the majority of garnet grains show well‐preserved Rayleigh‐style fractionated bell‐shaped Mn and Lu zoning profiles, and increasing Mg from core to rim. Pseudosection modelling indicates that garnet grew along a P–T path from ~470–525°C and ~2.4–2.6 GPa. The exceptional high‐Mn garnet core in one sample indicates an early growth during epidote–blueschist facies metamorphism at <460°C and <0.8 GPa. Therefore, the Lu–Hf dates of drilled cores record the early prograde garnet growth, whereas the Lu–Hf dates of rim‐rich fractions provide a maximum age for the end of garnet growth. The microsampling approach applied in this study can be broadly used in garnet‐bearing rocks, even those without extremely large garnet crystals, in an attempt to retrieve the early metamorphic timing recorded in older garnet cores. Given a proper selection of the drill bit size and a detailed crystal size distribution analysis, the cores of the mm‐sized garnet in most metamorphic rocks can be dated to yield critical constraints on the early timing of metamorphism. This study provides new crucial constraints on the timing of the initial subduction (before c. 469 Ma) and the ultimate closure (earlier than c. 452 Ma) of the fossil Qilian oceanic basin.     

Evidence for protracted prograde metamorphism followed by rapid exhumation of the Zermatt‐Saas Fee ophiolite

Major and trace‐element zoning in garnet, in combination with Rb–Sr, Sm–Nd and Lu–Hf geochronology, provide evidence for a protracted garnet growth history for the Zermatt‐Saas Fee (ZSF) ophiolite, western Alps. Four new Lu–Hf ages from Pfulwe (c. 52–46 Ma) and one from Chamois (c. 52 Ma) are very similar to a previously published Lu–Hf age from Lago di Cignana. Overall, the similarity of geochronological and garnet zoning patterns suggests that these three localities had a similar prograde tectonic history, commensurate with their similar structural position near the top of the ZSF. Samples from the lower part of the ZSF at Saas Fee and St. Jacques, however, produced much younger Lu–Hf ages (c. 41–38 Ma). Neither differences in whole‐rock geochemistry, which might produce distinct garnet growth histories, nor rare‐earth‐element zoning in garnet, can account for the age differences in the two suites. This suggests a much later prograde history for the lower part of the ZSF, supporting the idea that it was subducted diachronously. Such a model is consistent with changes in subduction vectors based on plate tectonic reconstructions, where early oblique subduction, which produced long prograde garnet growth, changed to more orthogonal subduction, which corresponds to shorter prograde garnet growth. Six new Rb–Sr phengite ages range from c. 42 to 39 Ma and, in combination with previously published Rb–Sr ages, constrain the timing of the transition from eclogite to upper greenschist facies P–T conditions. The proximity of the ZSF in the Saas Fee region to the underlying continental Monte Rosa unit and the similarity of peak‐metamorphic ages suggest these two units were linked for part of their tectonic history. This in turn indicates that the Monte Rosa may have been partly responsible for rapid exhumation of the ZSF unit.     

Pseudosection modelling and garnet Lu–Hf geochronology of HP amphibole schists constrain the closure of an ocean basin between the northern and southern Lhasa blocks,central Tibet

The P–T–t path of high‐P metamorphic rocks in subduction zones may reveal valuable information regarding the tectonic processes along convergent plate boundaries. Herein, we present a detailed petrological, pseudosection modelling and radiometric dating study of several amphibole schists of oceanic affinity from the Lhasa Block, Tibet. The amphibole schists experienced an overall clockwise P–T path that was marked by post‐P_max heating–decompression and subsequent isothermal decompression following the attainment of peak high‐P and low‐T conditions (~490°C and 1.6 GPa). Pseudosection modelling shows that the amphibole schists underwent water‐unsaturated conditions during prograde metamorphism, and the stability field of the assemblage extends to lower temperatures and higher pressures within the water‐unsaturated condition relative to water‐saturated model along the prograde path. The high‐P amphibole schists were highly reduced during retrograde metamorphism. Precise evaluation of the ferric iron conditions determined from the different compositions of epidote inclusions in garnet and matrix epidote is crucial for a true P–T estimate by garnet isopleth thermobarometry. Lu–Hf isotope analyses on garnet size separates from a garnet‐bearing amphibole schist yield four two‐point garnet–whole‐rock isochron ages from 228.2 ± 1.2 Ma to 224.3 ± 1.2 Ma. These Lu–Hf dates are interpreted to constrain the period of garnet growth and approximate the timing of prograde metamorphism because of the low peak metamorphic temperature of the rock and the well‐preserved Mn/Lu growth zoning in garnet. The majority of zircon U–Pb dates provide no constraints on the timing of metamorphism; however, two concordant U–Pb dates of 222.4 ± 3.9 Ma and 223.3 ± 4.2 Ma were obtained from narrow and uncommon metamorphic rims. Coexistence of zircon and sphene in the samples implies that the metamorphic zircon growth was likely assisted by retrogression of rutile to sphene during exhumation. The near coincident radiometric dates of zircon U–Pb and garnet Lu–Hf indicate rapid burial and exhumation of the amphibole schists, suggesting a closure time of c. 224–223 Ma for the fossil ocean basin between the northern and southern Lhasa blocks.     

Decoding polyphase migmatites using geochronology and phase equilibria modelling

In this study, in situ U–Pb monazite ages and Lu–Hf garnet geochronology are used to distinguish mineral parageneses developed during Devonian–Carboniferous and Cretaceous events in migmatitic paragneiss and orthogneiss from the Fosdick migmatite–granite complex in West Antarctica. SHRIMP U–Pb monazite ages define two dominant populations at 365–300 Ma (from cores of polychronic grains, dominantly from deeper structural levels in the central and western sectors of the complex) and 120–96 Ma (from rims of polychronic grains, dominantly from the central and western sectors of the complex, and from monochronic grains, mostly from shallower structural levels in the eastern sector of the complex). For five paragneisses and two orthogneisses, Lu–Hf garnet ages range from 116 to 111 Ma, c. 12–17 Ma older than published Sm–Nd garnet ages of 102–99 Ma from three of the same samples. Garnet grains in the analysed samples generally have Lu‐enriched rims relative to Lu‐depleted cores. By contrast, for three of the same samples, individual garnet grains have flat Sm concentrations consistent with high‐T diffusive resetting. Lutetium enrichment of garnet rims is interpreted to record the breakdown of a Lu‐rich accessory mineral during the final stage of garnet growth immediately prior to the metamorphic peak, and/or the preferential retention of Lu in garnet during breakdown to cordierite in the presence of melt concomitant with the initial stages of exhumation. Therefore, garnet is interpreted to be part of the Cretaceous mineral paragenesis and the Lu–Hf garnet ages are interpreted to record the timing of close‐to‐peak metamorphism for this event. For the Devonian–Carboniferous event, phase equilibria modelling of the metasedimentary protoliths to the paragneiss and a diatexite migmatite restrict the peak P–T conditions to 720–800 °C at 0.45–1.0 GPa. For the Cretaceous event, using both forward and inverse phase equilibria modelling of residual paragneiss and orthogneiss compositions, the P–T conditions after decompression are estimated to have been 850–880 °C at 0.65–0.80 GPa. These P–T conditions occurred between c. 106 and c. 96 Ma, determined from Y‐enriched rims on monazite that record the timing of garnet and biotite breakdown to cordierite in the presence of melt. The effects of this younger metamorphic event are dominant throughout the Fosdick complex.     

Protracted garnet growth in high‐P eclogite: constraints from multiple geochronology and P–T pseudosection

Understanding convergent margin processes requires determination of the onset and the termination of subduction, the duration of subduction‐zone metamorphism, and the subduction zone polarity. Garnet growth and intracrystalline zonation can be used to constrain the timing, duration and kinetics of tectonometamorphic processes. An eclogite from the Huwan shear zone in the Hong'an orogen was investigated with combined pseudosection analysis and multiple geochronologies. The pseudosection analysis illustrates that garnet growth is continuous and along an early near‐isothermal trajectory followed by a near‐isobaric heating path from 1.9 GPa/500 °C to 2.4 GPa/575 °C and subsequent near‐isothermal decompression. ⁴⁰Ar/³⁹Ar dating of an amphibole inclusion in garnet from the eclogite yielded an age of 310 ± 5 Ma, which is consistent with a U–Pb age of 305 ± 3 Ma for the metamorphic zircon within uncertainty. Garnet core and rim material produced Lu–Hf ages of 296.9 ± 3.8 and 256.9 ± 3.9 Ma respectively; the latter is consistent with its Sm–Nd age of 254.3 ± 4.6 Ma for the same aliquots. Similarly, limited zircon U–Pb ages of c. 257 Ma were obtained in zircon rims with garnet inclusions. These ages were interpreted to bracket the period of garnet growth and the difference of up to c. 40 Ma is best explained by protracted garnet growth. We propose that the rocks represent detachment of part of the downgoing slab and remained free of significant compression/decompression or heating/cooling close to the subduction channel, most likely underplating the mantle wedge, for a long time. These rocks were incorporated into the following subduction channel due to the successive entry of the buoyant materials, and exhumed at some time later than c. 254 Ma. The increasing observations of protracted garnet growth and long‐lived subduction in various orogens worldwide demand more sophisticated geodynamic models.     

Micro‐sampling Lu–Hf geochronology reveals episodic garnet growth and multiple high‐P metamorphic events

We present Lu–Hf dates from multiple growth zones within two large garnet porphyroblasts by a micro‐sampling method. A single eclogite sample, taken from the Huwan shear zone in the Hong'an orogen, was investigated with Lu–Hf geochronology by micro‐sampling and traditional bulk separation methods. The sample contains a few large garnet porphyroblasts up to several centimetres in diameter and a second major population of smaller (0.1–2.0 mm) garnet porphyroblasts, comprised of dark cores and pale rims. Elemental compositions and mineral inclusions in the garnet appear consistent with two garnet generations. Lu–Hf dates (c. 400–264 Ma) were determined from twelve micro‐sawed garnet sections from two large garnet porphyroblasts. These Lu–Hf dates overlap with age peaks defined by 115 SIMS zircon U–Pb analyses from the same eclogite sample. Bulk analyses of smaller garnet separates define a minimum date of c. 306 Ma and a maximum date of c. 252 Ma for the dark cores and pale rims respectively. These Lu–Hf dates were interpreted to bracket the period of garnet growth and the broad interval from c. 400 to 264 Ma is best explained by protracted and episodic garnet growth, which may require that these rocks experienced two subduction cycles that were initiated during the Devonian and terminated in the Triassic.     

Integrated garnet and zircon–titanite geochronology constrains the evolution of ultra‐high–pressure terranes: An example from the Sulu orogen

Dating ultra‐high–pressure (UHP) metamorphic rocks provides important timing constraints on deep subduction zone processes. Eclogites, deeply subducted rocks now exposed at the surface, undergo a wide range of metamorphic conditions (i.e. deep subduction and exhumation) and their mineralogy can preserve a detailed record of chronologic information of these dynamic processes. Here, we present an approach that integrates multiple radiogenic isotope systems in the same sample to provide a more complete timeline for the subduction–collision–exhumation processes, based on eclogites from the Dabie–Sulu orogenic belt in eastern China, one of the largest UHP terranes on Earth. In this study, we integrate garnet Lu–Hf and Sm–Nd ages with zircon and titanite U–Pb ages for three eclogite samples from the Sulu UHP terrane. We combine this age information with Zr‐in‐rutile temperature estimates, and relate these multiple chronometers to different P–T conditions. Two types of rutile, one present as inclusions in garnet and the other in the matrix, record the temperatures of UHP conditions and a hotter stage, subsequent to the peak pressure (‘hot exhumation') respectively. Garnet Lu–Hf ages (c. 238–235 Ma) record the initial prograde growth of garnet, while coupled Sm–Nd ages (c. 219–213 Ma) reflect cooling following hot exhumation. The maximum duration of UHP conditions is constrained by the age difference of these two systems in garnet (c. 235–220 Ma). Complementary zircon and titanite U–Pb ages of c. 235–230 Ma and c. 216–206 Ma provide further constraints on the timing of prograde metamorphism and the ‘cold exhumation' respectively. We demonstrate that timing of various metamorphic stages can thus be determined by employing complementary chronometers from the same samples. These age results, combined with published data from adjacent areas, show lateral diachroneity in the Dabie–Sulu orogeny. Three sub‐blocks are thus defined by progressively younger garnet ages: western Dabie (243–238 Ma), eastern Dabie–northern Sulu (238–235 Ma) and southern Sulu terranes (225–220 Ma), which possibly correlate to different crustal slices in the recently proposed subduction channel model. These observed lateral chronologic variations in a large UHP terrane can possibly be extended to other suture zones.     

Metamorphic petrology of a high‐T/low‐P granulite terrane (Damara belt,Namibia) – Constraints from pseudosection modelling and high‐precision Lu–Hf garnet‐whole rock dating

Migmatites comprise a minor volume of the high‐grade part of the Damara orogen of Namibia that is dominated by granite complexes and intercalated metasedimentary units. Migmatites of the Southern Central Zone of the Damara orogen consist of melanosomes with garnet+cordierite+biotite+K‐feldspar, and leucosomes, which are sometimes garnet‐ and cordierite‐bearing. Field evidence, petrographic observations, and pseudosection modelling suggest that, in contrast to other areas where intrusion of granitic magmas is more important, in situ partial melting of metasedimentary units was the main migmatite generation processes. Pseudosection modelling and thermobarometric calculations consistently indicate that the peak‐metamorphic grade throughout the area is in the granulite facies (~5 kbar at ~800°C). Cordierite coronas around garnet suggest some decompression from peak‐metamorphic conditions and rare andalusite records late, near‐isobaric cooling to <650°C at low pressures of ~3 kbar. The inferred clockwise P–T path is consistent with minor crustal thickening through continent–continent collision followed by limited post‐collisional exhumation and suggests that the granulite facies terrane of the Southern Central Zone of the Damara orogen formed initially in a metamorphic field gradient of ~35–40°C/km at medium pressures. New high‐precision Lu–Hf garnet‐whole rock dates are 530 ± 13 Ma, 522.0 ± 0.8 Ma, 520.8 ± 3.6 Ma, and 500.3 ± 4.3 Ma for the migmatites that record temperatures of ~800°C. This indicates that high‐grade metamorphism lasted for c. 20–30 Ma, which is compatible with previous estimates using Sm–Nd garnet‐whole rock systematics. In previous studies on Damara orogen migmatites where both Sm–Nd and Lu–Hf chronometers have been applied, the dates (c. 520–510 Ma) agree within their small uncertainties (0.6–0.8% for Sm–Nd and 0.1–0.2% for Lu–Hf). This implies rapid cooling after high‐grade conditions and, by implication, rapid exhumation at that time. The cause of the high geothermal gradient inferred from the metamorphic conditions is unknown but likely requires some extra heat that was probably added by intrusion of magmas from the lithospheric mantle, i.e., syenites that have been recently re‐dated at c. 545 Ma. Some granites derived from the lower crust at c. 545 Ma are the outcome rather than the cause of high‐T metamorphism. In addition, high contents of heat‐producing elements K, Th, and U may have raised peak temperatures by 150–200°C at the base of the crust, resulting in the widespread melting of fertile crustal rocks. The continuous gradation from centimetre‐scale leucosomes to decametre‐scale leucogranite sheets within the high‐grade metamorphic zone suggests that leucosome lenses coalesced to form larger bodies of anatectic leucogranites, thereby documenting a link between high‐grade regional metamorphism and Pan‐African magmatism. In view of the close association of the studied high‐T migmatites with hundreds of synmetamorphic high‐T granites that invaded the terrane as metre‐ to decametre‐wide sills and dykes, we postulate that crystallization of felsic lower crustal magma is, at least partly, responsible for heat supply. Late‐stage isobaric cooling of these granites may explain the occurrence of andalusite in some samples.     

Implications of garnet resorption for the Lu–Hf garnet geochronometer: an example from the contact aureole of the Makhavinekh Lake Pluton,Labrador

In the contact aureole of the Makhavinekh Lake Pluton (MLP), Labrador, garnet resorption caused redistribution of Lu and loss of Hf, creating spuriously young Lu–Hf garnet ages. Garnet grew during granulite facies regional metamorphism at 1860–1850 Ma. At 1322 Ma, garnet rims were replaced by coronas of cordierite and orthopyroxene during contact metamorphism. Garnet–ilmenite Lu–Hf geochronology using bulk‐garnet separates yields apparent ages that young from 1876 ± 21 Ma at 4025 m from the contact to 1396 ± 8 Ma at 450 m from the contact. Toward the contact, garnet crystals are progressively more resorbed. Concentrations of Lu measured by LA‐ICP‐MS along radial traverses on central sections through relict garnet decrease gently away from the cores but rise steeply within 50–200 μm of the edges of the relicts. Enrichments of Lu in rims of relict garnet demonstrates strong partitioning of Lu into garnet during resorption and modest intracrystalline diffusion. Hafnium distributions could not be measured, but considering the strong incompatibility of Hf with garnet, it is likely that nearly all Hf in resorbed portions of the garnet was lost from the crystals. Lu–Hf ages in the aureole are thus controlled predominantly by this retention of Lu and loss of Hf during garnet resorption. This deduction was tested with a simple numerical model in which the partial retention of Lu and loss of Hf is tracked as a population of garnet is resorbed. Assuming a spherical geometry for garnet porphyroblasts, Rayleigh fractionation is used to approximate initial Lu zoning profiles ranging from flat to steeply decreasing toward garnet rims. The model simulates: (i) Lu–Hf decay for a specified period before resorption; (ii) instantaneous resorption with retention of Lu and loss of Hf from the resorbed portion of the crystal and (iii) Lu–Hf decay during a specified period after resorption. Several parameters influence the modelled age, but garnet resorption and Lu retention are the primary factors. When all other parameters are held constant, larger amounts of resorption and higher degrees of Lu retention produce younger apparent ages (false ages). Similarly, flatter initial Lu profiles yield younger apparent ages as a consequence of the larger proportion of Lu and Hf that resides in the outer portions of the porphyroblast. The difference between the apparent and actual ages is greater if the duration of the pre‐resorption decay period is large relative to the post‐resorption decay period. Larger crystals in a Gaussian crystal‐size distribution (CSD) generally dominate the Lu–Hf budget and produce an older apparent age relative to the age of the mean crystal size. Compared to a symmetrical Gaussian CSD, positively skewed CSDs result in reduced resorption of large crystals and produce an older apparent age. Application of the model to the MLP aureole, positing growth at 1850 Ma and resorption at 1320 Ma, yields model ages that young from 1850 to 1374 Ma toward the contact, in good agreement with the apparent ages determined from geochronology.     

Time‐scale of deformation and intertectonic phases revealed by P–T–D–t relationships in the orogenic middle crust of the Orlica‐Śnieżnik Dome,Polish/Czech Central Sudetes

A section of the orogenic middle crust (Orlica‐?nie?nik Dome, Polish/Czech Central Sudetes) was examined to constrain the duration and significance of deformation (D) and intertectonic (I) phases. In the studied metasedimentary synform, three deformation events produced an initial subhorizontal foliation S1 (D₁), a subsequent subvertical foliation S2 (D₂) and a late subhorizontal axial planar cleavage S3 (D₃). The synform was intruded by pre‐, syn‐ and post‐D₂ granitoid sheets. Crystallization–deformation relationships in mica schist samples document I_1–2 garnet–staurolite growth, syn‐D₂ staurolite breakdown to garnet–biotite–sillimanite/andalusite, I_2–3 cordierite blastesis and late‐D₃ chlorite growth. Garnet porphyroblasts show a linear Mn–Ca decrease from the core to the inner rim, a zone of alternating Ca–Y‐ and P‐rich annuli in the inner rim, and a Ca‐poor outer rim. The Ca–Y‐rich annuli probably reflect the occurrence of the allanite‐to‐monazite transition at conditions of the staurolite isograd, whereas the Ca‐poor outer rim is ascribed to staurolite demise. The reconstructed P–T path, obtained by modelling the stability of parageneses and garnet zoning, documents near‐isobaric heating from ~4 kbar/485 °C to ~4.75 kbar/575 °C during I_1–2. This was followed by a progression to 4–5 kbar/580–625 °C and a subsequent pressure decrease to 3–4 kbar during D₂. Pressure decrease below 3 kbar is ascribed to I_2–3, whereas cooling below ~500 °C occurred during D₃. In the dated mica schist sample, garnet rims show strong Lu enrichment, oscillatory Lu zoning and a slight Ca increase. These features are also related to allanite breakdown coeval with staurolite appearance. As Lu‐rich garnet rims dominate the Lu–Hf budget, the 344 ± 3 Ma isochron age is ascribed to garnet crystallization at staurolite grade, near the end of I_1–2. For the dated sample of amphibole–biotite granitoid sheet, a Pb–Pb single zircon evaporation age of 353 ± 1 Ma is related to the onset of plutonic activity. The results suggest a possible Devonian age for D₁, and a Carboniferous burial‐exhumation cycle in mid‐crustal rocks that is broadly coeval with the exhumation of neighbouring HP rocks during D₂. In the light of published ages, a succession of telescoping stages with time spans decreasing from c. 10 to 2–3 Ma is proposed. The initially long period of tectonic quiescence (I_1–2 phase, c. 10 Ma) inferred in the middle crust contrasts with contemporaneous deformation at deeper levels and points to decoupled P–T–D histories within the orogenic wedge. An elevated gradient of ~30 °C km^?1 and assumed high heating rates of c. 20 °C Ma^?1 are explained by the protracted intrusion of granitoid sheets, with or without deformation, whereas fast vertical movements (2–3 Ma, D₂ phase) in the crust require the activity of deformation phases.     

Timing of high-pressure metamorphism and exhumation of the eclogite type-locality (Kupplerbrunn–Prickler Halt, Saualpe, south-eastern Austria): constraints from correlations of the Sm–Nd, Lu–Hf, U–Pb and Rb–Sr isotopic systems

Sm–Nd, Lu–Hf, Rb–Sr and SIMS U–Pb data are presented for meta‐gabbroic eclogites from the eclogite type‐locality ( Haüy, 1822 ) Kupplerbrunn–Prickler Halt and other areas of the Saualpe (SE Austria) and Pohorje Mountains (Slovenia). Mg‐rich eclogites derived from early gabbroic cumulates are kyanite‐ and zoisite rich, whereas eclogites with lower Mg contents contain clinozoisite ± kyanite. Calculated P–T conditions at the final stages of high‐pressure metamorphism are 2.2 ± 0.2 GPa at 630–740 °C. Kyanite‐rich eclogites did not yield geologically meaningful Sm–Nd ages due to incomplete Nd isotope equilibration, whereas Sm–Nd multifraction garnet–omphacite regression for a low‐Mg eclogite from Kupplerbrunn yields an age of 91.1 ± 1.3 Ma. The Sm–Nd age of 94.1 ± 0.8 Ma obtained from the Fe‐rich core fraction of this garnet dates the initial stages of garnet growth. Zircon that also crystallized at eclogite facies conditions gives a weighted mean U–Pb SIMS age of 88.4 ± 8.1 Ma. Lu–Hf isotope analysis of a kyanite–eclogite from Kupplerbrunn yields 88.4 ± 4.7 Ma for the garnet–omphacite pair. Two low‐Mg eclogites from the Gertrusk locality of the Saualpe yield a multimineral Sm–Nd age of 90.6 ± 1.0 Ma. A low‐Mg eclogite from the Pohorje Mountains (70 km to the SE) gives a garnet–whole‐rock Lu–Hf age of 93.3 ± 2.8 Ma. These new age data and published Sm–Nd ages of metasedimentary host rocks constrain the final stages of the eo‐Alpine high‐pressure event in the Saualpe–Pohorje part of the south‐easternmost Austroalpine nappe system suggesting that garnet growth in the high‐pressure assemblages started at c. 95–94 Ma and ceased at c. 90–88 Ma, probably at the final pressure peak. Zircon and amphibole crystallization was still possible during incipient isothermal decompression. Rapid exhumation of the high‐pressure rocks was induced by collision of the northern Apulian plate with parts of the Austroalpine microplate, following Jurassic closure of the Permo‐Triassic Meliata back‐arc basin.     

The timing of migmatization in the northern Arabian–Nubian Shield: Evidence for a juvenile sedimentary component in collision‐related batholiths

Collision‐related granitoid batholiths, like those of the Hercynian and Himalayan orogens, are mostly fed by magma derived from metasedimentary sources. However, in the late Neoproterozoic calcalkaline (CA) batholiths of the Arabian–Nubian Shield (ANS), which constitutes the northern half of the East African orogen, any sedimentary contribution is obscured by the juvenile character of the crust and the scarcity of migmatites. Here, we use paired in situ LASS‐ICP‐MS measurements of U–Th–Pb isotope ratios and REE contents of monazite and xenotime and SHRIMP‐RG analyses of separated zircon to demonstrate direct linkage between migmatites and granites in the northernmost ANS. Our results indicate a single prolonged period of monazite growth at 640–600 Ma, in metapelites, migmatites and peraluminous granites of three metamorphic suites: Abu‐Barqa (SW Jordan), Roded (S Israel) and Taba–Nuweiba (Sinai, Egypt). The distribution of monazite dates and age zoning in single monazite grains in migmatites suggest that peak thermal conditions, involving partial melting, prevailed for c. 10 Ma, from 620 to 610 Ma. REE abundances in monazite are well correlated with age, recording garnet growth and garnet breakdown in association with the prograde and retrograde stages of the melting reactions, respectively. Xenotime dates cluster at 600–580 Ma, recording retrogression to greenschist facies conditions as garnet continued to destabilize. Phase equilibrium modelling and mineral thermobarometry yield P–T conditions of ~650–680°C and 5–7 kbar, consistent with either water‐fluxed or muscovite‐breakdown melting. The expected melt production is 8–10 vol.%, allowing a melt connectivity network to form leading to melt segregation and extraction. U–Pb ages of zircon rims from leucosomes indicate crystallization of melt at 610 ± 10 Ma, coinciding with the emplacement of a vast volume of CA granites throughout the northern ANS, which were previously considered post‐collisional. Similar monazite ages (c. 620 Ma) retrieved from the amphibolite facies Elat schist indicate that migmatites are the result of widespread regional rather than local contact metamorphism, representing the climax of the East African orogenesis.     

Garnet re‐equilibration by coupled dissolution–reprecipitation: evidence from textural,major element and oxygen isotope zoning of ‘cloudy’ garnet

The analysis of texture, major element and oxygen isotope compositions of cloudy garnet crystals from a metapelite sampled on Ikaria Island (Greece) is used to assess the model of growth and re‐equilibration of these garnet crystals and to reconstruct the pressure–temperature–fluid history of the sample. Garnet crystals show complex textural and chemical zoning. Garnet cores (100–200 μm) are devoid of fluid inclusions. They are characterized by growth zoning demonstrated by a bell‐shaped profile of spessartine component (7–3 mol.%), an increase in grossular from 14 to 22 mol.% and δ¹⁸O values between 9.5 ± 0.3‰ and 10.4 ± 0.2‰. Garnet inner rims (90–130 μm) are fluid inclusion‐rich and show a decreasing grossular component from 22 to 5 mol.%. The trend of the spessartine component observed in the inner rim allows two domains to be distinguished. In contrast to domain I, where the spessartine content shows the same trend as in the core, the spessartine content of domain II increases outwards from 2 to 14 mol.%. The δ¹⁸O values decrease towards the margins of the crystals to a lowest value of 7.4 ± 0.2‰. The outer rims (<10 μm) are devoid of fluid inclusions and have the same chemical composition as the outermost part of domain II of the inner rim. Garnet crystals underwent a four‐stage history. Stage 1: garnet growth during the prograde path in a closed system for oxygen. Garnet cores are remnants of this growth stage. Stage 2: garnet re‐equilibration by coupled dissolution–reprecipitation at the temperature peak (630 < T < 650 °C). This causes the creation of porosity as the coupled dissolution–reprecipitation process allows chemical (Ca) and isotopic (O) exchange between garnet inner rims and the matrix. The formation of the outer rim is related to the closure of porosity. Stage 3: garnet mode decreases during the early retrograde path, but garnet is still a stable phase. The resulting garnet composition is characterized by an increasing Mn content in the inner rim’s domain II caused by intracrystalline diffusion. Stage 4: dissolution of garnet during the late retrograde path as garnet is not a stable phase anymore. This last stage forms corroded garnet. This study shows that coupled dissolution–reprecipitation is a possible re‐equilibration process for garnet in metamorphic rocks and that intra‐mineral porosity is an efficient pathway for chemical and isotopic exchange between garnet and the matrix, even for otherwise slow diffusing elements.     

Metamorphic evolution of a newly identified Mesoproterozoic oceanic slice in the Yuka terrane and its implications for a multi‐cyclic orogenic history of the North Qaidam UHPM belt

The North Qaidam Orogenic Belt (NQOB), lying at the northern margin of the Tibet Plateau, records two orogenic cycles: A Proterozoic cycle related to the amalgamation and breakup of the supercontinent Rodinia, and an Early Palaeozoic cycle including oceanic subduction and continental deep subduction. At present, the only information about the Proterozoic cycle is the concurrent c. 1,000–900 Ma magmatic and metamorphic events, which limited the understanding of the Proterozoic evolution of NQOB and the relationship between the Qaidam Block and other Rodinia fragments. In this study, a kyanite‐bearing eclogite was identified in Yuka terrane. It has positive‐slope chondrite‐normalized rare earth element distribution patterns, similar to present‐day N‐MORB. LA–ICP–MS zircon U–Pb dating obtained a protolith age of 1,273 Ma and an eclogite facies metamorphic age of 437 Ma, which is similar to the continental deep subduction age of the Yuka terrane. Zircon Lu–Hf analysis show that the magmatic zircon cores have high εHf(t) of 8.36–15.98 and T_DM1 of 1,450–1,131 Ma (M = 1,303 ± 55 Ma, consistent with its protolith age within error), indicating a juvenile crust protolith of the eclogite. The MORB‐like whole‐rock composition and zircon U–Pb and Lu–Hf analysis indicate that the protolith of the kyanite‐bearing eclogite was a Mesoproterozoic oceanic slice. P–T pseudosection analysis shows that the kyanite‐bearing eclogite experienced four metamorphic stages: (1) a prograde stage with the assemblage garnet+omphacite+talc+lawsonite+phengite+quartz at 22.4–23.2 kbar and 585°C; (2) a peak stage with the assemblage garnet+omphacite+lawsonite+phengite+coesite at 32.5 kbar and 670°C; (3) an early retrograde stage with the assemblage garnet+omphacite+kyanite+phengite+coesite/quartz±lawsonite at 27.1–30.0 kbar and 670–690°C; and (4) a late retrograde stage with the assemblage garnet+omphacite+epidote+hornblende+phengite+quartz at <18.0 kbar. The established clockwise P–T path is similar with other continental‐type eclogites in this area. On the basis of the geochemical and geochronological data, as well as the P–T path, we suggest that the protolith of the kyanite‐bearing eclogite was emplaced in the active margin of the Qaidam Block during the assembly of Rodinia and underwent continental deep subduction in the Early Palaeozoic. We conclude that (1) the Qaidam Block participated in the assembly of the Rodinia supercontinent. It was situated at or proximal to the margin of the supercontinent and probably close to India, east Antarctica and Tarim; and (2) both Mesoproterozoic and Early Palaeozoic oceanic crust slices occur in the NQOB. Thus, special caution is needed when using the metamorphic ages of oceanic affinity eclogites without protolith ages to constrain the evolution history of the North Qaidam UHPM belt.     

Modification of the Sm–Nd isotopic system in garnet induced by retrogressive fluids

Garnet, as a major constitutive mineral of eclogite, is important for Sm–Nd dating of eclogite due to its high Sm/Nd ratio and its stability during retrogression. However, a comprehensive study of the petrography, mineral chemistry, garnet water content, and Sm–Nd isotopic composition of eclogites from the Bixiling massif, Central Dabie Zone (CDZ), reveals significant modification of the Sm–Nd isotopic system in garnet as a result of retrogression. This problem constitutes a challenge for Sm–Nd dating of the Bixiling eclogites, with the Sm–Nd isochron ages of 218 ± 4 to 210 ± 9 Ma reported in the literature being younger than 226 ± 3 Ma, which is the generally accepted peak metamorphic age of the CDZ. Petrographic analysis reveals heterogeneity in colour within single fractured garnet grains. There are light‐pink garnet (Grt‐P) and red garnet (Grt‐R) types that possess distinct chemical compositions. Compared to Grt‐P, Grt‐R has higher Fe and andradrite contents but lower Al and grossular contents. Grt‐P also has lower water contents (15–35 ppm) than Grt‐R (34–65 ppm), which, together with the spatial association between Grt‐R and fractures, suggests that the colour change is related to fluid alteration. Grt‐P is an ultra‐high‐pressure (UHP) mineral, and Grt‐R is the product of the interaction between Grt‐P and a fluid during retrogression. Moreover, Grt‐R features lower Sm and Nd contents but higher Sm/Nd ratios than Grt‐P. The Sm–Nd isochrons defined by UHP minerals (Grt‐P+Omp+Rt or Grt‐P+Cpx+WR) from three eclogite samples yield consistent ages of 226.0 ± 3.8 Ma, 225.0 ± 3.9 Ma and 226.2 ± 6.9 Ma, which are identical to the peak metamorphic age of 226 ± 3 Ma for the CDZ. The retrogressed garnet (i.e., Grt‐R), omphacite and rutile, together define a pseudoisochron with younger ages of 218.9 ± 5.9 to 202.8 ± 4.8 Ma, which are geologically meaningless. The increase in the Sm/Nd ratio with constant or lower ¹⁴³Nd/¹⁴⁴Nd ratios during the transformation of Grt‐P to Grt‐R was probably the cause of these younger ages.     

Phase equilibria modelling and LASS monazite petrochronology: P–T–t constraints on the evolution of the Priest River core complex,northern Idaho

Phase equilibria modelling, laser‐ablation split‐stream (LASS)‐ICP‐MS petrochronology and garnet trace‐element geochemistry are integrated to constrain the P–T–t history of the footwall of the Priest River metamorphic core complex, northern Idaho. Metapelitic, migmatitic gneisses of the Hauser Lake Gneiss contain the peak assemblage garnet + sillimanite + biotite ± muscovite + plagioclase + K‐feldspar ± rutile ± ilmenite + quartz. Interpreted P–T paths predict maximum pressures and peak metamorphic temperatures of ~9.6–10.3 kbar and ~785–790 °C. Monazite and xenotime ²⁰⁸Pb/²³²Th dates from porphyroblast inclusions indicate that metamorphism occurred at c. 74–54 Ma. Dates from HREE‐depleted monazite formed during prograde growth constrain peak metamorphism at c. 64 Ma near the centre of the complex, while dates from HREE‐enriched monazite constrain the timing of garnet breakdown during near‐isothermal decompression at c. 60–57 Ma. Near‐isothermal decompression to ~5.0–4.4 kbar was followed by cooling and further decompression. The youngest, HREE‐enriched monazite records leucosome crystallization at mid‐crustal levels c. 54–44 Ma. The northernmost sample records regional metamorphism during the emplacement of the Selkirk igneous complex (c. 94–81 Ma), Cretaceous–Tertiary metamorphism and limited Eocene exhumation. Similarities between the Priest River complex and other complexes of the northern North American Cordillera suggest shared regional metamorphic and exhumation histories; however, in contrast to complexes to the north, the Priest River contains less partial melt and no evidence for diapiric exhumation. Improved constraints on metamorphism, deformation, anatexis and exhumation provide greater insight into the initiation and evolution of metamorphic core complexes in the northern Cordillera, and in similar tectonic settings elsewhere.     

Zircon U–Pb ages and Hf isotope compositions of migmatites from the North Qinling terrane and their geological implications

Migmatites are predominant in the North Qinling (NQ) orogen, but their formation ages are poorly constrained. This paper presents a combined study of cathodoluminescence imaging, U–Pb age, trace element and Hf isotopes of zircon in migmatites from the NQ unit. In the migmatites, most zircon grains occur as new, homogeneous crystals, while some are present as overgrowth rims around inherited cores. Morphological and trace element features suggest that the zircon crystals are metamorphic and formed during partial melting. The inherited cores have oscillatory zoning and yield U–Pb ages of c. 900 Ma, representing their protolith ages. The early Neoproterozoic protoliths probably formed in an active continental margin, being a response to the assembly of the supercontinent Rodinia. The migmatite zircon yields Hf model ages of 1911 ± 20 to 990 ± 22 Ma, indicating that the protoliths were derived from reworking of Palaeoproterozoic to Neoproterozoic crustal materials. The anatexis zircon yields formation ages ranging from 455 ± 5 to 420 ± 4 Ma, with a peak at c. 435 Ma. Combined with previous results, we suggest that the migmatization of the NQ terrane occurred at c. 455–400 Ma. The migmatization was c. 50 Ma later than the c. 490 Ma ultra‐high‐P (UHP) metamorphism, indicating that they occurred in two independent tectonic events. By contrast, the migmatization was coeval with the granulite facies metamorphism and the granitic magmatism in the NQ unit, which collectively argue for their formation due to the northward subduction of the Shangdan Ocean. UHP rocks were distributed mainly along the northern margin and occasionally in the inner part of the NQ unit, indicating that they were exhumed along the northern edge and detached from the basement by the subsequent migmatization process.     

P–T–t evolution of garnet amphibolites in the Wutai–Hengshan area,North China Craton: insights from phase equilibria and geochronology

Garnet amphibolites can provide valuable insights into geological processes of orogenic belts, but their metamorphic evolution is still poorly constrained. Garnet amphibolites from the Wutai–Hengshan area of the North China Craton mainly consist of garnet, hornblende, plagioclase, quartz, rutile and ilmenite, with or without titanite and epidote. Four samples selected in a south–north profile were studied by the pseudosection approach in order to elucidate the characteristics of their metamorphic evolution, and to better reveal the northwards prograde change in P–T conditions as established previously. For the sample from the lower Wutai Subgroup, garnet exhibits obvious two‐substage growth zoning characteristic of pyrope (X_py) increasing but grossular (X_gr) decreasing outwards in the core, and both X_py and X_gr increasing outwards in the rim. Phase modelling using thermocalc suggests that the garnet cores were formed by chlorite breakdown over 7–9 kbar at 530–600 °C, and rims grew from hornblende and epidote breakdown over 9.5–11.5 kbar at 600–670 °C. The isopleths of the minimum An in plagioclase and maximum X_py in garnet were used to constrain the peak P–T conditions of ~11.5 kbar/670 °C. The modelled peak assemblage garnet + hornblende + epidote+ plagioclase + rutile + quartz matches well the observed one. Plagioclase–hornblende coronae around garnet indicate post‐peak decompression and fluid ingress. For the samples from the south Hengshan Complex, the garnet zoning weaken gradually, reflecting modifications during decompression of the rocks. Using the same approach, the rocks are inferred to have suprasolidus peak conditions, increasing northwards from 11.5 kbar/745 °C, 12.5 kbar/780 °C to 13 kbar/800 °C. Their modelled peak assemblages involve diopside, garnet, hornblende, plagioclase, rutile and quartz, yet diopside is not observed petrographically. The post‐peak decompression is characterized by diopside + garnet + quartz + melt = hornblende + plagioclase, causing the diopside consumption and garnet compositions to be largely modified. Thus, the pesudosection approach is expected to provide better pressure results than conventional thermobarometry, because the later approach cannot be applied with confidence to rocks with multi‐generation assemblages. U–Pb dating of zircon in the Wutai sample records a protolith age of c. 2.50 Ga, and a metamorphic age of c. 1.95 Ga, while zircon in the Hengshan samples records metamorphic ages of c. 1.92 Ga. The c. 1.95 Ga is interpreted to represent the pre‐peak or peak metamorphic stages, and the ages of c. 1.92 Ga are assigned to represent the cooling stages. All rocks in the Wutai–Hengshan area share similar clockwise P–T morphologies. They may represent metamorphic products at different crustal depths in one orogenic event, which included a main thickening stage at c. 1.95 Ga followed by a prolonged uplift and cooling after 1.92 Ga.     

Petrology,geochemistry and P–T–t path of lawsonite‐bearing retrograded eclogites in the Changning–Menglian orogenic belt,southeast Tibetan Plateau

The Changning–Menglian orogenic belt (CMOB) in the southeastern Tibetan Plateau, is considered as the main suture zone marking the closure of the Palaeo‐Tethys Ocean between the Indochina and Sibumasu blocks. Here, we investigate the recently discovered retrograded eclogites from this suture zone in terms of their petrological, geochemical and geochronological features, with the aim of constraining the metamorphic evolution and protolith signature. Two types of metabasites are identified: retrograded eclogites and mafic schists. The igneous precursors of the retrograded eclogites exhibit rare earth element distribution patterns and trace element abundance similar to those of ocean island basalts, and are inferred to have been derived from a basaltic seamount in an intra‐oceanic tectonic setting. In contrast, the mafic schists show geochemical affinity to arc‐related volcanics with the enrichment of Rb, Th and U, and depletion of Nb, Ta, Zr, Hf and Ti, and their protoliths possibly formed at an active continental margin tectonic setting. Retrograded eclogites are characterized by peak metamorphic mineral assemblages of garnet, omphacite, white mica, lawsonite and rutile, and underwent five‐stage metamorphic evolution, including pre‐peak prograde stage (M₁) at 18–19 kbar and 400–420°C, peak lawsonite‐eclogite facies (M₂) at 24–26 kbar and 520–530°C, post‐peak epidote–eclogite facies decompression stage (M₃) at 13–18 kbar and 530–560°C, subsequent amphibolite facies retrogressive stage (M₄) at 8–10 kbar and 530–600°C, and late greenschist facies cooling stage (M₅) at 5–8 kbar and 480–490°C. Laser ablation inductively coupled plasma mass spectrometry (LA–ICP–MS) U–Pb spot analyses of zircon show two distinct age groups. The magmatic zircon from both the retrograded eclogite and mafic schist yielded protolith ages of 451 ± 3 Ma, which is consistent with the ages of Early Palaeozoic ophiolitic complexes and ocean island sequences in the CMOB reported in previous studies. In contrast, metamorphic zircon from the retrograded eclogite samples yielded consistent Triassic metamorphic ages of 246 ± 2 and 245 ± 2 Ma, which can be interpreted as the timing of closure of the Palaeo‐Tethys Ocean. The compatible peak metamorphic mineral assemblages, P–T–t paths and metamorphic ages, as well as the similar protolith signatures for the eclogites in the CMOB and Longmu Co–Shuanghu suture (LCSS) suggest that the two belts formed part of a cold oceanic subduction system in the Triassic. The main suture zone of the Palaeo‐Tethyan domain extends at least 1,500 km in length from the CMOB to the LCSS in the Tibetan Plateau. The identification of lawsonite‐bearing retrograded eclogites in the CMOB provides important insights into the tectonic framework and complex geological evolution of the Palaeo‐Tethys.