iBet uBet web content aggregator. Adding the entire web to your favor.
iBet uBet web content aggregator. Adding the entire web to your favor.



Link to original content: https://doi.org/10.3390/rs70100447
On Uncertainties of the Priestley-Taylor/LST-Fc Feature Space Method to Estimate Evapotranspiration: Case Study in an Arid/Semiarid Region in Northwest China
Next Article in Journal
Offshore Wind Resources Assessment from Multiple Satellite Data and WRF Modeling over South China Sea
Previous Article in Journal
Radiometric Non-Uniformity Characterization and Correction of Landsat 8 OLI Using Earth Imagery-Based Techniques
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On Uncertainties of the Priestley-Taylor/LST-Fc Feature Space Method to Estimate Evapotranspiration: Case Study in an Arid/Semiarid Region in Northwest China

1
State Key Laboratory of Remote Sensing Science, Institute of Remote Sensing and Digital Earth, Chinese Academy of Sciences, Beijing 100101, China
2
University of Chinese Academy of Sciences, Beijing 100049, China
3
Joint Center for Global Change Studies, Beijing 100875, China
*
Author to whom correspondence should be addressed.
Remote Sens. 2015, 7(1), 447-466; https://doi.org/10.3390/rs70100447
Submission received: 29 September 2014 / Accepted: 26 December 2014 / Published: 31 December 2014

Abstract

:
Accurate evapotranspiration (ET) estimation is very crucial for water resource management, particularly for the arid and semi-arid region. The remote sensing-based Priestley-Taylor method (RS-PT method) can estimate ET at regional scale, using the feature space of remotely sensed land surface temperature (LST) and vegetation index (VI). This study evaluates the RS-PT feature space method over an arid and semi-arid region in northwest China using satellite data from the moderate-resolution space-borne sensor Advanced Along-Track Scanning Radiometer (AATSR), the observations from the high-resolution airborne sensor Wide-angle Infrared Dual-mode line/area Array Scanner (WiDAS) and ground measurements of heat fluxes collected in summer 2008. The results show that the mean difference for latent heat flux (LE) estimates resulting from different domain sizes is 69.5 W/m2. When using high-resolution images from airborne measurements, the dry boundary is strongly affected by the pixels of impervious surfaces, which lead to a mean difference of 15.36 W/m2 for LE estimates. In addition, the physically based Surface Energy Balance Index (SEBI) model is used to analyze the accuracy of dry/wet boundaries in the RS-PT method. Compared with the SEBI-estimated relative evaporative fraction (Λr), the RS-PT method underestimated Λr by ~0.11. For the RS-PT method, the uncertainty in the determination of the dry/wet boundaries has a significant impact on the accuracy of the ET estimate, not only depending on the size of the area to build the feature space, but also on the land covers.

1. Introduction

Evapotranspiration (ET) is of great importance in modeling the water and energy interactions between land and atmosphere, and is a primary process in the energy and water exchange between hydrosphere, atmosphere and biosphere in terrestrial ecosystems. A precise estimation of regional ET is considered essential for many applications, such as global environmental change, regional ecosystem, water resource management and sustainable development of agriculture.
Traditional techniques provide point measurements, which usually are not representative of large regions due to the heterogeneity of the land surfaces and the dynamic nature of heat transfer. The development of remote sensing techniques provides a unique opportunity to estimate regional ET [1]. Many methods using remotely sensed data to acquire the spatially continuous ET have been proposed in the last few decades, including surface energy balance based models, empirical statistical methods, the traditional ET estimation approaches based on Penman-Monteith (PM) equation, methods based on land surface model and data assimilation techniques [2,3,4], all of which have been shown to provide acceptable results over certain range of land cover types and conditions [5,6].
Land surface ET is affected by numerous factors, including soil properties, weather condition, plant biophysics, and surface available energy, which have all been embedded in many physical based algorithms. However, difficulties in acquiring spatially representative ancillary measurements (e.g., near surface meteorological data, soil properties, and plant biophysical information) and in parameterization of various surface resistances, also have hindered the application of methods based on physical parameterization for ET estimation. As for the methods based on land surface model and data assimilation techniques, they need a lot of data and parameterizations as inputs and are relatively more computationally demanding than the remote sensing ET models. The empirical and semi-empirical statistical models, due to the lack of the physical foundation, need calibration to be applicable to local climate and ecosystems. Alternatively, methods built upon the triangle space between the land surface temperature (LST) and vegetation index (VI) e.g., Normalized Difference Vegetation Index (NDVI) and vegetation fractional cover (Fc) (referred to as the LST-Fc method) have been developed and widely used to estimate ET, evaporative fraction (EF) and soil moisture due to their simplicity and minimal request for inputs [7,8,9,10,11]. For instance, Jiang and Islam (1999) applied the triangle feature space method from the remotely sensed LST and NDVI to modify the coefficient of the Priestley-Taylor equation for actual ET estimates, which is called the RS-PT method. Such a method gives a clear physical meaning to the coefficient of the Priestley-Taylor equation and has been validated by several studies [12,13,14,15,16].
However, there are some uncertainties and limitations associated with the LST-Fc space method. Long et al. [17] assessed domain size dependencies of the LST-Fc feature space method using Landsat Thematic Mapper/Enhanced Thematic (TM/ETM) at the Soil Moisture-Atmosphere Coupling Experiment site in 2002. They indicated that the dry edge of the triangle method tends to move upward and the wet edge tends to move in the opposite direction when the domain size is increased, which could cause the a large uncertainties in EF estimates in terms of a mean absolute percentage difference (MAPD) up to 50%. Tian et al. [18] also quantitatively investigated variation of limiting edges due to the spatial domain size using ten Moderate Resolution Imaging Spectroradiometer (MODIS) observations in the Heihe River basin during growing season in 2009. Results showed that the extreme high surface temperature tends to increase and the extreme low surface temperature tends to decrease with the increase of domain size. Variation of intercept and slope of the dry edge would lead to the deviation in ET estimates with the maximum of 66 W/m2. Deviated from the previous studies, this study focuses primarily on offer a comprehensive view of the uncertainty in RS-PT method from a perspective of model physics and practical application through a study over arid/semiarid region in northwest China, including the preconditions of this methods and the influence of domain size as well as high-resolution scale effects. In addition, this study examines the accuracy of dry/wet boundaries from LST-Fc feature space by using a theoretical wet/dry edges derived from the physically based Surface Energy Balance Index (SEBI) model. Section 2 introduces fundamentals of RS-PT method and the theoretical dry/wet boundaries to assess dry/wet boundaries from LST-Fc feature space, followed by a description of study area and data collection in Section 3. Section 4 provides a systematic analysis of the uncertainties of RS-PT in application using moderate and high resolutions images, respectively, and a comparison between statistical and theoretical dry/wet boundaries. Conclusions and summaries are given in Section 5.

