# Mass balance of Glaciar Zongo, Bolivia, between 1956 and 2006, using glaciological, hydrological and geodetic methods

Alvaro Soruco, Christian Vincent, Bernard Francou, Pierre Ribstein, Thierry Berger, J.E. Sicart, Patrick Wagnon, Yves Arnaud, Vincent Favier, Y. Lejeune

► To cite this version:

Alvaro Soruco, Christian Vincent, Bernard Francou, Pierre Ribstein, Thierry Berger, et al.. Mass balance of Glaciar Zongo, Bolivia, between 1956 and 2006, using glaciological, hydrological and geodetic methods. Annals of Glaciology, 2009, 50 (50), pp.1-8. 10.3189/172756409787769799. insu-00413389

HAL Id: insu-00413389

https://insu.hal.science/insu-00413389

Submitted on 7 May 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Mass balance of Glaciar Zongo, Bolivia, between 1956 and 2006, using glaciological, hydrological and geodetic methods

A. SORUCO,<fnref n="1" />,<fnref n="2" /> C. VINCENT,<fnref n="3" /> B. FRANCOU,<fnref n="1" />,<fnref n="2" /> P. RIBSTEIN,<fnref n="4" /> T. BERGER,<fnref n="2" />
J.E. SICART,<fnref n="2" /> P. WAGNON,<fnref n="1" />,<fnref n="2" /> Y. ARNAUD,<fnref n="1" />,<fnref n="2" /> V. FAVIER,<fnref n="3" /> Y. LEJEUNE<fnref n="5" />

E-mail: alvaro.soruco@gmail.com

ABSTRACT. The longest continuous glaciological mass-balance time-series in the intertropical zone of South America goes back to 1991 on Glaciar Zongo, Bolivia. Photogrammetric and hydrological data have been used to (1) check the specific net balance over long periods and (2) extend the mass-balance time series over the last 50 years. These data reveal a bias in the glaciological mass balance which can be explained by the field-measurement sampling network. Our study shows a large temporal variability of the surface mass balances in the ablation area and reveals strong relationships between independent surface mass-balance data coming from selected ablation areas with numerous data. It demonstrates the very large contribution (80%) of low-elevation ranges (one-third of the surface) to the specific mass balance and, consequently, the importance of the reduction of the area of the tongue. With these new results, Glaciar Zongo offers the longest and most accurate mass-balance series in any Andean country. The dataset shows that Glaciar Zongo experienced a relatively steady state over the period 1956–75, with even a slight mass gain over 1963–75, and a rapid and continuous decrease since then.

## 1. INTRODUCTION

As shown in numerous studies, glacier mass balances are good indicators of climate change in tropical areas where historical climate data are rare (Dyurgerov, 2002; Coudrain and others, 2005; Solomon and others, 2007). According to the astronomical, thermal, and hygric delimitations proposed by Kaser (1999), approximately 99% of tropical glaciers are located in the highlands of the intertropical belt of South America. The others are located in East Africa and New Guinea (Hastenrath, 1984, 1995; Kaser and others, 1996; Kaser, 1999).

Moreover, mass-balance studies of glaciers in South America are essential to estimate the quantity of water available for future human consumption and hydropower generation. The assessment of water resources can help minimize the economic impact of the retreat of these glaciers (Vergara and others, 2007). Consequently, knowledge of the climatic parameters controlling glacier mass balances is of prime importance to understand the response of tropical glaciers to the specific climate forcing over the intertropical zone. For this purpose, long-term mass-balance data are needed.

In Bolivia, a long-term glaciological program has been carried out since 1991 by the Institut de Recherche pour le Développement (IRD), France, (Francou and others, 1995; Ribstein and others, 1995) and in the framework of the French glacier monitoring program GLACIOCLIM (http://www-lgge.obs.ujf-grenoble.fr/ServiceObs/index.htm). The objective of this program was to survey two glaciers: Glaciar Zongo and Glaciar Chacaltaya. Given the imminent disappearance of Glaciar Chacaltaya (Ramirez and others, 2001), another glacier monitoring program was started on Glaciar Charquini

Sur (Rabatel, 2005; Rabatel and others, 2006) in 2002. These glaciers were selected mainly because of the influence of ice melting on water resources for human consumption and hydropower generation for the city of La Paz.

Two independent methods have been used to assess these glacier mass balances using glaciological data (Francou and others, 1995) and hydrological data (Ribstein and others, 1995). In addition, numerous studies have been carried out to assess the relative importance of melting and sublimation in the energy balance (Wagnon and others, 1999a) in order to understand the relationships between energy balance and glacier runoff (Wagnon and others, 1999b) and to provide insight into the atmospheric control of the energy balance (Sicart and others, 2005). Moreover, hydrological investigations have been performed to compare mass-balance results from the glaciological and hydrological methods (Sicart and others, 2007).

Differences between the glaciological and hydrological methods were first pointed out by Francou and others (1995) and attributed to ice sublimation in the dry season. Later, this assumption was rejected by Sicart and others (2007), who assumed that sublimation was too small to explain the discrepancies. Sicart and others (2007) concluded that the hydrological method gives values that are approximately $0.6 \ \text{m w.e.} \cdot a^{-1}$ lower than the glaciological method. According to these authors, the large difference between the two methods was due to the catch deficiency of rain gauges and the lack of precipitation measurements at high altitudes around Glaciar Zongo.

However, neither method can estimate accurately the overall accumulation rate in the accumulation zone of the glacier. On the one hand, the glaciological mass balance is

<footnote n="1">¹IRD-Great Ice, Laboratoire de Glaciologie et Géophysique de l'Environnement du CNRS, BP 96,
38420 Saint-Martin-d'Hères Cedex, France</footnote>

<footnote n="2">²IRD-Great Ice, Maison des Sciences de l'Eau, BP 64501, 34394 Montpellier Cedex 5, France</footnote>

<footnote n="3">³Laboratoire de Glaciologie et Géophysique de l'Environnement du CNRS (associé à l'Université Joseph Fourier-Grenoble I),
54 rue Molière, BP 96, 38402 Saint-Martin-d'Hères Cedex, France</footnote>

