CO2-H2O mixtures in the geological sequestration of CO2. I. Assessment and calculation of mutual solubilities from 12 to 100°C and up to 600 bar

Evaluating the feasibility of CO₂ geologic sequestration requires the use of pressure-temperature-composition (P-T-X) data for mixtures of CO₂ and H₂O at moderate pressures and temperatures (typically below 500 bar and below 100°C). For this purpose, published experimental P-T-X data in this temperature and pressure range are reviewed. These data cover the two-phase region where a CO₂-rich phase (generally gas) and an H₂O-rich liquid coexist and are reported as the mutual solubilities of H₂O and CO₂ in the two coexisting phases. For the most part, mutual solubilities reported from various sources are in good agreement. In this paper, a noniterative procedure is presented to calculate the composition of the compressed CO₂ and liquid H₂O phases at equilibrium, based on equating chemical potentials and using the Redlich-Kwong equation of state to express departure from ideal behavior. The procedure is an extension of that used by King et al. (1992), covering a broader range of temperatures and experimental data than those authors, and is readily expandable to a nonideal liquid phase. The calculation method and formulation are kept as simple as possible to avoid degrading the performance of numerical models of water-CO₂ flows for which they are intended. The method is implemented in a computer routine, and inverse modeling is used to determine, simultaneously, (1) new Redlich-Kwong parameters for the CO₂-H₂O mixture, and (2) aqueous solubility constants for gaseous and liquid CO₂ as a function of temperature. In doing so, mutual solubilities of H₂O from 15 to 100°C and CO₂ from 12 to 110°C and up to 600 bar are generally reproduced within a few percent of experimental values. Fugacity coefficients of pure CO₂ are reproduced mostly within one percent of published reference data.     

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.     

Synthetic Fluid Inclusions in the Systems NaCl-H2O and NaCl-CO2-H2O: Dissolution Temperatures of Halite

Abstract. This study examined the effect of CO₂ on NaCl solubility in hydrothermal fluid, with the synthetic fluid inclusion technique. Fluid inclusions of 30–40 wt% NaCl and 5 mol % CO₂ were synthesized, and their halite dissolution temperatures, T_m(halite), were measured. The solubilities of NaCl in CO₂-bearing aqueous fluid were obtained at 160–320C under vapor-saturated pressures. The T_m(halite) value in aqueous fluid with 5 mol % CO₂ obtained in this study agrees with that of Schmidt et al. (1995), showing that 5 mol % CO₂ reduces the solubility of NaCl by about 1 wt%.
Calculation of magnetite solubility suggests that 5–10 mol % CO2 decreases magnetite solubility by 4.5–8.9 % relative to the magnetite solubility in CO2-free solution. Therefore, an increase of CO2 content in ore-forming solutions may cause deposition of iron minerals and produce ore deposits.     

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.     

Amorphous silica solubilities—I. Behavior in aqueous sodium nitrate solutions; 25–300°C, 0–6 molal

The solubility of amorphous silica was obtained in aqueous sodium nitrate solutions up to six molal and at temperatures from 25 to 300°C. It was expected that solubilities in aqueous sodium chloride solutions would be similar. At 25°C, the solubility of amorphous silica is lowered from that in water to 0.00086 m in 6.12 m sodium nitrate, or a decrease of 60%. At 300°C, the corresponding decrease is only 27% from a solubility of 0.0269 m in H₂O. From the change in solubility with temperature at a given constant molality of sodium nitrate, the molal heat of solution over the range, 100 to 300°C, increases from + 2.93 kcal mol^?1 in water to + 3.64 kcal mol^?1 in 6m sodium nitrate. The value approaches a constant of +3.8 kcal mol^?1 as sodium nitrate approaches saturation at 10.8 molal.     

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.     

A thermodynamic model for the solubility of barite and celestite in electrolyte solutions and seawater to 200°C and to 1 kbar

This paper describes a model for barite and celestite solubilities in the Na–K–Ca–Mg–Ba–Sr–Cl–SO₄–H₂O system to 200°C and to 1 kbar. It is based on Pitzer's ion interaction model for the thermodynamic properties of the aqueous phase and on values of the solubility products of the solids revised in this work. It is shown how barite and celestite solubilities in electrolyte solutions can be accurately predicted as a function of temperature and pressure from previously determined Pitzer's parameters. The equilibrium constant for the BaSO₄(aq) ion pair dissociation reaction is calculated from recently reported barite solubility in Na₂SO₄ solutions from 0 to 80°C. Pressure corrections are evaluated through partial molal volume calculations and are partially validated by comparing model predictions to measured barite and celestite solubilities in pure water to 1 kbar and in NaCl solutions to 500 bars. The model is then used to investigate the tendency of ion pairing of Ca, Sr and Ba with sulfate in seawater. Finally, the activity coefficient of aqueous barium sulfate in seawater is calculated for temperature, pressure and salinity values found in the ocean and compared to published values.     

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.     

