Remote Sensing Distribution of Actual Evapotranspiration over Qaidam Basin, an Arid Area in China

Evapotranspiration is a major loss flux of the water balance in arid and semi-arid areas. The estimation of actual evapotranspiration has significance for hydrological and environmental purposes. The Surface Energy Balance System (SEBS) algorithm was applied to estimate actual evapotranspiration in the Qaidam Basin and its eight hydrological sub-regions, Northwest China. There were 3,036 cloud-free and atmospherically corrected MODIS satellite images from 2001 to 2011 used in the SEBS algorithm to determine the actual evapotranspiration. The result indicated that the estimated annual actual evapotranspiration of the basin increased with time and the value varied from 72.7 to 182.3 mm. SEBS estimates were 7.5% and 14.1% of observed pan evaporation over the western and eastern areas, respectively. The variation of SEBS actual evapotranspiration is influenced by climate factors, vegetation, net radiation, land cover type and water table depth. The analysis of the evaporative behavior of different land cover types in the basin presented that water bodies, marsh, and farmland had relatively higher mean actual evapotranspiration though these land cover types make up less than 3.5% of the total basin. Bare soil has very low evapotranspiration and covered almost 60% of the study area. The actual evapotranspiration was observed to be decreased with an increase of water table depth. Overall, the SEBS algorithm proved to be useful and has potential for estimating spatial actual evapotranspiration on a regional scale.


Introduction
Evapotranspiration (ET) is the main part of water balance and plays an important role in the hydrological process.The surface ET on a regional scale can substantially influence the amount and spatial distribution of water resources.In an arid inland basin, ET is the major loss variable in water budget.It varies with land surface and local meteorological conditions, and it is difficult to be quantified [1].Therefore, quantitative estimation of surface ET is essential for understanding the hydrological cycle, water management, and sustainable utilization of water resources [2][3][4][5].
Qaidam Basin is an inland basin rich in mineral resources in northwestern China.However, the basin is a water shortage area due to a dry climate and fast economic development.Sustainable development of Qaidam Basin will mostly depend on the efficient utilization and strict management of water resources, which require a relatively accurate estimation of ET on a regional scale.
The ET estimation based on remote sensing techniques has developed rapidly since thermal infrared data was first used for evapotranspiration prediction [6].With the advantage of frequent updates and large area coverage, the remote sensing data has been used for ET estimation in many approaches [7][8][9].The satellite images can provide regional scale land surface information effectively and can be used to estimate ET through a remote sensing method without quantifying other hydrological parameters [10].Since ET cannot be measured directly by remote sensing, some models of indirect estimation of ET have been proposed for quantifying the energy balance parameters such as net radiation, soil heat flux, sensible heat flux, and then ET estimation [11][12][13][14][15]. Based on Gowda et al. [16] and Kalma et al. [17], the evapotranspiration estimation models using satellite data can be divided into two classes: (1) Semi-physical or physical approaches; (2) Empirical relationships for estimation of evapotranspiration.As Kustas and Anderson [18] stated, the physically based surface energy balance (SEB) models, especially SEBS, are less site-specific and require little subjective intervention, as opposed to techniques requiring the selection of hot and cold endpoints in the scene [14,19,20].Sun et al. [21,22] further put forward a simplified land surface model to calculate ET using the parameters estimated by a new stationarity-based method which couples the energy balance and water.The parameter estimated by the stationarity-based method can avoid calibration with field observation ET and can be applied broadly in a harsh area which has no continuous forcing (precipitation and radiation) and surface states (soil moisture and temperature).These physical models can also be served as a useful tool when time series of evapotranspiration is needed.
Physically based SEB models can be divided into single-source (e.g., SEBS, SEBAL) and two-source (e.g., TSEB) bulk transfer equations for sensible heat flux estimation [15,23].Many studies have indicated that both single-source and two-source models can provide good results for the surface energy balance partitioning at different scales [24][25][26].The single-source models were used relatively widely in hydro-meteorological fields due to their simplicity.For example, the SEBAL model can estimate actual ET by computing the different variables of the energy balance of the land surface [27][28][29].The modified single-source models such as SEBS can adjust the differences between aerodynamic and radiometric surface temperatures in partial vegetation cover by employing physical models of soil-canopy heat exchange instead of separation of canopy and soil as two source terms [30].Further, the modified surface energy budget technique can also be applied to map the distribution of lake evaporation and wetland ET [31].
The final aim of deriving the evapotranspiration by remote sensing is to reach good operational application under different land cover conditions [32,33].The SEBS model has been validated in operational utility of evapotranspiration deriving for sparsely vegetated and dry areas [34,35], and thus it was selected as the main method for regional evapotranspiration estimation in this study.
The main objective of this study is to estimate the ET of the Qaidam Basin in the Tibet Plateau and to discuss its controlling factors.The specific objectives are: (1) to estimate the ET a of the Qaidam Basin using MODIS land products; (2) to discuss the controlling factors of regional ET a , such as climate factors, vegetation, land cover types, and water table depth.

