Oceanic Subduction Research Articles

Origin of the Songpan–Garzê terrane, Tibetan Plateau: a perspective from the tectonic evolution of the Palaeo-Tethys Ocean

Abstract The Songpan–Garzê terrane is the largest Triassic remnant flysch basin on Earth and formed as the Palaeo-Tethys Ocean closed during the final amalgamation of the Pangaea supercontinent. However, the origin of the Songpan–Garzê terrane is highly controversial. A synthesis of the tectonic evolution of the Palaeo-Tethys Ocean and its branches surrounding the Songpan–Garzê terrane is presented, which clarifies the nature and relationships among the many Palaeo-Tethys sutures. Provenance analyses suggest that branches of the Palaeo-Tethys near the Songpan–Garzê terrane closed before the Early Triassic. In contrast, the main Palaeo-Tethys Ocean (Longmu Co-Shuanghu) did not close until the beginning of the Late Triassic. This study argues against the Songpan–Garzê terrane being a remnant ocean basin, and proposes that it was a back-arc basin of the main Palaeo-Tethys Ocean. It initially underwent extension by the combined effects of the main Palaeo-Tethys Ocean subduction and the Emeishan mantle plume in the Late Permian, and subsequently developed into a back-arc basin in the Triassic, into which huge turbiditie units were deposited derived from all surrounding terranes and orogens. The final closure of the main Palaeo-Tethys Ocean in the beginning of the Late Triassic and subsequent continent–continent collision led to basin inversion in the Late Triassic.

Metamorphic Remnants of the Variscan Orogeny across the Alps and Their Tectonic Significance

Lithospheric slices preserving pre-Alpine metamorphic imprints are widely described in the Alps. The Variscan parageneses recorded in continental, oceanic, and mantle rocks suggest a heterogeneous metamorphic evolution across the Alpine domains. In this contribution, we collect quantitative metamorphic imprints and ages of samples that document Variscan tectonometamorphic evolution from 420 to 290 Ma. Based on age distribution and metamorphic imprint, three main stages can be identified for the Variscan evolution of the Alpine region: Devonian (early Variscan), late Devonian–late Carboniferous (middle Variscan), and late Carboniferous–early Permian (late Variscan). The dominant metamorphic imprint during Devonian times was recorded under eclogite and HP granulite facies conditions in the Helvetic–Dauphinois–Provençal, Penninic, and eastern Austroalpine domains and under Ep-amphibolite facies conditions in the Southalpine domain. These metamorphic conditions correspond to a mean Franciscan-type metamorphic field gradient. During the late Devonian–late Carboniferous period, in the Helvetic–Dauphinois–Provençal and central Austroalpine domains, the dominant metamorphic imprint developed under eclogite and HP granulite facies conditions with a Franciscan field gradient. Amphibolite facies conditions dominated in the Penninic and Southalpine domains and corresponded to a Barrovian-type metamorphic field gradient. At the Carboniferous–Permian transition, the metamorphic imprints mainly developed under amphibolite-LP granulite facies conditions in all domains of the Alps, corresponding to a mean metamorphic field gradient at the transition between Barrovian and Abukuma (Buchan) types. This distribution of the metamorphic imprints suggests a pre-Alpine burial of oceanic and continental crust underneath a continental upper plate, in a scenario of single or multiple oceanic subductions preceding the continental collision. Both scenarios are discussed and revised considering the consistency of collected data and a comparison with numerical models. Finally, the distribution of Devonian to Triassic geothermal gradients agrees with a sequence of events that starts with subduction, continues with continental collision, and ends with the continental thinning announcing the Jurassic oceanization.

Lithospheric-scale dynamics during continental subduction: Evidence from a frozen-in plate interface

Abstract Continental subduction and collision are not merely follow-ups of oceanic subduction but mark the transition from lithospheric-scale deformation localized along the subduction interface to crustal-scale deformation distributed across the orogen. In order to unravel the processes typifying the dynamic changes from oceanic subduction to collision, we have characterized the pressure-temperature (P-T) and spatio-temporal evolution of rocks on each side of the tectonic contact (Briançonnais–Liguro-Piemont [Br-LP] contact) separating the subducted oceanic remnants from the subducted continental fragments along the Western Alps. Results indicate that the maximum temperature and pressure difference on each side of the contact is generally <30 °C and <0.3 GPa, evidencing that no significant metamorphic gap exists. The preservation of similar P-T conditions on both sides of the Br-LP contact is interpreted as resulting from offscraping of the Liguro-Piemont and later Briançonnais units at similar depths, as supported by the ~10 m.y. gap between peak burial ages of both zones. The similar depth range reached by the various units reflects systematic variations of slicing and mechanical coupling along the plate interface suggesting that (1) similar slicing mechanisms and strain localization prevailed during both oceanic and continental subduction and (2) the Br-LP contact represents a frozen-in subduction interface. The end of high-pressure and low-temperature metamorphism and continental subduction at ca. 33 Ma would mark the stalling of subduction interface dynamics and the onset of strain distribution across the plate interface and into the lower plate.

