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.     

Climatic dependence of stable carbon and oxygen isotope signals recorded in speleothems: From soil water to speleothem calcite

Understanding the relationship between stable isotope signals recorded in speleothems (δ¹³C and δ¹⁸O) and the isotopic composition of the carbonate species in the soil water is of great importance for their interpretation in terms of past climate variability. Here the evolution of the carbon isotope composition of soil water on its way down to the cave during dissolution of limestone is studied for both closed and open-closed conditions with respect to CO₂.The water entering the cave flows as a thin film towards the drip site. CO₂ degasses from this film within approx. 10 s by molecular diffusion. Subsequently, chemical and isotopic equilibrium is established on a time scale of several 10-100 s. The δ¹³C value of the drip water is mainly determined by the isotopic composition of soil CO₂. The evolution of the δ¹⁸O value of the carbonate species is determined by the long exchange time T_ex, between oxygen in carbonate and water of several 10,000 s. Even if the oxygen of the CO₂ in soil water is in isotopic equilibrium with that of the water, dissolution of limestone delivers oxygen with a different isotopic composition changing the δ¹⁸O value of the carbonate species. Consequently, the δ¹⁸O value of the rainwater will only be reflected in the drip water if it has stayed in the rock for a sufficiently long time.After the water has entered the cave, the carbon and oxygen isotope composition of the drip water may be altered by CO₂-exchange with the cave air. Exchange times, , of about 3000 s are derived. Thus, only drip water, which drips in less than 3000 s onto the stalagmite surface, is suitable to imprint climatic signals into speleothem calcite deposited from it.Precipitation of calcite proceeds with time constants, τ_p, of several 100 s. Different rate constants and equilibrium concentrations for the heavy and light isotopes, respectively, result in isotope fractionation during calcite precipitation. Since T_ex ? τ_p, exchange with the oxygen in the water can be neglected, and the isotopic evolution of carbon and oxygen proceed analogously. For drip intervals T_d < 0.1τ_p the isotopic compositions of both carbon and oxygen in the solution evolve linearly in time. The calcite precipitated at the apex of the stalagmite reflects the isotopic signal of the drip water.For long drip intervals, when calcite is deposited from a stagnant water film, long drip intervals may have a significant effect on the isotopic composition of the DIC. In this case, the isotopic composition of the calcite deposited at the apex must be determined by averaging over the drip interval. Such processes must be considered when speleothems are used as proxies of past climate variability.     

δC profiles along growth layers of stalagmites: Comparing theoretical and experimental results

The isotopic carbon ratio of a calcite-precipitating solution flowing as a water film on the surface of a stalagmite is determined by Rayleigh distillation. It can be calculated, when the -concentration of the solution at each surface point of the stalagmite and the fractionation factors are known. A stalagmite growth model based entirely on the physics of laminar flow and the well-known precipitation rates of a supersaturated solution of calcite, without any further assumptions, is employed to obtain the spatial distribution of the -concentration, which contributes more than 95% to the dissolved inorganic carbon (DIC). The δ¹³C profiles are calculated along the growth surface of a stalagmite for three cases: (A) isotopic equilibrium of both CO₂ outgassing and calcite precipitation; (B) outgassing of CO₂ is irreversible but calcite precipitation is in isotopic equilibrium. (C) Both CO₂ outgassing and calcite precipitation are irreversible. In all cases the isotopic shift δ¹³C increases from the apex along the distance on a growth surface. In cases A and B, calcite deposited at the apex is in isotopic equilibrium with the solution of the drip water. The difference between δ¹³C at the apex and the end of the growth layer is independent of the stalagmite’s radius, but depends on temperature. For case A, it is about half the value obtained for cases B and C. In case C, the isotopic composition of calcite at the apex equals that of the drip water, but further out it becomes practically identical with that of case B. The growth model has been applied to field data of stalagmite growth, where the thickness and the δ¹³C of calcite precipitated to a glass plate located on the top of a stalagmite have been measured as function of the distance from the drip point. The calculated data are in good agreement to the observed ones and indicate that deposition occurred most likely under conditions B, eventually also C. A sensitivity analysis has been performed, which shows that within the limits of observed external parameters, such as drip rates and partial pressure of carbon dioxide P_CO2 in the cave, the results remain valid.     

Controls on ostracod valve geochemistry: Part 2. Carbon and oxygen isotope compositions

