Search

The recovery of a 1.5 million yr long ice core from Antarctica represents a keystone of our understanding of Quaternary climate, the progression of glaciation over this time period and the role of greenhouse gas cycles in this progression. Here we tackle the question of where such ice may still be found in the Antarctic ice sheet. We can show that such old ice is most likely to exist in the plateau area of the East Antarctic ice sheet (EAIS) without stratigraphic disturbance and should be able to be recovered after careful presite selection studies. Based on a simple ice and heat flow model and glaciological observations, we conclude that positions in the vicinity of major domes and saddle position on the East Antarctic Plateau will most likely have such old ice in store and represent the best study areas for dedicated reconnaissance studies in the near future. In contrast to previous ice core drill site selections, however, we strongly suggest significantly reduced ice thickness to avoid bottom melting. For example for the geothermal heat flux and accumulation conditions at Dome C, an ice thickness lower than but close to about 2500 m would be required to find 1.5 Myr old ice (i.e., more than 700 m less than at the current EPICA Dome C drill site). Within this constraint, the resolution of an Oldest-Ice record and the distance of such old ice to the bedrock should be maximized to avoid ice flow disturbances, for example, by finding locations with minimum geothermal heat flux. As the geothermal heat flux is largely unknown for the EAIS, this parameter has to be carefully determined beforehand. In addition, detailed bedrock topography and ice flow history has to be reconstructed for candidates of an Oldest-Ice ice coring site. Finally, we argue strongly for rapid access drilling before any full, deep ice coring activity commences to bring datable samples to the surface and to allow an age check of the oldest ice.

Author(s): EPICA Community Members; EPICA Community Members; C. Barbante [1, 2]; J.-M. Barnola [3]; S. Becagli [4]; J. Beer [5]; M. Bigler [6, 7]; C. Boutron [3]; T. Blunier [6]; [...]

Author(s): O. Watanabe [1]; J. Jouzel (corresponding author) [2]; S. Johnsen [3]; F. Parrenin [2, 4]; H. Shoji [5]; N. Yoshida [6, 7] Recent ice core studies have raised the [...]

Documenting past changes in the East Antarctic surface mass balance is important to improve ice core chronologies and to constrain the ice-sheet contribution to global mean sea-level change. Here we reconstruct past changes in the ratio of surface mass balance (SMB ratio) between the EPICA Dome C (EDC) and Dome Fuji (DF) East Antarctica ice core sites, based on a precise volcanic synchronization of the two ice cores and on corrections for the vertical thinning of layers. During the past 216 000 a, this SMB ratio, denoted SMBEDC/SMBDF, varied between 0.7 and 1.1, being small during cold periods and large during warm periods. Our results therefore reveal larger amplitudes of changes in SMB at EDC compared with DF, consistent with previous results showing larger amplitudes of changes in water stable isotopes and estimated surface temperature at EDC compared with DF. Within the last glacial inception (Marine Isotope Stages, MIS-5c and MIS-5d), the SMB ratio deviates by up to 0.2 from what is expected based on differences in water stable isotope records. Moreover, the SMB ratio is constant throughout the late parts of the current and last interglacial periods, despite contrasting isotopic trends. [ABSTRACT FROM AUTHOR]

Two deep ice cores, Dome Fuji (DF) and EPICA Dome C (EDC), drilled at remote dome summits in Antarctica, were volcanically synchronized to improve our understanding of their chronologies. Within the past 216 kyr, 1401 volcanic tie points have been identified. DFO2006 is the chronology for the DF core that strictly follows O2 =N2 age constraints with interpolation using an ice flow model. AICC2012 is the chronology for five cores, including the EDC core, and is characterized by glaciological approaches combining ice flow modelling with various age markers. A precise comparison between the two chronologies was performed. The age differences between them are within 2 kyr, except at Marine Isotope Stage (MIS) 5. DFO2006 gives ages older than AICC2012, with peak values of 4.5 and 3.1 kyr at MIS 5d and MIS 5b, respectively. Accordingly, the ratios of duration (AICC2012 = DFO2006) range between 1.4 at MIS 5e and 0.7 at MIS 5a. When making a comparison with accurately dated speleothem records, the age of DFO2006 agrees well at MIS 5d, while the age of AICC2012 agrees well at MIS 5b, supporting their accuracy at these stages. In addition, we found that glaciological approaches tend to give chronologies with younger ages and with longer durations than age markers suggest at MIS 5d-6. Therefore, we hypothesize that the causes of the DFO2006-AICC2012 age differences at MIS 5 are (i) overestimation in surface mass balance at around MIS 5d-6 in the glaciological approach and (ii) an error in one of the O2 =N2 age constraints by ~3 kyr at MIS 5b. Overall, we improved our knowledge of the timing and duration of climatic stages at MIS 5. This new understanding will be incorporated into the production of the next common age scale. Additionally, we found that the deuterium signals of ice, δDice, at DF tends to lead the one at EDC, with the DF lead being more pronounced during cold periods. The lead of DF is by +710 years (maximum) at MIS 5d, -230 years (minimum) at MIS 7a and +60 to +126 years on average. [ABSTRACT FROM AUTHOR]

