Slab Subduction Research Articles

Two Mechanisms Generating Intraplate Volcanism in Southeast Asia: Insights From Mantle Transition Zone Discontinuities

AbstractDetermining the mechanisms responsible for intraplate volcanism ‐ such as slab devolatilization melting versus active mantle plumes ‐ remains a challenge. The greater South China Sea (SCS) region has experienced extensive intraplate Cenozoic volcanism across areas including Hainan, Southeast Indochina, northern Borneo, the northern SCS, and the post‐spreading SCS basin. The prevalence of volcanism distributed widely across this region prompts fundamental questions about the key geodynamic processes driving such diverse magmatic activities. In this study, we elucidate the mantle transition zone (MTZ) discontinuities in this region using SS precursors, which helps to overcome the sparse seismic coverage due to its predominantly oceanic setting. We collected over 16,000 high‐quality seismograms that sample the upper mantle and MTZ beneath this region from global earthquakes and stations. After correcting for the effects of shallow crustal variations and upper mantle heterogeneity on traveltimes of SS phases and their precursors, we unveil lateral variations in the MTZ boundaries (d410 and d660) and intricate features of the mid‐MTZ reflectors (S520S). Significant MTZ thinning and normal S520S waveforms beneath Hainan provide compelling evidence for mantle upwelling through the MTZ. Conversely, the evident splitting of S520S beneath the northern SCS, Southeast Indochina, and northern Borneo, all characterized by stagnant subducted slabs, indicates that the volcanism in these regions likely originated from a mechanism distinct from the active upwelling beneath Hainan. Dehydration melting attributed to devolatilizing stagnant slabs in the MTZ is a potential cause for Cenozoic volcanism in these regions.

Multiple Slabs and Complex Mantle Flows in the Molucca Sea Subduction Zone

AbstractThe Molucca Sea subduction zone is unique for its active divergent double‐subducted slab located in the north of the Java and Banda subduction zones. The spatial proximity of these subduction zones would cause a complex mantle flow field. To clarify the mantle dynamics, here we present P‐wave tilting‐axis anisotropic tomography obtained by inverting a large number of local and teleseismic travel‐time data recorded at 254 seismic stations in eastern Southeast Asia. Our results suggest that the mantle structure and dynamics of the western Molucca Sea subduction zone may be remotely controlled by the Java and Banda subduction zones. Mantle convection in the big mantle wedge west of the Molucca Sea subduction zone is possibly influenced by east‐west mantle flow associated with compression of the Indo‐Australian slab as well as north‐south mantle flow related to rollback of the Indo‐Australian slab. In contrast, the eastern Molucca Sea subduction zone is virtually unaffected by other subduction zones, probably due to the ongoing slab subduction there.

Slab break-off of the Kalamaili oceanic slab revealed by the latest Carboniferous mafic–ultramafic rocks in eastern North Tianshan (NW China)

Identifying the tectonic transition from oceanic subduction to collision is crucial for tracking the final stage evolution of ancient orogenic belts. In this study, we present new geochronological and geochemical data for the Mozbaysay mafic–ultramafic complex in the Balikun area, eastern North Tianshan of the southern Central Asian Orogenic Belt. This complex had intruded into the late Carboniferous volcano-sedimentary rocks and is comprised mainly of hornblende gabbro and lherzolite. Zircon U-Pb ages of the hornblende-gabbros reveal that this complex was emplaced at ca. 305 Ma. Geochemical analyses suggest these mafic–ultramafic rocks are characterized by slight enrichment of light rare earth elements (LREEs) and relatively depleted heavy rare earth elements (HREEs), resembling enriched mid-ocean ridge basalt (E-MORB). They also exhibit restricted (87Sr/86Sr)i ratios (0.702396–0.704295) and εNd(t) values (+7.0 to +9.1), indicative of a depleted mantle source with minimal crustal contamination. Incompatible element ratios (i.e., Nb/Ta, Zr/Hf, Rb/Nb, and Ba/Nb) suggest the involvement of subducted slab-derived aqueous fluids in their mantle source. These collectively indicate that the parental magmas of the Mozbaysay mafic–ultramafic rocks may have been generated by a mixed mantle source consisting of the E-MORB-like asthenospheric mantle, subcontinental lithospheric mantle (SCLM), and hydrous fluids from subducted slab. Furthermore, a slab break-off model is proposed to explain the generation of these latest Carboniferous mafic–ultramafic rocks. Integrating these findings with regional geological data, we propose that the tectonic transition from subduction (slab roll-back) to collision (slab break-off) along the Kalamaili suture zone occurred at ca. 305–300 Ma.

Mantle deformation in the highly oblique indo-burma subduction system inferred from shear wave splitting measurements

