Interhemispheric coupling, the West Antarctic Ice Sheet and warm Antarctic interglacials

Ice core evidence indicates that even though atmospheric CO2 concentrations did not exceed ~300 ppm at any point during the last 800 000 years, East Antarctica was at least ~3–4 °C warmer than preindustrial (CO2~280 ppm) in each of the last four interglacials. During the previous three interglacials, this anomalous warming was short lived (~3000 years) and apparently occurred before the completion of Northern Hemisphere deglaciation. Hereafter, we refer to these periods as "Warmer than Present Transients" (WPTs). We present a series of experiments to investigate the impact of deglacial meltwater on the Atlantic Meridional Overturning Circulation (AMOC) and Antarctic temperature. It is well known that a slowed AMOC would increase southern sea surface temperature (SST) through the bipolar seesaw and observational data suggests that the AMOC remained weak throughout the terminations preceding WPTs, strengthening rapidly at a time which coincides closely with peak Antarctic temperature. We present two 800 kyr transient simulations using the Intermediate Complexity model GENIE-1 which demonstrate that meltwater forcing generates transient southern warming that is consistent with the timing of WPTs, but is not sufficient (in this single parameterisation) to reproduce the magnitude of observed warmth. In order to investigate model and boundary condition uncertainty, we present three ensembles of transient GENIE-1 simulations across Termination II (135 000 to 124 000 BP) and three snapshot HadCM3 simulations at 130 000 BP. Only with consideration of the possible feedback of West Antarctic Ice Sheet (WAIS) retreat does it become possible to simulate the magnitude of observed warming.

Direct linking of Greenland and Antarctic ice cores at the Toba eruption (74 ka BP)

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 palaeo-environmental records other than ice cores into a precise climatic context.

An optimized multi-proxy, multi-site Antarctic ice and gas orbital chronology (AICC2012): 120--800 ka

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 δ18Oatm over Marine Isotope Stage (MIS) 11–12 on EDC and a complete δ18Oatm record of the TALDICE ice cores permit us to derive additional orbital gas age constraints. The coherency of the different orbitally deduced ages (from δ18Oatm, δO2/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.

A community-based geological reconstruction of Antarctic Ice Sheet deglaciation since the Last Glacial Maximum

A robust understanding of Antarctic Ice Sheet deglacial history since the Last Glacial Maximum is important in order to constrain ice sheet and glacial-isostatic adjustment models, and to explore the forcing mechanisms responsible for ice sheet retreat. Such understanding can be derived from a broad range of geological and glaciological datasets and recent decades have seen an upsurge in such data gathering around the continent and Sub-Antarctic islands. Here, we report a new synthesis of those datasets, based on an accompanying series of reviews of the geological data, organised by sector. We present a series of timeslice maps for 20 ka, 15 ka, 10 ka and 5 ka, including grounding line position and ice sheet thickness changes, along with a clear assessment of levels of confidence. The reconstruction shows that the Antarctic Ice sheet did not everywhere reach the continental shelf edge at its maximum, that initial retreat was asynchronous, and that the spatial pattern of deglaciation was highly variable, particularly on the inner shelf. The deglacial reconstruction is consistent with a moderate overall excess ice volume and with a relatively small Antarctic contribution to meltwater pulse 1a. We discuss key areas of uncertainty both around the continent and by time interval, and we highlight potential priorities for future work. The synthesis is intended to be a resource for the modelling and glacial geological community.

Phase relationships between orbital forcing and the composition of air trapped in Antarctic ice cores

Orbital tuning is central for ice core chronologies beyond annual layer counting, available back to 60 ka (i.e. thousands of years before 1950) for Greenland ice cores. While several complementary orbital tuning tools have recently been developed using δ¹⁸Oatm, δO₂⁄N₂ and air content with different orbital targets, quantifying their uncertainties remains a challenge. Indeed, the exact processes linking variations of these parameters, measured in the air trapped in ice, to their orbital targets are not yet fully understood. Here, we provide new series of δO₂∕N₂ and δ¹⁸Oatm data encompassing Marine Isotopic Stage (MIS) 5 (between 100 and 160 ka) and the oldest part (340–800 ka) of the East Antarctic EPICA Dome C (EDC) ice core. For the first time, the measurements over MIS 5 allow an inter-comparison of δO₂∕N₂ and δ¹⁸Oatm records from three East Antarctic ice core sites (EDC, Vostok and Dome F). This comparison highlights some site-specific δO₂∕N₂ variations. Such an observation, the evidence of a 100 ka periodicity in the δO₂∕N₂ signal and the difficulty to identify extrema and mid-slopes in δO2∕N2 increase the uncertainty associated with the use of δO₂∕N₂ as an orbital tuning tool, now calculated to be 3–4 ka. When combining records of δ¹⁸Oatm and δO₂∕N₂ from Vostok and EDC, we find a loss of orbital signature for these two parameters during periods of minimum eccentricity (∼ 400 ka, ∼ 720–800 ka). Our data set reveals a time-varying offset between δO₂∕N₂ and δ¹⁸Oatm records over the last 800 ka that we interpret as variations in the lagged response of δ¹⁸Oatm to precession. The largest offsets are identified during Terminations II, MIS 8 and MIS 16, corresponding to periods of destabilization of the Northern polar ice sheets. We therefore suggest that the occurrence of Heinrich–like events influences the response of δ¹⁸Oatm to precession.

Glacial–interglacial dynamics of Antarctic firn columns: comparison between simulations and ice core air-δ15N measurements

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.

Duration of Greenland Stadial 22 and ice-gas Δage from counting of annual layers in Greenland NGRIP ice core

High-resolution measurements of chemical impurities and methane concentrations in Greenland ice core samples from the early glacial period allow the extension of annual-layer counted chronologies and the improvement of gas age-ice age difference (Δage) essential to the synchronization of ice core records. We report high-resolution measurements of a 50 m section of the NorthGRIP ice core and corresponding annual layer thicknesses in order to constrain the duration of the Greenland Stadial 22 (GS-22) between Greenland Interstadials (GIs) 21 and 22, for which inconsistent durations and ages have been reported from Greenland and Antarctic ice core records as well as European speleothems. Depending on the chronology used, GS-22 occurred between approximately 89 (end of GI-22) and 83 kyr b2k (onset of GI-21). From annual layer counting, we find that GS-22 lasted between 2696 and 3092 years and was followed by a GI-21 pre-cursor event lasting between 331 and 369 yr. Our layer-based counting agrees with the duration of stadial 22 as determined from the NALPS speleothem record (3250 ± 526 yr) but not with that of the GICC05modelext chronology (2620 yr) or an alternative chronology based on gas-marker synchronization to EPICA Dronning Maud Land ice core. These results show that GICC05modelext overestimates accumulation and/or underestimates thinning in this early part of the last glacial period. We also revise the possible ranges of NorthGRIP Δdepth (5.49 to 5.85 m) and Δage (498 to 601 yr) at the warming onset of GI-21 as well as the Δage range at the onset of the GI-21 precursor warming (523 to 654 yr), observing that temperature (represented by the δ15N proxy) increases before CH4 concentration by no more than a few decades.