Redox state of iron in peralkaline rhyolitic glass/melt: X-ray absorption micro-spectroscopy experiments at high temperature

Assessing the ferric-ferrous ratio in magmas prior to eruption remains a challenging task. X-ray absorption near-edge structure (μXANES) spectra were collected at the iron K-edge in water-rich peralkaline silicic melt/glass inclusions trapped in quartz. These experiments were carried out between 800 and 20 °C. The chemical environment of iron was also determined in the naturally quenched samples (glass inclusions and matrix glass) and in the peralkaline rhyolitic reference glasses, with variable [Fe³⁺ / ∑Fe] ratios.In the reference glasses, both the intensity of the pre-peaks (Fe²⁺, Fe³⁺) and site geometry of iron change as the oxidation state increases. Fourfold-coordinated Fe³⁺ prevails in highly oxidised peralkaline silicic glasses, using alkalis for charge balance. The position of the pre-edge centroid of the 1s-3d transition correlates with the Fe³⁺ / ΣFe ratios that allowed calibration of the redox state of iron of our natural samples.At high temperatures, Fe²⁺ dominates in the pre-edge structure of melt inclusions. Upon cooling down to 20 °C, the intensity of the Fe³⁺ peak increases, the centroid position of the pre-edge features shifts by nearly 0.5 eV and the main edge moves slightly towards higher energies. The slower the cooling rate, the higher the ferric iron contribution. Iterative μXANES experiments performed on the same samples show that the process is reversible. However, this apparent oxidation of iron upon cooling is an artefact of changes in Fe coordination. It implies that the [Fe³⁺ / ΣFe] ratio of glassy samples, measured at 20 °C, may be overestimated by a factor > 1.7, and that this ratio cannot be reliably retrieved by probing naturally cooled glass inclusions, and most silicate glasses. High temperature μXANES experiments led first to an assessment of the ferric-ferrous ratio in the water-rich peralkaline melt in pre-eruptive magmatic conditions and second to the determination of the corresponding oxygen fugacity at 740 °C.     

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

Iron isotope fractionation during skarn-type alteration: Implications for metal source in the Han-Xing iron skarn deposit

The Han-Xing iron mineralization in the central North China Craton is a typical Fe skarn deposit associated with altered diorites. Here we report the Fe isotopic compositions of whole rocks and mineral separates from this deposit with a view to evaluate the Fe isotope fractionation during the formation of Fe skarn deposit, and to constrain the metal source. The Fe isotopes show a large variation both in whole rocks and mineral separates. Altered diorites show a wide range in δ⁵⁶Fe values (− 0.07‰ to + 0.21‰ relative to the Fe isotope standard IRMM-014) which positively correlate with their TFe₂O₃/TiO₂ ratios (Fe₂O₃ and FeO calculated as TFe₂O₃). The positive correlation indicates that heavy Fe isotopes were preferentially leached from diorites during the skarn-type alteration. Among the metallic minerals, pyrite and pyrrhotite are isotopically heavier (+ 0.12‰ to + 0.48‰) than the magnetite (+ 0.07‰ to + 0.21‰). Fe isotope fractionation between mineral pairs demonstrates that magnetite did not attain Fe isotopic equilibrium with pyrite and pyrrhotite, whereas pyrite and pyrrhotite might have attained isotopic equilibrium. Petrological observations and major element data also suggest that iron was leached from the diorites during the skarn-type alteration. If the leached iron provides the main Fe budget of the Han-Xing Fe skarn deposit, magnetite in ores would be isotopically heavier than the unaltered diorite. However, our results are in contrast with the magnetite being isotopically lighter than the unaltered diorite. This suggests that the major Fe source of the Han-Xing Fe skarn deposit is not from the leaching of diorites, and might be from magmatic fluid which is isotopically lighter than the silicate melt. Our data demonstrate that Fe isotopes can be used as important tracers in deciphering the metal source of Fe skarn deposits.     

High-temperature Si isotope fractionation between iron metal and silicate

The Si stable isotope fractionation between metal and silicate has been investigated experimentally at 1800, 2000, and 2200 °C. We find that there is a significant silicon stable isotope fractionation at high temperature between metal and silicate in agreement with Shahar et al. (2009). Further we find that this fractionation is insensitive to the structure and composition of the silicate as the fractionation between silicate melt and olivine is insignificant within the error of the analyses. The temperature-dependent silicon isotope fractionation is Δ³⁰Si_{silicate-metal} = 7.45 ± 0.41 × 10⁶/T². We also demonstrate the viability of using laser ablation MC-ICPMS as a tool for measuring silicon isotope ratios in high pressure and temperature experiments.     

