Explore the vast range of titles available.

Refreezing of meltwater in firn is a major component of Greenland ice-sheet's mass budget, but in situ observations are rare. Here, we compare the firn density and total ice layer thickness in the upper 15 m of 19 new and 27 previously published firn cores drilled at 15 locations in southwest Greenland (1850–2360 m a.s.l.) between 1989 and 2019. At all sites, ice layer thickness covaries with density over time and space. At the two sites with the earliest observations (1989 and 1998), bulk density increased by 15–18%, in the top 15 m over 28 and 21 years, respectively. However, following the extreme melt in 2012, elevation-detrended density using 30 cores from all sites decreased by 15 kg m−3 a−1 in the top 3.75 m between 2013 and 2019. In contrast, the lowest elevation site's density shows no trend. Thus, temporary build-up in firn pore space and meltwater infiltration capacity is possible despite the long-term increase in Greenland ice-sheet melting.

There is ample evidence for ice layers and lenses within glacial firn. The standard model for ice layer formation localizes the refreezing by perching of meltwater on pre‐existing discontinuities. Here we argue that even extreme melting events provide insufficient flux for this mechanism. Using a thermomechanical model we demonstrate a different mechanism of ice layer formation. After a melting event when the drying front catches up with the wetting front and arrests melt percolation, conductive heat loss freezes the remaining melt in place to form an ice layer. This model reproduces the depth of a new ice layer at the Dye‐2 site in Greenland. It provides a deeper insight into the interpretation of firn stratigraphy and past climate variability. It also improves the simulation of firn densification processes, a key source of uncertainty in assessing and attributing ice sheet mass balance based on satellite altimetry and gravimetry data. Plain Language Summary Firn covers a significant portion of Earth's glaciers and ice sheets. It can store surface meltwater and prevent runoff into the ocean. The widespread presence of ice layers embedded in firn formed by meltwater refreezing may prevent meltwater storage and contribute to sea level rise. However, current models of ice layer formation, originally developed for snow, do not seem to work in firn. This work presents a different mechanism for ice layer formation without invoking pre‐existing ice layers within the firn. Our model shows that the sequencing of ice layers formed by subsequent melting events depends on the overall heat added to the firn. Deeper layers occur in warmer, more porous firn during intense melt events in a warming climate. This insight enhances our understanding of firn layering and can help deduce past climate variations. Our model aids in understanding the density evolution of firn to reduce uncertainties in remote sensing data that determines the ice sheet mass loss and its contribution to global sea‐level rise. Key Points A new ice layer forms without meltwater perching when freezing localizes after a melt event as heat conduction dominates over advection Deeper ice layers form in warming climatic conditions whereas denser ice layers form near surface in net‐zero climatic conditions Results indicate the possibility of deducing past variability in climate from firn stratigraphy and vice versa

Firn simulations are essential for understanding Antarctic ice sheet mass change, as they enable us to convert satellite altimetry observed volume changes to mass changes and column thickness to ice thickness and to quantify the meltwater buffering capacity of firn. Here, we present and evaluate a simulation of the contemporary Antarctic firn layer using the updated semi-empirical IMAU Firn Densification Model (IMAU-FDM) for the period 1979–2020. We have improved previous fresh-snow density and firn compaction parameterizations and used updated atmospheric forcing. In addition, the model has been calibrated and evaluated using 112 firn core density observations across the ice sheet. We found that 62 % of the seasonal and 67 % of the decadal surface height variability are due to variations in firn air content rather than firn mass. Comparison of simulated surface elevation change with a previously published multi-mission altimetry product for the period 2003–2015 shows that performance of the updated model has improved, notably in Dronning Maud Land and Wilkes Land. However, a substantial trend difference (>10 cm yr−1) remains in the Antarctic Peninsula and Ellsworth Land, mainly caused by uncertainties in the spin-up forcing. By estimating previous climatic conditions from ice core data, these trend differences can be reduced by 38 %.

