Coupled ice sheet–climate modeling under glacial and pre-industrial boundary conditions

In the standard Paleoclimate Modelling Intercomparison Project (PMIP) experiments, the Last Glacial Maximum (LGM) is modeled in quasi-equilibrium with atmosphere–ocean–vegetation general circulation models (AOVGCMs) with prescribed ice sheets. This can lead to inconsistencies between the modeled climate and ice sheets. One way to avoid this problem would be to model the ice sheets explicitly. Here, we present the first results from coupled ice sheet–climate simulations for the pre-industrial times and the LGM. Our setup consists of the AOVGCM ECHAM5/MPIOM/LPJ bidirectionally coupled with the Parallel Ice Sheet Model (PISM) covering the Northern Hemisphere. The results of the pre-industrial and LGM simulations agree reasonably well with reconstructions and observations. This shows that the model system adequately represents large, non-linear climate perturbations. A large part of the drainage of the ice sheets occurs in ice streams. Most modeled ice stream systems show recurring surges as internal oscillations. The Hudson Strait Ice Stream surges with an ice volume equivalent to about 5 m sea level and a recurrence interval of about 7000 yr. This is in agreement with basic expectations for Heinrich events. Under LGM boundary conditions, different ice sheet configurations imply different locations of deep water formation