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.     

Combined iron and magnesium isotope geochemistry of pyroxenite xenoliths from Hannuoba,North China Craton: implications for mantle metasomatism

We present high-precision iron and magnesium isotopic data for diverse mantle pyroxenite xenoliths collected from Hannuoba, North China Craton and provide the first combined iron and magnesium isotopic study of such rocks. Compositionally, these xenoliths range from Cr-diopside pyroxenites and Al-augite pyroxenites to garnet-bearing pyroxenites and are taken as physical evidence for different episodes of melt injection. Our results show that both Cr-diopside pyroxenites and Al-augite pyroxenites of cumulate origin display narrow ranges in iron and magnesium isotopic compositions (δ⁵⁷Fe = ?0.01 to 0.09 with an average of 0.03 ± 0.08 (2SD, n = 6); δ²⁶Mg = ? 0.28 to ?0.25 with an average of ?0.26 ± 0.03 (2SD, n = 3), respectively). These values are identical to those in the normal upper mantle and show equilibrium inter-mineral iron and magnesium isotope fractionation between coexisting mantle minerals. In contrast, the garnet-bearing pyroxenites, which are products of reactions between peridotites and silicate melts from an ancient subducted oceanic slab, exhibit larger iron isotopic variations, with δ⁵⁷Fe ranging from 0.12 to 0.30. The δ⁵⁷Fe values of minerals in these garnet-bearing pyroxenites also vary widely (?0.25 to 0.08 in olivines, ?0.04 to 0.25 in orthopyroxenes, ?0.07 to 0.31 in clinopyroxenes, 0.07 to 0.48 in spinels and 0.31–0.42 in garnets). In addition, the garnet-bearing pyroxenite shows light δ²⁶Mg (?0.43) relative to the mantle. The δ²⁶Mg of minerals in the garnet-bearing pyroxenite range from ?0.35 for olivine and orthopyroxene, to ?0.34 for clinopyroxene, 0.04 for spinel and ?0.68 for garnet. These measured values stand in marked contrast to calculated equilibrium iron and magnesium isotope fractionation between coexisting mantle minerals at mantle temperatures derived from theory, indicating disequilibrium isotope fractionation. Notably, one phlogopite clinopyroxenite with an apparent later metasomatic overprint has the heaviest δ⁵⁷Fe (as high as 1.00) but the lightest δ²⁶Mg (as low as ?1.50) values of all investigated samples. Overall, there appears to be a negative co-variation between δ⁵⁷Fe and δ²⁶Mg in the Hannuoba garnet-bearing pyroxenite and in the phlogopite clinopyroxenite xenoliths and minerals therein. These features may reflect kinetic isotopic fractionation due to iron and magnesium inter-diffusion during melt–rock interaction. Such processes play an important role in producing inter-mineral iron and magnesium isotopic disequilibrium and local iron and magnesium isotopic heterogeneity in the subcontinental mantle.     

Iron isotope variations in spinel peridotite xenoliths from North China Craton: implications for mantle metasomatism

Iron isotopes, together with mineral elemental compositions of spinel peridotite xenoliths and clinopyroxenites from Hannuoba and Hebi Cenozoic alkaline basalts, were analyzed to investigate iron isotopic features of the lithospheric mantle beneath the North China Craton. The results show that the Hannuoba spinel peridotite xenoliths have small but distinguishable Fe isotopic variations. Overall variations in δ⁵⁷Fe are in a range of ?0.25 to 0.14‰ for olivine, ?0.17 to 0.17‰ for orthopyroxene, ?0.21 to 0.27‰ for clinopyroxene, and ?0.16 to 0.26‰ for spinel, respectively. Clinopyroxene has the heaviest iron isotopic ratio and olivine the lightest within individual sample. No clear linear relationships between the mineral pairs on “δ-δ” plot suggest that iron isotopes of mineral separates analyzed have been affected largely by some open system processes. The broadly negative correlations between mineral iron isotopes and metasomatic indexes such as spinel Cr#, (La/Yb)_N ratios of clinopyroxenes suggest that iron isotopic variations in different minerals and peridotites were probably produced by mantle metasomatism. The Hebi phlogopite-bearing lherzolite, which is significantly modified by metasomatic events, appears to be much heavier isotopically than clinopyroxene-poor lherzolite. This study further confirms previous conclusions that the lithospheric mantle has distinguishable and heterogeneous iron isotopic variations at the xenoliths scale. Mantle metasomatism is the most likely cause for the iron isotope variations in mantle peridotites.     

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.     

Erratum to: Iron isotope variations in spinel peridotite xenoliths from North China Craton: implications for mantle metasomatism

Iron isotopes, together with mineral elemental compositions of spinel peridotite xenoliths and clinopyroxenites from Hannuoba and Hebi Cenozoic alkaline basalts, were analyzed to investigate iron isotopic features of the lithospheric mantle beneath the North China Craton. The results show that the Hannuoba spinel peridotite xenoliths have small but distinguishable Fe isotopic variations. Overall variations in δ⁵⁷Fe are in a range of −0.25 to 0.14‰ for olivine, −0.17 to 0.17‰ for orthopyroxene, −0.21 to 0.27‰ for clinopyroxene, and −0.16 to 0.26‰ for spinel, respectively. Clinopyroxene has the heaviest iron isotopic ratio and olivine the lightest within individual sample. No clear linear relationships between the mineral pairs on “δ-δ” plot suggest that iron isotopes of mineral separates analyzed have been affected largely by some open system processes. The broadly negative correlations between mineral iron isotopes and metasomatic indexes such as spinel Cr#, (La/Yb)_N ratios of clinopyroxenes suggest that iron isotopic variations in different minerals and peridotites were probably produced by mantle metasomatism. The Hebi phlogopite-bearing lherzolite, which is significantly modified by metasomatic events, appears to be much heavier isotopically than clinopyroxene-poor lherzolite. This study further confirms previous conclusions that the lithospheric mantle has distinguishable and heterogeneous iron isotopic variations at the xenoliths scale. Mantle metasomatism is the most likely cause for the iron isotope variations in mantle peridotites.     

Li isotope fractionation in peridotites and mafic melts

We have measured the Li isotope ratios of a range of co-existing phases from peridotites and mafic magmas to investigate high-temperature fractionations of ⁷Li/⁶Li. The Li isotopic compositions of seven mantle peridotites, reconstructed from analyses of mineral separates, show little variation (δ⁷Li 3.2-4.9‰) despite a wide range in fertility and radiogenic isotopic compositions. The most fertile samples yield a best estimate of δ⁷Li ∼ 3.5‰ for the upper mantle. Bulk analyses of olivine separates from the xenoliths are typically ∼1.5‰ isotopically lighter than co-existing orthopyroxenes, suggestive of a small, high-temperature equilibrium isotope fractionation. On the other hand, bulk analyses of olivine phenocrysts and their host melts are isotopically indistinguishable. Given these observations, equilibrium mantle melting should generate melts with δ⁷Li little different from their sources (<0.5‰ lighter). In contrast to olivine and orthopyroxene, that dominate peridotite Li budgets, bulk clinopyroxene analyses are highly variable (δ⁷Li = 6.6‰ to −8.1‰). Phlogopite separated from a modally metasomatised xenolith yielded an extreme δ⁷Li of −18.9‰. Such large Li isotope variability is indicative of isotopic disequilibrium. This inference is strongly reinforced by in situ, secondary ion mass-spectrometry analyses which show Li isotope zonation in peridotite minerals. The simplest zoning patterns show isotopically light rims. This style of zoning is also observed in the phenocrysts of holocrystalline Hawaiian lavas. More dramatically, a single orthopyroxene crystal from a San Carlos xenolith shows a W-shaped Li isotope profile with a 40‰ range in δ⁷Li, close to the isotope variability seen in all terrestrial whole rock analyses. We attribute Li isotope zonation in mineral phases to diffusive fractionation of Li isotopes, within mineral phases and along melt pathways that pervade xenoliths. Given the high diffusivity of Li, the Li isotope profiles we observe can persist, at most, only a few years at magmatic temperatures. Our results thus highlight the potential of Li isotopes as a high-resolution geospeedometer of the final phases of magmatic activity and cooling.     

Lithium isotopic systematics of peridotite xenoliths from Hannuoba, North China Craton: Implications for melt-rock interaction in the considerably thinned lithospheric mantle

