Applicability and improvement of different evapotranspiration methods of reference crops in Jiangxi Province

To achieve accurate evaluation of evapotranspiration of reference crops (ET0) in Jiangxi, China, in the absence of systematic climatological data, with reference to the FAO-56 Penman–Monteith (P-M) equation, the Priestley-Taylor (P–T) method, the Makkink method, the Hargreaves-Samani (H–S) method, the Irmak-Allen (I-A) method, the Penman1948 (48PM) method, the Penman-Van Bavel (PVB) method, the Baier-Robertson (B-R) method, the improved Baier-Robertson (M-B-R) method, the Schendel (Sch) method, the Turc method, the Jensen-Haise (J-H) method, and the Brutsaert-Stricker (B-S) method were used to evaluate the daily climatological data collected by 26 weather stations in Jiangxi, China, and 17 weather stations in adjacent provinces. The results were compared with each other and parameter rate determination was conducted. The results indicated that the Turc method exhibited optimized applicability before parameter rate determination and the average root mean square error (RMSE) and the average normalized root mean square error (NRMSE) by this method were 0.39 mm/d and 0.157 mm, respectively. However, parameter rate determination led to negligible improvement in accuracy for this method. The Turc method could be directly applied in Jiangxi (except Nanchang). For special distribution of error after parameter rate determination, all methods exhibited significant errors in Northern Jiangxi. Herein, the 48PM method and the B-S method showed good applicability after parameter rate determination and RMSE and NRMSE of data by these methods ranged in 0.06 ~ 0.34 mm/d and 0.08 ~ 0.27, 8 ~ 27%, respectively, and their d-indices were close to 1. The annual over-estimations in weather stations in Jiangxi were below 30 mm. In the absence of data about relative humidity and wind speed, the P–T method was an appropriate simplified method for Jiangxi. In this case, α was slightly lower than the default value (1.05 ~ 1.18), RMSE was within 0.21 ~ 0.66 mm/d, and NRMSE was within 0.08 ~ 0.308 ~ 30%. Accuracy of RMSE, d-index, and NRMSE of data by the P–T method, the I-A method, and the PVB method was consistent with all stations, while that by the Mak method was slightly lower, which could be attributed to severe over-estimation in July and August. RMSE of the H–S method, the B-R method, the M-B-R method, the J-H method, and the Sch method were above 0.75 mm/d and these methods were not suitable for accurate evaluation of ET0 in Jiangxi, China. The annual ET0 was calculated by various methods (except the 48PM method and the B-S method) exhibited significant variation around 2003. This may be attributed to significant changes in certain meteorological factors over recent years. Evapotranspiration of reference crops (ET0) in Jiangxi was evaluated by 26 weather stations in Jiangxi, China, and 17 weather stations in adjacent provinces. The Turc method exhibited optimized applicability before parameter rate determination, and could be directly applied in Jiangxi (except Nanchang). In the absence of data about relative humidity and wind speed, the P–T method was an appropriate simplified method for Jiangxi. RMSE of the H–S method, the B-R method, the M-B-R method, the J-H method, and the Sch method were above 0.75 mm/d and these methods were not suitable for accurate evaluation of ET0 in Jiangxi, China. Evapotranspiration of reference crops (ET0) in Jiangxi was evaluated by 26 weather stations in Jiangxi, China, and 17 weather stations in adjacent provinces. The Turc method exhibited optimized applicability before parameter rate determination, and could be directly applied in Jiangxi (except Nanchang). In the absence of data about relative humidity and wind speed, the P–T method was an appropriate simplified method for Jiangxi. RMSE of the H–S method, the B-R method, the M-B-R method, the J-H method, and the Sch method were above 0.75 mm/d and these methods were not suitable for accurate evaluation of ET0 in Jiangxi, China.


Introduction
As basic data for determination of water demands of crops, the reference crop evapotranspiration is a key factor affecting the estimation accuracy of water demands of crops (Liu et al. 2006), as well as one of the key factors in irrigation system planning and hydrological model establishment (Allen et al. 1998;Wang et al. 2011). Owing to its wide applicability globally, the FAO Penman-Monteith (P-M) method was set by the United Nations Food Agriculture Organization (FAO) as the standard method for evaluation of ET 0 (Allen et al. 1998). However, the P-M method has a high requirement on the integrity of climatological data. Therefore, various simplified methods, including the temperature-based Hargreaves-Samani (H-S) method (Hargreaves & Allen 2003), the radiation-based IrMak method (Irmak, Irmak, Allen, & Jones), the Makkink method (Makkink 1957), and the Priestley-Taylor (P-T) method (Du et al. 2013), have been widely used. Additionally, the FAO 24 Penman radiation method (Zhang et al. 2015) and the Penman-Van Bavel (PVB) method (Yuan et al. 2014) have also attracted great attentions. Nevertheless, applications of these methods are limited by local climate and geographic conditions and many models cannot be directly applied. Many researchers have examined the applicability of different evapotranspiration models in different locations. For example, Djaman et al. (2015) assessed 16 ET0 models in the Senegal River Valley, and Muniandy et al. (2016) tested 26 ET0 models in Kluang, Malaysia, indicating that the temperature-based models tend to overestimate ET values. The Penman and Rohwer models perform better than the Penman-Monteith model. Akhavan et al. (2018) assessed and ranked different evapotranspiration models in daily corn evapotranspiration; the Hargreaves-M3 model (R MSE = 1.89 mm/day) in the temperature-based models, the Caprio (1974) model (R MSE = 1.99 mm/day) in the radiation-based models, and the Albrecht (1950) model (R MSE = 4.33 mm/day) in the mass transfer-based models were ranked first place. Paredes et al. (2018) assessed the PMT method (ETo PMT) and the HS equation (ETo HS), and the result showed that the PMT approach is more accurate than HS in the humid environments of Azores islands. Li et al. (2016) compared the applicability of six simplified methods, including the Penman (48PM) method, the PVB method, the H-S method, the P-T method, the I-A method, and the Mak method, in different climatic regions of Sichuan, China, suggesting that the H-S method should be used in Eastern Sichuan and the P-T method should be used in other regions. Zhang et al. (2015) evaluated the applicability of eight methods (including the Irmak method, the H-S method, the P-M method) in Xinxiang, Henan, and recommended the corrected Irmak method for this region. Wu et al. (2016) evaluated the applicability of various methods (including the H-S method) in Northwestern China and proposed alternative methods in the absence of climatological data. Nonetheless, reference evapotranspiration estimation is valuable when it is used in calculating actual evapotranspiration.
Jiangxi Province located in Southeastern China is in a central subtropical humid monsoon climatic region and is famous as one of the key rice producing areas in China owing to its rich hydrothermal resources. However, water resources and thermal resources may appear in different periods, resulting in droughts and floods. Indeed, summer droughts and/or floods are frequently observed, especially in July and August, resulting in severe threats to high and stable yields of rice. Additionally, Jiangxi shows varying geographical conditions and climates, but few weather stations (< 30) can provide long time meteorological data. Therefore, appropriate calculation methods for evapotranspiration of reference crops should be determined for the design of irrigation and drainage system. The farmland in this region is mostly managed by farmers, who usually lack higher education and have difficulty mastering the latest technology. Therefore, it is meaningful to establish a simpler ET0 model than any advanced methods. However, there is a lack of systematic comparison and feasibility evaluation of simple ET0 models for this region. For investigating the applicability of different calculation methods in Jiangxi and parameter localization, and eventually selection of optimized method, ET 0 has been simulated in daily scale by the 48PM method, the H-S method, the P-T method, the I-A method, the Mak method, the PVB method, the Baier-Robertson (B-R) method, the improved Baier-Robertson (M-B-R) method, the Schendel (Sch) method, the Turc method, the Jensen-Haise (J-H) method, and the Brutsaert-Stricker (B-S) method. The applicability of these methods was evaluated using the FAO56 P-M method as the reference and parameters of these methods were predetermined to select appropriate calculation method for evapotranspiration of reference crops in Jiangxi.

Background of testing region
Located at the south of the middle and lower reaches of Yangtse River, Jiangxi (24° 29′ N-30° 04′ N, 113° 34′ E-118° 28′ E) consists of Boyang Lake Plain (the north part), Jiangnan Hill (the south part), and mountains (the east and west part), as shown in Fig. 1. As a typical central subtropical humid monsoon climatic region, Jiangxi has average annual temperature of 16.3-19.5 °C (increased from north to south) and average annual precipitation of 1600 mm.

