Carbon dioxide and argon diffusion in silicate melts: Insights into the CO2 speciation in magmas

To investigate the influence of temperature and composition on the diffusivities of dissolved carbon dioxide and argon in silicate melts, diffusion experiments were performed at magmatic pressure and temperature conditions in (a) albite melts with excess Na₂O (0-8.6 wt%) and a constant Si/Al ratio of 3, and (b) albite₇₀quartz₃₀ to jadeite melts with decreasing SiO₂ content and a constant Na/Al ratio of 1. We obtained diffusion coefficients at 500 MPa and 1323-1673 K. In the fully polymerized system Ab₇₀Qz₃₀ - Jd, the change in composition only has a weak effect on bulk CO₂ diffusivity, but Ar diffusivity increases clearly with decreasing SiO₂ content. In the system Ab + Na₂O, bulk CO₂ and Ar diffusivity increase significantly with gradual depolymerisation. The relatively small change in composition on molar basis in the depolymerized system leads to a significantly larger change in diffusivities compared to the fully polymerized Ab₇₀Qz₃₀-Jd join. Within error, activation energies for bulk CO₂ and Ar diffusion in both systems are identical with decreasing silica content (Ab + Na₂O: 159 ± 25 kJ mol⁻¹ for bulk CO₂ and 130 ± 8 kJ mol⁻¹ for Ar; Ab₇₀Qz₃₀-Jd: 163 ± 16 kJ mol⁻¹ for bulk CO₂ and 148 ± 15 kJ mol⁻¹ for Ar) even though this results in depolymerisation in one system and not the other.Although there is a variation in CO₂ speciation with changing composition as observed in quenched glasses, it has previously established that this is not a true representation of the species present in the melt, with the ratio of molecular CO₂ to carbonate decreasing during quenching. Thus, diffusion coefficients for the individual CO₂ species cannot be directly derived by measuring molecular CO₂ and CO₃^2- concentration-distance profiles in the glasses. To obtain diffusivities of individual CO₂ species, we have made two assumptions that (1) inert Ar can be used as a proxy for molecular CO₂ diffusion characteristics as shown by our previous work and (2) the diffusivity of CO₃²⁻ can be calculated assuming it is identical to network forming components (Si⁴⁺ and Al³⁺). This is derived from viscosity data (Eyring eqn.) and suggests that CO₃²⁻ diffusion would be several orders of magnitude slower than molecular CO₂ diffusion.The systematics of measured bulk CO₂ diffusivity rates and comparison with the Ar proxy all suggest that the faster molecular CO₂ species is much more dominant in melts than measurements on resulting quenched glasses would suggest. This study has confirmed an observation of surprisingly consistent bulk CO₂ diffusivity across a range of natural compositions were Ar diffusivity significantly increases. This is consistent with an actual increase in molecular CO₂ mobility (similar to Ar) that is combined with an increase in the proportion of the slower carbonate in the melt.These results demonstrate that the CO₂ diffusion and speciation model provides an insight into the transport processes in the melt and is promising and an alternative tool to in situ speciation measurements at magmatic conditions, which at the moment are technically extremely difficult. We present the first high pressure high temperature in situ MIR spectra of a CO₂ bearing albitic glass/melt suggesting that molecular CO₂ is a stable species at high temperature, which is qualitatively consistent with the modelled CO₂ speciation data.     

CO2 solubility in dacitic melts equilibrated with H2O-CO2 fluids: Implications for modeling the solubility of CO2 in silicic melts

