A complete three-dimensional thermo-mechanical ice-shect model for the entire Antarctic ice sheet, including an ice shelf, grounding line-dynamics and isostatic bed adjustment, is employed to simulate the response of the ice sheet during the last glacial-interglacial cycle with respect to changing environmental conditions. To do this, the Vostok temperature signal is used to force changes in surface temperature and accumulation rate and sea level prescribed by a piecewise linear sawtooth function. Model calculations started at 160 ka B.P. In line with glacial geological evidence, the most pronounced fluctuations are found in the West Antarctic ice sheet and appear to be essentially controlled by changes in eustatic sea level. Grounding occurs more readily in the Weddell Sea than in the Ross Sea and, due to the long time scales involved, the ice sheet does not reach its full glacial extent until 16 ka B.p. The concomitant disintegration of the West Antarctic ice sheet is triggered by a rise in sea level and takes around 6000 years to complete. The ice sheet then halts close to the present state and no collapse takes place. This Holocene deglaciation appears to have added 6–8 million km3 of ice to the world oceans, corresponding with an Antarctic contribution to world-wide sea level of 12–15 m.
Reconstruction of changes in the Amundsen Sea and Bellingshausen Sea sector of the West Antarctic Ice Sheet since the Last Glacial Maximum