Ice cliffs are features along ice sheet margins, along tropical mountain glaciers, at termini of mountain glaciers and on debris-covered glacier tongues, that have received scattered attention in literature. They cover small relative areas of glacier or margin surface respectively, but have been involved in two apparent anomalies. On the one hand, they have been identified as potential hotspots of extreme melt rates on debris-covered tongues contributing to their relatively rapid ablation, compared to the surrounding glacier surface. On the other hand, they appear where the ice margin is stable (or temporarily advancing) even under conditions of negative mass balance. In this manuscript, we recapitulate why ice cliffs remain interesting features to investigate and what we know about them so far. We conclude by suggesting to further investigate their genesis and variable morphology and their potential as windows into past climates and processes.

Ice cliffs and ponds on debris-covered glaciers have received increased attention due to their role in amplifying local melt. However, very few studies have looked at these features on the catchment scale to determine their patterns and changes in space and time. We have compiled a detailed inventory of cliffs and ponds in the Langtang catchment, central Himalaya, from six high-resolution satellite orthoimages and DEMs between 2006 and 2015, and a historic orthophoto from 1974. Cliffs cover between 1.4% (± 0.4%) in the dry and 3.4% (± 0.9%) in the wet seasons and ponds between 0.6% (± 0.1%) and 1.6% (± 0.3%) of the total debris-covered tongues. We find large variations between seasons, as cliffs and ponds tend to grow in the wetter monsoon period, but there is no obvious trend in total area over the study period. The inventory further shows that cliffs are predominately north-facing irrespective of the glacier flow direction. Both cliffs and ponds appear in higher densities several hundred metres from the terminus in areas where tributaries reach the main glacier tongue. On the largest glacier in the catchment ~10% of all cliffs and ponds persisted over nearly a decade.

Mass losses originating from supraglacial ice cliffs at the lower tongues of debris-covered glaciers are a potentially large component of the mass balance, but have rarely been quantified. In this study, we develop a method to estimate ice cliff volume losses based on high-resolution topographic data derived from terrestrial and aerial photogrammetry. We apply our method to six cliffs monitored in May and October 2013 and 2014 using four different topographic datasets collected over the debris-covered Lirung Glacier of the Nepalese Himalayas. During the monsoon, the cliff mean backwasting rate was relatively consistent in 2013 (3.8 ± 0.3 cm w.e. d−1) and more heterogeneous among cliffs in 2014 (3.1 ± 0.7 cm w.e. d−1), and the geometric variations between cliffs are larger. Their mean backwasting rate is significantly lower in winter (October 2013–May 2014), at 1.0 ± 0.3 cm w.e. d−1. These results are consistent with estimates of cliff ablation from an energy-balance model developed in a previous study. The ice cliffs lose mass at rates six times higher than estimates of glacier-wide melt under debris, which seems to confirm that ice cliffs provide a large contribution to total glacier melt.

Air temperature is a key control of processes affecting snow and glaciers in high-elevation catchments, including melt, snowfall and sublimation. It is therefore a key input variable to models of land–surface–atmosphere interaction. Despite this importance, its spatial variability is poorly understood and simple assumptions are made to extrapolate it from point observations to the catchment scale. We use a dataset of 2.75 years of air temperature measurements (from May 2012 to November 2014) at a network of up to 27 locations in the Langtang River, Nepal, catchment to investigate air temperature seasonality and consistency between years. We use observations from high elevations and from the easternmost section of the basin to corroborate previous findings of shallow lapse rates. Seasonal variability is strong, with shallowest lapse rates during the monsoon season. Diurnal variability is also strong and should be taken into account since processes such as melt have a pronounced diurnal variability. Use of seasonal lapse rates seems crucial for glacio-hydrological modelling, but seasonal lapse rates seem stable over the 2–3 years investigated. Lateral variability at transects across valley is high and dominated by aspect, with south-facing sites being warmer than north-facing sites and deviations from the fitted lapse rates of up to several degrees. Local factors (e.g. topographic shading) can reduce or enhance this effect. The interplay of radiation, aspect and elevation should be further investigated with high-elevation transects.

Supraglacial ponds on debris-covered glaciers present a mechanism of atmosphere/glacier energy transfer that is poorly studied, and only conceptually included in mass-balance studies of debris-covered glaciers. This research advances previous efforts to develop a model of mass and energy balance for supraglacial ponds by applying a free-convection approach to account for energy exchanges at the subaqueous bare-ice surfaces. We develop the model using field data from a pond on Lirung Glacier, Nepal, that was monitored during the 2013 and 2014 monsoon periods. Sensitivity testing is performed for several key parameters, and alternative melt algorithms are compared with the model. The pond acts as a significant recipient of energy for the glacier system, and actively participates in the glacier’s hydrologic system during the monsoon. Melt rates are 2-4 cm d-1 (total of 98.5 m3 over the study period) for bare ice in contact with the pond, and <1 mmd-1 (total of 10.6m3) for the saturated debris zone. The majority of absorbed atmospheric energy leaves the pond system through englacial conduits, delivering sufficient energy to melt 2612 m3 additional ice over the study period (38.4 m3 d-1). Such melting might be expected to lead to subsidence of the glacier surface. Supraglacial ponds efficiently convey atmospheric energy to the glacier’s interior and rapidly promote the downwasting process.

Ice cliffs might be partly responsible for the high mass losses of debris-covered glaciers in the Hindu Kush-Karakoram-Himalaya region. The few existing models of cliff backwasting are point-scale models applied at few locations or assume cliffs to be planes with constant slope and aspect, a major simplification given the complex surfaces of most cliffs. We develop the first grid-based model of cliff backwasting for two cliffs on debris-covered Lirung Glacier, Nepal. The model includes an improved representation of shortwave and longwave radiation, and their interplay with the glacier topography. Shortwave radiation varies considerably across the two cliffs, mostly due to direct radiation. Diffuse radiation is the major shortwave component, as the direct component is strongly reduced by the cliffs’ aspect and slope through self-shading. Incoming longwave radiation is higher than the total incoming shortwave flux, due to radiation emitted by the surrounding terrain, which is 25% of the incoming flux. Melt is highly variable in space, suggesting that simple models provide inaccurate estimates of total melt volumes. Although only representing 0.09% of the glacier tongue area, the total melt at the two cliffs over the measurement period is 2313 and 8282 m3, 1.23% of the total melt simulated by a glacio-hydrological model for the glacier’s tongue.

Ice cliffs have been identified as a reason for higher ablation rates on debris-covered glaciers than are implied by the insulation effects of the debris. This study aims to improve our understanding of cliff backwasting, and the role of radiative fluxes in particular. An energy-balance model is forced with new data gathered in May and October 2013 on Lirung Glacier, Nepalese Himalaya. Observations show substantial variability in melt between cliffs, between locations on any cliff and between seasons. Using a high-resolution digital elevation model we calculate longwave fluxes incident to the cliff from surrounding terrain and include the effect of local shading on shortwave radiation. This is an advance over previous studies, that made simplified assumptions on cliff geometry and radiative fluxes. Measured melt rates varied between 3.25 and 8.6 cm d−1 in May and 0.18 and 1.34 cm d−1 in October. Model results reproduce the strong variability in space and time, suggesting considerable differences in radiative fluxes over one cliff. In October the model fails to reproduce stake readings, probably due to the lack of a refreezing component. Disregarding local topography can lead to overestimation of melt at the point scale by up to ∼9%.