Regional Estimation of Remotely Sensed Evapotranspiration Using the Surface Energy Balance-advection (seb-a) Method

Evapotranspiration (ET) is an essential part of the hydrological cycle and accurately estimating it plays a crucial role in water resource management. Surface energy balance (SEB) models are widely used to estimate regional ET with remote sensing. The presence of horizontal advection, however, perturbs the surface energy balance system and contributes to the uncertainty of energy influxes. Thus, it is vital to consider horizontal advection when applying SEB models to estimate ET. This study proposes an innovative and simplified approach, the surface energy balance-advection (SEB-A) method, which is based on the energy balance theory and also takes into account the horizontal advection to determine ET by remote sensing. The SEB-A method considers that the actual ET consists of two parts: the local ET that is regulated by the energy balance system and the exotic ET that arises from horizontal advection. To evaluate the SEB-A method, it was applied to the middle region of the Heihe River in China. Instantaneous ET for three days were acquired and assessed with ET measurements from eddy covariance (EC) systems. The results demonstrated that the ET estimates had a high accuracy, with a correlation coefficient (R 2) of 0.713, a mean average error (MAE) of 39.3 W/m 2 and a root mean square error (RMSE) of 54.6 W/m 2 between the estimates and corresponding measurements. Percent error was calculated to more rigorously assess the accuracy of these estimates, and it ranged from 0% to 35%, with over 80% of the locations within a 20% error. To better understand the SEB-A method, the relationship between the ET estimates and land use types was analyzed, and the results indicated that the ET estimates had spatial distributions that correlated with vegetation patterns and could well demonstrate the ET differences caused by different land use types. The sensitivity analysis suggested that the SEB-A method requested accurate estimation of the available energy, R n ´ G, but was less constrained with the difference between ground and air temperature, T 0 ´ T a–loc .


