Evolution of the stable carbon isotope composition of atmospheric CO₂ over the last glacial cycle

We present new δ¹³C measurements of atmospheric CO₂ covering the last glacial/interglacial cycle, complementing previous records covering Terminations I and II. Most prominent in the new record is a significant depletion in δ¹³C(atm) of 0.5‰ occurring during marine isotope stage (MIS) 4, followed by an enrichment of the same magnitude at the beginning of MIS 3. Such a significant excursion in the record is otherwise only observed at glacial terminations, suggesting that similar processes were at play, such as changing sea surface temperatures, changes in marine biological export in the Southern Ocean (SO) due to variations in aeolian iron fluxes, changes in the Atlantic meridional overturning circulation, upwelling of deep water in the SO, and long-term trends in terrestrial carbon storage. Based on previous modeling studies, we propose constraints on some of these processes during specific time intervals. The decrease in δ¹³C(atm) at the end of MIS 4 starting approximately 64 kyr B.P. was accompanied by increasing [CO₂]. This period is also marked by a decrease in aeolian iron flux to the SO, followed by an increase in SO upwelling during Heinrich event 6, indicating that it is likely that a large amount of δ¹³C-depleted carbon was transferred to the deep oceans previously, i.e., at the onset of MIS 4. Apart from the upwelling event at the end of MIS 4 (and potentially smaller events during Heinrich events in MIS 3), upwelling of deep water in the SO remained reduced until the last glacial termination, whereupon a second pulse of isotopically light carbon was released into the atmosphere.

A reconstruction of atmospheric carbon dioxide and its stable carbon isotopic composition from the penultimate glacial maximum to the last glacial inception

The reconstruction of the stable carbon isotope evolution in atmospheric CO2 (δ13Catm), as archived in Antarctic ice cores, bears the potential to disentangle the contributions of the different carbon cycle fluxes causing past CO2 variations. Here we present a new record of δ13Catm before, during and after the Marine Isotope Stage 5.5 (155 000 to 105 000 yr BP). The dataset is archived on the data repository PANGEA® (www.pangea.de) under 10.1594/PANGAEA.817041. The record was derived with a well established sublimation method using ice from the EPICA Dome C (EDC) and the Talos Dome ice cores in East Antarctica. We find a 0.4‰ shift to heavier values between the mean δ13Catm level in the Penultimate (~ 140 000 yr BP) and Last Glacial Maximum (~ 22 000 yr BP), which can be explained by either (i) changes in the isotopic composition or (ii) intensity of the carbon input fluxes to the combined ocean/atmosphere carbon reservoir or (iii) by long-term peat buildup. Our isotopic data suggest that the carbon cycle evolution along Termination II and the subsequent interglacial was controlled by essentially the same processes as during the last 24 000 yr, but with different phasing and magnitudes. Furthermore, a 5000 yr lag in the CO2 decline relative to EDC temperatures is confirmed during the glacial inception at the end of MIS5.5 (120 000 yr BP). Based on our isotopic data this lag can be explained by terrestrial carbon release and carbonate compensation.

Can the carbon isotopic composition of methane be reconstructed from multi-site firn air measurements?

Methane is a strong greenhouse gas and large uncertainties exist concerning the future evolution of its atmospheric abundance. Analyzing methane atmospheric mixing and stable isotope ratios in air trapped in polar ice sheets helps in reconstructing the evolution of its sources and sinks in the past. This is important to improve predictions of atmospheric CH4 mixing ratios in the future under the influence of a changing climate. The aim of this study is to assess whether past atmospheric δ13C(CH4) variations can be reliably reconstructed from firn air measurements. Isotope reconstructions obtained with a state of the art firn model from different individual sites show unexpectedly large discrepancies and are mutually inconsistent. We show that small changes in the diffusivity profiles at individual sites lead to strong differences in the firn fractionation, which can explain a large part of these discrepancies. Using slightly modified diffusivities for some sites, and neglecting samples for which the firn fractionation signals are strongest, a combined multi-site inversion can be performed, which returns an isotope reconstruction that is consistent with firn data. However, the isotope trends are lower than what has been concluded from Southern Hemisphere (SH) archived air samples and high-accumulation ice core data. We conclude that with the current datasets and understanding of firn air transport, a high precision reconstruction of δ13C of CH4 from firn air samples is not possible, because reconstructed atmospheric trends over the last 50 yr of 0.3–1.5 ‰ are of the same magnitude as inherent uncertainties in the method, which are the firn fractionation correction (up to ~2 ‰ at individual sites), the Kr isobaric interference (up to ~0.8 ‰, system dependent), inter-laboratory calibration offsets (~0.2 ‰) and uncertainties in past CH4 levels (~0.5 ‰).

The stable carbon isotopic composition of Daphnia ephippia in small, temperate lakes reflects in-lake methane availability

