Supraglacial meltwater lakes trigger ice-shelf break-up and modulate seasonal ice-sheet flow, and are thus agents by which warming is transmitted to the Antarctic and Greenland ice sheets. To characterize supraglacial lake variability we perform a comparative analysis of lake geometry and depth in two distinct regions, one on the pre-collapse (2002) Larsen B ice shelf, Antarctica, and the other in the ablation zone of Paakitsoq, a land-terminating region of the Greenland ice sheet. Compared to Paakitsoq, lakes on the Larsen B ice shelf cover a greater proportion of surface area (5.3% cf. 1%), but are shallower and more uniform in area. Other aspects of lake geometry (e.g. eccentricity, degree of convexity (solidity) and orientation) are relatively similar between the two regions. We attribute the notable difference in lake density and depth between ice-shelf and grounded ice to the fact that ice shelves have flatter surfaces and less distinct drainage basins. Ice shelves also possess more stimuli to small-scale, localized surface elevation variability, due to the various structural features that yield small variations in thickness and which float at different levels by Archimedes’ principle.

Monitoring the glacier mass balance of summer-accumulation-type Himalayan glaciers is critical to not only assess the impact of climate change on the volume of such glaciers but also predict the downstream water availability and the global sea-level change in future. To better understand the change in meteorological parameters related to glacier mass balance and runoff in a glacierized basin and to assess the highly heterogeneous glacier responses to climate change in the Nepal Himalaya and nearby ranges, the Cryosphere Monitoring Project (CMP) carries out meteorological observations in Langtang Valley and mass-balance measurements on Yala Glacier, a debris-free glacier in the same valley. A negative annual mass balance of –0.89m w.e. and the rising equilibrium-line altitude of Yala Glacier indicate a continuation of a secular trend toward more negative mass balances. Lower temperature lapse rate during the monsoon, the effect of convective precipitation associated with mesoscale thermal circulation in the local precipitation and the occurrence of distinct diurnal cycles of temperature and precipitation at different stations in the valley are other conclusions of this comprehensive scientific study initiated by CMP which aims to yield multi-year glaciological, hydrological and meteorological observations in the glacierized Langtang River basin.

Simulations of glacier evolution are needed to assess future changes in the runoff regime of mountain catchments. A simplified parameterized model is applied here to simulate future thickness changes and glacier retreat of Mer de Glace, French Alps. A normalized thickness change function describing the spatial distribution of surface-elevation changes as a function of elevation has been determined. The model reveals that under present climatic conditions Mer de Glace will continue to shrink dramatically in the coming decades, retreating by 1200 m between now and 2040. The method has certain limitations related to the uncertainties of the normalized function based on thickness change data. An error of 10% in the normalized function leads to uncertainties of 46%, 30% and 18% in Mer de Glace front, surface area and glacier-wide mass-balance changes respectively in 2040. Because the difference of the normalized function largely exceeds 10% from one glacier to another, even within a given glacier size class and elevation range, it would be very risky to extrapolate the normalized function to unmeasured glaciers. Consequently, the method is applicable only on glaciers where past surface elevation changes are well constrained.

To study the glaciological processes controlling the mass budget of Greenland’s peripheral glaciers and ice caps, field measurements were carried out on Qaanaaq ice cap, a 20 km long ice cap in northwestern Greenland. In the summer of 2012, we measured surface melt rate, ice flow velocity and ice thickness along a survey route spanning the ice margin (200m a.s.l.) to the ice-cap summit (1110m a.s.l.). Melt rates in the ablation area were clearly influenced by dark materials covering the ice surface, where degree-day factors varied from 5.44 mm w.e. K–1 d–1 on a clean surface to 8.26 mm w.e. K–1 d–1 in the dark regions. Ice velocity showed diurnal variations, indicating the presence of surface-meltwater induced basal sliding. Mean ice thickness along the survey route was 120 m, with a maximum thickness of 165 m. Ice velocity and temperature fields were computed using a thermomechanically coupled numerical glacier model. Modelled ice temperature, obtained by imposing estimated annual mean air temperature as the surface boundary condition, was substantially lower than implied by the observed ice velocity. This result suggests that the ice dynamics and thermodynamics of the ice cap are significantly influenced by heat transfer from meltwater and changing ice geometry.