<footnote n="4">⁴UMR Sisyphe, Université Pierre et Marie Curie, case 123, 4 Place Jussieu, 75252 Paris Cedex 05, France</footnote>

<footnote n="5">⁵Météo-France, 1441 rue de la Piscine, 38400 Saint-Martin-d'Hères Cedex, France</footnote>
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Contour map of Glaciar Zongo showing the surface topography in 1983 with 20 m contour intervals and the terminus of the glacier in 2006.</figcaption></figure>

calculated from only two or three snow pits. On the other hand, the highest rain gauge used to calculate the hydrological mass balance is only 5200 m a.s.l. and the accumulation rate inferred from these data can be questioned. Moreover, the accumulation zone corresponds to about 60% of the total area. Consequently, the uncertainty for both methods is large (Sicart and others, 2007) and these data cannot be used to accurately determine the mass balance of Glaciar Zongo.

The first purpose of our study is to assess the cumulative glacier mass balance from an independent method using photogrammetric data. This geodetic method has been widely used to check cumulative mass balances (Vincent, 2002; Bauder and others, 2007; Thibert and others, 2008). These new data make it possible to compare glacier mass balances from the glaciological method (1991–2006), the hydrological method (1974–2006) and the geodetic method (1956–2006). The second purpose is to extend the mass-balance time series into the past using both hydrological and photogrammetric data.

is temperate (Francou and others, 1995) and covers an area of 1.96 km² over a hydropower catchment area of 3.3 km². Its length is 3000 m and its width is 750 m. The glacier flows from 6100 to 4900 m a.s.l. The average ice-flow velocity is 20 m a⁻¹ between 4900 and 5200 m a.s.l.

Meteorological conditions in the intertropical belt have been thoroughly analyzed by Kaser and others (1996) and Kaser (2001), and two distinct areas can be distinguished: (i) inner tropics, with almost no seasonality in precipitation and (ii) outer tropics, with strong seasonality in precipitation. These meteorological differences have a strong impact on glacier surface mass balance. In the inner tropics, both accumulation and ablation occur throughout the year. In the outer tropics, accumulation and ice/snow melt occur mainly during the wet/summer season. During the dry/winter season, ice/snow melt is reduced and ice/snow sublimation is significant (Kaser, 2001; Favier and others, 2004). Glaciar Zongo is located in the outer tropics area (Wagnon and others, 1999a; Favier and others, 2004). Average precipitation, temperature and runoff are shown in Figure 2.

## 2. LOCATION

Glaciar Zongo is located in the Huayna Potosi region (Cordillera Real, Bolivia; 16° S, 68° W), 30 km north of La Paz, between the dry Altiplano plateau in the west and the wet Amazonian basin in the east (Fig. 1). Most of the precipitation comes from the Amazonian basin. The glacier

## 3. DATA AND METHODS

The mass balances of Glaciar Zongo have been reconstructed since 1956 using glaciological, hydrological and geodetic data. According to the wet/dry seasons, the hydrological year has been fixed between 1 September and 31 August (Ribstein and others, 1995). Only a short
description of each method is given since the methods have been thoroughly described in previous papers (Ribstein and others, 1995; Sicart and others, 2007).

