A synthetic fluid inclusion study of copper solubility in hydrothermal brines from 525 to 725 °C and 0.3 to 1.7 GPa

Fluid inclusions were synthesized in a piston-cylinder apparatus under mineral-buffered conditions over a range of Cl concentration (0.29 to 11.3 mol kg⁻¹), temperature (525 to 725 °C), and pressure (0.3 to 1.7 GPa). All fluids were buffered by the mineral assemblage native copper + cuprite + talc + quartz. In situ fluid composition was determined by analysing individual fluid inclusions by LA-ICPMS and independently analysing the quench solution. The solubility data provide basic information necessary to model the high temperature behaviour of Cu in magmatic-hydrothermal systems. Copper concentrations up to ∼15 wt% were measured at 630 °C and 0.34 GPa. These results give an upper limit for Cu in natural fluids and support field-based observations of similar high Cu concentrations in fluids at near-magmatic conditions. Experimental evidence indicates that Cu⁺ may form neutral chloride complexes with the general stoichiometry with n up to 4, though n ? 2 is typical for the majority of the experimental conditions. At high pressure (>∼0.5 GPa) there is evidence that hydroxide species, e.g., CuOH⁰, become increasingly important and may predominate over copper(I)-chloride complexes. The roles of fluid mixing, cooling and decompression in ore-forming environments are also discussed.     

A simple predictive model of quartz solubility in water-salt-CO2 systems at temperatures up to 1000 °C and pressures up to 1000 MPa

Knowledge of the solubility of quartz over a broad spectrum of aqueous fluid compositions and T-P conditions is essential to our understanding of water-rock interaction in the Earth’s crust. We propose an equation to compute the molality of aqueous silica, m_SiO2(aq), mol·(kg H₂O)⁻¹, in equilibrium with quartz and water-salt-CO₂ fluids, as follows:

     

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

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.     

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.     

An XAS study of the structure and thermodynamics of Cu(I) chloride complexes in brines up to high temperature (400 °C, 600 bar)

The transport and deposition of copper in saline hydrothermal fluids are controlled by the stability of copper(I) complexes with ligands such as chloride. Despite their role in the formation of most hydrothermal copper deposits, the nature and stability of Cu(I) chloride complexes in highly saline brines remains controversial. We present new X-ray absorption data (P = 600 bar, T = 25-400 °C, salinity up to 17.2 m Cl), which indicate that the linear (x = 1, 2) complexes are stable up to supercritical conditions. Distorted trigonal planar complexes predominate at room temperature and at high salinity (>3 m LiCl): subtle changes in the XANES spectrum with increasing salinity may reflect geometric distortions of this complex. Similar changes were observed in UV-Vis data [Liu, W., Brugger, J., McPhail, D.C., Spiccia, L., 2002. A spectrophotometric study of aqueous copper(I) chloride complexes in LiCl solutions between 100 °C and 250 °C. Geochim. Cosmochim. Acta66, 3615-3633], and were erroneously interpreted as a new species, . Our XAS data and ab-initio XANES calculations show that this tetrahedral species is not present to any significant degree in our solutions. The stability of the complexe decreases with increasing temperature; under supercritical conditions and in brines under magmatic-hydrothermal conditions (e.g., 15.58 m Cl, 400 °C, 600 bar), only the linear Cu(I) chloride complexes were observed. This result and the instability of the complex are also consistent with the recent ab-initio molecular dynamic calculations of Sherman [Sherman D. M.(2007) Complexation of Cu⁺ in hydrothermal NaCl brines: ab-initio molecular dynamics and energetics. Geochim. Cosmochim. Acta71, 714-722]. This study illustrates the power of the quantitative nature of XANES and EXAFS measurements for deciphering the speciation of weak transition metal complexes up to magmatic-hydrothermal conditions.The systematic XANES data are used to retrieve the formation constant for at 150 °C, which is in good agreement with the reinterpretation of the UV-Vis data of Liu et al. (Liu et al., 2002). At high temperatures (?400 °C), the solubility of chalcopyrite in equilibrium with hematite-magnetite-pyrite and K-feldspar-muscovite-quartz calculated with the new properties is lower than that calculated using the previous model, and the calculated solubilities are at the lower end of the range of values measured in brine inclusions from porphyry copper systems.     

An experimental study of the mechanism of the replacement of magnetite by pyrite up to 300 °C

