Modelling fractionation of stable isotopes in stalagmites

High resolution δ¹³C and δ¹⁸O profiles recorded in precisely dated speleothems are widely used proxies for the climate of the past. Both δ¹³C and δ¹⁸O depend on several climate related effects including meteorological processes, processes occurring in the soil zone above the cave and isotope fractionation processes occurring in the solution layer on the stalagmite surface. Here we model the latter using a stalagmite isotope and growth model and determine the relationship between the stable isotope values in speleothem calcite and cave parameters, such as temperature, drip interval, water p_CO2 and a mixing coefficient describing mixing processes between the solution layer and the impinging drop.The evolution of δ¹³C values is modelled as a Rayleigh distillation process and shows a pronounced dependence on the residence time of the solution on the stalagmite surface and the drip interval, respectively. The evolution of δ¹⁸O values, in contrast, is also influenced by buffering reactions between the bicarbonate in the solution and the drip water driving the δ¹⁸O value of the bicarbonate towards the value expected for equilibrium isotope fractionation between drip water and calcite. This attenuates the dependence of the δ¹⁸O values on drip interval. The temperature dependence of δ¹⁸O, however, is more pronounced than for δ¹³C and in a similar range as expected for fractionation under equilibrium conditions.We also investigate the isotopic enrichment of the δ¹³C and δ¹⁸O values along individual growth layers and, thus, the slopes expected for Hendy tests. The results show that a positive Hendy test is only possible if isotope fractionation occurred under disequilibrium conditions. However, a negative Hendy test does not exclude that isotope fractionation occurred under disequilibrium conditions. A more reliable indicator for disequilibrium fractionation is the enrichment of the δ¹³C values along an individual growth layer.     

Interpreting the Ca isotope record of marine biogenic carbonates

An 18 million year record of the Ca isotopic composition (δ^44/42Ca) of planktonic foraminiferans from ODP site 925, in the Atlantic, on the Ceara Rise, provides the opportunity for critical analysis of Ca isotope-based reconstructions of the Ca cycle. δ^44/42Ca in this record averages +0.37 ± 0.05 (1σ SD) and ranges from +0.21‰ to +0.52‰. The record is a good match to previously published Neogene Ca isotope records based on foraminiferans, but is not similar to the record based on bulk carbonates, which has values that are as much as 0.25‰ lower. Bulk carbonate and planktonic foraminiferans from core tops differ slightly in their δ^44/42Ca (i.e., by 0.06 ± 0.06‰ (n = 5)), while the difference between bulk carbonate and foraminiferan values further back in time is markedly larger, leaving open the question of the cause of the difference. Modeling the global Ca cycle from downcore variations in δ^44/42Ca by assuming fixed values for the isotopic composition of weathering inputs (δ^44/42Ca_w) and for isotope fractionation associated with the production of carbonate sediments (Δ_sed) results in unrealistically large variations in the total mass of Ca²⁺ in the oceans over the Neogene. Alternatively, variations of ±0.05‰ in the Ca isotope composition of weathering inputs or in the extent of fractionation of Ca isotopes during calcareous sediment formation could entirely account for variations in the Ca isotopic composition of marine carbonates. Ca isotope fractionation during continental weathering, such as has been recently observed, could easily result in variations in δ^44/42Ca_w of a few tenths of permil. Likewise a difference in the fractionation factors associated with aragonite versus calcite formation could drive shifts in Δ_sed of tenths of permil with shifts in the relative output of calcite and aragonite from the ocean. Until better constraints on variations in δ^44/42Ca_w and Δ_sed have been established, modeling the Ca²⁺ content of seawater from Ca isotope curves should be approached cautiously.     

Calcium isotopes in a proglacial weathering environment: Damma glacier, Switzerland

The biogeochemical cycling and isotopic fractionation of calcium during the initial stages of weathering were investigated in an alpine soil chronosequence (Damma glacier, Switzerland). This site has a homogeneous silicate lithology and minimal biological impacts due to sparse vegetation cover. Calcium isotopic compositions, obtained by TIMS using a ⁴³Ca-⁴⁶Ca double spike, were measured in the main Ca pools. During this very early stage of weathering, the young soils which have formed (δ^44/42Ca=+0.44‰) were indistinguishable to the rocks from which they were derived (δ^44/42Ca=+0.44‰) and stream water (δ^44/42Ca=+0.48‰) was also within error of the average rock. This lack of variation indicates that the dissolution of the bulk silicate rock does not strongly fractionate Ca isotopes. The only Ca pool which was strongly fractionated from bulk rock was vegetation, which exhibited an enrichment of light Ca isotopes. Significant Ca isotope fractionation between bulk rock and the dissolved flux of Ca is likely to only occur where the Ca biogeochemical cycle is dominated by secondary processes such as biological cycling, adsorption and secondary mineral precipitation.     

