Influence of physiology and climate on δD of leaf wax n-alkanes from C3 and C4 grasses

We measured hydrogen isotope compositions (δD) of high-molecular-weight n-alkanes (C₂₇-C₃₃) from grasses grown in greenhouses and collected from the US Great Plains. In both cases, n-alkanes from C₄ grasses are enriched in D by more than 20‰ relative to those from C₃ grasses. The apparent enrichment factor (ε_C29-GW) between C₂₉n-alkane and greenhouse water is −165 ± 12‰ for C₃ grasses and −140 ± 15‰ for C₄ grasses. For samples from the Great Plains, δD values of C₂₉n-alkanes range from −280 to −136‰, with values for C₄ grasses ca. 21‰ more positive than those for C₃ grasses from the same site. Differences in C₃ and C₄ grass n-alkane δD values are consistent with the shorter interveinal distance in C₄ grass leaves, and greater back-diffusion of enriched water from stomata to veins, than in C₃ grass leaves. Great Plains’ grass n-alkane isotopic ratios largely reflect precipitation δD values. However, the offset or apparent fractionation between n-alkanes and precipitation is not uniform and varies with annual precipitation and relative humidity, suggesting climatic controls on lipid δD values. The dryer sites exhibit smaller absolute apparent fractionation indicative of D-enrichment of source waters through transpiration and/or soil evaporation. To explore the relationship between climate and n-alkane δD values, we develop three models. (1) The ‘direct analog’ model estimates δD_C29 values simply by applying the apparent enrichment factors, ε_C29-GW, observed in greenhouse grasses to precipitation δD values from the Great Plains. (2) The ‘leaf-water’ model uses a Craig-Gordon model to estimate transpirational D-enrichment for both greenhouse and field sites. The transpiration-corrected enrichment factors between C₂₉ and bulk leaf-water, ε_C29-GW, calculated from the greenhouse samples (−181‰ for C₃ and −157‰ for C₄) are applied to estimate δD_C29 values relative to modeled bulk leaf-water δD values. (3) The ‘soil- and leaf-water’ model estimates the combined effects of soil evaporation, modeled by analogy with a flow-through lake, and transpiration on δD_C29 values. Predictions improve with the addition of the explicit consideration of transpiration and soil evaporation, indicating that they are both important processes in determining plant lipid δD values. D-enrichment caused by these evaporative processes is controlled by relative humidity, suggesting that important climatic information is recorded in leaf wax n-alkane δD values. Calibration studies such as this one provide a baseline for future studies of plant-water-deuterium systematics and form the foundation for interpretation of plant wax hydrogen isotope ratios as a paleo-aridity proxy.     

Carbon and hydrogen isotopic compositions of products of open-system catalytic hydrogenation of CO2: Implications for abiogenic hydrocarbons in Earth’s crust

This paper reports the isotope effects in an open-system Fischer-Tropsch type (FTT) synthesis, with implications for the origin of natural abiogenic hydrocarbons. The starting form of carbon was CO₂, with carbon and hydrogen isotopic compositions measured for products of catalytic hydrogenation of CO₂ on iron and cobalt catalysts (FTCO₂-Fe and FTCO₂-Co) at 350 and 245 °C, respectively, and 10 MPa. The carbon isotopic composition of the resulting saturated hydrocarbons (alkanes) as a function of carbon number shows a positive trend for both FTCO₂-Fe and FTCO₂-Co, with a fractionation of 2-4‰ and 3-6‰ between CH₄ and C₂H₆ over the Fe and Co catalysts, respectively. The unsaturated hydrocarbons (alkenes) do not show any trend. A strong kinetic isotope fractionation (>40‰) occurred between CO₂ and CH₄ in both experiments. The hydrogen isotope fractionation between alkanes appeared to be similar to that found in natural (thermogenic and biogenic) gases, with enrichment in deuterium of longer hydrocarbon chains; the dominant H/D fractionation occurred between CH₄ and C₂H₆. Alkenes in the products of the FTCO₂-Fe reaction are enriched in deuterium (∼50‰) and do not show any trend versus carbon number. We suggest that other than FTT reactions or a simple mixing are responsible for the occurrence of the inverse isotopic trends in both δ¹³C and δD found in light hydrocarbons in some terrestrial environments and meteorites.     

