# Glacier recession on Cerro Charquini (16° S), Bolivia, since the maximum of the Little Ice Age (17th century)

Antoine Rabatel, Abraham Machaca, Bernard Francou, Vincent Jomelli

► To cite this version:

Antoine Rabatel, Abraham Machaca, Bernard Francou, Vincent Jomelli. Glacier recession on Cerro Charquini (16° S), Bolivia, since the maximum of the Little Ice Age (17th century). Journal of Glaciology, 2006, 52 (176), pp.110 à 118. 10.3189/172756506781828917. insu-00375748

HAL Id: insu-00375748

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

Submitted on 23 Nov 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
# Glacier recession on Cerro Charquini (16° S), Bolivia, since the maximum of the Little Ice Age (17th century)

Antoine RABATEL,¹ Abraham MACHACA,² Bernard FRANCOU,³ Vincent JOMELLI⁴

¹Laboratoire de Glaciologie et Géophysique de l'Environnement, 54 rue Molière, BP 96, 38402 Saint-Martin-d'Hères Cedex, France
E-mail: rabatel@lggeobs.ujf-grenoble.fr

²Instituto Geografico Militar, Av. Saavedra, Estado Mayor de Miraflores, La Paz, Bolivia

³CP 9214, La Paz, Bolivia

⁴Maison des Sciences de l'Eau, BP 64501, 34394 Montpellier, France

**ABSTRACT.** Cerro Charquini, Bolivia (Cordillera Real, 5392 m.a.s.l.) was selected as a site to reconstruct glacier recession since the maximum of the Little Ice Age (LIA) in the central Andes. Five glaciers, located on differently exposed slopes, present comprehensive and well-preserved morainic systems attributed to former centuries. The moraines were dated by lichenometry and show a consistent organization on the different slopes. The past geometry of the glaciers was reconstructed using ground topography and aerophotogrammetry. Lichenometric dating shows that the LIA maximum occurred in the second half of the 17th century, after which the glaciers have receded nearly continuously. Over the last decades of the 20th century (1983–97), recession rates increased by a factor of four. On the northern and western slopes, glaciers receded more than on the southern and eastern slopes (by 78% and 65% of their LIA maximum area, respectively). The mean equilibrium-line altitude (ELA) rose by about 160 m between the LIA maximum and 1997. Recession rates were analysed in terms of climatic signal, suggesting that glacier recession since the LIA maximum was mainly due to a change in precipitation and that the 19th century may have been drier. For the 20th century, a temperature rise of about 0.6°C appears to be the main cause of glacier recession. Recent climatic conditions from 1983 to 1997 correspond to a mass deficit of about 1.36 m w.e. a⁻¹. If such conditions persist, the small glaciers below 5300 m a.s.l. in the Cordillera Real should disappear completely in the near future.

## 1. INTRODUCTION

Much evidence indicates that tropical glaciers were far more extensive in past centuries. This is especially true during the Little Ice Age (LIA), which is well documented for some mountain ranges in the Northern Hemisphere such as the European Alps and the Canadian Rockies (Le Roy Ladurie, 1967; Grove, 1988; Luckman, 2000), and also identified in the tropical Andes (Broggi, 1945; Hastenrath, 1981; Clapperton, 1983; Ames and Francou, 1995; Francou, 2004). In the Cordillera Real, many moraines attributed to the LIA exist below the present glacier snouts, and the maximum extent has been considered the most prominent phase of glacier advance of the entire Holocene (Gouze and others, 1986). However, few studies have yet been able to accurately date glacier fluctuations over this period (Solomina and others, in press). The large uncertainties involved in dating these glacier fluctuations are partly due to the dating method using ¹⁴C (Gouze and others, 1986; Seltzer, 1992), an isotope with atmospheric production rates that have varied over recent centuries, precluding dating within the last 300 years (Reimer and others, 2004). Information obtained from ice cores retrieved in the Andes has made it possible to roughly estimate the LIA's limits between the 15th and the late 19th centuries (Thompson and others, 1986). However more accurate dating of glacier evolution within this period requires other proxies.

Moraines are often used to date fluctuations and to reconstruct glacier areas and volumes with good accuracy. In the tropics, moraines are well preserved because reworking by frost penetration in the ground and snowmelt

is extremely limited. Moreover, time lags between glacier mass-balance variations and the dynamic response of the snouts are generally short, i.e. one to a few years (Francou and others, 2004), particularly for the small glaciers (<0.5 km²) that represent about 80% of the glacier area of the Cordillera Real, Bolivia (Jordan, 1991). Thus the moraines deposited by these small glaciers are believed to represent the form of the glaciers in an equilibrium state.

The aim of this paper is to describe the glacier evolution since the LIA maximum for this tropical Cordillera using moraines dated by lichenometry and to infer climate changes from these fluctuations.

## 2. STUDY AREA

Cerro Charquini (16°17'S, 68°09'W) is located 20 km northeast of La Paz, midway between Glaciar Zongo and Glaciar Chacaltaya where a permanent observation programme involving mass-balance, hydrological and energy-balance measurements has been run for more than a decade (Francou and others, 1995, 2003; Wagnon and others, 1999).

Despite the relatively low elevation of the summit (5392 m.a.s.l.), glaciers with similar morphology are present on the different slopes, allowing aspect effects to be analysed without interference from other factors. The Cerro Charquini glaciers are small (<0.5 km² in 1997) and, together with their moraines, are all located within the same 4430–5300 m a.s.l. elevation range (Fig. 1). Moraines extend 1.5 km and 200–450 m in elevation beyond the
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Cerro Charquini glaciers (Cordillera Real, Bolivia, 16° S). Although the ten moraines can be found on each glacier foreland, only the five biggest (M1, M3, M6, M8 and M9) are shown. Glacier limits in 1940, 1983 and 1997 are from stereophotogrammetric restitutions.</figcaption></figure>

