Chemical equilibrium model of solution behavior and solubility in the H-Na-K-OH-Cl-HSO4-SO4-H2O system to high concentration and temperature

This paper describes a chemical model that calculates solute and solvent activities and solid-liquid equilibria in the H-Na-K-OH-Cl-HSO₄-SO₄-H₂O system from dilute to high solution concentration within the 0° to 250°C temperature range. All binary and ternary subsystems are included in the model parameterization. The model is validated by comparing predictions with experimental data, primarily in higher order systems, not used in the parameterization process. Limitations of the model due to data insufficiencies are discussed.The Harvie and Weare (GCA44, 981, 1980) solubility modeling approach, incorporating their implementation of the concentration-dependent specific interaction equations of Pitzer (J. Phys. Chem.77, 268, 1973), is employed. This model expands the variable temperature Na-K-Cl-SO₄-H₂O model of Greenberg and Moller (GCA53, 2503, 1989) by including acid (H₂SO₄, HCl) and base (NaOH, KOH) species. Temperature functions for the chemical potentials of 5 acidic (sodium bisulfate, sodium sesquisulfate, mercallite, potassium sesquisulfate and misenite) and 6 basic (4 sodium hydroxide hydrates and 2 potassium hydroxide hydrates) solid phases are determined.     

A chemical equilibrium model of solution behavior and solubility in the H-Na-K-Ca-OH-Cl-HSO4-SO4-H2O system to high concentration and temperature

This paper describes a chemical model for calculating solute and solvent activities and solid-liquid equilibria in the H-Na-K-Ca-OH-Cl-HSO₄-SO₄-H₂O system from low to high solution concentration within the 0° to 250°C temperature range. The solubility modeling approach of Harvie and Weare (1980) and Harvie et al. (1984), and their implementation of the Pitzer (1973) specific interaction equations are employed. This model expands the variable temperature H-Na-K-OH-Cl-HSO₄-SO₄-H₂O model of Christov and Moller (2004) by evaluating the Ca²⁺-acid/base binary and ternary solution parameters and the chemical potentials of three basic calcium solid phases: Ca(OH)₂(s), CaCl₂.Ca(OH)₂.H₂O(s), and CaCl₂.3Ca(OH)₂.13H₂O and three calcium chloride hydrates (CaCl₂.nH₂O; n = 6, 4, 2). Comparisons of solubility and activity predictions with experimental data used and not used in model parameterization are given. Limitations of the model due to data insufficiencies are discussed.     

A thermodynamic model for calculating methane solubility, density and gas phase composition of methane-bearing aqueous fluids from 273 to 523 K and from 1 to 2000 bar

A thermodynamic model is presented to calculate methane solubility, liquid phase density and gas phase composition of the H₂O-CH₄ and H₂O-CH₄-NaCl systems from 273 to 523 K (possibly up to 573 K), from 1 to 2000 bar and from 0 to 6 mol kg⁻¹ of NaCl with experimental accuracy. By a more strict theoretical approach and using updated experimental data, this model made substantial improvements over previous models: (1) the accuracy of methane solubility in pure water in the temperature range between 273 and 283 K is increased from about 10% to about 5%, but confirms the accuracy of the Duan model [Duan Z., Moller N., Weare J.H., 1992a. Prediction of methane solubilities in natural waters to high ionic strength from 0 to 250 °C and from 0 to 1600 bar. Geochim. Cosmochim. Acta56, 1451-1460] above 283 K up to 2000 bar; (2) the accuracy of methane solubility in the NaCl aqueous solutions is increased from >12% to about 6% on average from 273 K and 1 bar to 523 K and 2000 bar; (3) this model is able to calculate water content in the gas phase and liquid phase density, which cannot be calculated by previous models; and (4) it covers a wider range of temperature and pressure space. With a simple approach, this model is extended to predict CH₄ solubility in other aqueous salt solutions containing Na⁺, K⁺, Mg²⁺, Ca²⁺, Cl⁻ and , such as seawater and geothermal brines, with excellent accuracy. This model is also able to calculate homogenization pressure of fluid inclusions (CH₄-H₂O-NaCl) and CH₄ solubility in water at gas-liquid-hydrate phase equilibrium. A computer code is developed for this model and can be downloaded from the website: www.geochem-model.org/programs.htm.     