The stable carbon and oxygen isotope compositions of fossil ostracods are powerful tools to estimate past environmental and climatic conditions. The basis for such interpretations is that the calcite of the valves reflects the isotopic composition of water and its temperature of formation. However, calcite of ostracods is known not to form in isotopic equilibrium with water and different species may have different offsets from inorganic precipitates of calcite formed under the same conditions. To estimate the fractionation during ostracod valve calcification, the oxygen and carbon isotope compositions of 15 species living in Lake Geneva were related to their autoecology and the environmental parameters measured during their growth. The results indicate that: (1) Oxygen isotope fractionation is similar for all species of Candoninae with an enrichment in ¹⁸O of more than 3‰ relative to equilibrium values for inorganic calcite. Oxygen isotope fractionation for Cytheroidea is less discriminative relative to the heavy oxygen, with enrichments in ¹⁸O for these species of 1.7 to 2.3‰. Oxygen isotope fractionations for Cyprididae are in-between those of Candoninae and Cytheroidea. The difference in oxygen isotope fractionation between ostracods and inorganic calcite has been interpreted as resulting from a vital effect. (2) Comparison with previous work suggests that oxygen isotope fractionation may depend on the total and relative ion content of water. (3) Carbon isotope compositions of ostracod valves are generally in equilibrium with DIC. The specimens’ δ¹³C values are mainly controlled by seasonal variations in δ¹³C_DIC of bottom water or variation thereof in sediment pore water. (4) Incomplete valve calcification has an effect on carbon and oxygen isotope compositions of ostracod valves. Preferential incorporation of at the beginning of valve calcification may explain this effect. (5) Results presented here as well as results from synthetic carbonate growth indicate that different growth rates or low pH within the calcification site cannot be the cause of oxygen isotope ‘vital effects’ in ostracods. Two mechanisms that might enrich the ¹⁸O of ostracod valves are deprotonation of that may also contribute to valve calcification, and effects comparable to salt effects with high concentrations of Ca and/or Mg within the calcification site that may also cause a higher temperature dependency of oxygen isotope fractionation.     

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     

Modelling δC and δO in the solution layer on stalagmite surfaces

We derive equations describing the evolution of the carbon and oxygen isotope composition of the bicarbonate in a calcite precipitating solution on the surface of a stalagmite using a classical Rayleigh approach. The combined effects of calcite precipitation, degassing of CO₂ and the buffering effect of the water reservoir are taken into account. Whereas δ¹³C shows a progressive increase to a final constant value, δ¹⁸O shows an initial isotopic enrichment, which exponentially decays due to the buffering effect of the water reservoir. The calculated evolution is significantly different for both carbon and oxygen isotopes than derived in a recent paper [Dreybrodt W. (2008) Evolution of the isotopic composition of carbon and oxygen in a calcite precipitating H₂O-CO₂-CaCO₃ solution and the related isotopic composition of calcite in stalagmites. Geochim. Cosmochim. Acta72, 4712-4724.].Furthermore, we discuss the isotopic evolution of the bicarbonate in the solution for long residence times on the stalagmite surface, i.e., for t→∞. The equilibrium isotope ratio of the bicarbonate is then determined by isotopic exchange between the cave atmosphere and the bicarbonate in the solution and can be calculated by equilibrium isotope fractionation. For strongly ventilated caves exchange with the cave atmosphere will result in higher δ¹³C and δ¹⁸O values than those observed in a pure Rayleigh distillation scenario, for sparsely ventilated caves it will result in lower δ¹³C and δ¹⁸O values.     

Modeling speleothem δC variability in a central Sierra Nevada cave using C and Sr/Sr

