Influence of Vegetation Cover on Regional Evapotranspiration in Semi-Arid Watersheds in Northwest China

Evapotranspiration (ET) is an important biospheric process whereby liquid water is vapor‐ ized from moist surfaces and from plant tissues [1]. It is one of the key processes in the hy‐ drological cycle and the energy balance of watersheds. Level of ET from a surface is primarily dependent on the availability of moisture on the surface and the amount of energy available to evaporate that moisture. To understand and quantify ET, three conceptual defi‐ nitions of ET exist in the scientific literature. Actual ET (AET) refers to the actual evapora‐ tion of liquid water from a surface under given atmospheric conditions. For the same atmospheric conditions, potential ET (PET) refers to the elevated evaporation of water when the amount of surface moisture is unlimited and vegetation conditions are ideal. Ideal con‐ ditions regarding vegetation is characterized by actively growing short vegetation covering a large surface area with unlimited supply of soil water [2, 3]. If the vegetation cover is standardized to grass or alfalfa, PET is considered as reference ET (i.e., ETo) [1]. Understand‐ ing spatiotemporal trends in ET is critical to interpreting eco-hydrometeorological processes of complex landscapes [4]. Specific factors affecting ET are divided into three main catego‐ ries, namely:


Introduction
Evapotranspiration (ET) is an important biospheric process whereby liquid water is vaporized from moist surfaces and from plant tissues [1]. It is one of the key processes in the hydrological cycle and the energy balance of watersheds. Level of ET from a surface is primarily dependent on the availability of moisture on the surface and the amount of energy available to evaporate that moisture. To understand and quantify ET, three conceptual definitions of ET exist in the scientific literature. Actual ET (AET) refers to the actual evaporation of liquid water from a surface under given atmospheric conditions. For the same atmospheric conditions, potential ET (PET) refers to the elevated evaporation of water when the amount of surface moisture is unlimited and vegetation conditions are ideal. Ideal conditions regarding vegetation is characterized by actively growing short vegetation covering a large surface area with unlimited supply of soil water [2,3]. If the vegetation cover is standardized to grass or alfalfa, PET is considered as reference ET (i.e., ET o ) [1]. Understanding spatiotemporal trends in ET is critical to interpreting eco-hydrometeorological processes of complex landscapes [4]. Specific factors affecting ET are divided into three main categories, namely: meteorological factors, including solar radiation, near-surface air temperature, wind velocity, and air humidity; ii. surface factors, including surface water and soil water content; and iii. plant factors, including rooting depth, leaf structure, and stomatal density and aperture [1,5,6].

Hydrological methods
Hydrological methods include the soil-water balance and weighing-lysimeter methods [18]. AET-determination from the soil-water balance is possible with a rearrangement of the water-budget equation, such that , AET P R D S = ---D (1) where P is the precipitation, R is the surface runoff, D is the soil drainage, and ∆S is the change in soil-water storage.
Each term in Eqn. (1) is given in unit catchment volume, in which proper delineation of catchment boundaries, distribution of precipitation and streamflow gauges, and accuracy of runoff are critical for the accurate assessment of AET [19].
Lysimeters are measuring devices with containers for soil and plants that estimate ET by measuring the difference in weight of water in individual containers [20]. Accuracy of AET with lysimeters rests on the precision of the weighing instrument and sampling frequency [18].

