CO2 concentration and stable isotope ratios of three Antarctic ice cores covering the period from 149.4 - 1.5 kyr before 1950

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

Table 2. Lead and Ba concentrations and Pb isotopic compositions measured in Law Dome ice cores

Lead isotopic compositions and Pb and Ba concentrations have been measured in ice cores from Law Dome, East Antarctica, covering the past 6500 years. 'Natural' background concentrations of Pb (ab. 0.4 pg/g) and Ba (ab. 1.3 pg/g) are observed until 1884 AD, after which increased Pb concentrations and lowered 206Pb/207Pb ratios indicate the influence of anthropogenic Pb. The isotopic composition of 'natural' Pb varies within the range 206Pb/207Pb=1.20-1.25 and 208Pb/207Pb=2.46-2.50, with an average rock and soil dust Pb contribution of 8-12%. A major pollution event is observed at Law Dome between 1884 and 1908 AD, elevating the Pb concentration four-fold and changing 206Pb/207Pb ratios in the ice to ab. 1.12. Based on Pb isotopic systematics and Pb emission statistics, this is attributed to Pb mined at Broken Hill and smelted at Broken Hill and Port Pirie, Australia. Anthropogenic Pb inputs are at their greatest from 1900 to 1910 and from ab. 1960 to ab. 1980. During the 20th century, Ba concentrations are consistently higher than 'natural' levels and are attributed to increased dust production, suggesting the influence of climate change and/or changes in land coverage with vegetation.

Stable carbon isotope ratios of carbon dioxide from EDC and Berkner Island ice cores for 40-50 ka BP

The stable carbon isotopic signature of carbon dioxide (d13CO2) measured in the air occlusions of polar ice provides important constraints on the carbon cycle in past climates. In order to exploit this information for previous glacial periods, one must use deep, clathrated ice, where the occluded air is preserved not in bubbles but in the form of air hydrates. Therefore, it must be established whether the original atmospheric d13CO2 signature can be reconstructed from clathrated ice. We present a comparative study using coeval bubbly ice from Berkner Island and ice from the bubble-clathrate transformation zone (BCTZ) of EPICA Dome C (EDC). In the EDC samples the gas is partitioned into clathrates and remaining bubbles as shown by erroneously low and scattered CO2 concentration values, presenting a worst-case test for d13CO2 reconstructions. Even so, the reconstructed atmospheric d13CO2 values show only slightly larger scatter. The difference to data from coeval bubbly ice is statistically significant. However, the 0.16 per mil magnitude of the offset is small for practical purposes, especially in light of uncertainty from non-uniform corrections for diffusion related fractionation that could contribute to the discrepancy. Our results are promising for palaeo-atmospheric studies of d13CO2 using a ball mill dry extraction technique below the BCTZ of ice cores, where gas is not subject to fractionation into microfractures and between clathrate and bubble reservoirs.

The consolidation and bond strength of rafted sea ice

The consolidation and bond strength of rafted sea ice were investigated through a series of experiments undertaken in the Ice Physics Laboratory at the UCL. To simulate a section of rafted sea ice, blocks of laboratory grown saline ice were stacked in an insulated tank with spacers between adjacent blocks to allow saline water to flood in. The rate of consolidation was then monitored using a combination of temperature readings recorded in the ice and liquid layer, salinity measurements of the liquid layer, and cores taken at specific times of interest. Two states of consolidation were observed: thermodynamic consolidation where the ice blocks were physically bonded but the bond strength was weak, and mechanical consolidation where the bond had reached full strength. Results showed that the rafted ice had physically bonded in less than a day, however it took many more days (6 to 30 depending on the environmental conditions) for the bond to reach maximum strength. Increasing the thickness of the ice, the salinity of the water and the inter-block gap size all increased the consolidation time. Once consolidated, ice cores were taken and sheared using the asymmetric four-point bending method to measure the strength of the bond between the ice blocks. These were then compared to the shear strength of solid ice blocks simulating level sea ice. Our results show that the shear strength of the bond between the rafted ice blocks is about 30% weaker than that of level ice.

A refined TALDICE-1a age scale from 55 to 112 ka before present for the Talos Dome ice core based on high-resolution methane measurements

