Diffuse gas emissions at the Ukinrek Maars,Alaska: Implications for magmatic degassing and volcanic monitoring

Diffuse CO₂ efflux near the Ukinrek Maars, two small volcanic craters that formed in 1977 in a remote part of the Alaska Peninsula, was investigated using accumulation chamber measurements. High CO₂ efflux, in many places exceeding 1000 g m⁻² d⁻¹, was found in conspicuous zones of plant damage or kill that cover 30,000–50,000 m² in area. Total diffuse CO₂ emission was estimated at 21–44 t d⁻¹. Gas vents 3-km away at The Gas Rocks produce 0.5 t d⁻¹ of CO₂ that probably derives from the Ukinrek Maars basalt based on similar δ¹³C values (∼−6‰), ³He/⁴He ratios (5.9–7.2 R_A), and CO₂/³He ratios (1–2 × 10⁹) in the two areas. A lower ³He/⁴He ratio (2.7 R_A) and much higher CO₂/³He ratio (9 × 10¹⁰) in gas from the nearest arc-front volcanic center (Mount Peulik/Ugashik) provide a useful comparison. The large diffuse CO₂ emission at Ukinrek has important implications for magmatic degassing, subsurface gas transport, and local toxicity hazards. Gas–water–rock interactions play a major role in the location, magnitude and chemistry of the emissions.     

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

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.     

The response of soil CO<Subscript>2</Subscript> efflux to desertification on alpine meadow in the Qinghai–Tibet Plateau

Assessing the global C budget requires a better understanding of the effect of environmental factors on soil CO₂ efflux from both experiments and theoretical research, especially in different desertified lands in the Qinghai–Tibet Plateau. Based on the enclosed chamber method, soil CO₂ efflux in four different desertified lands and one control [alpine meadow (AM)] were measured in June, August and September, 2008, respectively. Soil CO₂ efflux rates at the top, the middle, the bottom of a slope and the flat in front of the slope were obtained at Maduo County. The results showed that average daily soil CO₂ efflux rates were 3.72, 2.65, 2.68, 0.59 and 0.37 g m⁻² day⁻¹ in the AM, lightly (LDL), moderately (MDL), severely (SDL) and very severely desertified lands (VSDL) during the growing season, respectively. Soil CO₂ efflux decreased with the change of desertification. The response of soil CO₂ efflux to environmental factors was adequately described by the linear model; models accounted for 76, 65, 72, 59 and 71% of the variability on soil CO₂ efflux in the AM, LDL, MDL, SDL and VSDL, respectively. Any environmental factor, however, was insufficient to explain the soil CO₂ efflux; the common influence could perfectly reflect soil CO₂ efflux response to the desertification change.     

Stream geochemistry, chemical weathering and CO2 consumption potential of andesitic terrains, Dominica, Lesser Antilles

Recent studies of chemical weathering of andesitic-dacitic material on high-standing islands (HSIs) have shown these terrains have some of the highest observed rates of chemical weathering and associated CO₂ consumption yet reported. However, the paucity of stream gauge data in many of these terrains has limited determination of chemical weathering product fluxes. In July 2006 and March 2008, stream water samples were collected and manual stream gauging was performed in watersheds throughout the volcanic island of Dominica in the Lesser Antilles. Distinct wet and dry season solute concentrations reveal the importance of seasonal variations on the weathering signal. A cluster analysis of the stream geochemical data shows the importance of parent material age on the overall delivery of solutes. Observed Ca:Na, HCO₃:Na and Mg:Na ratios suggest crystallinity of the parent material may also play an important role in determining weathering fluxes. From total dissolved solids concentrations and mean annual discharge calculations we calculate chemical weathering yields of (6-106 t km⁻² a⁻¹), which are similar to those previously determined for basalt terrains. Silicate fluxes (3.1-55.4 t km⁻² a⁻¹) and associated CO₂ consumption (190-1575 × 10³ mol km⁻² a⁻¹) determined from our study are among the highest determined to date. The calculated chemical fluxes from our study confirm the weathering potential of andesitic-dacitic terrains and that additional studies of these terrains are warranted.     

