Isotopic fractionations associated with phosphoric acid digestion of carbonate minerals: Insights from first-principles theoretical modeling and clumped isotope measurements

Phosphoric acid digestion has been used for oxygen- and carbon-isotope analysis of carbonate minerals since 1950, and was recently established as a method for carbonate ‘clumped isotope’ analysis. The CO₂ recovered from this reaction has an oxygen isotope composition substantially different from reactant carbonate, by an amount that varies with temperature of reaction and carbonate chemistry. Here, we present a theoretical model of the kinetic isotope effects associated with phosphoric acid digestion of carbonates, based on structural arguments that the key step in the reaction is disproportionation of H₂CO₃ reaction intermediary. We test that model against previous experimental constraints on the magnitudes and temperature dependences of these oxygen isotope fractionations, and against new experimental determinations of the fractionation of ¹³C-¹⁸O-containing isotopologues (‘clumped’ isotopic species). Our model predicts that the isotope fractionations associated with phosphoric acid digestion of carbonates at 25 °C are 10.72‰, 0.220‰, 0.137‰, 0.593‰ for, respectively, ¹⁸O/¹⁶O ratios (1000 lnα^∗) and three indices that measure proportions of multiply-substituted isotopologues . We also predict that oxygen isotope fractionations follow the mass dependence exponent, λ of 0.5281 (where ). These predictions compare favorably to independent experimental constraints for phosphoric acid digestion of calcite, including our new data for fractionations of ¹³C-¹⁸O bonds (the measured change in Δ₄₇ = 0.23‰) during phosphoric acid digestion of calcite at 25 °C.We have also attempted to evaluate the effect of carbonate cation compositions on phosphoric acid digestion fractionations using cluster models in which disproportionating H₂CO₃ interacts with adjacent cations. These models underestimate the magnitude of isotope fractionations and so must be regarded as unsucsessful, but do reproduce the general trend of variations and temperature dependences of oxygen isotope acid digestion fractionations among different carbonate minerals. We suggest these results present a useful starting point for future, more sophisticated models of the reacting carbonate/acid interface. Examinations of these theoretical predictions and available experimental data suggest cation radius is the most important factor governing the variations of isotope fractionation among different carbonate minerals. We predict a negative correlation between acid digestion fractionation of oxygen isotopes and of ¹³C-¹⁸O doubly-substituted isotopologues, and use this relationship to estimate the acid digestion fractionation of for different carbonate minerals. Combined with previous theoretical evaluations of ¹³C-¹⁸O clumping effects in carbonate minerals, this enables us to predict the temperature calibration relationship for different carbonate clumped isotope thermometers (witherite, calcite, aragonite, dolomite and magnesite), and to compare these predictions with available experimental determinations. The success of our models in capturing several of the features of isotope fractionation during acid digestion supports our hypothesis that phosphoric acid digestion of carbonate minerals involves disproportionation of transition state structures containing H₂CO₃.     

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.     

An experimental study of oxygen isotope fractionation between inorganically precipitated aragonite and water at low temperatures

To determine oxygen isotope fractionation between aragonite and water, aragonite was slowly precipitated from Ca(HCO₃)₂ solution at 0 to 50°C in the presence of Mg²⁺ or SO₄²⁻. The phase compositions and morphologies of synthetic minerals were detected by X-ray diffraction (XRD) and scanning electron microscopy (SEM) techniques. The effects of aragonite precipitation rate and excess dissolved CO₂ gas in the initial Ca(HCO₃)₂ solution on oxygen isotope fractionation between aragonite and water were investigated. For the CaCO₃ minerals slowly precipitated by the CaCO₃ or NaHCO₃ dissolution method at 0 to 50°C, the XRD and SEM analyses show that the rate of aragonite precipitation increased with temperature. Correspondingly, oxygen isotope fractionations between aragonite and water deviated progressively farther from equilibrium. Additionally, an excess of dissolved CO₂ gas in the initial Ca(HCO₃)₂ solution results in an increase in apparent oxygen isotope fractionations. As a consequence, the experimentally determined oxygen isotope fractionations at 50°C indicate disequilibrium, whereas the relatively lower fractionation values obtained at 0 and 25°C from the solution with less dissolved CO₂ gas and low precipitation rates indicate a closer approach to equilibrium. Combining the lower values at 0 and 25°C with previous data derived from a two-step overgrowth technique at 50 and 70°C, a fractionation equation for the aragonite-water system at 0 to 70°C is obtained as follows:

     

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

