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.     

A high temperature equation of state for the H2O-CaCl2 and H2O-MgCl2 systems

An equation of state (EOS) is developed for salt-water systems in the high temperature range. As an example of the applications, this EOS is parameterized for the calculation of density, immiscibility, and the compositions of coexisting phases in the CaCl₂-H₂O and MgCl₂-H₂O systems from 523 to 973 K and from saturation pressure to 1500 bar. All available volumetric and phase equilibrium measurements of these binaries are well represented by this equation. This EOS is based on a Helmholtz free energy representation constructed from a reference system containing hard-sphere and polar contributions plus an empirical correction. For the temperature and pressure range in this study, the electrolyte solutes are assumed to be associated. The water molecules are modeled as hard spheres with point dipoles and the solute molecules, MgCl₂ and CaCl₂, as hard spheres with point quadrupoles. The free energy of the reference system is calculated from an analytical representation of the Helmholtz free energy of the hard-sphere contributions and perturbative estimates of the electrostatic contributions. The empirical correction used to account for deviations of the reference system predictions from measured data is based on a virial expansion. The formalism allows generalization to aqueous systems containing insoluble gases (CO₂, CH₄), alkali chlorides (NaCl, KCl), and alkaline earth chlorides (CaCl₂, MgCl₂). The program of this model is available as an electronic annex (see EA1 and EA2) and can also be downloaded at: http://www.geochem-model.org/programs.htm.     

Evaluation of carbon isotope fractionation of soil CO2 under an advective-diffusive regimen: A tool for computing the isotopic composition of unfractionated deep source

A physical model based on the advective-diffusion theory was developed in order to describe the mixing between a deep gas source and the atmosphere. The model was used to predict the isotopic fractionation of carbon in soil CO₂. Gas samples were collected at different depths in areas characterized by different geological settings and CO₂ fluxes. The relative theoretical and experimental isotopic profiles were compared and a good agreement was found. These profiles show how the isotopic composition of CO₂ changes through the upper few decimeters of soil and how the amount of the isotopic fractionation is strongly influenced by soil CO₂ flux. Finally, the model was used to derive the carbon isotopic composition of unfractioned deep CO₂ source for all the investigated sites.     

Equation of state of the H2O, CO2, and H2O-CO2 systems up to 10 GPa and 2573.15 K: Molecular dynamics simulations with ab initio potential surface

Based on our previous development of the molecular interaction potential for pure H₂O and CO₂ [Zhang, Z.G., Duan, Z.H. 2005a. Isothermal-isobaric molecular dynamics simulations of the PVT properties of water over wide range of temperatures and pressures. Phys. Earth Planet Interiors149, 335-354; Zhang, Z.G., Duan, Z.H. 2005b. An optimized molecular potential for carbon dioxide. J. Chem. Phys.122, 214507] and the ab initio potential surface across CO₂-H₂O molecules constructed in this study, we carried out more than one thousand molecular dynamics simulations of the PVTx properties of the CO₂-H₂O mixtures in the temperature-pressure range from 673.15 to 2573.15 K up to 10.0 GPa. Comparison with extensive experimental PVTx data indicates that the simulated results generally agree with experimental data within 2% in density, equivalent to experimental uncertainty. Even the data under the highest experimental temperature-pressure conditions (up to 1673 K and 1.94 GPa) are well predicted with the agreement within 1.0% in density, indicating that the high accuracy of the simulation is well retained as the temperature and pressure increase. The consistent and stable predictability of the simulation from low to high temperature-pressure and the fact that the molecular dynamics simulation resort to no experimental data but to ab initio molecular potential makes us convinced that the simulation results should be reliable up to at least 2573 K and 10 GPa with errors less than 2% in density. In order to integrate all the simulation results of this study and previous studies [Zhang and Duan, 2005a, 2005b] and the experimental data for the calculation of volumetric properties (volume, density, and excess volume), heat properties, and chemical properties (fugacity, activity, and possibly supercritical phase separation), an equation of state (EOS) is laboriously developed for the CO₂, H₂O, and CO₂-H₂O systems. This EOS reproduces all the experimental and simulated data covering a wide temperature and pressure range from 673.15 to 2573.15 K and from 0 to 10.0 GPa within experimental or simulation uncertainty.     

Oxidative dissolution potential of biogenic and abiogenic TcO2 in subsurface sediments