Quantifying inorganic sources of geochemical energy in hydrothermal ecosystems, Yellowstone National Park, USA

Combining analytical data from hot spring samples with thermodynamic calculations permits a quantitative assessment of the availability and ranking of various potential sources of inorganic chemical energy that may support microbial life in hydrothermal ecosystems. Yellowstone hot springs of diverse geochemical composition, ranging in pH from <2 to >9 were chosen for this study, and dozens of samples were collected during three field seasons. Field measurements of dissolved oxygen, nitrate, nitrite, total ammonia, total sulfide, alkalinity, and ferrous iron were combined with laboratory analyses of sulfate and other major ions from water samples, and carbon dioxide, hydrogen, methane, and carbon monoxide in gas samples to evaluate activity products for ∼300 coupled oxidation-reduction reactions. Comparison of activity products and independently calculated equilibrium constants leads to values of the chemical affinity for each of the reactions, which quantifies how far each reaction is from equilibrium. Affinities, in turn, show systematic behavior that is independent of temperature but can be correlated with pH of the hot springs as a proxy for the full spectrum of geochemical variability. Many affinities are slightly to somewhat dependent on pH, while others are dramatically influenced by changes in chemical composition. All reactions involving dissolved oxygen as the electron acceptor are potential energy sources in all hot spring samples collected, but the ranking of dominant electron donors changes from ferrous iron, and sulfur at high pH to carbon monoxide, hydrogen, and magnetite as pH decreases. There is a general trend of decreasing energy yields depending on electron acceptors that follows the sequence: O₂(aq) > NO₃⁻ ≈ NO₂⁻ ≈ S > pyrite ≈ SO₄⁻² ≈ CO(g) ≈ CO₂(g) at high pH, and O₂(aq) ≈ magnetite > hematite ≈ goethite > NO₃⁻ ≈ NO₂⁻ ≈ S ≈ pyrite ≈ SO₄⁻² at low pH. Many reactions that are favorable sources of chemical energy at one set of geochemical conditions fail to provide energy at other conditions, and vice versa. This results in energy profiles supplied by geochemical processes that provide fundamentally different foundations for chemotrophic microbial communities as composition changes.     

Stable isotopic evidence in support of active microbial methane cycling in low-temperature diffuse flow vents at 9°50′N East Pacific Rise

A unique dataset from paired low- and high-temperature vents at 9°50′N East Pacific Rise provides insight into the microbiological activity in low-temperature diffuse fluids. The stable carbon isotopic composition of CH₄ and CO₂ in 9°50′N hydrothermal fluids indicates microbial methane production, perhaps coupled with microbial methane consumption. Diffuse fluids are depleted in ¹³C by ∼10‰ in values of δ¹³C of CH₄, and by ∼0.55‰ in values of δ¹³C of CO₂, relative to the values of the high-temperature source fluid (δ¹³C of CH₄ =−20.1 ± 1.2‰, δ¹³C of CO₂ =−4.08 ± 0.15‰). Mixing of seawater or thermogenic sources cannot account for the depletions in ¹³C of both CH₄ and CO₂ at diffuse vents relative to adjacent high-temperature vents. The substrate utilization and ¹³C fractionation associated with the microbiological processes of methanogenesis and methane oxidation can explain observed steady-state CH₄ and CO₂ concentrations and carbon isotopic compositions. A mass-isotope numerical box model of these paired vent systems is consistent with the hypothesis that microbial methane cycling is active at diffuse vents at 9°50′N. The detectable ¹³C modification of fluid geochemistry by microbial metabolisms may provide a useful tool for detecting active methanogenesis.     

Methane microseepage from different sectors of the Yakela condensed gas field in Tarim Basin,Xinjiang, China

