Mantle Wedge Research Articles

Mineral chemistry and melt evolution of the mantle wedge peridotites in the Late Cretaceous Zagros Belt ophiolites (Iran): clues for the subduction initiation-induced forearc magmatism

We investigate in this paper mineral compositions and geochemical evolution of the mantle wedge peridotites preserved in the Late Cretaceous Zagros ophiolites of SW Iran. Mantle peridotites above subduction zones commonly experience distinct melting, depletion and refertilization processes as a result of the circulation of fluids derived from subducting slabs and flux melting. Our results reveal that the mantle wedge peridotites in the Zagros ophiolites are characterized mainly by residual and impregnated types. Residual peridotites resulted from early depletion and later refertilization processes, whereas impregnated peridotites developed due to episodic melt impregnations within and across the mantle. Mg#s and NiO contents, spinel Cr#, Mg#, and TiO 2 in olivines, Mg# and Al 2 O 3 contents of orthopyroxenes, and Mg#, TiO 2 and Al 2 O 3 contents in clinopyroxenes of dunites, harzburgites and lherzolites indicate the significant role of re-equilibration processes among different mineral phases and interactions with basaltic melts percolating within the host peridotites. The observed geochemical variations in the mineral chemistry of the Zagros peridotites reflect changes in magma chemistry and fluctuations in the degree of melt extraction and melt–rock interactions within the mantle peridotites. However, our data suggest that Mg–Fe distribution in the spinels of some dunites and harzburgites might also have resulted from subsolidus redistribution and exchange with surrounding olivine grains. Spinel and clinopyroxene phases in gabbroic rocks and ultramafic cumulates within the Zagros ophiolites also show significant variations in their compositions, suggesting that their magmas evolved from MORB-like to IAT, calc–alkaline and boninite suites, typical of subduction initiation-generated melts. Hence, the Zagros ophiolites present a case study of time-progressive melt evolution of the forearc oceanic lithosphere. Supplementary material: Datasets of the mineral chemistry and melt evolution of the mantle wedge peridotites in the Late Cretaceous Zagros Belt ophiolites (Iran) are available in tabular and figurative form at https://doi.org/10.6084/m9.figshare.c.7093528 . Thematic collection: This article is part of the Ophiolites, melanges and blueschists collection available at: https://www.lyellcollection.org/topic/collections/ophiolites-melanges-and-blueschists

Carbon-rich polyphasic inclusions in postcollisional mafic magmatic rocks from the Dabie Shan, China: Implications for the carbon cycle in continental subduction zones

Subduction is a fundamental geodynamic process that transfers carbon from Earth’s surface into the mantle. However, current understanding of the migration mechanisms, final storage region, and species involved in carbon recycling from continental crust remains limited. Here, we investigated the compositions of polyphasic inclusions and Mg isotopes in postcollisional mafic magmatic rocks from the Dabie Shan region of China. The main rock-forming minerals contained two distinct types of polyphasic inclusions, which displayed systematic differences in daughter mineral/gaseous phase assemblages, including host-like silicates ± carbonates (magnesite, dolomite, and calcite) + CH4 and carbonates + talc ± SiO2 (aqueous) + CH4, respectively. These inclusions indicate that carbon-rich silicate melts and carbon-rich magmatic fluids were trapped by host minerals during magmatic processes. The abundant carbonates and CH4 in both types of inclusions suggest that the mantle source of these postcollisional mafic magmatic rocks was rich in carbon, most likely existing in the forms of CO2 and CH4. Moreover, the studied postcollisional mafic magmatic rocks have mantle-like Mg isotope compositions, with δ26Mg values ranging from −0.23‰ to −0.16‰. The combined observations of polyphase inclusions and Mg isotopes indicate that a substantial carbon-rich mantle domain arose from the metasomatism of silicate melts derived from subducted continental slabs that had dissolved a certain quantity of CO2 and CH4. We proposed that continental subduction is an efficient pathway for transporting crustal carbon into an orogenic subcontinental lithospheric mantle wedge, where the recycled carbon can be stored for >100 m.y. and eventually released to the surface during postcollisional magmatism.

Detection of a recent large Hyuga-nada long-term slow slip event and estimation of its spatiotemporal slip distributions

