Supraglacial lakes (SGLs) form annually on the Greenland ice sheet 1,2 and, when they drain, their discharge enhances ice-sheet flow 3 by lubricating the base 4 and potentially by warming the ice 5 . Today, SGLs tend to form within the ablation zone, where enhanced lubrication is o set by e cient subglacial drainage 6,7 . However, it is not clear what impact a warming climate will have on this arrangement. Here, we use an SGL initiation and growth 8 model to show that lakes form at higher altitudes as temperatures rise, consistent with satellite observations 9 . Our simulations show that in southwest Greenland, SGLs spread 103 and 110 km further inland by the year 2060 under moderate (RCP 4.5) and extreme (RCP 8.5) climate change scenarios, respectively, leading to an estimated 48-53% increase in the area over which they are distributed across the ice sheet as a whole. Up to half of these new lakes may be large enough to drain, potentially delivering water and heat to the ice-sheet base in regions where subglacial drainage is ine cient. In such places, ice flow responds positively to increases in surface water delivered to the bed through enhanced basal lubrication 4,10,11 and warming of the ice 5 , and so the inland advance of SGLs should be considered in projections of ice-sheet change.The volume of water stored in SGLs on the surface of the Greenland ice sheet is determined by the presence of depressions in the local terrain 2 , by the amount of runo 8 (melt water plus rain minus refreezing in the snowpack) and by lake drainage 3 . It is estimated that 13% of Greenland's SGLs drain on timescales of the order of a few hours 12 , often by the creation of moulins as waterfilled fractures propagate through the full thickness of the ice sheet (termed hydro-fracture) 13 . SGLs act as a source of en-and subglacial water when they drain and afterwards, the moulin acts as a conduit allowing runo to pass between the ice-sheet surface and base 1,3 . Satellite and ground-based observations show a correlation between the degree of runo and the rate of ice motion 4,6,7 ; however, there are known spatial and temporal variations in the magnitude and sign of this relationship. For example, near the ice-sheet margin, lower annual ice speeds have been recorded in years of high melting 6,7 but further inland-at higher elevations-the reverse seems to be the case 4,11 . This dichotomy can be attributed to an abundance of melt water at the margin, enabling the evolution of e cient subglacial drainage early in the melt season 6,10 , and thicker ice and less water farther inland hindering the development of an e cient evacuation system 14,15 . In addition to their impact on basal sliding, draining SGLs, and moulins that persist post-drainage, can exert a local warming as relatively warm water passes through the colder