A glacier’s most fundamental altitude is the equilibrium-line altitude (ELA) because it divides the glacier into ablation and accumulation areas. The best parameterization of the ELA for glacier inventory is the balanced-budget ELA. We discuss direct estimation of balanced-budget ELA from mass-balance data for individual glaciers, and indirect estimation of balanced-budget ELA from simple topographic parameters available from the World Glacier Inventory (WGI), i.e. the area-median and maximum and minimum altitudes. Mass balance and ELA for individual glaciers are usually strongly correlated and we calculate balanced-budget ELA from the regression equation linking the two. We then compare balanced-budget ELA with area-median and mid-range altitudes for the 94 glaciers for which we have all the necessary data. The different ELA estimates agree well enough (±82 to ±125 m) to describe geographical variations in ELA and for application of glacier–climate models to glacier inventory data. Mid-range and area-median altitudes are already available for tens of thousands of glaciers in the current WGI and should be evaluated in future inventories.

A prototype snow albedo algorithm has been developed for the Moderate Resolution Imaging Spectroradiometer (MODIS). It complements existing MODIS products by providing albedo measurements for areas mapped as snow on a global daily basis by MODIS. Cloud detection and atmospheric correction are accomplished using existing MODIS products. Models of the bidirectional reflectance of snow created using a discrete-ordinate radiative transfer (DISORT) model are used to correct for anisotropic scattering effects over non-forested surfaces. Initial algorithm validation is undertaken through comparisons with broadband albedo measurements made at the U.S. National Oceanic and Atmospheric Administration (NOAA) Surface Radiation Budget Network (SURFRAD) site in Fort Peck, MT. In situ SURFRAD albedo measurements are compared to daily MODIS snow albedo retrievals for the period 21–26 November 2000 created from five narrow-to-broadband albedo conversion schemes. The prototype MODIS algorithm produces reasonable broadband albedo estimates. Maximum daily differences between the five MODIS broadband albedo retrievals and in situ albedo are 15%. Daily differences between the best MODIS broadband estimate and the measured SURFRAD albedo are 1–8%. However, no single conversion scheme consistently provides the closest albedo estimate. Further validation and algorithm development using data from North America and Greenland is ongoing.

The consequences of global warming on land ice masses are difficult to assess in detail, as two-dimensional glacier inventory data are still missing for many remote regions of the world. As the largest future temperature increase is expected to occur at high latitudes, the glaciers and ice caps in the Arctic will be particularly susceptible to the expected warming. This study demonstrates the possibilities of space-borne glacier inventorying at a remote site on Cumberland Peninsula, a part of Baffin Island in Arctic Canada, thereby providing glacier inventory data for this region. Our approach combines Landsat ETM+ and Terra ASTER satellite data, an ASTER-derived digital elevation model (DEM) and Geographic Information System-based processing. We used thresholded ratio images from ETM+ bands 3 and 5 and ASTER bands 3 and 4 for glacier mapping. Manual delineation of Little Ice Age trimlines and moraines has been applied to calculate area changes for 225 glaciers, yielding an average area loss of 11%. A size distribution has been obtained for 770 glaciers that is very different from that for Alpine glaciers. Numerous three-dimensional glacier parameters were derived from the ASTER DEM for a subset of 340 glaciers. The amount of working time required for the processing has been tracked, and resulted in 5 min per glacier, or 7 years for all estimated 160 000 glaciers worldwide.

This research presents a precise evaluation of the recession of Akshiirak and Ala Archa glaciers, Tien Shan, central Asia, based on data of geodetic surveys from 1861–69, aerial photographs from 1943, 1963, 1977 and 1981, 1:25000 scale topographic maps and SRTM and ASTER data from 2000–03. The Akshiirak glacierized massif in the central Tien Shan contains 178 glaciers covering 371.6 km2, and the Ala Archa glacier basin in the northern Tien Shan contains 48 glaciers covering 36.31 km2. The Tien Shan glaciers retreated as much as 3 km from the 1860s to 2003. Area shrinkage of Akshiirak and Ala Archa was 4.2% and 5.1%, respectively, from 1943 to 1977, and 8.7% and 10.6%, respectively, from 1977 to 2003. The volume of the Akshiirak glaciers was reduced by 3.566 km3 from 1943 to 1977 and 6.145 km3 from 1977 to 2000. The total reduction of the Tien Shan glaciers is 14.2% during the last 60 years (1943–2003). The northern and central Tien Shan have not experienced a significant precipitation increase during the last 100 years, but they have experienced an increase in summer air temperatures, especially observable since the 1970s, which accelerated the recession of the Tien Shan glaciers.

