Evolution Of Subduction Research Articles

Late Paleozoic–Early Mesozoic subduction history of Mudanjiang Ocean: New insights from greenschist in the Heilongjiang Complex, NE China

This paper presents zircon U-Pb ages, whole-rock major and trace element data, and Sr-Nd-Hf isotopic compositions for the greenschists from the Heilongjiang Complex in NE China. Based on the geochronological and geochemical characteristics, the greenschists from the Heilongjiang Complex can be categorized into three series: Early Permian tholeiitic series (275 ± 1 Ma and 273 ± 1 Ma), Middle Permian alkaline series (261 ± 1 Ma), and Early Jurassic alkaline series (179 ± 1 Ma). The Early Permian tholeiitic greenschist samples are characterized by remarkable negative Nb and Ta anomalies, and they exhibit relatively low whole-rock (87Sr/86Sr)i values (0.703896–0.704252), along with positive εNd(t) values (+7.66 to +8.65) and positive εHf(t) values (+0.92 to +11.51), implying that their protoliths originated from the partial melting of a depleted mantle wedge that has been metasomatized by subducted slab fluids. The Middle Permian alkaline greenschist samples have high Nb (17.3–19.0 ppm), TiO2 (1.40–1.64 wt%), and P2O5 (0.60–0.68 wt%) concentrations, and they are enriched in large ion lithophile elements and light rare earth elements, and depleted in high field strength elements and heavy rare earth elements, with high whole-rock (87Sr/86Sr)i values (0.709434–0.709949), negative εNd(t) values (−4.73 to −4.72) and varied εHf(t) values (−10.11 to +4.86), which compositionally align with those of Nb-enriched basalts. Based on the geochronological and geochemical characteristics, we propose that the protolith generation of Middle Permian alkaline greenschists was associated with the partial melting of a mantle wedge that had been metasomatized by subduction slab-derived melts. The Early Jurassic alkaline greenschist samples display OIB affinities, characterized by whole-rock (87Sr/86Sr)i values ranging from 0.705161 to 0.705926, εNd(t) values between +6.16 and +6.19, and εHf(t) values from +3.44 to +9.32, indicating the formation of the protolith for the Early Jurassic alkaline greenschist was related to the upwelling of the asthenosphere resulting from slab break-off. By integrating regional geological data, we can trace the subduction evolution of the Mudanjiang Ocean back to the Early Permian, and the slab break-off during the Early Jurassic likely corresponded to the early stage of the closure of the Mudanjiang Ocean.

The spatial–temporal evolution of Neotethyan subduction in central to Southeast Iran: Constraints from geochemistry, zircon U–Pb–Hf, and Sr–Nd–Pb isotopes

Igneous rocks in central Iran preserve proof of several episodes of Late Mesozoic–Cenozoic Neotethyan subduction, and record temporal and spatial variations from subduction to collisional geodynamic regimes that remain a subject of debate. We revisit magmatic rocks from the central Urumieh–Dokhtar magmatic arc (UDMA; 32° 30′ N to 36° 00′ N) and report new major and trace element analyses, whole-rock Sr–Nd–Pb-isotopic data and zircon U–Pb–Hf ages on plutonic rocks from the Kajan region. The data reveal magmatic pulses between 32 and 22 Ma that are represented by igneous rocks with geochemical and isotopic features resulting from melting of a metasomatized, enriched mantle. Modeling indicates the samples represent a mixture of 97% batch melt (from a source with 98.5% mantle melt +1.5% terrigenous sediment) with 3% upper continental crust. The oldest yielded age (i.e., 32 Ma) rules out a lull in magmatism in the central UDMA between 37 and 26 Ma and speaks against distinct flare-up magmatic episodes at ∼54–37 Ma and ∼ 20–5 Ma in the frontal arc. We hypothesize that magmatism records down-going slab tearing into two pieces beneath the central to SE UDMA, which then changes the geometry and buoyancy of subducting Neotethyan lithosphere. The down-going slab remained involved in central UDMA (i.e., rollback) during Early Miocene whilst arc retreat beneath the SE UDMA proceeded with detachment of subducted oceanic lithosphere (i.e., break-off). We hypothesize a diachronous collision in the UDMA, first in the southeastern segments, and subsequently in the central UDMA. These findings bear important implications for the geodynamic evolution of the Zagros orogen and run counter to a model suggesting that the collision initiated in the northwest to central and propagated progressively to the southeast along the Zagros suture zone.

