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

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

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

Gold and copper partitioning in magmatic-hydrothermal systems at 800 °C and 100 MPa

Porphyry-type ore deposits sometimes contain fluid inclusion compositions consistent with the partitioning of copper and gold into vapor relative to coexisting brine at the depositional stage. However, this has not been reproduced experimentally at magmatic conditions. In an attempt to determine the conditions under which copper and gold may partition preferentially into vapor relative to brine at temperatures above the solidus of granitic magmas, we performed experiments at 800 °C, 100 MPa, oxygen fugacity () buffered by Ni-NiO, and fixed at either 3.5 × 10⁻² by using intermediate solid solution-pyrrhotite, or 1.2 × 10⁻⁴ by using intermediate solid solution-pyrrhotite-bornite. The coexisting vapor (∼3 wt.% NaCl eq.) and brine (∼68 wt.% NaCl eq.) were composed initially of NaCl + KCl + HCl + H₂O, with starting HCl set to <1000 μg/g in the aqueous mixture. Synthetic vapor and brine fluid inclusions were trapped at run conditions and subsequently analyzed by laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS). Our experiments demonstrate that copper and gold partitioned strongly into the magmatic volatile phase(s) (MVP) (i.e., vapor or brine) relative to a silicate melt over the entire imposed range of . Nernst style partition coefficients between coexisting brine (b) and melt (m), D^b/m (±1σ), range from 3.6(±2.2) × 10¹ to 4(±2) × 10² for copper and from 1.2(±0.6) × 10² to 2.4(±2.4) × 10³ for gold. Partition coefficients between coexisting vapor (v) and melt, D^v/m range from 2.1 ± 0.7 to 18 ± 5 and 7(±3) × 10¹ to 1.6(±1.6) × 10² for copper and gold, respectively. Partition coefficients for all experiments between coexisting brine and vapor, D^b/v (±1σ), range from 7(±2) to 1.0(±0.4) × 10² and 1.7(±0.2) to 15(±2) for copper and gold, respectively. Observed average D^b/v at an of 1.2 × 10⁻⁴ were elevated, 95(±5) and 15 ± 1 for copper and gold, respectively, relative to those at the higher of 3.5 × 10⁻² where D^b/v were 10(±5) for copper and 7(±6) for gold. Thus, there is an inverse relationship between the and the D^b/v for both copper and gold with increasing resulting in a decrease in the D^b/v signifying increased importance of the vapor phase for copper and gold transport. This suggests that copper and gold may complex with volatile S-species as well as Cl-species at magmatic conditions, however, none of the experiments of our study at 800 °C and 100 MPa had a D^b/v ? 1. We did not directly determine speciation, but infer the existence of some metal-sulfur complexes based on the reported data. We suggest that copper and gold partition preferentially into the brine in most instances at or above the wet solidus. However, in most systems, the mass of vapor is greater than the mass of brine, and vapor transport of copper and gold may become more important in the magmatic environment at higher , lower , or near the critical point in a salt-water system. A D^b/v ? 1 at subsolidus hydrothermal conditions may also occur in response to changes in temperature, , , and/or acidity.Additionally, both copper and gold were observed to partition into intermediate solid solution and bornite much more strongly than into vapor, brine or silicate melt. This suggests that, although vapor and brine are both efficient at removing copper and gold from a silicate melt, the presence of Cu-Fe sulfides can sequester a substantial portion of the copper and gold contained within a silicate melt if the Cu-Fe sulfides are abundant.     

Experimental determination of the solubility of natural wollastonite in pure water up to pressures of 5 GPa and at temperatures of 400-800 °C

The solubility of natural, near-end-member wollastonite-I (>99.5% CaSiO₃) has been determined at temperatures from 400 to 800 °C and pressures between 0.8 and 5 GPa in piston-cylinder apparatus with the weight-loss method. Chemical analysis of quench products and optical monitoring in a hydrothermal diamond anvil cell demonstrates that no additional phases form during dissolution. Wollastonite-I, therefore, dissolves congruently in the pressure-temperature range investigated. The solubility of CaSiO₃ varies between 0.175 and 13.485 wt% and increases systematically with both temperature and pressure up to 3.0 GPa. Above 3.0 GPa wollastonite-I reacts rapidly to the high-pressure modification wollastonite-II. No obvious trends are evident in the solubility of wollastonite-II, with values between 1.93 and 10.61 wt%. The systematics of wollastonite-I solubility can be described well by a composite polynomial expression that leads to isothermal linear correlation with the density of water. The molality of dissolved wollastonite-I in pure water is then