We utilized shear wave splitting analysis of teleseismic SKS, SKKS, and PKS phases to infer upper mantle deformational fabrics across a substantial area of Southeast Asia, where splitting measurements were previously limited. We used newly available permanent and temporary broadband seismic networks deployed across the Indo-Burma subduction zone and the eastern Indochina peninsula. The resulting 492 well-constrained splitting and 654 null measurements from 185 stations reveal clear large-scale patterns in the mantle deformational fabrics in response to the highly oblique active subduction and a large transform plate boundary. We identified two distinct domains of mantle deformation fabrics in the western Burma microplate and the eastern Indochina peninsula. In the former, trench parallel N-S fast polarization directions with an average lag time (δt) of 1.9 s are observed beneath the Indo-Burman Ranges. We suggest the observed splitting is partly due to anisotropy in the sub-slab region and relates to shear induced by the north moving Indian plate. The lithospheric fabric within the Indo-Burman Ranges and underlying subducting slab fabric contribute to produce the observed average δt of 1.9 s. The δt value decreases to an average of 1.0 s towards the back-arc until we reach the dextral Sagaing fault. In the second domain, starting approximately 100 km east of the Sagaing fault, we observe a consistent E-W fast direction with an average δt of 1.10 s in the eastern Shan-Thai and Indochina blocks. We interpret the E-W fabric as due to the deformation associated with the westward spreading of the Hainan mantle plume, possibly driven by overriding plate motion. Low velocities in the shallow mantle and late Cenozoic intraplate volcanism in this region support the plume-driven asthenospheric flow model in the Indochina peninsula. The sudden transition of the fast polarization direction from N-S to E-W along the eastern edge of the Burma microplate indicates the Sagaing fault acts as a mantle flow boundary between the subduction dominated trench parallel flow to the west and plume induced asthenospheric flow to the east. We also observed no net splitting beneath the Bengal basin which is most likely due to the presence of frozen vertical fabric resulting from the Kerguelen plume activity during Early Cretaceous.

The petrogenesis of Cenozoic basalts from Daihai, western North China Craton: Constraints from 40Ar-39Ar chronology, major and trace elements, and Sr-Nd-Pb-Hf isotopes

The Daihai Cenozoic intraplate basalts are distributed in the western North China Craton (NCC), which is a part of the Cenozoic volcanic province in eastern China, and they are all alkaline basalts. The fine-scale determination of their mantle source region properties, partial melting mechanisms, and petrogenesis can provide crucial information for exploring lithospheric destruction and thinning in the western NCC. 40Ar-39Ar dating of potassium feldspar grains from the Daihai alkaline basalts yielded plateau ages of 18.22 ± 1.84 Ma and 26.86 ± 0.72 Ma, indicating that the Daihai basalts underwent multiple eruptive cycles during the Late Oligocene-Middle Miocene. These basalts exhibit ocean island basalt (OIB)-like geochemical features and were subjected to negligible crustal contamination. Moreover, basaltic magmas underwent intense fractional crystallization of olivine and clinopyroxene. The geochemical differences in the Daihai basalts were controlled by partial melting. The Daihai basaltic magmas were composed of at least two types (type I-enriched mantle and prevalent mantle) of mantle end-members that partially melted and then mixed, and the lithology of the mantle source region was predominantly peridotite. Under enriched mantle conditions, the mixing of garnet lherzolite partial melts (<1%) and spinel lherzolite partial melts (2-5%) can reasonably explain the elemental variations characteristic of the Daihai basalts. Most importantly, the melting depth and lithospheric thickness of the Daihai basalts were < 70 km or even close to 50 km, implying that the western NCC underwent lithospheric destruction and thinning during the Cenozoic. However, these effects were spatially heterogeneous. The most plausible genetic mechanism for the generation of the Daihai basalts was the coupled effects of subduction of the Pacific slab and subduction–collision of the Indo–Eurasian Plate since the Oligocene. The mantle flows generated by these two events convected, blocked and triggered upwelling mantle flows at the eastern margin of the Ordos Block. The upwelling mantle flows resulted in frequent magmatism in the region.

Stability of H2O‐Rich Fluid in the Deep Mantle Indicated by the MgO‐SiO2‐H2O Phase Relations at 23 GPa and 2,000 K

