Estimation and calibration of the water isotope differential diffusion length in ice core records

Palaeoclimatic information can be retrieved from the diffusion of the stable water isotope signal during firnification of snow. The diffusion length, a measure for the amount of diffusion a layer has experienced, depends on the firn temperature and the accumulation rate. We show that the estimation of the diffusion length using power spectral densities (PSDs) of the record of a single isotope species can be biased by uncertainties in spectral properties of the isotope signal prior to diffusion. By using a second water isotope and calculating the difference in diffusion lengths between the two isotopes, this problem is circumvented. We study the PSD method applied to two isotopes in detail and additionally present a new forward diffusion method for retrieving the differential diffusion length based on the Pearson correlation between the two isotope signals. The two methods are discussed and extensively tested on synthetic data which are generated in a Monte Carlo manner. We show that calibration of the PSD method with this synthetic data is necessary to be able to objectively determine the differential diffusion length. The correlation-based method proves to be a good alternative for the PSD method as it yields precision equal to or somewhat higher than the PSD method. The use of synthetic data also allows us to estimate the accuracy and precision of the two methods and to choose the best sampling strategy to obtain past temperatures with the required precision. In addition to application to synthetic data the two methods are tested on stable-isotope records from the EPICA (European Project for Ice Coring in Antarctica) ice core drilled in Dronning Maud Land, Antarctica, showing that reliable firn temperatures can be reconstructed with a typical uncertainty of 1.5 and 2 °C for the Holocene period and 2 and 2.5 °C for the last glacial period for the correlation and PSD method, respectively.

Triple isotope (δD, δ17O, δ18O) study on precipitation, drip water and speleothem fluid inclusions for a Western Central European cave (NW Switzerland)

Deuterium (δD) and oxygen (δ18O) isotopes are powerful tracers of the hydrological cycle and have been extensively used for paleoclimate reconstructions as they can provide information on past precipitation, temperature and atmospheric circulation. More recently, the use of δ17O excess derived from precise measurement of δ17O and δ18O gives new and additional insights in tracing the hydrological cycle whereas uncertainties surround this proxy. However, 17O excess could provide additional information on the atmospheric conditions at the moisture source as well as about fractionations associated with transport and site processes. In this paper we trace water stable isotopes (δD,δ17O and δ18O) along their path from precipitation to cave drip water and finally to speleothem fluid inclusions for Milandre cave in northwestern Switzerland. A two year-long daily resolved precipitation isotope record close to the cave site is compared to collected cave drip water (3 months average resolution) and fluid inclusions of modern and Holocene stalagmites. Amount weighted mean δD,δ18O and δ17O are -71.0‰, -9.9‰, -5.2‰ for precipitation, -60.3‰, -8.7‰, -4.6‰ for cave drip water and -61.3‰, -8.3‰, -4.7‰ for recent fluid inclusions respectively. Second order parameters have also been derived in precipitation and drip water and present similar values with 18 per meg for 17O excess whereas d-excess is 1.5‰ more negative in drip water. Furthermore, the atmospheric signal is shifted towards enriched values in the drip water and fluid inclusions (Δ of ~ + 10‰ for δD). The isotopic composition of cave drip water exhibits a weak seasonal signal which is shifted by around 8 - 10 months (groundwater residence time) when compared to the precipitation. Moreover, we carried out the first δ17O measurement in speleothem fluid inclusions, as well as the first comparison of the δ17 O behaviour from the meteoric water to the fluid inclusions entrapment in speleothems. This study on precipitation, drip water and fluid inclusions will be used as a speleothem proxy calibration for Milandre cave in order to reconstruct paleotemperatures and moisture source variations for Western Central Europe.

Diet and mobility of an Iron Age population in Switzerland; Stable carbon, nitrogen and sulphur isotope analysis of the human remains from Münsingen

The 220 abundantly equipped burials from the Late Iron Age cemetery of Münsingen (420 – 240 BC) marked a milestone for Iron Age research. The evident horizontal spread throughout the time of occupancy laid the foundation for the chronology system of the Late Iron Age. Today the skulls of 77 individuals and some postcranial bones are still preserved. The aim was to obtain information about nutrition, social stratification and migration of the individuals from Münsingen. Stable isotope ratios of carbon, nitrogen and sulphur were analysed. The results of 63 individuals show that all consumed C3 plants as staple food with significant differences between males and females in δ13C and δ15N values. The results indicate a gender restriction in access to animal protein. Stable isotope values of one male buried with weapons and meat as grave goods suggest a diet with more animal proteins than the other individuals. It is possible that he was privileged due to high status. Furthermore, the δ34S values indicate minor mobility. Assuming that the subadults represent the local signal of δ34S it is very likely that adults with enriched δ34S could have migrated to Münsingen at some point during their lives. This study presents stable isotope values of one of the most important Late Iron Age burial sites in Central Europe. The presented data provide new insight into diet, migration and social stratification of the population from Münsingen.

Can stable isotope analysis and tooth histology assist in detecting stillbirth? An approach using neonatal remains from the Roman settlement Studen Petinesca

Infant burials in Roman settlements are a common observation. Even though ancient authors provide information many questions remain uncertain. For instance, the burial ritual for stillbirth and infanticide neonates is not specifically mentioned. This study therefore aimed to investigate the application of stable nitrogen (δ15N) and carbon (δ13C) isotopes from neonatal bone collagen in differentiating between a breastfeeding signal and stillbirth or a short survival of less than ten days. For this purpose collagen of 11 human and 14 non-human bones from the Roman settlement Petinesca (1st - 3rd century AD, Switzerland) was extracted and analysed for δ15N and δ13C. Tooth histology was performed for the central incisor and canine of the right mandible in order to investigate the presence of a neonatal line. According to the length of the long bones the age varied between 8.5 lunar months to 2 months ex utero. The stable isotope results provided a breastfeeding signal for all except one individual where the breastfeeding signal was absent. The tooth histological analysis of this individual exhibited no neonatal line. It is concluded that stable isotope analysis could indicate stillbirth or a short survival after birth. The tooth histology confirmed the stable isotope results. Furthermore, this might indicate that the burial ritual did not differentiate between stillbirth and neonates, who died within the time span stated by ancient authors of up to 40 days of age or the appearance of teeth. However, for further justifications additional research is going to be conducted.

Oxygen isotope record of oceanic and high-pressure metasomatism: a P-T-time-fluid path for the Monviso eclogites (Italy)

Fluids are considered a fundamental agent for chemical exchanges between different rock types in the subduction system. Constraints on the sources and pathways of subduction fluids thus provide crucial information to reconstruct subduction processes. The Monviso ophiolitic sequence is composed of mafic, ultramafic and minor sediments that have been subducted to ~80 km depth. In this sequence, both localized fluid flow and channelized fluids along major shear zones have been documented. We investigate the timing and source of the fluids that affected the dominant mafic rocks using microscale U-Pb dating of zircon and oxygen isotope analysis of mineral zones (garnet, zircon and antigorite) in high pressure rocks with variable degree of metasomatic modification. In mafic eclogites, Jurassic zircon cores are the only mineralogical relicts of the protolith gabbros and retain δ18O values of 4.5–6 ‰, typical of mantle melts. Garnet and metamorphic zircon that grew during prograde to peak metamorphism display low δ18O values between 0.2 and 3.8 ‰, which are likely inherited from high-temperature alteration of the protolith on the sea floor. This is corroborated by δ18O values of 3.0 and 3.6 ‰ in antigorite from surrounding serpentinites. In metasomatised eclogites within the Lower Shear Zone, garnet rim formed at the metamorphic peak shows a shift to higher δ18O up to 6‰. The age of zircons in high-pressure veins and metasomatised eclogites constrains the timing of fluid flow at high pressure at around 45–46 Ma. Although the oxygen data do not contradict previous reports of interaction with serpentinite-derived fluids, the shift to isotopically heavier oxygen compositions requires contribution from sediment-derived fluids. The scarcity of metasediments in the Monviso sequence suggests that such fluids were concentrated and fluxed along the Lower Shear Zone in a sufficient amount to modify the oxygen composition of the eclogitic minerals.

Qualitative Distinction of Autotrophic and Heterotrophic Processes at the Leaf Level by Means of Triple Stable Isotope (C–O–H) Patterns