Micrometeorological methods
Micrometeorological methods estimate latent heat fluxes at potentially high temporal frequencies over homogeneous vegetation [18]. The Bowen ratio method estimates latent heat fluxes [i.e., λET (W m -2 ), where λ is the latent heat of vaporization expressed in J kg -1 ] from the ratio of surface energy to the Bowen ratio (non-dimensional), such that where R n is the net surface radiation (W m -2 ), G is the soil-heat flux (W m -2 ), and β (Bowen ratio) is the ratio to sensible heat and latent heat fluxes estimated from air temperature and water vapor pressure gradients taken within the same sample of air [21], or where γ is the psychometric constant (0.67 hPa K -1 , at sea level) and ∆T a (K) and ∆e a (hPa) are differences in air temperature and water vapor pressure taken at two separate levels [16].
The aerodynamic resistance method uses the atmospheric boundary-layer resistance to calculate the transfer of sensible heat (H; W m -2 ) from the surface to the air [22] and is formally expressed as Influence of Vegetation Cover on Regional Evapotranspiration in Semi-Arid Watersheds in Northwest where r a is the aerodynamic resistance (s m -1 ), ς and C p are the density and specific heat of the airmass involved (kg m -3 and J kg -1 K -1 , respectively), and T s and T a are the surface and air temperature (K). The aerodynamic resistance in Eqn. (4) is the most important and most difficult to define. Various methods have been proposed to estimate this resistance, many of which are summarized in [23]. Once H is determined, AET (i.e., AET = λET . λ -1 ) can be determined from given corresponding measurements of R n and G.
Eddy covariance facilitates determination of vertical fluxes of atmospheric gases (e.g., water vapor, carbon dioxide, and trace gases) and heat within the atmospheric boundary layer. Historical advancement of eddy-covariance concepts is summarized in [17]. Assuming that the vertical wind velocity is responsible for vertical fluxes of atmospheric gases and heat, there should be a high positive correlation between vertical wind velocities and individual fluxes [24]. The method monitors high-frequency fluctuations in vertical wind velocity (m s -1 ) and water vapor mixing ratios (kg kg -1 ), in estimating AET from the covariance of the two, i.e., where w' and q' are the instantaneous fluctuation (deviation from the interval mean) in vertical wind velocity and water vapor mixing ratio in a parcel of air.
The physiological approach determines AET for individual plants or groups of plants. Two of the more broadly used physiological methods, includes 1. the sap-flow method, and 2. the chamber method.
In the sap-flow method, sap flow in plants is assumed to be closely linked to plant transpiration and is quantified by applying heat pulses to the stem of plants and analyzing corresponding heat balances [18]. In a chamber setting, transpiring plants are enclosed in a transparent chamber and changes in within-chamber water vapor concentrations are quantified [25].
While AET-determination methods based on principles of micrometeorology and plant physiology provide suitable overall accuracy at point scales (within a few tens of metres), spatial interpolation of their results across entire landscapes is wholly inappropriate because of the inherent complexity of natural landscapes, particularly with respect to topography and landcover.

Regional estimates of AET from remote sensing-based data
For expansive landscapes, remote sensing (RS)-based methods have been gaining popularity during the past few decades with regard to estimating regional AET at daily, monthly, and annual time scales. Characterizing land-surface conditions with RS-based methods provides an important way of overcoming the difficulty of interpolating AET for complex landscapes [26]. Most methods using earth-observation data in approximating regional AET can be categorized into three main groups, i.e., methods based on: i. the surface-energy balance equation [27][28][29]; ii.
While all methods require an assessment of available net energy at the surface (i.e., R n -G) as primary input, they differ in the way they partition the energy into H and λET.

Surface energy balance
The surface-energy balance is based on the assumption that R n is equal to the sum of H, λET, and G; the energy required for photosynthesis is insignificant here as it accounts for < 1% of incoming solar radiation [35]. Based on a rewriting of the surface-energy balance, the expression of latent heat flux becomes: In Eqn. (8), R s ↓ and R s ↑ represent incoming and outgoing shortwave radiation and R L ↓ and R L ↑ , incoming and outgoing longwave radiation emitted by the atmosphere (including clouds, if present) and earth surface. Typically, R n and G (generally expressed as a fraction of R s ↓ ) are determined from sun-earth geometric relations and illumination angles imposed by variable terrain. Most of the functions and information for the determination of R n and G are available in present day geographic information systems (GIS) and digital terrain models [36][37][38][39]. Reflected and emitted radiative components are estimated from surface albedo, land surface and air temperatures, and emissivities derived from RS-data, e.g., as those identified in Table 1.
Most popular among surface-energy balance models are the surface-energy balance system (SEBS; [29]) and surface-energy balance algorithm for land (SEBAL; [27]). A comprehensive review of various surface-energy balance procedures appearing in the scientitific literature can be found in [40]. The SEBS model uses Monin-Obukhov similarity theory for the atmos-pheric surface layer to derive land surface physical parameters and roughness lengths in determining the evaporative fraction at the limiting ends of accessible water, i.e., dry and wet limits [29]. At the dry end, evaporation (ET dry ) is assumed to be zero due to the limitation in soil moisture and, as a result, sensible heat (H dry ) is maximum, equalling net available energy (i.e., H dry = R n -G). At the wet end, evaporation (ET wet ) is assumed to be at its potential rate (i.e. λET wet = λPET) and sensible heat (H wet ) is minimum. Relative evaporation fraction is defined as: where H is the actual sensible heat flux and determined by solving a system of non-linear equations proposed by Brutsaert [29,41]. AET is subsequently calculated with: where H wet is calculated using an equation proposed by Meneti [42] based on the Penman-Monteith combination equation [43].
The SEBAL model [27] was designed to estimate regional AET from RS-data consisting of surface temperature, surface reflectance, Normalized Difference Vegetation Index (NDVI), and their corresponding relations. The SEBAL model calculates H from a linear relationship between T s and T a derived from plotted distributions of wet and dry pixels.