Carbon isotopes in speleothems can vary in response to a number of complex processes active in cave systems that are both directly and indirectly related to climate. Progressing downward from the soil zone overlying the cave, these processes include soil respiration, fluid-rock interaction in the host limestone, degassing of CO₂ and precipitation of calcite upflow from the speleothem drip site, and calcite precipitation at the drip site. Here we develop a new approach to independently constrain the roles of water-rock interaction and soil processes in controlling stalagmite δ¹³C. This approach uses the dead carbon proportion (dcp) estimated from coupled ¹⁴C and ²³⁰Th/U measurements, in conjunction with Sr isotope analyses on stalagmite calcite from a central Sierra Nevada foothills cave in California, a region characterized by a highly seasonal Mediterranean-type climate, to determine the roles of water-rock interaction and soil processes in determining stalagmite δ¹³C. Increases in stalagmite dcp between 16.5 and 8.8 ka are coincident with decreased δ¹³C, indicating a varying yet substantial contribution from the soil organic matter (SOM) reservoir, likely due to significantly increased average age of SOM in the soil veneer above the cave during wet climatic intervals.We use geochemical and isotope mixing models to estimate the host-carbonate contribution throughout the δ¹³C time series and determine the degree of degassing and calcite precipitation that occurred prior to precipitation of stalagmite calcite. The degree of degassing and prior calcite precipitation we calculate varies systematically with other climate indicators, with less degassing and prior calcite precipitation occurring during wetter climatic intervals and more during drier intervals. Modeled δ¹³C values and degassing calculations suggest that some degree of prior calcite precipitation is necessary at all time intervals to explain measured stalagmite δ¹³C values, even during relatively wet intervals. These results illustrate the importance of constraining degassing and prior calcite precipitation in the interpretation of speleothem δ¹³C records, particularly those from caves that formed in seasonal semi-arid to arid environments.     

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.     

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.     

Stable Carbon Isotope Variations in Cave Percolation Waters and their Implications in Four Caves of Guizhou, China

Monitoring and sampling of main plants,soil CO2,soil water,bedrock,spring water,drip water and its corresponding speleothem were performed at four cave systems of Guizhou,Southwest China,from April 2003 to May 2004,in order to understand stable carbon isotope ratios variations of dissolved inorganic Carbon(DIC) in cave percolation waters(δ13CDIC) and their implications for paleoclimate.Stable carbon isotopic compositions and ions(Ca,Mg,Sr,SO4,Cl etc.) were measured for all samples.The results indicate that there are significant differences among the δ13CDIC values from inter-cave,even inter-drip of intra-cave in the four caves.The δ13CDIC values from the Liangfeng Cave(LFC) is lightest among the four caves,where vegetation type overlying the cave is primary forest dominated by tall trees with lighter average δ13C value(–29.9‰).And there are remarkable differences in δ13CDIC values of different drip waters in the Qixing Cave(QXC) and Jiangjun Cave(JJC),up to 6.9‰ and 7.8‰,respectively.Further analyses show that the δ13CDIC values in cave drip waters are not only controlled by vegetation biomass overlying the cave,but also hydro-geochemical processes.Therefore,accurate interpreting of δ13C recorded in speleothems cannot be guaranteed if these effects of the above mentioned factors are not taken into consideration.     

CO clumping in speleothems: Observations from natural caves and precipitation experiments

The oxygen isotope composition of speleothems is an important proxy of continental paleoenvironments, because of its sensitivity to variations in cave temperature and drip water δ¹⁸O. Interpreting speleothem δ¹⁸O records in terms of absolute paleotemperatures and δ¹⁸O values of paleo-precipitation requires quantitative separation of the effects of these two parameters, and correcting for possible kinetic isotope fractionation associated with precipitation of calcite out of thermodynamic equilibrium. Carbonate clumped-isotope thermometry, based on measurements of Δ₄₇ (a geochemical variable reflecting the statistical overabundance of ¹³C¹⁸O bonds in CO₂ evolved from phosphoric acid digestion of carbonate minerals), potentially provides a method for absolute speleothem paleotemperature reconstructions independent of drip water composition. Application of this new technique to karst records is currently limited by the scarcity of published clumped-isotope studies of modern speleothems. The only modern stalagmite reported so far in the literature yielded a lower Δ₄₇ value than expected for equilibrium precipitation, possibly due to kinetic isotope fractionation.Here we report Δ₄₇ values measured in natural speleothems from various cave settings, in carbonate produced by cave precipitation experiments, and in synthetic stalagmite analogs precipitated in controlled laboratory conditions designed to mimic natural cave processes. All samples yield lower Δ₄₇ and heavier δ¹⁸O values than predicted by experimental calibrations of thermodynamic equilibrium in inorganic calcite. The amplitudes of these isotopic disequilibria vary between samples, but there is clear correlation between the amount of Δ₄₇ disequilibrium and that of δ¹⁸O. Even pool carbonates believed to offer excellent conditions for equilibrium precipitation of calcite display out-of-equilibrium δ¹⁸O and Δ₄₇ values, probably inherited from prior degassing within the cave system.In addition to these modern observations, clumped-isotope analyses of a flowstone from Villars cave (France) offer evidence that the amount of disequilibrium affecting Δ₄₇ in a single speleothem can experience large variations at time scales of 10 kyr. Application of clumped-isotope thermometry to speleothem records calls for an improved physical understanding of DIC fractionation processes in karst waters, and for the resolution of important issues regarding equilibrium calibration of Δ₄₇ in inorganic carbonates.     

Chemical and isotopic characteristics of weathering and nitrogen release in non-glacial drainage waters on Arctic tundra

Soil-water interactions in coastal tundra soils are a potential source of nutrients for surrounding fjordal and coastal ecosystems. Changes in water chemistry and stable isotope composition from three streams in west Spitsbergen were examined to assess the sources and losses of nitrogen, sulfur and carbon in thin organic tundra soils overlying sediments. Studies were undertaken from snowmelt (mid June) through to the end of the summer (September) in both 2001 and 2002. Drainage water chemistry was dominated by the solution of Ca-Mg carbonates with δ¹³C values in the waters being uncharacteristically high (approx. −2‰ at the end of the season), reflecting a largely open system in which the CO₂ is derived equally from the atmosphere and plant/soil sources. Early melt waters had δ³⁴S values dominated by sea salt reflecting the close proximity to the ocean. However, as the season progressed the marine influence lessened. Extrapolation of the data suggests that the origin of non-sea salt δ³⁴S was the oxidation of reduced sulfur from coal particles in the subsoil. Concentrations of inorganic N in stream waters were generally very low. However, values were found to increase as the season progressed, possibly through increased microbial activity in the soil and the early senescence of tundra plants reducing demand. Dual isotope analysis of δ¹⁵N and δ¹⁸O suggested that the vast majority of snow-pack was assimilated by the soil microbial biomass before being released, recycled and lost to drainage waters. Organic N concentrations in drainage waters were generally equal to or greater than losses of inorganic N from tundra soils. The study demonstrated the effectiveness of stable isotope data for understanding biogeochemical cycling and soil-water interactions in tundra ecosystems. The implications of the results are discussed in relation to climate warming.     

Controls on dissolved inorganic carbon and δC in cave waters from DeSoto Caverns: Implications for speleothem δC assessments

Unraveling the factors controlling the carbon chemistry and transport of carbon within extant karst systems has important implications concerning the assessment of time-series δ¹³C records of speleothems. Here we report the results of a 3-year study of total dissolved inorganic carbon [DIC] and δ¹³C_DIC from cave waters at DeSoto Caverns (Southeastern USA) that offer valuable insight on carbon transport and the accompanied isotope fractionations from end-member sources to speleothems.[DIC] and δ¹³C_DIC values of cave waters range from 0.2 to 6.0 mM and 2.7 to −12.9 (‰ VPDB), respectively. [DIC] and δ¹³C_DIC of “seasonal drips” show seasonal, albeit noisy, variability and are inversely related (δ¹³C_DIC = −2.49[DIC] + 0.64, r² = 0.84). A shallow pool fed by multiple drips shows a bimodal δ¹³C_DIC distribution with an isotopically heavier mode during winter (−4‰ to −5‰ VPDB) relative to summer months (−9‰ to −10‰ VPDB). A multi-year trend of decreasing water availability during the study period is not reflected in a response of cave water carbon chemistry suggesting that rainfall amount may not be a significant controlling factor of the carbon chemistry. Coupled cave air winter ventilation/summer stagnation and varying CO₂ fluxes through the soil horizon and epikarst exert the strongest influence on seasonal [DIC] and δ¹³C_DIC variability. Measured values of high [DIC] and low δ¹³C_DIC from cave waters collected during the summer/early fall closely approximate isotopic equilibrium conditions. Conversely, low [DIC] and high δ¹³C_DIC values during winter/early months indicate kinetically enhanced isotopic fractionations within the cave waters. The kinetically enhanced isotopic fractionation of partitioned between degassed CO₂ and precipitated CaCO₃(1000lnα_{[(CO2-HCO3)+(CaCO3(AR)-HCO3)]/2}) is greater by about a factor of two (−6.7 ± 0.3‰) relative to the same isotopic fractionation under equilibrium conditions (−3.1‰).On the basis of ¹⁴C mass balance and paired ¹⁴C-U/Th measurements we estimate that on average about ∼23% of C delivered annually by the drips to the aragonite stalagmites is derived from ¹⁴C-dead dolomite cap while the remainder of ∼77% is derived from ¹⁴C-live biomass. δ¹³C measurements of aragonite (n = 12) sampled from the tips of active speleothems during the summer months are consistent with theoretical aragonite δ¹³C values calculated using the shallow pool summer/early fall data thus confirming the δ¹³C seasonality in both drips and coeval aragonite. δ¹³C values of an active stalagmite section spanning the last 200 years show a normal distribution with a mean of −7.1 ± 1.2‰ (n = 81) and a mode of −7‰ to −8‰ that are statistically indistinguishable from the annual mean and mode of all dripwaters. Thus secular time-series δ¹³C records of stalagmites at DeSoto Caverns with resolving power >10⁻¹ year will likely carry the imprints of drip annual means that record climate-driven δ¹³C seasonal biases.     

