Experimental studies of equilibrium iron isotope fractionation in ferric aquo-chloro complexes

Here we compare new experimental studies with theoretical predictions of equilibrium iron isotopic fractionation among aqueous ferric chloride complexes (Fe(H₂O)₆³⁺, FeCl(H₂O)₅²⁺, FeCl₂(H₂O)₄⁺, FeCl₃ (H₂O)₃, and FeCl₄^-), using the Fe-Cl-H₂O system as a simple, easily-modeled example of the larger variety of iron-ligand compounds, such as chlorides, sulfides, simple organic acids, and siderophores. Isotopic fractionation (⁵⁶Fe/⁵⁴Fe) among naturally occuring iron-bearing species at Earth surface temperatures (up to ∼3‰) is usually attributed to redox effects in the environment. However, theoretical modeling of reduced isotopic partition functions among iron-bearing species in solution also predicts fractionations of similar magnitude due to non-redox changes in speciation (i.e., ligand bond strength and coordination number). In the present study, fractionations are measured in a series of low pH ([H⁺] = 5 M) solutions of ferric chloride (total Fe = 0.0749 mol/L) at chlorinities ranging from 0.5 to 5.0 mol/L. Advantage is taken of the unique solubility of FeCl₄^- in immiscible diethyl ether to create a separate spectator phase, used to monitor changing fractionation in the aqueous solution. Δ⁵⁶Fe_aq-eth = δ⁵⁶Fe (total Fe remaining in aqueous phase)−δ⁵⁶Fe (FeCl₄^- in ether phase) is determined for each solution via MC-ICPMS analysis.Both experiments and theoretical calculations of Δ⁵⁶Fe_aq-eth show a downward trend with increasing chlorinity: Δ⁵⁶Fe_aq-eth is greatest at low chlorinity, where FeCl₂(H₂O)₄⁺ is the dominant species, and smallest at high chlorinity where FeCl₃(H₂O)₃ is dominant. The experimental Δ⁵⁶Fe_aq-eth ranges from 0.8‰ at [Cl^-] = 0.5 M to 0.0‰ at [Cl^-] = 5.0 M, a decrease in aqueous-ether fractionation of 0.8‰. This is very close to the theoretically predicted decreases in Δ⁵⁶Fe_aq-eth, which range from 1.0 to 0.7‰, depending on the ab initio model.The rate of isotopic exchange and attainment of equilibrium are shown using spiked reversal experiments in conjunction with the two-phase aqueous-ether system. Equilibrium under the experimental conditions is established within 30 min.The general agreement between theoretical predictions and experimental results points to substantial equilibrium isotopic fractionation among aqueous ferric chloride complexes and a decrease in ⁵⁶Fe/⁵⁴Fe as the Cl^-/Fe³⁺ ion ratio increases. The effects on isotopic fractionation shown by the modeling of this simple iron-ligand system imply that ligands present in an aqueous environment are potentially important drivers of fractionation, are indicative of possible fractionation effects due to other speciation effects (such as iron-sulfide systems or iron bonding with organic ligands), and must be considered when interpreting iron isotope fractionation in the geological record.     

Effects of changing solution chemistry on Fe/Fe isotope fractionation in aqueous Fe-Cl solutions

The range in ⁵⁶Fe/⁵⁴Fe isotopic compositions measured in naturally occurring iron-bearing species is greater than 5‰. Both theoretical modeling and experimental studies of equilibrium isotopic fractionation among iron-bearing species have shown that significant fractionations can be caused by differences in oxidation state (i.e., redox effects in the environment) as well as by bond partner and coordination number (i.e., nonredox effects due to speciation).To test the relative effects of redox vs. nonredox attributes on total Fe equilibrium isotopic fractionation, we measured changes, both experimentally and theoretically, in the isotopic composition of an Fe²⁺-Fe³⁺-Cl-H₂O solution as the chlorinity was varied. We made use of the unique solubility of FeCl₄⁻ in immiscible diethyl ether to create a separate spectator phase against which changes in the aqueous phase could be quantified. Our experiments showed a reduction in the redox isotopic fractionation between Fe²⁺- and Fe³⁺-bearing species from 3.4‰ at [Cl⁻] = 1.5 M to 2.4‰ at [Cl⁻] = 5.0 M, due to changes in speciation in the Fe-Cl solution. This experimental design was also used to demonstrate the attainment of isotopic equilibrium between the two phases, using a ⁵⁴Fe spike.To better understand speciation effects on redox fractionation, we created four new sets of ab initio models of the ferrous chloride complexes used in the experiments. These were combined with corresponding ab initio models for the ferric chloride complexes from previous work. At 20 °C, 1000 ln β (β = ⁵⁶Fe/⁵⁴Fe reduced partition function ratio relative to a dissociated Fe atom) values range from 6.39‰ to 5.42‰ for Fe(H₂O)₆²⁺, 5.98‰ to 5.34‰ for FeCl(H₂O)₅⁺, and 5.91‰ to 4.86‰ for FeCl₂(H₂O)₄, depending on the model. The theoretical models predict ferric-ferrous fractionation about half as large (depending on model) as the experimental results.Our results show (1) oxidation state is likely to be the dominant factor controlling equilibrium Fe isotope fractionation in solution and (2) nonredox attributes (such as ligands present in the aqueous solution, speciation and relative abundances, and ionic strength of the solution) can also have significant effects. Changes in the isotopic composition of an Fe-bearing solution will influence the resultant Fe isotopic signature of any precipitates.     

Calculating the partitioning of the isotopes of Mo between oxidic and sulfidic species in aqueous solution