log(m_woll)=2.2288-3418.23×T^-1+671386.84×T^-2+logρ_H2O×(5.4578+2359.11×T^-1).     

The solubility of molybdenum dioxide and trioxide in HCl-bearing water vapour at 350 °C and pressures up to 160 bars

The solubility and speciation of the assemblage MoO₂-MoO₃ in water vapour were investigated in experiments conducted at 350 °C, P_total from 59 to 160 bar and fHCl from 0 to 3.4 bar (0-2.0 mol%). Measured solubility at these conditions ranges from 22 to 2500 ppm (∑fMo from 4.4 × 10⁻⁴ to 6.5 × 10⁻² bar). The concentration of Mo in the vapour at fHCl below 0.1 bar is similar to that in pure water vapour, but increases by two orders of magnitude at fHCl above 0.1 bar. The fugacity of gaseous Mo species is independent of chloride concentration at fHCl below 0.1 bar, but increases with increasing fHCl above this pressure. The dominant Mo species at fHCl below 0.1 bar is interpreted to be the same as it is in pure water vapour, and to form as a result of the reaction

(A1)     

Controls on gold solubility in arc magmas: An experimental study at 1000 °C and 4 kbar

In order to (1) explain the worldwide association between epithermal gold-copper-molybdenum deposits and arc magmas and (2) test the hypothesis that adakitic magmas would be Au-specialized, we have determined the solubility of Au at 4 kbar and 1000 °C for three intermediate magmas (two adakites and one calc-alkaline composition) from the Philippines. The experiments were performed over a fO₂ range corresponding to reducing (∼NNO−1), moderately oxidizing (∼NNO+1.5) and strongly oxidizing (∼NNO+3) conditions as measured by solid Ni-Pd-O sensors. They were carried out in gold containers, the latter serving also as the source of gold, in presence of variable amounts of H₂O and, in a few additional experiments, of S. Concentrations of Au in glasses were determined by LA-ICPMS. Gold solubility in melt is very low (30-240 ppb) but increases with fO₂ in a way consistent with the dissolution of gold as both Au¹⁺ and Au³⁺ species. In the S-bearing experiments performed at ∼NNO−1, gold solubility reaches much higher values, from ∼1200 to 4300 ppb, and seems to correlate with melt S content. No systematic difference in gold solubility is observed between the adakitic and the non-adakitic compositions investigated. Oxygen fugacity and the sulfur concentration in melt are the main parameters controlling the incorporation and concentration of gold in magmas. Certain adakitic and non-adakitic magmas have high fO₂ and magmatic S concentrations favorable to the incorporation and transport of gold. Therefore, the cause of a particular association between some arc magmas and Au-Cu-Mo deposits needs to be searched in the origin of those specialized magmas by involvement of Au- and S-rich protoliths. The subducted slab, which contains metal-rich massive sulfides, may constitute a potentially favorable protolith for the genesis of magmas specialized with respect to gold.     

Microbial dissolution of calcite at T = 28 °C and ambient pCO2

This study used batch reactors to quantify the mechanisms and rates of calcite dissolution in the presence and absence of a single heterotrophic bacterial species (Burkholderia fungorum). Experiments were conducted at T = 28°C and ambient pCO₂ over time periods spanning either 21 or 35 days. Bacteria were supplied with minimal growth media containing either glucose or lactate as a C source, NH₄⁺ as an N source, and H₂PO₄⁻ as a P source. Combining stoichiometric equations for microbial growth with an equilibrium mass-balance model of the H₂O-CO₂-CaCO₃ system demonstrates that B. fungorum affected calcite dissolution by modifying pH and alkalinity during utilization of ionic N and C species. Uptake of NH₄⁺ decreased pH and alkalinity, whereas utilization of lactate, a negatively charged organic anion, increased pH and alkalinity. Calcite in biotic glucose-bearing reactors dissolved by simultaneous reaction with H₂CO₃ generated by dissolution of atmospheric CO₂ (H₂CO₃ + CaCO₃ → Ca²⁺ + 2HCO₃⁻) and H⁺ released during NH₄⁺ uptake (H⁺ + CaCO₃ → Ca²⁺ + HCO₃⁻). Reaction with H₂CO₃ and H⁺ supplied ∼45% and 55% of the total Ca²⁺ and ∼60% and 40% of the total HCO₃⁻, respectively. The net rate of microbial calcite dissolution in the presence of glucose and NH₄⁺ was ∼2-fold higher than that observed for abiotic control experiments where calcite dissolved only by reaction with H₂CO₃. In lactate bearing reactors, most H⁺ generated by NH₄⁺ uptake reacted with HCO₃⁻ produced by lactate oxidation to yield CO₂ and H₂O. Hence, calcite in biotic lactate-bearing reactors dissolved by reaction with H₂CO₃ at a net rate equivalent to that calculated for abiotic control experiments. This study suggests that conventional carbonate equilibria models can satisfactorily predict the bulk fluid chemistry resulting from microbe-calcite interactions, provided that the ionic forms and extent of utilization of N and C sources can be constrained. Because the solubility and dissolution rate of calcite inversely correlate with pH, heterotrophic microbial growth in the presence of nonionic organic matter and NH₄⁺ appears to have the greatest potential for enhancing calcite weathering relative to abiotic conditions.     

