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     

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.     

Intra-shell oxygen isotope ratios in the benthic foraminifera genus Amphistegina and the influence of seawater carbonate chemistry and temperature on this ratio

Using secondary ion mass spectrometry (SIMS) we looked at the natural variability in the oxygen isotope ratio of the shallow water, symbionts-bearing foraminiferan Amphistegina lobifera. Live foraminifera were collected in February 2005 in the Gulf of Eilat, Israel. Vertical section exposing the knob area of this species represents the growth history of this species from August 2004 to February 2005. SIMS profile at a resolution of ∼15 μm (representing about 2 weeks considering the size of the knob area and the life span of ≈6 months of this foraminifera species) yielded δ¹⁸O changes of ∼1.5‰ that are compatible with the known temperature changes for the Gulf of Eilat for this period (21-27 °C). Natural variability between primary and secondary calcite at the knob area were obtained on horizontal section of the upper knob area. This section is semi-tangential to the growth lines and exposes relatively wide belts of the primary calcite which could be analysed using the SIMS (beam size of 10 × 20 μm). The primary calcite δ¹⁸O value is on average more than 3‰ lower than the secondary calcite that represents the bulk of the skeleton (more than 95% by weight). A vertical profile at the knob was obtained by rastering an area of 50 × 50 μm at vertical steps of roughly 1 μm. The profile revealed a narrow zone of lower δ¹⁸O compared to the higher values above and below it. The difference between the lowest δ¹⁸O and the highest one was also close to 2‰. The δ¹⁸O in the margin - keel area of A. lobifera is also lower compared to the bulk secondary calcite. Specimens that were cultured in the laboratory at a constant temperature and inorganic carbon but at different pH have increased their CaCO₃ weight by roughly a factor of 8. Single specimen from each pH (ranging between 7.90 and 8.45) were investigated with the SIMS at the knob area. While there is some variability within each specimen (perhaps related to the primary calcite), the general trend was a decrease in δ¹⁸O with increasing pH (or CO₃²⁻ concentration), in agreement with previous studies on planktonic foraminifera. Some other specimens grown at different temperatures (between 21 and 33 °C) were also measured with the SIMS at the knob area. For each temperature, we observed also some variability, nevertheless the trend of −0.2‰/°C in δ¹⁸O is observed.     

Stable isotope variations in modern tropical speleothems: Evaluating equilibrium vs. kinetic isotope effects

Applications of speleothem calcite geochemistry in climate change studies require the evaluation of the accuracy and sensitivity of speleothem proxies to correctly infer paleoclimatic information. The present study of Harrison’s Cave, Barbados, uses the analysis of the modern climatology and groundwater system to evaluate controls on the C and O isotopic composition of modern speleothems. This new approach directly compares the δ¹⁸O and δ¹³C values of modern speleothems with the values for their corresponding drip waters in order to assess the degree to which isotopic equilibrium is achieved during calcite precipitation. If modern speleothems can be demonstrated to precipitate in isotopic equilibrium, then ancient speleothems, suitable for paleoclimatic studies, from the same cave environment may also have been precipitated in isotopic equilibrium. If modern speleothems are precipitated out of isotopic equilibrium, then the magnitude and direction of the C and O isotopic offsets may allow specific kinetic and/or equilibrium isotopic fractionation mechanisms to be identified.Carbon isotope values for the majority of modern speleothem samples from Harrison’s Cave fall within the range of equilibrium values predicted from the combined use of (1) calcite-water fractionation factors from the literature, (2) measured temperatures, and (3) measured δ¹³C values of the dissolved inorganic carbon of drip waters. Calcite samples range from ∼0.8‰ higher to ∼1.1‰ lower than predicted values. The ¹³C depletions are likely caused by kinetically driven departures in the fractionation between HCO₃⁻ (aq) and CaCO₃ from equilibrium conditions, caused by rapid calcite growth. ¹³C enrichments can be accounted for by Rayleigh distillation of the HCO₃⁻ (aq) reservoir during degassing of ¹³C-depleted CO₂.Modern speleothems from Harrison’s Cave are not in O isotopic equilibrium with their corresponding drip waters and are 0.2‰ to 2.3‰ enriched in ¹⁸O relative to equilibrium values. δ¹⁸O variations in modern calcite are likely controlled by kinetically driven changes in the fractionation between HCO₃⁻ (aq) and CaCO₃ from equilibrium conditions to nonequilibrium conditions, consistent with rapid calcite growth. In contrast to δ¹³C, δ¹⁸O values of modern calcite may not be affected by Rayleigh distillation during degassing because CO₂ hydration and hydroxylation reactions will buffer the O isotopic composition of the HCO₃⁻ (aq) reservoir. If the effects of Rayleigh distillation manifest themselves in the O isotopic system, they will result in ¹⁸O enrichment in the HCO₃⁻ (aq) reservoir and ultimately in the precipitated CaCO₃.     