Model for kinetic effects on calcium isotope fractionation (δCa) in inorganic aragonite and cultured planktonic foraminifera

The calcium isotope ratios (δ⁴⁴Ca = [(⁴⁴Ca/⁴⁰Ca)_sample/(⁴⁴Ca/⁴⁰Ca)_standard −1] · 1000) of Orbulina universa and of inorganically precipitated aragonite are positively correlated to temperature. The slopes of 0.019 and 0.015‰ °C⁻¹, respectively, are a factor of 13 and 16 times smaller than the previously determined fractionation from a second foraminifera, Globigerinoides sacculifer, having a slope of about 0.24‰ °C⁻¹. The observation that δ⁴⁴Ca is positively correlated to temperature is opposite in sign to the oxygen isotopic fractionation (δ¹⁸O) in calcium carbonate (CaCO₃). These observations are explained by a model which considers that Ca²⁺-ions forming ionic bonds are affected by kinetic fractionation only, whereas covalently bound atoms like oxygen are affected by kinetic and equilibrium fractionation. From thermodynamic consideration of kinetic isotope fractionation, it can be shown that the slope of the enrichment factor α(T) is mass-dependent. However, for O. universa and the inorganic precipitates, the calculated mass of about 520 ± 60 and 640 ± 70 amu (atomic mass units) is not compatible with the expected ion mass for ⁴⁰Ca and ⁴⁴Ca. To reconcile this discrepancy, we propose that Ca diffusion and δ⁴⁴Ca isotope fractionation at liquid/solid transitions involves Ca²⁺-aquocomplexes (Ca[H₂O]_n²⁺ · mH₂O) rather than pure Ca²⁺-ion diffusion. From our measurements we calculate that such a hypothesized Ca²⁺-aquocomplex correlates to a hydration number of up to 25 water molecules (490 amu). For O. universa we propose that their biologically mediated Ca isotope fractionation resembles fractionation during inorganic precipitation of CaCO₃ in seawater. To explain the different Ca isotope fractionation in O. universa and in G. sacculifer, we suggest that the latter species actively dehydrates the Ca²⁺-aquocomplex before calcification takes place. The very different temperature response of Ca isotopes in the two species suggests that the use of δ⁴⁴Ca as a temperature proxy will require careful study of species effects.     

Ca isotopes in carbonate sediment and pore fluid from ODP Site 807A: The Ca(aq)-calcite equilibrium fractionation factor and calcite recrystallization rates in Pleistocene sediments

The calcium isotopic compositions (δ⁴⁴Ca) of 30 high-purity nannofossil ooze and chalk and 7 pore fluid samples from ODP Site 807A (Ontong Java Plateau) are used in conjunction with numerical models to determine the equilibrium calcium isotope fractionation factor (α_s−f) between calcite and dissolved Ca²⁺ and the rates of post-depositional recrystallization in deep sea carbonate ooze. The value of α_s−f at equilibrium in the marine sedimentary section is 1.0000 ± 0.0001, which is significantly different from the value (0.9987 ± 0.0002) found in laboratory experiments of calcite precipitation and in the formation of biogenic calcite in the surface ocean. We hypothesize that this fractionation factor is relevant to calcite precipitation in any system at equilibrium and that this equilibrium fractionation factor has implications for the mechanisms responsible for Ca isotope fractionation during calcite precipitation. We describe a steady state model that offers a unified framework for explaining Ca isotope fractionation across the observed precipitation rate range of ∼14 orders of magnitude. The model attributes Ca isotope fractionation to the relative balance between the attachment and detachment fluxes at the calcite crystal surface. This model represents our hypothesis for the mechanism responsible for isotope fractionation during calcite precipitation. The Ca isotope data provide evidence that the bulk rate of calcite recrystallization in freshly-deposited carbonate ooze is 30-40%/Myr, and decreases with age to about 2%/Myr in 2-3 million year old sediment. The recrystallization rates determined from Ca isotopes for Pleistocene sediments are higher than those previously inferred from pore fluid Sr concentration and are consistent with rates derived for Late Pleistocene siliciclastic sediments using uranium isotopes. Combining our results for the equilibrium fractionation factor and recrystallization rates, we evaluate the effect of diagenesis on the Ca isotopic composition of marine carbonates at Site 807A. Since calcite precipitation rates in the sedimentary column are many orders of magnitude slower than laboratory experiments and the pore fluids are only slightly oversaturated with respect to calcite, the isotopic composition of diagenetic calcite is likely to reflect equilibrium precipitation. Accordingly, diagenesis produces a maximum shift in δ⁴⁴Ca of +0.15‰ for Site 807A sediments but will have a larger impact where sedimentation rates are low, seawater circulates through the sediment pile, or there are prolonged depositional hiatuses.     

