Application of iron and zinc isotopes to track the sources and mechanisms of metal loading in a mountain watershed

Here the hydrogeochemical constraints of a tracer dilution study are combined with Fe and Zn isotopic measurements to pinpoint metal loading sources and attenuation mechanisms in an alpine watershed impacted by acid mine drainage. In the tested mountain catchment, δ⁵⁶Fe and δ⁶⁶Zn isotopic signatures of filtered stream water samples varied by ∼3.5‰ and 0.4‰, respectively. The inherent differences in the aqueous geochemistry of Fe and Zn provided complimentary isotopic information. For example, variations in δ⁵⁶Fe were linked to redox and precipitation reactions occurring in the stream, while changes in δ⁶⁶Zn were indicative of conservative mixing of different Zn sources. Fen environments contributed distinctively light dissolved Fe (<−2.0‰) and isotopically heavy suspended Fe precipitates to the watershed, while Zn from the fen was isotopically heavy (>+0.4‰). Acidic drainage from mine wastes contributed heavier dissolved Fe (∼+0.5‰) and lighter Zn (∼+0.2‰) isotopes relative to the fen. Upwelling of Fe-rich groundwater near the mouth of the catchment was the major source of Fe (δ⁵⁶Fe ∼ 0‰) leaving the watershed in surface flow, while runoff from mining wastes was the major source of Zn. The results suggest that given a strong framework for interpretation, Fe and Zn isotopes are useful tools for identifying and tracking metal sources and attenuation mechanisms in mountain watersheds.     

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.     

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

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

Iron isotope fractionation in subterranean estuaries

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

Iron isotope variations in Holocene sediments of the Gotland Deep, Baltic Sea

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

SIMS analyses of silicon and oxygen isotope ratios for quartz from Archean and Paleoproterozoic banded iron formations

Banded iron formations (BIFs) are chemical marine sediments dominantly composed of alternating iron-rich (oxide, carbonate, sulfide) and silicon-rich (chert, jasper) layers. Isotope ratios of iron, carbon, and sulfur in BIF iron-bearing minerals are biosignatures that reflect microbial cycling for these elements in BIFs. While much attention has focused on iron, banded iron formations are equally banded silica formations. Thus, silicon isotope ratios for quartz can provide insight on the sources and cycling of silicon in BIFs. BIFs are banded by definition, and microlaminae, or sub-mm banding, are characteristic of many BIFs. In situ microanalysis including secondary ion mass spectrometry is well-suited for analyzing such small features. In this study we used a CAMECA IMS-1280 ion microprobe to obtain highly accurate (±0.3‰) and spatially resolved (∼10 μm spot size) analyses of silicon and oxygen isotope ratios for quartz from several well known BIFs: Isua, southwest Greenland (∼3.8 Ga); Hamersley Group, Western Australia (∼2.5 Ga); Transvaal Group, South Africa (∼2.5 Ga); and Biwabik Iron Formation, Minnesota, USA (∼1.9 Ga). Values of δ¹⁸O range from +7.9‰ to +27.5‰ and include the highest reported δ¹⁸O values for BIF quartz. Values of δ³⁰Si have a range of ∼5‰ from −3.7‰ to +1.2‰ and extend to the lowest δ³⁰Si values for Precambrian cherts. Isua BIF samples are homogeneous in δ¹⁸O to ±0.3‰ at mm- to cm-scale, but are heterogeneous in δ³⁰Si up to 3‰, similar to the range in δ³⁰Si found in BIFs that have not experienced high temperature metamorphism (up to 300 °C). Values of δ³⁰Si for quartz are homogeneous to ±0.3‰ in individual sub-mm laminae, but vary by up to 3‰ between multiple laminae over mm-to-cm of vertical banding. The scale of exchange for Si in quartz in BIFs is thus limited to the size of microlaminae, or less than ∼1 mm. We interpret differences in δ³⁰Si between microlaminae as preserved from primary deposition. Silicon in BIF quartz is mostly of marine hydrothermal origin (δ³⁰Si < −0.5‰) but silicon from continental weathering (δ³⁰Si ∼ 1‰) was an important source as early as 3.8 Ga.     