Technetium-99 (Tc) is an important fission product contaminant associated with sites of nuclear fuels reprocessing and geologic nuclear waste disposal. Tc is highly mobile in its most oxidized state and less mobile in the reduced form [Tc(IV)O₂·nH₂O]. Here we investigate the potential for oxidation of Tc(IV) that was heterogeneously reduced by reaction with biogenic Fe(II) in two sediments differing in mineralogy and aggregation state; unconsolidated Pliocene-age fluvial sediment from the upper Ringold (RG) Formation at the Hanford Site and a clay-rich saprolite from the Field Research Center (FRC) background site on the Oak Ridge Site. Both sediments contained Fe(III) and Mn(III/IV) as redox active phases, but FRC also contained mass-dominant Fe-phyllosilicates of different types. Shewanella putrefaciens CN32 reduced Mn(III/IV) oxides and generated Fe(II) that was reactive with Tc(VII) in heat-killed, bioreduced sediment. After bioreduction and heat-killing, biogenic Fe(II) in the FRC exceeded that in RG by a factor of two. More rapid reduction rates were observed in the RG that had lower biogenic Fe(II), and less particle aggregation. EXAFS measurements indicated that the primary reduction product was a TcO₂-like phase in both sediments. The biogenic redox product Tc(IV) oxidized rapidly and completely in RG when contacted with air. Oxidation, in contrast, was slow and incomplete in the FRC, in spite of similar molecular scale speciation of Tc compared to RG. X-ray microprobe, electron microprobe, X-ray absorption spectroscopy, and micro X-ray diffraction were applied to the whole sediment and isolated Tc-containing particles. These analyses revealed that non-oxidizable Tc(IV) in the FRC existed as complexes with octahedral Fe(III) within intra-grain domains of 50-100 μm-sized, Fe-containing micas presumptively identified as celadonite. The markedly slower oxidation rates in FRC as compared to RG were attributed to mass-transfer-limited migration of O₂ into intra-aggregate and intraparticle domains where Tc(IV) existed; and the formation of unique, oxidation-resistant, intragrain Tc(IV)-Fe(III) molecular species.     

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.     

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.     

The noble gas systematics of late-orogenic H2O-CO2 fluids, Mt Isa, Australia

The noble gases (He, Ne, Ar, Kr and Xe) are powerful geochemical tracers because they have distinctive isotopic compositions in the atmosphere, crust and mantle. This study illustrates how noble gases can be used to trace fluid origins in high-temperature metamorphic and mineralising environments; and at the same time provides new information on the composition of noble gases in deeper parts of the crust than have been sampled previously.We report data for H₂O and CO₂ fluid inclusions trapped at greenschist to amphibolite facies metamorphic conditions associated with three different styles of mineralisation and alteration in the Proterozoic Mt Isa Inlier of Australia. Sulphide fluid inclusions are dominated by crustal ⁴He. However, co-variations in fluid inclusion ²⁰Ne/²²Ne, ²¹Ne/²²Ne, ⁴⁰Ar/³⁶Ar and ¹³⁶Xe/¹³⁰Xe indicate noble gases were derived from three or more reservoirs. In most cases, the fluid inclusions elemental noble gas ratios (e.g. Ne/Xe) are close to the ranges expected in sedimentary and crystalline rocks. However, the elemental ratios have been modified in some of the samples providing evidence for independent pulses of CO_2, and interaction of CO₂ with high-salinity aqueous fluids.Compositional variation is attributed to mixing of: (i) magmatic fluids (or deeply sourced metamorphic fluids) characterised by basement-derived noble gases with ²⁰Ne/²²Ne ∼ 8.4, ²¹Ne/²²Ne ∼ 0.4, ⁴⁰Ar/³⁶Ar ∼ 40,000 and ¹³⁶Xe/¹³⁰Xe ∼ 8; (ii) basinal-metamorphic fluids with a narrow range of compositions including near-atmospheric values and (iii) noble gases derived from the meta-sedimentary host-rocks with ²⁰Ne/²²Ne ∼ 8-9.8, ²¹Ne/²²Ne < 0.1, ⁴⁰Ar/³⁶Ar < 2500 and ¹³⁶Xe/¹³⁰Xe ∼ 2.2.These data provide the strongest geochemical evidence available for the involvement of fluids from two distinct geochemical reservoirs in Mt Isa’s largest ore deposits. In addition the data show how noble gases in fluid inclusions can provide information on fluid origins, the composition of the crust’s major lithologies, fluid-rock interactions and fluid-fluid mixing or immiscibility processes.     

Numerical model to determine the composition of H2O-NaCl-CaCl2 fluid inclusions based on microthermometric and microanalytical data

