Simpliﬁed Priestley–Taylor Model to Estimate Land-Surface Latent Heat of Evapotranspiration from Incident Shortwave Radiation, Satellite Vegetation Index, and Air Relative Humidity

: An operational and accurate model for estimating global or regional terrestrial latent heat of evapotranspiration (ET) across different land-cover types from satellite data is crucial. Here, a simpliﬁed Priestley–Taylor (SPT) model was developed without surface net radiation (Rn) by combining incident shortwave radiation (Rs), satellite vegetation index, and air relative humidity (RH). Ground-measured ET for 2000–2009 collected by 100 global FLUXNET eddy covariance (EC) sites was used to calibrate and evaluate the SPT model. A series of cross-validations demonstrated the reasonable performance of the SPT model to estimate seasonal and spatial ET variability. The coefﬁcients of determination (R 2 ) of the estimated versus observed daily (monthly) ET ranged from 0.42 (0.58) ( p < 0.01) at shrubland (SHR) ﬂux sites to 0.81 (0.86) ( p < 0.01) at evergreen broadleaf forest (EBF) ﬂux sites. The SPT model was applied to estimate agricultural ET at high spatial resolution (16 m) from Chinese Gaofen (GF)-1 data and monitor long-term (1982–2018) ET variations in the Three-River Headwaters Region (TRHR) of mainland China using the Global LAnd-Surface Satellite (GLASS) normalized difference vegetation index (NDVI) product. The proposed SPT model without Rn provides an alternative model for estimating regional terrestrial ET across different land-cover types.


Introduction
Latent heat of evapotranspiration (ET) observations and simulations are crucial in monitoring the energy and water cycle among different land-surface ecosystems [1][2][3][4]. Especially in typical eco-environmentally vulnerable areas, ET is affected by serious surface dry status resulting in a food and water resources crisis as well as eco-environment deterioration [5][6][7][8]. It is tough to accurately simulate regional ET because of the heterogeneous surface and complicated biophysical properties [9]. Moreover, at large scale, sparse ground observations (e.g., eddy covariance flux towers) limit accurate estimation of spatiotemporal ET distribution.
Satellite remote sensing is the most important approach for acquiring regional water fluxes because it has improved regional-scale estimates of land-surface parameters linked to ET (e.g., surface net radiation, Rn; leaf area index, LAI; normalized difference vegetation index, NDVI; soil moisture, SM; land-cover types and albedo) [10][11][12]. Over the past 40 years, many satellite-derived ET approaches have been developed to simulate regional water fluxes [13][14][15][16][17][18][19]. Generally, these ET approaches include: (1) statistical The observed data contain half-hourly ET, sensible heat flux (H), Rs, Rn, air temperature (Ta), Tmin, RH and wind speed (WS). The daily and monthly variables were also linearly aggregated from the half-hour means. When the amount of missing daily data was more than 25% of the observations, the daily values were set as missing. The fixed Bowen ratio (BR) method designed by Twine et al. [37] was used to revise ET due to the unclosed energy budget.
where ET is the revised evapotranspiration; ETori and Hori are the original ET and the H derived from EC observations, respectively; and RBR refers to the energy closure ratio.
To develop the ET model, the MODerate resolution Imaging Spectroradiometer (MODIS) NDVI products were used and they have a spatial resolution of 1-km and temporal resolution of 8-day. The eight-day NDVI averages were interpolated to daily values based on the linear interpolation method. The observed data contain half-hourly ET, sensible heat flux (H), Rs, Rn, air temperature (Ta), Tmin, RH and wind speed (WS). The daily and monthly variables were also linearly aggregated from the half-hour means. When the amount of missing daily data was more than 25% of the observations, the daily values were set as missing. The fixed Bowen ratio (BR) method designed by Twine et al. [37] was used to revise ET due to the unclosed energy budget.
where ET is the revised evapotranspiration; ET ori and H ori are the original ET and the H derived from EC observations, respectively; and R BR refers to the energy closure ratio.
To develop the ET model, the MODerate resolution Imaging Spectroradiometer (MODIS) NDVI products were used and they have a spatial resolution of 1-km and temporal resolution of 8-day. The eight-day NDVI averages were interpolated to daily values based on the linear interpolation method.

