Numerical modeling the process of deep slab dehydration and magmatism

This study uses a 2D high-resolution thermo-mechanical coupled model to investigate the dynamic processes of deep plate hydration, dehydration, and subsequent magmatic activity in ocean-continent subduction zones. We reveal the pathways and temporal evolution of water transport to the deep mantle during the subduction process. Plate dehydration plays a critical role in triggering partial melting of the deep mantle and related magmatic activity. Our study shows significant differences in the volumes of melt produced at different depths, with dehydration reactions in deeper regions being weaker compared to shallower ones. It takes a longer time to reach the suitable P-T conditions for hydrous melting in the deep mantle. The results highlight the geophysical significance of water transport in deep subduction zones and its role in magmatic processes, particularly in the formation of magma chambers beneath continental plates.

An alternative to the igneous crust fluid + sediment melt paradigm for arc lava geochemistry

A long-standing paradigm of arc geochemistry is that the trace element compositions of arc lavas arise from two compositionally distinct slab components: an aqueous dehydration fluid from the subducting igneous ocean crust that transports "fluid-mobile" elements, such as barium (Ba), and a sediment melt that supplies thorium (Th) and the light rare earth elements. This two-component framework has been widely called upon to explain global geochemical trends as well as geochemical variations within individual arcs, such as the Marianas. Here, we show that this paradigm is inconsistent with mass balance, due to the low Ba contents of igneous ocean crust, and with experimental data, which show that aqueous fluids from the igneous oceanic crust would be too dilute to substantially affect arc compositions. Observations previously attributed to the sediment melt/igneous-crust-fluid hypothesis are better explained by diverse subducting sediment compositions coupled with ambient mantle wedge heterogeneity, both globally and for the Marianas.

A HIMU-like component in Mariana Convergent Margin magma sources during initial arc rifting revealed by melt inclusions

Compositions of island arc and back-arc basin basalts are often used to trace the recycling of subducted materials. However, the contribution of subducted components to the mantle source during initial arc rifting before back-arc basin spreading is not yet well constrained. The northernmost Mariana arc is ideal for studying this because the transition from rifting to back-arc spreading is happening here. Here we report major and trace element and Pb isotopic compositions of olivine-hosted melt inclusions from lavas erupted during initial rifting at 24°N (NSP-24) and compare them with those in active arc front at 21°N and mature back-arc basin at 18°N. NSP-24 high-K melt inclusions have highly radiogenic Pb compositions and are close to those of the HIMU end-member, suggesting the presence of this component in the magma source. The HIMU-like component may be stored in the over-riding plate and released into arc magma with rifting. HIMU-type seamounts may be subducted elsewhere beneath the Mariana arc, but obvious HIMU-type components appear only in the initial stages of arc rifting due to the low melting degree and being consumed during the process of back-arc spreading.

Ghost-arc geochemical anomaly at a spreading ridge caused by supersized flat subduction

The Southern Atlantic-Southwest Indian ridges (SASWIR) host mid-ocean ridge basalts with a residual subduction-related geochemical fingerprint (i.e., a ghost-arc signature) of unclear origin. Here, we show through an analysis of plate kinematic reconstructions and seismic tomography models that the SASWIR subduction-modified mantle source formed in the Jurassic close to the Georgia Islands slab (GI) and remained near-stationary in the mantle reference frame. In this analysis, the GI lies far inboard the Jurassic Patagonian-Antarctic Peninsula active margin. This was formerly attributed to a large-scale flat subduction event in the Late Triassic-Early Jurassic. We propose that during this flat slab stage, the subduction-modified mantle areas beneath the Mesozoic active margin and surrounding sutures zones may have been bulldozed inland by >2280 km. After the demise of the flat slab, this mantle anomaly remained near-stationary and was sampled by the Karoo mantle plume 183 Million years (Myr) ago and again since 55 Myr ago by the SASWIR. We refer to this process as asthenospheric anomaly telescoping. This study provides a hitherto unrecognized geodynamic effect of flat subduction, the viability of which we support through numerical modeling.