AbstractThe Earth's mantle contains significant amounts of water in the form of hydroxyl in hydrous minerals, nominally anhydrous minerals, and hydrous silicate melts. H2O fluid is thought to be present only in the shallow regions because it will always dissolve tens of weight percent of silicates by forming hydrous silicate melt in the deep mantle. Here I investigated the phase relations in the MgO‐SiO2‐H2O system by high‐pressure experiments at a pressure of 23 GPa and a temperature of 2,000 K, corresponding to the conditions at the bottom of the mantle transition zone and the topmost lower mantle. The experimental results indicate that hydrous melt can contain more than 90 wt.% of H2O, that is, it becomes H2O‐rich fluid when coexists only with stishovite. In contrast, silicate‐rich hydrous melt is formed when the system is enriched with MgO component. Therefore, H2O‐rich fluid may be stabilized in locally SiO2‐enriched rocks even at the topmost lower mantle, acting as a water source for the deep lower mantle by slab subduction. The H2O fluid also provide a possible cause for the occurrence of natural ice‐VII originated from 660 km depth.

Constraints on the subduction and closure of the Hegenshan Ocean: Magmatic and sedimentary records from central Inner Mongolia, China

Petrological, geochronological, and geochemical data from the volcano-sedimentary sequences, granitoids, and ophiolite relics of central Inner Mongolia, China, were used to reconstruct the subduction and final closure of the Hegenshan Ocean. Geochronological dating and compilation reveal four phases (ca. 360−355 Ma, 348−320 Ma, 320−310 Ma, and 310−275 Ma) of magmatism in the Uliastai continental margin. The ca. 356 Ma I-type Halatumiao granodiorite and Amanwusu ophiolite relics are subduction-related, and the Halatumiao granodiorite provides solid evidence of the northward subduction of the Hegenshan Ocean beneath the Uliastai continental margin at ca. 360−355 Ma. The ca. 348−320 Ma and 320−310 Ma volcanic rocks and granitoids constitute two linear magmatic belts roughly parallel to the Erenhot-Hegenshan ophiolite belt, which record two phases of continental arc magmatism in the Uliastai continental margin. Overall, the ca. 360−310 Ma arc magmatism shows landward migration and then oceanward migration in the Uliastai continental margin, which indicates advancing subduction and subsequent slab steepening of the Hegenshan Ocean. By contrast, the ca. 310−275 Ma magmatic rocks are dominated by I- and A-type felsic volcanic rocks, granites, and dikes, which are post-accretionary, extension-related, and pervasive in the Uliastai continental margin and Erenhot-Hegenshan ophiolite belt. A provenance shift was identified between the Benbatu and Amushan formations of the Amanwusu area of the Erenhot-Hegenshan ophiolite belt. The early detritus was derived from the early Paleozoic rocks in the Sonid Zuoqi arc belt, whereas the late detritus originated from the Early Carboniferous ophiolite relics in the Erenhot-Hegenshan ophiolite belt. The provenance shift and emplacement of pervasive extension-related magmatic rocks imply a Late Carboniferous closure of the Hegenshan Ocean. The Late Carboniferous oceanic closure event in the north of the southeast Central Asian Orogenic Belt is also evidenced by the transition of Hf isotopic composition of zircons dated between ca. 360−310 Ma and 310−275 Ma.

Identification of carbonated eclogitic source for Cenozoic intraplate alkaline basalts from western Central Asia and its geodynamic implication

Carbonated eclogite has been identified as a significant component in the formation of intraplate alkaline basalts, typically attributed to stagnated subducted oceanic slabs, such as the Pacific plate beneath eastern China. However, it is unclear whether this component exists in other tectonic settings. Here, we present new zircon UPb geochronology, mineral chemistry, and geochemical analyses of Cenozoic Tuoyun basanites and alkaline basalts from western Central Asia, where stagnated subducted oceanic slabs are absent, to elucidate the source properties and associated geodynamic processes. Olivine compositions (e.g., moderate Ni concentrations, low Mn/Fe and Ni/(Mg/Fe) ratios) and whole-rock chemistry (e.g., high Fe/Mn, Ca/Al, and Zr/Hf ratios, TiO2 contents, and weak anomalies of Zr-Hf-Ti) suggest a carbonated eclogite in the mantle source of the Tuoyun basanites. Trace elements (positive NbTa anomalies and absence of positive Eu anomaly) and depleted to slightly enriched Sr-Nd-Hf isotopic compositions indicate that the basanites originated from mantle sources composed of variable proportions of altered oceanic crust (AOC) and carbonates. Synthesizing regional geochemical, geophysical studies and tectonic history, we propose that the continuous India-Eurasia collision triggered mantle upwelling, leading to the formation of a widespread carbonated eclogitic mantle source beneath western Central Asia. Our findings suggest that carbonate-rich mantle sources may also be present in areas without stagnated subduction slabs.

Intermittent terrane arrival induces pulses of inland tectonic cycles

Pulsing volcanotectonic cycles characterized by short-lived (10 s Myr) switches in magmatic and tectonic patterns have been broadly identified in active margins. However, the specific mechanism that causes these switches remains ambiguous, i.e., whether the subduction continuity and/or terrane arrival (accretion, underthrusting, or subduction of buoyant continental/oceanic blocks with a thicker crust than their surrounding oceanic plate) plays a crucial role in controlling the observed volcanotectonic cycles remains controversial. Here, by modeling subduction processes involving the sequential arrival of buoyant terranes, we show that 1) in scenarios where the oceanic plate is weakly coupled with terranes, the entraining of the terrane into the subduction induces slab breakup to occur between the partially subducted terrane and its adjacent oceanic slab, with the terrane then rebounding and moving away from the trench. This evolution causes switches in the magmatic and tectonic patterns within the overriding plate. 2) In these models, the exposed terrane materials can preserve characteristic pressure-temperature-time trajectories, i.e., nearly isothermal compression to isothermal decompression and/or isobaric heating after partial exhumation. 3) slab breakup does not guarantee the occurrence of a trench jump; only when the terrane has a medium scale (∼300 km) will a new trench tend to develop prominently behind the rebounded terrane. 4) in scenarios where the composite slab (terrane and oceanic portions) resists yielding deformation and the terrane density is close to that of the oceanic plate (<0.6% density contrast), continuous subduction will occur. In this latter scenario, slab deformation (revealed by subduction angle and slab curvature), instead of slab breakup, will control the magmatic and tectonic patterns in the overriding plate. By further comparing model results with observations, we demonstrate that intermittent subduction interrupted by the subduction of terranes can be a tectonic driver for episodes of compression-to-extension transformations and magmatism (or piston-like volcanotectonic cycles) in two representative accretionary belts — with or without trench jumps (exemplars in south-central-Tibet and eastern-Mediterranean). In contrast, continuous subduction with strongly coupled upper and subducting plates could have contributed to similar cycles in the example without accretion (e.g., the Altiplano in the Central Andes). Therefore, over 10 s of Myr, the scale and frequency of terrane arrivals could essentially control the specific motion pattern of the subducting plate, creating the observed short-lived volcanotectonic switches at these three subtypes of active margins.

Fossil metasomatized and newly-accreted fertile lithospheric mantle volumes beneath the Bakony-Balaton Highland Volcanic Field (central Carpathian-Pannonian region)

This study focuses on upper mantle xenoliths from Sabar-hill (Bakony-Balaton Highland Volcanic Field), a newly discovered locality of mantle xenoliths, situated in the central Carpathian-Pannonian region (CPR). The investigated seven peridotite and one pyroxenite xenoliths record two episodes in the geochemical evolution of the local subcontinental lithospheric mantle. The moderately deformed Group I xenoliths reveal modal contents (≥9 vol% orthopyroxene), major element characteristics (depletion in Al2O3 and CaO in orthopyroxenes) and trace element features (e.g., enrichment in U, Pb and Sr in clinopyroxenes) typical in supra-subduction zones (i.e., mantle wedge) xenoliths. All these suggest that the studied xenoliths were formed by silica-rich fluid/melt - peridotitic wall rock interactions. The tectonic background to the formation of Group I xenoliths is likely linked to the subduction of an oceanic slab during the Mesozoic-Paleogene. This may have happened far from the current position of Sabar-hill, to where the lithosphere, including the metasomatized mantle volume, was transferred via plate extrusion. The less deformed, dominantly protogranular and porphyroclastic Group II xenoliths are depleted in light rare earth elements (LREE). In contrast, they are enriched in clinopyroxene (≥11 vol%), incompatible major elements (e.g., Al and Na) and heavy rare earth elements (HREE), by a minimum of four times compared to primitive mantle. All these suggest that Group II xenoliths represent a geochemically fertile mantle partition and the degree of partial melting affected this mantle section was limited (≤7% based on trace elements). Group II xenoliths have higher equilibrium temperatures (1088–1168 °C) compared to Group I xenoliths (1018–1089 °C). The lack of petrographic and geochemical evidence pointing to metasomatism-related annealing in Group II xenoliths suggests that they are derived from a greater depth. It follows that Group II xenoliths could represent prior asthenosphere volumes which recently accreted to the lithosphere. The lithospherization of the upper part of the asthenosphere may have initiated after the tectonic inversion of the region, when compressional tectonic regime beneath the CPR gradually overtook the former extensional one in the last ∼8 Ma. The high bulk structural hydroxyl content in both Group I (22–1699 ppm) and II (55–320 ppm) xenoliths is likely linked to the subduction events that took place in the CPR from the Mesozoic onwards. These events transported significant amounts of fluid (H2O) back to the upper mantle, affecting both the mantle wedge and the asthenosphere.

Reassessing early to mid-Neoproterozoic (∼885–720 Ma) tectonic evolution of the southwestern Jiangnan Orogen: A detrital zircon perspective