Using horizontal and vertical GNSS time series data from the GSI in sotheastern Kyushu from January 1, 2017 to June 30, 2022, we detected a recent long-term slow slip event (L-SSE) that occurred in the Hyuga-nada region, southwest Japan, and estimated its spatiotemporal slip distribution. We performed such analysis considering the piecewise evolution of slip over time with time windows of 0.8 years to ensure a good signal-to-noise ratio in the horizontal displacements in each time window. The results showed that a slip of more than 2 cm occurred in the central part of Miyazaki Prefecture from 2018.5 to 2019.3 (the units were 0.1 years = 36.5 days). Then, the amount of slip increased and expanded slightly to the southern part of Miyazaki Prefecture from 2019.3 to 2020.1, and the amount of slip reached a maximum of 3.9 cm in the subsequent period (2020.1–2020.9). A smaller slip occurred at almost the same location in the following 0.8 years. Therefore, the duration of the L-SSE was approximately 3.2 years from approximately 2018.5 to 2021.7. The annual average maximum slip rate was approximately 4.9 cm/yr (3.9 cm/0.8 yr) during the period from 2020.1 to 2020.9. The maximum total slip was estimated to be approximately 12.9 cm, and the equivalent release moment was 4.9 × 1019 Nm, corresponding to Mw7.1. Compared to previous L-SSEs that occurred in the Hyuga-nada region, the annual average maximum slip rate was relatively low, and the main total slip was estimated at almost the same location on the plate interface. On the other hand, the slip duration of this L-SSE was the longest, and the release moment and moment magnitude were the greatest. The total slip area of more than approximately 10 cm in the Hyuga-nada L-SSE estimated in this study almost overlapped with the afterslip area of the December 3, 1996 Hyuga-nada earthquake, with a depth range of approximately 30–40 km. Short-term SSEs also occurred in the same depth range. These coincidental same depth ranges of interplate seismic events stem from large thermal gradients. The large thermal gradients play an important role in narrowing the depth range of frictional parameters (a–b) at the plate interface. In addition, the pore pressure and normal stress are closely related to the critical stiffness. In particular, contrasting values of pore pressure and normal stress at the shallow and deep sides of the mantle wedge corner may also contribute to transient aseismic slip within the same depth range.

Neoarchean and Paleoproterozoic magmatism in the Longjiang area, NE China: Implications for the recycling history of supracrustal components in the eastern Central Asian Orogenic Belt

Substantial Phanerozoic crustal growth and reworking in the eastern Central Asian Orogenic Belt has partly obscured the Precambrian geological record. Accurate identification of such Precambrian magmatism and its petrogenesis is paramount to unveiling the early crustal evolution of the eastern Central Asian Orogenic Belt. Two outcrops of Neoarchean and Paleoproterozoic magmatism in the Longjiang area of NE China, eastern Central Asian Orogenic Belt, were recently discovered. The Kaoquantun outcrop consists of monzogranites that yield a zircon 207Pb/206Pb age of 2568 ± 8 Ma and are similar to peraluminous A1-type granites. They were derived from the partial melting of a hybrid crustal source comprising newly accreted Mesoarchean crust and abundant supracrustal metapelites in a post-collisional extensional setting. Such metapelites were initially part of the upper crust but were buried in the middle−lower crust due to crustal thickening prior to ca. 2.6 Ga. In the Zhanbeitun outcrop, basaltic andesites overlie monzogranites, and both have been intruded by syenogranites. The monzogranites and syenogranites yield zircon 207Pb/206Pb ages of 1881 ± 10 Ma and 1843 ± 3 Ma, respectively, whereas the basaltic andesites are inferred to have formed at 1.88−1.84 Ga. The 1.88 Ga monzogranites have relatively enriched zircon Hf-O isotopic compositions and were derived from a hybrid melt of mafic rocks in the lower crust and ancient metabasalts that had undergone low-temperature alteration at supracrustal depths. The 1.88−1.84 Ga basaltic andesites are Nb-rich and were generated by the partial melting of a mantle wedge that was metasomatized by subduction-derived fluids, recycled sediments, and slab-derived melts. The 1.84 Ga syenogranites are peraluminous and K-rich, and record a profound mixing and assimilation process involving melts derived from supracrustal sediments and the lower crust. Three stages of Paleoproterozoic magmatism emerged from subduction-related crustal thickening, slab breakoff, and collision, respectively. The zircon Hf isotopic data further indicate that the eastern Central Asian Orogenic Belt underwent at least one stage of proto-crustal growth at ca. 3.0 Ga, followed by multiple phases of crustal reworking during the Neoarchean and Paleoproterozoic, with dominant involvement of supracrustal components.