Study Area
Qaidam Basin stands between 34°40′N and 39°20′N and 90°00′-99°20′E located in the northern part of Tibet Plateau, China.The total area is 276,233 km 2 and the plain area covers 152,393 km 2 .The basin is typically classified as a continental arid climate, where the annual precipitation spatially varies between 16 and 190 mm from west to east and is concentrated from April to November.The annual pan evaporation ranges from 1,974 to 3,183 mm and significantly exceeds the annual precipitation.The annual mean air temperature is 4.77 °C in the center of the basin, while it is 1.53 °C at the boundary of the basin.The Qaidam Basin is a plateau basin and the elevation varies between 2,652 and 6,600 m (Figure 1).The amount of available surface runoff of Qaidam Basin is approximately 49.6 × 10 8 m 3 and the workable reserves of groundwater about 17 × 10 8 m 3 .However, the water resources are very unevenly distributed in space and time around Qaidam Basin.With the development of agriculture and industry, the water consumption increases and results in water shortage in irrigation.Considering the administrative divisions and the integrity of river systems, the Qaidam Basin is made up of eight hydrological sub-regions (Figure 1).The eight hydrological sub-regions and their areas are shown in Table 1.

Surface Energy Balance System Algorithm
The SEBS algorithm, developed by Su [15], can be used to determine turbulent heat fluxes by employing satellite and meteorological data.It consists of: (1) an estimation of a series of land surface physical parameters, such as emissivity, albedo, vegetation coverage etc. based on spectral reflectance and radiance [36]; (2) an extended model of roughness length estimation for heat transfer [25]; (3) an evaporative fraction estimation at limiting cases by energy balance.Three parts of information are required to be input in SEBS model.The first part is the land surface physical parameters, such as emissivity, albedo, temperature, vegetation coverage fraction, leaf area index and roughness height.The second part consists of climate factors of air temperature, air pressure, relative humidity and wind speed.The third part includes measured downward solar radiation and downward longwave radiation.
The basic surface energy balance equation can be expressed as: where R n is the net radiation, G 0 is the soil heat flux, H is the turbulent sensible heat flux, and LE is the turbulent latent heat.Net radiation, R n , can be calculated by: where α is the albedo, R swd and R lwd are the downward shortwave and longwave solar radiation (W•m −2 ), respectively, ε is the surface emissivity, σ is the Stefan-Bolzmann constant (5.67 ×10 −8 W•m −2 •K −4 ), and T 0 is the surface temperature (K).R swd and R lwd are estimated by meteorological measurement.α, ε and T 0 are physical parameters and can be obtained from remote sensing data.The soil heat flux, G 0 , which is energy entering into the soil, can be estimated as follows: where Γ c and Γ s are empirical coefficients.Γ s is 0.315 for most bare soil condition [37] and Γ c is commonly taken as 0.05 for full vegetation [38].The fractional canopy coverage f c would be used as an interpolation between these cases and it can be determined by the scaled Normalized Difference Vegetation Index (NDVI) as: where NDVI max and NDVI min are the NDVI for full vegetation (f c = 1) and bare soil (f c = 0), respectively.
In SEBS algorithm, the energy balance computation at limiting conditions would be used in order to derive the relative evaporation.At the dry limit, the sensible heat flux is at the maximum value, as well as evaporation would be zero because of limitation of soil moisture.The expression follows as: At the wet limit, the sensible heat flux can be taken as the minimum value, and the ET is determined by the available energy in the given surface and the atmospheric conditions.It follows from the equation as: Based on above process, the relative evaporation can be derived as: Finally, the evaporation fraction which is the energy used for the evapotranspiration process is estimated as: The basic formulation of SEBS algorithm is constituted by Equations (1-8) [15].The daily evaporation can be estimated based on the assumption that evaporative fraction is approximately equal to the instantaneous value.
where E daily is the actual evaporation (mm•d −1 ), L is the latent heat of vaporization (J•kg −1 ), is the daily net radiation and is ρ w the water density (kg•m −3 ).
Roughness terms are used to characterize the land surface condition.Roughness height of heat transfer, z oh , can be estimated from a model by Su et al. [25].kB −1 or excess resistance term is taken into account to determine momentum and heat transfer over different ground surfaces in one source model of SEBS [39].kB −1 is derived by: where f c is the fractional canopy coverage and f s is fraction of non-vegetated soil (1-f c ). C d is the drag coefficient of the foliage elements assumed to take the value of 0.2.C t and C t * are coefficients of heat transfer for a leaf and soil, respectively.n ec is the within-canopy wind speed profile extinction coefficient.z om is the roughness length for momentum transfer.The term of kB s −1 is used for bare soils and was formulated by Brutsaert [40].The roughness height, z oh is determined from: In this study, the daily ET values of the SEBS results were averaged for all pixels of the study area in each image and this mean value was considered to be the daily ET of the total Qaidam Basin.There was a total of 1,122,289 pixels in the study area.Because of cloud cover, there were only several cloud-free images available in one month to be used to estimate daily ET and the average ET of these days was represented the mean daily ET of this month.As a result, the monthly ET variables can be obtained based on mean daily ET by multiplying the number of days in a month and the summation yields the annual ET a of SEBS.

