A high resolution record of atmospheric carbon dioxid an its stable carb n isotopic composition from the penultimate glacial maximum to the glacial inception

Introduction Conclusions References


Introduction
Various processes are known to influence changes in the carbon distribution and its isotopic signature between the ocean, the atmosphere, terrestrial and marine organic carbon, reactive sediments and the lithosphere.Multiple processes operate simultaneously, and interact with each other non-linearly (K öhler et al., 2005;Brovkin et al., 2007;Sigman et al., 2010;Fischer et al., 2010;Tschumi et al., 2011) allowing for a wide range of possible scenarios to explain observed natural changes in atmospheric CO 2 .Thus, an unequivocal interpretation of past variations in the global carbon cycle is difficult.The stable carbon isotope signal of atmospheric CO 2 (δ 13 C atm ) represents a valuable tool to constrain processes affecting the global carbon cycle.Scrutinizing the potential processes and their contributions to the observed CO 2 variations, using longterm δ 13 C atm data sets represents an objective way to analyse the carbon cycle of the past.Antarctic ice cores allow for such reconstructions, however, the small changes in δ 13 C atm demand very high-precision measurements.Furthermore, due to the slow gas enclosure process, atmospheric changes are low-pass filtered in polar ice and fast changes in CO 2 and δ 13 C atm are damped (Joos and Spahni, 2008;K öhler et al., 2006, 2010, 2011), where the rate of damping is anti-correlated to the temperature and accumulation rate at the particular ice core site.Hence, a measurement precision significantly better than 0.1 ‰ is crucial for δ 13 C atm data, to provide useful constraints for global carbon cycle variability.Together with the small sample size available in polar ice cores, this represents a formidable analytical challenge.Apart from the requirement for high precision δ 13 C atm data, long-term coverage becomes crucial, since the balance of silicate and CaCO 3 weathering on land and of the preservation and dissolution response of CaCO 3 , opal and organic carbon in the deep ocean influence the CO 2 and δ 13 C atm dynamics over many thousand years (Elsig et al., 2009;Broecker and Clark, 2003).Since the time constants of these processes are comparable to the pacing of the major climate perturbations acting on orbital scales, neither the ocean's alkalinity budget nor the δ 13 C signature of the reservoirs necessarily reach a steady state.Previously published δ 13 C atm data focused mainly on time periods of distinct changes in atmospheric CO 2 including glacial terminations (Lourantou et al., 2010a,b;Schmitt et al., 2012;Smith et al., 1999) or the late rise of the CO 2 concentrations in the course of the Holocene (Inderm ühle et al., 1999;Elsig et al., 2009).
Here, we present a highly resolved δ 13 C atm record fulfilling the requirements in precision and covering the full Penultimate Glacial Maximum (PGM), the penultimate glacial/interglacial transition, the interglacial Marine Isotope Stage (MIS) 5.5, and the following glacial inception to an age of about 105 000 yr BP.Our record consists of Figures measurements from two Antarctic ice cores, the European Project for Ice Coring in Antarctica Dome C (EDC) and the Talos Dome ice core.The discussion of these data focuses on four main issues: (1) the critical evaluation of rapid changes in δ 13 C atm as found by Lourantou et al. (2010b) at the end of the penultimate glacial/interglacial transition; (2) the similarities and differences of this work to the previously documented glacial transition from the Last Glacial Maximum (LGM) to the Holocene (Schmitt et al., 2012); (3) the 5000 yr lag of the CO 2 decline with respect to decreasing Antarctic temperature during the glacial inception at the end of MIS5.5; and (4) long term changes of atmospheric δ 13 C atm .

Ice core details
The data presented here were obtained from the EDC ( 74• 39 S, 124 • 10 E, elevation 3240 m) and Talos Dome (72 • 49 S, 159 • 11 E, elevation 2315 m) ice cores.EDC samples were taken from depths between 1454 to 1887 m corresponding to the time interval between 104 600 and 154 500 yr BP on the EDC3 gas age scale (Parrenin et al., 2007;Loulergue et al., 2007).In total we measured EDC samples at 70 depth intervals, where 44 are single measurements and 26 represent the mean of 2 to 5 replicates from adjacent samples.To improve the robustness of the results and to increase overall temporal resolution, 13 depth intervals from the Talos Dome ice core were analysed, where 6 data points represent single measurements and 7 data points represent replicates.The Talos Dome samples cover a depth range between 1404 and 1444 m or a time interval from 126 000 to 154 000 yr BP on the TALDICE-1a gas age scale (where TALDICE-1a is synchronized to the EDC3 age scale via methane synchronization, Sch üpbach et al., 2011).At depths where replicate measurements were performed, the error bars in Fig. 1 represent the one σ standard deviation, while the mean reproducibility of the respective core was used for single measurements.The Introduction

Conclusions References
Tables Figures

Back Close
Full mean temporal resolution of our combined EDC and Talos Dome δ 13 C atm record is approximately 600 yr.

Method details
To release the air entrapped in Antarctic ice core samples, we used a sublimation extraction and an off-line coupled GC-IRMS system described in detail in Schmitt et al. (2011).The CO 2 mixing ratio was deduced by measuring the ratio of the total air volume released during the sublimation, and the CO 2 signal from the mass spectrometer.
The mean standard deviation of the CO 2 and δ 13 C atm for the Talos Dome samples measured in this study are approximately 1 ppm and 0.04 ‰, respectively.The δ 13 C atm precision of our method for Talos Dome is more than a factor of 2 better than other methodologies (Leuenberger et al., 1992;Lourantou et al., 2010a;Smith et al., 1999) and similar as previously measured data on the same ice core with the same method during Termination I (Schmitt et al., 2012).Reproducibility of CO 2 and δ 13 C atm of EDC ice are 3 ppm and 0.08 ‰, which is 1 ppm and 0.02 ‰ less precise than measurements on bubble ice in the same ice core using the same method (Elsig et al., 2009;Schmitt et al., 2011Schmitt et al., , 2012)).The origin of the lower precision in the deeper EDC ice is unknown but cannot be explained by methodological reasons as standard gas measurements were of usual precisions of 0.06 ‰ or better.The higher air content in the Talos Dome ice core with respect to EDC caused by the lower elevation of the drill site can also not account quantitatively for this difference.A possible change in impurity concentration could potentially add small contributions from in situ production of CO 2 , which can not be ruled out completely in Antarctic ice cores (Ahn et al., 2012).Note also that at depths below 1200 m in the EDC and 1100 m in the Talos Dome ice core, most air is stored in clathrates (L üthi et al., 2010).For dry extraction methods this could potentially affect the measurement results.However, since we use sublimation with 100 % extraction efficiency, the storage in bubbles or clathrates should have no effect on our analysis.Introduction

Conclusions References
Tables Figures