Solubility of NaCl in CO2 at high pressure and temperature: First experimental measurements

NaCl solubility in gaseous carbon dioxide has been measured in the pressure range from 30 to 70 MPa at 623 and 673 K. Our originally-designed high pressure apparatus allows in situ sampling of a portion of the fluid phase for chemical analysis. The results indicate that the solubility of NaCl increases with both temperature and pressure, and is about 4-5 orders of magnitude higher than saturated NaCl pressure values at the same temperature conditions (6.02 × 10⁻¹² at 623 K and 1.51 × 10⁻¹⁰ at 673 K). It is also 1-2 orders of magnitude greater than predictions according to the Equation of State of the ternary H₂O-CO₂-NaCl system by Duan, Moeller and Weare [Duan, Z., Moller, N., and Weare, J. H. (1995) Equation of state for the NaCl-H₂O-CO₂ system: prediction of phase equilibria and volumetric properties. Geochim. Cosmochim. Acta59, 2869] and has the opposite pressure dependence. The activity values of NaCl in the vapor phase, calculated from the experiments (with pure molten NaCl as a standard state in the vapor), have been fitted to the Darken Quadratic Formalism: , where, x_NaCl,v is mole the fraction of NaCl in the vapor phase, , , where P is the pressure in MPa and T the absolute temperature. Caution should be exerted while extrapolating this empirical equation far beyond the experimental P-T-compositional range.     

Comment on “Hydrolysis of neptunium(V) at variable temperatures (10-85 °C)” by L. Rao, T.G. Srinivasan, A.Yu. Garnov, P. Zanonato, P. Di Bernardo, and A. Bismondo

In a recent study [Rao, L., Srinivasan, T.G., Garnov, A.Yu., Zanonato, P., Di Bernardo, P., Bismondo, A., 2004. Hydrolysis of neptunium(V) at variable temperatures (10-85 °C). Geochim. Cosmochim. Acta68, 4821-4830.] the hydrolysis of Np(V) was investigated at 10-85 °C by absorption spectroscopy, potentiometry, and microcalorimetry along the titration of Np(V) solutions with tetramethylammonium hydroxide up to pH 10. However, there is strong evidence that the precautions to avoid competing reactions with carbonate were not sufficient and that the measured effects are not caused by the formation of Np(V) hydroxide complexes but primarily by the formation of Np(V) carbonate complexes. The reported equilibrium constants, enthalpies, entropies, and heat capacities for the complexes NpO₂OH(aq) and are severely in error and must not be used for the geochemical modeling of neptunium. If the hydrolysis constants reported by Rao et al. [Rao, L., Srinivasan, T.G., Garnov, A.Yu., Zanonato, P., Di Bernardo, P., Bismondo, A., 2004. Hydrolysis of neptunium(V) at variable temperatures (10-85 °C). Geochim. Cosmochim. Acta68, 4821-4830] are used to calculate neptunium solubilities in alkaline solutions relevant for nuclear waste repositories, the Np(V) concentrations are overestimated by orders of magnitude.     

An experimental study of the solubility and speciation of neodymium (III) fluoride in F-bearing aqueous solutions

