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.     

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 loss and nitrogen isotopic fractionation in biotite as a function of metamorphic grade in metapelites from western Maine, USA

Ammonium fixed in micas of metamorphic rocks is a sensitive indicator both of organic-inorganic interactions during diagenesis as well as of the devolatilization history and fluid/rock interaction during metamorphism. In this study, a collection of geochemically well-characterized biotite separates from a series of graphite-bearing Paleozoic greenschist- to upper amphibolite-facies metapelites, western Maine, USA, were analyzed for ammonium nitrogen () contents and isotopic composition (δ¹⁵N_NH4) using the HF-digestion distillation technique followed by the EA-IRMS technique. Biotite separates, sampled from 9 individual metamorphic zones, contain 3000 to 100 ppm with a wide range in δ¹⁵N from +1.6‰ to +9.1‰. Average contents in biotite show a distinct decrease from about 2750 ppm for the lowest metamorphic grade (∼500 °C) down to 218 ppm for the highest metamorphic grade (∼685 °C). Decreasing abundances in are inversely correlated in a linear fashion with increasing K⁺ in biotite as a function of metamorphic grade and are interpreted as a devolatilization effect. Despite expected increasing δ¹⁵N_NH4 values in biotite with nitrogen loss, a significant decrease from the Garnet Zones to the Staurolite Zones was found, followed by an increase to the Sillimanite Zones. This pattern for δ¹⁵N_NH4 values in biotite inversely correlates with Mg/(Mg + Fe) ratios in biotite and is discussed in the framework of isotopic fractionation due to different exchange processes between or , reflecting devolatilization history and redox conditions during metamorphism.     

Oxygen isotope fractionation between apatite-bound carbonate and water determined from controlled experiments with synthetic apatites precipitated at 10-37 °C

The oxygen isotope fractionation between the structural carbonate of inorganically precipitated hydroxyapatite (HAP) and water was determined in the range 10-37 °C. Values of 1000 ln α() are linearly correlated with inverse temperature (K) according to the following equation: 1000 ln α() = 25.19 (±0.53)·T⁻¹ − 56.47 (±1.81) (R² = 0.998). This fractionation equation has a slightly steeper slope than those already established between calcite and water ( [O’Neil et al., 1969] and [Kim and O’Neil, 1997]) even though measured fractionations are of comparable amplitude in the temperature range of these experimental studies. It is consequently observed that the oxygen isotope fractionation between apatite carbonate and phosphate increases from about 7.5‰ up to 9.1‰ with decreasing temperature from 37 °C to 10 °C. A compilation of δ¹⁸O values of both phosphate and carbonate from modern mammal teeth and bones confirms that both variables are linearly correlated, despite a significant scattering up to 3.5‰, with a slope close to 1 and an intercept corresponding to a 1000 ln α() value of 8.1‰. This apparent fractionation factor is slightly higher or close to the fractionation factor expected to be in the range 7-8‰ at the body temperature of mammals.     

Oxygen isotope fractionation factors involving cassiterite (SnO2): II. determination by direct isotope exchange between cassiterite and calcite

Direct oxygen isotope fractionation between cassiterite and calcite has been investigated experimentally at 15 kbar with temperature ranging from 800 to 1000°C. Combined with the quartz-calcite fractionation measured with the same technique (Clayton et al., 1989), the calcite-cassiterite and quartz-cassiterite oxygen isotope fractionations can be expressed as:

     

Experimental studies of oxygen isotope fractionation in the carbonic acid system at 15°, 25°, and 40°C