Secular Changes in the Occurrence of Subduction During the Archean

AbstractSubduction processes play a pivotal role in facilitating material exchange between the crust and mantle, contributing to the growth of continents. However, the onset and evolution of subduction remain hotly debated. Here, we developed a high‐dimensional machine learning (ML) model to use multiple compositional data (e.g., Nb/La, Nb, Ti, Nb/U, Pb/Nd, and Nb/Th) to distinguish arc‐type from non‐arc basalts worldwide, then applied this well‐trained ML model to identify and delineate the secular occurrence of Archean arc‐type basalts. Our findings indicate that basaltic arc magmatism can be traced back to only a few early Archean cratons. A noticeable increase in the prevalence of arc‐type basalts during the Mesoarchean and Neoarchean suggests the synchronous formation of Archean subduction across different cratons. This observation hints at transitioning from the predominantly stagnant geodynamic regime to the early, more mobile tectonic regime.

Inferring the Paleo‐Location of Proto‐Arc Magmas During Subduction Infancy in the Izu‐Bonin‐Mariana

AbstractSubduction zones are one of the principal drivers of the modern‐day plate tectonics, but the processes that develop as incipient subduction zones mature are yet to be understood. Finding modern analogs for different stages of subduction infancy remains an outstanding challenge to answer questions on how and when an oceanic subduction zone reaches a quasi‐steady state to become long‐lived. Here, we compare the southern Mariana intra‐oceanic arc, in which near‐trench spreading and infant arc magmatism developed ∼3–4 Myr ago, with the Izu‐Bonin‐Mariana (IBM) proto‐arc magmas that formed during subduction infancy ∼52 Myr ago. Building upon this comparison, we propose that the Southern Mariana may be considered as a modern analog to subduction infancy. The similarities between the Southern Mariana arc magmas and the IBM Eocene proto‐arc magmas suggest that the IBM proto‐arc crust might have formed within 80–90 km from the paleo‐trench (slab at < 90 km paleo‐depth), while the Eocene infant arc was likely located at ∼100 km from the paleo‐trench (∼100–125 km slab paleo‐depth). We use our constraints further to propose a new conceptual model of subduction evolution for IBM. In this model, development of the subduction channel during arc infancy facilitated downward slab penetration and intensified corner flow in the sub‐arc asthenosphere of the mantle wedge. As such, cooling and serpentinization of the fore‐arc mantle might be considered as the principal drivers of subduction maturation and stabilization over the IBM lifetime.

An Overview on the Composition and Age of Upper Crust of Proto‐Tethyan Lajishan Intra‐oceanic Arc, NE Tibet Plateau

AbstractIdentification and anatomy of oceanic arcs within ancient orogenic belt are significant for better understanding the tectonic framework and closure process of paleo‐ocean basin. This article summarizes the geological, geochemical, and geochronological characteristics of upper crust of Proto‐Tethyan Lajishan intra‐oceanic arc and provides new data to constrain the subduction evolution of the South Qilian Ocean. The intra‐oceanic arc volcanic rocks, including intermediate–mafic lava, breccia, tuff, and minor felsic rocks, are distributed along southern part of the Lajishan ophiolite belt. Geochemical and isotopic compositions indicate that the intermediate–mafic lava were originated from depleted mantle contaminated by sediment melts or hydrous fluids, whereas the felsic rocks were likely generated by partial melting of juvenile mafic crust in intra‐oceanic arc setting. Zircons from felsic rocks yield consistent and concordant ages ranging from 506 to 523 Ma, suggesting these volcanic rocks represent the relicts of upper crust of the Cambrian intra‐oceanic arc. Combined with the Cambrian forearc ophiolite and accretionary complex, we suggest that the Cambrian intra‐oceanic arc in the Lajishan ophiolite belt is belonging to the intra‐oceanic arc system which was generated by south‐directed subduction in the South Qilian Ocean at a relatively short interval between approximately 530 and 480 Ma.

