Isotope Fractionation Factors Research Articles

First-principles calculations of equilibrium barium isotope fractionation among silicate minerals

Barium isotopes have been increasingly applied to trace igneous processes, but the applications are hindered by the lack of equilibrium Ba isotope fractionation factors between silicate minerals. In this study, the reduced partition function ratios (RPFRs) of the major Ba-hosting silicate minerals, including muscovite, phlogopite, tremolite, richterite, microcline, albite, and epidote, as well as witherite and barite, are obtained using first-principles calculations. The effect of Ba concentration on the RPFR of those silicate minerals was examined, and no significant effect was observed in the natural range of Ba concentrations. Our results show a heavy Ba isotope enrichment order of tremolite > epidote > albite > witherite ∼ richterite ∼ microcline ∼ barite ∼ muscovite ∼ phlogopite. The RPFRs are controlled by the force constant of Ba and affected by the BaO bond length. The lattice environment of substituted ions controls the equilibrium Ba isotope fractionation factors of trace elements. For instance, the replacement of Ba2+ with Mg2+ or K+ in phlogopite introduces an approximately 0.16‰ offset of Ba isotope ratios (138Ba/134Ba) at 1000 K. In general, light Ba isotopes are preferentially enriched in crystalline positions with larger ion radii. Equilibrium Ba isotope fractionation among silicate minerals is up to 0.08‰ at 1000 K, which can be well resolved by the current analytical precision (∼0.05‰ for δ138/134Ba). This study provides the calculated equilibrium fractionation factors, which are the key to understanding Ba isotope fractionation in natural rocks.

Carbon isotope fractionation of di-(2-ethylhexyl)-phthalate during photosensitized degradation by •OH and [formula omitted] for characterization of reaction mechanisms

Compound specific isotope analysis (CSIA) has been used to explore the degradation pathways of di-(2-ethylhexyl)-phthalate (DEHP) in photosensitization experiments with H2O2, TiO2, Fe(III) and PMS. The reaction kinetics could be described by first order rate constants [k (h−1)]. The highest degradation rates of DEHP were observed in UV/TiO2 catalyzed reactions.Carbon isotope fractionation factors (εC) and metabolite patterns were examined to analyze the radical reaction. Similar carbon isotope fractionation factors have been found in experiments dominated by OH radical reactions (UV/H2O2, −1.2 ± 0.3‰; UV/TiO2, −1.2 ± 0.5‰). An inverse carbon isotope fractionation has been found for UV/Fe(III) (0.8 ± 0.2‰) and SO4− (0.6 ± 0.2‰ at pH 3) catalyzed reactions indicating that these reactions preferentially taking place at the H atom bond at tert-butyl carbon of the side chain. This hypothesis is further supported by an experiment with isobutyl acetate (εC = 1.8 ± 0.2‰ for UV/Fe(III) at pH 7; εC = -1.1 ± 0.2‰ for UV/PMS at pH 6 and εC = 0.9 ± 0.3‰ for UV/PMS at pH 3).Quenching experiments in the UV/PMS catalyzed reaction indicated a coexistence of SO4− and OH at pH 7. SO4− was the predominant radical at pH 3. OH was found to be the predominant radical in UV/H2O2 and UV/TiO2 reactions. CSIA indicated that the addition of OH to the aromatic ring of DEHP by the RAF pathway is the main reaction mechanism in UV/H2O2 and UV/TiO2 experiments. The SO4− attacks preferentially the C-H bond at the tert-butyl carbon position of the side chain of DEHP which causes an inverse carbon isotope effect. Thus, carbon isotope enrichment factors can be used to distinguish different radical reaction mechanisms with DEHP as a substrate.

Hydrogen isotope biogeochemistry of plant waxes in paired lake catchments

Studies of plants in modern catchment systems can serve to improve the level of understanding of sedimentary plant wax hydrogen isotope (δ2H) data by directly investigating the link between local climate conditions, plant source water, leaf water, and plant lipids for individual plant species. Here we present such an application to compare two lake catchments with different basin morphologies in Estonia. We sampled leaf and xylem water, as well as leaf waxes from the seven common plant species in each catchment, and soil water. We then measured the δ2H values of all waters, and of n-alkanes (δ2Hn–alk) from the plant waxes, as well as surface lake sediments. We applied a Péclet modified Craig-Gordon leaf water model using local precipitation isotope and climate data to characterize the δ2H values of the biosynthetic source water pool throughout the entire growth season. The data and model results provide a detailed view of how the input hydro-climatic signals from the precipitation δ2H values were modified by environmental and plant physiological conditions and ultimately by the biosynthetic isotope fractionation associated with n–alkane production by each plant species. We report both average apparent (εapp = –92 ± 21 ‰) and biosynthetic (εbio = –132 ± 19 ‰) hydrogen isotope fractionation factors of all species from the studied catchments. This information serves as a foundation for sedimentary organic geochemistry and paleoclimate studies, which allows for more direct and quantitative links to be made between sedimentary plant wax δ2H values and the climate signal contained in plant source water.

