Next Article in Journal
Evaluating Non-Composted Red Cotton Tree (Bombax ceiba) Sawdust Mixtures for Raising Okra (Abelmoschus esculentus (L.) Moench) in Pots
Previous Article in Journal
Short-term Response of Greenhouse Gas Emissions from Precision Fertilization on Barley
Previous Article in Special Issue
A Revised Equation of Water Application Efficiency in a Center Pivot System Used in Crop Rotation in No Tillage
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Advanced Hybrid Metaheuristic Machine Learning Models Application for Reference Crop Evapotranspiration Prediction

by
Rana Muhammad Adnan Ikram
1,
Reham R. Mostafa
2,
Zhihuan Chen
3,*,
Abu Reza Md. Towfiqul Islam
4,*,
Ozgur Kisi
5,6,
Alban Kuriqi
7,8 and
Mohammad Zounemat-Kermani
9
1
School of Economics and Statistics, Guangzhou University, Guangzhou 510006, China
2
Information Systems Department, Faculty of Computers and Information Sciences, Mansoura University, Mansoura 35516, Egypt
3
Engineering Research Center for Metallurgical Automation and Measurement Technology of Ministry of Education, Wuhan University of Science and Technology, Wuhan 431400, China
4
Department of Disaster Management, Begum Rokeya University, Rangpur 5400, Bangladesh
5
Department of Civil Engineering, Technical University of Lübeck, 23562 Lübeck, Germany
6
Civil Engineering Department, Ilia State University, 0162 Tbilisi, Georgia
7
CERIS, Instituto Superior Técnico, Universidade de Lisboa, 53089 Lisbon, Portugal
8
Civil Engineering Department, University for Business and Technology, 10000 Pristina, Kosovo
9
Department of Water Engineering, Shahid Bahonar University of Kerman, Kerman 93630, Iran
*
Authors to whom correspondence should be addressed.
Agronomy 2023, 13(1), 98; https://doi.org/10.3390/agronomy13010098
Submission received: 9 November 2022 / Revised: 19 December 2022 / Accepted: 24 December 2022 / Published: 28 December 2022
(This article belongs to the Special Issue Modernization and Optimization of Irrigation Systems)

Abstract

:
Hybrid metaheuristic algorithm (MA), an advanced tool in the artificial intelligence field, provides precise reference evapotranspiration (ETo) prediction that is highly important for water resource availability and hydrological studies. However, hybrid MAs are quite scarcely used to predict ETo in the existing literature. To this end, the prediction abilities of two support vector regression (SVR) models coupled with three types of MAs including particle swarm optimization (PSO), grey wolf optimization (GWO), and gravitational search algorithm (GSA) were studied and compared with single SVR and SVR-PSO in predicting monthly ETo using meteorological variables as inputs. Data obtained from Rajshahi, Bogra, and Rangpur stations in the humid region, northwestern Bangladesh, was used for this purpose as a case study. The prediction precision of the proposed models was trained and tested using nine input combinations and assessed using root mean square error (RMSE), mean absolute error (MAE), and Nash–Sutcliffe efficiency (NSE). The tested results revealed that the SVR-PSOGWO model outperformed the other applied soft computing models in predicting ETo in all input combinations, followed by the SVR-PSOGSA, SVR-PSO, and SVR. It was found that SVR-PSOGWO decreases the RMSE of SVR, SVR-PSO, and SVR-PSOGSA by 23%, 27%, 14%, 21%, 19%, and 5% in Rangpur and Bogra stations during the testing stage. The RMSE of the SVR, SVR-PSO, and SVR-PSOGSA reduced by 32%, 20%, and 3%, respectively, employing the SVR-PSOGWO for the Rajshahi Station. The proposed hybrid machine learning model has been recommended as a potential tool for monthly ETo prediction in a humid region and similar climatic regions worldwide.

1. Introduction

Strategies for saving water in agriculture and enhancing water use efficiency are vital in water-scarce regions. The increasing global shortage of water resources and high irrigation costs requires the development of precise water-saving irrigation strategies that can minimize water use in crop production [1]. Period information on the soil water status, crop water requirements, crop water stress status, and potential yield reduction under water-stressed conditions is crucial to optimize water and energy use and maximize profits. However, accurate information about the water balance is crucial for effective agriculture management [2]. One of the most important parameters needed to estimate the water balance in each ecosystem is evapotranspiration (ET). Evapotranspiration is an essential parameter of the hydrological cycle process in natural ecosystems that links land surface water and energy balance with the atmosphere [3]. ETo is expressed as the ET rate from a reference crop surface, where the reference crop surface is a theoretical grass or alfalfa with accurate and recognized characteristics. ETo plays a crucial role in water resource availability and stimulating the hydrological effect of climate change [4]. Accurate estimation of ETo is essential for climate change predictions, drought prediction and monitoring, evaluation of water resources scarcity/availability, assessing crop water requirements, and irrigation scheduling, among others [5]. The ETo rate and quantity change from region to region depending on meteorological factors such as temperature, solar radiation, air humidity, and wind speed [6]. ETo rate directly affects agriculture’s sustainable development; on the other side, crop water requirement is affected by environmental conditions during the growing season. Each crop has a specific growing season which depends on its special growing requirements and climate conditions [7]. Farmers often experiment with different agronomic practices, irrigation and water management techniques, and crop protection practices [8]. However, irrigation scheduling, which concerns the farmers’ decision on “how much” water to apply and “when” to irrigate to maximize crop productivity, may not be carried out correctly in many cases since it is a very complex decision-making process requiring specific knowledge on crop water requirements and water budget analysis; the latter one is highly influenced by ETo rate [9,10]. Therefore, recording continuous and long-term ETo measurements, especially in agricultural areas, is a consequential issue.
When continuous field observation is not possible due to their high cost, complicated installation of the instruments, and/or exhaustive maintenance, the estimation of ETo represents a feasible way of characterizing water variability and availability [11]. There are many methods to estimate ETo, such as the water balance, the energy balance, and several physical-based hydrological models [12]. Namely, ETo estimation methods may be segregated into five main categories: pan evaporation-based, mass-transfer-based, temperature-based, radiation-based, and combined [13]. Many researchers attempted to develop methods that consistently estimate ETo using various meteorological data [5]. Tejada et al. [14] developed models based on support vector machines (SVMs) and extreme learning machines (ELMs) for the estimation of daily ETo using different input combinations of meteorological data. Their findings indicate that the SVM and ELM models, with at least Tmax, Tmin, and Rs as inputs, provide the best daily ETo estimation.
Nevertheless, most of these methods are empirical approaches calibrated based on local field observations. It should be noted that in many reports, the accuracy of the outputs provided by these methods varied among them. It depended on many factors, such as data requirements, assumptions, complexity, and reliability [15].
The water balance approach has been used as a reference method to estimate ETo. It is worth mentioning that the FAO Penman–Monteith method is widely recognized as a standard reference method for ETo estimation. Nevertheless, this model and, generally, the water balance-based assessment methods require detailed and long-term meteorological data, which are not always available everywhere [5,12]. Thus, alternative methods for predicting ETo at different temporal and spatial scales should be developed, which are easily applied and require fewer input data without jeopardizing estimation accuracy. In this regard, especially during the last decades, artificial intelligence (AI) approaches are becoming very popular in estimating ETo [16]. Dou and Yang [17] compared daily ET estimated using extreme learning machine (ELM) and an adaptive neuro-fuzzy inference system (ANFIS) with observation data. They found that ELM and ANFIS models provided robust results and could even complement the traditional methods.
Similarly, Nourani et al. [18] estimated ETo using traditional and artificially-based algorithms; they also suggest that AI methods perform better than empirical methods. Antonopoulos and Antonopoulos [19] estimated ETo using artificial neural networks (ANN) and several empirical methods. They found that ANN could provide similar results with the empirical methods while using fewer input data. Although remarkable results were achieved using AI, most studies applied simple models, which have some drawbacks, such as overfitting and low performance [20]. Integrated or so-called hybrid AI approaches have been demonstrated to be more accurate in many hydrological processes’ computations [21].
Nevertheless, to our knowledge, the hybrid AI approaches are quite scarcely applied to estimate ETo [21]. For instance, Mehdizadeh [22] applied Multivariate Adaptive Regression Splines (MARS), and Gene Expression Programming (GEP) models combined with Autoregressive Conditional Heteroscedasticity (ARCH) to estimate ETo in regions characterized by different climate conditions. They concluded that both hybrid approaches improved ETo estimation significantly. Tikhamarine et al. [23] estimated ETo using three hybrid AI models, ANN, and several empirical methods. They found that the new hybrid AI model, i.e., support vector regression (SVR) integrated with grey wolf optimizer (SVR-GWO), performed better than ANN and provided robust results. Insightful findings were reported in other studies that applied hybrid AI to estimate ETo, such as in Tikhamarine et al. [24] and Seifi and Riahi [25]. Therefore, the primary objective of this work was to evaluate the predictive performance of three hybrid AI models in a study case in humid climate conditions. Secondly, this study intends to examine the prediction performance of novel hybrid methods for a data-scarce humid region. Namely, in this study, we evaluated the performance of the following models: support vector regression (SVR) which was integrated into particle swarm optimization (PSO); PSO which was integrated into grey wolf optimizer (GWO); and gravitational search algorithm (GSA) to improve the computation performance. Thus, to our knowledge, PSO-GWO and PSO-GSA are applied for the first time in this study to estimate ETo in data-scarce humid regions.

