Improving Estimation of Evapotranspiration under Water-Limited Conditions Based on SEBS and MODIS Data in Arid Regions

This study proposes a method for improving the estimation of surface turbulent fluxes in surface energy balance system (SEBS) model under water stress conditions using MODIS data. The normalized difference water index (NDWI) as an indicator of water stress is integrated into SEBS. To investigate the feasibility of the new approach, the desert-oasis region in the middle reaches of the Heihe River Basin (HRB) is selected as the study area. The proposed model is calibrated with meteorological and flux data over 2008–2011 at the Yingke station and is verified with data from 16 stations of the Heihe Watershed Allied Telemetry Experimental Research (HiWATER) project in 2012. The results show that soil moisture significantly affects evapotranspiration (ET) under water stress conditions in the study area. Adding the NDWI in SEBS can significantly improve the estimations of surface turbulent fluxes in water-limited regions, especially for spare vegetation cover area. The daily ET maps generated by the new model also show improvements in drylands with low ET values. This study demonstrates that integrating the NDWI into SEBS as an indicator of water stress is an effective way to improve the assessment of the regional ET in semi-arid and arid regions.


Introduction
Evapotranspiration (ET), which is the sum of soil surface evaporation and plant tissue transpiration, is an important energy exchange process of the soil-vegetation-atmosphere system that affects the redistribution of terrestrial surface radiation and precipitation [1].Accurate estimations of land surface ET are important for studies on water resources and agricultural management, flood and drought monitoring, weather forecasting, and climate change.Traditional methods for measuring ET, such as the Bowen ratio, eddy correlation systems, aerodynamics, and soil water balance, can only obtain ET over small and homogeneous regions, and fluxes at large spatial scales are difficult to estimate [2].For operational applications in agriculture, such as irrigation scheduling, water resource management, and drought monitoring, accurate regional variations in the surface fluxes, particularly ET, are difficult to determine.Satellite remote sensing (RS) is an attractive tool for obtaining land surface information related to the energy and water balances, such as the albedo, land surface temperature (LST), leaf area index (LAI), and surface soil moisture.Therefore, RS data are ideally suited for deriving spatially-continuous fields of ET [2].Over the last few decades, several methods have been proposed to estimate ET, which can be divided into three main categories: empirical approaches, physical models, and data assimilation [3].Empirical approaches often apply statistical regression to estimate ET using the satellite vegetation index, LST, LAI, and other meteorological data [4][5][6][7].Physical models can be classified into two types: surface conductance-based models and aerodynamic resistance-energy balance models [7].Surface conductance-based models were developed from the Penman-Monteith (PM) combination equation, where the required surface conductance is estimated from the remotely sensed LAI using a biophysical model [8][9][10].To overcome conductance parameterization problems, the Priestley-Taylor (PT) equation was developed, which simplified the PM equation by replacing the surface and aerodynamic resistance terms with an empirical multiplier α [11][12][13].Aerodynamic resistance-energy balance models are based on the bulk resistance equation for heat transfer, the difference between the land surface temperature and the air temperature, and the aerodynamic resistance to turbulence heat transport [14].The instantaneous latent heat flux can be estimated as a residual of the land surface energy balance equation.Depending on whether the soil and canopy are considered separately or together, these models can be divided into single-source and two-source models.Single-source models consider the land surface as a single system, while two-source models treat the land surface as a combination of soil and a canopy with variable turbulence fluxes.Examples of these models include the Surface Energy Balance Algorithm for Land (SEBAL) [15], the Simplified Surface Energy Balance Index (S-SEBI) [16], the Surface Energy Balance System (SEBS) [17], the satellite-based energy balance for Mapping Evapotranspiration with Internalized Calibration (METRIC) [18], and the Two-Source Energy Balance Model (TSEB) [19], which have been thoroughly discussed in a previous study [20].Data assimilation is such a methods that can be used to combine land surface models and observations to accurately estimate ET on a continuous spatial and temporal scale [21][22][23].Comparisons of different methods for estimating land surface fluxes at the local and global scales have demonstrated that the models have different degrees of success [24,25].
Evapotranspiration is controlled by the availability of moisture and energy fluxes at the surface [24].Under water-limited conditions, the available soil water limits ET.However, almost all empirical and physical models, including SEBS, lack a variable that directly estimates water stress.Several recent studies have integrated water stress into estimates of land surface evapotranspiration [14,31,[33][34][35].It has been found that water stress would underestimate the sensible heat flux and overestimate ET in SEBS models in semi-arid and arid regions [27,31,36].To deal with this limitation, Gokmen et al. [31] explicitly integrated soil moisture information as a water stress index through a modified definition of kB ´1, which is the parameter that summarizes the excess aerodynamic resistance to heat transfer compared with momentum transfer.However, the surface turbulent fluxes and ET are difficult to obtain at high spatial resolutions using this method due to the lack of soil moisture observations.To improve the estimates of surface turbulent flux and ET in SEBS under water stress conditions, we have proposed an alternative way to improve estimates of surface turbulent fluxes in SEBS under water stress conditions using MODIS data [37].In our previous study, the normalized difference water index (NDWI) is integrated into SEBS as an indicator of water stress based on linear, exponential, and S-curve equations, which can reduce kB ´1 as water stress increases and can further decrease the latent heat flux and ET as the sensible heat flux increases.One site (Yingke station from the Watershed Allied Telemetry Experimental Research (WATER) project during 2008-2009 [38]) was used to calibrate and validate with meteorological and EC flux observations in a previous study.The results show that the proposed scheme can significantly improve the land surface flux estimation in the Yingke station under water stress condition and the performance of revised SEBS based on linear, exponential, and S-curve equations are same.Wu et al. [39] also integrated the soil moisture information into SEBS for improving irrigated crop water consumption estimation.In order to ensure that the proposed scheme can be applied under a wide range of vegetation coverage and environmental conditions, in this study, the proposed scheme (only the S-curve equation) is calibrated with observations at the Yingke station during 2008-2011 and is validated by observations from 16 stations of the Heihe Watershed Allied Telemetry Experimental Research (HiWATER) project in 2012 [40].Only the S-curve equation is calibrated and validated on account of the conceptual framework for defining soil moisture regimes and the corresponding evapotranspiration regimes proposed by Seneviratne et al. [41] and the sigmoid function used by Gokmen et al. [31].Furthermore, the modified SEBS model (hereafter referred to as SEBS-NDWI) was used to calculate the regional land surface fluxes and daily evapotranspiration in the study area.

Study Area
The Heihe River Basin (HRB) is the second largest inland river basin in China, covering an area of approximately 143,000 km 2 .From upstream to downstream, the basin features glaciers, frozen soil, alpine meadows, forests, irrigated crops, riparian ecosystems, and deserts.The desert-oasis zone in the middle reach of the HRB, which is located between 100 ˝14 1 E-100 ˝32 1 E and 38 ˝42 1 N-39 ˝01 N, is chosen as the study area (Figure 1).The central part of the study area is a typical irrigated crop ecosystem, in which the main crops are maize and wheat.The region surrounding the oasis is mainly sandy desert, desert steppe, and the Gobi desert.The annual average temperature is 2.8-7.6 ˝C, and there are approximately 3000 to 4000 h of sunshine per year.The annual precipitation is approximately 100-250 mm, but the potential evaporation is as high as 1200-1800 mm.Comprehensive eco-hydrological experiments, such as the WATER [38] and HiWATER [40] projects, have successfully been performed in this region.Many satellite, airborne, and ground observations that are suitable for ET studies have been collected.

Field Measurements
The field measurements include meteorological data and flux data that were obtained during the WATER and HiWATER experiments.These stations are described in Table 1.
During the WATER experiment, the Yingke station was installed in the central area of the oasis, which is located within typical irrigated farmland of maize and wheat (Figure 1).An automatic meteorological station (AMS), eddy covariance (EC) system, and soil moisture and temperature measurement system (SMTMS) were installed at the station.A set of continuous observations were obtained from November 2007 to December 2011.In this study, the observations from 1 March to 31 October of 2008-2011 are used to calibrate the model parameters of the proposed SEBS-NDWI model.ecosystem, in which the main crops are maize and wheat.The region surrounding the oasis is mainly sandy desert, desert steppe, and the Gobi desert.The annual average temperature is 2.8-7.6 °C, and there are approximately 3000 to 4000 h of sunshine per year.The annual precipitation is approximately 100-250 mm, but the potential evaporation is as high as 1200-1800 mm.Comprehensive eco-hydrological experiments, such as the WATER [38] and HiWATER [40] projects, have successfully been performed in this region.Many satellite, airborne, and ground observations that are suitable for ET studies have been collected.During the HiWATER experiment, a multi-scale observational experiment on evapotranspiration over heterogeneous land surfaces (MUSOEXE) was conducted from May to September 2012 [40,42].Four AMSs were installed around the oasis (i.e., sandy desert, desert steppe, Gobi desert, and wetland stations); each station was supplemented with an EC system (Figure 1).An observation matrix composed of 17 EC towers and ordinary AMSs were located within the oasis (i.e., the kernel experimental area covered an area of 5.5 km ˆ5.5 km, Figure 1), where the land surface is heterogeneous and is dominated by seed corn, corn interplanted with spring wheat, vegetables, orchards, and residential areas (Figure 1) [40].Due to problems with the instruments, some of the observations at stations S3, S9, S10, S16, and S21 are missing, so the remaining observations from the 16 stations from 15 June to 15 September 2012, are used to validate the proposed SEBS-NDWI model.
The meteorological data (air temperature, humidity, wind speed, pressure, shortwave radiation, and longwave radiation) used in this study were obtained from AMS. Turbulent fluxes were measured by EC system with a sampling frequency 10-Hz.The raw high-frequency data were processed using EdiRe software developed by Edinburgh University, UK.Ground heat flux was Remote Sens. 2015, 7, 16795-16814 measured by soil heat flux plates buried at a depth of 5 cm and 6 cm, respectively, during the WATER experiment and HiWATER experiment.The meteorological data were recorded at 30 min interval during WATER experiment.The AMS in the HiWATER experiment recorded meteorological data at 10 min intervals, which were processed in 30 min interval.The turbulent fluxes from the EC system were processed at 30 min interval.The details of data processing can prefer to Liu et al. [43] and Liu et al. [44].In this study, we excluded the data with energy balance closure of less than 0.75.In addition, the Bowen ratio closure method was used to close the energy balance for turbulent fluxes from EC system [45].

MODIS Land Data Products
The MODIS land data products used in the study include surface reflectance in bands 1-7 and LST.These data are obtained from the online data pool at the NASA Land Processes Distributed Active Archive Center (LP DAAC), United States Geological Survey (USGS)/Earth Resources Observation and Science (EROS) Center, Sioux Falls, South Dakota.We choose the latest MOD09GA and MOD11A1 data, which provide bands 1-7 daily-gridded L2G Global 500 m reflectance values and the LST/E Daily L3 Global 1 km Grid product, respectively.We then use the MODIS Reprojection Tool (MRT) to reproject the original sinusoidal projection to a UTM projection with a WGS84 datum and a 500 m resolution.

Normalized Difference Water Index
The NDWI is defined as follows where ρ NIR and ρ SWIR represent the reflectances of the near-infrared and shortwave infrared bands, respectively, which correspond to band 2 (841-876 nm) and band 7 (2105-2135 nm) in MODIS images [46].A positive correlation has been found between the NDWI and vegetation water content [47], so the NDWI is a useful indictor of vegetation water stress.

Method
SEBS was developed to estimate surface turbulent fluxes and daily evapotranspiration using remote sensing data [17].
The details of the original SEBS formulation can be found in references [17,48].In this section, we first provide a short theoretical overview of SEBS and then describe how to integrate the NDWI into SEBS to address water stress in sparsely-vegetated semi-arid and arid regions.

Surface Energy Balance Terms
The surface energy balance is usually written as where R n is the net radiation, G is the soil heat flux, H is the sensible heat flux, and λE is the latent heat flux.In SEBS, λE is commonly calculated as the residual of the surface energy balance: The net radiation (R n ) is given by: R n " p1 ´αqR swd `εR lwd ´εδT 4 s where R swd and R lwd are the downward solar radiation and longwave radiation, respectively, which obtained from AMSs; α and ε are the albedo and emissivity of the land surface, respectively, δ is the Stefan-Boltzmann constant, and T s is the surface temperature.
The soil heat flux is usually controlled by R n and is affected by the fraction of canopy cover, which is parameterized as: G " R n rΓ c `p1 ´fc qpΓ s ´Γc qs (5) where f c is the fraction of canopy cover, Γ c " 0.05 for a full vegetation canopy, and Γ s " 0.315 for bare soil.

Calculation of Sensible Heat Flux
The sensible heat flux is obtained iteratively by solving the set of nonlinear equations shown below, which were derived in part from the Monin-Obukhov similarity theory: where Z is the height above the surface, µ is the wind speed, u ˚is the friction velocity, k = 0.4 is von Karman's constant, d 0 is the zero plane displacement, ρ is the air density, C p is the specific heat of air at a constant pressure, z 0m and z 0h are the roughness heights for momentum and heat transfer, respectively, ψ m and ψ h are the stability correction functions for momentum and heat transfer, respectively, θ a is the potential temperature at the height z, θ 0 is the potential temperature at the surface, L is the Obukhov length, θ v is the virtual temperature near the surface, and g is the acceleration due to gravity.

Calculation of Roughness Length for Heat Transfer
The roughness length for heat transfer can be derived from: where B ´1 is the inverse Stanton number, which is a dimensionless heat transfer coefficient.SEBS uses a physically-based model to calculate kB ´1 [48], which is expressed as: where f s " 1 ´fc is the fraction of soil coverage; C d is the drag coefficient of the foliage elements, which is assumed to have a value of 0.2; C t is the heat transfer coefficient of the leaf, which has values of 0.005N ď C t ď 0.075N (N is the number of sides of a leaf that participate in the heat exchange); h is the height of the vegetation canopy; n ec is the within-canopy wind speed profile extinction coefficient; uphq is the horizontal wind speed at the canopy height; C t is the heat transfer coefficient of the soil; and kB ´1 s is the value for the bare soil surface given by Brutsaert [49]: where Re ˚is the roughness Reynolds number.

Calculation of Daily ET
The daily ET is more applicable than the instantaneous ET and turbulent flux in hydrology and meteorology.Several methods have been developed to scale daily ET from instantaneous ET at the time of the satellite overpass.These methods are based on the assumption of diurnal self-preservation of the evaporative fraction (EF) [50][51][52][53].The EF is defined as the ratio of the actual evaporation to the available energy (i.e., the difference between R n and G): where Λ is the evaporative fraction.
In this study, we use the following scheme to calculate the daily ET.First, the evaporative fraction is evaluated according to SEBS [17]: where λE wet is the latent heat flux at the wet limit, and Λ r is the relative evaporative fraction: where H dry is the sensible heat flux at the dry limit.The estimation of λE wet and H dry is described in [17].
Second, the daily net radiation is derived from: where R n_daily is the daily net radiation, and the daily net shortwave radiation R ns_daily [54]: where W d is the daily solar radiance; a s , b s and the fraction of insolation n{N (0.7 in June, 0.67 in July, 0.65 in August, and 0.7 in September are used in this study [55]) are empirical coefficients used to perform the atmospheric correction; p 1 ρ 0 q 2 is the eccentricity correction; I 0 is the solar constant; ω 0 is the sunrise angle; and ϕ and δ are the latitude and the solar declination, respectively.Brunt [56] introduced an empirical parameterization of the daily net longwave radiation R nl_daily that is based on the air temperature, humidity and the fraction of insolation: where e d is the actual vapor pressure, and T a is the air temperature, which is obtained from AWSs.Depending on the climatic conditions, the empirical coefficients of a, b, c, and d can vary.In a continental climate, a, b, c, and d equal 0.39, 0.058, 0.1, and 0.9, respectively.The daily ET is finally given by: where ET daily is the daily ET (mm), and λ is the latent heat of vaporization.

Parameterization of Land Surface Parameters
The necessary parameters of SEBS can be divided into two categories: (1) surface parameters, which include the surface temperature T s , the surface albedo α, the surface emissivity ε, the fraction of canopy cover f c , the leaf area index LAI, and the height of the vegetation canopy h, based on MODIS data; and (2) the meteorological parameters, which include the air temperature T a , the barometric pressure p, the wind speed u, the specific humidity, the downward solar radiation R swd , and the downward longwave radiation R lwd , measured by the AMS.
The f c , LAI, and h are useful indicators for monitoring the condition of the land surface and plant health and are expressed as functions of the NDV I [57]: h " h min `ph max ´hmin q NDV I ´NDV I min NDV I max ´NDVI min ( 22) where NDVI min and NDVI max are the maximum and minimum NDV I values in the study area, respectively (0.05 and 0.87 were used in this study); h min and h max are the h values for bare soil and full vegetation, respectively (0.0012 m and 2 m were used in this study); and ρ NIR and represent the reflectance of the near-infrared and the red bands, respectively, which corresponds to band 2 (841-876 nm) and band 1 (620-670 nm) in MODIS images.The estimation approach of f c , LAI, and h were validated with ground measured data from HiWATER experiment in this study area.The determination coefficient (R 2 ) between estimated parameters (i.e., f c , LAI, and h) and measured data are 0.71, 0.67, and 0.73, respectively, which indicated that these functions can be used to estimate f c , LAI, and h in this region.
The surface albedo is used to estimate the amount of incident solar radiation that is reflected back to space; it greatly influences the surface energy budget and partitioning.Liang [58] proposed the narrowband-broadband conversion method for surface albedo using MODIS data: α=0.160ρ 1 +0.291ρ 2 +0.243ρ 3 +0.116ρ 4 +0.112ρ 5 +0.081ρ 7 -0.0015(24) where ρ i (i = 1-5, 7) represents the corresponding MODIS band surface reflectance values.The surface emissivity is calculated by weighing the vegetation emissivity ε c and soil emissivity ε s according to their proportions in a pixel [26,59,60]: where ε c = 0.985 for full vegetation and ε s = 0.96 for bare soil, and ă dε ą is a regional representative value that ranges from 0.00 to 0.02 (0.02 in this case).

kB ´1 Parameterization Scheme
In SEBS, provided that z 0m remains constant, a decrease in kB ´1 leads to decreased latent heat flux and ET as the sensible heat flux increases.Gokmen et al. [31] proposed an explicit integration of soil moisture information as a water stress index through a modified definition of kB ´1 in SEBS, which affects the exchange of sensible and latent heat between the canopy and the atmosphere.However, the surface turbulent fluxes and ET are difficult to obtain at high spatial resolutions using the method due to the lack of soil moisture observations.The NDWI was proposed by Gao [46] to detect the moisture conditions of vegetation canopies over large areas [4,47].Significant positive correlations have also been found between the NDWI and EF which is shown in Section 4.1.Therefore, the kB ´1 parameterization scheme [31] is adopted but the water stress index is replaced by NDWI in this study.The specific method used to modify the value of kB ´1 in SEBS is expressed as: where kB ´1 NDW I is the modified kB ´1 value, which indicates the level of water stress; f is a scaling factor, which is represented by a sigmoid function; and a, b, and c are the coefficients of the sigmoid function, which can be derived from an optimization (e.g., the Shuffled Complex Evolution Algorithm-University of Arizona (SCE-UA) [61] by reducing the error between the EC observations and the SEBS estimates of sensible heat flux.

Influence of Soil Moisture on Evaporative Fraction
Many authors have examined the influence of soil water availability on evapotranspiration by determining the relationship between soil moisture and evapotranspiration in different regions and soil types, land cover and scales.Positive correlations have been used when the soil water availability is insufficient [62][63][64].Consistent with earlier studies, a positive correlation was found between the near-surface soil moisture (SM) and the EF at the Yingke station during the WATER experiment in 2008-2011 (Figure 2a) and at sixteen stations of the HiWATER experiment in 2012 (Figure 2b).The SM is less than or equal to 0.3 m 3 /m 3 .The EF increases with an increase in SM until the SM is greater than 0.3 m 3 /m 3 .When the soil water availability is sufficiently high (i.e., SM is greater than 0.3 m 3 /m 3 ), ET is controlled by the available radiant energy, and the effect of soil moisture on ET is not significant.Remote Sens. 2015, 7 page-page In this case, we propose an explicit integration of NDWI as a water stress index into SEBS through a modified definition of 1 kB − .The scatter plots between the NDWI and EF at the Yingke station over 2008-2011 (Figure 3a) and at sixteen HiWATER stations in 2012 (Figure 3b) indicate significant positive correlations between the NDWI and EF, with determination coefficients of 0.62 and 0.75, respectively.In this case, we propose an explicit integration of NDWI as a water stress index into SEBS through a modified definition of kB ´1.The scatter plots between the NDWI and EF at the Yingke station over 2008-2011 (Figure 3a) and at sixteen HiWATER stations in 2012 (Figure 3b) indicate significant positive correlations between the NDWI and EF, with determination coefficients of 0.62 and 0.75, respectively.10 a 5 cm depth, and EF is the evaporative fraction of the clear-sky daytime averages of the EC measurements.
In this case, we propose an explicit integration of NDWI as a water stress index into SEBS through a modified definition of 1 kB − .The scatter plots between the NDWI and EF at the Yingke station over 2008-2011 (Figure 3a) and at sixteen HiWATER stations in 2012 (Figure 3b) indicate significant positive correlations between the NDWI and EF, with determination coefficients of 0.62 and 0.75, respectively.

Calibration of SEBS-NDWI
Generally, the spatial representativeness of EC is hundreds of meters, which is affected by the instruments heights of EC [3].In this study, estimations of the single pixel where the EC station located is taken to be compared with measurements from EC.To obtain the coefficients of the sigmoid function (Equation ( 27)), we calibrated SEBS-NDWI using EC flux observations at the Yingke station (2008-2011) according to the scheme proposed in Section 3.2.The coefficients of a , b , and c were −0.47, 0, and 8.97, respectively.During early March to early May, the Yingke station is near bare because maize is usually cultivated in later April and the vegetative emergence of maize begin in

Calibration of SEBS-NDWI
Generally, the spatial representativeness of EC is hundreds of meters, which is affected by the instruments heights of EC [3].In this study, estimations of the single pixel where the EC station located is taken to be compared with measurements from EC.To obtain the coefficients of the sigmoid function (Equation ( 27)), we calibrated SEBS-NDWI using EC flux observations at the Yingke station (2008-2011) according to the scheme proposed in Section 3.2.The coefficients of a, b, and c were ´0.47, 0, and 8.97, respectively.During early March to early May, the Yingke station is near bare because maize is usually cultivated in later April and the vegetative emergence of maize begin in early May.With the growth of maize, the surface is covered by vegetation gradually.After July, the maize grow up and enter into reproductive stage, the surface is covered with full vegetation.Maize harvesting in later September leaded to the Yingke station stand in spare vegetation areas.Owing to the seasonal diversity in vegetation coverage and environment conditions, the coefficients obtained by using Yingke Station are reasonable.
Furthermore, the observed surface turbulent fluxes from the EC system were compared with the modeled values of H and λE from the original SEBS (Figure 4a,c) and the SEBS-NDWI (Figure 4b,d).Consistent with earlier studies [27,31,33], the original SEBS clearly underestimated H (Figure 4a), particularly for low values of the NDWI; therefore, a shortage of soil water and a lack of vegetation cover were present.Underestimation of H leads to an overestimation of λE (Figure 4c).The estimates of H from the NDWI-integrated SEBS (H SEBS´NDW I ) are closer to the EC observations (H o ) for high values of H o compared with the estimates of H from the original SEBS (H SEBS ) (Figure 4b).Therefore, more EC and SEBS-NDWI point pairs are located near the 1:1 line in the scatter plot of λE (Figure 4d).The proposed scheme improves the underestimation of H when the NDWI is less than or equal to 0.28 (the median of the NDWI); it also works well when the NDWI is greater than 0.28 (Figure 4a,b).) (Figure 4b).Therefore, more EC and SEBS-NDWI point pairs are located near the 1:1 line in the scatter plot of E λ (Figure 4d).The proposed scheme improves the underestimation of H when the NDWI is less than or equal to 0.28 (the median of the NDWI); it also works well when the NDWI is greater than 0.28 (Figure 4a,b).To assess the effect of integrating the NDWI into SEBS, the RMSE and BIAS values for SEBS and SEBS-NDWI are compared based on the EC-observed surface turbulent fluxes (Table 2).The results show an overall improvement in the modeled values of H and E λ , with a 10 W/m 2 reduction in the RMSE and a 20 W/m 2 reduction in the BIAS when the NDWI is incorporated into SEBS.A significant improvement, with a large reduction in the RMSE and a 50% reduction in the BIAS, occurs when the NDWI is less than or equal to 0.28.When the NDWI is greater than 0.28, the RMSEs between the SEBS-estimated surface turbulent fluxes and the SEBS-NDWI-estimated surface turbulent fluxes are mostly less than 8 W/m 2 , which indicates that SEBS-NDWI results are similar to those of SEBS.
The integration of NDWI into SEBS can significantly improve the estimates of surface turbulent fluxes under water-limited conditions, but the estimates of surface turbulent fluxes slightly change when adequate soil water is available.Under near-bare soil and sparse vegetation areas, the NDWI may To assess the effect of integrating the NDWI into SEBS, the RMSE and BIAS values for SEBS and SEBS-NDWI are compared based on the EC-observed surface turbulent fluxes (Table 2).The results show an overall improvement in the modeled values of H and λE, with a 10 W/m 2 reduction in the RMSE and a 20 W/m 2 reduction in the BIAS when the NDWI is incorporated into SEBS.A significant improvement, with a large reduction in the RMSE and a 50% reduction in the BIAS, occurs when the NDWI is less than or equal to 0.28.When the NDWI is greater than 0.28, the RMSEs between the SEBS-estimated surface turbulent fluxes and the SEBS-NDWI-estimated surface turbulent fluxes are mostly less than 8 W/m 2 , which indicates that SEBS-NDWI results are similar to those of SEBS.The integration of NDWI into SEBS can significantly improve the estimates of surface turbulent fluxes under water-limited conditions, but the estimates of surface turbulent fluxes slightly change when adequate soil water is available.Under near-bare soil and sparse vegetation areas, the NDWI may potentially predicts soil moisture information, whereas it mainly index the moisture condition of vegetation rather than soil moisture content under the dense vegetation area.In the oasis area, the soil moisture content in the root zone of maize (depth at 10 cm-40 cm) is mostly among 0.2 m 3 /m 3 to 0.4 m 3 /m 3 , which is adequate for the consumption of maize.In other words, this area is under saturated moisture condition with higher NDWI value and vegetation moisture content and, thus, the evapotranspiration is controlled by the available energy rather than the soil moisture.This is the reason why the integration of NDWI into SEBS cannot refine the performance of the model when the pixels have high NDWI values.

Validation of SEBS-NDWI
The proposed scheme is validated with meteorological and flux data from 16 HiWATER stations in 2012.Figure 5 3).When the NDWI is less than or equal to 0.28, the BIAS is 166.6 W/m 2 for λE calculated with the original SEBS or 23.7 W/m 2 for λE calculated with SEBS-NDWI.Thus, integrating the NDWI into SEBS avoids overestimating λE (Table 3).Overall, SEBS-NDWI tends to outperform the original SEBS, particularly under water-limited conditions.

Validation of SEBS-NDWI
The proposed scheme is validated with meteorological and flux data from 16 HiWATER stations in 2012.Figure 5 compares the modeled and observed surface turbulent fluxes at the 16 stations, and Table 3 shows the error evaluation of the surface turbulent fluxes from the original SEBS and NDWIintegrated SEBS based on EC observations.The modeled values of H from the original SEBS (  3).
When the NDWI is less than or equal to 0.28, the BIAS is 166.6 W/m 2 for E λ calculated with the original SEBS or 23.7 W/m 2 for E λ calculated with SEBS-NDWI.Thus, integrating the NDWI into SEBS avoids overestimating E λ (Table 3).Overall, SEBS-NDWI tends to outperform the original SEBS, particularly under water-limited conditions.

Daily ET Mapping by SEBS-NDWI
The main aim of this study is to improve the mapping of regional surface turbulent fluxes and evapotranspiration in arid regions with SEBS by integrating the NDWI, which is an indicator of water stress.We calculated the instantaneous sensible heat flux (Figure 6), the instantaneous latent heat flux (Figure 7) and the daily evapotranspiration (Figure 8) in the study area using the original SEBS and the NDWI-integrated SEBS models on 24 June, 10 July, 2 August and 3 September 2012.
As shown in Figure 6, the spatial distribution of the sensible heat flux differs from the oasis to the desert areas.The instantaneous H modeled by the original SEBS in the oasis ranges from 0 W/m 2 to 100 W/m 2 , and the instantaneous H modeled by SEBS in the area surrounding the oasis (sandy desert, desert steppe, and the Gobi desert) is clearly greater than that in the oasis (100 W/m 2 to 200 W/m 2 ) (Figure 6a-d).Figure 5a,b show that the EC-observed values of H range from ´50 W/m 2 to 400 W/m 2 in the study area at the times of the satellite overpasses from June to September.Clearly, SEBS underestimates H when there is a shortage of soil water and a lack of vegetation coverage.In contrast, the estimates of H from SEBS-NDWI are more similar to the EC observations, which range from 150 W/m 2 to 400 W/m 2 (Figure 6e-g).The difference between the estimates of H on the four days is small, which indicates that no clear temporal variation occurs in the study area from June to September.The main aim of this study is to improve the mapping of regional surface turbulent fluxes and evapotranspiration in arid regions with SEBS by integrating the NDWI, which is an indicator of water stress.We calculated the instantaneous sensible heat flux (Figure 6), the instantaneous latent heat flux (Figure 7) and the daily evapotranspiration (Figure 8) in the study area using the original SEBS and the NDWI-integrated SEBS models on 24 June, 10 July, 2 August and 3 September 2012.
As shown in Figure 6, the spatial distribution of the sensible heat flux differs from the oasis to the desert areas.The instantaneous H modeled by the original SEBS in the oasis ranges from 0 W/m 2 to 100 W/m 2 , and the instantaneous H modeled by SEBS in the area surrounding the oasis (sandy desert, desert steppe, and the Gobi desert) is clearly greater than that in the oasis (100 W/m 2 to 200 W/m 2 ) (Figure 6a-d).Figure 5a,b show that the EC-observed values of H range from −50 W/m 2 to 400 W/m 2 in the study area at the times of the satellite overpasses from June to September.Clearly, SEBS underestimates H when there is a shortage of soil water and a lack of vegetation coverage.In contrast, the estimates of H from SEBS-NDWI are more similar to the EC observations, which range from 150 W/m 2 to 400 W/m 2 (Figure 6e-g).The difference between the estimates of H on the four days is small, which indicates that no clear temporal variation occurs in the study area from June to September.The spatial distribution of the instantaneous latent heat flux also differed between the oasis and desert areas (Figure 7).The estimates of E λ from SEBS range from 150 W/m 2 to 300 W/m 2 in the sandy desert, desert steppe, and Gobi desert, which lack soil water and vegetation cover (Figure 7a-d), while the EC-observed λE is less than 200 W/m 2 (Figure 5d).The E λ values modeled by SEBS are greater than the observed λE values from EC under water-limited conditions.The estimates of λE from SEBS-NDWI are consistent with the EC-observed λE in the sandy desert, desert steppe, and Gobi desert, ranging from 0 W/m 2 to 200 W/m 2 (Figure 7e-h).λE exhibits a clear temporal The spatial distribution of the instantaneous latent heat flux also differed between the oasis and desert areas (Figure 7).The estimates of λE from SEBS range from 150 W/m 2 to 300 W/m 2 in the sandy desert, desert steppe, and Gobi desert, which lack soil water and vegetation cover (Figure 7a-d), while the EC-observed λE is less than 200 W/m 2 (Figure 5d).The λE values modeled by SEBS are greater than the observed λE values from EC under water-limited conditions.The estimates of λE from SEBS-NDWI are consistent with the EC-observed λE in the sandy desert, desert steppe, and Gobi desert, ranging from 0 W/m 2 to 200 W/m 2 (Figure 7e-h).λE exhibits a clear temporal variation in the oasis from June to September due to the changes in the vegetation coverage and radiation.With increasing vegetation coverage from June to August, more surface radiation energy is used to exchange water via transpiration.The estimation of λE in the oasis on 3 September is much lower than that on the other three days because of the reduction in the net radiation.
The spatial distribution of the daily evapotranspiration is identical to that of the latent heat flux.The daily ET in the oasis is considerably greater than that in the surrounding area (Figure 8).The daily ET modeled by SEBS ranges from 1.5 mm/d to 4.0 mm/d in the sandy desert, desert steppe, and Gobi desert.In contrast, when the NDWI is integrated into SEBS as an indicator of water stress, the model-produced daily ET values are less than 3 mm/d in the sandy desert, desert steppe, and Gobi desert and are highly spatially variable from the oasis to the desert areas (Figure 8).Furthermore, the daily ET values that were modeled by SEBS-NDWI are greatly reduced in the arid and sparsely vegetated areas, while the high ET values for the oasis remain in the same range as those given by SEBS (Figure 8).A clear temporal variation in the daily ET occurs in the oasis from June to September due to the differences in vegetation coverage and solar radiance.The maximum values of the daily ET in the oasis are 5.88 mm/d, 6.17 mm/d, 6.41 mm/d, and 4.53 mm/d on 24 June, 10 July, 2 August, and 3 September, respectively (Figure 8e-h).In contrast, the temporal variation in the daily ET in the sandy desert, desert steppe, and Gobi desert is not apparent (Figure 8e-h).
Remote Sens. 2015, 7 page-page 14 variation in the oasis from June to September due to the changes in the vegetation coverage and radiation.With increasing vegetation coverage from June to August, more surface radiation energy is used to exchange water via transpiration.The estimation of λE in the oasis on 3 September is much lower than that on the other three days because of the reduction in the net radiation.The spatial distribution of the daily evapotranspiration is identical to that of the latent heat flux.The daily ET in the oasis is considerably greater than that in the surrounding area (Figure 8).The daily ET modeled by SEBS ranges from 1.5 mm/d to 4.0 mm/d in the sandy desert, desert steppe, and Gobi desert.In contrast, when the NDWI is integrated into SEBS as an indicator of water stress, the model-produced daily ET values are less than 3 mm/d in the sandy desert, desert steppe, and Gobi desert and are highly spatially variable from the oasis to the desert areas (Figure 8).Furthermore, the daily ET values that were modeled by SEBS-NDWI are greatly reduced in the arid and sparsely vegetated areas, while the high ET values for the oasis remain in the same range as those given by SEBS (Figure 8).A clear temporal variation in the daily ET occurs in the oasis from June to September due to the differences in vegetation coverage and solar radiance.The maximum values of the daily ET in the oasis are 5.88 mm/d, 6.17 mm/d, 6.41 mm/d, and 4.53 mm/d on 24 June, 10 July, 2 August, and 3 September, respectively (Figure 8e-h).In contrast, the temporal variation in the daily ET in the sandy desert, desert steppe, and Gobi desert is not apparent (Figure 8e-h).The spatial distribution of the daily evapotranspiration is identical to that of the latent heat flux.The daily ET in the oasis is considerably greater than that in the surrounding area (Figure 8).The daily ET modeled by SEBS ranges from 1.5 mm/d to 4.0 mm/d in the sandy desert, desert steppe, and Gobi desert.In contrast, when the NDWI is integrated into SEBS as an indicator of water stress, the model-produced daily ET values are less than 3 mm/d in the sandy desert, desert steppe, and Gobi desert and are highly spatially variable from the oasis to the desert areas (Figure 8).Furthermore, the daily ET values that were modeled by SEBS-NDWI are greatly reduced in the arid and sparsely vegetated areas, while the high ET values for the oasis remain in the same range as those given by SEBS (Figure 8).A clear temporal variation in the daily ET occurs in the oasis from June to September due to the differences in vegetation coverage and solar radiance.The maximum values of the daily ET in the oasis are 5.88 mm/d, 6.17 mm/d, 6.41 mm/d, and 4.53 mm/d on 24 June, 10 July, 2 August, and 3 September, respectively (Figure 8e-h).In contrast, the temporal variation in the daily ET in the sandy desert, desert steppe, and Gobi desert is not apparent (Figure 8e-h).To validate the estimated daily ET from the models, the modeled values of the daily ET from SEBS and SEBS-NDWI were compared with the observed daily ET (Figure 9).The original SEBS overestimated the observed low daily ET (Figure 9a) in villages, sandy deserts, desert steppes, and the Gobi Desert.The proposed model, which considers the influence of water stress on the evaporative fraction, improves the overestimation of the daily ET in these areas where the observed daily ET is less than 3 mm/d (Figure 9b).The RMSE of the estimated daily ET decreases from 1.37 mm/d to 1.09 m/d.The high daily ET values modeled by both SEBS and SEBS-NDWI are less than the observations (Figure 9); this result is likely due to underestimation of net radiation caused by mismatch of the spatial scale between observations and MODIS pixels (Figure 10).The observation towers in the oasis are located in center of farmland and the meteorological variables and turbulent fluxes are observed over homogenous crop or vegetable fields.However, the spatial resolution of MODIS products is 1 km, which is usually mixed by vegetation and bare soil in this region.The land surface temperature and albedo are usually higher for bare soil than those for vegetation covered soil, which will cause that surface temperature and albedo derived from MODIS are higher than observations from towers in vegetation covered farmland.Therefore, the estimated net radiation will decrease due to the negative correlation between net radiation and land surface temperature and albedo according to Equation (4).Additionally, the self-preservation assumption of the EF is not always valid under wet conditions.During the day time, the EF of farmlands curve is concave upward shows that the EF is higher in the early morning and later afternoon [65].Using the EF around noon to estimate daily ET may be another reason for the underestimation of daily ET.Furthermore, uncertainties of parameterization schemes related to vegetation also influence the estimation of daily ET.

15
To validate the estimated daily ET from the models, the modeled values of the daily ET from SEBS and SEBS-NDWI were compared with the observed daily ET (Figure 9).The original SEBS overestimated the observed low daily ET (Figure 9a) in villages, sandy deserts, desert steppes, and the Gobi Desert.The proposed model, which considers the influence of water stress on the evaporative fraction, improves the overestimation of the daily ET in these areas where the observed daily ET is less than 3 mm/d (Figure 9b).The RMSE of the estimated daily ET decreases from 1.37 mm/d to 1.09 m/d.The high daily ET values modeled by both SEBS and SEBS-NDWI are less than the observations (Figure 9); this result is likely due to underestimation of net radiation caused by mismatch of the spatial scale between observations and MODIS pixels (Figure 10).The observation towers in the oasis are located in center of farmland and the meteorological variables and turbulent fluxes are observed over homogenous crop or vegetable fields.However, the spatial resolution of MODIS products is 1 km, which is usually mixed by vegetation and bare soil in this region.The land surface temperature and albedo are usually higher for bare soil than those for vegetation covered soil, which will cause that surface temperature and albedo derived from MODIS are higher than observations from towers in vegetation covered farmland.Therefore, the estimated net radiation will decrease due to the negative correlation between net radiation and land surface temperature and albedo according to equation 4. Additionally, the self-preservation assumption of the EF is not always valid under wet conditions.During the day time, the EF of farmlands curve is concave upward shows that the EF is higher in the early morning and later afternoon [65].Using the EF around noon to estimate daily ET may be another reason for the underestimation of daily ET.Furthermore, uncertainties of parameterization schemes related to vegetation also influence the estimation of daily ET.

15
To validate the estimated daily ET from the models, the modeled values of the daily ET from SEBS and SEBS-NDWI were compared with the observed daily ET (Figure 9).The original SEBS overestimated the observed low daily ET (Figure 9a) in villages, sandy deserts, desert steppes, and the Gobi Desert.The proposed model, which considers the influence of water stress on the evaporative fraction, improves the overestimation of the daily ET in these areas where the observed daily ET is less than 3 mm/d (Figure 9b).The RMSE of the estimated daily ET decreases from 1.37 mm/d to 1.09 m/d.The high daily ET values modeled by both SEBS and SEBS-NDWI are less than the observations (Figure 9); this result is likely due to underestimation of net radiation caused by mismatch of the spatial scale between observations and MODIS pixels (Figure 10).The observation towers in the oasis are located in center of farmland and the meteorological variables and turbulent fluxes are observed over homogenous crop or vegetable fields.However, the spatial resolution of MODIS products is 1 km, which is usually mixed by vegetation and bare soil in this region.The land surface temperature and albedo are usually higher for bare soil than those for vegetation covered soil, which will cause that surface temperature and albedo derived from MODIS are higher than observations from towers in vegetation covered farmland.Therefore, the estimated net radiation will decrease due to the negative correlation between net radiation and land surface temperature and albedo according to equation 4. Additionally, the self-preservation assumption of the EF is not always valid under wet conditions.During the day time, the EF of farmlands curve is concave upward shows that the EF is higher in the early morning and later afternoon [65].Using the EF around noon to estimate daily ET may be another reason for the underestimation of daily ET.Furthermore, uncertainties of parameterization schemes related to vegetation also influence the estimation of daily ET.

Conclusions
ET is mainly controlled by the available radiant energy and soil water.A positive correlation exists between the near-surface soil moisture and the evaporative fraction in regions with different land cover and climates, including the study area.However, SEBS, which estimates evapotranspiration from the energy balance, does not explicitly consider soil moisture when calculating land surface fluxes and ET, but it assumes that the soil moisture information is incorporated into the land surface temperature state variable.This assumption works well when the available radiant energy is the limiting factor for ET, but the ET is overestimated under water-limited conditions in semi-arid and arid regions.
In this study, we proposed an explicit integration of the NDWI as a water stress index into SEBS through a modified definition of kB ´1.Significant positive correlations exist between the NDWI and EF, with correlation coefficients of 0.62 and 0.75, respectively, for the two sets of experimental data.Integrating the NDWI as an indicator of water stress into SEBS causes kB ´1 to decrease as the water stress increases and causes the latent heat flux and ET to decrease as the sensible heat flux increases.The modeled land surface fluxes improve when the NDWI is incorporated into SEBS, particularly for low values of the NDWI, where soil water and vegetation coverage are lacking.The proposed scheme not only significantly improves the estimation of surface turbulent fluxes under water-limited conditions but also provides good results when adequate soil water is available.The regional daily ET can also be reasonably estimated using only MODIS data products and meteorological data.The SEBS-NDWI method not only accurately captures the spatial distribution of the daily ET from the oasis to the desert areas but also represents the temporal variation in the daily ET in the oasis.However, the modeled values of the daily ET from SEBS and SEBS-NDWI are lower than the observations for high daily ET values because of the underestimation of the daily net radiation and the invalid self-preservation assumption of the EF.
The developed SEBS-NDWI scheme provides an effective method of estimating the regional ET in semi-arid and arid regions, and it only requires MODIS data products and meteorological data as inputs.The scheme effectively resolves the requirements of high spatial resolution soil moisture products in SEBS under water stress conditions [31].However, the estimation of ET in dense vegetation areas should be improved, particularly because of the uncertainties in the parameterization schemes of the land surface parameters and daily net radiation, which are usually related to the study region and the selected sensors.Additionally, the developed SEBS-NDWI does not take effect in dense vegetation areas, because canopy density also influences the correlation between NDWI and EF.In other words, NDWI and other water stress indices present a strong correlation with structural indices (Normalized Difference Vegetation Index, such as NDVI) related to canopy density [66].This correlation should be removed to capture the vegetation stress and it will be considered in future studies.

Figure 1 .
Figure 1.Spatial locations of the field observation stations in study area and land use map in the kernel experimental area.

Figure 1 .
Figure 1.Spatial locations of the field observation stations in study area and land use map in the kernel experimental area.

Figure 2 .
Figure 2. Box and whisker plots of near-surface soil moisture (SM) and evaporative fraction (EF) at (a) the Yingke station during the WATER experiment (2008-2011) and (b) sixteen stations of the HiWATER experiment (2012).SM denotes the daily average values of soil moisture measurements at a 5 cm depth, and EF is the evaporative fraction of the clear-sky daytime averages of the EC measurements.

Figure 2 .
Figure 2. Box and whisker plots of near-surface soil moisture (SM) and evaporative fraction (EF) at (a) the Yingke station during the WATER experiment (2008-2011) and (b) sixteen stations of the HiWATER experiment (2012).SM denotes the daily average values of soil moisture measurements at a 5 cm depth, and EF is the evaporative fraction of the clear-sky daytime averages of the EC measurements.

Figure 3 .
Figure 3. Relationships between the normalized difference water index (NDWI) and evaporative fraction (EF) at (a) the Yingke station during the WATER experiment (2008-2011) and (b) sixteen HiWATER stations (2012).EF is the evaporative fraction of the clear-sky daytime averages of the EC measurements.

Figure 3 .
Figure 3. Relationships between the normalized difference water index (NDWI) and evaporative fraction (EF) at (a) the Yingke station during the WATER experiment (2008-2011) and (b) sixteen HiWATER stations (2012).EF is the evaporative fraction of the clear-sky daytime averages of the EC measurements.

H
) for high values of o H compared with the estimates of H from the original SEBS ( SEBS H

Figure 4 .
Figure 4. Comparison of the modeled and observed surface turbulent fluxes at the Yingke station in 2008-2011.(a) and (b) show scatter plots of the measured and observed sensible heat fluxes from the original SEBS and the NDWI-integrated SEBS, respectively; (c) and (d) show scatter plots of the measured and observed latent heat fluxes from the original SEBS and the NDWI-integrated SEBS, respectively.

Figure 4 .
Figure 4. Comparison of the modeled and observed surface turbulent fluxes at the Yingke station in 2008-2011.(a) and (b) show scatter plots of the measured and observed sensible heat fluxes from the original SEBS and the NDWI-integrated SEBS, respectively; (c) and (d) show scatter plots of the measured and observed latent heat fluxes from the original SEBS and the NDWI-integrated SEBS, respectively.
compares the modeled and observed surface turbulent fluxes at the 16 stations, and Table 3 shows the error evaluation of the surface turbulent fluxes from the original SEBS and NDWI-integrated SEBS based on EC observations.The modeled values of H from the original SEBS (H SEBS ) are clearly underestimated when the observed values of H (H o ) are greater than 100 W/m 2 (Figure 5a), and the modeled values of H from the NDWI-integrated SEBS (H SEBS´NDW I ) are closer to H o for high values of H o (Figure 5c).When the NDWI is integrated into SEBS, the estimated values of λE decrease with as the estimated H increases for low NDWI values.The estimated values of λE remain nearly unchanged for high NDWI values.The modeled λE clearly improves when the NDWI is included in SEBS; the RMSE and BIAS decrease by approximately 33 W/m 2 (Table

2 (H
underestimated when the observed values of H ( o H ) are greater than 100 W/m Figure 5a), and the modeled values of H from the NDWI-integrated SEBS ( for high values of o H ( Figure 5c).When the NDWI is integrated into SEBS, the estimated values of E λ decrease with as the estimated H increases for low NDWI values.The estimated values of E λ remain nearly unchanged for high NDWI values.The modeled E λ clearly improves when the NDWI is included in SEBS; the RMSE and BIAS decrease by approximately 33 W/m 2 (Table

Figure 5 .
Figure 5.Comparison of modeled and observed surface turbulent fluxes at 16 stations in 2012.(a) and (b) show scatter plots of the measured and observed sensible heat fluxes from the original SEBS and the NDWI-integrated SEBS, respectively; (c) and (d) show scatter plots of the measured and observed latent heat fluxes from the original SEBS and the NDWI-integrated SEBS, respectively.

Figure 5 .
Figure 5.Comparison of modeled and observed surface turbulent fluxes at 16 stations in 2012.(a) and (b) show scatter plots of the measured and observed sensible heat fluxes from the original SEBS and the NDWI-integrated SEBS, respectively; (c) and (d) show scatter plots of the measured and observed latent heat fluxes from the original SEBS and the NDWI-integrated SEBS, respectively.

Figure 6 .
Figure 6.Instantaneous sensible heat flux maps of the study area in 2012.Figure 6 shows the instantaneous sensible heat flux maps obtained with the original SEBS for 24 June (a), 10 July (b), 2 August (c), and 3 September (d), respectively; and the instantaneous sensible heat flux maps obtained with the NDWI-integrated SEBS for 24 June (e), 10 July (f), 2 August (g), and 3 September (h), respectively.

Figure 6 .
Figure 6.Instantaneous sensible heat flux maps of the study area in 2012.Figure 6 shows the instantaneous sensible heat flux maps obtained with the original SEBS for 24 June (a), 10 July (b), 2 August (c), and 3 September (d), respectively; and the instantaneous sensible heat flux maps obtained with the NDWI-integrated SEBS for 24 June (e), 10 July (f), 2 August (g), and 3 September (h), respectively.

Figure 7 .
Figure 7. Instantaneous latent heat flux maps of the study area in 2012.Figure 7 show the instantaneous latent heat flux maps obtained with the original SEBS for 24 June (a), 10 July (b), 2 August (c), and 3 September (d), respectively; and the the instantaneous latent heat flux maps obtained with the NDWI-integrated SEBS for 24 June (e), 10 July (f), 2 August (g), and 3 September (h), respectively.

Figure 7
Figure 7. Instantaneous latent heat flux maps of the study area in 2012.Figure 7 show the instantaneous latent heat flux maps obtained with the original SEBS for 24 June (a), 10 July (b), 2 August (c), and 3 September (d), respectively; and the the instantaneous latent heat flux maps obtained with the NDWI-integrated SEBS for 24 June (e), 10 July (f), 2 August (g), and 3 September (h), respectively.

Figure 8 .
Figure 8. Daily evapotranspiration maps of the study area in 2012.Figure 8 show the daily evapotranspiration maps obtained with the original SEBS for 24 June (a), 10 July (b), 2 August(c), and 3 September (d), respectively; and the daily evapotranspiration maps obtained with the NDWIintegrated SEBS for 24 June (e), 10 July (f), 2 August (g), and 3 September (h), respectively.

Figure 7 .
Figure 7. Instantaneous latent heat flux maps of the study area in 2012.Figure 7 show the instantaneous latent heat flux maps obtained with the original SEBS for 24 June (a), 10 July (b), 2 August (c), and 3 September (d), respectively; and the the instantaneous latent heat flux maps obtained with the NDWI-integrated SEBS for 24 June (e), 10 July (f), 2 August (g), and 3 September (h), respectively.

Figure 7
Figure 7. Instantaneous latent heat flux maps of the study area in 2012.Figure 7 show the instantaneous latent heat flux maps obtained with the original SEBS for 24 June (a), 10 July (b), 2 August (c), and 3 September (d), respectively; and the the instantaneous latent heat flux maps obtained with the NDWI-integrated SEBS for 24 June (e), 10 July (f), 2 August (g), and 3 September (h), respectively.

Figure 7 .
Figure 7. Instantaneous latent heat flux maps of the study area in 2012.Figure 7 show the instantaneous latent heat flux maps obtained with the original SEBS for 24 June (a), 10 July (b), 2 August (c), and 3 September (d), respectively; and the the instantaneous latent heat flux maps obtained with the NDWI-integrated SEBS for 24 June (e), 10 July (f), 2 August (g), and 3 September (h), respectively.

Figure 7
Figure 7. Instantaneous latent heat flux maps of the study area in 2012.Figure 7 show the instantaneous latent heat flux maps obtained with the original SEBS for 24 June (a), 10 July (b), 2 August (c), and 3 September (d), respectively; and the the instantaneous latent heat flux maps obtained with the NDWI-integrated SEBS for 24 June (e), 10 July (f), 2 August (g), and 3 September (h), respectively.

Figure 8 .
Figure 8. Daily evapotranspiration maps of the study area in 2012.Figure 8 show the daily evapotranspiration maps obtained with the original SEBS for 24 June (a), 10 July (b), 2 August(c), and 3 September (d), respectively; and the daily evapotranspiration maps obtained with the NDWIintegrated SEBS for 24 June (e), 10 July (f), 2 August (g), and 3 September (h), respectively.

Figure 8 .
Figure 8. Daily evapotranspiration maps of the study area in 2012.Figure 8 show the daily evapotranspiration maps obtained with the original SEBS for 24 June (a), 10 July (b), 2 August(c), and 3 September (d), respectively; and the daily evapotranspiration maps obtained with the NDWI-integrated SEBS for 24 June (e), 10 July (f), 2 August (g), and 3 September (h), respectively.

Figure 9 .
Figure 9.Comparison of modeled and observed daily evapotranspiration at 16 stations on 24 June, 10 July, 2 August, and 3 September, 2012.(a) and (b) are scatter plots of the daily observations of evapotranspiration and the results from the original SEBS and NDWI-integrated SEBS.

Figure 10 .
Figure 10.Comparison of modeled daily net radiation from SEBS and daily net radiation observed from the 16 AMSs on 24 June, 10 July, 2 August, and 3 September, 2012.

Figure 9 .
Figure 9.Comparison of modeled and observed daily evapotranspiration at 16 stations on 24 June, 10 July, 2 August, and 3 September, 2012.(a) and (b) are scatter plots of the daily observations of evapotranspiration and the results from the original SEBS and NDWI-integrated SEBS.

Figure 9 .
Figure 9.Comparison of modeled and observed daily evapotranspiration at 16 stations on 24 June, 10 July, 2 August, and 3 September, 2012.(a) and (b) are scatter plots of the daily observations of evapotranspiration and the results from the original SEBS and NDWI-integrated SEBS.

Figure 10 .
Figure 10.Comparison of modeled daily net radiation from SEBS and daily net radiation observed from the 16 AMSs on 24 June, 10 July, 2 August, and 3 September, 2012.Figure 10.Comparison of modeled daily net radiation from SEBS and daily net radiation observed from the 16 AMSs on 24 June, 10 July, 2 August, and 3 September, 2012.

Figure 10 .
Figure 10.Comparison of modeled daily net radiation from SEBS and daily net radiation observed from the 16 AMSs on 24 June, 10 July, 2 August, and 3 September, 2012.Figure 10.Comparison of modeled daily net radiation from SEBS and daily net radiation observed from the 16 AMSs on 24 June, 10 July, 2 August, and 3 September, 2012.

Table 1 .
Description of observation stations in the WATER and HiWATER experiments.

Table 2 .
Error evaluations of the estimated sensible heat flux and latent heat flux from the original SEBS and the NDWI-integrated SEBS at the Yingke station (2008-2011).

Table 3 .
Error evaluations of the estimated sensible heat flux and latent heat flux from the original SEBS and the NDWI-integrated SEBS at 16 stations in 2012.