Iron and magnesium isotopic compositions of peridotite xenoliths from Eastern China

We present high-precision measurements of Mg and Fe isotopic compositions of olivine, orthopyroxene (opx), and clinopyroxene (cpx) for 18 lherzolite xenoliths from east central China and provide the first combined Fe and Mg isotopic study of the upper mantle. δ⁵⁶Fe in olivines varies from 0.18‰ to −0.22‰ with an average of −0.01 ± 0.18‰ (2SD, n = 18), opx from 0.24‰ to −0.22‰ with an average of 0.04 ± 0.20‰, and cpx from 0.24‰ to −0.16‰ with an average of 0.10 ± 0.19‰. δ²⁶Mg of olivines varies from −0.25‰ to −0.42‰ with an average of −0.34 ± 0.10‰ (2SD, n = 18), opx from −0.19‰ to −0.34‰ with an average of −0.25 ± 0.10‰, and cpx from −0.09‰ to −0.43‰ with an average of −0.24 ± 0.18‰. Although current precision (∼±0.06‰ for δ⁵⁶Fe; ±0.10‰ for δ²⁶Mg, 2SD) limits the ability to analytically distinguish inter-mineral isotopic fractionations, systematic behavior of inter-mineral fractionation for both Fe and Mg is statistically observed: Δ⁵⁶Fe_ol-cpx = −0.10 ± 0.12‰ (2SD, n = 18); Δ⁵⁶Fe_ol-opx = −0.05 ± 0.11‰; Δ²⁶Mg_ol-opx = −0.09 ± 0.12‰; Δ²⁶Mg_ol-cpx = −0.10 ± 0.15‰. Fe and Mg isotopic composition of bulk rocks were calculated based on the modes of olivine, opx, and cpx. The average δ⁵⁶Fe of peridotites in this study is 0.01 ± 0.17‰ (2SD, n = 18), similar to the values of chondrites but slightly lower than mid-ocean ridge basalts (MORB) and oceanic island basalts (OIB). The average δ²⁶Mg is −0.30 ± 0.09‰, indistinguishable from chondrites, MORB, and OIB. Our data support the conclusion that the bulk silicate Earth (BSE) has chondritic δ⁵⁶Fe and δ²⁶Mg.The origin of inter-mineral fractionations of Fe and Mg isotopic ratios remains debated. δ⁵⁶Fe between the main peridotite minerals shows positive linear correlations with slopes within error of unity, strongly suggesting intra-sample mineral-mineral Fe and Mg isotopic equilibrium. Because inter-mineral isotopic equilibrium should be reached earlier than major element equilibrium via chemical diffusion at mantle temperatures, Fe and Mg isotope ratios of coexisting minerals could be useful tools for justifying mineral thermometry and barometry on the basis of chemical equilibrium between minerals. Although most peridotites in this study exhibit a narrow range in δ⁵⁶Fe, the larger deviations from average δ⁵⁶Fe for three samples likely indicate changes due to metasomatic processes. Two samples show heavy δ⁵⁶Fe relative to the average and they also have high La/Yb and total Fe content, consistent with metasomatic reaction between peridotite and Fe-rich and isotopically heavy melt. The other sample has light δ⁵⁶Fe and slightly heavy δ²⁶Mg, which may reflect Fe-Mg inter-diffusion between peridotite and percolating melt.     

Exploration potential of Cu isotope fractionation in porphyry copper deposits

We examined the copper isotope ratio of primary high temperature Cu-sulfides, secondary low temperature Cu-sulfides (and Cu-oxides) as well as Fe-oxides in the leach cap, which represent the weathered remains of a spectrum of Cu mineralization, from nine porphyry copper deposits. Copper isotope ratios are reported as δ⁶⁵Cu‰ = ((⁶⁵Cu/⁶³Cu_sample/⁶⁵Cu/⁶³Cu_{NIST 976 standard}) − 1) ? 10³. Errors for all the analyses are ± 0.14‰ (determined by multiple analyses of the samples) and mass bias was corrected through standard-sample-standard bracketing. The overall isotopic variability measured in these samples range from − 16.96‰ to 9.98‰.     

The continental shelf benthic iron flux and its isotope composition

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