2. Case Study and Data Description

2.1. Case Study Area

The northwestern (NW) hydrological area in Bangladesh’s drought-prone, data-scarce region was chosen as the case study site to assess the accuracy of the aforementioned models ‘predictions. Bangladesh’s northwestern region accounts for 23.5% (34,515 km2) of the country’s overall geographical area. It is subdivided into 16 administrative districts, with Rajshahi, Bogra, and Rangpur selected as case studies for this study (Figure 1). The study area is within 23°47′ N to 25°50′ N latitude and 88°01′ E to 89°48′ E longitude and located on the west of the Brahmaputra River in Bangladesh and north of the Padma River [26]. The population density of the study area is 930 people per square kilometer. Several rivers traverse the region, which Plio-Pleistocene characterizes, and Holocene oxidized red sediments and soils [27]. Figure 1 depicts the northwestern region’s geography and the sites of meteorological stations. The case study region consists of three unique physiographic units: the Barind Tract in the Rajshahi district, the channel–floodplain complexes in the Borga district, and the Himalayan piedmont plain in the Rangpur district [28]. The height of the Barind Tract ranges between 11 and 48 m above mean sea level (AMSL).
In contrast, the channel–floodplain complexes are between 8 and 23 m, and the Himalayan piedmont plain is between 30 and 45 m. The subsurface lithology of the northwest hydrological zone varies considerably. The Barind Tract is primarily a drought-prone water-deficit zone, which often adversely affects agricultural crop output and the inhabitants’ way of life in the studied area [29]. In the district of Bogra, the predominant soil texture is silty loam and clay loam; in Rajshahi, clay loam; and in Rangpur, silty loam.
The northwest area has a dry winter and a humid monsoon season with an unpredictable rainfall pattern. In total, 82% of precipitation occurred between May and October (monsoon season), whereas the remaining 18% occurred between November and April (dry season) [30]. Rajshahi had the lowest average annual rainfall (1428 mm), while Rangpur had the most (2262 mm). Due to the continued effects of climate change, the northwestern region has been experiencing severe occurrences such as flooding and drought, which renders this region more vulnerable to natural disasters [31].
The Rajshahi district receives the least average annual precipitation (1427 mm), while the Rangpur district receives the most (2260 mm). In the Rajshahi district, the monthly mean relative humidity ranges from 62% in March to 87% in July, with a mean of 78%, while in the Rangpur district, the monthly mean relative humidity ranges from 70% in March to 90% in July, with a mean of 65% [32]. In the northwest area, the mean reference evapotranspiration is 1309 mm, which varies spatially and temporally, with the greatest value in Rajshahi and the lowest value in the Rangpur district. Approximately 43% of the reference evapotranspiration occurred during the dry period, whereas 57% occurred during the monsoon season. Rangpur district has a tropical climate that is hot, rainy, and humid [33]. It has a varied monsoon season with a mean annual temperature of 26 °C, ranging from 19 °C in January to 31 °C in August. Numerous monthly variations in meteorological indices indicate that the climate of the Rangpur area is continually evolving [33]. During the summer, some days in the Rajshahi district reach temperatures as high as 45 °C.
In comparison, the temperature decreases up to 5 °C during the winter. The annual mean maximum and lowest temperatures in the Bogra district are 34.6 °C and 11.9 °C, respectively, with a total annual precipitation of 1610 mm. The rate of average temperature increase is projected to be 5.39 °C by the end of the 21st century in the northwestern area, accompanied by a 0.66 mm drop in average precipitation [34]. The elevation of most parts of Bangladesh ranged from 1 to 37.5 m above the mean sea level, which makes for a mainly low-lying “delta-shaped” landform type. Bangladesh, on the other hand, has had trouble obtaining long-term and complete climate datasets, as well as dealing with similar problems caused by nature, such as a complex hydrogeologic and climate system, and by people, such as low economic growth, a lack of good data, and technological issues.

2.2. Data Sources and Quality Control

The Bangladesh Meteorological Department (BMD) operates just 43 meteorological stations around the country, with the majority of them located in the country’s south-eastern regions. BMD only has six meteorological sites in the country’s northwestern area (www.bmd.gov.bd, accessed on 23 December 2022). These six meteorological sites are dispersed irregularly. Because three out of six stations were constructed after 1990, these stations lack the necessary long-term data for determining ETo. When multiple climatic variables were needed, a daily meteorological record from fewer stations was available. Due to these limitations of long-term daily meteorological datasets, only three sites were selected to estimate monthly ETo from 1980 to 2017 for 37 years. These three sites have no data gaps or inconsistencies in the datasets. The three chosen meteorological stations of Rajshahi, Bogra, and Rangpur districts reflect the northwestern United States hydrological area. The recorded daily minimum and maximum temperature (°C), mean relative humidity (%), wind speed (Knots), and sunlight (hours per day−1) values for the three chosen stations were retrieved from the BMD. The BMD has followed World Meteorological Organization rules while collecting and documenting meteorological datasets (WMO) [35].
Before commencing the estimate of the ETo, all datasets were evaluated for quality. Careful quality control of the obtained datasets was undertaken by verifying that all parameters had positive values; for instance, Tmin is less than Tmax, and humidity is less than 100%. The homogeneity tests of the datasets were carried out using the one-way ANOVA test, and the findings indicated that all datasets are significant at a level of 95% (p < 0.05). All datasets were additionally validated by the BMD’s professional and competent quality assurance team.
Table 1 provides a concise geographical and meteorological description of the chosen sites. Meteorological stations cannot directly determine extraterrestrial radiation (Ra), wind speed at the height of 2 m (U2), or global solar radiation (Rs). For instance, the sunshine duration values (n/N) are 0.753 for Bogura, 0.7505 for Rajshahi, and 0.7192 for Rangpur, respectively. The actual height Z of U2 measurements is 18 m for Bogura, 20 m for Rajshahi, and 34 m for Rangpur. We approximated daily Rs, Ra, and U2 using the methods proposed by Allen et al. (1998) and the available meteorological data. Although Allen et al. [36] recommended a = 0.25 and b = 0.5 (or “as” and “bs” as used for this research), these values truly vary with location. The values of these constants in the radiation equation at various locations around the world over the year have already been established in the existing literature. Similar to our study, Adnan et al. [35], Salam and Islam [37], and Salam et al. [38] in Bangladesh, and Dabral et al. [39] in India, applied similar recommended values for their studies. In total, 80% of the dataset was utilized for training, while the remaining 20% was used for testing.

2.3. FAO56 Penman–Monteith Model (FAO56-PM Model)