C-O bonds in carbonate minerals: A new kind of paleothermometer

The abundance of the doubly substituted CO₂ isotopologue, ¹³C¹⁸O¹⁶O, in CO₂ produced by phosphoric acid digestion of synthetic, inorganic calcite and natural, biogenic aragonite is proportional to the concentration of ¹³C-¹⁸O bonds in reactant carbonate, and the concentration of these bonds is a function of the temperature of carbonate growth. This proportionality can be described between 1 and 50 °C by the function: Δ₄₇ = 0.0592 · 10⁶ · T⁻² − 0.02, where Δ₄₇ is the enrichment, in per mil, of ¹³C¹⁸O¹⁶O in CO₂ relative to the amount expected for a stochastic (random) distribution of isotopes among all CO₂ isotopologues, and T is the temperature in Kelvin. This relationship can be used for a new kind of carbonate paleothermometry, where the temperature-dependent property of interest is the state of ordering of ¹³C and ¹⁸O in the carbonate lattice (i.e., bound together vs. separated into different CO₃²⁻ units), and not the bulk δ¹⁸O or δ¹³C values. Current analytical methods limit precision of this thermometer to ca. ± 2 °C, 1σ. A key feature of this thermometer is that it is thermodynamically based, like the traditional carbonate-water paleothermometer, and so is suitable for interpolation and even modest extrapolation, yet is rigorously independent of the δ¹⁸O of water and δ¹³C of DIC from which carbonate grew. Thus, this technique can be applied to parts of the geological record where the stable isotope compositions of waters are unknown. Moreover, simultaneous determinations of Δ₄₇ and δ¹⁸O for carbonates will constrain the δ¹⁸O of water from which they grew.     

Calibration of the calcite-water oxygen-isotope geothermometer at Devils Hole, Nevada, a natural laboratory

The δ¹⁸O of ground water (−13.54 ± 0.05 ‰) and inorganically precipitated Holocene vein calcite (+14.56 ± 0.03 ‰) from Devils Hole cave #2 in southcentral Nevada yield an oxygen isotopic fractionation factor between calcite and water at 33.7 °C of 1.02849 ± 0.00013 (1000 ln α_{calcite-water} = 28.09 ± 0.13). Using the commonly accepted value of ∂(α_{calcite-water})/∂T of −0.00020 K⁻¹, this corresponds to a 1000 ln α_{calcite-water} value at 25 °C of 29.80, which differs substantially from the current accepted value of 28.3. Use of previously published oxygen isotopic fractionation factors would yield a calcite precipitation temperature in Devils Hole that is 8 °C lower than the measured ground water temperature. Alternatively, previously published fractionation factors would yield a δ¹⁸O of water, from which the calcite precipitated, that is too negative by 1.5 ‰ using a temperature of 33.7 °C. Several lines of evidence indicate that the geochemical environment of Devils Hole has been remarkably constant for at least 10 ka. Accordingly, a re-evaluation of calcite-water oxygen isotopic fractionation factor may be in order.Assuming the Devils Hole oxygen isotopic value of α_{calcite-water} represents thermodynamic equilibrium, many marine carbonates are precipitated with a δ¹⁸O value that is too low, apparently due to a kinetic isotopic fractionation that preferentially enriches ¹⁶O in the solid carbonate over ¹⁸O, feigning oxygen isotopic equilibrium.     

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.     

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

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

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.     

The record of temperature, wind velocity and air humidity in the δD and δO of water inclusions in synthetic and Messinian halites

Deuterium and oxygen isotope fractionations between liquid and vapor water were experimentally-determined during evaporation of a NaCl solution (35 g L⁻¹) as a function of water temperature and wind velocity. In the case of a null wind velocity, slopes of δD-δ¹⁸O trajectories of residual waters hyperbolically decrease with increasing water temperatures in the range 23-47 °C. For wind velocities ranging from 0.8 to 2.2 m s⁻¹, slopes of the δD-δ¹⁸O trajectories linearly increase with increasing wind velocity at a given water temperature. These experimental results can be modeled by using Rayleigh distillation equations taking into account wind-related kinetics effects. Deuterium and oxygen isotope compositions of water inclusions trapped by the precipitated halite crystals were determined by micro-equilibration techniques.These isotopic compositions accurately reflect those of the surrounding residual waters during halite growth. Isotopic compositions of water inclusions in twenty natural halites from the Messinian Realmonte mine in Sicily suggest precipitation temperatures of that match the homogenization temperatures obtained by microthermometry (median = 34 ± 5 °C). The similarity between the measured and experimental slopes of the δD-δ¹⁸O evaporation trajectories suggests that the effect of wind was negligible during the genesis of these halite deposits. Hydrogen and oxygen isotope compositions of water inclusions from Realmonte halite also define a linear trend whose extrapolation until intersection with the Mediterranean Meteoric Water Line allows the characterization of the water source with δD and δ¹⁸O values of −70 ± 10‰ and −11.5 ± 1.5‰, respectively. These results reveal that the huge amounts of salts deposited in Sicily result from the evaporation of seawater mixed with a dominant fraction (?50%) of meteoric waters most likely deriving from alpine fluvial discharge.     

O/O and D/H ratios of pedogenic kaolinite in a North American Cenomanian laterite: Paleoclimatic implications

Kaolinite, gibbsite and quartz are the dominant minerals in samples collected from two outcrops of a Cenomanian (∼95 Ma) laterite in southwestern Minnesota. A combination of measured yields and isotope ratios permitted mass balance calculations of the δD and δ¹⁸O values of the kaolinite in these samples. These calculations yielded kaolinite δD values of about −73‰ and δ¹⁸O values of about +18.7‰. The δD and δ¹⁸O values appear to preserve information on the ancient weathering system.If formed in hydrogen and oxygen isotope equilibrium with water characterized by the global meteoric water line (GMWL), the kaolinite δD and δ¹⁸O values indicate a crystallization temperature of 22 (±5) °C. A nominal paleotemperature of 22 °C implies a δ¹⁸O value for the corresponding water of −6.3‰. The combination of temperature and meteoric water δ¹⁸O values is consistent with relatively intense rainfall at that mid-paleolatitude location (∼40°N) on the eastern shore of the North American Western Interior Seaway. The inferred Cenomanian paleosol temperature of ∼22 °C is in general accord with published mid-Cretaceous continental mean annual temperatures (MAT) estimated from leaf margin analyses of fossil plants.When compared with results from a published GCM-based Cenomanian climate simulation which specifies a latitudinal sea surface temperature (SST) gradient that was either near modern or smaller-than-modern, the kaolinite paleotemperature of 22 °C is closer to the GCM-predicted MAT for a smaller equator-to-pole temperature difference in the mid-Cretaceous. Moreover, the warm, kaolinite-derived, mid-paleolatitude temperature of 22 °C is associated with proxy estimates of high concentrations of atmospheric CO₂ in the Cenomanian. The overall similarity of proxy and model results suggests that the general features of Cenomanian continental climate in that North American locale are probably being revealed.     

Calcite dissolution kinetics in Na-Ca-Mg-Cl brines

