Kinetic O effects in the hydrologic cycle: Indirect evidence and implications

The abundances of ¹⁸O and deuterium in the present and past hydrologic cycle have proven to be an important tool in Earth systems science. In contrast, the abundance of ¹⁷O in precipitation has thus far been assumed to carry no additional information to that of ¹⁸O. Here, we demonstrate, using known constraints on oxygen isotope abundances from the O₂ cycle and existing data about the natural abundance of ¹⁷O in water, that the relationship between the discrimination against ¹⁷O and ¹⁸O in water may vary. This relationship, presented here as θ = ln (¹⁷α)/ln (¹⁸α), is found to be 0.511 ± 0.005 for kinetic transport effects and 0.526 ± 0.001 for equilibrium effects, with very low temperature sensitivity. As a result, the ¹⁷Δ of precipitation is controlled primarily by kinetic effects during evaporation of the initial vapor and, in contrast to the deuterium excess, is independent of the temperature at the evaporation (and condensation) site. This makes ¹⁷Δ a unique tracer that complements ¹⁸O and deuterium, and may allow for a decoupling of changes in the temperature of the ocean, that serves as the vapor source, from changes in the relative humidity above it. In addition, the ¹⁷Δ of ice caps is influenced by the kinetic effects in ice formation, and therefore measurement of ice ¹⁷Δ can be used as an additional constraint for better understanding and parameterization of these effects.     

The triple isotopic composition of oxygen in leaf water

The isotopic composition of atmospheric O₂ depends on the rates of oxygen cycling in photosynthesis, respiration, photochemical reactions in the stratosphere and on δ¹⁷O and δ¹⁸O of ocean and leaf water. While most of the factors affecting δ¹⁷O and δ¹⁸O of air O₂ have been studied extensively in recent years, δ¹⁷O of leaf water—the substrate for all terrestrial photosynthesis—remained unknown. In order to understand the isotopic composition of atmospheric O₂ at present and in fossil air in ice cores, we studied leaf water in field experiments in Israel and in a European survey. We measured the difference in δ¹⁷O and δ¹⁸O between stem and leaf water, which is the result of isotope enrichment during transpiration. We calculated the slopes of the lines linking the isotopic compositions of stem and leaf water. The obtained slopes in ln(δ¹⁷O + 1) vs. ln(δ¹⁸O + 1) plots are characterized by very high precision (∼0.001) despite of relatively large differences between duplicates in both δ¹⁷O and δ¹⁸O (0.02-0.05‰). This is so because the errors in δ¹⁸O and δ¹⁷O are mass-dependent. The slope of the leaf transpiration process varied between 0.5111 ± 0.0013 and 0.5204 ± 0.0005, which is considerably smaller than the slope linking liquid water and vapor at equilibrium (0.529). We further found that the slope of the transpiration process decreases with atmospheric relative humidity (h) as 0.522-0.008 × h, for h in the range 0.3-1. This slope is neither influenced by the plant species, nor by the environmental conditions where plants grow nor does it show strong variations along long leaves.     

A general algorithm for the O abundance correction to C/C determinations from CO2 isotopologue measurements, including CO2 characterised by ‘mass-independent’ oxygen isotope distributions

It is widely recognised that a significant limitation to the ultimate precision of carbon stable isotope ratio measurements, as obtained from dual-inlet mass spectrometric measurements of CO₂ isotopologue ion abundances at m/z 44, 45, and 46, is the correction for interference from ¹⁷O-bearing molecular ions. Two long-established, alternative procedures for determining the magnitude of this correction are in widespread use (although only one has IAEA approval); their differences lead to small but potentially significant discrepancies in the magnitude of the resulting correction. Furthermore, neither approach was designed to accommodate oxygen three-isotope distributions which do not conform to terrestrial mass-dependent behaviour. Stratospheric CO₂, for example, contains a strongly ‘mass-independent’ oxygen isotope composition. A new strategy for determining the ¹⁷O-bearing ion correction is presented, for application where the oxygen three-isotope characteristics of the analyte CO₂ are accurately known (or assigned) in terms of the slope λ of the three-isotope fractionation line and the ordinate axis intercept 10³ ln(1 + k) on a 10³ ln(1 + δ¹⁷O) versus 10³ ln(1 + δ¹⁸O) plot. At the heart of the approach is the relationship between ¹⁷R, which is the ¹⁷O/¹⁶O ratio of the sample CO₂, and other assigned or empirically determined parameters needed for the δ¹³C evaluation:

     

Kinetic and equilibrium mass-dependent isotope fractionation laws in nature and their geochemical and cosmochemical significance