Back Close
Full Gas fractionation during the formation of clathrates in the ice, as observed for different gas species may potentially also result in δ 13 C atm variations on scales of centimetres as reported for O 2 /N 2 ratios (Kawamura et al., 2007) and CO 2 (Bereiter et al., 2009;L üthi et al., 2010), i.e. of equal scale as the length of our δ 13 C atm samples.Similar to the variations in CO 2 concentrations observed in EDC ice, potential δ 13 C variations due to clathrate fractionation may not be completely smoothed out by diffusion in cold EDC ice within the depth interval 200 m below the bubble clathrate transition zone (BCTZ).On the other hand small-scale CO 2 variations disappear quite rapidly below the BCTZ at Talos Dome (L üthi et al., 2010) and, thus, also a more efficient diffusional equilibration of potential δ 13 C variations is to be expected in this core.
We also compared our data with published data for Termination II by Lourantou et al. (2010b), which where also measured on the EDC ice core.A LGM value of −6.56 ‰ measured by Lourantou et al. (2010a) was about 0.16 ‰ lower than our sublimation measurements with a mean δ 13 C value of −6.40 ‰ (Schmitt et al., 2012).The offset between the Schmitt et al. (2012) and the Lourantou et al. (2010a) data in EDC the bubble ice for this time interval was systematic and was attributed to any method specific systematic fractionation.In contrast, our comparatively constant δ 13 C atm level of −6.8 ± 0.06 ‰ for the PGM derived on this clathrated ice appears to be in agreement with the PGM level of −6.8±0.1 ‰ from Lourantou et al. (2010b) (Fig. 1b) without invoking a systematic offset.However, given the combined scatter of the Lourantou et al. (2010b) and our data for the PGM time interval, an analytical offset on the order of 0.1 ‰ might still be hidden in the measurement error.

Gravitational correction
During the enclosure process in the ice, the air undergoes fractionation due to gravitational settling in the firn column, which is only dependent on the mass difference between 13 C and 13 C.The δ 13 C atm data are corrected for this gravitational fractionation using the stable isotope ratio of atmospheric N 2 (δ 15 N).Since N 2 has a much longer Introduction

Conclusions References
Tables Figures

Back Close
Full atmospheric lifetime than the time intervals considered here, δ 15 N changes only reflect the change in gravitational fractionation in the firn column (Dreyfus, 2008).Since δ 15 N data are not available for all EDC depth intervals, we used a linear regression between the measured water isotopes of the ice and measured δ 15 N.For the EDC ice core the correlation coefficient between δD and δ 15 N is R 2 = 0.88.The correction for gravitational fractionation in δ 13 C atm on ice from EDC is between −0.35 and −0.55 ‰, while Talos Dome samples are corrected by −0.35 to −0.42 ‰.The error introduced by the δ 15 N measurement itself is negligible because δ 15 N can be measured with precisions of 0.01 ‰, which is much smaller than the reproducibility of δ 13 C atm .However, the linear regression between δD and δ 15 N introduces an error between 0.04 and 0.05 ‰ at interglacial and glacial climate conditions, respectively.This additional uncertainty is affecting the overall values of our data, thus the absolute accuracy of the complete data set but not single point precision.

Monte Carlo cubic spline Average
Following the mathematical procedure as described in detail in Schmitt et al. (2012), a Monte Carlo cubic spline Average (MCA) was calculated to extract the robust δ 13 C atm trend information.We calculated 4000 cubic splines for the measured data.The cut off frequency (Enting, 1987) of the individual splines was chosen to roughly eliminate periods below 375 yr.Accordingly, high frequency measurement noise in the data is highly damped.Furthermore, a bootstrapping method was performed to detect possible outliers in the data set as described in Schmitt et al. (2012).However based on our bootstrapping outlier criterion (Schmitt et al., 2012), not a single point of our sublimation measurements was incompatible with the long-term trend, underlining the high quality of our δ 13 C atm data obtained with the sublimation method.Introduction

Conclusions References
Tables Figures

Back Close
Full In Fig. 1a we present our highly resolved CO 2 record between 155 000 and 105 000 yr BP from the EDC and Talos Dome ice cores.The stable carbon isotopic composition δ 13 C atm is shown in Fig. 1b.Both, the CO 2 concentration and δ 13 C atm , are compared to the EDC data of Lourantou et al. (2010b).The δD record (a proxy for local temperature) from the EDC ice core (Jouzel et al., 2007) is shown on a second ordinate in Fig. 1a to illustrate leads and lags between atmospheric CO 2 and temperature.Note, however, that the uncertainty in the ice age/gas age difference (∆age) must be taken into consideration: for instance ∆ age is 4200 and 2000 yr (Parrenin et al., 2007) at the PGM and the MIS5.5, respectively and relative errors in ∆age are about 20 % (Buiron et al., 2011).The data can be divided into six time periods of distinct characteristics for CO 2 , δ 13 C atm and δD.The first interval I (the end of the penultimate glacial), from 155 000 to 140 000 yr BP is characterised by fairly constant atmospheric values around 195 ppm CO 2 and −6.8 ‰ for δ 13 C atm , which is about 0.4 ‰ lower than during the LGM (Schmitt et al., 2011) (Fig. 2).At 140 000 yr BP our new CO 2 and the EDC δD data show minima of 193 ppm and −440 ‰, respectively.
Interval II (the penultimate glacial/interglacial transition from 140 000 to 128 500 yr BP) is characterised by rises in δD and CO 2 .Starting at the PGM, CO 2 gradually increases to peak values of 286 ppm at 128 500 yr BP.The Talos Dome data show an 8 ppm higher maximum than EDC at that time.The maximum in CO 2 at the onset of MIS5.5 lags behind the δD maximum by 500 to 1500 yr, which is still on the same order of magnitude than the uncertainty of about 400 yr in ∆ age at 128 500 yr BP (Parrenin et al., 2007, Buiron et al., 2011).No significant difference in the more slowly changing δ 13 C atm evolution is seen between the two ice cores.According to our new data a maximum atmospheric δ 13 C atm decrease of about 0.2 to −7.0 ‰ is observed between the PGM and 131 000 yr BP.Note that the temporal resolution is rather low in Introduction

Conclusions References
Tables Figures