Case Study I
To illustrate the performance of the SPT model for estimating agricultural ET at a high spatial resolution from Chinese GF-1 data, the Huailai agricultural region of China was selected as a study area ( Figure 2). The Huailai agricultural region (~7.8 km 2 ) is situated between 40.3425 • N-40.3621 • N and 115.7774 • E-115.8188 • E and located in the HR Basin in northern China. The land cover is dominated by agriculture and presents small portions of other land-cover types (e.g., forest, grass, water, and barren land). Among the crops grown, the summer crop is mainly dominated by maize. Two flux tower sites, namely EC1 site (40.3491 • N, 115.7880 • E) and EC2 site (40.3574 • N, 115.7923 • E), were in the northwestern parts of this study area. The main plant function type for the two sites is continuous maize land. EC data measured between 2013 and 2017 also include half-hourly ET, H, Rs, G, Rn, Ta, Tmin, RH, VPD (vapor pressure deficit) and WS. The same procedure as for the EC data for model development was used to process these data.
was selected as a study area ( Figure 2). The Huailai agricultural region (~7.8 km 2 ) is situated between 40.3425° N-40.3621° N and 115.7774° E-115.8188° E and located in the HR Basin in northern China. The land cover is dominated by agriculture and presents small portions of other land-cover types (e.g., forest, grass, water, and barren land). Among the crops grown, the summer crop is mainly dominated by maize. Two flux tower sites, namely EC1 site (40.3491° N, 115.7880° E) and EC2 site (40.3574° N, 115.7923° E), were in the northwestern parts of this study area. The main plant function type for the two sites is continuous maize land. EC data measured between 2013 and 2017 also include halfhourly ET, H, Rs, G, Rn, Ta, Tmin, RH, VPD (vapor pressure deficit) and WS. The same procedure as for the EC data for model development was used to process these data. Fifty-five Chinese GF-1 wide field view (WFV) images for 2013 to 2017 were used to calculate NDVI for ET estimation. The Chinese GF-1 satellite, launched on 26th April 2013, has produced 16 m images per 4 days with a continuous swath of approximately 800 km [38]. The GF-1 data have four bands. Band 1 is a blue band with a spectral range of 0.45-0.52 µm. Band 2 is a green band with a spectral range of 0.52-0.59 µm. Band 3 is a red band with a spectral range of 0.63-0.69 µm. Band 4 is a near-infrared band with a spectral range of 0.77-0.89 µm. The raw DN (digital number) values were released by the Chinese GF-1 data and ENVI software was used to conduct radiance calibration, atmospheric correction, and geometric correction. Fifty-five Chinese GF-1 wide field view (WFV) images for 2013 to 2017 were used to calculate NDVI for ET estimation. The Chinese GF-1 satellite, launched on 26th April 2013, has produced 16 m images per 4 days with a continuous swath of approximately 800 km [38]. The GF-1 data have four bands. Band 1 is a blue band with a spectral range of 0.45-0.52 µm. Band 2 is a green band with a spectral range of 0.52-0.59 µm. Band 3 is a red band with a spectral range of 0.63-0.69 µm. Band 4 is a near-infrared band with a spectral range of 0.77-0.89 µm. The raw DN (digital number) values were released by the Chinese GF-1 data and ENVI software was used to conduct radiance calibration, atmospheric correction, and geometric correction.