2. Methodology

2.1. Brief Overview of RS-PT Methods

2.1.1. The LST-Fc Feature Space

Although NDVI is commonly used as the VI in the RS-PT method, Fc is considered to be more appropriate than NDVI in the LST-VI diagram [15]. Therefore, we use Fc to construct the LST-VI triangle feature space. As shown in the LST-Fc diagram (Figure 1), point A represents the dry bare soil pixels with maximum LST and minimum Fc, point B represents well-watered soil at the potential ET with minimum LST. Point C represents well-watered vegetation pixels at potential ET and with minimum LST and maximum Fc. In Figure 1, the line AC is called the “dry edge”, along which Fc increases with decreasing LST under extreme dry surface conditions. The line BC is called the “wet edge”, along which the LST varies with the increase of Fc with potential ET (defined as the evaporation rate that would occur from a saturated surface with constant energy supply and constant atmospheric conditions). The wet and dry edges are derived statistically using linear regressions between the LST and Fc in the feature space diagram, for well-watered soil-vegetation canopy and completely dry soil conditions, respectively.
Figure 1. The concept of the LST-Fc triangular feature space.
Figure 1. The concept of the LST-Fc triangular feature space.
Remotesensing 07 00447 g001
In an ideal diagram, there is a point above the point C, which represents closed canopy pixels under water stress. The diagram would thus be a typical trapezoidal feature space [7]. However, in most practical applications using data from various remote sensors, this point is poorly defined or the position is very close to point C, so that a typical triangle shape in LST-Fc feather space is more often observed [8].

2.1.2. The RS-PT Method

The Priestley-Taylor equation, which has a clear physical meaning, has been widely used to calculate potential ET, provided advection can be neglected [19]:
L E = α [ ( R n G ) Δ Δ + γ ]
where LE (W/m2) is latent heat flux, α is the Priestley-Taylor (P-T) parameter, Rn and G (W/m2) are the net radiation flux and soil heat flux, respectively; Δ (Pa/K) is the slope of saturated vapor pressure with air temperature, and γ (Pa/K) is the psychometric constant. For well-watered soil-vegetation canopies, α is set as 1.26 [19]. Jiang and Islam proposed a modified version of the PT equation, which is to parameterize the P-T coefficient α based on the remotely sensed LST-VI feature space [12]. This method is referred to as the “remote sensing based Priestley-Taylor method” (RS-PT for short), and the symbol ϕ is used for the parameterized P-T coefficient, replacing the original coefficient α.
In the RS-PT model, the dry/wet boundaries in Figure 1 are straight lines and can be determined by a regression method applied to the scatter plot of LST and Fc. In this paper, the dry/wet boundaries are determined by an approach based on the scheme proposed by Tang et al., [15]. The value of ϕ for point A in Figure 1, corresponding to the driest soil with the lowest Fc, is set to 0. The value of ϕ for point B at the position of minimal Fc and minimal LST is set to 1.26. The value of ϕ along the wet edge is taken the same as that of point B, referred to ϕmax (= 1.26). The actual pixel-wise coefficient ϕ is derived using a two-steps interpolation scheme in the LST-Fc triangular space, as given below.
Firstly, the value of ϕ on the dry edge, referred as ϕmin,Fc, is linearly interpolated between 0 and 1.26 and changing with Fc,
ϕ min , F c = 1.26 F c
Secondly, the ϕ values under the same Fc are linearly interpolated between ϕmin,Fc and ϕmax, increasing with the decrease of LST [10,12],
ϕ = T max , F c T i T max , F c T min , F c ( ϕ max ϕ min , F c ) + ϕ min , F c
where Tmax,Fc and Tmin,Fc are LST values on the dry and wet edges, respectively, for any given FC. The Ti is the LST of pixel i (K).
The decrease of LST along the dry edge is due to the increase of the transpiration from vegetation with increasing Fc. For any given Fc, the ϕ increases linearly with LST, bounded by the ϕ values on the respective dry and wet edges (Equation (3)). In the LST-Fc triangle space, for a given FC, the relative size of LST is the indication of the soil moisture and the LST variation among different pixels is caused by the variation of soil moisture. Therefore, the LST-Fc triangle space implies that the LST is only affected by the variation of soil moisture and Fc, rather than by other factors, such as different atmospheric conditions or the topographical heterogeneity of the region [20,21,22].
From the above analysis, we can summarize the assumption and hypothesis in the RS-PT method as follows. (1) The dry/wet boundaries determined by the statistical method in the LST-Fc feature space are very close to the theoretical values; (2) For a given Fc, the evaporative fraction (EF) is a linear function of soil moisture availability; (3) The area of interest is under homogenous atmospheric forcing and with relatively small variations in terrain elevation.

2.2. Theoretical Dry/Wet Boundaries

