EPICA Dome C Research Articles

Local summer temperature changes over the past 440 ka revealed by the total air content in the Antarctic EPICA Dome C ice core

Abstract. Seasonal temperature reconstructions from ice cores are missing over glacial–interglacial timescales, preventing a good understanding of the driving factors of Antarctic past climate changes. Here the total air content (TAC) record from the Antarctic EPICA Dome C (EDC) ice core is analyzed over the last 440 ka (thousand years). While the water isotopic record, a tracer for annual mean surface temperature, exhibits a dominant ∼100 kyr cyclicity, the TAC record is associated with a dominant ∼40 kyr cyclicity. Our results show that the TAC record is anti-correlated with the mean insolation over the local astronomical summer half-year. They also show for the first time that it is highly anti-correlated with local summer temperature simulated with an Earth system model of intermediate complexity. We propose that (1) the local summer insolation controls the local summer temperature; (2) the latter, through the development of temperature gradients at the near-surface of the ice sheet (<2 m), is affecting the surface snow structure; and (3) those snow structure changes propagating down to the bottom of the firn through densification are eventually controlling the pore volume at the bubble close-off and consequently the TAC. Hence, our results suggest that the EDC TAC record could be used as a proxy for local summer temperature changes. Also, our new simulations show that the mean insolation over the local astronomical summer half-year is the primary driver of Antarctic summer surface temperature variations, while changes in atmospheric greenhouse gas (GHG) concentrations and Northern Hemisphere (NH) ice sheet configurations play a more important role in Antarctic annual surface temperature changes.

Stagnant ice and age modelling in the Dome C region, Antarctica

Abstract. The European Beyond EPICA project aims to extract a continuous ice core of up to 1.5 Ma, with a maximum age density of 20 kyr m−1 at Little Dome C (LDC). We present a 1D numerical model which calculates the age of the ice around Dome C. The model inverts for basal conditions and accounts either for melting or for a layer of stagnant ice above the bedrock. It is constrained by internal reflecting horizons traced in radargrams and dated using the EPICA Dome C (EDC) ice core age profile. We used three different radar datasets ranging from a 10 000 km2 airborne survey down to 5 km long ground-based radar transects over LDC. We find that stagnant ice exists in many places, including above the LDC relief where the new Beyond EPICA drill site (BELDC) is located. The modelled thickness of this layer of stagnant ice roughly corresponds to the thickness of the basal unit observed in one of the radar surveys and in the autonomous phase-sensitive radio-echo sounder (ApRES) dataset. At BELDC, the modelled stagnant ice thickness is 198±44 m and the modelled oldest age of ice is 1.45±0.16 Ma at a depth of 2494±30 m. This is very similar to all sites situated on the LDC relief, including that of the Million Year Ice Core project being conducted by the Australian Antarctic Division. The model was also applied to radar data in the area 10–15 km north of EDC (North Patch), where we find either a thin layer of stagnant ice (generally <60 m) or a negligible melt rate (<0.1 mm yr−1). The modelled maximum age at North Patch is over 2 Ma in most places, with ice at 1.5 Ma having a resolution of 9–12 kyr m−1, making it an exciting prospect for a future Oldest Ice drill site.

Effects of Last Glacial Maximum (LGM) sea surface temperature and sea ice extent on the isotope–temperature slope at polar ice core sites

