Skip to main content
×
×
Home

Information:

  • Access
  • Cited by 34
  • Cited by
    This article has been cited by the following publications. This list is generated based on data provided by CrossRef.

    Shukla, Tanuj Mehta, Manish Jaiswal, Manoj K. Srivastava, Pradeep Dobhal, D.P. Nainwal, H.C. and Singh, Atul K. 2018. Late Quaternary glaciation history of monsoon-dominated Dingad basin, central Himalaya, India. Quaternary Science Reviews, Vol. 181, Issue. , p. 43.

    Yue, Linwei Shen, Huanfeng Yu, Wei and Zhang, Liangpei 2018. Monitoring of Historical Glacier Recession in Yulong Mountain by the Integration of Multisource Remote Sensing Data. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, Vol. 11, Issue. 2, p. 388.

    Fountain, Andrew G. Glenn, Bryce and Basagic, Hassan J. 2017. The Geography of Glaciers and Perennial Snowfields in the American West. Arctic, Antarctic, and Alpine Research, Vol. 49, Issue. 3, p. 391.

    Swanger, Kate M. Lamp, Jennifer L. Winckler, Gisela Schaefer, Joerg M. and Marchant, David R. 2017. Glacier advance during Marine Isotope Stage 11 in the McMurdo Dry Valleys of Antarctica. Scientific Reports, Vol. 7, Issue. 1,

    DROLON, VANESSA MAISONGRANDE, PHILIPPE BERTHIER, ETIENNE SWINNEN, ELSE and HUSS, MATTHIAS 2016. Monitoring of seasonal glacier mass balance over the European Alps using low-resolution optical satellite images. Journal of Glaciology, Vol. 62, Issue. 235, p. 912.

    Funaki, Shota and Asaoka, Yoshihiro 2016. Long-term Change in Ablation Area of Tropical Glaciers by Landsat Data. Procedia Engineering, Vol. 154, Issue. , p. 168.

    Bisht, Pinkey Ali, S. Nawaz Shukla, Anil D. Negi, Sunil Sundriyal, Y.P. Yadava, M.G. and Juyal, Navin 2015. Chronology of late Quaternary glaciation and landform evolution in the upper Dhauliganga valley, (Trans Himalaya), Uttarakhand, India. Quaternary Science Reviews, Vol. 129, Issue. , p. 147.

    Ehlers, Todd A. Szameitat, Annika Enkelmann, Eva Yanites, Brian J. and Woodsworth, Glenn J. 2015. Identifying spatial variations in glacial catchment erosion with detrital thermochronology. Journal of Geophysical Research: Earth Surface, Vol. 120, Issue. 6, p. 1023.

    Carrivick, Jonathan L. Berry, Katie Geilhausen, Martin James, William H.m. Williams, Christopher Brown, Lee E. Rippin, David M. and Carver, Steve J. 2015. Decadal‐scale changes of the ödenwinkelkees, central austria, suggest increasing control of topography and evolution towards steady state. Geografiska Annaler: Series A, Physical Geography, Vol. 97, Issue. 3, p. 543.

    Ke, Linghong Ding, Xiaoli and Song, Chunqiao 2015. Heterogeneous changes of glaciers over the western Kunlun Mountains based on ICESat and Landsat-8 derived glacier inventory. Remote Sensing of Environment, Vol. 168, Issue. , p. 13.

    Nesje, A. 2014. Reference Module in Earth Systems and Environmental Sciences.

    Zhao, Liyun Ding, Ran and Moore, John C. 2014. Glacier volume and area change by 2050 in high mountain Asia. Global and Planetary Change, Vol. 122, Issue. , p. 197.

    Larsen, Christopher F. Motyka, Roman J. Arendt, Anthony A. Echelmeyer, Keith A. and Geissler, Paul E. 2014. 2009, State of the Cryosphere: Glaciers and Ice Sheets. p. 1.

    Sati, Sarswati Prakash Ali, Sheikh Nawaz Rana, Naresh Bhattacharya, Falguni Bhushan, Ravi Shukla, Anil Dutt Sundriyal, Yaspal and Juyal, Navin 2014. Timing and extent of Holocene glaciations in the monsoon dominated Dunagiri valley (Bangni glacier), Central Himalaya, India. Journal of Asian Earth Sciences, Vol. 91, Issue. , p. 125.

    Nawaz Ali, S. Biswas, R.H. Shukla, A.D. and Juyal, N. 2013. Chronology and climatic implications of Late Quaternary glaciations in the Goriganga valley, central Himalaya, India. Quaternary Science Reviews, Vol. 73, Issue. , p. 59.

    Nesje, A. 2013. Encyclopedia of Quaternary Science. p. 909.

    László, Péter Kern, Zoltán and Nagy, Balázs 2013. Late Pleistocene glaciers in the western Rodna Mountains, Romania. Quaternary International, Vol. 293, Issue. , p. 79.

    Chinn, T. Fitzharris, B.B. Willsman, A. and Salinger, M.J. 2012. Annual ice volume changes 1976–2008 for the New Zealand Southern Alps. Global and Planetary Change, Vol. 92-93, Issue. , p. 105.

    Arzhannikov, Sergei G. Braucher, Regis Jolivet, Marc Arzhannikova, Anastasia V. Vassallo, Riccardo Chauvet, Alain Bourlès, Didier and Chauvet, Frédéric 2012. History of late Pleistocene glaciations in the central Sayan-Tuva Upland (southern Siberia). Quaternary Science Reviews, Vol. 49, Issue. , p. 16.

    Bromley, Gordon R.M. Hall, Brenda L. Rademaker, Kurt M. Todd, Claire E. and Racovteanu, Adina E. 2011. Late Pleistocene snowline fluctuations at Nevado Coropuna (15°S), southern Peruvian Andes. Journal of Quaternary Science, Vol. 26, Issue. 3, p. 305.

    ×

