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.     

Application of sulphur isotope ratios to examine weaning patterns and freshwater fish consumption in Roman Oxfordshire, UK

This study investigates the application of sulphur isotope ratios (δ³⁴S) in combination with carbon (δ¹³C) and nitrogen (δ¹⁵N) ratios to understand the influence of environmental sulphur on the isotopic composition of archaeological human and faunal remains from Roman era sites in Oxfordshire, UK. Humans (n = 83), terrestrial animals (n = 11), and freshwater fish (n = 5) were analysed for their isotope values from four locations in the Thames River Valley, and a broad range of δ³⁴S values were found. The δ³⁴S values from the terrestrial animals were highly variable (−13.6‰ to +0.5‰), but the δ³⁴S values of the fish were clustered and ³⁴S-depleted (−20.9‰ to −17.3‰). The results of the faunal remains suggest that riverine sulphur influenced the terrestrial sulphur isotopic signatures. Terrestrial animals were possibly raised on the floodplains of the River Thames, where highly ³⁴S-depleted sulphur influenced the soil. The humans show the largest range of δ³⁴S values (−18.8‰ to +9.6‰) from any archaeological context to date. No differences in δ³⁴S values were found between the males (−7.8 ± 6.0‰) and females (−5.3 ± 6.8‰), but the females had a linear correlation (R² = 0.71; p < 0.0001) between their δ¹⁵N and δ³⁴S compositions. These δ³⁴S results suggest a greater dietary variability for the inhabitants of Roman Oxfordshire than previously thought, with some individuals eating solely terrestrial protein resources and others showing a diet almost exclusively based on freshwater protein such as fish. Such large dietary variability was not visible by analysing only the carbon and nitrogen isotope ratios, and this research represents the largest and most detailed application of δ³⁴S analysis to examine dietary practices (including breastfeeding and weaning patterns) during the Romano-British Period.     

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.     

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.     

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.     

Sulfur isotopic studies of continental flood basalts in the Noril’sk region: implications for the association between lavas and ore-bearing intrusions

Previous studies of both ore and non-ore-bearing intrusives in the Permo-Triassic flood basalts of the Siberian platform in the Noril’sk area have shown that high-grade Ni-Cu-platinum group elements (PGE) mineralization is associated with anomalously high δ³⁴S values of ∼8 to 12‰. In addition, several researchers have proposed that observed depletions in the Cu, Ni, and PGE content of basaltic lavas of the Nadezhdinsky (Nd) Formation are related to diffusional exchange with, and upgrading in metal tenor of, sulfides in the volcanic conduit system. Sulfur isotopic studies of the lavas at Noril’sk were initiated to determine if interaction with crustally derived sulfur in the conduit system was evident, and if the Nd lavas in particular were characterized by an anomalous isotopic signature. δ³⁴S values of the lavas range from −4.5 to 8.7‰ Vienna Cañon Diablo Troilite (VCDT), with S concentrations from <40 to 1373 ppm. The majority of δ³⁴S values range from 0 to 4‰, and are similar to those from S-poor intrusions in the Noril’sk area. Although textural data are not supportive of early sulfide saturation and the presence of immiscible sulfide droplets in the lavas, recrystallization may have erased expected mineralogical and textural evidence. Mineralogical data indicate that hydrothermal alteration of the lavas has occurred, but S redistribution has been restricted to localized areas and δ³⁴S values have not been affected. The relatively low S concentrations of the lavas are thought to be due in large part to degassing of the lavas in the shallow conduit system and during eruption. Our calculations are consistent with the premise that degassing of basaltic magmas at temperatures in excess of ∼900°C at QFM leads to only minor ³⁴S-depletion of sulfur remaining in the melt, and decreases in δ³⁴S values of less than 2‰ at 90% degassing. For this reason all lavas with δ³⁴S values in excess of ∼ 2‰ require a contribution of ³⁴S-enriched country rock sulfur. Because of the high S content and δ³⁴S value (∼ 16-20‰) of evaporites in the country rocks at Noril’sk, contamination of less than 0.5% is required to explain the most ³⁴S-enriched lavas. The Nd lavas have an average δ³⁴S of 2.9‰, but show no difference in S isotopic composition relative to the other lavas, suggesting that metal depletion involved only limited S transfer, or that exchange between mantle-derived S and S of crustal origin buffered δ³⁴S values to less than ∼5‰. Anomalously positive δ³⁴S values, similar to those of the ore-bearing intrusives in the Noril’sk region, are not consistently found in low-S rocks, either lavas or intrusives. Although the mechanism for the derivation of sulfide in the ore-bearing intrusions remain speculative, it is clear that the formation of sulfide ores characterized by high metal tenors proceeded only in the presence of sulfur of crustal origin.     