We present the results of an experimental study into the sulfidation of magnetite to form pyrite/marcasite under hydrothermal conditions (90-300 °C, vapor saturated pressures), a process associated with gold deposition in a number of ore deposits. The formation of pyrite/marcasite was studied as a function of reaction time, temperature, pH, sulfide concentration, solid-weight-to-fluid-volume ratio, and geometric surface area of magnetite in polytetrafluoroethylene-lined autoclaves (PTFE) and a titanium and stainless steel flow-through cell. Marcasite was formed only at pH_21°C <4 and was the dominant Fe disulfide at pH_21°C 1.11, while pyrite predominated at pH_21°C >2 and formed even under basic conditions (up to pH_21°C 12-13). Marcasite formation was favored at higher temperatures. Fine-grained pyrrhotite formed in the initial stage of the reaction together with pyrite in some experiments with large surface area of magnetite (grain size <125 μm). This pyrrhotite eventually gave way to pyrite. The transformation rate of magnetite to Fe disulfide increased with decreasing pH (at 120 °C; pH_120°C 0.96-4.42), and that rate of the transformation increased from 120 to 190 °C.Scanning electron microscope (SEM) imaging revealed that micro-pores (0.1-5 μm scale) existed at the reaction front between the parent magnetite and the product pyrite, and that the pyrite and/or marcasite were euhedral at pH_21°C <4 and anhedral at higher pH. The newly formed pyrite was micro-porous (0.1-5 μm); this micro-porosity facilitates fluid transport to the reaction interface between magnetite and pyrite, thus promoting the replacement reaction. The pyrite precipitated onto the parent magnetite was polycrystalline and did not preserve the crystallographic orientation of the magnetite. The pyrite precipitation was also observed on the PTFE liner, which is consistent with pyrite crystallizing from solution. The mechanism of the reaction is that of a dissolution-reprecipitation reaction with the precipitation of pyrite being the rate-limiting step relative to magnetite dissolution under mildly acidic conditions (e.g., pH_155°C 4.42).The experimental results are in good agreement with sulfide phase assemblage and textures reported from sulfidized Banded Iron Formations: pyrite, marcasite and pyrrhotite have been found to exist or co-exist in different sulfidized Banded Iron Formations, and the microtextures show no evidence of sub-μm-scale pseudomorphism of magnetite by pyrite.     

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

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.     

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

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.     

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)     

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.     

Heat capacities of aqueous sodium hydroxide/aluminate mixtures and prediction of the solubility constant of boehmite up to 300 °C

A modified commercial (Setaram C80) calorimeter has been used to measure the isobaric volumetric heat capacities of concentrated alkaline sodium aluminate solutions at ionic strengths from 1 to 6 mol kg⁻¹, with up to 40 mol.% substitution of hydroxide by aluminate, at temperatures from 50 to 300 °C and a pressure of 10 MPa. Apparent molar heat capacities for the mixtures, C_p?, derived from these data were found to depend linearly on the aluminate substitution level, i.e., they followed Young’s rule. These quantities were used to estimate the apparent molar heat capacities of pure, hypothetical sodium aluminate solutions, C_p? (‘NaAl(OH)₄’(aq)). Slopes of the Young’s rule plots were invariant with ionic strength at a given temperature but depended linearly on temperature. The heat capacities of ternary aqueous sodium hydroxide/aluminate mixtures could therefore be modelled using only two parameters in addition to those needed for the correlation of C_p? (NaOH(aq)) reported previously from these laboratories. An assessment of the standard thermodynamic quantities for boehmite, gibbsite and the aluminate ion yielded a set of recommended values that, together with the present heat capacity data, accurately predicts the solubility of gibbsite and boehmite at temperatures up to 300 °C.     

A potentiometric study of the stability of aqueous yttrium-acetate complexes from 25 to 175 °C and 1-1000 bar

The stability of yttrium-acetate (Y-Ac) complexes in aqueous solution was determined potentiometrically at temperatures 25-175 °C (at P_s) and pressures 1-1000 bar (at 25 and 75 °C). Measurements were performed using glass H⁺-selective electrodes in potentiometric cells with a liquid junction. The species YAc²⁺ and were found to dominate yttrium aqueous speciation in experimental solutions at 25-100 °C (log [Ac⁻] < −1.5, pH < 5.2), whereas at 125, 150 and 175 °C introduction of into the Y-Ac speciation model was necessary. The overall stability constants β_n were determined for the reaction

     

Ammonium in aqueous fluids to 600 °C, 1.3 GPa: A spectroscopic study on the effects on fluid properties, silica solubility, and K-feldspar to muscovite reactions