Polar ice cores provide exceptional archives of past environmental conditions. The dating of ice cores and the estimation of the age-scale uncertainty are essential to interpret the climate and environmental records that they contain. It is, however, a complex problem which involves different methods. Here, we present IceChrono1, a new probabilistic model integrating various sources of chronological information to produce a common and optimized chronology for several ice cores, as well as its uncertainty. IceChrono1 is based on the inversion of three quantities: the surface accumulation rate, the lock-in depth (LID) of air bubbles and the thinning function. The chronological information integrated into the model are models of the sedimentation process (accumulation of snow, densification of snow into ice and air trapping, ice flow), ice- and air-dated horizons, ice and air depth intervals with known durations, 1depth observations (depth shift between synchronous events recorded in the ice and in the air) and finally air and ice stratigraphic links in between ice cores. The optimization is formulated as a least squares problem, implying that all densities of probabilities are assumed to be Gaussian. It is numerically solved using the Levenberg-Marquardt algorithm and a numerical evaluation of the model's Jacobian. IceChrono follows an approach similar to that of the Datice model which was recently used to produce the AICC2012 (Antarctic ice core chronology) for four Antarctic ice cores and one Greenland ice core. IceChrono1 provides improvements and simplifications with respect to Datice from the mathematical, numerical and programming point of views. The capabilities of IceChrono1 are demonstrated on a case study similar to the AICC2012 dating experiment. We find results similar to those of Datice, within a few centuries, which is a confirmation of both IceChrono1 and Datice codes. We also test new functionalities with respect to the original version of Datice: observations as ice intervals with known durations, correlated observations, observations as air intervals with known durations and observations as mixed ice-air stratigraphic links. IceChrono1 is freely available under the General Public License v3 open source license. [ABSTRACT FROM AUTHOR]

A recent coherent chronology has been built for four Antarctic ice cores and the NorthGRIP (NGRIP) Greenland ice core (Antarctic Ice Core Chronology 2012, AICC2012) using a Bayesian approach for ice core dating (Datice). When building the AICC2012 chronology, and in order to prevent any confusion with official ice core chronology, the AICC2012 chronology for NGRIP was forced to fit exactly the GICC05 chronology based on layer counting. However, such a strong tuning did not satisfy the hypothesis of independence of background parameters and observations for the NGRIP core, as required by Datice. We present here the implementation in Datice of a new type of markers that is better suited for constraints deduced from layer counting: the duration constraints. Estimating the global error on chronology due to such markers is not straightforward and implies some assumption on the correlation between individual counting errors for each interval of duration. We validate this new methodological implementation by conducting twin experiments and a posteriori diagnostics on the NGRIP ice core. Several sensitivity tests on marker sampling and correlation between counting errors were performed to provide some guidelines when using such a method for future dating experiments. Finally, using these markers for NGRIP in a five-core dating exercise with Datice leads to new chronologies that do not differ by more than 410 years from AICC2012 for Antarctic ice cores and 150 years from GICC05 for NGRIP over the last 60 000 years. [ABSTRACT FROM AUTHOR]

Documenting past changes in the East Antarctic surface mass balance is important to improve ice core chronologies and to constrain the ice sheet contribution to global mean sea level. Here we reconstruct the past changes in the ratio of surface mass balance (SMB ratio) between the EPICA Dome C (EDC) and Dome Fuji (DF) East Antarctica ice core sites, based on a precise volcanic synchronisation of the two ice cores and on corrections for the vertical thinning of layers. During the past 216 000 years, this SMB ratio, denoted SMBEDC=SMBDF, varied between 0.7 and 1.1, decreasing during cold periods and increasing during warm periods. While past climatic changes have been depicted as homogeneous along the East Antarctic Plateau, our results reveal larger amplitudes of changes in SMB at EDC compared to DF, consistent with previous results showing larger amplitudes of changes in water stable isotopes and estimated surface temperature at EDC compared to DF. Within interglacial periods and during the last glacial inception (Marine Isotope Stages, MIS-5c and MIS-5d), the SMB ratio deviates by up to 30% from what is expected based on differences in water stable isotope records. Moreover, the SMB ratio is constant throughout the late parts of the current and last interglacial periods, despite contrasting isotopic trends. These SMB ratio changes not closely related to isotopic changes are one of the possible causes of the observed gaps between the ice core chronologies at DF and EDC. Such changes in SMB ratio may have been caused by (i) climatic processes related to changes in air mass trajectories and local climate, (ii) glaciological processes associated with relative elevation changes, or (iii) a combination of climatic and glaciological processes, such as the interaction between changes in accumulation and in the position of the domes. Our inferred SMB ratio history has important implications for ice sheet modeling (for which SMB is a boundary condition) or atmospheric modeling (our inferred SMB ratio could serve as a test). [ABSTRACT FROM AUTHOR]

The end of the Last Glacial Maximum (Termination I), roughly 20 thousand years ago (ka), was marked by cooling in the Northern Hemisphere, a weakening of the Asian monsoon, a rise in atmospheric CO2 concentrations and warming over Antarctica. The sequence of events associated with the previous glacial-interglacial transition (Termination II), roughly 136 ka, is less well constrained. Here we present high-resolution records of atmospheric CO2 concentrations and isotopic composition of N2-an atmospheric temperature proxy-from air bubbles in the EPICA Dome C ice core that span Termination II. We find that atmospheric CO2 concentrations and Antarctic temperature started increasing in phase around 136 ka, but in a second phase of Termination II, from 130.5 to 129 ka, the rise in atmospheric CO2 concentrations lagged that of Antarctic temperature unequivocally. We suggest that during this second phase, the intensification of the low-latitude hydrological cycle resulted in the development of a CO2 sink, which counteracted the CO2 outgassing from the Southern Hemisphere oceans over this period. [ABSTRACT FROM AUTHOR]