The influence of river-derived particles on estuarine and marine elemental cycles: Evidence from lithium isotopes

To examine the alteration of river-derived sediments through a large estuary and the implications for elemental cycling and global climate, this study analyses lithium (Li) isotopes and elemental concentrations (e.g., Li, Na, Mg, K, Ca, Fe and Al) of both the dissolved load and different phases of the sediment load (i.e., exchangeable, carbonate, oxide, clay and residue) in the Amazon estuary. The results show that river-derived sediments remove Li from the dissolved load, largely due to cation retention in secondary clays. By modelling the Li mass-balance and isotope fractionation, we estimate that the river-derived sediments gain 3–4 μg/g Li from the dissolved load in the Amazon estuary, with a Li isotope fractionation factor (αclay-solution) of approximately 0.975. Considering the whole Amazon estuary, the river-derived sediments remove around 3.6–4.8 × 109 g/yr of Li from the dissolved load. Specifically, around 1.0–1.7 × 108 g/yr of Li is removed from river water (∼1.8–3.0% of the dissolved Li discharge flux of the Amazon River) and around 3.5–4.7 × 109 g/yr of Li is removed from seawater, which represents a significant sink from the ocean. This estuarine Li sink is likely to be related to continental erosion rates; thus, continental weathering and erosion regimes could influence not only riverine Li input, but could also directly affect the Li sink, leading to a dual control on the Li budget and isotope composition in the ocean.

Temperature and organic carbon quality control the anaerobic carbon mineralization in peat profiles via modulating microbes: A case study of Changbai Mountain

Peatlands account for a significant fraction of the global carbon stock. However, the complex interplay of abiotic and biotic factors governing anaerobic carbon mineralization in response to warming remains unclear. In this study, peat sediments were collected from a typical northern peatland-Changbai Mountain to investigate the behavior and mechanism of anaerobic carbon mineralization in response to depth (0–200 cm) and temperature (5 °C, 15 °C and 20 °C), by integrating geochemical and microbial analysis. Several indices including humification indexes (HI), aromaticity, and water extractable organic carbon (WEOC) components were applied to evaluate carbon quality, while 16S rRNA sequencing was used to measure microbial composition. Regardless of temperature, degradations of carbon quality and associated reduction in microbial abundance as well as diversity resulted in a decrease in anaerobic carbon mineralization (both CO2 and CH4) towards greater depth. Warming either from 5 °C to 15 °C or 20 °C significantly increased anaerobic carbon mineralization in all depth profiles by improving carbon availability. Enhanced carbon availabilities were mediated by the change in microbial composition (p < 0.01) and an increase in metabolic activities, which was particularly evident in the enhanced β-glucosidase activity and microbial collaborations. A remarkable increase of over 10-fold in the relative abundance of the Geothrix genus was observed under warming. Overall, warming resulted in an enhanced contribution of CH4 emission and a higher ratio of hydrogenotrophic methanogenesis, as evidenced by carbon isotope fractionation factors. In addition, deep peat soils (>100 cm) with recalcitrant carbon demonstrated greater temperature sensitivity (Q10: ∼2.0) than shallow peat soils (Q10:∼1.2) when temperature increased from 15 °C to 20 °C. The findings of this study have significantly deepened our understanding for mechanisms of carbon quality and microbe-driven anaerobic carbon mineralization in peatlands under global warming.

Experimental and theoretical investigations of stable Sr isotope fractionation during its incorporation in aragonite

Strontium partitioning and isotope fractionation between aragonite and fluid have been determined experimentally at low values of the fluid saturation state (Ω) with respect to this mineral (1.1 ≤ Ωaragonite ≤ 2.2), and the measured isotope fractionation has been compared with the results of first-principles simulations. For the latter, Density Functional Theory (DFT) was used for estimation of the equilibrium Sr isotope fractionation between aragonite and Sr2+(aq). The obtained results suggest that, for values of Ωaragonite close to unity, the apparent distribution coefficient of Sr in aragonite (DSr,aragonite=XSrCO3XCaCO3(CaSr)fluid) exhibits values higher than one that rapidly decrease at increasing aragonite growth rate. Under equilibrium conditions (i.e. Ωaragonite=1) a DSr,aragoinite value of 2.7 can be extrapolated. Additionally, for aragonite growth rates rp ≤ 10−8.0±0.2 (mol/m2/s) the Sr isotope fractionation between aragonite and the fluid (i.e. Δ88/86Sraragonite-fluid) shows a constant value of −0.1±0.05‰, whereas it decreases to −0.40‰ when the growth rate increases to 10−7.7(mol/m2/s). The surface reaction kinetic model (SRKM) developed by DePaolo (2011) has been used to describe the dependence of DSr,aragonite and Δ88/86Sraragonite-fluid on mineral growth rate. In this model the best fit for DSr,aragonite and Δ88/86Sraragonite-fluid were 4 and −0.01‰, respectively, whereas the kinetic isotope fractionation factor for Δ88/86Sraragonite-fluid was −0.6‰. The results of first-principles calculations yield an equilibrium Sr isotope fractionation factor of −0.04‰ which is in excellent agreement with the experimental value of the present study. These results are the first experimental measurements of Sr isotope fractionation during inorganic aragonite precipitation as a function of growth rate and the first DFT calculations of Sr equilibrium fractionation in the aragonite-fluid system. The results of this study provide new insight into the mechanisms controlling stable Sr isotope composition in aragonite, which has implications for using Sr isotopes for paleo-reconstructions of natural archives, particularly those of abiogenic origin.

No biological effect on magnesium isotope fractionation during stromatolite growth

The growth and morphology of stromatolites have been previously linked to diverse microbial activities. However, the role of microbial metabolisms on carbonate formation during stromatolite growth remains controversial. Magnesium isotopes have been proposed to serve as a tracer of microbial carbonate formation, implying a potential biological isotope fractionation. To further elucidate whether Mg isotope fractionation is modified during microbial carbonate formation, this study reports Mg isotope compositions of Holocene stromatolites and pore waters from Lagoa Salgada, a coastal ephemeral lake in Brazil. The stromatolitic carbonates are composed of high-Mg calcites characterized by extremely positive inorganic δ13C values, up to +20‰, and variable δ26Mg values, ranging from −2.98‰ to −0.68‰. Multiple pieces of evidence consistently demonstrate changes in microbial metabolism resulting from ecosystem fluid chemistry evolution. However, the direction of Mg isotope fractionation associated with carbonate precipitation remained invariable despite the changes in microbial activity. Mineralogical features indicate that the stromatolitic carbonates formed via 2-D nucleation. An isotopic mass balance calculation based on the observed variations in δ26Mg values and Mg concentrations of associated pore waters argues for an Mg isotopic equilibrium fractionation factor of −2.60‰ associating with the carbonate formation at Lagoa Salgada, well-matching experimental values for abiotic calcite precipitation. Thus, the observed δ26Mg variability of stromatolitic carbonates is primarily controlled by the changes of physicochemical conditions in the ambient fluid. We infer that during the formation of stromatolites, a consortium of different microbial communities produces extracellular polymeric substances, which serve as a substrate for carbonate nucleation from ambient fluid without significantly affecting the fractionation of Mg-isotopes. Our findings shed light on the effects of microbial processes on carbonate formation during stromatolite growth and improve the current understanding of the Mg isotope record in microbial carbonates.

Equilibrium partitioning and isotopic fractionation of nitrogen between biotite, plagioclase, and K-feldspar during magmatic differentiation

A significant portion of the continental crust is composed of plutonic igneous rocks. However, little is known about the geochemical behaviour of N between the different minerals during magmatic differentiation. To provide new constraints for the behaviour of N during crust formation, we have characterised the geochemistry of nitrogen (N) in the compositionally zoned calc-alkaline pluton at Loch Doon, SW Scotland. We present N concentration and N isotope values for whole-rock data alongside biotite, plagioclase and K-feldspar mineral separates and assess the degree to which these data preserve equilibrium partitioning during magmatic differentiation. We show that whole rock likely inherited its N contents and δ15N signatures from the initial source composition and that this signature is homogenous at a pluton scale. Whilst the whole-rock data are best explained as crust-derived N in the source, the degree of homogenisation across a pluton scale is inconsistent with empirical N diffusivities, ruling out syn-emplacement crustal assimilation as the source of N. Instead, our data suggest a crustal signature inherited from depth associated with the Iapetus subduction zone. At a mineral scale, we find that N preferentially partitions into the feldspars over mica in this system in the order K-feldspar > plagioclase ≈ biotite > quartz, with average mineral–mineral distribution coefficients of DN plagioclase-biotite = 1.3 ± 0.6 and DN Kspar-biotite = 2.8 ± 0.6. Partitioning is accompanied by a large and near constant equilibrium isotope fractionation factor between biotite and both feldspars (averages are Δ15NPlag-Biotite = +7.8 ± 1.2‰ and Δ15NKspar-Biotite = +7.9 ± 1.0‰). In contrast, Δ15NKspar-Plagioclase closely approximates 0‰, where both minerals show δ15N values overlapping with the bulk rock δ15N values. These results show that mica crystallisation generates a 15N-depleted reservoir within plutonic rocks. Moreover, our dataset suggests that feldspars might be a more significant host of N in the igneous portion of Earth’s continental and oceanic crust than previously thought.

