A high-resolution model of the 100 ka ice-age cycle

Abstract
Significant improvements to the representation of climate forcing and mass-balance response in a coupled two-dimensional global energy balance climate model (EBM) and vertically integrated ice-sheet model (ISM) have led to the prediction of an ice-volume chronology for the most recent ice-age cycle of the Northern Hemisphere that is close to that inferred from the geological record. Most significant is that full glacial termination is delivered by the model without the need for new physical ingredients. In addition, a relatively close match is achieved between the Last Glacial Maximum (LGM) model ice topography and that of the recently-described ICE-4G reconstruction. These results suggest that large-scale climate system reorganization is not required to explain the main variations of the North American (NA) ice sheets over the last glacial cycle. Lack of sea-ice and marine-ice dynamics in the model leaves the situation over the Eurasian (EA) sector much more uncertain.The incorporation of a gravitationally self-consistent description of the glacial isostatic adjustment process demonstrates that the NA and EA bedrock responses can be adequately represented by simpler damped-relaxation models with characteristic time-scales of 3–5ka and 5 ka, respectively. These relaxation times agree with those independently inferred on the basis of postglacial relative sea-level histories.