Molybdenum isotope fractionation in pelagic euxinia: Evidence from the modern Black and Baltic Seas

Here we present stable isotope data for vertical profiles of dissolved molybdenum of the modern euxinic water columns of the Black Sea and two deeps of the Baltic Sea. Dissolved molybdenum in all water samples is depleted in salinity-normalized concentration and enriched in the heavy isotope (δ98Mo values up to + 2.9‰) compared to previously published isotope data of sedimentary molybdenum from the same range of water depths. Furthermore, δ98Mo values of all water samples from the Black Sea and anoxic deeps of the Baltic Sea are heavier than open ocean water. The observed isotope fractionation between sediments and the anoxic water column of the Black Sea are in line with the model of thiomolybdates that scavenge to particles under reducing conditions. An extrapolation to a theoretical pure MoS42− solution indicates a fractionation constant between MoS42− and authigenic solid Mo of 0.5 ± 0.3‰. Measured waters with all thiomolybdates coexisting in various proportions show larger but non-linear fractionation. The best explanation for our field observations is Mo scavenging by the thiomolybdates, dominantly — but not exclusively — present in the form of MoS42−. The Mo isotopic compositions of samples from the sediments and anoxic water column of the Baltic Sea are in overall agreement with those of the Black Sea at intermediate depth and corresponding sulphide concentrations. The more dynamic changes of redox conditions in the Baltic deeps complicate the Black Sea-derived relationship between thiomolybdates and Mo isotopic composition. In particular, the occasional flushing/mixing, of the deep waters, affects the corresponding water column and sedimentary data. δ98Mo values of the upper oxic waters of both basins are higher than predicted by mixing models based on salinity variations. The results can be explained by non-conservative behaviour of Mo under suboxic to anoxic conditions in the shallow bottom parts of the basin, most pronounced on the NW shelf of the Black Sea.

Isotope fractionation of selenium by biomethylation in microcosm incubations of soil

The natural abundance of stable Se isotopes in methylselenides reflects sources and formation conditions of methylselenides. We tested the effects of (i) different inorganic Se species spiked to soils and (ii) different soil samples on the extent of fungal biomethylation of Se and the Se isotope ratios (δ82/76Se) in methylselenides. Furthermore, we assessed the decrease of dissolved, bioavailable Se during three days of equilibration of the soils with Se-enriched solutions. We conducted closed microcosm experiments containing soil spiked with Se(IV) or Se(VI), a growth medium, and the fungus species Alternaria alternata for 11 d. The concentrations and isotope ratios of Se were determined in all components of the microcosm with multicollector ICP-MS. The equilibration of the spiked Se(IV) and Se(VI) for 3 d resulted in a decrease of dissolved, bioavailable Se concentrations by 32 to 44% and 8 to 14%, respectively. Very little isotope fractionation occurred during this phase, and it can be attributed to mixing of the added Se with the pre-existing Se in the soils and minor Se(IV) reduction in one experiment. In two of the incubated soils – moderately acidic roadside and garden soils – between 9.1 and 30% of the supplied Se(IV) and 1.7% of the supplied Se(VI) were methylated while in a strongly acidic forest soil no Se methylation occurred. The methylselenides derived from Se(IV) were strongly depleted in 82Se (δ82/76Se = − 3.3 to − 4.5‰) compared with the soil (0.16–0.45‰) and the added Se(IV) (0.20‰). The methylselenide yield of the incubations with Se(VI) was too small for isotope measurements. Our results demonstrate that Se source species and soil properties influence the extent of Se biomethylation and that the produced methylselenides contain isotopically light Se.

Magmatic–hydrothermal molybdenum isotope fractionation and its relevance to the igneous crustal signature

We analysed the Mo isotope composition of a comprehensive series of molybdenite samples from the porphyry- type Questa deposit (NM, USA), as well as one rhyolite and one granite sample, directly associated with the Mo mineralization. The δ98Mo of the molybdenites ranges between −0.48‰ and +0.40‰, with a median at −0.05‰. The median Mo isotope composition increases from early magmatic (−0.29‰) to hydrothermal (−0.05‰) breccia mineralization (median bulk breccia = −0.17‰) to late stockwork veining (+0.22‰). Moreover, variations of up to 0.34‰ are found between different molybdenite crystals within an individual hand specimen. The rhyolite sample with 0.12 μg g−1 Mo has δ98Mo = −0.57‰ and is lighter than all molybde- nites from the Questa deposit, interpreted to represent the igneous leftover after aqueous ore fluid exsolution. We recognize three Mo isotope fractionation processes that occur between about 700 and 350 °C, affecting the Mo iso- tope composition of magmatic–hydrothermal molybdenites. Δ1Mo: Minerals preferentially incorporate light Mo isotopes during progressive fractional crystallization in subvolcanic magma reservoirs, leaving behind a melt enriched in heavy Mo isotopes. Δ2Mo: Magmatic–hydrothermal fluids preferentially incorporate heavy Mo iso- topes upon fluid exsolution. Δ3Mo: Light Mo isotopes get preferentially incorporated in molybdenite during crys- tallization from an aqueous fluid, leaving behind a hydrothermal fluid that gets heavier with progressive molybdenite crystallization. The sum of all three fractionation processes produces molybdenites that record heavier δ98Mo compositions than their source magmas. This implies that the mean δ98Mo of molybdenites published so far (~0.4‰) likely represents a maximum value for the Mo isotope composition of Phanerozoic igneous upper crust.

