Mountain Lake Evaporation: A Comparative Study between Hourly Estimations Models and In Situ Measurements

: The quantiﬁcation of evaporation is very important for water resource management, determination of water availability and environmental modeling. Evaporation ( E ) is a key component of the hydrological cycle, which becomes more relevant under a changing climate scenario, where increases in temperature and E are projected on a regional scale. This study evaluated the hourly applicability of 21 methods to estimate E in a mountain lake. Models were grouped as combination, solar radiation-temperature and mass transfer methods. Estimated E by all models were compared to three measurement campaigns with an eddy covariance system in Lake Laja during the month of January (22 to 28), March (17 to 23) and November (09 to 15) of 2016. The average evaporated water during the measurement sessions were 3.40, 3.38 and 1.89 mm d − 1 , respectively. Best model performance was obtained with models whose main E principles are heat ﬂow and water vapor ﬂow. The best performance in this group was the Penman model using a calibrated wind function with a determination coe ﬃ cient ( R 2 ) of 0.91, Nash–Sutcli ﬀ e coe ﬃ cient (NS) of 0.9, and index of agreement (W) of 0.98. Comparing daytime and nighttime hours, it was found that the daytime E in the lake is best explained by the product of the vapor pressure deﬁcit and wind speed, while the nighttime E was explained by the ﬂow of heat in the water. The results highlight the importance of the analysis of diurnal dynamics of water ﬂux and energy stored in water to better understand the E in water bodies.


Introduction
The effects of climate change and rising temperatures throughout the world, and particularly in arid and semi-arid regions, are threatening to reduce available surface water through enhanced surface evaporation, especially in reservoirs [1]. Mountains and agricultural regions located in snowmelt-dominated Mediterranean climate watersheds have been identified as highly vulnerable to climate change. Temperature levels in mountain regions are expected to increase by 3-4 • C [2,3], thus increasing E rates in mountain lakes.
Evaporation is an important source of water loss in lakes and reservoirs in arid and semi-arid regions. The Mediterranean climate or dry summer climate exists in the Mediterranean basin, some parts of South Africa, California, Australia, and Chile [4]. Loss of water by evaporation can body but regulated in a controlled way by the El Toro Hydroelectric Plant. It has a volume of approximately 5600 Hm 3 , which is used for agricultural purposes and as the cornerstone of hydroelectric supply from the interconnected power system in Chile.
The measurement site was located on the southern shore of the lake (Figure 1), 4.5 km from the Antuco ski center (37°21′43.86″ S, 71°20′24.60″ O a 1360 m.a.s.l). To the south of the measurement site is the Antuco volcano, and to the north is a hilly area with little vegetation. The soils in the area adjacent to the lake are sandy and stony, not suitable for agriculture, livestock or forestry, and with limited use for wildlife, recreation or the protection of river basins. The site location was selected considering previous information on the direction of the prevailing winds in the lake. This previous analysis was considered to ensure that the fluxes of the surface energy balance measured by the EC system comes mainly from the body of water. Wind speed in the measurement time period reached 8 m s −1 , predominantly blowing in a direction from west to east.

Equipment and Measurements
An EC system was installed in the study area to measure the main components of the energy budget and the main environmental variables that affect E at 30 min intervals in January (from 22 to 28), March (from 17 to 23) and November (from 9 to 15) of 2016. Specific days during field campaigns were selected considering the weather forecast to find high probabilities of clear and sunny days during measurements.

Equipment and Measurements
An EC system was installed in the study area to measure the main components of the energy budget and the main environmental variables that affect E at 30 min intervals in January (from 22 to 28), March (from 17 to 23) and November (from 9 to 15) of 2016. Specific days during field campaigns were selected considering the weather forecast to find high probabilities of clear and sunny days during measurements. 1.5, 3 and 4.5 m depth, all attached to a marine rope with a weight and affixed with a buoy. Wind speed was measured with the same 3D sonic anemometer. Signals from these instruments were sampled every 5 s and stored as the mean value after 30 min in a data logger (CR3000, Campbell Scientic, Logan, UT, USA). Two 12 V batteries connected in parallel supplied the power.
To determine water heat flux (Q x ), the reduced form of the energy budget was used following [27,31,33,40,41]: where, Q x , Q * , λE and H s is in W m −2 .
In order to ensure that latent heat and sensible heat fluxes came mainly from the water body, a footprint analysis was performed using the method developed by Hsieh et al. [42].

