Next Article in Journal
Heavy Metal Contamination of Irrigation Water, Soil, and Vegetables and the Difference between Dry and Wet Seasons Near a Multi-Industry Zone in Bangladesh
Previous Article in Journal
Reliable Future Climatic Projections for Sustainable Hydro-Meteorological Assessments in the Western Lake Erie Basin
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Estimating Daily Dew Point Temperature Using Machine Learning Algorithms

1
Computer Science Department, College of Computer and Information Sciences, Al Imam Mohammad Ibn Saud Islamic University (IMSIU), Riyadh 11432, Saudi Arabia
2
Computer Science Department, Faculty of Applied Sciences, Taiz University, Taiz, Yemen
3
Department of Water Engineering, University of Tabriz, Tabriz 5166616471, Iran
4
School of the Built Environment, Oxford Brookes University, Oxford OX3 0BP, UK
5
Institute of Automation, Kando Kalman Faculty of Electrical Engineering, Obuda University, 1034 Budapest, Hungary
6
The Queensland University of Technology, Institute of Health and Biomedical Innovation, 60 Musk Avenue, Queensland 4059, Australia
7
Department for Management of Science and Technology Development, Ton Duc Thang University, Ho Chi Minh City, Vietnam
8
Faculty of Information Technology, Ton Duc Thang University, Ho Chi Minh City, Vietnam
9
Department of Civil and Environmental Engineering, Hong Kong Polytechnic University, Hong Kong, China
*
Author to whom correspondence should be addressed.
Water 2019, 11(3), 582; https://doi.org/10.3390/w11030582
Submission received: 4 February 2019 / Revised: 13 March 2019 / Accepted: 18 March 2019 / Published: 20 March 2019
(This article belongs to the Section Hydrology)

Abstract

:
In the current study, the ability of three data-driven methods of Gene Expression Programming (GEP), M5 model tree (M5), and Support Vector Regression (SVR) were investigated in order to model and estimate the dew point temperature (DPT) at Tabriz station, Iran. For this purpose, meteorological parameters of daily average temperature (T), relative humidity (RH), actual vapor pressure (Vp), wind speed (W), and sunshine hours (S) were obtained from the meteorological organization of East Azerbaijan province, Iran for the period 1998 to 2016. Following this, the methods mentioned above were examined by defining 15 different input combinations of meteorological parameters. Additionally, root mean square error (RMSE) and the coefficient of determination (R2) were implemented to analyze the accuracy of the proposed methods. The results showed that the GEP-10 method, using three input parameters of T, RH, and S, with RMSE of 0.96°, the SVR-5, using two input parameters of T and RH, with RMSE of 0.44, and M5-15, using five input parameters of T, RH, Vp, W, and S with RMSE of 0.37 present better performance in the estimation of the DPT. As a conclusion, the M5-15 is recommended as the most precise model in the estimation of DPT in comparison with other considered models. As a conclusion, the obtained results proved the high capability of proposed M5 models in DPT estimation.

1. Introduction

Dew point temperature (DPT) is defined as the temperature in which air becomes liquid water due to the high concentration of water molecules. Precise and accurate estimation of DPT has a significant role in solving agricultural problems, such as calculating the amount of available moisture in the air and estimating the near surface humidity [1]. DPT and relative humidity are commonly used to measure the air humidity level [2]. The DPT can also be used to estimate the temperature of crops considering glaciation [3]. Many studies have paid attention to the accurate estimation of DPT using regression methods. However, data-driven methods such as Gene Expression Programming (GEP) and Neuro-Fuzzy Inference System (ANFIS) have been developed to identify optimal functions and modeling for complex phenomena. In this regard, several studies have been carried out on the application of the mentioned methods in meteorological studies [4,5,6,7,8,9,10,11,12,13]. Shiri [2] compared the capabilities of the artificial neural network (ANN) and GEP to estimate the DPT using meteorological parameters at Seoul and Inchon stations, located in South Korea. They used two management scenarios: In the first scenario, the meteorological information of each station was used to estimate the DPT of the same station; in the second scenario, they used the meteorological information of adjacent stations. Their results showed that in both scenarios GEP was more accurate than ANN. Also, the application of the second scenario showed that GEP had more accurate results in estimating the DPT values of Seoul stations using Incheon station parameters. They also reported that the DPT values at Seoul Station could be estimated using the average temperature and relative humidity of the Incheon station with proper accuracy. Deka et al. [14] examined the ability of a support vector machine (SVM), ANN, and Extreme Learning Machine (ELM) to estimate DPT at two stations in Iran. They showed that the results of the ELM model were more similar to observed DPT at the two mentioned stations. In other research, Zounemat-Kermani [15] implemented two methods of multiple linear regression (MLR) and Levenberg–Marquardt algorithm (LMA) in the artificial neural network (LMA–ANN) in order to estimate DPT values at Ontario Station, Canada. The results of the LMA–ANN model had an appropriate match with observational data. Additionally, Jia et al. [16] investigated dew formation. For this purpose, they used meteorological data of average temperature, sunny hours, wind speed, saturated vapor pressure, relative humidity, and DPT values of three stations of Dagot, Pohang, and Ulsan, South Korea. They reported that the effects of sunny hours, wind speed, and saturated vapor pressure were lower than other parameters. Therefore, it was possible to estimate the DPT using average temperature and relative humidity. Attar et al. [17] used GEP, multivariate adaptive regression splines (MARS), and SVM models to estimate the DPT in arid regions of Iran. Using the meteorological data of 13 synoptic stations during the 55 years (1996 to 2014), and by defining 50 different scenarios. They concluded that the MARS model offers more accurate results than other studied models. In a similar study, Mehdizadeh et al. [18] estimated the DPT values in Tabriz and Urmia cities, in the northwest of Iran, using the GEP method. They defined three scenarios: A parameters-based scenario, a temperature-based scenario, and a periodicity-based scenario considered the meteorological parameters of minimum, maximum, and mean air temperature, actual vapor pressure, and atmospheric pressure. Their results showed that the actual vapor pressure is the most effective meteorological parameter in estimating the DPT in the study area.
Therefore, over the last decade, researchers have tried to estimate DPT values with suitable accuracy. For this reason, the main purpose of the current study was to implement three data-driven methods of GEP, M5, and SVR in order to improve the estimation accuracy and develop some mathematical formulations for obtaining precise estimations of DPT values using explicit formulations. To the best of our knowledge, the application of M5 has not been reported in the literature. In other words, the goals of the study were (i) evaluating the performance of the models above in the estimation of DPT, and (ii) investigating the role of climatic parameters estimation DPT values. The rest of the paper is structured as follows: Section 2 describes implemented methods, evaluation parameters, and characteristics of the study area. Additionally, Section 3 discussed the obtained results and, finally, the conclusion is presented in Section 4.

