An experimental study of crystalline basalt dissolution from 2 ? pH ? 11 and temperatures from 5 to 75 °C

Steady-state element release rates from crystalline basalt dissolution at far-from-equilibrium were measured at pH from 2 to 11 and temperatures from 5 to 75 °C in mixed-flow reactors. Steady-state Si and Ca release rates exhibit a U-shaped variation with pH where rates decrease with increasing pH at acid condition but increase with increasing pH at alkaline conditions. Silicon release rates from crystalline basalt are comparable to Si release rates from basaltic glass of the same chemical composition at low pH and temperatures ?25 °C but slower at alkaline pH and temperatures ?50 °C. In contrast, Mg and Fe release rates decrease continuously with increasing pH at all temperatures. This behaviour is interpreted to stem from the contrasting dissolution behaviours of the three major minerals comprising the basalt: plagioclase, pyroxene, and olivine. Calcium is primarily present in plagioclase, which exhibits a U-shaped dissolution rate dependence on pH. In contrast, Mg and Fe are contained in pyroxene and olivine, minerals whose dissolution rates decrease monotonically with pH. As a result, crystalline basalt preferentially releases Mg and Fe relative to Ca at acidic conditions. The injection of acidic CO₂-charged fluids into crystalline basaltic terrain may, therefore, favour the formation of Mg and Fe carbonates rather than calcite. Element release rates estimated from the sum of the volume fraction normalized dissolution rates of plagioclase, pyroxene, and olivine are within one order of magnitude of those measured in this study.     

Dissolution of basalts and peridotite in seawater, in the presence of ligands, and CO2: Implications for mineral sequestration of carbon dioxide

Steady-state silica release rates (r_Si) from basaltic glass and crystalline basalt of similar chemical composition as well as dunitic peridotite have been determined in far-from-equilibrium dissolution experiments at 25 °C and pH 3.6 in (a) artificial seawater solutions under 4 bar pCO₂, (b) varying ionic strength solutions, including acidified natural seawater, (c) acidified natural seawater of varying fluoride concentrations, and (d) acidified natural seawater of varying dissolved organic carbon concentrations. Glassy and crystalline basalts exhibit similar r_Si in solutions of varying ionic strength and cation concentrations. Rates of all solids are found to increase by 0.3-0.5 log units in the presence of a pCO₂ of 4 bar compared to CO₂ pressure of the atmosphere. At atmospheric CO₂ pressure, basaltic glass dissolution rates were most increased by the addition of fluoride to solution whereas crystalline basalt rates were most enhanced by the addition of organic ligands. In contrast, peridotite does not display any significant ligand-promoting effect, either in the presence of fluoride or organic acids. Most significantly, Si release rates from the basalts are found to be not more than 0.6 log units slower than corresponding rates of the peridotite at all conditions considered in this study. This difference becomes negligible in seawater suggesting that for the purposes of in-situ mineral sequestration, CO₂-charged seawater injected into basalt might be nearly as efficient as injection into peridotite.     

Generation of CO2-rich melts during basalt magma ascent and degassing

To test mechanisms of basaltic magma degassing, continuous decompressions of volatile-bearing (2.7–3.8 wt% H₂O, 600–1,300 ppm CO₂) Stromboli melts were performed from 250–200 to 50–25 MPa at 1,180–1,140 °C. Ascent rates were varied from 0.25 to ~1.5 m/s. Glasses after decompression show a wide range of textures, from totally bubble-free to bubble-rich, the latter with bubble number densities from 10⁴ to 10⁶ cm^?3, similar to Stromboli pumices. Vesicularities range from 0 to ~20 vol%. Final melt H₂O concentrations are homogeneous and always close to solubilities. In contrast, the rate of vesiculation controls the final melt CO₂ concentration. High vesicularity charges have glass CO₂ concentrations that follow theoretical equilibrium degassing paths, whereas glasses from low vesicularity charges show marked deviations from equilibrium, with CO₂ concentrations up to one order of magnitude higher than solubilities. FTIR profiles and maps reveal glass CO₂ concentration gradients near the gas–melt interface. Our results stress the importance of bubble nucleation and growth, and of volatile diffusivities, for basaltic melt degassing. Two characteristic distances, the gas interface distance (distance either between bubbles or to gas–melt interfaces) and the volatile diffusion distance, control the degassing process. Melts containing numerous and large bubbles have gas interface distances shorter than volatile diffusion distances, and degassing proceeds by equilibrium partitioning of CO₂ and H₂O between melt and gas bubbles. For melts where either bubble nucleation is inhibited or bubble growth is limited, gas interface distances are longer than volatile diffusion distances. Degassing proceeds by diffusive volatile transfer at the gas–melt interface and is kinetically limited by the diffusivities of volatiles in the melt. Our experiments show that CO₂-oversaturated melts can be generated as a result of magma decompression. They provide a new explanation for the occurrence of CO₂-rich natural basaltic glasses and open new perspectives for understanding explosive basaltic volcanism.     