Actions:

      • Send article to Kindle

        To send this article to your Kindle, first ensure no-reply@cambridge.org is added to your Approved Personal Document E-mail List under your Personal Document Settings on the Manage Your Content and Devices page of your Amazon account. Then enter the ‘name’ part of your Kindle email address below. Find out more about sending to your Kindle. Find out more about sending to your Kindle.

        Note you can select to send to either the @free.kindle.com or @kindle.com variations. ‘@free.kindle.com’ emails are free but can only be sent to your device when it is connected to wi-fi. ‘@kindle.com’ emails can be delivered even when you are not connected to wi-fi, but note that service fees apply.

        Find out more about the Kindle Personal Document Service.

        Map-based methods for estimating glacier equilibrium-line altitudes
        Available formats
        ×
        Send article to Dropbox

        To send this article to your Dropbox account, please select one or more formats and confirm that you agree to abide by our usage policies. If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your <service> account. Find out more about sending content to Dropbox.

        Map-based methods for estimating glacier equilibrium-line altitudes
        Available formats
        ×
        Send article to Google Drive

        To send this article to your Google Drive account, please select one or more formats and confirm that you agree to abide by our usage policies. If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your <service> account. Find out more about sending content to Google Drive.

        Map-based methods for estimating glacier equilibrium-line altitudes
        Available formats
        ×
Export citation

Abstract

We examine the validity of two methods for estimating glacier equilibrium-line altitudes (ELAs) from topographic maps. The ELA determined by contour inflection (the kinematic ELA) and the mean elevation of the glacier correlate extremely well with the ELA determined from mass-balance data (observed ELA). However, the range in glacier elevations above sea level is much larger than the variation in ELA, making this correlation unhelpful. The data were normalized and a reasonable correlation (r 2 = 0.59) was found between observed and kinematic ELA.The average of the normalized kinematic ELAs was consistently located down-glacier from the observed ELA, consistent with theory. The normalized mean elevation of the glacier exhibited no correlation and suggests that the toe–headwall altitude ratio is not a good approximation for the ELA. Kinematic waves had no effect on the position of the kinematic ELA. Therefore, topographic maps of glacier surfaces can be used to infer the position of the ELA and provide a method for estimating past ELAs from historic topographic maps.

Introduction

Field measurements of glacier mass change (e.g. Kasser, 1967,1973; Müller,1977; Haeberli,1985; Haeberli and Müller, 1988; Østrem and others, 1991; Haeberli and Hoelzle, 1993; Haeberli and others, 1998; Krimmel, 1999) are used to understand the response of glaciers to climatic variations (e.g. Holmlund,1996; Dyurgerov and Meier, 2000; McCabe and others, 2000). It is also important to understand the contribution of glacier melt to regional stream-flow variations (e.g. Meier, 1969; Braithwaite and Olesen, 1988) and to global sea-level change (Meier, 1984; Oerlemans and Fortuin, 1992; Dyurgerov and Meier, 1997). However, mass-balance studies are time-intensive and can be applied to only a small fraction of the glaciers in any given region. Unfortunately, methods to rapidly assess regional glacier change have been elusive. Monitoring changes in terminus position is an attractive method because glacier termini are clearly visible and easily monitored by aerial and ground-based photographic surveys (Gilbert, 1904; LaChapelle, 1962; Veatch, 1969; Aniya, 1988) and by satellite imagery (Williams and Ferrigno, 1993; Duncan and others, 1998; Paul, 2002). Inferring climate change or glacier mass change can be difficult because one must separate the dynamic response of the glacier from the changes in glacier mass (Nye, 1962). Recent advances in understanding the time-scale of glacier responses to climate change (Jóhannesson and others, 1989) have increased our capability to interpret the lag time between climate variation and terminus response. However, to track climate change through changes in terminus position is still problematic because of kinematic waves (Nye, 1965; Van de Wal and Oerlemans, 1995). To apply time-scale response methods also requires estimates of glacier thickness and some knowledge of mass change at the terminus, which are unknown for most glaciers.

