A glacierized terrain comprises different land covers, and their mapping using satellite data is challenged by their spectral similarity. We propose a hierarchical knowledge-based classification (HKBC) approach for differentiation of glacier terrain classes and mapping of glacier boundaries, using Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) imagery and Global Digital Elevation Model (GDEM). The methodology was tested over Kolahoi Glacier, Kashmir Himalaya. For the sequential extraction of various glacier terrain classes, several input layers were generated from the primary datasets by applying image-processing techniques. Noticeable differences in temperature and spectral response between supraglacial debris and periglacial debris facilitated the development of a thermal glacier mask and normalized-difference debris index, which together with slope enabled their differentiation. These and the other layers were then used in several discrete tests in HKBC, to map various glacier terrain classes. An ASTER visible near-infrared image and 42 field points were used to validate results. The proposed approach satisfactorily classified all the glacier terrain classes with an overall accuracy of 89%. The Z-test reveals that results obtained from HKBC are significantly (at 95% confidence level) better than those from a maximum likelihood classifier (MLC). Glacier boundaries obtained from HKBC were found to be plausibly better than those obtained from MLC and visual interpretation.

We investigated glacier changes in the Karatal river basin, the largest basin in Zhetysu (Dzhungar) Alatau, Kazakhstan, for the periods 1956–89, 1989–2001 and 2001–12, based on Landsat TM/ETM+ data analysis. In 1989, we found 243 glaciers with a total area of 142.8 km2; by 2012 these had shrunk to 214 glaciers with a total area of 109.3 km2, a decrease of 33.5 km2 over 23 years (1.02%a-1). This very high shrinkage rate is likely connected with a general trend of increasing temperatures, and small glaciers being situated at the relatively low altitude of the outer Zhetysu Alatau ranges. We also analyzed the shrinkage rate of glaciers based on their differences in size, altitude and aspect of slopes, as well as other topographic parameters, in four sub-basins where glacier shrinkage varied between 18% and 39%. Weather-station climate data showed a significant temperature increase and stable precipitation trends over the study period. We conclude that glacierized areas of the Karatal river basin are located in the most unfavorable conditions for glaciation, and as a result showed a higher shrinkage rate than other glacierized areas of the Tien Shan from 1956 to 2012.

The Qinghai-Tibetan Plateau (QTP) is characterized by a cold climate and a large number of lakes. The long ice season necessitates study of the widespread ice covers in the region. An unprecedented multidisciplinary field campaign was conducted on lake ice processes in the central QTP during the period 2019–13. The study lake generally froze up in late October or early November, and broke up in mid or late April, with a maximum ice thickness of 50–70 cm. The mass balances at both ice surface and bottom were measured continuously. Significant ice surface sublimation/ablation was detected and accounted for up to 40% of the whole ice thickness over the ice season. A simple heat-transfer model was developed for the surface ice loss. The calculated values were in good agreement with the observations. They also indicated that atmospheric conditions, including low air humidity and prevailing strong winds, are the primary drivers of the ice surface sublimation.

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.

An assessment of glacier shrinkage (reduction of area) for all of High Mountain Asia requires a complete compilation of measured rates of change and also a methodology for objective comparison of rates. I present a compilation from 155 publications reporting glacier area changes, and also a methodology that overcomes the main obstacles hindering comparison. Glacier areas are not always assigned uncertainties, and this problem is addressed with an error model derived from published estimates. The problem of discordant survey dates is addressed by interpolating measured areas to fixed dates at pentadal intervals. Interpolation error depends only incoherently on the time span between measurements, but strongly on glacier size: smaller glaciers, in addition to changing more rapidly on average, exhibit more variable rates of change. The overlapping boundaries of study regions are reconciled by mapping all of the information to a 0.5° geographical grid. When coupled with glacier area information from the Randolph Glacier Inventory, the widely observed inverse dependence of shrinkage rates on glacier size shows promise as a tool for treating incomplete spatial coverage. Over High Mountain Asia as a whole from 1960 to 2010, the unweighted average shrinkage rate is –0.57% a–1, but corrections for variable glacier size raise the average to –0.34% a–1, and filling unmeasured gridcells with rates based on size dependence alters the latter estimate to –0.40% a–1. The uncertainties in these rates are large. The Karakoram anomaly is found to be a zonal feature extending well to the east of the Karakoram proper.

