Next Article in Journal
Impact of Perceived Value and Community Attachment on Smart Renovation Participation Willingness for Sustainable Development of Old Urban Communities in China
Next Article in Special Issue
Use of RPA Images in the Mapping of the Chlorophyll Index of Coffee Plants
Previous Article in Journal
Introducing Adaptive Machine Learning Technique for Solving Short-Term Hydrothermal Scheduling with Prohibited Discharge Zones
Previous Article in Special Issue
Automatic Segmentation and Classification System for Foliar Diseases in Sunflower
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Linear Regression Machine Learning Algorithms for Estimating Reference Evapotranspiration Using Limited Climate Data

1
Institutes of Green Bio Science and Technology, Seoul National University, Pyeongchang-gun 25354, Korea
2
Department of Agricultural Engineering, Institute of Agriculture and Life Science, Gyeongsang National University, Jinju-si 52828, Korea
*
Authors to whom correspondence should be addressed.
Sustainability 2022, 14(18), 11674; https://doi.org/10.3390/su141811674
Submission received: 30 August 2022 / Revised: 11 September 2022 / Accepted: 14 September 2022 / Published: 16 September 2022

Abstract

:
A linear regression machine learning model to estimate the reference evapotranspiration based on temperature data for South Korea is developed in this study. FAO56 Penman–Monteith (FAO56 P–M) reference evapotranspiration calculated with meteorological data (1981–2021) obtained from sixty-two meteorological stations nationwide is used as the label. All study datasets provide daily, monthly, or annual values based on the average temperature, daily temperature difference, and extraterrestrial radiation. Multiple linear regression (MLR) and polynomial regression (PR) are applied as machine learning algorithms, and twelve models are tested using the training data. The results of the performance evaluation of the period from 2017 to 2021 show that the polynomial regression algorithm that learns the amount of extraterrestrial radiation achieves the best performance (the minimum root-mean-square errors of 0.72 mm/day, 11.3 mm/month, and 40.5 mm/year for daily, monthly, and annual scale, respectively). Compared to temperature-based empirical equations, such as Hargreaves, Blaney–Criddle, and Thornthwaite, the model trained using the polynomial regression algorithm achieves the highest coefficient of determination and lowest error with the reference evapotranspiration of the FAO56 Penman–Monteith equation when using all meteorological data. Thus, the proposed method is more effective than the empirical equations under the condition of insufficient meteorological data when estimating reference evapotranspiration.

1. Introduction

Evapotranspiration is an important factor in the hydrological cycle, and it is defined as the sum of the total amount of evaporation released into the atmosphere as water vapor from the surface of water or soil and the amount of transpiration released into the atmosphere through leaf surfaces [1,2,3]. Evapotranspiration is considered as the basic data for the hydrological cycle and rainfall-runoff model verification and development, and for climate and meteorological prediction models. It needs to be considered for the management of agricultural water, which includes the planning and design of hydrological facilities and irrigation plans [4,5,6,7,8]. Reliable water resource management and planning can only be achieved by accurately calculating evapotranspiration because evapotranspiration is sensitive to changes in the growth stage or meteorological conditions during crop production [1,7,9].
The primary meteorological factors affecting evapotranspiration include radiation, air temperature, humidity, and wind speed. Considering the reference crop (alfalfa or grass), the maximum amount of evapotranspiration affected only by meteorological factors in a state where the soil is completely saturated with moisture, i.e., potential evapotranspiration, is referred to as reference evapotranspiration (ET0), and the actual evapotranspiration value for each crop period is determined by the product of ET0 and crop coefficient for each crop-growth period. ET0 is calculated using only meteorological data according to the Food and Agriculture Organization (FAO) guidelines, which is an accurate model for use in any region [6,10]. Following the Penman, Blaney–Criddle, and Modified Penman equations, at present, the FAO recommends the FAO56 Penman–Monteith (FAO56 P–M) equation as the standard; furthermore, it is used in legal plans, such as the long-term comprehensive plan for water resources and the rationalization plan for agricultural water use in South Korea.
Despite the accuracy of the evapotranspiration calculation, the FAO56 P–M equation is difficult to apply to areas with insufficient or missing meteorological data [8,11]. This model is difficult to apply realistically in areas with a low meteorological observation density or short observation history because the installation and maintenance of a meteorological station that collects meteorological data is expensive and laborious [12,13,14,15]. Given these limitations, the Blaney–Criddle, Hargreaves, Thornthwaite, and Priestly Taylor equations, which can help estimate reference evapotranspiration values with minimal climatic parameters, such as average and daily temperature differences, have been used, even though they are not as accurate as the Penman–Monteith model [6,16,17,18,19]. However, the problems associated with low accuracy, such as overestimation compared to that using FAO56 P–M, have been identified as a limitation [1,16,19]. Meanwhile, techniques used for calculating the amount of evapotranspiration on a spatial continuum using the vegetation index and thermal imaging are being developed given the development of remote sensing technologies, such as drones and satellite images. However, there are still limitations in the scope of application because of the discrepancies in the ground observation values, low temporal resolutions, and lack of verification in a space other than the ground observation point [20].
Machine learning technology is useful for predicting time-series data, such as evapotranspiration, without the prior knowledge of physical processes. Various studies have attempted to estimate ET0 from a small amount of meteorological data. Bellido-Jiménez et al. [13] tested a machine learning technique to estimate ET0 from the temperature data of southern Spain, which is a semi-arid region. Among the testing algorithms, such as multilayer perceptron (MLP), generalized regression neural network, extreme learning machine (ELM), support vector machine (SVM), random forest (RF), and XGBoost, MLP and ELM that only use the highest and lowest temperatures and solar radiation data demonstrated the best performance, with a root-mean-square error (RMSE) of at least 0.67 mm. Mehdizadeh et al. [14] applied neuro-fuzzy inference algorithms, such as the adaptive neuro-fuzzy inference system (ANFIS), ANFIS shuffled frog-leaping algorithm, and ANFIS invasive-weed optimization in Iran, which has a dry climate, to calculate the evapotranspiration using only temperature data and an RMSE minimum value of up to 0.83 mm/day. Dimitriadou and Nikolakopoulos [21] tested the utility of multiple linear regression (MLR) for estimating the reference evapotranspiration (ET0) in the Peloponnese, Greece. Sixteen regression models were tested that combined seven factors, including sunshine hours (N), mean temperature (Tmean), solar radiation (Rs), net radiation (Rn), wind speed (u2), vapor pressure deficit (es–ea), and altitude (Z). As a result, the MLR5 model using N, Tmean, u2, es-ea, and Rn was evaluated as the best performance (RMSE 0.28 mm/day, Adjueted R2 = 98.1%). For the semi-arid regions of Brazil with a lack of meteorological data, Reis et al. [22] calculated evapotranspiration more accurately than the Hargreaves equation using artificial neural networks, multiple linear regression (MLR), and ELM algorithms using only temperature data. Dou and Yang [23] applied ANN, SVM, ANFIS, and ELM to four different grasslands, forests, arable land, and wetlands in the United States, Germany, Belgium, and Sweden, respectively, where the performance based on the temperature, sunlight, relative humidity, and soil temperature data was the most accurate in the forests and the least accurate in the wetlands. Granata [24] applied machine learning (RF, M5P, Bagging) for estimating crop evapotranspiration in Florida, which is a humid subtropical region, and he reported the best performance using all climate factors required for using the FAO56 P–M equation, as well as a sufficient performance with an RMSE of 0.35–0.40 mm/day using only three factors, such as average temperature, solar radiation, and relative humidity. Raza et al. [25] revealed that the SVM algorithm is the best model for estimating the reference evapotranspiration when all meteorological factors used in the FAO56 P–M equation are used, despite the different climatic zones. Ferreira et al. [26] tested the ANN, RF, XGBoost, and multivariate adaptive regression splines algorithms that calculated the ET0 using only temperature and humidity data, and they attempted to interpolate the prediction results by linear regression again to overcome the disadvantages of the meteorological conditions in Brazil. Choi et al. [27] tested the neural network model and MLR algorithm using all meteorological factors necessary for calculating the FAO56 P–M equation for two stations as independent variables, which reported that the performance of the MLR model was not inferior as long as there was no derivation of negative values. These aforementioned studies suggest that the appropriate algorithm for estimating the ET0 may differ based on the spatial and climatic characteristics of the study subjects and composition of training data owing to the characteristics of machine learning technology that reproduces empirical patterns.
Suitable methods with dependable accuracy are needed to estimate the ET0 due to its importance in meteorological, hydrological, and agricultural studies [28]. Especially recently, there has been an increase in the research and development for implementing smart agricultural water management; however, the meteorological data and evapotranspiration observations that support it rely on the observation network of the Korea Meteorological Administration with a low observation density. Although it is desirable to measure meteorological data and evapotranspiration directly on the farmland subject to water management, it is unrealistic in terms of cost and management to install the facilities required to measure all meteorological items necessary for calculating the ET0 on farmland pavements or in each neighborhood. There is no choice but to use distant point data through spatial interpolation or an area-weighted method. Even if all meteorological factors cannot be observed, calculating reliable ET0 by only observing minimum factors, such as air temperature, will allow building the observation system required for smart water management.
This study aims to develop a technology that enables the accurate estimation of ET0, even under conditions where it is difficult to estimate FAO56 P–M ET0 due to sparse meteorological data. Therefore, ET0 was estimated based on a machine learning algorithm using limited data (only temperature data), and its applicability was reviewed through a comparison with empirical equations, such as those of Hargreaves [8], Blaney–Criddle [29], and Thornthwaite [30]. The machine learning algorithm used simple and easy-to-understand linear regression algorithms [31].

