How well can satellite altimetry and firn models resolve Antarctic firn thickness variations?

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.