This study outlines a consistent methodology for identifying glacier surfaces from Landsat 5, 7 and 8 imagery that is applied to map all mainland North Asian glaciers, providing the first methodologically consistent and complete glacier inventory for the region ~2010. We identify 5065 glaciers covering a planimetric area of 2326 ± 186 km2, most of which is located in the Altai mountain subregion. The total glacier count is 15% higher, but the total glacier area is 32 ±11.6% lower, than the estimated glacier coverage provided in version 4.0 of the Randolph Glacier Inventory. We investigate the distribution of glacier size within North Asia and find that the majority of glaciers (82%) are smaller than 0.5 km2 but only account for a third of the total glacier area, with the largest 1 % (60 glaciers ≥ 5 km2 ) accounting for 28% of the total area. We present hypsometric characterizations of North Asian glaciers, largely substantiating existing findings that glaciers in this region are dominated by cold, relatively dry conditions. We provide a detailed assessment of errors and determine the uncertainty in our area estimate to be ±8.0%, with snow-cover uncertainty the largest contributing factor. Based on this assessment, the new glacier inventory presented here is more complete and of higher quality than other currently available data sources.

Glaciers in the Bugyai Kangri are located in a transition zone from southeast Tibet, where monsoonal temperate glaciers dominate, to inner Tibet, where continental glaciers dominate. Here we analyze glacier and glacial lake changes in this region using multi-year inventories based on Landsat images from 1981–2013. Results show that the total area of 141 glaciers in the region decreased by 30.44 ±0.89 km2 from 198.35 ±9.54 km2 (1980s) to 167.93 ±4.52 km2 (2010s). The annual area shrinkage rate (–0.48% a–1) is lower than that reported for southeastern Tibet but higher than that of inner Tibet. Both the number and total area of glacial lakes increased between 1981 and 2013. Among all lakes, proglacial lakes contribute most (~81 %) to the expansion. The total area of ten proglacial lakes increased by 150.3 ± 13.17% and of these ten lakes the four that expanded most sharply showed increased calving at their upper margins, resulting in more rapid retreat of lake-terminating glaciers than land-terminating glaciers. Owing to rapid calving, several lakes may undergo further growth in the near future, increasing the potential risk of glacial lake outburst floods.

High-resolution δ18O records from a Geladaindong mountain ice core spanning the period 1477-1982 were used to investigate past temperature variations in the Yangtze River source region of the central Tibetan Plateau (TP). Annual ice-core δ18O records were positively correlated with temperature data from nearby meteorological stations, suggesting that the δ18O record represented the air temperature in the region. A generally increasing temperature trend over the past 500 years was identified, with amplified warming during the 20th century. A colder stage, spanning before the 1850s, was found to represent the Little Ice Age with colder periods occurring during the 1470s–1500s, 1580s–1660s, 1700s–20s and 1770s–1840s. Compared with other temperature records from the TP and the Northern Hemisphere, the Geladaindong ice-core record suggested that the regional climate of the central TP experienced a stronger warming trend during the 20th century than other regions. In addition, a positive relationship between the Geladaindong δ18 O values and the North Atlantic Oscillation index, combined with a wavelet analysis of δ18 O records, indicated that there was a potential atmospheric teleconnection between the North Atlantic and the central TP.

Detailed studies on temporal changes of crevasses and their linkage with glacier dynamics are scarce in the Himalayan context. Observations of temporally changing surficial crevasse patterns and their orientations are suggestive of the processes that determine seasonal glacier flow characteristics. In the present study, on a Himalayan valley glacier, changing crevasse patterns and orientations were detected and mapped on Landsat 8 images in an automated procedure using the ratio of Thermal Infrared Sensor (TIRS) band 10 to Optical Land Imager (OLI) shortwave infrared (SWIR) band 6. The ratio was capable of mapping even crevasses falling under mountain shadows. Differential GPS observations suggested an average error of 3.65% and root-mean-square error of 6.32m in crevasse lengths. A year-round observation of these crevasses, coupled with field-based surface velocity measurements, provided some interesting interpretations of seasonal glacier dynamics.

This study presents a reconstruction of the seasonal mass balance of Glacier No. 354, located in the Akshiirak range, Kyrgyzstan, from 2003 to 2014. We use a distributed accumulation and temperature-index melt model driven by daily air temperature and precipitation from a nearby meteorological station. The model is calibrated with in situ measurements of the annual mass balance collected from 2011 to 2014. The snow-cover depletion pattern observed using satellite imagery provides additional information on the dynamics of mass change throughout the melting season. Two digital elevation models derived from high-resolution satellite stereo images acquired in 2003 and 2012 are used to calculate glacier volume change for the corresponding period. The geodetic mass change thus derived is used to validate the modelled cumulative glacier-wide balance. For the period 2003–12 we find a cumulative mass balance of –0.40±10mw.e.a-1. This result agrees well with the geodetic balance of –0.48±0.07mw.e.a-1over the same period.

