Next Article in Journal
Assessment of k-Nearest Neighbor and Random Forest Classifiers for Mapping Forest Fire Areas in Central Portugal Using Landsat-8, Sentinel-2, and Terra Imagery
Previous Article in Journal
A Tool for Pre-Operational Daily Mapping of Floods and Permanent Water Using Sentinel-1 Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Estimation of Evapotranspiration in Sparse Vegetation Areas by Applying an Optimized Two-Source Model

1
Institute of Forest Resource Information Techniques, Chinese Academy of Forestry, Beijing 100091, China
2
Key Laboratory of Forestry Remote Sensing and Information System, NFGA, Beijing 100091, China
3
National Academy of Forestry and Grassland Administration, Beijing 102600, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2021, 13(7), 1344; https://doi.org/10.3390/rs13071344
Submission received: 26 February 2021 / Revised: 23 March 2021 / Accepted: 27 March 2021 / Published: 1 April 2021

Abstract

:
Evapotranspiration (ET) is an important part of the water, carbon, and energy cycles in ecosystems, especially in the drylands. However, due to the particularity of sparse vegetation, the estimation accuracy of ET has been relatively low in the drylands. Therefore, based on the dry climate and sparse vegetation distribution characteristics of the drylands, this study optimized the core algorithms (canopy boundary resistance, aerodynamic resistance, and sparse vegetation coverage) and explored an ET estimation method in the Shuttleworth–Wallace two-layer model (SW model). Then, the Beijing–Tianjin sandstorm source region (BTSSR) was used as the study area to evaluate the applicability of the improved model in the drylands. Results show that: (1) The R2 value of the improved model results was increased by 1.4 and the RMSE was reduced by 1.9 mm, especially in extreme value regions of ET (maximum or minimum). (2) Regardless of the spatial distribution and seasonal changes of the ET (63–790 mm), the improved ET estimation model could accurately capture the differences. Furtherly, the different vegetation regions could stand for the different climate regions to a certain extent. The accuracy of the optimized model was higher in the semi-arid region (R2 = 0.92 and 0.93), while the improved model had the best improvement effect in the arid region, with R2 increasing by 0.12. (3) Precipitation was the decisive factor affecting vegetation transpiration and ET, with R2 value for both exceeding 0.9. The effect of vegetation coverage (VC) was less. This method is expected to provide a more accurate and adaptable model for the estimation of ET in the drylands.

Graphical Abstract

1. Introduction

Evapotranspiration (ET) plays a very important role in the ecosystem, connecting the global water, carbon, and energy cycles [1,2]. Water-limited ecosystems (including arid, semi-arid and sub-humid regions) account for more than 40% of the Earth’s land surface [3]. In these areas, 90% of the effective precipitation returns to the atmosphere through ET. This causes water resources to be the main limiting factor for economic development and ecological protection in arid areas [4,5]. Meanwhile, the drylands experience the frequent occurrence of climatic natural disasters, such as soil degradation and droughts. Early warnings of drought can be achieved through the dynamic monitoring of long-term ET [3,6,7,8,9]. Therefore, carrying out research on the accurate estimation of ET in drylands has great practical significance. Through this data, we can analyze the status of water profit, loss, and use efficiency. This can allow us to better guide agricultural production and ecosystem management [10,11].
The ET process is complex and observation sites are scarce; thus, it is difficult to precisely measure ET data obtained at large scales. However, remote sensing provides a cost-effective method to estimate regional or even global ET [12]. There are three main types of remote sensing methods used for estimating ET: The statistical empirical model method, the relationship between land surface temperature and vegetation index method, and the energy balance residual method [13,14,15,16]. The statistical empirical model method ignores the physical process of ET and relies heavily on training data. Therefore, it has a poor scalability to other areas and a low overall accuracy [17,18]. The relationship between land surface temperature and vegetation index method is based on establishing the “dry edge” and “wet edge” of regional ET. This method has certain physical significance and can improve the accuracy of ET estimation in a small area; however, this method still ignores the complex physical process of ET. It is difficult to obtain a theoretical “wet edge” in the drylands, which greatly affects the estimation accuracy of ET in these regions [19,20]. The energy balance residual method is used to simulate the physical process of ET, in which the physical meanings of the parameters are clear. It can estimate ET with a higher accuracy in different regions and is currently recognized as the best method for estimating ET [21,22].
There are three main categories in the energy balance residual method: single-layer models, two-layer models, and multi-layer models [23]. Single-layer models include the Penman–Monteith model [24,25,26], the improved Priestley–Taylor model [27,28], the Surface Energy Balance Algorithm for Land (SEBAL) [29,30], and the Surface Energy Balance System (SEBS) [31,32]. In a single-layer model, the vegetation is imagined as a “big leaf”, where all the material and energy exchanges between the atmosphere and the ground occur through this “big leaf”. Therefore, single-layer models can achieve good estimation results for areas with high vegetation coverage, but work poorly in the drylands with sparse vegetation [33]. Moreover, single-layer models can only calculate the ET and cannot obtain the vegetation transpiration and soil water evaporation separately. Although multi-layer models, such as the CO2 Energy Coupled Model (CECM) [34], have high accuracy in the estimation of ET, the physical processes are too complex to achieve ET estimation over a large area. This is because the model involves many variables and requires a large amount of known data as input. Two-layer models benefit from the advantages of both single- and multi-layer models. Taking into account the directional effect of radiation temperature, they clearly distinguish vegetation sources and soil sources in water and energy exchange, so the ET process was divided into two sources—vegetation and soil. Two-layer models can also estimate vegetation transpiration and soil water evaporation at the same time, which can improve the accuracy in estimating ET. Therefore, two-layer models provide one of the best methods for the estimation of ET using remote sensing data, especially for areas with sparse vegetation in the drylands [35].
The Shuttleworth–Wallace two-layer model (SW model) is a typical hybrid two-layer model based on the optimized Penman–Monteith (PM) model [36]. Due to its clear physical meaning and process, the SW model has relatively unique advantages for areas with sparse vegetation. It has been widely used in the remote sensing estimation of ET in the drylands [37,38,39]. Three two-layer models—the improved Penman–Monteith (PM) model, the Shuttleworth–Wallace (SW) two-layer model, and the two-patch (TP) model—have been compared for their adaptability to four ecosystems: deciduous forest, savanna, grassland, and farmland. The results showed that, compared with the other three ecosystems, all three models had relatively poor results in savanna areas [40].
In many previous studies on the estimation of ET in sparse vegetation areas [41,42,43], vegetation and soil have been used as a single layer in the two-layer models for parameter estimation and application research; however, the vegetation community within a single pixel often alternates between vegetation and soil. Remote sensing images with relatively low resolution cannot effectively capture such complex vegetation community characteristics within pixels, which is also the main reason for the poor simulation results of two-layer models in areas with sparse vegetation [44,45,46]. In view of such sparse vegetation characteristics, some studies have made adaptive improvements in some aspects [47,48]. The tree–grass ecosystems (TGEs) are a common ecological structure in the drylands. Vicente et al. [49] assumed that grassland and broadleaved trees were the dominant vegetation structure and cover separately in the growing period and dry period in semi-arid. On this basis, the leaf area index (LAI) input data of the original model were optimized to improve the model accuracy. However, the particularity of the distribution and ET process of the areas with sparse vegetation has not yet been fully grasped. Therefore, there are still large gaps in estimating ET of the drylands from biological process and sparse vegetation characteristics, including the canopy boundary resistance, aerodynamic resistance, vegetation coverage (VC), and so on.
The Beijing–Tianjin sandstorm source region (BTSSR) covers a wide area, and its climate and land-type changes are very complicated. Optimization of the ecological environment in this area can greatly improve the ecological environment in the Beijing–Tianjin–Hebei region, effectively curb the impact of meteorological disasters such as sandstorms, and slow down (or even improve) land degradation in central and northern China. The main objectives are: (1) Optimizing the core algorithm of the original SW model, in terms of the canopy boundary resistance, aerodynamic resistance, and regional ET estimation models, in order to adapt to the characteristics of dry climate and sparse vegetation in the drylands; and (2) Evaluating the applicability of the improved model in the drylands by analyzing the characteristics of spatial distribution and seasonal changes in vegetation transpiration, soil water evaporation, and ET in the BTSSR. This research is expected to improve the accuracy of the ET estimation in sparse vegetation areas and evaluate the water and energy cycle process of the ecosystem in the drylands.