Methane microseepage is the result of natural gas migration from subsurface hydrocarbon accumulations to the Earth’s surface, and it is quite common in commercial petroleum fields. While the role of microseepage as a pathfinder in petroleum exploration has been known for about 80 a, its significance as an atmospheric CH₄ source has only recently been studied, and flux data are currently available only in the USA and Europe. With the aim of increasing the global data-set and better understanding flux magnitudes and variabilities, microseepage is now being extensively studied in China. A static flux chamber method was recently applied to study microseepage emissions into the atmosphere in four different sectors of the Yakela condensed gas field in Tarim Basin, Xinjiang, China, and specifically in: (a) a faulted sector, across the Luntai fault systems; (b) an oil–water interface sector, at the northern margin of the field; (c) an oil–gas interface sector, in the middle of the field; (d) an external area, outside the northern gas field boundaries. The results show that positive CH₄ fluxes are pervasive in all sectors and therefore, only part of the CH₄ migrating from the deep oil–gas reservoirs is consumed in the soil by methanotrophic oxidation. The intensity of gas seepage seems to be controlled by subsurface geologic settings and lateral variabilities of natural gas pressure in the condensed gas field. The highest CH₄ fluxes, up to ∼14 mg m⁻² d⁻¹ (mean of 7.55 mg m⁻² d⁻¹) with higher spatial variability (standard deviation, σ: 2.58 mg m⁻² d⁻¹), occur in the Luntai fault sector. Merhane flux was lower in the oil–water area (mean of 0.53 mg m⁻² d⁻¹) and the external area (mean of 1.55 mg m⁻² d⁻¹), and at the intermediate level in the gas–oil sector (mean of 2.89 mg m⁻² d⁻¹). These values are consistent with microseepage data reported for petroleum basins in the USA and Europe. The build-up of methane concentration in the flux chambers is always coupled with an enrichment of ¹³C, from δ¹³C₁ of −46‰ to −42.5‰ (VPDB), which demonstrates that seeping methane is thermogenic, as that occurring in the deep Yakela reservoir. Daily variations of microseepage are very low, with minima in the afternoon, corresponding to higher soil temperature (and higher methanotrophic consumption), and maxima in the early morning (when soil temperatures are lowest). A preliminary and rough estimate of the total amount of CH₄ exhaled from the Yakela field is in the order of 10² tonnes a⁻¹.     

The fluvial geochemistry, contributions of silicate, carbonate and saline-alkaline components to chemical weathering flux and controlling parameters: Narmada River (Deccan Traps), India

The Narmada River in India is the largest west-flowing river into the Arabian Sea, draining through the Deccan Traps, one of the largest flood basalt provinces in the world. The fluvial geochemical characteristics and chemical weathering rates (CWR) for the mainstream and its major tributaries were determined using a composite dataset, which includes four phases of seasonal field (spot) samples (during 2003 and 2004) and a decade-long (1990-2000) fortnight time series (multiannual) data. Here, we demonstrate the influence of minor lithologies (carbonates and saline-alkaline soils) on basaltic signature, as reflected in sudden increases of Ca²⁺-Mg²⁺ and Na⁺ contents at many locations along the mainstream and in tributaries. Both spot and multiannual data corrected for non-geological contributions were used to calculate the CWR. The CWR for spot samples (CWR_spot) vary between 25 and 63 ton km⁻² year⁻¹, showing a reasonable correspondence with the CWR estimated for multiannual data (CWR_multi) at most study locations. The weathering rates of silicate (SilWR), carbonate (CarbWR) and evaporite (Sal-AlkWR) have contributed ∼38-58, 28-45 and 8-23%, respectively to the CWR_spot at different locations. The estimated SilWR (11-36 ton km⁻² year⁻¹) for the Narmada basin indicates that the previous studies on the North Deccan Rivers (Narmada-Tapti-Godavari) overestimated the silicate weathering rates and associated CO₂ consumption rates. The average annual CO₂ drawdown via silicate weathering calculated for the Narmada basin is ∼0.032 × 10¹² moles year⁻¹, suggesting that chemical weathering of the entire Deccan Trap basalts consumes approximately 2% (∼0.24 × 10¹² moles) of the annual global CO₂ drawdown. The present study also evaluates the influence of meteorological parameters (runoff and temperature) and physical weathering rates (PWR) in controlling the CWR at annual scale across the basin. The CWR and the SilWR show significant correlation with runoff and PWR. On the basis of observed wide temporal variations in the CWR and their close association with runoff, temperature and physical erosion, we propose that the CWR in the Narmada basin strongly depend on meteorological variability. At most locations, the total denudation rates (TDR) are dominated by physical erosion, whereas chemical weathering constitutes only a small part (<10%). Thus, the CWR to PWR ratio for the Narmada basin can be compared with high relief small river watersheds of Taiwan and New Zealand (1-5%) and large Himalayan Rivers such as the Brahmaputra and the Ganges (8-9%).     

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.     