Debris-covered glaciers play an important role in the high-altitude water cycle in the Himalaya, yet their dynamics are poorly understood, partly because of the difficult fieldwork conditions. In this study we therefore deploy an unmanned aerial vehicle (UAV) three times (May 2013, October 2013 and May 2014) over the debris-covered Lirung Glacier in Nepal. The acquired data are processed into orthomosaics and elevation models by a Structure from Motion workflow, and seasonal surface velocity is derived using frequency cross-correlation. In order to obtain optimal surface velocity products, the effects of different input data and correlator configurations are evaluated, which reveals that the orthomosaic as input paired with moderate correlator settings provides the best results. The glacier has considerable spatial and seasonal differences in surface velocity, with maximum summer and winter velocities 6 and 2.5 m a-1, respectively, in the upper part of the tongue, while the lower part is nearly stagnant. It is hypothesized that the higher velocities during summer are caused by basal sliding due to increased lubrication of the bed. We conclude that UAVs have great potential to quantify seasonal and annual variations in flow and can help to further our understanding of debris-covered glaciers.

With the exception of northern India, there are few, if any, consistent data records relating to avalanche activity in the high mountains of Asia. However, records do exist of avalanche fatalities in the region, contained in mountaineering expedition reports. In this paper, I review and analyze statistics of avalanche fatalities (both snow and ice) in the high mountains of Asia (Himalaya, Karakoram, Pamir, Hindu Kush, Tien Shan, Dazu Shan) from 1895 to 2014. The data are stratified according to accident cause, geographical region (Nepal-Tibet (Xizang), Pakistan, India, China, Central Asia), mountain range, personnel (hired or expedition members) and terrain. The character of the accidents is compared with data from North America and Europe. The data show that the important risk components are the temporal and spatial exposure probabilities. It is shown that human actions and decisions govern the pattern of fatal avalanches in the high mountains of Asia.

Surface digital elevation models (DEMs) and slope-related estimates of glacier thickness enable modelling of glacier-bed topographies over large ice-covered areas. Due to the erosive power of glaciers, such bed topographies can contain numerous overdeepenings, which when exposed following glacier retreat may fill with water and form new lakes. In this study, the bed overdeepenings for ~28 000 glaciers (40 775 km2) of the Himalaya-Karakoram region are modelled using GlabTop2 (Glacier Bed Topography model version 2), in which ice thickness is inferred from surface slope by parameterizing basal shear stress as a function of elevation range for each glacier. The modelled ice thicknesses are uncertain (±30%), but spatial patterns of ice thickness and bed elevation primarily depend on surface slopes as derived from the DEM and, hence, are more robust. About 16 000 overdeepenings larger than 104m2 were detected in the modelled glacier beds, covering an area of ~2200 km2 and having a volume of ~120km3 (3-4% of present-day glacier volume). About 5000 of these overdeepenings (1800 km2) have a volume larger than 106m3. The results presented here are useful for anticipating landscape evolution and potential future lake formation with associated opportunities (tourism, hydropower) and risks (lake outbursts).

We have compiled and analysed available records and data on the shrinkage of Satopanth (SPG) and Bhagirath Kharak (BKG) Glaciers, Uttarakhand, India, during the period 1936–2013. We estimate the mean retreat rates of the snouts of SPG and BKG for this period at 9.7±0.8m a–1 and 7.0±0.6m a–1 respectively. We have also revised the estimates of the area vacated during the period 1956–2013 to be 0.27 ± 0.05 km2 and 0.17 ± 0.04 km2 for SPG and BKG respectively, corresponding to front-averaged retreat rates of 5.7±0.6m a–1 and 6.0±0.9m a–1. The study revealed an average thinning of glacial ice in the lower ablation zone of SPG of 9 ± 11 m in the past 51 years. We observed that while the fronts of SPG and BKG depicted in the Survey of India topographic map published in 1962 are inconsistent with other available records, the elevation contours are consistent with them.