Kinetic isotope effect during reduction of iron from a silicate melt

Iron isotopic compositions measured in chondrules from various chondrites vary between δ⁵⁷Fe/⁵⁴Fe = +0.9‰ and −2.0‰, a larger range than for igneous rocks. Whether these compositions were inherited from chondrule precursors, resulted from the chondrule-forming process itself or were produced by later parent body alteration is as yet unclear. Since iron metal is a common phase in some chondrules, it is important to explore a possible link between the metal formation process and the observed iron isotope mass fractionation. In this experimental study we have heated a fayalite-rich composition under reducing conditions for heating times ranging from 2 min to 6 h. We performed chemical and iron isotope analyses of the product phases, iron metal and silicate glass. We demonstrated a lack of evaporation of Fe from the silicate melt in similar isothermal experiments performed under non-reducing conditions. Therefore, the measured isotopic mass fractionation in the glass, ranging between −0.32‰ and +3.0‰, is attributed to the reduction process. It is explained by the faster transport of lighter iron isotopes to the surface where reduction occurs, and is analogous to kinetic isotope fractionation observed in diffusion couples [Richter, F.M., Davis, A.M., Depaolo, D.J., Watson, E.B., 2003. Isotope fractionation by chemical diffusion between molten basalt and rhyolite. Geochim. Cosmochim. Acta67, 3905-3923]. The metal phase contains 90-99.8% of the Fe in the system and lacks significant isotopic mass fractionation, with values remaining similar to that of the starting material throughout. The maximum iron isotope mass fractionation in the glass was achieved within 1 h and was followed by an isotopic exchange and re-equilibration with the metal phase (incomplete at ∼6 h). This study demonstrates that reduction of silicates at high temperatures can trigger iron isotopic fractionation comparable in its bulk range to that observed in chondrules. Furthermore, if metal in Type I chondrules was formed by reduction of Fe silicate, our observed isotopic fractionations constrain chondrule formation times to approximately 60 min, consistent with previous work.     

Copper stable isotopes as tracers of metal-sulphide segregation and fractional crystallisation processes on iron meteorite parent bodies

We report high precision Cu isotope data coupled with Cu concentration measurements for metal, troilite and silicate fractions separated from magmatic and non-magmatic iron meteorites, analysed for Fe isotopes (δ⁵⁷Fe; permil deviation in ⁵⁷Fe/⁵⁴Fe relative to the pure iron standard IRMM-014) in an earlier study (Williams et al., 2006). The Cu isotope compositions (δ⁶⁵Cu; permil deviation in ⁶⁵Cu/⁶³Cu relative to the pure copper standard NIST 976) of both metals (δ⁶⁵Cu_M) and sulphides (δ⁶⁵Cu_FeS) span much wider ranges (−9.30 to 0.99‰ and −8.90 to 0.63‰, respectively) than reported previously. Metal-troilite fractionation factors (Δ⁶⁵Cu_M-FeS = δ⁶⁵Cu_M − δ⁶⁵Cu_FeS) are variable, ranging from −0.07 to 5.28‰, and cannot be explained by equilibrium stable isotope fractionation coupled with either mixing or reservoir effects, i.e. differences in the relative proportions of metal and sulphide in the meteorites. Strong negative correlations exist between troilite Cu and Fe (δ⁵⁷Fe_FeS) isotope compositions and between metal-troilite Cu and Fe (Δ⁵⁷Fe_M-FeS) isotope fractionation factors, for both magmatic and non-magmatic irons, which suggests that similar processes control isotopic variations in both systems. Clear linear arrays between δ⁶⁵Cu_FeS and δ⁵⁷Fe_FeS and calculated Cu metal-sulphide partition coefficients (D_Cu = [Cu]_metal/[Cu]_FeS) are also present. A strong negative correlation exists between Δ⁵⁷Fe_M-FeS and D_Cu; a more diffuse positive array is defined by Δ⁶⁵Cu_M-FeS and D_Cu. The value of D_Cu can be used to approximate the degree of Cu concentration equilibrium as experimental studies constrain the range of D_Cu between Fe metal and FeS at equilibrium to be in the range of 0.05-0.2; D_Cu values for the magmatic and non-magmatic irons studied here range from 0.34 to 1.11 and from 0.04 to 0.87, respectively. The irons with low D_Cu values (closer to Cu concentration equilibrium) display the largest Δ⁵⁷Fe_M-FeS and the lowest Δ⁶⁵Cu_M-FeS values, whereas the converse is observed in the irons with large values D_Cu that deviate most from Cu concentration equilibrium. The magnitudes of Cu and Fe isotope fractionation between metal and FeS in the most equilibrated samples are similar: 0.25 and 0.32‰/amu, respectively. As proposed in an earlier study (Williams et al., 2006) the range in Δ⁵⁷Fe_M-FeS values can be explained by incomplete Fe isotope equilibrium between metal and sulphide during cooling, where the most rapidly-cooled samples are furthest from isotopic equilibrium and display the smallest Δ⁵⁷Fe_M-FeS and largest D_Cu values. The range in Δ⁶⁵Cu_M-FeS, however, reflects the combined effects of partial isotopic equilibrium overprinting an initial kinetic signature produced by the diffusion of Cu from metal into exsolving sulphides and the faster diffusion of the lighter isotope. In this scenario, newly-exsolved sulphides initially have low Cu contents (i.e. high D_Cu) and extremely light δ⁶⁵Cu_FeS values; with progressive equilibrium and fractional crystallisation the Cu contents of the sulphides increase as their isotopic composition becomes less extreme and closer to the metal value. The correlation between Δ⁶⁵Cu_M-FeS and Δ⁵⁷Fe_M-FeS is therefore a product of the superimposed effects of kinetic fractionation of Cu and incomplete equilibrium between metal and sulphide for both isotope systems during cooling. The correlations between Δ⁶⁵Cu_M-FeS and Δ⁵⁷Fe_M-FeS are defined by both magmatic and non-magmatic irons record fractional crystallisation and cooling of metallic melts on their respective parent bodies as sulphur and chalcophile elements become excluded from crystallised solid iron and concentrated in the residual melt. Fractional crystallisation processes at shallow levels have been implicated in the two main classes of models for the origin of the non-magmatic iron meteorites; at (i) shallow levels in impact melt models and (ii) at much deeper levels in models where the non-magmatic irons represent metallic melts that crystallised within the interior of a disrupted and re-aggregated parent body. The presence of non-magmatic irons with a range of Fe and Cu isotope compositions, some of which record near-complete isotopic equilibrium implies crystallisation at a range of cooling rates and depths, which is most consistent with cooling within the interior of a meteorite parent body. Our data therefore lend support to models where the non-magmatic irons are metallic melts that crystallised in the interior of re-aggregated, partially differentiated parent bodies.     