The deep polar ice cores provide reference records commonly employed in global correlation of past climate events. However, temporal divergences reaching up to several thousand years (ka) exist between ice cores over the last climatic cycle. In this context, we are hereby introducing the Antarctic Ice Core Chronology 2012 (AICC2012), a new and coherent timescale developed for four Antarctic ice cores, namely Vostok, EPICA Dome C (EDC), EPICA Dronning Maud Land (EDML) and Talos Dome (TALDICE), alongside the Greenlandic NGRIP record. The AICC2012 timescale has been constructed using the Bayesian tool Datice (Lemieux-Dudon et al., 2010) that combines glaciological inputs and data constraints, including a wide range of relative and absolute gas and ice stratigraphic markers. We focus here on the last 120 ka, whereas the companion paper by Bazin et al. (2013) focuses on the interval 120-800 ka. Compared to previous timescales, AICC2012 presents an improved timing for the last glacial inception, respecting the glaciological constraints of all analyzed records. Moreover, with the addition of numerous new stratigraphic markers and improved calculation of the lock-in depth (LID) based on d15N data employed as the Datice background scenario, the AICC2012 presents a slightly improved timing for the bipolar sequence of events over Marine Isotope Stage 3 associated with the seesaw mechanism, with maximum differences of about 600 yr with respect to the previous Datice-derived chronology of Lemieux-Dudon et al. (2010), hereafter denoted LD2010. Our improved scenario confirms the regional differences for the millennial scale variability over the last glacial period: while the EDC isotopic record (events of triangular shape) displays peaks roughly at the same time as the NGRIP abrupt isotopic increases, the EDML isotopic record (events characterized by broader peaks or even extended periods of high isotope values) reached the isotopic maximum several centuries before. It is expected that the future contribution of both other long ice core records and other types of chronological constraints to the Datice tool will lead to further refinements in the ice core chronologies beyond the AICC2012 chronology. For the time being however, we recommend that AICC2012 be used as the preferred chronology for the Vostok, EDC, EDML and TALDICE ice core records, both over the last glacial cycle (this study), and beyond (following Bazin et al., 2013). The ages for NGRIP in AICC2012 are virtually identical to those of GICC05 for the last 60.2 ka, whereas the ages beyond are independent of those in GICC05modelext (as in the construction of AICC2012, the GICC05modelext was included only via the background scenarios and not as age markers). As such, where issues of phasing between Antarctic records included in AICC2012 and NGRIP are involved, the NGRIP ages in AICC2012 should therefore be taken to avoid introducing false offsets. However for issues involving only Greenland ice cores, there is not yet a strong basis to recommend superseding GICC05modelext as the recommended age scale for Greenland ice cores. [ABSTRACT FROM AUTHOR]

An accurate and coherent chronological framework is essential for the interpretation of climatic and environmental records obtained from deep polar ice cores. Until now, one common ice core age scale had been developed based on an inverse dating method (Datice), combining glaciological modelling with absolute and stratigraphic markers between 4 ice cores covering the last 50 ka (thousands of years before present) (Lemieux-Dudon et al., 2010). In this paper, together with the companion paper of Veres et al. (2013), we present an extension of this work back to 800 ka for the NGRIP, TALDICE, EDML, Vostok and EDC ice cores using an improved version of the Datice tool. The AICC2012 (Antarctic Ice Core Chronology 2012) chronology includes numerous new gas and ice stratigraphic links as well as improved evaluation of background and associated variance scenarios. This paper concentrates on the long timescales between 120-800 ka. In this framework, new measurements of d18Oatm over Marine Isotope Stage (MIS) 11-12 on EDC and a complete d18Oatm record of the TALDICE ice cores permit us to derive additional orbital gas age constraints. The coherency of the different orbitally deduced ages (from d18Oatm, dO2/N2 and air content) has been verified before implementation in AICC2012. The new chronology is now independent of other archives and shows only small differences, most of the time within the original uncertainty range calculated by Datice, when compared with the previous ice core reference age scale EDC3, the Dome F chronology, or using a comparison between speleothems and methane. For instance, the largest deviation between AICC2012 and EDC3 (5.4 ka) is obtained around MIS 12. Despite significant modifications of the chronological constraints around MIS 5, now independent of speleothem records in AICC2012, the date of Termination II is very close to the EDC3 one. [ABSTRACT FROM AUTHOR]

Correct estimation of the firn lock-in depth is essential for correctly linking gas and ice chronologies in ice core studies. Here, two approaches to constrain the firn depth evolution in Antarctica are presented over the last deglaciation: outputs of a firn densification model, and measurements of δ15N of N2 in air trapped in ice core, assuming that δ15N is only affected by gravitational fractionation in the firn column. Since the firn densification process is largely governed by surface temperature and accumulation rate, we have investigated four ice cores drilled in coastal (Berkner Island, BI, and James Ross Island, JRI) and semi-coastal (TALDICE and EPICA Dronning Maud Land, EDML) Antarctic regions. Combined with available ice core air-δ15N measurements from the EPICA Dome C (EDC) site, the studied regions encompass a large range of surface accumulation rates and temperature conditions.

Our δ15N profiles reveal a heterogeneous response of the firn structure to glacial–interglacial climatic changes. While firn densification simulations correctly predict TALDICE δ15N variations, they systematically fail to capture the large millennial-scale δ15N variations measured at BI and the δ15N glacial levels measured at JRI and EDML – a mismatch previously reported for central East Antarctic ice cores.

New constraints of the EDML gas–ice depth offset during the Laschamp event (~41 ka) and the last deglaciation do not favour the hypothesis of a large convective zone within the firn as the explanation of the glacial firn model–δ15N data mismatch for this site. While we could not conduct an in-depth study of the influence of impurities in snow for firnification from the existing datasets, our detailed comparison between the δ15N profiles and firn model simulations under different temperature and accumulation rate scenarios suggests that the role of accumulation rate may have been underestimated in the current description of firnification models. [ABSTRACT FROM AUTHOR]

The article examines the effect of atmospheric carbon dioxide (CO2) on global climatic changes, with a focus on CO2 variations relative to time and temperature. Topics include preserved CO2 records in Antarctic ice cores which indicate past temperature levels; a proposed relative age scale based on the concentration of CO2 and temperatures in Antarctica during deglacial warming periods; and trends in CO2 concentration which suggests a direct correlation between increasing temperatures and CO2 concentration.

The Toba eruption that occurred some 74 ka ago in Sumatra, Indonesia, is among the largest volcanic events on Earth over the last 2 million years. Tephra from this eruption has been spread over vast areas in Asia, where it constitutes a major time marker close to the Marine Isotope Stage 4/5 boundary. As yet, no tephra associated with Toba has been identified in Greenland or Antarctic ice cores. Based on new accurate dating of Toba tephra and on accurately dated European stalagmites, the Toba event is known to occur between the onsets of Greenland interstadials (GI) 19 and 20. Furthermore, the existing linking of Greenland and Antarctic ice cores by gas records and by the bipolar seesaw hypothesis suggests that the Antarctic counterpart is situated between Antarctic Isotope Maxima (AIM) 19 and 20 In this work we suggest a direct synchronization of Greenland (NGRIP) and Antarctic (EDML) ice cores at the Toba eruption based on matching of a pattern of bipolar volcanic spikes. Annual layer counting between volcanic spikes in both cores allows for a unique match. We first demonstrate this bipolar matching technique at the already synchronized Laschamp geomagnetic excursion (41 ka BP) before we apply it to the suggested Toba interval. The Toba synchronization pattern covers some 2000 yr in GI-20 and AIM- 19/20 and includes nine acidity peaks that are recognized in both ice cores. The suggested bipolar Toba synchronization has decadal precision. It thus allows a determination of the exact phasing of inter-hemispheric climate in a time interval of poorly constrained ice core records, and it allows for a discussion of the climatic impact of the Toba eruption in a global perspective. The bipolar linking gives no support for a long-term global cooling caused by the Toba eruption as Antarctica experiences a major warming shortly after the event. Furthermore, our bipolar match provides a way to place palaeoenvironmental records other than ice cores into a precise climatic context. [ABSTRACT FROM AUTHOR]

Correct estimate of the firn lock-in depth is essential for correctly linking gas and ice chronologies in ice cores studies. Here, two approaches to constrain the firn depth evolution in Antarctica are presented over the last deglaciation: output of a firn densification model and measurements of δ15 N of N2 in air trapped in ice core. Since the firn densification process is largely governed by surface temperature and accumulation rate, we have investigated four ice cores drilled in coastal (Berkner Island, BI, and James Ross Island, JRI) and semi coastal (TALDICE and EPICA Dronning Maud Land, EDML) Antarctic regions. Combined with available δ15 N measurements performed from the EPICA Dome C (EDC) site, the studied regions encompass a large range of surface accumulation rate and temperature conditions. While firn densification simulations are able to correctly represent most of the δ15N trends over the last deglaciation measured in the EDC, BI, TALDICE and EDML ice cores, they systematically fail to capture BI and EDML δ15N glacial levels, a mismatch previously seen for Central East Antarctic ice cores. Using empirical constraints of the EDML gas-ice depth offset during the Laschamp event (~41 ka), we can rule out the existence of a large convective zone as the explanation of the glacial firn model- δ15N data mismatch for this site. The good match between modelled and measured δ15N at TALDICE as well as the lack of any clear correlation between insoluble dust concentration in snow and δ15N records in the different ice cores suggest that past changes in loads of impurities are not the only main driver of glacial-interglacial changes in firn lock-in depth. We conclude that firn densification dynamics may instead be driven mostly by accumulation rate changes. The mismatch between modelled and measured δ15N may be due to inaccurate reconstruction of past accumulation rate or underestimated influence of accumulation rate in firnification models. [ABSTRACT FROM AUTHOR]

The deep polar ice cores provide reference records commonly employed in global correlation of past climate events. However, temporal divergences reaching up to several thousand years (ka) exist between ice cores over the last climatic cycle. In this context, we are hereby introducing the Antarctic Ice Core Chronology 2012 (AICC2012), a new and coherent timescale developed for four Antarctic ice cores, namely Vostok, EPICA Dome C (EDC), EPICA Dronning Maud Land (EDML) and Talos Dome (TALDICE), alongside the Greenlandic NGRIP record. The AICC2012 time scale has been constructed using the Bayesian tool Datice (Lemieux-Dudon et al., 2010) that combines glaciological inputs and data constraints, including a wide range of relative and absolute gas and ice stratigraphic markers. We focus here on the last 120 ka, whereas the companion paper by Bazin et al., (2012) focuses on the interval 120-800 ka. Compared to previous timescales, AICC2012 presents an improved timing for the last glacial inception respecting the glaciological constraints of all analyzed records. Moreover, with the addition of numerous new stratigraphic markers and improved calculation of the lock-in depth (LID) based on δ15 N data employed as the Datice background scenario, the AICC2012 presents a new timing for the bipolar sequence of events over Marine Isotope Stage 3 associated with the see-saw mechanism, with maximum differences of about 500 yr with respect to the previous Datice-derived chronology of Lemieux-Dudon et al. (2010), hereafter denoted LD2010. Our improved scenario confirms the regional differences for the millennial scale variability over the last glacial period: while the EDC isotopic record (events of triangular shape) displays peaks roughly at the same time as the NGRIP abrupt isotopic increases, the EDML isotopic record (events characterized by broader peaks or even extended periods of high isotope values) reached the isotopic maximum several centuries before. [ABSTRACT FROM AUTHOR]

An accurate and coherent chronological framework is essential for the interpretation of climatic and environmental records obtained from deep polar ice cores. Until now, one common ice core age scale has been developed based on an inverse dating method (Datice) combining glaciological modelling with absolute and stratigraphic markers between 4 ice cores covering the last 50 ka (thousand of years before present) (Lemieux- Dudon et al., 2010). In this paper, together with the companion paper of Veres et al. (2012), we present an extension of this work back to 800 ka for the NGRIP, TALDICE, EDML, Vostok and EDC ice cores using an improved version of the Datice tool. The AICC2012 (Antarctic Ice Core Chronology 2012) chronology includes numerous new gas and ice stratigraphic links as well as improved evaluation of background and associated variance scenarios. This paper concentrates on the long timescales between 120-800 ka. In this frame, new measurements of δ18 Oatm over Marine Isotope Stage (MIS) 11-12 on EDC and a complete δ18 Oatm record of the TALDICE ice cores permit us to derive new orbital gas age constraints. The coherency of the different orbitally deduced ages (from δ18 Oatm, δO2/N2 and air content) has been verified before implementation in AICC2012. The new chronology shows only small differences, well within the original uncertainty range, when compared with the previous ice core reference age scale EDC3. For instance, the duration of the last four interglacial periods is not affected by more than 5%. The largest deviation between AICC2012 and EDC3 (4.4 ka) is obtained around MIS 12. Despite significant modifications of the chronological constraints around MIS 5, now independent of speleothem records in AICC2012, the date of Termination II is very close to the EDC3 one. [ABSTRACT FROM AUTHOR]

The Toba eruption that occurred some 74 kyr ago in Sumatra, Indonesia, is among the largest volcanic events on Earth over the last 2 million years. Tephra from this eruption has been spread over vast areas in Asia where it constitutes a major time marker close to the Marine Isotope Stage 4/5 boundary. As yet, no tephra associated with Toba has been identified in Greenland or Antarctic ice cores. Based on new accurate dating of Toba tephra from Malaysia and on accurately dated European stalagmites the Toba event is known to occur between the onsets of Greenland Interstadials (GI) 19 and 20. Furthermore, the existing linking of Greenland and Antarctic ice cores by gas records and by the bipolar seesaw hypothesis suggests that the Antarctic counterpart is situated between Antarctic Isotope Maxima (AIM) 19 and 20. In this work we suggest a direct synchronization of Greenland (NGRIP) and Antarctic (EDML) ice cores at the Toba eruption based on matching of a pattern of bipolar volcanic spikes. Annual layer counting between volcanic spikes in both cores allows for a unique match. We first demonstrate this bipolar matching technique at the already synchronized Laschamp geomagnetic excursion (41 kyr BP) before we apply it to the suggested Toba interval. The Toba synchronization pattern covers some 2000 yr in GI-20 and AIM 19/20 and includes nine acidity peaks that are recognized in both ice cores. The suggested bipolar Toba synchronization has decadal precision. It thus allows a determination of the exact phasing of inter-hemispheric climate in a time interval of poorly constrained ice core records, and it allows for a discussion of the climatic impact of the Toba eruption in a global perspective. Furthermore, our bipolar match provides a way to place paleo-environmental records other than ice cores into a precise climatic context. [ABSTRACT FROM AUTHOR]

The main variations of ice volume of the last million years can be explained from orbital parameters by assuming climate oscillates between two states: glaciations and deglaciations (Parrenin and Paillard, 2003; Imbrie et al., 2011) (or terminations). An additional combination of ice volume and orbital parameters seems to form the trigger of a deglaciation, while only orbital parameters seem to play a role in the triggering of glaciations. Here we present an optimized conceptual model which realistically reproduce ice volume variations during the past million years and in particular the timing of the 11 canonical terminations. We show that our model looses sensitivity to initial conditions only after ∼ 200 kyr at maximum: the ice volume observations form a strong attractor. Both obliquity and precession seem necessary to reproduce all 11 terminations and both seem to play approximately the same role. More precisely, obliquity plays a fundamental role in the triggering of termination VI (∼530 kyr BP), while precession plays a fundamental role in the triggering of termination VIII (∼720 kyr ago). [ABSTRACT FROM AUTHOR]

The main variations of ice volume of the last million years can be explained from orbital parameters (Laskar et al., 2004) by assuming climate oscillates between two states: glaciations and deglaciations (Parrenin and Paillard, 2003; Imbrie et al., 2011) (or terminations). An additional combination of ice volume and orbital parameters seems to form the trigger of a deglaciation (Parrenin and Paillard, 2003), while only orbital parameters seem to play a role in the triggering of glaciations (Khodri et al., 2001). Here we present conceptual models, which realistically reproduce ice volume variations during the past million years and in particular, the timing of the 11 canonical terminations. Exploring our model's parameters exhaustively, we show that obliquity plays a fundamental role in the triggering of termination VI (~ 530 kyr BP), while precession plays a fundamental role in the triggering of termination VIII (~ 720 kyr ago). [ABSTRACT FROM AUTHOR]

This study aims at refining the synchronisation between the EPICA Dome C (EDC) and Vostok ice cores in the time interval 0-145 kyr BP by using the volcanic signatures. 102 common volcanic events were identified by using continuous electrical conductivity (ECM), di-electrical profiling (DEP) and sulfate measurements while trying to minimize the distortion of the glaciological chronologies. This is an update and a continuation of previous works performed over the 0-45 kyr interval that provided 56 tie points to the ice core chronologies (Udisti et al., 2004). This synchronisation will serve to establish Antarctic Ice Core Chronology 2012, the next synchronised Antarctic dating. A change of slope in the EDC-depth/Vostok-depth diagram is probably related to a change of accumulation regime as well as to a change of ice thickness upstream of the Lake Vostok, but we did not invoke any significant temporal change of surface accumulation at EDC relative to Vostok. No significant phase difference is detected between the EDC and Vostok isotopic records, but depth shifts between the Vostok 3G and 5G ice cores prevent from looking at this problem accurately. Three possible candidates for the Toba volcanic super-eruption ∼73 kyr ago are suggested in the Vostok and EDC volcanic records. Neither the ECM, DEP nor the sulfate fingerprints for these 3 events are significantly larger than many others in the records. [ABSTRACT FROM AUTHOR]

We compare a variety of methods for estimating the gas/ice depth offset (Δdepth) at EPICA Dome C (EDC, East Antarctica). (1) Purely based on modelling efforts, Δdepth can be estimated combining a firn densification with an ice flow model. Observations allow direct and indirect estimate of Δdepth. (2) The diffusive column height can be estimated from δ15N and converted to Δdepth using an ice flow model and assumptions about past average firn density and thickness of the convective zone. (3) Ice and gas synchronisation of the EDC ice core to the GRIP, EDML and TALDICE ice cores shifts the ice/gas offset problem into higher accumulation ice cores where it can be more accurately evaluated. (4) Finally, the bipolar seesaw hypothesis allows us to synchronise the ice isotopic record with the gas CH4 record, the later being taken as a proxy of Greenland temperature. The bipolar seesaw antiphase relationship is generally supported by the ice-gas cross synchronisation between EDC and the GRIP, EDML and TALDICE ice cores, which provide support for method 4. Applying the bipolar seesaw hypothesis to the deeper section of the EDC core confirms that the ice flow is complex and can help improving our reconstruction of the thinning function and thus of the EDC age scale. We confirm that method 1 overestimates the glacial Δdepth at EDC and we suggested that it is due to an overestimation of the glacial Close Off Depth by the firn densification model. In contrast we find that the glaciological models probably underestimate the Δdepth during termination II. Finally, we show that method 2 based on 15N data produces for the last deglaciation a Δdepth estimate which is in good agreement with methods 3 and 4. [ABSTRACT FROM AUTHOR]