The Neoproterozoic Jiangnan Orogen (JNO) marks the consolidation of the South China Block, sparking ongoing debates on its early to mid-Neoproterozoic tectonic evolution. In this study, we have gathered detrital zircon data from the southwestern JNO, presenting new detrital zircon U-Pb ages alongside trace element data from the mid-Neoproterozoic Banxi Group and Sinian System within this region. Our objective is to illuminate the early to mid-Neoproterozoic tectonic evolution within the southwestern JNO through a comprehensive detrital zircon analysis. The obtained dating results derived from detrital zircons reveal maximum depositional ages of ∼733 Ma and ∼713 Ma for Mobin Formation (Banxi Group) and Jiangkou Formation (Sinian System) sedimentary rocks, suggesting a proximal source, possibly originating from magmatic rocks within the southwestern JNO. Within the paleo-geographic context, it is plausible that these samples were deposited within the deep-water expanse of the Nanhua Rift Basin. Trace element analysis of detrital zircons within the southwestern JNO highlights distinctive fluctuations, delineating three crucial time points (∼885 Ma, ∼835 Ma, and ∼810 Ma): the beginning of early Neoproterozoic subduction of the South China Ocean towards the Yangtze Block beneath the southwestern JNO, the subduction slab rollback, and the Nanhua rifting initiation. The discernible shifts in detrital zircon age distributions within these sedimentary rocks serve as robust validation for the proposed model and crustal reworking. We then propose that the southwestern JNO was in a lull of magmatism before ∼885 Ma, then transitioned to a convergent setting around 885–835 Ma, subsequently subduction slab rollback after ∼835 Ma. The time of final assembly of the Yangtze and Cathaysia blocks in the southwestern JNO is no later than ∼810 Ma. This transformative period (∼810–720 Ma) witnessed the development of a rift basin concurrently with the accumulation of deep-water deposits within the southwestern JNO.

Implications of Late Triassic–Middle Jurassic detrital zircons from the southeastern margin of the South China Block for the Paleo-Pacific Plate subduction

There existed an early Mesozoic tectonic regime transformation in the South China Block (SCB). Although it was believed that this transformation was closely related to the subduction of Paleo-Pacific Plate, the process and initial time of the subduction of the Paleo-Pacific Plate have been controversial for a long time. Based on the published Upper Triassic and Middle Jurassic succession detrital zircons geochronology from the southeastern margin of the SCB, the newly obtained Late Triassic to Middle Jurassic detrital zircons Hf isotope data, and the available magmatic age data and Hf isotope data, this study discussed the subduction of the Paleo-Pacific Plate. The Hf isotope composition of zircon grains selected from the Early Mesozoic strata of the SCB, spanning from the Late Triassic to the Middle Jurassic, exhibits a consistent temporal trend with that of the Triassic–Middle Jurassic magmatic rocks. Both display a transition from negative to positive εHf(t) values, indicative of a gradually increasing contribution of mantle-derived materials from the Triassic to the Middle Jurassic. Zircon trace elements indicate that a magmatic arc appeared outside the southeastern margin of the SCB at 200–190 Ma and continued to develop into the Middle Jurassic, which may have been generated by the Paleo-Pacific Plate subduction. This study proposed that the Paleo-Pacific Plate subduction was confined to the southeastern coast of the SCB in the Late Triassic, and the subduction of the flat slab was halted by obstruction at the end of the Late Triassic. In the Early Jurassic, the Paleo-Pacific Plate began to subduction again, and arc-related magmatic rocks were formed along the coast of SCB. At the same time, due to the remote effect of the subduction of the Paleo-Pacific Plate, the weak tectonic belt existing in the Nanling area was reactivated, resulting in the Nanling area in the intraplate extensional setting. Subsequently, the continuous subduction of the Paleo-Pacific Plate led to the thickening of the crust along the SCB in the Middle Jurassic, and the Nanling area was still under in the extensional setting.

Anisotropic thermal conductivity of antigorite along slab subduction impacts seismicity of intermediate-depth earthquakes

Double seismic zones (DSZs) are a feature of some subducting slabs, where intermediate-depth earthquakes (~70–300 km) align along two separate planes. The upper seismic plane is generally attributed to dehydration embrittlement, whereas mechanisms forming the lower seismic plane are still debated. Thermal conductivity of slab minerals is expected to control the temperature evolution of subducting slabs, and therefore their seismicity. However, effects of the potential anisotropic thermal conductivity of layered serpentine minerals with crystal preferred orientation on slab’s thermal evolution remain poorly understood. Here we measure the lattice thermal conductivity of antigorite, a hydrous serpentine mineral, along its crystallographic b- and c-axis at relevant high pressure-temperature conditions of subduction. We find that antigorite’s thermal conductivity along the c-axis is ~3–4 folds smaller than the b-axis. Our numerical models further reveal that when the low-thermal-conductivity c-axis is aligned normal to the slab dip, antigorite’s strongly anisotropic thermal conductivity enables heating at the top portion of the slab, facilitating dehydration embrittlement that causes the seismicity in the upper plane of DSZs. Potentially, the antigorite’s thermal insulating effect also hinders the dissipation of frictional heat inside shear zones, promoting thermal runaway along serpentinized faults that could trigger intermediate-depth earthquakes.

