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.     

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 variations in Holocene sediments of the Gotland Deep, Baltic Sea

Holocene sediments from the Gotland Deep basin in the Baltic Sea were investigated for their Fe isotopic composition in order to assess the impact of changes in redox conditions and a transition from freshwater to brackish water on the isotope signature of iron. The sediments display variations in δ⁵⁶Fe (differences in the ⁵⁶Fe/⁵⁴Fe ratio relative to the IRMM-14 standard) from −0.27 ± 0.09‰ to +0.21 ± 0.08‰. Samples deposited in a mainly limnic environment with oxygenated bottom water have a mean δ⁵⁶Fe of +0.08 ± 0.13‰, which is identical to the mean Fe isotopic composition of igneous rocks and oxic marine sediments. In contrast, sediments that formed in brackish water under periodically euxinic conditions display significantly lighter Fe isotope signatures with a mean δ⁵⁶Fe of −0.14 ± 0.19‰. Negative correlations of the δ⁵⁶Fe values with the Fe/Al ratio and S content of the samples suggest that the isotopically light Fe in the periodically euxinic samples is associated with reactive Fe enrichments and sulfides. This is supported by analyses of pyrite separates from this unit that have a mean Fe isotopic composition of −1.06 ± 0.20‰ for δ⁵⁶Fe. The supply of additional Fe with a light Fe isotopic signature can be explained with the shelf to basin Fe shuttle model. According to the Fe shuttle model, oxides and benthic ferrous Fe that is derived from dissimilatory iron reduction from shelves is transported and accumulated in euxinic basins. The data furthermore suggest that the euxinic water has a negative dissolved δ⁵⁶Fe value of about −1.4‰ to −0.9‰. If negative Fe isotopic signatures are characteristic for euxinic sediment formation, widespread euxinia in the past might have shifted the Fe isotopic composition of dissolved Fe in the ocean towards more positive δ⁵⁶Fe values.     

Iron isotope fractionation in subterranean estuaries

Dissolved Fe concentrations in subterranean estuaries, like their river-seawater counterparts, are strongly controlled by non-conservative behavior during mixing of groundwater and seawater in coastal aquifers. Previous studies at a subterranean estuary of Waquoit Bay on Cape Cod, USA demonstrate extensive precipitation of groundwater-borne dissolved ferrous iron and subsequent accumulation of iron oxides onto subsurface sands. Waquoit Bay is thus an excellent natural laboratory to assess the mechanisms of Fe-isotope fractionation in redox-stratified environments and determine potential Fe-isotope signatures of groundwater sources to coastal seawater. Here, we report Fe isotope compositions of iron-coated sands and porewaters beneath the intertidal zone of Waquoit Bay. The distribution of pore water Fe shows two distinct sources of Fe: one residing in the upward rising plume of Fe-rich groundwater and the second in the salt-wedge zone of pore water. The groundwater source has high Fe(II) concentration consistent with anoxic conditions and yield δ⁵⁶Fe values between 0.3 and −1.3‰. In contrast, sediment porewaters occurring in the mixing zone of the subterranean estuary have very low δ⁵⁶Fe values down to −5‰. These low δ⁵⁶Fe values reflect Fe-redox cycling and result from the preferential retention of heavy Fe-isotopes onto newly formed Fe-oxyhydroxides. Analysis of Fe-oxides precipitated onto subsurface sands in two cores from the subterranean estuary revealed strong δ⁵⁶Fe and Fe concentration gradients over less than 2m, yielding an overall range of δ⁵⁶Fe values between −2 and 1.5‰. The relationship between Fe concentration and δ⁵⁶Fe of Fe-rich sands can be modeled by the progressive precipitation of Fe-oxides along fluid flow through the subterranean estuary. These results demonstrate that large-scale Fe isotope fractionation (up to 5‰) can occur in subterranean estuaries, which could lead to coastal seawater characterized by very low δ⁵⁶Fe values relative to river values.     

The origin and history of ordinary chondrites: A study by iron isotope measurements of metal grains from ordinary chondrites

