Search

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 ge

The British and Irish Ice Sheet (BIIS) was highly dynamic during the Late Quaternary, with considerable regional differences in the timing and extent of its change. This was reflected in equally variable offshore ice‐rafted debris (IRD) records. Here we reconcile these two records using the FRUGAL intermediate complexity iceberg–climate model, with varying BIIS catchment‐level iceberg fluxes, to simulate change in IRD origin and magnitude along the western European margin at 1000‐year time steps during the height of the last BIIS glaciation (31–6 ka bp). This modelled IRD variability is compared with existing IRD records from the deep ocean at five cores along this margin. There is general agreement of the temporal and spatial IRD variability between observations and model through this period. The Porcupine Bank off northwestern Ireland was confirmed by the modelling as a major dividing line between sites possessing exclusively northern or southern source regions for offshore IRD. During Heinrich events 1 and 2, the cores show evidence of a proportion of North American IRD, more particularly to the south of the British Isles. Modelling supports this southern bias for likely Heinrich impact, but also suggests North American IRD will only reach the British margin in unusual circumstances.

The findings of BRITICE-CHRONO Transect 2 through the North Sea Basin and eastern England are reported. We define ice-sheet marginal oscillation between -31 and 16 ka, with seven distinctive former ice-sheet limits (L1-7) constrained by Bayesian statistical analysis. The southernmost limit of the North Sea Lobe is recorded by the Bolders Bank Formation (L1; 25.8-24.6 ka). L2 represents ice-sheet oscillation and early retreat to the northern edge of the Dogger Bank (23.5-22.2 ka), with the Garret Hill Moraine in north Norfolk recording a significant regional readvance to L3 at 21.5-20.8 ka. Ice-marginal oscillations at similar to 26-21 ka resulted in L1, L2 and L3 being partially to totally overprinted. Ice-dammed lakes related to L1-3, including Lake Humber, are dated at 24.1-22.3 ka. Ice-sheet oscillation and retreat from L4 to L5 occurred between 19.7 and 17.3 ka, with grounding zone wedges marking an important transition from terrestrial to marine tidewater conditions, triggered by the opening of the Dogger Lake spillway between 19.9 and 17.5 ka. L6 relates to ice retreat under glacimarine conditions and final ice retreat into the Firth of Forth by 15.8 ka. L7 (-15 ka) represents an ice retreat from Bosies Bank into the Moray Firth. Copyright (C) 2021 The Authors Journal of Quaternary Science Published by John Wiley & Sons Ltd

Palaeo-ice sheets are important analogues for understanding contemporary ice sheets, offering a record of ice sheet behaviour that spans millennia. There are two main approaches to reconstructing palaeo-ice sheets. Empirical reconstructions use the available glacial geological and chronological evidence to estimate ice sheet extent and dynamics but lack direct consideration of ice physics. In contrast, numerically modelled simulations implement ice physics, but often lack direct quantitative comparison with empirical evidence. Despite being long identified as a fruitful scientific endeavour, few ice sheet reconstructions attempt to reconcile the empirical and model-based approaches. To achieve this goal, model-data comparison procedures are required. Here, we compare three numerically modelled simulations of the former British–Irish Ice Sheet with the following lines of evidence: (a) position and shape of former margin positions, recorded by moraines; (b) former ice-flow direction and flow-switching, recorded by flowsets of subglacial bedforms; and (c) the timing of ice-free conditions, recorded by geochronological data. These model–data comparisons provide a useful framework for quantifying the degree of fit between numerical model simulations and empirical constraints. Such tools are vital for reconciling numerical modelling and empirical evidence, the combination of which will lead to more robust palaeo-ice sheet reconstructions with greater explicative and ultimately predictive power.

Ice stream bed topography contains key evidence for the ways ice streams interact with, and are potentially controlled by, their beds. Here we present the first application of two–dimensional Fourier analysis to 22 marine and terrestrial topographies from 5 regions in Antarctica and Canada, with and without mega-scale glacial lineations (MSGLs). We find that the topography of MSGL-rich ice stream sedimentary beds is characterized by multiple, periodic wavelengths between 300 and 1200 m and amplitudes from decimeters to a few meters. This periodic topography is consistent with the idea that instability is a key element to the formation of MSGL bedforms. Dominant wavelengths vary among locations and, on one paleo ice stream bed, increase along the direction of ice flow by 1.7±0.52% km-1. We suggest that these changes are likely to reflect pattern evolution via downstream wavelength coarsening, even under potentially steady ice stream geometry and flow conditions. The amplitude of MSGLs is smaller than that of other fluvial and glacial topographies, but within the same order of magnitude. However, MSGLs are a striking component of ice stream beds because the topographic amplitude of features not aligned with ice flow is reduced by an order of magnitude relative to those oriented with the flow direction. This study represents the first attempt to automatically derive the spectral signatures of MSGLs. It highlights the plausibility of identifying these landform assemblages using automated techniques and provides a benchmark for numerical models of ice stream flow and subglacial landscape evolution.