Melt inclusion constraints on volatile systematics and degassing history of the 2014–2015 Holuhraun eruption,Iceland

The mass of volatiles emitted during volcanic eruptions is often estimated by comparing the volatile contents of undegassed melt inclusions, trapped in crystals at an early stage of magmatic evolution, with that of the degassed matrix glass. Here we present detailed characterisation of magmatic volatiles (H₂O, CO₂, S, Fl and Cl) of crystal-hosted melt and fluid inclusions from the 2014–2015 Holuhraun eruption of the Bárðarbunga volcanic system, Iceland. Based on the ratios of magmatic volatiles to similarly incompatible trace elements, the undegassed primary volatile contents of the Holuhraun parental melt are estimated at 1500–1700 ppm CO₂, 0.13–0.16 wt% H₂O, 60–80 ppm Cl, 130–240 ppm F and 500–800 ppm S. High-density fluid inclusions indicate onset of crystallisation at pressures?≥?0.4 GPa (~?12 km depth) promoting deep degassing of CO₂. Prior to the onset of degassing, the melt CO₂ content may have reached 3000–4000 ppm, with the total magmatic CO₂ budget estimated at  23–55 Mt. SO₂ release commenced at 0.12 GPa (~?3.6 km depth), eventually leading to entrapment of SO₂ vapour in low-density fluid inclusions. We calculate the syn-eruptive volatile release as 22.2 Mt of magmatic H₂O, 5.9–7.7 Mt CO₂, and 11.3 Mt of SO₂ over the course of the eruption; F and Cl release were insignificant. Melt inclusion constraints on syn-eruptive volatile release are similar to estimates made during in situ field monitoring, with the exception of H₂O, where field measurements may be heavily biased by the incorporation of meteoric water.     

CO2 solubility in Martian basalts and Martian atmospheric evolution

To understand possible volcanogenic fluxes of CO₂ to the Martian atmosphere, we investigated experimentally carbonate solubility in a synthetic melt based on the Adirondack-class Humphrey basalt at 1-2.5 GPa and 1400-1625 °C. Starting materials included both oxidized and reduced compositions, allowing a test of the effect of iron oxidation state on CO₂ solubility. CO₂ contents in experimental glasses were determined using Fourier transform infrared spectroscopy (FTIR) and Fe³⁺/Fe^T was measured by Mössbauer spectroscopy. The CO₂ contents of glasses show no dependence on Fe³⁺/Fe^T and range from 0.34 to 2.12 wt.%. For Humphrey basalt, analysis of glasses with gravimetrically-determined CO₂ contents allowed calibration of an integrated molar absorptivity of 81,500 ± 1500 L mol⁻¹ cm⁻² for the integrated area under the carbonate doublet at 1430 and 1520 cm⁻¹. The experimentally determined CO₂ solubilities allow calibration of the thermodynamic parameters governing dissolution of CO₂ vapor as carbonate in silicate melt, K_II, (Stolper and Holloway, 1988) as follows: , ΔV⁰ = 20.85 ± 0.91 cm³ mol⁻¹, and ΔH⁰ = −17.96 ± 10.2 kJ mol⁻¹. This relation, combined with the known thermodynamics of graphite oxidation, facilitates calculation of the CO₂ dissolved in magmas derived from graphite-saturated Martian basalt source regions as a function of P, T, and f_O2. For the source region for Humphrey, constrained by phase equilibria to be near 1350 °C and 1.2 GPa, the resulting CO₂ contents are 51 ppm at the iron-wüstite buffer (IW), and 510 ppm at one order of magnitude above IW (IW + 1). However, solubilities are expected to be greater for depolymerized partial melts similar to primitive shergottite Yamato 980459 (Y 980459). This, combined with hotter source temperatures (1540 °C and 1.2 GPa) could allow hot plume-like magmas similar to Y 980459 to dissolve 240 ppm CO₂ at IW and 0.24 wt.% of CO₂ at IW + 1. For expected magmatic fluxes over the last 4.5 Ga of Martian history, magmas similar to Humphrey would only produce 0.03 and 0.26 bars from sources at IW and IW + 1, respectively. On the other hand, more primitive magmas like Y 980459 could plausibly produce 0.12 and 1.2 bars at IW and IW + 1, respectively. Thus, if typical Martian volcanic activity was reduced and the melting conditions cool, then degassing of CO₂ to the atmosphere may not be sufficient to create greenhouse conditions required by observations of liquid surface water. However, if a significant fraction of Martian magmas derive from hot and primitive sources, as may have been true during the formation of Tharsis in the late Noachian, that are also slightly oxidized (IW + 1.2), then significant contribution of volcanogenic CO₂ to an early Martian greenhouse is plausible.     