The fractionation of the isotopes of Mo between different geological environments has recently been determined to high accuracy using mass spectrometry (Barling et al., 2001). Fractionation is observed between Mo in seawater, where it exists primarily in the form of the Mo(VI) anion molybdate, MoO₄⁻², and in oxic sediments, where the Mo is isotopically lighter than in sea water by ∼1.8‰ (in terms of the ⁹⁷Mo, ⁹⁵Mo isotope pair). EXAFS evidence exists for a five- or six-coordinate Mo environment in the Fe,Mn oxyhydroxides of ferromanganese nodules (Kuhn et al., 2003). In sediment regimes which are anoxic and sulfidic (sometimes referred to as euxinic), where the Mo(VI) is expected to exist as a sulfide, no fractionation is observed compared to seawater. This is presumably because of the stoichiometric conversion of the Mo from MoO₄⁻² to MoS₄⁻² (Erickson and Helz, 2000) and then to other sulfides. If the conversion is stoichiometrically complete, mass balance requires the same isotopic distribution in reactant and product. This is a result of the very high equilibrium constant for this reaction. Thus, to understand isotopic fractionation processes both the equilibrium constants for the isotopic fractionation reactions and the equilbrium constants for transformation of one chemical compound to another must be considered.We here present quantum mechanical calculations of the isotopic fractionation equilibrium constants for the isotopes ⁹²Mo and ¹⁰⁰Mo between MoO₄⁻², MoO₃(OH)⁻, MoO₂(OH)₂, MoO₃, MoO₃(OH₂)₃, MoS₄⁻² and a number of other oxidic and sulfidic complexes of Mo. The fractionation equilibrium constants are calculated directly from the computed vibrational, rotational and translational contributions to the free energy in the gas-phase using quantum methods. Calculated vibrational frequencies and ratios of frequencies for different isotopomers are first obtained using a number of different quantum methods and compared with available experimental data to establish the most reliable methodology.We have also calculated free energy changes in aqueous solution for a range of reactions of MoO₄⁻² and MoO₂(OH)₂ with H₂O and H₂S. We present evidence for the instability of the monomeric octahedral species Mo(OH)₆ commonly assumed to exist in acid solution and suggest highly distorted six-coordinate MoO₃(OH₂)₃ or three-coordinate MoO₃ as better representations of the species present. We have also calculated visible-UV absorption spectra to support our interpretation of the speciation. MoO₃ is calculated to be isotopically lighter than MoO₄⁻² by 1.6‰, consistent with the experimentally observed difference between sea water and oxic sediments. We explain the isotopic lightness of oxic sediments as arising from an intermediate step in which a three coordinate MoO₃ species is formed in aqueous solution, and subsequently attaches to the surface of a Fe,Mn oxyhydroxide mineral.     

Carbon and hydrogen isotope fractionation by moderately thermophilic methanogens

A series of laboratory studies were conducted to increase understanding of stable carbon (¹³C/¹²C) and hydrogen (D/H) isotope fractionation arising from methanogenesis by moderately thermophilic acetate- and hydrogen-consuming methanogens. Studies of the aceticlastic reaction were conducted with two closely related strains of Methanosaeta thermophila. Results demonstrate a carbon isotope fractionation of only 7‰ (α = 1.007) between the methyl position of acetate and the resulting methane. Methane formed by this process is enriched in ¹³C when compared with other natural sources of methane; the magnitude of this isotope effect raises the possibility that methane produced at elevated temperature by the aceticlastic reaction could be mistaken for thermogenic methane based on carbon isotopic content. Studies of H₂/CO₂ methanogenesis were conducted with Methanothermobacter marburgensis. The fractionation of carbon isotopes between CO₂ and CH₄ was found to range from 22 to 58‰ (1.023 ≤ α ≤ 1.064). Greater fractionation was associated with low levels of molecular hydrogen and steady-state metabolism. The fractionation of hydrogen isotopes between source H₂O and CH₄ was found to range from 127 to 275‰ (1.16 ≤ α ≤ 1.43). Fractionation was dependent on growth phase with greater fractionation associated with later growth stages. The maximum observed fractionation factor was 1.43, independent of the δD-H₂ supplied to the culture. Fractionation was positively correlated with temperature and/or metabolic rate. Results demonstrate significant variability in both hydrogen and carbon isotope fractionation during methanogenesis from H₂/CO₂. The relatively small fractionation associated with deuterium during H₂/CO₂ methanogenesis provides an explanation for the relatively enriched deuterium content of biogenic natural gas originating from a variety of thermal environments. Results from these experiments are used to develop a hypothesis that differential reversibility in the enzymatic steps of the H₂/CO₂ pathway gives rise to variability in the observed carbon isotope fractionation. Results are further used to constrain the overall efficiency of electron consumption by way of the hydrogenase system in M. marburgensis, which is calculated to be less than 55%.     

O and S isotopic composition of dissolved and attached oxidation products of pyrite by Acidithiobacillus ferrooxidans: Comparison with abiotic oxidations

