The World Glacier Monitoring Service (WGMS) collects and publishes mass-balance data of glaciers obtained by direct glaciological and geodetic methods as a contribution to the Global Terrestrial Network for Glaciers (GTN-G) which is part of the Global Climate/Terrestrial Observing Systems (GCOS/GTOS; GCOS 2004). The GTN-G monitoring strategy uses glacier observations in a system of tiers (cf. Reference Haeberli and HoelzleHaeberli and others, 2000; Haeberli, 2004). These tiers include extensive glacier mass-balance measurements within major climatic zones for improved process understanding and the calibration of numerical models (tier 2), as well as for the determination of regional volume changes within major mountain ranges using cost-saving methodologies (tier 3). Glacier-front variation and global inventories are further components of the monitoring strategy but are not discussed here. First surveys of accumulation and ablation of snow, firn and ice at individual stakes date back to the end of the 19th century and beginning of the 20th century, for example at Rhone glacier (Reference Larsen, Motyka, Arendt, Echelmeyer and GeisslerMercanton, 1916) and Silvretta glacier (Reference Hoelzle, Haeberli, Dischl and PeschkeHuss and others, 2008a). Annual glacier mass-balance measurements made by the direct glaciological method (cf. Reference OhmuraØstrem and Stanley, 1969; Reference OhmuraØstrem and Brugman, 1991), based on an extensive net of ablation stakes, snow pits and snow probing, were initiated in 1945 at Storglaciären, Sweden (Reference Hock and JensenHolmlund and others, 2005). Today, six decades of annual (and partially seasonal) mass-balance data are readily available from the WGMS and have been analyzed in detail by Reference DyurgerovDyurgerov (2002) and widely used for studies of glacier changes (e.g. Reference BraithwaiteBraithwaite, 2002; Reference Dyurgerov and MeierDyurgerov and Meier, 2005) and related questions from hydrology (e.g. Reference Braun, Weber and SchulzBraun and others, 2000), climate change (e.g. Reference MercantonOerlemans and Fortuin, 1992; Reference Francou, P. Wagnon, Mendoza and SicartFrancou and others, 2003; Reference Nakawo, Raymond and FountainOhmura, 2006) or sea-level variation studies (e.g. Reference Karlén, Pytte and ØstremKaser and others, 2006; Reference Østrem and BrugmanRaper and Braithwaite, 2006; Meier and others, 2007). However, due to the specific foci of these works, a sound and integrative discussion of the basic dataset and related issues is often missed out. Here we aim to give a review of the present monitoring network, a spatio-temporal analysis of the available data and discuss important issues related to the monitoring and interpretation of glacier mass-balance data.
Available Mass-Balance Data
The WGMS collects mass-balance data of glaciers annually from the preceding year through its collaboration network of national correspondents and principal investigators. This 1 year retention period allows the investigators time to properly analyze and publish their data before making them available to the scientific community and the wider public. Preliminary data are published annually on the WGMS website (www.wgms.ch) and every 2 years in the Glacier Mass Balance Bulletins (Reference ShumskiyWGMS, 2007, and earlier issues) as well as every 5 years, in full detail, in the Fluctuations of Glaciers (Reference Shumskiy and TolstiikovWGMS, 2008a, and earlier issues). All data are available digitally on request and free of charge. We aim to collect seasonal and annual mass-balance data that are measured according to the direct glaciological method (cf. Reference OhmuraØstrem and Stanley, 1969; Reference OhmuraØstrem and Brugman, 1991), i.e. based on a network of ablation stakes and snow pits distributed over the entire glacier, and inter-/extrapolated to the total area of the (same) glacier independently of meteorological or hydrological measurements. Ideally, these mass-balance measurements determined annually are combined with decadal volume-change assessments from geodetic surveys to reduce method-dependent uncertainties and systematic errors.
For the period 1946–2005, there are 3385 annual mass-balance results from 228 glaciers available through the WGMS (Fig. 1). Additional information on mass-balance vs altitude intervals has been reported for 45% of the annual observations and for 56% of the glaciers. A derived equilibrium-line altitude (ELA) and accumulation-area ratio (AAR) are available for 73% and 82% of the annual values and data series, respectively. Seasonal mass balances have been submitted for 45% of the observation years and for 63% of the glaciers. From the 228 available data series, 120 provide information from the 21st century whereas the remaining data series were interrupted in the past century. The average period of observations per series is 15 years, and 39 glaciers have more than 30 years of measurements. For the hydrological year 2006, preliminary mass-balance data have been reported from 97 glaciers. A temporal overview of the number of data series reported to the WGMS is given in Figure 2. The total number of series is shown together with the data series that are ongoing, continuous and consist of up to 15, 30, 45 and 61 observation years. Of the 228 data series, there are just 31 continuous measurement programmes reaching back to 1970, and only 12 back to 1960. Of these glaciers, the Glacier Mass Balance Bulletin (Reference ShumskiyWGMS, 2007) lists a set of 30 ‘reference’ glaciers with continuous measurement programmes back to 1976 and earlier (see Appendix).
Mass-Balance Series Vs The Global Glacier Distribution
About 90% of the mass-balance series come from the Northern Hemisphere and about 40% from Europe. An overview of the available mass-balance series in comparison with the global distribution of glaciers is given in Table 1. Most with the longest time series are found in the European Alps and Scandinavia, followed by North America and High Mountain Asia, with the earliest observations in the 1940s and 1950s, respectively. In the Canadian Arctic Archipelago and in High Mountain Asia, more than two-thirds of the series were interrupted in the last century. The only long-term observation series from the tropics, at Lewis Glacier (1979–96) on Mount Kenya, results in the high average observation length of the region class ‘Africa, New Guinea, Irian Jaya’. Northern Asia and Siberia, South America and the large ice masses around the two ice sheets in Greenland and Antarctica with few and only short-term data series are strongly underrepresented within the network (in respect of their ice cover).
The size and elevation distributions of the glaciers with available mass-balance observations and the glaciers with detailed information in the World Glacier Inventory (WGI; WGMS, 1989) are given in Tables 2 and 3, respectively. The total area covered by the 228 glaciers with mass-balance series is about 10 000km2. More than 75% of the glaciers have an area of 0.1–10.0km2, covering 2.5% of the total ice cover with mass-balance observations. About 90% of this ice cover comes from the 15 glaciers which are >100km2. The elevation of the terminus of 33%, 60% and 80% of these glaciers is below 1000, 2000 and 3000 ma.s.l., respectively. The glaciers of the WGI, with detailed area information, cover overall about 180 000km2 which corresponds to about 23% of the estimated global glacier cover (Reference Dyurgerov and MeierDyurgerov and Meier, 2005). About 80% of these glaciers are <1km2 and represent 7.5% of the inventoried ice cover, whereas 50% of this ice cover comes from the 235 glaciers >100km2. Half of the inventoried glaciers end between 3000 and 5000ma.s.l., and less than 10% reach below 1000ma.s.l.
Spatio-Temporal Analysis of Global and Regional Glacier Mass Changes
The 30 ‘reference’ glaciers (see WGMS, 2007) with (almost) continuous measurements since 1976 show an (arithmetic) average annual mass loss of 0.58mw.e. for the decade 1996–2005, which is more than twice the loss rate for the period 1986–95 (0.25mw.e.), and more than four times the rate for the period 1976–85 (0.14mw.e.). The preliminary value of –1.30mw.e. for 2006 (based on so far reported values from 27 reference glaciers) seems likely to become a new record loss surpassing the mass-balance values of the years 2003 (–1.24mw.e.), 2004 (–0.74mw.e.) and 1998 (–0.71mw.e.). The apparent advantage of analyzing only continuous mass-balance series comes with the disadvantage of having to deal with a small sample size and shorter data series. As a consequence, the available mass-balance values can be averaged annually, disregarding the completeness and length of the time series, but it is necessary to accept all the problems related to moving sample size (Reference Nakawo and RanaOhmura, 2004; Reference Cogley and AndersonCogley, 2005; Reference Dyurgerov and MeierDyurgerov and Meier, 2005). A third approach takes into account that the mean of all glaciers is strongly biased by the large proportion of Alpine and Scandinavian glaciers and, hence, calculates the annual mass balance from the (arithmetic) averaged regional mean values (cf. Haeberli and others, 2007). Other ways to correct for the spatial bias of the data are to use regional averages which are weighted by the corresponding glacier area or to apply polynomial spatial interpolation (e.g. Reference Cogley and AndersonCogley, 2005). Figure 3 shows the cumulative mean (annual) mass balances of: (1) the 30 reference glaciers (and the subsample before 1977); (2) all available glaciers; (3) all available glaciers excluding the reference glaciers; (4) the (subsample of) averages of the 12 macro-regions (Table 1); and (5) the (subsample of) area-weighted averages of the 12 macro-regions (Table 1). The samples used in approaches (1) and (3) are entirely independent; however, (2), (4) and (5) contain the sample from (1). From 1946 to 2006 (i.e. 61 years), the mean of all five approaches shows a cumulative ice loss of 21.2 mw.e., corresponding to an average annual ice loss of 0.35 mw.e. After six decades the difference between the five approaches, i.e. the range between (3) and (5), is 5.14 mw.e. This corresponds to an annual difference of 0.08 mw.e. and is in the same order of magnitude as the mean annual standard errors of the approaches (e.g. (1): 0.19mw.e.; (2) and (3): 0.12 mw.e.).
Regional and individual cumulative mass-balance curves have been shown and discussed in detail, for example by Reference Karlén, Pytte and ØstremKaser and others (2006, and references therein) and WGMS (2007). Here, we dissolve the global glacier mass balances, using approach (2), in six macro-regions and six decades (Fig. 4a–f). In addition to the global and regional mass balance, the standard deviation and the number of available observations are given in order to provide a measure for the variance and the observation density, respectively. For the first three decades (1946–75), mass-balance measurements are available only for latitudes higher than 30°N. Later, information becomes available from the regions between 30°N and 30° S. Latitudes below 30° S are covered with observations from South America after 1976, and from the northeastern side of the Antarctic Peninsula and New Zealand in the last decade. The regional data samples higher than 30°N consist of more than 100 observations per decade after 1976, whereas the six regions south of that exceed ten observations only in the case of South America in the last decade and, hence, are to be considered of limited significance.
In the decade 1946–55, observations of the five glaciers in the European Alps and two glaciers in Scandinavia indicate a mass loss of more than 0.4mw.e., whereas the mean of the two glaciers in North America features a zero balance over the last 3 years of that decade. In the following two decades (1956–75), North America, Europe and Northern Asia show a moderate ice loss, with a mean annual change of –0.10 to –0.25mw.e. After 1976, the decadal mass losses increase in North America and Northern Asia to mean annual mass balances of –0.65 and –0.43mw.e., respectively, in the last decade. The decadal mean balances in Europe remain moderately negative between 1976 and 1995, due to the temporal regain in mass of some glaciers, and become strongly negative (–0.72mw.e.) in the last decade. In South America, positive and close to steady-state balances are reported in the third decade (1976–85) from Quelccaya (PE) and Echaurren Norte (CL), respectively, followed by two decades of average ice losses with values between –0.26 and –0.90mw.e. from a sample of 12 glaciers (including four observations from Ventorillo (MX)). In Africa, negative mass balances have been reported from Lewis Glacier on Mount Kenya with observations between 1979 and 1996. The three negative decadal balances of Southern Asia come from 11 observations from AX010 and Rikha Samba, Nepal Himalaya, and from Changmekhangpu, Sikkim Himalaya. Furthermore, the average mass loss of Bahía del Diablo on Vega Island, Antarctica, (–0.25mw.e.) and the 1 year gain of Brewster (NZ) show up in the last decade.
Available dataset and ‘reference’ glacier concept
Unlike meteorological data, which are mainly measured, collected and made available through governmental agencies (GCOS, 2003), glacier mass-balance observations are usually carried out within scientific projects, collected within a cooperation network, and made available by the WGMS. As a consequence, the data collection runs with a very low funding level and depends fully on the cooperativeness of the individual investigators. An exception is Norway where for decades the majority of mass-balance observation programmes (over 40 series) have been run and made available by the Norwegian Water Resources and Energy Directorate. At first glance, the available mass-balance observations seem to be well distributed over the globe. However, looking at the length and continuity, or lack thereof, of the time series it is apparent that the observation network is strongly overrepresented by glaciers located in Europe and North America. Unfortunately, the regions hosting large proportions of the glacier cover of the Earth (e.g. the Canadian Archipelago, Patagonia, the Arctic Islands and the ice bodies around the two ice sheets in Greenland and Antarctica) are greatly underrepresented or even lack any long-term data series. In Asia a large number of observation series have been started but the vast majority were interrupted in the 1980s and 1990s.
With regard to the limited dataset, it is particularly unfortunate that data on mass balance vs altitude have been reported for only 45% of the annual measurements and ELA/AAR for only 73%. The seasonal balances were measured at 63% of the glaciers. Mass balance vs altitude and ELA/AAR can, at least in principle, be derived from the existing raw data and would provide valuable information about, for example, mass-balance gradients (e.g. Reference Dyurgerov and DwyerDyurgerov and Dwyer, 2001), mass turnover (e.g. Reference DyurgerovDyurgerov, 2002), the glacier climate regime (e.g. Reference Haeberli, Cihlar and BarryHoelzle and others, 2003) and climate sensitivity (e.g. Reference MercantonOerlemans and Fortuin, 1992). In order to improve the available seasonal dataset, the annual observation programmes need to be extended with seasonal field surveys, although the in situ determination of the winter/wet-season accumulation might be a major challenge. Once available, these data would provide insight into the processes behind the glacier fluctuations (e.g. Reference Dyurgerov and MeierDyurgerov and Meier, 1999; Reference Rabus and EchelmeyerSchöner and others, 2000; Reference Nakawo, Raymond and FountainOhmura, 2006) and would be of great value for model validation (e.g. Reference MeierOerlemans, 2001). The observation and standardized compilation of further parameters such as temperature and precipitation at the ELA, ice/firn temperatures or the remaining mean/maximum ice thickness would be of great value for many scientific questions, especially in view of the fast changes in nature.
Trend analyses ideally are based on long-term measurement series. It is for that purpose that the WGMS has introduced the concept of ‘reference’ glaciers. A mass-balance observation programme is considered to be a ‘reference’ series if it is ongoing, continuous and long-term (see Fig. 2), with reliable reporting of data and meta-information. At present, a set of 30 reference glaciers in nine mountain ranges (in five of the macro-regions; see Table 1) are listed in the Glacier Mass Balance Bulletin series (WGMS, 2007) with (almost) continuous observations since 1976, and for 11 of the reference glaciers observations reach back to 1960 or earlier (see Appendix). The set of reference glaciers is reviewed after each call-for-data and it will help to support the continuation of these unique data series. Further glaciers that could be considered next as reference glaciers include Baby and White (Canadian High Arctic), Peyto (Canadian Rocky Mountains), Helm (Canadian Coast Mountains) and Lemon Creek (US Coast Range).
Data provided by the WGMS, but also in general, are subject to errors and inaccuracies and, hence, have to be quality-checked against the related literature before being used in further analysis. The measurement and calculation of mass balances contain various sources of systematic and random errors that include: the accuracy of stake readings, snow probing and snow/firn density measurements (Reference OhmuraØstrem and Brugman, 1991; Reference Holmlund, Jansson and PetterssonJansson, 1999; Reference Ohmura, Sparling and HawkesworthØstrem and Haakensen, 1999), the distribution of the stake and pit network (Reference Kuhn, Dreiseitl, Hofinger, Markl, Span and G.Lliboutry, 1974; Reference CogleyCogley, 1999; Reference Fountain and VecchiaFountain and Vecchia, 1999), the method used to interpolate between the measurement points (Reference Ohmura, Sparling and HawkesworthØstrem and Brugman, 1991; Reference KääbKaser and others, 1996; Reference Funk, Morelli and StahelFunk and others, 1997; Reference Haeberli and HolzhauserHock and Jensen, 1999), the extrapolation for unprobed areas such as crevasse (Reference Jóhannesson, Raymond and WaddingtonKarlén, 1965; Jannson, 1999) or debris-covered zones (Reference LliboutryNakawo and Rana, 1999, Reference Meier and BahrNakawa and others, 2000), special problems related to internal accumulation (Reference O’Neel, Pfeffer, Krimmel and MeierRabus and Echelmeyer, 1998) and calving (Reference Arendt, Echelmeyer, Harrison, Lingle and ValentineArendt and others, 2002; Reference Østrem and StanleyRignot and others, 2003; Reference Oerlemans and FortuinO’Neel and others, 2005), as well as any changes in the glacier area that are usually not subject to annual measurement (Reference Elsberg, Harrison, Echelmeyer and KrimmelElsberg and others, 2001). Generally, the accuracy of the results from the summer balance and the ablation area exceeds that of the results from the winter balance and the accumulation zone, as the ablation processes spatially correlate much better than the snow accumulation (Reference Holmlund, Jansson and PetterssonJansson, 1999; Reference KrimmelKuhn and others, 1999). The dimensions of the various errors are glacier-specific, and their propagation is hard to quantify. Typical estimates of the annual mass-balance accuracy range lie between 0.1 mw.e. (Reference Holmlund, Jansson and PetterssonJansson, 1999) and 0.6mw.e. (Reference Funk, Morelli and StahelFunk and others, 1997). Further errors are introduced when comparing mass balances derived from different measurement systems (stratigraphic vs fixed-date systems) or from different hemispheres (i.e. the hydrological years of the two hemispheres are shifted by half a year).
In contrast to the in situ point measurements based on a changing reference (i.e. the previous year’s summer surface), geodetic methods have the advantage of providing volume changes over the entire glacier based on a non-changing reference (i.e. the surrounding bedrock), but encounter some problems with (fresh) snow in the accumulation area and density assumption for snow, firn and ice (Reference Kaser, Munari, Noggler, Oberschmied and ValentiniKrimmel, 1999). A combination of the direct glaciological method with decadal geodetic methods has proven to be an appropriate way to detect systematic biases (e.g. Reference Holmlund, Jansson and PetterssonJansson, 1999; Reference Kaser, Munari, Noggler, Oberschmied and ValentiniKrimmel, 1999; Reference KrimmelKuhn and others, 1999; Reference Bauder, Funk and HussBauder and others,2007). Several studies have shown that other remote-sensing methods are also useful for determining glacier volume changes over longer periods, for example those based on airborne laser altimetry (Reference Arendt, Echelmeyer, Harrison, Lingle and ValentineArendt and others, 2002) or the differencing of digital elevation models using SPOT5 (Systéme Probatoire pour l’Observation de la Terre) (Reference Berthier, Arnaud, Kumar, Ahmad, Wagnon and ChevallierBerthier and others, 2007), SRTM (Shuttle Radar Topography Mission) (Reference Østrem and StanleyRignot and others, 2003; Reference KuhnLarsen and others, 2007; Reference Paul and HaeberliSchiefer and others, 2007) or ASTER (Advanced Spaceborne Thermal Emission and Reflection Radiometer) and its optical stereo (Reference Huss, Bauder and FunkKääb, 2007). Differencing of digital elevation models from space-borne sensors is, at present, no replacement for such validation purposes, as the resolution and quality is still too low. However, it can already be used to assess the representativeness of available field measurements with respect to long-term glacier volume changes over large areas (e.g. Reference JanssonKääb, 2008).
Given the present state of knowledge, the manifold sources of errors, and the often missing meta-information, it is not possible to quantify the overall error of the available mass-balance data. As an objective measure of uncertainty we, hence, estimate the confidence interval of the here-presented global dataset to be in the magnitude of two standard errors of the reported annual mass-balance data, i.e. between 0.25 and 0.5 mw.e. but increasing with small sample sizes. With an average annual standard deviation of about 0.7mw.e., the annual signal of the mean mass balances is smaller than the regional variability. However, the change of the cumulative mean mass balances over extended periods is far beyond the estimated uncertainty and, hence, significant. In order to reduce, better understand and quantify these uncertainties, we reiterate the recommendation to validate and calibrate systematically the mass-balance data derived from annual/seasonal field observations with decadal mass/volume changes as assessed from geodetic methods, and to submit also the corresponding meta-information. It is one of the major shortcomings of the available datasets that this information is often missing and, as a consequence, it is difficult to reconstruct whether a mass-balance series is geodetically calibrated or not.
Global and regional glacier changes
The global average cumulative mass balance indicates a strong mass loss in the first decade after the start of the measurements in 1946 (for Europe), slowing down in the second decade (1956–65; based on observations above 30°N only), followed by a moderate ice loss between 1966 and 1985 (with data from the Southern Hemisphere only since 1976) and a subsequent acceleration until the present (2006). Over the six decades the global average ice loss has cumulated to about 21 mw.e., or an average mass balance of –0.35mw.e. a–1, a dramatic ice loss compared with the global average ice thickness which is estimated (by dividing volume by area) to be between 100 m (Reference Rignot, Rivera and CasassaSolomon and others, 2007) and about 180m (personal communication from A. Ohmura, 2008). The vast ice loss since the mid-20th century has already led to the disintegration of many glaciers within the observation network, including Lower Curties and Columbia 2057 (US), Chacaltaya (BO), Caréser (IT), Lewis (KE) and Urumqihe (CN), and presents one of the major challenges for glacier monitoring in the 21st century.
The cumulative mass loss is in the same order for all the sample/methods used for averaging. These methods are: 30 reference glaciers (and subsamples before 1976), all but the reference glaciers, all glaciers, arithmetic and area-weighted average of the regions. In addition, the global result corresponds well with the findings of several approaches as summarized by Reference Karlén, Pytte and ØstremKaser and others (2006) and published in Reference Rignot, Rivera and CasassaSolomon and others (2007). They state an average annual ice loss of 0.283 mw.e. between 1961 and 2004 (in this study, sampling approach (2) gave 0.286mw.e.). This is not surprising, as both approaches are based mainly on the dataset provided by the WGMS, but Reference Karlén, Pytte and ØstremKaser and others (2006) extended about 70 data series which have either not been reported to the WGMS or had to be rejected due to insufficient data quality or to their dependence on hydrological/meteorological measurements. For the last three decades (1976–2006) the strong cumulative changes calculated from a sample of 60–100 glaciers from 8 of the 12 glacierized macro-regions can be considered as representative of the global mass-balance signal, as the mass-balance variability in time is spatially correlated over distances of several hundred kilometres (Reference Kuhn, Markl, Kaser, Nickus, Obleitner and SchneiderLetréguilly and Reynaud, 1990; Reference Cogley, J.G. and AdamsCogley and Adams, 1998). Some reservation has to be made because of the lack of data series for the glaciers around Greenland and Antarctica as well as in New Zealand, and because the sample of glaciers with mass-balance observations might not be representative, with respect to characteristics (e.g. tidewater, debris cover), ice temperature and hypsometry, of the ice bodies that are major (potential) contributors to sea-level change. Going further back in time, the average is based only on observations from North America, Europe and Asia, increasingly dominated by the sample from Europe. Together with the less pronounced trends, the ‘global’ average mass balance from this small and regionally biased sample might not be representative any-more for the other regions. Approaches using regional (area-weighted) averaging might reduce this bias but cannot completely overcome the lack of observations.
The regional cumulative mass changes can be considered to be a good first-order index for the glacier changes in these areas over the corresponding period. However, further analysis might have to resample the series according to the focus of the study and have to consider the glacial characteristics of the different geographical regions (Reference Fountain and VecchiaFountain and Vecchia, 1999): mass-balance variability of temperate alpine glaciers in the Northern Hemisphere is mainly driven by accumulation and ablation during winter and summer seasons, respectively (Reference KrimmelKuhn and others, 1999). This is not valid for low-latitude glaciers where ablation occurs throughout the year and multiple accumulation seasons exist (Reference KääbKaser and Osmaston, 2002), monsoonal glaciers of the Himalaya where the accumulation and ablation seasons occur at about the same time (Reference Ageta and FujitaAgeta and Fujita, 1996), and high-altitude and polar glaciers where any season can be the accumulation season (Reference ChinnChinn, 1985). However, annual/net mass balances integrate these seasonal complexities and as such might not be relevant for some analysis. Strongly diverse mass-balance characteristics also exist between glaciers under dry-continental conditions and in maritime regions. The former can be found in northern Alaska, Arctic Canada, sub-Arctic Russia, parts of the Andes near the Atacama Desert and in many central Asian mountain chains. They show low mass turnover, equilibrium lines at high altitudes and firn/ice well below melting temperatures. The latter, glaciers in maritime regions, are found in Patagonia, Iceland, southeast Alaska, as well as in the coastal parts of Norway and New Zealand, and show high mass turnover, low equilibrium lines and firn/ice at melting temperatures (Reference Raper and BraithwaiteShumskiy, 1964; Reference Oerlemans and BalkemaOhmura and others, 1992; Haeberli and Burn, 2002).
For a sound analysis of global or regional glacier changes, results from the valuable but small sample of mass-balance measurements have to be complemented with other glacier observations, such as front variation measurements (available for about 1700 glaciers worldwide; WGMS, 2008, and earlier issues) or investigations of the glacier retreat from the trimlines of the Little Ice Age which can be found around the globe (Grove, 2004). Examples for such integrative analysis for entire mountain ranges are given by Reference Letréguilly and ReynaudMolnia (2007) for Alaska, by Reference CasassaCasassa and others (2007) for the Andean glaciers, by Reference KääbKaser and Osmaston (2002) for tropical glaciers, by Reference Andreassen, Elverøy, Kjøllmoen, Engeset and HaakensenAndreassen and others (2005) for Norway, by Reference SolomonZemp and others (2007b) for the European Alps, by Reference Kaser and OsmastonKotlyakov (2006) for Russia, and by Reference ChinnChinn (2001) for New Zealand, as well as by Reference Haeberli, Cihlar and BarryHoelzle and others (2003), Grove (2004), Reference VoldenZemp and others (2007a), and WGMS (2008b) for a global overview.
Glacier mass balance as a climate proxy
Glacier mass changes are used widely as a climate proxy in many environmental and climate change assessment reports (e.g. Reference Rignot, Rivera and CasassaSolomon and others, 2007). The mass change of glaciers, which are not influenced by thick debris covers, calving or ‘surge’ instabilities, is the direct, undelayed reaction of a glacier to climatic forcing. The mass-balance variability in time is well correlated spatially over larger distances (Reference Kuhn, Markl, Kaser, Nickus, Obleitner and SchneiderLetréguilly and Reynaud, 1990; Reference Cogley, J.G. and AdamsCogley and Adams, 1998) and with climatic parameters such as (seasonal) air temperature, precipitation and sunshine duration (Reference Kuhn, Dreiseitl, Hofinger, Markl, Span and G.Lliboutry, 1974; Reference Rabus and EchelmeyerSchöner and others, 2000). However, the glacier mass-balance change provides an integrative climatic signal, and the quantitative attribution of the forcing to individual meteorological parameters is not straightforward. The energy and mass balance at the glacier surface is influenced by changes in atmospheric conditions (solar radiation, air temperature, precipitation, wind, cloudiness, etc.). Air temperature thereby plays a predominant role, as it is related to the radiation balance, turbulent heat exchange and solid/liquid precipitation ratio (Reference Kaser, Cogley, Dyurgerov, Meier and OhmuraKuhn, 1981; Reference MolniaOhmura, 2001). The climatic sensitivity of a glacier not only depends on regional climate variability but also on local topographic effects, which can result in two adjacent glaciers featuring different specific mass-balance responses (Reference KotlyakovKuhn and others, 1985).
For a temperate glacier, an assumed step-change in climatic conditions would cause an initial mass-balance change followed by a return towards zero values, due to the adaptation of the size of the glacier (surface area) to the new climate (Reference Huss, Baude, Funk and HockJóhannesson and others, 1989; Reference Haeberli, Bamber and PayneHaeberli and Hoelzle, 1995). The observed trend of increasingly negative mass balance over reducing glacier surface areas thus leaves no doubt about the ongoing climatic forcing resulting from the change in climate and possible enhancement mechanisms such as mass-balance/altitude feedback, altered turbulent fluxes due to the size and existence of rock outcrops or changes in the surface albedo. The specific mass-balance data can be compared directly between different glaciers of any size and elevation range. These data series provide a combined hydrological and climatic signal. Glacier contribution to runoff can be calculated very easily by multiplying the specific mass balance with the corresponding glacier area, whereas a climatic interpretation needs to relate the mass changes to a glacier reference extent in order to derive a pure climate signal (cf. Reference Elsberg, Harrison, Echelmeyer and KrimmelElsberg and others, 2001; WGMS, 2007, and earlier issues).
In order to advance the present mass-balance monitoring, as part of the Global Climate Observing System, it is necessary to continue with the existing data series and to extend the network in respect of the global glacier distribution. Systematic use should be made of remote sensing and geo-informatics for assessing the representativeness of the available in situ observations (e.g. Reference Ohmura, Kasser and FunkPaul and Haeberli, 2008). Analytical or numerical modelling is needed to quantify the above-mentioned topographic effects as well as to attribute the glacier mass changes to individual meteorological or climate parameters (e.g. Reference Kaser, Cogley, Dyurgerov, Meier and OhmuraKuhn, 1981; Reference MeierOerlemans, 2001) and, in combination with measured and reconstructed glacier fluctuations, to compare the present mass changes with the (pre-)industrial variability (e.g. Reference Haeberli, Burn and SidleHaeberli and Holzhauser, 2003)
Glacier mass balance and sea-level changes
Measurements from the mass-balance monitoring network are used to estimate the contribution of glaciers to past, present and future sea-level changes (Reference Karlén, Pytte and ØstremKaser and others, 2006 and references therein; Reference Østrem and BrugmanRaper and Braithwaite, 2006). These estimates are hampered by the fact that: (a) no complete detailed inventory of the Earth’s glaciers exists; (b) the estimation of the overall ice volume of glaciers contains large uncertainties; (c) the spatial distribution of the available mass-balance series is disproportionate to the global ice cover; and (d) the small sample of mass-balance observations is (most probably) not representative for the entire sample of glaciers.
Most of the approaches use the regional ice extents of Reference Dyurgerov and MeierDyurgerov and Meier (2005 and earlier versions; mainly based on WGMS, 1989) as a baseline inventory to calculate the overall potential sea-level rise equivalent as well as sea-level changes. A detailed inventory, including information on glacier location, size and altitude extent, is only available for about 70 000 glaciers covering about 180 000km2. This corresponds to only 43% of the approximate total number and 23% of the overall glacier area based on rough estimates from Meier and Reference Bahr, Meier and PeckhamBahr (1996) and Reference Dyurgerov and MeierDyurgerov and Meier (2005), respectively. As a further uncertainty factor, the existing inventory contains no information on the proportion of ice below sea level. There are only a few glaciers where thickness measurements have been carried out (so far not compiled by the WGMS). Different approaches exist for estimating the overall ice volume (e.g. Reference Haeberli, Bamber and PayneHaeberli and Hoelzle, 1995; Reference Bahr, Meier and PeckhamBahr and others, 1997), but they all contain a number of uncertainties that could amount to 30–50% of the total ice volume. The latest assessment report of the Intergovernmental Panel on Climate Change (IPCC) (Reference Rignot, Rivera and CasassaSolomon and others, 2007) quotes the total area and corresponding potential sea-level rise as 510 000–540 000km2 and 150–370 mm, respectively. These estimates, as correctly noted in Reference Rignot, Rivera and CasassaSolomon and others (2007), do not include ice bodies around the ice sheets in Greenland (70 000km2 based on Reference Schöner, Auer and BöhmWeidick and Morrris, 1998) and Antarctica (169 000 km2 based on Reference RaupShumskiy, 1969) and, hence, might considerably underestimate the overall potential sea-level rise due to melting glaciers. As shown above, many of the regions with large ice covers, such as the Canadian Arctic, High Mountain Asia, South America and around the two ice sheets, are not represented by an adequate number of long-term mass-balance measurements. Mass-balance programmes require intensive fieldwork and are usually carried out on glaciers that are easy accessible, safe and not too large. Hence, these glaciers are neither representative of the glacier size distribution nor of the elevation distribution of all glaciers, at least when compared with the presented data of about 70 000 glaciers with detailed inventory information (the data for an exact comparison are not available).
The current first-order estimates of the contribution from glaciers to past, present and future sea-level changes can only be improved significantly by completing a detailed baseline inventory of the Earth’s glaciers as well as a review and enlargement of the available (measured) glacier thickness dataset. This would be needed to scale-up the few in situ series that we have to cover all glaciers. It is hoped that internationally coordinated efforts, such as the European Space Agency-funded GlobGlacier project (Reference Schiefer, Menounos and WheateVolden, 2007) or the Global Land Ice Measurements from Space initiative (Reference Østrem and HaakensenRaup and others, 2007), will make major steps in that direction. Furthermore, it is necessary to continue and extend the present mass-balance network in respect of the global distribution of the ice cover and to make systematic use of remote sensing and geo-informatics to assess the representativeness of the available in situ annual mass-balance series (Reference Ohmura, Kasser and FunkPaul and Haeberli, 2008) as well as of decadal volume changes in ice fields and ice caps that are too large for in situ measurements (e.g. Reference Østrem and StanleyRignot and others, 2003; Reference KuhnLarsen and others, 2007).
Conclusions and Consequences for The Monitoring of Glaciers
During the six decades of glacier mass-balance observation, the WGMS has compiled a dataset of more than 3400 annual mass-balance measurements from 228 glaciers worldwide. The collection and free availability of these data through a purely scientific collaboration network is a great success and, at the same time, one of the reasons why the present monitoring network is unevenly distributed in comparison with the global ice coverage. The 30 reference glaciers with continuous observation series since 1976 show an accelerated thinning, with mean annual ice losses of 0.14 m w.e. (1976–85), 0.25 m w.e. (1986–95) and 0.58mw.e. (1996–2005), which gives a total average ice-thickness reduction of about 10mw.e. The available data from the first three decades indicate strong mass losses as early as the 1940s and 1950s, followed by moderate mass losses until the end of the 1970s. The mass-balance data are widely used by the scientific community and represent one backbone of glacier research. Mass balance is recognized as an essential climate variable within the global climate-related observing systems and is, in effect, the largest non-steric contributor to the global sea-level rise at the turn of the century.
In view of the discrepancy between the relevance of glacier mass-balance data and the relatively small set of current long-term observations with a strong bias towards the Northern Hemisphere and Europe, we strongly recommend to:
continue the work on the glaciers with long-term measurements,
resume the interrupted long-term series,
replace vanishing glaciers by starting early with parallel observations on larger or higher-reaching glaciers,
extend the monitoring network to strategically important regions with few or no series,
validate and calibrate the field measurements with the results of geodetic methods (and clearly flag such data series),
make systematic use of remote sensing and geoinformatics for assessment of the representativeness of the field measurements and for (decadal) volume change analysis in mountain ranges lacking such data, and
make the data and related meta-information readily available to the scientific community and the wider public.
The potential dramatic changes as sketched for the 21st century by the IPCC report (Reference Rignot, Rivera and CasassaSolomon and others, 2007) require critical reflection and a rigorous implementation of the monitoring strategies for glaciers in order to face the challenges of the fast changes in nature and to bridge the gap between historical observation series and the new technologies.
We express our gratitude to the National Correspondents, Principal Investigators of the WGMS and all other individuals and institutions for the collection and free exchange of these important data over many years. Sincere thanks go to the field scientists and programme managers of the long-term data series: without their commitment to maintain the observation series through politically and financially difficult times, many of the current insights would not have been possible. We thank F. Paul for thoughtful remarks on a first draft of the manuscript, the two reviewers for constructive feedback, the scientific editor for providing guidance towards the final paper, and I. Woodhatch for polishing the English. This study has been mainly funded by the Swiss National Science Foundation and the Swiss Federal Office for the Environment.