Coupled phase and aqueous species equilibrium of the H2O-CO2-NaCl-CaCO3 system from 0 to 250 °C, 1 to 1000 bar with NaCl concentrations up to saturation of halite

A model is developed for the calculation of coupled phase and aqueous species equilibrium in the H₂O-CO₂-NaCl-CaCO₃ system from 0 to 250 °C, 1 to 1000 bar with NaCl concentrations up to saturation of halite. The vapor-liquid-solid (calcite, halite) equilibrium together with the chemical equilibrium of H⁺, Na⁺, Ca²⁺, , Ca(OH)⁺, OH⁻, Cl⁻, , , CO_2(aq) and CaCO_3(aq) in the aqueous liquid phase as a function of temperature, pressure, NaCl concentrations, CO_2(aq) concentrations can be calculated, with accuracy close to those of experiments in the stated T-P-m range, hence calcite solubility, CO₂ gas solubility, alkalinity and pH values can be accurately calculated. The merit and advantage of this model is its predictability, the model was generally not constructed by fitting experimental data.One of the focuses of this study is to predict calcite solubility, with accuracy consistent with the works in previous experimental studies. The resulted model reproduces the following: (1) as temperature increases, the calcite solubility decreases. For example, when temperature increases from 273 to 373 K, calcite solubility decreases by about 50%; (2) with the increase of pressure, calcite solubility increases. For example, at 373 K changing pressure from 10 to 500 bar may increase calcite solubility by as much as 30%; (3) dissolved CO₂ can increase calcite solubility substantially; (4) increasing concentration of NaCl up to 2 m will increase calcite solubility, but further increasing NaCl solubility beyond 2 m will decrease its solubility.The functionality of pH value, alkalinity, CO₂ gas solubility, and the concentrations of many aqueous species with temperature, pressure and NaCl_(aq) concentrations can be found from the application of this model. Online calculation is made available on www.geochem-model.org/models/h2o_co2_nacl_caco3/calc.php.     

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.     

Heat capacities of aqueous solutions of sodium hydroxide and water ionization up to 300 °C at 10 MPa

A commercial (Setaram C80) calorimeter has been modified to measure the heat capacities of highly caustic solutions at temperatures up to 300 °C and pressures up to 20 MPa. The improvements have allowed more accurate determination of the isobaric volumetric heat capacities of chemically aggressive liquids at high temperatures. Test measurements with aqueous solutions of sodium chloride showed a reproducibility of about ±0.1%, with an accuracy of ∼0.3% or better, over the whole temperature range. Heat capacities of aqueous solutions of sodium hydroxide at concentrations from 0.5 to 8 mol/kg were measured at temperatures from 50 to 300 °C and a pressure of 10 MPa. Apparent molar isobaric heat capacities of NaOH(aq) were calculated using densities determined previously for the same solutions by vibrating-tube densimetry. Standard state (infinite dilution) partial molar isobaric heat capacities of NaOH(aq) were obtained by extrapolation using an extended Redlich-Meyer equation. Values of the standard heat capacity change for the ionization of water up to 300 °C were derived by combining the present results with the literature data for HCl(aq) and NaCl(aq).     

A thermodynamic model for the prediction of phase equilibria and speciation in the H2O-CO2-NaCl-CaCO3-CaSO4 system from 0 to 250 °C, 1 to 1000 bar with NaCl concentrations up to halite saturation

