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

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.     

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

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.     

Assessment of grain-scale homogeneity and equilibration of carbon and oxygen isotope compositions of minerals in carbonate-bearing metamorphic rocks by ion microprobe

Nineteen samples of metamorphosed carbonate-bearing rocks were analyzed for carbon and oxygen isotope ratios by ion microprobe with a ∼5-15 μm spot, three from a regional terrain and 16 from five different contact aureoles. Contact metamorphic rocks further represent four groups: calc-silicate marble and hornfels (6), brucite marble (2), samples that contain a reaction front (4), and samples with a pervasive distribution of reactants and products of a decarbonation reaction (4). The average spot-to-spot reproducibility of standard calcite analyses is ±0.37‰ (2 standard deviations, SD) for δ¹⁸O and ±0.71‰ for δ¹³C. Ten or more measurements of a mineral in a sample that has uniform isotope composition within error of measurement can routinely return a weighted mean with a 95% confidence interval of 0.09-0.16‰ for δ¹⁸O and 0.10-0.29‰ for δ¹³C. Using a difference of >6SD as the criterion, only four of 19 analyzed samples exhibit significant intracrystalline and/or intercrystalline inhomogeneity in δ¹³C at the 100-500 μm scale, with differences within individual grains up to 3.7‰. Measurements are consistent with carbon isotope exchange equilibrium between calcite and dolomite in five of six analyzed samples at the same scale. Because of relatively slow carbon isotope diffusion in calcite and dolomite, differences in δ¹³C can survive intracrystalline homogenization by diffusion during cooling after peak metamorphism and likely represent the effects of prograde decarbonation and infiltration. All but 2 of 11 analyzed samples exhibit intracrystalline differences in δ¹⁸O (up to 9.4‰), intercrystalline inhomogeneity in δ¹⁸O (up to 12.5‰), and/or disequilibrium oxygen isotope fractionations among calcite-dolomite, calcite-quartz, and calcite-forsterite pairs at the 100-500 μm scale. Inhomogeneities in δ¹⁸O and δ¹³C are poorly correlated with only a single mineral (dolomite) in a single sample exhibiting both. Because of relatively rapid oxygen isotope diffusion in calcite, intracrystalline inhomogeneities in δ¹⁸O likely represent partial equilibration between calcite and fluid during retrograde metamorphism. Calcite is in oxygen isotope exchange equilibrium with forsterite in one of four analyzed samples, in equilibrium with dolomite in none of six analyzed samples, and in equilibrium with quartz in neither of two analyzed samples. There are no samples of contact metamorphic rock with analyzed reactants and products of an arrested metamorphic reaction that are in oxygen isotope equilibrium with each other. The degree of departure from equilibrium in analyzed samples is variable and is often related, at least in part, to alteration of δ¹⁸O of calcite during retrograde fluid-rock reaction. In situ sub-grain-scale carbon and oxygen isotope analyses of minerals are advisable in the common applications of stable isotope geochemistry to metamorphic petrology. Correlation of sub-mm scale stable isotope data with imaging will lead to improved understanding of reaction kinetics, reactive fluid flow, and thermal histories during metamorphism.     

Calibration of the carbonate ‘clumped isotope’ paleothermometer for otoliths

Paleothermometry is an essential tool for understanding past changes in climate. The ‘carbonate clumped isotope thermometer’ is a temperature proxy related to ordering of ¹³C and ¹⁸O in the carbonate lattice (based on measurements of ¹³C¹⁸O¹⁶O in CO₂ produced by acid digestion of carbonate). This thermometer has been previously calibrated for inorganic calcite and aragonitic corals [Ghosh P., Adkins J., Affek H., Balta B., Guo W. F., Schauble E. A., Schrag D., and Eiler J. M. (2006) C-13-O-18 bonds in carbonate minerals: a new kind of paleothermometer. Geochim. Cosmochim. Acta70 (6), 1439-1456]. Here we determine the relationship between growth temperatures of aragonitic fish otoliths and abundances of ¹³C¹⁸O¹⁶O produced by acid digestion of those otoliths. Our calibration is based on analyses of otoliths from six species from four genera of modern fish sampled from a latitudinal transect of the Atlantic Ocean between 54° S and 65° N, plus one species from the tropical western Pacific. The temperatures at which fish otoliths precipitated were estimated by the mean temperature in the waters in which they lived, averaged over their estimated lifetimes. Estimated growth temperatures of our samples vary between 2 and 25 °C. Our results show that the abundance of ¹³C¹⁸O¹⁶O in CO₂ produced by acid digestion of fish otolith aragonite is a function of growth temperature, following the relationship: , 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 closely approaches that previously documented for inorganic calcite and aragonitic coral (Ghosh et al., 2006).     

