Next Article in Journal / Special Issue
Complex Extreme Sea Levels Prediction Analysis: Karachi Coast Case Study
Previous Article in Journal
Pricing Constraint and the Complexity of IPO Timing in the Stock Market: A Dynamic Game Analysis
Previous Article in Special Issue
Simulation of S-Entropy Production during the Transport of Non-Electrolyte Solutions in the Double-Membrane System
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Reference Evapotranspiration Modeling Using New Heuristic Methods

1
College of Hydrology and Water Resources, State Key Laboratory of Hydrology-Water Resources and Hydraulic Engineering, Hohai University, Nanjing 210098, China
2
School of Information Science and Engineering, Wuhan University of Science and Technology, Wuhan 430081, China
3
School of Hydropower and Information Engineering, Huazhong University of Science & Technology, Wuhan 430074, China
4
Hubei Provincial Key Laboratory for Operation and Control of Cascaded Hydropower Station, China Three Gorges University, Yichang 443002, China
5
Department of Civil Engineering, School of Technology, Ilia State University, 0162 Tbilisi, Georgia
6
Department of Civil Engineering, Faculty of Engineering, University Malaya, Kuala Lumpur 50603, Malaysia
7
CERIS—Civil Engineering Research and Innovation for Sustainability, Instituto Superior Técnico—Universidade de Lisboa, 1049-001 Lisboa, Portugal
8
Department of Irrigation and Drainage, University of Agriculture, Faisalabad 38000, Pakistan
*
Authors to whom correspondence should be addressed.
Submission received: 7 April 2020 / Revised: 2 May 2020 / Accepted: 5 May 2020 / Published: 13 May 2020

Abstract

:
The study investigates the potential of two new machine learning methods, least-square support vector regression with a gravitational search algorithm (LSSVR-GSA) and the dynamic evolving neural-fuzzy inference system (DENFIS), for modeling reference evapotranspiration (ETo) using limited data. The results of the new methods are compared with the M5 model tree (M5RT) approach. Previous values of temperature data and extraterrestrial radiation information obtained from three stations, in China, are used as inputs to the models. The estimation exactness of the models is measured by three statistics: root mean square error, mean absolute error, and determination coefficient. According to the results, the temperature or extraterrestrial radiation-based LSSVR-GSA models perform superiorly to the DENFIS and M5RT models in terms of estimating monthly ETo. However, in some cases, a slight difference was found between the LSSVR-GSA and DENFIS methods. The results indicate that better prediction accuracy may be obtained using only extraterrestrial radiation information for all three methods. The prediction accuracy of the models is not generally improved by including periodicity information in the inputs. Using optimum air temperature and extraterrestrial radiation inputs together generally does not increase the accuracy of the applied methods in the estimation of monthly ETo.

1. Introduction

For irrigation management, agricultural processes, and hydrological cycles, the reference evapotranspiration (ETo) is considered as one of the fundamental variables that should be accurately estimated [1]. For example, it is essential to adequately estimate the ETo in order to figure out the amount, timing, schedule, and frequency for the agricultural procedure and irrigation management [2]. Additionally, nowadays under the expected increase of water scarcity and the reduction of food production, it is vital to get proper knowledge of the estimation of the ETo at different time increments in order to optimize the management of the available water resources, whether conventional or non-conventional water resources [3,4]. Estimation of ETo is fundamental due to the ETo key role in affecting the hydrological cycle [5]. Accurate estimation of ETo will help in understanding the impacts of human influences on the hydrological cycle and in improving water resource management [6].
Estimation of ETo is a challenging task due to complex interactions between meteorological and site-specific factors [7]. Generally, there is a need to standardize a single robust method to be used at a large scale (i.e., region, large basin) for reference evapotranspiration-ETo estimation, based on limited climatological and hydrological data. However, we are aware that such a method would still require further adjustment and calibration, considering the local conditions of the study site. In this context, even well-recognized methods in use so far, such as Penmen–Monteith (PM), require adjustment on many occasions [6,8]. Over the last 50 years, several scientists developed different methods to calculate the ETo; however, these methods have been subjected to local conditions and have several limitations to their applicability on a large scale. One of the most popular methods that have been validated and recommended by the Food and Agriculture Organization (FAO) is the PM method, and it is now recognized as FAO 56 PM. Although the FAO 56 PM method has been proved to be the most accurate one and achieved the most minimal errors, it required much information to be applied, which is mostly unavailable in several areas worldwide. Therefore, considering physical parameters and physically-based formulas might not be feasible and/or realistic for estimating the ETo.
Due to the complexity of estimating the ETo-based physical process, several machine learning methods were employed to estimate the ETo during the last two decades of the 21st century [8]. More recently, during the last two decades, artificial neural network (ANN) models have been examined for estimating the ETo, such as the multi-layer perceptron neural network (MLP-NN), fuzzy logic (FL), the adaptive neuro-fuzzy neural network (ANFIS), and least square support vector regression (LSSVR) [8,9,10,11]. The motivation for utilizing the ANN models is that these methods can provide high accuracy and robustness, are modeless, and can easily handle big data [12,13].
Although the ANN showed outstanding performance for estimating the ETo, there is a need to include several meteorological and hydrological variables in the model input in order to achieve relatively high estimating accuracy for the ETo utilizing ANN, [14,15]. Additionally, it has been reported in a few studies that the ETo prediction model using the ANN experienced a few challenges, such as over-fitting and the selection of the appropriate input variables based on the climate conditions of the study area [16].
Recently, the M5 model tree (M5RT) has been developed and considered as one of the most common methods used for simulation [17]. The process of the M5RT is designed to divide the data into sub-regions and consider the tree applying data of each sub-region. The input space should be divided into sub-regions, and then a fitted linear model could be used for each sub-region, [18]. There are a few studies which showed that the M5RT could achieve proper estimation accuracy but requires extensive sampling input data [19]. In general, M5RT is considered a straightforward procedure that could be worthwhile for estimation applications, especially with the availability of quite large amounts of data [17,20,21].
Furthermore, the ANFIS method has successfully shown a high ability to provide accurate predictions in different engineering applications, especially for the hydrological and climatological variables [22,23,24,25]. Even with a limited amount of available data of climate parameters, ANFIS could achieve accurate estimation for daily pan evaporation [24]. Further enhancement of the ANFIS model, including either grid partitioning (GP) or subtractive clustering (SC), has been developed by Sanikhani et al. [26]. It showed outstanding performance for estimating the daily evaporation; however, it required several input variables and a relatively large amount of data. More recent novel enhancement for the ANFIS model has been developed, which is defined as a dynamic evolving neural-fuzzy inference system (DENFIS) in order to be more suitable for dynamic time series prediction applications [27,28]. Conceptually, the DENFIS is designed to generate and update new fuzzy rules during the learning of the system. Such a procedure allows the DENFIS to calculate the desired output according to them-most activated rules, which have been chosen dynamically from the set of the fuzzy rules [27,28,29]. It could be noticed from the DENFIS structure and procedure and the nature of ETo that the DENFIS model has the potential to adequately and effectively estimate the ETo. To the knowledge of the authors, there is no study in the literature related to the application of DENFIS in modeling ETo.
More recently, LSSVR has been known as a useful tool for several engineering prediction applications. LSSVR has a simple and effective structure; however, it includes three unknown parameters that should be estimated and initialized at the beginning, which is considered one of the major disadvantages of the LSSVR. Kisi [19] showed that the more accurate the estimation of the LSSVR’s unknown parameters, the more accurate the estimation of the desired variable. However, the most common method to estimate these three unknown parameters is based on trial and error, as in many previous studies. One innovation for the development of a successful LSSVR model is to be integrated with a suitable optimization algorithm that could be able to search for the optimal estimation values for these three unknown parameters [19]. It has been reported that the gravitational search algorithm (GSA) is a useful optimization algorithm which is comparatively fast in convergence and free from trapping in local minima over the other optimization algorithms [30]. Therefore, integrating the LSSVR with (GSA) as an efficient optimizer for searching for the optimal values of the LSSVR’s unknown parameters could be a suitable solution to enhance the LSSVR’s ability to estimate the desired ETo. Available literature indicates that new methods are essential to enhancing ETo estimation preciseness and decreasing the model’s uncertainty. For this reason, novel hybrid heuristic soft computing methods were examined in this study for performing effective evapotranspiration modeling. A number of studies [31,32,33,34,35] have compared one or more of the hybrid heuristic models using several meteorological variables at different time steps, and the results obtained have shown that hybrid heuristic soft computing methods generally provided good accuracy. However, even though the aforementioned studies demonstrated successful applications of hybrid heuristic soft computing models, the most important point to note is that, generally speaking, these models do not use fewer independent input variables; instead, they require several inputs—i.e., U, RH, SH, Tmean, Tmax, and Tmin—and among them the U, RH, SH, and SR are the meteorological variables reported as the most significant factors influencing ETo and of primary importance. To the author’s knowledge, there is not any published study that investigates the potential of the LSSVR-GSA method in modeling ETo. First, we developed the hybrid of GSA and LSSVR for modeling ETo. Second, we used the already applied standalone soft computing method (M5RT) with only average temperature input, whereas most of the other studies in the literature used a mix of the Tmax and Tmin and several other meteorological variables during applications of such models. Third, we used extraterrestrial radiation which is easily obtained from the Julian date and latitude information as input to the models. Fourth, we demonstrated the effect of periodicity (month number) as an input variable for ETo estimation using limited climatic inputs.
In the current study, an attempt to develop a model to estimate the ETo based on a minimal number of input variables has been proposed utilizing advanced machine learning methods. LSSVR-GSA, DENFIS, and M5RT models have been developed to estimate the ETo considering only the temperature as an input variable. Data from three stations in the Jinsha river basin of China have been used in order to evaluate and examine all the proposed models. Additionally, in-depth, comprehensive analyses for the performance of all models have been accomplished in order to elucidate the logic and the cause behind each model’s estimation accuracy.

