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.     

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⁹.     

Concentrations,sources, and losses of H2O,CO2, and S in Kilauean basalt

Basaltic glasses included in olivine phenocrysts from Kilauea volcano contain concentrations of H₂O, CO₂, and S similar to glassy Kilauean basalt dredged from the deep sea floor and greater than vesicular, subaerial Kilauean basalt. Our result contrasts with earlier reports that inclusions of basaltic glass in phenocrysts have little or no H₂O and large ratios of $\text{CO}{}_2\text{H}{}_2\text{O}$. Our analysed inclusions of glass are larger than 100 micrometers thick and similar in chemical composition to the host glass surrounding the olivine crystals indicating that the trapped melts are representative of the bulk liquid from which the crystals grew. Crystallization of about 2–8% of olivine from the melts after they were trapped is indicated by slight departures from the experimentally established equilibrium distribution of Mg and Fe between olivine and liquid. The measured concentrations of CO₂ correspond to phenocryst crystallization pressures of about 1.3 kbar for a subaerial basalt and about 5 kbar for a submarine basalt, consistent with geophysical models of Kilauea volcano. The compositions of volcanic gas predicted from our analyses are consistent with restored compositions of actual Kilauean gases. The rate of sulfur emission predicted from our analyses is greater than the sulfur dioxide emission rate observed during repose, but probably consistent with total degassing including eruptive episodes. The concentrations of H₂O, K₂O, Cl, and P in parental Kilauean basalt can be derived from upper mantle phlogopitic mica, pargasitic amphibole and apatite with compositions close to those of natural primary minerals in ultramafic xenoliths from continental kimberlites, or solely from apatite and phlogopitic mica with $\text{H}{}_2\text{O}\text{K}{}_2\text{O}$ near 0.47 ± 0.03, slightly higher than the range of values reported. The amounts of phlogopitic mica and pargasitic amphibole contributing volatiles to Kilauean tholeiite is about 10 percent by mass of the parental liquid, or about 5% if the source does not include amphibole. In view of an estimated 20% of partial melting of mantle source rock to produce Kilauean tholeiites, there may be about 2 weight percent of mica plus amphibole in part of the mantle beneath Kilauea, or about 1 weight percent of phlogopitic mica if amphibole is absent.     

Platinum-group element systematics in Mid-Oceanic Ridge basaltic glasses from the Pacific, Atlantic, and Indian Oceans

The concentrations of Ir, Ru, Pt and Pd have been determined in 29 Mid-Oceanic Ridge basaltic (MORB) glasses from the Pacific (N = 7), the Atlantic (N = 10) and the Indian (N = 11) oceanic ridges and the Red Sea (N = 1) spreading centers. The effect of sulfide segregation during magmatic differentiation has been discussed with sample suites deriving from parental melts produced by high (16%) and low (6%) degrees of partial melting, respectively. Both sample suites define positive and distinct covariation trends in platinum-group elements (PGE) vs. Ni binary plots. The high-degree melting suite displays, for a given Ni content, systematically higher PGE contents relative to the low-degree melting suite. The mass fraction of sulfide segregated during crystallization (X_sulf), the achievement of equilibrium between sulfide melt and silicate melts (R_eff), and the respective proportions between fractional and batch crystallization processes (S_b) are key parameters for modeling the PGE partitioning behavior during S-saturated MORB differentiation. Regardless of the model chosen, similar sulfide melt/silicate melt partition coefficients for Ir, Ru, Pt and Pd are needed to model the sulfide segregation process, in agreement with experimental data. When corrected for the effect of magmatic differentiation, the PGE data display coherent variations with partial melting degrees. Iridium, Ru and Pt are found to be compatible in nonsulfide minerals whereas the Pd behaves as a purely chalcophile element. The calculated partition coefficients between mantle sulfides and silicate melts (assuming a PGE concentration in the oceanic mantle at ∼0.007 × CI-chondritic abundances) increase from Pd (∼10³) to Ir (∼10⁵). This contrasting behavior of PGE during S-saturated magmatic differentiation and mantle melting processes can be accounted for by assuming that Monosufide Solid Solution (Mss) controls the PGE budget in MORB melting residues whereas MORB differentiation processes involve Cu-Ni-rich sulfide melt segregation.     

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.     

