Mercury isotope fractionation during liquid-vapor evaporation experiments

Liquid-vapor mercury isotope fractionation was investigated under equilibrium and dynamic conditions. Equilibrium evaporation experiments were performed in a closed glass system under atmospheric pressure between 0 and 22 °C, where vapor above the liquid was sampled at chemical equilibrium. Dynamic evaporation experiments were conducted in a closed glass system under 10⁻⁵ bar vacuum conditions varying (1) the fraction of liquid Hg evaporated at 22 °C and (2) the temperature of evaporation (22-100 °C). Both, residual liquid and condensed vapor fractions were analyzed using stannous chloride CV-MC-ICP-MS.Equilibrium evaporation showed a constant liquid-vapor fractionation factor (α_202/198) of 1.00086 ± 0.00022 (2SD, n = 6) within the 0-22 °C range. The 22 °C dynamic evaporations experiments displayed Rayleigh distillation fractionation behavior with liquid-vapor α_202/198 = 1.0067 ± 0.0011 (2SD), calculated from both residual and condensed vapor fractions. Our results confirm historical data (1920s) from Brönsted, Mulliken and coworkers on mercury isotopes separation using evaporation experiments, for which recalculated δ²⁰²Hg′ showed a liquid-vapor α_202/198 of 1.0076 ± 0.0017 (2SD). This liquid-vapor α_202/198 is significantly different from the expected kinetic α_202/198 value ((202/198)^0.5 = 1.0101). A conceptual evaporation model of back condensation fluxes within a thin layer at the liquid-vapor interface was used to explain this discrepancy. The δ²⁰²Hg′ of condensed vapor fractions in the 22-100 °C temperature range experiments showed a negative linear relationship with 10⁶/T², explained by increasing rates of exchange within the layer with the increase in temperature.Evaporation experiments also resulted in non-mass-dependent fractionation (NMF) of odd ¹⁹⁹Hg and ²⁰¹Hg isotopes, expressed as Δ¹⁹⁹Hg′ and Δ²⁰¹Hg′, the deviation in ‰ from the mass fractionation relationship with even isotopes. Liquid-vapor equilibrium yielded Δ¹⁹⁹Hg′/Δ²⁰¹Hg′ relationship of 2.0 ± 0.6 (2SE), which is statistically not different from the one predicted for the nuclear field shift effect (Δ¹⁹⁹Hg/Δ²⁰¹Hg ≈ 2.47). On the other hand, evaporation under dynamic conditions at 22 °C led to negative anomalies in the residual liquid fractions that are balanced by positive anomalies in condensed vapors with lower Δ¹⁹⁹Hg′/Δ²⁰¹Hg′ ratios of 1.2 ± 0.4 (2SD). This suggests that either magnetic isotope effects may have occurred without radical chemistry or an unknown NMF process on odd isotopes operated during liquid mercury evaporation.     

高温下非传统稳定同位素分馏

过去十几年来,非传统稳定同位素地球化学在高温地质过程的研究中取得了的重大进展。多接收诱导耦合等离子质谱(MC-ICP-MS)的应用引发了稳定同位素分析方法的重大突破,使得精确测定重元素的同位素比值成为可能。本文总结了以Li、Fe和Mg同位素为代表的非传统稳定同位素在岩石地球化学研究中的应用。Li同位素目前被广泛地用于地幔地球化学、俯冲带物质再循环和变质作用的研究中,可以用来示踪岩浆的源区性质和扩散等动力学过程。不同价态的Fe在矿物熔体相之间的分配可以产生Fe同位素分馏,可以发生在地幔交代、部分熔融、分离结晶等过程中。岩浆岩的Mg同位素则大致反映其源区的特征,地幔的Mg同位素组成比较均一,这为研究低温地球化学过程中Mg同位素的分馏提供一个均一的背景。此外,Cl,Si,Cu,Ca,U等等同位素体系也具有广阔的应用前景。对同位素分馏机制的实验研究和理论模拟为理解非传统稳定同位素数据提供了必要的指导。实验表明,高温下具有不同的迁移速度的轻、重同位素可以产生显著的动力学同位素分馏,这一分馏可以在化学扩散、蒸发和凝华等过程中发生;同位素在矿物和熔体以及流体相中化学环境的差异使得不同相之间可以发生平衡分馏。而最近的硅酸盐岩浆的热扩散和热迁移实验则揭示了一种"新"的岩浆分异和同位素分馏机制。沿着温度梯度,硅酸盐岩浆可以发生显著的元素和同位素分异,湿的安山岩可以通过这种方式演变成花岗质成分,因此这个过程可能对陆壳的产生和演化有重大影响。如果温度梯度在岩浆作用中能长期存在,热扩散就可以产生稳定同位素的分馏,这一机制有别于传统的平衡和动力学同位素分馏。 而多个稳定同位素体系的正相关关系是示踪热迁移过程的最有力证据。在热扩散过程中,流体承载的物质的浓度和它的索瑞系数有关。但是这个系数对体系的很多参数非常敏感,变化极大,因此对热扩散效应的研究产生极大的困难。对热扩散实验的镁、钙和铁同位素测量表明,同位素比值的变化与体系的化学组成以及总温度无关,只和温度变化的幅度有关,这意味着即使元素的索瑞系数变化多端,某一元素的同位素之间的索瑞系数的差别总为常数。这一发现有助于简化对热扩散和索瑞系数这一基础物理问题的研究 。     