Natural fluids approximated by the H₂O-NaCl-CaCl₂ system are common in a wide range of geologic environments, including sedimentary basins associated with hydrocarbon occurrences and MVT deposits, submarine hydrothermal systems, and other metamorphic, magmatic and hydrothermal environments. We present a comprehensive numerical model and Microsoft® Excel©-based computer program to determine the compositions of fluid inclusions in the H₂O-NaCl-CaCl₂ system based on microthermometric and microanalytical data. The model consists of six polynomial correlation equations that describe liquid salinity as a function of NaCl/CaCl₂ ratio and melting temperature on each of the ice, hydrohalite, halite, antarcticite, CaCl₂·4H₂O and CaCl₂·2H₂O vapor-saturated liquidus surfaces. The cotectic and peritectic boundaries are determined from the intersections of the liquidus surfaces. The model is implicitly internally consistent and topologically correct.The model expands upon the compositional range of applicability and the data types that can be used for compositional determination. It reproduces experimental data for all compositions that lie within the H₂O-NaCl-CaCl₂·4H₂O compositional triangle in the H₂O-NaCl-CaCl₂ system and yields accurate reproductions of the H₂O-NaCl and H₂O-CaCl₂ binaries. Furthermore, in comparison to previously published models, the one presented here eliminates systematic errors, wavy isotherms and cotectic and peritectic curves with local “bumps.”     

CO2 solubility and speciation in intermediate (andesitic) melts: the role of H2O and composition

We determined total CO₂ solubilities in andesite melts with a range of compositions. Melts were equilibrated with excess C-O(-H) fluid at 1 GPa and 1300°C then quenched to glasses. Samples were analyzed using an electron microprobe for major elements, ion microprobe for C-O-H volatiles, and Fourier transform infrared spectroscopy for molecular H₂O, OH⁻, molecular CO₂, and CO₃²⁻. CO₂ solubility was determined in hydrous andesite glasses and we found that H₂O content has a strong influence on C-O speciation and total CO₂ solubility. In anhydrous andesite melts with ∼60 wt.% SiO₂, total CO₂ solubility is ∼0.3 wt.% at 1300°C and 1 GPa and total CO₂ solubility increases by about 0.06 wt.% per wt.% of total H₂O. As total H₂O increases from ∼0 to ∼3.4 wt.%, molecular CO₂ decreases (from 0.07 ± 0.01 wt.% to ∼0.01 wt.%) and CO₃²⁻ increases (from 0.24 ± 0.04 wt.% to 0.57 ± 0.09 wt.%). Molecular CO₂ increases as the calculated mole fraction of CO₂ in the fluid increases, showing Henrian behavior. In contrast, CO₃²⁻ decreases as the calculated mole fraction of CO₂ in the fluid increases, indicating that CO₃²⁻ solubility is strongly dependent on the availability of reactive oxygens in the melt. These findings have implications for CO₂ degassing. If substantial H₂O is present, total CO₂ solubility is higher and CO₂ will degas at relatively shallow levels compared to a drier melt. Total CO₂ solubility was also examined in andesitic glasses with additional Ca, K, or Mg and low H₂O contents (<1 wt.%). We found that total CO₂ solubility is negatively correlated with (Si + Al) cation mole fraction and positively correlated with cations with large Gibbs free energy of decarbonation or high charge-to-radius ratios (e.g., Ca). Combining our andesite data with data from the literature, we find that molecular CO₂ is more abundant in highly polymerized melts with high ionic porosities (>∼48.3%), and low nonbridging oxygen/tetrahedral oxygen (<∼0.3). Carbonate dominates most silicate melts and is most abundant in depolymerized melts with low ionic porosities, high nonbridging oxygen/tetrahedral oxygen (>∼0.3), and abundant cations with large Gibbs free energy of decarbonation or high charge-to-radius ratio. In natural silicate melt, the oxygens in the carbonate are likely associated with tetrahedral and network-modifying cations (including Ca, H, or H-bonds) or a combinations of those cations.     

TSR成因H2S的硫同位素分馏特征与机制

