Toward explaining the Holocene carbon dioxide and carbon isotope records: Results from transient ocean carbon cycle-climate simulations

[1] The Bern3D model was applied to quantify the mechanisms of carbon cycle changes during the Holocene (last 11,000 years). We rely on scenarios from the literature to prescribe the evolution of shallow water carbonate deposition and of land carbon inventory changes over the glacial termination (18,000 to 11,000 years ago) and the Holocene and modify these scenarios within uncertainties. Model results are consistent with Holocene records of atmospheric CO2 and δ13C as well as the spatiotemporal evolution of δ13C and carbonate ion concentration in the deep sea. Deposition of shallow water carbonate, carbonate compensation of land uptake during the glacial termination, land carbon uptake and release during the Holocene, and the response of the ocean-sediment system to marine changes during the termination contribute roughly equally to the reconstructed late Holocene pCO2 rise of 20 ppmv. The 5 ppmv early Holocene pCO2 decrease reflects terrestrial uptake largely compensated by carbonate deposition and ocean sediment responses. Additional small contributions arise from Holocene changes in sea surface temperature, ocean circulation, and export productivity. The Holocene pCO2 variations result from the subtle balance of forcings and processes acting on different timescales and partly in opposite direction as well as from memory effects associated with changes occurring during the termination. Different interglacial periods with different forcing histories are thus expected to yield different pCO2 evolutions as documented by ice cores.

Multiple tree-ring proxies (earlywood width, latewood width and δ13C) from pedunculate oak (Quercus robur L.), Hungary

The aim of this study was to analyse the effects of climatic factors (i.e. monthly mean temperature and total precipitation) on radial growth (earlywood width, latewood width, and total ringwidth) and on latewood stable carbon isotope composition in a pedunculate oak (Quercus robur L) stand in northeastern Hungary. Earlywood widths showed the weakest common variance and lack of statistically significant relationship to monthly precipitation and temperature. Latewood width showed the strongest common chronological signal. Correlation analysis with the monthly climate series pointed out the strongest positive/negative correlation with June precipitation for latewood width/stable carbon isotope ratio. These parameters shared the strongest climatic response also for seasonal scale since the highest correlation coefficients, 0.49 and -0.62 for latewood width and stable carbon isotope ratio, respectively, were obtained for both with a 10-month precipitation total (from previous November to current August of the growing season). A combined parameter, derived as difference between latewood width and stable carbon isotope indices showed improved statistical relationship compared to the hydroclimatic calibration target both for local and regional spatial scales. Spatial correlation analysis indicated that the hydroclimatic signal encoded in these moisture sensitive tree-ring parameters from Bakta Forest is expected to be representative for the northeastern Carpathians and for the large part of the Great Hungarian Plain. In addition, the hydroclimatic signal of latewood width chronology was compared to three independent records. Results showed that neither the strength nor the rank of the similarity of the local hydroclimate signals were stable throughout the past two centuries. Future palaeo(hydro)climatological efforts targeting the Carpathian(-Balkan) region are recommended to track carefully the spatial domains for which a given, local, proxy-derived hydroclimate reconstruction might provide useful information.

Stable Isotope and Trace Element Studies on Gladiators and Contemporary Romans from Ephesus (Turkey, 2nd and 3rd Ct. AD) - Implications for Differences in Diet

The gladiator cemetery discovered in Ephesus (Turkey) in 1993 dates to the 2nd and 3rd century AD. The aim of this study is to reconstruct diverse diet, social stratification, and migration of the inhabitants of Roman Ephesus and the distinct group of gladiators. Stable carbon, nitrogen, and sulphur isotope analysis were applied, and inorganic bone elements (strontium, calcium) were determined. In total, 53 individuals, including 22 gladiators, were analysed. All individuals consumed C3 plants like wheat and barley as staple food. A few individuals show indication of consumption of C4 plants. The δ13C values of one female from the gladiator cemetery and one gladiator differ from all other individuals. Their δ34S values indicate that they probably migrated from another geographical region or consumed different foods. The δ15N values are relatively low in comparison to other sites from Roman times. A probable cause for the depletion of 15N in Ephesus could be the frequent consumption of legumes. The Sr/Ca-ratios of the gladiators were significantly higher than the values of the contemporary Roman inhabitants. Since the Sr/Ca-ratio reflects the main Ca-supplier in the diet, the elevated values of the gladiators might suggest a frequent use of a plant ash beverage, as mentioned in ancient texts.

