OCO in Earth’s atmosphere

The chemistry and budgets of atmospheric gases are constrained by their bulk stable isotope compositions (e.g., δ¹³C values), which are based on mixing ratios of isotopologues containing one rare isotope (e.g., ¹⁶O¹³C¹⁶O). Atmospheric gases also have isotopologues containing two or more rare isotopes (e.g., ¹⁸O¹³C¹⁶O). These species have unique physical and chemical properties and could help constrain origins of atmospheric gases and expand the scope of stable isotope geochemistry generally. We present the first measurements of the abundance of ¹⁸O¹³C¹⁶O from natural and synthetic sources, discuss the factors influencing its natural distribution and, as an example of its applied use, demonstrate how its abundance constrains the sources of CO₂ in the Los Angeles basin. The concentration of ¹⁸O¹³C¹⁶O in air can be explained as a combination of ca. 1‰ enrichment (relative to the abundance expected if C and O isotopes are randomly distributed among all possible isotopologues) due to enhanced thermodynamic stability of this isotopologue during isotopic exchange with leaf and surface waters, ca. 0.1‰ depletion due to diffusion through leaf stomata, and subtle (ca. 0.05‰) dilution by ¹⁸O¹³C¹⁶O-poor anthropogenic CO₂. Some air samples are slightly (ca. 0.05‰) lower in ¹⁸O¹³C¹⁶O than can be explained by these factors alone. Our results suggest that ¹⁸O¹³C¹⁶O abundances should vary by up to ca. 0.2‰ with latitude and season, and might have measurable sensitivities to stomatal conductances of land plants. We suggest the greatest use of Δ₄₇ measurements will be to “leverage” interpretation of the δ¹⁸O of atmospheric CO₂.     

Isotopic fractionations associated with phosphoric acid digestion of carbonate minerals: Insights from first-principles theoretical modeling and clumped isotope measurements

Phosphoric acid digestion has been used for oxygen- and carbon-isotope analysis of carbonate minerals since 1950, and was recently established as a method for carbonate ‘clumped isotope’ analysis. The CO₂ recovered from this reaction has an oxygen isotope composition substantially different from reactant carbonate, by an amount that varies with temperature of reaction and carbonate chemistry. Here, we present a theoretical model of the kinetic isotope effects associated with phosphoric acid digestion of carbonates, based on structural arguments that the key step in the reaction is disproportionation of H₂CO₃ reaction intermediary. We test that model against previous experimental constraints on the magnitudes and temperature dependences of these oxygen isotope fractionations, and against new experimental determinations of the fractionation of ¹³C-¹⁸O-containing isotopologues (‘clumped’ isotopic species). Our model predicts that the isotope fractionations associated with phosphoric acid digestion of carbonates at 25 °C are 10.72‰, 0.220‰, 0.137‰, 0.593‰ for, respectively, ¹⁸O/¹⁶O ratios (1000 lnα^∗) and three indices that measure proportions of multiply-substituted isotopologues . We also predict that oxygen isotope fractionations follow the mass dependence exponent, λ of 0.5281 (where ). These predictions compare favorably to independent experimental constraints for phosphoric acid digestion of calcite, including our new data for fractionations of ¹³C-¹⁸O bonds (the measured change in Δ₄₇ = 0.23‰) during phosphoric acid digestion of calcite at 25 °C.We have also attempted to evaluate the effect of carbonate cation compositions on phosphoric acid digestion fractionations using cluster models in which disproportionating H₂CO₃ interacts with adjacent cations. These models underestimate the magnitude of isotope fractionations and so must be regarded as unsucsessful, but do reproduce the general trend of variations and temperature dependences of oxygen isotope acid digestion fractionations among different carbonate minerals. We suggest these results present a useful starting point for future, more sophisticated models of the reacting carbonate/acid interface. Examinations of these theoretical predictions and available experimental data suggest cation radius is the most important factor governing the variations of isotope fractionation among different carbonate minerals. We predict a negative correlation between acid digestion fractionation of oxygen isotopes and of ¹³C-¹⁸O doubly-substituted isotopologues, and use this relationship to estimate the acid digestion fractionation of for different carbonate minerals. Combined with previous theoretical evaluations of ¹³C-¹⁸O clumping effects in carbonate minerals, this enables us to predict the temperature calibration relationship for different carbonate clumped isotope thermometers (witherite, calcite, aragonite, dolomite and magnesite), and to compare these predictions with available experimental determinations. The success of our models in capturing several of the features of isotope fractionation during acid digestion supports our hypothesis that phosphoric acid digestion of carbonate minerals involves disproportionation of transition state structures containing H₂CO₃.     