Low stable carbon isotope fractionation by coccolithophore RubisCO

The ¹³C/¹²C ratio of carbon compounds is used to identify sources and sinks in the global carbon cycle. However, the relatively enriched ¹³C content observed for marine organic carbon remains enigmatic. The majority of oceanic carbon is fixed by algae and cyanobacteria via the Calvin-Benson-Bassham cycle, yet isotopic discrimination by the CO₂ fixation enzyme, RubisCO (ribulose 1,5-bisphosphate carboxylase/oxygenase), has only been measured for a single marine cyanobacterium. Different forms of RubisCO occur in different phytoplankton species (overall amino acid identity varying by as much as ∼75%) and thus may vary in the degree to which they fractionate carbon. Here we measured isotope discrimination by RubisCO from the coccolithophore Emiliania huxleyi, a cosmopolitan species used as a marine algal model.E. huxleyi RubisCO discriminated substantially less (ε = 11.1‰) against ¹³CO₂ than other RubisCO enzymes (18-29‰), despite having Michaelis-Menten kinetic parameters (K_CO2 = 72 μM; V_max = 0.66 μmol min⁻¹ mg⁻¹ protein) similar to those measured for RubisCO enzymes from different organisms. If widespread, decreased isotope discrimination of ¹³C by phytoplankton RubisCO may be a major factor influencing the enriched ¹³C content of marine organic carbon. This finding emphasizes the necessity of (a) determining ε values for RubisCOs of other marine phytoplankton and (b) re-evaluation of δ¹³C values from physiological, environmental, and geological studies.     

Late Pleistocene paleoclimates of North America as inferred from stable isotope studies of speleothems

Some speleothems (CaCO₃ cave deposits) can be demonstrated to have been formed in oxygen isotopic equilibrium with their parent seepage waters and thus a record of relative fluctuations in depositional temperature can be obtained from the δ¹⁸O variations in successive growth layers of such deposits. These temperature fluctuations reflect variations in the average annual air temperature at the surface above the cave, and therefore permit inference of past continental climate changes. Equilibrium deposits have been obtained from caves in San Luis Potosi, Bermuda, Kentucky, West Virginia, Iowa, and Alberta, ranging in age from 200,000 years BP to the present, as determined by ${}^{230}\text{Th}{}^{234}\text{>U}$ dating of the speleothems. The δ¹⁸O time curves for the six sites show the following synchronous climatic fluctuations: warm periods from 190 to 165 and from 120 to 100 Ka, at 60 and 10 Ka, and cold periods from 95 to 65 and from 55 to 20 Ka. The periods of thermal maxima correspond in time to the interglacial periods of the marine foraminiferal isotopic and faunal temperature records and to periods of high sea stand as observed from radiometric dating of raised coral reefs. Maxima and minima in insolation appear to be synchronous with this record as well.     

非传统稳定同位素分馏理论及计算