Experimental determination of oxygen isotope fractionations between CO2 vapor and soda-melilite melt

We report results of experiments constraining oxygen isotope fractionations between CO₂ vapor and Na-rich melilitic melt at 1 bar and 1250 and 1400°C. The fractionation factor constrained by bracketed experiments, 1000^.lnα_{CO2-Na melilitic melt}, is 2.65±0.25 ‰ (±2σ; n=92) at 1250°C and 2.16±0.16 ‰ (2σ; n=16) at 1400°C. These values are independent of Na content over the range investigated (7.5 to 13.0 wt. % Na₂O). We combine these data with the known reduced partition function ratio of CO₂ to obtain an equation describing the reduced partition function ratio of Na-rich melilite melt as a function of temperature. We also fit previously measured CO₂-melt or -glass fractionations to obtain temperature-dependent reduced partition function ratios for all experimentally studied melts and glasses (including silica, rhyolite, albite, anorthite, Na-rich melilite, and basalt). The systematics of these data suggest that reduced partition function ratios of silicate melts can be approximated either by using the Garlick index (a measure of the polymerization of the melt) or by describing melts as mixtures of normative minerals or equivalent melt compositions. These systematics suggest oxygen isotope fractionation between basalt and olivine at 1300°C of approximately 0.4 to 0.5‰, consistent with most (but not all) basalt glass-olivine fractionations measured in terrestrial and lunar basalts.     

Oxygen isotope fractionation between synthetic aragonite and water: Influence of temperature and Mg concentration

Aragonite was precipitated in the laboratory at 0, 5, 10, 25, and 40 °C to determine the temperature dependence of the equilibrium oxygen isotope fractionation between aragonite and water. Forced CO₂ degassing, passive CO₂ degassing, and constant addition methods were employed to precipitate aragonite from supersaturated solutions, but the resulting aragonite-water oxygen isotope fractionation was independent of the precipitation method. In addition, under the experimental conditions of this study, the effect of precipitation rate on the oxygen isotope fractionation between aragonite and water was almost within the analytical error of ±∼0.13‰ and thus insignificant. Because the presence of Mg²⁺ ions is required to nucleate and precipitate aragonite from Na-Ca-Cl-HCO₃ solutions under these experimental conditions, the influence of the total Mg²⁺ concentration (up to ∼0.9 molal) on the aragonite-water oxygen isotope fractionation was examined at 25 °C. No significant Mg²⁺ ion effect, or oxygen isotope salt effect, was detected up to 100 mmolal total Mg²⁺ but a noticeable isotope salt effect was observed at ∼0.9 molal total Mg²⁺.On the basis of results of the laboratory synthesis experiments, a new expression for the aragonite-water fractionation is proposed over the temperature range of 0-40 °C:

1000lnα_{aragonite-water}=17.88±0.13(10³/T)-31.14±0.46     

Silicon isotopic fractionation during adsorption of aqueous monosilicic acid onto iron oxide