In light of recent studies that show oxygen isotope fractionation in carbonate minerals to be a function of HCO₃⁻ and CO₃²⁻ concentrations, the oxygen isotope fractionation and exchange between water and components of the carbonic acid system (HCO₃⁻, CO₃²⁻, and CO_2(aq)) were investigated at 15°, 25°, and 40°C. To investigate oxygen isotope exchange between HCO₃⁻, CO₃²⁻, and H₂O, NaHCO₃ solutions were prepared and the pH was adjusted over a range of 2 to 12 by the addition of small amounts of HCl or NaOH. After thermal, chemical, and isotopic equilibrium was attained, BaCl₂ was added to the NaHCO₃ solutions. This resulted in immediate BaCO₃ precipitation; thus, recording the isotopic composition of the dissolved inorganic carbon (DIC). Data from experiments at 15°, 25°, and 40°C (1 atm) show that the oxygen isotope fractionation between HCO₃⁻ and H₂O as a function of temperature is governed by the equation:

     

Oxygen isotopic fractionation during inorganic calcite precipitation ― Effects of temperature, precipitation rate and pH

Stable oxygen isotopic fractionation during inorganic calcite precipitation was experimentally studied by spontaneous precipitation at various pH (8.3 < pH < 10.5), precipitation rates (1.8 < log R < 4.4 μmol m^− 2 h^− 1) and temperatures (5, 25, and 40 °C) using the CO₂ diffusion technique.The results show that the apparent stable oxygen isotopic fractionation factor between calcite and water (α_{calcite–water}) is affected by temperature, the pH of the solution, and the precipitation rate of calcite. Isotopic equilibrium is not maintained during spontaneous precipitation from the solution. Under isotopic non-equilibrium conditions, at a constant temperature and precipitation rate, apparent 1000lnα_{calcite–water} decreases with increasing pH of the solution. If the temperature and pH are held constant, apparent 1000lnα_{calcite–water} values decrease with elevated precipitation rates of calcite. At pH = 8.3, oxygen isotopic fractionation between inorganically precipitated calcite and water as a function of the precipitation rate (R) can be described by the expressions

at 5, 25, and 40 °C, respectively.The impact of precipitation rate on 1000lnα_{calcite–water} value in our experiments clearly indicates a kinetic effect on oxygen isotopic fractionation during calcite precipitation from aqueous solution, even if calcite precipitated slowly from aqueous solution at the given temperature range. Our results support Coplen's work [Coplen T. B. (2007) Calibration of the calcite–water oxygen isotope geothermometer at Devils Hole, Nevada, a natural laboratory. Geochim. Cosmochim. Acta 71, 3948–3957], which indicates that the equilibrium oxygen isotopic fractionation factor might be greater than the commonly accepted value.     

Dissolution rates of talc as a function of solution composition, pH and temperature

Steady-state talc dissolution rates, at far-from-equilibrium conditions, were measured as a function of aqueous silica and magnesium activity, pH from 1 to 10.6, and temperature from 25 to 150 °C. All rates were measured in mixed flow reactors and exhibited stoichiometric or close to stoichiometric dissolution. All measured rates at pH > 2 obtained at a fixed ionic strength of 0.02 M can be described to within experimental uncertainty using

     

Oxalate-promoted forsterite dissolution at low pH

We ran a series of 124 semi-batch reactor experiments to measure the dissolution rate of forsterite in solutions of nitric and oxalic acid solutions over a pH range of 0-7 and total oxalate concentrations between 0 and 0.35 m at 25 °C. We found that the empirical rate law for the dissolution of forsterite in these solutions is

     

An experimental study of Cobalt (II) complexation in Cl and H2S-bearing hydrothermal solutions

The speciation of cobalt (II) in Cl⁻ and H₂S-bearing solutions was investigated spectrophotometrically at temperatures of 200, 250, and 300 °C and a pressure of 100 bars, and by measuring the solubility of cobaltpentlandite at temperatures of 120-300 °C and variable pressures of H₂S. From the results of these experiments, it is evident that CoHS⁺ and predominate in the solutions except at 150 °C, for which the dominant chloride complex is CoCl₃⁻. The logarithms of the stability constant for CoHS⁺ show moderate variation with temperature, decreasing from 6.24 at 120 °C to 5.84 at 200 °C, and increasing to 6.52 at 300 °C. Formation constants for chloride species increase smoothly with temperature and at 300°C their logarithms reach 8.33 for , 6.44 for CoCl₃⁻, 4.94 to 5.36 for , and 2.42 for CoCl⁺. Calculations based on the composition of a model hydrothermal fluid (Ksp-Mu-Qz, KCl = 0.25 m, NaCl = 0.75 m, ΣS = 0.3 m) suggest that at temperatures ?200 °C, cobalt occurs dominantly as CoHS⁺, whereas at higher temperatures the dominant species is .     

The dissociation of carbonic acid in NaCl solutions as a function of concentration and temperature

Potentiometric measurements of the stoichiometric constants for the dissociation of carbonic acid in NaCl solutions ( and ) have been made as a function of molality (0-6 m) and temperature (0-50 °C). The results have been fitted to the equations

     

Molybdic acid ionisation under hydrothermal conditions to 300 °C

This UV spectrophotometric study was aimed at providing precise, experimentally derived thermodynamic data for the ionisation of molybdic acid (H₂MoO₄) from 30 to 300 °C and at equilibrium saturated vapour pressures. The determination of the equilibrium constants and associated thermodynamic parameters were facilitated by spectrophotometric measurements using a specially designed high temperature optical Ti-Pd flow-through cell with silica glass windows.The following van’t Hoff isochore equations describe the temperature dependence of the first and second ionisation constants of molybdic acid up to 300 °C:

     

Arsenous acid ionisation in aqueous solutions from 25 to 300 °C

The ultraviolet spectra of dilute, aqueous arsenic (III)-containing solutions have been measured from 25 to 300 °C at the saturated vapour pressure. From these measurements, the equilibrium constant was obtained for the reaction

     

A new value for the stable oxygen isotope fractionation between dissolved sulfate ion and water

Although the stable oxygen isotope fractionation between dissolved sulfate ion and H₂O (hereafter ) is of physico-chemical and biogeochemical significance, no experimental value has been established until present. The primary reason being that uncatalyzed oxygen exchange between and H₂O is extremely slow, taking ∼10⁵ years at room temperature. For lack of a better approach, values of 16‰ and 31‰ at 25 °C have been assumed in the past, based on theoretical ‘gas-phase’ calculations and extrapolation of laboratory results obtained at temperatures >75 °C that actually pertain to the bisulfate system. Here I use novel quantum-chemistry calculations, which take into account detailed solute-water interactions to establish a new value for of 23‰ at 25 °C. The results of the corresponding calculations for the bisulfate ion are in agreement with observations. The new theoretical values show that sediment -data, which reflect oxygen isotope equilibration between sulfate and ambient water during microbial sulfate reduction, are consistent with the abiotic equilibrium between and water.     

Rare earth element diffusion in natural enstatite

Chemical diffusion coefficients of La, Nd, Eu, Gd, and Yb in natural enstatite have been measured at 850-1250 °C and 1 atm. Anhydrous diffusion experiments were run in Pt capsules in air, or in sealed silica glass capsules under an iron-wüstite (IW) solid buffer. The sources of diffusant were pre-reacted mixtures of synthetic enstatite powder and microcrystalline rare-earth aluminate garnet. Rutherford Backscattering Spectrometry (RBS) was used to measure diffusion profiles. For Gd diffusion in air over the temperature range 1000-1250 °C, the following Arrhenius relation is found for diffusion normal to (210):

     

Copper partitioning in a melt-vapor-brine-magnetite-pyrrhotite assemblage

