Next Article in Journal
Variations of Stable Isotopes in Daily Precipitation in a Monsoon Region
Previous Article in Journal
Flood Frequency Estimation in Data-Sparse Wainganga Basin, India, Using Continuous Simulation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluation of Five Equations for Short-Term Reference Evapotranspiration Forecasting Using Public Temperature Forecasts for North China Plain

1
School of Water Conservation, North China University of Water Resources and Electric Power, Zhengzhou 450045, China
2
Department of Geology & Geophysics, Texas A&M University, College Station, TX 77843, USA
3
Department of Biological and Agricultural Engineering, Texas A&M University, College Station, TX 77843, USA
4
Metropolitan Solar Inc., Washington, DC 20032, USA
5
State Key Laboratory of Water Resources and Hydropower Engineering Science, Wuhan University, Wuhan 430072, China
*
Author to whom correspondence should be addressed.
Water 2022, 14(18), 2888; https://doi.org/10.3390/w14182888
Submission received: 21 June 2022 / Revised: 21 August 2022 / Accepted: 12 September 2022 / Published: 16 September 2022

Abstract

:
Accurate short-term forecasts of daily reference evapotranspiration (ET0) are essential for real-time irrigation scheduling. Many models rely on current and historical temperature data to estimate daily ET0. However, easily accessible temperature forecasts are relatively less reported in short-term ET0 forecasting. Furthermore, the accuracy of ET0 forecasting from different models varies locally and also across regions. We used five temperature-dependent models to forecast daily ET0 for a 7-day horizon in the North China Plain (NCP): the McCloud (MC), Hargreaves-Samani (HS), Blaney-Criddle (BC), Thornthwaite (TH), and reduced-set Penman–Monteith (RPM) models. Daily meteorological data collected between 1 January 2000 and 31 December 2014 at 17 weather stations in NCP to calibrate and validate the five ET0 models against the ASCE Penman–Monteith (ASCE-PM). Forecast temperatures for up to 7 d ahead for 1 January 2015–19 June 2021 were input to the five calibrated models to forecast ET0. The performance of the five models improved for forecasts at all stations after calibration. The calibrated RPM is the preferred choice for forecasting ET0 in NCP. In descending order of preference, the remaining models were ranked as HS, TH, BC, and MC. Sensitivity analysis showed that a change in maximum temperature influenced the accuracy of ET0 forecasting by the five models, especially RPM, HS, and TH, more than other variables. Meanwhile, the calibrated RPM and HS equations were better than the other models, and thus, these two equations were recommended for short-term ET0 forecasting in NCP.

1. Introduction

The North China Plain (NCP) is China’s most important grain-producing region, and it is facing severe resource constraints such as low carbon and inadequate water conservation [1]. Insufficient usable water resources and fierce competition for industrial and agricultural water have led to the severe overexploitation of groundwater in some areas [2]. Water-saving irrigation techniques and science-based irrigation decisions are important to ensure agricultural water supplies in the NCP [3,4].
Accurate daily reference evapotranspiration (ET0) forecasts are essential for real-time irrigation demand estimation [5]. The ET0 models available in the literature can be broadly classified as (1) fully physically based combination models that account for mass and energy conservation principles; (2) black-box models based on artificial neural networks, empirical relationships, and fuzzy and genetic algorithms; and (3) semi-physically based models dealing with either mass or energy conservation [6].
The fully physically based combination models mainly include the Penman–Monteith (ASCE-PM) [7], the Original Penman–Monteith [8], FAO-24 [9], Kimberly Penman [10,11,12], and FAO56 Penman–Monteith [13] models. The ASCE-PM model is recommended as the standard model for estimating ET0 and as a benchmark to evaluate other models [7,14]. However, the application of ASCE-PM is limited by the lack of required weather variables in many regions. Therefore, under the condition of the deficiency of climatological data, many studies have used machine learning models due to their excellent capability for tackling non-linear relationships between dependent and independent variables [15,16,17].
Many researchers have emphasized ET0 modeling using machine learning, such as multi-layer perceptron (MLP) [18,19], extreme learning machine (ELM) [20,21], support vector machines (SVMs) [22,23], support vector regression (SVR) [24], and extreme learning machines (ELMs). They have mainly emphasized estimating the current ET0 using various computing approaches, while less attention has been paid to machine learning codes to forecast futuristic ET0. The network was performed as a black box with no possibility of providing algebraic equations for rapid deployment. Traore et al. [19] forecasted daily ET0 based on gene-expression programming (GEP) with the limited public weather forecast in Jiangsu, China, and found that the GEP algorithm is a deployable tool for ET0 forecasts to anticipate decisions on short-term irrigation scheduling. They proposed a deployable model formulated into mathematic equations. The above methods take fewer variables as input and can produce precise results. However, for farmers and policymakers in rural areas, the machine learning methods are difficult to understand and employ. In addition, the soft computing methods are often not free and require high-performance computers. This makes it difficult to estimate ET0, particularly in remote areas where most irrigated fields are located, specifically in the developing and emerging world. For this reason, simpler empirical equations that require fewer climatic variables for ET0 calculation are utilized [25].
There are several simpler temperature-based equations such as Hargreaves, Hargreaves and Samani [26], Thornthwaite [27], and Blaney and Criddle [28], which were developed to be used with weather data. Many recent studies have investigated the feasibility of ET0 calculation based on forecast temperature. Trajkovic [29] found that a temperature-based radial basis function (RBF) network can be used to forecast ET0 when relative humidity, radiation, and wind speed data are unavailable. Lu et al. [30] compared three temperature-based estimation models with three radiation-based estimation models and recommended the Hamon model for use in the southeastern United States. The temperature-based Hargreaves–Samani (HS) model has been used to forecast ET0 in many studies and produced good results [31,32,33,34]. Other models, such as Thornthwaite (TH) and Blaney–Criddle (BC), can also be used to forecast ET0 after being calibrated [5,35,36,37].
Most of the above studies used current and historical weather data. Farmers and irrigation professionals make their decisions based on day-to-day water-balance estimation. Real-time irrigation demand estimation is challenging. Anticipating irrigation scheduling and water resources planning is a concern. Public weather forecast is easily accessible, and the forecast maximum and minimum temperatures are highly accurate [32,33,38,39]. Such temperature data could be used to accurately calculate and forecast ET0 and forecast irrigation demand in the near future [19,40]. Luo et al. [41] proposed using a locally calibrated HS model and temperature forecast for short-term ET0 forecasts and found that the method can provide daily ET0 forecasts with certain accuracy for real-time irrigation forecast. Xiong et al. [5] forecast ET0 with the BC model and temperature forecast, and the results indicate that the BC model can be an alternative and effective solution for forecast ET0 in east China. Chang et al. [35] used a modified daily TH equation and temperature forecasts for ET0 forecasts and found that TH produced good ET0 forecast results and can be a feasible solution for short-term ET0 forecasts. The temperature-based models are simple, convenient, and practical, and they are widely applicable in forecasting ET0 [38,42]. Although NCP has the greatest rate of agricultural irrigation in China, there has been little research into short-term ET0 forecasts for NCP. Some studies have limited their interest to one model or investigated models that were impractical for various reasons, so it seems that the real-time irrigation demand and irrigation water requirements forecasting in NCP may benefit from comparative research into five different temperature-based models for short-term ET0 forecasting. Most extant models require some form of correction for the particular area they are modeling, while the lack of model comparison caused confusion about whether this model is the best one for that region. It is therefore necessary to compare five models to determine the most accurate method for short-term ET0 forecasts in NCP.
The main objectives of this study were (1) to investigate the accuracy of short-time ET0 forecasting using five models with the public weather forecasts in NCP, (2) to produce five sets of calibrated parameters for temperature-based ET0 models and to recommend the best short-term ET0 model in NCP using the statistical indicators as criteria, and (3) to explore error source and the sensitivity of short-term ET0 forecasts to the percentage change in temperature forecast.