Foliar samples were harvested from two oaks, a beech, and a yew at the same site in order to trace the development of the leaves over an entire vegetation season. Cellulose yield and stable isotopic compositions (δ13C, δ18O, and δD) were analyzed on leaf cellulose. All parameters unequivocally define a juvenile and a mature period in the foliar expansion of each species. The accompanying shifts of the δ13C-values are in agreement with the transition from remobilized carbohydrates (juvenile period), to current photosynthates (mature phase). While the opponent seasonal trends of δ18O of blade and vein cellulose are in perfect agreement with the state-of-art mechanistic understanding, the lack of this discrepancy for δD, documented for the first time, is unexpected. For example, the offset range of 18 permil (oak veins) to 57 permil (oak blades) in δD may represent a process driven shift from autotrophic to heterotrophic processes. The shared pattern between blade and vein found for both oak and beech suggests an overwhelming metabolic isotope effect on δD that might be accompanied by proton transfer linked to the Calvin-cycle. These results provide strong evidence that hydrogen and oxygen are under different biochemical controls even at the leaf level.

Chlorine isotope composition in chlorofluorocarbons CFC-11, CFC-12 and CFC-113 in firn, stratospheric and tropospheric air

The stratospheric degradation of chlorofluorocarbons (CFCs) releases chlorine, which is a major contributor to the destruction of stratospheric ozone (O3). A recent study reported strong chlorine isotope fractionation during the breakdown of the most abundant CFC (CFC-12, CCl2F2, Laube et al., 2010a), similar to effects seen in nitrous oxide (N2O). Using air archives to obtain a long-term record of chlorine isotope ratios in CFCs could help to identify and quantify their sources and sinks. We analyse the three most abundant CFCs and show that CFC-11 (CCl3F) and CFC-113 (CClF2CCl2F) exhibit significant stratospheric chlorine isotope fractionation, in common with CFC-12. The apparent isotope fractionation (ϵapp) for mid- and high-latitude stratospheric samples are respectively −2.4 (0.5) and −2.3 (0.4) ‰ for CFC-11, −12.2 (1.6) and −6.8 (0.8) ‰ for CFC-12 and −3.5 (1.5) and −3.3 (1.2) ‰ for CFC-113, where the number in parentheses is the numerical value of the standard uncertainty expressed in per mil. Assuming a constant isotope composition of emissions, we calculate the expected trends in the tropospheric isotope signature of these gases based on their stratospheric 37Cl enrichment and stratosphere–troposphere exchange. We compare these projections to the long-term δ (37Cl) trends of all three CFCs, measured on background tropospheric samples from the Cape Grim air archive (Tasmania, 1978–2010) and tropospheric firn air samples from Greenland (North Greenland Eemian Ice Drilling (NEEM) site) and Antarctica (Fletcher Promontory site). From 1970 to the present day, projected trends agree with tropospheric measurements, suggesting that within analytical uncertainties, a constant average emission isotope delta (δ) is a compatible scenario. The measurement uncertainty is too high to determine whether the average emission isotope δ has been affected by changes in CFC manufacturing processes or not. Our study increases the suite of trace gases amenable to direct isotope ratio measurements in small air volumes (approximately 200 mL), using a single-detector gas chromatography–mass spectrometry (GC–MS) system.

A rapid and efficient ion-exchange chromatography for Lu–Hf, Sm–Nd, and Rb–Sr geochronology and the routine isotope analysis of sub-ng amounts of Hf by MC-ICP-MS