Kyagar glacier is located in the Chinese Karakoram mountains. The glacier tongue entirely blocks the riverbed in the upper Shaksgam valley and impounds a glacial lake, which was the source of several violent and disastrous glacial lake outburst floods (GLOFs). A GLOF early warning system was implemented between 2011 and 2013. We present an integrative analysis of the hazard potential of Kyagar lake, taking into account the ice flow dynamics of Kyagar glacier as well as the recent surface mass-balance response to climate change. Comparison of two high-resolution digital elevation models (DEMs) for the ice dam shows surface lowering rates of >5ma– 1 between 2002 and 2011, leading to a significant reduction in the maximum potential lake volume. However, two DEMs covering the entire glacier for the period 2000–10 indicate mass gains in its central part, and flow speed measurements show an acceleration in this region. This pattern of local ice-thickness changes combined with varying ice flow velocities is typical for surge-type glaciers. The velocity of the glacier surface and of the ice dam between 2011 and 2012 are analyzed at high temporal and spatial resolution, based on feature tracking of synthetic aperture radar (SAR) images.

We present a study of seasonal and interannual ice velocity changes at Polar Record Glacier, East Antarctica, using ERS-1/2, Envisat and PALSAR data with D-InSAR and intensity tracking. Ice flow showed seasonal variations at the front of the glacier tongue. Velocities in winter were 19% less than velocities during summer. No significant interannual changes were detected. Ice velocities in the grounding zone and grounded glacier did not show clear seasonal or interannual changes. The distribution of the seasonal variations suggests that the cause for the changes should be localized. Possible causes are seasonal sea-ice changes and iceberg blocking. Satellite images show that the sea ice surrounding Polar Record Glacier undergoes seasonal changes. Frozen sea ice in winter slowed the huge iceberg, and provided increased resistance to the glacier flow. The interaction between the glacier tongue, iceberg and sea ice significantly influences their flow pattern.

The Tianger area in the Tien Shan preserves many glacial landforms to which several methods of estimating ancient equilibrium-line altitudes (ELAs) can be applied. Thus it provides the opportunity to compare different methods and evaluate the correlation between the methods and their errors. This comparison provides guidance in the application of these methods in other mid-latitude alpine glaciated areas. In this paper we calculate characteristic ratios relating to the ELAs, including accumulation–area ratio (AAR), toe–headwall altitude ratio (THAR), terminus–summit altitude method (TSAM) and Höfer value, and their correlations with each other and the geomorphological characteristics. The AAR of a glacier is small when the slope of the bedrock beneath it is small and the area is large, and the AAR of a glacier with such characteristics is better correlated with its THAR, TSAM and Höfer value. Paleoclimate information derived from these characteristic ratios relating to ELAs is discussed. The glaciers and the glacial climate of the early and middle stage of the last glacial are more continental than those of the Last Glacial Maximum (LGM). Compared with those of the LGM on the Qinghai–Tibetan Plateau and surrounding mountains, the ELA characteristics of the early and middle stage of the last glacial in the Tianger area show more continental features.

We use remote-sensing and GIS technologies to monitor glacier changes in the Koshi River basin, central Himalaya. The results indicate that in 2009 there were 2061 glaciers in this region, with a total area of 3225 ±90.3 km2. This glacier population is divided into 1290 glaciers, with a total area of 1961 ±54.9 km2, on the north side of the Himalaya (NSH), and 771 glaciers, with a total area of 1264 ± 35.4 km2, on the south side of the Himalaya (SSH). From 1976 to 2009, glacier area in the basin decreased by about 19±5.6% (0.59±0.17%a–1). Glacier reduction was slightly faster on SSH (20.3 ±5.6%) than on NSH (18.8±5.6%). The maximum contribution to glacier area loss came from glaciers within the 1-5 km2 area interval, which accounted for 32% of total area loss between 1976 and 2009. The number of glaciers in the Koshi River catchment decreased by 145 between 1976 and 2009. Glacier area on SSH decreased at a rate of 6.2 ±3.2% (0.68 ±0.36% a–1), faster than on NSH, where the rate was 2.5 ±3.2% (0.27±0.36% a–1) during 2000-09. Based on records from Tingri weather station, we infer that temperature increase and precipitation decrease were the main causes of glacier thinning and retreat during the 1976-2000 period. Glacier retreat during the 2000-09 period appears to be controlled by temperature increase, since precipitation increase over this period did not offset ice losses to surface melting.

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.

We combine radio-echo sounding ice thickness data from the BEDMAP Project database and the PCMEGA (Prince Charles Mountains Expedition of Germany and Australia) dataset to generate a new ice thickness grid for the southern limit region of the Amery Ice Shelf, East Antarctica. We then reassess the mass balance of the central portion of the Lambert-Amery system, incorporating flow information derived from synthetic aperture radar interferometry (InSAR) and a modeled surface mass-balance dataset based on regional atmospheric modeling. Our analysis reveals that Mellor and Fisher Glaciers are approximately in balance to the level of our measurement uncertainty, while Lambert Glacier has a positive imbalance of 4.2 ±2.3 Gta1. The mass budget for the whole Lambert Glacier basin is approximately in balance, and the average basal melt rate in the downstream section of the ice shelf is 5.1 ± 3.0 m a-1. Our results differ substantially from other recent estimates using hydrostatically derived ice thickness data.