Horizontally forced initiation of the Izu-Bonin-Mariana subduction zone

The sparsity of a direct record for the moment of subduction zone initiation has led to various models describing the infancy and evolution of modern oceanic subduction systems. Recently, with increases in available samples and geochemical data for subduction zone initiation-to-mature-arc lavas, better constraints on subduction evolution are possible. Here, by systemically modeling the time-space pattern and geochemical characters of forearc magmas with forward numerical modeling, we attempt to search for a best-fit geodynamic scenario where Izu-Bonin-Mariana-type subduction tends to develop. Our modeling and geochemical constraints have identified a necessary and possibly transitory pre-subduction zone initiation trenchward contraction consistent with observed Izu-Bonin-Mariana forearc magma geochemistry. Our results also reveal a typical maturation process for Izu-Bonin-Mariana-type oceanic subductions, controlled by the pace of the upper plate’s rifting and solidification.

Origin of Early Cretaceous mafic volcanic rocks from the Erlian Basin west of the Great Xing’an Range of North China: Implications for the tectono-magmatic evolution of East Asia

Early Cretaceous intraplate volcanic rocks are widespread in NE Asia, but their origin remains controversial. This work presents zircon U-Pb ages, whole-rock element and Sr-Nd isotope data for mafic volcanic rocks from the Erlian Basin, a wide rift basin in NE Asia. There were two episodes of Early Cretaceous mafic volcanism in the Erlian Basin, and the eruptions show contrasting geochemical compositions. The early mafic volcanic rocks, with U-Pb ages of ca. 140−135 Ma, show slightly depleted Sr-Nd isotope compositions (ISr(t) = 0.7042−0.7052; εNd(t) = +0.82 to +3.0) and arc-like trace-element compositions, which are derived from subduction-related fluid/melt metasomatized lithosphere mantle. The late mafic volcanic rocks (dated at ca. 125 Ma) have enriched Sr-Nd isotopes (ISr(t) = 0.7055−0.7077; εNd(t) = −0.50 to −2.67) and oceanic-island basalt (OIB)-like trace-element compositions, revealing the metasomatism of melts from crustal materials and asthenosphere mantle. The two types of mafic volcanic rocks may record the interactions of the mantle and melts from the subducted paleo-Pacific oceanic slab at different depths. The landward-then-oceanward migration pattern of the Mesozoic volcanism from NE Asia can be explained by the flat subduction and subsequent slab roll-back of the Paleo-Pacific Ocean, consistent with migration patterns from the North China Craton and South China Block, implying similar Jurassic−Cretaceous subduction evolution along the entire East Asia margin. Some Late Jurassic to Early Cretaceous dates from east Mongolia and the southern margin of the Erlian Basin diverge from this trajectory. In combination with previous studies, we suggest that the Early Cretaceous pervasive intraplate volcanism in the Erlian Basin and adjacent areas of NE Asia mainly resulted from the slab roll-back of the Paleo-Pacific Ocean with a combined effect from the post-collision extension of the Mongol-Okhotsk orogen.

Continental subduction of Adria in the Apennines and relation with seismicity and hazard