Stable isotopes of nitrate (δ15N and δ18O) as functional indicators of nitrogen processing in constructed wetlands

The effectiveness of nitrogen removal in wetlands relies heavily on the biological processes that control its removal. Here, we used δ15N and δ18O of nitrate (NO3−) to assess the presence and the dominance of transformation processes of nitrogen in two urban water treatment wetlands in Victoria, Australia over two rainfall events. Laboratory incubation experiments were undertaken in both light and dark to measure the isotopic fractionation factor of nitrogen assimilation (by periphyton and algae) and benthic denitrification (using bare sediment). Highest isotopic fractionations were observed for nitrogen assimilation by algae and periphyton in the light, 15ε = −14.6 to −25 ‰ while the 15ε = −1.5 ‰ in bare sediment, consistent with that of benthic denitrification. Transect water samplings of the wetlands showed different rainfall patterns (discrete versus continuous) affect the removal capability of the wetlands. During the discrete event sampling, the observed 15ε of NO3− (an average of 3.0 to 4.3 ‰) within the wetland falls between the experimental 15ε of benthic denitrification and assimilation; coinciding with the decrease in NO3− concentrations, suggesting that both denitrification and assimilation were important removal pathways. Depletion of δ15N-NO3− throughout the whole wetland system also suggested the influence of water column nitrification during this time. In contrast, during continuous rain events, no fractionation effect was observed within the wetland and was consistent with limited NO3− removal. The difference in fractionation factors within the wetland during different sampling conditions suggested that nitrate removal was highly likely limited by changes in overall nutrient inputs, residence time and water temperature which impeded biological uptake or removal. These highlight that consideration of sampling condition is crucial when assessing the efficacy of a wetland in removing nitrogen.

Li isotopic variations of particulate non-silicate phases during estuarine water mixing

Lithium (Li) isotopes are widely used as a tracer of silicate weathering. However, Li isotopic compositions of non-silicate phases carried by sediments/suspended particulate matters (SPM) have not been systematically investigated, and their influences on the dissolved δ7Li variations remain poorly known. This study investigated Li geochemistry of different chemical fractions extracted by sequential leaching from SPM samples in the Changjiang (Yangtze River) Estuary. Our results demonstrate that Fe-Mn oxyhydroxide is the largest reservoir of non-silicate Li, accounting for ∼5% of Li in bulk SPM samples. The proportion of Li in exchangeable fraction is only 0.3 ± 0.1%, with δ7Li values ranging from 14.2‰ to 22.4‰. In response to seawater addition, the dissolved Na progressively replaces physically-adsorbed Li, and the chemically-adsorbed Li is isotopically fractionated with a factor of 0.9883. In the Changjiang Estuary, δ7Li of Fe-Mn oxyhydroxides responds rapidly to environmental variations. The isotope fractionation factors during Li uptake into authigenic Fe- and Mn-oxyhydroxides are 0.981 and 0.9735, respectively. Possibly due to short interaction time during SPM passing through the estuary, incorporation of Li into non-silicate phases cannot drive a detectable alteration of estuarine dissolved Li concentration. By contrast, this study highlights the possibility of oxyhydroxides affecting δ7Li values of river water. Further work may focus on determining Li isotopic behaviours of non-silicate phases and their applications on environment reconstruction.

Site-specific isotope effect: Insights from equilibrium magnesium isotope fractionation in mantle minerals