The solubility of CO₂ in dacitic melts equilibrated with H₂O-CO₂ fluids was experimentally investigated at 1250°C and 100 to 500 MPa. CO₂ is dissolved in dacitic glasses as molecular CO₂ and carbonate. The quantification of total CO₂ in the glasses by mid-infrared (MIR) spectroscopy is difficult because the weak carbonate bands at 1430 and 1530 cm⁻¹ can not be reliably separated from background features in the spectra. Furthermore, the ratio of CO_2,mol/carbonate in the quenched glasses strongly decreases with increasing water content. Due to the difficulties in quantifying CO₂ species concentrations from the MIR spectra we have measured total CO₂ contents of dacitic glasses by secondary ion mass spectrometry (SIMS).At all pressures, the dependence of CO₂ solubility in dacitic melts on x^fluid_CO2,total shows a strong positive deviation from linearity with almost constant CO₂ solubility at x_CO2fluid > 0.8 (maximum CO₂ solubility of 795 ± 41, 1376 ± 73 and 2949 ± 166 ppm at 100, 200 and 500 MPa, respectively), indicating that dissolved water strongly enhances the solubility of CO₂. A similar nonlinear variation of CO₂ solubility with x_CO2fluid has been observed for rhyolitic melts in which carbon dioxide is incorporated exclusively as molecular CO₂ (Tamic et al., 2001). We infer that water species in the melt do not only stabilize carbonate groups as has been suggested earlier but also CO₂ molecules.A thermodynamic model describing the dependence of the CO₂ solubility in hydrous rhyolitic and dacitic melts on T, P, f_CO2 and the mol fraction of water in the melt (x_water) has been developed. An exponential variation of the equilibrium constant K₁ with x_water is proposed to account for the nonlinear dependence of x_CO2,totalmelt on x_CO2fluid. The model reproduces the CO₂ solubility data for dacitic melts within ±14% relative and the data for rhyolitic melts within 10% relative in the pressure range 100-500 MPa (except for six outliers at low x_CO2fluid). Data obtained for rhyolitic melts at 75 MPa and 850°C show a stronger deviation from the model, suggesting a change in the solubility behavior of CO₂ at low pressures (a Henrian behavior of the CO₂ solubility is observed at low pressure and low H₂O concentrations in the melt). We recommend to use our model only in the pressure range 100-500 MPa and in the x_CO2fluid range 0.1-0.95. The thermodynamic modeling indicates that the partial molar volume of total CO₂ is much lower in rhyolitic melts (31.7 cm³/mol) than in dacitic melts (46.6 cm³/mol). The dissolution enthalpy for CO₂ in hydrous rhyolitic melts was found to be negligible. This result suggests that temperature is of minor importance for CO₂ solubility in silicic melts.     

CO2 solubility and speciation in intermediate (andesitic) melts: the role of H2O and composition

We determined total CO₂ solubilities in andesite melts with a range of compositions. Melts were equilibrated with excess C-O(-H) fluid at 1 GPa and 1300°C then quenched to glasses. Samples were analyzed using an electron microprobe for major elements, ion microprobe for C-O-H volatiles, and Fourier transform infrared spectroscopy for molecular H₂O, OH⁻, molecular CO₂, and CO₃²⁻. CO₂ solubility was determined in hydrous andesite glasses and we found that H₂O content has a strong influence on C-O speciation and total CO₂ solubility. In anhydrous andesite melts with ∼60 wt.% SiO₂, total CO₂ solubility is ∼0.3 wt.% at 1300°C and 1 GPa and total CO₂ solubility increases by about 0.06 wt.% per wt.% of total H₂O. As total H₂O increases from ∼0 to ∼3.4 wt.%, molecular CO₂ decreases (from 0.07 ± 0.01 wt.% to ∼0.01 wt.%) and CO₃²⁻ increases (from 0.24 ± 0.04 wt.% to 0.57 ± 0.09 wt.%). Molecular CO₂ increases as the calculated mole fraction of CO₂ in the fluid increases, showing Henrian behavior. In contrast, CO₃²⁻ decreases as the calculated mole fraction of CO₂ in the fluid increases, indicating that CO₃²⁻ solubility is strongly dependent on the availability of reactive oxygens in the melt. These findings have implications for CO₂ degassing. If substantial H₂O is present, total CO₂ solubility is higher and CO₂ will degas at relatively shallow levels compared to a drier melt. Total CO₂ solubility was also examined in andesitic glasses with additional Ca, K, or Mg and low H₂O contents (<1 wt.%). We found that total CO₂ solubility is negatively correlated with (Si + Al) cation mole fraction and positively correlated with cations with large Gibbs free energy of decarbonation or high charge-to-radius ratios (e.g., Ca). Combining our andesite data with data from the literature, we find that molecular CO₂ is more abundant in highly polymerized melts with high ionic porosities (>∼48.3%), and low nonbridging oxygen/tetrahedral oxygen (<∼0.3). Carbonate dominates most silicate melts and is most abundant in depolymerized melts with low ionic porosities, high nonbridging oxygen/tetrahedral oxygen (>∼0.3), and abundant cations with large Gibbs free energy of decarbonation or high charge-to-radius ratio. In natural silicate melt, the oxygens in the carbonate are likely associated with tetrahedral and network-modifying cations (including Ca, H, or H-bonds) or a combinations of those cations.     

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.     