Silicon isotopes in meteorites and planetary core formation

The silicon (Si) isotope compositions of 42 meteorite and terrestrial samples have been determined using MC-ICPMS with the aim of resolving the current debate over their compositions and the implications for core formation. No systematic δ³⁰Si differences are resolved between chondrites (δ³⁰Si = −0.49 ± 0.15‰, 2σ_SD) and achondrites (δ³⁰Si = −0.47 ± 0.11‰, 2σ_SD), although enstatite chondrites are consistently lighter (δ³⁰Si = −0.63 ± 0.07‰, 2σ_SD) in comparison to other meteorite groups. The data reported here for meteorites and terrestrial samples display an average difference Δ³⁰Si_{BSE−meteorite∗} = 0.15 ± 0.10‰, which is consistent within uncertainty with previous studies. No effect from sample heterogeneity, preparation, chemistry or mass spectrometry can be identified as responsible for the reported differences between current datasets. The heavier composition of the bulk silicate Earth is consistent with previous conclusions that Si partitioned into the metal phase during metal-silicate equilibration at the time of core formation. Fixing the temperature of core formation to the peridotite liquidus and using an appropriate metal silicate fractionation factor (ε ∼0.89), the Δ³⁰Si_{BSE−meteorite∗} value from this study indicates that the Earth core contains at least 2.5 and possibly up to 16.8 wt% Si.     

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.     

Diffusion induced Li isotopic fractionation during the cooling of magmatic rocks: The case of pyroxene phenocrysts from nakhlite meteorites

Ion-microprobe was used to measure Li abundances and isotopic compositions in pyroxenes from three Martian meteorites belonging to the nakhlite family. The profiles performed across augite crystals from Northwest Africa 817 show a large isotopic zoning from crystal cores (δ⁷Li ∼ 0‰) to rims (δ⁷Li ∼ +20‰) while Li abundances are almost constant (∼9.2 μg/g). Unlike NWA 817, the pyroxene studied in the Miller Range 03346 nakhlite shows a zoning in Li abundance, with concentrations increasing from ∼2.5 μg/g in the core to ∼9 μg/g in the rim. The augite rim (δ⁷Li = +7‰) is slightly enriched in ⁷Li with regard to the core (δ⁷Li = +4‰), but most of the isotopic variations observed occur at an intermediate position along the profile, where δ⁷Li falls down to ∼−11‰. In the case of Nakhla, Li concentrations in augite increase from cores (∼3.5 μg/g) to rims (∼6.5 μg/g), while the δ⁷Li variation is restricted (i.e., between δ⁷Li = +6.0 and +12.6‰). For the three meteorites the Li abundances were also measured in the groundmass, which was found to be enriched in lithium (∼10 μg/g). Conventional magmatic and post-magmatic processes such as alteration and fractional crystallization, fail to explain the dataset obtained on nakhlites. Degassing processes, which were previously proposed to explain the Li distribution in shergottite crystals, cannot result in the strong decoupling between Li abundances and isotopic composition observed in nakhlites. We suggest that the original magmatic Li distributions (concentrations and isotopic compositions) in nakhlites have been modified by diffusion of Li from the Li-rich groundmass towards the pyroxene crystals during sub-solidus cooling. Diffusion appears to have been efficient for NWA 817 and MIL 03346 but, apparently, did not produce a significant migration of Li in Nakhla, possibly because of the lower abundance of groundmass in the latter. Diffusion induced Li redistributions may also affect terrestrial porphyric rocks but very specific cooling rates are required to quench the diffusion profiles as observed in two of the present nakhlites.     

Variations of δSi and Ge/Si with weathering and biogenic input in tropical basaltic ash soils under monoculture