The behavior of ammonium, NH₄⁺, in aqueous systems was studied based on Raman spectroscopic experiments to 600 °C and about 1.3 GPa. Spectra obtained at ambient conditions revealed a strong reduction of the dynamic three-dimensional network of water with addition of ammonium chloride, particularly at small solute concentrations. The differential scattering cross section of the ν₁-NH₄⁺ Raman band in these solutions was found to be similar to that of salammoniac.The Raman band of silica monomers at ∼780 cm⁻¹ was present in all spectra of the fluid at high temperatures in hydrothermal diamond-anvil cell experiments with H₂O ± NH₄Cl and quartz or the assemblage quartz + kyanite + K-feldspar ± muscovite/tobelite. However, these spectra indicated that dissolved silica is less polymerized in ammonium chloride solutions than in comparable experiments with water. Quantification based on the normalized integrated intensity of the H₄SiO₄⁰ band showed that the silica solubility in experiments with H₂O + NH₄Cl was significantly lower than that in equimolal NaCl solutions. This suggests that ammonium causes a stronger decrease in the activity of water in chloridic solutions than sodium.The Raman spectra of the fluid also showed that a significant fraction of ammonium was converted to ammonia, NH₃, in all experiments at temperatures above 300 °C. This indicates a shift towards acidic conditions for experiments without a buffering mineral assemblage. The estimated pH of the fluid was ∼2 at 600 °C, 0.26 GPa, 6.6 m initial NH₄Cl, based on the ratio of the integrated ν₁-NH₃ and ν₁-NH₄⁺ intensities and the HCl⁰ dissociation constant. The NH₃/NH₄⁺ ratio increased with temperature and decreased with pressure. This implies that more ammonium should be retained in K-bearing minerals coexisting with chloridic fluids upon high-P low-T metamorphism. At 500 °C, 0.73 GPa, ammonium partitions preferentially into the fluid, as constrained from infrared spectroscopy on the muscovite and from mass balance.The conversion of K-feldspar to muscovite proceeded much faster in experiments with NH₄Cl solutions than in comparable experiments with water. This is interpreted as being caused by enhancement of the rate-limiting alumina solubility, suggesting complexation of Al with NH₄. Nucleation and growth of mica at the expense of K-feldspar and NH₄⁺/K⁺ exchange between fluid and K-feldspar occurred simultaneously, but incorporation of NH₄⁺ into K-feldspar was distinctly faster than K-feldspar consumption.     

Speciation of reduced C-O-H volatiles in coexisting fluids and silicate melts determined in-situ to ∼1.4 GPa and 800 °C

The structure of silicate melts in the system Na₂O·4SiO₂ saturated with reduced C-O-H volatile components and of coexisting silicate-saturated C-O-H solutions has been determined in a hydrothermal diamond anvil cell (HDAC) by using confocal microRaman and FTIR spectroscopy as structural probes. The experiments were conducted in-situ with the melt and fluid at high temperature (up to 800 °C) and pressure (up to 1435 MPa). Redox conditions in the HDAC were controlled with the reaction, Mo + H₂O = MoO₂ + H₂, which is slightly more reducing than the Fe + H₂O = FeO + H₂ buffer at 800 °C and less.The dominant species in the fluid are CH₄ + H₂O together with minor amounts of molecular H₂ and an undersaturated hydrocarbon species. In coexisting melt, CH₃ - groups linked to the silicate melt structure via Si-O-CH₃ bonding may dominate and possibly coexists with molecular CH₄. The abundance ratio of CH₃ - groups in melts relative to CH₄ in fluids increases from 0.01 to 0.07 between 500 and 800 °C. Carbon-bearing species in melts were not detected at temperatures and pressures below 400 °C and 730 MPa, respectively. A schematic solution mechanism is, Si-O-Si + CH₄?Si-O-CH₃+H-O-Si. This mechanism causes depolymerization of silicate melts. Solution of reduced (C-O-H) components will, therefore, affect melt properties in a manner resembling dissolved H₂O.     

Experimental study and modeling of fluid reaction paths in the quartz-kyanite ± muscovite-water system at 0.7 GPa in the 350-550 °C range: Implications for Al selective transfer during metamorphism

In order to quantify Al transfer in response to fluid-mineral equilibration under evolving metamorphic conditions, isobaric (0.7 GPa) experiments were conducted in the 350-550 °C range. Disequilibrium was induced (1) by holding initially pure water and natural minerals (kyanite + quartz ± muscovite enclosed in a perforated inner capsule) under isothermal conditions and (2) by stepwise temperature variations. In all experiments, secondary Al-bearing phases crystallized in the external tube of a “tube-in-tube” setup (SEM characterization); they are interpreted as witnesses of the evolution of the fluid composition (fluid reaction path). These reaction paths and the subsequent amount of secondary crystallizations were modeled using thermodynamic data from SUPCRT92 and estimates of both starting-mineral dissolution rates and elemental diffusion coefficients from the literature. A major result is that the amount of aluminum transferred to secondary phases is a thousand times larger than the calculated Al concentration in the fluid. Although the crystallization of Al-bearing phases was expected as a response to a temperature decrease, the stepwise temperature increase (20 °C/day) also led to aluminum transfer towards secondary phases. In the course of re-equilibration, the fluid first becomes saturated with respect to aluminosilicates and then reaches silica saturation, due to the low solubility of Al-minerals. Consequently, aluminosilicates partly recrystallize in response to a temperature increase. Crystallization of secondary Al-phases in the external tube implies that aqueous aluminum was efficiently transported from the inner capsule, even in the pure Al₂O₃-SiO₂-H₂O system. Therefore, mass balance calculations considering a constant Al reference frame, i.e., postulating Al immobility, should be regarded with caution.     

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.