The quantification of silicon isotopic fractionation by biotic and abiotic processes contributes to the understanding of the Si continental cycle. In soils, light Si isotopes are selectively taken up by plants, and concentrate in secondary clay-sized minerals. Si can readily be retrieved from soil solution through the specific adsorption of monosilicic acid () by iron oxides. Here, we report on the Si-isotopic fractionation during adsorption on synthesized ferrihydrite and goethite in batch experiment series designed as function of time (0-504 h) and initial concentration (ic) of Si in solution (0.21-1.80 mM), at 20 °C, constant pH (5.5) and ionic strength (1 mM). At various contact times, the δ²⁹Si vs. NBS28 compositions were determined in selected solutions (ic = 0.64 and 1.06 mM Si) by MC-ICP-MS in dry plasma mode with external Mg doping with an average precision of ±0.08‰ (±2σ_SEM). Per oxide mass, ferrihydrite (74-86% of initial Si loading) adsorbed more Si than goethite (37-69%) after 504 h of contact over the range of initial Si concentration 0.42-1.80 mM. Measured against its initial composition (δ²⁹Si = +0.01 ± 0.04‰ (±2σ_SD)), the remaining solution was systematically enriched in ²⁹Si, reaching maximum δ²⁹Si values of +0.70 ± 0.07‰ for ferrihydrite and +0.50 ± 0.08‰ for goethite for ic 1.06 mM. The progressive ²⁹Si enrichment of the solution fitted better a Rayleigh distillation path than a steady state model. The fractionation factor ²⁹ε (±1σ_SD) was estimated at −0.54 ± 0.03‰ for ferrihydrite and −0.81 ± 0.12‰ for goethite. Our data imply that the sorption of onto synthetic iron oxides produced a distinct Si-isotopic fractionation for the two types of oxide but in the same order than that generated by Si uptake by plants and diatoms. They further suggest that the concentration of light Si isotopes in the clay fraction of soils is partly due to sorption onto secondary clay-sized iron oxides.     

Controls on ostracod valve geochemistry: Part 2. Carbon and oxygen isotope compositions

The stable carbon and oxygen isotope compositions of fossil ostracods are powerful tools to estimate past environmental and climatic conditions. The basis for such interpretations is that the calcite of the valves reflects the isotopic composition of water and its temperature of formation. However, calcite of ostracods is known not to form in isotopic equilibrium with water and different species may have different offsets from inorganic precipitates of calcite formed under the same conditions. To estimate the fractionation during ostracod valve calcification, the oxygen and carbon isotope compositions of 15 species living in Lake Geneva were related to their autoecology and the environmental parameters measured during their growth. The results indicate that: (1) Oxygen isotope fractionation is similar for all species of Candoninae with an enrichment in ¹⁸O of more than 3‰ relative to equilibrium values for inorganic calcite. Oxygen isotope fractionation for Cytheroidea is less discriminative relative to the heavy oxygen, with enrichments in ¹⁸O for these species of 1.7 to 2.3‰. Oxygen isotope fractionations for Cyprididae are in-between those of Candoninae and Cytheroidea. The difference in oxygen isotope fractionation between ostracods and inorganic calcite has been interpreted as resulting from a vital effect. (2) Comparison with previous work suggests that oxygen isotope fractionation may depend on the total and relative ion content of water. (3) Carbon isotope compositions of ostracod valves are generally in equilibrium with DIC. The specimens’ δ¹³C values are mainly controlled by seasonal variations in δ¹³C_DIC of bottom water or variation thereof in sediment pore water. (4) Incomplete valve calcification has an effect on carbon and oxygen isotope compositions of ostracod valves. Preferential incorporation of at the beginning of valve calcification may explain this effect. (5) Results presented here as well as results from synthetic carbonate growth indicate that different growth rates or low pH within the calcification site cannot be the cause of oxygen isotope ‘vital effects’ in ostracods. Two mechanisms that might enrich the ¹⁸O of ostracod valves are deprotonation of that may also contribute to valve calcification, and effects comparable to salt effects with high concentrations of Ca and/or Mg within the calcification site that may also cause a higher temperature dependency of oxygen isotope fractionation.     

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.     

The origin of Zn isotope fractionation in sulfides