Abstract. Stable water isotopes in polar ice cores are widely used to reconstruct past temperature variations over several orbital climatic cycles. One way to calibrate the isotope–temperature relationship is to apply the present-day spatial relationship as a surrogate for the temporal one. However, this method leads to large uncertainties because several factors like the sea surface conditions or the origin and transport of water vapor influence the isotope–temperature temporal slope. In this study, we investigate how the sea surface temperature (SST), the sea ice extent, and the strength of the Atlantic Meridional Overturning Circulation (AMOC) affect these temporal slopes in Greenland and Antarctica for Last Glacial Maximum (LGM, ∼ 21 000 years ago) to preindustrial climate change. For that, we use the isotope-enabled atmosphere climate model ECHAM6-wiso, forced with a set of sea surface boundary condition datasets based on reconstructions (e.g., GLOMAP) or MIROC 4m simulation outputs. We found that the isotope–temperature temporal slopes in East Antarctic coastal areas are mainly controlled by the sea ice extent, while the sea surface temperature cooling affects the temporal slope values inland more. On the other hand, ECHAM6-wiso simulates the impact of sea ice extent on the EPICA Dome C (EDC) and Vostok sites through the contribution of water vapor from lower latitudes. Effects of sea surface boundary condition changes on modeled isotope–temperature temporal slopes are variable in West Antarctica. This is partly due to the transport of water vapor from the Southern Ocean to this area that can dampen the influence of local temperature on the changes in the isotopic composition of precipitation and snow. In the Greenland area, the isotope–temperature temporal slopes are influenced by the sea surface temperatures near the coasts of the continent. The greater the LGM cooling off the coast of southeastern Greenland, the greater the transport of water vapor from the North Atlantic, and the larger the temporal slopes. The presence or absence of sea ice very near the coast has a large influence in Baffin Bay and the Greenland Sea and influences the slopes at some inland ice core stations. The extent of the sea ice far south slightly influences the temporal slopes in Greenland through the transport of more depleted water vapor from lower latitudes to this area. The seasonal variations of sea ice distribution, especially its retreat in summer, influence the isotopic composition of the water vapor in this region and the modeled isotope–temperature temporal slopes in the eastern part of Greenland. A stronger LGM AMOC decreases LGM-to-preindustrial isotopic anomalies in precipitation in Greenland, degrading the isotopic model–data agreement. The AMOC strength modifies the temporal slopes over inner Greenland slightly and by a little on the coasts along the Greenland Sea where the changes in surface temperature and sea ice distribution due to the AMOC strength mainly occur.

Sub-millennial climate variability from high-resolution water isotopes in the EPICA Dome C ice core

Abstract. The EPICA Dome C (EDC) ice core provides the longest continuous climatic record, covering the last 800 000 years (800 kyr). A unique opportunity to investigate decadal to millennial variability during past glacial and interglacial periods is provided by the high-resolution water isotopic record (δ18O and δD) available for the EDC ice core. We present here a continuous compilation of the EDC water isotopic record at a sample resolution of 11 cm, which consists of 27 000 δ18O measurements and 7920 δD measurements (covering, respectively, 94 % and 27 % of the whole EDC record), including published and new measurements (2900 for both δ18O and δD) for the last 800 kyr. Here, we demonstrate that repeated water isotope measurements of the same EDC samples from different depth intervals obtained using different analytical methods are comparable within analytical uncertainty. We thus combine all available EDC water isotope measurements to generate a high-resolution (11 cm) dataset for the past 800 kyr. A frequency decomposition of the most complete δ18O record and a simple assessment of the possible influence of diffusion on the measured profile shows that the variability at the multi-decadal to multi-centennial timescale is higher during glacial than during interglacial periods and higher during early interglacial isotopic maxima than during the Holocene. This analysis shows as well that during interglacial periods characterized by a temperature optimum at the beginning, the multi-centennial variability is strongest over this temperature optimum.

Insight Into Provenance and Variability of Atmospheric Dust in Antarctic Ice Cores During the Late Pleistocene From Magnetic Measurements