The acidophilic iron-oxidizing bacterium, Acidithiobacillus ferrooxidans, plays a part in the pyrite oxidation process and has been widely studied in order to determine the kinetics of the reactions and the isotopic composition of dissolved product sulphates, but the details of the oxidation processes at the surface of pyrite are still poorly known. In this study, oxygen and sulphur isotopic compositions (δ¹⁸O and δ³⁴S) were analyzed for dissolved sulphates and water from experimental aerobic acidic (pH < 2) pyrite oxidation by A. ferrooxidans. The oxidation products attached to the pyrite surfaces were studied for their morphology (SEM), their chemistry (Raman spectroscopy) and for their δ¹⁸O (ion microprobe). They were compared to abiotically (Fe³⁺, H₂O₂, O₂) oxidized pyrite surface compounds in order to constrain the oxidation pathways and to look for the existence of potential biosignatures for this system.The pyrite dissolution evolved from non-stoichiometric (during the first days) to stoichiometric (with increasing time) resulting in dissolved sulphates having distinct δ¹⁸O (e.g. +11.0‰ and −2.0‰, respectively) and δ³⁴S (+4.5‰ and +2.8‰, respectively) values. The “oxidation layer” at the surface of pyrite is complex and made of iron oxides, sulphate, polysulphide, elemental sulphur and polythionates. Bio- and Fe³⁺-oxidation favour the development of monophased micrometric bumps made of hematite or sulphate while other abiotic oxidation processes result in more variable oxidation products. The δ¹⁸O of these oxidation products at the surface of oxidized pyrites are strongly variable (from ≈−40‰ to ≈+30‰) for all experiments.Isotopic fractionation between sulphates and pyrite, Δ³⁴S_SO4-pyrite, is equal to −1.3‰ and +0.4‰ for sulphates formed by stoichiometric and non-stoichiometric processes, respectively. These two values likely reflect either a S-S or a Fe-S bond breaking process. The Δ¹⁸O_SO4-H2O and Δ¹⁸O_SO4-O2 are estimated to be ≈+16‰ and ≈−25‰, respectively. These values are higher than previously published data and may reflect biological effects. The large δ¹⁸O heterogeneity measured at the surfaces of oxidized pyrites, whatever the oxidant, may be related (i) to the existence of local surface environments isolated from the solution in which the oxidation processes are different and (ii) to the stabilization at the pyrite surface of reaction intermediates that are not in isotopic equilibrium with the solution. Though the oxygen isotopic composition of surface oxidation products cannot be taken as a direct biosignature, the combined morphological, chemical and isotopic characterization of the surfaces of oxidized pyrites may furnish clues about a biological activity on a mineral surface.     

Isotopic compositions of oxygen, iron, chromium, and nickel in cosmic spherules: Toward a better comprehension of atmospheric entry heating effects

Large, correlated, mass-dependent enrichments in the heavier isotopes of O, Cr, Fe, and Ni are observed in type-I (metal/metal oxide) cosmic spherules collected from the deep sea. Limited intraparticle variability of oxygen isotope abundances, typically <5‰ in δ¹⁸O, indicates good mixing of the melts and supports the application of the Rayleigh equation for the calculation of fractional evaporative losses during atmospheric entry. Fractional losses for oxygen evaporation from wüstite, assuming a starting isotopic composition equal to that of air (δ¹⁸O = 23.5‰; δ¹⁷O = 11.8‰), are in the range 55%-77%, and are systematically smaller than evaporative losses calculated for Fe (69%-85%), Cr (81%-95%), and especially Ni (45%-99%). However, as δ¹⁸O values increase, fractional losses for oxygen approach those of Fe, Cr, and Ni indicating a shift in the evaporating species from metallic to oxidized forms as the spherules are progressively oxidized during entry heating. The observed unequal fractional losses of O and Fe can be reconciled by allowing for a kinetic isotope mass-dependent fractionation of atmospheric oxygen during the oxidation process and/or that some metallic Fe may have undergone Rayleigh evaporation before oxidation began.In situ measurements of oxygen isotopic abundances were also performed in 14 type-S (silicate) cosmic spherules, 13 from the Antarctic ice and one from the deep sea. Additional bulk Fe and Cr isotopic abundances were determined for two type-S deep-sea spherules. The isotopic fractionation of Cr isotopes suggest appreciable evaporative loss of Cr, perhaps as a sulfide. The oxygen isotopic compositions for the type-S spherules range from δ¹⁸O = −2‰ to + 27‰. The intraspherule isotopic variations are typically small, ∼5% relative, except for the less-heated porphyritic spherules which have preserved large isotopic heterogeneities in at least one case. A plot of δ¹⁷O vs. δ¹⁸O values for these spherules defines a broad parallelogram bounded at higher values of δ¹⁷O by the terrestrial fractionation line, and at lower values of δ¹⁷O by a line parallel to it and anchored near the isotopic composition of δ¹⁸O = −2.5‰ and δ¹⁷O = −5‰. Lack of independent evidence for substantial evaporative losses suggests that much of this variation reflects the starting isotopic composition of the precursor materials, which likely resembled CO, CM, or CI chondrites. However, the enrichments in heavy isotopes indicate that some mixing with atmospheric oxygen was probably involved during atmospheric entry for some of the spherules. Isotopic fractionation due to evaporation of incoming grain is not required to explain most of the oxygen isotopic data for type-S spherules. However spherules with barred olivine textures that are thought to have experienced a more intense heating than the porphyritic ones might have undergone some distillation. Two cosmic spherules, one classified as a radial pyroxene type and the other showing a glassy texture, show unfractionated oxygen isotopic abundances. They are probably chondrule fragments that survived atmospheric entry unmelted.Possible reasons type-I spherules show larger degrees of isotopic fractionation than type-S spherules include: a) the short duration of the heating pulse associated with the high volatile content of the type-S spherule precursors compared to type-I spherules; b) higher evaporation temperatures for at least a refractory portion of the silicates compared to that of iron metal or oxide; c) lower duration of heating of type-S spherules compared to type-I spherules as a consequence of their lower densities.     

A comparison of isotope fractionation of carbon and hydrogen from paddy field rice roots and soil bacterial enrichments during CO2/H2 methanogenesis