2. Methods

This study examined whether ETo can be satisfactorily estimated using a linear regression machine learning algorithm without using all the meteorological factors necessary for calculating the FAO56 P–M equation. As shown in Figure 1, twelve linear regression models were tested and the performance of the models was compared. A machine learning model was implemented using Python library’s Scikit-Learn, Pandas, Numpy, and Matplotlib. The regression model was used by changing the degree options to 1, 2, and 3 based on the MLR and PR models in the PolynomialFeatures function of the Scikit-Learn library. Furthermore, the SGDRegressor function of the Scikit-Learn library was used to apply the SGD method for minimizing the loss function. The maximum number of iterations was limited to 1000, and learning was performed at an adaptive learning rate; this was performed because the method using the adaptive learning rate had a greater learning efficiency than that when repeatedly using the learning rate directly designated by the user [19,31,32,33]. Each model was trained on a daily, monthly, and annual scale, respectively. The learning performance between twelve regression models was compared. In addition, the applicability of the proposed method was reviewed, compared with the ET0 experience equation calculated using only temperature data.

2.1. Linear Regression Machine Learning Models

The ET0 calculation method was tested by focusing on the most basic linear regression analysis algorithm among machine learning algorithms. MLR and polynomial regression (PR) algorithms were applied based on the composition of independent variables. Linear regression learns relatively quickly, has a high explanatory power, and has no significant difference in performance compared to other algorithms [31,34,35,36].
Linear regression is a representative regression algorithm that models a linear relationship between a dependent variable and one or more independent variables. The relationship between the independent and dependent variables is explained using the linear expression.
y   =   w 1 x 1 +   b
where y is the dependent variable (output value), x1 is the independent variable (input value), w1 is the weight for the independent variable, and b is the bias.
The MLR model is applied when one dependent variable and two or more independent variables are used as input data. This method is significant when each independent variable has a linear relationship with the dependent variable; however, the correlation between the independent variables is not high. The MLR model is expressed as
y   =   w 1 x 1   +   w 2 x 2 +   w 3 x 3 +     +   w n x n +   b
A PR model is a model that uses a linear model for learning nonlinear data; it is expressed as polynomials, such as quadratic and cubic equations, instead of the monomials of independent variables (Equation (2)). Quadratic and cubic polynomials were tested in this study.
y   =   w 2 x 1 2 +   w 1 x 1 +   b
In machine learning, linear regression is a learning model that identifies the parameter weights (w1, w2, …, wn) and deviation (b) that minimize the cost function (loss function) in a linear relationship. The cost function minimizes the loss by calculating the cost of the result predicted as an initial condition, i.e., repeating the process of calculating the loss, updating the model by modifying the parameters, and calculating the cost again, where the least-squares method to calculate the mean squared error (MSE) is often used. The least-squares method is a method of finding a parameter that minimizes the residual sum of squares between the target and predicted values, and is presented as follows:
MSE = 1 N   i   = 1 N ( y i , target   y i , prediction ) 2
where yi,target is the i-th target value (actual value), yi,prediction is the i-th predicted value, and N is the number of data.
Gradient descent is applied to find the optimal solution by gradually moving the learning rate and to find the parameter that minimizes the loss function. Gradient descent is simply used in machine learning to obtain the values of a function’s parameters that minimize a lose function as far as possible. In this study, stochastic gradient descent (SGD) was used among the gradient descent methods. SGD has received a considerable amount of attention just recently in the context of large-scale learning, and is suitable for regression problems with a high number of training samples (>10,000) [37,38]. Additionally, this method is more efficient than those using the entire batch, and it ensures more stable calculations than when using some randomized data [39].
Five-fold cross-validation was performed to use all independent variables for training and validation and to prevent overfitting or underfitting problems. The model was verified by dividing the training data into five equal parts, learning four parts (80%), and sequentially repeating the process of using the remaining one part (20%) as validation data five times.
The MSE, mean absolute error (MAE), RMSE, and coefficient of determination (R2) were used as indices to evaluate the input data composition of the machine learning model necessary for calculating ET0 and the performance according to the optimal model and for determining the best option.
MAE = 1 n   i   = 1 n | y i , target   y i , prediction |
MSE = 1 n   i   = 1 n ( y i , target   y i , prediction ) 2  
RMSE = MSE  
R 2 = SSE SST = 1 SSR SST = 1 1 n   i   = 1 n ( y i , target   y i , prediction ) 2 1 n   i   = 1 n ( y i , target y target ¯ ) 2 = 1 MSE var ( y target )
where SSE is the explained variation sum of squares, SSR is the unexplained variation sum of squares, SST is the total variation sum of squares (SSE + SSR), and n denotes the number of test data.
The coefficient of determination (R2) can be interpreted as the ratio of the variance of the predictable dependent variable in the independent variable. The coefficient of determination can take values between −∞ and 1 according to the mutual relation between the ground truth and prediction model [40]. The assumption that the R2 is greater than or equal to zero is based on the assumption that obtaining a negative R2 discards the regression calculation in use and utilizes the mean value. That is, when R2 is negative, it means an abnormal case in which the prediction performance is lower than the model predicted by the mean value. In particular, it can be negative when calculated on out-of-sample data (or for a linear regression without an intercept) [41].