The equilibrium Mg isotope fractionation factors (lnα) derived from both theoretical calculation and experimental studies are not consistent with the values from natural sample observation. Different cation sites, which have identical coordinate numbers but different bond strengths, contribute equally to lnα in theoretical calculations. These non-equivalent sites might complicate the inter-mineral equilibrium isotope fractionations. Here, we investigate reduced partition function ratios (103lnβ) of non-equivalent Mg cation sites (e.g., M1, M2, and M3) in forsterite, orthoenstatite, clinoenstatite, tremolite, omphacite, pyrope and spinel to explore the site-specific effect on the intra- and inter-mineral equilibrium isotope fractionation. Based on all-electron first-principles calculations, 103lnβ in different cation sites decreases in the following sequence of spinelM1 > tremoliteM3 > orthoenstatiteM1 > forsteriteM1 > omphaciteM1 > tremoliteM1 > tremoliteM2 > clinoenstatiteM1 > forsteriteM2 > orthoenstatiteM2 > clinoenstatiteM2 > pyropeM1. These findings suggest that cation sites with larger polyhedron volumes prefer light Mg isotopes to those with smaller polyhedron volumes inside minerals, and the different contribution of intra-mineral cation sites could result in significant variations of equilibrium Mg isotope fractionation between minerals. The range of lnα between minerals induced by specific sites is defined as the equilibrium isotope fractionation domain. The site-corrected lnα can provide better evaluation of equilibrium fractionation factor of Mg isotopes between spinel and olivine and clarify the equilibrium state between mantle minerals. In addition, we establish spinel-olivine and omphacite-pyrope Mg isotope geothermometers based on the site-specific effect. This study suggests that the site-specific effect could be prevalent in solid materials and thus important for other isotope systems.

One-atmosphere high-temperature CO–CO2–SO2 gas-mixing furnace: design, operation, and applications

Abstract. In this paper, we present a new design for a 1 atm gas-mixing furnace using the gas mixture CO–CO2–SO2. This furnace can simulate disequilibrium processes such as magmatic and volcanic degassing. Here, we present the technical aspects of the design. The furnace can sustain temperatures of up to 1650 ∘C and has a hot zone that spans 200 mm vertically, where the hotspot is determined to be ∼ 32 mm below the midpoint of the furnace enclosure. The four mass flow controllers are individually calibrated and accurate to within 0.8 % of the specified value. The fO2 is accurately reproduced in the furnace within ±0.002 log units, as calibrated by the Fe–FeO reaction across the iron–wüstite (IW) buffer at 1300 ∘C. The furnace can reliably simulate dynamic conditions, where the fO2 can be modulated at a maximum rate of 2.0 log units min−1 by varying the gas mixture. A delay of 40 s is observed to attain the fO2 calculated from the gas mixture, at the hotspot. A series of safety measures to protect the user from exposure to the toxic gases are detailed. In our experiments, the furnace is used to determine sulfur isotope fractionation factors among melt, sulfide, and the gas phase, within a magmatic context, using either crystals of olivine or silica glass tubes. The furnace has the potential to investigate various other dynamic high-temperature reactions occurring on Earth.

Isotopic fractionation of chlorine and potassium during chloride sublimation under lunar conditions

The Moon is depleted in volatile elements and compounds, and lunar samples exhibit a wide range of Cl isotopic compositions, which is believed to result from the volatilization of metal chlorides (e.g., NaCl, KCl, and FeCl2). However, the Cl isotopic fractionation behavior during volatilization is not well constrained, particularly for metal chlorides. Furthermore, the effect of metal chloride evaporation on metal isotopes is poorly known. In the present study, we performed NaCl and KCl sublimation experiments to study Cl and K isotopic fractionations at temperatures ranging from 923 K to 1061 K and at pressures of 7×10–5 bar to 1 bar in an N2 atmosphere. The isotope fractionation factors of 37/35Cl(αgas–solid) from NaCl sublimation experiments are 0.9985±0.0002, 0.9958±0.0004, and 0.99807±0.00004 at 1, 10–2, and 7×10–5 bar, respectively. Those of 41/39K(αgas–solid) and 37/35Cl(αgas–solid) from KCl sublimation experiments are 0.99884±0.00004 and 0.9988±0.0003 at 1 bar, 0.9977±0.0002 and 0.9972±0.0003 at 10–2 bar, and 0.9989±0.0002 and 0.9989±0.0001 at 7×10–5 bar, respectively. Chlorine and K isotopes fractionate more at 10–2 bar than at 7×10–5 bar and 1 bar. The saturation index in all the sublimation experiments was >95%, which resulted in near-equilibrium isotopic fractionation at the sublimation interface. Therefore, the isotopic fractionation was controlled by mass transfer processes in the gas and solid phases. The isotopic fractionation at 10–2 bar was controlled by the chemical diffusion of sublimated gas in an N2 atmosphere with almost no convection effect, (i.e., Pe number close to zero), whereas the isotopic fractionation at 1 bar was suppressed by atmospheric convection with a turbulence factor of 0.4±0.1 (i.e., Pe number >1). The extremely high sublimation rate and the very slow diffusion in the sublimating solid at 7×10–5 bar suppressed isotopic fractionations. Based on our experimental results, calculations using Cl/K and Na/K in lunar materials reveal that degassing of KCl contributed very little (<0.2‰) to the K isotopic fractionation (>0.58‰) during lunar magma ocean degassing. The Cl isotopic fractionation factor from lunar samples is similar to our results at 10–2 bar. This similarity of Cl isotope fractionation indicates that there may have been a transient atmosphere above the lunar magma ocean.