热化学硫酸盐还原反应(TSR)是深层碳酸盐岩油气藏中硫化氢的主要成因机制,目前已在全球发现了50多个TSR成因的大中型含硫化氢天然气田。通过对中国四川盆地含硫化氢气田硫化物的采集与同位素分析,结合全球含硫化氢天然气田硫同位素分析数据,研究了TSR过程中硫同位素的地球化学行为和分馏特征。研究发现,TSR成因的高含硫化氢天然气中,硫化氢与硫酸盐的硫同位素分馏值小于15‰,主要分布范围为2.5‰~13.82‰,平均在10‰。四川盆地海相层系膏岩的硫同位素值分布较宽,并呈现阶梯状变化,而硫化氢的硫同位素则呈现出相似的分布规律,表明各主要含硫化氢气田硫化氢中的硫来自于本层系的硫酸盐,TSR主要发生在各自的储集层中。四川盆地各气田TSR发生的温度条件相似,硫同位素分馏比较接近。TSR过程中硫同位素的分馏过程与硫酸盐本身硫同位素值的高低无关,而与TSR反应程度有关。TSR反应程度越高,硫化氢的硫同位素值与地层硫酸盐的硫同位素越相近。通过系统分析整理全球含硫化氢气田的硫化物硫同位素数据,并结合四川盆地地质条件和油气演化过程,揭示了TSR过程中硫同位素的分馏特征,并绘制出四川盆地和全球各时代硫化氢和石膏的硫同位素分布曲线图,为研究含油气盆地蒸发岩沉积演化和硫化氢成因提供了参考。     

CaO-MgO-Al2o3-SiO2 liquids: chemical and isotopic effects of Mg and Si evaporation in a closed system of solar composition

A method is shown for calculating vapor pressures over a CMAS droplet in a gas of any composition. It is applied to the problem of the evolution of the chemical and Mg and Si isotopic composition of a completely molten droplet having the composition of a likely refractory inclusion precursor during its evaporation into the complementary, i.e. modified solar, gas from which it originally condensed, a more realistic model than previous calculations in which the ambient gas is pure H_2(g). Because the loss rate of Mg is greater than that of Si, the vapor pressure of Mg_(g) falls and its ambient pressure rises faster than those of SiO_(g) during isothermal evaporation, causing the flux of Mg_(g) to approach zero faster and MgO to approach its equilibrium concentration sooner than SiO₂. As time passes, δ²⁵Mg and δ²⁹Si increase in the droplet and decrease in the ambient gas. The net flux of each isotope crossing the droplet/gas interface is the difference between its outgoing and incoming flux. δ²⁵Mg and δ²⁹Si of this instantaneous gas become higher, first overtaking their values in the ambient gas, causing them to increase with time, and later overtaking their values in the droplet itself, causing them to decrease with time, ultimately reaching their equilibrium values. If the system is cooling during evaporation and if mass transfer ceases at the solidus temperature, 1500 K, final MgO and SiO₂ contents of the droplet are slightly higher in modified solar gas than in pure H_2(g), and the difference increases with decreasing cooling rate and increasing ambient pressure. During cooling under some conditions, net fluxes of evaporating species become negative, causing reversal of the evaporation process into a condensation process, an increase in the MgO and/or SiO₂ content of the droplet with time, and an increase in their final concentrations with increasing ambient pressure and/or dust/gas ratio. At cooling rates <∼3 K/h, closed-system evaporation at P^tot ∼ 10⁻³ bar in a modified solar gas, or at lower pressure in systems with enhanced dust/gas ratio, can yield the same δ²⁵Mg in a residual CMAS droplet for vastly different evaporated fractions of Mg. The δ²⁵Mg of a refractory residue may thus be insufficient to determine the extent of Mg loss from its precursor. Evaporation of Mg into an Mg-bearing ambient gas causes δ²⁶Mg and δ²⁵Mg of the residual droplet to fall below values expected from Rayleigh fractionation for the amount of ²⁴Mg evaporated, with the degree of departure increasing with increasing fraction evaporated and ambient pressure of Mg. δ²⁶Mg and δ²⁵Mg do not depart proportionately from Rayleigh fractionation curves, with δ²⁵Mg being less than expected on the basis of δ²⁶Mg by up to ∼1.2‰. Such departures from Rayleigh fractionation could be used in principle to distinguish heavily from lightly evaporated residues with the same δ²⁵Mg.     

Kinetics of H2-O2-H2O redox equilibria and formation of metastable H2O2 under low temperature hydrothermal conditions