2. Materials and Methods

2.1. Case Study

The study uses monthly average temperatures (T), extraterrestrial radiation (Ra), and evapotranspiration (ETo) data of three meteorological stations, represented as 56004 (latitude 34.21 N, longitude 92.43 E), 56021 (latitude 34.13 N, longitude 95.78 E), and 56029 (latitude 33.01 N, longitude 97.01 E,) stations situated in the Jinsha river basin of China. The location of each station can be seen in Figure 1. These stations are operated by the China Meteorological Administration (CMA). Jinsha river basin consists of three provinces (Qinghai, Yunan, Sichuan, and Tibet Autonomous) covering the 473.2 × 103 km2 drainage area of China. The Jinsha river basin has a very vital role in regional and national economic development due to its contribution to irrigation, water supply, flood control, wood drift, tourism, and plentiful hydropower resources (58,060 MW). The mean annual rainfall in the Jinsha river basin is 750 mm, 90% of which occurs from May to October. The basin selected area belongs to humid warm temperate, which is the primary source of water, having annual mean rainfall of about 1200 mm with mean annual potential evapotranspiration of about 1350 mm.
In the study, the data periods of three meteorological stations from 1961 to 2012 were used. Available data were divided into three parts as training (50% of the whole data covering 1961–1986), validation (25% of the whole data covering 1986–1999), and testing (25% remaining part covering 2000–2012). The brief statistical characteristics of the user data are summed up in Table S1. Evapotranspiration data of three meteorological stations are calculated using the FAO-56 PM method. The detailed procedure of the FAO-56 PM method is described by Allen et al. [36] and She et al. [37].

2.2. Methods

2.2.1. Dynamic Evolving Neural-Fuzzy Inference System (DENFIS)

DENFIS is a neural-fuzzy modeling approach that evolves through gradual hybrid learning procedures resulting from the tuning of local elements regarding the new data entry to the model-specific categorization [27]. DENFIS is an extension of the evolving fuzzy neural network (EFuNN), which uses a clustering approach called the evolving clustering method (ECM) [38]. Dynamic features of the EFuNN make DENFIS useful for online adaptive systems also. ECM itself is an online fast-evolving clustering method that considers the maximum distance between points and the cluster’s center.
ECM is based on triangular membership functions, which serve as fuzzy sets and the utilization of an alternative weighting scheme for local learning of resulting parameters. Usually, the input and output neurons can be fuzzified by a fuzzy quantization approach [39]. The offspring fuzzy rules are consequently generated and updated while the system is in operation. The outputs resulting from DENFIS are computed through FIS-based activated fuzzy rules, which are dynamically chosen from a given fuzzy set of rules. Notably, the ECM used in the DENFIS interface utilizes a scatter partitioning of the input space to create fuzzy inference rules [38]. Uses of ECM improve the computational performance for the dynamic estimation of the cluster’s quantity in the given dataset by facilitating the finding of their current centers in the input space. A typical structure and detailed description of the EFuNN interface based on ECM from which the DENFIS approach is derived can be found at [27]. DENFIS uses a typical model architecture, which is based on the Takagi–Sugeno type fuzzy inference engine [34]. In this regard, DENFIS utilizes a model-based approach called the “lazy” learning approach. According to this approach, the fuzzy network estimates the position of each input vector in the attribute space and forms accordingly; meanwhile, a fuzzy inference system which predicts the output through a dynamic process is created during the incremental learning.
The classical “lazy” learning process of the fuzzy network utilizes a sample-based approach, wherein a small local model on demand is constructed regarding local samples taken from the closest query point. The learning process and governing equations used in the DENFIS model and applied in this study are described in detail in [39]. DENFIS models are being applied in different water resource-related areas and other disciplines as well. The use of DENFIS in rainfall-runoff modeling showed that results attained from the local learning model were significantly better than results achieved from physically-based models [38]. In another study related to evaporation modeling based on limited meteorological data, it was found that DENFIS models increased the accuracy of the outputs significantly [39]. Similar results were obtained in other studies as well [40]. DENFIS models were recently applied to predict solar radiation [41]; the results showed that DENFIS provides faster and much more accurate results compared to the other neuro-fuzzy models.