Social stratigraphy in Late Iron Age Switzerland: Stable carbon, nitrogen and sulphur isotope analysis of human remains from Münsingen

The Iron Age cemetery of Münsingen in Switzerland with 220 abundantly equipped burials marked a milestone for Iron Age research. The horizontal spread throughout the time of its occupancy laid the foundation for the chronology system of the Late Iron Age. Today, skulls of 77 individuals and some postcranial bones are still preserved. The aim of the study was to obtain information about diet, mobility and social stratification of the individuals. Stable isotope ratios of carbon, nitrogen and sulphur were analysed for 63 individuals. For all of them C3 plants were found to be the staple food. There are significant differences between males and females in δ13C and δ15N values. This points to a gender restriction in the access to animal protein with males probably having more access to meat and dairy products. Differences in δ15N values were also observed for different age classes. δ34S values indicate a terrestrial-based diet with no significant intake of marine or freshwater fish. Seven adults with enriched δ34S values might have immigrated to Münsingen, four of which were found in the oldest part of the cemetery. Furthermore, possible changes of the vegetation are indicated by the more positive stable carbon ratios in the later phases. The results lead to the suggestion that especially males buried with weapons might have played a special role in the Iron Age society.

The stable isotopic composition of Daphnia ephippia reflects changes in δ13C and δ18O values of food and water

The stable isotopic composition of fossil resting eggs (ephippia) of Daphnia spp. is being used to reconstruct past environmental conditions in lake ecosystems. However, the underlying assumption that the stable isotopic composition of the ephippia reflects the stable isotopic composition of the parent Daphnia, of their diet and of the environmental water have yet to be confirmed in a controlled experimental setting. We performed experiments with Daphnia pulicaria cultures, which included a control treatment conducted at 12 °C in filtered lake water and with a diet of fresh algae and three treatments in which we manipulated the stable carbon isotopic composition (δ13C value) of the algae, stable oxygen isotopic composition (δ18O value) of the water and the water temperature, respectively. The stable nitrogen isotopic composition (δ15N value) of the algae was similar for all treatments. At 12 °C, differences in algal δ13C values and in δ18O values of water were reflected in those of Daphnia. The differences between ephippia and Daphnia stable isotope ratios were similar in the different treatments (δ13C: +0.2 ± 0.4 ‰ (standard deviation); δ15N: −1.6 ± 0.4 ‰; δ18O: −0.9 ± 0.4 ‰), indicating that changes in dietary δ13C values and in δ18O values of water are passed on to these fossilizing structures. A higher water temperature (20 °C) resulted in lower δ13C values in Daphnia and ephippia than in the other treatments with the same food source and in a minor change in the difference between δ13C values of ephippia and Daphnia (to −1.3 ± 0.3 ‰). This may have been due to microbial processes or increased algal respiration rates in the experimental containers, which may not affect Daphnia in natural environments. There was no significant difference in the offset between δ18O and δ15N values of ephippia and Daphnia between the 12 and 20 °C treatments, but the δ18O values of Daphnia and ephippia were on average 1.2 ‰ lower at 20 °C than at 12 °C. We conclude that the stable isotopic composition of Daphnia ephippia provides information on that of the parent Daphnia and of the food and water they were exposed to, with small offsets between Daphnia and ephippia relative to variations in Daphnia stable isotopic composition reported from downcore studies. However, our experiments also indicate that temperature may have a minor influence on the δ13C, δ15N and δ18O values of Daphnia body tissue and ephippia. This aspect deserves attention in further controlled experiments.

Temperature and precipitation effects on δ13C depth profiles in SOM under temperate beech forests