Li concentrations and isotopic compositions of coexisting minerals (ol, opx, and cpx) from peridotite xenoliths entrained in the Hannuoba Tertiary basalts, North China Craton, provide insight into Li isotopic fractionation between mantle minerals during melt-rock interaction in the considerably thinned lithospheric mantle. Bulk analyses of mineral separates show significant enrichment of Li in cpx (2.4-3.6 ppm) relative to olivine (1.2-1.8 ppm), indicating that these peridotites have been affected by mantle metasomatism with mafic silicate melts. Bulk olivine separates (δ⁷Li ∼ +3.3‰ to +6.4‰) are isotopically heavier than coexisting pyroxenes (δ⁷Li ∼ −3.3‰ to −8.2‰ in cpx, and −4.0‰ to −6.7‰ in opx). Such large variation suggests Li elemental and isotopic disequilibrium. This conclusion is supported by results from in situ SIMS analyses of mineral grains where significant Li elemental and isotopic zonations exist. The olivine and opx have lower Li concentrations and heavier Li isotopes in the rims than in the cores. This reverse correlation of δ⁷Li with Li concentrations indicates diffusive fractionation of Li isotopes. However, the zoning patterns in coexisting cpx show isotopically heavier rims with higher Li abundances. This positive correlation between δ⁷Li and Li concentrations suggests a melt mixing trend. We attribute Li concentration and isotope zonation in minerals to the effects of two-stage diffusive fractionation coupled with melt-rock interaction. The earliest melts may have been derived from the subducted oceanic slab with low δ⁷Li values produced by isotopic fractionation during the dehydration of the seawater-altered slab. Melts at later stages were derived from the asthenosphere and interacted with the peridotites, producing the Li elemental and isotopic zoning in mineral grains. These data thus provide evidence for multiple-stage peridotite-melt interaction in the lithospheric mantle beneath the northern North China Craton.     

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

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

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

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

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

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.     

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.     

Coupled lithium- and iron isotope fractionation during magmatic differentiation

In this study, we investigated Fe and Li isotope fractionation between mineral separates of olivine pheno- and xenocrysts (including one clinopyroxyene phenocryst) and their basaltic hosts. Samples were collected from the Canary Islands (Teneriffa, La Palma) and some German volcanic regions (Vogelsberg, Westerwald and Hegau). All investigated bulk samples fall in a tight range of Li and Fe isotope compositions (δ⁵⁶Fe_wr = 0.06–0.17‰ and δ⁷Li_ma = 2.5–5.2‰, assuming δ⁷Li of the olivine-free matrix is virtually identical to that of the bulk sample for mass balance reasons). In contrast, olivine phenocrysts display highly variable, but generally light Fe and mostly light Li isotope compositions compared to their respective olivine-free basaltic matrix, which was considered to represent the melt (with δ⁵⁶Fe_ol = ? 0.24 to 0.14‰ and δ⁷Li_ol = ? 10.5 to + 6.5‰, respectively). Single olivine crystals from one sample display even a larger range of δ⁵⁶Fe_ol between ? 0.7 and + 0.1‰. One single clinopyroxene phenocryst displays the lightest Li isotope composition (δ⁷Li_cpx = ? 17.7‰), but no Fe isotope fractionation relative to melt. The olivine phenocrysts show variable Mg# and Ni (correlated in most cases) that range between 0.89 and 0.74 and between 300 and 3000 μg/g, respectively. These olivines likely grew by fractional crystallization in an evolving magma. One sample from the Vogelsberg volcano contained olivine xenocrysts (Mg# > 0.89 and Ni > 3000 μg/g), in addition to olivine phenocrysts. This sample displays the highest Li- and the second highest Fe-isotope fractionation between olivine and melt (Δ⁷Li_ol-melt = ? 13; Δ⁵⁶Fe_ol-melt = ? 0.29).Our data, i.e. the variable olivine- at constant whole rock and matrix isotope compositions, strongly indicate disequilibrium, i.e. kinetic Fe and Li isotope fractionation between olivine and melt (for Li also between cpx and melt) during fractional crystallization. Δ⁷Li_ol-melt is correlated with the Li partitioning between olivine and melt (i.e. with Li_ol/Li_melt), indicating Li isotope fractionation due to preferential (faster) diffusion of ⁶Li into olivine during fractional crystallization. Olivine with low Δ⁷Li_ol-melt, also have low Δ⁵⁶Fe_ol-melt, indicating that Fe isotope fractionation is also driven by diffusion of isotopically light Fe into olivine, potentially, as Fe–Mg inter-diffusion. The lowest Δ⁵⁶Fe_ol-melt (? 0.40) was observed in a sample from Westerwald (Germany) with abundant magnetite, indicating relatively oxidizing conditions during magma differentiation. This may have enhanced equilibrium Fe isotope fractionation between olivine and melt or fine dispersed magnetite in the basalt matrix may have shifted its Fe isotope composition towards higher δ⁵⁶Fe. The decoupling of Li- and Fe isotope fractionation in cpx is likely due to faster diffusion of Li relative to Fe in cpx, implying that the large investigated cpx phenocryst resided in the magma for only a short period of time which was sufficient for Li- but not for Fe diffusion. The absence of any equilibrium Fe isotope fractionation between the investigated cpx phenocryst and its basaltic host may be related to the similar Fe^3 +/Fe^2 + of cpx and melt. In contrast to cpx, the generally light Fe isotope composition of all investigated olivine separates implies the existence of equilibrium- (in addition to diffusion-driven) isotope fractionation between olivine and melt, on the order of 0.1‰.     