At both corners of the now stagnant Kamb Ice Stream (KIS, West Antarctica) outlet, shear margins of different ages confine wedge-shaped areas with relatively flat, smooth surfaces that stagnated before the main trunk of the ice stream. We identify these features as lakes or past lakes, and consider scenarios for their development in a regional history of ongoing adjustment to grounding-line position. We focus here on the centuries leading up to the recent stagnation of KIS, a time when its grounding line appears to have advanced >100km from an earlier upstream location. Starting from stagnation of Crary Ice Rise and changes in the grounding zone of Whillans Ice Stream, we trace feedbacks associated with local thickening, ice grounding and thickness transients that both advance the grounding line and leave remnant lakes in their wake. These lakes in turn promote the development of secondary margins that may appear as ‘margin jumps’ in the ice record.

δ18O measurements based on systematic sampling and isotopic modeling have been adopted to study the controls of stable isotopes in lake water in two lake basins (lakes Yamdrok-tso and Puma Yum-tso) at two different elevations on the southern Tibetan Plateau. Temporally, δ18O values in precipitation and lake water display a seasonal fluctuation in both lakes. Spatially, δ18O values in the two lake basins increase by 10% from the termini of glaciers to the lake shores, by ∽1% from the lake shores to the lake center and by 0.4% from the water surface to depth in these lakes. The clear annual δ18O variations indicate that lake water mixes sufficiently in a short time. Model results show that glacial meltwater and surface lake-water temperature are not the dominant factors in the balance process of stable isotopes in lake water. Equilibrium δ18O values decrease by 0.8% for Yamdrok-tso lake and 0.6% for Puma Yum-tso lake when glacial meltwater contributions to these lakes shrink by 60%. δ18O ratios increase rapidly during the initial stages and take a longer time to approach the equilibrium value.

In 2005, two ice cores with lengths of 58.7 and 57.6 m respectively to bedrock were recovered from the Miaoergou flat-topped glacier (43°03′19″ N, 94°19′21″ E; 4512 m a.s.l.), eastern Tien Shan. 210Pb dating of one of the ice cores (57.6 m) was performed, and an age of AD 1851 ± 6 at a depth of 35.2 mw.e. was determined. For the period AD 1851-2005, a mean annual net accumulation of 229 ± 7 mm w.e. a–1 was calculated. At the nearby oasis city of Hami (~80 km from the Miaoergou flat-topped glacier) the annual precipitation rate is 38 mm w.e. a–1, hence glacial meltwater is a major water supply for local residents. The surface activity concentration of 210Pbex was found to be ~400 mBq kg–1, which is higher than observed at other continental sites such as Belukha, Russia, and Tsambagarav, Mongolia, which have surface activity concentrations of 280 mBq kg–1. The 210Pb dating agrees well with the chronological sequence deduced from the annual-layer counting resulting from the seasonalities of d18O and trace metals for the period AD 1953-2005, and ^-activity horizons resulting from atmospheric nuclear testing during the period AD 1962-63. We conclude that 210Pb analysis is a suitable method for obtaining a continuous dating of the Miaoergou ice core for ~160 years, which can also be applied to other ice cores recovered from the mountains of western China.

Glacier surface melting can be described using energy-balance models. We conducted a surface energy budget experiment to quantify surface energy fluxes and to identify factors affecting glacial melt in the ablation zone of Laohugou glacier No. 12, western Qilian mountains. The surface energy budget was calculated based on data from an automatic weather station, and turbulent fluxes calculated using the bulk-aerodynamic approach were corrected using measurements from an eddy-covariance system. Simulated mass balances were validated by stake observations. Net shortwave radiation was the primary component of the surface energy balance (126Wm–2), followed by sensible heat flux. Net longwave radiation (–45Wm–2) and latent heat flux (–12.8 Wm–2) represented heat sinks. The bulk-aerodynamic method underestimated sensible and latent heat fluxes by 3.4 and 1.2 W m–2, respectively. The simulated total mass balance of –1703mmw.e. exceeded the observed total by 90 mm w.e. Daily positive accumulated temperature and low albedo were the main factors accelerating glacier melt. An uncertainty assessment showed that mass balance was very sensitive to albedo and varied by 36% when albedo changed by 0.1.