Penman-Monteith equation
To estimate regional ET using MODIS (or Moderate Resolution Imaging Spectroradiometer) data, Cleugh et al. [44] proposed the use of the Penman-Monteith equation, i.e., where ∆ is the slope of the saturation-vapor-pressure-to-temperature-curve, r a and r s are the aerodynamic (atmospheric) and surface resistances (s m -1 ) to the transfer of surface water vapor to the atmosphere, e s is saturation vapor pressure (hpa), and e a is actual vapor pres-sure (hPa). In their model, input to Eqn. (12) are based on MODIS-derived vegetation data and daily surface meteorological data. Surface conductance is estimated in the model as a function of MODIS LAI and NDVI.
The model was later modified by Mu et al. [5] by adding vapor pressure deficit (VPD) and minimum air temperature as two important factors influencing stomatal conductance. In this modification, EVI (enhanced vegetation index; [45] ) was used in replacing NDVI in the calculation of surface-vegetation fraction. Also, treatment of soil-water evaporation was incorporated. The model was further modified by Mu et al. [46] to include calculation of daytime and nighttime ET, soil-heat flux, and canopy conductance for both dry and wet foliar conditions. This variant of the model along with climate-station data worldwide were used in the production of a global product of AET (i.e., MOD16; [46]) that was subsequently validated against flux data acquired from 46 AmeriFlux sites distributed throughout the Americas.

Complementary-based methods
Complementary-based methods were first introduced by Bouchet [33], considering that whenever a well-watered surface dries, the decrease in AET is coupled with a corresponding increase in PET and, as a result 2 , where E w is defined as the value of potential evaporation, when AET and PET are equal for an unlimited moist surface [47].
The complementary relationship is based on the assumption that the energy required for evaporation is permanently available [48]. In such instance, drying of a wet surface would result in a reduction in AET and a corresponding reduction in energy consumption. The energy saved would lead to an increase in PET [49]. The complementary relationship was extended by Granger [34] for non-saturated surfaces and non-equal changes in AET and PET. He expressed these changes as Granger [34] defined PET, AET, and E w using a Dalton-type mass-transfer equation and introduced an alternative equation for the complementary relationship, i.e., Influence of Vegetation Cover on Regional Evapotranspiration in Semi-Arid Watersheds in Northwest China http://dx.doi.org/10.5772/52812 For non-saturated surfaces, Granger and Gray [50] proposed a relative evaporation fraction with the notion that wind velocity (u) has an equal impact on AET and PET, i.e., where e s s and e a s are the saturated and actual water vapor pressure at the surface (hPa), e a is the actual water vapor pressure of the air (hPa), and f(u) is a function of wind velocity.
where α is the Priestley-Taylor constant (i.e., 1.26) and u 2 is the wind velocity at 2-m above the ground or canopy surface (m s -1 ).
The CRAE model estimates AET from PET and E w , i.e.,

