Modeling monthly reference evapotranspiration process in Turkey: application of machine learning methods

BAYRAM S., ÇITAKOĞLU H.

ENVIRONMENTAL MONITORING AND ASSESSMENT, vol.195, no.1, 2023 (SCI-Expanded) identifier identifier identifier

  • Publication Type: Article / Article
  • Volume: 195 Issue: 1
  • Publication Date: 2023
  • Doi Number: 10.1007/s10661-022-10662-z
  • Journal Name: ENVIRONMENTAL MONITORING AND ASSESSMENT
  • Journal Indexes: Science Citation Index Expanded (SCI-EXPANDED), Scopus, ABI/INFORM, Aqualine, Aquatic Science & Fisheries Abstracts (ASFA), BIOSIS, CAB Abstracts, Compendex, EMBASE, Environment Index, Food Science & Technology Abstracts, Geobase, Greenfile, MEDLINE, Pollution Abstracts, Public Affairs Index, Veterinary Science Database, Civil Engineering Abstracts
  • Keywords: Reference evapotranspiration, Multi-gene genetic programming, M5 model trees, K-nearest neighbor algorithm, Turkiye
  • Erciyes University Affiliated: Yes

Abstract

In this study, the predictive power of three different machine learning (ML)-based approaches, namely, multi-gene genetic programming (MGGP), M5 model trees (M5Tree), and K-nearest neighbor algorithm (KNN), for long-term monthly reference evapotranspiration (ET0) prediction were investigated. The input data consist of monthly solar radiation (R-s), maximum air temperature (T-max), and wind speed (W-s) derived from 163 meteorological stations in Turkey. Different input combinations were created and analyzed. The model's performance was evaluated using criteria such as Nash-Sutcliffe efficiency, Kling-Gupta efficiency, relative root mean squared error, mean absolute percentage error, and determination coefficient. Moreover, Taylor, radar, and boxplot diagrams were created. It was determined that the MGGP model outperformed both the M5Tree and the KNN models. The equation obtained from the MGGP model, for the best-performed combination of R-s-T-max-W-s, was presented. The best weather conditions were obtained as 0.029 to 31.814 MJ/m(2), - 5.8 to 45.7 degrees C, and 0.140 to 5.086 m/s for R-s, T-max, and W-s, respectively. It was also found that the R-s was the most potent input variable for ET0 estimation while W-s was the weakest.