Recent climate changes have had a significant impact on the high-mountain glacial environment. Rapid melting of glaciers has resulted in the formation and expansion of moraine-dammed lakes, creating a potential danger from glacial lake outburst floods (GLOFs). Most lakes have formed during the second half of the 20th century. Glaciers in the Mount Everest (Sagamartha) region, Nepal, are retreating at an average rate of 10–59 ma–1. From 1976 to 2000, Lumding and Imja Glaciers retreated 42 and 34 ma–1, respectively, a rate that increased to 74 ma–1 for both glaciers from 2000 to 2007. During the past decade, Himalayan glaciers have generally been shrinking and retreating faster while moraine-dammed lakes have been proliferating. Although the number of lakes above 3500 m a.s.l. has decreased, the overall area of moraine-dammed lakes is increasing. Understanding the behaviour of glaciers and glacial lakes is a vital aspect of GLOF disaster management.

Global warming is causing an apparent rapid retreat of many glaciers worldwide. In addition to mass-balance investigation, the determination and monitoring of total glacial ice volume and ice-thickness distribution are important parameters for understanding the interactions between climate and the complex glacier system. Because of spatially irregular and sparse datasets, scaling of volume and ice-thickness distribution is often a challenging problem. This study focuses on two small (<2 km2) temperate glaciers in the Hohe Tauern (Eastern Alps) region of central Austria. The period 2003–04 saw the first use of ground-penetrating radar (GPR) to determine the total ice volume and ice-thickness distribution of the two glaciers. A centre frequency of 20 MHz was used in point measuring mode. Despite variable data quality, bedrock reflections up to depths of >100m were identified in the data. The acquired GPR data are irregularly distributed and the spatial density is too low to calculate reasonable bedrock topography with standard interpolation approaches. Thus one main focus of this study was to develop an appropriate interpolation technique. Eventually, kriging technique and a glacial mechanically based interpolation parameter were used. Mean calculated ice thicknesses for the two investigated glaciers are 40–50 m, with a maximum of 150–165 m. No direct validation data are available, so different considerations support the computed bedrock topography.

The influence of meteorological variables on snow/ice melting has been analyzed for two very contrasting months, in summer 2006, on Glacier de Saint-Sorlin, French Alps. July 2006 was the warmest July since 1950, and August 2006 was the coldest August since 1979. The total energy available for melting was just over half as much in August as in July, due to a sharp decrease in net shortwave radiation and in turbulent flux. This decrease of net shortwave radiation was mainly controlled by a strong increase in albedo responsible for an increase of reflected shortwave radiation, as well as by a reduction in incident shortwave radiation. During the two months, net longwave radiation remained almost unchanged. The mass balance computed from energy-balance modelling or with a degree-day approach was in good agreement with measured mass balance. Differences were attributed to space and time surface aspect variations which mainly controlled the observed mass balance.

Glacier surges tend to be initiated in relatively small regions, then propagate down-glacier, up-glacier and/or across-glacier. The processes controlling patterns and rates of surge propagation, however, are incompletely understood. In this paper, we focus on patterns of surge propagation in two confluent glaciers in Svalbard, and examine possible causes. One of these glaciers, Bakaninbreen, surged in 1985–95. The surge propagated ∽7 km down-glacier, but did not cross the medial moraine onto the other glacier, Paulabreen. When Paulabreen surged between 2003 and 2005, the surge wave travelled several km down-glacier, but its lateral boundary stayed very close to the medial moraine. The confluent glaciers formerly extended into a fjord, and bathymetric mapping and historical observations show that an active subglacial conduit has existed between Bakaninbreen and Paulabreen since at least the early 20th century. The existence of a persistent subglacial conduit below the medial moraine was confirmed when we entered and mapped a Nye channel at the confluence of Bakaninbreen and Paulabreen. We argue that the conduit acts as a barrier to surge propagation. If pressurized water below one branch of the glacier system reaches the conduit, water can be readily evacuated, preventing its propagation into the other branch.