.
w AET E PET = - PET is addressed by decomposing the Penman equation into its two main components: one component to address the energy balance at the surface and another, the transfer of water vapor from a moist surface, or and where e s p is saturation vapor pressure at T p and T p is the equilibrium temperature. The equilibrium temperature is defined as the temperature at which the energy balance [Eqn. (19)] and mass transfer [Eqn. (20)] give the same result for PET, when solved numerically. Parameter f T is a vapor-transfer coefficient considered constant for a specific pressure level and independent of wind velocity. Parameter f P is a function of atmospheric pressure (Ψ; hPa) and T p , such that E w is calculated with a regression equation, where b 1 and b 2 are coefficients determined by regression, ∆ TP is the slope of the saturationvapor-pressure-curve at T p , and R nTP is the net surface energy at T p , ε s is the surface emissivity (non-dimensional), and σ is the Stefan-Boltzmann constant (i.e., 5.67 × 10 -8 W m -2 K -4 ).
The two models were validated by Hobbins et al [53] using long term water-balance estimates from 120 watersheds distributed throughout the United States. Outcomes of this validation have revealed that the AA-model underestimated the annual AET by 10.6% of annual precipitation and the CRAE-model overestimated the same variable by 2.5% of annual precipitation. For highly arid watersheds, both models tended to overestimate AET. Validation of complementary-based models in other regions of the world, as well as comparing their results with direct observations of AET and PET strongly support the prevalence of the complementary relationship. Ramirez et al [54] found that complementary-based models tended to work well for temperate, humid regions of the world and less well for arid regions [55].  T T  AET  EF PET T T  -= = - The hypothetical temperature (T u ) is derived using the slope of the water vapor pressure curve at T d (i.e., Δ d ) and T s (Δ s ) and the relationship between water vapor pressures and corresponding temperatures, such that Combining Eqn.'s (15) and (23) and the Priestley-Taylor equation [32] with the definition of E w , the equation of AET then becomes These changes render the method less error prone and simpler regarding its data requirements. The method is numerically robust, producing errors < 15% of daily AET-measurements for heterogeneous landscapes of the Southern Great Plains of the United States [56].
Kalma et al. [57] summarized results of validation studies of 30 RS-based AET methods, where it was found that Venturini et al.'s [56] method produced the highest accuracy among the methods using MODIS-derived data. In the current study (Section 5.2), we use the Venturini et al. [56] method to calculate regional AET for the extremely complex landscape of westcentral Gansu, NW China.

Data requirements and accessibility
Data requirements in the calculation of AET depend on the method or model used in the calculation procedure. In general, input data for all methods can be grouped into four main categories: i. surface meteorological data, including near surface air temperature, wind velocity, and in-air water vapor pressure and humidity; ii.
radiative energy fluxes, including incoming shortwave ( R s ↓ ) and longwave radiation ( R L ↓ ) and their outgoing counterparts (i.e., R s ↑ and R L ↑ ); iii. surface attributes, including temperature, emissivity, albedo, and soil water content; and iv. vegetation attributes, including LAI, vegetation-cover density and extent, and stomatal density and aperture.
These data can be i. acquired by direct measurement in the field,

ii.
derived from other related variables by means of regression, or iii. derived from RS optical reflectance or thermal emission data.
While optical and thermal RS has been providing spatially-distributed values for the assessment of landscape AET, it poses several important challenges. RS-images from LANDSAT, ASTER, and SPOT-systems have high spatial resolution, but are acquired infrequently (i.e., every 16 days for LANDSAT and ASTER-systems and every 26 days for SPOT). In contrast, high temporal-acquisitions with geostationary satellites provide coarse spatial resolutions that are largely inadequate for regional-estimate of AET.
The MODIS-instrument is crucial to the production of important global image-products for land-and ocean-surface monitoring [58]. MODIS-instruments were launched to space by NASA (National Aeronautics and Space Administration) as part of the earth-observing system (EOS). MODIS Terra and Aqua image-products are available since February 24, 2000 and June 24, 2002, respectively [59]. Raw image-products have a temporal-acquisition interval of one or two days, covering 36 spectral bands between 0.405 and 14.385 μm [59,60]. Besides daily image scenes, several composite-products are also available at multiple-day intervals at different processing levels. Level-1 products provide the instrument data at full resolution for individual MODIS scenes, which can be exploited after converting to radiance values, surface reflectance, or brightness temperature. Level-2, 2G, and 3 images contain derived geo-biophysical parameters, while Level-4 images contain model-generated products. Level-2G and -3 image-products consist of aggregations of daily images to give 8 and 16-day cloud-free images. Table 1 lists a number of MODIS image-products valuable in estimating regional AET.  Table 1. Data needed in the calculation of regional AET and corresponding MODIS image-product sources and spatiotemporal resolutions.