Hydrogen isotope exchange kinetics between H2O and H4SiO4 from ab initio calculations

Hydrogen isotope exchange between water and orthosilicic acid (H₄SiO₄) was modeled using B3LYP calculations and classical transition-state theory. Configurations of 1, 2, 3 and 7 water molecules and H₄SiO₄ were used to investigate energetically viable reaction pathways. An upper-bound of 71 kJ/mol was assumed for the zero-point energy corrected barrier (ZPECB) because this is the experimentally determined activation energy for Si-O bond breaking (Rimstidt and Barnes, 1980) and ZPECB is expected to be close to this value. Long range solvation forces were accounted for using the integral equation formalism polarized continuum model (IEFPCM; Cancès et al., 1997). Primary and secondary isotope effects were computed by exchanging hydrogen atoms with deuterium. Results show that reaction mechanisms involving 3 and 7 water molecules have ZPECB of 34 to 38 kJ/mol, whereas those involving 1 and 2 water molecules have ZPECB in excess of the set upper-bound. The lower range of ZPECB with 3 or 7 water molecules is reasonable to explain rapid hydrogen isotope exchange with silicates. Rate constant calculations accounting for tunneling, anharmonicity and scaling factors indicate that the reaction is fast and equilibrium can be assumed under most geologic conditions.     

顶空气中氦氖氢气相色谱测试方法及在氦资源调查中的应用

渭河盆地具有良好的氦资源前景,在勘探调查中氦气含量的测定至关重要。笔者针对渭河盆地顶空气样品中的氦氖氢测试方法与应用进行研究探讨。本方法通过气密性进样装置解决进样过程中气体样品中含较多水汽与气体逸散等问题,提高结果准确性与稳定性。采用阀切换反吹技术有效避免无关组分干扰并保护色谱柱。经多次实验优化确定仪器分析条件：载气流速14 mL/min;检测器（TCD）温度：180℃;柱箱温度：40℃;阀切换时间：0.65 min。通过检出限、正确度和精密度等参数确定了方法的可行性。最后将该方法应用于渭河盆地某剖面的顶空气样品氦气的测试中,结果满意。证明该方法具有良好的应用前景。     

He,Ne, Ar stepwise crushing data on basalt glasses from different segments of Bouvet Triple Junction

Here we present the first data on He, Ne, Ar isotopic and elemental composition in fluid phases of tholeiitic chilled glasses from the Bouvet Triple Junction (BTJ). The chilled glasses from several dredging stations situated at different segments of BTJ have been investigated: Spiess Ridge, Mid Atlantic Ridge (MAR) and in a valley of the Southwest Indian Ridge (SWIR). The data allow to distinguish within BTJ three segments characterized by different geochemical behavior of He, Ne and Ar. MAR and Spiess samples contain MORB-like helium and neon while SWIR is characterized by addition of plume type He and Ne. The strong atmospheric contamination is typical of all segments, but for MAR it is less pronounced. The Ne-Ar isotope systematics suggests that the atmospheric component was most probably introduced into the mantle source of the fluids with fragments of oceanic crust/sediments.     

Equation of state of the H2O, CO2, and H2O-CO2 systems up to 10 GPa and 2573.15 K: Molecular dynamics simulations with ab initio potential surface