present snouts. As they have not yet been covered by vegetation, these moraines may be considered, a priori, as dating from the LIA. Moraines and bedrock are granodiorites and quartzites, rock types that favour growth of *Rhizocarpon* sp., a lichen commonly used to date moraines worldwide (Innes, 1985).

In the Cordillera Real, annual precipitation amounts to 800–1000 mm at above 4000 m a.s.l., which falls generally as snow on the glaciers. Precipitation comes with humid air masses from the northeast, according to a mechanism called the Amazonian monsoon, and presents a very contrasted seasonal pattern, with 80% falling from October to April (austral summer) and 65% from December to February. This marked seasonality of the climate plays an important role in the mass-balance fluctuations. During the austral summer, accounting for 90% of the annual mass-balance variability (Francou and others, 2003), the maximum net accumulation occurs in the upper part of the glaciers while ablation rates are highest in the lower part. Ablation is associated with a strong total radiation influx to the glacier surface (the sun is vertical and cloudiness is maximum from November to February), high humidity and low winds (which limit

sublimation) and a relatively high air temperature. Precipitation controls the energy balance at the glacier surface through the albedo (Wagnon and others, 1999). Every variation in the precipitation regime results in a change to the glacier surface albedo and consequently the ablation intensity (Sicart, 2002).

On a regional scale, the El Niño–Southern Oscillation (ENSO) controls the majority of the interannual variability of mass balance (Francou and others, 2003). Warm phases are generally accompanied by a precipitation deficit of 10–30% (i.e. less accumulation leading to a lower albedo), a decrease in cloudiness (i.e. more shortwave radiation reaching the glacier surface) and an increase in the sensible heat flux (air temperature anomaly in the +1° to +3°C range) (Wagnon and others, 2001). Consequently mass balance is very negative during El Niño years. An opposite situation occurs during the cold phases (La Niña), when mass balance shows near-equilibrium/positive values due to enhanced precipitation, colder air temperature and high cloudiness, which tend to decrease ablation strength and lower the ELA by 100–300 m compared with the warm phases. The climate during the LIA is believed to have been more extreme than at
<figure><img src="image_3.png" /><figcaption>Fig. 2. Glaciar Charquini Sur. Surface area reconstructed on the basis of the moraines until M10 and by aerial photography from 1940.</figcaption></figure>

present, and a direct effect of the ENSO variability during the last centuries on glacier fluctuations and moraine deposits has to be considered.

# 3. METHODS

## 3.1. Dating moraines by lichenometry

The dating of Charquini moraines by lichenometry with statistical analysis based on the extreme-value theory has been discussed elsewhere (Rabatel and others, 2005; Cooley and others, in press; Naveau and others, in press) and only the basic principles will be presented here.

The basic premise of lichenometry is that the diameter of the lichens growing on a surface is proportional to the length of time that the surface has been exposed to colonization and growth in a specific environmental context. Thus the largest lichens can be used to date the deposit, the stabilization or the exposure of surfaces.

All measurements were carried out with the *Rhizocarpon geographicum* s.l. lichen, which presents a circular growth pattern and has a lifetime longer than the period to be dated (Innes, 1985). Our dataset of lichenometric measurements is composed of 10 previously dated and 48 further surfaces. The surfaces of known age are presented in detail by Rabatel (2005). All were found in the eastern Cordillera (up to 4500 m a.s.l.), and so are expected to be representative of lichen growth rate on the glacier forelands. Eight date from the 20th century (between 1910 and 1965), one from the mid-18th century and one from the mid-17th century. Only the oldest one is based on a ¹⁴C dating. All the others are dated using documentary sources. Undated surfaces correspond to the main moraines of each of the glaciers to be dated.

In the current study, a new method based on the extreme-value theory called generalized extreme values (GEV; Katz and others, 2002) was used to analyse lichenometric data. Maximum lichen diameters and average surface ages are assumed to have two different distributions, a GEV and a normal distribution, respectively. This statement is based on

probability theory for maxima (extreme-value theory) and for averages (central limit theorem), respectively. These two distributions are linked through the relationship $\alpha = f(\mu)$, where the function $f$ represents the temporal variation of $\mu$ (for more details see the papers mentioned above).

## 3.2. Reconstructing glaciers from moraines and aerial photographs

Glacier fluctuations during this period were reconstructed on the basis of moraines. All the moraines of the LIA and the beginning of the 20th century were dated by lichenometry. For the second part of the 20th century, five aerial photograph pairs dated 1956, 1963, 1974, 1983 and 1997 (Instituto Geográfico Militar, Bolivia) were processed using classical (non-digital) photogrammetric restitution techniques. This information allowed us to reconstruct glacier areas and to calculate volume changes. Aerial pictures from 1940 could not provide a three-dimensional restitution because the flight parameters were not available. Nevertheless, it was possible to calculate the glacier area of 1940 by plotting the limits on the 1/10 000 map based on other aerial photographs. Restitution was also used to obtain a complete 1/10 000 mapping of Charquini on which we represent the ten main LIA morainic stages for each of the five glaciers (see Fig. 2 which shows an example, Glaciar Charquini Sur). In order to validate the results of the restitution, a complete ground topographic survey of each moraine was performed on Glaciar Charquini Sur. Results of both methods are consistent and allow us to be confident that the mapping of glaciers and moraines for the whole Charquini area is accurate.

## 3.3. ELA reconstruction

The equilibrium-line altitude (ELA) was reconstructed for each morainic stage and each pair of photographs using the accumulation-area ratio (AAR) method and a ratio of 0.65 determined for Glaciar Zongo when mass balance is in equilibrium (1991–2004 series, see Soruco and others, 2005). This 2:3 ratio was also observed on other south-facing glaciers of the sector (Rabatel and Mendoza, 2004) and thus it can be considered as a realistic approximation for AAR₀ on a regional scale. Nevertheless this ratio is valid only on glaciers that are in equilibrium, and may change significantly according to slope aspect on the scale of the Charquini area.

