Contrasting mechanisms and timescales of subduction and exhumation as recorded by Paleoproterozoic and late Paleozoic high-pressure granulites in the North China Craton

The North China Craton records multiple metamorphic events related to supercontinent assembly during the Paleoproterozoic, forming Columbia, and again during the late Paleozoic, forming Pangea. Here we show that the Paleoproterozoic high-pressure granulites (HPGs) formed from enriched mid-ocean ridge basalt protoliths and record a clockwise pressure-temperature-time (P–T–t) path with prograde metamorphism at 7.8–10.0 kbar and 780–820 °C, peak granulite-facies metamorphism at 12–12.3 kbar and ∼860–880 °C, and retrograde metamorphism at 8.7–9.1 kbar and 850–855 °C. Subduction initiated prior to 1.90 Ga, with final collision and orogeny at 1.88 Ga, followed by post-collision/exhumation at 1.80 Ga, defining a prolonged exhumation period (∼90 m.y.) that occurred at a slow velocity of ∼0.16 ± 0.08 mm/y. Late Paleozoic HPGs are normal mid-ocean ridge basalt type and record a near clockwise P–T–t path, with peak/post-peak amphibolitefacies metamorphism at 11.0–12.5 kbar and 860–890 °C, isothermal decompression to 7.2–7.5 kbar and 810–820 °C, and retrogression to 5.5–7.2 kbar and 805–850 °C. Subduction initiated earlier than ca. 340 Ma, exhumation and uplift initiated at 335–309 Ma and continued to 297–287 Ma. The exhumation was short-lived (∼50 m.y.) and relatively fast (0.38 ± 0.14 mm/y). When compared to granulite-facies metamorphism documented in many Paleoproterozoic HPGs, late Paleozoic HPGs appear to commonly form with an initial period of steep subduction leading to eclogite-facies metamorphism, with subsequent exhumation to middle/lower levels of the crust. Our results further reveal that the exhumation velocity for supercontinent collision was facilitated and duration shortened through time, and that the exhumation mechanism might have been controlled by subduction angle, compression pressure, and temperature

Cenozoic temporal variation of crustal thickness in the Urumieh-Dokhtar and Alborz magmatic belts, Iran

We present regional variations of whole-rock Sr/Y and (La/Yb)N ratios of magmatic rocks along the Cenozoic Urumieh-Dokhtar and Alborz magmatic belts, Iran. Both the magmatic belts are located at the north of the main Zagros Neo-Tethyan suture. The Urumieh-Dokhtar magmatic belt (UDMB), which trends NW-SE for 1000 km across Iran, was characterized by the intensive volcanism and plutonism, and defined the magmatic front (MF) of the Zagros orogenic belt. The Alborz magmatic belt (AMB) is situated to the north, and characterized by less intense magmatic activity. The Alborz magmatic belt was formed behind it in the rear-arc (RA) domain. A striking feature of the both magmatic belts is the transition from normal calc-alkaline arc magmatism during the Eocene–Oligocene to adakite-like calc-alkaline magmatism during the Middle to Late Miocene–Pliocene. The late-Cenozoic magmatism of the UDMB and AMB shows higher Sr/Y and (La/Yb)N. However, it should be noted that crustal thickening event is intensive in the UDMB than AMB during Late Cenozoic. Using the composition of the Lale-Zar zircons from the SE UDMB we determined the oxygen fugacity (fO2) during zircon crystallization to be between FMQ (fayalite–magnetite–quartz buffer) -0.69 to +2.41, whereas those of the Hashroud-Teckmdash-Gormolla zircons from NW AMB range from −1.22 to +5.99. The fO2 estimates suggest relatively more oxidized conditions for the Late Cenozoic igneous rocks of the AMB. Compiled data from the UDMB and AMB intrusions show an increase in average zircon crystallization temperatures with decreasing age. These outcomes have been interpreted in terms of variation of the crustal thickness, from 30 to 35 km during Eocene-Oligocene to 40–55 km during the middle-late Miocene. We propose the increase in crustal thickness is associated with the collision between the Arabian plate and Iran and subsequent convergence during the middle-late Miocene

Crustal architecture studies in the Iranian Cadomian arc: insights into source, timing and metallogeny