Inversion of marine magnetotelluric and controlled source electromagnetic data with site gap and low signal-to-noise ratio

SUMMARY Marine magnetotelluric (MT) and controlled source electromagnetic (CSEM) methods have been routinely applied to survey the crustal and upper mantle electrical resistivity structures beneath the sea floor. In practice, there are inevitably site gaps and contamination by noises in the collected data because of lost ocean bottom electromagnetic receivers, unusable data, and difficulties in deploying instruments near deep trenches. So far, it remains unclear to what degree those factors will lower the resolution and the credibility of marine MT and CSEM inversion models. In this paper, we investigate the individual and combined effects of site gaps and data noises on the inversion models through synthetic analyses based on a simple block resistivity model and a realistic resistivity structure derived from the Mariana Trench. The results suggest that data with a sufficiently high signal-to-noise ratio can reasonably recover the subseafloor structures in the area of data gap. The transverse electric mode and tipper data from the MT method are much more sensitive to the structure near the site gap. The joint inversion of MT and CSEM data would improve the model's resolution at the site gap area. The inversion of data with a relatively low signal-to-noise ratio, for example, 10 per cent, can recover the structures with few artefacts if there is no site gap. But if the site gap and noisy data are combined, even a joint inversion cannot correctly recover the burial depths and geometries of the anomalous bodies beneath the site gap where vertical strips are likely present. To improve the model's resolution and suppress inversion artefacts, we propose constraining part of the model with as much a prior information as possible. Specifically, for a survey in the subduction zone, we could reduce the penalties on the model's smoothness at the upper and low interfaces of the resistive subduction slab, or even fix the resistivity of the resistive slab with the help of other information, if any. The inversion models shown in this paper provide valuable references for the site design before marine MT and CSEM surveys as well as for interpreting real data inversion models that may be subject to the same biases introduced by the site gap and noise.

Tectonic evolution of the South Pamir Orogen: Insights from the Permian to cretaceous magmatism

The Pamir orogen originated through the amalgamation of various arc terranes. Despite undergoing numerous studies, the subduction polarity and time regarding the accretion processes remain inadequately constrained. To address these issues, we carried out a detailed study of zircon and apatite U-Pb ages, zircon Hf isotopes, and whole-rock geochemistry for the Permian to Cretaceous gabbro, volcanic rocks and granitoids from the Southern Pamir and the Rushan-Pshart suture. Three stages of magmatism are identified, including 279–266 Ma gabbro and andesite to dacite, 209–208 Ma granitoids, and 108–103 Ma gabbro and granitoids, respectively. The Permian gabbro and basaltic andesite to dacite in Pshart Range are enriched in light rare earth elements (LREE) and large ion lithophile elements (LILE), but are depleted in Ta, Nb, and Ti, suggesting formation from an enriched mantle source in a subduction zone. The Late Triassic I-type granitoids are typical arc calc-alkaline magmatic rocks with εHf(t) values varying from −8.7 to −3.5, which mainly sourced from the partial melting of mixed crustal sources predominately consisting mafic rocks with minor metasedimentary rocks. The Early Cretaceous Pshart granites are I-type granites with relatively low εHf(t) values ranging from −14.6 to −2.9, suggesting that they were crust-derived rocks. The Early Cretaceous Bazardara granites are identified as A2-type granites with the εHf(t) values ranging from −18.3 to +5.3, and were derived from mixing of mantle-derived and crust-derived magmas, or interaction between mantle melts and old crustal rocks. The early Cretaceous gabbro exhibits -MORB-like geochemical characteristics with flat REE pattern (La/YbN = 4.96), moderate enrichment of LILEs and weak depletion of Nb, Ta and Ti (Nb/LaPM = 0.94). Thus, the early Cretaceous A2-type granite, E-MORB-like gabbro and volcanic rocks were formed in an extensional setting. Our new data indicate that the southward subduction of the Rushan-Pshart oceanic slab started in the Early Permian and continued to the Late Triassic. Moreover, combining the previous data, we conclude that the Early Cretaceous granites and gabbro in the Southern Pamir formed in an extensional environment due to the southward rollback after northward near-flat subduction of the Neo-Tethyan oceanic lithosphere.

Neoarchean granitoid magmatism and geodynamic process in the northeastern North China craton