Sources of dehydration fluids underneath the Kamchatka arc

Fluids mediate the transport of subducted slab material and play a crucial role in the generation of arc magmas. However, the source of subduction-derived fluids remains debated. The Kamchatka arc is an ideal subduction zone to identify the source of fluids because the arc magmas are comparably mafic, their source appears to be essentially free of subducted sediment-derived components, and subducted Hawaii-Emperor Seamount Chain (HESC) is thought to contribute a substantial fluid flux to the Kamchatka magmas. Here we show that Tl isotope ratios are unique tracers of HESC contribution to Kamchatka arc magma sources. In conjunction with trace element ratios and literature data, we trace the progressive dehydration and melting of subducted HESC across the Kamchatka arc. In succession, serpentine (<100 km depth), lawsonite (100–250 km depth) and phengite (>250 km depth) break down and produce fluids that contribute to arc magmatism at the Eastern Volcanic Front (EVF), Central Kamchatka Depression (CKD), and Sredinny Ridge (SR), respectively. However, given the Tl-poor nature of serpentine and lawsonite fluids, simultaneous melting of subducted HESC is required to explain the HESC-like Tl isotope signatures observed in EVF and CKD lavas. In the absence of eclogitic crust melting processes in this region of the Kamchatka arc, we propose that progressive dehydration and melting of a HESC-dominated mélange offers the most compelling interpretation of the combined isotope and trace element data.

Fluids released from progressive breakdown of minerals at increasing pressure within a mélange may explain the trace element systematics and stable thallium isotope data of the Kamchatka arc lavas from volcanic front to back arc.

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.

Lithium, Oxygen and Magnesium Isotope Systematics of Volcanic Rocks in the Okinawa Trough: Implications for Plate Subduction Studies

Determining the influence of subduction input on back-arc basin magmatism is important for understanding material transfer and circulation in subduction zones. Although the mantle source of Okinawa Trough (OT) magmas is widely accepted to be modified by subducted components, the role of slab-derived fluids is poorly defined. Here, major element, trace element, and Li, O and Mg isotopic compositions of volcanic lavas from the middle OT (MOT) and southern OT (SOT) were analyzed. Compared with the MOT volcanic lavas, the T9-1 basaltic andesite from the SOT exhibited positive Pb anomalies, significantly lower Nd/Pb and Ce/Pb ratios, and higher Ba/La ratios, indicating that subducted sedimentary components affected SOT magma compositions. The δ7Li, δ18O, and δ26Mg values of the SOT basaltic andesite (−5.05‰ to 4.98‰, 4.83‰ to 5.80‰ and −0.16‰ to −0.09‰, respectively) differed from those of MOT volcanic lavas. Hence, the effect of the Philippine Sea Plate subduction component, (low δ7Li and δ18O and high δ26Mg) on magmas in the SOT was clearer than that in the MOT. This contrast likely appears because the amounts of fluids and/or melts derived from altered oceanic crust (AOC, lower δ18O) and/or subducted sediment (lower δ7Li, higher δ18O and δ26Mg) injected into magmas in the SOT are larger than those in the MOT and because the injection ratio between subducted AOC and sediment is always >1 in the OT. The distance between the subducting slab and overlying magma may play a significant role in controlling the differences in subduction components injected into magmas between the MOT and SOT.

Molybdenum isotopes unmask slab dehydration and melting beneath the Mariana arc

How serpentinites in the forearc mantle and subducted lithosphere become involved in enriching the subarc mantle source of arc magmas is controversial. Here we report molybdenum isotopes for primitive submarine lavas and serpentinites from active volcanoes and serpentinite mud volcanoes in the Mariana arc. These data, in combination with radiogenic isotopes and elemental ratios, allow development of a model whereby shallow, partially serpentinized and subducted forearc mantle transfers fluid and melt from the subducted slab into the subarc mantle. These entrained forearc mantle fragments are further metasomatized by slab fluids/melts derived from the dehydration of serpentinites in the subducted lithospheric slab. Multistage breakdown of serpentinites in the subduction channel ultimately releases fluids/melts that trigger Mariana volcanic front volcanism. Serpentinites dragged down from the forearc mantle are likely exhausted at >200 km depth, after which slab-derived serpentinites are responsible for generating slab melts.