2.2. Data Collection and Processing

In this study, daily data obtained from sixty-two meteorological stations were collected (Figure 2). The fifty-seven meteorological stations (excluding Jeju) had automated synoptic observation system (ASOS) data collected for more than 30 years, from 1981 to 2021. The five meteorological stations (Andong, Changwon, Taebaek, Bonghwa, and Cheorwon) collected data from 1983 to 2021, when continuous observation data existed.
The meteorological data between April and October (i.e., during the rice cultivation period) were limited to daily maximum temperature (DXT), daily minimum temperature (DNT), daily mean temperature (DMT), daily mean wind speed (DWS), daily relative humidity (DRH), daily solar radiation (DSR), or daily sunshine duration (DSD). From these collected data, 536,926 data were used as the research data by applying the average value of data collected from nearby observation stations. The daily ET0 was calculated using the FAO56 P–M equation as the label data for machine learning (Table 1). Based on the FAO56 P–M equation, the national average annual ET0 was 707.9 mm and the standard deviation was 59.6 mm. Figure 3 presents the spatiotemporal distribution of the average annual ET0 from April to October. ET0 was 711.1 mm/year from 1981 to 1990, 704.9 mm/year from 1991 to 2000, 706.3 mm/year from 2001 to 2010, and 720.0 mm/year from 2011 to 2021. In the past, ET0 was high in the southern region, but recently, it has been on the rise in the central and northern regions of Gyeonggi-do.

2.3. Learning Model and Training Dataset

The learning model was constructed by varying the conditions of the available meteorological items to search for the daily ET0 calculation model based on limited meteorological items. The daily ET0 based on the FAO56 P–M equation was set as a label, and 12 input variable datasets were created: 4 for each learning algorithm, with different combinations of temperature and radiation items (Table 2). Among the input variables, the diurnal temperature range (DTR) and daily extraterrestrial radiation (DER) were, respectively, calculated and processed from the daily maximum and minimum temperatures and the latitude and longitude coordinates of the meteorological station. Extraterrestrial radiation was included in the training data to consider the effect of solar radiation while minimizing the omission of training data because solar radiation was not observed at all the meteorological stations.
For each learning model, a total of 470,586 data from 1981 to 2016 collected from 62 meteorological stations nationwide were used as training data; the test performance was evaluated with a total of 66,340 data obtained from 2017 to 2021. For testing the learning model for monthly and annual ET0 calculations, a training dataset was added by processing the monthly and annual average values of the features in Table 2 for each observation station. In the monthly based training, a total of 17,563 training data and 2170 test data were used for 62 meteorological stations nationwide, whereas 2509 and 310 training and test data were used in the yearly based training.

3. Results and Discussion

3.1. Performance of Models in Daily ET0

The individual learning performances of the 12 models are presented in Table 3. The comparison of these performances indicates that the PR3_RR model trained as a cubic polynomial regression model using DMT, DTR, and DSR achieved the best results, with an MSE of 0.473 mm2/day2 (Table 3). For the test performance, all four performance indicators achieved the best results compared to the other models. Among the models that used only temperature data, the PR3_R and PR3_M models achieved the least errors, with RMSE = 0.888 mm/day and MAE = 0.724 mm/day. The PR algorithm showed an RMSE = 0.89 mm/day for the model that only learned temperature data without using extraterrestrial radiation data; this RMSE was less than that (RMSE = 0.94 mm/day) of the MLR algorithm. The model that learned extraterrestrial radiation together instead of using only temperature achieved a high performance in both training and test data regardless of the algorithm. For the R2 obtained with the daily ET0 using FAO56 P–M, all models learning DER showed a high explanatory power of 0.6 or more (Figure 4). Therefore, the method of using temperature and extraterrestrial radiation data together, and learning with the PR model instead of the MLR model, helped predict the daily ET0 more accurately than that when using only the temperature.
The error and regression fit of the monthly accumulated ET0 from 2017 to 2021 using the learning model were analyzed against the monthly ET0 (mm/month) calculated using FAO56 P–M for each national meteorological station (Figure 5). The PR3_RR and PR2_RR models achieved the least errors compared to the other ten models, and presented higher explanatory powers than the daily-based comparison with R2 = 0.78. As indicated by the daily-based comparison, models using extraterrestrial radiation data as learning data presented less errors than models only using temperature data (R2 was also higher than 0.75). The application of the PR algorithm was observed to be more appropriate than the MLR algorithm.
The comparison of the annual accumulated ET0 from April to October with the annual ET0 (mm/month) calculated using FAO56 P–M for each meteorological station in South Korea, as presented in Figure 6, indicated that the MLR_M model trained with only DXT and DNT with the MLR algorithm showed the best performance with an RMSE of 66.6 mm/year and an MAE of 51.5 mm/year. Unlike that of both the daily and monthly accumulated ET0, the annual accumulated ET0 presented the smallest error when using the learning result of an MLR algorithm. However, in all models, there was no correlation between the annual ET0 obtained using FAO56 P–M and an R2 of 0 or less. When R2 was less than 0, it indicates an abnormal case whose prediction performance was lower than the model predicted by the average value. Therefore, it was not appropriate to calculate the annual ET0 using daily learning results.