Investigation of datolite (CaB[SiO4/(OH)]) from basalts in the Northern Apennines ophiolites (Italy): Genetic implications

Datolite, ideally CaB[(OH/SiO₄)], from hydrothermal veins crosscutting pillow basalt in 10 different localities of the Northern Apennine ophiolites was investigated with regard to mineral chemistry and fluid inclusion microthermometry. Bulk analyses of datolite crystals show REE contents below chondritic, except for La and Ce. With respect to host rock, datolite is occasionally enriched in La, Rb, Cs, Be, and shows relatively high contents of chalcophile elements (Cu, Zn, Pb, Ni) when occurring in contact with sulfide-mineralized basalt. Volatiles escaped during the decomposition in the temperature range 600 and 700 °C. The main component is water. The temperature maximum of water release is different and frequently with a shoulder or a second maximum. Together with water, sulfur species as H₂S and SO₂ and traces of boron species escaped. The CO₂ release by the decomposition especially of datolite from Castellaro and Cinghi has a maximum in the range of 500-580 °C and is different from the decomposition of calcite. Together with CO₂ a boron species escaped. Chlorine does not detect. Two-phase (L+V) fluid inclusions texturally identifiable as primary and secondary were observed, yielding average homogenization-temperatures of 236 and 173 °C, respectively. Fluid inclusion cooling data yield calculated salinity in the range of 10-16 wt% NaCl equivalents, thus relatively higher compared with seawater. The results are compatible with those reported for fluids formed under diagenetic conditions, but differ from those observed in seafloor hydrothermal systems and/or emanating from magmas. Distribution of trace elements between datolite and host basalt indicates enrichment with respect to the host rock limited to a few elements such as La, Rb, Cs, Be, Ni, Cu, Zn and Pb. The lithophile elements can be hosted in the datolite lattice, whereas the chalcophile metals and Ni are probably carried in sub-microscopic inclusions.     

Experimental constraints on the mobility of Rhenium in silicate liquids

The volatization of Rhenium (Re) from melts of natural basalt, dacite and a synthetic composition in the CaO-MgO-Al₂O₃-SiO₂ system has been investigated at 0.1 MPa and 1250-1350 °C over a range of fO₂ conditions from log fO₂ = −10 to −0.68. Experiments were conducted using open top Pt crucibles doped with Re and Yb. Analysis of quenched glasses by laser ablation-inductively coupled plasma mass spectrometry (LA-ICP-MS) normal to the melt/gas interface showed concentration profiles for Re, to which a semi-infinite one-dimensional diffusion model could be applied to extract diffusion coefficients (D). The results show Re diffusivity in basalt at 1300 °C in air is log D_Re = −7.2 ± 0.3 cm²/s and increases to log D_Re = −6.6 ± 0.3 cm²/s when trace amounts of Cl were added to the starting material. At fO₂ conditions below the nickel-nickel oxide (NNO) buffer Re diffusivity decreases to and to in dacitic melt. In the CMAS composition, . The diffusivity of Re is comparable to Ar and CO₂ in basalt at 500 MPa favoring its release as a volatile. Our results support the contention that subaerial degassing is the cause of lower Re concentrations in arc-type and ocean island basalts compared to mid-ocean ridge basalts.     

Gold solubility in oxidized and reduced, water-saturated mafic melt

We have performed experiments to evaluate Au solubility in natural, water-saturated basaltic melts as a function of oxygen fugacity. Experiments were carried out at 1000 °C and 200 MPa, and oxygen fugacity was controlled at the fayalite-magnetite-quartz (FMQ) oxygen fugacity buffer and FMQ + 4. All experiments were saturated with a metal-chloride aqueous solution loaded initially as a 10 wt% NaCl eq. fluid. The stable phase assemblage at FMQ consists of basalt melt, olivine, clinopyroxene, a single-phase aqueous fluid, and metallic Au. The stable phase assemblage at FMQ + 4 consists of basalt melt, clinopyroxene, magnetite-spinel solid solution, a single-phase aqueous fluid, and metallic Au. Silicate glasses (i.e., quenched melt) and their contained crystalline material were analyzed by using both electron probe microanalysis (EPMA) and laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS). Measured Au concentrations in the quenched melt range from 4.8 μg g⁻¹ to 0.64 μg g⁻¹ at FMQ + 4, and 0.54 μg g⁻¹ to 0.1 μg g⁻¹ at FMQ. The measured solubility of Au in olivine and clinopyroxene was consistently below the LA-ICP-MS limit of detection (i.e., 0.1 μg g⁻¹). These melt solubility data place important limitations on the dissolved Au content of water-saturated, Cl- and S-bearing basaltic liquids at geologically relevant fO₂ values. The new data are compared to published, experimentally-determined values for Au solubility in dry and hydrous silicate liquids spanning the compositional range from basalt to rhyolite, and the effects of melt composition, oxygen fugacity, pressure and temperature are discussed.     