Variables MODIS image-products Spatiotemporal
Influence of Vegetation Cover on Regional Evapotranspiration in Semi-Arid Watersheds in Northwest China http://dx.doi.org/10.5772/52812

Land surface temperature
MODIS land surface temperatures (T s ; MOD11A2) are estimated from MODIS thermal infrared (TIR) data collected in near-cloud-free conditions with the application of a split-window algorithm to avoid misidentification with cloud-top conditions [61]. The algorithm uses MODIS-emissivity bands 31 and 32 in the calculation of T s [62]. Details of the calculationprocedure can be found in [63].

Vegetation indices
MODIS vegetation-products (MOD13Q1) provide two vegetation-indices (VI) both at 250-m and 16-day resolution [64]. NDVI is calculated as the normalized ratio of near infrared (NIR) and red bands. EVI (introduced earlier) was developed to improve the sensitivity of NDVI for high biomass regions [45]. NDVI and EVI are calculated from and where ρ's are wavelength-specific atmospherically-corrected or partially atmosphere-corrected (for Rayleigh and ozone absorption) surface reflectances; B is the gain factor; L is the canopy background adjustment which addresses nonlinear, differential NIR and red-radiant transfer through vegetative canopies; and C 1 , C 2 are the coefficients of aerosol resistance, which uses the blue band to correct for aerosol influences on the red band.

Land surface albedo
Land surface albedo is the fraction of shortwave radiation reflected in all directions and is a critical parameter in estimating surface net shortwave radiation [66]. MODIS-BRDF (bidirectional reflectance distribution function)-based albedo image-products combine atmospherically-corrected surface reflectance from multiple dates and sensors in creating 16-day, 1-km resolution images [67]. The MODIS product, MCD43B3, includes both black and white sky albedo acquired with the Terra and Aqua satellites [68], covering all seven spectral bands and three broad bands [60].

Air and dew point temperature
Air and dew point temperatures provided by the MODIS-MOD07 product (at 5-km resolution) is currently the best possible estimate of the two variables derived from RS methods [69]. MOD07 atmospheric profile information includes temperature and atmospheric humidity at 20 different vertical pressure levels between 5 to 1000 hPa [70].

