Influences of air CO<Subscript>2</Subscript> on hydrochemistry of drip water and implications for paleoclimate study in a stream-developed cave,SW China

Cave air CO₂ is a vital part of the cave environment. Most studies about cave air CO₂ variations are performed in caves with no streams; there are few studies to date regarding the relationship of cave air CO₂ variations and drip water hydrochemistry in underground stream–developed caves. To study the relationship of underground stream, drip water, and cave air CO₂, monthly and daily monitoring of air CO₂ and of underground stream and drip water was performed in Xueyu Cave from 2012 to 2013. The results revealed that there was marked seasonal variation of air CO₂ and stream hydrochemistry in the cave. Daily variations of cave air CO₂, and of stream and drip water hydrochemistry, were notable during continuous monitoring. A dilution effect was observed by analyzing hydrochemical variations in underground stream and drip water after rainfall. High cave air CO₂ along with low pH and low δ¹³C_DIC in stream and drip water indicated that air CO₂ was one of the dominant factors controlling stream and drip water hydrochemistry on a daily scale. On a seasonal scale, stream flows may promote increased cave air CO₂ in summer; in turn, the higher cave air CO₂ could inhibit degassing of drip water and make calcite δ¹³C more negative. Variation of calcite δ¹³C (precipitated from drip water) was in reverse of monthly temperature, soil CO₂, and cave air CO₂. Therefore, calcite δ¹³C in Xueyu Cave could be used to determine monthly changes outside the cave. However, considering the different precipitation rate of sediment in different seasons, it was difficult to use stalagmites to reconstruct environmental change on a seasonal scale.     

Carbon dioxide in Postojna Cave (Slovenia): spatial distribution,seasonal dynamics and evaluation of plausible sources and sinks

The CO₂ concentration of the air in Postojna Cave (400–7900 ppm) is found to be induced by CO₂ sources (human respiration contributing?~?20,000–58,000 ppm per breath, outgassing of dripwater and water seeping from the vadose zone/epikarst with a pCO₂ values of 5000–29,000 ppm, and underground Pivka River having pCO₂ at 2344–4266 ppm) and CO₂ dilution (inflow of outside air with a CO₂ concentration of?~?400 ppm). Measurements show that sinking Pivka River has the lowest CO₂ concentration among plausible CO₂ sources but still continuously exceeds the surrounding cave air CO₂ concentration. During the winter months, intensive ventilation reduces the cave air CO₂ concentration to outside levels (~?400 ppm), even in the centre of the cave system. CO₂ dilution is less pronounced in summer (CO_2(min)?≈?800 ppm), since the ventilation rate is not as strong as in winter and the outside air that enters the cave through breathing holes and fractures is enriched with soil CO₂. During spring and autumn, the daily alternation of the ventilation regime with a smaller rate of air exchange results in yearly cave air CO₂ peaks of up to?~?2400 ppm. Some dead-end passages can be much less affected by ventilation, resulting in a cave air CO₂ concentration of up to 7900 ppm. The strongest diurnal CO₂ peaks due to human respiration were recorded during the spring holidays (increase of up to 1300 ppm day^?1), compared to considerably smaller summer peaks despite peak visits (increase of?~?600 ppm day^?1).     

Intra‐annual variation of the calcite deposition rate of drip water in Shihua Cave,Beijing, China and its implications for palaeoclimatic reconstructions

Cai, B., Zhu, J., Ban, F. & Tan, M. 2011: Intra‐annual variation of the calcite deposition rate of drip water in Shihua Cave, Beijing, China and its implications for palaeoclimatic reconstructions. Boreas, Vol. 40, pp. 525–535. 10.1111/j.1502‐3885.2010.00201.x. ISSN 0300‐9483. Monthly in situ monitoring of the calcite deposition rate, drip‐water chemistry and surrounding cave environment was carried out at Shihua Cave, Beijing, China, through two hydrological years (from January 2006 to February 2008) to determine the seasonal variability and mechanisms of stalagmite growth in Shihua Cave. Calcite deposition rates exhibit significant intra‐annual variation, with the lowest values during the summer monsoonal rainy season (July–August) and peak values from autumn to spring. The temporal change in the calcite deposition rate is negatively correlated with the drip rate, cave‐air PCO₂ (CO₂ partial pressure) and Ca concentration, and positively correlated with the pH of the feeding drip water. The seasonal recharge regime of drip water is likely to be the primary control on the drip‐water quality and quantity, which, in turn, control the calcite deposition rate in Shihua Cave. During the summer rainy season, periodic and intense rainstorms increase the drip rate and cave‐air PCO₂, leading to drip water with a lower pH and saturation index of calcite, thereby reducing the calcite precipitation. It seems that the high cave‐air PCO₂ is the dominant control on the calcite deposition rate during the rainy season. Our previous study on the dissolved organic carbon of drip water concluded that the thin luminescent bands in stalagmite laminae from Shihua Cave form during the rainy season. The lower calcite deposition rate during the rainy season further supports this suggestion. The significant intra‐seasonal variability of the calcite deposition rate implies that the seasonal bias of δ¹⁸O of stalagmites should be considered when stalagmite δ¹⁸O is used as a high‐resolution palaeoclimatic archive.     