Iron and lithium isotope systematics of the Hekla volcano,Iceland — Evidence for Fe isotope fractionation during magma differentiation

In this study potential iron isotope fractionation by magmatic processes in the Earth's crust was systematically investigated. High precision iron isotope analyses by MC-ICP-MS were performed on a suite of rock samples representative for the volcanic evolution of the Hekla volcano, Iceland. The whole series of Hekla's rocks results from several processes. (i) Basaltic magmas rise and induce partial melting of meta-basalts in the lower part of the Icelandic crust. The resulting dacitic magma evolves to rhyolitic composition through crystal fractionation. During this differentiation the δ^56/54Fe_IRMM-014 values increase successively from 0.051 ± 0.021‰ for the primitive dacites to 0.168 ± 0.021‰ for the rhyolites. This increase can be described by a Rayleigh fractionation model using a constant bulk fractionation factor between all mineral phases (M) and the silicate liquid (L) of Δ^56/54Fe_M–L = ? 0.1‰. (ii) The basaltic magma itself differentiates by crystal fractionation to basaltic andesite composition. No Fe isotope fractionation was found in this series. All basalts and basaltic andesites have an average δ^56/54Fe_IRMM-014 value of 0.062 ± 0.042‰ (2SD, n = 9), identical to mean terrestrial basaltic values reported in previous studies. This observation is consistent with the limited removal of iron from the remaining silicate melt through crystal fractionation and small mineral-melt Fe isotope fractionation factors expected at temperatures in excess of 1050 °C. (iii) Andesites are produced by mixing of basaltic andesite with dacitic melts. The iron isotope composition of the andesites is matching that of the basaltic andesites and the less evolved dacites, in agreement with a mixing process. In the Hekla volcanic suite Li concentrations are positively correlated with indicators of magma differentiation. All Hekla rocks have δ⁷Li values typical for the upper mantle and demonstrate the absence of resolvable Li isotope fractionation during crystal fractionation. As a fluid-mobile trace element, Li concentrations and isotopes are a potential tracer of magma/fluid interaction. At Hekla, Li concentrations and isotope compositions do not indicate any extensive fluid exsolution. Hence, the heavy Fe isotope composition of the dacites and rhyolites can be predominately attributed to fractional crystallisation. Iron isotope analyses on single samples from other Icelandic volcanoes (Torfajökull, Vestmannaeyjar) confirm heavy Fe isotope enrichment in evolving magmas. Our results suggest that the iron isotope composition of highly evolved crust can be slightly modified by magmatic processes.     

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.     

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.     