Sedimentary basins can contain close to 20% by volume of pore fluids commonly classified as brines. These fluids can become undersaturated with respect to calcite as a result of migration, dispersive mixing, or anthropogenic injection of CO₂. This study measured calcite dissolution rates in geologically relevant Na-Ca-Mg-Cl synthetic brines (50-200 g L⁻¹ TDS). The dissolution rate dependency on brine composition, pCO₂ (0.1-1 bar), and temperature (25.0-82.5 °C) was modeled using the empirical rate equation

R=k(1-Ω)ⁿ     

Effects of aridity and vegetation on plant-wax δD in modern lake sediments

We analyzed the deuterium composition of individual plant-waxes in lake sediments from 28 watersheds that span a range of precipitation D/H, vegetation types and climates. The apparent isotopic fractionation (ε_a) between plant-wax n-alkanes and precipitation differs with watershed ecosystem type and structure, and decreases with increasing regional aridity as measured by enrichment of ²H and ¹⁸O associated with evaporation of lake waters. The most negative ε_a values represent signatures least affected by aridity; these values were −125 ± 5‰ for tropical evergreen and dry forests, −130‰ for a temperate broadleaf forest, −120 ± 9‰ for the high-altitude tropical páramo (herbs, shrubs and grasses), and −98 ± 6‰ for North American montane gymnosperm forests. Minimum ε_a values reflect ecosystem-dependent differences in leaf water enrichment and soil evaporation. Slopes of lipid/lake water isotopic enrichments differ slightly with ecosystem structure (i.e. open shrublands versus forests) and overall are quite small (slopes = 0-2), indicating low sensitivity of lipid δD variations to aridity compared with coexisting lake waters. This finding provides an approach for reconstructing ancient precipitation signatures based on plant-wax δD measurements and independent proxies for lake water changes with regional aridity. To illustrate this approach, we employ paired plant-wax δD and carbonate-δ¹⁸O measurements on lake sediments to estimate the isotopic composition of Miocene precipitation on the Tibetan plateau.     

Oxygen and carbon isotopic systematics of aragonite speleothems and water in Furong Cave, Chongqing, China

To understand oxygen and carbon stable isotopic characteristics of aragonite stalagmites and evaluate their applicability to paleoclimate, the isotopic compositions of active and fossil aragonite speleothems and water samples from an in situ multi-year (October 2005-July 2010) monitoring program in Furong Cave located in Chongqing of China have been examined. The observations during October 2005-June 2007 show that the meteoric water is well mixed in the overlying 300-500-m bedrock aquifer, reflected by relatively constant δ¹⁸O, ±0.11-0.14‰ (1σ), of drip waters in the cave, which represents the annual status of rainfall water. Active cave aragonite speleothems are at oxygen isotopic equilibrium with drip water and their δ¹⁸O values capture the surface-water oxygen isotopic signal. Aragonite-to-calcite transformation since the last glaciation is not noticeable in Furong stalagmites. Our multi-year field experiment approves that aragonite stalagmite δ¹⁸O records in this cave are suitable for paleoclimate reconstruction. With high U, 0.5-7.2 ppm, and low Th, 20-1270 ppt, the Furong aragonite stalagmites provide very precise chronology (as good as ±20s yrs (2σ)) of the climatic variations since the last deglaciation. The synchroneity of Chinese stalagmite δ¹⁸O records at the transition into the Bølling-Allerød (t-BA) and the Younger Dryas from Furong, Hulu and Dongge Caves supports the fidelity of the reconstructed East Asian monsoon evolution. However, the Furong record shows that the cold Older Dryas (OD) occurred at 14.0 thousand years ago, agreeing with Greenland ice core δ¹⁸O records but ∼200 yrs younger than that in the Hulu record. The OD age discrepancy between Chinese caves can be attributable to different regionally climatic/environmental conditions or chronological uncertainty of stalagmite proxy records, which is limited by changes in growth rate and subsampling intervals in absolute dating. Seasonal dissolved inorganic carbon δ¹³C variations of 2-3‰ in the drip water and 5-7‰ in the pool and spring waters are likely attributed to variable degrees of CO₂ degassing in winter and summer. The variable δ¹³C values of active deposits from −11‰ to 0‰ could be caused by kinetically mediated CO₂ degassing processes. The complicated nature of pre-deposition kinetic isotopic fractionation processes for carbon isotopes in speleothems at Furong Cave require further study before they can be interpreted in a paleoclimatic or paleoenvironmental context.     

