Correction of ground-water chemistry and carbon isotopic composition for effects of CO2 outgassing

Direct measurements on water samples from several CO₂-charged warm springs are significantly higher than values calculated from field pH and alkalinity (and other constituents). In addition, calcite saturation indices calculated from field pH and solution composition indicated supersaturation in samples which, on the basis of hydrogeologic concepts, should be near saturation or undersaturated. We attribute these discrepancies to uncertainties in field pH, resulting from CO₂ outgassing during pH measurement. Because samples for direct measurement can be taken with minimal disturbance to the water chemistry, we have used the measured to back calculate an estimate of the field pH and the carbon isotopic composition of the water before outgassing. By reconstructing water chemistry in this way, we find generally consistent grouping of $\text{δ}{}^{13}\text{C}$, pH, and degree of calcite saturation in samples taken from the same source at different times, an observation which we expect based on our understanding of the hydrogeology and geochemistry of the ground-water systems. This suggests that for very careful geochemical work, particularly on ground-waters much above ambient temperature, measurements may provide more information on the system and a better estimate of its state of saturation with respect to carbonate minerals than can field measurements of pH.     

A comparison of isotope fractionation of carbon and hydrogen from paddy field rice roots and soil bacterial enrichments during CO2/H2 methanogenesis

To better understand the isotope biogeochemistry of paddy field CH₄, we investigated carbon and hydrogen isotope fractionation during CO₂ reduction by a methanogenic community enriched from California paddy field soil and rice plants. Results from analyses of terminal restriction fragment length polymorphism (T-RFLP) and sequences of the archaeal small-subunit (SSU) rRNA-encoding genes (rDNA) showed a difference in methanogenic community structure between the soil (dominated by Methanobacteriaceae) and roots (dominated by Methanospirillaceae) which was essentially the same for sampling dates 15 and 99 days after flooding (DAF). CO₂/H₂ methanogenesis by these microbial communities produced CH₄ with different isotope ratios and fractionation factors (α factors). The carbon isotope α factors in an open system with a continuous supply of 0.5% H₂ were 1.050 ± 0.002 and 1.057 ± 0.001 for soil and root enrichment cultures at 15 DAF, and 1.052 ± 0.0.002 and 1.059 ± 0.002 for soil and root enrichment cultures at 99 DAF, respectively. These α factors are similar to, but distinct from values previously obtained from cultures of mesophilic methanogens and are larger than calculated values (1.045) for paddy soil. Fractionation of hydrogen isotopes was also studied in a closed system under 80% H₂. The difference in α factors between soil and root enrichment cultures remained clear. The hydrogen isotope fractionations between culture water and the product CH₄ were −327 ± 14‰ and −319 ± 18‰ for soil enrichments, and −389 ± 17‰ and −382 ± 21‰ for root enrichments at 15 DAF and 99 DAF, respectively.     

Stable carbon isotope composition and concentrations of CO2 and CH4 in the deep catotelm of a peat bog

Vertical profiles of concentration and C-isotopic composition of dissolved methane and carbon dioxide were observed over 26 months in the catotelm of a deep (6.5 m) peat bog in Switzerland. The dissolved concentrations of these gases increase with depth while CO₂ predominates over CH₄ (CO₂ ca. 5 times CH₄). This pattern can be reproduced by a reaction-advection-ebullition model, where CO₂ and CH₄ are formed in a ratio of 1:1. The less soluble methane is preferentially lost via outgassing (bubbles). The isotopic fractionation between CO₂ and CH₄ also increases with depth, with α_C values ranging from 1.045 to 1.075. The isotopic composition of the gases traces the passage of respiration-derived CO₂ (from the near surface) through a shallow zone with methanogenesis of low isotopic fractionation (splitting of fermentation-derived acetate). This solution then moves through the catotelm, where methanogenesis occurs by CO₂ reduction (large isotopic fractionation). In the upper part of the catotelm the C-13-depleted respiration-derived CO₂ pool buffers the isotopic composition of CO₂; the δ¹³C of CO₂ increases only slowly. At the same time strongly depleted CH₄ is formed as CO₂ reduction consumes the depleted CO₂. In the lower part of the catotelm, the respiration-derived CO₂ and shallow CH₄ become less important and CO₂ reduction is the dominant source of CO₂ and CH₄. Now, the δ¹³C values of both gases increase until equilibrium is reached with respect to the isotopic composition of the substrate. Thus, the δ¹³C values of methane reach a minimum at intermediate depth, and the deep methane has δ¹³C values comparable to shallow methane. A simple mixing model for the isotopic evolution is suggested. Only minor changes of the observed patterns of methanogenesis (in terms of concentration and isotopic composition) occur over the seasons. The most pronounced of these is a slightly higher rate of acetate splitting in spring.     

