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

The VI-Ts diagram determined by the scatter points of the vegetation index (VI) and surface temperature (Ts) has been widely applied in land surface studies. In the VI-Ts diagram, dry point is defined as a pixel with maximum Ts and minimum VI, while wet point is defined as a pixel with minimum Ts and maximum VI. If both dry and wet points can be obtained simultaneously, a triangular VI-Ts diagram can be readily defined. However, traditional methods cannot define an ideal VI-Ts diagram if there are no full ranges of land surface moisture and VI, such as during rainy season or in a period with a narrow VI range. In this study, a new method was proposed to define the VI-Ts diagram based on the subpixel vegetation and soil information, which was independent of the full ranges of land surface moisture and VI. In this method, a simple approach was firstly proposed to decompose Ts of a given pixel into two components, the surface temperatures of soil (Tsoil) and vegetation (Tveg), by means of Ts and VI information of neighboring pixels. The minimum Tveg and maximum Tsoil were then used to determine the wet and dry points respectively within a given sampling window. This method was tested over a 30 km × 30 km semiarid agricultural area in the North China Plain through 2003 using Advanced Spaceborne Thermal Emission Reflection Radiometer (ASTER) and MODerate-resolution Imaging Spectroradiometer (MODIS) data. The wet and dry points obtained from our proposed method and from a traditional method were compared with those obtained from ground data within the sampling window with the 30 km × 30 km size. Results show that Tsoil and Tveg can be obtained with acceptable accuracies, and that our proposed method can define reasonable VI-Ts diagrams over a semiarid agricultural region throughout the whole year, even for both cases of rainy season and narrow range of VI.


Introduction
The VI-Ts diagram determined by the scatter points of remotely sensed vegetation index (VI) and land surface temperature (Ts) has been widely used to retrieve information on the partitioning of available surface energy [1][2][3] and surface moisture status [4][5][6][7][8][9]. All these applications are under the condition of homogenous atmospheric forcing such as solar radiation and air temperature over a sampling window for defining a VI-Ts diagram, so the size of the sampling window cannot be too large. If the sampling window covers full ranges of land surface moisture (from dry to well-watered) and VI (from bare soil to closed canopy), the VI-Ts diagram typically represents a right triangle when canopy temperature is assumed to be equal to Ts of well-watered bare soil [2,8,10] (Figure 1). The triangular VI-Ts diagram has been widely applied in previous studies [2][3][4][5]. The key point in these applications is how to define an ideal VI-Ts diagram, while the key point in the definition of the VI-Ts diagram is how to determine a dry edge in the VI-Ts diagram. Two automatic methods were proposed to define the dry edge in previous studies [5,9,11]. However, these two traditional methods require enough pixels that cover full ranges of land surface moisture and VI. Practically, it is difficult to find enough ideal pixels within a limited sampling window, especially when using a satellite data with moderate / coarse resolution, such as 1 km MODerate-resolution Imaging Spectroradiometer (MODIS) data. Generally, a natural land surface at the 1 km scale is usually a mixture of vegetation and nonvegetation (water and soil). Since water and soil surfaces have similar properties of low VI value, water surface can be considered as soil surface saturated by water. Therefore, if surface temperature (Ts) of a mixture pixel can be separated as the component surface temperatures of soil and vegetation (T soil and T veg , defined as subpixel information in this paper), two extreme surface conditions, dry bare soil (dry point) and closed vegetation (wet point), will probably be easily found at the subpixel scale. From Figure 1, the VI-Ts diagram will be readily defined if the dry and wet points are determined, which is independent of full ranges of land surface moisture and VI. Figure 1. The concept of a triangular VI-Ts diagram. Point A is called a dry point, and Point C is called a wet point in the VI-Ts diagram. AC is named the dry edge, and BC is named the wet edge.
There are some available methods to estimate component surface temperatures in previous studies. Assuming that a 'mixed' pixel was composed of only two temperature fields, a 'target' temperature and a 'background' temperature, the surface temperatures of both 'target' and 'background' were determined using two equations based on the land surface thermal radiative energy balance of channels 3 and 4 of NOAA-6/AVHRR [12]. This multi-channel method needs the condition of the same spectral response, which is difficult to meet actually. A Bayesian-based method was proposed to estimate cloud-top subpixel brightness temperature using multi-channel infrared GOES radiometer data [13]. This method used a detailed optical model of the GOES multi-channel imaging system, which is complicated and difficult in practical applications. A genetic inverse algorithm and multi-angle thermal infrared data were used to retrieve the component temperatures of mixed pixels [14]. However, there are almost no sensors onboard the existing satellites to provide multi-angle information except the Along Track Scanning Radiometer (ATSR) launched by European Space Agency. The morning/noon data was used to retrieve thermal inertia information, and then obtain component temperatures [15]. This thermal inertia method depended on ground experiments. Because of the limitations of the above methods, some researchers were inclined to estimate the component temperatures by means of the remotely sensed vegetation index. A linear two-source model was developed to estimate evaporation fraction where soil surface temperature within a pixel was obtained using a VI-Ts diagram [2], however vegetation surface temperature within a sampling window for defining the VI-Ts diagram was assumed to be homogeneous. In this study, a new practical approach was developed to obtain the component surface temperatures (vegetation / soil) based on the spatial autocorrelation of land surface moisture [16][17][18], in which the variation of vegetation surface temperature due to the difference of water deficit was considered.
The MODIS and the Advanced Spaceborne Thermal Emission Reflection Radiometer (ASTER) are onboard the NASA's Earth Observation System (EOS)-Terra satellite launched in 1999, both of which can provide high quality observations of land surface. MODIS was designed to collect observational data over a wide range at moderate resolutions (250, 500 and 1,000 m) with almost daily coverage of the Earth (http://modis.gsfc.nasa.gov/). ASTER captures high spatial resolution data in 14 bands, from the visible (15 m) to the thermal infrared (90 m) wavelengths, and provides a capability of stereo viewing (30 m) for the digital elevation model creation (http://asterweb.jpl.nasa.gov/). As the "zoom lens" for Terra, ASTER data can be used by other Terra and space-borne instruments for validation and calibration. Since both MODIS and ASTER are on the same satellite, ASTER provides an opportunity to validate MODIS observational data.
The purposes of this study were: (1) to propose a new approach to obtain T soil and T veg within a given pixel; (2) to propose a new practical method to define a VI-Ts diagram using the information of vegetation and bare soil components within pixels; (3) to validate the proposed method by using ASTER, MODIS and ground-based data; (4) to compare the proposed method with the traditional method across a semiarid agricultural region in the North China Plain through 2003.

