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.

Recent progress has been made in quantifying snowmelt in the Himalaya. Although the conditions are favorable for refreezing, little is known about the spatial variability of meltwater refreezing, hindering a complete understanding of seasonal snowmelt dynamics. This study aims to improve our understanding about how refreezing varies in space and time. We simulated refreezing with the seNorge (v2.0) snow model for the Langtang catchment, Nepalese Himalaya, covering a 5-year period. Meteorological forcing data were derived from a unique elaborate network of meteorological stations and high-resolution meteorological simulations. The results show that the annual catchment average refreezing amounts to 122 mm w.e. (21% of the melt), and varies strongly in space depending on elevation and aspect. In addition, there is a seasonal altitudinal variability related to air temperature and snow depth, with most refreezing during the early melt season. Substantial intra-annual variability resulted from fluctuations in snowfall. Daily refreezing simulations decreased by 84% (annual catchment average of 19 mm w.e.) compared to hourly simulations, emphasizing the importance of using sub-daily time steps to capture melt–refreeze cycles. Climate sensitivity experiments revealed that refreezing is highly sensitive to changes in air temperature as a 2°C increase leads to a refreezing decrease of 35%.

We present changes of the ice margin in Northwest Greenland at the Eastern part of the Nunatarssuaq Ice Cap (NIC) over six decades. The ice margin in this area terminates as a near-vertical ice cliff of between 9 and 33 m thickness. During the years 1954–1957 and in 1965 multi-disciplinary studies were performed. We digitise and orthorectify material, that is often difficult to access, in order to use the historical data as an absolute starting point of our change assessment. We compare the cliff morphology of the mid-1950s and the mid-1960s with various time-steps between 1985 and 2017. The studied ice margin remained remarkably constant with very subtle changes of changing sign: rather slow advance rates are reported from the 1950s and 1960s that accelerated until 1985 and were followed by a general retreat until 2012 and a subsequent advance until 2017. Thickness changes are negative throughout the entire time-period, however, different rates of thinning are shown and there is a positive relationship with air temperature anomalies. Compared to similar elevations on the adjacent Greenland ice sheet, we find significantly weaker thinning rates at the NIC.

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.

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%.

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.

Estimates of melt from debris-covered glaciers require distributed estimates of meteorological variables and air temperature in particular. Meteorological data are scarce for this environment, and spatial variability of temperature over debris is poorly understood. Based on multiple measurements of air and surface temperature from three ablation seasons (2012–14) we investigate the variability of temperature over Lirung Glacier, Nepal, in order to reveal how air temperature is affected by the debris cover and improve ways to extrapolate it. We investigate how much on-glacier temperature deviates from that predicted from a valley lapse rate (LR), analyse on-glacier LRs and test regression models of air temperature and surface temperature. Air temperature over the debris-covered glacier tongue is much higher than what a valley LR would prescribe, so an extrapolation from off-glacier stations is not applicable. An on-glacier LR is clearly defined at night, with strong correlation, but not during the day, when the warming debris disrupts the elevation control. An alternative to derive daytime air temperature is to use a relationship between air and surface temperature, as previously suggested. We find strong variability during daytime that should be accounted for if these regressions are used for temperature extrapolation.

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.