Energetics of multicomponent diffusion in molten CaO-Al2O3-SiO2

The energetics of multicomponent diffusion in molten CaO-Al₂O₃-SiO₂ (CAS) were examined experimentally at 1440 to 1650°C and 0.5 to 2 GPa. Two melt compositions were investigated: a haplodacitic melt (25 wt.% CaO, 15% Al₂O₃, and 60% SiO₂) and a haplobasaltic melt (35% CaO, 20% Al₂O₃, and 45% SiO₂). Diffusion matrices were measured in a mass-fixed frame of reference with simple oxides as end-member components and Al₂O₃ as a dependent variable. Chemical diffusion in molten CAS shows clear evidence of diffusive coupling among the components. The diffusive flux of SiO₂ is significantly enhanced whenever there is a large CaO gradient that is oriented in a direction opposite to the SiO₂ gradient. This coupling effect is more pronounced in the haplodacitic melt and is likely to be significant in natural magmas of rhyolitic to andesitic compositions. The relative magnitude of coupled chemical diffusion is not very sensitive to changes in temperature and pressure.To a good approximation, the measured diffusion matrices follow well-defined Arrhenius relationships with pressure and reciprocal temperature. Typically, a change in temperature of 100°C results in a relative change in the elements of diffusion matrix of 50 to 100%, whereas a change in pressure of 1 GPa introduces a relative change in elements of diffusion matrix of 4 to 6% for the haplobasalt, and less than 5% for the haplodacite. At a pressure of 1 GPa, the ratios between the major and minor eigenvalues of the diffusion matrix λ₁/λ₂ are not very sensitive to temperature variations, with an average of 5.5 ± 0.2 for the haplobasalt and 3.7 ± 0.6 for the haplodacite. The activation energies for the major and minor eigenvalues of the diffusion matrix are 215 ± 12 and 240 ± 21 kJ mol⁻¹, respectively, for the haplodacite and 192 ± 8 and 217 ± 14 kJ mol⁻¹ for the haplobasalt. These values are comparable to the activation energies for self-diffusion of calcium and silicon at the same melt compositions and pressure. At a fixed temperature of 1500°C, the ratios λ₁/λ₂ increase with the increase of pressure, with λ₁/λ₂ varying from 2.5 to 4.1 (0.5 to 1.3 GPa) for the haplodacite and 4 to 6.5 (0.5 to 2.0 GPa) for the haplobasalt. The activation volumes for the major and minor eigenvalues of the diffusion matrix are 0.31 ± 0.44 and 2.3 ± 0.8 cm³ mol⁻¹, respectively, for the haplodacite and −1.48 ± 0.18 and −0.42 ± 0.24 cm³ mol⁻¹ for the haplobasalt. These values are quite different from the activation volumes for self-diffusion of calcium and silicon at the same melt compositions and temperature. These differences in activation volumes between the two melts likely result from a difference in the structure and thermodynamic properties of the melt between the two compositions (e.g., partial molar volume).Applications of the measured diffusion matrices to quartz crystal dissolution in molten CAS reveal that the activation energy and activation volume for quartz dissolution are almost identical to the activation energy and activation volume for diffusion of the minor or slower eigencomponent of the diffusion matrix. This suggests that the diffusion rate of slow eigencomponent is the rate-limiting factor in isothermal crystal dissolution, a conclusion that is likely to be valid for crystal growth and dissolution in natural magmas when diffusion in liquid is the rate-limiting factor.     

Behavior of molybdenum in alkali silicate melts: Effects of excess SiO2 and CO2