How the subducted oceanic lithosphere provides fluids and melts to flux the subarc mantle source of arc magmas is controversial. Here the authors use Mo and other isotopes to show serpentinites formed in both the forearc mantle and the subducted lithosphere contribute to generating arc magmas.

Building cratonic keels in Precambrian plate tectonics

The ancient cores of continents (cratons) are underlain by mantle keels-volumes of melt-depleted, mechanically resistant, buoyant and diamondiferous mantle up to 350 kilometres thick, which have remained isolated from the hotter and denser convecting mantle for more than two billion years. Mantle keels formed only in the Early Earth (approximately 1.5 to 3.5 billion years ago in the Precambrian eon); they have no modern analogues^1-4. Many keels show layering in terms of degree of melt depletion^5-7. The origin of such layered lithosphere remains unknown and may be indicative of a global tectonics mode (plate rather than plume tectonics) operating in the Early Earth. Here we investigate the possible origin of mantle keels using models of oceanic subduction followed by arc-continent collision at increased mantle temperatures (150-250 degrees Celsius higher than the present-day values). We demonstrate that after Archaean plate tectonics began, the hot, ductile, positively buoyant, melt-depleted sublithospheric mantle layer located under subducting oceanic plates was unable to subduct together with the slab. The moving slab left behind craton-scale emplacements of viscous protokeel beneath adjacent continental domains. Estimates of the thickness of this sublithospheric depleted mantle show that this mechanism was efficient at the time of the major statistical maxima of cratonic lithosphere ages. Subsequent conductive cooling of these protokeels would produce mantle keels with their low modern temperatures, which are suitable for diamond formation. Precambrian subduction of oceanic plates with highly depleted mantle is thus a prerequisite for the formation of thick layered lithosphere under the continents, which permitted their longevity and survival in subsequent plate tectonic processes.

Reply to 'Evidence for simple volcanic rifting not complex subduction initiation in the Laxmi Basin'

Recently, Pandey et al. proposed relict subduction initiation occurred along a passive margin in the northwest Indian Ocean, however, Clift et al. questioned their evidence for subduction initiation, suggesting that simpler rifting-related processes could more simply explain the available data. Here, Pandey et al. reply to Clift et al.’s comment, and argue that geochemical and isotope data for Laxmi basin lavas distinctly imply relict subduction initiation.

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.

Ancient recycled lower crust in the mantle source of recent Italian magmatism

Recycling of Earth’s crust through subduction and delamination contributes to mantle heterogeneity. Melt inclusions in early crystallised magmatic minerals record greater geochemical variability than host lavas and more fully reflect the heterogeneity of magma sources. To date, use of multiple isotope systems on small (< 300 μm) melt inclusions was hampered by analytical limitations. Here we report the first coupled Sr-Nd-Pb isotope data on individual melt inclusions from potassium-rich lavas from neighbouring Quaternary volcanoes in central Italy and infer the presence of a previously unidentified ancient lower crustal component in the mantle. We suggest derivation from Variscan or older basement included in the upper mantle by either delamination, sediment recycling, subduction erosion and/or slab detachment processes during Cenozoic subduction and collision of the western Mediterranean. The capability to determine isotope ratios in individual melt inclusions permits the detection of distinctive mantle contaminants and can provide insights into how geodynamic processes affect subduction recycling.

Recycling of Earth’s crust through subduction and delamination contributes to mantle heterogeneity. Here, the authors measure coupled Sr-Nd-Pb isotope compositions of melt inclusions in Italian potassium-rich lavas, they suggest their results indicate a potential ancient lower crustal component in the mantle source.

