Lake Evaporation in a Hyper-Arid Environment , Northwest of China — Measurement and Estimation

Lake evaporation is a critical component of the hydrological cycle. Quantifying lake evaporation in hyper-arid regions by measurement and estimation can both provide reliable potential evaporation (ET0) reference and promote a deeper understanding of the regional hydrological process and its response towards changing climate. We placed a floating E601 evaporation pan on East Juyan Lake, which is representative of arid regions’ terminal lakes, to measure daily evaporation and conducted simultaneous bankside synoptic observation during the growing season of 2013–2015. A semi-empirical evaporation model derived from Dalton model was parameterized and validated with measured data. The model was then used to estimate lake evaporation during 2002–2015. According to in situ measurements, maximum, minimum and mean lake evaporation were 8.1, 3.7 and 6.5 mm/day, and growing season evaporation was 1183.3 mm (~80% of the annual amount). Adding up non-growing season evaporation that we converted from φ20 pan evaporation at Ejina weather station, the annual mean lake evaporation, 1471.3 mm, was representative of lower Heihe River’s ET0. Model inter-comparison implied our model performed well both in simplicity and accuracy and has potential utilization in a data-sparse area. In 2002–2015, estimated mean daily evaporation was 6.5 mm/day and growing season evaporation was 1233.7 mm. Trend analysis of estimated evaporation proved the evaporation paradox’s existence in this hyper-arid region and validated complementary relationship theory’s adaptability.


