Iron isotope fractionation during hydrothermal ore deposition and alteration

Iron isotopes fractionate during hydrothermal processes. Therefore, the Fe isotope composition of ore-forming minerals characterizes either iron sources or fluid histories. The former potentially serves to distinguish between sedimentary, magmatic or metamorphic iron sources, and the latter allows the reconstruction of precipitation and redox processes. These processes take place during ore formation or alteration. The aim of this contribution is to investigate the suitability of this new isotope method as a probe of ore-related processes. For this purpose 51 samples of iron ores and iron mineral separates from the Schwarzwald region, southwest Germany, were analyzed for their iron isotope composition using multicollector ICP-MS. Further, the ore-forming and ore-altering processes were quantitatively modeled using reaction path calculations. The Schwarzwald mining district hosts mineralizations that formed discontinuously over almost 300 Ma of hydrothermal activity. Primary hematite, siderite and sulfides formed from mixing of meteoric fluids with deeper crustal brines. Later, these minerals were partly dissolved and oxidized, and secondary hematite, goethite and iron arsenates were precipitated. Two types of alteration products formed: (1) primary and high-temperature secondary Fe minerals formed between 120 and 300 °C, and (2) low-temperature secondary Fe minerals formed under supergene conditions (<100 °C). Measured iron isotope compositions are variable and cover a range in δ⁵⁶Fe between −2.3‰ and +1.3‰. Primary hematite (δ⁵⁶Fe: −0.5‰ to +0.5‰) precipitated by mixing oxidizing surface waters with a hydrothermal fluid that contained moderately light Fe (δ⁵⁶Fe: −0.5‰) leached from the crystalline basement. Occasional input of CO₂-rich waters resulted in precipitation of isotopically light siderite (δ⁵⁶Fe: −1.4 to −0.7‰). The difference between hematite and siderite is compatible with published Fe isotope fractionation factors. The observed range in isotopic compositions can be accounted for by variable fractions of Fe precipitating from the fluid. Therefore, both fluid processes and mass balance can be inferred from Fe isotopes. Supergene weathering of siderite by oxidizing surface waters led to replacement of isotopically light primary siderite by similarly light secondary hematite and goethite, respectively. Because this replacement entails quantitative transfer of iron from precursor mineral to product, no significant isotope fractionation is produced. Hence, Fe isotopes potentially serve to identify precursors in ore alteration products. Goethites from oolitic sedimentary iron ores were also analyzed. Their compositional range appears to indicate oxidative precipitation from relatively uniform Fe dissolved in coastal water. This comprehensive iron isotope study illustrates the potential of the new technique in deciphering ore formation and alteration processes. Isotope ratios are strongly dependent on and highly characteristic of fluid and precipitation histories. Therefore, they are less suitable to provide information on Fe sources. However, it will be possible to unravel the physico-chemical processes leading to the formation, dissolution and redeposition of ores in great detail.     

Iron isotopes in natural carbonate minerals determined by MC-ICP-MS with a Fe-Fe double spike

We have developed a method for iron isotope analysis by multiple-collector inductively coupled plasma mass spectrometry (MC-ICP-MS) using a ⁵⁸Fe-⁵⁴Fe double spike. A 20 min analysis produces mass-bias-corrected iron isotope data with an external reproducibility of ±0.05 (2 SD) on δ⁵⁶Fe, which represents a decrease in analysis time compared to sample-standard bracketing techniques. The estimation of external reproducibility is based on replicate analysis of the ETH hematite in-house standard. The double spike method has two advantages. First, matrix effects during MC-ICP-MS analysis are decreased with tests showing that accurate iron isotope data can, in some cases, be obtained even when matrix levels exceed iron concentration (Na/Fe, Mg/Fe, and Ca/Fe up to 5, 2, and 0.1, respectively). Because chemical separation reduces matrix/Fe to levels more than three orders of magnitude lower than this, measured Fe isotope compositions are unlikely to be compromised by matrix effects. Second, it is possible to spike samples before chemical purification, which enables any isotopic fractionation effect because of incomplete recovery of iron from a sample to be accounted for. This may be important where obtaining quantitative iron yields from samples is difficult, such as the extraction of dissolved iron from water samples. Fe isotope data on a set of standard reference materials (igneous rocks, ferromanganese nodules, sedimentary rocks, and ores) are presented, which are in agreement with previously published data considering analytical uncertainties. Mantle-derived standard rock samples that are the source of iron for surficial, (bio)geochemical cycling yield a mean δ⁵⁶Fe of 0.041 ± 0.11‰ (n = 8; 2 SD) with reference to IRMM-14. Hydrothermal and metamorphic calcium carbonate rocks with a relatively low iron content (100-4000 ppm) have δ⁵⁶Fe = −1.25 to −0.07‰. Structural Fe(II) in hydrothermal calcites has δ⁵⁶Fe = −1.25 to −0.27‰. The light iron in this range of carbonate minerals may reflect the iron isotope composition of the hydrothermal fluids from which the carbonate precipitated, or the presence of Fe(III) and/or organic material in the hydrothermal fluids during calcite precipitation.     