Oxygen and carbon isotopic systematics of aragonite speleothems and water in Furong Cave, Chongqing, China

To understand oxygen and carbon stable isotopic characteristics of aragonite stalagmites and evaluate their applicability to paleoclimate, the isotopic compositions of active and fossil aragonite speleothems and water samples from an in situ multi-year (October 2005-July 2010) monitoring program in Furong Cave located in Chongqing of China have been examined. The observations during October 2005-June 2007 show that the meteoric water is well mixed in the overlying 300-500-m bedrock aquifer, reflected by relatively constant δ¹⁸O, ±0.11-0.14‰ (1σ), of drip waters in the cave, which represents the annual status of rainfall water. Active cave aragonite speleothems are at oxygen isotopic equilibrium with drip water and their δ¹⁸O values capture the surface-water oxygen isotopic signal. Aragonite-to-calcite transformation since the last glaciation is not noticeable in Furong stalagmites. Our multi-year field experiment approves that aragonite stalagmite δ¹⁸O records in this cave are suitable for paleoclimate reconstruction. With high U, 0.5-7.2 ppm, and low Th, 20-1270 ppt, the Furong aragonite stalagmites provide very precise chronology (as good as ±20s yrs (2σ)) of the climatic variations since the last deglaciation. The synchroneity of Chinese stalagmite δ¹⁸O records at the transition into the Bølling-Allerød (t-BA) and the Younger Dryas from Furong, Hulu and Dongge Caves supports the fidelity of the reconstructed East Asian monsoon evolution. However, the Furong record shows that the cold Older Dryas (OD) occurred at 14.0 thousand years ago, agreeing with Greenland ice core δ¹⁸O records but ∼200 yrs younger than that in the Hulu record. The OD age discrepancy between Chinese caves can be attributable to different regionally climatic/environmental conditions or chronological uncertainty of stalagmite proxy records, which is limited by changes in growth rate and subsampling intervals in absolute dating. Seasonal dissolved inorganic carbon δ¹³C variations of 2-3‰ in the drip water and 5-7‰ in the pool and spring waters are likely attributed to variable degrees of CO₂ degassing in winter and summer. The variable δ¹³C values of active deposits from −11‰ to 0‰ could be caused by kinetically mediated CO₂ degassing processes. The complicated nature of pre-deposition kinetic isotopic fractionation processes for carbon isotopes in speleothems at Furong Cave require further study before they can be interpreted in a paleoclimatic or paleoenvironmental context.     

Sulfur isotope systematics of the 1982 El Chichón trachyandesite: an ion microprobe study

