Iron speciation and isotope fractionation during silicate weathering and soil formation in an alpine glacier forefield chronosequence

The chemical weathering of primary Fe-bearing minerals, such as biotite and chlorite, is a key step of soil formation and an important nutrient source for the establishment of plant and microbial life. The understanding of the relevant processes and the associated Fe isotope fractionation is therefore of major importance for the further development of stable Fe isotopes as a tracer of the biogeochemical Fe cycle in terrestrial environments. We investigated the Fe mineral transformations and associated Fe isotope fractionation in a soil chronosequence of the Swiss Alps covering 150 years of soil formation on granite. For this purpose, we combined for the first time stable Fe isotope analyses with synchrotron-based Fe-EXAFS spectroscopy, which allowed us to interpret changes in Fe isotopic composition of bulk soils, size fractions, and chemically separated Fe pools over time in terms of weathering processes. Bulk soils and rocks exhibited constant isotopic compositions along the chronosequence, whereas soil Fe pools in grain size fractions spanned a range of 0.4‰ in δ⁵⁶Fe. The clay fractions (<2 μm), in which newly formed Fe(III)-(hydr)oxides contributed up to 50% of the total Fe, were significantly enriched in light Fe isotopes, whereas the isotopic composition of silt and sand fractions, containing most of the soil Fe, remained in the range described by biotite/chlorite samples and bulk soils. Iron pools separated by a sequential extraction procedure covered a range of 0.8‰ in δ⁵⁶Fe. For all soils the lightest isotopic composition was observed in a 1 M NH₂OH-HCl-25% acetic acid extract, targeting poorly-crystalline Fe(III)-(hydr)oxides, compared with easily leachable Fe in primary phyllosilicates (0.5 M HCl extract) and Fe in residual silicates. The combination of the Fe isotope measurements with the speciation data obtained by Fe-EXAFS spectroscopy permitted to quantitatively relate the different isotope pools forming in the soils to the mineral weathering reactions which have taken place at the field site. A kinetic isotope effect during the Fe detachment from the phyllosilicates was identified as the dominant fractionation mechanism in young weathering environments, controlling not only the light isotope signature of secondary Fe(III)-(hydr)oxides but also significantly contributing to the isotope signature of plants. The present study further revealed that this kinetic fractionation effect can persist over considerable reaction advance during chemical weathering in field systems and is not only an initial transient phenomenon.     

Stable Cu isotope fractionation in soils during oxic weathering and podzolization

Copper stable isotope ratios are fractionated during various biogeochemical processes and may trace the fate of Cu during long-term pedogenetic processes. We assessed the effects of oxic weathering (formation of Cambisols) and podzolization on Cu isotope ratios (δ⁶⁵Cu). Two Cambisols (oxic weathered soils without strong vertical translocations of soil constituents) and two Podzols (soils showing vertical translocation of organic matter, Fe and Al) were analyzed for Cu concentrations, partitioning of Cu in seven fractions of a sequential extraction and δ⁶⁵Cu values in bulk soil. Cu concentrations in the studied soils were low (1.4-27.6 μg g⁻¹) and Cu was mainly associated with strongly bound Fe oxide- and silicate-associated forms. Bulk δ⁶⁵Cu values varied between −0.57‰ and 0.44‰ in all studied horizons. The O horizons had on average significantly lighter Cu isotope compositions (−0.21‰) than the A horizons (0.13‰) which can either be explained by Cu isotope fractionation during cycling through the plants or deposition of isotopically light Cu from the atmosphere. Oxic weathering without pronounced podzolization in both Cambisols and a weakly developed Podzol (Haplic Podzol 2) caused no significant isotope fractionation in the single profiles, while a slight tendency to lower δ⁶⁵Cu values with depth was visible in all four profiles. This is the opposite depth distribution of δ⁶⁵Cu values to that we observed in hydromorphic soils (soils which show indication of redox changes because of the influence of water saturation) in a previous study. In a more pronounced Podzol (Haplic Podzol 1), δ⁶⁵Cu values and Cu concentrations decreased from Ah to E horizons and increased again deeper in the soil. Humus-rich sections of the Bhs horizon had higher Cu concentrations (2.8 μg g⁻¹) and a higher δ⁶⁵Cu value (−0.18‰) than oxide-rich sections (1.9 μg g⁻¹, −0.35‰) suggesting Cu translocation between E and B horizons as organo-Cu complexes. The different depth distributions in oxic weathered and hydromorphic soils and the pronounced vertical differences in δ⁶⁵Cu values in Haplic Podzol 1 indicate a promising potential of δ⁶⁵Cu values to improve our knowledge of the fate of Cu during long-term pedogenetic processes.     

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

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