We measured saturation isothermal remanent magnetization (SIRM), coercivity of remanence (Hcr) and insoluble dust concentration (IDC) of 52 ice samples from Vostok and EPICA Dome-C ice cores (Antarctica) as a measure of magnetic properties of the aerosol dust entrapped in the ice. Samples range in age from marine isotopic stage (MIS) 7 to 19 in EPICA Dome-C ice core and from MIS 1 to 11 in Vostok ice core. Data from ice samples were compared with 86 samples from possible source areas (PSA) from East Antarctica, and 10 samples from South America and New Zealand. Previous results from MIS 1 to MIS 6 found that magnetic properties of aerosol dust could be divided in two distinct groups characterized by high-Hcr and low-SIRMdust for glacial samples, and low-Hcr and high-SIRMdust, for interglacial samples. The new data from older ice samples highlighted several discrepancies from this expectation with significant differences between Vostok and Dome-C sites. Magnetic properties of Antarctic PSA sample show a large variability, however PSA samples from Victoria Land and few other, have magnetic properties compatible with that of the glacial dust, or more precisely with samples characterized by high dust flux. The new data from Pleistocene ice and from PSA samples confirm the South American and Antarctic provenance of the largest atmospheric dust load typical of glacial stages. On the other hand, we did not found any PSA sample with properties compatible with the highly magnetic samples (mostly from interglacial stages), which are characterized by low IDC. These samples from the oldest and deepest part of the cores revealed a more complex picture than previously outlined from the analysis of MIS 1–6, and show unusual magnetic properties which can be tentatively attributed to chemical.

Compositions of Dust and Sea Salts in the Dome C and Dome Fuji Ice Cores From Last Glacial Maximum to Early Holocene Based on Ice‐Sublimation and Single‐Particle Measurements

AbstractWe analyzed the chemical compositions of dust and sea‐salt particles in the EPICA Dome C (EDC) ice core during 26–7 kyr BP using an ice‐sublimation technique and compared the results with existing data of the Dome Fuji (DF) ice core. Combined with ion concentration data, our data suggested similar sea‐salt fluxes in both cores and significantly lower dust flux in the EDC core. The differences in modal size and aspect ratio of dust particles between the two cores support the dominance of Patagonian source suggested by earlier works. The compositions of calcic dust showed major change at ~17 kyr BP, possibly reflecting a relative increase in dust transported via the upper troposphere. The calcium sulfate fraction was higher in the DF core than in the EDC core after ~17 kyr BP, suggesting that higher Patagonian dust contribution to the DF region. Abundant NaCl particles were found in the DF core in comparison with the EDC core from the LGM to early Holocene, possibly because of the high concentration of terrestrial dust in the DF core that reduced acid availability for sea‐salt modification. During the Holocene, the lower NaCl fraction and Cl−/Na+ ratio in the EDC core suggested that most Cl− was lost to the atmosphere from snow at Dome C, while it was preserved at Dome Fuji as NaCl and solid solution.

Unveiling the anatomy of Termination 3 using water and air isotopes in the Dome C ice core, East Antarctica

Each glacial – interglacial transition of the Quaternary occurs in a different orbital context leading to various timing for the deglaciation and sequence of high vs low latitudes events. Termination 3, 250 kyears before present (ka), is an unusual deglaciation in the context of the last 9 deglaciations recorded in the old EPICA Dome C (EDC) Antarctic ice core: it exhibits a three-phase sequence, two warming phases separated by a small cooling, the last phase suggesting a particularly rapid temperature increase. We present here new high resolution δ15N and deuterium excess (d-excess) data from the EDC ice core to provide a detailed temperature change estimate during this termination. Then, we combined the δD and δ18O to discuss the relationship between high and low latitude changes through the d-excess. We also provide the high vs low latitude sequence of events over this deglaciation without chronological uncertainty using low latitude ice core proxies. In agreement with previous studies based on speleothem analyses, we show that the first phase of Termination 3 (256–249 ka) is associated with small Heinrich like events linked to changes in ITCZ position, monsoon activity and teleconnections with Antarctica. In a context of minimum Northern Hemisphere insolation, this leads to a rather strong Antarctic warming, as observed in the δ15N record in contrast to the relatively small δD increase. The second warming phase occurs during the rise of the Northern hemisphere insolation, with a large Heinrich like event leading to the characteristic Antarctic warming observed in the δ15N and δD increase as for the more recent terminations.

On the use of δ18Oatm for ice core dating