Case Study II
The TRHR of China was also selected as another study area to assess the performance of the SPT model for monitoring long-term (1982-2018) ET variations using the Global LAnd-Surface Satellite (GLASS) NDVI product from Advanced Very High-Resolution Radiometer (AVHRR) data [39]. The TRHR (31.63 • N-36.33 • N, 89.52 • E-102.23 • E) is situated in southern parts of Qinghai Province (QH) of China and belongs to the hinterland of the Tibetan Plateau (TP) (Figure 2). TRHR is a famous region of the Lantsang, Yellow, and Yangtze rivers. Grassland (mainly alpine meadow) is the dominant ecosystem type of the TRHR. One grass flux tower site (Three-River Headwaters, TRH) (34.35 • N, 100.55 • E) was situated in the eastern parts of this study area. EC data measured between 2016 and 2017 also included half-hourly ET, H, Rs, G, Rn, Ta, Tmin, RH, VPD, and WS. The same procedure as for the EC data for model development was also used to process these variables.
To simulate the long-term evapotranspiration, the 8-day GLASS NDVI product with a spatial resolution of 0.05 degrees from 1982-2018 were used [39]. The 8-day NDVI data were linearly interpolated into daily values. The China Meteorological Forcing Dataset (CMFD) for daily Ta, Tmin, Rs, and RH for the period 1982-2018 were also used [40][41][42]. These meteorological variables (Ta, Tmin, Rs, and RH) were resampled into 0.05 • from a 0.1 • spatial resolution using linear interpolation.

Simplified Priestley-Taylor Model
The ET model was developed based on the simplified PT (SPT) model in which Rs was chosen to substitute available energy (Rn-G) since Rn can be reliably derived from Rs and other environmental variables [36] and G (soil heat flux) can be calculated using Rn [8,12,27]. In the SPT model, the ET is calculated as: where ∂ refers to PT parameter under wet surface conditions (1.26); ∆ refers to the slope for saturated vapor pressure (SVP) curve; γ stands for psychrometric constant. f(c) is the environmental constraint that is determined by meteorological and satellite variables [10,22]. It is tough to accurately calculate f(c) because of the large uncertainties in the characterization of the dynamics of the ET process using limited satellite and meteorological variables [4,43]. Because temperature and water are the key environmental constraints for controlling the evaporation fraction (EF), i.e., the ratio of ET to available energy, a linear combination of T min , NDVI, and RH was introduced to estimate f(c).
where α i (i = 0, 1, 2, 3) refers to the calibrated empirical coefficient (dimensionless) and varies among different land-cover types. In this model, T min represents the effects of temperature on EF, NDVI characterizes the vegetation transpiration influenced by temperature and water, and RH takes into account the contribution of SM on EF for a short period. The merits of the SPT model over many other physical evapotranspiration approaches are that it (1) improves operability for estimating long-term ET since it only depends on Rs, T min , NDVI and RH; (2) uses a linear combination of T min , NDVI and RH to avoid aerodynamic and surface resistance and to reduce biases of the required forcing datasets [1,27]; and (3) improves the accuracy of estimating ET by taking into account the differences in the calibrated coefficients of Equation (4) for a variety of land-cover types.

Cross-Validation
The empirical coefficients of Equation (4) were calibrated using a linear regression based on the measured data (EC observations and NDVI data) for 100 flux towers. To evaluate the estimated ET accuracy, the performance of the SPT model was assessed based on a ten-fold cross-validation approach that randomly divided the data into ten groups with about equal numbers for samples [44]. The simulated ET was independently evaluated for every of the ten groups using the calibrated empirical coefficients of Equation (4) based on the data of the remaining nine groups. The bias, root-mean-square error (RMSE), coefficients of determination (R 2 ), and p values of the simulated and observed ET were also used to illustrate the performance of the SPT model [45].
where N is the number of samples, S i is the simulated ET for sample i, M i is the groundobserved ET for sample i, S is the averaged S and M is the averaged M.
To assess the performance of the SPT model to simulate the spatial and temporal variations in ET, the ET variability was validated, including (1) seasonal variation and, (2) spatial (among-site) variation, and (3) annual anomalies, using estimated and groundmeasured data. The seasonal (daily and monthly) estimated and ground-measured ET was first validated. The spatial (among-site) variation was then validated by comparing the average ET of the estimated and ground-measured ET at every site for the whole period. Finally, the annual ET anomaly for each site was calculated by removing the multiyear average from the annual ET values.
To facilitate comparisons with other ET models, the dataset was randomly stratified into two different groups for about equal numbers for samples, and entire sites were assigned to each group [44]. Daily ET was validated using datasets from one group using the calibrated parameters of Equation (4) based on data from another group.