2. Materials and Methods

2.1. Study Area

The NCP is in the east of China, between 34°46′–40°25′ N and 112°30′–119°30′ E. It includes Beijing, Tianjin, the Hebei plain, the northern Henan plain, and the northern Shandong Plain north of the Yellow River [1]. The NCP is one of the three greatest plains in China. Its population is the largest in China, and it is an important grain-producing area. Agricultural irrigation is vital here. The NCP has a temperate and sub-tropical monsoon climate. It is the core area for economic development in north China, and water shortage is considered a severe threat there [2].
Figure 1 shows the locations of 17 weather stations (WS) in the NCP. Botou (BT), Tangshan (TS), and Xinxiang (XX) stations were selected as representative WS, and the characteristics of the WS are listed in Table 1. The latitudes of the WS are between 35.3° and 40.4° N, and the longitudes from 113.8° to 118.9°E, with their elevation ranging from 3.5 m to 77.3 m.

2.2. Data

Meteorological data, including daily minimum and maximum air temperature, average wind speed, average relative humidity, and sunshine duration, were obtained from 17 standard national weather stations (WS) from the China Meteorological Data Network (http://data.cma.gov.cn accessed on 30 June 2021). The 7 d ahead temperature forecast was collected from the China Weather Network (http://www.weather.com.cn accessed on 30 June 2021).
Daily meteorological data recorded from 1 January 2000 to 31 December 2014 at the 17 stations in NCP were used to calibrate and validate the five selected ET0 temperature-based forecast models against the ASCE Penman–Monteith (ASCE-PM) model. Forecast temperatures up to 7 d ahead for 1 January 2015–19 June 2021 were input to the five calibrated models to forecast ET0. Data used in the study are shown in Table 2.

2.3. Method

Five air temperature-based models, the McCloud (MC), Thornthwaite (TH), Blaney–Criddle (BC), Hargreaves–Samani (HS), and reduced-set Penman–Monteith (RPM) models, were selected for this study to calculate and forecast local ET0 [29,31,36,37,38,42,43]. Air-temperature-based models can accurately and reliably forecast daily ET0 in China [33,34,38,42]. Penman–Monteith (ASCE-PM) [7] has been recommended as the standard model for estimating ET0 and as a target to evaluate other models [7,14]. However, ASCE-PM requires several meteorological variables as input, and thus it is not practical for areas where the required daily meteorological input data is not available [5,29,36]. Figure 2 shows the workflow of the ET0 forecast in this study.

2.3.1. ASCE Penman–Monteith Model

The ASCE-Penman–Monteith (ASCE-PM) equation is [7]:
ET 0 , PM = 0.408 Δ R n G + γ Cn T + 273 U 2 e s e a Δ + γ 1 + CdU 2
where: ET 0 , PM is the daily ET0 (mm/d) calculated by the ASCE-PM equation; Rn is the net irradiance on the crop surface (MJ/m2/d); G is the soil heat flux density (MJ/m2/d); T is the air temperature at a height of 2 m (°C); U2 is the wind speed at a height of 2 m (m/s); e s is the saturated air vapor pressure (kPa); e a is the measured vapor pressure (kPa); Δ is the slope of the vapor pressure curve (kPa/°C); γ is the psychometric constant (kPa/°C).
If the wind speed at 2 m is not measured, the equation is empirically calculated as a function using measured sunshine and measured wind speed at 10 m [13] as follows:
              U 2 = U z 4.87 ln 67.8 z 5.42
R s = a s + b s n N R a
where: U 2 is the wind speed at 2 m height (m/s); U z is the measured wind speed at z m height (m/s); z is the measurement height (10 m); R s is the daily solar or shortwave irradiance (MJ/m2/d); n is the sunshine duration (h); N is the maximum possible duration of sunshine or daylight hours (h); Ra is the extraterrestrial radiation (MJ/m2/d); a s (0.25) and b s (0.5) are constants recommended in ASCE-EWRI [7]. Considering that there is short grass around the China Meteorological weather station and that ET0 was calculated on the daily time scale, the value of Cn and Cd are 900 and 0.34, respectively.

2.3.2. McCloud Model

McCloud [44] developed the McCloud (MC) equation for the ET0 estimation of land in Florida. The MC equation is:
ET 0 , MC = KW 1.8 T
where: ET 0 , MC is ET0 calculated by the MC equation (mm/d); K (0.254) and W (1.07) are both recommended default values [42]; T is the daily mean measured air temperature (°C) at a screen height of 2.0 m.

2.3.3. Hargreaves–Samani Model

The Hargreaves–Samani (HS) model [26] for ET0 estimation is expressed as:
ET 0 , HS = 1 λ C R a T max T min E T max + T min 2 + T
where: ET 0 , HS is the daily ET0 calculated by the HS equation (mm/d); T max and T min are the maximum and minimum air temperature (°C); Ra is the extraterrestrial irradiance (MJ/m2/d); C (0.0023) and E (0.5) are parameters with recommended values [26,33,34]. The HS model has been widely used in various climates, and many studies have found that the values of the parameters C and E vary from region to region. To ensure model accuracy, the local calibration of these two parameters is necessary [32,34].

2.3.4. Blaney–Criddle Model

Blaney and Criddle [28] developed the BC model to calculate daily ET0. The equation is:
ET 0 , BC = P 0.46 T m + 8.13  
where ET 0 , BC is the daily ET0 calculated by the BC model (mm/d), and Tm is the daily mean air temperature (°C) calculated according to daily T max and T min , as T m = T max + T min / 2 . The parameter p is the daily percentage of annual daytime hours (%).

2.3.5. Adjusted Thornthwaite Model

The Thornthwaite (TH) model was developed by Thornthwaite [27], and the adjusted TH model was developed by Pereira and Pruitt [45]. The adjusted TH model to calculate ET0 is:
ET 0 , TH 0   T ef < 0   16 · C · 10 · T ef I a   0   T ef 26   C · 415.85 + 32.24 T ef 0.43 T ef 2   T ef > 26  
with:
T ef = 0.5 K 3 Tmin max
I = i = 1 12 0.2 T i 1.514
a = 0.49 + 0.0179 I 0.0000771 I 2 + 0.000000675 I 3
where: ET 0 , TH is the daily ET0 calculated by the Thornthwaite equation (mm/d); I is a thermal index representing the local normal temperature regime; T i is the monthly average air temperature (°C), 0 if negative; the exponent a is a function of I; T ef is the effective daily temperature (°C); K is a calibration coefficient with a recommended value of 0.69 for calculating monthly ET0.

2.3.6. Reduced-Set Penman–Monteith Model

Allen et al. [13] proposed a modified version of the PM equation, the reduced-set Penman–Monteith (RPM) model that uses only T and U values as input data. Solar irradiance R S in the PM model and actual water vapor pressure e a are considered to be functions of temperature, and the wind speed at 2 m above the ground is assumed to be constant at 2 m/s because the most accurate RPM model is obtained by using 2 m/s wind speed as the input to the model [46]. The specific equation is:
R s = K T max T min 0.5 R a
e a = 0.6108 exp 17.27 T min T min
where: R S is the solar irradiance (MJ/m2/d); R a is the extraterrestrial irradiance (MJ/m2/d); T max and T min are the maximum and minimum air temperatures (°C); K is an adjustment coefficient: K = 0.16 for inland locations and K = 0.19 for coastal locations; e a is the actual vapor pressure (kPa). The RPM is used with these values to calculate ET0,RPM.

2.3.7. Calibration Method

In this study, calibration is the iterative process of comparing the ET0 by the five temperature-based models with ET0,P, revising the model if necessary, comparing again, until a model is accepted (validated). Although the HS equation appears to be applicable in any region or climate, its calculation accuracy actually varies according to region, and some studies have found the accuracy increased if C and E are calibrated for regional differences. The least squares method (LSM) is used to calibrate the values of the parameters C and E in the HS equation [29,31,33].
For the calibration of the other four temperature-based ET0 models (MC, BC, RPM, and TH), ET0 calculated by the ASCE-PM model (ET0,P) is considered to be the ground truth, and linear regression between ET0,P and ET0 calculated by the other models was used to correct the models for regional calibration. The linear regression model was:
ET 0 , P = bET 0 , T + a
where ET 0 , P is the ET0 calculated by ASCE-PM (mm/d); ET 0 , T is the ET0 estimated by the temperature-based models (mm/d); a and b are the correlation coefficients.

2.3.8. Statistical Analysis

We used several statistical indicators to test the correctness of the temperature-dependent models against the calculations of the ASCE-PM model. The indicators include accuracy (Acc), which is defined as the number of days in which the absolute errors of minimum and maximum temperature calculations are within ±2 °C as a percentage of the total number of samples [5,47]. The accuracy of ET0 forecasting was similarly defined as the number of days in which the absolute error of forecast ET0 was within ±1.5 mm/d as a percentage of the total number of samples [5]. Mean absolute error (MAE), root mean square error (RMSE), and the correlation coefficient (R) were also used for comparisons [5,41,42,46]. The equations for these indicators are:
Acc = N rK N fK × 100 %
MAE = i = 1 n E Ti E i n
RMSE = i = 1 n E Ti E i 2 n
R = i = 1 n E Ti E T ¯ E i E ¯ i = 1 n E Ti E T ¯ 2 · i = 1 n E i E ¯ 2
where: E Ti is the forecast temperature or the forecast ET0; E i is the measured air temperature or the calculated ET0 with the measured meteorological data based on ASCE-PM; i is the calculated sample ordinal number, I = 1, 2, 3…n; N is the total number of samples; E T ¯ is the mean value of the calculated sequence, and E ¯ is the mean value of the calculated sequence; K is T or E and is the error range for E Ti E i 2   ° C and E Ti E i 1.5   mm / d ; N rK is the number of days for which the forecasting is correct; N fK is the total number of days forecast.

3. Results and Discussion

3.1. Evaluation of Temperature Forecasts

A temperature-based ET0 model uses air temperature data as input to calculate ET0, so the accuracy of the air temperature forecast affects the accuracy of the ET0 forecast. We used Acc, MAE, RMSE, and R to evaluate the accuracy of the minimum and maximum temperature forecasts for the 17 WS in NCP. The values of the statistical indicators are shown in Figure 3. For all stations, the range of Acc for Tmin for 1 d ahead was 52–85% and for 7 d ahead was 42–60%. The corresponding ranges for R were 0.87–0.99 and 0.86–0.97, respectively. MAE and RMSE for T min were the least for the 1 d ahead and the greatest for 7 d ahead. The range of Acc for T max for 1 d ahead was 42.36–56.54% and for 7 d ahead was 38–49%. R for Tmax was >0.85 for all stations. The greatest values of RMSE and MAE for Tmax were for the 7 d ahead.
Figure 3 shows that the worst air temperature forecasts were for Lingxian WS (LX), but the forecast accuracy showed a similar pattern for the other WS. The T min   and T max forecast accuracy was similar for all WS; it decreased as the forecast period was extended. Table 3 shows the average values of the statistical indicators for all 17 WS for the seven forecast periods and the overall averages. Both Acc and R decreased and both MAE and RMSE increased as the forecast period increased. This indicates that, as expected, the forecast of T max and T min became less accurate as the forecast period increased. Average Acc for T min   was 52–72%. The average MAE for Tmax was 1.6 to 2.4 °C, and average RMSE was 2.2–3.1 °C. The average Acc for T max was 44–52%. The average MAE for Tmax was 2.5–2.9 °C, and the average RMSE was 3.3–3.8 °C. The R for Tmin was >0.96, and for T max , it was >0.94. Our results agreed with those of Yang et al. [42], who analyzed the 7 d temperature forecasts of 6 WS in China for 2012–2014 and found that R for both T max and T min was >0.92. We found that T min forecasts were more accurate than T max forecasts for NCP. This result is consistent with previous studies [5,38]. We concluded that the forecast accuracy of minimum and maximum air temperatures in public temperature forecasts were acceptable for use in the study of ET0 forecasts in NCP [34,41].

3.2. Calibration of Temperature-Based Models

Studies of various regions show that the ET0 for the same WS calculated by different methods can vary greatly [34,48].The ET0,P was used as the ground truth for the calibration of the temperature-based models. The calibrated parameters for each WS for the five models are shown in Table 4. The range of the calibrated coefficient for a model MC was 1.41 to 1.69, and the range of the calibrated coefficient b was 0.44 to 0.52. The corresponding ranges for BC were −0.69 to −0.02 and 0.69 to 0.86. The corresponding ranges for TH were 0.49 to 0.84 and 0.72 to 0.86. The corresponding ranges for RPM were −0.16 to 0.30 and 0.89 to 1.08. The parameters a and b for the calibrated RPM model were, respectively, close to zero and one, which indicates that the original RPM model was a good fit with the ASCE-PM. The recommended values of C and E for HS were 0.0023 and 0.50 [13]. We found that the calibrated parameter C was 0.001 or 0.002, agreeing well with the recommended value. The range of parameter E was 0.52 to 0.68, which is slightly greater than the recommended value of 0.50. These results show that the recommended values of C and E for HS are suitable for ET0 modeling NCP.
Figure 4 shows scatterplots for ET0 estimated by the five models and ASCE-PM for the calibration and validation periods for representative WS of BT, TS, and XX. These results show that, before regional calibration, ET0,MC was significantly greater than ET0,P and that ET0 estimated using RPM was closest to ET0,P. The ET0,HS and ET0,TH calculated before calibration were evenly distributed on both sides of y = x, which indicates that the calculated values of these two models were close to ET0,P. ET0,BC was generally greater than ET0,P. ET0 values calculated by the five calibrated models were closer to the line y = x than the uncalibrated models. Regional calibration increased the goodness-of-fit between all of the ET0 values calculated by the five models and ET0,P.
Figure 5 shows the variations in ET0 estimated by MC, RPM models, and ASCE-PM and Precipitation (Precip) for the three representative stations (for clarity, only the best model RPM and worst model MC are presented). After the calibration, although the maximum value of ET0 estimated by the five models was slightly less than the ground truth value, the estimated ET0 agreed well with ET0,P. That the maximum value of ET0,MC before calibration was much greater than ET0,P and that ET0,MC after calibration fitted ET0 P well. ET0,RPM fitted well with ET0,P both before and after calibration. However, the high value of ET0 estimated by the RPM model in the validation period is still lower than that of ET0,P, which is the main reason why the estimation accuracy of the RPM model is up to 98%. The estimation accuracy of the HS, TH, and BC models is between MC and RPM models. The low value of ET0 estimated by the BC model is closer to ET0,P than ET0,MC, making the estimation accuracy of the BC model in NCP 84% higher than that of the MC model. The comparison between the daily precipitation of representative stations with ET0,P shows that ET0,P is high on rainy days, which indicates that there is a positive correlation between precipitation and ET0,P. Liu et al. [49] compared the performance of several temperature models for north China and found that the ET0 calculated by HS well fitted FAO56-PM, which is consistent with our study. Awal et al. [50] evaluated the HS equation with FAO-ET0 at 88 stations of West Texas Mesonet and concluded that the monthly calibrated HS model yielded a performance (R2 = 0.82). Rodrigues and Braga [51] estimated ET0 during the irrigation season using nine temperature-based methods in a hot-summer Mediterranean climate. The results showed that the modified HS performed acceptably with a R2 higher than 0.78. Liu [52] evaluated six temperature-based PET methods, showing that the HS equation was the closeted to the PM calculation values (correlation coefficient = 0.85), followed by the Oudin and Hamon equations, whereas the revised Thornthwaite and Baier–Robertson equations performed relatively poorly. These agreed well with results in the study, and our studies performed better in some cases. Our results indicated that calibration increased the accuracy of all of the models’ ET0 estimation for the study area. Among them, the estimation accuracy of the modified MC model and the modified BC model is greatly improved, by 18% and 19% respectively, so that the calibrated models are suitable for practical use in the local area. The above results indicated that the five calibrated models can be used for ET0 forecasting.

3.3. ET0 Forecasts

Figure 6 shows a plot of the 1 d ahead, 4 d ahead, and 7 d ahead ET0 forecasts for three representatives WS in NCP for the period 1 January 2015 to 20 June 2021. These results show that forecast ET0 has the same trend as the measured ET0,P for all models. The minimum value of ET0,MC was greater than the minimum value of ET0,P, and there was an obvious difference between ET0,BC and ET0,P. ET0,RPM and ET0,HS were similar, and both were less than ET0,P. This is consistent with the previous findings that RPM and HS underestimate ET0,P [33,42]. Forecasts of ET0,BC were similar to those of ET0,MC, but the lower values of ET0,BC were closer to ET0,P. Predictions of ET0,TH were less accurate than the values of ET0,RPM and ET0,HS but more accurate than forecasts of ET0,MC and ET0,BC. The errors show that the forecasts of the five models differ from ET0,P, and it is evident that the accuracy of the ET0 forecasting is determined by the computational model used.
Table 5 presents the average values of the four statistical indicators of ET0 forecast accuracy for the 17 WS in the 1–7 d ahead forecasting period in NCP. The Acc and R decreased, and the MAE and RMSE increased as the forecast time increased. The Acc for MC for the 7 d ahead forecast was >75%, and R for the same forecast was greater than 0.65. The Acc for the 7 d ahead forecast was greater than 80%, and the R for the same forecast was greater than 0.79 for the other four models. The MAE and RMSE for 7 d ahead forecasts were greater for MC than for the other four models. The greatest average values of MAE (0.1 mm/d) and RMSE (0.2 mm/d) were for MC forecasts. The changes in the four statistical indicators show that the ET0 forecasts of each model become less accurate as the forecast time increases, which is consistent with the trend of temperature forecasts. However, the ET0 forecast accuracy of the five models’ maximum decreased by 3% less than the decrease in temperature forecasts.
With the above analysis, the highest average forecast accuracy of the RPM model is 85%, and the lowest average forecast accuracy of the MC model is 75%. The reason for the greater accuracy of RPM may be that the model is conceptually similar to ASCE-PM. In practice, the forecast of ET0 using RPM is fundamental to irrigation decision-making, a view that has been expressed in other studies [42]. ET0,HS was second in accuracy, and the originally uncalibrated HS model’s average forecast accuracy was 84%, which indicates that the architecture of the model is suitable for ET0 forecasts. The forecast accuracy of TH and BC ranked third and fourth. Xiong et al. [5] forecast ET0 using a monthly calibrated BC equation using temperature forecasts in east China, and the accuracy of the ET0 forecasting was acceptable (MAE varied between 0.73 and 0.82 mm/d, and R ranged from 0.74 to 0.90). According to the study of Chang et al. [35], a modified daily Thornthwaite (TH) equation can be used to forecast ET0 with temperature forecasts and found that the RMSE of forecasting ranged from 0.86 to 1.01 mm/d. The results were slightly better than this study, and it may be explained that the performance of the annually calibrated BC and TH model for ET0 forecast is further improved after monthly calibration. Yang et al. [42] studied representative meteorological stations in different climate zones in China and found that the calibrated TH equation can better forecast the daily ET0 in different climate zones in China. It agrees well with this study that BC and MC were the fourth and worst choices for daily ET0 forecasts in the NCP climate zone, compared to RPM, HS, and TH. The reason for the worse performance of BC and MC may be related to the purposes for which the models were created; the BC equation was developed to use U.S. national experimental data, and the MC equation was first used to calculate ET0 for grass and golf courses in Florida. For a while, the calibrated RPM and HS models were better than the other models, and thus these two models were recommended for short-term ET0 forecasting in NCP.
Daily ET0 prediction using the public weather forecast information can be adopted for an irrigation decision support system and real-time water allocation [41,53]. Zhang et al. [54] proposed a methodology to forecast short-term daily ETc using the ‘Kc-ETo’ approach and public weather forecasts. Considering a basic database for the short-term predicted ETc of different crops in the whole growth stages and irrigation system engineering information, the irrigation decision support system can make efficient decisions for water allocation with the help of real-time, in-field moisture detection and automatic weather forecast acquisition [55]. It uses technical methods to reduce requirements of a user’s specialized knowledge and can take a user’s managerial experience into account. It can also maximize their productivity and improve irrigation water-use efficiency [56,57].

3.4. Sensitivity and Error Analysis

The daily average minimum and maximum air temperatures for representative WS for the past 40 years were used as baseline values, and sensitivity analysis was performed by varying 5 to 20% in 5% from these base values. The amount of increase or decrease in ET0 with a 5, 10, 15, and 20% increase or decrease in maximum and minimum air temperature was calculated [35]. The percentage changes in ET0 forecasts for representative WS versus the percentage changes in maximum and minimum air temperatures are shown in Figure 7.
The percentage changes in ET0 for the five temperature-based models varied linearly with the percentage changes in Tmax and Tmin but at a greater rate, which indicates that Tmax has greater influence on ET0 forecast than Tmin. Other studies have also found that ET0 is most sensitive to maximum air temperature [35,58]. The preceding analysis of forecast temperature accuracy showed that forecast Tmin was more accurate than forecast Tmax, so the increased accuracy of the Tmax forecast will influence the accuracy of the five models.
The ET0,MC and EC0,BC decreased as Tmin decreased and increased as Tmin increased, and the change trend is consistent with that of Tmax. The ET0,HS, ET0,Th, and ET0,RPM behave oppositely, and HS, TS, and RPM were more sensitive to Tmin than BC and MC. ET0,HS, ET0,TH, and ET0,RPM changed greatly as Tmax changed. Although the change trend was the same as for ET0,BC and ET0,MC, the range of the change was greater. Figure 7 shows that the percentage change in ET0,MC varied from −4% to 5% with the percentage change in Tmin and from −7% to 9% with the percentage change in Tmax, which was similar to the behavior of ET0,BC. The percentage change in ET0,RPM varied from −9% to 8% with the percentage change in Tmin and from −29% to 32% with the percentage change in Tmax. The percentage change in ET0,HS varied from −8% to 7% with the percentage change in Tmin and from −29% to 29% with the percentage change in Tmax. The percentage change in ET0,TH varied from −4% to 5% with the percentage change in Tmin and from −25% to 29% with the percentage change in Tmax. This indicates that HS, TH, and RPM are much more sensitive to temperature than BC and MC.
Change in temperature had the greatest influence on RPM, followed in descending order by HS, TH, MC, and BC. This order corresponds almost exactly to the ET0 forecast accuracy order of the five models for the NCP. Our results are consistent with those of Yang et al. [42], who found that, when the same input data is used, differences in forecast accuracy clearly reflect the quality of the model architecture. Figure 7g shows the influence of variables on the ASCE-PM model, including Tmax, Tmin, wind speed (U2), mean relative humidity (RHmean), and sunshine hours (n). These results show Tmax has a greater impact on ET0,PM than Tmin. According to the study of G. T. Patle and D. K. Singh [58], the more sensitive the model is to the maximum temperature, the better the ET0 forecast performance would be. The ET0 forecast error factors might be divided into two categories. The errors in temperature forecasts directly influenced the ET0 forecast accuracy, and the effects of wind and relative humidity are not taken into account when forecasting ET0 [51]. With an increase in the temperature forecast accuracy, especially in maximum temperature, the model performance can be improved.

4. Conclusions

We used five temperature-based models, MC, HS, BC, TH, and RPM, to make daily forecasts of ET0 for a 1 to 7 d ahead period for 1 January 2015 to 20 June 2021 using public weather forecasts from 17 WS in the NCP. We investigated the accuracy of temperature forecasts and analyzed the sensitivity of each model to air temperature. The main conclusions are as follows.
Forecasts of Tmin were more accurate than forecasts of Tmax, and forecast temperature accuracy decreased by 14% as the forecast period increased. The accuracy of ET0 forecasts by the five models similarly decreased by 3% as the period of the forecast increased. The forecast accuracy of the five calibrated models was all greater than 75%. All five calibrated temperature-based models can be used for short-term ET0 forecasts for NCP, but RPM, HS, and TH are, in that order, more accurate than BC or MC. Therefore, the calibrated RPM and HS models were better than the other models and thus these two equations were recommended for short-term ET0 forecasting in NCP.
Sensitivity analysis showed that Tmax had greater influence on the ET0 forecast than Tmin. Changes in temperature had the greatest influence on RPM, followed in descending order by HS, TH, MC, and BC. The ET0 forecast error factors might be in two categories. The errors in temperature forecasts directly influenced the ET0 forecast accuracy, while other weather factors not considered in the temperature-based model also caused ET0 forecast errors. With an increase in the temperature forecast accuracy, especially in maximum temperature, the model performance can be improved.
The five calibrated models described herein can be used to forecast crop water and irrigation water requirements over a 7-day forecast horizon. Public weather forecast data at the provincial level can now be directly used for ET0 forecasts to be used as an irrigation scheduling decision support tool for direct field applications in crop cultivation in the NCP. Although we have taken the results provided by the ASCE-PM equation as the benchmark values, it should be noted that this is also a calculation model, which may be inadequate for some cases. It would be therefore much more reliable to attempt to improve five equations based on actual lysimeter measurements.

Author Contributions

Conceptualization, J.G.; data curation, X.Z. and J.Z.; funding acquisition, L.Z.; methodology, X.Z. and S.T.; project administration, L.Z.; supervision, G.F. and Y.L.; validation, J.Z. and S.T.; writing—original draft, X.Z.; writing—review and editing, L.Z., J.G., G.F. and Y.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financially supported by the National Natural Science Foundation of China (NSFC 51909092), the Key Science and Technology Program of Henan Province, China (222102110058), and the Key Scientific Research Projects of Higher Education Institutions in Henan Province (21A570003).

Data Availability Statement

Meteorological data recorded at 17 weather stations on the North China Plain were downloaded from the China Meteorological Data Network (http://data.cma.gov.cn accessed on 30 June 2021). The 7 d ahead temperature forecast data for the North China Plain was downloaded from the China Weather Network (http://www.weather.com.cn accessed on 30 June 2021).

Acknowledgments

We thank China Scholarship Council for providing financial support to L.Z. We thank the China Meteorological Data Sharing Service System for providing the observed meteorological data.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. He, C.; Wang, Y.Q.; Yu, W.B.; Kou, Y.H.; Yves, B.N.d.; Zhao, X.; Zhang, H.L. Comprehensive Analysis of Resource Utilization Efficiency under Different Tillage Systems in North China Plain. J. Clean. Prod. 2022, 347, 131289. [Google Scholar] [CrossRef]
  2. Yang, G.; Li, S.; Wang, H.; Wang, L. Study on Agricultural Cultivation Development Layout Based on the Matching Characteristic of Water and Land Resources in North China Plain. Agric. Water Manag. 2022, 259, 107272. [Google Scholar] [CrossRef]
  3. Koch, J.; Zhang, W.; Martinsen, G.; He, X.; Stisen, S. Estimating Net Irrigation Across the North China Plain Through Dual Modeling of Evapotranspiration. Water Resour. Res. 2020, 56, e2020WR027413. [Google Scholar] [CrossRef]
  4. Zhang, Z.; Lin, A.; Zhao, L.; Zhao, B. Attribution of Local Land Surface Temperature Variations Response to Irrigation over the North China Plain. Sci. Total Environ. 2022, 826, 154104. [Google Scholar] [CrossRef] [PubMed]
  5. Xiong, Y.; Luo, Y.; Wang, Y.; Traore, S.; Xu, J.; Jiao, X.; Fipps, G. Forecasting Daily Reference Evapotranspiration Using the Blaney-Criddle Model and Temperature Forecasts. Arch. Agron. Soil Sci. 2015, 62, 790–805. [Google Scholar] [CrossRef]
  6. Srivastava, A.; Sahoo, B.; Raghuwanshi, N.S.; Singh, R. Evaluation of Variable-Infiltration Capacity Model and MODIS-Terra Satellite-Derived Grid-Scale Evapotranspiration Estimates in a River Basin with Tropical Monsoon-Type Climatology. J. Irrig. Drain. Eng. 2017, 143, 1. [Google Scholar] [CrossRef]
  7. Allen, R.G.; Walter, I.A.; Elliot, R.L.; Howell, T.A.; Su, D.I.; Jensen, M.E.; Snyder, R.L. The ASCE Standardized Reference Evapotranspiration Equation; American Society of Civil Engineers: Reston, VA, USA, 2005. [Google Scholar] [CrossRef]
  8. Monteith, J.L. Evaporation and Environment. In Symposia of the Society for Experimental Biology; Cambridge University Press: Cambridge, UK, 1965; Volume 19, pp. 205–234. [Google Scholar]
  9. Doorenbos, J.; Pruitt, W.O. Guidelines for Prediction of Crop Water Requirements; FAO Irrigation and Drainage Paper No. 24; FAO: Rome, Italy, 1977. [Google Scholar]
  10. Wright, J.L. Derivation of alfalfa and grass reference evapotranspiration. In Evapotranspiration and Irrigation Scheduling; ASAE: San Antonio, TX, USA, 1996; pp. 133–140. [Google Scholar]
  11. Wright, J.L.; Jensen, M.E. Peak Water Requirements of Crops in Southern Idaho. J. Irrig. Drain. Div. 1972, 96, 193–201. [Google Scholar] [CrossRef]
  12. Wright, J.L. New Evapotranspiration Crop Coefficients. J. Irrig. Drain. Div. 1982, 108, 57–74. [Google Scholar] [CrossRef]
  13. 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]
  14. Allen, R.G.; Smith, M.; Pereira, L.S.; Perrier, A. An Update for the Calculation of Reference Evapotranspiration. ICID Bull. 1994, 43, 35–92. [Google Scholar]
  15. Feng, Y.; Peng, Y.; Cui, N.; Gong, D.; Zhang, K. Modeling Reference Evapotranspiration Using Extreme Learning Machine and Generalized Regression Neural Network Only with Temperature Data. Comput. Electron. Agric. 2017, 136, 71–78. [Google Scholar] [CrossRef]
  16. Landeras, G.; Ortiz-Barredo, A.; López, J.J. Forecasting Weekly Evapotranspiration with ARIMA and Artificial Neural Network Models. J. Irrig. Drain. Eng. 2009, 135, 323–334. [Google Scholar] [CrossRef]
  17. Yassin, M.A.; Alazba, A.A.; Mattar, M.A. Artificial Neural Networks versus Gene Expression Programming for Estimating Reference Evapotranspiration in Arid Climate. Agric. Water Manag. 2016, 163, 110–124. [Google Scholar] [CrossRef]
  18. Ladlani, I.; Houichi, L.; Djemili, L.; Heddam, S.; Belouz, K. Estimation of Daily Reference Evapotranspiration (ET0) in the North of Algeria Using Adaptive Neuro-Fuzzy Inference System (ANFIS) and Multiple Linear Regression (MLR) Models: A Comparative Study. Arab. J. Sci. Eng. 2014, 39, 5959–5969. [Google Scholar] [CrossRef]
  19. Traore, S.; Luo, Y.; Fipps, G. Deployment of Artificial Neural Network for Short-Term Forecasting of Evapotranspiration Using Public Weather Forecast Restricted Messages. Agric. Water Manag. 2015, 163, 363–379. [Google Scholar] [CrossRef]
  20. Abdullah, S.S.; Malek, M.A.; Abdullah, N.S.; Kisi, O.; Yap, K.S. Extreme Learning Machines: A New Approach for Prediction of Reference Evapotranspiration. J. Hydrol. 2015, 527, 184–195. [Google Scholar] [CrossRef]
  21. Gocic, M.; Petković, D.; Shamshirband, S.; Kamsin, A. Comparative Analysis of Reference Evapotranspiration Equations Modelling by Extreme Learning Machine. Comput. Electron. Agric. 2016, 127, 56–63. [Google Scholar] [CrossRef]
  22. Fan, J.; Yue, W.; Wu, L.; Zhang, F.; Cai, H.; Wang, X.; Lu, X.; Xiang, Y. Evaluation of SVM, ELM and Four Tree-Based Ensemble Models for Predicting Daily Reference Evapotranspiration Using Limited Meteorological Data in Different Climates of China. Agric. For. Meteorol. 263 2018, 225–241. [Google Scholar] [CrossRef]
  23. Ferreira, L.B.; da Cunha, F.F.; de Oliveira, R.A.; Fernandes Filho, E.I. Estimation of Reference Evapotranspiration in Brazil with Limited Meteorological Data Using ANN and SVM—A New Approach. J. Hydrol. 2019, 572, 556–570. [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. Droogers, P.; Allen, R.G. Estimating Reference Evapotranspiration under Inaccurate Data Conditions. Irrig. Drain. Syst. 2002, 16, 33–45. [Google Scholar] [CrossRef]
  26. Hargreaves, G.H.; Samani, Z.A. Reference Crop Evapotranspiration From Ambient Air Temperature. Pap.-Am. Soc. Agric. Eng. 1985, 1, 96–99. [Google Scholar]
  27. Thornthwaite, C.W. An Approach toward a Rational Classification of Climate. Geogr. Rev. 1948, 38, 55–94. [Google Scholar] [CrossRef]
  28. Blaney, H.F.; Criddle, W.D. Determining Consumptive Use and Irrigation Water Requirements; US Department of Agriculture: Washington, DC, USA, 1962.
  29. Trajkovic, S. Temperature-Based Approaches for Estimating Reference Evapotranspiration. J. Irrig. Drain. Eng. 2005, 131, 316–323. [Google Scholar] [CrossRef]
  30. Lu, J.; Sun, G.; Mcnulty, S.; Amatya, D.M. A Comparison of Six Potential Evapotranspiration Methods for Regional Use in the Southeastern United States. J. Am. Water Resour. Assoc. 2005, 41, 621–633. [Google Scholar] [CrossRef]
  31. Zhang, Q.; Duan, A.; Wang, G.; Shiva, J.K.; Shen, X.; Cai, H. Middle and Short Term Forecasting Models for Reference Evapotranspiration Based on Daily Weather Forecast. Nongye Jixie Xuebao/Trans. Chin. Soc. Agric. Mach. 2015, 46, 107–114. [Google Scholar] [CrossRef]
  32. Ren, X.; Qu, Z.; Martins, D.S.; Paredes, P.; Pereira, L.S. Daily Reference Evapotranspiration for Hyper-Arid to Moist Sub-Humid Climates in Inner Mongolia, China: I. Assessing Temperature Methods and Spatial Variability. Water Resour. Manag. 2016, 30, 3769–3791. [Google Scholar] [CrossRef]
  33. Feng, Y.; Jia, Y.; Cui, N.; Zhao, L.; Li, C.; Gong, D. Calibration of Hargreaves Model for Reference Evapotranspiration Estimation in Sichuan Basin of Southwest China. Agric. Water Manag. 2017, 181, 1–9. [Google Scholar] [CrossRef]
  34. Zhu, X.; Luo, T.; Luo, Y.; Yang, Y.; Guo, L.; Luo, H.; Fang, C.; Cui, Y. Calibration and Validation of the Hargreaves-Samani Model for Reference Evapotranspiration Estimation in China. Irrig. Drain. 2019, 68, 822–836. [Google Scholar] [CrossRef]
  35. Chang, X.; Wang, S.; Gao, Z.; Luo, Y.; Chen, H. Forecast of Daily Reference Evapotranspiration Using a Modified Daily Thornthwaite Equation and Temperature Forecasts. Irrig. Drain. 2018, 68, 297–317. [Google Scholar] [CrossRef]
  36. Chen, D.; Gao, G.; Xu, C.Y.; Guo, J.; Ren, G. Comparison of the Thornthwaite Method and Pan Data with the Standard Penman-Monteith Estimates of Reference Evapotranspiration in China. Clim. Res. 2005, 28, 123–132. [Google Scholar] [CrossRef]
  37. Rahimikhoob, A.; Hosseinzadeh, M. Assessment of Blaney-Criddle Equation for Calculating Reference Evapotranspiration with NOAA/AVHRR Data. Water Resour. Manag. 2014, 28, 3365–3375. [Google Scholar] [CrossRef]
  38. Luo, Y.; Li, S.; Peng, S.; Wang, W.; Jiao, X.; Jiang, Y.; Gu, H. Forecasting Reference Crop Evapotranspiration Based on Temperature Forecast and Hargreaves—Samani Equation. J. Drain. Irrig. Mach. Eng. 2013, 31, 987–992. [Google Scholar]
  39. Wang, W.; Xing, W.; Shao, Q. How Large Are Uncertainties in Future Projection of Reference Evapotranspiration through Different Approaches? J. Hydrol. 2015, 524, 696–700. [Google Scholar] [CrossRef]
  40. Liu, B.; Liu, M.; Cui, Y.; Shao, D.; Mao, Z.; Zhang, L.; Khan, S.; Luo, Y. Assessing Forecasting Performance of Daily Reference Evapotranspiration Using Public Weather Forecast and Numerical Weather Prediction. J. Hydrol. 2020, 590, 125547. [Google Scholar] [CrossRef]
  41. Luo, Y.; Chang, X.; Peng, S.; Khan, S.; Wang, W.; Zheng, Q.; Cai, X. Short-Term Forecasting of Daily Reference Evapotranspiration Using the Hargreaves-Samani Model and Temperature Forecasts. Agric. Water Manag. 2014, 136, 42–51. [Google Scholar] [CrossRef]
  42. Yang, Y.; Luo, Y.; Wu, C.; Zheng, H.; Zhang, L.; Cui, Y.; Sun, N.; Wang, L. Evaluation of Six Equations for Daily Reference Evapotranspiration Estimating Using Public Weather Forecast Message for Different Climate Regions across China. Agric. Water Manag. 2019, 222, 386–399. [Google Scholar] [CrossRef]
  43. Mohawesh, O.E. Spatio-Temporal Calibration of Blaney-Criddle Equation in Arid and Semiarid Environment. Water Resour. Manag. 2010, 24, 2187–2201. [Google Scholar] [CrossRef]
  44. McCloud, D.E. Water requirements of field crops in Florida as influenced by climate. Proc. Soil Sci. Soc. Fla. 1995, 15, 165–172. [Google Scholar]
  45. Pereira, A.R.; Pruitt, W.O. Adaptation of the Thornthwaite Scheme for Estimating Daily Reference Evapotranspiration. Agric. Water Manag. 2004, 66, 251–257. [Google Scholar] [CrossRef]
  46. Yang, Y.; Cui, Y.; Bai, K.; Luo, T.; Dai, J.; Wang, W.; Luo, Y. Short-Term Forecasting of Daily Reference Evapotranspiration Using the Reduced-Set Penman-Monteith Model and Public Weather Forecasts. Agric. Water Manag. 2019, 211, 70–80. [Google Scholar] [CrossRef]
  47. China Meteorological Administration (CMA). Quality Inspection of Medium-Short-Term Weather Forecast; Meteorology Press: Beijing, China, 2005. (In Chinese) [Google Scholar]
  48. Irmak, S.; Allen, R.G.; Whitty, E.B. Daily Grass and Alfalfa-Reference Evapotranspiration Estimates and Alfalfa-to-Grass Evapotranspiration Ratios in Florida. J. Irrig. Drain. Eng. 2003, 129, 360–370. [Google Scholar] [CrossRef]
  49. Liu, X.; Li, Y.; Wang, Q. Evaluation on Several Temprature-Based Methods for Estimating Reference Crop Evapotranspiration. Trans. CSAE 2006, 22, 12. [Google Scholar]
  50. Awal, R.; Habibi, H.; Fares, A.; Deb, S. Estimating Reference Crop Evapotranspiration under Limited Climate Data in West Texas. J. Hydrol. Reg. Stud. 2020, 28, 100677. [Google Scholar] [CrossRef]
  51. 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]
  52. Liu, Z. Estimating Land Evapotranspiration from Potential Evapotranspiration Constrained by Soil Water at Daily Scale. Sci. Total Environ. 2022, 834, 155327. [Google Scholar] [CrossRef]
  53. Cai, J.B.; Liu, Y.; Xu, D.; Paredes, P.; Pereira, L.S. Simulation of the Soil Water Balance of Wheat Using Daily Weather Forecast Messages to Estimate the Reference Evapotranspiration. Hydrol. Earth Syst. Sci. 2009, 13, 1045–1059. [Google Scholar] [CrossRef]
  54. Zhang, L.; Cui, Y.; Xiang, Z.; Zheng, S.; Traore, S.; Luo, Y. Short-Term Forecasting of Daily Crop Evapotranspiration Using the ‘Kc-ETo’ Approach and Public Weather Forecasts. Arch. Agron. Soil Sci. 2018, 64, 903–915. [Google Scholar] [CrossRef]
  55. Wang, W.; Cui, Y.; Luo, Y.; Li, Z.; Tan, J. Web-Based Decision Support System for Canal Irrigation Management. Comput. Electron. Agric. 2019, 161, 312–321. [Google Scholar] [CrossRef]
  56. Kim, Y.; Evans, R.G. Software Design for Wireless Sensor-Based Site-Specific Irrigation. Comput. Electron. Agric. 2009, 66, 159–165. [Google Scholar] [CrossRef]
  57. Chauhan, Y.S.; Wright, G.C.; Holzworth, D.; Rachaputi, R.C.N.; Payero, J.O. AQUAMAN: A Web-Based Decision Support System for Irrigation Scheduling in Peanuts. Irrig. Sci. 2013, 31, 271–283. [Google Scholar] [CrossRef]
  58. Patle, G.T.; Singh, D.K. Sensitivity of Annual and Seasonal Reference Crop Evapotranspiration to Principal Climatic Variables. J. Earth Syst. Sci. 2015, 124, 819–828. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Location of the 17 weather stations used in this study. The red sites Xinxiang (XX), Botou (BT), and Tangshan (TS) are representative of the north–south and east–west geographic locations.
Figure 1. Location of the 17 weather stations used in this study. The red sites Xinxiang (XX), Botou (BT), and Tangshan (TS) are representative of the north–south and east–west geographic locations.
Water 14 02888 g001
Figure 2. Workflow of the ET0 forecasting in this study.
Figure 2. Workflow of the ET0 forecasting in this study.
Water 14 02888 g002
Figure 3. Statistical indices for daily minimum temperature forecasts: (a) Accuracy, (b) MAE, (c) RMSE, and (d) R; for daily maximum temperature forecasts: (e) Accuracy, (f) MAE, (g) RMSE, and (h) R at 17 stations.
Figure 3. Statistical indices for daily minimum temperature forecasts: (a) Accuracy, (b) MAE, (c) RMSE, and (d) R; for daily maximum temperature forecasts: (e) Accuracy, (f) MAE, (g) RMSE, and (h) R at 17 stations.
Water 14 02888 g003
Figure 4. Scatterplots for ET0 calculated by the temperature-based models and the PM model for original, calibrated, and validated models in stations BT, TS, and XX; (ac) MC model, (df) HS model, (gi) BC model, (jl) TH model, (mo) RPM model; Fit-O, Fit-C, and Fit-V are respectively the best fit lines for the original, calibrated, and validated models.
Figure 4. Scatterplots for ET0 calculated by the temperature-based models and the PM model for original, calibrated, and validated models in stations BT, TS, and XX; (ac) MC model, (df) HS model, (gi) BC model, (jl) TH model, (mo) RPM model; Fit-O, Fit-C, and Fit-V are respectively the best fit lines for the original, calibrated, and validated models.
Water 14 02888 g004
Figure 5. Variation in ET0 estimated by the MC and RPM temperature-based models and PM and Precipitation (Precip) for stations BT, TS, and XX; (ac) MC model and the (df) RPM model.
Figure 5. Variation in ET0 estimated by the MC and RPM temperature-based models and PM and Precipitation (Precip) for stations BT, TS, and XX; (ac) MC model and the (df) RPM model.
Water 14 02888 g005
Figure 6. Comparisons of ET0 forecasts between MC and RPM temperature-based models and PM for BT, TS, and XX; (ac) 1 d ahead, (df) 4 d ahead, and (gi) 7 d ahead.
Figure 6. Comparisons of ET0 forecasts between MC and RPM temperature-based models and PM for BT, TS, and XX; (ac) 1 d ahead, (df) 4 d ahead, and (gi) 7 d ahead.
Water 14 02888 g006
Figure 7. Percentage change in ET0 forecasts with respect to the percentage change in temperature for stations BT, TS, and XX; Tmin (a,c,e), Tmax (b,d,f); percentage change in ET0,PM with respect to the percentage change in climatic variables (g).
Figure 7. Percentage change in ET0 forecasts with respect to the percentage change in temperature for stations BT, TS, and XX; Tmin (a,c,e), Tmax (b,d,f); percentage change in ET0,PM with respect to the percentage change in climatic variables (g).
Water 14 02888 g007
Table 1. Characteristics of weather stations in the study areas.
Table 1. Characteristics of weather stations in the study areas.
No.WMO No.StationProvinceLongitudeLatitudeElevation (m)
153,798Xingtai (XT)Hebei114.5° E37.1° N77.3
254,518Bazhou (BZ)Hebei116.4° E39.2° N8.9
354,534Tangshan (TS)Hebei118.1° E39.7° N23.2
454,539Laoting (LT)Hebei118.9° E39.4° N8.5
554,602Baoding (BG)Hebei115.5° E38.7° N16.8
654,606Raoyang (RY)Hebei115.7° E38.2° N19.0
754,618Botou (BT)Hebei116.6° E38.1° N13.2
854,624Huanghua (HH)Hebei117.3° E38.4° N4.5
954,705Nangong (NG)Hebei115.4° E37.4° N27.4
1053,986Xinxiang (XX)Henan113.9° E35.3° N73.2
1154,416Miyun (MY)Beijing116.9° E40.4° N71.8
1254,511Beijing (BJ)Beijing116.5° E39.8° N31.3
1354,525Baodi (BI)Tianjin117.3° E39.7° N5.1
1454,527Tianjin (TJ)Tianjin117.1° E39.1° N3.5
1554,715Lingxian (LX)Shandong116.6° E37.3° N18.6
1654,725Huimin (HM)Shandong117.5° E37.5° N11.7
1754,808Shenxian (SX)Shandong115.6° E36.2° N37.8
Table 2. Data used in the study.
Table 2. Data used in the study.
Data TypeDataTime PeriodData Source
Meteorological dataTmax, Tmin, RHmean, SH, U21 January 1970–20 June 2021http://data.cma.gov.cn
Temperature forecast Tmax, Tmin1 January 2015–19 June 2021http://www.weather.com.cn
Note: Tmax = daily maximum air temperature; Tmin = daily minimum air temperature; RHmean = average relative humidity; SH = sunshine duration; U2 = wind speed at 2 m.
Table 3. Average values of statistical indicators for temperature forecasts.
Table 3. Average values of statistical indicators for temperature forecasts.
Lead Time (d)TminTmax
ACC (%)MAE (°C)RMSE (°C)RACC (%)MAE (°C)RMSE (°C)R
1721.62.20.98522.53.30.95
2681.82.30.97512.63.50.95
3661.92.40.97502.63.40.95
4612.02.70.97492.73.50.95
5582.12.70.97472.83.60.95
6552.32.90.96462.83.60.94
7522.43.10.96442.93.80.94
M622.02.60.97492.73.50.95
Table 4. Calibrated parameters for the five models.
Table 4. Calibrated parameters for the five models.
StationMCHSBCTHRPM
abCEababab
XT1.450.450.0020.56−0.690.820.620.790.031.01
XX1.540.440.0020.54−0.550.790.680.770.080.97
MY1.550.480.0010.64−0.020.690.570.730.140.89
BJ1.690.490.0020.54−0.180.770.840.780.301.00
BZ1.460.490.0020.61−0.490.800.550.81−0.0021.02
BI1.550.480.0020.60−0.090.700.660.720.180.90
TJ1.580.460.0020.55−0.260.750.730.770.160.99
TS1.500.500.0020.57−0.300.760.600.800.061.02
LT1.420.510.0020.52−0.260.730.580.800.071.02
BG1.410.460.0020.63−0.550.790.560.790.011.00
RY1.470.520.0010.64−0.510.820.490.82−0.081.01
BT1.530.500.0020.64−0.520.830.580.83−0.041.05
HH1.610.490.0020.61−0.450.830.660.840.061.08
NG1.520.520.0020.64−0.600.860.490.86−0.161.07
LX1.480.510.0020.58−0.550.820.520.82−0.051.01
HM1.610.490.0010.68−0.390.800.620.820.011.03
SX1.470.460.0010.62−0.480.770.530.77−0.040.96
Table 5. Average values of the statistical indicators of five temperature-based models’ ET0 forecast accuracy for 17 stations in 1–7 d ahead forecasts for the NCP.
Table 5. Average values of the statistical indicators of five temperature-based models’ ET0 forecast accuracy for 17 stations in 1–7 d ahead forecasts for the NCP.
Statistical IndicatorLead TimeMCHSBCTHRPM
Acc (%)17587818487
27586818386
37586818386
47585818285
57585818285
67584818284
77584818184
M7585818285
MAE (mm/d)11.110.740.870.830.75
21.110.760.880.850.77
31.110.760.880.850.76
41.120.780.880.860.78
51.120.790.880.870.79
61.120.800.890.880.80
71.130.810.890.880.81
M1.120.780.880.860.78
RMSE (mm/d)11.421.011.141.101.00
21.431.031.161.121.03
31.431.031.151.121.02
41.431.051.161.131.05
51.441.061.161.151.06
61.451.081.161.161.07
71.451.101.171.171.09
M1.441.051.161.131.05
R10.670.860.800.820.86
20.660.850.800.810.85
30.660.850.800.810.85
40.660.840.800.810.84
50.660.840.790.800.84
60.650.840.790.800.83
70.650.830.790.790.83
M0.660.840.800.810.84
Note: M stands for the mean value of the statistical indicators of 1–7 d ahead ET0 forecasts by five temperature-based models for 17 stations in the NCP.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Zhang, L.; Zhao, X.; Ge, J.; Zhang, J.; Traore, S.; Fipps, G.; Luo, Y. Evaluation of Five Equations for Short-Term Reference Evapotranspiration Forecasting Using Public Temperature Forecasts for North China Plain. Water 2022, 14, 2888. https://doi.org/10.3390/w14182888

AMA Style

Zhang L, Zhao X, Ge J, Zhang J, Traore S, Fipps G, Luo Y. Evaluation of Five Equations for Short-Term Reference Evapotranspiration Forecasting Using Public Temperature Forecasts for North China Plain. Water. 2022; 14(18):2888. https://doi.org/10.3390/w14182888

Chicago/Turabian Style

Zhang, Lei, Xin Zhao, Jiankun Ge, Jiaqi Zhang, Seydou Traore, Guy Fipps, and Yufeng Luo. 2022. "Evaluation of Five Equations for Short-Term Reference Evapotranspiration Forecasting Using Public Temperature Forecasts for North China Plain" Water 14, no. 18: 2888. https://doi.org/10.3390/w14182888

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