Search

Elevation changes of the Antarctic Ice Sheet (AIS) related to surface mass balance and firn processes vary strongly in space and time. Their subdecadal natural variability is large and hampers the detection of long-term climate trends. Firn models or satellite altimetry observations are typically used to investigate such firn thickness changes. However, there is a large spread among firn models. Further, they do not fully explain observed firn thickness changes, especially on smaller spatial scales. Reconciled firn thickness variations will facilitate the detection of long-term trends from satellite altimetry; the resolution of the spatial patterns of such trends; and, hence, their attribution to the underlying mechanisms. This study has two objectives. First, we quantify interannual Antarctic firn thickness variations on a 10 km grid scale. Second, we characterise errors in both the altimetry products and firn models. To achieve this, we jointly analyse satellite altimetry and firn modelling results in time and space. We use the timing of firn thickness variations from firn models and the satellite-observed amplitude of these variations to generate a combined product (“adjusted firn thickness variations”) over the AIS for 1992–2017. The combined product characterises spatially resolved variations better than either firn models alone or altimetry alone. It provides a higher resolution and a more precise spatial distribution of the variations compared to model-only solutions and eliminates most of the altimetry errors compared to altimetry-only solutions. Relative uncertainties in basin-mean time series of the adjusted firn thickness variations range from 20 % to 108 %. At the grid cell level, relative uncertainties are higher, with median values per basin in the range of 54 % to 186 %. This is due to the uncertainties in the large and very dry areas of central East Antarctica, especially over large megadune fields, where the low signal-to-noise ratio poses a challenge for both models and altimetry to resolve firn thickness variations. A large part of the variance in the altimetric time series is not explained by the adjusted firn thickness variations. Analysis of the altimetric residuals indicate that they contain firn model errors, such as firn signals not captured by the models, and altimetry errors, such as time-variable radar penetration effects and errors in intermission calibration. This highlights the need for improvements in firn modelling and altimetry analysis.

The Antarctic firn layer provides pore space in which an estimated 94 % to 96 % of the surface melt refreezes or is retained as liquid water. Future depletion of firn pore space by increased surface melt, densification and formation of low-permeability ice slabs can potentially lead to meltwater ponding, hydrofracturing and ice-shelf disintegration. Here, we investigate the 21st-century evolution of total firn air content (FAC) and accessible FAC (i.e. the pore space that meltwater can reach) across Antarctic ice shelves. We use the semi-empirical IMAU Firn Densification Model (IMAU-FDM) with an updated dynamical densification expression to cope with changing climate forcing. The firn model is forced by general circulation model output of the Community Earth System Model version 2 (CESM2) for three climate emission scenarios (SSP1-2.6, SSP2-4.5 and SSP5-8.5), dynamically downscaled to a 27 km horizontal resolution by the Regional Atmospheric Climate Model version 2.3p2 (RACMO2.3p2). To estimate accessible FAC, we prescribe a relationship between ice-slab thickness and permeability. In our simulations, ice shelves on the Antarctic Peninsula and the Roi Baudouin Ice Shelf in Dronning Maud Land are particularly vulnerable to total FAC depletion (> 50 % decrease by 2100), even for low-emission (SSP1-2.6) and intermediate-emission (SSP2-4.5) scenarios. In the high-emission (SSP5-8.5) scenario in particular, the formation of ice slabs further reduces accessible FAC on ice shelves with low accumulation rates (current rates of < 500 mmw.e.yr-1), including many East Antarctic ice shelves and the Filchner–Ronne, Ross, Pine Island and Larsen C ice shelves. These results underline the potentially large vulnerability of low-accumulation ice shelves to firn air depletion through ice-slab formation.

Abstract Despite in‐situ observations of perennial firn aquifers (PFAs) at specific locations of the Antarctic ice sheet, a comprehensive continent‐wide mapping of PFA distribution is currently lacking. We present an estimate of their distribution across Antarctica in the form of a probability assessment using a Monte Carlo technique. Our approach involves a novel methodology that combines observations from Sentinel‐1 and Advanced SCATterometer (ASCAT) with output from a regional climate model. To evaluate our method, we conduct an extensive comparison with Operation Ice Bridge observations from the Greenland Ice Sheet. Application to Antarctica reveals high PFA probabilities in the Antarctic Peninsula (AP), particularly along its northern, northwestern, and western coastlines, as well as on the Wilkins, Müller, and George VI ice shelves. Outside the AP, PFA probability is low, except for some locations with marginally higher probabilities, such as on the Abbot, Totten, and Shackleton ice shelves.

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 %.