In this study, SEBI algorithm proposed by Menenti and Choudhury [23] is used to obtain the theoretical-based dry/wet boundaries in the LST-Fc feature space. A brief description of SEBI is given below.
The relative evaporation fraction Λr, defined as the ratio of actual ET (λE) to potential ET (λEp), is expressed as:
Λ r = λ E / λ E P = 1 S E B I
where SEBI, according to Menenti and Choudhury (1993), is an index of standardized difference of temperature difference between the land surface and the air scaled by its dry and wet conditions and given below:
S E B I = ( T s T a ) r a ( T s , w T a ) r a,w ( T s , d T a ) r a,d ( T s , w T a ) r a,w
where Ts (K) is the land surface temperature, Ta (K) is the air temperature at reference height, Ts,w is the temperature at the wet-limit, Ts,d is the temperature at the dry-limit.
The (Ts,d – Ta) in Equation (5) is the surface-air temperature difference corresponding to the extreme dry condition:
T s , d T a = ( r a , d ) R n G ρ C p
where ρ (kg∙m−3) is the air density; Cp (= 1005 J∙kg−1∙K−1) is the specific heat of air at constant pressure.
The (Ts,w − Ta) in Equation (5) is the surface-air temperature difference for wet condition:
T s , w T a = ( r a , w ) R n G ρ C p 1 γ ( e s a t e ) 1 + Δ γ
where Δ (Pa/K) is the slope of e with air temperature and γ (Pa/K) is the psychometric constant, as before; e (Pa) is actual vapor pressure near surface; esat (Pa) is the saturated vapor pressure at mean air temperature. The ra (s∙m−1) in Equation (5) is the aerodynamic resistance and expressed as:
r a = 1 k u * [ ln ( z d 0 z 0 h ) ψ h ( z d 0 L ) + ψ h ( z 0 h L ) ]
where k is the von Karman constant (k = 0.4), z (m) is the reference height, zoh (m) is the roughness length for heat transfer, d0 (m) is the zero plane displacement height, u * (m∙s−1) is the friction velocity, Ψh is the stability correction function for turbulent heat flux.
The aerodynamic resistance for the theoretical boundaries, i.e., ra,d for zero evaporation (dry extreme) and ra,w for potential evaporation (wet extreme) are respectively defined as
r a , d = 1 k u * [ ln ( z d 0 z 0 h ) ψ h ( z d 0 L d ) + ψ h ( Z 0 h L d ) ]
r a , w = 1 k u * [ ln ( z d 0 z 0 h ) ψ h ( z d 0 L w ) + ψ h ( Z 0 h L w ) ]
The L (m) in Equation (8) is the Monin-Obukhov length and defined as,
L = ρ C p u * 3 θ v k g H
where g (m∙s−2) is the acceleration due to gravity, θv (K) is the potential virtual temperature near the surface.
In Equations (9) and (10), the Monin-Obukhov length for the two extreme limits are expressed, respectively, as,
L w = ρ a u * 3 k g ( H w c p T a + 0.608 R n G λ )
L d = ρ a u * 3 k g · ( R n G c p T a )
In Equations (5), (9) and (10), the subscript “w” denotes the variables (e.g., Ts,w and Lw) corresponding to the well-watered condition, where evapotranspiration take its potential value. The subscript “d” denotes the dry condition when evaporation becomes zero due to the limitation of water availability. Assuming that the air temperature at reference height does not change under two extreme conditions, Ts,d and Ts,w can be obtained using Equations (6) and (7). For each pixel (under a given Fc), the Ts,d and Ts,w can be treated as the theoretical dry and wet edges, respectively.
First, the theoretical dry/wet boundaries are estimated from SEBI, while the statistical boundaries are got from LST-Fc feature space method according to Section 2.1.1. Second, the Ts,w-FC and Ts,d-FC scatter plots are constructed just like the LST-FC triangular feature space (Figure 1), and same coordinate system is shared by these three scatter plots. Finally, the positions of the obtained two sets of boundaries are compared in the same frame, and the accuracy of the boundaries from RS-PT thus is evaluated. In addition, to quantitatively compare the discrepancy by dry/wet boundaries determining strategies, the Λr is also calculated for the RS-PT result according to Equation (4), in which the potential ET is computed by Penman-Monteith equation [24,25].

3. Study Area and Data

The study area is located in the Heihe River basin between 37.5°–39.6°N, 98.1°–100.6°E in northwest China, and exhibits a continental arid and semiarid climate. The study area encompasses the Qilian Mountains and the upper and middle reaches of the Heihe River basin. The land cover, which is provided by Heihe Plan Science Data Center [26], show that this area is dominated by sparse cover grassland, irrigated cropland, forest and desert as shown in Figure 2.
Figure 2. The location and the land cover of the study area. The small green rectangle in the map on the left is the Zhangye oasis.
Figure 2. The location and the land cover of the study area. The small green rectangle in the map on the left is the Zhangye oasis.
Remotesensing 07 00447 g002
Remote sensing data, including visible, near infrared and thermal infrared spectral domains, in two different spatial resolutions, were collected on 7 July 2008. One dataset is at 1 km spatial resolution from the Advanced Along-Track Scanning Radiometer (AATSR) on-board the Envisat satellite launched by the European Space Agency (ESA) and covers the entire study area indicated by the black rectangle in Figure 2. The second dataset is airborne remotely sensed data with 7.5 m spatial resolution, only covering the central part of the middle reaches of the Heihe River basin, as indicated by the green rectangle. The latter is acquired by the Wide-angle Infrared Dual-mode line/area Array Scanner (WiDAS) sensor during the Watershed Allied Telemetry Experiment (WATER) conducted in the Heihe River basin in 2008 [27,28].
To avoid the effect of measurements angular difference on the model behavior, only data with nadir view are adopted for LST, Fc and NDVI products. In addition, the MODIS albedo product, i.e., MCD43B3, which are available from the MODIS website [29], is also used in this study.
The meteorological forcing data over the AATSR domain, including wind speed, air temperature, air specific humidity, air pressure, downward shortwave and long wave radiations at 10 km spatial resolution and three-hourly temporal interval has been obtained from the China Meteorological Forcing Dataset developed by Data Assimilation and Modeling Center for Tibetan Multi-spheres, Institute of Tibetan Plateau Research, Chinese Academy of Sciences.
The flux measurements on 7 July 2008 have been collected from eddy covariance system (EC) at the Yingke oasis station [30]. More detailed information about the methods of data processing and quality control of EC flux measurements is given in [31]. The near surface meteorological variables over the WiDAS domain are from the automatic weather station (AWS) at Yingke oasis station.

4. Results and Discussion

