Ice-contact proglacial lakes affect ice dynamics and the transition of glacier termini from land- to lake-terminating has been shown to cause ice flow acceleration. In recent decades, the number and size of Greenlandic ice-marginal lakes has increased, highlighting the need to further understand these lake-terminating ice-margins as their influence on ice sheet mass balance increases. Here, time series of satellite-derived observations of ice velocity, surface elevation, and terminus position were generated at a lake-terminating outlet glacier, Isortuarsuup Sermia, and the nearby land-terminating Kangaasarsuup Sermia in south-west Greenland. At Isortuarsuup Sermia, annual surface velocity at the terminus increased by a factor of 2.5 to 214 ± 4 m yr−1 (2013–2021), with the magnitude of this acceleration declining with distance up-glacier. Meanwhile, near-terminus surface elevation changed at a rate of −2.3 ± 1.1 m yr−1 (2012–2021). Conversely, velocity change at Kangaasarsuup Sermia was minimal, while surface elevation change was approximately half at comparable elevations (−1.2 ± 0.3 m yr−1). We attribute these dynamic differences to thinning at Isortuarsuup Sermia and subsequent retreat from a stabilising sublacustrine moraine, and emphasise the potential of proglacial lakes to enhance future rates of mass loss from the Greenland Ice Sheet.

The Arctic is a hotspot for climate warming, making it crucial to quantify the sea level rise contribution from its ice masses. Novaya Zemlya's ice caps are the largest glacier complex in Europe and are a major contributor to contemporary sea level rise. Here we show that Novaya Zemlya outlet glaciers on the Barents Sea coast respond rapidly and consistently to oceanic forcing at annual timescales, likely due to their exposure to Atlantic Water variability. Glaciers on the Kara Sea show more variable response, likely reflecting their reduced exposure to Atlantic Water. Data demonstrate that the pause in glacier retreat previously observed on Novaya Zemlya between 2013 and 2015 has not persisted and that these changes correspond to ocean temperature variability on the Barents Sea coast. We document a marked shift to warmer air and ocean temperatures, and reduced sea ice concentrations from 2005 onwards. Although we identify ocean warming as the primary trigger for glacier retreat, we suggest that multi-year thinning, driven by the shift towards warmer air temperatures since 2005, pre-conditioned Novaya Zemlya's glaciers to retreat. Despite commonality in the timing of outlet glacier retreat, the magnitude is highly variable during rapid retreat phases, which we attribute to glacier-specific factors.

Supraglacial debris cover regulates the melt rates of many glaciers in mountainous regions around the world, thereby modifying the availability and quality of downstream water resources. However, the influence of supraglacial debris is often poorly represented within glaciological models, due to the absence of a technique to provide high-precision, spatially continuous measurements of debris thickness. Here, we use high-resolution UAV-derived thermal imagery, in conjunction with local meteorological data, visible UAV imagery and vertically profiled debris temperature time series, to model the spatially distributed debris thickness across a portion of Llaca Glacier in the Cordillera Blanca of Peru. Based on our results, we simulate daily sub-debris melt rates over a 3-month period during 2019. We demonstrate that, by effectively calibrating the radiometric thermal imagery and accounting for temporal and spatial variations in meteorological variables during UAV surveys, thermal UAV data can be used to more precisely represent the highly heterogeneous patterns of debris thickness and sub-debris melt on debris-covered glaciers. Additionally, our results indicate a mean sub-debris melt rate nearly three times greater than the mean melt rate simulated from satellite-derived debris thicknesses, emphasising the importance of acquiring further high-precision debris thickness data for the purposes of investigating glacier-scale melt processes, calibrating regional melt models and improving the accuracy of runoff predictions.

Greenland's future contribution to sea-level rise is strongly dependent on the extent to which dynamic perturbations, originating at the margin, can drive increased ice flow within the ice-sheet interior. However, reported observations of ice dynamical change at distances >~50 km from the margin have a very low spatial and temporal resolution. Consequently, the likely response of the ice-sheet's interior to future oceanic and atmospheric warming is poorly constrained. Through combining GPS and satellite-image-derived ice velocity measurements, we measure multi-decadal (1993–1997 to 2014–2018) velocity change at 45 inland sites, encompassing all regions of the ice sheet. We observe an almost ubiquitous acceleration inland of tidewater glaciers in west Greenland, consistent with acceleration and retreat at glacier termini, suggesting that terminus perturbations have propagated considerable distances (>100 km) inland. In contrast, outside of Kangerlussuaq, we observe no acceleration inland of tidewater glaciers in east Greenland despite terminus retreat and near-terminus acceleration, and suggest propagation may be limited by the influence of basal topography and ice geometry. This pattern of inland dynamical change indicates that Greenland's future contribution to sea-level will be spatially complex and will depend on the capacity for dynamic changes at individual outlet glacier termini to propagate inland.

Frontal ablation from tidewater glaciers is a major component of the total mass loss from the Greenland ice sheet. It remains unclear, however, how changes in atmospheric and oceanic temperatures translate into changes in frontal ablation, in part due to sparse observations at sufficiently high spatial and temporal resolution. We present high-frequency time-lapse imagery (photos every 30 min) of iceberg calving and meltwater plumes at Kangiata Nunaata Sermia (KNS), southwest Greenland, during June–October 2017, alongside satellite-derived ice velocities and modelled subglacial discharge. Early in the melt season, we infer a subglacial hydrological network with multiple outlets that would theoretically distribute discharge and enhance undercutting by submarine melt, an inference supported by our observations of terminus-wide calving during this period. During the melt season, we infer hydraulic evolution to a relatively more channelised subglacial drainage configuration, based on meltwater plume visibility indicating focused emergence of subglacial water; these observations coincide with a reduction in terminus-wide calving and transition to an incised planform terminus geometry. We suggest that temporal variations in subglacial discharge and near-terminus subglacial hydraulic efficiency exert considerable influence on calving and frontal ablation at KNS.

The increasingly negative mass balance of the Greenland ice sheet (GrIS) over the last ~25 years has been associated with enhanced surface melt and increased ice loss from marine-terminating outlet glaciers. Accelerated retreat during 2000–2010 was concentrated in the southeast and northwest sectors of the ice sheet; however, there was considerable spatial and temporal variability in the timing and magnitude of retreat both within and between these regions. This behaviour has yet to be quantified and compared for all glaciers in both regions. Furthermore, it is unclear whether retreat has continued after 2010 in the northwest, and whether the documented slowdown in the southeast post-2005 has been sustained. Here, we compare spatial and temporal patterns of frontal change in the northwest and southeast GrIS, for the period 2000–2015. Our results show near-ubiquitous retreat of outlet glaciers across both regions for the study period; however, the timing and magnitude of inter-annual frontal position change is largely asynchronous. We also find that since 2010, there is continued terminus retreat in the northwest, which contrasts with considerable inter-annual variability in the southeast. Analysis of the role of glacier-specific factors demonstrates that fjord and bed geometry are important controls on the timing and magnitude of glacier retreat.

The physical properties of eight ice cores recovered from along a flowline at Glacier de Tsanfleuron, Switzerland, have led to the identification of three distinctive internal zones. We use variations in the bulk ionic chemistry of these zones to approximate their relative liquid-water concentrations and ice viscosities. Results suggest that relative bulk water concentration and ice softness vary by over an order of magnitude between the zones. Implications of this variability for predictions of the glacier’s response to climate change are evaluated by incorporating these relative softnesses into a multi-layered (two-dimensional) model of ice flow. Model output is compared with that from an identical model constrained with a spatially uniform ice viscosity under advance and retreat modelling scenarios. The former scenario is used to tune viscosity by growing a glacier to its present long-section geometry, resulting in best-fit ice hardness values of 1.2 a1 bar–3 for the englacial ice in the multi-layered model and 7.0 a–1 bar–3 for all of the ice in the single-layered model. Both result in close approximations to the current long profile, yielding rms deviationsbetween measured and modelled ice thicknesses that are 5 5 m. In contrast, a single-layered model constrained with a hardness of 1.2 a–1 bar–3 overestimates the current measured long-section area by 31%, having a rms ice-thickness error of 15.0 m. Under the retreat modelling scenario, which gauges the response of the glacier to an imposed 75 m rise in equilibrium-line altitude (ELA), the multi-layered model predicts a long-section area reduction of 78%. This contrasts with a reduction of 64% for the single-layered model (hardness = 7.0 a–1 bar–3) and 85% for the single-layered model (hardness = 1.2 a–1 bar–3).

Records of surface motion, englacial tilt and repeat inclinometry are used to determine patterns of surface, internal and basal motion across the tongue of Haut Glacier d’Arolla, Switzerland, over temporal scales ranging from days to months. Findings are interpreted with reference to contemporaneous measurements of subglacial water pressures, and prior knowledge of the glacier’s subglacial drainage-system structure. Long-term inclinometry results show pronounced extrusion flow over a subglacial drainage axis, with basal velocities up to twice those measured at the glacier surface. Deformation profiles are more conventional away from the drainage axis, with basal velocities ∼60–70% of surface velocities. Comparison of long-term tilt rates from repeat inclinometry and englacial tiltmeters shows close correspondence. Englacial tiltmeter data are used to reconstruct internal velocity profiles and to split surface velocities into internal deformation and basal motion contributions over spring, summer and autumn/winter periods. Although, spatial patterns of surface movement are similar between periods, patterns of internal and basal motion are not. Results are interpreted in terms of the location of sticky and slippery spots, with temporally changing patterns of basal drag reflecting changing distributions of water pressure.

Channels incised into bedrock, or Nye channels, often form an important component of subglacial drainage at temperate glaciers, and their structure exerts control over patterns and rates of (a) channel erosion, (b) water flow-velocity and (c) water pressure. The latter, in turn, exerts a strong control over basal traction and, thus, ice dynamics. In order to investigate these controls, it is necessary to quantify detailed flow processes in subglacial Nye channels. However, it is effectively impossible to acquire such measurements from fully pressurized, subglacial channels. To solve this problem, we here apply a three-dimensional, finite-volume solution of the Reynolds averaged Navier– Stokes (RANS) equations with a one-equation mixing-length turbulence closure to simulate flow in a 3 m long section of an active Nye channel located in the immediate foreground of Glacier de Tsanfleuron, Switzerland. Numerical model output permits high-resolution visualization of water flow through the channel reach, and enables evaluation of the experimental manipulation of the pressure field adopted across the overlying ice lid. This yields an increased theoretical understanding of the hydraulic behaviour of Nye channels, and, in the future, of their effect on glacier drainage, geomorphology and ice dynamics.

Measurements of surface dynamics on polythermal John Evans Glacier, Nunavut, Canada, over two winter periods and every 7–10 days throughout two melt seasons (June–July 2000, 2001) provide new insight into spatio-temporal patterns of High Arctic glacier dynamics. In the lower ablation zone, mean annual surface velocities are 10–21 m a–1, but peak velocities up to 50% higher are attained during late June/early July. In the upper ablation zone and lower accumulation zone, mean annual surface velocities are typically 10–18 m a–1, and peak velocities up to 40% higher occur during late July. In the upper accumulation zone, mean annual surface velocities are 2–9 m a–1, and motion in mid- to late July exceeds this by up to 10%. Rapid drainage of ponded supraglacial water in the upper ablation zone to an initially distributed subglacial drainage system in mid-June may force excess surface motion in the warm-based lower glacier. The data indicate that the duration of the velocity response may be related to the rate of channelization of the basal drainage, and the velocity response may be transmitted up-glacier by longitudinal coupling. An increase in surface velocities in the middle glacier in late July occurs in conjunction with the opening of two further moulins in the accumulation zone.

The mass balance of polythermal ice masses is critically dependent on the proportion of surface-generated meltwater that subsequently refreezes in the snowpack and firn. In order to quantify this effect and to characterize its spatial variability, we measured near-surface (<10 m) snow and firn densities at an elevation of ~1945ma.s.l. in the percolation zone of the Greenland ice sheet in spring and autumn 2004. Results indicate that local snowpack depth above the previous end-of-summer 2003 melt surface increased by ±5% (7.6 cm) from spring to autumn while, over the same period, snowpack density increased by >26%, resulting in a 32% increase in net accumulation. This ‘seasonal densification’ increased at lower elevations, rising to 47% 10 km closer to the ice-sheet margin at 1860ma.s.l. Density/depth profiles from nine sites within 1 km2 at ~1945ma.s.l. reveal complex stratigraphies that change over short spatial scales and seasonally. We conclude that estimates of mass-balance change cannot be calculated solely from observed changes in surface elevation, but that near-surface densification must also be considered. However, predicting spatial and temporal variations in densification may not be straightforward. Further, the development of complex firn-density profiles both masks discernible annual layers in the near-surface firn and ice stratigraphy and is likely to introduce error into radar-derived estimates of surface elevation.

An indirect methodology for determining the distribution of mass balance at high spatial resolution using remote sensing and ice-flow modelling is presented. The method, based on the mass-continuity equation, requires two datasets collected over the desired monitoring interval: (i) the spatial pattern of glacier surface-elevation change, and (ii) the mass-flux divergence field. At Haut Glacier d’Arolla, Valais, Switzerland, the mass-balance distribution between September 1992 and September 1993 is calculated at 20 m resolution from the difference between the pattern of surface-elevation change derived from analytical photogrammetry and the mass-flux divergence field determined from three-dimensional, numerical flow modelling constrained by surface-velocity measurements. The resultant pattern of mass balance is almost totally negative, showing a strong dependence on elevation, but with large localized departures. The computed distribution of mass balance compares well (R2 = 0.91) with mass-balance measurements made at stakes installed along the glacier centre line over the same period. Despite the highly optimized nature of the flow-modelling effort employed in this study, the good agreement indicates the potential this method has as a strategy for deriving high spatial and temporal-resolution estimates of mass balance.

Interferometric synthetic aperture radar (InSAR) observations of ice-shelf flow contain ocean-tide and atmospheric-pressure signals. A model-based correction can be applied, but this method is limited by its dependency upon model accuracy, which in remote regions can be uncertain. Here we describe a method to determine two-dimensional ice-shelf flow vectors independently of model predictions of tide and atmospheric pressure, by stacking conventional and multiple aperture InSAR (MAI) observations of the Dotson Ice Shelf, West Antarctica. In this way we synthesize a longer observation period, which enhances long-period (flow) displacement signals, relative to rapidly varying (tide and atmospheric pressure) signals and noise. We estimate the error associated with each component of the velocity field to be ~22 ma-1, which could be further reduced if more images were available to stack. With the upcoming launch of several satellite missions, offering the prospect of regular short-repeat SAR acquisitions, this study demonstrates that stacking can improve estimates of ice-shelf flow velocity.

The longitudinal pattern of surface velocity of a large, predominantly cold, polythermal glacier (John Evans Glacier, Ellesmere Island, Canada) was measured over summer and winter periods. In the accumulation and upper ablation areas, where ice is predominantly cold-based, summer velocities were slightly higher than overwinter velocities. Predicted velocities due to ice deformation alone in these areas closely matched these observations in the winter, with limited basal motion likely in the summer. In the lower ablation area, where ice is likely warm-based, measured summer velocities were up to double overwinter velocities. Predicted ice deformation could not account for all of these measured velocities in either summer or winter. This suggests that basal motion occurs throughout the year over at least part of the lower ablation area. This finding is supported by radio-echo sounding, subglacial drainage reconstructions and analyses of early-summer meltwater chemistry, which suggest that subglacial water is present throughout the year in this region. In summer, basal motion may account for up to 75% of the total surface velocity throughout the lower ablation area. The inferred rate of basal motion increases sharply directly below a set of moulins by which most surface meltwater reaches the glacier bed.

Three early-melt-season high-velocity events (or “spring events”) occurred on Haut Glacier d’Arolla, Switzerland, during the melt seasons of 1998 and 1999. The events involve enhanced glacier velocity during periods of rapidly increasing bulk discharge in the proglacial stream and high subglacial water pressures. However, differences in spatial patterns of surface velocity, internal ice deformation rates, the spatial extent of high subglacial water pressures and in rates of subglacial sediment deformation suggest different hydrological and mechanical controls. The data from two of the events suggest widespread ice–bed decoupling, particularly along a subglacial drainage axis creating the highest rates of basal motion and “plug flow” in the overlying ice. The other event showed evidence of less extensive ice–bed decoupling and sliding along the drainage axis with more mechanical support for ice overburden transferred to areas adjacent to decoupled areas. We suggest that: (1) plug flow may be a common feature on glaciers experiencing locally induced reductions in basal drag; (2) under certain circumstances, enhanced surface motion may be due in part to non-locally forced enhanced bed deformation; and (3) subglacial sediment deformation is confined to a depth of the order of centimetres to decimetres.

The surface motion of Haut Glacier d’Arolla, Switzerland, was monitored at a high spatial and temporal resolution. Data are analyzed to calculate surface velocities, surface strain rates and the components of the glacier force budget before, during and after an early melt season speed-up or “spring event”. We investigate the extent to which variations in glacier motion can be attributed to hydrologically induced local forcing or to non-local forcing transmitted via horizontal stress gradients. Enhanced glacier motion is dependent on a change in the spatial distribution of areas of high drag across the glacier.

The surface velocity of a predominantly cold polythermal glacier (John Evans Glacier, Ellesmere Island, Canada) varies significantly on both seasonal and shorter time-scales. Seasonal variations reflect the penetration of supraglacial water to the glacier bed through significant thicknesses of cold ice. Shorter-term events are associated with periods of rapidly increasing water inputs to the subglacial drainage system. Early-season short-term events immediately follow the establishment of a drainage connection between glacier surface and glacier bed, and coincide with the onset of subglacial outflow at the terminus. A mid-season short-term event occurred as surface melting resumed following cold weather, and may have been facilitated by partial closure of subglacial channels during this cold period. There is a close association between the timing and spatial distribution of horizontal and vertical velocity anomalies, the temporal pattern of surface water input to the glacier, and the formation, seasonal evolution and distribution of subglacial drainage pathways. These factors presumably control the occurrence of highwater-pressure events and water storage at the glacier bed. The observed coupling between surface water inputs and glacier velocity may allow predominantly cold polythermal glaciers to respond rapidly to climate-induced changes in surface melting.

Dye-tracer experiments undertaken over two summer melt seasons at polythermal John Evans Glacier, Ellesmere Island, Canada, were designed to investigate the character of the subglacial drainage system and its evolution over a melt season. In both summers, dye injections were conducted at several moulins and traced to a single subglacial outflow. Tracer breakthrough curves suggest that supraglacial meltwater initially encounters a distributed subglacial drainage system in late June. The subsequent development and maintenance of a channelled subglacial network are dependent upon sustained high rates of surface melting maintaining high supraglacial inputs. In a consistently warm summer (2000), subglacial drainage became rapidly and persistently channelled. In a cooler summer (2001), distributed subglacial drainage predominated. These observations confirm that supraglacial meltwater can access the bed of a High Arctic glacier in summer, and induce significant structural evolution of the subglacial drainage system. They do not support the view that subglacial drainage systems beneath polythermal glaciers are always poorly developed. They do suggest that the effects on ice flow of surface water penetration to the bed of predominantly cold glaciers may be short-lived.

The transfer of surface-generated meltwater to the subglacial drainage system through full ice thickness crevassing may lead to accelerated glacier velocities, with implications for ice motion under future climatic scenarios. Accurate predictions of where surface meltwater accesses the ice/bed interface are therefore needed in fully coupled hydrodynamic ice-sheet models. We present a spatially distributed modelling routine for predicting the location and timing of delivery of surface-derived meltwater to the ice/bed interface through moulins and supraglacial lake drainage. The model is explained as it is applied to the Croker Bay glacial catchment of Devon Ice Cap, Canada. The formation of moulins, drainage of lakes, and the transfer of meltwater through the full ice thickness are modelled for the 2004 and 2006 ablation seasons. Through this case study we assess the model’s sensitivity to degree-day factors, fracture toughness, tensile strength and crevasse width, and confirm that parameters influencing the rate at which water fills a crevasse are the most significant controls on the ability of a crevasse to reach the bed. Increased surface melt production, therefore, has the potential to significantly influence the spatial and temporal transfer of meltwater through surface-to-bed connections in a warmer climate

This paper describes the physical and chemical properties of the snow- pack on John Evans Glacier, Ellesmere Island, Canadian Arctic Archipelago, and investigates the controls on snowpack solute concentrations and atmospheric deposition. The snowpack contains three layers that are traceable across the whole glacier. These represent fall accumulation that has been metamorphosed to depth hoar, winter accumulation mixed with snow reworked by wind from the underlying depth hoar, and spring accumulation mixed with wind-reworked snow. The seasonal cycle in snow chemistry closely reflects changes in the composition of the atmospheric aerosol at Alert, with some modification of NO3− concentrations by post-depositional processes. Mean water-weighted solute concentrations in the snowpack are largely independent of accumulation, while atmospheric deposition tends to increase with accumulation. This suggests that, for most species, wet deposition is the dominant depositional process throughout the year. However, concentrations of Ca2+ and K+ increase with both accumulation and elevation, implying an enhanced input from dry deposition of soil dust above 800 m elevation. Concentrations of SO42− are inversely related to accumulation, especially in the winter layer, suggesting a significant input from non-precipitating events, such as dry deposition or riming, during this period of very limited snowfall.