3.2. Performance of Models in Monthly ET0

The monthly average data of the independent variables for all twelve learning models presented in Table 3 were prepared to learn about monthly evapotranspiration for testing the application of linear regression machine learning technology for estimating monthly ET0. The monthly sum of the daily ET0 values calculated by FAO56 P–M was used as the label, and the monthly average of DXT (MDXT), monthly average of DNT (MDNT), monthly average of DTR (MDTR), and monthly average of DER (MDER) were prepared as the training data.
Although the PR3_RR model was the best in daily learning, in monthly learning, the PR3_MR model presented a slightly better performance than PR3_RR in all phases of training, validation, and testing (Table 4 and Figure 7). The test performance of the PR3_MR model was evaluated as MSE = 127.69 mm2/month2, RMSE = 11.30 mm/month, MAE = 8.75 mm/month, and R2 = 0.81. For the model that only learned temperature data by excluding extraterrestrial radiation, PR3_R using MDMT and MDTR as input data achieved the best results, with MSE = 19.28 mm2/month2 and MAE = 16.47 mm/month. However, there was a considerable difference in the performances from the models learning MDER together. Similar to daily learning, the performance was better when using PR rather than MLR, and when using extraterrestrial radiation data simultaneously with temperature data rather than only temperature data. Furthermore, R2 presented greater explanatory power for monthly learning than daily learning. Instead of accumulating the daily ET0, as presented in Figure 5, the approach of predicting the monthly ET0 was better in terms of error and coefficient of determination.

3.3. Performance of Models in Annual ET0

The learning model for estimating the annual ET0 was tested using the accumulated ET0 for 214 days from April 1 to October 31 calculated using FAO56 P–M as the label. From the daily meteorological data, the annual average of DXT (ADXT), annual average of the DNT (ADNT), annual average of the DTR (ADTR), and annual average of the DER (ADER) were prepared as the input data.
As with monthly learning, the PR3_MR model performed the best (Table 5). The test performance of the PR3_MR model was evaluated as MSE = 1639.7 mm2/year2, RMSE = 40.49 mm/year, MAE = 32.26 mm/year, and R2 = 0.49. Among the models that only learned temperature data, PR3_M using ADXT and ADNT as input data achieved the best results, with RMSE = 41.32 mm/year and MAE = 32.94 mm/year. Similar to daily and monthly learning, the performance was better when using PR rather than MLR, and when using extraterrestrial radiation data simultaneously with temperature data rather than only temperature data. However, R2 = 0.5 or less in all models, and the explanatory power was lower than that of the daily or monthly learning. Furthermore, compared with the results presented in Figure 6, the value predicted by learning the annual meteorological data was more appropriate than the value obtained by accumulating daily ET0 estimates.

3.4. Comparison with Temperature-Based Empirical Methods

Methods, such as Hargreaves, Blaney–Criddle, and Thornthwaite, which can estimate ET0 using temperature data alone, are used when there are insufficient observed meteorological data or when the reliability of the data is insufficient. These methods are recommended when calculating ET0 for at least ten days, and the accuracy of calculating ET0 decreases depending on the climate. However, these methods are still widely used in situations with limited meteorological data because they only require temperature data for the calculations [8,29]. Monthly and annual ET0 were calculated using the Hargreaves, Blaney–Criddle, and Thornthwaite equations using only temperature data obtained from 2017 to 2021 for sixty-two meteorological stations nationwide to analyze the usefulness of linear regression machine learning technology. For the Hargreaves equation, which can also calculate daily evapotranspiration, the daily ET0 was calculated and compared with the result obtained by machine learning [42].
The results of the PR3_RR and PR3_R models were compared with the Hargreaves equation’s daily evapotranspiration calculation results (Table 6). Daily ET0 according to FAO56 P–M used as the reference had maximum and minimum values of 9.80 and 0.76 mm/day, respectively; the range of the daily ET0 by Hargreaves and the PR3_R and PR3_RR learning models was narrower. PR3_RR had the maximum and minimum values of 6.73 and 1.05 mm/day, respectively, in the narrowest range. The average of the daily ET0 using the PR3_R and PR3_RR learning models was 3.37 mm/day, which was similar to the 3.38 mm/day of the FAO56 P–M. However, as shown in Figure 8, Hargreaves achieved a tendency of overestimation with 3.93 mm/day, compared to that of FAO56 P–M and PR3_RR. Comparing the error with FAO56 P–M, the Hargreaves equation presented a greater error with RMSE = 0.95 mm/day compared to PR3_R (0.89 mm/day) and PR3_RR (0.72 mm/day). Therefore, the PR learning model achieved an appropriate level of accuracy as an alternative to the Hargreaves equation when estimating the daily ET0.
The results of PR3_MR and PR3_R and those of the Hargreaves, Blaney–Criddle, and Thornthwaite equations were compared with the monthly ET0 by FAO56 P–M (Table 6). The national average of monthly evapotranspiration by FAO56 P–M was about 103.3 mm/month; the maximum value was 171.55 mm/month, and the minimum value was 42.25 mm/month. The averages of the Hargreaves, Blaney–Criddle, and Thornthwaite equations were 120.27 mm/month, 157.03 mm/month, and 105.76 mm/month, respectively. Figure 9 shows that the temperature-based empirical equation overestimates the monthly ET0 compared to the FAO56 P–M. Such results are attributed to the limitations of the empirical equations, as reported in the previous studies conducted by [1,43,44]. The PR3_R and PRPR3_MR models, to which machine learning was applied, presented similar statistics as that of FAO56 P–M with averages of 103.03 and 102.96 mm/month. For the error against FAO56 P–M, PR3_MR achieved the lowest RMSE of 11.30 mm/month and MAE of 8.75 mm/month. The error of the Blaney–Criddle equation, which was widely used in the design of agricultural reservoirs in the past when meteorological data were not sufficient, was the highest at RMSE = 57.52 mm/month and MAE = 53.77 mm/month. In addition to the Thornthwaite equation, R2 was also less than 0, which indicated no correlation at all. The Hargreaves equation showed a better error than the Blaney–Criddle and Thornthwaite equations; however, it achieved a larger error than the RMSE of 19.28 mm/month and MAE of 16.47 mm/month of PR3_R, which learned the monthly average temperature and diurnal temperature range. Thus, the machine learning method was observed to be more useful than the empirical equations when calculating the monthly ET0.
For the annual evapotranspiration values, the PR machine learning method showed a lower error against FAO56 P–M than the empirical equations. The national average of annual ET0 by FAO56 P–M was 722.94 mm/year, whereas the Hargreaves and Blaney–Criddle equations overestimated it to 841.91 and 1099.24 mm/year, respectively. For Thornthwaite, the average value was 740.34 mm/year and the maximum value was 831.50 mm/year, which was similar to the values for FAO56 P–M. The averages of PR3_MR and PR3_M, which showed the best performances in annual learning, were 723.89 and 723.55 mm/year, respectively; the maximum values were 791.01 and 792.02 mm/year, respectively, which suggest statistics closer to FAO56 P–M than those of the Hargreaves or Blaney–Criddle equations (Figure 10).
For the error with FAO56 P–M, the Hargreaves equation had a high error rate with an RMSE of 149.00 mm/year and R2 less than 0. This meant that there was no correlation, and this was an abnormal case where the prediction performance was lower than the model predicted by the mean value. The RMSE of the Blaney–Criddle equation was excessive at 378.75 mm/year. However, the Thornthwaite equation had an RMSE of 44.44 mm/year with FAO56 P–M, which was the same for PR3_M and PR3_MR, with RMSEs of 41.32 and 40.49 mm/year, respectively. The Hargreaves or Blaney–Criddle equations were not appropriate for annual ET0, and Thornthwaite showed the best error among the empirical equations. Similar to the daily and monthly results, the polynomial regression machine learning showed the closest results to the FAO56 P–M in annual reference evapotranspiration, and it was appropriate as an alternative to the temperature-based empirical equation.