This research used the FAO56-PM equation to estimate daily ETo. Allen et al. [36] proposed this conventional and widely utilized model. The following Equation (1) states the original statistical form of the FAO56-PM model:
ETo = 0.408 Δ R n   G + γ 900 T + 273 U 2 e s   e a Δ + γ 1 + 0.34 U 2
where, ETo represents the reference evapotranspiration (mmd−1), R n is the net radiation at the crop surface (MJm−2 d−1). Allen et al. [38] recommended G (soil heat flux density) =0. FAO’s 56 study describes in full the processes for ETo estimate.
Rn is calculated by Equations (2)–(11):
Rn = Rns − Rnl
Rns = (1 − α)Rs
Rs = [ as + bs n N ]   Ra
Ra = 24   60 π   Gsc   dr   [ ω ssin ( ϕ ) sin ( δ ) + cos ( ϕ ) cos ( δ ) sin ( ω s ) ]
dr = 1   0.033 cos 2 π 365 J
δ = 0.409   sin 2 π 365 J 1.39
ωs= arccos [−tan (ϕ) tan (δ)]
Radians = π/180(decimaldegrees)
Rnl = σ [ T max k 4 +   T min K 4 2 ] ( 0.34 0.14 e a ) [ 1.35 R s R so 0.35 ]
Rso = ( 0.75 + 2 × 10 5 Z )   Ra
U2 is calculated from the following Equation (12), recommended by Allen et al. [36],
U 2 = U z 4.87 In 67.8 z 5.42
where, Rns stands for the net solar or shortwave radiation (MJ m−2 d−1), Rnl is the net outgoing longwave radiation (MJ m−2 d−1), Rs is the global solar or shortwave radiation (MJ m−2 d−1), N and n are, respectively, maximum and actual possible sunshine duration, Ra is extraterrestrial radiation (MJ m−2 d−1), Gsc is solar constant (0.0820 MJ m−2 min−1), dr is inverse relative distance Earth–Sun, ωs is the sunset hour angle (rad), ϕ is latitude (rad), δ is the solar declination (rad), J is the number of the day in the year between 1 (1 January) and 365 or 366 (31 December), σ is Stefan–Boltzmann constant (4.903 × 10−9 MJ K−4 m−2 d−1), α is albedo (α = 0.23), Tmaxk and Tmin k are, respectively, maximum and minimum absolute temperatures during 24 h, and Rso is clear sky solar radiation (MJ m−2 d−1). Allen et al. [36] recommended 0.25 for as and 0.50 for bs. Uz is measured wind speed at Zm above the ground surface (ms−2), and z is the respective station elevation above sea level (m).
According to Allen et al. [36], saturation vapor pressure(es), actual vapor pressure (ea), slope vapor pressure curve (∆), and psychrometric constant (γ) are calculated by the following Equations (13)–(19), respectively:
e s = e 0 T m a x + e 0 T m i n 2
e 0 T m a x = 0.6108   e x p 17.27 T m a x T m a x + 237.3
e 0 T m i n = 0.6108   e x p 17.27 T m i n T m i n + 237.3
e a = H r m e a n 100 e 0 T m a x + e 0 T m i n 2
Δ = 4098 0.6108   e x p 17.27   T T + 237.3 T + 237.3 2
γ = C p P ε λ = 0.665   ×   10 3 P
P = 101.3 293 0.0065 Z 293 5.26
where, es represents the mean saturation vapor pressure (kPa), e 0 T m a x   and   e 0 T m i n represent the saturation vapor pressure at the maximum and lowest temperatures, Ea represents the actual vapor pressure function (kPa), and Hr represents the mean relative humidity. Tave, Tmax, and Tmin are the average, maximum, and minimum air temperatures, respectively, in °C, and exp [ ] is 2.7183 (i.e., the base of natural logarithm) raised to the power [3]. P is atmospheric pressure (kPa), λ is the latent heat of vaporization (2.45 MJ kg−1), Cp is the specific heat at constant pressure (1.013×10−3 MJ kg−1 °C−1), ε is s the ratio of water vapor molecular weight to dry molecular air weight (0.622).

3. Methods

3.1. Machine Learning Model (SVR)

This study employs support vector regression (SVR) as a prevailing and well-known machine learning model to predict ETo values. It is worth noting that the SVR has already been successfully applied and reported in several hydrological modeling studies [40,41,42].
Consider a training set such as T = x i , y i , i = 1 , 2 , , l ; the SVR model aims to map the initial data in a higher dimensional feature space. Therefore, the SVR needs to construct a decision function, f ( x ) , to simulate the nonlinear relationship between the input vector, x i , and the target value, yi. Assuming M as the order of the polynomial of the hyper-surface, the estimated regression function can be written as below:
f ( x j ) = w , x + b = j = 1 M w j x j + b
where . , . means the dot product of two parameters and w j is the coefficient vector. In Figure 2, the geometrical illustration of the linear form of the SVR method is shown. SVR tries to find the widest tube, ε SVR , between the support vectors around the hyper-surface, which is the range between f ( x )   +   ε and f ( x ) ε ( ε is the permitted error threshold, see Figure 2). Therefore, the objective function, 1 2 w 2 , subject to y j f x j ε should be minimized (here y is the target value). For the data outside of the tube, the boundary should be optimized. Thus, two slack factors, ξ > 0 and ξ * > 0 , are introduced and applied for the optimization problem [43].
Min Object   function = 1 2 w 2 + C i = 1 M ξ + ξ * s . t .   y j - f x j ε + ξ * f x j y j ε + ξ ,     j = 1 , 2 , , M
where the constant C > 0 is the penalty coefficient. It should be noted that in real cases, most data cannot be appropriately separated by a linear hyperplane (hyper-surface).
Thus, kernel functions are used to cope with this problem to map the available data from low-dimensional feature space to a higher response space where a linear separation works properly. In this sense, the radial basis function, RBF, has been used in this study.
k x , x i = e x p γ x x i 2
where γ denotes the kernel’s hyperparameter. In a standard SVR model, tuning parameters (e.g., the values for C ,   ε ,   and   γ ) are optimized according to mathematical methodologies such as the least square optimization method. However, based on the conclusion of several researchers [44,45], embedding heuristic algorithms for optimization procedures can improve the accuracy of the standard models. Thus, we have applied integrative SVR models based on individual and hybrid heuristic algorithms. Detailed information regarding the heuristic algorithms and the integration procedure is explained in the following sections.

3.2. Heuristic Optimization Methods (PSO, GWO, and GSA)

In the current research, three types of heuristic algorithms were used either as an individual (e.g., the particle swarm optimization, PSO) or as hybrid algorithms (the grey wolf optimizer, GWO and the gravitational search algorithm, GSA) to optimize the standard SVR model.
The PSO is a well-known population-based optimization algorithm inspired by the social behavior of animals’ behavior for searching food or immigration, such as fish schooling or bird flocking [46]. The PSO method generates the initial population randomly through the search space. During the training process, the best location of each agent (particle) is kept in the algorithm’s memory. Hence, in each iteration, particles in the swarm (the group/herd of the animals) would update their positions based on the following equations [46]:
x n + 1 i = x n i + v n + 1 i
v n + 1 i = ω v n i + c 1 r 1 p n i x n i + c 2 r 2 p n g x n i
where x and v are the position and velocity vectors, respectively. In addition, i stands for the particle, n is the iteration number in the epoch. ω denotes the inertia weight parameter, and r1 and r2 represent the two random numbers between zero and unity (rand (0,1)). Pi is the best position achieved by the ith particle. At the same time, Pg refers to the best position information in the swarm. Similar to the other iteration algorithms, the computation would continue until reaching the final epoch or converging to a predetermined stopping criterion.
Similar to the PSO, the GWO algorithm is another nature-based method inspired by the leadership hierarchy of grey wolf packs as the apex predators. In the GWO, the agents are rated in four types according to the four types of wolves in the leadership hierarchy: alpha, beta, delta, and omega. However, the alpha agents introduce the solution within the searching space. Beta and delta agents represent the second and third-best solutions, and omega agents are the solution candidates. In other words, the training process for seeking the best solution (hunting for the wolves) is carried out by the alpha, beta, and delta wolves. The GWO algorithm generally follows three main steps for searching for the best solution for the response parameter. These steps include [47]:
(i)
Tracking, chasing, and approaching the prey;
(ii)
Pursuing, encircling, and harassing the prey;
(iii)
And finally, getting close to the prey and attacking.
Unlike the previously mentioned PSO and GWO algorithms, the GSA has not considered a nature-based optimization algorithm. However, it is a swarm optimization heuristic algorithm based on Newton’s gravitational law between objects. In this sense, each object with a specific mass exerts a force on the other objects around it. Naturally, the other surrounding objects impose the same force on the object of interest in a mutual way. This process is known as the gravitation force among objects. The gravitational force between two objects complies with two main characteristics: (i) the mass of the objects and (ii) the distance between them [48]. Assuming two objects with a distance equal to R having the masses of M1 and M2, the induced gravitational force (F) between them can be calculated as the following equation:
F = G M 1 M 2 R 2
In the above equation, G denotes the gravitational constant.
The GSA aims to find the best solution according to the movement of particles (objects). Having the initial population. The position and velocity of the particles are updated in compliance with the resultant force (F, see Equation (6)) and acceleration (a) associated with each particle along with the G factor. Considering a particle, its velocity and position vectors would be updated in the next iteration as follows:
x n + 1 i = x n i + v n + 1 i
v n + 1 i = r i v n i + a n i
During the solution process, the particles are attracted and move towards those with greater inertia mass (M). Detailed information regarding the GSA methodology can be found by Rashedi et al. [49].

3.3. Hybrid Optimization Methods (PSOGWO and PSOGSA)

In this study, we have used two hybrid heuristic algorithms, including the PSOGWO (the combination of the particle search swarm optimization, PSO, and the grey wolf optimizer, GWO) and the PSOGSA (the combination of the PSO, and Gravitational Search Algorithm, GSA). In addition to the abovementioned hybrid algorithms, the integrative SVR-PSO model is also applied to evaluate the performances of SVR-PSOGWO and SVR-PSOGSA. In the following, explanations regarding the hybrid optimization methods are given.
-
PSOGWO
Here, the developed hybrid PSOGWO algorithm works on the concepts of the general operation of the original PSO and GWO algorithms (Algorithm 1). It should be noted that the PSO is a robust nature-based algorithm, and it can be successfully utilized for several simulating problems. However, it has been declared that one can reduce the possibility of trapping in local minima by attaching it to the GWO algorithm [50]. In the original PSO algorithm, some particles are allocated random positions that enhance the risk of falling in a local minimum. Attaching the GWO algorithm to the PSO would direct the randomly positioned particles to the improved positions specified by the GWO algorithm. This procedure improves the precision of the PSO algorithm; nonetheless, it encounters the shortcoming of longer running time (computational time). In the following, the pseudo-code for the PSOGWO algorithm is explained. Detailed information regarding the hybridization procedure of the PSOGWO is available at Şenel et al. [50].
Algorithm 1. PSOGWO
•   Setting up parameters
  Epoch: the number of iterations (either set by the user or reached according to the other types of stopping criteria)
  SP: Initial swarm population number (particles in the PSO algorithm)
