We present in situ firn temperatures from the extreme 2012 melt season in the southwestern lower accumulation area of the Greenland ice sheet. The upper 2.5 m of snow and firn was temperate during the melt season, when vertical meltwater percolation was inefficient due to a ~5.5 m thick ice layer underlying the temperate firn. Meltwater percolation and refreezing beneath 2.5 m depth only occurred after the melt season. Deviations from temperatures predicted by pure conductivity suggest that meltwater refroze in discrete bands at depths of 2.0–2.5, 5.0–6.0 and 8.0–9.0 m. While we find no indication of meltwater percolation below 9 m depth or complete filling of pore volume above, firn at 10 and 15 m depth was respectively 4.2–4.5°C and 1.7°C higher than in a conductivity-only simulation. Even though meltwater percolation in 2012 was inefficient, firn between 2 and 15 m depth the following winter was on average 4.7°C warmer due to meltwater refreezing. Our observations also suggest that the 2012 firn conditions were preconditioned by two warm summers and ice layer formation in 2010 and 2011. Overall, firn temperatures during the years 2009–13 increased by 0.6°C.

The retention and release of liquid water in glacierized basins was modelled with a conceptual, semi-distributed model of the water and ice balance designed for long-term averages with monthly resolution for 100 m elevation bands. Here we present the components of the liquid water balance of 86 mostly glacierized basins on either side of the main Alpine divide between 10 and 13°E in the period 1998–2006 and compare them with the records of 30 basins monitored from 1970 to 1997. Basin average of liquid water retention has maxima in excess of 100 mm per month in May, often followed by maximum release when the retaining snow matrix melts. Glacier storage peaks in August partly due to ice melt and the ensuing filling of the englacial reservoirs and partly on account of a precipitation maximum. These two components combined to a common maximum of storage in summer in the first period 1970–97 and developed two distinct maxima in the warmer period 1998–2006. A further maximum of liquid water storage that was often found in October is most likely due to a peak in precipitation in the southern part of the study region.

Glacier naledi are extrusive ice masses that appear in front of glaciers as a consequence of refreezing of meltwater seepage during the accumulation season. These structures provide a unique opportunity to understand subglacial drainage activity during the accumulation season; however, only few detailed studies have previously focused on their characteristics. Here, we investigated glacier-derived naled assemblages in the proglacial zone of the polythermal glacier Werenskioldbreen (27.4 km2) in SW Svalbard. We determined the spatial distribution of naledi using ground penetrating radar surveys. The main subglacial drainage pattern was related to a channel under the medial moraine, and three sources are linked to a distributed subglacial drainage network. The relation between atmospherically-corrected (Ca2+ + Mg2+) and (SO4 2−) in sub-naled waters was closely related to sulphide oxidation coupled with carbonate dissolution (r = 0.99; slope = 1.6). This is consistent with the local lithology, which is dominated by schist containing carbonates. We also found high carbonate saturation indices in pale white ice layers within the naled. We conclude that sulphide oxidation coupled with carbonate dissolution is the dominant chemical weathering process in the subglacial drainage system of Werenskioldbreen during the accumulation season.

Supraglacial streams are a significant part of the glacial hydrological system and important for understanding the connection between glacial hydrology and glacier dynamics. Here we determine the factors that influence the development of step-pool formation and pulsating flow in a supraglacial stream on Bylot Island, Nunavut. Results show that during the second week of a 2-week study, multiple successive rainfall events occurred, stream temperature increased and ablation decreased; which also caused stream discharge to decrease. In addition, the stream, which flowed over a 13 m high waterfall off the front of Fountain Glacier, rapidly formed 21 step-pools and began to pulsate. The pulsating phenomenon involved the complete stoppage of flow over the waterfall and the subsequent restart between 8 and 20 s later. Pulsating flow resulted from rapid changes in the streambed morphology. In particular, the formation of the step-pool sequence was caused by helical flow around meander bends and hydrologically induced slippage along transverse shear planes, evidenced by observations of high-pressure artesian flow from transverse fractures. Contrary to previous literature, this study shows that high discharge is not necessarily the cause of step-pool formation and pulsating flow within supraglacial streams.