Back Close
Full this time interval, which compromises the robustness of the magnitude of the δ 13 C atm drop.At 131 000 yr BP δ 13 C atm starts to rise from −7.0 to −6.7 ‰ and continues to do so until 128 500 yr BP.To illustrate differences in the CO 2 records across Termination II data from Lourantou et al. (2010b), Petit et al. (1999) and from this study are plotted in Fig. 3. Before 128 500 yr BP our new CO 2 data agree well with the data published in Petit et al. (1999) Fischer et al. (1999) and Lourantou et al. (2010b).However, between 128 000 to 126 000 yr BP our data suggest very constant CO 2 concentrations of about 275 ppm, whereas Petit et al. (1999) and Lourantou et al. (2010b) report CO 2 fluctuations between 265 and 280 ppm.Small overall offsets between CO 2 records from clathrate ice of the same ice cores along the same time intervals are most probably due to the specific air extraction techniques used in different studies.Air release after mechanical destruction of the ice matrix as used by Petit et al. (1999) and Lourantou et al. (2010b) suffers from incomplete extraction as release of CO 2 from decomposing air clathrates is slower than for other gases such as O 2 and N 2 (L üthi et al., 2010).Moreover, larger clathrates may be depleted in CO 2 compared to smaller clathrates but are preferentially released during dry extraction (Bereiter et al., 2012).In contrast sublimation of the ice sample, as used in our study, ensures a 100 % effective extraction of the air and, therefore, an unbiased CO 2 reconstruction.
In the case of the CO 2 offset of 8 ppm around 128 000 yr BP of our measurements derived for Talos Dome and EDC (Fig. 3), both records were measured with the sublimation technique.Hence, other potential processes have to be considered to explain this difference: (i) dating inconsistencies between the records; (ii) in situ production of CO 2 ; (iii) differences in gravitational enrichment; and (iv) a different damping of fast CO 2 changes connected to the different gas enclosure characteristics.Sch üpbach Introduction

Conclusions References
Tables Figures

Back Close
Full i) between 129 000 and 126 000 yr BP, the broad shape of the early CO 2 peak from the Talos Dome matches the timing of the peak reconstructed from the EDC ice core, but can not explain the higher CO 2 concentrations.In situ production of CO 2 (process ii) in the Talos Dome ice core is unlikely to account for the offset since the δ 13 C atm data from both ice cores agree within their uncertainties for that time interval.A strong depletion in the measured δ 13 C would be expected (Francey et al., 1997) from CO 2 in situ production due to microbial activity.Only if the CO 2 in situ production results from inorganic, acid-carbonate chemical reactions the effect on δ 13 C could be smaller and may remain undetected in our δ 13 C atm data.However, so far there is no indication of such carbonate reactions in the Talos Dome ice core.Correcting both ice cores for gravitational settling and thus enrichment of CO 2 (process iii), the offset increases by 1 ppm.Thus, a damping of relative short-lived CO 2 changes (process iv), which is different in the two ice cores is likely to account for the offset.Accumulation rates are higher in the Talos Dome (Stenni et al., 2011)  Interval IV starts at 120 000 yr BP and lasts until 115 000 yr BP.This interval is characterised by rather constant CO 2 concentrations around 276 ppm, whereas EDC temperature decreases by 3 K (Jouzel et al., 2007) as derived from the decline of δD from −390 to −408 ‰ .At the same time, atmospheric δ 13 C atm drops from −6.55 by about −0.2 to −6.75 ‰ .
Finally at 115 000 yr BP, the beginning of interval V, CO 2 drops by 40 ppm in 10 000 yr to 240 ppm.During the CO 2 decline starting at 115 000 yr BP, δ 13 C atm increases again to −6.55 ‰ .Note that the Antarctic temperature decrease during interval IV leads the CO 2 decrease in interval V by up to 5000 yr on the EDC3 age scale as previously observed in the Vostok ice core (Fischer et al., 1999;Cuffey and Vimeux, 2001).This lag is significantly larger than the ∆ age uncertainty at this time.
Between 108 500 to 105 000 yr BP (interval VI) our new CO 2 and δ 13 C atm data cover the Antarctic Isotope Maxima (AIM) event number 24 (Capron et al., 2010).EDC temperature rises by 4 K, according to an increase of δD from −434 to −409 ‰ .At the same time, our CO 2 data show a clear peak of 20 ppm, which is reproduced in the EDML record (Bereiter et al., 2012).While we also find the pronounced CO 2 maximum at the onset of MIS5.5, the strong drop in δ 13 C atm is not visible in our data.We performed high resolution sampling on both, the EDC and Talos Dome ice cores, around the time interval of this large δ 13 C atm excursion but no trace of such a drop could be found (Fig. 1b).Such a fast change in δ 13 C atm and the connected changes in the carbon cycle have not been found in any other pre-industrial time interval previously studied for δ 13 C atm .
We discuss in the following whether such an anomaly in δ 13 C atm can potentially be explained by natural carbon cycle fluctuations.To put such high anomalies into perspective to variations potentially observed in the carbon cycle we briefly report on the human perturbation of the carbon cycle of the recent past.The unprecedented anthropogenic CO 2 emissions from fossil fuel burning and land use, which amount to more than 500 Gt C in the last 250 yr (Boden et al., 2011;Houghton, 2010), has led to a rise of more than 100 ppm in atmospheric CO 2 (Etheridge et al., 1996;Keeling et al., 2009) and a δ 13 C atm decrease of −1.9 ‰ from pre-industrial conditions of −6.3 to −8.2 ‰ in 2010 AD (Francey et al., 1999;Keeling et al., 2010).Thus, a 0.7 ‰ drop in 200 yr during Termination II accompanied by a CO 2 increase of not more than 20 ppm represents an extraordinary perturbation.Moreover, the true atmospheric signal in CO 2 and δ 13 C atm was even larger than the signal recorded in the ice core.The amplitude of a CO 2 peak recorded in the ice core depends on the age distribution during gas enclosure and the rate of CO 2 emissions.öAccumulation rates in Antarctica during MIS5.5 were higher than at present day (Parrenin et al., 2007), therefore the gas age distribution during gas enclosure was narrower at that time than for recent times.For the Holocene the mean width of this age distribution was calculate to be about 200 yr (Joos and Spahni, 2008).In the EDC ice Introduction

Conclusions References
Tables Figures

Back Close
Full core this leads to a damping to about 40 % of the true atmospheric δ 13 C atm variations for periodicities of 600 yr (K öhler et al., 2010).Thus, neglecting any effect of uncertain emission rates we estimate that the true atmospheric signal underlying the data of Lourantou et al. (2010b) around 128 000 yr BP would have been as much as −1.6 ‰ .If such a negative anomaly in δ 13 C atm represents a true climatic signal, one has to assume drastic environmental changes such as a sea-level change during Termination II leading to a rapid sea-level rise of 20-40 m within 1000 yr and a release of carbon formerly stored on coastal areas (Siddall et al., 2006;Kopp et al., 2009).A similar flooding hypothesis has recently been proposed to explain the rapid jump in CO 2 around 14 600 yr BP during melt-water pulse 1A (K öhler et al., 2011).Such a scenario would release terrestrial carbon with isotopic light 13 C in line with a large negative anomaly in δ 13 C atm .However, a box modelling approach (K öhler et al., 2011), shows that a true atmospheric anomaly in δ 13 C atm of −0.87, −1.26 or −1.62 ‰, caused by the injection of 200, 300 or 400 Gt C within 100 yr would also lead to an atmospheric CO 2 rise of 47, 73 and 99 ppm, respectively, which is not seen in our record.
The recovery in δ 13 C atm between 127 900 and 127 700 yr BP from −7.55 to −6.8 ‰, i.e. 0.75 ‰ in 200 yr, implies even higher rates of change, and is even more difficult to explain on the basis of natural carbon cycle fluctuations.
To explain the dramatic excursion in the Lourantou et al. (2010b) data, we must rather consider (i) unidentified analytical problems (e.g.contamination) or (ii) fractionation due to an incomplete gas release during the dry extraction.In summary, we conclude that rapid variations in δ 13 C atm within a few centuries appear to be unlikely to have occurred at any time in our record in line with the rather slow evolution of CO 2 concentration in the same time interval.