## 3.4. Ice-volume and mass-balance reconstruction

For the ten LIA stages, glacier volumes were reconstructed by contour line interpolation based on the moraine height and the bedrock morphology. Since depth to the bedrock below the current glaciers was unknown, we could only calculate volume fluctuations between the different stages. For the 20th century, the volume variations computed from the photogrammetric restitution of each pair of photographs were considered to be accurate. After dating the moraines and aerial photographs, it was possible to calculate the mean annual mass balance for each period, given by the following equation:

$$ \overline{B_{(T_i, T_{i+1})}} = \rho \frac{(V_{i+1} - V_i)}{(T_{i+1} - T_i) \times \frac{(S_{i+1} + S_i)}{2}} $$

where $\overline{B_{(T_i, T_{i+1})}}$ is the mean annual mass balance for the period $T_i$ to $T_{i+1}$, $T_i$ is the glacier stage date (moraine or aerial
<table><caption>**Table 1.** Moraine dating by lichenometry with the associated error (years), surface area (km²) and ELA (m a.s.l.) for the LIA stages and six 20th-century aerial photographs of the Charquini glaciers. On Glaciar Charquini Sureste blocks were too scarce on the M2 and M5 moraines to allow lichenometric measurements</caption><thead><tr><th rowspan="2">Moraine</th><th colspan="3">Charquini Sur</th><th colspan="3">Charquini Sureste</th><th colspan="3">Charquini Noreste</th><th colspan="3">Charquini Norte</th><th colspan="3">Charquini Oeste</th></tr><tr><th>Date</th><th>Surface area</th><th>ELA</th><th>Date</th><th>Surface area</th><th>ELA</th><th>Date</th><th>Surface area</th><th>ELA</th><th>Date</th><th>Surface area</th><th>ELA</th><th>Date</th><th>Surface area</th><th>ELA</th></tr></thead><tbody><tr><td>LIA max</td><td>1686 ± 14</td><td>1.22</td><td>4930</td><td>1664 ± 14</td><td>1.41</td><td>4815</td><td>1662 ± 14</td><td>1.04</td><td>4870</td><td>1663 ± 14</td><td>1.06</td><td>4990</td><td>1663 ± 14</td><td>0.83</td><td>4910</td></tr><tr><td>M2</td><td>1703 ± 12</td><td>1.16</td><td></td><td>1700 ± 12</td><td>1.02</td><td></td><td>1700 ± 12</td><td>1.02</td><td></td><td>1706 ± 12</td><td>1.04</td><td></td><td>1700 ± 12</td><td>0.78</td><td></td></tr><tr><td>M3</td><td>1734 ± 12</td><td>1.16</td><td>4935</td><td>1736 ± 12</td><td>1.38</td><td>4820</td><td>1740 ± 12</td><td>0.99</td><td>4885</td><td>1740 ± 12</td><td>1.02</td><td>4995</td><td>1739 ± 12</td><td>0.70</td><td>4915</td></tr><tr><td>M4</td><td>1765 ± 10</td><td>1.11</td><td></td><td>1755 ± 10</td><td>1.34</td><td></td><td>1758 ± 10</td><td>0.96</td><td></td><td>1755 ± 10</td><td>0.96</td><td></td><td>1755 ± 10</td><td>0.66</td><td></td></tr><tr><td>M5</td><td>1802 ± 10</td><td>1.08</td><td></td><td>1755 ± 10</td><td>1.34</td><td></td><td>1767 ± 10</td><td>0.92</td><td></td><td>1769 ± 10</td><td>0.91</td><td></td><td>1763 ± 10</td><td>0.62</td><td></td></tr><tr><td>M6</td><td>1808 ± 10</td><td>1.04</td><td>4965</td><td>1792 ± 10</td><td>1.22</td><td>4820</td><td>1794 ± 10</td><td>0.90</td><td>4910</td><td>1794 ± 10</td><td>0.88</td><td>5010</td><td>1791 ± 10</td><td>0.58</td><td>4930</td></tr><tr><td>M7</td><td>1825 ± 10</td><td>0.98</td><td></td><td>1819 ± 10</td><td>1.21</td><td></td><td>1817 ± 10</td><td>0.85</td><td></td><td>1817 ± 10</td><td>0.86</td><td></td><td>1815 ± 10</td><td>0.47</td><td></td></tr><tr><td>M8</td><td>1843 ± 9</td><td>0.95</td><td>4985</td><td>1849 ± 9</td><td>1.20</td><td>4825</td><td>1848 ± 9</td><td>0.82</td><td>4920</td><td>1847 ± 9</td><td>0.85</td><td>5020</td><td>1852 ± 9</td><td>0.46</td><td>4950</td></tr><tr><td>M9</td><td>1871 ± 9</td><td>0.85</td><td>4995</td><td>1868 ± 9</td><td>1.18</td><td>4830</td><td>1864 ± 9</td><td>0.78</td><td>4940</td><td>1870 ± 9</td><td>0.83</td><td>5030</td><td>1873 ± 9</td><td>0.41</td><td>4965</td></tr><tr><td>M10</td><td>1912 ± 9</td><td>0.78</td><td>5025</td><td>1909 ± 9</td><td>0.94</td><td>4905</td><td>1905 ± 9</td><td>0.62</td><td>4950</td><td>1910 ± 9</td><td>0.66</td><td>5070</td><td>1907 ± 9</td><td>0.31</td><td>4980</td></tr><tr><td>Aerial photographs</td><td>1940</td><td>0.75</td><td>5030</td><td>1940</td><td>0.84</td><td>4925</td><td>1940</td><td>0.57</td><td>4965</td><td>1940</td><td>0.53</td><td>5110</td><td>1940</td><td>0.26</td><td>4990</td></tr><tr><td></td><td>1956</td><td>0.71</td><td>5035</td><td>1956</td><td>0.79</td><td>4935</td><td>1956</td><td>0.51</td><td>4985</td><td>1956</td><td>0.42</td><td>5120</td><td>1956</td><td>0.24</td><td>5000</td></tr><tr><td></td><td>1963</td><td>0.68</td><td>5040</td><td>1963</td><td>0.73</td><td>4940</td><td>1963</td><td>0.51</td><td>5000</td><td>1963</td><td>0.39</td><td>5130</td><td>1963</td><td>0.21</td><td>5010</td></tr><tr><td></td><td>1974</td><td>0.65</td><td>5055</td><td>1974</td><td>0.69</td><td>4945</td><td>1974</td><td>0.48</td><td>5005</td><td>1974</td><td>0.35</td><td>5140</td><td>1974</td><td>0.18</td><td>5015</td></tr><tr><td></td><td>1983</td><td>0.59</td><td>5075</td><td>1983</td><td>0.63</td><td>4950</td><td>1983</td><td>0.44</td><td>5010</td><td>1983</td><td>0.31</td><td>5145</td><td>1983</td><td>0.14</td><td>5020</td></tr><tr><td></td><td>1997</td><td>0.49</td><td>5095</td><td>1997</td><td>0.52</td><td>4960</td><td>1997</td><td>0.36</td><td>5060</td><td>1997</td><td>0.23</td><td>5165</td><td>1997</td><td>0.11</td><td>5025</td></tr></tbody></table>

photograph), $V_i$, $S_i$ are the volume and surface area for glacier stage $i$, and $\rho$ is the ratio of ice density to water density (fixed at 0.9).