The solubility of neodymium (III) fluoride was investigated at temperatures of 150, 200 and 250 °C, saturated water vapor pressure, and a total fluoride concentration (HF°_aq + F⁻) ranging from 2.0 × 10⁻³ to 0.23 mol/l. The results of the experiments show that Nd³⁺ and NdF²⁺ are the dominant species in solution at the temperatures investigated and were used to derive formation constants for NdF²⁺ and a solubility product for NdF₃. The solubility product of NdF₃(logK_sp=loga_Nd³⁺+3loga_F^-) is −24.4 ± 0.2, −22.8 ± 0.1, and −21.5 ± 0.2 at 250, 200 and 150 °C, respectively, and the formation constant of NdF²⁺(logβ=loga_NdF²⁺-loga_Nd³⁺-loga_F^-) is 6.8 ± 0.1, 6.2 ± 0.1, and 5.5 ± 0.2 at 250, 200 and 150 °C, respectively. The results of this study show that published theoretical predictions significantly overestimate the stability of NdF²⁺ and the solubility of NdF₃.The potential impact of the results on natural systems was evaluated for a hypothetical fluid with a composition similar to that responsible for REE mineralization in the Capitan pluton, New Mexico. In contrast to results obtained using the theoretical predictions of Haas [Haas J. R., Shock E. L., and Sassani D. C. (1995) Rare earth elements in hydrothermal systems: estimates of standard partial molal thermodynamic properties of aqueous complexes of the rare earth elements at high pressures and temperatures. Geochim. Cosmochim. Acta59, 4329-4350.], which indicate that NdF²⁺ is the dominant species in solution, calculations employing the data presented in this paper and previously published experimental data for chloride and sulfate species [Migdisov A. A., and Williams-Jones A. E. (2002) A spectrophotometric study of neodymium(III) complexation in chloride solutions. Geochim. Cosmochim. Acta66, 4311-4323; Migdisov A. A., Reukov V. V., and Williams-Jones A. E. (2006) A spectrophotometric study of neodymium(III) complexation in sulfate solutions at elevated temperatures. Geochim. Cosmochim. Acta70, 983-992.] show that neodymium chloride species predominate and that neodymium fluoride species are relatively unimportant. This suggests that accepted models for REE deposits that invoke fluoride complexation as the method of hydrothermal REE transport may need to be re-evaluated.     

The prediction of mineral solubilities in natural waters: The Na-K-Mg-Ca-H-Cl-SO4-OH-HCO3-CO3-CO2-H2O system to high ionic strengths at 25°C

The mineral solubility model of Harvie and Weare (1980) is extended to the eight component system, Na-K-Mg-Ca-H-Cl-SO₄-OH-HCO₃-CO₃-CO₂-H₂O at 25°C to high concentrations. The model is based on the semi-empirical equations of Pitzer (1973) and co-workers for the thermodynamics of aqueous electrolyte solutions. The model is parameterized using many of the available isopiestic, electromotive force, and solubility data available for many of the subsystems. The predictive abilities of the model are demonstrated by comparison to experimental data in systems more complex than those used in parameterization. The essential features of a chemical model for aqueous electrolyte solutions and the relationship between pH and the equilibrium properties of a solution are discussed.     

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.     

Residence times for protons bound to three oxygen sites in the AlO4Al12(OH)24(H2O)12 polyoxocation

