From cyclic ice streaming to Heinrich-like events: the grow-and-surge instability in the Parallel Ice Sheet Model

Published: Oct 1, 2017 by The PISM Authors

Title From cyclic ice streaming to Heinrich-like events: the grow-and-surge instability in the Parallel Ice Sheet Model
Authors J. Feldmann and A. Levermann
Venue The Cryosphere

When thermodynamics are coupled to basal sliding through the production of liquid water in a modeled basal till layer, PISM may exhibit cyclic, surge-like behavior. This paper studies such ice-discharge instabilities in PISM for a buttressed ice-sheet-stream-shelf system. When present the simulated surge cycles, which are believed to be fundamentally physical, have 1000 to 5000 year periods in the studied parameter ranges. This work extends prior work on cyclic PISM behavior (Bueler & Brown 2009, van Pelt & Oerlemans 2012) to marine ice sheet conditions, but the 1 km resolution model here is significantly more complete than in those studies. It includes stable0.7 version implementations of enthalpy-based thermodynamics, an improved till-layer hydrology model (Bueler & van Pelt 2015), a sub-grid scheme at the grounding line (Feldman et al 2014), and a buttressed ice shelf. The surface conditions are all constant, independent even of surface elevation, so the cyclic behavior arises purely internally in the ice and the basal layer. This paper includes a new and careful analysis of cyclic behavior, identifying three competing time-scales controlling such cycles. Results suggest that ice sheets of medium thickness may be more susceptible to surging than thin or thick ones. The mechanisms studied here may play a role in Heinrich events of the late-LGM Laurentide Ice Sheet and/or in ice-stream shutdown and reactivation observed in West Antarctica. This is the definitive study of the PISM parameter space for modeled ice stream-shelf systems.


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.