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.     

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.     

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.     

The measurement of sulfate mineral solubilities in the Na-K-Ca-Cl-SO4-H2O system at temperatures of 100, 150 and 200°C

At T > 100°C development of thermodynamic models suffers from missing experimental data, particularly for solubilities of sulfate minerals in mixed solutions. Solubilities in Na⁺-K⁺-Ca²⁺-Cl⁻-SO₄²⁻/H₂O subsystems were investigated at 150, 200°C and at selected compositions at 100°C. The apparatus used to examine solid-liquid phase equilibria under hydrothermal conditions has been described.In the system NaCl-CaSO₄-H₂O the missing anhydrite (CaSO₄) solubilities at high NaCl concentrations up to halite saturation have been determined. In the system Na₂SO₄-CaSO₄-H₂O the observed glauberite (Na₂SO₄ · CaSO₄) solubility is higher than that predicted by the high temperature model of Greenberg and Møller (1989), especially at 200°C. At high salt concentrations, solubilities of both anhydrite and glauberite increase with increasing temperature. Stability fields of the minerals syngenite (K₂SO₄ · CaSO₄ · H₂O) and goergeyite (K₂SO₄ · 5 CaSO₄ · H₂O) were determined, and a new phase was found at 200°C in the K₂SO₄-CaSO₄-H₂O system. Chemical and single crystal structure analysis give the formula K₂SO₄ · CaSO₄. The structure is isostructural with palmierite (K₂SO₄ · PbSO₄). The glaserite (“3 K₂SO₄ · Na₂SO₄”) appears as solid solution in the system Na₂SO₄-K₂SO₄-H₂O. Its solubility and stoichiometry was determined as a function of solution composition.     

An isopiestic study of aqueous NaBr and KBr at 50 °C: Chemical equilibrium model of solution behavior and solubility in the NaBr-H2O, KBr-H2O and Na-K-Br-H2O systems to high concentration and temperature

The isopiestic method has been used to determine the osmotic coefficients of the binary solutions NaBr-H₂O (from 0.745 to 5.953 mol kg⁻¹) and KBr-H₂O (from 0.741 to 5.683 mol kg⁻¹) at the temperature t = 50 °C. Sodium chloride solutions have been used as isopiestic reference standards. The isopiestic results obtained have been combined with all other experimental thermodynamic quantities available in literature (osmotic coefficients, water activities, bromide mineral’s solubilities) to construct a chemical model that calculates solute and solvent activities and solid-liquid equilibria in the NaBr-H₂O, KBr-H₂O and Na-K-Br-H₂O systems from dilute to high solution concentration within the 0-300 °C temperature range. The Harvie and Weare [Harvie C., and Weare J. (1980) The prediction of mineral solubilities in naturalwaters: the Na-K-Mg-Ca-Cl-SO₄-H₂O system from zero to high concentration at 25 °C. Geochim. Cosmochim. Acta44, 981-997] solubility modeling approach, incorporating their implementation of the concentration-dependent specific interaction equations of Pitzer [Pitzer K. (1973) Thermodynamics of electrolytes. I. Theoretical basis and general equations. J. Phys. Chem.77, 268-277] is employed. The model for binary systems is validated by comparing activity coefficient predictions with those given in literature, and not used in the parameterization process. Limitations of the mixed solutions model due to data insufficiencies are discussed. This model expands the variable temperature sodium-potassium model of Greenberg and Moller [Greenberg J., and Moller N. (1989) The prediction of mineral solubilities in natural waters: a chemical equilibrium model for the Na-K-Ca-Cl-SO₄-H₂O system to high concentration from 0 to 250 °C. Geochim. Cosmochim. Acta53, 2503-2518] by evaluating Br⁻ pure electrolyte and mixing solution parameters and the chemical potentials of three bromide solid phases: NaBr-2H₂O (cr), NaBr (cr) and KBr (cr).     

Experimental determination and analysis of the solubility of corundum in 0.1-molal CaCl2 solutions between 400 and 600°C at 0.6 to 2.0 kbar