Negative δO values in Allan Hills 84001 carbonate: Possible evidence for water precipitation on mars

The Martian meteorite ALH84001 contains ∼1% by weight of carbonate formed by secondary processes on the Martian surface or in the shallow subsurface. The major form of this carbonate is chemically and isotopically zoned rosettes which have been well documented elsewhere. This study concentrates upon carbonate regions ∼200 μm across which possess previously unobserved magnesium rich inner cores, interpreted here as rosette fragments, surrounded by a later stage cement containing rare Ca-rich carbonates (up to Ca₈₁Mg₀₇Fe₀₄Mn₀₇) intimately associated with feldspar. High spatial resolution ion probe analyses of Ca-rich carbonate surrounding rosette fragments have δ¹⁸O_V-SMOW values as low as −10‰. These values are not compatible with deposition from a global Martian atmosphere invoked to explain ALH84001 rosettes. The range of δ¹⁸O values are also incompatible with a fluid that has equilibrated with the Martian crust at high temperature or from remobilisation of carbonate of rosette isotopic composition. At Martian atmospheric temperatures, the small CO₂(gas)-CO₂(ice) fractionation makes meteoric CO₂ an unlikely source for −10‰ carbonates. In contrast, closed system Rayleigh fractionation of H₂O can generate δ¹⁸O_H2O −30‰, as observed at high latitudes on Earth. We suggest that atmospheric transport and precipitation of H₂O in a similar fashion to that on Earth provides a source of suitably ¹⁸O depleted water for generation of carbonate with δ¹⁸O_V-SMOW = −10‰.     

Determination of carbon fractionation factor between aqueous carbonate and CO2(g) in two-direction isotope equilibration

The experiments were conducted in the open CO₂ system to find out the equilibrium fractionation between the carbonate ion and CO₂(g). The existence of isotopic equilibrium was checked using the two-direction approach by passing the CO₂−N₂ gases with different δ¹³C compositions (− 1.5‰ and − 23‰) through the carbonate solution with δ¹³C = − 4.2‰. The Δ_{CO3T²⁻−CO2(g)} equilibrium fractionation is given as 6.03 ± 0.17‰ at 25 °C. Discussion is provided about the significance of carbonate complexing in determination of Δ_{CO3T²⁻−CO2(g)} and Δ_{HCO3T⁻−CO2(g)} fractionations. Finally, an isotope numerical model of flow and kinetics of hydration and dehydroxylation is built to predict the isotopic behaviour of the system with time.     

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.     

Defining an absolute reference frame for ‘clumped’ isotope studies of CO2

