The formation and rapid drainage of supraglacial lakes (SGL) influences the mass balance and dynamics of the Greenland Ice Sheet (GrIS). Although SGLs are expected to spread inland during the 21st century due to atmospheric warming, less is known about their future spatial distribution and volume. We use GrIS surface elevation model and regional climate model outputs to show that at the end of the 21st century (2070–2099) approximately 9.8 ± 3.9 km3 (+113% compared to 1980‐2009) and 12.6 ± 5 km3 (+174%) of meltwater could be stored in SGLs under moderate and high representative concentration pathways (RCP 4.5 and 8.5), respectively. The largest increase is expected in the northeastern sector of the GrIS (191% in RCP 4.5 and 320% in RCP 8.5), whereas in west Greenland, where the most SGLs are currently observed, the future increase will be relatively moderate (55% in RCP 4.5 and 68% in RCP 8.5).

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.

Predicting the amplitude and distribution of surface undulations on ice sheets and glaciers is useful because of their influence on surface mass and energy balance, atmospheric boundary layer processes, and supraglacial meltwater routing. We develop an approximate method of calculating the surface elevation response due to spatial perturbations in basal topography and slipperiness, on two‐dimensional flow sections whose thickness, surface slope, and basal slip ratio vary along flow. Our main result is an integral expressing nonuniform transfer of basal variability to the surface. It uses published Fourier transfer functions derived through perturbing plane‐slab Stokes flow but circumvents the need to subwindow the spatial domain to estimate the response. We test the method on ice flow synthesized by a finite‐element model of Stokes flow with constant viscosity and known basal topography and slipperiness perturbations; in this case, it predicts the observed size and shape of the surface undulations well, capturing more than 90% of their variance. Application of the method to the central flowline of Columbia Glacier, Alaska, and a flowline on the Greenland Ice Sheet ending on Nordenskiöld Glacier, using knowledge of the approximate bed topography and ignoring the unknown slipperiness forcing, yields less faithful prediction of their surface undulations (40–50% of their variance) but demonstrates the method's potential to reproduce their qualitative features. We discuss the factors limiting the method's performance on real flowlines.

Ice flow can transfer variations in basal topography and basal slipperiness to the ice surface. Recent developments in this theory have made it possible to conduct numerical experiments to predict mesoscale surface topographical undulations and surface relief on an ice sheet-scale. Focussing here on the contemporary Greenland Ice Sheet (GrIS), we demonstrate that the theory can be used to predict the surface relief of the ice sheet from bed topography, ice thickness and basal slip ratio datasets. In certain regions of the GrIS our approach overestimates, while in others underestimates, the observed surface relief. The magnitude and spatial pattern of these mismatches correspond with the theory's limitations and known uncertainties in the bed topography and basal slip ratio datasets. Our prediction experiment establishes that the first-order control on GrIS surface relief is basal topography modulated by ice thickness, surface slope and basal slip ratio. Additional analyses show that the surface relief, which is controlled by the bed-to-surface transfer of basal topography, preconditions the large scale spatial structure of surface drainage, with other factors such as surface runoff modulating the actual drainage system through influencing the temporal evolution of meltwater features. It follows that the spatial structure of surface drainage depends strongly on the transfer of basal topography to the ice surface. These findings represent an important step toward investigating and understanding the net long-term (>10 2 years) effect of surface drainage on ice sheet mass balance and dynamics during deglaciation events.

Mt. Ojos del Salado (6893 m a.s.l.) lies within the Andean Arid Diagonal, on the Chilean-Argentinean border. Due to the extremely arid climate, surface ice is not widespread on Mt. Ojos del Salado and at similar high-altitude massifs in the region, although ice-bearing permafrost might be present. However, the thermal regime of the ground is relatively unknown in the region, especially outside of rock glaciers at high elevations north of 30°S. To study ground thermal regimes, in-situ shallow ground temperature and snow coverage from satellite imagery have been surveyed for four years (2012-2016) at six sites between the elevations of 4200-6893 m a.s. l. on Mt. Ojos del Salado (27°07′S, 68°32′W). According to the ground temperature and snow coverage data at the six monitoring sites, the presence of permafrost is unlikely below 4550 m a.s.l. but likely above 5250 m a.s.l. on Mt. Ojos del Salado.In addition, the active layer becomes extremely thin around 6750 m a.s.l.