Abstract. The last glacial period (LGP; ca.110–10 ka BP) was marked by the existence of two types of abrupt climatic changes, Dansgaard-Oeschger (DO) and Heinrich (H) events. Although the mechanisms behind these are not fully understood, it is generally accepted that the presence of ice sheets played an important role in their occurrence. While an important effort has been made to investigate the dynamics and evolution of the Laurentide Ice Sheet (LIS) during this period, the Eurasian Ice Sheet (EIS) has not received much attention, in particular from a modeling perspective. However, meltwater discharge from this and other ice sheets surrounding the Nordic Seas is often implied as a potential cause of ocean instabilities that lead to glacial abrupt climate changes. Thus, a better understanding of its variations during the LGP is important to understand its role in glacial abrupt climate changes. Here we investigate the response of the EIS to millennial-scale climate variability during the LGP. We use a hybrid, three-dimensional, thermomechanical ice-sheet model that includes ice shelves and ice streams. The model is forced offline through a novel perturbative approach that includes the effect of both atmospheric and oceanic variations and provides a more realistic treatment of millennial-scale climatic variability than conventional methods. Our results show that the EIS responds with enhanced ice discharge in phase with interstadial warming in the North Atlantic when forced with surface ocean temperatures. Conversely, when subsurface ocean temperatures are used, enhanced ice discharge occurs both during stadials and at the beginning of the interstadials. Separating the atmospheric and oceanic effects demonstrates the major role of the ocean in controlling the dynamics of the EIS on millennial time scales. While the atmospheric forcing alone is only able to produce modest iceberg discharges, warming of the ocean leads to higher rates of iceberg discharges as a result of relatively strong basal melting at the margins of the ice sheet. Together with previous work, our results provide a consistent explanation for the response of the LIS and the EIS to glacial abrupt climate changes, and highlight the need for stronger constraints on the local North Atlantic behavior in order to improve our understanding of the ice sheet's glacial dynamics.
Millennial-scale iceberg discharges in the Irminger Basin during the Last Glacial Period: Relationship with the Heinrich events and environmental settings