Temperature dependence of sulfide and sulfate solubility in olivine-saturated basaltic magmas

The sulfur concentration at pyrrhotite- and anhydrite-saturation in primitive hydrous basaltic melt of the 2001-2002 eruption of Mt. Etna was determined at 200 MPa, T = 1050-1250 °C and at log fO₂ from FMQ to FMQ+2.2 (FMQ is Fayalite-Magnetite-Quartz oxygen buffer). At 1050 °C Au sample containers were used. A double-capsule technique, using a single crystal olivine sample container closed with an olivine piston, embedded in a sealed Au₈₀Pd₂₀ capsule, was developed to perform experiments in S-bearing hydrous basaltic systems at T > 1050 °C. Pyrrhotite is found to be a stable phase coexisting with melt at FMQ-FMQ+0.3, whereas anhydrite is stable at FMQ+1.4-FMQ+2.2. The S concentration in the melt increases almost linearly from 0.12 ± 0.01 to 0.39 ± 0.02 wt.% S at FeS-saturation and from 0.74 ± 0.01 to 1.08 ± 0.04 wt.% S at anhydrite-saturation with T ranging from 1050-1250 °C. The relationships between S concentration at pyrrhotite and/or anhydrite saturation, MgO content of the olivine-saturated melt, T, and log fO₂ observed in this study and from previous data are used to develop an empirical model for estimating the magmatic T and fO₂ from the S and MgO concentrations of H₂O-bearing olivine-saturated basaltic melts. The model can also be used to determine maximum S concentrations, if fO₂ and MgO content of the melt are known. The application of the model to compositions of melt inclusions in olivines from Mt. Etna indicates that the most primitive magmas trapped in inclusions might have been stored at log fO₂ slightly higher than FMQ+1 and at T = 1100-1150 °C, whereas more evolved melts could have been trapped at T ? 1100 °C. These values are in a good agreement with the estimates obtained by other independent methods reported in the literature.     

Time-dependent changes of the electrical conductivity of basaltic melts with redox state

The electrical conductivity of basaltic melts has been measured in real-time after fO₂ step-changes in order to investigate redox kinetics. Experimental investigations were performed at 1 atm in a vertical furnace between 1200 and 1400 °C using air, pure CO₂ or CO/CO₂ gas mixtures to buffer oxygen fugacity in the range 10⁻⁸ to 0.2 bars. Ferric/ferrous ratios were determined by wet chemical titrations. A small but detectable effect of fO₂ on the electrical conductivity is observed. The more reduced the melt, the higher the conductivity. A modified Arrhenian equation accounts for both T and fO₂ effects on the electrical conductivity. We show that time-dependent changes in electrical conductivity following fO₂ step-changes monitor the rate of Fe²⁺/Fe³⁺ changes. The conductivity change with time corresponds to a diffusion-limited process in the case of reduction in CO-CO₂ gas mixtures and oxidation in air. However, a reaction at the gas-melt interface probably rate limits oxidation of the melt under pure CO₂. Reduction and oxidation rates are similar and both increase with temperature. Those rates range from 10⁻⁹ to 10⁻⁸ m²/s for the temperature interval 1200-1400 °C and show activation energy of about 200 kJ/mol. The redox mechanism that best explains our results involves a cooperative motion of cations and oxygen, allowing such fast oxidation-reduction rates.     

Experimental determination of oxygen isotope fractionations between CO2 vapor and soda-melilite melt

We report results of experiments constraining oxygen isotope fractionations between CO₂ vapor and Na-rich melilitic melt at 1 bar and 1250 and 1400°C. The fractionation factor constrained by bracketed experiments, 1000^.lnα_{CO2-Na melilitic melt}, is 2.65±0.25 ‰ (±2σ; n=92) at 1250°C and 2.16±0.16 ‰ (2σ; n=16) at 1400°C. These values are independent of Na content over the range investigated (7.5 to 13.0 wt. % Na₂O). We combine these data with the known reduced partition function ratio of CO₂ to obtain an equation describing the reduced partition function ratio of Na-rich melilite melt as a function of temperature. We also fit previously measured CO₂-melt or -glass fractionations to obtain temperature-dependent reduced partition function ratios for all experimentally studied melts and glasses (including silica, rhyolite, albite, anorthite, Na-rich melilite, and basalt). The systematics of these data suggest that reduced partition function ratios of silicate melts can be approximated either by using the Garlick index (a measure of the polymerization of the melt) or by describing melts as mixtures of normative minerals or equivalent melt compositions. These systematics suggest oxygen isotope fractionation between basalt and olivine at 1300°C of approximately 0.4 to 0.5‰, consistent with most (but not all) basalt glass-olivine fractionations measured in terrestrial and lunar basalts.     