The largest dischargers of ice in Greenland are glaciers that terminate in the ocean and melt in contact with sea water. Studies of ice-sheet/ocean interactions have mostly focused on melting beneath near-horizontal floating ice shelves. For tidewater glaciers, melting instead takes place along the vertical face of the calving front. Here we modify the Massachusetts Institute of Technology general circulation model (MITgcm) to include ice melting from a calving face with the freshwater outflow at the glacier grounding line. We use the model to predict melt rates and their sensitivity to ocean thermal forcing and to subglacial discharge. We find that melt rates increase with approximately the one-third power of the subglacial water flux, and increase linearly with ocean thermal forcing. Our simulations indicate that, consistent with limited field data, melting ceases when subglacial discharge is shut off, and reaches several meters per day when subglacial discharge is high in the summer. These results are a first step toward a more realistic representation of subglacial discharge and of ocean thermal forcing on the subaqueous melting of tidewater glaciers in a numerical ocean model. Our results illustrate that the ice-front melting process is both complex and strongly time-dependent.

Recent observational and modeling studies indicate that the Arctic sea-ice cover is undergoing significant climate-induced changes, affecting both its extent and thickness. The thickness or, more precisely, the mass balance of the ice cover is a key climate-change indicator since it is an integrator of both the surface heat budget and the ocean heat flux. Accordingly, efforts are underway to develop and deploy in situ observing systems which, when combined with satellite remote-sensing information and numerical models, can effectively monitor and attribute changes in the mass balance of the Arctic sea-ice cover. As part of this effort, we have developed an autonomous ice mass-balance buoy (IMB), which is equipped with sensors to measure snow accumulation and ablation, ice growth and melt, and internal ice temperature, plus a satellite transmitter. The IMB is unique in its ability to determine whether changes in the thickness of the ice cover occur at the top or bottom of the ice cover, and hence provide insight into the driving forces behind the change. Since 2000, IMBs have been deployed each spring from the North Pole Environmental Observatory and in several other areas, including a few in the Beaufort Sea and Central Basin. At this point, the collective time series is too short to draw significant and specific conclusions regarding interannual and regional variability in ice mass balance. Comparisons of available data indicate that ice surface ablation is greater in the Beaufort region (67–80 cm), relative to the North Pole (0–30 cm), consistent with a longer period of melt in the more southerly location. Ablation at the bottom of the ice (22 cm), maximum ice thickness (235 cm) and maximum snow depth (28 cm) were comparable in the two regions.

Temperature records from slations on the west roast of the Antarctic Peninsula show a very high level of interannual variability and, over the last 50 years, larger warming trends than are seen elsewhere in Antarctica. in this paper we investigate the role of atmospheric circulation variability and sea-ice extent variations in driving these changes. Owing to a lack of independent data, the reliability of Antarctic atmospheric analyses produced in the 1950s and 1960s cannot be readily established, but examination of the available data suggests that there has been an increase in the northerly component of the circulation over the Peninsula since the late 1950s. Few observations of sea-ice extent are available prior to 1973, but the limited data available indicate that the ice edge to the west of the Peninsula lay to the north of recently observed extremes during the very cold conditions prevailing in the late 1950s. The ultimate cause of the atmospheric-circulation changes remains to be determined and may lie outside the Antarctic region.

