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. Supraglacial lakes represent an ephemeral storage buffer for meltwater runoff and lead to significant, yet shortlived, episodes of ice-flow acceleration by decanting large meltwater and energy fluxes into the ice sheet's hydrological system. Here, a methodology for calculating lake volume is used to quantify storage and drainage across Russell Glacier, West Greenland, between 2002 and 2012. Using 502 MODIS scenes, water volume at ∼ 200 seasonally occurring lakes was derived using a depth-reflectance relationship, which was independently calibrated and field validated against lake bathymetry. The inland expansion of lakes is strongly correlated with air temperature: during the record melt years of 2010 and 2012, lakes formed and drained earlier, attaining their maximum volume 38 and 20 days earlier than the 11 yr mean, as well as occupying a greater area and forming at higher elevations (> 1800 m) than previously. Despite occupying under 2 % of the study area, lakes delay the transmission of up to 7-13 % of the bulk meltwater discharged. Although the results are subject to an observational bias caused by periods of cloud cover, we estimate that across Russell Glacier, 28 % of supraglacial lakes drain rapidly (< 4 days). Clustering of such events in space and time suggests a synoptic trigger mechanism. Further, we find no evidence to support a unifying critical size or depth-dependent drainage threshold.

Subglacial hydrology plays a key role in many glaciological processes, including ice dynamics via the modulation of basal sliding. Owing to the lack of an overarching theory, however, a variety of model approximations exist to represent the subglacial drainage system. The Subglacial Hydrology Model Intercomparison Project (SHMIP) provides a set of synthetic experiments to compare existing and future models. We present the results from 13 participating models with a focus on effective pressure and discharge. For many applications (e.g. steady states and annual variations, low input scenarios) a simple model, such as an inefficient-system-only model, a flowline or lumped model, or a porous-layer model provides results comparable to those of more complex models. However, when studying short term (e.g. diurnal) variations of the water pressure, the use of a two-dimensional model incorporating physical representations of both efficient and inefficient drainage systems yields results that are significantly different from those of simpler models and should be preferentially applied. The results also emphasise the role of water storage in the response of water pressure to transient recharge. Finally, we find that the localisation of moulins has a limited impact except in regions of sparse moulin density.

Amplitude-versus-angle (AVA) analysis of a seismic reflection line, imaged 13km from Russell Glacier terminus, near the western margin of the Greenland ice sheet (GrIS), suggests the presence of sediment at the bed. The analysis was complicated by the lack of identifiable multiples in the data due to a highly irregular and crevassed ice surface, rendering deeper seismic returns noisy. A modified technique for AVA processing of glacial seismic data using forward modelling with primary reflection amplitudes and simulated multiple amplitudes is presented here. Our analysis demonstrates that AVA analysis can be applied to areas with noisy seismic returns and indicates that sediment underlies the seismic study site. Our data are inconsistent with the common assumption that the GrIS is underlain only by hard bedrock, but consistent with the presence of subglacial sediment with porosity between 30% and 40%. As analysis and modelling of ice-sheet dynamics requires a sound knowledge of the underlying basal materials, subglacial sediment should be taken into account when considering ice dynamics in this region of the GrIS.Peer reviewe

Ice shelf basal channels cause transverse fractures that can be exacerbated by surface rivers, culminating in calving.