Open to closed system transition traced through the TDIC isotopic signature at the aquifer recharge stage, implications for groundwater C dating

¹⁴C dating models are limited when considering recent groundwater for which the carbon isotopic signature of the total dissolved inorganic carbon (TDIC) is mainly acquired in the unsaturated zone. Reducing the uncertainties of dating thus implies a better identification of the processes controlling the carbon isotopic composition of the TDIC during groundwater recharge. Geochemical interactions between gas, water and carbonates in the unsaturated zone were investigated for two aquifers (the carbonate-free Fontainebleau sands and carbonate-bearing Astian sands, France) in order to identify the respective roles of CO₂ and carbonates on the carbon isotopic signatures of the TDIC; this analysis is usually approached using open or closed system terms. Under fully open system conditions, the seasonality of the ¹³C values in the soil CO₂ can lead to important uncertainties regarding the so-called “initial ¹⁴C activity” used in ¹⁴C correction models. In a carbonate-bearing unsaturated zone such as in the Astian aquifer, we show that an approach based on fully open or closed system conditions is not appropriate. Although the chemical saturation between water and calcite occurs rapidly within the first metre of the unsaturated zone, the carbon isotopic contents (δ¹³C) of the CO₂ and the TDIC evolve downward, impacted by the dissolution-precipitation of the carbonates. In this study, we propose a numerical approach to describe this evolution. The δ¹³C and the A¹⁴C (radiocarbon activity) of the TDIC at the base of the carbonate-bearing unsaturated zone depends on (i) the δ¹³C and the A¹⁴C of the TDIC in the soil determined by the soil CO_2, (ii) the water’s residence time in the unsaturated zone and (iii) the carbonate precipitation-dissolution fluxes. In this type of situation, the carbonate δ¹³C-A¹⁴C evolutions indicate the presence of secondary calcite and permit the calculation of its accretion flux, equal to . More generally, for other sites under temperate climate and with similar properties to the Astian sands site, this approach allows for a reliable determination of the carbon isotopic composition at the base of the unsaturated zone as the indispensable “input function” data of the carbon cycle into the aquifer.     

Carbonate biomineralization induced by soil bacterium Bacillus megaterium

Biogenic carbonates spawned from microbial activities are common occurrences in soils. Here, we investigate the carbonate biomineralization mediated by the bacterium Bacillus megaterium, a dominant strain separated from a loess profile in China. Upon completing bacterial cultivation, the ensuring products are centrifuged, and the resultant supernatant and the concentrated bacterial sludge as well as the un-separated culture are added separately into a Ca-CO₃ containing solution for crystallization experiments. Results of XRD and SEM analysis indicate that calcite is the dominant mineral phase formed when the bacteria are present. When the supernatant alone is used, however, a significant portion of vaterite is also precipitated. Experimental results further reveal that the bacteria have a strong tendency to colonize the center area of the calcite {1 0  4} faces. Observed crystal morphology suggests that the bacterial colony may promote the growth normal to each individual {1 0  4} face of calcite when the cell concentration is high, but may retard it or even cause dissolution of the immediate substrate surfaces when the concentration is low. SEM images taken at earlier stages of the crystallization experiments demonstrate the nucleation of calcite on the bacterial cell walls but do not show obvious morphological changes on the nanometer- to submicron-sized nuclei. δ¹³C measurements unveil that the crystals grown in the presence of bacteria are further enriched in the heavy carbon isotope, implying that the bacterial metabolism may not be the carbon sources for the mineralization. Based upon these findings, we propose a mechanism for the B. megaterium mediated calcite mineralization and conclude that the whole process involves epi- and inter-cellular growth in the local microenvironments whose conditions may be controlled by cell sequestration and proton pumping during bacterial respiration.     