prob: possibility rate (set by the user)
•   Hybrid procedure
  Initializing particles in the solution space
   FOR i = 1 to Epoch
      FOR j = 1 to SP
      Run PSO (updating the x and v vectors)
      Evaluating the fitness values
      Updating Pg (memorizing the best values of the swarm)
      IF rand (0,1) < prob then (to avoid trapping in local minima) THEN
      Run GWO
          Evaluating the fitness of all wolves
          Updating the positions of the Alpha, Beta, and delta wolves
          Calculating the mean of the position of three best (α, β, δ) wolves
          Returning updated values for the particles in the PSO algorithm
      END IF
   END FOR
END FOR
-
PSOGSA
In developing the hybrid version of PSOGSA (Algorithm 1), we tried to cope with the main shortcomings of the individual heuristic algorithms. According to the concepts of forming the GSA, particles move towards the one with greater mass. This procedure makes the algorithm a good candidate for conducting an efficient exploitation phase. Nevertheless, this upside characteristic of the GSA might cause weakness in properly searching the whole response domain, i.e., the deficiency in the global exploring potential. On the other hand, as mentioned earlier, the PSO algorithm has a strong global exploring ability (i.e., exploration phase).
Nevertheless, according to some reports [51], its exploiting phase is not as sufficient as it should be. Therefore, similar to the hybrid PSOGWO algorithm, the hybridization process has been carried out to overcome this drawback using the GSA (instead of the GWO) so that one can evaluate the precision and potential of these two hybrid methods [52]. In the following, the pseudo-code developed for the PSOGSA algorithm is explained. Further information for the applied hybridization procedure of the PSOGSA is available at Song et al. [53].
Algorithm 2. PSOGSA
•   Setting up initial values and parameters
    Epoch: the number of iterations; SP: Initial swarm population number; prob: possi bility rate
•   Hybrid procedure
Initializing particles in the solution space
   FOR i = 1 to Epoch
      FOR j = 1 to SP
      Run PSO
      Evaluating the fitness values of the particles updating Pg
      IF rand(0,1) < prob then (to avoid trapping in local minima) THEN
      Run GSA
          Computing the resultant force (F) and the acceleration (a)
          Updating values for the velocity and positions (Pi)
          Returning updated values for the particles in the PSO algorithm
      END IF
   END FOR
END FOR

3.4. Performance Evaluation

Two novel SVR methods combined with hybrid PSO-GWO and PSO-GSA meta-heuristic algorithms are compared with the single SVR and SVR-PSO method in estimating monthly ET0 using climatic data involving Tmin, Tmax, Ra, Rs, U2, and HR. The following statistics were utilized for assessing the implemented methods:
      R M S E : R o o t   M e a n   S q u a r e   E r r o r = 1 N i = 1 N [ ( E T 0 ) i ( E T e ) i ] 2    
M A E : M e a n   A b s o l u t e   E r r o r = 1 N i = 1 N | ( E T 0 ) i ( E T e ) i |
N S E : N a s h S u t c l i f f e   E f f i c i e n c y = 1 i = 1 N [ ( E T 0 ) i ( E T e ) i ] 2 i = 1 N [ ( E T 0 ) i E T ¯ 0 ] 2 , < N S E 1        
where E T 0 ,   E T e ,   E T ¯ 0 are FAO56-PM ET0, estimated an average FAO56-PM ET0, respectively, and N indicates the data quantity. Distinct values were attempted to reach the optimal values reported in Table 2. This table also represents each algorithm’s population number, iterations, and several runs. Table 3 lists the input combinations considered in this study and the corresponding variables.

4. Results and Discussion

Four AI models with nine different input combinations (Table 3) are applied to three selected climatic stations. Table 4 sums up the training and testing results of the implemented methods in estimating the ETo of Bogra Station. In all methods, full data offer the highest accuracy, and the SVR-PSOGWO acts better than the other models; improvement in RMSE is 27%, 21%, and 5% in the testing stage compared to SVR, SVR-PSO, and SVR-PSOGSA, respectively. Other models are ranked from best to worst as SVR-PSOGSA > SVR-PSO > SVR in the estimation of monthly ETo.
Training efficiencies of the SVR-PSOGWO and SVR-GSA are almost equal, and they could approximate the phenomenon better compared to SVR and SVR-PSO; viz. NSE increases from 0.988 to 0.995, and RMSE decreases from 0.073/0.071 to 0.042/0.043. Another important information derived from the results is that the second input case, including Tmin, Tmax, and Ra, also offers good accuracy in estimating monthly ETo for this station. This might carry importance in practical applications because measuring temperatures is very easy, and Ra can be directly calculated from the Julian date.
Table 5 lists the accuracy of the single and hybrid SVR-based models for various input cases of Rajshahi Station. In contrast to the previous station, the methods produced the best accuracy in different input cases. The SVR-PSOGSA with Tmin, Tmax, Ra, Rs, and U2 acted better than the other models; improvements in RMSE by 32%, 20%, and 3% in the testing stage compared to the SVR, SVR-PSO, and SVR-PSOGWO, respectively. The SVR-PSOGWO, SVR-PSO, and SVR follow the best accuracy of SVR-PSOGSA in estimating monthly ETo. In the training (simulation) stage, however, the SVR-PSOGWO is better fitted to data than SVR-PSOGSA. In this station, the second input case (Tmin, Tmax, Ra) may also be another alternative in estimating monthly ETo when the other climatic data (Rs, U2, and HR) are missing. The other important thing is a slight difference between the input cases V and IX for the SVR-PSOGWO. This input case (Tmin, Tmax, Ra, Rs) also offers good accuracy in all methods. In this station, Rs seems more effective on ETo than Bogra Station, which can be seen from the differences between II and V input cases.
The SVR-based methods’ training and testing results are summarized in Table 6 for Rangpur Station. The SVR, SVR-PSO, and SVR-PSOGWO offered the best accuracy for the IX input case, while the SVR-PSOGSA had the best outcomes for V inputs. However, there is a slight difference between inputs IX and V for this method. The SVR-PSOGWO with the Tmin, Tmax, Ra, Rs, U2, and HR acted better than the other models; improvements in RMSE by 23%, 14%, and 19% in the testing stage compared to SVR, SVR-PSO, and SVR-PSOGSA, respectively. The other models are ranked as SVR-PSOGSA > SVR-PSO > SVR in the estimation of monthly ETo. However, variations in the predictive capabilities between two novel MA models using the combinations of three input climatic variables plans (Tmax, Tmin, and Ra) were slight (RMSE < 10%). As evident from the first part of (columns 3–6) Table 6, the training (simulation) accuracy of the SVR-PSOGWO is also better than the other models; NSE increases from 0.976/0.988/0.990 to 0.998, and RMSE decreases from 0.139/0.099/0.089 to 0.041 by applying the SVR-PSOGWO compared to SVR/SVR-PSO/SVR-PSOGSA.
Table 4, Table 5 and Table 6 clearly show that the SVR-PSOGWO generally offered better accuracy than the other methods in estimating monthly ETo. In the second station, the SVR-PSOGSA performed superior to the SVR-PSOGWO, but the difference was marginal. In addition, the SVR-PSOGSA produced inferior results compared to SVR-PSO in the third station. It is also clear that the single SVR offered the worst outcomes in all stations. This matter indicates the necessity of hybrid metaheuristic algorithms in the training of the SVR method in ETo prediction. The input combination of Tmin, Tmax, and Rs (V input case) also offered good accuracy in ETo estimation.
Figure 3, Figure 4 and Figure 5 reveal the time variation diagrams of the FAO56-PM ETo and estimated ETo by the optimal SVR-based models for the Bogra, Rajshahi, and Rangpur stations, respectively. It is clear from the figures that the estimates of the SVR-based hybrid models are closer to the FAO56-PM ETo values than the single SVR model.
The scatter diagrams of the estimated ETo are illustrated in Figures S1–S3 (see Supplementary Materials) for the three stations. All graphs clearly show that the SVR-PSOGWO has less scattered estimates with the highest R2 values (0.9334, 0.9508, 0.9501) for the Bogra, Rajshahi, and Rangpur stations, respectively. At the same time, the single SVR provides the most scattered estimations.
Similar to our study, Granata [54] found that the SVR was the weakest model in the humid tropical region of Florida, USA. Shiri [55] stated that hybrid tree-based methods were effectively used for monthly ETo forecasting in southern Iran. Similarly, Huang et al. [56] argued that SVM models showed the worst performance with the lowest increases in testing RMSE from 4.1% to 37.3%. Moreover, the metaheuristic algorithms, e.g., SVR-PSOGWO and SVR-PSOGSA models, generally perform better for classification issues but lower for regression issues because they cannot give static outputs. In such a case, these two novel MA tools cannot produce better prediction precision, which leads to an overfitting issue when noisy testing data are employed for estimating ETo. Overall, the SVR-PSOGWO and SVR-PSOGSA offer better accuracy than the SVR-PSO method in monthly ETo estimation. The main advantages of these hybrid algorithms (GWO and GSA) are incorporated to improve the exploitation ability of PSO, as PSO is known for better exploration but lags in exploitation. However, the two novels’ MA models showed stability with a satisfactory % increase with the lowest increases in the testing RMSE stage.
Salam and Islam [37] found similar results in Bangladesh’s subhumid tropical region. On the other hand, the SVR and SVR-PSO models have many hyper-variables, which need to be carefully tuned for monthly ETo purposes. Our findings suggest that the two newly novel metaheuristic algorithm, with much development upon the traditional SVR model, are a promising tool for preventing the difficulty of the overfitting issue. Thus, K-fold cross-validation is needed to select the optimal parameter and validate the better stability of the soft computing methods. These deserve further investigation.
Khosravi et al. [57] investigated the accuracy of four ANFIS including single ANFIS, hybrid ANFIS models tuned with differential evolution (ANFIS-DE), genetic algorithm (ANFIS-GA), and imperialistic competitive algorithm (ANFIS-ICA), and five data mining models including M5P, random forest (RF), random tree (RT), reduced error runing tree and Kstar models in predicting monthly ETo of Baghdad and Mosul (Iraq) using different meteorological input combinations involving sunshine hours (n), maximum and minimum temperature (Tmax and Tmin), wind speed (U2), rainfall (P), and relative humidity (HR). For the Bagdat Station, all the above-mentioned models except RT provided the best accuracy for the input combination of n, Tmax, and RH, while for the Mosul station, the input combinations of (Tmax, Tmin, n, HR, and U2) and (Tmax, Tmin, n, HR, U2, P) produced best ETo predictions. Khosravi et al. [57,58] reported that the best model accuracy varies concerning the type of the machine learning or data mining methods. This can be explained by the fact that each method has a different structure and different calculation processes.
Reasonably, many soft computing methods generally consider a sole dataset for validating for testing, which may mislead or give partly strong conclusions because the performance capabilities are solely described to an exact distribution of the main dataset and may not be identical [16,30]. The definition of the data combination included in each data combination is important for evaluating the method’s performance. Therefore, the performance of any machine learning method is favorably based on the data distributions and climatic areas. For example, an analogous model showed various performances in various cited works under different climate conditions [59]. Similar to this study, some scholars, e.g., Adnan et al. [60], reported that the hybrid model outperformed other state-of-the-art models. Thus, our work suggests that the SVR-based metaheuristic algorithm can be a likely soft computing method for good accuracy and consistency for monthly ETo prediction using limited datasets in the data-scarce humid region of north-western Bangladesh and is applicable for like climatic settings around the world. However, further study is needed to assess the performance of potential SVR-based metaheuristic models at different time durations, including hourly or daily basis or in other analogous areas, worldwide with humid tropical climate conditions