Introduction
Evapotranspiration (ET), which is a combined process of evaporation and transpiration, is crucial to the hydrological cycle.It links energy and water exchanges among the hydrosphere, atmosphere and biosphere [1][2][3][4][5].ET is the largest outgoing water flux from the earth's surface; most precipitation is lost in the form of ET with the percentage varying from region to region globally.In semi-arid regions, the amount of ET is roughly equivalent to the precipitation [6].Therefore, the accurate quantification of ET is critical to developing a greater understanding of local and global energy and water cycles, global climate change, land atmosphere interactions and ecosystem processes, and beneficial to improving applications in many fields, such as water resources management, drought monitoring and assessment [7], irrigation scheduling, hydrological modeling, weather forecasts and so on [8][9][10][11][12][13].Ground-based observations of ET have some advantages but the required ground instrumentation is available in a limited number of places and can only obtain localized values of ET, which often cannot satisfy the needs of regional-scale studies.In addition, the employment of instrumentations is always time consuming, labor intensive and sometimes subject to instrument failure [9].
Recent developments in remote sensing have made it possible to acquire crucial variables for characterizing land surface interactions.Over the last two decades, the requirement of quantifying ET at a regional scale, together with the recent advances in satellite remote sensing technology, has led to many studies on mapping large-scale ET [9,11,14,15].The surface energy balance algorithm based on the energy balance equation (SEB; Equation ( 1)), is one of the most widely used approaches for obtaining regional estimates of remotely sensed ET at multiple temporal and spatial scales [9].Energy and mass are exchanged between the land surface and the atmosphere [16,17].When the heat fluxes that are transported by horizontal advection and consumed by photosynthetic vegetation are not considered, the one-dimensional form of the surface energy balance equation on an instantaneous time scale can be expressed as: R n " G `H `LE where R n (W/m 2 ) is the net surface radiation, G (W/m 2 ) is the soil heat flux, H (W/m 2 ) is the sensible heat flux, and LE (W/m 2 ) is the latent heat flux.In Equation (1), sensible heat flux can be calculated using the remotely sensed surface temperature, air temperature and a series of conductance formulas; estimates of R n and G are available using remote sensing or ground data, and hence, the latent heat flux is taken as a residual term of Equation (1) [16,17].
Based on how the sinks or sources of heat fluxes are parameterized between the earth's surface and the atmosphere, SEB models can be divided into two categories: the single-source models and two-source models, or multi-source models [18].Single-source SEB models consider vegetation and soil to be a single unit and often utilize land surface temperature (LST) and empirical resistance corrections to derive the latent heat flux [19].Since they are highly dependent on empirical corrections, single-source SEB models are inappropriate for application over areas with partial vegetation cover, where significant errors would be found [20][21][22].To overcome the limitations associated with single-source SEB models in partially vegetated areas [22], it is necessary to apply a two-source SEB model that separately simulates the heat and water exchange and the interactions between soil and atmosphere and between vegetation and atmosphere.Thus, in two-source SEB models, the component H is partitioned between the soil and vegetation, and the component LE is treated as the sum of evaporation from the soil surface and transpiration from vegetation.
From the aforementioned review of SEB models, we know that both single-source and two-source models are based on the energy balance equation (Equation ( 1)) that is founded on the premise of ignoring the energy transported by horizontal advection and consumed by photosynthesis [16,17].Nevertheless, the presence of horizontal advection always characterizes exchanges in the near surface layer and contributes to the uncertainty of energy influxes [23,24].Neglecting systematic sampling errors, systematic instrument bias, low and high frequency loss of turbulent fluxes and energy sinks, horizontal advection of heat and water vapor is regarded as the primary reason for the non-closure of the surface energy balance [25][26][27].To a certain region, uncertainties from horizontal advection on the surface energy balance have two aspects.One is the horizontal advection itself, which develops on heterogeneous surfaces and exerts profound influence on the energy balance equation.This has been reported in many previous studies [28,29].Therefore, horizontal advection should be discussed as one reason for energy imbalances and regarded as a necessary item to add to the energy balance equation.The other uncertainty refers to the ground-based meteorological variables of air temperature, humidity and other auxiliary surface measurements inputted to the energy balance equation.These meteorological variables are actually measured and thus characterize the practical circumstance rather than describe the environment under the condition of energy balance closure, which means that the actually observed meteorological variables have already contained effects from the horizontal advection and would affect the energy balance equation when they are used as inputs.In conclusion, horizontal advection affects the energy balance equation directly through energy inputs or outputs, whereas ground-based observations impact the energy balance equation indirectly by including the effects of advection.
In view of the prevalence of the horizontal advection and its influence on the energy balance system, it is of great necessity to take into account the horizontal advection when applying SEB models to estimate ET.The objective of this paper is to propose an innovative method, the surface energy balance-advection (SEB-A) method, which is based on the energy balance equation and considers the horizontal advection as well to derive ET.The SEB-A method assumes that the energy consumed in the actual ET consists of two parts; one part is provided by the energy balance system and the other part is brought in by the horizontal advection.The former, called local ET, can be obtained by the modified energy balance equation (Equation ( 3)), and the calculation method of it is similar to previous studies that are on the SEB models; the latter, called exotic ET, can be calculated with the help of ground-based EC observations.Section 2 describes the datasets used to assess the method.Section 3.1 presents a description on the theory of calculating the local ET that is regulated by the energy balance system, and Section 3.2 introduces approaches for deriving the exotic ET that is driven by the horizontal advection.Section 4 shows the results of the SEB-A method that was applied on a regional scale by using variables retrieved from ASTER imagery and HJ-1A/B imagery and measured at meteorological and EC stations.Section 5 discusses the advantages and limitations of the SEB-A method.Section 6 concludes the work.

Study Area
The study area is located in Zhangye, Gansu, China, and it is from 38.83 ˝N to 38.93 ˝N in latitude and from 100.32 ˝E to 100.42 ˝E in longitude, as shown in Figure 1.It belongs to the middle reaches of Heihe River that is the second longest inland river in China.The annual mean air temperature is 7.0 ˝C.The climate is dry and rainless, with an annual mean precipitation of 124.9 mm and a potential evaporation of 2000 mm.The ground is flat and cultivated with maize, vegetable crops, orchards and woodlands.Some built-up areas (villages) are spread throughout the study area.Every year from May to September, staple crops are grown, including seed corns, vegetables and fruit trees.
To estimate the ET, remote sensing data that include ASTER images and HJ-1A/B images and ground-based observations that include meteorological variables and surface fluxes are required.