The glaciological method estimates the loss or gain in mass over the whole glacier surface from field measurements obtained from stakes in the ablation area and snow pits/drilled cores in the accumulation area (Fig. 1). The specific mass balance of the glacier is obtained using the traditional glaciological method (Paterson, 1994). Annual mass balances have been calculated with this method over Glaciar Zongo since 1991, using 14 stakes for monthly measurements and 3 snow pits/drilled cores for yearly measurements averaged over the ablation and accumulation zones respectively (http://glaciog.ujf-grenoble.fr/ServiceObs/ SiteWebAndes/index.htm). Provided that snow accumulation can occur at any time during the year in the ablation area, snow height measurements are required together with stake emergence measurements. Unfortunately, snow height in the ablation area was not always measured between 1991 and 2001, which could lead to greater uncertainties in annual mass balance over this period. For these measurements in the ablation area, the snow and ice density were assumed to be 400 and 900 kg m⁻³ respectively (Dyurgerov, 2002). In the accumulation area, snow-pit/drilled-core measurements were carried out at least once a year at the end of the hydrologic year. The snow density was measured for each pit. Until now, the surface area used to calculate the specific mass balance of Glaciar Zongo came from a 1983 map (Fig. 1) without taking into account the decrease in surface area (Francou and others, 1995; Sicart and others, 2007). Consequently, it was assumed that the influence of area reduction since 1983 was not significant. As can be seen below, this assumption is not supported by photogrammetric data.

For these reasons, the 15 annual mass balances of Glaciar Zongo (1991–2006) were recalculated for our study. The main changes involve the ablation area. First, the snow heights were roughly estimated in the ablation area for the 1991–2001 period when the snow layer persisted at the end of the hydrological year. These estimates were made using monthly mass-balance measurements (stake emergence measurements). No change was made in the accumulation area. Second, the assumption of constant surface area was rejected and the specific mass balance was computed using the actual surface areas. For this purpose, areas for different altitude ranges were calculated each year using detailed maps from photogrammetric data (see below). These two errors have been corrected, leading to changes in mass balance of 0.11 and 0.23 m w.e. a⁻¹ respectively over the period 1991–2006.

The second method of obtaining mass balance is based on hydrological data. The hydrological mass balance is obtained from precipitation, glacier runoff and sublimation. The hydrological mass balance ($b_{n,H}$) is given by $b_{n,H} = P - 1/S_g [D - (S - S_g)c_e P] - E$ (Sicart and others, 2007), where $S_g$ is the glacier surface area (m²), $S$ is the total surface area of the basin (m²), $c_e$ is the runoff coefficient, $P$ is the precipitation (m a⁻¹), $D$ is the runoff at the outlet of the basin (m³ a⁻¹) and $E$ is the sublimation of the glacier (m a⁻¹). Precipitation and the runoff of Glaciar Zongo were measured over two periods. From 1973 to 1991, precipitation (rain gauge P4750, at 4750 m a.s.l., 1 km away from the glacier tongue with an opening area of 314 cm²) and runoff were measured twice a day (at 0700 and 1500 h) by the Bolivian

<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Climatic conditions of the outer tropical Glaciar Zongo. The black curve shows the mean monthly runoff (1973–2006), the grey curve shows the mean monthly precipitation (P4750; 1971–2006) and the dashed curve shows the mean monthly temperature (El Alto station; 1968–2006).</figcaption></figure>

Hydropower Company (Cobee). When the glaciological program started in 1991, a network of five additional rain gauges (opening area of 2000 cm²) was set up and was read monthly around Glaciar Zongo together with an automated stream gauge at the outlet of the Zongo basin. Comparison between P4750 and the rain gauge network around the glacier shows that rain gauge P4750 underestimates precipitation around the glacier by 20% (Francou and others, 1995) and over the glacier by 29% (Sicart and others, 2007). Additional precipitation measurements, from a Geonor pluviometer set up in 2003, reveal that P4750 underestimates precipitation by 37%. To estimate the hydrological mass balance of Glaciar Zongo, we used the daily rain gauge (P4750) to assess the accumulation rate between 1973 and 2006 and a positive correction of 37%. In addition, a runoff coefficient (ratio of runoff to precipitation) of 0.8 was used for the surface area not covered by the glacier (Ribstein and others, 1995, 1999; Sicart and others, 2007). The water flow is deduced from the stream gauge measurements between 1973 and 1991 and from the automated stream gauge measurements between 1991 and 2006. To estimate monthly sublimation ($E$), we used monthly values of US National Centers for Environmental Prediction (NCEP)-US National Center for Atmospheric Research (NCAR) re-analysis data (Kalnay and others, 1996) from the gridcell closest to Glaciar Zongo (16° S, 68° W) at 500 hPa (around 5500 m a.s.l.) and an empirical relationship (Favier and others, 2008) given by $E = \alpha(q - q_s)\nu$ (mm w.e. month⁻¹) where $q$, $q_s$ and $\nu$ are the re-analyzed NCEP-NCAR data for specific humidity (g kg⁻¹), saturated specific humidity (g kg⁻¹) and horizontal wind speed (m s⁻¹) respectively, and $\alpha$ is a coefficient for homogeneity of the empirical relationship (Favier and others, 2008). This coefficient has been calibrated with monthly sublimation data measured in the field from automatic weather stations by Wagnon and others (2001) (September 1996–August 1998) and Sicart (2002) (September 1999–August 2000) except for a data gap in March 2000). The correlation $r^2$ over the 35 month period is 0.71, and yields a value for $\alpha$ of 0.694. The catchment and glacier areas were properly determined from photogrammetric data.

