Slab Subduction Research Articles

Comparison of Sn-related granitoids in subduction and collision settings by accessory mineral geochemistry: A case study in the Tengchong-Lianghe tin belt, SW China

Magmatic-hydrothermal Sn deposits can develop in both subduction and collision settings. However, a systematic comparison of the geochemical fingerprints and petrogenetic evolution of Sn-related granitoids from both settings is still lacking. The Tengchong-Lianghe tin belt hosts numerous Sn deposits but of different metallogenic ages, such as the Xiaolonghe and Lailishan Sn deposits, which provide an ideal study area for addressing this issue. Zircon U–Pb dating suggests that the Xiaolonghe and Lailishan Sn-related granitoids were emplaced at 74.2–75.5 and 50.0–51.1 Ma, respectively. Given the initial India-Asia collision commenced at ca. 65–60 Ma, these ages, consistent with available Sn metallogenic ages, suggest the formation of the Xiaolonghe and Lailishan Sn deposits in subduction and collisional settings, respectively. New zircon and apatite compositions, and apatite Sr isotopic data, combined with previously reported whole-rock geochemical results, all indicate that the Xiaolonghe and Lailishan Sn-related granitoids share similar magma sources, evolution processes, redox states, and volatile contents. The presence of hornblende, high apatite NdN/NdN* (mostly > 1), and whole-rock evolutionary trends, together indicate that they represent fractionated I-type granites. Distinctly elevated apatite 87Sr/86Sr ratios (0.70988–0.71430) suggest that they originated from an ancient lower crust. Subsequently, they underwent an extensive fractional crystallization dominated by feldspar, as revealed by whole-rock, apatite, and zircon elemental variation trends. The low zircon CeN/CeN* (mainly < 200) and ΔFMQ (–1.9 to 1.7), together with elevated apatite F content (2.30–3.69 wt%), reveal that they crystallized from reduced and F-rich magmas. The relatively high whole-rock zircon saturation temperatures (mainly > 800 ℃) and Ba/Pb ratios imply that they were produced at a higher temperature by biotite-dehydration melting, which requires additional heat from the mantle. The specific mechanism that triggers mantle upwelling (oceanic slab subduction or break-off) could be the most significant difference between Sn-related granitoids formed during subduction and collision processes.

Zircon and apatite U-Pb geochronology of the Paleoproterozoic (Eburnean) basement and late Neoproterozoic (Pan-African) metamorphism and magmatism from Port-Béni, Armorican Massif (France)

By re-examining the historical outcrops of Port-Béni located in the Trégor unit of the North Armorican Cadomian belt, the present work delivers four new ages that provide additional constraints on the Proterozoic history of northern Brittany. It is established that granitic, porphyritic rocks crystallized at the end of the Rhyacian (Paleoproterozoic), 2038 ± 12 Ma ago, before being transformed into orthogneisses at a late Neoproterozoic (Ediacaran) age of 621 ± 2 Ma, which is a minimum age, given the retrograde alteration these rocks underwent. The age of ca. 1.8 Ga previously proposed for the protolith of the Port-Béni orthogneiss should be discarded, and these two new ages are consistent with most of those yielded so far by the other Icartian (i.e., Eburnean) basement relics from the Armorican Massif. The gneissic basement was then intruded and disrupted into xenoliths by a granodioritic magma that crystallized 604.5 ± 2.0 Ma ago. This age, slightly younger than previously thought, corresponds to the emplacement age of one of the main units of the North Trégor batholith − the Pleubian-Talbert unit −, part of the Trégor volcano-plutonic complex, which may have built up over a longer period than that indicated by the uncertainty associated with this age. Caution should be exercised in extrapolating this age to that of the whole complex. Finally, doleritic dykes, possibly resulting in two swarms previously thought to be Paleozoic in age, have crosscut this complex. One of the latest yielded an age of 597 ± 15 Ma, indicating that the Trégor doleritic dyking episodes also occurred during the late Neoproterozoic, in between ca. 605 Ma and ca. 580 Ma. As the doleritic dykes are of tholeiitic composition, which distinguishes them from the earlier calc-alkaline magmas, they suggest that the intra-arc extension, documented in the southern, adjacent Saint-Brieuc unit of the belt, also affected the Trégor unit. They may likely have fed northern equivalents of the lava flows from the Paimpol Formation (exposed in between the Saint-Brieuc and the Trégor units), when magma production became moderately influenced by the Cadomian (i.e., Pan-African) subduction and mostly dominated by extension, possibly as a result of a steepening of a north-dipping subduction slab. Indeed, a re-examination of the available geochemical and geochronological data in the light of our new results documents that arc-magma production moved progressively from north (Trégor unit) to south (Saint-Brieuc unit) over time, in the interval 605-580 Ma.