Metamorphic and geochronological constraints on the final assembly of West Gondwana: an Ediacaran-Cambrian collisional setting in the northern Araguaia Orogen, central Brazil

ABSTRACT The Araguaia Orogen, located in Central Brazil, represents an Ediacaran-Cambrian collisional orogen associated with the final assembly of West Gondwana. Two contrasting tectonic models have been proposed for the evolution of this orogen, involving different histories for the Amazonian palaeocontinent drift: (i) subduction of a large ocean and final collision of the Amazonian palaeocontinent with Gondwana between ~ 560 and 500 Ma and (ii) intracontinental orogeny after early amalgamation of the Amazonian palaeocontinent at ca. 620 Ma. Detailed tectono-metamorphic studies in this key locality for the evolution of Gondwana are also still scarce. Hence, in this study we present a detailed metamorphic investigation of two metapelite samples (schists) and one garnet amphibolite aiming to unravel the pressure-temperature-time (P-T-t) paths for the higher grade rocks exposed in this orogen. Pressure-temperature estimates based on thermodynamic modelling using isochemical phase diagrams in the MnNCKFMASHTO system for the schists constrained a single and continuous prograde metamorphic path (burial and heating) from ~ 4–6 kbar and ~ 450–500°C to ~ 7–8 kbar and ~ 580–650°C. The near-peak metamorphic conditions of the garnet amphibolite were estimated from classic geothermobarometry at 7.5 ± 0.8 kbar and 552 ± 40°C. Lutetium-hafnium garnet dating constrained an average age for the prograde metamorphism at around 537–533 Ma. These results show that the Araguaia Orogen records one single metamorphic event during the late Ediacaran-early Cambrian, and likely represents part of a collisional system developed between the Amazonian Craton and Parnaíba Block after the closure of the Clymene Ocean.

Light Mo Isotopes of Post‐Collisional Ultrapotassic Rocks in Southern Tibet Derived From Subducted Indian Continental Crust

AbstractRecycling of molybdenum isotopes in continental subduction zones remains debated. In this contribution, we re‐visit the Mo isotope compositions of the Sailipu post‐collisional ultrapotassic rocks in the Himalaya‐southern Tibet orogen. These ultrapotassic rocks have very varying δ98/95Mo values of −0.66 to −0.07‰ and Mo/Ce ratios of 0.0008–0.005, which are lower than those of mid‐ocean ridge basalts (MORB; δ98/95Mo = −0.20 ± 0.06‰, and Mo/Ce = 0.03) and oceanic subduction‐related (i.e., mantle source involving fluids, residual slab, or oceanic sediments) magmatic rocks (e.g., modern arc lavas, Cenozoic OIB‐type basalts in eastern China and the central Mariana Trough basalts in the back‐arc basin, syn‐collisional andesitic rocks in southern Tibet). Combined with the light Mo isotopes of the Himalayan schists and gneisses, we suggest that the light Mo isotopic signature of the Sailipu ultrapotassic rocks is derived from subducted Indian continental crust. This is consistent with the extremely low δ11B (−17.4 to −9.7‰) and B/Nb (0.16–1) values and enriched Sr‐Nd‐Pb isotopes of the Sailipu ultrapotassic rocks. Thus, this study reveals the recycling of light Mo‐B isotopes during continental subduction and demonstrates that Mo‐B isotopes can effectively distinguish between continental and oceanic subduction.

Tectonic evolution of the Middle-Late Permian orogenic belt in the eastern part of the CAOB: Implications from the magmatism in the Changchun-Kaiyuan area

AbstractVarious magmatisms during the subduction-collision process are crucial to reveal the long-term tectonic evolution of the eastern Central Asian Orogenic Belt. In this paper, we present major and trace elements of whole-rock, zircon U-Pb dating and Hf isotope of the Shanmen pluton. Results imply that the Shanmen pluton consists of quartz diorite and mylonitic granite, with zircon U-Pb ages of 263.7–259.6 Ma. The studied quartz diorite contains high Sr/Y (51.19–90.87) and (La/Yb)N (7.82–13.62) ratios, and belongs to adakitic rocks. Coupled with the positive εHf(t) values of +5.71 to +12.8 with no obvious Eu anomaly, we propose that quartz diorite is the product of the interaction between different degrees of slab melt and the overlying mantle wedge. In contrast, the mylonitic granite has lower MgO (0.28 wt% – 0.47 wt%) contents and positive εHf(t) values of +7.79 to +10.15, indicating an affinity with I-type granite originated by partial melting of the intermediate-basic lower crust. The geochemical characteristics and lithological assemblages, along with the Permian magmatic rocks in the Changchun-Kaiyuan area displaying arc rocks affinity, propose their formation is related to the southward subduction of the Paleo-Asian Ocean (PAO). Based on this study and previous evidence, we lean towards adopting a middle-late Permian slab break-off model, wherein the PAO did not close until the late Permian.