Based on our previous development of the molecular interaction potential for pure H₂O and CO₂ [Zhang, Z.G., Duan, Z.H. 2005a. Isothermal-isobaric molecular dynamics simulations of the PVT properties of water over wide range of temperatures and pressures. Phys. Earth Planet Interiors149, 335-354; Zhang, Z.G., Duan, Z.H. 2005b. An optimized molecular potential for carbon dioxide. J. Chem. Phys.122, 214507] and the ab initio potential surface across CO₂-H₂O molecules constructed in this study, we carried out more than one thousand molecular dynamics simulations of the PVTx properties of the CO₂-H₂O mixtures in the temperature-pressure range from 673.15 to 2573.15 K up to 10.0 GPa. Comparison with extensive experimental PVTx data indicates that the simulated results generally agree with experimental data within 2% in density, equivalent to experimental uncertainty. Even the data under the highest experimental temperature-pressure conditions (up to 1673 K and 1.94 GPa) are well predicted with the agreement within 1.0% in density, indicating that the high accuracy of the simulation is well retained as the temperature and pressure increase. The consistent and stable predictability of the simulation from low to high temperature-pressure and the fact that the molecular dynamics simulation resort to no experimental data but to ab initio molecular potential makes us convinced that the simulation results should be reliable up to at least 2573 K and 10 GPa with errors less than 2% in density. In order to integrate all the simulation results of this study and previous studies [Zhang and Duan, 2005a, 2005b] and the experimental data for the calculation of volumetric properties (volume, density, and excess volume), heat properties, and chemical properties (fugacity, activity, and possibly supercritical phase separation), an equation of state (EOS) is laboriously developed for the CO₂, H₂O, and CO₂-H₂O systems. This EOS reproduces all the experimental and simulated data covering a wide temperature and pressure range from 673.15 to 2573.15 K and from 0 to 10.0 GPa within experimental or simulation uncertainty.     

Molecular H2O in armenite, BaCa2Al6Si9O30·2H2O, and epididymite, Na2Be2Si6O15·H2O: Heat capacity, entropy and local-bonding behavior of confined H2O in microporous silicates

Armenite, ideal formula BaCa₂Al₆Si₉O₃₀·2H₂O, and its dehydrated analog BaCa₂Al₆Si₉O₃₀ and epididymite, ideal formula Na₂Be₂Si₆O₁₅·H₂O, and its dehydrated analog Na₂Be₂Si₆O₁₅ were studied by low-temperature relaxation calorimetry between 5 and 300 K to determine the heat capacity, C_p, behavior of their confined H₂O. Differential thermal analysis and thermogravimetry measurements, FTIR spectroscopy, electron microprobe analysis and powder Rietveld refinements were undertaken to characterize the phases and the local environment around the H₂O molecule.The determined structural formula for armenite is Ba_0.88(0.01)Ca_1.99(0.02)Na_0.04(0.01)Al_5.89(0.03)Si_9.12(0.02)O₃₀·2H₂O and for epididymite Na_1.88(0.03)K_0.05(0.004)Na_0.01(0.004)Be_2.02(0.008)Si_6.00(0.01)O₁₅·H₂O. The infrared (IR) spectra give information on the nature of the H₂O molecules in the natural phases via their H₂O stretching and bending vibrations, which in the case of epididymite only could be assigned. The powder X-ray diffraction data show that armenite and its dehydrated analog have similar structures, whereas in the case of epididymite there are structural differences between the natural and dehydrated phases. This is also reflected in the lattice IR mode behavior, as observed for the natural phases and the H₂O-free phases. The standard entropy at 298 K for armenite is S° = 795.7 ± 6.2 J/mol K and its dehydrated analog is S° = 737.0 ± 6.2 J/mol K. For epididymite S° = 425.7 ± 4.1 J/mol K was obtained and its dehydrated analog has S° = 372.5 ± 5.0 J/mol K. The heat capacity and entropy of dehydration at 298 K are Δ = 3.4 J/mol K and ΔS^rxn = 319.1 J/mol K and Δ = −14.3 J/mol K and ΔS^rxn = 135.7 J/mol K for armenite and epididymite, respectively. The H₂O molecules in both phases appear to be ordered. They are held in place via an ion-dipole interaction between the H₂O molecule and a Ca cation in the case of armenite and a Na cation in epididymite and through hydrogen-bonding between the H₂O molecule and oxygen atoms of the respective silicate frameworks. Of the three different H₂O phases ice, liquid water and steam, the C_p behavior of confined H₂O in both armenite and epididymite is most similar to that of ice, but there are differences between the two silicates and from the C_p behavior of ice. Hydrogen-bonding behavior and its relation to the entropy of confined H₂O at 298 K is analyzed for various microporous silicates.The entropy of confined H₂O at 298 K in various silicates increases approximately linearly with increasing average wavenumber of the OH-stretching vibrations. The interpretation is that decreased hydrogen-bonding strength between a H₂O molecule and the silicate framework, as well as weak ion-dipole interactions, results in increased entropy of H₂O. This results in increased amplitudes of external H₂O vibrations, especially translations of the molecule, and they contribute strongly to the entropy of confined H₂O at T < 298 K.     