Benchmarking the simulation of Cr isotope fractionation

A benchmark problem set consisting of four problem levels was developed for the simulation of Cr isotope fractionation in 1D and 2D domains. The benchmark is based on a recent field study where Cr(VI) reduction and accompanying Cr isotope fractionation occurs abiotically by an aqueous reaction with dissolved Fe 2+ (Wanner et al., 2012., Appl. Geochem., 27, 644–662). The problem set includes simulation of the major processes affecting the Cr isotopic composition such as the dissolution of various Cr(VI) bearing minerals, fractionation during abiotic aqueous Cr(VI) reduction, and non-fractionating precipitation of Cr(III) as sparingly soluble Cr-hydroxide. Accuracy of the presented solutions was ensured by running the problems with four well-established reactive transport modeling codes: TOUGHREACT, MIN3P, CRUNCHFLOW, and FLOTRAN. Results were also compared with an analytical Rayleigh-type fractionation model. An additional constraint on the correctness of the results was obtained by comparing output from the problem levels simulating Cr isotope fractionation with the corresponding ones only simulating bulk concentrations. For all problem levels, model to model comparisons showed excellent agreement, suggesting that for the tested geochemical processes any code is capable of accurately simulating the fate of individual Cr isotopes.

Constraints on barium isotope fractionation during aragonite precipitation by corals

We present a barium (Ba) isotope fractionation study of marine biogenic carbonates (aragonitic corals). The major aim is to provide first constraints on the Ba isotope fractionation between modern surface sea water and coral skele- ton. Mediterranean surface sea water was found to be enriched in the heavy Ba isotopes compared to previously reported values for marine open ocean authi- genic and terrestrial minerals. In aquarium experiments with a continuous sup- ply of Mediterranean surface water, the Ba isotopic composition of the bulk sample originating from cultured, aragonitic scleractinian corals (d137/134Ba between +0.16 +/- 0.12permil and +0.41 +/-0.12permil) were isotopically identical or lighter than that of the ambient Mediterranean surface sea water (d137/134Ba = +0.42 +/- 0.07permil, 2SD), which corresponds to an empirical maximum value of Ba isotope fractionation of D137/134Bacoral-seawater = -0.26 +/- 0.14permil at 25°C. This maximum Ba isotope fractionation is close and identical in direction to previous results from inorganic precipitation experiments with aragonite- structured pure BaCO3 (witherite). The variability in measured Ba concentrations of the cultured corals is at odds with a uniform distribution coefficient, DBa/Ca, thus indicating stronger vital effects on isotope than element discrimination. This observation supports the hypothesis that the Ba isotopic compositions of these corals do not result from simple equilibrium between the skeleton and the bulk sea water. Complementary coral samples from natural settings (tropical shallow-water corals from the Bahamas and Florida and cold- water corals from the Norwegian continental shelf) show an even wider range in d137/134Ba values (+0.14 +/- 0.08permil and +0.77 +/- 0.11permil), most probably due to additional spatial and/or temporal sea water heterogeneity, as indicated by recent publications.

Experimental determination of barium isotope fractionation during diffusion and adsorption processes at low temperatures

Variations in barium (Ba) stable isotope abundances measured in low and high temperature environments have recently received increasing attention. The actual processes controlling Ba isotope fractionation, however, remain mostly elusive. In this study, we present the first experimental approach to quantify the contribution of diffusion and adsorption on mass- dependent Ba isotope fractionation during transport of aqueous Ba2+ ions through a porous medium. Experiments have been carried out in which a BaCl2 solution of known isotopic composition diffused through u-shaped glass tubes filled with silica hydrogel at 10 C and 25 C for up to 201 days. The diffused Ba was highly fractionated by up to -2.15‰ in d137/134Ba, despite the low relative difference in atomic mass. The time-dependent isotope fractionation can be successfully reproduced by a diffusive transport model accounting for mass-dependent differences in the effective diffusivities of the Ba isotope species (D137Ba/D134Ba = (m134/m137)^b). Values of b extracted from the transport model were in the range of 0.010–0.011. Independently conducted batch experiments revealed that adsorption of Ba onto the surface of silica hydrogel favoured the heavier Ba isotopes (a = 1.00015 ± 0.00008). The contribution of adsorption on the overall isotope fractionation in the diffusion experiments, however, was found to be small. Our results contribute to the understanding of Ba isotope fractionation pro- cesses, which is crucial for interpreting natural isotope variations and the assessment of Ba isotope ratios as geochemical proxies.