Heavy oxygen recycled into the lithospheric mantle

Magmas in volcanic arcs have geochemical and isotopic signatures that can be related to mantle metasomatism due to fluids and melts released by the down-going oceanic crust and overlying sediments, which modify the chemistry and mineralogy of the mantle wedge. However, the effectiveness of subduction-related metasomatic processes is difficult to evaluate because the composition of arc magmas is often overprinted by interactions with crustal lithologies occurring during magma ascent and emplacement. Here, we show unequivocal evidence for recycling of continental crust components into the mantle. Veined peridotite xenoliths sampled from Tallante monogenetic volcanoes in the Betic Cordillera (southern Spain) provide insights for mantle domains that reacted with Si-rich melts derived by partial melting of subducted crustal material. Felsic veins crosscutting peridotite and the surrounding orthopyroxene-rich metasomatic aureoles show the highest ¹⁸O/¹⁶O ratios measured to date in upper mantle assemblages worldwide. The anomalously high oxygen isotope compositions, coupled with very high ⁸⁷Sr/⁸⁶Sr values, imply the continental crust origin of the injected melts. Isotopic anomalies are progressively attenuated in peridotite away from the veins, showing ¹⁸O isotope variations well correlated with the amount of newly formed orthopyroxene. Diffusion may also affect the isotope ratios of mantle rocks undergoing crustal metasomatism due to the relaxation of ¹⁸O isotope anomalies to normal mantle values through time. Overall, the data define an O isotope “benchmark” allowing discrimination between mantle sources that attained re-equilibration after metasomatism (>5 Myr) and those affected by more recent subduction-derived enrichment processes.

The causes of spatiotemporal variations in erupted fluxes and compositions along a volcanic arc

Decades of study on volcanic arcs have provided insight into the overarching processes that control magmatism, and how these processes manifest at individual volcanoes. However, the causes of ubiquitous and dramatic intra-arc variations in volcanic flux and composition remain largely unresolved. Investigating such arc-scale issues requires greater quantitative comparison of geophysical and geochemical data, linked through sets of common intensive variables. To work towards these goals, we use observed lava compositions to estimate the heat budget associated with Quaternary volcanism in the Cascades Arc and compare this to the heat required to produce the observed geophysical properties of the crust. Here we show that along-strike volcanic variability in the Quaternary Cascades Arc is primarily related to variations in the flux of basalt into the crust, rather than variations in their crustal storage history. This approach shows promise for studying other large-scale frontier geologic problems in volcanic arcs.

The primary causes of dramatic variations in volcanic flux and composition along strike in subduction zones remain largely unknown. Here we use a promising new approach to show that along-strike volcanic variability in the Quaternary Cascades Arc is primarily due to variations in the flux of basalt into the base of the crust, rather than crustal magma storage.

Episodic zircon age spectra mimic fluctuations in subduction

Decades of geochronological work have shown the temporal distribution of zircon ages to be episodic on billion-year timescales and seemingly coincident with the lifecycle of supercontinents, but the physical processes behind this episodicity remain contentious. The dominant, end-member models of fluctuating magmatic productivity versus selective preservation of zircon during times of continental assembly have important and very different implications for long-term, global-scale phenomena, including the history of crustal growth, the initiation and evolution of plate tectonics, and the tempo of mantle outgassing over billions of years. Consideration of this episodicity has largely focused on the Precambrian, but here we analyze a large collection of Phanerozoic zircon ages in the context of global, full-plate tectonic models that extend back to the mid-Paleozoic. We scrutinize two long-lived and relatively simple active margins, and show that along both, a relationship between the regional subduction flux and zircon age distribution is evident. In both cases, zircon age peaks correspond to intervals of high subduction flux with a ~10–30 Ma time lag (zircons trailing subduction), illuminating a possibly intrinsic delay in the subduction-related magmatic system. We also show that subduction fluxes provide a stronger correlation to zircon age distributions than subduction lengths do, implying that convergence rates play a significant role in regulating the volume of melting in subduction-related magmatic systems, and thus crustal growth.