Lithium isotopic disequilibrium of minerals in the spinel lherzolite xenoliths from Boeun,Korea

Spinel lherzolite xenoliths found in Boeun, Korea, have protogranular to porphyroclastic textures and are enclosed in a Miocene alkali basalt. The lithium concentration and isotopic compositions of olivine, clinopyroxene, and orthopyroxene separates from the spinel lherzolite, and whole rocks of the spinel lherzolites and alkali basalt were determined by inductively coupled plasma mass spectrometry (ICP-MS) and multi-collector inductively coupled plasma mass spectrometry (MC-ICP-MS). The lithium concentrations of the olivines and orthopyroxenes range from 2.2 to 5.0 ppm and from 2.1 to 6.4 ppm, respectively. In contrast, the clinopyroxenes have larger lithium concentrations, from 2.0 to 8.4 ppm, which reflect their preferential lithium enrichment. The lithium isotopic compositions (δ⁷Li) of olivines (-5.4 to + 3.5‰), orthopyroxenes (-11.4 to -0.1‰), and clinopyroxenes (-14.4 to -4.7‰) range far beyond the average mantle composition of + 4 ± 2‰. The lithium isotopic composition of the host rock, alkali basalt (3.4‰), is within the range of the intraplate and oceanic island basalts. The spinel lherzolites from Boeun exhibits strong elemental and isotopic disequilibria due to the different lithium and lithium isotope diffusion velocities in the olivine, orthopyroxene, and clinopyroxene minerals after eruption and magma cooling.     

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.     

Iron isotopic systematics of UHP eclogites respond to oxidizing fluid during exhumation

The iron stable isotope compositions (δ⁵⁶Fe) and iron valence states of ultrahigh‐pressure eclogites from Bixiling in the Dabie orogen belt, China, were measured to trace the changes of geochemical conditions during vertical transportation of earth materials, for example, oxygen fugacity. The bulk Fe³⁺/ΣFe ratios of retrograde eclogites, determined by Mössbauer spectroscopy, are consistently higher than those of fresh eclogites, suggesting oxidation during retrograde metamorphism and fluid infiltration. The studied eclogites (five samples) display limited mid‐ocean ridge basalts (MORB)‐like (~0.10‰) δ⁵⁶Fe values, which are indistinguishable from their protoliths, that is, gabbro cumulates formed through differentiation of mantle‐derived basaltic magma. This suggests that Fe isotope fractionation during continental subduction is limited. Garnet separates display limited δ⁵⁶Fe variation ranging from ?0.08 ± 0.07 ‰ to 0.02 ± 0.07‰, whereas coexisting omphacite displays a large variation of δ⁵⁶Fe values from 0.15 ± 0.07‰ to 0.47 ± 0.07‰. Omphacite also has highly variable Fe³⁺/ΣFe ratios from 0.367 ± 0.025 to 0.598 ± 0.024, indicating modification after peak metamorphism. Omphacite from retrograde eclogites has elevated Fe³⁺/ΣFe ratios (0.54–0.60) compared to that from fresh eclogites (~0.37), whereas garnet displays a narrow range of ferric iron content with Fe³⁺/ΣFe ratios from 0.039 ± 0.013 to 0.065 ± 0.022. The homogenous δ⁵⁶Fe values and Fe³⁺/ΣFe ratios of garnet suggest that it survived the retrograde metamorphism and preserved its Fe‐isotopic features and ferric contents of peak metamorphism. Because of similar diffusion rates of Fe and Mg in garnet and omphacite, and constant Δ²⁶Mg_{omphacite‐garnet} values (1.14 ± 0.04‰), equilibrium iron isotope fractionation between garnet and omphacite was probably achieved during peak metamorphism. Elevated Fe³⁺/ΣFe ratios of omphacite from retrograde eclogites and variant Δ⁵⁶Fe_{omphacite‐garnet} values of the studied eclogites (0.13 ± 0.10‰ to 0.48 ± 0.10‰) indicate that oxidized geofluid infiltration resulted in the elevation of δ⁵⁶Fe values of omphacite during retrograde metamorphism.     

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

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

Magmatic modification of the uppermost mantle beneath the Basin and Range to Colorado Plateau Transition Zone; Evidence from xenoliths, Wikieup, Arizona