A thermodynamic model is developed for the calculation of both phase and speciation equilibrium in the H₂O-CO₂-NaCl-CaCO₃-CaSO₄ system from 0 to 250 °C, and from 1 to 1000 bar with NaCl concentrations up to the saturation of halite. The vapor-liquid-solid (calcite, gypsum, anhydrite and halite) equilibrium together with the chemical equilibrium of H⁺,Na⁺,Ca²⁺, , , and CaSO_4(aq) in the aqueous liquid phase as a function of temperature, pressure and salt concentrations can be calculated with accuracy close to the experimental results.Based on this model validated from experimental data, it can be seen that temperature, pressure and salinity all have significant effects on pH, alkalinity and speciations of aqueous solutions and on the solubility of calcite, halite, anhydrite and gypsum. The solubility of anhydrite and gypsum will decrease as temperature increases (e.g. the solubility will decrease by 90% from 360 K to 460 K). The increase of pressure may increase the solubility of sulphate minerals (e.g. gypsum solubility increases by about 20-40% from vapor pressure to 600 bar). Addition of NaCl to the solution may increase mineral solubility up to about 3 molality of NaCl, adding more NaCl beyond that may slightly decrease its solubility. Dissolved CO₂ in solution may decrease the solubility of minerals. The influence of dissolved calcite on the solubility of gypsum and anhydrite can be ignored, but dissolved gypsum or anhydrite has a big influence on the calcite solubility. Online calculation is made available on www.geochem-model.org/model.     

Oxygen and hydrogen isotope fractionation in serpentine-water and talc-water systems from 250 to 450 °C, 50 MPa

Oxygen and hydrogen isotope fractionation factors in the talc-water and serpentine-water systems have been determined by laboratory experiment from 250 to 450 °C at 50 MPa using the partial exchange technique. Talc was synthesized from brucite + quartz, resulting in nearly 100% exchange during reaction at 350 and 450 °C. For serpentine, D-H exchange was much more rapid than ¹⁸O-¹⁶O exchange when natural chrysotile fibers were employed in the initial charge. In experiments with lizardite as the starting charge, recrystallization to chrysotile enhanced the rate of ¹⁸O-¹⁶O exchange with the coexisting aqueous phase. Oxygen isotope fractionation factors in both the talc-water and serpentine-water systems decrease with increasing temperature and can be described from 250 to 450 °C by the relationships: 1000 ln  = 11.70 × 10⁶/T² − 25.49 × 10³/T + 12.48 and 1000 ln  = 3.49 × 10⁶/T² − 9.48 where T is temperature in Kelvin. Over the same temperature interval at 50 MPa, talc-water D-H fractionation is only weakly dependent on temperature, similar to brucite and chlorite, and can be described by the equation: 1000 ln = 10.88 × 10⁶/T² − 41.52 × 10³/T + 5.61 where T is temperature in Kelvin. Our D-H serpentine-water fractionation factors calibrated by experiment decrease with temperature and form a consistent trend with fractionation factors derived from lower temperature field calibrations. By regression of these data, we have refined and extended the D-H fractionation curve from 25 to 450 °C, 50 MPa as follows: 1000 ln  = 3.436 × 10⁶/T² − 34.736 × 10³/T + 21.67 where T is temperature in Kelvin. These new data should improve the application of D-H and ¹⁸O-¹⁶O isotopes to constrain the temperature and origin of hydrothermal fluids responsible for serpentine formation in a variety of geologic settings.     

Solubilities of corundum, wollastonite and quartz in H2O-NaCl solutions at 800 °C and 10 kbar: Interaction of simple minerals with brines at high pressure and temperature