Chondrules and chondrites provide unique insights into early solar system origin and history, and iron plays a critical role in defining the properties of these objects. In order to understand the processes that formed chondrules and chondrites, and introduced isotopic fractionation of iron isotopes, we measured stable iron isotope ratios ⁵⁶Fe/⁵⁴Fe and ⁵⁷Fe/⁵⁴Fe in metal grains separated from 18 ordinary chondrites, of classes H, L and LL, ranging from petrographic types 3-6 using multi-collector inductively coupled plasma mass spectrometry. The δ⁵⁶Fe values range from −0.06 ± 0.01 to +0.30 ± 0.04‰ and δ⁵⁷Fe values are −0.09 ± 0.02 to +0.55 ± 0.05‰ (relative to IRMM-014 iron isotope standard). Where comparisons are possible, these data are in good agreement with published data. We found no systematic difference between falls and finds, suggesting that terrestrial weathering effects are not important in controlling the isotopic fractionations in our samples. We did find a trend in the ⁵⁶Fe/⁵⁴Fe and ⁵⁷Fe/⁵⁴Fe isotopic ratios along the series H, L and LL, with LL being isotopically heavier than H chondrites by ∼0.3‰ suggesting that redox processes are fractionating the isotopes. The ⁵⁶Fe/⁵⁴Fe and ⁵⁷Fe/⁵⁴Fe ratios also increase with increasing petrologic type, which again could reflect redox changes during metamorphism and also a temperature dependant fractionation as meteorites cooled. Metal separated from chondrites is isotopically heavier by ∼0.31‰ in δ⁵⁶Fe than chondrules from the same class, while bulk and matrix samples plot between chondrules and metal. Thus, as with so many chondrite properties, the bulk values appear to reflect the proportion of chondrules (more precisely the proportion of certain types of chondrule) to metal, whereas chondrule properties are largely determined by the redox conditions during chondrule formation. The chondrite assemblages we now observe were, therefore, formed as a closed system.     

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

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.     

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.     

The isotopic effects of electron transfer: An explanation for Fe isotope fractionation in nature

Isotope fractionation of electroplated Fe was measured as a function of applied electrochemical potential. As plating voltage was varied from −0.9 V to 2.0 V, the isotopic signature of the electroplated iron became depleted in heavy Fe, with δ⁵⁶Fe values (relative to IRMM-14) ranging from −0.18(±0.02) to −2.290(±0.006) ‰, and corresponding δ⁵⁷Fe values of −0.247(±0.014) and −3.354(±0.019) ‰. This study demonstrates that there is a voltage-dependent isotope fractionation associated with the reduction of iron. We show that Marcus’s theory for the kinetics of electron transfer can be extended to include the isotope effects of electron transfer, and that the extended theory accounts for the voltage dependence of Fe isotope fractionation. The magnitude of the electrochemically-induced fractionation is similar to that of Fe reduction by certain bacteria, suggesting that similar electrochemical processes may be responsible for biogeochemical Fe isotope effects. Charge transfer is a fundamental physicochemical process involving Fe as well as other transition metals with multiple isotopes. Partitioning of isotopes among elements with varying redox states holds promise as a tool in a wide range of the Earth and environmental sciences, biology, and industry.     

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.     

Iron isotope fractionation in oxic soils by mineral weathering and podzolization

Stable iron isotope ratios in three soils (two Podzols and one Cambisol) were measured by MC-ICPMS to investigate iron isotope fractionation during pedogenic iron transformation and translocation processes under oxic conditions. Podzolization is a soil forming process in which iron oxides are dissolved and iron is translocated and enriched in the subsoil under the influence of organic ligands. The Cambisol was studied for comparison, representing a soil formed by chemical weathering without significant translocation of iron. A three-step sequential extraction procedure was used to separate operationally-defined iron mineral pools (i.e., poorly-crystalline iron oxides, crystalline iron oxides, silicate-bound iron) from the soil samples. Iron isotope ratios of total soil digests were compared with those of the separated iron mineral pools. Mass balance calculations demonstrated excellent agreement between results of sequential extractions and total soil digestions. Systematic variations in the iron isotope signature were found in the Podzol profiles. An enrichment of light iron isotopes of about 0.6‰ in δ⁵⁷Fe was found in total soil digests of the illuvial Bh horizons which can be explained by preferential translocation of light iron isotopes. The separated iron mineral pools revealed a wide range of δ⁵⁷Fe values spanning more than 3‰ in the Podzol profiles. Strong enrichments of heavy iron isotopes in silicate-bound iron constituting the residue of weathering processes, indicated the preferential transformation of light iron isotopes during weathering. Iron isotope fractionation during podzolization is probably linked to the ligand-controlled iron translocation processes. Comparison of iron isotope data from eluvial and illuvial horizons of the Podzol profiles revealed that some iron must have been leached out of the profile. However, uncertainties in the initial iron content and iron isotopic composition of the parent materials prevented thorough mass balance calculations of iron fluxes within the profiles. In contrast to the Podzol profiles, the Cambisol profile displayed uniform δ⁵⁷Fe values across soil depth and showed only a small enrichment of light iron isotopes of about 0.4‰ in the poorly-crystalline iron oxide pool extracted by 0.5 M HCl. This work demonstrates that significant iron isotope fractionations can occur during pedogenesis in oxic environments under the influence of organic ligands. Our findings provide new insights into fractionation mechanisms of iron isotopes and will help in the development of stable iron isotopes as tracers for biogeochemical iron cycling in nature.     