Comparison of sequence of events around Termination II and Termination I
From our new data we calculate a δ 13 C atm Monte Carlo Average (MCA) using the method described in Schmitt et al. (2012) to illustrate the significant millennial changes Introduction

Conclusions References
Tables Figures

Back Close
Full in δ 13 C atm .The MCAs for the last 24 000 yr and for 155 000-105 000 yr BP are shown in Fig. 2. To compare Termination II with Termination I, the x-axes are aligned, such that the first CO 2 maxima at the end of Termination I and II match.Thus, the LGM time slice, corresponding to 22 000 yr BP (Clark et al., 2009), coincides with the PGM at 140 000 yr BP.Also illustrated in the same figure is the atmospheric CO 2 and EDC δD evolution (Jouzel et al., 2007) within the same time frames.Within the uncertainty, δ 13 C atm is roughly constant around −6.8 ‰ during the PGM between 155 000 and 140 000 yr BP.This is consistent with the carbon cycle being in dynamic equilibrium during both the LGM (Schmitt et al., 2012) and the PGM, however the average PGM δ 13 C atm value is 0.4 ‰ higher than for the LGM (see Sect. 4.4 for further discussion of this difference).The δ 13 C atm evolution is characterised by a sharp synchronous incline at the start of the CO 2 rise in Termination I (Fig. 2).From 17 000 to 16 000 yr BP a strong δ 13 C atm decrease to −6.7 ‰ is observed, which is accompanied by rising CO 2 concentrations, warmer Antarctic temperatures and a strong drop in aeolian dust deposition in the Southern Ocean region.Note that dust fluxes are higher during the LGM than at the onset of the CO 2 rise (Fig. 2).
During Termination II the stable carbon isotope signature decreases by about 0.2 ‰ from about 136 000 to 130 000 yr BP.As temporal resolution as well as precision of the δ 13 C atm data are comparably low in that time interval it is hard to tell when exactly δ 13 C atm starts to drop and how large the maximum drop is.However, we can say that the major increase in CO 2 (and δD) starts at 136 000 yr BP, while the first strong decrease in aeolian dust deposition (Fig. 2) occurs already at 138 000 yr BP (Lambert et al., 2012).Also it becomes clear that the overall magnitudes of the δ  Farell andPrell, 1989, Hodell et al., 2011).This region is exclusively flushed with low alkalinity water originating from Circumpolar Deep Water (Hodell et al., 2001(Hodell et al., , 2003;;Rickaby et al., 2010).The Atlantic Basin style is represented by ODP1090 drilled at a depth of 3700 m, where northern sourced waters of high alkalinity control the preservation history (Hodell et al., 2001(Hodell et al., , 2003;;Rickaby et al., 2010).Hence, ODP1090 is mainly influenced by changes in ocean circulation strength and stratification of North Atlantic Deep Water (NADW), whereas preservation and dissolution events characteristic for the deep Indo-Pacific Ocean are recorded in ODP1089.
The CaCO 3 preservation histories of ODP1089 and ODP1090 are also consistent with carbonate ion reconstructions from Broecker et al. (1999) and Barker and Elderfield (2002).Nevertheless, trends are only qualitatively covered, as changes in weight percent do not linearly correlate with the sedimentation rate.

CPD Introduction Conclusions References
Tables Figures

Back Close
Full The CaCO 3 preservation histories of the Atlantic and the Pacific Oceans show opposing behaviour during interglacials and the glacial inception but the same trend during glacials and the deglaciations.Increasing carbonate ion concentration in the deep ocean due to ocean reorganization within a glacial termination may have triggered a deepening of the carbonate compensation depth by about 700 m (Farrel and Prell, 1989;Broecker et al., 1993), which allows for enhanced CaCO 3 preservation in the Atlantic and the Pacific Oceans.A strengthening of deep ocean ventilation during the termination is also supported by deep sea δ 13 C reconstructions of benthic sediment cores (Oliver et al., 2010).Those show increasing deep ocean δ 13 C values throughout the LGM and Termination I, indicating a decrease of the δ 13 C surface to deep ocean gradient (Hodell et al., 2003).
While an intensified deep ocean circulation at the start of the deglaciation helps to explain the rising CaCO 3 preservation during the terminations, it falls short to explain the fact that preservation increased already 6000 yr before the deglaciation starts.In fact the preservation records (Atlantic and Pacific) are already half way up from their minimum values when the deglaciation starts, pointing to internal processes of the alkalinity budget of the ocean rather than circulation to be responsible for this longterm trend in the preservation history.
Note also that the preservation state in the deep ocean was already on a higher level during the PGM than for the LGM, while at the end of Termination II the magnitude of preservation was the same as after Termination I. Additionally, deep ocean δ 13 C DIC (Oliver et al., 2010) started to rise at least 10 000 yr before the PGM, while δ 13 C DIC in the deep ocean decreased for 10 000 yr before the LGM.
As shown by transient simulations over the past 130 000 yr with a 3-d dynamic ocean model, a multitude of physical and biogeochemical processes influences the spatio-termporal evolution of carbonate preservation, δ 13 C of dissolved inorganic carbon, atmospheric CO 2 , δ 13 C atm , and other proxies on millennial and even longer time scales.This is a consequence of ocean-sediment interactions and the ∼ 100 000 yr time scales associated with the weathering-burial cycle of organic carbon, opal, and Introduction

Conclusions References
Tables Figures

Back Close
Full The major δ 13 C atm increases in Fig. 2 starting at around 130 000 and 12 000 yr BP, respectively, can be largely explained by the regrowth of the terrestrial biosphere (K öhler et al., 2005;K öhler and Fischer, 2006;Elsig et al., 2009;Kleinen et al., 2010;Menviel and Joos, 2011) as well as warming SST, where a SST change of +1 K induces a δ 13 C atm fractionation by about +0.1 ‰ (Zhang et al., 1995).The δ 13 C atm signal induced from terrestrial carbon build-up is a mixture of a signal composed of C4 (δ 13 C ∼ −10 to −14 ‰) and C3 (δ 13 C ∼ −22 to −30 ‰) plant type origin (O' Leary, 1981, Ehleringer et al., 1993).As CO 2 and temperature increase during glacial terminations, many terrestrial ecosystems store more carbon on land.Additionally, large areas once covered by ice sheets transform to boreal forests and peat lands with high carbon stocks.This net gain in terrestrial carbon is assumed to be partially reduced by carbon release due to shrinking permafrost extent (Zimov et al., 2009).
The CO 2 evolution at the end of the terminations is similar in that the maximum values are reached already at the very beginning of the interglacial.We attribute the early CO 2 maximum mainly to the opposing effects of the ongoing process of carbonate compensation in response to the net glacial/interglacial transfer of carbon from the ocean/atmosphere system to the terrestrial biosphere, as well as the continuing regrowth of the terrestrial biosphere during the early interglacial (Fischer et al., 1999, Inderm ühle et al., 1999).Note however that CO 2 variations result from a subtle balance of forcings and processes acting on different time scales and partly in opposite direction as well as from memory effects associated with changes occurring during the termination (Menviel and Joos, 2011).Different interglacial periods with different forcing histories are thus expected to yield different CO 2 and δ 13 C atm evolutions as documented by our measurements.In this respect, Termination I is special since it is characterised by the intermittent cooling during the YD event, which was also accompanied by an intermittent decline in terrestrial carbon storage (K öhler et al., 2005).This

Conclusions References
Tables Figures

Back Close
Full late anomaly in terrestrial carbon storage during Termination I would also lead to a delay of the carbonate compensation signal in CO 2 in the Holocene compared to MIS5.5.Accordingly, we attribute the increasing CO 2 trend later in the Holocene and the lack of such a trend during MIS 5.5 to the effect of carbonate compensation in response to the occurrence of the YD cold anomaly during Termination I.
Both during Termination I and the subsequent early Holocene on the one side and Termination II plus the early MIS5.5 on the other side, δ 13   marine records used for this correction are given on their individual age scales, which are also different from the EDC3/TALDICE-1a ages scale used for our ice core records.This could lead to erroneous δ 13 C atm changes in the SST corrected version on a millennial time scale.Absolute magnitudes and evolutions in our SST corrected δ 13 C atm record highly depend on past temperature reconstructions, which suffer from regional differences and potential dating inconsistencies, as well as model specific assumptions and calculation limitations.Accordingly, the temperature correction is shown only to illustrate the physical effect of changing SST on δ 13 C atm and we refrain from quantitatively interpreting the millennial variations in the SST corrected δ 13 C atm evolution in Fig. 2.
In summary, taking into account CO 2 and δ 13 C atm data, essentially the same processes have been active during Termination II as in Termination I but with different strengths, different relative timing and from different starting conditions.This sequence of events is comprised of (i) the up-welling of isotopically light water from the deep ocean starting at the onset of the termination accompanied by a decline in aeolian Fe fertilization; (ii) the increase in SST leading to higher CO 2 and δ 13 C atm ; (iii) the regrowth of the terrestrial biosphere starting in the second half of the transition, which ceases only a few thousand years into the following interglacial; (iv) the long-term effect of carbonate compensation in response to the transfer of carbon into the terrestrial biosphere.To explain the differences in the temporal evolution of carbon cycle changes during both Terminations and the subsequent interglacials the different temperature evolution during the terminations (such as the lack of a Younger Dryas cold spell in Termination II) have to be considered.Furthermore, other processes, such as changes in marine productivity or in the re-mineralisation depth of organic matter, nutrient release from shelves, shallow water carbonate compensation or changes in ocean salinity, have likely contributed to the overall evolutions of atmospheric CO 2 and δ 13 C.
A quantitative attribution of carbon cycle changes during Termination I and II cannot be achieved from our data only, however, our new δ 13 C atm record provides crucial data constraints for future carbon cycle model experiments.Introduction

Conclusions References
Tables Figures

Back Close
Full At 120 000 yr BP, Antarctic temperature begins to decline in Fig. 2 (Jouzel et al., 2007), whereas atmospheric CO 2 concentration remains at a constant value of 276 ppm for another 5000 yr.One may expect that the net uptake of CO 2 by the ocean should follow Antarctic temperature based on the relatively high correlation between these two parameter over the past 800 000 yr (Fischer et al., 2010).Given the near constancy of atmospheric CO 2 during the early inception, individual processes that increase or lower atmospheric CO 2 must balance each other.Fischer et al. (1999) attributed the delayed response of CO 2 mainly to the terrestrial release of carbon into the atmosphere, which counterbalances a transfer of CO 2 to the deep ocean.In contrast Hodell et al. (2001) attributed the continuously high CO 2 levels to a delayed release of carbon from the ocean due to carbonate compensation induced from terrestrial biosphere build-up and ocean reorganization during Termination II and the early MIS5.5.The main difference between the two theories is that a release from biomass requires a net change in δ 13 C atm , while the carbon isotopic signature of CO 2 released by carbonate compensation hardly affects atmospheric δ 13 C atm , as the released carbon originates from ocean DIC with δ 13 C between +2 and +3 ‰, which is similar to surface waters in equilibrium with δ 13 C atm .Note however, that besides these two, other processes, e.g.changes in ocean circulation or marine productivity act simultaneously and that a quantitative understanding demands a complete carbon cycle consideration.
Our new data show δ 13 C atm to be declining by 0.2 ‰ between 120 000 and 115 000 yr BP, thus taken at face value supporting the biosphere hypothesis.However, as SST decreases in parallel the solubility pump drags down CO 2 and also induces an atmospheric fractionation to 13 C depleted values; δ 13 C atm is reduced by 0.1 ‰ as mean SST cools by 1 K (Mook, 1986, Zhang et al., 1995), hence a 0.2 ‰ decrease in δ 13 C atm as observed in our record would require a mean SST cooling by 2 K at that time.MARGO Project Members (2011) reported a global mean SST increase of only 1.9 ± 1.8 K, from Introduction

Conclusions References
Tables Figures

Back Close
Full the LGM to the Holocene; thus, it is unlikely that the 0.2 ‰ drop in δ 13 C atm can be completely attributed to the much more subtle cooling between 120 000 and 115 000 yr BP.Moreover, the solubility effect would imply a reduction in atmospheric CO 2 concentrations as well, which is not observed.Assuming that the entire δ 13 C atm effect is explained by the solubility effect, the carbon transfer from the atmosphere to the surface ocean, has to be compensated by some other process as atmospheric CO 2 concentrations remain unchanged.This could be due to carbon release from the terrestrial biosphere, which favourable would add to the δ 13 C atm decrease, thus, would relax the required SST cooling constraint.In summary, three mechanisms have to be considered: Cooling SST draws down both, atmospheric CO 2 and δ 13 C atm .As the terrestrial biosphere expanded during Termination II, it is reasonable to assume that the terrestrial carbon storage shrinks to some degree when temperatures decrease at the glacial inception.This release of terrestrial carbon induces further δ 13 C atm lowering.In total the reduction of CO 2 due to the solubility effect has to be balanced by the release of terrestrial carbon into the atmosphere and/or carbon release from carbonate compensation while at the same time closing the δ 13 C budget.At the moment we cannot quantitatively decipher how much each of the processes discussed above contributed.A precise and ice core synchronized global SST reconstruction for this time period will be essential to better constrain the SST effect on δ 13 C atm .Also global higher resolution δ 13 C compilations for the deep and surface ocean over this time period could help to answer this question as a release of terrestrial carbon leads to a net decrease in δ 13 C in the ocean/atmosphere system.

Understanding the 0.4 ‰ offset in δ 13 C atm between the PGM and LGM
δ 13 C atm during the LGM was on average about 0.4 ‰ heavier than during the PGM.At the same time an enrichment of 0.2 to 0.4 ‰ for the mean ocean δ 13 C DIC change is recorded in many marine sediment cores (e.g.Oliver et al., 2010;Hoogakker et al., 2007).The change of the combined atmosphere/ocean system to 13 C enriched values Introduction

Conclusions References
Tables Figures

Back Close
Full suggests an overall shift of the isotopic signature of the complete carbon stock.An overall δ 13 C shift of the combined ocean plus atmosphere carbon storage has to be induced by fluctuations in the carbon fluxes entering and leaving the ocean/atmosphere system.In contrast such an overall shift cannot be explained by ocean ventilation changes, which would enrich the surface ocean and deplete the deep ocean in δ 13 C DIC .
To discuss an overall long term δ 13 C drift of the ocean/atmosphere system at orbital time scales the carbon fluxes into and out of the ocean/atmosphere system have to be considered.The main ingoing carbon fluxes are volcanic CO 2 emission at subduction zones and mid ocean ridges, (0.07-0.08 Gt C yr −1 , Sano and Williams, 1996), carbonate weathering on the continents (0.1 to 0.2 Gt C yr −1 , Archer et al., 1998;Zeebe and Westbroek, 2003) and mineralization of old sedimentary organic carbon on continents and shelves (0.06-0.12 Gt yr −1 , Zeebe and Caldeira, 2008).The main outgoing carbon fluxes are burial of carbonates in marine sediments and coral reefs as well as burial of organic matter on the ocean floor and on the continents.An important example of such partially buried carbon reservoirs on continents is the organic carbon stored in the deep permafrost and yedoma in high northern latitudes.The estimated amount of organic carbon in this region is large (1700 Gt C, Tarnocai et al., 2009), yet it is unclear when and at what rates this carbon was accumulated in the past and under what circumstances this carbon becomes remineralized and reintroduced into the ocean/atmosphere reservoir.Since all the carbon fluxes mentioned above are small compared to the annual exchange fluxes between biosphere/atmosphere and atmosphere/ocean, they lead to an estimated mean residence time of carbon in the ocean/atmosphere system of around 100 000 to 150 000 yr.In other words, after one glacial cycle 50-70 % of the carbon is replaced by new carbon.On the one hand this effect is able to produce a drift in the mean δ 13 C of the atmosphere/ocean system, on the other hand an imprinted δ 13 C signal (for instance caused by a growth in the terrestrial biosphere carbon storage) becomes gradually lost on these time scales.A long term δ 13 C drift of the ocean/atmosphere system can result either from a change in these fluxes or in their isotopic δ 13 C signature or a combination of both.

CPD Introduction Conclusions References
Tables Figures

Back Close
Full In the following we want to discuss the following alternatives to explain the observed 0.4 ‰ offset in δ 13 C atm between the PGM and LGM: scenario (1) invokes differences in carbon fluxes and the corresponding isotopic signatures by continental weathering and sedimentation at the ocean floor; scenario (2) invokes changes in volcanic outgassing of CO 2 ; and scenario (3) a net shift in the terrestrial carbon storage between the PGM and LGM (Lourantou et al., 2010b).Zeebe and Caldeira (2008) examined the effect of a long-term imbalance in carbon input and output fluxes on atmospheric CO 2 .They studied the hypothesis that carbon input into the ocean/atmosphere system by volcanism and net organic matter oxidation is balanced by silicate rock weathering and subsequent carbonate mineral burial.They argue that the mean decrease in atmospheric CO 2 of 22 ppm during the last 650 000 yr is due to an imbalance between carbon addition and removal with a net flux removing carbon from the ocean/atmosphere system.Interpolating over the last 150 000 yr, this might account for the 5 to 10 ppm concentration decline between the PGM and LGM.Stable carbon isotope signatures were not considered yet in the study by Zeebe and Caldeira (2008).
Carbonate weathering causes a carbon flux entering the atmosphere-ocean reservoir of about 0.1 to 0.2 Gt C possibly affecting δ 13 C.Note that even if carbon fluxes are totally balanced a change in the isotopic signature of the ocean/atmosphere system is possible, if ingoing and outgoing carbon fluxes do not carry the same isotopic signature.E.g. the source could change as different rocks are weathered and with changing sources also the carbon isotopic composition of the carbonate could be altered (scenario (1)).Sedimentary carbonate rocks that are subject to continental weathering carry carbon isotopic signatures between 2 and 3 ‰ (Rubinson and Clayton, 1969, Emrich et al., 1970, Romanek et al., 1992), which would allow for a long-term trend in δ 13 C atm .Unfortunately, there is not sufficient observational evidence available to reconstruct a change in the δ 13 C signature of the carbonate weathering.
The second possible scenario (2) concentrates on the volcanic carbon release.According to Sano and Williams (1996) volcanic outgassing of CO 2 releases about 0.07 to

Conclusions References
Tables Figures

Back Close
Full 0.08 Gt C yr −1 .These inputs are assumed to be compensated by sedimentary carbonate burial in the ocean over time scales of thousands of years, which is controlled by the adaptation of the carbonate compensation depth in the ocean in response to changing riverine carbonate and alkalinity fluxes.Changing the ratio of the volcanic carbon flux via subduction along the continental crust relative to the release from mid-ocean-ridge basalt may shift the mean δ 13 C signature of the ocean-atmosphere carbon stock.Carbon released from mid-ocean-ridge regions (δ 13 C ∼ −6.5 ‰) shows 13 C depleted CO 2 values relative to subduction areas by −1 to −2 ‰ (Sano and Williams, 1996).Alternatively, a net global reduction of the volcanic carbon release results in enrichment of the mean δ 13 C input signal as the isotopic signature of the riverine input becomes more important.Hence, the PGM to LGM offset could also be induced by changed volcanic emissions.Roth and Joos (2011) attribute only minor δ 13 C atm changes of 0.06 ‰ to an enhanced volcanic carbon release during Termination I.However, on time scales of a complete glacial cycle continuing changes in emission rates or isotopic signature of volcanic carbon fluxes might potentially alter δ 13 C atm more strongly and induce long-term drifts.Note however, that continuous changes in emission rate require other processes balancing CO 2 concentrations.
In the third scenario (3) a larger terrestrial carbon pool at the LGM compared to the PGM could cause a net shift in δ 13 C.This could be due to a net change in the mid and low latitude biomass or, more likely, by carbon locked away in permafrost soils.
Studies on northern permafrost regions concerning amount and age of stored organic carbon suggested permafrost soil pool sizes of up to 1700 Gt C (Tarnocai et al., 2009), with 88 % being perennially frozen.Permafrost in the central Yukon region in Canada revealed ages older than 740 000 yr BP (Froese et al., 2008), showing that a long-term accumulation of peat is possible.Also Reyes et al. (2010) claimed the persistence of deep discontinuous permafrost during the last 200 000 yr.Such organic carbon pools represent a possibility to remove carbon from the ocean-atmosphere system on orbital time scales.

Conclusions References
Tables Figures

Back Close
Full A simple mass balance approach is performed in Eqs. ( 1) and ( 2) to quantify the necessary amount of carbon that has to be removed from the ocean-atmosphere carbon stock by organic carbon in permafrost to alter atmospheric δ 13 C atm by ∆δ atm = 0.4 ‰ as observed in our study and mean ocean δ 13 C by ∆δ DIC = 0.2 ‰ to 0.4 ‰ (Oliver et al., 2010).Note, that we allow different values for ∆δ atm and ∆δ DIC as internal ocean-atmosphere changes can potentially alter the δ 13 C difference between the ocean and atmosphere reservoirs.
(∆δ ocean + δ PGM ocean ) Changes of Total Organic Carbon, rain ratio and marine production are not considered in the following, essentially assuming that these processes acted in a comparable way both during the last and the penultimate glacial.Furthermore, the isotopic signature of formed biomass ∆C terr is assumed to remain constant at δ terr = −28 ‰, representative of the C3 plant metabolism in high latitudes (Novak et al., 1999;Dioumaeva et al., 2003).Hence this calculation represents a minimum flux scenario for a terrestrial carbon source of distinct isotopic signature.To quantify the atmospheric carbon reservoir during the PGM C PGM atm , as well as the change between PGM and LGM ∆C atm , atmospheric CO 2 concentration is converted to GtC.One ppm corresponds to 2.123 Gt C.
The size of the LGM ocean carbon stock C PGM ocean is set to 39 400 Gt C (K öhler et al., 2005).Our new δ 13 C atm data show stable carbon isotope values of δ PGM atm = −6.8‰ at the PGM.The ocean value is set to δ PGM ocean = −0.5 ‰ (Oliver et al., 2010).Using this mass balance approach we can quantify the amount of carbon required to explain a 0.4 ‰ offset.However, in this mass balance scenario we assume an instantaneous shift from the PGM to the LGM time window.In reality both time intervals are about 120 000 yr apart.Accordingly, part of the carbon anomaly will be removed over time by the continuous replacement of carbon in the biosphere/ocean/atmosphere 2039 Introduction

Conclusions References
Tables Figures

Back Close
Full system by the carbonate weathering and volcanic carbon input flux on the one hand and sediment burial on the other.Summing up these fluxes between 0.274 to 0.344 Gt C enter and leave the ocean-atmosphere system during one year and the complete ocean-atmosphere carbon pool is recycled once in less than 150 000 yr. Thus, after 120 000 yr only a small fraction of the δ 13 C anomaly is preserved.Approximations using the Bern3D climate model of intermediate complexity lead to a preservation of 25 % of any δ 13 C anomaly after 100 000 yr (Schneiter, 2011).For a maximum estimate we assume that this isotopic dilution process is switched off, i.e. 100 % of the anomaly are preserved.
Results of the mass balance above, which account for the isotopic dilution by the sedimentary turn-over of carbon, are shown in Fig. 4. The green curve in the centre is calculated for changes in oceanic ∆δ ocean = 0.3 ‰, whereas the black envelopes result from ∆δ ocean = 0.2 ‰ and 0.4 ‰, respectively.The errors for the isotopic value of the deep ocean δ PGM ocean and the ocean carbon storage C PGM ocean are negligible compared to the unknown rate of signal dilution.
In the case of 25 % signal perservation, 1000 to 2000 Gt C have to be transferred from the ocean-atmosphere system to the terrestrial carbon pool to explain the 0.4 ‰ offset in δ 13 C atm between LGM and PGM as shown in away such terrestrial carbon on time scales of 10 000 yr and more.From our new data we also conclude that the major part of such a possible increase in terrestrial carbon storage did not occur during MIS5.5 as δ 13 C atm still was 0.2 ‰ lower at 100 000 yr BP compared to the LGM (0.4 ‰ in the SST corrected records).Hence, the major part of such as net biomass accumulation must have occurred after MIS5.5 but before the LGM.
In summary, our new data confirm the PGM to LGM shift as reported in Lourantou et al. (2010b), but show a larger magnitude.Based on current knowledge it is not possible to unambiguously exclude one of the three explanation approaches of (1) changes in carbon isotope fluxes of carbonate weathering and sedimentation at the sea floor, (2) varying contributions from volcanic outgassing or (3) glacial peat build-up on land that would be conserved till present days.In contrast the required net growth of terrestrial biomass outside of permafrost regions is probably too large to explain the offset.

Conclusions
Using ice core measurements of atmospheric δ 13 C atm , we are able to better constrain the possible pathways of the global carbon cycle evolution in the past.First order interpretation of our new data for the penultimate glacial/interglacial transition shows that the same processes influence the global carbon cycle at that time as already reported for the last 24 000 yr, however with different magnitude, timing and starting conditions.
The early MIS5.5 is dominated by a build-up of land biosphere and increased preservation of calcium carbonate in the ocean.In this respect it is interesting to note that the different temporal evolution of CO 2 in the Holocene and in MIS5.5 does not require an extraordinary anthropogenic source of carbon into the atmosphere (Ruddiman et al., 2003), but can be explained by the different temporal evolution of the carbon cycle the long-term memory of the carbon cycle through carbonate compensation have to be considered to explain the increasing CO 2 levels during the Holocene (Elsig et al., 2009), while such an increase is not observed during MIS5.5.The appearance of a minimum in δ 13 C atm during the late MIS5.5, likely induced by a terrestrial carbon release and colder SST, can explain stable CO 2 levels for many thousand years, when Antarctic temperatures already declined.However, the amount of carbon released from the biosphere at that time cannot be quantified with certainty, since carbonate compensation may have also contributed to CO 2 emissions into the atmosphere.
A reasons for an observed offset between the mean values in δ 13 C atm of the current and last interglacials cannot unambiguously be assigned.A build-up of permafrost by 1000 to 2000 Gt C with a C3 plant isotopic signature is one possible scenario to explain the observed shifts between LGM and PGM in atmospheric CO 2 and δ 13 C atm , while changes in terrestrial carbon stocks other than in permafrost are unlikely to account for the observed δ 13 C change.Furthermore, a temporal change in the isotopic ratio of carbonate weathering between the end of the previous interglacial and the last glacial maximum or fluctuations in the contributions from volcanic input fluxes may provide alternative explanations.A future continuous record of δ 13 C changes over the entire last glacial cycle connecting the MIS5.5 with the LGM will be key to better understand these long-term changes in carbon (isotope) imbalances.Furthermore extensive three dimensional dynamical ocean modelling on a regional and global scale including all relevant carbon cycle processes and dynamic input/output fluxes of carbonate and silicate weathering, volcanic emission and CaCO 3 fluxes into the sediments, are crucial as carbon cycle fluxes are not decoupled from each other.Introduction

Conclusions References
Tables Figures

Back Close
Full Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | than the EDC(Parrenin et al., 2007)   ice core.Therefore fast fluctuations in atmospheric CO 2 are less damped in the Talos Dome than in the EDC ice core.Using recent (interglacial) bubble enclosure characteristics for EDC and Talos Dome (Joos andSpahni, 2008) a CO 2 variation of 20 ppm with a wavelength of 1000 yr (similar to the CO 2 change observed in the Talos Dome ice core at the beginning of MIS 5.5) would be damped to 13 ppm in the EDC ice core in line with our observations.In summary, a combination of both an age scale problem and different damping characteristics of the analysed ice cores can explain the CO 2 mismatch during this time interval.At the onset of time interval III, 128 500 yr BP, CO 2 concentrations as well as EDC δD slightly decline to 276 ppm and −387 ‰, respectively.From 126 000 yr BP onwards, CO 2 concentrations as well as δD remain constant at about 276 ppm and −387 ‰Discussion Paper | Discussion Paper | Discussion Paper | rapid δ 13 C atm variations previously found during Termination II In Fig. 1 our new δ 13 C atm and CO 2 data are compared to published measurements from the same EDC ice core covering Termination II (Lourantou et al., 2010b).CO 2 measurements of both studies basically show the same picture: fairly constant CO 2 concentrations during the PGM, which strongly rose during the Termination and ended with CO 2 peak concentrations of 290 ppm at 128 500 yr BP.However, δ 13 C atm measurements from Lourantou et al. (2010b) show deviations of −0.8 to +0.2 ‰ relative to our δ 13 C atm record.The most striking difference to the MCA is observed at 128 000 yr BP.Lourantou et al. (2010b) found a sharp negative excursion in δ 13 C atm of about Discussion Paper | Discussion Paper | Discussion Paper | −0.66 ‰ in 1000 yr, which is in concert with the early CO 2 peak at 128 500 yr BP.This excursion in their data set is represented by three data points at 128 500, 128 200 and 127 900 yr BP with δ 13 C atm values of −7.25, −7.4 and −7.55 ‰, and CO 2 concentrations of 290, 290 and 287 ppm, respectively.At 127 700 yr BP the stable carbon isotopic signal rises again to −6.8 ‰ in δ 13 C atm and 277 ppm in CO 2 within only 200 yr.
Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | 13 C atm decline at the start of Termination I and II are slightly different: After the LGM δ 13 C atm drops by about 0.3 ‰, whereas along Termination II a maximum decline of 0.2 ‰ is found.Interestingly the terminations are also characterised by generally enhanced preservation of CaCO 3 as reconstructed from the sediment cores ODP1089 and ODP1090, drilled in the Cape Basin, located in the Atlantic sector of the Southern Ocean (Hodell Discussion Paper | Discussion Paper | Discussion Paper | et al., 2001) (Fig. 2).ODP1089 was drilled at 4600 m depth, recording the Pacific Basin style CaCO 3 development ( ,Anderson et al. (2009),Tschumi et al. (2011) andSchmitt et al. (2012), a strengthening of the Southern Ocean mixing, transporting isotopically light carbon from the deep ocean to the surface, was most likely responsible for the decreasing δ 13 C atm at the start of the deglaciation in accordance with the CO 2 and CaCO 3 preservation increases.Most likely, this carbon had been previously accumulated in the deep ocean due to decreased ventilation of Antarctic Bottom Waters during glacial times.Also a decrease in iron fertilisation of the surface productivity in the Southern Ocean synchronous to the decline in aeolian dust deposition at the onset of the termination may have contributed to the early δ 13 C atm decrease, however model experiments show that the change in the biological pump due to increasing iron limitation at that time can only explain about half of the δ 13 C atm decline(K öhler et al., 2005, 2010)  observed in our δ 13 C atm record.Note also in Fig.2that CO 2 increases significantly only after dust fluxes have already declined from the glacial maximum values to about 4-6 mg m −2 yr at the start of the termination.
Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | calcium carbonate.A fully consistent interpretation of the proxy records over Termination I and II would thus require knowledge of the previous evolution of the Earth System.
Discussion Paper | Discussion Paper | Discussion Paper | δ 13 C atm corrections are derived from approximated SST evolutions (shown on a second ordinate) which represent the area weighted mean from the individual changes in the Equatorial Pacific, Equatorial Atlantic, Southern Ocean, North Pacific and North Atlantic Ocean boxes (K öhler and Fischer, 2006).The SST corrected δ 13 C atm evolution along MIS5.5 supports the theory of an additional terrestrial build-up until 123 000 yr BP, as a net δ 13 C rise of about 0.1 to 0.2 ‰ remains after correcting for SST changes.Note that the Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | 4.3 Lag of the CO 2 decrease behind the temperature drop during the last glacial inception Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper |

Fig. 4 .
Such a range may still be realistic in the framework of the permafrost reconstructions byFroese et al. (2008), Tarnocai et al. (2009) and Reyes et al. (2010).In contrast if we assume a change in biomass in low and middle latitudes with a typical δ 13 C terr signature of −22 ‰ for the LGM(Francois et al., 1998), thus significantly heavier than carbon locked in permafrost, the total amount of additional stored terrestrial carbon required to explain the δ 13 C anomaly in the PGM for 25 % preservation increases to more than 2000 Gt C (not shown).This is on the same order as the total carbon stored in living biomass and soils today and is approximately a factor of 4 larger than the net change in the terrestrial carbon storage from the LGM to the Holocene.Accordingly, such an increase in nonpermafrost terrestrial carbon is most likely too large and unlikely to be able to account for the δ 13 C anomaly.Hence, permafrost represent the most likely candidate to lock FiguresBack Close FullScreen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | during the preceding deglaciation.The unique Younger Dryas climate event during Termination I, which has no comparable counterpart during Termination II, together with Discussion Paper | Discussion Paper | Discussion Paper | Screen / Esc Printer-friendly Version Interactive Discussion Discussion Paper | Discussion Paper | Discussion Paper | Discussion Paper | IPEV and PNRA (at Dome C) and AWI (at Dronning Maud Land).This work is also a contribution to the EU Integrated Programme "Past4Future: Climate change -Learning from the past climat" within the European Union's Seventh Framework programme (FP7/2007-2013) under grant agreement no.243908.We thank A. Landais, D. Rodriguez, E. Capron and G. Dreyfus for the contribution of δ 15 N data.We acknowledge financial support by the Swiss NSF, the DFG priority program INTERDYNAMIK and the German climate program DEKLIM.We thank J. Chappellaz for many discussions on δ 13 C atm in ice cores and E.W. Wolff for discussions on the impact of climate variations during glacial/interglacial terminations on carbonate compensation effects in the subsequent warm periods.This is EPICA publication no.xxx and TALDICE publication no.XXX.