Non-biological fractionation of stable Ca isotopes in soils of the Atacama Desert, Chile

We measured Ca stable isotope ratios (δ^44/40Ca) in an ancient (2 My), hyperarid soil where the primary source of mobile Ca is atmospheric deposition. Most of the Ca in the upper meter of this soil (3.5 kmol m⁻²) is present as sulfates (2.5 kmol m⁻²), and to a lesser extent carbonates (0.4 kmol m⁻²). In aqueous extracts of variably hydrated calcium sulfate minerals, δ^44/40Ca_E values (vs. bulk Earth) increase with depth (1.4 m) from a minimum of −1.91‰ to a maximum of +0.59‰. The trend in carbonate-δ^44/40Ca in the top six horizons resembles that of sulfate-δ^44/40Ca, but with values 0.1-0.6‰ higher. The range of observed Ca isotope values in this soil is about half that of δ^44/40Ca values observed on Earth. Linear correlation among δ^44/40Ca, δ³⁴S and δ¹⁸O values indicates either (a) a simultaneous change in atmospheric input values for all three elements over time, or (b) isotopic fractionation of all three elements during downward transport. We present evidence that the latter is the primary cause of the isotopic variation that we observe. Sulfate-δ³⁴S values are positively correlated with sulfate-δ¹⁸O values (R² = 0.78) and negatively correlated with sulfate δ^44/40Ca_E values (R² = 0.70). If constant fractionation and conservation of mass with downward transport are assumed, these relationships indicate a δ^44/40Ca fractionation factor of −0.4‰ in CaSO₄. The overall depth trend in Ca isotopes is reproduced by a model of isotopic fractionation during downward Ca transport that considers small and infrequent but regularly recurring rainfall events. Near surface low Ca isotope values are reproduced by a Rayleigh model derived from measured Ca concentrations and the Ca fractionation factor predicted by the relationship with S isotopes. This indicates that the primary mechanism of stable isotope fractionation in CaSO₄ is incremental and effectively irreversible removal of an isotopically enriched dissolved phase by downward transport during small rainfall events.     

Quadruple sulfur isotope analysis of ca. 3.5 Ga Dresser Formation: New evidence for microbial sulfate reduction in the early Archean

Multiple sulfur isotope system is a powerful new tracer for atmospheric, volcanic, and biological influences on sulfur cycles in the anoxic early Earth. Here, we report high-precision quadruple sulfur isotope analyses (³²S/³³S/³⁴S/³⁶S) of barite, pyrite in barite, and sulfides in related hydrothermal and igneous rocks occurring in the ca. 3.5 Ga Dresser Formation, Western Australia. Our results indicate that observed isotopic variations are mainly controlled by mixing of mass-dependently (MD) and non-mass-dependently fractionated (non-MD) sulfur reservoirs. Based on the quadruple sulfur isotope systematics (δ³⁴S-Δ³³S-Δ³⁶S) for these minerals, four end-member sulfur reservoirs have been recognized: (1) non-MD sulfate (δ³⁴S = −5 ± 2‰; Δ³³S = −3 ± 1‰); (2) MD sulfate (δ³⁴S = +10 ± 3‰); (3) non-MD sulfur (δ³⁴S > +6‰; Δ³³S > +4‰); and (4) igneous MD sulfur (δ³⁴S = Δ³³S = 0‰). The first and third components show a clear non-MD signatures, thus probably represent sulfate and sulfur aerosol inputs. The MD sulfate component (2) is enriched in ³⁴S (+10 ± 3‰) and may have originated from microbial and/or abiotic disproportionation of volcanic S or SO₂. Our results reconfirm that the Dresser barites contain small amounts of pyrite depleted in ³⁴S by 15-22‰ relative to the host barite. These barite-pyrite pairs exhibit a mass-dependent relationship of δ³³S/δ³⁴S with slope less than 0.512, which is consistent with that expected for microbial sulfate reduction and is significantly different from that of equilibrium fractionation (0.515). The barite-pyrite pairs also show up to 1‰ difference in Δ³⁶S values and steep Δ³⁶S/Δ³³S slopes, which deviate from the main Archean array (Δ³⁶S/Δ³³S = −0.9) and are comparable to isotope effects exhibited by sulfate reducing microbes (Δ³⁶S/Δ³³S = −5 to −11). These new lines of evidence support the existence of sulfate reducers at ca. 3.5 Ga, whereas microbial sulfur disproportionation may have been more limited than recently suggested.     