BRYOCARB: A process-based model of thallose liverwort carbon isotope fractionation in response to CO2, O2, light and temperature

Evidence from laboratory experiments indicates that fractionation against the heavy stable isotope of carbon (Δ¹³C) by bryophytes (liverworts and mosses) is strongly dependent on atmospheric CO₂. This physiological response may therefore provide the basis for developing a new terrestrial CO₂ proxy [Fletcher, B.J., Beerling, D.J., Brentnall, S.J., Royer, D.L., 2005. Fossil bryophytes as recorders of ancient CO₂ levels: experimental evidence and a Cretaceous case study. Global Biogeochem. Cycles19, GB3012]. Here, we establish a theoretical basis for the proxy by developing an extended model of bryophyte carbon isotope fractionation (BRYOCARB) that integrates the biochemical theory of photosynthetic CO₂ assimilation with controls on CO₂ supply by diffusion from the atmosphere. The BRYOCARB model is evaluated against measurements of the response of liverwort photosynthesis and Δ¹³C to variations in atmospheric O₂, temperature and irradiance at different CO₂ concentrations. We show that the bryophyte proxy is at least as sensitive to variations in atmosphere CO₂ as the two other leading carbon isotope-based approaches to estimating palaeo-CO₂ levels (δ¹³C of phytoplankton and of paleosols). Mathematical inversion of BRYOCARB provides a mechanistic means of estimating atmospheric CO₂ levels from fossil bryophyte carbon that can explicitly account for the effects of past differences in O₂ and climate.     

Temperature dependence of oxygen isotope fractionation of CO2 from magnesite-phosphoric acid reaction

New experimental results are reported on oxygen isotope fractionation factors, α_T, between the δ¹⁸O compositions of carbon dioxide liberated by phosphoric acid in the temperature interval of 323 to 373K and that of total oxygen from a natural magnesite (MgCO₃). These results are distinctly different from some previously published mutually inconsistent data, and can be expressed as a linear relationship: 10³ lnα_T = [{(6.845 ± 0.475)∗10⁵}/T²] + (4.22 ± 0.08), where 10³ lnα_T refers to fractionation at different temperatures T in Kelvin.Fractionation factors have also been determined at 323 and 368K on a natural calcite. The results on calcite are in excellent agreement with previously published data and can be written as:10³ lnα_T = [{(5.608 ± 0.151)∗10⁵}/T²] + (3.89 ± 0.08).The combined results on magnesite and calcite yield a computed value of α = 1.01117 for dolomite at 298K, assuming equal proportion of 0.5 mole of magnesium and calcium in dolomite, the previously reported experimental values being 1.01109 and 1.01110.     

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.     

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.     

The yield and isotopic composition of radiolytic H2, a potential energy source for the deep subsurface biosphere