Although, the kinetic reactivity of a mineral surface is determined, in part, by the rates of exchange of surface-bound oxygens and protons with bulk solution, there are no elementary rate data for minerals. However, such kinetic measurements can be made on dissolved polynuclear clusters, and here we report lifetimes for protons bound to three oxygen sites on the AlO₄Al₁₂(OH)₂₄(H₂O)₁₂⁷⁺ (Al₁₃) molecule, which is a model for aluminum-hydroxide solids in water. Proton lifetimes were measured using ¹H NMR at pH ∼ 5 in both aqueous and mixed solvents. The ¹H NMR peak for protons on bound waters (η-H₂O) lies near 8 ppm in a 2.5:1 mixture of H₂O/acetone-d₆ and broadens over the temperature range −20 to −5 °C. Extrapolated to 298 K, the lifetime of a proton on a η-H₂O is τ²⁹⁸ ∼ 0.0002 s, which is surprisingly close to the lifetime of an oxygen in the η-H₂O (∼0.0009 s), but in the same general range as lifetimes for protons on fully protonated monomer ions of trivalent metals (e.g., Al(H₂O)₆³⁺). The lifetime is reduced somewhat by acid addition, indicating that there is a contribution from the partly deprotonated Al₁₃ molecule in addition to the fully protonated Al₁₃ at self-buffered pH conditions. Proton lifetimes on the two distinct sets of hydroxyls bridging two Al(III) (μ₂-OH) differ substantially and are much shorter than the lifetime of an oxygen at these sites. The average lifetimes for hydroxyl protons were measured in a 2:1 mixture of H₂O/dmso-d₆ over the temperature range 3.7-95.2 °C. The lifetime of a hydrogen on one of the μ₂-OH was also measured in D₂O. The τ²⁹⁸ values are ∼0.013 and ∼0.2 s in the H₂O/dmso-d₆ solution and the τ²⁹⁸ value for the μ₂-OH detectable in D₂O is τ²⁹⁸ ∼ 0.013 s. The ¹H NMR peak for the more reactive μ₂-OH broadens slightly with acid addition, indicating a contribution from an exchange pathway that involves a proton or hydronium ion. These data indicate that surface protons on minerals will equilibrate with near-surface waters on the diffusional time scale.     

Silicon isotopic fractionation of CAI-like vacuum evaporation residues

Calcium-, aluminum-rich inclusions (CAIs) are often enriched in the heavy isotopes of magnesium and silicon relative to bulk solar system materials. It is likely that these isotopic enrichments resulted from evaporative mass loss of magnesium and silicon from early solar system condensates while they were molten during one or more high-temperature reheating events. Quantitative interpretation of these enrichments requires laboratory determinations of the evaporation kinetics and associated isotopic fractionation effects for these elements. The experimental data for the kinetics of evaporation of magnesium and silicon and the evaporative isotopic fractionation of magnesium is reasonably complete for Type B CAI liquids (Richter F. M., Davis A. M., Ebel D. S., and Hashimoto A. (2002) Elemental and isotopic fractionation of Type B CAIs: experiments, theoretical considerations, and constraints on their thermal evolution. Geochim. Cosmochim. Acta66, 521-540; Richter F. M., Janney P. E., Mendybaev R. A., Davis A. M., and Wadhwa M. (2007a) Elemental and isotopic fractionation of Type B CAI-like liquids by evaporation. Geochim. Cosmochim. Acta71, 5544-5564.). However, the isotopic fractionation factor for silicon evaporating from such liquids has not been as extensively studied. Here we report new ion microprobe silicon isotopic measurements of residual glass from partial evaporation of Type B CAI liquids into vacuum. The silicon isotopic fractionation is reported as a kinetic fractionation factor, α_Si, corresponding to the ratio of the silicon isotopic composition of the evaporation flux to that of the residual silicate liquid. For CAI-like melts, we find that α_Si = 0.98985 ± 0.00044 (2σ) for ²⁹Si/²⁸Si with no resolvable variation with temperature over the temperature range of the experiments, 1600-1900 °C. This value is different from what has been reported for evaporation of liquid Mg₂SiO₄ (Davis A. M., Hashimoto A., Clayton R. N., and Mayeda T. K. (1990) Isotope mass fractionation during evaporation of Mg₂SiO₄. Nature347, 655-658.) and of a melt with CI chondritic proportions of the major elements (Wang J., Davis A. M., Clayton R. N., Mayeda T. K., and Hashimoto A. (2001) Chemical and isotopic fractionation during the evaporation of the FeO-MgO-SiO₂-CaO-Al₂O₃-TiO₂-REE melt system. Geochim. Cosmochim. Acta65, 479-494.). There appears to be some compositional control on α_Si, whereas no compositional effects have been reported for α_Mg. We use the values of α_Si and α_Mg, to calculate the chemical compositions of the unevaporated precursors of a number of isotopically fractionated CAIs from CV chondrites whose chemical compositions and magnesium and silicon isotopic compositions have been previously measured.     