Data source
The daily climatological data (including hours of sunshine, average temperature, maximum temperature, minimum temperature, relative humidity, and wind speed at 2 m) collected from 26 weather stations in Jiangxi and 17 weather stations in provinces adjacent to Jiangxi (Fig. 1) in this study were obtained from National Meteorological Science data sharing Center of China and the monthly and annual ET0 were calculated based on the daily data. In order to compare the accuracy of different ET0 estimation methods, the quality control of meteorological data was carried out: when one or more of the daily meteorological data was missing, all the data of that day were deleted.

Methods
Twelve ET0 estimation methods were used to calculate reference crop evapotranspiration, and compared with ET0 estimation methods in FAO56P-M (Table 1).

Statistical evaluation and data processing
Four evaluation indices, including the root mean square error (R MSE ), consistency coefficient (d-index), normalized root mean square error (N RMSE ), and average deviation (M BE ), commonly used for evaluation of reference crop evapotranspiration were used and the equations are as follows: where n refers to the sample size, C refers to the observed values, P refers to the simulated values, and C refers to the average observed value.
Parameter redetermination of different reference crop evapotranspiration calculation methods was achieved using the lsqcurvefit function in MATLAB 2014a. To simplify the models, the parameters of all models were redetermined twice. After the first round, the average of parameters collected by different stations was set as the value of this specific if the coefficient of variation (CV) was below 0.1. If not, the parameters were redetermined again. Figure 2 compares the daily scale statistical results of the 12 evapotranspiration calculation methods and the P-M method under default conditions. As observed, applicability of these methods varied significantly. The Turc method exhibited best applicability, followed by the P-T method. Indeed, R MSE of data by the Turc method and the P-T method were significantly lower than those of other methods, while d-index of data of the two methods were significantly higher than those of other methods. The medians of N RMSE of data obtained by the Turc method and the P-T method were 0.16 and 0.20, respectively. According to MBE, deviations of under-estimation and over-estimation were relatively small for both methods, indicating better performance of the Turc method. Nevertheless, outliers in box plots demonstrated that accuracy of data collected by certain stations by the Turc method was lower than that by the P-T method. The average R MSE of data by other 10 methods exceeded 0.75 mm/d and these methods were not applicable in Jiangxi. Especially, the 48PM method and the B-S method had identical requirements on the integrity of climatological data as the P-M method, but their accuracies were poor (the accuracy of the 48PM method is slightly higher than that of the B-S method). For the 48PM method, median of R MSE was approximately 2 mm/d, d-index was below 0.9, N RMSE was significantly higher than 0.4, and MBE was over-estimated by 1.5 mm/d. For the B-S method, most R MSE were above 6 mm/d, d-index was below 0.6, N RMSE was significantly higher than 0.4, and MBE was over-estimated by 4 mm/d. In summary, the methods other than the Turc method and the P-T method could not be directly applied in Jiangxi before localized redetermination. ET 0 refers to the reference crop evapotranspiration (mm/d), R n refers to the net radiation (MJ/(m 2 ·d)), R s refers to the overall surface radiation(MJ/(m 2 ·d)), G refers to the soil heat flux density (MJ/(m 2 ·d)), T refers to the average temperature at 2 m (℃), u 2 refers to the wind speed at 2 m (m/s), e s refers to the saturation vapor pressure (kPa), e a refers to the actual vapor pressure (kPa), Δ refers to the slope of temperature-vapor pressure curve (kPa/℃), γ refers to the hydrometer constant (kPa/℃), λ refers to the latent heat constant of vaporization(2.45 MJ·kg −1 ), T max refers to the daily high temperature (℃), T min refers to the daily low temperature (℃), T mean refers to the daily average temperature (℃), and RH refers to the relative humidity (%)

