Effects of spatial discretization in ice-sheet modelling using the shallow-ice approximation

Author： van den Berg J. van de Wal R.S.W. Oerlemans J.

Publisher： International Glaciological Society

ISSN： 1727-5652

Source： Journal of Glaciology, Vol.52, Iss.176, 2006-01, pp. : 89-98

Disclaimer: Any content in publications that violate the sovereignty, the constitution or regulations of the PRC is not accepted or approved by CNPIEC.

Previous Menu Next

Abstract

This paper assesses a two-dimensional, vertically integrated ice model for its numerical properties in the calculation of ice-sheet evolution on a sloping bed using the shallow-ice approximation. We discuss the influence of initial conditions and individual model parameters on the model's numerical behaviour, with emphasis on varying spatial discretizations. The modelling results suffer badly from numerical problems. They show a strong dependence on gridcell size and we conclude that the widely used gridcell spacing of 20 km is too coarse. The numerical errors are small in each single time-step, but increase non-linearly over time and with volume change, as a result of feedback of the mass balance with height. We propose a new method for the calculation of the surface gradient near the margin, which improves the results significantly. Furthermore, we show that we may use dimension analysis as a tool to explain in which situations numerical problems are to be expected.