Estimates of Daily Evapotranspiration in the Source Region of the Yellow River Combining Visible/Near-Infrared and Microwave Remote Sensing

: The spatial variation of surface net radiation, soil heat ﬂux, sensible heat ﬂux, and latent heat ﬂux at different times of the day over the northern Tibetan Plateau were estimated using the Surface Energy Balance System algorithm, data from the FY-2G geostationary meteorological satellite, and microwave data from the FY-3C polar-orbiting meteorological satellite. In addition, the evaporative fraction was analyzed, and the total evapotranspiration (ET) was obtained by the effective evaporative fraction to avoid the error from accumulation. The hourly change of latent heat ﬂux presented a sound unimodal diurnal variation. The results showed the regional ET ranged between 2.0 and 4.0 mm over the Source Region of the Yellow River. The conditional expectations of surface energy components during the experimental period of the study area were statistically analyzed, and the correspondence between different surface temperatures and the effective energy distribution was examined. The effective energy distribution of the surface changed signiﬁcantly with the increase in temperature; in particular, when the surface temperature exceeded 290 K, the effective energy was mainly used for surface ET. The aim of this study was to avoid the use of surface meteorological observations that are not readily available over large areas, and the ﬁndings lay a foundation for the commercialization of land surface evapotranspiration.


Introduction
The variation of regional evapotranspiration is of significance for the land surface process and the energy exchanges in the hydrosphere, atmosphere, and biosphere [1][2][3]. Due to the development of observation methods and understanding of the turbulence process in the atmospheric boundary layer, the method of calculating land surface evapotranspiration has gradually developed from the traditional spot-measurements-based method to regional evapotranspiration estimation based on the technology of large-scale satellite remote sensing [4].
Since the mid-1970s, numerous scholars have undertaken a large amount of research on the accuracy of evapotranspiration retrieved by remote sensing [5,6]. According to different theories, remote sensing evapotranspiration models can be roughly divided into three categories. First category is based on the Penman-Monteith algorithm with an improvement of its parameters, using remote sensing data and ground observation data to estimate the actual evapotranspiration [7]. The second type is a statistical empirical model, which attempts to establish a correlation between evapotranspiration and vegetation index (VI) or surface-air temperature difference [8,9], or derives the two-dimensional R n = (1 − α)Q + L ↓ −L ↑ (1) where, α is the surface albedo, determined by MODIS observation data; Q is the total solar radiation, including direct radiation, scattered radiation, and reflected radiation [19].
L↓ is the downward long-wave radiation, which can be obtained by: where, ε a is the air effective emissivity, obtained from actual water vapor pressure and air temperature [20]; σ is the Stefan-Boltzman constant, which is 5.669 × 10 −8 Wm −2 k −4 ; T a is the measured air temperature. L↑ is the upward long-wave radiation, which can be calculated from the formula: where, ε is the land surface emissivity, obtained by the methods of Chen et al. [19]; T S is the surface skin temperature, which can be obtained using the Black Body Temperature (TBB) data observed by the Visible Infrared Spin Scan Radiometer (VISSR) on FY-2G and the MODerate resolution atmospheric TRANsmission model (MODTRAN). Specifically, the band response functions of MODTRAN and VISSR are adopted to calculate the atmospheric thermal infrared transmittance and atmospheric thermal infrared radiance. Then, brightness temperature from VISSR thermal infrared band-1 and band-2 are chosen in the surface radiation transfer model. Finally, the results can be derived using an iterative algorithm.

Net Radiation under Cloud Cover
The SRYR is located in the northeast of the Qinghai-Tibet Plateau. The development of convective clouds over the plateau in the summer afternoon limits visible and nearinfrared remote sensing. Therefore, a simplified parameterization scheme of surface radiation flux [21] was adopted under cloud-cover conditions. This scheme assumes that the long-wave radiations are approximately equal under cloud cover.
(1) Atmospheric attenuation The attenuation effect of the atmosphere on solar radiation can be presented by the transmittance τ: where, Q E is the solar radiation that reaches the surface after atmospheric attenuation; Q 0 is the total solar radiation at the top of the atmosphere.
(2) The energy received by the satellite sensor The energy received by the satellite sensor can be divided into two parts: one is the energy of solar radiation that has not been attenuated by the atmosphere, and the other is the energy that reaches the surface after atmospheric attenuation and is reflected [21].
where, α TOA is the albedo of the top cloud, α TOA Q 0 is the energy received by the sensor; α s is the surface albedo; Q 0 τ α s τ is the energy passing to the surface after atmospheric attenuation, and reflected and attenuated by the surface; (1 − τ)Q 0 is the energy reflected at the cloud top. Using Equation (5), the formula of atmospheric transmittance can be given by the following: where α TOA and α s are taken from the observations of the geostationary and polar-orbiting satellites, respectively. The hourly temporal scale transmittance can be derived from Equation (6), and the net radiation calculation formula under the cloudy condition is: It should be noted that the atmospheric transmittance obtained from the above formula is not the complete atmospheric transmittance of solar radiation, but the total equivalent extinction effect of the earth-atmosphere system on the radiation transfer path from outside the atmosphere to the ground surface. It represents the interaction between absorption and scattering effects. For satellite sensor observations that are missing or with obviously wrong values, we established the relationship between solar short-wave radiation and solar radiation before the atmospheric effect to calculate the atmospheric transmittance.

Soil Heat Flux
The formula for soil heat flux is given by: where, Γ c and Γ s represent the ratios between G 0 and R n with 100% vegetation cover and bare soil area. Vegetation coverage follows f v = NDV I−NDV I min NDV I max −NDV I min , where NDVI min is the normalized difference vegetation index of the bare soil surface, and NDVI max is the normalized difference vegetation index when the area is fully covered with vegetation [21].

Calculation Method of Sensible Heat Flux
An iterative method is adopted to calculate the sensible heat flux and the specific calculation formula is as follows: where Z is the height of the reference surface, u is the wind speed, u * is the friction wind speed, d o is the zero-plane displacement height, and Z om and Z oh denote momentum and thermal roughness lengths, respectively. ψ m and ψ h are the stability correction functions for dynamic and thermodynamic transfer, respectively; θ o and θ a are the potential temperature of the surface air and the potential temperature of the reference height, respectively; H is the sensible heat flux; L is the Obukhov length; k is the von Kármán constant; ρ is the air density, C p is the gas constant; g is the acceleration of gravity; and θ v is the virtual potential temperature near the surface. More parameters can be seen in Su [14]. The Leaf Area Index (LAI) needed for the calculation of Z oh and Z om can be derived from the method of Chen et al. [19]. In this research, the 37 GHz vertically polarized brightness temperature from the MWRI on the FY-3C polar-orbiting meteorological satellite was used to retrieve the surface temperature under cloud-cover conditions at 16:00 [22].

Evaporation Fraction
The evaporation fraction is the ratio of actual ET to effective energy (the difference between net radiation and soil heat flux) and can be presented as where EF is the evaporation fraction, LE is the latent heat flux, R n is the net radiation flux, and G is the soil heat flux. The SEBS model defines the limit cases as deficient and sufficient [23]. Since evaporation fraction is stable during the whole day, if we know the net radiation and soil heat flux at a certain time, then the latent heat flux at that time can be calculated using the latent heat of vaporization (2.49 × 10 6 Wm −2 mm −1 ).

Validation
Ground observations from the regular meteorological station and eddy covariance system from the SRYR Alpine Climate Station (102 • 08 E, 33 • 53 N, 3443 m) were used to validate the instantaneous heat fluxes that were estimated from remote sensing retrievals from 1-30 August 2019. The statistical analysis using the root mean square error is also presented in Section 4.2.
Since the soil heat flux observation was at 10 cm, we used the Averaging Soil Thermocouple Probe (TCAV) method to elicit the surface soil heat flux through soil heat storage, which was obtained from the average soil temperature [24]. The daily average soil volumetric water content of 0-10 cm at the observation site was 0.39, and the average porosity of shallow soil was 0.50.