2.2.2. Least Squares Support Vector Regression (LSSVR)

Least squares support vector regression was initially proposed by [42]; it is a version of the support vector regression (SVR) algorithm improved by adopting the loss function differently from SVR and minimizing the square error. In the LSSVR interface model, the inequality constraints have been replaced with equality constraints by transforming quadratic programming issues into a linear equation to overcome the calculation issue of the large-scale dataset [43]. SVR itself is well known, and a robust algorithm applied in different areas mainly for data regression and also to overcome the over-fitting issue. LSSVR differs from SVR, mainly due to the types of functions used by each of them. LSSVR handles square errors instead of non-negative errors, and it uses equality constraints instead of inequality constraints, as opposed to the conventional SVR [44].
Similarly, to the DENFIS, LSSVR performs the training procedure in terms of solving a linear system instead of using a quadratic programming problem, which thus significantly improves the computation time and the accuracy of the model learning ability [45]. Figure 2 illustrates the structure of a typical LSSVR model. A detailed explanation about LSSVR applications, parameters, and respective governing equations is given by Yuan et.al [44]. LSSVR is widely applied in solving challenging problems in different areas. Wu and Peng [43] used LSSVR to improve the prediction performance of the potential wind power; they found that LSSVR can enhance the accuracy of up to 20% compared to other single or hybrid models. A similar conclusion was obtained by Lu et al. [46] who applied LSSVR for short-term forecasting of wind power. The application of the LSSVR has shown promising results in the prediction of nanoparticles as well [47].

2.2.3. Gravitational Search Algorithm (GSA)

GSA, which is based on the law of gravity and motion proposed initially by Rashedi et al. [30], represents one of the most effective optimization algorithms compared to other evolutionary algorithms. In GSA, each nod is characterized by four parameters: position, inertial mass, gravitational mass, and velocity. The location of each node corresponds to a solution of the problem, while the gravitational and inertia masses for each node were obtained utilizing a fitness function [48]. The location of the nodes can be expressed as follows:
X i = ( x i 1 , , x i k , , x i s )   i = 1 , 2 , , np
where x i k represents the location of the ith node for the kth dimension. The mass of each node is computed after calculating the fitness of the given population as follows:
X i = f i t i t w o r s t   t b e s t t w o r s t   t
M i t = m i t j = 1 N m j t
where fiti(t) and Mi(t) represent the fitness value and mass of the ith node at time t, respectively, whereas best(t) and worst(t) are the minimum fitness value and maximum fitness value, respectively, the gravitational acceleration of the node i, is calculated as follows: firstly, the force exerted by a large node on the node i is computed.
F i j k t = G c t M i t . M j t || x i t , x j t || 2 + ε .   x i k t x j k t
where Mi(t) and Mj(t) are the passive and active gravitational mass, respectively, corresponding to nodes i and j at the t generation; Gc(t) and ε are gravitational and small constants; xk i (t) and xk j (t) indicate the positions of the kth dimensions of nodes i and j at the t generation; and ‖xi(t), xj(t)‖ is the Euclidean distance between nodes i and j. The total gravitational acceleration of the ith node was calculated using the law of motion as follows:
a i k t = j = 1 ,   j i N r a n d .   F i j k t M i t
where a i k t   represents the gravitational acceleration of the node i in the kth dimension, and the rand represents a random variable with uniform distribution within the interval (0,1). The total gravitational force exerted on the node i in the kth dimension is calculated as follow:
F i j k t = M i t x a i k t = j = 1 ,   j i N r a n d .   F i j k t
Afterward, the speed and location of each node were updated as follows:
v i k t + 1 = r a n d   v i k t + a i k t
x i k t + 1 = x i k t + v i k t + 1
As highlighted above also, GSA utilizes the gravitational force as the direct form with which to communicate the cooperation of the nodes. The heavy nodes in GSA are processed, infer reasonable solutions, and move more gradually than lighter ones. Thus, the GSA searches for the ideal solution by appropriately calibrating the inertia and gravitational masses of nodes where every node provides a specific solution.

2.2.4. HLGSA (Hybrid LSSVR-GSA)

The combination of LSSVR with the gravitational search algorithm (GSA) has been reported to improve the forecasting accuracy and computational performance [43] significantly. To our knowledge, there are only a few studies applied to hybrid models, such as combining GSA with other ANN models. Yuan et al. [44] applied LSSVM–GSA to predict day-ahead wind power; they concluded that GSA provided higher accuracy and improved the computational duration. Similar conclusions were reported by Lu at al. [46] and Sun at al. [47] in the case of PM2.5 concentration prediction. The process of the evapotranspiration prediction model HLGSA using the hybrid of LSSVR and GSA methods consists of the following steps:
  • Firstly, divide meteorological datasets into training, validation, and testing parts.
  • Second, select the RBF kernel function and initial parameters for the HLGSA method to build the initial LSSVR model. The initial values of the parameters are set as follows: the range of penalty factor γ is 0.1 to 1500, the range of RBF parameter σ2 is 0.001 to 10, the iteration range is 20–30, the number of particles can be set up to 30, and constant alpha was found to perform better in range of 12–18, whereas initial gravitational constant G0 was found to perform better in the range from 102 to 120.
  • Third, compute the particle fitness value of each node. The RMSE is selected as the fitness function in this study. The fitness function of this method is defined as follows:
  • Fourth, chose the best parameter combination through GSA to obtain the optimal values for the LSSVR parameters.
  • Fifth, utilize the new combination of parameters to reconstruct the LSSVR if it does not meet the stopping criterion.
  • Last, the optimal LSSVR model for forecasting evapotranspiration was built based on the typical parameter values.

2.2.5. M5 Model Tree

M5 model tree (M5RT) is a novel neural approach developed initially by Quinlan [49] to help in solving continuous class learning problems. This model is based on a binary tree that serves as the backbone for the model itself. An M5 model tree consists of a linear regression function applied to terminal leaf nodes to describe the relationships of independent and dependent variables. Since the M5 model tree is a quantitative data focus model, this feature increases the accuracy, and as a result, the importance compared to other models [50]. In general, the principle of the tree-based model consists of “divide-and-conquer” approach for constructing a relationship between independent and dependent variables or input and output parameters. They can also be used for qualitative and quantitative data assessment [51]. All model trees can effectively learn and succeed in tasks with very high dimensionality. The main advantage of the M5 model tree compared to regression trees and other models is that it is much smaller than regression trees [52]. Additionally, the decision strength is evident, and the regression functions do not usually contain abundant variables, which in some cases can reduce the computational performance.
M5 model tree consists of two steps. The first step has to do with the splitting of the datasets into subsets. Splitting process: it often creates an excessive tree-like structure that leads to overfitting. In the second step, the overrun tree is trimmed; then trimmed subtrees are replaced with linear regression functions (Figure 3). M5 model tree has found many applications in water resource-related areas. Rahimikhoob et al. [53] applied the M5 model tree to compute evapotranspiration (ETo) in a semi-arid region. They found that the M5 model tree provides more robust results compared to other conventional or empirical methods. However, for the arid zone, Rahimikhoob [54] found that ANN estimated ET0 better than the M5 model tree.
Nevertheless, both models performed in agreement for that study case, and the results were very similar to the results obtained from the FAO56-PM method. Kisi and Kilic [52] estimated ETo considering several empirical and neural-fuzzy approaches, including the M5 model tree. They found that the M5 model tree showed better accuracy, particularly comparing to the empirical models.

3. Application and Results

The potential of new machine learning methods, LSSVR-GSA and DENFIS, was investigated by applying them to temperature (T) and extraterrestrial radiation (Ra) data of three stations located in China, and the results were compared with the well-known M5RT method. The models’ accuracy was evaluated based on three statistical criteria: root mean square error (RMSE), mean absolute error (MAE), and determination coefficient (R2). The RMSE and MAE can be defined as:
MAE =     i = 1 N E T i m E T i p N
RMSE = i = 1 N E T i m E T i p 2 N
where N is the number of datapoints used, E T ip is predicted ETo, and E T im is ETo calculated by FAO 56 PM.
Table S2 sums up the validation and test results of the and M5RT models concerning different input combinations comprising previous values of temperature (T) and extraterrestrial radiation (Ra) for the first station (56004). In Table S2, Tt and Rat refer to the mean air temperature and extraterrestrial radiation at time t and vice versa. In the tables, t represents the month which ETo needs to be predicted, whereas t-1 represents the previous month.
It is apparent from Table S2 that the LSSVR-GSA has slightly better accuracy than the DENFIS, and they both perform superiorly to the M5RT in modeling monthly ETo. There are not large differences between the optimal T-based and Ra-based models; the RMSE differences are 12.6%, 12.9%, and 5.9% for the LSSVR-GSA, DENFIS, and M5RT, respectively. It is worth noting that Ra-based models only use Julian date and do not require temperature measurement.
In Table 1, the model results are presented for the optimal (best) inputs, including periodicity input (α), which indicates the month number of the output. In the table, Opt T and Opt Ra refer the optimum T and optimum Ra inputs, which provided the best accuracy in the test stage (see Table S2); Opt T based models for LSSVR-GSA, DENFIS, and M5RT are found as Tt, Tt-1, Tt-2, Tt, Tt-1, Tt-2, Tt-3, Tt, Tt-1, and Tt-2 respectively. For the optimal Ra input combination, Rat, Rat-1, Rat-2, and Rat-3 input combination provided the best results for DENFIS and LSSVR-GSA models, whereas, Rat, Rat-1, Rat-2 provided best results for the M5RT model. As clearly observed from Table 1, the models’ performances are also examined by considering both optimum temperature and extraterrestrial radiation inputs. As evident from the table, including periodicity input marginally improves the models’ accuracy, and combining optimal T and optimal Ra inputs worsens the performances of the three methods in prediction of monthly ETo.
Model estimates of each method with Opt T and Opt Ra inputs in the test period are compared in Figure 4a,b in the form of hydrograph for the first station (56021). It can be seen from the zoomed sections that the LSSVR-GSA is closer to the FAO 56 PM ETo than the DENFIS and M5RT models. Figure 5 makes a scatterplot comparison of the model estimates. As seen from the graphs, LSSVR-GSA has less scattered points, closely followed by the DENFIS model, for both temperature and extraterrestrial radiation inputs.
Validation and test statistics of the three applied methods are reported in Table S3 for the prediction of monthly ETo of the second station (56021). It is observable from Table S3 that the temperature-based LSSVR-GSA models outperform the corresponding DENFIS and M5RT models. The best LSSVR-GSA model has lower RMSE (0.230 mm) and MAE (0.179 mm) and higher R2 (0.956) compared to the other two best models. In the case of extraterrestrial radiation input, however, LSSVR-GSA has slightly better accuracy than the DENFIS and M5RT models. It is worth noting here that Ra-based models perform superiorly regarding the T based models, except LSSVR-GSA, so that there is a marginal difference between R-based and T-based models. This result caries practical importance because R-based models only use Julian’s date, and they can predict monthly ETo without climatic data. For station 56021, the optimal T inputs were found Tt, Tt-1, and Tt-2, for LSSVR-GSA and M5RT models, whereas they were Tt, Tt-1, Tt-2, and Tt-3 for DENFIS model.
The optimal models with the periodicity information are compared in Table 2 for the second station (56021). In this station, also adding periodicity information slightly improves prediction accuracy, and in some cases, it negatively affects the models’ accuracies (e.g., see DENFIS with Opt Ra input and M5RT with Opt T and Opt T, Opt Ra inputs). Here, also combining optimal temperature and extraterrestrial radiation inputs does not improve the models’ exactness compared to optimal T and/or optimal Ra input-based models.
Figure 6a,b illustrate the ETo estimates of three methods with Opt T and Opt Ra inputs in the test period for the second station (56021). As can be seen, the LSSVR-GSA estimates are closer to the FAO 56 PM ETo than the DENFIS and M5TR models. It can be said that the fluctuations of the model estimates are very close to each other in the case of Opt Ra input (Figure 6b). This confirms the statistics provided in Table 2. Scatterplot comparison of the model estimates is shown in Figure 7. It is evident from the graphs that the LSSVR-GSA has less scattered estimates compared to the other two models and is closely followed by the DENFIS. In the case of Ra input, the distributions of all three models are similar to each other.
Table S4 gives the validation and test statistics of the LSSVR-GSA, DENFIS, and M5RT models for prediction of monthly ETo of the third station (56029) using previous values of temperature and extraterrestrial radiation to find their optimal input combinations. In this station, the LSSVR-GSA method performs better than the DENFIS and M5RT methods; however, there is a slight difference among the Ra-based methods, similarly to the previous station (56021). Among the all models, the LSSVR-GSA with Tt, Tt-1, Tt-2, and Tt-3 inputs has the lowest RMSE (0.230 mm) and MAE (0.172 mm) and the highest R2 (0.954) in monthly ETo prediction.
Table 3 compares the optimal T- and R-based models, together with periodicity input (α). For the LSSVR-GSA method including α generally, it worsens the prediction accuracy, while the performances of DENFIS models are considerably improved by importing periodicity input; the most considerable improvement is 25% for the model with Opt T input. Combining optimum T and optimum Ra inputs also considerably improve the prediction accuracy of DENFIS; improvements are 22% and 10% compared to Opt T and Opt Ra inputs. For this station, it can be concluded that the DENFIS model with Opt T, α input performs better than the other models.
Figure 8a,b provide the ETo estimates of three methods with Opt T and Opt Ra inputs in the test period for the third station (56029). In this station, the difference between LSSVR and the other two methods is more clearly observed, mainly for the optimum temperature inputs. Figure 9 compares the models’ estimates concerning their scatterplot distributions. The LSSVR-GSA has the estimates which are closer to the FAO 56 PM ETo with less scattered distribution compared to DENFIS and M5RT models. In the case of Ra input, the other two models and are closely followed by the DENFIS. In the case of Ra input, however, the distributions of all three models are very close to each other.