Using the Moderate Resolution Imaging Spectroradiometer (MODIS) Level 1 radiance Swath Data (MOD02QKM) with a spatial resolution of 250 m, we derive snowlines during July–September 2001–12 for several mountain ranges distributed across the Tibetan Plateau (TP). Radiance bands 1 and 2 are projected to the study area and processed automatically. The discrimination between snow and ice is done using a k-mean cluster analysis and the snowlines are delineated based on a fixed percentile of the snow-cover altitude. The highest transient snowline altitude is then taken as a proxy for the equilibrium-line altitude (ELA). In the absence of measured glaciological, meteorological or hydrological data, our ELA time series enable better understanding of atmosphere-cryosphere couplings on the TP. Interannual ELA variability is linked to local and remote climate indices using a correlation analysis. Southerly flow and higher temperatures are linked with a higher ELA in most regions. Eastern and Trans-Himalayan sites show positive correlations between winter temperatures and ELA. As winter temperatures are substantially below zero, this suggests an enhancement of winter sublimation as opposed to a reduction in accumulation. It appears that large-scale atmospheric forcing has varying and sometimes opposite influences on the annual ELA in different regions on the TP.

Understanding the hydrology of glacierized catchments is an important step in assessing the vulnerability of water resources to a changing climate. While there have been increased efforts recently to understand the dynamics of Asia's cryosphere, glacier melt dynamics and hydrograph separation of river discharge are open questions. A multi-year, multi-seasonal dataset of water chemistry from the Langtang Valley, Nepal, is used to explore water sources and flow paths that contribute to Langtang River discharge. Differences in hydrochemistry of samples from debris-free Khimsung Glacier and debris-covered Lirung Glacier demonstrate the effect of debris cover on glacier outflow. Additional data show seasonal transitions in the composition of Langtang River discharge. End-member mixing analysis (EMMA) using geochemical and isotopic tracers suggests that reacted meltwater contributes the majority of flow during most of the year, with the exception of the summer when unreacted meltwater and precipitation dominate streamflow. We hypothesize our dataset is missing characteristic monsoon water and utilize a Late May river sample as a proxy for precipitation-influenced groundwater in the EMMA. Results offer insight into the plausibility of flow sources and pathways in the basin.

Despite being a key location for paleoglaciological research in north-central Asia, with the largest number of modern and Pleistocene glaciers, and in the transition zone between the humid Russian Altai and dry Gobi Altai, little is known about the precise extent and timing of Holocene and late Pleistocene glaciations in western Mongolia. Here we present detailed information on the distribution of modern and late Holocene glaciers, and new results addressing the geomorphological differentiation and numerical dating (by optically stimulated luminescence, OSL) of Pleistocene glacial sequences in these areas. For the Mongolian Altai, geochronological results suggest large ice advances correlative to marine isotope stages (MIS) 4 and 2. This is in contrast to results from the Khangai mountains, central Mongolia, showing that significant ice advances additionally occurred during MIS3. During the Pleistocene, glacial equilibrium-line altitudes (ELAs) were ~500 to >1000m lower in the more humid portion of the Russian and western Mongolian Altai, compared to 300-600 m in the drier ranges of the eastern Mongolian Altai. Pleistocene ELAs in the Khangai mountains were depressed by 700-1000 m, suggesting more humid conditions at times of major glaciation than in the eastern Mongolian Altai. This paleo-ELA pattern reveals that the precipitation gradient from the drier to the more humid regions was more pronounced during glacial times than at present.

A long-standing problem in avalanche dynamics is to model the flow of a mixed flowing/powder avalanche. Here we use the thermodynamic concept of configurational energy to describe the blow-out of air from the avalanche core. Configurational energy is the mean potential energy associated with the location of snow and ice particles in the avalanche core. As such, configurational energy determines the avalanche flow density. Expansion of the particle ensemble reduces the flow density and leads to the intake of air. Compression of the particle ensemble causes the blow-out of the intaken air, now laden with ice dust. Once formed, the cloud moves independently of the flowing avalanche with the initial momentum acquired in the core. Configurational energy changes in the avalanche core are therefore intimately related to the formation of the powder suspension cloud. In this paper we use the concept of configurational energy to predict the mass of air taken into and blown out of the core. This requires calculating the dispersive pressure arising from random particle movements and configuration changes related to the expansion and collapse of the flowing particle ensemble. The ice avalanche that struck the Everest base camp on 25 April 2015 is simulated using the proposed concept.

A temperature reconstruction in the glacierized Himalaya over the past centuries using glacial length fluctuation records is challenging due to the abundance of debris-covered glaciers and a scarcity of glacial length fluctuation data. Using idealized flowline model simulations, we show that supraglacial debris cover significantly alters the length fluctuations only when the debris cover is very thick. An expanded database of length fluctuation records for 43 glaciers in the Himalaya and Karakoram is compiled and a standard linear inversion procedure is applied to a subset of 34 glaciers in this database. The reconstructed temperature anomaly during 1860-2010 indicates a continued warming of the region with a total temperature change of ~1.6 K. A close resemblance of the regional temperature anomaly to global trends is seen.

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.