The composition of Archean granitoid rocks changed from predominantly tonalite-trondhjemite-granodiorite (TTG) gneisses in the early Archean (4−3 Ga) to diversified granitoid rock assemblages in the late Archean (3.0−2.5 Ga), marking a crucial transformation in the geodynamic processes of early Earth. However, the reason for this major transition remains enigmatic because the petrogenetic features of different granitoid assemblages and their crust-mantle interactions during different periods are poorly understood. We use variations in the spatial-temporal distribution, lithological association, chemical composition, and petrogenesis of Neoarchean (2.7−2.5 Ga) granitoids and inferred correlative crust-mantle interactions in the Eastern Liaoning Range (ELR) of the northeastern North China craton to explore this geodynamic transition. The early Neoarchean (ca. 2.7 Ga) ELR granitoids were dominated by TTG gneisses, and the late Neoarchean (2.6−2.5 Ga) ELR granitoid typology and compositions became more complex, changing into TTGs and more K2O-rich granitoid rocks. The TTGs can be subdivided into a high-Ca group and a low-Ca group: The 2.71−2.68 Ga high-Ca group TTG magma originated from partial melting of subducted juvenile oceanic crust, and the low-Ca group TTG magma was derived from fractionation crystallization of the high-Ca group TTG magma. The chemical composition of the magmatic sources played a dominant role on the 2.60−2.50 Ga TTG magmatism: the high-Ca and low-Ca group TTG magmas came from low-K mafic rocks and tonalites, respectively. The 2.58−2.49 Ga K2O-rich granitoids can be divided into three petrogenetic series: (1) The high-Ca-Mg group K2O-rich granitoid magma originated from partial melting of high-K mafic rocks, (2) the low-Ca-Mg group K2O-rich granitoid magma was derived from partial melting of sedimentary rocks, and (3) the transition group K2O-rich granitoid magma was sourced from metagreywackes. The 2.71−2.68 Ga TTGs were generated in an island arc belt, and subducted slab melting and subsequent magmatic differentiation were the dominant mechanisms of the TTG magmatism. The 2.60−2.50 Ga diversified granitoids were formed in the oceanic-continental subduction process under the active continental margin; the complicated oceanic slab subduction and arc-arc and arc-continent collisions contributed to the diversity of late Neoarchean granitoid magmatism.

K- and Na-rich davemaoite inclusion in diamond is not inherited from deeply subducted oceanic crusts

Deep potassium (K) and sodium (Na) cycles can affect the physiochemical properties and geodynamic processes of the Earth's interior, to some extent. Recently reported natural K- and Na-rich davemaoite inclusion in superdeep diamond (Tschauner et al., 2021), which is still controversial (Walter et al. 2022a; Tschauner et al., 2022), motivated us to reinvestigate the recycling and partitioning behaviours of K and Na in the subducting oceanic crust under deep-mantle P-T conditions. In this study, we experimentally determined the partitioning of Na and K, as well as other major elements, in a natural phengite-rich eclogite system under 5–30 GPa and 850–1550 °C, relevant to P-T conditions from the deep upper mantle to the uppermost lower mantle. We found that the produced davemaoite with the perovskite structure contained minimal amounts of other elements, including K and Na. Therefore, we conclude that K- and Na-rich davemaoite inclusion in the superdeep diamond cannot be inherited from deeply subducted oceanic crusts. In addition, our experimental results further support previous arguments overestimating the K and Na, and other elemental contents in natural davemaoite, which may further indicate that the reported diamond “inclusion” with abnormally high heteroatoms (Tschauner et al., 2021) is not a davemaoite from the lower mantle (Walter et al. 2022a). Finally, we combined the high-pressure phase relationship in the upper continental crust, oceanic basalt, and pyrolite systems from previous studies, and constructed a geodynamical model for the deep K and Na cycles along with the subduction of the oceanic slabs into the uppermost lower mantle.

The influence of viscous slab rheology on numerical models of subduction

