From precipitation to ice core: On the importance of surface processes for stable water-isotope records in East Antarctica

Stable water-isotope records from Antarctic ice cores allow the reconstruction of past temperature variability. However, accurate interpretation of the isotopic signal requires comprehensive understanding of the processes leading to its archiving in snow and ice, which can be documented by in situ measurements

On the similarity and apparent cycles of isotopic variations in East Antarctic snow and ice cores

Oxygen and hydrogen isotope ratios in polar ice typically show variations over a large range of timescales. Since the isotope ratios are interpreted as a proxy for atmospheric temperatures, their variations can provide essential information about the natural climate variability and cycles. Nowadays high-resolution isotope samplings corresponding to depth intervals below or around the local accumulation of snow per year are routinely performed, and observed variations in the isotopic composition at a given site have frequently been interpreted as the reflection of the seasonal cycle in temperature and also to indicate multi-year quasi-periodic climatic cycles. However, studies from strongly different accumulation conditions in East Antarctica reported similar isotopic variability and comparable apparent cycles in isotope profiles with typical wavelengths of around 20 cm, which is inconsistent with a climatically driven origin. Here we show, based on spectral analysis, that these features do not correspond to truly or quasi-periodic cycles. In addition, the typical wavelengths increase with depth for most East Antarctic sites, which is inconsistent with the effect of burial and compression on a climatic cyclic signal. We explain these results by isotopic diffusion acting on a noise-dominated isotope signal. The firn diffusion length is rather stable across the Antarctic Plateau, leading to similar power spectral densities of the isotopic variations, and increases with depth in the near-surface firn. Since the first moments of the spectral density govern the characteristic spacing of the extrema of a time series – a fundamental relationship known as Rice’s law – the similar isotope spectra in turn imply similar average distances between the isotopic minima and maxima that get larger with increasing depth. Our results bear important implications for the interpretation of isotope records in terms of cyclical climate variability. They underline that simply counting isotopic extrema is not sufficient to detect periodicities, instead robust spectral analyses have to be applied in order to differentiate between true climate cycles and the apparent cycles created in the diffusion process. This has consequences for the dating of ice-core records, which is often based on or underpinned by counting isotopic maxima, but also for the detection and interpretation of quasi-periodic climate phenomena on longer timescales. Finally, the general implications of our findings are not restricted to ice cores but likely also apply to other paleo-climate archives, as other smoothing processes, e.g. the bioturbational smoothing of proxy records from marine sediments, might lead to similar apparent cycles

Spatial variability in isotopic composition of surface snow along the East Antarctic International Ice Sheet Traverse (EAIIST)

The isotopic signal of oxygen and hydrogen, archived in the Antarctic ice sheet through snow precipitation, is an important proxy of climatic conditions. This signal depends on several parameters such as local temperature, climatic conditions in the moisture source areas and air mass pathways. Moreover, the isotopic composition may be affected by spatial variability induced by the interactions of the snow surface with the overlying atmosphere along the direction of the prevailing winds. In regions where the snow accumulation is very low, interactions between the atmosphere and the snow surface could modify the pristine signal through isotopic exchanges, sublimation processes and mechanical mixing originated from wind action. The EAIIST (East Antarctic International Ice Sheet Traverse) traverse, that took place during the 2019-2020 Antarctic field season, starting from Concordia Station towards the South Pole, provides a perfect path of study. Along the EAIIST traverse, areas with homogeneous accumulation rates can be compared to areas influenced by wind scouring and mega-dunes formation. Extremely low accumulation and wind-surface snow interaction observed in these areas could be representative of glacial period conditions in the Antarctic Plateau. Here we present the isotopic composition (dD and d18O) of surface (a few cm of depth), bulk (1 m depth) and snow pit (2 m depth) samples along the EAIIST traverse to evaluate the parameters explaining the spatial variability of this proxy. The dD, d18O and, the deuterium excess will be evaluated with respect to geographical features (elevation, latitude, distance from the coast, slope) and climatic conditions (temperature, accumulation, wind speed and direction). Wind action is expected to play a major role in explaining the isotopic spatial variability in these areas. Understanding the spatial variability in the deposition process, which strongly decreases the ratio between signal and noise, is essential to better interpret high-resolution isotopic profiles from firn and ice cores, collected along the EAIIST traverse, which will be analyzed soon