Daphnia can ingest methane-oxidizing bacteria and incorporate methanogenic carbon into their biomass, leading to low stable carbon isotope ratios (expressed as δ13C values) of their tissue. Therefore, δ13C analysis of Daphnia resting eggs (ephippia) in lake sediment records can potentially be used to reconstruct past in-lake availability of methane (CH4). However, detailed multilake studies demonstrating that δ13C values of recently deposited Daphnia ephippia (δ13Cephippia) are systematically related to in-lake CH4 concentrations (CH4aq) are still missing. We measured δ13Cephippia from surface sediments of 15 small lakes in Europe, and compared these values with late-summer CH4aq. δ13Cephippia ranged from −51.6‰ to −25.9‰, and was strongly correlated with CH4aq in the surface water and above the sediment (r −0.73 and −0.77, respectively), whereas a negative rather than the expected positive correlation was found with δ13C values of carbon dioxide (CO2) (r −0.54), and no correlation was observed with CO2aq. At eight sites, offsets between δ13 CCO2 and δ13Cephippia exceeded offsets between δ13 CCO2 and δ13Calgae reported in literature. δ13Cephippia was positively correlated with δ13C values of sedimentary organic matter (r 0.54), but up to 20.7‰ lower in all except one of the lakes (average −6.1‰). We conclude that incorporation of methanogenic carbon prior to ephippia formation must have been widespread by Daphnia in our study lakes, especially those with high CH4aq. Our results suggest a systematic relationship between δ13Cephippia values and CH4aq in small temperate lakes, and that δ13Cephippia analysis on sediment records may provide insights into past changes in in-lake CH4aq.

Deep ocean ventilation, carbon isotopes, marine sedimentation and the deglacial CO2 rise

The link between the atmospheric CO2 level and the ventilation state of the deep ocean is an important building block of the key hypotheses put forth to explain glacial-interglacial CO2 fluctuations. In this study, we systematically examine the sensitivity of atmospheric CO2 and its carbon isotope composition to changes in deep ocean ventilation, the ocean carbon pumps, and sediment formation in a global three-dimensional ocean-sediment carbon cycle model. Our results provide support for the hypothesis that a break up of Southern Ocean stratification and invigorated deep ocean ventilation were the dominant drivers for the early deglacial CO2 rise of ~35 ppm between the Last Glacial Maximum and 14.6 ka BP. Another rise of 10 ppm until the end of the Holocene is attributed to carbonate compensation responding to the early deglacial change in ocean circulation. Our reasoning is based on a multi-proxy analysis which indicates that an acceleration of deep ocean ventilation during the early deglaciation is not only consistent with recorded atmospheric CO2 but also with the reconstructed opal sedimentation peak in the Southern Ocean at around 16 ka BP, the record of atmospheric δ13CCO2, and the reconstructed changes in the Pacific CaCO3 saturation horizon.

Stable N isotope composition of nitrate reflects N transformations during the passage of water through a montane rain forest in Ecuador

Knowledge of the fate of deposited N in the possibly N-limited, highly biodiverse north Andean forests is important because of the possible effects of N inputs on plant performance and species composition. We analyzed concentrations and fluxes of NO3 −–N, NH4 +–N and dissolved organic N (DON) in rainfall, throughfall, litter leachate, mineral soil solutions (0.15–0.30 m depths) and stream water in a montane forest in Ecuador during four consecutive quarters and used the natural 15N abundance in NO3 − during the passage of rain water through the ecosystem and bulk δ15N values in soil to detect N transformations. Depletion of 15N in NO3 − and increased NO3 −–N fluxes during the passage through the canopy and the organic layer indicated nitrification in these compartments. During leaching from the organic layer to mineral soil and stream, NO3 − concentrations progressively decreased and were enriched in 15N but did not reach the δ15N values of solid phase organic matter (δ15N = 5.6–6.7‰). This suggested a combination of nitrification and denitrification in mineral soil. In the wettest quarter, the δ15N value of NO3 − in litter leachate was smaller (δ15N = −1.58‰) than in the other quarters (δ15N = −9.38 ± SE 0.46‰) probably because of reduced mineralization and associated fractionation against 15N. Nitrogen isotope fractionation of NO3 − between litter leachate and stream water was smaller in the wettest period than in the other periods probably because of a higher rate of denitrification and continuous dilution by isotopically lighter NO3 −–N from throughfall and nitrification in the organic layer during the wettest period. The stable N isotope composition of NO3 − gave valuable indications of N transformations during the passage of water through the forest ecosystem from rainfall to the stream.

Deep ocean ventilation, carbon isotopes, marine sedimentation and the deglacial CO2 rise