4. Conclusions

The potentials of two new machine learning methods, least-square support vector regression with a gravitational search algorithm and the dynamic evolving neural-fuzzy inference system, for modeling reference evapotranspiration using only temperature data, were investigated in the study, and estimates of the new methods were compared with the well-known M5 model tree method. Previous values of air temperature and extraterrestrial radiation were tried as inputs to the models to estimate monthly ETo. Three stations from China were used as case studies. The results obtained from the applications provided the following conclusions:
  • In all three stations, the temperature or extra-terrestrial radiation-based LSSVR-GSA models performed superiorly to the DENFIS and M5RT models in estimating monthly ETo. In some cases, especially for the Ra based models, however, a slight difference was found between the LSSVR-GSA and DENFIS methods, while the M5RT provided the worst estimates.
  • The results revealed that the only extra-terrestrial radiation input might provide better prediction accuracy for all the three methods, and this implies that the monthly ETo can be easily calculated with only Julian date or without temperature information.
  • Importing periodicity information to the model’s inputs generally improved the prediction accuracy, and the accuracy of DENFIS was considerably increased in the third station (56029).
  • Combining optimum air temperature and extra-terrestrial radiation inputs generally did not increase the accuracy of the methods in terms of estimation of monthly ETo. In one station (third station), however, the combination of both inputs improved the accuracy of DENFIS methods by 22% and 10% (RMSE) compared to optimum T and optimum Ra inputs, respectively. In this station, this method performed better than the LSSVR-GSA and M5RT.
The results of this study recommend the use of the LSSVR-GSA model as an efficient tool for estimating monthly ETo using only temperature or extraterrestrial radiation. This is very important in practical applications, especially for the developing countries, where some climatic data may be missing or absent because of technical reasons or lack of opportunities.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/1099-4300/22/5/547/s1. Table S1: The statistical parameters of the applied data. Table S2: Validation and test statistics of the models for monthly ET0 prediction—Station 1 (56004). Table S3: Validation and test statistics of the models for monthly ET0 prediction—Station 2 (56021). Table S4: Validation and test statistics of the models for monthly ET0 prediction—Station 3 (56029).

Author Contributions

Formal analysis, R.M.A. and O.K.; methodology, O.K. and M.I.; writing—original draft, R.M.A., A.K. and A.E.-S.; supervision, Z.C. and X.Y. All authors have read and agreed to the published version of the manuscript.

Funding

Open Fund of the Hubei Provincial Key Laboratory for Operation and Control of Cascaded Hydropower Station in China Three Gorges University: No. 2019KJX02.

Acknowledgments

This work was supported by the Open Fund of the Hubei Provincial Key Laboratory for Operation and Control of Cascaded Hydropower Station in China, Three Gorges University (2019KJX02). Alban Kuriqi was supported by a Ph.D. scholarship granted by Fundação para a Ciência e a Tecnologia, I.P. (FCT), Portugal, under the Ph.D. Program FLUVIO–River Restoration and Management, grant number: PD/BD/114558/2016.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Gavili, S.; Sanikhani, H.; Kisi, O.; Mahmoudi, M.H. Evaluation of several soft computing methods in monthly evapotranspiration modelling. Meteorol. Appl. 2018, 25, 128–138. [Google Scholar] [CrossRef] [Green Version]
  2. Sanikhani, H.; Kisi, O.; Maroufpoor, E.; Yaseen, Z.M. Temperature-based modeling of reference evapotranspiration using several artificial intelligence models: Application of different modeling scenarios. Theor. Appl. Climatol. 2019, 135, 449–462. [Google Scholar] [CrossRef]
  3. Almorox, J.; Senatore, A.; Quej, V.H.; Mendicino, G. Worldwide assessment of the Penman–Monteith temperature approach for the estimation of monthly reference evapotranspiration. Theor. Appl. Climatol. 2018, 131, 693–703. [Google Scholar] [CrossRef]
  4. Kumar, M.; Raghuwanshi, N.S.; Singh, R.; Wallender, W.W.; Pruitt, W.O. Estimating evapotranspiration using artificial neural network. J. Irrig. Drain. Eng. 2002, 128, 224–233. [Google Scholar] [CrossRef]
  5. Tie, Q.; Hu, H.; Tian, F.; Holbrook, N.M. Comparing different methods for determining forest evapotranspiration and its components at multiple temporal scales. Sci. Total Environ. 2018, 633, 12–29. [Google Scholar] [CrossRef]
  6. Li, G.; Zhang, F.; Jing, Y.; Liu, Y.; Sun, G. Response of evapotranspiration to changes in land use and land cover and climate in China during 2001–2013. Sci. Total Environ. 2017, 596, 256–265. [Google Scholar] [CrossRef]
  7. Gharbia, S.S.; Smullen, T.; Gill, L.; Johnston, P.; Pilla, F. Spatially distributed potential evapotranspiration modeling and climate projections. Sci. Total Environ. 2018, 633, 571–592. [Google Scholar] [CrossRef]
  8. Keshtegar, B.; Kisi, O.; Arab, H.G.; Zounemat-Kermani, M. Subset modeling basis ANFIS for prediction of the reference evapotranspiration. Water Resour. Res. 2018, 32, 1101–1116. [Google Scholar]
  9. Malik, A.; Kumar, A.; Kisi, O. Daily pan evaporation estimation using heuristic methods with gamma test. J. Irrig. Drain. Eng. 2018, 144, 04018023. [Google Scholar] [CrossRef]
  10. Malik, A.; Kumar, A. Pan evaporation simulation based on daily meteorological data using soft computing techniques and multiple linear regression. Water Resour. Res. 2015, 29, 1859–1872. [Google Scholar] [CrossRef]
  11. Adnan, R.M.; Malik, A.; Kumar, A.; Parmar, K.S.; Kisi, O. Pan evaporation modeling by three different neuro-fuzzy intelligent systems using climatic inputs. Arab. J. Geosci. 2019, 12, 606. [Google Scholar] [CrossRef]
  12. Yuan, X.; Chen, C.; Lei, X.; Yuan, Y.; Adnan, R.M. Monthly runoff forecasting based on LSTM–ALO model. Stoch. Env. Res. Risk A. 2018, 32, 2199–2212. [Google Scholar] [CrossRef]
  13. Kisi, O.; Shiri, J.; Karimi, S.; Adnan, R.M. Three different adaptive neuro fuzzy computing techniques for forecasting long-period daily streamflows. In Big Data in Engineering Applications; Springer: Singapore, Singapore, 2018; pp. 303–321. [Google Scholar]
  14. Guven, A.; Kişi, Ö. Daily pan evaporation modeling using linear genetic programming technique. Irrig. Sci. 2011, 29, 135–145. [Google Scholar] [CrossRef]
  15. Kim, S.; Seo, Y.; Singh, V.P. Assessment of pan evaporation modeling using bootstrap resampling and soft computing methods. J. Comput. Civ. Eng. 2013, 29, 04014063. [Google Scholar] [CrossRef]
  16. El-Shafie, A.; Alsulami, H.M.; Jahanbani, H.; Najah, A. Multi-lead ahead prediction model of reference evapotranspiration utilizing ANN with ensemble procedure. Stoch. Env. Res. Risk A. 2013, 27, 1423–1440. [Google Scholar] [CrossRef]
  17. Chen, Q.; Mynett, A.E. Predicting Phaeocystis globosa bloom in Dutch coastal waters by decision trees and nonlinear piecewise regression. Ecol. Model. 2004, 176, 277–290. [Google Scholar] [CrossRef]
  18. Adnan, R.M.; Yuan, X.; Kisi, O.; Anam, R. Improving accuracy of river flow forecasting using LSSVR with gravitational search algorithm. Adv. Meteorol. 2017. [Google Scholar] [CrossRef]
  19. Kisi, O. Pan evaporation modeling using least square support vector machine, multivariate adaptive regression splines and M5 model tree. J. Hydrol. 2015, 528, 312–320. [Google Scholar] [CrossRef]
  20. Adnan, R.M.; Liang, Z.; Yuan, X.; Kisi, O.; Akhlaq, M.; Li, B. Comparison of LSSVR, M5RT, NF-GP, and NF-SC Models for Predictions of Hourly Wind Speed and Wind Power Based on Cross-Validation. Energies 2019, 12, 329. [Google Scholar] [CrossRef] [Green Version]
  21. Adnan, R.M.; Yuan, X.; Kisi, O.; Adnan, M.; Mehmood, A. Stream Flow Forecasting of Poorly Gauged Mountainous Watershed by Least Square Support Vector Machine, Fuzzy Genetic Algorithm and M5 Model Tree Using Climatic Data from Nearby Station. Water Resour. Manag. 2018, 32, 4469–4486. [Google Scholar] [CrossRef]
  22. Terzi, Ö.; Erol Keskin, M.; Dilek Taylan, E. Estimating evaporation using ANFIS. J. Irrig. Drain. Eng. 2006, 132, 503–507. [Google Scholar] [CrossRef]
  23. Moghaddamnia, A.; Gousheh, M.G.; Piri, J.; Amin, S.; Han, D. Evaporation estimation using artificial neural networks and adaptive neuro-fuzzy inference system techniques. Adv. Water Resour. 2009, 32, 88–97. [Google Scholar] [CrossRef]
  24. Shiri, J.; Dierickx, W.; Pour-Ali Baba, A.; Neamati, S.; Ghorbani, M.A. Estimating daily pan evaporation from climatic data of the State of Illinois, USA using adaptive neuro-fuzzy inference system (ANFIS) and artificial neural network (ANN). Hydrol. Res. 2011, 42, 491–502. [Google Scholar] [CrossRef]
  25. Muhammad Adnan, R.; Yuan, X.; Kisi, O.; Yuan, Y.; Tayyab, M.; Lei, X. Application of soft computing models in streamflow forecasting. In Proceedings of the Institution of Civil Engineers-Water Management; Thomas Telford Ltd.: London, UK, 2019; Volume 172, pp. 123–134. [Google Scholar]
  26. Sanikhani, H.; Kisi, O. River flow estimation and forecasting by using two different adaptive neuro-fuzzy approaches. Water Resour. Manag. 2012, 26, 1715–1729. [Google Scholar] [CrossRef]
  27. Heddam, S.; Dechemi, N. A new approach based on the dynamic evolving neural-fuzzy inference system (DENFIS) for modelling coagulant dosage (Dos): Case study of water treatment plant of Algeria. Desalin. Water Treat. 2015, 53, 1045–1053. [Google Scholar]
  28. Kasabov, N.K.; Song, Q. DENFIS: Dynamic evolving neural-fuzzy inference system and its application for time-series prediction. IEEE T. Fuzzy Sys. 2002, 10, 144–154. [Google Scholar] [CrossRef] [Green Version]
  29. Adnan, R.M.; Liang, Z.; El-Shafie, A.; Zounemat-Kermani, M.; Kisi, O. Prediction of Suspended Sediment Load Using Data-Driven Models. Water 2019, 11, 2060. [Google Scholar] [CrossRef] [Green Version]
  30. Rashedi, E.; Nezamabadi-Pour, H.; Saryazdi, S. GSA: A gravitational search algorithm. Inf. Sci. 2009, 179, 2232–2248. [Google Scholar] [CrossRef]
  31. Tao, H.; Diop, L.; Bodian, A.; Djaman, K.; Ndiaye, P.M.; Yaseen, Z.M. Reference evapotranspiration prediction using hybridized fuzzy model with firefly algorithm: Regional case study in Burkina Faso. Agric. Water Manag. 2018, 208, 140–151. [Google Scholar] [CrossRef]
  32. Adnan, R.M.; Liang, Z.; Trajkovic, S.; Zounemat-Kermani, M.; Li, B.; Kisi, O. Daily streamflow prediction using optimally pruned extreme learning machine. J. Hydrol. 2019, 577, 123981. [Google Scholar] [CrossRef]
  33. Mosavi, A.; Edalatifar, M. A hybrid neuro-fuzzy algorithm for prediction of reference evapotranspiration. In International Conference on Global Research and Education; Springer: Cham, Switzerland, 2018; pp. 235–243. [Google Scholar]
  34. Wu, L.; Zhou, H.; Ma, X.; Fan, J.; Zhang, F. Daily reference evapotranspiration prediction based on hybridized extreme learning machine model with bio-inspired optimization algorithms: Application in contrasting climates of China. J. Hydrol. 2019, 577, 123960. [Google Scholar] [CrossRef]
  35. Huang, G.; Wu, L.; Ma, X.; Zhang, W.; Fan, J.; Yu, X.; Zeng, W.; Zhou, H. Evaluation of CatBoost method for prediction of reference evapotranspiration in humid regions. J. Hydrol. 2019, 574, 1029–1041. [Google Scholar] [CrossRef]
  36. Allen, R.G.; Pereira, L.S.; Raes, D.; Smith, M. Crop Evapotranspiration-Guidelines for Computing Crop Water Requirements; FAO Irrigation and Drainage Paper 56; FAO: Rome, Italy, 1998; Volume 300. [Google Scholar]
  37. She, D.; Xia, J.; Zhang, Y. Changes in reference evapotranspiration and its driving factors in the middle reaches of Yellow River Basin, China. Sci. Total Environ. 2017, 607, 1151–1162. [Google Scholar] [CrossRef] [PubMed]
  38. Talei, A.; Chua, L.H.C.; Quek, C.; Jansson, P.E. Runoff forecasting using a Takagi–Sugeno neuro-fuzzy model with online learning. J. Hydrol. 2013, 488, 17–32. [Google Scholar] [CrossRef]
  39. Kisi, O.; Mansouri, I.; Hu, J.W. A New Method for Evaporation Modeling: Dynamic Evolving Neural-Fuzzy Inference System. Adv. Meteorol. 2017. [Google Scholar] [CrossRef] [Green Version]
  40. Heddam, S.; Watts, M.J.; Houichi, L.; Djemili, L.; Sebbar, A. Evolving connectionist systems (ECoSs): A new approach for modeling daily reference evapotranspiration (ET 0). Environ. Monit. Assess. 2018, 190, 516. [Google Scholar] [CrossRef]
  41. Kisi, O.; Heddam, S.; Yaseen, Z.M. The implementation of univariable scheme-based air temperature for solar radiation prediction: New development of dynamic evolving neural-fuzzy inference system model. Appl. Energy 2019, 241, 184–195. [Google Scholar] [CrossRef]
  42. Suykens, J.A.; Vandewalle, J. Recurrent least squares support vector machines. In IEEE Transactions on Circuits and Systems I: Fundamental Theory and Applications; IEEE: Piscataway, NJ, USA, 2000; Volume 47, pp. 1109–1114. [Google Scholar]
  43. Wu, Q.; Peng, C. Wind power grid connected capacity prediction using LSSVM optimized by the bat algorithm. Energies 2015, 8, 14346–14360. [Google Scholar] [CrossRef] [Green Version]
  44. Yuan, X.; Chen, C.; Yuan, Y.; Huang, Y.; Tan, Q. Short-term wind power prediction based on LSSVM–GSA model. Energy Conver. Manag. 2015, 101, 393–401. [Google Scholar] [CrossRef]
  45. Gu, Y.; Zhao, W.; Wu, Z. Online adaptive least squares support vector machine and its application in utility boiler combustion optimization systems. J. Process Control 2011, 21, 1040–1048. [Google Scholar] [CrossRef]
  46. Lu, P.; Ye, L.; Sun, B.; Zhang, C.; Zhao, Y.; Teng, J. A new hybrid prediction method of ultra-short-term wind power forecasting based on EEMD-PE and LSSVM optimized by the GSA. Energies 2018, 11, 697. [Google Scholar] [CrossRef] [Green Version]
  47. Sun, W.; Sun, J. Daily PM2. 5 concentration prediction based on principal component analysis and LSSVM optimized by cuckoo search algorithm. J. Environ. Manag. 2017, 188, 144–152. [Google Scholar] [CrossRef] [PubMed]
  48. Sahu, R.K.; Panda, S.; Padhan, S. A novel hybrid gravitational search and pattern search algorithm for load frequency control of nonlinear power system. Appl. Soft Comput. 2015, 29, 310–327. [Google Scholar] [CrossRef]
  49. Quinlan, J.R. Learning with continuous classes. In Proceedings of the 5th Australian Joint Conference on Artificial Intelligence; World Scientific: Singapore, Singapore, 1992; pp. 343–348. [Google Scholar]
  50. Adnan, R.M.; Liang, Z.; Heddam, S.; Zounemat-Kermani, M.; Kisi, O.; Li, B. Least square support vector machine and multivariate adaptive regression splines for streamflow prediction in mountainous basin using hydro-meteorological data as inputs. J. Hydrol. 2019, 124371. [Google Scholar] [CrossRef]
  51. Alipour, A.; Yarahmadi, J.; Mahdavi, M. Comparative study of M5 model tree and artificial neural network in estimating reference evapotranspiration using MODIS products. J. Climatol. 2014. [Google Scholar] [CrossRef] [Green Version]
  52. Kisi, O.; Kilic, Y. An investigation on generalization ability of artificial neural networks and M5 model tree in modeling reference evapotranspiration. Theor. Appl. Climatol. 2016, 126, 413–425. [Google Scholar] [CrossRef]
  53. Rahimikhoob, A.; Asadi, M.; Mashal, M. A comparison between conventional and M5 model tree methods for converting pan evaporation to reference evapotranspiration for semi-arid region. Water Resour. Manag. 2013, 27, 4815–4826. [Google Scholar] [CrossRef]
  54. Rahimikhoob, A. Comparison between M5 model tree and neural networks for estimating reference evapotranspiration in an arid environment. Water Resour. Manag. 2014, 28, 657–669. [Google Scholar] [CrossRef]
