Mass balance, flow and subglacial processes of a modelled Younger Dryas ice cap in Scotland

We use an empirically validated high-resolution three-dimensional ice-sheet model to investigate the mass-balance regime, flow mechanisms and subglacial characteristics of a simulated Younger Dryas Stadial ice cap in Scotland, and compare the resulting model forecasts with geological evidence. Input data for the model are basal topography, a temperature forcing derived from GRIP δ18O fluctuations and a precipitation distribution interpolated from modern data. The model employs a positive-degree-day scheme to calculate net mass balance within a domain of 112 500 km2, which, under the imposed climate, gives rise to an elongate ice cap along the axis of the western Scottish Highlands. At its maximum, the ice cap is dynamically and thermally zoned, reflecting topographic and climatic controls, respectively. In order to link these palaeoglaciological conditions to geological interpretations, we calculate the relative balance between sliding and creep within the simulated ice cap, forecast areas of the ice cap with the greatest capacity for basal erosion and predict the likely pattern of subglacial drainage. We conclude that ice flow in central areas of the ice cap is largely due to internal deformation, and is associated with geological evidence of landscape preservation. Conversely, the distribution of streamlined landforms is linked to faster-flowing ice whose velocity is predominantly the result of basal sliding. The geometry of the main ice mass focuses subglacial erosion in the mid-sections of topographic troughs, and produces glaciohydraulic gradients that favour subglacial drainage through low-order arterial routes.