Mineral solution equilibria—II. An experimental study of mineral solubilities and the thermodynamic properties of aqueous CaCl2 in the system CaO-SiO2-H2O-HCl

Speciation of aqueous calcium chloride and the solubility of wollastonite represented by the reaction wollastonite + 2HCl° → CaCl₂° + quartz + H₂O were experimentally investigated at 1 and 2 kbar in the range 425–600°C using rapid-quench hydrothermal techniques and a modified Ag + AgCl buffer technique (Frantz and Popp, 1979). Variation in the measured concentration in HCl° as a function of total dissolved calcium was used to identify associated aqueous CaCl₂° as the predominant calcium species in the fluid at temperatures above 500°C at 2 kbar. The data were used to calculate the equilibrium constant for the above reaction as a function of temperature and pressure, from which the difference in Gibbs free energy of formation between CaCl₂° and HCl° at 1 and 2 kbar, 450°–600°C was calculated. Solubility constants for minerals in the system MgO-CaO-SiO₂-H₂O-HCl-CO₂ were calculated using the data from this study and from Frantz and Popp (1979). Calculated mineral solubilities were used to calculate the solution compositions and solid alteration products resulting from interactions of a Ca-Mg silicate mineral (diopside) with hydrothermal solutions containing a range of different total chloride concentrations. High total chloride (2.0 m) in the solution results in Si-Mg enrichment in the solids and Ca enrichment in the fluid, whereas low total chloride (0.008 m) results in Mg enrichment in the solids and Ca-Si enrichment in the fluid.     

Melting relationships in the systems CaOCO2 and MgOCO2 to 33 kilobars

The melting temperatures of calcite and magnesite in the presence of excess CO₂ have been measured using Ag₂C₂O₄ in sealed capsules m a piston-cylinder apparatus. At 27 kbar, 11.5 wt % CO₂ dissolves in molten CaCO₂, depressing the freezing temperature from 1610 to 1505°C; and 6.5 wt % CO₂ dissolves in molten MgCO₃, depressing the freezing temperature from 1590 to 1510°C. The eutectic between calcite and lime was located at 1385°C at 27 kbar. These and other new results, combined with previously published data, permit completion of PT diagrams for the systems CaO-CO₂ and MgO-CO₂ from 1 bar to 35 kbar. The dissociation curve for each carbonate terminates at an invariant point where melting begins, at 40 bars and 1230°C for CaO-CO₂ and 23 kbar and 1550°C for MgO-CO₂ The differences between the two systems are explained by the different solubilities of CO₂ in the invariant liquids consequent upon the large pressure difference between the locations of these two invariant points. The results show that the temperatures for the beginning of melting of carbonates in the asthenosphere are lowered by about 100°C in the presence of CO₂.     

Phase equilibria in the system CO2-H2O I: New equilibrium relations at low temperatures