The third and independent method used to assess glacier mass balance comes from geodetic data. It can be used to calculate the volume change and, consequently, the cumu-
<figure><img src="image_4.png" /><figcaption>**Fig. 3.** Thickness variations (m) of Glaciar Zongo over the periods 1956–63, 1963–75, 1975–83, 1983–97, 1997–2006 and 1956–2006. The color scale is the same in each case. The white areas corresponding to shadowed or low-contrast areas were not measured. The surface areas of Glaciar Zongo are 2.29, 2.31, 2.26, 2.23, 2.12 and 1.96 km², corresponding to 1956, 1963, 1975, 1983, 1997 and 2006 photogrammetric flights respectively.</figcaption></figure>

lative mass balance of the glacier using digital photogram-
metric data. Six photogrammetric flights have been carried
out over Glaciar Zongo on 20 May 1956, 21 June 1963, 23
July 1975, 20 June 1983, 20 April 1997 and 13 July 2006 by
the National Service of Aerophotogrammetry, Bolivia (SNA).
The resulting photographs were scanned with a resolution of
14 µm using a photogrammetric scanner (Vexcel UltraScan
5000). In order to orient the images, a stereo preparation
was performed in 2006 using the geodetic differential global
positioning system (GPS) on ground-control points. Since the
photographs had been taken earlier, obvious features were
used as control points. Moreover, as the relative orientation
of the images between two dates is more important than the
absolute ground orientation, at least 25 tie points were
added per stereo pair to improve consistency (Kruck, 2002).
The photogrammetric restitution was performed using a
Leica digital system (Orima and LPS®) following an irregular
grid: every 20 m on flat surfaces of the glacier and as low as
every 5 m in bumpy areas. Approximately 10 000 points
were measured for each aerial campaign, over the whole
area of the glacier (average area 2.2 km²). The average
measurement accuracy is 7.5 µm for photograph coordinates
and 0.64 m for ground-control point coordinates. Some
areas with low contrast or shadows were not measured
(shown in white in Fig. 3)

The digital elevation models (DEMs) were constructed
using Surfer software with a grid of 10 × 10 m and a
minimum curvature algorithm. The algorithm employed for
DEM construction has a weak impact on the whole glacier

volume (Thibert and others, 2008). The subtraction of the
DEMs yields variations of surface elevation that must be
converted to water equivalent using the ice density, i.e.
900 kg m⁻³ (see section 5).

4. RESULTS

The geodetic method gives the long-term cumulative mass
balance. Based on photogrammetric measurements, digital
maps were plotted for 1956, 1963, 1975, 1983, 1997 and
2006, and thickness variations were calculated from these
data over the corresponding periods (Fig. 3). Some global
features can be observed. Over the first period (1956–63),
changes in volume are small, while the second period
(1963–75) displays a positive volume variation in the highest
part of the glacier. The period 1975–2006 shows a strong
change of volume in the ablation zone. Over the whole
period 1956–2006, a gradual loss of volume is observed in
the ablation zone, with no significant changes in the upper
part of the glacier. Finally, the cumulative mass balance of
Glaciar Zongo calculated over the last 50 years using
photogrammetric and hydrological methods is displayed in
Figure 4 and Table 1. The figure reveals a striking agreement
between independent data obtained by the hydrological and
photogrammetric methods. The cumulative glaciological
mass balances are also plotted in Figure 4 and shown in
Table 1. However, as explained in the next section, these
glaciological results have been adjusted according to the
photogrammetric data.
<figure><img src="image_5.png" /><figcaption>**Fig. 4.** Zongo cumulative mass balance (mm w.e.) calculated by the geodetic method (black triangles), the hydrological method (grey curve) and the glaciological method (black curve). The hydrological mass balances were obtained from annual data between 1974 and 2006. The glaciological mass balances were calculated using yearly field measurements and adjusted according to photogrammetric data. The origin of the hydrological mass balance is the geodetic value in 1975. The hydrological mass balance envelops correspondents to one standard deviation (thin grey curves).</figcaption></figure>

Glaciar Zongo lost 16.32 m w.e. in mass between 1956 and 2006 and 15% of its surface area between 1956 and 2006. It seems that between 1956 and 1975 the glacier was near steady-state conditions. Conversely, the periods 1975–83 and 1997–2006 show strong negative mass-balance rates of –0.91 and –0.84 m w.e. a⁻¹ respectively.

## 5. DISCUSSION

The mass balance of Glaciar Zongo has been assessed using three independent methods. Uncertainties related to each method are analyzed below.

The comparison between glaciological and hydrological mass balance reveals a large discrepancy over the period of the common data, i.e. 1991–2006. The cumulative mass balances resulting from the two methods diverge by 7.4 m w.e. over the 15 years of observations.