MODIS Data
The MODIS satellite data were downloaded from the website of NASA's Earth Observing System Data and Information System.The satellite and solar zenith and azimuth angles were "unpacked" by MODIS Swath re-projection tool (MRTSwath) which is a standard operation.MODIS surface reflectance products (MOD09A1 product) and MODIS land surface temperature data (MOD11A2 product) were derived by monitoring the Earth's surface.MOD9 and MOD11 are based on eight-day composites and have spatial resolution of 500 m and 1,000 m, respectively.1-7 bands of MOD09 and 31, 32 bands of MOD11 were selected.In this study, 3,036 cloud-free MODIS images covered the Qaidam Basin in different months of the 11 years from 2001 to 2011 and were used to estimate the evapotranspiration in regional scale.The numbers of available image composites for total basin were 46 per year after mosaicking.All the images were resampled to 500 m resolution with the Nearest Neighborhood method to match spatial resolution of spectral reflectance bands.
MODIS NDVI data (MOD13A1), which is generated by daily MODIS Level-2G surface reflectance and has been corrected for ozone absorption, aerosols and molecular scattering, can depict spatial and temporal variation of vegetation activities.With 250 m spatial resolution, the MODIS NDVI data has been used in a wide range of researches, such as drought monitoring, global vegetation variation, agricultural and hydrologic modeling [41][42][43][44][45][46][47][48].Since June, July, August and September are the most productive seasons for vegetation growth during a year in Northwest China, 88 MODIS NDVI images of 16-day composites of June, July, August and September from the year 2001 to 2011 were used for illustrating the long-term variation of vegetation in the Qaidam Basin.

GLDAS
GLDAS (the Global Land Data Assimilation System) drives multiple, offline (not coupled to the atmosphere) land surface models, integrates a huge quantity of observation based data, and executes globally at high resolutions (2.5° to 1 km), enabled by the Land Information System (LIS) [49].The temporal resolution of the GLDAS products is three-hourly, and monthly products are also generated through temporal averaging of the three-hourly products.Currently, GLDAS drives four land surface models (LSMs): Mosaic, Noah, the Community Land Model (CLM), and the Variable Infiltration Capacity (VIC).GLDAS bands 19-22 to specific humidity (kg•kg −1 ), wind speed (m•s −1 ), air temperature (K), and air pressure (Pa), respectively, were selected in this study.In total, 132 Sets of Monthly GLDAS products at a spatial resolution of 0.25 degree over the Qaidam Basin from the year 2001 to 2011 were used to estimate evapotranspiration.The GLDAS data was resampled and interpolated to have the same spatial resolution as the ready-processed MODIS data.

Meteorological Observation Data
The corresponding meteorological dataset of air temperature, precipitation, relative humidity, wind speed and pan observed evaporation to the remote sensing data in 8 meteorological stations (Mangya, Lenghu, Xiaozaohuo, Dachaidan, Delingha, Golmud, Nuomuhong, Dulan, Wulan) was used in this study (Figure 1).The monthly average precipitation, air temperature, relatively humidity and wind speed of these eight stations during the period of 2001 to 2011 were available in Qaidam Basin.However, only monthly pan evaporation from 2006 to 2009 data were collected and subsequently measured by E601 dish.The maximum elevation is 3,950 m for Wulan station with a minimum of 2,767 m in Xiaozaohuo station.The yearly net radiation of Golmud station from 2001 to 2011 was used in analysis of its impact on ET.
Conversion coefficient of evaporation was defined as the ratio of actual evapotranspiration to observed pan evaporation in this study.The regional actual evapotranspiration plays a very important role in the hydrological cycle and groundwater modeling.However, due to inconvenient transportation and a lower population in Northwestern China, there is little field measurement data of evaporation except the pan observations in meteorological stations.Therefore, the conversion coefficient of evaporation has a significant effect on estimation of actual evapotranspiration and depends mostly on the empirical formulation in these areas.In this study, the correlation between the SEBS actual ET a and the pan observations was established to find the conversion coefficient of evaporation.

Water Table Depth
Water table depth is a factor that can influence actual evapotranspiration at a site.The relationship between ET a and water table depth reflects the evaporation behavior to groundwater depth.Field measurements of groundwater level in the Qaidam Basin were investigated at WTM and DXS sub-regions during June 2011.Groundwater levels in most of the sites were measured in wells as shown in Figure 2. Consequently, the number of groundwater level measurements in the WTM sub-region is 133 with a number of 107 in the DXS sub-region.The contour maps of water table depth in these two areas were constructed by interpolating field measurements to the same resolution grid as the MODIS data.In most parts of the DXS area, the water table depth is less than 8 m, and the depth in the WTM area is smaller than 3 m.Therefore, 3 m and 10 m can be considered the limit of the water table depth in these two areas, respectively.
For correlation analysis with ET a , the contour map of the water table depth was graded to the same 500 m resolution grid as the MODIS data using the arithmetic mean.The new water table depth map matches the ET a map.For each grid cell, a pair of ET a and water table depth values can be read.There are in total 4,833 and 8,244 pairs of ET a and water table depth in WTM area and DXS area, respectively.
The range of water table depth from 0 to 8 m in DXS was divided into a total of 80 intervals with 0.1 m in each interval.In each interval the corresponding daily actual evapotranspiration values in June 2011 were averaged.The mean values represent the general conditions of the actual evapotranspiration in different water table depth.The same above program was also applied in WTM area.

