Remote sensing is commonly used to monitor supraglacial lakes on the Greenland Ice Sheet (GrIS); however, most satellite records must trade off higher spatial resolution for higher temporal resolution (e.g. MODIS) or vice versa (e.g. Landsat). Here, we overcome this issue by developing and applying a dual-sensor method that can monitor changes to lake areas and volumes at high spatial resolution (10–30 m) with a frequent revisit time ( ∼ 3 days). We achieve this by mosaicking imagery from the Landsat 8 Operational Land Imager (OLI) with imagery from the recently launched Sentinel-2 Multispectral Instrument (MSI) for a ∼ 12 000 km² area of West Greenland in the 2016 melt season. First, we validate a physically based method for calculating lake depths with Sentinel-2 by comparing measurements against those derived from the available contemporaneous Landsat 8 imagery; we find close correspondence between the two sets of values (R² = 0.841; RMSE = 0.555 m). This provides us with the methodological basis for automatically calculating lake areas, depths, and volumes from all available Landsat 8 and Sentinel-2 images. These automatic methods are incorporated into an algorithm for Fully Automated Supraglacial lake Tracking at Enhanced Resolution (FASTER). The FASTER algorithm produces time series showing lake evolution during the 2016 melt season, including automated rapid ( ≤ 4 day) lake-drainage identification. With the dual Sentinel-2–Landsat 8 record, we identify 184 rapidly draining lakes, many more than identified with either imagery collection alone (93 with Sentinel-2; 66 with Landsat 8), due to their inferior temporal resolution, or would be possible with MODIS, due to its omission of small lakes < 0.125 km². Finally, we estimate the water volumes drained into the GrIS during rapid-lake-drainage events and, by analysing downscaled regional climate-model (RACMO2.3p2) run-off data, the water quantity that enters the GrIS via the moulins opened by such events. We find that during the lake-drainage events alone, the water drained by small lakes ( < 0.125 km²) is only 5.1 % of the total water volume drained by all lakes. However, considering the total water volume entering the GrIS after lake drainage, the moulins opened by small lakes deliver 61.5 % of the total water volume delivered via the moulins opened by large and small lakes; this is because there are more small lakes, allowing more moulins to open, and because small lakes are found at lower elevations than large lakes, where run-off is higher. These findings suggest that small lakes should be included in future remote-sensing and modelling work.

Uncertainty remains about how the surface hydrology of the Greenland ice sheet influences its subglacial drainage system, affecting basal water pressures and ice velocities, particularly over intraseasonal and interseasonal timescales. Here we apply a high spatial (200 m) and temporal (1 h) resolution subglacial hydrological model to a marginal (extending ~25 km inland), land‐terminating, ~200 km² domain in the Paakitsoq region, West Greenland. The model is based on that by Hewitt (2013) but adapted for use with both real topographic boundary conditions and calibrated modeled water inputs. The inputs consist of moulin hydrographs, calculated by a surface routing and lake‐filling/draining model, which is forced with distributed runoff from a surface energy‐balance model. Results suggest that the areal density of lake‐bottom moulins and their timing of opening during the melt season strongly affects subglacial drainage system development. A higher moulin density causes an earlier onset of subglacial channelization (i.e., water transport through channels rather than the distributed sheet), which becomes relatively widespread across the bed, whereas a lower moulin density results in a later onset of channelization that becomes less widespread across the bed. In turn, moulin density has a strong control on spatial and temporal variations in subglacial water pressures, which will influence basal sliding rates, and thus ice motion. The density of active surface‐to‐bed connections should be considered alongside surface melt intensity and extent in future predictions of the ice sheet's dynamics.

Glaciers in High Mountain Asia, many of which exhibit surface debris, contain the largest volume of ice outside of the polar regions. Many contain supraglacial pond networks that enhance melt rates locally, but no large‐scale assessment of their impact on melt rates exists. Here we use surface energy balance modeling forced using locally measured meteorological data and monthly satellite‐derived pond distributions to estimate the total melt enhancement for the four main glaciers within the 400‐km² Langtang catchment, Nepal, for a 6‐month period in 2014. Ponds account for 0.20 ± 0.03 m/year of surface melt, representing a local melt enhancement of a factor of 14 ± 3 compared with the debris‐covered area, and equivalent to 12.5 ± 2.0% of total catchment ice loss. Given the prevalence of supraglacial ponds across the region, our results suggest that effective incorporation of melt enhancement by ponds is essential for accurate predictions of future mass balance change in the region.

Surface runoff at the margin of the Greenland Ice Sheet (GrIS) drains to the ice-sheet bed, leading to enhanced summer ice flow. Ice velocities show a pattern of early summer acceleration followed by mid-summer deceleration due to evolution of the subglacial hydrology system in response to meltwater forcing. Modelling the integrated hydrological–ice dynamics system to reproduce measured velocities at the ice margin remains a key challenge for validating the present understanding of the system and constraining the impact of increasing surface runoff rates on dynamic ice mass loss from the GrIS. Here we show that a multi-component model incorporating supraglacial, subglacial, and ice dynamic components applied to a land-terminating catchment in western Greenland produces modelled velocities which are in reasonable agreement with those observed in GPS records for three melt seasons of varying melt intensities. This provides numerical support for the hypothesis that the subglacial system develops analogously to alpine glaciers and supports recent model formulations capturing the transition between distributed and channelized states. The model shows the growth of efficient conduit-based drainage up-glacier from the ice sheet margin, which develops more extensively, and further inland, as melt intensity increases. This suggests current trends of decadal-timescale slowdown of ice velocities in the ablation zone may continue in the near future. The model results also show a strong scaling between average summer velocities and melt season intensity, particularly in the upper ablation area. Assuming winter velocities are not impacted by channelization, our model suggests an upper bound of a 25 % increase in annual surface velocities as surface melt increases to 4× present levels.

Surface meltwater on ice shelves can exist as slush, it can pond in lakes or crevasses, or it can flow in surface streams and rivers. The collapse of the Larsen B Ice Shelf in 2002 has been attributed to the sudden drainage of ∼3000 surface lakes and has highlighted the potential for surface water to cause ice-shelf instability. Surface meltwater systems have been identified across numerous Antarctic ice shelves, although the extent to which these systems impact ice-shelf instability is poorly constrained. To better understand the role of surface meltwater systems on ice shelves, it is important to track their seasonal development, monitoring the fluctuations in surface water volume and the transfer of water across ice-shelf surfaces. Here, we use Landsat 8 and Sentinel-2 imagery to track surface meltwater across the Nivlisen Ice Shelf in the 2016–2017 melt season. We develop the Fully Automated Supraglacial-Water Tracking algorithm for Ice Shelves (FASTISh) and use it to identify and track the development of 1598 water bodies, which we classify as either circular or linear. The total volume of surface meltwater peaks on 26 January 2017 at 5.5×10⁷ m³. At this time, 63 % of the total volume is held within two linear surface meltwater systems, which are up to 27 km long, are orientated along the ice shelf's north–south axis, and follow the surface slope. Over the course of the melt season, they appear to migrate away from the grounding line, while growing in size and enveloping smaller water bodies. This suggests there is large-scale lateral water transfer through the surface meltwater system and the firn pack towards the ice-shelf front during the summer.