Annals of Glaciology 46 2007 97 Fast computation of a viscoelastic deformable Earth model for ice-sheet simulations
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...
Main Authors: | , |
---|---|
Other Authors: | |
Format: | Text |
Language: | English |
Subjects: | |
Online Access: | http://citeseerx.ist.psu.edu/viewdoc/summary?doi=10.1.1.429.7792 http://www.igsoc.org/annals/46/a46a130.pdf |
Summary: | 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. |
---|