Effects of magnesium ions on near-equilibrium calcite dissolution: Step kinetics and morphology

Dissolution kinetics at the aqueous solution-calcite interface at 50 °C were investigated using in situ atomic force microscopy (AFM) to reveal the influence of magnesium concentration and solution saturation state on calcite dissolution kinetics and surface morphology. Under near-equilibrium conditions, dissolved Mg²⁺ displayed negligible inhibitory effects on calcite dissolution even at concentrations of . Upon the introduction of , the solution saturation state with respect to calcite, , acted as a “switch” for magnesium inhibition whereby no significant changes in step kinetics were observed at Ω_calcite<0.2, whereas a sudden inhibition from Mg²⁺ was activated at Ω_calcite?0.2. The presence of the Ω-switch in dissolution kinetics indicates the presence of critical undersaturation in accordance with thermodynamic principles. The etch pits formed in solutions with exhibited a unique distorted rhombic profile, different from those formed in Mg-free solutions and in de-ionized water. Such unique etch pit morphology may be associated with the anisotropy in net detachment rates of counter-propagating kink sites upon the addition of Mg²⁺.     

Isotopic fractionations associated with phosphoric acid digestion of carbonate minerals: Insights from first-principles theoretical modeling and clumped isotope measurements

Phosphoric acid digestion has been used for oxygen- and carbon-isotope analysis of carbonate minerals since 1950, and was recently established as a method for carbonate ‘clumped isotope’ analysis. The CO₂ recovered from this reaction has an oxygen isotope composition substantially different from reactant carbonate, by an amount that varies with temperature of reaction and carbonate chemistry. Here, we present a theoretical model of the kinetic isotope effects associated with phosphoric acid digestion of carbonates, based on structural arguments that the key step in the reaction is disproportionation of H₂CO₃ reaction intermediary. We test that model against previous experimental constraints on the magnitudes and temperature dependences of these oxygen isotope fractionations, and against new experimental determinations of the fractionation of ¹³C-¹⁸O-containing isotopologues (‘clumped’ isotopic species). Our model predicts that the isotope fractionations associated with phosphoric acid digestion of carbonates at 25 °C are 10.72‰, 0.220‰, 0.137‰, 0.593‰ for, respectively, ¹⁸O/¹⁶O ratios (1000 lnα^∗) and three indices that measure proportions of multiply-substituted isotopologues . We also predict that oxygen isotope fractionations follow the mass dependence exponent, λ of 0.5281 (where ). These predictions compare favorably to independent experimental constraints for phosphoric acid digestion of calcite, including our new data for fractionations of ¹³C-¹⁸O bonds (the measured change in Δ₄₇ = 0.23‰) during phosphoric acid digestion of calcite at 25 °C.We have also attempted to evaluate the effect of carbonate cation compositions on phosphoric acid digestion fractionations using cluster models in which disproportionating H₂CO₃ interacts with adjacent cations. These models underestimate the magnitude of isotope fractionations and so must be regarded as unsucsessful, but do reproduce the general trend of variations and temperature dependences of oxygen isotope acid digestion fractionations among different carbonate minerals. We suggest these results present a useful starting point for future, more sophisticated models of the reacting carbonate/acid interface. Examinations of these theoretical predictions and available experimental data suggest cation radius is the most important factor governing the variations of isotope fractionation among different carbonate minerals. We predict a negative correlation between acid digestion fractionation of oxygen isotopes and of ¹³C-¹⁸O doubly-substituted isotopologues, and use this relationship to estimate the acid digestion fractionation of for different carbonate minerals. Combined with previous theoretical evaluations of ¹³C-¹⁸O clumping effects in carbonate minerals, this enables us to predict the temperature calibration relationship for different carbonate clumped isotope thermometers (witherite, calcite, aragonite, dolomite and magnesite), and to compare these predictions with available experimental determinations. The success of our models in capturing several of the features of isotope fractionation during acid digestion supports our hypothesis that phosphoric acid digestion of carbonate minerals involves disproportionation of transition state structures containing H₂CO₃.     

Reconstructing palaeoprecipitation from an active cave flowstone