The continental shelf benthic iron flux and its isotope composition

Benthic iron fluxes from sites along the Oregon-California continental shelf determined using in situ benthic chambers, range from less than 10 μmol m⁻² d⁻¹ to values in excess of ∼300 μmol m⁻² d⁻¹. These fluxes are generally greater than previously published iron fluxes for continental shelves contiguous with the open ocean (as opposed to marginal seas, bays, or estuaries) with the highest fluxes measured in the regions around the high-sediment discharge Eel River and the Umpqua River. These benthic iron fluxes do not covary with organic carbon oxidation rates in any systematic fashion, but rather seem to respond to variations in bottom water oxygen and benthic oxygen demand. We hypothesize that the highest rates of benthic iron efflux are driven, in part, by the greater availability of reactive iron deposited along these river systems as compared to other more typical continental margin settings. Bioirrigation likely plays an important role in the benthic Fe flux in these systems as well. However, the influence of bottom water oxygen concentrations on the iron flux is significant, and there appears to be a threshold in dissolved oxygen (∼60-80 μM), below which sediment-ocean iron exchange is enhanced. The isotope composition of this shelf-derived benthic iron is enriched in the lighter isotopes, and appears to change by ∼3‰ (δ⁵⁶Fe) during the course of a benthic chamber experiment with a mean isotope composition of −2.7 ± 1.1‰ (2 SD, n = 9) by the end of the experiment. This average value is slightly heavier than those from two high benthic Fe flux restricted basins from the California Borderland region where δ⁵⁶Fe is −3.4 ± 0.4‰ (2 SD, n = 3). These light iron isotope compositions support previous ideas, based on sediment porewater analyses, suggesting that sedimentary iron reduction fractionates iron isotopes and produces an isotopically light iron pool that is transferred to the ocean water column. In sum, our data suggest that continental shelves may export a higher efflux of iron than previously hypothesized, with the likelihood that along river-dominated margins, the benthic iron flux could well be orders of magnitude larger than non-river dominated shelves. The close proximity of the continental shelf benthos to the productive surface ocean means that this flux is likely to be essential for maintaining ecosystem micronutrient supply.     

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.     

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.     

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.     

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.     

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 Isotopic Compositions of Forearc Mantle Peridotites Constrained by Ophiolites from the North Qilian Orogen,Northern Tibet

Previous studies on iron isotope compositions of subduction zone magmas have revealed significant and complex variations that have great bearings on petrogenetic processes in the mantle wedge, e.g., partial melting, fluid metasomatism and redox state. However, interpretations for the fractionations are highly debatable and lack direct constraints from mantle wedge peridotites. This study presents iron isotope compositions for whole rocks and mineral separates in fresh forearc peridotites from the Yushigou ophiolite, North Qilian orogen in northern Tibet. Major and trace element compositions of whole rock and mineral indicate that the peridotites are highly depleted forearc peridotites overprinted by melt metasomatism, in contrast to the long‐holding opinion that the peridotites are derived from mid‐oceanic ridges. The minerals fall on a line with a slope of ～1 on the plot of δ⁵⁶Fe vs. δ⁵⁶Fe, indicating isotope equilibrium between minerals. δ⁵⁶Fe fractionation between olivine and orthopyroxene is within the range of 0～0.05, while fractionation between olivine and spinel is about 0.05～0.10. The fractionation trend between olivine and spinel is opposite to previous theoretical and experimental constraints, which may be due to substantial Cr substitution into the spinel. This indicates that negative correlations between spinel Cr#, fO₂ and spinel δ⁵⁶Fe in previous studies are probably a reflection of gradual Cr enrichment in spinel during melt extraction, and spinel δ⁵⁶Fe values are not a proxy for oxygen fugacity. Whole rock δ⁵⁶Fe values are well correlated with mineral δ⁵⁶Fe values, varying from overlapping with depleted mantle to slightly lower than depleted mantle. Therefore, variations in iron isotope compositions of subduction zone magmas are probably due to combined effect of source heterogeneity and partial melting fractionation.     

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.