The subduction of continental lithosphere is a complex process because the buoyancy of the crust is higher than the oceanic and should resist sinking into the mantle. Anyway, studies on the Alpine-Himalayan collision system indicate that a large portion of the continental crust is subducted, while some material is accreted in the orogens. The Apennine is a perfect case for studying how such processes evolve, thanks to high quality seismic images that illuminate a critical depth range not commonly resolved in many collisional settings. In this paper, we show the structure of the Apennines orogen, as jointly revealed by seismicity and deep structure from regional and teleseismic tomography and receiver function profiles. The westward subducting Adria lithosphere is well defined along the orogen showing a mid-crustal delamination. Seismicity within the underthrusting lower crust and velocity anomalies in the mantle wedge highlight how the subduction evolution is entangled with the liberation of fluids. The eclogitization of subducted material enhances the fluid release into the wedge, the delamination and retreat of the Adria plate. This delamination/subduction generates a coupled compression and extension system that migrates eastward following the retreat of the lithosphere, with broad sets of normal faults that invert or interfere with pre-existing compressional structures all over the roof plate. The sparseness and non-ubiquity of intermediate depth earthquakes along the subduction panel suggest that the brittle response of the subducting crust is governed by its different composition and fluid content. Therefore, the lower crust composition appears essential in conditioning the evolution of continental subduction.

Pre‐Subduction Architecture Controls Coherent Underplating During Subduction and Exhumation (Nevado‐Filábride Complex, Southern Spain)

AbstractThe interplay between structural and metamorphic processes operating along the deep plate interface in subduction zones remains elusive as much of the geologic record is recycled into the mantle. In some cases, metamorphosed subducted rocks are underplated and exhumed to the surface, providing critical constraints on structural processes and the rheological evolution of subduction interfaces at convergent margins. One such exhumed high‐pressure/low‐temperature subduction complex is the Cenozoic Nevado‐Filábride Complex (NFC) in Southern Spain. This study presents new data from the NFC that elucidate the syn‐metamorphic deformation, stacking, and underplating of continental slivers along the subduction interface. The structurally lowest NFC dominantly comprises lithologically monotonous Paleozoic metamorphic basement rocks recorded by apatite U‐Pb ages and shows no evidence for large‐scale internal duplications suggesting it behaved as a coherent basement succession during subduction. In contrast, structurally higher levels of the NFC are characterized by the stacking of older‐on younger coherent slices and distinctly different metamorphic ages. These relationships document syn‐subduction structural repetitions and tectonic stacking of imbricate thin slivers (∼100s m) during subduction underplating. Structurally higher levels of the NFC exhibit both Eocene and Miocene metamorphic zircon rims and apatite ages, along with microstructures indicative of relatively higher temperature metamorphism. Large‐scale underplating and antiformal stacking of slivers in the subduction channel can provide buoyancy forces to underplate and assist exhumation. We demonstrate that the presubduction stratigraphic architecture is a key control on the style and timing of deformation and metamorphism, facilitating coherent subduction underplating.

Tectonic evolution of circum-Rodinia subduction: Evidence from Neoproterozoic A-type granitic magmatism in the Central Tianshan Block, northwest China

Neoproterozoic igneous rocks associated with the assembly and configuration of the Rodinia supercontinent are widely distributed in the Central Tianshan Block, southwestern Central Asian Orogenic Belt (CAOB), northwest China. The Central Tianshan Block is an important geological unit in tectonic reconstructions of the CAOB and Rodinia. However, the early Neoproterozoic tectonic framework and evolution of the Central Tianshan Block are poorly understood. Here we present a geochronological, mineral and whole-rock geochemical, and Nd–Hf isotopic study of early Neoproterozoic augen granites from the Bingdaban area. In situ zircon U–Pb dating shows that the augen granites have crystallization ages of 1002–992 Ma. The augen granites have high SiO2 (70.46–76.24 wt%) and K2O (3.36–4.96 wt%) contents, and FeOT/(FeOT + MgO) ratios (0.76–0.84), low Mg# values (29–39), slightly enriched light rare earth element patterns, enriched large-ion lithophile elements, depleted high-field-strength elements (e.g., Nb, Ta, and Ti), and high 10,000 × Ga/Al ratios (2.3–8.3) and Zr + Nb + Ce + Y contents (256–408 ppm). The augen granites have geochemical features typical of A2-type granites, variable εNd(t) values of –1.64 to +6.87 and Nd model ages of 1.78–1.08 Ga, and zircon εHf(t) values of −1.96 to +6.42 with Hf model ages of 1.98–1.48 Ga. These data indicate the granites had a juvenile crustal-dominated source, with subordinate input of ancient crustal components, and formed in a subduction-related extensional setting. The extension might have been related to slab break-off beneath Rodinia. The early Neoproterozoic igneous rocks in the Central Tianshan and adjacent microcontinental blocks in the southwestern CAOB formed at or near an active continental margin, possibly part of the circum-Rodinian subduction zone during or after the assembly of Rodinia. Similar Neoproterozoic magmatism in the Central Tianshan and Yili blocks suggests that they had a comparable crustal evolution and a close tectonic and paleogeographic relationship during the Neoproterozoic. However, significant age and isotopic differences characterize the Central Tianshan–Yili blocks and Tarim Craton, suggesting that the Central Tianshan–Yili blocks were not separated from the Tarim Craton during the early Neoproterozoic.