Statistical results of different evapotranspiration methods
Penman1948 (48PM) [12] ET 0 = 1 Δ Δ+ R n − G + a Δ+ 1 + bu 2 e s − e a a = 6.43, b = 0.536 Irmak-Allen (I-A) [6] ET 0 = a + bR n + cT mean Jensen-Haise (J-H) [15] ET Figure 2 compares the redetermined daily scale statistical results of 12 evapotranspiration calculation methods and the P-M method. As observed, the accuracy of all evapotranspiration calculation methods was significantly improved after redetermination and the amplitudes varied significantly. Herein, the B-S method showed the greatest improvement (minimum R MSE , all medians < 0.1 mm/d, d-index ≈ 1, N RMSE < 10%) and optimized applicability after model redetermination (close to that of the P-M method). Additionally, R MSE of data by the 48PM method was low and the N RMSE was below 10% (except Nanchang, which may be attributed to over-estimation, according to MBE). The R MSE , d-index, and N RMSE of data by the P-T method, the I-A method, and the PVB method were close to each other. The accuracy of the P-T method, the I-A method, and the PVB method was lower than that of the B-S method and the 48PM method while higher than that of other methods. Herein, under-estimation of the P-T method was 0 ~ 1.4 mm/d, that of the PVB method was below 0.3 mm/d, and that of the I-A method was negligible. In terms of the climatological data required, the PVB method has the highest requirement on the integrity of data set: the data of relative humidity and wind speed (not required in the P-T method and the I-A method) were required for calculation of actual vapor pressures and canopy resistances. However, the PVB method exhibited no significant advantages in accuracy over the P-T method and the I-A method. Therefore, the PVB method was not applicable in Jiangxi.

Redetermined statistical results of different evapotranspiration calculation methods
The Mak method shared a similar structure with the P-T method. R MSE and N RMSE of the Mak method was slightly higher than that of the P-T method, the I-A method, and the PVB method; d-index of the Mak method was slightly lower than that of the P-T method, the I-A method, and the PVB method, while MBE of the Mak method indicated neither over-estimation nor under-estimation. In summary, the Mak method was slightly inferior to the P-T method in terms of applicability. Considering the fact that the requirements on climatological data by the Mak method and the P-T method were consistent, the P-T method was preferred in this area.
The accuracy of the Turc method and the Mak method was consistent, while parameter redetermination brought negligible improvements on the accuracy of the Turc method. For the Turc method, the defaults of R MSE and Fig. 2 Statistic results of difference reference crop evapotranspiration calculation methods under default parameters N RMSE were 0.389 mm/d and 0.157 mm, respectively, while the redetermined R MSE and N RMSE were 0.375 mm/d and 0.151 mm, respectively. In other words, the Turc method could be directly applied in Jiangxi (except Nanchang, where R MSE = 0.972 mm/d) without parameter redetermination. The average R MSE of data by the H-S method before and after parameter redetermination were 0.75 mm/d and 0.70 mm/d, respectively, indicating negligible accuracy improvement by parameter redetermination.
The average R MSE of data by the B-R method and the M-B-R method before parameter redetermination were 1.49 mm/d and 0.91 mm/d, respectively, while the ratedetermined average R MSE of data by the B-R method and the M-B-R method were 0.85 mm/d and 0.55 mm/d, respectively. This demonstrated significant accuracy improvement by parameter redetermination for the B-R method and the M-B-R method. Herein, the B-R method required temperature data but its accuracy was lower than that of the H-S method, indicating that the B-R method was not applicable in Jiangxi. In the presence of data about relative humidity, the M-B-R method showed better accuracy compared with the H-S method. However, the M-B-R method required five empirical parameters and this may be a major factor limiting the model application. The accuracy of the Sch method was still relatively low after parameter redetermination: the average R MSE was approximately 0.8 mm/d, the d-index ranged in 0.85 ~ 0.95, and the average N RMSE was above 0.3. Hence, the Sch method was not applicable in Jiangxi.
The J-H method requires sunshine data and this increases the model complexity. However, R MSE of the J-H method was increased by 0.3 mm/d compared with the I-A model, which required similar data with the J-H method, and accuracy improvement over the H-S model, which requires temperature data only, was limited. Therefore, the J-H method showed poor applicability and was not applicable in Jiangxi.

Redetermined parameters of different evapotranspiration calculation methods
The parameters were determined based on the results by the P-M method, and the linear regressions of different evapotranspiration calculation methods to the P-M method were established using the least square method. To simplify calculation methods, the secondary regression was applied for some methods. Specifically, if CV is relatively small after primary regression, average of values collected by different stations is defined as the value of this specific parameter and other parameters are using the least square method. Table 2 summarizes redetermined parameters of different evapotranspiration calculation methods. The redetermined parameter a of the 48PM method ranged in 0.81 ~ 1.21 in Jiangxi and was lower than the default value (6.43), demonstrating that the contribution to the reference crop evapotranspiration by the aerodynamic term would be over-estimated when this method was applied in Jiangxi. Indeed, the result was over-estimated by 1 ~ 2.2 mm/d. After redetermination, under-estimations of most stations were below 0.05 mm/d while few were around 0.18 mm/d. Future studies may divide Jiangxi into several districts according to the climate and align parameters within the specific district. In this way, applicability of a model may be improved without affecting its accuracy ( Table 2).
The redetermined parameter a of the H-S model ranged in 0.0009 ~ 0.0022 in Jiangxi and was lower than the default value (0.0023); parameter b of the rate-determined H-S model ranged in 0.62 ~ 0.92 and was higher than the default value (0.5), indicating that the daily temperature range contributed significantly to the reference crop evapotranspiration in Jiangxi; parameter c of the rate-determined H-S model ranged in 7.59 ~ 28.5 and the highest and second highest value appeared in the Lushan Mountain and the Jinggangshan Mountain, indicating significant effects of altitude on this parameter (Table 2).
Parameter a of the P-T method ranged in 1.05 ~ 1.18 and was lower than the default value (1.26), indicating that the contribution on the reference crop evapotranspira1tion by the aerodynamic term was relatively low when this method was applied in Jiangxi. Specifically, the contribution by the aerodynamic term was less than 20% of that by the radiation term without considering errors (Table 2).
In this study, after parameter redetermination of the I-A method, parameters b and c were set as 0.32 and 0.015, respectively, and parameter a ranged in − 0.69 ~ − 0.35. Similar to the I-A method, parameter b of the Mak method ranged in − 0.37 ~ − 0.08 with parameter a being set to be 0.66. Additionally, terms in the I-A method and the Mak method are constant and could only be used for correction of overall deviations of the ultimate results. As a result, MBE of the I-A method and the Mak method was 0 (see Fig. 3d).
Parameters a, b, and c of the PVB method varied significantly within − 378.6 ~ 681.9, − 113.7 ~ 60.2, and − 3.43 ~ 2.63, respectively. Hence, the esp calculated based on average temperature was highly unreliable, resulting in poor model accuracy ( Table 2).
The B-R method has four parameters. According to the calculation equation structure, coefficients of the high temperature term, the daily temperature range term, and the extraterrestrial radiation term were above 0.1. Despite their specific units, contributions of these three parameters to reference crop evapotranspiration were in similar scale. Also, redetermined coefficients of the extraterrestrial radiation term were basically less than half of the high temperature term, the daily temperature range term for all stations, indicating the dominant role of temperature in application of this model in Jiangxi. Contributions of the extraterrestrial radiation terms to reference crop evapotranspiration were small and the accuracy was relatively low (Table 2).
By introducing the difference term of the saturation vapor pressure and the actual vapor pressure into the B-R method, the M-B-R method has five parameters and average coefficient of the extraterrestrial radiation term was 20% higher than that of the daily temperature range term (Table 3). In other words, unlike the B-R method, the contribution of radiation on reference crop evapotranspiration was significant in the M-B-R method and its accuracy was significantly higher than that of the B-R method.
The parameters of the Sch method ranged in 10.0 ~ 11.4 and were lower than the default value (16) ( Table 3) However, the accuracy of this method was low, suggesting that this method neglects some key factors although the effects of temperature and relative humidity have been considered. As a result, applicability of the Sch method in Jiangxi was poor. For the Turc method, parameter b ranged in 15.65 ~ 17.47 if parameter a was set to be 0.018. Considering the fact that redetermination led to limited improvement on the model accuracy, it was concluded that the Turc method could be directly applied in Jiangxi. Parameters a, b, and c of the J-H method ranged in − 0.34 ~ 1.03, − 0.31 ~ 4.53, and − 0.14 ~ 2.12, respectively, and the parameter values varied significantly. This demonstrated that the J-H method was not reliable, and thus not applicable in Jiangxi (Table 3).
With parameter a of the B-S method set to be 0.69, parameters b and c ranged in − 2.36 ~ − 0.97 and 0.24 ~ 1.22, respectively, and the contribution of the wind speed on the aerodynamic term was adjusted by these two parameters. Moreover, as the aerodynamic term in this method was negative, parameter b should be negative to bring positive contribution to the aerodynamic term (Table 3).  Figure 4 shows deviations of monthly values by different evapotranspiration methods after redetermination. As observed, the monthly variation was less than 3 mm for the 48PM method, indicating good applicability. For the H-S method, deviations in April to August ranged in 4 ~ 10 mm and errors would be exacerbated in calculations of water demands of crops as crop coefficients were usually larger than 1 in this period. The absolute deviation of the P-T method was relatively high in October to December. However, this period was not the major water consumption period for crops; the P-T method was applicable in Jiangxi. The monthly absolute deviations of the I-A method were highly consistent (< 4 mm even in April to August), indicating the good applicability of the I-A method in Jiangxi. The absolute deviation of the Mak method was relatively high in July and August (under-estimation was 8.9 mm for July and 5.9 mm for August), demonstrating that the Mak method was inferior to the P-T method in terms of applicability in Jiangxi. The monthly absolute deviations of the PVB method were highly consistent (< 3 mm in May to September). The B-R method exhibited over-estimations of 6 mm in May to July and severe under-estimations in July and August. Compared with the B-R method, the M-B-R method exhibited significant reductions of monthly absolute deviations, although underestimations of 5 mm were observed in July and August. The monthly absolute deviations of the Sch method were above 5 mm and maximized in July to November. Redetermination by month may enhance the applicability of this model. The monthly absolute deviations of the Turc method were highly consistent (< 5 mm) and maximized deviations appeared in January and July. The J-H method exhibited under-estimations of 10 mm in October to April and over-estimations of 5 mm in July and August, demonstrating that the J-H method was not applicable in Jiangxi. The monthly absolute deviations of the B-S method were small, demonstrating good applicability of this method in Jiangxi.