综述了稳定同位素平衡和动力学分馏的一些主要的理论方法。首先介绍了平衡分馏的理论方法,从平衡分馏的核心理论-Bigeleisen Mayer公式(或称Urey模型)以及对它的一些高级能量校正开始,介绍了基于路径积分的分子动力学方法和蒙特卡罗方法对同位素非谐效应的处理,介绍了压力效应的理论计算方法,最后介绍了核体积效应及其理论计算方法,并强调核体积效应是未来重金属同位素研究的重要部分。另外,综述了同位素动力学分馏的主要理论和计算方法,首先介绍了在稳态下因温度梯度引发的同位素分馏,着重介绍了基于局部热力学平衡理论的计算方法和最新结果;然后从如何使用过渡态理论计算化学反应导致的同位素动力学分馏出发,深入介绍了磁同位素效应和由于核体积效应引发的异常同位素动力学分馏;然后对低温下矿物生长的同位素动力学理论模型进行了介绍,尤其是仔细介绍了其中DePaolo的表面动力学模型,并对由吸附和共沉降过程产生的动力学分馏也进行了介绍。最后详细介绍了在气化过程、固体中、高温硅酸洋熔体中由于扩散引起的同位素动力学分馏,以及如何对这些过程进行理论处理,强调了开展固体矿物扩散理论计算的重要性。     

Magnesium stable isotope fractionation in marine biogenic calcite and aragonite

This survey of magnesium stable isotope compositions in marine biogenic aragonite and calcite includes samples from corals, sclerosponges, benthic porcelaneous and planktonic perforate foraminifera, coccolith oozes, red algae, and an echinoid and brachiopod test. The analyses were carried out using MC-ICP-MS with an external repeatability of ±0.22‰ (2SD for δ²⁶Mg; n = 37), obtained from a coral reference sample (JCp-1).Magnesium isotope fractionation in calcitic corals and sclerosponges agrees with published data for calcitic speleothems with an average Δ²⁶Mg_{calcite-seawater} = −2.6 ± 0.3‰ that appears to be weakly related to temperature. With one exception (Vaceletia spp.), aragonitic corals and sclerosponges also display uniform Mg isotope fractionations relative to seawater with Δ²⁶Mg_{biogenic aragonite-seawater} = −0.9 ± 0.2.Magnesium isotopes in high-Mg calcites from red algae, echinoids and perhaps some porcelaneous foraminifera as well as in all low-Mg calcites (perforate foraminifera, coccoliths and brachiopods) display significant biological influences. For planktonic foraminifera, the Mg isotope data is consistent with the fixation of Mg by organic material under equilibrium conditions, but appears to be inconsistent with Mg removal from vacuoles. Our preferred model, however, suggests that planktonic foraminifera synthesize biomolecules that increase the energetic barrier for Mg incorporation. In this model, the need to remove large quantities of Mg from vacuole solutions is avoided. For the high-Mg calcites from echinoids, the precipitation of amorphous calcium carbonate may be responsible for their weaker Mg isotope fractionation.Disregarding superimposed biological effects, it appears that cation light isotope enrichments in CaCO₃ principally result from a chemical kinetic isotope effect, related to the incorporation of cations at kink sites. In this model, the systematics of cation isotope fractionations in CaCO₃ relate to the activation energy required for cation incorporation, which probably reflects the dehydration of the cation and the crystal surface and bond formation at the incorporation site. This kinetic incorporation model predicts (i) no intrinsic dependence on growth rate, unless significant back reaction upon slow growth reduces the isotope fractionation towards that characteristic for equilibrium isotope partitioning (this may be observed for Ca isotopes in calcites), (ii) a small decrease of isotope fractionation with increasing temperature that may be amplified if higher temperatures promote back reaction and (iii) a sensitivity to changes in the activation barrier caused by additives such as anions or biomolecules or by the initial formation of amorphous CaCO₃.     

A global model of mass independent mercury stable isotope fractionation

