Description: In order to interpret the paleoclimatic record stored in the air enclosed in polar ice cores, it is crucial to understand the fundamental lock-in process. Within the porous firn, bubbles are sealed continuously until the respective horizontal layer reaches a critical porosity. Present-day firn air models use a postulated temperature dependence of this value as the only parameter to adjust to the surrounding conditions of individual sites. However, no direct measurements of the firn microstructure could confirm these assumptions. Here we show that the critical porosity is a climate-independent constant by providing an extensive data set of micrometer-resolution 3-D X-ray computer tomographic measurements for ice cores representing different extremes of the temperature and accumulation ranges. We demonstrate why indirect measurements suggest a climatic dependence and substantiate our observations by applying percolation theory as a theoretical framework for bubble trapping. The incorporation of our results significantly influences the dating of trace gas records, changing gas-age–ice-age differences by up to more than 1000 years. This may further help resolve inconsistencies, such as differences between East Antarctic δ15N records (as a proxy for firn height) and model results. We expect our findings to be the basis for improved firn air and densification models, leading to lower dating uncertainties. The reduced coupling of proxies and surrounding conditions may allow for more sophisticated reinterpretations of trace gas records in terms of paleoclimatic changes and will benefit the development of new proxies, such as the air content as a marker of local insolation.

Description: Surface mass balances of polar ice sheets are essential to estimate the contribution of ice sheets to sea level rise. Uncertain snow and firn densities lead to significant uncertainties in surface mass balances, especially in the interior regions of the ice sheets, such as the East Antarctic Plateau (EAP). Robust field measurements of surface snow density are sparse and challenging due to local noise. Here, we present a snow density dataset from an overland traverse in austral summer 2016/17 on the Dronning Maud Land plateau. The sampling strategy using 1 m carbon fiber tubes covered various spatial scales, as well as a high-resolution study in a trench at 79∘ S, 30∘ E. The 1 m snow density has been derived volumetrically, and vertical snow profiles have been measured using a core-scale microfocus X-ray computer tomograph. With an error of less than 2 %, our method provides higher precision than other sampling devices of smaller volume. With four spatially independent snow profiles per location, we reduce the local noise and derive a representative 1 m snow density with an error of the mean of less than 1.5 %. Assessing sampling methods used in previous studies, we find the highest horizontal variability in density in the upper 0.3 m and therefore recommend the 1 m snow density as a robust measure of surface snow density in future studies. The average 1 m snow density across the EAP is 355 kg m−3, which we identify as representative surface snow density between Kohnen Station and Dome Fuji. We cannot detect a temporal trend caused by the temperature increase over the last 2 decades. A difference of more than 10 % to the density of 320 kg m−3 suggested by a semiempirical firn model for the same region indicates the necessity for further calibration of surface snow density parameterizations. Our data provide a solid baseline for tuning the surface snow density parameterizations for regions with low accumulation and low temperatures like the EAP.

Description: Along a traverse through North Greenland in May 2015 we collected snow cores up to 2 m depth and analyzed their density and water isotopic composition. A new sampling technique and an adapted algorithm for comparing data sets from different sites and aligning stratigraphic features are presented. We find good agreement of the density layering in the snowpack over hundreds of kilometers, which allows the construction of a representative density profile. The results are supported by an empirical statistical density model, which is used to generate sets of random profiles and validate the applied methods. Furthermore we are able to calculate annual accumulation rates, align melt layers and observe isotopic temperatures in the area back to 2010. Distinct relations of δ18O with both accumulation rate and density are deduced. Inter alia the depths of the 2012 melt layers and high-resolution densities are provided for applications in remote sensing.

Description: One of the great challenges in glaciology is the ability to estimate the bulk ice anisotropy in ice sheets and glaciers, which is needed to improve our understanding of ice-sheet dynamics. We investigate the effect of crystal anisotropy on seismic velocities in glacier ice and revisit the framework which is based on fabric eigenvalues to derive approximate seismic velocities by exploiting the assumed symmetry. In contrast to previous studies, we calculate the seismic velocities using the exact c axis angles describing the orientations of the crystal ensemble in an ice-core sample. We apply this approach to fabric data sets from an alpine and a polar ice core. Our results provide a quantitative evaluation of the earlier approximative eigenvalue framework. For near-vertical incidence our results differ by up to 135ms−1 for P-wave and 200ms−1 for S-wave velocity compared to the earlier framework (estimated 1% difference in average P-wave velocity at the bedrock for the short alpine ice core). We quantify the influence of shear-wave splitting at the bedrock as 45ms−1 for the alpine ice core and 59ms−1 for the polar ice core. At non-vertical incidence we obtain differences of up to 185ms−1 for P-wave and 280ms−1 for S-wave velocities. Additionally, our findings highlight the variation in seismic velocity at non-vertical incidence as a function of the horizontal azimuth of the seismic plane, which can be significant for non-symmetric orientation distributions and results in a strong azimuth-dependent shear-wave splitting of max. 281ms−1 at some depths. For a given incidence angle and depth we estimated changes in phase velocity of almost 200ms−1 for P wave and more than 200ms−1 for S wave and shear-wave splitting under a rotating seismic plane. We assess for the first time the change in seismic anisotropy that can be expected on a short spatial (vertical) scale in a glacier due to strong variability in crystal-orientation fabric (±50ms−1 per 10cm). Our investigation of seismic anisotropy based on ice-core data contributes to advancing the interpretation of seismic data, with respect to extracting bulk information about crystal anisotropy, without having to drill an ice core and with special regard to future applications employing ultrasonic sounding.

Description: 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 ten common IRHs in the measured RES data and the synthetic trace. Then we 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 time scale. 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.