Arc-like magmas generated by mélange-peridotite interaction in the mantle wedge

The mechanisms of transfer of crustal material from the subducting slab to the overlying mantle wedge are still debated. Mélange rocks, formed by mixing of sediments, oceanic crust, and ultramafics along the slab-mantle interface, are predicted to ascend as diapirs from the slab-top and transfer their compositional signatures to the source region of arc magmas. However, the compositions of melts that result from the interaction of mélanges with a peridotite wedge remain unknown. Here we present experimental evidence that melting of peridotite hybridized by mélanges produces melts that carry the major and trace element abundances observed in natural arc magmas. We propose that differences in nature and relative contributions of mélanges hybridizing the mantle produce a range of primary arc magmas, from tholeiitic to calc-alkaline. Thus, assimilation of mélanges into the wedge may play a key role in transferring subduction signatures from the slab to the source of arc magmas.

Mélange rocks are predicted to form at the slab-mantle interface in most subduction zones, but their role in arc magmatism is still debated. Here, the authors show that melting of peridotite hybridized by mélange rocks produces melts that carry the major and trace element abundances of natural arc magmas.

Source and fractionation controls on subduction-related plutons and dike swarms in southern Patagonia (Torres del Paine area) and the low Nb/Ta of upper crustal igneous rocks

The subduction system in southern Patagonia provides direct evidence for the variability of the position of an active continental arc with respect to the subducting plate through time, but the consequences on the arc magmatic record are less well studied. Here we present a geochemical and geochronological study on small plutons and dykes from the upper crust of the southern Patagonian Andes at ~ 51°S, which formed as a result of the subduction of the Nazca and Antarctic plates beneath the South American continent. In situ U-Pb geochronology on zircons and bulk rock geochemical data of plutonic and dyke rocks are used to constrain the magmatic evolution of the retro-arc over the last 30 Ma. We demonstrate that these combined U-Pb and geochemical data for magmatic rocks track the temporal and spatial migration of the active arc, and associated retro-arc magmatism. Our dataset indicates that the rear-arc area is characterized by small volumes of alkaline basaltic magmas at 29-30 Ma that are characterized by low La/Nb and Th/Nb ratios with negligible arc signatures. Subsequent progressive eastward migration of the active arc culminated with the emplacement of calc-alkaline plutons and dikes ~ 17-16 Ma with elevated La/Nb and Th/Nb ratios and typical subduction signatures constraining the easternmost position of the southern Patagonian batholith at that time. Geochemical data on the post-16 Ma igneous rocks including the Torres del Paine laccolith indicate an evolution to transitional K-rich calc-alkaline magmatism at 12.5 ± 0.2 Ma. We show that trace element ratios such as Nb/Ta and Dy/Yb systematically decrease with increasing SiO₂, for both the 17-16 Ma calc-alkaline and the 12-13 Ma K-rich transitional magmatism. In contrast, Th/Nb and La/Nb monitor the changes in the source composition of these magmas. We suggest that the transition from the common calc-alkaline to K-rich transitional magmatism involves a change in the source component, while the trace element ratios, such as Nb/Ta and Dy/Yb, of derivative higher silica content liquids are controlled by similar fractionating mineral assemblages. Analysis of a global compilation of Nb/Ta ratios of arc magmatic rocks and simple geochemical models indicate that amphibole and variable amounts of biotite exert a major control on the low Dy/Yb and Nb/Ta of derivative granitic liquids. Lastly, we suggest that the low Nb/Ta ratio of silica-rich magmas is a natural consequence of biotite fractionation and that alternative models such as amphibolite melting in subduction zones and diffusive fractionation are not required to explain the Nb/Ta ratio of the upper continental crust.