C-O isotope signatures and ‘clumped isotope’ thermometry in foraminifera and coccoliths

Accurate constraints on past ocean temperatures and compositions are critical for documenting climate change and resolving its causes. Most proxies for temperature are not thermodynamically based, appear to be subject to biological processes, require regional calibrations, and/or are influenced by fluid composition. As a result, their interpretation becomes uncertain when they are applied in settings not necessarily resembling those in which they were empirically calibrated. Independent proxies for past temperature could provide an important means of testing and/or expanding on existing reconstructions. Here we report measurements of abundances of stable isotopologues of calcitic and aragonitic benthic and planktic foraminifera and coccoliths, relate those abundances to independently estimated growth temperatures, and discuss the possible scope of equilibrium and kinetic isotope effects. The proportions of ¹³C-¹⁸O bonds in these samples exhibits a temperature dependence that is generally similar to that previously been reported for inorganic calcite and other biologically precipitated carbonate-containing minerals (apatite from fish, reptile, and mammal teeth; calcitic brachiopods and molluscs; aragonitic coral and mollusks). Most species that exhibit non-equilibrium ¹⁸O/¹⁶O (δ¹⁸O) and ¹³C/¹²C (δ¹³C) ratios are characterized by ¹³C-¹⁸O bond abundances that are similar to inorganic calcite and are generally indistinguishable from apparent equilibrium, with possible exceptions among benthic foraminiferal samples from the Arctic Ocean where temperatures are near-freezing. Observed isotope ratios in biogenic carbonates can be explained if carbonate minerals generally preserve a state of ordering that reflects the extent of isotopic equilibration of the dissolved inorganic carbon species.     

Effects of speciation on equilibrium fractionations and rates of oxygen isotope exchange between (PO4)aq and H2O

The effects of phosphate speciation on both rates of isotopic exchange and oxygen isotope equilibrium fractionation factors between aqueous phosphate and water were examined over the temperature range 70 to 180°C. Exchange between phosphate and water is much faster at low pH than at high pH, an observation that is similar to what has been observed in the analogous sulfate-water system. Oxygen isotope fractionations between protonated species like H₃PO₄ and H₂PO₄⁻ that are dominant at relatively low pH and species like PO₄³⁻ and ion pairs like KHPO₄⁻ that are dominant at relatively high pH, range between 5 and 8‰ at the temperatures of the experiments. In aqueous phosphate systems at equilibrium, ¹⁸O/¹⁶O ratios increase with increasing degree of protonation of phosphate. This effect can be explained in part by the relative magnitudes of the dissociation constants of the protonated species. Under equilibrium conditions, carbonate in solution or in solid phases concentrates ¹⁸O relative to orthophosphate in solution or in solid phases at all temperatures, supporting the traditional view that biogenic phosphate is precipitated in near oxygen isotope equilibrium with body/ambient aqueous fluids with no attendant vital effects.     

Oxygen- and carbon-isotopic composition of gypsum-calcite-dolomite crystals and metamorphic marble assemblages

Gypsum and celestine crystals coexisting with those of calcite and dolomite in marble geodes from the Apuan Alps (Italy) display oxygen-isotopic ratios of between +25.8 and +27.8‰ relative to SMOW, with an average value of +27.2‰.The δ ¹⁸O and δ ¹³C values of calcite-dolomite-marble assemblages fall in the range of marine limestone, the dolomites being enriched by about 1‰ both for oxygen and carbon isotopes.The homogeneity of gypsum and celestine δ ¹⁸O values seems to reflect an isotopic equilibrium of sulfate ions with environmental water, prior to their precipitation and at a temperature higher than 40°C.The calcite-dolomite pairs are not in isotopic equilibrium, and their oxygen-isotopic fractionations cannot be used as indicators of their deposition temperatures.     

Reformulated O correction of mass spectrometric stable isotope measurements in carbon dioxide and a critical appraisal of historic ‘absolute’ carbon and oxygen isotope ratios

