The natural fluctuations of firn densification and their effect on the geodetic determination of ice sheet mass balance

A one-dimensional, numerical model of time-evolving firn densification was used to simulate the response of the density profile through an ice sheet to changes in the temperature, density and accumulation rate at the surface. The equilibrium response of the model was compared with ice-core density profiles from Byrd, Antarctica and Site 2, Greenland, and the model predicted the density to within 10% of both cores.. The response of the model to step-wise changes and random fluctuations in the surface boundary conditions was investigated. The standard deviation of elevation changes as a function of observation interval was computed. These changes were found to be small in comparison with the magnitude of present uncertainties in the mass balances of the Antarctic and Greenland Ice Sheets. It was concluded that, in the dry snow zones, natural variability in the densification will not prevent the geodetic determination of ice sheet mass balance from improving upon current estimates. Uncertainty in the constitutive equation for snow and firn is the dominant source of error in the calculations.