The most important step in the RS-PT method is to determine the coefficient ϕ, for which the crucial point is to get the dry/wet boundaries in the LST-Fc diagram. Therefore, the discussion of this paper will focus on the identification of the dry/wet boundaries in the LST-Fc diagram and the subsequent estimation of coefficient ϕ.

4.1. Evaluation of the RS-PT Method with Moderate Resolution Data from AATSR

To evaluate the uncertainties in estimating the dry/wet boundaries and in turn their impact on the coefficient ϕ, the evaluation experiment is done by dividing the study area into three nested windows of different sizes, as shown in Figure 3:
-
Window A: the largest sub-region with elevations ranging from 1500 m to 5000 m and land cover consisting of forest land, cropland, grassland and unused land (e.g., Gobi and desert);
-
Window B: the intermediate sub-region, which encompasses the middle reaches of the Heihe River basin with elevations from 1500 m to 2180 m and consist of cropland and unused land;
-
Window C: the target area to be evaluated, which is flat with elevations ranging from 1500 m to 1730 m and land cover is dominated by cropland, located in the middle reach of Heihe River basin.
The sub-region C in Figure 3 covering mainly part of the middle reach of the Heihe River Basin was chosen as the target area to evaluate the RS-PT method. The pixel-wise ET is calculated for this target area using RS-PT method (Equations (1), (2) and (3)) with the coefficient ϕ calculated from LST-Fc feature space constructed from the three different image windows A, B and C. Figure 4 shows the spatial distribution of the coefficient ϕ and the corresponding LST-Fc feature space. The dry/wet boundaries, determined using statistical method, are also given in each diagram.
Figure 3. The spatial distribution of (A) Digital Elevation Model (DEM) (unit: m); (B) Fc and (C) LST (unit: K) in the entire study area. The three rectangles (A), (B), and (C) indicate the nested windows to build the LST-Fc feature space.
Figure 3. The spatial distribution of (A) Digital Elevation Model (DEM) (unit: m); (B) Fc and (C) LST (unit: K) in the entire study area. The three rectangles (A), (B), and (C) indicate the nested windows to build the LST-Fc feature space.
Remotesensing 07 00447 g003
The LST-Fc diagram of region A is nearly a standard triangle feature space (Figure 4A). Due to the elevation effect, the LST in mountain area is lower than that of cropland in the study area as shown in Figure 3C. As a consequence, many wet pixels with low LST along the wet edge might not necessarily be well-watered pixels but can be pixels at higher elevations. Although the shape of the LST-Fc feature space in Figure 4A seems to be perfect, the mis-defined well-watered pixels in the mountainous areas may cause errors in calculating the coefficient ϕ.
For the areas B and C, to remove the elevation effect, pixels with elevations greater than 1750 m were first masked out before constructing the LST-Fc diagrams. As shown in Figure 4B,C, the shapes of the LST-Fc diagrams from sub-regions B and C are obviously different. We used the scatter plot to evaluate the effect of study area size on the coefficient ϕ estimation. The P-T coefficient ϕ based on the image windows B and C as well as the resulting ET are compared and presented in Figure 5.
As shown in Figure 5, the coefficient ϕ derived from LST-Fc space of window C is larger than that from the LST-Fc space of window B, with mean difference of 0.13 and a standard deviation (STD) of 0.27 for the difference. The overestimation is more obvious for the pixels with higher ϕ values. For the pixels with ϕ < 0.2, i.e., near to the dry bare soil, the value ϕ from area B is lower than that derived from area C. The comparison of latent heat flux in Figure 5B gives similar result: the latent heat flux estimated using LST-Fc space of window C is larger than that from LST-Fc space of window B with a mean discrepancy of 69.5 W/m2 and STD of 114.7 W/m2. In theory, the value of ϕ for a given pixel depends on the relative position of the pixel between the dry and wet edges in the LST-Fc feature space. With different spatial coverage of the study area, the respective statistical boundaries also change in the varying LST-Fc feature space (e.g., Figure 4B,C). Thus, the relative position of any pixel changes with the varying wet or dry edge. This leads to a varying value ϕ with the change of size and location of the windows that determine the LST-Fc feature space as shown in Figure 4. From the above analysis, we can draw the conclusion that the difference both in the LST-Fc feature spaces and the coefficient ϕ between window B and window C is completely caused by the different size of the two windows and a resulting different range of land surface properties (e.g., Fc, LST and soil moisture).
Figure 4. The LST-Fc feature space constructed respectively using: (A) image window A; (B) image window B; and (C) image window C in Figure 3; (D,E,F)—the coefficient ϕ of the target area estimated using the three windows A, B and C in Figure 3 for LST-Fc feature space, respectively.
Figure 4. The LST-Fc feature space constructed respectively using: (A) image window A; (B) image window B; and (C) image window C in Figure 3; (D,E,F)—the coefficient ϕ of the target area estimated using the three windows A, B and C in Figure 3 for LST-Fc feature space, respectively.
Remotesensing 07 00447 g004
Figure 5. Comparison of estimated values for target area (window C) using LST-Fc feature space constructed from windows B and C respectively: (A) coefficient ϕ; and (B) latent heat flux.
Figure 5. Comparison of estimated values for target area (window C) using LST-Fc feature space constructed from windows B and C respectively: (A) coefficient ϕ; and (B) latent heat flux.
Remotesensing 07 00447 g005
It should be noted that the size of image window C is relative small compared with image window B, thus the range of soil moisture is narrow and the number of soil pixels with lower LST in window C is rather rare. Therefore, it is questionable that the wet boundary can be derived correctly from the LST-Fc diagram feature space constructed using the image window C.
The results show that the LST-Fc diagram does not only depend on the size of the area to build the LST-Fc feature space, but also on the ranges of the land surface Fc and soil moisture found in the area. The varying size of the sample window may lead to varying LST-Fc scatter plots and dry/wet boundaries, in turn, affect the coefficient ϕ. In addition, the area to construct the LST-Fc diagram should meet the precondition, which is to have a relatively small variation in elevation.

4.2. Evaluation of RS-PT Method with High Resolution Data from WiDAS