Source and Evolution of Subduction–Related Hot Springs Discharged in Tengchong Geothermal Field, Southwest China: Constrained by Stable H, O, and Mg Isotopes

The hydrothermal system plays a crucial role in material and energy cycling between the lithosphere and hydrosphere. In general, seafloor hydrothermal systems are one of important Mg sinks, but the situation may not be the same as it is in terrestrial hydrothermal systems. In addition, the behavior of Mg isotopes during hydrothermal circulation is still unclear. Thus, in this study, we determined the Mg isotopic compositions of the hydrothermal fluids discharged in the Tengchong region to understand better the fate of Mg in the continental hydrothermal system. The δ2H and δ18O values of the Tengchong hydrothermal fluids indicate that the recharge water sources are primary from meteoric water and influenced by the evaporation process. In contrast, the subduction–related volcanic water input is limited, except in for the Rehai area. The Mg in most of the samples is contributed by percolated meteoric water. The measured δ26Mg values range from –0.969 to 0.173‰, which are enriched in light Mg compared to the volcanic rocks of Tengchong. Combined with the precipitation dissolution of carbonate, we calculated the δ26Mg value for the endmember fluid before precipitation, which shows that the process of carbonate precipitation changes the Mg isotope of the fluid, substantially. The Shiqiang (SQ) vent is unique among all of the samples, characterized by an extremely a high δ26Mg value and Mg concentration, and it is estimated that it could have been mixed with an upper crustal material. This also reveals the diversity of the hydrothermal fluid material sources in the subduction zone.

How do pre-existing weak zones and rheological layering of the continental lithosphere influence the development and evolution of intra-continental subduction?

• The depth of a preexisting weak zone is one of the key factors in initiating intra-continental subduction. • Crustal internal friction angle plays an important role in altering the deformation type of a continental collision system. • The reactivation of the suture may be an auxiliary condition for the evolution of the south-dipping continental subduction underneath the North Pamir. Intra-continental subduction is of special importance for studying the formation of intra-continental orogens, crust-mantle structural evolution, and the far-field effects of continental collision, whose mechanism is still a matter of discussion. In this work, we investigated the role of pre-existing weak zones and the continental lithospheric rheological layering in the formation and evolution of the intra-continental subduction based on a 2D finite element numerical technique. The model results indicate that the deeper the intra-continental weak zone is and the faster the convergence velocity is, the more likely it is to develop into a new intra-continental subduction. Altering the rheological strength of the overriding plate may not have a substantial impact on the intra-continental subduction mode when the depth of the pre-existing weak zone is larger than half of the lithospheric thickness. In contrast, the lithospheric rheological strength is closely related to the continental collision system’s deformation style: Models with a weaker overriding plate are inclined to delaminate continuously under collision, whereas a strong overriding plate results in the subducting plate’s roll-back. The reactivation of the suture that runs deep into the lithosphere as a result of the Indian-Asian continental collision could be one of the crucial factors controlling the formation of the south-dipping subduction under the North Pamir.

Re Variation Triggered from the Paleo-Pacific Plate Evolution: Constrains from Mo Polymetallic Deposits in Zhejiang Province, South China Mo Province