# 4. RESULTS

## 4.1. LIA chronology on Charquini glaciers

Only the principal results are presented here. Details about the interpretation of moraine dating can be found in Rabatel and others (2005). On the five glaciers the ten most important moraines were considered to determine the chronology of glacier evolution since the LIA maximum. The ten moraines can be found on all five slopes and show the same relative position and morphology (Fig. 1). Five moraines are always prominent (M1, M3, M6, M8 and M9) whereas the others (M2, M4, M5, M7 and M10) are less marked. The morphologies of the prominent moraines are believed to reflect phases of stabilization or small advances, whereas the others represent only short stops in a continuous recession process. By dating the moraines (Table 1), we found that the prominent outer moraine M1, which corresponds to the LIA maximum, was formed during the second half of the 17th century, i.e. between 1648 and 1700, taking into account the five slopes and the margins of error. The moraines M2, M3, M4 and M5, formed during the 18th century, reflect a gradual recession of glaciers, only interrupted by an advance which deposited moraine M3, which dated from 1722–52 depending on the glacier. In the late 18th to early 19th centuries (1781–1818), M6 attests to another advance. Then, during the first half of the 19th century, the glaciers retreated slightly and did not display an advance which could remove the older moraines. The long distance, more than 250 m on average, between the last moraine of the 19th century, M9 (about 1870), and M10 (about 1910), suggests an acceleration of the retreat during the late 19th century.

## 4.2. Ice surface area variations

Throughout the period extending from the mid-17th to the late 20th century, all the glaciers retreated over a distance of 950–1400 m. However, Glaciares Norte and Oeste lost much more (>78% of the maximum LIA area) than Glaciares Sur, Sureste and Noreste (<65%; see Table 1). Expressed in terms of annual means for the five glaciers, the area reduction was $0.15\% a^{-1}$ between M1 and M9, increasing to $0.30\% a^{-1}$ between M9 and M10. Afterwards, the retreat slowed during the first part of the 20th century ($0.18\% a^{-1}$) then accelerated after the 1940s, reaching $0.35\% a^{-1}$. Since the beginning of the 1980s, the retreat rates ($0.52\% a^{-1}$) have been the highest recorded over the last four centuries (Fig. 3). Between the LIA maximum and the late 19th century, the greatest retreats occurred on Glaciares Sur and Oeste. These glaciers lost 30–50% of the area of the LIA maximum while the others lost only 16–25%. Nevertheless, during the 20th century, Glaciares Norte, Noreste and Sureste retreated more than the others (40–57% vs 29–36%).

## 4.3. ELA fluctuations

Between the maximum LIA extent and the late 20th century, the ELA rise determined geometrically using the AAR method can be estimated at $158 \pm 30$ m (mean for the five glaciers) ranging from 115 m (for Glaciar Oeste) to 190 m (for Glaciar Noreste). The ELA rose from 4815–4990 m a.s.l. during the morainic stage M1 to 4960–5165 m a.s.l. during the last decades of the 20th century (Table 1).

Figure 4 shows that between the morainic stages M1 (~1665) and M3 (~1735) the ELA rise is very limited, less than 10 m. Its rise between M1 and M9 (~1870) is of the order of 50 m (ranging from 15 m on Glaciar Sureste to 70 m on Glaciar Noreste).

## 4.4. Volume variations and mass-balance changes

Figure 5 shows the mean annual mass balance, reconstructed for the following four periods: LIA maximum to late 19th
<figure><img src="image_4.png" /><figcaption>**Fig. 3.** Surface area lost for each period for the five Charquini glaciers. The area lost is given as a percentage of the initial area.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Fig. 4.** Fluctuation of the ELA since the LIA maximum. The zero represents the average altitude of the ELA over the whole period.</figcaption></figure>