Review of the systematics of CO2–H2O fluid inclusions

Aqueous solutions that contain volatile (gas) components are one of the most important types of fluid in the Earth's crust. The record that such fluids have left in the form of fluid inclusions in minerals provides a wealth of insight into the geochemical and petrologic processes in which the fluids participated. This article reviews the systematics of CO₂–H₂O fluid inclusions as a starting point for interpreting the chemically more complex systems. The phase relations of the binary are described with respect to a qualitative P–T–X model, and isoplethic–isochoric paths through this model are used to explain the equilibrium and non-equilibrium behaviour of fluid inclusions during microthermometric heating and cooling. The P–T–X framework is then used to discuss the various modes of fluid inclusion entrapment, and how the resulting assemblage textures can be used to interpret the P–T conditions, phase states, and evolution paths of the parent solutions. Finally, quantitative methods are reviewed by which bulk molar volume and composition of CO₂–H₂O fluid inclusions can be determined from microthermometric observations of phase transitions.     

The microdetermination of H2O,CO2, and SO2 in glass using a 1280°C microscope vacuum heating stage,cryopumping, and vapor pressure measurements from 77 to 273 K

A new microscope vacuum heating stage and gas analyzer has been developed for measurement of H₂O, CO₂, SO₂, and noncondensable gas (H₂, CO, N₂, Ar, CH₄, etc.) evolved from samples, particularly natural glass, at temperatures up to 1280°C. The gas evolved upon heating to 1280δC is collected in a liquid nitrogen cold trap. Gas components are identified by the characteristic vapor pressure and temperature ranges over which solid and vapor are in equilibrium during sublimation of individual components. The masses of CO₂, SO₂, and H₂O derived from samples and blanks are calculated using the ideal gas law, the molecular weights of the components, and the gauge constant (i.e. the ratio of the number of moles of a gas to its partial pressure in the constant volume). Results obtained by repeated determinations of H₂O, CO₂, and SO₂ evolved from a submarine basaltic glass from Kilauea volcano, Hawaii, (average sample mass = 3 × 10^?3 g) gave probable errors for the determinations of H₂O (0.23%), CO₂ (0.025%), and S (0.071%) equal to 4, 10, and 8% respectively, of the concentrations. Determinations of H₂O in smaller samples of H₂O-poor basaltic pumice show a linear proportionality (0.063%) between the measured H₂O and the sample mass over the range 0.1 × 10^?6 to 1.7 × 10^?6 g H₂O. Comparisons of H₂O determinations by this technique with those obtained by Penfield, gas chromatic, microcoulometric, and vacuum fusion techniques used elsewhere show reasonably good agreement. Determinations of SO₂ by this technique agree reasonably well X-ray fluorescence and electron microprobe determinations of sulfur. Determinations of CO₂ by the present technique are reproducible but cannot be compared directly to measurements made in other labs because of differences in samples analyzed. The principle advantages of this analytical technique are the very small sample required, the simultaneous determination of H₂O, CO₂, SO₂ and noncondensable gas, the avoidance of calibration procedures dependent on chemical standards, and the visual observations that can be made during sample outgassing.     

Evolution of the chemistry of Fe bearing waters during CO2 degassing

The rates of Fe(II) oxidation and precipitation from groundwater are highly pH dependent. Elevated levels of dissolved CO₂ can depress pH and cause difficulty in removing dissolved Fe and associated metals during treatment of ferruginous water. This paper demonstrates interdependent changes in pH, dissolved inorganic C species, and Fe(II) oxidation rates that occur as a result of the removal (degassing) of CO₂ during aeration of waters discharged from abandoned coal mines. The results of field monitoring of aeration cascades at a treatment facility as well as batchwise aeration experiments conducted using net alkaline and net acidic waters in the UK are combined with geochemical modelling to demonstrate the spatial and temporal evolution of the discharge water chemistry. The aeration cascades removed approximately 67% of the dissolved CO₂ initially present but varying the design did not affect the concentration of Fe(II) leaving the treatment ponds. Continued removal of the residual CO₂ by mechanical aeration increased pH by as much as 2 units and resulted in large increases in the rates of Fe(II) oxidation and precipitation. Effective exsolution of CO₂ led to a reduction in the required lime dose for removal of remaining Fe(II), a very important factor with regard to increasing the sustainability of treatment practices. An important ancillary finding for passive treatment is that varying the design of the cascades had little impact on the rate of CO₂ removal at the flow rates measured.     

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.     

Oxygen sites in hydrous aluminosilicate glasses: the role of Al-O-Al and H2O

We describe here high-field ¹⁷O magic-angle-spinning (MAS) and triple-quantum MAS (3QMAS) NMR spectra for several alkali silicate and Na, K, and Ca aluminosilicate glasses containing up to 10 wt.% water. The H₂O site appears to have a large quadrupolar coupling constant, and its chemical shift increases from Na- to K- glasses, suggesting significant cation-H₂O interactions. In ¹⁷O one-pulse MAS and 3QMAS and ²⁷Al one-pulse NMR experiments, major differences were seen between spectra for anhydrous and hydrous calcium aluminosilicate glasses. The changes in the ¹⁷O MAS spectra can be explained by the addition of an H₂O peak and to the disappearance of an Al-O-Al peak from the ¹⁷O NMR spectrum for the hydrous glass. The ²⁷Al results are consistent with this interpretation.     

Experimental evidence for carbonate precipitation and CO2 degassing during sea ice formation

Chemical and stable carbon isotopic modifications during the freezing of artificial seawater were measured in four 4 m³ tank incubations. Three of the four incubations were inoculated with a nonaxenic Antarctic diatom culture. The 18 days of freezing resulted in 25 to 27 cm thick ice sheets overlying the residual seawater. The ice phase was characterized by a decrease in temperature from −1.9 to −2.2°C in the under-ice seawater down to −6.7°C in the upper 4 cm of the ice sheet, with a concurrent increase in the salinity of the under-ice seawater and brine inclusions of the ice sheet as a result of physical concentration of major dissolved salts by expulsion from the solid ice matrix. Measurements of pH, total dissolved inorganic carbon (C_T) and its stable isotopic composition (δ¹³C_T) all exhibited changes, which suggest minimal effect by biological activity during the experiment. A systematic drop in pH and salinity-normalized C_T by up to 0.37 pH_SWS units and 376 μmol C kg⁻¹ respectively at the lowest temperature and highest salinity part of the ice sheet were coupled with an equally systematic ¹³C enrichment of the C_T. Calculations based on the direct pH and C_T measurements indicated a steady increase in the in situ concentration of dissolved carbon dioxide (CO₂(aq)) with time and increasing salinity within the ice sheet, partly due to changes in the dissociation constants of carbonic acid in the low temperature-high salinity range within sea ice. The combined effects of temperature and salinity on the solubility of CO₂ over the range of conditions encountered during this study was a slight net decrease in the equilibrium CO₂(aq) concentration as a result of the salting-out overriding the increase in solubility with decreasing temperature. Hence, the increase in the in situ CO₂(aq) concentration lead to saturation or supersaturation of the brine inclusions in the ice sheet with respect to atmospheric pCO₂ (≈3.5 × 10⁻⁴ atm). When all physico-chemical processes are considered, we expect CO₂ degassing and carbonate mineral precipitation from the brine inclusions of the ice sheet, which were saturated or highly supersaturated with respect to both the anhydrous (calcite, aragonite, vaterite) and hydrated (ikaite) carbonate minerals.     

