CO2 content beneath northern Iceland and the variability of mantle carbon

© The Author(s), 2017. This article is distributed under the terms of the Creative Commons Attribution License. The definitive version was published in Geology 46 (2018): 55-58, doi:10.1130/G39413.1.Primitive basalt melt inclusions from Borgarhraun, northern Iceland, display large correlated variations in CO2 and nonvolatile incompatible trace elements (ITEs) such as Nb, Th, Rb, and Ba. The average CO2/ITE ratios of the Borgarhraun melt inclusion population are precisely determined (e.g., CO2/Nb = 391 ± 16; 2σM [two standard errors of the mean], n = 161). These data, along with published data on five other populations of undegassed mid-oceanic ridge basalt (MORB) glasses and melt inclusions, demonstrate that upper mantle CO2/Ba and CO2/Rb are nearly homogeneous, while CO2/Nb and CO2/Th are broadly correlated with long-term indices of mantle heterogeneity reflected in Nd isotopes (143Nd/144Nd) in five of the six regions of the upper mantle examined thus far. Our results suggest that heterogeneous carbon contents of the upper mantle are long-lived features, and that average carbon abundances of the mantle sources of Atlantic MORB are higher by a factor of two than those of Pacific MORB. This observation is correlated with a similar distinction in water contents and trace elements characteristic of subduction fluids (Ba, Rb). We suggest that the upper mantle beneath the younger Atlantic Ocean basin contains components of hydrated and carbonated subduction-modified mantle from prior episodes of Iapetus subduction that were entrained and mixed into the upper mantle during opening of the Atlantic Ocean basin.Maclennan is supported by Natural Environment Research Council grant NE/M000427/1. This research was supported by the Carnegie Institution of Washington

CO2 content beneath northern Iceland and the variability of mantle carbon

Primitive basalt melt inclusions from Borgarhraun, northern Iceland, display large correlated variations in CO2 and nonvolatile incompatible trace elements (ITEs) such as Nb, Th, Rb, and Ba. The average CO2/ITE ratios of the Borgarhraun melt inclusion population are precisely determined (e.g., CO2/Nb = 391 ± 16; 2σM [two standard errors of the mean], n = 161). These data, along with published data on five other populations of undegassed mid-oceanic ridge basalt (MORB) glasses and melt inclusions, demonstrate that upper mantle CO2/Ba and CO2/Rb are nearly homogeneous, while CO2/Nb and CO2/Th are broadly correlated with long-term indices of mantle heterogeneity reflected in Nd isotopes (143Nd/144Nd) in five of the six regions of the upper mantle examined thus far. Our results suggest that heterogeneous carbon contents of the upper mantle are long-lived features, and that average carbon abundances of the mantle sources of Atlantic MORB are higher by a factor of two than those of Pacific MORB. This observation is correlated with a similar distinction in water contents and trace elements characteristic of subduction fluids (Ba, Rb). We suggest that the upper mantle beneath the younger Atlantic Ocean basin contains components of hydrated and carbonated subduction-modified mantle from prior episodes of Iapetus subduction that were entrained and mixed into the upper mantle during opening of the Atlantic Ocean basin

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

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

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

The competing effects of sulfide saturation versus degassing on the behavior of the chalcophile elements during the differentiation of hydrous melts

There is a lack of consensus regarding the roles of sulfide saturation versus volatile degassing on the partitioning of Cu and Ag during differentiation and eruption of convergent margin magmas. Because of their oxidized character, volatile-rich magmas from the Eastern Manus Back-arc Basin (EMBB) only reach sulfide saturation following magnetite-driven reduction of the melt: the so-called “magnetite crisis.” If sulfide saturation typically precedes volatile saturation, the magnetite crisis will limit the proportion of Cu and Ag that can partition from the melt into an exsolving volatile-rich phase, which may contribute to the sporadic occurrence of magmatic-hydrothermal ore deposits at convergent margins. However, it is unclear whether the magnetite crisis is a common or rare event during differentiation of volatile-rich magmas. We report major and trace element data for submarine volcanic glasses from the Tonga arc-proximal Valu Fa Ridge (VFR; SW Pacific). Cu-Se-Ag systematics of samples erupting at the southern VFR suggest magnetite fractionation-triggered sulfide saturation. The similarity in chalcophile element systematics of the southern VFR and EMBB samples is unlikely to be coincidental, and may indicate that the magnetite crisis is a common event during differentiation of hydrous melts. However, unlike many convergent margin magmas, it is unlikely that the evolving VFR and EMBB were saturated in a S-bearing volatile phase prior to magnetite fractionation. Hence, the metal-depleting magnetite crisis may be restricted to back-arc basin magmas that do not degas volatiles prior to magnetite fractionation and potentially convergent margin magmas fractionating at high pressures in the continental crust

Asymmetric shallow mantle structure beneath the Hawaiian Swell—evidence from Rayleigh waves recorded by the PLUME network