This study aims at refining the synchronisation between the EPICA Dome C (EDC) and Vostok ice cores in the time interval 0-145 kyr BP by using the volcanic signatures. 111 common volcanic events were identified by using continuous electrical conductivity (ECM), di-electrical profiling (DEP) and sulfate measurements while trying to minimize the distortion of the glaciological chronologies. This is an update and a continuation of previous works performed over the 0-45 kyr interval which provided 56 tie points to the ice core chronologies (Udisti et al., 2004). This synchronisation will serve for the establishment of the next synchronised Antarctic dating. A change of slope in the EDC-depth/Vostok-depth diagram is probably related to a change of accumulation regime as well as to a change of ice thickness upstream of the Vostok lake, but we did not invoke any significant temporal change of surface accumulation at EDC relative to Vostok. A significant phase difference is detected between the EDC and Vostok isotopic records during the 95-120 kyr interval, but not during Termination II. Three possible candidates for the Toba volcanic super-eruption ∼73 kyr ago are suggested in the Vostok and EDC volcanic records. However the ECM, DEP and sulfate fingerprints for these three events are not significantly larger than many others in the records. [ABSTRACT FROM AUTHOR]

A new deep ice core drilling program, TALDICE, has been successfully handled by a European team at Talos Dome, in the Ross Sea sector of East Antarctica, down to 1620m depth. Using stratigraphic markers and a new inverse method, we produce the first official chronology of the ice core, called TALDICE-1. We show that it notably improves an a priori chronology resulting from a one-dimensional ice flow model. It is in agreement with a posteriori controls of the resulting accumulation rate and thinning function along the core. An absolute uncertainty of only 300 yr is obtained over the course of the last deglaciation. This uncertainty remains lower than 600 yr over Marine Isotope Stage 3, back to 50 kyr BP. The phasing of the TALDICE ice core climate record with respect to the central East Antarctic plateau and Greenland records can thus be determined with a precision allowing for a discussion of the mechanisms at work at submillennial time scales. [ABSTRACT FROM AUTHOR]

The Greenland Ice Core Chronology 2005 (GICC05) is a time scale based on annual layer counting of high-resolution records from Greenland ice cores. Whereas the Holocene part of the time scale is based on various records from the DYE-3, the GRIP, and the NorthGRIP ice cores, the glacial part is solely based on NorthGRIP records. Here we present an 18 ka extension of the time scale such that GICC05 continuously covers the past 60 ka. The new section of the time scale places the onset of Greenland Interstadial 12 (GI-12) at 46.9±1.0 ka b2k (before year AD 2000), the North Atlantic Ash Zone II layer in GI-15 at 55.4±1.2 ka b2k, and the onset of GI-17 at 59.4±1.3 ka b2k. The error estimates are derived from the accumulated number of uncertain annual layers. In the 40-60 ka interval, the new time scale has a discrepancy with the Meese-Sowers GISP2 time scale of up to 2.4 ka. Assuming that the Greenland climatic events are synchronous with those seen in the Chinese Hulu Cave speleothem record, GICC05 compares well to the time scale of that record with absolute age differences of less than 800 years throughout the 60 ka period. The new time scale is generally in close agreement with other independently dated records and reference horizons, such as the Laschamp geomagnetic excursion, the French Villars Cave and the Austrian Kleegruben Cave speleothem records, suggesting high accuracy of both event durations and absolute age estimates. [ABSTRACT FROM AUTHOR]

The article presents a response by F. Parrenin on the paper "Ice Thinning, Upstream Advection, and Non-Climactic Biases for 89% of the EDML Ice Core From a Nested Model of the Antarctic Ice Sheet" by P. Huybrechts et al. that was published in the previous issue in the journal. The editor notes that while the paper is well written, details on the models and sensitivity experiments are missing. He also questions the clarity of the coupling schemes. He suggests showing some snapshots of the LSM surface for several typical climates.

The article presents a response by F. Parrenin et al. to a comment on their paper "The EDC3 Chronology for the EPICA Dome C Ice Core" that was published in a previous issue of the journal. The authors of the paper offer the new EDC3 time scale for the EPICA Dome-C deep ice core. They provide a useful discussion on the different types of accuracies involved, such as absolute, relative and the duration of climatic events.

The article comments on the paper "Ice flow modelling at EPICA Dome C and Dome Fuji, East Antarctica" by F. Parrenin et al. The authors, in the submitted manuscript, summarized theoretical issues and provided further development of their previous work on the application of a 1-D glaciological model, first implemented to construct a chronology for the Vostok ice core.

The Greenland Ice Core Chronology 2005 (GICC05) is a time scale based on annual layer counting of high-resolution records from Greenland ice cores. Whereas the Holocene part of the time scale is based on various records from the DYE-3, the GRIP, and the NorthGRIP ice cores, the glacial part is solely based on NorthGRIP records. Here we present an 18 kyr extension of the time scale such that GICC05 continuously covers the past 60 kyr. The new section of the time scale places the onset of Greenland Interstadial 12 (GI-12) at 46.9±1.0 kyr b2k (before year AD 2000), the North Atlantic Ash Zone 2 layer in GI-15 at 55.4±1.2 kyr b2k, and the onset of GI-17 at 59.4±1.3 kyr b2k. The error estimates are derived from the accumulated number of uncertain annual layers and can be regarded as 1σ uncertainties. In the 40-60 kyr interval the new time scale has a discrepancy with the Meese-Sowers GISP2 time scale of up to 2.4 kyr, whereas GICC05 compares well to the dating of the Hulu Cave record with absolute age differences of less than 800 years throughout the 60 kyr period. The new time scale is generally in close agreement with other independently dated records and reference horizons, such as the Laschamp geomagnetic excursion and the Kleegruben speleothem record from the Austrian Alps, suggesting high accuracy of both event durations and absolute age estimates. [ABSTRACT FROM AUTHOR]

