The last W elsh I ce C ap: Part 1 – Modelling its evolution, sensitivity and associated climate

A high‐resolution, three‐dimensional, thermomechanical ice‐flow model is used to investigate the glaciodynamics of the L ast G lacial M aximum W elsh I ce C ap – a large, independent ice centre of the B ritish– I rish I ce S heet. The model uses higher‐order physics to solve longitudinal stresses, and is coupled to climate via a distributed, positive degree‐day mass‐balance scheme. A suite of model experiments driven by the GISP 2 δ 18 O curve was initiated from a climatic optimum at 38.3 ka BP through to the D evensian/ H olocene boundary to identify an icecap configuration compatible with available empirical evidence. An enhanced cooling from present of 11.85° C and strong precipitation suppression are required between 27.4 and 23.5 ka BP for the modelled icecap to attain well‐established empirical limits, a scenario probably associated with H einrich E vent‐2 and the potential collapse of thermohaline circulation in the N orth A tlantic. The experiments indicate ice‐dispersal centres located in N orth and M id W ales, the latter being essential for forcing ice southwards of the B recon B eacons during the Last Glacial Maximum. Deglaciation of the Welsh Ice Cap was relatively rapid, occurring within one millennium. Dynamic stability is governed largely by the dominance and vigour with which fast‐flowing outlet glaciers drain the icecap interior, which in turn are linked to variations in the climatic forcing. The distribution of permanently cold‐based ice across the uplands and summits indicates the probable preservation of relict landscapes in these areas throughout the full glacial cycle.