An ice-flow model is used to simulate the front variations of the Pasterze glacier, Austria. The model deals explicitly with the ice flux from sub-streams and tributaries to the main ice stream. The dynamic calibration method adopted successfully calibrates the ice-flow model under a non-steady-state condition. Despite the complexity of the glacier geometry, the ice dynamics of the Pasterze are adequately simulated.

Results of the sensitivity experiments suggest that the Pasterze glacier has been in a non-steady state most of the time and has a response time of 34–50 years.

Projections of the behaviour of the Pasterze in the next 100 years are made under various climate scenarios. Results suggest that the Pasterze will undergo a substantial retreat if there is future warming. The glacier is likely to retreat 2–5 km by the year 2100. The ice volume could be reduced by 24–63% by the end of the 21st century.

A numerical simulation of the flow of the cold glacier of Dôme du Goûter (4304 m, Mont Blanc, France) is presented. Owing to the large thickness of the firn layer, the simulation was done by using a rheological model for porous ice derived from a model for ceramic sintering and adapted to fit available data on in situ measured density profiles and firn mechanical behaviour. The flow calculation was made under the assumptions of axisymmetric geometry and stationary conditions, by solving a coupled problem. For a given density field, the velocities were obtained by the finite-element method. Then the integration of the mass-conservation equation along the streamlines derived from this velocity field gave the corresponding stationary densities. The results of the numerical simulation, besides the velocity and density fields, are the age of the ice along the streamlines. They are compared with observation and field data.

Finsterwalderbreen is a 35 km2, polythermal glacier in southern Spitsbergen, which last surged at around the start of the 20th century, and is still in the quiescent phase. Surface elevations measured several times since 1898 show strong thinning and retreat at the front, and gradual build-up in the accumulation area. Present-day annual velocities increase from about 1 m a–1 at the snout to a maximum of 13 m a–1 near the equilibrium line, then drop to 5 m a–1 at the bergschrund. Summer velocities are higher than the annual average, and winter velocities are lower (10–60% of summer velocities, where both are measured over several months). Velocities measured along the glacier centre line in 1950–52 are higher than those measured in 1994–95 in the upper basin, but are lower in the ablation area, which may be attributable to changes in surface profile on recovery from the surge. The measured modern ice flux (≈2.5 x 10–3 km3 a–1 w.e.) is only about 60% of the flux required to maintain balance (≈4.1 x 10–3 km3 a–1 w.e.), as determined from continuing mass-balance studies. Multi-frequency radar profiles show that most of the glacier is at the pressure-melting point at the base. This is consistent with hydrological studies showing high suspended-sediment loads, and hydro-chemical evidence of high sulphate levels, both indicating a well-developed basal drainage system during the summer which probably accounts for the seasonal velocity variations. Since the glacier still has a temperate bed, and appears to be building up in the accumulation area, it may be developing towards another surge.

A field of vectors showing the average velocity of the surging glacier Osbornebreen, Svalbard, was determined by comparing sequential SPOT (Système pour l’Observation de la Terre) and Landsat thematic mapper images. Crevasses which developed during the initial phase of the surge in the winter of 1986–87 were tracked using a fast Fourier chip cross-correlation technique. A digital elevation model (DEM) was developed using digital photogrammetry on aerial photographs from 1990. This new DEM was compared with a map drawn in 1966. The velocity field could be almost entirely determined with 1 month separation of the images, but only partly determined with images 1 year apart, due to changes of the crevasse pattern. The velocity field is similar to that found for Kronebreen, a continuously fast-moving tidewater glacier. No distinct zones of compressive flow were present and the data gave no evidence of a compression zone/surge front traveling downstream. The velocity field, the rapid advance of the terminus and the development of transverse crevasses in the upper accumulation area within a 6 month period may indicate that the surge developed as a zone of extension starting near the terminus and propagating quickly upstream. The crevasse pattern in the images is therefore interpreted to be the result of the extension zone traveling upstream, and, as the whole glacier starts to slide, the crevasse pattern alters according to the bedrock topography.

Radio-echo soundings in four different frequency bands ranging from 30 to 1000 MHz were compared with temperature measurements in boreholes in the accumulation area and ablation area of Finsterwalderbreen (77°26′ N, 15°15′ E), southern Spitsbergen. Finsterwalderbreen is a polythermal surge-type glacier in the quiescent phase after its last surge around AD 1900. The objective of the study was to investigate the relation between internal echos and the glacier ice temperature to map the overall thermal structure of the glacier. The thermal structure is important for ice flow velocities and hydrology of glaciers, and it also affects their ability to surge. At the borehole site in the accumulation area (three boreholes within a range of 60 m), a change in the relative amplitude of the reflected signal is detected in the 320–370 and 600–650 MHz bands at 52–55 m depth. The high-resolution temperature measurements with 2 m intervals show that the transition zone between cold and temperate ice corresponds to the change in the relative amplitude on the 320–370 and 600–650 MHz bandwidth data. The overall thermal structure of the glacier was mapped based on the radar sounding. The radar results show (a) that the glacier is at the pressure-melting point over most of its bed except within 500–700 m of the terminus, and (b) that there is an upper cold ice layer of variable thickness (25–170 m) underlain by temperate ice. This thermal structure is confirmed by the thermistor-instrumented access holes to the bed in both the accumulation and ablation zones of the glacier. The variations in the thermal structure in lower parts of the accumulation area are explained by superimposed ice and ice layers that cause variations in the downward heat transfer by refreezing of meltwater.