To better understand the isotope biogeochemistry of paddy field CH₄, we investigated carbon and hydrogen isotope fractionation during CO₂ reduction by a methanogenic community enriched from California paddy field soil and rice plants. Results from analyses of terminal restriction fragment length polymorphism (T-RFLP) and sequences of the archaeal small-subunit (SSU) rRNA-encoding genes (rDNA) showed a difference in methanogenic community structure between the soil (dominated by Methanobacteriaceae) and roots (dominated by Methanospirillaceae) which was essentially the same for sampling dates 15 and 99 days after flooding (DAF). CO₂/H₂ methanogenesis by these microbial communities produced CH₄ with different isotope ratios and fractionation factors (α factors). The carbon isotope α factors in an open system with a continuous supply of 0.5% H₂ were 1.050 ± 0.002 and 1.057 ± 0.001 for soil and root enrichment cultures at 15 DAF, and 1.052 ± 0.0.002 and 1.059 ± 0.002 for soil and root enrichment cultures at 99 DAF, respectively. These α factors are similar to, but distinct from values previously obtained from cultures of mesophilic methanogens and are larger than calculated values (1.045) for paddy soil. Fractionation of hydrogen isotopes was also studied in a closed system under 80% H₂. The difference in α factors between soil and root enrichment cultures remained clear. The hydrogen isotope fractionations between culture water and the product CH₄ were −327 ± 14‰ and −319 ± 18‰ for soil enrichments, and −389 ± 17‰ and −382 ± 21‰ for root enrichments at 15 DAF and 99 DAF, respectively.     

Negative δO values in Allan Hills 84001 carbonate: Possible evidence for water precipitation on mars

The Martian meteorite ALH84001 contains ∼1% by weight of carbonate formed by secondary processes on the Martian surface or in the shallow subsurface. The major form of this carbonate is chemically and isotopically zoned rosettes which have been well documented elsewhere. This study concentrates upon carbonate regions ∼200 μm across which possess previously unobserved magnesium rich inner cores, interpreted here as rosette fragments, surrounded by a later stage cement containing rare Ca-rich carbonates (up to Ca₈₁Mg₀₇Fe₀₄Mn₀₇) intimately associated with feldspar. High spatial resolution ion probe analyses of Ca-rich carbonate surrounding rosette fragments have δ¹⁸O_V-SMOW values as low as −10‰. These values are not compatible with deposition from a global Martian atmosphere invoked to explain ALH84001 rosettes. The range of δ¹⁸O values are also incompatible with a fluid that has equilibrated with the Martian crust at high temperature or from remobilisation of carbonate of rosette isotopic composition. At Martian atmospheric temperatures, the small CO₂(gas)-CO₂(ice) fractionation makes meteoric CO₂ an unlikely source for −10‰ carbonates. In contrast, closed system Rayleigh fractionation of H₂O can generate δ¹⁸O_H2O −30‰, as observed at high latitudes on Earth. We suggest that atmospheric transport and precipitation of H₂O in a similar fashion to that on Earth provides a source of suitably ¹⁸O depleted water for generation of carbonate with δ¹⁸O_V-SMOW = −10‰.     

Kinetic nitrogen isotope fractionation associated with thermal decomposition of NH3: Experimental results and potential applications to trace the origin of N2 in natural gas and hydrothermal systems

Ammonia (NH₃) is the major intermediate phase in the pathway of nitrogen (N) transfer from the fixed N phases (e.g., in crustal material) to free N₂ (e.g., in natural gas reservoirs and volcanic gases). Yet the N isotopic behavior during these N-cycling processes remains poorly known. In an attempt to contribute to the understanding of N cycling using N isotopes, we carried out laboratory experiments to investigate the N isotopic effect associated with thermal decomposition of ammonia (2NH₃ → N₂ + 3H₂). Pure NH₃ (with initial δ¹⁵N_NH3 of ∼ −2‰, relative to air standard) was sealed into quartz tubes and thermally decomposed at 600, 700 or 800 °C from 2 hours to 500 days. With the progress of the reaction, the δ¹⁵N of the remaining NH₃ and the accumulated N₂ increased from −2 to +35‰ and from −20 to −2‰, respectively. The differences of the N-isotope fractionations at the three temperatures are not significant. Modeling using the Rayleigh distillation model yielded similar kinetic N-isotope fractionation factors (α_N2-NH3) of 0.983 ± 0.002 for 600, 700 and 800 °C. Applied to geological settings, this significant isotope discrimination (∼17‰) associated with partial decomposition of NH₃/NH₄⁺ from crustal sources (δ¹⁵N_average ∼ +6.3‰) can produce mantle-like (i.e. ∼ −5‰) or even lower δ¹⁵N values of N₂. This may explain the large variation of δ¹⁵N (−20 to +30‰) of N₂ in natural gas reservoirs. It can also possibly explain the extreme ¹⁵N-depletion of N₂ in some volcanic gases. This possibility has to be carefully considered when using N isotopes to trace geological N cycling across subduction zones by analysis of volcanic N₂.     

Non-biological fractionation of stable Ca isotopes in soils of the Atacama Desert, Chile

We measured Ca stable isotope ratios (δ^44/40Ca) in an ancient (2 My), hyperarid soil where the primary source of mobile Ca is atmospheric deposition. Most of the Ca in the upper meter of this soil (3.5 kmol m⁻²) is present as sulfates (2.5 kmol m⁻²), and to a lesser extent carbonates (0.4 kmol m⁻²). In aqueous extracts of variably hydrated calcium sulfate minerals, δ^44/40Ca_E values (vs. bulk Earth) increase with depth (1.4 m) from a minimum of −1.91‰ to a maximum of +0.59‰. The trend in carbonate-δ^44/40Ca in the top six horizons resembles that of sulfate-δ^44/40Ca, but with values 0.1-0.6‰ higher. The range of observed Ca isotope values in this soil is about half that of δ^44/40Ca values observed on Earth. Linear correlation among δ^44/40Ca, δ³⁴S and δ¹⁸O values indicates either (a) a simultaneous change in atmospheric input values for all three elements over time, or (b) isotopic fractionation of all three elements during downward transport. We present evidence that the latter is the primary cause of the isotopic variation that we observe. Sulfate-δ³⁴S values are positively correlated with sulfate-δ¹⁸O values (R² = 0.78) and negatively correlated with sulfate δ^44/40Ca_E values (R² = 0.70). If constant fractionation and conservation of mass with downward transport are assumed, these relationships indicate a δ^44/40Ca fractionation factor of −0.4‰ in CaSO₄. The overall depth trend in Ca isotopes is reproduced by a model of isotopic fractionation during downward Ca transport that considers small and infrequent but regularly recurring rainfall events. Near surface low Ca isotope values are reproduced by a Rayleigh model derived from measured Ca concentrations and the Ca fractionation factor predicted by the relationship with S isotopes. This indicates that the primary mechanism of stable isotope fractionation in CaSO₄ is incremental and effectively irreversible removal of an isotopically enriched dissolved phase by downward transport during small rainfall events.     

Carbon isotope effects in the open-system Fischer-Tropsch synthesis

Carbon isotopic composition was measured for products of the Fischer-Tropsch synthesis: catalytic reaction between CO and H₂ to produce CO, CO₂, light hydrocarbons C1-C4 and “oil” fraction. Hydrogen isotopes were also measured in the oil fraction and the produced water. Experimental runs were conducted in the flow-through reactor at 260-310 °C and 30 bar using the synthesis gas composed of 5N₂ + 3H₂ + 2CO, on Fe-catalyst mixed with ZSM-5 synthetic zeolite. In the two of seven runs a Fe + Co-catalyst was used that gives a lower yield of unsaturated hydrocarbons in reaction products. The isotopic effects depended on the conversion of the carbon monoxide. Under steady-state conditions (CO conversion more than 90%) a strong kinetic fractionation was observed between CO and CO₂ (∼−10‰) and CO and hydrocarbons (∼+38‰). At low conversion a clear “inverse” isotopic trend of the depletion in ¹³C of longer hydrocarbon chains was observed. On average, Δ₁₂ = δ¹³C(CH₄) − δ¹³C(C₂H₆) correlates well with the CO conversion: the C₂H₆ is ∼6‰ isotopically lighter than CH₄ at low conversion and ∼2‰ heavier at steady-state regime. Under steady-state conditions there almost no difference was observed in the isotopic composition of methane and ethane and higher hydrocarbons. The chemical composition of light hydrocarbons in the products of flow-through, dynamic FTS is different from that found in the static FTS-type experiments with Fe-catalyst, but isotopic effects are similar. Our results suggest that the isotopic distribution of carbon found in so-called “abiogenic” hydrocarbons from some natural gases (δ¹³C1 > δ¹³C2 > δ¹³C3  >?) is somewhat similar to that at low conversion of CO, but do not resemble the distribution characteristic for the high conversion products, at least, on Fe-catalyst. Other processes (a simple mixing of two or more endmembers) or other P-T conditions of the carbon reduction could be responsible for the “inverse” isotopic trend found in meteorites and some natural gases.     

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

Application of the Fe isotope system to studies of natural rocks and fluids requires precise knowledge of equilibrium Fe isotope fractionation factors among various aqueous Fe species and minerals. These are difficult to obtain at the low temperatures at which Fe isotope fractionation is expected to be largest and requires careful distinction between kinetic and equilibrium isotope effects. A detailed investigation of Fe isotope fractionation between [Fe^III(H₂O)₆]³⁺ and hematite at 98°C allows the equilibrium ⁵⁶Fe/⁵⁴Fe fractionation to be inferred, which we estimate at 10³lnα_{Fe(III)-hematite} = −0.10 ± 0.20‰. We also infer that the slope of Fe(III)-hematite fractionation is modest relative to 10⁶/T², which would imply that this fractionation remains close to zero at lower temperatures. These results indicate that Fe isotope compositions of hematite may closely approximate those of the fluids from which they precipitated if equilibrium isotopic fractionation is assumed, allowing inference of δ⁵⁶Fe values of ancient fluids from the rock record. The equilibrium Fe(III)-hematite fractionation factor determined in this study is significantly smaller than that obtained from the reduced partition function ratios calculated for [Fe^III(H₂O)₆]³⁺ and hematite based on vibrational frequencies and Mössbauer shifts by [Polyakov 1997] and [Polyakov and Mineev 2000], and Schauble et al. (2001), highlighting the importance of experimental calibration of Fe isotope fractionation factors. In contrast to the long-term (up to 203 d) experiments, short-term experiments indicate that kinetic isotope effects dominate during rapid precipitation of ferric oxides. Precipitation of hematite over ∼12 h produces a kinetic isotope fractionation where 10³lnα_{Fe(III)-hematite} = +1.32 ± 0.12‰. Precipitation under nonequilibrium conditions, however, can be recognized through stepwise dissolution in concentrated acids. As expected, our results demonstrate that dissolution by itself does not measurably fractionate Fe isotopes.     

Zn isotopic fractionation caused by sorption on goethite and 2-Lines ferrihydrite