Deciphering formation processes of banded iron formations from the Transvaal and the Hamersley successions by combined Si and Fe isotope analysis using UV femtosecond laser ablation

To investigate the genesis of BIFs, we have determined the Fe and Si isotope composition of coexisting mineral phases in samples from the ∼2.5 billion year old Kuruman Iron Formation (Transvaal Supergroup, South Africa) and Dales Gorges Member of the Brockman Iron Formation (Hamersley Group, Australia) by UV femtosecond laser ablation coupled to a MC-ICP-MS. Chert yields a total range of δ³⁰Si between −1.3‰ and −0.8‰, but the Si isotope compositions are uniform in each core section examined. This uniformity suggests that Si precipitated from well-mixed seawater far removed from its sources such as hydrothermal vents or continental drainage. The Fe isotope composition of Fe-bearing mineral phases is much more heterogeneous compared to Si with δ⁵⁶Fe values of −2.2‰ to 0‰. This heterogeneity is likely due to variable degrees of partial Fe(II) oxidation in surface waters, precipitation of different mineral phases and post-depositional Fe redistribution. Magnetite exhibits negative δ⁵⁶Fe values, which can be attributed to a variety of diagenetic pathways: the light Fe isotope composition was inherited from the Fe(III) precursor, heavy Fe(II) was lost by abiotic reduction of the Fe(III) precursor or light Fe(II) was gained from external fluids. Micrometer-scale heterogeneities of δ⁵⁶Fe in Fe oxides are attributed to variable degrees of Fe(II) oxidation or to isotope exchange upon Fe(II) adsorption within the water column and to Fe redistribution during diagenesis. Diagenetic Fe(III) reduction caused by oxidation of organic matter and Fe redistribution is supported by the C isotope composition of a carbonate-rich sample containing primary siderite. These carbonates yield δ¹³C values of ∼−10‰, which hints at a mixed carbon source in the seawater of both organic and inorganic carbon. The ancient seawater composition is estimated to have a minimum range in δ⁵⁶Fe of −0.8‰ to 0‰, assuming that hematite and siderite have preserved their primary Fe isotope signature. The long-term near-zero Fe isotope composition of the Hamersley and Transvaal BIFs is in balance with the assumed composition of the Fe sources. The negative Fe isotope composition of the investigated BIF samples, however, indicates either a perturbation of the steady state, or they have to be balanced spatially by deposition of isotopically heavy Fe. In the case of Si, the negative Si isotope signature of these BIFs stands in marked contrast to the assumed source composition. The deviation from potential source composition requires a complementary sink of isotopically heavy Si in order to maintain steady state in the basin. Perturbing the steady state by extraordinary hydrothermal activity or continental weathering in contrast would have led to precipitation of light Si isotopes from seawater. Combining an explanation for both elements, a likely scenario is a steady state ocean basin with two sinks. When all published Fe isotope records including BIFs, microbial carbonates, shales and sedimentary pyrites, are considered, a complementary sink for heavy Fe isotopes must have existed in Precambrian ocean basins. This Fe sink could have been pelagic sediments, which however are not preserved. For Si, such a complementary sink for heavy Si isotopes might have been provided by other chert deposits within the basin.     

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.     

Iron isotope fractionation during skarn-type metallogeny: A case study of Xinqiao Cu–S–Fe–Au deposit in the Middle–Lower Yangtze valley

Fe isotope compositions of mineral separates and bulk samples from Xinqiao Cu–S–Fe–Au skarn type deposit were investigated. An overall variation in δ⁵⁷Fe values from − 1.22‰ to + 0.73‰ has been observed, which shows some regularity. The δ⁵⁷Fe values of endoskarn and the earliest formed Fe-mineral phase magnetite are ca.1.2‰ and ca. 0.3‰ lower, respectively, relative to the quartz–monzodiorite stock, indicating that fluid exsolved from the stock is enriched in light Fe isotopes. Moreover, spatial and temporal variations in δ⁵⁷Fe values are observed, which suggest iron isotope fractionation during fluid evolution. Precipitation of Fe-bearing minerals results in the Fe isotope composition of residual fluids evolving with time. Precipitation of Fe (III) minerals incorporating heavy iron isotopes preferentially leaves the remaining fluid enriched in light isotopes, while precipitation of Fe (II) minerals preferentially taking-up light iron isotopes, and makes the Fe isotopic composition of the fluid progressively heavier. The regularity of Fe isotope variations occurred during fluid exsolution and evolution indicates that the dominant Fe source of Xinqiao deposit is magmatic. Overall, this study demonstrates that Fe isotope composition has great potential in unraveling ore-forming processes, as well as constraining the metal sources of ore deposits.     

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.     

Accurate measurement of silver isotopic compositions in geological materials including low Pd/Ag meteorites

Very precise silver (Ag) isotopic compositions have been determined for a number of terrestrial rocks, and high and low Pd/Ag meteorites by utilizing multicollector inductively coupled plasma mass spectrometry (MC-ICPMS). The meteorites include primitive chondrites, the Group IAB iron meteorites Canyon Diablo and Toluca, and the Group IIIAB iron meteorite Grant. Silver isotopic measurements are primarily of interest because ¹⁰⁷Ag was produced by decay of the short-lived radionuclide ¹⁰⁷Pd during the formation of the solar system and hence the Pd-Ag chronometer has set constraints on the timing of early planetesimal formation. A 2σ precision of ±0.05‰ can be obtained for analyses of standard solutions when Ag isotopic ratios are normalized to Pd, to correct for instrumental mass discrimination, and to bracketing standards. Caution must be exercised when making Ag isotopic measurements because isotopic artifacts can be generated in the laboratory and during mass spectrometry. The external reproducibility for geological samples based on replicate analyses of rocks is ±0.2‰ (2σ).All chondrites analyzed have similar Ag isotopic compositions that do not differ significantly (>0.3‰) from the ‘terrestrial’ value of the NIST SRM 978a Ag isotope standard. Hence, they show no evidence of excess ¹⁰⁷Ag derived from ¹⁰⁷Pd decay or, of stable Ag isotope fractionation associated with volatile element depletion within the accretion disk or from parent body metamorphism. The Group IAB iron meteorite samples analyzed show evidence of complex behavior and disturbance of Ag isotope systematics. Therefore, care must be taken when using this group of iron meteorites to obtain chronological information based on the Pd-Ag decay scheme.     

Grain-scale iron isotopic distribution of pyrite from Precambrian shallow marine carbonate revealed by a femtosecond laser ablation multicollector ICP-MS technique: Possible proxy for the redox state of ancient seawater

The redox state of Precambrian shallow seas has been linked with material cycle and evolution of the photosynthesis-based ecosystem. Iron is a redox-sensitive element and exists as a soluble Fe(II) species or insoluble Fe(III) species on Earth’s surface. Previous studies have shown that the iron isotopic ratio of marine sedimentary minerals is useful for understanding the ocean redox state, although the redox state of the Archean shallow sea is poorly known. This is partly because the conventional bulk isotope analytical technique has often been used, wherein the iron isotopic record may be dampened by the presence of isotopically different iron-bearing minerals within the same sample. Here we report a microscale iron isotopic ratio of individual pyrite grains in shallow marine stromatolitic carbonates over geological time using a newly developed, near-infrared femtosecond laser ablation multicollector ICP-MS technique (NIR-fs-LA-MC-ICP-MS).We have determined that the grain-scale iron isotopic distribution of pyrite from coeval samples shows a bimodal (2.7 and 2.3 Ga) or unimodal pattern (2.9, 2.6, and 0.7 Ga). In particular, pyrite from the 2.7 Ga Fortescue Group shows a unique bimodal distribution with highly positive (+1.0‰ defined as Type 1) and negative δ⁵⁶Fe values (−1.8‰ defined as Type 2). Type 1 and 2 pyrites occasionally occur within different siliceous layers in the same rock specimen. Layer-scale iron isotopic heterogeneity indicates that the iron isotopic ratios of the two types of pyrite are not homogenized by diagenesis after deposition. Some cubic pyrites have a core with a positive δ⁵⁶Fe value (1‰) and a rim with a crustal δ⁵⁶Fe value (0‰). The observed isotopic zoning suggests that the positive δ⁵⁶Fe value is a primary signature at the time of stromatolite formation, while secondary pyrite precipitated during diagenesis.The positive δ⁵⁶Fe value of Type 1 and the large iron isotopic difference between Type 1 and 2 (2.8‰.) suggest partial Fe(II) oxidation in the 2.7-Ga shallow sea, i.e., pyritization of ⁵⁶Fe-enriched ferric oxyhydroxide (Type 1) and ⁵⁶Fe depleted Fe²⁺aq in seawater (Type 2). Type 2 pyrite was probably not produced by microbial iron redox cycling during diagenesis because this scenario requires a higher abundance of pyrite with δ⁵⁶Fe of 0‰ than of −1.8‰. Consequently, the degree of Fe(II) oxidation in the 2.7-Ga shallow sea can be estimated by a Fe²⁺aq steady-state model. The model calculation shows that half the Fe²⁺aq influx was oxidized in the seawater. This implies that O₂ produced by photosynthesis would have been completely consumed by oxidation of the Fe²⁺aq influx. Grain-scale iron isotopic distribution of pyrite could be a useful index for reconstructing the redox state of the Archean shallow sea.     

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.     