2. Materials and Methods

2.1. Study Area

The BTSSR (36°52′–46°44′N, 105°10′–121°E) extends from Wulatehou Banner in Inner Mongolia in the west, to Aluhorqin Banner in Inner Mongolia in the east, to Dingbian County in Shaanxi in the south, and Dongwuzhumuqin Banner in Inner Mongolia in the north, covering 138 counties (banners, cities and districts) of six provinces (autonomous regions and municipalities): Beijing, Tianjin, Hebei, Shanxi, Inner Mongolia and Shaanxi. It has a total area of 710,500 square kilometers and a desertification area of 226,900 square kilometers, accounting for 31.9% of the total area of the region. A schematic diagram of the study area is shown in Figure 1.
The BTSSR landform consists of three regions: Plains, mountains and plateau. The area includes the Mu Us Sandy Land, the Kubuqi Desert, the Otindag Sandy Land and the Horqin Sandy Land. The western arid areas are mainly deserts and the Gobi, with bare land and low vegetation cover. The BTSSR includes a sub-humid arid zone, a semi-arid region, and an arid region, from southeast to northwest. The average annual precipitation decreases from about 450 mm in the southeast to about 150 mm in the northwest. The vegetation in the area is distributed sequentially, from the southeast to the northwest, as forest grassland, typical grassland, desert grassland, steppe desert, and desert, with large areas of sand vegetation distributed on the sandy land. According to the results of China’s vegetation zone division, the study area includes eight vegetation regions; namely, (1) A northern temperate meadow steppe subzone; (2) A temperate shrub, tree and semi-desert subzone; (3) A typical steppe subzone in the northern temperate zone; (4) A warm–temperate northern deciduous oak forest zone; (5) A temperate southern desert steppe subzone; (6) A typical steppe subzone in the southern temperate zone; (7) A temperate southern forest (meadow) grassland zone; and (8) A temperate shrub, semi-shrub and desert subzone. Their distribution is shown in Figure 2.

2.2. Data