2. Study Area

The Tabriz synoptic station is one of the oldest meteorological stations belonging to the Iranian Meteorological Organization, located in East Azerbaijan province with a latitude of 38° 05′ N and longitude of 46° 17′ E and an elevation of 1364 m above sea level. In the current study, the DPT values of Tabriz station in the period of 1998 to 2016 were utilized to evaluate the precision of the considered models. The geographic location of the study area is shown in Figure 1.

3. Materials and Methods

3.1. Support Vector Regression (SVR)

SVM, which is established by statistical learning theory, has been broadly applied for identifying complex patterns of different environmental phenomena. Moreover, SVR, as one of the types of SVM, has been used previously for regression problems [20,21]. The basics of the support vector regression are presented briefly below.
If the available data can be linearly split, they can be distinguished using cloud-computing system data. In some cases, data cannot be linearly separated. In these cases, the data is mapped to a larger dimensional space, after which, the data are separated. Figure 2B shows the non-linear mapping of φ (0): Rn → Rnh, in which, the nonlinear training data are mapped to a higher-dimensional space. Following that, a linear relationship is obtained between the input and output data in the converted space. The general form of the linear relation is as follows:
f ( x ) = W T φ ( x ) + b
where, f(x) is the estimated variable, WT is the transpose of the vector of coefficients and b is a constant coefficient. The SVM is based on minimizing the value of empirical error.
R e m p ( f ) = 1 N i = 1 N Θ ε ( y i , W T φ ( x ) + b )
where, Θ ε (y, f(x)) is the ε-insensitive error function defined as Equation (3):
Θ ε ( y , f ( x ) ) = { | f ( x ) y | ε                     i f   | f ( x ) y | ε 0                                                                                   o t h e r w i s e            
In addition, the function Θ ε (y, f(x)) is used to find the optimal separator plain in a high dimension space (which may have infinitive dimensions). In a space with a high dimension, an optimal separator plate will maximize the distance between training data. Following this, the SVR model minimizes the general error concerning the constraints. These limitations presented in Equations (4) to (8).
min w , b , ξ , ξ R ε ( W , ξ , ξ ) = 1 2 W T W + C i = 1 N ( ξ i + ξ i )  
where, ξi is an error greater than –ε, ξi* is the error greater than , and c is a constant.
y i W T φ ( x i ) b ε + ξ i ,   i = 1 , 2 , , N
y i + W T φ ( x i ) + b ε + ξ i ,   i = 1 , 2 , , N
  ξ i 0 ,   i = 1 , 2 , ,
  ξ i 0 ,   i = 1 , 2 , , N
After obtaining the optimal separator plate, the vector of the coefficients w obtained as Equation (9).
W = i = 1 N ( β i β i ) φ ( x i   )  
in which, βi* and βi are calculated by applying quadratic programming and Lagrange coefficients, respectively. Finally, the SVR function is obtained as Equation (10) in a two-dimensional space, as follows:
f ( x ) = i = 1 N ( β i β i   ) K ( X i , X ) + b  
K (Xi, Xj) is called the Kernel function and is equal to the inner product of the two vectors of Xi and Xj in a space with a high dimension. In the SVR method, several Kernel functions are used, including the polynomial Kernel function, the normalized polynomial Kernel function, the radial-base Kernel function, and the Pearson Kernel function [22].

3.2. Gene Expression Programming

GEP, which is a variant of genetic programming (GP), is a generalized genetic algorithm. GEP is considered to be a circular method based on Darwin’s theory of evolution. GEP at the beginning of the process does not take into account the functional relationship and can optimize the structure of the model and its components [23]. Unlike the genetic algorithm, GEP acts on the tree structure of formulas rather than a series of binary numbers. The tree structures created from the set of functions (mathematical operators used in formulas) and terminals (problem variables and constant numbers). Before the implementation of GEP, the following factors are determined in the following:
  • Terminal set (problem variables, randomized constant numbers),
  • The mathematical operators used in formulas,
  • Select the fitness function (RMSE, MSE, MAE,…) to measure the fitness of the formulas,
  • Select the parameters controlling the implementation of the program (population size, the probability associated with the use of genetic operators and other details related to the implementation of the program),
  • The completion benchmark and the presentation of the results of the program implementation (the number of new population production, the determination of the specified amount for the fitness of the formulas if the fitness level is equal to or greater than that value stopped) [24]. The outlines of the mentioned steps are shown in Figure 3. Moreover, the parameters used in the implementation of the GEP presented in Table 1.

3.3. M5 Model Tree

The M5 model tree is a subset of the machine learning and data mining methods developed by Quinlan [25]. Data mining refers to the process of searching and discovering various models, summarizing and obtaining quantities from the collections. Learning machine and data mining methods have the ability to discover data patterns, semi-automatically. The main reasons for using a model tree are as follows:
  • The model tree is directly related to estimative variables; therefore, the results of the model are easy to understand.
  • Model trees are non-parametric, and there is no user intervention on them.
  • The output of the model has a high degree of accuracy that can be compared to other models.
The structure of the model trees includes roots, branches, nodes, and leaves. The nodes are represented by a circle and the branches represent the connection between the nodes [26]. The generation of the model tree structure consists of different steps of creating a tree and pruning it. In the first step, an inferential algorithm or division criterion is used for the production of the tree. The decision criterion for the M5 is the standard deviation of the class values, which is calculated as a quantity of error to a node and calculates the expected reduction in this error as the result of the test of each attribute in that node. The standard deviation ratio (SDR) is calculated as follows (Equation (11)):
S D R = s d ( T ) | T I | | T | s d ( T I )
In which, T is a collection of input samples to each node, Ti is a subset of the samples that have the i output of the potential series and sd denotes the standard deviation [27]. As a result of the branching process, the data in the child nodes has a lower standard deviation than the parent node and is purer. After maximizing all possible branches, the M5 selects an attribute that maximizes the expected reduction. This division forms a large quasi-tree structure, which causes over-fitting. To overcome the mentioned problem, the tree should be pruned. This is done by replacing a sub-tree with a leaf; therefore, the second step in the design of a model tree is to prune the grown tree and to replace the sub-trees with linear regression functions. This method for creating the model tree divides the space of the input parameters into smaller areas or sub-areas. In each area, a linear regression model is fitted. After obtaining a linear model, the simplification of the model can minimize the estimation error by deleting the model parameters [28].

3.4. Evaluation Criteria

Error values between observed and estimated data were determined by RMSE and coefficient of determination (R2). The RMSE was used to evaluate the accuracy of the estimations. The consistent estimated values of the model lead to the minimization of the statistical index. Furthermore, R2 is a statistical tool for determining the type and degree of the relationship of a variable with other variables. This coefficient varies from 0 to 1; when there is no relationship between two variables, its value is equal to zero [29,30]. Furthermore, Taylor diagrams [31] were used to check the accuracy of the mentioned models. It is noteworthy that Taylor suggested a diagram, in which measured parameters and some characteristics of the model are summed up, coincidentally. Surprisingly, Taylor diagrams utilize several points on a polar plot for comparing the accuracy of measured and estimated values. In these diagrams, the coefficient of determination and normalized standard deviation are represented by an azimuth angle and radial distances from the base point, respectively [31,32].

4. Results and Discussion

In order to reach the research objectives, daily average temperature, relative humidity, actual vapor pressure, wind speed, and sunny hours at the Tabriz synoptic station were collected from the Meteorological Organization of East Azerbaijan province, Iran during the period 1998 to 2016. The statistical characteristics of the implemented data are presented in Table 2. There is no basic way of separating training and testing data. For example, the study of Kurup and Dudani [33] used a total of 63% of their data for model development, whereas Samadianfard et al. [4] and Samadianfard et al. [7] used 67% of total data, and Deo et al. [34] used 70% of total data to develop their models. Thus, to develop the studied GEP, M5, and SVR models for estimation DPT, we divided the data into training (67%) and testing (33%). Therefore, the accuracy of the models in estimating DPT evaluated through Taylor diagrams. Additionally, the effects of considered meteorological parameters were inspected by defining 15 different input combinations (Table 3).
After performing the computations for different input combinations, the accuracy of the considered models was determined in the testing phase based on the statistical criteria (Equations (9) and (10)) and Taylor diagrams. The obtained results are presented in Table 4.
As can be seen in Table 4, GEP-10 with RMSE of 0.96 degrees and R2 equal to 0.902 with the parameters of T, RH, and S shows better performance compared to GEP models. However, SVR-6 with RMSE of 0.44 degree and R2 of 0.996 presents more accurate estimation compared to the SVR models. Furthermore, the best estimation of the DPT, based on M5 models, was related to M5-15 with RMSE of 0.37 degree and R2 of 0.996 and using all considered meteorological parameters as the input. In other words, a comprehensive comparison between the mentioned models exhibited that M5-15 had the best performance in estimation DPT values by using input combinations of T, S, RH, W, Vp. After selecting the most accurate models for estimation DPT values, the time series plots and scatterplots are finalized and illustrated in the Figure 4 and Figure 5.
It can be comprehended from Table 4 and Figure 4 that the estimation accuracy of the M5-15 was higher than the GEP-10 and SVR-6. The above-mentioned conclusion, regarding the high accuracy of the M5-15 model in estimation the DPT for Tabriz station, can be deduced from Figure 5. In this figure, it can be seen that the distribution of the points around the bisector line in the M5-15 model was less than the corresponding points of GEP-10 and SVR-6.
Furthermore, Taylor charts were used to examine the standard deviation and correlation values among estimated and measured DPT values for the GEP, M5, and SVR models with different input parameters. Taylor diagrams for models mentioned above are shown in Figure 6. The length of the space from the reference point (a green color point) to each point is defined as centered RMSE [31]. Therefore, the most accurate model has a minimum distance between the green point and its corresponding point. According to Figure 6, M5-15 (a blue color point) offered the most accurate estimations of DPT values at Tabriz station.
One of the advantages of GEP and the M5 models, in comparison with other data-driven methods, is their ability to provide explicit relationships to calculate the output parameter. Therefore, for the current study, Equation (12) was obtained for estimation DPT values using GEP-10 as the most accurate GEP model.
T d e w = R H 9.7 9.7 e ( S R H ) + T ( e ( 3.4 R H ) 2 ) 3 + R H 9.7 9.7 S R H
Additionally, the list of linear equations (presented in Table 5) was the outcome of M5-15 estimation of the DPT values using meteorological parameters of T, S, RH, W, Vp.
As previously mentioned, Deka et al. [14] used SVM, ANN, and ELM and implemented meteorological parameters of minimum, maximum, and average temperatures, relative humidity, atmospheric pressure, water vapor pressure, sunny hours, and solar radiation in order to estimate the DPT in two cities of Kerman province, Iran. The minimum RMSE reported in the mentioned study was 0.49, related to the ELM method by using minimum temperature and water vapor pressure data as input parameters. In the present study, with the application of the M5 and applying meteorological parameters of average temperature, relative humidity, actual vapor pressure, wind speed, and sunny hours, the RMSE decreased to 0.37, indicating the high accuracy of the M5 model tree for estimation DPT values. The output of the M5 was a simple linear relationship that can be used to calculate the DPT values easily, while the ELM does not have such a capability. Additionally, in another study by Baghban et al. [35] the maximum estimation accuracy of the DPT was reported using an SVM model with an RMSE value of 0.4.
Furthermore, the precision of the M5-15 in the current study was more than the accuracy of the proposed SVM method by Baghban et al. [35]. To investigate the influence of input parameters on the DPT estimation, the RMSE and R2 were utilized for different groupings of input variables. For this purpose, all utilized models, including GEP, M5, and SVR were selected for sensitivity analysis (Table 6). Each model confirmed the extent to which the eliminated variable would affect the model accuracy. As shown in Table 6, the precision of all models decreased if each of T, RH, Vp, W, and S input parameters were removed from the modeling. Furthermore, it can be comprehended that T had the greatest effect in increasing the prediction accuracy. In other words, eliminating T caused a sharp increase in RMSE values in all studied models.

5. Conclusions

In the current study, three data-driven methods including GEP, M5, and SVR were used to estimate DPT values at Tabriz synoptic station, Iran. For this purpose, the meteorological parameters were collected from the Meteorological Organization of East Azerbaijan province in the period 1998 to 2016. Also, 15 different input combinations were defined to study the effect of meteorological parameters on the estimation of DPT values. The results of this study revealed that the SVR-6, using two input parameters of T and RH, and GEP-10 using three parameters of T, RH, and S, had appropriate performance in the estimation of DPT values. Furthermore, the overall analysis of the studied methods showed that the M5-15 using five parameters of T, S, RH, W, and Vp had the best performance in the estimation of DPT values at Tabriz station in comparison with all considered models with different input combinations. To conclude, M5-15 is proposed as the most accurate method for the estimation of DPT values at the Tabriz synoptic station, Iran.

Author Contributions

The authors have made equal contributions.

Funding

This research received no external funding.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Famiglietti, C.A.; Fisher, J.B.; Halverson, G.; Borbas, E.E. Global Validation of MODIS Near-Surface Air and Dew Point Temperatures. Geophys. Res. Lett. 2018, 45, 7772–7780. [Google Scholar] [CrossRef]
  2. Shiri, J. Prediction vs. estimation of dewpoint temperature: Assessing GEP, MARS and RF models. Hydrol. Res. 2018. [Google Scholar] [CrossRef]
  3. Ali, H.; Fowler, H.J.; Mishra, V. Global observational evidence of strong linkage between dew point temperature and precipitation extremes. Geophys. Res. Lett. 2018, 45, 12320–12330. [Google Scholar] [CrossRef]
  4. Samadianfard, S.; Delirhasannia, R.; Kisi, O.; Agirre-Basurko, E. Comparative analysis of ozone level prediction models using gene expression programming and multiple linear regression. GEOFIZIKA 2013, 30, 43–74. [Google Scholar]
  5. Mosavi, A.; Ozturk, P.; Chau, K.W. Flood prediction using machine learning models: Literature review. Water 2018, 10, 1536. [Google Scholar] [CrossRef]
  6. Samadianfard, S.; Nazemi, A.H.; Sadraddini, A.A. M5 model tree and gene expression programming based modeling of sandy soil water movement under surface drip irrigation. Agric. Sci. Dev. 2014, 3, 178–190. [Google Scholar]
  7. Samadianfard, S.; Sattari, M.T.; Kisi, O.; Kazemi, H. Determining flow friction factor in irrigation pipes using data mining and artificial intelligence approaches. Appl. Artif. Intell. 2014, 28, 793–813. [Google Scholar] [CrossRef]
  8. Dehghani, M.; Riahi-Madvar, H.; Hooshyaripor, F.; Mosavi, A.; Shamshirband, S.; Zavadskas, E.K.; Chau, K.W. Prediction of Hydropower Generation Using Grey Wolf Optimization Adaptive Neuro-Fuzzy Inference System. Energies 2019, 12, 289. [Google Scholar] [CrossRef]
  9. Lee, O.; Kim, S. Estimation of Future Probable Maximum Precipitation in Korea Using Multiple Regional Climate Models. Water 2018, 10, 637. [Google Scholar] [CrossRef]
  10. Jabbari, A.; Bae, D.H. Application of Artificial Neural Networks for Accuracy Enhancements of Real-Time Flood Forecasting in the Imjin Basin. Water 2018, 10, 1626. [Google Scholar] [CrossRef]
  11. Samadianfard, S.; Asadi, E.; Jarhan, S.; Kazemi, H.; Kheshtgar, S.; Kisi, O.; Sajjadi, S.; Abdul Manaf, A. Wavelet neural networks and gene expression programming models to estimate short-term soil temperature at different depths. Soil Tillage Res. 2018, 175, 37–50. [Google Scholar] [CrossRef]
  12. Nie, J.; Liu, J.; Li, N.; Meng, X. Dew point measurement using dual quartz crystal resonator sensor. Sens. Actuators B Chem. 2017, 246, 792–799. [Google Scholar] [CrossRef]
  13. Shamshirband, S.; Jafari Nodoushan, E.; Adolf, J.E.; Abdul Manaf, A.; Mosavi, A.; Chau, K.W. Ensemble models with uncertainty analysis for multi-day ahead forecasting of chlorophyll a concentration in coastal waters. Eng. Appl. Comput. Fluid Mech. 2019, 13, 91–101. [Google Scholar] [CrossRef]
  14. Deka, P.C.; Patil, A.P.; Yeswanth Kumar, P.; Naganna, S.R. Estimation of dew point temperature using SVM and ELM for humid and semi-arid regions of India. ISH J. Hydraul. Eng. 2018, 24, 190–197. [Google Scholar] [CrossRef]
  15. Zounemat-Kermani, M. Hourly predictive Levenberg–Marquardt ANN and multi linear regression models for predicting of dew point temperature. Meteorol. Atmos. Phys. 2012, 117, 181–192. [Google Scholar] [CrossRef]
  16. Jia, Z.; Wang, Z.; Wang, H. Characteristics of Dew Formation in the Semi-Arid Loess Plateau of Central Shaanxi Province, China. Water 2019, 11, 126. [Google Scholar] [CrossRef]
  17. Attar, N.; Khalili, K.; Behmanesh, J.; Khanmohammadi, N. On the reliability of soft computing methods in the estimation of dew point temperature: The case of arid regions of Iran. Comput. Electron. Agric. 2018, 153, 334–346. [Google Scholar] [CrossRef]
  18. Mehdizadeh, S.; Behmanesh, J.; Khalili, K. Application of gene expression programming to estimate daily dew point tempreture. Appl. Therm. Eng. 2017, 112, 1097–1107. [Google Scholar] [CrossRef]
  19. Google Earth Data. Available online: https://earth.google.com/web/@32.00526119,53.69582824,2905.05963411a,3214415.88153899d,35y,0h,0t,0r (accessed on 3 February 2018).
  20. Vapnik, V.; Golowich, S.; Smola, A. Support vector method for function approximation regression estimation, and signal processing. In Advances in Neural Information Processing Systems; MIT Press: Cambridge, MA, USA, 1996; Volume 9. [Google Scholar]
  21. Suykens, J.A.K.; Van Gestel, T.; Brabanter, J.; De Moor, B.; Vandewalle, J. Least Squares Support Vector Machines; World Scientific: Singapore, 2002. [Google Scholar]
  22. Sihag, P.; Jain, P.; Kumar, M. Modelling of impact of water quality on recharging rate of storm water filter system using various kernel function based regression. Model. Earth Syst. Environ. 2018, 4, 61–68. [Google Scholar] [CrossRef]
  23. Sette, S.; Boullart, L. Genetic programming: Principles and applications. Eng. Appl. Artif. Intell. 2001, 14, 727–736. [Google Scholar] [CrossRef]
  24. Aytek, A.; Kisi, O. A genetic programming approach to suspended sediment modeling. J. Hydrol. 2008, 351, 288–298. [Google Scholar] [CrossRef]
  25. Quinlan, J.R. Learning with continuous classes. In Proceedings of the 5th Australian Joint Conference on Artificial Intelligence, Hobart, Tasmania, 16–18 November 1992; World Scientific: Singapore, 1992; pp. 343–348. [Google Scholar]
  26. Pal, M. M5 model tree for land cover classification. Int. J. Remote Sens. 2006, 27, 825–831. [Google Scholar] [CrossRef]
  27. Najafzadeh, M.; Shiri, J.; Sadeghi, G.; Ghaemi, A. Prediction of the friction factor in pipes using model tree. ISH J. Hydraul. Eng. 2018, 24, 9–15. [Google Scholar] [CrossRef]
  28. Witten, I.H.; Frank, E. Data Mining: Practical Machine Learning Tools and Techniques with Java Implementations; Morgan Kaufmann: San Francisco, CA, USA, 2005. [Google Scholar]
  29. Alizadeh, Z.; Yazdi, J.; Kim, J.H.; Al-Shamiri, A.K. Assessment of Machine Learning Techniques for Monthly Flow Prediction. Water 2018, 10, 1676. [Google Scholar] [CrossRef]
  30. Sudheer, K.P.; Gosain, A.K.; Rangan, D.M.; Saheb, S.M. Modeling evaporation using an artificial neural network algorithm. Hydrol. Process. 2003, 16, 3189–3202. [Google Scholar] [CrossRef]
  31. Taylor, K.E. Summarizing multiple aspects of model performance in a single diagram. J. Geophys. Res. Atmos. 2001, 106, 7183–7192. [Google Scholar] [CrossRef]
  32. Choubin, B.; Moradi, E.; Golshan, M.; Adamowski, J.; Sajedi-Hosseini, F.; Mosavi, A. An Ensemble prediction of flood susceptibility using multivariate discriminant analysis, classification and regression trees, and support vector machines. Sci. Total Environ. 2019, 651, 2087–2096. [Google Scholar] [CrossRef] [PubMed]
  33. Kurup, P.U.; Dudani, N.K. Neural networks for profiling stress history of clays from PCPT data. J. Geotech. Geoenviron. Eng. 2014, 128, 569–579. [Google Scholar] [CrossRef]
  34. Deo, R.C.; Ghorbani, M.A.; Samadianfard, S.; Maraseni, T.; Bilgili, M.; Biazar, M. Multi-layer perceptron hybrid model integrated with the firefly optimizer algorithm for windspeed prediction of target site using a limited set of neighboring reference station data. Renew. Energy 2018, 116, 309–323. [Google Scholar] [CrossRef]
  35. Baghban, A.; Bahadori, M.; Rozyn, J.; Lee, M.; Abbas, A.; Bahadori, A. Estimation of air dew point temperature using computational intelligence schemes. Appl. Therm. Eng. 2016, 93, 1043–1052. [Google Scholar] [CrossRef]
Figure 1. Location of the study area [19].
Figure 1. Location of the study area [19].
Water 11 00582 g001
Figure 2. (A) Initial Space. (B) Feature space. (C) Insensitive error function.
Figure 2. (A) Initial Space. (B) Feature space. (C) Insensitive error function.
Water 11 00582 g002
Figure 3. The general form of the initial steps of gene expression programming (GEP) (Alvisi, Mascellani, Franchini, and Bardossy, 2005).
Figure 3. The general form of the initial steps of gene expression programming (GEP) (Alvisi, Mascellani, Franchini, and Bardossy, 2005).
Water 11 00582 g003
Figure 4. Observed and estimated values of dew point temperature (DPT) with the best models including (GEP-10), (M5-15), and (SVR-6).
Figure 4. Observed and estimated values of dew point temperature (DPT) with the best models including (GEP-10), (M5-15), and (SVR-6).
Water 11 00582 g004
Figure 5. The scatter plots of observed and estimated DPT values with most precise models including (GEP-10), (M5-15), and (SVR-6).
Figure 5. The scatter plots of observed and estimated DPT values with most precise models including (GEP-10), (M5-15), and (SVR-6).
Water 11 00582 g005
Figure 6. Taylor diagrams of estimated DPT values in the test period.
Figure 6. Taylor diagrams of estimated DPT values in the test period.
Water 11 00582 g006
Table 1. Parameters used in the GEP method.
Table 1. Parameters used in the GEP method.
ParameterQuantity
Functions used + ,   ,   × ,   ÷ ,     ,   ln ( x ) , e x p , r , S i n , C o s , A r c t a n
Number of chromosomes30
Number of genes3
Linking functionSum
Jump speed0
Mutation rate0.044
Inversion rate version0.1
One-point recombination rate T0.3
Two-point recombination rate two points0.3
Gene recombination rate the gene0.1
Gene transposition rate0.1
Table 2. Statistical characteristics of the meteorological data.
Table 2. Statistical characteristics of the meteorological data.
CCSkewnessStandard DeviationMaxMinMeanParameter
0.59−0.1310.2634.0−15.013.3Tavg (°C)
0.120.2417.4596.010.050.0RH (%)
0.010.134.33880.0848.3864.3Vp (kpa)
0.210.861.5713.00.003.40W (m/s)
0.23−0.713.7814.00.007.90S (h)
Table 3. Different combinations of input parameters in the estimation of dew point temperature.
Table 3. Different combinations of input parameters in the estimation of dew point temperature.
NumberInput ParametersNumberInput Parameters
1T9T, S
2RH10T, S, RH
3Vp11T, S, Vp
4W12T, S, W
5S13T, S, RH, Vp
6T, RH14T, S, RH, W
7T, Vp15T, S, RH, W, Vp
8T, W
Table 4. Evaluation of the performance of GEP, M5 model tree (M5), and support vector regression (SVR) models in the testing phase.
Table 4. Evaluation of the performance of GEP, M5 model tree (M5), and support vector regression (SVR) models in the testing phase.
ScenariosGEPM5SVR
RMSE (Degree)R2RMSE (Degree)R2RMSE (Degree)R2
13.400.7193.360.7273.370.724
26.200.0876.110.0926.150.102
35.740.2415.580.2435.590.247
45.900.1685.850.1735.980.158
55.200.4035.770.1885.760.187
61.560.9350.400.9960.440.996
73.440.7143.340.7313.330.731
83.500.7013.300.7343.300.736
93.180.7512.980.7873.000.783
100.960.9020.400.9960.460.994
113.100.7602.960.7882.990.784
123.210.7482.900.7952.910.796
132.570.8400.380.9960.540.994
141.050.9740.380.9960.470.994
152.600.8350.370.9960.550.989
Table 5. Obtained Equation from the M5 model tree for scenario No. 15.
Table 5. Obtained Equation from the M5 model tree for scenario No. 15.
Obtained Equation from the M5 Model TreeConditions of Input
T d e w RHT
T dew = 0.9196 T + 0.2162 RH 0.0032 EA + 0.0798 W + 0.0124 S 17.7731 RH ≤ 65.5T ≤ −7.95
T dew = 0.9311 T + 0.189 RH 0.0032 EA + 0.0399 W + 0.0127 S 15.7956 RH > 65.5T ≤ −7.95
T dew = 1.0015 T + 0.248 RH 0.0055 EA + 0.0959 W + 0.0316 S 17.3586 RH ≤ 61.5−7.95 < T ≤ −5.05
T dew = 0.915 T + 0.1943 RH 0.0043 EA + 0.0266 W + 0.0049 S 15.237 61.5 < RH ≤ 73.5−7.95 < T ≤ −5.05
T dew = 0.9199 T + 0.1804 RH + 0.0016 EA + 0.0293 W 0.0011 S 19.2341 RH > 73.5−7.95 < T ≤ −5.05
T dew = 0.8274 T + 0.2887 RH 0.0005 EA + 0.076 W 0.0087 S 23.8787 RH ≤ 50.5−5.05 < T ≤ −0.15
T dew = 0.9152 T + 0.2308 RH 0.0005 EA + 0.0213 W 0.0276 S 20.454 50.5 < RH ≤ 61.5−5.05 < T ≤ −0.15
T dew = 0.9186 T + 0.1929 RH + 0.0033 EA + 0.0302 W 0.0131 S 21.4432 61.5 < RH ≤ 74.5−5.05 < T ≤ −0.15
T dew = 0.9565 T + 0.1694 RH 0.001 EA + 0.0095 W 0.0093 S 15.8692 RH > 74.5−5.05 < T ≤ −0.15
T dew = 0.8804 T + 0.3409 RH + 0.0005 EA + 0.0036 W 0.0043 S 26.7889 RH ≤ 48.5−0.15 < T ≤ 5.75
T dew = 0.9219 T + 0.2566 RH + 0.0111 EA + 0.0036 W 0.030 S 31.9987 48.5 < RH ≤ 61.5-0.15 < T ≤ 5.75
T dew = 0.9029 T + 0.3429 RH + 0.0108 EA + 0.0036 W 0.0556 S 35.5285 RH ≤ 49.55.75 < T ≤ 9.95
T dew = 0.9108 T + 0.2683 RH + 0.0189 EA + 0.0036 W 0.0193 S 39.3861 49.5 < RH ≤ 61.55.75 < T ≤ 9.95
T dew = 0.8932 T + 0.2097 RH 0.0008 EA + 0.0033 W 0.0173 S 18.9323 61.5 < RH ≤ 73.5−0.15 < T ≤ 1.55
T dew = 0.917 T + 0.2395 RH + 0.0204 EA 0.0033 W + 0.0002 39.2768 61.5 < RH ≤ 65.51.55 < T ≤ 4.05
T dew = 0.9644 T + 0.206 RH 0.0003 EA + 0.0033 W + 0.0002 S 19.2912 65.5 < RH ≤ 73.51.55 < T ≤ 4.05
T dew = 0.9634 T + 0.1762 RH 0.0017 EA + 0.0033 W 0.0078 S 15.844 RH > 73.51.55 < T ≤ 4.05
T dew = 0.9456 T + 0.2034 RH + 0.0097 EA + 0.0033 W 0.0182 S 27.6514 RH > 73.54.05 < T ≤ 9.95
T dew = 0.8883 T + 0.4574 RH + 0.0347 EA + 0.0041 W 0.0601 S 60.2489 RH ≤ 39.59.95 < T ≤ 15.35
T dew = 0.8800 T + 0.3631 RH + 0.0346 EA + 0.0041 W 0.0337 S 56.861 39.5 < RH ≤ 46.59.95 < T ≤ 15.35
T dew = 0.9395 T + 0.6002 RH + 0.0095 EA + 0.0041 W 0.0105 S 43.7243 RH < 28.515.35 < T ≤ 20.85
T dew = 0.9251 T + 0.4472 RH + 0.0400 EA + 0.0041 W 0.0454 S 65.1975 28.5 < RH ≤ 35.515.35 < T ≤ 20.85
T dew = 0.9034 T + 0.3602 RH + 0.0358 EA + 0.0334 W 0.0234 S 58.3174 35.5 < RH ≤ 46.515.35 < T ≤ 20.85
T dew = 0.9004 T + 0.2864 RH + 0.0449 EA + 0.0053 W 0.0361 S 62.6688 46.5 < RH ≤ 59.5T ≤ 14.45
T dew = 0.9328 T + 0.2231 RH + 0.0151 EA + 0.0053 W 0.0233 S 33.5626 RH > 59.5T ≤ 14.45
T dew = 0.9333 T + 0.2794 RH + 0.0292 EA + 0.0649 W 0.0196 S 49.493 All values14.45 < T ≤ 20.85
T dew = 0.8695 T + 0.6839 RH + 0.0749 EA + 0.0808 W 0.0566 S 100.5357 RH ≤ 25.5T > 20.85
T dew = 0.8514 T + 0.4688 RH + 0.0167 EA + 0.1031 W 0.0543 S 44.4641 25.5 < RH ≤ 36.5T > 20.85
T dew = 0.8677 T + 0.3523 RH 0.0016 EA + 0.0998 W + 0.0019 S 25.4284 RH > 36.5T > 20.85
Table 6. Effect of removing input variables on the utilized model’s accuracy for predicting DPT values.
Table 6. Effect of removing input variables on the utilized model’s accuracy for predicting DPT values.
ModelInput ParametersGEPM5SVR
RMSE (Degree)R2RMSE (Degree)R2RMSE (Degree)R2
1All2.600.8350.370.9960.550.989
2Remove T5.430.2274.180.1734.650.169
3Remove S2.580.8472.730.7533.160.980
4Remove RH3.180.7523.720.6893.830.342
5Remove W2.580.9302.680.8432.630.863
6Remove Vp2.780.6422.910.5412.310.763

Share and Cite

MDPI and ACS Style

Qasem, S.N.; Samadianfard, S.; Sadri Nahand, H.; Mosavi, A.; Shamshirband, S.; Chau, K.-w. Estimating Daily Dew Point Temperature Using Machine Learning Algorithms. Water 2019, 11, 582. https://doi.org/10.3390/w11030582

AMA Style

Qasem SN, Samadianfard S, Sadri Nahand H, Mosavi A, Shamshirband S, Chau K-w. Estimating Daily Dew Point Temperature Using Machine Learning Algorithms. Water. 2019; 11(3):582. https://doi.org/10.3390/w11030582

Chicago/Turabian Style

Qasem, Sultan Noman, Saeed Samadianfard, Hamed Sadri Nahand, Amir Mosavi, Shahaboddin Shamshirband, and Kwok-wing Chau. 2019. "Estimating Daily Dew Point Temperature Using Machine Learning Algorithms" Water 11, no. 3: 582. https://doi.org/10.3390/w11030582

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