Evaporation Estimation Methods
Twenty-one models were used to estimate evaporation on an hourly basis, as indicated in the following groups: Combination equations: The combined method generally includes energy transfer and aerodynamic components. For practical applications, the formulation requires measurement of radiation fluxes on the water body, water surface temperature, and meteorological data such as wind speed, air temperature, and dew point temperature at a given height. Solar radiation and temperature: This group consists of four models, and has the advantage of requiring fewer physical parameters, such as solar radiation, air temperature, atmospheric pressure, and relative humidity. However, as in all the other methods, measurements should be representative of the conditions of the water body.
Mass transfer: This group consists of three models. Mass transfer is based on physical principles. Evaporation calculated from this method is considered proportional to a function of the mean wind speed, and the difference between air vapor pressure at the evaporation surface and at some level above the surface of the water body. Table 1 shows the equations involved for each of the 21 models. The last 7 models of the mass transfer group were proposed by Singh and Xu [29], and require the estimation of the parameters a, b and c. The parameters were estimated using a method based on the least squares criterion between measured and estimated evaporation. For empirical models, the values of the parameters need to be calibrated from a sample of input and output values of the model, and from an objective function whose value must be minimized. One of the simplest objectives is a linear regression model that uses an input variable x, and an output variable y. The model to be calibrated is an equation of type y' = A + BX in which parameters A and B can minimize the objective function n i=1 (y − y) 2 , where n is the number of data pairs.

Mass transfer group
Wind speed at 2 m above water surface (m s −1 ) (km hr −1 for Trivett). The Penman, Brutsaert-Stricker and deBruin models consider the following wind function: where a and b are adjustment constants, which will be calibrated to maximize the function of the coefficient of determination (R 2 ).

Results
The components of the water surface energy budget and the environmental variables that affect evaporation were measured for seven days in January (from 22 to 28), March (from 17 to 23) and November (from 9 to 15) of 2016. The data were filtered, eliminating the values that were attributable to noise in the data, sensor failure or rainy periods, resulting in a total of 747 half hour data (259, 255 and 233 in January, March and November, respectively). Figure 2 shows the environmental variables measured. Table 3 (Figure 2e) ranged between 10% and 90%, reaching its maximum and minimum values in morning and afternoon hours, respectively, showing a behavior contrary to the vapor pressure deficit (VPD) curve. Maximum values for relative humidity were recorded in January. The maximum temperature difference between water surface (T 0 ) and air temperature (Figure 2f) occurred between six and eight in the morning, reaching up to 14 • C. The daily variation of water surface temperature (Figure 2d) did not exceed ±2 • C, reaching its maximum during the afternoon hours. In November, there was a failure in the sensor that recorded water surface temperature, so there are no data for the 315 day of the year (DOY) and half of the day 316.

Environmental Components
Water 2020, 12, x FOR PEER REVIEW 8 of 25 did not exceed ±2 °C, reaching its maximum during the afternoon hours. In November, there was a failure in the sensor that recorded water surface temperature, so there are no data for the 315 day of the year (DOY) and half of the day 316. is wind speed at 2 m above water surface (a), is the vapor pressure deficit (c), is the relative humidity (e), is the air temperature (b), is water surface temperature (d) and, ∆T is the difference between water lake temperature and air temperature (f).  . Measured environmental variables. U 2 is wind speed at 2 m above water surface (a), VPD is the vapor pressure deficit (c), RH is the relative humidity (e), T a is the air temperature (b), T 0 is water surface temperature (d) and, ∆T is the difference between water lake temperature and air temperature (f).  Figure 3 shows the measured solar radiation (Q s ). The maximum values were recorded between 2 and 3 p.m., reaching 1200 W m −2 , 900 W m −2 and 1050 W m −2 in January, March and November, respectively.
Water 2020, 12, x FOR PEER REVIEW 9 of 25 Figure 3 shows the measured solar radiation ( ). The maximum values were recorded between 2 and 3 p.m., reaching 1200 W m −2 , 900 W m −2 and 1050 W m −2 in January, March and November, respectively.

Footprint
The footprint analysis shows that 90% of the fluxes, under unstable conditions, came from an area 500 m away from the station, while the most important contributions came from less than 100 m of distance. The fluxes in the study area predominantly come from the northwest, this assures that all the surface energy balance measurements mainly came from the surface of interest. Figure 4 shows the flux direction under unstable conditions (north west direction), distance (Xpeak) to the point of maximum influence from the measurement point (generally less than 100 m) and the flux footprint distribution function.

Footprint
The footprint analysis shows that 90% of the fluxes, under unstable conditions, came from an area 500 m away from the station, while the most important contributions came from less than 100 m of distance. The fluxes in the study area predominantly come from the northwest, this assures that all the surface energy balance measurements mainly came from the surface of interest. Figure 4 shows the flux direction under unstable conditions (north west direction), distance (Xpeak) to the point of maximum influence from the measurement point (generally less than 100 m) and the flux footprint distribution function.   Figure 5 shows the daily curves for the surface energy budget components. Table 4 shows the daily maximum and mean values of all components. The net radiation reached a maximum of 862 W m −2 at 2:30 p.m. in January, recording negative values at nighttime. Furthermore, 73% of the observations of the sensible heat flux were positive, indicating that there is heat transfer from the lake  Figure 5 shows the daily curves for the surface energy budget components. Table 4 shows the daily maximum and mean values of all components. The net radiation reached a maximum of 862 W m −2 at 2:30 p.m. in January, recording negative values at nighttime. Furthermore, 73% of the observations of the sensible heat flux were positive, indicating that there is heat transfer from the lake to the atmosphere. Negative values of sensible heat flux were observed during the afternoon and evening hours, indicating that there is heat transfer from the atmosphere to the lake. Water heat fluxes were calculated according to Equation (1) 22  184  91  722  207  37  2  591  113  23  217  117  862  181  56  2  674  62  24  183  104  221  6  92  57  99  −155  25  159  102  534  78  55  40  345  −64  26  178  120  669  139  60  42  472  −23  27  189  113  508  59  87  46  315  −100  78  209  121  608  100  128  38  324  −58  79  182  114  593  109  97  14  479  −20  80  205  99  602  114  41  3  556  12  81  202  85  599  113  130  24  473  4  82  191  71  580  108  145  10  549  28  315  104  45  665  190  157  12  559  133  316  129  56  686  190  54  −2  587  136  317  117  61  704  202  42  6  576  134  318  103  63  606  151  62  22  458  67 The following variables were correlated with the measured latent heat flux to identify those that most significantly affect E: the relationships between latent heat flux (measured every 30 min) and wind speed above the lake surface (Figure 6a  shows that the lowest relationship is given by  and , followed by  and * . It can be observed that λE is more influenced by the product of and 2 (R 2 = 0.43).  The following variables were correlated with the measured latent heat flux to identify those that most significantly affect E: the relationships between latent heat flux (measured every 30 min) and wind speed above the lake surface (Figure 6a Figure 6 shows that the lowest relationship is given by λE and Q x , followed by λE and Q * . It can be observed that λE is more influenced by the product of VPD and U 2 (R 2 = 0.43). During these campaigns, the average evaporated water depths reached 3.40, 3.38 and 1.89 mm per day in January, March and November, respectively. Table 5 shows the evaporated water depths for the days measured, and its equivalence as E flow (m 3 s −1 ), which was calculated based on the area flooded by the lake (approximately 70 km 2 ). The lowest E rate was observed in November, while the highest value was recorded in January.  Figure 7 shows the volume of daily water evaporated according to the elevation of the water level of the lake for different depths of evaporated water. This indicates that the higher the elevation During these campaigns, the average evaporated water depths reached 3.40, 3.38 and 1.89 mm per day in January, March and November, respectively. Table 5 shows the evaporated water depths for the days measured, and its equivalence as E flow (m 3 s −1 ), which was calculated based on the area flooded by the lake (approximately 70 km 2 ). The lowest E rate was observed in November, while the highest value was recorded in January.  Figure 7 shows the volume of daily water evaporated according to the elevation of the water level of the lake for different depths of evaporated water. This indicates that the higher the elevation of the lake, the greater the volume of water lost by evaporation. An important change occurs between elevation 1360 and 1370 m because one of its branches is active, which implies a significant jump in the curves of evaporated water volume.

DOY
Water 2020, 12, x FOR PEER REVIEW 13 of 25 of the lake, the greater the volume of water lost by evaporation. An important change occurs between elevation 1360 and 1370 m because one of its branches is active, which implies a significant jump in the curves of evaporated water volume.
On the other hand, if the lake were at its maximum level (1390 m), at an evaporation rate of 4 mm d −1 , the volume of water evaporated could reach 480,000 m 3 d −1 .  As observed in Figure 8, the Penman and Penman-Monteith models show similar relationships (R 2 ), indicating the similarity between the models. The same occurs with the Priestley-Taylor and deBrauin-Keijman models. The best relationship between measured and estimated evaporation in this group was obtained with the Priestley-Taylor model. On the other hand, if the lake were at its maximum level (1390 m), at an evaporation rate of 4 mm d −1 , the volume of water evaporated could reach 480,000 m 3 d −1 .     The wind functions of the Penman, deBruin, Brutsaert-Stricker models were adjusted according to equation (2). The adjustment was performed to maximize the function of the coefficient of determination (R 2 ) as:

Model Comparisons with Measurements
With the adjustment, R 2 changed from 0.84 to 0.91 in Penman, 0.43 to 0.44 in deBruin and 0.14 to 0.84 in Brutsaert-Stricker.
As for the solar radiation and temperature group, the data were filtered only considering the daytime ranges where solar radiation exceeded 100 W m −2 .
All models in this group overestimate E ( Figure 10). The efficiency indicator (NS) for all models in this group is negative, indicating that the mean of the observed data has a higher predictive power than the model. A similarity of measured and estimated E relationships is observed between the Stephens-Stewart and Jensen-Haise models. This is explained by the similarity in their equations. Based on this, the best model of this group is the Stephens-Stewart model because it presents the lowest overestimation. As observed in Figure 8, the Penman and Penman-Monteith models show similar relationships (R 2 ), indicating the similarity between the models. The same occurs with the Priestley-Taylor and deBrauin-Keijman models. The best relationship between measured and estimated evaporation in this group was obtained with the Priestley-Taylor model.
The wind functions of the Penman, deBruin, Brutsaert-Stricker models were adjusted according to equation (2). The adjustment was performed to maximize the function of the coefficient of determination (R 2 ) as: With the adjustment, R 2 changed from 0.84 to 0.91 in Penman, 0.43 to 0.44 in deBruin and 0.14 to 0.84 in Brutsaert-Stricker.
As for the solar radiation and temperature group, the data were filtered only considering the daytime ranges where solar radiation exceeded 100 W m −2 .
All models in this group overestimate E (Figure 10). The efficiency indicator (NS) for all models in this group is negative, indicating that the mean of the observed data has a higher predictive power than the model. A similarity of measured and estimated E relationships is observed between the Stephens-Stewart and Jensen-Haise models. This is explained by the similarity in their equations. Based on this, the best model of this group is the Stephens-Stewart model because it presents the lowest overestimation.
Similarly, the mass transfer models present an acceptable linear relationship, same as most of their indicators, except for the NS indicator. As mentioned before, water surface temperature data of two days were not recorded and, therefore, these models were evaluated with less data (n = 659). Figure 11 shows the relationship between estimated values with three models of the mass transfer group and measured values. The best linear relationship is given by the Ryan-Harleman model. However, this does not mean that it is the best model, since it presents an overestimation of E values (72%). Of these three models, the best is Trivett's as it has the best indicators and the least overestimation. The parameters a, b and c of the generalized form of the mass transfer model developed by Sing and Xu [29], were adjusted through a least-squares fitting Table 6. The best adjustments are given by the models B, D, E and G ( Figure 12).
According to the analyzed indicators, the G model is the best model of the generalized form of the mass transfer models. It is important to note that the D model requires the estimation of one parameter less, which is an advantage when processing the data. In addition, the difference obtained in the NS index was not significant.  The parameters a, b and c of the generalized form of the mass transfer model developed by Sing and Xu [29], were adjusted through a least-squares fitting Table 6. The best adjustments are given by the models B, D, E and G ( Figure 12).  As many of the equations for estimating E have constants or coefficients that were determined by those who developed or modified them to meet their particular needs, adjustment factors ( ) for each of the hourly models were obtained with the least squares model. The models to estimate E are expressed as follows: = ( ) * (6) where ( ) is the E function of each of the models evaluated. Table 7 shows the statistical indicators evaluated for conditions prior to and after the application of KLL. For both conditions, it is observed that the best model is Penman including the wind function ( ) (Penman WF). For the models of the group of solar radiation and temperature, the best model was the Stephens-Stewart model even after applying . Finally, the best model of the mass transfer models was Trivett, both prior to and after the application of .  According to the analyzed indicators, the G model is the best model of the generalized form of the mass transfer models. It is important to note that the D model requires the estimation of one parameter less, which is an advantage when processing the data. In addition, the difference obtained in the NS index was not significant.
As many of the equations for estimating E have constants or coefficients that were determined by those who developed or modified them to meet their particular needs, adjustment factors (K LL ) for each of the hourly models were obtained with the least squares model. The models to estimate E are expressed as follows: where f (E) is the E function of each of the models evaluated. Table 7 shows the statistical indicators evaluated for conditions prior to and after the application of K LL . For both conditions, it is observed that the best model is Penman including the wind function f (U) (Penman WF). For the models of the group of solar radiation and temperature, the best model was the Stephens-Stewart model even after applying K LL . Finally, the best model of the mass transfer models was Trivett, both prior to and after the application of K LL . Table 7. Statistical indicators prior to and after the application of adjustment parameters for the Laja Lake (K LL ) (* indicates that the displayed factor is the one obtained after applying K LL ).

Measurements
In the southern hemisphere, January and March correspond to summer months. January is the month with the highest evaporative demand in the year, while November corresponds to a spring month. Days 22 and 23 were clear days, and days 24 to 27 were cloudy, which caused a drop in the air temperature recorded on those days compared to previous days ( Figure 2). The variation in water surface temperature between a cloudy and a clear day for the period between days 22 and 27 reached an average of 0.8 • C. Between days 24 and 27, it was observed that the water surface temperature was higher than the air temperature, indicating that there was heat transfer from the lake surface to the near atmosphere. The decrease in air temperature caused a decrease in e a , which reduced the capacity of the atmosphere to store water vapor, as reflected in the decrease in E around those days. The water surface temperature was lower in November, considering snow-rain driven regime. The decrease in water temperature in November can be associated with the beginning of thaws.
Air temperature dropped sharply during nighttime due to the absence of solar radiation. Water temperature had an average daytime variation of 1.3, 1.6 and 2.9 • C in January, March and November, respectively. Therefore, water temperature was higher than that of air during the night, causing heat transfer from the lake water to the atmosphere. For the night period, water heat flux was negative. There was a positive but not significant correlation between ∆T and Hs (R 2 = 0.22).
The latent heat flux for a water surface is expected to be high and even greater than that of a crop. However, the effect of water depth is considerable and influences the heat storage capacity of the water body. The maximum latent heat flux (217 W m −2 ) was recorded on a clear day in January at 3:00 p.m., a similar value was found in [10]. This peak coincides with the highest wind speed recorded during the measurement time period (9.1 m s −1 ) but occurred later than the maximum solar radiation (2:30 p.m.). Daytime measurements of latent heat flux had a similar pattern, reaching the highest values in the afternoon hours. The λE reaches average values of 108, 98 and 56 W m −2 during the months of January, March and November, respectively. Miranda et al. [11] measured E using an EC system in a Mediterranean climate reservoir (Alqueba, Portugal) during summer months, they found for the hottest month a λE value of 113 W m −2 similar to those found for the Laja Lake.
Nighttime E rates (between 8:00 p.m. and 7:00 a.m.) corresponded to an average of 36.9, 40.7 and 26.4% of the daytime evaporation rates recorded in January, March and November, respectively. This can be explained by the energy release from the water body. Similar rates were found by Liu et al. [41] at the Ross Barnett dam, Mississippi. This lake has an area of 133 km 2 and water depths between 4-8 m. Measurements with an EC system showed that nighttime evaporation represents 45% of total E. The heat flux in the water presented the lowest correlation when compared to latent heat. However, when analyzing the daytime and nighttime ranges separately (Figure 13), it can be observed that nighttime measurements showed a significant correlation (R 2 = 0.64) with λE, while variables were inversely proportional (between 8:00 a.m. and 7:00 p.m.), i.e., as water heat flux decreases, latent heat flux increases. As for the daytime range (between 7:30 a.m. and 7:30 p.m.), the relationship between latent heat flux and heat flux in water was not significant (R 2 = 0.01). On the other hand, when analyzing the other components for the daytime range, it is observed that E is better correlated during these hours of the day, due to the product between vapor pressure deficit and wind speed (R 2 = 0.51).
Water 2020, 12, x FOR PEER REVIEW 19 of 25 correlated during these hours of the day, due to the product between vapor pressure deficit and wind speed (R 2 = 0.51). Large and high-altitude lakes have a high capacity to store energy. In summer, most of the solar radiation is absorbed by the water body and the greatest heat storage capacity is recorded in late summer and early fall. Afterwards, the stored energy begins to release. In late fall and early winter, when solar radiation is small, this energy is released to the atmosphere as latent and sensible heat fluxes; the lowest stored heat is recorded in spring, when the heating process of the water body starts again [43]. This explains the decrease in the nighttime E rate recorded in November (26.4%).
In the Laja Lake, there is a low relationship between latent heat flux and net radiation. This is mainly due to the depth of the lake; short-wave radiation penetrates the lake surface where the energy is absorbed. Granger and Hedstrom [32] studied this relationship for a 25 m deep lake with hourly measurements and found an R 2 coefficient of 0.01, an even lower value than that in the present study.
The relationship between wind speed and latent heat flux was greater than the relationship between net radiation and evaporation fluxes. Same as in [32], it is observed that E in the Laja Lake is controlled by wind speed above the lake, similar observations were found in [10]. In addition, there is a directly proportional relationship between the latent heat flux and the product between vapor pressure deficit and wind speed. This occurs because evaporation rate depends directly on vapor pressure deficit (VPD) and is enhanced by higher wind speeds.
During measurement in January and March, evaporated water reached an average of 3.4 mm per day, equivalent to 266 million L. If we consider a daily water consumption of 200 L per person, the evaporated volume is equivalent to the consumption of 3644 people in a year. By expressing E in terms of flux, E flux averaged 2.76 m 3 s −1 for the measurement days in January and March, and 1.52 m 3 s −1 for those in November. Considering that the Laja Lake receives contributions from the Polcura River with flow rates that averaged 3.47 m 3 s −1 , 2.34 m 3 s −1 and 8,9 m 3 s −1 in January, March and November, respectively, it can be indicated that at least half of the water entering through the Polcura tunnel is lost by E.

Models
The evaluation of the combined group models showed that the best estimate was given by the Priestley-Taylor model. This model shows similarity in its equation with the deBrauin-Keijman Large and high-altitude lakes have a high capacity to store energy. In summer, most of the solar radiation is absorbed by the water body and the greatest heat storage capacity is recorded in late summer and early fall. Afterwards, the stored energy begins to release. In late fall and early winter, when solar radiation is small, this energy is released to the atmosphere as latent and sensible heat fluxes; the lowest stored heat is recorded in spring, when the heating process of the water body starts again [43]. This explains the decrease in the nighttime E rate recorded in November (26.4%).
In the Laja Lake, there is a low relationship between latent heat flux and net radiation. This is mainly due to the depth of the lake; short-wave radiation penetrates the lake surface where the energy is absorbed. Granger and Hedstrom [32] studied this relationship for a 25 m deep lake with hourly measurements and found an R 2 coefficient of 0.01, an even lower value than that in the present study.
The relationship between wind speed and latent heat flux was greater than the relationship between net radiation and evaporation fluxes. Same as in [32], it is observed that E in the Laja Lake is controlled by wind speed above the lake, similar observations were found in [10]. In addition, there is a directly proportional relationship between the latent heat flux and the product between vapor pressure deficit and wind speed. This occurs because evaporation rate depends directly on vapor pressure deficit (VPD) and is enhanced by higher wind speeds.
During measurement in January and March, evaporated water reached an average of 3.4 mm per day, equivalent to 266 million L. If we consider a daily water consumption of 200 L per person, the evaporated volume is equivalent to the consumption of 3644 people in a year. By expressing E in terms of flux, E flux averaged 2.76 m 3 s −1 for the measurement days in January and March, and 1.52 m 3 s −1 for those in November. Considering that the Laja Lake receives contributions from the Polcura River with flow rates that averaged 3.47 m 3 s −1 , 2.34 m 3 s −1 and 8,9 m 3 s −1 in January, March and November, respectively, it can be indicated that at least half of the water entering through the Polcura tunnel is lost by E.

Models
The evaluation of the combined group models showed that the best estimate was given by the Priestley-Taylor model. This model shows similarity in its equation with the deBrauin-Keijman model. Both models incorporate the variable γ, which is dependent on elevation. However, the deBrauin-Keijman model was developed at sea level, so the Priestley-Taylor model was more suitable for a mountain lake. Similar results were obtained by Wang et al. [31]. The Penman-Kimberly model overestimates E in the Laja Lake. This model is based on Penman and presents a modification of the wind function. Therefore, the wind function proposed by the Kimberly model was not adequate to estimate the hourly evaporation in the Laja Lake.
One of the lowest yields was given by the deBruin model. The author combined the Penman and Priestley-Taylor equations, and obtained a formula that predicts E. The model does not include net radiation, but the drying power of air (a function of the wind associated with vapor pressure gradient), and an energy component represented by the latent heat. deBruin's experimental data showed that the model performed well in a shallow body of water and in periods greater than 10 days [43]. The poor performance of this model in the Laja Lake can be associated with the exclusion of Q * and Q x in its formula. As abovementioned, Q x plays an important role in nighttime E.
The Brutsaert-Stricker model also presented a low fit with the measured data. Same as the models in its group, it bases its estimates on a radiative and an aerodynamic component. Brutsaert-Stricker [44] determined that the model was relatively insensitive to the choice of wind function, which would result in a poor performance of the model when applied to the Laja Lake. The wind function was adjusted for the Penman, deBruin and Brutsaert-Stricker models in order to maximize R 2 . The adjusted R 2 increased from 0.84 to 0.91; from 0.43 to 0.44; from 0.14 to 0.84 in Penman, deBruin and Brutsaert-Stricker, respectively. The adjustment of this last model resulted in the elimination of the aerodynamic component (a and b = 0 in the wind function).
The models of the solar radiation and temperature group overestimated E in the lake. These models base predictions on solar radiation and temperature, excluding net radiation, water heat flux, and the aerodynamic component. The overestimation can be explained by the effect of not including lake depth since part of the short-wave radiation is absorbed at the water surface, but the fraction not absorbed by the surface can penetrate deeper layers [41]. Due to the high heat capacity of the water, an important amount of energy is stored in the water during spring and summer months. Normally the water temperature increases and the heat storage peaks occurs in late summer or early fall [41,45].
The Ryan-Harleman, Trivett and Quinn models have adjustment constants calibrated for daily estimates. Due to this, it was necessary to divide the output of these models into 24 parts to obtain mm h −1 . The results obtained show that the Trivett model had the best fit. Additionally, the least squares method was applied to determine the constants a, b and c in the generalized form of the mass transfer models proposed by Sing and Xu [29]. Once the constants were determined and the model applied, it was obtained that the best fit was given by the G model. However, it is important to note that the D model presents similar fit indicators and requires fewer parameters. Although the models in this group tend to overestimate evaporation, they produced a good fit with the measured data. This occurred due to the influence of the wind on the E of the lake. Figure 14 shows the distribution of the model groups with respect to the indicators used. For RMSE versus R 2 , the closer to zero on the abscissa axis and to one on the ordinate axis, the better fit of the model. As for NS and W versus R 2 , a better fit is close to one on both axes. The models of the combined group had the best fit, followed by the mass transfer models. Conversely, the models based on solar radiation and temperature exhibited the worst performance.
This occurred due to the influence of the wind on the E of the lake. Figure 14 shows the distribution of the model groups with respect to the indicators used. For RMSE versus R 2 , the closer to zero on the abscissa axis and to one on the ordinate axis, the better fit of the model. As for NS and W versus R 2 , a better fit is close to one on both axes. The models of the combined group had the best fit, followed by the mass transfer models. Conversely, the models based on solar radiation and temperature exhibited the worst performance.  Figure 15 shows maximum, minimum, and average absolute errors between modeled and observed values for hourly estimates. It is observed that the mean value of most of the combination models and the mass transfer group was found within an error range of ± 0.04 mm h −1 , while the  Figure 15 shows maximum, minimum, and average absolute errors between modeled and observed values for hourly estimates. It is observed that the mean value of most of the combination models and the mass transfer group was found within an error range of ± 0.04 mm h −1 , while the models of the solar radiation and temperature group present a significant overestimation. The mean values of the absolute errors of the generalized form of mass transfer models exhibited the smallest distance to the observed value.  The use of the adjustment factors (KLL) resulted in a general improvement in all the models for the RMSE, W and NS indicators. However, the best fit was obtained using the combined group, followed by the mass transfer and the temperature and solar radiation group. Rosenberry et al. [30] evaluated 15 models to estimate evaporation in a small mountain lake on a daily scale, and determined that the best comparisons with the measured E were achieved with the methods that include available energy and aerodynamic components. Similar results were found in [31]. Figure 16 shows maximum, minimum and average absolute errors after applying the adjustment factor KLL. The mean value of the absolute errors improved, achieving an error range of ±0.02 mm h −1 for all the models. As in the previous case, the maximum amplitude of the error is found in the group of the solar radiation and temperature models, while the minimum amplitudes are achieved with the combined group. This shows the effectiveness of the use of an adjustment constant for the particular conditions of the lake. The use of the adjustment factors (K LL ) resulted in a general improvement in all the models for the RMSE, W and NS indicators. However, the best fit was obtained using the combined group, followed by the mass transfer and the temperature and solar radiation group. Rosenberry et al. [30] evaluated 15 models to estimate evaporation in a small mountain lake on a daily scale, and determined that the best comparisons with the measured E were achieved with the methods that include available energy and aerodynamic components. Similar results were found in [31]. Figure 16 shows maximum, minimum and average absolute errors after applying the adjustment factor K LL . The mean value of the absolute errors improved, achieving an error range of ±0.02 mm h −1 for all the models. As in the previous case, the maximum amplitude of the error is found in the group of the solar radiation and temperature models, while the minimum amplitudes are achieved with the combined group. This shows the effectiveness of the use of an adjustment constant for the particular conditions of the lake.

Conclusions
The main goals of this study were to evaluate the performance of 21 models to estimate hourly E rates in water bodies and to identify the main environmental variables that control its diurnal dynamics. Hourly E and environmental variables that affect it were measured in a mountain lake for 18 days using an EC system. Average E rates measured during three field campaigns reached 3.40, 3.38 and 1.89 mm d −1 in January, March and November 2016, respectively.
Best model performance was obtained with the models that include the available energy and aerodynamic components (combined group). The best fit was obtained with the Priestley-Taylor equation; however, the Penman model showed the best performance when the wind function was adjusted. The models based on solar radiation and air temperature exhibited the worst performance because the aerodynamic and available energy components were not considered, the Stephens-Stewart model produced the best estimates of the group. The mass transfer models presented a good agreement with the E measured. The greatest advantage of this group is that it does not require water heat flux measurements, the G model had the best performance in this group.
In order to adjust the results to the particular conditions of the Laja Lake, an adjustment factor (KLL) was applied to all models, with the exception of those proposed by Sing and Xu [29] as they present their own adjustment parameters. After applying the adjustment factor, the Penman model (with the wind function) presented the best performance. In general, all the models improved their indicators after applying the KLL.
Nighttime E represented 36.9 and 40.7% of the total daily rates in January and March, respectively, which is explained by the energy release from the water body to the atmosphere. The water heat flux during nighttime accounted for 64% of E. In November, the nighttime component reduced E rates to 26.4%, mainly due to the energy loss from the lake in winter. Daytime E is aerodynamically driven in the Laja Lake and hourly fluxes are strongly correlated with the product of wind speed above the lake and the vapor pressure deficit.
The analysis of the nighttime component of latent heat flux and the better correlation with wind speed found during daytime highlights the importance of the analysis of diurnal dynamics of water flux and energy stored in water. Therefore, hourly measurements of latent heat flux, net radiation, sensible heat flux, and water heat flux are recommended to better understand the E variability and energy balance components of water bodies.

Conclusions
The main goals of this study were to evaluate the performance of 21 models to estimate hourly E rates in water bodies and to identify the main environmental variables that control its diurnal dynamics. Hourly E and environmental variables that affect it were measured in a mountain lake for 18 days using an EC system. Average E rates measured during three field campaigns reached 3.40, 3.38 and 1.89 mm d −1 in January, March and November 2016, respectively.
Best model performance was obtained with the models that include the available energy and aerodynamic components (combined group). The best fit was obtained with the Priestley-Taylor equation; however, the Penman model showed the best performance when the wind function was adjusted. The models based on solar radiation and air temperature exhibited the worst performance because the aerodynamic and available energy components were not considered, the Stephens-Stewart model produced the best estimates of the group. The mass transfer models presented a good agreement with the E measured. The greatest advantage of this group is that it does not require water heat flux measurements, the G model had the best performance in this group.
In order to adjust the results to the particular conditions of the Laja Lake, an adjustment factor (K LL ) was applied to all models, with the exception of those proposed by Sing and Xu [29] as they present their own adjustment parameters. After applying the adjustment factor, the Penman model (with the wind function) presented the best performance. In general, all the models improved their indicators after applying the K LL .
Nighttime E represented 36.9 and 40.7% of the total daily rates in January and March, respectively, which is explained by the energy release from the water body to the atmosphere. The water heat flux during nighttime accounted for 64% of E. In November, the nighttime component reduced E rates to 26.4%, mainly due to the energy loss from the lake in winter. Daytime E is aerodynamically driven in the Laja Lake and hourly fluxes are strongly correlated with the product of wind speed above the lake and the vapor pressure deficit.
The analysis of the nighttime component of latent heat flux and the better correlation with wind speed found during daytime highlights the importance of the analysis of diurnal dynamics of water flux and energy stored in water. Therefore, hourly measurements of latent heat flux, net radiation, sensible heat flux, and water heat flux are recommended to better understand the E variability and energy balance components of water bodies.

Funding:
This research was supported by Chilean government through the project CRHIAM (ANID/FONDAP/15130015)