Introduction
During the past 25 years, many remote-sensing techniques have been applied to the study of glaciers, including measurements of (1) absolute ice thickness by radio-echo sounding from surface and airborne platforms, (2) changes in surface elevation over time with aerial photogrammetric methods and by geodetic airborne or spaceborne radar and laser altimetry, (3) delineation of the surface expression of glacier facies with satellite sensors, and (4) measurements of the fluctuations of the termini of valley glaciers and outlet glaciers (and margins) from ice fields, ice caps and ice sheets (Williams and Hall, in press). Often, satellite and airborne remote-sensing studies of glacier fluctuations do not have the benefit of ground observations to confirm the measurements. In this paper, however, we analyze Landsat images acquired between 1973 and 1992, in combination with field observations made during the same interval, to compare the measurement of changes in the termini of outlet glaciers and margin of Vatnajökull, Iceland (figs 1 and 2). Variations of Iceland’s glaciers are reported annually in Jökull.
The Earth’s smaller glaciers have generally been receding during the 20th century (Wood, 1988; Oerlemans, 1994). This has been accompanied by an estimated 1–3 mm a−1 rise in sea level (Paterson, 1994), perhaps as much as one-third of which has been caused by reduction in glacier-ice volume primarily from southeastern Alaska (Meier, 1984). The terminus-position changes of small glaciers often occur on a decadal time-scale, because smaller glaciers react more rapidly to climate change, according to Haeberli (1995) and Haeberli and Hoelzle (1995). Termini of Icelandic glaciers react to climate on a decadal time-scale regardless of their length (Sigurðsson and Jánsson, 1995). Sustained recession of a glacier terminus is usually indicative of a negative glacier mass balance, although it is possible for a glacier to increase in volume while its front is receding (Meier, 1984). In general, however, glaciers with a negative mass balance will retreat over time, although “large and relatively flat glaciers have a tendency of downwasting rather than retreating in case of rapid warming” (personal communication from W. Haeberli, 1996; cf. Haeberli and Hoelzle, 1995).
An ice cap is a dome-shaped mass of ice, < 50 000 km2 in area, that completely covers the terrain except for a few internal and marginal nunataks; it usually has one or more lobate or valley outlet glaciers. Ice caps are vulnerable to climate change, particularly those in the temperate areas of the Northern Hemisphere, because even small changes in the elevation of the equilibrium line can markedly affect the relationship of the size of the accumulation area to the ablation area. There is a marked decrease in the accumulation area in response to upward movement of the equilibrium line and a concomitant increase in the ablation area; as a result, total ice-cap area can shrink rapidly (e.g. Ok ice cap in Iceland (Williams, 1986a)), the shrinkage including recession of its outlet glaciers.
Background
Vatnajökull (Fig. 2), Iceland’s largest glacier (Björnsson, 1980a, 1980b; Williams, 1983), is the Earth’s most studied ice cap and represents a classical glaciological field site on the basis of the work in the late 18th century by the pioneer glaciologist S. Pálsson (Thórarinsson, 1960), the glaciological work in the late 19tlt century/early 20th century by Thoroddsen (1892) and pioneering glaciological research in the late 1930s and early 1940s by Ahlmann and Thórarinsson (1937–43). For the past century, Vatnajökull has been the focus of a vast array of glaciological studies by scientists from many nations, and much of the research has been published in Jökull and Geografiska Annaler.
For the past 45 years, Vatnajökull, the largest (8100 km2) of Iceland’s 13 ice caps > 10 km2 in area (Fig. 1), has been the subject of numerous remote-sensing investigations, using sensors mounted on ground, airborne and satellite platforms. A large variety of sensors have been used to (1) determine its thickness (and subglacier topography) (Björnsson, 1988), (2) measure its area (Williams, 1987), (3) delineate the surface expression of subglacier volcanoes and surface manifestations of subglacier volcanic and/or geothermal activity (Thórarinsson and others, 1974), (4) map glacier facies (Williams and others, 1991), (5) measure changes in surface elevation (Garvin and Williams, 1993), and (6) monitor and measure fluctuations of the margins of the outlet glaciers (Hall and others, 1992). For the 19 year period considered in this paper, most of its outlet glaciers measured were in recession, many as a result of post-surge downwasting.
A large ice cap such as Vatnajökull, with its 38 named and several unnamed outlet glaciers (Fig. 2), does not function as a unified mass of ice. Each of its outlet glaciers occupies its own hydrological basin (Björnsson, 1988); each responds to its unique subglacier topographic configuration and state of mass balance. In addition, many of Vatnajökull’s outlet glaciers surge periodically (Thórarinsson, 1969) or are affected by jökulhlaups caused by either subglacier volcanic/geothennal activity (Björnsson, 1988) or failure of ice-dammed lakes (Thórarinsson, 1974).
History of fluctuation measurements of Vatnajokull’s outlet glaciers
In the 1930s, systematic ground observations of the annual fluctuations of some of the termini of glaciers in Iceland, including some of the outlet glaciers of Vatnajökull, were initiated by Eypórsson (1931, 1949, 1963) and Eiriksson (1932). Since the publication of the first volume of Jökull in 1951, annual compilations of glacier fluctuations in Iceland (“Jöklabreytingar”) have been published by a succession of Icelandic glaciologists: J. Eypórsson (vols 1–16 (1951–66)), S. Rist (vols 17–37 (1967–87)), and O. Sigurðsson (vol. 38– (1988–present)). At the present time, the frontal positions of 19 named outlet glaciers from Vatnajökull are monitored on the ground (at 23 observation stations) generally on an annual basis, by trained lay persons who belong to the Ice-land Glaciological Society, other occasional observers and glaciologists. Ground-based measurements are made with a tape measure from a fixed point (e.g. cairn, steel pole) with an error of ± 1 m; in extreme cases, where it is difficult to determine the exact frontal position of an outlet glacier, the error can be as much as ± 5 m.
Landsat sensors
Landsat sensors have been operating since 23 July 1972, when the first multispectral scanner (MSS) and return-beam vidicon (RBV) instruments were launched on the ERTS 1 (later renamed the Landsat 1) satellite. The MSS sensor recorded images in four spectral bands, band 4 (0.5–0.6 μm); band 5 (0.6–-0.7 μm), band 6 (0.7–0.8 μm) and band 7 (0.8–1.1 μm), that covered parts of the visible and near-infrared parts of the electromagnetic spectrum at a picture clement (pixel) resolution of 79 m. An improved RBV sensor (0.505–0.750 μm panchromatic band) was flown on Landsat 3 (launched on 5 March 1978) and had a pixel resolution of 30 m. The thematic mapper (TM) sensor was first flown, along with the MSS sensor, on the Landsat 4 satellite (launched on 16 July 1978). The TM acquires images in seven spectral bands, band 1 (0.45–0.52 μm), band 2 (0.52–0.6 μm), band 3 (0.63–0.69 μm), band 4 (0.76–0.9 μm), band 5 (1.55–1.75 μm), band 6 (10.4–12.5 μm), and band 7 (2.08–2.35 μm), that cover parts of the visible, near-infrared, mid-infrared and thermal infrared spectrum. TM bands 1–5 and 7 have a pixel resolution of 30 m; TM band 6 has a pixel resolution of 120 m. The first three Landsat satellites (1–3) had orbital repeat cycles of 18 d; Landsats 4 and 5, however, image the same area of the Earth’s surface every 16 d, cloud cover permitting.
Glacier variation from satellite measurements
Landsat images were recognized soon after launch of the ERTS 1 spacecraft to be of value to investigations of glaciers (Williams and others, 1974; Krimmel and Meier, 1975; èstrem, 1975). In fact, Vatnajòkull was, on 8 September 1972, the second Icelandic glacier to be imaged by Landsat; the first was Mýrdalsjökull on 7 September 1972. During the past 24 years, Landsat images have been used to measure glacier variations using sequential images (Williams and others, 1974; Knight and others, 1987; Lucchitta and others, 1993); to carry out preliminary glacier inventories (Williams, 1986a); and to compile a satellite image atlas of glaciers of the world (Williams and Ferrigno, 1994). Landsat images have also been used to measure spectral reflectance of glacier surfaces and to determine if such reflectivity differences can be correlated with the surface expression of glacier fades (Hall and others, 1987; Williams and others, 1991), as defined by Benson (1962) and Benson and Motyka (1979) from Benson’s work on the Greenland ice sheet. Landsat images have also provided an important source of information for obtaining a global perspective of glaciers (Williams, 1986b; Williams and Hall, 1993).
Changes in glaciers at a number of locations throughout the world have been studied using Landsat data. Baseline information, using Landsat MSS images acquired in the mid-1970s, on glacier area and terminus position for all of the Earth’s glacierized regions is in preparation as an eleven volume U.S. Geological Survey Professional Paper 1386 A-K (Williams and Ferrigno, 1994); five of the eleven volumes have been published. Other work has shown that glacier recession in many areas, such as Alaska, Iceland and Austria, can be measured to an accuracy of approximately one Landsat pixel (Hall and others, 1992; Bayr and others, 1994). For example, the Pasterze glacier in the Austrian Alps receded >90 m between 1984 and 1990. Some of the outlet glaciers of Vatnajökull (Fig. 2) receded 1 km or more between 1973 and 1992 according to Landsat images. The tidewater Muir Glacier in Glacier Bay, Alaska, receded >7 km from 1973 to 1992, as measured from satellite images (Hall and others, 1995a).
Satellite images can be used to measure changes in glacier margins on a time-lapse basis, if the glacier margin can be spectrally discriminated from the terminal moraines and/or sandur (glacial outwash) deposits. In the current study, Landsat MSS, Landsat 3 RBV and TM images, spanning a 19 year period, were co-registered to measure changes in the margins of Vatnajökull. To carry out the comparative analysis between satellite-derived and ground-based measurements, only optimum Landsat images were used, with no cloud cover on one-half or more of the glacier and acquisition during the latter part of the melt season (to minimize the amount of marginal snowpack). Six Landsat images were used that fit these criteria (table 1): one MSS (22 September 1973), two Landsat 3 RBV (9 August 1978; and three TM (26 July 1987, 7 September 1991 and 19 October 1992). The cloud-free 22 September 1973 MSS image was used as the baseline image.
Methodology Of Measuring Change With Landsat Images
Registration of multiple images
In order to register Landsat images, one image is selected to be the baseline image. When images of different pixel resolution arc registered, the resolution of the registered images becomes that of the image with the least pixel resolution. Thus, when the MSS and TM images were registered for the study of Vatnajökull, the pixel resolution of the resulting data is 79 m.
To register images, first a set of “control points” is selected using image-processing software. Control points are fixed features that arc common to both images, such as a shoreline of a stable lake or a rock outcrop. Greater registration accuracy is generally achieved as the number of control points is increased. In order to register the Landsat images for this project, 175 control points were located and used to register each image. A review of control points for registration of the 1973 MSS and 1987 TM images found no discernible difference; a comparison of the 1973 MSS and 1992 TM images, however, found a one-pixel misalignment of registration on the northeastern margin of Vatnajökull that affected satellite-derived measurements of Brúarjokull and Eyjabakkajökull. Therefore, the 1973–92 measurement error for these two outlet glaciers is ± 191 m instead of 112 m (see “Measurement of changes between image dates” in next section). Also, depending on the type of printer used, the image geometry of computer-generated prints may be distorted.
Delineation of the glacier margin
Satellite images of a glacier permit the entire margin to be analyzed for changes, not just the 23 fixed survey stations for Vatnajökull from which marginal fluctuations are generally monitored and measured annually (see arrows in Figure 2). Once two different images are registered, one can delineate the changes in margins of ice caps and the termini of outlet glaciers. Delineating the exact position ofau ice-cap margin or terminus of an outlet glacier can be straight forward, if there is a sharp spectral difference between glacier ice and the surrounding terrain, On the computer monitor, a mouse-driven cursor is used to digitize a line that traces the margin of the ice cap; the glaciologist, preferably one with “local knowledge”, must decide where the precise margin is. Situations in which ambiguity exists as to where the margin is will be discussed in the next section. Various color lines can be used to outline the changing positions of the margins of ice caps and outlet glaciers. Glacier margins imaged in different years can then be compared (Fig. 3). One-can also “flicker” the registered images on the computer monitor to observe changes.
Analysis
Problems associated with delineating margins
Complicating factors in the delineation of the margin of Vatnajökull included the following: (1) lack of the third dimension in the Landsat images, (2) spectral similarities of supraglacial morainic debris and the terminal moraine and/or sandur (not usually a problem with an advancing outlet glacier, often a problem with a receding and down-wasting margin), (3) marginal snowpack, and (4) changes in spectral reflectivity caused by low solar elevation angle of illumination of the scene. To some degree, the delineation problem can be reduced by image enhancement, the presence of ice-margin (proglacial) lakes, location of the origin of glaciofluvial outwash channels and use of images that have solar-elevation angles of ≥25°.
Lack of a third dimension
Except in image-overlap areas of orbitally adjacent images, satellite images provide only a two-dimensional view of glaciers and the surrounding terrain, so that spectral differentiation becomes the key to accurately delineating the glacier margin. Stereoscopic images would help delineate the margin where there is a marked topographic difference between the glacier, terminal moraine and surrounding terrain.
Spectral similarity
When a glacier front has been stationary, receding or down-wasting, the entrained rock debris that is being carried by the glacier will appear on the glacier surface as the ice melts. This debris has the same spectral signature as the surrounding bedrock and surrounding glacial landforms, such as sandur and moraines. Under such conditions of spectral similarity, it is difficult to locate the position of the terminus of a receding glacier (Hall and others, 1992). Therefore, one must use other evidence to be able to delineate the terminus. In figures 3 and 4, ice-margin lakes are visible that define where the debris-covered glacier margin is located. Ice-margin lakes form at favorable locations along ice-cap and outlet glacier margins, wherever meltwater is impounded by the glacier on one side, sloping terrain exists in front of the glacier and a morainic dam is present. In some cases, tributary valleys are dammed by an outlet glacier forming a distributary, such as the small western lobe of Skeiðarárjökull, and may cause a lake to form (e.g. GrÖnalón; figs 2 and 4). Most such lakes change their shape seasonally, and some may be the source of periodic rapid-drainage events (jökulhlaups), such as Grænalón (Fig. 4). Advancing glaciers often have “clean” termini that are markedly dissimilar spectrally from the surrounding morainic and outwash deposits. Glacier termini that end in a lake or the ocean also have a large spectral dissimilarity between ice and water (Fig. 4). Under such circumstances the delineation of the glacier margin is simple.
Marginal snowpack
Snowpack can also cause problems in accurately delineating the margin of a glacier, by either masking terminal or lateral moraines or being difficult to differentiate from glacier ice. Under such conditions, it is easy to produce a false measurement of the glacier margin and glacier area. The best time of year to measure glacier change is at the end of the melt season, after the previous year’s snowpack has melted or reduced to a minimum. The time of minimum snow cover provides the opportunity to achieve good spectral contrast between the glacier margin and the surrounding terrain.
Low solar elevation angle of illumination of scene
The 19 October 1992 TM image, with a solar elevation angle of 14°. initially caused problems in accurately delineating the margin of outlet glaciers in southwestern Vatnajökull; the area of downwasting and recession of Tungnaárjökull between 1973 and 1992 was misidentibed, causing the margin to be delineated in the wrong location. A re-examination and reanalysis of the 1973–92 image pairs produced an accurate position of the margin.
Time-lapse interval
It is ideal if satellite data can be acquired at regular intervals, so that one can measure the changes in a glacier. If the interval between images is too short, there may not be measurable changes. On the other hand, if the time interval is too long, interim changes may be missed; therefore, the observed changes may be misinterpreted. Eyjabakkajökull (Fig. 3) surged about 0.5–lkm in 1972 and an additional 1 km in 1973, after being in a quiescent state since its last surge in 1938. Considerable recession had occurred between 1973 and 1992 (Fig. 3), so images taken prior to the 1972–73 surge that are compared with 1992 images would “miss” some of the 1.5–2 km advance.
Measurement of changes between image dates
Co-registration of sequential Landsat images and delineation of the terminus of an outlet glacier on the co-registered images permitted measurements to be made of the amount of advance or recession between two or more dates with an accuracy of ± 112 m for MSS images and for comparison with changes delineated on Landsat 3 RBV and TM data. The accuracy of Landsat 3 RBV and TM or TM and TM data is ± 42 m.
The ± 112 m (and ± 42 mj error is calculated as follows. A single measurement of glacier margin position will fall within or at the edge of a single 79 m pixel, on MSS 30 m pixel on Landsat 3 RBV and TM images) giving a range of ± 79 m for the width of a pixel to either side (±30 m for RBV and TM data). The second measurement of glacier margin will be the same error. The error of a series of measurements of similar quantities is where E is the error of measurement (79 m) and n is the number of measurements in the series (2). Therefore, the error is 79 m × = ± 112 m (or ± 42 m for RBV and TM data).
Using PCI EASI/PACETM software on a Silicon GraphicsTM workstation, the number of pixels can be counted between the two positions of the margin or terminus and multiplied by the size of the pixel (79 m in the case of superimposed images, one being an MSS image) to calculate the distance that the outlet margin or outlet glacier terminus has receded or advanced between the two dates. Table 2 summarizes these calculations for 1973–87, 1973–78. 1973–91, 1973–92, and 1987–92 registered image pairs. Most of the measurements of glacier fluctuation were done with the 1973–87 and 1973–92 image pairs.
Results
Changes in position of outlet glaciers
Changes of the margins of many of the outlet glaciers of Vatnajökull were measured using Landsat MSS and TM data from 1973–92. While most of the margins of the ire cap receded, some margins advanced and others were stable. The largest changes during the 19 year period occurred in the large lobate outlet glaciers along the southwestern, western and northern margins of Vatnajökull. Beginning at the southwestern margin of Vatnajökull at the Tungnaárjökull outlet glacier and proceeding in a counter-clockwise direction, we will summarize the changes measured using Landsat data and compare these measurements with ground-based observation? (table 2).
Comparison of satellite-derived and ground-based fluctuations
Table 2 summarizes satellite-derived fluctuation measurements on 11 of Vatnajòkull’s 38 named outlet glaciers (Fig. 2) and compares six of them with ground-based observations. Two outlet glaciers, Fláajökull and Kverkjökull, were measured on 1973 Landsat MSS and 1987 and 1991 Landsat TM images (table 1), but the changes were smaller than the pixel resolution of the MSS image. Two outlet glaciers from the Öræfajökull interior ice cap, Skaftafellsjökull and Svinafellsjökull,
were measured on the 1973 MSS and 1987 TM images, but the fluctuation was also too small to be resolved by the 79 m pixel resolution of the MSS image. Total ground-based measurements of fluctuation for the period 1973–87 were +110 ± 1 m for Skaftafellsjökull and +11 ± 1 m for Svinafellsjòknll, each too small a variation to be discerned on the MSS image.
Five outlet glaciers of Vatnajökull on its southwestern, southeastern, northwestern and western margins do not have systematic observations from established ground-based stations. In a counter-clockwise direction, they are: Skaftárjökull, Heinabergsjökull, Dyngjujökull, Köldukvslarjökull and Sylgjujökull Consequently, only the satellite-derived measurements of fluctuations for these five outlet glaciers are given in Table 2.
Three lobate outlet glaciers in southwestern Vatnajökull Tungnaárjòkull, Skaftárjökull and Siðujökull; Fig. 2) surged in 1994; they had all been in recession before 1973 and during the period 1973–92. Tungnaárjökull had the greatest recession of all of the lobate outlet glaciers in southwestern Vatnajökull. On 1973 MSS and 1987 TM images, it receded −1140 ± 112 m; ground-based observations during the same period showed a recession of −1130 ± 1 m. On 1973 MSS and 1992 TM images, it receded −1413 ± 112 m in the center near the ground-observation station, but only −474 + 112 m on the southern margin of the lobe. Ground measurements during the same period were −1380 ± 1m. Skaftárjökull had receded −485 ± 112 m on 1973 MSS and 1987 TM images and −425 ± 112 m on 1973 MSS and 1992 TM images; there is no ground-based observation station to corroborate these measurements.
On 1973 MSS and 1987 TM images, Siàujökull had receded −527 ± 112 m in its west-central part, −485 ± 112 m in its central part, and −513 ± 112 m in its east-central part. During the same period of time, ground observations showed a recession of −643 ± 1 m on its eastern part. On 1973 MSS and 1992 TM images, Siðujökull had receded −1030 ± 112 m; ground observations showed a recession of −1021 ± 1 m for the same period. Siðujökull last surged in 1963–64 at least 0.5 km, although field evidence suggests that the total surge was probably much more.
Skeiðarárjökull (Fig. 4), on the basis of 1973 MSS and 1987 TM images, had advanced +257 ± 112 m on its western lobe and had receded −150 ± 112 m in its central part. Comparable ground observations of the western lobe during the 1973–87 period showed an advance of +244 ± 1 m. Skeiðarárjökull began a surge in 1991, so the Landsat images that postdate 1991 will show an advance. In the eastern part of its margin, 1973 MSS and 1992 TM images show an advance of +474 ± 112 m; comparable ground observations showed an advance of +520 ± 1 m. 1987 TM and 1992 TM images show an advance of its terminus of +480 ± 42 m in the west-central part and +223 ± 42 m in the eastern part. Comparable ground observations for this latter period showed an advance of +276 ± 1 m. Changes in the western part of the lobe on 1987 and 1992 TM images were too small to measure on the Landsat images.
The post-surge (surge in 1972–73, for a total advance of 1.5–2 km determined from ground observations) recession of Eyjabakkajökull has been −790 ± 191 m on 1973 MSS and 1991 TM images, −921 ± 191 m on 1973 MSS and 1992 TM images (Fig. 3). There have been no ground measurements of the terminus of Eyjabakkajökull since 1985. On the 1973 MSS image its terminus was 2095 ± 79 m wide; on the 1992 TM image it was 1387 ± 30 m wide, reflecting the rapid downwasting and recession of the stagnant terminus following the 1972–73 surge.
Following a surge in 1963–64 of +8 km, Brúarjökull has been actively receding. Ground observations for the period 1977–88 show that its broad terminus had receded −1048 ± 5m. No ground observations were made during the 1972–77 and 1988–92 periods, although a linear extrapolation of downwasting and recession would give −2096 ± 5 m for the period 1972–92. Although Brúarjökull presents a difficult problem in delineating the glacier ice and terminal moraine contact, because of massive down-wasting all along its stagnant terminus, measurements were made on 1973 MSS and 1992 TM images. The terminus had, during this period of time, receded −1975 ± 191 m at the approximate position of the ground-observation station.
On 1973 MSS and 1991 TM images, Dyngjujökull’s margin had advanced +790 ± 112 m on the east-northeastern, +632 ± 112 m on the eastern, and +474 ± 112 m on the east-southeastern parts of its eastern margin; on the 1973 MSS and 1992 TM images its cast-southeastern margin advanced +553 ± 112 m (Fig. 5). On Köldukvislarjökull, on 1973 MSS and 1978 RBV images, it had advanced from +290 ± 112 in to +330 ± 112 m on its southern margin and +215 ± 112 m on its northern margin; on 1973 MSS and 1992 TM images its western margin had advanced between +285 ± 112 m on the north and +883 ± 112 m on the south-central part of the lobe. On 1973 MSS and 1992 TM images, Sylgjujökull receded −395 ± 112 m in the center of the lobe to between −553 ± 112m and −711 ± 112 m on its southwestern lobe.
Conclusions
In the specific case of Iceland, Landsat images permit the monitoring of fluctuations along most of Vatnajökull’s margin, not just at the 23 ground-observation stations (19 individual outlet glaciers) monitored on an annual basis al an accuracy of about 112 m with MSS images and 42 m for Landsat 3 RBV or TM images. Within the margins of error, we round good correspondence between the satellite-derived and ground-based measurements of six outlet glaciers of Vatnajökull during the period 1973–92. One aspect of our comparison of measurements of glacier fluctuations made from Landsat images vs ground-based observations is the critical importance of having “local knowledge” of glaciers to support the analysis of Landsat images.
Based on these results, we have confidence in the accuracy of the measurement of fluctuation of five additional outlet glaciers for which ground observations are not available. For four other outlet glaciers for which ground-observation stations exist, the ground-based measurements were ≤110 ± 1 m; no fluctuation could be discerned with satellite-derived measurements of these four outlet glaciers.
For the larger glaciers, such as ice fields, ice caps and ice sheets, the capability now exists to begin to build an archive of long-term data sets of the fluctuations of large glacier margins and area on a global basis, using glaciological data based on analysis of current satellite sensor data. If the World Glacier Monitoring service expands its area of interest from small glaciers to the larger glaciers (e.g. ice fields, ice caps, ice sheets and the outlet glaciers that comprise each), satellite images will be an important source of data for monitoring and measuring changes in glacier margins on a regional (and global) basis, in association with glaciologists who have “local knowledge” of glaciers of their specific region, and in the building of an international archive of long-term data sets. Improvements in the pixel resolution of imaging sensors, use of radar imagers, the annual acquisition of satellite images and acquisition of geodetically accurate laser altimetry from Earth-orbiting satellites will eventually provide considerably more data of value to glaciologists who are monitoring and measuring changes in glacier margins, area and volume (Williams and Hall, in press).