Calcium-isotope fractionation in selected modern and ancient marine carbonates

The calcium-isotope composition (δ^44/42Ca) was analyzed in modern, Cretaceous and Carboniferous marine skeletal carbonates as well as in bioclasts, non-skeletal components, and diagenetic cements of Cretaceous and Carboniferous limestones. In order to gain insight in Ca²⁺_aq-CaCO₃-isotope fractionation mechanisms in marine carbonates, splits of samples were analyzed for Sr, Mg, Fe, and Mn concentrations and for their oxygen and carbon isotopic composition. Biological carbonates generally have lower δ^44/42Ca values than inorganic marine cements, and there appears to be no fractionation between seawater and marine inorganic calcite. A kinetic isotope effect related to precipitation rate is considered to control the overall discrimination against ⁴⁴Ca in biological carbonates when compared to inorganic precipitates. This is supported by a well-defined correlation of the δ^44/42Ca values with Sr concentrations in Cretaceous limestones that contain biological carbonates at various stages of marine diagenetic alteration. No significant temperature dependence of Ca-isotope fractionation was found in shells of Cretaceous rudist bivalves that have recorded large seasonal temperature variations as derived from δ¹⁸O values and Mg concentrations. The reconstruction of secular variations in the δ^44/42Ca value of seawater from well preserved skeletal calcite is compromised by a broad range of variation found in both modern and Cretaceous biological carbonates, independent of chemical composition or mineralogy. Despite these variations that may be due to still unidentified biological fractionation mechanisms, the δ^44/42Ca values of Cretaceous skeletal calcite suggest that the δ^44/42Ca value of Cretaceous seawater was 0.3-0.4‰ lower than that of the modern ocean.     

Non-biological fractionation of stable Ca isotopes in soils of the Atacama Desert, Chile

We measured Ca stable isotope ratios (δ^44/40Ca) in an ancient (2 My), hyperarid soil where the primary source of mobile Ca is atmospheric deposition. Most of the Ca in the upper meter of this soil (3.5 kmol m⁻²) is present as sulfates (2.5 kmol m⁻²), and to a lesser extent carbonates (0.4 kmol m⁻²). In aqueous extracts of variably hydrated calcium sulfate minerals, δ^44/40Ca_E values (vs. bulk Earth) increase with depth (1.4 m) from a minimum of −1.91‰ to a maximum of +0.59‰. The trend in carbonate-δ^44/40Ca in the top six horizons resembles that of sulfate-δ^44/40Ca, but with values 0.1-0.6‰ higher. The range of observed Ca isotope values in this soil is about half that of δ^44/40Ca values observed on Earth. Linear correlation among δ^44/40Ca, δ³⁴S and δ¹⁸O values indicates either (a) a simultaneous change in atmospheric input values for all three elements over time, or (b) isotopic fractionation of all three elements during downward transport. We present evidence that the latter is the primary cause of the isotopic variation that we observe. Sulfate-δ³⁴S values are positively correlated with sulfate-δ¹⁸O values (R² = 0.78) and negatively correlated with sulfate δ^44/40Ca_E values (R² = 0.70). If constant fractionation and conservation of mass with downward transport are assumed, these relationships indicate a δ^44/40Ca fractionation factor of −0.4‰ in CaSO₄. The overall depth trend in Ca isotopes is reproduced by a model of isotopic fractionation during downward Ca transport that considers small and infrequent but regularly recurring rainfall events. Near surface low Ca isotope values are reproduced by a Rayleigh model derived from measured Ca concentrations and the Ca fractionation factor predicted by the relationship with S isotopes. This indicates that the primary mechanism of stable isotope fractionation in CaSO₄ is incremental and effectively irreversible removal of an isotopically enriched dissolved phase by downward transport during small rainfall events.     

