Throughout Africa, good quality water is being recognized as a finite resource that is in short supply-limiting development in many cases  . The Blue Nile basin is one of areas in Africa of extreme water scarcity: nearly 150 million people are dependent on the Blue Nile river with a discharge of less than 85 billion m3   . Understandably, careful assessment of the use of water needs to be made because in the future, it will be increasingly important to match the supply (i.e., excess rainfall) in the Ethiopian highlands with the demand by irrigation systems, industry, and population in the basin. Estimating reference evaporation provides fundamental information on water abstraction in the highlands and downstream when used in combination with crop coefficients  .
Estimating reference evaporation (ET0) has taken many forms ranging from direct measurement to indirect methods employing meteorological measurements. One of the direct measurements is the Class A evaporation pan. The pan evaporation can be related to reference evaporation with coefficients developed by Pereira et al.  and Allen et al.  . The disadvantage of Class A evaporation pans are only their availability at a few locations in developing countries. Indirect methods are, therefore, a good alternative and range from locally developed, empirical relationships to physically based energy- and mass-transfer models  . The meteorological data required vary with the type of indirect methods. The Penman and Priestley and Taylor methods require most meteorological data and the Penman is generally considered as the most accurate one      , but as for the Class A pan data, the required data is only collected at few meteorological stations. Thornthwaite Mather, which is used frequently in Ethiopia, needs fewer climatic data (temperature and sunshine hour). The Hargraves, Blaney-Criddle, and Enku (locally developed) need only temperature to estimate the reference evaporation. Finally, methods that use the measurements of a Piche evaporimeter to calculate the reference evaporation were developed by  and  .
Reference evaporation is defined by Meyer  as the rate of evaporation from a hypothetical crop with an assumed crop height of 12 cm and fixed canopy resistance of 70 s m−1 and an albedo of 0.23 and would give the same evaporation rate as envisioned by Smith et al.  of an extensive surface of green actively growing completely shading the ground and not short of water. It is of course difficult to replicate the conditions that define the reference evaporation. Several methods have been used to validate the direct and indirect methods. In the literature, these validation methods are divided in four categories. The first category uses the evaporation measured with Class A pan    . Other studies employ the ET0 of FAO-Penman Monteith for validation        . Limited experimental studies have used the ET0 calculated from Piche data for validation  . The final method used for validation of calculated reference evaporation are lysimeters     but these are expensive to construct and operate and therefore not available in most developing countries.
Agronomists and hydrologists have little guidance   to choose among the many available methods to estimate reference evaporation in the Ethiopian highlands where the density of weather station is low   and often only measure temperature and precipitation. The objective of this study is, therefore, to select the most accurate method for estimating reference evaporation for the Ethiopian highlands. The Bahir Dar meteorological station was chosen for this study because as the only station in the Ethiopian highlands, all the meteorological variables were measured required for calculating the reference evaporation for all the methods considered. The data recorded on a daily basis included maximum and minimum temperature, relative humidity, wind speed, and sunshine hours and evaporation measured with the Piche and Class A pan.
2. Material and Methods
2.1. Site Description and Data Set
The Ethiopian highlands is the region in Ethiopia above 1500 m and covers a total area of 537,000 km2 (43% of Ethiopia)   . It constitutes more than half of all the highland areas of Africa  . Most of the area of the highlands is cultivated land and land degradation is a major threat. One of the main rivers is the Blue Nile with a watershed of 180,000 km2. It has a monsoonal climate with rainfall varying between 800 and 3000 mm per year and evapotranspiration between 1400 and 1681 per year  . The major rain phase for the highlands of Ethiopia is between June and September.
One of the difficulties in the Ethiopian highlands is that pan evaporation data is measured only at the Bahir Dar Station for long period. It is located 2.2 km from Lake with latitude of 11˚35'59" and longitude of 37˚21'36". The elevation is 1805 m (Figure S1 in the Supplementary Material). The station has more than 30 years of data but Pan Evaporation started only in 2005 (Table S1 in the Supplementary Material list the collected data). Daily maximum and minimum temperature, relative humidity, sunshine hour, wind speed at 2 meter, pan and Piche evaporation data were obtained for the 11 years from 2005 to 2015 from Bahir Dar Meteorology Directorate.
2.2. Reference Evaporation: Description of the Selected Methods
FAO-Penman Monteith, Priestley & Taylor, Hargraves, Thonthwaite, Blaney-Criddle, Enku, Pan-Allen, Pan-Pereira, Piche-Stanhill and Piche-Adam & Ahmed methods were selected. The first four are used locally by experts and by researchers for hydrological modeling and in computations of irrigation water requirement. Piche-Stanhill and Piche-Adam & Ahmed methods were included because of availability Piche readings in several weather stations in the Ethiopian highlands. Pan-Allen and Pan-Pereira methods were used for validation of the rest of the empirical models. Since the Penman has not been validated for the Ethiopian highlands, the Class A pan data using coefficients proposed by Pereira et al.  and Allen et al.  are used as the direct measure of reference evaporation against which the calculated values are compared.
1) FAO-Penman Monteith Method (FAO-PEN)
FAO-Penman Monteith method computes reference evapotranspiration (ET0) from meteorological data. It has been recommended as the standard method for computation of the reference evapotranspiration. The method considers all parameters that govern energy exchange and corresponding latent heat flux and requires air temperature, relative humidity, sunshine intensity/hour, wind speed data and elevation as input parameters.
The modified Penman-Monteith equation can be written as  :
where ET0 reference evapotranspiration [mm day−1], Rn net radiation at the crop surface [MJ m−2 day−1], G soil heat flux density [MJ m−2 day−1], T mean daily air temperature at 2 m height [˚C], u2 wind speed at 2 m height [m s−1], es saturation vapor pressure [kPa], ea actual vapor pressure [kPa], es − ea saturation vapor pressure deficit [kPa], Δ slope vapor pressure curve [kPa ˚C−1], psychrometric constant [kPa ˚C−1]. Further details how these parameters are calculated are given in the Supplementary Material Section.
2) Priestley & Taylor Method (PT)
The Priestley & Taylor form of the equation is    :
where; rs is surface resistance [s m−1] which is 70 s m−1 for the grass reference surface, ra is aerodynamic resistance which helps to determine the transfer of heat and water vapor from the evaporating surface into the air above the canopy. The aerodynamic resistance ra [s m−1] for the grass reference surface becomes:
where u2 is the wind speed [m s−1] at 2 m.
3) Hargraves Method (HAR)
This method requires daily maximum and minimum air temperature, and extraterrestrial solar radiation, which is computed from the latitude of the study site. The Hargreaves and Samani  equation is defined as follows:
where Tmin daily minimum temperature [˚C], Tmax saturation vapor pressure at daily maximum temperature and Ra extraterrestrial radiation [MJ m−2 day−1].
4) Thornthwaite Method (TH)
This formula is based mainly on temperature with an adjustment factor for the number of daylight hours. The potential evapotranspiration, calculated on a monthly basis, is  :
where the subscript m indicates the months 1, 2, 3, …, 12, Tm is the monthly mean temperature (˚C), I is the heat index for the year, given by:
where, im is monthly values of heat index and a is empirical coefficient defined as
where, Nm is the monthly adjustment factor related to hours of daylight, ET0 reference evaporation [mm/day], d duration of daylight [hr].
5) Blaney-Criddle Method (BC)
The Blaney-Criddle equation has been used to calculate the reference crop evapotranspiration ET0 when pan evaporation is not available. It can be expressed as 
where ET0 is reference crop evapotranspiration (mm/day) averaged over the month Tmean = mean daily temperature (˚C), and p is mean daily percentage of annual daytime hours and varies between 0.26 and 0.29 for Bahir Dar.
6) Enku’s Simple Temperature Method (ENKU)
The new simple empirical temperature method, which was named by “Enku’s simple temperature method”, was developed  and tested in Ethiopia with Penman Montieth reference Evaporation and the Piche Evaporimeter. The equation is:
where ET0 is the reference evapotranspiration (mm day−1); n = 2.5 for the Lake Tana area; k is the coefficient, which is calibrated for local conditions ranging from about 600 for lower mean annual maximum temperature areas to 1300 for higher mean annual maximum temperature areas. The coefficient, k, was found by Enku and Melesse  as for combined wet and dry conditions or for the dry phase, and for the rain phase, where Tmm (C) is the long term daily mean maximum temperature for the seasons under consideration.
7) Pan Methods
The evaporation rate from Class A Pans filled with water is easily obtained. In the absence of rain, the amount of water evaporated during a period (mm/day) corresponds with the decrease in water depth. Pans provide a measurement of the integrated effect of radiation, wind, temperature and humidity on the evaporation from an open water surface  .
There are two types of Pan Models to estimate reference evaporation.
a) Pereira Model (PAN-P)
According to Pereira et al.  , the reference evaporation, ET0 from pan data calculated with:
where; ET0 reference evapotranspiration [mm/day], K1 pan coefficients [-], EPan Class A pan evaporation [mm/day].
b) Allen Model (PAN-A)
As Allen et al.  , ET0 from pan data calculated with:
The pan coefficient K2 for Class A pan with green fetch can be found as  :
where F is fetch or distance of the reference grass (m) and RH is the relative humidity in percent.
8) Piche Methods
The Piche evaporimeter is a type of atmometer used to measure the rate of evaporation from a wet disc of absorbent paper. It is used mainly in hot, dry climates where water loss through evaporation must be observed regularly. Since the results are dependent on wind speed past the disc, as well as the wet bulb saturation deficit, it is almost essential to expose the evaporimeter inside a meteorological screen    .
There are two types of Piche models used in this study to estimate evapotranspiration.
a) Stanhill Method (PI-S)
Stanhill  suggested that it may be possible to estimate the second term in Penman’s equation from available sheltered Piche evaporation data (Epi) as follows.
where a slope and b is intercept of a linear relation with the aerodynamic term of Penman equation:
where f(u) is an aerodynamic wind function and (es − ea) is the difference between saturated vapour pressure and actual vapour pressure in hPa, evaluated at mean air temperature and at 2 m above the ground or water surface.
The equivalent to the wind function of Penman  is
where u is the wind speed in m s−1 at 2 m elevation, and a, and b, are empirical coefficients. Penman  , Penman  suggested using values equivalent to 1 and 0.537 for a, and b, respectively for a short grass cover when u is measured in m s−1.
b) Adam & Ahmed Method (PI-ADAH)
Adam and Ahmed  showed that the ratio of Penman estimated and Piche evaporation have relation with relative humidity exponentially. That is:
where a and b are constants that can be find from the exponential relation and the RH humidity in percent
2.3. Methods of Data Analysis
The collected meteorological data for Bahir Dar (2005-2015) was of good quality with less than 5% of missing data. A simple arithmetic mean was used to determine the missing data. Reference evaporation was computed at daily and monthly time step using the 10 methods listed above after data preparation.
In addition, we compared the calculated reference evaporation with the ten methods using three techniques: visual inspection, descriptive statistics and statistical methods to test methods of efficiency and their standard error. The visual inspection of plotted ET0 prediction methods reveals whether the calculated reference evaporation is in agreement with either other methods or an outlier. Descriptive statistics like maximum, absolute minimum, median value, total range of values, standard deviation, and coefficient of variation were used to compare the ET0 prediction methods.
In this study, every single model was correlated with the pan ET0s to assess the model performances and standard error. Pearson’s correlation and root mean square error (RMSE) equations were used for this purpose.
The formula for the Pearson product moment correlation coefficient, r, is:
The equation for the standard error of the predicted y is:
where and are the sample means.
A sensitivity analysis was carried out to evaluate the percentage response of calculated ET0 to selected weather variables. Changes of model outputs and their variability induced by change in weather variables were evaluated. The change in variables was analyzed for values 25% above and 25% below the mean value.
Mann Kendall trend test was used to indicate whether there are trends in the ET0 computed using the ten methods. Non-parametric Mann Kendall Trend Test is useful to examine the temporal variation trend. It is based on the significance of differences, not directly on the random values. Therefore the trend that’s been determined is less affected by outliers  . The nonparametric Mann-Kendall trend test has been applied in many studies to identify whether monotonic trends exist in hydro-meteorological data such as temperature, rainfall and stream flow  . The Mann Kendall Trend Test, S is calculated by using the equation below   .
where xj and xi are the sequential data value and j greater than i, N is the length of the data set.
As indicated in Mann  and Kendall and Stuart  , when N ≥ 8, the distribution of S approaches the Gaussian form with mean E(S) = 0 and variance Var(S) given by:
where: ti is the number of ties of length m.
The statistic S is then standardized (Z), and its significance can be estimated from the normal cumulative distribution function.
The positive Z value indicates an increasing trend while a negative Z value indicates a decreasing trend. When testing two sided trends at a selected level of significance a, the null hypothesis (H0) of no trend is rejected if the absolute value of Z is greater than Za/2 where a represents the chosen significance level (5% with Z0.025 = 1.96).
Reference Evaporation: Calculated Values
The reference evaporation, ET0, was calculated for a daily and monthly time step. In Figure 1, as an example, the 2005 ET0 for the ten methods are depicted. Values calculated for the other years are shown in tabular form in Table S2 in the supplementary material. The minimum reference evaporation calculated was 0 mm/day (Thornthwaite) and 16.3 mm/day (Hargraves and Priestley-Taylor). Blaney-Criddle has the smallest range (the difference in daily maximum and minimum values) while Priestley-Taylor has the largest range of values. The FAO-Penman and Enku methods have acceptable daily ET0 ranges between 1.5 and 6.5 mm/day    .
Seasonally, the reference evaporation of each of the methods is greater during the dry phase from February to May when the temperature is high and there are few clouds than during the rain phase from June to September when it is cloudy and less warm (Figure 1, Figure 2). In addition, the reference evaporation decreases in December and January for some of the methods (Figure 1, Figure 2) due to cloudiness caused by of easterlies carrying moisture from Arabian Sea.
In equatorial countries like Ethiopia, the difference in sunshine duration between months is minimal and the difference in ET0 is due to cloud cover mainly and temperature and relatively humidity secondly. This is different from temperate climates where the day length is the primary factor that determines the magnitude of the reference evaporation during the year.
The methods by Enku, Class A pan (both), FAO-Pennman Monteith and the Adam and Ahmed Piche evaporimeter had all the same annual average reference evaporation of 1460 mm y−1 (Figure 3). The predictions of the Blaney-Criddle and Thornthwaite are slightly more at 1754 mm y−1 (Figure 3). The last three
Figure 1. Daily plot of reference potential evaporation (ET0) of the ten methods for 2005 for the Bahir Dar meteorological station. The following abbreviations are used: FAO-PEN is FAO-Penman Monteith (Equation (1)), PT is Priestley & Taylor (Equation (2)), HAR is Hargraves (Equation (5)), TH is Thornthwaite (Equation (9)), BC is Blaney-Criddle (Equation (10)), ENKU is the Enku’s simple temperature method (Equation (11)), PAN-P is the Pan-Pereira (Equation (12)), PAN-A is the Pan-Allen (Equation (14)), PI-S is the Piche-Stanhill (Equation (16)) and PI-ADAH is the Piche-Adam and Ahmed (Equation (20)).
Figure 2. Monthly plot of ET0 of the ten methods for the Bahir Dar meteorological station from 2005 to 2015. Abbreviations of the reference evaporation methods are listed under Figure 1.
Figure 3. Annual reference evaporation, ET0, of the ten methods for the Bahir Dar meteorological station from 2005 to 2015. Abbreviations of the reference evaporation methods are listed under Figure 1.
methods also over predicted pan evaporation in the Ethiopian highlands: the Stanhill Piche evaporimeter method predicted on the average 1571 mm y−1, Priestley & Taylor methods as 1838 mm y−1 and Hargraves almost twice the pan evaporation ET0 at 2792 mm y−1. In addition, Figure 3 showed that the year to year variation in reference evaporation was small, indicating that the amount of rainfall and soil moisture have little or no effect on the loss of water in the atmosphere.
4.1. Reference Evaporation: Comparing the Methods
Figure 4. Long-term mean monthly evapotranspiration (2005-2015) of Bahir Dar for ten different prediction techniques. Abbreviations of the reference evaporation methods are explained Figure 1.
(2005-2015) are depicted. As in Figure 2 the Hargraves (HR), the Priestley & Taylor (PT) and the Piche-Stanhill method (PI-S) over predicts the reference evaporation greatly and will not be discussed in any detail. These three methods should not be applied in the Ethiopian highlands to calculate evaporation.
While all methods have the smallest reference evaporation during the rain phase (Figure 4), the Blaney-Criddle (blue solid line) is the exception. At the end of dry phase, the ET0 calculated with the Blaney-Criddle method is in the same order as other indirect measurement that predicts realistically ET0 such as the FAO-Penman and the Enku methods, but severely over predicts during the rain phase. It has also the smallest variability compared to all other methods (Figure 4). The reason is that the reference evaporation is calculated as linear function of the average monthly temperature (Equation (10)) which in the countries near the equator varies little during the year and thus the ET0 values are within a narrow band as well (Figure 4, Figure 5). By replacing the constant values for the whole year in the Blaney-Criddle equation (Equation (10)), by monthly varying constants and calibrating these, the model will fit much better  . The modified Blaney-Criddle Equation becomes in this way similar to the Enku method (Equation (11)). Interestingly, in Figure 4 has shown that during the rain phase starting in May (when the first rains fell) through September (when the rains ended), the reference evaporation, ET0, calculated with the three indirect methods (the FAO-Penman Monteith, the Thornthwaite and the Enku methods) compared well with the direct measurement of both pan evaporation methods and the Piche evaporimeter using coefficients proposed by Adam and Ahmed  .
During the rain phase when precipitation exceeds evaporation, the condition of a well-watered surface on which the reference evaporation is based is similar to that in the Ethiopian highlands. So, during the wet phase the direct and indirect
Figure 5. Box plots of daily ET0 prediction methods for the Bahir Dar weather station. Abbreviations of the reference evaporation methods are explained Figure 1.
measurements should agree and they do for the FAO-Penman Monteith, the Thornthwaite and the Enku methods with the two pan methods and one of the Piche methods (Figure 4). During the dry phase, the reference evaporation for these seven methods started to deviate. Initially after the rains stopped and the surface of soil was still wet, only the Thorntwaite method predicted higher ET0 values while the remaining five agreed well. In January to April when the landscape is very dry and rivers are without water, the direct measurements of the ET0 are much greater than the calculated values using the meteorological data (Figure 4). The reason for the deviation can be explained with the complementary relationship of actual evaporation and apparent (or pan) evaporation originally introduced by Bouchet  and refined by Brutsaert and Stricker  . The theory is as follows: When the soil is dry and the evaporation is less than the reference evaporation, the energy that under wet conditions is used for evaporation is converted to the sensible heat and causes an increase in the evaporative demand of the atmosphere. The reference evaporation rates derived from the pan during the dry season are therefore greater than it would be during the wet season under otherwise similar atmospheric conditions. The FAO-Penman method predicts reference evaporation rates of a well-watered irrigation field independent of the condition of the landscape. Hence, the discrepancy between the direct and indirect methods under extreme dry conditions in the period is from January through April (Figure 4).
The question whether the direct or indirect measurement of the reference evaporation is more appropriate during the dry phase is a mood point because the evaporation for most of the landscape, where the soil is dry, is limited by the soil and not by the atmosphere. Only for irrigated fields and lakes, the reference evaporation determines the rate of evaporation during the dry phase. As pointed out by Bouchet  the upward wind site of an irrigated field, the evaporation rate is similar to the pan measurement while in the remaining of the area the evaporation rate is equal to reference evaporation calculated with the Penman Monteith.
In addition to plotting the averaged reference evaporation, ET0, in Figure 5, it is also of interesting to investigate the variation in the reference evaporation for each of the 10 methods. Therefore, the maximum, minimum, mean, lower and upper quartile of the long term daily reference evaporation for the period are shown in Figure 5 from 2005 to the end of 2015. Figure 5 indicates that the Blaney-Criddle method has the least variation because it is only dependent of the average temperature that (despite what the Ethiopians claim) varies little throughout the year. The FAO-Penman Monteith has only a slightly larger variation in ET0 than the Blaney-Criddle method because it calculates likewise the ET0 based on meteorological variables that vary little throughout the year except for the cloud cover which is the reason of the additional spread in ET0 values. The Enku method that depends on temperature has a larger variation than the Blaney-Criddle, because besides temperature it is dependent on a few other fitted functions so that it can simulate the lower ET0 during the month with clouds. The two pan measurements have a large variation, because they are based on measured data. The measured ET0 values have a large spread because of the sensible heat during the dry phase that increases the pan evaporation but is not included in any of the indirect methods.
In the Pearson correlation statistics, we looked the performance of eight reference evaporation methods with the pan evaporation methods (Table 1). As expected form the discussion above, the FAO-Penman Monteith and Enku simple maximum temperature methods were highly correlated with correlation coefficient of respectively 0.91 and 0.93 (monthly), and 0.64 and 0.69 (daily). Thornthwaite method holds the third place of good ET0 monthly estimator with correlation of 0.89 for both pan ET0 methods on a monthly basis.
Rácz et al.  stated that methods having smallest RMSE have the lowest systematic error in predicting ET0. Table 2 confirms our findings above that the Penman, Enku, Piche and Pan are most similar because they have the smallest RMSE (Table 2). It is surprising that the Blaney-Criddle model has such a small RMSE with the pan on daily time step. It is quirk of how the daily ET0 values are calculated.
Table 1. Pearson’s correlation between Pan ET0 and ET0 of prediction methods (2005- 2015).
Table 2. Root mean square error (RMSE) between modeled daily ET0 and pan ET0 (2005- 2015).
Figure 6. Relative change in reference potential evaporation of Bahir Dar station for the change of temperature between the year 2005 and 2015. Abbreviations of the reference evaporation methods are explained Figure 1.
4.2. Sensitivity Analysis
In the sensitivity analysis, we looked the relative effect on reference evaporation of a 25% change in temperature, sunshine hours, relative humidity and wind speed (Figure 6, Figure S2). As expected, an increase in temperature increased reference evaporation (Figure 6). For the indirect methods (Enku, Thornthwaite and Hargreaves methods), the reference evaporation increased the same percentage or more as the temperature. The direct methods were the least sensitive to the change in temperature. Changing sunshine hours, relative humidity and wind speed did not change or increased the reference evaporation with the exception that an increase in relatively humidity for the FAO-Penman Monteith and an increase in winds Speed for both pan methods decreased the reference evaporation (Figure S2).
4.3. Trend Analysis
The trend analysis results of the reference potential evaporation obtained by applying the Mann-Kendall method is shown in Table 3. To our surprise despite our short 11-year record, there were strong trends for the majority of the 10
Table 3. Mann Kendall trend test statistics for daily ET0 of Bahir Dar station from 2005 to 2015.
Figure 7. Time series plots of ET0 for Bahir Dar station from 2005 to 2015: (a) FAO-Penman Monteith; (b) Pan-Pereira.
methods. Only the direct measurement of the ET0 with the two pan methods did not change over the eleven year period (i.e., accept the null hypothesis that there is not a trend in the ET0 values, Table 3, Figure 7(b)). Methods that calculated the reference evaporation from meteorological data (e.g., FAO-Penman Monteith (FAO-PEN), Hargraves (HAR), Thornthwaite (TH), Blaney-Criddle (PT), Enku (ENKU)) increased with time (Table 3, Figure 7(a) and Figure S3). The Priestly and Taylor methods that was ill suited for the Ethiopian highlands was an exception. Thus, since these methods were directly based on measured meteorological data (two only on temperature) it clearly indicates that the climate is changing in the Ethiopian Highlands. In other words, since three of the indirect methods were primarily based on temperature, it is getting warmer fast in Bahir Dar! Finally, the Piche methods that not directly depend on the measured standard indicated that the reference evaporation was decreasing.
This paradox of increasing reference evaporation calculated by indirect methods and either not changing (both pan measurements) or decreasing (Piche measurement) by direct measurement has been noted before in a slightly different context by Brutsaert and Parlange  as one of the first. They noted that pan evaporation rates were generally decreasing and related that to increased terrestrial evaporation because more runoff occurring with climate change. In terms of the complementary relationship introduced earlier, they explain that more of the incoming energy is used for evaporation and consequently less is converted to sensible heat. This in turn reduces pan evaporation. For the Bahir Dar meteorological station, the explanation might be slightly different. Temperatures in Bahir Dar are increasing either due to climate change or because of rapid urbanization. The dry season has been always so dry that all evaporative energy was converted to sensible heat and increasing temperatures did not affect pan evaporation. The Piche evaporimeter measurements were likely not affected either by the higher temperatures but the relatively humidity decreased due to the higher temperature decreasing the ET0 value as can be seen from Equation (20).
Ten methods to predict the reference evaporation were tested for the Ethiopian highlands. The Priestley and Taylor, Hargreaves, and Piche-Stanhill methods over predicted ET0 and should not be used without recalibration. In addition, the Blaney-Criddle method over predicted the reference evaporation during the rain phase. The reference evaporation calculated with the FAO-Penman and Enku methods resembles most closely the direct measurements with the Class A pan evaporation using confident introduced by Pareira and Allen. Thorthwaite’s monthly ET0 model performed well too.
A significant increasing trend in calculated reference evaporation using meteorological variables was found, indicating that temperatures were increasing during the past 11 years at the Bahir Dar weather station. The direct measurement by the Class A pan did not show this trend. More research is needed to research whether other parts of the Ethiopian highlands show similar trends in reference evaporation.
Our recommendation is that the FAO-Penman Monteith is recommended for locations where the input data are available. Otherwise, the Enku method using maximum daily temperature is best for estimating the reference evaporation.
Authors gratefully acknowledge the Ethiopian National Meteorological Agency (NMA)―Bahir Dar Meteorology Directorate for providing the meteorological datasets of Bahir Dar principal station. The frank discussion on evaporation methods with Professor Brutsaert was very helpful.
Supplementary Material S1: Parameters of FAO-Penman Monteith Equation
The slope of the relationship between saturation vapor pressure and temperature, Δ, is
where, T air temperature [˚C]. The actual vapor pressure can be determined from the difference between the dry and wet bulb temperatures, ea
where ea actual vapor pressure [kPa], saturation vapor pressure at daily minimum temperature [kPa], saturation vapor pressure at daily maximum temperature [kPa], RHmean is the mean relative humidity, defined as the average between RHmax and RHmin.
As saturation vapor pressure is related to air temperature, it can be calculated from the air temperature. The relationship is expressed by:
where saturation vapour pressure at the air temperature T [kPa], T air
Figure S1. Location of Bahir Dar principal meteorological station.
Table S1. Data used and overview of ET0 methods calculation.
temperature [˚C], exp[..] 2.7183 (base of natural logarithm) raised to the power [..].
Saturation vapor pressure as a function of air temperature, es
Soil heat flux (G):
where Tmonth,i mean air temperature of month i [˚C] and Tmonth,i−1 mean air temperature of previous month [˚C].
The psychrometric constant, , is given by:
where psychrometric constant [kPa ˚C−1], P atmospheric pressure [kPa], λ latent heat of vaporization, 2.45 [MJ kg−1], cp specific heat at constant pressure, 1.013 ´ 10−3 [MJ kg−1 ˚C−1], ratio molecular weight of water vapor/dry air = 0.622.
The atmospheric pressure, P, is the pressure exerted by the weight of the earth’s atmosphere. Evaporation at high altitudes is promoted due to low atmospheric pressure as expressed in the psychrometric constant. The effect is, however, small and in the calculation procedures, the average value for a location is sufficient. A simplification of the ideal gas law, assuming 20˚C for a standard atmosphere, can be employed to calculate P:
where P atmospheric pressure [kPa], and Z elevation above sea level [m]
The net radiation (Rn) is the difference between the incoming net shortwave radiation (Rns) and the outgoing net long wave radiation (Rnl):
where Rns net solar or shortwave radiation [MJ m−2 day−1], α albedo or canopy reflection coefficient, which is 0.23 for the hypothetical grass reference crop [dimensionless], Rs the incoming solar radiation [MJ m−2 day−1].
If the solar radiation, Rs, is not measured, it can be calculated with the Angstrom formula which relates solar radiation to extraterrestrial radiation and relative sunshine duration:
where Rs solar or shortwave radiation [MJ m−2 day−1], n actual duration of sunshine [hour], N maximum possible duration of sunshine or daylight hours [hour], n/N relative sunshine duration [-], Ra extraterrestrial radiation [MJ m−2 day−1], as regression constant, expressing the fraction of extraterrestrial radiation reaching the earth on overcast days (n = 0), as + bs fraction of extraterrestrial radiation reaching the earth on clear days (n = N). If no actual solar radiation data are available and no calibration has been carried out for improved as and bs parameters, the values as = 0.25 and bs = 0.50 are recommended.
The daylight hours, N, are given by:
where ωs is the sunset hour angle in radians
where Ra extraterrestrial radiation in the hour (or shorter) period [MJ m−2 hour−1], Gsc solar constant = 0.0820 MJ m−2 min−1, dr inverse relative distance Earth-Sun, δ solar declination [rad], latitude [rad], ωs sunset hour angle [rad].
The inverse relative distance Earth-Sun, dr, and the solar declination, δ, are given by:
Table S2. Results overview of potential evaporation methods between 2005 and 2015 for Bahir Dar meteorological station. Abbreviations of the reference evaporation methods are explained Figure 1.
where; J is the number of the day in the year between 1 (1 January) and 365 or 366 (31 December). The sunset hour angle, ωs, is given by:
Net long wave radiation (Rnl):
where Rnl net outgoing longwave radiation [MJ m−2 day−1], σ Stefan-Boltzmann constant [4.903 ´ 10−9 MJ K−4 m−2 day−1], Tmax,K maximum absolute temperature during the 24-hour period [K = ˚C + 273.16], Tmin,K minimum absolute temperature during the 24-hour period [K = ˚C + 273.16], ea actual vapour pressure [kPa], Rs/Rso relative shortwave radiation (limited to ≤ 1.0), Rs measured or calculated solar radiation [MJ m−2 day−1], Rso calculated clear-sky radiation [MJ m−2 day−1].
where; Z station elevation above sea level [m].
Figure S2. Relative change in reference potential evaporation of Bahir Dar station for the change of sunshine hour (a), relative humidity (b) and wind speed (c) between the year 2005 and 2015. Abbreviations of the reference evaporation methods are explained Figure 1.
Figure S3. Daily time series plots of ET0 for (a) Priestley & Taylor; (b) Hargraves; (c) Blaney-Criddle; (d) Enku; (e) Pan-Allen; (f) Piche-Stanhill; and (g) Piche-Adam and Ahmed.
 Smith, M., Segeren, A., Santos Pereira, L., Perrier, A. and Allen, R. (1991) Report on the Expert Consultation on Procedures for Revision of FAO Guidelines for Prediction of Crop Water Requirements. Rome, 28-31 May 1990.
 Djaman, K., Tabari, H., Balde, A.B., Diop, L., Futakuchi, K. and Irmak, S. (2016) Analyses, Calibration and Validation of Evapotranspiration Models to Predict Grass-Reference Evapotranspiration in the Senegal River Delta. Journal of Hydrology: Regional Studies, 8, 82-94.
 Gotardo, J.T., Rodrigues, L.N. and Gomes, B.M. (2016) Comparison of Methods for Estimating Reference Evapotranspiration: An Approach to the Management of Water Resources within an Experimental Basin in the Brazilian Cerrado. Engenharia Agrícola, 36, 1016-1026.
 Pandey, P.K., Dabral, P.P. and Pandey, V. (2016) Evaluation of Reference Evapotranspiration Methods for the Northeastern Region of India. International Soil and Water Conservation Research, 4, 52-63.
 Rácz, C., Nagy, J. and Dobos, A.C. (2013) Comparison of Several Methods for Calculation of Reference Evapotranspiration. Acta Silvatica et Lignaria Hungarica, 9, 9-24.
 Hashem, A., Engel, B., Bralts, V., Radwan, S. and Rashad, M. (2016) Performance Evaluation and Development of Daily Reference Evapotranspiration Model. Irrigation & Drainage Systems Engineering, 5, 1-6.
 Bogawski, P. and Bednorz, E. (2014) Comparison and Validation of Selected Evapotranspiration Models for Conditions in Poland (Central Europe). Water Resources Management, 28, 5021-5038.
 Djaman, K., Balde, A.B., Sow, A., Muller, B., Irmak, S., N’diaye, M.K., Manneh, B., Moukoumbi, Y.D., Futakuchi, K. and Saito, K. (2015) Evaluation of Sixteen Reference Evapotranspiration Methods under Sahelian Conditions in the Senegal River Valley. Journal of Hydrology: Regional Studies, 3, 139-159.
 Čadro, S., Uzunović, M., Žurovec, J. and Žurovec, O. (2017) Validation and Calibration of Various Reference Evapotranspiration Alternative Methods under Bosnia and Herzegovina Climate Conditions. International Soil and Water Conservation Research.
 Xing, Z., Chow, L., Meng, F.-R., Rees, H.W., Steve, L. and Monteith, J. (2008) Validating Evapotranspiraiton Equations using Bowen Ratio in New Brunswick, Maritime, Canada. Sensors, 8, 412-428. https://doi.org/10.3390/s8010412
 Douglas, E.M., Jacobs, J.M., Sumner, D.M. and Ray, R.L. (2009) A Comparison of Models for Estimating Potential Evapotranspiration for Florida Land Cover Types. Journal of Hydrology, 373, 366-376.
 Razzaghi, F. and Sepaskhah, A.R. (2010) Assessment of Nine Different Equations for ETo Estimation using Lysimeter Data in a Semi-Arid Environment. Archives of Agronomy and Soil Science, 56, 1-12. https://doi.org/10.1080/03650340902829180
 Razzaghi, F. and Sepaskhah, A.R. (2012) Calibration and Validation of Four Common ET0 Estimation Equations by Lysimeter Data in a Semi-Arid Environment. Archives of Agronomy and Soil Science, 58, 303-319.
 Worqlul, A.W., Ayana, E.K., Maathuis, B.H., Macalister, C., Philpot, W.D., Leyton, J.M.O. and Steenhuis, T.S. (2017) Performance of Bias Corrected MPEG Rainfall Estimate for Rainfall-Runoff Simulation in the Upper Blue Nile Basin, Ethiopia. Journal of Hydrology.
 Worqlul, A.W., Yen, H., Collick, A.S., Tilahun, S.A., Langan, S. and Steenhuis, T.S. (2017) Evaluation of CFSR, TMPA 3B42 and Ground-Based Rainfall Data as Input for Hydrological Models, in Data-Scarce Regions: The Upper Blue Nile Basin, Ethiopia. Catena, 152, 242-251.
 Nyssen, J., Vandenreyken, H., Poesen, J., Moeyersons, J., Deckers, J., Haile, M., Salles, C. and Govers, G. (2005) Rainfall Erosivity and Variability in the Northern Ethiopian Highlands. Journal of Hydrology, 311, 172-187.
 Priestley, C. and Taylor, R. (1972) On the Assessment of Surface Heat Flux and Evaporation using Large-Scale Parameters. Monthly Weather Review, 100, 81-92.
 Penman, H.L. (1948) Natural Evaporation from Open Water, Bare Soil and Grass. Proceedings of the Royal Society of London A: Mathematical, Physical and Engineering Sciences, The Royal Society, 120-145.
 Suhaila, J., Deni, S.M., Zin, W.Z.W. and Jemain, A.A. (2010) Trends in Peninsular Malaysia Rainfall Data during the Southwest Monsoon and Northeast Monsoon Seasons: 1975-2004. Sains Malaysiana, 39, 533-542.
 Sulaiman, N.H., Kamarudin, M.K.A., Mustafa, A.D., Amran, M.A., Azaman, F., Abidin, I.Z. and Hairoma, N. (2015) Trend Analysis of Pahang River Using Non-Parametric Analysis: Mann Kendall’s Trend Test. Malaysian Journal of Analytical Sciences, 19, 1327-1334.