Although highly dispersed, critical Re metal has attracted lots of attention from geoscientists, the controlling factors of Re-content variation are not completely understood, especially with regards to the genetic relationship between Re-bearing Mo polymetallic deposits and plate subduction evolution. It is well documented that the South China Mo Province, in Zhejiang Province, is characterized by multi-stage Mo polymetallic mineralization associated with Paleo-Pacific plate subduction. The Xianlin Mo(Cu)–Fe deposit occurs in Western Zhejiang as porphyry mineralization or skarn mineralization between the granodiorite and limestone. Zircon U–Pb analysis of the ore-forming granodiorite yields a Concordia age of 150.8 ± 1.1 Ma. Six molybdenite samples have relatively high Re contents (128.9~155.7 ppm) and deliver a weighted mean model age of 149.6 ± 1.3 Ma. These geochronological data suggest the Xianlin polymetallic mineralization was genetically related to the granodiorite in the Late Jurassic. Moreover, a new compilation of reliable Re contents and Re–Os isotope age data in Zhejiang Province indicates a decreasing trend in Re contents of molybdenite from the Jurassic Fe-/Cu-dominated Mo mineralization stage to the Cretaceous PbZn-enriched Mo mineralization stage in the South China Mo Province. Based on previously proposed models relating tectonic, magmatic, and hydrothermal processes, it is suggested that the Jurassic Re-enriched Mo mineralization, associated with I-type granitoids, formed in a compressive setting during the low-angle subduction of the Paleo-Pacific slab, whilst the Cretaceous Re-poorer Mo/Mo–Pb–Zn mineralization, related to both I- and A-type granitoids, formed in an extensional back-arc setting triggered by the rollback of the Paleo-Pacific slab.

How do pre-existing weak zones and rheological layering of the continental lithosphere influence the development and evolution of intra-continental subduction?

How do pre-existing weak zones and rheological layering of the continental lithosphere influence the development and evolution of intra-continental subduction?

Subduction Evolution Controlled Himalayan Orogenesis: Implications from 3-D Subduction Modeling

Himalayan orogenesis remains enigmatic in terms of Tibetan Plateau geodynamics originating from the Cenozoic India–Eurasian continental collision. India underthrusts below Tibet to the Yarlung–Tsangpo suture, which has been identified as the northernmost boundary for underplating. However, the way in which the historical evolution of continental subduction induces plateau uplift and the way it controls the variation in uplift between outboard and inboard areas is still unclear. To interpret the evolutionary mechanisms involved in the Himalayan growth history, we constructed different 3-D dynamic models at important stages to address these questions related to the formation of the Himalayas on the basis of paleoenthalpy evidence encoded in fossil leaves from recently documented assemblages in southern Tibet. The results show that (1) the effect of crustal thickening was the predominant factor in the early evolution from the Paleocene to the early Eocene, which resulted in a moderate growth rate. (2) The consecutive slab break-off eastward from the western syntaxis and the associated slab rebound significantly accelerated orogenesis from the late Eocene to the Oligocene. The upwelling asthenospheric flow was a key control of increasing crustal buoyancy, which resulted in the fastest growth of the Himalayas during the early Miocene. (3) Thereafter, the gradually enhanced monsoon and surface erosion during accompanying the increasing mountain height resulted in a slowdown of the orogenic rate, which counterbalanced the buoyant force produced by asthenospheric flow driving continuous Himalayan growth.

Provenance and tectonic setting of the Sumdo Formation in the Lhasa Terrane, Tibet: Implications for early subduction evolution of the Sumdo Paleo–Tethys Ocean