Electrical conductivity of antigorite–olivine aggregates under high temperature and pressure: Implications for the water content in mantle wedges

Antigorite–olivine aggregates were synthesized to serve as representative lithologies of subduction zones, which often suffered serpentinization due to the hydration of the mantle wedge. The electrical conductivities (EC) of antigorite–olivine aggregates with different antigorite proportions before and after dehydration at 2.0 GPa were measured. An increase in electrical conductivity prior to antigorite dehydration was due to Fe and Mg interdiffusion with an activation enthalpy of 130–200 kJ/mol. The electrical conductivity of the aggregates increased to 0.1 S/m due to antigorite dehydration. The electrical conductivity of 0.1 S/m derived from the laboratory experiments is consistent with geophysical observations. Combined cube model and modified Archie's law model with previous studies, the volume fractions of aqueous fluids in the mantle wedge of cold and hot subduction zones were estimated.

An enhanced chemical maturation of juvenile arc crust recorded by the Urad Houqi intrusive suite in Langshan arc, Central Asia

Continental arc magmatism is crucial in producing juvenile continental crust of andesitic to dacitic composition, yet its impact on the formation of the modern mature continental crust, which is more enriched in potassic and granitic components, remains poorly understood. We examined a suite comprised of norite gabbro, gabbroic diorite, high-Mg diorite, dioritic enclave, and medium-K to high-K granites (the Urad Houqi suite) from the Langshan arc in the southern Central Asian Orogenic Belt (CAOB). Zircon U-Pb geochronology indicates that the Urad Houqi suite was formed ca. 278−268 Ma by a continuous magmatic event, marking the final episode of the late Paleozoic continental arc magmatism in the southern CAOB induced by the subduction of the Paleo-Asian Ocean. The suite exhibits a more potassic and fertile composition compared to typical continental arc magmas, and shows a continuous increasing K2O/Na2O ratio and incompatible element concentrations, including K, Rb, Th, and U, with increasing SiO2 content. The geochemical evolution of the suite offers valuable insights into the maturation process of the juvenile arc crust. By combining zircon Hf and bulk-rock Sr-Nd-Pb isotope studies, we demonstrated that the various lithologies within the Urad Houqi suite recorded a three-stage enhancement of the geochemical fertility of the juvenile Langshan arc crust: (1) generation of the initially fertile primary melt from the mantle wedge metasomatized by the recycled subducted sediments; (2) trans-crustal open-system hybridization processes, including peritectic reaction and magma mixing, that further fertilized the derived melts in addition to fractional crystallization; and (3) diverse petrogenetic processes that contributed to the fertility of the granitic magmas, including polybaric fractionation of the primary magma, anatexis of ancient crust, and extraction of high-silica melts from the shallow mush system. The three-stage maturation of the juvenile Langshan arc crust was typically coupled with the accelerated subduction of the Paleo-Asian Ocean since the earliest Permian. Our study finds that the Urad Houqi suite can serve as an example of modern continental crust maturation at continental arc settings.

Regional-Scale Paleoproterozoic Heating Event on Archean Acasta Gneisses in Slave Province, Canada: Insights from K–Ar and 40Ar/39Ar Chronology