A puzzling phenomenon of ice-stream flow is the stick–slip motion displayed by Whillans Ice Stream (WIS), West Antarctica. In this study we test the hypothesis that the WIS stick–slip motion has features similar to those of other known stick–slip systems, and thus might be of the same origin. To do so, we adapt a simple mechanical model widely used in seismology to study classic stick–slip behavior observed in tectonic faults, in which the difference between static and dynamic friction allows for the generation and spatial propagation of abrupt slip events. We show how spatial variability in friction properties, as well as a periodic forcing intended to mimic the effect of tides, can reproduce the observed duration and periodicity of stick–slip motion in an ice stream. An intriguing aspect of the association of WIS with mechanical stick–slip oscillators is that the onset of stick–slip cycling from a condition of permanent slip appears to be associated with the reduction in overall speed of WIS. If this association is true, then stick–slip behavior of WIS is a transitional phase of behavior associated with the ice stream's recent deceleration.

During three summer field seasons (2003, 2005 and 2006) we carried out radio-echo sounding measurements with a 5MHz (central frequency) ice-penetrating radar on McCall Glacier, Arctic Alaska, USA, along the central flowline and 17 cross-profiles. Two-way travel time was, after migration, converted to ice thickness, which, in combination with a recent digital elevation model of the surface of the glaciated area, resulted in a detailed map of the bed topography. This reveals a complex basal topography in the confluence area of the different glacial cirques. The pattern of subglacial water flow following the hydraulic potential gradient was calculated for the whole glacier area and shows a confluence of subglacial water downstream from the confluence of the glacier cirques. From the ice-thickness map the total ice volume was estimated as slightly less than 0.5 km3. Bed reflection power (BRP) was determined for the glacier after correction for ice-thickness dependence. Results reveal a clear relationship between the BRP pattern and basal sliding anomalies along the central flowline.

Numerous glaciological data have been obtained from measurements carried out on Glacier d’Argentière, Mont Blanc area, France, since the beginning of the 20th century. Moreover, data on annual mass balance, ice-flow velocity, thickness variation and length fluctuation have been obtained from yearly measurements performed since 1975. This dataset provides an excellent opportunity to analyze the relationships between surface mass balance and dynamic response over time periods during which net mass balance changed from positive to negative. Following a positive specific-net-balance period between 1960 and 1981, the ablation zone experienced a large increase in thickness and ice-flow velocities. Conversely, the highly negative specific-net-balance period since 1982 has led to strong thinning, deceleration and retreat of the tongue. The response of these observed dynamics to surface mass balance is analyzed from ice-flux calculations performed on three transverse cross-sections. Our results reveal that the ice fluxes are largely accommodated by ice-flow velocities. Velocity fluctuations are synchronous over the entire area studied. In the largest part of the glacier, no compressing/extending flow change has been observed over the last 30 years and thickness changes are solely driven by surface mass-balance changes. However, on the tongue of the glacier, thickness changes do not depend on surface mass balance but are mainly driven by changes in the longitudinal strain rate.

The ice cover of the Nianchu river basin, southern Tibetan Plateau, was mapped for 2005, using a SPOT5 scene, and for 1990 and 2000 from Landsat TM/ETM. Digital elevation models (DEMs) were generated from 1 : 50 000 scale topographical maps. The results show that in 2005 there were 136 glaciers in this basin, with a total area of 224 km2. Of these, 37 glaciers had an area >1 km2 and 10 were larger than 5 km2; the average snout altitude was 5608m a.s.l. A comparison of outlines from the last 15 years shows that most glaciers have decreased in size; none have advanced. From 1990 to 2005, Xiaquepu glacier No. 56 and Shimozongpu glacier No. 38 retreated 310 and 560 m, respectively. The mean reduction in glacier surface area was 5% (10 a)–1 while the area of glacial lakes expanded by 10%; nine new lakes formed in this basin over the 15 year period. Because air temperatures here have been increasing, while precipitation has remained steady, glacier retreat is considered to be related to rising temperature.

