Comparison of Three Theoretical Methods for Determining Dry and Wet Edges of the LST / FVC Space : Revisit of Method Physics

Land surface temperature and fractional vegetation coverage (LST/FVC) space is a classical model for estimating evapotranspiration, soil moisture, and drought monitoring based on remote sensing. One of the key issues in its utilization is to determine its boundaries, i.e., the dry and wet edges. In this study, we revisited and compared three methods that were presented by Moran et al. (1994), Long et al. (2012), and Sun (2016) for calculating the dry and wet edges theoretically. Results demonstrated that: (1) for the dry edge, the Sun method is equal to the Long method and they have greater vegetation temperature than that of the Moran method. (2) With respect to the wet edge, there are greater differences among the three methods. Generally, Long’s wet edge is a horizontal line equaling air temperature. Sun’s wet edge is an oblique line and is higher than that of the Long’s. Moran’s wet edge intersects them with a higher soil temperature and a lower vegetation temperature. (3) The Sun and Long methods are simpler in calculation and can circumvent some complex parameters as compared with the Moran method. Moreover, they outperformed the Moran method in a comparison of estimating latent heat flux (LE), where determination coefficients varied between 0.45 ~ 0.66 (Sun), 0.47 ~ 0.68 (Long), and 0.39 ~ 0.57 (Moran) among field stations.


Introduction
Heat energy and water vapor that are required for atmospheric motion come mainly from the land surface.Moreover, land surface heat and momentum fluxes also determine the strength and stability of turbulence and diffusion in the boundary layer as well as control the wind, temperature, and humidity changes.Consequently, land surface heat and water budget play a significant role in regional and global climate.Remote sensing is very promising in studying land surface heat and water budget since it provides a relatively cheap and rapid way to obtain up-to-date information of land surface evapotranspiration (ET), temperature, and soil moisture over a large geographical area [1].
Many methods for estimating ET and soil moisture based on remote sensing data were developed, such as SEBI (Surface Energy Balance Index) [2], SEBAL (Surface Energy Balance Algorithm for Land) [3], METRIC (Mapping EvapoTranspiration at high Resolution with Internalized Calibration) [4], a thermal inertia model for soil moisture [5], and so on.Reviews about these methods can be found [6][7][8].However, developing a method that can avoid complex parameterization of aerodynamic and surface resistances for water and heat transfer and can enable the estimation using a remote sensing image without additional information is still a challenge [9][10][11].Currently, the models based on Land Surface Temperature and Fractional Vegetation Coverage (LST/FVC) space show the prospect of a Remote Sens. 2017, 9, 528 2 of 15 wide application to meet the challenge and, therefore, are of continuous interest to the remote sensing community [11][12][13][14].One key advantage of such LST/FVC space-based methods is their simplicity in terms of their implementation, as well as the requirement of easily obtained parameters from remote sensing data.Moreover, those methods also include all the advantages of both the optical and thermal infrared methods [8] such as providing fine spatial and temporal resolution, use of relatively mature technology, long historical data, etc.In addition, the LST/FVC space is also a classical model for drought monitoring based on remote sensing.Temperature Vegetation Drought Index (TVDI) [15] and the Vegetation Temperature Condition Index (VTCI) [16] are two typical drought indexes based on the LST/FVC model.
However, one fatal drawback of such LST/FVC space-based methods is their determination of theoretical boundaries of LST/FVC space, i.e., the dry and wet edges.Generally, there are three ways to solve this issue: visual recognition, automatic fitting, and theoretical calculation.The visual recognition way is easy to implement but relies on expert's experience, thus, it is too subjective with large uncertainties.The automatic fitting way determines the boundaries directly from remote sensing images and reduces the subjectivity and uncertainty of the visual recognition method [10].However, this method relies on the study area generating the LST/FVC spaces.In order to achieve a satisfactory effect, it requires the study area to have areas of extreme dry and wet as well as areas with full vegetation and bare surfaces.In addition, the above two ways generally suffer from all the characteristics of an empirically derived methodology, e.g., lack of transferability to other regions, fine-tuning, and weakness to describe physical processes.Conversely, theoretical calculation derives the boundaries from an energy balance equation of land surface and, thus, allows the determination of dry and wet edges independently from the study area [17][18][19][20].Although some parameters that are not easy to obtain by remote sensing are still necessary, the theoretical calculation method facilitates understanding the physical meanings and influence factors of dry and wet edges.Moreover, it facilitates finding a way in the future for combining the advantages of the above methods such as the objectivity and independency of the theoretical way and the simplicity and parameter accessibility of the automatic fitting way.In summary, as an objective and independent method, the theoretical calculation method is worthy of further study.
Currently, there are three typical methods belonging to the theoretical calculation method.They were, respectively, proposed by Moran et al. in 1994 [21], Long et al. in 2012 [18], and Sun Hao in 2016 [19].In this study, the three methods are called Moran method, Long method, and Sun method for short.The three typical methods were presented with different deduction processes and expressions, which means they are hard to understand and they are limited in their application and further improvement.In addition, there is very limited research that compares them in determining the dry and wet edges or in estimating latent heat flux (LE).Motivated by the significance and potential of theoretical calculation, the aims of this study are to compare the three typical theoretical calculation methods, revisiting their basic physics, and clarify their differences and relations.

