Spatial and Temporal Variations of Potential Evapotranspiration in the Loess Plateau of China During 1960–2017

Potential evapotranspiration (ET0) is an integral component of the hydrological cycle and the global energy balance, and its long-term variation is of much concern in climate change studies. The Loess Plateau is an important area of agricultural civilization and water resources research. This study analyzed the spatial and temporal evolution processes and influential parameters of ET0 at 70 stations in different topographical areas of the Chinese Loess Plateau (CLP). Using the Mann–Kendall trend, Cross wavelet transform, and the ArcGIS platform, the ET0 of each station was quantified using the Penman–Monteith equation, and the effects of climatic factors on ET0 were assessed by analyzing the correlation coefficients and contribution rates of the climatic factors. The results showed that: (1) the overall trend of the ET0 in different terrains of the Loess Plateau is consistent, however, the ET0 values differ; the hill region (HR) has the highest ET0, followed by the valley region (VR), and the mountain region (MR) has the lowest, and ET0 changes differ between seasons. (2) Spatial distribution characteristics of multiyear mean ET0 in the study are as follows: the ET0 values in mountain and hilly areas are decreasing from west to east, and the higher mean annual ET0 value in the VR is mainly concentrated in the eastern CLP. (3) In the past 58 years, the annual mean and the seasonal ET0 of the region showed increasing trends, however, differences in different terrains were obvious. (4) ET0 has significant correlations with El Niño–Southern Oscillation (ENSO), Pacific–North American teleconnection (PNA), and Atlantic Multidecadal Oscillation (AMO). The resonance period of ET0 and ENSO was 3–6 a, mainly in 1976–1985. The mean coherence phase angle was close to 360◦, indicating that ET0 lags behind PNA by approximately 2–6 a; ET0 has a very strong positive correlation with AMO. (5) Relative humidity (RH) is the main influencing factor of ET0 change in the Loess Plateau. Temperature (T) variation has the highest contribution rate (42%) to the regional ET0 variation in the entire CLP. We should pay more attention to the variation of evaporation under future climate change, especially temperature change.


Introduction
Potential evapotranspiration (ET 0 ) is an important parameter of hydrological processes and is closely related to regional precipitation and energy spatial distribution [1,2]. With the changing global climate, rising temperatures will further result in increasing rates of water circulation which further change regional climate patterns [3][4][5]. The global temperature has a significant increasing place of origin of world agricultural civilization, the fragile ecological environment in the CLP is very sensitive to climate change. Therefore, it is necessary to improve our understanding of the temporal and spatial characteristics of regional ET 0 in the CLP. In this study, data from 70 weather stations with three topography types in the CLP during 1960-2017 were used to analyze spatial and temporal ET 0 variations. The main objectives of this study were to: (1) analyze the temporal and spatial variation characteristics of the ET 0 with different topographies in the CLP; (2) determine the main controlling factors of the variations of regional ET 0 with different topographies in the CLP; and (3) determine the relationship between regional ET 0 of the CLP and larger-scale circulation parameters. The results of this study will be of great significance to improve our understanding of regional hydrological processes and regional water resource management.

Study Area
In the current study, the CLP was defined as the area between 33 • 43'-41 • 16' N and 100 • 54'-114 • 33' E [36] and has three geomorphic units: the mountain region (MR, including Qin Ling Mountains in the southern CLP, Wu Shao Ling Mountains in the western CLP, and Taihang Mountains and Luliang Mountains in the eastern CLP), the valley region (VR, including Hetao Basin in the northern CLP, Weihe Basin in the southern CLP, Yinchuan Basin in the western CLP, and Fenhe Basin in the eastern CLP), and the hill region (HR, mainly concentrated in the inner part of the CLP) ( Figure 1). The area is approximately 64 x 10 4 km 2 , which is approximately 6.67% of the total area of China. The terrain of the CLP is high in the southwest region and low in the southeast. Based on long-term observation data from 70 national meteorological observatory stations in and around the CLP, we found the mean annual temperatures ranging from 3.6 to 14.3 • C and annual precipitation ranging from 150 to 750 mm. The soil type is mainly brown soil and cinnamon soil. Vegetation types are diverse and are mainly grassland, woodland (the categories of trees in the CLP mainly include Chinese catalpa, cypress, and Chinese prickly ash), and arable land [37]. The CLP has a long history of agricultural cultivation and its main types of crop producation are wheat and millet [38]. of the temporal and spatial characteristics of regional ET0 in the CLP. In this study, data from 70 weather stations with three topography types in the CLP during 1960-2017 were used to analyze spatial and temporal ET0 variations. The main objectives of this study were to: (1) analyze the temporal and spatial variation characteristics of the ET0 with different topographies in the CLP; (2) determine the main controlling factors of the variations of regional ET0 with different topographies in the CLP; and (3) determine the relationship between regional ET0 of the CLP and larger-scale circulation parameters. The results of this study will be of great significance to improve our understanding of regional hydrological processes and regional water resource management.