The mass-dependent fractionation laws that describe the partitioning of isotopes are different for kinetic and equilibrium reactions. These laws are characterized by the exponent relating the fractionation factors for two isotope ratios such that α_2/1 = α_3/1^β. The exponent β for equilibrium exchange is (1/m₁ − 1/m₂)/(1/m₁ − 1/m₃), where m_i are the atomic masses and m₁ < m₂ < m₃. For kinetic fractionation, the masses used to evaluate β depend upon the isotopic species in motion. Reduced masses apply for breaking bonds whereas molecular or atomic masses apply for transport processes. In each case the functional form of the kinetic β is ln(M₁/M₂)/ln(M₁/M₃), where M_i are the reduced, molecular, or atomic masses. New high-precision Mg isotope ratio data confirm that the distinct equilibrium and kinetic fractionation laws can be resolved for changes in isotope ratios of only 3‰ per amu. The variability in mass-dependent fractionation laws is sufficient to explain the negative Δ¹⁷O of tropospheric O₂ relative to rocks and differences in Δ¹⁷O between carbonate, hydroxyl, and anhydrous silicate in Martian meteorites. (For simplicity, we use integer amu values for masses when evaluating β throughout this paper.)     

The non-mass-dependent oxygen isotope effect in the electrodissociation of carbon dioxide: A step toward understanding NoMaD chemistry

A non-mass dependent (NoMaD) oxygen isotope effect is demonstrated in the dissociation of CO₂ similar to that observed in the electrosynthesis of ozone. The molecular oxygen produced carries the signature of two separate isotopic fractionation processes; a mass-dependent fractionation probably due to CO₂ + O isotopic exchange, and a secondary NoMaD fractionation (δ¹⁷O = 0.97 ± 0.09δ¹⁸O, with the O₂ depleted in ¹⁷O and ¹⁸O). It is suggested that the effect is due to either the formation or relaxation of ozone in an excited electronic state. This represents the latest advance in the understanding of chemical NoMaD effects which may be essential to the explanation of non-mass-dependent fractionations observed in meteorites.     

Carbon isotopic fractionation during the CO2 exchange process between air and sea water under equilibrium and kinetic conditions

Carbon isotopic fractionation during the air/sea exchange process is not fully understood at present. Information on the equilibrium and kinetic fractionation factors is an essential requirement, together with the value of the CO₂ partial pressure, for understanding the carbon cycle in the atmosphere and marine environments. Using a specially designed countercurrent equilibrator system, the fractionation factors between gaseous CO₂ and dissolved inorganic carbon in sea water were determined under both kinetic and equilibrium conditions. The following results were obtained: kinetic fractionation factor for air to sea (α_as) is 0.998 at 288.2 K; kinetic fractionation factor for sea to air (α_sa) is 0.990; equilibrium fractionation factor (α_eq) is 0.991 at pH = 8.3 and 288.2 K. From these results, the carbon isotopic ratio of CO₂ passed through the air/ sea interface is estimated to be about ?10 %. for air to sea and ?8 %. for sea to air when CO₂ exchange takes place between air (δ¹³C = ?8 %.) and surface sea water (δ¹³C = 2 %.) at 288.2 K.     

Determining the O/O ratio of water using a water-CO2 equilibration method: Application to glacial-interglacial changes in O-excess from the Dome Fuji ice core, Antarctica

Recent studies show that oxygen three isotope measurement (¹⁶O, ¹⁷O, and ¹⁸O) of water provides additional information for investigating the hydrological cycle and paleoclimate. For determining the ¹⁸O/¹⁶O value of water, a conventional CO₂-water equilibration method involves measurement of the ratios of CO₂ isotopologues which were equilibrated with water. However, this long-established technique was not intended to measure the ¹⁷O/¹⁶O ratio, primarily because the historic ion correction scheme does not allow for possible deviations from a fixed (and mass-dependent) relationship between ¹⁷O/¹⁶O and ¹⁸O/¹⁶O isotope ratios. Here, we propose an improved method for obtaining the ¹⁷O/¹⁶O isotope ratio of fresh water by the equilibration method and measurement of the 45/44 CO₂ ion abundance ratio. Equations which we formulated for ¹⁷O/¹⁶O measurement have two features: first, instead of absolute isotope ratio (R), all equations are formulated in δ values, measured by isotope ratio mass spectrometry. Second, we include two “assigned” δ values of water standards in the equations, because the δ¹⁸O are commonly measured against two working standards to normalize the span of the δ scale. This approach clarifies that the contribution from ¹⁷O (¹²C¹⁶O¹⁷O⁺) to the molecular ion current at mass-to-charge ratio m/z 45 signal depends not on the absolute ¹³C/¹²C ratio, but on the relative δ¹³C differences between the working standards and the sample. The pH value of water affects δ¹⁷O estimation because δ¹³C of CO₂ was changed in the water-CO₂ system. We reevaluated this effect using a set of equations, which explicitly includes CO₂ partial pressure effect on pH value. Our new estimation of pH effect is significantly smaller than previously reported value, but it does not alter the main conclusions in the previous study. The method was verified by δ¹⁷O measurements of an international standard reference water (GISP) provided by the IAEA. We applied the method to investigate ¹⁷O-excess of the ice core drilled at the Dome Fuji station, Antarctica. A total of 1320 samples from a 130 m section around Marine Isotope Stage 9.3 (∼330,000 years before present) were measured. The error of a measurement for δ¹⁷O is 0.175‰ and that of ¹⁷O-excess is 184 per meg. Although these analytical uncertainties hampered accurate estimation of the changes in ¹⁷O-excess, the averaged data indicate that ¹⁷O-excess around MIS 9.3 was higher than during the subsequent glacial period. This approach can be applied only to fresh water samples, and additional improvements will be needed to measure samples which contains significant amount of carbonate minerals.     