In soils, silicon released by mineral weathering can be retrieved from soil solution through clay formation, Si adsorption onto secondary oxides and plant uptake, thereby impacting the Si-isotopic signature and Ge/Si ratio of dissolved Si (DSi) exported to rivers. Here we use these proxies to study the contribution of biogenic Si (BSi) in a soil-plant system involving basaltic ash soils differing in weathering degree under intensive banana cropping. δ³⁰Si and Ge/Si ratios were determined in bulk soils (<2 mm), sand (50-2000 μm), silt (2-50 μm), amorphous Si (ASi, 2-50 μm) and clay (<2 μm) fractions: δ³⁰Si by MC-ICP-MS Nu Plasma in medium resolution, operating in dry plasma with Mg doping (δ³⁰Si vs. NBS28 ± 0.12‰ ± 2σ_SD), Ge/Si computed after determination of Ge and Si concentrations by HR-ICP-MS and ICP-AES, respectively. Components of the ASi fraction were quantified by microscopic counting (phytoliths, diatoms, ashes). Compared to fresh ash (δ³⁰Si = −0.38‰; Ge/Si = 2.21 μmol mol⁻¹), soil clay fractions (<2 μm) were enriched in light Si isotopes and Ge: with increasing weathering degree, δ³⁰Si decreased from −1.19 to −2.37‰ and Ge/Si increased from 4.10 to 5.25 μmol mol⁻¹. Sand and silt fractions displayed δ³⁰Si values close to fresh ash (−0.33‰) or higher due to saharian dust quartz deposition, whose contribution was evaluated by isotopic mass balance calculation. Si-isotopic signatures of bulk soils (<2 mm) were strongly governed by the relative proportions of primary and secondary minerals: the bulk soil Si-isotopic budget could be closed indicating that all the phases involved were identified. Microscopic counting highlighted a surface accumulation of banana phytoliths and a stable phytolith pool from previous forested vegetation. δ³⁰Si and Ge/Si values of clay fractions in poorly developed volcanic soils, isotopically heavier and Ge-depleted in surface horizons, support the occurrence of a DSi source from banana phytolith dissolution, available for Si sequestration in clay-sized secondary minerals (clay minerals formation and Si adsorption onto Fe-oxide). In the soil-plant system, δ³⁰Si and Ge/Si are thus highly relevant to trace weathering and input of DSi from phytoliths in secondary minerals, although not quantifying the net input of BSi to DSi.     

The origins and behaviour of carbon in a major semi-arid river,the Murray River,Australia, as constrained by carbon isotopes and hydrochemistry

δ¹³C values of dissolved inorganic C (DIC), dissolved organic C (DOC), and particulate organic C (POC) together with δ¹⁸O and δ²H values of water, δ³⁴S values of dissolved SO₄, and major ion concentrations were measured in the Murray River and its tributaries between November 2005 and April 2007 to constrain the origins and behaviour of riverine C. δ¹³C_DIC values in the Murray River vary between −9.5 and −4.7‰ with a range of <3‰ within any sampling round. δ¹³C_DIC values of the tributaries are −11.0‰ to −5.1‰. DIC concentrations of the Murray River increase from ∼25 mg/L in the middle and upper reaches of the river to 45–55 mg/L in the lower reaches. However, the mass ratio of DIC as a proportion of the total dissolved solids (TDS) decreases from ∼0.6–0.7 in the headwaters to ∼0.2–0.3 in the lower reaches of the river, with similar downstream changes in DIC/Cl ratios. This precludes simple evaporative concentration of DIC and is interpreted as the river evading CO₂; this interpretation is consistent with pCO₂ values that are in the range 550–11,200 ppm volume (ppmv), which are far higher than those in equilibrium with the atmosphere (∼360 ppmv). The δ¹³C_DIC values are similar to those that would be produced by the weathering of marine limestone (δ¹³C ∼ 0‰). However, the lack of marine limestones cropping out in the Murray–Darling Basin and the relatively uniform δ¹³C_DIC values of the Murray River (even in upland reaches where the dominant rock types are metamorphosed silicates and granites) make this unlikely. Rather the high pCO₂ values and δ¹³C_DIC values are best explained by a combination of mineralisation of low δ¹³C organic C and evasion to the atmosphere. The rate of these two processes may attain near steady state and control both DIC concentrations and δ¹³C values.     

Indoor and outdoor urban atmospheric CO2: Stable carbon isotope constraints on mixing and mass balance

