Evaluation of Eleven Reference Evapotranspiration Models in Semiarid Conditions ()
1. Introduction
Mean annual rainfall across the West African Sahel, varies from 100 mm to 600 mm covering 1 to 6 months with the decreasing trends from the southern limit to the northern limit [1] . The extremely high air temperature and the low relative humidity increase the Evapotranspiration demand. Evapotranspiration is the main source of water losses and one of the important parameters under the hydrological, agricultural and environmental studies [2] . Under the Sahelo-Saharan climate conditions, water resource is one of the most factors limiting food production where insufficient and erratic rainfall limits rainfed agriculture in high risk production systems in terms of water stress. The Sahel is characterized by semi-arid climate under which effective water management is, therefore, critical. Accurate crop evapotranspiration estimation is critical for the water resources management under agricultural, hydrological and environmental processes. Management of water resources in the limited available water resources environment like the west African Sahel merit primordial attention for resource sustainability and improving water use efficiency [3] [4] [5] . Crop evapotranspiration is mostly estimated by the indirect method combining crop reference evapotranspiration (ETo) with crop coefficients [6] beside the direct measurement though lysimeters [7] [8] . Different methods for estimating ETo have been developed for different parts of the globe using climatic variables such as temperature, radiation, and combined variables [9] - [14] . While most of developed ETo equations have relative accuracy and adaptability to the local environment different from where they were developed, the Penman-Monteith ETo equation was demonstrated and shown to be the most accurate and adapted to all climatic conditions [12] [15] - [20] . The applicability of the Penman-Monteith ETo equation is constrained by the number and the non-availability of the full climatic dataset (temperature, relative humidity, solar radiation, and wind speed) mostly the developing countries like the sub-Saharan African countries. Scientists have to evaluate the adaptability of simple ETo equations to the local climatic conditions. One of the Valiantzas ETo equations was shown suitable ETo estimation across Burkina Faso [4] [21] , Tanzania, Kenya and Uganda [20] [22] . Djaman et al. [3] reported the Trabert, Mahringer, Penman 1948, Albrecht, and two of Valiantzas’ equations to perform well under the semiarid condition in Senegal River Delta and their calibration to the local climatic condition had improved daily ETo [5] . Tabari et al. [23] reported suitability of very few ETo equation to the Iranian environment while [24] indicated good performance of the calibrated forms of two of the Valiantzas equations to the Pilbara region of Western Australia.
While different equations have been tested and calibrated for different regions and sub-regions under different climatic conditions including in Canada [25] , in Iran [26] [27] [28] , in China [29] , in Poland [30] , in Florida (USA) [31] , in Southeast Australia [32] , and in Senegal [3] , limited studies were conducted on the adaptability of reference evapotranspiration equation to the African Sahelo- Saharan region. Therefore, it is critical to assess the adaptability of ETo models and improve them through proper calibration to the local climate conditions. Thus, the objectives of this study were to: 1) evaluate the FAO-PM ETo equation under limited data condition, 2) evaluate nine temperature reference evapotranspiration equations with comparison to the FAO-PM method, and 3) calibrate and validate the Abtew ETo equation [33] under the Malian semiarid and arid climate conditions.
2. Materials and Methods
2.1. Data and Study Area
Climatic data were collected at eleven weather stations across Mali (West Africa) for the period of 1990-2012. Maximum temperature (Tmax), minimum temperature (Tmin), maximum relative humidity (RHmax), minimum relative humidity (RHmin), wind speed (u2), and solar radiation (Rs) data were collected from Bougouni, Hombori, Koutiala, Mopti, Nara, San, Segou, Senou, Sikasso, and Tessalit (Table 1). Mali has three climatic zones: the Sudanian zone with 700 to 1000 mm of annual precipitation, the Sahelian zone which receives 200 to 400 mm of precipitation, and the Saharan zone with little or no rain. Mali is among the hottest countries in the world. Most of Mali receives negligible rainfall and droughts are very frequent. The rainy season covers generally late June to early December in the southern area. The vast northern desert part of Mali has a hot desert climate with long, extremely hot summers and scarce rainfall which decreases northwards. The central area has a hot semi-arid climate with very high temperatures year-round, a long, intense dry season and a brief, irregular rainy season.
2.2. Reference Evapotranspiration Equations
- Penman-Monteith equation (FAO-PM)
Daily grass reference evapotranspiration (ETo-Ref) was calculated using the Penman-Monteith (FAO-PM) equation [12] . The Penman-Monteith reference evapotranspiration equation for grass surface is:
Table 1. Geographic coordinates and the climatic zones of the eleven weather stations.
(1)
where all variables are defined as described in [22] . All parameters necessary for computing ETo were computed according to the procedure developed in FAO-56 by [12] .
The FAO-PM ETo model was evaluated under limited climatic data conditions for its performance and suitability under similar conditions as it is almost the case in most of the developing countries where climatic data record is not consistent and with missing variables [3] [34] [35] [36] [37] .
Thus, the following acronyms were used for ETo estimated using: 1) FAO-56 Penman Monteith equation with full data set is referred as FAO-PM ETo, 2) ETo-Rs when Rs is missing, 3) ETo-Tmin when RH is missing, 4) ETo-um when u2 is missing, 5) ETo-RsTmin when Rs and RH are missing, 6) ETo-Rsum when Rs and u2 are missing, 7) ETo-umTmin when u2 and RH are missing, 8) ETo-RsumTmin when Rs, RH and u2 are missing.
- Jensen and Haise [6] method:
(2)
- Hansen [38] method:
(3)
- Abtew [33] method 1: Abtew
(4)
- Calibrated Christiansen [33] method:
(5)
- Droogers and Allen [39] method: Dr-Al
(6)
- Hargreaves and Allen [40] method: Harg
(7)
- Irmak [41] method:
(8)
- Tabari [23] method 1: Tabari 1
(9)
- Tabari [23] method 2: Tabari 2
(10)
where, Tmax, Tmin and Tmean are daily maximum, minimum and mean air temperature (˚C), respectively; Rs is solar radiation, Ra is extraterrestrial radiation (MJm−2・day−1), K is dimensionless coefficient to be determined for each location.
- Multi-model ensemble (MME)
Simple multi-model ensemble (MME) was constructed by combining all the nine the individual ensemble of ETo models with equal weights [42] [43] . Multi-model ensemble outperforms single models in their skill due to error cancellation and the nonlinearity of the skill metrics applied. Hagedorn et al. [43] reported that it is not usually possible to identify a “best” or a “poorest” model from a set of models, as their individual strengths and flaws typically vary with location and initialization time. In this method, MME forecasts are generated by simply pooling together the participating simple models, with all ensemble members having equal weight [43] . We assume that this procedure might solve the quantification of all aspects of all simple models uncertainties with increasing model performance as demonstrated by [44] [45] [46] [47] .
2.3. Calibration of the Abtew ETo Equation
To calibrate the ETo equations, a linear regression relationship between daily PM-ETo and daily ETo estimates by the Abtew equation was determined and the calibration coefficients were then obtained by multiplying the slope of a regression line between ETo estimate by an ETo equation and the PM-ETo by its inverse to bring the slope of the regression line to the unity. And, the opposite value of the intercept was added to the new regression relationship to minimize the new intercept (as close to zero as possible). The dependent variable was ETo estimated by the PM-ETo and the independent variable was the ETo estimations by the Abtew equation. Therefore, the calibration processes tend to have a new regression relationship with a slope as unity and intercept as zero. The Abtew ETo equation was calibrated and validated using all eleven weather stations. The data from 1990 to 2003 were used for the equation calibration and data from 2004 to 2012 were used for the validation. This partitioning is due to the need of more data for training the equation as suggested by [28] [48] .
2.4. Evaluation Criteria
Simple linear regression was used for models comparison with reference to the FAO-PM. Root mean squared error (RMSE), relative error (RE), mean bias error (MBE), and the absolute mean error (AME) were also used for model evaluation and calculated as follow:
(11)
(12)
(13)
(14)
where, Ei is the estimated ETo with FAO-PM under limited data and the temperature ETo models; and Oi is ETo estimated with FAO-PM model with full dataset, at the ith data point and n is the total number of data points.
3. Results and Discussion
3.1. Performance of the FAO-PM Equation under Limited Data
The FAO-PM equation performed relatively well under limited data conditions in the semiarid and arid conditions in Mali. Under solar radiation missing, the FAO-PM (ETo-Rs) overall underestimated the daily ETo with the best performance at Nara in the Sahelian zone and Tessalit in the Sahara zone. The regression slope between FAO-PM ETo and ETo-Rs was very high and varied from 0.934 to 0.999 and the R2 ranged from 0.68 to 0.96 (Table 2). ETo underestimation
Table 2. Comparison between FAO-PM ETo computed from full Data set and FAO-PM ETo computed with limited data when Rs, RH, and u2 are missing.
was revealed through the MBE that average −0.18 mm/day and the MAE average of 0.38 mm/day (Table 3). Under missing RH conditions and when the actual vapor pressure is estimated with Tmin, ETo-Tmin basically underestimated daily ETo at lower rate of daily evapotranspiration less than 6 mm/day (Figure 1). The largest underestimation was observes at Nara and Tessalit when the MBE was −0.66 and −0.69 mm/day, and the MAE was 0.66 and 0.69 mm/day, respectively (Table 3). The least RE was observed at Bougouni (7.1%), the largest RE (12%) was observed at Nara and the average RE was 10% that showed the applicability of the method in the case of missing RH data. Large overestimation of daily ETo between 7% and 14% that averaged 10% (Figure 1) was obtained when the global average wind speed of 2 m/s was used (Table 2) with RMSE that varied from 0.51 and 0.85 mm/day, RE from 7.7% to 15.5%, and MBE from 0.44 and 0.79 mm/day (Table 3). However, very high R2 varying from 0.87 and 1.0 was obtained between FAO-PM ETo with full data set and ETo-um. When Rs and RH data are missing, ETo-RsTmin underestimated the daily ETo with RMSE ranging from 0.78 to 1.08 mm/day, MBE from −0.58 to −1.01 mm/day, high RE that varied from 12.5% to 19.6% (Table 3). ETo overestimation average of 7% was observed in the case of missing Rs and wind speed with no geographical specificity. In this case, RE was as high as 17.6% at Bougouni and average 13.9% (Table 3). Under missing wind speed and RH data, ETo-umTmin had relatively low RMSE that varied from 0.20 to 0.56 mm/day and low RE always less than 10% and averaging 6.7%, MBE average of −0.25 mm/day, and MAE average of 0.30 mm/day. The lowest RE was obtained at Bougouni (4%) when the null MBE was observed and the lowest MAE of 0.14 mm/day. When the RS, RH, and u2 are missing the FAO-PM has the poorest performance at all locations across Mali with large RMSE varying from 0.64 to 0.84 mm/day, high RE averaging 13%, and MBE varying from −0.75 to −0.29 mm/day, and MAE averaging 0.59 mm/day (Table 3). Therefore it is not recommend using the FAO-PM equation under this condition across the semiarid and arid conditions across Mali.
The results of this study are in agreement with previous research under similar climatic conditions. Similar results were reported by [4] during their study across Burkina Faso. Rojas and Thepadia [37] reported better results of the FAO-PM ETo equation when using wind speed data from a neighboring site in northeast Louisiana compared to the adoption of the global average wind speed of 2 m/s with a mean ratio of 0.98 and MAE of 0.56 mm/day. Trajkovic and Kolakovic [8] reported that the discrepancies between ETo under full data set and ETo under limited data set increased with increasing number of estimates weather parameters. Popova et al. [49] reported that the ETo estimation by FAO-PM under limited data provided accurate estimates of ETo with small standard errors of estimates. In contrast, [50] showed small differences in terms of MBE varying −0.22 to 0.25 mm/day and small RMSE varying 0.06 - 0.50 mm/day when they comparing ETo-RS and FAO-PM ETo with full data set in Korea. The results are close enough to the results of [51] who indicated that under missing Rs data, temperatures could be used to derive Rs for ETo estimation
Table 3. Statistical indices for the evaluation of FAO-PM ETo computed with limited data when Rs, RH, and u2 are missing.
Figure 1. Relationship between the FAO-PM daily ETo computed with full data and the daily ETo by the FAO-PM under limited data conditions.
in the Mediterranean environment. They reported very high R2 ≥ 0.98 using this method and RMSE varying from 0.42 to 0.71 mm/day. Wang et al. [52] indicated accurate estimates of ETo when RH is missing in Malawi. Jabloun and Sahli [34] also reported similar results under semiarid conditions in Tunisia. These results corroborated the finding of [4] for their study under semiarid climate in Burkina Faso. In Southern Ontario, Canada, ETo-Tmin overestimated ETo up to 12% as reported by [35] . Kwon and Choi [50] reported large RMSE of 0.6 - 0.73 mm/day under missing RH data in Korea. Under missing wind speed data, [34] reported close to unity regression slopes and very high R2 (>0.96) in Tunisia. [4] [8] [34] [37] reported that the use of global wing speed average of 2 m/s should be replaced by the local average wind speed data that provided better estimates of ETo. Under missing RH and u2, the results of this study are in agreement with [35] who reported that under missing RH and u2, FAO-PM method showed good estimation of daily ETo in the Southern Ontario, with RMSE < 0.53 mm/day. However, when three climatic variables were mission, the results showed the poorest performance of the FAO-PM model as reported by [34] [50] [53] .
3.2. Evaluation of the Selected ETo Equations
The tested ETo equations had different performance at the eleven weather stations in Mali. The Jensen-Haise equation systematically overestimated ETo at all sites (Table 3) with RMSE ranging from 1.6 to 2.0 mm/day, the RE within the range of 26.3% to 37.6% averaging 32.1%, and average AME of 1.55 mm/day (Table 5). The highest overestimation was recorded at Bougouni, Senou and Sikasso, all under the Sudanian climate. The Hansen, Christiansen, and Irmak ETo equations had similar performance and slightly underestimated ETo (Table 4). The regression slope varied from 0.79 to 0.94, from 0.76 to 0.92, and from 0.79 to 0.95 for the Hansen, Christiansen, and Imak equations, respectively, and the RMSE varied from 0.52 to 1.44 mm/day, from 0.63 to 1.66 mm/day, and from 0.45 to 1.43 mm/day for the respective ETo equations (Table 5). The relative error averaged 17%, 19%, and 15% under the Hansen, Christiansen, and Imak equations while the MBE averaged −0.75, −0.85, and −0.69 mm/day for the respective equations. These three ETo equations had the best performance at Bougouni under the Sudanian climate and the worst performance at Nara under the Sahelian climate.
The Abtew dimensionless coefficient K varied with locations and was 58.67, 57.75, and 57.55 at Bougouni, Senou, Sikasso, respectively, in the Soudanian zone; 57.54, 55.58, 55.08, 55.96, 56.94, 55.37 at Koutiala, Mopti, Nara, Niono, San, and Segou, respectively, for the Sahelian Zone; and 57.97 and 56.31 at Hombori and Tessalit, respectively, for the Sahara zone. Overall, there was not particular correlation between K and the geographical coordinates of the weather stations. However, good correlation between K values and the latitudes of the weather stations was found only for the Sahelian zone (K = 30.864 × Altitude − 1438.5 with R2 = 0.63) and this relation could be introduced into the original Abtew equation for regionalization under the Sahelian conditions. The Abtew equation showed the best performance among the ETo equations, slightly better than the [39] (Dr-Al) and [40] (Harg.) equations (Table 4 and Table 5). The regression slopes varied from 0.98 to 1.03, from 0.91 to 0.99, and from 0.90 to 1.07 for the Abtew, Dr-Al and Harg equations with average R2 of 0.86, 0.74 and 0.79, and
Table 4. Comparison between FAO-PM ETo computed from full data set and ETo computed with selected models.
MBE of −0.05, −0.25 and −0.05 mm/day for the respective ETo equations. The Abtew equation obtained the least relative errors that averaged 9.83% and AME of 0.41 mm/day. They showed better performance under the Sudanian and Sahelian semiarid climates. The Tabari 1 and Tabari 2 ETo equations did not show good performance under the Sudano-Sahelo semiarid and Saharan arid conditions in Mali. Both equations systematically underestimated daily ETo across the study area with RMSE that averaged 1.53 and 1.54 mm/day, MBE of −1.42 and −1.41 mm/day, and AME of 1.42 and 1.42 mm/day, respectively. The Tabari ETo equations showed the best performance with the lowest RMSE of 1.04 and 1.01 mm/day and the lowest AME of 0.93 and 0.86 mm/day at Bougouni in the Sudanian semiarid climate zone (Table 5). It can be deducted that the Abtew ETo equation was revealed the best among the nine temperature ETo equations which obtained the best evaluation indices across the study area. Therefore, it
Table 5. Statistical indices summary for the reference evapotranspiration equations evaluation.
could be used as the most adapted and specific ETo for the Sudanian, Sahelian, and Saharan climate in Mali. However, with 10% of RE, it should be adjusted to the Malian climate to improve its performance using proper calibration to the local conditions.
The results of this study showed better performance of the Hargreaves and Abtew under the semiarid and arid climate in Mali than under the Western Australian semiarid climate where [24] reported RMSE of 0.94 mm/day, MBE of −0.38 mm/day and RE of 17% for the Hargreaves equation while the Abtew equation induced 14% of relative error, RMSE of 1.014 mm/day and MBE of −0.4 mm/day. Contradictory, the Irmak’s model performed well with the lowest value of MBE of 0.27 mm/day at Pantnagarin India among twenty ETo models [54] . The performance of Jensen-Haise equation at Tessalit is in agreement with [55] who reported large discrepancies when using this equation under extremely arid condition like the Saharan arid climate at Tessalit. Similar to the results of this study, the Jensen-Haise equation presented also the highest RMSE of 1.63 mm/day with a relative error over 40% in the State of Rio de Janeiro, Southeast of Brazil [56] . However, [57] reported that the Jensen-Haise equation was the best among 23 ETo methods evaluated under the extremely arid climate conditions in the central Saudi Arabia. Kingston et al. [58] reported the uncertainty of the Jensen-Haise that provided the highest estimate of ETo at 20˚N, but the lowest ETo between 50˚S - 60˚S and [59] reported that Hansen equation was one of the best two performing ones with the least average monthly error in Greece. While the Hargreaves equation showed overall good performance under the semiarid and arid climates in Mali similar to the results of [60] in eastern arid and semiarid regions of Iran, it overestimated ETo under humid climate in northeast Louisiana’s [37] . Jensen-Haise model showed inaccurate estimation of ETo at California wit RMSD of 4.5 mm/day, and 2.36 mm/day at Bushland (Texas) and Davis (California), respectively, [19] . Under semiarid conditions in the Southern Spain, the Hargreaves equation achieved large values of the MBE ranging from 0.74 to 1.13 mm/day and RMSE from 0.46 to 1.65 mm/day representing large maxima under- and overestimation of 24.5% and 22.5%, respectively [61] . Sabziparvar and Mirgaloybayat [61] reported good performance of Irmak equation with comparison to FAO-PM model at coastal sites in Iran with low altitude and high relative humidity.
3.3. Simple Multi-Model Ensemble (MME)
The simple multi-model ensembles (MME) achieved better performance than every single equation of the group of nine ETo equations (Figure 2 and Figure 3). The regression slope between FAO-PM ETo and MMEs ETO estimates varied from 0.858 to 0.979 with very high R2 greater than 0.79 (Figure 3, Table 4). While the RMSE varied from 0.37 to 2.29 mm/day for all temperature methods, it ranged between 0.31 and 1.02 mm/day for the MMEs (Table 5). The MMES also achieved the lowest RE within the range of 6.3% - 16.2% while RE varied from 8.49 to 36.61 % for the group of the equations. The Abtew, Hargreaves and Droogers-Allen equations had performed better than the MMEs. The Jensen-Haise, Chriatiansen, Tabari 1 and Tabari 2 equations were revealed the worst compared the MMES (Table 5). Overall, the MMEs method improved the accuracy of the daily ETo estimation across Mali and could be used when a selection of some ETo equation should be used for Eto estimation under inaccurate and or limited data conditions. The MME method was successful applied by
Figure 2. Relationship between the daily FAO-PM ETo estimates and the daily ETo computed by the ETo models under evaluation.
Figure 3. Relationship between the daily FAO-PM ETo estimates and the daily MME ETo estimates at all weather stations (1990-2012).
[62] and [63] who reported more than 10% reduction in model uncertainty under was achieved under MME when estimating irrigation water requirement through multi-model ensemble. Wang et al. [64] also reported the outperformance of the multiple models in a nowcast system for the monitoring of the current state of soil water.
3.4. Calibration and Validation of the Abtew Equation
The Abtew ETo equation was revealed the best among the selected equations with a regression slope between the FAO-PM ETo estimates and the original Abtew ETo of 0.9898 and high R2 of 0.86 (Figure 4(a)) All weather stations combined, the simple linear regression between the daily FAO-PM ETo estimates and the calibrated Abtew ETo estimates for the 1990-2003 period is presented in Figure 4(b). The regression showed the good fitness of the calibration with regression slope of 1.0035 close to unity and R2 = 0.89 (Figure 4(b)). The validation of the calibrated Abtew equation is presented in Figure 4(c). With a regression slope of almost unity (1.0288) and the R2 equal to 0.92, the calibrated equation showed good performance and can be used for ETo estimation under the Sudanian, Sahelian, and Saharan climates in Mali. Further, the calibration process improved the RMSE of the ETo estimates from an average of 0.55 mm/day to 0.31 mm/day representing 41% improvement. The highest improvement of 63% by calibration was achieved at Tessalit while the lowest improvement of 30% was achieved at Sikasso (Table 6). The calibrated Abtew
Figure 4. Relationship between (a) the original Abtew ETo estimates and the FAO-PM ETo estimates for the 1990-2012 period; (b) Calibrated Abtew ETo estimates and the FAO-PM ETo estimates for the 1990-2003.
Table 6. Improvement of the RMSE of the calibrated Abtew ETo equation.
equation performed as well during the calibration period as the validation period with the RMSE of the ETo estimated was 0.32 mm/day. The results of this study corroborated with the findings of [65] who reported that the Abtew model showed the best overall performance with respect to the data from all available climate stations of Central Greece, and [29] in Gansu Province, northwest China. Djaman et al. [5] reported the good fitness of the Abtew equation with FAO-PM equation under the semiarid climate in Tanzania and Kenya. Xu et al. [48] also have concluded that the simple Abtew equation can be used in the state of Vaud in Switzerland when other meteorological data except solar radiation are not available. The Calibrated Abtew equation to be used under the semiarid and arid climates in Mali climate is recommended for the study area.
4. Summary and Conclusion
The FAO-PM ETo equations using missing climatic data and nine temperature reference evapotranspiration methods were evaluated for their accuracy relative to FAO-PM equation under the Sudano-Sahelo-Saharan climate across Mali for the period of 1990-2012. The results showed that under limited data conditions, the FAO-PM equation showed good performance when solar radiation (Rs), relative humidity (RH), and wind speed (u2) are lacking solely and when both relative humidity and wind speed data are missing with RMSE varying lower than 0.58 mm/day and average relative error of 6.7%. The Abtew ETo equation that requires solar radiation and maximum temperature showed the best performance across all three climatic zones in Mali. The Jensen-Haise equation systematically overestimated ETo while the Tabari’s two equations underestimated ETo. The Irmak, Hargreaves, Hansen, Christiansen performed relatively well in the study area. With 41% improvement of the performance of the best performing Abtew equation, a new form of the Abtew equation is recommended for ETo estimation across Mali and similar climatic conditions. Also in the case of non-available u2 and RH data, the simplified forms of the FAO-PM equation is recommended for use for reasonable ETo estimation across Mali.