As the WiDSA is an airborne remote sensing sensor, the spatial coverage of WiDSA image is much smaller compared with that of satellite AATSR image. In this study, the image of WiDAS covers part of Zhangye Oasis, which is dominated by irrigated maize crops with village buildings and roads distributed in between. The WiDAS image was acquired on the same day during growing season, i.e., 7 July 2008.
Because the spatial resolution of WiDAS image is very high, the non-vegetated area (e.g., house roof, roads, pool, etc.) could be easily identified by visualization. As the pixels of water body are characterized by lower LST and strong evaporation, they can be treated as well-watered “soil” pixels. The pixels of impervious surface have a heating effect through absorbing solar energy rapidly. Due to the lower heat capacity with respect to that of bare soil, the temperature of an impervious surface is usually higher than that of a dry bare soil. At the coarse resolution of an AATSR image, the influence of impervious surface pixels on the LST-Fc diagram is not obvious due to the mixture of different elements in a pixel. At high resolution, however, these pixels of impervious surfaces appear explicitly around point A in Figure 1, corresponding to the maximum LST and minimal Fc. These pixels do not satisfy the definition of the dry pixel in Section 2.1. We refer to these pixels with impervious surfaces as “erroneous dry pixels” (fake dry pixels) and the derived dry edge based on the fake dry pixels is referred to as an erroneous dry edge (fake dry edge).
The results of the estimated ϕ and LST-Fc feature space measured by WiDAS are shown in Figure 6. The LST-Fc feature space in Figure 6A, referred to as original feature space, is the result constructed normally with all pixels in the image. The feature space in Figure 6B corresponds to the refined result obtained using refined feature space in Figure 6D, that is constructed with the residual pixels after removing the erroneous dry pixels.
Figure 6. The LST-Fc feature space constructed respectively by: (A) using the WiDAS image, denoted as original feature space; (B) using WiDAS image data without erroneous dry pixels, denoted as refined feature space. (C,D) the results of coefficient ϕ estimated from LST-Fc feature space (A) and (B) respectively.
Figure 6. The LST-Fc feature space constructed respectively by: (A) using the WiDAS image, denoted as original feature space; (B) using WiDAS image data without erroneous dry pixels, denoted as refined feature space. (C,D) the results of coefficient ϕ estimated from LST-Fc feature space (A) and (B) respectively.
Remotesensing 07 00447 g006
The maize was irrigated several times during the growing season. The irrigation schedule was not the same for different maize fields, this result in a relative broad range of the soil moisture for vegetation-covered pixels over this area. As a consequence, the diagram of WiDAS shows a typical trapezoid shape (Figure 6A,B). The comparison between the original feature space (Figure 6A) and the refined feature space (Figure 6B) shows that the position of the erroneous dry edge is higher than that of the real dry edge in the LST-Fc diagram. This would cause errors in the calculation of P-T coefficient ϕ and ET. Analogous to Figure 5, the scatterplots of original and refined results (Figure 7) were used to evaluate the influence of erroneous dry pixels on the estimation of the P-T coefficient ϕ and of ET. It is obvious that the values of coefficient ϕ derived from the original feature space are higher than that from the refined feature space due to its erroneous dry edge, and it is more obvious for the pixels with relatively lower ϕ. The difference in coefficient ϕ between the original feature space (Figure 6C) and the refined one (Figure 6D) ranges from −0.01 to 0.46 with a mean of 0.02 and STD of 0.05, and the corresponding difference in LE varies from −35 W/m2 to 167 W/m2 with a mean of 15.36 W/m2 and a STD of 23.65 W/m2. The results in Figure 7 demonstrate that the erroneous dry pixels cause an overestimation of both the coefficient ϕ and of the latent heat flux.
Figure 7. Comparison of original and refined result obtained from original and refined feature space: (A) for coefficient ϕ, and (B) for latent heat flux.
Figure 7. Comparison of original and refined result obtained from original and refined feature space: (A) for coefficient ϕ, and (B) for latent heat flux.
Remotesensing 07 00447 g007

4.3. Assessment of the RS-PT Boundaries by SEBI