The Jalal Abad magmatic rocks, situated at the southern edge of the Saghand-Bafgh-Zarand district, include a thick pile of Cadomian extrusive and pyroclastic units intruded by younger granitoid stocks. New zircon U–Pb ages show eruptions at ∼552 Ma, followed by emplacement of granodiorite at ∼537 Ma. The Jalal Abad magmatic rocks have typical high-K and shoshonitic signatures, and are characterized by enrichment in large-ion lithophile elements (LILEs) and depletion in high-field-strength elements (HFSE). Zircon ɛHf(t) from the Jalal Abad magmatic rocks ranges from +3.9 to −3.9 for volcanic rocks and −1.2 to +8.1 for granodiorite. Zircon δ18O values for the Jalal Abad are variable from +5.1 to +8.8‰, progressively higher than those of mantle-derived melts. The whole-rock ɛNd(t) values range between −7.7 to −7.4 for granodiorite, −4.6 to −3 for volcanic rocks and −6.2 to −8.2 for ignimbrites/tuff. The whole-rock Nd and zircon Hf crustal model ages (TDMC) for the Jalal Abad magmatic rocks range between 0.8 and 2.3 Ga. All of the Jalal Abad magmatic rocks have quite similar trace element patterns, and slightly different whole-rock Nd and zircon Hf isotopic composition, indicating the involvement of the thick continental crust during the formation of these rocks. Modeling of zircon Hf–O data, bulk-rock trace elements, and Sr–Nd isotopes suggest the magmas were generated by interaction of mantle–derived melts with thick continental crust through assimilation/fractional crystallization (AFC) processes. However, crustal architecture studies in the Iranian Cadomian arcs show that AFC processes were more important during the Mesoarchean–Early Neoproterozoic (3000–1000 Ma), whereas juvenile magmas became increasingly important to the Cadomian (600–500 Ma) magmatism. Early Cambrian intrusive magmas seemingly intruded sedimentary sequences in the study region and provided magmatic constituents and a heat source for hydrothermal processes and mineralization

Metasomatism of the continental crust and its impact on surface uplift: Insights from reactive‐transport modelling

High-elevation, low-relief continental plateaus are major topographic features and profoundly influence atmospheric circulation, sediment transport and storage, and biodiversity. Although orogenic surface-uplift mechanisms for modern continental plateaus near known plate margins like Tibet are well-characterized, they cannot account for examples in intracontinental settings like the Colorado Plateau. In contrast to canonical plate-tectonic uplift mechanisms, broad-scale hydration-induced metasomatism of the lower crust has been suggested to reduce its density and increase its buoyancy sufficiently to contribute to isostatic uplift. However, the relationships between key petrophysical properties in these environments are not fully quantified, which limits application of this model. Here, we develop a series of petrological models that describe the petrological and topographic effects of fluid–rock interaction in non-deforming continental crust of varying composition. We apply an open-system petrological modelling framework that utilizes reactive-transport calculations to determine the spatial and temporal scales over which mineralogic transformations take place compared with the magnitude of infiltration of aqueous fluids derived from devolatilization of subducting oceanic lithosphere. The buoyancy effect of hydration-induced de-densification is most significant for metabasic lower crust, intermediate for metapelitic crust, and minimal for granodioritic crust. We apply these results to a case study of the ~2 km-high Colorado Plateau and demonstrate that under ideal conditions, hydration of its lower–middle crust by infiltrating aqueous fluids released by the Farallon slab during Cenozoic low-angle subduction could have uplifted the plateau surface by a maximum of ~1 km over 16 Myr. However, realistically, although hydration likely has a measurable effect on surface tectonics, the uplift of orogenic plateaus is likely dominantly controlled by other factors, such as lithospheric delamination

Quantifying geological uncertainty in metamorphic phase equilibria modelling; a Monte Carlo assessment and implications for tectonic interpretations