Variations of Li and Mg isotope ratios in bulk chondrites and mantle xenoliths

We present whole rock Li and Mg isotope analyses of 33 ultramafic xenoliths from the terrestrial mantle, which we compare with analyses of 30 (mostly chondritic) meteorites. The accuracy of our new Mg isotope ratio measurement protocol is substantiated by a combination of standard addition experiments, the absence of mass independent effects in terrestrial samples and our obtaining identical values for rock standards using two different separation chemistries and three different mass-spectrometric introduction systems. Carbonaceous, ordinary and enstatite chondrites have irresolvable mean stable Mg isotopic compositions (δ²⁵Mg = −0.14 ± 0.06; δ²⁶Mg = −0.27 ± 0.12‰, 2SD), but our enstatite chondrite samples have lighter δ⁷Li (by up to ∼3‰) than our mean carbonaceous and ordinary chondrites (3.0 ± 1.5‰, 2SD), possibly as a result of spallation in the early solar system. Measurements of equilibrated, fertile peridotites give mean values of δ⁷Li = 3.5 ± 0.5‰, δ²⁵Mg = −0.10 ± 0.03‰ and δ²⁶Mg = −0.21 ± 0.07‰. We believe these values provide a useful estimate of the primitive mantle and they are within error of our average of bulk carbonaceous and ordinary chondrites. A fuller range of fresh, terrestrial, ultramafic samples, covering a variety of geological histories, show a broad positive correlation between bulk δ⁷Li and δ²⁶Mg, which vary from −3.7‰ to +14.5‰, and −0.36‰ to + 0.06‰, respectively. Values of δ⁷Li and δ²⁶Mg lower than our estimate of primitive mantle are strongly linked to kinetic isotope fractionation, occurring during transport of the mantle xenoliths. We suggest Mg and Li diffusion into the xenoliths is coupled to H loss from nominally anhydrous minerals following degassing. Diffusion models suggest that the co-variation of Mg and Li isotopes requires comparable diffusivities of Li and Mg in olivine. The isotopically lightest samples require ∼5-10 years of diffusive ingress, which we interpret as a time since volatile loss in the host magma. Xenoliths erupted in pyroclastic flows appear to have retained their mantle isotope ratios, likely as a result of little prior degassing in these explosive events. High δ⁷Li, coupled with high [Li], in rapidly cooled arc peridotites may indicate that these samples represent fragments of mantle wedge that has been metasomatised by heavy, slab-derived fluids. If such material is typically stirred back into the convecting mantle, it may account for the heavy δ⁷Li seen in some oceanic basalts.     

Stable isotope and petrologic evidence for open-system degassing during the climactic and pre-climactic eruptions of Mt. Mazama, Crater Lake, Oregon

Evaluation of the extent of volatile element recycling in convergent margin volcanism requires delineating likely source(s) of magmatic volatiles through stable isotopic characterization of sulfur, hydrogen and oxygen in erupted tephra with appropriate assessment of modification by degassing. The climactic eruption of Mt. Mazama ejected approximately 50 km³ of rhyodacitic magma into the atmosphere and resulted in formation of a 10-km diameter caldera now occupied by Crater Lake, Oregon (lat. 43°N, long. 122°W). Isotopic compositions of whole-rocks, matrix glasses and minerals from Mt. Mazama climactic, pre-climactic and postcaldera tephra were determined to identify the likely source(s) of H₂O and S. Integration of stable isotopic data with petrologic data from melt inclusions has allowed for estimation of pre-eruptive dissolved volatile concentrations and placed constraints on the extent, conditions and style of degassing.Sulfur isotope analyses of climactic rhyodacitic whole rocks yield δ³⁴S values of 2.8-14.8‰ with corresponding matrix glass values of 2.4-13.2‰. δ³⁴S tends to increase with stratigraphic height through climactic eruptive units, consistent with open-system degassing. Dissolved sulfur concentrations in melt inclusions (MIs) from pre-climactic and climactic rhyodacitic pumices varies from 80 to 330 ppm, with highest concentrations in inclusions with 4.8-5.2 wt% H₂O (by FTIR). Up to 50% of the initial S may have been lost through pre-eruptive degassing at depths of 4-5 km. Ion microprobe analyses of pyrrhotite in climactic rhyodacitic tephra and andesitic scoria indicate a range in δ³⁴S from −0.4‰ to 5.8‰ and from −0.1‰ to 3.5‰, respectively. Initial δ³⁴S values of rhyodacitic and andesitic magmas were likely near the mantle value of 0‰. Hydrogen isotope (δD) and total H₂O analyses of rhyodacitic obsidian (and vitrophyre) from the climactic fall deposit yielded values οf −103 to −53‰ and 0.23-1.74 wt%, respectively. Values of δD and wt% H₂O of obsidian decrease towards the top of the fall deposit. Samples with depleted δD, and mantle δ¹⁸O values, have elevated δ³⁴S values consistent with open-system degassing. These results imply that more mantle-derived sulfur is degassed to the Earth’s atmosphere/hydrosphere through convergent margin volcanism than previously attributed. Magmatic degassing can modify initial isotopic compositions of sulfur by >14‰ (to δ³⁴S values of 14‰ or more here) and hydrogen isotopic compositions by 90‰ (to δD values of −127‰ in this case).     

The oxygen-isotope composition of chondrules and isolated forsterite and olivine grains from the Tagish Lake carbonaceous chondrite

The oxygen-isotope compositions (obtained by laser fluorination) of hand-picked separates of isolated forsterite, isolated olivine and chondrules from the Tagish Lake carbonaceous chondrite describe a line (δ¹⁷O = 0.95 * δ¹⁸O − 3.24; R² = 0.99) similar to the trend known for chondrules from other carbonaceous chondrites. The isolated forsterite grains (Fo_99.6-99.8; δ¹⁸O = −7.2‰ to −5.5‰; δ¹⁷O = −9.6‰ to −8.2‰) are more ¹⁶O-rich than the isolated olivine grains (Fo_39.6-86.8; δ¹⁸O = 3.1‰ to 5.1‰; δ¹⁷O = −0.3‰ to 2.2‰), and have chemical and isotopic characteristics typical of refractory forsterite. Chondrules contain olivine (Fo_97.2-99.8) with oxygen-isotope compositions (δ¹⁸O = −5.2‰ to 5.9‰; δ¹⁷O = −8.1‰ to 1.2‰) that overlap those of isolated forsterite and isolated olivine. An inverse relationship exists between the Δ¹⁷O values and Fo contents of Tagish Lake isolated forsterite and chondrules; the chondrules likely underwent greater exchange with ¹⁶O-poor nebular gases than the forsterite. The oxygen-isotope compositions of the isolated olivine grains describe a trend with a steeper slope (1.1 ± 0.1, R² = 0.94) than the carbonaceous chondrite anhydrous mineral line (CCAM_slope = 0.95). The isolated olivine may have crystallized from an evolving melt that exchanged with ¹⁶O-poor gases of somewhat different composition than those which affected the chondrules and isolated forsterite. The primordial components of the Tagish Lake meteorite formed under conditions similar to other carbonaceous chondrite meteorite groups, especially CMs. Its alteration history has its closest affinities to CI carbonaceous chondrites.     

Stable sulfur and oxygen isotope ratios of the Świętokrzyski National Park spring waters generated by natural and anthropogenic factors (south-central Poland)

This paper presents the results of an isotopic study of spring waters in ?wi?tokrzyski (Holy Cross Mountain) National Park (?NP), south-central Poland. The δ³⁴S_V-CDT and δ¹⁸O_V-SMOW of soluble sulfates (n = 40) varied from 0.5‰ to 18.1‰ and from 3.5‰ to 12.2‰, respectively. The average δ³⁴S values are closely similar to those of rainwater, soils and rocks (comprising scattered pyrite). This suggests that soluble sulfates in the springs originated from mixing of recent and historic deposition, sulfates derived from pyrite oxidation, and CS-mineralization in soils and debris. An additional anthropogenic sulfur input (inorganic fertilizer) occurs in the water of spring S-61 located in the ?wi?tokrzyski National Park buffer zone. The δ¹⁸O_V-SMOW of spring waters (n = 4) were in the range of −10.6‰ to −10.2‰ indicating that they are derived from vadose groundwater in ?NP. This was the first isotope study of spring waters in the national parks of Poland. It enabled the determination of sulfur pathways and discrimination between natural and anthropogenic sources of this element in a relatively pristine area.     

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.     

Sulfur isotope fractionation during the May 2003 eruption of Anatahan volcano, Mariana Islands: Implications for sulfur sources and plume processes