The remote sensing data used were Moderate Resolution Imaging Spectroradiometer (MODIS, http://earthdata.nasa.gov, accessed on 25 February 2021) images, including MOD13A2, MOD15A2, MOD11A1, MCD43A1, MCD43B3 and MCD12Q1C. Table 1 shows the specific data parameter information. Preprocessing was carried out using the Google Earth Engine (GEE, https://code.earthengine.google.com, accessed on 25 February 2021) platform, including cropping and stitching actions. The remote sensing data were mainly used as the parameters of the SW two-layer model for calculation.
The meteorological data were obtained from the China Meteorological Data Network (http://data.cma.cn, accessed on 25 February 2021). They include the daily value data of meteorological stations from 365 days in 2019, including average temperature, relative humidity, precipitation, average wind speed, and soil moisture content. The station data were differentiated using the Kriging space difference method, in order to form surface data.
Digital elevation model (DEM) data (with a spatial resolution of 90 m) were obtained from the NASA Earth Observing System Data and Information System (EOSDIS, https://earthdata.nasa.gov/eosdis, accessed on 25 February 2021), for topographic correction when estimating the net surface radiation. According to the vectorized 1:1 million Chinese soil data provided by the Harmonized World Soil Database (HWSD) version 1.1, the saturated water content, field capacity, and wilting point of different soil types were calculated for the estimation of soil surface resistance and canopy resistance.
The correlation analysis data came from the Global Land Data Assimilation System (GLDAS). Bai et al. evaluated the accuracy and adaptability of GLDAS ET products in various parts of China, and their results showed that the ET error in China was 0.44 mm/year [50]. As the GLDAS ET products have a high accuracy in China, they were used as reference data to analyze the authenticity of the ET results of this study.

2.3. Methods

In this study, we optimized the original SW model, based on the characteristics of sparse vegetation and environmental differences in the drylands, in order to improve its ET estimation accuracy. The optimization model was mainly used due to the limitations of the original SW model, in terms of canopy boundary resistance, aerodynamic resistance, vegetation coverage (VC), and leaf area index (LAI). Figure 3 shows the research flow chart. Vegetation transpiration and soil water evaporation were calculated separately, where both required canopy boundary resistance and aerodynamic resistance as input parameters. The VC and LAI directly affect the net surface radiation, canopy boundary resistance, and aerodynamic resistance.

2.3.1. SW Model

In order to simulate the energy change in water flowing in different media (e.g., from soil water to plant water) and different state transitions (e.g., gaseous to liquid), the SW model introduces the concept of a resistance system. There are five types of resistance: aerodynamic resistance from canopy to reference height ( r a a ; s·m−1), aerodynamic resistance from the surface to the canopy height ( r a s ; s·m−1), canopy boundary resistance ( r a c ; s·m−1), soil surface resistance ( r s s ; s·m−1), and canopy stomatal resistance ( r s c ; s·m−1). These resistances are closely related to the physiological characteristics of the vegetation (e.g., stomatal conductance, photosynthesis, and so on), the atmospheric environment of the vegetation (e.g., solar radiation, atmospheric temperature, water vapor pressure, and so on), and the process of water and heat transfer between soil, vegetation, and the atmosphere.
The ET ( λ ET ) of the ecosystem can be decomposed into vegetation transpiration ( λ T ) and soil water evaporation ( λ E ) through the calculation of energy conversion by the resistance system:
λ ET = λ E + λ T = C s   P M s + C c P M c ,
P M c = Δ R + ρ C p D Δ r a c R s / r a a + r a c Δ + γ 1 + r s c r a a + r a c ,
P M s = Δ R + ρ C p D Δ r a s ( R R s ) / r a a + r a s Δ + γ 1 + r s s r a a + r a s ,
where P M c is the latent heat flux of the canopy (W·m−2), P M s is the latent heat flux of the soil covered by the canopy (W·m−2), C s and C c are the distribution coefficients, D is the saturated vapor pressure difference (KPa), Δ is the change rate of saturated water vapor pressure difference with temperature, γ is the psychrometer constant (0.067 kPa·K−1), R is the total effective energy (W·m−2), R s is the effective ground energy (W·m−2), ρ is the air density, and C p is the constant pressure specific heat.
The calculation formulae of the distribution coefficients C s and C c are as follows:
C c = 1 1 + A c A a A s A c + A a ,
C s = 1 1 + A s A a A c A s + A a
where:
A a = Δ + γ r a a ,
A c = Δ + γ r a c + γ r s c ,
A s = Δ + γ r a s + γ r s s .
The calculation formulae of the total effective energy R and ground effective energy R s are as follows:
R = R n G ,
R s = R n s G ,
R n s = R n exp C · LAI ,
where R n and R n s are the net radiant flux (W·m−2) of the canopy and the ground, respectively; G is the soil heat flux (W·m−2); and LAI is the leaf area index. According to the Lambert–Beer law, R n s is closely related to the leaf area index and vegetation type, and C is the attenuation coefficient of energy in the vegetation canopy community. In this study, different attenuation coefficients are substituted for different land-cover types.

2.3.2. Adaptability Optimization of SW Model in the Drylands

(1)
Optimization of Canopy Boundary Resistance r a c
The calculation of canopy boundary resistance in the original model simply divides the entire pixel into two parts—vegetation and bare soil—and the canopy boundary curve is the boundary between these two parts. This is an extreme situation. All vegetation gathers together and affects each other, so moisture is not easily lost to the outside air, resulting in increased canopy boundary resistance. In the drylands, the vegetation grows sparsely, and have little influence on each other. Therefore, the resistance of the canopy boundary is lower. If the calculation is still carried out according to the original model, there will be deviations (Figure 4). The left side of Figure 4 shows the accumulation of different vegetation within the pixel. The original models unified these, according to the ideal assumption on the right. Obviously, the two situations on the left differ, in terms of their canopy boundary resistance.
Even under the same VC and LAI, due to the difference in the degree of leaf aggregation within the pixel, the canopy boundary resistance will still vary greatly. Therefore, we creatively established the escaping index (EI), in order to revise the canopy boundary resistance in this situation. The biggest influencing factor of the EI is the canopy structure and vegetation community structure. Under the same VC and LAI, the denser the canopy structure and vegetation community structure are, the lower the EI is; meanwhile, the sparser the canopy structure and vegetation community structure are, the higher the EI is. Thus, the EI can be used to optimize the canopy boundary resistance of the original model in order to improve the calculation accuracy. Due to the differing degrees of vegetation accumulation, the shadow occlusion in different directions is different, which leads to different reflectivity in different directions. Therefore, the bidirectional reflectance distribution function (BRDF) contains vegetation aggregation and dispersion information. The large-area EI can be obtained based on remote sensing data and through the normalized difference hot spots and dark spots (NDHD) [51,52]:
EI = A θ s × NDHD + B θ s ,
where A and B are simulation coefficients, which vary with the solar zenith angle (SZA) (θs) and crown shape, and the definition of the NDHD is:
NDHD = ρ h ρ d ρ h + ρ d ,
where ρ h and ρ d represent the reflectivities of the hot and dark spots in the main plane, respectively. A hot spot is the peak of backscattering when the solar radiation is in the same direction as the viewing direction and there is no shadow. A dark spot is in the forward scattering direction away from a hot spot, where the largest shadow is observed, and corresponds to the smallest reflectivity in the main plane.
The ρ h and ρ d of each pixel are calculated from the BRDF product (MCD43) of the MODIS data. This product uses the RossThick-LiSparse-Reciprocal (RTLSR) model to calculate the anisotropic surface reflectance of each pixel [53,54]. The core is represented by three kernel functions: isotropic scattering, volume scattering and geometric optical surface scattering. The values of   ρ h and ρ d are related to the vegetation coverage. However, there are many low-coverage vegetation areas in the study area and, so, it is difficult to obtain hot spots without background interference. Therefore, for areas with low vegetation coverage, the relative observation angle needs to be increased. In this study, when the vegetation coverage is greater than 25%, the relative observation angle is the image observation angle; when the vegetation coverage is less than 25%, the relative observation angle is 60° [49].
The calculation formulae for the aerodynamic resistance inside the vegetation from leaf resistance ( r b ) to canopy resistance ( r a c ) are as follows [55,56,57]:
r a c = r b 2 L ,
r b = 100 / n w / u h 0.5 1 exp n 2 1 ,
where L is the LAI, w is the leaf width, and u h is the wind speed at the top of the canopy. In summary, the optimized formula for calculating the canopy boundary resistance is:
r a c = E I × r a c = E I 100 / n w / u h 0.5 1 exp n 2 1 2 L .
(2)
Optimization of Aerodynamic Resistance r a a and r a s
The aerodynamic resistance can be mainly divided into two parts: Closed vegetation and bare soil. The calculation formulae for the aerodynamic resistance of the closed vegetation area are shown in Equations (17) and (18):
r a s a = ln h r d z 0 k 2 · u h · e L 1 e L z 0 + d h L · h d ,
r a a a = ln h r d z 0 k 2 · u l n z d h d + h e L 1 z 0 + d h 1 L · h d ,
where   h r is the reference height, z 0 is the roughness, d is the zero-plane displacement, L is the LAI, h is the canopy height, u is the wind speed, and k is the Kaman constant.
The calculation formulae for the aerodynamic resistance of bare soil are shown in Equations (19) and (20):
r a s 0 = 1 k 2 · u ln h r z 0 l n z 0 + d z 0 ,
r a a 0 = 1 k 2 · u ln 2 h r z 0 r a s .
The aerodynamic resistance revision in the original SW model is based on the LAI [37,54]. When LAI > 4, it is assumed that the vegetation is completely closed, while when LAI < 4, it is assumed that there is a linear correlation between the aerodynamic resistance and the LAI. However, there are significant differences among different types of vegetation, regarding leaf type, crown structure, and tree shape. Therefore, there will be great errors in dividing closed vegetation or not when using a single LAI. In response to this problem, in this paper, we used the porosity to quantify the canopy closure of vegetation. Each pixel was divided completely into closed vegetation and bare soil, based on the porosity. Then, the aerodynamic resistance was calculated for the closed vegetation and bare soil, respectively. The formulae are as follows:
r a s = 1 P · r a s a + P · r a s 0 ,  
r a a = 1 P · r a a a + P · r a a 0 .
The calculation formula for porosity is based on the Lambert–Beer formula:
P θ = exp L G Ω / cos θ ,
where L is the LAI of the pixel, G is the leaf correlation coefficient, and Ω is the aggregation index.
(3)
Optimization of the ET Based on Poisson Distribution
The original SW model uses a single VC in each pixel to calculate the vegetation transpiration and soil water evaporation. However, in areas with sparse vegetation, the VC changes dramatically and different VCs may coexist within a 500 m pixel. For areas with sparse vegetation, vegetation growth can be regarded as random, and the value of each sub-pixel is a binary value of vegetation or bare soil. Therefore, from the perspective of the entire pixel, it can approximate a Poisson distribution. This can introduce the influence of different vegetation coverage into the model, thereby better characterizing areas with sparse vegetation by using a Poisson distribution to optimize the model. In this study, we divide the entire pixel into nine parts, each with a probability of 1–9 (VC × 10), where the VC × 10 of the entire pixel is used as the mean value of the Poisson distribution. The vegetation transpiration and soil water evaporation of a pixel are the sum of the vegetation transpiration and soil water evaporation of all sub-pixels. In this way, different vegetation coverage information can be fully utilized, and the characteristics of vegetation transpiration and soil water evaporation under different levels of vegetation coverage can also be expressed. The calculation formulae are as follows:
P X i = e VC m e a n · VC m e a n X i X i !   ,
E c = i = 1 n P   X i E c , i ,
E s = i = 1 n P   X i E s , i ,
where   X i   is the VC × 10 of sub-pixel i (with i ranging from 1 to 9), E c is the vegetation transpiration of the pixel, E c , i is the vegetation transpiration of sub-pixel i, E s is the soil water evaporation of the pixel, and E s , i is the soil water evaporation of sub-pixel i.
Similarly, the LAI can also be approximated as a Poisson distribution. The method is the same as that of VC; the formulae are as follows:
P X i = e LAI m e a n · LAI m e a n X i X i !   ,
E c = i = 1 n P   X i E c , i ,
E s = i = 1 n P   X i E s , i ,
where   X i   is the LAI of sub-pixel i (with i ranging from 1 to 9), E c is the vegetation transpiration of the pixel, E c , i is the vegetation transpiration of sub-pixel i, E s is the soil water evaporation of the pixel, and E s , i is the soil water evaporation of sub-pixel i.

2.3.3. Relative Validation of Results

As it is difficult to determine the relative true value of the ET in a large region, we adopted the cross-analysis method, which is a common method used in the analysis of ET estimation results [58,59,60]. This method uses a reliable surface ET product, which has been tested against a large number of absolute true values [61]. In addition, the cross-analysis method can also obtain the correlation results on a surface scale. Therefore, it can show the spatial distribution, which is important in order for us to grasp the improved model effect in different regions. The GLDAS ET products were used as the reference data in this study, and the root mean square error (RMSE) and R-squared (R2) were used to analyze the results. The calculation formulae are as follows:
RMSE = 1 n i = 1 n y i y ^ i 2 ,
R 2 = 1 i = 1 n y ^ i y i 2 i = 1 n y i y i ¯ 2 ,
where   y ^ i and y i are the reference data value and model-estimated value, respectively; y i ¯   is the average value of the reference data; and n is the total number of samples.

2.3.4. Changes Analysis of ET in the Vegetation Regions

The different vegetation regions could stand for the different climate regions to a certain extent. In order to better study the temporal and spatial characteristics of the ET and to explore the heterogeneity of ET in different vegetation regions, five statistics were calculated including minimum (min), maximum (max), mean (mean), median (med), and standard deviation (SD) of the eight vegetation regions and the entire study area. The minimum, maximum, and mean values were used to represent the numerical distribution of the vegetation regions, while the mean–median difference and standard deviation were used to characterize the heterogeneity of ET in this vegetation zone and to characterize the deviation when taking the mean value as the ET in the whole area.

3. Results and Analysis

3.1. Validation Results of the Model and the Reference Data

According to the analysis method introduced in Section 2.3.3, the effect of the optimization model is shown in Figure 5. Compared with the estimation results of the original model, different parts of the optimization model had different effects on the correlation coefficient. Among them, the R2 of the original model was 0.68 and the RMSE was 8.2 mm. After the three-part optimization algorithm was carried out at the same time, the R2 increased by 0.14 and the RMSE decreased by 1.9 mm. This indicates that the optimized model can effectively improve the estimation correlation of the original model. From the change in the scattered point distribution, it can be seen that the shape of the optimized model in Figure 5 was slenderer, which indicates that the optimized results were more stable and the results were more accurate. In particular, the optimization model was significantly optimized for the high-value area underestimated in the original model, while the low-value area was also significantly improved. Therefore, this indicates that the model has a significant improvement in extreme value regions (maximum or minimum) and has strong adaptability to the drylands.

3.2. Temporal and Spatial Distribution of ET in the BTSSR

According to the optimized SW model, the annual soil water evaporation, vegetation transpiration, and ET in the BTSSR in 2019 were obtained, as shown in Figure 6.
Because of the discontinuity and heterogeneity of vegetation in the drylands, the spatial distribution of the soil water evaporation, vegetation transpiration, and ET showed obvious differences. The improved model could accurately capture the differences. In the BTSSR, the highest value (261 mm) of soil water evaporation was observed in the northwest of the study area. The main land-cover type is Gobi in this area, with a low VC. The limited precipitation mainly dissipates into the atmosphere, in the form of soil water evaporation. On the other hand, the lowest value (19 mm) of soil water evaporation was in the southeast of the study area. The main land-cover type is forest in this area, with relatively closed VC. The precipitation is mainly dissipated in the form of vegetation transpiration. The regional distribution of vegetation transpiration showed an opposite trend to the regional distribution of soil water evaporation, being lower (>11 mm) in the northwest and higher (<573 mm) in the southeast. The highest values (790 mm) for ET were in the central and southern parts of the study area, while the lowest values (63 mm) were in the northwestern parts. It can be seen, from Figure 6, that no matter whether soil water evaporation, vegetation transpiration, or ET was considered, there were obvious spatial differences and gradual changes. Additionally, the higher the VC, the higher the vegetation transpiration, the higher the ET, and the lower the soil water evaporation.
On the time scale, we analyzed the spatial distribution of soil water evaporation, vegetation transpiration, and ET over the four seasons (January–March, April–June, July–September and October–December) in key areas of the study area. The selection factor for the key areas was to have as many vegetation regions as possible. The changes in precipitation and vegetation distribution were also complicated, and the types of land-cover were also as diverse as possible. The results are shown in Figure 7.
From the seasonal changes in Figure 7, we can see that the vegetation transpiration, soil water evaporation, and ET were significantly higher in the second and third quarters than in the first and fourth quarters, with a magnitude difference of about 10 times. In addition, vegetation transpiration accounted for a large proportion of the ET—exceeding 80% in most periods and regions—and the proportions in the second and third quarters were significantly higher than those in the first and fourth quarters. The seasonal spatial distribution of the soil water evaporation had little change, and there was no obvious change in the high- and low-value areas. Due to the influence of perennial vegetation, the high-value areas of vegetation transpiration in the first and fourth quarters disappeared in the second and third quarters, and were replaced by areas with better precipitation conditions and a higher vegetation coverage. As vegetation transpiration accounts for a relatively high proportion of the ET, the seasonal spatial variation in ET was similar to that of vegetation transpiration. This demonstrates that the precipitation conditions and annual maximum VC have a greater impact on the regional vegetation transpiration and ET. This is also the main reason why precipitation and VC were selected as influencing factors for the following analysis. In general, regardless of spatial or seasonal distribution, the improved model could accurately capture the differences in the drylands, which reflects the robustness of the improved model.

3.3. Temporal and Spatial Characteristics of ET in Different Vegetation Regions

The results of the five statistics are shown in Figure 8, in which the numbers 1–8 correspond to the vegetation regions in Figure 2 and the category “All” indicates the statistics referring to the entire study area.
Among the results, although the annual maximum occurred in the typical steppe subzone in the southern temperate zone, the mean and median of this area were lower than in the warm-temperate northern deciduous oak forest zone, while the mean-median differences and standard deviation were higher than those in this zone. Therefore, these results show that the overall ET was relatively high and the heterogeneity was low in the warm-temperate northern deciduous oak forest zone. The regional heterogeneity was relatively high and the distribution of ET was uneven in the typical steppe subzone in the southern temperate zone. On the other hand, the highest heterogeneity and the most uneven distribution of ET in this area was in the typical steppe subzone in the northern temperate zone. The smallest of the five statistics were observed in the temperate shrub, semi-shrub, and desert subzone, where the mean–median difference was also the lowest. This indicates that the annual ET in this area was not only scarce, but also very evenly distributed.
On the time scale, Figure 9 shows the daily average changes in vegetation transpiration and soil water evaporation in the eight vegetation regions in 2019, and the ET could be calculated by adding the two. The northern temperate meadow steppe subzone reached 5.31 mm on 21 July, which was the highest regional daily average value; however, the seasonal variation in ET in this area was large. This is mainly because this area is a typical grassland area, where there is more annual vegetation. This led to the ET from July to August being higher and being lower in other months. The opposite was true in the warm-temperate northern deciduous oak forest zone, where there are many arbor forests and the seasonal changes in ET are not very dramatic. By analyzing the daily average changes in the temperate shrub, semi-shrub, and desert subzone, it can be seen that, although the annual ET in this area was relatively low, there were still high values of daily averages (e.g., exceeding 4 mm) on certain dates, due to the influence of precipitation. The spatial distribution of the annual ET in the northern temperate meadow steppe subzone was relatively uniform, but with large seasonal variation, indicating that the spatial and temporal characteristics of the region were not consistent.
The different vegetation regions could stand for different climate regions to a certain extent. The comparative analysis of the eight vegetation regions on the daily scale is shown in Figure 10. The red dots and line mean the comparative result between the original model with reference data, while the green ones mean the comparative result between the improved model with reference data. Among them, the northern temperate meadow steppe subzone and the warm-temperate northern deciduous oak forest zone are in the semi-arid region, while the temperate shrub, semi-shrub, and desert subzone is in the arid region, and other vegetation regions belong to the transitional areas between the arid and semi-arid regions. By analysis, the accuracy of the optimized model was higher in semi-arid regions than in other regions, with R2 values exceeding 0.9. In the arid region, the model had the best improvement effect, with R2 increasing by 0.12. On the daily scale, the regions with the low ET (in the drylands) were the most difficult to estimate ET; However, the optimized model had the best improvement effect.

3.4. Influencing Factors on ET

From the input parameters of the SW model, it can be seen that there were three main types of factors affecting ET: The vegetation itself, soil topographic factors, and climate conditions. The relevant analysis diagram of influencing factors is shown in Figure 11.
After the analysis, it can be seen that the annual precipitation had a strong influence on the vegetation transpiration and ET, with R2 exceeding 0.9. This indicates that, in drylands, annual precipitation is the biggest influencing factor for vegetation transpiration and ET. The annual maximum NDVI (Normalized Difference Vegetation Index) also had a certain impact on ET, with an R2 of 0.52. This was because areas with a higher VC generally have higher rainfall and ET; however, these correlations are indirect effects. The annual maximum NDVI had a greater impact on vegetation transpiration than ET. This was because the VC affects the vegetation transpiration to a certain extent, and vegetation transpiration accounts for a far greater proportion of ET than soil water evaporation. However, in drylands, the vegetation is more drought-tolerant and the ET may be quite different under the same VC. Therefore, in general, the annual precipitation has a strong effect on the ET and the effect of VC is less in this area.

4. Conclusions

Taking into account the complex characteristics of sparse vegetation in the drylands, we optimized the canopy boundary resistance, aerodynamic resistance, VC, and LAI of the SW model. Then, we comprehensively analyzed the spatial distribution and temporal changes of the vegetation transpiration, soil water evaporation, and ET in the BTSSR in 2019. Results demonstrated that: (1) The R2 value of the improved model results was increased by 1.4 and the RMSE was reduced by 1.9 mm, especially in extreme value regions of ET (maximum or minimum). This indicates that the optimized model can effectively improve the estimation accuracy. (2) In general, the higher the VC, the higher the vegetation transpiration, the higher the ET, and the lower the soil water evaporation. Additionally, the seasonal spatial variation in soil water evaporation did not change much, while the vegetation transpiration and ET changed along with changes in precipitation and VC. Regardless of the spatial distribution and seasonal changes of the ET (63–790 mm), the improved ET estimation model could accurately capture the differences in the drylands. Furtherly, the different vegetation regions could stand for different climate regions to a certain extent. By analysis, the accuracy of the optimized model was higher in semi-arid regions (R2 = 0.92 and 0.93). In the arid region, the model had the best improvement effect, with R2 increasing by 0.12. (3) In terms of influencing factors, precipitation was the decisive factor affecting vegetation transpiration and ET, with R2 value for both exceeding 0.9. The effect of VC on the vegetation transpiration and ET in the region was less. The results show that the optimized SW model has high accuracy and applicability for the estimation of vegetation transpiration, soil water evaporation, and ET in drylands.

5. Discussion

Based on the climate and sparse vegetation characteristics of drylands, this study used the BTSSR as the study area and optimized the core algorithms, such as canopy boundary resistance, aerodynamic resistance, and sparse vegetation coverage, in the SW model. This method achieved the aim of further improving the estimation ability for regional ET in the drylands.
(1)
The estimation results of the improved model are based on the BTSSR, and it is very site-specific. It may have the applicability in other drylands of the world which have the similar characteristics of climate and vegetation, but this needs to be verified. Although the field data of regional ET is lacking to estimate the improved model, the spatial and temporal trends of ET in many previous studies were consistent with the results of this paper [62,63,64].
(2)
When calculating the canopy boundary resistance, the improved model assumes that the vegetation source and the soil source are on the same plane, and ignores the water and energy exchange in the vertical direction. This may lead to underestimating the results. In addition, the improved model assumes that there is indeed an area of bare soil within the pixel, and makes it more accurate in the sparse vegetation. As for the areas with continuous vegetation, it may be a little difficult to parameterize the model, and there is need to judge in advance whether the pixel is a continuous vegetation pixel.
(3)
Considering the other factors simplified in this paper, the SW model may further improve the estimation accuracy of ET. First, the precipitation interception was not calculated, which may have caused the ET to seem higher. Second, we ignored the impact of snowmelt. There are freezing periods, ranging from the north to south in the study area, based on the geographical location. Although there is less precipitation in winter, the area is still affected by snowmelt. Zhao et al. found that snowmelt can cause the SW model to estimate vegetation transpiration to be too high and soil water evaporation to be too low [65]. Third, when discussing the analysis of factors affecting regional ET, we only considered the two most important factors—precipitation and VC—in this study. The next step should further analyze the vegetation itself, as well as atmospheric and surface environmental factors, in order to make the results more accurate.

Author Contributions

Conceptualization, C.L., Z.L., Z.G. and B.S.; methodology, C.L.; software, C.L.; validation, C.L., Z.L., Z.G. and B.S.; formal analysis, C.L.; investigation, C.L., B.S.; resources, Z.L. and Z.G.; data curation, Z.L. and Z.G.; writing—original draft preparation, C.L.; writing—review and editing, Z.G., B.S.; visualization, C.L.; supervision, Z.L., Z.G. and B.S.; project administration, Z.L., Z.G. and B.S.; funding acquisition, Z.L., Z.G. and B.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Fundamental Research Funds for the Central Non-profit Research Institution of CAF (No. CAFYBB2019ZB004) and Commissioned project by State Forestry and Grassland Administration.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chen, Y. Comparison of satellite-based evapotranspiration models over terrestrial ecosystems in China. Remote Sens. Environ. 2014, 140, 279–293. [Google Scholar] [CrossRef]
  2. Jung, M. Recent decline in the global land evapotranspiration trend due to limited moisture supply. Nature 2010, 467, 951–954. [Google Scholar] [CrossRef] [PubMed]
  3. Liu, Q.; Yang, Z.; Cui, B. Spatial and temporal variability of annual precipitation during 1961–2006 in Yellow River Basin, China. J. Hydrol. 2008, 361, 330–338. [Google Scholar] [CrossRef]
  4. Sun, G. A general predictive model for estimating monthly ecosystem evapotranspiration. Ecohydrology 2011, 4, 245–255. [Google Scholar] [CrossRef]
  5. Wen, X. Support-vector-machine-based models for modeling daily reference evapotranspiration with limited climatic data in extreme arid regions. Water Resour. Manag. 2015, 29, 3195–3209. [Google Scholar] [CrossRef]
  6. Feng, S.; Fu, Q. Expansion of global drylands under a warming climate. Atmos. Chem. Phys. 2013, 13, 10081–10094. [Google Scholar] [CrossRef] [Green Version]
  7. Zhang, Q. Conversion features of evapotranspiration responding to climate warming in transitional climate regions in northern China. Clim. Dyn. 2019, 52, 3891–3903. [Google Scholar] [CrossRef]
  8. Shadmani, M.; Marofi, S.; Roknian, M. Trend analysis in reference evapotranspiration using Mann-Kendall and Spearman’s Rho tests in arid regions of Iran. Water Resour. Manag. 2012, 26, 211–224. [Google Scholar] [CrossRef] [Green Version]
  9. Espelta, J.M. “Mediterranean Oak Woodland Working Landscapes: Dehesas of Spain and Ranchlands of California” de Pablo Campos et al., 2013. Rev. Ecosist. 2014, 23, 79. [Google Scholar] [CrossRef] [Green Version]
  10. Du, T. Effects of distinguishing vegetation types on the estimates of remotely sensed evapotranspiration in arid regions. Remote Sens. 2019, 11, 2856. [Google Scholar] [CrossRef] [Green Version]
  11. Cleverly, J.R. Riparian ecohydrology: Regulation of water flux from the ground to the atmosphere in the Middle Rio Grande, New Mexico. Hydrol. Process. Int. J. 2006, 20, 3207–3225. [Google Scholar] [CrossRef] [Green Version]
  12. Zhang, K.; Kimball, J.S.; Running, S.W. A review of remote sensing based actual evapotranspiration estimation. Wiley Interdiscip. Rev. Water 2016, 3, 834–853. [Google Scholar] [CrossRef]
  13. Khosa, F.V. Evaluation of modeled actual evapotranspiration estimates from a land surface, empirical and satellite-based models using in situ observations from a South African semi-arid savanna ecosystem. Agric. For. Meteorol. 2019, 279, 107706. [Google Scholar] [CrossRef]
  14. Knipper, K.R. evapotranspiration estimates derived using thermal-based satellite remote sensing and data fusion for irrigation management in California vineyards. Irrig. Sci. 2019, 37, 431–449. [Google Scholar] [CrossRef]
  15. Li, Y. Mapping daily evapotranspiration based on spatiotemporal fusion of ASTER and MODIS images over irrigated agricultural areas in the Heihe River Basin, Northwest China. Agric. For. Meteorol. 2017, 244, 82–97. [Google Scholar] [CrossRef]
  16. Kalma, J.D.; McVicar, T.R.; McCabe, M.F. Estimating land surface evaporation: A review of methods using remotely sensed surface temperature data. Surv. Geophys. 2008, 29, 421–469. [Google Scholar] [CrossRef]
  17. Granata, F. evapotranspiration evaluation models based on machine learning algorithms—A comparative study. Agric. Water Manag. 2019, 217, 303–315. [Google Scholar] [CrossRef]
  18. Pan, S. Evaluation of global terrestrial evapotranspiration by state-of-the-art approaches in remote sensing, machine learning, and land surface modeling. Hydrol. Earth Syst. Sci. 2020, 24, 1485–1509. [Google Scholar] [CrossRef] [Green Version]
  19. Tang, R.; Li, Z.-L.; Tang, B. An application of the Ts–VI triangle method with enhanced edges determination for evapotranspiration estimation from MODIS data in arid and semi-arid regions: Implementation and validation. Remote Sens. Environ. 2010, 114, 540–551. [Google Scholar] [CrossRef]
  20. Vinukollu, R.K. Multi-model, multi-sensor estimates of global evapotranspiration: Climatology, uncertainties and trends. Hydrol. Process. 2011, 25, 3993–4010. [Google Scholar] [CrossRef]
  21. Chen, J.M.; Liu, J. Evolution of evapotranspiration models using thermal and shortwave remote sensing data. Remote Sens. Environ. 2020, 237, 111594. [Google Scholar] [CrossRef]
  22. Aboutalebi, M. Incorporation of unmanned aerial vehicle (UAV) point cloud products into remote sensing evapotranspiration models. Remote Sens. 2020, 12, 50. [Google Scholar] [CrossRef] [Green Version]
  23. Liou, Y.-A.; Kar, S.K. evapotranspiration estimation with remote sensing and various surface energy balance algorithms—A review. Energies 2014, 7, 2821–2849. [Google Scholar] [CrossRef] [Green Version]
  24. Hua, D. Uncertainty assessment of potential evapotranspiration in arid areas, as estimated by the Penman-Monteith method. J. Arid Land 2020, 12, 166–180. [Google Scholar] [CrossRef] [Green Version]
  25. McColl, K.A. Practical and Theoretical Benefits of an Alternative to the Penman-Monteith evapotranspiration Equation. Water Resour. Res. 2020, 56, e2020WR027106. [Google Scholar] [CrossRef]
  26. Yang, Y. Short-term forecasting of daily reference evapotranspiration using the reduced-set Penman-Monteith model and public weather forecasts. Agric. Water Manag. 2019, 211, 70–80. [Google Scholar] [CrossRef]
  27. Yao, Y. Differences in estimating terrestrial water flux from three satellite-based Priestley-Taylor algorithms. Int. J. Appl. Earth Obs. Geoinf. 2017, 56, 1–12. [Google Scholar] [CrossRef]
  28. Hao, Y.; Baik, J.; Choi, M. Developing a soil water index-based Priestley–Taylor algorithm for estimating evapotranspiration over East Asia and Australia. Agric. For. Meteorol. 2019, 279, 107760. [Google Scholar] [CrossRef]
  29. Glenn, E.P. Integrating remote sensing and ground methods to estimate evapotranspiration. Crit. Rev. Plant Sci. 2007, 26, 139–168. [Google Scholar] [CrossRef]
  30. Ochege, F.U. Mapping evapotranspiration variability over a complex oasis-desert ecosystem based on automated calibration of Landsat 7 ETM+ data in SEBAL. GISci. Remote Sens. 2019, 56, 1305–1332. [Google Scholar] [CrossRef]
  31. Zhao, J. Higher temporal evapotranspiration estimation with improved SEBS model from geostationary meteorological satellite data. Sci. Rep. 2019, 9, 1–15. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Wang, D. Improving meteorological input for surface energy balance system utilizing mesoscale weather research and forecasting model for estimating daily actual evapotranspiration. Water 2020, 12, 9. [Google Scholar] [CrossRef] [Green Version]
  33. Stannard, D.I. Comparison of Penman-Monteith, Shuttleworth-Wallace, and modified Priestley-Taylor evapotranspiration models for wildland vegetation in semiarid rangeland. Water Resour. Res. 1993, 29, 1379–1392. [Google Scholar] [CrossRef]
  34. Zhan, X.; Kustas, W. A coupled model of land surface CO2 and energy fluxes using remote sensing data. Agric. For. Meteorol. 2001, 107, 131–152. [Google Scholar] [CrossRef]
  35. Gao, G.; Zhang, X.; Yu, T. evapotranspiration of a Populus euphratica forest during the growing season in an extremely arid region of northwest China using the Shuttleworth–Wallace model. J. For. Res. 2016, 27, 879–887. [Google Scholar] [CrossRef]
  36. Shuttleworth, W.J.; Wallace, J. Evaporation from sparse crops-an energy combination theory. Q. J. R. Meteorol. Soc. 1985, 111, 839–855. [Google Scholar] [CrossRef]
  37. Huang, S. Modeling evapotranspiration for cucumber plants based on the Shuttleworth-Wallace model in a Venlo-type greenhouse. Agric. Water Manag. 2020, 228, 105861. [Google Scholar] [CrossRef]
  38. Chen, H.; Huang, J.J.; McBean, E. Partitioning of daily evapotranspiration using a modified shuttleworth-wallace model, random Forest and support vector regression, for a cabbage farmland. Agric. Water Manag. 2020, 228, 105923. [Google Scholar] [CrossRef]
  39. Jiang, Z.-Y. Revealing the spatio-temporal variability of evapotranspiration and its components based on an improved Shuttleworth-Wallace model in the Yellow River Basin. J. Environ. Manag. 2020, 262, 110310. [Google Scholar] [CrossRef] [PubMed]
  40. Lu, H. A hybrid dual-source model of estimating evapotranspiration over different ecosystems and implications for satellite-based approaches. Remote Sens. 2014, 6, 8359–8386. [Google Scholar] [CrossRef] [Green Version]
  41. Bao, Y.; Daun, L.; Liu, T. Simulation of evapotranspiration and its components for the mobile dune using an improved dual-source model in semi-arid regions. J. Hydrol. 2021, 592, 125796. [Google Scholar] [CrossRef]
  42. Delogu, E.; Boulet, G.; Olioso, A. Evaluation of the SPARSE dual-source model for predicting water stress and evapotranspiration from thermal infrared data over multiple crops and climates. Remote Sens. 2018, 10, 1806. [Google Scholar] [CrossRef] [Green Version]
  43. Wu, B.; Zhu, W.; Yan, N. Regional Actual Evapotranspiration Estimation with Land and Meteorological Variables Derived from Multi-Source Satellite Data. Remote Sens. 2020, 12, 332. [Google Scholar] [CrossRef] [Green Version]
  44. Carpintero, E. Remote-Sensing-Based Water Balance for Monitoring of evapotranspiration and Water Stress of a Mediterranean Oak–Grass Savanna. Water 2020, 12, 1418. [Google Scholar] [CrossRef]
  45. Li, X. A simple and objective method to partition evapotranspiration into transpiration and evaporation at eddy-covariance sites. Agric. For. Meteorol. 2019, 265, 171–182. [Google Scholar] [CrossRef]
  46. Turner, J.S. Termites as Mediators of the Water Economy of Arid Savanna Ecosystems, in Dryland Ecohydrology; Springer: Berlin/Heidelberg, Germany, 2019; pp. 401–414. [Google Scholar]
  47. Bellvert, J. Feasibility of using the two-source energy balance model (TSEB) with Sentinel-2 and Sentinel-3 images to analyze the spatio-temporal variability of vine water status in a vineyard. Remote Sens. 2020, 12, 2299. [Google Scholar] [CrossRef]
  48. Cha, M.; Li, M.; Wang, X. Estimation of Seasonal evapotranspiration for Crops in Arid Regions Using Multisource Remote Sensing Images. Remote Sens. 2020, 12, 2398. [Google Scholar] [CrossRef]
  49. Burchard-Levine, V. Seasonal adaptation of the thermal-based two-source energy balance model for estimating evapotranspiration in a semiarid tree-grass ecosystem. Remote Sens. 2020, 12, 904. [Google Scholar] [CrossRef] [Green Version]
  50. Bai, P.; Liu, X. Intercomparison and evaluation of three global high-resolution evapotranspiration products across China. J. Hydrol. 2018, 566, 743–755. [Google Scholar] [CrossRef]
  51. Wei, S. Global 500 m clumping index product derived from MODIS BRDF data (2001–2017). Remote Sens. Environ. 2019, 232, 111296. [Google Scholar] [CrossRef]
  52. Wei, S.; Fang, H. Estimation of canopy clumping index from MISR and MODIS sensors using the normalized difference hotspot and darkspot (NDHD) method: The influence of BRDF models and solar zenith angle. Remote Sens. Environ. 2016, 187, 476–491. [Google Scholar] [CrossRef]
  53. Lucht, W.; Roujean, J.L. Considerations in the parametric modeling of BRDF and albedo from multiangular satellite sensor observations. Remote Sens. Rev. 2000, 18, 343–379. [Google Scholar] [CrossRef]
  54. Schaaf, C.B. First operational BRDF, albedo nadir reflectance products from MODIS. Remote Sens. Environ. 2002, 83, 135–148. [Google Scholar] [CrossRef] [Green Version]
  55. Li, X. Modeling cherry orchard evapotranspiration based on an improved dual-source model. Agric. Water Manag. 2010, 98, 12–18. [Google Scholar] [CrossRef]
  56. Shuttleworth, W.J.; Gurney, R.J. The theoretical relationship between foliage temperature and canopy resistance in sparse crops. Q. J. R. Meteorol. Soc. 1990, 116, 497–519. [Google Scholar] [CrossRef]
  57. Hu, Z. Modeling and partitioning of regional evapotranspiration using a satellite-driven water-carbon coupling model. Remote Sens. 2017, 9, 54. [Google Scholar] [CrossRef] [Green Version]
  58. Xu, T. Evaluation of twelve evapotranspiration products from machine learning, remote sensing and land surface models over conterminous United States. J. Hydrol. 2019, 578, 124105. [Google Scholar] [CrossRef]
  59. Lu, X. Partitioning of evapotranspiration using a stable isotope technique in an arid and high temperature agricultural production system. Agric. Water Manag. 2017, 179, 103–109. [Google Scholar] [CrossRef] [Green Version]
  60. Tadesse, H.K. Evaluating evapotranspiration estimation methods in APEX model for dryland cropping systems in a semi-arid region. Agric. Water Manag. 2018, 206, 217–228. [Google Scholar] [CrossRef]
  61. Jia, Z. Validation of remotely sensed evapotranspiration over the Hai River Basin, China. J. Geophys. Res. Atmos. 2012, 117. [Google Scholar] [CrossRef]
  62. Shan, N. Spatiotemporal trends of reference evapotranspiration and its driving factors in the Beijing–Tianjin Sand Source Control. Project Region., China. Agric. For. Meteorol. 2015, 200, 322–333. [Google Scholar] [CrossRef]
  63. Han, J. Effects of different land use types on potential evapotranspiration in the Beijing-Tianjin-Hebei region, North China. J. Geogr. Sci. 2019, 29, 922–934. [Google Scholar] [CrossRef] [Green Version]
  64. Bian, Y. Spatial distribution of potential evapotranspiration trends in the Inner Mongolia Autonomous Region (1971–2016). Theor. Appl. Climatol. 2020, 140, 1161–1169. [Google Scholar] [CrossRef]
  65. Zhao, P. Comparison of dual crop coefficient method and Shuttleworth–Wallace model in evapotranspiration partitioning in a vineyard of northwest China. Agric. Water Manag. 2015, 160, 41–56. [Google Scholar] [CrossRef]
Figure 1. Location of the study area.
Figure 1. Location of the study area.
Remotesensing 13 01344 g001
Figure 2. Schematic diagram of vegetation regions in the study area.
Figure 2. Schematic diagram of vegetation regions in the study area.
Remotesensing 13 01344 g002
Figure 3. Research flow chart.
Figure 3. Research flow chart.
Remotesensing 13 01344 g003
Figure 4. Schematic diagram of canopy boundary resistance calculation. The whole is a pixel, and each cell is a sub-pixel. Here, the green is vegetation, and the gray is bare soil. (a): The actual distribution of vegetation and bare soil inside the pixel. (b): The simplified distribution of vegetation and bare soil inside the pixel used in previous studies.
Figure 4. Schematic diagram of canopy boundary resistance calculation. The whole is a pixel, and each cell is a sub-pixel. Here, the green is vegetation, and the gray is bare soil. (a): The actual distribution of vegetation and bare soil inside the pixel. (b): The simplified distribution of vegetation and bare soil inside the pixel used in previous studies.
Remotesensing 13 01344 g004
Figure 5. Comparative analysis of the observed results and modeled results.
Figure 5. Comparative analysis of the observed results and modeled results.
Remotesensing 13 01344 g005
Figure 6. Spatial distribution of soil water evaporation (a), vegetation transpiration (b), and evapotranspiration (ET, c) in the study area.
Figure 6. Spatial distribution of soil water evaporation (a), vegetation transpiration (b), and evapotranspiration (ET, c) in the study area.
Remotesensing 13 01344 g006
Figure 7. Spatial distribution of vegetation transpiration, soil water evaporation, and ET in the four seasons in key areas.
Figure 7. Spatial distribution of vegetation transpiration, soil water evaporation, and ET in the four seasons in key areas.
Remotesensing 13 01344 g007
Figure 8. Basic statistics of the vegetation regions and the entire study area.
Figure 8. Basic statistics of the vegetation regions and the entire study area.
Remotesensing 13 01344 g008
Figure 9. Daily scale changes in vegetation transpiration, soil water evaporation in vegetation regions.
Figure 9. Daily scale changes in vegetation transpiration, soil water evaporation in vegetation regions.
Remotesensing 13 01344 g009
Figure 10. The comparative analysis of the eight vegetation regions on the daily scale.
Figure 10. The comparative analysis of the eight vegetation regions on the daily scale.
Remotesensing 13 01344 g010
Figure 11. Analysis of influencing factors on ET and vegetation transpiration.
Figure 11. Analysis of influencing factors on ET and vegetation transpiration.
Remotesensing 13 01344 g011
Table 1. The Moderate Resolution Imaging Spectroradiometer (MODIS) data products used in the paper.
Table 1. The Moderate Resolution Imaging Spectroradiometer (MODIS) data products used in the paper.
Short NameData ProductsSpatial ResolutionTime Resolution
MOD13A2Vegetation Indices1 km16 Day
MOD15A2Leaf Area Index500 m8 Day
MOD11A1Land Surface Temperature and Emissivity1 kmDaily
MCD43A1Bidirectional Reflectance Distribution Function, BRDF500 mDaily
MCD43B3Albedo1 km16 Day
MCD12Q1CLand Cover Type500 mYearly
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, C.; Li, Z.; Gao, Z.; Sun, B. Estimation of Evapotranspiration in Sparse Vegetation Areas by Applying an Optimized Two-Source Model. Remote Sens. 2021, 13, 1344. https://doi.org/10.3390/rs13071344

AMA Style

Li C, Li Z, Gao Z, Sun B. Estimation of Evapotranspiration in Sparse Vegetation Areas by Applying an Optimized Two-Source Model. Remote Sensing. 2021; 13(7):1344. https://doi.org/10.3390/rs13071344

Chicago/Turabian Style

Li, Changlong, Zengyuan Li, Zhihai Gao, and Bin Sun. 2021. "Estimation of Evapotranspiration in Sparse Vegetation Areas by Applying an Optimized Two-Source Model" Remote Sensing 13, no. 7: 1344. https://doi.org/10.3390/rs13071344

APA Style

Li, C., Li, Z., Gao, Z., & Sun, B. (2021). Estimation of Evapotranspiration in Sparse Vegetation Areas by Applying an Optimized Two-Source Model. Remote Sensing, 13(7), 1344. https://doi.org/10.3390/rs13071344

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop