Establishing Coupled Models for Estimating Daily Dew Point Temperature Using Nature-Inspired Optimization Algorithms

: Potential of a classic adaptive neuro-fuzzy inference system (ANFIS) was evaluated in the current study for estimating the daily dew point temperature (Tdew). The study area consists of two stations located in Iran, namely the Rasht and Urmia. The daily Tdew time series of the studied stations were modeled through the other effective variables comprising minimum air temperature (Tmin), extraterrestrial radiation (Ra), vapor pressure deﬁcit (VPD), sunshine duration (n), and relative humidity (RH). The correlation coefﬁcients between the input and output parameters were utilized to determine the most effective inputs. Furthermore, novel hybrid models were proposed in this study in order to increase the estimation accuracy of Tdew. For this purpose, two optimization algorithms named bee colony optimization (BCO) and dragonﬂy algorithm (DFA) were coupled on the classic ANFIS. It was concluded that the hybrid models (i.e., ANFIS-BCO and ANFIS-DFA) demonstrated better performances compared to the classic ANFIS. The full-input pattern of the coupled models, speciﬁcally the ANFIS-DFA, was found to present the most accurate results for both the selected stations. Therefore, the developed hybrid models can be proposed as alternatives to the classic ANFIS to accurately estimate the daily Tdew.