Study Area
In the current study, the CLP was defined as the area between 33°43'-41°16' N and 100°54'-114°33' E [36] and has three geomorphic units: the mountain region (MR, including Qin Ling Mountains in the southern CLP, Wu Shao Ling Mountains in the western CLP, and Taihang Mountains and Luliang Mountains in the eastern CLP), the valley region (VR, including Hetao Basin in the northern CLP, Weihe Basin in the southern CLP, Yinchuan Basin in the western CLP, and Fenhe Basin in the eastern CLP), and the hill region (HR, mainly concentrated in the inner part of the CLP) ( Figure 1). The area is approximately 64 x 10 4 km 2 , which is approximately 6.67% of the total area of China. The terrain of the CLP is high in the southwest region and low in the southeast. Based on longterm observation data from 70 national meteorological observatory stations in and around the CLP, we found the mean annual temperatures ranging from 3.6 to 14.3 °C and annual precipitation ranging from 150 to 750 mm. The soil type is mainly brown soil and cinnamon soil. Vegetation types are diverse and are mainly grassland, woodland (the categories of trees in the CLP mainly include Chinese catalpa, cypress, and Chinese prickly ash), and arable land [37]. The CLP has a long history of agricultural cultivation and its main types of crop producation are wheat and millet [38].

Data
This study used continuous daily series data from 1960 to 2017 at 70 national meteorological observatory stations in and around the CLP (Table 1). Data included wind speed (WS) (m/s), sunshine duration (SD) (h), daily mean temperature (T) ( • C), daily maximum temperature (Tmax) ( • C), daily minimum temperature (Tmin) ( • C), and relative humidity (RH) (%). Data were provided by the China Meteorological Administration (http://data.cma.cn/site/index.html). Large-scale circulation parameters data were downloaded from the National Oceanic and Atmospheric Administration (http://www.esrl.noaa.gov/psd/data/climateindices). Four seasons were divided in this study as: spring (March to May), summer (June to August), autumn (September to November), and winter (December to February).

ET 0 Calculation Method
In this study, the FAO Penman-Monteith method [13] was used to estimate ET 0 variations. The FAO-PM equation to calculate ET 0 (mm.day −1 ) is given as: where ET 0 is the daily potential evapotranspiration rate, mm·d −1 ; ∆ is the slope of the saturation vapor pressure curve, KP a · • C −1 ; Rn is the net radiation, MJ·m −2 ·d −1 ; G is the soil heat flux density, MJ·m −2 ·d −1 ; γ is the psychometric constant, KP a · • C −1 ; u 2 is the wind speed at the height of 2 m, m·s −1 ; T is the temperature at the height of 2 m, • C; e s is the saturated water vapor pressure; e a is the actual water vapor pressure, and (e s − e a ) is the saturated water vapor pressure differential value, KP a .

Mann-Kendall Trend Test
The Mann-Kendall (MK) test was used to detect regional ET 0 trends and hydrological imbalances [39,40]. For a given data series, the statistical value S and the standardized test statistics Z MK were calculated as follows [41]: Var(S) = 1 18 and where X i and X j are the values of the year i and j of the time series, n is the length of the time series data, t p is the binding value corresponding to the pth number, and Z MK is the variation trend of the time series data. If Z MK > 0, the time series data have an increasing trend over time; if the contrary, the data have a decreasing trend. If |Z MK | > Z (1−a/2) , then the null hypothesis that there is a significant trend in the time series data was rejected. The Z (1−a/2) values can be found in the standard normal distribution chart. When taking a = 5% as the significance level, the corresponding value of Z (1−a/2) is 1.96; when taking a = 1% as the significance level, the corresponding value of Z (1−a/2) is 2.58.

Wavelet Transform Coherence and Contribution Rate Analysis
Wavelet transform coherence is a new signal analysis technology combining wavelet transform and cross-spectrum analysis. It can study the correlation between two time series in the time-frequency domain from multiple time scales [42]. Wavelet transform coherence was introduced by Hudgins et al. [43] among others. In hydrology, it has been used in rainfall-runoff cross-analysis [44].
Here, W XY n (s) indicates the cross wavelet power spectrum, S is the smoothing function, c 1 , c 2 are the normalized constants, and (0.6s) indicates a rectangular function. In the formula, Y is the standardized value of the ET 0 , X 1 , X 2 · · · X n are standardized values of the meteorological parameters, a 1 , a 2 · · · a n are regression coefficients after sequence normalization, g 1 is the relative contribution rate of X n variation to Y variation. Y = a 1 X 1 + a 2 X 2 + a 3 X 3 + · · · + a n X n (9) Figure 2 shows the seasonal and annual mean ET 0 over the entire CLP and the three topographic regions from 1960 to 2017. The annual mean ET 0 was 1074 mm during 1960-2017 in the CLP, and the maximum and minimum values were 1181 mm and 987 mm, respectively. Comparing the ET 0 values of the three topographic regions (MR, HR, and VR), the mean annual ET 0 value of the HR (1078 mm) was larger than that of the MR (933 mm) and the VR (1021 mm). Over recent decades, the mean annual ET 0 value of the CLP displayed an upward trend of 7.7 mm/10a. There was a decreasing trend during the 1960s and 1970s. During this period, the ET 0 value exhibited a decreasing trend of fluctuations from the 1960s (1074 mm) to the 1970s (1052 mm), and then, a relatively stable trend of regional ET 0 during 1960-2017 was higher than that for winter. A larger temporal variation can be observed in the mean ET0 values for the spring and autumn. For spring (Figure 2b), the ET0 exhibited a relatively consistent value from 1960 to 1990, after which time a remarkable increase can be seen since the 1990s. A similar variation appeared in the mean ET0 for the summer and autumn. A slightly decreasing trend was displayed in the mean ET0 values for the summer and autumn during the period of the 1960s to the mid-1980s. From the mid-1980s, the mean ET0 for the summer and autumn showed increasing trends. The temporal variation of the mean ET0 in winter showed a similar trend to that in the spring. Regarding the three topographic regions (MR, HR, and VR), the mean ET0 values for the three regions exhibited similar temporal variations in all seasons except for the summer. In summer, the mean ET0 of the MR showed a slightly increasing trend (Figure 2c).  Figure 3 shows the spatial patterns of the seasonal and annual ET0 values mean from 1960 to 2017. During 1960-2017, the annual mean ET0 of the CLP was 1044 mm, ranging from 495 to 1467mm. The highest value of ET0 was at the Xiangning station (1467 mm), and the lowest value was at the Wutaishan station (495 mm). The ET0 values in the southern and eastern areas of the CLP were always higher than those in the western part of the CLP (Figure 3), which may be related to the different During the study period, the mean ET 0 values in spring, summer, autumn, and winter were 250, 445, 247, and 97mm, respectively. Comparing the four seasons, the mean ET 0 value for summer during 1960-2017 was higher than that for winter. A larger temporal variation can be observed in the mean ET 0 values for the spring and autumn. For spring (Figure 2b), the ET 0 exhibited a relatively consistent value from 1960 to 1990, after which time a remarkable increase can be seen since the 1990s. A similar variation appeared in the mean ET 0 for the summer and autumn. A slightly decreasing trend was displayed in the mean ET 0 values for the summer and autumn during the period of the 1960s to the mid-1980s. From the mid-1980s, the mean ET 0 for the summer and autumn showed increasing trends. The temporal variation of the mean ET 0 in winter showed a similar trend to that in the spring. Regarding the three topographic regions (MR, HR, and VR), the mean ET 0 values for the three regions exhibited similar temporal variations in all seasons except for the summer. In summer, the mean ET 0 of the MR showed a slightly increasing trend (Figure 2c). The highest value of ET 0 was at the Xiangning station (1467 mm), and the lowest value was at the Wutaishan station (495 mm). The ET 0 values in the southern and eastern areas of the CLP were always higher than those in the western part of the CLP (Figure 3), which may be related to the different spatial distributions of the climatic conditions (e.g., T, WS, RH, and SD) and the complicated topography in the CLP. In spring, the maximum mean ET 0 value occurred at the Xiangning station (388 mm) and the minimum mean ET 0 value was observed at the Wutaishan station (100 mm). The mean ET0 values of the VR were higher than those of the MR and HR. In summer, the Xiangning station (578 mm) displayed the maximum average ET 0 value, and the Wutaishan station (219 mm) displayed the minimum average ET 0 value. The average ET 0 values of the VR were higher than those of the HR and MR. As for the autumn, there was a minimum average ET 0 value observed at the Wutaishan station (133 mm), and a maximum average ET 0 value occurred at the Xiangning station (347 mm). In winter, the Xianging station (155 mm), located in the southeastern part of the Loess Plateau, displayed the maximum average ET 0 value, and a minimum average ET 0 value occurred at the Wutaishan station (44 mm). The average ET 0 values of the MR were lower than those of the HR and VR.

Spatial Distribution of Et 0 in Different Terrain Regions and Seasons
Sustainability 2019, 11, x FOR PEER REVIEW 7 of 14 spatial distributions of the climatic conditions (e.g., T, WS, RH, and SD) and the complicated topography in the CLP. In spring, the maximum mean ET0 value occurred at the Xiangning station (388 mm) and the minimum mean ET0 value was observed at the Wutaishan station (100 mm). The mean ET0 values of the VR were higher than those of the MR and HR. In summer, the Xiangning station (578 mm) displayed the maximum average ET0 value, and the Wutaishan station (219 mm) displayed the minimum average ET0 value. The average ET0 values of the VR were higher than those of the HR and MR. As for the autumn, there was a minimum average ET0 value observed at the Wutaishan station (133 mm), and a maximum average ET0 value occurred at the Xiangning station (347 mm). In winter, the Xianging station (155 mm), located in the southeastern part of the Loess Plateau, displayed the maximum average ET0 value, and a minimum average ET0 value occurred at the Wutaishan station (44 mm). The average ET0 values of the MR were lower than those of the HR and VR. In addition, the ET0 of the VR had higher values throughout the entire year, especially in spring and summer. The VR region is an important food-growing area for the CLP, and the higher ET0 values in summer and autumn may result in the regional water shortages being more serious. Therefore, water management should pay more attention to the higher ET0 values of the VR during the spring and summer.  In addition, the ET 0 of the VR had higher values throughout the entire year, especially in spring and summer. The VR region is an important food-growing area for the CLP, and the higher ET 0 values in summer and autumn may result in the regional water shortages being more serious. Therefore, water management should pay more attention to the higher ET 0 values of the VR during the spring and summer.

Spatial Variation of the Trends of Annual and Seasonal Et 0 Series-MK Tests
Nonparametric Mann-Kendall tests were performed for 70 meteorological stations to analyze trends of annual and seasonal ET 0 series in the CLP during 1960-2017 at the 5% significance level (Figure 4). The annual mean ET 0 series of the majority of stations displayed an increasing trend during 1960-2017. Approximately 93% of stations exhibited a statistically significant increasing trend in the annual mean ET 0 series. The annual ET 0 of the HR regions had a stronger significant increasing trend than that of the VR and MR regions. In the HR region, the most significant increasing trend appeared at the Wuzhong station (6.47 mm/a), and the slightly increasing trend of the annual mean ET 0 was observed at the Xintai (1.91 mm/a) and Tongxin (1.61 mm/a) stations. In the VR region, the most significant increasing trend occurred at the Xiangning station (6.39 mm/a), and the increasing trend of the Mengjin station (1.46 mm/a) was not obvious. In the MR region, the annual mean ET 0 series had overall significant increasing trends, except for Wushaoling station. However, the increase trend of the Wushaoling station (0.13 mm/a) was not obvious.   Figure 5 shows the wavelet transform coherence (WTC) results for the ET0 and the large-scale circulation patterns (for example, NAO, Atlantic Multidecadal Oscillation (AMO), Pacific-North American teleconnection (PNA), etc.). For ET0 (Figure 5a), the WTC with NAO had very good In spring, the majority of stations showed statistically significant increasing trends, with the largest increase at the Xiangning station (6.57 mm/a), and regional magnitudes declined from the VR regions to the MR regions. In summer, 59 stations showed increasing trends, while decreasing trends occurred in the majority of stations in the MR and HR. In autumn, more than 90% of the total stations showed an increasing trend during 1960-2017, and only six stations (WuQiaoLing, SanMenXia, MengJin, XiangNing, YongJi, YuLin) had statistically significant decreases. In winter, the ET 0 exhibited increasing trends at 69 stations, mostly in the MR and HR regions. A decreasing trend was observed at Hequ station, which was located on the eastern edge of the Loess Plateau.

Wavelet Transform Coherence of Et0 and the Large-Scale Circulation Patterns
Overall, on a temporal scale, an increasing trend existed in the CLP throughout almost the entire year, especially in spring. On a spatial scale, the seasonal ET 0 in the majority of the CLP exhibited increasing trends, and decreasing trends mainly occurred in the HR and VR regions.  Figure 5c shows that in the cone of influence (COI) region, the resonance period of ET 0 and El Niño-Southern Oscillation (ENSO) passing the 95% red noise test was 3-6 a; there was a significant resonance relationship at this frequency band, and the oscillation cohesion was the strongest. The coherence phase angle between ET0 and ENSO was −60 • , showing the ENSO lags behind the ET 0 by around 0.5-1 year band between 1978 and 1984. This phase change was relatively stable. There was significant coherence power and a consistent phase angle for the period 1983-1993 associated with signals in the 2-6 years band, and in this period the phase angle was 360 • , which indicated the ET 0 lags about 2-6 a behind PNA (Figure 5d). Similar to the PNA, a significant coherence power and consistent phase appeared in the period 1985-1994 associated with signals in the 2-5 years band and the coherence phase angle of ET 0 and AMO was −45 • , which means that AMO lags behind ET 0 by around 0.25-0.63 a. The ET0 and WPI showed significant resonance relationships between 6 and 8 years in 1968-1976, reflecting that the ET 0 was ahead of WPI; another significant resonance relationship was 7-10 a, and the coherence phase angle of ET0 and WPI was 10 • , indicating that ET 0 lags behind WPI by around 0.19-0.28 a (Figure 5f). Sustainability 2019, 11, x FOR PEER REVIEW 10 of 14 lags behind ET0 by around 0.25-0.63 a. The ET0 and WPI showed significant resonance relationships between 6 and 8 years in 1968-1976, reflecting that the ET0 was ahead of WPI; another significant resonance relationship was 7-10 a, and the coherence phase angle of ET0 and WPI was 10°, indicating that ET0 lags behind WPI by around 0.19-0.28 a (Figure 5f).  Table 2 shows the correlation coefficients between the ET0 and the meteorological parameters. In the CLP, the regional RH variation had a significant association with the variation of the ET0 (i.e., the correlation coefficient is 0.76). In addition, there is an obvious correlation between the variation of regional ET0 and Tmax/T. Regarding the three geomorphic regions, higher correlation coefficients between the regional ET0 variation and Tmax variation were observed at the MR (0.77) and VR (0.84) regions. For the HR region, the regional ET0 displayed the most remarkable (i.e., passed the 0.1 level test) association with the variation of RH, followed by the variations of Tmax and T. According to the results in Table 2, the main meteorological controlling parameters for the regional variation of ET0 in the CLP were RH, T, and Tmax, which indicated that regional ET0 is more sensitive to the variations in regional relative humidity and air temperature.

Relationships between Regional Et0 and Climatic Factors
However, the influences of meteorological controlling parameters on regional ET0 variations have spatial variation among the three topographic regions. There is a remarkable correlation  Table 2 shows the correlation coefficients between the ET 0 and the meteorological parameters. In the CLP, the regional RH variation had a significant association with the variation of the ET 0 (i.e., the correlation coefficient is 0.76). In addition, there is an obvious correlation between the variation of regional ET 0 and Tmax/T. Regarding the three geomorphic regions, higher correlation coefficients between the regional ET 0 variation and Tmax variation were observed at the MR (0.77) and VR (0.84) regions. For the HR region, the regional ET 0 displayed the most remarkable (i.e., passed the 0.1 level test) association with the variation of RH, followed by the variations of Tmax and T. According to the results in Table 2, the main meteorological controlling parameters for the regional variation of ET 0 in the CLP were RH, T, and Tmax, which indicated that regional ET 0 is more sensitive to the variations in regional relative humidity and air temperature.

Relationships between Regional Et0 and Climatic Factors
However, the influences of meteorological controlling parameters on regional ET 0 variations have spatial variation among the three topographic regions. There is a remarkable correlation relationship between the variation of regional T and the regional ET 0 , and variations in regional RH and T may cause regional ET 0 variability. Besides, the WS plays an important role in regional ET 0 variation, whereas the SD has no significant association with the variation of regional ET 0 . To assess the contributions of the main meteorological controlling parameters to the regional ET 0 variation in the entire CLP, the contribution rates were analyzed using the method contribution rate analysis. The results (Table 3) show the T variation had the highest contribution rate (42%). The variations of RH, WS, and SD made smaller contributions to the ET 0 in the CLP. As for the MR region, the largest contribution rate for the regional ET 0 variation was T variation (40%), followed by variations of Tmax (26%), Tmin (26%), SD (4%), RH (4%), and WS (1%). As for the VR region, the contribution rate pattern had a similar order to that of the MR region. As for the ET 0 variation in the HR region, the variation of Tmax displayed the highest contribution rate (39%). The contribution rates of the regional WS, SD, and RH were also relatively low. In short, we can conclude that the variation of regional ET 0 in the CLP was mainly controlled by regional air temperature variations, while the influences of the regional WS, SD, and RH on regional ET 0 variation were relatively weak.

Comparison with Other Study Areas
Recently, some researchers have paid increasing attention to assessing regional ET 0 variations. Due to the different environments in various regions, temporal and spatial variations of the ET 0 in different regions have been proven to have significant differences (Table 4). In a study of the ET 0 for the whole of China, researchers found that the overall ET 0 variations in China showed a downward trend [45]. This same variation was also found in the northwestern part of China [28]. However, in the southwestern part of China [46], the Wei River basin [47], and the Yangtze River [48], the ET 0 had increasing trends. The trend of the CLP is the opposite to the trend of the whole of China. Comparing the ET 0 trend in the three topographic regions (MR, HR, and VR), the trend of the ET 0 in the VR (1.49 mm/a) was larger than in the MR (1.28 mm/a) and HR (1.09 mm/a).

Conclusions
During the period of 1960-2017, the mean annual ET 0 of the CLP was 1074 mm, and the maximum and minimum values were 1181 mm and 987 mm, respectively. The trends of ET 0 in three different topographies were generally the same, however, the ET 0 values differed. The ET 0 in the HR was the highest of the three regions, followed by the VR, and the ET 0 in the MR was the lowest. The magnitude of the ET 0 varied in different seasons; in spring, summer, and autumn, the magnitudes of the ET 0 were consistent with the annual ET 0 variation, which was the largest in the HR, the second-largest in the VR, and the smallest in the MR. In winter, the ET 0 values of the HR and MR were very similar, while the ET 0 of the VR was the smallest.
The spatial distribution of the annual ET 0 in the study area varied in different topographies. The ET 0 value in the MR decreased from west to east, and the distribution of the ET 0 in the HR was consistent with the MR. The larger ET 0 values in the VR were mainly concentrated in the eastern CLP, and the rest were mostly distributed in Inner Mongolia. In the four seasons of comparative analysis, the larger values of the MR were distributed in the western region, and the ET 0 distributions in the HR and VR differed. The annual mean ET 0 series of the majority of stations displayed increasing trends during 1960-2017. Approximately 93% of stations exhibited a statistically significant increasing trend for the annual mean ET 0 series. The annual ET 0 of the HR regions had a stronger significant increasing trend than that of the VR and MR regions.
The wavelet transform coherence results of ET 0 and large-scale circulation factors showed that ET 0 had significant correlation with ENSO, PNA, and AMO. The resonance period of ET 0 and ENSO was 3-6 a, mainly in 1976-1985. The resonance period of ET 0 and PNA was mainly in 1983-1993; the average coherence phase angle was close to 360 • , indicating that the ET 0 lags behind the PNA by about 2-6 a; ET 0 had a very good positive correlation with AMO, and its significant resonance period was 2-5 a.
The RH, T, and Tmax are important meteorological parameters for the variation of control of regional ET 0 . In the CLP, RH is the main influencing parameter of the ET 0 variation; however, in the MR and VR, the correlation coefficients of Tmax were the largest (0.77 and 0.84, respectively). The correlation coefficient of the RH in the HR was 0.71, which is the main controlling parameter. The temperature variation had the highest contribution rate (42%) for the regional ET 0 variation in the entire CLP. As an important controlling factor for the regional energy balance, temporal-spatial variation of evaporation usually resulted in the regional water cycle rate increasing and further affected the regional water-heat balance and ecological security. We should pay more attention to the variation of evaporation under climate change, especially temperature change.
Author Contributions: Conceptualization, C.S.; methodology, Z.Z. and Y.W.; resources, C.S.; writing, original draft preparation, C.S. and Z.Z.; writing, review and editing, C.S. and W.C.; supervision, W.C.; funding acquisition, C.S. All authors have read and agreed to the published version of the manuscript.