Deep ice core chronologies have been improved over the past years through the addition of new age constraints. However, dating methods are still associated with large uncertainties for ice cores from the East Antarctic plateau where layer counting is not possible. Indeed, an uncertainty up to 6 ka is associated with AICC2012 chronology of EPICA Dome C (EDC) ice core, which mostly arises from uncertainty on the delay between changes recorded in δ18Oatm and in June 21st insolation variations at 65°N used for ice core orbital dating. Consequently, we need to enhance the knowledge of this delay to improve ice core chronologies.We present new high-resolution EDC δ18Oatm record (153–374 ka) and δO2/N2 measurements (163–332 ka) performed on well-stored ice to provide continuous records of δ18Oatm and δO2/N2 between 100 and 800 ka. The comparison of δ18Oatm with the δ18Ocalcite from East Asian speleothems shows that both signals present similar orbital and millennial variabilities, which may represent shifts in the InterTropical Convergence Zone position, themselves associated with Heinrich events. We thus propose to use the δ18Ocalcite as target for δ18Oatm orbital dating. Such a tuning method improves the ice core chronology of the last glacial inception compared to AICC2012 by reconciling NGRIP and mid-latitude climatic records. It is especially marked during Dansgaard-Oeschger 25 where the proposed chronology is 2.2 ka older than AICC2012. This δ18Oatm – δ18Ocalcite alignment method applied between 100 and 640 ka improves the EDC ice core chronology, especially over MIS 11, and leads to lower ice age uncertainties compared to AICC2012.

Is there 1.5-million-year-old ice near Dome C, Antarctica?

Abstract. Ice sheets provide exceptional archives of past changes in polar climate, regional environment and global atmospheric composition. The oldest dated deep ice core drilled in Antarctica has been retrieved at EPICA Dome C (EDC), reaching ∼ 800 000 years. Obtaining an older paleoclimatic record from Antarctica is one of the greatest challenges of the ice core community. Here, we use internal isochrones, identified from airborne radar coupled to ice-flow modelling to estimate the age of basal ice along transects in the Dome C area. Three glaciological properties are inferred from isochrones: surface accumulation rate, geothermal flux and the exponent of the Lliboutry velocity profile. We find that old ice (> 1.5 Myr, 1.5 million years) likely exists in two regions: one ∼ 40 km south-west of Dome C along the ice divide to Vostok, close to a secondary dome that we name Little Dome C (LDC), and a second region named North Patch (NP) located 10–30 km north-east of Dome C, in a region where the geothermal flux is apparently relatively low. Our work demonstrates the value of combining radar observations with ice flow modelling to accurately represent the true nature of ice flow, and understand the formation of ice-sheet architecture, in the centre of large ice sheets.

Comparison of measurements from different radio-echo sounding systems and synchronization with the ice core at Dome C, Antarctica

Abstract. We present a compilation of radio-echo sounding (RES) measurements of five radar systems (AWI, BAS, CReSIS, INGV and UTIG) around the EPICA Dome C (EDC) drill site, East Antarctica. The aim of our study is to investigate the differences of the various systems in their resolution of internal reflection horizons (IRHs) and bed topography, penetration depth and capacity of imaging the basal layer. We address the questions of the compatibility of existing radar data for common interpretation and the suitability of the individual systems for reconnaissance surveys. We find that the most distinct IRHs and IRH patterns can be identified and transferred between most data sets. Considerable differences between the RES systems exist in range resolution and depiction of the bottom-most region. Considering both aspects, which we judge as crucial factors in the search for old ice, the CReSIS and the UTIG systems are the most suitable ones. In addition to the RES data set comparison we calculate a synthetic radar trace from EDC density and conductivity profiles. We identify 10 common IRHs in the measured RES data and the synthetic trace. We then conduct a sensitivity study for which we remove certain peaks from the input conductivity profile. As a result the respective reflections disappear from the modeled radar trace. In this way, we establish a depth conversion of the measured travel times of the IRHs. Furthermore, we use these sensitivity studies to investigate the cause of observed reflections. The identified IRHs are assigned ages from the EDC's timescale. Due to the isochronous character of these conductivity-caused IRHs, they are a means to extend the Dome C age structure by tracing the IRHs along the RES profiles.