Genesis of fumarolic emissions as inferred by isotope mass balances: CO2 and water at Vulcano Island, Italy

We have developed a quantitative model of CO₂ and H₂O isotopic mixing between magmatic and hydrothermal gases for the fumarolic emissions of the La Fossa crater (Vulcano Island, Italy). On the basis of isotope balance equations, the model takes into account the isotope equilibrium between H₂O and CO₂ and extends the recent model of chemical and energy two-end-member mixing by Nuccio et al. (1999). As a result, the H₂O and CO₂ content and the δD, δ¹⁸O, and δ¹³C isotope compositions for both magmatic and hydrothermal end-members have been assessed. Low contributions of meteoric steam, added at a shallow depth, have been also recognized and quantified in the fumaroles throughout the period from 1988 to 1998. Nonequilibrium oxygen isotope exchange also seems to be occurring between ascending gases and wall rocks along some fumarolic conduits.The δ¹³C_CO2 of the magmatic gases varies around −3 to 1‰ vs. Peedee belemnite (PDB), following a perfect synchronism with the variations of the CO₂ concentration in the magmatic gases. This suggests a process of isotope fractionation because of vapor exsolution caused by magma depressurization. The hydrogen isotopes in the magmatic gases (−1 to −‰ vs. standard mean ocean water [SMOW]), as well as the above δ¹³C_CO2 value, are coherent with a convergent tectonic setting of magma generation, where the local mantle is widely contaminated by fluids released from the subducted slab. Magma contamination in the crust probably amplifies this effect.The computed isotope composition of carbon and hydrogen in the hydrothermal vapors has been used to calculate the δD and δ¹³C of the entire hydrothermal system, including mixed H₂O-CO₂ vapor, liquid water, and dissolved carbon. We have computed values of about 10‰ vs. SMOW for water and −2 to −6.5‰ vs. PDB for CO₂. On these grounds, we think that Mediterranean marine water (δD_H2O ≈ 10‰) feeds the hydrothermal system. It infiltrates at depth throughout the local rocks, reaching oxygen isotope equilibrium at high temperatures. Interaction processes between magmatic gases and the evolving seawater also seem to occur, causing the dissolution of isotopically fractionated aqueous CO₂ and providing the source for hydrothermal carbon. These results have important implications concerning fluid circulation beneath Vulcano and address the more convenient routine of geochemical surveillance.     

Gas seepage from Tokamachi mud volcanoes, onshore Niigata Basin (Japan): Origin, post-genetic alterations and CH4-CO2 fluxes

Methane and CO₂ emissions from the two most active mud volcanoes in central Japan, Murono and Kamou (Tokamachi City, Niigata Basin), were measured in from both craters or vents (macro-seepage) and invisible exhalation from the soil (mini- and microseepage). Molecular and isotopic compositions of the released gases were also determined. Gas is thermogenic (δ¹³C_CH4 from −32.9‰ to −36.2‰), likely associated with oil, and enrichments of ¹³C in CO₂ (δ¹³C_CO2 up to +28.3‰) and propane (δ¹³C_C3H8 up to −8.6‰) suggest subsurface petroleum biodegradation. Gas source and post-genetic alteration processes did not change from 2004 to 2010. Methane flux ranged within the orders of magnitude of 10¹-10⁴ g m⁻² d⁻¹ in macro-seeps, and up to 446 g m⁻² d⁻¹ from diffuse seepage. Positive CH₄ fluxes from dry soil were widespread throughout the investigated areas. Total CH₄ emission from Murono and Kamou were estimated to be at least 20 and 3.7 ton a⁻¹, respectively, of which more than half was from invisible seepage surrounding the mud volcano vents. At the macro-seeps, CO₂ fluxes were directly proportional to CH₄ fluxes, and the volumetric ratios between CH₄ flux and CO₂ flux were similar to the compositional CH₄/CO₂ volume ratio. Macro-seep flux data, in addition to those of other 13 mud volcanoes, supported the hypothesis that molecular fractionation (increase of the “Bernard ratio” C₁/(C₂ + C₃)) is inversely proportional to gas migration fluxes. The CH₄ “emission factor” (total measured output divided by investigated seepage area) was similar to that derived in other mud volcanoes of the same size and activity. The updated global “emission-factor” data-set, now including 27 mud volcanoes from different countries, suggests that previous estimates of global CH₄ emission from mud volcanoes may be significantly underestimated.     

