[1] High-resolution density profiles of 16 firn cores from Greenland and Antarctica are investigated in order to improve our understanding of the densification of layered polar firn. A vertical resolution of 1-5 mm enables us to study the detailed densification processes and the evolution of the layering and the resulting variability in density with increasing depth. The densification of layered firn is important for the process of air enclosure in ice and is connected with the observed formation of a nondiffusive zone. Our findings show the following. (1) Mean density profiles, obtained from high-resolution measurements, only partly show clear transitions in densification rate at densities of 550, 730, or 820-840 kg/m 3 , as they are commonly used in literature. (2) The density variability, induced by the layering, shows a similar pattern at all sites: high variabilities at the surface, a rapid drop to a relative minimum in variability at mean density of 600-650 kg/m 3 , followed by a second relative maximum. (3) This leads to increased variability at densities of the firn-ice transition for most of the sites. (4) The variability at the surface decreases with increasing mean annual temperature and accumulation rate, whereas the variability at the firn-ice transition increases. We can exclude a change in local climate conditions as an explanation for the density variability since the firn cores in this study cover a broad range in mean annual temperature, accumulation rate, and age. Overall, high-resolution density profiles deliver a more complex picture of compaction of polar firn as a layered granular medium than has been obtained from mean density profiles in the past.

Editor: P. DeMenocalKeywords: polar firn density densification impurity density variability Understanding polar firn densification is crucial for reconstructing the age of greenhouse gas concentrations extracted from ice cores, and for the interpretation of air in ice as a dating tool or as a climate proxy. Firn densification is generally modeled as a steady burial and sintering process of defined layers, where the structure of the layering is maintained along the whole firn and ice column. However, available high-resolution density data, as well as firn air samples, question this picture and point to a lack of understanding of firn densification. Based on analysis of high-resolution density and calcium concentration records from Antarctic and Greenland ice cores, we show for the first time that also impurities may have a significant impact on the densification. Analysis of firn cores shows a correlation between density and the calcium ion (Ca++) concentration, and this correlation increases with depth. The existence of this relationship is independent of the local climatic conditions at the core sites analyzed. The strong positive correlation between the density and the logarithm of Ca++ concentration indicates that impurities induce softening and lead to faster densification over a wide range of concentrations. In one core, the impurity effect manifests itself so strongly that the density develops a seasonal cycle closely following the seasonal cycle of Ca++. Our results clearly show that the structure of the firn layering changes with depth and suggest that the increased variability in density observed in deep firn, recently described as a universal feature of polar firn, may arise from the influence of Ca++ and/or other impurities. The impurity effect is likely to have direct implications on our understanding of glacial firn densification and on glacial gas age estimates.

The Sentinel Application Platform (SNAP) architecture facilitates Earth Observation data processing. In this work, we present results from a new Snow Processor for SNAP. We also describe physical principles behind the developed snow property retrieval technique based on the analysis of Ocean and Land Colour Instrument (OLCI) onboard Sentinel-3A/B measurements over clean and polluted snow fields. Using OLCI spectral reflectance measurements in the range 400–1020 nm, we derived important snow properties such as spectral and broadband albedo, snow specific surface area, snow extent and grain size on a spatial grid of 300 m. The algorithm also incorporated cloud screening and atmospheric correction procedures over snow surfaces. We present validation results using ground measurements from Antarctica, the Greenland ice sheet and the French Alps. We find the spectral albedo retrieved with accuracy of better than 3% on average, making our retrievals sufficient for a variety of applications. Broadband albedo is retrieved with the average accuracy of about 5% over snow. Therefore, the uncertainties of satellite retrievals are close to experimental errors of ground measurements. The retrieved surface grain size shows good agreement with ground observations. Snow specific surface area observations are also consistent with our OLCI retrievals. We present snow albedo and grain size mapping over the inland ice sheet of Greenland for areas including dry snow, melted/melting snow and impurity rich bare ice. The algorithm can be applied to OLCI Sentinel-3 measurements providing an opportunity for creation of long-term snow property records essential for climate monitoring and data assimilation studies—especially in the Arctic region, where we face rapid environmental changes including reduction of snow/ice extent and, therefore, planetary albedo.

No abstract