Isotope fractionation of Zn between aqueous sulfide, chloride, and carbonate species (Zn²⁺, Zn(HS)₂, , , ZnS(HS)⁻, ZnCl⁺, ZnCl₂, , and ZnCO₃) was investigated using ab initio methods. Only little fractionation is found between the sulfide species, whereas carbonates are up to 1‰ heavier than the parent solution. At pH > 3 and under atmospheric-like CO₂ pressures, isotope fractionation of Zn sulfides precipitated from sulfidic solutions is affected by aqueous sulfide species and the δ⁶⁶Zn of sulfides reflect these in the parent solutions. Under high P_CO2 conditions, carbonate species become abundant. In high P_CO2 conditions of hydrothermal solutions, Zn precipitated as sulfides is isotopically nearly unfractionated with respect to a low-pH parent fluid. In contrast, negative δ⁶⁶Zn down to at least −0.6‰ can be expected in sulfides precipitated from solutions with pH > 9. Zinc isotopes in sulfides and rocks therefore represent a potential indicator of mid to high pH in ancient hydrothermal fluids.     

Iron isotope fractionation and atom exchange during sorption of ferrous iron to mineral surfaces

The application of stable Fe isotopes as a tracer of the biogeochemical Fe cycle necessitates a mechanistic knowledge of natural fractionation processes. We studied the equilibrium Fe isotope fractionation upon sorption of Fe(II) to aluminum oxide (γ-Al₂O₃), goethite (α-FeOOH), quartz (α-SiO₂), and goethite-loaded quartz in batch experiments, and performed continuous-flow column experiments to study the extent of equilibrium and kinetic Fe isotope fractionation during reactive transport of Fe(II) through pure and goethite-loaded quartz sand. In addition, batch and column experiments were used to quantify the coupled electron transfer-atom exchange between dissolved Fe(II) (Fe(II)_aq) and structural Fe(III) of goethite. All experiments were conducted under strictly anoxic conditions at pH 7.2 in 20 mM MOPS (3-(N-morpholino)-propanesulfonic acid) buffer and 23 °C. Iron isotope ratios were measured by high-resolution MC-ICP-MS. Isotope data were analyzed with isotope fractionation models. In batch systems, we observed significant Fe isotope fractionation upon equilibrium sorption of Fe(II) to all sorbents tested, except for aluminum oxide. The equilibrium enrichment factor, , of the Fe(II)_sorb-Fe(II)_aq couple was 0.85 ± 0.10‰ (±2σ) for quartz and 0.85 ± 0.08‰ (±2σ) for goethite-loaded quartz. In the goethite system, the sorption-induced isotope fractionation was superimposed by atom exchange, leading to a δ^56/54Fe shift in solution towards the isotopic composition of the goethite. Without consideration of atom exchange, the equilibrium enrichment factor was 2.01 ± 0.08‰ (±2σ), but decreased to 0.73 ± 0.24‰ (±2σ) when atom exchange was taken into account. The amount of structural Fe in goethite that equilibrated isotopically with Fe(II)_aq via atom exchange was equivalent to one atomic Fe layer of the mineral surface (∼3% of goethite-Fe). Column experiments showed significant Fe isotope fractionation with δ^56/54Fe(II)_aq spanning a range of 1.00‰ and 1.65‰ for pure and goethite-loaded quartz, respectively. Reactive transport of Fe(II) under non-steady state conditions led to complex, non-monotonous Fe isotope trends that could be explained by a combination of kinetic and equilibrium isotope enrichment factors. Our results demonstrate that in abiotic anoxic systems with near-neutral pH, sorption of Fe(II) to mineral surfaces, even to supposedly non-reactive minerals such as quartz, induces significant Fe isotope fractionation. Therefore we expect Fe isotope signatures in natural systems with changing concentration gradients of Fe(II)_aq to be affected by sorption.     

Experimental studies of oxygen isotope fractionation between rhodochrosite (MnCO3) and water at low temperatures

