Experimental constraints on Fe isotope fractionation during magnetite and Fe carbonate formation coupled to dissimilatory hydrous ferric oxide reduction

Iron isotope fractionation between aqueous Fe(II) and biogenic magnetite and Fe carbonates produced during reduction of hydrous ferric oxide (HFO) by Shewanella putrefaciens, Shewanella algae, and Geobacter sulfurreducens in laboratory experiments is a function of Fe(III) reduction rates and pathways by which biogenic minerals are formed. High Fe(III) reduction rates produced ⁵⁶Fe/⁵⁴Fe ratios for Fe(II)_aq that are 2-3‰ lower than the HFO substrate, reflecting a kinetic isotope fractionation that was associated with rapid sorption of Fe(II) to HFO. In long-term experiments at low Fe(III) reduction rates, the Fe(II)_aq-magnetite fractionation is −1.3‰, and this is interpreted to be the equilibrium fractionation factor at 22°C in the biologic reduction systems studied here. In experiments where Fe carbonate was the major ferrous product of HFO reduction, the estimated equilibrium Fe(II)_aq-Fe carbonate fractionations were ca. 0.0‰ for siderite (FeCO₃) and ca. +0.9‰ for Ca-substituted siderite (Ca_0.15Fe_0.85CO₃) at 22°C. Formation of precursor phases such as amorphous nonmagnetic, noncarbonate Fe(II) solids are important in the pathways to formation of biogenic magnetite or siderite, particularly at high Fe(III) reduction rates, and these solids may have ⁵⁶Fe/⁵⁴Fe ratios that are up to 1‰ lower than Fe(II)_aq. Under low Fe(III) reduction rates, where equilibrium is likely to be attained, it appears that both sorbed Fe(II) and amorphous Fe(II)(s) components have isotopic compositions that are similar to those of Fe(II)_aq.The relative order of δ⁵⁶Fe values for these biogenic minerals and aqueous Fe(II) is: magnetite > siderite ≈ Fe(II)_aq > Ca-bearing Fe carbonate, and this is similar to that observed for minerals from natural samples such as Banded Iron Formations (BIFs). Where magnetite from BIFs has δ⁵⁶Fe >0‰, the calculated δ⁵⁶Fe value for aqueous Fe(II) suggests a source from midocean ridge (MOR) hydrothermal fluids. In contrast, magnetite from BIFs that has δ⁵⁶Fe ≤0‰ apparently requires formation from aqueous Fe(II) that had very low δ⁵⁶Fe values. Based on this experimental study, formation of low-δ⁵⁶Fe Fe(II)_aq in nonsulfidic systems seems most likely to have been produced by dissimilatory reduction of ferric oxides by Fe(III)-reducing bacteria.     

Oxygen and iron isotope constraints on near-surface fractionation effects and the composition of lunar mare basalt source regions

Oxygen and iron isotope analyses of low-Ti and high-Ti mare basalts are presented to constrain their petrogenesis and to assess stable isotope variations within lunar mantle sources. An internally-consistent dataset of oxygen isotope compositions of mare basalts encompasses five types of low-Ti basalts from the Apollo 12 and 15 missions and eight types of high-Ti basalts from the Apollo 11 and 17 missions. High-precision whole-rock δ¹⁸O values (referenced to VSMOW) of low-Ti and high-Ti basalts correlate with major-element compositions (Mg#, TiO₂, Al₂O₃). The observed oxygen isotope variations within low-Ti and high-Ti basalts are consistent with crystal fractionation and match the results of mass-balance models assuming equilibrium crystallization. Whole-rock δ⁵⁶Fe values (referenced to IRMM-014) of high-Ti and low-Ti basalts range from 0.134‰ to 0.217‰ and 0.038‰ to 0.104‰, respectively. Iron isotope compositions of both low-Ti and high-Ti basalts do not correlate with indices of crystal fractionation, possibly owing to small mineral-melt iron fractionation factors anticipated under lunar reducing conditions.The δ¹⁸O and δ⁵⁶Fe values of low-Ti and the least differentiated high-Ti mare basalts are negatively correlated, which reflects their different mantle source characteristics (e.g., the presence or absence of ilmenite). The average δ⁵⁶Fe values of low-Ti basalts (0.073 ± 0.018‰, n = 8) and high-Ti basalts (0.191 ± 0.020‰, n = 7) may directly record that of their parent mantle sources. Oxygen isotope compositions of mantle sources of low-Ti and high-Ti basalts are calculated using existing models of lunar magma ocean crystallization and mixing, the estimated equilibrium mantle olivine δ¹⁸O value, and equilibrium oxygen-fractionation between olivine and other mineral phases. The differences between the calculated whole-rock δ¹⁸O values for source regions, 5.57‰ for low-Ti and 5.30‰ for high-Ti mare basalt mantle source regions, are solely a function of the assumed source mineralogy. The oxygen and iron isotope compositions of lunar upper mantle can be approximated using these mantle source values. The δ¹⁸O and δ⁵⁶Fe values of the lunar upper mantle are estimated to be 5.5 ± 0.2‰ (2σ) and 0.085 ± 0.040‰ (2σ), respectively. The oxygen isotope composition of lunar upper mantle is identical to the current estimate of Earth’s upper mantle (5.5 ± 0.2‰), and the iron isotope composition of the lunar upper mantle overlaps within uncertainty of estimates for the terrestrial upper mantle (0.044 ± 0.030‰).     

Iron isotope fractionation during microbially stimulated Fe(II) oxidation and Fe(III) precipitation

Interpretation of the origins of iron-bearing minerals preserved in modern and ancient rocks based on measured iron isotope ratios depends on our ability to distinguish between biological and non-biological iron isotope fractionation processes. In this study, we compared ⁵⁶Fe/⁵⁴Fe ratios of coexisting aqueous iron (Fe(II)_aq, Fe(III)_aq) and iron oxyhydroxide precipitates (Fe(III)_ppt) resulting from the oxidation of ferrous iron under experimental conditions at low pH (<3). Experiments were carried out using both pure cultures of Acidothiobacillus ferrooxidans and sterile controls to assess possible biological overprinting of non-biological fractionation, and both SO₄²⁻ and Cl⁻ salts as Fe(II) sources to determine possible ionic/speciation effects that may be associated with oxidation/precipitation reactions. In addition, a series of ferric iron precipitation experiments were performed at pH ranging from 1.9 to 3.5 to determine if different precipitation rates cause differences in the isotopic composition of the iron oxyhydroxides. During microbially stimulated Fe(II) oxidation in both the sulfate and chloride systems, ⁵⁶Fe/⁵⁴Fe ratios of residual Fe(II)_aq sampled in a time series evolved along an apparent Rayleigh trend characterized by a fractionation factor α_{Fe(III)aq-Fe(II)aq} ∼ 1.0022. This fractionation factor was significantly less than that measured in our sterile control experiments (∼1.0034) and that predicted for isotopic equilibrium between Fe(II)_aq and Fe(III)_aq (∼1.0029), and thus might be interpreted to reflect a biological isotope effect. However, in our biological experiments the measured difference in ⁵⁶Fe/⁵⁴Fe ratios between Fe(III)_aq, isolated as a solid by the addition of NaOH to the final solution at each time point under N₂-atmosphere, and Fe(II)_aq was in most cases and on average close to 2.9‰ (α_{Fe(III)aq-Fe(II)aq} ∼ 1.0029), consistent with isotopic equilibrium between Fe(II)_aq and Fe(III)_aq. The ferric iron precipitation experiments revealed that ⁵⁶Fe/⁵⁴Fe ratios of Fe(III)_aq were generally equal to or greater than those of Fe(III)_ppt, and isotopic fractionation between these phases decreased with increasing precipitation rate and decreasing grain size. Considered together, the data confirm that the iron isotope variations observed in our microbial experiments are primarily controlled by non-biological equilibrium and kinetic factors, a result that aids our ability to interpret present-day iron cycling processes but further complicates our ability to use iron isotopes alone to identify biological processing in the rock record.     

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.     

Iron isotope fractionation between pyrite (FeS2), hematite (Fe2O3) and siderite (FeCO3): A first-principles density functional theory study

In addition to equilibrium isotopic fractionation factors experimentally derived, theoretical predictions are needed for interpreting isotopic compositions measured on natural samples because they allow exploring more easily a broader range of temperature and composition. For iron isotopes, only aqueous species were studied by first-principles methods and the combination of these data with those obtained by different methods for minerals leads to discrepancies between theoretical and experimental isotopic fractionation factors. In this paper, equilibrium iron isotope fractionation factors for the common minerals pyrite, hematite, and siderite were determined as a function of temperature, using first-principles methods based on the density functional theory (DFT). In these minerals belonging to the sulfide, oxide and carbonate class, iron is present under two different oxidation states and is involved in contrasted types of interatomic bonds. Equilibrium fractionation factors calculated between hematite and siderite compare well with the one estimated from experimental data (ln α⁵⁷Fe/⁵⁴Fe = 4.59 ± 0.30‰ and 5.46 ± 0.63‰ at 20 °C for theoretical and experimental data, respectively) while those for Fe(III)_aq-hematite and Fe(II)_aq-siderite are significantly higher that experimental values. This suggests that the absolute values of the reduced partition functions (β-factors) of aqueous species are not accurate enough to be combined with those calculated for minerals. When compared to previous predictions derived from Mössbauer or INRXS data [Polyakov V. B., Clayton R. N., Horita J. and Mineev S. D. (2007) Equilibrium iron isotope fractionation factors of minerals: reevaluation from the data of nuclear inelastic resonant X-ray scattering and Mössbauer spectroscopy. Geochim. Cosmochim. Acta71, 3833-3846], our iron β-factors are in good agreement for siderite and hematite while a discrepancy is observed for pyrite. However, the detailed investigation of the structural, electronic and vibrational properties of pyrite as well as the study of sulfur isotope fractionation between pyrite and two other sulfides (sphalerite and galena) indicate that DFT-derived β-factors of pyrite are as accurate as for hematite and siderite. We thus suggest that experimental vibrational density of states of pyrite should be re-examined.     

Kinetic fractionation of Fe isotopes during transport through a porous quartz-sand column

Sorption and desorption processes are an important part of biological and geochemical metallic isotope cycles. Here, we address the dynamic aspects of metallic isotopic fractionation in a theoretical and experimental study of Fe sorption and desorption during the transport of aqueous Fe(III) through a quartz-sand matrix. Transport equations describing the behavior of sorbing isotopic species in a water saturated homogeneous porous medium are presented; isotopic fractionation of the system (Δ_{sorbedmetal-soln}) being defined in terms of two parameters: (i) an equilibrium fractionation factor, α_e; and (ii) a kinetic sorption factor, α₁. These equations are applied in a numerical model that simulates the sorption-desorption of Fe isotopes during injection of a Fe(III) solution pulse into a quartz matrix at pH 0-2 and explores the effects of the kinetic and equilibrium parameters on the Fe-isotope evolution of porewater. The kinetic transport theory is applied to a series of experiments in which pulses of Na and Fe(III) chloride solutions were injected into a porous sand grain column. Fractionation factors of α_e = 1.0003 ± 0.0001 and α₁ = 0.9997 ± 0.0004 yielded the best fit between the transport model and the Fe concentration and δ⁵⁶Fe data. The equilibrium fractionation (Δ⁵⁶Fe_{sorbedFe-soln}) of 0.3‰ is comparable with values deduced for adsorption of metallic cations on iron and manganese oxide surfaces and suggests that sandstone aquifers will fractionate metallic isotopes during sorption-desorption reactions. The ability of the equilibrium fractionation factor to describe a natural system, however, depends on the proximity to equilibrium, which is determined by the relative time scales of mass transfer and chemical reaction; low fluid transport rates should produce a system that is less dependent on kinetic effects. The results of this study are applicable to Fe-isotope fractionation in clastic sediments formed in highly acidic conditions; such conditions may have existed on Mars where acidic oxidizing ground and surface waters may have been responsible for clastic sedimentation and metallic element transport.     

Evidence for equilibrium iron isotope fractionation by nitrate-reducing iron(II)-oxidizing bacteria

Iron isotope fractionations produced during chemical and biological Fe(II) oxidation are sensitive to the proportions and nature of dissolved and solid-phase Fe species present, as well as the extent of isotopic exchange between precipitates and aqueous Fe. Iron isotopes therefore potentially constrain the mechanisms and pathways of Fe redox transformations in modern and ancient environments. In the present study, we followed in batch experiments Fe isotope fractionations between Fe(II)_aq and Fe(III) oxide/hydroxide precipitates produced by the Fe(III) mineral encrusting, nitrate-reducing, Fe(II)-oxidizing Acidovorax sp. strain BoFeN1. Isotopic fractionation in ⁵⁶Fe/⁵⁴Fe approached that expected for equilibrium conditions, assuming an equilibrium Δ⁵⁶Fe_{Fe(OH)3-Fe(II)aq} fractionation factor of +3.0‰. Previous studies have shown that Fe(II) oxidation by this Acidovorax strain occurs in the periplasm, and we propose that Fe isotope equilibrium is maintained through redox cycling via coupled electron and atom exchange between Fe(II)_aq and Fe(III) precipitates in the contained environment of the periplasm. In addition to the apparent equilibrium isotopic fractionation, these experiments also record the kinetic effects of initial rapid oxidation, and possible phase transformations of the Fe(III) precipitates. Attainment of Fe isotope equilibrium between Fe(III) oxide/hydroxide precipitates and Fe(II)_aq by neutrophilic, Fe(II)-oxidizing bacteria or through abiologic Fe(II)_aq oxidation is generally not expected or observed, because the poor solubility of their metabolic product, i.e. Fe(III), usually leads to rapid precipitation of Fe(III) minerals, and hence expression of a kinetic fractionation upon precipitation; in the absence of redox cycling between Fe(II)_aq and precipitate, kinetic isotope fractionations are likely to be retained. These results highlight the distinct Fe isotope fractionations that are produced by different pathways of biological and abiological Fe(II) oxidation.     

Iron isotope fractionation by Fe(II)-oxidizing photoautotrophic bacteria

Photoautotrophic bacteria that oxidize ferrous iron (Fe[II]) under anaerobic conditions are thought to be ancient in origin, and the ferric (hydr)oxide mineral products of their metabolism are likely to be preserved in ancient rocks. Here, two enrichment cultures of Fe(II)-oxidizing photoautotrophs and a culture of the genus Thiodictyon were studied with respect to their ability to fractionate Fe isotopes. Fe isotope fractionations produced by both the enrichment cultures and the Thiodictyon culture were relatively constant at early stages of the reaction progress, where the ⁵⁶Fe/⁵⁴Fe ratios of poorly crystalline hydrous ferric oxide (HFO) metabolic products were enriched in the heavier isotope relative to aqueous ferrous iron (Fe[II]_aq) by ∼1.5 ± 0.2‰. This fractionation appears to be independent of the rate of photoautotrophic Fe(II)-oxidation, and is comparable to that observed for Fe isotope fractionation by dissimilatory Fe(III)-reducing bacteria. Although there remain a number of uncertainties regarding how the overall measured isotopic fractionation is produced, the most likely mechanisms include (1) an equilibrium effect produced by biological ligands, or (2) a kinetic effect produced by precipitation of HFO overlaid upon equilibrium exchange between Fe(II) and Fe(III) species. The fractionation we observe is similar in direction to that measured for abiotic oxidation of Fe(II)_aq by molecular oxygen. This suggests that the use of Fe isotopes to identify phototrophic Fe(II)-oxidation in the rock record may only be possible during time periods in Earth’s history when independent evidence exists for low ambient oxygen contents.     

First experimental determination of iron isotope fractionation between hematite and aqueous solution at hydrothermal conditions

Although iron isotopes provide a new powerful tool for tracing a variety of geochemical processes, the unambiguous interpretation of iron isotope ratios in natural systems and the development of predictive theoretical models require accurate data on equilibrium isotope fractionation between fluids and minerals. We investigated Fe isotope fractionation between hematite (Fe₂O₃) and aqueous acidic NaCl fluids via hematite dissolution and precipitation experiments at temperatures from 200 to 450 °C and pressures from saturated vapor pressure (P_sat) to 600 bar. Precipitation experiments at 200 °C and P_sat from aqueous solution, in which Fe aqueous speciation is dominated by ferric iron (Fe^III) chloride complexes, show no detectable Fe isotope fractionation between hematite and fluid, Δ⁵⁷Fe_{fluid-hematite} = δ⁵⁷Fe_fluid − δ⁵⁷Fe_hematite = 0.01 ± 0.08‰ (2 × standard error, 2SE). In contrast, experiments at 300 °C and P_sat, where ferrous iron chloride species (FeCl₂ and FeCl⁺) dominate in the fluid, yield significant fluid enrichment in the light isotope, with identical values of Δ⁵⁷Fe_{fluid-hematite} = −0.54 ± 0.15‰ (2SE) both for dissolution and precipitation runs. Hematite dissolution experiments at 450 °C and 600 bar, in which Fe speciation is also dominated by ferrous chloride species, yield Δ⁵⁷Fe_{fluid-hematite} values close to zero within errors, 0.15 ± 0.17‰ (2SE). In most experiments, chemical, redox, and isotopic equilibrium was attained, as shown by constancy over time of total dissolved Fe concentrations, aqueous Fe^II and Fe^III fractions, and Fe isotope ratios in solution, and identical Δ⁵⁷Fe values from dissolution and precipitation runs. Our measured equilibrium Δ⁵⁷Fe_{fluid-hematite} values at different temperatures, fluid compositions and iron redox state are within the range of fractionations in the system fluid-hematite estimated using reported theoretical β-factors for hematite and aqueous Fe species and the distribution of Fe aqueous complexes in solution. These theoretical predictions are however affected by large discrepancies among different studies, typically ±1‰ for the Δ⁵⁷Fe _{Fe(aq)-hematite} value at 200 °C. Our data may thus help to refine theoretical models for β-factors of aqueous iron species. This study provides the first experimental calibration of Fe isotope fractionation in the system hematite-saline aqueous fluid at elevated temperatures; it demonstrates the importance of redox control on Fe isotope fractionation at hydrothermal conditions.     

Iron and magnesium isotopic compositions of peridotite xenoliths from Eastern China

We present high-precision measurements of Mg and Fe isotopic compositions of olivine, orthopyroxene (opx), and clinopyroxene (cpx) for 18 lherzolite xenoliths from east central China and provide the first combined Fe and Mg isotopic study of the upper mantle. δ⁵⁶Fe in olivines varies from 0.18‰ to −0.22‰ with an average of −0.01 ± 0.18‰ (2SD, n = 18), opx from 0.24‰ to −0.22‰ with an average of 0.04 ± 0.20‰, and cpx from 0.24‰ to −0.16‰ with an average of 0.10 ± 0.19‰. δ²⁶Mg of olivines varies from −0.25‰ to −0.42‰ with an average of −0.34 ± 0.10‰ (2SD, n = 18), opx from −0.19‰ to −0.34‰ with an average of −0.25 ± 0.10‰, and cpx from −0.09‰ to −0.43‰ with an average of −0.24 ± 0.18‰. Although current precision (∼±0.06‰ for δ⁵⁶Fe; ±0.10‰ for δ²⁶Mg, 2SD) limits the ability to analytically distinguish inter-mineral isotopic fractionations, systematic behavior of inter-mineral fractionation for both Fe and Mg is statistically observed: Δ⁵⁶Fe_ol-cpx = −0.10 ± 0.12‰ (2SD, n = 18); Δ⁵⁶Fe_ol-opx = −0.05 ± 0.11‰; Δ²⁶Mg_ol-opx = −0.09 ± 0.12‰; Δ²⁶Mg_ol-cpx = −0.10 ± 0.15‰. Fe and Mg isotopic composition of bulk rocks were calculated based on the modes of olivine, opx, and cpx. The average δ⁵⁶Fe of peridotites in this study is 0.01 ± 0.17‰ (2SD, n = 18), similar to the values of chondrites but slightly lower than mid-ocean ridge basalts (MORB) and oceanic island basalts (OIB). The average δ²⁶Mg is −0.30 ± 0.09‰, indistinguishable from chondrites, MORB, and OIB. Our data support the conclusion that the bulk silicate Earth (BSE) has chondritic δ⁵⁶Fe and δ²⁶Mg.The origin of inter-mineral fractionations of Fe and Mg isotopic ratios remains debated. δ⁵⁶Fe between the main peridotite minerals shows positive linear correlations with slopes within error of unity, strongly suggesting intra-sample mineral-mineral Fe and Mg isotopic equilibrium. Because inter-mineral isotopic equilibrium should be reached earlier than major element equilibrium via chemical diffusion at mantle temperatures, Fe and Mg isotope ratios of coexisting minerals could be useful tools for justifying mineral thermometry and barometry on the basis of chemical equilibrium between minerals. Although most peridotites in this study exhibit a narrow range in δ⁵⁶Fe, the larger deviations from average δ⁵⁶Fe for three samples likely indicate changes due to metasomatic processes. Two samples show heavy δ⁵⁶Fe relative to the average and they also have high La/Yb and total Fe content, consistent with metasomatic reaction between peridotite and Fe-rich and isotopically heavy melt. The other sample has light δ⁵⁶Fe and slightly heavy δ²⁶Mg, which may reflect Fe-Mg inter-diffusion between peridotite and percolating melt.     

Copper isotope fractionation in acid mine drainage

We measured the Cu isotopic composition of primary minerals and stream water affected by acid mine drainage in a mineralized watershed (Colorado, USA). The δ⁶⁵Cu values (based on ⁶⁵Cu/⁶³Cu) of enargite (δ⁶⁵Cu = −0.01 ± 0.10‰; 2σ) and chalcopyrite (δ⁶⁵Cu = 0.16 ± 0.10‰) are within the range of reported values for terrestrial primary Cu sulfides (−1‰ < δ⁶⁵Cu < 1‰). These mineral samples show lower δ⁶⁵Cu values than stream waters (1.38‰ ? δ⁶⁵Cu ? 1.69‰). The average isotopic fractionation (Δ_aq-min = δ⁶⁵Cu_aq − δ⁶⁵Cu_min, where the latter is measured on mineral samples from the field system), equals 1.43 ± 0.14‰ and 1.60 ± 0.14‰ for chalcopyrite and enargite, respectively. To interpret this field survey, we leached chalcopyrite and enargite in batch experiments and found that, as in the field, the leachate is enriched in ⁶⁵Cu relative to chalcopyrite (1.37 ± 0.14‰) and enargite (0.98 ± 0.14‰) when microorganisms are absent. Leaching of minerals in the presence of Acidithiobacillus ferrooxidans results in smaller average fractionation in the opposite direction for chalcopyrite (Δ_aq-min^o=-0.57±0.14‰, where min^o refers to the starting mineral) and no apparent fractionation for enargite (Δ_aq-min^o=0.14±0.14‰). Abiotic fractionation is attributed to preferential oxidation of ⁶⁵Cu⁺ at the interface of the isotopically homogeneous mineral and the surface oxidized layer, followed by solubilization. When microorganisms are present, the abiotic fractionation is most likely not seen due to preferential association of ⁶⁵Cu_aq with A. ferrooxidans cells and related precipitates. In the biotic experiments, Cu was observed under TEM to occur in precipitates around bacteria and in intracellular polyphosphate granules. Thus, the values of δ⁶⁵Cu in the field and laboratory systems are presumably determined by the balance of Cu released abiotically and Cu that interacts with cells and related precipitates. Such isotopic signatures resulting from Cu sulfide dissolution should be useful for acid mine drainage remediation and ore prospecting purposes.     

Copper isotopic zonation in the Northparkes porphyry Cu-Au deposit, SE Australia

Significant, systematic Cu isotopic variations have been found in the Northparkes porphyry Cu-Au deposit, NSW, Australia, which is an orthomagmatic porphyry Cu deposit. Copper isotope ratios have been measured in sulfide minerals (chalcopyrite and bornite) by both solution and laser ablation multi-collector inductively coupled plasma-mass spectrometry (MC-ICP-MS). The results from both methods show a variation in δ⁶⁵Cu of hypogene sulfide minerals of greater than 1‰ (relative to NIST976). Significantly, the results from four drill holes through two separate ore bodies show strikingly similar patterns of Cu isotope variation. The patterns are characterized by a sharp down-hole decrease from up to 0.8‰ (0.29 ± 0.56‰, 1σ, n = 20) in the low-grade peripheral alteration zones (phyllic-propylitic alteration zone) to a low of ∼−0.4‰ (−0.25 ± 0.36‰, 1σ, n = 30) at the margins of the most mineralized zones (Cu grade >1 wt%). In the high-grade cores of the systems, the compositions are more consistent at around 0.2‰ (0.19 ± 0.14‰, 1σ, n = 40). The Cu isotopic zonation may be explained by isotope fractionation of Cu between vapor, solution and sulfides at high temperature, during boiling and sulfide precipitation processes. Sulfur isotopes also show an isotopically light shell at the margins of the high-grade ore zones, but these are displaced from the low δ⁶⁵Cu shells, such that there is no correlation between the Cu and S isotope signatures. Fe isotope data do not show any discernable variation along the drill core. This work demonstrates that Cu isotopes show a large response to high-temperature porphyry mineralizing processes, and that they may act as a vector to buried mineralization.     

The experimental calibration of the iron isotope fractionation factor between pyrrhotite and peralkaline rhyolitic melt

A first experimental study was conducted to determine the equilibrium iron isotope fractionation between pyrrhotite and silicate melt at magmatic conditions. Experiments were performed in an internally heated gas pressure vessel at 500 MPa and temperatures between 840 and 1000 °C for 120-168 h. Three different types of experiments were conducted and after phase separation the iron isotope composition of the run products was measured by MC-ICP-MS. (i) Kinetic experiments using ⁵⁷Fe-enriched glass and natural pyrrhotite revealed that a close approach to equilibrium is attained already after 48 h. (ii) Isotope exchange experiments—using mixtures of hydrous peralkaline rhyolitic glass powder (∼4 wt% H₂O) and natural pyrrhotites (Fe_1 − xS) as starting materials— and (iii) crystallisation experiments, in which pyrrhotite was formed by reaction between elemental sulphur and rhyolitic melt, consistently showed that pyrrhotite preferentially incorporates light iron. No temperature dependence of the fractionation factor was found between 840 and 1000 °C, within experimental and analytical precision. An average fractionation factor of Δ ⁵⁶Fe/⁵⁴Fe_{pyrrhotite-melt} = −0. 35 ± 0.04‰ (2SE, n = 13) was determined for this temperature range. Predictions of Fe isotope fractionation between FeS and ferric iron-dominated silicate minerals are consistent with our experimental results, indicating that the marked contrast in both ligand and redox state of iron control the isotope fractionation between pyrrhotite and silicate melt. Consequently, the fractionation factor determined in this study is representative for the specific Fe²⁺/ΣFe ratio of our peralkaline rhyolitic melt of 0.38 ± 0.02. At higher Fe²⁺/ΣFe ratios a smaller fractionation factor is expected. Further investigation on Fe isotope fractionation between other mineral phases and silicate melts is needed, but the presented experimental results already suggest that even at high temperatures resolvable variations in the Fe isotope composition can be generated by equilibrium isotope fractionation in natural magmatic systems.     

Fractionation of Cu and Zn isotopes during adsorption onto amorphous Fe(III) oxyhydroxide: Experimental mixing of acid rock drainage and ambient river water

Fractionation of Cu and Zn isotopes during adsorption onto amorphous ferric oxyhydroxide is examined in experimental mixtures of metal-rich acid rock drainage and relatively pure river water and during batch adsorption experiments using synthetic ferrihydrite. A diverse set of Cu- and Zn-bearing solutions was examined, including natural waters, complex synthetic acid rock drainage, and simple NaNO₃ electrolyte. Metal adsorption data are combined with isotopic measurements of dissolved Cu (⁶⁵Cu/⁶³Cu) and Zn (⁶⁶Zn/⁶⁴Zn) in each of the experiments. Fractionation of Cu and Zn isotopes occurs during adsorption of the metal onto amorphous ferric oxyhydroxide. The adsorption data are modeled successfully using the diffuse double layer model in PHREEQC. The isotopic data are best described by a closed system, equilibrium exchange model. The fractionation factors (α_soln-solid) are 0.99927 ± 0.00008 for Cu and 0.99948 ± 0.00004 for Zn or, alternately, the separation factors (Δ_soln-solid) are −0.73 ± 0.08‰ for Cu and −0.52 ± 0.04‰ for Zn. These factors indicate that the heavier isotope preferentially adsorbs onto the oxyhydroxide surface, which is consistent with shorter metal-oxygen bonds and lower coordination number for the metal at the surface relative to the aqueous ion. Fractionation of Cu isotopes also is greater than that for Zn isotopes. Limited isotopic data for adsorption of Cu, Fe(II), and Zn onto amorphous ferric oxyhydroxide suggest that isotopic fractionation is related to the intrinsic equilibrium constants that define aqueous metal interactions with oxyhydroxide surface sites. Greater isotopic fractionation occurs with stronger metal binding by the oxyhydroxide with Cu > Zn > Fe(II).     

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.     

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.     

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.     

The effect of early diagenesis on the Fe isotope compositions of porewaters and authigenic minerals in continental margin sediments

Iron isotope compositions in marine pore fluids and sedimentary solid phases were measured at two sites along the California continental margin, where isotope compositions range from δ⁵⁶Fe = −3.0‰ to +0.4‰. At one site near Monterey Canyon off central California, organic matter oxidation likely proceeds through a number of diagenetic pathways that include significant dissimilatory iron reduction (DIR) and bacterial sulfate reduction, whereas at our other site in the Santa Barbara basin DIR appears to be comparatively small, and production of sulfides (FeS and pyrite) was extensive. The largest range in Fe isotope compositions is observed for Fe(II)_aq in porewaters, which generally have the lowest δ⁵⁶Fe values (minimum: −3.0‰) near the sediment surface, and increase with burial depth. δ⁵⁶Fe values for FeS inferred from HCl extractions vary between ∼−0.4‰ and +0.4‰, but pyrite is similar at both stations, where an average δ⁵⁶Fe value of −0.8 ± 0.2‰ was measured. We interpret variations in dissolved Fe isotope compositions to be best explained by open-system behavior that involves extensive recycling of Feflux. This study is the first to examine Fe isotope variations in modern marine sediments, and the results show that Fe isotopes in the various reactive Fe pools undergo isotopic fractionation during early diagenesis. Importantly, processes dominated by sulfide formation produce high-δ⁵⁶Fe values for porewaters, whereas the opposite occurs when Fe(III)-oxides are present and DIR is a major pathway of organic carbon respiration. Because shelf pore fluids may carry a negative δ⁵⁶Fe signature it is possible that the Fe isotope composition of ocean water reflects a significant contribution of shelf-derived iron to the open ocean. Such a signature would be an important means for tracing iron sources to the ocean and water mass circulation.     

Calcium isotope (δCa) fractionation along hydrothermal pathways, Logatchev field (Mid-Atlantic Ridge, 14°45′N)

We investigate the Logatchev Hydrothermal Field at the Mid-Atlantic Ridge, 14°45′N to constrain the calcium isotope hydrothermal flux into the ocean. During the transformation of seawater to a hydrothermal solution, the Ca concentration of pristine seawater ([Ca]_SW) increases from about 10 mM to about 32 mM in the hydrothermal fluid endmember ([Ca]_HydEnd) and thereby adopts a δ^44/40Ca_HydEnd of −0.95 ± 0.07‰ relative to seawater (SW) and a ⁸⁷Sr/⁸⁶Sr isotope ratio of 0.7034(4). We demonstrate that δ^44/40Ca_HydEnd is higher than that of the bedrock at the Logatchev field. From mass balance calculations, we deduce a δ^44/40Ca of −1.17 ± 0.04‰ (SW) for the host-rocks in the reaction zone and −1.45 ± 0.05‰ (SW) for the isotopic composition of the entire hydrothermal cell of the Logatchev field. The values are isotopically lighter than the currently assumed δ^44/40Ca for Bulk Earth of −0.92 ± 0.18‰ (SW) [Skulan J., DePaolo D. J. and Owens T. L. (1997) Biological control of calcium isotopic abundances in the global calcium cycle. Geochim. Cosmochim. Acta61,(12) 2505-2510] and challenge previous assumptions of no Ca isotope fractionation between hydrothermal fluid and the oceanic crust [Zhu P. and Macdougall J. D. (1998) Calcium isotopes in the marine environment and the oceanic calcium cycle. Geochim. Cosmochim. Acta62,(10) 1691-1698; Schmitt A. -D., Chabeaux F. and Stille P. (2003) The calcium riverine and hydrothermal isotopic fluxes and the oceanic calcium mass balance. Earth Planet. Sci. Lett. 6731, 1-16]. Here we propose that Ca isotope fractionation along the fluid flow pathway of the Logatchev field occurs during the precipitation of anhydrite. Two anhydrite samples from the Logatchev Hydrothermal Field show an average fractionation of about Δ^44/40Ca = −0.5‰ relative to their assumed parental solutions. Ca isotope ratios in aragonites from carbonate veins from ODP drill cores indicate aragonite precipitation directly from seawater at low temperatures with an average δ^44/40Ca of −1.54 ± 0.08‰ (SW). The relatively large fractionation between the aragonite precipitates and seawater in combination with their frequent abundance in weathered mafic and ultramafic rocks suggest a reconsideration of the marine Ca isotope budget, in particular with regard to ocean crust alteration.     

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.