Study area description
The study area consists of the Shiyang and Hei River watersheds in westcentral Gansu, NW China (Fig. 1). The Shiyang River originates from the Qilian Mountains and flows northwestward before terminating in the Minqin-lake district [71]. The total basin area is approximately 49,500 km 2 . Elevation in the Shiyang River basin varies from 1,284 to 5,161 m above mean sea level (AMSL), with an average elevation of 1,871 m AMSL. The Hei River watershed, with a land surface area of approximately 128,000 km 2 , is the second largest inland river basin in NW China [72]. The Hei River watershed includes the Zhangye watershed, with a total land area of about 31,100 km 2 . Elevation in the Zhangye watershed varies from 1,287 to 5,045 m AMSL, with an average elevation of 2,679 m AMSL (Fig. 1).
The natural landscape of the study area comprises of mountains, oases, and deserts, all interacting with each other [73]. The oases play an important role in the sustainability of the region's overall ecological and socio-economic integrity [74]. The area overlaps four distinct ecoregions [75]. The northern part, noted for its arid to semi-arid conditions, includes portion of the Badain Jaran and Tengger deserts and oases in the southwest portions of the Alashan Plateau. Liangzhou, at the south, and Minqin, at the north, are two important oases in the Shiyang River watershed [71]. Zhangye is the main oasis in the Zhangye watershed. Spring wheat is the main food crop grown in the oases, which is usually supported by irrigation [76]. In the deserts, salt-tolerant, xerophytic shrub species, i.e., saxaul (Haloxylon ammodendron) and Reaumuria soongolica [77] are common.
Locally-generated rainfall in the oases is normally insufficient (< 170 mm yr -1 ) to support agriculture in the oases [7]. The main source of water to the oases is the runoff from the Qilian Mountains generated from snowmelt and in-mountain precipitation transported by the Shiyang and Hei River systems [78,79]. Glacial meltwater contributes to about 3.8% and 8.3% of total runoff in the Shiyang and Hei Rivers, respectively [78]. The meltwater usually flows during the spring-summer period due to the warming of the mountain glaciers and previous snow-season's snow cover [74].
A primary source of water in the rivers during the summer is orographic precipitation [80] formed in the Qilian Mountains [81]. Annual precipitation in the watersheds is < 80 mm in the desert-portion of the watershed and > 800 mm in the Qilian Mountains, proper. Annual PET in the deserts ranges between 2,000 to 2,600 mm and 700 to 1,200 mm in the mountains [82]. Most of the precipitation occurs during June-August. About 94% of the water delivered from the mountains is through surface runoff. Average annual runoff delivered by the Shiyang River is about 15.8 ×10 8 m 3 and about 37.7 × 10 8 m 3 by the Hei River [82]. Presence of the Qinghai-Tibetan Plateau at the south of the study area blocks the northward passage of southwest monsoonal precipitation and the westerly airflow on the northern side of the plateau interferes with southerly airflow from reaching the region [83]. Dry northwesterly winds during summer generated from the Azores high pressure system and cold dry northerly winds during the winter generated from the Siberian high pressure system limit the outside contribution of moisture to this area [84]. Cyclic exchange of atmospheric and surface and shallow subsurface water in the region forms a nearly perfect close system [74]. ET in the series of oases at the base of the Qilian Mountains is shown to play an important role in the recycling of water in the region and seasonal evolution of snow cover in the Qilian Mountains [7,85].

Landcover classification
The MODIS annual global landcover map currently available (as of 2012) is produced from seven spectral maps, BRDF-adjusted reflectance, T s , EVI, and an application of supervised classification using ground data from 1860 field sites [86]. Assessments of the product have shown that this map is not entirely realistic for zones of transition or for mountainous regions [87]. Improved landcover definition at regional or local scales with supervised classification usually involves much greater amounts of ground data that are normally available for most regions. Recently, decision tree-based classification has been applied to RS-data and has been shown to produce better results than other classification systems based on maximum likelihood or unsupervised clustering and labeling [88]. One advantage of decision tree-based classification is that it is able to use local knowledge of vegetation characteristics together with other pertinent data, such as terrain characteristics. In the current study, we use chronological-sequences of MODIS-based EVI and digital terrain information (e.g., aspect, elevation) to classify landcover with decision trees.
Vegetation distribution in the study area has a unique preferential association with elevation, slope, and slope direction [89]. North-facing slopes of the Qilian Mountains support alpine meadow at elevations between 2,500 to 3,300 m AMSL. At elevations above 3,300 m AMSL, deciduous shrubs represent the most dominant vegetation type. Isolated patches of conifer forests in the Qilian Mountains are classified as a separate ecoregion [90] found at elevations between 2,500 m to 3,300 m AMSL. Vegetation density and seasonal vegetation growth vary as a function of vegetation type and, consequently, landcover.
Based on vegetation site preferences, the study area is subdivided into four main elevation zones defined by elevations: i. < 2,500; ii. between 2,500 to 3,300; iii. between 3,300 to 3,900; and iv. > 3,900 m AMSL.
Different landcover types in these elevation zones were then identified based on EVI and terrain attributes, in particular slope orientation (aspect). Landcover types and their discrimination are summarized in Table 2.  Table 2. Landcover definition as a function of elevation zone, EVI, and slope orientation (aspect). Ten landcover maps were generated for 2000-2009 using the classification standards summarized in Table 2. From these maps, a final landcover composite (LCOV dom for all image pixels; Fig. 2) was then created based on a pixel-level, landcover-dominance evaluation, i.e., where LCOV i and LCOV dom represent landcover at the pixel-level for individual years (2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008)(2009)) and the dominant landcover over the same ten-year time period, respectively.