Study Area and Ground Data Collection
The North China Plain (NCP) is one of main crop regions in China. The region displays a typical continental monsoon climate. The yearly mean air temperature is 13.1 °C, and the annual precipitation is about 610 mm, of which about 70% falls between June and August. Therefore, the NCP is zoned as a semiarid agricultural region. Our study area (30 km × 30 km) locates at the center of the NCP ( Figure 2). The light, temperature and water conditions support a 1-year 2-harvest cropping system (winter wheat (Oct.  Regular meteorological data recorded at the time when the EOS-Terra satellite overpassed our study area were collected from a flux station in the YES through 2003, including air temperature, humidity, wind speed, precipitation, downward and upward shortwave solar radiation, and downward and upward long-wave radiation [19]. The bulk temperature of an infinitely thick vegetation canopy is close to ambient air temperature [10], so observed air temperature can be used to validate Ts of wet points (Ts_wet, it corresponds to the minimum T veg ) in this study. Assuming that the latent heat flux is 0 over the surface of dry point, the Ts of dry point (Ts_dry, it corresponds to the maximum T soil ) can be obtained based on the energy balance over the dry bare soil surface [6,[20][21], a p a excess a n T C where R n and G are net radiation and soil heat flux in W/m 2 , ρ a and C p are the density and heat capacity of air in kg/m 3 and J/kg/°C, r a and r excess are aerodynamic resistance and excess resistance in s/m, and T a is air temperature in °C. The methods of r a and r excess calculations are the same to those in previous literatures [6,20,21]. In this study, Ts_dry estimated based on ground data was used to validate those derived from remote sensing data.

Satellite Data Collection and Processing
Two level-2 data products of ASTER (projection: UTM-50N; datum: WGS-84) over our study area on May 9 2003 were collected from the Japanese Ground Data System (http://www.gds.aster. ersdac.or.jp/gds_www2002/index_e.html): AST07 (the atmospheric corrected surface reflectance, resolution: 15 m) and AST08 (the surface temperature qualitatively assessed by cloud mask, resolution: 90 m). The AST08 is produced using the Temperature Emissivity Separation (TES) algorithm that yields accuracies around 0.01 for surface emissivity and 1 °C for radiometric temperature, respectively [22]. At the same observational time of ASTER (about 11:00 am of local standard time on May 9 2003), two MODIS data products in version 5 (projection: Sinusoidal; datum: WGS-84; resolution: 1 km) were also collected from the EOS data gateway (http://edcimswww.cr.usgs.gov/pub/ imswelcome/): MOD11 (the daily land surface temperature and emissivity masked by clouds) and MOD09 (the atmospheric corrected surface reflectance). The MOD11 product has been validated in South America, and results showed that the accuracy was better than 1 °C in the range from -10 °C to 50 °C [23]. Both MODIS and ASTER products are georegistered in their making processes. The accuracy of MODIS geolocation approximates 50 m at the nadir [24]. The hand-to-hand registration accuracy of ASTER is better than 0.2 pixels [25]. The reflectances of red and near-infrared bands were used to calculate the Normalized Difference Vegetation Index (MODIS-NDVI and ASTER-NDVI). ASTER-Ts was resampled to the resolution of 1 km from 90 m by averaging pixels in order to compare the system difference between MODIS and ASTER sensors. The projection of MODIS data was transferred from Sinusoidal to UTM-50N. The statistics of above datasets were listed in Table 1. It is found that both maximum and range of MODIS-NDVI are greater than those of ASTER-NDVI, which is consistent with the report in [26]. Due to the pixel-average scaling effect [27,28], 1 km-ASTER-Ts has a larger minimum and a smaller maximum values, and thus a narrow range compared to the 90 m-ASTER-Ts. The comparison of 1 km-ASTER-Ts with MODIS-Ts reveals that Ts observed by ASTER is a little larger than that by MODIS ( Figure  3a). The root mean square error (RMSE) of the difference between 1 km-ASTER-Ts and MODIS-Ts is 3.61 °C. From Table 1, only the minimum 1 km-ASTER-Ts is close to the minimum MODIS-Ts, and other statistical items show larger differences between 1 km-ASTER-Ts and MODIS-Ts. This is caused by the difference of their respective retrieval algorithms [22][23]27]. In order to make MODIS-Ts comparable to 1 km-ASTER-Ts, 1 km-ASTER-Ts was normalized based on their relationship in Figure 3a: where 1km-ASTER-Ts-N is the normalized 1 km-ASTER-Ts. The RMSE of difference between the normalized 1 km-ASTER-Ts and MODIS-Ts is reduced to 0.98 °C (Figure 3b). The 90 m-ASTER-Ts was also normalized using equation (2) to remove the effects caused by the difference between their respective algorithms, and the spatial variability and scaling issues.   Figure 4. Firstly, pure vegetation and soil pixels were identified using their respective 15 m-ASTER-NDVI thresholds, and then were resampled from 15 m to 90 m. Based on a histogram analysis and a viewing-identification on 15 m-ASTER-NDVI, it was found that NDVI of pure soil pixels was less than 0.20, while NDVI of pure vegetation pixels was greater than 0.7. These values are consistent with those proposed in a previous study [29]. Then, Ts of pure vegetation and soil pixels with 90 m resolution were averaged within a 1 km × 1 km pixel.
In order to further compare the VI-Ts diagrams defined respectively using both MODIS pixel and subpixel information, daily MODIS-Ts and MODIS-NDVI on other days through 2003 were also collected. Although only 13 cloud-free (100%) datasets were collected over our study area, these datasets could still show seasonal variations of the VI-Ts diagram ( Table 2). Here, MODIS-NDVI was calculated using the MOD09 reflectance product. Where i, j <= 121.
If i or j <1, null is set as the pixel Surface temperature of pure soil within a pixel with 1 km resolution

Method for Estimating the Component Surface Temperatures of Vegetation and Soil
It is well known that environmental factors such as land surface moisture have a characteristic of the spatial autocorrelation [16][17]. Zhang et al. developed a two-layer remote sensing model for land surface flux estimations under the assumption that land surface moisture status of soil and vegetation within a pixel are the same or similar because of the same water supply sources (e.g.; from rain, irrigation, or upgoing ground water) [21]. In this study, the assumption was enlarged from one pixel to a moving window (3 pixels × 3 pixels). Thus, all these 9 pixels can be assumed to be on an isoline of land surface moisture in the VI-Ts diagram [4][5][6][7][8][9]. Therefore, T soil and T veg within this given pixel can be obtained by extending this isoline to the vertical lines of f veg = 0 and f veg = 1 ( Figure 5 where s is the slope of the isoline of land surface moisture, and f veg is vegetation cover fraction. It can be calculated using a second-order scaled NDVI [30][31][32] where NDVI max and NDVI min are the NDVIs of full vegetation (f veg = 1) and bare soil (f veg = 0). Based on the analysis of MODIS-NDVI data during the whole growing lifecycle of winter wheat in our study area in 2003, NDVI max was given as 0.85, and NDVI min was given as 0.20. The isoline in Figure 5 can be obtained from 9 pixels by a linear regression. As for m × n Ts image, T soil and T veg images in the size of (m-2) × (n-2) can be obtained using equations (3) and (4). Then, the size of sampling window should be determined to define the VI-Ts diagram. In this study, the size of study area is 30 km × 30 km (30 pixels × 30 pixels for 1 km-MODIS data). Hence the 28 pixels × 28 pixels area is considered as a sampling window. Within this sampling window, dry and wet points respectively correspond to soil surface with minimum moisture and vegetation surface with maximum moisture, therefore, the maximum T soil and minimum T veg can be used to determine dry and wet points, then the VI-Ts diagram can be defined only using dry and wet points. In this study, a 2 nd scaled NDVI was used to define the VI-Ts diagram (f veg -Ts) instead of NDVI. ASTER-Ts and ASTER-NDVI data over our study area on May 9 2003 were used to compare the NDVI-Ts, 1 st order NDVI-Ts and 2 nd order NDVI-Ts diagrams. Results show that the use of the 2 nd order NDVI instead of NDVI in a VI-Ts diagram does not change the triangular form.

Obtaining MODIS-T veg and MODIS-T soil using the proposed approach
MODIS-T veg and MODIS-T soil were obtained using the approach proposed in Section 3 across the whole year of 2003. The accuracies of MODIS-T veg and MODIS-T soil depend on the isoline slope of land surface moisture based on the linear regression from 9 neighboring pixels. R 2 is a measure of goodness-of-fit of linear regression. From Table 3, it is found that the means of R 2 over the 28 pixels × 28 pixels area are all larger than 0.5 across the whole year of 2003. Ts of pure soil and vegetation were also obtained from the normalized 90 m-ASTER-Ts according to the procedure shown in Figure 4. From Table 4, it is shown that MODIS-T veg and MODIS-T soil agree well with 90 m-ASTER-T veg -N and 90 m-ASTER-T soil -N, respectively. The differences of statistics between MODIS-T veg and 90 m-ASTER-T veg -N are less than 1 °C. Compared with 90 m-ASTER-T soil -N, however, MODIS-T soil has a narrower range because of the inherent bias between MODIS-Ts and ASTER-Ts although ASTER-Ts has been normalized. These accuracies of MODIS-T veg and MODIS-T soil estimations by the proposed approach are similar to those reported in [14,33], especially for the estimation of MODIS-T veg .

Validation using ground data across the whole year of 2003
It is hard to measure T soil and T veg at subpixel scales in fields, so it is difficult to validate retrieved T soil and T veg pixel by pixel. Another strategy was used to validate MODIS-T veg and MODIS-T soil obtained using the proposed approach in this study. The minimum T veg is used to determine wet point, and the maximum T soil is used to determine dry point in this study. Therefore, Ts_wet (wet point) and Ts_dry (dry point) from ground data (Section 2.1) can be reversely used to validate MODIS-T veg and MODIS-T soil to a certain extent. From Figure 6, results show that the minimum MODIS-T veg and maximum MODIS-T soil are respectively close to the ground-based Ts_wet and Ts_dry across the whole year of 2003. Both R 2 in Figure 6 are very close to 1. Their respective RMSEs are 1.28 °C and 1.16 °C.

Comparisons based on the Normalized 90 m-ASTER-Ts on May 9 2003
Three triangular VI-Ts diagrams defined by different methods are shown in Figure 7. In Figure 7a, the maximum 90 m-ASTER-T soil -N and minimum 90 m-ASTER-T veg -N were used to determine dry and wet points, then to define a triangular VI-Ts diagram on May 9, 2003. This VI-Ts diagram based on the normalized 90 m-ASTER-Ts can be considered as a "true" VI-Ts diagram in this study because of its fine spatial resolution. The triangular VI-Ts diagram in Figure 7b was defined based on the information within MODIS pixels using the proposed method in this study. In Figure 7c, the triangular VI-Ts diagram was defined using the traditional method based on the information of MODIS pixels. The traditional method is a combination of two previous automatic methods [5,9,11]. Firstly, the pixels with maximum temperature within each small VI interval are selected in the VI-Ts scatter plot. Then, some selected pixels far from the dry edge are excluded. Finally, the dry edge is defined by a linear regression on the available selected pixels.
The VI-Ts diagram in Figure 7b is similar to that in Figure 7a, but the VI-Ts diagram in Figure 7c is far from that in Figure 7a. The wet point in Figure 7b is very close to that in Figure 7a because the difference between the minimum 90 m-ASTER-T veg -N and the minimum MODIS-T veg is less than 0.2 °C. However, the dry point in Figure 7b is a little lower than that in Figure 7a because the Ts_dry in  Figures 7a and 7b resulted from the Ts bias between MODIS and ASTER is ignored, the VI-Ts diagrams are very close between two figures. However, the Ts_dry in Figure 7c (26.84 °C) is a little lower, and the Ts_wet (21.56 °C) is a little higher by comparing with Figure 7a. This results in a VI-Ts triangle far from the "true" VI-Ts diagram in Figure 7a. The reason of higher Ts_wet and lower Ts_dry in Figure 7c is due to few pixels for determining the dry edge within the range between 0.5 and 1.0 of vegetation cover fraction. This weakness cannot be avoided by the traditional method because it depends on both the full ranges of vegetation cover fraction and land surface moisture. From the comparisons in Figure 7, it is indicated that the proposed method has a similar capability to define the VI-Ts diagram using the coarse-resolution MODIS data just like using the fine-resolution ASTER data.  The proposed and traditional methods were also tested using other 12 MODIS datasets on cloudless days through 2003. Results are shown in Figure 6 and Figure 8 (including results on May 9, 2003), and the VI-Ts diagrams are shown in Figure 9. Both RMSEs of Ts_dry and Ts_wet by the proposed method are less than 1.28 °C, but for the traditional method, the RMSE of Ts_wet is 6.11 °C, and the RMSE of Ts_dry is 2.40 °C. From Figure 9, it is obviously shown that there are larger errors of Ts_dry by the traditional method on Jul. 26 (5.64 °C for absolute bias to the ground-based Ts_dry) and Sep. 21 (4.81 °C for absolute bias to the ground-based Ts_dry). The period from July to September is the rainy season in the NCP. There were 63.5 mm and 11 mm rainfalls in the half month just before Jul. 26 and Sep. 21 in 2003, respectively ( Figure 10). Both soil and vegetation had been watered before the EOS-Terra satellite overpassed on Jul. 26 and Sep. 21, so the Ts of bare soil pixels were close to the Ts of vegetation pixels, and the ranges of Ts were only 4.5 °C and 6.68 °C over the study area, respectively. Therefore, it was very difficult to find the true dry point using the traditional method at the 1-km pixel scale in the rainy season, although the concrete, asphalt or some bare soil surfaces within pixels had become dry after the sunrise. In contrast, the proposed method could obtain the dry points with good accuracies based on subpixel information in the rainy season. The absolute biases to the ground-based Ts_dry were reduced to 0.88 °C and 2.09 °C on Jul. 26 and Sep. 21, respectively.  Figure 9. The VI-Ts diagrams defined by the proposed (black) and traditional (red) methods using MODIS data across 2003. The horizontal axis is vegetation cover fraction, and the vertical axis is surface temperature in °C. From Figure 9, there are large errors of Ts_wet by the traditional method mainly on Mar. 27 (13.82 °C for absolute bias to the ground-based Ts_wet, before the re-growing of winter wheat), Jun. 28 (11.79 °C for absolute bias to the ground-based Ts_wet, after the harvest of winter wheat) and Dec. 26 (10.71 °C for absolute bias to the ground-based Ts_wet, in winter).The ranges of MODIS-NDVI are all narrower than 0.5 on these three days, so less information on Ts of vegetation pixels contributes to dry edges, and which results in dry edges and Ts_wet with large errors. However, the proposed method can obtain wet points using the minimum T veg within a sampling window even in the case of narrow range of NDVI. The absolute biases were reduced to 0.75 °C, 2.51 °C and 1.23 °C on Mar. 27, Jun. 28 and Dec. 26, respectively.

Issues about the Bias between MODIS-Ts and ASTER-Ts
MODIS-Ts and ASTER-Ts retrievals were compared by considering the spatial variability of the 11 × 11 aggregated ASTER pixel (90 m) values inside the corresponding MODIS pixel (1 km) in a previous study [27]. 2.5 °C was taken as a given threshold. If the standard deviation of radiometric temperature was larger than 2.5 °C, the MODIS and ASTER pixels were not considered, and then their results showed that selected MODIS-Ts and scaled ASTER-Ts retrievals were in a good agreement. In order to compare all MODIS-Ts and ASTER-Ts over our study area, the threshold of Ts standard deviation was not used to exclude pixels. From Figure 3a, MODIS-Ts is lower than 1 km-ASTER-Ts. The first reason comes from their respective algorithms [22,23,27], and the second reason is probably that the spatial variability and scaling issues are not considered [27,28]. In order to remove these effects, ASTER-Ts was normalized in this study. MODIS-Ts is then relatively consistent with the normalized 1 km-ASTER-Ts (Figure 3b), but MODIS-Ts is still a little lower than 1 km-ASTER-Ts-N, especially for pixels with higher Ts (see points in the upper of scatter plot in Figure 3b). This results in a small shift of dry point in Figure 7. If this shift is igonred, the VI-Ts diagrams are very close between Figure 7a and 7b.

Issues about the Approach of Obtaining T veg and T soil
In the proposed method of defining the VI-Ts diagram, the key process is to obtain T veg and T soil , and then to determine dry and wet points. In this study, T veg and T soil are obtained under the assumption of homogenous land surface moisture for neighboring pixels. Result from [18] show that soil moisture has a spatial correlation within a scale of 6,400 m × 6,400 m in an agricultural landscape. This spatial distribution of land surface moisture is mainly controlled by land cover/use [34]. If the 3 × 3 moving window in Figure 5 covers the same land cover / use, the assumption of homogenous land surface moisture for neighboring pixels is valid. In practice, the 3 × 3 window may covers the boundary between different land covers / uses, for example, the boundary between water and land surfaces, and the boundary between city and surrounding areas. When this case occurs, a land cover / use data such as MOD12 (http://modis-land.gsfc.nasa.gov/landcover.htm) is allowed to exclude pixels whose land types are different to that of the pixel (i, j) within the 3 × 3 window. When the excluded pixels are more than 3, the values of T veg and T soil within the pixel (i, j) are set as nulls. Actually, there were a lot of 3 × 3 moving windows in the images used in this study, and thus a window with the minimum soil water content can be obtained. It is considered that the surface temperature of soil within this window will be closest to the dry point defined in Figure 1.

Issues about Obtaining Wet and Dry Points
The minimum T veg (wet point) can be considered as the "true" value of vegetation without water deficit. In this study, we assume that at least vegetation without water deficit can always be found at subpixel scales. Results show that this assumption was valid in our study area. Wet points were obtained with a RMSE error of 1.28 °C through the whole year of 2003. However, true wet points cannot be easily found in some arid regions because well-watered vegetation surface may not exist even at subpixel scales. Wet points will be overestimated in this case.
However, T soil may be not the "true" value of soil but of other components such as roads, house roofs and water surface. In this study, dry point is defined as a surface without vegetation cover and evaporation. Water surface can be excluded because of its low Ts and strong evaporation, but the impervious surfaces such as roads, house roofs satisfy the definition of dry point. It is well known that impervious surfaces have a heating effect through absorbing and holding solar energy, so Ts of them are usually higher than Ts of vegetation, water and general soil. Analysis of ASTER images for examining the relationship between urban thermal features and biophysical descriptors in Indianapolis (Indiana, USA) showed that Ts of impervious surfaces in urban areas were relatively higher than that of soil surface in agricultural areas [35][36]. This was also validated using ASTER images in our study area, and results show that Ts of impervious surfaces in urban areas are 2-6 °C higher than that of soil surface in agricultural areas on May 9, 2003. Further analysis on the VI-Ts diagram from 90 m-ASTER-Ts and 90 m-ASTER-NDVI shows that some scatter points near the position of dry point in the VI-Ts diagram correspond to urban areas. Therefore, there may be several dry points in a VI-Ts diagram, especially in or near urban areas. If several dry points occur in the same VI-Ts diagram (Figure 11), the dry point defined by our proposed method must be the point (dry point 1) with the maximum Ts (e.g.; roads, house roofs), not the true dry point (dry point 3, e.g.; bare soil). In this case, the true dry point will be overestimated. However, this error (RMSE of 1.16 °C for dry point estimations) from this reason is not significant because our study area is an agricultural landscape. In order to identify the true dry point in or near urban areas, the technology of spectral mixture analysis [35][36] is suggested to distinguish bare soil and impervious surfaces (e.g.; roads and house roofs) in further studies. In practices, alternatively, dry point can be determined using the average of several T soil near the maximum T soil in order to reduce the effects from impervious surfaces. Sometimes, dry point can also be underestimated in some humid regions because bare soil surface may not exist even at subpixel scales. False dry edges

Issues about the Definition of the VI-Ts Diagram
To define an ideal VI-Ts diagram, the traditional method requires the full and continuous ranges of land surface moisture and VI [5,9,11]. It needs enough heterogeneity of land surface moisture and vegetation cover within a limited sampling window. However, this condition is difficult to meet at coarse-resolution scales such as 1 km scale for MODIS data. For example, narrow ranges of VI resulted in large errors of Ts_wet by the traditional method on Mar. 27, Jun. 28 and Dec. 26 in 2003 in this study. The advantage of our proposed method is independent of the land surface heterogeneity, but only depends on two extreme land surface conditions at subpixel scales, dry bare soil and wellwatered vegetation. If dry and wet points can be found within the sampling window, a VI-Ts diagram will be easily defined. This is also a limitation of our proposed method. The proposed method will be helpless to define a VI-Ts diagram if the land cover is entirely homogenous within the sampling window, such as continuous bare soil (e.g.; North African and Mid-Asian deserts) and vegetation (e.g.; tropical and boreal forests) regions. In these regions, all pixels within the sample window will cluster near dry or wet points (no wet point in desert regions, no dry point in forest regions). In applications, if such a cluster occurs in the VI-Ts diagram, it will also contribute to studies about land surface energy balance. The cluster near dry point indicates almost no latent heat partitioned from land surface available energy, while the cluster near wet point indicates almost no sensible heat partitioned from land surface available energy. Even if the land surface is a mixture of bare soil and vegetation, the sampling window should still cover dry bare soil. Vegetation such as forest, grass and crop can be easily found from a mixed land surface. Dry bare soil is also easily found in arid or semiarid areas. However, if the land surface is fully wet across the entire sampling window, such as in humid areas, it is impossible to find dry soil even at subpixel scale. In this case, the VI-Ts scatter plot looks like a horizontal line (no dry point in humid areas), and Ts of all pixels within the sampling window are close to air temperature [10]. Almost all available energy of land surface is transferred to latent heat in such a case. In brief, our proposed method cannot define a triangular VI-Ts diagram without dry and wet points. Our study area is a typical semiarid agricultural landscape (a mixture of bare soil and crop), so our proposed method can obtain dry and wet points with good accuracies through the whole year in this study.
Several previous studies reported that a positive relationship of VI-Ts may occur in lowtemperature or low-radiation regions [37][38][39][40]. If the slope of dry edge in a VI-Ts diagram is positive, the VI-Ts diagram will not contribute to the estimations of land surface moisture status or evapotranspiration from land surface because there is almost no evapotranspiration in the cases of lowtemperature or low-radiation. Thus the positive relationship of VI-Ts was seldom mentioned in applications.

Conclusions
The traditional method cannot always define a reasonable VI-Ts diagram within a limited sampling window. It is invalid in two cases of rainy season and narrow range of VI. In this study, therefore, a new method was proposed to define a VI-Ts diagram using dry and wet points from the subpixel vegetation and soil information. In this method, a new approach was proposed to obtain T veg and T soil within a given pixel by means of neighboring pixels. The VI-Ts diagram defined by the proposed method is close to that defined using the ASTER-Ts of pure vegetation and soil. The proposed method can obtain "true" wet and dry points with respective average accuracies of 1.28 °C and 1.16 °C across the whole year of 2003 in our study area. This indicates that our proposed method can define "true" VI-Ts diagrams in our study area across the whole year, even in two cases of rainy season and narrow range of vegetation index. In applications, two limitations of the proposed method should be noted. A triangular VI-Ts diagram cannot be obtained if either dry point or wet point is absent. The other is that some false dry points may be involved in or near urban areas.