We use the Shreve hydraulic potential equation as a simplified approach to investigate potential subglacial lake locations and meltwater drainage pathways beneath the Antarctic and Greenland ice sheets. We validate the method by demonstrating its ability to recall the locations of >60% of the known subglacial lakes beneath the Antarctic Ice Sheet. This is despite uncertainty in the ice-sheet bed elevation and our simplified modelling approach. However, we predict many more lakes than are observed. Hence we suggest that thousands of subglacial lakes remain to be found. Applying our technique to the Greenland Ice Sheet, where very few subglacial lakes have so far been observed, recalls 1607 potential lake locations, covering 1.2% of the bed. Our results will therefore provide suitable targets for geophysical surveys aimed at identifying lakes beneath Greenland. We also apply the technique to modelled past ice-sheet configurations and find that during deglaciation both ice sheets likely had more subglacial lakes at their beds. These lakes, inherited from past ice-sheet configurations, would not form under current surface conditions, but are able to persist, suggesting a retreating ice-sheet will have many more subglacial lakes than advancing ones. We also investigate subglacial drainage pathways of the present-day and former Greenland and Antarctic ice sheets. Key sectors of the ice sheets, such as the Siple Coast (Antarctica) and NE Greenland Ice Stream system, are suggested to have been susceptible to subglacial drainage switching. We discuss how our results impact our understanding of meltwater drainage, basal lubrication and ice-stream formation.

Cosmogenic-nuclide concentrations in subglacial bedrock cores show that the West Antarctic Ice Sheet (WAIS) at a site between Thwaites and Pope glaciers was at least 35 m thinner than present in the past several thousand years and then subsequently thickened. This is important because of concern that present thinning and grounding line retreat at these and nearby glaciers in the Amundsen Sea Embayment may irreversibly lead to deglaciation of significant portions of the WAIS, with decimeter- to meter-scale sea level rise within decades to centuries. A past episode of ice sheet thinning that took place in a similar, although not identical, climate was not irreversible. We propose that the past thinning–thickening cycle was due to a glacioisostatic rebound feedback, similar to that invoked as a possible stabilizing mechanism for current grounding line retreat, in which isostatic uplift caused by Early Holocene thinning led to relative sea level fall favoring grounding line advance.

Widespread existing geological records from above the modern ice sheet surface and outboard of the current ice margin show that the Antarctic Ice Sheet (AIS) was much more extensive at the Last Glacial Maximum (∼ 20 ka) than at present. However, whether it was ever smaller than present during the last few millennia, and (if so) by how much, is known only for a few locations because direct evidence lies within or beneath the ice sheet, which is challenging to access. Here, we describe how retreat and readvance (henceforth “readvance”) of AIS grounding lines during the Holocene could be detected and quantified using subglacial bedrock, subglacial sediments, marine sediment cores, relative sea-level (RSL) records, geodetic observations, radar data, and ice cores. Of these, only subglacial bedrock and subglacial sediments can provide direct evidence for readvance. Marine archives are of limited utility because readvance commonly covers evidence of earlier retreat. Nevertheless, stratigraphic transitions documenting change in environment may provide support for direct evidence from subglacial records, as can the presence of transgressions in RSL records, and isostatic subsidence. With independent age control, ice structure revealed by radar can be used to infer past changes in ice flow and geometry, and therefore potential readvance. Since ice cores capture changes in surface mass balance, elevation, and atmospheric and oceanic circulation that are known to drive grounding line migration, they also have potential for identifying readvance. A multidisciplinary approach is likely to provide the strongest evidence for or against a smaller-than-present AIS in the Holocene.

The rapidly retreating Thwaites and Pine Island glaciers together dominate present-day ice loss from the West Antarctic Ice Sheet and are implicated in runaway deglaciation scenarios. Knowledge of whether these glaciers were substantially smaller in the mid-Holocene and subsequently recovered to their present extents is important for assessing whether current ice recession is irreversible. Here we reconstruct relative sea-level change from radiocarbon-dated raised beaches at sites immediately seawards of these glaciers, allowing us to examine the response of the earth to loading and unloading of ice in the Amundsen Sea region. We find that relative sea level fell steadily over the past 5.5 kyr without rate changes that would characterize large-scale ice re-expansion. Moreover, current bedrock uplift rates are an order of magnitude greater than the rate of long-term relative sea-level fall, suggesting a change in regional crustal unloading and implying that the present deglaciation may be unprecedented in the past ~5.5 kyr. While we cannot preclude minor grounding-line fluctuations, our data are explained most easily by early Holocene deglaciation followed by relatively stable ice positions until recent times and imply that Thwaites and Pine Island glaciers have not been substantially smaller than present during the past 5.5 kyr.