ALBERT

Description: The last glacial period is characterized by a number of millennial climate events that have been identified in both Greenland and Antarctic ice cores and that are abrupt in Greenland climate records. The mechanisms governing this climate variability remain a puzzle that requires a precise synchronization of ice cores from the two hemispheres to be resolved. Previously, Greenland and Antarctic ice cores have been synchronized primarily via their common records of gas concentrations or isotopes from the trapped air and via cosmogenic isotopes measured on the ice. In this work, we apply ice core volcanic proxies and annual layer counting to identify large volcanic eruptions that have left a signature in both Greenland and Antarctica. Generally, no tephra is associated with those eruptions in the ice cores, so the source of the eruptions cannot be identified. Instead, we identify and match sequences of volcanic eruptions with bipolar distribution of sulfate, i.e. unique patterns of volcanic events separated by the same number of years at the two poles. Using this approach, we pinpoint 82 large bipolar volcanic eruptions throughout the second half of the last glacial period (12–60 ka). This improved ice core synchronization is applied to determine the bipolar phasing of abrupt climate change events at decadal-scale precision. In response to Greenland abrupt climatic transitions, we find a response in the Antarctic water isotope signals (δ18O and deuterium excess) that is both more immediate and more abrupt than that found with previous gas-based interpolar synchronizations, providing additional support for our volcanic framework. On average, the Antarctic bipolar seesaw climate response lags the midpoint of Greenland abrupt δ18O transitions by 122±24 years. The time difference between Antarctic signals in deuterium excess and δ18O, which likewise informs the time needed to propagate the signal as described by the theory of the bipolar seesaw but is less sensitive to synchronization errors, suggests an Antarctic δ18O lag behind Greenland of 152±37 years. These estimates are shorter than the 200 years suggested by earlier gas-based synchronizations. As before, we find variations in the timing and duration between the response at different sites and for different events suggesting an interaction of oceanic and atmospheric teleconnection patterns as well as internal climate variability.

Description: 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.

Description: 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.

Description: We investigate for the first time through continuous measurements the loss and alteration of past atmospheric information from air trapping mechanisms under low accumulation conditions. Methane concentration changes were measured over the Dansgaard-Oeschger event 17 (D0-17, ~ 60,000 yrBP) 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, 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 which aims at minimizing the bias in the overall signal. Once the layered-trapping induced anomalies are removed from the record, the 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 and surprisingly, the degree of smoothing appears 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 usually expected. We also developed a method to extract the gas age distributions in ice layers that can be applied even for sites without firn-air measurements. 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 accumulation conditions.

Description: Using new high-resolution 10Be measurements in the NGRIP, EDML and Vostok ice cores, together with previously published data from EDC, we present an improved synchronization between Greenland and Antarctic ice cores during the Laschamp geomagnetic excursion  ∼  41 kyr ago. We estimate the precision of this synchronization to be ±20 years, an order of magnitude better than previous work. We discuss the implications of this new synchronization for making improved estimates of the depth difference between ice and enclosed gas of the same age (Δdepth), difference between age of ice and enclosed gas at the same depth (Δage) in the EDC and EDML ice cores, spectral properties of the 10Be profiles and phasing between Dansgaard–Oeschger-10 (in NGRIP) and AIM-10 (in EDML and EDC).