Study Area and Observation Data
The source region of the Yellow River is located in the northeastern part of the Qinghai-Tibet Plateau (95.5 • E-103.5 • E, 31.5 • N-36.5 • N), with an area of about 1.22 × 10 5 km 2 , accounting for 16% of the total area of the Yellow River Basin. This area, mostly with an altitude of above 4000 m, is mainly covered with alpine meadows (Figure 1). The source region of the Yellow River is located in the "arid-semi-arid" transition zone between the Qinghai-Tibet Plateau and the Loess Plateau. It has a typical continental alpine climate with an average annual temperature of 2.9 • C. Precipitation is mainly concentrated in the period from May to September, accounting for more than 80% of the total, which is 560 mm per year . The annual evaporation is 920 mm, gradually decreasing from southeast to northwest [21].
Remote Sens. 2020, 12, x 5 of 15 In this research, the 37 GHz vertically polarized brightness temperature from the MWRI on the FY-3C polar-orbiting meteorological satellite was used to retrieve the surface temperature under cloud-cover conditions at 16:00 [22].

Evaporation Fraction
The evaporation fraction is the ratio of actual ET to effective energy (the difference between net radiation and soil heat flux) and can be presented as where EF is the evaporation fraction, LE is the latent heat flux, Rn is the net radiation flux, and G is the soil heat flux. The SEBS model defines the limit cases as deficient and sufficient [23]. Since evaporation fraction is stable during the whole day, if we know the net radiation and soil heat flux at a certain time, then the latent heat flux at that time can be calculated using the latent heat of vaporization (2.49 × 10 6 Wm −2 mm −1 ).

Validation
Ground observations from the regular meteorological station and eddy covariance system from the SRYR Alpine Climate Station (102°08′E, 33°53′N, 3443 m) were used to validate the instantaneous heat fluxes that were estimated from remote sensing retrievals from 1-30 August 2019. The statistical analysis using the root mean square error is also presented in Section 4.2.
Since the soil heat flux observation was at 10 cm, we used the Averaging Soil Thermocouple Probe (TCAV) method to elicit the surface soil heat flux through soil heat storage, which was obtained from the average soil temperature [24]. The daily average soil volumetric water content of 0-10 cm at the observation site was 0.39, and the average porosity of shallow soil was 0.50.