Speleothem record of the last 180 ka in Villars cave (SW France): Investigation of a large δ18O shift between MIS6 and MIS5

The Vil-car-1 flowstone core from Villars cave (SW France) provides one of the first European speleothem records extending back to 180 ka, based on U–Th TIMS and MC-ICP-MS measurements. The core offers a continuous record of Termination II and the Last Interglacial. The penultimate deglaciation is characterized by a prominent 5‰ depletion in calcite δ¹⁸O. Determining which specific environmental factors controlled such a large oxygen isotopic shift offers the opportunity to assess the impact of various factors influencing δ¹⁸O variations in speleothem calcite.Oxygen isotope analyses of fluid inclusions indicate that drip water δ¹⁸O remained within a very narrow range of ±1‰ from Late MIS6 to the MIS5 δ¹⁸O optimum. The possibility of such a stable behaviour is supported by simple calculations of various effects influencing seepage water δ¹⁸O.Although this could suggest that the isotopic shift in calcite is mainly driven by temperature increase, attempts to quantify the temperature shift from Late MIS6 to the MIS5 δ¹⁸O optimum by assuming an equilibrium relationship between calcite and fluid inclusion δ¹⁸O yield unreasonably high estimates of ～20 °C warming and Late MIS6 cave temperatures below 0 °C; this suggests that the flowstone calcite precipitated out of thermodynamic equilibrium at this site.Using a method proposed by Guo et al. (submitted for publication) combining clumped isotope measurements, fluid inclusion and modern calcite δ¹⁸O analyses, it is possible to quantitatively correct for isotopic disequilibrium and estimate absolute paleotemperatures. Although the precision of these absolute temperature reconstructions is limited by analytical uncertainties, the temperature rise between Late MIS6 and the MIS5 optimum can be robustly constrained between 13.2 ± 2.6 and 14.6 ± 2.6 °C (1σ), consistent with existing estimates from Western Europe pollen and sea-surface temperature records.     

Characterization of calcium isotopes in natural and synthetic barite

The mineral barite (BaSO₄) accommodates calcium in its crystal lattice, providing an archive of Ca-isotopes in the highly stable sulfate mineral. Holocene marine (pelagic) barite samples from the major ocean basins are isotopically indistinguishable from each other (δ^44/40Ca = −2.01 ± 0.15‰) but are different from hydrothermal and cold seep barite samples (δ^44/40Ca = −4.13 to −2.72‰). Laboratory precipitated (synthetic) barite samples are more depleted in the heavy Ca-isotopes than pelagic marine barite and span a range of Ca-isotope compositions, Δ^44/40Ca = −3.42 to −2.40‰. Temperature, saturation state, , and aCa²⁺/aBa²⁺ each influence the fractionation of Ca-isotopes in synthetic barite; however, the fractionation in marine barite samples is not strongly related to any measured environmental parameter. First-principles lattice dynamical modeling predicts that at equilibrium Ca-substituted barite will have much lower ⁴⁴Ca/⁴⁰Ca than calcite, by −9‰ at 0 °C and −8‰ at 25 °C. Based on this model, none of the measured barite samples appear to be in isotopic equilibrium with their parent solutions, although as predicted they do record lower δ^44/40Ca values than seawater and calcite. Kinetic fractionation processes therefore most likely control the extent of isotopic fractionation exhibited in barite. Potential fractionation mechanisms include factors influencing Ca²⁺ substitution for Ba²⁺ in barite (e.g. ionic strength and trace element concentration of the solution, competing complexation reactions, precipitation or growth rate, temperature, pressure, and saturation state) as well as nucleation and crystal growth rates. These factors should be considered when investigating controls on isotopic fractionation of Ca²⁺ and other elements in inorganic and biogenic minerals.     

