Search

We investigate the lock-in zone (LIZ) of the East Greenland Ice Core Project (EastGRIP) region, northeastern Greenland, in detail. We present results from the firn air-pumping campaign of the S6 borehole, forward modeling, and a novel technique for finding the lock-in depth (LID, the top of the LIZ) based on the visual stratigraphy of the EastGRIP ice core. The findings in this work help to deepen our knowledge of how atmospheric gases are trapped in ice cores. CO2, δ15N, and CH4 data suggest that the LID lies around 58 to 61 m depth. With the pixel value intensity and bright-spot analysis based on visual stratigraphy, we can pinpoint a change in ice properties to exactly 58.3 m depth, which we define as the optical lock-in depth (OLID). This visual change in ice properties is caused by the formation of rounded and enclosed air bubbles that alter the measured refraction of the light pathways. The results for the LID and OLID agree accurately on the depth. We furthermore use the visual stratigraphy images to obtain information on the sharpness of the open- to closed-porosity transition. Combining traditional methods with the independent optical method presented here, we can now better constrain the bubble closure processes in polar firn.

While the processes controlling pore closure are broadly understood, the physical mechanisms driving the associated elemental fractionation remains ambiguous. Previous studies have shown that the pore closure process leads to a depletion in small-sized molecules (e.g. H2, O2, Ar, Ne, He) in ice core bubbles relative to larger-sized molecules like N2. This size-dependent fractionation, identified using ice core δ(O2/N2) records, exhibits a clear anti-correlation with local summer solstice insolation, making δ(O2/N2) a valuable ice core dating tool. Mechanisms controlling this relationship are attributed to the physical properties of deep firn. In this study, we compile δ(O2/N2) records from 15 polar ice cores and show a new additional link between δ(O2/N2) and local surface temperature and/or accumulation rate. Using the Crocus snowpack model, we perform sensitivity tests to identify the response of near-surface snow properties to changes in insolation intensity, accumulation rate, and air temperature. These tests support a mechanism linked to firn grain size, such that the larger the grain size for a given density, the stronger the pore closure fractionation and, hence, the lower the δ(O2/N2) values archived in the ice. Based on both snowpack model outputs and data compilation, our findings suggest that local accumulation rate and temperature should be considered when interpreting δ(O2/N2) as a local insolation proxy.

The EPICA (European Project for Ice Coring in Antarctica) Dome C (EDC) ice core drilling in East Antarctica reaches a depth of 3260 m. The reference EDC chronology, the AICC2012 (Antarctic Ice Core Chronology 2012), provides an age vs. depth relationship covering the last 800 kyr (thousands of years), with an absolute uncertainty rising up to 8000 years at the bottom of the ice core. The origins of this relatively large uncertainty are twofold: (1) the δ18Oatm, δO2/N2 and total air content (TAC) records are poorly resolved and show large gaps over the last 800 kyr, and (2) large uncertainties are associated with their orbital targets. Here, we present new highly resolved δ18Oatm, δO2/N2 and δ15N measurements for the EDC ice core covering the last five glacial–interglacial transitions; a new low-resolution TAC record over the period 440–800 ka BP (ka: 1000 years before 1950); and novel absolute 81Kr ages. We have compiled chronological and glaciological information including novel orbital age markers from new data on the EDC ice core as well as accurate firn modeling estimates in a Bayesian dating tool to construct the new AICC2023 chronology. For the first time, three orbital tools are used simultaneously. Hence, it is possible to observe that they are consistent with each other and with the other age markers over most of the last 800 kyr (70 %). This, in turn, gives us confidence in the new AICC2023 chronology. The average uncertainty in the ice chronology is reduced from 1700 to 900 years in AICC2023 over the last 800 kyr (1σ). The new timescale diverges from AICC2012 and suggests age shifts reaching 3800 years towards older ages over marine isotope stages (MISs) 5, 11 and 19. But the coherency between the new AICC2023 timescale and independent chronologies of other archives (Italian Lacustrine succession from Sulmona Basin, Dome Fuji ice core and northern Alpine speleothems) is improved by 1000 to 2000 years over these time intervals.