High precision iron isotope measurements of meteoritic material by cold plasma ICP-MS

The first cold plasma ICP-MS (inductively coupled plasma mass spectrometer) Fe isotope study is described. Application of this technique to the analyses of Fe isotopes in a number of meteorites is also reported. The measurement technique relies on reduced temperature operation of the ICP source to eliminate pervasive molecular interferences from Ar complexes associated with conventional ICP-MS. Instrumental mass bias corrections are performed by sample-standard bracketing and using Cu as an external mass bias drift monitor. Repeated measurements of a terrestrial basalt reference sample indicate an external reproducibility of ± 0.06 ‰ for δ⁵⁶Fe and ± 0.25 ‰ for δ⁵⁸Fe (1 σ). The measured iron isotopic compositions of various bulk meteorites, including irons, chondrites and pallasites are identical, within error, to the composition of our terrestrial basalt reference sample suggesting that iron mass fractionation during planet formation and differentiation was non-existent. Iron isotope compositions measured for eight chondrules from the unequilibrated ordinary chondrite Tieschitz range from −0.5 ‰ < δ⁵⁶Fe_chondrules < 0.0 ‰ relative to the terrestrial/meteorite average. Mechanisms for fractionating iron in these chondrules are discussed.     

Micro-scale tracing of Fe and Si isotope signatures in banded iron formation using femtosecond laser ablation

We have detected micrometre-scale differences in Fe and Si stable isotope ratios between coexisting minerals and between layers of banded iron formation (BIF) using an UV femtosecond laser ablation system connected to a MC-ICP-MS. In the magnetite–carbonate–chert BIF from the Archean Old Wanderer Formation in the Shurugwi Greenstone Belt (Zimbabwe), magnetite shows neither intra- nor inter-layer trends giving overall uniform δ⁵⁶Fe values of 0.9‰, but exhibits intra-crystal zonation. Bulk iron carbonates are also relatively uniform at near-zero values, however, their individual δ⁵⁶Fe value is highly composition-dependent: both siderite and ankerite and mixtures between both are present, and δ⁵⁶Fe end member values are 0.4‰ for siderite and −0.7‰ for ankerite. The data suggest either an early diagenetic origin of magnetite and iron carbonates by the reaction of organic matter with ferric oxyhydroxides catalysed by Fe(III)-reducing bacteria; or more likely an abiotic reaction of organic carbon and Fe(III) during low-grade metamorphism. Si isotope composition of the Old Wanderer BIF also shows significant variations with δ³⁰Si values that range between −1.0‰ and −2.6‰ for bulk layers. These isotope compositions suggest rapid precipitation of the silicate phases from hydrothermal-rich waters. Interestingly, Fe and Si isotope compositions of bulk layers are covariant and are interpreted as largely primary signatures. Moreover, the changes of Fe and Si isotope signatures between bulk layers directly reflect the upwelling dynamics of hydrothermal-rich water which govern the rates of Fe and Si precipitation and therefore also the development of layering. During periods of low hydrothermal activity, precipitation of only small amounts of ferric oxyhydroxide was followed by complete reduction with organic carbon during diagenesis resulting in carbonate–chert layers. During periods of intensive hydrothermal activity, precipitation rates of ferric oxyhydroxide were high, and subsequent diagenesis triggered only partial reduction, forming magnetite–carbonate–chert layers. We are confident that our micro-analytical technique is able to detect both the solute flux history into the sedimentary BIF precursor, and the BIF’s diagenetic history from the comparison between coexisting minerals and their predicted fractionation factors.     