Introduction
Dew point temperature (Tdew) as a weather parameter can be described as a temperature in which the air is saturated with water. In other words, the air is cooled under constant pressure to become saturated [1]. Tdew, along with relative humidity, can be considered as an index to recognize the content of air moisture [2]. The knowledge and availability of Tdew with reliable precision plays remarkable roles in various research areas such as agriculture, climatology, hydrology, and agronomy [3][4][5][6].
Tdew is an essential parameter for the estimation of evaporation and evapotranspiration [7]. Dew is an important factor for plant survival, specifically in arid regions with infrequent amounts of precipitation [8]. Tdew and wet-bulb temperature could be applied to compute the ambient temperature, which makes it possible to protect crops against frosts in which plants may acquire damage [9][10][11][12]. Climate change in each area can also be recognized by assessing Tdew trends over time [11].
Unlike weather parameters such as air temperature, relative humidity, precipitation, etc. that are usually recorded at many meteorological stations, Tdew information may not be measured and is therefore unavailable. Hence, alternative approaches should be taken into consideration to estimate reliable Tdew data. The use of artificial intelligence (AI)-based models has recently attracted attention in estimating Tdew, of which some are briefly reviewed below.
Hydrology 2022, 9, 9 2 of 15 The hourly Tdew data of 20 stations located in Georgia, United States of America, were estimated by Shank et al. [10] through artificial neural networks (ANN). Other weather variables were used as inputs to feed into the ANN. It was concluded that it could be considered as a dependable technique if Tdew data are not available. Zounemat-Kermani [13] tested the potential of a multiple linear regression (MLR) and feed-forward type of ANN in hourly Tdew estimation of a station located in Ontario, Canada. It was reported that the feed-forward ANN showed superiority over the conventional MLR. Nadig et al. [14] evaluated the performance of ANN in the estimation of hourly air and dew temperatures in Georgia, USA. Mohammadi et al. [11] compared the accuracy of three AI-based models in the estimation of daily Tdew of two stations in Iran. An extreme learning machine (ELM) was the best-performing model at both stations when comparing its performance with a support vector machine (SVM) and ANN. Kim et al. [15] estimated the time series of daily Tdew of two stations located in California, USA, utilizing the generalized regression neural networks (GRNN) and multi-layer perceptron (MLP). They found that these models could be considered as promising techniques to accurately estimate daily Tdew. Baghban et al. [16] used a genetic algorithm (GA) for optimizing the ANFIS and least square support vector machine (LSSVM) parameters in estimating Tdew and concluded that the GA could increase the performance of classic models. The classic ANFIS was applied by Mohammadi et al. [1] to estimate the daily Tdew of two locations in Iran, using other available weather variables, and they reported the reliable capability of ANFIS. Amirmojahedi et al. [17] recommended a coupled method by hybridizing the wavelet transform (WT) and an ELM for estimating the daily Tdew of a station in Iran. Its performance was also compared with the classic SVM and ANN. The best estimates of the daily Tdew were observed via the hybrid WT-ELM. Mehdizadeh et al. [6] assessed the performance of gene expression programming (GEP) when estimating daily Tdew of two sites in Iran and reported its dependable potential for estimating Tdew by applying other available weather parameters. The AI-based models, including the GEP, SVM, and M5 model tree, were used by Qasem et al. [18] in the estimation of the daily Tdew of Tabriz site, Iran. The M5 model tree outperformed the other models used. Dong et al. [12] applied seven classic AI-based models and three hybrid models, including an ELM optimized by three optimization algorithms, namely Bat, particle swarm optimization (PSO), and GA. The Yangling station was considered as a case study. They verified the abilities of mentioned optimizers to increase the potential of ELM when estimating the daily Tdew. Besides, the overall results of other studies such as Bickici Arikan et al. [19], Golizadeh Akhlaghi et al. [20], and Alizamir et al. [21] verified the capabilities of AI-based approaches for estimating Tdew.
In recent years, many efforts have been made to propose and implement the hybrid models as alternatives to the classic AI models when modeling meteorological, hydrological, and environmental variables [22][23][24][25]. It has also been well documented that the hybrid models can outperform the classic AI ones in estimating the times series of the investigated target variables [26]. In recent years, bio-inspired optimizers have received remarkable attention to developing the coupled methods created via coupling the AI and bio-inspired optimization techniques. Optimization is a vital engineering issue, and hence, various algorithms have been introduced with metaheuristics that have recently had many applications in engineering sciences [27]. Lately, the application of nature-inspired optimization algorithms has proven to be a significant solution for reducing the complexity of hydrological modeling [28]. These hybrid metaheuristic models can automatically provide powerful parameter estimation of the artificial intelligence models and enhance hydrological modeling performance; various nature-inspired algorithms have been employed for overcoming the limitation of the machine learning models for hydrological forecasting tasks [29][30][31].
The current study applied two nature-inspired optimization algorithms (BCO and DFA) for changing the attributes of the ANFIS model, such as learning rate and weights, in order to have a capable training and to reduce the modeling error. In addition, the ability of coupled ANFIS-BCO and ANFIS-DFA approaches were compared with the classical ANFIS for simulating dew point temperature to reach an efficient simulated result. The proposed models reach the optimal conditions when a suitability occurs between weights of ANFIS and the BCO and DFA, or the maximum number of iterations satisfies. This study seeks to estimate daily Tdew of two sites in Iran, through the classic adaptive neuro-fuzzy inference system (ANFIS) and novel hybrid models via coupling the ANFIS with the bio-inspired optimization algorithms consisting of the bee colony optimization (BCO), and dragonfly algorithm (DFA). The chief purpose of the present study is to enhance the estimation performance of Tdew via hybrid models. The effective variables of Tdew were recognized utilizing the correlation analysis. There are several hybrid models in the literature on the hybridization of AI and bio-inspired optimizers to estimate Tdew. Additionally, no studies have been reported in previous studies on proposing and testing the efficiency of coupled ANFIS-BCO and ANFIS-DFA techniques when estimating Tdew. Hence, this study is the first attempt to develop the aforementioned hybrid models.

Study Region and Data Used
Two stations were selected in Iran, Rasht and Urmia, as the study regions to achieve the objectives of this study (i.e., estimating the daily Tdew via other available weather data). The Rasht station in northern Iran includes a humid climate classification, having a latitude of 37 • 19 N, longitude of 49 • 37 E, and an elevation of −8.6 m. Furthermore, the Urmia station, with a semi-arid climate, is located in northwestern Iran. Its latitude, longitude, and elevation are 37 • 40 N, 45 • 03 E, and 1328.0 m, respectively. Figure 1 illustrates the geographical locations of the study locations on the map of Iran.
The current study applied two nature-inspired optimization algorithms (BCO and DFA) for changing the attributes of the ANFIS model, such as learning rate and weights, in order to have a capable training and to reduce the modeling error. In addition, the ability of coupled ANFIS-BCO and ANFIS-DFA approaches were compared with the classical ANFIS for simulating dew point temperature to reach an efficient simulated result. The proposed models reach the optimal conditions when a suitability occurs between weights of ANFIS and the BCO and DFA, or the maximum number of iterations satisfies. This study seeks to estimate daily Tdew of two sites in Iran, through the classic adaptive neuro-fuzzy inference system (ANFIS) and novel hybrid models via coupling the ANFIS with the bio-inspired optimization algorithms consisting of the bee colony optimization (BCO), and dragonfly algorithm (DFA). The chief purpose of the present study is to enhance the estimation performance of Tdew via hybrid models. The effective variables of Tdew were recognized utilizing the correlation analysis. There are several hybrid models in the literature on the hybridization of AI and bio-inspired optimizers to estimate Tdew. Additionally, no studies have been reported in previous studies on proposing and testing the efficiency of coupled ANFIS-BCO and ANFIS-DFA techniques when estimating Tdew. Hence, this study is the first attempt to develop the aforementioned hybrid models.

Study Region and Data Used
Two stations were selected in Iran, Rasht and Urmia, as the study regions to achieve the objectives of this study (i.e., estimating the daily Tdew via other available weather data). The Rasht station in northern Iran includes a humid climate classification, having a latitude of 37°19′ N, longitude of 49°37′ E, and an elevation of −8.6 m. Furthermore, the Urmia station, with a semi-arid climate, is located in northwestern Iran. Its latitude, longitude, and elevation are 37°40′ N, 45°03′ E, and 1328.0 m, respectively. Figure 1 illustrates the geographical locations of the study locations on the map of Iran. The required data for this study consist of dew point temperature (Tdew), minimum air temperature (Tmin), maximum air temperature (Tmax), mean air temperature (Tmean), soil temperature at 5 cm depth (Ts), sunshine duration (n), relative humidity The required data for this study consist of dew point temperature (Tdew), minimum air temperature (Tmin), maximum air temperature (Tmax), mean air temperature (Tmean), soil temperature at 5 cm depth (Ts), sunshine duration (n), relative humidity (RH), wind speed (Ws) and precipitation (P), which were collected from the Iran Meteorological Organization (IMO) during 2006-2017. Table 1 summarizes some of the statistical properties of the data used for both stations. Ra and VPD indicate the extraterrestrial radiation and vapor pressure deficit, respectively, which were calculated using the theoretical relationships proposed by Allen et al. [31]. The daily time series of Tdew for both the training and testing Hydrology 2022, 9, 9 4 of 15 periods at the studied stations are demonstrated in Figure 2. As can be seen, there are approximately similar trends for the daily Tdew time series of both the training and testing periods for the study locations. extraterrestrial radiation and vapor pressure deficit, respectively, which were calculated using the theoretical relationships proposed by Allen et al. [31]. The daily time series of Tdew for both the training and testing periods at the studied stations are demonstrated in Figure 2. As can be seen, there are approximately similar trends for the daily Tdew time series of both the training and testing periods for the study locations.   The data were split into training and testing datasets. The first nine-year data from 2006 to 2014 (i.e., 75% of the entire data) were applied when training models were developed. Moreover, the last three-year data between 2015 and 2017 (i.e., 25% of the entire data) were used for testing the models.

Adaptive Neuro-Fuzzy Inference System (ANFIS)
ANFIS as a machine learning approach is used as the base model in the current study. ANFIS used is according to the Takagi-Sugeno fuzzy system method, and this technique was proposed in 1993 by Jang et al. [32]. Considering ANFIS as an integrated neural network and fuzzy logic concept, it can take advantage of both of them in a single frame. The appropriate structure of the ANFIS technique is selected according to the input data, membership degree, input and output membership functions. During the training stage, input values can be closer to measured values by adjusting membership degree parameters according to the acceptable error rate. The ANFIS approach utilizes neural network learning algorithms and fuzzy logic for handling the nonlinear relationship between input and output values, and it has excellent potential in time series simulation and classification issues [33]. Numerical information can be extracted from a formally formulated rule according to the advantages of fuzzy rules, and the ANFIS can simply analyze the complex transformation of human intelligence into fuzzy systems. Its inference system is in accordance with the set of fuzzy if-then laws that can be learned to approximate nonlinear functions. Therefore, ANFIS is considered as a popular estimator tool in engineering fields. In this study, Fuzzy Cluster Means (FCM) method was used for fuzzy clustering based on the k-means clustering algorithm and the Backpropagation method was used for finding the best optimization method in classic ANFIS [34]. The fuzzy c-means (FCM) clustering by extracting a set of rules used for generating the fuzzy inference system (FIS) was applied in the current study for implementing the ANFIS model, and the optimal parameters of the ANFIS model are provided in Table 2. The BCO was recently introduced as an optimization technique that simulates the behavior of a bee colony [35]. A bee colony, a swarm of birds, or bees are simple examples of a demographic system. Another example is the collective intelligence of the bee colony around the hive. The BCO is an optimization algorithm based on the intelligent behavior of the bee population. In this mathematical model, artificial bee colonies have three kinds of bees. By the effort of worker bees, food-collecting occurs, and they bring it to the apiary from another food source. Determining the value of a resource of food is part of the duty of observer bees, and they patrol between worker bees to understand the worth of each food source. The food source can be considered as a state in the search space for a solution to the optimization problem in the BCO algorithm, and at the beginning, the number of food sources are considered as the number of bees in the apiary. The value of the objective function can determine the quality of food resources at each level of the algorithm running according to the amount of fitness [36].

Dragonfly Algorithm (DFA)
A new bio-inspired optimization technique is proposed, called the Dragonfly algorithm, as introduced by Mirjalili [37]. The basic idea of the DFA is inspired by the behavior of dragonflies in the nature. There are two basic stages in optimization: exploration, and exploitation. They are modeled on (teamwork) leading, searching for food, and avoiding enemies in the static and dynamic modes. The basic concepts of the Dragonfly algorithm include: (i) Separation: refers to avoiding one person dealing with other neighbors. (ii) Alignment: indicates the speed adjustment of people concerning other neighbors. (iii) Cohesion: refers to the tendency of individuals toward the center of neighborly imitation. The main goal of every mass and group is to maintain survival; thus, everyone must be attracted to food sources and distracted from enemies [37].

Models' Development
In the current study, BCO and DFA approaches were employed for enhancing the precision of the ordinary ANFIS model by optimizing its parameters. In the ordinary ANFIS model, parameters are nominated by a trial-and-error method, and sometimes this mentioned traditional method cannot select real optimal parameters of ANFIS. This is a reason for having more errors in the modeling process by ordinary ANFIS (because the ordinary ANFIS in some complex conditions cannot find optimal values by a traditional trial-and-error method). By hybridizing the nature-inspired optimization algorithms with an ordinary ANFIS model, the new proposed boosted hybrid models (namely ANFIS-BCO and ANFIS-DFA) can find the optimal values better than standalone ANFIS. The BCO and DFA algorithms were applied for this goal, and they try to find optimal parameters and then establish their updated membership functions based on the minimum error in Tdew estimation. Substantially, the suggested method can upgrade the abilities of the ordinary ANFIS model for time series prediction as well daily Tdew simulations. Then, the weight of the ANFIS model can be updated according to the gained values from the nature-inspired algorithms. The maximum number of iterations for running both optimization algorithms were considered as 500, and the number of search agents was considered as 40 in DFA. In addition, for running BCO, population size (colony size), number of decision variables, and acceleration coefficient upper bound were all considered as 100, 5, and 1, respectively. Figure 3 shows the schematic diagram of the modeling process by the classic ANFIS and hybrid models in this research. exploration, and exploitation. They are modeled on (teamwork) leading, searching for food, and avoiding enemies in the static and dynamic modes. The basic concepts of the Dragonfly algorithm include: (i) Separation: refers to avoiding one person dealing with other neighbors. (ii) Alignment: indicates the speed adjustment of people concerning other neighbors. (iii) Cohesion: refers to the tendency of individuals toward the center of neighborly imitation. The main goal of every mass and group is to maintain survival; thus, everyone must be attracted to food sources and distracted from enemies [37].

Models' Development
In the current study, BCO and DFA approaches were employed for enhancing the precision of the ordinary ANFIS model by optimizing its parameters. In the ordinary ANFIS model, parameters are nominated by a trial-and-error method, and sometimes this mentioned traditional method cannot select real optimal parameters of ANFIS. This is a reason for having more errors in the modeling process by ordinary ANFIS (because the ordinary ANFIS in some complex conditions cannot find optimal values by a traditional trial-and-error method). By hybridizing the nature-inspired optimization algorithms with an ordinary ANFIS model, the new proposed boosted hybrid models (namely ANFIS-BCO and ANFIS-DFA) can find the optimal values better than standalone ANFIS. The BCO and DFA algorithms were applied for this goal, and they try to find optimal parameters and then establish their updated membership functions based on the minimum error in Tdew estimation. Substantially, the suggested method can upgrade the abilities of the ordinary ANFIS model for time series prediction as well daily Tdew simulations. Then, the weight of the ANFIS model can be updated according to the gained values from the nature-inspired algorithms. The maximum number of iterations for running both optimization algorithms were considered as 500, and the number of search agents was considered as 40 in DFA. In addition, for running BCO, population size (colony size), number of decision variables, and acceleration coefficient upper bound were all considered as 100, 5, and 1, respectively. Figure 3 shows the schematic diagram of the modeling process by the classic ANFIS and hybrid models in this research.

Error Metrics Used to Evaluate the Models' Performance
Four statistical measures, namely root mean square error (RMSE), relative RMSE (RRMSE), mean absolute error (MAE), and coefficient of determination (R 2 ), were utilized in the current work to evaluate the models' performances in estimating the daily Tdew of the studied sites. They can be formulated as the following equations [38,39]: where T m,i and T e,i indicate the ith daily measured and estimated Tdew, respectively, T m and T e denote the mean of measured and estimated values of daily Tdew, respectively, and N is the total number of measured data used to train and test the models separately. It is obvious that lower values achieved for the RMSE, RRMSE, and MAE, and higher values for the R 2 , illustrate a better potential of each model in estimating daily Tdew.

Results and Discussions
The values of daily Tdew for two locations selected in Iran were first estimated through the classic ANFIS. To determine the significance level of the weather variables on Tdew estimation, correlation coefficients existing among the input predictors and output variables were taken into consideration. The hybrid models were then proposed to improve the daily Tdew estimation using the hybridization of classic ANFIS with two types of bio-inspired optimization algorithms consisting of BCO and DFA.

Performance Investigation of the Classic and Hybrid Models Proposed
To better feed the models developed in this study, correlation coefficients of the inputs and target variables were calculated as graphically illustrated in Figure 4. According to this radar graph, temperature components (i.e., Tmin, Tmax, Tmean, and Ts) illustrate the highest correlations with Tdew. Among them, Tmin followed by Tmean relatively show higher correlations than those of the Ts and Tmax. After the temperature parameters, the highest correlations are related to Ra and VPD for both the Rasht and Urmia stations. After the variables mentioned above, n is a high-correlated variable with Tdew. In addition, there exist weakest correlations for Tdew with the Ws and P variables. High-correlated variables with the Tdew were applied to define the input combinations when estimating the daily Tdew. In this context, Tmin from the temperature-based components, and Ra, VPD, n, and RH were utilized as effective variables. The input combinations considered to implement the classic and hybrid models are tabulated in Table 3. As is apparent, the first five models represent the sole-input combinations.
Hydrology 2022, 9, 9 8 of 16 components, and Ra, VPD, n, and RH were utilized as effective variables. The input combinations considered to implement the classic and hybrid models are tabulated in Table 3. As is apparent, the first five models represent the sole-input combinations.  The statistical results of the classic ANFIS, including the RMSE, RRMSE, MAE, and R 2 metrics at the Rasht and Urmia stations are provided in the first parts of Tables 4 and 5, respectively. As expected, Tmin had the best accuracy of Tdew estimation at both stations. Ra followed by VPD demonstrated higher performances than the n and RH variables. It can be therefore said that the input variables with the higher correlation coefficients with the target variable (i.e., Tdew in this study) could provide reasonable  The statistical results of the classic ANFIS, including the RMSE, RRMSE, MAE, and R 2 metrics at the Rasht and Urmia stations are provided in the first parts of Tables 4 and 5, respectively. As expected, Tmin had the best accuracy of Tdew estimation at both stations. Ra followed by VPD demonstrated higher performances than the n and RH variables. It can be therefore said that the input variables with the higher correlation coefficients with the target variable (i.e., Tdew in this study) could provide reasonable estimation accuracy. The high-performance variable, i.e., Tmin was assumed as a fixed variable, and other variables were added into Tmin for defining the input combinations with more inputs. For example, the second-best performing variable (Ra) was added to Tmin to develop the input combination with two inputs. Besides, the last models in Table 3 illustrate the full-input combination. The accuracy of ANFIS6 models at both stations is somewhat similar to that of the ANFIS1 models, which means that adding Ra to the Tmin in ANFIS6 could not improve the performance of ANFIS1. However, considering the Tmin, Ra, and VPD in the ANFIS7 model of Urmia station considerably increased the accuracy of Tdew estimation, while no significant difference was observed in the performance of ANFIS6 and ANFIS7 models developed at the Rasht station. Finally, the ANFIS8 (four-input) and ANFIS9 (full-input) models performed better than the models with fewer inputs at both the studied locations; however, the ANFIS9 was the best-performing classic ANFIS with the highest degree of accuracy (Rasht station; RMSE = 1.  Table 4. Statistical metrics of the RMSE, RRMSE, MAE, and R 2 for the classical and hybrid models developed at the Rasht station.

Models
Training Testing  As already mentioned, the main purpose of this study is to increase the estimation accuracy of daily Tdew via the novel hybrid models. In doing so, two different types of optimization techniques, including BCO and DFA, were applied to develop and propose the hybrid ANFIS-BCO and ANFIS-DFA models. Error statistics of hybrid models are given in the second and third sections of Tables 4 and 5. The same outcomes achieved for the classic ANFIS were also obtained for the coupled techniques implemented, i.e., the better estimation accuracy of the Tdew at both the stations through the Tmin followed by Ra and VPD compared to the n and RH variables, improving the performance of all the hybrid ANFIS1-BCO and ANFIS1-DFA models with the inclusion of other effective variables (i.e., Ra, VPD, n, RH) when estimating the Tdew and the most accurate Tdew estimates considering the full-input combination.
The full-input combination considers the other weather variables such as Tmin, Ra, VPD, n, and RH to estimate the target variable (Tdew). Among these variables, Tmin, n, and RH are usually recorded at most weather stations. Furthermore, Ra and VPD are calculated by the other available parameters. Ra is computed by the Julian day and latitude of the station, while the Tmin, Tmax, Tmean, and RH are needed to calculate the VPD. It can therefore be concluded that all the hybrid models proposed in this study, especially the ANFIS-DFA, can be used for estimating the daily Tdew by utilizing the predictors as mentioned above with a dependable degree of precision.

Performance Comparison of Classic and Coupled Models Proposed
A performance evaluation of the classic ANFIS and hybrid models developed in the current study focusing on error metrics of RMSE, RRMSE, MAE, and R 2 , tabulated in Tables 4 and 5, clearly illustrates that the novel coupled methods outperformed the classic ANFIS to estimate the daily Tdew time series of the study sites due to the lower values for the RMSE, RRMSE, and MAE, and higher values for the R 2 criterion. This can be verified, taking into consideration the ability of the new hybrid models in Tdew estimation. This study shows that the optimization method can help to solve engineering problems, especially BCO and DFA, which can improve the ability of ANFIS as an estimator tool for estimating Tdew at different climates by different input variables.
Regarding the capabilities of all the models developed for estimating daily Tdew at both locations with different climates, it can be apparently observed that all the classic and coupled approaches present better accuracies in the Urmia station than those of the Rasht station under sole-input combinations. Conversely, input combinations with more inputs, from two to five, in all the classic and coupled models show superior outcomes at the Rasht station compared to the Urmia station. Different climates can therefore produce different results for the models. Climate difference at the two stations studied has affected the results of the models. The Rasht station has a humid climate due to its proximity to the world's largest lake (Caspian Sea), and this climatic complexity has made the estimation models at this station require more input variables.
Generally, among the hybrid models, the hybrid ANFIS-BCO under all the defined input combinations presented lower estimation accuracy. Meanwhile, the hybrid ANFIS-DFA performed better than the ANFIS and ANFIS-BCO and is therefore the best-performing hybrid model to accurately estimate the daily Tdew at both the Rasht and Urmia stations. It can be justified considering the point of the reliable capability of the DFA algorithm to find the optimal parameters of the ANFIS model. Also, the ability of DFA for finding any possible optimal solution helped to have a better performance in Tdew estimation. The results demonstrate DFA's applicability and robustness to solving real problems with unknown search spaces, and these results are consistent with the results of Mirjalili [37]. To schematically investigate the estimation potential of classic ANFIS and proposed hybrid techniques, scatter graphs for the measured and estimated Tdew data via any of the models were prepared and are shown in Figure 5 for the studied stations. To this, the best-performing classic ANFIS (i.e., ANFIS9) and hybrid models (i.e., ANFIS9-BCO and ANFIS9-DFA) were employed during the testing period. The dashed lines in the scatter plots indicate the 1:1 exact line. In general, the values of R 2 are closer to 1 in the hybrid models, specifically the hybrid ANFIS-DFA, than those of the classic ANFIS, which denotes the dependable performance of the hybrid approaches when estimating daily Tdew.
Besides the scatter plots, Taylor diagrams were also depicted for the superior classic and hybrid models (i.e., ANFIS9, ANFIS9-DFA, and ANFIS9-BCO) at both the studied locations, as shown in Figure 6. This diagram is made up of two semicircles, one for negative correlations and the other for positive correlations. Three statistics comprising the standard deviation, correlation coefficient, and centered RMSE can easily be observed in these diagrams. The measured point (i.e., black point) is defined, taking into account the standard deviation, which is depicted on the horizontal axis. The shorter distance related to each model with the measured point indicates the better performance of that model. In this context, ANFIS9-DFA (i.e., green points in Figure 6) is the best-performing model. Figure 5. Scatter graphs of the measured and estimated Tdew through the best-performing classical and hybrid models at the studied stations during the testing period. Figure 5. Scatter graphs of the measured and estimated Tdew through the best-performing classical and hybrid models at the studied stations during the testing period.

Conclusions
This study estimated the daily Tdew time series of Rasht and Urmia stations located in Iran by the classic ANFIS and its hybridized form with two bio-inspired optimization techniques, namely BCO and DFA. Other available weather variables were used as the predictors to feed the models developed. The impact of considered variables on Tdew was investigated via correlation analysis. The weather variables, including Tmin, Ra, VPD, n, and RH were used as inputs, and therefore nine different input combinations were defined. Based on the outcomes achieved, Tmin followed by Ra were the most effective

Conclusions
This study estimated the daily Tdew time series of Rasht and Urmia stations located in Iran by the classic ANFIS and its hybridized form with two bio-inspired optimization techniques, namely BCO and DFA. Other available weather variables were used as the predictors to feed the models developed. The impact of considered variables on Tdew was investigated via correlation analysis. The weather variables, including Tmin, Ra, VPD, n, and RH were used as inputs, and therefore nine different input combinations were defined. Based on the outcomes achieved, Tmin followed by Ra were the most effective variables to accurately estimate daily Tdew. Moreover, the n and RH variables showed the least impact on Tdew. Accordingly, all the classic and hybrid models under sole-input combinations presented the highest and lowest accuracy regarding the application of Tmin and RH, respectively. The performance of the methods developed in the present study was improved with the inclusion of more inputs such that the ANFIS9, ANFIS9-BCO, and ANFIS9-DFA with full inputs illustrated the best performance in each class of classic and hybrid models. For both the studied stations, the hybrid models led to better estimates of the daily Tdew in comparison to the classic ANFIS; however, the ANFIS-DFA was the best hybrid paradigm. It was concluded that the coupled methods surpassed the classic model used in this study (i.e., ANFIS) for the estimation of daily Tdew. There are some advantages in the optimization of the proposed methods, such as searching in response to global and local spaces and checking all possible responses around a potential, optimal candidate. In addition, such metaheuristic algorithms are capable, and small numbers of iterations are needed for solving the problem. Then, finding an optimal response using a strong search can make sure the ANFIS is trained by a capable method and that a minimum error occurs during the Tdew simulation process. It is therefore suggested that more hybrid models could be proposed through the hybridization of other classic models, such as SVR, MLP, etc., with different types of bio-inspired optimization techniques, including evolution-based algorithms, ecology-based algorithms, and swarm-intelligence-based algorithms for the precise estimation of a Tdew time series. Besides, the estimation accuracy of hybrid models implemented could be tested by estimating other temperature variables consisting of air temperature and soil temperature, as well as other hydrological and agro-meteorological variables such as evaporation, evapotranspiration, precipitation, and so on.