Carbon monoxide (CO) is a naturally occurring atmospheric trace gas, a regulated pollutant, and one of the main components determining the oxidative capacity of the atmosphere. Evaluating climate–chemistry models under different conditions than today and constraining past CO sources requires a reliable record of atmospheric CO mixing ratios ([CO]) that includes data since preindustrial times. Here, we report the first continuous record of atmospheric [CO] for Southern Hemisphere (SH) high latitudes over the past 3 millennia. Our continuous record is a composite of three high-resolution Antarctic ice core gas records and firn air measurements from seven Antarctic locations. The ice core gas [CO] records were measured by continuous flow analysis (CFA), using an optical feedback cavity-enhanced absorption spectrometer (OF-CEAS), achieving excellent external precision (2.8–8.8 ppb; 2σ) and consistently low blanks (ranging from 4.1±1.2 to 7.4±1.4 ppb), thus enabling paleo-atmospheric interpretations. Six new firn air [CO] Antarctic datasets collected between 1993 and 2016 CE at the DE08-2, DSSW19K, DSSW20K, South Pole, Aurora Basin North (ABN), and Lock-In sites (and one previously published firn CO dataset at Berkner) were used to reconstruct the atmospheric history of CO from ∼1897 CE, using inverse modeling that incorporates the influence of gas transport in firn. Excellent consistency was observed between the youngest ice core gas [CO] and the [CO] from the base of the firn and between the recent firn [CO] and atmospheric [CO] measurements at Mawson station (eastern Antarctica), yielding a consistent and contiguous record of CO across these different archives. Our Antarctic [CO] record is relatively stable from −835 to 1500 CE, with mixing ratios within a 30–45 ppb range (2σ). There is a ∼5 ppb decrease in [CO] to a minimum at around 1700 CE during the Little Ice Age. CO mixing ratios then increase over time to reach a maximum of ∼54 ppb by ∼1985 CE. Most of the industrial period [CO] growth occurred between about 1940 to 1985 CE, after which there was an overall [CO] decrease, as observed in Greenland firn air and later at atmospheric monitoring sites and attributed partly to reduced CO emissions from combustion sources. Our Antarctic ice core gas CO observations differ from previously published records in two key aspects. First, our mixing ratios are significantly lower than reported previously, suggesting that previous studies underestimated blank contributions. Second, our new CO record does not show a maximum in the late 1800s. The absence of a [CO] peak around the turn of the century argues against there being a peak in Southern Hemisphere biomass burning at this time, which is in agreement with (i) other paleofire proxies such as ethane or acetylene and (ii) conclusions reached by paleofire modeling. The combined ice core and firn air [CO] history, spanning −835 to 1992 CE, extended to the present by the Mawson atmospheric record, provides a useful benchmark for future atmospheric chemistry modeling studies.

In the 1990s, closed and open porosity volumes of firn samples were measured by J.-M. Barnola using the technique of gas pycnometry, on firn from three different polar sites. They are the basis of a parameterization of closed porosity in polar firn, first introduced in Goujon et al. (2003) and used in several firn physics models (e.g., Buizert et al., 2012). However, these data and their processing have not been published in their own right yet. In this short article, we detail how they were processed by J.-M. Barnola and how the closed porosity parameterization was obtained. We show that the original data processing only partially accounts for the presence of reopened bubbles in the samples. Since the proper correction to apply for this effect is hard to estimate, we also processed the data without including a correction for reopened bubbles. Finally, we made these pycnometry data available in order to be used by the glaciology community, notably for the study of polar ice formation and of the composition of gas records in ice cores. They are hosted on the PANGAEA database: https://doi.org/10.1594/PANGAEA.907678 (Fourteau et al., 2019a).