We have used ground-penetrating radar (GPR) to observe englacial structural control upon the development of an esker formed during a high-magnitude outburst flood (jökulhlaup). The surge-type Skeiðarárjökull, an outlet glacier of the Vatnajökull ice cap, Iceland, is a frequent source of jökulhlaups. The rising-stage waters of the November 1996 jökulhlaup travelled through a dense network of interconnected fractures that perforated the margin of the glacier. Subsequent discharge focused upon a small number of conduit outlets. Recent ice-marginal retreat has exposed a large englacial esker associated with one of these outlets. We investigated structural controls on esker genesis in April 2006, by collecting >2.5km of GPR profiles on the glacier surface up-glacier of where the esker ridge has been exposed by meltout. In lines closest to the exposed esker ridge, we interpret areas of englacial horizons up to ~30m wide and ~10–15m high as an up-glacier continuation of the esker sediments. High-amplitude, dipping horizons define the base of esker materials across many lines. Similar dipping surfaces deeper in the profiles suggest that: (1) the dipping surfaces beneath the esker are englacial tephera bands; (2) floodwaters were initially discharged along structurally controlled englacial surfaces (tephra bands); (3) the rapid increase in discharge resulted in hydrofracturing; (4) establishment of preferential flow paths resulted in conduit development along the tephra bands due to localized excavation of surrounding glacier ice; and (5) sedimentation took place within the new accommodation space to form the englacial structure melting out to produce the esker.

A new glacier inventory of Iran, compiled according to GLIMS guidelines through the use of photogrammetry and remote sensing supported by fieldwork, provides the first comprehensive study of its mountain glaciers. The glaciers are found in five main areas: two in the higher elevations of the Alborz mountain range (Damavand and Takhte–Soleiman regions), two on the Zardkuh and Oshtorankuh mountain chain in the Zagros mountain range and one in the Sabalan Mountains in northwest Iran. Several important glacier attributes, including minimum and maximum height of ice, area and maximum length and width, together with glacier extent, were successfully extracted using aerial and satellite imagery. Thereafter a comprehensive glacier database was established in a GIS environment.

Water stored in and released from supraglacial lakes is an important factor when considering the seasonal and long-term evolution of the Greenland ice sheet. Here we use a radiative transfer model to estimate changes in the depth and volume of a supraglacial lake on the surface of Jakobshavn Isbræ, West Greenland, between 2002 and 2005. When compared to estimates of the lake depth determined from airborne lidar observations, we estimate that the root-mean-square departure of the modelled lake depths was 0.3 m during cloud-free conditions. The maximum lake area, depth and volume were 3.4 km2, 9.6 ±1.0 m and (18.6±3.7)×106 m3, respectively. When sequenced according to the number of positive degree-days (PDDs) accumulated prior to each image, we observe that the lake volume evolves in three distinct phases. At the start of the melting season, the rate of filling is slow; after approximately 80 PDDs the rate of filling increases by a factor ∽3, and after approximately 125 PDDs the lake drains rapidly. We estimate that the lake drains at a minimum rate of (2.66±0.53)×106 m3 d–1 over a 6 day period.

Glacier mass-balance models that employ the degree-day method of melt modeling are most commonly driven by surface air temperatures that have been downscaled over the area of interest, using digital elevation models and assuming a constant free air lapse rate that is often taken to be the moist adiabatic lapse rate (MALR: –6.5°Ckm–1). Air-temperature lapse rates measured over melting glacier surface are, however, consistently less steep than free air values and have been shown to vary systematically with lower-tropospheric temperatures. In this study, the implications of including a variable near-surface lapse rate in a 26 year (1980–2006) degree-day model simulation of the surface mass balance of Devon Ice Cap, Nunavut, Canada, are examined and compared with estimates derived from surface air temperatures downscaled using a constant near-surface lapse rate equal to the measured summer mean (–4.9°Ckm–1) and the MALR. Our results show that degree-day models are highly sensitive to the choice of lapse rate. When compared with 23 years of surface mass-balance measurements from the northwest sector of the ice cap, model estimates are significantly better when surface air temperatures are downscaled using a modeled daily lapse rate rather than a constant lapse equal to either the summer mean or the MALR.