century, late 19th to early 20th century, early 20th century to 1983 and 1983–97. This information leads to the same conclusions as the reconstructed areas. Between the LIA maximum and the late 19th century, ice loss was low, with a mean mass balance estimated at about –0.10 m w.e. $a^{-1}$. Between M9 and M10 (~1870 and ~1910), the mass balance of the five glaciers was, on average, four times more negative (–0.40 m w.e. $a^{-1}$). For the ~1910–1983 period, deficits were less pronounced, with an average of about 0.25 m w.e. $a^{-1}$. The last two decades of the 20th century marked a shift in the deficits, which increased to 1.36 m w.e. $a^{-1}$ (average for the five Charquini glaciers). This clearly represents the largest mass loss recorded over the past four centuries, and it is also for the shortest time-scale. In addition, we observe that Glaciares Oeste and Norte present the strongest recession for recent years (1983–97) with losses as high as 3 and 2.3 m w.e. $a^{-1}$, respectively. However, Glaciar Oeste is believed to be less representative than the others because of its smaller size (<0.15 km$^2$) and because its low elevation puts it completely in the ablation zone. The three other glaciers facing east and south experienced clearly less negative mass balances that did not exceed –0.66 m w.e. $a^{-1}$.

# 5. DISCUSSION

## 5.1. Comparison with other glaciers in the central Andes

The chronology obtained for the Charquini glaciers indicates that recession of glaciers in the Bolivian Andes since their LIA maximum is in good agreement with that of glaciers in the Peruvian Cordillera Blanca (Solomina and others, in press). After the LIA maximum, in the second half of the 17th century, the Cordillera Blanca glaciers receded nearly continuously, with only minor readvances, each less than that which deposited the outer moraine M1. This recession occurred even during the early 19th century, when many Northern Hemisphere glaciers grew and sometimes

exceeded their 17th-century maximum (Grove, 1988). Surface area losses since the LIA maximum on Charquini glaciers are consistent with those obtained on the near-by south-facing Glaciar Chacaltaya (0.035 km$^2$, 5395–5150 m a.s.l.), which in 1998 had lost 80% of its maximum LIA area (Ramírez and others, 2001). Similarly Hastenrath and Ames (1995a, b) found that Glaciar Yanamarey (0.8 km$^2$, 5200–4600 m a.s.l., facing southwest) in the Cordillera Blanca lost about 53% of its surface area between the LIA maximum (undated) and the early 1980s, consistent with Glaciar Charquini Sur which lost almost the same area (51%) over the same period. The accelerated retreat and very negative mean annual mass balance during the M9–M10 period could be the consequence of a significant climate change. Both support the idea that the LIA ended between about 1870 and 1910. This result is consistent with evidence of rapid glacier decrease in other places of the tropical Andes in Peru and Ecuador, as reported by voyagers and scientists of the late 19th century (Broggi, 1945; Hastenrath, 1981; Ames and Francou, 1995; Francou, 2004).

Several studies of 20th-century recession have been performed in Peru (Cordillera Blanca and Cordillera de Vilcanota) by Broggi (1945), Kinzl (1969), Brecher and Thompson (1993), Hastenrath and Ames (1995a, b) Kaser and Georges (1997) and Thompson and others (2000). In the Cordillera Blanca, glacier retreat accelerated after the 1890s (Kinzl, 1969). Then it slowed down during the first half of the 20th century, with a small but marked readvance in the 1920s. This event was followed by another significant retreat in the 1930s–1940s (Broggi, 1945; Kaser and Georges, 1997; Georges, 2004). During 1950–1970, glaciers retreated very slowly (Hastenrath and Ames, 1995b), followed by a general acceleration (Ames and Francou, 1995; Kaser and Georges, 1997). For the Cordillera de Vilcanota, Brecher and Thompson (1993) and Thompson and others (2000) noted that the Qori Kalis, a glacier issuing from the Quelccaya ice cap, receded between 1963 and 1998. The retreat increased after 1980. The area lost was 8.4 times greater between 1983 and 1991 than during the 1963–78 period, and the retreat
accelerated dramatically during the 1990s. For the Bolivian eastern Cordillera, the recession of the Charquini glaciers can be summarized in four stages: (1) a major retreat starting in the late 19th century; (2) a relative slowdown during the 1910s/1930s; (3) an acceleration since the 1940s; and (4) a very strong recession during the 1980s/1990s. Thus, the consistency of glacier fluctuations observed throughout the tropical Andes between the late 19th and the 20th centuries highlights that glaciers appear to respond to the same climatic signal on a regional scale.

## 5.2. Climatic interpretation

### From the mid-17th to the late 19th century

Whatever their aspect, the Charquini glaciers have evolved in a similar way since their LIA maximum and their retreat has been practically continuous since the second half of the 17th century. The presence of many morainic ridges shows that the retreat was interrupted by several periods of stabilization or small advance, but these advances were never large enough to bury the moraines deposited previously. Being similar in size, morphology and elevation, the main parameter which may determine fluctuations of different magnitudes for these glaciers is their aspect, which controls precipitation amounts, cloudiness and radiative fluxes. Thus, variations from one glacier to another can be interpreted in terms of climate influence.

According to the sensitivity analysis carried out by Hastenrath and Ames (1995a) on Glaciar Yanamarey in Peru and assuming melting to be the only ablation process (with a latent heat of melting $L_m = 33 \times 10^4 \text{ J kg}^{-1}$), the Charquini glaciers' mean annual balance value between the LIA maximum and the late 19th century of -0.10 m w.e. translates to increasing energy absorbed by the glacier surface of about $1 \text{ W m}^{-2}$. This excess energy could have been the consequence of a rise in temperature of $0.15^\circ\text{C}$, an additional humidity of $0.05 \text{ g kg}^{-1}$, a decrease in cloudiness of about 0.1/10, or a combination of these variables. The increase in incoming solar radiation caused by a decrease in cloudiness can explain the strong sensitivity of glaciers to a change in this parameter. However, the resulting decrease in incoming longwave radiation can offset this effect (Ambach, 1974). It is thus difficult to determine the exact influence of changes in cloudiness on the glacier surface energy balance. A change in cloudiness can also affect precipitation and thus the glacier surface albedo (assuming that precipitation falls in mainly solid form on the glacier) but such a small change is hardly quantifiable. However, these two parameters (less cloudiness plus less precipitation) are linked and even with a slight shift they may sum their effects to increase melting at the glacier surface. If the glacier retreat between the LIA maximum and the late 19th century was the result of a $0.15^\circ\text{C}$ rise in temperature, the resulting ELA shift according to Kaser's model linking ELA to climate variables (Kaser, 2001) would have been only 20 m, inconsistent with the 50 m found on Charquini glaciers. In the model of Kaser (2001) a change of the ELA can be interpreted in terms of changes in accumulation, energy balance and temperature or any combination of these variables. Considering change in accumulation only, this 50 m ELA shift between the LIA maximum and the late 19th century (of which 42 m was between M3 (~1735) and M9 (~1870)) could have been the consequence of a decrease in annual accumulation of about 340 mm. Taking into account both temperature and precipitation changes, with a temperature rise of $0.10^\circ\text{C}$

