Energy-Water Balance and Ecosystem Response To Climate Change in Southwest China
Energy-Water Balance and Ecosystem Response To Climate Change in Southwest China
Energy-Water Balance and Ecosystem Response To Climate Change in Southwest China
Abstract
It is important to highlight energy-water balance and ecosystem response to climate
changes. The change of water-energy balance and ecosystem due to climate change will
affect the regional ecological and human living significantly, especially in Southwest
China which is an ecologically fragile area. This chapter presents the retrieval
methodology of parameters (reconstruction of vegetation index, land cover semiautomatic classification, a time series reconstruction of land surface temperature based
on Kalman filter and precipitation interpolation based on thin plate smoothing splines),
time-series analysis methodology (land cover change, vegetation succession and
drought index) and correlate analysis methodology (correlation coefficient and
principal component analysis). Then, based on the above method, remote sensing data
were integrated, a time series analysis on a 30-year data was used to illustrate the waterenergy balance and ecosystem variability in Southwest China. The result showed that
energy-water balance and ecosystem (ecosystem structures, vegetation and droughts)
have severe response to climate change.
Keywords: energy-water balance, climate change, ecosystem, droughts estimation,
vegetation index, time series analysis, land cover change
1. Methodology
Southwest China consists of the municipality of Chongqing and the four provinces of Guang
xi, Yunnan, Guizhou and Sichuan, 2107N to 3414N, 9730E to 11205E (Figure 1). South
west China, which was considered to be rich in water resources, has suffered climate changes
2016 The Author(s). Licensee InTech. This chapter is distributed under the terms of the Creative Commons
Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution,
and reproduction in any medium, provided the original work is properly cited.
in the last decades. It is particularly important to highlight energy-water balance and ecosys
tem response to climate changes in Southwest China using the remote sensing techniques and
Geographic Information System(GIS). This chapter introduced the methodologies for retriev
ing parameters of land cover, vegetation index, gridded precipitation, land surface tempera
ture, which can be used as the factors of climate change and ecosystem. Furthermore, the climate
change response of droughts estimation methods was presented. The time series analysis of
water-energy balance and ecosystem variability of Southwest China in recent 30 years showed
that ecosystem structures, vegetation and droughts were affected by the change of climate. The
following paragraphs describe the recent achievements in more detail.
d
l=
d
i
-1
i
(1)
EVI = l i EVI i
i =1
(2)
where EVI is EVI value of poor quality pixel after interpolation and EVIi is the EVI value of a
good pixel. The continuous initial EVI time series is recorded as EVI0.
Each pixel in the 163 images in seven years was studied using time series analyzing methods.
Based on the EVI0 data, we firstly generated a time series curve of each pixel and dealt with
noisy signal for each curve using Savitzky-Golay filter, and then extracted a long-term change
trend (EVItr). After that, the weight of each point, which was in the 163 samples time series
( W i ) data, was computed and the EVI time series were traversed.
1.1.2. Ecosystem (land cover): land cover semi-automatic classification from multispectral remote
sensing imagery
Through the conversion of forests and grasslands to croplands and pastures, humans have
affected the exchange of energy, water and carbon between the atmosphere and the land
surface. As the key input to ecosystem researches, various studies focused on land cover
mapping [11]. Whereas, land cover data at large scale is hard to approached except remote
sensing [12]. A variety of approaches was used to map land cover based on remote sensing
data [1317], including visual interpretation classification, unsupervised clustering coupled
with extensive ancillary data and manual labeling of clusters, supervised classification, expert
system classification, artificial intelligence neural network classification, and decision tree
classification. However, the accuracy and the efficiency of land cover classification is not
guaranteed and the land cover classifications are arbitrary using these method. For example,
supervised classification methods require selecting training samples which rely on substantial
expertise and human participation, so the result of land cover classification is influenced
greatly, and it is impossible to classify land cover automatically with these methods. The
algorithms such as neural network classification and fuzzy logic classification are difficult to
understand and apply widely because of highly complicated in their algorithm basis. The
construction of the decision tree and the assignment of thresholds for each sub-nodes are the
key problem of decision tree classification, and they heavily depends on human experience
and varies spatially and temporally. In order to solve these problems to improve the accuracy
and efficiency of classification, Jiang et al. proposed an efficient automatic landscape classifi
cation approach using prior accurate land-cover data as the background experience [18, 19].
This method consists of two steps: (1) semi-automatically detecting land cover changed pixels
from satellite images compared with prior land cover map; (2) semi-automatically classifying
the land cover of changed pixels based on pattern recognition and changed rules.
Automatic collection of training samples: In this method, pure pixels of land cover were
extracted automatically with an accurate previous land cover dataset as prior knowledge. The
interiors of individual land cover areas and larger patches are considered to be more ecolog
ically stabile areas. Based on accumulation area threshold (Pa), the samples of different land
covers sorted in descending order of their area. The accumulation area threshold is calculated
based on the percent of largest patches of specific land cover categories occupied in its total
area.
After spatial buffer analysis, the joint region of different land cover were discarded. The buffer
area are obtained using the different distance to all patches because of patches of vary areas.
The distance of buffer analysis is
buffer
Ab
A
, d<0
(3)
where Pbuffer is area threshold for buffer analysis, Ab d is the buffer area of the patch with a
distance of d, with d negative, and A is the area of the patch.
Establishment of three-dimensional feature space: The data in all spectral bands of each land
cover class, which were extracted from the region of interest, were processed by principal
component statistical analysis. The first three principal components for orthogonal decompo
sition was selected to construct the three-dimensional feature space of different land cover
classes.
j =1
P ji - MP ji
s ji
<c
(4)
where P is the principal component (PC), MP and are mean and standard deviation of PC,
respectively, c is the radius of the three-dimensional feature space.
Change detection and classification of changed pixels: Combining the satellite images and
early land cover maps, the spectral data of the images were extracted according to accurate
early land cover maps. Based on spatial intersect analysis with corresponding three-dimen
sional feature space, the pixels with spectral data outside the ellipsoid were detected as
changed pixels. After obtaining the changed land cover pixels, the satellites images and threedimensional feature space were used to classify the changed area of land cover by calculating
the minimum spectral distance. For each changed pixel, the spectral data of all bands were
input to the formula for all land cover classes in three-dimensional feature space to calculate
the minimum spectral distance dmi .
To express the rules of changed land cover, we proposed a drag coefficient of changed land
cover (r). Combining dmi and r, we determined the final land cover classification of changed
pixels. The minimum distance of the land cover classification based on changed rules is defined
as:
dL
where dL
mi is
mi
= p * d mi
ij
(5)
the minimum distance of the mth pixel to the ith land cover class based on the
changed rules, and rij is the drag coefficient of the ith land cover class that changed to the jth
land cover class.
The class information would be contaminated by adjacent class codes, so the post-classification
results was modified to solve the problem. Any group of pixels, which was <3 pixels in size,
was identified as noise in dilation operation.
1.1.3. Energy parameter: land surface temperature
The parameter of temperature of surface energy balances has be important index for studies
of ecosystem response to climate change. As land cover mapping, the land surface temperature
retrieval based on remote sensing has been the primary method for gaining temperature at
large scale. This is due to the accessibility and convenience of satellite remote sensed infor
mation. However, the surface energy parameters retrieved from remote-sensing data are often
interrupted on the spatial and temporal scales by clouds, aerosols, solar elevation angle and
bidirectional reflection. so, it is always difficult to obtain complete dataset for a large region.
In addition, the accuracy of surface parameter retrieval was affected with various degrees by
an indirect retrieval method and the instantaneous features of monitoring [2, 20]. To reduce
such impacts, a time-composite method is generally adopted. For surface energy parameters,
the time series fitting and noise-removal methods commonly used to include the mean diurnal
variation and nonlinear regression methods. Compared with nonlinear regression method, insitu measurement data are necessary for the mean diurnal variation method [10, 21]. However,
the modeled results often cannot represent the actual situation under significant changes of
environmental conditions [20, 22]. In recent years, data-assimilation methods have been
adopted in the reconstruction of time series data to solve these problems. In this paper, a time
series reconstruction method based on an ensemble Kalman filter was established, which is
can be applied in various studies[19,2326]. It focused on evaluating the state of a discretetime controlled process. The state Xk can be expressed as [27]:
X =f
k
k , k -1
k -1
+ G k , k -1W k -1
(6)
k ,k 1and k ,k 1are
the noise
T
E V k = 0, E V kV j = R kd kj
T
E
=0
WV j
(7)
where Qk is the nonnegative covariance matrix of W k and Rk is the positive covariance matrix
of V k ; kj is the function of Kronecker-.
= 0.5 * cov(randn(1, N ))
(8)
= 0.5 * cov(randn(1, N ))
(9)
where X k 1 is defined to be our a priori state estimate at step k, given knowledge of the process
X
and the a posteriori state as follows:
k , k -1
=f
k , k -1
k -1
(10)
X =X
k
k , k -1
+ K k Z k - H k X k , k - 1
(11)
The calculation process of a Kalman filter is a constant forecast correction process. Based on
time-update and observation-update, the values were reconstructed with minimum variance
by compared to in situ values.
1.1.4. Water balance: precipitation interpolation based on thin plate smoothing splines
Precipitation data are one of the important input elements of ecological mechanism model,
and it play an extremely important role in simulating and researching ecosystem changes at
regional or global scale. It also a key index in global change researches, which is a direct
parameter for drought estimation [28]. The precision of precipitation data may affect the
simulation and prediction precision. At present, the precipitation data acquisition mainly rely
on the long-term observation of the weather stations, but scarcity and uneven distribution of
weather stations reduce the accurate of spatial precipitation data.
The precipitation measurements from the WMO stations are point measurements. Spatial
interpolation technique for meteorological elements is often used to obtain the meteorological
data of every position in the scope which is covered by the meteorological stations. Several
methods have been developed to interpolate these point data to a real estimation of rainfall.
The software ANUSPLIN, developed by Hutchinson, Australian National University, was
applied to generate gridded precipitation data [29]. ANUSPLIN is a kind of analysis and
interpolation tool for multivariable data using local thin plate spline function. The model
required to enter the location of the meteorological site, elevation and other ancillary data. Its
statistical analyzes and diagnose multivariable data to realize spatial interpolation function.
In order to fit to datasets distributed across an unlimited number of climate station locations,
interpolating methods of ANUSPLIN were applied, which use thin plate smoothing splines
for spatial interpolation with a third parameter of elevation [3036]. The main advantage of
thin plate splines is that splines do not require prior estimation of spatial auto-covariance
structure, which was difficult to estimate and validate. The partial spline observational model
for n data values zi at positions x is given by setting [37]:
= f
( x ) b j ( x ) + e (i = 1,..., n; j = 1,...p)
i
j =1
(12)
E
where = (,
(ee ) = V s
T
(13)
the errors i are uncorrelated if V is diagonal and correlated otherwise. The function f and the
parameters j are estimated by minimizing
( z- g ) V
T
where z = (z,
-1
( z - g ) + pJ ( f )
(14)
g = g(x ) = f (x ) + b j (x )
i
where
j =1
(15)
10
n
S =
i =1
( DS
S ) 100
i- j
t 100%
(16)
where S is the comprehensive land cover dynamic degree in t period, S i j is the total area of
land cover change of ith type converted into jth type from monitor beginning to end. Si is the
area of ith type when the monitor started.
K = (S
i
it 2
- S it 1
) S (t t ) 100%
it 1
2 1
(17)
K i is the single land cover dynamic degree of ith type from t1 to t2 time, Sit 2and Sit 1
were areas of ith type in t2and t1time separately.
where
S
=
Sij ...21
S n1
S
S
12
22
...
n2
S
S
...
2n
... ...
... S nn
...
1n
(18)
(19)
where x represents the grid cell and t represents the period in which NPP is accumulated, for
example, a month. APAR is determined by the fraction of photo synthetically active radiation
(FPAR) and the total solar surface radiation (SOL) (MJ m2) [57] as
11
12
(20)
where the constant 0.5 represents the ratio of the total solar radiation (with a wavelength range
of 0.40.7 m) used by the vegetation [58].
LUE is calculated as the product of maximum light use efficiency, and its temperature and
moisture stressors [56] as
LUE ( x , t ) = T e 1 ( x , t ) T e 2 ( x , t ) W e ( x , t ) e max
(21)
where LUE (x, t) represents the actual light use efficiency, maxthe maximum light use efficien
cy, and the value for grass (0.604 g/MJ), simulated by Running based on BIOME-BGC model
[59], was used here; T 1(x, t)and T 2(x, t)are temperature scalars, and W (x, t) is the moisture
stress coefficient. T 1(x, t), T 2(x, t) and W (x, t)were computed at every location at each time
step. T 1(x, t)and T 2(x, t)are calculated as [56, 57].
T e ( x,t ) =
2
where
{1 + exp 0.2 (T
T opt is
opt
opt
( x)
2
1.1814
opt
( x ) - 10 - T ( x , t )}
{1 + exp 0.3 ( -T
opt
(22)
( x ) - 10 - T ( x , t ) )
(23)
maximum normal differential vegetation index (NDVI). T is the monthly mean temperature;
T 2( x ) = 1, when T = Topt ; it decreases to 0.5 when T is 10C above or 13C below T opt .
W (x, t) reflects the effect of water condition, and it generally increases when available water
increases. Atmospheric vapor pressure deficit reflects air humidity, which affects transpiration
and then the LUE [60]. Therefore, there are currently studies using vapor pressure deficit (D
in kPa) to calculate the moisture stress coefficient [61, 62], computed as [49]
W e ( x , t ) = (1.2e
( -0.35 D )
) - 0.2,
- 273
- 273
D = 0.611 exp (17.27 T s
) - exp (17.27 T d
)
36
Ts
T d - 36
(24)
where
T d is
Ts
T s T d < 0, D = 0. Td was derived from Guo Jies regression model for Sichuan province based
on Yang Jingmeis findings of a significant linear relationship between dew point temperature
and the logarithm of total perceptible water [63, 64] as follows:
Ln(U ) = 1.8084 + 0.0735Td
(25)
13
14
PAI = ( P - P ) / P 100%
(26)
PAI (%)
Month
Season
Year
None
[40, )
[25, )
[15, )
Mild
[60, 40)
[50, 25)
[30, 15)
Moderate
[80, 60)
[70, 50)
[40, 30)
Severe
[95, 80)
[80, 70)
[45, 40)
Extreme
( , 95)
( , 80)
( , 45)
Because of difficulty to evaluate absolute droughts, the meteorological drought level was used
to estimate the drought severity. Table 1 shows the relationship between the PAI and mete
orological drought level in three temporal scales of month season and year [80].
The drought frequency (DF) were defined as:
DF = n N 100%
(27)
(X - X)(Y - Y)
(X - X) (Y - Y)
n
xy
i=1
i=1
i=1
(28)
where r xy is the correlation coefficient of variable X and variable Y; n is the number of sample;
is the mean of variable X; Y
is the mean of variable Y.
X
The significance test of correlation coefficient generally uses the t-test. The statistic calculation
formula is:
xy
r n- 2
1- r
2
(29)
15
16
Figure 4 shows the change of mean temperature of Sichuan Province in four seasons from 1982
to 2010. It turns out the spatial distribution pattern changes of mean temperature at different
seasons are not significant. The changes of mean temperature in summer ranged from 6.59
to 22.06C, the change ranged from 6.59 to 22.06C in spring, from 6.59 to 22.06C in autumn,
and from 6.59 to 22.06C in winter. The mean temperature variations in winner is biggest
(30.746C), and it in autumn is smallest (25.27C).
Figure 5 shows the monthly mean temperature changes of time series from 1982 to 2010. As
can be seen from the Figure 5, the monthly mean temperature display cyclic variations under
this time series. It offered upgrade firstly than descending latter tendency. In June, July, August
and September, the temperature was higher and reaches a maximum in July and August. In
January, the monthly mean temperature is lower. The changes of every year were basically
similar, and they have certain periodicity with an obvious fluctuation. The fluctuation of the
highest temperature of each year is very small. A maximum of the highest temperature appears
in the June 2006 was 23.65CA minimum of the highest temperature has an obvious fluctuation
ranged from 1.23(1984) to 4.28(1987).
Figure 6 shows annual changes of the mean temperature from 1982 to 2010. The mean
temperature fluctuated between 12.48 and 13.90C. The minimum of the mean temperature
17
18
occurred in 1992 and the maximum occurred in 2006. From 1982 to 2010, the annual mean
temperature was rising obviously which was consistent with global warming trends.
5.1.2. Precipitation
Taking Sichuan province as an example, we selected the monthly mean precipitation data of
nine weather stations from 1982 to 2010. As can be seen from Figure 7, the spatial distribution
of annual average precipitation was uneven in 29 years with descending from the east to the
west. Taking Songpan Country-Li Country-Kangding City-Jiulong Country-Yanyuan Country
as a boundary, the annual precipitation of the east of the boundary was abundant, and it is
more than 900 mm/a. Especially, the annual mean precipitation of Yaan is more than 1200
mm/a and ranked first in the whole province. The annual mean precipitation of most areas in
the west is not more than 800 mm. Among them, there is a very little rain in Shiqu country,
less than 600 m/a.
As shown in Figure 8, the annual average precipitation was concentrated mostly in summer,
and there are obvious differences in the spatial distributions in different seasons. The spring
precipitation mainly concentrated in the central and eastern regions in Sichuan, and the
maximum precipitation was 311.55 mm. The precipitation in summer moves westward
concentrating in the central regions with the maximum precipitation of 859.68 mm. The scope
of autumn rainfall mainly concentrated in the central and southern Sichuan with uneven
spatial distribution. The maximum precipitation was 340.02 mm. The precipitation of most
regions in winter was low, and decreased from the west to the east.
19
20
Monthly mean precipitation changing trends in every years were basically the same, namely,
it offered upgrade firstly than descending latter tendency (Figure 9). In June, July, August and
September, the precipitation was higher and reaches a maximum in July and August. In
January, the monthly mean precipitation is lower. The precipitation change has certain
periodicity with an obvious fluctuation from 1982 to 2010. The maximum of annual highest
precipitation was 294.49 mm in 1984, and the minimum is 154.03 mm in 2006. The annual
lowest precipitation has no obvious fluctuation ranged from 2.86 to 15.71 mm. The monthly
mean precipitation had not a visible downtrend by 0.11 mm per decade.
Figure 10 shows the changes of the annual precipitation from 1982 to 2010. The annual
precipitation showed obvious decreasing trend. The maximum of the annual precipitation
occurred in 1989 (1131.17 mm), and the minimum occurred in 2006 (812.84 mm). The precip
itation changed in a wide range and had obvious fluctuations.
across the other region. The grassland, including low coverage grassland, medium coverage
grassland and high coverage grassland, were distributed mainly in Sichuan and Guizhou. The
woodland were mainly distributed at Yunnan and Guangxi.
The time series analysis on land cover indicated the whole change is not obvious with certain
changes on land use types during 20 years. High coverage grassland, woodland, urban land
and construction land increase, others decrease.
Figure 11. Distribution of land cover in 1980, 1990, 1995, 2000 and 2005.
21
22
the boundary of Sichuan and Chongqing Province. In general, NDVI in Southwest China
decreased from 2004 to 2010, especially in Sichuan Province.
Figure 12. Distribution of NDVI in 1980, 1990, 1995, 2000 and 2005.
SC suffered from droughts in a large area in January, February, March, October, November
and December. As shown in Figure 13, Yunnan Province and Guangxi Province are severely
drought-afflicted areas. A drought pattern also appeared in the monthly variability of DF over
time. From January to March, central and eastern Yunnan Province, southwestern Sichuan
Province, and southern Guangxi Province experienced droughts (DF >40%) over large area. In
May, the drought area rapidly narrowed, and the area of DF >40% was located in a limited
extension of Yunnan Province. From June to September, the DF was low with a slightly
increase. The distribution of droughts was spread from east to west. In October, droughts
began to occur again in a large area, and the DF of eastern Guangxi Province and the DF of
northwestern Yunnan Province exceeded 40 and 30%, respectively. From November, the
drought area increased rapidly. Half of study area is at a high drought risk. The DF of the entire
Yunnan Province in western SC is >40%. In December, the DF of most of the area of Yunnan
Province exceeded 50%, except in the limited extension of the central portion, which meaning
that these regions suffered droughts biyearly.
23
24
Figure 14. Spatial distribution of annual DF in Southwest China from 1961 to 2012.
To illustrate the spatial distribution of drought variability, the annual DF of SC from 1961 to
2012 was calculated. As shown in Figure 14, the annual droughts DF were scattered, in contrast
to the monthly scale. More than 61% of the area in SC had a relatively high DF (>15%), that is,
all of SC has been a region of high drought risk for more than half century. The southern and
eastern mountain zones have high drought frequency. In contrast, Sichuan Basin, which
occupies a large portion of the study area with relatively flat and fertile grounds, suffered
fewer droughts events.
5.3.2. Ecosystem response: droughts analysis based on vegetation index
The vegetation index data were 1 km, 16 days composited MODIS EVI (MOD13A12), which
were downloaded from the NASA EOS Data Gateway (EDG) (http://modis.gsfc.nasa.gov/
index.php). To be consistent with monthly GRACE estimates, the days EVI product was
recomposited using the maximum value composite (MVC). The composited monthly EVI data
were then linearly interpolated to a 1 1 grid.
Droughts are common in Southwest China, and several episodes of potential severe droughts
were detected using temporal variability analysis of the GRACE TWS change. To evaluate
drought events in Southwest China from 2003 to 2013, three indicators were taken into account:
the PAI, the AVI and the annual cycle removed TWS (TWSI). To compute the TWSI, monthly
averaged GRACE TWS change data for 10 years were removed at each pixel. Using the
spatiotemporal variability analysis described above for the TWS, the precipitation, and the
EVI, we chose representative pixels from three regions for drought analysis. Figure 15 shows
the time series of the three droughts indicators between 2006 and 2012 when droughts were
likely to occur in Southwest China.
Figure 15. Time series of the TWSI, the PAI, and the AVI for selected pixels.
The correlation coefficients between EVI and precipitation were all more than 0.84 from 2003
to 2013, whereas they were 0.05, 0.07, 0.08 and 0.16 between PAI and AVI, respectively. The
low correlation coefficients between PAI and AVI imply that these two drought indicators
predict different water resources deficit conditions that accompany droughts. However,
correlation analysis showed that TWSI present low correlation with EVI, precipitation and
corresponding droughts indicators. Figure 15 shows that the occurrence and release of drought
AVI lagged behind PAI for 13 months and droughts of AVI were more severe than PAI. This
is because of the delay in recharge of surface and soil water from rainfall and the vegetation
growth. Both indicators mainly reflected water depletion in surface water and shallow soil
water. PAI and AVI were both invalid for detecting droughts. This was most obvious during
the period from the end of 2011 to the beginning of 2012 in Figure 15, when the time series for
PAI and AVI had no extreme droughts and a normal and gentle amplitude fluctuation.
However, TWSI showed that Southwest China experienced great water resource decreases
during this period, which may have caused an extreme drought in 2006. Considering the
annual cycle of precipitation and vegetation growth and their relation to shallow water, the
TWS change mainly contributed to the discharge of groundwater to surface water, which
implied a drought risk.
25
26
6. Conclusion
Southwest China is an ecologically fragile area with more than 242 million populations. The
change of water-energy balance and ecosystem due to climate change will affect the regional
ecological and human living significantly. Our study of Southwest China in past 30 years
shows that ecosystem structures destroyed by shrinking of water ecosystem, forest ecosystem
and grass ecosystem, together with the precipitation reduction and temperature rise. In
addition, the climate changes also affected the artificial ecosystems such as crops land, which
resulting in food production mainly caused by frequently occurred severe droughts. Our
research showed that, even in the region with abundant water resources of Southwest China,
energy-water balance and ecosystem have severe response to climate change, which is
significant to human productions and activities of daily livings.
Author details
Jianhua Wang1*, Yaohuan Hang2, Dong Jiang2 and Xiaoyang Song3
*Address all correspondence to: jiangd@igsnrr.ac.cn
1 State Key Laboratory of Simulation and Regulation of Water Cycle in River Basin, China
Institute of Water Resources and Hydropower Research, Beijing, China
2 Institute of Geographical Sciences and Natural Resources Research, Chinese Academy of
Sciences, Beijing, China
3 China University of Mining and Technology, Beijing, China
References
[1] Stephenson N L. Climatic control of vegetation distributionthe role of the waterbalance. American Naturalist, 1990, 135(5): 649670.
[2] Francis A P, Currie D J. A globally consistent richness-climate relationship for angio
sperms. American Naturalist, 2003, 161(4): 523536.
[3] Zhao Y-S, Chen D-M, Yang L-M. The application principle and method of remote
sensing. Science Press, Beijing, 2003.
[4] Tucker C J, Townshend J R G, Goff T E. African land-cover classification using satellite
data. Science, 1985, 227(4685): 369375.
27
28
[18] Jiang D, Huang Y H, Zhuang D F, et al. A simple semi-automatic approach for land
cover classification from multispectral remote sensing imagery. Plos One, 2012,
7(9):e45889e45889.
[19] Fu J Y, Jiang D, Huang Y H, et al. A Kalman filter-based method for reconstructing
GMS-5 global solar radiation by introduction of in situ data. Energies, 2013, 6(6): 2804
2818.
[20] Neteler M. Estimating daily land surface temperatures in mountainous environments
by reconstructed MODIS LST data. Remote Sensing, 2010, 2(1): 333351.
[21] Huang Y H, Wang J H, J D, et al. Reconstruction of MODIS-EVI time-series data with
S-G filter. Geomatics and Information Science of Wuhan University, 2009, 34(12): 1440
1443.
[22] Ziwei X U, Liu S, Tongren X U, et al. Comparison of the gap filling methods of evapo
transpiration measured by eddy covariance system. Advances in Earth Science, 2009,
24(4): 372382.
[23] Aravequia J A, Szunyogh I, Fertig E J, et al. Evaluation of a strategy for the assimilation
of satellite radiance observations with the local ensemble transform Kalman filter.
Monthly Weather Review, 2011, 139(6): 19321951.
[24] Huang C L, Li X, Lu L. Retrieving soil temperature profile by assimilating MODIS LST
products with ensemble Kalman filter. Remote Sensing of Environment, 2008, 112(4):
13201336.
[25] Rosema A. Using meteosat for operational evapotranspiration and biomass monitoring
in the Sahel region. Remote Sensing of Environment, 1993, 46(1): 2744.
[26] Sakov P, Evensen G, Bertino L. Asynchronous data assimilation with the EnKF. Tellus
Series a-Dynamic Meteorology and Oceanography, 2010, 62(1): 2429.
[27] Welch G, Bishop G. An introduction to the Kalman filter. University of North Carolina
at Chapel Hill, Springer, Heidelberg, Berlin, 1995, vol. 785: 127132. ISBN
978-3-642-02287-6.
[28] Wang J H, Jiang D, Huang Y H, et al. Drought analysis of the Haihe River Basin based
on GRACE terrestrial water storage. Scientific World Journal, 2014, 2014:578372578372.
[29] Huang Y, Xu C, Yang H, et al. Temporal and spatial variability of droughts in Southwest
China from 1961 to 2012. Sustainability, 2015, 7(10): 1359713609.
[30] Fleming M D, Chapin F S, Cramer W, et al. Geographic patterns and dynamics of
Alaskan climate interpolated from a sparse station record. Global Change Biology,
2000, 6(S1): 4958.
[31] He B, L A, Wu J, et al. Drought hazard assessment and spatial characteristics analysis
in China. Journal of Geographical Sciences, 2011, 21(2): 235249.
[32] Hong X, Guo S, Xiong L, et al. Spatial and temporal analysis of drought using entropybased standardized precipitation index: a case study in Poyang Lake basin, China.
Theoretical and Applied Climatology, 2015, 122(34): 543556.
[33] Hutchinson M F. Interpolation of rainfall data with thin plate smoothing splinesPart
I: two dimensional smoothing of data with short range correlation. Journal of Geo
graphic Information and Decision Analysis, 1998, 2(2): 153167.
[34] Liu Z, Xingxiu Y U, Wang S, et al. Comparative analysis of three covariates methods
in thin-plate smoothing splines for interpolating precipitation. Progress in Geography,
2012, 31(1): 5662.
[35] Price D T, Mckenney D W, Nalder I A, et al. A comparison of two statistical methods
for spatial interpolation of Canadian monthly mean climate data. Agricultural and
Forest Meteorology, 2000, 101(s 23): 8194.
[36] Qian Y L, Lv H Q, Zhang Y H. Application and assessment of spatial interpolation
method on daily meteorological elements based on ANUSPLIN software. Journal of
Meteorology and Environment, 2010, 26(2):714.
[37] Hutchinson M F. Interpolation of rainfall data with thin plate smoothing splinesPart
I: two dimensional smoothing of data with short range correlation. Journal of Geo
graphic Information and Decision Analysis, 1998, 2(2): 139151.
[38] Jiyuan L. Study on spatial-temporal feature of modern land use change in China: using
remote sensing techniques. Quaternary Sciences, 2000, 20(3): 11.
[39] Crabtree R, Potter C, Mullen R, et al. A modeling and spatio-temporal analysis
framework for monitoring environmental change using NPP as an ecosystem indicator.
Remote Sensing of Environment, 2009, 113(7): 14861496.
[40] Ramakrishna N, Michael W, Peter T, et al. Recent trends in hydrologic balance have
enhanced the carbon sink in the United States. Geophysical Research Letters, 2002,
29(10): 106-1106-4.
[41] Lieth H. Modeling the primary productivity of the world. Springer, Heidelberg, Berlin,
1975: 237263.
[42] Melillo J M, Mcguire A D, Kicklighter D W, et al. Global climate change and terrestrial
net primary production. Nature International Weekly Journal of Science, 1993,
363(6426): 234240.
[43] Golubyatnikov L L, Denisenko E A. Modeling the values of net primary production for
the zonal vegetation of European Russia. Izvestiia Akademii Nauk Seriia Biologiche
skaia, 2001, 28(3): 353361.
[44] Wang P, Xie D, Zhou Y, et al. Estimation of net primary productivity using a processbased model in Gansu Province, Northwest China. Environmental Earth Sciences, 2014,
71(2): 647658.
29
30
[45] Bala G, Joshi J, Chaturvedi R K, et al. Trends and variability of AVHRR-derived NPP
in India. Remote Sensing, 2013, 5(2): 810829.
[46] Cramer W, Kicklighter D W, Bondeau A, et al. Comparing global models of terrestrial
net primary productivity (NPP): overview and key results. Global Change Biology,
1999, 5(S1): 115.
[47] Kale M P, Roy P S. Net primary productivity estimation and its relationship with tree
diversity for tropical dry deciduous forests of central India. Biodiversity and Conser
vation, 2012, 21(5): 11991214.
[48] Potter C S, Rerson J T, Field C B, et al. Terrestrial ecosystem production: a process model
based on global satellite and surface data. Global Biogeochemical Cycles, 1993, 7(4):
811841.
[49] Prince S D, Goward S N. Global primary production: a remote sensing approach.
Journal of Biogeography, 1995, 22(4/5): 815835.
[50] Cao M, Prince S D, Small J, et al. Remotely sensed interannual variations and trends in
terrestrial net primary productivity 19812000. Ecosystems, 2004, 7(3): 233242.
[51] Nayak R K, Patel N R, Dadhwal V K. Estimation and analysis of terrestrial net primary
productivity over India by remote-sensing-driven terrestrial biosphere model. Envi
ronmental Monitoring and Assessment, 2010, 170(14): 195213.
[52] Potter C, Klooster S, Myneni R, et al. Continental-scale comparisons of terrestrial carbon
sinks estimated from satellite data and ecosystem modeling 19821998. Global and
Planetary Change, 2003, 39(s 34): 201213.
[53] Zhang Y, Wei Q, Zhou C, et al. Spatial and temporal variability in the net primary
production of alpine grassland on the Tibetan Plateau since 1982. Journal of Geograph
ical Sciences, 2014, 24(2): 269287.
[54] Bian J, Li A, Deng W. Estimation and analysis of net primary Productivity of Ruoergai
Wetland in China for the recent 10 years based on remote sensing. Procedia Environ
mental Sciences, 2010, 2(6): 288301.
[55] Mccallum I, Wagner W, Schmullius C, et al. Satellite-based terrestrial production
efficiency modeling. Carbon Balance and Management, 2009, 4(9): 114.
[56] Field C B, Randerson J T, Malmstrm C M. Global net primary production: combining
ecology and remote sensing. Remote Sensing of Environment, 1995, 51(1): 7488.
[57] Piao S L, Fang J Y, Guo Q H. Application of casa model to the estimation of Chinese
terrestrial net primary productivity. Acta Phytoecologica Sinica, 2001, 25:603608
[58] Mccree K J. Photosynthetically active radiation. Springer, Heidelberg, Berlin, 1981: 41
55.
[59] Running S W, Thornton P E, Nemani R, et al. Global terrestrial gross and net primary
productivity from the earth observing system. Springer, New York, 2000: 4457.
[60] Jiang L P, Qin Z H, Xie W, et al. A research of net primary productivity model of
grassland based on MODIS data. Chinese Journal of Grassland, 2006, 28(6): 7276.
[61] Chen L J, Liu G H, Feng X F. Estimation of net primary productivity of terrestrial
vegetation in China by remote sensing. Acta Botanica Sinica, 2001, 43(11): 11911198.
[62] Chen L J, Liu G H, Hui Guo L I. Estimating net primary productivity of terrestrial
vegetation in China using remote sensing. Journal of Remote Sensing, 2002, 6(2): 129
135.
[63] Jie G, Ping L G. Climatic characteristics of precipitable water vapor and relations to
surface water vapor column in Sichuan and Chongqing region. Journal of Natural
Resources, 2009, 24(2): 344350.
[64] Yang J, Qiu J. A method for estimating precipitable water and effective water vapor
content from ground humidity parameters. Chinese Journal of Atmospheric Sciences,
2002, 26(1): 922.
[65] Keyantash J, Dracup J A. The quantification of drought: an evaluation of drought
indices. Bulletin of the American Meteorological Society, 2002, 83(8): 11671180.
[66] Lu X, Zhuang Q. Areal changes of land ecosystems in the Alaskan Yukon River Basin
from 1984 to 2008. Environmental Research Letters, 2011, 6(3): 329346.
[67] Raziei T, Saghafian B, Paulo A A, et al. Spatial patterns and temporal variability of
drought in Western Iran. Water Resources Management, 2009, 23(3): 439455.
[68] Glantz D A W, H. M. Understanding: the drought phenomenon: the role of definitions.
Water International, 1985, 10(3): 111120.
[69] Wang J, Jiang D, Huang Y, et al. Drought analysis of the Haihe River Basin based on
GRACE terrestrial water storage. Advances in Meteorology, 2014, 2014: 578372578372.
[70] Lu X, Zhuang Q. Evaluating evapotranspiration and water-use efficiency of terrestrial
ecosystems in the conterminous United States using MODIS and AmeriFlux data.
Remote Sensing of Environment, 2010, 114(9): 19241939.
[71] Huang Y, Wang J, Dong J, et al. Surface water deficiency zoning of China based on
surface water deficit index (SWDI). Water Resources, 2014, 41(4): 481487.
[72] Vicente-Serrano S M. Differences in spatial patterns of drought on different time scales:
an analysis of the Iberian Peninsula. Water Resources Management, 2006, 20(1): 3760.
[73] Li S, Xiong L, Dong L, et al. Effects of the Three Gorges Reservoir on the hydrological
droughts at the downstream Yichang station during 20032011. Hydrological Process
es, 2013, 27(26): 39813993.
31
32