Zn isotopic fractionation caused by sorption on 2-Lines ferrihydrite (Fh2L) and goethite was investigated to assess the role of reactions at the Fe-oxyhydroxide/water interface in changes of the isotopic distribution of Zn. Since sorption reactions are ubiquitous in Earth’s surface environments, it is important to evaluate their influence on the isotopic distribution of Zn before it can be used to track and quantify contributions of various sources and/or biogeochemical processes involving this element. Our results show that Zn isotopes are fractionated upon sorption on Fe-oxyhydroxides with an enrichment of the heavy isotopes present on the solid’s surface. This fractionation appears to proceed through an equilibrium mechanism and yields different (Δ^66/64Zn)_{sorbed-aqueous} values for Zn sorption on goethite [(Δ^66/64Zn)_{sorbed-aqueous} around +0.29‰] and Fh2L [(Δ^66/64Zn)_{sorbed-aqueous} around +0.53‰]. These different magnitudes of Zn fractionation are related to structural differences between Zn complexes existing on the surface of goethite (octahedrally coordinated Zn by oxygen atoms) and Fh2L (tetrahedrally coordinated Zn by oxygen atoms), as evidenced by Extended X-ray Absorption Fine Structure (EXAFS) spectroscopy and CD-MUSIC modeling. These results show the importance of accounting for reactions at the Fe-oxyhydroxide/water interface when dealing with the isotopic distribution of Zn at the Earth’s surface. Considering the large range of other possible sorbents (Mn or Al oxides, phyllosilicates, carbonates, biologic surfaces, etc.) and the importance of reactions at sorbent/water interfaces for other non-traditional stable isotopes (i.e. Cr, Fe, Ni and Cu) that are increasingly used in environmental studies, these results emphasize the need for further experimental studies that are needed to quantify the isotopic fractionation of these elements possibly accompanying their sorption.     

Model for kinetic effects on calcium isotope fractionation (δCa) in inorganic aragonite and cultured planktonic foraminifera

The calcium isotope ratios (δ⁴⁴Ca = [(⁴⁴Ca/⁴⁰Ca)_sample/(⁴⁴Ca/⁴⁰Ca)_standard −1] · 1000) of Orbulina universa and of inorganically precipitated aragonite are positively correlated to temperature. The slopes of 0.019 and 0.015‰ °C⁻¹, respectively, are a factor of 13 and 16 times smaller than the previously determined fractionation from a second foraminifera, Globigerinoides sacculifer, having a slope of about 0.24‰ °C⁻¹. The observation that δ⁴⁴Ca is positively correlated to temperature is opposite in sign to the oxygen isotopic fractionation (δ¹⁸O) in calcium carbonate (CaCO₃). These observations are explained by a model which considers that Ca²⁺-ions forming ionic bonds are affected by kinetic fractionation only, whereas covalently bound atoms like oxygen are affected by kinetic and equilibrium fractionation. From thermodynamic consideration of kinetic isotope fractionation, it can be shown that the slope of the enrichment factor α(T) is mass-dependent. However, for O. universa and the inorganic precipitates, the calculated mass of about 520 ± 60 and 640 ± 70 amu (atomic mass units) is not compatible with the expected ion mass for ⁴⁰Ca and ⁴⁴Ca. To reconcile this discrepancy, we propose that Ca diffusion and δ⁴⁴Ca isotope fractionation at liquid/solid transitions involves Ca²⁺-aquocomplexes (Ca[H₂O]_n²⁺ · mH₂O) rather than pure Ca²⁺-ion diffusion. From our measurements we calculate that such a hypothesized Ca²⁺-aquocomplex correlates to a hydration number of up to 25 water molecules (490 amu). For O. universa we propose that their biologically mediated Ca isotope fractionation resembles fractionation during inorganic precipitation of CaCO₃ in seawater. To explain the different Ca isotope fractionation in O. universa and in G. sacculifer, we suggest that the latter species actively dehydrates the Ca²⁺-aquocomplex before calcification takes place. The very different temperature response of Ca isotopes in the two species suggests that the use of δ⁴⁴Ca as a temperature proxy will require careful study of species effects.     

Calibration of the calcite-water oxygen-isotope geothermometer at Devils Hole, Nevada, a natural laboratory

The δ¹⁸O of ground water (−13.54 ± 0.05 ‰) and inorganically precipitated Holocene vein calcite (+14.56 ± 0.03 ‰) from Devils Hole cave #2 in southcentral Nevada yield an oxygen isotopic fractionation factor between calcite and water at 33.7 °C of 1.02849 ± 0.00013 (1000 ln α_{calcite-water} = 28.09 ± 0.13). Using the commonly accepted value of ∂(α_{calcite-water})/∂T of −0.00020 K⁻¹, this corresponds to a 1000 ln α_{calcite-water} value at 25 °C of 29.80, which differs substantially from the current accepted value of 28.3. Use of previously published oxygen isotopic fractionation factors would yield a calcite precipitation temperature in Devils Hole that is 8 °C lower than the measured ground water temperature. Alternatively, previously published fractionation factors would yield a δ¹⁸O of water, from which the calcite precipitated, that is too negative by 1.5 ‰ using a temperature of 33.7 °C. Several lines of evidence indicate that the geochemical environment of Devils Hole has been remarkably constant for at least 10 ka. Accordingly, a re-evaluation of calcite-water oxygen isotopic fractionation factor may be in order.Assuming the Devils Hole oxygen isotopic value of α_{calcite-water} represents thermodynamic equilibrium, many marine carbonates are precipitated with a δ¹⁸O value that is too low, apparently due to a kinetic isotopic fractionation that preferentially enriches ¹⁶O in the solid carbonate over ¹⁸O, feigning oxygen isotopic equilibrium.     

Stable isotope geochemistry of ground and surface waters associated with undisturbed massive sulfide deposits; constraints on origin of waters and water-rock reactions