Carbon dioxide degassing and thermal energy release in the Monte Amiata volcanic-geothermal area (Italy)

The quaternary volcanic complex of Mount Amiata is located in southern Tuscany (Italy) and represents the most recent manifestation of the Tuscan Magmatic Province. The region is characterised by a large thermal anomaly and by the presence of numerous CO₂-rich gas emissions and geothermal features, mainly located at the periphery of the volcanic complex. Two geothermal systems are located, at increasing depths, in the carbonate and metamorphic formations beneath the volcanic complex. The shallow volcanic aquifer is separated from the deep geothermal systems by a low permeability unit (Ligurian Unit). A measured CO₂ discharge through soils of 1.8 × 10⁹ mol a⁻¹ shows that large amounts of CO₂ move from the deep reservoir to the surface. A large range in δ¹³C_TDIC (−21.07 to +3.65) characterises the waters circulating in the aquifers of the region and the mass and isotopic balance of TDIC allows distinguishing a discharge of 0.3 × 10⁹ mol a⁻¹ of deeply sourced CO₂ in spring waters. The total natural CO₂ discharge (2.1 × 10⁹ mol a⁻¹) is slightly less than minimum CO₂ output estimated by an indirect method (2.8 × 10⁹ mol a⁻¹), but present-day release of 5.8 × 10⁹ mol a⁻¹ CO₂ from deep geothermal wells may have reduced natural CO₂ discharge. The heat transported by groundwater, computed considering the increase in temperature from the infiltration area to the discharge from springs, is of the same order of magnitude, or higher, than the regional conductive heat flow (>200 mW m⁻²) and reaches extremely high values (up to 2700 mW m⁻²) in the north-eastern part of the study area. Heat transfer occurs mainly by conductive heating in the volcanic aquifer and by uprising gas and vapor along fault zones and in those areas where low permeability cover is lacking. The comparison of CO₂ flux, heat flow and geological setting shows that near surface geology and hydrogeological setting play a central role in determining CO₂ degassing and heat transfer patterns.     

Can pelagic net heterotrophy account for carbon fluxes from eastern Canadian lakes?

Lakes worldwide are commonly oversaturated with CO₂, however the source of this CO₂ oversaturation is not well understood. To examine the magnitude of the C flux to the atmosphere and determine if an excess of respiration (R) over gross primary production (GPP) is sufficient to account for this C flux, metabolic parameters and stable isotopes of dissolved O₂ and C were measured in 23 Québec lakes. All of the lakes sampled were oversaturated with CO₂ over the sampling period, on average 221 ± 25%. However, little evidence was found to conclude that this CO₂ oversaturation was the result of an excess of pelagic R over GPP. In lakes Croche and à l’Ours, where CO₂ flux, R and GPP were measured weekly, the annual difference between pelagic GPP and R, or net primary production (NPP), was not sufficient to account for the size of the CO₂ flux to the atmosphere. In Lac Croche average annual NPP was 14.4 mg C m⁻² d⁻¹ while the average annual flux of CO₂ to the atmosphere was 34 mg C m⁻² d⁻¹. In Lac à l’Ours average annual NPP was −9.1 mg C m⁻² d⁻¹ while the average annual flux of CO₂ to the atmosphere was 55 mg C m⁻² d⁻¹. In all of the lakes sampled, O₂ saturation averaged 104.0 ± 1.7% during the ice-free season and the isotopic composition of dissolved O₂ (δ¹⁸O_DO) was 22.9 ± 0.3‰, lower than atmospheric values and indicative of net autotrophy. Carbon evasion was not a function of R, nor did the isotopic signature of dissolved CO₂ in the lakes present evidence of excess R over GPP. External inputs of C must therefore subsidize the lake to explain the continued CO₂ oversaturation. The isotopic composition of dissolved inorganic C (δ¹³C_DIC) indicates that the CO₂ oversaturation cannot be attributed to in situ aerobic respiration. δ¹³C_DIC reveals a source of excess C enriched in ¹³C, which may be accounted for by anaerobic sediment respiration or groundwater inputs followed by kinetic isotope fractionation during degassing under open system conditions.     