4. Conclusions

A linear regression machine learning model that can estimate the reference evapotranspiration under limited conditions of meteorological data was tested, and its reliability was reviewed by comparing it to FAO56 P–M. The applicability of the proposed method was evaluated by comparing it to temperature-based empirical equations, such as Hargreaves, Blaney–Criddle, and Thornthwaite. To this end, meteorological data obtained from sixty-two meteorological stations across the country collected over 41 years (1981 to 2021) were used by applying multiple linear regression (MLR) and polynomial regression (PR) algorithms using Python-based libraries, such as Scikit-Learn. A total of twelve learning models were tested based on the composition of input datasets, such as maximum temperature, minimum temperature, average temperature, diurnal temperature range, and extraterrestrial radiation. Using four performance indicators, such as MSE, RMSE, MAE, and R2, the daily ET0, monthly accumulated ET0, and annual accumulated ET0 by FAO56 P–M from 2017 to 2021, respectively, were used as controls, and the test performance of each learning model was evaluated on a daily, monthly, and annual scale.
The cubic PR machine learning model (PR3_MR or PR3_RR) showed the best performance in the daily ET0, monthly ET0, and annual ET0 prediction tests. Compared to models that did not learn extraterrestrial radiation, models that learned the DTR or average temperature and extraterrestrial radiation together showed less errors in all tests. Furthermore, the PR algorithm showed a better performance than the MLR algorithm in all time scales (for example, the best R2 were daily 0.62, 0.69; monthly 0.77, 0.81; annual 0.40, 0.49 MLR and PR algorithms, respectively) and predicted an ET0 closer to FAO56 P–M than the temperature-based empirical equations. Therefore, if only the temperature data were measured, they can be learned using the cubic PR machine learning model along with the extraterrestrial radiation data (calculated based on temperature-observation-point-location information). This proposed method can calculate the ET0 of the smallest error compared to the FAO56 P–M, and can substitute the empirical equation.
This study proposed a linear regression learning model that has sufficient applicability compared to the empirical equations and can be applied to all of South Korea. However, inland and coastal areas are different, and there may be errors in predicting the ET0 according to the spatiotemporal deviation of learning data. Thus, it is necessary to study the optimal model for each region by examining other training data and learning models considering regional characteristics.

Author Contributions