Paleozoic rapid denudation of the Proterozoic Jinchuan Ni-Cu-PGE sulfide deposit, NW China: Insights from single-grain zircon (U-Th)/He thermochronology

The Proterozoic Jinchuan Ni-Cu-PGE sulfide deposit is one of the three largest magmatic sulfide deposits on the earth which was formed in either the Mesoproterozoic or Neoproterozoic. After mineralization of the giant Jinchuan Ni-Cu-PGE orebody, orogenic movements occurred frequently and strongly in a large-scale area. However, less is known on how the deposit denudated after mineralization yet. In this study, we report a new finding of the rapid uplift and denudation event during Late Paleozoic which is supported by single-grain zircon (U-Th)/He dating of 7 samples with elevations ranging between 1453 m and 2180 m collected from a vertical profile in the Jinchuan area. The numerical dating results cluster within 341.08 ± 14.63 – 262.60 ± 29.10 Ma for 7 samples with elevations ranging between 2180 and 1453 m, and 154.74 ± 11.54 Ma for the lherzolite with elevation of 1220 m. Confined with the published 40Ar-39Ar data, the inverse modelling results support that the Jinchuan deposit underwent a rapid cooling at ∼320 Ma. The numerical dating and inverse modelling results obtained in this study support that the Jinchuan deposit and adjacent area experienced a rapid uplift during 341 – 262 Ma (Late Paleozoic). The rapid denudation during the Late Paleozoic is likely a response to the southward subduction of the Paleo-Asian Ocean, inferred from the regional orogenic movements and magmatism in the adjacent areas. Besides, the zircon (U-Th)/He age of 154.74 ± 11.54 Ma obtained from lherzolite supports a cooling event which is to be further studied.

Volcanism Pacing Slumping Gravity Flow Deposits during the Late Carboniferous in the Southern Margin of the Junggar Basin, China

AbstractDeep‐water gravity depositional processes and evolution in arc systems have become topics of intense research focus in recent years. This study discusses the co‐evolution of volcanism and deep‐water gravity flow deposits at the southern margin of the Junggar Basin, based on petrology, geochronology and geochemical analyses. The results show that a massive collapse of unstable sediments from the slope was triggered by volcanism, resulting in the formation of slumping gravity flows. The occurrence of volcanic beds in the slump deposits confirm that synchronous volcanism likely affected sediment instability, triggering gravity flows. The Th/Yb, Ta/Yb and Th/Ta elemental ratios, U‐Pb ages of detrital zircons and paleocurrent directions indicate that the North Tianshan (NTS) island arc represents the provenance of the Qianxia Formation. Moreover, statistical data on the pyroclastic components in the gravity flow deposits reveal an intensity index of volcanism, indicating that volcanism is strongly related to gravity flow deposits, especially in terms of the type and distribution of the deposits. A model for volcanically‐triggered deep‐water gravity flow deposits is established, in order to provide a more in‐depth understanding of the co‐evolution of volcanism and gravity flow deposits within the depositional setting of the late Paleozoic NTS oceanic subduction margin in the Junggar Basin.

Geochemistry and zircon U-Pb-Hf isotopes of Paleozoic granitoids along the Solonker suture zone in Inner Mongolia, China: Constraints on bidirectional subduction and closure timing of the Paleo-Asian Ocean

The closure of the Paleo-Asian Ocean resulted in formation of the Central Asian Orogenic Belt (CAOB) accompanied by considerable Phanerozoic juvenile crustal growth. As an international frontier topic, whether the Solonker suture zone is the site of final closure of the Paleo-Asian Ocean and records the terminal evolution of the CAOB in Inner Mongolia is still in controversial. This study presents the geochemistry and zircon U-Pb-Hf isotopes of granitoids along the Solonker suture zone to constrain the subduction and closure processes of the Paleo-Asian Ocean. Granitoids from the Xilinhot area in the northern accretionary zone comprise diorites and granodiorites with LA-ICPMS zircon U-Pb ages of 316 ± 2 Ma, 302 ± 2 Ma, 285 ± 2 Ma and 243 ± 4 Ma, positive εHf(t) values of +2.39 to +13.31 and TDM2 model ages of 481–1191 Ma. The ca. 316–285 Ma granitoids were formed by fractional crystallization of regional coeval calc-alkaline basaltic magma in an arc-related setting and the ca. 243 Ma granitoid was generated mainly from partial melting of dominant basaltic source rocks and minor metagraywackes in the northward subduction-related setting. The granitoids from the Linxi area in the southern accretionary zone include muscovite monzogranites and monzogranites with LA-ICPMS zircon U-Pb ages of 525 ± 6 Ma, 274 ± 2 Ma and 261 ± 4 Ma. The ca. 525 Ma muscovite monzogranites have negative zircon εHf(t) values of −19.36 to −13.61 with TDM2 model ages of 2348–2711 Ma. They were generated from partial melting of Neoarchean to Paleoproterozoic supracrustal materials. The ca. 274–261 Ma monzogranites show negative zircon εHf(t) values ranging from −10.48 to −0.66 with TDM2 model ages of 1335–1942 Ma. They were derived from remelting of Paleoproterozoic to Mesoproterozoic supracrustal materials in the southward subduction zone along the northern active margin of the North China Craton. Combining our results with existing data, we infer that the early Paleozoic southward subduction of the Paleo-Asian Ocean started before ca. 525 Ma, and the northward subduction zone lasted to ca. 243 Ma, which implies the bidirectional subduction of the ocean lasted at least until ca. 243 Ma. Our work suggests that along the Solonker suture zone, the regional amalgamation may have occurred after ca. 243 Ma, which closed the Paleo-Asian Ocean.