Gas is trapped in polar ice sheets at ~50-120 m below the surface and is therefore younger than the surrounding ice. Firn densification models are used to evaluate this ice age-gas age difference (Δage) in the past. However, such models need to be validated by data, in particular for periods colder than present day on the East Antarctic plateau. Here we bring new constraints to test a firn densification model applied to the EPICA Dome C (EDC) site for the last 50 kyr, by linking the EDC ice core to the EPICA Dronning Maud Land (EDML) ice core, both in the ice phase (using volcanic horizons) and in the gas phase (using rapid methane variations). We also use the structured 10Be peak, occurring 41 kyr before present (BP) and due to the low geomagnetic field associated with the Laschamp event, to experimentally estimate the Δage during this event. Our results seem to reveal an overestimate of the Δage by the firn densification model during the last glacial period at EDC. Tests with different accumulation rates and temperature scenarios do not entirely resolve this discrepancy. Although the exact reasons for the Δage overestimate at the two EPICA sites remain unknown at this stage, we conclude that current densification model simulations have deficits under glacial climatic conditions. Whatever the cause of the Δage overestimate, our finding suggests that the phase relationship between CO2 and EDC temperature previously inferred for the start of the last deglaciation (lag of CO2 by 800±600 yr) seems to be overestimated. [ABSTRACT FROM AUTHOR]

The EPICA (European Project for Ice Coring in Antarctica) Dome C drilling in East Antarctica has now been completed to a depth of 3260 m, at only a few meters above bedrock. Here we present the new EDC3 chronology, which is based on the use of 1) a snow accumulation and mechanical flow model, and 2) a set of independent age markers along the core. These are obtained by pattern matching of recorded parameters to either absolutely dated paleoclimatic records, or to insolation variations. We show that this new time scale is in excellent agreement with the Dome Fuji and Vostok ice core time scales back to 100 kyr within 1 kyr. Discrepancies larger than 3 kyr arise during MIS 5.4, 5.5 and 6, which points to anomalies in either snow accumulation or mechanical flow during these time periods. We estimate that EDC3 gives accurate event durations within 20% (2δ) back to MIS 11 and accurate absolute ages with a maximum uncertainty of 6 kyr back to 800 kyr. [ABSTRACT FROM AUTHOR]

A chronology called EDML1 has been developed for the EPICA ice core from Dronning Maud Land (EDML). EDMLI is closely interlinked with EDC3, the new chronology for the EPICA ice core from Dome-C (EDC) through a stratigraphic match between EDML and EDC that consists of 322 volcanic match points over the last 128 ka. The EDC3 chronology comprises a glaciological model at EDC, which is constrained and later selectively tuned using primary dating information from EDC as well as from EDML, the latter being transferred using the tight stratigraphic link between the two cores. Finally, EDMLI was built by exporting EDC3 to EDML. For ages younger than 41 ka BP the new synchronized time scale EDML1/EDC3 is based on dated volcanic events and on a match to the Greenlandic ice core chronology GICC05 via 10Be and methane. The internal consistency between EDML1 and EDC3 is estimated to be typically ~6 years and always less than 450 years over the last 128 ka (always less than 130 years over the last 60 ka), which reflects an unprecedented synchrony of time scales. EDML1 ends at 150 ka BP (2417 m depth) because the match between EDML and EDC becomes ambiguous further down. This hints at a complex ice flow history for the deepest 350 m of the EDML ice core. [ABSTRACT FROM AUTHOR]

A common time scale for the EPICA ice cores from Dome C (EDC) and Dronning Maud Land (EDML) has been established. Since the EDML core was not drilled nn a dome, the development of the EDMLI time scale for the EPICA ice core drilled in Dronning Maud Land was based on the creation of a detailed stratigraphic link between EDML and EDC. which was dated by a simpler ID ice-flow model. The synchronisation between the two EPICA ice cores was done through the identification of several common volcanic signatures. This paper describes the rigorous method, using the signature of volcanic sulfate, which was employed for the last 52 kyr of the record. We estimated the discrepancies between the modelled EDC and EDML glaciological age scales during the studied period, by evaluating the ratio R of the apparent duration of temporal intervals between pairs of isochrones. On average R ranges between 0.8 and 1.2 corresponding to an uncertainty of up to 20% in the estimate of the time duration in at least one of the two ice cores. Significant deviations of R up to 1.4-1.5 are observed between 18 and 28 kyr before present (BP), where present is defined as 1950. At this stage our approach does not allow us unequivocally to find out which of the models is affected by errors, but assuming that the thinning function at both sites and accumulation history at Dome C (which was drilled on a dome) are correct, this anomaly can be ascribed to a complex spatial accumulation variability (which may be different in the past compared to the present day) upstream of the EDML core. [ABSTRACT FROM AUTHOR]

While there are no indications of mixing back to 800 000 years in the EPICA Dome C ice core record, comparison with marine sediment records shows significant differences in the timing and duration of events prior to stage 11 (∼430ka, thousands of years before 1950). A relationship between the isotopic composition of atmospheric oxygen (δ18O of O2, noted δ18Oatm) and daily northern hemisphere summer insolation has been observed for the youngest four climate cycles. Here we use this relationship with new δ18O of O2 measurements to show that anomalous flow in the bottom 500 m of the core distorts the duration of events by up to a factor of 2. By tuning δ18Oatm to orbital precession we derive a corrected thinning function and present a revised age scale for the interval corresponding to Marine Isotope Stages 11-20 in the EPICA Dome C ice core. Uncertainty in the phasing of δ18Oatm with respect to insolation variations in the precession band limits the accuracy of this new agescale to ±6 kyr (thousand of years). The previously reported ∼ 30 kyr duration of interglacial stage 11 is unchanged. In contrast, the duration of stage 15.1 is reduced by a factor of 2, from 31 to 16 kyr. [ABSTRACT FROM AUTHOR]