Bacterial physiology highlighted by the δ13C fractionation of bacteriohopanetetrol isomers

Lipid biomarkers, such as the various bacteriohopanetetrol (BHT) isomers studied here, are useful tools in tracing bacterially mediated nitrogen and carbon cycle processes affecting greenhouse gas emissions, including the anaerobic oxidation of ammonia. Three BHT isomers occur commonly in the environment. By gas chromatography, BHT-34S elutes first; it is produced by numerous bacteria. The two later eluting isomers are more constrained in their origin. The marine anammox bacteria ‘Ca. Scalindua’ is the only known producer of a BHT isomer of unknown stereochemistry (BHT-x), making BHT-x a diagnostic biomarker in anoxic marine settings. The BHT-34R isomer is produced by three freshwater aerobic heterotrophic producers (Frankia spp., Acetobacter pasteurianus, and Komagataeibacter xylinus), a freshwater serine-cycle (Type II) methanotroph (Methylocella palustris), and the freshwater anammox ‘Ca. Brocadia’, which makes the detection of freshwater anammox using BHT-34R more complicated. We investigated whether the source of BHT-34R in freshwater environments could be ascertained via its δ13C value. We used conventional on-column gas chromatography-combustion-isotope ratio mass spectrometry (GC-C-IRMS) (as opposed to high temperature GC-C-IRMS) to determine the δ13C composition of acetylated BHT isomers in cultured bacteria and bacterial enrichments. We combined these with bulk biomass and substrate δ13C compositions to establish carbon isotopic fractionation factors. The two anammox genera had large fractionation factors from dissolved inorganic carbon (DIC) to biomass (Δ13Cbiomass – DIC = –43.8 to –26.4 ‰) and to BHTs (Δ13CBHT – DIC = –53.8 to –38.2 ‰), which clearly distinguished them from the freshwater aerobic heterotrophic producers (Δ13Cbiomass – substrate = –2.3 to –0.1 ‰; Δ13CBHT – substrate = –12.8 to 5.2 ‰). Methylocella assimilated mainly carbon from DIC, rather than from methane, into its biomass and BHT, and previous work suggested this assimilation comes with relatively small fractionation. Thus, in peatlands, the BHT δ13C values of Methylocella would not reflect the low δ13C values of biogenic methane. Consequently, the presence of BHT-34R with low δ13C values relates to ‘Ca. Brocadia’ and presents a novel tool to trace anammox in freshwater environments.

First-principles calculations of equilibrium Ga isotope fractionations between several important Ga-bearing minerals and aqueous solutions

This study predicts the equilibrium isotope fractionation factors for some important Ga-bearing species, including major minerals, aqueous solutions and gas phase systems. Equilibrium isotope fractionations of Ga were investigated by using the first-principles quantum chemistry method at the B3LYP/6-311+G(d) level. The 103ln(RPFR) values of orthoclase, albite, quartz, kaolinite, forsterite, montmorillonite, gibbsite, cassiterite, aragonite, sphalerite and calcite were calculated with the volume variable cluster model. The 103ln(RPFR)s of these minerals decrease in the following order: orthoclase > albite > quartz > kaolinite > forsterite > montmorillonite > gibbsite > cassiterite > aragonite > sphalerite > calcite. The solvation effect of Ga3+-bearing aqueous species is modeled by the water-droplet method, and the 103ln(RPFR)s of Ga3+-bearing aqueous species decrease in the following order: [Ga(OH)4]− > [Ga(OH)3] > [Ga(OH)]2+ > [Ga(OH)2]+ > [Ga(H2O)6]3+. The calculation results show that equilibrium isotope fractionations of Ga between different minerals, solutions and gas phases are appreciable. Among minerals, Ga isotope fractionation exhibits the largest value between orthoclase and calcite. Ga isotopic fractionation factor between these two minerals can reach 3.18 per mil at 100 °C. Ga isotope fractionations between Ga-bearing aqueous species and minerals are important for obtaining information about the different geochemical processes, such as surficial geochemistry. This study has provided important Ga isotope fractionation factors.

Zinc Stable Isotope Fractionation Mechanisms during Adsorption on and Substitution in Iron (Hydr)oxides.