Mass-spectrometric stable isotope measurements of CO₂ use molecular ion currents at mass-to-charge ratios m/z 44, 45 and 46 to derive the elemental isotope ratios n(¹³C)/n(¹²C) and n(¹⁸O)/n(¹⁶O), abbreviated ¹³C/¹²C and ¹⁸O/¹⁶O, relative to a reference. The ion currents have to be corrected for the contribution of ¹⁷O-bearing isotopologues, the so-called ‘¹⁷O correction’. The magnitude of this correction depends on the calibrated isotope ratios of the reference. Isotope ratio calibrations are difficult and are therefore a matter of debate. Here, I provide a comprehensive evaluation of the existing ¹³C/¹²C (¹³R), ¹⁷O/¹⁶O (¹⁷R) and ¹⁸O/¹⁶O (¹⁸R) calibrations of the reference material Vienna Standard Mean Ocean Water (VSMOW) and CO₂ generated from the reference material Vienna Pee Dee Belemnite (VPDB) by reaction with 100% H₃PO₄ at 25 °C (VPDB-CO₂). I find , ¹⁸R_VSMOW/10⁻⁶ = 2005.20 ± 0.45, ¹³R_VPDB-CO2/10^-6= 11124 ± 45, and ¹⁸R_VPDB-CO2/10^-6=2088.37±0.90. I also rephrase the calculation scheme for the ¹⁷O correction completely in terms of relative isotope ratio differences (δ values). This reveals that only ratios of isotope ratios (namely, ¹⁷R/¹³R and ¹³R¹⁷R/¹⁸R) are required for the ¹⁷O correction. These can be, and have been, measured on conventional stable isotope mass spectrometers. I then show that the remaining error for these ratios of isotope ratios can lead to significant uncertainty in the derived relative ¹³C/¹²C difference, but not for¹⁸O/¹⁶O. Even though inter-laboratory differences can be corrected for by a common ‘ratio assumption set’ and/or normalisation, the ultimate accuracy of the ¹⁷O correction is hereby limited. Errors of similar magnitude can be introduced by the assumed mass-dependent relationship between ¹⁷O/¹⁶O and ¹⁸O/¹⁶O isotope ratios. For highest accuracy in the ¹³C/¹²C ratio, independent triple oxygen isotope measurements are required. Finally, I propose an experiment that allows direct measurement of ¹³R¹⁷R/¹⁸R.     

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.     

Isotopic and chemical intra-shell variations in low-Mg calcite of rudist bivalves (Mollusca-Hippuritacea): disequilibrium fractionations and late Cretaceous seasonality

 Isotopic (δ¹³C, δ¹⁸O) and elemental (Mg, Sr, Mn, Fe) compositions were analysed in sclerochronological profiles of several shells of late Cretaceous rudist bivalves from Greece, Turkey, Somalia and the Arabian Peninsula. The preservation of original compositions of low-Mg calcite of outer shell layers is indicated by constant and high Sr, generally low Fe and Mn, and the preservation of fibrous-prismatic ultrastructures. Cyclic variations in δ¹⁸O and Mg are interpreted to reflect seasonal temperature/salinity cycles and, thus, annual growth increments. In shells of Torreites, amplitudes of correlated δ¹³C and δ¹⁸O cycles cannot be related to reasonable palaeotemperatures or salinity. This isotopic pattern reflects vital fractionations of an extent which is unknown from modern bivalves. In contrast, almost identical ranges and amplitudes of δ¹⁸O cycles are observed in 13 shells of five species from Santonian-Campanian localities in central Greece and northern Turkey, suggesting that seasonal variations in environmental conditions were recorded without significant vital fractionations. The effect of seasonal salinity changes on δ¹⁸O of the shells is evaluated, and mean palaeotemperatures are constrained within the range of 30–32.5  °C. The annual range of temperature was estimated to be 7  °C, assuming a constant salinity. This agrees with other isotopic proxies of Late Cretaceous palaeotemperatures, and with global circulation models which predict higher low-latitude sea-surface temperatures than the present ones. Received: 12 February 1998 / Accepted: 24 May 1999     

Equilibrium thermodynamics of multiply substituted isotopologues of molecular gases