We present a revised approach for standardizing and reporting analyses of multiply substituted isotopologues of CO₂ (i.e., ‘clumped’ isotopic species, especially the mass-47 isotopologues). Our approach standardizes such data to an absolute reference frame based on theoretical predictions of the abundances of multiply-substituted isotopologues in gaseous CO₂ at thermodynamic equilibrium. This reference frame is preferred over an inter-laboratory calibration of carbonates because it enables all laboratories measuring mass 47 CO₂ to use a common scale that is tied directly to theoretical predictions of clumping in CO₂, regardless of the laboratory’s primary research field (carbonate thermometry or CO₂ biogeochemistry); it explicitly accounts for mass spectrometric artifacts rather than convolving (and potentially confusing) them with chemical fractionations associated with sample preparation; and it is based on a thermodynamic equilibrium that can be experimentally established in any suitably equipped laboratory using commonly available materials.By analyzing CO₂ gases that have been subjected to established laboratory procedures known to promote isotopic equilibrium (i.e., heated gases and water-equilibrated CO₂), and by reference to thermodynamic predictions of equilibrium isotopic distributions, it is possible to construct an empirical transfer function that is applicable to data with unknown clumped isotope signatures. This transfer function empirically accounts for the fragmentation and recombination reactions that occur in electron impact ionization sources and other mass spectrometric artifacts. We describe the protocol necessary to construct such a reference frame, the method for converting gases with unknown clumped isotope compositions to this reference frame, and suggest a protocol for ensuring that all reported isotopic compositions (e.g., Δ₄₇ values; [Eiler and Schauble, 2004] and [Eiler, 2007]) can be compared among different laboratories and instruments, independent of laboratory-specific analytical or methodological differences. We then discuss the use of intra-laboratory secondary reference frames (e.g., based on carbonate standards) that can be more easily used to track the evolution of each laboratory’s empirical transfer function. Finally, we show inter-laboratory reproducibility on the order of ±0.010 (1σ) for four carbonate standards, and present revised paleotemperature scales that should be used to convert carbonate clumped isotope signatures to temperature when using the absolute reference frame described here. Even when using the reference frame, small discrepancies remain between two previously published synthetic carbonate calibrations. We discuss possible reasons for these discrepancies, and highlight the need for additional low temperature (<15 °C) synthetic carbonate experiments.     

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.     

Effects of changing solution chemistry on Fe/Fe isotope fractionation in aqueous Fe-Cl solutions

The range in ⁵⁶Fe/⁵⁴Fe isotopic compositions measured in naturally occurring iron-bearing species is greater than 5‰. Both theoretical modeling and experimental studies of equilibrium isotopic fractionation among iron-bearing species have shown that significant fractionations can be caused by differences in oxidation state (i.e., redox effects in the environment) as well as by bond partner and coordination number (i.e., nonredox effects due to speciation).To test the relative effects of redox vs. nonredox attributes on total Fe equilibrium isotopic fractionation, we measured changes, both experimentally and theoretically, in the isotopic composition of an Fe²⁺-Fe³⁺-Cl-H₂O solution as the chlorinity was varied. We made use of the unique solubility of FeCl₄⁻ in immiscible diethyl ether to create a separate spectator phase against which changes in the aqueous phase could be quantified. Our experiments showed a reduction in the redox isotopic fractionation between Fe²⁺- and Fe³⁺-bearing species from 3.4‰ at [Cl⁻] = 1.5 M to 2.4‰ at [Cl⁻] = 5.0 M, due to changes in speciation in the Fe-Cl solution. This experimental design was also used to demonstrate the attainment of isotopic equilibrium between the two phases, using a ⁵⁴Fe spike.To better understand speciation effects on redox fractionation, we created four new sets of ab initio models of the ferrous chloride complexes used in the experiments. These were combined with corresponding ab initio models for the ferric chloride complexes from previous work. At 20 °C, 1000 ln β (β = ⁵⁶Fe/⁵⁴Fe reduced partition function ratio relative to a dissociated Fe atom) values range from 6.39‰ to 5.42‰ for Fe(H₂O)₆²⁺, 5.98‰ to 5.34‰ for FeCl(H₂O)₅⁺, and 5.91‰ to 4.86‰ for FeCl₂(H₂O)₄, depending on the model. The theoretical models predict ferric-ferrous fractionation about half as large (depending on model) as the experimental results.Our results show (1) oxidation state is likely to be the dominant factor controlling equilibrium Fe isotope fractionation in solution and (2) nonredox attributes (such as ligands present in the aqueous solution, speciation and relative abundances, and ionic strength of the solution) can also have significant effects. Changes in the isotopic composition of an Fe-bearing solution will influence the resultant Fe isotopic signature of any precipitates.     

Seasonal and diurnal variations of COO in air: Initial observations from Pasadena, CA