Mass independent fractionation (MIF) of stable isotopes associated with terrestrial geochemical processes was first observed in the 1980s for oxygen and in the 1990s for sulfur isotopes. Recently mercury (Hg) was added to this shortlist when positive odd Hg isotope anomalies were observed in biological tissues. Experimental work identified photoreduction of aquatic inorganic divalent Hg^II and photodegradation of monomethylmercury species as plausible MIF inducing reactions. Observations of continental receptors of atmospheric Hg deposition such as peat, lichens, soils and, indirectly, coal have shown predominantly negative MIF. This has led to the suggestion that atmospheric Hg has negative MIF signatures and that these are the compliment of positive Hg MIF in the aquatic environment. Recent observations on atmospheric vapor phase Hg⁰ and Hg^II in wet precipitation reveal zero and positive Hg MIF respectively and are in contradiction with a simple aquatic Hg^II photoreduction scenario as the origin for global Hg MIF observations.This study presents a synthesis of all terrestrial Hg MIF observations, and these are integrated in a one-dimensional coupled continent-ocean-atmosphere model of the global Hg cycle. The model illustrates how Hg MIF signatures propagate through the various Earth surface reservoirs. The scenario in which marine photoreduction is the main MIF inducing process results in negative atmospheric Δ¹⁹⁹Hg and positive ocean Δ¹⁹⁹Hg of −0.5‰ and +0.25‰, yet does not explain atmospheric Hg⁰ and Hg^II wet precipitation observations. Alternative model scenarios that presume in-cloud aerosol Hg^II photoreduction and continental Hg^II photoreduction at soil, snow and vegetation surfaces to display MIF are necessary to explain the ensemble of natural observations. The model based approach is a first step in understanding Hg MIF at a global scale and the eventual incorporation of Hg stable isotope information in detailed global mercury chemistry and transport models.     

Determination of tin equilibrium isotope fractionation factors from synchrotron radiation experiments

A method of determination of the reduced isotopic partition function ratio (β-factor) from the partial density of state (PDOS) obtained by inelastic nuclear resonant X-ray scattering (INRXS) in synchrotron radiation experiments has been established. The method has been demonstrated by the example of tin isotopes. The tin β-factors for CaSnO₃, SnO₂, SnO have been computed from the INRXS-derived PDOSs.

     

Redox-driven stable isotope fractionation in transition metals: Application to Zn electroplating

Redox processes are ubiquitous in Earth science and are often associated with large isotope fractionations. In a previous study, voltage-dependent amplification of stable isotope fractionation was observed for an Fe reduction process. Here, we describe experiments showing a similar effect for a second transition metal, zinc. After electrochemical reduction, the composition of plated Zn metal is enriched in the light isotope (⁶⁴Zn) with respect to the Zn²⁺ leftover in solution, with a voltage-dependent fractionation factor. Results from voltage-dependent electroplating experiments are in good agreement with a second data set following equilibrium fractional isotope evolution of Zn isotopes during an electroplating process which stepwise removes most of the Zn from the aqueous reservoir. Taken together, the results indicate a voltage-dependent isotope fractionation (in permil) of ⁶⁶Zn with respect to ⁶⁴Zn to be equal to −3.45 to 1.71 V. The negative slope trend is in contrast with previously published results on iron isotope fractionation during electroplating which shows a positive slope. These results are interpreted using an extension of Marcus theory, which predicts isotope fractionations as a function of driving force in an electrochemical system. Taken together with observations of natural fractionation of redox-sensitive and non redox-active elements, our modified Marcus theory provides a framework for quantitatively predicting transition metal isotope geochemical signatures during environmentally relevant redox processes in terms of simple energetic parameters.     

Nuclear field shift effects on stable isotope fractionation: a review

An anomalous isotope effect exists in many heavy element isotope systems (e.g., Sr, Gd, Zn, U). This effect used to be called the “odd–even isotope effect” because the odd mass number isotopes behave differently from the even mass number isotopes. This mass-independent isotope fractionation driving force, which originates from the difference in the ground-state electronic energies caused by differences in nuclear size and shape, is currently denoted as the nuclear field shift effect (NFSE). It is found that the NFSE can drive isotope fractionation of some heavy elements (e.g., Hg, Tl, U) to an astonishing degree, far more than the magnitude caused by the conventional mass-dependent effect (MDE). For light elements, the MDE is the dominant factor in isotope fractionation, while the NFSE is neglectable. Furthermore, the MDE and the NFSE both decrease as temperatures increase, though at different rates. The MDE decreases rapidly with a factor of 1/T², while the NFSE decreases slowly with a factor of 1/T. As a result, even at high temperatures, the NFSE is still significant for many heavy element isotope systems. In this review paper, we begin with an introduction of the basic concept of the NSFE, including its history and recent progress, and follow with the potential implications of the inclusion of the NFSE into the kinetic isotope fractionation effect (KIE) and heavy isotope geochronology.     

