Abstract. Estimates of snow and firn density are required for satellite-altimetry-based
retrievals of ice sheet mass balance that rely on volume-to-mass
conversions. Therefore, biases and errors in presently used density models
confound assessments of ice sheet mass balance and by extension ice sheet
contribution to sea level rise. Despite this importance, most contemporary
firn densification models rely on simplified semi-empirical methods, which are
partially reflected by significant modeled density errors when compared to
observations. In this study, we present a new drifting-snow compaction scheme
that we have implemented into SNOWPACK, a physics-based land surface snow
model. We show that our new scheme improves existing versions of SNOWPACK
by increasing simulated near-surface (defined as the top 10 m) density
to be more in line with observations (near-surface bias reduction from −44.9
to −5.4 kg m−3). Furthermore, we demonstrate high-quality
simulation of near-surface Antarctic snow and firn density at 122 observed
density profiles across the Antarctic ice sheet, as indicated by reduced model
biases throughout most of the near-surface firn column when compared to two
semi-empirical firn densification models (SNOWPACK mean
bias=-9.7 kg m−3, IMAU-FDM mean bias=-32.5 kg m−3, GSFC-FDM mean bias=15.5 kg m−3). Notably, our analysis is restricted to the
near surface where firn density is most variable due to accumulation and
compaction variability driven by synoptic weather and seasonal climate
variability. Additionally, the GSFC-FDM exhibits lower mean density bias from
7–10 m (SNOWPACK bias=-22.5 kg m−3, GSFC-FDM
bias=10.6 kg m−3) and throughout the entire
near surface at high-accumulation sites (SNOWPACK bias=-31.4 kg m−3, GSFC-FDM bias=-4.7 kg m−3).
However, we found that the performance of SNOWPACK did not degrade when
applied to sites that were not included in the calibration of semi-empirical
models. This suggests that SNOWPACK may possibly better represent firn
properties in locations without extensive observations and under future
climate scenarios, when firn properties are expected to diverge from their
present state.