The effect of sulfur on the partitioning of Cu in a melt-vapor-brine ± magnetite ± pyrrhotite assemblage has been quantified at 800 °C, 140 MPa, f_O2 = nickel-nickel oxide (NNO), logf_S2=-3.0 (i.e., on the magnetite-pyrrhotite curve at NNO), logf_H2S=-1.3 and logf_SO2=-1. All experiments were vapor + brine saturated. Vapor and brine fluid inclusions were trapped in silicate glass and self-healed quartz fractures. Vapor and brine are dominated by NaCl, KCl and HCl in the S-free runs and NaCl, KCl and FeCl₂ in S-bearing runs. Pyrrhotite served as the source of sulfur in S-bearing experiments. The composition of fluid inclusions, glass and crystals were quantified by laser-ablation inductively coupled plasma mass spectrometry. Major element, chlorine and sulfur concentrations in glass were quantified by using electron probe microanalysis. Calculated Nernst-type partition coefficients (±2σ) for Cu between melt-vapor, melt-brine and vapor-brine are , , and , respectively, in the S-free system. The partition coefficients (±2σ) for Cu between melt-vapor, melt-brine and vapor-brine are , , and , respectively, in the S-bearing system. Apparent equilibrium constants (±1σ) describing Cu and Na exchange between vapor and melt and brine and melt were also calculated. The values of are 34 ± 21 and 128 ± 29 in the S-free and S-bearing runs, respectively. The values of are 33 ± 22 and60 ± 5 in the S-free and S-bearing runs, respectively. The data presented here indicate that the presence of sulfur increases the mass transfer of Cu into vapor from silicate melt. Further, the nearly threefold increase in suggests that Cu may be transported as both a chloride and sulfide complex in magmatic vapor, in agreement with hypotheses based on data from natural systems. Most significantly, the data demonstrate that the presence of sulfur enhances the partitioning of Cu from melt into magmatic volatile phases.     

Sr/Ca and Ca/Ca fractionation during inorganic calcite formation: I. Sr incorporation

Partitioning of strontium during spontaneous calcite formation was experimentally studied using an advanced CO₂-diffusion technique. Results at different precipitation rates and T = 5, 25, and 40 °C show that at constant temperature Sr incorporation into calcite is controlled by the precipitation rate (R in μmol/m²/h) according to the individual expressions

     

Complexation of Cu in Hydrothermal NaCl Brines: Ab initio molecular dynamics and energetics

Chloride complexation of Cu⁺ controls the solubility of copper(I) oxide and sulfide ore minerals in hydrothermal and diagenetic fluids. Solubility measurements and optical spectra of high temperature CuCl solutions have been interpreted as indicating the formation of CuCl, , and complexes. However, no other monovalent cation forms tri- and tetrachloro complexes. EXAFS spectra of high temperature Cu-Cl solutions, moreover, appear to show only CuCl and complexes at T > 100 °C. To reconcile these results, I investigated the nature and stability of Cu-Cl complexes using ab initio cluster calculations and ab initio (Car-Parrinello) molecular dynamics simulations for CuCl-NaCl-H₂O systems at 25 to 450 °C. Ab initio molecular dynamic simulations of 1 m CuCl in a 4 m Cl solution give a stable complex at 25 °C over 4 ps but show that the third Cl⁻ is weakly bound. When the temperature is increased along the liquid-vapour saturation curve to 125 °C, the complex dissociates into and Cl⁻; only forms at 325 °C and 1 kbar. Even in a 15.6 m Cl brine at 450 °C, only the complex forms over a 4 ps simulation run.Cluster calculations with a static dielectric continuum solvation field (COSMO) were used in an attempt directly estimate free energies of complex formation in aqueous solution. Consistent with the MD simulations, the complex is slightly stable at 25 °C but decreases in stability with decreasing dielectric constant (ε). The complex is predicted to be unstable at 25 °C and becomes increasingly unstable with decreasing dielectric constant. In hydrothermal fluids (ε < 30) both the and complexes are unstable to dissociation into and Cl⁻.The results obtained here are at odds with recent equations of state that predict and complexes are the predominant species in hydrothermal brines. In contrast, I predict that only complexes will be significant at T > 125 °C, even in NaCl-saturated brines. The high-temperature (T > 125 °C) optical spectra of CuCl solutions and solubility measurements of Cu minerals in Cl-brines need to be reinterpreted in terms of only the CuCl and complexes.     

