Simulations of Northern Hemisphere ice-sheet retreat: sensitivity to physical mechanisms involved during the Last Deglaciation

A 3-dimensional thermomechanical ice-sheet model is used to simulate the evolution of the geometry of Northern Hemisphere ice sheets through the Last Deglaciation. The ice-sheet model is forced by a time-evolving climatology provided by the linear interpolation through time of climate snapshots simulated by the LMD5.3 atmospheric general circulation model (AGCM) at different periods of the Last Deglaciation (21, 15, 9, 6 and 0 kyr BP). The AGCM is driven by insolation, atmospheric CO2 content, ice-sheet configuration and sea surface temperatures. Although our approach is able to produce the complete continental ice retreat, our simulated deglaciation presents a phase-lag with reconstructions based on observational evidences. This suggests that physical mechanisms related to climate forcing and/or ice-sheet internal dynamics are not properly represented. The influence of millennial-scale forcing, feedback mechanisms between ice-sheet elevation and surface mass balance and parameterization of the ice flow is also tested through a set of sensitivity experiments. The rapid variability has a strong impact on the evolution of the ice volume because of nonlinear effects in temperature-mass balance relationships. Fennoscandia appears to be strongly sensitive to the small-scale ice-sheet instability. Both ice sheets are to some extent sensitive to an increased basal sliding. (C) 2001 Elsevier Science Ltd. All rights reserved.