IngentaConnect Fast computation of a viscoelastic deformable Earth model for ice...

Authors: Bueler, E.d.; Lingle, Craig S.; Brown, Jed

Source: Annals of Glaciology, Volume 46, Number 1, October 2007 , pp. 97-105(9)

Publisher: International Glaciological Society

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.

The full text electronic article is available for purchase. You will be able to download the full text electronic article after payment.

$32.46 plus tax Refund Policy

Fast computation of a viscoelastic deformable Earth model for ice-sheet simulations