0.9 Ma oxygen isotope stratigraphy for a shallow-water sedimentary transect across three IODP 317 sites in the Canterbury Bight of the southwest Pacific Ocean

Sedimentary records in shallow-water environments provide a unique opportunity to further our understanding of regional relative sea level changes in relation to global climate change. Here we present a 0.9Ma oxygen isotope stratigraphy for a shallow-water sedimentary transect across three IODP 317 sites in the Canterbury Bight of the southwest Pacific Ocean. The sites, located on the eastern margin of the South Island of New Zealand, include a continental slope site (U1352) and two continental shelf sites (U1354 and U1351). We first generated benthic foraminifers (Nonionella flemingi) δ18O records for the three sites and a planktonic (Globigerina bulloides) record for U1352B. An initial chronological framework for the benthic δ18O record of U1352B was constructed using two accelerator mass spectrometry radiocarbon dates and four biostratigraphic events. Then a refined age model was established by correlating the U1352B benthic δ18O record with the EDC (EPICA Dome C) δD record on the AICC2012 time-scale and LR04 benthic δ18O stack. Although U1354B and U1351B have lower sedimentation rates, their benthic δ18O records correlate well with that of U1352B. To ensure the accuracy of the chronostratigraphic framework established, we also analyzed characteristics of sedimentary grain size and planktonic and benthic δ18O values. In accord with the adjacent sites, the results show that melting of the Southern Alps glaciers resulting from the warming climate during the Marine Isotope Stage (MIS) 11 and MIS 5.5 led to increased fresh water delivery, with massive terrigenous deposits. The warm sea surface temperature during MIS 7 is related to Subtropical Front migration, which led to strong current activity and coarsened grain size. Meanwhile, records of the benthic δ18O, sedimentation rate, and >63μm coarse fraction contents of Site U1352 all indicate that MIS 20 was indeed a colder interval compared to subsequent glacial times.

Climate dependent contrast in surface mass balance in East Antarctica over the past 216 ka

ABSTRACTDocumenting 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.

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

Abstract. 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 δ18Oatm, δO2⁄N2 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 δO2∕N2 and δ18Oatm 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 δO2∕N2 and δ18Oatm records from three East Antarctic ice core sites (EDC, Vostok and Dome F). This comparison highlights some site-specific δO2∕N2 variations. Such an observation, the evidence of a 100 ka periodicity in the δO2∕N2 signal and the difficulty to identify extrema and mid-slopes in δO2∕N2 increase the uncertainty associated with the use of δO2∕N2 as an orbital tuning tool, now calculated to be 3–4 ka. When combining records of δ18Oatm and δO2∕N2 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 δO2∕N2 and δ18Oatm records over the last 800 ka that we interpret as variations in the lagged response of δ18Oatm 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 δ18Oatm to precession.

Volcanic synchronization of Dome Fuji and Dome C Antarctic deep ice cores over the past 216 kyr

Abstract. 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.

No evidence for planetary influence on solar activity 330 000 years ago

Context. Abreu et al. (2012, A&A. 548, A88) have recently compared the periodicities in a 14 C – 10 Be proxy record of solar variability during the Holocene and found a strong similarity with the periodicities predicted on the basis of a model of the time-dependent torque exerted by the planets on the sun’s tachocline. If verified, this effect would represent a dramatic advance not only in the basic understanding of the Sun’s variable activity, but also in the potential influence of this variability on the Earth’s climate. Cameron and Schussler (2013, A&A. 557, A83) have seriously criticized the statistical treatment used by Abreu et al. to test the significance of the coincidences between the periodicities of their model with the Holocene proxy record.Aims. If the Abreu et al. hypothesis is correct, it should be possible to find the same periodicities in the records of cosmogenic nuclides at earlier times.Methods. We present here a high-resolution record of 10 Be in the EPICA Dome C (EDC) ice core from Antarctica during the Marine Interglacial Stage 9.3 (MIS 9.3), 325–336 kyr ago, and investigate its spectral properties.Results. We find very limited similarity with the periodicities seen in the proxy record of solar variability during the Holocene, or with that of the model of Abreu et al.Conclusions. We find no support for the hypothesis of a planetary influence on solar activity, and raise the question of whether the centennial periodicities of solar activity observed during the Holocene are representative of solar activity variability in general.