The analysis in the previous section shows that the dry/wet boundaries derived from LST-Fc feature space have significant impact on the accuracy of LE estimates. In practice, we believe that the dry/wet boundaries derived from LST-Fc feature space are reasonable if the study area covers the full range of soil moisture and Fc conditions according to the assumption that the dry/wet boundaries determined by regression in the LST-Fc feature space are very close to the theoretical values (see Section 2). Therefore, we need to assess whether the dry/wet boundaries derived from the LST-VI feature space are reasonable and accurate.
In order to investigate the uncertainty associated with the statistical dry/wet boundaries in LST-Fc feature space, the Ts,d and Ts,w calculated by the SEBI model are used to investigate the position of the theoretical dry/wet boundaries. The results are presented in Figure 8. In the LST-Fc feature space, there is only one pair of values Ts,d and Ts,w for all the pixels at the same Fc value, that are on the dry/wet boundaries, respectively. However, according to the SEBI concept, the Ts,d and Ts,w are defined for each pixel, say pixel-dependent. Particularly, difference in aerodynamic resistance will also lead to different dry and wet edges for the same Fc value. That means, for each pixel there would be a corresponding Ts,d and Ts,w. That is the reason why the scatter plots of the Fc and Ts,d do not form a straight line but rather an area as shown in Figure 8, as do the scatter points of Fc and Ts,w.
Figure 8. Comparison between the position of theoretical dry/wet pixels from SEBI and the dry/wet boundaries derived from LST-Fc diagram using the moderate resolution AATSR image Window B. Red points represent the Ts,d and blue points represent the Ts,w. The cyan line indicates the dry edge and purple line indicates the wet edge derived by from LST-Fc feature space.
Figure 8. Comparison between the position of theoretical dry/wet pixels from SEBI and the dry/wet boundaries derived from LST-Fc diagram using the moderate resolution AATSR image Window B. Red points represent the Ts,d and blue points represent the Ts,w. The cyan line indicates the dry edge and purple line indicates the wet edge derived by from LST-Fc feature space.
Remotesensing 07 00447 g008
Comparing the cluster of the theoretical (SEBI) wet pixels and the position of RS-PT wet edge in Figure 8, the latter seems reasonable. However, there is a significant difference between the feature-space dry edge and the cluster of the theoretical (SEBI) dry pixels (Figure 8). The feature-space dry edge is much lower than the theoretical (SEBI) dry edge. Similar to the coefficient ϕ, the value of Λr depends on the relative position of the pixel value of surface temperature between the two extremes Ts,d and Ts,w in the LST-Fc scatter plot. As presented in Figure 8, the dry edge of the LST-Fc feature space is lower than the theoretical (SEBI) Ts, d, this is to say that Λr derived from the LST-Fc feature space would be lower than that deduced from SEBI. The spatial pattern and frequency curves of Λr derived from SEBI and RS-PT method over the target area are shown in Figure 9, demonstrating that, compared to SEBI, the RS-PT method underestimates surface wetness (i.e., smaller Λr) for most pixels in Zhangye Oasis. The mean discrepancy of Λr between SEBI and RS-PT is 0.11 with a STD of 0.32. The results in Figure 9 suggest that there are large uncertainties associated with the statistical boundaries from the LST-Fc feature space.
We also derived the theoretical dry/wet pixels for the WiDAS image Window B as shown in Figure 10. As described in Section 4.2, there is abundant irrigation to the maize fields in the growing season. For this reason, there should be many farmland pixels evaporating at potential ET on 7 July 2008. Therefore, the wet edge derived from the RS-PT method using LST-Fc feature space is very close to the theoretical Ts,w (Figure 10). However, the feature-space dry edge is still far from the theoretical (SEBI) dry edge. The comparison of Λr derived from SEBI and RS-PT is shown in Figure 11. Compared to Λr calculated from SEBI, the Λr from the RS-PT method is underestimated by 0.13 with a STD of 0.04. That is because the feature space dry edge is lower than the theoretical (SEBI) dry edge (Figure 10), as a consequence, the Λr derived from RS-PT method is lower than that from SEBI over the WiDAS image.
Figure 9. The spatial pattern and frequency curves of Λr derived from different methods over the Zhangye Oasis, using the same data as in Figure 8: (A) Λr is the result from the statistical dry/wet boundaries based on LST-Fc feature space (Λr-RS-PT); (B) Λr is calculated based on SEBI (Λr-SEBI); and (C) the frequency curves of Λr for (A) and (B).
Figure 9. The spatial pattern and frequency curves of Λr derived from different methods over the Zhangye Oasis, using the same data as in Figure 8: (A) Λr is the result from the statistical dry/wet boundaries based on LST-Fc feature space (Λr-RS-PT); (B) Λr is calculated based on SEBI (Λr-SEBI); and (C) the frequency curves of Λr for (A) and (B).
Remotesensing 07 00447 g009
Figure 10. Comparison between the position of theoretical dry/wet boundaries (pixels) from SEBI and the dry/wet boundaries derived from LST-Fc diagram using the high resolution WiDAS image. The labels are identical to those in Figure 8.
Figure 10. Comparison between the position of theoretical dry/wet boundaries (pixels) from SEBI and the dry/wet boundaries derived from LST-Fc diagram using the high resolution WiDAS image. The labels are identical to those in Figure 8.
Remotesensing 07 00447 g010
Figure 11. The spatial pattern and frequency curves of Λr from different methods using WiDAS data over Zhangye Oasis on 7 July 2008: (A) Λr is obtained using the LST-Fc feature-space (Λr_RS-PT); (B) Λr_SEBI is obtained based on SEBI (Λr-SEBI); and (C) the frequency curves of (A) and (B).
Figure 11. The spatial pattern and frequency curves of Λr from different methods using WiDAS data over Zhangye Oasis on 7 July 2008: (A) Λr is obtained using the LST-Fc feature-space (Λr_RS-PT); (B) Λr_SEBI is obtained based on SEBI (Λr-SEBI); and (C) the frequency curves of (A) and (B).
Remotesensing 07 00447 g011
Note that the Λr derived from RS-PT shows quasi Gaussian distribution but the SEBI-estimated Λr show asymmetric distribution peaked around 0.56 in the frequency curves (Figure 9C). The main reason for this is that the Λr associated with parameter ϕ derived from RS-PT is calculated using interpolation between 0 and 1; on the contrary the SEBI-estimated Λr completely depends on the environmental situation for each pixel. Considering the effect of advection on the oasis, which may cause the failure of energy balance closure, the ground-measured Λr without energy balance closure correction was adopted in current study. The ground-measured Λr from eddy covariance (EC) measurements at the Yingke station was obtained as reference value for the maize field in the whole study area [30]. Since it was the growing season of maize in July, the Λr value reached 0.97 at the sensor overpass time in current study. Therefore, it is reasonable that the value of the SEBI-estimated Λr = 0.96 for crop pixels.
From the analysis in Section 4.2, the erroneous dry points would lead to higher dry edge. However, compared with the theoretical dry pixels in Figure 10, the erroneous dry pixels are closer to the theoretical dry edge in this case. Therefore, the original feature space in Figure 6B is not reasonable, but may nevertheless be closer to the actual value than the refined feature space for this case. This apparent contradiction is still caused by the limitation of the method for the dry/wet boundaries depending on LST-Fc scatter plot in the RS-PT model.

4.4. Summary of the Issues about RS-PT Method

4.4.1. Issues about the Uncertainty in RS-PT Application

The varying size of the area of interest used to construct the LST-Fc feature space causes the different LST-Fc feature spaces and dry/wet boundaries, and therefore leading to uncertainties in the coefficient ϕ and ET estimates. The low resolution and large spatial coverage AATSR image captured a broader range of soil moisture and vegetation cover, which is essential to construct a reasonable feature space. However, the inhomogeneous atmospheric conditions and the difference of elevation over AATSR image may cause errors in the estimate on ET. Contrarily, the small area covered by the airborne WiDAS can be safely assumed to have homogenous atmospheric conditions and small elevation variation, but the range of soil wetness and Fc is likely to be too narrow. The latter may cause errors in the estimate of the dry and wet edges using the LST-Fc feature space method.