Measurements of basal water pressure from 15 boreholes located at both local (tens of meters) and regional (kilometers) length scales were used to elucidate the pressure/sliding relationship during an autumn rapid motion event on Bench Glacier, Alaska, USA. The 8 day event had two distinct phases, each with a ten-fold speed-up with respect to winter velocity. The water pressure in all 15 boreholes varied synchronously during the speed-up. The first phase of rapid sliding began after a peak in basal water pressure and continued while the pressure was elevated and stable, or decreasing. The second phase of rapid sliding occurred when the basal water pressure was low but increasing, and terminated before the pressure peaked. Pressure and velocity do not appear unrelated, but the pressure/sliding relationship was not consistently linked to increasing, decreasing or a critical water pressure. The pressure variations and sliding accelerations are a response to a warm rainstorm, although equally large input events occurred in weeks prior with no apparent response. Drainage system evolution therefore appears to play a key role in both the acceleration and the pressure/velocity relationship. Basal cavity dynamics are likely responsible for three episodes of reverse (up-valley) motion observed after enhanced sliding.

According to the division into subareas of water-resource distribution in China, and based on the Glacier Inventory of China (GIC), China’s total glacier water storage is 5040.2 x 109m3, 33.0% of which is distributed in the southwest drainage basins and 64.1% in the northwest inland drainage basins, forming enormous solid reservoirs with plentiful freshwater storage. Glacier change in China is estimated for the periods from the Little Ice Age (LIA) maximum to about 1960 and from 1960 to 1995. The relative glacierized area loss is 23% and 8.9% respectively for maritime glaciers and 15% and 4.9% respectively for continental glaciers. The normal annual glacier meltwater runoff is estimated at 60 465 x 106m3 by the climate parameter temperature-index method, 38.7% of which is distributed in the northwest inland drainage basins, and at 61 574 x 106m3 by the glacier system temperature-index method, 41.5% of which is distributed in the northwest inland drainage basins. Simulation of glacier meltwater runoff under the temperature change ratio of 0.03 K a–1 by the glacier system model in west China between 1980 and 2000 indicates that the total glacier meltwater runoff increment is 10.8% overall, 14.3% in the inland drainage basins of northwest China and 9.0% in the outflow drainage basins of southwest China.

In this study, long-term series of winter mass balances from two neighbouring glaciers in the southern Oetztal Alps, Austria, i.e. Hintereisferner and Vernagtferner, are analyzed with respect to the methods used in their determination. For this purpose, (1) some basic data of field surveys are presented, (2) the influence of different temporal systems is discussed, and (3) the profile, contour and a ‘model’ method based on energy-balance ablation modelling and measured net mass balance are discussed with respect to the reliability of the resulting series. The main findings of the investigations are: (1) The winter mass-balance series for Hintereisferner and Vernagtferner as determined with all applied methods result in a reliable climatologic average of 1000±100mmw.e. (2) When using the profile method, different spatial integration approaches are quite sensitive to the altitudinal coverage and the spatial pattern of observations. (3) The error of the model method occurs randomly, whereas contour as well as profile-method errors are more systematic. (4) Filtered time series from the two glaciers show similar tendencies for the last three decades.

Ice-thickness measurements are needed to calculate fluxes through fast-flowing outlet glaciers in Greenland, Alaska, Patagonia and Antarctica. However, relatively high attenuation of radio waves by dielectric absorption and volume scattering from englacial water hampers detection of the bed through warm deep ice. In the past we have had success measuring ice thickness of temperate glaciers using a ground-based monopulse radar system operating at low frequencies (2 MHz). Here we adapt the same system to operate from an airplane. Test flights over Bering Glacier, Alaska, USA, detected the bed through ice up to 1250m thick. Flights across the Seward–Malaspina Glacier system, Alaska, resolved the ice thickness of Malaspina Glacier, but strong hyperbolic-shaped returns obscured the bed echo through the Seward throat. It is likely that this clutter in the signal was caused by off-nadir returns from chaotic surface crevasses that are ubiquitous in the throat region.