Fe isotopic fractionation during mineral dissolution with and without bacteria

Fe released into solution is isotopically lighter (enriched in the lighter isotope) than hornblende starting material when dissolution occurs in the presence of the siderophore desferrioxamine mesylate (DFAM). In contrast, Fe released from goethite dissolving in the presence of DFAM is isotopically unchanged. Furthermore, Δ⁵⁶Fe_{solution-hornblende} for Fe released to solution in the presence of ligands varies with the affinity of the ligand for Fe. The extent of isotopic fractionation of Fe released from hornblende also increases when experiments are agitated continuously. The Fe isotope fractionation observed during hornblende dissolution with organic ligands is attributed predominantly to retention of ⁵⁶Fe in an altered surface layer, while the lack of isotopic fractionation during goethite dissolution in DFAM is consistent with the lack of an altered layer. When a siderophore-producing soil bacterium is added to the system (without added organic ligands), Fe released to solution from both hornblende and goethite differs isotopically from Fe in the bulk mineral: Δ⁵⁶Fe_{solution-starting material} = −0.56 ± 0.19 (hornblende) and −1.44 ± 0.16 (goethite). Increased isotopic fractionation is attributed in this case to the fact that as bacterial respiration depletes the system in oxygen and aqueous Fe is reduced, equilibration between aqueous ferrous and ferric iron creates a pool of isotopically heavy ferric iron that is assimilated by bacterial cells. Adsorption of isotopically heavy ferrous iron (Fe(II) enriched in the heavier isotope) or precipitation of isotopically heavy Fe minerals may also contribute to observed fractionations.To test whether these Fe isotope signatures are recorded in natural systems, we also investigated extractions of samples of soils from which the bacteria were isolated. These extractions show variability in the isotopic signatures of exchangeable Fe and Fe oxyhydroxide fractions from one soil sample to another, but exchangeable Fe is observed to be lighter than Fe in soil Fe oxyhydroxides and hornblende. This observation is consistent with isotopically light Fe-organic complexes in soil pore water derived from the Fe-silicate starting materials in the presence of growing microorganisms, as documented in experiments reported here. The contributions from phenomena including organic ligand-promoted nonstoichiometric dissolution of Fe silicates, uptake of ferric iron by organisms, adsorption of isotopically heavy ferrous iron, and precipitation of iron minerals should create complex isotopic signatures in soils. Better understanding of these processes and the timescales over which they contribute to fractionation is needed.     

Iron isotopes constrain biologic and abiologic processes in banded iron formation genesis

The voluminous 2.5 Ga banded iron formations (BIFs) from the Hamersley Basin (Australia) and Transvaal Craton (South Africa) record an extensive period of Fe redox cycling. The major Fe-bearing minerals in the Hamersley-Transvaal BIFs, magnetite and siderite, did not form in Fe isotope equilibrium, but instead reflect distinct formation pathways. The near-zero average δ⁵⁶Fe values for magnetite record a strong inheritance from Fe³⁺ oxide/hydroxide precursors that formed in the upper water column through complete or near-complete oxidation. Transformation of the Fe³⁺ oxide/hydroxide precursors to magnetite occurred through several diagenetic processes that produced a range of δ⁵⁶Fe values: (1) addition of marine hydrothermal , (2) complete reduction by bacterial dissimilatory iron reduction (DIR), and (3) interaction with excess that had low δ⁵⁶Fe values and was produced by DIR. Most siderite has slightly negative δ⁵⁶Fe values of ∼ −0.5‰ that indicate equilibrium with Late Archean seawater, although some very negative δ⁵⁶Fe values may record DIR. Support for an important role of DIR in siderite formation in BIFs comes from previously published C isotope data on siderite, which may be explained as a mixture of C from bacterial and seawater sources.Several factors likely contributed to the important role that DIR played in BIF formation, including high rates of ferric oxide/hydroxide formation in the upper water column, delivery of organic carbon produced by photosynthesis, and low clastic input. We infer that DIR-driven Fe redox cycling was much more important at this time than in modern marine systems. The low pyrite contents of magnetite- and siderite-facies BIFs suggests that bacterial sulfate reduction was minor, at least in the environments of BIF formation, and the absence of sulfide was important in preserving magnetite and siderite in the BIFs, minerals that are poorly preserved in the modern marine record. The paucity of negative δ⁵⁶Fe values in older (Early Archean) and younger (Early Proterozoic) BIFs suggests that the extensive 2.5 Ga Hamersley-Transvaal BIFs may record a period of maximum expansion of DIR in Earth’s history.     

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.     