To try to explain this discrepancy, glaciological surface mass-balance data from stake and pit measurements were analyzed. Three areas of 4 ha were selected in the ablation area at 5050, 5100 and 5150 m a.s.l., and three others in the accumulation area at 5550, 5650 and 5750 m a.s.l. These areas correspond to surface mass-balance field measurements covering the entire period. In each area, centered mass balances were calculated by subtracting the average mass balance from the annual values (Fig. 5a). The figure shows striking consistency for the data obtained from the ablation area, despite the high variability. The relationships between surface mass balance in the ablation area are strong, with correlation coefficients between r² = 0.92 and r² = 0.86. This confirms the validity of the surface mass-balance data in this area. Figure 5a also shows the large variability difference between the accumulation and ablation areas. The standard deviations are high for the elevation ranges of the ablation area (σ = 1.6, 1.5 and 1.4 m w.e. a⁻¹ for elevations 5050, 5100 and 5150 m respectively, over the period 1991–2006) and are low for the elevation ranges of the accumulation area (σ = 0.5, 0.4 and 0.4 m w.e. a⁻¹ for elevations 5550, 5650 and 5750 m respectively, over the period 1991–2006). Figure 5a suggests that specific mass-balance change is mainly controlled by the ablation area. Over the whole glacier area, the standard deviation ranges roughly from 1.6 m w.e. a⁻¹ in the ablation area to 0.4 m w.e. a⁻¹ in the accumulation area. Assuming constant surface areas, the

variance of the specific mass balance is given by σ² = Σ (σᵢ² Sᵢ² / S²), where Sᵢ and σᵢ are the surface area and surface mass-balance standard deviation of elevation range i, and S is the surface area of the glacier (i = 1, ...). Over the

<table><caption>**Table 1.** Zongo cumulative mass balance (m w.e.) calculated by the hydrological, glaciological and geodetic methods</caption><thead><tr><th>Period</th><th>Cumulative hydrological mass balance</th><th>Adjusted cumulative glaciological mass balance</th><th>Period</th><th>Cumulative geodetic mass balance</th></tr></thead><tbody><tr><td>1973/74</td><td>3.60</td><td></td><td></td><td></td></tr><tr><td>1974/75</td><td>4.35</td><td></td><td></td><td></td></tr><tr><td>1975/76</td><td>4.35</td><td></td><td></td><td></td></tr><tr><td>1976/77</td><td>3.24</td><td></td><td>1956–63</td><td>0.90</td></tr><tr><td>1977/78</td><td>2.72</td><td></td><td></td><td></td></tr><tr><td>1978/79</td><td>2.33</td><td></td><td></td><td></td></tr><tr><td>1979/80</td><td>1.25</td><td></td><td></td><td></td></tr><tr><td>1980/81</td><td>0.16</td><td></td><td></td><td></td></tr><tr><td>1981/82</td><td>-0.77</td><td></td><td></td><td></td></tr><tr><td>1982/83</td><td>-3.06</td><td></td><td>1963–75</td><td>3.86</td></tr><tr><td>1983/84</td><td>-3.78</td><td></td><td></td><td></td></tr><tr><td>1984/85</td><td>-3.86</td><td></td><td></td><td></td></tr><tr><td>1985/86</td><td>-1.66</td><td></td><td></td><td></td></tr><tr><td>1986/87</td><td>-2.14</td><td></td><td></td><td></td></tr><tr><td>1987/88</td><td>-3.02</td><td></td><td></td><td></td></tr><tr><td>1988/89</td><td>-4.27</td><td></td><td></td><td></td></tr><tr><td>1989/90</td><td>-5.34</td><td></td><td>1975–83</td><td>-3.45</td></tr><tr><td>1990/91</td><td>-5.78</td><td></td><td></td><td></td></tr><tr><td>1991/92</td><td>-7.13</td><td></td><td></td><td></td></tr><tr><td>1992/93</td><td>-6.99</td><td></td><td></td><td></td></tr><tr><td>1993/94</td><td>-7.07</td><td></td><td></td><td></td></tr><tr><td>1994/95</td><td>-8.05</td><td></td><td></td><td></td></tr><tr><td>1995/96</td><td>-9.04</td><td></td><td>1983–97</td><td>-8.81</td></tr><tr><td>1996/97</td><td>-9.03</td><td></td><td></td><td></td></tr><tr><td>1997/98</td><td>-11.58</td><td>-11.28</td><td></td><td></td></tr><tr><td>1998/99</td><td>-13.71</td><td>-12.18</td><td></td><td></td></tr><tr><td>1999/00</td><td>-14.16</td><td>-12.52</td><td></td><td></td></tr><tr><td>2000/01</td><td>-14.24</td><td>-12.31</td><td></td><td></td></tr><tr><td>2001/02</td><td>-14.97</td><td>-12.80</td><td></td><td></td></tr><tr><td>2002/03</td><td>-15.06</td><td>-13.47</td><td>1997–2006</td><td>-16.32</td></tr><tr><td>2003/04</td><td>-15.55</td><td>-14.30</td><td></td><td></td></tr><tr><td>2004/05</td><td>-16.69</td><td>-16.00</td><td></td><td></td></tr><tr><td>2005/06</td><td>-16.87</td><td>-16.38</td><td></td><td></td></tr></tbody></table>
<figure><img src="image_6.png" /><figcaption>**Fig. 5.** (a) Centered annual surface mass balance (m w.e. a⁻¹) for different elevation ranges; and (b) centered specific mass balance (m w.e. a⁻¹) obtained from glaciological data and from hydrological data.</figcaption></figure>