Sulfur isotope compositions of pumice and adsorbed volatiles on ash from the first historical eruption of Anatahan volcano (Mariana arc) are presented in order to constrain the sources of sulfur erupted during the period 10-21 May, 2003. The isotopic composition of S extracted from erupted pumice has a narrow range, from δ³⁴S_V-CDT +2.6‰ to +3.2‰, while the composition of sulfur adsorbed onto ash has a larger range (+2.8‰ to +5.3‰). Fractionation modeling for closed and open system scenarios suggests that degassing of SO₂ raised the δ³⁴S_V-CDT value of S dissolved in the melt from an initial composition of between +1.6‰ and +2.6‰ for closed-system degassing, or between −0.5‰ and +1.5‰ for open-system degassing, however closed-system degassing is the preferred model. The calculated values for the initial composition of the magma represent a MORB-like (δ³⁴S_V-CDT ∼ 0‰) mantle source with limited contamination by subducted seawater sulfate (δ³⁴S_V-CDT +21‰). Modeling also suggests that the δ³⁴S_V-CDT value of SO₂ gas in closed-system equilibrium with the degassed magma was between +0.9‰ and +2.5‰. The δ³⁴S_V-CDT value of sulfate adsorbed onto ash in the eruption plume (+2.8‰ to +5.1‰) is consistent with sulfate formation by oxidation of magmatic SO₂ in the eruption column. The sulfur isotope composition of sulfate adsorbed to ash changes from lower δ³⁴S values for ash erupted early in the eruption to higher δ³⁴S values for ash erupted later in the eruption. We interpret the temporal/stratigraphic change in sulfate isotopic composition to primarily reflect a change in the isotopic composition of magmatic SO₂ released from the progressively degassing magma and is attributed to the expulsion of an accumulated gas phase at the beginning of the eruption. More efficient oxidation of magmatic SO₂ gas to sulfate in the early water-rich eruption plume probably contributed to the change in S isotope compositions observed in the ash leachates.     

Copper isotope fractionation in acid mine drainage

We measured the Cu isotopic composition of primary minerals and stream water affected by acid mine drainage in a mineralized watershed (Colorado, USA). The δ⁶⁵Cu values (based on ⁶⁵Cu/⁶³Cu) of enargite (δ⁶⁵Cu = −0.01 ± 0.10‰; 2σ) and chalcopyrite (δ⁶⁵Cu = 0.16 ± 0.10‰) are within the range of reported values for terrestrial primary Cu sulfides (−1‰ < δ⁶⁵Cu < 1‰). These mineral samples show lower δ⁶⁵Cu values than stream waters (1.38‰ ? δ⁶⁵Cu ? 1.69‰). The average isotopic fractionation (Δ_aq-min = δ⁶⁵Cu_aq − δ⁶⁵Cu_min, where the latter is measured on mineral samples from the field system), equals 1.43 ± 0.14‰ and 1.60 ± 0.14‰ for chalcopyrite and enargite, respectively. To interpret this field survey, we leached chalcopyrite and enargite in batch experiments and found that, as in the field, the leachate is enriched in ⁶⁵Cu relative to chalcopyrite (1.37 ± 0.14‰) and enargite (0.98 ± 0.14‰) when microorganisms are absent. Leaching of minerals in the presence of Acidithiobacillus ferrooxidans results in smaller average fractionation in the opposite direction for chalcopyrite (Δ_aq-min^o=-0.57±0.14‰, where min^o refers to the starting mineral) and no apparent fractionation for enargite (Δ_aq-min^o=0.14±0.14‰). Abiotic fractionation is attributed to preferential oxidation of ⁶⁵Cu⁺ at the interface of the isotopically homogeneous mineral and the surface oxidized layer, followed by solubilization. When microorganisms are present, the abiotic fractionation is most likely not seen due to preferential association of ⁶⁵Cu_aq with A. ferrooxidans cells and related precipitates. In the biotic experiments, Cu was observed under TEM to occur in precipitates around bacteria and in intracellular polyphosphate granules. Thus, the values of δ⁶⁵Cu in the field and laboratory systems are presumably determined by the balance of Cu released abiotically and Cu that interacts with cells and related precipitates. Such isotopic signatures resulting from Cu sulfide dissolution should be useful for acid mine drainage remediation and ore prospecting purposes.     

Stable isotope fractionation during bacterial sulfate reduction is controlled by reoxidation of intermediates