Raman spectroscopic measurements of CO2 density: Experimental calibration with high-pressure optical cell (HPOC) and fused silica capillary capsule (FSCC) with application to fluid inclusion observations

Raman spectroscopy is a powerful method for the determination of CO₂ densities in fluid inclusions, especially for those with small size and/or low fluid density. The relationship between CO₂ Fermi diad split (Δ, cm⁻¹) and CO₂ density (ρ, g/cm³) has been documented by several previous studies. However, significant discrepancies exist among these studies mainly because of inconsistent calibration procedures and lack of measurements for CO₂ fluids having densities between 0.21 and 0.75 g/cm³, where liquid and vapor phases coexist near room temperature.In this study, a high-pressure optical cell and fused silica capillary capsules were used to prepare pure CO₂ samples with densities between 0.0472 and 1.0060 g/cm³. The measured CO₂ Fermi diad splits were calibrated with two well established Raman bands of benzonitrile at 1192.6 and 1598.9 cm⁻¹. The relationship between the CO₂ Fermi diad split and density can be represented by: ρ = 47513.64243 − 1374.824414 × Δ + 13.25586152 × Δ² − 0.04258891551 × Δ³ (r² = 0.99835, σ = 0.0253 g/cm³), and this relationship was tested by synthetic fluid inclusions and natural CO₂-rich fluid inclusions. The effects of temperature and the presence of H₂O and CH₄ on this relationship were also examined.     

Dissolution kinetics of diopside as a function of solution saturation state: Macroscopic measurements and implications for modeling of geological storage of CO2

Measurements of the dissolution rate of diopside (r) were carried out as a function of the Gibbs free energy of the dissolution reaction (ΔG_r) in a continuously stirred flow-through reactor at 90 °C and pH_90 °C = 5.05. The overall relation between r and ΔG_r was determined over a free energy range of −130.9 < ΔG_r < −47.0 kJ mo1⁻¹. The data define a highly non-linear, sigmoidal relation between r and ΔG_r. At far-from-equilibrium conditions (ΔG_r ? −76.2 kJ mo1⁻¹), a rate plateau is observed. In this free energy range, the rates of dissolution are constant, independent of [Ca], [Mg] and [Si] concentrations, and independent of ΔG_r. A sharp decrease of the dissolution rate (∼1 order of magnitude) occurs in the transition ΔG_r region defined by −76.2 < ΔG_r ? −61.5 kJ mo1⁻¹. Dissolution closer to equilibrium (ΔG_r > −61.5 kJ mo1⁻¹) is characterised by a much weaker inverse dependence of the rates on ΔG_r. Modeling the experimental r-ΔG_r data with a simple classical transition state theory (TST) law as implemented in most available geochemical codes is found inappropriate. An evaluation of the consequences of the use of geochemical codes where the r-ΔG_r relation is based on basic TST was carried out and applied to carbonation reactions of diopside, which, among other reactions with Ca- and Mg-bearing minerals, are considered as a promising process for the solid state sequestration of CO₂ over long time spans. In order to take into account the actual experimental r-ΔG_r relation in the geochemical code that we used, a new module has been developed. It reveals a dramatic overestimation of the carbonation rate when using a TST-based geochemical code. This points out that simulations of water-rock-CO₂ interactions performed with classical geochemical codes should be evaluated with great caution.     