Cave air control on dripwater geochemistry, Obir Caves (Austria): Implications for speleothem deposition in dynamically ventilated caves

There are very few process studies that demonstrate the annual variation in cave environments depositing speleothems. Accordingly, we initiated a monitoring program at the Obir Caves, an Austrian dripstone cave system characterized by a seasonally changing air flow that results in a predictable pattern of high pCO₂ during summer and low pCO₂ in winter. Although similar seasonal changes in soil pCO₂ occur, they are not directly connected with the changes in the subsurface since the dripwaters are fed from a well-mixed source showing little seasonal variation. Cold season flushing by relatively CO₂-poor air enhances degassing of CO₂ in the cave and leads to a high degree of supersaturation of dripwater with regard to calcite. Forced calcite deposition during the cold season also gives rise to a pronounced pattern of synchronous seasonal variations in electrical conductivity, alkalinity, pH, Ca and δ¹³C_DIC which parallel variations recorded in δ¹³C_{cave air}. Chemical components unaffected by calcite precipitation (e.g., δD, δ¹⁸O, SiO₂, SO₄) lack a seasonal signal attesting to a long residence in the karst aquifer. Modeling shows that degassing of CO₂ from seepage waters results in kinetically-enhanced C isotopic fractionation, which contrasts with the equilibrium degassing shown from the Soreq cave in Israel. The Obir Caves may serve as a case example of a dripstone cave whose seepage waters (and speleothems) show intra-annual geochemical variability that is primarily due to chemical modification of the groundwater by a dynamic, bidirectional subsurface air circulation.     

Cave air and hydrological controls on prior calcite precipitation and stalagmite growth rates: Implications for palaeoclimate reconstructions using speleothems

Hourly resolved cave air P_CO2 and cave drip water hydrochemical data illustrate that calcite deposition on stalagmites can be modulated by prior calcite precipitation (PCP) on extremely short timescales. A very clear second-order covariation between cave air P_CO2 and drip water Ca²⁺ concentrations during the winter months demonstrates the effects of degassing-induced PCP on drip water chemistry. Estimating the strength of the cave air P_CO2 control on PCP is possible because the PCP signal is so clear; at our drip site a one ppm shift in Ca²⁺ concentrations requires a P_CO2 shift of between 333 and 667 ppm. This value will undoubtedly vary from site to site, depending on drip water flow rate, residence time, drip water-cave air P_CO2 differential, and availability of low P_CO2 void spaces in the vadose zone above the cave. High-resolution cave environmental measurements were used to model calcite deposition on one stalagmite in Crag Cave, SW Ireland, and modelled growth over the study period (222 μm over 171 days) is extremely similar to the amount of actual calcite growth (240 μm) over the same time interval, strongly suggesting that equations used to estimate stalagmite growth rates are valid. Although cave air P_CO2 appears to control drip water hydrochemistry in the winter, drip water dilution caused by rain events may have played a larger role during the summer, as evidenced by a series of sudden drops in Ca²⁺ concentrations (dilution) followed by much more gradual increases in drip water Ca²⁺ concentrations (slow addition of diffuse water). This research demonstrates that PCP on stalactites, cave ceilings, and void spaces within the karst above the cave partially controls drip water chemistry, and that thorough characterisation of this process at individual caves is necessary to most accurately interpret climate records from those sites.     

Speleothem calcite farmed in situ: Modern calibration of δO and δC paleoclimate proxies in a continuously-monitored natural cave system

Understanding the relationships between speleothem stable isotopes (δ¹³C δ¹⁸O) and in situ cave forcing mechanisms is important to interpreting ancient stalagmite paleoclimate records. Cave studies have demonstrated that the δ¹⁸O of inorganically precipitated (low temperature) speleothem calcite is systematically heavier than the δ¹⁸O of laboratory-grown calcite for a given temperature. To understand this apparent offset, rainwater, cave drip water, groundwater, and modern naturally precipitated calcite (farmed in situ) were grown at multiple locations inside Hollow Ridge Cave in Marianna, Florida. High resolution micrometeorological, air chemistry time series and ventilation regimes were also monitored continuously at two locations inside the cave, supplemented with periodic bi-monthly air gas grab sample transects throughout the cave.Cave air chemistry and isotope monitoring reveal density-driven airflow pathways through Hollow Ridge Cave at velocities of up to 1.2 m s⁻¹ in winter and 0.4 m s⁻¹ in summer. Hollow Ridge Cave displays a strong ventilation gradient in the front of the cave near the entrances, resulting in cave air that is a mixture of soil gas and atmospheric CO₂. A clear relationship is found between calcite δ¹³C and cave air ventilation rates estimated by proxies pCO₂ and ²²²Rn. Calcite δ¹³C decreased linearly with distance from the front entrance to the interior of the cave during all seasons, with a maximum entrance-to-interior gradient of Δδ¹³C_CaCO3 = −7‰. A whole-cave “Hendy test” at multiple contemporaneous farming sites reveals that ventilation induces a +1.9 ± 0.96‰ δ¹³C offset between calcite precipitated in a ventilation flow path and calcite precipitated on the edge or out of flow paths. This interpretation of the “Hendy test” has implications for interpreting δ¹³C records in ancient speleothems. Calcite δ¹³C_CaCO3 may be a proxy not only for atmospheric CO₂ or overlying vegetation shifts but also for changes in cave ventilation due to dissolution fissures and ceiling collapse creating and plugging ventilation windows.Farmed calcite δ¹⁸O was found to exhibit a +0.82 ± 0.24‰ offset from values predicted by both theoretical calculations and laboratory-grown inorganic calcite. Unlike δ¹³C_CaCO3, oxygen isotopes showed no ventilation effects, i.e. Δδ¹⁸O_CaCO3 appears to be a function of growth temperature only although we cannot rule out a small effect of (unmeasured) gradients in relative humidity (evaporation) accompanying ventilation. Our results support the findings of other cave investigators that water-calcite fractionation factors observed in speleothem calcite are higher that those measured in laboratory experiments. Cave and laboratory calcite precipitates may differ mainly in the complex effects of kinetic isotope fractionation. Combining our data with other recent speleothem studies, we find a new empirical relationship for cave-specific water-calcite oxygen isotope fractionation across a range of temperatures and cave environments:

1000lnα=16.1(10³T^-1)-24.6     

Processes controlling dripwater hydrochemistry variations in Xueyu Cave,SW China: implications for speleothem palaeoclimate signal interpretations

Cave dripwater hydrochemistry responds to environmental changes, both within and outside of the cave, and thereby conveys this information to any stalagmite fed by the drips. As stalagmites are important archives of climate proxy information, understanding how dripwater hydrochemistry responds to environmental forcing is critical. However, despite the large number of speleothems in SW China, the response of dripwater to regional climate variability is not yet adequately understood. A 3‐year study of three drip sites in Xueyu Cave, Chongqing Municipality, SW China, revealed the most important mechanisms controlling dripwater chemical variability. The principal chemical indices (pH, specific conductivity, Ca²⁺, Mg²⁺, Sr²⁺ and ) in collected dripwaters and the local climate data were analysed in this study. The principal controls on the hydrochemistry were found to be the external climate and its changes, groundwater residence time, cave ventilation and prior calcite precipitation (PCP) processes. Dripwater hydrochemistry showed strongly coherent seasonal patterns despite the fact that all sites are Ca–HCO₃ type waters and supersaturated with calcite. Seasonal changes in dripwater hydrochemistry were influenced by the soil and vadose zone CO₂ content as well as groundwater residence time in the upper karst zone. Cave‐air CO₂ seasonal variations were consistent with changes in dripwater PCO₂ and cave ventilation. Trace element ratios (Mg/Ca and Sr/Ca) of dripwater were controlled by PCP processes. Seasonal variations in dripwater Mg/Ca and Sr/Ca ratios in Xueyu Cave showed inverse changes with the Asia Monsoon Index during the monitoring period, reflecting the seasonal climate changes that may have been recorded in the speleothems. Based on a linear regression of PCO₂ and the Ca²⁺ data in the cold–dry winter season, a 130‐ppm shift in cave‐air PCO₂ results in a 1‐ppm shift in dripwater Ca²⁺ concentration in Xueyu Cave. This study illustrates the importance of understanding factors controlling the changes in the composition of dripwater before it reaches the speleothem.     

The Role of CaCO<Subscript>3</Subscript> Reactions in the Contemporary Oceanic CO<Subscript>2</Subscript> Cycle

The present analysis adjusts previous estimates of global ocean CaCO₃ production rates substantially upward, to 133 × 10¹² mol yr^?1 plankton production and 42 × 10¹² mol yr^?1 shelf benthos production. The plankton adjustment is consistent with recent satellite-based estimates; the benthos adjustment includes primarily an upward adjustment of CaCO₃ production on so-called carbonate-poor sedimentary shelves and secondarily pays greater attention to high CaCO₃ mass (calcimass) and turnover of shelf communities on temperate and polar shelves. Estimated CaCO₃ sediment accumulation rates remain about the same as they have been for some years: ~20 × 10¹² mol yr^?1 on shelves and 11 × 10¹² mol yr^?1 in the deep ocean. The differences between production and accumulation of calcareous materials call for dissolution of ~22 × 10¹² mol yr^?1 (~50 %) of shelf benthonic carbonate production and 122 × 10¹² mol yr^?1 (>90 %) of planktonic production. Most CaCO₃ production, whether planktonic or benthonic, is assumed to take place in water depths of <100 m, while most dissolution is assumed to occur below this depth. The molar ratio of CO₂ release to CaCO₃ precipitation (CO₂↑/CaCO₃↓) is <1.0 and varies with depth. This ratio, Ψ, is presently about 0.66 in surface seawater and 0.85 in ocean waters deeper than about 1000 m. The net flux of CO₂ associated with CaCO₃ reactions in the global ocean in late preindustrial time is estimated to be an apparent influx from the atmosphere to the ocean, of +7 × 10¹² mol C yr^?1, at a time scale of 10²–10³ years. The CaCO₃-mediated influx of CO₂ is approximately offset by CO₂ release from organic C oxidation in the water column. Continuing ocean acidification will have effects on CaCO₃ and organic C metabolic responses to the oceanic inorganic C cycle, although those responses remain poorly quantified.     

