The western Nyainqentanglha Range is located in the south-eastern centre of the Tibetan Plateau. Its north-western slopes drain into Lake Nam Co. The region is of special interest for glacio-climatological research as it is influenced by both the continental climate of Central Asia and the Indian Monsoon system, and situated at the transition zone between temperate and subcontinental glaciers. A glacier inventory for the whole mountain range was generated for the year around 2001 using automated remote sensing and GIS techniques based on Landsat ETM+ and SRTM3 DEM data. Glacier change analysis was based on data from Hexagon KH-9 and Landsat MSS (both 1976), Metric Camera (1984), and Landsat TM/ETM+ (1991, 2001, 2005, 2009). Manual adjustment was especially necessary for delineating the debris-covered glaciers and the glaciers on the panchromatic Hexagon data. In the years around 2001 the whole mountain range contained about 960 glaciers covering an area of 795.6 ± 22.3 km² while the ice in the drainage basin of Nam Co covered 198.1 ± 5.6 km². The median elevation of the glaciers was about 5800 m with the majority terminating around 5600 m. Five glaciers with debris-covered tongues terminated lower than 5200 m. The glacier area decreased by −6.1 ± 3% between 1976 and 2001. This is less than reported in previous studies based on the 1970s topographic maps and Landsat data from 2000. Glaciers continued to shrink during the period 2001–2009. No advancing glaciers were detected. Detailed length measurements for five glaciers indicated a retreat of around 10 m per year (1976–2009). Ice cover is higher south-east of the mountain ridge which reflects the windward direction to the monsoon. The temperature increase during the ablation period was probably the main driver of glacier wastage, but the complex glacier-climate interactions need further investigation.

Glacier changes are a vivid example of how environmental systems react to a changing climate. Distributed surface mass balance models, which translate the meteorological conditions on glaciers into local melting rates, help to attribute and detect glacier mass and volume responses to changes in the climate drivers. A well-calibrated model is a suitable test bed for sensitivity, detection, and attribution analyses for many scientific applications and often serves as a tool for quantifying the inherent uncertainties. Here, we present the open-source COupled Snowpack and Ice surface energy and mass balance model in PYthon (COSIPY), which provides a flexible and user-friendly framework for modeling distributed snow and glacier mass changes. The model has a modular structure so that the exchange of routines or parameterizations of physical processes is possible with little effort for the user. The framework consists of a computational kernel, which forms the runtime environment and takes care of the initialization, the input–output routines, and the parallelization, as well as the grid and data structures. This structure offers maximum flexibility without having to worry about the internal numerical flow. The adaptive subsurface scheme allows an efficient and fast calculation of the otherwise computationally demanding fundamental equations. The surface energy balance scheme uses established standard parameterizations for radiation as well as for the energy exchange between atmosphere and surface. The schemes are coupled by solving both surface energy balance and subsurface fluxes iteratively such that consistent surface skin temperature is returned at the interface. COSIPY uses a one-dimensional approach limited to the vertical fluxes of energy and matter but neglects any lateral processes. Accordingly, the model can be easily set up in parallel computational environments for calculating both energy balance and climatic surface mass balance of glacier surfaces based on flexible horizontal grids and with varying temporal resolution. The model is made available on a freely accessible site and can be used for non-profit purposes. Scientists are encouraged to actively participate in the extension and improvement of the model code.

Climate variables that control the annual cycle of the surface energy and mass balance on Zhadang glacier in the central Tibetan Plateau were examined over a 2 year period using a physically based energy-balance model forced by routine meteorological data. The modelled results agree with measured values of albedo, incoming longwave radiation, surface temperature and surface level of the glacier. For the whole observation period, the radiation component dominated (82%) the total surface energy heat fluxes. This was followed by turbulent sensible (10%) and latent heat (6%) fluxes. Subsurface heat flux represented a very minor proportion (2%) of the total heat flux. The sensitivity of specific mass balance was examined by perturbations of temperature (±1 K), relative humidity (±20%) and precipitation (±20%). The results indicate that the specific mass balance is more sensitive to changes in precipitation than to other variables. The main seasonal variations in the energy balance were in the two radiation components (net shortwave radiation and net longwave radiation) and these controlled whether surface melting occurred. A dramatic difference in summer mass balance between 2010 and 2011 indicates that the glacier surface mass balance was closely related to precipitation seasonality and form (proportion of snowfall and rainfall).

Bedrock overdeepenings exposed by continued glacial retreat can store precipitation and meltwater, potentially leading to the formation of new proglacial lakes. These lakes may pose threats of glacial lake outburst floods (GLOFs) in high mountain areas, particularly if new lakes form in geomorphological setups prone to triggering events such as landslides or moraine collapses. We present the first complete inventory for future glacial lakes in High Mountain Asia by computing the subglacial bedrock for ~100 000 glaciers and estimating overdeepening area, volume and impact hazard for the larger potential lakes. We detect 25 285 overdeepenings larger than 10⁴m²with a volume of 99.1 ± 28.6 km³covering an area of 2683 ± 773.8 km². For the 2700 overdeepenings larger than 10⁵m², we assess the lake predisposition for mass-movement impacts that could trigger a GLOF by estimating the hazard of material detaching from surrounding slopes. Our findings indicate a shift in lake area, volume and GLOF hazard from the southwestern Himalayan region toward the Karakoram. The results of this study can be used for anticipating emerging threats and potentials connected to glacial lakes and as a basis for further studies at suspected GLOF hazard hotspots.

We present a new open-source, collaborative “COupled Snowpack and Ice surface energy and MAss balance model” (COSIMA) that is evaluated for Zhadang glacier, Tibetan Plateau. The model is calibrated, run, and validated based on in situ measurements and atmospheric model data from the High Asia Refined analysis (HAR) over the period April 2009 to June 2012. Results for the model runs forced by both in situ measurements and HAR agree well with observations of various atmospheric, glaciological, surface, and subsurface parameters on the glacier. A time-lapse camera system next to the glacier provides a 3-year image time series of the mean transient snow line altitude and the snow cover pattern, which is used for the spatial and temporal validation of the model. The model output corresponds very well to the observed temporal and spatial snow cover variability. The model is then run for a 10-year period of October 2001 to September 2011 forced with HAR data. In general, the radiation components dominate...