From July to November 2009, concentrations of CO₂ in 78 samples of ambient air collected in 18 different interior spaces on a university campus in Dallas, Texas (USA) ranged from 386 to 1980 ppm. Corresponding δ¹³C values varied from −8.9‰ to −19.4‰. The CO₂ from 22 samples of outdoor air (also collected on campus) had a more limited range of concentrations from 385 to 447 ppm (avg. = 408 ppm), while δ¹³C values varied from −10.1‰ to −8.4‰ (avg.=-9.0‰). In contrast to ambient indoor and outdoor air, the concentrations of CO₂ exhaled by 38 different individuals ranged from 38,300 to 76,200 ppm (avg. = 55,100 ppm), while δ¹³C values ranged from −24.8‰ to −17.7‰ (avg. = −21.8‰). The residence times of the total air in the interior spaces of this study appear to have been on the order of 10 min with relatively rapid approaches (∼30 min) to steady-state concentrations of ambient CO₂ gas. Collectively, the δ¹³C values of the indoor CO₂ samples were linearly correlated with the reciprocal of CO₂ concentration, exhibiting an intercept of −21.8‰, with r² = 0.99 and p < 0.001 (n = 78). This high degree of linearity for CO₂ data representing 18 interior spaces (with varying numbers of occupants), and the coincidence of the intercept (−21.8‰) with the average δ¹³C value for human-exhaled CO₂ demonstrates simple mixing between two inputs: (1) outdoor CO₂ introduced to the interior spaces by ventilation systems, and (2) CO₂ exhaled by human occupants of those spaces. If such simple binary mixing is a common feature of interior spaces, it suggests that the intercept of a mixing line defined by two data points (CO₂ input from the local ventilation system and CO₂ in the ambient air of the room) could be a reasonable estimate of the average δ¹³C value of the CO₂ exhaled by the human occupants. Thus, such indoor spaces appear to constitute effective “sample vessels” for collection of CO₂ that can be used to determine the average proportions of C₃ and C₄-derived C in the diets of the occupants. For the various groups occupying the rooms sampled in this study, C₄-derived C appears to have constituted ∼40% of the average diet.     

Carbon isotope effects in the open-system Fischer-Tropsch synthesis

Carbon isotopic composition was measured for products of the Fischer-Tropsch synthesis: catalytic reaction between CO and H₂ to produce CO, CO₂, light hydrocarbons C1-C4 and “oil” fraction. Hydrogen isotopes were also measured in the oil fraction and the produced water. Experimental runs were conducted in the flow-through reactor at 260-310 °C and 30 bar using the synthesis gas composed of 5N₂ + 3H₂ + 2CO, on Fe-catalyst mixed with ZSM-5 synthetic zeolite. In the two of seven runs a Fe + Co-catalyst was used that gives a lower yield of unsaturated hydrocarbons in reaction products. The isotopic effects depended on the conversion of the carbon monoxide. Under steady-state conditions (CO conversion more than 90%) a strong kinetic fractionation was observed between CO and CO₂ (∼−10‰) and CO and hydrocarbons (∼+38‰). At low conversion a clear “inverse” isotopic trend of the depletion in ¹³C of longer hydrocarbon chains was observed. On average, Δ₁₂ = δ¹³C(CH₄) − δ¹³C(C₂H₆) correlates well with the CO conversion: the C₂H₆ is ∼6‰ isotopically lighter than CH₄ at low conversion and ∼2‰ heavier at steady-state regime. Under steady-state conditions there almost no difference was observed in the isotopic composition of methane and ethane and higher hydrocarbons. The chemical composition of light hydrocarbons in the products of flow-through, dynamic FTS is different from that found in the static FTS-type experiments with Fe-catalyst, but isotopic effects are similar. Our results suggest that the isotopic distribution of carbon found in so-called “abiogenic” hydrocarbons from some natural gases (δ¹³C1 > δ¹³C2 > δ¹³C3  >?) is somewhat similar to that at low conversion of CO, but do not resemble the distribution characteristic for the high conversion products, at least, on Fe-catalyst. Other processes (a simple mixing of two or more endmembers) or other P-T conditions of the carbon reduction could be responsible for the “inverse” isotopic trend found in meteorites and some natural gases.     

Oxygen and hydrogen isotopes of rainfall and dripwater at DeSoto Caverns (Alabama, USA): Key to understanding past variability of moisture transport from the Gulf of Mexico