Human impact in a tourist karstic cave (Aracena, Spain)

 Human intrusion on the Cave of Marvels (southwestern Spain) has produced a series of effects on the water (fall in the level of the pools due to pumping from nearby wells), the air (increased temperature and CO₂ concentration as well as decreased relative humidity) and the rock. In addition, plant colonization, favored by the lighting system, has irreversibly altered numerous speleothems. The processes of degradation are especially intense in the sectors with less air volume and limited ventilation. The analysis of the cave deposits by scanning electron microscopy and thin section analysis revealed that floral pollution constitutes one of the most aggressive agents against the calcite and aragonite precipitates, being responsible for biochemical and biophysical degradation of the first order. Received: 18 December 1995 · Accepted 10 September 1996     

Effect of injected CO2 on geochemical alteration of the Altmark gas reservoir in Germany

Capturing CO₂ from point sources and storing it in geologic formations is a potential option for allaying the CO₂ level in the atmosphere. In order to evaluate the effect of geological storage of CO₂ on rock-water interaction, batch experiments were performed on sandstone samples taken from the Altmark reservoir, Germany, under in situ conditions of 125 °C and 50 bar CO₂ partial pressure. Two sets of experiments were performed on pulverized sample material placed inside a closed batch reactor in (a) CO₂ saturated and (b) CO₂ free environment for 5, 9 and 14 days. A 3M NaCl brine was used in both cases to mimic the reservoir formation water. For the “CO₂ free” environment, Ar was used as a pressure medium. The sandstone was mainly composed of quartz, feldspars, anhydrite, calcite, illite and chlorite minerals. Chemical analyses of the liquid phase suggested dissolution of both calcite and anhydrite in both cases. However, dissolution of calcite was more pronounced in the presence of CO₂. In addition, the presence of CO₂ enhanced dissolution of feldspar minerals. Solid phase analysis by X-ray diffraction and Mössbauer spectroscopy did not show any secondary mineral precipitation. Moreover, Mössbauer analysis did not show any evidence of significant changes in redox conditions. Calculations of total dissolved solids’ concentrations indicated that the extent of mineral dissolution was enhanced by a factor of approximately 1.5 during the injection of CO₂, which might improve the injectivity and storage capacity of the targeted reservoir. The experimental data provide a basis for numerical simulations to evaluate the effect of injected CO₂ on long-term geochemical alteration at reservoir scale.     

Microclimatic characterization of a karstic cave: human impact on microenvironmental parameters of a prehistoric rock art cave (Candamo Cave, northern Spain)

The Candamo Cave contains an important group of paleolithic paintings which have been seriously deteriorated due to mass tourism. In this work, an analysis was carried out of different climatic parameters (CO₂, temperature, humidity, ²²²Rn) during annual cycles with the cave closed to the public and during an experimental period of controlled visits. The effect of visits on the geochemical characteristics of karstic water was also analyzed together with the cave ventilation. The natural variations in the cave air CO₂ were above 3000 ppm, the increase produced through visits was only 100–110 ppm and since the humidity is almost permanently at saturation point, the critical parameter which limits the visitor capacity becomes air temperature. The temperature changes during the annual cycle are of the order of 1  °C in the external part and less than 0.5  °C in the internal part of the cave and a maximum increase of 0.13  °C was observed during the period of the visits. The ²²²Rn and CO₂ concentration minimums in the summer period (July–October) show that this is the most propitious time for visits, since the greatest ventilation is produced in the cave at this time and, therefore, the greatest capacity for recovery. The geochemistry of the water, on the other hand, indicated that this is the period of the year in which processes of wall corrosion can be most easily introduced, although this would be of limited magnitude. The visitor capacity calculated was 29 visitors/day. Received: 29 August 1996 · Accepted: 23 June 1997     

CO<Subscript>2</Subscript> flux variation and its contribution area in the debris-covered area of Koxkar Glacier,Mt. Tianshan in China