Frontal oscillations since the beginning of the 20th century are known at Glaciar Perito Moreno, an eastward outlet glacier of Hielo Patagónico Sur (southern Patagonia ice field). In 1900, the calving front was located about 1 km from the opposite bank. From 1935 to 1988, ruptures of ice-dams occurred at intervals of 1–5 years. Although this glacier has thus oscillated, it can be regarded as having been in a rather stable condition during the last half-century. Ice thickness in the ablation area has also remained unchanged from 1990 to 1996. The near-steady behavior of Glaciar Perito Moreno may be attributed to a regulating effect of the calving rate, namely, a decrease in the ablation amount due to calving with a retreat of the glacier.

Using 12 m long ablation poles, ice-flow velocities at the ablation area were measured several times in 1993 and 1994. The velocity in the early summer (November) was found to be slightly larger than the annual mean. It is concluded that basal sliding is significant throughout the year at this temperate glacier, with large fluctuations within a short period.

The continentality index is a good measure of the nature of the climate in a region, as it reflects not only the temperature but also the large-scale circulation. It correlates well with glacier mass-balance parameters. The climate along the west–east transect slightly north of the Arctic Circle across the Scandinavian Caledonides is governed by the prevailing westerlies; however, during winter the eastern part of the Caledonides is influenced by weather systems from the southeast. The differences in continentality meant by temperature and precipitation have a major impact on the response times of glaciers. The climatic change in this area has been dominated by increased summer mean temperatures (1910–20) and increased maritime influence since the 1980s. The slower-reacting glaciers on the Swedish side of the mountains are still adapting to the temperature increase at the beginning of the 20th century, and the increase in maritime influence gradually becomes less important towards the west. Thus, differences in the behaviour of glacier fronts along the west–east transect mirror differences in continentality.

As part of a study of ice dynamics, 25 boreholes were drilled with high-pressure hot water to the base of Haut Glacier d’Arolla, Switzerland. The boreholes were distributed across a half-section of the glacier, with closest spacing towards the glacier margin. The interior structure of the glacier was investigated by lowering a miniature color video camera down 11 boreholes, and the glacier bed was observed in 3 boreholes.

The video showed distinct changes in ice conditions with depth in the glacier, including ice foliations and changes in ice-bubble content and color. A basal ice layer of varying thickness was recorded in two of the three boreholes in which the glacier bed was observed. This basal ice layer was characterized by relatively sediment-rich, coarse-clear ice, and was thickest in a zone of high water-pressure fluctuations. The available evidence suggests that localized freezing at the glacier bed caused by variations in water pressure is the primary source of this basal ice layer. The near-surface transition between finegrained and coarse-bubbly ice, and between coarse-bubbly and coarse-clear ice, appeared to correlate with the depths of influence of diurnal and annual cold waves, respectively. Two types of ice foliation were identified, consisting of distinct planar coarse-clear ice, or more pervasive coarse-bubbly or fine-grained ice. An origin from crevasse closure or the transposition of sedimentary stratification provides the most likely explanation for these.

Glaciers of both the Arctic and mid-latitude mountain systems within Eurasia have retreated intensively during the past century. Measured and reconstructed glacier mass balances show that glacier retreat began around the 1880s. The mean annual mass-balance value for 1880–1990 was −480 mm a−1 for glaciers with maritime climatic conditions, and −140 mm a−1 for continental glaciers. It can be concluded that warming in the Caucasus occurred during at least the last 60 years, according to the distribution of crystal sizes in an ice core from the Dzhantugan firn plateau. Temperatures measured in 1962 at 20 m on the Gregoriev ice cap, Tien Shan, were −4.2°C while in 1990 they were −2°C, a warming of 2.2°C over 28 years. Changes in the isotopic composition of glacier ice during the 20th century indicate recent and continuing warming in different regions of Eurasia. The δ18O records reveal an enrichment at the Gregoriev ice cap during the last 50 years, while surface temperatures at the Tien Shan meteorological station have increased 0.5°C since 1930.

Combined measurements of meltwater discharge from the portal and of water level in a borehole drilled to the bed of Findelengletscher, Switzerland, were obtained during the later part of the 1993 ablation season. A severe storm, lasting from 22 through 24 September, produced at least 130 mm of precipitation over the glacier, largely as rain. The combined hydrological records indicate periods during which the basal drainage system became constricted and water storage in the glacier increased, as well as phases of channel growth. During the storm, water pressure generally increased as water backed up in the drainage network. Abrupt, temporary falls in borehole water level were accompanied by pulses in portal discharge. On 24 September, whilst borehole water level continued to rise, water started to escape under pressure with a resultant increase in discharge. As the drainage network expanded, a large amount of debris was flushed from a wide area of the bed. Progressive growth in channel capacity as discharge increased enabled stored water to drain and borehole water level to fall rapidly. Possible relationships between observed borehole water levels and water pressures in subglacial channels are influenced by hydraulic conditions at the base of the hole, distance between the hole and a channel, and the nature of the substrate.

Airborne scanning laser altimetry is a relatively new technique for remote sensing of ground elevation. A laser ranger is scanned across a swath beneath the aircraft, producing a two-dimensional distribution of elevations when combined with data on aircraft position and orientation. Smooth snow-covered glaciers are ideal surfaces for laser scanning since they are highly reflective. A new prototype laser system is described together with results from Hardangerjökulen, Norway. An analysis of the data shows that noise levels are very low at around 2 cm, and that repeatability between overlapping swaths is approximately ± 10 cm. This is consistent with an absolute accuracy of 15 cm or better from manufacturer’s and other measurements. Swath widths of over 1 km are attainable, at lowing complete coverage of small to medium-sized glaciers using parallel flight tracks. The high accuracy and dense, even coverage (about 20 000 points per km2) gives good-quality derived products such as DEMs and enables reliable measurement of glacier volume change. Scanning laser altimetry has many advantages over photogrammetry, which was previously the only remote-sensing method of measuring elevations over large areas but which performs poorly over snow-covered glaciers.

Velocity measurements carried out on Hintereisferner, Central Alps, Austria, provide the unique opportunity to study 100 years of ice dynamics of this glacier. During this time, three periods of accelerated flow occurred, around 1920, in 1940 and in the 1970s; but only around 1920 did the acceleration actually lead to an advance of about 60 m. The velocity increased from 30 m year−1 in 1914 to more than 120 m year−1 in 1919, and doubled during the accelerations of 1940 and 1980. In the course of the third event, the velocity increase spread over a period of more than a decade (1965–79) with a comparatively low maximum. These velocity changes cannot be explained by increased deformation velocity due to increased ice thickness alone.

Time series of the velocities at various locations along the glacier are given for the entire period, and an attempt was made to construct a time series of the velocity at a point 2 km from the strongly retreating front. The flow divergence was about 0.1 per year in the lowest 2 km, and emergence velocities reached 5 myear−1.

South Cascade Glacier, Washington, U.S.A., is one of the most extensively studied glaciers in the Western Hemisphere. In addition to mass-balance measurements, which date to 1958, numerous hydrological investigations have been carried out during the last three decades, and repeated ice-thickness determinations have been made using a variety of techniques. In the late 1960s, the basal topography was initially determined by gravitimetric methods. In the mid-1970s some of the fist depth measurements using radar on temperate ice were made. The basal topography was remapped soon after from a series of point radar measurements and boreholes drilled to the glacier bottom. During the 1990s, the ice thickness was remapped using digital recording of continuous profiles that obtained over 5000 ice-thickness measurements. Profiles have been corrected for the finite beamwidth of the antenna radiation pattern and reflections in steep terrain, resulting in a significantly improved depiction of the basal surface and internal structures. The map based on our recent radar profiles confirms the large-scale features of the basal topography previously depicted and reveals more structural detail. A bright reflector was detected at the base of the glacier and could be traced in adjacent profiles. Comparison with results from water-level measurements in boreholes drilled to the bed indicates that the reflector is a subglacial conduit.

The paper describes a non-linear dynamic thermomechanical model of the snow–firn development in the high mountain environment during the cold season. The model allows the estimation of the thermal and mechanical development of snow in an inhomogeneous stratified snowpack and in the upper part of the firn layer, the development of the snow–firn structure and transformation of snow into firn. The non-linear mathematical model consists of a heat problem, water-vapor diffusion problem, densification of the snow–firn layer, structural equations and strength conditions of the metamorphosed snow–firn depending on its temperature, structural parameters, density and rate of sublimation. The model includes conditions of snow transformation into firn and its structural changes. These equations are non-linear and coupled. The iterative finite-difference numerical method was used for the calculations and involved a full mathematical model. Numerical experiments were implemented using the complete model as well as real. The numerical computations were made for various meteorological data, temperature gradients and in the wide range of structural parameters. The calculation results were found to be in good agreement with the results of observations.