A laser fluorination method for oxygen isotope analysis of biogenic silica and a new oxygen isotope calibration of modern diatoms in freshwater environments

An empirical calibration for the oxygen isotope fractionation between biogenic silica and water was determined for diatom frustules sampled from living diatom communities in the Jemez Mountains of northern New Mexico, USA. Over a temperature range from 5.1 to 37.8 °C, the silica-water fractionation is defined by the equation 1000 ln α_{(silica-water)} = 2.39(±0.13) × 10⁶T⁻² + 4.23(±1.49). This relationship is in close agreement with other published silica-water fractionation factors for laboratory cultured diatom samples; however, it is as much as 8‰ lower than equilibrium quartz-water fractionations and 3-4‰ lower than observed silica-water fractionations in diatomaceous silica collected from sediment traps and sediment cores. There are three possible explanations for the disparate silica-water fractionation factors observed in diatom silica: (1) silica does not precipitate in equilibrium with ambient water, (2) silica does precipitate in equilibrium with ambient water, but the silica-water fractionation factor for diatom silica is considerably less than the equilibrium fractionation factor for quartz-water, or (3) silica precipitation is influenced by a ‘vital’ effect, where the δ¹⁸O value of the water inside the diatom cell walls is lower than the δ¹⁸O values of ambient water.Post-mortem loss of organic material results in an alteration or ‘maturation’ of diatom silica in which silica reequilibrates with a silica-water fractionation closer to the equilibrium quartz-water fractionation. Alteration is likely to occur rapidly after the diatom frustule loses its organic coating, either as it settles through the water column or at the sediment-water interface; δ¹⁸O values recorded by paleo-diatom silica therefore do not record growing conditions but more likely record conditions at the sediment-water interface. In the case of lacustrine environments, where the bottom water remains at a nearly constant 4 °C, the reequilibration of diatom silica with bottom conditions could reduce or remove the conflating effects of temperature on δ¹⁸O values recorded by paleo-diatom silica and provide direct information on the δ¹⁸O value of the lake water.     

Iron and magnesium isotopic compositions of peridotite xenoliths from Eastern China

We present high-precision measurements of Mg and Fe isotopic compositions of olivine, orthopyroxene (opx), and clinopyroxene (cpx) for 18 lherzolite xenoliths from east central China and provide the first combined Fe and Mg isotopic study of the upper mantle. δ⁵⁶Fe in olivines varies from 0.18‰ to −0.22‰ with an average of −0.01 ± 0.18‰ (2SD, n = 18), opx from 0.24‰ to −0.22‰ with an average of 0.04 ± 0.20‰, and cpx from 0.24‰ to −0.16‰ with an average of 0.10 ± 0.19‰. δ²⁶Mg of olivines varies from −0.25‰ to −0.42‰ with an average of −0.34 ± 0.10‰ (2SD, n = 18), opx from −0.19‰ to −0.34‰ with an average of −0.25 ± 0.10‰, and cpx from −0.09‰ to −0.43‰ with an average of −0.24 ± 0.18‰. Although current precision (∼±0.06‰ for δ⁵⁶Fe; ±0.10‰ for δ²⁶Mg, 2SD) limits the ability to analytically distinguish inter-mineral isotopic fractionations, systematic behavior of inter-mineral fractionation for both Fe and Mg is statistically observed: Δ⁵⁶Fe_ol-cpx = −0.10 ± 0.12‰ (2SD, n = 18); Δ⁵⁶Fe_ol-opx = −0.05 ± 0.11‰; Δ²⁶Mg_ol-opx = −0.09 ± 0.12‰; Δ²⁶Mg_ol-cpx = −0.10 ± 0.15‰. Fe and Mg isotopic composition of bulk rocks were calculated based on the modes of olivine, opx, and cpx. The average δ⁵⁶Fe of peridotites in this study is 0.01 ± 0.17‰ (2SD, n = 18), similar to the values of chondrites but slightly lower than mid-ocean ridge basalts (MORB) and oceanic island basalts (OIB). The average δ²⁶Mg is −0.30 ± 0.09‰, indistinguishable from chondrites, MORB, and OIB. Our data support the conclusion that the bulk silicate Earth (BSE) has chondritic δ⁵⁶Fe and δ²⁶Mg.The origin of inter-mineral fractionations of Fe and Mg isotopic ratios remains debated. δ⁵⁶Fe between the main peridotite minerals shows positive linear correlations with slopes within error of unity, strongly suggesting intra-sample mineral-mineral Fe and Mg isotopic equilibrium. Because inter-mineral isotopic equilibrium should be reached earlier than major element equilibrium via chemical diffusion at mantle temperatures, Fe and Mg isotope ratios of coexisting minerals could be useful tools for justifying mineral thermometry and barometry on the basis of chemical equilibrium between minerals. Although most peridotites in this study exhibit a narrow range in δ⁵⁶Fe, the larger deviations from average δ⁵⁶Fe for three samples likely indicate changes due to metasomatic processes. Two samples show heavy δ⁵⁶Fe relative to the average and they also have high La/Yb and total Fe content, consistent with metasomatic reaction between peridotite and Fe-rich and isotopically heavy melt. The other sample has light δ⁵⁶Fe and slightly heavy δ²⁶Mg, which may reflect Fe-Mg inter-diffusion between peridotite and percolating melt.     