Constraining subducting slab viscosity with topography and gravity fields in free-surface mantle convection models

Subducted slabs provide the primary driving force for mantle convection and plate tectonics. The strength of the slab directly controls the transfer of forces between it and the lithospheric plates at the surface. However, slab viscosity constrained by surface topography and gravity data is significantly lower than that suggested by mineral physics laboratory experiments. It is unclear whether the free-slip top boundary condition used by many studies affects the inverted slab viscosity with gravity or geoid data. In this study, we constrain the subducting slab viscosity structure by comparing the computed topography and gravity in slab subduction models with a free surface with observations. The free-surface model results support relatively weak slabs (20–120 times more viscous than the upper mantle) at the bending region, consistent with previous studies with a free-slip top boundary. The viscosity of the slab below the bending region barely affects the surface topography and gravity field, and both strong and weak slabs fit the observed topography and gravity field, suggesting that extra independent observations are needed to constrain the deep slab viscosity. We investigate the comprehensive relations between subduction interface viscosity, surface topography and gravity anomaly, and trench motion. Models with trench advance have significantly low topography and gravity above the volcanic arc, contradicting subduction zone observations. Together with present trench motion observations and previous studies, we support the idea that the trench retreats under normal single-slab subduction conditions.

Big Mantle Wedge and Intraplate Volcanism in Alaska: Insight From Anisotropic Tomography

AbstractWe determine high‐resolution tomographic models of isotropic P‐wave velocity (Vp) and tilting‐axis anisotropy of the Alaska subduction zone using a large number of local and teleseismic data recorded at many portable and permanent network stations in and around Alaska. We find a flat high‐Vp slab in the mantle transition zone (410–670 km depths) beneath western Alaska, which is connected with the subducting Pacific slab at 0–410 km depths, suggesting that a big mantle wedge has formed under western Alaska. Our tilting‐axis anisotropy model reveals complex mantle flows in the asthenosphere. Corner flow in the mantle wedge above the subducting Pacific slab and toroidal flow in the big mantle wedge are revealed, which may cause the Cenozoic intraplate volcanoes in western Alaska and the Bering Sea. In central Alaska, the mantle wedge beneath the Denali volcanic gap is characterized by high‐Vp and subhorizontal fast velocity directions normal to the volcanic arc, which may reflect a remnant of the subducted Yakutat slab. In SE Alaska, the shallow subduction of the Wrangell slab is visible above 150 km depth, and hot mantle upwelling through the Wrangell‐Yakutat slab gap may contribute to the Wrangell volcanic field.

Differential Uplift Triggered Basin‐And‐Range System: Evidence From Low‐Temperature Thermochronology in Eastern NE Asia

AbstractSince the Mesozoic, eastern NE Asia has experienced multiple tectonic events, resulting in a complex structure and forming one of the world's largest Meso‐Cenozoic lacustrine basin systems. Presently, basin evolution models require further elucidation regarding the simultaneous generation of diverse rift basins and the potential impact stemming from the closure of the Mudanjiang Ocean, whose oceanic closure demarcated the boundary between the Songliao Basin and the eastern basins, raises questions about its influence on the development of the basin‐and‐range system. To address these questions, we augment new low‐temperature thermochronological data on basement highs separating the eastern NE Asia basins to investigate the shallow‐deep coupling process of tectonic evolution since the Mesozoic. The new cooling age pattern shows non‐overlapping and spatial differences among major basement highs. Inverse thermal modeling revealed five‐stage cooling episodes among the basement highs, but with different onset and cooling rates of each episode, indicating a significant differential uplift mode. A major reburial stage was identified throughout eastern NE Asia during the mid‐Cretaceous. Compiling cooling age patterns and inverse thermal modeling, we note the existence of a proto‐basin covering an area much larger than the previously contemplated “Pan‐Sanjiang” Basin. In general, our study indicates the final closure of the Mudanjiang Ocean occurred at ca. 150–140 Ma. Since the Early Cretaceous, with changes in the subduction direction, two‐stage flat slab subduction of the Paleo‐Pacific plate and the consequent subduction of the Pacific plate co‐dominated the basements' differential uplift and the formation of the eastern NE Asia basin‐and‐range framework.