Comparison to the PT-JPL Model
The revised Priestley-Taylor (PT) model from the Jet Propulsion Laboratory (PT-JPL) was developed based on the original PT equation by introducing ecophysiological constraints (NDVI and LAI) and atmospheric constraints (RH and VPD) to estimate terrestrial ET [1,10]. The validation for 16 global EC sites illustrates that the PT-JPL model accounts for 90% of ET variations, and comparisons with multiple models indicate that this model is the most promising ET method model [30,46]. An obvious difference between the SPT model and PT-JPL is that the SPT model uses ground measurements to calibrate the coefficients of the ecophysiological constraints to develop f(c) for different land-cover types.

Model Parameterization
Long-term (2000-2009) ground measurements and the corresponding NDVI datasets for different land-cover types were analyzed at 100 flux tower sites to detect the driving variables for determining terrestrial ET. Table 1 illustrates the correlation coefficients (r) between ET and Rs, Tmin, RH, and MODIS NDVI for different land-cover types. For all land-cover types, Rs presented the highest R 2 with ET (0.20 ≤ R 2 ≤ 0.56) and Tmin presented the second highest (0.17 ≤ R 2 ≤ 0.53), demonstrating that both Rs and Tmin are the dominant variables driving terrestrial ET. The NDVI was also highly correlated with ET, with R 2 varying from 0.16 to 0.50 for most land-cover types (except for WET), which illustrates that the NDVI is an important variable controlling the variation of terrestrial ET. The NDVI characterizes vegetation amounts and reflects vegetation transpiration by coupling CO 2 assimilation with vegetation photosynthesis [47]. The R 2 values between RH and ET were relatively low (less than 0.11) for all land-cover types due to the seasonal variability of RH. To parameterize f(c) of the SPT model, correlation coefficients (r) between groundmeasured f(c) and the corresponding Tmin, RH, and MODIS NDVI from all EC sites covering different terrestrial biomes were analyzed. Table 2 summarizes the f(c) correlation coefficients with Tmin, RH, and MODIS NDVI for different land-cover types. For most land-cover types, ground-measured f(c) illustrated the highest correlations with the NDVI (0.12 ≤ R 2 ≤ 0.46), Tmin (0.11 ≤ R 2 ≤ 0.32), and RH (0.10 ≤ R 2 ≤ 0.30). However, for the WET sites, f(c) had no correlation with the NDVI because the surface water supply was saturated and water evaporation mainly occurred in wetlands. Similarly, there were no obvious correlations between f(c) and Tmin for the SHR, CRO, and GRA sites, which might partially be attributed to the lower sensitivity of surface conductivity to Tmin in crops, grasses and shrubs [4,48].