The effects of excess SiO₂ and CO₂ on the solubility of molybdenite in hydrous sodium disilicate melts were experimentally determined at 680 bars and 650°C. The molybdenite solubility decreases with increasing SiO₂ and CO₂. Under the experimental conditions, the MoS₂ content of the vapor-saturated liquid decreases from 10 wt.% to 2.5 wt.% at SiO₂ saturation. In the presence of CO₂, the solubility decreases to 4.6 wt.% MoS₂ and becomes negligible at high P_CO2. These results are explained as deriving from the increased polymerization and hence decreased NBO/Si ratio of the melt with increasing SiO₂ content and CO₂, respectively. Sulfur dissolves principally as SO₄^?2 at the relatively high fo₂ of the experiments. Consequently, the effect of sulfur is to lower the Mo solubility by effectively decreasing the NBO/Si ratio of the melt. Sulfur saturation is, therefore, likely to be a limiting factor in the Mo content of alkali silicate melts because of the chalcophile affinities of molybdenum.     

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.     

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.     

Oxygen isotope exchange between H2O and CO2 at elevated CO2 pressures: Implications for monitoring of geological CO2 storage

Traditionally, the application of stable isotopes in Carbon Capture and Storage (CCS) projects has focused on δ¹³C values of CO₂ to trace the migration of injected CO₂ in the subsurface. More recently the use of δ¹⁸O values of both CO₂ and reservoir fluids has been proposed as a method for quantifying in situ CO₂ reservoir saturations due to O isotope exchange between CO₂ and H₂O and subsequent changes in δ¹⁸O_H2O values in the presence of high concentrations of CO₂. To verify that O isotope exchange between CO₂ and H₂O reaches equilibrium within days, and that δ¹⁸O_H2O values indeed change predictably due to the presence of CO₂, a laboratory study was conducted during which the isotope composition of H₂O, CO₂, and dissolved inorganic C (DIC) was determined at representative reservoir conditions (50 °C and up to 19 MPa) and varying CO₂ pressures. Conditions typical for the Pembina Cardium CO₂ Monitoring Pilot in Alberta (Canada) were chosen for the experiments. Results obtained showed that δ¹⁸O values of CO₂ were on average 36.4 ± 2.2‰ (1σ, n = 15) higher than those of water at all pressures up to and including reservoir pressure (19 MPa), in excellent agreement with the theoretically predicted isotope enrichment factor of 35.5‰ for the experimental temperatures of 50 °C. By using ¹⁸O enriched water for the experiments it was demonstrated that changes in the δ¹⁸O values of water were predictably related to the fraction of O in the system sourced from CO₂ in excellent agreement with theoretical predictions. Since the fraction of O sourced from CO₂ is related to the total volumetric saturation of CO₂ and water as a fraction of the total volume of the system, it is concluded that changes in δ¹⁸O values of reservoir fluids can be used to calculate reservoir saturations of CO₂ in CCS settings given that the δ¹⁸O values of CO₂ and water are sufficiently distinct.     

CO2 in haplo-phonolite melt: solubility, speciation and carbonate complexation

CO₂ solubility was measured in a synthetic iron-free phonolite (haplo-phonolite) by equilibrating melt with excess CO₂ fluid in a piston cylinder apparatus for a range of pressures (1.0- 2.5 GPa) and temperatures (1300 to 1550°C). The quenched glasses were then analysed using a bulk carbon analytical method (LECO). The measured solubilities are between 0.65 and 2.77 wt.% for the range of conditions studied and show a negative correlation with temperature as reported for most other silicate melt compositions.A range of carbonate species are present within the glass, as well as minor amounts of molecular CO₂. FTIR and NMR analyses suggest that carbonate is present as both ‘network’ and ‘depolymerised’ units as shown for relatively highly polymerised compositions in the model of Brooker et al. (2001b). The bulk CO₂ analyses were used to calibrate the IR extinction coefficient for the carbonate groups. However, the results show that the values obtained for the glasses vary with the melt equilibration conditions, presumably because the ratio of the different carbonate species changes as a complex function of run pressure, temperature and quench rate. Thus the use of IR may not be a reliable method for the quantification of dissolved CO₂ concentrations in natural glasses of ‘intermediate’ composition.     

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.     

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.     

Kimberlite petrogenesis: Insights from clinopyroxene-melt partitioning experiments at 6 GPa in the CaO-MgO-Al2O3-SiO2-CO2 system