The production rate and isotopic composition of H₂ derived from radiolytic reactions in H₂O were measured to assess the importance of radiolytic H₂ in subsurface environments and to determine whether its isotopic signature can be used as a diagnostic tool. Saline and pure, aerobic and anaerobic water samples with pH values of 4, 7, and 10 were irradiated in sealed vials at room temperature with an artificial γ source, and the H₂ abundance in the headspace and its isotopic composition were measured. The H₂ concentrations were observed to increase linearly with dosage at a rate of 0.40 ± 0.04 molecules (100 eV)⁻¹ within the dosage range of 900 to 3500 Gray (Gy; Gy = 1 J Kg⁻¹) with no indication of a maximum limit on H₂ concentration. At ∼2000 Gy, the H₂ concentration varied only by 16% across the experimental range of pH, salinity, and O₂. Based upon this measured yield and H₂ yields for α and β particles, a radiolytic H₂ production rate of 10⁻⁹ to 10⁻⁴ nM s⁻¹ was estimated for the range of radioactive element concentrations and porosities typical of crustal rocks. The δD of H₂ was independent of the dosage, pH (except for pH 4), salinity, and O₂ and yielded an αD_H2O-H2 of 2.05 ± 0.07 (αD_H2O-H2 = (D/H)_H2O to (D/H)_H2), slightly less than predicted radiolytic models. Although this radiolytic fractionation value is significantly heavier than that of equilibrium isotopic exchange between H₂ and H₂O, the isotopic exchange rate between H₂ and H₂O will erase the heavy δD of radiolytic H₂ if the age of the groundwater is greater than ∼10³ to 10⁴ yr. The millimolar concentrations of H₂ observed in the groundwater of several Precambrian Shields are consistent with radiolysis of water that has resided in the subsurface for a few million years. These concentrations are well above those required to support H₂-utilizing microorganisms and to inhibit H₂-producing, fermentative microorganisms.     

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.     

Impact of climate and CO2 on a millennium-long tree-ring carbon isotope record

We present one millennium-long (1171-year), and three 100 year long annually resolved δ¹³C tree-ring chronologies from ecologically varying Juniperus stands in the Karakorum Mountains (northern Pakistan), and evaluate their response to climatic and atmospheric CO₂ changes. All δ¹³C records show a gradual decrease since the beginning of the 19th century, which is commonly associated with a depletion of atmospheric δ¹³C due to fossil fuel burning. Climate calibration of high-frequency δ¹³C variations indicates a pronounced summer temperature signal for all sites. The low-frequency component of the same records, however, deviates from long-term temperature trends, even after correction for changes in anthropogenic CO₂. We hypothesize that these high-elevation trees show a response to both climate and elevated atmospheric CO₂ concentration and the latter might explain the offset with target temperature data. We applied several corrections to tree-ring δ¹³C records, considering a range of potential CO₂ discrimination changes over the past 150 years and calculated the goodness of fit with the target via calibration/verification tests (R², residual trend, and Durbin-Watson statistics). These tests revealed that at our sites, carbon isotope fixation on longer timescales is affected by increasing atmospheric CO₂ concentrations at a discrimination rate of about 0.012‰/ppmv. Although this statistically derived value may be site related, our findings have implications for the interpretation of any long-term trends in climate reconstructions using tree-ring δ¹³C, as we demonstrate with our millennium-long δ¹³C Karakorum record. While we find indications for warmth during the Medieval Warm Period (higher than today’s mean summer temperature), we also show that the low-frequency temperature pattern critically depends on the correction applied. Patterns of long-term climate variation, including the Medieval Warm Period, the Little Ice Age, and 20th century warmth are most similar to existing evidence when a strong influence of increased atmospheric CO₂ on plant physiology is assumed.     

Experimental determination of carbon isotope fractionation between CaCO3 and graphite

Carbon isotopic exchange between graphite and three polymorphs of CaCO₃ was investigated at temperatures of 600-1400 °C and at pressures from 1.4 to 2.3 GPa. Fractionation factors at all temperatures were determined by the partial exchange treatment of Northrop and Clayton (1966).Graphite starting material for the majority of the experiments was milled in water for 20-25 h, producing aggregates of nanosheets. The sheets range in width from 50 to 1000 nm and in thickness from 20 to 30 nm, and they retain hexagonal symmetry.Isotopic exchange appears to be the sum of surface exchange and interior exchange. At 1100-1400 °C, interior exchange exceeded surface exchange, probably by a combination of grain growth, as determined by increase in crystallite size, recrystallization, as observed in FESEM images, and diffusion. In some runs at 1200 and 1400 °C with an isotopic contrast between the initial graphite and calcite of close to 50‰, equilibrium fractionation was actually overstepped due to a kinetic effect. A weighted regression of fractionation factors from the high-temperature runs yields the line of equilibrium interior exchange:

     

Kinetic isotopic fractionation during carbonate dissolution in laboratory experiments: Implications for detection of microbial CO2 signatures using δC-DIC

Laboratory experiments on reagent-grade calcium carbonate and carbonate rich glacial sediments demonstrate previously unreported kinetic fractionation of carbon isotopes during the initial hydrolysis and early stages of carbonate dissolution driven by atmospheric CO₂. There is preferential dissolution of Ca¹²CO₃ during hydrolysis, resulting in δ¹³C-DIC values that are significantly lighter isotopically than the bulk carbonate. The fractionation factor for this kinetic isotopic effect is defined as ε_carb. ε_carb is greater on average for glacial sediments (−17.4‰) than for calcium carbonate (−7.8‰) for the < 63 μm size fraction, a sediment concentration of 5 g L⁻¹ and closed system conditions at 5°C. This difference is most likely due to the preferential dissolution of highly reactive ultra-fine particles with damaged surfaces that are common in subglacial sediments. The kinetic isotopic fractionation has a greater impact on δ¹³C-DIC at higher CaCO₃:water ratios and is significant during at least the first 6 h of carbonate dissolution driven by atmospheric CO₂ at sediment concentrations of 5 g L⁻¹. Atmospheric CO₂ dissolving into solution following carbonate hydrolysis does not exhibit any significant equilibrium isotopic fractionation for at least ∼ 6 h after the start of the experiment at 5°C. This is considerably longer than previously reported in the literature. Thus, kinetic fractionation processes will likely dominate the δ¹³C-DIC signal in natural environments where rock:water contact times are short <6-24 h (e.g., glacial systems, headwaters in fluvial catchments) and there is an excess of carbonate in the sediments. It will be difficult apply conventional isotope mass balance techniques in these types of environment to identify microbial CO₂ signatures in DIC from δ¹³C-DIC data.     

Argon and CO2 on the race track in silicate melts: A tool for the development of a CO2 speciation and diffusion model

We have analysed the kinetics of Argon and CO₂ diffusion in simplified iron free rhyolitic to hawaiitic melts using the diffusion couple technique. The concentration distance profiles of Ar and CO₂ were measured with electron microprobe analysis and Fourier Transform Infrared Spectroscopy, respectively. Error functions were fitted to the symmetrical concentration distance profiles to extract the diffusion coefficients.In the temperature range 1373 to 1773 K the activation energies for Ar diffusion range from 169 ± 20 to 257 ± 62 kJ mol⁻¹. Ar diffusivity increases exponentially with the degree of depolymerisation. In contrast, the mobility of total CO₂, that is identical to Ar mobility in rhyolitic melt, keeps constant with changing bulk composition from rhyolite to hawaiite. CO₂ speciation at 1623 K and 500 MPa was modeled for the range of compositions studied using the diffusion data of Ar and total CO₂ in combination with network former diffusion calculated from viscosity data. Within error this model is in excellent agreement with CO₂ speciation data extrapolated from temperatures near the glass transition temperature for dacitic melt composition. This model shows that even in highly depolymerised hawaiitic and tholeiitic melts molecular CO₂ is a stable species and contributes 70 to 80% to the total CO₂ diffusion, respectively.     

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.     

