Abstract. We present detailed records of lake discharge, ice motion and passive seismicity capturing the behaviour and processes preceding, during and following the rapid drainage of a ∼ 4 km 2 supraglacial lake through 1.1-km-thick ice on the western margin of the Greenland Ice Sheet. Peak discharge of 3300 m 3 s −1 coincident with maximal rates of vertical uplift indicates that surface water accessed the ice-bed interface causing widespread hydraulic separation and enhanced basal motion. The differential motion of four global positioning system (GPS) receivers located around the lake record the opening and closure of the fractures through which the lake drained. We hypothesise that the majority of discharge occurred through a ∼ 3-km-long fracture with a peak width averaged across its wetted length of ∼ 0.4 m. We argue that the fracture's kilometre-scale length allowed rapid discharge to be achieved by combining reasonable water velocities with sub-metre fracture widths. These observations add to the currently limited knowledge of in situ supraglacial lake drainage events, which rapidly deliver large volumes of water to the ice-bed interface.

Abstract. We present a multi-temporal digital inventory of Svalbard glaciers with the most recent from the late 2000s containing 33 775 km 2 of glaciers covering 57 % of the total land area of the archipelago. At present, 68 % of the glacierized area of Svalbard drains through tidewater glaciers that have a total terminus width of ∼ 740 km. The glacierized area over the entire archipelago has decreased by an average of 80 km 2 a −1 over the past ∼ 30 yr, representing a reduction of 7 %. For a sample of ∼ 400 glaciers (10 000 km 2 ) in the south and west of Spitsbergen, three digital inventories are available from the 1930/60s, 1990 and 2007 from which we calculate average changes during 2 epochs. In the more recent epoch, the terminus retreat was larger than in the earlier epoch, while area shrinkage was smaller. The contrasting pattern may be explained by the decreased lateral wastage of the glacier tongues. Retreat rates for individual glaciers show a mix of accelerating and decelerating trends, reflecting the large spatial variability of glacier types and climatic/dynamic response times in Svalbard. Lastly, retreat rates estimated by dividing glacier area changes by the tongue width are larger than centerline retreat due to a more encompassing frontal change estimate with inclusion of lateral area loss.

Meltwater drainage across the surface of the Greenland ice sheet (GrIS) is well constrained by measurements and modeling, yet despite its critical role, knowledge of its transit through the subglacial environment remains limited. Here we present a subglacial hydrological analysis of a land-terminating sector of the GrIS at unprecedented resolution that predicts the routing of surface-derived meltwater once it has entered the basal drainage system. Our analysis indicates the probable existence of small subglacial lakes that remain undetectable by methods using surface elevation change or radar techniques. Furthermore, the analysis suggests transient behavior with rapid switching of subglacial drainage between competing catchments driven by seasonal changes in the basal water pressure. Our findings provide a cautionary note that should be considered in studies that attempt to relate and infer future response from surface temperature, melt, and runoff from point measurements and/or modeling with measurements of proglacial discharge and ice dynamics. Hydraulic potential analysis enables subglacial drainage to be estimated from ice thickness and basal topography on spatial scales where the main control of the drainage is the geometry of the ice sheet [Shreve, 1972]. This method has been applied in Greenland [e.g.

Abstract. The climate in Svalbard is undergoing amplified change compared to the global mean. This has major implications for runoff from glaciers and seasonal snow on land. We use a coupled energy balance–subsurface model, forced with downscaled regional climate model fields, and apply it to both glacier-covered and land areas in Svalbard. This generates a long-term (1957–2018) distributed dataset of climatic mass balance (CMB) for the glaciers, snow conditions, and runoff with a 1 km×1 km spatial and 3-hourly temporal resolution. Observational data including stake measurements, automatic weather station data, and subsurface data across Svalbard are used for model calibration and validation. We find a weakly positive mean net CMB (+0.09 m w.e. a−1) over the simulation period, which only fractionally compensates for mass loss through calving. Pronounced warming and a small precipitation increase lead to a spatial-mean negative net CMB trend (−0.06 m w.e. a−1 decade−1), and an increase in the equilibrium line altitude (ELA) by 17 m decade−1, with the largest changes in southern and central Svalbard. The retreating ELA in turn causes firn air volume to decrease by 4 % decade−1, which in combination with winter warming induces a substantial reduction of refreezing in both glacier-covered and land areas (average −4 % decade−1). A combination of increased melt and reduced refreezing causes glacier runoff (average 34.3 Gt a−1) to double over the simulation period, while discharge from land (average 10.6 Gt a−1) remains nearly unchanged. As a result, the relative contribution of land runoff to total runoff drops from 30 % to 20 % during 1957–2018. Seasonal snow on land and in glacier ablation zones is found to arrive later in autumn (+1.4 d decade−1), while no significant changes occurred on the date of snow disappearance in spring–summer. Altogether, the output of the simulation provides an extensive dataset that may be of use in a wide range of applications ranging from runoff modelling to ecosystem studies.

We present a multi-temporal digital inventory of Svalbard glaciers with the most recent from the late 2000s containing 33 775 km² of glaciers, or 57% of the total land area of the archipelago. At present, 68% of the glaciated area of Svalbard drains through tidewater glaciers that have a summed terminus width of ~ 740 km. The glaciated area over the entire archipelago has decreased by an average of 80 km² a⁻¹ over the past ~ 30 yr, representing a reduction of 7%. For a sample of ~ 400 glaciers (10 000 km²) in the south and west of Spitsbergen, three digital inventories are available from 1930/60s, 1990 and 2007 from which we calculate average changes during 2 epochs. In the more recent epoch, the terminus retreat was larger than in the earlier epoch while area shrinkage was smaller. The contrasting pattern may be explained by the decreased lateral wastage of the glacier tongues. Temporal retreat rates for individual glaciers show a mix of accelerating and decelerating trends, reflecting the large spatial variability of glacier types and climatic/dynamic response times in Svalbard. Last, retreat rates estimated by dividing glacier area changes by the tongue width are larger than centerline retreat due to a more encompassing frontal change estimate with inclusion of lateral area loss