Projections of evapotranspiration form the basis of future runoff and water availability assessment in a climate change context. The scarcity of data or insufficiency of time/funds compels the application of simple reference evapotranspiration (ETo) methods requiring less meteorological inputs for ETo projections which adds uncertainty to the projected changes. This study investigates the bias in ETo climate change signals derived from seven simple temperature- and radiation-based methods (Blaney–Criddle, Hargreaves–Samani, Schendel, Makkink, Turc, Jensen–Haise, Tabari) compared with that from the standard Penman–Monteith FAO 56 method on the basis of 12 general circulation model (GCM) outputs from the Coupled Model Intercomparison Project Phase 5 for central Belgium for four future greenhouse gas scenarios (RCP2.6, RCP4.5, RCP6.0, RCP8.5). The results show the lack of conformity on the amount of ETo changes between the simple and standard methods, with biases of over 100% for some simple methods. The uncertainty affiliated with ETo methods for monthly ETo changes is smaller but of comparable magnitude to GCM uncertainty, which is usually the major source of uncertainty, and larger for daily extreme ETo changes. This emphasizes the imperative of addressing the uncertainty associated with ETo methods for quantifying the hydrological response to climate change.
Evapotranspiration (ET) is a main flux term in the water cycle and an important process in the climate system, thereby playing a key role in water and energy balance. It is also vital for climate change impact studies, drought prediction and monitoring, agricultural water requirement planning and effective management of water resources. Specifically for hydrology, reference ET (ETo) is a necessary component for computing actual ET for a watershed as a function of soil moisture deficit and a key input parameter for hydrological water-balance models (Bergström et al. 2001; Obeysekera 2013; Tabari & Hosseinzadeh Talaee 2014; Meng et al. 2016). ETo is, moreover, used to assess the hydrological impact of climate change at catchment scales (Taye et al. 2011; Ntegeka et al. 2014; Vansteenkiste et al. 2014; Rudd & Kay 2016).
Climate change is one of the most critical global problems of our time (IPCC 2007), which may influence watershed hydrology and water availability by changing the patterns of precipitation and ET (King et al. 2015). Future changes in ET rate together with changes in precipitation patterns will affect runoff and water balance and consequently water infrastructure design (Prudhomme & Williamson 2013; Mohor et al. 2015; Koedyk & Kingston 2016).
There are numerous methods to estimate ETo based on meteorological data ranging from simple temperature (e.g., Hargreaves–Samani: Hargreaves & Samani 1985) or radiation (e.g., Jensen–Haise: Jensen & Haise 1963) based methods to complex Penman family methods (e.g., Penman: Penman 1948; Kimberly–Penman: Wright 1982). Among them, the Penman–Monteith FAO 56 (PMF-56) method was strongly recommended by the United Nations Food and Agriculture Organizations (FAO) and the American Society of Civil Engineers (ASCE) as a standard method for ETo estimation (Allen et al. 1998; ASCE-EWRI 2005), because the PMF-56 method incorporates most of the meteorological variables controlling the ET process. Nevertheless, to simplify the process of ETo estimation due to limited time and/or budget and in some cases because of unavailability of the full meteorological data required for the standard PMF-56 method, simple, empirical equations have been used in many climate change studies for this purpose (e.g., Immerzeel et al. 2012; Bastola 2013; Lu et al. 2013; Capell et al. 2013; Kopytkovskiy et al. 2015). The simple equations relate ETo to an easily measured process variable such as temperature and do not consider the effect of ‘advective’ variables (i.e., wind speed and vapor pressure) (Donohue et al. 2010). However, ETo is primarily a function of four key weather variables of radiation, wind speed, humidity (or vapor pressure), air temperature and these variables must be considered for a better understanding of ETo changes under climate change conditions (McVicar et al. 2012). The ETo changes projected by simple approaches may differ substantially from those by the standard method that may add uncertainty to derived climate change signals (CCS) and future water balance projections (Kingston et al. 2009; Bae et al. 2011; McAfee 2013; Thompson et al. 2014; Wang et al. 2015). Therefore, the existence of potential systematic bias or uncertainty in the derivation of ETo CCS by using simple ETo methods needs to be considered. The uncertainty in ETo CCS by simple methods is, however, largely unknown. This research gap is addressed in this paper. More specifically, the bias in ETo CCS obtained from seven simple methods (Blaney–Criddle, Hargreaves–Samani, Schendel, Makkink, Turc, Jensen–Haise, Tabari) compared with those from the PMF-56 method as a benchmark is investigated under different time resolutions, i.e., seasonal, monthly and daily.
The general circulation model-Coupled Model Intercomparison Project Phase 5 (GCM-CMIP5) ensemble considered for this study consists of 12 GCMs (listed in Table S1 with their respective resolution; Table S1 is available with the online version of this paper). All the meteorological data needed for ETo estimation by the PMF-56 method (which include those needed for seven simple methods) were obtained from 264 CMIP5 GCM runs (44 runs for each of six meteorological variables including maximum, minimum, and mean temperature, radiation, wind speed, and humidity). In fact, all the GCMs providing the complete set of meteorological inputs for ETo estimation were selected. Although some authors recommend 50-year periods (Bonell & Bruijnzeel 2005), in this study the period 1961–1990 was applied as the control period. This period is recommended by WMO (2007) and commonly applied in climate change studies. The period 2071–2100 was considered as the scenario period. Four 21st century scenarios for future greenhouse gas concentrations including RCP2.6, RCP4.5, RCP6.0, and RCP8.5 were considered, where RCP8.5 represents the highest concentrations and RCP2.6 corresponds to the lowest concentrations (Moss et al. 2010). The analyses were done for the GCM grid cell covering central Belgium.
After extracting all the meteorological data from the GCMs, ETo is calculated by different methods, that are briefly explained in the next section. Daily ETo values are averaged to obtain monthly ETo. Afterwards, CCS are derived for both daily and monthly ETo as the difference between the values for the scenario period (2071–2100) and those for the control period (1961–1990). The changes in daily ETo are shown by means of box-and-whisker plots for summer and winter seasons. The calculation of these results involves the following steps:
For each ETo model, separate the ETo results for summer (JJA: June, July, August) and winter (DJF: December, January, February) seasons.
Combine CCS of all GCMs.
Compute 5%, 25%, 50%, 75%, and 95% percentiles of the signals (%) and show the results in a box-and-whisker plot.
High, mean, and low climate scenarios for monthly or daily ETo change are later defined as 95%, 50%, and 5% quantiles of the relative changes. For the monthly ETo, the changes are considered as such; for the daily ETo, the changes are considered as changes in the quantiles (percentiles). As the changes in high ET rates are of vital importance for drought monitoring and agricultural water planning and management, the changes in high daily ETo quantiles (percentiles) are also investigated.
In this study, the ETo CCS derived from the PMF-56 method are compared with the ones from three temperature-based (i.e., Blaney–Criddle, Hargreaves–Samani and Schendel) and four radiation-based (i.e., Makkink, Turc, Jensen–Haise, and Tabari) methods.
Blaney–Criddle (Blaney & Criddle 1950)
Hargreaves–Samani (Hargreaves & Samani 1985)
Schendel (Schendel 1967)
Jensen–Haise (Jensen & Haise 1963)
Makkink (Makkink 1957)
Turc (Turc 1961)
Tabari (Tabari et al. 2013)
Bias in the mean scenario of seasonal ETo calculated based on the simple ETo methods is presented in Table 1. As the results indicate, the performance of the simple methods changes seasonally. The Hargreaves–Samani, Makkink, and Tabari methods underestimate the mean scenario and the rest of the methods overestimate it. In winter, the best performance in estimation of the mean scenario for seasonal ETo is produced by the Hargreaves–Samani, Tabari and Jensen–Haise methods, while the Makkink and Tabari methods perform the best in spring. The mean scenario for summer ETo is underestimated by most of the simple methods, among which Turc and Schendel show the smallest biases. For autumn, the Tabari and Makkink methods more closely reproduce the PMF-56-driven scenario. Overall, among all the simple methods, Schendel has the worst performance in estimation of monthly and seasonal ETo changes.
The results obtained indicate that the ETo change signals projected by the simple methods broadly follow the pattern of those by the PMF-56 method. However, the magnitude of the changes differs substantially between the simple and standard method. The main reason for such difference is the inclusion of more meteorological variables controlling the ET process in the PMF-56 method. Moreover, simple ETo approaches with the same meteorological inputs present different CCS. This difference stems from the difference in empirical formulation for each approach and the difference in climatic regions for which they were originally developed. Let us consider the drawback of the Hargreaves–Samani method as the most widely used simple approach in hydrological impact studies of climate change. First of all, the Hargreaves–Samani method uses extraterrestrial radiation rather than solar radiation, implying the use of maximum possible radiation and neglecting atmospheric transmissivity. The existence of high moisture content in the atmosphere in humid regions influences transmissivity or sky clearness and leads to an increase of the attenuation of solar radiation at the surface (Fontenot 2004). The other important variable in humid regions is atmospheric moisture, which is not considered in the Hargreaves–Samani method. Atmospheric moisture (directly correlated with vapor pressure) has an indirect relationship with ET, and in fact, ET decreases as atmospheric moisture increases (McKenney & Rosenberg 1993; Tabari 2010; Moratiel et al. 2010). These two important factors result in bias in the estimates of the Hargreaves–Samani method. The calibration of the method to local conditions may decrease this bias (Tabari & Hosseinzadeh Talaee 2011).
The findings show that choosing an alternative method to the PMF-56 method with less data demand is difficult. The performance of the simple methods varies remarkably by time scale and season. Taking the mean scenario into account, the Hargreaves–Samani, Makkink, Turc and Tabari methods are the best choices for seasonal and monthly ETo climate change studies in winter, spring, summer, and autumn, respectively, while the Turc, Jensen–Haise and Hargreaves–Samani methods perform best for investigation of the changes in high ETo percentiles. This suggests using different simple methods for different time scales and seasons in the absence of full meteorological data for the PMF-56 method application.
In order to investigate the importance of the uncertainty associated with ETo methods, it is compared with GCM uncertainty which is considered as the main source of uncertainty in climate change studies (Wilby & Harris 2006; Kay et al. 2009; Chen et al. 2011; Exbrayat et al. 2014) and also that of future greenhouse gas scenarios (RCP). For the monthly ETo change and also taking the full range of daily ETo change into consideration, the uncertainty in GCMs is larger than that in the ETo methods and RCPs (Table 2). For both mentioned cases, the ETo uncertainty is of comparable magnitude to the GCM uncertainty. As for the changes in high daily ETo percentiles, the uncertainty related to the ETo methods is the dominant source of uncertainty for both summer and winter seasons, while RCP has the lowest contribution to the total uncertainty in ETo changes. Generally, the uncertainty analysis shows that the uncertainty in ETo methods is noticeable for all the cases considered. This shows the importance of addressing the uncertainty associated with ETo methods next to that in GCMs and future greenhouse gas scenarios as is traditionally done.
|Uncertainty||Monthly ETo||Daily ETo percentile changes (full range of||Daily ETo percentile changes (>90%)|
|Uncertainty||Monthly ETo||Daily ETo percentile changes (full range of||Daily ETo percentile changes (>90%)|
CONCLUSIONS AND RECOMMENDATIONS
In this study, the difference between the ETo CCS derived from the standard PMF-56 method and those from seven simple methods (Blaney–Criddle, Hargreaves–Samani, Schendel, Makkink, Turc, Jensen–Haise, Tabari) was investigated. The meteorological data required for ETo estimation by the methods were obtained from 12 CMIP5 GCMs for control (1961–1990) and scenario (2071–2100) periods. The results show that the standard and the simple methods generally agree on ETo increase for all time scales, but disagree on the magnitude of the increase. Biases of over 100% are found for CCS extracted from some simple methods such as Turc and Schendel.
A comparison between the uncertainties related to ETo methods and GCMs shows that the uncertainty associated with ETo methods is of comparable magnitude to GCM uncertainty as the main source of uncertainty, suggesting that the uncertainty in the choice of ETo methods must not be neglected in such studies. Use of further advanced methods for ETo estimation, which is a function of time scale and season (as the results of this study show), may reduce the uncertainty associated with ETo methods. Comparing the ETo uncertainty with the other sources of uncertainty, such as GCM boundary conditions, was impossible in this study because of the limited number of GCM runs providing all meteorological inputs for the ETo methods. This type of comparison can be made for precipitation as another important input for hydrological modeling in future studies. The uncertainty related to the choice of a downscaling method was not addressed in this study. Thus, further work is needed to assess the degree to which the choice of the statistical method to transfer the climate change signal to the hydrological model inputs, i.e., the downscaling method, contributes to the total uncertainty. Moreover, the effect of the different sources of uncertainties affiliated with precipitation and ETo may be explored for hydrological impact projections.
To conclude, considering the large bias of the simple ETo methods, a fully physically based formulation of ETo like the PMF-56 method is recommended for projecting ET CCS. This, moreover, requires that the climate modeling centers produce and provide model outputs for all meteorological variables governing the evaporative process. In the case of data unavailability, the simple methods must be used with caution. The bias in the results of the simple methods depends on the time scale and season, as is shown in this study for central Belgium. This bias obviously may differ for other regions of the world.
The authors acknowledge the World Climate Research Programme's Working Group on Coupled Modelling, and the climate modeling centers involved for the CMIP5 model outputs.