Several recent studies from both Greenland and Antarctica have reported significant changes in the water isotopic composition of near‐surface snow between precipitation events. These changes have been linked to isotopic exchange with atmospheric water vapor and sublimation‐induced fractionation, but the processes are poorly constrained by observations. Understanding and quantifying these processes are crucial to both the interpretation of ice core climate proxies and the formulation of isotope‐enabled general circulation models. Here, we present continuous measurements of the water isotopic composition in surface snow and atmospheric vapor together with near‐surface atmospheric turbulence and snow‐air latent and sensible heat fluxes, obtained at the East Greenland Ice‐Core Project drilling site in summer 2016. For two 4‐day‐long time periods, significant diurnal variations in atmospheric water isotopologues are observed. A model is developed to explore the impact of this variability on the surface snow isotopic composition. Our model suggests that the snow isotopic composition in the upper subcentimeter of the snow exhibits a diurnal variation with amplitudes in δ 18 O and δD of ~2.5‰ and ~13‰, respectively. As comparison, such changes correspond to 10–20% of the magnitude of seasonal changes in interior Greenland snow pack isotopes and of the change across a glacial‐interglacial transition. Importantly, our observation and model results suggest, that sublimation‐induced fractionation needs to be included in simulations of exchanges between the vapor and the snow surface on diurnal timescales during summer cloud‐free conditions in northeast Greenland.

Abstract. Stable isotope ratios δ 18 O and δD in polar ice provide a wealth of information about past climate evolution. Snow-pit studies allow us to relate observed weather and climate conditions to the measured isotope variations in the snow. They therefore offer the possibility to test our understanding of how isotope signals are formed and stored in firn and ice. As δ 18 O and δD in the snowfall are strongly correlated to air temperature, isotopes in the near-surface snow are thought to record the seasonal cycle at a given site. Accordingly, the number of seasonal cycles observed over a given depth should depend on the accumulation rate of snow. However, snow-pit studies from different accumulation conditions in East Antarctica reported similar isotopic variability and comparable apparent cycles in the δ 18 O and δD profiles with typical wavelengths of ∼ 20 cm. These observations are unexpected as the accumulation rates strongly differ between the sites, ranging from 20 to 80 mm w.e. yr −1 (∼ 6-21 cm of snow per year). Various mechanisms have been proposed to explain the isotopic variations individually at each site; however, none of these are consistent with the similarity of the different profiles independent of the local accumulation conditions.Here, we systematically analyse the properties and origins of δ 18 O and δD variations in high-resolution firn profiles from eight East Antarctic sites. First, we confirm the suggested cycle length (mean distance between peaks) of ∼ 20 cm by counting the isotopic maxima. Spectral analysis further shows a strong similarity between the sites but indicates no dominant periodic features. Furthermore, the apparent cycle length increases with depth for most East Antarctic sites, which is inconsistent with burial and compression of a regular seasonal cycle. We show that these results can be explained by isotopic diffusion acting on a noise-dominated isotope signal. The firn diffusion length is rather stable across the Antarctic Plateau and thus leads to similar power spectral densities of the isotopic variations. This in turn implies a similar distance between isotopic maxima in the firn profiles.Our results explain a large set of observations discussed in the literature, providing a simple explanation for the interpretation of apparent cycles in shallow isotope records, without invoking complex mechanisms. Finally, the results underline previous suggestions that isotope signals in single ice cores from low-accumulation regions have a small signal-to-noise ratio and thus likely do not allow the reconstruction of interannual to decadal climate variations.

The density of firn is an important property for monitoring and modeling the ice sheets as well as to model the pore close‐off and thus to interpret ice core‐based greenhouse gas records. One feature, which is still in debate, is the potential existence of an annual cycle of firn density in low‐accumulation regions. Several studies describe or assume seasonally successive density layers, horizontally evenly distributed, as seen in radar data. On the other hand, high‐resolution density measurements on firn cores in Antarctica and Greenland show no clear seasonal cycle in the top few meters. A major caveat of most existing snow‐pit and firn‐core‐based studies is that they represent one vertical profile from a laterally heterogeneous density field. To overcome this, we created an extensive data set of horizontal and vertical density data at Kohnen Station, Dronning Maud Land, on the East Antarctic Plateau. We drilled and analyzed three 90 m long firn cores as well as 143 one‐meter‐long vertical profiles from two elongated snow trenches to obtain a two‐dimensional view of the density variations. The analysis of the 45 m wide and 1 m deep density fields reveals a seasonal cycle in density. However, the seasonality is overprinted by strong stratigraphic noise, making it invisible when analyzing single firn cores. Our density data set extends the view from the local ice core perspective to a hundred meter scale and thus supports linking spatially integrating methods such as radar and seismic studies to ice and firn cores.

[1] Recent advances in three-dimensional (3D) imaging of snow and firn combined with numerical modeling of flow through complex geometries have greatly improved the ability to predict permeability values based on microstructure. In this work, we combined 3D reconstructions of polar firn microstructure obtained from microcomputed tomography (mCT) and a 3D lattice-Boltzmann (LB) model of air flow. We compared the modeled results to measurements of permeability for polar firn with a wide range of grain and pore-scale characteristics. The results show good agreement between permeability measurements and calculated permeability values from the LB model over a wide range of sample types. The LB model is better at predicting measured permeability values than traditional empirical equations for polar firn.