skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: A subduction influence on ocean ridge basalts outside the Pacific subduction shield
Abstract The plate tectonic cycle produces chemically distinct mid-ocean ridge basalts and arc volcanics, with the latter enriched in elements such as Ba, Rb, Th, Sr and Pb and depleted in Nb owing to the water-rich flux from the subducted slab. Basalts from back-arc basins, with intermediate compositions, show that such a slab flux can be transported behind the volcanic front of the arc and incorporated into mantle flow. Hence it is puzzling why melts of subduction-modified mantle have rarely been recognized in mid-ocean ridge basalts. Here we report the first mid-ocean ridge basalt samples with distinct arc signatures, akin to back-arc basin basalts, from the Arctic Gakkel Ridge. A new high precision dataset for 576 Gakkel samples suggests a pervasive subduction influence in this region. This influence can also be identified in Atlantic and Indian mid-ocean ridge basalts but is nearly absent in Pacific mid-ocean ridge basalts. Such a hemispheric-scale upper mantle heterogeneity reflects subduction modification of the asthenospheric mantle which is incorporated into mantle flow, and whose geographical distribution is controlled dominantly by a “subduction shield” that has surrounded the Pacific Ocean for 180 Myr. Simple modeling suggests that a slab flux equivalent to ~13% of the output at arcs is incorporated into the convecting upper mantle.  more » « less
Award ID(s):
1714892
PAR ID:
10354652
Author(s) / Creator(s):
; ; ; ; ;
Date Published:
Journal Name:
Nature Communications
Volume:
12
Issue:
1
ISSN:
2041-1723
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Delineation of geochemically distinct domains in Earth’s mantle is essential for understanding large-scale mantle convective flow and dynamics. Previous studies identify possible long-lived (>60 million-year) mantle isotopic domains (i.e. Antarctic-Zealandia, Pacific and Indian) near the Philippine Sea and western Pacific. Here we compile published basalt geochemistry of the Philippine Sea and surroundings and add new Mo isotopic and water content data for Gagua Ridge lavas, northwestern Philippine Sea, to distinguish slab-derived components during subduction. The water content, trace element, and Mo-Sr-Nd isotope compositions of Gagua Ridge arc lavas suggest that slab fluids and sediment melts are responsible for element recycling to the arc. The Philippine Sea basalts show both Indian and Zealandia-Antarctic Pb isotopic signatures; restoration of the basalt locations within a plate reconstruction shows the far-travelled Philippine Sea traversed these mantle domains. We establish the Indian mantle domain eastern boundary at ~120°E under SE Asia and the Indian Ocean. The Antarctic-Zealandia mantle domain lies south of ~10°N within the SW Pacific and has mostly remained in oceanic realms since ~400 Ma with only limited continental material input. 
    more » « less
  2. Abstract The Tonga‐Kermadec subduction zone exhibits the fastest observed trench retreat and convergence near its northern end. However, a paradox exists: despite the rapid trench retreat, the Tonga slab maintains a relatively steep dip angle above 400 km depth. The slab turns flat around 400 km, then steepening again until encountering a stagnant segment near 670 km. Despite its significance for understanding slab dynamics, no existing numerical model has successfully demonstrated how such a distinct slab morphology can be generated under the fast convergence. Here we run subduction models that successfully reproduce the slab geometries while incorporating the observed subduction rate. We use a hybrid velocity boundary condition, imposing velocities on the arc and subducting plate while allowing the overriding plate to respond freely. This approach is crucial for achieving a good match between the modeled and observed Tonga slab. The results explain how the detailed slab structure is highly sensitive to physical parameters including the seafloor age and the mantle viscosity. Notably, a nonlinear rheology, where dislocation creep reduces upper mantle viscosity under strong mantle flow, is essential. The weakened upper mantle allows for a faster slab sinking rate, which explains the large dip angle. Our findings highlight the utilizing rheological parameters that lead to extreme viscosity variations within numerical models to achieve an accurate representation of complex subduction systems like the Tonga‐Kermadec zone. Our study opens new avenues for further study of ocean‐ocean subduction systems, advancing our understanding of their role in shaping regional and global tectonics. 
    more » « less
  3. Abstract We investigated the state of the arc background mantle (i.e. mantle wedge without slab component) by means of olivine CaO and its Cr-spinel inclusions in a series of high-Mg# volcanic rocks from the Quaternary Trans-Mexican Volcanic Belt. Olivine CaO was paired with the Cr# [molar Cr/(Cr + Al) *100] of Cr-spinel inclusions, and 337 olivine+Cr-spinel pairs were obtained from 33 calc-alkaline, high-K and OIB-type arc front volcanic rocks, and three monogenetic rear-arc basalts that lack subduction signatures. Olivine+Cr-spinels display coherent elemental and He–O isotopic systematics that contrast with the compositional diversity of the bulk rocks. All arc front olivines have low CaO (0.135 ± 0.029 wt %) relative to rear-arc olivines which have the higher CaO (0.248 ± 0.028 wt %) of olivines from mid-ocean ridge basalts. Olivine 3He/4He–δ18O isotope systematics confirm that the olivine+Cr-spinels are not, or negligibly, affected by crustal basement contamination, and thus preserve compositional characteristics of primary arc magmas. Variations in melt H2O contents in the arc front series and the decoupling of olivine CaO and Ni are inconsistent with controls on the olivine CaO by melt water and/or secondary mantle pyroxenites. Instead, we propose that low olivine CaO reflects the typical low melt CaO of high-Mg# arc magmas erupting through thick crust. We interpret the inverse correlation of olivine CaO and Cr-spinel Cr# over a broad range of Cr# (~10–70) as co-variations of CaO, Al and Cr of their (near) primary host melts, which derived from a mantle that has been variably depleted by slab-flux driven serial melt extraction. Our results obviate the need for advecting depleted residual mantle from rear- and back-arc region, but do not upset the larger underlying global variations of melt CaO high-Mg# arc magmas worldwide, despite leading to considerable regional variations of melt CaO at the arc front of the Trans-Mexican Volcanic Belt. 
    more » « less
  4. Abstract The extensive fast seismic anomalies in the mantle transition zone beneath East Asia are often interpreted as stagnant Pacific slabs, and a reason for the widespread tectonics since the Mesozoic. Previous hypotheses for their formation mostly emphasize vertical resistances to slab penetration or trench retreat. In this study, we investigate the origin of these stagnant slabs using global‐scale thermal‐chemical models with data‐assimilation. We find that subduction of the Izanagi‐Pacific mid‐ocean ridge marked the transition of mantle flow beneath western Pacific from being surface‐driven Couette‐type flow to pressure‐driven Poiseuille‐type flow, a result previously unrealized. This Cenozoic westward mantle wind driven by the pressure gradient independently explains seismic anisotropy in the region. We conclude that the mantle wind is the dominant mechanism for the formation of stagnant slabs by advecting them westward while the pressure gradient holds them in the transition zone. 
    more » « less
  5. Abstract Subduction zones are fundamental features of Earth's mantle convection and plate tectonics, but mantle flow and pressure around slabs are poorly understood because of the lack of direct observational constraints on subsurface flow. To characterize the linkages between slabs and mantle flow, we integrate high‐resolution representations of Earth's lithosphere and slabs into a suite of global mantle convection models to produce physically plausible present‐day flow fields for Earth's mantle. We find that subduction zones containing wide, thick, and long slabs dominate regional mantle flow in the neighboring regions and this flow conforms to patterns predicted by simpler regional subduction models. These subduction zones, such as Kuril‐Japan‐Izu‐Bonin‐Mariana, feature prismatic poloidal flow coupled to the downgoing slab that rotates toward toroidal slab‐parallel flow near the slab edge. However, other subduction zones, such as Sumatra, deviate from this pattern because of the competing influence of other slabs or longer‐wavelength mantle flow, showing that upper mantle flow can link separate subduction zones and how flow at subduction zones is influenced by broader scale mantle flow. We find that the non‐linear dislocation creep reduces the coupling between slab motion and asthenospheric flow and increases the occurrence of non‐ideal flow, in line with inferences derived from seismological constraints on mantle anisotropy. 
    more » « less