Glaciation and Deglaciation Mechanisms in a Coupled 2-dimensional Climate-ice-sheet Model

A two-dimensional model which links the atmosphere, the mixed layer of the ocean, the sea ice, the continents, the ice sheets and their underlying bedrock has been used to test the Milankovitch theory over the last glacial-interglacial cycle. It was found that the orbital variations alone can induce, in such a system, feed-backs sufficient to generate the low-frequency part of the climatic variations over the last 122 kyear. These simulated variations at the astronomical time-scale are broadly in agreement with ice volume and sea-level reconstructions independently obtained from geological data. Imperfections in the simulated climate were the insufficient southward extent of the ice sheets and the too small hemispheric cooling during the last glacial maximum. These deficiencies were partly remedied in a further experiment (Gallee and others, in press) by using the time-dependent CO2 atmospheric concentration given by the Vostok ice core in addition to the astronomical forcing. For this second experiment, the main mechanisms and feedbacks responsible for the glaciation and the deglaciation in the model are discussed here.