Description: Continental relative humidity (RH) is a key climate parameter, but there is a lack of quantitative RH proxies suitable for climate model–data comparisons. Recently, a combination of climate chamber and natural transect calibrations have laid the groundwork for examining the robustness of the triple oxygen isotope composition (δ′18O and 17O-excess) of phytoliths, that can preserve in sediments, as a new proxy for past changes in RH. However, it was recommended that besides RH, additional factors that may impact δ′18O and 17O-excess of plant water and phytoliths be examined. Here, the effects of grass leaf length, leaf development stage and day–night alternations are addressed from growth chamber experiments. The triple oxygen isotope compositions of leaf water and phytoliths of the grass species F. arundinacea are analysed. Evolution of the leaf water δ′18O and 17O-excess along the leaf length can be modelled using a string-of-lakes approach to which an unevaporated–evaporated mixing equation must be added. We show that for phytoliths to record this evolution, a kinetic fractionation between leaf water and silica, increasing from the base to the apex, must be assumed. Despite the isotope heterogeneity of leaf water along the leaf length, the bulk leaf phytolith δ′18O and 17O-excess values can be estimated from the Craig and Gordon model and a mean leaf water–phytolith fractionation exponent (λPhyto-LW) of 0.521. In addition to not being leaf length dependent, δ′18O and 17O-excess of grass phytoliths are expected to be impacted only very slightly by the stem vs. leaf biomass ratio. Our experiment additionally shows that because a lot of silica polymerises in grasses when the leaf reaches senescence (58 % of leaf phytoliths in mass), RH prevailing during the start of senescence should be considered in addition to RH prevailing during leaf growth when interpreting the 17O-excess of grass bulk phytoliths. Although under the study conditions 17O-excessPhyto do not vary significantly from constant day to day–night conditions, additional monitoring at low RH conditions should be done before drawing any generalisable conclusions. Overall, this study strengthens the reliability of the 17O-excess of phytoliths to be used as a proxy of RH. If future studies show that the mean value of 0.521 used for the grass leaf water–phytolith fractionation exponent λPhyto-LW is not climate dependent, then grassland leaf water 17O-excess obtained from grassland phytolith 17O-excess would inform on isotope signals of several soil–plant-atmosphere processes.

Description: The isotopic compositions of oxygen and hydrogen in ice cores are invaluable tools for the reconstruction of past climate variations. Used alone, they give insights into the variations of the local temperature, whereas taken together they can provide information on the climatic conditions at the point of origin of the moisture. However, recent analyses of snow from shallow pits indicate that the climatic signal can become erased in very low accumulation regions, due to local processes of snow reworking. The signal-to-noise ratio decreases and the climatic signal can then only be retrieved using stacks of several snow pits. Obviously, the signal is not completely lost at this stage, otherwise it would be impossible to extract valuable climate information from ice cores as has been done, for instance, for the last glaciation. To better understand how the climatic signal is passed from the precipitation to the snow, we present here results from varied snow samples from East Antarctica. First, we look at the relationship between isotopes and temperature from a geographical point of view, using results from three traverses across Antarctica, to see how the relationship is built up through the distillation process. We also take advantage of these measures to see how second-order parameters (d-excess and 17O-excess) are related to δ18O and how they are controlled. d-excess increases in the interior of the continent (i.e., when δ18O decreases), due to the distillation process, whereas 17O-excess decreases in remote areas, due to kinetic fractionation at low temperature. In both cases, these changes are associated with the loss of original information regarding the source. Then, we look at the same relationships in precipitation samples collected over 1 year at Dome C and Vostok, as well as in surface snow at Dome C. We note that the slope of the δ18O vs. temperature (T) relationship decreases in these samples compared to those from the traverses, and thus caution is advocated when using spatial slopes for past climate reconstruction. The second-order parameters behave in the same way in the precipitation as in the surface snow from traverses, indicating that similar processes are active and that their interpretation in terms of source climatic parameters is strongly complicated by local temperature effects in East Antarctica. Finally we check if the same relationships between δ18O and second-order parameters are also found in the snow from four snow pits. While the d-excess remains opposed to δ18O in most snow pits, the 17O-excess is no longer positively correlated to δ18O and even shows anti-correlation to δ18O at Vostok. This may be due to a stratospheric influence at this site and/or to post-deposition processes.