Exploring paleoproterozoic metamorphism in the western Río de la Plata craton from drill-core evidence

Rhyacian basement is not exposed in the western part of the Río de la Plata craton, but deep borehole samples (depth of 1046–3340 m) provide information on its composition and metamorphic evolution. An amphibole schist sample from the Camilo Aldao borehole underwent two metamorphic stages: M1 and M2. The P-T conditions of M1 are not constrained whereas the M2 main metamorphic stage is characterized by conditions of 6.1–7.4 kbar and 575–605 °C. In contrast, an olivine gabbronorite sample from the Santiago Temple borehole records crystallization conditions of 1.7 ± 0.6 kbar and 1180 ± 10 °C based on clinopyroxene composition, and late magmatic to subsolidus low-pressure granulite facies metamorphism at 2.3 ± 1.2 kbar and 749 ± 85 °C based on multi-equilibrium thermobarometry. The Rhyacian evolution of the Río de la Plata craton is characterized by at least three thermo-tectonic orogenic stages: T1 (ca. 2200–2100 Ma), T2 (ca. 2100–2040 Ma), and T3 (ca. 2040–2010 Ma) that correspond respectively to: 1) early oceanic subduction and juvenile continental arc magmatism; 2) subsequent continental collision (clockwise intermediate P/T paths of metamorphism), and 3) final low-pressure metamorphism coeval with post-orogenic magmatism. The Camilo Aldao amphibole schist is representative of the T2 orogenic stage while the Santiago Temple gabbronorite is representative of T3 magmatism and metamorphism.

The competition between anthropogenic aerosol and greenhouse gas climate forcing is revealed by North Pacific water-mass changes.

Modeled water-mass changes in the North Pacific thermocline, both in the subsurface and at the surface, reveal the impact of the competition between anthropogenic aerosols (AAs) and greenhouse gases (GHGs) over the past 6 decades. The AA effect overwhelms the GHG effect during 1950-1985 in driving salinity changes on density surfaces, while after 1985 the GHG effect dominates. These subsurface water-mass changes are traced back to changes at the surface, of which ~70% stems from the migration of density surface outcrops, equatorward due to regional cooling by AAs and subsequent poleward due to warming by GHGs. Ocean subduction connects these surface outcrop changes to the main thermocline. Both observations and models reveal this transition in climate forcing around 1985 and highlight the important role of AA climate forcing on our oceans' water masses.

Mesozoic–Cenozoic cooling, exhumation and tectonic implications of Chaqiabeishan–Shaliuquan Li[sbnd]Be ore district in the northeastern Qinghai–Tibet Plateau

The Mesozoic-Cenozoic exhumation and deformation of the northeastern Qinghai–Tibetan Plateau generated extensive regional rugged mountains. We aimed to explore the prolonged kinematic evolutionary history since the Mesozoic. We collected Mesozoic magmatic rock samples from the Chaqiabeishan-Shaliuquan lithium‑beryllium ore district in the Northern Qaidam tectonic belt (NQTB) to constrain the cooling exhumation history using apatite fission track (AFT) thermochronology. Twenty samples yielded AFT ages from 111 to 11 Ma, including four age groups:111–103, 89–75, 71–47, and 30–11 Ma, with mean track lengths of 12.5–11.2 μm. Four thermal history phases of 120–90, 90–50, 50–20, and 20–0 Ma were established based on the inverse thermal model and isostatic adjustment of the lithosphere. Additionally, the average cooling rates of the four phases were 1.67 ± 0.17, 1.25 ± 0.13, 0.33, and 1.75 ± 0.5 °C/Myr, with a corresponding exhumation amount of 1.15 ± 0.07 km, 0.98 ± 0.06 km, 0.86 ± 0.05 km, and 1.09 ± 0.06 km respectively, and the maximum exhumation was 6.021–6.787 km since 120 Ma. In the first phase, the far-field effects of Mid-Tethys Ocean subduction and closure drove the rapid exhumation of the research area. The closure of the Neo-Tethys Oceans and the continuing collision of the India-Asia plate disturbed the northeastern plateau, inducing regional backlash fault events in the NQTB, demonstrated in the second phase of thermal history. The third phase of the stable tectonic setting is consistent with previous research results, suggesting a similar cooling history for the NQTB and Qilian Mountains, evidenced by the widespread absence of Eocene to Miocene cooling signals. The rapid exhumation and reactivation of the NQTB and Qilian thrust fault system since 20 Ma represents the final stage of growth and expansion of the northeastern Qinghai–Tibet Plateau.