安徽新桥矿床矿相学与Fe同位素特征及其对矿床成因的制约

对安徽新桥矿床进行系统的野外地质调查和矿相学研究发现,层状矿体中的胶状黄铁矿交代矽卡岩磁铁矿矿体,为探讨层状硫化物矿床是早期沉积成因还是岩浆热液成因提供了新的地质约束。对铜陵矿集区内的新桥矿床层状菱铁矿矿体和凤凰山矽卡岩型矿体中的菱铁矿开展了Fe同位素组成的对比研究,结果显示：新桥矿床菱铁矿与典型低温热液脉型矿床和沉积铁矿中的菱铁矿在Fe同位素组成特征上有所不同,而与凤凰山矽卡岩型矿床中的菱铁矿更为接近;新桥矿床中胶状黄铁矿和菱铁矿相对于磁铁矿富集Fe的轻同位素,表明磁铁矿不是过去认为的由胶状黄铁矿和菱铁矿矿胚层经热液改造形成,而是与典型的岩浆热液有关。新桥矿区层状硫化物矿体和矽卡岩型矿体中,近岩体矽卡岩和最早形成的金属矿物磁铁矿比岩体更为富集Fe的轻同位素,而赋矿围岩比岩体更为富集Fe的重同位素。同时,不同矿化阶段形成的含铁矿物和不同空间位置的硫化物中的Fe同位素组成呈现出时空分带现象,Fe同位素组成的时空演化特征与流体出溶、流体演化非常一致,并且符合同位素分馏的基本理论,表明层状硫化物矿体和矽卡岩型矿体具有相同的成矿物质来源,为同一流体体系演化的产物。新桥矿区岩相学的研究结果和Fe同位素组成特征均表明,新桥层状硫化物矿床不是海西期喷流沉积成矿作用的产物,而是燕山期热液成矿作用的产物,为一个典型的热液成因矿床。     

Iron isotope systematics in estuaries: The case of North River, Massachusetts (USA)

Recent studies have suggested that rivers may present an isotopically light Fe source to the oceans. Since the input of dissolved iron from river water is generally controlled by flocculation processes that occur during estuarine mixing, it is important to investigate potential fractionation of Fe-isotopes during this process. In this study, we investigate the influence of the flocculation of Fe-rich colloids on the iron isotope composition of pristine estuarine waters and suspended particles. The samples were collected along a salinity gradient from the fresh water to the ocean in the North River estuary (MA, USA). Estuarine samples were filtered at 0.22 μm and the iron isotope composition of the two fractions (dissolved and particles) were analyzed using high-resolution MC-ICP-MS after chemical purification. Dissolved iron results show positive δ⁵⁶Fe values (with an average of 0.43 ± 0.04‰) relative to the IRMM-14 standard and do not display any relationships with salinity or with percentage of colloid flocculation. The iron isotopic composition of the particles suspended in fresh water is characterized by more negative δ⁵⁶Fe values than for dissolved Fe and correlate with the percentage of Fe flocculation. Particulate δ⁵⁶Fe values vary from −0.09‰ at no flocculation to ∼0.1‰ at the flocculation maximum, which reflect mixing effects between river-borne particles, lithogenic particles derived from coastal seawaters and newly precipitated colloids. Since the process of flocculation produces minimal Fe-isotope fractionation in the dissolved Fe pool, we suggest that the pristine iron isotope composition of fresh water is preserved during estuarine mixing and that the value of the global riverine source into the ocean can be identified from the fresh water values. However, this study also suggests that δ⁵⁶Fe composition of rivers can also be characterized by more positive δ⁵⁶Fe values (up to 0.3‰) relative to the crust than previously reported. In order to improve our current understanding of the oceanic iron isotope cycling, further work is now required to determine the processes controlling the fractionation of Fe-isotopes during continental run-off.     

Controls on iron-isotope fractionation in organic-rich sediments (Kimmeridge Clay, Upper Jurassic, Southern England)