Application of the Pitzer ion interaction model to isopiestic data for the Fe2(SO4)3-H2SO4-H2O system at 298.15 and 323.15 K

Recent isopiestic studies of the Fe₂(SO₄)₃-H₂SO₄-H₂O system at 298.15 K are represented with an extended version of Pitzer’s ion interaction model. The model represents osmotic coefficients for aqueous {(1 − y)Fe₂(SO₄)₃ + yH₂SO₄} mixtures from 0.45 to 3.0 m at 298.15 K and 0.0435 ? y ? 0.9370. In addition, a slightly less accurate representation of a more extended molality range to 5.47 m extends over the same y values, translating to a maximum ionic strength of 45 m. Recent isopiestic data for the system at 323.15 K are represented with the extended Pitzer model over a limited range in molality and solute fraction. These datasets are also represented with the usual “3-parameter” version of Pitzer’s model so that it may be incorporated in geochemical modeling software, but is a slightly less accurate representation of thermodynamic properties for this system. Comparisons made between our ion interaction model and available solubility data display partial agreement for rhomboclase and significant discrepancy for ferricopiapite. The comparisons highlight uncertainty remaining for solubility predictions in this system as well as the need for additional solubility measurements for Fe³⁺-bearing sulfate minerals. The resulting Pitzer ion interaction models provide an important step toward an accurate and comprehensive representation of thermodynamic properties in this geochemically important system.     

Decomposition reactions of magnesium sulfate hydrates and phase equilibria in the MgSO4-H2O and Na-Mg-Cl-SO4-H2O systems with implications for Mars

We report new measurements of equilibrium relative humidities for stable and metastable hydration-dehydration equilibria involving several magnesium sulfates in the MgSO₄·nH₂O series. We also report a comprehensive thermodynamic treatment of the system including solution properties and experimental data from the published literature, i.e. solubilities, heat capacities and additional decomposition humidities. While for some magnesium sulfate hydrates solubility data in the binary system MgSO₄-H₂O are sparse, there is a reasonable database of solubility measurements of these hydrates in the ternary MgCl₂-MgSO₄-H₂O and the quaternary reciprocal Na⁺-Mg²⁺-Cl⁻-SO₄²⁻-H₂O systems. To make these data suitable for the determination of solubility products, we parameterized a Pitzer ion interaction model for the calculation of activity coefficients and water activities in mixed solutions of these systems and report the ion interaction parameters for the Na⁺-Mg²⁺-Cl⁻-SO₄²⁻-H₂O system. The model predicted solubilities in the reciprocal system are in very good agreement with experimental data. Using all available experimental data and the solution model an updated phase diagram of the MgSO₄-H₂O system covering the whole temperature range from about 170 to 473 K is established. This treatment includes MgSO₄·H₂O (kieserite), MgSO₄·4H₂O (starkeyite), MgSO₄·5H₂O (pentahydrite), MgSO₄·6H₂O (hexahydrite), MgSO₄·7H₂O (epsomite) and MgSO₄·11H₂O (meridianiite). It is shown that only kieserite, hexahydrite, epsomite and meridianiite show fields of stable existence while starkeyite and pentahydrite are always metastable. Due to sluggish kinetics of kieserite formation, however, there is a rather extended field of metastable existence of starkeyite which makes this solid a major product in dehydration reactions. The model predicted behavior of the magnesium sulfates is in excellent agreement with observations reported in the literature under terrestrial temperature and relative humidity conditions. We also discuss the implications of the new phase diagram for sulfates on Mars.     

Iron isotope fractionation between pyrite (FeS2), hematite (Fe2O3) and siderite (FeCO3): A first-principles density functional theory study

In addition to equilibrium isotopic fractionation factors experimentally derived, theoretical predictions are needed for interpreting isotopic compositions measured on natural samples because they allow exploring more easily a broader range of temperature and composition. For iron isotopes, only aqueous species were studied by first-principles methods and the combination of these data with those obtained by different methods for minerals leads to discrepancies between theoretical and experimental isotopic fractionation factors. In this paper, equilibrium iron isotope fractionation factors for the common minerals pyrite, hematite, and siderite were determined as a function of temperature, using first-principles methods based on the density functional theory (DFT). In these minerals belonging to the sulfide, oxide and carbonate class, iron is present under two different oxidation states and is involved in contrasted types of interatomic bonds. Equilibrium fractionation factors calculated between hematite and siderite compare well with the one estimated from experimental data (ln α⁵⁷Fe/⁵⁴Fe = 4.59 ± 0.30‰ and 5.46 ± 0.63‰ at 20 °C for theoretical and experimental data, respectively) while those for Fe(III)_aq-hematite and Fe(II)_aq-siderite are significantly higher that experimental values. This suggests that the absolute values of the reduced partition functions (β-factors) of aqueous species are not accurate enough to be combined with those calculated for minerals. When compared to previous predictions derived from Mössbauer or INRXS data [Polyakov V. B., Clayton R. N., Horita J. and Mineev S. D. (2007) Equilibrium iron isotope fractionation factors of minerals: reevaluation from the data of nuclear inelastic resonant X-ray scattering and Mössbauer spectroscopy. Geochim. Cosmochim. Acta71, 3833-3846], our iron β-factors are in good agreement for siderite and hematite while a discrepancy is observed for pyrite. However, the detailed investigation of the structural, electronic and vibrational properties of pyrite as well as the study of sulfur isotope fractionation between pyrite and two other sulfides (sphalerite and galena) indicate that DFT-derived β-factors of pyrite are as accurate as for hematite and siderite. We thus suggest that experimental vibrational density of states of pyrite should be re-examined.     

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

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.     

Europium retention onto clay minerals from 25 to 150 °C: Experimental measurements, spectroscopic features and sorption modelling

The sorption of Eu(III) onto kaolinite and montmorillonite was investigated up to 150 °C. The clays were purified samples, saturated with Na in the case of montmorillonite. Batch experiments were conducted at 25, 40, 80 and 150 °C in 0.5 M NaClO₄ solutions to measure the distribution coefficients (Kd) of Eu as a trace element (<10⁻⁶ mol/L) between the solution and kaolinite. For the Na-montmorillonite, we used Kd results from a previous study [Tertre, E., Berger, G., Castet, S., Loubet, M., Giffaut, E., 2005. Experimental study of adsorption of Ni²⁺, Cs⁺ and Ln³⁺ onto Na-montmorillonite up to 150 °C. Geochim. Cosmochim. Acta69, 4937-4948] obtained under exactly the same conditions. The number and nature of the Eu species sorbed onto both clay minerals were investigated by time resolved laser fluorescence spectroscopy (TRLFS) in specific experiments in the same temperature range. We identified a unique inner-sphere complex linked to the aluminol sites in both clays, assumed to be AlOEu²⁺ at the edge of the particles, and a second exchangeable outer-sphere complex for montmorillonite, probably in an interlayer position. The Kd values were used to adjust the parameters of a surface complexation model (DLM: diffuse layer model) from 25 to 150 °C. The number of Eu complexes and the stoichiometry of reactions were constrained by TRLFS. The acidity constants of the amphoteric aluminol sites were taken from another study [Tertre, E., Castet, S., Berger, G., Loubet, M., Giffaut, E. Acid/base surface chemistry of kaolinite and Na-montmorillonite at 25 and 60 °C: experimental study and modelling. Geochim. Cosmochim. Acta, in press], which integrates the influence of the negative structural charge of clays on the acid/base properties of edge sites as a function of temperature and ionic strength. The results of the modelling show that the observed shift of the sorption edge towards low pH with increasing temperature results solely from the contribution of the AlOEu²⁺ edge complexes. Finally, we successfully tested the performance of our model by confronting the predictions with experimental Kd data. We used our own data obtained at lower ionic strength (previous study) or higher suspension density and higher starting concentration (TRLFS runs, this study), as well as published data from other experimental studies [Bradbury, M.H., Baeyens, B., 2002. Sorption of Eu on Na and Ca-montmorillonite: experimental investigations and modeling with cation exchange and surface complexation. Geochim. Cosmochim. Acta66, 2325-2334; Kowal-Fouchard, A., 2002. Etude des mécanismes de rétention des ions U(IV) et Eu(III) sur les argiles: influence des silicates. Ph.D. Thesis, Université Paris Sud, France, 330p].     

The activity coefficients of Fe(III) hydroxide complexes in NaCl and NaClO4 solutions

The osmotic coefficients of FeCl₃ at 25 °C from 0.15 to 1.7 m [Rumyantsev et al., Z. Phys. Chem., 218, 1089-1127, 2004] have been used to determine the Pitzer parameters (β⁽⁰⁾, β⁽¹⁾ and C^?) for FeCl₃. Since the differences in the Pitzer coefficients of rare earths in NaCl and NaClO₄ are small, the values of Fe(ClO₄)₃ have been estimated using the differences between La(ClO₄)₃ and LaCl₃. The Pitzer coefficients for FeCl₃ combined with enthalpy and heat capacity data for the rare earths can be used to estimate the activity coefficients of Fe³⁺ in NaCl over a wide range of temperatures (0 to 50 °C) and ionic strength (0 to 6 m).The activity coefficients of Fe³⁺ in NaCl and NaClO₄ solutions have been used to determine the activity coefficients of Fe(OH)²⁺ in these solutions from the measured first hydrolysis constants of Fe³⁺ [Byrne et al., Mar. Chem., 97, 34-48, 2005]. The activity coefficients of , Fe(OH)₃ and from 0 to 50 °C have also been determined from the solubility measurements of Fe(III) in NaCl solutions [Liu and Millero, Geochim. Cosmochim Acta, 63, 3487-3497, 1999]. These activity coefficients have been fitted to the Pitzer equations. These results can be used to estimate the speciation of Fe(III) with OH⁻ in natural waters with high concentrations of NaCl from 0 to 50 °C.     

Experimental investigation of the solubility of albite and jadeite in H2O, with paragonite + quartz at 500 and 600 °C, and 1-2.25 GPa

The solubilities of the assemblages albite + paragonite + quartz and jadeite + paragonite + quartz in H₂O were determined at 500 and 600 °C, 1.0-2.25 GPa, using hydrothermal piston-cylinder methods. The three minerals are isobarically and isothermally invariant in the presence of H₂O, so fluid composition is uniquely determined at each pressure and temperature. A phase-bracketing approach was used to achieve accurate solubility determinations. Albite + quartz and jadeite + quartz dissolve incongruently in H₂O, yielding residual paragonite which could not be retrieved and weighed. Solution composition fixed by the three-mineral assemblage at a given pressure and temperature was therefore bracketed by adding NaSi₃O_6.5 glass in successive experiments, until no paragonite was observed in run products. Solubilities derived from experiments bounding the appearance of paragonite thus constrain the equilibrium fluid composition. Results indicate that, at a given pressure, Na, Al, and Si concentrations are higher at 600 °C than at 500 °C. At both 500 and 600 °C, solubilities of all three elements increase with pressure in the albite stability field, to a maximum at the jadeite-albite-quartz equilibrium. In the jadeite stability field, element concentrations decline with continued pressure increase. At the solubility maximum, Na, Al, and Si concentrations are, respectively, 0.16, 0.05, and 0.48 molal at 500 °C, and 0.45, 0.27, and 1.56 molal at 600 °C. Bulk solubilities are 3.3 and 10.3 wt% oxides, respectively. Observed element concentrations are everywhere greater than those predicted from extrapolated thermodynamic data for simple ions, monomers, ion pairs, and the silica dimer. The measurements therefore require the presence of additional, polymerized Na-Al-Si-bearing species in the solutions. The excess solubility is >50% at all conditions, indicating that polymeric structures are the predominant solutes in the P-T region studied. The solubility patterns likely arise from combination of the large solid volume change associated with the albite-jadeite-quartz equilibrium and the rise in Na-Al-Si polymerization with approach to the hydrothermal melting curves of albite + quartz and jadeite + quartz. Our results indicate that polymerization of Na-Al-Si solutes is a fundamental aspect of fluid-rock interaction at high pressure. In addition, the data suggest that high-pressure metamorphic isograds can impose unexpected controls on metasomatic mass transfer, that significant metasomatic mass transfer prior to melting should be considered in migmatitic terranes, and that polymeric complexes may be an important transport agent in subduction zones.     