In situ study of the reaction phase A plus high-P clinoenstatite to forsterite plus water at reduced water activity

Abstract. We examined the reaction phase A plus high-P clinoenstatite to forsterite plus water (Reaction R1) by means of in situ X-ray diffraction measurements with the large volume press at the synchrotron PETRA III, Hamburg. Contrary to the study of Lathe et al. (2022), in which all experiments on Reaction (R1) were performed at a water activity of 1, the reversed experiments presented in this study were performed at reduced water activity with mole fractions of about XH2O= XCO2=0.5. The intention of this investigation was to test the observation made by Perrillat et al. (2005), which was that dehydration reactions are kinetically faster at reduced than under water-saturated conditions. The position of Reaction (R1) at the reduced conditions was determined by reversal brackets at 9.1 and 9.5 GPa (630 and 700 ∘C), at 9.7 and 10.0 GPa (725 and 700 ∘C), at 9.8 and 10.2 GPa (675 and 750 ∘C), and at 10.5 GPa (675 and 740 ∘C). Additionally, we performed two offline experiments with brackets at 10.0 and 10.6 GPa (750 and 700 ∘C, respectively) that are in agreement with the results of the in situ experiments. We do not observe any “intermediate” precursor phase in our experiments. The equilibrium of Reaction (R1) is shifted by about 100 ∘C to lower temperature compared to the results under water-saturated conditions. Thus, at a water activity (aH2O) below 1 the phase A plus clinoenstatite dehydration reaction can only occur in extremely cold subduction slabs. The kinetics of Reaction (R1) dehydration at reduced water activity is slower than that determined previously by Lathe et al. (2022) under water-saturated conditions. Thus, the above-mentioned hypothesis of Perrillat et al. (2005) could not be confirmed. However, in both of our studies on Reaction (R1), the newly formed dehydration product forsterite was of nanometer size, which supports earlier experimental observations, which is that product phases of dehydration reactions are generally very fine-grained and might promote the concept that intermediate-depth earthquakes in subduction zones are initiated by mechanical instabilities from extremely fine-grained materials formed during dehydration reactions.

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.

Potential of porphyry copper mineralization in the Triassic Yunxian-Jinggu arc, in the Sanjiang orogen, SW China

Porphyry copper deposits are commonly associated with porphyries generated in arc environments. Assessing the contemporaneous arc volcanic rocks provides crucial insights into the potential for porphyry copper mineralization. The Triassic Yidun arc in the northern part of the Sanjiang orogen in southwest China, which formed in association with the subduction of the Paleo-Tethys Oceanic slab, hosts several porphyry copper deposits. However, the potential for porphyry copper mineralization in the Yunxian-Jinggu volcanic arc in the southern part of the Sanjiang orogen, which was also linked to the subduction of the Paleo-Tethys Oceanic slab, remains unclear. This study reports zircon U–Pb ages, trace element compositions and Hf isotopes, and whole-rock element compositions and Sr–Nd–Pb isotopes of the basaltic-andesitic-dacitic rocks in the Yunxian-Jinggu volcanic arc. Zircon SHRIMP U–Pb dating of the Wenyu basalts and andesites yielded ages of 244.3 ± 1.9 Ma and 242.5 ± 2.5 Ma, respectively. The basalts have arc geochemical affinities with εHf (t) values of 2.3 to 9.3, εNd(t) values of –0.6 to 0.2, and (87Sr/86Sr)i values of 0.7065 to 0.7068. They have high Ba concentrations (∼ 550 ppm), highly variable Ba/La ratios (0.01 to 0.10), and a narrow range of Th/Yb ratios (1.2 to 1.9). These characteristics suggest that the basalts were derived from mantle metasomatized by fluids likely released from the subducted Paleo-Tethys slab. The Wenyu andesites, with εHf (t) values range from 1.6 to 7.0, εNd (t) value of –0.4, and (87Sr/86Sr)i value of 0.7031, along with distinct Pb isotopes, were probably produced through basaltic magmas undergoing crustal fractional crystallization. The Minle dacites, sharing similar εHf (t) value (–1.6 to 10.7) but lower εNd (t) values (–3.1 to –4.9) compared to Wenyu volcanics, are interpreted as derived from partial melting of basaltic magmas with contamination from ancient crustal components. The zircon Eu/Eu* ratios (mostly < 0.4), zircon ΔFMQ value (–0.6 ± 1.1), and bulk-rock Sr/Y ratios (< 10.3) for the Yunxian-Jinggu volcanic rocks indicate that co-genetic intrusions, if present, would likely be less hydrous and less oxidized, making them less favorable for porphyry copper mineralization.