Solubilities of corundum (Al₂O₃) and wollastonite (CaSiO₃) were measured in H₂O-NaCl solutions at 800 °C and 10 kbar and NaCl concentrations up to halite saturation by weight-loss methods. Additional data on quartz solubility at a single NaCl concentration were obtained as a supplement to previous work. Single crystals of synthetic corundum, natural wollastonite or natural quartz were equilibrated with H₂O and NaCl at pressure (P) and temperature (T) in a piston-cylinder apparatus with NaCl pressure medium and graphite heater sleeves. The three minerals show fundamentally different dissolution behavior. Corundum solubility undergoes large enhancement with NaCl concentration, rising rapidly from Al₂O₃ molality (m_Al2O3) of 0.0013(1) (1σ error) in pure H₂O and then leveling off to a maximum of ∼0.015 at halite saturation (X_NaCl ≈ 0.58, where X is mole fraction). Solubility enhancement relative to that in pure H₂O, , passes through a maximum at X_NaCl ≈ 0.15 and then declines towards halite saturation. Quenched fluids have neutral pH at 25 °C. Wollastonite has low solubility in pure H₂O at this P and T(m_CaSiO3=0.0167(6)). It undergoes great enhancement, with a maximum solubility relative to that in H₂O at X_NaCl ≈ 0.33, and solubility >0.5 molal at halite saturation. Solute silica is 2.5 times higher than at quartz saturation in the system H₂O-NaCl-SiO₂, and quenched fluids are very basic (pH 11). Quartz shows monotonically decreasing solubility from m_SiO2=1.248 in pure H₂O to 0.202 at halite saturation. Quenched fluids are pH neutral. A simple ideal-mixing model for quartz-saturated solutions that requires as input only the solubility and speciation of silica in pure H₂O reproduces the data and indicates that hydrogen bonding of molecular H₂O to dissolved silica species is thermodynamically negligible. The maxima in for corundum and wollastonite indicate that the solute products include hydrates and Na⁺ and/or Cl⁻ species produced by molar ratios of reactant H₂O to NaCl of 6:1 and 2:1, respectively. Our results imply that quite simple mechanisms may exist in the dissolution of common rock-forming minerals in saline fluids at high P and T and allow assessment of the interaction of simple, congruently soluble rock-forming minerals with brines associated with deep-crustal metamorphism.     

The system H2O-NaCl. Part I: Correlation formulae for phase relations in temperature-pressure-composition space from 0 to 1000 °C, 0 to 5000 bar, and 0 to 1 XNaCl

Realistic simulations of fluid flow in geologic systems have severely been hampered by the lack of a consistent formulation for fluid properties for binary salt-water fluids over the temperature-pressure-composition ranges encountered in the Earth’s crust. As the first of two companion studies, a set of correlations describing the phase stability relations in the system H₂O-NaCl is developed. Pure water is described by the IAPS-84 equation of state. New correlations comprise the vapor pressure of halite and molten NaCl, the NaCl melting curve, the composition of halite-saturated liquid and vapor, the pressure of vapor + liquid + halite coexistence, the temperature-pressure and temperature-composition relations for the critical curve, and the compositions of liquid and vapor on the vapor + liquid coexistence surface. The correlations yield accurate values for temperatures from 0 to 1000 °C, pressures from 0 to 5000 bar, and compositions from 0 to 1 X_NaCl (mole fraction of NaCl). To facilitate their use in fluid flow simulations, the correlations are entirely formulated as functions of temperature, pressure and composition.     

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.     

Molecular dynamics simulation of the CH4 and CH4-H2O systems up to 10 GPa and 2573 K

Based on our previous study of the intermolecular potential for pure H₂O and the strict evaluation of the competitive potential models for pure CH₄ and the ab initio fitting potential surface across CH₄-H₂O molecules in this study, we carried out more than two thousand molecular dynamics simulations for the PVTx properties of pure CH₄ and the CH₄-H₂O mixtures up to 2573 K and 10 GPa. Comparison of 1941 simulations with experimental PVT data for pure CH₄ shows an average deviation of 0.96% and a maximum deviation of 2.82%. The comparison of the results of 519 simulations of the mixtures with the experimental measurements reveals that the PVTx properties of the CH₄-H₂O mixtures generally agree with the extensive experimental data with an average deviation of 0.83% and 4% in maximum, which is equivalent to the experimental uncertainty. Moreover, the maximum deviation between the experimental data and the simulation results decreases to about 2% as temperature and pressure increase, indicating that the high accuracy of the simulation is well retained in the high temperature and pressure region.After the validation of the simulation method and the intermolecular potential models, we systematically simulated the PVTx properties of this binary system from 673 K and 0.05 GPa to 2573 K and 10 GPa. In order to integrate all the simulation results and the experimental data for the calculation of thermodynamic properties, an equation of state (EOS) is developed for the CH₄-H₂O system covering 673-2573 K and 0.01-10 GPa. Isochores for compositions <4 mol% CH₄ up to 773 K and 600 MPa are also determined in this paper. The program for the EOS can be downloaded from www.geochem-model.org/programs.htm.     

The partitioning of molybdenum(VI) between aqueous liquid and vapour at temperatures up to 370 °C