Rhodochrosite crystals were precipitated from Na-Mn-Cl-HCO₃ parent solutions following passive, forced and combined passive-to-forced CO₂ degassing methods. Forced and combined passive-to-forced CO₂ degassing produced rhodochrosite crystals with a small non-equilibrium oxygen isotope effect whereas passive CO₂ degassing protocols yielded rhodochrosite in apparent isotopic equilibrium with water. On the basis of the apparent equilibrium isotopic data, a new temperature-dependent relation is proposed for the oxygen isotope fractionation between rhodochrosite and water between 10 and 40 °C:

1000lnα_{rhodochrosite-water}=17.84±0.18(10³/T)-30.24±0.62     

Abundance of mass 47 CO2 in urban air, car exhaust, and human breath

Atmospheric carbon dioxide is widely studied using records of CO₂ mixing ratio, δ¹³C and δ¹⁸O. However, the number and variability of sources and sinks prevents these alone from uniquely defining the budget. Carbon dioxide having a mass of 47 u (principally ¹³C¹⁸O¹⁶O) provides an additional constraint. In particular, the mass 47 anomaly (Δ₄₇) can distinguish between CO₂ produced by high temperature combustion processes vs. low temperature respiratory processes. Δ₄₇ is defined as the abundance of mass 47 isotopologues in excess of that expected for a random distribution of isotopes, where random distribution means that the abundance of an isotopologue is the product of abundances of the isotopes it is composed of and is calculated based on the measured ¹³C and ¹⁸O values. In this study, we estimate the δ¹³C (vs. VPDB), δ¹⁸O (vs. VSMOW), δ47, and Δ₄₇ values of CO₂ from car exhaust and from human breath, by constructing ‘Keeling plots’ using samples that are mixtures of ambient air and CO₂ from these sources. δ47 is defined as , where is the R⁴⁷ value for a hypothetical CO₂ whose δ¹³C_VPDB = 0, δ¹⁸O_VSMOW = 0, and Δ₄₇ = 0. Ambient air in Pasadena, CA, where this study was conducted, varied in [CO₂] from 383 to 404 μmol mol⁻¹, in δ¹³C and δ¹⁸O from −9.2 to −10.2‰ and from 40.6 to 41.9‰, respectively, in δ47 from 32.5 to 33.9‰, and in Δ₄₇ from 0.73 to 0.96‰. Air sampled at varying distances from a car exhaust pipe was enriched in a combustion source having a composition, as determined by a ‘Keeling plot’ intercept, of −24.4 ± 0.2‰ for δ¹³C (similar to the δ¹³C of local gasoline), δ¹⁸O of 29.9 ± 0.4‰, δ47 of 6.6 ± 0.6‰, and Δ₄₇ of 0.41 ± 0.03‰. Both δ¹⁸O and Δ₄₇ values of the car exhaust end-member are consistent with that expected for thermodynamic equilibrium at∼200 °C between CO₂ and water generated by combustion of gasoline-air mixtures. Samples of CO₂ from human breath were found to have δ¹³C and δ¹⁸O values broadly similar to those of car exhaust-air mixtures, −22.3 ± 0.2 and 34.3 ± 0.3‰, respectively, and δ47 of 13.4 ± 0.4‰. Δ₄₇ in human breath was 0.76  ± 0.03‰, similar to that of ambient Pasadena air and higher than that of the car exhaust signature.     

Carbon isotopic fractionation between CO2 vapour,silicate and carbonate melts: an experimental study to 30 kbar

