A numerical comparison of approximations to the Stokes equations used in ice sheet and glacier modeling

Hindmarsh, R.C.A. ORCID: https://orcid.org/0000-0003-1633-2416. 2004 A numerical comparison of approximations to the Stokes equations used in ice sheet and glacier modeling. Journal of Geophysical Research, 109 (F1), F01012. 15, pp. 10.1029/2003JF000065

Full text not available from this repository. (Request a copy)

Abstract/Summary

A computational analysis of the accuracy of different approximations to the Stokes equations for momentum balance used in ice sheet modeling is performed by solving a particular tractable form of the equations appropriate for small perturbations of the ice surface, describing the uniform flow of ice with a Glen rheology on an infinitely long and broad section. The approximants comprise the shallow ice approximation and various schemes for incorporating longitudinal stresses and, in one case, the horizontal gradient of the horizontal plane shear stresses. The simplifications lead to a vertically one-dimensional numerical problem, whose solution can be computed rapidly. The relaxation rate of perturbations as well as other response descriptors for the stable full system and approximants are compared. Compared with the shallow ice approximation, the inclusion of longitudinal stresses increases accuracy at shorter wavelengths, but accuracy is poor at wavelengths around or less than the ice sheet thickness. Even though analysis shows that the horizontal gradients of the horizontal plane shear stresses are of similar magnitude to longitudinal stress effects, computations show, in agreement with glaciological belief, that longitudinal stress effects are more significant and need to be corrected for first in practice. Two schemes, a multilayer scheme and a one-layer scheme, are particularly good and should be investigated further in cases where perturbations from uniformity are large. Some other apparently plausible approximations introduce nonphysical instabilities. New schemes need to be assessed in the way described in this paper before being used in real ice sheet models.

Item Type:	Publication - Article
Digital Object Identifier (DOI):	10.1029/2003JF000065
Programmes:	BAS Programmes > Antarctic Science in the Global Context (2000-2005) > Global Interactions of the Antarctic Ice Sheet
ISSN:	0148-0227
Additional Keywords:	ice sheets, mechanics, mathematical geophysics
NORA Subject Terms:	Glaciology Earth Sciences Mathematics
Date made live:	18 Jan 2012 10:47 +0 (UTC)
URI:	https://nora.nerc.ac.uk/id/eprint/12192

Actions (login required)

View Item

Document Downloads

More statistics for this item...