Late Cretaceous basalts in the eastern South China Block: The result of slab tearing during retreating subduction of the Paleo-Pacific plate

Cretaceous arc-like, ocean island basalt (OIB)-like, and transitional basalts erupted on the eastern South China Block (SCB), and their petrogenesis was possibly related to the subduction of the Paleo-Pacific plate. However, the detailed petrogenesis in terms of slab subduction, magma sources, and the spatial distribution of the different types in the coastal and inland regions is poorly understood. Here we conduct a comparative study on the Late Cretaceous basalts from both coastal and inland regions. Our compiled age data demonstrate that the basalts in the coastal regions erupted in the period of 104–81 Ma. They show low Nb/La (<0.8), Nb/Y (<0.7), and other typical arc-like signatures, e.g., negative Nb-Ta anomalies, variable enrichment of LILE and LREE, and more radiogenic Sr-Nd isotopic compositions. In contrast, the basalts in the inland regions erupted at 86–71 Ma. They are characterized by high Nb/La (>0.9), slightly negative to positive Nb-Ta anomalies; they mostly possess positive εNd(t) values and geochemical features transitional to OIB-like basalts. The systematic geochemical differences may imply the spatial heterogeneity in their magma sources. The three end-member mixing model suggests that the isotopic signatures of arc-like basalts largely resulted from partial melting of the mantle wedge metasomatized by the subducted sediments, while those of transitional basalts likely resulted from depleted mantle source metasomatized by melts from altered oceanic crust (AOC) with minor sediments. We suggest that the retreating of the Paleo-Pacific plate promoted the melting of sediment-hybridized mantle wedge, and then triggered the opening of a slab window and asthenosphere-slab interactions.

Speculations on the Paleozoic legacy of Gondwana amalgamation

Using Gondwana as an example, we show how the geological record can be interrogated to detect significant changes in mantle convection patterns at critical junctures in Earth’s evolution. Evidence of major changes in mantle circulation in the aftermath of late Neoproterozoic-early Paleozoic Gondwana assembly is provided by widespread (i) plume-related magmatism around Gondwana’s periphery, (ii) ironstone deposits related to mantle plume-ocean ridge interaction and enhanced hydrothermal activity, and (iii) super-mature clastic deposits that reflect epeirogenic uplift triggered by mantle upwelling beneath Gondwana combined with deep tropical weathering.In our model, Gondwana assembled above a region of mantle downwelling in which subducted slabs between the converging Gondwanan continents descended to the core-mantle boundary. Renewed subduction along Gondwana’s periphery yielded early arc magmas. But as downwelling beneath Gondwana evolved into upwelling as a result of the ponding of subducted slabs at the base of the mantle, mantle plumes rose from the margins of the nascent upwelling to interact with the edges of Gondwana, where they penetrated the peripheral subduction zones via slab windows, tears and transform faults to generate voluminous calc-alkalic crustal melts in hydrated arc regions and A-type magmas in dry back-arc regions. The plumes also underplated oceanic lithosphere and interacted with adjacent ocean ridges, thereby enhancing hydrothermal activity and the flux of bioessential nutrients, leading to the recurrence of marine iron-rich sedimentary rocks in the geological record. At the same time, upwelling beneath a tropical to equatorial Gondwana led to epeirogenic uplift, deep weathering and erosion, resulting in the production of widespread super-mature clastic deposits. We contend that major changes in mantle convection patterns were encoded into the geological record of Gondwana assembly, influenced global-scale mantle convection patterns, and should be incorporated into geodynamic models for the assembly of Pangea.