Description: The oldest ice core records are obtained from the East Antarctic plateau. Water stable isotopes records are key for reconstructions of past climatic conditions both over the ice sheet and at the evaporation source. The accuracy of such climate reconstructions crucially depends on the knowledge of all the processes affecting the water vapour, precipitation and snow isotopic composition. Atmospheric fractionation processes are well understood and can be integrated in Rayleigh distillation and complex isotope enabled climate models. However, a comprehensive quantitative understanding of processes potentially altering the snow isotopic composition after the deposition is still missing, especially for exchanges between vapour and snow. In low accumulation sites such as found on the East Antarctic Plateau, these poorly constrained processes are especially likely to play a significant role. This limits the interpretation of isotopic composition from ice core records, specifically at short time scales. Here, we combine observations of isotopic composition in the vapour, the precipitation, the surface snow and the buried snow from various sites of the East Antarctic Plateau. At the seasonal scale, we highlight a significant impact of metamorphism on surface snow isotopic signal compared to the initial precipitation isotopic signal. In particular, in summer, exchanges of water molecules between vapour and snow are driven by the sublimation/condensation cycles at the diurnal scale. Using highly resolved isotopic composition profiles from pits in five East Antarctic sites, we identify a common 20 cm cycle which cannot be attributed to the seasonal variability of precipitation. Altogether, the smaller range of isotopic compositions observed in the buried and in the surface snow compared to the precipitation, and also the reduced slope between surface snow isotopic composition and temperature compared to precipitation, constitute evidences of post-deposition processes affecting the variability of the isotopic composition in the snow pack. To reproduce these processes in snow-models is crucial to understand the link between snow isotopic composition and climatic conditions and to improve the interpretation of isotopic composition as a paleoclimate proxy.

Description: The last interglacial period (LIG, ∼ 129–116 thousand years ago) provides the most recent case study of multimillennial polar warming above the preindustrial level and a response of the Greenland and Antarctic ice sheets to this warming, as well as a test bed for climate and ice sheet models. Past changes in Greenland ice sheet thickness and surface temperature during this period were recently derived from the North Greenland Eemian Ice Drilling (NEEM) ice core records, northwest Greenland. The NEEM paradox has emerged from an estimated large local warming above the preindustrial level (7.5 ± 1.8 °C at the deposition site 126 kyr ago without correction for any overall ice sheet altitude changes between the LIG and the preindustrial period) based on water isotopes, together with limited local ice thinning, suggesting more resilience of the real Greenland ice sheet than shown in some ice sheet models. Here, we provide an independent assessment of the average LIG Greenland surface warming using ice core air isotopic composition (δ15N) and relationships between accumulation rate and temperature. The LIG surface temperature at the upstream NEEM deposition site without ice sheet altitude correction is estimated to be warmer by +8.5 ± 2.5 °C compared to the preindustrial period. This temperature estimate is consistent with the 7.5 ± 1.8 °C warming initially determined from NEEM water isotopes but at the upper end of the preindustrial period to LIG temperature difference of +5.2 ± 2.3 °C obtained at the NGRIP (North Greenland Ice Core Project) site by the same method. Climate simulations performed with present-day ice sheet topography lead in general to a warming smaller than reconstructed, but sensitivity tests show that larger amplitudes (up to 5 °C) are produced in response to prescribed changes in sea ice extent and ice sheet topography.

Description: The transformation of snow into ice is a complex phenomenon 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 that the surrounding ice. The resulting gas-ice age difference is essential to document 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 in several sites of East Antarctica, predicting larger firn thickness during the Last Glacial Maximum, whereas δ15N suggests a reduced firn thickness compared to the Holocene. We present here modifications of the LGGE firn densification model, which significantly reduce the model-data mismatch for the gas trapping depth evolution over the last deglaciation, while preserving the good agreement between measured and modelled modern firn density profiles. In particular, we introduce a dependency of the activation energy to temperature and impurities in the firn densification rate calculation. The temperature influence reflects the existence 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. However, the inclusion of impurities effects deteriorates the agreement between modelled and measured δ15N evolution in Greenland and Antarctic sites with high accumulation.