Influence of the enzyme dissimilatory sulfite reductase on stable isotope fractionation during sulfate reduction

The stable isotopes of sulfate are often used as a tool to assess bacterial sulfate reduction on the macro scale. However, the mechanisms of stable isotope fractionation of sulfur and oxygen at the enzymatic level are not yet fully understood. In batch experiments with water enriched in ¹⁸O we investigated the effect of different nitrite concentrations on sulfur isotope fractionation by Desulfovibrio desulfuricans.With increasing nitrite concentrations, we found sulfur isotope enrichment factors ranging from −11.2 ± 1.8‰ to −22.5 ± 3.2‰. Furthermore, the δ¹⁸O values in the remaining sulfate increased from approximately 50-120‰ when ¹⁸O-enriched water was supplied. Since ¹⁸O-exchange with ambient water does not take place in sulfate, but rather in intermediates of the sulfate reduction pathway (e.g. ), we suggest that nitrite affects the steady-state concentration and the extent of reoxidation of the metabolic intermediate sulfite to sulfate during sulfate reduction. Given that nitrite is known to inhibit the production of the enzyme dissimilatory sulfite reductase, our results suggest that the activity of the dissimilatory sulfite reductase regulates the kinetic isotope fractionation of sulfur and oxygen during bacterial sulfate reduction. Our novel results also imply that isotope fractionation during bacterial sulfate reduction strongly depends on the cell internal enzymatic regulation rather than on the physico-chemical features of the individual enzymes.     

Climatic dependence of stable carbon and oxygen isotope signals recorded in speleothems: From soil water to speleothem calcite

Understanding the relationship between stable isotope signals recorded in speleothems (δ¹³C and δ¹⁸O) and the isotopic composition of the carbonate species in the soil water is of great importance for their interpretation in terms of past climate variability. Here the evolution of the carbon isotope composition of soil water on its way down to the cave during dissolution of limestone is studied for both closed and open-closed conditions with respect to CO₂.The water entering the cave flows as a thin film towards the drip site. CO₂ degasses from this film within approx. 10 s by molecular diffusion. Subsequently, chemical and isotopic equilibrium is established on a time scale of several 10-100 s. The δ¹³C value of the drip water is mainly determined by the isotopic composition of soil CO₂. The evolution of the δ¹⁸O value of the carbonate species is determined by the long exchange time T_ex, between oxygen in carbonate and water of several 10,000 s. Even if the oxygen of the CO₂ in soil water is in isotopic equilibrium with that of the water, dissolution of limestone delivers oxygen with a different isotopic composition changing the δ¹⁸O value of the carbonate species. Consequently, the δ¹⁸O value of the rainwater will only be reflected in the drip water if it has stayed in the rock for a sufficiently long time.After the water has entered the cave, the carbon and oxygen isotope composition of the drip water may be altered by CO₂-exchange with the cave air. Exchange times, , of about 3000 s are derived. Thus, only drip water, which drips in less than 3000 s onto the stalagmite surface, is suitable to imprint climatic signals into speleothem calcite deposited from it.Precipitation of calcite proceeds with time constants, τ_p, of several 100 s. Different rate constants and equilibrium concentrations for the heavy and light isotopes, respectively, result in isotope fractionation during calcite precipitation. Since T_ex ? τ_p, exchange with the oxygen in the water can be neglected, and the isotopic evolution of carbon and oxygen proceed analogously. For drip intervals T_d < 0.1τ_p the isotopic compositions of both carbon and oxygen in the solution evolve linearly in time. The calcite precipitated at the apex of the stalagmite reflects the isotopic signal of the drip water.For long drip intervals, when calcite is deposited from a stagnant water film, long drip intervals may have a significant effect on the isotopic composition of the DIC. In this case, the isotopic composition of the calcite deposited at the apex must be determined by averaging over the drip interval. Such processes must be considered when speleothems are used as proxies of past climate variability.     

Quantification of carbon flow from stable isotope fractionation in rice field soils with different organic matter content