The carbon isotopic fractionation between CO₂ vapour and sodamelilite (NaCaAlSi₂O₇) melt over a range of pressures and temperatures has been investigated using solid-media piston-cylinder high pressure apparatus. Ag₂C₂O₄ was the source of CO₂ and experimental oxygen fugacity was buffered at hematite-magnetite by the double capsule technique. The abundance and isotopic composition of carbon dissolved in sodamelilite (SM) glass were determined by stepped heating and the ¹³C of coexisting vapour was determined directly by capsule piercing. CO₂ solubility in SM displays a complex behavior with temperature. At pressures up to 10 kbars CO₂ dissolves in SM to form carbonate ion complexes and the solubility data suggest slight negative temperature dependence. Above 20 kbars CO₂ reacts with SM to form immiscible Na-rich silicate and Ca-rich carbonate melts and CO₂ solubility in Na-enriched silicate melt rises with increasing temperature above the liquidus. Measured values for carbon isotopic fractionation between CO₂ vapour and carbonate ions dissoived in sodamelilite melt at 1200°–1400° C and 5–30 kbars average 2.4±0.2, favouring¹³C enrichment in CO₂ vapour. The results are maxima and are independent of pressure and temperature. Similar values of 2 are obtained for the carbon isotopic fractionation between CO₂ vapour and carbonate melts at 1300°–1400° C and 20–30 kbars.     

An experimental determination of chlorine isotope fractionation in acid systems and applications to volcanic fumaroles

The δ³⁷Cl values of volcanic fumarole gases and bubbling springs were measured from the Central American and the Kurile arcs. Low temperature gas samples from the Central American arc have δ³⁷Cl values generally between −2 and 2‰, whereas high-temperature fumaroles (>100 °C) range from 4 to 12‰, with several outliers. This is in contrast to the high-temperature fumaroles from the Kurile island Kudryavy which have slightly positive δ³⁷Cl values, averaging 0.8‰ (±0.6, 1σ), and from our previous work on Izu and Mariana arc samples in which the δ³⁷Cl values of fumarole and ash samples are similar to each other and negative. Assuming that the source for the high-T Central American fumaroles has typical subduction δ³⁷Cl values (−2.5 to 1‰), then there must be a large Cl isotope fractionation in the near-surface fumarolic system. The most likely fractionation mechanism for the high δ³⁷Cl values is between Cl⁻_aq − HCl(g), but published theoretical fractionation for this pair is only ∼1.5‰, insufficient to explain the large range of values observed in the fumaroles. Three experiments were undertaken in order to identify a process that could cause the wide range of δ³⁷Cl values observed in the high-temperature fumaroles. Results are the following: (1) A sub-boiling equilibration experiment between aqueous chloride and HCl gas had , in agreement with the theoretical calculations. (2) Evaporation of HCl(g) from hydrochloric acid at room temperature had fractionation in the opposite sense, with a . (3) A ‘synthetic fumarole’ gave large positive fractionations up to 9‰, with ³⁷Cl strongly partitioned into the vapor phase. The ‘fumarole’ experiments were made by bubbling dry air through boiling hydrochloric acid in an Erlenmeyer flask, and collecting the evolved HCl(g) in a second ‘downstream’ flask filled with distilled water. This extreme enrichment is likely due to a distillation process in which ³⁷Cl-enriched HCl(g) is stripped from the hydrochloric acid followed by a significant fraction of the light HCl(g) redissolving into the constantly condensing water vapor on the walls of the first flask. This distillation experiment creates a Cl isotope fractionation that is the same order of magnitude as observed in the high-temperature fumaroles in Central America. These results suggest that there must be a H₂O liquid-vapor region in the sub-surface fumarole conduit where light Cl is stripped from the HCl gas as it passes through the fumarole. Similar ³⁷Cl enrichments are expected in fossil epithermal boiling systems.     

Oxygen and sulfur isotope systematics of sulfate produced by bacterial and abiotic oxidation of pyrite

