Journal cover Journal topic
Hydrology and Earth System Sciences An interactive open-access journal of the European Geosciences Union
Journal topic
Hydrol. Earth Syst. Sci., 23, 657-668, 2019
https://doi.org/10.5194/hess-23-657-2019
Hydrol. Earth Syst. Sci., 23, 657-668, 2019
https://doi.org/10.5194/hess-23-657-2019

Review article 06 Feb 2019

Review article | 06 Feb 2019

# Modelling Lake Titicaca's daily and monthly evaporation

Modelling Lake Titicaca's daily and monthly evaporation
Ramiro Pillco Zolá1, Lars Bengtsson2, Ronny Berndtsson2, Belen Martí-Cardona3, Frederic Satgé4, Franck Timouk5, Marie-Paule Bonnet6, Luis Mollericon1, Cesar Gamarra7, and José Pasapera7 Ramiro Pillco Zolá et al.
• 1Instituto de Hidráulica e Hidrología, Universidad Mayor de San Andrés, La Paz, Bolivia
• 2Division of Water Resources Engineering and Center for Middle Eastern Studies, Lund University, Lund, Sweden
• 3Department of Civil and Environmental Engineering, University of Surrey, Guildford, UK
• 4CNES, UMR HydroSciences, Univeristy of Montpellier, Place E. Bataillon, 34395 Montpellier CEDEX 5, France
• 5Laboratoire GET UMR5563, CNRS, IRD, Université Paul Sabatier, OMP, Toulouse, France
• 6IRD, UMR Espace-Dev, Maison de la télédétection, 500 Rue JF Breton, 34093 Montpellier CEDEX 5, France
• 7IMARPE, Instituto del Mar del Perú, Puno, Peru
Abstract

Lake Titicaca is a crucial water resource in the central part of the Andean mountain range, and it is one of the lakes most affected by climate warming. Since surface evaporation explains most of the lake's water losses, reliable estimates are paramount to the prediction of global warming impacts on Lake Titicaca and to the region's water resource planning and adaptation to climate change. Evaporation estimates were done in the past at monthly time steps and using the four methods as follows: water balance, heat balance, and the mass transfer and Penman's equations. The obtained annual evaporation values showed significant dispersion. This study used new, daily frequency hydro-meteorological measurements. Evaporation losses were calculated following the mentioned methods using both daily records and their monthly averages to assess the impact of higher temporal resolution data in the evaporation estimates. Changes in the lake heat storage needed for the heat balance method were estimated based on the morning water surface temperature, because convection during nights results in a well-mixed top layer every morning over a constant temperature depth. We found that the most reliable method for determining the annual lake evaporation was the heat balance approach, although the Penman equation allows for an easier implementation based on generally available meteorological parameters. The mean annual lake evaporation was found to be 1700 mm year−1. This value is considered an upper limit of the annual evaporation, since the main study period was abnormally warm. The obtained upper limit lowers by 200 mm year−1, the highest evaporation estimation obtained previously, thus reducing the uncertainty in the actual value. Regarding the evaporation estimates using daily and monthly averages, these resulted in minor differences for all methodologies.

1 Introduction

Lake Titicaca, the largest freshwater lake in South America, is located in the endorheic Andean mountain range plateau Altiplano, straddling the border between Peru and Bolivia (Fig. 1). The lake plays an essential role in shaping the semiarid Altiplano climate; feeding the downstream Desaguadero River and Lake Poopó (Pillco and Bengtsson, 2006; Abarca-del-Río et al., 2012); and supplying the inhabitants with water resources for domestic, agricultural, and industrial use (Revollo, 2001). Anthropogenic pressure on the Altiplano water resources has increased during the last decades due to population growth and increased evapotranspiration losses (FAO, 2011; Canedo et al., 2016; Satgé et al., 2017) as well as to industrial pollution (UNEP, 1996; CMLT, 2014). The challenge of managing water resources in the Altiplano Basin is further exacerbated by climate conditions; annual rainfall is highly variable (Garreaud et al., 2003), while warming in this region exceeds the average global trend (López-Moreno et al., 2015), which is expected to intensify the evaporation from the lake surface and the evapotranspiration losses from the whole basin. The combined impact of these pressures becomes evident at the downstream end of the system, where Lake Poopó is situated. In recent years this lake suffered extreme water shortages, including its complete drying out in December 2015 (Satgé et al., 2017).

Figure 1Lake Titicaca and the TDPS system within the Altiplano.

Lake Titicaca has a large surface area of about 8500 km2 on average. Over a certain water surface level, the lake spills out at the south-eastern end and feeds the Desaguadero River. However, the major water output from Lake Titicaca is due to evaporation, which accounts for approximately 90 % of the losses (Roche et al., 1992; Pouyaud, 1993; Talbi et al., 1999; Delclaux et al., 2007). In recent years, Lake Titicaca's level dropped below the outlet threshold for some periods. Thus, a small increase in evaporation or decrease in precipitation may turn the lake into a closed system with no outflow.

Since evaporation dominates the water balance in Lake Titicaca, it is essential to improve the knowledge of the lake's evaporation. This is especially important in light of anthropogenic pressure and due to the evident strong global warming that this region experiences. Previous studies of Lake Titicaca's evaporation have all been based on monthly meteorological observations. Due to the importance of lake evaporation, detailed calculations using daily as well as monthly observations may be necessary. Consequently, this paper investigates different methods for calculating evaporation using both daily and monthly data; in addition we discussed the possibility for the appropriate evaporation models at both timescales to be used on the study the climatic functioning and sensitivity. The main problem with Lake Titicaca's evaporation estimation is the lack of high-resolution temporal data. Taking into account only the mass transfer models for different timescale, Singh and Xu (1997) calculated monthly evaporation. However, doing the same calculations on a daily basis could give radically different results. For both timescales, the evaporation estimation could be more sensitive to vapour pressure. On the other hand, random errors in input data could have a significant effect on evaporation estimation at a monthly scale rather than at a daily scale (Singh et al., 1997).