Sulfur isotopic compositions were determined by ion microprobe for 36 spots on anhydrite crystals in trachyandesitic pumices erupted from El Chichón Volcano in 1982. Individual anhydrite crystals are homogeneous in δ³⁴S, within the ±1‰ (2σ) uncertainty of the method, but crystal-to-crystal variations are large (+2.5 to +10.9‰). The mean δ³⁴S for anhydrite (+6.4 ± 2.1‰, 1σ) is significantly lower than earlier results for bulk anhydrite separates (+9.0 to +9.2‰). The difference between the mean δ³⁴S values in these two populations may reflect a grain-size effect, with heavier sulfur concentrated in smaller anhydrite crystals, few of which were analyzed by ion microprobe. Variations in δ³⁴S show no correlation with complex textures in anhydrite revealed by cathodoluminescence color. Ion-microprobe analyses of δ³⁴S were also obtained on six ovoid-shaped inclusions of pyrrhotite, chalcopyrite, and/or intermediate sulfide solid solution hosted by silicate or oxide crystals, interpreted to be magmatic (δ³⁴S = −0.1 to +2.7‰; mean +0.7‰), and on four irregularly shaped multiphase sulfide fragments in the matrix, interpreted as xenocrystic, which range widely in δ³⁴S (−3.7 to +5.5‰). We evaluate four different mixing scenarios involving (1) magmatic anhydrite and sedimentary sulfate, (2) magmatic anhydrite and hydrothermal anhydrite, and anhydrite and coexisting sulfide crystals precipitated in different domains of a common magma reservoir that were affected by (3) different degrees of degassing or (4) different degrees of crustal sulfur contamination. The model involving physical contamination of sedimentary sulfate is considered untenable. The other three models are considered to be viable, but none of them can explain all observations. The results of this study and other recent investigations prompt a re-evaluation of the sulfur budget for the 1982 El Chichón eruption. We estimate that 2.2 × 10¹³ g of S was emitted, and that 58 wt.% of the sulfur was present as anhydrite prior to eruption, with the remainder in a vapor phase, with H₂S/SO₂ ≈ 9. The bulk magmatic δ³⁴S value for the 1982 El Chichón trachyandesite is estimated as +4.1 to +5.8‰, typical of the relatively heavy sulfur isotopic compositions that characterize subduction-related magmas.     

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.     

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.     

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.     

O and S isotopic composition of dissolved and attached oxidation products of pyrite by Acidithiobacillus ferrooxidans: Comparison with abiotic oxidations

The acidophilic iron-oxidizing bacterium, Acidithiobacillus ferrooxidans, plays a part in the pyrite oxidation process and has been widely studied in order to determine the kinetics of the reactions and the isotopic composition of dissolved product sulphates, but the details of the oxidation processes at the surface of pyrite are still poorly known. In this study, oxygen and sulphur isotopic compositions (δ¹⁸O and δ³⁴S) were analyzed for dissolved sulphates and water from experimental aerobic acidic (pH < 2) pyrite oxidation by A. ferrooxidans. The oxidation products attached to the pyrite surfaces were studied for their morphology (SEM), their chemistry (Raman spectroscopy) and for their δ¹⁸O (ion microprobe). They were compared to abiotically (Fe³⁺, H₂O₂, O₂) oxidized pyrite surface compounds in order to constrain the oxidation pathways and to look for the existence of potential biosignatures for this system.The pyrite dissolution evolved from non-stoichiometric (during the first days) to stoichiometric (with increasing time) resulting in dissolved sulphates having distinct δ¹⁸O (e.g. +11.0‰ and −2.0‰, respectively) and δ³⁴S (+4.5‰ and +2.8‰, respectively) values. The “oxidation layer” at the surface of pyrite is complex and made of iron oxides, sulphate, polysulphide, elemental sulphur and polythionates. Bio- and Fe³⁺-oxidation favour the development of monophased micrometric bumps made of hematite or sulphate while other abiotic oxidation processes result in more variable oxidation products. The δ¹⁸O of these oxidation products at the surface of oxidized pyrites are strongly variable (from ≈−40‰ to ≈+30‰) for all experiments.Isotopic fractionation between sulphates and pyrite, Δ³⁴S_SO4-pyrite, is equal to −1.3‰ and +0.4‰ for sulphates formed by stoichiometric and non-stoichiometric processes, respectively. These two values likely reflect either a S-S or a Fe-S bond breaking process. The Δ¹⁸O_SO4-H2O and Δ¹⁸O_SO4-O2 are estimated to be ≈+16‰ and ≈−25‰, respectively. These values are higher than previously published data and may reflect biological effects. The large δ¹⁸O heterogeneity measured at the surfaces of oxidized pyrites, whatever the oxidant, may be related (i) to the existence of local surface environments isolated from the solution in which the oxidation processes are different and (ii) to the stabilization at the pyrite surface of reaction intermediates that are not in isotopic equilibrium with the solution. Though the oxygen isotopic composition of surface oxidation products cannot be taken as a direct biosignature, the combined morphological, chemical and isotopic characterization of the surfaces of oxidized pyrites may furnish clues about a biological activity on a mineral surface.     

Paleoclimate record in the Nubian Sandstone Aquifer,Sinai Peninsula,Egypt

Sixteen groundwater samples collected from production wells tapping Lower Cretaceous Nubian Sandstone and fractured basement aquifers in Sinai were analyzed for their stable isotopic compositions, dissolved noble gas concentrations (recharge temperatures), tritium activities, and ¹⁴C abundances. Results define two groups of samples: Group I has older ages, lower recharge temperatures, and depleted isotopic compositions (adjusted ¹⁴C model age: 24,000–31,000 yr BP; δ¹⁸O: − 9.59‰ to − 6.53‰; δ²H: − 72.9‰ to − 42.9‰; < 1 TU; and recharge T: 17.5–22.0°C) compared to Group II (adjusted ¹⁴C model age: 700–4700 yr BP; δ¹⁸O: − 5.89‰ to − 4.84‰; δ²H: − 34.5‰ to − 24.1‰; < 1 to 2.78 TU; and recharge T: 20.6–26.2°C). Group II samples have isotopic compositions similar to those of average modern rainfall, with larger d-excess values than Group I waters, and locally measurable tritium activity (up to 2.8 TU). These observations are consistent with (1) the Nubian Aquifer being largely recharged prior to and/or during the Last Glacial Maximum (represented by Group I), possibly through the intensification of paleowesterlies; and (2) continued sporadic recharge during the relatively dry and warmer interglacial period (represented by Group II) under conditions similar to those of the present.     

Mineralogical and isotopic record of biotic and abiotic diagenesis of the Callovian-Oxfordian clayey formation of Bure (France)

The Callovian-Oxfordian (COx) clayey unit is being studied in the Eastern part of the Paris Basin at depths between 400 and 500 m depth to assess of its suitability for nuclear waste disposal. The present study combines new mineralogical and isotopic data to describe the sedimentary history of the COx unit. Petrologic study provided evidence of the following diagenetic mineral sequence: (1) framboidal pyrite and micritic calcite, (2) iron-rich euhedral carbonates (ankerite, sideroplesite) and glauconite (3) limpid calcite and dolomite and celestite infilling residual porosity in bioclasts and cracks, (4) chalcedony, (5) quartz/calcite. Pyrite in bioturbations shows a wide range of δ³⁴S (−38‰ to +34.5‰), providing evidence of bacterial sulphate reduction processes in changing sedimentation conditions. The most negative values (−38‰ to −22‰), measured in the lower part of the COx unit indicate precipitation of pyrite in a marine environment with a continuous sulphate supply. The most positive pyrite δ³⁴S values (−14‰ up to +34.5‰) in the upper part of the COx unit indicate pyrite precipitation in a closed system. Celestite δ³⁴S values reflect the last evolutionary stage of the system when bacterial activity ended; however its deposition cannot be possible without sulphate supply due to carbonate bioclast dissolution. The ⁸⁷Sr/⁸⁶Sr ratio of celestite (0.706872-0.707040) is consistent with deposition from Jurassic marine-derived waters. Carbon and oxygen isotopic compositions of bulk calcite and dolomite are consistent with marine carbonates. Siderite, only present in the maximum clay zone, has chemical composition and δ¹⁸O consistent with a marine environment. Its δ¹³C is however lower than those of marine carbonates, suggesting a contribution of ¹³C-depleted carbon from degradation of organic matter. δ¹⁸O values of diagenetic chalcedony range between +27‰ and +31‰, suggesting precipitation from marine-derived pore waters. Late calcite crosscutting a vein filled with chalcedony and celestite, and late euhedral quartz in a limestone from the top of the formation have lower δ¹⁸O values (∼+19‰), suggesting that they precipitated from meteoric fluids, isotopically close to present-day pore waters of the formation. Finally, the study illustrates the transition from very active, biotic diagenesis to abiotic diagenesis. This transition appears to be driven by compaction of the sediment, which inhibited movement of bacterial cells by reduction of porosity and pore sizes, rather than a lack of inorganic carbon or sulphates.     

Environmental memory and a possible seasonal bias in the stable isotope composition of (U-Th)/He-dated goethite from the Canadian Arctic

Goethite (Ax-2) from Axel Heiberg Island (∼80°N) on the margin of the Arctic Ocean is the dominant mineral in a sample of “petrified” Eocene wood, but U, Th, and He measurements suggest that the goethite (α-FeOOH) crystallized in the latest Miocene/Pliocene (ca. 5.5 to 2.8 Ma). Measured δD and δ¹⁸O values of Ax-2 are −221 (±6)‰ and −9.6 (±0.5)‰, respectively. The inferred δD and δ¹⁸O values of the ancient water were about −139‰ and −18.6‰, respectively, with a calculated temperature of crystallization of 3 (±5)°C, which compares with the modern summer (J-J-A) temperature of 3 °C and contrasts with a modern MAT of −19 °C. Published results from various biological proxies on nearby Ellesmere Island indicate a Pliocene (∼4 Ma) MAT of either −6 or −0.4 °C and corresponding seasonal amplitudes of about 18 or 13 °C. A conductive heat flow model suggests that a temperature of 3 °C could represent goethite crystallization at depths of ∼100-200 cm (for MAT = −6 °C) or ∼250-450 cm (for MAT = −0.4 °C) over seasonally restricted intervals of time.The δ¹⁸O value of the Ax-2 water (−18.6‰) is more positive than the modern J-J-A precipitation (−22‰). In combination, the paleotemperatures and δ¹⁸O values of ancient waters (from Ax-2 and published results from three Eocene or Pliocene proxy sites on Axel Heiberg and Ellesmere Islands) are consistent with a warm season bias in those isotopic proxies. The results are also consistent with higher proportions of J-J-A precipitation in the annual total. If so, this emphasizes the importance of seasonality at high latitudes even in times of warmer global climates, and suggests that the Arctic hydrologic cycle, as expressed in the seasonal distribution and isotopic composition of precipitation (perhaps modified by a warmer Arctic Ocean), differed from modern.The δ¹³C value of the Fe(CO₃)OH component in the Ax-2 goethite is +6.6‰, which is much more positive than expected if crystallizing goethite incorporated CO₂ derived primarily from oxidation of relict Eocene wood with δ¹³C values of about −24‰. This apparent paradox may be resolved if the goethite is a product of oxidation of ¹³C-rich siderite, which had previously replaced wood in an Eocene methanogenic burial environment. Thus, the goethite retains a carbon isotope “memory” of a diagenetic Eocene event, but a δD and δ¹⁸O record of the latest Miocene/Pliocene Arctic climate.     

Stable isotope geochemistry of ground and surface waters associated with undisturbed massive sulfide deposits; constraints on origin of waters and water-rock reactions

Stable isotopes (H, O, C) were determined for ground and surface waters collected from two relatively undisturbed massive sulfide deposits (Halfmile Lake and Restigouche) in the Bathurst Mining Camp (BMC), New Brunswick, Canada. Additional waters from active and inactive mines in the BMC were also collected. Oxygen and hydrogen isotopes of surface and shallow groundwaters from both the Halfmile Lake and Restigouche deposits are remarkably uniform (− 13 to − 14‰ and − 85 to − 95‰ for δ¹⁸O_VSMOW and δ²H_VSMOW, respectively). These values are lighter than predicted for northern New Brunswick and, combined with elevated deuterium excess values, suggest that recharge waters are dominated by winter precipitation, recharged during spring melting. Deeper groundwaters from the Restigouche deposit, and from active and inactive mines have heavier δ¹⁸O_VSMOW ratios (up to − 10.8‰) than shallow groundwaters suggesting recharge under warmer climate or mixing with Shield-type brines. Some of the co-variation in Cl concentrations and δ¹⁸O_VSMOW ratios can be explained by mixing between saline and shallow recharge water end-members. Carbon isotopic compositions of dissolved inorganic carbon (DIC) are variable, ranging from − 15 to − 5‰ δ¹³C_VPDB for most ground and surface waters. Much of the variation in the carbon isotopes is consistent with closed system groundwater evolution involving soil zone CO₂ and fracture zone carbonate minerals (calcite, dolomite and siderite; average = − 6.5‰ δ¹³C_VPDB). The DIC of saline Restigouche deposit groundwater is isotopically heavy (∼+ 12‰ δ¹³C_VPDB), indicating carbon isotopic fractionation from methanogenesis via CO₂ reduction, consistent with the lack of dissolved sulfate in these waters and the observation of CH₄-degassing during sampling.     

Micro-scale S isotope studies of the Kharaelakh intrusion, Noril’sk region, Siberia: Constraints on the genesis of coexisting anhydrite and sulfide minerals

The coexistence of magmatic anhydrite and sulfide minerals in non-arc-related mafic magmas has only rarely been documented. Likewise the S isotope fractionation between sulfate and sulfide in mafic rocks has infrequently been measured. In the Kharaelakh intrusion associated with the world-famous Noril’sk ore district in Siberia coexisting magmatic anhydrite and sulfide minerals have been identified. Sulfur isotope compositions of the anhydrite-sulfide assemblages have been measured via both ion microprobe and conventional analyses to help elucidate the origin of the anhydrite-sulfide pairs. Magmatic anhydrite and chalcopyrite are characterized by δ³⁴S values between 18.8‰ and 22.8‰, and 9.3‰ and 13.2‰, respectfully. Coexisting anhydrite and chalcopyrite show Δ values that fall between 8.5‰ and 11.9‰. Anhydrite in the Kharaelakh intrusion is most readily explained by the assimilation of sulfate from country rocks; partial reduction to sulfide led to mixing between sulfate-derived sulfide and sulfide of mantle origin. The variable anhydrite and sulfide δ³⁴S values are a function of differing degrees of sulfate reduction, variable mixing of sulfate-derived and mantle sulfide, incomplete isotopic homogenization of the magma, and a lack of uniform attainment of isotopic equilibrium during subsolidus cooling. The δ³⁴S values of sulfide minerals have changed much less with cooling than have anhydrite values due in large part to the high sulfide/sulfate ratio. Variations in both sulfide and anhydrite δ³⁴S values indicate that isotopically distinct domains existed on a centimeter scale. Late stage hydrothermal anhydrite and pyrite also occur associated with Ca-rich hydrous alteration assemblages (e.g., thomsonite, prehnite, pectolite, epidote, xonotlite). δ³⁴S values of secondary hydrothermal anhydrite and pyrite determined by conventional analyses are in the same range as those of the magmatic minerals. Anhydrite-pyrite Δ values are in the 9.1-10.1‰ range, and are smaller than anticipated for the low temperatures indicated by the silicate alteration assemblages. The small Δ values are suggestive of either sulfate-sulfide isotopic disequilibrium or closure of the system to further exchange between ∼550 and 600 °C. Our results confirm the importance of the assimilation of externally derived sulfur in the generation of the elevated δ³⁴S values in the Kharaelakh intrusion, but highlight the sulfur isotopic variability that may occur in magmatic systems. In addition, our results confirm the need for more precise experimental determination of sulfate-sulfide sulfur isotope fractionation factors in high-T systems.     

Fractionation of multiple sulfur isotopes during phototrophic oxidation of sulfide and elemental sulfur by a green sulfur bacterium

We present multiple sulfur isotope measurements of sulfur compounds associated with the oxidation of H₂S and S⁰ by the anoxygenic phototrophic S-oxidizing bacterium Chlorobium tepidum. Discrimination between ³⁴S and ³²S was +1.8 ± 0.5‰ during the oxidation of H₂S to S⁰, and −1.9 ± 0.8‰ during the oxidation of S⁰ to , consistent with previous studies. The accompanying Δ³³S and Δ³⁶S values of sulfide, elemental sulfur, and sulfate formed during these experiments were very small, less than 0.1‰ for Δ³³S and 0.9‰ for Δ³⁶S, supporting mass conservation principles. Examination of these isotope effects within a framework of the metabolic pathways for S oxidation suggests that the observed effects are due to the flow of sulfur through the metabolisms, rather than abiotic equilibrium isotope exchange alone, as previously suggested. The metabolic network comparison also indicates that these metabolisms work to express some isotope effects (between sulfide, polysulfides, and elemental sulfur in the periplasm) and suppress others (kinetic isotope effects related to pathways for oxidation of sulfide to sulfate via the same enzymes involved in sulfate reduction acting in reverse). Additionally, utilizing fractionation factors for phototrophic S oxidation calculated from our experiments and for other oxidation processes calculated from the literature (chemotrophic and inorganic S oxidation), we constructed a set of ecosystem-scale sulfur isotope box models to examine the isotopic consequences of including sulfide oxidation pathways in a model system. These models demonstrate how the small δ³⁴S effects associated with S oxidation combined with large δ³⁴S effects associated with sulfate reduction (by SRP) and sulfur disproportionation (by SDP) can produce large (and measurable) effects in the Δ³³S of sulfur reservoirs. Specifically, redistribution of material along the pathways for sulfide oxidation diminishes the net isotope effect of SRP and SDP, and can mask the isotopic signal for sulfur disproportionation if significant recycling of S intermediates occurs. We show that the different sulfide oxidation processes produce different isotopic fields for identical proportions of oxidation, and discuss the ecological implications of these results to interpreting minor S isotope patterns in modern systems and in the geologic record.     

Constraining atmospheric oxygen and seawater sulfate concentrations during Paleoproterozoic glaciation: In situ sulfur three-isotope microanalysis of pyrite from the Turee Creek Group, Western Australia

Previous efforts to constrain the timing of Paleoproterozoic atmospheric oxygenation have documented the disappearance of large, mass-independent sulfur isotope fractionation and an increase in mass-dependent sulfur isotope fractionation associated with multiple glaciations. At least one of these glacial events is preserved in diamictites of the ∼2.4 Ga Meteorite Bore Member of the Kungarra Formation, Turee Creek Group, Western Australia. Outcrop exposures of this unit show the transition from the Boolgeeda Iron Formation of the upper Hamersley Group into clastic, glaciomarine sedimentary rocks of the Turee Creek Group. Here we report in situ multiple sulfur isotope and elemental abundance measurements of sedimentary pyrite at high spatial resolution, as well as the occurrence of detrital pyrite in the Meteorite Bore Member. The 15.3‰ range of Δ³³S in one sample containing detrital pyrite (−3.6‰ to 11.7‰) is larger than previously reported worldwide, and there is evidence for mass-independent sulfur isotope fractionation in authigenic pyrite throughout the section (Δ³³S from −0.8‰ to 1.0‰). The 90‰ range in δ³⁴S observed (−45.5‰ to 46.4‰) strongly suggests microbial sulfate reduction under non-sulfate limiting conditions, indicating significant oxidative weathering of sulfides on the continents. Multiple generations of pyrite are preserved, typically represented by primary cores with low δ³⁴S (<−20‰) overgrown by euhedral rims with higher δ³⁴S (4-7‰) and enrichments in As, Ni, and Co. The preservation of extremely sharp sulfur isotope gradients (30‰/<4 μm) implies limited sulfur diffusion and provides time and temperature constraints on the metamorphic history of the Meteorite Bore Member. Together, these results suggest that the Meteorite Bore Member was deposited during the final stages of the “Great Oxidation Event,” when pO₂ first became sufficiently high to permit pervasive oxidative weathering of continental sulfides, yet remained low enough to permit the production and preservation of mass-independent sulfur isotope fractionation.     

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.