Graphical analysis of free-energy relationships involving binary quadruple points and their associated univariant equilibria in the system CO₂-H₂O suggests the presence of at least 2 previously unrecognized quadruple points and a degenerate binary invariant point involving an azeotrope between CO₂-rich gas and liquid. Thermodynamic data extracted from the equilibrium involving clathrate (hydrate), gas, and ice (H = G+I) are employed along with published data to calculate the P-T range of the 3-ice equilibrium curve, S+I = H, where S is solid CO₂. This equilibrium curve intersects the H = G+I curve approximately where the latter curve intersects the S+H = G curve, thus confirming the existence of one of the inferred quadruple points involving the phases S, G, H, and I. Recognition of some binary equilibria probably have been hampered by extremely low mutual solubilities of CO₂ and H₂O in the fluids phases which, for example, render the S+H = G virtually indistinguishable from the CO₂-sublimation curve.To make the published portion of the L(liquid CO₂)-G-H equilibrium “connect” with the other new quadruple point involving S, L, G, and H, it is necessary to change the sense of the equilibrium from L = G+H at higher pressures to L+H = G at lower pressures by positing a L = G azeotrope at very low concentrations of H₂O. At the low-pressure origin of the azeotrope, which is only a few bars above the CO₂-triple point, the azeotrope curve intersects the 3-phase curve tangentially, creating a degenerate invariant point at which the 3-phase equilibrium changes from L+H = G at lower pressures to L = G+H at higher pressures. The azeotrope curve is offset at slightly lower temperature from the L = G+H curve until the 3-phase equilibrium terminates at the quadruple point involving G, L, H, and W (water). With further increase in pressure the azeotrope curve tracks the L = G+W equilibrium and apparently terminates at a critical end point in close proximity to critical endpoints for the CO₂-saturation curve and the L = G+W curve.Thermodynamic data for clathrate extracted from the slope of the H = G+I curve are consistent with a solid-state phase transformation in CO₂-clathrate between 235 and 255 K. Published work shows that the type-I clathrate phase, whose atomic structure is a framework of water molecules with CO₂ molecules situated in large “guest” sites within the framework, is variable in composition with ∼1 guest site vacancy per unit cell at the high-temperature limit of its stability; the number of water molecules, however, remains constant. The formula (CO₂)_8-y·46H₂O, where y is the number of vacancies per unit cell, is in keeping with the atomic structure, whereas the traditional formula, CO₂·nH₂O, where n (hydration number) = 5.75, is misleading.Ambient P-T conditions in the Antarctic and Greenland ice sheets are compatible with sequestering large amounts of carbon as liquid CO₂ and/or clathrate.     

Solubility of grossular, Ca3Al2Si3O12, in H2O-NaCl solutions at 800 °C and 10 kbar, and the stability of garnet in the system CaSiO3-Al2O3-H2O-NaCl

The solubility and stability of synthetic grossular were determined at 800 °C and 10 kbar in NaCl-H₂O solutions over a large range of salinity. The measurements were made by evaluating the weight losses of grossular, corundum, and wollastonite crystals equilibrated with fluid for up to one week in Pt capsules and a piston-cylinder apparatus. Grossular dissolves congruently over the entire salinity range and displays a large solubility increase of 0.0053 to 0.132 molal Ca₃Al₂Si₃O₁₂ with increasing NaCl mole fraction (X_NaCl) from 0 to 0.4. There is thus a solubility enhancement 25 times the pure H₂O value over the investigated range, indicating strong solute interaction with NaCl. The Ca₃Al₂Si₃O₁₂ mole fraction versus NaCl mole fraction curve has a broad plateau between X_NaCl = 0.2 and 0.4, indicating that the solute products are hydrous; the enhancement effect of NaCl interaction is eventually overtaken by the destabilizing effect of lowering H₂O activity. In this respect, the solubility behavior of grossular in NaCl solutions is similar to that of corundum and wollastonite. There is a substantial field of stability of grossular at 800 °C and 10 kbar in the system CaSiO₃-Al₂O₃-H₂O-NaCl. At high Al₂O₃/CaSiO₃ bulk compositions the grossular + fluid field is limited by the appearance of corundum. Zoisite appears metastably with corundum in initially pure H₂O, but disappears once grossular is nucleated. At X_NaCl = 0.3, however, zoisite is stable with corundum and fluid; this is the only departure from the quaternary system encountered in this study. Corundum solubility is very high in solutions containing both NaCl and CaSiO₃: Al₂O₃ molality increases from 0.0013 in initially pure H₂O to near 0.15 at X_NaCl = 0.4 in CaSiO₃-saturated solutions, a >100-fold enhancement. In contrast, addition of Al₂O₃ to wollastonite-saturated NaCl solutions increases CaSiO₃ molality by only 12%. This suggests that at high pH (quench pH is 11-12), the stability of solute Ca chloride and Na-Al ± Si complexes account for high Al₂O₃ solubility, and that Ca-Al ± Si complexes are minor. The high solubility and basic dissolution reaction of grossular suggest that Al may be a very mobile component in calcareous rocks in the deep crust and upper mantle when migrating saline solutions are present.     

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

Carbonation of alkaline paper mill waste to reduce CO2 greenhouse gas emissions into the atmosphere

