Control of slab tears and slab flat wedging on volcanism in the Alaska subduction zone

Multistage plate subduction plays a crucial role in magmatism; however, the mechanisms by which deep geodynamic processes govern volcanism in the Alaska subduction zone remain controversial. Using numerous travel-time data from several seismic arrays, we constructed high-resolution tomographic models to investigate the velocity structure of the Pacific Plate and Yakutat slab. Our tomographic results revealed high-velocity anomalies in the Pacific Plate and Yakutat slab, while the low-velocity areas within the Pacific Plate were identified as slab tears. We suggest that the Pacific Plate transitioned from oblique subduction along the Aleutian volcano chain to lower-angle subduction beneath the Pacific-Yakutat Plate interaction zone, forming two slab tears that enhance hot asthenosphere materials upwelling. The partial melting of the mantle wedge induced by Pacific slab dehydration and the concurrent upwelling of mantle materials jointly drove volcanism in the transition zone. Conversely, the flat subduction of the Yakutat slab into the mantle wedge overlying the Pacific slab effectively hindered the upwelling of hot hybrid materials, cooling the Pacific mantle wedge. These results offer a new perspective on the influence of slab dynamics on volcanic and magmatic processes in the region and represent an advancement in our understanding compared to previous studies, which did not resolve the tears within the slab or their geodynamic implications at this level of detail.

Melting of subducted slab dictates trace element recycling in global arcs

Arc magma acquires continental crust-like trace element signatures through selective recycling of incompatible elements from the subducted slab. The long-standing model of element recycling through aqueous fluid from altered oceanic crust (AOC) and sediment melt has been challenged by the resurgence of mélange diapir (a mix of AOC, sediment, and serpentinite) and saline aqueous fluid models. Here, we present experimental data for near-solidus sediment melts and a framework for calculating trace element concentrations in subduction fluids from metamorphosed sediment and oceanic crust. We observe that variation of element ratios in global primitive arc basalts is comparable with that of sediment and/or oceanic crustal melt, rather than (saline) aqueous fluid or mélange melt. In particular, the systematic correlation of element ratios in arc basalt corresponds to element fractionation in slab melt with temperature and therefore follows a power function. Our findings suggest that slab melt is primarily responsible for element recycling to the arc.

Geochronology, Geochemistry, and Geodynamic Relationship of the Mafic Dykes and Granites in the Qianlishan Complex, South China

The Qianlishan complex, located in Hunan Province of South China, is closely associated with intense W-dominated polymetallic mineralization. The Qianlishan complex is composed of three phases: the main-phase porphyritic and equigranular granites, granite porphyry, and mafic dykes. Geochronologically, the zircon U-Pb dating results show that the porphyritic and equigranular granites have ages of approximately 159 and 158 Ma, respectively, similar to those of mafic dykes (approximately 158 Ma), while the granite porphyry was formed later at approximately 145 Ma. Geochemically, the mafic dykes are characterized by calc-alkaline high-Mg andesite (HMA) with high MgO, TiO2, Mg#, and CA/TH index. They exhibit significantly depleted εNd(t) and εHf(t) with high Ba/La, La/Nb, and (La/Yb)N, indicating that they formed from mixing melts of depleted asthenospheric mantle and metasomatized subcontinental lithospheric mantle (SCLM). The main-phase granites are peraluminous and are characterized by high SiO2, low (La/Yb)N ratios, and relative depletion in Ba, Sr, Ti, and Eu. They also display negative correlations between La, Ce, Y, and Rb contents, suggesting that they are highly fractionated S-type granites. Furthermore, they show high εNd(t) and εHf(t), CaO/Na2O ratios, HREE, and Y contents, indicating that they were produced by parental melting of ancient basement mixed with mantle-derived components. In contrast, the granite porphyry shows A-type signature granites, with higher εNd(t) and εHf(t) and CaO/Na2O ratios than the main-phase granites but similar Zr/Nb and Zr/Hf ratios to the mafic dykes, suggesting that they are the products of partial melting of a hybrid source with ancient basement and the mafic dykes. We thus infer that the slab roll-back led to generation of Qianlishan back-arc basalt and HMA and further triggered the formation of the Qianlishan granite.

Ophiolitic Pyroxenites Record Boninite Percolation in Subduction Zone Mantle

The peridotite section of supra-subduction zone ophiolites is often crosscut by pyroxenite veins, reflecting the variety of melts that percolate through the mantle wedge, react, and eventually crystallize in the shallow lithospheric mantle. Understanding the nature of parental melts and the timing of formation of these pyroxenites provides unique constraints on melt infiltration processes that may occur in active subduction zones. This study deciphers the processes of orthopyroxenite and clinopyroxenite formation in the Josephine ophiolite (USA), using new trace and major element analyses of pyroxenite minerals, closure temperatures, elemental profiles, diffusion modeling, and equilibrium melt calculations. We show that multiple melt percolation events are required to explain the variable chemistry of peridotite-hosted pyroxenite veins, consistent with previous observations in the xenolith record. We argue that the Josephine ophiolite evolved in conditions intermediate between back-arc and sub-arc. Clinopyroxenites formed at an early stage of ophiolite formation from percolation of high-Ca boninites. Several million years later, and shortly before exhumation, orthopyroxenites formed through remelting of the Josephine harzburgites through percolation of ultra-depleted low-Ca boninites. Thus, we support the hypothesis that multiple types of boninites can be created at different stages of arc formation and that ophiolitic pyroxenites uniquely record the timing of boninite percolation in subduction zone mantle.

Evidence for ancient fractional melting, cryptic refertilization and rapid exhumation of Tethyan mantle (Civrari Ophiolite, NW Italy)

Western Tethyan peridotites exposed in the European Alps show limited amounts of partial melting and mostly fertile compositions. Here we investigate the Civrari Ophiolite (northwestern Italy), which is composed of depleted spinel-harzburgites and serpentinites associated with MOR-type gabbros and basalts. The ultramafic rocks are unique amongst western Tethyan peridotites, showing homogeneous residual compositions after ~ 15% near-fractional melting, lack of pervasive melt percolation and mineral compositions that indicate high-temperature equilibration ≥ 1200 °C. Clinopyroxene chemistry records some of the lowest abundances of Na₂O, Ce, and Zr/Hf amongst abyssal peridotites worldwide, suggesting that most abyssal peridotites have been affected by variable degrees of melt retention upon melting or cryptic melt percolation. Locally, cryptic MORB-like melt migration in Civrari peridotites produced orthopyroxene + plagioclase intergrowth around reacted clinopyroxene. These clinopyroxene preserve micron-scale chemical zoning indicating rapid cooling after melt crystallization. ¹⁴³Nd/¹⁴⁴Nd isotopic data indicate that Civrari mantle rocks, gabbros, and basalts are not in isotopic equilibrium. Civrari spinel-peridotites represent a highly radiogenic endmember amongst Western Tethys depleted spinel-peridotites, which together form a partial melting errochron of 273 Ma ± 24 Ma. Ancient near-fractional melting and cryptic melt-rock reaction cause variations in radiogenic εNd and εHf, leading to isotopic heterogeneity of Western Tethys mantle rocks. Such inherited signatures in mantle rocks are most likely to be preserved along (ultra-)slow-spreading systems and ocean-continent transition zones.

Geochemical evidence for mélange melting in global arcs

In subduction zones, sediments and hydrothermally altered oceanic crust, which together form part of the subducting slab, contribute to the chemical composition of lavas erupted at the surface to form volcanic arcs. Transport of this material from the slab to the overlying mantle wedge is thought to involve discreet melts and fluids that are released from various portions of the slab. We use a meta-analysis of geochemical data from eight globally representative arcs to show that melts and fluids from individual slab components cannot be responsible for the formation of arc lavas. Instead, the data are compatible with models that first invoke physical mixing of slab components and the mantle wedge, widely referred to as high-pressure mélange, before arc magmas are generated.

Advent of Continents: A New Hypothesis

The straightforward but unexpected relationship presented here relates crustal thickness to magma type in the Izu-Ogasawara (Bonin) and Aleutian oceanic arcs. Volcanoes along the southern segment of the Izu-Ogasawara arc and the western Aleutian arc (west of Adak) are underlain by thin crust (10–20 km). In contrast those along the northern segment of the Izu-Ogasawara arc and eastern Aleutian arc are underlain by crust ~35 km thick. Interestingly, andesite magmas dominate eruptive products from the former volcanoes and mostly basaltic lavas erupt from the latter. According to the hypothesis presented here, rising mantle diapirs stall near the base of the oceanic crust at depths controlled by the thickness of the overlying crust. Where the crust is thin, melting occurs at relatively low pressures in the mantle wedge producing andesitic magmas. Where the crust is thick, melting pressures are higher and only basaltic magmas tend to be produced. The implications of this hypothesis are: (1) the rate of continental crust accumulation, which is andesitic in composition, would have been greatest soon after subduction initiated on Earth, when most crust was thin; and (2) most andesite magmas erupted on continental crust could be recycled from “primary” andesite originally produced in oceanic arcs.

Crustal thickness control on Sr/Y signatures of recent arc magmas: an Earth scale perspective

Arc magmas originate in subduction zones as partial melts of the mantle, induced by aqueous fluids/melts liberated by the subducted slab. Subsequently, they rise through and evolve within the overriding plate crust. Aside from broadly similar features that distinguish them from magmas of other geodynamic settings (e.g., mid-ocean ridges, intraplate), arc magmas display variably high Sr/Y values. Elucidating the debated origin of high Sr/Y signatures in arc magmas, whether due to mantle-source, slab melting or intracrustal processes, is instrumental for models of crustal growth and ore genesis. Here, using a statistical treatment of >23000 whole rock geochemical data, I show that average Sr/Y values and degree of maturation (MgO depletion at peak Sr/Y values) of 19 out of 22 Pliocene-Quaternary arcs correlate positively with arc thickness. This suggests that crustal thickness exerts a first order control on the Sr/Y variability of arc magmas through the stabilization or destabilization of mineral phases that fractionate Sr (plagioclase) and Y (amphibole ± garnet). In fact, the stability of these mineral phases is function of the pressure at which magma evolves, which depends on crustal thickness. The data presented show also that high Sr/Y Pliocene-Quaternary intermediate-felsic arc rocks have a distinct origin from their Archean counterparts.

Implications of estimated magmatic additions and recycling losses at the subduction zones of accretionary (non-collisional) and collisional (suturing) orogens

Arc magmatism at subduction zones (SZs) most voluminously supplies juvenile igneous material to build rafts of continental and intra-oceanic or island arc (CIA) crust. Return or recycling of accumulated CIA material to the mantle is also most vigorous at SZs. Recycling is effected by the processes of sediment subduction, subduction erosion, and detachment and sinking of deeply underthrust sectors of CIA crust. Long-term (>10–20 Ma) rates of additions and losses can be estimated from observational data gathered where oceanic crust underruns modern, long-running (Cenozoic to mid-Mesozoic) ocean-margin subduction zones (OMSZs, e.g. Aleutian and South America SZs). Long-term rates can also be observationally assessed at Mesozoic and older crust-suturing subduction zone (CSSZs) where thick bodies of CIA crust collided in tectonic contact (e.g. Wopmay and Appalachian orogens, India and SE Asia). At modern OMSZs arc magmatic additions at intra-oceanic arcs and at continental margins are globally estimated atc. 1.5 AU andc. 1.0 AU, respectively (1 AU, or Armstrong Unit,=1 km3a−1of solid material). During collisional suturing at fossil CSSZs, global arc magmatic addition is estimated at 0.2 AU. This assessment presumes that in the past the global length of crustal collision zones averagedc. 6000 km, which is one-half that under way since the early Tertiary. The average long-term rate of arc magmatic additions extracted from modern OMSZs and older CSSZs is thus evaluated at 2.7 AU. Crustal recycling at Mesozoic and younger OMSZs is assessed atc. 60 km3Ma−1km−1(c. 60% by subduction erosion). The corresponding global recycling rate isc. 2.5 AU. At CSSZs of Mesozoic, Palaeozoic and Proterozoic age, the combined upper and lower plate losses of CIA crust via subduction erosion, sediment subduction, and lower plate crustal detachment and sinking are assessed far less securely atc. 115 km3Ma−1km−1. At a global length of 6000 km, recycling at CSSZs is accordinglyc. 0.7 AU. The collective loss of CIA crust estimated for modern OMSZs and for older CSSZs is thus estimated atc. 3.2 AU. SZ additions (+2.7 AU) and subtractions (−3.2 AU) are similar. Because many uncertainties and assumptions are involved in assessing and applying them to the deep past, the net growth of CIA crust during at least Phanerozoic time is viewed as effectively nil. With increasing uncertainty, the long-term balance can be applied to the Proterozoic, but not before the initiation of the present style of subduction atc. 3 Ga. Allowing that since this time a rounded-down rate of recycling of 3 AU is applicable, a startlingly high volume of CIA crust equal to that existing now has been recycled to the mantle. Although the recycled volume (c. 9×109km3) is small (c. 1%) compared with that of the mantle, it is large enough to impart to the mantle the signature of recycled CIA crust. Because subduction zones are not spatially fixed, and their average global lengths have episodically been less or greater than at present, recycling must have contributed significantly to creating recognized heterogeneities in mantle geochemistry.