Several firn/ice cores were recovered from the Siberian Altai (Belukha plateau), central Tien Shan (Inilchek glacier) and the Tibetan Plateau (Zuoqiupu glacier, Bomi) from 1998 to 2003. The comparison analyses of stable-isotope/geochemistry records obtained from these firn/ice cores identified the physical links controlling the climate-related signals at the seasonal-scale variability. The core data related to physical stratigraphy, meteorology and synoptic atmospheric dynamics were the basis for calibration, validation and clustering of the relationships between the firn-/ice-core isotope/ geochemistry and snow accumulation, air temperature and precipitation origin. The mean annual accumulation (in water equivalent) was 106 gcm−2 a−1 at Inilchek glacier, 69 gcm−2 a−1 at Belukha and 196 g cm−2 a−1 at Zuoqiupu. The slopes in regression lines between the δ18O ice-core records and air temperature were found to be positive for the Tien Shan and Altai glaciers and negative for southeastern Tibet, where heavy amounts of isotopically depleted precipitation occur during summer monsoons. The technique of coupling synoptic climatology and meteorological data with δ18O and d-excess in firn-core records was developed to determine climate-related signals and to identify the origin of moisture. In Altai, two-thirds of accumulation from 1984 to 2001 was formed from oceanic precipitation, and the rest of the precipitation was recycled over Aral–Caspian sources. In the Tien Shan, 87% of snow accumulation forms by precipitation originating from the Aral–Caspian closed basin, the eastern Mediterranean and Black Seas, and 13% from the North Atlantic.

This paper reports changes in supraglacial debris cover and supra-/proglacial lake development associated with recent glacier retreat (1985–2000) in the central Caucasus Mountains, Russia. Satellite imagery (Landsat TM and ETM+) was used to map the surface area and supraglacial debris cover on six neighbouring glaciers in the Adylsu valley through a process of manual digitizing on a false-colour composite of bands 5, 4, 3 (red, green, blue). The distribution and surface area of supraglacial and proglacial lakes was digitized for a larger area, which extended to the whole Landsat scene. We also compare our satellite interpretations to field observations in the Adylsu valley. Supraglacial debris cover ranges from <5% to >25% on individual glaciers, but glacier retreat between 1985 and 2000 resulted in a 3–6% increase in the proportion of each glacier covered by debris. The only exception to this trend was a very small glacier where debris cover did not change significantly and remote mapping proved more difficult. The increase in debris cover is characterized by a progressive up-glacier migration, which we suggest is being driven by focused ablation (and therefore glacier thinning) at the up-glacier limit of the debris cover, resulting in the progressive exposure of englacial debris. Glacier retreat has also been accompanied by an increase in the number of proglacial and supraglacial lakes in our study area, from 16 in 1985 to 24 in 2000, representing a 57% increase in their cumulative surface area. These lakes appear to be impounded by relatively recently lateral and terminal moraines and by debris deposits on the surface of the glacier. The changes in glacier surface characteristics reported here are likely to exert a profound influence on glacier mass balance and their future response to climate change. They may also increase the likelihood of glacier-related hazards (lake outbursts, debris slides), and future monitoring is recommended.

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

Most of the central Himalayan glaciers have surface debris layers of variable thickness, which greatly affect the ablation rate. An attempt has been made to relate debris-cover thickness to glacier surface melting. Thirty stakes were used to calculate ablation for debris-covered and clean ice of Dokriani Glacier (7 km2) from 2009/10 to 2012/13. Our study revealed significant altitude-wise difference in the rate of clean and debris-covered ice melting. We found a high correlation (R2 = 0.92) between mean annual clean-ice ablation and altitude, and a very low correlation (R2 = 0.14) between debris-covered ice melting and altitude. Debris-covered ice ablation varies with variation in debris thickness from 1 to 40 cm; ablation was maximum under debris thicknesses of 1–6 cm and minimum under 40 cm. Even a small debris-cover thickness (1–2 cm) reduces ice melting as compared to that of clean ice on an annual basis. Overall, debris-covered ice ablation during the study period was observed to be 37% less than clean-ice ablation. Strong downwasting was also observed in the Dokriani Glacier ablation area, with average annual ablation of 1.82 m w.e. a–1 in a similar period. Our study suggests that a thinning glacier rapidly becomes debris-covered over the ablation area, reducing the rate of ice loss.