5. Conclusions

The present study investigated the ability of two novel SVR methods merged with hybrid PSOGSA and PSOGWO meta-heuristic algorithms in estimating ETo using climatic data as inputs. Having a better evaluation, the SVR-PSOGSA and SVR-PSOGWO methods were also compared with single SVR and SVR-PSO methods. The outcomes provided the following conclusions:
(i)
Monthly discharge, Tmin, Tmax, Ra, Rs, U2, and HR data from three stations were used for assessing the above-mentioned methods. Based on the root mean square error, mean absolute error, Nash–Sutcliffe efficiency and determination coefficient and graphical methods, the SVR–PSOGWO was superior to the other methods, followed by the SVR–PSOGSA, SVR–PSO, and SVR. This implies the necessity of hybrid metaheuristic algorithms in SVR training.
(ii)
It was observed that the input combination involving whole climatic data generally produced the best accuracy. The SVR–PSOGWO with Tmin, Tmax, Ra, Rs, and U2 inputs improved the accuracy of single SVR by 27%, 32%, and 23% for Bogra, Rajshahi, and Rangpur stations with respect to root mean square errors in the testing stage, respectively. The second input combination comprising Tmin, Tmax, and Ra also provided good accuracy (NSE ranges from 0.808 to 0.897). The models with this input combination can be a good alternative when other climatic data are unavailable. The viability of the presented hybrid metaheuristic algorithms can be assessed for improving other machine learning methods such as extreme leaning machine, neural networks, or neuro-fuzzy systems in future studies.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/agronomy13010098/s1. Figure S1. Scatterplots of the observed and predicted ET different SVR based models in the test period of Bogra Station. Figure S2. Scatterplots of the observed and predicted ET different SVR based models in the test period of Rajshahi Station. Figure S3. Scatterplots of the observed and predicted ET different SVR based models in the test period of Rangpur Station.

Author Contributions

Conceptualization: R.M.A.I., R.R.M. and O.K.; formal analysis: R.R.M.; validation: R.M.A.I., O.K., Z.C. and M.Z.-K.; supervision: Z.C. and O.K.; writing—original draft: R.M.A.I., Z.C., A.R.M.T.I., O.K., R.R.M., A.K. and M.Z.-K.; visualization: R.M.A.I. and A.R.M.T.I.; investigation: A.R.M.T.I., O.K., A.K. and Z.C., writing—review and editing: R.M.A.I., Z.C., O.K., R.R.M., A.K. and M.Z.-K. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study will be available upon interesting request from the corresponding author.

Conflicts of Interest

There is no conflict of interest in this study.

References

  1. Bozkurt Çolak, Y.; Yazar, A.; Alghory, A.; Tekin, S. Evaluation of crop water stress index and leaf water potential for differentially irrigated quinoa with surface and subsurface drip systems. Irrig. Sci. 2021, 39, 81–100. [Google Scholar] [CrossRef]
  2. Huang, D.; Wang, J.; Khayatnezhad, M. Estimation of Actual Evapotranspiration Using Soil Moisture Balance and Remote Sensing. Iran. J. Sci. Technol. Trans. Civ. Eng. 2021, 45, 2779–2786. [Google Scholar] [CrossRef]
  3. Ghiat, I.; Mackey, H.R.; Al-Ansari, T. A Review of Evapotranspiration Measurement Models, Techniques and Methods for Open and Closed Agricultural Field Applications. Water 2021, 13, 2523. [Google Scholar] [CrossRef]
  4. Lugato, E.; Lavallee, J.M.; Haddix, M.L.; Panagos, P.; Cotrufo, M.F. Different climate sensitivity of particulate and mineral-associated soil organic matter. Nat. Geosci. 2021, 14, 295–300. [Google Scholar] [CrossRef]
  5. Shah, S.; Duan, Z.; Song, X.; Li, R.; Mao, H.; Liu, J.; Ma, T.; Wang, M. Evaluating the added value of multi-variable calibration of SWAT with remotely sensed evapotranspiration data for improving hydrological modeling. J. Hydrol. 2021, 603, 127046. [Google Scholar] [CrossRef]
  6. Fan, G.; Sarabandi, A.; Yaghoobzadeh, M. Evaluating the climate change effects on temperature, precipitation and evapotranspiration in eastern Iran using CMPI5. Water Supply 2021, 21, 4316–4327. [Google Scholar] [CrossRef]
  7. Machakaire, A.T.B.; Steyn, J.M.; Franke, A.C. Assessing evapotranspiration and crop coefficients of potato in a semi-arid climate using Eddy Covariance techniques. Agric. Water Manag. 2021, 255, 107029. [Google Scholar] [CrossRef]
  8. Makwana, J.J.; Tiwari, M.K.; Deora, B.S. Development and comparison of artificial intelligence models for estimating daily reference evapotranspiration from limited input variables. Smart Agric. Technol. 2023, 3, 100115. [Google Scholar] [CrossRef]
  9. Rodrigues, G.C.; Braga, R.P. Estimation of Reference Evapotranspiration during the Irrigation Season Using Nine Temperature-Based Methods in a Hot-Summer Mediterranean Climate. Agriculture 2021, 11, 124. [Google Scholar] [CrossRef]
  10. Samadi, S.Z. Assessing the sensitivity of SWAT physical parameters to potential evapotranspiration estimation methods over a coastal plain watershed in the southeastern United States. Hydrol. Res. 2016, 48, 395–415. [Google Scholar] [CrossRef]
  11. Gisolo, D.; Previati, M.; Bevilacqua, I.; Canone, D.; Boetti, M.; Dematteis, N.; Balocco, J.; Ferrari, S.; Gentile, A.; N’Sassila, M.; et al. A calibration free radiation driven model for estimating actual evapotranspiration of mountain grasslands (CLIME-MG). J. Hydrol. 2022, 610, 127948. [Google Scholar] [CrossRef]
  12. Abeysiriwardana, H.D.; Muttil, N.; Rathnayake, U. A Comparative Study of Potential Evapotranspiration Estimation by Three Methods with FAO Penman–Monteith Method across Sri Lanka. Hydrology 2022, 9, 206. [Google Scholar] [CrossRef]
  13. Wanniarachchi, S.; Sarukkalige, R. A Review on Evapotranspiration Estimation in Agricultural Water Management: Past, Present, and Future. Hydrology 2022, 9, 123. [Google Scholar] [CrossRef]
  14. Tejada, A.T.; Ella, V.B.; Lampayan, R.M.; Reaño, C.E. Modeling Reference Crop Evapotranspiration Using Support Vector Machine (SVM) and Extreme Learning Machine (ELM) in Region IV-A, Philippines. Water 2022, 14, 754. [Google Scholar] [CrossRef]
  15. Abedi-Koupai, J.; Dorafshan, M.-M.; Javadi, A.; Ostad-Ali-Askari, K. Estimating potential reference evapotranspiration using time series models (case study: Synoptic station of Tabriz in northwestern Iran). Appl. Water Sci. 2022, 12, 212. [Google Scholar] [CrossRef]
  16. Chia, M.Y.; Huang, Y.F.; Koo, C.H. Improving reference evapotranspiration estimation using novel inter-model ensemble approaches. Comput. Electron. Agric. 2021, 187, 106227. [Google Scholar] [CrossRef]
  17. Dou, X.; Yang, Y. Evapotranspiration estimation using four different machine learning approaches in different terrestrial ecosystems. Comput. Electron. Agric. 2018, 148, 95–106. [Google Scholar] [CrossRef]
  18. Nourani, V.; Elkiran, G.; Abdullahi, J. Multi-station artificial intelligence based ensemble modeling of reference evapotranspiration using pan evaporation measurements. J. Hydrol. 2019, 577, 123958. [Google Scholar] [CrossRef]
  19. Antonopoulos, V.Z.; Antonopoulos, A.V. Daily reference evapotranspiration estimates by artificial neural networks technique and empirical equations using limited input climate variables. Comput. Electron. Agric. 2017, 132, 86–96. [Google Scholar] [CrossRef]
  20. Zhang, G.; Liu, B.; Zhu, T.; Zhou, A.; Zhou, W. Visual privacy attacks and defenses in deep learning: A survey. Artif. Intell. Rev. 2022, 55, 4347–4401. [Google Scholar] [CrossRef]
  21. Emadi, A.; Zamanzad-Ghavidel, S.; Fazeli, S.; Zarei, S.; Rashid-Niaghi, A. Multivariate modeling of pan evaporation in monthly temporal resolution using a hybrid evolutionary data-driven method (case study: Urmia Lake and Gavkhouni basins). Environ. Monit. Assess. 2021, 193, 355. [Google Scholar] [CrossRef] [PubMed]
  22. Mehdizadeh, S. Estimation of daily reference evapotranspiration (ETo) using artificial intelligence methods: Offering a new approach for lagged ETo data-based modeling. J. Hydrol. 2018, 559, 794–812. [Google Scholar] [CrossRef]
  23. Tikhamarine, Y.; Malik, A.; Souag-Gamane, D.; Kisi, O. Artificial intelligence models versus empirical equations for modeling monthly reference evapotranspiration. Environ. Sci. Pollut. Res. 2020, 27, 30001–30019. [Google Scholar] [CrossRef] [PubMed]
  24. Tikhamarine, Y.; Malik, A.; Kumar, A.; Souag-Gamane, D.; Kisi, O. Estimation of monthly reference evapotranspiration using novel hybrid machine learning approaches. Hydrol. Sci. J. 2019, 64, 1824–1842. [Google Scholar] [CrossRef]
  25. Seifi, A.; Riahi, H. Estimating daily reference evapotranspiration using hybrid gamma test-least square support vector machine, gamma test-ANN, and gamma test-ANFIS models in an arid area of Iran. J. Water Clim. Chang. 2018, 11, 217–240. [Google Scholar] [CrossRef]
  26. Hossain, M.A.; Rahman, M.M.; Hasan, S.S. Application of combined drought index to assess meteorological drought in the south western region of Bangladesh. Phys. Chem. Earth Parts A/B/C 2020, 120, 102946. [Google Scholar] [CrossRef]
  27. Mojid, M.A.; Aktar, S.; Mainuddin, M. Rainfall-induced recharge-dynamics of heavily exploited aquifers—A case study in the North-West region of Bangladesh. Groundw. Sustain. Dev. 2021, 15, 100665. [Google Scholar] [CrossRef]
  28. Sumon, K.A.; Rashid, H.; Peeters, E.T.H.M.; Bosma, R.H.; Van den Brink, P.J. Environmental monitoring and risk assessment of organophosphate pesticides in aquatic ecosystems of north-west Bangladesh. Chemosphere 2018, 206, 92–100. [Google Scholar] [CrossRef]
  29. Jamal, M.R.; Kristiansen, P.; Kabir, M.J.; Kumar, L.; Lobry de Bruyn, L. Trajectories of cropping system intensification under changing environment in south-west coastal Bangladesh. Int. J. Agric. Sustain. 2022, 20, 722–742. [Google Scholar] [CrossRef]
  30. Dewan, A.; Shahid, S.; Bhuian, M.H.; Hossain, S.M.J.; Nashwan, M.S.; Chung, E.-S.; Hassan, Q.K.; Asaduzzaman, M. Developing a high-resolution gridded rainfall product for Bangladesh during 1901–2018. Sci. Data 2022, 9, 471. [Google Scholar] [CrossRef]
  31. Salehie, O.; Ismail, T.; Shahid, S.; Ahmed, K.; Adarsh, S.; Asaduzzaman, M.; Dewan, A. Ranking of gridded precipitation datasets by merging compromise programming and global performance index: A case study of the Amu Darya basin. Theor. Appl. Climatol. 2021, 144, 985–999. [Google Scholar] [CrossRef]
  32. Pour, S.H.; Shahid, S.; Chung, E.-S.; Wang, X.-J. Model output statistics downscaling using support vector machine for the projection of spatial and temporal changes in rainfall of Bangladesh. Atmos. Res. 2018, 213, 149–162. [Google Scholar] [CrossRef]
  33. Wahiduzzaman, M.; Luo, J.-J. A statistical analysis on the contribution of El Niño–Southern Oscillation to the rainfall and temperature over Bangladesh. Meteorol. Atmos. Phys. 2021, 133, 55–68. [Google Scholar] [CrossRef]
  34. Rouf, M.A.; Uddin, M.K.; Debsarma, S.K.; Rahman, M.M. Climate of Bangladesh: An Analysis of Northwestern and Southwestern Part Using High Resolution Atmosphere-Ocean General Circulation Model (AOGCM). Agriculturists 2012, 9, 143–154. [Google Scholar] [CrossRef] [Green Version]
  35. Adnan, R.M.; Mostafa, R.R.; Islam, A.R.M.T.; Kisi, O.; Kuriqi, A.; Heddam, S. Estimating reference evapotranspiration using hybrid adaptive fuzzy inferencing coupled with heuristic algorithms. Comput. Electron. Agric. 2021, 191, 106541. [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, p. D05109. [Google Scholar]
  37. Salam, R.; Islam, A.R.M.T. Potential of RT, Bagging and RS ensemble learning algorithms for reference evapotranspiration prediction using climatic data-limited humid region in Bangladesh. J. Hydrol. 2020, 590, 125241. [Google Scholar] [CrossRef]
  38. Salam, R.; Islam, A.R.M.T.; Pham, Q.B.; Dehghani, M.; Al Ansari, N.; Linh, N.T.T. The optimal alternative for quantifying reference evapotranspiration in climatic sub-regions of Bangladesh. Sci. Rep. 2020, 10, 20171. [Google Scholar] [CrossRef]
  39. Dabral, P.P.; Mor, N.; Jhajharia, D. Time series modelling of monthly reference evapotranspiration for Bikaner, Rajasthan (India). Indian J. Soil Conserv. 2018, 46, 42–51. [Google Scholar]
  40. Aliku, O.; Oshunsanya, S.O.; Aiyelari, E.A. Estimation of crop evapotranspiration of Okra using drainage Lysimeters under dry season conditions. Sci. Afr. 2022, 16, e01189. [Google Scholar] [CrossRef]
  41. Moazenzadeh, R.; Mohammadi, B.; Shamshirband, S.; Chau, K.-W. Coupling a firefly algorithm with support vector regression to predict evaporation in northern Iran. Eng. Appl. Comput. Fluid Mech. 2018, 12, 584–597. [Google Scholar] [CrossRef] [Green Version]
  42. Adnan, R.M.; Kisi, O.; Mostafa, R.R.; Ahmed, A.N.; El-Shafie, A. The potential of a novel support vector machine trained with modified mayfly optimization algorithm for streamflow prediction. Hydrol. Sci. J. 2022, 67, 161–174. [Google Scholar] [CrossRef]
  43. Sreedhara, B.M.; Rao, M.; Mandal, S. Application of an evolutionary technique (PSO–SVM) and ANFIS in clear-water scour depth prediction around bridge piers. Neural Comput. Appl. 2019, 31, 7335–7349. [Google Scholar] [CrossRef]
  44. Ikram, R.M.A.; Dai, H.-L.; Ewees, A.A.; Shiri, J.; Kisi, O.; Zounemat-Kermani, M. Application of improved version of multi verse optimizer algorithm for modeling solar radiation. Energy Rep. 2022, 8, 12063–12080. [Google Scholar] [CrossRef]
  45. Meng, E.; Huang, S.; Huang, Q.; Fang, W.; Wang, H.; Leng, G.; Wang, L.; Liang, H. A Hybrid VMD-SVM Model for Practical Streamflow Prediction Using an Innovative Input Selection Framework. Water Resour. Manag. 2021, 35, 1321–1337. [Google Scholar] [CrossRef]
  46. Huang, W.; Liu, H.; Zhang, Y.; Mi, R.; Tong, C.; Xiao, W.; Shuai, B. Railway dangerous goods transportation system risk identification: Comparisons among SVM, PSO-SVM, GA-SVM and GS-SVM. Appl. Soft Comput. 2021, 109, 107541. [Google Scholar] [CrossRef]
  47. Mirjalili, S.; Mirjalili, S.M.; Lewis, A. Grey Wolf Optimizer. Adv. Eng. Softw. 2014, 69, 46–61. [Google Scholar] [CrossRef] [Green Version]
  48. Mittal, H.; Tripathi, A.; Pandey, A.C.; Pal, R. Gravitational search algorithm: A comprehensive analysis of recent variants. Multimed. Tools Appl. 2021, 80, 7581–7608. [Google Scholar] [CrossRef]
  49. Rashedi, E.; Nezamabadi-pour, H.; Saryazdi, S. GSA: A Gravitational Search Algorithm. Inf. Sci. 2009, 179, 2232–2248. [Google Scholar] [CrossRef]
  50. Şenel, F.A.; Gökçe, F.; Yüksel, A.S.; Yiğit, T. A novel hybrid PSO–GWO algorithm for optimization problems. Eng. Comput. 2019, 35, 1359–1373. [Google Scholar] [CrossRef]
  51. Song, B.; Wang, Z.; Zou, L. An improved PSO algorithm for smooth path planning of mobile robots using continuous high-degree Bezier curve. Appl. Soft Comput. 2021, 100, 106960. [Google Scholar] [CrossRef]
  52. Eappen, G.; Shankar, T. Hybrid PSO-GSA for energy efficient spectrum sensing in cognitive radio network. Phys. Commun. 2020, 40, 101091. [Google Scholar] [CrossRef]
  53. Song, B.; Xiao, Y.; Xu, L. Design of fuzzy PI controller for brushless DC motor based on PSO–GSA algorithm. Syst. Sci. Control Eng. 2020, 8, 67–77. [Google Scholar] [CrossRef]
  54. Granata, F. Evapotranspiration evaluation models based on machine learning algorithms—A comparative study. Agric. Water Manag. 2019, 217, 303–315. [Google Scholar] [CrossRef]
  55. Shiri, J. Improving the performance of the mass transfer-based reference evapotranspiration estimation approaches through a coupled wavelet-random forest methodology. J. Hydrol. 2018, 561, 737–750. [Google Scholar] [CrossRef]
  56. 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]
  57. Khosravi, K.; Daggupati, P.; Alami, M.T.; Awadh, S.M.; Ghareb, M.I.; Panahi, M.; Pham, B.T.; Rezaie, F.; Qi, C.; Yaseen, Z.M. Meteorological data mining and hybrid data-intelligence models for reference evaporation simulation: A case study in Iraq. Comput. Electron. Agric. 2019, 167, 105041. [Google Scholar] [CrossRef]
  58. Khosravi, K.; Mao, L.; Kisi, O.; Yaseen, Z.M.; Shahid, S. Quantifying hourly suspended sediment load using data mining models: Case study of a glacierized Andean catchment in Chile. J. Hydrol. 2018, 567, 165–179. [Google Scholar] [CrossRef]
  59. Zounemat-Kermani, M.; Kisi, O.; Piri, J.; Mahdavi-Meymand, A. Assessment of artificial intelligence–based models and metaheuristic algorithms in modeling evaporation. J. Hydrol. Eng. 2019, 24, 04019033. [Google Scholar] [CrossRef]
  60. 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. 2020, 586, 124371. [Google Scholar] [CrossRef]
Figure 1. Location of selected stations.
Figure 1. Location of selected stations.
Agronomy 13 00098 g001
Figure 2. A schematic diagram of the two-dimensional feature space mapped data using a linear SVR.
Figure 2. A schematic diagram of the two-dimensional feature space mapped data using a linear SVR.
Agronomy 13 00098 g002
Figure 3. Time variation graphs of the FAO56-PM Eto and predicted ETo by different SVR-based models in the test period of Bogra Station.
Figure 3. Time variation graphs of the FAO56-PM Eto and predicted ETo by different SVR-based models in the test period of Bogra Station.
Agronomy 13 00098 g003
Figure 4. Time variation graphs of the FAO56-PM Eto and predicted ETo by different SVR-based models in the test period of Rajshahi Station.
Figure 4. Time variation graphs of the FAO56-PM Eto and predicted ETo by different SVR-based models in the test period of Rajshahi Station.
Agronomy 13 00098 g004
Figure 5. Time variation graphs of the FAO56-PM Eto and predicted ETo by different SVR-based models in the test period of Rangpur Station.
Figure 5. Time variation graphs of the FAO56-PM Eto and predicted ETo by different SVR-based models in the test period of Rangpur Station.
Agronomy 13 00098 g005
Table 1. The geographical sites and daily mean values of three selected in situ observational datasets in northwest Bangladesh.
Table 1. The geographical sites and daily mean values of three selected in situ observational datasets in northwest Bangladesh.
StationsLatitude
(N)
Longitude
(E)
Altitude
(m)
Tmax
(°C)
Tmin
(°C)
Rs
(MJm−2d−1)
Ra
(MJm−2d−1)
U2
(ms−1)
Hr
(%)
ETo
(mmd−1)
Bogura24.8589.3717.9029.9121.0416.6932.841.0678.143.69
Rajshahi24.3788.719.5030.1120.5617.2532.971.0078.183.78
Rangpur25.7389.2732.6128.9620.2516.6032.631.0380.263.53
Table 2. Parameters setting of each optimization algorithm.
Table 2. Parameters setting of each optimization algorithm.
SVR C 10
γ 0.1
ε 0.01
Kernel typeRadial bias function (RBF)
PSOCognitive component ( c 1 )2
Social component ( c 2 )2
Inertia weight0.2–0.9
GWO a decreased from 2 to 0
GSAInitial gravitational constant G 0 100
Search parameter α 20
PSOGWOAs in both PSO and GWO
PSOGSAAs in both PSO and GSA
All algorithmsPopulation25
Number of iterations100
Number of runs for each algorithm8
Table 3. The input combinations used for model development.
Table 3. The input combinations used for model development.
Input CombinationsVariables
(i)Tmin, Tmax
(ii)Tmin, Tmax, Ra
(iii)Tmin, Tmax, Rs
(iv)Tmin, Tmax, U2
(v)Tmin, Tmax, Ra, Rs
(vi)Tmin, Tmax, Rs, U2
(vii)Tmin, Tmax, Ra, U2
(viii)Tmin, Tmax, Ra, Rs, U2
(ix)Tmin, Tmax, Ra, Rs, U2, HR
Table 4. The results of Station 1 (Bogra) using SVR-based models.
Table 4. The results of Station 1 (Bogra) using SVR-based models.
ModelsInput
Combinations
TrainingTesting
RMSEMAENSER2RMSEMAENSER2
SVRI0.5080.4030.6820.6820.5970.5110.5840.692
II0.3960.3100.8070.8070.3900.3240.8230.834
III0.2090.1540.9460.9460.4320.3160.7820.875
IV0.3950.3070.8080.8210.4700.3320.7430.749
V0.1900.1360.9550.9550.3260.1980.8760.898
VI0.1720.1280.9630.9630.4410.3060.7730.866
VII0.3110.2400.8810.8810.3790.2690.8330.863
VIII0.1470.1070.9740.9740.3730.2500.8380.891
IX0.1000.0730.9880.9880.3520.2320.8560.912
SVR-PSOI0.4110.3380.7920.8030.4980.4190.7100.743
II0.3080.2320.8830.8840.3640.3050.8450.871
III0.1820.1330.9590.9600.3980.2900.8150.888
IV0.3170.2400.8770.8770.4300.3120.7840.792
V0.1530.1090.9710.9710.3530.2440.8550.909
VI0.1520.1120.9720.9720.4260.2960.7880.876
VII0.2410.1900.9290.9290.3460.2630.8600.883
VIII0.1270.0960.9800.9800.4200.3180.7940.905
IX0.0970.0710.9880.9890.3350.2220.8690.923
SVR-
PSOGSA
I0.3690.2930.8320.8330.4900.3910.7200.761
II0.2380.1830.9300.9300.3680.2910.8420.916
III0.1310.0980.9790.9790.3090.2130.8880.901
IV0.2810.2210.9030.9030.4200.3050.7940.805
V0.1270.0930.9800.9800.2830.1810.9070.919
VI0.1180.0830.9830.9830.4720.3670.7400.890
VII0.2150.1700.9430.9430.3040.2160.8920.893
VIII0.1090.0760.9850.9850.3690.2470.8410.897
IX0.0610.0430.9950.9950.2920.1780.9000.927
SVR-
PSOGWO
I0.3160.2470.8770.8770.5120.3880.6940.782
II0.2330.1770.9330.9330.2980.2410.8970.931
III0.1480.1090.9730.9730.3300.2190.8730.898
IV0.2830.2220.9020.9020.3900.2780.8230.832
V0.1060.0780.9860.9860.3060.2080.8910.927
VI0.1220.0860.9820.9820.4460.3430.7680.892
VII0.2150.1680.9430.9430.3050.2330.8920.922
VIII0.1130.0820.9840.9840.3010.2120.8940.929
IX0.0610.0420.9950.9950.2770.1670.9110.933
Table 5. The results of Station 2 (Rajshahi) using SVR-based models.
Table 5. The results of Station 2 (Rajshahi) using SVR-based models.
ModelsInput CombinationsTrainingTesting
RMSEMAENSER2RMSEMAENSER2
SVRI0.3450.2660.8920.8920.5500.4140.7170.856
II0.3070.2310.9140.9140.4540.3200.8070.883
III0.2500.1760.9430.9430.3580.2860.8800.908
IV0.2810.2250.9280.9280.3950.3150.8540.873
V0.2130.1440.9590.9590.3190.2260.9090.919
VI0.3240.2580.9050.9050.4540.3630.8070.843
VII0.2930.2360.9220.9220.3930.3180.8550.864
VIII0.2610.2010.9380.9470.3660.2770.8750.903
IX0.3230.2280.9050.9210.3270.2390.9060.913
SVR-PSOI0.3060.2310.9150.9150.5270.4000.7400.872
II0.2600.1910.9390.9390.4310.3360.8260.906
III0.2260.1510.9540.9550.3460.2710.8880.918
IV0.2540.1990.9410.9420.3920.2930.8570.881
V0.2080.1430.9610.9610.2900.2260.9210.936
VI0.1820.1380.9700.9700.3150.2460.9070.920
VII0.2660.2090.9360.9360.3450.2610.8890.911
VIII0.2060.1520.9610.9620.3150.2470.9070.922
IX0.2270.1730.9530.9530.2980.2300.9170.929
SVR-
PSOGSA
I0.2760.2060.9310.9310.5250.4120.7420.875
II0.2450.1820.9460.9460.3790.3050.8650.928
III0.1860.1210.9690.9690.3020.2100.9150.932
IV0.2230.1750.9550.9550.3770.2790.8670.893
V0.1920.1210.9670.9670.2710.1990.9310.945
VI0.1530.1110.9790.9790.3020.2190.9150.928
VII0.2180.1660.9570.9570.3160.2350.9070.920
VIII0.1340.0960.9840.9840.2410.1440.9460.947
IX0.0920.0520.9920.9920.2520.1470.9410.944
SVR-
PSOGWO
I0.2640.1970.9370.9370.4960.3850.7700.884
II0.2430.1770.9460.9460.3890.3160.8590.939
III0.1940.1270.9660.9660.3120.2330.9090.931
IV0.1990.1520.9640.9640.3170.2490.9060.916
V0.1800.1120.9710.9710.2550.1880.9390.943
VI0.1300.0880.9850.9850.2740.1840.9300.936
VII0.1990.1480.9640.9640.2980.2280.9170.933
VIII0.1110.0710.9890.9890.2700.1910.9320.936
IX0.0820.0470.9940.9940.2480.1450.9430.950
Table 6. The results of Station 3 (Rangpur) using SVR-based models.
Table 6. The results of Station 3 (Rangpur) using SVR-based models.
ModelsInput CombinationsTrainingTesting
RMSEMAENSER2RMSEMAENSER2
SVRI0.3690.2760.8310.8310.5160.4120.6510.773
II0.2810.1800.9020.9020.3900.3160.8000.887
III0.2400.1690.9290.9290.3520.2300.8380.878
IV0.3480.2540.8500.8500.4460.3310.7390.787
V0.2050.1330.9480.9480.3920.2640.7980.879
VI0.2060.1490.9480.9480.3390.2220.8500.882
VII0.2970.1930.8900.8900.3230.2510.8630.876
VIII0.1710.1110.9640.9640.3060.1790.8770.890
IX0.1390.0860.9760.9760.2460.1430.9200.923
SVR-PSOI0.3450.2560.8520.8530.5180.4130.6490.773
II0.2500.1950.9220.9240.3830.3060.8080.885
III0.2280.1590.9360.9380.3200.2310.8660.896
IV0.2860.2340.8980.9050.4160.3250.7730.793
V0.1520.1090.9710.9710.3170.2110.8680.899
VI0.1670.1130.9660.9660.3230.2170.8630.884
VII0.2340.1660.9320.9320.3180.2450.8680.891
VIII0.1350.1000.9770.9770.2990.1840.8830.909
IX0.0990.0640.9880.9880.2280.1450.9320.936
SVR-
PSOGSA
I0.2530.1980.9210.9210.4900.3830.6850.793
II0.2350.1830.9320.9320.3890.3160.8020.888
III0.1630.1200.9670.9670.2900.2090.8890.909
IV0.2550.1980.9190.9190.3970.3090.7930.804
V0.0890.0620.9900.9900.2380.1690.9260.945
VI0.1060.0780.9860.9860.2960.1920.8850.904
VII0.1490.1090.9730.9730.3230.2510.8630.895
VIII0.1220.0840.9810.9810.2620.1880.9100.942
IX0.0980.0720.9880.9880.2420.1770.9230.943
SVR-
PSOGWO
I0.2340.1780.9320.9320.5240.4070.6400.795
II0.1790.1380.9600.9600.3850.3140.8050.890
III0.1500.1090.9720.9720.2940.2100.8860.918
IV0.2040.1570.9480.9480.3910.2950.8000.825
V0.0860.0590.9910.9910.2430.1840.9220.939
VI0.1000.0710.9880.9880.3420.2370.8470.898
VII0.1720.1280.9630.9630.3250.2470.8610.905
VIII0.1060.0700.9860.9860.2780.2030.8990.939
IX0.0410.0290.9980.9980.2000.1320.9480.951
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ikram, R.M.A.; Mostafa, R.R.; Chen, Z.; Islam, A.R.M.T.; Kisi, O.; Kuriqi, A.; Zounemat-Kermani, M. Advanced Hybrid Metaheuristic Machine Learning Models Application for Reference Crop Evapotranspiration Prediction. Agronomy 2023, 13, 98. https://doi.org/10.3390/agronomy13010098

AMA Style

Ikram RMA, Mostafa RR, Chen Z, Islam ARMT, Kisi O, Kuriqi A, Zounemat-Kermani M. Advanced Hybrid Metaheuristic Machine Learning Models Application for Reference Crop Evapotranspiration Prediction. Agronomy. 2023; 13(1):98. https://doi.org/10.3390/agronomy13010098

Chicago/Turabian Style

Ikram, Rana Muhammad Adnan, Reham R. Mostafa, Zhihuan Chen, Abu Reza Md. Towfiqul Islam, Ozgur Kisi, Alban Kuriqi, and Mohammad Zounemat-Kermani. 2023. "Advanced Hybrid Metaheuristic Machine Learning Models Application for Reference Crop Evapotranspiration Prediction" Agronomy 13, no. 1: 98. https://doi.org/10.3390/agronomy13010098

APA Style

Ikram, R. M. A., Mostafa, R. R., Chen, Z., Islam, A. R. M. T., Kisi, O., Kuriqi, A., & Zounemat-Kermani, M. (2023). Advanced Hybrid Metaheuristic Machine Learning Models Application for Reference Crop Evapotranspiration Prediction. Agronomy, 13(1), 98. https://doi.org/10.3390/agronomy13010098

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