The lithospheric cycle of ore-forming ingredients and the resultant recurrent ‘flare-ups’ of porphyry Cu deposits: Examples from the Yidun Arc, eastern Tibet

The porphyry Cu deposits with distinct mineralisation ages frequently occur as recurrent ‘flare-ups’ in a narrow region; however, the reasons controlling their occurrence are still debated. The eastern Yidun Arc in eastern Tibet produced a cluster of porphyry Cu deposits of the late Triassic and late Cretaceous ages, making it the ideal location for conducting this research. As a result, we use the late Triassic Xuejiping Cu–(Au) and late Cretaceous Hongshan–Hongniu Cu–Mo deposits as research objects, compiling new and published geochronological, elemental, and Sr–Nd–Hf isotopic data from their ore-forming porphyries to shed light on the above issue. The zircon U–Pb ages of the Xuejiping and Hongshan–Hongniu porphyries are 217.8 – 217.9 Ma and 78.7 Ma, respectively, and display adakitic characteristics. The Xuejiping porphyries present enriched Sr–Nd–Hf isotopes (0.70529 – 0.70590, −3.30 – −1.09, 0.40 – 2.20, respectively), low (Gd/Yb)N (2.24 – 3.20) and Nb/Ta (12.8 – 15.9) ratios, relatively high Nb (14.9 – 20.0 ppm) and Zr/Sm (23.4 – 41.5) values, indicating their origin in the partial melting of enriched mantle that is metasomatised by slab-melts and fluids in subduction zone. Furthermore, the Hongshan–Hongniu porphyry is adakitic, like the late Triassic igneous rocks in this area. Its Sr–Nd–Hf isotopes (0.69873 – 0.70996, −10.7 – −6.72, −14.9 – −0.50, respectively) are close to the isotopic evolution line of the late Triassic Xuejiping intrusions and show a trend towards the regional crystalline basement, which signify that the Hongshan–Hongniu porphyry mainly originated from the late Triassic juvenile crust, with some contributions from crystalline basement. One highly fractionated sample (i.e. HZK17-11-76) has extremely evolved Sr–Nd isotopes (0.70996 and −10.7, respectively), indicative of assimilation of sedimentary rocks. The new genetic models of the late Triassic and late Cretaceous porphyries in the Geza Arc show that the westward subduction of the Garzê–Litang oceanic crust in the late Triassic triggered the occurrence of the highly oxidised and Cu–S–H2O rich adakitic melts in the middle–lower crust, which evolved to amphibole cumulates and porphyry Cu–(Au) deposits (i.e. Xuejiping) in the upper crust. In the late Cretaceous, the extensional environment triggered the re-melting of late Triassic juvenile cumulates and Mo-rich crystalline basement, resulting in Cu–Mo–S–H2O-rich, highly oxidised magma. Through fractional crystallisation and assimilation of Mo-rich sedimentary rocks, the magma eventually transformed into Hongshan–Hongniu porphyry-(skarn) type Cu–Mo deposits. In summary, we conclude that the periodic recycling of the ore-forming ingredients on a lithospheric scale from slab subduction to post-subduction periods controls the recurrent ‘flare-ups’ of PCDs in this confined area.

Petrogenesis of Oligocene Granitoids and its Implications on Beryllium Mineralization at Pagele, Nianqingtanggula Mountain (Central Tibet)