Slave Province in Canada is an Archean granite–supracrustal terrane at the northwestern corner of the Canadian Shield. It is bordered by the Thelon–Taltson orogen (2.0 to 1.9 Ga) to the southeast and the Wopmay orogen (1.9 to 1.8 Ga) to the west. Acasta gneisses, exposed in the westernmost Slave Province, and the Wopmay rocks, located close to the gneisses, were systematically collected for K–Ar and laser step-heating 40Ar/39Ar single-crystal analyses of the biotite and amphibole. The K–Ar biotite ages of the four Wopmay samples range from 1816 ± 18 Ma to 1854 ± 26 Ma. The 40Ar/39Ar biotite analyses of the three Wopmay samples yield plateau ages of 1826 ± 21 Ma, 1886 ± 13 Ma, and 1870 ± 18 Ma. These ages fall within the reported U–Pb zircon age range of the Wopmay orogen. The K–Ar biotite ages of the fifteen Acasta gneisses range from 1779 ± 25 Ma to 1877 ± 26 Ma, except for one younger sample (1711 ± 25 Ma). The 40Ar/39Ar analyses of the biotite crystals from three samples give the plateau ages of 1877 ± 8 Ma, 1935 ± 14 Ma, and 1951 ± 11 Ma. The K–Ar amphibole ages from twelve samples range from 1949 ± 19 Ma to 1685 ± 25 Ma. Two samples of them give ages older than the zircon U-Pb age of Hepburn plutons. The 40Ar/39Ar analyses of the amphibole crystals show varied age relations. The two samples give plateau ages of 1814 ± 22 Ma and 1964 ± 12 Ma. Some samples exhibit apparent old ages of ~2000 Ma in the middle temperature fractions. These old fractions result from the amphibole crystals, originally formed in the Archean, being affected by the thermal events during the Wopmay orogeny but not fully resetting. These observations suggest that the K–Ar system ages of the biotite and amphibole in the Archean Acasta gneiss were rejuvenated during the Paleoproterozoic ages. The Discussion explores the possibility that the heat source rejuvenating the K–Ar system ages may have arisen due to asthenospheric extrusion into the wedge mantle, a process likely triggered by subduction rollback.

Oligocene melting of subducted mélange and its mantle dynamics in northeast Asia

Melting of subducted mélange can potentially transport mass from the slab-mantle interface to the mantle wedge in subduction zones. The mélange diapir model was primarily proposed from the results of laboratory experiments and thermodynamic modeling. However, the melting mechanisms of mélange diapirs in subduction zones remain unclear. To further constrain the mantle dynamics of a mélange diapir, we studied Oligocene alkaline intermediate rocks on the northeast Asian continental margin. We report whole-rock geochemical and Sr-Nd-Pb-Mg-Zn isotope data and show that these rocks formed by partial melting of mélange. We conclude that a diapir was the mechanism for Oligocene melting of the mélange. We also identified younger rocks formed by melting of mélange in the eastern part of northeast Asia, implying an eastward shift in such magmatism since the Oligocene. Our results and the tectonic setting indicate that melting of mélange diapirs occurred preferentially during tectonic transitions, such as the formation of a back-arc basin triggered by trench-perpendicular mantle flow. The low-viscosity mantle with an incompressible stress field triggered melting of the mélange diapirs. Interactions occurred between the mélange diapirs and carbonated peridotites, constraining the depth of mélange-mantle interactions to the asthenosphere, which is deeper than the depth inferred in previous studies.

Multilayer anisotropy along the Alaska-Aleutians Subduction zone

SUMMARY Increasing evidence from seismic methods shows that anisotropy within subduction zones should consist of multiple layers. To test this, we calculate and model shear wave splitting across the Alaska-Aleutians Subduction Zone (AASZ), where previous studies have argued for separate layers of anisotropy in the subslab, slab and mantle wedge. We present an updated teleseismic splitting catalogue along the span of the AASZ, which has many broad-band seismometers recently upgraded to three components. Splitting observations are sparse in the Western Aleutians, and fast directions are oriented generally trench parallel. There are significantly more splitting measurements further east along the AASZ. We identify six regions in the Central and Eastern Aleutians, Alaskan Peninsula and Cook Inlet with a high density of splits suitable for multilayered anisotropy analyses. These regions were tested for multilayer anisotropy, and for five of the six regions we favour multiple layers over a single layer of anisotropy. We find that the optimal setup for our models is one with a dipping middle layer oriented parallel to palaeospreading. A prominent feature of our modelling is that fast directions above and below the dipping layer are generally oriented parallel to the strike of the slab. Additionally, we lay out a framework for robust and statistically reliable multilayer shear wave splitting modelling.

Volatiles in the mantle transition zone and their effects on big mantle wedge systems.

Volatiles in the mantle transition zone and their effects on big mantle wedge systems.

Ultramafic Rocks from the Sanbagawa Belt: Records of Mantle Wedge Processes