The budget of atmospheric CO₂ is widely studied using records of temporal and spatial variations of concentrations, δ¹³C and δ¹⁸O values. However, the number and diversity of sources and sinks prevents these alone from fully constraining the budget. Molecules containing two rare isotopes can serve as additional tracers and potentially provide additional, independent, constraints. We present data documenting seasonal and diurnal variations of CO₂ having a mass of 47 u (mostly ¹³C¹⁸O¹⁶O) in air from Pasadena, CA. We report these data using the ‘mass 47 anomaly’ (Δ₄₇), which is defined as the deviation of R⁴⁷(=[47]/[44]) from that expected for a random distribution of isotopologues. Between February 2004 and December 2005, Δ₄₇ showed a seasonal pattern that differed significantly from that expected based on thermodynamic equilibrium. During the year 2004 Δ₄₇ was 0.76‰ in winter, increased to 0.87‰ in summer and gradually decreased through the autumn to 0.81‰ at the end of the year. Δ₄₇ then increased again through the winter and spring of 2005 to 0.97‰ in summer followed by a decrease to 0.88‰ at the end of 2005. The seasonal variations cannot be accounted for by variations in the relative contribution of local fossil fuel sources. Diurnal variations were the combined effect of both fuel combustion and respiration having Δ₄₇ values of 0.41‰ and ca. 0.77‰, respectively. The seasonal cycle may be interpreted as a competition between low Δ₄₇ values in respiration and higher Δ₄₇ values resulting from CO₂-water exchange in photosynthesis.     

C-O bonds in carbonate minerals: A new kind of paleothermometer

The abundance of the doubly substituted CO₂ isotopologue, ¹³C¹⁸O¹⁶O, in CO₂ produced by phosphoric acid digestion of synthetic, inorganic calcite and natural, biogenic aragonite is proportional to the concentration of ¹³C-¹⁸O bonds in reactant carbonate, and the concentration of these bonds is a function of the temperature of carbonate growth. This proportionality can be described between 1 and 50 °C by the function: Δ₄₇ = 0.0592 · 10⁶ · T⁻² − 0.02, where Δ₄₇ is the enrichment, in per mil, of ¹³C¹⁸O¹⁶O in CO₂ relative to the amount expected for a stochastic (random) distribution of isotopes among all CO₂ isotopologues, and T is the temperature in Kelvin. This relationship can be used for a new kind of carbonate paleothermometry, where the temperature-dependent property of interest is the state of ordering of ¹³C and ¹⁸O in the carbonate lattice (i.e., bound together vs. separated into different CO₃²⁻ units), and not the bulk δ¹⁸O or δ¹³C values. Current analytical methods limit precision of this thermometer to ca. ± 2 °C, 1σ. A key feature of this thermometer is that it is thermodynamically based, like the traditional carbonate-water paleothermometer, and so is suitable for interpolation and even modest extrapolation, yet is rigorously independent of the δ¹⁸O of water and δ¹³C of DIC from which carbonate grew. Thus, this technique can be applied to parts of the geological record where the stable isotope compositions of waters are unknown. Moreover, simultaneous determinations of Δ₄₇ and δ¹⁸O for carbonates will constrain the δ¹⁸O of water from which they grew.     

Oxygen isotopic compositions of IVA iron meteorites: implications for the thermal evolution derived from in situ ultraviolet laser microprobe analyses

Oxygen isotopic compositions of silicate inclusions in IVA iron meteorites have been measured with an in situ UV laser microprobe technique. The homogeneity of oxygen isotopic compositions within and among individual mineral grains has also been examined. Oxygen isotope fractionations between coexisting mineral pairs were utilized in oxygen isotope thermometry. Our measured Δ¹⁷O values, ranging from 0.97 to 1.25‰, are characteristic of a single reservoir and fully confirm the oxygen isotopic similarity between IVA irons and L/LL chondrites. Steinbach and São João Nepomuceno, containing inclusions of two silicate minerals in mutual contact, exhibit a mass-dependent fractionation of ¹⁸O/¹⁶O between tridymite and bronzite with apparent oxygen isotopic heterogeneity. The SiO₂-bearing member, Gibeon, gives homogeneous oxygen isotopic compositions without detectable fractionation of ¹⁸O/¹⁶O between tridymite and quartz. Oxygen isotope equilibrium temperatures are estimated for coexisting tridymite and bronzite in the same sample slabs or clusters in Steinbach and São João Nepomuceno. The fractionations of ¹⁸O/¹⁶O between bronzite and tridymite range from 1.6 to 2.3‰ in different sample slabs or clusters. On the basis of the closure temperature concept, cooling rates are estimated at approximately 20 to 1000°C/Myr between 800 and 1000°C, a range of temperatures not accessible to other cooling rate methods. Using the Fast Grain Boundary diffusion model, we have demonstrated that significant oxygen heterogeneity both in tridymite and bronzite is probably due to isotope exchange during cooling between minerals with various grain sizes and mineral abundances in different regions of the samples. The new estimates of cooling rate by oxygen isotope thermometry refine previous cooling curves of IVA irons and support the breakup-reassembly model for the IVA parent body.     