咸水层CO2地质封存地下利用空间评估方法研究

咸水层CO2地质封存技术是我国实现碳中和目标的重要支撑技术,也是一项深部地下空间开发利用技术。咸水层CO2地质封存工程利用的深部地下空间,需要在确定CO2羽流、扰动边界和经济因素“三级边界”的基础上进行综合评估。以我国唯一的深部咸水层CO2地质封存项目——国家能源集团鄂尔多斯碳捕集与封存(Carbon Capture and Storage, CCS)示范工程为实例,基于封存场地储层CO2羽流监测以及扰动边界的推断预测结果综合评估,认为示范工程平面上4个1'×1'经纬度范围可作为地下利用空间平面边界,垂向上以纸坊组顶界(深度约958 m)为地下封存体顶部边界,以深度2 800 m为底板封隔层底界。提出的咸水层CO2地质封存地下利用空间评估方法,能够为未来封存工程地下利用空间审批与监管提供一定参考,但也需要进一步结合已有法律法规及规模化封存工程实践完善提升。     

GEOCARBSULF: A combined model for Phanerozoic atmospheric O2 and CO2

A model for the combined long-term cycles of carbon and sulfur has been constructed which combines all the factors modifying weathering and degassing of the GEOCARB III model [Berner R.A., Kothavala Z., 2001. GEOCARB III: a revised model of atmospheric CO₂ over Phanerozoic time. Am. J. Sci. 301, 182-204] for CO₂ with rapid recycling and oxygen dependent carbon and sulfur isotope fractionation of an isotope mass balance model for O₂ [Berner R.A., 2001. Modeling atmospheric O₂ over Phanerozoic time. Geochim. Cosmochim. Acta65, 685-694]. New isotopic data for both carbon and sulfur are used and new feedbacks are created by combining the models. Sensitivity analysis is done by determining (1) the effect on weathering rates of using rapid recycling (rapid recycling treats carbon and sulfur weathering in terms of young rapidly weathering rocks and older more slowly weathering rocks); (2) the effect on O₂ of using different initial starting conditions; (3) the effect on O₂ of using different data for carbon isotope fractionation during photosynthesis and alternative values of oceanic δ¹³C for the past 200 million years; (4) the effect on sulfur isotope fractionation and on O₂ of varying the size of O₂ feedback during sedimentary pyrite formation; (5) the effect on O₂ of varying the dependence of organic matter and pyrite weathering on tectonic uplift plus erosion, and the degree of exposure of coastal lands by sea level change; (6) the effect on CO₂ of adding the variability of volcanic rock weathering over time [Berner, R.A., 2006. Inclusion of the weathering of volcanic rocks in the GEOCARBSULF model. Am. J. Sci.306 (in press)]. Results show a similar trend of atmospheric CO₂ over the Phanerozoic to the results of GEOCARB III, but with some differences during the early Paleozoic and, for variable volcanic rock weathering, lower CO₂ values during the Mesozoic. Atmospheric oxygen shows a major broad late Paleozoic peak with a maximum value of about 30% O₂ in the Permian, a secondary less-broad peak centered near the Silurian/Devonian boundary, variation between 15% and 20% O₂ during the Cambrian and Ordovician, a very sharp drop from 30% to 15% O₂ at the Permo-Triassic boundary, and a more-or less continuous rise in O₂ from the late Triassic to the present.     

Indoor and outdoor urban atmospheric CO2: Stable carbon isotope constraints on mixing and mass balance

