In order to improve the understanding of variations of speleothem δ18O and δ13C values in the context of palaeoclimate research, it is important to quantify the isotope fractionation processes ...influencing the δ18O and δ13C values of stalagmites. Here we present an extended version of speleothem stable isotope model that accounts for evaporation and condensation effects during precipitation of calcite. The ISOLUTION-model allows to calculate the effect of evaporation on δ18Ocalcite and δ13Ccalcite values in dependence on relative humidity and wind velocity. Our results reveal that evaporation may have a significant effect on δ18Ocalcite and δ13Ccalcite due to the loss of H2O from the solution layer and the related increase of the Ca2+ concentration. This leads to higher precipitation rates and, consequently, larger isotope fractionation effects. The effect on speleothem δ18Ocalcite and δ13Ccalcite can be as high as changes caused by variations of temperature, drip interval, pCO2.drip and pCO2.cave. Furthermore, the change of pCO2.cave due to ventilation effects, may also affect speleothem δ18Ocalcite and δ13Ccalcite values because of the effect of pCO2.cave on the equilibrium concentrations of the CO2–H2O–CaCO3-system and, thus, the calcite precipitation rate. The variability of other cave parameters (temperature, drip interval, pCO2.drip, pCO2.cave) can have a comparable effect on δ18Ocalcite and δ13Ccalcite values depending on the magnitude of variability. We validated the ISOLUTION-model by modelling the δ18Ocalcite values from a cave-analogue experiment and successfully reproduced the values for high temperatures.
Here we present novel cave-analogue experiments directly investigating stable carbon and oxygen isotope fractionation between the major involved species of the carbonate system (HCO3−, CO2, CaCO3 and ...H2O). In these experiments, which were performed under controlled conditions inside a climate box, a thin film of solution flew down an inclined marble or glass plate. After different distances of flow and, thus, residence times on the plate, pH, electrical conductivity, supersaturation with respect to calcite, precipitation rate as well as the δ18O and δ13C values of the dissolved inorganic carbon (DIC) and the precipitated CaCO3 were obtained.
Progressive precipitation of CaCO3 along the plate is accompanied by degassing of CO2 and stable isotope fractionation, and the system is driven out of isotope equilibrium. We observe a strong enrichment of the δ13C values with increasing residence time and a smaller enrichment in δ18O. The temporal evolution of the δ18O and δ13C values of both the DIC and the precipitated CaCO3 can be explained by a Rayleigh fractionation model, but the observed enrichment in δ13C values is much larger than expected based on isotope equilibrium fractionation factors.
Our setup enables to determine the fractionation between CaCO3 and HCO3−, i.e., εCaCO3/HCO3−. Carbon isotope fractionation, 13εCaCO3/HCO3−, is strongly negative for all experiments and much lower than equilibrium isotope fractionation (0–1‰). In addition, 13εCaCO3/HCO3− decreases with increasing residence time on the plate, and thus decreasing supersaturation with respect to calcite. Thus, isotope fractionation depends on precipitation rate and consequently occurs under kinetic conditions. This is in contrast to previous studies, which found no rate-dependence and no or even a positive carbon isotope fractionation between CaCO3 and HCO3−. Oxygen isotope fractionation, 18εCaCO3/HCO3−, is also negative and dependent on precipitation rate. Since no literature values for 18εCaCO3/HCO3− are available, we calculated 18εCaCO3/HCO3− using equilibrium oxygen isotope fractionation factors between water and calcite and water and HCO3−, respectively. At the beginning of the plate, the fractionation is in agreement with the fractionation calculated using fractionation factors determined in cave systems.
The observed fractionation between CaCO3 and water, 1000ln18α, is also in good agreement with the values determined in cave systems and shows a very similar temperature dependence 1000ln18α=16.516±1.267∗103T−26.141±4.356. However, with progressive precipitation of CaCO3 along the plate, the system is forced out of isotope equilibrium with the water, and 1000ln18α increases.
The large, negative, rate-dependent isotope fractionations observed in this study suggest that precipitation of speleothem calcite is strongly kinetically controlled and may, thus, have a large effect on speleothem δ18O and δ13C values. Since these values may erroneously be interpreted as reflecting changes in past temperature, precipitation and/or vegetation density, these results have important implications for paleoclimate reconstructions from speleothems.
Abstract
Surface temperature is a fundamental parameter of Earth’s climate. Its evolution through time is commonly reconstructed using the oxygen isotope and the clumped isotope compositions of ...carbonate archives. However, reaction kinetics involved in the precipitation of carbonates can introduce inaccuracies in the derived temperatures. Here, we show that dual clumped isotope analyses, i.e., simultaneous ∆
47
and ∆
48
measurements on the single carbonate phase, can identify the origin and quantify the extent of these kinetic biases. Our results verify theoretical predictions and evidence that the isotopic disequilibrium commonly observed in speleothems and scleractinian coral skeletons is inherited from the dissolved inorganic carbon pool of their parent solutions. Further, we show that dual clumped isotope thermometry can achieve reliable palaeotemperature reconstructions, devoid of kinetic bias. Analysis of a belemnite rostrum implies that it precipitated near isotopic equilibrium and confirms the warmer-than-present temperatures during the Early Cretaceous at southern high latitudes.
Stable oxygen and carbon isotope ratios (613C and 618O) are the most applied climate and environmental proxies in speleothems allowing to infer past changes in cave drip water 613C and 618O related ...to climate and environmental variations from above the cave. However, disequilibrium isotope fractionation processes can modify 613C and 618O values in speleothems, which is in most cases difficult to estimate due to inter-dependencies on various cave specific parameter. To better understand the effect of these disequilibrium isotope fractionation processes proxy system models were developed in recent years, such as the ISOLUTION model. Here the code of the ISOLUTION model is made available for the public and the speleothem community to be applied to research questions that arise from e.g. monitoring programs that investigate 613C and 618O values of in situ calcite precipitates on watch glasses or modern speleothem calcite, respectively. Another application of the ISOLUTION model is to investigate the dependence of calcite 613C and 618O on the variation of one or multiple cave specific parameter, such as cave air temperature, drip interval, cave air pCO2, Ca2+ concentration of the drip water as well as on relative humidity and wind velocity. This allows to quantitatively estimate the effect of disequilibrium isotope fractionation processes in individual caves and drip sites on speleothem 613C and 618O values for modern and past climates and may help to further elucidate the complex interplay of kinetic and disequilibrium isotope fractionation.
The high spatial resolution of laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) microanalysis of speleothems and biogenic calcium carbonates makes it possible to improve our ...understanding of past climatic conditions. However, there are analytical problems associated with this method, such as molecular interferences, elemental fractionation, and mass-load-dependent matrix effects, when using NIST silicate reference glasses for calibration.
Our study with a sector-field ICP-MS shows that many masses are affected by interferences, such as 24Mg+ by 48Ca++ or 31P+ by 15N16O+. Elemental fractionation and mass-load-dependent matrix effects have been detected for both 213 and 193nm Nd:YAG lasers. They are small for refractory lithophile elements, in particular for the 193nm laser (<5%). For chalcophile/siderophile trace elements with low boiling points, however, these effects are large (up to 40% and 20% for the 213 and 193nm lasers, respectively), and variable for the NIST glasses and carbonates.
In this paper, a protocol for precise and accurate LA-ICP-MS analysis of calcium carbonates is established. Isotopes for interference-free measurements at low (M/ΔM ~300) and medium (M/ΔM ~4000) mass resolution have been identified. The NIST glasses are suitable for calibration of refractory lithophile element concentrations. However, matrix-matched calibration with carbonate reference materials, such as USGS MACS-1 or MACS-3, is necessary for accurate analysis of chalcophile/siderophile elements with low boiling points. Analytical reproducibility (RSD) is a factor of ~2 better using the 193nm laser than the 213nm laser. As shown by the analysis of MACS-1 and MACS-3, LA-ICP-MS analyses yield results that agree with the reference values within relative uncertainties of ca. 5–10% at the 95% confidence level.
We have applied LA-ICP-MS for the determination of trace element concentrations in calcite and aragonite layers of a stalagmite and found large variations for Mg, Zn, Sr, and U. In ostracod shells, the concentrations of some trace elements (e.g., Sr and Ba) vary significantly, indicating the potential for paleoclimate research.
► We identify problems of LA-ICP-MS for microanalysis of calcium carbonates. ► Accurate data are obtained by a careful choice of isotopes and reference materials. ► LA-ICP-MS analyses of carbonates can be used to study past climate conditions.
We compared the suitability of two skeletal materials of the Atlantic brain coral Diploria strigosa for 230Th/U-dating: the commonly used bulk material comprising all skeletal elements and the denser ...theca wall material. Eight fossil corals of presumably Last Interglacial age from Bonaire, southern Caribbean Sea, were investigated, and several sub-samples were dated from each coral. For four corals, both the ages and the activity ratios of the bulk material and theca wall agree within uncertainty. Three corals show significantly older ages for their bulk material than for their theca wall material as well as substantially elevated 232Th content and (230Th/238U) ratios. The bulk material samples of another coral show younger ages and lower (230Th/238U) ratios than the corresponding theca wall samples. This coral also contains a considerable amount of 232Th.
The application of the available open-system models developed to account for post-depositional diagenetic effects in corals shows that none of the models can successfully be applied to the Bonaire corals. The most likely explanation for this observation is that the assumptions of the models are not fulfilled by our data set.
Comparison of the theca wall and bulk material data enables us to obtain information about the open-system processes that affected the corals. The corals showing apparently older ages for their bulk material were probably affected by contamination with a secondary (detrital) phase. The most likely source of the detrital material is carbonate sand. The higher (230Th/232Th) ratio of this material implies that detrital contamination would have a much stronger impact on the ages than a contaminant with a bulk Earth (230Th/232Th) ratio and that the threshold for the commonly applied 232Th reliability criterion would be much lower than the generally used value of 1ngg−1. The coral showing apparently younger ages for its bulk material was probably influenced by more than one diagenetic process. A potential scenario is a combination of detrital contamination and U addition by secondary pore infillings.
Our results show that the dense theca wall material of D. strigosa is generally less affected by post-depositional open-system behaviour and better suited for 230Th/U-dating than the bulk material. This is also obvious from the fact that all ages of theca wall material reflect a Last Interglacial origin (∼125ka), whereas the bulk material samples are either substantially older or younger. However, for some corals, the 230Th/U-ages and activity ratios of the bulk material and the theca wall samples are similar. This shows that strictly reliable 230Th/U-ages can also be obtained from bulk material samples of exceptionally well-preserved corals. However, the bulk material samples more frequently show elevated activity ratios and ages than the corresponding theca wall samples. Our findings should be generally applicable to brain corals (Mussidae) that are found in tropical oceans worldwide and may enable reliable 230Th/U-dating of fossil corals with similar skeletal architecture, even if their bulk skeleton is altered by diagenesis.
The 230Th/U-ages we consider reliable (120–130ka), along with a recently published age of 118ka, provide the first comprehensive dating of the elevated lower reef terrace at Bonaire (118–130ka), which is in agreement in timing and duration with other Last Interglacial records.
In situ laser ablation analyses rely on the microanalytical homogeneity of reference materials (RMs) and a similar matrix and mass fraction between unknown samples and RMs to obtain reliable results. ...Suitable carbonate and phosphate RMs for determination of Sr isotope ratios in such materials are limited. Thus, we determined 87Sr/86Sr ratios of several carbonate (JCt‐1, JCp‐1, MACS‐1, MACS‐3) and phosphate (MAPS‐4, MAPS‐5, NIST SRM 1400, NIST SRM 1486) international RMs using dissolved samples and two different multi‐collector inductively coupled plasma‐mass spectrometers (MC‐ICP‐MS). Our Sr isotope data are in agreement with published data and have an improved measurement precision for some RMs. For MACS‐1, we present the first 87Sr/86Sr value. We tested the suitability of these materials for microanalytical analyses by LA‐MC‐ICP‐MS, with two different laser ablation systems: a conventional nanosecond laser and a state‐of‐the‐art femtosecond laser. We investigated the RMs micro‐homogeneity and compared the data with our solution data. Both laser ablation systems yielded identical 87Sr/86Sr ratios within uncertainty to the solution data for RMs with low interferences of REEs. Therefore, these carbonate and phosphate RMs can be used to achieve accurate and precise results for in situ Sr isotope investigations by LA‐MC‐ICP‐MS of similar materials.
Key Points
Sr isotope data for calcium carbonate and phosphate reference materials.
Comparison of laser ablation and solution MC‐ICP‐MS values with literature data.
Applicability of high‐precision data for validation and monitoring purposes for laser ablation work on carbonate and phosphate samples.
Abstract
The Late Cretaceous orogeny followed by the Eocene collision of the Adriatic with the European plate dissected the Northern Calcareous Alps (NCA) by a number of well-studied strike-slip ...fault systems accommodating N-S shortening and E-W stretching. However, the post-Miocene fault activity is poorly constrained due to lack of Neogene faulted sediments, and glacial erosion of geomorphic indicators. Using the protected environment of caves, we fill the knowledge gap in the post-Miocene evolution of the NCA by paleostress analysis of 172 reactivated faults that offset passages in 28 caves near major faults. Constrained maximum age of caves, our results indicate that the NCA have been subjected to N to NE trending compression since Pliocene. Faulted speleothems dated with
230
Th/U method, indicate that the recorded present-day stress state did not significantly change during the last 0.5 Ma. In contrast to the previously proposed post-Miocene N-S extension of NCA, but in agreement with what was observed in Vienna and Pannonian basins, we conclude that the eastward extrusion resulting from N-S convergence has continued despite a distinct slowdown of plate tectonic velocities in the late Miocene. The N-S extension affected only the Alpine front during Pliocene Molasse basin inversion, while at the scale of the Alpine orogen the NCA underwent successive N-S shortening and E-W stretching.
Abstract
Hells Bells are underwater secondary carbonates discovered in sinkholes (cenotes) southeast of Cancun on the north-eastern Yucatán peninsula, Mexico. These authigenic calcite precipitates, ...reaching a length of up to 4 m, most likely grow in the pelagic redoxcline. Here we report on detailed
230
Th/U-dating and in-depth geochemical and stable isotope analyses of specimens from cenotes El Zapote, Maravilla and Tortugas. Hells Bells developed since at least ~ 8000 years ago, with active growth until present day. Initial (
234
U/
238
U) activity ratios (δ
234
U
0
) in Hells Bells calcite decreas from 55 to 15‰ as sea level converges toward its present state. The temporal evolution of the geochemistry and isotope composition of Hells Bells calcites thus appears to be closely linked to sea-level rise and reflects changing hydrological conditions (desalinization) of the aquifer. We suggest that decelerated leaching of excess
234
U from the previously unsaturated bedrock traces Holocene relative sea-level rise. Considering this proxy, the resulting mean sea-level reconstruction contains half as much scatter, i.e. improves by a factor of two, when compared to previously published work for the period between 8 and 4 ky BP.