High-resolution mineral dust and sea ice proxy records from the Talos Dome ice core

Abstract. In this study we report on new non-sea salt calcium (nssCa2+, mineral dust proxy) and sea salt sodium (ssNa+, sea ice proxy) records along the East Antarctic Talos Dome deep ice core in centennial resolution reaching back 150 thousand years (ka) before present. During glacial conditions nssCa2+ fluxes in Talos Dome are strongly related to temperature as has been observed before in other deep Antarctic ice core records, and has been associated with synchronous changes in the main source region (southern South America) during climate variations in the last glacial. However, during warmer climate conditions Talos Dome mineral dust input is clearly elevated compared to other records mainly due to the contribution of additional local dust sources in the Ross Sea area. Based on a simple transport model, we compare nssCa2+ fluxes of different East Antarctic ice cores. From this multi-site comparison we conclude that changes in transport efficiency or atmospheric lifetime of dust particles do have a minor effect compared to source strength changes on the large-scale concentration changes observed in Antarctic ice cores during climate variations of the past 150 ka. Our transport model applied on ice core data is further validated by climate model data. The availability of multiple East Antarctic nssCa2+ records also allows for a revision of a former estimate on the atmospheric CO2 sensitivity to reduced dust induced iron fertilisation in the Southern Ocean during the transition from the Last Glacial Maximum to the Holocene (T1). While a former estimate based on the EPICA Dome C (EDC) record only suggested 20 ppm, we find that reduced dust induced iron fertilisation in the Southern Ocean may be responsible for up to 40 ppm of the total atmospheric CO2 increase during T1. During the last interglacial, ssNa+ levels of EDC and EPICA Dronning Maud Land (EDML) are only half of the Holocene levels, in line with higher temperatures during that period, indicating much reduced sea ice extent in the Atlantic as well as the Indian Ocean sector of the Southern Ocean. In contrast, Holocene ssNa+ flux in Talos Dome is about the same as during the last interglacial, indicating that there was similar ice cover present in the Ross Sea area during MIS 5.5 as during the Holocene.

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

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

The Antarctic ice core chronology (AICC2012): an optimized multi-parameter and multi-site dating approach for the last 120 thousand years

Abstract. 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 δ15N 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.

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

Abstract. 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.

A 250 ka oxygen isotope record from diatoms at Lake El'gygytgyn, far east Russian Arctic

Abstract. In 2003 sediment core Lz1024 was drilled at Lake El'gygytgyn, far east Russian Arctic, in an area of the Northern Hemisphere which has not been glaciated for the last 3.6 Ma. Biogenic silica was used for analysing the oxygen isotope composition (δ18Odiatom) in the upper 13 m long section dating back about 250 ka with samples dominated by one taxa in the <10 μm fraction (Cyclotella ocellata). Downcore variations in δ18O values show that glacial-interglacial cycles are present throughout the core and δ18Odiatom-values are mainly controlled by δ18Oprecipitation. Changes reflect the Holocene Thermal Maximum, the Last Glacial Maximum (LGM) and the interglacial periods corresponding to MIS 5.5 and MIS 7 with a peak-to-peak amplitude between LGM and MIS 5.5 of Δ18O = 5.3‰. This corresponds to a mean annual air temperature difference of about 9 °C. Our record is the first continuous δ18Odiatom record from an Arctic lake sediment core directly responding to precipitation and dating back more than 250 ka and correlates well with the stacked marine δ18O LR04 (r = 0.58) and δD EPICA Dome-C record (r = 0.69). With δ18O results indicating strong links to both marine and ice-core records, records from Lake El'gygytgyn can be used to further investigate the sensitivity of the Arctic climate to both past and future global climatic changes.