Kinetics of H2-O2-H2O redox equilibria and formation of metastable H2O2 under low temperature hydrothermal conditions

Hydrothermal experiments were conducted to evaluate the kinetics of H_2(aq) oxidation in the homogeneous H₂-O₂-H₂O system at conditions reflecting subsurface/near-seafloor hydrothermal environments (55-250 °C and 242-497 bar). The kinetics of the water-forming reaction that controls the fundamental equilibrium between dissolved H_2(aq) and O_2(aq), are expected to impose significant constraints on the redox gradients that develop when mixing occurs between oxygenated seawater and high-temperature anoxic vent fluid at near-seafloor conditions. Experimental data indicate that, indeed, the kinetics of H_2(aq)-O_2(aq) equilibrium become slower with decreasing temperature, allowing excess H_2(aq) to remain in solution. Sluggish reaction rates of H_2(aq) oxidation suggest that active microbial populations in near-seafloor and subsurface environments could potentially utilize both H_2(aq) and O_2(aq), even at temperatures lower than 40 °C due to H_2(aq) persistence in the seawater/vent fluid mixtures. For these H₂-O₂ disequilibrium conditions, redox gradients along the seawater/hydrothermal fluid mixing interface are not sharp and microbially-mediated H_2(aq) oxidation coupled with a lack of other electron acceptors (e.g. nitrate) could provide an important energy source available at low-temperature diffuse flow vent sites.More importantly, when H_2(aq)-O_2(aq) disequilibrium conditions apply, formation of metastable hydrogen peroxide is observed. The yield of H₂O_2(aq) synthesis appears to be enhanced under conditions of elevated H_2(aq)/O_2(aq) molar ratios that correspond to abundant H_2(aq) concentrations. Formation of metastable H₂O₂ is expected to affect the distribution of dissolved organic carbon (DOC) owing to the existence of an additional strong oxidizing agent. Oxidation of magnetite and/or Fe⁺⁺ by hydrogen peroxide could also induce formation of metastable hydroxyl radicals (•OH) through Fenton-type reactions, further broadening the implications of hydrogen peroxide in hydrothermal environments.     

Iron weathering products in a CO2 + (H2O or H2O2) atmosphere: Implications for weathering processes on the surface of Mars

Various iron-bearing primary phases and rocks have been weathered experimentally to simulate possible present and past weathering processes occurring on Mars. We used magnetite, monoclinic and hexagonal pyrrhotites, and metallic iron as it is suggested that meteoritic input to the martian surface may account for an important source of reduced iron. The phases were weathered in two different atmospheres: one composed of CO₂ + H₂O, to model the present and primary martian atmosphere, and a CO₂ + H₂O + H₂O₂ atmosphere to simulate the effect of strong oxidizing agents. Experiments were conducted at room temperature and a pressure of 0.75 atm. Magnetite is the only stable phase in the experiments and is thus likely to be released on the surface of Mars from primary rocks during weathering processes. Siderite, elemental sulfur, ferrous sulfates and ferric (oxy)hydroxides (goethite and lepidocrocite) are the main products in a water-bearing atmosphere, depending on the substrate. In the peroxide atmosphere, weathering products are dominated by ferric sulfates and goethite. A kinetic model was then developed for iron weathering in a water atmosphere, using the shrinking core model (SCM). This model includes competition between chemical reaction and diffusion of reactants through porous layers of secondary products. The results indicate that for short time scales, the mechanism is dominated by a chemical reaction with second order kinetics (k = 7.75 × 10⁻⁵ g⁻¹/h), whereas for longer time scales, the mechanism is diffusion-controlled (De_A = 2.71 × 10⁻¹⁰ m²/h). The results indicate that a primary CO₂- and H₂O-rich atmosphere should favour sulfur, ferrous phases such as siderite or Fe²⁺-sulfates, associated with ferric (oxy)hydroxides (goethite and lepidocrocite). Further evolution to more oxidizing conditions may have forced these precursors to evolve into ferric sulfates and goethite/hematite.