Introduction
Evaporation links the atmosphere and land surface system in the hydrological cycle [1][2][3].Water loss from evaporation is recognized as the largest component of the regional hydrological cycle in arid regions, especially for lakes or reservoirs [4][5][6].Open water evaporation, which includes lake evaporation, can also represent regional potential evaporation (PET) [7] and is often used as an important parameter in many methods to estimate terrestrial evapotranspiration (ET) [8][9][10].However, in remote areas, it is difficult to conduct in situ measurement to obtain reliable evaporation data due to tough conditions; many ET studies have thus been based on empirical or estimated values rather than measured PET data-an approximation that increases uncertainty and has some limitations [11][12][13][14][15][16].Quantifying evaporation is important for hydrological research on groundwater modeling and large-scale hydrological cycle simulation, and also helps us gain a better understanding of hydrological responses to climate change in extreme weather conditions [17][18][19].Thus, it is critical to provide fundamental and reliable open water evaporation reference for relevant ET research in validation, improvement of remote sensing methods to estimate ET, and seeking for practical methods for open water evaporation estimation in arid regions.
In situ measurement and model estimation are two main approaches for quantifying open water evaporation.The in situ measurement method originated in E. Halley's research in 1687, and different types of measuring instruments were subsequently developed in the following studies [20][21][22].The World Meteorological Organization suggests the following as standard equipment for measuring open water evaporation [23]: a 20 m 2 evaporation tank, a GGI-3000 from the Soviet Union [24,25], and a Class A evaporator from the USA [22,[26][27][28].The 20 m 2 evaporation tank is found only at some evaporation stations due to high expenses in building and maintaining the equipment, whereas the GGI-3000 is found only in former Soviet Union (FSU) countries, Russia and China [29].Compared with the evaporation pans previously mentioned, the E601 evaporation pan (i.e., modified GGI-3000 pan) is practical, cost-effective and commonly used by Chinese researchers [10,[30][31][32][33][34][35].One of the characteristics of the E601 pan is that it has four arc water troughs 20 cm wide, which comprise a water circle to reduce the effects of turbulence generated by the pan (particularly by the rim of the pan) while ensuring accuracy [29,36].However, for convenience and simplicity, researchers often locate the E601 pan on land rather than on the water, which may inaccurately represent open water evaporation due to different underlying conditions; therefore, evaporation data collected in this way must be converted to actual open water evaporation using a coefficient, which increases uncertainty in measured data [15,16,29].Researchers also applied the eddy covariance (EC) method, which is theoretically based on the correlation of vertical wind speed and air moisture content fluctuation and an accurate way to measure open water evaporation with fine accuracy [37][38][39][40][41][42][43].Although the EC method and pan method can both provide convincing and reasonable lake evaporation data [41,44], we cannot ignore that in a remote area with tough residence conditions, it is very hard to carry out in situ experiments, especially since EC equipment's durability and reliability may be affected by the extreme temperature (up to 41 • C and down to −36 • C) and strong wind (up to 16.5 m/s) in this hyper-arid study area [10], while the E601 pan, with three modifications, can meet our need of obtaining the actual lake evaporation with satisfying adaptability and affordable cost.Another approach for quantifying lake evaporation is model estimation; many studies have shown the reliability and practicality of this method, including the Dalton model [45], Bowen ratio-energy balance (BREB) method [46], Penman model [47], Priestley-Taylor model [48], water budget method [49], complementary relationship method [50][51][52] and energy balance method [53,54], etc. Table 1 shows a summary of relevant literature on open water evaporation researches with measurement and estimation methods [4,7,9,41,53,[55][56][57][58][59][60].
As discussed above, the modified floating E601 pan measurement is a more adaptive, reliable and affordable method to reveal actual lake evaporation in remote and tough areas compared with the EC method.Moreover, the floating E601 evaporation pan on the lake behaves better in revealing actual lake evaporation than evaporation pans on land [10,28,57].However, in extreme weather conditions (i.e., extreme dryness, hotness and coldness, or very high wind speed), due to the complexity and difficulty in instrument maintenance, this method is rarely used and it is not suitable for long-term in situ measurement.Meanwhile, evaporation models require a great deal of input data, which limits its utility in data-sparse regions.Concerned with such shortcomings, in this study, we directly measured daily lake evaporation using a floating E601 pan on the lake for the first time and sought a practical and low-input evaporation model to quantify open water evaporation in hyper-arid regions (Aridity Index < 0.05, AI = P/PET and P is mean annual precipitation [61]), taking East Juyan Lake (EJL) as an example, which is located in Ejina Delta, a hyper-arid and remote region in the northwest of China.We also obtained long time series of lake evaporation and added solid and more convincing proof of the evaporation paradox's existence in this region.In a wetland at the U.S. Geological Survey Norman Landfill Research Site in Oklahoma, USA; an area of 8800 m 2 .
In situ measurement; Monthly evaporation data from Class A pan and floating pan are used.
1. Floating pan can better simulate actual physical conditions on the water surface than evaporation pan on land.2. Floating pan to land pan ratios varied in different months, ranging from 0.69 to 0.87.
Bowen ratio; sensible heat (H) and basic synoptic data are used; Eddy covariance method is used.
1.A new method to estimate wet surface evaporation using easily obtained data was proposed based on Bowen ratio equation, avoiding the term of G in the energy budget that is hard to measure on wet surfaces.2. Validation by eddy covariance values indicated excellent accuracy when applied to LG.

Granger and Hedstrom (2011)
Three lakes with fetch distances ranging from 150 m to 11,000 m in Canada.
Direct measurement of eddy covariance method; wind speed, air temperature and humidity data are measured on the lake.
1. Wind speed is the most significant factor governing lake evaporation, followed by land-water contrast of temperature and vapor pressure, while net radiation bears no relationship with lake evaporation in short time periods.2. Relationships were developed between the hourly rates of lake evaporation and the following significant variables and parameters, and a versatile model for estimating hourly lake evaporation rate is proposed with good accuracy.3. Improvements can be made to the traditional mass transfer and DYRESM models for many applications, including modelling the effects small lakes have on regional weather.
Aerodynamic method and energy balance method; Synoptic factors are measured at the station located near the center of the lake.
1.The primary factors determining lake evaporation are solar radiation, near surface air and water temperatures, relative humidity, wind speed and near surface atmospheric stability.
2. Lake evaporation during exceptionally hot summer days was lower than normal when wind speed was much lower.
3. The weakening of the permanent synoptic Etesian winds and the descent of the marine inversion to a height of the topographic ridge are responsible for the reduction of the Mediterranean Sea Breeze and evaporation.In this research, we plan to utilize the in situ measurement method and propose a daily lake evaporation model with bankside synoptic data to measure as well as estimate the daily evaporation of EJL during the growing season (ice-free, from April to September) by three steps: (1) obtaining daily evaporation data of EJL from the floating E601 evaporation pan; (2) establishing a daily lake evaporation model for EJL with measured synoptic data; and (3) utilizing the proposed method to estimate EJL's daily evaporation during 2002-2015.

Study Area
East Juyan Lake (101 • 11 31"-101 • 19 07" E, 42 • 15 25"-42 • 20 00" N) is located ~60 km north of Ejina weather station, which is a hyper-arid region in the lower reach of the Heihe River-the second largest inland river in northwest China.Ejina Delta is characterized by a continental climate that is extremely hot in the summer and cold in the winter.The maximum daily temperature was 41 • C (in July) and the minimum daily temperature was −36 • C (in January).There are strong winds, especially in spring and autumn, with maximum wind speed of 16.5 m/s [62][63][64].According to data from the Ejina weather station from 1959 to 2015, the mean annual temperature, precipitation, pan evaporation, relative humidity and wind speed were 9.0 • C, 36.6 mm, 3755 mm (φ20 evaporation pan, which is a copper pan with a diameter of 20 cm and height of 10 cm), 33.7% and 3.3 m/s, respectively.The growing season in the delta is from April to September, during which time it is ice free.Ejina Delta' s dominant landscape is the Gobi Desert (Figure 1), with seasonal runoff in the east river and west river [65,66].
Water 2016, 8, 527 5 of 21 In this research, we plan to utilize the in situ measurement method and propose a daily lake evaporation model with bankside synoptic data to measure as well as estimate the daily evaporation of EJL during the growing season (ice-free, from April to September) by three steps: (1) obtaining daily evaporation data of EJL from the floating E601 evaporation pan; (2) establishing a daily lake evaporation model for EJL with measured synoptic data; and (3) utilizing the proposed method to estimate EJL's daily evaporation during 2002-2015.

Study Area
East Juyan Lake (101°11′31″-101°19′07″ E, 42°15′25″-42°20′00″ N) is located ~60 km north of Ejina weather station, which is a hyper-arid region in the lower reach of the Heihe River-the second largest inland river in northwest China.Ejina Delta is characterized by a continental climate that is extremely hot in the summer and cold in the winter.The maximum daily temperature was 41 °C (in July) and the minimum daily temperature was −36 °C (in January).There are strong winds, especially in spring and autumn, with maximum wind speed of 16.5 m/s [62][63][64].According to data from the Ejina weather station from 1959 to 2015, the mean annual temperature, precipitation, pan evaporation, relative humidity and wind speed were 9.0 °C, 36.6 mm, 3755 mm (φ20 evaporation pan, which is a copper pan with a diameter of 20 cm and height of 10 cm), 33.7% and 3.3 m/s, respectively.The growing season in the delta is from April to September, during which time it is ice free.Ejina Delta' s dominant landscape is the Gobi Desert (Figure 1), with seasonal runoff in the east river and west river [65,66].EJL, surrounded by the Gobi Desert, is a shallow lake lying at the end of the east river with the mean water depth of 1.8 m.According to remote sensing data, the mean area of EJL during the past decade was ~40 km 2 , with a maximum area of 59.6 km 2 [13].It must be noted that EJL dried up in 1992 due to decreasing river runoff in the lower Heihe River [67,68].Since 2002, EJL has been recharged under the official ecological water conveyance project (EWCP), which was launched to transport water from the upper and middle reaches of the Heihe River to the lower reach-Ejina Delta-to mitigate ecological problems such as lake shrinkage, groundwater decline and vegetation degradation in EJL and the associated riparian zone [69,70].Annually, ~45 million m 3 water is conveyed to EJL.

Data Sources
In this study, synoptic data and evaporation data are used.Synoptic data includes data from the Ejina weather station and EJL bankside weather station that we have set up.Evaporation data includes E601 pan evaporation data at Ejina weather station (on land) and floating E601 pan evaporation on EJL.
The Ejina weather station is maintained by the local weather bureau and is 60 km away from EJL (Figure 1).Because data from this station may not reflect the actual climate conditions of the lake, we set up the East Juyan Lake bankside weather station 20 m away from the lake (Figure 2) to automatically observe precipitation, air temperature, wind speed (at 2 m), relative humidity, solar radiation and atmospheric pressure in 5 min intervals.The data are transferred by the DTU (Data Transfer Unit) online in 60 min intervals.Daily synoptic data were auto-recorded from 25 August 2013 to 30 September 2015 in the growing season, except for occasional gaps due to equipment being out of order.Daily synoptic data from the Ejina weather station are thus available accordingly.
EJL, surrounded by the Gobi Desert, is a shallow lake lying at the end of the east river with the mean water depth of 1.8 m.According to remote sensing data, the mean area of EJL during the past decade was ~40 km 2 , with a maximum area of 59.6 km 2 [13].It must be noted that EJL dried up in 1992 due to decreasing river runoff in the lower Heihe River [67,68].Since 2002, EJL has been recharged under the official ecological water conveyance project (EWCP), which was launched to transport water from the upper and middle reaches of the Heihe River to the lower reach-Ejina Delta-to mitigate ecological problems such as lake shrinkage, groundwater decline and vegetation degradation in EJL and the associated riparian zone [69,70].Annually, ~45 million m 3 water is conveyed to EJL.

Data Sources
In this study, synoptic data and evaporation data are used.Synoptic data includes data from the Ejina weather station and EJL bankside weather station that we have set up.Evaporation data includes E601 pan evaporation data at Ejina weather station (on land) and floating E601 pan evaporation on EJL.
The Ejina weather station is maintained by the local weather bureau and is 60 km away from EJL (Figure 1).Because data from this station may not reflect the actual climate conditions of the lake, we set up the East Juyan Lake bankside weather station 20 m away from the lake (Figure 2) to automatically observe precipitation, air temperature, wind speed (at 2 m), relative humidity, solar radiation and atmospheric pressure in 5 min intervals.The data are transferred by the DTU (Data Transfer Unit) online in 60 min intervals.Daily synoptic data were auto-recorded from 25 August 2013 to 30 September 2015 in the growing season, except for occasional gaps due to equipment being out of order.Daily synoptic data from the Ejina weather station are thus available accordingly.The evaporation data include E601 pan evaporation data from Ejina weather station from 2002 to 2015, and the floating E601 evaporation pan (Figure 2) from 25 August 2013 to 30 September 2015.The pan is located on the lake approximately 150 m away from the bank (minimum water depth >1.2 m, 101 • 15 42" E and 42 • 19 45" N) and was maintained weekly.The sensor can measure the water level of the E601 pan in 5 min intervals with 0.1 mm precision, and when the water level of the pan reaches a low value due to evaporation, the pump will automatically add water to maintain a reasonable water level (the distance between water surface to top of the pan's edge is less than 15 cm).The floating E601 pan water level data from the sensor yielded 239 days of pan evaporation data (except for periods of ice or equipment's out of work).

Data Quality Control
In case of occasional strong wind and high waves, we have designed an anti-wave steel boat that can prevent pan inflow and water level fluctuation.The E601 pan is floating in the center of the boat with arms stretching out for 5 m in four directions, and the boat is fixed by four steel tubes in the lake to allow only vertical displacement, making the big boat a relatively steady platform whereby the water inside the boat is almost still and the water level data is hardly affected by high waves (Figure 2c).In addition, the E601 pan is floating on the water that is freely connected to the water outside of the boat.The pan is thus actually floating on the lake inside the boat and the steel plates (40 cm in height) around the main body of the boat.Due to a distance of 50 cm from the steel plate to the edge of the E601 pan, pan inflow can be effectively prevented (Figure 2d).For unpredicted wind and high waves, pan inflow can be identified and eliminated after quality control of processing high measuring frequency water level data.Such measures mentioned above can ensure the reliability of the pan evaporation data we obtained.
We calculated the daily evaporation by considering the starting time water level, the ending time water level, and the upward (water out of the pan) and downward (water into the pan) water level jump during a day.Water jump is caused by pumping water, occasional precipitation, splashed water etc., and can be easily identified with the help of water level data in high frequency (5 min steps), since water jump can be 0.3-0.8mm within one or two time steps, whereas the evaporation rate scarcely exceeds 0.2 mm within 10 min.

Brief Description of the E601 Pan
The widely used E601 pan is made of fiber-reinforced polymer and is of cylindrical design with a surface area of 3000 cm 2 and depth of 60 cm.It is set in the soil with its rim 30 cm above the ground and with a cone-shaped or a cylindrical-shaped bottom, and there is a needle ruler attached at the rim of the pan to measure water level change caused by evaporation [71].Apart from the identical designs and dimensions with GGI-3000 pans, there are two modifications in the E601 pan: (1) when installing, the E601 pan is supported 22.5 cm higher than the GGI-3000 pan, making the rim 30 cm above the ground; (2) there is a water ring around the pan consisting of four arc water troughs 20 cm in width, giving E601 pan the advantage of reducing the edge effects of turbulence generated by the rim of the pan over the Russian GGI-3000 pan [36].
E601 pan is a frequently used way to measure evaporation, despite its remaining shortcomings.In this study, concerned with the disadvantages of the traditional E601 pan, we made three major modifications: (1) In order to have the pan in the lake rather than on land, we designed a big steel boat with four 5-m-long arms, placing the pan inside this boat's main body on the lake for less interference from waves (Figure 2d).Consequently, the measured pan evaporation is more representative of actual lake evaporation than data from evaporation pans on land; (2) We replaced the needle ruler with an auto-recorded water level sensor to measure evaporation, the sensor set at the rim is of high accuracy (0.1 mm precision) and can record water level data automatically, avoiding errors caused during manual operation and reading, thereby, ensuring water level data accuracy; (3) We also set up a pump that can automatically add water into the pan until it reaches a given water level at midnight, if the water was reduced to a certain level during the daily evaporation, thereby diminishing possible sudden heat changes from the added water during daytime and simplifying the measurement process.
In short, the evaporation pan we used can represent actual lake evaporation and ensure fine data accuracy without intensive maintenance during the experiment.

Derivation of the Model
Lake evaporation process is governed by a number of factors, and vapor pressure deficit (VPD), air temperature (T a ) and wind speed (u) are the three main synoptic factors [45,60,72].Dalton's traditional method to estimate open water evaporation has low-input demand.However, after applying Dalton's equation to calculate EJL's daily evaporation, we found the calculation a little complicated and the modeling accuracy too low (R 2 is below 0.5).Thus, to improve Dalton method's modeling accuracy and adaptability, we derived a semi-empirical lake evaporation model based on the Dalton model, in which the original form can be expressed as: where E is the open water evaporation, mm; f (u) is the wind speed function and e s − e a is VPD, kPa; e s is the saturation vapor pressure at the temperature of water surface (T s ) and e a is the vapor pressure at air temperature (T a , at 2 m).kPa are empirically expressed as: e a = RH 100 e(T s ) where RH is relative humidity, %, T s is evaporation surface temperature and T a is air temperature at 2 m, • C. Due to the lack of accessible surface temperature data, we use measured air temperature and assume that the temperature gradient between water surface and the measuring level of T a is zero [59,73], which implies e s (T s ) = e s (T a ), and e a can be expressed empirically as the function of T a , thus obtaining: where g(RH) and h(T a ) are function of relative humidity and air temperature, yielding: Wind speed has a positive relationship with E in transporting evaporated vapor to neighborhood, and the function can be expressed as: High RH restrains the evaporation process as moisture accumulates in the air during the evaporation process, creating a negative relationship with E; high T a intensifies the escaping process of Water 2016, 8, 527 9 of 21 liquid water modular on the evaporation surface, so there is a positive relationship with evaporation.The g(RH) and h(T a ) function can be generalized as: where a, b, c, d, m and n are undermined parameters and α, β, γ are often valued according to local synoptic conditions.The power of α often ranges from 0.5 to 1.0, and according to the facts that measured wind speed data bears no good linear relationship with daily evaporation and the wind speed on EJL is quite high, we choose the value to be 0.5; β reflects the relative humidity's sensitivity regarding evaporation and ranges from 0.5 to 2.0, so we choose the value to be 1.5 considering that atmospheric capacity of vapor restrains the process of lake evaporation in arid regions considerably; γ ranges from 0.5 to 2.0 depending on different occasions, and considering that air temperature during the growing season ranges from −10 • C to 40 • C; the commonly used Equation (2) [74] can therefore be replaced by Equation ( 8) with the power of 1.0 for simplicity.
Thus, we generated a multivariate nonlinear evaporation model consisting of T a , u and RH.Parameterization and validation can be done with measured daily data from bankside weather station and floating evaporation pan.
where σ is the Stephan-Boltzman constant =4.933 × 10 −9 (MJ/m 2 /K 4 /day), α is reflection coefficient or albedo and for open water surface α = 0.08, R S is solar radiation and R SO is clear sky radiation (W•m −2 ), calculated as Water 2016, 8, 527 10 of 21 For East Juyan Lake, the calibrated value of a, b are 0.45 and 0.39, Z is elevation and R A is upper (outer limits of the atmosphere) solar radiation, W•m −2 .
The other meteorological parameters that appear in the two methods are calculated as follows: 3. Results

Time Series of East Juyan Lake's Daily Evaporation
After quality control, we processed the water level data and obtained 239 daily evaporation data values (mainly from August to September in 2013, May to August in 2014 and April to September in 2015) (Figure 3).
For East Juyan Lake, the calibrated value of a, b are 0.45 and 0.39, Z is elevation and RA is upper (outer limits of the atmosphere) solar radiation, W•m −2 .
The other meteorological parameters that appear in the two methods are calculated as follows:  

Time Series of East Juyan Lake's Daily Evaporation
After quality control, we processed the water level data and obtained 239 daily evaporation data values (mainly from August to September in 2013, May to August in 2014 and April to September in 2015) (Figure 3).The daily evaporation data cover most days from April to September in 2013-2015, along with synoptic data measured at the EJL bankside weather station.The mean daily lake evaporation during the growing season is 6.5 mm, the maximum and minimum daily evaporation are 8.1 mm (in August) and 3.7 mm (in September).
To obtain monthly evaporation of EJL, since data are occasionally unavailable in a month, we assume that if there are 15 or more measured daily evaporation scattering in each month, the mean daily evaporation of these available data is representative of the mean daily evaporation of this month, so the monthly evaporation can be calculated.Using the available data, we obtained the monthly lake evaporation of EJL (Table 3), from which we conclude that the mean total lake evaporation during annual growing season is 1183.3mm.The daily evaporation data cover most days from April to September in 2013-2015, along with synoptic data measured at the EJL bankside weather station.The mean daily lake evaporation during the growing season is 6.5 mm, the maximum and minimum daily evaporation are 8.1 mm (in August) and 3.7 mm (in September).
To obtain monthly evaporation of EJL, since data are occasionally unavailable in a month, we assume that if there are 15 or more measured daily evaporation scattering in each month, the mean daily evaporation of these available data is representative of the mean daily evaporation of this month, so the monthly evaporation can be calculated.Using the available data, we obtained the monthly lake evaporation of EJL (Table 3), from which we conclude that the mean total lake evaporation during annual growing season is 1183.3mm.As described before, the semi-empirical daily lake evaporation for EJL is expressed as: In deciding the parameters of a, b, c, d, m and n, we used u, RH and T a data measured at the EJL bankside weather station and floating E601 pan evaporation data of EJL from 2013 to 2014 (77 days of available data).After parameterization, the model can be expressed as: E 0 = 0.0345 + 0.002u 0.5 42.6824 − 0.0122RH 1.5 (2.66 + 0.08T a ) Comparing with measured daily lake evaporation data, we see good accuracy in estimating daily evaporation of EJL using the proposed model, and the mean absolute error (MAE) is 0.52 mm, root-mean-square error (RMSE) is 0.63 mm, R 2 is 0.84 and the mean relative error (MRE) is 8.6% (Table 4).
Validation is also conducted.We used measured bankside synoptic data (T, W and RH) from April to August in 2015 to estimate the daily EJL evaporation with the proposed model, and compare with measured data from the floating E601 pan (114 days of available data) (Figure 4).Meanwhile, we compared the obtained monthly lake evaporation with estimated monthly evaporation, which is calculated from the daily evaporation estimation.As described before, the semi-empirical daily lake evaporation for EJL is expressed as: In deciding the parameters of a, b, c, d, m and n, we used u, RH and Ta data measured at the EJL bankside weather station and floating E601 pan evaporation data of EJL from 2013 to 2014 (77 days of available data).After parameterization, the model can be expressed as: Comparing with measured daily lake evaporation data, we see good accuracy in estimating daily evaporation of EJL using the proposed model, and the mean absolute error (MAE) is 0.52 mm, root-mean-square error (RMSE) is 0.63 mm, R 2 is 0.84 and the mean relative error (MRE) is 8.6% (Table 4).
Validation is also conducted.We used measured bankside synoptic data (T, W and RH) from April to August in 2015 to estimate the daily EJL evaporation with the proposed model, and compare with measured data from the floating E601 pan (114 days of available data) (Figure 4).Meanwhile, we compared the obtained monthly lake evaporation with estimated monthly evaporation, which is calculated from the daily evaporation estimation.

Inter-Comparison of Evaporation Methods
The deBruin-Keijman method and Brutsaert-Stricker method are believed to estimate daily open water evaporation with fine accuracy [56], so we applied them to calculate EJL's daily evaporation and compare the results with measured evaporation (Emea) (Figure 5 and Table 5).
Water 2016, 8, 527 12 of 21 The deBruin-Keijman method and Brutsaert-Stricker method are believed to estimate daily open water evaporation with fine accuracy [56], so we applied them to calculate EJL's daily evaporation and compare the results with measured evaporation (Emea) (Figure 5 and Table 5).Our proposed model tends to yield better estimations of daily evaporation than the d-K and B-S methods, and d-K behaves better than B-S in this hyper-arid environment, with the MAE and RMSE being 0.81 mm and 0.96 mm.In general, our proposed model can estimate daily evaporation of EJL with moderate accuracy, while the d-K method tends to underestimate daily evaporation and the B-S method tends to overestimate (Figure 5).

Estimation of EJL's Evaporation during 2002-2015
Synoptic data (u, RH and Ta) are needed to estimate time series of EJL's evaporation from 2002 to 2015.Input synoptic data near the lake helps to estimate the evaporation more accurately [58].However, the EJL bankside weather station was set up in 2013 and gathered only three years of synoptic data.Because the Ejina weather station is rich in accessible data , we expanded EJL's synoptic data series based on the correlativity (Table 6, Figure 6) of available data (2013-2015) from these two stations to get longer time series of synoptic data (2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015), with which we estimated EJL's daily evaporation.Although the Ejina weather station is 60 km away from EJL, air temperature, wind speed and relative humidity measured at these two stations are well correlated and therefore a longer series of EJL bankside synoptic data can be obtained using the correlations.During 2002-2015, air  Our proposed model tends to yield better estimations of daily evaporation than the d-K and B-S methods, and d-K behaves better than B-S in this hyper-arid environment, with the MAE and RMSE being 0.81 mm and 0.96 mm.In general, our proposed model can estimate daily evaporation of EJL with moderate accuracy, while the d-K method tends to underestimate daily evaporation and the B-S method tends to overestimate (Figure 5).

Estimation of EJL's Evaporation during 2002-2015
Synoptic data (u, RH and T a ) are needed to estimate time series of EJL's evaporation from 2002 to 2015.Input synoptic data near the lake helps to estimate the evaporation more accurately [58].However, the EJL bankside weather station was set up in 2013 and gathered only three years of synoptic data.Because the Ejina weather station is rich in accessible data , we expanded EJL's synoptic data series based on the correlativity (Table 6, Figure 6) of available data (2013-2015) from these two stations to get longer time series of synoptic data (2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)(2010)(2011)(2012)(2013)(2014)(2015), with which we estimated EJL's daily evaporation.temperature, mean annual wind speed and relative humidity have a trend of 0.013 °C/year, −0.036 m/s/year and −0.46%/year, respectively (Figure 7).For the growing season of 2002-2015 the estimated mean, maximum and minimum daily lake evaporation are 6.5 mm, 7.8 mm and 3.4 mm; the estimated mean, maximum and minimum total lake evaporation are 1233.7 mm, 1190.7 mm (in 2003) and 1258.3 mm (in 2009), and there is an increasing trend in annual total evaporation of 0.83 mm/year.Generally, the maximum monthly lake evaporation appears in July or August, and the minimum monthly lake evaporation appears in April, which both accord with the temporal distribution of in situ measured monthly lake evaporation (Table 3).

Floating Pan Evaporation's Representativeness and Reliability
During the growing season of 2013-2015, EJL's total measured evaporation was 1183.3 mm.To obtain the annual open water evaporation, we had to consider the evaporation during the non-growing season (October to March, the next year), when it is the ice period and measuring instruments cannot work.Therefore, we converted the φ20 evaporation pan data from the Ejina weather station with the empirical conversion coefficient of 0.48 [77] to obtain the evaporation in this period.For 2013-2015, the mean converted total evaporation during the non-growing season was ~288.0 mm, making the mean annual open water evaporation (ET0) 1471.3 mm.Lake evaporation during the growing season took up ~80% of the annual amount, while it was only ~20% during the non-growing season.Lake evaporation during the growing season was in situ measured with high frequency and high reliability; although we could not neglect uncertainty of pan evaporation conversion method used to get the open water evaporation in non-growing season, its influence on the annual lake evaporation was quite small, considering its small portion (20%) and the error of conversion method (within 15% [78,79]).Consequently, we are confident that the combined annual evaporation, 1471.3 mm, is quite convincing as the annual lake evaporation and standard representation of the lower Heihe River's annual ET0.Although the Ejina weather station is 60 km away from EJL, air temperature, wind speed and relative humidity measured at these two stations are well correlated and therefore a longer series of EJL bankside synoptic data can be obtained using the correlations.During 2002-2015, air temperature, mean annual wind speed and relative humidity have a trend of 0.013 • C/year, −0.036 m/s/year and −0.46%/year, respectively (Figure 7).For the growing season of 2002-2015 the estimated mean, maximum and minimum daily lake evaporation are 6.5 mm, 7.8 mm and 3.4 mm; the estimated mean, maximum and minimum total lake evaporation are 1233.7 mm, 1190.7 mm (in 2003) and 1258.3 mm (in 2009), and there is an increasing trend in annual total evaporation of 0.83 mm/year.Generally, the maximum monthly lake evaporation appears in July or August, and the minimum monthly lake evaporation appears in April, which both accord with the temporal distribution of in situ measured monthly lake evaporation (Table 3).

Floating Pan Evaporation's Representativeness and Reliability
During the growing season of 2013-2015, EJL's total measured evaporation was 1183.3 mm.To obtain the annual open water evaporation, we had to consider the evaporation during the non-growing season (October to March, the next year), when it is the ice period and measuring instruments cannot work.Therefore, we converted the φ20 evaporation pan data from the Ejina weather station with the empirical conversion coefficient of 0.48 [77] to obtain the evaporation in this period.For 2013-2015, the mean converted total evaporation during the non-growing season was ~288.0 mm, making the mean annual open water evaporation (ET0) 1471.3 mm.Lake evaporation during the growing season took up ~80% of the annual amount, while it was only ~20% during the non-growing season.Lake evaporation during the growing season was in situ measured with high frequency and high reliability; although we could not neglect uncertainty of pan evaporation conversion method used to get the open water evaporation in non-growing season, its influence on the annual lake evaporation was quite small, considering its small portion (20%) and the error of conversion method (within 15% [78,79]).Consequently, we are confident that the combined annual evaporation, 1471.3 mm, is quite convincing as the annual lake evaporation and standard representation of the lower Heihe River's annual ET0.For the growing season of 2002-2015 the estimated mean, maximum and minimum daily lake evaporation are 6.5 mm, 7.8 mm and 3.4 mm; the estimated mean, maximum and minimum total lake evaporation are 1233.7 mm, 1190.7 mm (in 2003) and 1258.3 mm (in 2009), and there is an increasing trend in annual total evaporation of 0.83 mm/year.Generally, the maximum monthly lake evaporation appears in July or August, and the minimum monthly lake evaporation appears in April, which both accord with the temporal distribution of in situ measured monthly lake evaporation (Table 3).

Floating Pan Evaporation's Representativeness and Reliability
During the growing season of 2013-2015, EJL's total measured evaporation was 1183.3 mm.To obtain the annual open water evaporation, we had to consider the evaporation during the non-growing season (October to March, the next year), when it is the ice period and measuring instruments cannot work.Therefore, we converted the φ20 evaporation pan data from the Ejina weather station with the empirical conversion coefficient of 0.48 [77] to obtain the evaporation in this period.For 2013-2015, the mean converted total evaporation during the non-growing season was ~288.0 mm, making the mean annual open water evaporation (ET 0 ) 1471.3 mm.Lake evaporation during the growing season took up ~80% of the annual amount, while it was only ~20% during the non-growing season.Lake evaporation during the growing season was in situ measured with high frequency and high reliability; although we could not neglect uncertainty of pan evaporation conversion method used to get the open water evaporation in non-growing season, its influence on the annual lake evaporation was quite small, considering its small portion (20%) and the error of conversion method (within 15% [78,79]).Consequently, we are confident that the combined annual evaporation, 1471.3 mm, is quite convincing as the annual lake evaporation and standard representation of the lower Heihe River's annual ET 0 .
Open water evaporation, reference evapotranspiration and pan evaporation are used to represent regional ET 0 [52,[79][80][81].We converted E601 and φ20 pan evaporation at Ejina weather station with the conversion coefficient of 0.76 and 0.48 (referred to a research based on similar arid climate settings in Tarim River Basin [77]), the results demonstrating that the lower Heihe River's ET 0 represented by pan evaporation was 1565.9 mm.According to recent studies, calculated reference evapotranspiration in the lower Heihe River using the Food and Agriculture Organization (FAO) Penman-Monteith method were 1631 mm by Li [82] in 2013 and 1553 mm by Du [83] in 2016.Table 7 shows there is acceptable difference (less than 10%) when using open water evaporation, E601 pan evaporation and reference evapotranspiration to represent ET 0 of the lower Heihe River.

Model Sensitivity and Reliability
The d-K and B-S methods tended to demonstrate larger uncertainty in estimating EJL's daily evaporation than the proposed model.The reason may be that these two methods have considered both aerodynamic and energy terms in evaporation, which requires more input variables, but in a data-sparse region like the Ejina Delta, some parameters such as R n , Q x and VPD are unavailable, and hence have to be calculated using empirical methods, which increases the uncertainty due to these variables' high sensitivity [84,85].On the contrary, the proposed model involves basic and accessible synoptic variables (T a , u and RH) and is with moderate accuracy.The proposed model considered three governing synoptic factors as model inputs, and they accorded well with Granger and Hedstrom's finding that wind speed, air temperature and relative humidity are the most significant factors governing evaporation rate, after doing in situ evaporation measurement in three lakes [60].Therefore, the proper model variables we selected could ensure the proposed model's reliability.Moreover, this model (Equation ( 20)), derived from the Dalton model, balances the simplicity of input data and accuracy of evaporation estimation well, and we are confident that, after parameterization with local input data, it is suitable for lakes in China's vast arid and data-sparse regions.
To find the most sensitive variables in the model, we analyzed their sensitivity coefficients, according to EJL's synoptic data, and estimated lake evaporation during the growing season of 2002-2015 (Figure 8).The sensitivity coefficient (S v , dimensionless) vividly represents the changes in actual lake evaporation caused by any perturbation of the governing synoptic factors [86], and is defined as [87]: where ∆V i is the change of synoptic variables (T a , u and RH), and ∆E is the change of evaporation.For 2002-2015, the mean S v of T a , u and RH are 0.67, 0.51 and −0.36, respectively, indicating that air temperature has the highest sensitivity, followed by wind speed and relative humidity (negative).In Granger and Hedstrom's research, wind speed was the most significant factor in lake evaporation [60].The climatic settings of the two sites are different.Crean Lake is a small lake located in Canada with mild weather conditions (annual precipitation greater than 500 mm), while EJL is located in a hyper-arid region with extreme weather conditions (annual precipitation of 36.6 mm).The environmental factors (Ta, u and atmospheric turbulence) governing evaporation of water bodies can be substantially distinctive over small water bodies than over larger lakes [3,58].Secondly, air temperature and wind speed are all very significant factors in open water evaporation; for EJL, a shallow lake (depth of ~2 m), compared with the deep Crean Lake (depth of ~20 m), the water is easily heated thoroughly (mean daily air temperature of 21.6 °C during the growing season) and, with the help of high wind speed (mean value of 3.4 m/s) on the lake, air temperature shows even greater significance in lake evaporation than Crean Lake with its obvious temperature gradients of different water levels.Therefore, Ta is the most sensitive factor, followed by u in EJL.

Complementary Relationship Theory's Adaptability in the Hyper-Arid Region
The complementary relationship (CR) theory, used by Brutsaert to explain the evaporation paradox, suggests that in non-humid regions, the increase of temperature results in the increase of actual evaporation and the decrease of observed pan evaporation [50].According to Figures 7a and  9b, the variation trends of Ta and E for the growing season of 2002-2015 were 0.013 °C/year and 0.83 mm/year, indicating that there was an increasing trend of lake evaporation (representative of regional actual evaporation) along with the increase of air temperature, whereas, annual E601 pan evaporation (during the growing season) at Ejina weather station decreased over the last 14 years (Figure 10).In Granger and Hedstrom's research, wind speed was the most significant factor in lake evaporation [60].The climatic settings of the two sites are different.Crean Lake is a small lake located in Canada with mild weather conditions (annual precipitation greater than 500 mm), while EJL is located in a hyper-arid region with extreme weather conditions (annual precipitation of 36.6 mm).The environmental factors (T a , u and atmospheric turbulence) governing evaporation of water bodies can be substantially distinctive over small water bodies than over larger lakes [3,58].Secondly, air temperature and wind speed are all very significant factors in open water evaporation; for EJL, a shallow lake (depth of ~2 m), compared with the deep Crean Lake (depth of ~20 m), the water is easily heated thoroughly (mean daily air temperature of 21.6 • C during the growing season) and, with the help of high wind speed (mean value of 3.4 m/s) on the lake, air temperature shows even greater significance in lake evaporation than Crean Lake with its obvious temperature gradients of different water levels.Therefore, T a is the most sensitive factor, followed by u in EJL.

Complementary Relationship Theory's Adaptability in the Hyper-Arid Region
The complementary relationship (CR) theory, used by Brutsaert to explain the evaporation paradox, suggests that in non-humid regions, the increase of temperature results in the increase of actual evaporation and the decrease of observed pan evaporation [50].According to Figures 7a and  9b, the variation trends of T a and E for the growing season of 2002-2015 were 0.013 • C/year and 0.83 mm/year, indicating that there was an increasing trend of lake evaporation (representative of regional actual evaporation) along with the increase of air temperature, whereas, annual E601 pan evaporation (during the growing season) at Ejina weather station decreased over the last 14 years (Figure 10).In Granger and Hedstrom's research, wind speed was the most significant factor in lake evaporation [60].The climatic settings of the two sites are different.Crean Lake is a small lake located in Canada with mild weather conditions (annual precipitation greater than 500 mm), while EJL is located in a hyper-arid region with extreme weather conditions (annual precipitation of 36.6 mm).The environmental factors (Ta, u and atmospheric turbulence) governing evaporation of water bodies can be substantially distinctive over small water bodies than over larger lakes [3,58].Secondly, air temperature and wind speed are all very significant factors in open water evaporation; for EJL, a shallow lake (depth of ~2 m), compared with the deep Crean Lake (depth of ~20 m), the water is easily heated thoroughly (mean daily air temperature of 21.6 °C during the growing season) and, with the help of high wind speed (mean value of 3.4 m/s) on the lake, air temperature shows even greater significance in lake evaporation than Crean Lake with its obvious temperature gradients of different water levels.Therefore, Ta is the most sensitive factor, followed by u in EJL.

Complementary Relationship Theory's Adaptability in the Hyper-Arid Region
The complementary relationship (CR) theory, used by Brutsaert to explain the evaporation paradox, suggests that in non-humid regions, the increase of temperature results in the increase of actual evaporation and the decrease of observed pan evaporation [50].According to Figures 7a and  9b, the variation trends of Ta and E for the growing season of 2002-2015 were 0.013 °C/year and 0.83 mm/year, indicating that there was an increasing trend of lake evaporation (representative of regional actual evaporation) along with the increase of air temperature, whereas, annual E601 pan evaporation (during the growing season) at Ejina weather station decreased over the last 14 years (Figure 10).The evaporation paradox suggests that with a warming climate, a declining rather than increasing trend of observed pan evaporation (i.e., apparent potential evaporation) during the past 50 years has been witnessed worldwide, especially in the Northern Hemisphere [50,88,89].However, the "pan evaporation paradox" is a manifestation of the CR theory, as major declines of pan evaporation were witnessed in many sites in the US, Canada, Mexico, Iran, Thailand and India [90][91][92][93][94][95][96], many of which are in arid environmental settings.Similar proof for the evaporation paradox was reported in China's northwestern arid regions [33,[97][98][99].Set against a global background of the evaporation paradox, our research results-i.e., in situ measurements reliably estimated using our proposed model-indicate an increase of actual evaporation and decrease of observed pan evaporation in accordance with the available evidence; furthermore, it adds concrete proof of the existence of the evaporation paradox, as well as to the adaptability of the CR theory in China's hyper-arid region.

Conclusions
Quantifying lake evaporation is a theoretical and practical issue, especially in a hyper-arid region, as it offers important and reliable reference for relevant research.Through fine and fundamental evaporation measurement on East Juyan Lake (EJL) in the growing season of 2013-2015, EJL's measured mean total evaporation was 1183.3 mm, and the maximum, minimum, mean daily lake evaporation were 8.1, 3.7, 6.5 mm, respectively.The measured annual open water evaporation was 1471.3 mm; along with E601 pan evaporation and the calculated reference evapotranspiration, these measurements well represent the lower Heihe River's ET0 with little differences.The measurement of lake evaporation in growing season was solid and reliable, while measured data in non-growing season is lacking, so we converted the data from φ20 pan evaporation at Ejina weather station.Although such a conversion method has unneglectable errors, the results can still be favorable (less than 3% of uncertainty) since non-growing season evaporation only has a small portion (~20%) of annual evaporation; therefore, the measured open water evaporation as potential evaporation (ET0) is convincing.Moreover, this is an improvement to finding a more reliable and practical method to quantify lake evaporation (sublimation) during the non-growing season.
Derived from the Dalton model, parameterized and validated by in situ data, in addition to being compared with acknowledged evaporation models, the proposed daily lake evaporation model (Equation ( 21)) has the advantages not only of simplicity, reliability and accuracy, but also of potential utilization for other lakes in arid and data-sparse regions.Furthermore, apart from the governing synoptic factors considered in this study, lake area, water depth and water quality are also important factors in lake evaporation [8,100,101] that could be used for modifying the proposed model for wider application in arid regions worldwide.The evaporation paradox suggests that with a warming climate, a declining rather than increasing trend of observed pan evaporation (i.e., apparent potential evaporation) during the past 50 years has been witnessed worldwide, especially in the Northern Hemisphere [50,88,89].However, the "pan evaporation paradox" is a manifestation of the CR theory, as major declines of pan evaporation were witnessed in many sites in the US, Canada, Mexico, Iran, Thailand and India [90][91][92][93][94][95][96], many of which are in arid environmental settings.Similar proof for the evaporation paradox was reported in China's northwestern arid regions [33,[97][98][99].Set against a global background of the evaporation paradox, our research results-i.e., in situ measurements reliably estimated using our proposed model-indicate an increase of actual evaporation and decrease of observed pan evaporation in accordance with the available evidence; furthermore, it adds concrete proof of the existence of the evaporation paradox, as well as to the adaptability of the CR theory in China's hyper-arid region.

Conclusions
Quantifying lake evaporation is a theoretical and practical issue, especially in a hyper-arid region, as it offers important and reliable reference for relevant research.Through fine and fundamental evaporation measurement on East Juyan Lake (EJL) in the growing season of 2013-2015, EJL's measured mean total evaporation was 1183.3 mm, and the maximum, minimum, mean daily lake evaporation were 8.1, 3.7, 6.5 mm, respectively.The measured annual open water evaporation was 1471.3 mm; along with E601 pan evaporation and the calculated reference evapotranspiration, these measurements well represent the lower Heihe River's ET 0 with little differences.The measurement of lake evaporation in growing season was solid and reliable, while measured data in non-growing season is lacking, so we converted the data from φ20 pan evaporation at Ejina weather station.Although such a conversion method has unneglectable errors, the results can still be favorable (less than 3% of uncertainty) since non-growing season evaporation only has a small portion (~20%) of annual evaporation; therefore, the measured open water evaporation as potential evaporation (ET 0 ) is convincing.Moreover, this is an improvement to finding a more reliable and practical method to quantify lake evaporation (sublimation) during the non-growing season.
Derived from the Dalton model, parameterized and validated by in situ data, in addition to being compared with acknowledged evaporation models, the proposed daily lake evaporation model (Equation ( 21)) has the advantages not only of simplicity, reliability and accuracy, but also of potential utilization for other lakes in arid and data-sparse regions.Furthermore, apart from the governing synoptic factors considered in this study, lake area, water depth and water quality are also important factors in lake evaporation [8,100,101] that could be used for modifying the proposed model for wider application in arid regions worldwide.
Lastly, we utilized the model to obtain EJL's evaporation during 2002-2015; the mean growing season evaporation was 1233.7 mm (4% larger than measured value).The increasing trends of air temperature, actual lake evaporation and E601 pan evaporation proved the evaporation paradox's existence in this hyper-arid region and supported the complementary relationship (CR) theory's adaptability in this region.In conclusion, the results of this study gave us reliable ET 0 reference for related researches in hyper-arid regions and deepened our understanding of the regional hydrological cycle and its response towards climate change.

Figure 1 .
Figure 1.Location of the Ejina Delta and East Juyan Lake.

Figure 1 .
Figure 1.Location of the Ejina Delta and East Juyan Lake.

Figure 3 .
Figure 3.Time series of measured daily evaporation of EJL.

Figure 3 .
Figure 3.Time series of measured daily evaporation of EJL.

Figure 4 .
Figure 4. Calibration (a) and validation (b) of the daily evaporation model.

Figure 4 .
Figure 4. Calibration (a) and validation (b) of the daily evaporation model.

Figure 5 .
Figure 5. Comparisons of results by the three methods and measurement.

Figure 5 .
Figure 5. Comparisons of results by the three methods and measurement.

Table 1 .
Summary of relevant literature on open water evaporation studies based on in situ measurement or model estimation (the current paper is added for completeness).Keijman, Priestley-Taylor and Penman equation estimated monthly evaporation that agreed most closely with energy budget values.2. Input data measured near the lake is essential to apply these methods, and if only distant data are available, Jensen-Haise and Makkink methods are preferable.
1. Penman equation is simplified for routine hydrologic applications and a more simplified version without wind speed is given and tested.2.The new open water evaporation formulas were adapted for calculating reference crop ET.Rosenberry et al. (2007)Mirror Lake in New Hampshire, USA; 0.15 km 2 ; annual precipitation and evaporation is 1220 and 490 mm. 15 evaporation methods are applied; BREB method is chosen as the standard; In situ measured data are used.1.Priestley-Taylor, deBruin-Keijman and Penman method values compared most favorably with BREB's (within 20% difference).2. Methods may be unduly sensitive to wind speed, when wind speed is neither substantially larger nor smaller.3. Temperature-only methods compared remarkably well with BREB values.4. Results of this research also serve well in other lakes with a similar physical and climatic setting.Masoner et al. (2008)

Table 1 .
Cont.Various modelling methods were used to estimate hourly latent heat fluxes and the theoretical mass transfer model performed the best, the Granger and Hedstrom the worst.2. Estimates by the DYRESM model tended to be greater than measured values.

Table 2 .
Selected evaporation method for comparison with the proposed model.

Table 5 .
Inter-comparison of statistical indexes of evaporation methods.

Table 6 .
Correlativity between synoptic data at Ejina station and EJL station.
Notes: In each expression, X stands for data measured at Ejina weather station and Y stands for EJL bankside data.

Table 5 .
Inter-comparison of statistical indexes of evaporation methods.

Table 6 .
Correlativity between synoptic data at Ejina station and EJL station.
Notes: In each expression, X stands for data measured at Ejina weather station and Y stands for EJL bankside data.

Table 7 .
Different method to represent lower Heihe River's annual ET 0 .