To better understand reaction pathways of pyrite oxidation and biogeochemical controls on δ¹⁸O and δ³⁴S values of the generated sulfate in acid mine drainage (AMD) and other natural environments, we conducted a series of pyrite oxidation experiments in the laboratory. Our biological and abiotic experiments were conducted under aerobic conditions by using O₂ as an oxidizing agent and under anaerobic conditions by using dissolved Fe(III)_aq as an oxidant with varying δ¹⁸O_H2O values in the presence and absence of Acidithiobacillus ferrooxidans. In addition, aerobic biological experiments were designed as short- and long-term experiments where the final pH was controlled at ∼2.7 and 2.2, respectively. Due to the slower kinetics of abiotic sulfide oxidation, the aerobic abiotic experiments were only conducted as long term with a final pH of ∼2.7. The δ³⁴S_SO4 values from both the biological and abiotic anaerobic experiments indicated a small but significant sulfur isotope fractionation (∼−0.7‰) in contrast to no significant fractionation observed from any of the aerobic experiments. Relative percentages of the incorporation of water-derived oxygen and dissolved oxygen (O₂) to sulfate were estimated, in addition to the oxygen isotope fractionation between sulfate and water, and dissolved oxygen. As expected, during the biological and abiotic anaerobic experiments all of the sulfate oxygen was derived from water. The percentage incorporation of water-derived oxygen into sulfate during the oxidation experiments by O₂ varied with longer incubation and lower pH, but not due to the presence or absence of bacteria. These percentages were estimated as 85%, 92% and 87% from the short-term biological, long-term biological and abiotic control experiments, respectively. An oxygen isotope fractionation effect between sulfate and water (ε¹⁸O_SO4-H2O) of ∼3.5‰ was determined for the anaerobic (biological and abiotic) experiments. This measured value was then used to estimate the oxygen isotope fractionation effects between sulfate and dissolved oxygen in the aerobic experiments which were −10.0‰, −10.8‰, and −9.8‰ for the short-term biological, long-term biological and abiotic control experiments, respectively. Based on the similarity between δ¹⁸O_SO4 values in the biological and abiotic experiments, it is suggested that δ¹⁸O_SO4 values cannot be used to distinguish biological and abiotic mechanisms of pyrite oxidation. The results presented here suggest that Fe(III)_aq is the primary oxidant for pyrite at pH < 3, even in the presence of dissolved oxygen, and that the main oxygen source of sulfate is water-oxygen under both aerobic and anaerobic conditions.     

Experimental investigation of the effects of temperature and ionic strength on Mo isotope fractionation during adsorption to manganese oxides

The Mo stable isotope system is being applied to study changes in ocean redox. Such applications implicitly assume that Mo isotope fractionation in aqueous systems is relatively insensitive to frequently changing environmental variables such as temperature (T) and ionic strength (I). A major driver of fractionation is the adsorption of Mo to Mn oxyhydroxide surfaces [Barling J. and Anbar A. D. (2004) Molybdenum isotope fractionation during adsorption by manganese oxides. Earth Planet. Sci. Lett.217(3-4), 315-329]. Here, we report the results of experiments that determine the extent to which Mo isotope fractionation during adsorption of Mo to the Mn oxyhydroxide mineral birnessite is sensitive to T and I. The results are compared to new predictions from quantum chemical computations. We measured fractionation from 1 to 50 °C at I = 0.1 m and found that Δ^97/95Mo_{dissolved-adsorbed} varies from 1.9‰ to 1.6‰ over this temperature range. Experiments were also performed at 25 °C in synthetic seawater (I = 0.7); fractionation at this condition was the same within analytical error as in low ionic strength experiments. These findings confirm that the Mo isotope fractionation during adsorption to Mn oxyhydroxides is relatively insensitive to variations and T and I over environmentally relevant ranges. To relate these findings to potential mechanisms of Mo isotope fractionation, we also report results for density functional theory computations of the fractionation between and various possible structures of molybdic acid as a function of temperature. Because no plausible species fractionates from with a magnitude matching the experiments, we are left with three possibilities to explain the fractionation: (1) solvation effects on the vibrational frequencies of aqueous species considered thus far are significant, such that our calculations in vacuo yield inaccurate fractionations; (2) a trace aqueous species not yet considered fractionates from and then adsorbs to birnessite; or (3) a surface complex not present in solution forms on birnessite in which Mo is not tetrahedrally coordinated. Our findings help validate assumptions underlying paleoceanographic applications of the Mo isotope system and also lead us closer to understanding the mechanism of isotope fractionation during adsorption of Mo to Mn oxyhydroxides.     

Sulfur isotope signatures in gypsiferous sediments of the Estancia and Tularosa Basins as indicators of sulfate sources, hydrological processes, and microbial activity

In order to reconstruct paleo-environmental conditions for the saline playa lakes of the Rio Grande Rift, we investigated sediment sulfate sources using sulfur isotope compositions of dissolved ions in modern surface water, groundwater, and precipitated in the form of gypsum sediments deposited during the Pleistocene and Holocene in the Tularosa and Estancia Basins. The major sulfate sources are Lower and Middle Permian marine evaporites (δ³⁴S of 10.9-14.4‰), but the diverse physiography of the Tularosa Basin led to a complex drainage system which contributed sulfates from various sources depending on the climate at the time of sedimentation. As inferred from sulfur isotope mass balance constraints, weathering of sulfides of magmatic/hydrothermal and sedimentary origin associated with climate oscillations during Last Glacial Maximum contributed about 35-50% of the sulfates and led to deposition of gypsum with δ³⁴S values of −1.2‰ to 2.2‰ which are substantially lower than Permian evaporates. In the Estancia Basin, microbial sulfate reduction appears to overprint sulfur isotopic signatures that might elucidate past groundwater flows. A Rayleigh distillation model indicates that about 3-18% of sulfates from an inorganic groundwater pool (δ³⁴S of 12.6-13.8‰) have been metabolized by bacteria and preserved as partially to fully reduced sulfur-bearing minerals species (elemental sulfur, monosulfides, disulfides) with distinctly negative δ³⁴S values (−42.3‰ to −20.3‰) compared to co-existing gypsum (−3.8‰ to 22.4‰). For the Tularosa Basin microbial sulfate reduction had negligible effect on δ³⁴S value of the gypsiferous sediments most likely because of higher annual temperatures (15-33 °C) and lower organic carbon content (median 0.09%) in those sediments leading to more efficient oxidation of H₂S and/or smaller rates of sulfate reduction compared to the saline playas of the Estancia Basin (5-28 °C; median 0.46% of organic carbon).The White Sands region of the Tularosa Basin is frequently posited as a hydrothermal analogue for Mars. High temperatures of groundwater (33.3 °C) and high δ¹⁸O(H₂O) values (1.1‰) in White Sands, however, are controlled predominantly by seasonal evaporation rather than the modern influx of hydrothermal fluids. Nevertheless, it is possible that some of the geochemical processes in White Sands, such as sulfide weathering during climate oscillations and upwelling of highly mineralized waters, might be considered as valid terrestrial analogues for the sulfate cycle in places such as Meridiani Planum on Mars.     

Exchange and fractionation of Mg isotopes between epsomite and saturated MgSO4 solution

The equilibrium Mg isotope fractionation factor between epsomite and aqueous MgSO₄ solution has been measured using the three isotope method in recrystallization experiments conducted at 7, 20, and 40 °C. Complete or near-complete isotopic exchange was achieved within 14 days in all experiments. The Mg isotope exchange rate between epsomite and MgSO₄ solution is dependent on the temperature, epsomite seed crystal grain size, and experimental agitation method. The Mg isotope fractionation factors (Δ²⁶Mg_eps-sol) at 7, 20, and 40 °C are 0.63 ± 0.07‰, 0.58 ± 0.16‰, and 0.56 ± 0.03‰, respectively. These values are indistinguishable within error, indicating that the Mg isotope composition of epsomite is relatively insensitive to temperature. The magnitude of the isotope fractionation factor (Δ²⁶Mg_eps-sol = ca. 0.6‰ between 7 and 40 °C) indicates that significant Mg isotope variations can be produced in evaporite sequences, and Mg isotopes may therefore, constrain the degree of closed-system behavior, paleo-humidity, and hydrological history of evaporative environments.