The link between the atmospheric CO2 level and the ventilation state of the deep ocean is an important building block of the key hypotheses put forth to explain glacial-interglacial CO2 fluctuations. In this study, we systematically examine the sensitivity of atmospheric CO2 and its carbon isotope composition to changes in deep ocean ventilation, the ocean carbon pumps, and sediment formation in a global 3-D ocean-sediment carbon cycle model. Our results provide support for the hypothesis that a break up of Southern Ocean stratification and invigorated deep ocean ventilation were the dominant drivers for the early deglacial CO2 rise of ~35 ppm between the Last Glacial Maximum and 14.6 ka BP. Another rise of 10 ppm until the end of the Holocene is attributed to carbonate compensation responding to the early deglacial change in ocean circulation. Our reasoning is based on a multi-proxy analysis which indicates that an acceleration of deep ocean ventilation during early deglaciation is not only consistent with recorded atmospheric CO2 but also with the reconstructed opal sedimentation peak in the Southern Ocean at around 16 ka BP, the record of atmospheric δ13CCO2, and the reconstructed changes in the Pacific CaCO3 saturation horizon.

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.

Highly accurate isotope composition measurements by a miniature laser ablation mass spectrometer designed for in situ investigations on planetary surface

An experimental procedure for precise and accurate measurements of isotope abundances by a miniature laser ablation mass spectrometer for space research is described. The measurements were conducted on different untreated NIST standards and galena samples by applying pulsed UV laser radiation (266 nm, 3 ns and 20 Hz) for ablation, atomisation, and ionisation of the sample material. Mass spectra of released ions are measured by a reflectron-type time-of-flight mass analyser. A computer controlled performance optimiser was used to operate the system at maximum ion transmission and mass resolution. At optimal experimental conditions, the best relative accuracy and precision achieved for Pb isotope compositions are at the per mill level and were obtained in a range of applied laser irradiances and a defined number of accumulated spectra. A similar relative accuracy and precision was achieved in the study of Pb isotope compositions in terrestrial galena samples. The results for the galena samples are similar to those obtained with a thermal ionisation mass spectrometer (TIMS). The studies of the isotope composition of other elements yielded relative accuracy and precision at the per mill level too, with characteristic instrument parameters for each element. The relative accuracy and precision of the measurements is degrading with lower element/isotope concentration in a sample. For the elements with abundances below 100 ppm these values drop to the percent level. Depending on the isotopic abundances of Pb in minerals, 207Pb/206Pb ages with accuracy in the range of tens of millions of years can be achieved.

Coupling of LMS with a fs-laser ablation ion source: elemental and isotope composition measurements

Mass spectrometric analysis of elemental and isotopic compositions of several NIST standards is performed by a miniature laser ablation/ionisation reflectron-type time-of-flight mass spectrometer (LMS) using a fs-laser ablation ion source (775 nm, 190 fs, 1 kHz). The results of the mass spectrometric studies indicate that in a defined range of laser irradiance (fluence) and for a certain number of accumulations of single laser shot spectra, the measurements of isotope abundances can be conducted with a measurement accuracy at the per mill level and at the per cent level for isotope concentrations higher and lower than 100 ppm, respectively. Also the elemental analysis can be performed with a good accuracy. The LMS instrument combined with a fs-laser ablation ion source exhibits similar detection efficiency for both metallic and non-metallic elements. Relative sensitivity coefficients were determined and found to be close to one, which is of considerable importance for the development of standard-less instruments. Negligible thermal effects, sample damage and excellent characteristics of the fs-laser beam are thought to be the main reason for substantial improvement of the instrumental performance compared to other laser ablation mass spectrometers.

On the interference of Kr during carbon isotope analysis of methane using continuous-flow combustion–isotope ratio mass spectrometry

Stable carbon isotope analysis of methane (delta C-13 of CH4) on atmospheric samples is one key method to constrain the current and past atmospheric CH4 budget. A frequently applied measurement technique is gas chromatography (GC) isotope ratio mass spectrometry (IRMS) coupled to a combustion-preconcentration unit. This report shows that the atmospheric trace gas krypton (Kr) can severely interfere during the mass spectrometric measurement, leading to significant biases in delta C-13 of CH4, if krypton is not sufficiently separated during the analysis. According to our experiments, the krypton interference is likely composed of two individual effects, with the lateral tailing of the doubly charged Kr-86 peak affecting the neighbouring m/z 44 and partially the m/z 45 Faraday cups. Additionally, a broad signal affecting m/z 45 and especially m/z 46 is assumed to result from scattered ions of singly charged krypton. The introduced bias in the measured isotope ratios is dependent on the chromatographic separation, the krypton-to-CH4 mixing ratio in the sample, the focusing of the mass spectrometer as well as the detector configuration and can amount to up to several per mil in delta C-13. Apart from technical solutions to avoid this interference, we present correction routines to a posteriori remove the bias.