Corundum (α-Al₂O₃) solubility was measured in 0.1-molal CaCl₂ solutions from 400 to 600°C between 0.6 and 2.0 kbar. The Al molality at 2 kbar increases from 3.1 × 10⁻⁴ at 400°C to 12.7 × 10⁻⁴ at 600°C. At 1 kbar, the solubility increases from 1.5 × 10⁻⁴m at 400°C to 3.4 × 10⁻⁴m at 600°C. These molalities are somewhat less than corundum solubility in pure H₂O (Walther, 1997) at 400°C but somewhat greater at 600°C. The distribution of species was computed considering the Al species Al(OH)₃⁰ and Al(OH)₄⁻, consistent with the solubility of corundum in pure H₂O of Walther (1997) and association constants reported in the literature. The calculated solubility was greater than that measured except at 600°C and 2.0 kbar, indicating that neutral-charged species interactions are probably important.A Setchénow model for neutral species resulted in poor fitting of the measured values at 1.0 kbar. This suggests that Al(OH)₃⁰ has a greater stability relative to Al(OH)₄⁻ than given by the models of Pokrovskii and Helgeson (1995) or Diakonov et al. (1996). The significantly lower Al molalities in CaCl₂ relative to those in NaCl solutions at the same concentration confirm the suggestions of Walther (2001) and others that NaAl(OH)₄⁰ rather than an Al-Cl complex must be significant in supercritical NaCl solutions to give the observed increase in corundum solubility with increasing NaCl concentrations.     

Isopiestic measurement of the osmotic and activity coefficients for the NaOH-NaAl(OH)4-H2O system at 313.2 K

By using a specially designed and constructed isopiestic apparatus, we measured the osmotic coefficients at 313.2 K for the NaOH-NaAl(OH)₄-H₂O system with the total alkali molality, m_NaOHT (m_NaOH + m_NaAl[OH]4), from 0.05 mol/kg H₂O to 12 mol/kg H₂O and α_K (m_NaOHT/m_NaAl(OH)4) from 1.64 to 5.53. The mean standard deviation of the measurements is 0.0038. Several sets of the Pitzer model parameters for NaOH-NaAl(OH)₄-H₂O system were then obtained by regressing the measured osmotic coefficients with the Pitzer model and the Pitzer model parameters for NaOH(aq). One set of the results is as follows: β⁽⁰⁾_NaOH: 0.08669, β⁽¹⁾_NaOH: 0.31446, β⁽²⁾_NaOH: −0.00007367, C^Φ_NaOH: 0.003180, β⁽⁰⁾_NaAl(OH)4: 0.03507, β⁽¹⁾_NaAl(OH)4: 0.02401, C^Φ_NaAl(OH)4: −0.001066, θ_{OH⁻Al(OH)4⁻}: 0.08177, Ψ_Na⁺OH⁻_Al(OH)4⁻: −0.01162. The mean standard difference between the calculated and the measured osmotic coefficients is 0.0088. With the obtained Pitzer model parameters, we calculated the values of K = (γ_{NaAl(OH)4,cal²} · m_{Al(OH)4⁻,exp})/(γ_NaOH,cal² · m_OH⁻,exp) for the gibbsite solubility. The results show that the obtained Pitzer model parameters are reliable, and the relative error of the calculated activity coefficients should be < 2.1%. We also compared the calculated gibbsite solubility data among several activity coefficients models over a range of m_NaOHT at various temperatures. The comparison indicates that our activity coefficients model may be approximately applied in the ranges of temperature from 298.2 to 323.2 K and m_NaOHT from 0 to 8 mol/kg H₂O. We also calculated the stoichiometric activity coefficients of NaOH and NaAl(OH)₄ and the activity of H₂O for the NaOH-NaAl(OH)₄-H₂O system, and these calculations establish their variations with m_NaOHT and α_K. These variations imply that the strengths of the repulsive interactions among various anions are in the following sequence: Al(OH)₄⁻-Al(OH)₄⁻ < Al(OH)₄⁻-OH⁻ < OH⁻-OH⁻, and the attractive interaction between Al(OH)₄⁻ and H₂O is weaker than that between OH⁻ and H₂O.     

Wollastonite solubility and free energy of supercritical aqueous CaCl2

The solubility of wollastonite in chloride solutions has been measured between 750° C and 850° C at 2kb pressure and at 800°C at 1 kb pressure. The fugacity of HCl at P and T was controlled by combining Ag+AgCl internally with magnetite +hematite+water as an external hydrogen buffer. Silica activity was fixed by the presence of quartz. The calcium and chloride contents of the solutions were measured after quench. CaCl₂ appears to be the major solute species.Using these compositional data, the Gibbs Free Energy of formation for CaCl₂ was calculated, assuming that the activity coefficient quotient is near unity. The results were combined with quartz solubility data and calculations made by Joesten (1974) to approximate the compositions of CO₂-H₂O rich fluids which prevailed during the metasomatic formation of mono-minerallic zones in calc-silicate nodules from a contact aureole in the Christmas Mountains, Texas. The largest predicted gradients in H₄SiO₄ and CaCl₂ occur across the wollastonite zone.     

Thermodynamic modeling of binary CH4-H2O fluid inclusions

This work reports the application of thermodynamic models, including equations of state, to binary (salt-free) CH₄-H₂O fluid inclusions. A general method is presented to calculate the compositions of CH₄-H₂O inclusions using the phase volume fractions and dissolution temperatures of CH₄ hydrate. To calculate the homogenization pressures and isolines of the CH₄-H₂O inclusions, an improved activity-fugacity model is developed to predict the vapor-liquid phase equilibrium. The phase equilibrium model can predict methane solubility in the liquid phase and water content in the vapor phase from 273 to 623 K and from 1 to 1000 bar (up to 2000 bar for the liquid phase), within or close to experimental uncertainties. Compared to reliable experimental phase equilibrium data, the average deviation of the water content in the vapor phase and methane solubility in the liquid phase is 4.29% and 3.63%, respectively. In the near-critical region, the predicted composition deviations increase to over 10%. The vapor-liquid phase equilibrium model together with the updated volumetric model of homogenous (single-phase) CH₄-H₂O fluid mixtures (Mao S., Duan Z., Hu J. and Zhang D. (2010) A model for single-phase PVTx properties of CO₂-CH₄-C₂H₆-N₂-H₂O-NaCl fluid mixtures from 273 to 1273 K and from 1 to 5000 bar. Chem. Geol.275, 148-160), is applied to calculate the isolines, homogenization pressures, homogenization volumes, and isochores at specified homogenization temperatures and compositions. Online calculation is on the website: http://www.geochem-model.org/.     

Experimental determination of the solubility constant for magnesium chloride hydroxide hydrate (Mg3Cl(OH)5·4H2O, phase 5) at room temperature, and its importance to nuclear waste isolation in geological repositories in salt formations

In this study, the solubility constant of magnesium chloride hydroxide hydrate, Mg₃Cl(OH)₅·4H₂O, termed as phase 5, is determined from a series of solubility experiments in MgCl₂-NaCl solutions. The solubility constant in logarithmic units at 25 °C for the following reaction,

Mg₃Cl(OH)₅·4H₂O+5H⁺=3Mg²⁺+9H₂O(l)+Cl^-     

Solubility of Fe-ettringite (Ca6[Fe(OH)6]2(SO4)3 · 26H2O)

The solubility of Fe-ettringite (Ca₆[Fe(OH)₆]₂(SO₄)₃ · 26H₂O) was measured in a series of precipitation and dissolution experiments at 20 °C and at pH-values between 11.0 and 14.0 using synthesised material. A time-series study showed that equilibrium was reached within 180 days of ageing. After equilibrating, the solid phases were analysed by XRD and TGA while the aqueous solutions were analysed by ICP-OES (calcium, sulphur) and ICP-MS (iron). Fe-ettringite was found to be stable up to pH 13.0. At higher pH-values Fe-monosulphate (Ca₄[Fe(OH)₆]₂(SO₄) · 6H₂O) and Fe-monocarbonate (Ca₄[Fe(OH)₆]₂(CO₃) · 6H₂O) are formed. The solubilities of these hydrates at 25 °C are:      

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.     

Solid-liquid equilibria of Mg(OH)2(cr) and Mg2(OH)3Cl·4H2O(cr) in the system Mg-Na-H-OH-Cl-H2O at 25°C

The solubility of crystalline Mg(OH)₂(cr) was determined by measuring the equilibrium H⁺ concentration in water, 0.01-2.7 m MgCl₂, 0.1-5.6 m NaCl, and in mixtures of 0.5 and 5.0 m NaCl containing 0.01-0.05 m MgCl₂. In MgCl₂ solutions above 2 molal, magnesium hydroxide converted into hydrated magnesium oxychloride. The solid-liquid equilibrium of Mg₂(OH)₃Cl·4H₂O(cr) was studied in 2.1-5.2 m MgCl₂. Using known ion interaction Pitzer coefficients for the system Mg-Na-H-OH-Cl-H₂O (25°C), the following equilibrium constants at I = 0 are calculated:

     