The tectonic setting of the Tangjia–Sumdo metamorphic belt, which separates the North Lhasa Block (NLB) and South Lhasa Block (SLB), remains controversial. Moreover, the sedimentary evolution of the Tangjia–Sumdo metamorphic belt has been little studied, particularly during the early subduction evolution of the Sumdo Paleo-Tethys Ocean (SPTO). Here we present petrological, geochemical, and zircon U-Pb geochronological data on the clastic rocks and metabasaltic interlayers collected from the Sumdo Formation in the Sumdo area. The Sumdo Formation consists of metamorphosed sandstone, graywacke, phyllite, and interlayered metabasalts. Detrital zircon U-Pb dating of the metasandstone and quartz schist identified a youngest zircon age group of 389–316 Ma, which have negative εHf(t) values (−20.5 to −1.5) and old Hf crustal model ages (TCDM = 1450–2654 Ma). The zircon trace element features and Hf isotopic compositions show that the sedimentary detritus was sourced from the NLB, which indicates that this terrane experienced Late Devonian–early Carboniferous magmatism. The ages of the youngest detrital zircons (YC2σ age of 324 ± 13 Ma) indicate that the Sumdo Formation formed mainly during the late Carboniferous–early Permian. Combined with geochemical data of the metabasalts and other regional geological evidence, we conclude that the Sumdo Formation formed in an initial forearc basin. The Sumdo Formation was likely deposited along an active continental margin according to the in–situ suture model, where northward subduction of the Sumdo Paleo–Tethys oceanic lithosphere occurred beneath the NLB. As such, the Sumdo Formation provides a sedimentary record of the early subduction evolution of the SPTO.

Short-lived intra-oceanic arc-trench system in the North Qaidam belt (NW China) reveals complex evolution of the Proto-Tethyan Ocean

Abstract Recognition of any intra-oceanic arc-trench system (IOAS) could provide invaluable information on the tectonic framework and geodynamic evolution of the vanished ocean basin. The Tanjianshan Complex and mafic-ultramafic rocks along the North Qaidam ultra-high pressure metamorphic belt in NW China record the subduction process of the Proto-Tethyan Ocean. Four lithotectonic units, including island arc, ophiolite, forearc basin, and accretionary complex, are recognized based on detailed field investigation. They rest on the northern margin of the Qaidam block and occur as allochthons in fault contact with underlying high-grade metamorphic rocks. The ophiolite unit mainly consists of ultramafic rocks, 527–506 Ma gabbro, 515–506 Ma plagiogranite, dolerite, and massive lava. High-Cr spinels in serpentinite, dolerite with forearc basalt affinity, and boninitic lava collectively indicate a forearc setting. The accretionary complex, exposed to the south of the ophiolite complex and island arc, is highly disrupted and contains repeated slices of basalt, 495–486 Ma tuff, chert, limestone, and mélange. Tuffs with positive zircon εHf(t) values indicate derivation from a nearby juvenile island arc. These lithotectonic units, as well as the back-arc basin, are interpreted to constitute a Cambrian IOAS that formed during the northward subduction of the Proto-Tethyan Ocean. Combined with regional geology, we propose a new geodynamic model involving short-lived Mariana-type subduction and prolonged Andean-type subduction to account for the complex evolution of the Proto-Tethyan Ocean. The reconstruction of a relatively complete IOAS from the North Qaidam belt not only reveals a systematic evolution of intra-oceanic subduction but also advances our understanding of the subduction and accretion history of the Proto-Tethyan Ocean.

Geodynamic evaluation of the pacific tectonic model for chortis block evolution using 3D numerical models of subduction

Past position and subsequent tectonic evolution of the Chortis Block represents one of the major uncertainties related to the geodynamic evolution of southern Mexico and northern Central America. This uncertainty has led to different hypotheses and evolutionary models based on tectonic, geological and paleographic studies that attempt to explain the origin and paleoposition of the Chortis Block. Among these models, a challenging hypothesis suggests the Chortis Block originates from a west-southwest position in the Pacific rather than attached to the Mexican margin. Additionally, this model proposes a high rate (200 km/20 Myr) of trench erosion along the Mexican margin just north of the present-day Chortis Block. In the Pacific model, the Chortis Block is being pushed eastward by the subduction movement of the Farallon and Cocos oceanic plates, and is finally captured by the Caribbean plate. In this study, the evolution of subduction along the Mexican and Central America margins is investigated numerically taking into account kinematic reconstruction of the Pacific model. The main goal is to investigate validity of the Pacific model for Chortis Block evolution by comparing our predicted slab geometries with present-day tomographic images. The kinematic 3D numerical models are integrated 45 Myr in time and use a mantle reference frame for tectonic plate evolution of the study area. These models generate subducting plate geometries that are only partial consistent with what is observed from seismic tomography, they reproduce realistic slab geometries only for the region located beneath the North America Plate, where the Cocos slab is flat at shallow depths and also stagnates along the mantle transition zone. On the other hand, for the region located beneath the Chortis Block the slab penetrates vertically through the upper mantle and even curls back in the upper portion of the lower mantle, a geometry that is inconsistent with seismic tomography images which shows the presence of east-dipping slab that tends to stagnate at 1000 km depth. This study indicate that the Pacific model is not an appropriate hypothesis for explaining the Chortis Block origin, and future studies should focus on alternative tectonic models.

