The coupled ice sheet-Earth system model Bern3D v3.0

Understanding climate variability from millennial to glacial–interglacial time scales remains challenging due to the complex and nonlinear feedbacks between ice, ocean, sediments, biosphere, and atmosphere. Complex climate models generally struggle to dynamically and comprehensively simulate such long time periods as a result of the large computational costs. Here, we therefore coupled a dynamical ice sheet model to the Bern3D Earth system model of intermediate complexity, which allows for simulating multiple glacial–interglacial cycles. The performance of the model is first validated against modern observations and its response to abrupt perturbations, such as atmospheric CO2 changes and North Atlantic freshwater hosing, is investigated. To further test the fully coupled model, the climate evolution over the entire last glacial cycle is explored in a transient simulation forced by variations in the orbital configuration and greenhousegases and aerosols. The model simulates global mean surface temperature in fair agreement with reconstructions, exhibiting a gradual cooling trend since the last interglacial that is interrupted by two more rapid cooling events during the early Marine Isotope Stage (MIS) 4 and Last Glacial Maximum (LGM). Simulated Northern Hemispheric ice sheets show pronounced variability on orbital time scales, and ice volume more than doubles from MIS3 to the LGM in good agreement with recent sea level reconstructions. At the LGM, the Atlantic overturning has a strength of about 14 Sv (1 Sv ; 10 6 m 3 s 21 ), which is a reduction by about one-quarter compared to the preindustrial. We thus demonstrate that the new coupled model is able to simulate large-scale aspects of glacial–interglacial cycles