period 1991–2006, the surface mass-balance standard deviation is estimated to be −0.75 m w.e. a⁻¹, although it cannot be determined properly over several elevation ranges. The values ($\sigma_t^2 S_t^2$)/($\sigma^2 S^2$) can be used to assess the relative contribution of each elevation range $i$ to the specific mass-balance change for 1 year. The very large contribution of low-elevation ranges to the specific mass balance is apparent: over the period 1991–2006, more than 80% of the variance comes from elevation ranges below 5300 m, representing roughly one-third of the total surface area. The contribution of the tongue changes in time with the area changes. For instance, the mass-balance changes between 4900 and 5200 m accounted for about 73% of the specific mass-balance change in 1991, although the surface area was 25% of the glacier. This contribution fell to 24% in 2006 following the area reduction. The centered specific mass balances are plotted in Figure 5b, using glaciological and hydrological data. Despite large unexplained discrepancies in 1999, fluctuations of the two sets of mass-balance data are consistent, with a correlation coefficient of $r^2 = 0.80$.

However, the difference between the average glaciological and hydrological mass balance remains. Provided that surface mass-balance measurements are reliable and that surface area changes have been properly taken into account, the discrepancy should come from the sampling network. It can be concluded that the field measurement sites are not representative of the entire area of the glacier and that the linear interpolation between surface mass-balance measurements from different elevation ranges (5250–5550 m) is not valid for calculating the specific mass balance. Consequently, to determine the specific mass balance, the glaciological data have been adjusted according to the photogrammetric data between 1997 and 2006. For this purpose, the average volumetric mass balance over this period has been added to the centered specific glaciological mass balance calculated over the same period (Fig. 4).

Error analysis of photogrammetric results was performed using the method of Thibert and others (2008) and is shown in Table 2. Additional uncertainties result from snow/ice density uncertainties and unmeasured areas resulting from lack of stereoscopic measurements. The overall error in the cumulative mass balance obtained from photogrammetric data is estimated at 2.0 m w.e. over the entire period 1956–2006. The uncertainty of the hydrological method is difficult to assess accurately. As shown by Sicart and others (2007), the runoff discharge error is small, i.e. 5% in the range 0–600 L s⁻¹ (Ribstein and others, 1995). Additionally, provided that the portion of basin area free of ice is reduced, the influence of the runoff coefficient is small. The maximum changes in mass balance due to changes in the runoff coefficient from 0.5 to 1 are −0.20 and +0.15 m w.e. (Sicart and others, 2007). The average sublimation rate calculated over the last 50 years over Glaciar Zongo is 171 mm w.e. a⁻¹. The uncertainty of this parameter is difficult to assess and it can be assumed that it is less than 0.1 m w.e. a⁻¹. Therefore, the uncertainty of the hydrological mass balance depends mainly on the accumulation uncertainties. As mentioned by Sicart and others (2007), no clear relationship can be found between altitude and precipitation. The field glaciological measurements show that the accumulation rate is low both in the ablation zone and in the accumulation zone (about 1 m w.e. a⁻¹). It is assumed that the uncertainty in the accumulation rate calculated from precipitation should not exceed 0.2 m w.e. a⁻¹. As a result, the uncertainty of annual hydrological mass balance is 0.3 m w.e. a⁻¹ maximum and the overall error in the cumulative mass balance from hydrological data is estimated at 1.68 m w.e. over the period 1975–2006 (Fig. 4). The very good agreement with photogrammetric data confirms that glacier mass balance obtained from hydrological data is reliable.

## 6. CONCLUSIONS

These glacier cumulative mass-balance results over the last 50 years are unique for this region. Hydrological data have made it possible to reconstruct annual glacier mass balance since 1974.

Only two other glacier mass-balance series have been investigated in the Cordillera Real over a similar period: those for the Charquini glaciers located 5 km from Glaciar Zongo, and Glaciar Chacaltaya 9 km away. The cumulative mass balance of the Charquini glaciers between 1956 and 1997 was assessed by Rabatel (2005) using photogrammetric data. Glaciar Charquini Sur (0.36 km²; the glacier flows approximately from 5280 to 4985 m a.s.l.), with similar ex-position, shows similar fluctuations over the whole period, with a strong negative mass balance of −0.79 m w.e. a⁻¹ between 1975 and 1983.