One appealing alternative to measuring terminus position is to find the equilibrium-line altitude (ELA). The position of the ELA is controlled entirely by climatic processes (Kuhn, 1981; Ohmura and others, 1992; Seltzer, 1994; Fountain and others, 1999). Up-glacier from this “line”, the accumulation zone gains more mass annually than it loses. Down-glacier, in the ablation zone, the glacier loses more mass than it gains. At the ELA, the annual mass change is zero (Paterson,1994). The most accurate method to determine the ELA is by field measurements (e.g. Østrem and Brugman, 1991). One approach is to contour the point measurements of mass balance on a map to determine the glacier’s ELA for that year. Such field methods are slow and expensive and therefore are only applied to a few glaciers. These constraints have led to a variety of remote monitoring techniques to find the position of a glacier’s ELA. Assuming that the snowline at the end of the summer is a good approximation of the ELA, then aerial and ground-based photography are attractive methods to monitor inferred changes in ELA (LaChapelle, 1962). This method is difficult in practice because the timing of the photography is critical. Observations have to be collected just prior to the start of the accumulation season, and early snowfall or late warming events may require repeated imaging. Also, poor weather, common to mountainous regions, can delay image acquisition for weeks.This method cannot easily be applied to environments where accumulation and ablation seasons do not neatly fit into the climatic seasons of winter and summer, respectively. In the dry valleys of Antarctica, for example, snowfall occurs year-round on all surfaces and distinct accumulation/ablation seasons do not exist (Fountain and others,1998). In the Himalaya, the accumulation and ablation seasons occur simultaneously during the monsoonal season (Ageta and Higuchi, 1984), and in the equatorial regions multiple seasons of accumulation and ablation exist (Kaser,1999).

An alternative method for determining ELA is based on the inflection of elevation contour lines on a topographic map of a glacier (Østrem, 1966; Porter, 1975). This method was first introduced by Hess (1904) to determine firn-line locations, and in the European scientific literature it is sometimes referred to as the “Hess method”. The accumulation zone of the glacier is bowl-shaped, because snow accumulating on the margins of the glacier advects toward the center. Below the ELA, in the ablation zone, the contour lines are convex, because mass is lost from all sides and ice is advected away from the glacier center towards the margins. The transition or inflection from a concave to convex contour is a relatively flat surface that should be close to the location of the long-term averageELA (Fig.1). If this is true, then the topography of the glacier surface can be used to infer the average location of the ELA. Because the ELA position determined by the contour method results in part from glacier motion, we term this ELA the “kinematic ELA” to distinguish it from the year-to-year ELA determined by observations of glacier surface mass changes, or “observed ELA”.

Fig. 1. Contour map of South Cascade Glacier, Washington, U.S.A., from the United States Geologic Survey Dome Peak 7.5′ quadrangle. According to the map, the glacier was field-checked in1965.

This paper examines the correspondence between the kinematic and observed ELA. To our knowledge, the “kinematic” (“Hess”) method has not been empirically tested in this way. If these two values correspond, historical maps of glaciers can be used to estimate ELA change since the turn of the 20th century, when reliable topographic maps of alpine regions became available, and modern maps can be used to estimate ELAs for glaciers on which mass-balance studies are not performed.We also evaluate the correlation between ELA and mean glacier elevation. The normalized form of mean glacier elevation, which we use, is similar to the toe–headwall altitude ratio (THAR), a procedure for estimating ELA in the glacial geologic literature (Meierding, 1982; Hawkins, 1985; Torsnes and others,1993).

Table 1. The glaciers used in this study

Methods

