Mantle melting as a function of water content beneath back-arc basins

Subduction zone magmas are characterized by high concentrations of H_(2)O, presumably derived from the subducted plate and ultimately responsible for melting at this tectonic setting. Previous studies of the role of water during mantle melting beneath back-arc basins found positive correlations between the H_(2)O concentration of the mantle (H_(2)O_o ) and the extent of melting (F), in contrast to the negative correlations observed at mid-ocean ridges. Here we examine data compiled from six back-arc basins and three mid-ocean ridge regions. We use TiO_2 as a proxy for F, then use F to calculate H_(2)O_o from measured H_(2)O concentrations of submarine basalts. Back-arc basins record up to 0.5 wt % H_(2)O or more in their mantle sources and define positive, approximately linear correlations between H_(2)O_o and F that vary regionally in slope and intercept. Ridge-like mantle potential temperatures at back-arc basins, constrained from Na-Fe systematics (1350°–1500°C), correlate with variations in axial depth and wet melt productivity (∼30–80% F/wt % H_(2)O_o ). Water concentrations in back-arc mantle sources increase toward the trench, and back-arc spreading segments with the highest mean H_(2)O_o are at anomalously shallow water depths, consistent with increases in crustal thickness and total melt production resulting from high H_(2)O. These results contrast with those from ridges, which record low H_(2)O_o (<0.05 wt %) and broadly negative correlations between H_(2)O_o and F that result from purely passive melting and efficient melt focusing, where water and melt distribution are governed by the solid flow field. Back-arc basin spreading combines ridge-like adiabatic melting with nonadiabatic mantle melting paths that may be independent of the solid flow field and derive from the H_(2)O supply from the subducting plate. These factors combine significant quantitative and qualitative differences in the integrated influence of water on melting phenomena in back-arc basin and mid-ocean ridge settings

Subducting carbon

A hidden carbon cycle exists inside Earth. Every year, megatons of carbon disappear into subduction zones, affecting atmospheric carbon dioxide and oxygen over Earth’s history. Here we discuss the processes that move carbon towards subduction zones and transform it into fluids, magmas, volcanic gases and diamonds. The carbon dioxide emitted from arc volcanoes is largely recycled from subducted microfossils, organic remains and carbonate precipitates. The type of carbon input and the efficiency with which carbon is remobilized in the subduction zone vary greatly around the globe, with every convergent margin providing a natural laboratory for tracing subducting carbon

NanoSIMS results from olivine-hosted melt embayments: Magma ascent rate during explosive basaltic eruptions

The explosivity of volcanic eruptions is governed in part by the rate at which magma ascends and degasses. Because the time scales of eruptive processes can be exceptionally fast relative to standard geochronometers, magma ascent rate remains difficult to quantify. Here we use as a chronometer concentration gradients of volatile species along open melt embayments within olivine crystals. Continuous degassing of the external melt during magma ascent results in diffusion of volatile species from embayment interiors to the bubble located at their outlets. The novel aspect of this study is the measurement of concentration gradients in five volatile elements (CO2, H2O, S, Cl, F) at fine-scale (5–10 μm) using the NanoSIMS. The wide range in diffusivity and solubility of these different volatiles provides multiple constraints on ascent timescales over a range of depths. We focus on four 100–200 μm, olivine-hosted embayments erupted on October 17, 1974 during the sub-Plinian eruption of Volcán de Fuego. H2O, CO2, and S all decrease toward the embayment outlet bubble, while F and Cl increase or remain roughly constant. Compared to an extensive melt inclusion suite from the same day of the eruption, the embayments have lost both H2O and CO2 throughout the entire length of the embayment. We fit the profiles with a 1-D numerical diffusion model that allows varying diffusivities and external melt concentrations as a function of pressure. Assuming a constant decompression rate from the magma storage region at approximately 220 MPa to the surface, H2O, CO2 and S profiles for all embayments can be fit with a relatively narrow range in decompression rates of 0.3–0.5 MPa/s, equivalent to 11–17 m/s ascent velocity and an 8 to 12 minute duration of magma ascent from ~ 10 km depth. A two stage decompression model takes advantage of the different depth ranges over which CO2 and H2O degas, and produces good fits given an initial stage of slow decompression (0.05–0.3 MPa/s) at high pressure (< 145 MPa), with similar decompression rates to the single-stage model for the shallower stage. The magma ascent rates reported here are among the first for explosive basaltic eruptions and demonstrate the potential of the embayment method for quantifying magmatic timescales associated with eruptions of different vigor

Thermochemical evolution of the sub-arc mantle due to back-arc spreading

We present the results of a series of numerical geodynamic experiments designed to characterize the thermal and compositional evolution of the sub-arc mantle in response to spreading in the back-arc. We find large changes in both the temperature and composition of the sub-arc mantle with time as the BASC migrates away from the arc. In particular, the sub-arc mantle becomes increasingly more depleted with time following the onset of spreading, as mantle that has experienced decompression melting and melt extraction beneath the BASC is gradually drawn beneath the arc plate by slab-induced corner flow. The rate at which this depletion increases during the ~2 Myr immediately following the onset of spreading is controlled by the spreading rate at the BASC, with faster spreading leading to a more rapid increase in depletion. Following this initial period, depletion within the sub-arc mantle continues to increase at a somewhat slower pace. During this phase, the rate at which depletion increases is chiefly dictated by the subduction rate, with faster subduction leading to a more rapid increase in depletion beneath the arc. Depletion within the sub-arc mantle is also found to increase with increasing mantle potential temperature, decreasing age of the overriding plate, and decreasing distance between the initial location of the BASC and the arc. Predicted changes in the depletion of the sub-arc mantle with time are shown to be consistent with observations of systematic along-strike geochemical variations within a portion of the Tonga Arc adjacent to the Eastern Lau Spreading Center

Mantle Melting as a Function of Water Content beneath the Mariana Arc

Subduction zone magmas are characterized by high concentrations of pre-eruptive H_2O, presumably as a result of an H_2Oflux originating from the dehydrating, subducting slab. The extent of mantle melting increases as a function of increasing water content beneath back-arc basins and is predicted to increase in a similar manner beneath arc volcanoes. Here, we present new data for olivine-hosted, basaltic melt inclusions from the Mariana arc that reveal pre-eruptive H_2O contents of ~1•5-6•0 wt %, which are up to three times higher than concentrations reported for the Mariana Trough back-arc basin. Major element systematics of arc and back-arc basin basalts indicate that the back-arc basin melting regime does not simply mix with wet, arc-derived melts to produce the observed range of back-arc magmatic H_2O concentrations. Simple melting models reveal that the trend of increasing extents of melting with increasing H_2O concentrations of the mantle source identified in the Mariana Trough generally extends beneath the Mariana volcanic front to higher mantle water contents and higher extents of melting. In detail, however, each Mariana volcano may define a distinct relationship between extent of melting and the H_2O content of the mantle source. We develop a revised parameterization of hydrous melting, incorporating terms for variable pressure and mantle fertility, to describe the distinct relationships shown by each arc volcano. This model is used in combination with thermobarometry constraints to show that hydrous melts equilibrate at greater depths (34-87 km) and temperatures (>1300°C) beneath the Mariana arc than beneath the back-arc basin (21-37 km), although both magma types can form from a mantle of similar potential temperature (~1350°C).The difference lies in where the melts form and equilibrate. Arc melts are dominated by those that equilibrate within the hot core of the mantle wedge, whereas back-arc melts are dominated by those that equilibrate within the shallow zone of decompression melting beneath the spreading center. Despite higher absolute melting temperatures (>1300°C), Mariana arc melts reflect lower melt productivity as a result of wet melting conditions and a more refractory mantle source

Reconciling mantle attenuation-temperature relationships from seismology, petrology, and laboratory measurements