Author Posting. © The Author(s), 2011. This article is posted here by permission of John Wiley & Sons for personal use, not for redistribution. The definitive version was published in Geophysical Journal International 187 (2011): 1725–1742, doi:10.1111/j.1365-246X.2011.05238.x.We present models of the 3-D shear velocity structure of the lithosphere and asthenosphere beneath the Hawaiian hotspot and surrounding region. The models are derived from long-period Rayleigh-wave phase velocities that were obtained from the analysis of seismic recordings collected during two year-long deployments for the Hawaiian Plume-Lithosphere Undersea Mantle Experiment. For this experiment, broad-band seismic sensors were deployed at nearly 70 seafloor sites as well as 10 sites on the Hawaiian Islands. Our seismic images result from a two-step inversion of path-averaged dispersion curves using the two-station method. The images reveal an asymmetry in shear velocity structure with respect to the island chain, most notably in the lower lithosphere at depths of 60 km and greater, and in the asthenosphere. An elongated, 100-km-wide and 300-km-long low-velocity anomaly reaches to depths of at least 140 km. At depths of 60 km and shallower, the lowest velocities are found near the northern end of the island of Hawaii. No major velocity anomalies are found to the south or southeast of Hawaii, at any depth. The low-velocity anomaly in the asthenosphere is consistent with an excess temperature of 200–250 °C and partial melt at the level of a few percent by volume, if we assume that compositional variations as a result of melt extraction play a minor role. We also image small-scale low-velocity anomalies within the lithosphere that may be associated with the volcanic fields surrounding the Hawaiian Islands.This research was financed by the National Science Foundation under grants OCE-00-02470 and OCE-00-02819. Markee was partly sponsored by a SIO graduate student fellowship

MPI-DING reference glasses for in situ microanalysis: New reference values for element concentrations and isotope ratios

We present new analytical data of major and trace elements for the geological MPI-DING glasses KL2-G, ML3B-G, StHs6/80-G, GOR128-G, GOR132-G, BM90/21-G, T1-G, and ATHO-G. Different analytical methods were used to obtain a large spectrum of major and trace element data, in particular, EPMA, SIMS, LA-ICPMS, and isotope dilution by TIMS and ICPMS. Altogether, more than 60 qualified geochemical laboratories worldwide contributed to the analyses, allowing us to present new reference and information values and their uncertainties (at 95% confidence level) for up to 74 elements. We complied with the recommendations for the certification of geological reference materials by the International Association of Geoanalysts (IAG). The reference values were derived from the results of 16 independent techniques, including definitive (isotope dilution) and comparative bulk (e.g., INAA, ICPMS, SSMS) and microanalytical (e.g., LA-ICPMS, SIMS, EPMA) methods. Agreement between two or more independent methods and the use of definitive methods provided traceability to the fullest extent possible. We also present new and recently published data for the isotopic compositions of H, B, Li, O, Ca, Sr, Nd, Hf, and Pb. The results were mainly obtained by high-precision bulk techniques, such as TIMS and MC-ICPMS. In addition, LA-ICPMS and SIMS isotope data of B, Li, and Pb are presented. Copyright 2006 by the American Geophysical Union

Lithosphere versus asthenosphere mantle sources at the Big Pine Volcanic Field, California

[1] Here we report the first measurements of the H2O content of magmas and mantle xenoliths from the Big Pine Volcanic Field (BPVF), California, in order to constrain the melting process in the mantle, and the role of asthenospheric and lithospheric sources in this westernmost region of the Basin and Range Province, western USA. Melt inclusions trapped in primitive olivines (Fo82–90) record surprisingly high H2O contents (1.5 to 3.0 wt.%), while lithospheric mantle xenoliths record low H2O concentrations (whole rock <75 ppm). Estimates of the oxidation state of BPVF magmas, based on V partitioning in olivine, are also high (FMQ +1.0 to +1.5). Pressures and temperatures of equilibration of the BPVF melts indicate a shift over time, from higher melting temperatures (∼1320°C) and pressures (∼2 GPa) for magmas that are >500 ka, to cooler (∼1220°C) and shallower melting (∼1 GPa) conditions in younger magmas. The estimated depth of melting correlates strongly with some trace element ratios in the magmas (e.g., Ce/Pb, Ba/La), with deeper melts having values closer to upper mantle asthenosphere values, and shallower melts having values more typical of subduction zone magmas. This geochemical stratification is consistent with seismic observations of a shallow lithosphere-asthenosphere boundary (∼55 km depth). Combined trace element and cryoscopic melting models yield self-consistent estimates for the degree of melting (∼5%) and source H2O concentration (∼1000 ppm). We suggest two possible geodynamic models to explain small-scale convection necessary for magma generation. The first is related to the Isabella seismic anomaly, either a remnant of the Farallon Plate or foundered lithosphere. The second scenario is related to slow extension of the lithosphere

MPI-Ding reference glasses for in situ microanalysis: New reference values for element concentrations and isotope ratios

We present new analytical data of major and trace elements for the geological MPI-DING glasses KL2-G, ML3B-G, StHs6/80-G, GOR128-G, GOR132-G, BM90/21-G, T1-G, and ATHO-G. Different analytical methods were used to obtain a large spectrum of major and trace element data, in particular, EPMA, SIMS, LA-ICPMS, and isotope dilution by TIMS and ICPMS. Altogether, more than 60 qualified geochemical laboratories worldwide contributed to the analyses, allowing us to present new reference and information values and their uncertainties (at 95% confidence level) for up to 74 elements. We complied with the recommendations for the certification of geological reference materials by the International Association of Geoanalysts (IAG). The reference values were derived from the results of 16 independent techniques, including definitive (isotope dilution) and comparative bulk (e.g., INAA, ICPMS, SSMS) and microanalytical (e.g., LA-ICPMS, SIMS, EPMA) methods. Agreement between two or more independent methods and the use of definitive methods provided traceability to the fullest extent possible. We also present new and recently published data for the isotopic compositions of H, B, Li, O, Ca, Sr, Nd, Hf, and Pb. The results were mainly obtained by high-precision bulk techniques, such as TIMS and MC-ICPMS. In addition, LA-ICPMS and SIMS isotope data of B, Li, and Pb are presented