The configuration, sensitivity and rapid retreat of the Late Weichselian Icelandic ice sheet

The fragmentary glacial-geological record across the Icelandic continental shelf has hampered reconstruction of the volume, extent and chronology of the Late Weichselian ice sheet particularly in key offshore zones. Marine geophysical data collected over the last two decades reveal that the ice sheet likely attained a continental shelf-break position in all sectors during the Last Glacial Maximum, though its precise timing and configuration remains largely unknown. Within this context, we review the available empirical evidence and use a well-constrained three-dimensional thermomechanical model to investigate the drivers of an extensive Late Weichselian Icelandic ice-sheet, its sensitivity to environmental forcing, and phases of deglaciation. Our reconstruction attains the continental shelf break across all sectors with a total ice volume of 5.96 × 10⁵ km³ with high precipitation rates being critical to forcing extensive ice sheet flow offshore. Due to its location astride an active mantle plume, a relatively fast and dynamic ice sheet with a low aspect ratio is maintained. Our results reveal that once initial ice-sheet retreat was triggered through climate warming at 21.8 ka BP, marine deglaciation was rapid and accomplished in all sectors within c. 5 ka at a mean rate of 71 Gt of mass loss per year. This rate of ice wastage is comparable to contemporary rates observed for the West Antarctic ice sheet. The ice sheet subsequently stabilised on shallow pinning points across the near shelf for two millennia, but abrupt atmospheric warming during the Bølling Interstadial forced a second, dramatic collapse of the ice sheet onshore with a net wastage of 221 Gt a^-1 over 750 years, analogous to contemporary Greenland rates of mass loss. Geothermal conditions impart a significant control on the ice sheet's transient response, particularly during phases of rapid retreat. Insights from this study suggests that large sectors of contemporary ice sheets overlying geothermally active regions, such as Siple Coast, Antarctica, and NE Greenland, have the potential to experience rapid phases of mass loss and deglaciation once initial retreat is initiated.