AbstractThe newly discovered Oligocene granitoids (33.1–28.7 Ma) at Pagele are magmatic rocks related to beryllium mineralization during the India‐Asia late‐collisional stage. This discovery provides an ideal example to study the late‐collisional orogeny and beryllium prospecting in the Lhasa terrane. The Oligocene granitoids include porphyritic granodiorite, Stage I, II and III granites, and granitic pegmatite. Geochemical analysis shows that the porphyritic granodiorite is characterized by high SiO2, K2O, total ΣREE contents, and (La/Yb)N ratios; while the latter two by higher SiO2, lower ΣREE and (La/Yb)N ratios. Notably, the granitic pegmatite has extremely high Y/Ho, low K/Rb and Zr/Hf, and distinct REE tetrad effect (1.14–1.21). This study suggests that the porphyritic granodiorite may be derived from partial melting of beryllium‐rich materials composed of Lhasa ancient crust (70–80%) and enriched Lhasa lithospheric mantle (20–30%) under the tearing subduction of Indian slab. The three‐stage granites and granitic pegmatite, which contain higher beryllium contents or beryls, were likely generated by highly fractionation of the porphyritic granodioritic magma or other homologous magma. Considering the possible genetic and spatial link between Indian slab tear and rifts, we suggest that highly‐fractionated granites in rifting systems represent important Be prospecting targets in the Lhasa terrane.

The role of high oxygen fugacity on genesis of the late Cenozoic Abaga basalts in Xilin Gol League, Inner Mongolia, China

ABSTRACT The Abaga basalts, which erupted within an intraplate background in the late Cenozoic, are situated in eastern central Inner Mongolia. Previous studies have explored the genesis of these basalts, which were influenced by the subduction of slab. Although trace element ratios demonstrate that altered oceanic crust (AOC) and sediments have contributed to the genesis of the basalts, limited information about the detailed percentages of mixing among sediments, AOC, and depleted mantle (DM) based on Sr‒Nd‒Pb isotopic values has been reported. In this study, we modelled the possible mixing percentages among AOC, DM, and sediments. The calculated results illustrate that a proportion of 80% DM, 16% AOC, and 4% sediments is appropriate for interpreting the Sr‒Nd‒Pb isotopic variation of the basalts. In addition, the associated major and trace elements of olivine and the partition coefficient of V between olivine and melt are used to estimate the crystallization temperatures of the olivine phenocrysts and oxygen fugacity of the basalts. The calculated results display that the olivine mainly crystallized between 1248 and 1293°C (standard error estimate: 51°C), and the fO2 values of the Abaga basalts range from FMQ + 0.14 to FMQ + 1.46 (FMQ, fayalite–magnetite–quartz buffer). The estimated oxygen fugacity values plot in the range of island arc basalts (from FMQ to FMQ + 2) but are higher than those of both mid-ocean ridge basalt (MORB) (from FMQ-0.4 to FMQ + 0.5) and Abaga xenoliths (from FMQ-1.55 to FMQ + 0.53), which indicates that the basaltic melts were more oxidized than those of MORB and the Abaga lithospheric mantle. Furthermore, the primary magma and the corresponding mantle source were more oxidized than MORB. The Ba/Th and Ce/Pb ratios of the basalt present positive correlations with fO2 values, demonstrating that the addition of sediments and oceanic crust into the mantle source may have contributed to increasing the fO2 values of the basalts.

Middle to Late Triassic felsic and mafic magmatism in the eastern West Qinling Orogen, Central China: A record of the Paleo-Tethyan oceanic crust advance to retreat

The tectonic setting of Indosinian magmatism in West Qinling Orogen and the tectonic evolution of the Paleo-Tethys remain hotly debated. Here we present SHRIMP zircon UPb dating, major-trace elemental and Sr-Nd-Pb-Hf isotopic data for four granitic and gabbroic plutons across the eastern West Qinling Orogen, from north to south including Guanshan, Fengxian, Mishuling and Miba. Our new data indicate that the Fengxian granitoids, emplaced at 238 Ma in a continental arc setting, show adakitic affinity and were formed by partial melting of subducted sediments. The southern Guanshan granitoids, emplaced at 235–229 Ma in a continental arc setting, also show adakitic affinity but were generated by partial melting of thickened (>50 km) lower crust. The northern Guanshan granites and gabbros were emplaced at 217–216 Ma in a back-arc setting, with the former generated by partial melting of normal (40–50 km) lower crust and the latter derived from asthenospheric mantle. The Mishuling granitic pluton contains abundant mafic microgranular enclaves, both emplaced at 214 Ma in a back-arc setting with the host granitoids formed by partial melting of normal (40–50 km) lower crust and the enclave magmas derived from asthenospheric-lithospheric mantle. The Miba granitoids, emplaced at 216 Ma in a continental arc setting, show adakitic affinity and were generated by partial melting of thickened (>50 km) lower crust. We suggest that the Paleo-Tethyan oceanic crust had subducted near the Shangdan suture by 238 Ma and successively advanced northward, approaching southeastern Qilian Orogen at 235–229 Ma. Then slab rollback initiated and gradually migrated southward. The slab subduction did not cease before 214 Ma.