Bacterial sulfate reduction is one of the most important respiration processes in anoxic habitats and is often assessed by analyzing the results of stable isotope fractionation. However, stable isotope fractionation is supposed to be influenced by the reduction rate and other parameters, such as temperature. We studied here the mechanistic basics of observed differences in stable isotope fractionation during bacterial sulfate reduction. Batch experiments with four sulfate-reducing strains (Desulfovibrio desulfuricans, Desulfobacca acetoxidans, Desulfonatronovibrio hydrogenovorans, and strain TRM1) were performed. These microorganisms metabolize different carbon sources (lactate, acetate, formate, and toluene) and showed broad variations in their sulfur isotope enrichment factors. We performed a series of experiments on isotope exchange of ¹⁸O between residual sulfate and ambient water. Batch experiments were conducted with ¹⁸O-enriched (δ¹⁸O_water = +700‰) and depleted water (δ¹⁸O_water = −40‰), respectively, and the stable ¹⁸O isotope shift in the residual sulfate was followed. For Desulfovibrio desulfuricans and Desulfonatronovibrio hydrogenovorans, which are both characterized by low sulfur isotope fractionation (ε_S > −13.2‰), δ¹⁸O values in the remaining sulfate increased by only 50‰ during growth when ¹⁸O-enriched water was used for the growth medium. In contrast, with Desulfobacca acetoxidans and strain TRM1 (ε_S < −22.7‰) the residual sulfate showed an increase of the sulfate δ¹⁸O close to the values of the enriched water of +700‰. In the experiments with δ¹⁸O-depleted water, the oxygen isotope values in the residual sulfate stayed fairly constant for strains Desulfovibrio desulfuricans, Desulfobacca acetoxidans and Desulfonatronovibrio hydrogenovorans. However, strain TRM1, which exhibits the lowest sulfur isotope fractionation factor (ε_S < −38.7‰) showed slightly decreasing δ¹⁸O values.Our results give strong evidence that the oxygen atoms of sulfate exchange with water during sulfate reduction. However, this neither takes place in the sulfate itself nor during formation of APS (adenosine-5′-phosphosulfate), but rather in intermediates of the sulfate reduction pathway. These may in turn be partially reoxidized to form sulfate. This reoxidation leads to an incorporation of oxygen from water into the “recycled” sulfate changing the overall ¹⁸O isotopic composition of the remaining sulfate fraction. Our study shows that such incorporation of ¹⁸O is correlated with the stable isotope enrichment factor for sulfur measured during sulfate reduction. The reoxidation of intermediates of the sulfate reduction pathway does also strongly influence the sulfur stable isotope enrichment factor. This aforesaid reoxidation is probably dependent on the metabolic conversion of the substrate and therefore also influences the stable isotope fractionation factor indirectly in a rate dependent manner. However, this effect is only indirect. The sulfur isotope enrichment factors for the kinetic reactions themselves are probably not rate dependent.     

Cl,N, C isotopic analysis of common agro-chemicals for identifying non-point source agricultural contaminants

The isotopic compositions of commercially available herbicides were analyzed to determine their respective ¹⁵N, ¹³C and ³⁷Cl signatures for the purposes of developing a discrete tool for tracing and identifying non-point source contaminants in agricultural watersheds. Findings demonstrate that of the agrochemicals evaluated, chlorine stable isotopes signatures range between δ³⁷Cl = −4.55‰ and +3.40‰, whereas most naturally occurring chlorine stable isotopes signatures, including those of road salt, sewage sludge and fertilizers, vary in a narrow range about the Standard Mean Ocean Chloride (SMOC) between −2.00‰ and +1.00‰. Nitrogen stable isotope values varied widely from δ¹⁵N = −10.86‰ to +1.44‰ and carbon stable isotope analysis gave an observed range between δ¹³C = −37.13‰ and −21.35‰ for the entire suite of agro-chemicals analyzed. When nitrogen, carbon and chlorine stable isotope analyses were compared in a cross-correlation analysis, statistically independent isotopic signatures exist suggesting a new potential tracer tool for identifying herbicides in the environment.     

Oxygen and iron isotope constraints on near-surface fractionation effects and the composition of lunar mare basalt source regions