Spatial and Temporal Distribution of Actual Evapotranspiration by SEBS
The spatial distribution of mean daily ET a by SEBS model in eight sub-regions of Qaidam Basin is illustrated in Figure 3 for 11 years.The mean estimated daily ET a over the total basin varies from 0 to 3.2 mm•d −1 .The areas with higher mean daily ET a are lakes and vegetated areas and the mean ET a in these areas ranges from 1.7 to 2.5 mm•d −1 .The western and northern parts of Qaidam Basin are dominated by desert where the daily ET a varies from 0 to 0.7 mm•d −1 .The eastern and southern parts of Qaidam Basin are mainly covered by shrubs and alpine meadows, which have relatively higher actual evapotranspiration.The SEBS modeled annual ET a in Qaidam Basin and its eight hydrological sub-regions from 2001 to 2011 are shown in Table 2.
In Table 2 it can be noted that the maximum estimated annual ET a by SEBS is 266.8 mm in 2011 over DXS sub-regions, and the lowest ET a estimate is 139.0 mm in WTM areas.The areas with higher estimated ET a in Qaidam Basin are GLM and QDL sub-regions where the east parts are located and the annual ET a of 2011 are 227.4mm, 239.3 mm, respectively.The lower annual ET a areas are MLH and YCD sub-regions in west parts of Qaidam Basin and the annual ET a values are 145.9mm, 165.7 mm, respectively.It can be seen that the annual actual ET a of all eight sub-regions increased with time and approximately achieved its biggest increase from the year 2005.Although the minimum and maximum values of actual ET a were different, the rate of increase was more or less the the same.Compared with the other sub-regions, the actual ET a of MLH, WTM and YCD (the red, cyan and green lines) sub-regions in the western part of the basin were the smallest because of sparse vegetation and very little precipitation.QDL, BDL, DXS and GLM sub-regions (the black, purple, orange and brown lines) in the eastern part of Qaidam Basin presented relative high vegetation cover and thus had higher values of actual ET a .

Conversion Coefficient of Evaporation
The observed pan evaporation of seven meteorological stations (Lenghu, Xiaozaohuo, Dachaidan, Golmud, Nuomuhong, Wulan and Dulan) from 2006 to 2009 are available and are used for estimation of the conversion coefficient of evaporation.Since the annual precipitation increases from the western to the eastern area, which is same as the variation of estimated ET a by SEBS, the Qaidam Basin is divided into two parts for estimation of conversion coefficient between SEBS ET a and pan observations.Lenghu, Dachaidan and Xiaozaohuo stations were selected for west area as Golmud, Nuomuhong, Wulan and Dulan stations were considered as the analysis in eastern regions.
Figure 5a presents the correlation of yearly estimated SEBS ET a with corresponding observed pan evaporation of Lenghu, Dachaidan and Xiaozaohuo stations in the western part of Qaidam Basin.The R-square is 0.69 and the evaporation coefficient is 0.075.It means the SEBS ET a is 7.5% of observed pan evaporation.The correlation of yearly SEBS ET a with corresponding observed pan evaporation of Golmud, Nuomuhong, Dulan and Wulan stations in the eastern area of Qaidam Basin is illustrated in Figure 5b.The evaporation coefficient is 0.141 and the R-square is 0.72.

Effect of Precipitation, NDVI and Net Radiation on Actual Evapotranspiration
The increase of annual actual evapotranspiration can be understood by considering the variation about the annual precipitation, vegetation cover and net radiation.Figure 6 shows the increase of annual precipitation, mean NDVI and net radiation across Qaidam Basin.The annual precipitation ranged from 60 to 133 mm, with the maximum being 154 mm in 2010 (Figure 6a).The big increase of precipitation from 2001 to 2005 resulted in a spike of actual evapotranspiration in 2005.NDVI is an index for monitoring the vegetation cover and it was increased gently after 2002 in Qaidam Basin.It varied between 0.127 and 0.165 with a maximum of 0.172 in 2010 (Figure 6b).The net radiation is the energy driver for ET change and the yearly net radiation ranged from 1,840 to 2,149 W•J•m −2 with the maximum of 2,343 W•J•m −2 in 2010 (Figure 6c).The combined effect of precipitation, vegetation and net radiation resulted in increases of actual evapotranspiration.However, precipitation and vegetation are the dominant constraint factors of actual evapotranspiration in the Qaidam Basin.The ET a is very small in the dry and cool winter season (November, December, January and February) and it is mostly concentrated in spring, summer and autumn seasons.Higher precipitation and vegetation cover results in higher ET a in wet and hot seasons in the basin.