Several drill cores were obtained from a laminated, actively forming flowstone from a shallow cave in Austria. Highly resolved petrographic and geochemical analyses combined with multi‐annual cave monitoring reveal a distinct sensitivity of flowstone growth and composition to late Holocene meteoric precipitation. The regular sub‐millimetre‐scale lamination consists of thicker, translucent laminae and thinner (organic) inclusion‐rich laminae. There is also a macroscopic millimetre‐scale banding of darker and lighter bands comprising several laminae. Stable isotope analyses of drill cores and modern calcite precipitates show a pronounced positive covariation of δ¹³C and δ¹⁸O values indicative of kinetic isotope effects. Comparing the isotope values with petrography shows gradual changes across several of the annual laminae, i.e. changes of several per mille on a multi‐annual to decadal timescale. The stable isotope and trace‐element composition, as well as the flowstone petrography, are mainly controlled by the variable drip‐water discharge controlling the water‐film thickness and water residence time on the flowstone surface and consequently the intensity of CO₂‐degassing, kinetic isotope enrichment and concomitant calcite precipitation. Drill core PFU6 provides an isotope record of the last ca. 3000 years at near‐annual resolution. A distinct phase of low C and O isotope values – interpreted as increased discharge and hence higher meteoric precipitation – occurred from ca. 300 to 140 a b2k (second half of the Little Ice Age) and another wet interval occurred around 700 a, corresponding to reported Medieval glacier advances. The Roman Warm Period was also dominated by relatively wet conditions, although significant decadal variability prevailed. Increased precipitation further characterized the intervals from ca. 2480 to 2430 and 2950 to 2770 a. Dry conditions persisted during the Medieval Climate Anomaly, although this trend towards reduced precipitation started earlier. The highest C isotope values of the last 2 ka are recorded around 750 a and another dry phase is centred at 1480 a. This new record shows that inter‐annual to decadal oscillations are a dominant mode of variation during the last 3 ka in the Alps. Copyright © 2011 John Wiley & Sons, Ltd.     

Carbon isotope fractionation in synthetic magnesian calcite

Mg-calcite was precipitated at 25 °C in closed system, free-drift experiments, from solutions containing NaHCO₃, CaCl₂ and MgCl₂. The carbon stable isotope composition of bulk solid and solution were analyzed from subsamples collected during time course experiments of 24 h duration. Considering only the Mg-content and δ¹³C values for the bulk solid, the carbon isotope fractionation factor for the Mg-calcite-HCO_3(aq)⁻ system (as ) increased with average mol percentage of Mg (X_Mg) in the solid at a rate of (0.024 ± 0.011) per mol% MgCO₃. Extrapolation of this relationship to the pure calcite end member yields a value of 0.82 ± 0.09, which is similar to published values for the calcite-HCO_3(aq)⁻ system. Although did not vary for precipitation rates that ranged from 10^3.21 to 10^4.60 μmol m⁻² h⁻¹, it was not possible to hold Mg-content of the solid constant, so kinetic effect on 10³ ln α could not be evaluated from these experiments.     

贵州七星洞系统中水文地球化学特征对滴水δ~(13)C_(DIC)的影响及其意义

2003年4月至2004年5月,笔者对贵州七星洞(QXD)进行了较为详细的监测,逐月采集了土壤水和洞穴滴水等样品,分别测定了样品的稳定碳同位素组成和水文地球化学参数。结果显示,9个滴水点同期的溶解无机碳同位素值(δ13CDIC)之间存在着大的差异,最大达6.9‰;δ13CDIC值偏重的Ⅰ组滴水,其Ca、Sr、HCO3、电导率(EC)和方解石饱和指数(SIC)等水文地球化学指标偏小,而Mg/Ca比值偏大,Ⅱ组则相反;不同滴水点的δ13CDIC值分别与相应滴水的Ca、Sr、HCO3、EC、Mg/Ca和SIC等水文地球化学指标之间存在较好的相关关系。进一步分析表明,与土壤水δ13CDIC平均值-9.9‰相比,Ⅰ组1#、2#、6#、7#和8#滴水点δ13CDIC值偏重4.5‰～5.7‰,主要是由大量的基岩溶解以及前期方解石沉积(prior calcite precipitation,PCP)共同作用的结果;Ⅱ组3#、4#、5#和9#滴水点δ13CDIC值偏重0.6‰～1.6‰,受基岩溶解和PCP过程影响较小,尤其是9#滴水点受影响最小。因此,若不考虑基岩溶解和PCP过程影响作用,将会极大地影响洞穴化学沉积物碳同位素记录的准确解释。