A large-scale numerical model for computing isochrone geometry

Hindmarsh, Richard C.A. ORCID:; Leysinger Vieli, Gwendolyn J.-M.C.; Parrenin, Frédéric. 2009 A large-scale numerical model for computing isochrone geometry. Annals of Glaciology, 50 (51). 130-140.

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


A finite-difference model for the calculation of radar layer geometries in large ice masses is presented. Balance velocities are used as coefficients in the age equation and in the heat equation. Solution of the heat equation allows prediction of sliding areas and computation of basal melt rates. Vertical distributions of velocity are parameterized using shape functions. These can be set uniformly, or allowed to vary in space according to the distribution of sliding. The vertical coordinate can either be uniformly distributed within the thickness of the ice, or be uniformly distributed within the flux. The finite-difference scheme results in a large set of linear equations. These are solved using a nested factorization preconditioned conjugate gradient scheme. The convergence properties of some other iteration solution schemes are studied. The output is computations of age and temperature assuming steady state, in large ice masses at high resolution. Age calculations are used to generate isochrones which show the best fit to observed layers. Comparisons with analytical solutions are made, and the influence of the order of the finite-difference approximation and the choice of vertical coordinate on solution accuracy is considered.

Item Type: Publication - Article
Digital Object Identifier (DOI):
Programmes: BAS Programmes > Polar Science for Planet Earth (2009 - ) > Ice Sheets
ISSN: 0260-3055
NORA Subject Terms: Glaciology
Earth Sciences
Date made live: 02 Nov 2010 14:08 +0 (UTC)

Actions (login required)

View Item View Item

Document Downloads

Downloads for past 30 days

Downloads per month over past year

More statistics for this item...