Lake Titicaca's surface water is cold, with a temperature that remains 12–17 C throughout the year, and below 40 m depth the temperature is almost constant (Richerson et al., 1977). The water is usually warmer than the air during the daytime, which means that the air immediately above the lake is unstable. The air temperature shows large diurnal variations, often exceeding 15 C in summer. At an average terrain elevation above 4000 m a.s.l. the solar radiation is strong and the atmospheric pressure is low, which means that the ratio between sensible and latent heat flux (Bowen ratio) is lower than at sea level. To determine the evaporation rate using the aerodynamic mass transfer approach, the atmospheric vapour pressure and surface temperature must be known. Furthermore, a wind function must be used, because the atmosphere over Lake Titicaca is unstable most of the time. This means that the wind function may be different from the function used for most other lakes.

It can generally be assumed that during a year the lake water temperature returns to the value at the beginning of the year. Thus, for the heat balance, it is sufficient to know the annual net radiation, provided that the sensible heat flux can be estimated from the constant Bowen ratio. When using the method for shorter time periods, the time variation of the lake water temperature profile must be known. The heat balance approach and the aerodynamic method can be combined. The Penman method is such a combined approach. A wind function must also be included in this approach.

## Previous investigations

One of the first evaporation studies for Lake Titicaca applied the water balance method using measurements for the period 1956–1973 (Carmouze et al., 1977) and estimated a mean annual lake evaporation of 1550 mm year−1. Taylor and Aquize (1984) applied a bulk transfer approach for a shorter period and determined the lake evaporation to be 1350 mm year−1. The largest reported annual evaporation is from using the energy balance approach. Richerson et al. (1977) found the lake evaporation to be equal to 1900 mm year−1. Later, Carmouze (1992) used the same approach and found the lake evaporation to be 1720 mm year−1. Using observations for the period 1965–1983 and the water balance method, Pouyaud et al. (1993) found the mean annual evaporation to be equal to about 1600 mm year−1. Thus, the mean annual evaporation has been estimated in the range 1350–1900 mm year−1. While the precipitation can vary much from year to year, the large range of calculated annual evaporation, 1350 to 1900 mm year−1, is likely to be the result of uncertainties in the evaporation estimations and in the temporal resolution of the measurements. Recently, Delclaux et al. (2007) studied the evaporation from Lake Titicaca using in situ pan-evaporation measurements, energy balance, mass transfer and the Penman methods. They concluded that the mean annual evaporation may be about 1650 mm year−1, with low seasonal variation between 135 mm in July (winter) and 165 mm in November (summer).

Table 1Climatological and hydrological components of Lake Titicaca for 1966–2011.

This study applies the methodologies mentioned above using the frequent and accurate hydro-meteorological measurements acquired at Lake Titicaca in 2015 and 2016, with the aim of reducing the uncertainty in evaporation estimates and evaluating the effect of using daily records instead of monthly ones.

2 Study area

Lake Titicaca is a unique biosphere due to its large depth and volume, high elevation, and tropical latitude. It is located in the northern part of the Peruvian–Bolivian Altiplano, between latitudes of 1545 S and 6925 W. It is surrounded by the eastern and western Andean Mountains. The total Lake Titicaca basin area, including the lake itself, is close to 57 000 km2, with a mean elevation higher than 4000 m a.s.l. The outlet sill is at 3807 m a.s.l. The lake volume is about 903 km3, with a corresponding mean depth of 105 m (Boulange and Aquize, 1981; Wirrmann, 1992). The only outlet is the Desaguadero River, which ends in the shallow Lake Poopó. The modern Lake Titicaca consists of the major Titicaca lake (Lake Chuquito), which is 284 m at the deepest point, and the smaller Titicaca lake (Lake Huiñamarca). The latter lake represents 1200 km2, with a maximum depth of 35 m below the spill level. The threshold between the two lake basins is 19 m below the spill level (see Fig. 2). The lake is described by Dejoux and Iltis (1992) and in the Encyclopedia of Lakes and Reservoirs edited by Bengtsson et al. (2012).

Figure 2Lake Titicaca basin with sub-basins and major and smaller lakes.

## 2.1 Hydrology

The Lake Titicaca watershed is a part of the TDPS system (Titicaca, Desaguadero, Poopó and Salares) within the Altiplano (Revollo, 2001). Lake Poopó is considered a terminal lake, with only one discharge event into the downstream Coipasa salt pan occurring in the last century (Pillco and Bengtsson, 2006). The basin of Lake Titicaca itself includes the sub-basins Katari, Coata, Huancané, Huaycho, Ilave, Illpa, Keka-Achacachi, Ramis and Suchez. The largest is the Ramis River basin, with an area of 15 000 km2, representing 30 % of the total basin (Fig. 2). The mean flow of the Ramis River for the period 1965–2011 was 72 m3 s−1. The mean outflow from Lake Titicaca through the Desaguadero River for the same period was 35 m3 s−1. During those 50 years, Lake Titicaca experienced large changes in water level, with a mean close to 3808.1 m a.s.l., which is about 1 m above the outlet threshold (Pillco and Bengtsson, 2006). From its low to high water level, the Lake Titicaca water surface area might change from a minimum of 7000 to a maximum of 9000 km2.

## 2.2 Climate

The northern part of the Altiplano is semiarid, while the southern part, including the biggest salt pans in the world, is arid (TDPS, 1993). The climate is further characterized by a short wet season (December–March) and a long dry season (April–November; Garreaud et al., 2003). The average precipitation over the Lake Titicaca basin is about 800 mm year−1, out of which more than 70 % fall during the wet season (Garreaud et al., 2003). Over the lake, annual precipitation is assumed to vary from 1200 mm year−1 in the central part to 800 mm year−1 along the shores (TDPS, 1993). January is the wettest (about 180 mm month−1) and July the driest (less than 10 mm) month. In the central and southern parts of the Altiplano, the total annual precipitation is about 350 mm (Roche et al., 1992; Pillco and Bengtsson, 2006) and is less than 200 mm over the Salares in the southernmost areas (Satgé et al., 2016). The seasonal variability of precipitation in the basin is related to changes in the upper troposphere circulation. During the Austral summer, an upper-level cyclone is established to the south-east of the central Andes. The Bolivian high brings easterly winds and allows influx of moisture from the central continent over the plateau during periods, intensifying the precipitation (Garreaud, 1999; Vuille et al., 2000).

The daily air temperature over Lake Titicaca is rather constant throughout the year, usually varying between 7 and 12 C but sometimes up to 20 C in summer. The Titicaca region is more humid than the more southern parts of the Altiplano. The relative humidity varies between 52 % to 68 % as a monthly average, with diurnal variation between 33 % and 80 %. According to Carmouze (1992), the dominant wind on the lake is in the north-west to south-east direction, with mean monthly wind velocity close to 2 m s−1, rarely reaching 5 m s−1 at the daily time step. The general climate and hydrology are summarized in Table 1. The total river inflow was estimated through a representative area approach based on the Ramis River discharge.

3 Methods

## Lake evaporation models

Four evaporation estimation conventional methods were applied in this study: water balance, energy balance, mass transfer and the Penman method. These approaches have previously been used by other researchers to estimate Lake Titicaca's evaporation at a monthly time step (Carmouze, 1992; Pouyaud, 1993; Delclaux et al., 2007). The methods are briefly described as follows.

## 3.1 Energy balance approach

The energy balance approach (Maidment, 1993) which comes from the integral energy balance equation of a reference volume at the air–water interface, the evaporation component in terms of latent heat flux is

$\begin{array}{}\text{(1)}& \mathit{\lambda }E=\frac{{R}_{\mathrm{n}}-{Q}_{\mathrm{heat}}}{\mathrm{1}+\mathit{\beta }},\end{array}$

where λ is the latent heat vaporization (J kg−1), E is evaporation rate (mm day−1), Rn is net radiation (W m−2), Qheat is heat storage within the water (W m−2) and β is the Bowen ratio (Bowen, 1926):

$\begin{array}{}\text{(2)}& & \mathit{\beta }=\mathit{\gamma }\frac{{T}_{\mathrm{w}}-{T}_{\mathrm{a}}}{{e}_{\mathrm{w}}-{e}_{\mathrm{a}}},\text{(3)}& & \mathit{\gamma }=\frac{{c}_{\mathrm{p}}{p}_{\mathrm{a}}}{\mathrm{0.622}\mathit{\lambda }},\end{array}$

where γ is the psychrometric constant (mbar C−1); Tw and Ta are the surface water and air temperature (C), respectively; ew and ea are the water surface and air vapour pressures (Pa), respectively; cp is the specific heat capacity (J kg−1C−1); and Pa is the atmospheric pressure (kg Pa). The psychrometric constant, and thus also the Bowen ratio, are lower at this high altitude than at sea level. The net radiation is the sum of net short-wave and net long-wave radiation. The net short-wave radiation is Rs (1− albedo), where Rs is the solar radiation reaching the lake (W m−2). The atmospheric long-wave radiation as well as the back radiation are computed from Stefan's law. The emissivity of the water is well known and was set to 0.98, and the emissivity atmosphere must be known. The emissivity of the atmosphere depends on humidity, temperature and cloudiness. The atmospheric emissivity accounting for clouds was proposed by Crawford and Duchon (1999):

$\begin{array}{}\text{(4)}& & {\mathit{\epsilon }}_{e}=\left(\mathrm{1}-s\right)+s{\mathit{\epsilon }}_{o}\left({T}_{\mathrm{a}},{e}_{\mathrm{a}}\right),\text{(5)}& & s=\frac{{R}_{\mathrm{s}}}{{R}_{\mathrm{s},o}},\text{(6)}& & {R}_{\mathrm{s},o}={R}_{\mathrm{a}}{e}^{\left(\frac{-\mathrm{0.0018}{P}_{\mathrm{a}}}{{K}_{\mathrm{t}}\mathrm{sin}{\mathit{\varphi }}_{\mathrm{24}}}\right)},\text{(7)}& & {\mathit{\epsilon }}_{o}=\mathrm{1.18}{\left(\frac{{e}_{\mathrm{a}}}{{T}_{\mathrm{a}}}\right)}^{\frac{\mathrm{1}}{\mathrm{7}}},\end{array}$

where s is the proxy cloudiness defined as the ratio of measured incoming solar radiation Rs (W m−2) to the solar radiation received for the clear-sky conditions Rs,o (W m−2), and εo is the emissivity in the clear-sky condition, which is determined from the vapour pressure ea expressed in hPa and Ta temperature in Kelvin. The Φ24 is the mean daily sun angle. The constant 1.18 describes the attenuation defined for the region according to Lhomme et al. (2007). The extraterrestrial solar radiation Ra (W m2) is determined as a function of local latitude and altitude and time of year, using the turbidity coefficient Kt=0.85.

The energy equation is fairly easy to use over a full year, since the lake water usually returns to its initial state when computations were started or when Qheat equals zero (W m−2). When using the approach over shorter periods the variation of the water temperature in the lake must be accounted for. In Eq. (1), the change of heat storage is included. From temperature water profiling observations, it was assumed that the water temperature below 40 m does not change from month to month. The temperature, Tw, above this mixing depth, hmix (m), changes but remains almost homothermal after convective mixing during the night (Richerson et al., 1977), which also is corroborated by our own field investigations. Thus, the change of heat content can be estimated from measured surface temperature:

$\begin{array}{}\text{(8)}& {Q}_{\mathrm{heat}}=\mathit{\rho }{c}_{\mathrm{p}}\frac{{V}_{\mathrm{mix}}}{{A}_{\mathrm{lake}}}\frac{\partial {T}_{\mathrm{w}}}{\partial t},\end{array}$

where ρ is density of water (kg m−3), cp is the specific heat capacity of water and Vmix is the volume above the mixing depth (km−3). Carmouze et al. (1992) introduced the concept of the exchange of heat between surface and deep water using the energy balance concept. The results of Carmouze et al. (1992) were compared to the calculation results in the present study.

## 3.2 Mass transfer approach

The mass transfer aerodynamic approach is used in various models based on Dalton's law (Dalton, 1802). The latent heat transfer is related to the vapour pressure deficit. Most often a linear wind function is used (e.g. Carmouze et al., 1992):

$\begin{array}{}\text{(9)}& E=\left(a+bU\right)\left({e}_{\mathrm{w}}-{e}_{\mathrm{a}}\right),\end{array}$

where E is the evaporation rate, U is wind velocity (m s−1) and ewea is the vapour pressure deficit (mbar). The parameter a accounts for unstable atmospheric conditions. Carmouze et al. (1992) used a=0.17 (mm mbar−1 day−1) and b=0.30 (mm mbar−1 s m−1).

## 3.3 Penman approach

The Penman equation is a combination of energy balance and mass transfer used for evaluating open water evaporation (Penman, 1948):

$\begin{array}{}\text{(10)}& E=\frac{\mathrm{\Delta }}{\mathrm{\Delta }+\mathit{\gamma }}\frac{\left({R}_{\mathrm{n}}-{Q}_{\mathrm{heat}}\right)}{\mathit{\lambda }\mathit{\rho }}+\frac{\mathit{\gamma }}{\mathrm{\Delta }+\mathit{\gamma }}c\left({a}^{\prime }+{b}^{\prime }U\right)\left({e}_{\mathrm{s}}-{e}_{\mathrm{a}}\right),\end{array}$

where E is open water evaporation. The slope of the water pressure–temperature curve is denoted by Δ(K Pa s−1), and esea is the saturation deficit of the air (K Pa); here ea is dependent on the relative humidity (%). Delclaux et al. (2007) applied the Penman equation to Lake Titicaca using ${a}^{\prime }=\mathrm{0.26}$, ${b}^{\prime }=\mathrm{0.14}$ and c=1 after optimizing (mm day−1 mbar).

## 3.4 Lake water balance model

The water balance approach was applied to calculate water levels in Lake Titicaca in a previous study by Pillco and Bengtsson (2007). The water balance is

$\begin{array}{}\text{(11)}& {A}_{\mathrm{lake}}\frac{\partial h}{\partial t}=\left(P-E\right){A}_{\mathrm{lake}}+{Q}_{\mathrm{in}}-{Q}_{\mathrm{out}},\end{array}$

where $\partial h/\partial t$ represents change in water depth; P is precipitation on the lake (mm); E is evaporation from open water (mm day−1); Alake is water surface of the lake (km−2), which is a function of depth; Qin is inflow to the lake; and Qout represents the outflow from the lake (m3 s−1). Computations were carried out at a monthly timescale for two periods, one for 1966–2011 and another for 2015–2016. As already pointed out, the most reliable method of computing evaporation over long periods is probably the water balance method. However, the computation only can be general, since the inflow to Lake Titicaca is not measured in all rivers.

## 3.5 Possible errors when using monthly averages

The evaporation during individual days is not important for the water balance but is only important over longer periods like months. However, since the equations for calculating evaporation are not linear, the monthly evaporation computed from monthly mean meteorological data may differ from what is found when data with a higher time resolution are used. In the aerodynamic approach the wind speed is multiplied by the vapour deficit. The energy balance approach includes the Bowen ratio, which may differ from day to day and can even be negative for certain periods. If high atmospheric vapour pressure is related to strong winds, the aerodynamic equation using monthly means can yield lower evaporation estimates than when daily values are used. This is further discussed below. The Bowen ratio changes during a month. When the net radiation is large, the air temperature is likely to be rather high but is not necessarily related to high vapour pressure. For this situation, the Bowen ratio is relatively high, and the computed evaporation is higher than it would have been using a constant monthly Bowen ratio. This means that when using monthly averages, the computed evaporation will tend to be low.

4 Instrumentation and data

For this study, hydro-meteorological parameters and water surface temperature were measured near Lake Titicaca's bank for 24 consecutive months (2015–2016). Vertical lake temperature profiles were also acquired periodically. Observations were taken at 15 min intervals. These records were averaged to daily and monthly values. A Campbell Scientific research-grade automatic weather station (AWS) was installed at the Isla de la Luna (latitude 160159${}^{\prime \prime }$, longitude 690401${}^{\prime \prime }$), near the shore of Lake Titicaca (Fig. 3). The AWS was equipped with a rain gauge sensor, a CS215 probe for measuring relative humidity and air temperature, an A100R vector anemometer and W200P wind vane to measure wind speed, and a Skye SP1100 pyranometer for solar radiation measurement. The surface water temperature was taken from Juli, Puno (latitude 161258${}^{\prime \prime }$, longitude 692731${}^{\prime \prime }$), at a distance of 42 km from Isla de la Luna. A handheld thermometer was used to measure water surface temperature at 8 h intervals at approximately 60 m from the shoreline. An increased daily surface recorded temperature is representative of heat storage changes.

Hydrological data, such as inflow to the lake, were observed at the outlet of the Ramis River. The outflow through the Desaguadero River was observed at Aguallamaya. This is 40 km downstream of the lake outlet. However, there are only a few tributaries between the lake and this point that may contribute to the data uncertainty. The water level was observed at Huatajata at the daily time step, shown as depth in Table 1. Additional lake water temperature soundings were carried out close to Isla de la Luna (latitude 163000${}^{\prime \prime }$, longitude 691510${}^{\prime \prime }$) for specific days during the summer, spring and winter of the study period, using the Hydrolab DS5 multiparameter data sonde. The sounding reached a maximum lake depth of 95 m, with a water temperature recording for each 5 cm at the surface and each 0.5 m below 1 m depth.

Long-term monthly temperature and wind observations from 1960 onwards were available from the Copacabana weather station mentioned above (Fig. 3; Table 1). The El Alto station observations, 50 km from Copacabana, were used to fill 2.5 % of the missing wind data for the period. The monthly precipitation on the lake was determined using the rain gauge at Copacabana and Puno on the lake shore. The total inflow from all rivers was estimated from a representative area approach assuming the specific run-off to be the same for all rivers entering into Lake Titicaca. The long-term outflow from the lake was measured at the outlet of the lake and treated by Gutiérrez and Molina (2014).

Figure 3Location of observation points.

Tables 2 and 3 summarize hydrological and meteorological measurements used in this study. Subscripts for vapour pressure are “w” for water, “s” for saturated air vapour pressure and “a” for actual air vapour pressure. The computed variables required for evaporation calculations are given in Table 4.

Table 2Monthly mean of hydrological variables observed during the 2015–2016 period.

Table 3Monthly averages of main climatic variables observed during the 2015–2016 period.

Table 4Monthly average parameters for evaporation calculations for the 2015–2016 period.

Short-wave radiation was measured, while long-wave radiation was computed as described above. The average for all components is shown in Fig. 4. The radiation budget is positive every day, with a mean of about 150 Wm−2, varying from 100 in winter to 200 Wm−2 in summer.

Figure 4Monthly average radiative budget for Lake Titicaca in 2015–2016.

5 Results and discussion

## 5.1 Monthly data

Detailed energy balance computations over the period 2015–2016 should give good estimates of the total lake evaporation for that period. After 24 months the lake surface temperature at Puno more or less returned to the temperature at the beginning of 2015. When applying this method over the 2 years of study the mean annual lake evaporation is 1700 mm year−1. When computing the evaporation month by month, the change of heat storage was considered in the way previously described. The mixing depth was set to 40 m. The change of the heat storage is shown in Fig. 5. The values suggested by Carmouze et al. (1992) are shown for comparison. The calculated monthly heat storage agrees well with the Carmouze estimates.

Figure 5Change of heat storage in 2015–2016.

The computed monthly evaporation using monthly average data and the energy balance method was somewhat higher in 2016 than in 2015, 1725 mm year−1 as compared to 1680 mm year−1. The small gap of evaporation between 2 consecutive years is mainly explained by the warmer season that occurred in autumn of 2015; otherwise the evaporation was fairly evenly distributed over the year, being about 140 mm month−1, with somewhat lower evaporation rates from July to September (see in Fig. 6).

Figure 6Monthly evaporation computed using energy balance approach.

From the energy balance and the water balance methods, the annual evaporation from Lake Titicaca was estimated in the range of about 1700 mm year−1. The monthly variation depends on the change of heat storage, therefore the calculated evaporation may be high one month and low the following month. When using the mass transfer approach, similar annual evaporation to that from the energy balance approach may be anticipated when applying the approaches over 2 full years. This may be the case even though there may be differences when comparing monthly calculations. However, when the coefficients suggested by Carmouze et al. (1992) were used, the evaporation was much higher than 1720 mm year−1, which was found from the energy balance method. A good fit for the total evaporation was found using the coefficients a=0.17 mbar and b=0.155 mm mbar−1 s m−1.

Figure 7Comparison of monthly evaporation computed by energy balance and mass transfer method for 2015.

The monthly evaporation computed by mass transfer over the 2 years is compared with the energy balance calculations in Fig. 7 for 2015 and in Fig. 8 for 2016. The computed annual evaporation by the last method was 1700 mm year−1 in 2015 and 1675 mm year−1 in 2016. Consequently, for individual years the two methods gave similar results and similar seasonal trends. This is expected, since the coefficients in the mass transfer equation were chosen to fit over the 2-year period. For individual months, there are deviations. However, it is not possible to note any systematic differences related to different seasons of the year. The largest difference between the two methods for an individual month was about 30 mm.

Figure 8Comparison of monthly evaporation computed by energy balance and mass transfer method for 2016.

A summary and comparison of all investigated methods for the study period are shown in Fig. 9 and Table 5. As seen from these, the evaporation calculated from water balances differs from the three other methods. The water balance method yielded 1672 mm year−1 as the average for the tow year. As well, the same method gave the largest standard deviation. The variation in mean annual evaporation was 1633–1711 mm year−1. Since the lake water level can be observed at best with a resolution in centimetres, individual monthly evaporation becomes highly uncertain, like in February 2016 and May 2016. Thus, calculated annual evaporation is better performed using the three other methods. In general, the mass transfer, energy balance and Penman methods gave a similar monthly variation as described above; only the evaporation by Penman gave a smaller rate than the rest of models, 1621 mm year−1, while energy balance gave the highest evaporation rate, 1701 mm year−1.

Figure 9Monthly actual evaporation calculated by the four methods for the period from January 2015 to December 2016.

Table 5Descriptive statistics of monthly evaporation calculated by the four methods for the period from January 2015 to December 2016.

Water balances were computed for the long-term period and continuous available data in 1966–2011 as well. During the computation the Alake parameter in the model was assumed to be constantly equal to 8800 km2. The resulting annual lake evaporation for the period 1966–2011 was about 1600 mm year−1. The mean water balance components for the last period at monthly mean values are shown in Table 6. When we used the water balance approach the computed evaporation internally varied greatly from month to month and year to year. This is an indication that some hydrological input data are uncertain, like the precipitation, which can be improved from one point observation at lake shore by the satellite observation. In any case, water balance computations over a long-term period should give a reasonable estimate of mean lake evaporation.

Table 6Monthly mean water balances for Lake Titicaca for the period 1966–2011.

## 5.2 Using daily data

When calculating evaporation using daily data, it was found that there are large differences between the methods and ignoring the water balance method. The maximum daily evaporation using the mass transfer method was 12 mm day−1. Neither the energy balance nor the Penman method gave a higher evaporation than 8 mm day−1. There was poor agreement between the computed daily evaporation computed by mass transfer and the corresponding results using the other two methods.

When the mass transfer approach is used, it is straight-forward to determine the daily evaporation. Using the energy balance, the change of heat storage in the lake must be determined with high resolution. Detailed water temperature measurements are not available, however. Instead, it was assumed that the temperature changes at a steady rate through individual months. August is an example, since the temperature for the whole month changed very little (0.2 C). The computed daily evaporation is shown for August 2015 in Fig. 10. There were 2 days with average winds exceeding 6 ms−1. Consequently, the evaporation was high during these days, when the mass transfer approach was used.

Figure 10Daily evaporation computed by the mass transfer (shaded staples) and energy balance (filled staples) method.

When annual evaporation was determined using daily data instead of monthly mean data, there was hardly any difference for the mass transfer method. As indicated above it is not possible to use the energy balance method with short time resolution when temperature changes have to be taken in to account from day to day. However, the evaporation can be computed while neglecting the heat change, keeping the Bowen ratio constant throughout a month and changing the Bowen ratio day by day. In this case, it was found that evaporation increased by about 2 %. The conclusion, considering the many uncertainties involved in estimating evaporation, is that it is sufficient to use monthly means when estimating evaporation.

The evaporation computed with the Penman equation falls between what was found by the energy balance and the mass transfer approach, being somewhat closer to the energy balance than to the mass transfer results. Since monthly means are sufficient for computing evaporation with the two above methods, mean values are also sufficient when using the Penman method.

It has to be noted that Lake Titicaca's near-bank surface temperatures have been observed to be warmer than the lake surface's average during daytime using satellite thermal imagery, as reported in other lakes (e.g. Marti-Cardona et al., 2008). According to this observation, the temperatures acquired for this study are likely to be an overestimation.

The spatial distribution of Lake Titicaca's surface temperature and its impact on the evaporation losses is currently under analysis. However, for the energy balance method, daily changes rather than absolute temperatures were used, which are considered to be reasonable approximations of the heat storage changes. Over the larger period of air temperatures observed at the Copacabana weather station (1966–2016), the particular months in 2015–2016 have been characterized by the strongest El Niño dry phenomena during the last 50 years (http://www.ciifen.org, last access: 10 January 2017); in comparison to the rest of the years, the air temperatures recorded were higher than the average of 10 C and close to 20 C at the daily time step. Then the rates of evaporation found might express, somehow, the indicated warmer period.

6 Conclusions

Due to uncertainty of most observed data such as river inflow to Lake Titicaca, and mainly the discharge data, it might not be easy to improve the water balance results; thus it is suggested that the most reliable method of determining the lake evaporation is using the heat balance approach. To estimate the lake evaporation using this method, heat storage changes must be known. Since convection from the surface layer is intense during nights, resulting in a well-mixed top layer every morning, it is possible to determine the change of heat storage from the measured morning surface temperature. The lake evaporation is fairly uniformly distributed over the year, with lows between July and September. The mean annual evaporation is about 1700 mm year−1, and the mean monthly evaporation is 141.8 mm month−1. When using the mass transfer approach, the required coefficients in the aerodynamic equation was set so that the mean annual evaporation agreed with that obtained from the heat balance calculations. These coefficients were found to be lower than coefficients used in previous studies. Also, when using the mass transfer approach, the evaporation was found to be lowest in July–September.

However, for the purpose of assessing climate change effects on Lake Titicaca's evaporation, the practical approach, rather than the two empirical models, might be the Penman equation due to available observed data for this lake and the integral behaviour of the equation. Also in comparison with the two models proposed in Delclaux et al. (2007) for modelling the lake evaporation, the first model only depends on the solar radiation data, and, additionally, the second one depends on the air temperature factor; thus both models cannot be applied broadly. In the Penman model based on the adjusted wind coefficient, the mean annual evaporation is 1620 mm year−1, and the mean monthly is 135 mm month−1. So far, monthly evaporation computed using daily data and monthly means resulted in minor differences. The most practical model for use at the daily scale might be the mass transfer and the Penman models in comparison to the energy balance approach, which is of highly demanded observed data. Particularly the Penman equation at the daily temporal scale might correctly be applied for the climate change assessment at this altitude. Nonetheless, according to spatial available data from remote sensing, the evaporation equations used at daily and monthly scales could be applied from now on to improve the spatial pattern of the lake evaporation. Since we had really extreme single warmer days during the period 2015–2016 due to the El Niño phenomenon, higher daily rates of evaporation must be expected; therefore the application of the models at both timescales for the study period we believe that was found the upper limits of yearly evaporation.

Data availability
Data availability.

All data can now be freely accessed through requests to rpillco@umsa.edu.bo.

Author contributions
Author contributions.

RPZ coordinated the research and was directly involved in all steps, from fieldwork to proofreading. He built the database, computed the evaporation for the different models, analyzed the results, and prepared figures and tables. LB contributed to the conceptual approach and structure of the paper. He supervised and validated the calculations and contributed to the writing of the objectives and the scientific background of the paper. RB revised and validated the calculations and collaborated on the writing, mainly for Sect. 5. BMC contributed to the discussion and analysis of results and writing, especially for the Abstract and Sects. 1 and 6. FS collaborated on the analysis of results, particularly on the interpretation of some evaporative models. He prepared the figures depicting maps. FT was in charge of the installation and maintenance of the gauging stations and data quality assurance. MPB improved the conceptual approach of the paper. She also helped to obtain funding for the field data acquisition. LM assisted in the paper drafting and building the database base. He also prepared the chart figures. CG and JP facilitated meteorological records from Peruvian gauging stations. They contributed to the paper structure and to the content of Sects. 1 and 5.

Competing interests
Competing interests.

The authors declare that they have no conflict of interest.

Special issue statement
Special issue statement.

Acknowledgements
Acknowledgements.

We would like to express our sincere appreciation to the HASM, Research Programme – Hydrology of Altiplano from Space to Modeling at GET-IRD and IHH-UMSA (Instituto de Hidráulica e Hidrología, UMSA, Bolivia), financed by the TOSCA-CNES (Centre National d'Etudes Spatiales). We would like to thank SENAMHI-Bolivia (Servicio Nacional de Hidrometeorología de Bolivia) for providing long-term climatic data. Thanks also to the IMARPE-Perú (Instituto para el Mar del Perú/Puno) for providing additional hydrological data as well as surface water temperatures of Lake Titicaca. In addition, our acknowledgment is directed to the project Fortalecimiento de Planes Locales de Intervención y Adaptación al Cambio Climático en el Altiplano Boliviano at Agua Sustentable-Bolivia for providing the Lake Titicaca discharge data. Finally, we thank the programme BABEL Erasmus EU for providing economic assistance and completing this work in Sweden.

Reviewed by: two anonymous referees

References

Abarca-Del-Río, R., Crétaux, J. F., Berge-Nguyen, M., and Maisogrande, P.: Does Lake Titicaca still control the Lake Poopó system water levels? An investigation using satellite altimetry and MODIS data (2000–2009), Remote Sens. Lett., 3, 707–714, https://doi.org/10.1080/01431161.2012.667884, 2012.

Bengtsson, L., Herschy, R. W., and Fairbridge, R. W.: Encyclopedia of Lakes and Reservoirs, Springer Dordrecht, ISBN-13: 978-1-4020-5616-1, 2012.

Boulange, B. and Aquize, J. E.: Morphologie, hydrographie et climatologie du lac Titicaca et de son basin versant, Rev. Hydrobiol. Trop., 14, 269–287, 1981.

Bowen, I. S.: The ratio of the heat losses by conduction and by evaporation from any water surface, Phys. Rev., 27, 779–787, 1926.

Canedo, C., Pillco, Z. R., and Berndtsson, R.: Role of hydrological studies for the development of the TDPS system, Water, 8, 1–14, https://doi.org/10.3390/w8040144, 2016.

Carmouze, J. P.: The energy balance in the Lake Titicaca, in: A Synthesis of Limnological Knowledge, edited by: Dejoux, C. and Iltis, A., Kluwer, 131–146, 1992.

Carmouze, J. P., Arce, C., and Quintanilla, J.: La régulation hydrique des lacs Titicaca et Poopo, Cahiers ORSTOM, Série Hydrobiologique, XI, 269–283, 1977.

CMLT: Comisión multisectorial para la preservación y recuperación ambiental del Lago Titicaca y sus afluentes, Estado de la calidad ambiental de la cuenca del Lago Titicaca ámbito peruano, CMLT, 152 pp., 2014.

Crawford, T. and Duchon, C. E.: An improved parametrization for estimating effective atmospheric emissivity for use in calculating daytime downwelling long-wave radiation, J. Appl. Meteorol., 38, 474–480, https://doi.org/10.1175/1520-0450(1999)038<0474:AIPFEE>2.0.CO;2, 1999.

Dalton, J.: Experimental essays on the constitution of mixed gases, Manchester Literary and Philosophical Society Memo, 5, 535–602, 1802.

Dejoux, C. and Iltis, A. (Eds.): Lake Titicaca, A synthesis of Limnological Knowledge, Springer Science, Kluwer Academic Publisher, 1992.

Delclaux, F., Coudrain, A., and Condom, T.: Evaporation estimation on Lake Titicaca: a synthesis review and modeling, Hydrol. Process., 21, 1664–1677, https://doi.org/10.1002/hyp.6360, 2007.

FAO: The state of the world's land and water resources for the food and agriculture (SOLAW) – Managing system at risk, Food and organization of the United Nation, Rome and Earthscan, London, ISBN-13: 978-92-5-106614-0, 2011.

Garreaud, R.: Multiscale Analysis of the Summertime Precipitation over the Central Andes, Mon. Weather Rev., 127, 901–921, https://doi.org/10.1175/1520-0493(1999)127<0901:MAOTSP>2.0.CO;2, 1999.

Garreaud, R., Vuille, M., and Clement, A. C.: The climate of the Altiplano: observed current conditions and mechanism of past changes, PALAEO, 194, 5–22, https://doi.org/10.1016/S0031-0182(03)00269-4, 2003.

Gutiérrez, C. B. and Molina, C. J.: Final report: Balance hídrico y escenarios de cambio climático en el lago Titicaca, Project: Fortalecimiento de Planes Locales de inversión y adaptación al cambio climático en el Altiplano boliviano, IDRC-CRDI, Agua Sustentable, La Paz, Bolivia, https://doi.org/10.13140/RG.2.2.24230.22083, 2014.

Lhomme, J. P., Vacher, J. J., and Rocheteau, A.: Estimating downward long-wave radiation on the Andean Altiplano, Agr. Forest Meteorol., 145, 139–148, https://doi.org/10.1016/j.agrformet.2007.04.007, 2007.

López-Moreno, J. I., Morán-Tejada, S. M., Vicente-Serrano, J., Bazo, C., Azorin-Molina, J. I., Revuelto, A., Sánchez-Lorenzo, F., Navarro-Serrano, E., Aguilar, E., and Chura, O.: Recent temperature variability and changes in the Altiplano of the Bolivia and Perú, Int. J. Climatol, 36, 1773–1796, 2015.

Maidment, D. R.: Handbook of Hydrology, Mc-GRAW-HILL, INC., ISBN-13: 978-0-07-039732-3, 1993.

Marti-Cardona, B., Streissberg, T. E., Schladow, S. G., and Hook, S. J.: Relating fish kills to upwellings and wind patterns at the Salton Sea, Hydrobiologia, 604, 85–95, 2008.

Penman, H. L.: Natural Evaporation from open water, Bare Soil and Grass, P. Roy. Soc. A-Math. Phy., A193, 120–145, https://doi.org/10.1098/rspa.1948.0037, 1948.

Pillco, Z. R. and Bengtsson, L.: Long-term and extreme water levels variations of the shallow Lake Poopó, Bolivia, Hydrolog. Sci. J., 51, 98–114, https://doi.org/10.1623/hysj.51.1.98, 2006.

Pillco, Z. R. and Bengtsson, L.: Doctoral theses: Response of Bolivian Altiplano Lakes to Seasonal and Annual Climate Variation, Lund University Publications, 1, 107 pp., 2007.

Pouyaud, B.: Lago Titicaca, Programa estudios TDPS: misión d'evaluation sur la prise en compte de l'evaporation, GIE Hydro Consult International, GIE ORSTOM-EDF, 1993.

Revollo, M. M.: Management issues in the Lake Titicaca and Lake Poopó System: Importance of development a water budget, Lakes and Reservoirs: Research and Management, 6, 225–229, https://doi.org/10.1046/j.1440-1770.2001.00151.x, 2001.

Richerson, P. J., Widmer, C., and Kittel, T.: The Limnology of the Lake Titicaca (Perú-Bolivia), a Large, High Altitude Tropical Lake, Institute of Ecology Publication, 14, 82 pp., 1977.

Roche, M. A., Bouges, J., and Mattos, R.: Climatology and hydrology of the Lake Titicaca basin, in: Lake Titicaca a Synthesis of Limnological Knowledge, edited by: Dejoux, C. and Iltis, A., Kluwer Academic Publisher, 63–88, 1992.

Satgé, F., Espinoza, R., Pillco, Z. R., Roig, H., Timouk, F., Molina, J., Garnier, J., Calmant, S., Seyler, F., and Bonnet, M.-P.: Role of climate variability and human activity on Poopó Lake droughts between 1990 and 2015 assessed using Remote Sensing Data, Remote Sens., 9, 1–17, https://doi.org/10.3390/rs9030218, 2017.

Singh, V. P. and Xu, C.-Y.: Sensitivity of mass transfer-based evaporation equations to errors in daily and monthly input data, Hydrol. Process., 11, 1465–1473, 1997.

Talbi, A., Coudrain, A., Ribstein, P., and Pouyaud, B.: Calcul de la pluie sur le bassin-versant du lac Titicaca pendant l'Holocène, Comptes Rendus de l'Accadémie des Sciences, Paris, Sciences de la Terre et des Planètes, 329, 197–203, 1999.

Taylor, M. and Aquize, E. A.: Climatological energy Budget of the Lake Titicaca (Peru/Bolivia), Internationale Vereinigung für Theoretische und Angewandte Limnologie: Verhandlungen, 22, 1246–1251, https://doi.org/10.1080/03680770.1983.11897478, 1984.

TDPS: Climatología del sistema de los lagos Titicaca, Desaguadero, Poopó y Salares Coipasa y Uyuni (TDPS), Comisión de comunidades de europeas, Repúblicas del Perú y Bolivia, convenios ALA/86/03 y ALA/87/23, LP, Bolivia, 1993.

UNEP.: Diagnóstico Ambiental del Sistema Titicaca-Desaguadero-Poopó-Salar de Coipasa (TDPS) Bolivia-Perú, División de Aguas Continentales del Programa de las Naciones Unidas para el Medio Ambiente y Autoridad Binacional del Sistema TDPS, Washington, DC, USA, p. 192, 1996.

Vuille, M., Radley, R. S., and Keiming, F.: Interannual climate variability in the Central Andes and its relation to tropical Pacific and Atlantic forcing, J. Geophys. Res.-Atmos., 105, 12447–12460, https://doi.org/10.1029/2000JD900134, 2000.

Wirrmann, D.: Morphology and bathymetry, in: Lake Titicaca, In a synthesis of Limnological Knowledge, edited by: Dejoux, C. and Iltis, A., Monographiae Biologicae, 68, Kluwer Academic Publisher, Dordrecht, 16–22, 1992.