Oxygen isotope signature of Pi regeneration from organic compounds by phosphomonoesterases and photooxidation

Phosphomonoesters are an important source of dissolved inorganic orthophosphate (PO₄ or P_i), the preferred form of P utilized by microbiota and aquatic plants in marine and freshwater ecosystems. Two enzymes involved in phosphomonoester metabolism and commonly detected in natural waters (alkaline phosphatase and 5′-nucleotidase) have been studied to determine the oxygen isotope signature of P_i-regeneration from phosphomonoesters by enzymatic degradation. Oxygen (O) isotope ratios of water and P_i released from phosphomonoesters during enzyme hydrolysis experiments demonstrate that released P_i incorporates one oxygen atom from water. The isotopic fractionation between this incorporated water O and ambient water O is −30 (±8)‰ for alkaline phosphatase and −10 (±1)‰ for 5′-nucleotidase, with very weak dependence on temperature. The result of these enzyme-specific isotopic fractionations at one of the four O sites in PO₄ is that the δ¹⁸O value of P_i regenerated by 5′-nucleotidase is 5‰ higher than P_i regenerated by alkaline phosphatase from the same phosphomonoester substrate. The δ¹⁸O value of regenerated P_i also reflects inheritance of 75% of O from the phosphomonoester substrate, thus making the δ¹⁸O of phosphomonoester-derived P_i a potential tracer of organophosphorous compound sources. Phosphomonoesterase-regenerated P_i has a distinct phosphate oxygen isotope signature that is different and distinguishable from that of biologically recycled and subsequently equilibrated P_i and P_i regenerated from photooxidation of organic matter. The δ¹⁸O value of regenerated P_i will correlate positively with the δ¹⁸O value of bulk water and the fractionation, α, between regenerated P_i and water, α_{regen Pi-water}, should converge toward equilibrium α_Pi-water values with increased biological cycling of P_i.     

O and S isotopic composition of dissolved and attached oxidation products of pyrite by Acidithiobacillus ferrooxidans: Comparison with abiotic oxidations

The acidophilic iron-oxidizing bacterium, Acidithiobacillus ferrooxidans, plays a part in the pyrite oxidation process and has been widely studied in order to determine the kinetics of the reactions and the isotopic composition of dissolved product sulphates, but the details of the oxidation processes at the surface of pyrite are still poorly known. In this study, oxygen and sulphur isotopic compositions (δ¹⁸O and δ³⁴S) were analyzed for dissolved sulphates and water from experimental aerobic acidic (pH < 2) pyrite oxidation by A. ferrooxidans. The oxidation products attached to the pyrite surfaces were studied for their morphology (SEM), their chemistry (Raman spectroscopy) and for their δ¹⁸O (ion microprobe). They were compared to abiotically (Fe³⁺, H₂O₂, O₂) oxidized pyrite surface compounds in order to constrain the oxidation pathways and to look for the existence of potential biosignatures for this system.The pyrite dissolution evolved from non-stoichiometric (during the first days) to stoichiometric (with increasing time) resulting in dissolved sulphates having distinct δ¹⁸O (e.g. +11.0‰ and −2.0‰, respectively) and δ³⁴S (+4.5‰ and +2.8‰, respectively) values. The “oxidation layer” at the surface of pyrite is complex and made of iron oxides, sulphate, polysulphide, elemental sulphur and polythionates. Bio- and Fe³⁺-oxidation favour the development of monophased micrometric bumps made of hematite or sulphate while other abiotic oxidation processes result in more variable oxidation products. The δ¹⁸O of these oxidation products at the surface of oxidized pyrites are strongly variable (from ≈−40‰ to ≈+30‰) for all experiments.Isotopic fractionation between sulphates and pyrite, Δ³⁴S_SO4-pyrite, is equal to −1.3‰ and +0.4‰ for sulphates formed by stoichiometric and non-stoichiometric processes, respectively. These two values likely reflect either a S-S or a Fe-S bond breaking process. The Δ¹⁸O_SO4-H2O and Δ¹⁸O_SO4-O2 are estimated to be ≈+16‰ and ≈−25‰, respectively. These values are higher than previously published data and may reflect biological effects. The large δ¹⁸O heterogeneity measured at the surfaces of oxidized pyrites, whatever the oxidant, may be related (i) to the existence of local surface environments isolated from the solution in which the oxidation processes are different and (ii) to the stabilization at the pyrite surface of reaction intermediates that are not in isotopic equilibrium with the solution. Though the oxygen isotopic composition of surface oxidation products cannot be taken as a direct biosignature, the combined morphological, chemical and isotopic characterization of the surfaces of oxidized pyrites may furnish clues about a biological activity on a mineral surface.     

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.     