Pseudosection modelling is rapidly becoming an essential part of a petrologist’s toolkit and often forms the basis of interpreting the tectonothermal evolution of a rock sample, outcrop, or geological region. Of the several factors that can affect the accuracy and precision of such calculated phase diagrams, “geological” uncertainty related to natural petrographic variation at the hand sample- and/or thin section-scale is rarely considered. Such uncertainty influences the sample’s bulk composition, which is the primary control on its equilibrium phase relationships and thus the interpreted pressure–temperature (P–T) conditions of formation. Two case study examples—a garnet–cordierite granofels and a garnet–staurolite–kyanite schist—are used to compare the relative importance that geological uncertainty has on bulk compositions determined via (1) X-ray fluorescence (XRF) or (2) point counting techniques. We show that only minor mineralogical variation at the thin-section scale propagates through the phase equilibria modelling procedure and affects the absolute P–T conditions at which key assemblages are stable. Absolute displacements of equilibria can approach ±1 kbar for only a moderate degree of modal proportion uncertainty, thus being essentially similar to the magnitudes reported for analytical uncertainties in conventional thermobarometry. Bulk compositions determined from multiple thin sections of a heterogeneous garnet–staurolite–kyanite schist show a wide range in major-element oxides, owing to notable variation in mineral proportions. Pseudosections constructed for individual point count-derived bulks accurately reproduce this variability on a case-by-case basis, though averaged proportions do not correlate with those calculated at equivalent peak P–T conditions for a whole-rock XRF-derived bulk composition. The main discrepancies relate to varying proportions of matrix phases (primarily mica) relative to porphyroblasts (primarily staurolite and kyanite), indicating that point counting preserves small-scale petrographic features that are otherwise averaged out in XRF analysis of a larger sample. Careful consideration of the size of the equilibration volume, the constituents that comprise the effective bulk composition, and the best technique to employ for its determination based on rock type and petrographic character, offer the best chance to produce trustworthy data from pseudosection analysis.RMP acknowledges a NERC postgraduate grant (reference number NE/H524781/1) for funding analytical work performed at the University of Oxford, UK

Editorial: Seeing convergent margin processes through metamorphism

Plate convergence can induce large-scale metamorphism and magmatism, reshape large parts of continental margins, and subsequently change regional climate and biodiversity. Metamorphic rocks in orogenic belts commonly record different metamorphic evolutions and temporal-spatial distributions at the regional scale, which are strongly influenced by convergent processes through time. In some cases, ultrahigh-pressure (UHP) and ultrahigh-temperature (UHT) metamorphic rocks are observed at both ancient and young convergent plate margins, marking the operation of extreme tectonism in the regime of plate tectonics. This Research Topic aims to understand how regional metamorphism operated at convergent plate margins through the study of field and petrographic observations, geochemical and petrological analysis, high-pressure experiments, and thermodynamic modeling. The scope is to gather new ideas and interpretations on the structure and processes of convergent plate margins

Quantifying geological uncertainty in metamorphic phase equilibria modelling; a Monte Carlo assessment and implications for tectonic interpretations

Pseudosection modelling is rapidly becoming an essential part of a petrologist's toolkit and often forms the basis of interpreting the tectonothermal evolution of a rock sample, outcrop, or geological region. Of the several factors that can affect the accuracy and precision of such calculated phase diagrams, “geological” uncertainty related to natural petrographic variation at the hand sample- and/or thin section-scale is rarely considered. Such uncertainty influences the sample's bulk composition, which is the primary control on its equilibrium phase relationships and thus the interpreted pressure–temperature (P–T) conditions of formation. Two case study examples—a garnet–cordierite granofels and a garnet–staurolite–kyanite schist—are used to compare the relative importance that geological uncertainty has on bulk compositions determined via (1) X-ray fluorescence (XRF) or (2) point counting techniques. We show that only minor mineralogical variation at the thin-section scale propagates through the phase equilibria modelling procedure and affects the absolute P–T conditions at which key assemblages are stable. Absolute displacements of equilibria can approach ±1 kbar for only a moderate degree of modal proportion uncertainty, thus being essentially similar to the magnitudes reported for analytical uncertainties in conventional thermobarometry. Bulk compositions determined from multiple thin sections of a heterogeneous garnet–staurolite–kyanite schist show a wide range in major-element oxides, owing to notable variation in mineral proportions. Pseudosections constructed for individual point count-derived bulks accurately reproduce this variability on a case-by-case basis, though averaged proportions do not correlate with those calculated at equivalent peak P–T conditions for a whole-rock XRF-derived bulk composition. The main discrepancies relate to varying proportions of matrix phases (primarily mica) relative to porphyroblasts (primarily staurolite and kyanite), indicating that point counting preserves small-scale petrographic features that are otherwise averaged out in XRF analysis of a larger sample. Careful consideration of the size of the equilibration volume, the constituents that comprise the effective bulk composition, and the best technique to employ for its determination based on rock type and petrographic character, offer the best chance to produce trustworthy data from pseudosection analysis