The Southeast and the US Gulf Coast in particular are notably lacking isotope data in the water cycle despite the fact that moisture transport from the Gulf of Mexico (GOM) has a considerable influence on both regional and continental rainfall patterns. This study reports time-series of oxygen and hydrogen isotopes acquired over a 3-year period (2005-2008) from GOM-derived rainfall, cave dripwater and shallow groundwaters, and offers valuable insights on the links between factors controlling regional rainfall and the ubiquitous karst hydrology.Amount-weighted mean monthly rainwater δ¹⁸O and δD values in Tuscaloosa, Alabama range from −1.5 to −8.3‰ and −1.2 to −49.5‰, respectively, and show mean seasonal amplitudes of ∼4‰ (δ¹⁸O) and ∼25‰ (δD). In comparison d-excess values display large seasonal amplitudes of 10-20‰ resulting from differences in the degree of evaporation from falling raindrops between summer and winter months, and correlate well with the coeval air temperature (r² = 0.59; p < 0.05). Deviations of the Gulf Coast Meteoric Water Line (GCMWL) slope and d-excess from the global meteoric water line (GMWL) are attributed to different rates of evaporation after condensation, and to humidity contrast between the cloud boundary layer and the surrounding atmosphere in the vapor source region, respectively. Rainfall amounts declined during the study interval from an excessive “wet” year, ascribed to six tropical storms incursions during an unusually active hurricane season in 2005, to an onset of a regional drought during 2007-2008 with monthly rainfall amounts substantially below normal values (30-year monthly means). An interannual trend of ¹⁸O and ²H enrichments is discerned from 2005 to 2008 (1.4‰ and 11.6‰, respectively) coeval with the decline in rainfall amounts.Dripwater samples from nearby DeSoto Caverns show weak δ¹⁸O and δD seasonal variations and record only 20% and 51% of the ¹⁸O and ²H enrichments, respectively, discerned in the rainwater 3-year time-trends. The seasonal and interannual amplitude attenuations in the dripwaters are attributed to a relatively thick overlying bedrock (∼30-40 m) and a relatively large, well-mixed, epikarst-storage reservoir. Residence time of water in the cave’s epikarst is estimated to be 1-3 months based on high-resolution flow-rate data.Our investigation suggests that global atmospheric circulation patterns (ENSO and Bermuda High) likely govern the interannual δ¹⁸O and δD isotope trends discerned in the water cycle compartments but much longer time-series are required to confirm our conjectures. The results of this study form a solid basis for future acquisition and interpretation of climate proxy records from regional speleothems.     

Nature of volatile depletion and genetic relationships in enstatite chondrites and aubrites inferred from Zn isotopes