Sr/Ca and Ca/Ca fractionation during inorganic calcite formation: II. Ca isotopes

Ca isotope fractionation during inorganic calcite formation was experimentally studied by spontaneous precipitation at various precipitation rates (1.8 < log R < 4.4 μmol/m²/h) and temperatures (5, 25, and 40 °C) with traces of Sr using the CO₂ diffusion technique.Results show that in analogy to Sr/Ca [see Tang J., Köhler S. J. and Dietzel M. (2008) Sr²⁺/Ca²⁺ and ⁴⁴Ca/⁴⁰Ca fractionation during inorganic calcite formation: I. Sr incorporation. Geochim. Cosmochim. Acta] the ⁴⁴Ca/⁴⁰Ca fractionation during calcite formation can be followed by the Surface Entrapment Model (SEMO). According to the SEMO calculations at isotopic equilibrium no fractionation occurs (i.e., the fractionation coefficient α_calcite-aq = (⁴⁴Ca/⁴⁰Ca)_s/(⁴⁴Ca/⁴⁰Ca)_aq = 1 and Δ^44/40Ca_calcite-aq = 0‰), whereas at disequilibrium ⁴⁴Ca is fractionated in a primary surface layer (i.e., the surface entrapment factor of ⁴⁴Ca, F_44Ca < 1). As a crystal grows at disequilibrium, the surface-depleted ⁴⁴Ca is entrapped into the newly formed crystal lattice. ⁴⁴Ca depletion in calcite can be counteracted by ion diffusion within the surface region. Our experimental results show elevated ⁴⁴Ca fractionation in calcite grown at high precipitation rates due to limited time for Ca isotope re-equilibration by ion diffusion. Elevated temperature results in an increase of ⁴⁴Ca ion diffusion and less ⁴⁴Ca fractionation in the surface region. Thus, it is predicted from the SEMO that an increase in temperature results in less ⁴⁴Ca fractionation and the impact of precipitation rate on ⁴⁴Ca fractionation is reduced.A highly significant positive linear relationship between absolute ⁴⁴Ca/⁴⁰Ca fractionation and the apparent Sr distribution coefficient during calcite formation according to the equation

Δ^44/40Ca_calcite-aq=(−1.90±0.26)·logD_Sr−2.83±0.28     

Stable isotope profiles (Ca,O, C) through modern brachiopod shells of T. septentrionalis and G. vitreus: Implications for calcium isotope paleo-ocean chemistry

Brachiopod shells are widely used as an archive to reconstruct elemental and isotopic composition of seawater. Studies, focused on oxygen and carbon isotopes over the last decades, are increasingly extending to the emerging calcium isotope system. To date, only little attention has been paid to test the reliability of fossil brachiopods on their modern counterparts.In this context, the present study investigates two modern brachiopods, Terebratulina septentrionalis (eastern Canada, 5–30 m depth, 7.1 °C seasonal temperature variation, two-layer shell) and Gryphus vitreus (northern Mediterranean, 200 m depth, constant all-year round temperature, three-layer shell). Both species were sampled along the ontogenetic growth direction and calcium, oxygen, and carbon isotopes as well as elemental concentration were measured. Calcium isotopes were analyzed on TIMS. The elemental composition was analyzed by LA-ICP-MS and ICP-AES.The results indicate an intra-specimen δ^44/40Ca variation ranging from 0.16 to 0.33‰, pointing to a fairly homogenous distribution of calcium isotopes in brachiopod shells. However, in the light of the suggested 0.7‰ increase in calcium isotopes over the Phanerozoic such intra-specimen variations constrain ocean reconstruction. δ^44/40Ca values of T. septentrionalis do not seem to be affected by growth rate. Calcium isotopic values of G. vitreus are heavy in the central part of the shell and trend towards lighter values in peripheral areas approaching the maximum isotopic composition of T. septentrionalis. The maximum inter-species δ^44/40Ca difference of 0.62‰ between T. septentrionalis and G. vitreus indicates that care should be taken when using different taxa, species with different strontium content or brachiopods with specialized shell structure, such as G. vitreus, for ocean water reconstruction in terms of Ca isotopic composition. T. septentrionalis may record Ca isotopic fractionation related to seasonal seawater temperature variations in its shell but this is difficult to resolve at the current analytical precision. Average δ¹⁸O-derived temperatures of the two investigated species are close to on-site measured temperatures.     