Oxygen isotope compositions of phosphate from arvicoline teeth and Quaternary climatic changes, Gigny, French Jura

Oxygen isotope compositions of biogenic phosphates from mammals are widely used as proxies of the isotopic compositions of meteoric waters that are roughly linearly related to the air temperature at high- and mid-latitudes. An oxygen isotope fractionation equation was determined by using present-day European arvicoline (rodents) tooth phosphate: δ¹⁸O_p = 20.98(±0.59) + 0.572(±0.065) δ¹⁸O_w. This fractionation equation was applied to the Late Pleistocene karstic sequence of Gigny, French Jura. Comparison between the oxygen isotope compositions of arvicoline tooth phosphate and Greenland ice core records suggests to reconsider the previously established hypothetical chronology of the sequence. According to the δ¹⁸O value of meteoric water-mean air temperature relationships, the δ¹⁸O value of arvicoline teeth records variations in mean air temperatures that range from 0° to 15°C.     

Model for kinetic effects on calcium isotope fractionation (δCa) in inorganic aragonite and cultured planktonic foraminifera

The calcium isotope ratios (δ⁴⁴Ca = [(⁴⁴Ca/⁴⁰Ca)_sample/(⁴⁴Ca/⁴⁰Ca)_standard −1] · 1000) of Orbulina universa and of inorganically precipitated aragonite are positively correlated to temperature. The slopes of 0.019 and 0.015‰ °C⁻¹, respectively, are a factor of 13 and 16 times smaller than the previously determined fractionation from a second foraminifera, Globigerinoides sacculifer, having a slope of about 0.24‰ °C⁻¹. The observation that δ⁴⁴Ca is positively correlated to temperature is opposite in sign to the oxygen isotopic fractionation (δ¹⁸O) in calcium carbonate (CaCO₃). These observations are explained by a model which considers that Ca²⁺-ions forming ionic bonds are affected by kinetic fractionation only, whereas covalently bound atoms like oxygen are affected by kinetic and equilibrium fractionation. From thermodynamic consideration of kinetic isotope fractionation, it can be shown that the slope of the enrichment factor α(T) is mass-dependent. However, for O. universa and the inorganic precipitates, the calculated mass of about 520 ± 60 and 640 ± 70 amu (atomic mass units) is not compatible with the expected ion mass for ⁴⁰Ca and ⁴⁴Ca. To reconcile this discrepancy, we propose that Ca diffusion and δ⁴⁴Ca isotope fractionation at liquid/solid transitions involves Ca²⁺-aquocomplexes (Ca[H₂O]_n²⁺ · mH₂O) rather than pure Ca²⁺-ion diffusion. From our measurements we calculate that such a hypothesized Ca²⁺-aquocomplex correlates to a hydration number of up to 25 water molecules (490 amu). For O. universa we propose that their biologically mediated Ca isotope fractionation resembles fractionation during inorganic precipitation of CaCO₃ in seawater. To explain the different Ca isotope fractionation in O. universa and in G. sacculifer, we suggest that the latter species actively dehydrates the Ca²⁺-aquocomplex before calcification takes place. The very different temperature response of Ca isotopes in the two species suggests that the use of δ⁴⁴Ca as a temperature proxy will require careful study of species effects.     

Oxygen isotope fractionation in synthetic magnesian calcite