Calcite solubility and speciation in supercritical NaCl-HCl aqueous fluids

The solubility of calcite in NaCl-H₂O and in HCl-H₂O fluids was measured using an extraction-quench hydrothermal apparatus. Experiments were conducted at 2 kbar, between 400° C and 600° C. Measurements in NaCl-H₂O were conducted in two ways: 1) at constant pressure and NaCl concentration, as a function of temperature; and 2) at constant pressure and temperature, as a function of NaCl concentration. In both the NaCl-H₂O and the HCl-H₂O systems, the solubility of calcite increases with increasing chlorine concentrations. For example, the log calcium molality in equilibrium with calcite increases from –3.75 at 2 kbar and 500° C, in pure H₂O to –3.10 at 2 kbar and 500° C at log NaCl molality=–1.67. At fixed pressure and NaCl molality, the solubility of calcite is almost constant from 400° C to 550° C, but increases somewhat at higher temperatures. The results can be used to determine the dominant calcium species in the experimental solutions as a function of NaCl concentration and to obtain values for the second dissociation constant of CaCl_2(aq). At 2 kbar, 400° C, 500° C, and 600° C, we calculate values for the log of the dissociation constant of CaCl⁺ of –2.1, –3.2, and –4.3, respectively. The 400° C and 500° C values are consistent with those obtained by Frantz and Marshall (1982) using electrical conductance techniques. However, our 600° C value is 0.8 log units higher than that reported by Frantz and Marshall. The calcite solubilities in the NaCl-H₂O and HCl-H₂O systems are inconsistent with the solubilities of calcite in pure H₂O reported by Walther and Long (1986). They are, however, consistent with the measurements of calcite solubilities in pure H₂O presented in this study. These results allow for the calculation of the solubilities of calcium silicates and carbonates in fluids that contain CO₂ and NaCl.     

Structure of H2O-saturated peralkaline aluminosilicate melt and coexisting aluminosilicate-saturated aqueous fluid determined in-situ to 800 °C and ∼800 MPa

The structure of H₂O-saturated silicate melts and of silicate-saturated aqueous solutions, as well as that of supercritical silicate-rich aqueous liquids, has been characterized in-situ while the sample was at high temperature (to 800 °C) and pressure (up to 796 MPa). Structural information was obtained with confocal microRaman and with FTIR spectroscopy. Two Al-bearing glasses compositionally along the join Na₂O•4SiO₂-Na₂O•4(NaAl)O₂-H₂O (5 and 10 mol% Al₂O₃, denoted NA5 and NA10) were used as starting materials. Fluids and melts were examined along pressure-temperature trajectories of isochores of H₂O at nominal densities (from PVT properties of pure H₂O) of 0.85 g/cm³ (NA10 experiments) and 0.86 g/cm³ (NA5 experiments) with the aluminosilicate + H₂O sample contained in an externally-heated, Ir-gasketed hydrothermal diamond anvil cell.Molecular H₂O (H₂O°) and OH groups that form bonds with cations exist in all three phases. The OH/H₂O° ratio is positively correlated with temperature and pressure (and, therefore, fugacity of H₂O, f_H2O) with (OH/H₂O°)^melt > (OH/H₂O°)^fluid at all pressures and temperatures. Structural units of Q³, Q², Q¹, and Q⁰ type occur together in fluids, in melts, and, when outside the two-phase melt + fluid boundary, in single-phase liquids. The abundance of Q⁰ and Q¹ increases and Q² and Q³ decrease with f_H2O. Therefore, the NBO/T (nonbridging oxygen per tetrahedrally coordination cations), of melt is a positive function of f_H2O. The NBO/T of silicate in coexisting aqueous fluid, although greater than in melt, is less sensitive to f_H2O.The melt structural data are used to describe relationships between activity of H₂O and melting phase relations of silicate systems at high pressure and temperature. The data were also combined with available partial molar configurational heat capacity of Qⁿ-species in melts to illustrate how these quantities can be employed to estimate relationships between heat capacity of melts and their H₂O content.     

Solubility and solution mechanism of H2O in alkali silicate melts and glasses at high pressure and temperature