Mantle wedge domains beneath the forearc Moho are unique regions of Earth’s interior where mantle encounters subducting oceanic plates. Crystal-plastic deformation and fluid-induced reactions in the supra-subduction mantle control global material circulation, arc volcanism, and seismicity within subduction zones. The Sanbagawa metamorphic belt contains numerous ultramafic blocks in its higher-grade zones, some of which likely originated as lower crustal arc cumulates that were subsequently incorporated into the mantle wedge and transported to the slab–mantle interface by mantle flow. Properties of these ultramafic rocks provide a valuable opportunity to understand the dynamic processes of the mantle wedge up to 80 km depth, including mantle flow, hydration/dehydration, and fluid–rock interactions near the slab–mantle interface of a warm subduction zone.

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.

The Role of Slab Remnants in Modulating Free Subduction Dynamics: A 3‐D Spherical Numerical Study

AbstractSeismic tomography of Earth's mantle images abundant slab remnants, often located in close proximity to active subduction systems. The impact of such remnants on the dynamics of subduction remains underexplored. Here, we use simulations of multi‐material free subduction in a 3‐D spherical shell geometry to examine the interaction between visco‐plastic slabs and remnants that are positioned above, within and below the mantle transition zone. Depending on their size, negatively buoyant remnants can set up mantle flow of similar strength and length scales as that due to active subduction. As such, we find that remnants located within a few hundred km from a slab tip can locally enhance sinking by up to a factor 2. Remnant location influences trench motion: the trench advances toward a remnant positioned in the mantle wedge region, whereas remnants in the sub‐slab region enhance trench retreat. These motions aid in rotating the subducting slab and remnant toward each other, reducing the distance between them, and further enhancing the positive interaction of their mantle flow fields. In this process, the trench develops along‐strike variations in shape that are dependent on the remnant's location. Slab‐remnant interactions may explain the poor correlation between subducting plate velocities and subducting plate age found in recent plate tectonic reconstructions. Our results imply that slab‐remnant interactions affect the evolution of subducting slabs and trench geometry. Remnant‐induced downwelling may also anchor and sustain subduction systems, facilitate subduction initiation, and contribute to plate reorganization events.

The Lithosphere and Upper Mantle of the Western‐Central Mediterranean Region From Integrated Geophysical‐Geochemical Modeling

AbstractThis study integrates geophysical‐geochemical data to investigate the thermochemical structure of the lithosphere and sublithospheric mantle, along the Southern Tyrrhenian Basin, Apennines, Adriatic Sea, Dinarides, and Carpathians‐Balkanides. We present the lithospheric structure of the Adria microplate and the two opposing mantle slabs along its NE and SW margins. The modeling shows the presence of two asthenospheric mantle wedges aligning with the Apenninic and Dinaric continental mantle slab rollback, along with cold (−200°C) sublithospheric anomalies beneath Adria's NE and SW margins. In the northern Adria region, the lithosphere undergoes synchronous thinning in the Tyrrhenian domain and thickening toward the forefront of the northern Apennines. This is associated with the northeastward rollback of the SW Adriatic slab, leading to subsequent delamination of the continental mantle. In the southern Adria region, the complex deep structure results from the variably oriented lithospheric slabs, and nearly 90‐degree shift of the tectonic grain between the southern Apennines and the Calabrian Arc. At the SW Adria margin, beneath the northern Apennines, the thermal sublithospheric anomaly is attached to the shallower lithosphere, while a slab gap is modeled in the southern Apennines. One possibility is that the gap is due to a recent horizontal slab tear. Along the NE margin of Adria, the thermal anomaly penetrates to depths of about 200 km in the northern Dinarides and 280 km in the southern Dinarides, shallower than the SW Adria anomaly, which extends to at least 400 km depth.

Triassic oceanic–continental transition of the Ganzi–Litang Ocean basin: Evidence from high-Mg andesites and Nb-enriched basalts in the Litang area, Eastern Tibet