A precise synchronization of different climate records is indispensable for a correct dynamical interpretation of paleoclimatic data. A chronology for the TALDICE ice core from the Ross Sea sector of East Antarctica has recently been presented based on methane synchronization with Greenland and the EDC ice cores and δ18Oice synchronization with EDC in the bottom part (TALDICE-1). Using new high-resolution methane data obtained with a continuous flow analysis technique, we present a refined age scale for the age interval from 55–112 thousand years (ka) before present, where TALDICE is synchronized with EDC. New and more precise tie points reduce the uncertainties of the age scale from up to 1900 yr in TALDICE-1 to below 1100 yr over most of the refined interval and shift the Talos Dome dating to significantly younger ages during the onset of Marine Isotope Stage 3. Thus, discussions of climate dynamics at sub-millennial time scales are now possible back to 110 ka, in particular during the inception of the last ice age. Calcium data of EDC and TALDICE are compared to show the impact of the refinement to the synchronization of the two ice cores not only for the gas but also for the ice age scale.

Simultaneous stable isotope analysis of methane and nitrous oxide on ice core samples

Methane and nitrous oxide are important greenhouse gases which show a strong increase in atmospheric mixing ratios since pre-industrial time as well as large variations during past climate changes. The understanding of their biogeochemical cycles can be improved using stable isotope analysis. However, high-precision isotope measurements on air trapped in ice cores are challenging because of the high susceptibility to contamination and fractionation. Here, we present a dry extraction system for combined CH4 and N2O stable isotope analysis from ice core air, using an ice grating device. The system allows simultaneous analysis of δD(CH4) or δ13C(CH4), together with δ15N(N2O), δ18O(N2O) and δ15N(NO+ fragment) on a single ice core sample, using two isotope mass spectrometry systems. The optimum quantity of ice for analysis is about 600 g with typical "Holocene" mixing ratios for CH4 and N2O. In this case, the reproducibility (1σ ) is 2.1‰ for δD(CH4), 0.18‰ for δ13C(CH4), 0.51‰ for δ15N(N2O), 0.69‰ for δ18O(N2O) and 1.12‰ for δ15N(NO+ fragment). For smaller amounts of ice the standard deviation increases, particularly for N2O isotopologues. For both gases, small-scale intercalibrations using air and/or ice samples have been carried out in collaboration with other institutes that are currently involved in isotope measurements of ice core air. Significant differences are shown between the calibration scales, but those offsets are consistent and can therefore be corrected for.

Simulating the temperature and precipitation signal in an Alpine ice core

Accumulation and delta O-18 data from Alpine ice cores provide information on past temperature and precipitation. However, their correlation with seasonal or annual mean temperature and precipitation at nearby sites is often low. This is partly due to the irregular sampling of the atmosphere by the ice core (i.e. ice cores almost only record precipitation events and not dry periods) and the possible incongruity between annual layers and calendar years. Using daily meteorological data from a nearby station and reanalyses, we replicate the ice core from the Grenzgletscher (Switzerland, 4200m a.s.l.) on a sample-by-sample basis by calculating precipitation-weighted temperature (PWT) over short intervals. Over the last 15 yr of the ice core record, accumulation and delta O-18 variations can be well reproduced on a sub-seasonal scale. This allows a wiggle-matching approach for defining quasi-annual layers, resulting in high correlations between measured quasi-annual delta O-18 and PWT. Further back in time, the agreement deteriorates. Nevertheless, we find significant correlations over the entire length of the record (1938-1993) of ice core delta O-18 with PWT, but not with annual mean temperature. This is due to the low correlations between PWT and annual mean temperature, a characteristic which in ERA-Interim reanalysis is also found for many other continental mid-to-high-latitude regions. The fact that meteorologically very different years can lead to similar combinations of PWT and accumulation poses limitations to the use of delta O-18 from Alpine ice cores for temperature reconstructions. Rather than for reconstructing annual mean temperature, delta O-18 from Alpine ice cores should be used to reconstruct PWT over quasi-annual periods. This variable is reproducible in reanalysis or climate model data and could thus be assimilated into conventional climate models.

Microgram-Level Radiocarbon Determination of Carbonaceous Particles in Firn and Ice Samples: Pretreatment and Oc/Ec Separation

Carbonaceous particles that comprise organic carbon (OC) and elemental carbon (EC) are of increasing interest in climate research because of their influence on the radiation balance of the Earth. The radiocarbon determination of particulate OC and EC extracted from ice cores provides a powerful tool to reconstruct the long-term natural and anthropogenic emissions of carbonaceous particles. However, this C-14-based source apportionment method has not been applied for the firn section, which is the uppermost part of Alpine glaciers with a typical thickness of up to 50 m. In contrast to glacier ice, firn samples are more easily contaminated through drilling and handling operations. In this study, an alternative decontamination method for firn samples consisting of chiselling off the outer parts instead of rinsing them was developed and verified. The obtained procedural blank of 2.8 +/- 0.8 mu g C for OC is a factor of 2 higher compared to the rinsing method used for ice, but still relatively low compared to the typical OC concentration in firn samples from Alpine glaciers. The EC blank of 0.3 +/- 0.1 mu g C is similar for both methods. For separation of OC and EC for subsequent C-14 analysis, a thermal-optical method instead of the purely thermal method was applied for the first time to firn and ice samples, resulting in a reduced uncertainty of both the mass and C-14 determination. OC and EC concentrations as well as their corresponding fraction of modern for firn and ice samples from Fiescherhorn and Jungfraujoch agree well with published results, validating the new method.