Isotopologues of molecular gases containing more than one rare isotope (multiply substituted isotopologues) can be analyzed with high precision (1σ <0.1‰), despite their low natural abundances (∼ ppm to ppt in air), and can constrain geochemical budgets of natural systems. We derive a method for calculating abundances of all such species in a thermodynamically equilibrated population of isotopologues, and present results of these calculations for O₂, CO, N₂, NO, CO₂, and N₂O between 1000 and 193 to 77 K. In most cases, multiply substituted isotopologues are predicted to be enriched relative to stochastic (random) distributions by ca. 1 to 2‰ at earth-surface temperatures. This deviation, defined as Δ_i for isotopologue i, generally increases linearly with 1/T at temperatures ≤ 500 K. An exception is N₂O, which shows complex temperature dependences and 10’s of per-mill enrichments or depletions of abundances for some isotopologues. These calculations provide a basis for discriminating between fractionations controlled by equilibrium thermodynamics and other sorts of isotopic fractionations in the budgets of atmospheric gases. Moreover, because abundances of multiply substituted isotopologues in thermodynamically equilibrated populations of molecules vary systematically with temperature, they can be used as geothermometers. Such thermometers are unusual in that they involve homogeneous rather than heterogeneous equilibria (e.g., isotopic distribution in gaseous CO₂ alone, rather than difference in isotopic composition between CO₂ and coexisting water). Also, multiple independent thermometers exist for all molecules having more than one multiply substituted isotopologue (e.g., thermometers based on abundances of ¹⁸O¹³C¹⁶O and ¹⁸O¹²C¹⁸O are independent); thus, temperatures estimated by this method can be tested for internal consistency.     

Environmental setting, (micro)morphologies and stable C–O isotope composition of cold climate carbonate precipitates—a review and evaluation of their potential as paleoclimatic proxies

Given the growing interest in carbonate deposits from polar regions as paleoclimatic proxies, this review paper first provides a classification of the various types of cold-climate carbonate precipitates followed by a summary of the ¹³C and ¹⁸O composition of the carbonate deposits and parent water from which the carbonates precipitated. The cold-climate carbonate precipitates were classified into three broad categories: powders, crusts and speleothem. The carbonate powders include those that precipitated in relation to aufeis aggradation (cryogenic aufeis calcite) and in relation to the growth of various annual/perennial ice formations in freezing caves (cryptocrystalline calcite and calcite pearls). The carbonate crusts can be further subdivided based on their lithic environment; those that precipitated on the upper surface of bedrock/clasts (i.e. subglacially precipitated calcite and evaporative calcite crusts); those that are located on the underside of clasts (i.e. pedogenic carbonates); and those that precipitated in rock outcrop fissures (i.e. endostromatolites). The cold-climate carbonate precipitates have a highly variable isotopic composition with δ¹⁸O values ranging between −6.5‰ and 28‰ VSMOW and δ¹³C values in the −10–20‰ VPDB range. However, each type of carbonate precipitates has a specific δ¹³C and δ¹⁸O range, suggesting that their environmental setting and the mechanism by which they formed controls their ¹³C and ¹⁸O signature. It was found that carbonate deposits that precipitated under equilibrium physico-chemical conditions had a δ¹³C value that is in equilibrium with that of the parent water, while its δ¹⁸O composition was more variable, as it is in part controlled by the temperature of reaction and by the δ¹⁸O and calcite saturation state of the parent water. By contrast, the δ¹⁸O composition of biologically precipitated carbonate deposits (endostromatolites) reflect that of the parent water, while its δ¹³C composition was enriched over that of the parent water due to bacterial methanogenesis. In the case of kinetically precipitated carbonate deposits, the δ¹⁸O and δ¹³C values are out-of-equilibrium relative to that of the parent water due to the faster rate of reaction.     

Isotopic partitioning between scallop shell calcite and seawater: effect of shell growth rate

The relationship between molluscan shell growth rate and skeletal δ¹⁸O and δ¹³C was investigated in a detailed field study for the scallop, Pecten maximus. Seasonal variation in shell growth rate was found to be a governing factor influencing shell δ¹⁸O and δ¹³C. At low shell growth rates, shell δ¹⁸O were more positive (of the order +0.4‰) and δ¹³C more negative (up to −2‰) as compared with predicted values for precipitation of inorganic calcite in isotopic equilibrium with seawater. The deviations in δ¹⁸O were hypothesized as reflecting possible differences in solution carbonate chemistry at the site of mineralization in the extrapallial fluid as compared with that of the external seawater medium. The deviations in shell δ¹³C were consistent with incorporation of isotopically depleted respiratory ¹³C (i.e., a metabolic effect). A trend toward more depleted shell δ¹⁸O and δ¹³C values occurred at higher shell growth rates, with negative δ¹⁸O values as compared with predicted equilibrium at shell growth rates above 0.13 mm per day. These simultaneous negative deviations in skeletal δ¹⁸O and δ¹³C were interpreted as resulting from a kinetic effect. The implications for environmental reconstruction from molluscan isotopic records are discussed in light of a model of isotopic behavior based on the findings of the study.     