Lithium isotope systematics in a forested granitic catchment (Strengbach, Vosges Mountains, France)

Over the last decade it has become apparent that Li isotopes may be a good proxy to trace silicate weathering. However, the exact mechanisms which drive the behaviour of Li isotopes in surface environments are not totally understood and there is a need to better calibrate and characterize this proxy. In this study, we analysed the Li concentrations and isotopic compositions in the various surface reservoirs (soils, rocks, waters and plants) of a small forested granitic catchment located in the Vosges Mountains (Strengbach catchment, France, OHGE http://ohge.u-strasbg.fr). Li fluxes were calculated in both soil profiles and at the basin scale and it was found that even in this forested basin, atmospheric inputs and litter fall represented a minor flux compared to input derived from the weathering of rocks and soil minerals (which together represent a minimum of 70% of dissolved Li). Li isotope ratios in soil pore waters show large depth dependent variations. Average dissolved δ⁷Li decreases from −1.1‰ to −14.4‰ between 0 and −30 cm, but is +30.7‰ at −60 cm. This range of Li isotopic compositions is very large and it encompasses almost the entire range of terrestrial Li isotope compositions that have been previously reported. We interpret these variations to result from both the dissolution and precipitation of secondary phases. Large isotopic variations were also measured in the springs and stream waters, with δ⁷Li varying from +5.3‰ to +19.6‰. δ⁷Li increases from the top to the bottom of the basin and also covaries with discharge at the outlet. These variations are interpreted to reflect isotopic fractionations occurring during secondary phase precipitation along the water pathway through the rocks. We suggest that the dissolved δ⁷Li increases with increasing residence time of waters through the rocks, and so with increasing time of interaction between waters and solids. A dissolution precipitation model was used to fit the dissolved Li isotopic compositions. It was found that the isotopic compositions of springs and stream waters are explicable by an isotopic fractionation of −5‰ to −14‰ (best fit −10.8‰), in agreement with Li incorporation into clay. In soil solutions, it was found that isotopic fractionation during secondary precipitation is larger (at least −23‰), suggesting a major role for different secondary phases, such as iron oxides that maybe incorporate Li with a higher isotope fractionation.     

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

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

Impact of soil weathering degree on silicon isotopic fractionation during adsorption onto iron oxides in basaltic ash soils, Cameroon

The sequestration of silicon in soil clay-sized iron oxides may affect the terrestrial cycle of Si. Iron oxides indeed specifically adsorb aqueous monosilicic acid (H₄SiO₄⁰), thereby influencing Si concentration in soil solution. Here we study the impact of H₄SiO₄⁰ adsorption on the fractionation of Si isotopes in basaltic ash soils differing in weathering degree (from two weathering sequences, Cameroon), hence in clay and Fe-oxide contents, and evaluate the potential isotopic impact on dissolved Si in surrounding Cameroon rivers. Adsorption was measured in batch experiment series designed as function of time (0-72 h) and initial concentration (ic) of Si in solution (0.61-1.18 mM) at 20 °C, constant pH (5.5) and ionic strength (1 mM). After various soil-solution contact times, the δ³⁰Si vs. NBS28 compositions were determined in selected solutions by MC-ICP-MS (Nu Plasma) in medium resolution, operating in dry plasma with Mg doping with an average precision of ±0.15‰ (±2σ_SEM). The quantitative adsorption of H₄SiO₄⁰ by soil Fe-oxides left a solution depleted in light Si isotopes, which confirms previous study on synthetic Fe-oxides. Measured against its initial composition (δ³⁰Si = +0.02 ± 0.07‰ (±2σ_SD)), the solutions were systematically enriched in ³⁰Si reaching maximum δ³⁰Si values ranging between +0.16‰ and +0.95‰ after 72 h contact time. The enrichment of the solution in heavy isotopes increased with increasing values of three parameters: soil weathering degree, iron oxide content, and proportion of short-range ordered Fe-oxide. The Si-isotopic signature of the solution was partly influenced by Si release, possibly through mineral dissolution and Si desorption from oxide surfaces, depending on soil type, highlighting the complex pattern of natural soils. Surrounding Cameroon rivers displayed a mean Si-isotopic signature of +1.19‰. Our data imply that in natural environments, H₄SiO₄⁰ adsorption by soil clay-sized Fe-oxides at least partly impacts the Si-isotopic signature of the soil solution exported to water streams.     

Natural fractionation of U/U

The isotopic composition of U in nature is generally assumed to be invariant. Here, we report variations of the ²³⁸U/²³⁵U isotope ratio in natural samples (basalts, granites, seawater, corals, black shales, suboxic sediments, ferromanganese crusts/nodules and BIFs) of ∼1.3‰, exceeding by far the analytical precision of our method (≈0.06‰, 2SD). U isotopes were analyzed with MC-ICP-MS using a mixed ²³⁶U-²³³U isotopic tracer (double spike) to correct for isotope fractionation during sample purification and instrumental mass bias. The largest isotope variations found in our survey are between oxidized and reduced depositional environments, with seawater and suboxic sediments falling in between. Light U isotope compositions (relative to SRM-950a) were observed for manganese crusts from the Atlantic and Pacific oceans, which display δ²³⁸U of −0.54‰ to −0.62‰ and for three of four analyzed Banded Iron Formations, which have δ²³⁸U of −0.89‰, −0.72‰ and −0.70‰, respectively. High δ²³⁸U values are observed for black shales from the Black Sea (unit-I and unit-II) and three Kupferschiefer samples (Germany), which display δ²³⁸U of −0.06‰ to +0.43‰. Also, suboxic sediments have slightly elevated δ²³⁸U (−0.41‰ to −0.16‰) compared to seawater, which has δ²³⁸U of −0.41 ± 0.03‰. Granites define a range of δ²³⁸U between −0.20‰ and −0.46‰, but all analyzed basalts are identical within uncertainties and slightly lighter than seawater (δ²³⁸U = −0.29‰).Our findings imply that U isotope fractionation occurs in both oxic (manganese crusts) and suboxic to euxinic environments with opposite directions. In the first case, we hypothesize that this fractionation results from adsorption of U to ferromanganese oxides, as is the case for Mo and possibly Tl isotopes. In the second case, reduction of soluble U^VI to insoluble U^IV probably results in fractionation toward heavy U isotope compositions relative to seawater. These findings imply that variable ocean redox conditions through geological time should result in variations of the seawater U isotope compositions, which may be recorded in sediments or fossils. Thus, U isotopes might be a promising novel geochemical tracer for paleo-redox conditions and the redox evolution on Earth. The discovery that ²³⁸U/²³⁵U varies in nature also has implications for the precision and accuracy of U-Pb dating. The total observed range in U isotope compositions would produce variations in ²⁰⁷Pb/²⁰⁶Pb ages of young U-bearing minerals of up to 3 Ma, and up to 2 Ma for minerals that are 3 billion years old.     

Calcium isotopes in a proglacial weathering environment: Damma glacier, Switzerland

The biogeochemical cycling and isotopic fractionation of calcium during the initial stages of weathering were investigated in an alpine soil chronosequence (Damma glacier, Switzerland). This site has a homogeneous silicate lithology and minimal biological impacts due to sparse vegetation cover. Calcium isotopic compositions, obtained by TIMS using a ⁴³Ca-⁴⁶Ca double spike, were measured in the main Ca pools. During this very early stage of weathering, the young soils which have formed (δ^44/42Ca=+0.44‰) were indistinguishable to the rocks from which they were derived (δ^44/42Ca=+0.44‰) and stream water (δ^44/42Ca=+0.48‰) was also within error of the average rock. This lack of variation indicates that the dissolution of the bulk silicate rock does not strongly fractionate Ca isotopes. The only Ca pool which was strongly fractionated from bulk rock was vegetation, which exhibited an enrichment of light Ca isotopes. Significant Ca isotope fractionation between bulk rock and the dissolved flux of Ca is likely to only occur where the Ca biogeochemical cycle is dominated by secondary processes such as biological cycling, adsorption and secondary mineral precipitation.     

Natural variations of δSi ratios during progressive basalt weathering, Hawaiian Islands

Silicon stable isotopes can be used to trace the biogeochemical pathways of Si as it moves from its continental sources to its sink in ocean sediments. Along the way, Si is incorporated into clay minerals, taken up by plants where it forms plant opal, and leached into rivers, the major land-to-ocean conduit. Compared to igneous rocks, the waters that drain continents are enriched in heavy Si isotopes, but the mechanisms that control fractionation have not been elucidated. We studied Si isotope fractionation along a 4 million yr basaltic soil chronosequence on the Hawaiian Islands. Using the natural context of these samples in combination with laboratory experiments, we demonstrate that the isotopic composition of dissolved Si in weathering systems is determined by the combined effects of rock disintegration, clay mineral neosynthesis, and Si biocycling. Weathering preferentially releases ²⁸Si into solution, whereas secondary mineral formation preferentially removes ²⁸Si from solution. In humid environments, leached soils have lost large amounts of this soluble Si, thus creating a net loss of ³⁰Si from the entire soil system. As soils develop and greater fractions of Si reside in neoformed clay minerals, δ³⁰Si_{bulk soil} values change progressively toward more negative values; basalt δ³⁰Si values are about −0.5‰, but older soils have δ³⁰Si values up to −2.5‰. The difference between the solid and solution δ³⁰Si values remains more or less constant with progressive weathering, and therefore, soil water from older soils has a more negative δ³⁰Si composition. In the upper horizons of the Hawaiian soils, this weathering-driven δ³⁰Si shift is modified by the addition of unweathered primary minerals via dust, carrying δ³⁰Si values of about −0.5‰, and by biocycling of Si via plants, producing negative δ³⁰Si values in phytoliths and positive δ³⁰Si values in soil solutions derived from upper horizons. Due to the high concentrations of dissolved Si in these near-surface layers, rivers have more positive δ³⁰Si values than predicted based on the weathering status of the lower horizons. When combined with published δ³⁰Si values from large rivers worldwide, we find that the results from Hawaii point to weathering control of Si isotopes delivered to the oceans, and thus, to an important continent-ocean linkage that warrants further investigation.     

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

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

Adsorption as a cause for iron isotope fractionation in reduced groundwater

Iron isotopes were used to investigate iron transformation processes during an in situ field experiment for removal of dissolved Fe from reduced groundwater. This experiment provided a unique setting for exploring Fe isotope fractionation in a natural system. Oxygen-containing water was injected at a test well into an aquifer containing Fe(II)-rich reduced water, leading to oxidation of Fe(II) and precipitation of Fe(III)(hydr)oxides. Subsequently, groundwater was extracted from the same well over a time period much longer than the injection time. Since the surrounding water is rich in Fe(II), the Fe(II) concentration in the extracted water increased over time. The increase was strongly retarded in comparison to a conservative tracer added to the injected solution, indicating that adsorption of Fe(II) onto the newly formed Fe(III)(hydr)oxides occurred. A series of injection-extraction (push-pull) cycles were performed at the same well. The δ⁵⁷Fe/⁵⁴Fe of pre-experiment background groundwater (−0.57 ± 0.17 ‰) was lighter than the sediment leach of Fe(III) (−0.24 ± 0.08 ‰), probably due to slight fractionation (only ∼0.3 ‰) during microbial mediated reductive dissolution of Fe(III)(hydr)oxides present in the aquifer. During the experiment, Fe(II) was adsorbed from native groundwater drawn into the oxidized zone and onto Fe(III)(hydr)oxides producing a very light groundwater component with δ⁵⁷Fe/⁵⁴Fe as low as −4 ‰, indicating that heavier Fe(II) is preferentially adsorbed to the newly formed Fe(III)(hydr)oxides surfaces. Iron concentrations increased with time of extraction, and δ⁵⁷Fe/⁵⁴Fe linearly correlated with Fe concentrations (R² = 0.95). This pattern was reproducible over five individual cycles, indicating that the same process occurs during repeated injection/extraction cycles. We present a reactive transport model to explain the observed abiotic fractionation due to adsorption of Fe(II) on Fe(III)(hydr)oxides. The fractionation is probably caused by isotopic differences in the equilibrium sorption constants of the various isotopes (K_ads) and not by sorption kinetics. A fractionation factor α_57/54 of 1.001 fits the observed fractionation.     

Ca isotope cycling in a forested ecosystem

Reports of large Ca isotope fractionations between trees and soils prompted this study of a Boreal forest ecosystem near La Ronge, Saskatchewan, to improve understanding of this phenomenon. The results on five tree species (black spruce, trembling aspen, white spruce, jack pine, balsam poplar) confirm that nutrient Ca uptake by plants favors the light isotopes, thus driving residual Ca in plant available soil pools towards enrichment in the heavy isotopes. Substantial within-tree fraction occurs in tissues formed along the transpiration stream, with low δ⁴⁴Ca values in fine roots (2 mm), intermediate values in stemwood, and high values in foliage. Separation factors between different plant tissues are similar between species, but the initial fractionation step in the tips of the fine roots is species specific, and/or sensitive to the local soil environment. Soil water δ⁴⁴Ca values appear to increase with depth to at least 35 cm below the top of the forest floor, which is close to the deepest level of fine roots. The heavy plant fractionated signature of Ca in the finely rooted upper soils filters downward where it is retained on ion exchange sites, leached into groundwater, and discharged into surface waters.The relationship between Ca uptake by tree fine roots and the pattern of δ⁴⁴Ca enrichment with soil depth was modeled for two Ca pools: the forest floor (litter) and the underlying (upper B) mineral soil. Six study plots were investigated along two hillside toposequences trending upwards from a first order stream. We used allometric equations describing the Ca distribution in boreal tree species to calculate weighted average δ⁴⁴Ca values for the stands in each plot and estimate Ca uptake rates. The δ⁴⁴Ca value of precipitation was measured, and soil weathering signatures deduced, by acid leaching of lower B mineral soils. Steady state equations were used to derive a set of model Ca fluxes and fractionation factors for each plot. The model reproduces the increase in δ⁴⁴Ca with depth found in forest floor and upper B soil waters. Transient model runs show that the forest Ca cycle is sensitive to changes in plant Ca uptake rate, such as would occur during ontogeny or disturbance. Accordingly, secular records of δ⁴⁴Ca in tree ring cellulose have the potential to monitor changes in the forest Ca cycle through time, thus providing a new tool for evaluating natural and anthropogenic impacts on forest health. Another model run shows that by changing the size of the isotope fractionation factor and adjusting for differences in forest productivity, that the range in Ca isotope fractionation in forested ecosystems reported in the literature, thus far, is reproduced. As a quantitative tool, the Ca cycling model produces a reasonable set of relative Ca fluxes for the La Ronge site, consistent with Environment Canada’s measurements for wet deposition in the region and simulated Ca release from soil mineral weathering using the PROFILE model. But the sensitivity of the model is limited by the small range of fractionation observed in this boreal shield setting of ∼1‰, which limits accuracy. If the model were applied to a site with a greater range in δ⁴⁴Ca values among the principal Ca fluxes, it is capable of producing robust and reliable estimations of Ca fluxes that are otherwise difficult to measure in forested ecosystems.     

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

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

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

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

Interpreting the Ca isotope record of marine biogenic carbonates

An 18 million year record of the Ca isotopic composition (δ^44/42Ca) of planktonic foraminiferans from ODP site 925, in the Atlantic, on the Ceara Rise, provides the opportunity for critical analysis of Ca isotope-based reconstructions of the Ca cycle. δ^44/42Ca in this record averages +0.37 ± 0.05 (1σ SD) and ranges from +0.21‰ to +0.52‰. The record is a good match to previously published Neogene Ca isotope records based on foraminiferans, but is not similar to the record based on bulk carbonates, which has values that are as much as 0.25‰ lower. Bulk carbonate and planktonic foraminiferans from core tops differ slightly in their δ^44/42Ca (i.e., by 0.06 ± 0.06‰ (n = 5)), while the difference between bulk carbonate and foraminiferan values further back in time is markedly larger, leaving open the question of the cause of the difference. Modeling the global Ca cycle from downcore variations in δ^44/42Ca by assuming fixed values for the isotopic composition of weathering inputs (δ^44/42Ca_w) and for isotope fractionation associated with the production of carbonate sediments (Δ_sed) results in unrealistically large variations in the total mass of Ca²⁺ in the oceans over the Neogene. Alternatively, variations of ±0.05‰ in the Ca isotope composition of weathering inputs or in the extent of fractionation of Ca isotopes during calcareous sediment formation could entirely account for variations in the Ca isotopic composition of marine carbonates. Ca isotope fractionation during continental weathering, such as has been recently observed, could easily result in variations in δ^44/42Ca_w of a few tenths of permil. Likewise a difference in the fractionation factors associated with aragonite versus calcite formation could drive shifts in Δ_sed of tenths of permil with shifts in the relative output of calcite and aragonite from the ocean. Until better constraints on variations in δ^44/42Ca_w and Δ_sed have been established, modeling the Ca²⁺ content of seawater from Ca isotope curves should be approached cautiously.     

Germanium isotope measurements of high-temperature geothermal fluids using double-spike hydride generation MC-ICP-MS

We present a double-spike isotope dilution MC-ICP-MS technique for the determination of germanium (Ge) isotope fractionation. Using this technique we determined Ge isotope compositions of geothermal spring fluids, a Columbia River Basalt sample, and an in-house diatom standard. Our technique uses a ⁷³Ge/⁷⁰Ge double spike in combination with hydride generation for Ge extraction from the sample matrix. Fractionation is determined on the ⁷⁴Ge/⁷²Ge mass ratio. The double spike allows us to effectively correct analytical isotope fractionation. Our external standard reproducibility is 0.4‰ (2 SD) over the course of several months. The minimum quantity of Ge needed for isotope analysis is approximately 2 ng. Consistent with previous work on geothermal fluids, Ge in the geothermal spring samples presented here is enriched over Si as compared to low temperature weathering signatures. This observation is typically interpreted as Ge exclusion during silicate mineral precipitation (e.g., quartz). Our isotope results indicate that the analyzed high temperature fluids fractionate Ge isotopes with a range in δ⁷⁴Ge between −0.4‰ and −1.4‰ relative to a Columbia River basalt. We cautiously interpret the observed fractionation as preferential removal of heavy Ge isotopes out of solution during cooling of the hydrothermal fluid and subsequent precipitation of quartz.     

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.     

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

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

Iron isotope fractionation during hydrothermal ore deposition and alteration

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

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

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