Petrological evidence of the oceanic–continental transition of intra-oceanic arcs enables the identification and reconstruction of the transformation of ocean basins into continents. In this study, we present new elemental and Sr–Nd isotopic compositions and zircon U–Pb ages for Triassic andesites and basalts in the Litang area, Eastern Tibet. Andesites in Cunge and basalts in Juewu contain zircon grains with U–Pb ages of 246 ± 2 and 248 ± 2 Ma, respectively, indicating a Middle Triassic origin. Andesites in Cunge have high SiO2 contents (53.1–56.5 wt%), MgO contents (5.15–7.13 wt%), and Mg# values (54–61), similar to Sanukite-type high-Mg andesites. Basalts in Juewu have high TiO2 (1.43–1.80 wt%) and Nb (12.5–20.2 ppm) contents, and high Nb/U and (Nb/La)N values, typical of Nb-enriched basalts. The andesites and basalts were both characterized by low (87Sr/86Sr)i values (0.7033–0.7044) and positive εNd(t) values (+4.2–+4.8). Therefore, interactions between melt from the subducting slab and mantle wedge peridotite formed andesites in Cunge, whereas the partial melting of mantle peridotite metasomatized by an adakitic melt (produced by the melted slab) formed basalts in Juewu. The findings of this study indicate that the formation of the andesites and basalts in the Litang area were related to intra-oceanic subduction magmatism, which was the product of the westward subduction of the Ganzi–Litang Ocean basin in the Middle Triassic.

Recycling of crustal Os in nascent subduction zones revealed by ophiolitic boninites

Olivine-rich lavas in the Margi area, Cyprus, are part of the boninitic Upper Pillow Lava series of the Troodos ophiolite. Re-Os isotope compositions of both whole rock and olivine phenocrysts of these lavas have been measured to determine the initial 187Os/188Os isotope compositions of the primary boninitic magma. The Margi lavas have high and variable whole rock Re-Os isotope compositions (187Re/188Os: 1.060–8.613; 187Os/188Os: 0.1307–0.1599), with (87Sr/88Sr)90Ma (0.705936 ± 0.001021), εNd(90 Ma) (5.6 ± 2.5), and δ18O (8.59 ± 0.38 ‰). In contrast, olivine phenocrysts, which host spinel and melt inclusions, have low and constant 187Os/188Os from 0.1301 to 0.1336 and low 187Re/188Os from 0.1350 to 1.557. Initial 187Os/188Os of these olivine phenocrysts (0.1300 ± 0.0041) is lower than that of whole rock (0.1373 ± 0.0135). These facts indicate that the olivine phenocrysts faithfully record initial 187Os/188Os of the primary boninitic magmas, which are comparable to the 187Os/188Os of primitive mantle, whereas whole rock Re-Os isotopes have been affected by secondary processes after olivine crystallization. Further geochemical modeling indicates that, although Os can be mobilized during slab dehydration, the contribution of crustal Os is insufficient to significantly raise the 187Os/188Os of the mantle wedge in a nascent subduction zone. Overall, our findings provide new insights about Os recycling in the early stages of subduction.

Melting behavior of impure limestone under H2O-poor conditions: Implications for the contribution of carbonate-rich sediments to arc magmatic carbon output

Sedimentary carbon may be transferred from subducting slabs to the overlying mantle wedge as a component of sediment-derived melts and contribute to carbon output in subduction zones. To investigate the melting behavior of subducting carbonate-rich sediments at shallow depths and evaluate the consequent efficiency of carbon recycling, this study presents the phase relations of H2O-poor impure limestone at 1.5–4.5 GPa and 750–1300 °C. The onset of melting occurs between 750 and 800 °C at 1.5–2.5 GPa, between 950 and 1000 °C at 3.5 GPa, and between 1000 and 1050 °C at 4.5 GPa. The melting experiments yield silicate melts with CO2 contents below 5 wt% and H2O contents below 10 wt%, exhibiting compositional evolution from rhyolitic to dacitic and trachytic as the temperature increases. Combining these results with previous experimental studies on various sediment systems, we conclude that the position of the solidus and melt compositions are influenced by multiple factors, including pressure, bulk composition, and volatile concentrations. Calculations demonstrate that the subduction of carbonate-rich sediment leads to a near-surface recycling efficiency for carbon below 5%, notably lower than the recycling efficiency of water (>75%). This study also indicates that sediment melting in subarc regions is restricted to either hot subduction zones or to subducted slabs at subarc depths exceeding 150 km, which are uncommon among global subduction zones on the modern Earth.However, sediment diapirism may be an efficient way for subducting sedimentary carbon to migrate into source regions of arc magmas. Reaction experiments involving hydrous impure limestone and lherzolite were conducted at 2.5 GPa and 1000 °C, as well as 4.5 GPa and 1100 °C to constrain the melting behavior of carbonate-rich sediment diapirs in the mantle wedge. Silicate melts and a thin websterite reaction layer are generated at 2.5 GPa, while at 4.5 GPa carbonate melts are produced, resulting in the thorough metasomatism of lherzolite. These findings highlight distinct and diverse contributions that carbonate-rich sediment diapirs make to compositional heterogeneity in the mantle wedge at different depths, and also to carbon degassing in magmatic arcs.