Meteorological Factors and Actual SEBS Evapotranspiration
Monthly precipitation, air temperature, relative humidity and wind speed were averaged, respectively, over eight meteorological stations (Lenghu, Dachaidan, Xiaozaohuo, Golmud, Nuomuhong, Delingha, Wulan and Dulan) in Qaidam Basin. Figure 7 shows the correlation of monthly SEBS ET a with corresponding monthly precipitation, air temperature, relative humidity and wind speed from 2001 to 2011.The result indicated that SEBS ET a has a positive correlation with precipitation and air temperature; the R-square is 0.656 and 0.675, respectively.On the other hand, the relative humidity and wind speed have little correlation with SEBS ET a ; the R-square is 0.453 and 0.076, respectively.Monthly ET (mm) SEBS is highly sensitive to temperature in all input meteorological conditions [39].Su [15] and Badola [50] also reported that SEBS is sensitive to the gradient of surface-air temperature.With the way SEBS is formulated, air temperature can affect the calculated value of sensible heat flux and further has a propagation effect on determination of latent heat flux.

Vegetation and Actual SEBS Evapotranspiration
Generally, the vegetation of Qaidam Basin is sparse and the vegetation cover fraction is very low.Vegetation is dominated by desert plants.As the depth to groundwater variation, low water table plants such as Artemisia desertorum will be changed into high water table plants such as Phragmites australis.Correspondingly, low salt tolerance vegetation such as Artemisia desertorum and Phragmites australis will be replaced by high salt tolerance vegetation such as Tamarix ramossissima Ledeb and Achnatherum splendens with the variation of soil salt content.MODIS NDVI of June, July, August and September data were averaged to represent the vegetation condition of the Qaidam Basin. Figure 8 presents the variation of annual SEBS ET a and mean NDVI during the period of 2001-2011.The result indicated that the annual actual evapotranspiration was positively correlated with NDVI and the R-square is 0.67.

Land Cover Type and Actual SEBS Evapotranspiration
The land cover map used in this research is from a combination of MODIS NDVI of 2011 and field investigation (Figure 9).It is used in the analysis of the actual evapotranspiration of different land cover types in the Qaidam Basin.The variation of actual evapotranspiration across the basin is an index of the water consumption of existing land cover types.The long-term change of evapotranspiration or water use of the different land cover types in Qaidam Basin is analyzed by monthly SEBS ET a and the land cover map.The long-term variation of monthly SEBS ET a for different land cover types in Qaidam Basin are as shown in Figure 10.It can be noted in Figure 10 that the mean monthly actual evapotranspiration of all land cover types were increased with time.The monthly ET a of water bodies, mash, mixture of farmland/trees and shrubs had a relatively substantial increase since 2005.However, the ET a of bare soil had no significant change from 2001 to 2011.The variation of mean monthly ET a for different land cover types can be better understood by increases of precipitation and vegetation cover.A section across different land cover types from north to south in the basin was selected to present and validate the ET a variation (Figure 11).The two figures (Figure 10 and 11) also showed that water bodies have a higher mean actual evapotranspiration than expected and this can be explained by high moisture availability for water bodies.In terms of the relationship between areal coverage of land cover type and monthly SEBS ET a , the small land covered areas have higher mean actual evapotranspiration.Water bodies covering 1.5% of the basin have a higher mean monthly actual evapotranspiration than 110 mm.Mosaic of trees/farmlands and marsh that covered 1.1% and 0.9% of the basin also have relatively high mean actual evapotranspiration rates of 52 mm and 90 mm, respectively.On the contrary, land cover types that covered large parts of the basin have low mean daily actual evapotranspiration.Bare soil has lowest monthly actual evapotranspiration of 15 mm but the land cover type covered 59.5% of Qaidam Basin.Shrubs covered 24.1% of the basin has relatively low mean monthly actual evapotranspiration rates of 34 mm.