We collected topographic maps of glaciers for which observed ELAs have been recorded by the World Glacier Monitoring Service (WGMS). The 40 glaciers we examined are listed in Table 1. Our dataset is biased towards glaciers in the Northern Hemisphere, but that should not change our results. We selected the contour that best represented the inflectionbetween the regions of surface concavity (accumulation zone) and convexity (ablation zone). The altitude represented by this contour is defined as the kinematic ELA. In some cases it is difficult to determine the specific contour defining the inflection of the glacier surface because the contour line itself may contain concave and convex segments. In these situations the line that shared the most symmetrical distribution of concavity and convexity was selected. The error of the selection of the contour line is, at minimum, the error of the contour map, typically half of the contour interval. In those cases where several adjacent contour lines seemed suitable, we chose the average. The date of the glacier surface depicted was determined from either the date of photography used to make the map or the date the surface was surveyed. If a “field check” was performed on the glacier, the map’s date was set to that of the fieldwork. Information on the date of aerial photography or field checking is typically included in the map legend.

The highest and lowest elevations of each glacier were recorded from the maps and checked against values reported to theWGMS. We used them to calculate the mean elevation of each glacier (maximum elevation plus minimum elevation divided by two). This is not the area-weighted mean discussed by Braithwaite and Müller (1980), rather it corresponds to their “E CR”. Kurowski (1891) introduced this method of estimating a glacier’s firn line, and Porter (1975) suggested the mean elevation of a glacier is a good estimate of the ELA, but neither author presents data to support this assertion.

ELA data were obtained from the WGMS reports (Kässer, 1967, 1973; Müller, 1977; Haeberli, 1985; Haeberli and Müller, 1988; Haeberli and Hoelzle, 1993; Haeberli and others, 1998). We averaged the observed ELAs reported for the years 1965–95, and compared this average with the kinematic ELA for a date within this time period. We did not calculate steady-state ELAs because most glaciers are not in a steady state.We chose to average the whole record, rather than select a smaller time interval, both because the kinematic ELA results from long-term climate trends and because we did not want abnormally high or low ELAs to bias the average.Typically, the maps were prepared at the beginning of the mass-balance program and precluded averaging the observed ELA record for the years preceding the map date, as might be suggested by response-time theories (Jóhannesson and others,1989).

Results

Table 2 summarizes the results of our kinematic ELA estimates and the observed ELAs and mean elevation.The correlation between the kinematic and observed ELAs is very good (Fig. 2a), with a correlation coefficient, r 2, of 0.99 (Table 3). The slope of the least-squares fit is 1.05, only slightly steeper than a one-to-one ratio. The intercept of the equation is +10.8 m, indicating the observed ELA is a bit higher than the kinematic ELA. The relation between mean glacier altitude and observed ELA (Fig. 2b) is similar to Figure 2a, with a slope slightly smaller (0.99) and a larger intercept, 46 m, with the observed ELA higher than the mean altitude. Correlations were also calculated between the observed ELA for only that year in which the map was made, and the kinematic ELA and mean altitude (Fig. 2c and d). Results were similar to those of the previous findings, but with lower correlation coefficients and larger intercepts, as one might expect from the year-to-year variability of the observed ELA relative to the mean ELA.

Fig. 2. Relations between the observed ELA and the kinematic ELA and mean elevation: (a) kinematic ELA vs mean observed ELA; (b) mean glacier elevation vs mean observed ELA; (c) kinematic ELA vs same-year observed ELA; (d) mean glacier elevation vs same-year observed ELA.

Table 2. The data used in this study

Table 3. Regression equations for the studied ELA relationships corresponding to Figures 2 and 3

The correlation between kinematic and observed ELAs should be good, since the two values are constrained to the altitude range of each glacier, while the range itself varies greatly among the glaciers. Our sample is also dominated by relatively small alpine glaciers, on which many mass-balance studies are based. To remove the spurious correlation induced by the wide range of altitudes, we normalized the data by the altitude range of the glacier,

(1)

where E is the normalized value, subscript i indicates either the kinematic (k) or observed values (o), ELA is observed or kinematic ELA, z 0 is the minimum glacier altitude, and, z m is the maximum glacier altitude. The normalization limits the altitude of each ELA to a range of 0–1. The mean elevation for all glaciers becomes 0.5. Averages of the normalized data showed that the kinematic ELA was about 20% lower than the observed ELA. An F test on the normalized data showed that the kinematic and observed ELAs were not statistically different. The mean elevation data cannot be analyzed in the same way against the observed ELAs, as this would compare a single point (0.5) with a variable dataset. As can be seen in Figure 3c and d, there is no correlation between the normalized mean elevation and normalized observed ELA.