Kinetic and equilibrium Fe isotope fractionation between aqueous Fe(II) and Fe(III)

Equilibrium and kinetic Fe isotope fractionation between aqueous ferrous and ferric species measured over a range of chloride concentrations (0, 11, 110 mM Cl⁻) and at two temperatures (0 and 22°C) indicate that Fe isotope fractionation is a function of temperature, but independent of chloride contents over the range studied. Using ⁵⁷Fe-enriched tracer experiments the kinetics of isotopic exchange can be fit by a second-order rate equation, or a first-order equation with respect to both ferrous and ferric iron. The exchange is rapid at 22°C, ∼60-80% complete within 5 seconds, whereas at 0°C, exchange rates are about an order of magnitude slower. Isotopic exchange rates vary with chloride contents, where ferrous-ferric isotope exchange rates were ∼25 to 40% slower in the 11 mM HCl solution compared to the 0 mM Cl⁻ (∼10 mM HNO₃) solutions; isotope exchange rates are comparable in the 0 and 110 mM Cl⁻ solutions.The average measured equilibrium isotope fractionations, Δ_{Fe(III)-Fe(II)}, in 0, 11, and 111 mM Cl⁻ solutions at 22°C are identical within experimental error at +2.76±0.09, +2.87±0.22, and +2.76±0.06 ‰, respectively. This is very similar to the value measured by Johnson et al. (2002a) in dilute HCl solutions. At 0°C, the average measured Δ_{Fe(III)-Fe(II)} fractionations are +3.25±0.38, +3.51±0.14 and +3.56±0.16 ‰ for 0, 11, and 111 mM Cl⁻ solutions. Assessment of the effects of partial re-equilibration on isotope fractionation during species separation suggests that the measured isotope fractionations are on average too low by ∼0.20 ‰ and ∼0.13 ‰ for the 22°C and 0°C experiments, respectively. Using corrected fractionation factors, we can define the temperature dependence of the isotope fractionation from 0°C to 22°C as: where the isotopic fractionation is independent of Cl⁻ contents over the range used in these experiments. These results confirm that the Fe(III)-Fe(II) fractionation is approximately half that predicted from spectroscopic data, and suggests that, at least in moderate Cl⁻ contents, the isotopic fractionation is relatively insensitive to Fe-Cl speciation.     

Carbon solubility in core melts in a shallow magma ocean environment and distribution of carbon between the Earth’s core and the mantle

The solubility of carbon in Fe and Fe-5.2 wt.% Ni melts, saturated with graphite, determined by electron microprobe analysis of quenched metal melts was 5.8 ± 0.1 wt.% at 2000 °C, 6.7 ± 0.2 wt.% at 2200 °C, and 7.4 ± 0.2 wt.% at 2410 °C at 2 GPa, conditions relevant for core/mantle differentiation in a shallow magma ocean. These solubilities are slightly lower than low-pressure literature values and significantly beneath calculated values for even higher pressures [e.g., Wood B. J. (1993) Carbon in the core. Earth Planet. Sci. Lett.117, 593-607]. The trend of C solubility versus temperature for Fe-5.2 wt.% Ni melt, within analytical uncertainties, is similar to or slightly lower (∼0.2-0.4 wt.%) than that of pure Fe. Carbon content of core melts and residual mantle silicates derived from equilibrium batch or fractional segregation of core liquids and their comparison with our solubility data and carbon content estimate of the present day mantle, respectively, constrain the partition coefficient of carbon between silicate and metallic melts, in a magma ocean. For the entire range of possible bulk Earth carbon content from chondritic to subchondritic values, of 10⁻⁴ to 1 is derived. But for ∼1000 ppm bulk Earth carbon, is between 10⁻² and 1. Using the complete range of possible for a magma ocean at ∼2200 °C, we predict maximum carbon content of the Earth’s core to be ∼6-7 wt.% and a preferred value of 0.25 ± 0.15 wt.% for a bulk Earth carbon concentration of ∼1000 ppm.