A Comparative Study of Potential Evapotranspiration Estimation by Eight Methods with FAO Penman – Monteith Method in Southwestern China

Potential evapotranspiration (PET) is crucial for water resources assessment. In this regard, the FAO (Food and Agriculture Organization)–Penman–Monteith method (PM) is commonly recognized as a standard method for PET estimation. However, due to requirement of detailed meteorological data, the application of PM is often constrained in many regions. Under such circumstances, an alternative method with similar efficiency to that of PM needs to be identified. In this study, three radiation-based methods, Makkink (Mak), Abtew (Abt), and Priestley–Taylor (PT), and five temperature-based methods, Hargreaves–Samani (HS), Thornthwaite (Tho), Hamon (Ham), Linacre (Lin), and Blaney–Criddle (BC), were compared with PM at yearly and seasonal scale, using long-term (50 years) data from 90 meteorology stations in southwest China. Indicators, viz. (videlicet) Nash–Sutcliffe efficiency (NSE), relative error (Re), normalized root mean squared error (NRMSE), and coefficient of determination (R2) were used to evaluate the performance of PET estimations by the above-mentioned eight methods. The results showed that the performance of the methods in PET estimation varied among regions; HS, PT, and Abt overestimated PET, while others underestimated. In Sichuan basin, Mak, Abt and HS yielded similar estimations to that of PM, while, in Yun-Gui plateau, Abt, Mak, HS, and PT showed better performances. Mak performed the best in the east Tibetan Plateau at yearly and seasonal scale, while HS showed a good performance in summer and autumn. In the arid river valley, HS, Mak, and Abt performed better than the others. On the other hand, Tho, Ham, Lin, and BC could not be used to estimate PET in some regions. In general, radiation-based methods for PET estimation performed better than temperature-based methods among the selected methods in the study area. Among the radiation-based methods, Mak performed the best, while HS showed the best performance among the temperature-based methods.


Introduction
Evapotranspiration (ET) plays an important role in maintaining water balance of terrestrial ecosystem.Accurate assessment of evapotranspiration is essential for efficient irrigation management, water resources management, crop production, environmental assessment, ecosystem modelers and solar energy system [1][2][3][4].Potential evapotranspiration (PET) has been commonly applied to calculate Water 2017, 9, 734 3 of 18 (Figure 1).Southwestern China is located in the first and second stairway of Chinese terrain with a complicated topography [39].There are four geomorphic units [40], viz.(a) Tibetan Plateau with an average attitude about 4000 m and stronger solar radiation, less rainfall and low temperature; (b) the Hengduan Mountains, consisting of a series of north-south oriented mountain ranges with altitudes of 4000-5000 m and major rivers, in which the temperature has obvious vertical variation, the annual average temperature on plateau surface is 14-16 • C while it is above 20 • C in valley bottom, and the wet season is between May and October; (c) Yun-Gui plateau with altitudes of 1800-1900 m, which has a subtropical monsoon climate with large temperature difference among seasons and the rainfall is concentrate in April to October; and (d) Sichuan basin with an elevation range of 300-700 m and annual average precipitation and temperature of 1000-1300 mm and 16-18 • C, respectively.Generally, southwestern China is a typical monsoonal climate region, which includes southwest monsoon and southeast monsoon.The precipitation and air temperature is daedal in spatial distribution while the dry and wet seasons are obvious.According to Kottek and coworkers' [41] classification theory about climate, southwestern China is Cmb (C means warm temperate; m means monsoonal; b means warm summer).
average attitude about 4000 m and stronger solar radiation, less rainfall and low temperature; (b) the Hengduan Mountains, consisting of a series of north-south oriented mountain ranges with altitudes of 4000-5000 m and major rivers, in which the temperature has obvious vertical variation, the annual average temperature on plateau surface is 14-16 C while it is above 20 C in valley bottom, and the wet season is between May and October; (c) Yun-Gui plateau with altitudes of 1800-1900 m, which has a subtropical monsoon climate with large temperature difference among seasons and the rainfall is concentrate in April to October; and (d) Sichuan basin with an elevation range of 300-700 m and annual average precipitation and temperature of 1000-1300 mm and 16-18C, respectively.Generally, southwestern China is a typical monsoonal climate region, which includes southwest monsoon and southeast monsoon.The precipitation and air temperature is daedal in spatial distribution while the dry and wet seasons are obvious.According to Kottek and coworkers' [41] classification theory about climate, southwestern China is Cmb (C means warm temperate; m means monsoonal; b means warm summer).