A distributed surface energy-balance study was performed to determine sub-debris ablation across a large part of Baltoro glacier, a wide debris-covered glacier in the Karakoram range, Pakistan. The study area is ~124km2. The study aimed primarily at analyzing the influence of debris thickness on the melt distribution. The spatial distribution of the physical and thermal characteristics of the debris was calculated from remote-sensing (ASTER image) and field data. Meteorological data from an automatic weather station at Urdukas (4022ma.s.l.), located adjacent to Baltoro glacier on a lateral moraine, were used to calculate the spatial distribution of energy available for melting during the period 1–15 July 2004. The model performance was evaluated by comparisons with field measurements for the same period. The model is reliable in predicting ablation over wide debris-covered areas. It underestimates melt rates over highly crevassed areas and water ponds with a high variability of the debris thickness distribution in the vicinity, and over areas with very low debris thickness (<0.03 m). We also examined the spatial distribution of the energy-balance components (global radiation and surface temperature) over the study area. The results allow us to quantify, for the study period, a meltwater production of 0.058 km3.

Ablation of glacier ice has been observed with artificial debris layers prepared with Ottawa sand (ASTM C-109) ranging from 0.01 to 0.1 m thick. Data on external variables observed during the experiments and determination of physical constants of the debris layers have allowed the testing of a proposed simple model. Theoretical predictions compare favourably with the observations. Discussion is extended to a proposal for a simple method by which ablation under a debris layer could be estimated even if the thermal conductivity or thermal resistance of the material were unknown.

Beginning in 1996, aerial photographs were taken for a new Austrian glacier inventory, resulting in digital elevation models (DEMs) and digital orthoimages. An earlier inventory of the Austrian glaciers containing the original aerial photographs and glacier maps and a manual evaluation of various glacier parameters as of 1969 has been re-evaluated at the present state of the art. The two inventories provide the basis for the comparison of glacier reactions over a period of 29 years. In general a reduction of glacier area is observed for almost all Austrian glaciers between 1969 and 1998. The overall reduction in ice-covered area is 17%. The glacier volume change calculated from the DEMs amounts to about 5 km3. This is almost 22% of the ice volume in Austria in 1969, estimated from a volume–area relation. Changes of individual glaciers, however, show a wide variability, depending on their size and physiographic setting.

This study presents a reconstruction of the mass balance (MB) of Chhota Shigri glacier, Western Himalaya, India, and discusses the regional climatic drivers responsible for its evolution since 1969. The MB is reconstructed by a temperature-index and an accumulation model using daily air-temperature and precipitation records from the nearest meteorological station, at Bhuntar Observatory. The only adjusted parameter is the altitudinal precipitation gradient. The model is calibrated against 10 years of annual altitudinal MB measurements between 2002 and 2012 and decadal cumulative MBs between 1988 and 2010. Three periods were distinguished in the MB series. Periods I (1969-85) and III (2001-12) show significant mass loss at MB rates of -0.36±0.36 and -0.57±0.36mw.e.a-1 respectively, whereas period II (1986-2000) exhibits steady-state conditions with average MBs of -0.01 ±0.36mw.e.a–1. The comparison among these three periods suggests that winter precipitation and summer temperature are almost equally important drivers controlling the MB pattern of Chhota Shigri glacier at decadal scale. The sensitivity of the modelled glacier-wide MB to temperature is -0.52 m w.e. a–1 °C–1 whereas the sensitivity to precipitation is calculated as 0.16mw.e.a-1 for a 10% change.

The importance of glaciers in mainland Norway for runoff is reflected in the extensive glacier measurement record. Mass balance has been measured for 42 glaciers. Length (or front-position) records exist for about 60 glaciers, and nearly half of these are presently measured. The mass-balance and front-position data have been analyzed with respect to spatial and temporal variations. The maritime glaciers with a large annual mass turnover have had a mass surplus between 1962 and 2000. In contrast, the continental glaciers with smaller summer and winter balances had a mass deficit over the same period. Since 2001 all monitored glaciers have had a marked mass deficit. The Norwegian glaciers have all retreated during the 20th century. However, both local and regional variations have been observed. Advances were recorded around 1910, around 1930, in the second half of the 1970s and around 1990. This last advance stopped in most glaciers at the turn of the century.