The cumulative mass balance of Glaciar Chacaltaya (0.007 km²; the glacier flows approximately from 5325 to 5225 m a.s.l.) has been roughly reconstructed between 1963 and 1998 using photogrammetric and stake-measurements data (Ramirez and others, 2001). The mean mass balance is very negative, −0.75 m w.e. a⁻¹ over the whole period. Over the period 1983–98 this series shows a stronger negative mass balance of −0.96 m w.e. a⁻¹. However, this glacier will disappear in the near future, and its specific mass balance was mainly controlled by the high surface-area change over recent decades.
<table><caption>**Table 2.** Error analysis in the volumetric mass-balance measurements (see Thibert and others, 2008, for details)</caption><thead><tr><th rowspan="2">Random error in DEM</th><th colspan="6">Year</th></tr><tr><th>1956</th><th>1963</th><th>1975</th><th>1983</th><th>1997</th><th>2006</th></tr></thead><tbody><tr><td>Surface roughness error for a single point, σ<sub>rgh</sub> (m)</td><td>0.35</td><td>0.35</td><td>0.35</td><td>0.35</td><td>0.35</td><td>0.35</td></tr><tr><td>Planimetric stereoscopic error (single point), σ<sub>stereo.x</sub> (m)</td><td>1.58</td><td>1.21</td><td>2.69</td><td>2.04</td><td>2.01</td><td>1.38</td></tr><tr><td>Altimetric stereoscopic error (single point), σ<sub>stereo.z</sub> (m)</td><td>1.65</td><td>1.21</td><td>3.01</td><td>2.89</td><td>2.24</td><td>3.5</td></tr><tr><td>Mean altimetric error for the DEM, σ'<sub>z</sub> (m)</td><td>0.02</td><td>0.02</td><td>0.03</td><td>0.03</td><td>0.02</td><td>0.03</td></tr><tr><td rowspan="3">Orientation errors (ground-control points), σ<sub>x</sub>, σ<sub>y</sub>, σ<sub>z</sub> (m)</td><td>0.55</td><td>0.75</td><td>0.92</td><td>0.75</td><td>0.72</td><td>0.45</td></tr><tr><td>0.40</td><td>0.67</td><td>0.56</td><td>0.97</td><td>0.71</td><td>0.48</td></tr><tr><td>0.78</td><td>0.39</td><td>0.69</td><td>0.65</td><td>0.67</td><td>0.43</td></tr><tr><td>Altimetric orientation error, σ<sub>a.o</sub> (m)</td><td>0.8</td><td>0.48</td><td>0.75</td><td>0.73</td><td>0.72</td><td>0.47</td></tr></tbody></table>

As shown in other studies of tropical glaciers (Kaser, 1999), the vertical budget gradient is high (~2 m w.e. (100 m)⁻¹). Furthermore, our study reveals strong temporal variability of the surface mass balance in the ablation area, more than twice the mean variability.

Mass-balance fluctuations derived from hydrological, glaciological and photogrammetric data agree well. However, the glaciological method reveals a strong bias in specific mass balance, as many areas cannot be sampled because of difficult access. The glaciological mass balance has therefore been corrected using photogrammetric results. The photogrammetric method is therefore needed to provide a reliable average specific mass balance over the whole period. Generally, the hydrological method depends strongly on the correction of the measured precipitation in the catchment area. However, for Glaciar Zongo it seems that accumulation changes are low and do not significantly influence hydrological mass-balance changes. Comparison between annual surface mass balances observed from stake measurements reveals very similar temporal changes and proves that surface mass-balance changes are strongly related to meteorological conditions. On the whole, the glacier was close to steady-state conditions between 1956 and 1963, between 1963 and 1975 it gained mass, and after 1975 the specific mass balance is negative. The strong temporal variability of the surface mass balance in the ablation area leads to a large influence of surface area reduction. The relative contribution of the tongue (between 5200 m and the terminus) to the specific mass balance has decreased from 73% in 1991 to 24% in 2006.

As a result, the surface area of this glacier converges to steady state faster than glaciers in other regions of the world.

## ACKNOWLEDGEMENTS