One-dimensional (1-D) ice flow models are used to construct the age scales at the Dome C and Dome Fuji drilling sites (East Antarctica). The poorly constrained glaciological parameters at each site are recovered by fitting independent age markers identified within each core. We reconstruct past accumulation rates, that are larger than those modelled using the classical vapour saturation pressure relationship during glacial periods by up to a factor 1.5. During the Early Holocene, changes in reconstructed accumulation are not linearly related to changes in ice isotopic composition. A simple model of past elevation changes is developed and shows an amplitude variation of 110-120m at both sites. We suggest that there is basal melting at Dome C (0.56±0.19 mm/yr). The reconstructed velocity profile is highly non-linear at both sites, which suggests complex ice flow effects. This induces a non-linear thinning function in both drilling sites, which is also characterized by bumps corresponding to variations in ice thickness with time. [ABSTRACT FROM AUTHOR]

While there are no indications of mixing back to 800 000 years in the EPICA Dome C ice core record, comparison with marine sediment records shows significant differences in the timing and duration of events prior to stage 11 (∼430 ka, thousand of years before 1950). A relationship between the isotopic composition of atmospheric oxygen (δ18O of O2, noted δ18Oatm) and daily northern hemisphere summer insolation has been observed for the youngest four climate cycles. Here we use this relationship with new δ18O of O2 measurements to show that anomalous flow in the bottom 500m of the core distorts the duration of events by up to a factor of 2. By tuning ∼18Oatm to orbital precession we derive a corrected thinning function and present a revised age scale for the interval corresponding to Marine Isotope Stages 11-20 in the EPICA Dome C ice core. Uncertainty in the phasing of δ18Oatm with respect to insolation variations in the precession band limits the accuracy of this new agescale to ±6 kyr (thousand of years). The previously reported ∼30 kyr duration of interglacial stage 11 is unchanged. In contrast, the duration of stage 15.1 is reduced by a factor of 2, from 31 to 16 kyr. [ABSTRACT FROM AUTHOR]

1-D ice flow models are used to construct the age scales at the Dome C and Dome Fuji drilling sites (East Antarctica). The poorly constrained glaciological parameters at each sites are recovered by fitting independent age markers identified on each core. We construct past accumulation rates, that are larger than those modeled using the classical vapour saturation pressure relationship during glacial periods by up to a factor 1.5. During the Early Holocence, changes in reconstructed accumulation are not linearly related to changes in ice isotopic composition. A simple model of past elevation changes is developed and shows and amplitude variation of 110-120 m at both sites. We suggest that there is basal melting at Dome C (0.56±0.19 mm/yr). The reconstructed velocity profile is highly non linear at both sites, which suggests complex ice flow effects. This induces a non linear thinning function in both drilling sites, which is also characterized by bumps corresponding to variations in surface elevation with time. [ABSTRACT FROM AUTHOR]

The study of the distribution of the crystallographic orientations (the fabric) along ice cores supplies information on the past and current ice flows of ice-sheets. Beside the usually observed formation of a vertical single maximum fabric, the EPICA Dome Concordia ice core (EDC) shows an abrupt and unexpected strengthening of its fabric during termination II around 1750 m depth. Such strengthenings were already observed for sites located on an ice-sheet. This suggests that horizontal shear could occur along the EDC core. Moreover, the change in the fabric leads to a modification of the viscosity between neighbouring ice layers. Through the use of an anisotropic ice flow model, we quantify the change in viscosity and investigate its implication on ice flow and dating. [ABSTRACT FROM AUTHOR]

Ice cores provide unique archives of past climate and environmental changes based only on physical processes. Quantitative temperature reconstructions are essential for the comparison between ice core records and climate models. Several methods have been developed to reconstruct past local temperatures from deep ice cores. Here we first analyse the long term fluctuations of temperature as depicted in the long Antarctic record from EPICA Dome C. The long term imprint of obliquity changes in the EPICA Dome C record is highlighted and compared to simulations conducted with the ECBILT-CLIO intermediate complexity climate model. We discuss the comparison between the current interglacial period and the long interglacial corresponding to marine isotopic stage 11, ∼400 kyr BP. Previous studies had focused on the role of precession and the thresholds required to induce glacial inceptions. We suggest that, due to the low eccentricity configuration of MIS 11 and the Holocene, the effect of precession on the incoming solar radiation is damped and that changes in obliquity must be taken into account. The EPICA Dome C alignment of terminations I and VI published in 2004 corresponds to a phasing of the obliquity signals. A conjunction of low obliquity and minimum northern hemisphere summer insolation is not found in the next tens of thousand years, supporting the idea of an unusually long interglacial ahead. As a second point relevant for future climate change, we discuss the magnitude and rate of change of past temperatures reconstructed from Greenland (NorthGRIP) and Antarctic (Dome C) ice cores. Past episodes of temperatures above the present-day values by up to 5°C are recorded at both locations during the penultimate interglacial period. The polar warming simulated by coupled climate models forced by a CO2increase of 1% per year is compared to ice-core-based temperature reconstructions. In Antarctica, the CO2-induced warming lies clearly beyond the natural rhythm of temperature fluctuations. In Greenland, the CO2-induced warming is as fast or faster than the most rapid temperature shifts of the last ice age. The magnitude of polar temperature change in response to a quadrupling of atmospheric CO2 is comparable to the magnitude of the polar temperature change from the Last Glacial Maximum to present-day. When forced by prescribed changes in ice sheet reconstructions and CO2 changes, climate models systematically underestimate the glacial-interglacial polar temperature change. [ABSTRACT FROM AUTHOR]