Fig. 3. The relation between the normalized values of the kinematic and mean observed ELAs and the mean elevation and mean observed ELAs: (a) normalized kinematic ELA vs mean observed ELA with all glaciers included; (b) normalized kinematic vs mean observed ELA for glaciers with >15 years of data; (c) normalized mean elevation vs mean observed ELA for all glaciers; (d) normalized mean elevation vs mean observed ELA for glaciers with >15 years of observed ELA data.

The relation between the normalized observed and kinematic ELA was poor (Fig. 3a; r 2 = 0.1). Concerned that short time series of ELAs could bias the relation, we replotted our data after removing all glaciers with observation records of <15 years. The results for these long-record glaciers are much better (Fig. 3b; Table 3: r 2 = 0.59).

Glaciers with kinematic waves produce large variations in terminus position that might strongly affect the application of the proposed method. We compared the kinematic ELA with the mean glacier elevation for a series of 12 maps of Nisqually Glacier, Mount Rainier, Washington, dating from 1913 to 1994. Nisqually Glacier is well known for its kinematic waves (Meier, 1962), but its mass balance has never been measured due to the difficulty of the glacier terrain. Our results show that the mean elevation is far more variable than the kinematic ELA (Fig. 4). The change in kinematic ELA is coincident with extensive glacier retreat from 1913 to 1940 (T. Nylen, 2003). Since 1940 the ELA has been relatively stable and, relative to the errors in the method, the ELA can be considered constant. We believe the pre-1940 rise in ELA is real and is perhaps an artifact of the Little Ice Age in this region. After 1940, the variation in the mean elevation is large compared to the kinematic ELA and is a result of kinematic waves which temporarily lengthen the glacier. In addition, defining the terminus of such glaciers is problematic because of the large portions of “dead” ice that accompany large advances from the kinematic waves (T. Nylen, unpublished information).

Fig. 4. Estimates of the ELA with time using the mean elevation and topographic methods for Nisqually Glacier which regularly experiences kinematic waves. No kinematic ELA was recorded for the years 1925–45 because the maps from those years only included the ablation area.

Discussion

The correlation between the kinematic and observed ELA is good, but we were initially puzzled by the offset (20% of the total glacier elevation in Fig. 3b) between the normalized kinematic and observed ELAs. Hooke (1998) points out that an offset should be expected because glaciers move down slopes. For a steady-state glacier, the net balance (b n) is directly related to the emergence (vertical) velocity (w s) at the surface. For a point on a surface with slope α and horizontal velocity u s,

(2)

At the observed ELA, b n = 0, and since u s is positive and the slope is negative, the emergence velocity w s is negative, or into the glacier. Thus, at the observed ELA, mass is being advected into the glacier and the kinematic ELA of zero emergence velocity is located down-glacier. Alternatively, where w s is equal to zero (kinematic ELA), the balance is negative, placing the kinematic ELA lower on the glacier than the observed ELA, in the ablation zone, as shown by our results. Another possible cause of this offset is that the map was made before field observations were collected and the ELA has risen in elevation since that time.

The mean elevation of a glacier correlates well with the observed ELA (Fig. 2b and d) when elevations above sea level are used. This is spurious because the relation is biased by the large range in glacier elevations. Moreover, as previously described, the normalized ELA shows great scatter about the normalized mean elevation, fixed at 0.5 (Fig. 3c and d). Normalized mean elevation is similar to the THAR method (Meierding, 1982; Hawkins, 1985; Torsnes and others, 1993). The THAR method is used principally in glacial geology to determine a former glacier’s ELA based on valley morphology. The difference in elevation between the former glacier’s terminus, or toe, (usually determined from moraine location) and the elevation of the former glacier’s headwall is multiplied by the THAR to estimate the elevation of the former equilibrium line relative to the terminus. ATHAR of 0.4 means that the ELA is 40% of the glacier’s elevation range higher than the terminus. Both Meierding (1982) and Torsnes and others (1993) compared the THAR of 0.4 with the accumulation–area ratio (AAR) of 0.6, an equilibrium value for “normal” alpine glaciers (Meier and Post, 1962), and found a good correlation.

For the small mountain glaciers used in this study we do not observe a constant THAR, as the normalized ELA is quite variable from glacier to glacier. This brings into question the validity of the THAR method for small alpine glaciers. However, it may be useful in the broadest of contexts because the average of the normalized observed ELAs is 0.55, close to the normalized value of glacier elevation (THAR) of 0.5.