An experimental study of the solubility of molybdenum in H2O and KCl-H2O solutions from 500 °C to 800 °C, and 150 to 300 MPa

The solubility of molybdenum (Mo) was determined at temperatures from 500 °C to 800 °C and 150 to 300 MPa in KCl-H₂O and pure H₂O solutions in cold-seal experiments. The solutions were trapped as synthetic fluid inclusions in quartz at experimental conditions, and analyzed by laser ablation inductively coupled plasma mass spectrometry (LA ICPMS).Mo solubilities of 1.6 wt% in the case of KCl-bearing aqueous solutions and up to 0.8 wt% in pure H₂O were found. Mo solubility is temperature dependent, but not pressure dependent over the investigated range, and correlates positively with salinity (KCl concentration). Molar ratios of ∼1 for Mo/Cl and Mo/K are derived based on our data. In combination with results of synchrotron X-ray absorption spectroscopy of individual fluid inclusions, it is suggested that Mo-oxo-chloride complexes are present at high salinity (>20 wt% KCl) and ion pairs at moderate to low salinity (<11 wt% KCl) in KCl-H₂O aqueous solutions. Similarly, in the pure H₂O experiments molybdic acid is the dominant species in aqueous solution. The results of these hydrothermal Mo experiments fit with earlier studies conducted at lower temperatures and indicate that high Mo concentrations can be transported in aqueous solutions. Therefore, the Mo concentration in aqueous fluids seems not to be the limiting factor for ore formation, whereas precipitation processes and the availability of sulfur appear to be the main controlling factors in the formation of molybdenite (MoS₂).     

New data on the kinetics of the dissolution of alkali feldspars at 200°C in CO2 charged water

The kinetics of the dissolution of a (Na. K) sanidine in CO₂-charged water, were studied experimentally at 200°C as a function of the surface area s and of the time t. The molalities of Na⁺ and K⁺ have been plotted against the product st. As for albite and adularia (Lagache, Bull. Soc. Fr. Minéral. Cristallogr. 88, 223–253, 1965), the rates of dissolution of Na² and K⁺ are smooth continuous functions of st, which implies that they are controlled by the composition of the solution.The comparison between the dissolution of pure sodic or potassic feldspars and that of an intermediate feldspar shows that the sanidine dissolves as if it were composed of albite and K-feldspar grains in the proportions corresponding to its composition.Theoretical considerations presented by Helgeson (Geochim. Cosmochim. Acta35, 421–169, 1971; The Feldspars, pp. 184–217, 1972) and Pa?es (Geochim. Cosmochim. Acta37, 2641–2663, 1973) are reviewed: both had suggested, arguing from my first experimental results, that the dissolution could be described by a process of diffusional mass transfer through a surface layer of reaction products.The present experiments do not agree with such an interpretation of the mechanism of dissolution.