Density functional theory predictions of equilibrium isotope fractionation of iron due to redox changes and organic complexation

We present molecular orbital/density functional theory (MO/DFT) calculations that predict a greater isotopic fractionation in redox reactions than in reactions involving ligand exchange. The predicted fractionation factors, reported as 1000·ln(^56-54α), associated with equilibrium between Fe-organic and Fe-H₂O species were <1.6‰ in vacuo and <1.2‰ in solution when the oxidation state of the system was held constant. These fractionation factors were significantly smaller than those predicted for equilibrium between different oxidation states of Fe, for which 1000·ln(^56-54α) was >2.7‰ in vacuo and >2.2‰ in solution when the bound ligands were unchanged. The predicted ⁵⁶Fe/⁵⁴Fe ratio was greater in complexes containing Fe³⁺ and in complexes with shorter Fe-O bond lengths; both of these trends follow previous theoretical results. Our predictions also agree with previous experimental measurements that suggest that the largest biological fractionations will be associated with processes that change the oxidation state of Fe, and that identification of biologically controlled Fe isotope fractionation may be difficult when abiotic redox fractionations are present in the system. The models studied here also have important implications for future theoretical isotope calculations, because we have discovered the necessity of using vibrational frequencies instead of reduced masses when predicting reduced partition functions in aqueous-phase species.     

Carbon isotopic fractionation during a mesocosm bloom experiment dominated by Emiliania huxleyi: Effects of CO2 concentration and primary production

We investigated the effect of CO₂ and primary production on the carbon isotopic fractionation of alkenones and particulate organic matter (POC) during a natural phytoplankton bloom dominated by the coccolithophore Emiliania huxleyi. In nine semi-closed mesocosms (∼11 m³ each), three different CO₂ partial pressures (pCO₂) in triplicate represented glacial (∼180 ppmv CO₂), present (∼380 ppmv CO₂), and year 2100 (∼710 ppmv CO₂) CO₂ conditions. The largest shift in alkenone isotopic composition (4-5‰) occurred during the exponential growth phase, regardless of the CO₂ concentration in the respective treatment. Despite the difference of ∼500 ppmv, the influence of pCO₂ on isotopic fractionation was marginal (1-2‰). During the stationary phase, E. huxleyi continued to produce alkenones, accumulating cellular concentrations almost four times higher than those of exponentially dividing cells. Our isotope data indicate that, while alkenone production was maintained, the interaction of carbon source and cellular uptake dynamics by E. huxleyi reached a steady state. During stationary phase, we further observed a remarkable increase in the difference between δ¹³C of bulk organic matter and of alkenones spanning 7-12‰. We suggest that this phenomenon is caused mainly by a combination of extracellular release of ¹³C-enriched polysaccharides and subsequent particle aggregation induced by the production of transparent exopolymer particles (TEP).     

Oxygen and carbon isotopic systematics of aragonite speleothems and water in Furong Cave, Chongqing, China