Figure 1. Study area.
Figure 1. Study area.
Entropy 22 00547 g001
Figure 2. LSSVR model for evapotranspiration (ETo) modeling.
Figure 2. LSSVR model for evapotranspiration (ETo) modeling.
Entropy 22 00547 g002
Figure 3. Schematic view of M5 model tree (a) structure and (b) splitting data space into sub-regions.
Figure 3. Schematic view of M5 model tree (a) structure and (b) splitting data space into sub-regions.
Entropy 22 00547 g003
Figure 4. (a) Time variation graphs of the FAO 56 PM and estimated ETo by LSSVR-GSA, DENFIS, and M5RT in the test period of Station 56004 using optimal T inputs. (b) Time variation graphs of the FAO 56 PM and estimated ETo by LSSVR-GSA, DENFIS, and M5RT in the test period of Station 56004 using optimal Ra inputs.
Figure 4. (a) Time variation graphs of the FAO 56 PM and estimated ETo by LSSVR-GSA, DENFIS, and M5RT in the test period of Station 56004 using optimal T inputs. (b) Time variation graphs of the FAO 56 PM and estimated ETo by LSSVR-GSA, DENFIS, and M5RT in the test period of Station 56004 using optimal Ra inputs.
Entropy 22 00547 g004aEntropy 22 00547 g004b
Figure 5. Time scatterplots of the observed and estimated ETo (by LSSVR-GSA, DENFIS, and M5RT) in the test period of Station 56029 using optimal (a) T and (b) Ra inputs.
Figure 5. Time scatterplots of the observed and estimated ETo (by LSSVR-GSA, DENFIS, and M5RT) in the test period of Station 56029 using optimal (a) T and (b) Ra inputs.
Entropy 22 00547 g005
Figure 6. (a) Time variation graphs of the FAO 56 PM and estimated ETo (by LSSVR-GSA, DENFIS, and M5RT) in the test period of Station 56021 using optimal T inputs. (b) Time variation graphs of the FAO 56 PM and estimated ETo by LSSVR-GSA, DENFIS, and M5RT in the test period of Station 56021 using optimal Ra inputs.
Figure 6. (a) Time variation graphs of the FAO 56 PM and estimated ETo (by LSSVR-GSA, DENFIS, and M5RT) in the test period of Station 56021 using optimal T inputs. (b) Time variation graphs of the FAO 56 PM and estimated ETo by LSSVR-GSA, DENFIS, and M5RT in the test period of Station 56021 using optimal Ra inputs.
Entropy 22 00547 g006
Figure 7. Time scatterplots of the observed and estimated ETo (by LSSVR-GSA, DENFIS, and M5RT) in the test period of Station 56021 using optimal (a) T and (b) Ra inputs.
Figure 7. Time scatterplots of the observed and estimated ETo (by LSSVR-GSA, DENFIS, and M5RT) in the test period of Station 56021 using optimal (a) T and (b) Ra inputs.
Entropy 22 00547 g007
Figure 8. (a)Time variation graphs of the FAO 56 PM and estimated ETo (by LSSVR-GSA, DENFIS, and M5RT) in the test period of Station 56029 using optimal T inputs. (b) Time variation graphs of the FAO 56 PM and estimated ETo by LSSVR-GSA, DENFIS, and M5RT in the test period of Station 56029 using optimal Ra inputs.
Figure 8. (a)Time variation graphs of the FAO 56 PM and estimated ETo (by LSSVR-GSA, DENFIS, and M5RT) in the test period of Station 56029 using optimal T inputs. (b) Time variation graphs of the FAO 56 PM and estimated ETo by LSSVR-GSA, DENFIS, and M5RT in the test period of Station 56029 using optimal Ra inputs.
Entropy 22 00547 g008aEntropy 22 00547 g008b
Figure 9. Time scatterplots of the observed and estimated ETo (by LSSVR-GSA, DENFIS, and M5RT) in the test period of Station 56021 using optimal (a) T and (b) Ra inputs.
Figure 9. Time scatterplots of the observed and estimated ETo (by LSSVR-GSA, DENFIS, and M5RT) in the test period of Station 56021 using optimal (a) T and (b) Ra inputs.
Entropy 22 00547 g009
Table 1. Validation and test statistics of the models in monthly ETo prediction using the optimal inputs of T, Ra, and α—Station 1 (56004).
Table 1. Validation and test statistics of the models in monthly ETo prediction using the optimal inputs of T, Ra, and α—Station 1 (56004).
Model InputsValidation PeriodTest Period
RMSEMAER2RMSEMAER2
LSSVM-GSA
Opt T0.2460.1840.9580.2460.1820.953
Opt T, α0.2430.1800.9590.2440.1800.955
Opt Ra0.3410.2460.9200.2770.2190.935
Opt Ra, α0.3390.2440.9220.2740.2170.936
Opt T, Opt Ra0.2850.2100.9400.3010.2260.930
Opt T, Opt Ra, α0.2840.2090.9420.2990.2250.932
DENFIS
Opt T0.2710.1980.9470.2490.1900.950
Opt T, α0.2630.1930.9500.2460.1850.952
Opt Ra0.3420.2480.9190.2810.2210.934
Opt Ra, α0.3410.2500.9200.2800.2210.935
Opt T, Opt Ra0.2880.2120.9390.3050.2290.927
Opt T, Opt Ra, α0.2860.2110.9400.3040.2280.928
M5RT
Opt T0.3220.2210.9230.3050.2320.928
Opt T, α0.3040.2150.9310.3230.2390.931
Opt Ra0.3500.2660.9170.3230.2670.932
Opt Ra, α0.3480.2640.9180.3220.2640.933
Opt T, Opt Ra0.3040.2190.9310.3090.2310.923
Opt T, Opt Ra, α0.3010.2150.9330.3060.2300.925
Table 2. Validation and test statistics of the models in monthly ETo prediction using the optimal inputs of T, Ra, and α—Station 2 (56021).
Table 2. Validation and test statistics of the models in monthly ETo prediction using the optimal inputs of T, Ra, and α—Station 2 (56021).
Model InputsValidation PeriodTest Period
RMSEMAER2RMSEMAER2
LSSVM-GSA
Opt T0.2350.1610.9540.2300.1790.956
Opt T, α0.2330.1580.9550.2280.1760.961
Opt Ra0.2760.1950.9330.2360.1760.949
Opt Ra, α0.2730.1920.9350.2340.1750.950
Opt T, Opt Ra0.2450.1700.9440.2380.1950.940
Opt T, Opt Ra, α0.2360.1640.9500.2320.1820.952
DENFIS
Opt T0.2420.1620.9510.3010.2530.945
Opt T, α0.2400.1620.9520.2270.1710.961
Opt Ra0.2860.2090.9320.2410.1860.947
Opt Ra, α0.2880.2130.9300.2680.2070.942
Opt T, Opt Ra0.2760.1870.9440.2840.2150.938
Opt T, Opt Ra, α0.2600.1780.9450.2690.2080.940
M5RT
Opt T0.2650.1890.9400.3100.2270.921
Opt T, α0.3030.2200.9290.3410.2600.908
Opt Ra0.2890.2110.9290.2510.1910.942
Opt Ra, α0.2810.2050.9310.2500.1850.943
Opt T, Opt Ra0.3000.2160.9300.3360.2580.905
Opt T, Opt Ra, α0.3060.2200.9250.3390.2610.901
Table 3. Validation and test statistics of the models in monthly ETo prediction using the optimal inputs of T, Ra, and α—Station 3 (56029).
Table 3. Validation and test statistics of the models in monthly ETo prediction using the optimal inputs of T, Ra, and α—Station 3 (56029).
Model InputsValidation PeriodTest Period
RMSEMAER2RMSEMAER2
LSSVM-GSA
Opt T0.2020.1600.9680.2300.1720.954
Opt T, α0.1990.1570.9700.2910.2400.957
Opt Ra0.2300.1790.9570.2620.2050.947
Opt Ra, α0.2280.1760.9580.2600.2020.949
Opt T, Opt Ra0.2080.1610.9660.2980.2430.959
Opt T, Opt Ra, α0.2080.1610.9660.3010.2460.959
DENFIS
Opt T0.2170.1600.9640.2910.2390.951
Opt T, α0.2220.1660.9620.2180.1650.959
Opt Ra0.2310.1780.9570.2680.2090.944
Opt Ra, α0.2420.1900.9570.2520.1920.947
Opt T, Opt Ra0.2260.1690.9600.2260.1670.950
Opt T, Opt Ra, α0.2210.1660.9610.2190.1580.952
M5RT
Opt T0.2760.2090.9380.3170.2240.910
Opt T, α0.2700.2050.9420.3150.2200.912
Opt Ra0.2320.1930.9540.2720.2150.940
Opt Ra, α0.2280.1800.9560.2690.2100.943
Opt T, Opt Ra0.2630.1950.9440.3020.2130.920
Opt T, Opt Ra, α0.2600.1930.9460.2990.2100.922

Share and Cite

MDPI and ACS Style

Muhammad Adnan, R.; Chen, Z.; Yuan, X.; Kisi, O.; El-Shafie, A.; Kuriqi, A.; Ikram, M. Reference Evapotranspiration Modeling Using New Heuristic Methods. Entropy 2020, 22, 547. https://0-doi-org.brum.beds.ac.uk/10.3390/e22050547

AMA Style

Muhammad Adnan R, Chen Z, Yuan X, Kisi O, El-Shafie A, Kuriqi A, Ikram M. Reference Evapotranspiration Modeling Using New Heuristic Methods. Entropy. 2020; 22(5):547. https://0-doi-org.brum.beds.ac.uk/10.3390/e22050547

Chicago/Turabian Style

Muhammad Adnan, Rana, Zhihuan Chen, Xiaohui Yuan, Ozgur Kisi, Ahmed El-Shafie, Alban Kuriqi, and Misbah Ikram. 2020. "Reference Evapotranspiration Modeling Using New Heuristic Methods" Entropy 22, no. 5: 547. https://0-doi-org.brum.beds.ac.uk/10.3390/e22050547

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop