Ice-contact lakes modify glacier geometry and dynamics by shifting the majority of mass loss from the ice surface to the terminus. Lake-terminating glaciers are known to experience greater thinning rates and higher velocities than land-terminating glaciers, but the controls on variability in surface elevation change and ice flow between lake-terminating glaciers in different regions remain poorly explored. We combined existing datasets of glacier velocity, surface elevation change and glacial lake area to characterise the evolution of 352 lake-terminating and land-terminating glaciers within three Himalayan sub-regions between 2000 and 2019. These analyses show that the influence of ice-contact lakes propagates up-glacier across only the lowermost 30% of the hypsometric distribution, even where lakes are well established. We find that ice-contact lakes only affect glacier behaviour when the lakes reach an advanced evolutionary stage; most clearly manifested in the Eastern Himalaya by statistically robust differences in glacier-wide surface elevation change between lake-terminating (–0.68 ± 0.05 m a–1) and land-terminating (–0.54 ± 0.04 m a–1) glaciers. These differences are driven by the presence of a greater number of well-developed ice-contact lakes in the Eastern Himalaya compared to in the Western and Central Himalaya, resulting from greater mass loss rates to date.

The debris that covers the ablation areas of high-elevation debris-covered glaciers contributes to the distinctive features and processes occurring both on and within such glaciers. Despite recent advances, knowledge of the subsurface environments of high-elevation debris-covered glaciers is still extremely limited. In particular, targeted field-based data are needed to parameterise and refine the projections of these glaciers in numerical models. Here, we outline the current understanding of the internal properties of high-elevation debris-covered glaciers based on direct field-based methods and suggest future research directions for field-based studies.

Characterising the structures within glaciers can give unique insight into ice motion processes. On debris-covered glaciers, traditional structural glaciological mapping is challenging because the lower glacier is hidden by the supraglacial debris layer. Here, we use high-resolution optical televiewer (OPTV) image logs from four boreholes drilled into Khumbu Glacier, Nepal, to overcome this limitation and investigate englacial structural features within a Himalayan debris-covered glacier. The OPTV logs show structural features that are up to an order of magnitude thinner than those observed at the glacier surface and reveal five structural units: (I) primary stratification of ice; (II) debris-rich planes that conform with the primary stratification; (III) water-healed crevasse traces; (IV) healed crevasse traces; and (V) steeply dipping planes of basally derived fine sediment near the glacier terminus. The OPTV logs also reveal that the primary stratification both decreases in dip with depth (by up to 56° over 20 m) and rotates with depth (by up to 100° over 20 m) towards parallelism with the proximal lateral moraine. This transformation and the presence of relict layers of basally derived sediment raised into an englacial position – possibly involving thrusting – near the glacier's now stagnant terminus reveal a previously more dynamic glacier regime.

Glaciers in the Southern Patagonia Icefield (SPI) have been shrinking in recent decades, but due to a lack of field observations, understanding of the drivers of ablation is limited. We present a distributed surface energy balance model, forced with meteorological observations from a west–east transect located in the north of the SPI. Between October 2015 and June 2016, humid and warm on-glacier conditions prevailed on the western side compared to dry and cold conditions on the eastern side. Controls of ablation differ along the transect, although at glacier-wide scale sensible heat (mean of 72 W m−2 to the west and 51 W m−2 to the east) and net shortwave radiation (mean of 54 W m−2 to the west and 52 W m−2 to the east) provided the main energy sources. Net longwave radiation was an energy sink, while latent heat was the most spatially variable flux, being an energy sink in the east (−4 W m−2) and a source in the west (20 W m−2). Ablation was high, but at comparable elevations, it was greater to the west. These results provide new insights into the spatial variability of energy-balance fluxes and their control over the ablation of Patagonian glaciers.

Rock debris covers ~30% of glacier ablation areas in the Central Himalaya and modifies the impact of atmospheric conditions on mass balance. The thermal properties of supraglacial debris are diurnally variable but remain poorly constrained for monsoon-influenced glaciers over the timescale of the ablation season. We measured vertical debris profile temperatures at 12 sites on four glaciers in the Everest region with debris thickness ranging from 0.08 to 2.8 m. Typically, the length of the ice ablation season beneath supraglacial debris was 160 days (15 May to 22 October)—a month longer than the monsoon season. Debris temperature gradients were approximately linear (r2 > 0.83), measured as −40°C m–1 where debris was up to 0.1 m thick, −20°C m–1 for debris 0.1–0.5 m thick, and −4°C m–1 for debris greater than 0.5 m thick. Our results demonstrate that the influence of supraglacial debris on the temperature of the underlying ice surface, and therefore melt, is stable at a seasonal timescale and can be estimated from near-surface temperature. These results have the potential to greatly improve the representation of ablation in calculations of debris-covered glacier mass balance and projections of their response to climate change.

Turbulent fluxes make a substantial and growing contribution to the energy balance of ice surfaces globally, but are poorly constrained owing to challenges in estimating the aerodynamic roughness length (z0). Here, we used structure from motion (SfM) photogrammetry and terrestrial laser scanning (TLS) surveys to make plot-scale 2-D and 3-D microtopographic estimations of z0 and upscale these to map z0 across an ablating mountain glacier. At plot scales, we found spatial variability in z0 estimates of over two orders of magnitude with unpredictable z0 trajectories, even when classified into ice surface types. TLS-derived surface roughness exhibited strong relationships with plot-scale SfM z0 estimates. At the glacier scale, a consistent increase in z0 of ~0.1 mm d−1 was observed. Space-for-time substitution based on time since surface ice was exposed by snow melt confirmed this gradual increase in z0 over 60 d. These measurements permit us to propose a scale-dependent temporal z0 evolution model where unpredictable variability at the plot scale gives way to more predictable changes of z0 at the glacier scale. This model provides a critical step towards deriving spatially and temporally distributed representations of z0 that are currently lacking in the parameterisation of distributed glacier surface energy balance models.

Glaciers in the eastern Pamir have reportedly been gaining mass during recent decades, even though glaciers in most other regions in High Mountain Asia have been in recession. Questions still remain about whether the trend is strengthening or weakening, and how far the positive balances extend into the eastern Pamir. To address these gaps, we use three different digital elevation models to reconstruct glacier surface elevation changes over two periods (2000–09 and 2000–15/16). We characterize the eastern Pamir as a zone of transition from positive to negative mass balance with the boundary lying at the northern end of Kongur Tagh, and find that glaciers situated at higher elevations are those with the most positive balances. Most (67% of 55) glaciers displayed a net mass gain since the 21st century. This led to an increasing regional geodetic glacier mass balance from −0.06 ± 0.16 m w.e. a−1 in 2000–09 to 0.06 ± 0.04 m w.e. a−1 in 2000–15/16. Surge-type glaciers, which are prevalent in the eastern Pamir, showed fluctuations in mass balance on an individual scale during and after surges, but no statistical difference compared to non-surge-type glaciers when aggregated across the region.

While hot-water drilling is a well-established technique used to access the subsurface of ice masses, drilling into high-elevation (≳ 4000 m a.s.l.) debris-covered glaciers faces specific challenges. First, restricted transport capacity limits individual equipment items to a volume and mass that can be slung by small helicopters. Second, low atmospheric oxygen and pressure reduces the effectiveness of combustion, limiting a system's ability to pump and heat water. Third, thick supraglacial debris, which is both highly uneven and unstable, inhibits direct access to the ice surface, hinders the manoeuvring of equipment and limits secure sites for equipment placement. Fourth, englacial debris can slow the drilling rate such that continued drilling becomes impracticable and/or boreholes deviate substantially from vertical. Because of these challenges, field-based englacial and subglacial data required to calibrate numerical models of high-elevation debris-covered glaciers are scarce or absent. Here, we summarise our experiences of hot-water drilling over two field seasons (2017–2018) at the debris-covered Khumbu Glacier, Nepal, where we melted 27 boreholes up to 192 m length, at elevations between 4900 and 5200 m a.s.l. We describe the drilling equipment and operation, evaluate the effectiveness of our approach and suggest equipment and methodological adaptations for future use.

Fine resolution topographic data derived from methods such as Structure from Motion (SfM) and Multi-View Stereo (MVS) have the potential to provide detailed observations of geomorphological change, but have thus far been limited by the logistical constraints of conducting repeat surveys in the field. Here, we present the results from an automated time-lapse camera array, deployed around an ice-marginal lake on the western margin of the Greenland ice sheet. Fifteen cameras acquired imagery three-times per day over a 426 day period, yielding a dataset of ~19 000 images. From these data we derived 18 point clouds of the ice-margin across a range of seasons and successfully identified calving events (ranging from 234 to 1475 m2 in area and 815–8725 m3 in volume) induced by ice cliff undercutting at the waterline and the collapse of spalling flakes. Low ambient light levels, locally reflective surfaces and the large survey range hindered analysis of smaller scale ice-margin dynamics. Nevertheless, this study demonstrates that an integrated SfM-MVS and time-lapse approach can be employed to generate long-term 3-D topographic datasets and thus quantify ice-margin dynamics at a fine spatio-temporal scale. This approach provides a template for future studies of geomorphological change.

Measurements of glacier ice cliff evolution are sparse, but where they do exist, they indicate that such areas of exposed ice contribute a disproportionate amount of melt to the glacier ablation budget. We used Structure from Motion photogrammetry with Multi-View Stereo to derive 3-D point clouds for nine ice cliffs on Khumbu Glacier, Nepal (in November 2015, May 2016 and October 2016). By differencing these clouds, we could quantify the magnitude, seasonality and spatial variability of ice cliff retreat. Mean retreat rates of 0.30–1.49 cm d−1 were observed during the winter interval (November 2015–May 2016) and 0.74–5.18 cm d−1 were observed during the summer (May 2016–October 2016). Four ice cliffs, which all featured supraglacial ponds, persisted over the full study period. In contrast, ice cliffs without a pond or with a steep back-slope degraded over the same period. The rate of thermo-erosional undercutting was over double that of subaerial retreat. Overall, 3-D topographic differencing allowed an improved process-based understanding of cliff evolution and cliff-pond coupling, which will become increasingly important for monitoring and modelling the evolution of thinning debris-covered glaciers.

We present satellite-derived velocity patterns for the two contrasting melt seasons of 2009–10 across Russell Glacier catchment, a western, land-terminating sector of the Greenland ice sheet which encompasses the K(angerlussuaq)-transect. Results highlight great spatial heterogeneity in flow, indicating that structural controls such as bedrock geometry govern ice discharge into individual outlet troughs. Results also reveal strong seasonal flow variability extending 57 km up-glacier to 1200 m elevation, with the largest acceleration (100% over 11 days) occurring within 10 km of the margin coincident with spring melt. By late July 2010, 2 weeks before peak melt and runoff, 48 % of the 2400 km2 catchment had slowed to less than the winter mean. This observation supports the hypothesis that the subglacial hydrological system evolves from an inefficient distributed to an efficient drainage system, regulating flow dynamics. Despite this, the cumulative surface flux over the record melt year of 2010 was still greater compared with the perturbation over the average melt year of 2009. This study supports the proposition that local surface meltwater runoff couples to basal hydrology driving ice-sheet dynamics, and although the effect is nonlinear, our observations indicate that greater meltwater runoff yields increased net flux over this sector of the ice sheet.

Relationships between stable isotopes (δD–δ18O), ice facies and glacier structures have hitherto gone untested in the mid-latitude maritime glaciers of the Southern Hemisphere. Here, we present δD–δ18O values as part of a broader study of the structural glaciology of Fox Glacier, New Zealand. We analyzed 94 samples of δD–δ18O from a range of ice facies to investigate whether isotopes have potential for structural glaciological studies of a rapidly deforming glacier. The δD–δ18O measurements were aided by structural mapping and imagery from terminus time-lapse cameras. The current retreat phase was preceded by an advance of 1 km between 1984 and 2009, with the isotopic sampling and analysis undertaken at the end of that advance (2010/11). Stable isotopes from debris-bearing shear planes near the terminus, interpreted as thrust faults, are isotopically enriched compared with the surrounding ice. When plotted on co-isotopic diagrams (δD–δ18O), ice sampled from the shear planes appears to show a subtle, but distinctive isotopic signal compared with the surrounding clean ice on the lower glacier. Hence, stable isotopes (δD–δ18O) have potential within the structural glaciology field, but larger sample numbers than reported here may be required to establish isotopic contrasts between a broad range of ice facies and glacier structures.