4.4.2. Issues about the Uncertainty on the Dry/Wet Boundaries

From the analysis in Section 4.3, we may conclude that the assumptions and preconditions must be met to derive reasonable wet and dry boundaries from LST-Fc feature space in the RS-PT method. The impact of the erroneous dry pixels must be must be paid enough attentions in application, particularly when the high-resolution images used include urban areas. Compared with the theoretical dry/wet boundaries from SEBI, it is relative difficult to get correct and accurate dry/wet boundaries just using the LST-Fc feature space in general, especially for the dry edge. Taking into account the effect of domain size effect on the wet and dry boundaries in the LST-Fc feature space in practical applications, one should pay much attention to the accuracy of the dry/wet boundaries when applying the RS-PT method for ET estimation.

5. Conclusions

This study evaluates the remote sensing-based Priestley-Taylor method (RS-PT method) over an arid and semi-arid region in northwest China from a comprehensive view, including the preconditions of this methods and the influence of domain size as well as high resolution scale effects, and examines the accuracy of dry/wet boundaries from land surface temperature-fractional vegetation cover (LST-Fc) feature space by using theoretical dry/wet boundaries. The uncertainties of the RS-PT method caused by domain size are analyzed using the moderate-resolution remotely sensed data from the Advanced Along-Track Scanning Radiometer (AATSR). The results show that when the research area is heterogeneous surface with the mountainous areas, the pixel along the wet edge might not necessarily be well-watered pixels but the pixels at higher elevations. The mean difference for Priestley-Taylor (P-T) coefficient (ϕ) estimates resulted from different domain sizes is 0.13 with a standard deviation (STD) of 0.27 and is 69.5 W/m2 with a STD of 114 W/m2 for latent heat flux (LE) estimates. When the high resolution images from airborne Wide-angle Infrared Dual-mode line/area Array Scanner (WiDAS) are used, it can be concluded that impervious surfaces leading to the erroneous dry pixels strongly affected the dry edge and resulted in the mean difference of 0.02 with a STD of 0.05 for ϕ estimates and of 15.36 W/m2 with a STD of 23.65 W/m2 for LE estimates. Therefore, the erroneous dry pixels from impervious surfaces must be paid enough attentions when the high-resolution images are used. The statistical dry/wet boundaries in LST-Fc feature space are also compared with the theoretical dry/wet boundaries derived from the physically based Surface Energy Balance Index (SEBI) model. In general, the statistical feature-space dry edge is lower than the SEBI-based theoretical dry edge, and the RS-PT method underestimates relative evaporative fraction (Λr) by ~0.11 compared with the SEBI-estimated values.
The applicability of the RS-PT method is restricted by the pre-conditions of homogeneous atmospheric conditions and a large range of land surface cover and soil moisture conditions. ET estimates from the RS-PT method strongly depend on the size of the area to build the LST- Fc feature space. However, how large the research area is the best scale for constructing LST- Fc feature space is still unclear. This work only analyzed the uncertainties of the LST-Fc method applied to the arid and semi-arid condition. More works need to be done to further quantify the uncertainties associated with dry/wet boundaries in the LST-Fc feature space method over different land cover and climate regions.

Acknowledgments

This work is jointly supported by the project of the National High Technology Research and Development Program of China (Grant No. 2012AA12A304), the CAS/SAFEA International Partnership Program for Creative Research Teams (Grant No. KZZD-EW-TZ-09).

Author Contributions