We have conducted experiments to evaluate the vapour-liquid fractionation of Mo(VI) in the system MoO₃-NH₃-H₂O at 300-370 °C and saturated vapour pressure, using a two-chamber autoclave that allows separate trapping of the vapour and liquid. The measured total Mo concentrations in each phase were used to calculate a distribution coefficient, , which increases as the density of the vapour approaches that of the liquid, and is greater than one for pH ? 4. Molybdenum speciation in the vapour is described by a single complex, MoO₃H₂O. By contrast, thermodynamic modeling of the distribution of Mo species in the liquid indicates that bimolybdate (HMoO₄⁻) is the dominant aqueous species at the conditions of our experiments, and that molybdate (MoO₄²⁻) and molybdic acid (H₂MoO₄⁰) are present in smaller quantities. As vapour-liquid fractionation occurs between neutral species, it is governed by the reaction H₂MoO₄⁰(aq) = MoO₃ · H₂O(g). Fractionation is therefore controlled by the concentration of H₂MoO₄⁰ in the liquid, which increases with increasing temperature and decreasing pH. Owing to the pH dependence of , it cannot be used to describe Mo fractionation in aqueous vapour-liquid systems with compositions different than those of this study. We have therefore calculated a composition-independent (Henry’s Law) constant, , for each experimental point, using the measured total Mo concentration in the vapour and the modeled concentration of H₂MoO₄⁰ in the liquid. This constant may be applied to aqueous vapour-liquid systems of known liquid composition to estimate the concentration of Mo in a vapour for which little chemical information is available, and thereby supplement the available fractionation data for natural porphyry-forming systems. The results of this study demonstrate that at conditions typical of natural porphyry ore-forming systems, a significant amount of molybdenum fractionates into the vapour over the liquid, and the vapour may transport quantities of Mo in excess of that in the liquid at pH conditions below those of the muscovite-microcline reaction boundary.     

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

Zinc complexation in aqueous sulfide solutions: Determination of the stoichiometry and stability of complexes via ZnS(cr) solubility measurements at 100 °C and 150 bars

The solubility of ZnS_(cr) was measured at 100 °C, 150 bars in sulfide solutions as a function of sulfur concentration (m(S_total) = 0.02-0.15) and acidity (pH_t = 2-11). The experiments were conducted using a Ti flow-through hydrothermal reactor enabling the sampling of large volumes of solutions at experimental conditions, with the subsequent concentration and determination of trace quantities of Zn. Prior to the experiments, a long-term in situ conditioning of the solid phase was performed in order to attain the reproducible Zn concentrations (i.e. solubilities). The ZnS_(cr) solubility product was monitored in the course of the experiment. The following species were found to account for Zn speciation in solution: Zn²⁺ (pH_t < 3), (pH_t 3-4.5), (pH_t 5-8), and ZnS(HS)⁻ (pH_t > 8) (pH_t predominance regions are given for m(S_total) = 0.1). Solubility data collected in this study at pH_t > 3 were combined with the ZnS_(cr) solubility product determined at lower pH to yield the following equilibrium constants (t = 100 °C, P = 150 bars):

     

Experimental measurement of the solubility of bismuth phases in water vapor from 220 °C to 300 °C: Implications for ore formation

Preliminary measurements were carried out of the solubility of the O_2-buffering assemblage bismuth + bismite (Bi₂O₃) in aqueous liquid–vapor and vapor-only systems at temperatures of 220, 250 and 300 °C. All experiments were carried out in Ti reaction vessels and were designed such that the Bi solids were contained in a silica tube that prevented contact with liquid water at any time during the experiment. Two blank (no Bi solids present) liquid–vapor experiments at 220 °C yielded Bi concentrations (±1σ) in the condensed liquid of 0.22 ± 0.02 mg/L, whereas the solubility measurements at this temperature yielded an average value of approximately 6 ± 9 mg/L, with replicate experiments ranging from 0.3 to 26 mg/L. Although the 6 mg/L value is associated with a considerable degree of uncertainty, the experiments do indicate transport of Bi through the vapor phase. Measured Bi concentrations in the condensed liquid at 250 °C were in the same range as those at 220 °C, whereas those at 300 °C were significantly lower (i.e., all below the blank value). Vapor-only experiments necessarily contained much smaller initial volumes of water, thereby making the results more susceptible to contamination. Single blank runs at 220 and 300 °C yielded Bi concentrations of 82 and 16 mg/L, respectively. Measured concentrations (±1σ) of Bi in the vapor-only solubility experiments at 220 °C were 235 ± 78 mg/L for an initial water volume of 0.5 mL, and at 300 °C were 56 ± 30 mg/L and 33 ± 21 for initial water volumes of 1 and 2 mL, respectively, suggesting strong preferential partitioning of Bi into the vapor. The results indicate a negative dependence of Bi solubility on temperature, but are inconclusive with respect to the dependence of Bi solubility on water density or fugacity.