Pathways and regulation of carbon, sulfur and energy transfer in marine sediments overlying methane gas hydrates on the Opouawe Bank (New Zealand)

This study combines sediment geochemical analysis, in situ benthic lander deployments and numerical modeling to quantify the biogeochemical cycles of carbon and sulfur and the associated rates of Gibbs energy production at a novel methane seep. The benthic ecosystem is dominated by a dense population of tube-building ampharetid polychaetes and conspicuous microbial mats were unusually absent. A 1D numerical reaction-transport model, which allows for the explicit growth of sulfide and methane oxidizing microorganisms, was tuned to the geochemical data using a fluid advection velocity of 14 cm yr⁻¹. The fluids provide a deep source of dissolved hydrogen sulfide and methane to the sediment with fluxes equal to 4.1 and 18.2 mmol m⁻² d⁻¹, respectively. Chemosynthetic biomass production in the subsurface sediment is estimated to be 2.8 mmol m⁻² d⁻¹ of C biomass. However, carbon and oxygen budgets indicate that chemosynthetic organisms living directly above or on the surface sediment have the potential to produce 12.3 mmol m⁻² d⁻¹ of C biomass. This autochthonous carbon source meets the ampharetid respiratory carbon demand of 23.2 mmol m⁻² d⁻¹ to within a factor of 2. By contrast, the contribution of photosynthetically-fixed carbon sources to ampharetid nutrition is minor (3.3 mmol m⁻² d⁻¹ of C). The data strongly suggest that mixing of labile autochthonous microbial detritus below the oxic layer sustains high measured rates of sulfate reduction in the uppermost 2 cm of the sulfidic sediment (100-200 nmol cm⁻³ d⁻¹). Similar rates have been reported in the literature for other seeps, from which we conclude that autochthonous organic matter is an important substrate for sulfate reducing bacteria in these sediment layers. A system-scale energy budget based on the chemosynthetic reaction pathways reveals that up to 8.3 kJ m⁻² d⁻¹ or 96 mW m⁻² of catabolic (Gibbs) energy is dissipated at the seep through oxidation reactions. The microorganisms mediating sulfide oxidation and anaerobic oxidation of methane (AOM) produce 95% and 2% of this energy flux, respectively. The low power output by AOM is due to strong bioenergetic constraints imposed on the reaction rate by the composition of the chemical environment. These constraints provide a high potential for dissolved methane efflux from the sediment (12.0 mmol m⁻² d⁻¹) and indicates a much lower efficiency of (dissolved) methane sequestration by AOM at seeps than considered previously. Nonetheless, AOM is able to consume a third of the ascending methane flux (5.9 mmol m⁻² d⁻¹ of CH₄) with a high efficiency of energy expenditure (35 mmol CH₄ kJ⁻¹). It is further proposed that bioenergetic limitation of AOM provides an explanation for the non-zero sulfate concentrations below the AOM zone observed here and in other active and passive margin sediments.     

The role of sulfate groups in controlling CaCO3 polymorphism

The nucleation and growth of CaCO₃ phases from aqueous solutions with SO₄²⁻:CO₃²⁻ ratios from 0 to 1.62 and a pH of ∼10.9 were studied experimentally in batch reactors at 25 °C. The mineralogy, morphology and composition of the precipitates were characterized by X-ray diffraction, Fourier transform infrared spectroscopy, scanning electron microscopy and microanalyses. The solids recovered after short reaction times (5 min to 1 h) consisted of a mixture of calcite and vaterite, with a S content that linearly correlates with the SO₄²⁻:CO₃²⁻ ratio in the aqueous solution. The solvent-mediated transformation of vaterite to calcite subsequently occurred. After 24 h of equilibration, calcite was the only phase present in the precipitate for aqueous solutions with SO₄²⁻:CO₃²⁻ ? 1. For SO₄²⁻:CO₃²⁻ > 1, vaterite persisted as a major phase for a longer time (>250 h for SO₄²⁻:CO₃²⁻ = 1.62). To study the role of sulfate in stabilizing vaterite, we performed a molecular simulation of the substitution of sulfate for carbonate groups into the crystal structure of vaterite, aragonite and calcite. The results obtained show that the incorporation of small amounts (<3 mole%) of sulfate is energetically favorable in the vaterite structure, unfavorable in calcite and very unfavorable in aragonite. The computer modeling provided thermodynamic information, which, combined with kinetic arguments, allowed us to put forward a plausible explanation for the observed crystallization behavior.     

Major ion chemistry of the Yarlung Tsangpo-Brahmaputra river: Chemical weathering, erosion, and CO2 consumption in the southern Tibetan plateau and eastern syntaxis of the Himalaya

The Yarlung Tsangpo-Brahmaputra river drains a large portion of the Himalaya and southern Tibetan plateau, including the eastern Himalayan syntaxis, one of the most tectonically active regions on the globe. We measured the solute chemistry of 161 streams and major tributaries of the Tsangpo-Brahmaputra to examine the effect of tectonic, climatic, and geologic factors on chemical weathering rates. Specifically, we quantify chemical weathering fluxes and CO₂ consumption by silicate weathering in southern Tibet and the eastern syntaxis of the Himalaya, examine the major chemical weathering reactions in the tributaries of the Tsangpo-Brahmaputra, and determine the total weathering flux from carbonate and silicate weathering processes in this region. We show that high precipitation, rapid tectonic uplift, steep channel slopes, and high stream power generate high rates of chemical weathering in the eastern syntaxis. The total dissolved solids (TDS) flux from the this area is greater than 520 tons km⁻² yr⁻¹ and the silicate cation flux more than 34 tons km⁻² yr⁻¹. In total, chemical weathering in this area consumes 15.2 × 10⁵ mol CO₂ km⁻² yr⁻¹, which is twice the Brahmaputra average. These data show that 15-20% of the total CO₂ consumption by silicate weathering in the Brahmaputra catchment is derived from only 4% of the total land area of the basin. Hot springs and evaporite weathering provide significant contributions to dissolved Na⁺ and Cl⁻ fluxes throughout southern Tibet, comprising more than 50% of all Na⁺ in some stream systems. Carbonate weathering generates 80-90% of all dissolved Ca²⁺ and Mg²⁺ cations in much of the Yarlung Tsangpo catchment.     

Regional and temporal variations in CO2/He,He/He and δC along the North Anatolian Fault Zone,Turkey

New He and C relative abundance, isotope and concentration results from nine geothermal locations situated along an 800-km transect of the North Anatolian Fault Zone (NAFZ), Turkey, that were monitored during the period November 2001–November 2004, are reported. The geothermal waters were collected every 3–6 months to study possible links between temporal geochemical variations and seismic activity along the NAFZ. At the nine sample locations, the He isotope ratios range from 0.24 to 2.3R_A, δ¹³C values range from −4.5 to +5.8‰, and CO₂/³He ratios range from 5 × 10⁹ to 5 × 10¹⁴. The following geochemical observations are noted: (1) the highest ³He/⁴He ratios are found near the Galatean volcanic region, in the central section of the NAFZ, (2) at each of the nine sample locations, the ³He/⁴He ratios are generally constant; however, CO₂/³He ratios and He contents both show one order of magnitude variability, and δ¹³C values show up to ∼4‰ variability, and (3) at all locations (except Re?adiye), δ¹³C values show positive correlations with CO₂ contents. The results indicate that at least three processes are necessary to account for the geochemical variations: (1) binary mixing between crustal and mantle-derived volatiles can explain the general characteristics of ³He/⁴He ratios, δ¹³C values, and CO₂/³He ratios at the nine sample locations; (2) preferential degassing of He from the geothermal waters is responsible for variations in CO₂/³He values and He contents at each sample location; and (3) CO₂ dissolution followed by calcite precipitation is responsible for variations in CO₂ contents and δ¹³C values at most locations. For each of the geochemical parameters, anomalies are defined in the temporal record by values that fall outside two standard deviations of average values at each specific location. Geochemical anomalies that may be related to seismic activity are recorded on June 28, 2004 at Yalova, where a M = 4.2 earthquake occurred 43 days earlier at 15 km distance from the sample location, and on April 7, 2003 at Efteni, where a M = 4.0 earthquake occurred 44 days later at a distance of 12 km. At both locations, the sampling periods containing geochemical anomalies were preceded by an increase in M ? 3 earthquakes occurring within 60 days and less than 40 km distance.     

Carbon dioxide and argon diffusion in silicate melts: Insights into the CO2 speciation in magmas

To investigate the influence of temperature and composition on the diffusivities of dissolved carbon dioxide and argon in silicate melts, diffusion experiments were performed at magmatic pressure and temperature conditions in (a) albite melts with excess Na₂O (0-8.6 wt%) and a constant Si/Al ratio of 3, and (b) albite₇₀quartz₃₀ to jadeite melts with decreasing SiO₂ content and a constant Na/Al ratio of 1. We obtained diffusion coefficients at 500 MPa and 1323-1673 K. In the fully polymerized system Ab₇₀Qz₃₀ - Jd, the change in composition only has a weak effect on bulk CO₂ diffusivity, but Ar diffusivity increases clearly with decreasing SiO₂ content. In the system Ab + Na₂O, bulk CO₂ and Ar diffusivity increase significantly with gradual depolymerisation. The relatively small change in composition on molar basis in the depolymerized system leads to a significantly larger change in diffusivities compared to the fully polymerized Ab₇₀Qz₃₀-Jd join. Within error, activation energies for bulk CO₂ and Ar diffusion in both systems are identical with decreasing silica content (Ab + Na₂O: 159 ± 25 kJ mol⁻¹ for bulk CO₂ and 130 ± 8 kJ mol⁻¹ for Ar; Ab₇₀Qz₃₀-Jd: 163 ± 16 kJ mol⁻¹ for bulk CO₂ and 148 ± 15 kJ mol⁻¹ for Ar) even though this results in depolymerisation in one system and not the other.Although there is a variation in CO₂ speciation with changing composition as observed in quenched glasses, it has previously established that this is not a true representation of the species present in the melt, with the ratio of molecular CO₂ to carbonate decreasing during quenching. Thus, diffusion coefficients for the individual CO₂ species cannot be directly derived by measuring molecular CO₂ and CO₃^2- concentration-distance profiles in the glasses. To obtain diffusivities of individual CO₂ species, we have made two assumptions that (1) inert Ar can be used as a proxy for molecular CO₂ diffusion characteristics as shown by our previous work and (2) the diffusivity of CO₃²⁻ can be calculated assuming it is identical to network forming components (Si⁴⁺ and Al³⁺). This is derived from viscosity data (Eyring eqn.) and suggests that CO₃²⁻ diffusion would be several orders of magnitude slower than molecular CO₂ diffusion.The systematics of measured bulk CO₂ diffusivity rates and comparison with the Ar proxy all suggest that the faster molecular CO₂ species is much more dominant in melts than measurements on resulting quenched glasses would suggest. This study has confirmed an observation of surprisingly consistent bulk CO₂ diffusivity across a range of natural compositions were Ar diffusivity significantly increases. This is consistent with an actual increase in molecular CO₂ mobility (similar to Ar) that is combined with an increase in the proportion of the slower carbonate in the melt.These results demonstrate that the CO₂ diffusion and speciation model provides an insight into the transport processes in the melt and is promising and an alternative tool to in situ speciation measurements at magmatic conditions, which at the moment are technically extremely difficult. We present the first high pressure high temperature in situ MIR spectra of a CO₂ bearing albitic glass/melt suggesting that molecular CO₂ is a stable species at high temperature, which is qualitatively consistent with the modelled CO₂ speciation data.     

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.