All authors contributed to the design and development of this manuscript. Conceptualization and methodology, M.-W.J. and S.-J.K.; data collection, S.-J.B.; implementation of machine learning algorithms, S.-J.K.; writing—original draft preparation, M.-W.J.; writing—review and editing, M.-W.J. and S.-J.B.; project administration, M.-W.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Research Foundation of Korea (NRF), grant funded by the Korea government (MSIT) (No. NRF-2019R1A2C1010125).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kim, S.J.; Kim, M.I.; Lim, C.H.; Lee, W.K.; Kim, B.J. Applicability analysis of FAO56 Penman-Monteith methodology for estimating potential evapotranspiration in Andong dam watershed using limited meteorological data. J. Clim. Change Res. 2017, 8, 125–143. [Google Scholar] [CrossRef]
  2. Dimitriadou, S.; Nikolakopoulos, K.G. Artificial Neural Networks for the Prediction of the Reference Evapotranspiration of the Peloponnese Peninsula, Greece. Water 2022, 14, 2027. [Google Scholar] [CrossRef]
  3. Qu, W.; Jin, Z.; Zhang, Q.; Gao, Y.; Zhang, P.; Chen, P. Estimation of evapotranspiration in the yellow river basin from 2002 to 2020 based on GRACE and GRACE follow-on observations. Remote Sens. 2022, 14, 730. [Google Scholar] [CrossRef]
  4. Dimitriadou, S.; Nikolakopoulos, K.G. Evapotranspiration trends and interactions in light of the anthropogenic footprint and the climate crisis: A review. Hydrology 2021, 8, 163. [Google Scholar] [CrossRef]
  5. Ajjur, S.B.; Al-Ghamdi, S.G. Evapotranspiration and water availability response to climate change in the middle east and north Africa. Clim. Change 2021, 166, 28. [Google Scholar] [CrossRef]
  6. 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]
  7. Wanniarachchi, S.; Sarukkalige, R. A review on evapotranspiration estimation in agricultural water management: Past, present, and future. Hydrology 2022, 9, 123. [Google Scholar] [CrossRef]
  8. 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]
  9. Feng, K.; Tian, J. Forecasting reference evapotranspiration using data mining and limited climatic data. Eur. J. Remote Sens. 2021, 54, 363–371. [Google Scholar] [CrossRef]
  10. Berti, A.; Tardivo, G.; Chiaudani, A.; Rech, F.; Borin, M. Assessing reference evapotranspiration by the Hargreaves method in north-eastern Italy. Agric. Water Manag. 2014, 140, 20–25. [Google Scholar] [CrossRef]
  11. Seo, H.; Kim, J.; Park, H.; Kim, Y. Estimating evapotranspiration with the complementary relationship at Fluxnet sites over Asia. KSCE J. Civil. Environ. Eng. Res. 2017, 37, 303–310. [Google Scholar] [CrossRef] [Green Version]
  12. Eccel, E. Estimating air humidity from temperature and precipitation measures for modelling applications. Meteorol. Appl. 2012, 19, 118–128. [Google Scholar] [CrossRef]
  13. Bellido-Jiménez, J.A.; Estévez, J.; García-Marín, A.P. New machine learning approaches to improve reference evapotranspiration estimates using intra-daily temperature-based variables in a semi-arid region of Spain. Agric. Water Manag. 2021, 245, 106558. [Google Scholar] [CrossRef]
  14. Mehdizadeh, S.; Mohammadi, B.; Pham, Q.B.; Duan, Z. Development of boosted machine learning models for estimating daily reference evapotranspiration and comparison with empirical approaches. Water 2021, 13, 3489. [Google Scholar] [CrossRef]
  15. Dlouhá, D.; Dubovský, V.; Pospíšil, L. Optimal calibration of evaporation models against Penman–Monteith equation. Water 2021, 13, 1484. [Google Scholar] [CrossRef]
  16. Lang, D.; Zheng, J.; Shi, J.; Liao, F.; Ma, X.; Wang, W.; Chen, W.; Zhang, M. A comparative study of potential evapotranspiration estimation by eight methods with FAO Penman–Monteith method in southwestern China. Water 2017, 9, 734. [Google Scholar] [CrossRef]
  17. Elferchichi, A.; Giorgio, G.A.; Lamaddalena, N.; Ragosta, M.; Telesca, V. Variability of temperature and Its impact on reference evapotranspiration: The test case of the Apulia region (southern Italy). Sustainability 2017, 9, 2337. [Google Scholar] [CrossRef]
  18. Gentilucci, M.; Bufalini, M.; Materazzi, M.; Barbieri, M.; Aringoli, D.; Farabollini, P.; Pambianchi, G. Calculation of potential evapotranspiration and calibration of the Hargreaves equation using geostatistical methods over the last 10 years in central Italy. Geosciences 2021, 11, 348. [Google Scholar] [CrossRef]
  19. Aschale, T.M.; Sciuto, G.; Peres, D.J.; Gullotta, A.; Cancelliere, A. Evaluation of reference evapotranspiration estimation methods for the assessment of hydrological impacts of photovoltaic power plants in Mediterranean climates. Water 2022, 14, 2268. [Google Scholar] [CrossRef]
  20. Lee, S.; Kim, K.; Kim, Y.; Kim, J.; Park, S.; Yun, Y.; Kim, N.; Lee, Y. Deep learning-based estimation and mapping of evapotranspiration in cropland using local weather prediction model and satellite data. J. Kor. Cartograph. Assoc. 2018, 18, 105–116. [Google Scholar] [CrossRef]
  21. Dimitriadou, S.; Nikolakopoulos, K.G. Multiple linear regression models with limited data for the prediction of reference evapotranspiration of the Peloponnese, Greece. Hydrology 2022, 9, 124. [Google Scholar] [CrossRef]
  22. Reis, M.M.; da Silva, A.J.; Junior, J.Z.; Santos, L.D.T.; Azevedo, A.M.; Lopes, É.M.G. Empirical and learning machine approaches to estimating reference evapotranspiration based on temperature data. Comput. Electron. Agric. 2019, 165, 104937. [Google Scholar] [CrossRef]
  23. 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]
  24. Granata, F. Evapotranspiration evaluation models based on machine learning algorithms—A comparative study. Agric. Water Manag. 2019, 217, 303–315. [Google Scholar] [CrossRef]
  25. Raza, A.; Shoaib, M.; Faiz, M.A.; Baig, F.; Khan, M.M.; Ullah, M.K.; Zubair, M. Comparative assessment of reference evapotranspiration estimation using conventional method and machine learning algorithms in four climatic regions. Pure Appl. Geophys. 2020, 177, 4479–4508. [Google Scholar] [CrossRef]
  26. Ferreira, L.B.; da Cunha, F.F.; Filho, E.I.F. Exploring machine learning and multi-task learning to estimate meteorological data and reference evapotranspiration across Brazil. Agric. Water Manag. 2022, 259, 107281. [Google Scholar] [CrossRef]
  27. Choi, Y.; Kim, M.; O’Shaughnessy, S.; Jeon, J.; Kim, Y.; Song, W.J. Comparison of artificial neural network and empirical models to determine daily reference evapotranspiration. J. Kor. Soc. Agric. Eng. 2018, 60, 43–54. [Google Scholar] [CrossRef]
  28. Ahmadi, F.; Mehdizadeh, S.; Mohammadi, B.; Pham, Q.B.; Doan, T.N.C.; Vo, N.D. Application of an artificial intelligence technique enhanced with intelligent water drops for monthly reference evapotranspiration estimation. Agric. Water Manag. 2021, 244, 106622. [Google Scholar] [CrossRef]
  29. Brouwer, C.; Heibloem, M. Irrigation Water Management: Irrigation Water Needs; Training Manual; FAO: Rome, Italy, 1986; Volume 3. [Google Scholar]
  30. Thornthwaite, C.W. An approach toward a rational classification of climate. Geogr. Rev. 1948, 38, 55–94. [Google Scholar] [CrossRef]
  31. Li, Y.; Ren, X.; Zhao, F.; Yang, S. A Zeroth-order adaptive learning rate method to reduce cost of hyperparameter tuning for deep learning. Appl. Sci. 2021, 11, 10184. [Google Scholar] [CrossRef]
  32. Tang, M.; Huang, Z.; Yuan, Y.; Wang, C.; Peng, Y. A bounded scheduling method for adaptive gradient methods. Appl. Sci. 2019, 9, 3569. [Google Scholar] [CrossRef] [Green Version]
  33. Yi, D.; Ji, S.; Park, J. An adaptive optimization method based on learning rate schedule for neural networks. Appl. Sci. 2021, 11, 850. [Google Scholar] [CrossRef]
  34. Mallikarjuna, P.; Jyothy, S.A.; Sekhar Reddy, K.C. Daily reference evapotranspiration estimation using linear regression and ANN models. J. Inst. Eng. India Ser. A 2013, 93, 215–221. [Google Scholar] [CrossRef]
  35. Yirga, S.A. Modelling reference evapotranspiration for Megecha catchment by multiple linear regression. Modeling Earth Syst. Environ. 2019, 5, 471–477. [Google Scholar] [CrossRef]
  36. Mosre, J.; Suárez, F. Actual evapotranspiration estimates in arid cold regions using machine learning algorithms with in situ and Remote Sensing Data. Water 2021, 13, 870. [Google Scholar] [CrossRef]
  37. Scikit-Learn User Guide 1.1.2. Available online: https://scikit-learn.org/stable/user_guide.html (accessed on 1 April 2022).
  38. Zhou, J.; Wei, W.; Zhang, R.; Zheng, Z. Damped newton stochastic gradient descent method for neural networks training. Mathematics 2021, 9, 1533. [Google Scholar] [CrossRef]
  39. Kang, M.J. Comparison of gradient descent for deep learning. J. Korea Acad.-Ind. Coop. Soc. 2020, 21, 189–194. [Google Scholar] [CrossRef]
  40. Chicco, D.; Warrens, M.J.; Jurman, G. The coefficient of determination R-squared is more informative than SMAPE, MAE, MAPE, MSE and RMSE in regression analysis evaluation. PeerJ Comput. Sci. 2021, 7, e623. [Google Scholar] [CrossRef]
  41. Jansen, S. Machine Learning for Algorithmic Trading: Predictive Models to Extract Signals from Market and Alternative Data for Systematic Trading Strategies with Python, 2nd ed.; Packt Publishing: Birmingham, UK, 2020; ISBN 978-1839217715. [Google Scholar]
  42. Farmer, W.; Strzepek, K.; Schlosser, C.A.; Droogers, P.; Gao, X. A Method for Calculating Reference Evapotranspiration on Daily Time Scales; MIT Joint Program on the Science and Policy of Global Change; Massachusetts Institute of Technology: Cambridge, MA, USA, 2011. [Google Scholar]
  43. Lee, K.H.; Park, J.H. Calibration of the Hargreaves equation for the reference evapotranspiration estimation on a nation-wide scale. KSCE J. Civil. Environ. Eng. Res. 2008, 28, 675–681. [Google Scholar]
  44. Kim, C.G.; Lee, J.; Lee, J.E.; Kim, H. Evaluation of improvement effect on the spatial-temporal correction of several reference evapotranspiration methods. J. Korea Water Res. Assoc. 2020, 53, 701–715. [Google Scholar] [CrossRef]