Calibration of oxygen isotope fractionation and calcite-corundum thermometry in emery at Naxos,Greece

Corundum (Crn), including sapphire, occurs in emery pods surrounded by marble on the island of Naxos, Greece. The emery formed from bauxite deposited in karst that was metamorphosed to 400–700°C at 20–15 Ma. Many of these rocks initially appeared well suited for refractory accessory mineral (RAM) thermometry, which uses oxygen isotope fractionation between a RAM – corundum – and a modally dominant phase with faster diffusion of oxygen – calcite (Cc) – to determine peak metamorphic temperatures. However, previous attempts at oxygen isotope thermometry were confounded by highly variable fractionations (Δ¹⁸O) measured at mm-scale and the uncertain calibration of Δ¹⁸O(Cc-Crn) versus temperature. Secondary ion mass spectrometry (SIMS) permits in situ analysis of δ¹⁸O in corundum and calcite at the 10-μm scale in adjacent grains where textures suggest peak metamorphic equilibrium was attained. SIMS analyses of adjacent mineral pairs in eight rocks yield values of Δ(Cc-Crn) that systematically decrease from 7.2 to 2.9‰ at higher metamorphic grade. Pairing these data with independent temperature estimates from mineral isograds yields an empirical calibration of 1,000 lnα(Cc-Crn) = 2.72 ± 0.3 × 10⁶/T² (T in K). The new fractionations (2.7‰ at 1,000 K) are significantly smaller than those calculated from the modified increment method (6.5‰ at 1,000 K; Zheng, Geochimica et Cosmochimica Acta, 1991, 55:2299–2307; Zheng, Mineral Mag, 1994, 58A:1000–1001), which yield unreasonably high temperatures of 630 to 1,140°C when applied to the new Naxos data. The new calibration of Δ(Cc-Crn) can be combined with published fractionations to calculate A-factors for corundum versus a range of 14 other minerals. These new fractionation factors can be used for thermometry or to constrain the genesis of corundum. A compilation of gem corundum δ¹⁸O values shows that many igneous sapphires, including important deposits of basalt-associated sapphire, are mildly elevated in δ¹⁸O relative to the calculated range in equilibrium with mantle values (4.4–5.7‰) and formed from evolved magmas.     

Stable isotope fractionation between mollusc shells and marine waters from Martinique Island

Forty-nine aragonitic and calcitic shells from 14 species of marine tropical molluscs (Bivalvia, Gastropoda, Polyplacophora) and ambient waters from Martinique have been analyzed for their carbon and oxygen isotope compositions. Mineralogy of shells was systematically determined by Raman spectroscopy that reveals composite shell structures and early processes of diagenetic alteration. In mangrove, brackish waters result from the mixing between 89±1% of seawater and 11±1% of freshwater, a hydrological budget quantified by both oxygen isotope and salinity mass balance calculations. Mollusc shells from the mangrove environment (S=31‰; δ¹⁸O=0.5‰) are characterized by mean δ¹³C values (−1.2‰) lower than those (+2.6‰) living in the open sea (S=35‰; δ¹⁸O=1‰). These low carbon isotope compositions result from the oxidation of organic matter into bicarbonate ions used in the building of mollusc shells. The oxygen isotope compositions of the studied mollusc species are mainly controlled by the temperature and composition of seawater whereas the role of the so-called “vital effects” is negligible. Contrasting with carbon isotopes, variability in the δ¹⁸O values among and within species of mollusc shells is very low (1σ=0.15) for a given littoral environment. Using ambient temperatures of seawater (28-30 °C), oxygen isotope fractionations between all studied living species and environmental waters match those extrapolated from the fractionation equation established for molluscs by Grossman and Ku [Chem. Geol., Isot. Geosci. Sect. 59 (1986) 59] in the range 3-20 °C. By analyzing calcite and aragonite layers from the same shell or by comparing shells from different species living in the same environment, there is no evidence that oxygen isotope fractionation between aragonite and water differs from that between calcite and water. On the basis of these results, we conclude that the oxygen isotope compositions of shells from most fossil mollusc species are suitable to estimate past seawater temperatures at any paleolatitude.     

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.     