Carbonate xenoliths in La Palma: Carbonatite or alteration product?

Carbonate xenoliths containing olivine and rimmed by kaersutitic amphibole were collected in basaltic rocks of the Basal Complex of La Palma. The mineralogical composition and microscopic appearance may suggest a relationship with carbonatites in general, thus a major element, trace element and stable isotope study was conducted to investigate the origin of the carbonate formation. Based on electron microprobe analyses, the carbonate is calcite with up to 6.3 wt% MgO and 7.2 wt% SiO₂. The elevated SiO₂ content may suggest a melt origin for the carbonate. However, the C and O isotope compositions of the carbonate xenoliths (δ¹³C and δ¹⁸O around −1‰ and 13‰, respectively) are similar to those of calcite veins and amygdales in basaltic rocks of the Basal Complexes of La Palma and Fuerteventura and are interpreted as produced by fluid degassing and metasomatism by CO₂-H₂O fluid derived from mobilization of sedimentary material. Trace element contents determined by laser-ablation ICP-MS analyses support the assumed origin, thus, the relationship with carbonatitic melts can be excluded. Based on trace element compositions, the amphibole surrounding the xenoliths is not related genetically to the carbonate. The elevated SiO₂ content of the calcite can be attributed to submicron relics of pyroxene, thus, the use of this feature as an evidence for melt origin is questionable.     

Carbon dioxide in silicate melts: A molecular dynamics simulation study

We have performed a series of molecular dynamics simulations aimed at the evaluation of the solubility of CO₂ in silicate melts of natural composition (from felsic to ultramafic). In making in contact within the simulation cell a supercritical CO₂ phase with a silicate melt of a given composition, we have been able to evaluate (i) the solubility of CO₂ in the P-T range 1473-2273 K and 20-150 kbar, (ii) the density change experienced by the CO₂-bearing melt, (iii) the respective concentrations of CO₂ and species in the melt, (iv) the lifetime and the diffusivity of these species and (v) the structure of the melt around the carbonate groups. The main results are the following:(1) The solubility of CO₂ increases markedly with the pressure in the three investigated melts (a rhyolite, a mid-ocean ridge basalt and a kimberlite) from about ∼2 wt% CO₂ at 20 kbar to ∼25 wt% at 100 kbar and 2273 K. The solubility is found to be weakly dependent on the melt composition (as far as the present compositions are concerned) and it is only at very high pressure (above ∼100 kbar) that a clear hierarchy between solubilities occurs (rhyolite < MORB < kimberlite). Furthermore at a given pressure the calculated solubility is negatively correlated with the temperature.(2) In CO₂-saturated melts, the proportion of carbonate ions is positively correlated with the pressure at isothermal condition and is negatively correlated with the temperature at isobaric condition (and vice versa for molecular CO₂). Furthermore, at fixed (P, T) conditions the proportion of carbonate ions is higher in CO₂-undersaturated melts than in the CO₂-saturated melt. Although the proportion of molecular CO₂ decreases when the degree of depolymerization of the melt increases, it is still significant in CO₂-saturated basic and ultrabasic compositions at high temperatures. This finding is at variance with experimental data on CO₂-bearing glasses which show no evidence of molecular CO₂ as soon as the degree of depolymerization of the melt is high (e.g. basalt). These conflicting results can be reconciled with each other by noticing that a simple low temperature extrapolation of the simulation data predicts that the proportion of molecular CO₂ in basaltic melts might be negligible in the glass at room temperature.(3) The carbonate ions are found to be transient species in the liquid phase, with a lifetime increasing exponentially with the inverse of the temperature. Contrarily to a usual assumption, the diffusivity of carbonate ions into the liquid silicate is not vanishingly small with respect to that of CO₂ molecules: in MORB they differ from each other by a factor of ∼6 at 1473 K and only a factor of ∼2 at 2273 K. Although the bulk diffusivity of CO₂ is governed primarily by the diffusivity of CO₂ molecules, the carbonate ions contribute significantly to the diffusivity of CO₂ in depolymerized melts.(4) Concerning the structure of the CO₂-bearing silicate melt, the carbonate ions are found to be preferentially associated with NBO’s of the melt, with an affinity for NBOs which exceeds that for BOs by almost one order of magnitude. This result explains why the concentration in carbonate ions is positively correlated with the degree of depolymerization of the melt and diminishes drastically in fully polymerized melts where the number of NBO’s is close to zero. Furthermore, the network modifier cations are not randomly distributed in the close vicinity of carbonate groups but exhibit a preferential ordering which depends at once on the nature of the cation and on the melt composition. However at the high temperatures investigated here, there is no evidence of long lived complexes between carbonate groups and metal cations.     