Isotopic evolutionary track of water due to interaction with rocks and its use for tracing water cycle through the lithosphere

Hydrogen and oxygen isotopic compositions of natural spa waters from hot and/or mineral-rich springs (including flowing/pumping wells) potentially provide important clues about water cycle in Earth's interior, because some of the waters are deep-seated and non-meteoric. However, the isotopic compositions of these waters are typically altered by contamination with meteoric water and are thus poorly known. To reconstruct their original (i.e., before these waters mix with meteoric water) isotope compositions, here we develop an isotopic evolution model of the lithosphere-enclosed waters and derive the “ocean-origin lithospheric water curve” (OLWC; δ2H = 60 / (δ18O – 11)) as a typical evolutionary track starting from ocean water. Then we apply them to spa waters from central Japan, a globally unique geotectonic region where double subduction of slabs occurs. The OLWC explains well observed isotope data of subseafloor pore waters, submarine mud volcano pore waters, coastal oil-field brines, and subduction zone volcanic steam. Mixing analysis indicates that lithospheric water accounts for 19–95 % of each spa water we investigated. The nonlinear relationships of the reconstructed water δ18O and Cl or Li concentrations justify our hypothesis that the lithospheric waters are undergoing evolution from ocean water to magmatic water by rock-water isotopic exchange and that they are not simply mixtures of the two. A comparison of the reconstruction results with model predictions under slab subduction settings suggests that most spa waters in the southern to central part of the study region contain progressively evolving lithospheric waters released from the subducting Philippine Sea slab, and that a few at western edge of the region contain waters released from the Pacific slab. The others in the northern and eastern part of the region contain waters at the early-to-mid evolutionary stage, which are discharged from Tertiary marine sedimentary rocks (including tuffs). All these waters were enclosed in the lithosphere for more than 1–5 Ma. Our approach provides useful tools for tracing sources of non-meteoric water and highlights water cycle through lithosphere for more than one million years.

Upper mantle melt caused by a subducted slab in the Indian-Eurasian continental subduction zone

A low-velocity layer atop the mantle transition zone has been extensively observed worldwide. In subduction zones, this layer is widely explained as partial melting triggered by slab subduction on a regional or global scale. However, direct observational evidence is still absent, and the response of the layer to slab subduction is not well known. Here, we image the seismic velocity around the mantle transition zone by matching synthetic and observed triplicated seismic P and sP waveforms in the Indian–Eurasian continental subduction zone. Our observations reveal a laterally varied low-velocity layer atop the mantle transition zone beneath the Hindu Kush, where a subducted slab extends to the mantle transition zone. The geometric morphology of the low-velocity layer documents that it is a partially molten layer induced by the subducted slab on a regional scale. Interestingly, our observations also support that the layer has a low viscosity. The decreased viscosity may facilitate slab motion in the deep mantle, contributing to slab stretching, tearing and break-off and its resultant rare recurring large intermediate-depth earthquakes in an intracontinental setting.

Evaluation of Empirical Ground-Motion Models for the 2022 New Zealand National Seismic Hazard Model Revision

ABSTRACT This article presents an evaluation of empirical ground-motion models (GMMs) for active shallow crustal, subduction interface, and subduction slab earthquakes using a recently developed New Zealand (NZ) ground-motion database for the 2022 New Zealand National Seismic Hazard Model revision. This study considers both NZ-specific and global models, which require evaluation to inform of their applicability in an NZ context. A quantitative comparison between the models is conducted based on intensity measure residuals and a mixed-effects regression framework. The results are subsequently investigated to assess how the models are performing in terms of overall accuracy and precision, as well as to identify the presence of any biases in the model predictions when applied to NZ data. Many models showed reasonable performance and could be considered appropriate for inclusion within suites of models to properly represent ground-motion predictions and epistemic uncertainty. In general, the recent models that are NZ-specific or developed on large international databases performed the best. This evaluation of models helped inform suitable GMMs for the ground-motion characterization model logic tree. In addition, spatial trends in systematic site-to-site residuals to the west of the Taupō Volcanic Zone demonstrated the need for backarc attenuation modifications for slab earthquakes.