This study explores the fractionation of iron isotopes (⁵⁷Fe/⁵⁴Fe) in an organic-rich mudstone succession, focusing on core and outcrop material sampled from the Upper Jurassic Kimmeridge Clay Formation type locality in south Dorset, UK. The organic-rich environments recorded by the succession provide an excellent setting for an investigation of the mechanisms by which iron isotopes are partitioned among mineral phases during biogeochemical sedimentary processes.Two main types of iron-bearing assemblage are defined in the core material: mudstones with calcite ± pyrite ± siderite mineralogy, and ferroan dolomite (dolostone) bands. A cyclic data distribution is apparent, which reflects variations in isotopic composition from a lower range of δ⁵⁷Fe values associated with the pyrite/siderite mudstone samples to the generally higher values of the adjacent dolostone samples. Most pyrite/siderite mudstones vary between −0.4 and 0.1‰ while dolostones range between −0.1 and 0.5‰, although in very organic-rich shale samples below 360 m core depth higher δ⁵⁷Fe values are noted. Pyrite nodules and pyritized ammonites from the type exposure yield δ⁵⁷Fe values of −0.3 to −0.45‰. A fractionation model consistent with the δ⁵⁷Fe variations relates the lower δ⁵⁷Fe pyrite and siderite ± pyrite mudstones values to the production of isotopically depleted Fe(II) during biogenic reduction of the isotopically heavier lithogenic Fe(III) oxides. A consequence of this reductive dissolution is that a ⁵⁷Fe-enriched iron species must be produced that potentially becomes available for the formation of the higher δ⁵⁷Fe dolostones. An isotopic profile across a dolostone band reveals distinct zonal variations in δ⁵⁷Fe, characterized by two peaks, respectively located above and below the central part of the band, and decoupling of the isotopic composition from the iron content. This form of isotopic zoning is shown to be consistent with a one-dimensional model of diffusional-chromatographic Fe-isotope exchange between dolomite and isotopically enriched pore water. An alternative mechanism envisages the infiltration of dissolved ferrous iron from variable (high and low) δ⁵⁷Fe sources during coprecipitation of Fe(II) ion with dolomite. The study provides clear evidence that iron isotopes are cycled during the formation and diagenesis of organic carbon-rich sediments.     

Matrix-independent Fe isotope ratio determination in silicates using UV femtosecond laser ablation

UV femtosecond laser ablation coupled to MC-ICP-MS provides a promising in situ tool to investigate elemental and isotope ratios by non-matrix-matched calibration. In this study, we investigate Fe isotope composition in siliceous matrices including biotite, hornblende, garnet, fayalite and forsterite (San Carlos Olivine), and an oceanic Fe–Mn crust using the iron reference material IRMM-014 for calibration. To test the accuracy of the laser ablation data, Fe isotope compositions were obtained independently by solution ICP-MS after chromatographic separation of Fe. Sample materials with low Cr content, i.e. biotite, hornblende, fayalite and the Fe–Mn crust, reveal δ^56/54Fe and δ^57/54Fe values that agree with those from solution ICP-MS data within the measured precision. For high Cr concentration (⁵⁴Cr/⁵⁴Fe >0.0001), i.e. in the garnet and forsterite sample, δ^56/54Fe and δ^57/54Fe values were derived from ⁵⁷Fe/⁵⁶Fe ratios as correction of the isobaric interference of ⁵⁴Cr on ⁵⁴Fe is unsatisfactory. This approach provides accurate results for both minerals. Moreover, the garnet crystal exhibits isotopic zonation with differences of 0.3‰ in δ^56/54Fe showing that substantial heterogeneities exist in high-temperature metamorphic minerals. Multiple analyses of homogeneous sample materials reveal a repeatability of 0.1‰ (2 SD) for δ^56/54Fe and 0.2‰ (2 SD) for δ^57/54Fe, respectively. This study adds to the observations of Horn et al. (2006) who have shown that the determination of Fe isotope ratios in various matrices including iron alloys, iron oxides and hydroxides, iron sulfide and iron carbonates can be performed with high accuracy and precision at high spatial resolution using UV femtosecond laser ablation ICP-MS. These results demonstrate that femtosecond laser ablation ICP-MS is a largely matrix-independent method, which provides a substantial advantage over commonly employed nanosecond laser ablation systems.     

Inter-mineral Fe isotope variations in mantle-derived rocks and implications for the Fe geochemical cycle