Experimental studies of oxygen isotope fractionation between rhodochrosite (MnCO3) and water at low temperatures

Rhodochrosite crystals were precipitated from Na-Mn-Cl-HCO₃ parent solutions following passive, forced and combined passive-to-forced CO₂ degassing methods. Forced and combined passive-to-forced CO₂ degassing produced rhodochrosite crystals with a small non-equilibrium oxygen isotope effect whereas passive CO₂ degassing protocols yielded rhodochrosite in apparent isotopic equilibrium with water. On the basis of the apparent equilibrium isotopic data, a new temperature-dependent relation is proposed for the oxygen isotope fractionation between rhodochrosite and water between 10 and 40 °C:

1000lnα_{rhodochrosite-water}=17.84±0.18(10³/T)-30.24±0.62     

Effect of oxygen fugacity on the H2O storage capacity of forsterite in the carbon-saturated systems

High pressure experiments have been performed in the systems Mg₂SiO₄-C-O-H and Mg₂SiO₄-K₂CO₃-C at 6.3 GPa and 1200 to 1600 °C using a split-sphere multi-anvil apparatus. In the Mg₂SiO₄-C-O-H system the composition of fluid was modeled by adding different amounts of water and stearic acid. The fO₂ was controlled by the Mo-MoO₂ or Fe-FeO oxygen buffers. Several experiments in the Mg₂SiO₄-C-O-H system and all experiments in the Mg₂SiO₄-K₂CO₃-C system have been conducted without buffering the fO₂. Forsterite in the system Mg₂SiO₄-K₂CO₃-C does not reveal OH absorption bands in the IR spectra, while forsterite coexisting with carbon-bearing fluid and silicate melt at logfO₂ from FMQ-2 to FMQ-5 (from 2 to 5 log units below fayalite-magnetite-quartz oxygen buffer) contains 800-1850 wt. ppm H₂O. The maximum concentrations were detected at 1400 °C and FMQ-3.5. We observed an increase in the solidus temperature in the system Mg₂SiO₄-C-O-H from 1200 to above 1600 °C with log fO₂ decreasing from FMQ-2 to FMQ-5. The increase of the solidus temperature and the broadening of the stability field of the H₂O-H₂-CH₄ subsolidus fluid phase at 1400-1600 °C explain the high H₂O storage capacity of forsterite relative to that crystallized from carbon-free, oxidized, hydrous, silicic melt. At temperatures above 1400 °C liquidus forsterite precipitated along with diamond from oxidized (FMQ-1) carbonate-silicate melt and from silicate melt dissolving the moderately reduced C-O-H fluid (from FMQ-2 to FMQ-3.5). Formation of diamond was not detected under ultra-reduced conditions (FMQ-5) at 1200-1600 °C. Olivine co-precipitating with diamond from dry carbonate-silicate or hydrous-silicic fluid/melt can provide information on the H₂O contents and speciation of the diamond-forming media in the mantle. The conditions for minimum post-crystallization alteration of olivine and its hydrogen content are discussed.     

Basalt contamination by continental crust: Some experiments and models

Chemical interaction between molten basalt and felsic minerals of the continental crust (quartz, K-feldspar, and oligoclase) was examined in static and dynamic experiments at 1,200°–1,400° C. Under circumstances of continuous stirring at 1,400°, -quartz dissolves in tholeiite melt at a rate of 3.3×10^–6 g/s per cm² of contact area; at 1,300°, the solution rate is 1.5×10^–6 g/cm{cm2}s. The feldspars are molten at the experimental conditions, and interact with contacting basalt melt by diffusion in the liquid state. This is a complex process characterized by rapid initial diffusion of alkalies to establish a distribution between felsic melt and basalt similar to that observed in cases of actual two-liquid equilibrium (both alkalies reach concentrations in the felsic melt 1.5–3 times those in the basalt). Alkali diffusion may be uphill or downhill, depending on which direction of net flux is required to produce a two-liquid type distribution. Once this distribution is attained, subsequent diffusion of all melt species is slow and apparently limited by the diffusivity of SiO₂, which is 10^–9-l0^–10 cm²/s at 1,200° C. Interdiffusion experiments involving molten basalt and synthetic granite confirm the behavior illustrated by the feldspar/basalt results, and give similar SiO₂ diffusivities.The solution rates and interdiffusion data can be used to model basalt contamination processes likely to occur in the continental crust. For the restricted case of solid quartzitic xenoliths, the uptake of SiO₂ in a well-mixed basalt magma is quite fast: appreciable SiO₂ contamination may occur over exposure times of only days to years. If basalt magma induces local melting of crustal rocks, the assimilation process becomes one of liquid-state interdiffusion. In this case, the varying diffusivities of ions and their differing preferences for silicic relative to basaltic melts can produce marked selective contamination effects. Selective contamination of ascending basaltic magmas is particularly likely in the case of K₂O, which may be introduced in substantial amounts even when other elements remain unaffected. The Na₂O content of mantle-derived magmas is buffered against contamination by crustal materials, and K₂O is buffered against further increases once it reaches a level of 1–1.5 wt.%.     