<figure><img src="image_6.png" /><figcaption>**Fig. 5.** Mean annual mass balance reconstructed for each period for the five Charquini glaciers. The value for Glaciar Chacaltaya is shown for the 1983–1997 period for comparison.</figcaption></figure>

(i.e. two-thirds of the sensitivity analysis from Hastenrath and Ames) the decrease in accumulation between M1 (~1665) and M9 (~1870) becomes 210 mm. This change represents a decrease of 20–25% of the present accumulation amount. This is consistent with the results of Thompson and others (1986) based on the Quelccaya ice core that showed an accumulation deficit of about 20% of the present average amounts during the 1720–1860 period.

Furthermore, as noted in section 4, Glaciar es Oeste and Sur exhibited the greatest retreat and the most negative mass balance between the LIA maximum and the late 19th century, a trend which was particularly marked during the 19th century. Precipitation and cloudiness coming from the northeast are crucial for mass balance (section 2). Lower accumulation on the leeward slopes than on the northeast-facing slopes (Rabatel and Mendoza, 2004) leads to a reduced albedo on Glaciar es Sur and Oeste. Since cloudiness follows the same pattern, these slopes experience a higher shortwave radiation influx which, combined with the reduced albedo, causes greater energy absorption by the glacier surface. In the tropics, the energy available to melt snow and ice is particularly high on southern slopes during summer when the sun is at its highest (Francou and others, 1995). The recent energy balance measured year-round on Glaciar Zongo by Wagnon and others (1999, 2001) shows that the consequences of a precipitation/cloudiness deficit during the austral summer are considerable for the glacier. Therefore, considering the importance of precipitation deficits on glacier mass balance nowadays, the assumption that the mid-18th to late 19th centuries could have been dry is reasonable.

Strengthening the findings of Thompson and others (1986) about a dry 19th century, Valero-Garces and others (2003) have pointed out a precipitation deficit during the late 19th to early 20th century from palaeohydrological reconstructions based on sedimentological, geochemical and isotopic records in the Bolivian Altiplano. This evidence is also consistent with the historical ENSO records (Quinn and Neal, 1992; Ortlieb, 2002) which revealed that the last
decades of the 19th century were marked by strong and long-lasting ENSO warm events. As pointed out in section 2, melting accelerates dramatically during the ENSO warm phases. The El Niño occurrence during this period is thus consistent with the accelerated retreat observed on the Charquini glaciers (section 4.2 and 4.4) and throughout the central Andes.

## The 20th century

For the 20th century, a 100 year record (from 1898) of precipitation in La Paz (San Calixto station, 3660 m a.s.l., 20 km south of the Charquini massif) recently analysed by Gioda and others (in press) did not display any significant trend over this period. In consequence, the estimated $-0.4 \text{ m w.e. } a^{-1}$ mean annual mass balance for all the Charquini glaciers over the 20th century, corresponding to an increasing energy at the glacier surface of $4 \text{ W m}^{-2}$, could have resulted principally from a rise in temperature of $0.57^{\circ}\text{C}$ according to the mass-balance sensitivity analysis. According to Kaser's model, the 110 m rise of the ELA (mean of the five glaciers) corresponds to a rise of $0.6^{\circ}\text{C}$, assuming only a change in temperature. Note that both results of air-temperature rise are consistent with the $0.55^{\circ}\text{C}$ warming generally estimated in the literature over the same period (Lean and others, 1995). Thus, in contrast to the mid-17th to late 19th centuries, the recession of glaciers during the 20th century appears to have been the result of increasing air temperature.

Consistent with this hypothesis is the more severe retreat throughout the 20th century observed on Glaciares Norte and Oeste. Since the snouts there were 150 m lower at the end of the 19th century than on southern and western slopes, a rise in temperature, which controls the altitude of the snow/rain limit, might have led to higher ablation rates and consequently to a more rapid retreat of these tongues.

During the last two decades of the 20th century, the acceleration of glacier retreat is here, as elsewhere in the tropical Andes, coincident with the Pacific climate shift of 1976 after which the El Niño frequency and intensity increased significantly (Francou and others, 2003, 2004). The mean annual mass balance over the 1983–97 period was about $-0.51 \text{ m w.e. } a^{-1}$ for Glaciares Sur, Noreste and Sureste and as high as $-2.64 \text{ m w.e. } a^{-1}$ for Glaciares Norte and Oeste. This represents an increase in energy available at the glacier surface of about $5 \text{ W m}^{-2}$ for Glaciares Sur, Sureste and Noreste and $28 \text{ W m}^{-2}$ for Glaciares Norte and Oeste. Depending on the slope, the possible causes of the rapid glacier retreat in the late 20th century appear to be a rise of $0.8–3^{\circ}\text{C}$ in air temperature (enhancing sensible heat transfer to the glacier), an increase in specific humidity of $0.3–1.5 \text{ kg}^{-1}$ (reducing sublimation and increasing melting), a decrease in cloudiness of less than 0.5/10 to 1/10 (increasing solar radiation) or any combination of these processes.

## 5.3. The imbalance of Charquini glaciers