Annual values by improved evapotranspiration methods
In the background of global climate change, model applicability may vary significantly with time. Figure 5 shows the model applicability in Jiangxi with interannual deviation, which facilitated the study. The annual ET 0 calculated by various methods exhibited significant variation around Therefore, the J-H method was not applicable in Jiangxi. The B-S method exhibited annual under-estimation in 0 ~ 10 mm, indicating that the overall deviation of the B-S method was minimum among that of all 12 methods. The long-term ET 0 curve of the B-S method was highly consistent with that of the P-M method.   Overall, the B-S method and the 48P-M method show optimized accuracy (maximized N RMSE < 5%), followed by the P-T method, the I-A method, the Mak method, the PVB method, and the Turc method (N RMSE = 0.1 ~ 0.2). In all these methods, overall surface radiation or net radiation has been considered and led to significant improvement for the accuracy. Among temperature-based methods, the J-H method was slightly better than the H-S method, while N RMSE of data by other methods ranged in 0.3 ~ 0.4, indicating that applicability of these methods was limited.

Discussion
The applicability of simplified ET0 models exhibits obvious regional differences due to various controlling factors in ET0 for different climates (Feng et al. 2016;Salam et al. 2020). Our study area, Jiangxi Province, has a subtropical monsoon climate, characterized by a high precipitation amount and a high spatiotemporal variability of precipitation and solar radiation. Previous studies have demonstrated soil radiation dominates ET0 in Jiangxi (Fan et al. 2018). Compared with other radiation-based models, the Turc model performed the best. The B-R, M-B-R, and J-H models have large errors because the default parameter values are not suitable for the local region.
The 48PM method has been widely applied globally (Li et al. 2016;Matsui & Osawa 2015;Xu et al. 2010), but it is limited by severe errors before parameter redetermination. For instance, R MSE of data by the 48PM method was between 1.18 and 2.13 mm/d (Li et al. 2016) in Sichuan, above 1.1 mm/d in Kunshan, Jiangsu (Xu et al. 2010), and about 0.5 mm/d in the middle and lower reaches of Yangtse River (Jia et al. 2016). Nevertheless, applicability of the B-S method has not been studied in China. In this study, R MSE of data by the 48PM method before parameter redetermination ranged in 1.25 ~ 3 mm/d in Jiangxi and the error was mainly attributed to over-estimation of the aerodynamic term. After parameter redetermination, R MSE of data by the 48PM method ranged in 0.06 ~ 0.34 mm/d, and the annual over-estimation of data by the 48PM method was below 30 mm and decreased since 2003. On the other hand, R MSE of data by the B-S method ranged in 0.08 ~ 0.27 mm and the annual over-estimation of data by the B-S method was below 20 mm. With comparable accuracy to the P-M method, the 48PM method and the B-S method involve no information about wind speed in the radiation term; thus, they are more suitable for investigations of contributions by the radiation term and the aerodynamic term to reference crop evapotranspiration and the mechanism behind. Tabari (2010) evaluated applicability of the Makkink method, the Turc method, the H-S method, and the P-T method in Iran, and proposed that the Turc method exhibited optimized applicability. Qin et al. (2016) reported that the Turc method showed good applicability in Qinhuai District only in September to November and no significant advantages in the upper reach of the Yellow River (Du et al. 2013;Li 2012). In this study, average R MSE and N RMSE of data by the Turc method redetermination were 0.39 mm/d and 0.157 m, respectively, while redetermination brought no significant enhancement. In summary, the Turc method could be directly applied in Jiangxi.
The P-T method is a calculation method for reference crop evapotranspiration commonly applied in humid regions. In this method, the aerodynamic term was neglected and the contribution by this term was covered by multiplying the radiation term by a coefficient larger than 1 (the default value is 1.26). Basically, ET 0 would be over-estimated in humid regions and under-estimated in arid regions (Estã Vez et al. 2010;Szilagyi 2014). In this study, Parameter α of the P-T method ranged in 1.05 ~ 1.18 and the method presents excellent performance in April to August. Therefore, the P-T method is an alternative method in absence of data about wind speed and relative humidity. Shared similar structure with the P-T method, the Mak method showed advantages in accuracy over the P-T method in Northwestern China (Wang et al. 2009;Wu et al. 2016). However, in this study, the Mak method was inferior to the P-T method in terms of accuracy and applicability, and it exhibited under-estimation by 6 mm and 9 mm in July and August, respectively.
It should be pointed out that the results by the P-M method were used as the reference data in this study. However, these data were not perfectly aligned with practical evapotranspiration data, and the reliable evaluations of applicability of different calculation methods for crop evapotranspiration in regional scale need further efforts. Bormann evaluated performances of 18 simplified evapotranspiration models in the background of climate changes (Bormann 2011). The results demonstrated that ET 0 is sensitive to meteorological factors to different degrees. The annual ET 0 calculated by different methods (except the 48PM method and the B-S method) exhibited significant variations around 2003. This may be attributed to significant changes in certain meteorological factors over recent years. The 48PM method and the B-S method shared identical climatological data with the P-M method, and the changes could be reflected by annual ET 0 . Also, it indicated that fixed empirical parameters in empirical methods could not reflect effects by key meteorological factors in long time simulations in the background of climate changes. The intensive studies on the causes and key meteorological factors of ET 0 in Jiangxi Province would facilitate the selection of appropriate empirical models, thus improving their applicability and reducing the risks.

Conclusions
The N RMSE was maximized in Northern Jiangxi for all methods, while minimized N RMSE by different models appeared in different regions. Performances of irradiation-based methods were better than those of temperature-based ones.
The 48PM method and the B-S method showed good applicability after parameter redetermination: their R MSE and N RMSE ranged in 0.06 ~ 0.34 mm/d and 0.08 ~ 0.27 mm, respectively, and their annual over-estimations were below 30 mm. R MSE and N RMSE of the Turc method were 0.389 mm/d and 0.157 mm, respectively, before parameter redetermination and the values did not vary significantly after parameter redetermination. The Turc method can be used in Jiangxi (except Nanchang). In the absence of relative humidity and wind speed, the P-T method was an appropriate simplified method for Jiangxi. In this case, α was between the default value 1.05 ~ 1.18, R MSE was within 0.21 ~ 0.66 mm/d, and N RMSE was within 0.08 ~ 0.30 m. Owing to severe monthly and annual deviations, the H-S method, the B-R method, and the M-B-R method showed poor accuracy.
In the background of climate changes, applications of empirical methods are challenged. The intensive studies of reasons on the causes and key meteorological factors of ET 0 in Jiangxi would facilitate selection of appropriate empirical models, thus improving their applicability and reducing the risks.
Code availability Not applicable.

Declarations
Ethics approval and consent to participate Not applicable.

Consent for publication Not applicable.
Competing Financial interests The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

Financial interests
The authors declare the following financial interests/personal relationships which may be considered as potential competing interests.

Conflict of interest
The authors declare no competing interests.