Study Area and Observation Data
The source region of the Yellow River is located in the northeastern part of the Qinghai-Tibet Plateau (95.5°E-103.5°E, 31.5°N-36.5°N), with an area of about 1.22 × 10 5 km 2 , accounting for 16% of the total area of the Yellow River Basin. This area, mostly with an altitude of above 4000 m, is mainly covered with alpine meadows (Figure 1). The source region of the Yellow River is located in the "arid-semi-arid" transition zone between the Qinghai-Tibet Plateau and the Loess Plateau. It has a typical continental alpine climate with an average annual temperature of 2.9 °C. Precipitation is mainly concentrated in the period from May to September, accounting for more than 80% of the total, which is 560 mm per year . The annual evaporation is 920 mm, gradually decreasing from southeast to northwest [21].  air temperature, air pressure, relative air humidity, surface temperature, wind speed at 10 m, precipitation, and sunshine hours. An eddy covariance system was erected at 3.02 m, with a frequency of 10 Hz, observing short-wave radiation flux and long-wave radiation flux (upward and downward), sensible heat flux, latent heat flux, soil heat flux (two layers: −10 cm and −20 cm). The system was set to output average flux every 10 min, and the results were finally averaged to obtain the hourly temporal scale data.
FY-2G and FY-3C are the geostationary meteorological and the polar-orbiting meteorological satellites, respectively, developed and operated by China. The Visible Infrared Spin Scan Radiometer carried by the FY-2G is the main earth observing instrument on the FY series satellites of the China Earth Observing System. It provides five channels of data from 1 July 2015. The resolution is 1.25 km for the visible channel, and 5 km for infrared and water vapor channels. The MicroWave Radiation Imager (MWRI) on FY-3C can be used for monitoring and research of soil moisture, snow cover, and ocean salinity. The sensor contains five frequencies-10.65, 18.7, 23.8, 36.5, and 89 GHz-and the corresponding spatial resolutions of sub-satellite points are 51 km, 85 km, 30 × 50 km, 27 × 45 km, 18 × 30 km, and 9 × 15 km. The data of the FY series satellites used in the study were all downloaded from the Fengyun Satellite Data Center (http://satellite.nsmc.org.cn/portalsite/default.aspx), including the 1-h average thermal infrared radiation brightness temperature and visible wide-band albedo, and MWRI 36.5 GHz microwave radiation brightness temperature. For specific processing steps, refer to the relevant literature [25].
In addition, the remote sensing data used in this study also included surface albedo, NDVI, LAI, and DEM. Surface albedo, NDVI, and LAI were obtained from MODIS (http://modis.gsfc.nasa.gov/). NDVI was adopted by the maximum value composite from an 8-day daily MODIS-NDVI, which eliminated the adverse effects of clouds, atmosphere, and satellite zenith angles. The DEM was obtained from the Shuttle Radar Topography Mission (SRTM) (http://srtm.csi.cgiar.org/), and its spatial resolution was 90 m. After the data or products were resampled to the nominal projection format of FY-2G, the final results were converted to the 0.078 • longitude and latitude coordinate system to complete the drawing and statistical analysis. Figure 2 shows the surface albedo, NDVI, and LAI in the SRYR retrieved by satellite remote sensing during the study period. The distribution figure shows that the surface albedo, NDVI, and LAI in the study area were highly consistent with its land use. Since the study area was dominated by alpine meadow grasslands, the surface albedo of the entire area was generally below 0.25. The minimum albedo (<0.10) occurred mainly in the west of the study area, indicating that there were large lakes, whereas the maximum value (>0.25) was in the northwest, indicating there was snow or permafrost tundra. The NDVI value in the area was between 0.45 and 0.85. Due to the distribution of lakes, certain areas also showed low values between 0.10 and 0.40. The value of LAI was between 0.9 and 2.1, indicating that although the distribution of vegetation presented certain spatial differences, there was no fundamental difference in general, and LAI had a more obvious spatial variation range than NDVI. The altitude in the region was around 3200-5200 m, and the location of the observation site was relatively flat. It should be noted that this research did not examine the surface ET of snow and lakes.

Results of Energy Flux
Because the underlying surface of the study area was relatively flat and due to the footprint effect, the surface observation was of spatial representativeness from tens of meters to several kilometers. That means the energy flux from the surface observation can be approximately taken as the surface energy flux of the entire study area. To verify the satellite remote sensing algorithm and the estimated results, the field observations were compared with the estimated results. Figure 3 shows the daily variation of the near-surface energy flux in the SRYR retrieved by satellite remote sensing on August 18th, 2019. The figure shows the hourly net radiation, soil heat flux, sensible heat flux, and latent heat flux from 11:00 to 16:00 Beijing time. The solar radiation absorbed by the ground after sunrise was greater than the expended radiation. At 11:00, the net radiation was about 400.0 Wm −2 , the soil heat flux was about 80.0 Wm −2 , the sensible heat flux was about 120.0 Wm −2 , and the latent heat flux value was about 200.0 Wm −2 ; At midday, when the solar radiation strengthened, the net radiation continued to increase at 12:00 and 13:00, and the net radiation, soil heat flux, and sensible heat flux all reached the daily maximum at 12:00 and 13:00, which were about 650.0 Wm

Results of Energy Flux
Because the underlying surface of the study area was relatively flat and due to the footprint effect, the surface observation was of spatial representativeness from tens of meters to several kilometers. That means the energy flux from the surface observation can be approximately taken as the surface energy flux of the entire study area. To verify the satellite remote sensing algorithm and the estimated results, the field observations were compared with the estimated results. Figure   From Figure 3, it can be seen that as the plateau convective cloud started to d at 15:00, the near-infrared band of the geostationary satellite was affected by cloud thus, the flux could not be retrieved in most of the area. However, the MWRI on the orbiting satellite passed through the SRYR at around 15:40. Figure 4 shows the late flux retrieved from the FY-3C MWRI and FY-2G visible/near-infrared data. The m data over the cloud cover area from FY-2G was replaced with data from FY-3C. As in the figure, in areas with no cloud cover, the spatial distribution of high and low retrieved by the two methods were basically consistent. Cloudy areas were main centrated in the north and south of the study area, and the latent heat flux was abou From Figure 3, it can be seen that as the plateau convective cloud started to develop at 15:00, the near-infrared band of the geostationary satellite was affected by cloud cover, thus, the flux could not be retrieved in most of the area. However, the MWRI on the polar-orbiting satellite passed through the SRYR at around 15:40. Figure 4 shows the latent heat flux retrieved from the FY-3C MWRI and FY-2G visible/near-infrared data. The missing data over the cloud cover area from FY-2G was replaced with data from FY-3C. As shown in the figure, in areas with no cloud cover, the spatial distribution of high and low values retrieved by the two methods were basically consistent. Cloudy areas were mainly concentrated in the north and south of the study area, and the latent heat flux was about 200.0 Wm −2 (Figure 4b). However, the latent heat flux obtained from the visible and near-infrared band displayed sharp jumps in the cloudy area and its edge (Figure 4a), indicating that the results of the geostationary satellite obviously affected the retrieval results. The latent heat flux retrieved from the microwave was more continuous and smoother in the same area, and there were no jumps in edge regions (Figure 4b).
Remote Sens. 2020, 12, x 10 of 15 Wm −2 (Figure 4b). However, the latent heat flux obtained from the visible and near-infrared band displayed sharp jumps in the cloudy area and its edge (Figure 4a), indicating that the results of the geostationary satellite obviously affected the retrieval results. The latent heat flux retrieved from the microwave was more continuous and smoother in the same area, and there were no jumps in edge regions (Figure 4b). There was no significant difference in the spatial distribution of the retrieved latent heat fluxes between the two data sources. The two methods had specific advantages and disadvantages. Although the surface temperature retrieved by the FY-3C microwave was not affected by clouds, its resolution was about 25 km. The surface temperature retrieval needed to be spatially interpolated before being used, and the accuracy was slightly lower than that of the thermal infrared remote sensing. The polar-orbiting satellites observed variables instantaneously, whereas the geostationary satellites obtained one-hour average values.

Verification
Due to the footprint effect in the turbulence observations, the value observed by the eddy covariance system can not only be regarded as a value at a certain point but is of spatial representativeness from tens of meters to several kilometers. The time scale of the remote sensing was one hour, so the hourly observation data from the eddy covariance system were collected during the experiment to verify the estimated results. Figure 5 shows the comparisons between hourly surface energy flux measurements and the estimates on August 18th, 2019, including net radiation flux, soil heat flux, sensible heat flux, and latent heat flux. The net radiation flux and sensible heat flux estimated by satellite remote sensing were compared with the ground measurements, and the root mean square error were 108.31 Wm −2 and 68 Wm −2 respectively; the root mean square error of the estimated soil heat flux and that obtained by the TCAV method was 14.28 Wm −2 and the root mean square error of the estimated latent heat flux and the ground measurement was 45.64 Wm −2 . The errors were all within a reasonable range.  There was no significant difference in the spatial distribution of the retrieved latent heat fluxes between the two data sources. The two methods had specific advantages and disadvantages. Although the surface temperature retrieved by the FY-3C microwave was not affected by clouds, its resolution was about 25 km. The surface temperature retrieval needed to be spatially interpolated before being used, and the accuracy was slightly lower than that of the thermal infrared remote sensing. The polar-orbiting satellites observed variables instantaneously, whereas the geostationary satellites obtained one-hour average values.

Verification
Due to the footprint effect in the turbulence observations, the value observed by the eddy covariance system can not only be regarded as a value at a certain point but is of spatial representativeness from tens of meters to several kilometers. The time scale of the remote sensing was one hour, so the hourly observation data from the eddy covariance system were collected during the experiment to verify the estimated results. Figure 5 shows the comparisons between hourly surface energy flux measurements and the estimates on 18 August 2019, including net radiation flux, soil heat flux, sensible heat flux, and latent heat flux. The net radiation flux and sensible heat flux estimated by satellite remote sensing were compared with the ground measurements, and the root mean square error were 108.31 Wm −2 and 68 Wm −2 respectively; the root mean square error of the estimated soil heat flux and that obtained by the TCAV method was 14.28 Wm −2 and the root mean square error of the estimated latent heat flux and the ground measurement was 45.64 Wm −2 . The errors were all within a reasonable range. mates on August 18th, 2019, including net radiation flux, soil heat flux, sensible heat flux, and latent heat flux. The net radiation flux and sensible heat flux estimated by satellite remote sensing were compared with the ground measurements, and the root mean square error were 108.31 Wm −2 and 68 Wm −2 respectively; the root mean square error of the estimated soil heat flux and that obtained by the TCAV method was 14.28 Wm −2 and the root mean square error of the estimated latent heat flux and the ground measurement was 45.64 Wm −2 . The errors were all within a reasonable range.

The Relationship between the Conditional Expectation and the Temperature
Surface temperature is the most direct feedback of the surface to meteorological conditions and is one of the key factors affecting the net radiation distribution on the surface. The surface temperature and energy distribution showed certain relationships at the spatial scale due to the climate and land cover of the study area. The spatially conditional expectation of the surface energy components and surface temperature retrieved by remote sensing was derived and analyzed after smoothing (Figure 6a). The trend line was smoothed. The results show that the distribution of net radiation generally displayed a similar trend to that of sensible heat. From a starting point of 275 K, with the increase in temperature, the conditional expectations of net radiation and latent heat flux increased, until stabilizing at 298-301 K, and decreasing after 302 K. Sensible heat flux showed an upward trend from 275 to 288 K, but decreased slowly after 288 K. In general, the expected values of sensible heat, net radiation and latent heat flux from 275 to 280 K were consistent, showing an upward trend. However, starting from 290 K, the expected value of sensible heat decreased slightly, while the expected latent heat increased. The expected value of soil heat flux increased slowly. Therefore, in mid-August, the daily surface sensible heat flux in the SRYR over the Qinghai-Tibet Plateau rose with the increase in temperature below 288 K. The sensible heat flux changed slightly within the temperature range of 280-285 K, and the latent heat flux increased significantly between 280 and 298 K.

The Relationship between the Conditional Expectation and the Temperature
Surface temperature is the most direct feedback of the surface to meteorological conditions and is one of the key factors affecting the net radiation distribution on the surface. The surface temperature and energy distribution showed certain relationships at the spatial scale due to the climate and land cover of the study area. The spatially conditional expectation of the surface energy components and surface temperature retrieved by remote sensing was derived and analyzed after smoothing (Figure 6a). The trend line was smoothed. The results show that the distribution of net radiation generally displayed a similar trend to that of sensible heat. From a starting point of 275 K, with the increase in temperature, the conditional expectations of net radiation and latent heat flux increased, until stabilizing at 298-301 K, and decreasing after 302 K. Sensible heat flux showed an upward trend from 275 to 288 K, but decreased slowly after 288 K. In general, the expected values of sensible heat, net radiation and latent heat flux from 275 to 280 K were consistent, showing an upward trend. However, starting from 290 K, the expected value of sensible heat decreased slightly, while the expected latent heat increased. The expected value of soil heat flux increased slowly. Therefore, in mid-August, the daily surface sensible heat flux in the SRYR over the Qinghai-Tibet Plateau rose with the increase in temperature below 288 K. The sensible heat flux changed slightly within the temperature range of 280-285 K, and the latent heat flux increased significantly between 280 and 298 K.
values of sensible heat, net radiation and latent heat flux from 275 to 280 K were consistent, showing an upward trend. However, starting from 290 K, the expected value of sensible heat decreased slightly, while the expected latent heat increased. The expected value of soil heat flux increased slowly. Therefore, in mid-August, the daily surface sensible heat flux in the SRYR over the Qinghai-Tibet Plateau rose with the increase in temperature below 288 K. The sensible heat flux changed slightly within the temperature range of 280-285 K, and the latent heat flux increased significantly between 280 and 298 K. The above analysis shows that under the abundant surface soil moisture conditions, with the increase in temperature and net radiation, the effective energy distribution of the surface changed significantly with the increase in temperature; in particular, when the surface temperature exceeded 290 K, the effective energy was mainly used for surface ET. This also indicates that, when the surface temperature exceeded 290 K, its retrieval accuracy had a significant impact on the ET retrieval, i.e., the higher the surface temperature, the greater the ET retrieval error.

The Relationship between Evaporation Fraction, Latent Heat Flux, and Daily ET
The evaporation fraction represented the energy exchange between the lower atmosphere and the near-surface at a regional scale, that is, the distribution of effective energy between surface sensible heat and latent heat flux. Under ideal conditions, when the soil was at its dry limit and the soil moisture was not available for evaporation, the latent heat flux was approximately zero and the sensible heat flux reached its maximum. In the wet limit environment, in which the soil moisture was sufficient, the latent heat flux reached its maximum. The First ISLSCP Field Experiment (FIFE) indicated that the evaporation fraction remained relatively stable under typical clear sky conditions, with the exception of mornings and evenings, so the ET throughout the day could be derived from the evaporation fraction at a certain moment and the surface effective energy of the whole day [26,27]. The observation from geostationary satellites made it possible to obtain the hourly scale evaporation fraction by remote sensing and provided a data source for verification of the hourly scale evaporation fraction using the eddy covariance observation system. Figure 6b presents the retrieved and observed latent heat flux and evaporation fraction at an hourly scale. Although the evaporation fractions of the two were similar, with values between 0.51 and 0.65, the fitting straight line showed different trends. The observed evaporation fraction showed a more obvious trend of decreasing with the increase in latent heat flux, whereas the retrieved evaporation fraction displayed a slight increase. The two lines intersected when the latent heat flux was 210.19 Wm −2 . However, due to the lack of accurate model-driven and surface verification data-sets, which leads to errors in the model forcing field and uncertainty of the model parameters, model calibration and error analysis based on land surface observations was complicated. More importantly, SEBS is a closure energy model, which retrieves surface ET based on the energy residual method, whereas the eddy covariance system used by the station for surface flux observation is not closed, however, and is still widely used for surface ET observation.
The evaporation fraction was evaluated using the observations from the eddy covariance system and the retrievals from the remote sensing data. At 12:00, the observed evaporation fraction was the closest to the estimated evaporation fraction. Hence, 12:00 was defined as the analysis time. The retrieved evaporation fraction at this time was 0.57, which was considered to be the effective evaporation fraction of the day. The daily ET calculated based on this value was 4.85mm, and the value of the total ET derived from the evaporation fraction at other times was within the range of 4.04-5.25 mm. Compared with the cumulative value, the daily ET derived from the effective evaporation fraction was closer to the observed value with the same input forcing and the model physics schemes. Therefore, the evaporation fraction at 12:00 was used to obtain the daily ET and this method was able to eliminate the retrieval error at a certain moment without destroying the physical mechanism of the retrieval model. The method has the potential to become an effective means of retrieving ET using geostationary satellite observations.

Daily ET
The spatial distribution of daily ET on 18 August 2019 in the SRYR was obtained by collecting the evaporation fraction and accumulating the hourly latent heat flux (Figure 7). The regional ET was between 2.0 and 4.0 mm. The low-value area appeared in the north of the study area with ET below 2.5 mm. The maximum daily ET was about 4.0 mm in the southeast of the SRYR, indicating that the surface soil water, vegetation, water available for ET, and meteorological conditions were all conducive to ET. In contrast, the northern part of the study area is located at the northern foot of Tanggula Mountain, which is not favorable for energy transfer from the surface to the atmosphere due to the high altitude and harsh weather conditions.

Conclusions
In this research, the spatial distribution of energy fluxes and daily ET at different times of the day in the SRYR were estimated using the SEBS algorithm and ground meteorological observation data, combined with the data of visible/near-infrared and microwave observations from the Fengyun series of meteorological satellites. Considering the influence of convective clouds during the afternoon above the plateau, which prevents visible and near-infrared observation in cloud-cover conditions, we used the MWRI brightness temperature to estimate the surface temperature. The hourly change of latent heat flux presented a sound unimodal diurnal variation. As the solar radiation increased after sunrise, the latent heat flux gradually increased from 180.0 Wm −2 at 11:00 to 200.0 Wm −2 at 12:00. The latent heat flux was about 280.0 Wm −2 at 13:00. At 14:00, the value reached its daily maximum of around 350.0 Wm −2 , and in some areas reached 400.0 Wm −2 . As the solar altitude, radiation, and air saturation deficit weakened in the afternoon, at 16:00, we adopted the surface temperature retrieved from the microwave data and found the latent heat flux in the cloud-free area was about 280.0 Wm −2 , and the cloud-covered area 150.0 Wm −2 . Converting the latent heat flux to the actual land surface ET, it was found that the daily ET in most of the study area was between 2.0 mm and 4.0 mm. The results were within a reasonable error range and lay a foundation for realizing the commercialization of land surface ET in the future.
The focus of this research was to investigate the hourly retrieval of ET using satellite remote sensing data. Due to the limitations of ground observation data and the overpass time of polar-orbiting satellites, the accuracy of the algorithm under cloud cover was not further analyzed. In addition, cloud-free and cloud-covered represent theoretical conditions, and cloud total amount (CTA) was obtained using the radiation transfer equation and the data from the visible and near-infrared channels. This could be used to determine

Conclusions
In this research, the spatial distribution of energy fluxes and daily ET at different times of the day in the SRYR were estimated using the SEBS algorithm and ground meteorological observation data, combined with the data of visible/near-infrared and microwave observations from the Fengyun series of meteorological satellites. Considering the influence of convective clouds during the afternoon above the plateau, which prevents visible and nearinfrared observation in cloud-cover conditions, we used the MWRI brightness temperature to estimate the surface temperature. The hourly change of latent heat flux presented a sound unimodal diurnal variation. As the solar radiation increased after sunrise, the latent heat flux gradually increased from 180.0 Wm −2 at 11:00 to 200.0 Wm −2 at 12:00. The latent heat flux was about 280.0 Wm −2 at 13:00. At 14:00, the value reached its daily maximum of around 350.0 Wm −2 , and in some areas reached 400.0 Wm −2 . As the solar altitude, radiation, and air saturation deficit weakened in the afternoon, at 16:00, we adopted the surface temperature retrieved from the microwave data and found the latent heat flux in the cloud-free area was about 280.0 Wm −2 , and the cloud-covered area 150.0 Wm −2 . Converting the latent heat flux to the actual land surface ET, it was found that the daily ET in most of the study area was between 2.0 mm and 4.0 mm. The results were within a reasonable error range and lay a foundation for realizing the commercialization of land surface ET in the future.
The focus of this research was to investigate the hourly retrieval of ET using satellite remote sensing data. Due to the limitations of ground observation data and the overpass time of polar-orbiting satellites, the accuracy of the algorithm under cloud cover was not further analyzed. In addition, cloud-free and cloud-covered represent theoretical conditions, and cloud total amount (CTA) was obtained using the radiation transfer equation and the data from the visible and near-infrared channels. This could be used to determine whether the study area at the regional scale was under cloud cover, but the accuracy of determining whether a single spot was affected by cloud cover was yet to be considered. In the next step, it will be necessary to take manual field observation experiments, that is, to record the cloud amount at the observation spot for each moment, and analyze the difference between the measured ET and the estimated results under cloud cover. The focus will be to improve the estimation accuracy of the model under cloud cover for a more reliable estimate of regional ET under all weather conditions.