Rice fields are an important source for the greenhouse gas methane produced by acetoclastic and hydrogenotrophic methanogenesis. Fractionation of ¹³C/¹²C can in principle be used to quantify the relative contribution of these pathways, but our knowledge of isotopic fractionation during reduction of CO₂ and turnover of acetate in different methanogenic environments is still scarce. We therefore measured δ¹³C signatures in two types of anoxic Italian rice field soils, one with high and one with low degradable organic matter (OM) content. Both soils were incubated in the presence and absence of methyl fluoride, a specific inhibitor of acetoclastic methanogenesis. Optimization of methyl fluoride concentration resulted in complete inhibition of acetoclastic methanogenesis. CH₄ was then exclusively produced by hydrogenotrophic methanogenesis, allowing determination of the isotopic signatures and fractionation factors specific for this methanogenic pathway. Acetate, which was then no longer consumed, accumulated and was used for determination of the isotopic signature of the fermentatively produced acetate (both total acetate and methyl carbon of acetate). Hence, all isotopic signatures, including fractionation factors were determined for the methanogenic soil. These data, were then used for computation of the relative contribution of the two methanogenic pathways. In the high OM soil, the contribution of acetoclastic methanogenesis to total CH₄ production increased simultaneously with decreasing acetate concentration. In the low OM soil, methanogenesis from H₂/CO₂ was clearly greater than theoretically expected. Furthermore, isotope fractionation of hydrogenotrophic methanogenesis indicated that the in situ energy status of methanogens strongly depended on the availability of organic carbon in the rice field soil system. Collectively, our data show that the study of isotopic fractionation in methanogenic environments allows a deeper insight into the ongoing processes, which may be quite different in the same ecosystem with different content of degradable OM.     

A new value for the stable oxygen isotope fractionation between dissolved sulfate ion and water

Although the stable oxygen isotope fractionation between dissolved sulfate ion and H₂O (hereafter ) is of physico-chemical and biogeochemical significance, no experimental value has been established until present. The primary reason being that uncatalyzed oxygen exchange between and H₂O is extremely slow, taking ∼10⁵ years at room temperature. For lack of a better approach, values of 16‰ and 31‰ at 25 °C have been assumed in the past, based on theoretical ‘gas-phase’ calculations and extrapolation of laboratory results obtained at temperatures >75 °C that actually pertain to the bisulfate system. Here I use novel quantum-chemistry calculations, which take into account detailed solute-water interactions to establish a new value for of 23‰ at 25 °C. The results of the corresponding calculations for the bisulfate ion are in agreement with observations. The new theoretical values show that sediment -data, which reflect oxygen isotope equilibration between sulfate and ambient water during microbial sulfate reduction, are consistent with the abiotic equilibrium between and water.     

石英一水一盐体系氧同位素分馏作用

温度为180—550℃,盐度(wt.％)分别为0、5、25和40条件下,在高压釜内完成了由硅胶合成石英的氧同位素分馏作用实验研究,目的是了解:①盐同位素效应;②△t值对同位素分馏的影响;③温度与同位素分馏系数的关系。研究资料表明:低温条件下矿物和纯水之间同位素平衡作用不可能发生;影响含氧矿物(初)之间氧同位素平衡速率的因素包括盐度、△t值大小和温度等;我们的研究还表明,盐度对同位素分馏作用同系数无影响,即不存在所谓的“同位素盐效应”。在180—550℃温度范围内,不同盐度条件下获得的石英-水氧同位素分馏实验方程为:10001nα_(石英-水)=3.306×10~5T~(-2)—2.71。     

Hydration in solution is critical for stable oxygen isotope fractionation between carbonate ion and water

Various isotope studies require accurate fractionation factors (α’s) between different chemical compounds in thermodynamic equilibrium. Although numerous isotope systems involve aqueous solutions, the conventional theory is formulated for the gas-phase and predicts incorrect α’s for many compounds dissolved in water. Here I show that quantum-chemistry calculations, which take into account solute–water interactions, accurately predict, for instance, oxygen isotope fractionation between dissolved and H₂O (hereafter ). Simple force field and quantum-chemistry calculations for the ‘gas-phase’ ion predict (15‰) at 25 °C. However, based on -clusters with up to 22 H₂O molecules, I calculate a value of 25‰, which agrees with the experimental value of 24.5 ± 0.5‰. Effects of geometry and anharmonicity on the calculated α were also examined. The calculations reveal the critical role of hydration in solution, which is ignored in the gas-phase theory. The approach presented provides an adequate framework for calculating fractionation factors involving dissolved compounds; it may also be used to predict α’s that cannot (or have not yet been) determined experimentally.     