Climatic and environmental controls on speleothem oxygen-isotope values

Variations in speleothem oxygen-isotope values (δ¹⁸O) result from a complicated interplay of environmental controls and processes in the ocean, atmosphere, soil zone, epikarst, and cave system. As such, the controls on speleothem δ¹⁸O values are extremely complex. An understanding of the processes that control equilibrium and kinetic fractionation of oxygen isotopes in water and carbonate species is essential for the proper interpretation of speleothem δ¹⁸O as paleoclimate and paleoenvironmental proxies, and is best complemented by study of site-specific cave processes such as infiltration, flow routing, drip seasonality and saturation state, and cave microclimate, among others. This review is a process-based summary of the multiple controls on δ¹⁸O in the atmosphere, soil, epikarst, and speleothem calcite, illustrated with case studies. Primary controls of δ¹⁸O in the atmosphere include temperature and relative humidity through their role in the multiple isotope “effects”. Variability and modifications of water δ¹⁸O values in the soil and epikarst zones are dominated by evaporation, mixing, and infiltration of source waters. The isotopically effective recharge into a cave system consists of those waters that participate in precipitation of CaCO₃, resulting in calcite deposition rates which may be biased to time periods with optimal dripwater saturation state. Recent modeling, experimental, and observational data yield insight into the significance of kinetic fractionation between dissolved carbonate phases and solid CaCO₃, and have implications for the ‘Hendy’ test. To assist interpretation of speleothem δ¹⁸O time series, quantitative and semi-quantitative δ¹⁸O-climate calibrations are discussed with an emphasis on some of the difficulties inherent in using modern spatial and temporal isotope gradients to interpret speleothems as paleoclimate proxy records. Finally, several case studies of globally significant speleothem paleoclimate records are discussed that show the utility of δ¹⁸O to reconstruct past climate changes in regions that have been typically poorly represented in paleoclimate records, such as tropical and subtropical terrestrial locations. The new approach to speleothem paleoclimatology emphasizes climate teleconnections between regions and attribution of forcing mechanisms. Such investigations allow paleoclimatologists to infer regional to global-scale climate dynamics.     

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.     

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.     

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.     

Calcium isotope fractionation in calcite and aragonite

Calcium isotope fractionation was measured on skeletal aragonite and calcite from different marine biota and on inorganic calcite. Precipitation temperatures ranged from 0 to 28°C. Calcium isotope fractionation shows a temperature dependence in accordance with previous observations: 1000 · ln(α_cc) = −1.4 + 0.021 · T (°C) for calcite and 1000 · ln(α_ar) = −1.9 + 0.017 · T (°C) for aragonite. Within uncertainty the temperature slopes are identical for the two polymorphs. However, at all temperatures calcium isotopes are more fractionated in aragonite than in calcite. The offset in δ^44/40Ca is about 0.6‰. The underlying mechanism for this offset may be related to the different coordination numbers and bond strengths of the calcium ions in calcite and aragonite crystals, or to different Ca reaction behavior at the solid-liquid interface. Recently, the observed temperature dependence of the Ca isotope fractionation was explained quantitatively by the temperature control on precipitation rates of calcium carbonates in an experimental setting (Lemarchand et al., 2004). We show that this mechanism can in principle also be applied to CaCO₃ precipitation in natural environments in normal marine settings. Following this model, Ca isotope fractionation in marine Ca carbonates is primarily controlled by precipitation rates. On the other hand the larger Ca isotope fractionation of aragonite compared to calcite can not be explained by different precipitation rates. The rate control model of Ca isotope fractionation predicts a strong dependence of the Ca isotopic composition of carbonates on ambient CO₃²⁻ concentration. While this model is in general accordance with our observations in marine carbonates, cultured specimens of the planktic foraminifer Orbulina universa show no dependence of Ca-isotope fractionation on the ambient CO₃²⁻ concentration. The latter observation implies that the carbonate chemistry in the calcifying vesicles of the foraminifer is independent from the ambient carbonate ion concentration of the surrounding water.     

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.     