Oxygen and iron isotope analyses of low-Ti and high-Ti mare basalts are presented to constrain their petrogenesis and to assess stable isotope variations within lunar mantle sources. An internally-consistent dataset of oxygen isotope compositions of mare basalts encompasses five types of low-Ti basalts from the Apollo 12 and 15 missions and eight types of high-Ti basalts from the Apollo 11 and 17 missions. High-precision whole-rock δ¹⁸O values (referenced to VSMOW) of low-Ti and high-Ti basalts correlate with major-element compositions (Mg#, TiO₂, Al₂O₃). The observed oxygen isotope variations within low-Ti and high-Ti basalts are consistent with crystal fractionation and match the results of mass-balance models assuming equilibrium crystallization. Whole-rock δ⁵⁶Fe values (referenced to IRMM-014) of high-Ti and low-Ti basalts range from 0.134‰ to 0.217‰ and 0.038‰ to 0.104‰, respectively. Iron isotope compositions of both low-Ti and high-Ti basalts do not correlate with indices of crystal fractionation, possibly owing to small mineral-melt iron fractionation factors anticipated under lunar reducing conditions.The δ¹⁸O and δ⁵⁶Fe values of low-Ti and the least differentiated high-Ti mare basalts are negatively correlated, which reflects their different mantle source characteristics (e.g., the presence or absence of ilmenite). The average δ⁵⁶Fe values of low-Ti basalts (0.073 ± 0.018‰, n = 8) and high-Ti basalts (0.191 ± 0.020‰, n = 7) may directly record that of their parent mantle sources. Oxygen isotope compositions of mantle sources of low-Ti and high-Ti basalts are calculated using existing models of lunar magma ocean crystallization and mixing, the estimated equilibrium mantle olivine δ¹⁸O value, and equilibrium oxygen-fractionation between olivine and other mineral phases. The differences between the calculated whole-rock δ¹⁸O values for source regions, 5.57‰ for low-Ti and 5.30‰ for high-Ti mare basalt mantle source regions, are solely a function of the assumed source mineralogy. The oxygen and iron isotope compositions of lunar upper mantle can be approximated using these mantle source values. The δ¹⁸O and δ⁵⁶Fe values of the lunar upper mantle are estimated to be 5.5 ± 0.2‰ (2σ) and 0.085 ± 0.040‰ (2σ), respectively. The oxygen isotope composition of lunar upper mantle is identical to the current estimate of Earth’s upper mantle (5.5 ± 0.2‰), and the iron isotope composition of the lunar upper mantle overlaps within uncertainty of estimates for the terrestrial upper mantle (0.044 ± 0.030‰).     

Limitations on the climatic and ecological signals provided by the δC values of phytoliths from a C4 North American prairie grass

Silica phytoliths, which are deposits of opal-A that precipitate in the intra- and intercellular spaces of plant tissues during transpiration, commonly contain small amounts of occluded organic matter. In this paper, we investigate whether the δ¹³C values of phytoliths from a C₄ grass, Calamovilfa longifolia, vary in response to climatic variables that can affect the carbon-isotope composition of plant tissues. There is no significant correlation (r² < 0.3) between climate variables and the δ¹³C values of C. longifolia tissues (average δ¹³C_tissue = −13.1 ± 0.6 ‰; n = 70) across the North American prairies. However, plant tissue δ¹³C values are lower for grasses collected in populated areas where the δ¹³C value of atmospheric CO₂ is expected to be lower because of fossil fuel burning. Phytolith δ¹³C values are more variable (δ¹³C = −27.3 to −23.0‰; average = −25.1 ± 1.3‰; n = 34) and more sensitive to changes in aridity than whole tissue δ¹³C values. The strongest correlations are obtained between the δ¹³C values of stem or sheath phytoliths and humidity (r² = 0.3), latitude (r² = 0.4) and amount of precipitation (r² = 0.5). However, use of these relationships is limited by the wide spread in δ¹³C values of phytoliths from different plant tissues at the same location. We have been unable to infer any relationship between δ¹³C values of phytoliths and expected variations in the δ¹³C values of atmospheric CO₂. The C. longifolia phytoliths are depleted of ¹³C relative to tissue carbon by 10-14‰. This means that the phytoliths examined in this study have carbon isotopic compositions within the range reported previously for phytoliths from C₃ plants. This observation may further limit the usefulness of soil-phytolith assemblage δ¹³C values for identifying shifts in grassland C₃:C₄ ratios.     

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.     

Note on the isotopic geochemistry of fossil-lacustrine tufas in carbonate plateau—A study from Dungul region (SW Egypt)

Lithological, chemical, and stable isotope data are used to characterize lacustrine tufas dating back to pre-late Miocene and later unknown times, capping different surfaces of a Tertiary carbonate (Sinn el-Kedab) plateau in Dungul region in the currently hyperarid south-western Egypt. These deposits are composed mostly of calcium carbonate, some magnesium carbonate and clastic particles plus minor amounts of organic matter. They have a wide range of (Mg/Ca)_molar ratios, from 0.03 to 0.3. The bulk-tufa carbonate has characteristic isotope compositions: (δ¹³C_mean = −2.49 ± 0.99‰; δ¹⁸O_mean = −9.43 ± 1.40‰). The δ¹³C values are consistent with a small input from C4 vegetation or thinner soils in the recharge area of the tufa-depositing systems. The δ¹⁸O values are typical of fresh water carbonates. Covariation between δ¹³C and δ¹⁸O values probably is a reflection of climatic conditions such as aridity. The tufas studied are isotopically similar to the underlying diagenetic marine chalks, marls and limestones (δ¹³C_mean = −2.06 ± 0.84‰; δ¹⁸O_mean = −10.06 ± 1.39‰). The similarity has been attributed to common meteoric water signatures. This raises large uncertainties in using tufas (Mg/Ca)_molar, δ¹³C and δ¹⁸O records as proxies of paleoclimatic change and suggests that intrinsic compositional differences in material sources within the plateau may mask climatic changes in the records.     

Sulfur isotope signatures in gypsiferous sediments of the Estancia and Tularosa Basins as indicators of sulfate sources, hydrological processes, and microbial activity

In order to reconstruct paleo-environmental conditions for the saline playa lakes of the Rio Grande Rift, we investigated sediment sulfate sources using sulfur isotope compositions of dissolved ions in modern surface water, groundwater, and precipitated in the form of gypsum sediments deposited during the Pleistocene and Holocene in the Tularosa and Estancia Basins. The major sulfate sources are Lower and Middle Permian marine evaporites (δ³⁴S of 10.9-14.4‰), but the diverse physiography of the Tularosa Basin led to a complex drainage system which contributed sulfates from various sources depending on the climate at the time of sedimentation. As inferred from sulfur isotope mass balance constraints, weathering of sulfides of magmatic/hydrothermal and sedimentary origin associated with climate oscillations during Last Glacial Maximum contributed about 35-50% of the sulfates and led to deposition of gypsum with δ³⁴S values of −1.2‰ to 2.2‰ which are substantially lower than Permian evaporates. In the Estancia Basin, microbial sulfate reduction appears to overprint sulfur isotopic signatures that might elucidate past groundwater flows. A Rayleigh distillation model indicates that about 3-18% of sulfates from an inorganic groundwater pool (δ³⁴S of 12.6-13.8‰) have been metabolized by bacteria and preserved as partially to fully reduced sulfur-bearing minerals species (elemental sulfur, monosulfides, disulfides) with distinctly negative δ³⁴S values (−42.3‰ to −20.3‰) compared to co-existing gypsum (−3.8‰ to 22.4‰). For the Tularosa Basin microbial sulfate reduction had negligible effect on δ³⁴S value of the gypsiferous sediments most likely because of higher annual temperatures (15-33 °C) and lower organic carbon content (median 0.09%) in those sediments leading to more efficient oxidation of H₂S and/or smaller rates of sulfate reduction compared to the saline playas of the Estancia Basin (5-28 °C; median 0.46% of organic carbon).The White Sands region of the Tularosa Basin is frequently posited as a hydrothermal analogue for Mars. High temperatures of groundwater (33.3 °C) and high δ¹⁸O(H₂O) values (1.1‰) in White Sands, however, are controlled predominantly by seasonal evaporation rather than the modern influx of hydrothermal fluids. Nevertheless, it is possible that some of the geochemical processes in White Sands, such as sulfide weathering during climate oscillations and upwelling of highly mineralized waters, might be considered as valid terrestrial analogues for the sulfate cycle in places such as Meridiani Planum on Mars.     

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.     

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.     

The record of temperature, wind velocity and air humidity in the δD and δO of water inclusions in synthetic and Messinian halites

Deuterium and oxygen isotope fractionations between liquid and vapor water were experimentally-determined during evaporation of a NaCl solution (35 g L⁻¹) as a function of water temperature and wind velocity. In the case of a null wind velocity, slopes of δD-δ¹⁸O trajectories of residual waters hyperbolically decrease with increasing water temperatures in the range 23-47 °C. For wind velocities ranging from 0.8 to 2.2 m s⁻¹, slopes of the δD-δ¹⁸O trajectories linearly increase with increasing wind velocity at a given water temperature. These experimental results can be modeled by using Rayleigh distillation equations taking into account wind-related kinetics effects. Deuterium and oxygen isotope compositions of water inclusions trapped by the precipitated halite crystals were determined by micro-equilibration techniques.These isotopic compositions accurately reflect those of the surrounding residual waters during halite growth. Isotopic compositions of water inclusions in twenty natural halites from the Messinian Realmonte mine in Sicily suggest precipitation temperatures of that match the homogenization temperatures obtained by microthermometry (median = 34 ± 5 °C). The similarity between the measured and experimental slopes of the δD-δ¹⁸O evaporation trajectories suggests that the effect of wind was negligible during the genesis of these halite deposits. Hydrogen and oxygen isotope compositions of water inclusions from Realmonte halite also define a linear trend whose extrapolation until intersection with the Mediterranean Meteoric Water Line allows the characterization of the water source with δD and δ¹⁸O values of −70 ± 10‰ and −11.5 ± 1.5‰, respectively. These results reveal that the huge amounts of salts deposited in Sicily result from the evaporation of seawater mixed with a dominant fraction (?50%) of meteoric waters most likely deriving from alpine fluvial discharge.     

Sulfur and carbon isotope records from 1700 to 800 Ma carbonates of the Jixian section, northern China: Implications for secular isotope variations in Proterozoic seawater and relationships to global supercontinental events

This study is a comprehensive, stable isotope survey of the marine carbonate-dominated, upper Paleo- to lower Neoproterozoic stratigraphy of Jixian County, China. Carbonate-associated sulfate (CAS) was extracted and measured for δ³⁴S_CAS using the same samples analyzed for δ¹³C_carbonate. This integrated proxy approach is a step towards a more comprehensive picture of secular variation in the composition of Proterozoic seawater. We specifically sampled marine carbonate intervals from the lower section of the Chuanlinggou Formation, Changcheng Group (ca. 1700 Ma) to the top of the Jingeryu Formation, Qingbaikou Group (ca. 800 Ma). δ¹³C_carbonate values are mostly negative in the upper Paleoproterozoic Changcheng Group, with an ascending trend from −3‰ to 0‰. We observed variation of approximately 0 ± 1‰ in the Mesoproterozoic Jixian Group, and positive values of +2 ± 2‰ characterize the lower Neoproterozoic Qingbaikou Group. Stratigraphic variations in δ³⁴S_CAS are more remarkable in their ranges and magnitudes, including conspicuously high values exceeding +30‰ in the three intervals at ca. 1700 Ma, 1300-1100 Ma, and 1000-900 Ma. In the Changcheng Group, δ³⁴S_CAS values are typically higher than +25‰, with only a few values of less than +15‰. In contrast, most of the data spanning from the Mesoproterozoic Tieling Formation of the Jixian Group to the lower Neoproterozoic Jingeryu Formation of the Qingbaikou Group are highly variable between +10‰ and +25‰, with some values exceeding +25‰.In the late Paleoproterozoic (1700-1600 Ma), a >10‰ decrease in δ³⁴S_CAS and ∼3‰ increase in δ¹³C_carbonate are coincident with, and likely related to, the breakup of Columbia, a supercontinent that predated Rodinia. Carbon and sulfur isotope data from the Mesoproterozoic, when global tectonic activity was comparatively weaker, fall mostly in the ranges of +15 ± 10‰ and 0 ± 1‰, respectively, but fluctuations of >20‰ for δ³⁴S_CAS and >3‰ for the δ¹³C_carbonate at ca. 1450-1400 Ma may reflect subduction and large-scale magmatic activity in island arcs marking the end of Columbia breakup. From the late Mesoproterozoic (ca. 1300-1100 Ma) to the early Neoproterozoic (ca. 800 Ma), the δ¹³C and δ³⁴S of seawater increased gradually with increasing variability. Most impressive areδ³⁴S_CAS values that exceed +30‰ in two intervals at ca. 1300-1100 Ma and ca. 1000-900 Ma, which may reflect the assembly and early breakup of Rodinia. Although gaps in the record remain, and studies of even higher resolution are warranted, our results suggest that changes in paleoceanographic conditions linked to global tectonics strongly influenced the biogeochemical cycles of C and S. Furthermore, periods of the Proterozoic previously noted for their isotopic invariability show clear isotopic expressions of this tectonic activity.