Hydrogen isotopic fractionation in lipid biosynthesis by H2-consuming Desulfobacterium autotrophicum

We report hydrogen isotopic fractionations between water and fatty acids of the sulfate-reducing bacterium Desulfobacterium autotrophicum. Pure cultures were grown in waters with deuterium (D) contents that were systematically varied near the level of natural abundance (−37‰ ? δD ? 993‰). H₂ of constant hydrogen isotope (D/H) ratio was supplied to the cultures. The D/H ratios of water, H₂, and specific fatty acids were measured by isotope-ratio mass spectrometry. The results demonstrate that D. autotrophicum catalyzes hydrogen isotopic exchange between water and H₂, and this reaction is conclusively shown to approach isotopic equilibrium. In addition, variation in the D/H ratio of growth water accounts for all variation in the hydrogen isotopic composition of fatty acids. The D/H ratios of fatty acids from cultures grown on H₂/CO₂ are compared with those from a separate set of cultures grown on D-enriched formate, an alternative electron donor. This comparison rules out H₂ as a significant source of fatty acid hydrogen. Grown on either H₂/CO₂ or formate, D. autotrophicum produces fatty acids in which all hydrogen originates from water. For specific fatty acids, biosynthetic fractionation factors are mostly in the range 0.60 ? α_FA-water ? 0.70; the 18:0 fatty acid exhibits a lower fractionation factor of 0.52. The data show that α_FA-water generally increases with length of the carbon chain from C₁₄ to C₁₇ among both saturated and unsaturated fatty acids. These results indicate a net fractionation associated with fatty acid biosynthesis in D. autotrophicum that is slightly smaller than in another H₂-consuming bacterium (Sporomusa sp.), but much greater than in most photoautotrophs.     

Hydrogen isotopes in volcanic plumes: Tracers for remote temperature sensing of fumaroles

In high-temperature volcanic fumaroles (>400 °C), the isotopic composition of molecular hydrogen (H₂) reaches equilibrium with that of the fumarolic H₂O. In this study, we used this hydrogen isotope exchange equilibrium of fumarolic H₂ as a tracer for the remote temperature at volcanic fumaroles. In this remote sensing, we deduced the hydrogen isotopic composition (δD value) of fumarolic H₂ from those in the volcanic plume. To ascertain that we can estimate the δD value of fumarolic H₂ from those in a volcanic plume, we estimated the values in three fumaroles with outlet temperatures of 630 °C (Tarumae), 203 °C (Kuju), and 107 °C (E-san). For this we measured the concentration and δD value of H₂ in each volcanic plume, along with those determined directly at each fumarole. The average and maximum mixing ratios of fumarolic H₂ within a plume’s total H₂ were 97% and 99% (at Tarumae), 89% and 96% (at Kuju), and 97% and 99% (at E-san). We found a linear relationship between the depletion in the δD values of H₂, with the reciprocal of H₂ concentration. Furthermore, the estimated end-member δD value for each H₂-enriched component (−260 ± 30‰ vs. VSMOW in Tarumae, −509 ± 23‰ in Kuju, and −437 ± 14‰ in E-san) coincided well with those observed at each fumarole (−247.0 ± 0.6‰ in Tarumae, −527.7 ± 10.1‰ in Kuju, and −432.1 ± 2.5‰ in E-san). Moreover, the calculated isotopic temperatures at the fumaroles agreed to within 20 °C with the observed outlet temperature at Tarumae and Kuju. We deduced that the δD value of the fumarolic H₂ was quenched within the volcanic plume. This enabled us to remotely estimate these in the fumarole, and thus the outlet temperature of fumaroles, at least for those having the outlet temperatures more than 400 °C. By applying this methodology to the volcanic plume emitted from the Crater 1 of Mt. Naka-dake (the volcano Aso) where direct measurement on fumaroles was impractical, we estimated that the δD value of the fumarolic H₂ to be −172 ± 16‰ and the outlet temperature to be 868 ± 97 °C. The remote temperature sensing using hydrogen isotopes developed in this study is widely applicable to many volcanic systems.     