A strong word of caution: our methods to estimate the ELA from maps or average glacier altitude and THAR apply only to alpine glaciers that terminate “normally” on dry land. Like the application of the AAR, these two methods do not apply to glaciers that terminate in marine or lacustrine waters or to glaciers that terminate on cliffs. In these situations, the relation between mass balance and altitude is strongly biased by calving at the glacier terminus. The calving flux is largely controlled by processes at or near the terminus and, for marine and lacustrine systems, is largely independent of annual climate (Brown and others, 1982; see papers inVan der Veen, 1997).

For “normal” glaciers the kinematic ELA will provide a reasonable approximation of a glacier’s average ELA at far less expense than a traditional mass-balance study. Recent advances in remote-sensing techniques for mapping glacier surfaces (Echelmeyer and others, 1996; Allen, 1998; Favey and others,1999; Kääb and Funk,1999; Hubbard and others, 2000) will make the kinematic ELA method even more practical.

Conclusions

The “kinematic” method of estimating the ELA, by selecting the topographic contour representing the inflection between the concave and convex surfaces on a glacier, is a reasonable approximation for the actual, long-term, ELA acquired through traditional mass-balance studies. The kinematic ELA occurs lower on the glacier than the observed ELA, consistent with theoretical models of glacier flow and with a climate-driven rise in the position of the ELA during the period of observed ELA data collection. Distortions to the glacier surface by kinematic waves do not seem to affect the results of this method. Application of this method is limited to those “normal” alpine glaciers that do not terminate in cliffs or in water.

The mean elevation of a glacier does not provide a reasonable approximation for the observed ELA but may be useful as the most general of indices. The THAR method, an approach common to glacial geological studies, may be applicable in a broad sense but should not be used with individual glaciers. Our work suggests an appropriateTHAR for modern alpine glaciers (using the highest elevation on each glacier as the headwall altitude) is 0.55, much higher than the value of 0.4 commonly used.

These methods allow the examination of the spatial and temporal variation in ELA starting in the early 20th century when reliable topographic maps of alpine regions became available. With increasing access to airborne and satellite imagery, practical application for constructing and revising topographic maps will increase in time. By applying the kinematic method to these topographic maps, we can remotely infer ELAs for many glaciers over large regions and track changes in ELA through time. This will increase our understanding of variations in glacier activity in relation to climatic variations and local topographic settings.

Acknowledgements

We are grateful to G. Østrem and M. Dyurgerov for kind assistance in locating glacier maps for our study. G. Østrem was also very helpful in tracking down some of the original Hess material. In that effort we also owe a debt to J. O. Hagen and L. Braun. T. Nylen provided numerous maps of Nisqually Glacier from his M.S. thesis and we appreciate his willingness to share his results with us. J. Cook was helpful in the early stages of this project. Part of this work was supported by U.S. National Science Foundation grant OPP-9810219 which is gratefully acknowledged. Reviews by C. Waythomas, G. Østrem and an anonymous reviewer helped strengthen the manuscript.

References