Archean continental crust formed by magma hybridization and voluminous partial melting

Archean (4.0–2.5 Ga) tonalite–trondhjemite–granodiorite (TTG) terranes represent fragments of Earth’s first continents that formed via high-grade metamorphism and partial melting of hydrated basaltic crust. While a range of geodynamic regimes can explain the production of TTG magmas, the processes by which they separated from their source and acquired distinctive geochemical signatures remain uncertain. This limits our understanding of how the continental crust internally differentiates, which in turn controls its potential for long-term stabilization as cratonic nuclei. Here, we show via petrological modeling that hydrous Archean mafic crust metamorphosed in a non-plate tectonic regime produces individual pulses of magma with major-, minor-, and trace-element signatures resembling—but not always matching—natural Archean TTGs. Critically, magma hybridization due to co-mingling and accumulation of multiple melt fractions during ascent through the overlying crust eliminates geochemical discrepancies identified when assuming that TTGs formed via crystallization of discrete melt pulses. We posit that much Archean continental crust is made of hybrid magmas that represent up to ~ 40 vol% of partial melts produced along thermal gradients of 50–100 °C/kbar, characteristic of overthickened mafic Archean crust at the head of a mantle plume, crustal overturns, or lithospheric peels

HT-MP metamorphism in Central Qilian Block, NE Tibet Plateau: implications on the tectonic evolution of the Qilian Orogen

The metamorphism of the Central Qilian Block in the northeastern Tibetan Plateau records a complete tectonic history of the Qilian Orogen. Here, results of structural measurements, geochronological data, and thermobarometry for metamorphic rocks in the Central Qilian Block were presented, which trace the tectonic evolution of the Qilian Orogen. New U–Pb dating of detrital zircon from one paragneiss shows a main age population between c. 1500 and c. 1250 Ma, with the youngest age of 1085 Ma. This unit was intruded by an orthogneiss which has a U–Pb weighted mean age of 920 ± 18 Ma. Together with c. 1200–1000 Ma ages from inherited zircon cores in amphibolites, these results indicate that the protoliths of the Huangyuan Group were formed during the Mesoproterozoic and Neoproterozoic. Rims of these zircons obtain tightly constrained and concordant ages ranging from c. 459 to c. 427 Ma, with a weighted mean age of c. 450 Ma. Phase equilibrium modelling and conventional thermobarometry jointly indicate high-temperature/medium-pressure HT-MP metamorphism along a clockwise pressure–temperature (P–T) path at c. 450 Ma, passing through prograde conditions of 7.8–8.0 kbar and 620–650°C to peak conditions of ~7 kbar and ~780–800°C. Together with documented widespread metamorphism, magmatism, and ductile shear belts, these new results reassert that the Central Qilian Block experienced a three-stage tectono-metamorphic evolution during the Early Palaeozoic and the HT-MP metamorphism of Huangyuan Group was developed by a continental collision geodynamic setting with coeval mafic and granitic magmatism

Intraoceanic subduction system within the Neo-Tethys: evidence from Late Cretaceous arc magmatic rocks of the eastern Himalaya

The tectonic evolution of the Neo-Tethys Ocean remains highly controversial, with several models existing in the community that conflict with each other. Here, we present new geochronologic and geochemical data for orthogneisses and amphibolites from the Greater Himalayan Sequence, eastern Himalayan orogen, which indicate that these rocks have Cenozoic metamorphic ages (∼52–3 Ma), but were derived from Late Cretaceous (∼89 Ma) magmas with arc-like and depleted mantle geochemical signatures. Considering that northern India was a passive continental margin during the Mesozoic, and the previously reported Late Cretaceous magmatic rocks in the eastern Himalaya formed in a continental rifting setting, we suggest that the studied Late Cretaceous arc-type magmatic rocks formed in an intraoceanic arc setting within the Neo-Tethys, and accreted onto the passive margin of the Indian continent prior to the terminal continental collision. When combined with the existence of Late Mesozoic and intraoceanic arc-type magmatic rocks in the western Himalaya, we suggest that a huge Late Cretaceous subduction system operated within the eastern Neo-Tethys Ocean. This study supports two subduction zones having been responsible for the consumption and closure of the Neo-Tethys basin, and a two-stage collision history between India, Asia, and the intermediate island arc system. Our data therefore provide important constraints on the evolution of the Neo-Tethys Ocean and India-Asia collisional orogeny in southern Tibet