Investigating the evolution of major Northern Hemisphere ice sheets during the last glacial-interglacial cycle

A 2.5-dimensional climate model of intermediate complexity, CLIMBER-2, fully coupled with the GREMLINS 3-D thermo-mechanical ice sheet model is used to simulate the evolution of major Northern Hemisphere ice sheets during the last glacial-interglacial cycle and to investigate the ice sheets responses to both insolation and atmospheric CO<sub>2</sub> concentration. This model reproduces the main phases of advance and retreat of Northern Hemisphere ice sheets during the last glacial cycle, although the amplitude of these variations is less pronounced than those based on sea level reconstructions. At the last glacial maximum, the simulated ice volume is 52.5×10<sup>15</sup> m<sup>3</sup> and the spatial distribution of both the American and Eurasian ice complexes is in reasonable agreement with observations, with the exception of the marine parts of these former ice sheets. <br> A set of sensitivity studies has also been performed to assess the sensitivity of the Northern Hemisphere ice sheets to both insolation and atmospheric CO<sub>2</sub>. Our results suggest that the decrease of summer insolation is the main factor responsible for the early build up of the North American ice sheet around 120 kyr BP, in agreement with benthic foraminifera δ<sup>18</sup>O signals. In contrast, low insolation and low atmospheric CO<sub>2</sub> concentration are both necessary to trigger a long-lasting glaciation over Eurasia