The solubility behavior of H₂O in melts in the system Na₂O-SiO₂-H₂O was determined by locating the univariant phase boundary, melt = melt + vapor in the 0.8-2 GPa and 1000°-1300°C pressure and temperature range, respectively. The NBO/Si-range of the melts (0.25-1) was chosen to cover that of most natural magmatic liquids. The H₂O solubility in melts in the system Na₂O-SiO₂-H₂O (X_H2O) ranges between 18 and 45 mol% (O = 1) with (∂X_H2O/∂P)_T∼14-18 mol% H₂O/GPa. The (∂X_H2O/∂P)_T is negatively correlated with NBO/Si (= Na/Si) of the melt. The (∂X_H2O/∂T)_P is in the −0.03 to +0.05 mol% H₂O/°C range, and is negatively correlated with NBO/Si. The [∂X_H2O/∂(NBO/Si)]_P,T is in the −3 to −8 mol% H₂O/(NBO/Si) range. Melts with NBO/Si similar to basaltic liquids (∼0.6-∼1.0) show (∂X_H2O/∂T)_P<0, whereas more polymerized melts exhibit (∂X_H2O/∂T)_P>0. Complete miscibility between hydrous melt and aqueous fluid occurs in the 0.8-2 GPa pressure range for melts with NBO/Si ≤0.5 at T >1100°C. Miscibility occurs at lower pressure the more polymerized the melt.     

Solubility of KFe(CrO4)2·2H2O at 4–75°C

The solubility of KFe(CrO₄)₂·2H₂O, a precipitate recently identified in a Cr(VI)-contaminated soil, was studied in dissolution and precipitation experiments. Ten dissolution experiments were conducted at 4–75°C and initial pH values between 0.8 and 1.2 using synthetic KFe(CrO₄)₂·2H₂O. Four precipitation experiments were conducted at 25°C with final pH values between 0.16 and 1.39. The log K_SP for the reaction

相似文献   

Phase relations in the CH4-H2O-NaCl system at 2 kbar, 300 to 600°C as determined using synthetic fluid inclusions

Synthesis of fluid inclusions in the CH₄-H₂O-NaCl system was accomplished by subjecting fractured quartz or fluorite, along with known quantities of CH₄, H₂O, and NaCl, to a pressure of 2 kbar and temperatures of 300, 400, 500, or 600°C, in sealed Au capsules. Under the elevated P-T conditions, some of the fractures healed, trapping fluids as inclusions. Microthermometric measurements conducted on the fluid inclusions show that at 2 kbar and 400 to 600°C, there are very broad regions of fluid unmixing in the CH₄-H₂O-NaCl system. For those bulk fluid compositions that lie in the two-phase (i.e., immiscible fluids) field, the high-density phase is enriched in NaCl, whereas the low-density phase is enriched in CH₄. For any given bulk composition, the degree of NaCl enrichment in the high-density phase increases, whereas the degree of CH₄ enrichment in the low-density phase decreases, as temperature increases from 400 to 600°C. Our experimental constraints on the size of the two-phase field are generally consistent with results generated using the equation-of-state GEOFLUIDS (available at http://geotherm.ucsd.edu/geofluids/). However, when comparing the compositions of coexisting immiscible fluids, as determined experimentally vs. calculated using GEOFLUIDS, we find that some relatively small but probably significant differences exist between our experiments and this equation of state.     

A high temperature equation of state for the H2O-CaCl2 and H2O-MgCl2 systems

An equation of state (EOS) is developed for salt-water systems in the high temperature range. As an example of the applications, this EOS is parameterized for the calculation of density, immiscibility, and the compositions of coexisting phases in the CaCl₂-H₂O and MgCl₂-H₂O systems from 523 to 973 K and from saturation pressure to 1500 bar. All available volumetric and phase equilibrium measurements of these binaries are well represented by this equation. This EOS is based on a Helmholtz free energy representation constructed from a reference system containing hard-sphere and polar contributions plus an empirical correction. For the temperature and pressure range in this study, the electrolyte solutes are assumed to be associated. The water molecules are modeled as hard spheres with point dipoles and the solute molecules, MgCl₂ and CaCl₂, as hard spheres with point quadrupoles. The free energy of the reference system is calculated from an analytical representation of the Helmholtz free energy of the hard-sphere contributions and perturbative estimates of the electrostatic contributions. The empirical correction used to account for deviations of the reference system predictions from measured data is based on a virial expansion. The formalism allows generalization to aqueous systems containing insoluble gases (CO₂, CH₄), alkali chlorides (NaCl, KCl), and alkaline earth chlorides (CaCl₂, MgCl₂). The program of this model is available as an electronic annex (see EA1 and EA2) and can also be downloaded at: http://www.geochem-model.org/programs.htm.     

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.