Redox equilibria of iron and silicate melt structure: Implications for olivine/melt element partitioning

Olivine/melt partitioning of ΣFe, Fe²⁺, Mg²⁺, Ca²⁺, Mn²⁺, Co²⁺, and Ni²⁺ has been determined in the systems CaO-MgO-FeO-Fe₂O₃-SiO₂ (FD) and CaO-MgO-FeO-Fe₂O₃-Al₂O₃-SiO₂ (FDA3) as a function of oxygen fugacity (f_O2) at 0.1 MPa pressure. Total iron oxide content of the starting materials was ∼20 wt%. The f_O2 was to used to control the Fe³⁺/ΣFe (ΣFe: total iron) of the melts. The Fe³⁺/ΣFe and structural roles of Fe²⁺ and Fe³⁺ were determined with ⁵⁷Fe resonant absorption Mössbauer spectroscopy. Changes in melt polymerization, NBO/T, as a function of f_O2 was estimated from the Mössbauer data and existing melt structure information. It varies by ∼100% in melts coexisting with olivine in the FDA3 system and by about 300% in the FD system in the Fe³⁺/ΣFe range of the experiments (0.805-0.092). The partition coefficients ( in olivine/wt% in melt) are systematic functions of f_O2 and, therefore, NBO/T of the melt. There is a -minimum in the FDA3 system at NBO/T-values corresponding to intermediate Fe³⁺/ΣFe (0.34-0.44). In the Al-free system, FD, where the NBO/T values of melts range between ∼1 and ∼2.9, the partition coefficients are positively correlated with NBO/T (decreasing Fe³⁺/ΣFe). These relationships are explained by consideration of solution behavior in the melts governed by Qⁿ-unit distribution and structural changes of the divalent cations in the melts (coordination number, complexing with Fe³⁺, and distortion of the polyhedra).     

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

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

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

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

Structure and Composition Effects on the Oxygen Isotope Fractionation in Silicate Melts

The influence of melt composition and structure on the oxygen isotope fractionation was studied for the multicomponent (SiO₂ ± TiO₂ + Al₂O₃ ± Fe₂O₃ + MgO ± CaO) system at 1500°C and 1 atm. The experiments show that significant oxygen isotope effects can be observed in silicate melts even at such high temperature. It is shown that the ability of silicate melt to concentrate ¹⁸O isotope is mainly determined by its structure. In particular, an increase of the NBO/T ratio in the experimental glasses from 0.11 to 1.34 is accompanied by a systematic change of oxygen isotope difference between melt and internal standard by values from–0.85 to +1.29‰. The obtained data are described by the model based on mass-balance equations and the inferred existence of O⁰, O^–, and O^2– (bridging, non-bridging, and free oxygen) ions in the melts. An application of the model requires the intra-structure isotope fractionation between bridging and non-bridging oxygens. Calculations show that the intra-structure isotope fractionation in our experiments is equal to 4.2 ± 1.0‰. To describe the obtained oxygen isotope effects at the melts relatively to temperature and fraction of non-bridging oxygen a general equation was proposed.     

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

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

The structural behavior of ferric and ferrous iron in aluminosilicate glass near meta-aluminosilicate joins

Iron-57 resonant absorption Mössbauer spectroscopy was used to describe the redox relations and structural roles of Fe³⁺ and Fe²⁺ in meta-aluminosilicate glasses. Melts were formed at 1500 °C in equilibrium with air and quenched to glass in liquid H₂O with quenching rates exceeding 200 °C/s. The aluminosilicate compositions were NaAlSi₂O₆, Ca_0.5AlSi₂O₆, and Mg_0.5AlSi₂O₆. Iron oxide was added in the form of Fe₂O₃, NaFeO₂, CaFe₂O₄, and MgFe₂O₄ with total iron oxide content in the range ∼0.9 to ∼5.6 mol% as Fe₂O₃. The Mössbauer spectra, which were deconvoluted by assuming Gaussian distributions of the hyperfine field, are consistent with one absorption doublet of Fe²⁺ and one of Fe³⁺. From the area ratios of the Fe²⁺ and Fe³⁺ absorption doublets, with corrections for differences in recoil-fractions of Fe³⁺ and Fe²⁺, the Fe³⁺/ΣFe is positively correlated with increasing total iron content and with decreasing ionization potential of the alkali and alkaline earth cation. There is a distribution of hyperfine parameters from the Mössbauer spectra of these glasses. The maximum in the isomer shift distribution function of Fe³⁺, δ_Fe³⁺, ranges from about 0.25 to 0.49 mm/s (at 298 K relative to Fe metal) with the quadrupole splitting maximum, Δ_Fe³⁺, ranging from ∼1.2 to ∼1.6 mm/s. Both δ_Fe³⁺ and δ_Fe²⁺ are negatively correlated with total iron oxide content and Fe³⁺/ΣFe. The dominant oxygen coordination number Fe³⁺ changes from 4 to 6 with decreasing Fe³⁺/ΣFe. The distortion of the Fe³⁺-O polyhedra of the quenched melts (glasses) decreases as the Fe³⁺/ΣFe increases. These polyhedra do, however, coexist with lesser proportions of polyhedra with different oxygen coordination numbers. The δ_Fe²⁺ and Δ_Fe²⁺ distribution maxima at 298 K range from ∼0.95 to 1.15 mm/s and 1.9 to 2.0 mm/s, respectively, and decrease with increasing Fe³⁺/ΣFe. We suggest that these hyperfine parameter values for the most part are more consistent with Fe²⁺ in a range of coordination states from 4- to 6-fold. The lower δ_Fe²⁺-values for the most oxidized melts are consistent with a larger proportion of Fe²⁺ in 4-fold coordination compared with more reduced glasses and melts.     