Calcium isotope ratios in animal and human bone

Calcium isotopes in tissues are thought to be influenced by an individual’s diet, reflecting parameters such as trophic level and dairy consumption, but this has not been carefully assessed. We report the calcium isotope ratios (δ^44/42Ca) of modern and archaeological animal and human bone (n = 216). Modern sheep raised at the same location show 0.14 ± 0.08‰ higher δ^44/42Ca in females than in males, which we attribute to lactation by the ewes. In the archaeological bone samples the calcium isotope ratios of the herbivorous fauna vary by location. At a single site, the archaeological fauna do not show a trophic level effect. Humans have lower δ^44/42Ca than the mean site fauna by 0.22 ± 0.22‰, and the humans have a greater δ^44/42Ca range than the animals. No effect of sex or age on the calcium isotope ratios was found, and intra-individual skeletal δ^44/42Ca variability is negligible. We rule out dairy consumption as the main cause of the lower human δ^44/42Ca, based on results from sites pre-dating animal domestication and dairy availability, and suggest instead that individual physiology and calcium intake may be important in determining bone calcium isotope ratios.     

Controls on calcium isotope fractionation in cultured planktic foraminifera, Globigerinoides ruber and Globigerinella siphonifera

Specimens of two species of planktic foraminifera, Globigerinoides ruber and Globigerinella siphonifera, were grown under controlled laboratory conditions at a range of temperatures (18-31 °C), salinities (32-44 psu) and pH levels (7.9-8.4). The shells were examined for their calcium isotope compositions (δ^44/40Ca) and strontium to calcium ratios (Sr/Ca) using Thermal Ionization Mass Spectrometry and Inductively Coupled Plasma Mass Spectrometry. Although the total variation in δ^44/40Ca (∼0.3‰) in the studied species is on the same order as the external reproducibility, the data set reveals some apparent trends that are controlled by more than one environmental parameter. There is a well-defined inverse linear relationship between δ^44/40Ca and Sr/Ca in all experiments, suggesting similar controls on these proxies in foraminiferal calcite independent of species. Analogous to recent results from inorganically precipitated calcite, we suggest that Ca isotope fractionation and Sr partitioning in planktic foraminifera are mainly controlled by precipitation kinetics. This postulation provides us with a unique tool to calculate precipitation rates and draws support from the observation that Sr/Ca ratios are positively correlated with average growth rates. At 25 °C water temperature, precipitation rates in G. siphonifera and G. ruber are calculated to be on the order of 2000 and 3000 μmol/m²/h, respectively. The lower δ^44/40Ca observed at ?29 °C in both species is consistent with increased precipitation rates at high water temperatures. Salinity response of δ^44/40Ca (and Sr/Ca) in G. siphonifera implies that this species has the highest precipitation rates at the salinity of its natural habitat, whereas increasing salinities appear to trigger higher precipitation rates in G. ruber. Isotope effects that cannot be explained by precipitation rate in planktic foraminifera can be explained by a biological control, related to a vacuolar pathway for supply of ions during biomineralization and a pH regulation mechanism in these vacuoles. In case of an additional pathway via cross-membrane transport, supplying light Ca for calcification, the δ^44/40Ca of the reservoir is constrained as −0.2‰ relative to seawater. Using a Rayleigh distillation model, we calculate that calcification occurs in a semi-open system, where less than half of the Ca supplied by vacuolization is utilized for calcite precipitation. Our findings are relevant for interpreting paleo-proxy data on δ^44/40Ca and Sr/Ca in foraminifera as well as understanding their biomineralization processes.     

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 stalagmite growth and δC as a function of drip interval and temperature

Growth and stable isotope composition of stalagmites are affected by climate changes. To understand the underlying mechanisms, we developed a time dependent multi-box model that describes stalagmite growth and stable isotope fractionation of carbon under disequilibrium conditions. The model takes variations of the drip interval, temperature and the amount of mixing between the impinging drops and the solution on top of the stalagmite into account, which allows to quantify the influence of these parameters. To calculate the variations of δ¹³C, the multi-box model was calibrated by comparison with an existing growth model. The results show that drip interval, temperature and the mixing coefficient do have a significant influence on δ¹³C. However, considering the higher natural variability of the drip interval, this parameter might have the largest influence.     

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.