Study-area AET variation for different landcover types
Monthly AET for the study area was calculated using the complementary method of Venturini et al. [56] for the period of 2000-2009. Computed AET was compared with ET o calculated from pan evaporation data (Fig. 3) corrected with season-specific coefficients reported in [91].  Fig. 3 show that modeled AET has very high positive correlation with ET o for both watersheds (R 2 > 0.78; Fig.  3a for the Shiyang and Fig. 3b for the Hei River watershed). The line and error graphs in Fig.  3 are generated from AET-values extracted at the location of four meteorological stations within the Liangzhou, Minqin, and Zhangye oases (Fig. 1). The lines represent mean total monthly AET based on 2000-2009 data; error bars represent standard deviation of corresponding values. Modeled AET display seasonal patterns similar to those expressed in ET o , but at substantially reduced levels (Fig. 3). This discrepancy can be rationalized by three important realities, i.e., iii. pan evaporation coefficients are based on seasonal averages independent of location.
Spatially-distributed average total monthly AET was calculated by averaging monthly AETimages generated for the 2000-2009 period (Fig. 4). Yearly total growing-season AET was calculated by summing AET from April-October of each year (Fig. 5). From Fig. 4, it is clear that AET in the watershed is very low (≤ 25 mm month -1 ) during winter (January-March and October-December periods) and progressively higher in summer (> 75 mm month -1 ). AET reaches its maximum during the June-August period and begins to decrease prior the start of winter. Both monthly and total growing-season AET reveals greatest AET in the oases and low-to-mid-slope positions of the Qilian Mountains and lowest in the deserts (not shown) and high-elevation portions of the Qilian Mountains. Within-year variation in AET also appears within the same elevation bands due to changes in vegetation type. Comparison of chronological-series of monthly AET for three landcover types is shown in Fig. 6. For most years, forest landcover is shown to contribute the most to AET in the Shiyang River watershed; crop cover, however, contributes the most in the Hei River watershed (Fig. 6). For the two watersheds, areas sparse of vegetation (as defined in Table 2) consistently contribute the least to AET during the growing seasons of 2000-2009.   Influence of Vegetation Cover on Regional Evapotranspiration in Semi-Arid Watersheds in Northwest China http://dx.doi.org/10.5772/52812

Vegetation influence on AET
Using EVI as an indicator of vegetation density and vitality, monthly AET was compared with same-month EVI extracted at 220 randomly distributed points across the study area for different landcover types. Scattergraphs using all points and partitioned according to landcover type (Fig. 7) show that independent-evaluations of AET have an overall positive correlation with same-month EVI; the strength of correlation (R 2 ), however, varies with landcover type. Strongest correlation occurs for cropland areas and weakest, for areas sparse of vegetation.

Conclusions
This chapter provides i. a summary of methodologies available to estimate AET at point to regional scales, and ii. case-study calculations of AET for two watersheds in NW China based on the complementary relationship of Venturini et al. [56].
The complementary relationship of Venturini et al. is considered appropriate for complex landscapes of NW China, because of the method's independence from wind velocity, its overall accurracy, and its ability to regionalize AET-calculations with assistance of RS-data as input. Because of the frequent acquisition of MODIS data (primary input datatype to the calculation of AET), AET-calculations can be updated frequently. For NW China, AET-calculations at 250-m resolution were carried out on a monthly interval over a ten-year time period (2000-2009). Based on a landcover map generated from decision-tree classification and landcover-dominance analysis, regional AET was partitioned along three vegetation-dominated landcover types. Forest and crop landcover types were shown to contribute the most to AET across the study area, particularly in lowland areas. Areas of sparse vegetation (among the three landcover types) contributed the least to regional AET. This supports the view that the state and abundance of vegetation (defined here by EVI), particularly in the lowlands and low-to-mid-slope positions of the Qilian Mountains, have an important influence on regional AET and on the water budget of the study area.