Model Validation and Comparison
The SPT model was validated using ten-fold cross-validation. Validations of the simulated and observed daily (monthly) ET at all sites illustrated the good performance of the SPT model for estimating seasonal ET variability ( Figure 3). For different land-cover types, the performance of the SPT model differs greatly. The R 2 of the estimated versus observed daily (monthly) ET ranged from 0.42 (0.58) (p < 0.01) at all SHR tower sites to 0.81 (0.86) (p < 0.01) at all EBF tower sites. The RMSE varied from 15.8 (11.7) W/m 2 at all EBF tower sites to 28.2 (23.9) W/m 2 at all DBF tower sites. Figure 3 also demonstrates the good ability of the SPT model for estimating the spatial variation. The RMSE of the site-averaged estimated ET varied from 6.7 W/m 2 at all EBF tower sites to 18 all CRO tower sites, and R 2 (99% confidence) varied from 0.22 at all SHR tower sites to 0.97 at all EBF tower sites. Overall, the SPT model accurately captured seasonal and spatial variations in ET. The SPT model also accurately yielded long-term variations in ET at all sites. The R 2 (99% confidence) of the estimated and observed annual ET anomalies varied from 0.20 at all SHR tower sites and 0.77 at all MIF tower sites, and the largest RMSE was less than 12 W/m 2 at all DNF tower sites, which was partially because of missing observed ET datasets caused by worse climatological conditions [18,46].    The estimated daily ET using the SPT model was also compared with those for the PT-JPL model. Table 3 shows the statistics (R 2 , RMSE and bias) of the comparison between the estimated and ground-observed daily ET from the first group of data using the second group of samples to calibrate the SPT model. For the forest EC sites, the RMSE for the estimated daily ET using the SPT model and PT-JPL model was 22.3 W/m 2 and 31.2 W/m 2 , respectively. R 2 of the estimated daily ET using the SPT model was higher than that for the PT-JPL model. For the crop, grass, and wetland sites, the SPT model had better performance than the PT-JPL model, with RMSEs less than 30 W/m 2 and R 2 values more than 0.85 (p < 0.01). For the SAW and SHR sites, the R 2 of the estimated ET using the SPT model was much higher and the RMSE was less than 25 W/m 2 compared with the PT-JPL model. Overall, the R 2 for the SPT model increased by about 0.05 (p < 0.01) and the RMSE decreased by about 9 W/m 2 at most forest flux tower sites. Similarly, the R 2 for the SPT model increased by about 0.10 (p < 0.01) and the RMSE decreased by about 3 W/m 2 at the CRO, GRA, WET, SAW, and SHR sites. This improvement over PT-JPL and other PT models is to calibrate the SPT model using ET ground measurements, with which the ground-measured ET was also used for validation. It is clear that the SPT model has better performance for estimating ET than that of the PT-JPL model. Therefore, compared with other PT models, the proposed SPT model has provided a better fit to EC observations. Table 3. Summary of a series of statistics (bias, RMSE, and R 2 ) of the comparison between the ground-measured and estimated average daily ET using the two ET models of the first group. The second group samples were considered to be training datasets to calibrate coefficients of the proposed ET model. CRO: cropland; GRA: grassland; DBF: deciduous broadleaf forest; MIF: mixed forest; DNF: deciduous needleleaf forest; SAW: savanna; EBF: evergreen broadleaf forest; SHR: shrubland; ENF: evergreen needleleaf forest; WET: wetland. The ET validation for each land-cover types agreed, in general, with results documented in the literature [30,34,49]. Ershadi et al. [30] evaluated four ET models using ground-measured ET from 20 FLUXNET tower sites and found that the best model is PT-JPL model with the RMSE of 66 W/m 2 . Vinukollu et al. [34] assessed three processbased ET models from 12 flux tower sites across the US and reported that these models yielded monthly ET with mean correlation of 0.54-0.57. Yebra et al. [49] compared estimates of ET produced with six different vegetation indices from the MODIS and three contrasting estimation approaches using ground-measured ET at 16 FLUXNET sites and documented that the PM-Gs (surface conductance) approach provided the lowest RMSE (average RMSE = 38 W/m 2 ), and highest R 2 (average R 2 = 0.72). The validation accuracy of the estimated ET using the SPT model is higher than the above reported results.

