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.     

Ammonium in thermal waters of Yellowstone National Park: Processes affecting speciation and isotope fractionation

Dissolved inorganic nitrogen, largely in reduced form (), has been documented in thermal waters throughout Yellowstone National Park, with concentrations ranging from a few micromolar along the Firehole River to millimolar concentrations at Washburn Hot Springs. Indirect evidence from rock nitrogen analyses and previous work on organic compounds associated with Washburn Hot Springs and the Mirror Plateau indicate multiple sources for thermal water NH₄(T), including Mesozoic marine sedimentary rocks, Eocene lacustrine deposits, and glacial deposits. A positive correlation between NH₄(T) concentration and δ¹⁸O of thermal water indicates that boiling is an important mechanism for increasing concentrations of NH₄(T) and other solutes in some areas. The isotopic composition of dissolved NH₄(T) is highly variable (δ¹⁵N = −6‰ to +30‰) and is positively correlated with pH values. In comparison to likely δ¹⁵N values of nitrogen source materials (+1‰ to +7‰), high δ¹⁵N values in hot springs with pH >5 are attributed to isotope fractionation associated with loss by volatilization. NH₄(T) in springs with low pH typically is relatively unfractionated, except for some acid springs with negative δ¹⁵N values that are attributed to condensation. NH₄(T) concentration and isotopic variations were evident spatially (between springs) and temporally (in individual springs). These variations are likely to be reflected in biomass and sediments associated with the hot springs and outflows. Elevated NH₄(T) concentrations can persist for 10s to 1000s of meters in surface waters draining hot spring areas before being completely assimilated or oxidized.     

Ammonium loss and nitrogen isotopic fractionation in biotite as a function of metamorphic grade in metapelites from western Maine, USA

Ammonium fixed in micas of metamorphic rocks is a sensitive indicator both of organic-inorganic interactions during diagenesis as well as of the devolatilization history and fluid/rock interaction during metamorphism. In this study, a collection of geochemically well-characterized biotite separates from a series of graphite-bearing Paleozoic greenschist- to upper amphibolite-facies metapelites, western Maine, USA, were analyzed for ammonium nitrogen () contents and isotopic composition (δ¹⁵N_NH4) using the HF-digestion distillation technique followed by the EA-IRMS technique. Biotite separates, sampled from 9 individual metamorphic zones, contain 3000 to 100 ppm with a wide range in δ¹⁵N from +1.6‰ to +9.1‰. Average contents in biotite show a distinct decrease from about 2750 ppm for the lowest metamorphic grade (∼500 °C) down to 218 ppm for the highest metamorphic grade (∼685 °C). Decreasing abundances in are inversely correlated in a linear fashion with increasing K⁺ in biotite as a function of metamorphic grade and are interpreted as a devolatilization effect. Despite expected increasing δ¹⁵N_NH4 values in biotite with nitrogen loss, a significant decrease from the Garnet Zones to the Staurolite Zones was found, followed by an increase to the Sillimanite Zones. This pattern for δ¹⁵N_NH4 values in biotite inversely correlates with Mg/(Mg + Fe) ratios in biotite and is discussed in the framework of isotopic fractionation due to different exchange processes between or , reflecting devolatilization history and redox conditions during metamorphism.     

Modeling the effects of bond environment on equilibrium iron isotope fractionation in ferric aquo-chloro complexes

Four or five sets of ab initio models, including Unrestricted Hartree Fock (UHF) and hybrid Density Functional Theory (DFT) are calculated for each species in a series of aqueous ferric aquo-chloro complexes: , , , FeCl₃(H₂O)₃, FeCl₃(H₂O)₂, , FeCl₅H₂O²⁻, , ) in order to determine the relative isotopic fractionation among the complexes, to compare the results of different models for the same complexes, to examine factors that influence the magnitude of the isotopic fractionation, and to compare bond-partner-driven fractionation with redox-driven fractionation.Relative to , all models show a nearly linear decrease in ⁵⁶Fe/⁵⁴Fe as the number of Cl⁻ ions per Fe³⁺ ion increases, with slopes of −0.8‰ to −1.0‰ per Cl⁻ at 20 °C. At 20 °C, 1000 ln β (β = ⁵⁶Fe/⁵⁴Fe reduced partition function ratio relative to a dissociated Fe atom) values range from 8.93‰ to 9.73‰ for , 8.04-9.12‰ for , 7.61-8.73‰ for , 7.14-8.25‰ for , and 3.09-4.41‰ for . The fractionation between and ranges from 1.5‰ to 2.6‰, depending on the model; this is comparable in magnitude to fractionation effects due to Fe³⁺/Fe²⁺ redox reactions. β values from the UHF models are consistently higher than those from the hybrid DFT models.Isotopic fractionation is shown to be sensitive to differences in ligand bond stiffness (above), coordination number, bond length, and the frequency of the asymmetric Fe-X stretching vibrational mode, as predicted by previous theoretical studies. Complexes with smaller coordination numbers have higher 1000 ln β (7.46‰, 5.25‰, and 3.48‰ for , ,, respectively, from the B3LYP/6-31G(d) model). Species with the same number of chlorides but fewer waters also show the effect of coordination number on 1000 ln β: (7.46‰ vs. 7.05‰ for FeCl₃(H₂O)₂ vs. FeCl₃(H₂O)₃ and 5.25‰ vs. 4.94‰ for vs. FeCl₅H₂O²⁻ with the B3LYP/6-31G(d) model). As more Fe-Cl bonds substitute for Fe-OH₂ bonds (with a resulting decrease in β), the lengths of the Fe-Cl bonds and the Fe-O bonds increase.Preliminary modeling of shows an Fe³⁺/Fe²⁺ fractionation of 3.2‰ for the B3LYP/6-31G(d) model, in agreement with previous studies. The addition of an explicit outer hydration sphere of 12 H₂O molecules to models of improves agreement with measured vibrational frequencies and bond lengths; 1000 ln β increases by 0.8-1.0‰. An additional hydration sphere around increases 1000 ln β by only 0.1‰.Isotopic fractionations predicted for this simple system imply that ligands present in an aqueous iron environment are potentially important drivers of fractionation, and suggest that significant fractionation effects are likely in other aqueous systems containing sulfides or organic ligands. Fractionation effects due to both speciation and redox must be considered when interpreting iron isotope fractionations in the geological record.     

Unsteady diagenetic processes and sulfur biogeochemistry in tropical deltaic muds: Implications for oceanic isotope cycles and the sedimentary record

Sedimentary S cycling is usually conceptualized and interpreted within the context of steadily accreting (1-D) transport-reaction regimes. Unsteady processes, however, are common in many sedimentary systems and can result in dramatically different S reaction balances and diagenetic products than steady conditions. Globally important common examples include tropical deltaic topset and inner shelf muds such as those extending from the Amazon River ∼1600 km along the Guianas coast of South America. These deposits are characterized by episodic reworking of the surface seabed over vertical depths of ∼0.1-3 m. Reworked surface sediments act as unsteady, suboxic batch reactors, unconformably overlying relict anoxic, often methanic deposits, and have diagenetic properties largely decoupled from net accumulation of sediment. Despite well-oxygenated water and an abundant reactive organic matter supply, physical disturbance inhibits macrofauna, and benthic communities are dominated by microbial biomass across immense areas. In the surficial suboxic layer, molecular biological analyses, tracer experiments, sediment C/S/Fe compositions, and δ³⁴S, δ¹⁸O of pore water indicate close coupling of anaerobic C, S, and Fe cycles. δ¹⁸O- can increase by 2-3‰ during anaerobic recycling without net change in δ³⁴S-, demonstrating reduction coupled to complete anaerobic reoxidation to and a δ¹⁸O- reduction + reoxidation fractionation factor?12‰ (summed magnitudes). S reoxidation must be coupled to Fe-oxide reduction, contributing to high dissolved Fe²⁺ (∼1 mM) and Fe mobilization-export. The reworking of Amazon-Guianas shelf muds alone may isotopically alter δ¹⁸O- equivalent in mass to?25% of the annual riverine delivery of to the global ocean. Unsteady conditions result in preservation of unusually heavy δ³⁴S isotopic compositions of residual Cr reducible S, ranging from 0‰ to >30‰ in physically reworked deposits. In contrast, bioturbated facies adjacent to physically reworked regions accumulate isotopically light S (δ³⁴S to −20‰) in otherwise similar decomposition regimes. The isotopic patterns of both physically and biologically reworked regions can be simulated with simple diagenetic models. Heavy S isotopic signatures are largely a consequence of unsteady diffusion and progressive anaerobic burndown into underlying deposits, whereas isotopically depleted bioturbated deposits predominantly reflect biogenic diffusive scaling and isotopic distillation/diffusive pumping associated with reoxidation in burrow walls immediately adjacent to reduced zones. The S isotopic transition from unsteady physically controlled regions of the Amazon delta moving laterally into bioturbated facies mimics the transition of S isotopic patterns temporally in the geologic record during the rise of bioturbation. No special role for S disproportionation is required to explain these differences. The potential role of unsteady, suboxic diagenesis and dynamic reworking of sediments has been largely ignored in models of the evolution of surficial elemental cycling and interpretations of the geologic record.     

The origin of Zn isotope fractionation in sulfides

Isotope fractionation of Zn between aqueous sulfide, chloride, and carbonate species (Zn²⁺, Zn(HS)₂, , , ZnS(HS)⁻, ZnCl⁺, ZnCl₂, , and ZnCO₃) was investigated using ab initio methods. Only little fractionation is found between the sulfide species, whereas carbonates are up to 1‰ heavier than the parent solution. At pH > 3 and under atmospheric-like CO₂ pressures, isotope fractionation of Zn sulfides precipitated from sulfidic solutions is affected by aqueous sulfide species and the δ⁶⁶Zn of sulfides reflect these in the parent solutions. Under high P_CO2 conditions, carbonate species become abundant. In high P_CO2 conditions of hydrothermal solutions, Zn precipitated as sulfides is isotopically nearly unfractionated with respect to a low-pH parent fluid. In contrast, negative δ⁶⁶Zn down to at least −0.6‰ can be expected in sulfides precipitated from solutions with pH > 9. Zinc isotopes in sulfides and rocks therefore represent a potential indicator of mid to high pH in ancient hydrothermal fluids.     

Silicon isotopic fractionation during adsorption of aqueous monosilicic acid onto iron oxide

The quantification of silicon isotopic fractionation by biotic and abiotic processes contributes to the understanding of the Si continental cycle. In soils, light Si isotopes are selectively taken up by plants, and concentrate in secondary clay-sized minerals. Si can readily be retrieved from soil solution through the specific adsorption of monosilicic acid () by iron oxides. Here, we report on the Si-isotopic fractionation during adsorption on synthesized ferrihydrite and goethite in batch experiment series designed as function of time (0-504 h) and initial concentration (ic) of Si in solution (0.21-1.80 mM), at 20 °C, constant pH (5.5) and ionic strength (1 mM). At various contact times, the δ²⁹Si vs. NBS28 compositions were determined in selected solutions (ic = 0.64 and 1.06 mM Si) by MC-ICP-MS in dry plasma mode with external Mg doping with an average precision of ±0.08‰ (±2σ_SEM). Per oxide mass, ferrihydrite (74-86% of initial Si loading) adsorbed more Si than goethite (37-69%) after 504 h of contact over the range of initial Si concentration 0.42-1.80 mM. Measured against its initial composition (δ²⁹Si = +0.01 ± 0.04‰ (±2σ_SD)), the remaining solution was systematically enriched in ²⁹Si, reaching maximum δ²⁹Si values of +0.70 ± 0.07‰ for ferrihydrite and +0.50 ± 0.08‰ for goethite for ic 1.06 mM. The progressive ²⁹Si enrichment of the solution fitted better a Rayleigh distillation path than a steady state model. The fractionation factor ²⁹ε (±1σ_SD) was estimated at −0.54 ± 0.03‰ for ferrihydrite and −0.81 ± 0.12‰ for goethite. Our data imply that the sorption of onto synthetic iron oxides produced a distinct Si-isotopic fractionation for the two types of oxide but in the same order than that generated by Si uptake by plants and diatoms. They further suggest that the concentration of light Si isotopes in the clay fraction of soils is partly due to sorption onto secondary clay-sized iron oxides.     

Seasonal δS variations in two high elevation snow pits measured by S-S double spike thermal ionization mass spectrometry

δ³⁴S and sulfate concentrations were determined in snow pit samples using a thermal ionization mass spectrometric technique capable of 0.2‰ accuracy and requires ≈5 μg (0.16 μmol) natural S. The technique utilizes a ³³S-³⁶S double spike for instrumental mass fractionation correction, and has been applied to snow pit samples collected from the Inilchek Glacier, Kyrgyzstan and from Summit, Greenland. These δ³⁴S determinations provide the first high-resolution seasonal data for these sites, and are used to estimate seasonal sulfate sources. Deuterium (δD) and oxygen (δ¹⁸O) isotope data show that the Inilchek and Summit snow pit samples represent precipitation over ≈20 months.The δ³⁴S values for the Inilchek ranged from +2.6 ± 0.4‰ to +7.6 ± 0.4‰ on sample sizes ranging from 0.3 to 1.8 μmol S. δ³⁴S values for Greenland ranged from +3.6 ± 0.7‰ to +13.3 ± 5‰ for sample sizes ranging from 0.05 to 0.29 μmol S. The concentration ranged from 92.6 ± 0.4 to 1049 ± 4 ng/g for the Inilchek and 18 ± 9 to 93 ± 6 ng/g for the Greenland snow pit. Anthropogenic sulfate dominates throughout the sampled time interval for both sites based on mass balance considerations. Additionally, both sites exhibit a seasonal signature in both δ³⁴S and concentration. The thermal ionization mass spectrometric technique has three advantages compared to gas source isotopic methods: (1) sample size requirements of this technique are 10-fold less permitting access to the higher resolution S isotope record of low concentration snow and ice, (2) the double spike technique permits δ³⁴S and S concentration to be determined simultaneously, and (3) the double spike is an internal standard.     

Oxygen isotope fractionation between apatite-bound carbonate and water determined from controlled experiments with synthetic apatites precipitated at 10-37 °C

The oxygen isotope fractionation between the structural carbonate of inorganically precipitated hydroxyapatite (HAP) and water was determined in the range 10-37 °C. Values of 1000 ln α() are linearly correlated with inverse temperature (K) according to the following equation: 1000 ln α() = 25.19 (±0.53)·T⁻¹ − 56.47 (±1.81) (R² = 0.998). This fractionation equation has a slightly steeper slope than those already established between calcite and water ( [O’Neil et al., 1969] and [Kim and O’Neil, 1997]) even though measured fractionations are of comparable amplitude in the temperature range of these experimental studies. It is consequently observed that the oxygen isotope fractionation between apatite carbonate and phosphate increases from about 7.5‰ up to 9.1‰ with decreasing temperature from 37 °C to 10 °C. A compilation of δ¹⁸O values of both phosphate and carbonate from modern mammal teeth and bones confirms that both variables are linearly correlated, despite a significant scattering up to 3.5‰, with a slope close to 1 and an intercept corresponding to a 1000 ln α() value of 8.1‰. This apparent fractionation factor is slightly higher or close to the fractionation factor expected to be in the range 7-8‰ at the body temperature of mammals.     

An experimental determination of chlorine isotope fractionation in acid systems and applications to volcanic fumaroles

The δ³⁷Cl values of volcanic fumarole gases and bubbling springs were measured from the Central American and the Kurile arcs. Low temperature gas samples from the Central American arc have δ³⁷Cl values generally between −2 and 2‰, whereas high-temperature fumaroles (>100 °C) range from 4 to 12‰, with several outliers. This is in contrast to the high-temperature fumaroles from the Kurile island Kudryavy which have slightly positive δ³⁷Cl values, averaging 0.8‰ (±0.6, 1σ), and from our previous work on Izu and Mariana arc samples in which the δ³⁷Cl values of fumarole and ash samples are similar to each other and negative. Assuming that the source for the high-T Central American fumaroles has typical subduction δ³⁷Cl values (−2.5 to 1‰), then there must be a large Cl isotope fractionation in the near-surface fumarolic system. The most likely fractionation mechanism for the high δ³⁷Cl values is between Cl⁻_aq − HCl(g), but published theoretical fractionation for this pair is only ∼1.5‰, insufficient to explain the large range of values observed in the fumaroles. Three experiments were undertaken in order to identify a process that could cause the wide range of δ³⁷Cl values observed in the high-temperature fumaroles. Results are the following: (1) A sub-boiling equilibration experiment between aqueous chloride and HCl gas had , in agreement with the theoretical calculations. (2) Evaporation of HCl(g) from hydrochloric acid at room temperature had fractionation in the opposite sense, with a . (3) A ‘synthetic fumarole’ gave large positive fractionations up to 9‰, with ³⁷Cl strongly partitioned into the vapor phase. The ‘fumarole’ experiments were made by bubbling dry air through boiling hydrochloric acid in an Erlenmeyer flask, and collecting the evolved HCl(g) in a second ‘downstream’ flask filled with distilled water. This extreme enrichment is likely due to a distillation process in which ³⁷Cl-enriched HCl(g) is stripped from the hydrochloric acid followed by a significant fraction of the light HCl(g) redissolving into the constantly condensing water vapor on the walls of the first flask. This distillation experiment creates a Cl isotope fractionation that is the same order of magnitude as observed in the high-temperature fumaroles in Central America. These results suggest that there must be a H₂O liquid-vapor region in the sub-surface fumarole conduit where light Cl is stripped from the HCl gas as it passes through the fumarole. Similar ³⁷Cl enrichments are expected in fossil epithermal boiling systems.