Thermodynamics of a seasonal supraglacial lake were examined based on field data from three summers. At maximum, the lake body consisted of an upper layer with thin ice on top, and a lower layer with slush, hard ice and sediment at the bottom. Sublimation from the upper ice surface averaged to 0.7 mm d−1, and melting in the ice interior averaged to 9.1 mm d−1 during summer. Albedo was on average 0.6 and light attenuation coefficient was ~1 m−1. Averaged over December and January, and over 3 different years, we found that the net solar heating was 137 W m−2, while the losses averaged to 62 W m−2 for the longwave radiation, 16 W m−2 for the sensible heat flux, 24 W m−2 for the latent heat flux and 3 W m−2 for the bottom flux. The depth scale is determined by the light attenuation distance and thermal diffusion coefficient, and the net liquid water volume ranged from 0.5 to 1.0 m in different years. The potential winter growth is more than summer melting, and thus the lake freezes up completely in winter in the present climate.

GPS campaigns on glaciers during jökulhlaups show how subglacial floods affect glacier motion and shed light on the dynamics of such floods. Three such campaigns have been carried out on southern and western Vatnajökull, southeast Iceland, over known jökulhlaup paths. Two slowly rising jökulhlaups from Grímsvötn and two rapidly rising jökulhlaups from the western and eastern Skaftá cauldrons were captured in these campaigns, with maximum discharge ranging from 240 to 3300 m3 s−1. Glacier surface movements measured in these campaigns are presented along with the corresponding discharge curves. The measurements are interpreted as indicating: (1) initiation of rapidly rising jökulhlaups with a propagating subglacial pressure wave, (2) decreased glacier basal friction during jökulhlaups, (3) subglacial accumulation of water in slowly rising jökulhlaups and (4) lifting of the glacier caused by subglacial water pressure exceeding overburden in both rapidly and slowly rising jökulhlaups. The latter two observations are inconsistent with assumptions that are typically made in theoretical and numerical modelling of jökulhlaups. Both viscous and elastic deformation of the glacier as well as turbulent hydraulic fracture at the ice/bedrock interface are important in the dynamics of the subglacial pressure wave at the front of rapidly rising jökulhlaups.

A suite of surface and basal measurements during and after borehole drilling is used to perform in situ investigation of the local basal drainage system and pressure forcing in western Greenland. Drill and borehole water temperature were monitored during borehole drilling, which was performed with dyed hot water. After drilling, borehole water pressure and basal dye concentration were measured concurrently with positions in a GPS strain diamond at the surface. Water pressure exhibited diurnal changes in antiphase with velocity. Dye monitoring in the borehole revealed stagnant basal water for nearly 2 weeks. The interpretation of initial connection to an isolated basal cavity is corroborated by the thermal signature of borehole water during hot water drilling. Measurement-based estimates of cavity size are on the order of cubic meters, and analysis indicates that small changes in its volume could induce the observed pressure variations. It is found that longitudinal coupling effects are unable to force necessary volume changes at the site. Sliding-driven basal cavity opening and elastic uplift from load transfer are plausible mechanisms controlling pressure variations. Elastic uplift requires forcing from a hydraulically connected reach, which observations suggest must be relatively small and in close proximity to the isolated cavity.

We use daily surface velocities measured over several weeks in 2007 and 2008 on a slowly surging glacier in Yukon, Canada, to examine the ordinary melt-season dynamics in the context of the ongoing surge. Horizontal velocities within and just below the ~1.5 km-long zone of fastest flow, where the surge is occurring, are often correlated during intervals of low melt. This correlation breaks down during melt events, with the lower reaches of the fast-flow zone responding first. Velocity variability in this lower reach is most highly correlated with melt; velocities above and below appear to respond at least as strongly to the velocity variations of this reach as to local melt. GPS height records are suggestive of ice/bed separation occurring in the fast-flow zone but not below it, pointing to a hydrological cause for the short-term flow variability in the surging region. Independent velocity measurements over 6 years show a maximum July flow anomaly coincident with the location most responsive to melt. Results from a simple model of dashpots and frictional elements lend support to the hypothesis that this zone partly drives the dynamics of the ice above and below it. We speculate that the slow surge may enhance glacier sensitivity to melt-season processes, including short-term summer sliding events.

The subglacial hydrologic system exerts strong controls on the dynamics of the overlying ice, yet the parameters that govern the evolution of this system are not widely known or observable. To gain a better understanding of these parameters, we invert a spatially averaged model of subglacial hydrology from observations of ice surface velocity and outlet stream discharge at Kennicott Glacier, Wrangell Mountains, AK, USA. To identify independent parameters, we formally non-dimensionalize the forward model. After specifying suitable prior distributions, we use a Markov-chain Monte Carlo algorithm to sample from the distribution of parameter values conditioned on the available data. This procedure gives us not only the most probable parameter values, but also a rigorous estimate of their covariance structure. We find that the opening of cavities due to sliding over basal topography and turbulent melting are of a similar magnitude during periods of large input flux, though turbulent melting also exhibits the greatest uncertainty. We also find that both the storage of water in the englacial system and the exchange of water between englacial and subglacial systems are necessary in order to explain both surface velocity observations and the relative attenuation in the amplitude of diurnal signals between input and output flux observations.

Antarctica Bedmap2 datasets are used to calculate subglacial hydraulic potential and the area, depth and volume of hydraulic potential sinks. There are over 32 000 contiguous sinks, which can be thought of as predicted lakes. Patterns of subglacial melt are modelled with a balanced ice flux flow model, and water fluxes are cumulated along predicted flow pathways to quantify steady-state fluxes from the main basin outlets and from known subglacial lakes. The total flux from the continent is ~21 km3 a−1. Byrd Glacier has the greatest basin flux of ~2.7 km3 a−1. Fluxes from subglacial lakes range from ~1 × 10−4 to ~1.5 km3 a−1. Lake turnover times are calculated from their volumes and fluxes, and have median values of ~100 a for known ‘active’ lakes and ~500 a for other lakes. Recurrence intervals of a 0.25 km3 flood range from ~2 months to ~2000 a (median ≈130 a) for known ‘active’ lakes and from ~2 to ~2400 a (median ≈ 360 a) for other lakes. Thus, several lakes that have recently been observed to fill and drain may not do so again for many centuries; and several lakes that have not, so far, been observed to fill and drain have the potential to do so, even at annual to decadal timescales.

Lakes beneath the Antarctic Ice Sheet are known to decrease traction at the ice base and therefore can have a great impact on ice dynamics. However, the total extent of Antarctic subglacial lakes is still unknown. We address this issue by combining modeling and remote-sensing strategies to predict potential lake locations using the general hydraulic potential equation. We are able to reproduce the majority of known lakes, as well as predict the existence of many new and so far undetected potential lakes. To validate our predictions, we analyzed ice-penetrating radar profiles from radio-echo sounding flights acquired over 1994–2013 in Dronning Maud Land, East Antarctica, and this led to the identification of 31 new subglacial lakes. Based on these findings, we estimate the total number of Antarctic subglacial lakes to be ~1300, a factor of three higher than the total number of lakes discovered to date. We estimate that only ~30% of all Antarctic subglacial lakes and ~65% of the total estimated lake-covered area have been discovered, and that lakes account for 0.6% of the Antarctic ice/bed interface.