Identification of Moho Discontinuity Depth and Subduction Slab in Bengkulu and South Sumatra Region using Receiver Function Method: A Preliminary Result

Bengkulu and South Sumatra regions are areas with high seismic activity due to faults and subduction zones. This study aims to analyze the depth of the Moho discontinuity layer and subduction slab under 4 BMKG seismic stations that form a perpendicular line of trench in the Bengkulu and South Sumatra areas. The distribution of stations used are Enggano Station (EGSI) in the front arc islands, Manna Station (MNAI) in the front arc basin, Pagar Alam Station (PKSI) in the volcanic-arc zone, and Lahat Station (LHSI) in the proto-volcanic zone. We selected teleseismic earthquake data with a distance of 30°-90° from the stations and magnitudes above 6 (M>6). For the identification of Moho and slab depths, we used analysis of receiver functions with iterative time domain deconvolution and migration of receiver functions with the AK-135 velocity model. The Oceanic Moho layer and the subduction slab under the EGSI Station were identified at depths of 34 km and 23 km respectively, under the MNAI Station the Moho layer and the subduction slab were identified at a depth of 18 km and 87 km respectively, under the PKSI Station the Moho layer and the subduction slab were identified at a depth of 34 km and 129 km respectively, and under the LHSI Station the Moho layer and subduction slab were identified at depths of 35 km and 175 km respectively.

Depth Estimation of Moho Discontinuity Layer at 5 BMKG Seismic Stations in East Java Using Receiver Function Method

East Java is one of the areas with a high level of seismicity due to the existence of a subduction zone in the south of Java Island. In addition, there are also active faults and active volcanos in East Java. This study aims to determine the depth of the Moho Discontinuity layer and subduction slab using receiver function method based on the iterative time domain deconvolution. This study was conducted using teleseismic earthquake data with an epicenter distance of 30°–90° from the receiver station and magnitude more than 6 (M≥6). 5 BMKG seismic stations that form a straight line are used, including those in the Bawean Arc (BWJI), Rembang Zone (BAJI), Kendeng Zone (SIJM), Modern Mountain Arc (PPJI), Southern Mountain Zone (GEJI). The depth of the Moho Discontinuity layer at BWJI Station was observed in ranges from 36–38 km, at BAJI Station it ranges from 39–40 km, at SIJM Station it ranges from 39–40 km, at PPJI Station it ranges from 46–48 km, and at GEJI Station it ranges from 33–36 km. In general, the Moho Discontinuity layer in the mountainous region is deeper due to the isostasy effect.

The Crustal Magmatic Structure Beneath the Denali Volcanic Gap Imaged by a Dense Linear Seismic Array

AbstractThe crustal structure in south‐central Alaska has been influenced by terrane accretion, flat slab subduction, and a modern strike‐slip fault system. Within the active subduction system, the presence of the Denali Volcanic Gap (DVG), a ∼400 km region separating the active volcanism of the Aleutian Arc to the west and the Wrangell volcanoes to the east, remains enigmatic. To better understand the regional tectonics and the nature of the volcanic gap, we deployed a month‐long north‐south linear geophone array of 306 stations with an interstation distance of 1 km across the Alaska Range. By calculating multi‐component noise cross‐correlation and jointly inverting Rayleigh wave phase velocity and ellipticity across the array, we construct a 2‐D shear wave velocity model along the transect down to ∼16 km depth. In the shallow crust, we observe low‐velocity structures associated with sedimentary basins and image the Denali fault as a narrow localized low‐velocity anomaly extending to at least 12 km depth. About 12 km, below the fold and thrust fault system in the northern flank of the Alaska Range, we observe a prominent low‐velocity zone with more than 15% velocity reduction. Our velocity model is consistent with known geological features and reveals a previously unknown low‐velocity zone that we interpret as a magmatic feature. Based on this feature's spatial relationship to the Buzzard Creek and Jumbo Dome volcanoes and the location above the subducting Pacific Plate, we interpret the low‐velocity zone as a previously unknown subduction‐related crustal magma reservoir located beneath the DVG.