During the formation and development of glacial meltwater runoff, hydrochemical erosion is abundant, especially the hydrolysis of K/Na feldspar and carbonates, which can consume H⁺ in the water, promote the formation of bicarbonate by dissolving atmospheric CO₂, and affect the regional carbon cycle. From July 21, 2015, to July 18, 2017, the CO₂ concentration and flux were observed by the eddy covariance (EC) method in the relatively flat and open moraine cover area of Koxkar Glacier in western Mt. Tianshan, China. We found that: (1) atmospheric CO₂ fluxes ranged from ??408.95 to 81.58 mmol m^?2 day^?1 (average ? 58.68 mmol m^?2 day^?1), suggesting that the study area is a significant carbon sink, (2) the CO₂ flux footprint contribution areas were primarily within 150 m of the EC station, averaging total contribution rates of 93.30%, 91.39%, and 90.17% of the CO₂ flux in the snow accumulation, snow melting, and glacial melting periods, respectively. Therefore, the contribution areas with significant influences on CO₂ flux observed at EC stations were concentrated, demonstrating that grassland CO₂ flux around the glaciers had little effect at the EC stations, (3) in the predominant wind direction, under stable daytime atmospheric stratification, the measurement of CO₂ flux, as interpreted by the Agroscope Reckenholz Tanikon footprint tool, was 79.09% ± 1.84% in the contribution area. This was slightly more than seen at night, but significantly lower than the average under unstable atmospheric stratification across the three periods of interest (89%). The average distance of the farthest point of the flux footprint under steady state atmospheric conditions was 202.61?±?69.33 m, markedly greater than that under non-steady state conditions (68.55?±?10.34 m). This also indicates that the CO₂ flux observed using EC was affected primarily by hydrochemical erosion reactions in the glacier area, (4) a good negative correlation was found between net glacier exchange (NGE) of CO₂ and air temperature on precipitation-free days. Strong ice and snow ablation could promote hydrochemical reactions of soluble substances in the debris area and accelerated sinking of atmospheric CO₂. Precipitation events might reduce snow and ice melting, driven by reduced regional temperatures. However, a connection between NGE and precipitation, when less than 8.8 mm per day, was not obvious. When precipitation was greater than 8.8 mm per day, NGE decreased with increasing precipitation, (5) graphically, the slope of NGE, related to daily runoff, followed a trend: snow melting period?>?snow accumulation period?>?early glacial ablation period?>?late glacier ablation period?>?dramatic glacier ablation period. The slope was relatively large during snow melting, likely because of CO₂ sinking caused by water–rock interactions. The chemical reaction during elution in the snow layer might also promote atmospheric CO₂ drawdown. At the same time, the damping effect of snow cover and the almost-closed glacier hydrographic channel inhibited the formation of regional runoff, possibly providing sufficient time for the chemical reaction, thus promoting further CO₂ drawdown.     

Experimental and numerical simulation study of the mineral sequestration mechanism of the Shiqianfeng saline aquifers in the Ordos Basin,Northwest China

This study focused on typical injection layers of deep saline aquifers in the Shiqianfeng Formation used in the Carbon Capture and Sequestration Demonstration Projects in the Ordos Basin, Northwest China. The study employed experiments and numerical simulations to investigate the mechanism of CO₂ mineral sequestration in these deep saline aquifers. The experimental results showed that the dissolved minerals are plagioclase, hematite, illite–smectite mixed layer clay and illite, whereas the precipitated minerals are quartz (at 55, and 70 °C) and kaolinite (at 70 °C). There are rare carbonate mineral precipitations at the experimental time scale, while the precipitation of quartz as a product of the dissolution of silicate minerals and some intermediate minerals rich in K and Mg that transform to clay minerals, reveals the possibility of carbonate precipitation at the longer time scale. These results are consistent with some results previously reported in the literature. We calibrated the kinetic parameters of mineral dissolution and precipitation by these experimental results and then simulated the CO₂ mineral sequestration under deep saline aquifer conditions. The simulation results showed that the dissolved minerals are albite, anorthite and minor hematite, whereas the precipitated minerals are calcite, kaolinite and smectite at 55 and 70 °C. The geochemical reaction of illite is more complex. At 55 °C, illite is dissolved at the relatively lag time and transformed to dawsonite; at 70 °C, illite is precipitated in the early reaction period and then transformed to kaolinite. Based on this research, sequestrated CO₂ minerals, which are mainly related to the temperature of deep saline aquifers in Shiqianfeng Fm., are calcite and dawsonite at lower temperature, and calcite at higher temperature. The simulation results also establish that calcite could precipitate over a time scale of thousands of years, and the higher the temperature the sooner such a process would occur due to increased reaction rates. These characteristics are conducive, not only to the earlier occurrence of mineral sequestration, but also increase the sequestration capacity of the same mineral components. For a sequestration period of 10,000 years, we determined that the mineral sequestration capacity is 0.786 kg/m³ at 55 °C, and 2.180 kg/m³ at 70 °C. Furthermore, the occurrence of mineral sequestration indirectly increases the solubility of CO₂ in the early reaction period, but this decreases with the increase in temperature.     

Carbon mass-balance modelling and carbon isotope exchange processes in dynamic caves

Diverse interpretations have been made of carbon isotope time series in speleothems, reflecting multiple potential controls. Here we study the dynamics of ¹³C and ¹²C cycling in a particularly well-constrained site to improve our understanding of processes affecting speleothem δ¹³C values. The small, tubular Grotta di Ernesto cave (NE Italy) hosts annually-laminated speleothem archives of climatic and environmental changes. Temperature, air pressure, pCO₂, dissolved inorganic carbon (DIC) and their C isotopic compositions were monitored for up to five years in soil water and gas, cave dripwater and cave air. Mass-balance models were constructed for CO₂ concentrations and tested against the carbon isotope data. Air advection forces winter pCO₂ to drop in the cave air to ca. 500 ppm from a summer peak of ca. 1500 ppm, with a rate of air exchange between cave and free atmosphere of approximately 0.4 days. The process of cave ventilation forces degassing of CO₂ from the dripwater, prior to any calcite precipitation onto the stalagmites. This phase of degassing causes kinetic isotope fractionation, i.e. ¹³C-enrichment of dripwater whose δ¹³C_DIC values are already higher (by about 1‰) than those of soil water due to dissolution of the carbonate rock. A subsequent systematic shift to even higher δ¹³C values, from −11.5‰ in the cave drips to about −8‰ calculated for the solution film on top of stalagmites, is related to degassing on the stalagmite top and equilibration with the cave air. Mass-balance modelling of C fluxes reveals that a very small percentage of isotopically depleted cave air CO₂ evolves from the first phase of dripwater degassing, and shifts the winter cave air composition toward slightly more depleted values than those calculated for equilibrium. The systematic ¹³C-enrichment from the soil to the stalagmites at Grotta di Ernesto is independent of drip rate, and forced by the difference in pCO₂ between cave water and cave air. This implies that speleothem δ¹³C values may not be simply interpreted either in terms of hydrology or soil processes.     

Suitability of Existing Numerical Model Codes and Thermodynamic Databases for the Prognosis of Calcite Dissolution Processes in Near-Surface Sediments Due to a CO2 Leakage Investigated by Column Experiments

Among the risks of CO₂ storage is the potential of CO₂ leakage into overlaying formations and near-surface potable aquifers. Through a leakage, the CO₂ can intrude into protected groundwater resources, which can lead to groundwater acidification followed by potential mobilisation of heavy metals and other trace metals through mineral dissolution or ion exchange processes. The prediction of pH buffer reactions in the formations overlaying a CO₂ storage site is essential for assessing the impact of CO₂ leakages in terms of trace metal mobilisation. For buffering the pH-value, calcite dissolution is one of the most important mechanisms. Although calcite dissolution has been studied for decades, experiments conducted under elevated CO₂ partial pressures are rare. Here, the first study for column experiments is presented applying CO₂ partial pressures from 6 to 43 bars and realising a near-natural flow regime. Geochemical calculations of calcite dissolution kinetics were conducted using PHREEQC together with different thermodynamic databases. Applying calcite surface areas, which were previously acquired by N₂-BET or calculated based on grain diameters, respectively, to the rate laws according to Plummer et al. (Am J Sci 278:179–216, doi:10.2475/ajs.278.2.179, 1978) or Palandri and Kharaka (US Geol Surv Open file Rep 2004–1068:71, 2004) in the numerical simulations led to an overestimation of the calcite dissolution rate by up to three orders of magnitude compared to the results of the column experiments. Only reduction of the calcite surface area in the simulations as a fitting procedure allowed reproducing the experimental results. A reason may be that the diffusion boundary layer (DBL), which depends on the groundwater flow velocity and develops at the calcite grain surface separating it from the bulk of the solution, has to be regarded: The DBL leads to a decrease in the calcite dissolution rate under natural laminar flow conditions compared to turbulent mixing in traditional batch experiments. However, varying the rate constants by three orders of magnitudes in a field scale PHREEQC model simulating a CO₂ leakage produced minor variations in the pH buffering through calcite dissolution. This justifies the use of equilibrium models when calculating the calcite dissolution in CO₂ leakage scenarios for porous aquifers and slow or moderate groundwater flow velocities. However, the selection of the thermodynamic database has an impact on the dissolved calcium concentration, leading to an uncertainty in the simulation results. The resulting uncertainty, which applies also to the calculated propagation of an aquifer zone depleted in calcite through dissolution, seems negligible for shallow aquifers of approximately 60 m depth, but amounts to 35 % of the calcium concentration for aquifers at a depth of approximately 400 m.     

Origin and palaeoenvironmental significance of lamination in stalagmites from Katerloch Cave,Austria

The origin and environmental dependencies of lamination in stalagmites from Katerloch, common in speleothems from other cave sites, are examined in detail. Petrographic observations and chemical analyses (including isotopes) of stalagmites and modern calcite were combined with multi‐annual cave monitoring. All investigated stalagmites are composed of low‐Mg calcite and show white, porous laminae and typically thinner, translucent dense laminae. The binary lamination pattern results from changes in the calcite fabric: white, porous laminae are characterized by a high porosity and abundant fluid inclusions and also by enhanced vertical growth and thinning towards the flanks. Translucent, dense laminae exhibit a compact fabric and constant thickness of individual growth layers. U‐Th dating supports an annual origin of the lamination and the seasonally changing intensity of cave ventilation provides a robust explanation for the observed relationships between lamination, stable C isotopic compositions and trace elements (Mg, Sr and Ba). The seasonally variable air exchange, driven by temperature contrasts between the cave interior and outside atmosphere, modulates the rate and amount of CO₂ degassing from the drip water and affects the hydrochemistry and consequently the fabric of the precipitating calcite. Although cave air composition and drip rate are both major variables in controlling CO₂ degassing from the drip water, the seasonally changing ventilation in Katerloch exerts the primary control and the results suggest a secondary (amplifying/attenuating) influence of the drip rate. Drip rate, however, might be the controlling parameter for lamina development at cave sites experiencing only small seasonal cave air exchange. Importantly, the seasonally variable composition of drip water does not reflect the seasonal cycle of processes in the soil zone, but results from exchange with the cave atmosphere. The alternating porous and dense calcite fabric is the expression of a variable degree of lateral coalescence of smaller crystallites forming large columnar crystals. The white, porous laminae represent partial coalescence and form during the warm season: low calcite δ¹³C values are linked to low δ¹³C values of cave air and drip water during that time. This observation corresponds to times of reduced cave ventilation, high pCO₂ of cave air, low drip water pH, lower calcite supersaturation and typically high drip rates. In contrast, the translucent, dense laminae represent more or less complete lateral coalescence (inclusion‐free) during the cold season (high calcite, drip water and cave air δ¹³C values), i.e. times of enhanced cave ventilation, low cave air pCO₂, increased drip water pH, relatively high calcite supersaturation and typically low drip rates. In essence, the relative development of the two lamina types reflects changes in the seasonality of external air temperature and precipitation, with a strong control of the winter air temperature on the intensity of cave‐air exchange. Thick translucent, dense laminae are favoured by long, cold and wet winters and such conditions may be related closely to the North Atlantic Oscillation mode (weak westerlies) and enhanced Mediterranean cyclone activity during the cold season. Studies of speleothem lamination can thus help to better understand (and quantify) the role of seasonality changes, for example, during rapid climate events.     

Temperature-driven meltwater production and hydrochemical variations at a glaciated alpine karst aquifer: implication for the atmospheric CO<Subscript>2</Subscript> sink under global warming

About two hydrological years of continuous data of discharge, temperature, electrical conductivity and pH have been recorded at the Glarey spring in the Tsanfleuron glaciated karst area in the Swiss Alps, to understand how glaciated karst aquifer systems respond hydrochemically to diurnal and seasonal recharge variations, and how calcite dissolution by glacial meltwater contributes to the atmospheric CO₂ sink. A thermodynamic model was used to link the continuous data to monthly water quality data allowing the calculation of CO₂ partial pressures and calcite saturation indexes. The results show diurnal and seasonal hydrochemical variations controlled chiefly by air temperature, the latter influencing karst aquifer recharge by ice and snowmelt. Karst process-related atmospheric CO₂ sinks were more than four times higher in the melting season than those in the freezing season. This finding has implication for understanding the atmospheric CO₂ sink in glaciated carbonate rock terrains: the carbon sink will increase with increasing runoff caused by global warming, i.e., carbonate weathering provides a negative feedback for anthropogenic CO₂ release. However, this is a transient regulation effect that is most efficient when glacial meltwater production is highest, which in turn depends on the future climatic evolution.     

A multiple cave deposit assessment of suitability of speleothem isotopes for reconstructing palaeo‐vegetation and palaeo‐temperature

The suitability of speleothems for interpreting palaeoclimate is typically determined by using either the Hendy Test, overlapping analysis or long‐term cave environment monitoring. However, in many cases, these methods are not applicable, because a speleothem lacks clearly traceable layers for the Hendy Test, it is difficult to obtain an overlapping speleothem nearby, or long‐term cave monitoring is impractical. The authors propose a multiple cave deposit approach to assess the suitability of speleothems for palaeoclimate study. Speleothems collected from two sites within Raccoon Mountain Cave, Tennessee (USA) exhibit remarkable spatial variation (δ¹³C: ?10·3‰ to ?2·2‰) over a relatively short distance (ca 260 m). Drip water δ¹⁸O values exhibit a seasonal precipitation signal at Site 1 and an annual signal at Site 2. Combining field observations, water isotope analysis and trace‐element data, the authors propose that the speleothem formation at Site 1 and Site 2 tapped distinct sources of CO₂: (i) CO₂ derived from overlying soils for Site 1; and (ii) limestone dissolved inorganic carbon induced by ground water dissolution for Site 2. Using fresh cave deposits (modern speleothem) δ¹³C (100% C3 vegetation) as an analogue, a simple model was developed to estimate land surface vegetation for speleothems. The speleothem formation temperature estimated using fresh cave deposit δ¹⁸O values generally reflects the mean annual temperature in this region. This study indicates that spatial variations in carbon isotopes could be caused by different carbon sources dominating in different parts of the cave, which should be taken into consideration by researchers when using speleothem δ¹³C values to reconstruct temporal palaeo‐vegetation changes. This study demonstrates a practical sampling strategy for verifying suitability of speleothems for palaeo‐vegetation and palaeo‐temperature reconstructions by analysing multiple cave deposits, especially for cases in which the Hendy Test, parallel sampling and long‐term monitoring of cave environment are not feasible.     

Dissolution sequestration mechanism of CO<Subscript>2</Subscript> at the Shiqianfeng saline aquifer in the Ordos Basin,northwestern China

This study focused on the target injection layers of deep saline aquifers in the Shiqianfeng Fm. in the Carbon Capture and Sequestration (CCS) Demonstration Projects in the Ordos Basin, northwestern China. The study employed a combination method of experiments and numerical simulation to investigate the dissolution mechanism and impact factors of CO₂ in these saline aquifers. The results showed (1) CO₂ solubility in different types of water chemistry were shown in ascending order: MgCl₂-type water < CaCl₂-type water < Na₂SO₄-type water < NaCl-type water < Na₂CO₃-type water < distilled water. These results were consistent with the calculated results undertaken by TOUGHREACT with about 5% margin of error. CO₂ solubility of Shiqianfeng Fm. saline was 1.05 mol/L; (2) compared with distilled water, the more complex the water’s chemical composition, the greater the increase in HCO₃ ^?concentration. While the water’s composition was relatively simple, the tested water’s HCO₃ ^?concentrations were in close accord with the calculated value undertaken by the TOUGHREACT code, and the more complex the water’s composition, the poorer the agreement was, probably due to the complex and unstable HCO₃ ^? complicating matters when in an aqueous solution system including both tested HCO₃ ^?concentration and calculated HCO₃ ^?concentration; (3) the CO₂ solubility in the saline at the temperature conditions of 55 °C and 70 °C were 1.17 and 1.02 mol/L. When compared with the calculated value of 1.20 and 1.05 mol/L, they were almost the same with only 1 and 3% margin of error; concentrations of HCO₃ ^? were 402.73 mg/L (0.007 mol/L) and 385.65 mg/L (0.006 mol/L), while the simulation results were 132.16 mg/L (0.002 mol/L) and 128.52 mg/L (0.002 mol/L). From the contrast between the tested data and the calculated data undertaken by the TOUGHREACT code, it was shown that TOUGHRACT code could better simulate the interaction between saline and CO₂ in the dissolution sequestration capacity. Therefore, TOUGHREACT code could be used for the inter-process prediction of CO₂ long-term geological storage of CO₂; (4) The Ca²⁺ concentration and SO₄ ^2?concentration in saline water had less effect on the solubility of CO₂ and HCO₃ ^?concentration. In addition, TDS and pH values of saline affected not only the solubility of CO₂, but also the conversion of CO₂ to HCO₃ ^? due to that they can affect the activity and acid-base balance. So in fact, we just need to consider that the TDS and pH values are main impact factors in the dissolution sequestration capacity of CO₂ geological sequestration in deep saline aquifers.     

The effect of alkalis and polymerization on the solubility of H<Subscript>2</Subscript>O and CO<Subscript>2</Subscript> in alkali-rich silicate melts

The effect of alkalis on the solubility of H₂O and CO₂ in alkali-rich silicate melts was investigated at 500 MPa and 1,250 °C in the systems with H₂O/(H₂O + CO₂) ratio varying from 0 to 1. Using a synthetic analog of phonotephritic magma from Alban Hills (AH1) as a base composition, the Na/(Na + K) ratio was varied from 0.28 (AH1) to 0.60 (AH2) and 0.85 (AH3) at roughly constant total alkali content. The obtained results were compared with the data for shoshonitic and latitic melts having similar total alkali content but different structural characteristics, e.g., NBO/T parameter (the ratio of non-bridging oxygens over tetrahedrally coordinated cations), as those of the AH compositions. Little variation was observed in H₂O solubility (melt equilibrated with pure H₂O fluid) for the whole compositional range in this study with values ranging between 9.7 and 10.2 wt. As previously shown, the maximum CO₂ content in melts equilibrated with CO₂-rich fluids increases strongly with the NBO/T from 0.29 wt % for latite (NBO/T = 0.17) to 0.45 wt % for shoshonite (NBO/T = 0.38) to 0.90 wt % for AH2 (NBO/T = 0.55). The highest CO₂ contents determined for AH3 and AH1 are 1.18 ± 0.05 wt % and 0.86 ± 0.12 wt %, respectively, indicating that Na is promoting carbonate incorporation stronger than potassium. At near constant NBO/T, CO₂ solubility increases from 0.86 ± 0.12 wt % in AH1 [Na/(Na + K)] = 0.28, to 1.18 ± 0.05 wt % in AH3 [Na/(Na + K)] = 0.85, suggesting that Na favors CO₂ solubility on an equimolar basis. An empirical equation is proposed to predict the maximum CO₂ solubility at 500 MPa and 1,100–1,300 °C in various silicate melts as a function of the NBO/T, (Na + K)/∑cations and Na/(Na + K) parameters: \({\text{wt}}\% \;{\text{CO}}_{2} = - 0.246 + 0.014\exp \left( {6.995 \cdot \frac{\text{NBO}}{T}} \right) + 3.150 \cdot \frac{{{\text{Na}} + {\text{K}}}}{{\varSigma {\text{cations}}}} + 0.222 \cdot \frac{\text{Na}}{{{\text{Na}} + {\text{K}}}}.\) This model is valid for melt compositions with NBO/T between 0.0 and 0.6, (Na + K)/∑cation between 0.08 and 0.36 and Na/(Na + K) ratio from 0.25 to 0.95 at oxygen fugacities around the quartz–fayalite–magnetite buffer and above.