The Asian–Pacific Oscillation (APO) is a climate pattern related to the thermal differences between the Asian continent and the north Pacific. We present a 1000 year record of the major ions in a 108.83 m ice core from east Rongbuk glacier (28°01′ N, 86°58′ E; 6518ma.s.l.) on the northeast slope of Qomolangma (Mount Everest), and discuss its relationship with a 993 year time series of the summer APO index that has been reconstructed from Beijing’s summer temperature (recorded by a stalagmite) and the Pacific Decadal Oscillation (PDO) index. Empirical orthogonal function (EOF) analysis shows that crustal major ions (Mg2+, Ca2+, SO42– and NO3–) are highly loaded on EOF1, whereas marine major ions (Cl– and Na+) are highly loaded on EOF2. Higher EOF1 is associated with lower upper-troposphere temperature (UTT) over Mongolia, corresponding to a lower APO index and higher surface pressure. Higher EOF2 is associated with higher UTT over Mongolia, corresponding to a higher APO index and lower surface pressure. The higher burden of major ions and higher summer APO index during the period AD 1000–1323 may be due to drought and warm climate of the major-ion source regions. From AD1323 to 1900, EOF1 and EOF2 do not show a consistent correlation with summer APO index, indicating the complex mechanisms of ion transport over the southern Tibetan Plateau (TP) during this period. After AD 1900, the summer APO index is correlated negatively with EOF1 and positively with EOF2, indicating that ion transport over the southern TP during this period is influenced significantly by the APO. Our examinations of the PDO index and major-ion record show that higher PDO corresponds to higher EOF1 and lower EOF2. This suggests that the correlation between the recorded major-ion concentration and the APO index originates in the seesaw relationship in temperature between the Asian continent and the north Pacific.

Blue-ice areas (BIAs) and their geographical distribution in Antarctica were mapped using Landsat-7 ETM+ images with 15 m spatial resolution obtained during the 1999–2003 austral summers and covering the area north of 82.5° S, and a snow grain-size image of the MODIS-based Mosaic of Antarctica (MOA) dataset with 125 m grid spacing acquired during the 2003/04 austral summer from 82.5°S to the South Pole. A map of BIAs was created with algorithms of thresholds based on band ratio and reflectance for ETM+ data and thresholds based on snow grain size for the MOA dataset. The underlying principle is that blue ice can be separated from snow or rock by their spectral discrepancies and by different grain sizes of snow and ice. We estimate the total area of BIAs in Antarctica during the data acquisition period is 234 549 km2, or 1.67% of the area of the continent. Blue ice is scattered widely over the continent but is generally located in coastal or mountainous regions. The BIA dataset presented in this study is the first map covering the entire Antarctic continent sourced solely from ETM+ and MODIS data. This dataset can potentially benefit other studies in glaciology, meteorology, climatology and paleoclimate, meteorite collection and airstrip site selection.

Alberta’s Bow River has its headwaters in the glaciated eastern slopes of the Canadian Rockies and is a major source of water in southern Alberta. Glacial retreat, declining snowpacks and increased water demand are all expected in the coming century, yet there are relatively few studies focusing on quantifying glacial meltwater in the Bow River. We develop a new radiation-temperature melt model for modelling distributed glacier mass balance and runoff in the Bow River basin. The model reflects physical processes through the incorporation of near-surface air temperature and absorbed radiation, while avoiding problems of collinearity through the use of a radiation-decorrelated temperature index. The model is calibrated at Haig Glacier in the southern portion of the basin and validated at Haig and Peyto Glaciers. Application of the model to the entire Bow River basin for 2000-09 shows glacier ice melt is equivalent to 3% of annual discharge in Calgary on average. Modelled ice melt in August is equal to 8-20% of the August Bow River discharge in Calgary. This emphasizes the importance of glacier runoff to late-summer streamflow in the region, particularly in warm, dry years.

The ice cover on the Qinghai–Tibetan Plateau plays an important role in the environmental and ecological systems. We analyze the in situ measurements of ice growth and examine the thermal diffusivity of thermokarst lake ice in the Beiluhe basin. We evaluate numerically the change of thermal diffusivity of thermokarst lake ice with changing ice temperature using an optimal control model. In a higher ice temperature regime (–3 to 0°C), the thermal diffusivity of thermokarst lake ice decreases exponentially with increasing ice temperature, and approaches the thermal diffusivity value of fresh water near the freezing-point temperature. In a lower ice temperature regime (–15 to –3°C), the thermal diffusivity increases slowly with decreasing ice temperature.