Visible and infrared satellite images reveal numerous lineations on the Siple Coast region of West Antarctica. We used 5 MHz ice-penetrating radar to probe the interior and the bed of the ice sheet beneath a lineation at the boundary between Engelhardt Ice Ridge and flat-ice terrain to the south of the Kamb Ice Stream (KIS) outlet. Results show curved reflectors that emerge from the bed beneath 600 m thick ice. The tops of the reflectors extend about 100m into the ice above the bed, where they become almost horizontal. Apparent reflectivity of the horizontal section is about 20 dB less than that of the bed. We conclude that the likely cause of such strong reflection is sea water that was accreted into basal crevasses when the flat-ice terrain was floating. Internal layers are warped downward just downslope from the basal reflectors. It is thought that the downwarping was caused by localized basal melting in the past. The spatial pattern of downwarping suggests that localized basal melting was stronger on the north side than on the south side of KIS; apparently ice/ocean interactions on the two sides of KIS were different.

This paper describes a new compilation of both direct and geodetic mass-balance measurements, develops a procedure to reduce diverse balance measurements over different time-spans to common time-spans, and presents updated estimates of global average balance of small glaciers based on the enlarged compilation. Although geodetic measurements are fewer than direct measurements, they cover four times as many balance years. Direct and geodetic measurements are unbiased with respect to one another, but differences are often substantial. The statistical procedure can be understood by imagining that an n-year balance measurement is an average of a series of 1 year measurements. The series is hypothetical but we can calculate the uncertainty of each of its elements if, in addition to its measured average, we can also estimate its standard deviation. For this claim to be valid, the annual series must be stationary and normally distributed with independent (roughly, uncorrelated) elements, for which there is reasonable evidence. The need to know the standard deviation means that annual direct measurements from a nearby glacier, or equally reliable information about variability, are indispensable. Given this information, the new methodology results in moderately more negative balances. This is probably because tidewater glaciers are better represented in the geodetic data. In any case, the most recent published estimate of global average balance, 0.8–1.0mma–1 of sea-level equivalent for 2001–04, is now increased substantially to 1.1–1.4 mma–1 for 2001–05.

Results of a new glacier inventory of the upper Huasco valley, which lies within the arid Norte Chico zone of the Chilean Andes, are presented for 2004. Despite the high altitude, the glaciation in this region is limited in extent and is not classical mountain glaciation, which poses difficulties in completing standard inventory attribute tables. Small cornice-style ridgeline features constitute a large number of the non-transient ice bodies identified, and glaciers with surface areas <0.1 km2 comprise 18% of the glacierized area and 3% of the water resource stored as glacier ice within the Huasco valley. Rock glaciers are an important component of the cryosphere, comprising 12% of the total water volume stored in glacial features. Changes in glacier area over the last ~50 years are in line with those for glaciers in central Chile despite the contrasting climate conditions. Projections of glacier area change based on glacier hypsometry and zero isotherm shifts predicted using the PRECIS regional model temperature change for IPCC scenario B2 conditions suggest that the survival of 65% of glacier area and 77% of active rock-glacier area will be threatened under forecast conditions for the end of the 21st century.

Using the finite-element code Elmer, we show that the full Stokes modeling of the ice-sheet/ice-shelf transition we propose can give consistent predictions of grounding-line migration. Like other marine ice-sheet models our approach is highly sensitive to the chosen mesh resolution. However, with a grid size down to <5 km in the vicinity of the grounding line, predictions start to be robust because: (1) whatever the grid size (<5 km) the steady-state grounding-line position is sensibly the same (6 km standard deviation), and (2) with a grid-size refinement in the vicinity of the grounding line (200 m), the steady-state solution is independent of the applied perturbation in fluidity, provided this perturbation remains monotonic.