Mo isotopes record recycling of anoxic sediment in a Paleo-oceanic subduction zone

Mafic arc igneous rocks with oceanic crust-derived fluids in their mantle sources generally show heavy Mo isotope compositions. However, these rocks may exhibit variable Mo isotope compositions due to involvement of seafloor sediments under different redox conditions. Thus, it is intriguing how Mo isotope recycles with sediment subduction and how Mo isotope fractionation occurs during this process. Here we report Mo isotope data for Paleozoic mafic arc igneous rocks from the Qinling-Tongbai orogenic belt in Central China. The Qinling mafic rocks have δ98Mo values of −0.03 to 0.40‰, and the Tongbai mafic rocks have δ98Mo values of −0.26 to 0.17‰. These values are considerably higher than the normal mantle values of −0.20 ± 0.01‰. In combination with low Mo/Ce ratios, enriched SrNd isotope compositions, and high melt-mobile incompatible trace element contents for these rocks, the high δ98Mo values suggest that anoxic sediment-derived melts with heavy Mo isotopes were incorporated by oceanic subduction into the mantle sources. The high δ98Mo values of the mafic rocks indicate the limited Mo isotope fractionation during subduction of the anoxic sediment. This inference is further verified by modeling of Mo isotope fractionation under reducing conditions in the subduction zone. In this regard, Mo isotopes can be used to identify the recycling of sediments in different types and thus provide geochemical constraints on the redox conditions of oceanic subduction zones.

Circum-cratonic mantle archives the cumulative effects of plume and convergence events