Oxygen isotope study of metamorphic and granitic rocks of the Yanai district in the Ryoke belt,Japan

¹⁸O/¹⁶O ratios have been obtained for 134 whole-rocks and minerals from metamorphic and granitic rocks of the Yanai district in the Ryoke belt, Southwest Japan. The ¹⁸O/¹⁶O ratios of pelitic rocks of the marginal metamorphic zone decrease progressively with increasing metamorphic grade. In the gneiss-granite complex (zone of migmatite [1]), the most characteristic feature of the rocks is that oxygen isotopic homogenization proceeds on both local and regional scales in parallel with “granitization” or chemical homogenization. Granitic rocks of various origin are fairly uniform in isotopic composition with δ ¹⁸O of quartz of 12 to 14‰ (SMOW) and δ ¹⁸O of biotite of 7 to 9‰ and are about 3 to 4‰ enriched in ¹⁸O compared to other Cretaceous granites of non-metamorphic terranes in Japan. The high ¹⁸O/¹⁶O ratios of granitic rocks of this district were discussed in relation to the ¹⁸O-depletion in metasediments. Oxygen isotopic fractionations among coexisting minerals from various rock-types of the gneiss-granite complex indicate that these minerals were formed under near isotopic equilibrium at a temperature of about 600 to 700° C. Some abnormal fractionations of quartz-biotite pairs also were obtained for rocks which had undergone a progressive ¹⁸O-depletion or ¹⁸O-enrichment. This is due to high resistivity of quartz and contrastive susceptibility of biotite to isotopic exchange during metamorphism and “granitization”.     

Diagenesis and geochemistry of porites corals from Papua New Guinea: Implications for paleoclimate reconstruction

Coral proxy records of sea surface temperature (SST) and hydrological balance have become important tools in the field of tropical paleoclimatology. However, coral aragonite is subject to post-depositional diagenetic alteration in both the marine and vadose environments. To understand the impact of diagenesis on coral climate proxies, two mid-Holocene Porites corals from raised reefs on Muschu Island, Papua New Guinea, were analysed for Sr/Ca, δ¹⁸O, and δ¹³C along transects from 100% aragonite to 100% calcite. Thin-section analysis showed a characteristic vadose zone diagenetic sequence, beginning with leaching of primary aragonite and fine calcite overgrowths, transitional to calcite void filling and neomorphic, fabric selective replacement of the coral skeleton. Average calcite Sr/Ca and δ¹⁸O values were lower than those for coral aragonite, decreasing from 0.0088 to 0.0021 and −5.2 to −8.1‰, respectively. The relatively low Sr/Ca of the secondary calcite reflects the Sr/Ca of dissolving phases and the large difference between aragonite and calcite Sr/Ca partition coefficients. The decrease in δ¹⁸O of calcite relative to coral aragonite is a function of the δ¹⁸O of precipitation. Carbon-isotope ratios in secondary calcite are variable, though generally lower relative to aragonite, ranging from −2.5 to −10.4%. The variability of δ¹³C in secondary calcite reflects the amount of soil CO₂ contributing ¹³C-depleted carbon to the precipitating fluids. Diagenesis has a greater impact on Sr/Ca than on δ¹⁸O; the calcite compositions reported here convert to SST anomalies of 115°C and 14°C, respectively. Based on calcite Sr/Ca compositions in this study and in the literature, the sensitivity of coral Sr/Ca-SST to vadose-zone calcite diagenesis is 1.1 to 1.5°C per percent calcite. In contrast, the rate of change in coral δ¹⁸O-SST is relatively small (−0.2 to 0.2°C per percent calcite). We show that large shifts in δ¹⁸O, reported for mid-Holocene and Last Interglacial corals with warmer than present Sr/Ca-SSTs, cannot be caused by calcite diagenesis. Low-level calcite diagenesis can be detected through X-ray diffraction techniques, thin section analysis, and high spatial resolution sampling of the coral skeleton and thus should not impede the production of accurate coral paleoclimate reconstructions.