Dry and Wet Edges
Figure 1 shows the schematic diagram of the dry and wet edges in a LST/FVC space.Points A and D represent the soil and vegetation components with the maximum water stress, respectively.Points B and C represent the soil and vegetation components with saturated water supply, respectively.Lines within the LST/FVC space are isopiestic lines of soil moisture availability [11,18].Line BC is the wet edge and AC is the dry edge.Most research supposes the dry and wet edges as a straight line [11,13].Following this assumption, to determine the dry and wet edges we only require the determination of LST over points A (T s max ), B (T s min ), C (T c min ), and D (T c max ).shortwave radiation (W/m 2 ); a ε is atmosphere emissivity which can be determined using the method provided in [22]; σ is a constant with a value of 5.67 × 10 −8 W/m 2 /K 4 ; n is a fraction coefficient (unitless); ρ is the air density (kg/m 3 );  For bare soil and full-cover vegetation canopy components, Formula Error!Reference source not found.has different expressions.For bare soil, we have the following formula: where the subscript s represents the parameters of the bare soil component; For the full-cover vegetation canopy, there is the following formula: According to the physical meanings of the dry and wet edges, their theoretical expressions can be derived from the Energy Balance Equation: where R n is the net downward radiative flux (W/m 2 ); G is the downward ground heat flux into the subsurface (W/m 2 ); H is the upward surface sensible heat flux (W/m 2 ); and LE is the upward surface latent heat flux (W/m 2 ).Generally, R n , G, and H are determined by the following equation: where α and ε are land surface albedo and emissivity (unitless); S d is the down-welling shortwave radiation (W/m 2 ); ε a is atmosphere emissivity which can be determined using the method provided in [22]; σ is a constant with a value of 5.67 × 10 −8 W/m 2 /K 4 ; n is a fraction coefficient (unitless); ρ is the air density (kg/m 3 ); c p is the specific heat of air at constant pressure (J/kg/K); r a is the aerodynamic resistance (s/m); and T a is the air temperature (K).
For bare soil and full-cover vegetation canopy components, Formula (2) has different expressions.For bare soil, we have the following formula: where the subscript s represents the parameters of the bare soil component; T s is the bare soil surface temperature; r a s is aerodynamic resistance above the bare soil surface.
For the full-cover vegetation canopy, there is the following formula: where the subscript c represents the parameters of the vegetation canopy component; T c is the vegetation canopy surface temperature; r a c is the aerodynamic resistance above the vegetation canopy.
In terms of n s and n c , default values of 0.35 and 0 are generally adopted [18,19].

Moran Method
In the Moran method, the following formula was utilized to determine LE of a full-cover canopy: where γ is psychrometric constant (kPa/K); r c is canopy resistance (s/m) to vapor transport.e c * is saturated vapor pressure with respect to canopy radiometric temperature T c ; e a is vapor pressure of the air.Combining the Formulas (1), ( 2) and ( 5), there is the following equation for fully vegetated surface: where e a * is saturated vapor pressure of the air and ∆ represents the slope of the saturated vapor pressure-temperature relation, i.e., ∆ = (e c * − e a * )/(T c − T a ).
For a full-cover and well-watered vegetation canopy, r c = r c m , and for a full-cover vegetation with no available water, r c = r c x .Therefore, T c at the dry and wet edges (T c max and T c min ) can be determined by: where r c m and r c x can be obtained from measurements of stomatal resistance and the leaf area index (LAI), i.e., r c m = r sm /LAI and r c x = r sx /LAI; r sm and r sx are minimum and maximum stomatal resistance, respectively, and they were generally set as r sm = 25-100 s/m and r sx = 1000-1500 s/m [21].In this study, r sm and r sx are set as 25 s/m and 1500 s/m, respectively, representing the minimum and maximum values in the above value range.
For saturated bare soil, r c = 0 and for dry bare soil, r c = ∞.Therefore, bare soil radiometric temperature at the dry and wet edges (T s max and T s min ) can be determined by: Formulas ( 7) and ( 8) are the theoretical expressions of the original Moran method, where a single value of R n (measured on-site) was used [21].However, Formulas (3) and ( 4) indicate that the extreme component temperatures, i.e., T c max , T min c , T s max , and T min s are inputs of the calculation of R n .By introducing the Formulas (3) and (4) to calculate (R n − G) s and (R n − G) c , we improved the original Moran method and finally acquired the following expressions:

Long Method
The Long method assumes that LE equals to zero on the dry edge and H equals to zero on the wet edge.On the wet edge, the assumption of H equaling to zero establishes the equation that T c min = T s min = T a .On the dry edge, the assumption of LE equaling to zero enables the following equation established: Combining Equations ( 1), (2), and (11), we obtain the following equation: For bare soil radiometric temperature T s at the dry and wet edges (T s max and T min s ), their expressions can then be determined by: For canopy radiometric temperature T c at the dry and wet edges (T c max and T min c ), their expressions are:

Sun Method
In the Sun method, the following formula was used to determine LE: where φ is a complex effective Priestley-Taylor's parameter (unitless); the other variables have the same meanings as in the above equations.Combining the Formulas (1), ( 2) and ( 15), we obtained: Finally, radiometric temperatures of full-cover vegetation at the dry and wet edges (T max c and T c min ) can be determined by: Radiometric temperatures of bare soil at the dry and wet edges (T max s and T min s ) can be determined by: where φ min and φ max are set as 0 and 1.26 according to previous research [10,[23][24][25][26][27].∆ and γ can be determined by some empirical equations such as in [23].r s a and r c a can be calculated using the equations provided in [18].

Study Area and Materials
In 2012, there was an observation experiment in the Heihe River Basin, China.This observation experiment was called Heihe Watershed Allied Telemetry Experimental Research (HiWATER).It was a comprehensive experiment, which integrated the observations by satellite and airborne remote sensing as well as ground stations.Multiscale Observation Experiment on Evapotranspiration (MUSOEXE) is a subproject of HiWATER.In the MUSOEXE, several ground-based observation sites were established in Zhangye oasis located in the middle reaches of the Heihe River Basin.They were selected according to the crop structure, shelterbelt, residential area, soil moisture, and irrigation status.At each site, an automatic weather station and an eddy covariance (EC) system were installed to measure the surface fluxes of momentum, energy, and water vapor.
In this study, the observations by the EC systems and weather stations at five sites in the corn field were obtained from the "Heihe Plan Science Data Center, National Natural Science Foundation of China" (http://www.heihedata.org).These five sites are No. 10 to No. 14. Figure 2 shows the geographical location of the observation sites over oasis surfaces with a false-color Landsat-7 image of 10 July 2012 to show the underlying background.Several meteorological parameters were measured by the observation sites listed in Figure 2 from June to September in 2012.Six parameters among them are eligible for our study including near surface air temperature T a ( • C), air relative humidity RH (%), solar down welling shortwave radiation S d (W/m 2 ), friction velocity µ * (m/s), sensible heat flux H (W/m 2 ), and latent heat flux LE (W/m 2 ).The H and LE measured using the EC system have been corrected for closure using the Bowen ratio closure method.More information about the HiWATER and MUSOEXE can be found in [28,29].
Additionally, four Moderate Resolution Imaging Spectroradiometer (MODIS) products covering the same period from June to September in 2012 were collected in this study.They are Land Surface Temperature/Emissivity Daily L3 Global 1km products (MOD11A1 and MYD11A1) and Vegetation Indices 16-Day L3 Global 1km product (MOD13A2 and MYD13A2).The criterions "Pixel produced", "good data quality", "average emissivity error ≤ 0.02", and "average LST error ≤ 2 K" were used to control LST data.Additionally, the vegetation index (VI) usefulness parameter in the MOD13A2 and MYD13A2 products Quality Assessment Science Data Sets was used to ensure the data quality of normalized difference vegetation index (NDVI).The NDVI data with a VI usefulness from 0000 (highest quality) to 1100 (lowest quality) were retained in this study, while the others were excluded.Daily NDVI during a 16-day interval were assumed to be invariant as done in [12,30].
The LST and NDVI were extracted directly from the pixel in which the observation station was located, because the land cover type in the study area is approximately homogeneous.Then, we matched the MODIS data with the in-situ observations by comparing the view time.The matched values were then utilized to compare the three theoretical methods for estimating the dry and wet edges.In the comparison, ε s and ε c were set as 0.95 and 0.98; α c and α s were set as 0.18 and 0.24; and the vegetation height (h c ) was set as 1.0 m according to the field observations over the study area [31].To calculate r c m and r c x from r sm and r sx , maximum possible LAI is required [21].Since the vegetation under each site is maize, a maximum possible LAI of 5.0 was selected in the Moran method [32].FVC was calculated from NDVI using a non-linear relationship, i.e., FVC = [(NDVI − NDVI min )/(NDVI max − NDVI min )] 2 where NDVI min and NDVI max were set as 0.2 and 0.86 [19].
were established in Zhangye oasis located in the middle reaches of the Heihe River Basin.They were selected according to the crop structure, shelterbelt, residential area, soil moisture, and irrigation status.At each site, an automatic weather station and an eddy covariance (EC) system were installed to measure the surface fluxes of momentum, energy, and water vapor.
In this study, the observations by the EC systems and weather stations at five sites in the corn field were obtained from the "Heihe Plan Science Data Center, National Natural Science Foundation of China" (http://www.heihedata.org).These five sites are No. 10 to No. 14. Figure 2 shows the geographical location of the observation sites over oasis surfaces with a false-color Landsat-7 image of 10 July 2012 to show the underlying background.Several meteorological parameters were measured by the observation sites listed in Figure 2   Additionally, four Moderate Resolution Imaging Spectroradiometer (MODIS) products covering the same period from June to September in 2012 were collected in this study.They are Land Surface Temperature/Emissivity Daily L3 Global 1km products (MOD11A1 and MYD11A1) and Vegetation Indices 16-Day L3 Global 1km product (MOD13A2 and MYD13A2).The criterions "Pixel produced", "good data quality", "average emissivity error ≤ 0.02", and "average LST error ≤ 2 K" were used to control LST data.Additionally, the vegetation index (VI) usefulness parameter in the MOD13A2 and MYD13A2 products Quality Assessment Science Data Sets was used to ensure the data quality of normalized difference vegetation index (NDVI).The NDVI data with a VI usefulness from 0000 (highest quality) to 1100 (lowest quality) were retained in this study, while the others were excluded.Daily NDVI during a 16-day interval were assumed to be invariant as done in [12,30].
The LST and NDVI were extracted directly from the pixel in which the observation station was located, because the land cover type in the study area is approximately homogeneous.Then, we matched the MODIS data with the in-situ observations by comparing the view time.The matched values were then utilized to compare the three theoretical methods for estimating the dry and wet edges.In the comparison, s ε and c ε were set as 0.95 and 0.98; c α and s α were set as 0.18 and

Estimation of LE
Furthermore, we compared the three theoretical methods for estimating LE by determining the dry and wet edges in a Two-source Model for estimating Evaporative Fraction (TMEF) [19].The LE was estimated based on the same TMEF model and the same input parameters but with different boundaries of the LST/FVC space.For convenience, the maximum and minimum canopy resistance for dry vegetation in the Moran method were set as infinity and zero.
Figure 5 shows the scatter plots between estimated and observed LE.The scatter plots indicate that estimated LE by the Sun method and that by the Long method are comparable, and they have greater differences with the estimated LE by the Moran method.The estimations of LE by the Sun and Long methods show better agreement with the in-situ observations of LE.Two specific evaluation indicators, determination coefficient (R 2 ) and Root-Mean-Square-Error (RMSE), are also presented in Figure 5.The R 2 values by the Sun and Long method are greater than those by the Moran method, and the RMSE values by the Sun and Long method are less than those by the Moran method.R 2 and RMSE values by the Sun and Long methods are very close.Taking the station of No. 11, for example, the R 2 and RMSE are 0.66 and 69.85 for the Sun method, respectively, 0.68 and 68.50 for the Long method, respectively, and 0.57 and 78.86 for the Moran method, respectively.The same trend can be found at the sites of No. 10, No. 12, No. 13, and No. 14.The R 2 and RMSE values at No. 14 are 0.57 and 75.69 for the Sun method, respectively, 0.59 and 75.83 for the Long method, respectively, and 0.46 and 88.13 for the Moran method, respectively.In summary, the Long and Sun methods are comparable with each other in this study.They are more appropriate than the Moran method in determining the dry and wet edges for estimating LE.

Differences and Relations
The three theoretical calculation methods are all derived from the land surface energy balance equation.Although they were presented in different forms, we found they have similar mathematical expressions by revisiting their method physics.The prime difference among them is the determination of LE.The Moran method requires canopy resistance to determine LE.It infers the dry and wet edges through the adjustment of canopy resistance.For the bare soil components, canopy resistances on the dry and wet edges are equal to infinity and zero values.For the vegetation components, canopy resistances on the dry and wet edges are approximated to some empirical values.The Long method directly assumes LE equaling to zero on the dry edge and assumes LE equaling to available energy (i.e., H equaling to zero) on the wet edge.The Sun method introduces the Priestley-Taylor method to determine LE and calculates the dry and wet edges by adjusting the Priestley-Taylor parameter.On the dry and wet edges, the Priestley-Taylor parameter is respectively set as zero and an empirical maximum.
The differences in calculating LE led to different boundaries of the LST/FVC space.For determining the dry edge, the three theoretical methods are almost the same.That is because LE on the dry edge is considered to approximate zero for all the three methods, especially when the canopy resistance of vegetation and the Priestley-Taylor parameter on the dry edge are, respectively, set as infinity and zero.With respect to the wet edge, the three theoretical methods show greater difference.Generally, Long's wet edge is a horizontal line equaling to T a .Sun's wet edge is an oblique line and higher than Long's wet edge.Moran's wet edge intersects Long's and Sun's wet edges with a higher T min s and a lower T min c .

Suggestions
Among the three methods, the Sun method and the Long method are suggested in this study.Firstly, Table 2 shows some general and special parameters of the three theoretical methods.This table and the calculation equations of the three methods indicate that the Moran method is the most complex one and requires the most input parameters.Particularly, it requires some parameters that are hard to obtain, such as water vapor deficit and canopy resistance.The Sun and Long methods are much simpler than the Moran method.Moreover, they circumvent some complex and barely accessible parameters (i.e., e * a − e a , r x c , r m c ).Secondly, we compared the three theoretical methods for estimating LE using the same TMEF model and the same input parameters.Results demonstrated that the Sun and Long methods outperform the Moran method.Thirdly, on the Moran's wet edge T min c is generally less than T a (see Figure 3), which implies that the sensible heat flux of vegetation components with saturated water supply is generally negative and its latent heat flux is greater than the available energy.

Special Parameters General Parameters
Moran method e Further, a sensitivity analysis was performed in this study.We compared the relative change of a new estimated T min c or T min s to its initial value with the relative change of a new input variable to its initial value.The selected parameters are S d , µ * , α s , α c , and T a with initial values of 798.80 (W/m 2 ), 0.25 (m/s), 0.24, 0.18, and 285.82 (K), respectively.Results are presented in Figures 6 and 7.It can be found that a slight perturbation in the input parameters would bring greater error in the Moran method than in the Sun and Long methods except for T a .A 20% perturbation in S d causes more than a 1 K error in the estimation of T min c and almost a 4 K error in estimating T min s by the Moran method.In contrast, a 20% perturbation in S d causes less than a 1 K error in T min c and less than a 3 K error in T min s by the Sun method.The same phenomena can be found for the perturbations in µ * , α s , and α c .As for the perturbation in T a , approximate sensitivity was found between the Moran and Sun methods.This sensitivity analysis provides another argument for our suggestion.by the Sun method.The same phenomena can be found for the perturbations in µ * , s α , and c α .As for the perturbation in a T , approximate sensitivity was found between the Moran and Sun methods.This sensitivity analysis provides another argument for our suggestion.α s , and α c change from −20% to 20%.Additionally, determining the wet edge temperature as T a in the Long method implies that the EF on the wet edge is a constant value of 1, but EF is controlled not only by T a but also by wind speed [23].Yang et al. (2015) indicated that the underestimation of temperatures for the wet edge is a main reason for underestimating EF in the hybrid dual-source scheme and trapezoid framework-based evapotranspiration model (HTEM) model [13].The wet edge by the Sun method is higher than that by the Long method.Therefore, the Sun method is preferred to determine the dry and wet edges of the LST/FVC space.

Conclusions
In this study, we revisited and compared three theoretical methods for determining the dry and wet edges of LST/FVC space, i.e., Moran, Long, and Sun methods.Results are concluded here: (a) For determining the dry edge, the Sun method is equal to the Long method.They are a little bit different from the Moran method.However, if the canopy resistance of dry full-covered vegetation is set as infinity, all of the three method are identical.At that time, LE on the dry edge is equal to zero.
(b) For determining the wet edge, the three methods have greater differences.Long's wet edge is a horizontal line equaling air temperature.Sun's wet edge is an oblique line and is higher than that of Long's.Moran's wet edge intersects with the wet edges of Long and Sun methods with a higher soil component temperature and a lower vegetation component temperature.
(c) Among the three methods, the Sun and Long methods are suggested for estimating LE or monitoring variation of soil moisture, because they are simpler in calculation; they can circumvent some complex parameters; and they had better performance in estimating LE.We are more inclined to the Sun method since Long's wet edge is only correlated with air temperature.
From the Moran method to the Sun and Long methods, complex and barely accessible parameters decrease and performance in estimating LE is not necessarily lost.This implies that it is possible to determine the dry and wet edges using easily accessible parameters or parameters obtained from the image itself in the future.For example, air temperature is an essential parameter in the Sun and Long methods.There are already some methods for estimating the air temperature using satellite images [14, 30,33].Consequently, in future work, it will be worthwhile to improve further the Sun and Long methods, making them only rely on easily accessible parameters. d

pc
is the specific heat of air at constant pressure (J/kg/K); a r is the aerodynamic resistance (s/m); and a T is the air temperature (K).

Figure 1 .
Figure 1.Diagram of the dry and wet edges in Land surface temperature and fractional vegetation coverage (LST/FVC) space [19].

sT
is the bare soil surface temperature; s a r is aerodynamic resistance above the bare soil surface.

Figure 1 .
Figure 1.Diagram of the dry and wet edges in Land surface temperature and fractional vegetation coverage (LST/FVC) space [19].
from June to September in 2012.Six parameters among them are eligible for our study including near surface air temperature a T (°C), air relative humidity RH (%), solar down welling shortwave radiation d S (W/m 2 ), friction velocity μ * (m/s), sensible heat flux H (W/m 2 ), and latent heat flux LE (W/m 2 ).The H and LE measured using the EC system have been corrected for closure using the Bowen ratio closure method.More information about the HiWATER and MUSOEXE can be found in[28,29].

Figure 2 .
Figure 2. Geographic location of the stations.

Figure 2 .
Figure 2. Geographic location of the stations.

Figure 3
Figure 3 presents the estimation of T max s , T max

Figure 3 .
Figure 3. LST components on dry and wet edges estimated by the three theoretical methods on different sites.(a,b) are on No. 10 site; (c,d) are on No. 11 site; (e,f) are on No. 12 site.Figure 3. LST components on dry and wet edges estimated by the three theoretical methods on different sites.(a,b) are on No. 10 site; (c,d) are on No. 11 site; (e,f) are on No. 12 site.

Figure 4 .
Figure 4. Dry and wet edges estimated by the three theoretical methods on sites of (a) No. 10, (b) No. 11, (c) No. 12, and (d) No. 13 where the maximum and minimum canopy resistance were calculated using empirical parameters.

Figure 4 .
Figure 4. Dry and wet edges estimated by the three theoretical methods on sites of (a) No. 10, (b) No. 11, (c) No. 12, and (d) No. 13 where the maximum and minimum canopy resistance were calculated using empirical parameters.

Figure 4 .
Figure 4. Dry and wet edges estimated by the three theoretical methods on sites of (a) No. 10, (b) No. 11, (c) No. 12, and (d) No. 13 where the maximum and minimum canopy resistance were calculated using empirical parameters.
causes more than a 1 K error in the estimation of c T and almost a 4 K error in estimating s T by the Moran method.In contrast, a 20% perturbation in d S causes less than a 1 K error in min c T and less than a 3 K error in min s T

Figure 6 .Figure 6 .
Figure 6.Sensitivity analysis of the min c T by the three theoretical methods to input parameters

Figure 7 .
Figure 7. Sensitivity analysis of the m in s T by the three theoretical methods to input parameters

Figure 7 .
Figure 7. Sensitivity analysis of the T min s by the three theoretical methods to input parameters (a) S d , (b) µ * , (c) α s , and (d) T a

Table 1 .
Dry and wet edges estimated by the three theoretical methods on sites of No. 10~No.14 where the maximum and minimum canopy resistance were set as infinity and zero.FVC = fractional vegetation coverage.

Table 2 .
Input parameters of the three theoretical methods.