Data needed in the calculation of regional AET and corresponding MODIS image-product sources and spatiotemporal resolutions.
Evapotranspiration (ET) is an important biospheric process whereby liquid water is vaporized from moist surfaces and from plant tissues . 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., ETo) . Understanding spatiotemporal trends in ET is critical to interpreting eco-hydrometeorological processes of complex landscapes . 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;
surface factors, including surface water and soil water content; and
Spatiotemporal variability in these factors determines the variability of ET across large areas. This variability increases as the complexity of the underlying terrain and vegetation cover intensifies.
Back and forth exchange of water vapor and liquid water from oases at the base of the Qilian Mountains (Northwest China) and from the Qilian Mountains to oases as surface and shallow subsurface flow has been previously shown to be a potentially significant mechanism in the long term maintenance of oases in westcentral Gansu, NW China . In general, direct precipitation (rain + snow) to the foothills and oases at the base of the Qilian Mountains is inadequate to maintain vegetation in these areas. Maintenance of oases vegetation is shown to be dependent on the surface water flowing from the higher regions of the mountain range, where direct precipitation and snowmelt are greatest. Supply of atmospheric moisture that leads to the formation of precipitation in the mountains is sustained by seasonal ET at the base of the mountain range.
Influence of vegetation and landcover on local and regional climate is well documented in the scientific literature [7-9]. Changes in landcover result in changes in surface albedo, surface roughness, leaf area index (LAI), stomatal conductance, rooting depth, and soil texture and structure [10-12]. Changes in surface vegetation have been known to impact
the partitioning of net radiative energy into sensible and latent heat fluxes, and
the formation of convective rainfall by affecting the state of the convective boundary layer .
Agricultural landuse affects climate by modifying the physiological attributes (e.g., canopy conductance) of the land .
2. AET-determination methods
2.1. Point measurements
Suitability of AET-measurement methods depend on the reasons for assessment and on the spatiotemporal dimensions of the problem. Measurement of AET is mainly done indirectly by measuring the effect of AET on the water and energy balance. Rose and Sharma  categorized AET-measurement methods into three broad categories, i.e.,
hydrological approaches, based on the residual of the water-budget equation;
plant-physiological approaches at the scale of individual plants or groups of plants by collecting sap-flow or chamber measurements.
2.1.1. Hydrological methods
Hydrological methods include the soil-water balance and weighing-lysimeter methods . AET-determination from the soil-water balance is possible with a rearrangement of the water-budget equation, such that
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 .
Lysimeters are measuring devices with containers for soil and plants that estimate ET by measuring the difference in weight of water in individual containers . Accuracy of AET with lysimeters rests on the precision of the weighing instrument and sampling frequency .
2.1.2. Micrometeorological methods
Micrometeorological methods estimate latent heat fluxes at potentially high temporal frequencies over homogeneous vegetation . 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 Rn 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 , or
where γ is the psychometric constant (0.67 hPa K-1, at sea level) and ∆Ta (K) and ∆ea (hPa) are differences in air temperature and water vapor pressure taken at two separate levels .
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  and is formally expressed as
where ra is the aerodynamic resistance (s m-1), ς and Cp are the density and specific heat of the airmass involved (kg m-3 and J kg-1 K-1, respectively), and Ts and Ta 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 . Once H is determined, AET (i.e., AET = λET.λ-1) can be determined from
given corresponding measurements of Rn 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 . 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 . 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
the sap-flow method, and
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 . In a chamber setting, transpiring plants are enclosed in a transparent chamber and changes in within-chamber water vapor concentrations are quantified .
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.
2.2. 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 . Most methods using earth-observation data in approximating regional AET can be categorized into three main groups, i.e., methods based on:
While all methods require an assessment of available net energy at the surface (i.e., Rn-G) as primary input, they differ in the way they partition the energy into H and λET.
2.2.1. Surface energy balance
The surface-energy balance is based on the assumption that Rn 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 . Based on a rewriting of the surface-energy balance, the expression of latent heat flux becomes:
In Eqn. (8), and represent incoming and outgoing shortwave radiation and and, incoming and outgoing longwave radiation emitted by the atmosphere (including clouds, if present) and earth surface. Typically, Rn and G (generally expressed as a fraction of) are determined from sun-earth geometric relations and illumination angles imposed by variable terrain. Most of the functions and information for the determination of Rn and G are available in present day geographic information systems (GIS) and digital terrain models [36-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; ) and surface-energy balance algorithm for land (SEBAL; ). A comprehensive review of various surface-energy balance procedures appearing in the scientitific literature can be found in . The SEBS model uses Monin-Obukhov similarity theory for the atmospheric 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 . At the dry end, evaporation (ETdry) is assumed to be zero due to the limitation in soil moisture and, as a result, sensible heat (Hdry) is maximum, equalling net available energy (i.e., Hdry = Rn - G). At the wet end, evaporation (ETwet) is assumed to be at its potential rate (i.e. λETwet= λPET) and sensible heat (Hwet) is minimum. Relative evaporation fraction is defined as:
The SEBAL model  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 Ts and Ta derived from plotted distributions of wet and dry pixels.
2.2.2. Penman-Monteith equation
To estimate regional ET using MODIS (or Moderate Resolution Imaging Spectroradiometer) data, Cleugh et al.  proposed the use of the Penman-Monteith equation, i.e.,
where ∆ is the slope of the saturation-vapor-pressure-to-temperature-curve, ra and rs are the aerodynamic (atmospheric) and surface resistances (s m-1) to the transfer of surface water vapor to the atmosphere, es is saturation vapor pressure (hpa), and ea is actual vapor pressure (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.  by adding vapor pressure deficit (VPD) and minimum air temperature as two important factors influencing stomatal conductance. In this modification, EVI (enhanced vegetation index;  ) 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.  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; ) that was subsequently validated against flux data acquired from 46 AmeriFlux sites distributed throughout the Americas.
2.2.3. Complementary-based methods
Complementary-based methods were first introduced by Bouchet , considering that whenever a well-watered surface dries, the decrease in AET is coupled with a corresponding increase in PET and, as a result
where Ew is defined as the value of potential evaporation, when AET and PET are equal for an unlimited moist surface .
The complementary relationship is based on the assumption that the energy required for evaporation is permanently available . 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 . The complementary relationship was extended by Granger  for non-saturated surfaces and non-equal changes in AET and PET. He expressed these changes as
Granger  defined PET, AET, and Ew using a Dalton-type mass-transfer equation and introduced an alternative equation for the complementary relationship, i.e.,
For non-saturated surfaces, Granger and Gray  proposed a relative evaporation fraction with the notion that wind velocity (u) has an equal impact on AET and PET, i.e.,
where and are the saturated and actual water vapor pressure at the surface (hPa), ea is the actual water vapor pressure of the air (hPa), and f(u) is a function of wind velocity.
Two of the most widely used models based on the complementary relationship are
the advection-aridity (AA) model of Brutsaert and Stricker  and
the complementary relationship aerial evaporation (CRAE) model of Morton .
The AA-model combines
the Priestley-Taylor equation ,
the Penman equation ,
an empirical wind-velocity function for relative evaporation , and
the complementary relationship, giving rise to
where α is the Priestley-Taylor constant (i.e., 1.26) and u2 is the wind velocity at 2-m above the ground or canopy surface (m s-1).
The CRAE model estimates AET from PET and Ew, i.e.,
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
whereis saturation vapor pressure at Tp and Tp 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 fT is a vapor-transfer coefficient considered constant for a specific pressure level and independent of wind velocity. Parameter fP is a function of atmospheric pressure (Ψ; hPa) and Tp, such that
Ew is calculated with a regression equation,
where b1 and b2 are coefficients determined by regression, ∆TP is the slope of the saturation-vapor-pressure-curve at Tp, and RnTP is the net surface energy at Tp, ε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  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  found that complementary-based models tended to work well for temperate, humid regions of the world and less well for arid regions .
One important constraint of conventional complementary-based methods is the requirement of wind velocity in the calculation of relative evaporation, which is not always available or reliable. To overcome this limitation, Venturini et al  modified the relative evaporation fraction (EFr) of Eqn. (16) by replacing the expression of water vapor pressure with the corresponding temperature used in calculating it. In the case of actual surface water vapor pressure, a hypothetical temperature is used. This hypothetical temperature (Tu) is defined as the temperature at which the surface becomes saturated without changing water vapor pressure; this is equivalent to the definition of dew point temperature for air (Td; K). Using these temperatures, Venturini et al.  modified Eqn. (16) to the following:
The hypothetical temperature (Tu) is derived using the slope of the water vapor pressure curve at Td (i.e., Δd) and Ts (Δs) and the relationship between water vapor pressures and corresponding temperatures, such that
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 . Kalma et al.  summarized results of validation studies of 30 RS-based AET methods, where it was found that Venturini et al.'s  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.  method to calculate regional AET for the extremely complex landscape of westcentral Gansu, NW China.
3. 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:
surface meteorological data, including near surface air temperature, wind velocity, and in-air water vapor pressure and humidity;
radiative energy fluxes, including incoming shortwave () and longwave radiation () and their outgoing counterparts (i.e., and);
surface attributes, including temperature, emissivity, albedo, and soil water content; and
vegetation attributes, including LAI, vegetation-cover density and extent, and stomatal density and aperture.
These data can be
acquired by direct measurement in the field,
derived from other related variables by means of regression, or
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 . 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 . 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.
|Variables||MODIS image-products||Spatiotemporal Resolution|
|Land surface temperature (Ts)|
Land surface emissivity (εs)
|MODIS land surface temperature and emissivity data (MOD11A2)||1000||8|
|Normalized difference vegetation index (NDVI)|
Enhanced vegetation index (EVI)
|MODIS vegetation indices (MOD13Q1)||250||16|
|Land surface albedo (As)||MODIS products combined with BRDF-albedo products||1000||16|
|Air temperature (Ta)|
Dew point temperature (Td)
|MODIS atmospheric profile data (MOD07)||5000||1|
3.1. Land surface temperature
MODIS land surface temperatures (Ts; 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 . The algorithm uses MODIS-emissivity bands 31 and 32 in the calculation of Ts . Details of the calculation-procedure can be found in .
3.2. Vegetation indices
MODIS vegetation-products (MOD13Q1) provide two vegetation-indices (VI) both at 250-m and 16-day resolution . 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 . NDVI and EVI are calculated from
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 C1, C2 are the coefficients of aerosol resistance, which uses the blue band to correct for aerosol influences on the red band.
Coefficients adopted in the EVI-algorithm are L=1, C1=6, C2 = 7.5, and G= 2.5 . While NDVI shows some problems with dense forests , both vegetation indices provide similar accuracy when used for cropland delineation .
3.3. 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 . 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 . The MODIS product, MCD43B3, includes both black and white sky albedo acquired with the Terra and Aqua satellites , covering all seven spectral bands and three broad bands .
3.4. 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 . MOD07 atmospheric profile information includes temperature and atmospheric humidity at 20 different vertical pressure levels between 5 to 1000 hPa .
4. 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 . The total basin area is approximately 49,500 km2. 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 km2, is the second largest inland river basin in NW China . The Hei River watershed includes the Zhangye watershed, with a total land area of about 31,100 km2. 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 . The oases play an important role in the sustainability of the region’s overall ecological and socio-economic integrity . The area overlaps four distinct ecoregions . 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 . 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 . In the deserts, salt-tolerant, xerophytic shrub species, i.e., saxaul (Haloxylon ammodendron) and Reaumuria soongolica  are common.
Locally-generated rainfall in the oases is normally insufficient (< 170 mm yr-1) to support agriculture in the oases . 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 . The meltwater usually flows during the spring-summer period due to the warming of the mountain glaciers and previous snow-season’s snow cover .
A primary source of water in the rivers during the summer is orographic precipitation  formed in the Qilian Mountains . 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 . 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 ×108 m3 and about 37.7 × 108 m3 by the Hei River .
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 . 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 . Cyclic exchange of atmospheric and surface and shallow subsurface water in the region forms a nearly perfect close system . 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].
5. Seasonal and annual variation of AET for different landcover types
5.1. Landcover classification
The MODIS annual global landcover map currently available (as of 2012) is produced from seven spectral maps, BRDF-adjusted reflectance, Ts, EVI, and an application of supervised classification using ground data from 1860 field sites . Assessments of the product have shown that this map is not entirely realistic for zones of transition or for mountainous regions . 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 . 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 . 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  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:
between 2,500 to 3,300;
between 3,300 to 3,900; and
> 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.
|Elevation zone||Landcover||Discrimination criteria|
(< 2,500 m AMSL)
|Desert||Mean growing-season EVI < 0.1130|
|Crop||Maximum growing-season EVI > 0.27 and mean growing-season EVI < 0.1130|
|Dense grass||Maximum growing-season EVI > 0.27, but different from cropland|
|Sparse grass and/or shrub||Mean growing-season EVI between 0.113-0.27|
(2,500-3,300 m AMSL)
|Alpine meadow||Maximum growing-season EVI > 0.27 on north-facing slopes|
|Coniferous forest||Maximum growing-season EVI > 0.27 on other than north-facing slopes|
|Sparse grass and/or shrub||Mean growing-season EVI between 0.113 – 0.27|
|Bare land||Mean growing-season EVI < 0.113|
(3,300-3,900 m AMSL)
|Deciduous shrub||Maximum growing-season EVI > 0.27|
|Bare land||Mean growing-season EVI < 0.113|
(> 3,900 m AMSL)
|Sparse shrub||Mean growing-season EVI between 0.113-0.27|
|Snow and/or ice||Mean growing-season EVI < 0.113|
Ten landcover maps were generated for 2000-2009 using the classification standards summarized in Table 2. From these maps, a final landcover composite (LCOVdom for all image pixels; Fig. 2) was then created based on a pixel-level, landcover-dominance evaluation, i.e.,
where LCOVi and LCOVdom represent landcover at the pixel-level for individual years (2000-2009) and the dominant landcover over the same ten-year time period, respectively.
5.2. Study-area AET variation for different landcover types
Monthly AET for the study area was calculated using the complementary method of Venturini et al.  for the period of 2000-2009. Computed AET was compared with ETo calculated from pan evaporation data (Fig. 3) corrected with season-specific coefficients reported in . Pan evaporation coefficients were calculated by relating 50 years of pan evaporation data collected at 580 climate stations distributed across China to ETo calculated for the same stations with the FAO Penman-Monteith equation. Pan evaporation coefficients by Chen et al.  varied from 0.45 to 0.54 from spring to winter. The scattergraphs in Fig. 3 show that modeled AET has very high positive correlation with ETo for both watersheds (R2> 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 ETo, but at substantially reduced levels (Fig. 3). This discrepancy can be rationalized by three important realities, i.e.,
modeled AET represents the average conditions of individual image pixels (covering 62,500 m2) and not points;
surface moisture and vegetation conditions and water-vapor-transfer resistances are openly different from those driving ETo (i.e., open water vs. soil-plant environments); and
pan evaporation coefficients are based on seasonal averages independent of location.
Spatially-distributed average total monthly AET was calculated by averaging monthly AET-images 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.
5.3. 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 (R2), however, varies with landcover type. Strongest correlation occurs for cropland areas and weakest, for areas sparse of vegetation.
This chapter provides
a summary of methodologies available to estimate AET at point to regional scales, and
case-study calculations of AET for two watersheds in NW China based on the complementary relationship of Venturini et al. .
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.