Ageta, Y. and Higuchi, K.. 1984. Estimation of mass balance components of a summer-accumulation type glacier in the Nepal Himalaya. Geogr. Ann., 66A(3), 249255.
Allen, T. R. 1998. Topographic context of glaciers and perennial snowfields, Glacier National Park, Montana. Geomorphology, 21(3–4), 207216.
Aniya, M. 1988. Glacier inventory for the Northern Patagonia Icefield, Chile, and variations 1944/45 to 1985/86. Arct. Alp. Res., 20(2), 179–2.
Braithwaite, R. J. and Müller, F.. 1980. On the parameterization of glacier equilibrium line altitude. International Association of Hydrological Sciences Publication 126 (Riederalp Workshop 1978 — World Glacier Inventory), 263–271.
Braithwaite, R. J. and Olesen, O. B.. 1988. Effect of glaciers on annual runoff, Johan Dahl Land, south Greenland. J. Glaciol., 34(117), 200207.
Brown, C. S., Meier, M. F. and Post, A.. 1982. Calving speed of Alaska tidewater glaciers, with application to Columbia Glacier. U.S.Geol. Surv. Prof. Pap. 1258-C.
Duncan, C. C., Klein, A. J., Masek, J.G. and Isacks, B. L.. 1998. Comparison of late Pleistocene andmodern glacier extents in central Nepal based on digital elevation data and satellite imagery. Quat. Res., 49(3), 241254.
Dyurgerov, M. B. and Meier, M. F.. 1997. Year-to-year fluctuations of global mass balance of small glaciers and their contribution to sea-level changes. Arct. Alp. Res., 29(4), 392402.
Dyurgerov, M. B. and Meier, M. F.. 2000. Twentieth century climate change: evidence from small glaciers. Proc.Natl. Acad. Sci.U.S.A., 97(4), 14061411.
Echelmeyer, K. A. and 8 others.1996. Airborne surface profiling of glaciers: a case-study in Alaska. J. Glaciol., 42(142), 538547.
Favey, E., Geiger, A., Gudmundsson, G.H. and Wehr, A.. 1999. Evaluating the potential of an airborne laser-scanning system for measuring volume changes of glaciers. Geogr. Ann., 81A(4), 555561.
Fountain, A. G., Dana, G. L., Lewis, K. J., Vaughn, B. H. and McKnight, D. M.. 1998. Glaciers of the McMurdo Dry Valleys, southern Victoria Land, Antarctica. In Priscu, J. C., ed. Ecosystem dynamics in a polar desert: the McMurdo DryValleys, Antarctica. Washington, DC, American Geophysical Union, 6575. (Antarctic Research Series 72.)
Fountain, A. G., Lewis, K. J. and Doran, P.T.. 1999. Spatial climatic variation and its control on glacier equilibrium line altitude inTaylorValley, Antarctica. Global Planet. Change, 22(1–4), 110.
Gilbert, G. K. 1904. Variations of Sierra glaciers. Sierra Club Bulletin, 5(1), 2025.
Haeberli, W., comp. 1985. Fluctuations of glaciers 1975–1980 (Vol. IV). Paris, International Commission on Snow and Ice of the International Association of Hydrological Sciences/UNESCO.
Haeberli, W. and Hoelzle, M., comps.1993. Fluctuations of glaciers1985–1990(Vol. VI). Wallingford, Oxon, IAHS Press; Nairobi, UNEP; Paris, UNESCO.
Haeberli, W. and Müller, P., comps. 1988. Fluctuations of glaciers 1980–1985(Vol. V). Wallingford, Oxon, IAHS Press; Nairobi, UNEP; Paris, UNESCO.
Haeberli, W., Hoelzle, M., Suter, S. and Frauenfelder, R., comps. 1998. Fluctuations of glaciers 1990–1995with addendas from earlier years (Vol.VII). Wallingford, Oxon, IAHS Press; Nairobi, UNEP; Paris, UNESCO.
Hawkins, F. F. 1985. Equilibrium-line altitudes and paleoenvironment in the Merchants Bay area, Baffin Island, N.W.T., Canada. J. Glaciol., 31(109), 205213.
Hess, H. 1904. Die Gletscher. Braunschweig, Vieweg & Sohn.
Holmlund, P. 1996. Mass balance studies in northern Sweden. Z. Gletscherkd. Glazialgeol., 31, Part 1, 1995, 105114.
Hooke, R. LeB. 1998. Principles of glacier mechanics. Upper Saddle River, NJ, Prentice Hall.
Hubbard, A. and 6 others. 2000. Glacier mass-balance determination by remote sensing and high-resolution modelling. J. Glaciol., 46(154), 491498.
Jóhannesson, T., Raymond, C. and Waddington, E.. 1989. Time-scale for adjustment of glaciers to changes in mass balance. J. Glaciol., 35(121), 355–2.
Kääb, A. and Funk, M.. 1999. Modelling mass balance using photogrammetric and geophysical data: a pilot study at Griesgletscher, Swiss Alps. J. Glaciol., 45(151), 575583.
Kaser, G. 1999. A review of the modern fluctuations of tropical glaciers. Global Planet. Change, 22(1–4), 93103.
Kasser, P. 1967. Fluctuations of glaciers 1959–1965 [Vol. I]. Paris, International Commission of Snow and Ice of the International Association of Scientific Hydrology/UNESCO.
Kasser, P., comp. 1973. Fluctuations of glaciers 1965–1970 [Vol. II]. Paris, International Commission on Snow and Ice of the International Association of Hydrological Sciences/UNESCO.
Krimmel, R. M. 1999. Water, ice, meteorological and speed measurements at South Cascade Glacier, Washington, 1998 balance year. U.S. Geol. Surv. Water-Resour. Invest. Rep. 99-4049.
Kuhn, M. 1981. Climate and glaciers. International Association of Hydrological Sciences Publication 131 (Symposium at Canberra 1979 — Sea Level, Ice and Climatic Change), 320.
Kurowski, L. 1891. Die Hohe der Schneegrenze mit besonderer Berucksichtigung der Finsteraarhorn-Gruppe. Geographische Abhandlungen, Berlin Universitat, 5(1), 119160.
LaChapelle, E. 1962. Assessing glacier mass budgets by reconnaissance aerial photography. J. Glaciol., 4(33), 290297.
McCabe, G. J., Fountain, A. G. and Dyurgerov, M.. 2000.Variability in winter mass balance of Northern hemisphere glaciers and relations with atmospheric circulation. Arct. Antarct. Alp. Res., 32(1), 6472.
Meier, M. F. 1962.The kinematic wave on Nisqually Glacier, Washington. J. Geophys. Res., 67(2), 886.
Meier, M. F. 1969. Glaciers and water supply. American Water Works Association, Journal, 61(1), 812.
Meier, M. F. 1984. Contribution of small glaciers to global sea level. Science, 226(4681), 14181421.
Meier, M. F. and Post, A. S.. 1962. Recent variations in mass net budgets of glaciers in western North America. International Association of Scientific Hydrology Publication 58 (Symposium at Obergurgl 1962 — Variations of the Regime of Existing Glaciers), 6377.
Meierding, T. C. 1982. Late Pleistocene glacial equilibrium-line altitudes in the Colorado Front Range: a comparison of methods. Quat. Res., 18(3), 289310.
Müller, F., comp. 1977. Fluctuations of glaciers 1970–1975 (Vol. III). Paris, International Commission on Snow and Ice of the International Association of Hydrological Sciences/UNESCO.
Nye, J. F. 1962. Implications of mass balance studies. J. Glaciol., 4(33), 264–2.
Nye, J. F. 1965. A numerical method of inferring the budget history of a glacier from its advance and retreat. J. Glaciol., 5(41), 589607.
Nylen, T. 2003. Spatial and temporal variations of glaciers on Mt. Rainier between 1913 and 1994. (M.S. thesis, Portland State University.)
Oerlemans, J. and Fortuin, J. P. F.. 1992. Sensitivity of glaciers and small ice caps to greenhouse warming. Science, 258(5079), 115117.
Ohmura, A., Kasser, P. and Funk, M.. 1992. Climate at the equilibrium line of glaciers. J. Glaciol., 38(130), 397411.
Østrem, G. 1966. The height of the glaciation limit in southern British Columbia and Alberta. Geogr. Ann., 48A(3), 126138.
Østrem, G. and Brugman, M.. 1991. Glacier mass-balance measurements. A manual for field and office work. Saskatoon, Sask., Environment Canada. National Hydrology Research Institute. (NHRI Science Report 4.)
Østrem, G., Haakensen, N., Kjöllmoen, B., Laumann, T. and Wold, B.. 1991. Massebalansemalinger pa norske breer 1988 og 1989: resultater fra feltarbeid og beregninger. Norges Vassdrags- og Elektrisitetsvesen. Vassdragsdirektoratet, Publikasjon 11.
Paterson, W. S. B. 1994.The physics of glaciers. Third edition. Oxford, etc., Elsevier.
Paul, F. 1997. Changes of glacierarea in the Austrian Alps between 1973 and1992 derived from Landsat data. Hamburg, Max-Planck-Institut für Meteorologie. (Report 242.)
Paul, F. 2002. Changes in glacier area in Tyrol, Austria, between 1969 and 1992 derived from Landsat TM and Austrian glacier inventory data. Int. J. Remote Sensing, 23(4), 787799.
Porter, S. C. 1975. Equilibrium-line altitudes of Late Quaternary glaciers in the Southern Alps, New Zealand. Quat. Res., 5(1), 2747.
Seltzer, G. O. 1994. Climatic interpretation of alpine snowline variations on millennial time scales. Quat. Res., 41(2), 154159.
Torsnes, I., Rye, N. and Nesje, A.. 1993. Modern and Little Ice Age equilibrium-line altitudes on outlet valley glaciers from Jostedalsbreen, western Norway: an evaluation of different approaches to their calculation. Arct. Alp. Res., 25(2), 106116.
Van der Veen, C. J., ed. 1997. Calving glaciers: report of aWorkshop, February 28–March 2,1997, Columbus, OH. Byrd Polar Res. Cent. Rep. 15.
Van de Wal, R. S.W. and Oerlemans, J.. 1995. Response of valley glaciers to climate change and kinematic waves: a study with a numerical ice-flow model. J. Glaciol., 41(137), 142152.
Veatch, F. M. 1969. Analysis of a 24-year photographic record of Nisqually Glacier, Mount Rainier National Park, Washington. U.S. Geol. Surv. Prof. Pap. 631.
Williams, R. S. Jr and Ferrigno, J. G.. 1993. Glaciers of Europe. Satellite image atlas of glaciers of the world. U.S. Geol. Surv. Prof. Pap. 1386-E.