Iron isotope fractionation in subterranean estuaries

Dissolved Fe concentrations in subterranean estuaries, like their river-seawater counterparts, are strongly controlled by non-conservative behavior during mixing of groundwater and seawater in coastal aquifers. Previous studies at a subterranean estuary of Waquoit Bay on Cape Cod, USA demonstrate extensive precipitation of groundwater-borne dissolved ferrous iron and subsequent accumulation of iron oxides onto subsurface sands. Waquoit Bay is thus an excellent natural laboratory to assess the mechanisms of Fe-isotope fractionation in redox-stratified environments and determine potential Fe-isotope signatures of groundwater sources to coastal seawater. Here, we report Fe isotope compositions of iron-coated sands and porewaters beneath the intertidal zone of Waquoit Bay. The distribution of pore water Fe shows two distinct sources of Fe: one residing in the upward rising plume of Fe-rich groundwater and the second in the salt-wedge zone of pore water. The groundwater source has high Fe(II) concentration consistent with anoxic conditions and yield δ⁵⁶Fe values between 0.3 and −1.3‰. In contrast, sediment porewaters occurring in the mixing zone of the subterranean estuary have very low δ⁵⁶Fe values down to −5‰. These low δ⁵⁶Fe values reflect Fe-redox cycling and result from the preferential retention of heavy Fe-isotopes onto newly formed Fe-oxyhydroxides. Analysis of Fe-oxides precipitated onto subsurface sands in two cores from the subterranean estuary revealed strong δ⁵⁶Fe and Fe concentration gradients over less than 2m, yielding an overall range of δ⁵⁶Fe values between −2 and 1.5‰. The relationship between Fe concentration and δ⁵⁶Fe of Fe-rich sands can be modeled by the progressive precipitation of Fe-oxides along fluid flow through the subterranean estuary. These results demonstrate that large-scale Fe isotope fractionation (up to 5‰) can occur in subterranean estuaries, which could lead to coastal seawater characterized by very low δ⁵⁶Fe values relative to river values.     

Stable isotope studies of fluid inclusions in speleothems and their paleoclimatic significance

Fluid inclusions found trapped in speleothems (cave deposited travertine) are interpreted as samples of seepage water from which enclosing calcium carbonate was deposited. The inclusions are assumed to have preserved their D/H ratios since the time of deposition. Initial ¹⁸O/¹⁶O ratios can be inferred from δD because rain- and snow-derived seepage waters fall on the meteoric water line (δD = 8δ¹⁸O + 10). Estimates of temperature of deposition of the carbonate can be calculated from inclusion D/H ratios and δ¹⁸O of enclosing calcite in Pleistocene speleothems. For most speleothems investigated (0–200,000 yr old) δ¹⁸O of calcite appears to have decreased with increasing temperature of deposition indicating that the dominant cause of climate-dependent change in δ¹⁸O of calcite was the change in K_cw, the isotope fractionation equilibrium constant, with temperature; δ¹⁸O of meteoric precipitation generally increased with increasing temperature, but not sufficiently to compensate for the decrease in K_cw.     

Calculation of sulfur isotope fractionation in sulfides

The increment method has been successfully applied to calculate thermodynamic isotope fractionation factors of oxygen in silicates, oxides, carbonates, and sulfates. In this paper, we modified the increment method to calculate thermodynamic isotope fractionation factors of sulfur in sulfides, based on chemical features of sulfur-metal bonds and crystal features of sulfide minerals. To approximate the bond strength of sulfides, a new constant, known as the Madelung constant, was introduced. The increment method was then extended to calculate the reduced partition function ratios of sphalerite, chalcopyrite, galena, pyrrhotite, greenockite, bornite, cubanite, sulvanite, and violarite, as well as the isotope fractionation factors between them over the temperature range from 0 to 1000 °C. The order of ³⁴S enrichment in these nine minerals is pyrrhotite > greenockite > sphalerite > chalcopyrite > cubanite > sulvanite > bornite > violarite > galena. Our improved method constitutes another model for calculating the thermodynamic isotope fractionation factors of sulfur in sulfides of geochemical interest.