The evolution of surface topography of glaciers in the Swiss Alps is well documented with high-resolution aerial photographs repeatedly recorded since the 1960s and further back in time with topographic maps including elevation contour lines first surveyed in the mid-19th century. In order to quantify and interpret glacier changes in the Swiss Alps, time series of volume changes over the last 100–150 years have been collected. The available datasets provide a detailed spatial resolution for the retreat period since the end of the Little Ice Age. The spatial distribution as well as temporal variations of the thickness change were analyzed. A significant ice loss since the end of the 19th century was observed in the ablation area, while the changes in the accumulation area were small. We found moderate negative secular rates until the 1960s, followed by steady to positive rates for about two decades and strong ice loss starting in the 1980s which has lasted until the present. An evaluation of 19 glaciers revealed a total ice volume loss of about 13km3 since the 1870s, of which 8.7 km3 occurred since the 1920s and 3.5 km3 since 1980. Decadal mean net balance rates for the periods 1920–60, 1960–80 and 1980–present are –0.29, –0.03 and –0.53ma–1w.e., respectively.

Determination of the spatial snow-depth distribution is important in potential avalanche-starting zones, both for avalanche prediction and for the dimensioning of permanent protection measures. Knowledge of the spatial distribution of snow is needed in order to validate snow depths computed from snowpack and snowdrift models. The inaccessibility of alpine terrain and the acute danger of avalanches complicate snow-depth measurements (e.g. when probes are used), so the possibility of measuring the snowpack using terrestrial laser scanning (TLS) was tested. The results obtained were compared to those of tachymetry and manual snow probing. Laser measurements were taken using the long-range laser profile measuring system Riegl LPM-i800HA. The wavelength used by the laser was 0.9 μm (near-infrared). The accuracy was typically within 30 mm. The highest point resolution was 30 mm when measured from a distance of 100 m. Tachymetry measurements were carried out using Leica TCRP1201 systems. Snowpack depths measured by the tachymeter were also used. The datasets captured by tachymetry were used as reference models to compare the three different methods (TLS, tachymetry and snow probing). This is the first time that the accuracy of TLS systems in snowy and alpine weather conditions has been quantified. The relative accuracy between the three measurement methods is bounded by a maximum offset of ±8 cm. Between TLS and the tachymeter the standard deviation is 1σ = 2 cm, and between manual probing and TLS it is up to 1σ = 10 cm, for maximum distances for the TLS and tachymeter of 300 m.

In order to extend diagnoses of recent sea-ice variations beyond the past few decades, a century-scale digital dataset of Arctic sea-ice coverage has been compiled. For recent decades, the compilation utilizes satellite-derived hemispheric datasets. Regional datasets based primarily on ship reports and aerial reconnaissance are the primary inputs for the earlier part of the 20th century. While the various datasets contain some discrepancies, they capture the same general variations during their period of overlap. The outstanding feature of the time series of total hemispheric ice extent is a decrease that has accelerated during the past several decades. The decrease is greatest in summer and weakest in winter, contrary to the seasonality of the greenhouse changes projected by most global climate models. The primary spatial modes of sea-ice variability diagnosed in terms of empirical orthogonal functions, also show a strong seasonality. The first winter mode is dominated by an opposition of anomalies in the western and eastern North Atlantic, corresponding to the well-documented North Atlantic Oscillation. The primary summer mode depicts an anomaly of the same sign over nearly the entire Arctic and captures the recent trend of sea-ice coverage.

The degree-day factor (DDF) is an important parameter for the degree-day model, which is a widely used method for ice- and snowmelt computation. Spatial variations of the DDF greatly affect the accuracy of snow- and ice-melt modelling. This study analyzes the spatial variability of DDFs obtained from observed glaciers in different regions of western China. The results clearly show that the DDF for a single glacier is subject to significant small-scale variations, and the factor for maritime glaciers is higher than that for subcontinental and extremely continental glaciers. In western China the factors increase gradually from northwest to southeast. In general, the regional patterns of DDFs are detectable on the glaciers due to the unique climatic environment and heat budget of the Tibetan Plateau and surrounding regions. Low DDFs can be expected for cold-dry areas, whereas high DDFs can be expected for warm-wet areas in western China. Depending on spatial variation of the characteristics of DDFs and the meteorological data, we can provide gridded degree-day models for non-monitored glaciers to reconstruct gridded historical glacier mass-balance series in western China.