From July to November 2009, concentrations of CO₂ in 78 samples of ambient air collected in 18 different interior spaces on a university campus in Dallas, Texas (USA) ranged from 386 to 1980 ppm. Corresponding δ¹³C values varied from −8.9‰ to −19.4‰. The CO₂ from 22 samples of outdoor air (also collected on campus) had a more limited range of concentrations from 385 to 447 ppm (avg. = 408 ppm), while δ¹³C values varied from −10.1‰ to −8.4‰ (avg.=-9.0‰). In contrast to ambient indoor and outdoor air, the concentrations of CO₂ exhaled by 38 different individuals ranged from 38,300 to 76,200 ppm (avg. = 55,100 ppm), while δ¹³C values ranged from −24.8‰ to −17.7‰ (avg. = −21.8‰). The residence times of the total air in the interior spaces of this study appear to have been on the order of 10 min with relatively rapid approaches (∼30 min) to steady-state concentrations of ambient CO₂ gas. Collectively, the δ¹³C values of the indoor CO₂ samples were linearly correlated with the reciprocal of CO₂ concentration, exhibiting an intercept of −21.8‰, with r² = 0.99 and p < 0.001 (n = 78). This high degree of linearity for CO₂ data representing 18 interior spaces (with varying numbers of occupants), and the coincidence of the intercept (−21.8‰) with the average δ¹³C value for human-exhaled CO₂ demonstrates simple mixing between two inputs: (1) outdoor CO₂ introduced to the interior spaces by ventilation systems, and (2) CO₂ exhaled by human occupants of those spaces. If such simple binary mixing is a common feature of interior spaces, it suggests that the intercept of a mixing line defined by two data points (CO₂ input from the local ventilation system and CO₂ in the ambient air of the room) could be a reasonable estimate of the average δ¹³C value of the CO₂ exhaled by the human occupants. Thus, such indoor spaces appear to constitute effective “sample vessels” for collection of CO₂ that can be used to determine the average proportions of C₃ and C₄-derived C in the diets of the occupants. For the various groups occupying the rooms sampled in this study, C₄-derived C appears to have constituted ∼40% of the average diet.     

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

Evolution of the isotopic composition of carbon and oxygen in a calcite precipitating H2O-CO2-CaCO3 solution and the related isotopic composition of calcite in stalagmites

The isotopic composition of carbon and oxygen in a calcite precipitating CO₂-H₂O-CaCO₃ solution is preserved in the calcite precipitated. For the interpretation of isotopic proxies from stalagmites knowledge of the evolution of δ¹³C and δ¹⁸O in the solution during precipitation is required. A system of differential equations is presented from which this evolution can be derived. Both, irreversible loss of carbon and oxygen from the solution with precipitation time τ and exchange of oxygen in the carbonates with the oxygen in the water with exchange time T are considered. For carbon, where no exchange is active, a modified equation of Rayleigh-distillation is found, which takes into account that precipitation stops at c_eq, the saturation concentration of DIC with respect to calcite, and that c_eq as well as the precipitation time τ is slightly different for the heavy and the light isotope. This, however, requires introducing a new parameter γ = (A_eq/B_eq)/(A₀/B₀), which has to be determined experimentally. (A_eq/B_eq) is the isotopic ratio for the heavy (A) and the light isotope (B) at both chemical and isotopic equilibrium and (A₀/B₀) is the initial isotopic ratio of the solution. In the case of oxygen, where exchange is present, the isotopic shifts are reduced with increasing values of the precipitation time τ. For τ ? T the solution stays in isotopic equilibrium with the oxygen in the water during the entire time in which precipitation is active. The isotopic ratios in a calcite precipitating solution R(t)/R₀ = (1 + δ(t)/1000) for carbon are plotted versus those of oxygen. R₀ is the isotopic ratio at time t = 0, when precipitation starts and δ(t) the isotopic shift in the solution after time t. These show positive correlations for the first 50% of calcite, which can precipitate. Their slopes increase with increasing values of τ and they closely resemble Hendy-tests performed along growth layers of stalagmites. Our results show that stalagmites, which grow by high supply of water with drip times less than 50 s, exhibit positive correlations between δ¹³C and δ¹⁸O along a growth layer. But in spite of this the isotopic composition of oxygen in the solution at the apex is in isotopic equilibrium with the oxygen in the water, and therefore also that of calcite deposited at the apex.