The Zn isotope fingerprint is widely used as a proxy of various environmental geochemical processes, so it is crucial to determine which are the mechanisms responsible for isotopic fractionation. Iron (Fe) (hydr)oxides greatly control the cycling and fate and thus isotope fractionation factors of Zn in terrestrial environments. Here, Zn isotope fractionation and related mechanisms during adsorption on and substitution in three FeOOH polymorphs are explored. Results demonstrate that heavy Zn isotopes are preferentially enriched onto solids, with almost similar isotopic offsets (Δ66/64Znsolid-solution = 0.25-0.36‰) for goethite, lepidocrocite, and feroxyhyte. This is consistent with the same average Zn-O bond lengths for adsorbed Zn on these solids as revealed by Zn K-edge X-ray absorption fine structure spectroscopy. In contrast, at an initial Zn/Fe molar ratio of 0.02, incorporation of Zn into goethite and lepidocrocite by substituting for lattice Fe preferentially sequesters light Zn isotopes with Δ66/64Znsubstituted-stocksolution of -1.52 ± 0.09‰ and -1.18 ± 0.15‰, while Zn-substituted feroxyhyte (0.06 ± 0.11‰) indicates almost no isotope fractionation. This is closely related to the different crystal nucleation and growth rates during the Zn-doped FeOOH formation processes. These results provide direct experimental evidence of incorporation of isotopically light Zn into Fe (hydr)oxides and improve our understanding of Zn isotope fractionation mechanisms during mineral-solution interface processes.

Boron isotope fractionation between B(OH)3 and B(OH)4− in aqueous solution: A theoretical investigation beyond the harmonic and Born−Oppenheimer approximations

Boron isotope fractionation between boric acid and borate ion in aqueous solution (hereafter, α3–4) is a key coefficient for reconstructing paleo-pH and paleo-pCO2 records, and a valuable parameter for understanding boron isotope geochemistry related to mineral-water interactions. Although boron isotopes have wide geologic application, the effect of temperature and ion pairing on boron isotope fractionation is poorly described. Moreover, conventional density functional theory (DFT) calculations do not provide an accurate estimate of the boron isotope fractionation factor. Here, we provide a new strategy to accurately calculate α3–4 values in aqueous solution by evaluating harmonic frequencies, higher-order energy terms, and solvation effects. Using benchmark coupled cluster (CCSD(T)) calculations, our 25 °C results for α3–4 ranges from 1.0259 to 1.0275, which is in good agreement with recent experimental data. Solvation effects reduced the α3–4 value by ∼6‰, while high-order corrections raised the α3–4 value by ∼2‰. Our calculations evaluating the effect of aqueous ion pairs suggest that the α3–4 values in seawater and pure water should be almost identical. Our results provide a reasonable estimation of the partition function ratios of dissolved boron species, which will benefit studies of boron isotope fractionation in aqueous environments.

Subducted carbonates not required: Deep mantle melting explains stable Ca isotopes in kimberlite magmas

Ca isotope geochemistry has great potential for improving our understanding of magmatic systems and for tracing the deep Earth carbon cycle. There are still many open questions, however, regarding the proper application of this relatively novel proxy to the study of mantle-derived magmas, including (i) the possible effects of pressure on mineral-melt Ca isotope fractionation factors, and (ii) the potential for Ca isotopes to be used as tracers of recycled marine carbonates in mantle-derived magmas. Kimberlites are mantle-derived melts that are highly enriched in CO2 and are the deepest-sourced magmas (>200 km depth) known to erupt at Earth’s surface, providing an excellent opportunity to explore these questions. We present Ca isotope data combined with detailed petrographic observations, bulk-carbonate C-O isotope data, and bulk-rock major element analyses, for a suite of 23 well-characterized kimberlite samples from their type-locality (Kimberley, South Africa). These kimberlites have abundant previous evidence for recycled surface materials in their mantle source, including low S isotope and moderately radiogenic Sr isotope compositions, yet display only limited variations in their Ca isotope compositions (δ44CaBSE of −0.08‰ to −0.27‰), with an average of −0.17 ± 0.02‰ (2SE, n = 21). This composition is indistinguishable from average carbonatites [−0.19 ± 0.03‰ (2SE, n = 106)] and OIB from recent studies [−0.16 ± 0.01‰ (2SE, n = 41)], and slightly lower than average MORB [−0.11 ± 0.02‰ (2SE, n = 31)]. Although our samples display a wide range of emplacement styles, alteration conditions, extents of magmatic differentiation, and degrees of mantle-cargo entrainment (i.e., xenocryst accumulation), we find no correlations between Ca isotopes and any of these factors. Instead, we find that low-degree partial melting of the likely kimberlite source lithology (i.e., carbon-bearing garnet lherzolite) yields modelled melt δ44CaBSE values ranging between −0.12‰ and −0.16‰ (at 1400–1500 °C), in agreement with the measured Ca isotope compositions of the Kimberley kimberlites. This observation, and the lack of heavy carbon isotope signatures in the examined samples, indicates that kimberlites do not require subducted carbonates in their mantle sources, despite their very high CO2 contents. Although several recent studies have suggested that equilibrium mineral-melt Ca isotope fractionation factors (e.g., 1000lnαgrt-melt) could be significantly different at higher pressures (i.e., due to pressure-induced changes in CaO bond lengths and coordinations), our models successfully reproduce the kimberlite data using pressure-independent predictions for mineral-melt fractionations. It remains possible, however, that differences in isotopic fractionation due to the peculiar composition of kimberlite melts (e.g., high CO2, low SiO2) are effectively cancelled out by competing pressure effects, and future work independently targeting these factors will be especially important for our understanding of Ca isotope fractionation in mantle-derived melts and the Ca isotope systematics of Earth’s mantle.