In this experimental study, we examine the mineral-melt partitioning of major and trace elements between clinopyroxene and CO₂-rich kimberlitic melts at a pressure of 6 GPa and temperatures of 1410°C and 1430°C. The melts produced contain ∼ 28 wt% dissolved CO₂, and are saturated with olivine and clinopyroxene. To assess the effects of temperature, crystal and melt compositions on trace element partitioning, experiments were performed in the model CaO-MgO-Al₂O₃-SiO₂-CO₂ system. Our results reveal that all the elements studied, except Al, Mg, Si, and Ga, are incompatible in clinopyroxene. Partition coefficients show a considerable range in magnitude, from ∼ 10⁻³ for D_U and D_Ba to ∼ 2.5 for D_Si. The two experimental runs show similar overall partitioning patterns with the D values being lower at 1430°C. Rare earth elements display a wide range of partition coefficients, D_La (0.012-0.026) being approximately one order of magnitude lower than D_Lu (0.18-0.23). Partition coefficients for the 2+ and 3+ cations entering the M2-site exhibit a near-parabolic dependence on radius of the incorporated cations as predicted from the lattice strain model. This underlines the contribution made by the crystal structure toward controlling the distribution of trace elements. Using data obtained in this study combined with that in the published literature, we also discuss the effects that other important parameters, namely, melt composition, pressure, and temperature, could have on partitioning.Our partition coefficients have been used to model the generation of the Group I (GI) kimberlites from South Africa. The numerical modeling shows that kimberlitic melts can be produced by ∼0.5% melting of a MORB-type depleted source that has been enriched by small-degree melts originating from a similar depleted source. This result suggests that the source of GI kimberlites may be located at the lithosphere-asthenosphere transition. Percolation of small degree melts from the asthenosphere would essentially create a metasomatic horizon near the bottom of the non-convecting sublithospheric mantle. Accumulation of such small degree melts together with the presence of volatiles and conductive heating would trigger melting of the ambient mantle and subsequently lead to eruption of kimberlitic melts. Additionally, our model shows that the GI source can be generated by metasomatism of a 2 Ga old MORB source ca. 1 Ga ago. Assuming that MORB-type mantle is the most depleted source of magmas on earth, then this is the oldest age at which the GI source could have existed. However, this age most likely reflects the average age of a series of metasomatic events than that of a single event.     

A new view of the He-Ar-CO2 degassing at mid-ocean ridges: Homogeneous composition of magmas from the upper mantle

Deep-sea exploration is rapidly improving our understanding of volatiles geochemistry in mid-ocean-ridge igneous products. It is also placing greater constraints on degassing processes of the Earth’s mantle, with the result that degassing models based on vapour-melt equilibrium are no longer able to explain the increasing number of data. In fact, such models force to postulate an upper mantle strongly heterogeneous at any scale, and cannot account for the widespread carbon supersaturation of the recovered igneous products. Here we review the global He-Ar-CO₂ dataset of fluid inclusions in mid-ocean-ridge glasses using the framework of advanced modelling of multicomponent bubble growth in magmas. We display that non-equilibrium fractionations among He, Ar and CO₂, driven by their different diffusivities in silicate melts, are common in most of the natural conditions of magma decompression and their signature strongly depends on pressure of degassing. Due to the comparable Ar and CO₂ diffusivity, magma degassing at low pressure fractionates both the He/Ar and He/CO₂ ratio by a similar extent, while the slower CO₂ diffusion at high pressure causes early kinetic effects on Ar/CO₂ ratio and dramatically changes the degassing path. On this ground, the very different geochemical signatures among suites of data coming from different ridge segments mainly depend on the depth of the magma chamber where the melt was stored. Besides, the variations inside a single suite highlight variable ascent speed and cooling rate of the emplaced lava. The large variations in both the He/CO₂ and Ar/CO₂ ratios at almost constant He/Ar, displayed in glasses coming from the Mid-Atlantic Ridge 24-30°N segment and the Rodriguez Triple Junction, are therefore interpreted as a high-pressure signature. In contrast, the simultaneous increase in both He/CO₂ and He/Ar of the East Pacific Rise, Pito Seamount and South-East Indian Ridge data sets suggests the dominance of low-pressure fractionation, implying that the shallow magma chambers are at a lower depth than those of the Mid-Atlantic Ridge 24-30°N and Rodriguez Triple Junction. Our conclusions support the presence of a relationship between spreading rate and depth of high-temperature zones below ridges, and are consistent with the depth of magma chambers as suggested from seismic studies. Non-equilibrium degassing explains the volatile systematics of mid-ocean-ridge basalts by starting from a single mantle-derived magma, dispensing with the supposed need for heterogeneities in abundance ratios of volatiles in the mantle below oceanic ridges.     