Upper mantle xenoliths from Wikieup, AZ, provide abundant evidence for magmatic modification of the uppermost mantle beneath the Transition Zone between the Colorado Plateau and the southern Basin and Range province. Upper mantle lithologies in this xenolith suite are represented by spinel peridotite, wehrlite, plagioclase peridotite, and Al-augite group pyroxenites. Isotopic data for these xenoliths yield relatively uniform values and suggest a common petrogenesis. Al-augite-bearing gabbro and pyroxenite xenoliths from this locality are interpreted to have formed by crystal fractionation processes from parent alkali basalts similar to the Wikieup host basalt. Mineral and whole rock compositions show consistent trends of increasing incompatible element contents (Fe, Al, Ca, Na, K, LIL, and LREE), and decreasing compatible element contents (Mg, Cr, Ni) from spinel peridotite to wehrlite to plagioclase peridotite to the host basalt composition. These compositional trends are interpreted as resulting from varying degrees of magma-mantle wall rock interaction as ascending mafic magmas infiltrated upper mantle peridotite. Small degrees of melt infiltration resulted in slightly modified spinel peridotite compositions while moderate degrees metasomatized spinel peridotite to wehrlite, and the highest degrees metasomatized it to plagioclase peridotite. Whole rock compositions and clinopyroxene, plagioclase, and whole rock isotopic data suggest that the infiltrating magmas were the same as those from which the gabbros and pyroxenites crystallized, and that they were alkalic in composition, similar to the Wikieup host alkali olivine basalts. Relatively uniform ¹⁴³Nd/¹⁴⁴Nd for the mineral separates and whole rocks in spite of the significantly wide range in their ¹⁴⁷Sm/¹⁴⁴Nd (0.71–0.23 in clinopyroxene) suggests that the Wikieup xenoliths including gabbro, pyroxenite, peridotite, wehrlite, and plagioclase peridotite, are all relatively young rocks formed or metasomatized by a relatively recent magmatic episode. Received: 21 May 1996 / Accepted: 23 December 1996     

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

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

Fe isotope systematics of coexisting amphibole and pyroxene in the alkaline igneous rock suite of the Ilímaussaq Complex,South Greenland

The Ilímaussaq intrusion, South Greenland, provides an exceptional test case for investigating the changes of stable Fe isotope fractionation of solidus phases with changes in the Fe³⁺/∑Fe ratio of an evolving melt. The intrusion comprises a sequence of four melt batches that were fed from the same parental alkali basaltic magma. Differentiation produced cumulate rocks that range from augite syenite (phase I) over peralkaline granite (phase II) to agpaitic syenites (phases IIIa and IIIb). Fe³⁺/∑Fe ratios in amphiboles increase substantially from phase I to phase II and III rocks and mark a major change in the parental magma composition from augite syenites to peralkaline granites and agpaitic syenites. Before this transition, olivine, clinopyroxene, and amphibole in augite syenite, the most primitive rock type in the Ilímaussaq Complex, have a uniform Fe isotope composition that is identical to that of the bulk of igneous crustal rocks and approximated by the average isotopic composition of basalts (δ^56/54Fe_IRMM-014 = 0.072 ± 0.046‰). After the transition, amphiboles in the peralkaline granites and agpaitic syenites yield significantly heavier Fe isotope compositions with δ^56/54Fe_IRMM-014 values ranging from 0.123 to 0.237‰. Contamination of the Ilímaussaq magma by ongoing crustal assimilation as cause for this increase can be excluded on the grounds of Nd isotope data. Large-scale metasomatic overprint with an external fluid can also be dismissed based on amphibole O and Li isotope systematics. Rather, the increase towards heavy Fe isotope compositions most likely reflects the change in chemical compositions of amphiboles (calcic in augite syenite to sodic in the agpaitic syenites) and their Fe³⁺/ΣFe ratios that mirror changes in the chemical composition of the melt and its oxygen fugacity. A sensitive adjustment of equilibrium Fe isotope fractionation factors to amphibole ferric/ferrous ratios is also supported by beta-factors calculated from Mössbauer spetroscopy data. Comparison of the measured isotope fractionation between clinopyroxene and amphibole with that predicted from Mössbauer data reveal Fe isotope systematics close to equilibrium in augite syenites but Fe isotopic disequilibrium between these two phases in phase IIIa agpaitic syenites. These results are in agreement with O and Li isotope systematics. While amphiboles in all Ilímaussaq lithologies crystallized at temperatures between 650 and 850 °C, textural evidence reveals later clinopyroxene crystallization at temperatures as low as 300–400 °C. Therefore, isotopic equilibrium at crystallization conditions between these two phases can not be expected, but importantly, subsolidus reequilibration can also be dismissed.