D/H ratios of atmospheric H2 in urban air: results using new methods for analysis of nano-molar H2 samples

We present the results of a study of the concentration and D/H ratio of molecular hydrogen from air in the Los Angeles Basin and adjacent San Gabriel Mountains. These data define a mixing relationship in dimensions of D/H ratio vs. 1/(H₂) which constrains the δD_VSMOW of unpolluted winter air in this region to be ca. +100 to +125 ‰ and that of urban H₂ sources to be ca. −270 ‰. This study makes use of a new method for measuring the deuterium content of molecular hydrogen in small samples (∼100 to 500 cc) of air, which we describe in detail. The method consists of an off-line separation of H₂ followed by introduction to the mass spectrometer in a continuous flow of He. Off-line, all components of an atmospheric gas sample, with the exception of He, H₂, and Ne are condensed by exposure to a cold-trap held at 30 Kelvin. This separation is followed by cryo-transfer of non-condensable gases to a small volume molecular sieve finger, with assist from a mercury piston pump. At the mass spectrometer, the sample is put in line with a continuous flow of He where it is focused on to an additional column of molecular sieve before subsequent introduction into the ion source. Analyses of DH/H₂ ratio have accuracy and precision of ±4 to 7 per mil. Comparison of sample peak area to peak areas of standards of known size allows for determination of H₂ concentration with accuracy and precision of ∼±5%, relative. The method reduces sample size and processing time by several orders of magnitude compared to previous methods, allowing for sampling at proportionately higher spatial and temporal resolution.     

Exchange and fractionation of Mg isotopes between epsomite and saturated MgSO4 solution

The equilibrium Mg isotope fractionation factor between epsomite and aqueous MgSO₄ solution has been measured using the three isotope method in recrystallization experiments conducted at 7, 20, and 40 °C. Complete or near-complete isotopic exchange was achieved within 14 days in all experiments. The Mg isotope exchange rate between epsomite and MgSO₄ solution is dependent on the temperature, epsomite seed crystal grain size, and experimental agitation method. The Mg isotope fractionation factors (Δ²⁶Mg_eps-sol) at 7, 20, and 40 °C are 0.63 ± 0.07‰, 0.58 ± 0.16‰, and 0.56 ± 0.03‰, respectively. These values are indistinguishable within error, indicating that the Mg isotope composition of epsomite is relatively insensitive to temperature. The magnitude of the isotope fractionation factor (Δ²⁶Mg_eps-sol = ca. 0.6‰ between 7 and 40 °C) indicates that significant Mg isotope variations can be produced in evaporite sequences, and Mg isotopes may therefore, constrain the degree of closed-system behavior, paleo-humidity, and hydrological history of evaporative environments.     

Effects of speciation on equilibrium fractionations and rates of oxygen isotope exchange between (PO4)aq and H2O

The effects of phosphate speciation on both rates of isotopic exchange and oxygen isotope equilibrium fractionation factors between aqueous phosphate and water were examined over the temperature range 70 to 180°C. Exchange between phosphate and water is much faster at low pH than at high pH, an observation that is similar to what has been observed in the analogous sulfate-water system. Oxygen isotope fractionations between protonated species like H₃PO₄ and H₂PO₄⁻ that are dominant at relatively low pH and species like PO₄³⁻ and ion pairs like KHPO₄⁻ that are dominant at relatively high pH, range between 5 and 8‰ at the temperatures of the experiments. In aqueous phosphate systems at equilibrium, ¹⁸O/¹⁶O ratios increase with increasing degree of protonation of phosphate. This effect can be explained in part by the relative magnitudes of the dissociation constants of the protonated species. Under equilibrium conditions, carbonate in solution or in solid phases concentrates ¹⁸O relative to orthophosphate in solution or in solid phases at all temperatures, supporting the traditional view that biogenic phosphate is precipitated in near oxygen isotope equilibrium with body/ambient aqueous fluids with no attendant vital effects.     

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₂.     

Different types of methane monooxygenases produce similar carbon and hydrogen isotope fractionation patterns during methane oxidation