Evolution of cold dry snow and firn plays important roles in glaciology; however, the physical formulation of a densification law is still an active research topic. We forced eight firn-densification models and one seasonal-snow model in six different experiments by imposing step changes in temperature and accumulation-rate boundary conditions; all of the boundary conditions were chosen to simulate firn densification in cold, dry environments. While the intended application of the participating models varies, they are describing the same physical system and should in principle yield the same solutions. The firn models all produce plausible depth-density profiles, but the model outputs in both steady state and transient modes differ for quantities that are of interest in ice core and altimetry research. These differences demonstrate that firn-densification models are incorrectly or incompletely representing physical processes. We quantitatively characterize the differences among the results from the various models. For example, we find depth-integrated porosity is unlikely to be inferred with confidence from a firn model to better than 2 m in steady state at a specific site with known accumulation rate and temperature. Firn Model Intercomparison Experiment can provide a benchmark of results for future models, provide a basis to quantify model uncertainties and guide future directions of firn-densification modeling.

The density structure of firn has implications for hydrological and climate modelling, and ice-shelf stability. The structure of firn can be evaluated from depth models of seismic velocity, widely obtained with Herglotz–Wiechert inversion (HWI), an approach that considers the slowness of refracted seismic arrivals. However, HWI is strictly appropriate only for steady-state firn profiles and the inversion accuracy can be compromised where firn contains ice layers. In these cases, full waveform inversion (FWI) may yield more success than HWI. FWI extends HWI capabilities by considering the full seismic waveform and incorporates reflected arrivals. Using synthetic firn density profiles, assuming both steady- and non-steady-state accumulation, we show that FWI outperforms HWI for detecting ice slab boundaries (5–80 m thick, 5–80 m deep) and velocity anomalies within firn. FWI can detect slabs thicker than one wavelength (here, 20 m, assuming a maximum frequency of 60 Hz) but requires the starting velocity model to be accurate to ±2.5%. We recommend for field practice that the shallowest layers of velocity models are constrained with ground-truth data. Nonetheless, FWI shows advantages over established methods, and should be considered when the characterisation of firn ice slabs is the goal of the seismic survey.

Processes governing meltwater penetration into cold firn remain poorly constrained. Here, in situ experiments are used to develop a grain-scale model to investigate physical limitations on meltwater infiltration in firn. At two sites in Greenland, drilling pumped water into cold firn to >75 m depth, and the thermo-hydrologic evolution of the firn column was measured. Rather than filling all available pore space, the water formed perched aquifers with downward penetration halted by thermal and density conditions. The two sites formed deep aquifers at ~40 m depth and at densities considerably less than the air pore close-off density (~725 kg m−3 at −18°C, and ~750 kg m−3 at −14°C), demonstrating that some pore space at depth remains inaccessible. A geometric grain-scale model of firn is constructed to quantify the limits of a descending fully saturated wetting front in cold firn. Agreement between the model and field data implies the model includes the first-order effects of water and heat flow in a firn lattice. The model constrains the relative importance of firn density, temperature and grain/pore size in inhibiting wetting front migration. Results imply that deep infiltration, including that which leads to firn aquifer formation, does not have access to all available firn pore space.

Constraining the onset of dynamic recrystallization (DRX) and its effects on the mechanical properties of firn is crucial for firn densification modeling. To that end, samples from a depth of 13 m in a Summit, Greenland (72°35′N, 38°25′W) firn core were subjected to creep tests at −14°C and 0.21 MPa compressive stress to strains of 7%, 12%, 18%, and 29%. Microstructural analyses using thin‐section imaging and microcomputed x‐ray tomography (micro‐CT) revealed smaller grain sizes, reduced specific surface area and connectivity, and increased density in relation to reduced porosity as the strain increases. These results show that DRX occurs in firn under creep, with strain‐induced boundary migration (SIBM) and nucleation and growth starting at ∼7%. DRX leads to elongated grains, reduced grain size, and the development of a preferred crystallographic orientation, indicating that DRX occurs by both SIBM and nucleation and growth. Plain Language Summary Firn is multi‐year snow that undergoes densification due to the load from the snow overburden and from sintering. Understanding firn densification is important for interpreting ice core records, predicting ice sheet mass balance and elevation changes, and studying climate change effects. Previous densification models focused on accumulation rate and temperature, overlooking the role of recrystallization. To address this gap, compression tests were performed on Greenland firn samples from a depth of 13 m. The deformation resulted in reduced median grain size, preferred crystallographic orientation, and increased density. Our findings indicate that dynamic recrystallization starts when the firn is subjected to a strain of about 7% through boundary migration of old grains, around which new stress‐free grains also start to form. Key Points Dynamic Recrystallization occurs in firn through strain‐induced boundary migration, and nucleation and growth Average grain size in firn decreases under constant temperature and compressive stress in firn