Equilibrium isotope fractionation factors of H exchange between steam and soil clay fractions.

Steam equilibration overcomes the problem of the traditional measurements of H isotope compositions, which leave an arbitrary amount of adsorbed water in the sample, by controlling for the entire exchangeable H pool, including adsorbed water and hydroxyl-H. However, the use of steam equilibration to determine nonexchangeable stable H isotope compositions in environmental media (expressed as δ2 Hn values) by mathematically eliminating the influence of exchangeable H after sample equilibration with waters of known H-isotopic composition requires the knowledge of the equilibrium isotope fractionation factor between steam-H and exchangeable H of the sample (αex-w ), which is frequently unknown. We developed a new method to determine the αex-w values for clay minerals, topsoil clay fractions, and mica by manipulating the contributions of exchangeable H to the total H pool via different degrees of post-equilibration sample drying. We measured the δ2 H values of steam-equilibrated mineral and soil samples using elemental analyzer-pyrolysis-isotope ratio mass spectrometry. The αex-w values of seven clay minerals ranged from 1.071 to 1.140, and those of 19 topsoil clay fractions ranged from 0.885 to 1.216. The αex-w value of USGS57 biotite, USGS58 muscovite, and of cellulose was 0.965, 0.871, and 1.175, respectively. The method did not work for kaolinite, because its small exchangeable H pool did not respond to the selected drying conditions. Structurally different mineral groups such as two- and three-layer clay minerals or mica showed systematically different αex-w values. The αex-w value of the topsoil clay fractions correlated with the soil clay content (r=0.63, P= 0.004), the local mean annual temperature (r = 0.68, P= 0.001), and the δ2 H values of local precipitation (r = 0.72, P< 0.001), likely to reflect the different clay mineralogy under different weathering regimes. Our new αex-w determination method yielded realistic results in line with the few previously published values for cellulose. The determined αex-w values were similar to the widely assumed values of 1.00-1.08 in the literature, suggesting that the adoption of one of these values in steam equilibration approaches is appropriate.

An investigation of zinc isotope fractionation in cacao (Theobroma cacao L.) and comparison of zinc and cadmium isotope compositions in hydroponic plant systems under high cadmium stress

This study aims to establish whether zinc (Zn) and cadmium (Cd) share similar physiological mechanisms for uptake and translocation in cacao plants (Theobroma cacao L.). Multiple-collector ICP-MS was used to determine the Zn stable isotope compositions in the roots, stems and leaves of 19 diverse cacao genotypes grown in hydroponics with 20 µmol L−1 CdCl2. Additional plants of one genotype were grown in hydroponic solutions containing lower Cd concentrations (0 and 5 µmol L−1 added CdCl2). Regardless of the Cd concentration used in the exposures, the Zn stable isotope compositions show the same systematic patterns in plant organs, with δ66Znroot > δ66Znstem > δ66Znleaf (δ66Zn denotes relative differences in 66Zn/64Zn ratios in parts per thousand). The mean Zn stable isotope fractionation between the plants and the hydroponic solutions was ε66Znuptake = –1.15 ± 0.36‰ (2SD), indicating preferential uptake of isotopically light Zn by plants from the hydroponic solution. The mean stable isotope fractionation factor associated with translocation of Zn from roots to shoots, ε66Znseq-mob = + 0.52 ± 0.36‰ (2SD), shows that isotopically heavy Zn is preferentially sequestered in the cacao roots, whilst isotopically light Zn is mobilised to the leaves. A comparison with the Cd stable isotope compositions of the same plants shows that both isotopically light Zn and Cd are preferentially taken up by cacao plants. In contrast to Zn, however, the cacao roots retain isotopically light Cd and transfer isotopically heavy Cd to the leaves.