Enrichment of 13C in SOM with soil depth is related to interacting processes influenced by temperature and precipitation. Our objectives were to derive climate effects on patterns of vertical δ13C values of soil organic matter (SOM) while minimizing the effect of confounding variables. We investigated vertical changes in δ13C values of SOM in 1-cm depth intervals in silvicultural mature beech (Fagus sylvatica L.) forest ecosystems in northern Rhineland-Palatinate across gradients of MAT (7.9 to 9.7 °C mean annual temperature) and MAP (607 to 1085 mm mean annual precipitation) in winter 2011. Forest stands (n = 10) were chosen based on data sets provided by the Rhineland-Palatinate Forest Administration so that variations in these gradients occurred while other environmental factors like physico-chemical soil properties, tree species, stand age, exposition and precipitation (for the temperature gradient) or temperature (for the precipitation gradient) did not differ among study sites. From litter down to the mineral soil at 10 cm depth, soil organic carbon (SOC) content decreased (47.5 ± SE 0.1% to 2.5 ± 0.1%) while the δ13C values increased (− 29.4 ± 0.1‰ to − 26.1 ± 0.1‰). Litter of sites under higher MAP/lower MAT had lower δ13C values which was in line with literature data on climate driven plant physiological process. To compare the dimension of the vertical 13C enrichment, δ13C values were regressed linearly against log-transformed carbon contents yielding absolute values of these slopes (beta). Beta values ranged between 0.6 and 4.5 (range of r from − 0.7 to − 1.0; p < 0.01). Due to an assumed decay continuum and similar variations of δ13C values in litter and in 10 cm depth, we conclude that effects on isotope composition in the Oi layer continue vertically and therefore, δ13C values in litter do not solely control beta values. Beta values decreased with increasing MAT (r = − 0.83; p < 0.05). Reduced soil moisture and therefore both, reduced microbial activity and reduced downward transport of microbial cycled DOM (=13C enriched) might be responsible for less pronounced δ13C depth profiles in case of high temperatures. Greater C:N ratios (lower degradability) of the litter under higher temperatures likely contributed to these depth trends. Beta values increased with increasing MAP (r = 0.73; p < 0.05). We found decreasing C:N ratios in the mineral soil that possibly indicates higher decomposition under higher precipitation. Exclusion of the organic layers from linear regressions indicated a stronger impact of MAP on the development of δ13C depth profiles. Our results confirm temperature and precipitation effects on δ13C depth profiles and indicate stronger 13C enrichment under lower MAT/higher MAP. Therefore, time series of vertical δ13C depth profiles might provide insights into climate change effects.

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.

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.

Eclogitic diamonds from variable crustal protoliths in the northeastern Siberian craton: trace elements and coupled δ13C – δ18O signatures in diamonds and garnet inclusions

Diamonds of eclogitic assemblages are dominant in the placer diamond deposits of the northeastern Siberian platform. In this study we present new trace elements and stable isotopes (δ13C and δ18O) data for alluvial diamonds and their garnet inclusions from this locality. Cr-rich garnets of peridotitic affinity in the studied diamonds have a narrow range of δ18O values from 5.7‰ to 6.2‰, which is largely overlapping with the accepted mantle range. This narrow range suggests that the garnet inclusions showing different REE patterns and little variations in oxygen isotopes may have formed by different processes involving fluid/melts that, however, were in oxygen isotopic equilibrium with the mantle. The trace element composition of the eclogitic garnet inclusions supports a crustal origin for at least the high-Ca garnets, which show flat HREE patterns and in some cases a positive Eu-anomaly. High-Ca eclogitic garnets generally show heavier oxygen isotope compositions (δ18O 6.5–9.6‰) than what is observed in low-Ca garnets (δ18O 5.7–7.4‰). The variability in oxygen isotopes and trace elements is suggested to be inherited from contrasting crustal protoliths. The relationship between the high δ18O values of inclusions and the low δ13C values of the host diamonds implies that the high-Ca garnet inclusions were derived from intensely hydrated (e.g., δ18O > 7‰) and typically oxidised basaltic rock close to the seawater interface, and that the carbon for diamonds was closely associated with this protolith.

Real-time analysis of δ13C- and δD-CH4 in ambient air with laser spectroscopy: method development and first intercomparison results

In situ and simultaneous measurement of the three most abundant isotopologues of methane using mid-infrared laser absorption spectroscopy is demonstrated. A field-deployable, autonomous platform is realized by coupling a compact quantum cascade laser absorption spectrometer (QCLAS) to a preconcentration unit, called trace gas extractor (TREX). This unit enhances CH4 mole fractions by a factor of up to 500 above ambient levels and quantitatively separates interfering trace gases such as N2O and CO2. The analytical precision of the QCLAS isotope measurement on the preconcentrated (750 ppm, parts-per-million, µmole mole−1) methane is 0.1 and 0.5 ‰ for δ13C- and δD-CH4 at 10 min averaging time. Based on repeated measurements of compressed air during a 2-week intercomparison campaign, the repeatability of the TREX–QCLAS was determined to be 0.19 and 1.9 ‰ for δ13C and δD-CH4, respectively. In this intercomparison campaign the new in situ technique is compared to isotope-ratio mass spectrometry (IRMS) based on glass flask and bag sampling and real time CH4 isotope analysis by two commercially available laser spectrometers. Both laser-based analyzers were limited to methane mole fraction and δ13C-CH4 analysis, and only one of them, a cavity ring down spectrometer, was capable to deliver meaningful data for the isotopic composition. After correcting for scale offsets, the average difference between TREX–QCLAS data and bag/flask sampling–IRMS values are within the extended WMO compatibility goals of 0.2 and 5 ‰ for δ13C- and δD-CH4, respectively. This also displays the potential to improve the interlaboratory compatibility based on the analysis of a reference air sample with accurately determined isotopic composition.

Northern Mediterranean climate since the Middle Pleistocene: a 637 ka stable isotope record from Lake Ohrid (Albania/Macedonia)

Lake Ohrid (Macedonia/Albania) is an ancient lake with unique biodiversity and a site of global significance for investigating the influence of climate, geological, and tectonic events on the generation of endemic populations. Here, we present oxygen (δ18O) and carbon (δ13C) isotope data from carbonate over the upper 243 m of a composite core profile recovered as part of the Scientific Collaboration on Past Speciation Conditions in Lake Ohrid (SCOPSCO) project. The investigated sediment succession covers the past ca. 637 ka. Previous studies on short cores from the lake (up to 15 m, < 140 ka) have indicated the total inorganic carbon (TIC) content of sediments to be highly sensitive to climate change over the last glacial–interglacial cycle. Sediments corresponding to warmer periods contain abundant endogenic calcite; however, an overall low TIC content in glacial sediments is punctuated by discrete bands of early diagenetic authigenic siderite. Isotope measurements on endogenic calcite (δ18Oc and δ13Cc) reveal variations both between and within interglacials that suggest the lake has been subject to palaeoenvironmental change on orbital and millennial timescales. We also measured isotope ratios from authigenic siderite (δ18Os and δ13Cs) and, with the oxygen isotope composition of calcite and siderite, reconstruct δ18O of lake water (δ18Olw) over the last 637 ka. Interglacials have higher δ18Olw values when compared to glacial periods most likely due to changes in evaporation, summer temperature, the proportion of winter precipitation (snowfall), and inflow from adjacent Lake Prespa. The isotope stratigraphy suggests Lake Ohrid experienced a period of general stability from marine isotope stage (MIS) 15 to MIS 13, highlighting MIS 14 as a particularly warm glacial. Climate conditions became progressively wetter during MIS 11 and MIS 9. Interglacial periods after MIS 9 are characterised by increasingly evaporated and drier conditions through MIS 7, MIS 5, and the Holocene. Our results provide new evidence for long-term climate change in the northern Mediterranean region, which will form the basis to better understand the influence of major environmental events on biological evolution within Lake Ohrid.

Stable isotope response to lake eutrophication: Calibration of a high-resolution lacustrine sequence from Baldeggersee, Switzerland

Stable isotope analyses of discrete seasonal layers from a 108-yr annually laminated freeze-core from Baldeg-gersee, a small, eutrophic lake in central Switzerland, provide information on the climatological and environmental factors, including lake eutrophication, that control oxygen and carbon isotopic composition of epilimnic biologically induced calcite precipitate. During the last 100 yr, Baldeggersee has undergone major increases in productivity and eutrophication in response to nutrient loading from agriculture and industrialization in the lake's watershed. Calibration of the isotopic signal in Baldeggersee to historical limnological data quantitatively links evidence of isotopic depletion in the sedimented calcite to trophic state of the lake. δ18O values from the spring/summer “light” sediment layers steadily diverged to more depleted values in response to historical eutrophication: measured δ18O values were up to 21.5‰ more negative than calculated equilibrium δ18O values. Evidence for 13C depletion in the calcite, relative to equilibrium values, is more difficult to ascertain because of an overall dominance of isotopic enrichment in the dissolved inorganic pool as productivity in Baldeggersee increases. A positive association exists between the degree of oxygen-18 depletion and the calcite crystal size. Thus, large amorphous calcite grains can be used as a proxy for recognizing apparent isotopic nonequilibrium in sediment sequences from highly productive lacustrine environments from all geologic time scales. In contrast to the light layers, the oxygen isotopic composition of the calcite in the late summer/fall “dark” sediment layers is unaffected by the apparent isotope nonequilibrium. Oxygen and carbon isotope values from the dark laminae in the Baldeggersee sediment therefore provide environmental and climatological proxies that can be calibrated with known environmental and regional climate data for the last century.

Benthic foraminiferal stable isotope record from ODP Site 138-849 on the East Pacific Rise

Benthic foraminifer and delta13C data from Site 849, on the west flank of the East Pacific Rise (0°11 'N, 110°31'W; 3851 m), give relatively continuous records of deep Pacific Ocean stable isotope variations between 0 and 5 Ma. The mean sample spacing is 4 k.y. Most analyses are from Cibicides wuellerstorfi, but isotopic offsets relative to Uvigerina peregrina appear roughly constant. Because of its location west of the East Pacific Rise, Site 849 yields a suitable record of mean Pacific Ocean delta13C, which approximates a global oceanic signal. The ~100-k.y.-period climate cycle, which is prevalent in delta18O does not dominate the long-term delta13C record. For delta13C, variations in the ~400- and 41-k.y. periods are more important. Phase lags of delta13C relative to ice volume in the 41- and 23-k.y. bands are consistent with delta13C as a measure of organic biomass. A model-calculated exponential response time of 1-2 k.y. is appropriate for carbon stored in soils and shallow sediments responding to glacial-interglacial climate change. Oceanic delta13C leads ice volume slightly in the 100-k.y. band, and this suggests another process such as changes in continental weathering to modulate mean river delta13C at long periods. The delta13C record from Site 849 diverges from that of Site 677 in the Panama Basin mostly because of decay of 13C-depleted organic carbon in the relatively isolated Panama Basin. North Atlantic to Pacific delta13C differences calculated using published data from Sites 607 and 849 reveal variations in Pliocene deep water within the range of those of the late Quaternary. Maximum delta13C contrast between these sites, which presumably reflects maximum influx of high-delta13C northern source water into the deep North Atlantic Ocean, occurred between 1.3 and 2.1 Ma, well after the initiation of Northern Hemisphere glaciation. Export of high-delta13C North Atlantic Deep Water from the Atlantic to the circumpolar Antarctic, as recorded by published delta13C data from Subantarctic Site 704, appears unrelated to the North Atlantic-Pacific delta13C contrast. To account for this observation, we suggest that deep-water formation in the North Atlantic reflects northern source characteristics, whereas export of this water into the circumpolar Antarctic reflects Southern Hemisphere wind forcing. Neither process appears directly linked to ice-volume variations.

Great Barrier Reef coral element/Ca and stable isotope data and U-Th ages from IODP Expedition 325

Tropical south-western Pacific temperatures are of vital importance to the Great Barrier Reef (GBR), but the role of sea surface temperatures (SSTs) in the growth of the GBR since the Last Glacial Maximum remains largely unknown. Here we present records of Sr/Ca and d18O for Last Glacial Maximum and deglacial corals that show a considerably steeper meridional SST gradient than the present day in the central GBR. We find a 1-2 °C larger temperature decrease between 17° and 20°S about 20,000 to 13,000 years ago. The result is best explained by the northward expansion of cooler subtropical waters due to a weakening of the South Pacific gyre and East Australian Current. Our findings indicate that the GBR experienced substantial meridional temperature change during the last deglaciation, and serve to explain anomalous deglacial drying of northeastern Australia. Overall, the GBR developed through significant SST change and may be more resilient than previously thought.