Mg-bearing calcite was precipitated at 25°C in closed system free-drift experiments from solutions containing NaHCO₃, CaCl₂ and MgCl₂. The chemical and isotope composition of the solution and precipitate were investigated during time course experiments of 24-h duration. Monohydrocalcite and calcite precipitated early in the experiments (<8 h), while Mg-calcite was the predominant precipitate (>95%) thereafter. Solid collected at the end of the experiments displayed compositional zoning from pure calcite in crystal cores to up to 23 mol% MgCO₃ in the rims. Smaller excursions in Mg were superimposed on this chemical record, which is characteristic of oscillatory zoning observed in synthetic and natural solid-solution carbonates of differing solubility. Magnesium also altered the predominant morphology of crystals over time from the {104} to {100} and {110} growth forms.The oxygen isotope fractionation factor for the magnesian-calcite-water system (as 10³lnα_Mg-cl-H2O) displayed a strong dependence on the mol% MgCO₃ in the solid phase, but quantification of the relationship was difficult due to the heterogeneous nature of the precipitate. Considering only the Mg-content and δ¹⁸O values for the bulk solid, 10³lnα_Mg-cl-H2O increased at a rate of 0.17 ± 0.02 per mol% MgCO₃; this value is a factor of three higher than the single previous estimate (Tarutani T., Clayton R.N., and Mayeda T. K. (1969) The effect of polymorphims and magnesium substitution on oxygen isotope fractionation between calcium carbonate and water. Geochim. Cosmochim. Acta 33, 987-996). Nevertheless, extrapolation of our relationship to the pure calcite end member yielded a value of 27.9 ± 0.02, which is similar in magnitude to published values for the calcite-water system. Although no kinetic effect was observed on 10³lnα_Mg-cl-H2O for precipitation rates that ranged from 10^3.21 to 10^4.60 μmol · m⁻² · h⁻¹, it was impossible to disentangle the potential effect(s) of precipitation rate and Mg-content on 10³lnα_Mg-cl-H2O due to the heterogeneous nature of the solid.The results of this study suggest that paleotemperatures inferred from the δ¹⁸O values of high magnesian calcite (>10 mol% MgCO₃) may be significantly underestimated. Also, the results underscore the need for additional experiments to accurately characterize the effect of Mg coprecipitation on the isotope systematics of calcite from a chemically homogeneous precipitate or a heterogeneous material that is analyzed at the scale of chemical and isotopic zonation.     

Kinetic and equilibrium Fe isotope fractionation between aqueous Fe(III) and hematite

Application of the Fe isotope system to studies of natural rocks and fluids requires precise knowledge of equilibrium Fe isotope fractionation factors among various aqueous Fe species and minerals. These are difficult to obtain at the low temperatures at which Fe isotope fractionation is expected to be largest and requires careful distinction between kinetic and equilibrium isotope effects. A detailed investigation of Fe isotope fractionation between [Fe^III(H₂O)₆]³⁺ and hematite at 98°C allows the equilibrium ⁵⁶Fe/⁵⁴Fe fractionation to be inferred, which we estimate at 10³lnα_{Fe(III)-hematite} = −0.10 ± 0.20‰. We also infer that the slope of Fe(III)-hematite fractionation is modest relative to 10⁶/T², which would imply that this fractionation remains close to zero at lower temperatures. These results indicate that Fe isotope compositions of hematite may closely approximate those of the fluids from which they precipitated if equilibrium isotopic fractionation is assumed, allowing inference of δ⁵⁶Fe values of ancient fluids from the rock record. The equilibrium Fe(III)-hematite fractionation factor determined in this study is significantly smaller than that obtained from the reduced partition function ratios calculated for [Fe^III(H₂O)₆]³⁺ and hematite based on vibrational frequencies and Mössbauer shifts by [Polyakov 1997] and [Polyakov and Mineev 2000], and Schauble et al. (2001), highlighting the importance of experimental calibration of Fe isotope fractionation factors. In contrast to the long-term (up to 203 d) experiments, short-term experiments indicate that kinetic isotope effects dominate during rapid precipitation of ferric oxides. Precipitation of hematite over ∼12 h produces a kinetic isotope fractionation where 10³lnα_{Fe(III)-hematite} = +1.32 ± 0.12‰. Precipitation under nonequilibrium conditions, however, can be recognized through stepwise dissolution in concentrated acids. As expected, our results demonstrate that dissolution by itself does not measurably fractionate Fe isotopes.     

Kinetic oxygen isotope effects during dissimilatory sulfate reduction: A combined theoretical and experimental approach

Kinetic isotope effects related to the breaking of chemical bonds drive sulfur isotope fractionation during dissimilatory sulfate reduction (DSR), whereas oxygen isotope fractionation during DSR is dominated by exchange between intercellular sulfur intermediates and water. We use a simplified biochemical model for DSR to explore how a kinetic oxygen isotope effect may be expressed. We then explore these relationships in light of evolving sulfur and oxygen isotope compositions (δ³⁴S_SO4 and δ¹⁸O_SO4) during batch culture growth of twelve strains of sulfate-reducing bacteria. Cultured under conditions to optimize growth and with identical δ¹⁸O_H2O and initial δ¹⁸O_SO4, all strains show ³⁴S enrichment, whereas only six strains show significant ¹⁸O enrichment. The remaining six show no (or minimal) change in δ¹⁸O_SO4 over the growth of the bacteria. We use these experimental and theoretical results to address three questions: (i) which sulfur intermediates exchange oxygen isotopes with water, (ii) what is the kinetic oxygen isotope effect related to the reduction of adenosine phosphosulfate (APS) to sulfite (SO₃²⁻), (iii) does a kinetic oxygen isotope effect impact the apparent oxygen isotope equilibrium values? We conclude that oxygen isotope exchange between water and a sulfur intermediate likely occurs downstream of APS and that our data constrain the kinetic oxygen isotope fractionation for the reduction of APS to sulfite to be smaller than 4‰. This small oxygen isotope effect impacts the apparent oxygen isotope equilibrium as controlled by the extent to which APS reduction is rate-limiting.     

N/N and O/O stable isotope ratios of nitrous oxide produced during denitrification in temperate forest soils

Anaerobic incubations of upland and wetland temperate forest soils from the same watershed were conducted under different moisture and temperature conditions. Rates of nitrous oxide (N₂O) production by denitrification of nitrate () and the stable isotopic composition of the N₂O (δ¹⁵N, δ¹⁸O) were measured. In all soils, N₂O production increased with elevated temperature and soil moisture. At each temperature and moisture level, the rate of N₂O production in the wetland soil was greater than in the upland soil. The ¹⁵N isotope effect (ε) _{(product − substrate)} ranged from −20‰ to −29‰. These results are consistent with other published estimates of ¹⁵N fractionation from both single species culture experiments and soil incubation studies from different ecosystems.A series of incubations were conducted with ¹⁸O-enriched water (H₂O) to determine if significant oxygen exchange (O-exchange) occurred between H₂O and N₂O precursors during denitrification. The exchange of H₂O-O with nitrite () and/or nitric oxide (NO) oxygen has been documented in single organism culture studies but has not been demonstrated in soils prior to this study. The fraction of N₂O-O derived from H₂O-O was confined to a strikingly narrow range that differed between soil types. H₂O-O incorporation into N₂O produced from upland and wetland soils was 86% to 94% and 64% to 70%, respectively. Neither the temperature, soil moisture, nor the rate of N₂O production influenced the magnitude of O-exchange. With the exception of one treatment, the net ¹⁸O isotope effect (ε_net) _{(product-substrate)} ranged from +37‰ to +43‰.Most previous studies that have reported ¹⁸O isotope effects for denitrification of to N₂O have failed to account for the effect of oxygen exchange with H₂O. When high amounts of O-exchange occur after fractionation during reductive O-loss, the ¹⁸O-enrichment is effectively lost or diminished and δ¹⁸O-N₂O values will be largely dictated by δ¹⁸O-H₂O values and subsequent fractionation. The process and extent of O-exchange, combined with the magnitude of oxygen isotope fractionation at each reduction step, appear to be the dominant controls on the observed oxygen isotope effect. In these experiments, significant oxygen isotope fractionation was observed to occur after the majority of water O-exchange. Due to the importance of O-exchange, the net oxygen isotope effect for N₂O production in soils can only be determined using δ¹⁸O-H₂O addition experiments with δ¹⁸O-H₂O close to natural abundance.The results of this study support the continued use of δ¹⁵N-N₂O analysis to fingerprint N₂O produced from the denitrification of . The utilization of ¹⁸O/¹⁶O ratios of N₂O to study N₂O production pathways in soil environments is complicated by oxygen exchange with water, which is not usually quantified in field studies. The oxygen isotope fractionation observed in this study was confined to a narrow range, and there was a clear difference in water O-exchange between soil types regardless of temperature, soil moisture, and N₂O production rate. This suggests that ¹⁸O/¹⁶O ratios of N₂O may be useful in characterizing the actively denitrifying microbial community.     

Oxygen isotope exchange between H2O and CO2 at elevated CO2 pressures: Implications for monitoring of geological CO2 storage

Traditionally, the application of stable isotopes in Carbon Capture and Storage (CCS) projects has focused on δ¹³C values of CO₂ to trace the migration of injected CO₂ in the subsurface. More recently the use of δ¹⁸O values of both CO₂ and reservoir fluids has been proposed as a method for quantifying in situ CO₂ reservoir saturations due to O isotope exchange between CO₂ and H₂O and subsequent changes in δ¹⁸O_H2O values in the presence of high concentrations of CO₂. To verify that O isotope exchange between CO₂ and H₂O reaches equilibrium within days, and that δ¹⁸O_H2O values indeed change predictably due to the presence of CO₂, a laboratory study was conducted during which the isotope composition of H₂O, CO₂, and dissolved inorganic C (DIC) was determined at representative reservoir conditions (50 °C and up to 19 MPa) and varying CO₂ pressures. Conditions typical for the Pembina Cardium CO₂ Monitoring Pilot in Alberta (Canada) were chosen for the experiments. Results obtained showed that δ¹⁸O values of CO₂ were on average 36.4 ± 2.2‰ (1σ, n = 15) higher than those of water at all pressures up to and including reservoir pressure (19 MPa), in excellent agreement with the theoretically predicted isotope enrichment factor of 35.5‰ for the experimental temperatures of 50 °C. By using ¹⁸O enriched water for the experiments it was demonstrated that changes in the δ¹⁸O values of water were predictably related to the fraction of O in the system sourced from CO₂ in excellent agreement with theoretical predictions. Since the fraction of O sourced from CO₂ is related to the total volumetric saturation of CO₂ and water as a fraction of the total volume of the system, it is concluded that changes in δ¹⁸O values of reservoir fluids can be used to calculate reservoir saturations of CO₂ in CCS settings given that the δ¹⁸O values of CO₂ and water are sufficiently distinct.     

Role of basicity and tetrahedral speciation in controlling the thermodynamic properties of silicate liquids, part 1: the system CaO-MgO-Al2O3-SiO2

Activity coefficients of oxide components in the system CaO-MgO-Al2O3-SiO2 (CMAS) were calculated with the model of Berman (Berman R. G., “A thermodynamic model for multicomponent melts with application to the system CaO-MgO-Al2O3-SiO2,” Ph.D. dissertation, University of British Columbia, 1983) and used to explore large-scale relationships among these variables and between them and the liquid composition. On the basis of Berman’s model, the natural logarithm of the activity coefficient of MgO, ln(γMgOLiq), and ln(γMgOLiq/γSiO2Liq) are nearly linear functions of ln(γCaOLiq). All three of these variables are simple functions of the optical basicity Λ with which they display minima near Λ ∼ 0.54 that are generated by liquids with low ratios of nonbridging to tetrahedral oxygens (NBO/T) (<0.3) and a mole fraction ratio, XSiO2Liq/XAl2O3Liq, in the range 4 to 20. Variations in ln(γCaOLiq) at constant Λ near the minimum are due mostly to liquids with (XCaOLiq + XMgOLiq)/XAl2O3Liq < 1. The correlations with optical basicity imply that the electron donor power is an important factor in determining the thermodynamic properties of aluminosilicate liquids.For a constant NBO/T, ln(γCaOLiq/γAl2O3Liq) and ln(γMgOLiqγAl2O3Liq) form curves in terms of XSiO2Liq/XAl2O3Liq. The same liquids that generate minima in the Λ plots are also associated with minima in ln(γCaOLiqγAl2O3Liq) and ln(γMgOLiqγAl2O3Liq) as a function of XSiO2Liq/XAl2O3Liq. In addition, there are maxima or sharp changes in slope for NBO/T > 0.3, which occur for XSiO2Liq/XAl2O3Liq ranging from ∼0 to ∼6 and increase with increasing NBO/T. The systematic variations in activity coefficients as a function of composition and optical basicity reflect underlying shifts in speciation as the composition of the liquid is changed. On the basis of correlations among the activity coefficients, it is likely that the use of CaO, an exchange component such as SiMg−1 and two of MgO, CaAl2O4, or MgAl2O4 would yield significant savings in the number of parameters required to model the excess free energy surface of liquids over large portions of CMAS relative to the use of oxide end members.Systematic behavior of thermodynamic properties extends to small amounts of other elements dissolved in otherwise CMAS liquids. For example, ln(XFe2+Liq/XFe3+Liq) at constant oxygen fugacity is linearly correlated with ln(γCaOLiq). Similarly, ln(CS), where CS is the sulfide capacity is linearly correlated at constant temperature with each of the optical basicity, ln(aCaOLiq) and ln(γCaOLiq), although the correlation for the latter breaks down for low values of Λ. The well-known systematic behavior of sulfide capacity as a function of optical basicity for systems inside as well as outside CMAS suggests that ln(γCaSLiq) is also a simple function of optical basicity and that the relationships observed among the activity coefficients in CMAS may hold for more complex systems.     

Abundance of mass 47 CO2 in urban air, car exhaust, and human breath

Atmospheric carbon dioxide is widely studied using records of CO₂ mixing ratio, δ¹³C and δ¹⁸O. However, the number and variability of sources and sinks prevents these alone from uniquely defining the budget. Carbon dioxide having a mass of 47 u (principally ¹³C¹⁸O¹⁶O) provides an additional constraint. In particular, the mass 47 anomaly (Δ₄₇) can distinguish between CO₂ produced by high temperature combustion processes vs. low temperature respiratory processes. Δ₄₇ is defined as the abundance of mass 47 isotopologues in excess of that expected for a random distribution of isotopes, where random distribution means that the abundance of an isotopologue is the product of abundances of the isotopes it is composed of and is calculated based on the measured ¹³C and ¹⁸O values. In this study, we estimate the δ¹³C (vs. VPDB), δ¹⁸O (vs. VSMOW), δ47, and Δ₄₇ values of CO₂ from car exhaust and from human breath, by constructing ‘Keeling plots’ using samples that are mixtures of ambient air and CO₂ from these sources. δ47 is defined as , where is the R⁴⁷ value for a hypothetical CO₂ whose δ¹³C_VPDB = 0, δ¹⁸O_VSMOW = 0, and Δ₄₇ = 0. Ambient air in Pasadena, CA, where this study was conducted, varied in [CO₂] from 383 to 404 μmol mol⁻¹, in δ¹³C and δ¹⁸O from −9.2 to −10.2‰ and from 40.6 to 41.9‰, respectively, in δ47 from 32.5 to 33.9‰, and in Δ₄₇ from 0.73 to 0.96‰. Air sampled at varying distances from a car exhaust pipe was enriched in a combustion source having a composition, as determined by a ‘Keeling plot’ intercept, of −24.4 ± 0.2‰ for δ¹³C (similar to the δ¹³C of local gasoline), δ¹⁸O of 29.9 ± 0.4‰, δ47 of 6.6 ± 0.6‰, and Δ₄₇ of 0.41 ± 0.03‰. Both δ¹⁸O and Δ₄₇ values of the car exhaust end-member are consistent with that expected for thermodynamic equilibrium at∼200 °C between CO₂ and water generated by combustion of gasoline-air mixtures. Samples of CO₂ from human breath were found to have δ¹³C and δ¹⁸O values broadly similar to those of car exhaust-air mixtures, −22.3 ± 0.2 and 34.3 ± 0.3‰, respectively, and δ47 of 13.4 ± 0.4‰. Δ₄₇ in human breath was 0.76  ± 0.03‰, similar to that of ambient Pasadena air and higher than that of the car exhaust signature.     

Reconstruction of seasonal temperature variability in the tropical Pacific Ocean from the shell of the scallop, Comptopallium radula

We investigated the oxygen isotope composition (δ¹⁸O) of shell striae from juvenile Comptopallium radula (Mollusca; Pectinidae) specimens collected live in New Caledonia. Bottom-water temperature and salinity were monitored in-situ throughout the study period. External shell striae form with a 2-day periodicity in this scallop, making it possible to estimate the date of precipitation for each calcite sample collected along a growth transect. The oxygen isotope composition of shell calcite (δ¹⁸O_{shell calcite}) measured at almost weekly resolution on calcite accreted between August 2002 and July 2003 accurately tracks bottom-water temperatures. A new empirical paleotemperature equation for this scallop species relates temperature and δ¹⁸O_{shell calcite}:

t(°C)=20.00(±0.61)-3.66(±0.39)×(δ¹⁸O_{shell calcite VPDB}-δ¹⁸O_{water VSMOW})     

A model for oxygen and sulfur isotope fractionation in sulfate during bacterial sulfate reduction processes

We present a model of bacterial sulfate reduction that includes equations describing the fractionation relationship between the sulfur and the oxygen isotope composition of residual sulfate (δ³⁴S_{SO4_residual}, δ¹⁸O_{SO4_residual}) and the amount of residual sulfate. The model is based exclusively on oxygen isotope exchange between cell-internal sulfur compounds and ambient water as the dominating mechanism controlling oxygen isotope fractionation processes. We show that our model explains δ³⁴S_{SO4_residual} vs. δ¹⁸O_{SO4_residual} patterns observed from natural environments and from laboratory experiments, whereas other models, favoring kinetic isotope fractionation processes as dominant process, fail to explain many (but not all) observed δ³⁴S_{SO4_residual} vs. δ¹⁸O_{SO4_residual} patterns. Moreover, we show that a “typical” δ³⁴S_{SO4_residual} vs. δ¹⁸O_{SO4_residual} slope does not exist. We postulate that measurements of δ³⁴S_{SO4_residual} and δ¹⁸O_{SO4_residual} can be used as a tool to determine cell-specific sulfate reduction rates, oxygen isotope exchange rates, and equilibrium oxygen isotope exchange factors. Data from culture experiments are used to determine the range of sulfur isotope fractionation factors in which a simplified set of equations can be used. Numerical examples demonstrate the application of the equations. We postulate that, during denitrification, the oxygen isotope effects in residual nitrate are also the result of oxygen isotope exchange with ambient water. Consequently, the equations for the relationship between δ³⁴S_{SO4_residual}, δ¹⁸O_{SO4_residual}, and the amount of residual sulfate could be modified and used to calculate the fractionation-relationship between δ¹⁵N_{NO3_residual}, δ¹⁸O_{NO3_residual}, and the amount of residual nitrate during denitrification.     

Deconvolving the δO seawater component from subseasonal coral δO and Sr/Ca at Rarotonga in the southwestern subtropical Pacific for the period 1726 to 1997

To reconstruct oceanographic variations in the subtropical South Pacific, 271-year long subseasonal time series of Sr/Ca and δ¹⁸O were generated from a coral growing at Rarotonga (21.5°S, 159.5°W). In this case, coral Sr/Ca appears to be an excellent proxy for sea surface temperature (SST) and coral δ¹⁸O is a function of both SST and seawater δ¹⁸O composition (δ¹⁸O_sw). Here, we focus on extracting the δ¹⁸O_sw signal from these proxy records. A method is presented assuming that coral Sr/Ca is solely a function of SST and that coral δ¹⁸O is a function of both SST and δ¹⁸O_sw. This method separates the effects of δ¹⁸O_sw from SST by breaking the instantaneous changes of coral δ¹⁸O into separate contributions by instantaneous SST and δ¹⁸O_sw changes, respectively. The results show that on average δ¹⁸O_sw at Rarotonga explains ∼39% of the variance in δ¹⁸O and that variations in SST explains the remaining ∼61% of δ¹⁸O variance. Reconstructed δ¹⁸O_sw shows systematic increases in summer months (December-February) consistent with the regional pattern of variations in precipitation and evaporation. The δ¹⁸O_sw also shows a positive linear correlation with satellite-derived estimated salinity for the period 1980 to 1997 (r = 0.72). This linear correlation between reconstructed δ¹⁸O_sw and salinity makes it possible to use the reconstructed δ¹⁸O_sw to estimate the past interannual and decadal salinity changes in this region. Comparisons of coral δ¹⁸O and δ¹⁸O_sw at Rarotonga with the Pacific decadal oscillation index suggest that the decadal and interdecadal salinity and SST variability at Rarotonga appears to be related to basin-scale decadal variability in the Pacific.     

Theoretical investigation of iron isotope fractionation between Fe(H2O)6 and Fe(H2O)6: Implications for iron stable isotope geochemistry

The magnitude of equilibrium iron isotope fractionation between Fe(H₂O)₆³⁺ and Fe(H₂O)₆²⁺ is calculated using density functional theory (DFT) and compared to prior theoretical and experimental results. DFT is a quantum chemical approach that permits a priori estimation of all vibrational modes and frequencies of these complexes and the effects of isotopic substitution. This information is used to calculate reduced partition function ratios of the complexes (10³ · ln(β)), and hence, the equilibrium isotope fractionation factor (10³ · ln(α)). Solvent effects are considered using the polarization continuum model (PCM). DFT calculations predict fractionations of several per mil in ⁵⁶Fe/⁵⁴Fe favoring partitioning of heavy isotopes in the ferric complex. Quantitatively, 10³ · ln(α) predicted at 22°C, ∼ 3 ‰, agrees with experimental determinations but is roughly half the size predicted by prior theoretical results using the Modified Urey-Bradley Force Field (MUBFF) model. Similar comparisons are seen at other temperatures. MUBFF makes a number of simplifying assumptions about molecular geometry and requires as input IR spectroscopic data. The difference between DFT and MUBFF results is primarily due to the difference between the DFT-predicted frequency for the ν₄ mode (O-Fe-O deformation) of Fe(H₂O)₆³⁺ and spectroscopic determinations of this frequency used as input for MUBFF models (185-190 cm⁻¹ vs. 304 cm⁻¹, respectively). Hence, DFT-PCM estimates of 10³ · ln(β) for this complex are ∼ 20% smaller than MUBFF estimates. The DFT derived values can be used to refine predictions of equilibrium fractionation between ferric minerals and dissolved ferric iron, important for the interpretation of Fe isotope variations in ancient sediments. Our findings increase confidence in experimental determinations of the Fe(H₂O)₆³⁺ − Fe(H₂O)₆²⁺ fractionation factor and demonstrate the utility of DFT for applications in “heavy” stable isotope geochemistry.