We measured the methane mixing ratios of enclosed air in five ice core sections drilled on the East Antarctic Plateau. Our work aims to study two effects that alter the recorded gas concentrations in ice cores: layered gas trapping artifacts and firn smoothing. Layered gas trapping artifacts are due to the heterogeneous nature of polar firn, where some strata might close early and trap abnormally old gases that appear as spurious values during measurements. The smoothing is due to the combined effects of diffusive mixing in the firn and the progressive closure of bubbles at the bottom of the firn. Consequently, the gases trapped in a given ice layer span a distribution of ages. This means that the gas concentration in an ice layer is the average value over a certain period of time, which removes the fast variability from the record. Here, we focus on the study of East Antarctic Plateau ice cores, as these low-accumulation ice cores are particularly affected by both layering and smoothing. We use high-resolution methane data to test a simple trapping model reproducing the layered gas trapping artifacts for different accumulation conditions typical of the East Antarctic Plateau. We also use the high-resolution methane measurements to estimate the gas age distributions of the enclosed air in the five newly measured ice core sections. It appears that for accumulations below 2 cm ice equivalent yr−1 the gas records experience nearly the same degree of smoothing. We therefore propose to use a single gas age distribution to represent the firn smoothing observed in the glacial ice cores of the East Antarctic Plateau. Finally, we used the layered gas trapping model and the estimation of glacial firn smoothing to quantify their potential impacts on a hypothetical 1.5-million-year-old ice core from the East Antarctic Plateau. Our results indicate that layering artifacts are no longer individually resolved in the case of very thinned ice near the bedrock. They nonetheless contribute to slight biases of the measured signal (less than 10 ppbv and 0.5 ppmv in the case of methane using our currently established continuous CH4 analysis and carbon dioxide, respectively). However, these biases are small compared to the dampening experienced by the record due to firn smoothing.

We study a firn and ice core drilled at the new “Lock-In” site in East Antarctica, located 136 km away from Concordia station towards Dumont d'Urville. High-resolution chemical and physical measurements were performed on the core, with a particular focus on the trapping zone of the firn where air bubbles are formed. We measured the air content in the ice, closed and open porous volumes in the firn, firn density, firn liquid conductivity, major ion concentrations, and methane concentrations in the ice. The closed and open porosity volumes of firn samples were obtained using the two independent methods of pycnometry and tomography, which yield similar results. The measured increase in the closed porosity with density is used to estimate the air content trapped in the ice with the aid of a simple gas-trapping model. Results show a discrepancy, with the model trapping too much air. Experimental errors have been considered but do not explain the discrepancy between the model and the observations. The model and data can be reconciled with the introduction of a reduced compression of the closed porosity compared to the open porosity. Yet, it is not clear if this limited compression of closed pores is the actual mechanism responsible for the low amount of air in the ice. High-resolution density measurements reveal the presence of strong layering, manifesting itself as centimeter-scale variations. Despite this heterogeneous stratification, all layers, including the ones that are especially dense or less dense compared to their surroundings, display similar pore morphology and closed porosity as a function of density. This implies that all layers close in a similar way, even though some close in advance or later compared to the bulk firn. Investigation of the chemistry data suggests that in the trapping zone, the observed stratification is partly related to the presence of chemical impurities.

Understanding the slow densification process of polar firn into ice is essential in order to constrain the age difference between the ice matrix and entrapped gases. The progressive microstructure evolution of the firn column with depth leads to pore closure and gas entrapment. Air transport models in the firn usually include a closed porosity profile based on available data. Pycnometry or melting–refreezing techniques have been used to obtain the ratio of closed to total porosity and air content in closed pores, respectively. X-ray-computed tomography is complementary to these methods, as it enables one to obtain the full pore network in 3-D. This study takes advantage of this nondestructive technique to discuss the morphological evolution of pores on four different Antarctic sites. The computation of refined geometrical parameters for the very cold polar sites Dome C and Lock In (the two Antarctic plateau sites studied here) provides new information that could be used in further studies. The comparison of these two sites shows a more tortuous pore network at Lock In than at Dome C, which should result in older gas ages in deep firn at Lock In. A comprehensive estimation of the different errors related to X-ray tomography and to the sample variability has been performed. The procedure described here may be used as a guideline for further experimental characterization of firn samples. We show that the closed-to-total porosity ratio, which is classically used for the detection of pore closure, is strongly affected by the sample size, the image reconstruction, and spatial heterogeneities. In this work, we introduce an alternative parameter, the connectivity index, which is practically independent of sample size and image acquisition conditions, and that accurately predicts the close-off depth and density. Its strength also lies in its simple computation, without any assumption of the pore status (open or close). The close-off prediction is obtained for Dome C and Lock In, without any further numerical simulations on images (e.g., by permeability or diffusivity calculations).