The Great Ice glaciological program is supported by Institut de Recherche pour le Développement (IRD) and by the French Observation de Recherche en Environnement (ORE) program GLACIOCLIM (http://www-lgge.ujf-grenoble.fr/ServiceObs/index.htm). We thank everyone who took part in the extensive field measurements on Glaciar Zongo. We are grateful for the assistance of Instituto de Hidráulica e Hidrología (IHH) and Compañia Boliviana de Energia Eléctrica (COBEE) in La Paz. We thank the US National Oceanic and Atmospheric Administration (NOAA)/Office of Oceanic and Atmospheric Research (OAR)/Earth and Space Research Laboratory (ESRL) Physical Sciences Division (PSD), Boulder, Colorado, USA, for providing the NCEP re-

analysis data from their website at http://www.cdc.noaa.gov/. We are grateful to the anonymous reviewers whose comments improved the quality of the manuscript.

## REFERENCES

Bauder, A., M. Funk and M. Huss. 2007. Ice-volume changes of selected glaciers in the Swiss Alps since the end of the 19th century. *Ann. Glaciol.*, **46**, 145–149.

Coudrain, A., B. Francou and Z.W. Kundzewicz. 2005. Glacier shrinkage in the Andes and consequences for water resources – editorial. *Hydrol. Sci. J.*, **50**(6), 925–932.

Dyurgerov, M. 2002. *Glacier mass balance and regime: data of measurements and analysis*. Boulder, CO, University of Colorado. Institute of Arctic and Alpine Research. (INSTAAR Occasional Paper 55.)

Favier, V., P. Wagnon and P. Ribstein. 2004. Glaciers of the outer and inner tropics: a different behaviour but common response to climatic forcing. *Geophys. Res. Lett.*, **31**(16), L16403. (10.1029/2004GL020654.)

Favier, V. and 8 others. 2008. Evidence of groundwater flow on Antizana ice-covered volcano, Ecuador. *Hydrol. Sci. J.*, **53**(1), 278–291.

Francou, B., P. Ribstein, R. Saravia and E. Tiriau. 1995. Monthly balance and water discharge of an inter-tropical glacier: Zongo Glacier, Cordillera Real, Bolivia, 16° S. *J. Glaciol.*, **41**(137), 61–67.

Hastenrath, S. 1984. *The glaciers of equatorial East Africa*. Dordrecht, etc., D. Reidel Publishing Co.

Hastenrath, S. 1995. Glacier recession on Mount Kenya in the context of the global tropics. *Bull. Inst. Fr. Étud. Andin. [Lima]*, **24**(3), 633–638.

Kalnay, E. and 21 others. 1996. The NCEP/NCAR 40-year reanalysis project. *Bull. Am. Meteorol. Soc.*, **77**(3), 437–471.

Kaser, G. 1999. A review of the modern fluctuations of tropical glaciers. *Global Planet. Change*, **22**(1–4), 93–103.

Kaser, G. 2001. *Glacier-climate interaction at low latitudes*. *J. Glaciol.*, **47**(157), 195–204.

Kaser, G., S. Hastenrath and A. Ames. 1996. Mass balance profiles on tropical glaciers. *Z. Gletscherkd. Glazialgeol.*, **32**, 75–81.

Kruck, E. 2002. *BINGO–F: bundle adjustment for engineering applications. Version 4.0*. Aalen, Gesellschaft für Industrie-photogrammetrie.

Paterson, W.S.B. 1994. *The physics of glaciers. Third edition*. Oxford, etc., Elsevier.

Rabatel, A. 2005. Chronologie et interprétation paléoclimatique des fluctuations des glaciers dans les Andes de Bolivie (16° S) depuis le maximum du Petit Age Glaciaire (17ème siècle). (PhD thesis, Université Joseph Fourier.)

Rabatel, A., A. Machaca, B. Francou and V. Jomelli. 2006. Glacier recession on Cerro Charquini (16° S), Bolivia, since the
maximum of the Little Ice Age (17th century). *J. Glaciol.*, **52**(176), 110–118.

Ramírez, E. and 8 others. 2001. Small glaciers disappearing in the tropical Andes: a case-study in Bolivia: Glaciar Chacaltaya (16° S). *J. Glaciol.*, **47**(157), 187–194.

Ribstein, P., E. Tiriau, B. Francou and R. Saravia. 1995. Tropical climate and glacier hydrology: a case study in Bolivia. *J. Hydrol.*, **165**(1–4), 221–234.

Ribstein, P., B. Pouyaud, J.E. Sicart, P. Wagnon, E. Ramírez and B. Francou. 1999. Variabilité climatique et fonctionnement hydrologique d'un glacier tropical. *In Comité National Français de Géodésie et Géophysique, Rapport Quadriennal, 1995–1998, Section 6*. Paris, Comité National Français de Géodésie et Géophysique, 279–287.

Sicart, J.E. 2002. Contribution à l'étude des lux d'energie, du bilan de masse et du debit du fonte d'un glacier tropical: le Zongo, Bolivie. (PhD thesis, Université de Paris.)

Sicart, J.E., P. Wagnon and P. Ribstein. 2005. Atmospheric controls of the heat balance of Zongo Glacier (16°S, Bolivia). *J. Geophys. Res.*, **110**(D12), D12106. (10.1029/2004JD005732.)

Sicart, J.E., P. Ribstein, B. Francou, B. Pouyaud and T. Condom. 2007. Glacier mass balance of tropical Glaciar Zongo, Bolivia, comparing hydrological and glaciological methods. *Global Planet. Change*, **59**(1–4), 27–36.

Solomon, S. and 7 others, eds. 2007. *Climate change 2007: the physical science basis. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change*. Cambridge, etc., Cambridge University Press.

Thibert, E., R. Blanc, C. Vincent and N. Eckert. 2008. Glaciological and volumetric mass-balance measurements: error analysis over 51 years for Glacier de Sarennes, French Alps. *J. Glaciol.*, **54**(186), 522–532.

Vergara, W. and 6 others. 2007. Economic impacts of rapid glacier retreat in the Andes. *Eos*, **88**(25), 261, 264.

Vincent, C. 2002. Influence of climate change over the 20th century on four French glacier mass balances. *J. Geophys. Res.*, **107**(D19), 4375. (10.1029/2001JD000832.)

Wagnon, P., P. Ribstein, B. Francou and B. Pouyaud. 1999a. Annual cycle of energy balance of Zongo Glacier, Cordillera Real, Bolivia. *J. Geophys. Res.*, **104**(D4), 3907–3924.

Wagnon, P., P. Ribstein, G. Kaser and P. Berton. 1999b. Energy balance and runoff seasonality of a Bolivian glacier. *Global Planet. Change*, **22**(1–4), 49–58.

Wagnon, P., P. Ribstein, B. Francou and J.E. Sicart. 2001. Anomalous heat and mass budget of Glaciar Zongo, Bolivia, during the 1997/98 El Niño year. *J. Glaciol.*, **47**(156) 21–28.