Zhansheng Li wrote the paper; Li Jia supervised the study, reviewed and edited the manuscript; Jing Lu reviewed and edited the manuscript. All authors read and approved the manuscript. We thank Ronald Hutjes for final language corrections.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. 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]
  2. Xu, T.; Liang, S.; Liu, S. Estimating turbulent fluxes through assimilation of geostationary operational environmental satellites data using ensemble kalman filter. J. Geophys. Res. Atmos. 2011, 116. [Google Scholar] [CrossRef]
  3. Caparrini, F.; Castelli, F.; Entekhabi, D. Estimation of surface turbulent fluxes through assimilation of radiometric surface temperature sequences. J. Hydrometeorol. 2004, 5, 145–159. [Google Scholar] [CrossRef]
  4. Sun, J.; Salvucci, G.D.; Entekhabi, D. Estimates of evapotranspiration from modis and AMSR-E land surface temperature and moisture over the southern Great Plains. Remote Sens. Environ. 2012, 127, 44–59. [Google Scholar] [CrossRef]
  5. Li, Z.L.; Tang, R.; Wan, Z.; Bi, Y.; Zhou, C.; Tang, B.; Yan, G.; Zhang, X. A review of current methodologies for regional evapotranspiration estimation from remotely sensed data. Sensors 2009, 9, 3801–3853. [Google Scholar] [CrossRef] [PubMed]
  6. Jia, Z.; Liu, S.; Xu, Z.; Chen, Y.; Zhu, M. Validation of remotely sensed evapotranspiration over the Hai River Basin, China. J. Geophys. Res. 2012, 117, D13113. [Google Scholar] [CrossRef]
  7. Moran, M.S.; Clarke, T.R.; Inoue, Y.; Vidal, A. Estimating crop water-deficit using the relation between surface-air temperature and spectral vegetation index. Remote Sens. Environ. 1994, 49, 246–263. [Google Scholar] [CrossRef]
  8. Price, J.C. Using spatial context in satellite data to infer regional scale evapotranspiration. IEEE Trans. Geosci. Remote Sens. 1990, 28, 940–948. [Google Scholar] [CrossRef]
  9. Jackson, R.D.; Idso, S.B.; Reginato, R.J.; Pinter, P.J. Canopy temperature as a crop water-stress indicator. Water Resour. Res. 1981, 17, 1133–1138. [Google Scholar] [CrossRef]
  10. Stisen, S.; Sandholt, I.; Nørgaard, A.; Fensholt, R.; Jensen, K.H. Combining the triangle method with thermal inertia to estimate regional evapotranspiration—Applied to msg-seviri data in the Senegal River Basin. Remote Sens. Environ. 2008, 112, 1242–1255. [Google Scholar] [CrossRef]
  11. Sun, Z.; Wang, Q.; Matsushita, B.; Fukushima, T.; Ouyang, Z.; Watanabe, M. A new method to define the VI-TS diagram using subpixel vegetation and soil information: A case study over a semiarid agricultural region in the north China Plain. Sensors 2008, 8, 6260–6279. [Google Scholar] [CrossRef] [Green Version]
  12. Jiang, L.; Islam, S. A methodology for estimation of surface evapotranspiration over large areas using remote sensing observations. Geophys. Res. Lett. 1999, 26, 2773–2776. [Google Scholar] [CrossRef]
  13. Jiang, L.; Islam, S. Estimation of surface evaporation map over southern great plains using remote sensing data. Water Resour. Res. 2001, 37, 329–340. [Google Scholar] [CrossRef]
  14. Jiang, L.; Islam, S. An intercomparison of regional latent heat flux estimation using remote sensing data. Int. J. Remote Sens. 2003, 24, 2221–2236. [Google Scholar] [CrossRef]
  15. 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]
  16. Yao, Y.; Qin, Q.; Ghulam, A.; Liu, S.; Zhao, S.; Xu, Z.; Dong, H. Simple method to determine the priestley–taylor parameter for evapotranspiration estimation using Albedo-VI triangular space from MODIS data. J. Appl. Remote Sens. 2011, 5, 053505:1–053505:16. [Google Scholar] [CrossRef]
  17. Long, D.; Singh, V.P.; Scanlon, B.R. Deriving theoretical boundaries to address scale dependencies of triangle models for evapotranspiration estimation. J. Geophys. Res. Atmos. 2012, 117. [Google Scholar] [CrossRef]
  18. Tian, J.; Su, H.; Sun, X.; Chen, S.; He, H.; Zhao, L. Impact of the spatial domain size on the performance of the TS-VI triangle method in terrestrial evapotranspiration estimation. Remote Sens. 2013, 5, 1998–2013. [Google Scholar] [CrossRef]
  19. Priestley, C.; Taylor, R. On the assessment of surface heat flux and evaporation using large-scale parameters. Mon Weather Rev. 1972, 100, 81–92. [Google Scholar] [CrossRef]
  20. Carlson, T.N.; Gillies, R.R.; Schmugge, T.J. An interpretation of methodologies for indirect measurement of soil-water content. Agr. For. Meteorol. 1995, 77, 191–205. [Google Scholar] [CrossRef]
  21. Nemani, R.; Pierce, L.; Running, S.; Goward, S. Developing satellite-derived estimates of surface moisture status. J. Appl. Meteorol. 1993, 32, 548–557. [Google Scholar] [CrossRef]
  22. Goward, S.N.; Cruickshanks, G.D.; Hope, A.S. Observed relation between thermal emission and reflected spectral radiance of a complex vegetated landscape. Remote Sens. Environ. 1985, 18, 137–146. [Google Scholar] [CrossRef]
  23. Menenti, M.; Choudhury, B.J. Parameterization of land-surface evaporation by means of location dependent potential evaporation and surface-temperature range. In Exchange Processes at the Land Surface for a Range of Space and Time Scales; Bolle, H.J., Feddes, R.A., Kalma, J.D., Eds.; IAHS Publication: Washington, DC, USA, 1993; pp. 561–568. [Google Scholar]
  24. Penman, H.L. Natural evaporation from open water, bare soil and grass. Proc. R. Soc. Lond. Ser. A. 1948, 193, 120–145. [Google Scholar] [CrossRef]
  25. Monteith, J. Evaporation and Environment. Available online: http://www.unc.edu/courses/2007fall/geog/801/001/www/ET/Monteith65.pdf (accessed on 29 December 2014).
  26. Multi-Source Integrated Chinese Land Cover Map. Available online: http://westdc.westgis.ac.cn (accessed on 5 March 2013).
  27. Water: Dataset of Airborne Widas Mission in the Zhangye-Yingke-Huazhaizi Flight Zone on 7 July 2008. Available online: http://westdc.westgis.ac.cn (accessed on 5 March 2013).
  28. Li, X.; Li, X.W.; Li, Z.Y.; Ma, M.G.; Wang, J.; Xiao, Q.; Liu, Q.; Che, T.; Chen, E.; Yan, G.J.; et al. Watershed allied telemetry experimental research. J. Geophys. Res. Atmos. (1984–2012) 2009, 114. [Google Scholar] [CrossRef]
  29. Reverb. Available online: http://reverb.echo.nasa.gov (accessed on 5 March 2013).
  30. WATER: Dataset of eddy covariance observations at the Yingke oasis station. Available online: http://westdc.westgis.ac.cn (accessed on 5 March 2013).
  31. Liu, S.; Xu, Z.; Wang, W.; Jia, Z.; Zhu, M.; Bai, J.; Wang, J. A comparison of eddy-covariance and large aperture scintillometer measurements with respect to the energy balance closure problem. Hydrol. Earth Syst. Sci. 2011, 15, 1291–1306. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Li, Z.; Jia, L.; Lu, J. On Uncertainties of the Priestley-Taylor/LST-Fc Feature Space Method to Estimate Evapotranspiration: Case Study in an Arid/Semiarid Region in Northwest China. Remote Sens. 2015, 7, 447-466. https://doi.org/10.3390/rs70100447

AMA Style

Li Z, Jia L, Lu J. On Uncertainties of the Priestley-Taylor/LST-Fc Feature Space Method to Estimate Evapotranspiration: Case Study in an Arid/Semiarid Region in Northwest China. Remote Sensing. 2015; 7(1):447-466. https://doi.org/10.3390/rs70100447

Chicago/Turabian Style

Li, Zhansheng, Li Jia, and Jing Lu. 2015. "On Uncertainties of the Priestley-Taylor/LST-Fc Feature Space Method to Estimate Evapotranspiration: Case Study in an Arid/Semiarid Region in Northwest China" Remote Sensing 7, no. 1: 447-466. https://doi.org/10.3390/rs70100447

Article Metrics

Back to TopTop