Sulfur disproportionation in deep COHS slab fluids drives mantle wedge oxidation.

Sulfur degassed at volcanic arcs calls for dissolved sulfate ions (S6+) released by subduction-zone fluids, oxidizing (in association with carbon) the subarc mantle, but sulfur speciation in subduction fluids at subarc depths remains unclear. We apply electrolytic fluid thermodynamics to model the dissolution behavior of pyrite in metacarbonate sediments as a function of P, T and rock redox state up to 4.3 gigapascals and 730°C. At subarc depth and the redox conditions of the fayalite-magnetite-quartz oxygen buffer, pyrite dissolution releases oxidized sulfur in fluids by disproportionation into sulfate, bisulfite, and sulfide species. These findings indicate that oxidized, sulfur-rich carbon-oxygen-hydrogen-sulfur (COHS) fluids form within subducting slabs at depths greater than 100 kilometers independent from slab redox state and that sulfur can be more effective than the concomitantly dissolved carbon at oxidizing the mantle wedge, especially when carbonates are stable.

Local-S shear wave splitting along the length of the Alaska–Aleutian subduction zone

SUMMARY The Alaska–Aleutian subduction zone represents an ideal location to study dynamics within a mantle wedge. The subduction system spans several thousand kilometres, is characterized by a slab edge, and has ample seismicity. Additionally, the majority of islands along the arc house broad-band seismic instruments. We examine shear wave splitting of local-S phases originating along the length of the subduction zone. We have dense measurement spacing in two regions, the central Aleutians and beneath Alaska. Beneath Alaska, we observe a rotation in fast splitting directions near the edge of the subducting slab. Fast directions change from roughly trench perpendicular away from the slab edge to trench parallel near the boundary. This is indicative of toroidal flow around the edge of the subducting Alaska slab. In the central Aleutians, local-S splitting is primarily oriented parallel to, or oblique to, the strike of the trench. The local-S measurements, however, exhibit a depth dependence where deeper events show more consistently trench-parallel directions indicating prevalent trench-parallel mantle flow. Our local-S shear wave splitting results suggest trench-parallel orientation are likely present along much of the subduction zone excited by the slab edge, but that additional complexities exist along strike.

Petrogenesis and tectonic implications of late Permian andesites from Rongzharinian, northern Qiangtang, Tibet

At convergent margins, high-magnesian andesitic rocks (HMAs) can offer vital information about the formation of the continental crust and the interactions between the deep crust and mantle. This study focuses on new identified HMAs from the Rongzharinian area, northern Qiangtang Terrane (NQT). We report the major, trace elemental, Sr-Nd-Hf isotopic data, and zircon U-Pb age of volcanic rocks with the aim of understanding their petrogenesis and tectonic setting. A concordant age of 257.7 ± 1.5 Ma for the Rongzharinian high-Mg andesites is determined by zircon U-Pb dating using LA-ICP-MS techniques. The Rongzharinian HMAs have high MgO contents (5.88–7.39 wt%), Mg# values (63–68), and Cr (143–185 ppm) and Ni (48–66 ppm) concentrations, close to primitive HMAs. They underwent variable degrees of fractionation of amphibole and minor clinopyroxene and negligible crustal contamination during magma evolution. REE modeling results show that the Rongzharinian HMAs can be generated by 4 %–7% melting of hydrous spinel lherzolite at relatively shallow depths. Geochemical data indicate that the Rongzharinian HMAs were dominantly derived from the mixed mantle wedge source composed of N-MORB-like and E-MORB-like components. Based on the regional tectonic evolution, we propose that the Late Permian rollback of the Paleo-Tethyan slab may induce upwelling of the asthenosphere mantle, resulting in melting of the mantle source and the generation of the Rongzharinian high-Mg andesites.