Error sources in basal yield stress inversions for Jakobshavn Isbræ, Greenland, derived from residual patterns of misfit to observations

Published: Nov 1, 2017 by The PISM Authors

Title Error sources in basal yield stress inversions for Jakobshavn Isbræ, Greenland, derived from residual patterns of misfit to observations
Authors M. Habermann and others
Venue J. Glaciol.

Glaciologists know little about the conditions at the beds of glaciers because these places are so hard to observe directly. Nonetheless most PISM runs need and use a model of the basal conditions. This includes a model of the basal topography, derived from ice-penetrating radar (though necessarily more than the “raw” reflections), and a model of the relation between basal sliding velocity and shear stress (e.g. a power law with uncertain parameters). Inverse modeling can examine the relationship between available surface observations and the conditions at the base. This paper assesses the intermediate quantities in an inverse modeling exercise. Such a calculation generates a (forward) model result for which a norm of the difference between the observations and the model result is minimized in some sense. This difference itself, the residual, contains information which is lost in computing its norm. These authors examine the residual, and other quantities as well, for meaning. There are several types of errors in glaciological inverse modeling which generates basal shear stress. These include observation errors (in the surface elevation and velocity), geometry errors, ice temperature errors, and “model errors”. The last are particularly interesting in the PISM context; a question addressed by this paper is “what are the manifestations of the simplifications in PISM in the sense of residual features?” A particular focus is on the difference between upstream and near-terminus regions of a deep-trough outlet glacier in the Greenland ice sheet.


Latest news

PISM 2.0 is out

PISM developers have been hard at work to bring you a brand new version of PISM, packed with new features. After years of development, PISM finally includes a Blatter solver, warranting a new major version: PISM 2.0.

Version 1.2

We are pleased to announce the release of the Parallel Ice Sheet Model (PISM) v1.2.

MPI-M Hamburg, Germany: open postdoc for coupled atmosphere-ocean-ice sheet model

The Max Planck Institute for Meteorology (MPI-M) contributes to the BMBF project “From the Last Interglacial to the Anthropocene: Modeling a Complete Glacial Cycle” (PalMod,, which aims at simulating the climate from the peak of the last interglacial up to the present using comprehensive Earth System Models. Phase II of this project has an open position Postdoctoral Scientist (W073). The successful candidate will be part of a local team performing and analysing long-term transient simulations covering the last glacial and the transition into the Holocene with an interactively coupled atmosphere-ocean-ice sheet model. Additionally, the candidate will contribute to the continued development of this model. The model system consists of the MPI-Earth system model, the ice sheet model PISM, and the solid-earth model VILMA.