The development and improvement of MC-ICP-MS instruments have fueled the growth of Lu–Hf geochronology over the last two decades, but some limitations remain. Here, we present improvements in chemical separation and mass spectrometry that allow accurate and precise measurements of 176Hf/177Hf and 176Lu/177Hf in high-Lu/Hf samples (e.g., garnet and apatite), as well as for samples containing sub-nanogram quantities of Hf. When such samples are spiked, correcting for the isobaric interference of 176Lu on 176Hf is not always possible if the separation of Lu and Hf is insufficient. To improve the purification of Hf, the high field strength elements (HFSE, including Hf) are first separated from the rare earth elements (REE, including Lu) on a first-stage cation column modified after Patchett and Tatsumoto (Contrib. Mineral. Petrol., 1980, 75, 263–267). Hafnium is further purified on an Ln-Spec column adapted from the procedures of Münker et al. (Geochem., Geophys., Geosyst., 2001, DOI: 10.1029/2001gc000183) and Wimpenny et al. (Anal. Chem., 2013, 85, 11258–11264) typically resulting in Lu/Hf < 0.0001, Zr/Hf < 1, and Ti/Hf < 0.1. In addition, Sm–Nd and Rb–Sr separations can easily be added to the described two-stage ion-exchange procedure for Lu–Hf. The isotopic compositions are measured on a Thermo Scientific Neptune Plus MC-ICP-MS equipped with three 1012 Ω resistors. Multiple 176Hf/177Hf measurements of international reference rocks yield a precision of 5–20 ppm for solutions containing 40 ppb of Hf, and 50–180 ppm for 1 ppb solutions (=0.5 ng sample Hf 0.5 in ml). The routine analysis of sub-ng amounts of Hf will facilitate Lu–Hf dating of low-concentration samples.

Geochronology and isotope geochemistry of Eocene dykes intruding the Ladakh Batholith

In order to determine the extent and timing of dyke formation in the Ladakh Batholith we examined about 30 mostly andesitic dykes intruding the Ladakh batholith in a ca. 50 km wide area to the west of Leh (NW India). The dykes in the east of the area trend E-NE and those in the west trend N-NW. The difference in orientation is also evident in the petrography and isotopic signatures. The eastern dykes contain corroded quartz xenocrysts and show negative ε0(Nd) and positive ε0(Sr) values, where as the western dykes do not contain quartz xenocrysts and exhibit positive ε0(Nd) and near-zero ε0(Sr) values. The variability in Sr-Nd isotopes (ε0(Nd) = 3.6 to −9.6, ε0(Sr) = 0.4 to 143) and the quartz xenocrysts can best be explained by (differing degrees of) crustal assimilation of the parent magma of the dykes. Separated minerals from five dykes were dated by 40Ar-39Ar incremental heating: amphibole ages range between 50 and 54 Ma, and one biotite dated both by Rb-Sr and by 40Ar-39Ar gave an age of 45 Ma. One dated pseudotachylyte sample attests to brittle faulting at ca. 54 Ma. The combination of structural field evidence with petrographic, isotopic and geochronological analyses demonstrates that the dykes did not form from a single, progressively differentiating magma chamber, despite having formed in the same tectonic setting around the same time, and that processes such as crustal assimilation and magma mixing/mingling also played a significant role in magma petrogenesis.

Constraints on barium isotope fractionation during aragonite precipitation by corals

We present a barium (Ba) isotope fractionation study of marine biogenic carbonates (aragonitic corals). The major aim is to provide first constraints on the Ba isotope fractionation between modern surface sea water and coral skele- ton. Mediterranean surface sea water was found to be enriched in the heavy Ba isotopes compared to previously reported values for marine open ocean authi- genic and terrestrial minerals. In aquarium experiments with a continuous sup- ply of Mediterranean surface water, the Ba isotopic composition of the bulk sample originating from cultured, aragonitic scleractinian corals (d137/134Ba between +0.16 +/- 0.12permil and +0.41 +/-0.12permil) were isotopically identical or lighter than that of the ambient Mediterranean surface sea water (d137/134Ba = +0.42 +/- 0.07permil, 2SD), which corresponds to an empirical maximum value of Ba isotope fractionation of D137/134Bacoral-seawater = -0.26 +/- 0.14permil at 25°C. This maximum Ba isotope fractionation is close and identical in direction to previous results from inorganic precipitation experiments with aragonite- structured pure BaCO3 (witherite). The variability in measured Ba concentrations of the cultured corals is at odds with a uniform distribution coefficient, DBa/Ca, thus indicating stronger vital effects on isotope than element discrimination. This observation supports the hypothesis that the Ba isotopic compositions of these corals do not result from simple equilibrium between the skeleton and the bulk sea water. Complementary coral samples from natural settings (tropical shallow-water corals from the Bahamas and Florida and cold- water corals from the Norwegian continental shelf) show an even wider range in d137/134Ba values (+0.14 +/- 0.08permil and +0.77 +/- 0.11permil), most probably due to additional spatial and/or temporal sea water heterogeneity, as indicated by recent publications.