Ice-contact proglacial lakes are generally absent from numerical model simulations of glacier evolution, and their effects on ice dynamics and on rates of deglaciation remain poorly quantified. Using the BISICLES ice flow model, we analyzed the effects of an ice-contact lake on the Pukaki Glacier, New Zealand, during recession from the Last Glacial Maximum. The ice-contact lake produced a maximum effect on grounding line recession >4 times further and on ice velocities up to 8 times faster, compared to simulations of a land-terminating glacier forced by the same climate. The lake contributed up to 82% of cumulative grounding line recession and 87% of ice velocity during the first 300 years of the simulations, but those values decreased to just 6% and 37%, respectively, after 5,000 years. Numerical models that ignore lake interactions will, therefore, misrepresent the rate of recession especially during the transition of a land-terminating to a lake-terminating environment. Plain Language Summary Lakes form at the margins of glaciers as meltwater accumulates against hillsides and behind ridges of glacier debris. Lakes at a glacier terminus are known to affect its behavior. However, glaciers terminating into such lakes are usually absent from computer simulations, and the effects of these lakes on the rates of deglaciation and on glacier behavior are poorly quantified. In this study, we tested the effect of a lake on glacier recession under two different scenarios; a land-terminating versus a lake-terminating glacier. We used an ice flow model called BISICLES and applied it to what was once the Pukaki Glacier in New Zealand during the end of the last ice age. We found that the presence of a lake caused the glacier to recede more than 4 times further and it accelerated ice flow by up to 8 times when compared to the same glacier that terminated on land under the same climate. Our simulated lake processes predominantly influenced the glacier over decades to centuries rather than over millennia. We suggest, therefore, that simulations of glacier evolution ignoring glacial lakes will likely misrepresent the timing and rate of recession, especially during the transition from a land-terminating to a lake-terminating environment. Recent observations, both field-based (e.g., Watson et al., 2020) and via satellite imagery (e.g., King et al., 2019), have highlighted the spatiotemporal frequency and magnitude of changes in glacial lakes and the associated glaciers that feed meltwater to them. However, field-based measurements are limited to

The BRITICE‐CHRONO consortium of researchers undertook a dating programme to constrain the timing of advance, maximum extent and retreat of the British–Irish Ice Sheet between 31 000 and 15 000 years before present. The dating campaign across Ireland and Britain and their continental shelves, and across the North Sea included 1500 days of field investigation yielding 18 000 km of marine geophysical data, 377 cores of sea floor sediments, and geomorphological and stratigraphical information at 121 sites on land; generating 690 new geochronometric ages. These findings are reported in 28 publications including synthesis into eight transect reconstructions. Here we build ice sheet‐wide reconstructions consistent with these findings and using retreat patterns and dates for the inter‐transect areas. Two reconstructions are presented, a wholly empirical version and a version that combines modelling with the new empirical evidence. Palaeoglaciological maps of ice extent, thickness, velocity, and flow geometry at thousand‐year timesteps are presented. The maximum ice volume of 1.8 m sea level equivalent occurred at 23 ka. A larger extent than previously defined is found and widespread advance of ice to the continental shelf break is confirmed during the last glacial. Asynchrony occurred in the timing of maximum extent and onset of retreat, ranging from 30 to 22 ka. The tipping point of deglaciation at 22 ka was triggered by ice stream retreat and saddle collapses. Analysis of retreat rates leads us to accept our hypothesis that the marine‐influenced sectors collapsed rapidly. First order controls on ice‐sheet demise were glacio‐isostatic loading triggering retreat of marine sectors, aided by glaciological instabilities and then climate warming finished off the smaller, terrestrial ice sheet. Overprinted on this signal were second order controls arising from variations in trough topographies and with sector‐scale ice geometric readjustments arising from dispositions in the geography of the landscape. These second order controls produced a stepped deglaciation. The retreat of the British–Irish Ice Sheet is now the world’s most well‐constrained and a valuable data‐rich environment for improving ice‐sheet modelling.

Abstract. Uncertainties in future sea level projections are dominated by our limited understanding of the dynamical processes that control instabilities of marine ice sheets. The last deglaciation of the British–Irish Ice Sheet offers a valuable example to examine these processes. The Minch Ice Stream, which drained a large proportion of ice from the northwest sector of the British–Irish Ice Sheet during the last deglaciation, is constrained with abundant empirical data which can be used to inform, validate, and analyse numerical ice sheet simulations. We use BISICLES, a higher-order ice sheet model, to examine the dynamical processes that controlled the retreat of the Minch Ice Stream. We perform simplified experiments of the retreat of this ice stream under an idealised climate forcing to isolate the effect of marine ice sheet processes, simulating retreat from the continental shelf under constant “warm” surface mass balance and sub-ice-shelf melt. The model simulates a slowdown of retreat as the ice stream becomes laterally confined at the mouth of the Minch strait between mainland Scotland and the Isle of Lewis, resulting in a marine setting similar to many large tidewater glaciers in Greenland and Antarctica. At this stage of the simulation, the presence of an ice shelf becomes a more important control on grounded ice volume, providing buttressing to upstream ice. Subsequently, the presence of a reverse slope inside the Minch strait produces an acceleration in retreat, leading to a “collapsed” state, even when the climate returns to the initial “cold” conditions. Our simulations demonstrate the importance of the marine ice sheet instability and ice shelf buttressing during the deglaciation of parts of the British–Irish Ice Sheet. We conclude that geological data could be applied to further constrain these processes in ice sheet models used for projecting the future of contemporary ice sheets.

No abstract

This is a repository copy of Recent progress on combining geomorphological and geochronological data with ice sheet modelling, demonstrated using the last British-Irish Ice Sheet.