We investigate for the first time the loss and alteration of past atmospheric information from air trapping mechanisms under low-accumulation conditions through continuous CH4 (and CO) measurements. Methane concentration changes were measured over the Dansgaard–Oeschger event 17 (DO-17, ∼ 60 000 yr BP) in the Antarctic Vostok 4G-2 ice core. Measurements were performed using continuous-flow analysis combined with laser spectroscopy. The results highlight many anomalous layers at the centimeter scale that are unevenly distributed along the ice core. The anomalous methane mixing ratios differ from those in the immediate surrounding layers by up to 50 ppbv. This phenomenon can be theoretically reproduced by a simple layered trapping model, creating very localized gas age scale inversions. We propose a method for cleaning the record of anomalous values that aims at minimizing the bias in the overall signal. Once the layered-trapping-induced anomalies are removed from the record, DO-17 appears to be smoother than its equivalent record from the high-accumulation WAIS Divide ice core. This is expected due to the slower sinking and densification speeds of firn layers at lower accumulation. However, the degree of smoothing appears surprisingly similar between modern and DO-17 conditions at Vostok. This suggests that glacial records of trace gases from low-accumulation sites in the East Antarctic plateau can provide a better time resolution of past atmospheric composition changes than previously expected. We also developed a numerical method to extract the gas age distributions in ice layers after the removal of the anomalous layers based on comparison with a weakly smoothed record. It is particularly adapted for the conditions of the East Antarctic plateau, as it helps to characterize smoothing for a large range of very low-temperature and low-accumulation conditions.

The transformation of snow into ice is a complex phenomenon that is difficult to model. Depending on surface temperature and accumulation rate, it may take several decades to millennia for air to be entrapped in ice. The air is thus always younger than the surrounding ice. The resulting gas–ice age difference is essential to documenting the phasing between CO2 and temperature changes, especially during deglaciations. The air trapping depth can be inferred in the past using a firn densification model, or using δ15N of air measured in ice cores. All firn densification models applied to deglaciations show a large disagreement with δ15N measurements at several sites in East Antarctica, predicting larger firn thickness during the Last Glacial Maximum, whereas δ15N suggests a reduced firn thickness compared to the Holocene. Here we present modifications of the LGGE firn densification model, which significantly reduce the model–data mismatch for the gas trapping depth evolution over the last deglaciation at the coldest sites in East Antarctica (Vostok, Dome C), while preserving the good agreement between measured and modelled modern firn density profiles. In particular, we introduce a dependency of the creep factor on temperature and impurities in the firn densification rate calculation. The temperature influence intends to reflect the dominance of different mechanisms for firn compaction at different temperatures. We show that both the new temperature parameterization and the influence of impurities contribute to the increased agreement between modelled and measured δ15N evolution during the last deglaciation at sites with low temperature and low accumulation rate, such as Dome C or Vostok. We find that a very low sensitivity of the densification rate to temperature has to be used in the coldest conditions. The inclusion of impurity effects improves the agreement between modelled and measured δ15N at cold East Antarctic sites during the last deglaciation, but deteriorates the agreement between modelled and measured δ15N evolution at Greenland and Antarctic sites with high accumulation unless threshold effects are taken into account. We thus do not provide a definite solution to the firnification at very cold Antarctic sites but propose potential pathways for future studies.

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.

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.

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.

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.

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.