Understanding firn densification is essential for interpreting ice core records, predicting ice sheet mass balance, elevation changes and future sea-level rise. Current models of firn densification on the Antarctic ice sheet (AIS), such as the Herron and Langway (1980) model are either simple semi-empirical models that rely on sparse climatic data and surface density observations or complex physics-based models that rely on poorly understood physics. In this work, we introduce a deep learning technique to study firn densification on the AIS. Our model, FirnLearn, evaluated on 225 cores, shows an average root-mean-square error of 31 kg m−3 and explained variance of 91%. We use the model to generate surface density and the depths to the $550\\,\\mathrm{kg\\,m}^{-3}$ and $830\\,\\mathrm{kg\\,m}^{-3}$ density horizons across the AIS to assess spatial variability. Comparisons with the Herron and Langway (1980) model at ten locations with different climate conditions demonstrate that FirnLearn more accurately predicts density profiles in the second stage of densification and complete density profiles without direct surface density observations. This work establishes deep learning as a promising tool for understanding firn processes and advancing towards a universally applicable firn model.

Accurate modeling of firn densification is necessary for ice core interpretation and assessing the mass balance of glaciers and ice sheets. In this paper, we revisit the nonlinear-viscous firn rheology introduced by Gagliardini and Meyssonnier (1997) that allows multidimensional firn densification problems to be posed, subject to arbitrary stress and temperature fields. First, we extend the calibration of the coefficient functions that control firn compressibility and viscosity to five additional Greenlandic sites, showing that the original calibration is not universally valid. Next, we demonstrate that the transient collapse of a Greenlandic firn tunnel can be reproduced in a cross-section model, but that anomalous warm summer temperatures during 2012–14 reduce confidence in attempts to independently validate the rheology. Finally, we show that the rheology can explain the increased densification rate and varying bubble close-off depth observed across the shear margins of the Northeast Greenland Ice Stream. Although we suggest more work is needed to constrain the near-surface compressibility and viscosity functions of the rheology, our results strengthen the empirical grounding of the rheology for future use, such as modeling horizontal firn density variations over ice sheets for mass-loss estimates or estimating ice-gas age differences in ice cores subject to complex strain histories.

The firn layer covers 98% of Antarctica's ice sheets, protecting underlying glacial ice from the external environment. Accurate measurement of firn properties is essential for assessing cryosphere mass balance and climate change impacts. Characterizing firn structure through core sampling is expensive and logistically challenging. Seismic surveys, which translate seismic velocities into firn densities, offer an efficient alternative. This study employs Distributed Acoustic Sensing technology to transform an existing fiber‐optic cable near the South Pole into a multichannel, low‐maintenance, continuously interrogated seismic array. The data resolve 16 seismic wave propagation modes at frequencies up to 100 Hz that constrain P and S wave velocities as functions of depth. Using co‐located geophones for ambient noise interferometry, we resolve very weak radial anisotropy. Leveraging nearby SPICEcore firn density data, we find prior empirical density‐velocity relationships underestimate firn air content by over 15%. We present a new empirical relationship for the South Pole region. Plain Language Summary Firn, the layer of compacted snow merging into glacial ice covering Antarctica, acts as an insulating blanket that mitigates environmental perturbations to the polar ice sheet. Understanding the density and seismic characteristics of the firn layer helps scientists better infer its properties and variation, including factors relevant to glacial stability and sea level change. Firn density is the major uncertainty source for measuring ice sheet mass changes via satellite and airborne sensing. Traditional methods of assessing firn density involve drilling or snow pit analyses and are expensive and time‐consuming. We utilize the rapidly developing technology of Distributed Acoustic Sensing to transform a data communication cable near the South Pole into a dense array of seismic sensors, allowing us to noninvasively estimate firn properties by studying seismic waves propagating in the firn to assess its physical properties. Our findings suggest that previous parameterizations overestimate firn density by over 5% and underestimate its air content by over 15% and highlight the value of seismology for advancing glaciological and polar region's climate research. Key Points Distributed Acoustic Sensing repurposes an 8 km fiber‐optic cable at the South Pole into a dense seismic array Gathered data resolve 16 dispersion modes at frequencies up to 100 Hz that constrain P‐ and S‐wave velocities in the firn layer Previous density‐velocity empirical relations overestimate the dry firn density at South Pole