A model for O/O variations in CO2 evolved from goethite during the solid-state α-FeOOH to α-Fe2O3 phase transition

“Plateau” δ¹⁸O values of CO₂ that evolved from the Fe(CO₃)OH component during isothermal vacuum dehydrations (200-230 °C) of 18 natural goethites range from 8.2 to 28.1‰. In contrast, the measured δ¹⁸O values of the goethite structural oxygen range from −11.3 to 1.7‰. The results of this study indicate that the apparent oxygen isotope fractionation factor (¹⁸α_app) between plateau CO₂ and initial goethite is systematically related to the rate of isothermal vacuum dehydration. The nonlinear correlation and the magnitudes of the ¹⁸α_app values are predicted by a relatively simple mass balance model with the following assumptions: (1) the rate of isothermal vacuum dehydration of goethite (for the interval from 0 to ∼60 to 80% loss of structural hydroxyl hydrogen) can be reasonably well represented by first-order kinetics and (2) isotopic exchange between evolving H₂O vapor and solid occurs only in successive, local transition states. The generally good correspondence between the model predictions and the experimental data seems to validate these assumptions. Thus, the ¹⁸O/¹⁶O ratios of the evolved CO₂ can act as probes into the transient processes operating at the molecular level during the solid-state goethite-to-hematite phase transition. For example, the activation energy for the rate constant associated with the transition state, oxygen isotopic exchange between solid and H₂O vapor, is tentatively estimated as 28 ± 11 KJ/mol. Such knowledge may be of consequence in understanding the significance of ¹⁸O/¹⁶O ratios in hematites from some natural environments (e.g., Mars?).Kinetic data and δ¹⁸O values of CO₂ are routinely obtained in the course of measurements of the abundance and δ¹³C values of the Fe(CO₃)OH in goethite. The observed correlation between ¹⁸α_app and dehydration rates suggests that plateau δ¹⁸O values of evolved CO₂ may provide complementary estimates of the δ¹⁸O values of total goethite structural oxygen (O, OH, CO₂) with an overall precision of about ±1‰. However, because of isotopic exchange during the dehydration process, δ¹⁸O values of the evolved CO₂ do not reflect the original δ¹⁸O values of the CO₂ that was occluded as Fe(CO₃)OH in goethite.     

H2O diffusion in dacitic and andesitic melts

The diffusion of water in dacitic and andesitic melts was investigated at temperatures of 1458 to 1858 K and pressures between 0.5 and 1.5 GPa using the diffusion couple technique. Pairs of nominally dry glasses and hydrous glasses containing between 1.5 and 6.3 wt.% dissolved H₂O were heated for 60 to 480 s in a piston cylinder apparatus. Concentration profiles of hydrous species (OH groups and H₂O molecules) and total water (C_H2Ot = sum of OH and H₂O) were measured along the cylindrical axis of the diffusion sample using IR microspectroscopy. Electron microprobe traverses show no significant change in relative proportions of anhydrous components along H₂O profiles, indicating that our data can be treated as effective binary interdiffusion between H₂O and the rest of the silicate melt. Bulk water diffusivity (D_H2Ot) was derived from profiles of total water using a modified Boltzmann-Matano method as well as using fittings assuming a functional relationship between D_H2Ot and C_H2Ot. In dacitic melts D_H2Ot is proportional to C_H2Ot up to 6 wt.%. In andesitic melts the dependence of D_H2Ot on C_H2Ot is less pronounced. A pressure effect on water diffusivity could not be resolved for either dacitic or andesitic melt in the range 0.5 to 1.5 GPa. Combining our results with previous studies on water diffusion in rhyolite and basalt show that for a given water content D_H2Ot increases monotonically with increasing melt depolymerization at temperatures >1500 K. Assuming an Arrhenian behavior in the whole compositional range, the following formulation was derived to estimate D_H2Ot (m²/s) at 1 wt.% H₂O_t in melts with rhyolitic to andesitic composition as a function of T (K), P (MPa) and S (wt.% SiO₂):

     

