Mass-conserving subglacial hydrology in the Parallel Ice Sheet Model v0.6

Published: Jun 1, 2015 by The PISM Authors

Thickness of subglacial water in "routing" model.

Thickness of subglacial water in “routing” model.

Title Mass-conserving subglacial hydrology in the Parallel Ice Sheet Model v0.6
Authors E. Bueler and W. van Pelt
Venue Geoscientific Model Development

This paper describes and tests a major extension of PISM, introduced in v0.6 and fully-supported in v0.7, a two-horizontal-dimension subglacial hydrology model which combines till with a distributed system of water-filled, linked cavities. This sub-model accomplishes three specific goals: (1) conservation of the mass of water, (2) simulation of spatially- and temporally-variable basal shear stress from physical mechanisms based on a minimal number of free parameters, and (3) convergence under grid refinement. Besides a broad approach to the source of the model equations and detailed attention to the implementation and testing of the numerics, this paper demonstrates the model at scale by modeling the whole Greenland ice sheet at 2 km horizontal resolution, with one million nodes in the hydrology grid. But the model is far from complete. It both takes a very conservative approach to coupling hydrology to a model for basal shear, and it does not include the physics which determines the location and evolution of subglacial conduits.


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.