Annals of Glaciology 46 2007 97 Fast computation of a viscoelastic deformable Earth model for ice-sheet simulations

• Main Authors: Craig S. Lingle, Jed Brown
• Other Authors: The Pennsylvania State University CiteSeerX Archives
• Format: Text
• Language: English
• Subjects: Ice Sheet
• Online Access: http://citeseerx.ist.psu.edu/viewdoc/summary?doi=10.1.1.429.7792; http://www.igsoc.org/annals/46/a46a130.pdf

ABSTRACT. The model used by Lingle and Clark (1985) to approximate the deformation of the Earth under a single ice stream is adapted to the purposes of continent-scale ice-sheet simulation. The model combines a layered elastic spherical Earth (Farrell, 1972) with a viscous half-space overlain by an elastic plate lithosphere (Cathles, 1975). For the half-space model we identify a new mathematical formulation, essentially a time-dependent partial differential equation, which generalizes and improves upon the standard elastic plate lithosphere with relaxing asthenosphere model widely used in ice-sheet simulation. The new formulation allows a significantly faster numerical strategy, a spectral collocation method based directly on the fast Fourier transform. We verify this method by comparing to an integral formula for a disk load. We also demonstrate that the magnitudes of numerical errors made in approximating coupled ice-flow/Earth-deformation systems are significantly smaller than pairwise differences between several Earth models. Our implementation of the Lingle and Clark (1985) model offers important features of spherical, layered, self-gravitating, viscoelastic Earth models without the computational expense. 1.