Hydrothermal experiments were conducted to evaluate the kinetics of H_2(aq) oxidation in the homogeneous H₂-O₂-H₂O system at conditions reflecting subsurface/near-seafloor hydrothermal environments (55-250 °C and 242-497 bar). The kinetics of the water-forming reaction that controls the fundamental equilibrium between dissolved H_2(aq) and O_2(aq), are expected to impose significant constraints on the redox gradients that develop when mixing occurs between oxygenated seawater and high-temperature anoxic vent fluid at near-seafloor conditions. Experimental data indicate that, indeed, the kinetics of H_2(aq)-O_2(aq) equilibrium become slower with decreasing temperature, allowing excess H_2(aq) to remain in solution. Sluggish reaction rates of H_2(aq) oxidation suggest that active microbial populations in near-seafloor and subsurface environments could potentially utilize both H_2(aq) and O_2(aq), even at temperatures lower than 40 °C due to H_2(aq) persistence in the seawater/vent fluid mixtures. For these H₂-O₂ disequilibrium conditions, redox gradients along the seawater/hydrothermal fluid mixing interface are not sharp and microbially-mediated H_2(aq) oxidation coupled with a lack of other electron acceptors (e.g. nitrate) could provide an important energy source available at low-temperature diffuse flow vent sites.More importantly, when H_2(aq)-O_2(aq) disequilibrium conditions apply, formation of metastable hydrogen peroxide is observed. The yield of H₂O_2(aq) synthesis appears to be enhanced under conditions of elevated H_2(aq)/O_2(aq) molar ratios that correspond to abundant H_2(aq) concentrations. Formation of metastable H₂O₂ is expected to affect the distribution of dissolved organic carbon (DOC) owing to the existence of an additional strong oxidizing agent. Oxidation of magnetite and/or Fe⁺⁺ by hydrogen peroxide could also induce formation of metastable hydroxyl radicals (•OH) through Fenton-type reactions, further broadening the implications of hydrogen peroxide in hydrothermal environments.     

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.     

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.     

CO2 solubility in dacitic melts equilibrated with H2O-CO2 fluids: Implications for modeling the solubility of CO2 in silicic melts

The solubility of CO₂ in dacitic melts equilibrated with H₂O-CO₂ fluids was experimentally investigated at 1250°C and 100 to 500 MPa. CO₂ is dissolved in dacitic glasses as molecular CO₂ and carbonate. The quantification of total CO₂ in the glasses by mid-infrared (MIR) spectroscopy is difficult because the weak carbonate bands at 1430 and 1530 cm⁻¹ can not be reliably separated from background features in the spectra. Furthermore, the ratio of CO_2,mol/carbonate in the quenched glasses strongly decreases with increasing water content. Due to the difficulties in quantifying CO₂ species concentrations from the MIR spectra we have measured total CO₂ contents of dacitic glasses by secondary ion mass spectrometry (SIMS).At all pressures, the dependence of CO₂ solubility in dacitic melts on x^fluid_CO2,total shows a strong positive deviation from linearity with almost constant CO₂ solubility at x_CO2fluid > 0.8 (maximum CO₂ solubility of 795 ± 41, 1376 ± 73 and 2949 ± 166 ppm at 100, 200 and 500 MPa, respectively), indicating that dissolved water strongly enhances the solubility of CO₂. A similar nonlinear variation of CO₂ solubility with x_CO2fluid has been observed for rhyolitic melts in which carbon dioxide is incorporated exclusively as molecular CO₂ (Tamic et al., 2001). We infer that water species in the melt do not only stabilize carbonate groups as has been suggested earlier but also CO₂ molecules.A thermodynamic model describing the dependence of the CO₂ solubility in hydrous rhyolitic and dacitic melts on T, P, f_CO2 and the mol fraction of water in the melt (x_water) has been developed. An exponential variation of the equilibrium constant K₁ with x_water is proposed to account for the nonlinear dependence of x_CO2,totalmelt on x_CO2fluid. The model reproduces the CO₂ solubility data for dacitic melts within ±14% relative and the data for rhyolitic melts within 10% relative in the pressure range 100-500 MPa (except for six outliers at low x_CO2fluid). Data obtained for rhyolitic melts at 75 MPa and 850°C show a stronger deviation from the model, suggesting a change in the solubility behavior of CO₂ at low pressures (a Henrian behavior of the CO₂ solubility is observed at low pressure and low H₂O concentrations in the melt). We recommend to use our model only in the pressure range 100-500 MPa and in the x_CO2fluid range 0.1-0.95. The thermodynamic modeling indicates that the partial molar volume of total CO₂ is much lower in rhyolitic melts (31.7 cm³/mol) than in dacitic melts (46.6 cm³/mol). The dissolution enthalpy for CO₂ in hydrous rhyolitic melts was found to be negligible. This result suggests that temperature is of minor importance for CO₂ solubility in silicic melts.