The assumption that moraine deposits represent glaciers in an equilibrium state is valid for small glaciers. It is therefore appropriate to determine the ELAs 'geometrically' from glacier areas between the LIA maximum and the first part of the 20th century. However, the method is spurious when glaciers are in a constant and rapid recession and hence markedly unbalanced, as occurred during the last decades of the 20th century. In this situation, the snout fluctuations

reflect ablation conditions in the lower zone rather than oscillations of the ELA, which are generally located at the top of the glacier or above. In that case a comparison of the ELAs calculated using the AAR method with those determined from field observations (ablation stakes) shows how unbalanced these glaciers are in relation to the present climate and how rapid their shrinkage will be in the near future (Ramírez and others, 2001). From ablation stakes during the 1991–2004 period, the ELA$_0$ is observed on Glaciar Zongo (facing south to east) and Glaciar Chacaltaya (facing south) at 5215 m a.s.l. (Soruco and others, 2005). Calculated for glaciers of similar aspect on the Cerro Charquini using the AAR method for 1997, the ELA is found at 4960–5095 m a.s.l. That is a difference between direct measurements and the AAR calculation of about 120 m for the southern slopes and about 250 m for the southeastern slopes. This average difference of 175 m represents the imbalance of the Charquini glaciers under current climatic conditions. Due to the small gap existing between their mean ELA and the maximum elevation of their upper reaches, these glaciers could recover equilibrium in the future only by a drastic reduction in their area. This implies that the Charquini glaciers, and in general all the glaciers of the Cordillera Real extending below $\sim 5300 \text{ m a.s.l.}$, are unbalanced under the present climate and could disappear in the near future.

## 6. CONCLUSIONS

A new chronology and a quantification of the glacier recession since the LIA maximum have been presented in this study from the Charquini area in the Bolivian Andes. Surface area, ELA, volume and mass balance were reconstructed using well-preserved moraines dated by lichenometry and using photogrammetry for the last six decades of the 20th century.

The LIA culminated during the second half of the 17th century and glaciers subsequently retreated nearly continuously. Retreat slowed down during the early 19th century but no important readvances occurred.

The major acceleration of glacier recession observed between M9 and M10 suggests the LIA ended between about 1870 and 1910, consistent with evidence found throughout the Andes and in many places in the Northern Hemisphere.

The retreat between the LIA maximum and the late 19th century was probably associated with a significant decrease in precipitation and cloudiness, particularly during the 19th century.

In contrast, a sensitivity analysis shows that glacier retreat during the 20th century has been principally due to a rise in air temperature of about $0.57^{\circ}\text{C}$. In order to reverse the present deficit and to balance these glaciers, a decrease in energy available for melting ice of about $28 \text{ W m}^{-2}$ would be required. This exceeds the current climatic variability.

The AAR method suggests the ELA might have been on average 158 m lower during the second half of the 17th century (LIA maximum) than during the last decades of the 20th century, assuming that current glacier geometry corresponds to a steady state of the glaciers.
The ELA measured on nearby Glaciar Zongo for the 1991–2004 period (i.e. 5215 m a.s.l.) differs from the ELA geometric computation based on the AAR method by 175 m; this difference shows that the current glacier extent does not correspond to a balanced situation. It can be deduced that the Charquini glaciers, as well as all the small glaciers of the Cordillera Real below 5300 m a.s.l., are not in equilibrium with the present climate and that they could disappear in the next few decades. This confirms the assertion from Glaciar Chacaltaya (Ramírez and others, 2001).

## ACKNOWLEDGEMENTS

We thank all those who took part in the field measurements (lichenometry and ground topography), the French Eclipse program (INSU-CNRS) managed by V. Masson-Delmotte (LSCE, CEA-CNRS), and P. Naveau (LSCE-CNRS) and D. Grancher (LGP-CNRS) for analysis using the Lichen Power® software. We also thank R. Naruse, editor of this paper, and two referees, H. Kerschner and B. Mark, for their helpful comments and suggestions.

## REFERENCES

Ambach, W. 1974. The influence of cloudiness on the net radiation balance of a snow surface with high albedo. *J. Glaciol.*, **13**(67), 73-84.

Ames, A. and B. Francou. 1995. Cordillera Blanca. Glacaires en la historia. *Bull. Inst. fr. étud. andin.*, **24**(1), 37-64.

Brecher, H.H. and L.G. Thompson. 1993. Measurement of the retreat of Qori Kalis glacier in the tropical Andes of Peru by terrestrial photogrammetry. *Photogramm. Eng. Rem. Sens.*, **59**(6), 1017-1022.

Broggi, J.A. 1945. La desglaciation actual de los Andes del Perú. *Bol. mus. hist. nat.* "Javier Prado", **9**(34-35), 222-248.

Clapperton, C.M. 1983. The glaciation of the Andes. *Quat. Sci. Rev.*, **2**(2-3), 83-155.

Cooley, D., P. Naveau and V. Jomelli. In press. A Bayesian hierarchical extreme value model for lichenometry. *Environmetrics*.

Francou, B. 2004. Andes del Ecuador: los glacaires en la época de los viajeros, siglos XVIII a XX. In *Dollfus, O.*, J.P. Deler and E. Mesclier, eds. *Los Andes el reto del espacio mundo andino homenaje*. Lima, Instituto Francés de Estudios Andinos/Instituto de Ecología Política, 137-152.

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.

Francou, B., M. Vuille, P. Wagnon, J. Mendoza and J.E. Sicart. 2003. Tropical climate change recorded by a glacier in the central Andes during the last decades of the twentieth century: Chacaltaya, Bolivia, 16°S. *J. Geophys. Res.*, **108**(D5), 4154. (10.1029/2002JD002959.)

Francou, B., M. Vuille, V. Favier and B. Cáceres. 2004. New evidence for an ENSO impact on low-latitude glaciers: Antizana 15, Andes of Ecuador, 0°28'S. *J. Geophys. Res.*, **109**(D18), D18106. (10.1029/2003JD004484.)