Solubility of palladium in picritic melts: 1. The effect of iron

In order to improve our understanding of HSE geochemistry, we evaluate the effect of Fe on the solubility of Pd in silicate melts. To date, experimentally determined Pd solubilities in silicate melt are only available for Fe-free anorthite-diopside eutectic compositions. Here we report experiments to study the solubility of Pd in a natural picritic melt as a function of pO₂ at 1300 °C in a one atm furnace. Palladium concentrations in the run products were determined by laser-ablation-ICP-MS. Palladium increases from 1.07 ± 0.26 ppm at FMQ-2, to 306 ± 19 ppm at FMQ+6.6. At a relative pO₂ of FMQ the slope in log Pd concentration vs. log pO₂ space increases considerably, and Pd concentrations are elevated over those established for AnDi melt compositions. In the same pO₂ range, ferric iron significantly increases relative to ferrous iron. Furthermore, at constant pO₂ (FMQ+0.5) Pd concentrations significantly increase with increasing X_FeO-total in the melt. Therefore, we consider ferric Fe to promote the formation of Pd²⁺ enhancing the solubility of Pd in the picrite melt significantly.The presence of FeO in the silicate melt has proven to be an important melt compositional parameter, and should be included and systematically investigated in future experimental studies, since most natural compositions have substantial FeO contents.     

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

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

Redox-controlled iron isotope fractionation during magmatic differentiation: an example from the Red Hill intrusion, S. Tasmania

This study presents accurate and precise iron isotopic data for 16 co-magmatic rocks and 6 pyroxene–magnetite pairs from the classic, tholeiitic Red Hill sill in southern Tasmania. The intrusion exhibits a vertical continuum of compositions created by in situ fractional crystallisation of a single injection of magma in a closed igneous system and, as such, constitutes a natural laboratory amenable to determining the causes of Fe isotope fractionation in magmatic rocks. Early fractionation of pyroxenes and plagioclase, under conditions closed to oxygen exchange, gives rise to an iron enrichment trend and an increase in $ f_{{{\text{O}}_{2} }} $ of the melt relative to the Fayalite–Magnetite–Quartz (FMQ) buffer. Enrichment in Fe³⁺/ΣFe_melt is mirrored by δ⁵⁷Fe, where ^VIFe²⁺-bearing pyroxenes partition ⁵⁷Fe-depleted iron, defining an equilibrium pyroxene-melt fractionation factor of $ \Updelta^{57} {\text{Fe}}_{{{\text{px}} - {\text{melt}}}} \le - 0.25\,\permille \times 10^{6} /T^{2} $ . Upon magnetite saturation, the $ f_{{{\text{O}}_{2} }} $ and δ⁵⁷Fe of the melt fall, commensurate with the sequestration of the oxidised, ⁵⁷Fe-enriched iron into magnetite, quantified as $ \Updelta^{57} {\text{Fe}}_{{{\text{mtn}} - {\text{melt}}}} = + 0.20\,\permille \times 10^{6} /T^{2} $ . Pyroxene–magnetite pairs reveal an equilibrium fractionation factor of $ \Updelta^{57} {\text{Fe}}_{{{\text{mtn}} - {\text{px}}}} \approx + 0.30\,\permille $ at 900–1,000?°C. Iron isotopes in differentiated magmas suggest that they may act as an indicator of their oxidation state and tectonic setting.