Caprock interaction with CO2: A laboratory study of reactivity of shale with supercritical CO2 and brine

Crushed rock from two caprock samples, a carbonate-rich shale and a clay-rich shale, were reacted with a mixture of brine and supercritical CO₂ (CO₂–brine) in a laboratory batch reactor, at different temperature and pressure conditions. The samples were cored from a proposed underground CO₂ storage site near the town of Longyearbyen in Svalbard. The reacting fluid was a mixture of 1 M NaCl solution and CO₂ (110 bar) and the water/rock ratio was 20:1. Carbon dioxide was injected into the reactors after the solution had been bubbled with N₂, in order to mimic O₂-depleted natural storage conditions. A control reaction was also run on the clay-rich shale sample, where the crushed rock was reacted with brine (CO₂-free brine) at the same experimental conditions. A total of 8 batch reaction experiments were run at temperatures ranging from 80 to 250 °C and total pressures of 110 bar (∼40 bar for the control experiment). The experiments lasted 1–5 weeks.Fluid analysis showed that the aqueous concentration of major elements (i.e. Ca, Mg, Fe, K, Al) and SiO₂ increased in all experiments. Release rates of Fe and SiO₂ were more pronounced in solutions reacted with CO₂–brine as compared to those reacted with CO₂-free brine. For samples reacted with the CO₂–brine, lower temperature reactions (80 °C) released much more Fe and SiO₂ than higher temperature reactions (150–250 °C). Analysis by SEM and XRD of reacted solids also revealed changes in mineralogical compositions. The carbonate-rich shale was more reactive at 250 °C, as revealed by the dissolution of plagioclase and clay minerals (illite and chlorite), dissolution and re-precipitation of carbonates, and the formation of smectite. Carbon dioxide was also permanently sequestered as calcite in the same sample. The clay-rich shale reacted with CO₂–brine did not show major mineralogical alteration. However, a significant amount of analcime was formed in the clay-rich shale reacted with CO₂-free brine; while no trace of analcime was observed in either of the samples reacted with CO₂–brine.     

Na2O solubility in CaO-MgO-SiO2 melts

The sodium solubility in silicate melts in the CaO-MgO-SiO₂ (CMS) system at 1400 °C has been measured by using a closed thermochemical reactor designed to control alkali metal activity. In this reactor, Na_(g) evaporation from a Na₂O-xSiO₂ melt imposes an alkali metal vapor pressure in equilibrium with the molten silicate samples. Because of equilibrium conditions in the reactor, the activity of sodium-metal oxide in the molten samples is the same as that of the source, i.e., aNa₂O_(sample) = aNa₂O_(source). This design also allows to determine the sodium oxide activity coefficient in the samples. Thirty-three different CMS compositions were studied. The results show that the amount of sodium entering from the gas phase (i.e., Na₂O solubility) is strongly sensitive to silica content of the melt and, to a lesser extent, the relative amounts of CaO and MgO. Despite the large range of tested melt compositions (0 < CaO and MgO < 40; 40 < SiO₂ < 100; in wt%), we found that Na₂O solubility is conveniently modeled as a linear function of the optical basicity (Λ) calculated on a Na-free basis melt composition. In our experiments, γNa₂O_(sample) ranges from 7 × 10⁻⁷ to 5 × 10⁻⁶, indicating a strongly non-ideal behavior of Na₂O solubility in the studied CMS melts (γNa₂O_(sample) ? 1). In addition to showing the effect of sodium on phase relationships in the CMS system, this Na₂O solubility study brings valuable new constraints on how melt structure controls the solubility of Na in the CMS silicate melts. Our results suggest that Na₂O addition causes depolymerization of the melt by preferential breaking of Si-O-Si bonds of the most polymerized tetrahedral sites, mainly Q⁴.     