Recent developments of global positioning systems (GPS) have provided opportunities for rapid and frequent glacier mapping. The speed and accuracy of GPS techniques make them particularly suitable for repeated glacier mapping. Using two receivers in differential mode provides very high accuracy. In July 1995, a kinematic differential GPS survey of Austre Okstindbreen, one of the glaciers in the Norwegian national programme of mass-balance studies, provided three-dimensional positions of 2228 points in less than 6.5 h. Handling the data in a geographic information system permitted construction of a triangular irregular network digital terrain model (TIN DTM), which could be compared with TIN DTMs constructed from 1981 and 1993 surveys based on aerial photogrammetry and electronic distance measurement, respectively. A TIN DTM has advantages over the more usual grid-based DTM, which incorporates many interpolated values, although interpolation also is necessary if contours are to be derived from a TIN model. The aspect and gradient of triangular facets in the 1995 DTM have been used to produce a map which has considerable potential for detailed energy-balance studies of the glacier’s accumulation area.

Long-term observations of surface velocities and strain rates at the Norwegian glacier Austre Okstindbreen revealed both temporal and spatial variations. During a period of 6 years, the amount of ice passing through a cross-section slightly below the mean equilibrium-line altitude (1250 m) was some 30% less than the amount which accumulated above the equilibrium line. The mean horizontal component of surface velocity at the centre of the cross-section was of the order of 45–50 m a−1, whilst the thinner marginal ice moved less rapidly. At an altitude of about 1230–1200 m, surface velocities generally increased as the ice entered a steep icefall. In the lower part of the icefall, mean surface velocities again were of the order of 50 m a−1. From there, a general decrease down-glacier was evident, and longitudinal compression along the curving centre line of flow was accompanied by lateral extension. The contribution of internal deformation to surface flow at the lower part of the glacier, which was less than 150 m thick, is likely to have been relatively small, and between-year variations of the horizontal component of surface flow which affected a large area probably were a response to changes of basal sliding rates, reflecting variations of mass balance and water availability.

Results from analysis of kinematic global positioning system (GPS) surveying in Antarctica are examined to evaluate the feasibility of this technique for monitoring small glaciers. The experiment focuses on GPS field methods and a comparison of profiles surveyed in different years. Kinematic GPS has proved to be a useful method of monitoring glacier surfaces.

A series of boreholes were drilled with high-pressure hot water across a section of Haut Glacier d’Arolla, Switzerland, in summer 1995. Twenty-three of the boreholes were profiled with a digital inclinometer soon aller drilling, and 14 were re-profiled up to 6 weeks later to determine changes in the longitudinal shape of boreholes with time. In addition to the main surveys, three boreholes were surveyed 14 times each to assess the accuracy and reproducibility of inclinometry measurements. These repeat surveys suggest that caution is needed in the interpretation of short-term borehole displacement measurements, and that the reoccupation of boreholes from one year to the next may be a better way to determine patterns of internal deformation and basal sliding. The annual scale may also have advantages in providing more long-term insight into glaciological processes than short term (single season) measurements.

The paper describes the characteristics and behavior of unstable glaciers in the Pamirs, central Asia, over the last 30 years. These glaciers are divided into three groups: glaciers with registered surges; glaciers with signs of current instability; glaciers with signs of instability in the recent past. Examples are given to illustrate the nature of glacier fluctuations at a time of general glacier retreat in the Pamirs. A description is given of the evolution of Medvezhiy surging glacier, which was the subject of detailed annual observation from 1963 to 1991. Comparative characteristics of two building-up stages of this glacier, 1964–72 and 1974–88, are produced. On the basis of aerotopographic monitoring and with the aid of computer analysis, the kinematics of the glacier surface during the 1988–89 surge are traced in detail and the role of microsurges in the course of glacier build-up and the glacier surge is established.

As on all Peruvian Cordilleras, the glaciers on the eastern slopes of the Cordillera Blanca extend to generally lower elevations than those on the western slopes. The mountain range of Santa Cruz–Pucahirca possesses the largest east–west extension within the Cordillera Blanca. A significant retreat of the glaciers between two quasi-stationary situations around 1930 and 1950 was reconstructed from air photographs. The derived ELAs as well as the ∆ELA1930–50 show spatial diversities. The pattern of the ELAs is caused by differences in both accumulation and effective global radiation. The change in ELAs1930–50 is partly due to a spatially uniform increase in air temperature. The remaining rise of the equilibrium line, which varies within different parts of the investigation area, has to be related to changes in precipitation and effective global radiation. Both correspond to changes in air humidity which is suggested to be an important factor for tropical glacier fluctuations. A model of superposed typical tropical circulation patterns of different scales and ELA–climate model-based approaches is presented.