Land-Cover Types
To estimate the regional terrestrial ET, the parameters of Equation (4) were recalibrated based on ground-observed meteorological datasets and the MODIS NDVI datasets from all 100 EC flux tower sites. Table 4 demonstrates the Equation (4) coefficients for different land-cover types using a linear regression from MODIS NDVI, ground-observed Ta, Tmin, and RH. The SPT model is sufficiently representative to estimate regional ET due to the different terrestrial biomes and their locations of 100 EC sites. The calibrated SPT model was applied to estimate agricultural ET of the Huailai region of China at high spatial resolution derived from GF-1 data during 2013-2017. Figure 4 illustrates the good correspondence between the estimated and ground-observed daily ET at the EC1 and EC2 sites. RMSE of estimated versus observed ET was 23.6 W/m 2 and 24.4 W/m 2 , while that for the bias was −11.6 W/m 2 and −1.2 W/m 2 . In addition, the R 2 (99% confidence) was 0.84 and 0.81 at the EC1 and EC2 sites, respectively. These good results indicate that the ET estimation of the SPT model is reliable in the Huailai agricultural region of China.   , the estimated ET ranged from 0 to 80 W/m 2 , and the highes curred in the wetland regions of the top left part of the study area, which w with the spatial patterns in the spring. Overall, the estimated ET tended tow values during the period of summer because high vegetation coverage occu region, while low ET values occurred in late fall after crops were harvested grass and forest declined. Few previous studies have focused on regional L at relatively high spatial resolutions using GF-1 data [46]. The SPT model s acquires finer spatial resolution ET using GF-1 NDVI data.  study area, which was consistent with the spatial patterns in the spring. Overall, the estimated ET tended towards high values during the period of summer because high vegetation coverage occurred in this region, while low ET values occurred in late fall after crops were harvested and other grass and forest declined. Few previous studies have focused on regional LE estimation at relatively high spatial resolutions using GF-1 data [46]. The SPT model successfully acquires finer spatial resolution ET using GF-1 NDVI data.
with the spatial patterns in the spring. Overall, the estimated ET tended towards high values during the period of summer because high vegetation coverage occurred in this region, while low ET values occurred in late fall after crops were harvested and other grass and forest declined. Few previous studies have focused on regional LE estimation at relatively high spatial resolutions using GF-1 data [46]. The SPT model successfully acquires finer spatial resolution ET using GF-1 NDVI data. The SPT model driven by CMFD meteorological data and GLASS NDVI products was also applied to generate the daily terrestrial ET of TRHR of China at a 0.05° spatial resolution between 1982 and 2018. This daily ET product was also validated based on ground-observed EC data collected from one TRH grass site. Figure 6 demonstrates the statistical comparison between the estimated and ground-observed ET at the TRH site. The R 2 for the TRH site was approximately 0.67 at the 99% level of confidence, the bias The SPT model driven by CMFD meteorological data and GLASS NDVI products was also applied to generate the daily terrestrial ET of TRHR of China at a 0.05 • spatial resolution between 1982 and 2018. This daily ET product was also validated based on ground-observed EC data collected from one TRH grass site. Figure 6 demonstrates the statistical comparison between the estimated and ground-observed ET at the TRH site. The R 2 for the TRH site was approximately 0.67 at the 99% level of confidence, the bias was 2.0 W/m 2 and the RMSE was 14.9 W/m 2 . The slightly positive bias may be caused by the limitations of the SPT model and the unclosed energy balance problem in EC data. However, Seguin et al. [50] pointed out that the ET estimation accuracy for application can be accepted as less than 50 W/m 2 . The accuracy of the results meets this basic requirement. However, Seguin et al. [50] pointed out that the ET estimation accuracy for application can be accepted as less than 50 W/m 2 . The accuracy of the results meets this basic requirement.  Figure 7 presents the spatial pattern of the multiyear (between 1982 and 2018) averaged annual ET in the TRHR. A large ET occurred in the eastern areas of the TRHR, while a small ET occurred across the western regions in the TRHR. The average annual ET of TRHR was 200 mm/year. Overall, the estimated annual ET anomalies of TRHR showed a significantly increasing trend of 0.3 W/m 2 per decade (p < 0.05) from 1982-2018, and they coincided with widespread regional grass greening [51] (Figure 8). Regionally, the ET in the western areas of the TRHR accounted for more than 75% of the variation in the whole regional ET, and the positive annual ET in this region was mainly attributed to the regional vegetation (mainly grass) greening caused by global warning (0.065 °C/decade). The pattern of rising ET matched by an increasing NDVI indicated that increases in the NDVI (vegetation greening) represented the main contributor to the increasing ET trend over the western areas of the TRHR between 1982 and 2018. In contrast, the large-scale droughts characterized by the decreasing RH led to the regional ET decline in the eastern areas of the TRHR (Figure 9). More than 70% of the eastern region of the TRHR had positive correlations between annual ET and RH. The strong consistency of the patterns in terrestrial ET and RH trends illustrated that large-scale droughts were the main mechanism leading to the decreasing ET trend in the eastern areas of the TRHR from 1982-2018.  Figure 7 presents the spatial pattern of the multiyear (between 1982 and 2018) averaged annual ET in the TRHR. A large ET occurred in the eastern areas of the TRHR, while a small ET occurred across the western regions in the TRHR. The average annual ET of TRHR was 200 mm/year. Overall, the estimated annual ET anomalies of TRHR showed a significantly increasing trend of 0.3 W/m 2 per decade (p < 0.05) from 1982-2018, and they coincided with widespread regional grass greening [51] (Figure 8). Regionally, the ET in the western areas of the TRHR accounted for more than 75% of the variation in the whole regional ET, and the positive annual ET in this region was mainly attributed to the regional vegetation (mainly grass) greening caused by global warning (0.065 • C/decade). The pattern of rising ET matched by an increasing NDVI indicated that increases in the NDVI (vegetation greening) represented the main contributor to the increasing ET trend over the western areas of the TRHR between 1982 and 2018. In contrast, the large-scale droughts characterized by the decreasing RH led to the regional ET decline in the eastern areas of the TRHR (Figure 9). More than 70% of the eastern region of the TRHR had positive correlations between annual ET and RH. The strong consistency of the patterns in terrestrial ET and RH trends illustrated that large-scale droughts were the main mechanism leading to the decreasing ET trend in the eastern areas of the TRHR from 1982-2018. NDVI (vegetation greening) represented the main contributor to the increasing ET over the western areas of the TRHR between 1982 and 2018. In contrast, the large droughts characterized by the decreasing RH led to the regional ET decline in the e areas of the TRHR (Figure 9). More than 70% of the eastern region of the TRHR had tive correlations between annual ET and RH. The strong consistency of the patte terrestrial ET and RH trends illustrated that large-scale droughts were the main m nism leading to the decreasing ET trend in the eastern areas of the TRHR from 1982