O and H diffusion in uraninite: Implications for fluid-uraninite interactions, nuclear waste disposal, and nuclear forensics

Diffusion coefficients for oxygen and hydrogen were determined from a series of natural uraninite-H₂O experiments between 50 and 700 °C. Under hydrous conditions there are two diffusion mechanisms: (1) an initial extremely fast-path diffusion mechanism that overprinted the oxygen isotopic composition of the entire crystals regardless of temperature and (2) a slower volume-diffusive mechanism dominated by defect clusters that displace or eject nearest neighbor oxygen atoms to form two interstitial sites and two partial vacancies, and by vacancy migration. Using the volume diffusion coefficients in the temperature range of 400-600 °C, diffusion coefficients for oxygen can be represented by D = 1.90e⁻⁵ exp (−123,382 J/RT) cm²/s and for temperatures between 100 and 300 °C the diffusion coefficients can be represented by D = 1.95e⁻¹⁰ exp (−62484 J/RT) cm²/s, where the activation energies for uraninite are 123.4 and 62.5 kJ/mol, respectively. Hydrogen diffusion in uraninite appears to be controlled by similar mechanisms as oxygen. Using the volume diffusion coefficients for temperatures between 50 and 700 °C, diffusion coefficients for hydrogen can be represented by D = 9.28e⁻⁶ exp (−156,528 J/RT) cm²/s for temperatures between 450 and 700 °C and D = 1.39e⁻¹⁴ exp (−34518 J/RT) cm²/s for temperatures between 50 and 400 °C, where the activation energies for uraninite are 156.5 and 34.5 kJ/mol, respectively.Results from these new experiments have implications for isotopic exchange during natural UO₂-water interactions. The exceptionally low δ¹⁸O values of natural uraninites (i.e. −32‰ to −19.5‰) from unconformity-type uranium deposits in Saskatchewan, in conjunction with theoretical and experimental uraninite-water and UO₃-water fractionation factors, suggest that primary uranium mineralization is not in oxygen isotopic equilibrium with coeval clay and silicate minerals. The low δ¹⁸O values have been interpreted as resulting from the low temperature overprinting of primary uranium mineralization in the presence of relatively modern meteoric fluids having δ¹⁸O values of ca. −18‰, despite petrographic and U-Pb isotope data that indicate limited alteration. Our data show that the anomalously low oxygen isotopic composition of the uraninite from the Athabasca Basin can be due to meteoric water overprinting under reducing conditions, and meteoric water or groundwater can significantly affect the oxygen isotopic composition of spent nuclear fuel in a geologic repository, with minimal change to the chemical composition or texture. Moreover, the rather fast oxygen and hydrogen diffusion coefficients for uraninite, especially at low temperatures, suggest that oxygen and hydrogen diffusion may impart characteristic isotopic signals that can be used to track the route of fissile material.     

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.