The spatial variability in isotopic composition of surface snow and snowpits on the East Antarctic Ice Sheet

The water isotope composition of snow precipitations, archived in the Antarctic ice sheet every year, is an important proxy of climatic conditions. This signal depends on several parameters such as local temperature, altitude, moisture source areas and air mass pathways. However, especially in areas where snow accumulation is very low (as on the East Antarctic Plateau), the isotopic composition is affected by additional spatial variability induced by the interactions between the atmosphere and snow surface, and the pristine signal may be modified through isotopic exchanges, sublimation processes and mechanical mixing originated from wind action. Here, we present the isotopic composition (D and 18O) and the second-order parameter d-excess of surface snow and snowpit samples collected during the Italian-French campaign in Antarctica (2019-2020). The sampling sites cover the area from Dumont D'Urville to Concordia Station and from Concordia Station towards the South Pole (EAIIST – East Antarctic International Ice Sheet Traverse). These data, compared with a previous dataset of Antarctic surface snow isotopic composition (Masson-Delmotte et al. 2008), are analyzed to determine the variability of the spatial relationship between precipitation isotopic composition and local temperature in relation to geographical parameters (latitude, distance from the coast and elevation). The interpretation of these factors determining the isotope signature is the base to better define the amount of the effects caused by subsequent interaction between atmosphere and surface snow, and by the wind action. Understanding the spatial variability of this proxy, which strongly decreases the signal-to-noise ratio, could permit to improve the use of the “isotopic thermometer” to quantify past changes in temperature based on the stable isotopic record of deep ice cores

Chronology reconstruction for the disturbed bottom section of the GISP2 and the GRIP ice cores: Implications for Termination II in Greenland

We have reconstructed chronology for the disturbed bottom parts of the GRIP and GISP2 ice cores using the combined paleoatmospheric records of CH₄ concentration and δ¹⁸O[subscript atm] in the trapped gases. Our reconstructed ages for basal ice samples are based on comparison of published measurements of CH₄ and δ¹⁸O[subscript atm] from the disturbed section of the GRIP and GISP2 cores with the same properties in the Vostok ice core. NGRIP δ¹⁸O[subscript ice] values are also used to constrain the chronology during the end of marine isotope stage 5e. For each sample, we assign an age that represents the unique or most probable time of gas trapping, given its gas composition. Of 157 samples with CH₄ and δ¹⁸O[subscript atm] data, 10 give unique ages. Twenty-five newly measured values of the triple isotope composition of O₂ from the disturbed section of the GISP2 core add a third time-dependent gas property that agrees with our reconstruction. Our reconstruction supports earlier conclusions of Landais et al. (2003) that the disturbed section primarily includes ice from the last interglacial (MIS 5e) and the penultimate glacial period (MIS 6). The oldest ice in the basal layer of GISP2 and GRIP has an age ≥237 ka. The climate history we derive suggests that the last interglacial at Summit, Greenland, around 127 ka was slightly warmer than the current interglacial period. Reduction of various ion concentrations in ice and thickening of the ice sheet during Termination II was similar to that in Termination I.Keywords: chronology, Termination II, ice cor

How warm was Greenland during the last interglacial period?

The last interglacial period (LIG, ~ 129–116 thousand years ago) provides the most recent case study for multi-millennial polar warming above pre-industrial level and a respective 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 NEEM ice core records, North-West Greenland. The NEEM paradox has emerged from an estimated large local warming above pre-industrial level (7.5 ± 1.8 °C at the deposition site 126 ka ago without correction for any overall ice sheet altitude changes between the LIG and pre-industrial) 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 +7 to +11 °C (+8 °C being the most likely estimate according to constraints on past accumulation rate) compared to the pre-industrial period. This temperature estimate is consistent with the 7.5 ± 1.8 °C warming initially determined from NEEM water isotopes. Moreover, we show that under such warm temperatures, melting of snow probably led to a significant firn shrinking by ~ 15 m. Climate simulations performed with present day ice sheet topography lead to much smaller warming but larger amplitudes (up to 5 °C) can be obtained from changes in sea ice extent and ice sheet topography. Still, ice sheet simulations forced by 5 °C surface warming lead to large ice sheet decay that are not compatible with existing data. Our new, independent temperature constrain therefore reinforces the NEEM paradox