Halogen diffusion in a basaltic melt

The diffusion of the halogens fluorine, chlorine and bromine was measured in a hawaiitic melt from Mt. Etna at 500 MPa and 1.0 GPa, 1250 to 1450 °C at anhydrous conditions; the diffusion of F and Cl in the melt was also studied with about 3 wt% of dissolved water. Experiments were performed using the diffusion-couple technique in a piston cylinder. Most experiments were performed with only one halogen diffusing between the halogen-enriched and halogen-poor halves of the diffusion couple, but a few experiments with a mixture of halogens (F, Cl and Br) were also performed in order to investigate the possibility of interactions between the halogens during diffusion. Fluorine and chlorine diffusivity show a very similar behavior, slightly diverging at low temperature. Bromine diffusion is a factor of about 2-5 lower than the other halogens in this study. Diffusion coefficients for fluorine range between 2.3 × 10⁻¹¹ and 1.4 × 10⁻¹⁰ m² s⁻¹, for chlorine between 1.1 × 10⁻¹¹ and 1.3 × 10⁻¹⁰ and for bromine between 9.4 × 10⁻¹² and 6.8 × 10⁻¹¹ m² s⁻¹. No pressure effect was detected at the conditions investigated. In experiments involving mixed halogens, the diffusivities appear to decrease slightly (by a factor of ∼3), and are more uniform among the three elements. However, activation energies for diffusion do not appear to differ between experiments with individual halogens or when they are all mixed together. The effect of water increases the diffusion coefficients of F and Cl by no more than a factor of 3 compared to the anhydrous melt (D_F = 4.0 × 10⁻¹¹ to 1.6 × 10⁻¹⁰ m² s⁻¹; D_Cl = 3.0 × 10⁻¹¹ to 1.9 × 10⁻¹⁰ m² s⁻¹). Comparing our results to the diffusion coefficients of other volatiles in nominally dry basaltic melts, halogen diffusivities are about one order of magnitude lower than H₂O, similar to CO₂, and a factor of ∼5 higher than S. The contrasting volatile diffusivities may affect the variable extent of volatile degassing upon melt depressurization and vesiculation, and can help our understanding of the compositions of rapidly grown magmatic bubbles.     

CO2-water-basalt interaction. Numerical simulation of low temperature CO2 sequestration into basalts

The interaction between CO₂-rich waters and basaltic glass was studied using reaction path modeling in order to get insight into the water-rock reaction process including secondary mineral composition, water chemistry and mass transfer as a function of CO₂ concentration and reaction progress (ξ). The calculations were carried out at 25-90 °C and pCO₂ to 30 bars and the results were compared to recent experimental observations and natural systems. A thermodynamic dataset was compiled from 25 to 300 °C in order to simulate mineral saturations relevant to basalt alteration in CO₂-rich environment including revised key aqueous species for mineral dissolution reactions and apparent Gibbs energies for clay and carbonate solid solutions observed to form in nature. The dissolution of basaltic glass in CO₂-rich waters was found to be incongruent with the overall water composition and secondary mineral formation depending on reaction progress and pH. Under mildly acid conditions in CO₂ enriched waters (pH <6.5), SiO₂ and simple Al-Si minerals, Ca-Mg-Fe smectites and Ca-Mg-Fe carbonates predominated. Iron, Al and Si were immobile whereas the Mg and Ca mobility depended on the mass of carbonate formed and water pH. Upon quantitative CO₂ mineralization, the pH increased to >8 resulting in Ca-Mg-Fe smectite, zeolites and calcite formation, reducing the mobility of most dissolved elements. The dominant factor determining the reaction path of basalt alteration and the associated element mobility was the pH of the water. In turn, the pH value was determined by the concentration of CO₂ and extent of reaction. The composition of the carbonates depended on the mobility of Ca, Mg and Fe. At pH <6.5, Fe was in the ferrous oxidation state resulting in the formation of Fe-rich carbonates with the incorporation of Ca and Mg. At pH >8, the mobility of Fe and Mg was limited due to the formation of clays whereas Ca was incorporated into calcite, zeolites and clays. Competing reactions between clays (Ca-Fe smectites) and carbonates at low pH, and zeolites and clays (Mg-Fe smectites) and carbonates at high pH, controlled the availability of Ca, Mg and Fe, playing a key role for low temperature CO₂ mineralization and sequestration into basalts. Several problems of the present model point to the need of improvement in future work. The determinant factors linking time to low temperature reaction path modeling may not only be controlled by the primary dissolving phase, which presents challenges concerning non-stoichiometric dissolution, the leached layer model and reactive surface area, but may include secondary mineral precipitation kinetics as rate limiting step for specific reactions such as retrieved from the present reaction path study.     