Marine-terminating outlet glaciers discharge most of the Greenland ice sheet's mass through frontal ablation and meltwater runoff. While calving can be estimated by in situ and remote sensing observations, submarine melting and subglacial meltwater transport are more challenging to quantify. Here we investigate the subglacial hydrology of Rink Isbræ, a fast-flowing West Greenland tidewater glacier, using time-lapse photography, modeled runoff estimates and daily satellite imagery from 2007 to 2011. We find that sediment plumes appear episodically at four distinct locations across the terminus, and last between 2 h and 17 d. This suggests short-term variability in discharge and the existence of persistent pathways. The seasonal onset of sediment plumes occurs before supraglacial lake drainages, shortly after the onset of runoff, and only after the wintertime ice mélange has begun disintegrating. Plumes were also visible after the cessation of runoff (23 ± 5 d), which is indicative of subglacial storage. The lack of either a seasonal velocity change or a correspondence between meltwater availability and plume occurrence suggests that the subglacial system persists in a state of inefficient drainage. Subglacial hydrology at tidewater glaciers is of critical importance in understanding dynamics at the ice front.

Fláajökull is a non-surging outlet glacier draining the south-eastern part of the Vatnajökull, southeast Iceland. Fláajökull was stationary or advanced slightly between 1966 and 1995 and formed a prominent end moraine. Glacial retreat since then has revealed a cluster of 15 drumlins. This study focuses on the morphology and sedimentology of the drumlins. They are 100–600 m long, 40–130 m wide, and have cores of glaciofluvial sediment or till. The drumlins are draped by ~1 m thick, massive subglacial traction till. The glacier forefield is characterized by a number of arcuate and saw-tooth, terminal and recessional moraine ridges, overridden moraines with fluted surfaces, and glaciofluvial outwash. Some of the drumlins extend towards the 1995 end moraine but terminate abruptly at the moraine and are not observed in front of it. This suggests that they were formed sub-marginally during the 1966–1995 terminal position. The sedimentary structure of the drumlins is best explained by the sticky spot model. Dating and dendrochronological analyses of birch logs found on the surface of one of the drumlins indicate that the valley was forested about 2100 calendar year BP, after which the glacier started to reform, possibly due to an abrupt change in climate.

The hypothesized link between glacier surging and bedrock geology motivates this study of the suspended sediment size distributions (SSSD) from surge-type and non-surge-type glaciers. We analyze SSSDs from proglacial streams in 20 individual basins comprising various fractions of metasedimentary (MS) and felsic plutonic rocks. We compare the size distributions by performing tests of significance on the distribution statistics, and a principal component analysis on discrete grain sizes. We find that surge-type and non-surge-type glaciers underlain solely by MS rocks have significantly different SSSDs, while surge-type glaciers as a whole have remarkably similar SSSDs, regardless of the underlying bedrock geology. These observations hint at a relationship between sediment characteristics and glacier surging, though causation in either direction cannot be established without additional data.

Virkisjökull is a rapidly retreating outlet glacier draining the western flanks of Öræfajökull in SE Iceland. Since 2011 there have been continuous measurements of flow in the proglacial meltwater channel and regular campaigns to sample stable isotopes δ2H and δ18O from the river, ice, moraine springs and groundwater. The stable isotopes provide reliable end members for glacial ice and shallow groundwater. Analysis of data from 2011 to 2014 indicates that although ice and snowmelt dominate summer riverflow (mean 5.3–7.9 m3 s−1), significant flow is also observed in winter (mean 1.6–2.4 m3 s−1) due primarily to ongoing glacier icemelt. The stable isotope data demonstrate that the influence of groundwater discharge from moraines and the sandur aquifer increases during winter and forms a small (15–20%) consistent source of baseflow to the river. The similarity of hydrological response across seasons reflects a highly efficient glacial drainage system, which makes use of a series of permanent englacial channels within active and buried ice throughout the year. The study has shown that the development of an efficient year round drainage network within the lower part of the glacier has been coincident with the stagnation and subsequent rapid retreat of the glacier.