We determined the stable carbon and hydrogen isotope fractionation factors for methane oxidation under oxic conditions using strains with known degradation pathways. The aerobic oxidation of methane can be initiated by two different forms of enzymes known as methane monooxygenases (MMO). The expression of these enzymes is type-specific and dependent upon the adjusted copper concentration in the medium (or environment). In this study, the expression of either the soluble MMO or the particulate MMO was supported by adjusting the copper concentrations in the growth medium. Taxonomically different aerobic methanotrophic strains, mainly belonging to the alpha- and gamma- classes of Proteobacteria, produced methane isotope enrichment factors (ε_bulk) ranging from −14.8 to −27.9‰ for carbon, and from −110.0 to −231.5‰ for hydrogen. The ratio of hydrogen versus carbon discrimination (Λ = (α_H⁻¹ − 1)/(α_C⁻¹ − 1) ≈ Δ(δ²H)/Δ(δ¹³C)) were similar for all tested cultures, and are also identical to values calculated from previously published enrichment factors for aerobic and anaerobic methane degradation. In contrast, Λ-values for the abiotic oxidation of methane with OH radicals (this process is considered as the main removal process for methane from the atmosphere) were significantly higher than the values derived from biotic oxidation. Due to the low variability of microbial methane isotope fractionation patterns, we propose that combined carbon and hydrogen isotope fractionation analyses can be used to monitor and assess the occurrence of microbial methane oxidation in marine or terrestrial environments. However, it is not possible to distinguish distinct aerobic or anaerobic methane-oxidation pathways by this approach.     

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.     

Sulfur isotope fractionation during growth of sulfate-reducing bacteria on various carbon sources

Stable sulfur isotope fractionation during microbial sulfate reduction is a potential tool to estimate sulfate reduction rates at field sites. However, little is known about the influence of the utilized carbon source on the magnitude of sulfur isotope fractionation. To investigate this effect, both a pure culture (strain PRTOL1) and enrichment cultures from a petroleum hydrocarbon (PHC)-contaminated aquifer were used and grown in batch cultures on various carbon sources with an initial sulfate concentration of 1 mmol/L. As sole carbon sources the PHC components naphthalene, 1,3,5-trimethylbenzene, and heating oil (enrichment culture) and the organic acids acetate, pyruvate, benzoate, and 3-phenylpropionate (enrichment culture and PRTOL1) were used. Sulfate reduction rates of all cultures ranged from 6 ± 1 nmol cm⁻³ d⁻¹ (enrichment culture grown on 1,3,5-trimethylbenzene) to 280 ± 6 nmol cm⁻³ d⁻¹ (enrichment culture grown on pyruvate). Cell-specific sulfate reduction rates ranged from 1.1 × 10⁻¹⁴ mol cell⁻¹ d⁻¹ (PRTOL1 grown on pyruvate) to 1.5 × 10⁻¹³ mol cell⁻¹ d⁻¹ (PRTOL1 grown on acetate). Sulfur isotope enrichment factors (ε) for the enrichment culture ranged from 16.1‰ (3-phenylpropionate) to 34.5‰ (1,3,5-trimethylbenzene) and for PRTOL1 from 30.0‰ (benzoate) to 36.0‰ (pyruvate). Cultures of PRTOL1 always showed higher ε values than the enrichment culture when grown on the same carbon source due to culture-specific properties. Higher ε values were obtained when the enrichment culture was grown on PHC components than on organic acids. No relationship between ε values and cell-specific sulfate reduction rate existed when all data were combined. When comparing the magnitude of ε values determined in this laboratory study with ε values measured at contaminated and uncontaminated field sites, it becomes evident that a multitude of factors influences ε values at field sites and complicates their interpretation. The results of this study help us assess some of the general parameters that govern the magnitude of ε in sulfate-reducing environments.     

Carbon and hydrogen isotope fractionation by moderately thermophilic methanogens

A series of laboratory studies were conducted to increase understanding of stable carbon (¹³C/¹²C) and hydrogen (D/H) isotope fractionation arising from methanogenesis by moderately thermophilic acetate- and hydrogen-consuming methanogens. Studies of the aceticlastic reaction were conducted with two closely related strains of Methanosaeta thermophila. Results demonstrate a carbon isotope fractionation of only 7‰ (α = 1.007) between the methyl position of acetate and the resulting methane. Methane formed by this process is enriched in ¹³C when compared with other natural sources of methane; the magnitude of this isotope effect raises the possibility that methane produced at elevated temperature by the aceticlastic reaction could be mistaken for thermogenic methane based on carbon isotopic content. Studies of H₂/CO₂ methanogenesis were conducted with Methanothermobacter marburgensis. The fractionation of carbon isotopes between CO₂ and CH₄ was found to range from 22 to 58‰ (1.023 ≤ α ≤ 1.064). Greater fractionation was associated with low levels of molecular hydrogen and steady-state metabolism. The fractionation of hydrogen isotopes between source H₂O and CH₄ was found to range from 127 to 275‰ (1.16 ≤ α ≤ 1.43). Fractionation was dependent on growth phase with greater fractionation associated with later growth stages. The maximum observed fractionation factor was 1.43, independent of the δD-H₂ supplied to the culture. Fractionation was positively correlated with temperature and/or metabolic rate. Results demonstrate significant variability in both hydrogen and carbon isotope fractionation during methanogenesis from H₂/CO₂. The relatively small fractionation associated with deuterium during H₂/CO₂ methanogenesis provides an explanation for the relatively enriched deuterium content of biogenic natural gas originating from a variety of thermal environments. Results from these experiments are used to develop a hypothesis that differential reversibility in the enzymatic steps of the H₂/CO₂ pathway gives rise to variability in the observed carbon isotope fractionation. Results are further used to constrain the overall efficiency of electron consumption by way of the hydrogenase system in M. marburgensis, which is calculated to be less than 55%.     

Incremental vacuum dehydration-decarbonation experiments on a natural gibbsite (α-Al(OH3)): CO2 abundance and δC values

Incremental vacuum dehydration-decarbonation experiments were performed at 190°C on chemically “cleaned” aliquots of a gibbsite-dominated, Eocene-age bauxite sample with evolution of CO₂ and H₂O. “Plateau” F (CO₂/H₂O ratios) and δ¹³C values of the CO₂ derived from gibbsite were attained over the dehydration interval, X_v(H₂) = 0.16 to 0.67 (i.e., 16 to 67% breakdown of gibbsite). The plateau value of F for gibbsite was 0.0043 ± 0.0003, while the corresponding δ¹³C value of evolved CO₂ was −16.0‰±0.4‰. Additional experiments on chemically cleaned aliquots included (1) treatment with a solution of 0.3M Na-Citrate + 0.1M Na-Dithionite and (2) an exchange experiment with 0.1 bar of ¹³C-depleted CO₂ (−46‰) at 105°C for 64.5 h. Neither of these additional treatments resulted in a measurable perturbation of plateau values of F or δ¹³C for CO₂ evolved from gibbsite during dehydroxylation. These results support published work on Holocene samples which suggested that CO₂ occluded in gibbsite may preserve information on δ¹³C values of CO₂ in ancient terrestrial systems. The plateau values of F observed in the Eocene gibbsite indicate that it may be possible to experimentally calibrate a relationship between the concentration of CO₂ occluded in gibbsite and CO₂ in the environment at the time of crystallization. Such a calibration would significantly enhance the value of gibbsite as a source of information on ancient oxidized carbon systems.     

Surface heat flow and CO2 emissions within the Ohaaki hydrothermal field,Taupo Volcanic Zone,New Zealand