Water Table Depth and Actual SEBS Evapotranspiration
The scatter plots ET a against water table depth are presented in Figure 12, and the red line refers to the variation of average daily ET a with water table depth.In Figure 12a, daily ET a values vary between 0 and 3 mm and water table depth range from 0 to 3 m, a dark cloud of data points concentrate in the area of 1 m < water table depth < 2 m.There is declining trend of daily ET a with increasing water table depth.This agrees well with field investigation: in the WTM area where groundwater is shallow, the density of shrubs is larger, and more than 90% of the WTM area is covered by shrubs dominated by P. australis.The root zone of P. australis is shallow (less than 2 m) and the shrub is highly dependent on soil water near the land surface.The average peak value of ET a is 1.5 mm•d −1 and the groundwater depth of 1.2 m corresponded to the P. australis.The average ET a decreases gradually when the groundwater depth is greater than 1.5 m and the mean ET a of this area is 1.3 mm•d −1 .The relationship between ET a and water table depth in DXS area is more complex because this zone is a mixture of farmlands and shrub lands.Two dark clouds of data point concentration occur: one at the water table depth between 0 and 2 m and another between 4 and 7 m.The first may correspond to the farmland.The crops are irrigated and have higher ET a values.The second cloud of data may reflect on shrubs dominated by A. desertorum.The average ET a is in declining trend with increasing groundwater depth.However, the average ET a varies little when the water table depth is larger than 4 m and the value is around 1.5 mm•d −1 which is close to the mean daily ET a of DXS area.Since the recharge of groundwater is small, the ET a mainly comes from soil water in the areas of water table depth larger than 4 m.It is deduced that the threshold depth of groundwater evaporation of the rainy season in the DXS area is about 4 m.The surface ET is roughly equal to soil water evaporation (1.5 mm•d −1 ) when the groundwater depth is larger than 4 m.

Conclusions
The Surface Energy Balance System (SEBS) algorithm was used in this research to estimate spatial actual evapotranspiration in Qaidam Basin and its eight hydrological sub-regions.In general, the actual evapotranspiration of east part was higher than west part of the basin.SEBS actual evapotranspiration estimates were evaluated in two stations and it was on average 7.5% and 14.1% of observed pan evaporation in west and east areas respectively.It was discussed that the difference could be due to variation of precipitation, temperature, roughness and heterogeneity parameters.
In terms of controlling factors, the east area was observed to be high actual evapotranspiration was discussed to be due the high vegetation coverage and precipitation.It was indicated that the land cover types with high evapotranspiration were water bodies, marsh and trees/farmland.These three land cover types covered less than 3.5% of the basin.Land cover types with lower evaporative behavior were sparse shrubs and bare soil, they covered almost 84% of the total study area.The impact of water table depth on actual evapotranspiration is more complex.Generally, the actual evapotranspiration is decreased with increase of water table depth.
Overall, SEBS was proved to be a useful tool for estimation of actual evapotranspiration in basin area.The SEBS actual evapotranspiration can be very helpful for regional hydrological assessment, water resources management and planning in arid and semi-arid areas.However, the more incisive analysis of effect of water table depth on actual evapotranspiration is needed in the near future.

Figure 1 .
Figure 1.Location of Qaidam Basin and eight hydrological sub-regions.The black points represent the meteorological stations and the black solid lines are the boundary of hydrological sub-regions.

Figure 2 .
Figure 2. Distribution of groundwater observation wells in WTM and DXS sub-regions.

Figure 4
Figure 4 presents the increase of SEBS annual ET a in mm during the period 2001-2011 over Qaidam Basin and the eight sub-regions.Overall, the annual SEBS ET a of the basin is increased with time and the value varies from 72.7 to 182.3 mm with the rate of 6.25 mm per year (the bold black line).

Figure 4 .
Figure 4. Variation of annual SEBS ET a of eight hydrological sub-regions and the total basin during the period of 2001-2011.

Figure 6 .
Figure 6.Variation of (a) annual precipitation, (b) NDVI and (c) net radiation during the period of 2001-2011 in Qaidam Basin.

Figure 8 .
Figure 8. Variation of annual SEBS ET a and mean NDVI from 2001 to 2011.

Figure 9 .
Figure 9. Land cover and vegetation types of Qaidam Basin in 2011.The black solid line represents a section from north to south.

Figure 10 .
Figure 10.Monthly SEBS ET a of different land cover types in Qaidam Basin.

Figure 11 .
Figure 11.Section of monthly ET a variation for different land cover types.

Figure 12 .
Figure 12.The scatter plot of daily ET a against water table depth in (a) WTM and (b) DXS in June 2011.

Table 1 .
Areas and abbreviations of hydrological sub-regions in the Qaidam Basin.

Table 2 .
Annual SEBS ET a (mm) of total basin and eight hydrological sub-regions from 2001 to 2011.