Model Performance
Model validation at 100 global EC sites demonstrated that the proposed SPT for estimating terrestrial ET was reliable across different terrestrial biomes. Howev SPT model had large differences in estimating ET across different land-cover types the same environmental conditions. For example, the SPT model accounted for mo Figure 9. Spatial pattern of the RH trends for 1982-2018 over the TRHR. The black solid dots refer to grids with 95% confidence.

Model Performance
Model validation at 100 global EC sites demonstrated that the proposed SPT model for estimating terrestrial ET was reliable across different terrestrial biomes. However, the SPT model had large differences in estimating ET across different land-cover types under the same environmental conditions. For example, the SPT model accounted for more than 75% of the ET variability for DBF, which may be associated with the strong seasonal variations in vegetation LAI in these biomes and the sensitivity of the NDVI to red reflectance and LAI changes [27,49,52]. In contrast, the SPT model only explained 50% of the ET variability for SHR and ENF because vegetation seasonal changes were weak and the NDVI yielded a poor performance for modeling vegetation transpiration [22,52,53]. This finding is consistent with previous studies documenting that the vegetation conductance of both ENF is half that of DBF [4,53]. In the SPT model, environmental factors, including Tmin, NDVI and RH, are considered to be the important regulators for downscaling the potential ET to actual ET. However, at the CRO, GRA, and SHR sites, Table 1 shows that the contribution of Tmin to f(c) is limited because the short vegetation is mainly controlled by water availability rather than Tmin [4,26]. Similarly, NDVI has no influence on f(c) at WET sites because the wetland is saturated with water [4,19]. The SPT model explains their differences using the same environmental factors with different calibration coefficients.
Although the SPT model only uses a simple linear formulation to parameterize f(c), it has lower errors in the inputs and has a better performance than that of the PT-JPL model. Importantly, the SPT model avoids Rn estimation directly using the Rs product as an important input and considers the differences in the f(c) calculation at different land-cover types, while the PT-JPL model ignores the classifications among land-cover types [10,27].
The simulation accuracy of the proposed model is affected by the errors in EC observations, the biases of forcing data, the scaling effects of different data and the structure of the model. Generally, EC observations have a bias of approximately 10-25% [54]; however, the EC technique is considered a good method to measure ET. Furthermore, EC suffers from energy imbalance problems because it does not capture large eddies, and the annual R BR for approximately 200 global FLUXNET sites is only 0.8 [55]. Although the fixed BR method was used to correct the ET, the uncertainty from EC observations remains indistinct [37,54]. Additionally, errors in meteorological variables (Tmin, Ta and RH) will also lead to errors in ET estimation. Recent studies have presented errors for CMFD datasets compared to ground observations [40][41][42]. Thus, the biases of meteorological data will reduce the accuracy of ET estimations.
Generally, the footprints of EC observations are less than 300 m [9]. However, the spatial resolutions of the MODIS (~1 km), GLASS (~5 km) and GF-1 (~16 m) NDVI products do not match the footprints of the EC observations. The mismatch among these datasets will cause large errors in the ET estimates. The structure of the SPT model also affects the ET simulation errors because the SPT model does not include SM and only uses RH to characterize moisture constraints, which may cause large uncertainties in ET estimation in semiarid and arid regions [16,22].

Implication for Regional Water Resources Assessment
Accurately quantifying the regional terrestrial ET is important for agricultural water resource management. Especially for agricultural regions, satellites with high spatial resolution (e.g., GF-1) have provided spatially and temporally frequent terrestrial observations for estimating field ET. Compared with LST, the NDVI changes slowly under a relatively short period to characterize SM conditions. However, the NDVI cannot reflect soil evaporation for bare soil [56]. Here, RH was used as a substitute for SM to estimate soil evaporation. Bouchet [57] reported a complementary theory that surface SM is associated with the atmospheric evaporation demand. Subsequently, Fisher et al. [10] developed a SM metric (RH VPD ) to calculate soil evaporation. Similarly, Yan and Shugart [58] also found that RH is a better indicator than SM to estimate terrestrial ET, especially for late crop harvest periods. This good theoretical hypothesis illustrates the ability of the proposed SPT model to accurately estimate agricultural ET and has provided a usefully diagnostic assessment for field water consumption.
Monitoring long-term ET variations in ecologically fragile regions is critical for understanding regional hydrological dynamics and performing water resource assessments. In this study, the proposed SPT model has provided evidence for vegetation greening caused by global warming, which is important for controlling long-term ET variations in the TRHR of China. Previous studies have reported that climatic change (Ta and RH variations) has increased vegetation greening in northern high latitudes [59,60]. Similarly, the increase in ET is caused by the greening of vegetation, which is consistent with previous studies [2,61]. The results in this study also support this conclusion (Figures 8 and 9). Other mechanisms that account for increasing ET seem to be less important due to the weak correlations between ET changes and other meteorological variables [62,63]. Further studies should highlight the combinations of satellite ET and ground measurements to investigate the impacts of human activities on hydrological processes.

Conclusions
A simplified Priestley-Taylor (SPT) ET model was developed to estimate global terrestrial ET across multiple biomes by combining incident shortwave radiation, a satellite vegetation index, and air relative humidity. The model coefficients of the environmental parameters are calibrated using the satellite-derived NDVI and ground-observed Ta, Tmin, RH, and ET. The proposed model was also validated based on a ten-fold cross-validation method across different plant function types using global FLUXNET observations from 100 global EC sites.
The validation results showed that the SPT model accurately captures seasonal and spatial variations in ET. For multiple biomes, the performance of the simplified PT model differed greatly. The R 2 of the estimated versus observed daily (monthly) ET ranged from 0.42 (0.58) (p < 0.01) at all shrubland (SHR) flux tower sites to 0.81 (0.86) (p < 0.01) at all EBF flux tower sites. The RMSE of the site-averaged estimated ET varied from 6.7 W/m 2 at all EBF flux tower sites to 18.2 W/m 2 at all cropland (CRO) flux tower sites. The performance of the SPT model was good for simulating interannual ET variation, with an R 2 of approximately 0.48 (p < 0.01). When compared with the PT-JPL model, the SPT model performed better at different regional scales. The SPT model was also applied to estimate agricultural ET at high spatial resolution from Chinese GF-1 data and to monitor long-term (1982-2018) ET variations in the TRHR of China. The SPT model presented here provides a bridge between EC flux data and regional applications.