Partial Melting of Subducted Sediments Produced Early Mesozoic Calc-alkaline Lamprophyres from Northern Guangxi Province, South China

There is growing agreement that subducted sediments recycled into the deep mantle could make a significant contribution to the generation of various mantle-derived rocks. However, solid evidence and examples to support this view are few, and whether or not the subducted sediments can act as the dominating material source for the magma is unclear. Here, we report a comprehensive geochemical study that demonstrates that the newly identified Early Mesozoic calc-alkaline lamprophyres in the northern Guangxi Province, southeastern Yangtze Block in South China were likely derived in large part from the partial melting of the subducted terrigenous sediments in the deep mantle. The investigated lamprophyres are SiO₂-rich minettes, characterized by moderate TFeO and MgO and high Mg^# (>70). The multi-element pattern shows a typical crustal-like signature, such as enrichments in large-ion lithophile elements (LILE) and light rare earth elements (LREE) with troughs in Nb-Ta, Ti and Eu and peaks in Th-U and Pb. These rocks also show sediment-like ratios of Nb/U, Nb/Th and Ce/Pb, together with extremely radiogenic ⁸⁷Sr/⁸⁶Sr (0.71499–0.71919), unradiogenic ¹⁴³Nd/¹⁴⁴Nd (0.51188–0.51195) and radiogenic ²⁰⁷Pb/²⁰⁴Pb (15.701–15.718) isotopic compositions.

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.

Continental versus oceanic subduction zones

Subduction zones are tectonic expressions of convergent plate margins, where crustal rocks descend into and interact with the overlying mantle wedge. They are the geodynamic system that produces mafic arc volcanics above oceanic subduction zones but high- to ultrahigh-pressure metamorphic rocks in continental subduction zones. While the metamorphic rocks provide petrological records of orogenic processes when descending crustal rocks undergo dehydration and anataxis at forearc to subarc depths beneath the mantle wedge, the arc volcanics provide geochemical records of the mass transfer from the subducting slab to the mantle wedge in this period though the mantle wedge becomes partially melted at a later time. Whereas the mantle wedge overlying the subducting oceanic slab is of asthenospheric origin, that overlying the descending continental slab is of lithospheric origin, being ancient beneath cratons but juvenile beneath marginal arcs. In either case, the mantle wedge base is cooled down during the slab–wedge coupled subduction. Metamorphic dehydration is prominent during subduction of crustal rocks, giving rise to aqueous solutions that are enriched in fluid-mobile incompatible elements. Once the subducting slab is decoupled from the mantle wedge, the slab–mantle interface is heated by lateral incursion of the asthenospheric mantle to allow dehydration melting of rocks in the descending slab surface and the metasomatized mantle wedge base, respectively. Therefore, the tectonic regime of subduction zones changes in both time and space with respect to their structures, inputs, processes and products. Ophiolites record the tectonic conversion from seafloor spreading to oceanic subduction beneath continental margin, whereas ultrahigh-temperature metamorphic events mark the tectonic conversion from compression to extension in orogens.

Magnesium isotope geochemistry in arc volcanism

Incorporation of subducted slab in arc volcanism plays an important role in producing the geochemical and isotopic variations in arc lavas. The mechanism and process by which the slab materials are incorporated, however, are still uncertain. Here, we report, to our knowledge, the first set of Mg isotopic data for a suite of arc lava samples from Martinique Island in the Lesser Antilles arc, which displays one of the most extreme geochemical and isotopic ranges, although the origin of this variability is still highly debated. We find the δ(26)Mg of the Martinique Island lavas varies from -0.25 to -0.10, in contrast to the narrow range that characterizes the mantle (-0.25 ± 0.04, 2 SD). These high δ(26)Mg values suggest the incorporation of isotopically heavy Mg from the subducted slab. The large contrast in MgO content between peridotite, basalt, and sediment makes direct mixing between sediment and peridotite, or assimilation by arc crust sediment, unlikely to be the main mechanism to modify Mg isotopes. Instead, the heavy Mg isotopic signature of the Martinique arc lavas requires that the overall composition of the mantle wedge is buffered and modified by the preferential addition of heavy Mg isotopes from fluids released from the altered subducted slab during fluid-mantle interaction. This, in turn, suggests transfer of a large amount of fluid-mobile elements from the subducting slab to the mantle wedge and makes Mg isotopes an excellent tracer of deep fluid migration.

The zircon archive of continent formation through time

The strong resilience of the mineral zircon and its ability to host a wealth of isotopic information make it the best deep-time archive of Earth's continental crust. Zircon is found in most felsic igneous rocks, can be precisely dated and can fingerprint magmatic sources; thus, it has been widely used to document the formation and evolution of continental crust, from pluton- to global-scale. Here, we present a review of major contributions that zircon studies have made in terms of understanding key questions involving the formation of the continents. These include the conditions of continent formation on early Earth, the onset of plate tectonics and subduction, the rate of crustal growth through time and the governing balance of continental addition v. continental loss, and the role of preservation bias in the zircon record.

Supplementary material:

A compilation used in this study of previously published detrital zircon U-Pb-Hf isotope data are available at http://www.geolsoc.org.uk/SUP18791

Subduction of fracture zones controls mantle melting and geochemical signature above slabs

For some volcanic arcs, the geochemistry of volcanic rocks erupting above subducted oceanic fracture zones is consistent with higher than normal fluid inputs to arc magma sources. Here we use enrichment of boron (B/Zr) in volcanic arc lavas as a proxy to evaluate relative along-strike inputs of slab-derived fluids in the Aleutian, Andean, Cascades and Trans-Mexican arcs. Significant B/Zr spikes coincide with subduction of prominent fracture zones in the relatively cool Aleutian and Andean subduction zones where fracture zone subduction locally enhances fluid introduction beneath volcanic arcs. Geodynamic models of subduction have not previously considered how fracture zones may influence the melt and fluid distribution above slabs. Using high-resolution three-dimensional coupled petrological-thermomechanical numerical simulations of subduction, we show that enhanced production of slab-derived fluids and mantle wedge melts concentrate in areas where fracture zones are subducted, resulting in significant along-arc variability in magma source compositions and processes.

Subduction of the oceanic Hikurangi Plateau and its impact on the Kermadec arc

Large igneous province subduction is a rare process on Earth. A modern example is the subduction of the oceanic Hikurangi Plateau beneath the southern Kermadec arc, offshore New Zealand. This segment of the arc has the largest total lava volume erupted and the highest volcano density of the entire Kermadec arc. Here we show that Kermadec arc lavas south of ~32°S have elevated Pb and Sr and low Nd isotope ratios, which argues, together with increasing seafloor depth, forearc retreat and crustal thinning, for initial Hikurangi Plateau-Kermadec arc collision ~250 km north of its present position. The combined data set indicates that a much larger portion of the Hikurangi Plateau (the missing Ontong Java Nui piece) than previously believed has already been subducted. Oblique plate convergence caused southward migration of the thickened and buoyant oceanic plateau crust, creating a buoyant 'Hikurangi' mélange beneath the Moho that interacts with ascending arc melts.

A 'hidden' 18O-enriched reservoir in the sub-arc mantle

Plate subduction continuously transports crustal materials with high-δ(18)O values down to the mantle wedge, where mantle peridotites are expected to achieve the high-δ(18)O features. Elevated δ(18)O values relative to the upper mantle value have been reported for magmas from some subduction zones. However, peridotites with δ(18)O values significantly higher than the well-defined upper mantle values have never been observed from modern subduction zones. Here we present in-situ oxygen isotope data of olivine crystals in Sailipu mantle xenoliths from South Tibet, which have been subjected to a long history of Tethyan subduction before the India-Asia collision. Our data identify for the first time a metasomatized mantle that, interpreted as the sub-arc lithospheric mantle, shows anomalously enriched oxygen isotopes (δ(18)O = +8.03 ± 0.28 ‰). Such a high-δ(18)O mantle commonly does not contribute significantly to typical island arc basalts. However, partial melting or contamination of such a high-δ(18)O mantle is feasible to account for the high-δ(18)O signatures in arc basalts.

Separation of supercritical slab-fluids to form aqueous fluid and melt components in subduction zone magmatism

Subduction-zone magmatism is triggered by the addition of H(2)O-rich slab-derived components: aqueous fluid, hydrous partial melts, or supercritical fluids from the subducting slab. Geochemical analyses of island arc basalts suggest two slab-derived signatures of a melt and a fluid. These two liquids unite to a supercritical fluid under pressure and temperature conditions beyond a critical endpoint. We ascertain critical endpoints between aqueous fluids and sediment or high-Mg andesite (HMA) melts located, respectively, at 83-km and 92-km depths by using an in situ observation technique. These depths are within the mantle wedge underlying volcanic fronts, which are formed 90 to 200 km above subducting slabs. These data suggest that sediment-derived supercritical fluids, which are fed to the mantle wedge from the subducting slab, react with mantle peridotite to form HMA supercritical fluids. Such HMA supercritical fluids separate into aqueous fluids and HMA melts at 92 km depth during ascent. The aqueous fluids are fluxed into the asthenospheric mantle to form arc basalts, which are locally associated with HMAs in hot subduction zones. The separated HMA melts retain their composition in limited equilibrium with the surrounding mantle. Alternatively, they equilibrate with the surrounding mantle and change the major element chemistry to basaltic composition. However, trace element signatures of sediment-derived supercritical fluids remain more in the melt-derived magma than in the fluid-induced magma, which inherits only fluid-mobile elements from the sediment-derived supercritical fluids. Separation of slab-derived supercritical fluids into melts and aqueous fluids can elucidate the two slab-derived components observed in subduction zone magma chemistry.

Slab melting versus slab dehydration in subduction-zone magmatism

The second critical endpoint in the basalt-H(2)O system was directly determined by a high-pressure and high-temperature X-ray radiography technique. We found that the second critical endpoint occurs at around 3.4 GPa and 770 °C (corresponding to a depth of approximately 100 km in a subducting slab), which is much shallower than the previously estimated conditions. Our results indicate that the melting temperature of the subducting oceanic crust can no longer be defined beyond this critical condition and that the fluid released from subducting oceanic crust at depths greater than 100 km under volcanic arcs is supercritical fluid rather than aqueous fluid and/or hydrous melts. The position of the second critical endpoint explains why there is a limitation to the slab depth at which adakitic magmas are produced, as well as the origin of across-arc geochemical variations of trace elements in volcanic rocks in subduction zones.

Trace element signature of subduction-zone fluids, melts and supercritical liquids at 120–180 km depth

Fluids and melts liberated from subducting oceanic crust recycle lithophile elements back into the mantle wedge, facilitate melting and ultimately lead to prolific subduction-zone arc volcanism. The nature and composition of the mobile phases generated in the subducting slab at high pressures have, however, remained largely unknown. Here we report direct LA-ICPMS measurements of the composition of fluids and melts equilibrated with a basaltic eclogite at pressures equivalent to depths in the Earth of 120-180 km and temperatures of 700-1,200 degrees C. The resultant liquid/mineral partition coefficients constrain the recycling rates of key elements. The dichotomy of dehydration versus melting at 120 km depth is expressed through contrasting behaviour of many trace elements (U/Th, Sr, Ba, Be and the light rare-earth elements). At pressures equivalent to 180 km depth, however, a supercritical liquid with melt-like solubilities for the investigated trace elements is observed, even at low temperatures. This mobilizes most of the key trace elements (except the heavy rare-earth elements, Y and Sc) and thus limits fluid-phase transfer of geochemical signatures in subduction zones to pressures less than 6 GPa.