Sulfur isotope fractionation during the May 2003 eruption of Anatahan volcano, Mariana Islands: Implications for sulfur sources and plume processes

Sulfur isotope compositions of pumice and adsorbed volatiles on ash from the first historical eruption of Anatahan volcano (Mariana arc) are presented in order to constrain the sources of sulfur erupted during the period 10-21 May, 2003. The isotopic composition of S extracted from erupted pumice has a narrow range, from δ³⁴S_V-CDT +2.6‰ to +3.2‰, while the composition of sulfur adsorbed onto ash has a larger range (+2.8‰ to +5.3‰). Fractionation modeling for closed and open system scenarios suggests that degassing of SO₂ raised the δ³⁴S_V-CDT value of S dissolved in the melt from an initial composition of between +1.6‰ and +2.6‰ for closed-system degassing, or between −0.5‰ and +1.5‰ for open-system degassing, however closed-system degassing is the preferred model. The calculated values for the initial composition of the magma represent a MORB-like (δ³⁴S_V-CDT ∼ 0‰) mantle source with limited contamination by subducted seawater sulfate (δ³⁴S_V-CDT +21‰). Modeling also suggests that the δ³⁴S_V-CDT value of SO₂ gas in closed-system equilibrium with the degassed magma was between +0.9‰ and +2.5‰. The δ³⁴S_V-CDT value of sulfate adsorbed onto ash in the eruption plume (+2.8‰ to +5.1‰) is consistent with sulfate formation by oxidation of magmatic SO₂ in the eruption column. The sulfur isotope composition of sulfate adsorbed to ash changes from lower δ³⁴S values for ash erupted early in the eruption to higher δ³⁴S values for ash erupted later in the eruption. We interpret the temporal/stratigraphic change in sulfate isotopic composition to primarily reflect a change in the isotopic composition of magmatic SO₂ released from the progressively degassing magma and is attributed to the expulsion of an accumulated gas phase at the beginning of the eruption. More efficient oxidation of magmatic SO₂ gas to sulfate in the early water-rich eruption plume probably contributed to the change in S isotope compositions observed in the ash leachates.     

Sources, degassing, and contamination of CO2, H2O, He, Ne, and Ar in basaltic glasses from Kolbeinsey Ridge, North Atlantic

New volatile data (CO₂, H₂O, He, Ne, and Ar) are presented for 24 submarine basaltic glasses from the Kolbeinsey Ridge, Tjörnes Fracture Zone and Mohns Ridge, North Atlantic. Low CO₂ and He contents indicate that magmas were strongly outgassed with the extent of degassing increasing toward the south, as expected from shallower ridge depths. Ne and Ar are significantly more abundant in the southernmost glasses than predicted for degassed melt. The strong atmospheric isotopic signal associated with this excess Ne and Ar suggests syn- or posteruptive contamination by air. Degassing, by itself, cannot generate the large variations in δ¹³C values of dissolved CO₂ or coupled CO₂-Ar variations. This suggests that δ¹³C values were also affected by some other processes, most probably melt-crust interaction. Modelling indicates that degassing had a negligible influence on water owing to its higher solubility in basaltic melt than the other volatiles. Low H₂O contents in the glasses reflect melting of a mantle source that is not water-rich relative to the source of N-MORB.Before eruption, Kolbeinsey Ridge melts contained ∼400 ppm CO₂ with δ¹³C of −6‰, 0.1 to 0.35 wt.% H₂O, ³He/⁴He ∼11 R_A, and CO₂/³He of ∼2 × 10⁹. We model restored volatile characteristics and find homogeneous compositions in the source of Kolbeinsey Ridge magmas. Relative to the MORB-source, He and Ne are mildly fractionated while the ⁴⁰Ar/³⁶Ar may be low. The ³He/⁴He ratios in Tjörnes Fracture Zone glasses are slightly higher (13.6 R_A) than on Kolbeinsey Ridge, suggesting a greater contribution of Icelandic mantle from the south, but the lack of ³He/⁴He variation along the Kolbeinsey Ridge is inconsistent with active dispersal of Icelandic mantle beyond the Tjörnes Fracture Zone.