The CO2-He-Ar-H2O systematics of the manus back-arc basin: resolving source composition from degassing and contamination effects

We report new stepped heating He, Ar, CO₂ and water data on a petrogenetically diverse suite of lavas from the Manus back-arc basin, where a plume component has previously been identified. The aim of this study is to evaluate the superimposed effects of degassing and contamination in order to identify mantle source characteristics. CO₂ abundances and carbon isotopes in both the vesicle ([CO₂] up to 180ppm; δ¹³C as low as -33.6 ‰) and glass ([CO₂] up to 270ppm; δ¹³C as low as -34.3 ‰) phases reveal that samples have been modified by varying degrees of degassing. High water concentration samples (back-arc basin basalts (BABB) and arc type samples) show the highest degrees of degassing (i.e. lower δ¹³C values and lower CO₂ contents). The results are modelled for both the glass and vesicle phases using batch and fractional degassing models. Parental melt compositions can be constrained to show the following CO₂ concentration trend: arc-type > BABB s.r. (southern rift) > MORB-2, E-MORB, X-BABB (extreme BABB), BABB > MORB-1 and MORB-smt. ⁴He/⁴⁰Ar∗ ratios of samples (14.6-1100) are consistent with residual volatiles from a degassed source. Variations in CO₂/³He values are likely due to degassing, followed by contamination from a crustal source (either the subducting Solomon Sea Plate or the pre-existing crust through which the lavas erupt), as evidenced by high K₂O/TiO₂ ratios and low δ¹³C. The CO₂/³He of the Manus plume is best estimated by the MORB-smt and MORB-1 samples at 3.1 ± 0.6 x 10⁹. This value is similar to previous estimates of plume CO₂/³He values, which are either equal to or slightly greater than the upper mantle average of 2 x 10⁹.     

The effect of CO2 on the speciation of RbBr in solution at temperatures to 579 °C and pressures to 0.26 GPa

Carbon dioxide- and salt-bearing solutions are common in granulite, ore-forming and magmatic environments. The presence of CO₂ affects mineral solubilities, fluid miscibility, and viscosity and wetting properties, and is expected to affect salt speciation. EXAFS measurements of RbBr-H₂O-CO₂ fluids contained in corundum-osed synthetic fluid inclusions (SFLINCs) have been used to investigate the effect of CO₂ on salt speciation at temperatures to 579 °C and pressures to around 0.26 GPa.Forward modelling indicates that solute dehydration is difficult to distinguish from up to around 40% of Rb-Br ion-pairing, so results refer to the total number of nearest neighbours, which are likely to be mostly O present in waters of hydration, but may also include Br, if ion pairing is present. Additionally, results relate to the number of well-ordered neighbours in the first shell, because nearest neighbours with a high degree of disorder may be present but contribute minimally to the EXAFS signal. Analysis of the EXAFS results at the Rb edge for the CO₂-free solution is consistent with previous work and shows that the number of nearest neighbours for Rb in CO₂-free solutions decreases from 6 ± 0.6 to 1.4 ± 0.1 as temperature increases from 20 to 534 °C. The decrease is accompanied by a decrease in Rb-x bondlengths of 0.05 Å, where x is the first shell scatterer. Results for the CO₂-bearing solution are different to those for the CO₂-free solution. The number of nearest neighbours is 16 and 22% less than for the CO₂-bearing solution at 312 and 445 °C respectively. Changes in the numbers of nearest neighbours correlate well with calculated changes in the bulk solution dielectric constant; CO₂-bearing and CO₂-free solutions lie on the same trend, which suggests that it may be possible to calculate the number of nearest neighbours from dielectric constant. Rb-x bondlengths for the CO₂-bearing solution are statistically indistinguishable to those for the CO₂-free inclusions. Results for Br are worse quality than for Rb so EXAFS analysis could not be completed, however XANES spectra for CO₂-free and CO₂-bearing solutions are consistent with solute dehydration similar to that recorded by the Rb spectra. The conclusions of this study provide support for the notion that CO₂ has a fundamental effect on the mechanics of solubility, and that these effects should be incorporated into conceptual and quantitative thermodynamic models.