Stable isotopes (H, O, C) were determined for ground and surface waters collected from two relatively undisturbed massive sulfide deposits (Halfmile Lake and Restigouche) in the Bathurst Mining Camp (BMC), New Brunswick, Canada. Additional waters from active and inactive mines in the BMC were also collected. Oxygen and hydrogen isotopes of surface and shallow groundwaters from both the Halfmile Lake and Restigouche deposits are remarkably uniform (− 13 to − 14‰ and − 85 to − 95‰ for δ¹⁸O_VSMOW and δ²H_VSMOW, respectively). These values are lighter than predicted for northern New Brunswick and, combined with elevated deuterium excess values, suggest that recharge waters are dominated by winter precipitation, recharged during spring melting. Deeper groundwaters from the Restigouche deposit, and from active and inactive mines have heavier δ¹⁸O_VSMOW ratios (up to − 10.8‰) than shallow groundwaters suggesting recharge under warmer climate or mixing with Shield-type brines. Some of the co-variation in Cl concentrations and δ¹⁸O_VSMOW ratios can be explained by mixing between saline and shallow recharge water end-members. Carbon isotopic compositions of dissolved inorganic carbon (DIC) are variable, ranging from − 15 to − 5‰ δ¹³C_VPDB for most ground and surface waters. Much of the variation in the carbon isotopes is consistent with closed system groundwater evolution involving soil zone CO₂ and fracture zone carbonate minerals (calcite, dolomite and siderite; average = − 6.5‰ δ¹³C_VPDB). The DIC of saline Restigouche deposit groundwater is isotopically heavy (∼+ 12‰ δ¹³C_VPDB), indicating carbon isotopic fractionation from methanogenesis via CO₂ reduction, consistent with the lack of dissolved sulfate in these waters and the observation of CH₄-degassing during sampling.     

Interpreting the Ca isotope record of marine biogenic carbonates

An 18 million year record of the Ca isotopic composition (δ^44/42Ca) of planktonic foraminiferans from ODP site 925, in the Atlantic, on the Ceara Rise, provides the opportunity for critical analysis of Ca isotope-based reconstructions of the Ca cycle. δ^44/42Ca in this record averages +0.37 ± 0.05 (1σ SD) and ranges from +0.21‰ to +0.52‰. The record is a good match to previously published Neogene Ca isotope records based on foraminiferans, but is not similar to the record based on bulk carbonates, which has values that are as much as 0.25‰ lower. Bulk carbonate and planktonic foraminiferans from core tops differ slightly in their δ^44/42Ca (i.e., by 0.06 ± 0.06‰ (n = 5)), while the difference between bulk carbonate and foraminiferan values further back in time is markedly larger, leaving open the question of the cause of the difference. Modeling the global Ca cycle from downcore variations in δ^44/42Ca by assuming fixed values for the isotopic composition of weathering inputs (δ^44/42Ca_w) and for isotope fractionation associated with the production of carbonate sediments (Δ_sed) results in unrealistically large variations in the total mass of Ca²⁺ in the oceans over the Neogene. Alternatively, variations of ±0.05‰ in the Ca isotope composition of weathering inputs or in the extent of fractionation of Ca isotopes during calcareous sediment formation could entirely account for variations in the Ca isotopic composition of marine carbonates. Ca isotope fractionation during continental weathering, such as has been recently observed, could easily result in variations in δ^44/42Ca_w of a few tenths of permil. Likewise a difference in the fractionation factors associated with aragonite versus calcite formation could drive shifts in Δ_sed of tenths of permil with shifts in the relative output of calcite and aragonite from the ocean. Until better constraints on variations in δ^44/42Ca_w and Δ_sed have been established, modeling the Ca²⁺ content of seawater from Ca isotope curves should be approached cautiously.     

Impact of soil weathering degree on silicon isotopic fractionation during adsorption onto iron oxides in basaltic ash soils, Cameroon

The sequestration of silicon in soil clay-sized iron oxides may affect the terrestrial cycle of Si. Iron oxides indeed specifically adsorb aqueous monosilicic acid (H₄SiO₄⁰), thereby influencing Si concentration in soil solution. Here we study the impact of H₄SiO₄⁰ adsorption on the fractionation of Si isotopes in basaltic ash soils differing in weathering degree (from two weathering sequences, Cameroon), hence in clay and Fe-oxide contents, and evaluate the potential isotopic impact on dissolved Si in surrounding Cameroon rivers. Adsorption was measured in batch experiment series designed as function of time (0-72 h) and initial concentration (ic) of Si in solution (0.61-1.18 mM) at 20 °C, constant pH (5.5) and ionic strength (1 mM). After various soil-solution contact times, the δ³⁰Si vs. NBS28 compositions were determined in selected solutions by MC-ICP-MS (Nu Plasma) in medium resolution, operating in dry plasma with Mg doping with an average precision of ±0.15‰ (±2σ_SEM). The quantitative adsorption of H₄SiO₄⁰ by soil Fe-oxides left a solution depleted in light Si isotopes, which confirms previous study on synthetic Fe-oxides. Measured against its initial composition (δ³⁰Si = +0.02 ± 0.07‰ (±2σ_SD)), the solutions were systematically enriched in ³⁰Si reaching maximum δ³⁰Si values ranging between +0.16‰ and +0.95‰ after 72 h contact time. The enrichment of the solution in heavy isotopes increased with increasing values of three parameters: soil weathering degree, iron oxide content, and proportion of short-range ordered Fe-oxide. The Si-isotopic signature of the solution was partly influenced by Si release, possibly through mineral dissolution and Si desorption from oxide surfaces, depending on soil type, highlighting the complex pattern of natural soils. Surrounding Cameroon rivers displayed a mean Si-isotopic signature of +1.19‰. Our data imply that in natural environments, H₄SiO₄⁰ adsorption by soil clay-sized Fe-oxides at least partly impacts the Si-isotopic signature of the soil solution exported to water streams.     

Sources and biological fractionation of Silicon isotopes in the Eastern Equatorial Pacific

Silicon isotopes in dissolved silicic acid were measured in the upper four kilometers between 4°N and 3°S latitude at 110°W longitude in the eastern Equatorial Pacific. Silicon isotopes became progressively heavier with silicic acid depletion of surface water as expected from biological fractionation. The value of ε estimated by applying a steady-state isotope fractionation model to data from all stations between 4°N and 3°S was −0.77 ± 0.12‰ (std. err.). When the analysis was restricted to those stations whose temperature and salinity profiles indicated that they were directly influenced by upwelling of the Equatorial Undercurrent (EUC), the resulting value of ε was −1.08 ± 0.27‰ (std. err.) similar to the value established in culture studies (−1.1‰). When the non steady state Rayleigh model was applied to the same restricted data set the resulting value of ε was significantly more positive, −0.61 ± 0.16‰ (std. err.). To the extent that the equatorial system approximates a steady state these results support a value of −1.1‰ for the fractionation factor for isotopes of Si in the sea. Without the assumption of steady state the value of ε can only be constrained to be between −0.6 and −1.1‰. Silicic acid in Equatorial Pacific Deep Water below 2000 m had a near constant δ³⁰Si of +1.32 ± 0.05‰. That value is significantly more positive than obtained for North Pacific Deep Water at similar depths at stations to the northwest of our study area (0.9-1.0‰) and it is slightly less positive than new measures of the δ³⁰Si of silicic acid from the silicic acid plume centered over the Cascadia basin in the Northeast Pacific (Si(OH)₄ > 180  μM, δ³⁰Si = +1.46 ± 0.12‰ (SD, n = 4). We show that the data from the equator and Cascadia basin fit a general trend of increasing δ³⁰Si(OH)₄ with increasing silicic acid concentration in the deep sea, but that the isotope values from the Northeast Pacific are anomalously light. The observed level of variation in the silicon isotope composition of deep waters from this single ocean basin is considerably larger than that predicted by current models based on fractionation during opal formation with no isotope effect during dissolution. Confirmation of such high variability in deep water δ³⁰Si(OH)₄ within individual ocean basins will require reassessment of the mechanisms controlling the distribution of isotopes of silicon in the sea.     

Rare sulfur and triple oxygen isotope geochemistry of volcanogenic sulfate aerosols

We present analyses of stable isotopic ratios ¹⁷O/¹⁶O, ¹⁸O/¹⁶O, ³⁴S/³²S, and ³³S/³²S, ³⁶S/³²S in sulfate leached from volcanic ash of a series of well known, large and small volcanic eruptions. We consider eruptions of Mt. St. Helens (Washington, 1980, ∼1 km³), Mt. Spurr (Alaska, 1953, <1 km³), Gjalp (Iceland, 1996, 1998, <1 km³), Pinatubo (Phillipines, 1991, 10 km³), Bishop tuff (Long Valley, California, 0.76 Ma, 750 km³), Lower Bandelier tuff (Toledo Caldera, New Mexico, 1.61 Ma, 600 km³), and Lava Creek and Huckleberry Ridge tuffs (Yellowstone, Wyoming, 0.64 Ma, 1000 km³ and 2.04 Ma 2500 km³, respectively). This list covers much of the diversity of sizes and the character of silicic volcanic eruptions. Particular emphasis is paid to the Lava Creek tuff for which we present wide geographic sample coverage.This global dataset spans a significant range in δ³⁴S, δ¹⁸O, and Δ¹⁷O of sulfate (29‰, 30‰, and 3.3‰, respectively) with oxygen isotopes recording mass-independent (Δ¹⁷O > 0.2‰) and sulfur isotopes exhibiting mass-dependent behavior. Products of large eruptions account for most of‘ these isotopic ranges. Sulfate with Δ¹⁷O > 0.2‰ is present as 1-10 μm gypsum crystals on distal ash particles and records the isotopic signature of stratospheric photochemical reactions. Sediments that embed ash layers do not contain sulfate or contain little sulfate with Δ¹⁷O near 0‰, suggesting that the observed sulfate in ash is of volcanic origin.Mass-dependent fractionation of sulfur isotopic ratios suggests that sulfate-forming reactions did not involve photolysis of SO₂, like that inferred for pre-2.3 Ga sulfates from Archean sediments or Antarctic ice-core sulfate associated with few dated eruptions. Even though the sulfate sulfur isotopic compositions reflect mass-dependent processes, the products of caldera-forming eruptions display a large δ³⁴S range and exhibit fractionation relationships that do not follow the expected equilibrium slopes of 0.515 and 1.90 for ³³S/³²S vs. ³⁴S/³²S and ³⁶S/³²S vs. ³⁴S/³²S, respectively. The data presented here are consistent with modification of a chemical mass-dependent fractionation of sulfur isotopes in the volcanic plume by either a kinetic gas phase reaction of volcanic SO₂ with OH and/or a Rayleigh processes involving a residual Rayleigh reactant—volcanic SO₂ gas, rather than a Rayleigh product. These results may also imply at least two removal pathways for SO₂ in volcanic plumes.Above-zero Δ¹⁷O values and their positive correlation with δ¹⁸O in sulfate can be explained by oxidation by high-δ¹⁸O and high-Δ¹⁷O compounds such as ozone and radicals such as OH that result from ozone break down. Large caldera-forming eruptions have the highest Δ¹⁷O values, and the largest range of δ¹⁸O, which can be explained by stratospheric reaction with ozone-derived OH radicals. These results suggest that massive eruptions are capable of causing a temporary depletion of the ozone layer. Such depletion may be many times that of the measured 3-8% depletion following 1991 Pinatubo eruption, if the amount of sulfur dioxide released scales with the amount of ozone depletion.