Carbon dioxide emissions and heat flow have been determined from the Ohaaki hydrothermal field, Taupo Volcanic Zone (TVZ), New Zealand following 20 a of production (116 MW_e). Soil CO₂ degassing was quantified with 2663 CO₂ flux measurements using the accumulation chamber method, and 2563 soil temperatures were measured and converted to equivalent heat flow (W m⁻²) using published soil temperature heat flow functions. Both CO₂ flux and heat flow were analysed statistically and then modelled using 500 sequential Gaussian simulations. Forty subsoil CO₂ gas samples were also analysed for stable C isotopes. Following 20 a of production, current CO₂ emissions equated to 111 ± 6.7 T/d. Observed heat flow was 70 ± 6.4 MW, compared with a pre-production value of 122 MW. This 52 MW reduction in surface heat flow is due to production-induced drying up of all alkali–Cl outflows (61.5 MW) and steam-heated pools (8.6 MW) within the Ohaaki West thermal area (OHW). The drying up of all alkali–Cl outflows at Ohaaki means that the soil zone is now the major natural pathway of heat release from the high-temperature reservoir. On the other hand, a net gain in thermal ground heat flow of 18 MW (from 25 MW to 43.3 ± 5 MW) at OHW is associated with permeability increases resulting from surface unit fracturing by production-induced ground subsidence. The Ohaaki East (OHE) thermal area showed no change in distribution of shallow and deep soil temperature contours despite 20 a of production, with an observed heat flow of 26.7 ± 3 MW and a CO₂ emission rate of 39 ± 3 T/d. The negligible change in the thermal status of the OHE thermal area is attributed to the low permeability of the reservoir beneath this area, which has limited production (mass extraction) and sheltered the area from the pressure decline within the main reservoir. Chemistry suggests that although alkali–Cl outflows once contributed significantly to the natural surface heat flow (∼50%) they contributed little (<1%) to pre-production CO₂ emissions due to the loss of >99% of the original CO₂ content due to depressurisation and boiling as the fluids ascended to the surface. Consequently, the soil has persisted as the major (99%) pathway of CO₂ release to the atmosphere from the high temperature reservoir at Ohaaki. The CO₂ flux and heat flow surveys indicate that despite 20 a of production the variability in location, spatial extent and magnitude of CO₂ flux remains consistent with established geochemical and geophysical models of the Ohaaki Field. At both OHW and OHE carbon isotopic analyses of soil gas indicate a two-stage fractionation process for moderate-flux (>60 g m⁻² d⁻¹) sites; boiling during fluid ascent within the underlying reservoir and isotopic enrichment as CO₂ diffuses through porous media of the soil zone. For high-flux sites (>300 g m⁻² d⁻¹), the δ¹³CO₂ signature (−7.4 ± 0.3‰ OHW and −6.5 ± 0.6‰ OHE) is unaffected by near-surface (soil zone) fractionation processes and reflects the composition of the boiled magmatic CO₂ source for each respective upflow. Flux thresholds of <30 g m⁻² d⁻¹ for purely diffusive gas transport, between 30 and 300 g m⁻² d⁻¹ for combined diffusive–advective transport, and ?300 g m⁻² d⁻¹ for purely advective gas transport at Ohaaki were assigned. δ¹³CO₂ values and cumulative probability plots of CO₂ flux data both identified a threshold of ∼15 g m⁻² d⁻¹ by which background (atmospheric and soil respired) CO₂ may be differentiated from hydrothermal CO₂.     

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.     

Graphite-bearing CO2-fluid inclusions in granulites: Insights on graphite precipitation and carbon isotope evolution