Reconstructing Jurassic‐Cretaceous Intra‐Oceanic Subduction Evolution in the Northwestern Panthalassa Ocean Using Ocean Plate Stratigraphy From Hokkaido, Japan

AbstractPlate reconstructions of the Panthalassa Ocean typically portray a simple system of diverging plates surrounded by active margins, yet geological and seismic tomographic records demonstrate that intra‐oceanic subduction existed. Here, we reconstruct the plate tectonic evolution of the pre‐Cretaceous intra‐oceanic Oku‐Niikappu island arc, remnants of which are exposed on Hokkaido, Japan. This arc formed at a Jurassic subduction zone separating two oceanic plates: the Izanagi Plate and the here proposed ‘Izanami’ Plate. The Oku‐Niikappu arc was previously shown to have gone extinct in an intra‐oceanic setting, was subsequently (hyper)extended, and overlain by Berriasian cherts. The extinct arc remained on the Panthalassa ocean floor for ∼45 Myr until its ∼100 Ma accretion to Hokkaido, revealing an original position far from the continental margin and likely above the previously identified Telkhinia slabs. We show that arc extinction coincided with a northwestern Panthalassa plate reorganization recorded by a ∼30° change in spreading direction, and that extinction and subsequent extension of the arc is straightforwardly explained by subduction of the Izanami‐Pacific ridge followed by continued divergence between the Izanagi and Pacific plates. From our reconstruction, it follows that the outer zone of Japan, to which the accretionary complex in which the Oku‐Niikappu Complex resides belongs, was separated from the inner zone by a back‐arc basin during the Early to mid‐Cretaceous. This study illustrates the value of accretionary orogens in the development of plate reconstructions of lost oceanic plates and ancient continental margins, particularly when combined with seismic tomographic and marine geophysical data sets.

On exhumation velocities of high-pressure units based on insights from chemical zoning in garnet (Tianshan, NW China)

Exhumed high-pressure rocks formed during subduction provide insight into global geodynamic and plate tectonic processes. Their metamorphic history is partially encoded in their mineralogy and mineral chemistry, which can record the P-T evolution and rates of subduction and exhumation. Here, we study a garnet-glaucophane-albite schist from the SW Tianshan in China in which complexly zoned garnet displays a distinct and sharp oscillatory Mn-zonation and irregular, jagged Ca distributions with peak widths of ∼5 μm. Combined Raman quartz-in-garnet barometry and carbonaceous matter thermometry suggests a P-T path that reached ∼2.5 GPa and was followed by ∼0.9 GPa of nearly-isothermal decompression at temperatures of 530–540°C, associated with final garnet growth. This P-T path is largely confirmed by thermodynamic modeling. Three different approaches to modeling coupled Ca, Mn, Mg and Fe diffusion in garnet were used to estimate the duration between establishment of high-resolution Mn, Fe and Ca zoning during growth and cooling below nominal closure conditions upon exhumation. All models agree best with conditions of ∼530°C maintained for around 100 kyr, though slight differences are found depending on published diffusion coefficients. For the constrained exhumation P-T path, we infer that the studied sample experienced exhumation velocities in the range of 40–45 cm/yr from ∼80 km to ∼50 km depth. We propose an exhumation model involving diapiric Stokes flow along the plate interface and low viscosities around the exhumed unit that may facilitate such high exhumation velocities.