Abstract. Numerical models of subduction commonly use diffusion and dislocation creep laws from laboratory deformation experiments to determine the rheology of the lithosphere. The specific implementation of these laws varies from study to study, and the impacts of this variation on model behavior have not been thoroughly explored. We run simplified 2D numerical models of free subduction in SULEC, with viscoplastic slabs following (1) a diffusion creep law, (2) a dislocation creep law, and (3) both simultaneously, as well as several variations of model 3 with reduced resistance to bending. We compare the results of these models to a model with a constant-viscosity slab to determine the impact of the implementation of different lithospheric flow laws on subduction dynamics. In creep-governed models, higher subduction velocity causes a longer effective slab length, increasing slab pull and asthenospheric drag, which, in turn, affect subduction velocity. Numerical and analogue models implementing constant-viscosity slabs lack this feedback but still capture morphological patterns observed in more complex models. Dislocation creep is the primary deformation mechanism throughout the subducting lithosphere in our models. However, both diffusion creep and dislocation creep predict very high viscosities in the cold core of the slab. At the trench, the effective viscosity is lowered by plastic failure, rendering effective slab thickness the primary control on bending resistance and subduction velocity. However, at depth, plastic failure is not active, and the viscosity cap is reached in significant portions of the slab. The resulting high slab stiffness causes the subducting plate to curl under itself at the mantle transition zone, affecting patterns in subduction velocity, slab dip, and trench migration over time. Peierls creep and localized grain size reduction likely limit the stress and viscosity in the cores of real slabs. Numerical models implementing only power-law creep and neglecting Peierls creep are likely to overestimate the stiffness of subducting lithosphere, which may impact model results in a variety of respects.

Altered spinels act as a mirror of multi-episodic fluid metasomatism in the forearc mantle: An example from the Minhe ophiolite in Qilian Orogen, NW China

Forearc mantle peridotites commonly undergo complex melt/fluid metasomatism, which inevitably has a considerable influence on their geochemical compositions at the whole-rock scale. Serpentinites predominantly comprise the Minhe ophiolite, which is a fragment of the Proto-Tethys Ocean that outcrops in the Qilian Orogen in NW China. These serpentinites contain zoned spinels that vary compositionally from the centers to the margins, with a spinel core overgrown by ferritchromite layer and magnetite at the outermost rim. The mineral chemistry of the spinel cores (i.e., Cr#, Mg#, Cr2O3, and Al2O3) suggests that the protolith of the serpentinites was an oceanic mantle remnant in a forearc setting. Moreover, the ferritchromite exhibited increases in Fe, Mn, Sc, and Ni contents and decreases in Mg, Cr, Al, V, Zn, and Ga contents. Considering the linear correlation among elements (e.g., Co, Ni, and Zn) and the absence of mineral inclusions, the ferritchromite likely formed through ionic substitutions induced by slab fluids. The magnetite was enriched in Cr, Ni, V, Zn, Co, and Cu, but did not have a magmatic origin owing to low V concentrations. All evidence indicates that the peridotites underwent at least two episodes of fluidic metasomatism in the subduction zone. During early slab subduction, aqueous fluids released from the dehydration of the subducted slab directly serpentinized the forearc mantle peridotites, thereby altering the olivine and orthopyroxene to lizardite/chrysotile. During this process, certain mobile elements (e.g., Fe, Ni, and Mn) in the fluids and/or olivine serpentinization were incorporated into existing spinels, favoring ferritchromite growth. The second fluidic episode was dependent on prograde metamorphism at greater depths in the subduction channel, which not only produced antigorite via the recrystallization of lizardite and chrysotile, but also drove the fluids towards a more oxidizing state, thereby facilitating the growth of Cr-magnetites. Subsolidus re-equilibration controlled by fluid infiltration further promoted the inter-diffusion of elements (e.g., Cr, V, Co, Ni, Cu, and Zn) in the zoned spinels.

Mantle Flow Induced by the Interplay of Downgoing Slabs Revealed by Seismic Anisotropy Beneath the Sula Block in Eastern Indonesia

AbstractThe North Sulawesi subduction zone is characterized by southward subduction of the Celebes Sea slab to a depth of ∼250 km, mainly overlying the Sangihe slab that subducts west from the Molucca Sea and penetrates the mantle transition zone. The palaeo‐subducted Sula slab dips northward and partially underlies both the Sangihe and Celebes Sea slabs. Adjacent subduction zones with horizontal overlapping subducting slabs in the upper mantle have unclear dynamic interactions. An extensive strike‐slip fault forms the western boundary of the active North Sulawesi subduction zone, providing an ideal setting to study mantle flow between overlapping slabs. We use local S‐wave and teleseismic S and SK(K)S waveform splitting analysis to measure seismic anisotropy in the northern Sulawesi region. Our observations reveal typical mantle wedge corner flow within the Sangihe subduction system. In the Gulf of Tomini, the observed trench‐oblique fast‐axis orientations above the Celebes Sea slab are likely a consequence of the interaction between two subducting slabs. The southernmost measurement with an E–W‐trending fast direction in the mantle wedge might be related to the subduction of the Sula slab. Furthermore, fault‐parallel fast‐axis orientations of anisotropy near the southern segment of the Palu‐Koro fault are attributed to large‐scale shearing across this lithospheric‐scale strike‐slip fault system. Overall, our observations suggest that the strain caused by lithospheric and asthenospheric deformation is mainly confined within the microplate, displaying a restricted flow pattern and localized effects due to the size of the plate boundaries, such as the Palu‐Koro fault.