Graphite in deep crustal enderbitic (orthopyroxene + garnet + plagioclase + quartz) granulites (740°C, 8.9 kb) of Nilgiri hills, southern India were investigated for their spectroscopic and isotopic characteristics. Four types of graphite crystals were identified. The first type (Gr_I), which is interstitial to other mineral grains, can be grouped into two subtypes, Gr_IA and Gr_IB. Gr_IA is either irregular in shape or deformed, and rough textured with average δ¹³C values of −12.7 ± 0.4‰ (n = 3). A later generation of interstitial graphite (Gr_IB) shows polygonal crystal shapes and highly reflecting smooth surface features. These graphite grains are more common and have δ¹³C values of −11.9 ± 0.3‰ (n = 14). Both subtypes show well-defined Raman shifts suggesting a highly crystalline nature. Cores of interstitial graphite grains have, on average, lower δ¹³C values by ∼0.5‰ compared to that of the rim. The second type of graphite (Gr_II) occurs as solid inclusions in silicate minerals, commonly forming regular hexagonal crystals with a slightly disordered structure. The third type of graphite (Gr_III) is associated with solid inclusions (up to 100 μm) that have decrepitation halos of numerous small (<15 μm) satellite fluid inclusions of pure CO₂ with varying density (1.105 to 0.75 g/cm³). The fourth type of graphite (Gr_IV) is found as daughter crystals within primary type CO₂-fluid inclusions in garnet and quartz. These fluid inclusions have a range of densities (1.05 to 0.90 g/cm³), but in general are significantly less dense than graphite-free primary, pure CO₂ fluid inclusions (1.12 g/cm³). Raman spectral characteristics of graphite inside fluid inclusions suggest graphite crystallization at low temperature (∼ 500°C). The precipitation of graphite probably occurred during the isobaric cooling of CO₂-rich peak metamorphic fluid as a result of oxyexsolution of oxide phases. The oxyexsolution process is evidenced by the magnetite-ilmenite granular exsolution textures and the systematic presence of numerous micron-sized rutile and other oxide inclusions in association with fluid inclusions within garnet, plagioclase, and quartz.The carbon isotope compositions of coexisting CO₂ (in fluid inclusions) and graphite show a fractionation (α₂CO−gr) of ∼6‰ in garnet, consistent with the existing theoretical estimates of α₂CO−gr at 800°C. A subsequent generation of CO₂ inclusions trapped in matrix quartz and quartz segregation have higher δ¹³C values, −4‰ and −2.9‰ respectively. Graphite in quartz segregations also has higher δ¹³C values (−9.8‰) than those in enderbite (−12.7‰). Micro-graphite crystals included in garnet, quartz (enderbite), and quartz (segregation) have average δ¹³C values of −11.1, −10.4, and −8.7‰ respectively, indicating progressive enrichment in ¹³C with a decrease in temperature of recrystallization of respective minerals. This progressive enrichment is also observed in carbon isotope compositions of fluid inclusion CO₂, suggesting isotopic equilibrium during graphite precipitation from CO₂ fluids. Thus, the carbon isotope record preserved in these rocks by the interstitial graphite, CO₂ fluid in enderbite, graphite microcrystals, graphite in quartz segregation, and CO₂ fluid in quartz segregation, suggests a temperature-controlled isotopic evolution. This evolution is in accordance with a closed system Rayleigh-type graphite precipitation process which progressively enriched residual CO₂ in ¹³C.     

A model for O/O variations in CO2 evolved from goethite during the solid-state α-FeOOH to α-Fe2O3 phase transition

“Plateau” δ¹⁸O values of CO₂ that evolved from the Fe(CO₃)OH component during isothermal vacuum dehydrations (200-230 °C) of 18 natural goethites range from 8.2 to 28.1‰. In contrast, the measured δ¹⁸O values of the goethite structural oxygen range from −11.3 to 1.7‰. The results of this study indicate that the apparent oxygen isotope fractionation factor (¹⁸α_app) between plateau CO₂ and initial goethite is systematically related to the rate of isothermal vacuum dehydration. The nonlinear correlation and the magnitudes of the ¹⁸α_app values are predicted by a relatively simple mass balance model with the following assumptions: (1) the rate of isothermal vacuum dehydration of goethite (for the interval from 0 to ∼60 to 80% loss of structural hydroxyl hydrogen) can be reasonably well represented by first-order kinetics and (2) isotopic exchange between evolving H₂O vapor and solid occurs only in successive, local transition states. The generally good correspondence between the model predictions and the experimental data seems to validate these assumptions. Thus, the ¹⁸O/¹⁶O ratios of the evolved CO₂ can act as probes into the transient processes operating at the molecular level during the solid-state goethite-to-hematite phase transition. For example, the activation energy for the rate constant associated with the transition state, oxygen isotopic exchange between solid and H₂O vapor, is tentatively estimated as 28 ± 11 KJ/mol. Such knowledge may be of consequence in understanding the significance of ¹⁸O/¹⁶O ratios in hematites from some natural environments (e.g., Mars?).Kinetic data and δ¹⁸O values of CO₂ are routinely obtained in the course of measurements of the abundance and δ¹³C values of the Fe(CO₃)OH in goethite. The observed correlation between ¹⁸α_app and dehydration rates suggests that plateau δ¹⁸O values of evolved CO₂ may provide complementary estimates of the δ¹⁸O values of total goethite structural oxygen (O, OH, CO₂) with an overall precision of about ±1‰. However, because of isotopic exchange during the dehydration process, δ¹⁸O values of the evolved CO₂ do not reflect the original δ¹⁸O values of the CO₂ that was occluded as Fe(CO₃)OH in goethite.     

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.