Data Sources
Because of the inconsistency in data length and data integrity, we selected 90 comparable stations with complete daily meteorological data from 1962 to 2013, and all the stations are national standard stations.The modern nation-wide network of weather observing stations in China began operation in the 1950s.The meteorological data were downloaded from the National Climate Center, China Meteorological Administration (CMA) [42] and included daily precipitation, relative humidity, average temperature, maximum temperature, minimum temperature, wind speed, etc.To detect the best model for a region, the study region was divided into four sub-regions (Figure 1): Sichuan basin (19 stations), located in the northeast of the study region; Yun-Gui plateau

Data Sources
Because of the inconsistency in data length and data integrity, we selected 90 comparable stations with complete daily meteorological data from 1962 to 2013, and all the stations are national standard stations.The modern nation-wide network of weather observing stations in China began operation in the 1950s.The meteorological data were downloaded from the National Climate Center, China Meteorological Administration (CMA) [42] and included daily precipitation, relative humidity, average temperature, maximum temperature, minimum temperature, wind speed, etc.
To detect the best model for a region, the study region was divided into four sub-regions (Figure 1): Sichuan basin (19 stations), located in the northeast of the study region; Yun-Gui plateau (39 stations), located in the south of the study region; the eastern margin of the Tibetan Plateau (11 stations), located in the northwest of the study region [43]; and the arid river valley region (21 stations), which belongs to the eastern margin of the Tibetan Plateau, which is made a separate region because of special climate characteristics [44], such as a small annual but a large diurnal variation of air temperature, abundant precipitation, and significantly dry and humid seasons.Due to high evaporation, the soil is deficient in water for almost the whole year, especially during the growth season.This area is also provided with abundant sunshine and intense solar radiation.The annual and seasonal values of PET were calculated in the above-mentioned four sub-regions.The seasonal scale was divided into spring (March-May), summer (June-August), autumn (September-November), and winter (December-February).All the calculated values were based on daily or monthly meteorological data; however, the present study has accumulated them into seasonal and annual values.

Daily Based methods
The FAO-PM method as given by FAO Irrigation and Drainage Paper No. 56 [10] as where PET PM is the potential evapotranspiration (mm/d); ∆ is the slope of the saturation vapor pressure function (kPa/ • C); R n is the net radiation (MJ/m 2 /day) (MJ means megajoule), which was estimated from total incoming solar radiation measurements following the procedure of Allen et al. [10]; G is the soil heat flux density (MJ/m 2 /day), which was considered as null for daily estimates; γ is the psychometric constant (kPa/ • C); T mean is the daily average temperature ( • C), which is the average value of the sum of maximum and minimum temperature; u 2 is the wind speed at 2 m height (m/s); e s is the vapor pressure of the air at saturation (kPa); and e a is the actual vapor pressure (kPa).
Hargreaves and Samni [30] proposed several equations for calculating daily potential evapotranspiration.One of the equations is which, when compared with modified Penman, showed reasonable accuracy to estimate the reference crop evapotranspiration [30].In Equation ( 2), R a is the extraterrestrial solar radiation (MJ/m 2 /day), which was measured following the procedure of Allen et al. [10]; T max and T min are the maximum and minimum temperatures ( • C), respectively; and T is the mean air temperature ( • C).
Hamon [32] computed ET in millimeters per day as where Ld is the daytime length (h); RHOSAT is the saturated vapor density (g/m 3 ); and KPEC is calibration coefficient, which we set as 1.2 in the current study [11].Priestley and Taylor [34] proposed a simplified version of PET equation as where α is the calibration constant; α = 1.26; and λ is the latent heat of vaporization (MJ/kg), which we set as 2.45 in the present study [45].
Water 2017, 9, 734 5 of 18 Linacre [33] developed an equation including air and dew point temperature where A is elevation of the station (m); ϕ is the latitude of the station ( • C); and T d is dew point temperature ( • C).
Makkink [35] estimated PET in millimeters per day as where R s is solar radiation (MJ/m 2 /day), which was measured following the procedure of Allen et al. [10], and relates surface shortwave radiation to extraterrestrial radiation and daily sunshine duration: where a s and b s are regression constants, the recommended values a s = 0.25 and b s = 0.50 were used in this study; n is daily sunshine duration (h); and N is daylight hours (h).
Abtew [29] used a simple model that estimates ET from solar radiation as follows where k is a dimensionless coefficient = 0.53.

Monthly Based Methods
Thornthwaite [28] reported a PET equation based on monthly meteorological data as where T is monthly average temperature ( • C); H is annual heat index; and A is constant.Brouwer and Heibloem [31] set Blaney-Criddle equation in their study as where q is daily percentage of annual daytime hours; and T m is the mean monthly temperature ( • C).

Statistical Analysis
To assess the best method to evaluate PET, four performance indicators, viz.Nash-Sutcliffe efficiency (NSE), relative error (Re), normalized root mean squared error (NRMSE) and linear regression [46,47], were used.The value of NSE, Re and NRMSE were calculated by the following equations and linear regression was obtained using statistical software, viz.SigmaPlot 12.5 (https://systatsoftware.com/; Systat software Inc., San Jose, CA, USA).
Water 2017, 9, 734 PET PM (13) where n is the number of PET value sample; PET PM is the PET value calculated by PM (mm); and PET o is the PET value calculated by other methods (mm).Good performance of the methods are required to satisfy all the following conditions: 0 ≤ NSE ≤ 1; 0 ≤ |Re| ≤ 0.2; 0 ≤ NRMSE ≤ 0.2; and 0.8 ≤ R 2 ≤ 1.
Percent bias (PBIAS) was used to measure the average tendency of the simulated data to be larger or smaller than standard estimated data [48], which is calculated as below: PBIAS is expressed as percentage, and positive values show model underestimation bias, while negative values show model overestimation bias.

PET Estimated by FAO-PM Method in the Four Sub-Regions
Annual and seasonal PET values estimated by FAO-PM method in the four sub-regions are shown in Table 1.Yun-Gui plateau displayed the maximum value in yearly, autumn and winter, whereas the arid river valley region showed the maximum value in other two seasons.Sichuan basin showed the minimum values at all time except summer, whereas the eastern margin of the Tibetan Plateau showed the minimum value for summer.Overestimation and underestimation errors in the four sub-regions of all time scales are shown in Table 2. PBIAS suggested HS, PT, and Abt overestimate PET in all regions, while the five other methods underestimate PET.However, HS showed underestimation in autumn and winter in the eastern margin of the Tibetan Plateau, which were 0.37% and 2.87%, respectively.According to the PBIAS values, HS and Mak ranged from −21.63% to 2.87% and from 6.44% to 15.97%, respectively, and they showed better performance than other methods in all regions and all time scales.The PET value of Tho was much lower than PM due to the lowest PBIAS 67.02% and the highest PBIAS value even up to 92.18%.

Sichuan Basin
To understand the relationship between PM and the eight other methods, NSE, Re, and NRMSE values for one year and four seasons are listed in Table 3. Yearly, NSE values ranged from −58.724 to −1.040.Results for Abt and Mak were −1.040 and −1.857, respectively, while the others were much smaller than zero.Compared to Re, Abt, Mak and HS performed better with values 0.110, −0.155 and 0.188, respectively.With respect to the NRMSE values, Abt performed the best with the lowest value of 0.134.Considering the coefficient of determination (R 2 ) of the linear regression between PM and the eight other methods (Figure 2e), Ham performed the best with the value of 0.925 followed by PT (0.920) and Mak (0.899).The results indicated that no methods performed well yearly, but Mak and Abt were narrowly acceptable.The seasonal statistical analysis between the measured and simulated PET are shown in Table 3 and Figure 2.For spring, Abt performed the best with the highest value of NSE (0.664), the lowest value of |Re| (0.075), the lowest value of NRMSE (0.129), and a higher value of R 2 (0.836).Mak was in the second place with NSE of 0.477, R 2 of 0.982, |Re| of 0.152 and NRMSE of 0.161.Although the R 2 value of PT was close to 1, the results of other statistical analysis showed poor performance.It indicated that PT is unacceptable for PET simulation for spring.2b.Mak performed the second best with the values of NSE = 0.691, |Re| = 0.160, NRMSE = 0.170, and R 2 = 0.992.
For autumn, the applicability of each method was not similar to that for yearly, spring, and summer.Mak performed the best with the highest value of NSE (0.630), the lowest value of |Re| (0.155), the lowest value of NRMSE (0.164), and the second highest value of R 2 (0.980) (Table 3 and Figure 2c).HS performed the second best with NSE = 0.552, |Re| = 0.155, NRMSE = 0.181, and R 2 = 0.934.Although PT owned the highest R 2 (0.981), other statistical analysis of PT appeared poor which indicated that PT performed not well in autumn.
For winter, Mak was the best estimation method with the highest value of NSE (0.886), the smallest value of |Re| (0.148), the smallest value of NRMSE (0.173), and the value of R 2 = 0.7797, which indicated that this method was acceptable and available.

Yun-Gui Plateau
As seen in Table 3 and Figure 3e, for year, the NSE values of Abt and Mak were 0.334 and 0.138, respectively, and the Re and NRMSE values were 0.110 and −0.136, and 0.123 and 0.140, respectively.It indicated that these two methods estimated PET well.The values of R 2 of PT, Mak, Ham, and Abt indicated that these methods were efficient (R 2 ＞ 0.850).The overall results showed that Abt was the best method to perform year PET in this region, followed by Mak.

Yun-Gui Plateau
As seen in Table 3 and Figure 3e, for year, the NSE values of Abt and Mak were 0.334 and 0.138, respectively, and the Re and NRMSE values were 0.110 and −0.136, and 0.123 and 0.140, respectively.It indicated that these two methods estimated PET well.The values of R 2 of PT, Mak, Ham, and Abt indicated that these methods were efficient (R 2 >0.850).The overall results showed that Abt was the best method to perform year PET in this region, followed by Mak.
For spring, the NSE values of Abt, HS and Mak were above 0.2, especially for Abt whose NSE was 0.845 (Table 4 and Figure 3a).The Re and NRMSE of Abt, Mak and HS were lower.The Re value ranged from −0.140 to 0.185, and the NRMSE was in the range 0.127-0.217.As for Mak, PT, Abt and HS, the R 2 values were above 0.93, which indicated that these four methods performed well.These results showed that Abt was the best estimate method to simulate PET for the region in spring, followed by Mak and HS.Although PT owned the highest value of R 2 (0.989), its other statistical parameter values were not useful, indicating that this estimation model is not acceptable and applicable.NSE values of four methods performed well among which the values of Mak and HS were 0.886 and 0.803, respectively.For these methods, the Re values ranged from −0.126 to 0.185, while NRMSE ranged from 0.138 to 0.201.The values of R 2 were above 0.80 except for Tho method; the highest value of R 2 was shown by PT.For winter, Mak was the best choice followed by HS and Abt.For summer (Table 4 and Figure 3b), the NSE values were above 0.10 for Abt, Mak and HS.Abt showed the highest NSE value of 0.805.The Re values ranged from −0.145 to 0.168 for Abt, Mak and HS, while the NRMSE values were between 0.089 and 0.190.Note that the Re and NRMSE values of Abt were close to zero.The R 2 value of PT was the highest (0.979) followed by Mak (0.978), both close to 1. Similar for year and spring, Abt displayed the best performance followed by Mak and HS.
As for autumn (Table 4 and Figure 3c), the NSE values of Mak, HS, and Abt were higher than 0.58.The Re and NRMSE values were low for Mak, HS and Abt.The values of Re ranged from −0.132 to 0.136, while the values of NRMSE ranged from 0.144 to 0.180.The values of R 2 were higher than 0.830 for Mak, HS and Abt.These results indicated that Mak, HS and Abt were all acceptable and applicable, among which Mak was the best.
NSE values of four methods performed well among which the values of Mak and HS were 0.886 and 0.803, respectively.For these methods, the Re values ranged from −0.126 to 0.185, while NRMSE ranged from 0.138 to 0.201.The values of R 2 were above 0.80 except for Tho method; the highest value of R 2 was shown by PT.For winter, Mak was the best choice followed by HS and Abt.

The Eastern Margin of the Tibetan Plateau
A comparative evaluation of error and regression analysis indicated that Mak performed the best for yearly scale with the highest NSE (0.238), the second lowest |Re| (0.095), the lowest NRMSE (0.120) and high R 2 value of 0.983 (Table 5 and Figure 4e).The other methods were not quite suitable for yearly scale in this region.
The spring results were similar to those of yearly; only Mak was acceptable and applicable with the highest NSE (0.049), the lowest |Re| (0.100) and NRMSE (0.111), and the second-highest R 2 (0.830) (Table 5 and Figure 4a).The NSE values were not good for spring and yearly scale because, not only was the value of Mak close to 0 (lying above 0), but the value of the seven other methods were also very low, lying below 0. The R 2 values were also low; most of them were below 0.9.
The Mak and HS performed the best and second best, respectively, in summer and autumn (Table 5 and Figure 4b,c).In summer, Mak showed the highest NSE (0.747) and second-highest R 2 (0.983), the lowest |Re| (0.123) and the lowest NRMSE (0.142).HS also could be used to estimate PET with the second-highest NSE (0.351), high R 2 (0.924), the second-lowest |Re| (0.133) and the second-lowest NRMSE (0.227).In autumn, the NSE values of Mak and HS were above 0.6, which generated good simulation results, 0.849 and 0.665, respectively.The Re values of Mak and HS were closed to 0, especially HS whose Re value was −0.004.Mak had the lowest NRMSE (0.135), while HS owned the second-lowest value (0.201).Five methods had R 2 values above 0.8, and PT owned the highest value of 0.977.These results showed that Mak was the best method to estimate PET followed by HS in summer and autumn.In winter (Table 5 and Figure 4e), the NSE values of Mak and HS were 0.860 and 0.738, respectively, both of them close to 1.The Re values of Mak (−0.064) and HS (−0.029) were very close to zero.The NRMSE values of Mak and HS were 0.090 and 0.123, respectively.The highest R 2 value was 0.960, which belonged to PT, while Mak showed the second-highest value of 0.934; however, HS showed a low R 2 value (0.766), which was below 0.8.These results indicated that Mak was the best method to estimate PET in winter in this region.

Arid River Valley Region
For yearly scale (Table 6 and Figure 5e), only NSE values of two methods were above 0, viz.HS (0.433) and Mak (0.344).The Re value of HS was close to zero, while for Mak, it was −0.120.As for NRMSE value, HS had the lowest value (0.121) and Mak owned the second-lowest (0.131).Four methods had R 2 values above 0.8; PT had the highest value while HS had the lowest value among these four methods.These results indicated that HS was the best and Mak was the second-best for yearly scale in this region.In winter (Table 5 and Figure 4e), the NSE values of Mak and HS were 0.860 and 0.738, respectively, both of them close to 1.The Re values of Mak (−0.064) and HS (−0.029) were very close to zero.The NRMSE values of Mak and HS were 0.090 and 0.123, respectively.The highest R 2 value was 0.960, which belonged to PT, while Mak showed the second-highest value of 0.934; however, HS showed a low R 2 value (0.766), which was below 0.8.These results indicated that Mak was the best method to estimate PET in winter in this region.

Arid River Valley Region
For yearly scale (Table 6 and Figure 5e), only NSE values of two methods were above 0, viz.HS (0.433) and Mak (0.344).The Re value of HS was close to zero, while for Mak, it was −0.120.As for NRMSE value, HS had the lowest value (0.121) and Mak owned the second-lowest (0.131).Four methods had R 2 values above 0.8; PT had the highest value while HS had the lowest value among these four methods.These results indicated that HS was the best and Mak was the second-best for yearly scale in this region.The spring simulation result was contrary to the yearly results (Table 6 and Figure 5a).Mak, with the highest NSE (0.600) and R 2 (0.964), the lowest NRMSE (0.129) and the second-lowest |Re| (0.117), was the best performing method.The second was HS, which had the lowest |Re| (0.114), the second-highest NSE (0.449), the second-lowest NRMSE (0.152) and high R 2 (0.872).
For summer (Table 6 and Figure 5b), the highest NSE and high R 2 values was obtained for Abt, about 0.677 and 0.893, respectively, as well as the lowest |Re| and NRMSE, about 0.097 and 0.134, For summer (Table 6 and Figure 5b), the highest NSE and high R 2 values was obtained for Abt, about 0.677 and 0.893, respectively, as well as the lowest |Re| and NRMSE, about 0.097 and 0.134, respectively, indicating that Abt performed the best.Mak was the second with the highest R 2 values, about 0.980, the second-highest NSE (0.628), and low |Re| (0.134) and NRMSE (0.144).The third was HS with high NSE (0.548) and R 2 (0.921), and low |Re| (0.130) and NRMSE (0.150).

Discussion
This paper compared eight different methods for the evaluation of PET with PM in southwestern China.We applied three radiation-based methods (PT, Mak, and Abt) and five temperature-based methods (Tho, Ham, Lin, HS, and BC), all of which showed large differences in their applicability.The radiation-based methods showed better performance than the temperature-based methods in the evaluation of PET in the study area.This is in accordance with findings of previous studies in other regions [25,49,50].
Among the radiation-based methods, Mak showed the best performance followed by the Abt and PT.These three methods were established in a humid climate condition, which are suitable for regions with abundant rainfall like our study area.This is probably one of the reasons they can be used to estimate PET in the study area.Mak showed less errors and high correlation comparable with PM among the radiation-based methods used in the present study.Lu et al. [11] suggested Mak with poor performance in the southeastern United States, whereas Tukimat et al. [49] reported that Mak performed the best among the studied methods in Malaysia.The reason for these different results may be that the same method yielded different results under different climatic conditions and geographical environments.PT was acceptable and applicable to evaluate PET during the winter in Yun-Gui plateau because of its third place among the three methods which can evaluate PET in this situation.Wei and Menzel [51] believed PT to be the most suitable method for global application in the estimation of PET.However, Abt only needs two parameters to estimate ET compared to PT which requires five parameters.Therefore, Abt can be used more easily than PT in estimating PET in the study area.
Among the temperature-based methods, only HS can be used to estimate PET, but it was not the best simulation method in the study area.HS can be used on all time scales in arid river valley region, possibly because it was established in an arid and semi-arid climate condition, and the present study area has a similar climate condition to some extent, i.e., the arid river valley region.Meanwhile, the input parameters of the HS included extraterrestrial solar radiation while other temperature-based methods exclude this, which might be the reason that only HS can be used to estimate PET in the study region.Chen et al. [13], Lu et al. [11] and Nikam et al. [52] found that Tho performed worse in many regions.In the present study, Tho also showed the worst performance among all the selected methods.This is probably because Tho is only based on temperature, i.e.only single input parameter.Another reason might be that Tho was established in a humid climate and the underlying surface was a valley, and PET is determined by other variables such as humidity, wind speed, vapor pressure and solar radiation.However, Tho was found to be a suitable method for PET estimation in Malaysia [48,53].Besides, other temperature-based methods were also found to be unsatisfactory in PET estimation in the study area due to the high error analysis and poor correlation.
The four statistical indicators showed various simulation effects.On a seasonal scale, NSE showed simulation methods match FAO-PM method better than annual except for the eastern margin of the Tibetan Plateau.The greatest disadvantage of NSE was that calculations of squared values were required to show the relationships between the measured and the simulated values.As a result, larger values were strongly overestimated while the lower values were neglected [54,55].The three other indicators presented no special features in annual and seasonal scales.Overall, the simulation effects of the four statistical indicators in the eastern margin of the Tibetan Plateau were poorer than the other three regions.The high altitude, changeable weather and complex terrain may be responsible for such results.PBIAS suggested that Abt, PT, and HS overestimate PET, while the other five models typically underestimate PET.
All the selected methods comparable with PM showed different simulation results, in which some methods viz.Mak and Abt led to good estimation, some viz.PT and HS were acceptable and applicable, while some viz.Ham, Lin, Tho and BC were completely unusable to estimate PET.The input parameters may be the main factors that affected their performance.Different methods need different input parameters.Another important reason is that each method was established on a specific underlying surface and climate conditions; however, the study area has complex topography and changeable weather.The heterogeneity problem in the climate data we used may be influenced by climate or non-climate factors [56,57] (including the relocation of weather stations and observation methods, urbanization, and agricultural irrigation).

Conclusions
From the above analysis, it can be concluded that the radiation-based methods were more suitable for estimating PET in southwest China.Among the radiation-based methods, Mak performed the best, Abt was second, and PT was third.Among the temperature-based methods, HS was the only one which was acceptable and applicable, while the others showed a poor performance in evaluating PET, especially Tho.Besides, the PET values of HS, PT, and Abt were lager than PM, while Ham, Lin, Mak, Tho, and BC were smaller than PM.
For Sichuan basin, Mak performed the best for year, autumn and winter PET estimations, while Abt was the most suitable for spring and summer.For the Yun-Gui plateau, Abt was perfect for estimating in year, spring and summer PET, while, Mak performed the best for autumn and winter.In the eastern margin of the Tibetan Plateau, Mak was the most suitable for all time scales.On the other hand, for the arid river valley region, HS ranked first in year and autumn PET estimations, Mak performed the best for spring and winter, and Abt was perfect for summer.
To some extent, it is expected that the conclusions of this study can be used in regions with similar topography and climatic conditions in the world.If a region is featured with low latitude, warm and moist climate similar to Sichuan basin, or if a region is a plateau having monsoon climate with low latitude, Mak and Abt can be used to estimate PET.Researchers may be able to use Mak in regions that are similar to the eastern margin of the Tibetan Plateau.If a region has abundant rainfall but high evaporation, and intensive solar radiation similar to the arid river valley region in this study, Mak and HS may be suitable for PET estimating.A region whose geographic feature is complex and climate varies greatly similar to southwest China, Mak can be recommended.We hope that the present study will be helpful in general to select the appropriate methods according to the availability of meteorological data.

Figure 1 .
Figure 1.Location of the weather stations in southwestern China.(Thered lines were dividing lines which divided the study region into three zones; and  in red represents the weather stations in arid river valley region which belongs to one of the three parts.)

Figure 1 .
Figure 1.Location of the weather stations in southwestern China.(The red lines were dividing lines which divided the study region into three zones; and ∆ in red represents the weather stations in arid river valley region which belongs to one of the three parts.)

Figure 4 .
Figure 4. Linear regression of measured and simulation values in the eastern margin of the Tibetan Plateau: (a) spring; (b) summer; (c) autumn; (d) winter; and (e) annual.

Figure 5 .
Figure 5. Linear regression of measured and simulation values in the arid river valley region: (a) spring; (b) summer; (c) autumn; (d) winter; and (e) annual.

Figure 5 .
Figure 5. Linear regression of measured and simulation values in the arid river valley region: (a) spring; (b) summer; (c) autumn; (d) winter; and (e) annual.

Table 1 .
Annual and seasonal PET estimated by FAO-PM method in the four sub-regions.

Table 2 .
Annual and seasonal PBIAS (%) in the four sub-regions a .
a : SCB means Sichuan basin; YGP means Yun-Gui plateau; ETP means the eastern margin of the Tibetan Plateau; ARV means arid river valley region

Table 3 .
Accuracy assessment of PET estimation in Sichuan basin.Abt performed the best with the highest value of NSE (0.937), the lowest values of |Re| (0.007), the lowest value of NRMSE (0.076,) and high value of R 2 (0.964), as shown in

Table 3 and
Figure

Table 4 .
Accuracy assessment of PET estimation in the Yun-Gui plateau.

Table 5 .
Accuracy assessment of PET estimation in the eastern margin of the Tibetan Plateau.

Table 6 .
Accuracy assessment of PET estimation in the arid river valley region.