Figure 1. Procedure of ET0 calculation using MLR and PR algorithms.
Figure 1. Procedure of ET0 calculation using MLR and PR algorithms.
Sustainability 14 11674 g001
Figure 2. Distribution of sixty-two meteorological stations in South Korea.
Figure 2. Distribution of sixty-two meteorological stations in South Korea.
Sustainability 14 11674 g002
Figure 3. Spatiotemporal distribution of annual ETo.
Figure 3. Spatiotemporal distribution of annual ETo.
Sustainability 14 11674 g003
Figure 4. Scatter plots of the daily ET0 calculated by FAO56 P–M and predicted by MLR and PR algorithms.
Figure 4. Scatter plots of the daily ET0 calculated by FAO56 P–M and predicted by MLR and PR algorithms.
Sustainability 14 11674 g004
Figure 5. Comparison of monthly accumulated ET0 using different models.
Figure 5. Comparison of monthly accumulated ET0 using different models.
Sustainability 14 11674 g005
Figure 6. Comparison of annual accumulated ET0 using different models.
Figure 6. Comparison of annual accumulated ET0 using different models.
Sustainability 14 11674 g006
Figure 7. Scatter plots of monthly ET0 calculated by FAO56 P–M and predicted by MLR and PR algorithms.
Figure 7. Scatter plots of monthly ET0 calculated by FAO56 P–M and predicted by MLR and PR algorithms.
Sustainability 14 11674 g007
Figure 8. Daily ET0 by Hargreaves compared with FAO56 P–M and PR3_RR.
Figure 8. Daily ET0 by Hargreaves compared with FAO56 P–M and PR3_RR.
Sustainability 14 11674 g008
Figure 9. Monthly ET0 by empirical equations compared with FAO56 P–M and PR3_MR.
Figure 9. Monthly ET0 by empirical equations compared with FAO56 P–M and PR3_MR.
Sustainability 14 11674 g009
Figure 10. Annual ET0 by empirical equations compared with FAO56 P–M and PR3_MR.
Figure 10. Annual ET0 by empirical equations compared with FAO56 P–M and PR3_MR.
Sustainability 14 11674 g010
Table 1. Nationwide average statistics of the collected climate data and calculated ET0 for the irrigation period from April to October.
Table 1. Nationwide average statistics of the collected climate data and calculated ET0 for the irrigation period from April to October.
StatisticsDaily Climate DataET0
by FAO56 P–M (mm/Year)
DXT (°C/Day)DNT (°C/Day)DRH (%/Day)DWS (m/s)DSR (MJ/m2/Day)DSD (h/Day)
Maximum41.030.9100.017.915.714.3969.0
Minimum–1.5–12.412.10.01.50.0423.7
Average24.814.572.31.97.76.4707.9
Standard deviation5.46.613.41.22.04.159.6
DXT: daily maximum temperature, DNT: daily minimum temperature, DRH: daily relative humidity, DWS: daily mean wind speed, DSR: daily solar radiation, DSD: daily sunshine duration, FAO56 P–M: FAO56 Penman–Monteith.
Table 2. Twelve models of daily-based machine learning datasets.
Table 2. Twelve models of daily-based machine learning datasets.
ModelsAlgorithmsLabelFeatures
MLR_MMultiple linear regression (MLR)Daily ET0 by FAO56 P–MDXT, DNT
MLR_MRDXT, DNT, DER
MLR_RDMT, DTR
MLR_RRDMT, DTR, DER
PR2_MQuadratic polynomial regression (PR)Daily ET0 by FAO56 P–MDXT, DNT
PR2_MRDXT, DNT, DER
PR2_RDMT, DTR
PR2_RRDMT, DTR, DER
PR3_MCubic polynomial regression (PR)Daily ET0 by FAO56 P–MDXT, DNT
PR3_MRDXT, DNT, DER
PR3_RDMT, DTR
PR3_RRDMT, DTR, DER
DXT: daily maximum temperature, DNT: daily minimum temperature, DER: daily extraterrestrial radiation, DMT: daily mean temperature, DTR: daily diurnal temperature range, and FAO56 P–M: FAO56 Penman–Monteith.
Table 3. Performance of the applied models in training and test phases for daily ET0 predictions.
Table 3. Performance of the applied models in training and test phases for daily ET0 predictions.
ModelsTrainingTest
MSE5-Fold MSEMSERMSEMAER2
MLR_M0.8020.8090.8760.9360.7750.485
MLR_MR0.5860.5950.6420.8010.6430.623
MLR_R0.8010.8090.8740.9350.7730.486
MLR_RR0.5860.5950.6410.8000.6420.623
PR2_M0.7380.7460.7960.8920.7280.532
PR2_MR0.4950.5040.5470.7400.5790.678
PR2_R0.7330.7410.7930.8910.7280.534
PR2_RR0.4770.4850.5280.7270.5610.690
PR3_M0.7320.7400.7890.8880.7240.536
PR3_MR0.4800.4880.5290.7270.5640.689
PR3_R0.7300.7380.7880.8880.7240.537
PR3_RR0.4730.4820.5210.7220.5570.694
MSE: mean squared error, 5-fold MSE: Average value of MSE with K-fold cross-validation, RMSE: root-mean-square error, MAE: mean absolute error, and R2: coefficient of determination. Bold font: best results.
Table 4. Performance of applied models in training and test phases for monthly ET0 predictions.
Table 4. Performance of applied models in training and test phases for monthly ET0 predictions.
ModelsTrainingTest
MSE5-Fold MSEMSERMSEMAER2
MLR_M333.21337.51401.3120.0317.120.40
MLR_MR127.40131.67154.0712.419.600.77
MLR_R333.20337.59400.5820.0117.100.40
MLR_RR127.40131.77153.7712.409.590.77
PR2_M321.98326.43384.6919.6116.700.43
PR2_MR109.06113.47135.8411.669.030.80
PR2_R320.28325.01379.7619.4916.640.43
PR2_RR109.04113.76136.2611.679.050.80
PR3_M320.76325.21381.8919.5416.650.43
PR3_MR100.77104.94127.6911.308.750.81
PR3_R317.44321.85371.7519.2816.470.45
PR3_RR101.84106.45128.8411.358.810.81
MSE: mean squared error, 5-fold MSE: average value of MSE with K-fold cross-validation, RMSE: root-mean-square error, MAE: mean absolute error, and R2: coefficient of determination. Bold font: best results.
Table 5. Performance of the applied models in training and test phases for annual ET0 predictions.
Table 5. Performance of the applied models in training and test phases for annual ET0 predictions.
ModelsTrainingTest
MSE5-Fold MSEMSERMSEMAER2
MLR_M1837.71907.51921.043.8334.490.40
MLR_MR1836.41918.91924.043.8634.510.40
MLR_R1838.01907.11926.143.8934.530.40
MLR_RR1836.71918.61931.143.9434.570.39
PR2_M1792.41866.01787.542.2833.650.44
PR2_MR1776.01888.01720.541.4833.070.46
PR2_R1786.71866.21758.341.9333.310.45
PR2_RR1769.11887.71697.041.2032.730.47
PR3_M1780.31853.41707.441.3232.940.46
PR3_MR1761.11887.21639.740.4932.260.49
PR3_R1781.61860.91727.441.5633.160.46
PR3_RR1761.81892.61663.940.7932.450.48
MSE: mean squared error, 5-fold MSE: average value of MSE with K-fold cross-validation, RMSE: root-mean-square error, MAE: mean absolute error, and R2: coefficient of determination. Bold font: best results.
Table 6. Comparison of the daily, monthly, and annual ET0 from 2017 to 2021 for all sixty-two meteorological stations by ET0 estimation methods.
Table 6. Comparison of the daily, monthly, and annual ET0 from 2017 to 2021 for all sixty-two meteorological stations by ET0 estimation methods.
ModelsStatisticsErrors
Max.Min.Avg.Std. devRMSEMAER2
By daily parameters
FAO56 P–M9.800.763.381.30---
Hargreaves7.870.663.931.250.950.790.47
PR3_R7.621.353.370.890.890.720.54
PR3_RR6.731.053.371.030.720.560.69
By monthly parameters
FAO56 P–M171.5542.25103.2825.89---
Hargreaves190.3455.96120.2728.4922.3919.730.25
Blaney–Criddle211.7087.67157.0330.0457.5253.77−3.94
Thornthwaite192.9322.42105.7642.8932.3226.93−0.56
PR3_R158.6168.45103.0315.8419.2816.470.45
PR3_MR149.3150.74102.9621.4311.308.750.81
By annual parameters
FAO56 P–M859.98573.33722.9456.46
Hargreaves1016.38643.33841.9171.67149.00134.99−5.96
Blaney–Criddle1162.08932.431099.2436.88378.75376.30−44.0
Thornthwaite831.50570.26740.3443.0244.4437.200.38
PR3_M792.02614.77723.5532.5441.3232.940.46
PR3_MR791.01611.96723.8932.7440.4932.260.49
Max.: maximum, Min.: minimum, Avg.: average, Std.dev: standard deviation, RMSE: root-mean-square error, MAE: mean absolute error, and R2: coefficient of determination.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kim, S.-J.; Bae, S.-J.; Jang, M.-W. Linear Regression Machine Learning Algorithms for Estimating Reference Evapotranspiration Using Limited Climate Data. Sustainability 2022, 14, 11674. https://doi.org/10.3390/su141811674

AMA Style

Kim S-J, Bae S-J, Jang M-W. Linear Regression Machine Learning Algorithms for Estimating Reference Evapotranspiration Using Limited Climate Data. Sustainability. 2022; 14(18):11674. https://doi.org/10.3390/su141811674

Chicago/Turabian Style

Kim, Soo-Jin, Seung-Jong Bae, and Min-Won Jang. 2022. "Linear Regression Machine Learning Algorithms for Estimating Reference Evapotranspiration Using Limited Climate Data" Sustainability 14, no. 18: 11674. https://doi.org/10.3390/su141811674

APA Style

Kim, S. -J., Bae, S. -J., & Jang, M. -W. (2022). Linear Regression Machine Learning Algorithms for Estimating Reference Evapotranspiration Using Limited Climate Data. Sustainability, 14(18), 11674. https://doi.org/10.3390/su141811674

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