To understand oxygen and carbon stable isotopic characteristics of aragonite stalagmites and evaluate their applicability to paleoclimate, the isotopic compositions of active and fossil aragonite speleothems and water samples from an in situ multi-year (October 2005-July 2010) monitoring program in Furong Cave located in Chongqing of China have been examined. The observations during October 2005-June 2007 show that the meteoric water is well mixed in the overlying 300-500-m bedrock aquifer, reflected by relatively constant δ¹⁸O, ±0.11-0.14‰ (1σ), of drip waters in the cave, which represents the annual status of rainfall water. Active cave aragonite speleothems are at oxygen isotopic equilibrium with drip water and their δ¹⁸O values capture the surface-water oxygen isotopic signal. Aragonite-to-calcite transformation since the last glaciation is not noticeable in Furong stalagmites. Our multi-year field experiment approves that aragonite stalagmite δ¹⁸O records in this cave are suitable for paleoclimate reconstruction. With high U, 0.5-7.2 ppm, and low Th, 20-1270 ppt, the Furong aragonite stalagmites provide very precise chronology (as good as ±20s yrs (2σ)) of the climatic variations since the last deglaciation. The synchroneity of Chinese stalagmite δ¹⁸O records at the transition into the Bølling-Allerød (t-BA) and the Younger Dryas from Furong, Hulu and Dongge Caves supports the fidelity of the reconstructed East Asian monsoon evolution. However, the Furong record shows that the cold Older Dryas (OD) occurred at 14.0 thousand years ago, agreeing with Greenland ice core δ¹⁸O records but ∼200 yrs younger than that in the Hulu record. The OD age discrepancy between Chinese caves can be attributable to different regionally climatic/environmental conditions or chronological uncertainty of stalagmite proxy records, which is limited by changes in growth rate and subsampling intervals in absolute dating. Seasonal dissolved inorganic carbon δ¹³C variations of 2-3‰ in the drip water and 5-7‰ in the pool and spring waters are likely attributed to variable degrees of CO₂ degassing in winter and summer. The variable δ¹³C values of active deposits from −11‰ to 0‰ could be caused by kinetically mediated CO₂ degassing processes. The complicated nature of pre-deposition kinetic isotopic fractionation processes for carbon isotopes in speleothems at Furong Cave require further study before they can be interpreted in a paleoclimatic or paleoenvironmental context.     

Kinetic nitrogen isotope fractionation associated with thermal decomposition of NH3: Experimental results and potential applications to trace the origin of N2 in natural gas and hydrothermal systems

Ammonia (NH₃) is the major intermediate phase in the pathway of nitrogen (N) transfer from the fixed N phases (e.g., in crustal material) to free N₂ (e.g., in natural gas reservoirs and volcanic gases). Yet the N isotopic behavior during these N-cycling processes remains poorly known. In an attempt to contribute to the understanding of N cycling using N isotopes, we carried out laboratory experiments to investigate the N isotopic effect associated with thermal decomposition of ammonia (2NH₃ → N₂ + 3H₂). Pure NH₃ (with initial δ¹⁵N_NH3 of ∼ −2‰, relative to air standard) was sealed into quartz tubes and thermally decomposed at 600, 700 or 800 °C from 2 hours to 500 days. With the progress of the reaction, the δ¹⁵N of the remaining NH₃ and the accumulated N₂ increased from −2 to +35‰ and from −20 to −2‰, respectively. The differences of the N-isotope fractionations at the three temperatures are not significant. Modeling using the Rayleigh distillation model yielded similar kinetic N-isotope fractionation factors (α_N2-NH3) of 0.983 ± 0.002 for 600, 700 and 800 °C. Applied to geological settings, this significant isotope discrimination (∼17‰) associated with partial decomposition of NH₃/NH₄⁺ from crustal sources (δ¹⁵N_average ∼ +6.3‰) can produce mantle-like (i.e. ∼ −5‰) or even lower δ¹⁵N values of N₂. This may explain the large variation of δ¹⁵N (−20 to +30‰) of N₂ in natural gas reservoirs. It can also possibly explain the extreme ¹⁵N-depletion of N₂ in some volcanic gases. This possibility has to be carefully considered when using N isotopes to trace geological N cycling across subduction zones by analysis of volcanic N₂.     

Experimental determination of oxygen isotope fractionations between CO2 vapor and soda-melilite melt

We report results of experiments constraining oxygen isotope fractionations between CO₂ vapor and Na-rich melilitic melt at 1 bar and 1250 and 1400°C. The fractionation factor constrained by bracketed experiments, 1000^.lnα_{CO2-Na melilitic melt}, is 2.65±0.25 ‰ (±2σ; n=92) at 1250°C and 2.16±0.16 ‰ (2σ; n=16) at 1400°C. These values are independent of Na content over the range investigated (7.5 to 13.0 wt. % Na₂O). We combine these data with the known reduced partition function ratio of CO₂ to obtain an equation describing the reduced partition function ratio of Na-rich melilite melt as a function of temperature. We also fit previously measured CO₂-melt or -glass fractionations to obtain temperature-dependent reduced partition function ratios for all experimentally studied melts and glasses (including silica, rhyolite, albite, anorthite, Na-rich melilite, and basalt). The systematics of these data suggest that reduced partition function ratios of silicate melts can be approximated either by using the Garlick index (a measure of the polymerization of the melt) or by describing melts as mixtures of normative minerals or equivalent melt compositions. These systematics suggest oxygen isotope fractionation between basalt and olivine at 1300°C of approximately 0.4 to 0.5‰, consistent with most (but not all) basalt glass-olivine fractionations measured in terrestrial and lunar basalts.     

The isotopic composition of nitrate produced from nitrification in a hardwood forest floor

Dual isotopic analysis of nitrate (¹⁵N/¹⁴N and ¹⁸O/¹⁶O) is increasingly used to investigate the environmental impacts of human-induced elevated atmospheric nitrate deposition. In forested ecosystems, the nitrate found in surface water and groundwater can originate from two sources: (1) atmospheric deposition, and (2) nitrate produced from nitrification in forest soils (microbial nitrate). Application of the dual nitrate isotope technique for determining the relative importance of nitrate sources in forested catchments requires knowledge of the isotopic composition of microbial nitrate. We excluded precipitation inputs to three zero-tension lysimeters installed below the F-horizon (Oe) at the Turkey Lakes Watershed (TLW) in order to measure the isotopic composition of microbial nitrate produced in situ. To our knowledge, this is the first in situ study of the isotopic composition of microbial nitrate in forest soils. Over a 2-week period, nitrate produced by nitrification was periodically flushed to the lysimeters by watering the area with a nitrogen-free solution. Nitrate produced in the forest floor had δ¹⁸O values ranging from +3.1‰ to +10.1‰ with a mean of +5.2‰. These values were only slightly higher than from the expected value of +1.0‰ calculated for chemolithoautotrophic nitrification, which depends on the δ¹⁸O of available O₂ and H₂O. In addition to nitrate, we also collected soil gas to determine if soil respiration and O₂ diffusion affected soil gas δ¹⁸O-O₂, which is typically assumed to be identical to atmospheric O₂ (+23.5‰) when calculating microbial nitrate δ¹⁸O values. No significant difference in δ¹⁸O-O₂ from atmospheric O₂ was found in forest soils to a depth of 55 cm, and therefore ¹⁸O-enrichment of soil gas O₂ could not explain the modest enrichment of nitrate ¹⁸O. Evaporative ¹⁸O-enrichment of soil water available to nitrifiers in the forest floor is a plausible mechanism for slightly elevated nitrate δ¹⁸O values. However, the observed nitrate δ¹⁸O values could also be explained by a minor contribution of nitrate from heterotrophic nitrifiers. The δ¹⁵N of nitrate produced ranged from −10.4 to −7.3‰ and, as expected, was depleted in ¹⁵N relative to soil organic nitrogen. Microbial nitrate produced in the forest floor was also significantly depleted in ¹⁵N relative to microbial nitrate exported in groundwater and headwater streams at the TLW. We hypothesize that ¹⁵N-depleted forest floor nitrate is not detected in groundwaters largely because of: (1) the immobilization of forest floor nitrate in the mineral soil and (2) the mixing of the remaining forest floor nitrate with nitrate generated in the mineral soil, which is expected to have higher δ¹⁵N values. This study demonstrates that current methods of calculating a priori the δ¹⁸O of microbial nitrate provide a reasonable value for nitrate produced by nitrification at the TLW.