The global warming of Earth’s near-surface, air and oceans in recent decades is a direct consequence of anthropogenic emission of greenhouse gases into the atmosphere such as CO₂, CH₄, N₂O and CFCs. The CO₂ emissions contribute approximately 60% to this climate change. This study investigates experimentally the aqueous carbonation mechanisms of an alkaline paper mill waste containing about 55 wt% portlandite (Ca(OH)₂) as a possible mineralogical CO₂ sequestration process. The overall carbonation reaction includes the following steps: (1) Ca release from portlandite dissolution, (2) CO₂ dissolution in water and (3) CaCO₃ precipitation. This CO₂ sequestration mechanism was supported by geochemical modelling of final solutions using PHREEQC software, and observations by scanning electron microscope and X-ray diffraction of final reaction products. According to the experimental protocol, the system proposed would favour the total capture of approx. 218 kg of CO₂ into stable calcite/ton of paper waste, independently of initial CO₂ pressure. The final product from the carbonation process is a calcite (ca. 100 wt%)-water dispersion. Indeed, the total captured CO₂ mineralized as calcite could be stored in degraded soils or even used for diverse industrial applications. This result demonstrates the possibility of using the alkaline liquid–solid waste for CO₂ mitigation and reduction of greenhouse effect gases into the atmosphere.     

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.     

An empirical solvus for CO2-H2O-2.6 wt% salt

The solvus in the system CO₂-H₂O-2.6 wt% NaCl-equivalent was determined by measuring temperature of homogenization in fluid inclusions which contained variable $\text{CO}{}_2\text{H}{}_2\text{O}$ but the same amount of salt dissolved in the aqueous phase at room temperature. The critical point of the solvus is at 340 ± 5°C, at pressures between 1 and 2 kbar; this is about 65°C higher than for the pure CO₂-H₂O system. The solvus is assymetrical, with a steeper H₂O-rich limb and with the critical point at mole fraction of water between 0.65 and 0.8.     

Combined microthermometric and Raman spectroscopic technique to determine the salinity of H2O-CO2-NaCl fluid inclusions based on clathrate melting

Fluid inclusions approximated by the system H₂O-CO₂-NaCl are common in many geologic environments. In order to apply microthermometric data from these inclusions to infer P-T (pressure-temperature) trapping conditions, the composition of the inclusions, including the salinity, must be known. Normally, salinities of aqueous inclusions are determined from ice-melting temperatures obtained during microthermometry. However, when CO₂-bearing aqueous fluid inclusions are cooled they often form a hydrate that incorporates H₂O into the structure, and salinities estimated from ice-melting temperatures are therefore higher than the actual salinity. A technique that combines data from Raman spectroscopic and microthermometric analyses of individual inclusions was developed to determine the salinity of CO₂-bearing aqueous inclusions based on measured clathrate melting temperatures and CO₂ pressures obtained from Raman analyses. In this study, the pressure within inclusions was determined using Raman spectroscopy based on the splitting of the Fermi diad of CO₂, measured at the clathrate melting temperature. The CO₂ densities (and pressures) predicted by the equation developed in this study are in relatively good agreement with previously published equations, except for very low densities and correspondingly low pressures. The combined Raman spectroscopy - microthermometry technique thus provides both the temperature and the pressure in the inclusion at clathrate melting. For inclusions in which the clathrate melts in the presence of CO₂ liquid, the salinity can be determined with a precision of a few tenths of a wt% NaCl, whereas for inclusions in which clathrate melts in the presence of CO₂ vapor the salinity error may be a few wt% NaCl. Applying the method to synthetic fluid inclusions with known salinity suggests that the technique is valid for determining salinity of H₂O-CO₂-NaCl fluid inclusions in which clathrate melts in the presence of liquid CO₂ only or vapor CO₂ only.     

Solubility of the buffer assemblage pyrite + pyrrhotite + magnetite in NaCl solutions from 200 to 350°C

In the design of hydrothermal solubility studies it is important that the system be completely defined chemically. If the solubilities of minerals containing m metallic elements are to be determined in hydrothermal NaCl solutions, the phase rule requires that a total of m + 6 independent intensive parameters be controlled or measured in order to determine completely the system.In this study the solubility of the univariant assemblage pyrite + pyrrhotite + magnetite has been determined in vapor saturated hydrothermal solutions from 200 to 350°C for NaCl concentrations ranging from 0.0 to 5.0 molal. At any temperature, oxygen and sulfur fugacities were buffered by the chosen assemblage. System pH was determined from excess CO₂ partial pressures and computed ionic equilibria. Equilibrium constants were calculated by regression analysis of solubility data. The results show that more than 10 ppm of each mineral can dissolve in typical hydrothermal solutions under geologically realistic conditions. Solubilities were best represented by the species Fe²⁺ and FeCl⁺ at 200 and 250°C; Fe²⁺, FeCl⁺ and FeCl₂⁰ at 300°C; and Fe²⁺ and FeCl₂⁰ at 350°C. Ore deposition would occur by lowering temperature, diluting chloride concentration, or by raising pH through wall rock alteration reactions.