Remote Sensing Data
The ASTER is a sensor that collects multispectral visible, near-infrared and thermal infrared images.It is a special remote sensing tool intended to monitor land surface energy balance, hydrological processes and climate [30].ASTER images with a high spatial resolution of 15 m in both the visible and near-infrared wave bands (0.52-0.86 µm) were used in this study to estimate the vegetation fraction and the surface emissivity.The LST products, which had a spatial resolution of 90 m, were acquired from the thermal infrared wave band (8.125-11.65 µm) of ASTER images and provided by the Heihe Plan Science Data Center [31,32].The albedo in the study area was calculated from images acquired from the HJ-1A/B, which is a small satellite constellation that is used for monitoring and forecasting environmental changes and disasters.The HJ-1A/B was launched by China in 2008.Two CCDs (0.45-0.89 µm) with a spatial resolution of 30 m are carried by each satellite.

Ground Observations
In the study area, there are 17 stations (Figure 1) equipped with Automatic Meteorological Station Systems and EC Systems that are both automatic data acquisition systems.The Automatic Meteorological Station System can acquire meteorological variables, including air temperature, humidity, air pressure, solar radiation, wind speed, wind direction and so on, at sample intervals of 10 min and one minute.The EC System can monitor surface fluxes such as sensible heat flux, latent heat flux and so on, at a sample interval of 30 min.
All data mentioned above were provided by the HiWATER experiment which is an watershed-scale and eco-hydrological experiment designed from an interdisciplinary perspective to address problems that include heterogeneity, scaling, uncertainty and closing of the water cycle at the watershed scale [33].The experiment was performed in the Heihe River Basin.Users can attain the data by submitting an application to the data center, Cold and Arid Regions Science Data, at Lanzhou (http://westdc.westgis.ac.cn/).
provided by the Heihe Plan Science Data Center [31,32].The albedo in the study area was calculated from images acquired from the HJ-1A/B, which is a small satellite constellation that is used for monitoring and forecasting environmental changes and disasters.The HJ-1A/B was launched by China in 2008.Two CCDs (0.45-0.89 µm) with a spatial resolution of 30 m are carried by each satellite.

Ground Observations
In the study area, there are 17 stations (Figure 1) equipped with Automatic Meteorological Station Systems and EC Systems that are both automatic data acquisition systems.The Automatic Meteorological Station System can acquire meteorological variables, including air temperature, humidity, air pressure, solar radiation, wind speed, wind direction and so on, at sample intervals of 10 min and one minute.The EC System can monitor surface fluxes such as sensible heat flux, latent heat flux and so on, at a sample interval of 30 min.
All data mentioned above were provided by the HiWATER experiment which is an watershed-scale and eco-hydrological experiment designed from an interdisciplinary perspective to address problems that include heterogeneity, scaling, uncertainty and closing of the water cycle at the watershed scale [33].The experiment was performed in the Heihe River Basin.Users can attain the data by submitting an application to the data center, Cold and Arid Regions Science Data, at Lanzhou (http://westdc.westgis.ac.cn/).

Methodology
In this section, the SEB-A method is described in detail.According to previous studies, horizontal advection is the main reason for surface energy imbalances [25][26][27].As a result, without regard to the internal generating mechanism of ET, the actual ET ( ) is directly composed of two parts: the local ET ( ) dominated by the energy balance system, and the exotic ET ( ) given rise to by horizontal advection.This is inspired by the research of Zhang et al. [34].
With no advection or when there is a surface energy balance closure, the actual ET ( ) is equal to the local ET ( ) and can be calculated by the energy balance equation.In reality, however, the horizontal advection would perturb the energy balance system and influence the energy balance equation.In this case, the actual ET ( ) can be regarded as the sum of the local ET ( ) and the exotic ET ( ).The basic idea can be described as: the local driving force (the energy balance system) produces a local ET, and the exotic driving force (horizontal advection) changes the local ET by adding an exotic ET.The actual ET can be expressed as:

Methodology
In this section, the SEB-A method is described in detail.According to previous studies, horizontal advection is the main reason for surface energy imbalances [25][26][27].As a result, without regard to the internal generating mechanism of ET, the actual ET (LE act ) is directly composed of two parts: the local ET (LE loc ) dominated by the energy balance system, and the exotic ET (LE exo ) given rise to by horizontal advection.This is inspired by the research of Zhang et al. [34].
With no advection or when there is a surface energy balance closure, the actual ET (LE act ) is equal to the local ET (LE loc ) and can be calculated by the energy balance equation.In reality, however, the horizontal advection would perturb the energy balance system and influence the energy balance equation.In this case, the actual ET (LE act ) can be regarded as the sum of the local ET (LE loc ) and the exotic ET (LE exo ).The basic idea can be described as: the local driving force (the energy balance system) produces a local ET, and the exotic driving force (horizontal advection) changes the local ET by adding an exotic ET.The actual ET can be expressed as: where LE act is the actual ET, LE loc is the local ET, and LE exo is the exotic ET.Three items are measured in energy units, W/m 2 .Figure 2 shows the principle of the SEB-A method.
where is the actual ET, is the local ET, and is the exotic ET.Three items are measured in energy units, W/m 2 .
Figure 2 shows the principle of the SEB-A method.Because contains both and , to estimate , it is necessary to calculate the two parts: and .

Figure 2.
The principle of the surface energy balance-advection (SEB-A) method.

Obtaining the Local ET,
The can be determined by using the energy balance equation.Under the condition of energy balance closure, that is to say, there is no horizontal advection, the energy exchange between the land surface and atmosphere can be described by the modified energy balance equation (Equation ( 3)).Here, we also ignore the energy consumed by photosynthesis.
where (W/m 2 ) and (W/m 2 ) are the sensible heat flux and the latent heat flux driven by the energy balance system, respectively.
To calculate , we need to know the sensible heat flux, , as well.It is given as Equation (4): where is the air density, (1004 J/(Kg•K)) is the specific heat at constant pressure of air, (s/m) is the aerodynamic resistance, (°C) is the aerodynamic temperature that is usually substituted with LST in applications [35,36], and (°C) is the surface air temperature when there is no advection.Thus, can be derived as: Since there is no horizontal advection, the variables in Equation ( 5) should be acquired under the condition of energy balance closure.Among these variables, , , and LST are thought to remain the same, regardless of whether the energy balance is closed or not on an instantaneous time scale, whereas the other two variables, and -, are specific and need to be calculated in the case of energy balance closure.

Obtaining the Local ET,LE loc
The LE loc can be determined by using the energy balance equation.Under the condition of energy balance closure, that is to say, there is no horizontal advection, the energy exchange between the land surface and atmosphere can be described by the modified energy balance equation (Equation ( 3)).Here, we also ignore the energy consumed by photosynthesis.
where H loc (W/m 2 ) and LE loc (W/m 2 ) are the sensible heat flux and the latent heat flux driven by the energy balance system, respectively.To calculate LE loc , we need to know the sensible heat flux, H loc , as well.It is given as Equation ( 4): where ρ is the air density, C p (1004 J/(Kg¨K)) is the specific heat at constant pressure of air, r a-loc (s/m) is the aerodynamic resistance, T 0 ( ˝C) is the aerodynamic temperature that is usually substituted with LST in applications [35,36], and T a-loc ( ˝C) is the surface air temperature when there is no advection.Thus, LE loc can be derived as: Since there is no horizontal advection, the variables in Equation ( 5) should be acquired under the condition of energy balance closure.Among these variables, R n , G, ρC p and LST are thought to remain the same, regardless of whether the energy balance is closed or not on an instantaneous time scale, whereas the other two variables, r a-loc and T a-loc , are specific and need to be calculated in the case of energy balance closure.

Calculation of r a-loc
The measurement of the special r a-loc was through a laboratory experiment.In the experiment, the observational object was bare soil with a saturated soil water content.Parameters for the calculation including the radiometric soil surface temperature (T s ), actual surface vapor pressure (e a ) and evaporation (LE), were measured with no wind.The formula used for deriving r a-loc is written as: where γ is the psychrometric constant (Pa/K), r s (s/m) is the surface water vapor resistance and equals zero when the soil is fully saturated, e s (Pa) is the surface saturated vapor pressure, and LE (g) is the evaporation from the soil and measured by the weighing method.
In the experiment, we set five different initial air conditions with different T a and RH, and we calculated r a-loc once per hour from 08:00 to 18:00 at the heights of 1 m and 0.5 m, respectively.Tables 1 and 2 show the calculation results of r a-loc at the heights of 0.5 m and 1 m.In the tables, T a is the air temperature, RH is the air relative humidity, and av. is the average value of r a-loc .The r a-loc varied around 150 s/m, and therefore, a value of r a-loc , 150 s/m, was determined and used to compute LE loc .Because T a-loc is specific under the condition of energy balance closure, it can be calculated by the energy balance equation (Equation ( 3)).Together with the sensible heat flux equation (Equation ( 4)) and the Bowen ratio, β, equation (Equation ( 7)), T a-loc can be derived as follows: T a-loc " LST ´β pR n ´Gq β `1 In Equation ( 8), the Bowen ratio, β, is determined according to an empirical equation presented by Zhang et al. [37].The equation is written as: where A is an empirical coefficient that equals 0.66, P (W¨m ´2¨s 0.5 ¨˝C) is the simplified thermal inertia, and P max and P min are, respectively, the maximum and minimum thermal inertia that are determined from the trapezoid space formed by the vegetation index (VI) and LST.The simplified thermal inertia P is defined as: where R n (W/m 2 ) is the average net radiation during t 1 (s) and t 2 (s), T 02 ( ˝C) is the surface temperature when the satellite is passing, and T 01 ( ˝C) is the daily minimum surface temperature.t 2 and t 1 are the occurring times of T 02 and T 01 , respectively.

Calculation of R n and G
Net radiation, R n , can be derived by Equation (11); soil heat flux, G, can be estimated by an empirical equation (Equation (12)) that is used in some previous studies [34].
R n " p1 ´αq S 0 `Rld ´εs σT 0 4 (11) where α is the surface albedo, S 0 (W/m 2 ) is the downward shortwave radiation, R ld (W/m 2 ) is the downward longwave radiation, ε s is the surface emissivity, σ (5.67 ˆ10 ´8 W¨m ´2¨K ´4) is the Stefan-Boltzmann constant, and f v is the fraction of vegetated surface.

Obtaining the Exotic ET, LE exo
Referring to LE exo , it requires ground-based observations of latent heat flux to determine its value.For each EC station inputted to the method, we know its local ET, which can be estimated according to Section 3.1 and its actual ET, which is measured by EC stations; then the exotic ET can be determined by Equation (2).However, the attained exotic ET, LE exo , is ground-based and needs to be extended to a regional scale from a point scale.In the study, we used a spatial interpolation method, the inverse distance weighting (IDW), to interpolate the ground-based LE exo to acquire the spatial distribution results.We assume that LE exo is affected by the intensity of horizontal advection, which moves to other places through diffusion.Since the IDW interpolates according to the distance between two inputted objectives, it is appropriate to use it to acquire the spatial distribution of LE exo .The process of interpolation is described as follows.
If the pixel covering the inputted EC station is expressed as Pixel i with LE exopiq , the exotic ET of other pixels (Pixel x, LE exopxq ) around Pixel i can be calculated by the IDW.Formulas can be written as: where LE exopxq (W/m 2 ) is the exotic ET of any pixel to be calculated, i is the number of the inputted EC stations, w piq is the weight, d i (m) is the distance between the Pixel i and Pixel x, and k is the specified power exponent.

Results
Three clear sky days, 10 July, 11 August and 27 August 2012, were chosen to conduct the study.The spatially derived ET on three days was arrived at the moment of the ASTER satellite overpass, as shown in Figure 3a-c, respectively.To evaluate the results estimated by the SEB-A method, a result validation, an error analysis, a relationship analysis between ET estimates and land use status and a sensitivity analysis were performed.
Remote Sens.2016, 8, 644; doi:10.3390/rs80806448 of 14 where ( ) (W/m 2 ) is the exotic ET of any pixel to be calculated, is the number of the inputted EC stations, ( ) is the weight, (m) is the distance between the Pixel and Pixel x, and is the specified power exponent.

Results
Three clear sky days, 10 July, 11 August and 27 August 2012, were chosen to conduct the study.The spatially derived ET on three days was arrived at the moment of the ASTER satellite overpass, as shown in Figure 3a-c, respectively.To evaluate the results estimated by the SEB-A method, a result validation, an error analysis, a relationship analysis between ET estimates and land use status and a sensitivity analysis were performed.

Validation of ET Estimates
The results were validated against measurements from ground-based EC stations.For each of the three days, nine EC stations were used as inputs to the SEB-A method to calculate the exotic ET, and the other eight EC stations were used to verify the estimates.Because three days' results were tested together, a total of twenty-four points were conducted for the validation.The comparison between estimates obtained by the SEB-A method and measurements from EC stations is shown in Figure 4.It shows that the estimates have an R 2 of 0.713, a MAE of 39.3 W/m 2 and an RMSE of 54.6 W/m 2 .Except for three points that have absolute errors higher than 100 W/m 2 and two points that have absolute errors between 60 W/m 2 and 70 W/m 2 , the other points have absolute errors lower than 20 W/m 2 .Validation results are similar to previous studies.Kalma made an assessment of some 30 published validations that showed an RMSE value of approximately 50 W/m 2 and relative errors of 15%-30% [15].The correlation analysis indicated that the simulated results were quite consistent with the available measurements, illustrating that ET estimates obtained by the SEB-A method were close to the 1:1 line and had high accuracy.

Validation of ET Estimates
The results were validated against measurements from ground-based EC stations.For each of the three days, nine EC stations were used as inputs to the SEB-A method to calculate the exotic ET, and the other eight EC stations were used to verify the estimates.Because three days' results were tested together, a total of twenty-four points were conducted for the validation.The comparison between estimates obtained by the SEB-A method and measurements from EC stations is shown in Figure 4.It shows that the estimates have an R 2 of 0.713, a MAE of 39.3 W/m 2 and an RMSE of 54.6 W/m 2 .Except for three points that have absolute errors higher than 100 W/m 2 and two points that have absolute errors between 60 W/m 2 and 70 W/m 2 , the other points have absolute errors lower than 20 W/m 2 .Validation results are similar to previous studies.Kalma made an assessment of some 30 published validations that showed an RMSE value of approximately 50 W/m 2 and relative errors of 15%-30% [15].The correlation analysis indicated that the simulated results were quite consistent with the available measurements, illustrating that ET estimates obtained by the SEB-A method were close to the 1:1 line and had high accuracy.

Error Analysis
To further evaluate the ET estimates, the frequency distribution of the percent errors was analyzed.According to numerical analysis and error analysis theory, the percent error is more scientifically rigorous than the absolute error when making assessments on the accuracy of estimates [38].The percent error is the absolute error divided by the magnitude of the exact value, and then expressed in terms of per 100.The expression can be written as: where is the percent error, ∆ is the absolute error and L is the exact value.Figure 5 displays the frequency distribution of the percent errors on ET estimates, which varies from 0% to 35%.An accuracy of the percent errors within 20% was achieved at more than 80% of study locations, and over 40% of the locations had percent errors lower than 5%.

Error Analysis
To further evaluate the ET estimates, the frequency distribution of the percent errors was analyzed.According to numerical analysis and error analysis theory, the percent error is more scientifically rigorous than the absolute error when making assessments on the accuracy of estimates [38].The percent error is the absolute error divided by the magnitude of the exact value, and then expressed in terms of per 100.The expression can be written as: where δ is the percent error, ∆ is the absolute error and L is the exact value.Figure 5 displays the frequency distribution of the percent errors on ET estimates, which varies from 0% to 35%.An accuracy of the percent errors within 20% was achieved at more than 80% of study locations, and over 40% of the locations had percent errors lower than 5%.

Relationship between ET Estimates and Land Use Status
Since spatial distributions of ET estimates were obtained (Figure 3), we can conduct deep analyses by combining the land use types (Figure 1).The land use type within the study area, which was based on the aerial remote sensing data of the Compact Airborne Spectrographic Imager (CASI) and used the Support Vector Machine (SVM) for classification, was also provided by the Heihe Plan Science Data Center [39].
There were five land uses types in the study area, including maize growing, village habitats, vegetable crops, orchards and woodlands.For each land use type, we performed statistical computation and gave an average ET.Table 3 shows the average ET of each land use type on the three days.Maize had the maximum ET, and the average ETs on the three days were 528 W/m 2 , 389 W/m 2 and 425 W/m 2 , respectively.Vegetable crops, orchards and woodlands had an average ET between 400 W/m 2 and 500 W/m 2 on 10 July and an average ET of approximately 300 W/m 2 on the other two days.Villages showed the lowest ET, with an average value of 330 W/m 2 on 10 July and less than 200 W/m 2 on the other two days.The spatial distribution of ET was intimately linked with vegetation patterns.In general, high ET is in accordance with dense vegetation cover, whereas low ET agrees with sparse vegetation cover.It indicates that the SEB-A method can well demonstrate ET differences caused by land use types.This agrees with previous suggestions that estimating ET at a regional level requires the incorporation of heterogeneous surface conditions [40].We used the normalized difference vegetation index (NDVI) as an indicator to describe the vegetation characterization and gave phenology of maize, which is the staple crop in the study area.The summarized NDVIs of maize, vegetables, orchards and woodland, and the phenology of maize, are shown in Tables 4 and 5, respectively.In July and August, maize is in its vigorous growth season (Table 5) with frequent irrigation, and hence, the vegetation cover (Table 4) is dense and the ET is high.Because the orchards, vegetation (pepper, celery and cauliflower) and woodlands had sparse vegetation cover compared to the maize during this time, their associated ETs are lower than those of maize.As for villages, they had largely bare and solidified underlying surfaces, and thus their ETs were the lowest.These results are consistent with the facts.

Relationship between ET Estimates and Land Use Status
Since spatial distributions of ET estimates were obtained (Figure 3), we can conduct deep analyses by combining the land use types (Figure 1).The land use type within the study area, which was based on the aerial remote sensing data of the Compact Airborne Spectrographic Imager (CASI) and used the Support Vector Machine (SVM) for classification, was also provided by the Heihe Plan Science Data Center [39].
There were five land uses types in the study area, including maize growing, village habitats, vegetable crops, orchards and woodlands.For each land use type, we performed statistical computation and gave an average ET.Table 3 shows the average ET of each land use type on the three days.Maize had the maximum ET, and the average ETs on the three days were 528 W/m 2 , 389 W/m 2 and 425 W/m 2 , respectively.Vegetable crops, orchards and woodlands had an average ET between 400 W/m 2 and 500 W/m 2 on 10 July and an average ET of approximately 300 W/m 2 on the other two days.Villages showed the lowest ET, with an average value of 330 W/m 2 on 10 July and less than 200 W/m 2 on the other two days.The spatial distribution of ET was intimately linked with vegetation patterns.In general, high ET is in accordance with dense vegetation cover, whereas low ET agrees with sparse vegetation cover.It indicates that the SEB-A method can well demonstrate ET differences caused by land use types.This agrees with previous suggestions that estimating ET at a regional level requires the incorporation of heterogeneous surface conditions [40].We used the normalized difference vegetation index (NDVI) as an indicator to describe the vegetation characterization and gave phenology of maize, which is the staple crop in the study area.The summarized NDVIs of maize, vegetables, orchards and woodland, and the phenology of maize, are shown in Tables 4 and 5, respectively.In July and August, maize is in its vigorous growth season (Table 5) with frequent irrigation, and hence, the vegetation cover (Table 4) is dense and the ET is high.Because the orchards, vegetation (pepper, celery and cauliflower) and woodlands had sparse vegetation cover compared to the maize during this time, their associated ETs are lower than those of maize.As for villages, they had largely bare and solidified underlying surfaces, and thus their ETs were the lowest.These results are consistent with the facts.

Sensitivity Analysis
The sensitivity analysis plays a significant role in understanding the contribution of individual variables to the model and can give suggestions on applying the model better.In the study, we followed the procedures used in a prior study by Masson.et al. [24] to conduct the sensitivity analysis.Table 6 gives the initial and conventional values of variables as references to examine the sensitivity of the SEB-A method.The direct inputs of Equation ( 5), Rn ´G and LST ´Ta-loc , were analyzed (the LE exo was ignored here because it is added to Equation (2) directly, and the sensitivity of it is similar to R n ´G).Perturbations of Rn ´G were specified as [´20%, 20%], with a variation step of 10%.Perturbations of LST ´Ta-loc were also specified as [´20%, 20%], with a variation step of 10%.Tables 7 and 8 show the deviations in ET computed after perturbing Rn ´G and LST ´Ta-loc from their initial values individually.Rn ´G is positively correlated with ET, but LST ´Ta-loc is negatively correlated with ET.When comparing the sensitivity of Rn ´G with that of LST ´Ta-loc , Rn ´G is found to be more sensitive than LST ´Ta´oc to ET.A 10% increase in Rn ´G could result in a 30 W/m 2 increase in ET estimates, whereas a 20% increase in LST ´Ta-loc only results in a 17.3 W/m 2 decrease in ET estimates.The sensitivity results indicate that the SEB-A method requests accurate estimates of Rn ´G, but it is less constrained with LST ´Ta-loc .

Discussion
The SEB-A method ignores the inner mechanisms of ET; rather, it only concerns itself with the final estimates and thus the principle of it is clear to understand.Besides the principle, the SEB-A approach has another advantage, which uses the exact value of air aerodynamic resistance (150 s/m).Because the air aerodynamic resistance required by the method is under the condition of no advection, its value is exact, which will avoid complex computational processing of it in practice.In previous studies, the air aerodynamic resistance was difficult to determine.In addition, the surface air temperature, T a-loc , needed is under the condition of no advection and can be estimated by remote sensing, which is another situation where the SEB-A method has an edge over previous methods that usually use interpolated air temperature from meteorological stations.Furthermore, the sensitivity analyses indicate that the SEB-A method is insensitive to the temperature difference between the air and the ground.This is promising because determining LST and T a-loc is complex and may result in transferring errors.The insensitivity of LST ´Ta-loc improves the applicability of the SEB-A method.
Despite the advantages of the SEB-A method, it still carries some uncertainties.On the one hand, the approach requires plenty of ground-based observations to calculate the exotic ET driven by horizontal advection; therefore, the lack of sufficient ground-based observations would lead to inaccurate ET estimates.On the other hand, because the exotic ET is calculated by using EC observations, it is ground-based and needs to be extended from a point scale to a regional scale by the IDW.The IDW, however, interpolates according to the distance between inputted objectives with different values (see Equations ( 13) and ( 14)).As a result, the interpolation results may display in the form of circular high-value or low-value areas that do not exist in practice.Thus, the IDW may be inappropriate to handle the exotic ET.It is essential to find out the diffusion mechanism of the exotic ET and propose an appropriate way to address it.
The SEB-A method provides a way to take into account the horizontal advection when applying the energy balance theory and gives avenues for quantifying the exotic ET produced by the horizontal advection.

Conclusions
The study proposes the SEB-A method, which is innovative in that it considers the horizontal advection when applying the energy balance theory to estimate ET through remote sensing.The advantage of the SEB-A method is that it simplifies the ET process and regards the ET as regulated by two different driving forces: the local driving force and the exotic driving force.The local driving force is generated by radiation effects and can be described by the energy balance equation, whereas the exotic driving force is caused by the horizontal advection and can be expressed indirectly through ground-based EC observations.The idea makes the ET retrieval much simpler for the exact air aerodynamic resistance and for the quantification of the exotic ET.Model application and validation indicated that ET estimates from the SEB-A method had high accuracy, with an R 2 of 0.713, a MAE of 39.3 W/m 2 and an RMSE of 54.6 W/m 2 between estimates and measurements.In addition, ET estimates had percent errors within 20% at more than 80% of the locations and could well demonstrate ET differences caused by heterogeneous underlying surfaces.The sensitivity results showed that the SEB-A method requested accurate estimation of Rn ´G, but it was insensitive to the difference between surface and air temperatures.
Since the influence of horizontal advection on ET is ubiquitous in practice, it is necessary to consider the exotic ET caused by the horizontal advection.The SEB-A method puts forward the way of utilizing ground-based EC measurements to quantify the exotic ET, which is verified with high accuracy.If remote sensing data and ground-based measurements including meteorological variables and surface fluxes are available, the SEB-A method can perform well.However, the method cannot work well when the ground-based measurements is insufficient.

Figure 1 .
Figure 1.Location of the study area including land use status and 17 ground-based stations.

Figure 1 .
Figure 1.Location of the study area including land use status and 17 ground-based stations.

Figure 2 .
Figure 2. The principle of the surface energy balance-advection (SEB-A) method.

Figure 4 .
Figure 4. Comparison between ET observations from EC stations and ET estimates by the SEB-A method.

Figure 4 .
Figure 4. Comparison between ET observations from EC stations and ET estimates by the SEB-A method.

Figure 5 .
Figure 5.The frequency distribution of percent errors on ET estimates.

Figure 5 .
Figure 5.The frequency distribution of percent errors on ET estimates.
Because LE act contains both LE loc and LE exo , to estimate LE act , it is necessary to calculate the two parts: LE loc and LE exo .

Table 1 .
r a-loc (s/m) measured at a height of 0.5 m.

Table 2 .
r a-loc (s/m) measured at a height of 1 m.

Table 3 .
ET estimates and the exotic ET retrieved by the SEB-A method on five different land use types on 10 July, 11 August and 27 August 2012.

Table 3 .
ET estimates and the exotic ET retrieved by the SEB-A method on five different land use types on 10 July, 11 August and 27 August 2012.

Table 5 .
The phenology of maize in the study area in 2012.

Table 6 .
The preset initial values of variables in the SEB-A method for the sensitivity analysis.

Table 7 .
Deviations in ET after perturbing Rn ´G.

Table 8 .
Deviations in ET after perturbing LST ´Ta-loc .