Iron isotope and major- and minor-element compositions of coexisting olivine, clinopyroxene, and orthopyroxene from eight spinel peridotite mantle xenoliths; olivine, magnetite, amphibole, and biotite from four andesitic volcanic rocks; and garnet and clinopyroxene from seven garnet peridotite and eclogites have been measured to evaluate if inter-mineral Fe isotope fractionation occurs in high-temperature igneous and metamorphic minerals and if isotopic fractionation is related to equilibrium Fe isotope partitioning or a result of open-system behavior. There is no measurable fractionation between silicate minerals and magnetite in andesitic volcanic rocks, nor between olivine and orthopyroxene in spinel peridotite mantle xenoliths. There are some inter-mineral differences (up to 0.2‰ in ⁵⁶Fe/⁵⁴Fe) in the Fe isotope composition of coexisting olivine and clinopyroxene in spinel peridotites. The Fe isotope fractionation observed between clinopyroxene and olivine appears to be a result of open-system behavior based on a positive correlation between the Δ⁵⁶Fe_{clinopyroxene-olivine} fractionation and the δ⁵⁶Fe value of clinopyroxene and olivine. There is also a significant difference in the isotopic compositions of garnet and clinopyroxene in garnet peridotites and eclogites, where the average Δ⁵⁶Fe_{clinopyroxene-garnet} fractionation is +0.32 ± 0.07‰ for six of the seven samples. The one sample that has a lower Δ⁵⁶Fe_{clinopyroxene-garnet} fractionation of 0.08‰ has a low Ca content in garnet, which may reflect some crystal chemical control on Fe isotope fractionation. The Fe isotope variability in mantle-derived minerals is interpreted to reflect subduction of isotopically variable oceanic crust, followed by transport through metasomatic fluids. Isotopic variability in the mantle might also occur during crystal fractionation of basaltic magmas within the mantle if garnet is a liquidus phase. The isotopic variations in the mantle are apparently homogenized during melting processes, producing homogenous Fe isotope compositions during crust formation.     

Subsurface processes at the lucky strike hydrothermal field, Mid-Atlantic ridge: evidence from sulfur, selenium, and iron isotopes

At Lucky Strike near the Azores Triple Junction, the seafloor setting of the hydrothermal field in a caldera system with abundant low-permeability layers of cemented breccia, provides a unique opportunity to study the influence of subsurface geological conditions on the hydrothermal fluid evolution. Coupled analyses of S isotopes performed in conjunction with Se and Fe isotopes have been applied for the first time to the study of seafloor hydrothermal systems. These data provide a tool for resolving the different abiotic and potential biotic near-surface hydrothermal reactions. The δ³⁴S (between 1.5‰ and 4.6‰) and Se values (between 213 and 1640 ppm) of chalcopyrite suggest a high temperature end-member hydrothermal fluid with a dual source of sulfur: sulfur that was leached from basaltic rocks, and sulfur derived from the reduction of seawater sulfate. In contrast, pyrite and marcasite generally have lower δ³⁴S within the range of magmatic values (0 ± 1‰) and are characterized by low concentrations of Se (<50 ppm). For ⁸²Se/⁷⁶Se ratios, the δ⁸²Se values range from basaltic values of near −1.5‰ to −7‰. The large range and highly negative values of hydrothermal deposits observed cannot be explained by simple mixing between Se leached from igneous rock and Se derived from seawater. We interpret the Se isotope signature to be a result of leaching and mixing of a fractionated Se source located beneath hydrothermal chimneys in the hydrothermal fluid. At Lucky Strike we consider two sources for S and Se: (1) the “end-member” hydrothermal fluid with basaltic Se isotopic values (−1.5‰) and typical S isotope hydrothermal values of 1.5‰; (2) a fractionated source hosted in subsurface environment with negative δ³⁴S values, probably from bacterial reduction of seawater sulfate and negative δ⁸²Se values possibly derived from inorganic reduction of Se oxyanions. Fluid trapped in the subsurface environment is conductively cooled and has restricted mixing and provide favorable conditions for subsurface microbial activity which is potentially recorded by S isotopes. Fe isotope systematic reveals that Se-rich high temperature samples have δ⁵⁷Fe values close to basaltic values (∼0‰) whereas Se-depleted samples precipitated at medium to low temperature are systematically lighter (δ⁵⁷Fe values between −1 to −3‰). An important implication of our finding is that light Fe isotope composition down to −3.2‰ may be explained entirely by abiotic fractionation, in which a reservoir effect during sulfide precipitation was able to produce highly fractionated compositions.