Enstatite meteorites include the undifferentiated enstatite chondrites and the differentiated enstatite achondrites (aubrites). They are the most reduced group of all meteorites. The oxygen isotope compositions of both enstatite chondrites and aubrites plot along the terrestrial mass fractionation line, which suggests some genetic links between these meteorites and the Earth as well.For this study, we measured the Zn isotopic composition of 25 samples from the following groups: aubrites (main group and Shallowater), EL chondrites, EH chondrites and Happy Canyon (impact-melt breccia). We also analyzed the Zn isotopic composition and elemental abundance in separated phases (metal, silicates, and sulfides) of the EH4, EL3, and EL6 chondrites. The different groups of meteorites are isotopically distinct and give the following values (‰): aubrite main group (−7.08 < δ⁶⁶Zn < −0.37); EH3 chondrites (0.15 < δ⁶⁶Zn < 0.31); EH4 chondrites (0.15 < δ⁶⁶Zn < 0.27); EH5 chondrites (δ⁶⁶Zn = 0.27 ± 0.09; n = 1); EL3 chondrites (0.01 < δ⁶⁶Zn < 0.63); the Shallowater aubrite (1.48 < δ⁶⁶Zn < 2.36); EL6 chondrites (2.26 < δ⁶⁶Zn < 7.35); and the impact-melt enstatite chondrite Happy Canyon (δ⁶⁶Zn = 0.37).The aubrite Peña Blanca Spring (δ⁶⁶Zn = −7.04‰) and the EL6 North West Forrest (δ⁶⁶Zn = 7.35‰) are the isotopically lightest and heaviest samples, respectively, known so far in the Solar System. In comparison, the range of Zn isotopic composition of chondrites and terrestrial samples (−1.5 < δ⁶⁶Zn < 1‰) is much smaller ( [Luck et al., 2005] and [Herzog et al., 2009]).EH and EL3 chondrites have the same Zn isotopic composition as the Earth, which is another example of the isotopic similarity between Earth and enstatite chondrites. The Zn isotopic composition and abundance strongly support that the origin of the volatile element depletion between EL3 and EL6 chondrites is due to volatilization, probably during thermal metamorphism. Aubrites show strong elemental depletion in Zn compared to both EH and EL chondrites and they are enriched in light isotopes (δ⁶⁶Zn down to −7.04‰). This is the opposite of what would be expected if Zn elemental depletion was due to evaporation, assuming the aubrites started with an enstatite chondrite-like Zn isotopic composition. Evaporation is therefore not responsible for volatile loss from aubrites. On Earth, Zn isotopes fractionate very little during igneous processes, while differentiated meteorites show only minimal Zn isotopic variability. It is therefore very unlikely that igneous processes can account for the large isotopic fractionation of Zn in aubrites. Condensation of an isotopically light vapor best explains Zn depletion and isotopically light Zn in these puzzling rocks. Mass balance suggests that this isotopically light vapor carries Zn lost by the EL6 parent body during thermal metamorphism and that aubrites evolved from an EL6-like parent body. Finally, Zn isotopes suggest that Shallowater and aubrites originate from distinct parent bodies.     

Oxygen- and magnesium-isotope compositions of calcium-aluminum-rich inclusions from Rumuruti (R) chondrites

We report oxygen- and magnesium-isotope compositions of Ca,Al-rich inclusions (CAIs) from several Rumuruti (R) chondrites measured in situ using a Cameca ims-1280 ion microprobe. On a three-isotope oxygen diagram, δ¹⁷O vs. δ¹⁸O, compositions of individual minerals in most R CAIs analyzed fall along a slope-1 line. Based on the variations of Δ¹⁷O values (Δ¹⁷O = δ¹⁷O − 0.52 × δ¹⁸O) within individual inclusions, the R CAIs are divided into (i) ¹⁶O-rich (Δ¹⁷O ∼ −23-26‰), (ii) uniformly ¹⁶O-depleted (Δ¹⁷O ∼ −2‰), and (iii) isotopically heterogeneous (Δ¹⁷O ranges from −25‰ to +5‰). One of the hibonite-rich CAIs, H030/L, has an intermediate Δ¹⁷O value of −12‰ and a highly fractionated composition (δ¹⁸O ∼ +47‰). We infer that like most CAIs in other chondrite groups, the R CAIs formed in an ¹⁶O-rich gaseous reservoir. The uniformly ¹⁶O-depleted and isotopically heterogeneous CAIs subsequently experienced oxygen-isotope exchange during remelting in an ¹⁶O-depleted nebular gas, possibly during R chondrite chondrule formation, and/or during fluid-assisted thermal metamorphism on the R chondrite parent asteroid.Three hibonite-bearing CAIs and one spinel-plagioclase-rich inclusion were analyzed for magnesium-isotope compositions. The CAI with the highly fractionated oxygen isotopes, H030/L, shows a resolvable excess of ²⁶Mg (²⁶Mg^∗) corresponding to an initial ²⁶Al/²⁷Al ratio of ∼7 × 10⁻⁷. Three other CAIs show no resolvable excess of ²⁶Mg (²⁶Mg^∗). The absence of ²⁶Mg^∗ in the spinel-plagioclase-rich CAI from a metamorphosed R chondrite NWA 753 (R3.9) could have resulted from metamorphic resetting. Two other hibonite-bearing CAIs occur in the R chondrites (NWA 1476 and NWA 2446), which appear to have experienced only minor degrees of thermal metamorphism. These inclusions could have formed from precursors with lower than canonical ²⁶Al/²⁷Al ratio.