Seismic attenuation measurements provide a powerful tool for sampling mantle properties. Laboratory experiments provide calibrations at seismic frequencies and mantle temperatures for dry melt-free rocks, but require ∼10²−10³ extrapolations in grain size to mantle conditions; also, the effects of water and melt are not well understood. At the same time, body wave attenuation measured from dense broadband arrays provides reliable estimates of shear wave attenuation (Q_S⁻¹), affording an opportunity for calibration. We reanalyze seismic data sets that sample arc and back-arc mantle in Central America, the Marianas, and the Lau Basin, confirming very high attenuation (Q_S ∼ 25–80) at 1 Hz and depths of 50–100 km. At each of these sites, independent petrological studies constrain the temperature and water content where basaltic magmas last equilibrated with the mantle, 1300–1450°C. The Q_S measurements correlate inversely with the petrologically inferred temperatures, as expected. However, dry attenuation models predict Q_S too high by a factor of 1.5–5. Modifying models to include effects of H₂O and rheology-dependent grain size shows that the effects of water-enhanced dissipation and water-enhanced grain growth nearly cancel, so H₂O effects are modest. Therefore, high H₂O in the arc source region cannot explain the low Q_S, nor in the back arc where lavas show modest water content. Most likely, the high attenuation reflects the presence of melt, and some models of melt effects come close to reproducing observations. Overall, body wave Q_S can be reconciled with petrologic and laboratory inferences of mantle conditions if melt has a strong influence beneath arcs and back arcs

Global variations in H_2O/Ce: 2. Relationships to arc magma geochemistry and volatile fluxes

We compiled a data set of 100 primitive arc magma compositions from melt inclusion and whole rock analyses to compare volatile contents, slab tracers, and calculated subduction component compositions between 18 subduction zone segments spanning the global range in slab thermal structure. The average primitive magma H_2O content in our data set is 3.3 ± 1.2 wt.% (1 s.d.) for melts erupted within 50 km of the volcanic front. While there is a wide range of volatile contents in magmas within individual arcs, the highest values occur in magmas erupted from vents along the volcanic front, where the subducting slab is located 104 ± 29 km (avg ± 1 s.d.) beneath the surface. This observation, coupled with positive correlations between H_2O, Cl, S, and B contents and predictions from geodynamic models, provides strong evidence for the active supply of volatile-rich slab-derived components from the subducting oceanic plate beneath volcanic arcs. We also show that temperature-sensitive ratios (e.g., H_2O/Ce) for both primitive arc magmas and calculated subduction components are similar and display monotonic behavior with slab thermal parameter. Furthermore, calculated subduction component compositions have higher trace element to H_2O ratios in arcs with hotter slabs (lower thermal parameter), suggesting that hydrous melts of differing compositions are added beneath different arcs. Finally, we present new volatile outflux estimates for Central Cascades magmatism and then compare these to estimates for the Central American and Kamchatka-Kurile arcs to create a combined data set spanning a large range in slab thermal parameter

Melting during late-stage rifting in Afar is hot and deep

Investigations of a variety of continental rifts and margins worldwide have revealed that a considerable volume of melt can intrude into the crust during continental breakup, modifying its composition and thermal structure. However, it is unclear whether the cause of voluminous melt production at volcanic rifts is primarily increased mantle temperature or plate thinning. Also disputed is the extent to which plate stretching or thinning is uniform or varies with depth with the entire continental lithospheric mantle potentially being removed before plate rupture. Here we show that the extensive magmatism during rifting along the southern Red Sea rift in Afar, a unique region of sub-aerial transition from continental to oceanic rifting, is driven by deep melting of hotter-than-normal asthenosphere. Petrogenetic modelling shows that melts are predominantly generated at depths greater than 80 kilometres, implying the existence of a thick upper thermo-mechanical boundary layer in a rift system approaching the point of plate rupture. Numerical modelling of rift development shows that when breakup occurs at the slow extension rates observed in Afar, the survival of a thick plate is an inevitable consequence of conductive cooling of the lithosphere, even when the underlying asthenosphere is hot. Sustained magmatic activity during rifting in Afar thus requires persistently high mantle temperatures, which would allow melting at high pressure beneath the thick plate. If extensive plate thinning does occur during breakup it must do so abruptly at a late stage, immediately before the formation of the new ocean basin

Tracking along-arc sediment inputs to the Aleutian arc using thallium isotopes.