Georges, C. 2004. 20th century glacier fluctuations in the tropical Cordillera Blanca, Peru. *Arct. Antarct. Alp. Res.*, **36**(1), 100-107.

Gioda, A., Y. L'Hote, L.A. Drake, J. Ronchail and B. Pouyaud. In press. Analyse et variabilité temporelle d'une longue série de pluies des Andes en relation avec l'Oscillation Australe (La Paz, 3658 m, 1891-2000). In *Demarée, G.*, ed. *TCMH-2001*. Brussels, Royal Meteorological Institute of Belgium.

Gouze, P., J. Argollo, J.F. Saliege and M. Servant. 1986. Interprétation paléoclimatique des oscillations des glaciers au cours des

20 derniers millénaires dans les régions tropicales; exemple des Andes boliviennes. *CR Acad. Sci. (Paris)*, **303**(Série II), 219-224.

Grove, J.M. 1988. *The Little Ice Age*. London and New York, Methuen.

Hastenrath, S. 1981. *The glaciation of the Ecuadorian Andes*. Rotterdam, A.A. Balkema.

Hastenrath, S. and A. Ames. 1995a. Diagnosing the imbalance of Yanamarey Glacier in the Cordillera Blanca of Peru. *J. Geophys. Res.*, **100**(D3), 5105-5112.

Hastenrath, S. and A. Ames. 1995b. Recession of Yanamarey Glacier in Cordillera Blanca, Peru, during the 20th century. *J. Glaciol.*, **41**(137), 191-196.

Innes, J.L. 1985. Lichenometry. *Prog. Phys. Geog.*, **9**(2), 187-254.

Jordan, E. 1991. *Die Gletscher der bolivianischen Anden: eine photogrammetrisch-kartographische Bestandsaufnahme der Gletscher Boliviens als Grundlage für klimatische Deutungen und Potential für die wirtschaftliche Nutzung*. Stuttgart, Franz Steiner Verlag. (Erdwissenschaftliche Forschung 23.) [In German with English summary.]

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

Kaser, G. and C. Georges. 1997. Changes of the equilibrium-line altitude in the tropical Cordillera Blanca, Peru, 1930-50, and their spatial variations. *Ann. Glaciol.*, **24**, 344-349.

Katz, R.W., M. Parlange and P. Naveau. 2002. Extremes in hydrology. *Adv. Water Resour.*, **25**(8), 1287-1304.

Kinzl, H. 1969. La glaciacion actual y pleistocenica en los Andes centrales. *Bol. Soc. Geog. Lima*, **89**, 89-100.

Lean, J., J. Beer and R. Bradley. 1995. Reconstruction of solar irradiance since 1610: implications for climate change. *Geophys. Res. Lett.*, **22**(23), 3195-3198.

Le Roy Ladurie, E. 1967. *Histoire du climate depuis l'an Mil*. Paris, Flammarion.

Luckman, B.H. 2000. The Little Ice Age in the Canadian Rockies. *Geomorphology*, **32**, 357-384.

Naveau, P., V. Jomelli, D. Cooley and A. Rabatel. In press. Modelling uncertainties in lichenometry studies with an application: the tropical Andes (Charquini Glacier, Bolivia). *Arct. Antarct. Alp. Res.*

Ortlieb, L. 2002. Manifestations historiques du phénomène El Niño en Amérique du Sud depuis 16ème siècle. *Houille Blanche*, **6/7**, 115-120.

Quinn, W.H. and V.T. Neal. 1992. The historical record of El Niño events. In Bradley, R.S. and P.D. Jones, eds. *Climate since A.D. 1500*. London, etc., Routledge, 623-648.

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, Grenoble.)

Rabatel, A. and J. Mendoza. 2004. Glacaires Zongo, Chacaltaya y Charquini Sur (16°S, Bolivia). *Mediciones meteorológicas, hidrológicas y glaciológicas, año hidrológico 2002-2003*. La Paz, Great-Ice.

Rabatel, A., V. Jomelli, P. Naveau, B. Francou and D. Grancher. 2005. Dating of Little Ice Age glacier fluctuations in the tropical Andes: Charquini glaciers, Bolivia, 16°S. *CR Geoscience*, **337**(15), 1311-1322.

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.

Reimer, P.J. and 27 others. 2004. IntCal04 terrestrial radiocarbon age calibration, 0-26 cal kyr BP. *Radiocarbon*, **46**, 1029-1058.

Seltzer, G.O. 1992. Late Quaternary glaciation of the Cordillera Real, Bolivia. *J. Quat. Sci.*, **7**(2), 87-98.

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.)

Solomina, O., V. Jomelli, G. Kaser, A. Ames and B. Pouyaud. In press. Little Ice Age moraines in the Cordillera Blanca: lichenometric data replication. *Global Planet. Change*.
Soruco, A. and 13 others. 2005. Glaciares Zongo-Chacaltaya-Charquini Sur – Bolivia 16°S. *Mediciones glaciológicas, hidrológicas y meteorológicas, año hidrológico 2003-2004*. La Pas, Great-Ice.

Thompson, L.G., E. Mosley-Thompson, W. Dansgaard and P.M. Grootes. 1986. The Little Ice Age as recorded in the stratigraphy of the tropical Quelccaya ice cap. *Science*, **234**(4774), 361–364.

Thompson, L.G., E. Mosley-Thompson and K.A. Henderson. 2000. Ice-core palaeoclimate records in tropical South America since the Last Glacial Maximum. *J. Quat. Sci.*, **15**(4), 377–394.

Valero-Garces, B.L., A. Delgado-Huertas, A. Navas, L. Edwards, A. Schwalb and N. Ratto. 2003. Patterns of regional hydrological variability in central-Southern Altiplano (18°-26° S) lakes during the last 500 years. *Palaeogeogr., Palaeoclimatol., Palaeoecol.*, **194**, 319-338.

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

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.

MS received 14 December 2004 and accepted in revised form 22 October 2005