C-14 Measurements of Ice Samples from the Juvfonne Ice Tunnel, Jotunheimen, Southern Norway-Validation of a C-14 Dating Technique for Glacier Ice

Establishing precise age-depth relationships of high-alpine ice cores is essential in order to deduce conclusive paleoclimatic information from these archives. Radiocarbon dating of carbonaceous aerosol particles incorporated in such glaciers is a promising tool to gain absolute ages, especially from the deepest parts where conventional methods are commonly inapplicable. In this study, we present a new validation for a published C-14 dating method for ice cores. Previously C-14-dated horizons of organic material from the Juvfonne ice patch in central southern Norway (61.676 degrees N, 8.354 degrees E) were used as reference dates for adjacent ice layers, which were C-14 dated based on their particulate organic carbon (POC) fraction. Multiple measurements were carried out on 3 sampling locations within the ice patch featuring modern to multimillennial ice. The ages obtained from the analyzed samples were in agreement with the given age estimates. In addition to previous validation work, this independent verification gives further confidence that the investigated method provides the actual age of the ice.

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.

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.

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

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.

Spatial gradients of temperature, accumulation and delta18O-ice in Greenland over a series of Dansgaard-Oeschger events

Air and water stable isotope measurements from four Greenland deep ice cores (GRIP, GISP2, NGRIP and NEEM) are investigated over a series of Dansgaard–Oeschger events (DO 8, 9 and 10), which are representative of glacial millennial scale variability. Combined with firn modeling, air isotope data allow us to quantify abrupt temperature increases for each drill site (1σ = 0.6 °C for NEEM, GRIP and GISP2, 1.5 °C for NGRIP). Our data show that the magnitude of stadial–interstadial temperature increase is up to 2 °C larger in central and North Greenland than in northwest Greenland: i.e., for DO 8, a magnitude of +8.8 °C is inferred, which is significantly smaller than the +11.1 °C inferred at GISP2. The same spatial pattern is seen for accumulation increases. This pattern is coherent with climate simulations in response to reduced sea-ice extent in the Nordic seas. The temporal water isotope (δ18O)–temperature relationship varies between 0.3 and 0.6 (±0.08) ‰ °C−1 and is systematically larger at NEEM, possibly due to limited changes in precipitation seasonality compared to GISP2, GRIP or NGRIP. The gas age−ice age difference of warming events represented in water and air isotopes can only be modeled when assuming a 26% (NGRIP) to 40% (GRIP) lower accumulation than that derived from a Dansgaard–Johnsen ice flow model.

Topographic Control of Regional Accumulation Rate Variability at South Pole and Implications for Ice-Core Interpretation

Snow-accumulation rates are known to be sensitive to local changes in ice-sheet surface slope because of the effect of katabatic winds. These topographic effects can be preserved in ice cores that are collected at non-ice-divide locations. The trajectory of an ice-core site at South Pole is reconstructed using measurements of ice-sheet motion to show that snow was probably deposited at places of different surface slope during the past 1000 years. Recent accumulation rates, derived from shallow firn cores, vary along this trajectory according to surface topography, so that on a relatively steep flank mean annual accumulation is similar to 18% smaller than on a nearby topographic depression. These modern accumulation rates are used to reinterpret the cause of accumulation rate variability with time in the long ice-core record as an ice-dynamics effect and not a climate-change signal. The results highlight the importance of conducting ancillary ice-dynamics measurements as part of ice-coring programs so that topographic effects can be deconvolved from potential climate signals.

Summer Temperature Trend Over the Past Two Millennia Using Air Content in Himalayan Ice

Two Himalayan ice cores display a factor-two decreasing trend of air content over the past two millennia, in contrast to the relatively stable values in Greenland and Antarctica ice cores over the same period. Because the air content can be related with the relative frequency and intensity of melt phenomena, its variations along the Himalayan ice cores provide an indication of summer temperature trend. Our reconstruction point toward an unprecedented warming trend in the 20th century but does not depict the usual trends associated with "Medieval Warm Period" (MWP), or "Little Ice Age" (LIA).