The cumulative effects of heating, various melts/fluids, and tectonic stress from plume and convergent events can produce complex geochemical signatures in circum-cratonic mantle domains. The lithospheric mantle beneath the northwestern Tarim Craton witnessed the Tarim plume and multi-stage convergent events, thus providing a good opportunity to investigate the depletion and enrichment of circum-craton mantle driven by their cumulative effects. Here, we present the in-situ mineral major/trace elements, clinopyroxene Sr isotopic and whole-rock analyses of platinum group elements and Re-Os isotopes for peridotite xenoliths from Piqiang in the northwestern Tarim Craton. Combined with peridotite xenoliths from Tuoyun (South Tianshan Belt) and Xikeer (inner Tarim), we established a lithospheric mantle profile to constrain the formation and subsequent evolution of the circum-Tarim mantle in response to the subduction of the archipelagic ocean and the impingement of the Tarim plume. The Piqiang peridotites have moderately refractory compositions (Mg# in Ol: 90.1–90.9; average values) and can be divided into Group A lherzolites (subgroups A1, A2 and A3) and Group B wehrlite. Group A1 lherzolites contain copper-poor sulfides and clinopyroxenes depleted in incompatible elements, and are interpreted as melting residua that experienced negligible metasomatism. Group A2 lherzolites have disequilibrium textures and are depleted in Nb and Ta with high Ti/Eu but low (La/Yb)N ratios in clinopyroxene, demonstrating recent volatile-rich melt metasomatism. Carbonatitic metasomatism from a mixed crust-mantle source is recorded in the Group A3 lherzolites, as their clinopyroxenes have highly variable 87Sr/86Sr ratios (0.7028–0.7057) with low Ti/Eu and high (La/Yb)N. Overall, Group A lherzolites display high initial 187Os/188Os ratios (0.1198–0.1271, excluding suprachondritic values), yielding the Neoproterozoic peak (∼650 Ma) TRD eruption model ages and the minimum estimate for the age of lithospheric mantle formation of 1.1 Ga. Their high-Na clinopyroxenes, low-Ca orthopyroxenes, and low melting pressures (<3 GPa) suggest that the initial melt extraction of these peridotites was the result of lateral mantle accretion triggered by multi-stage subduction of the Neoproterozoic Pan-Rodina Ocean. In this case, the greater fertility of the Tuoyun peridotites compared to those from Piqiang defines a trend of lateral accretion, while the low whole-rock Al2O3 but high 187Os/188Os ratios of the Xikeer peridotites point to a second melt extraction event. This melt extraction event is interpreted as reflecting low-flux plume melts that heated and melted the overlying lithospheric mantle at the edge of the early Permian Tarim plume. Meanwhile, the Piqiang peridotites, further from the center of the Tarim plume, underwent carbonatitic metasomatism originating from the lower-flux plume melts that involved the high-density components from recycled oceanic crustal materials. The mantle melt activity and faulting related to the India-Eurasia collision provided the reactive melts to form the Group B wehrlite and prepared migration pathways for volatile-rich melts to metasomatize Group A2 lherzolites. Collectively, the accretion, re-extraction, and metasomatic modification of circum-craton mantle reflected a response to the cumulative effects of the plume and convergence events. Our study reveals that the mechanical, thermal, and chemical effects from different dynamic events interacted with each another and worked together to affect the circum-craton mantle evolution.

The geochemistry and geochronology of garnet from the Jinchanghe Fe-Cu-Pb-Zn deposit, Baoshan Block, western Yunnan

The Jinchanghe polymetallic Fe-Cu-Pb-Zn deposit is located in the northern region of the Baoshan Block in western Yunnan. Using electron probe micro-analysis and laser ablation inductively coupled plasma mass spectrometry, we analyzed the major, trace, and rare earth elemental contents of garnets from the Jinchanghe deposit. Moreover, U-Pb dating was performed using the laser ablation sector field inductively coupled plasma mass spectrometry (LA-SF-ICP-MS) method to determine the direct mineralization age of the Jinchanghe deposit. Petrographic analyses revealed three garnet generations: early (GrtI), middle (GrtII), and late (GrtIII). GrtI (And55.69-64.47Gro30.28-41.21Ura + Pyr + Spe + Alm2.35-5.49) was characterized by irregular shapes and weak optical properties, and primarily comprised Fe-rich grossular. GrtII (And53.69-96.17Gro2.78-45.07Ura + Pyr + Spe + Alm0.59-2.12) was characterized by andradite, in which GrtIIa was a homogeneous body, whereas GrtIIb shows exhibited oscillatory zoning with abnormal interference colors. Fine-veined andradite aggregates were predominant in GrtIII (And70.51-90.33Gro6.34-27.80Ura + Pyr + Spe + Alm1.03-3.32). The GrtI, GrtIIa, and the rims of GrtIIb were depleted in heavy rare earth elements and enriched in light rare earth elements, with positive Eu anomalies, indicating that these garnets crystallized from mildly acidic magmatic hydrothermal fluids. However, the cores of GrtIIb had negative Eu anomalies and may have crystallized from a relatively neutral Cl-rich fluid. GrtIII was severely depleted in light rare earth elements with positive Eu anomalies, indicating the presence of a weakly acidic environment. Accordingly, owing to a decrease in temperature and an increase in oxygen fugacity, the pH of the Jinchanghe deposit may have transitioned from an early weakly acidic to neutral, then to a late weakly acidic environment. Furthermore, the ages of GrtII and GrtIII ranged from 520 to 496 Ma, which places the formation of the Jinchanghe deposit during the Late Cambrian. The deposit is closely related to the arc magmatic rocks produced by the subduction of the Proto-Tethys Ocean.

Coexisting Nb-enriched basalts and arc volcanic rocks in the northern Qiangtang Terrane, China: Implications for the effects of ambient mantle on subduction zone magmatism

The petrogenesis of Nb-enriched basalts (NEBs) is controversial. NEBs are generally considered to form by melting of either a mixed enriched and depleted mantle source or a mantle wedge metasomatized by adakitic melts. Here we present geochronological, petrological, geochemical, and Sr-Nd-Os isotope data for Early Permian volcanic rocks from the Tuotuohe area in the northern Qiangtang Terrane, China. The studied volcanic rocks can be divided into two groups: arc basalts and basaltic andesites in the lower sequence and NEBs in the upper sequence. Zircon grains from one basaltic andesite sample in the lower sequence yielded a concordia age of 302.0 ± 1.1 Ma. The geochemical and isotopic characteristics suggest that the NEBs were derived from a mixed source consisting of enriched mid-oceanic ridge basalts (MORB)- to oceanic-island basalt-like ambient mantle and slab-derived fluids, whereas the arc basalts and basaltic andesites originated from a source made of the MORB-type ambient mantle modified by slab-derived melts. We propose for the first time that the Nb enrichment of the NEBs was most likely inherited from high-Nb ambient mantle and is independent of the addition of slab-derived components. In addition to slab-derived components, the nature of the ambient mantle is also an important factor controlling the diversity of arc magmatism. Our new findings show that the Late Carboniferous−Early Permian magmatism in the northern Qiangtang Terrane was related to northward subduction of the Paleo-Tethys Ocean, suggesting that subduction might have occurred in at least the Early Permian rather than in the previously thought Middle Permian.

Source and Migration of Fluids in a Meso-Tethyan Subduction Zone: Fluid Inclusion Study of Syn-Mélange Veins from the Mugagangri Accretionary Complex

The Mugagangri Group (MG), located at the southern margin of the Qiangtang terrane in Tibet, is a crucial research target for understanding the subduction and accretion history of the Meso-Tethys Ocean. Extensional crack-seal veins restricted within sandstone blocks from the broken formation in the MG (Gaize) formed synchronously in the mélange formation. The primary inclusions trapped in the veins recorded multiple pieces of information during the formation of the accretionary wedge. To precisely constrain the MG subduction–accretion processes, we investigated the trapping temperature, salinity, density, and composition of the fluid inclusions within the crack-seal veins derived from the broken formation in the MG (Gaize). The primary inclusions indicate that the crack was sealed at ~151–178 °C. The salinity of the primary inclusions exhibited a well-defined average of 3.3 ± 0.7 wt% NaCl equivalent, slightly lower than the average of seawater (3.5 wt%). There were no nonpolar gases, and only H2O (low salinity) was detectable in the primary inclusions. These characteristics suggest that the syn-mélange fluids were a type of pore fluid in the shallow subduction zone, with the principal source being pore water from sediments overlying the oceanic crust. Because of mineral dehydration and compaction, the pore fluids became more diluted with H2O and fluid overpressure owing to a pore fluid pressure that was greater than the hydrostatic pressure. Subsequently, the creation of cracks through hydraulic fracturing provided a novel pathway for the flow of fluids which, in turn, contributed to the décollement step-down and underthrusting processes. These fractures acted as conduits for fluid movement and played a crucial role in facilitating these peculiar occurrences of quartz veins. The depth (~5 km) and temperature estimates of the fluid expulsion align with the conditions of the décollement step-down, thereby leading to the trapping of fluids within the sandstone blocks and their subsequent underplating to the accretionary complex. In our preferred model, such syn-mélange fluids have the potential to provide valuable constraints on the subduction–accretion processes occurring in other accretionary complexes.

Early Paleozoic oblique convergence from subduction to collision: Insights from timing and structural style of the transpressional dextral shear zone in the Qilian orogen, northern Tibet of China

Transpressional shear zones commonly occur in ancient and modern convergent plate boundaries to accommodate oblique plate convergence. The early Paleozoic Qilian orogen in northeastern Tibet records the subduction of Proto-Tethyan Ocean lithosphere and the accretion-collision of various magmatic arcs and continental terranes. This study focused on the Datong ductile shear zone, which represents the central part of the WNW-ESE−striking ductile shear zone along the northern margin of the Qilian block in the Qilian orogen. This structure bears key information about the evolution of oblique convergence during the early Paleozoic orogeny. The kinematics and timing of the Datong ductile shear zone were investigated via field-based, microstructural, and mica 40Ar/39Ar dating analyses. Mesostructural and microstructural data showed predominantly dextral strike-slip shearing within the Datong ductile shear zone. Microstructural features and quartz c-axis crystallographic preferred orientation patterns indicated that dextral ductile shearing occurred under lower-amphibolite-facies conditions (∼500−550 °C and ∼5.6 kbar) within the shear zone. Microstructures of quartz showed subgrain rotation (SGR) and grain boundary migration (GBM), suggesting dislocation creep−dominated deformation. A strain rate of 10−12 s−1 and a differential stress of 25−39 MPa were estimated by the rheological flow law and quartz paleopiezometry. Finite strain measurements indicated that all deformed rocks of the Datong ductile shear zone exhibit a weakly oblate ellipsoid near the plane strain. Kinematic vorticity (ranging 0.47−0.83) analysis suggested the coexistence of simple shear and pure shear strains within the Datong ductile shear zone, indicating a transpressional setting. Biotite and muscovite 40Ar/39Ar data showed that transpressional shearing deformation started in the Ordovician (before 453 Ma) and lasted to the Silurian (ca. 430 Ma). Our new data combined with regional geological data show that the deformation type, kinematics, and dynamics of the Datong ductile shear zone were controlled by the southward oblique subduction of the Paleo-Qilian Ocean (Proto-Tethyan Ocean) and the following oblique collision between the Qilian block and the Alxa block. The intensive transpressional deformation along the northern Qilian block may reflect strong coupling between the subducting Paleo-Qilian oceanic slab and the overriding Qilian block as well as a high degree of convergence obliquity during the ongoing early Paleozoic convergence.

Exhumation of ultrahigh-pressure eclogite and metarodingite within UHP serpentinites from the Chinese southwestern Tianshan

Extensive research over the past two decades suggests that serpentinite loses its ability to effectively facilitate the exhumation of high-density high-pressure (HP) eclogites at depths of 60–80 km under average subduction zone geotherms. However, at Changawuzi, the western termination of Chinese southwestern Tianshan, eclogite and eclogitic metarodingite are commonly enclosed within ultrahigh-pressure Ti-chondrodite serpentinites (510–530 °C, 37 ± 7 kbar). This metamorphic mafic–ultramafic association, which has undergone a complex evolution from oceanic alteration to deep subduction and subsequent exhumation, presents a unique opportunity to study the incompletely understood mechanisms leading to the exhumation of ultrahigh-pressure (UHP) eclogite in UHP serpentinite.The Changawuzi serpentinite wraps three different rock types: retrograded eclogite, Rodingite I and Rodingite II. Field observations and petrological analyses indicate that the precursor of Rodingite I intruded as vein into the serpentinite protolith prior to subduction, sharing the evolution towards ultrahigh-pressure conditions. According to its composition, the serpentinite originated from a depleted oceanic mantle peridotite, with only minor enrichment in fluid-mobile elements (FME) during oceanic serpentinization. The eclogites have N-MORB like bulk compositions, with lower Si but higher Mg content compared to other eclogites in the Tianshan area, and appear to have been relatively unaffected by oceanic alteration processes. The Changawuzi serpentinite and eclogites likely evolved within a slab setting with limited interaction with seawater. Petrological studies and phase equilibria modeling suggest that the eclogites experienced peak UHP metamorphism at 28–34 kbar and around 460–465 °C. Density measurements yielded relatively low densities of 2.66–2.97 g/cm3 for the host serpentinites, but 3.27–3.63 g/cm3 for Rodingite I. The similarity in the low temperature/ultrahigh pressure conditions, as well as the geochemical characteristics shared by the serpentinite and eclogites, further support the connection between the Changawuzi serpentinites and the interior of the subducting slab. According to our data, two types of UHP eclogitic mafic rocks exhumated with UHP serpentinite. (1) High density Rodingite I veins are embedded within serpentinite and share a common evolution and exhumation history with the host serpentinite. (2) The eclogite and serpentinite tectonically detached from the interior of the subducting slab at different depths and ascended to the plate interface above the slab (subduction channel). These two rock types experienced separate exhumation processes within the subduction channel during the early stages of exhumation. Subsequently, the buoyant serpentinite captured the denser retrograded eclogite at a depth of approximately 80 km, and together they underwent further retrograde metamorphism as they ascended towards the surface. During a later exhumation stage, external fluid induced a second serpentinization process, leading to the formation of Rodingite II within the subduction channel. This study presents a novel exhumation mechanism specific to the Tianshan eclogites and highlights the significant role of deeply subducted serpentinite in the exhumation of eclogites and eclogitic metarodingites from sub-arc depths within an oceanic subduction zone.

Zircon and apatite fingerprints on the magma fertility for the Jurassic igneous rocks in the southern Lhasa subterrane, Tibet

The products of long-lasting magmatism (ca. 170 Myr) related to the northward subduction of the Neo-Tethys Ocean are preserved in the southern Lhasa subterrane (SLT). However, only a few large porphyry Cu–Au deposits related to these igneous rocks have been discovered in the Xiongcun ore concentration area. Here, we present zircon and apatite data on Jurassic fertile suites (FS), barren suites (BS), and recently proposed prospecting target suites (PTS) in the SLT, intending to present a case study of the regional-scale porphyry metallogenic potential. The zircon geochemical data with Rayleigh fractionation modeling show that the zircon geochemical characteristics of the FS and PTS were affected by high proportions of amphibole and titanite crystallization, whereas the zircon geochemistry of the BS displays the effects of plagioclase fractionation or mixing with the least evolved melts. The zircon Ce4+/Ce3+ and Eu/Eu* values cannot be used to trace the oxygen fugacity of the magma due to prior or simultaneous mineral crystallization. However, the zircon Ce/Ce* and Ce/Nd ratios may be used as semiquantitative measures of the magma redox conditions. This result is similar to that calculated by the new magmatic oxybarometer using trace elements in zircon, showing that the FS (ΔFMQ = 1.21 ± 0.42) and PTS (ΔFMQ = 1.29 ± 0.51) have higher oxygen fugacity levels than the BS (ΔFMQ = 0.40 ± 0.88). The melt S contents estimated from magmatic apatites are lower in the FS and PTS than in the BS, which is attributed to prior anhydrite crystallization under high oxygen fugacity and low temperature conditions. The higher Cl in the melt of the FS (average 0.098 wt%) than in the melts of the PTS (average 0.045 wt%) and BS (average 0.04 wt%) shows that elevated Cl contents in melts are likely important for forming porphyry Cu–Au deposits. The combined geochemical characteristics of zircon and apatite suggest that a high oxygen fugacity (ΔFMQ > 1.2), prior or concurrent anhydrite crystallization, and high zircon Eu/Eu* (>0.6), Ce4+/Ce3+, and (Yb/Dy)N values can be used to refine potential porphyry Cu–Au targets on a regional scale in the SLT. High melt Cl contents and apatite MnO and Pb contents can effectively indicate fertile suites at the deposit scale.