Sediment transport from the subducted slab to the mantle wedge is an important process in understanding the chemical and physical conditions of arc magma generation. The Aleutian arc offers an excellent opportunity to study sediment transport processes because the subducted sediment flux varies systematically along strike (Kelemen et al., 2003) and many lavas exhibit unambiguous signatures of sediment addition to the sub-arc mantle (Morris et al., 1990). However, the exact sediment contribution to Aleutian lavas and how these sediments are transported from the slab to the surface are still debated. Thallium (Tl) isotope ratios have great potential to distinguish sediment fluxes in subduction zones because pelagic sediments and low-temperature altered oceanic crust are highly enriched in Tl and display heavy and light Tl isotope compositions, respectively, compared with the upper mantle and continental crust. Here, we investigate the Tl isotope composition of lavas covering almost the entire Aleutian arc a well as sediments outboard of both the eastern (DSDP Sites 178 and 183) and central (ODP Hole 886C) portions of the arc. Sediment Tl isotope compositions change systematically from lighter in the Eastern to heavier in the Central Aleutians reflecting a larger proportion of pelagic sediments when distal from the North American continent. Lavas in the Eastern and Central Aleutians mirror this systematic change to heavier Tl isotope compositions to the west, which shows that the subducted sediment composition is directly translated to the arc east of Kanaga Island. Moreover, quantitative mixing models of Tl and Pb, Sr and Nd isotopes reveal that bulk sediment transfer of ∼0.6–1.0% by weight in the Eastern Aleutians and ∼0.2–0.6% by weight in the Central Aleutians can account for all four isotope systems. Bulk mixing models, however, require that fractionation of trace element ratios like Ce/Pb, Cs/Tl, and Sr/Nd in the Central and Eastern Aleutians occurs after the sediment component was mixed with the mantle wedge. Models of Sr and Nd isotopes that involve sediment melting require either high degrees of sediment melting (>50%), in which case trace element ratios like Ce/Pb, Cs/Tl, and Sr/Nd of Aleutian lavas need to be produced after mixing with the mantle, or significant fluid additions from the underlying oceanic crust with Sr and Nd isotope compositions indistinguishable from the mantle wedge as well as high Sr/Nd ratios similar to that of low (<20%) degree sediment melts. Thallium isotope data from Western Aleutian lavas exhibit compositions slightly lighter than the upper mantle, which implies a negligible sediment flux at this location and probably involvement of low-temperature altered oceanic crust in the generation of these lavas. In general, the lightest Tl isotope compositions are observed for the highest Sr/Y ratios and most unradiogenic Sr and Pb isotope compositions, which is broadly consistent with derivation of these lavas via melting of eclogitized altered oceanic crust

Volatile loss from melt inclusions in pyroclasts of differing sizes

We have investigated the loss of H2O from olivine-hosted melt inclusions (MIs) by designing an experiment using tephra samples that cooled at different rates owing to their different sizes: ash, lapilli, and bomb samples that were deposited on the same day (10/17/74) of the sub-Plinian eruption of Volcán de Fuego in Guatemala. Ion microprobe, laser ablation-ICPMS, and electron probe analyses show that MIs from ash and lapilli record the highest H2O contents, up to 4.4 wt%. On the other hand, MIs from bombs indicate up to 30 % lower H2O contents (loss of ~1 wt% H2O) and 10 % post-entrapment crystallization of olivine. This evidence is consistent with the longer cooling time available for a bomb-sized clast, up to 10 min for a 3–4-cm radius bomb, assuming conductive cooling and the fastest H diffusivities measured in olivine (D~10−9 to 10−10 m2/s). On the other hand, several lines of evidence point to some water loss prior to eruption, during magma ascent and degassing in the conduit. Thus, results point to both slower post-eruptive cooling and slower magma ascent affecting MIs from bombs, leading to H2O loss over the timescale of minutes to hours. The important implication of this study is that a significant portion of the published data on H2O concentrations in olivine-hosted MIs may reflect unrecognized H2O loss via diffusion. This work highlights the importance of reporting clast and MI sizes in order to assess diffusive effects and the potential benefit of using water loss as a chronometer of magma ascent