Analysis of long-term terrestrial water storage variations in the Yangtze River basin

In this study, we analyze 32 yr of terrestrial water storage (TWS) data obtained from the Interim Reanalysis Data (ERA-Interim) and Noah model from the Global Land Data Assimilation System (GLDAS-Noah) for the period 1979 to 2010. The accuracy of these datasets is validated using 26 yr (1979–2004) of runoff data from the Yichang gauging station and comparing them with 32 yr of independent precipitation data obtained from the Global Precipitation Climatology Centre Full Data Reanalysis Version 6 (GPCC) and NOAA’s PRECipitation REConstruction over Land (PREC/L). Spatial and temporal analysis of the TWS data shows that TWS in the Yangtze River basin has decreased significantly since the year 1998. The driest period in the basin occurred between 2005 and 2010, and particularly in the middle and lower Yangtze reaches. The TWS figures changed abruptly to persistently high negative anomalies in the middle and lower Yangtze reaches in 2004. The year 2006 is identified as major inflection point, at which the system starts exhibiting a persistent decrease in TWS. Comparing these TWS trends with independent precipitation datasets shows that the recent decrease in TWS can be attributed mainly to a decrease in the amount of precipitation. Our findings are based on observations and modeling datasets and confirm previous results based on gauging station datasets.


Introduction
Terrestrial water storage (TWS) is determined by all physical phases of water stored above and below the surface of the Earth, including soil moisture, snow and ice, canopy water storage, groundwater, etc.As a key component of terrestrial and global hydrological cycles, TWS strongly influences water, energy, and biogeochemical fluxes, thereby playing a major role in the Earth's climate system (Famiglietti, 2004).TWS is not only an indicator of the Earth's climate variability, but also affects various components of the Earth's hydrological cycle (Niu et al., 2006).Soil moisture plays a key role in both the water and energy cycles through its impact on the energy partitioning at the surface, and soil moisture also has links with the biogeochemical cycle via plant transpiration and photosynthesis (Seneviratne et al., 2010).Snow cover has a strong influence on the onset of the summer monsoon and runoff production in spring (Ding et al., 2009).Therefore, the spatial and temporal variability in TWS due to climate change and human-induced impacts both form important components in the water and energy cycles, and should be taken into account in river basin management.
From a historical perspective, there is limited information about the TWS distribution in time and space, as TWS is not routinely assessed like other hydrometeorological measurements.Isolated datasets are available for only a few regions and rarely for periods of more than several years.Moreover, the in situ observations are point measurements, and not always representative for larger spatial domains (Famiglietti et al., 2008;van der Velde et al., 2008).Fortunately, progress in satellite remote sensing and corresponding retrieval techniques enables large scale monitoring of land surface bio-geophysical properties (e.g.soil moisture and temperature).This may potentially improve our understanding of the spatially heterogeneous hydrometeorological processes.Advances in microwave remote sensing have demonstrated their use in providing large-scale soil moisture information, resulting in satellite missions specifically dedicated to soil moisture (Entekhabi et al., 2010).Microwave observations can, however, only provide information on the top few centimeters of the soil.In addition, Tapley et al. (2004a,b) and others have shown that, using measurements of the Earth's gravity field, terrestrial water storage change (TWSC) may be inferred on a monthly scale.The first space mission that employs this technology is the Gravity Recovery and Climate Experiment (GRACE) launched on 17 March 2002.
Data assimilation products such as Interim Reanalysis Data (ERA-Interim) and Global Land Data Assimilation System (GLDAS) combine the virtues of in situ data, remotely sensed observations, and modeling.The models in these systems simulate the main components of TWS and, by fusing these components with other data sources, reduce uncertainties in the hydrological interpretations.These systems have been extensively applied in TWS and related studies, and have, for example, been utilized in regional, continental, and global TWS variation analysis (Chen et al., 2005;Seneviratne et al., 2004;Syed et al., 2008).As well, these systems offer long-term records of data, making them suitable for long-term analysis, while remotely sensed data and in situ observations are most likely time limited.
In this study, we focus on the analysis of the long-term variation in TWS of the Yangtze River basin.The Yangtze River, the longest river in China, forms one of the world's top ten rivers basins as far as water shortage is concerned.This shortage is caused by intensive human water use, and is despite the large volume of runoff the river basin receives (Wong et al., 2007).During the past three decades, the Yangtze River basin has experienced fundamental changes, e.g. a marked increase in temperature, population growth, economic development, water consumption, as well as the dam construction.The Three Gorges Dam (TGD) is the largest hydroelectric dam, and has created the largest manmade lake (more than 600 km 2 of former land) in the world.Such sizeable land use changes alter many factors, such as albedo, regional climate, and the hydrological cycle.In recent years, the basin has experienced an increasing trend in the frequency of extreme events, i.e. low runoff in drought years, and floods during intense rainfall (Dai et al., 2008;IPCC, 2001).A better understanding of the changes occurring in the Yangtze River basin and its hydrological state variables is thus important.However, previous work has mostly focused on the interaction between runoff, precipitation, and evapotranspiration in the basin, while little attention has been paid to space-time variability in TWS and its response to climate change and human activity.
In this paper, we examine the spatial and temporal variation in TWS in the Yangtze River basin, with the aim to improve our understanding of the water cycle and aid management of the water resources.The specific objectives of this paper are (1) to use an ERA-Interim dataset to estimate the TWS and TWSC for the period 1979 to 2010, (2) to assess the accuracy of the TWS estimation from ERA-Interim and GLDAS datasets in the basin, (3) to examine the climatology of the spatial pattern of TWS in the basin, and (4) to detect trends and abrupt changes, as well as their possible causes.

Study area
The Yangtze River basin is located in the subtropical zone in China.The river originates in the Qinghai-Tibetan Plateau and flows 6300 km eastwards to the sea.The upper Yangtze reaches, the headwaters, extend from the westernmost point, at Tuotuohe, to Yichang.The middle reaches extend from Yichang to Hukou, and the lower reaches extend from Hukou to the river mouth near Shanghai (Fig. 1).Cuntan forms the entrance to the Three Gorges Dam (TGD), which extends more than 600 km along the mainstream of the Yangtze River.The Three Gorges Dam was constructed 37 km upstream from Yichang for multiple purposes: energy generation, flood control, and water supply.

ERA-Interim
The ERA-Interim reanalysis dataset contains physical data of atmosphere and surface analyses covering the period from 1979 to the present based on the ECMWF Integrated Forecast System (IFS) release Cy31r2 (Berrisford et al., 2011;Simmons et al., 2006).The reanalysis incorporates a forecast model with three fully coupled components for atmosphere, land surface and ocean waves, and assimilates various types of observations, including satellite and ground based measurements.It uses the Tiled ECMWF Scheme for Surface Exchanges over Land (Viterbo et al., 1995) to simulate heat and water exchanges between land and atmosphere.The TESSEL model structure includes four soil layers (0-7, 7-28, 28-100, and 100-289 cm) for each type of vegetation scheme and each type of snow scheme.As the latest global atmospheric reanalysis produced by ECMWF, it has been confirmed that the performance of this system is substantially improved in certain key aspects (the representation of the hydrological cycle, the quality of the stratospheric circulation, and the consistency in time of the reanalyzed fields) compared to EAR-40 (Dee et al., 2011).
The monthly means, based on daily means of volumetric soil water in the four layers, as well as snow depth and sis data for 1994 to 1999, NOAA/GDAS atmospheric analysis fields for 2000, and a combination of NOAA/GDAS atmospheric analysis fields, spatially and temporally disaggregated NOAA Climate Prediction Center Merged Analysis of Precipitation (CMAP) fields, and observation-based downward shortwave and long wave radiation fields, using the method of the Air Force Weather Agency's AGRicultural METeorological modeling system(AGRMET), for the period 2001 to the present (Rui, 2011).
In this study, we use the model output produced by the Noah land surface scheme.These data are available from 1979 to the present at 3-hourly intervals.The Noah soil moisture profile includes four layers, namely 10, 30, 60, and 100 cm, from the soil surface down.The monthly products generated through temporal averaging of the 3-hourly products of soil moisture content(SMC) and snow water equivalent (SWE) during the period January, 1979 to December, 2010, are used in this study.Data generated by the GLDAS-Noah is publicly available on http://ldas.gsfc.nasa.gov/gldas/.

Field data
The monthly river discharge at the Yichang Gauging Station (Fig. 1) has been recorded over the period January 1979 to December 2004.This data set is used to validate the ERA-Interim land GLDAS-Noah outputs for the Yangtze River Basin.

GPCC
The Global Precipitation Climatology Centre (GPCC) offers gauge-based, gridded, monthly precipitation data sets for the global land surface from 1901 to 2010.The GPCC Full Data Reanalysis Version 6 with a spatial resolution of 1.0 • , which is fully independent of the precipitation data from ERA-200 Interim and GLDAS-Noah, is used in this study.It uses the complete GPCC station database (ca.67.200 stations with at least 10 years of data) available at the time of analysis and is therefore recommended for use in global and regional water balance studies, the calibration/validation of remote sensing 205 based rainfall estimations, and the verification of numerical models (Schneider et al., 2011).

PREC/L
NOAA's PRECipitation REConstruction over Land (PREC/L), a further gauge-based dataset regarding monthly 210 precipitation over land (Chen., 2002), is included in this study (spatial resolution of 2.5 • ).PREC/L is based on different collections of gauge data than are used for the GPCC, and draws on a large number of stations until 1990s.It has been documented by Wang et al. (2008) to have a 215 very high anomaly correlation coefficient (root mean square error) with the observation.

225
The RL05 L3-land data are based on the RL05 spherical harmonics from the CSR, the Jet Propulsion Laboratory (JPL) and the German Research Centre for Geosciences (GFZ),

GLDAS-Noah
The Global Land Data Assimilation System (GLDAS) supplies users with a model output of state-of-the-art land surface schemes created with atmospheric variables that originate from various data sources.The model has been forced by multiple datasets: bias-corrected ECMWF Reanalysis data for the time period 1979 to 1993, bias-corrected National Center for Atmospheric Research (NCAR) Reanalysis data for 1994 to 1999, NOAA/GDAS atmospheric analysis fields for 2000 and a combination of NOAA/GDAS atmospheric analysis fields, spatially and temporally disaggregated NOAA Climate Prediction Center Merged Analysis of Precipitation (CMAP) fields, and observation-based downward shortwave and longwave radiation fields, using the method of the Air Force Weather Agency's AGRicultural METeorological modeling system (AGRMET), for the period 2001 to the present (Rui, 2011).
In this study, we use the model output produced by the Noah land surface scheme.These data are available from 1979 to the present at 3 h intervals.The Noah soil moisture profile includes four layers, namely 10, 30, 60, and 100 cm, from the soil surface down.The monthly products generated through temporal averaging of the 3 h interval products of soil moisture content (SMC) and snow water equivalent (SWE) during the period January 1979 to December 2010 are used in this study.Data generated by the GLDAS-Noah is publicly available on http://ldas.gsfc.nasa.gov/gldas/.

Field data
The monthly river discharge at the Yichang Gauging Station (Fig. 1) has been recorded over the period January 1979 to December 2004.This dataset is used to validate the ERA-Interim land GLDAS-Noah outputs for the Yangtze River basin.

GPCC
The Global Precipitation Climatology Centre (GPCC) offers gauge-based, gridded, monthly precipitation datasets for the global land surface from 1901 to 2010.The GPCC Full Data Reanalysis Version 6 with a spatial resolution of 1.0 • , which is fully independent of the precipitation data from ERA-Interim and GLDAS-Noah, is used in this study.It uses the complete GPCC station database (ca.67 200 stations with at least 10 yr of data) available at the time of analysis and is therefore recommended for use in global and regional water balance studies, the calibration/validation of remote sensing based rainfall estimations, and the verification of numerical models (Schneider et al., 2011).

PREC/L
NOAA's PRECipitation REConstruction over Land (PREC/L), a further gauge-based dataset regarding monthly precipitation over land (Chen et al., 2002), is included in this study (spatial resolution of 2.5 • ).PREC/L is based on different collections of gauge data than are used for the GPCC, and draws on a large number of stations until the 1990s.It has been documented by Wang et al. (2008) (Bettadpur et al., 2012a,b) have proved that RL05 is more accurate than previously released GRACE products.The RL05 L3-land data are based on the RL05 spherical harmonics from the CSR, the Jet Propulsion Laboratory (JPL) and the German Research Centre for Geosciences (GFZ), and have additional, post processing steps, summarized on http://grace.jpl.nasa.gov/data/gracemonthlymassgridsland/.

Water storage estimation
TWS is generally defined as all phases of water stored above and below the surface of the Earth: soil moisture, canopy water storage, snow water equivalent and ground water, surface water storage, etc.Our analysis of storage is, however, limited to the total soil moisture column (TSM) and snow water equivalent (SWE) and does not give a complete description of the lateral and vertical distribution of water storage unless surface and groundwater components are added to the land model used here.We also neglect canopy water storage (CWS), although this is included in the GLDAS-Noah simulation.The reason is that CWS in the Yangtze River basin is negligible in comparison with soil moisture (Zhong et al., 2010).Therefore, TWS is expressed as Eq. ( 1), where N is an index representing the month of the year.
The monthly change in terrestrial water storage (TWSC N ) can be calculated at each pixel as follows: This method elicits promising results and also compares well with the Gravity Recovery and Climate Experiment (GRACE) estimation and the monthly basin-scale terrestrial water balance approach from flux variables (Chen et al., 2009;Rodell et al., 2004;Syed et al., 2008).The ERA-Interim soil profile includes four layers of 7, 21, 72, and 189 cm depth (forming a total of 289 cm), while the Noah soil profile includes four layers of 10, 30, 60, and 100 cm (200 cm in total).In order to be able to compare the TWS information obtained from both these datasets, we only considered the first 200 cm of soil in both cases.

Statistical analysis
Trend analyses involve linear regression and the nonparametric Mann-Kendall (MK) test (Mann, 1945;Kendall, 1975).A linear regression model is used to compute the annual trend in the TWS for each pixel.The MK test is a rankbased procedure and is applied to detect the significance of the trends.The MK test statistics are given by where The σ term is given by 1/18 n(n − 1)(n − 5) − t t (t − 1)(2t + 5) and X j and X i are the sequential data values, n is the dataset record length, t is the extent of any given tie (the number of annual maxima in a given tie), and is the summation of all ties.Positive and negative values of Z indicate increasing and decreasing trends, respectively.The statistic Z follows a normal distribution N (0,1) (Burn and Hag Elnur, 2002;Yang et al., 2010).To analyze whether the trend is stationary in the TWS anomalies, the Mann-Kendall-Sneyers (MKS) test (Sneyers, 1975) is also applied.This test, a sequential version of the MK test, enables not only detection of significant trends, but also approximation of the transition point in the temporal behavior of a series.Let x 1 , . .., x n be the data points.For each element x i , the number n i of element x j preceeds it (j < i) such that x j < x i is computed.Under the null hypothesis (no trend), the test statistic t k = k i=1 n i is normally distributed, with the mean and variance given by 2) be the normalized variable, which is the forward sequence.This principle can be usefully extended to the backward sequence UB k , which is calculated using the same equation but with a reversed series of data.When points in the forward series are outside the confidence interval, this indicates the detection of a significantly increasing (UF > 0) or a significantly decreasing (UF < 0) trend.If an intersection occurs between UF and UB within the confidence interval, this indicates an inflection (Li et al., 2004(Li et al., , 2007;;Moraes et al., 1998).Figure 2 shows that the ERA-Interim modeled runoff ts the observed values better than the GLDAS-Noah modled runoff does, for the period between 1979 and 2004.0.87 and 4.19mm, respectively, while for GLDAS-Noah , they are 0.68 and 14.58mm, respectively.Note that the runoff is consistently underestimated by GLDAS-Noah, which is also confirmed by Zaitchik et al. (2010).GLDAS-Noah outputs show errors in 1996 and 1997.

335
Apparently, ERA-Interim datasets show higher accuracy and reliability in the Yangtze River Basin.
To explore the quality of these datasets further and as precipitation arguably forms the most critical input into an accurate TWS, precipitation estimates of ERA-Interim and 340 GLDAS-Noah are compared with products from the GPCC and PREC/L, which are derived more directly from observations.The spatially averaged time series of standardized annual anomalies have been computed and compared for these four data sets.The result (see Fig. 3) shows a notable error 345    To explore the quality of these datasets further and as precipitation arguably forms the most critical input into an accurate TWS, precipitation estimates of ERA-Interim and 340 GLDAS-Noah are compared with products from the GPCC and PREC/L, which are derived more directly from observations.The spatially averaged time series of standardized annual anomalies have been computed and compared for these four data sets.The result (see Fig. 3) shows a notable error 345 in 1996 concerning GLDAS-Noah.GPCC and PREC/L fit very well (their R-squared value is 0.86).Generally speaking, ERA-Interim precipitation fits PREC/L and GPCC bet-  G-O , RMSE G-O ), they are 0.68 and 14.58 mm, respectively.Note that the runoff is consistently underestimated by GLDAS-Noah, which is also confirmed by Zaitchik et al. (2010).GLDAS-Noah outputs show errors in 1996 and 1997.Apparently, ERA-Interim datasets show higher accuracy and reliability for the Yangtze River basin.
To explore the quality of these datasets further and as precipitation arguably forms the most critical input into an accurate TWS, precipitation estimates of ERA-Interim and GLDAS-Noah are compared with products from the GPCC and PREC/L, which are derived more directly from observations.The spatially averaged time series of standardized annual anomalies have been computed and compared for these four datasets.The result (see Fig. 3) shows a notable error in 1996 concerning GLDAS-Noah.GPCC and PREC/L fit very well (their R-squared value is 0.86).Generally speaking, ERA-Interim precipitation fits PREC/L and GPCC better than GLDAS-Noah does.The R-squared between ERA-Interim and PREC/L (R 2 E-P ) and between ERA-Interim and GPCC (R 2 E-G ) are 0.49 and 0.66, respectively, while the Rsquared between GLDAS-Noah and PREC/L (R 2 G-P ) and between GLDAS-Noah and GPCC (R 2 G-G ) are 0.18 and 0.13, respectively.ERA-Interim generally shows good agreement with GPCC and PREC/L; however, a small shift has become apparent in the past decade, which we will discuss later.
In situ measurements of soil moisture are invaluable for the calibration and validation of a land surface model and satellite-based soil moisture retrieval.Unfortunately, there is a very low sampling rate with only 1 sample being available in the Yangtze River basin from the International Soil Moisture Network (ISMN) (Dorigo et al., 2011).However, the error structures of the ERA-Interim and GLDAS-Noah soil moisture products have been estimated using the triple collocation technique by Dorigo et al. (2010) and Scipal et al. (2008).ERA-Interim reanalyzed soil moisture is characterized by a relatively low mean global error of 0.018 m 3 m −3 (Dorigo et al., 2010), which is fairly consistent with the average error (a mean global error of 0.020 m 3 m −3 ) obtained by Scipal et al. (2008) by applying the triple collocation model to three satellite-based and model-based soil moisture products.It is found that the errors of soil moisture estimates in the Yangtze River basin are at an intermediate level.This can also be confirmed by the high correlation with ASCAT retrievals for the years 2007 and 2008 (Dorigo et al., 2010) and ERS-2 retrievals for the years 1998, 1999, and 2000(Scipal et al., 2008)).In addition, Liu et al. (2011) has shown that there is a high correlation coefficient (R) between GLDAS-Noah and ASCAT retrievals for the Yangtze River basin in 2007.It has been firmly proven that active microwave satellite-based (e.g.ASCAT) retrievals result in smaller errors in moderately to densely vegetated areas (e.g. the Yangtze River basin) than passive microwave products do (Liu et al., 2011).Therefore, the high correlation between ERA-Interim, or GLDAS-Noah, and active microwave satellite-based soil moisture retrievals provides some confidence in the EAR-Interim and GLDAS-Noah soil moisture qualities in the Yangtze River basin.
Other components such as surface water and groundwater form a large proportion of the TWS.To assess their impact on the matchup, we compared TWS products derived from ERA-Interim to those derived from GRACE observations (reprocessed Release-05, GRACE RL05) for a sevenyear period (2004)(2005)(2006)(2007)(2008)(2009)(2010). Figure 4 shows that the magnitude of the spatially averaged TWS anomalies from these two datasets (ERA-Interim and GRACE RL05) is similar and exhibits the same variation, with a coefficient of determination as high as 0.79.This means that the ERA-Interim product on TWS over a soil depth of 2 m is representative for the GRACE observations that are affected by water storage fluctuations in the entire air-land column, including surface water and groundwater.

Climatology
Annual standardized anomalies are calculated by using the monthly values subtracted with the annual mean and divided by the annual standard deviation of the annual mean, which can be expressed as follows: with where A ij is the annual TWS standardized anomaly in the ith month of the j th year, TWS ij is the TWS in the ith month of the j th year, TWS j is the mean TWS of the all months in the j th year, and σ j is the standardized deviation of all months in the j th year.The spatial distribution of TWS and TWSC climatological annual standardized anomalies derived from ERA-Interim and GLDAS-Noah are shown in Figs. 5 and 6, respectively.The performance of ERA-Interim and GLDAS-Noah elicits large variability, to a large extent due to insufficient physical descriptions of land-surface processes (Niu et al., 2006;Zeng et al., 2008).Nevertheless, the spatial patterns reveal strong consistencies.
After suffering the dry season (December-February), the southeast corner of the Yangtze basin starts to become wet during March to May, due mainly to the south China rainfall belt extension and the mean precipitation increase in the lower basin (Ding and Chan, 2005;Qian et al., 2002).High, positive TWS standardized anomalies emerge in most of the Yangtze basin during June to October (Fig. 5), and show a large increase in July compared to June (Fig. 6), corresponding with the intensive precipitation observed along the whole Yangtze River from mid-June to mid-July (Ding, 1992) sharply, while they stay quite high and positive in the upper Yangtze reaches till October, mainly due to the continuous rainy season from mid-June to mid-September.
This striking consistency between TWS and the rainfall pattern is not unexpected.That higher precipitation leads to higher soil moisture can generally be considered predictable, though there are a few exceptions.For instance, in the case of intense precipitation with rates beyond the infiltration rate, or precipitation over very wet or saturated areas, the rainfall anomalies will result in runoff anomalies rather than soil moisture (Dunne et al., 1978;Horton et al., 1933).Nevertheless, except in the extreme cases, there is an obvious and direct response of soil moisture to precipitation.On the other hand, the feedback, via the return path from soil moisture through evapotranspiration to precipitation, can also play an important role in the TWS variability, though a weaker one (Seneviratne et al., 2010).Abundant previous research (Dirmeyer, 2011;Jung et al., 2010;Wei et al., 2012) shows that the Yangtze River Basin is dominated by wet soil moisture regimes, where soil moisture does not mainly control the variability in evapotranspiration and has only a minor impact on the change in rainfall.Dirmeyer (2011) also confirms that soil moisture neither typically provides feedback to the atmosphere nor has a damping effect on climate variability.Thus, it is reasonable to speculate that the TWS variability in the Yangtze River Basin is mainly controlled by large-scale atmospheric circulations, as is also established by Wei et al. (2012).Moreover, as displayed in Fig. 5, the Yangtze River Basin suffered the highest TWS anomalies during June-July, which implies a high flood risk during this period, since runoff is sensitive to soil moisture content under wet soil regimes.When soil moisture is very high and soil becomes saturated, high precipitation variability may lead to high runoff variability, which cannot be damped by soil moisture storage (Seneviratne et al., 2010).dle and lower Yangtze reaches, which may be explained b large-scale circulation and heterogeneous land-surface con ditions.The upper Yangtze reaches are mainly influence by the South Asian (or Indian) summer monsoon, and th 495 middle and lower Yangtze reaches are controlled by the Eas Asian summer monsoon (Ding et al., 2005).The seasona process of the Asian summer monsoon plays a crucial rol in heat and moisture transport and the hydrological cycle Related rainfall systems perform differently in the uppe 500 reaches than in the middle and lower reaches (Qian et al 2002).Since the topography in the upper Yangtze reache is totally different from that in the middle and lower reache (Fig. 1), land-surface heterogeneities in temperature are ex pected (Giorgi et al., 1997;Salama et al., 2012).The lan 505 cover and hydrological conditions differ in the two area (Piao et al., 2010).The inhomogeneous surface results i heterogeneity in surface energy partitioning, which in tur has an impact on land-atmosphere interactions (Brunsell e al., 2011;Ma et al., 2008).Therefore, different response 510 are expected by land-surface systems in the upper Yangtz reaches to in the middle-lower reaches.For example, so moisture exerts a significant positive control on the max mum and mean temperature in the middle-lower reaches du ing summer, while no significant control is elicited in the up 515 per reaches.Furthermore, while soil moisture and precip tation are positively coupled in the upper Yangtze reache this coupling is negative for the middle and lower Yangtz reaches (Zhang et al., 2011).2010) pointed out that the major El Nin event in 1998 was followed by changes in the behaviou of some land water cycle components, such as the SMC (around 107 • E) for ERA-Interim, which is different from the negative GLDAS-Noah value in the upper area.According to the climatological rainfall differences between May and June and between June and July (Qian et al., 2002, Fig. 4), the increased rainfall in June compared to May appears in the Plateau and in southwest China, with the center lying in the upper Yangtze reaches.

TWS trend analysis
Another area of increased rainfall is located along the eastern coastland with its center in the lower Yangtze reaches.There is no obvious increase in rainfall in the central part of China.In July, the increased rainfall has migrated to the north of the lower Yangtze River basin, while rainfall is steadily increasing in the upper parts.This pattern of change in precipitation from May to July resembles the TWS pattern derived from ERA-Interim (Fig. 5).After July, the TWS anomalies in the middle and lower Yangtze reaches decrease sharply, while they stay quite high and positive in the upper Yangtze reaches till October, mainly due to the continuous rainy season from mid-June to mid-September.
This striking consistency between TWS and the rainfall pattern is not unexpected.That higher precipitation leads to higher soil moisture can generally be considered predictable, though there are a few exceptions.For instance, in the case of intense precipitation with rates beyond the infiltration rate, or precipitation over very wet or saturated areas, the rainfall anomalies will result in runoff anomalies rather than soil moisture (Dunne, 1978;Horton et al., 1933).Nevertheless, except in the extreme cases, there is an obvious and direct response of soil moisture to precipitation.On the other hand, the feedback, via the return path from soil moisture through evapotranspiration to precipitation, can also play an important role in the TWS variability, though a weaker one (Seneviratne et al., 2010).Abundant previous research (Dirmeyer, 2011;Jung et al., 2010;Wei et al., 2012) shows that the Yangtze River basin is dominated by wet soil moisture regimes, where soil moisture does not mainly control the variability in evapotranspiration and has only a minor impact on the change in rainfall.Dirmeyer (2011) also con-firms that soil moisture neither typically provides feedback to the atmosphere nor has a damping effect on climate variability.Thus, it is reasonable to speculate that the TWS variability in the Yangtze River basin is mainly controlled by large-scale atmospheric circulations, as is also established by Wei et al. (2012).Moreover, as displayed in Fig. 5, the Yangtze River basin suffered the highest TWS anomalies during June-July, which implies a high flood risk during this period, since runoff is sensitive to soil moisture content under wet soil regimes.When soil moisture is very high and soil becomes saturated, high precipitation variability may lead to high runoff variability, which cannot be damped by soil moisture storage (Seneviratne et al., 2010).
It should be recognized that the TWS pattern in the upper Yangtze reaches is completely different from that in the middle and lower Yangtze reaches, which may be explained by large-scale circulation and heterogeneous land-surface conditions.The upper Yangtze reaches are mainly influenced by the South Asian (or Indian) summer monsoon, and the middle and lower Yangtze reaches are controlled by the East Asian summer monsoon (Ding and Chan, 2005).The seasonal process of the Asian summer monsoon plays a crucial role in heat and moisture transport and the hydrological cycle.Related rainfall systems perform differently in the upper reaches than in the middle and lower reaches (Qian et al., 2002).Since the topography in the upper Yangtze reaches is totally different from that in the middle and lower reaches (Fig. 1), land-surface heterogeneities in temperature are expected (Giorgi et al., 1997;Salama et al., 2012).The land cover and hydrological conditions differ in the two areas (Piao et al., 2010).The inhomogeneous surface results in heterogeneity in surface energy partitioning, which in turn has an impact on land-atmosphere interactions (Brunsell et al., 2011;Ma et al., 2008).Therefore, different responses are expected by land-surface systems in the upper Yangtze reaches than in the middle/lower reaches.For example, soil moisture exerts a significant positive control on the maximum and mean temperature in the middle/lower reaches during summer, while no significant control is elicited in the upper reaches.Furthermore, while soil moisture and precipitation are positively coupled in the upper Yangtze reaches, this coupling is negative for the middle and lower Yangtze reaches (Zhang et al., 2011).Jung et al. (2010) pointed out that the major El Nino event in 1998 was followed by changes in the behavior of some land water cycle components, such as the SMC.Additionally, GLDAS-Noah outputs show errors in 1996 and 1997.Therefore, we only show the results of the ERA-Interim dataset and separate the whole study period into two parts, from January 1979 to December 1997, and from January 1998 to December 2010.In Fig. 7, the ERA-Interim dataset shows decreasing TWS trends over large parts of the Yangtze River basin between 1998 and 2010, which match the descending trend in the SMC from microwave satellite observations between 1998 to 2008 (Jung et al., 2010).It shows significantly decreasing trends (most of which surpass the 95 %, while some even surpass the 99 % confidence level) in the middle and lower reaches, with a maximum of −3.93 mm yr −1 .The upper reaches suffer milder decreases and even insignificant trends in some parts during the period 1998 to 2010.Between 1979 and 1997, it provides insignificant trends for most regions of the basin.This result indicates that the Yangtze River basin is drying up, the conclusion also reached by a new WWF study (http://www.asianscientist.com/topnews/yangtze-river-basin-is-drying-up-wwf-china-2012/).

TWS trend analysis
Monthly standardized anomalies of the TWS are calculated by the monthly TWS minus the corresponding monthly value of the annual cycle, and then divided by the standard deviation of the values of the same months within the period 1979 to 2010, in order to eliminate the influence of interannual variability for intra-annual analysis.Monthly standardized anomalies of the TWS can be expressed as follows: with where M ij is the monthly TWS standardized anomaly in the ith month of the j th year.The subscripts i and j represent the ith month and j th year, respectively; TWS i is the TWS of the ith month averaged over all the years; σ i is the standardized deviation of ith month TWS over all the years.
As discussed before, the ERA-Interim dataset has better accuracy, at least in the upper reaches of the Yangtze River, while GLDAS-Noah outputs show notable errors in 1996 and 1998.Therefore, only the ERA-Interim dataset is used in this section.
The MKS test is applied to detect the transition points in the temporal behavior of TWS standardized anomalies based on the annual mean, the wet season mean and the dry season mean, respectively.The definition of wet season and dry season is based on the precipitation climatology of the Yangtze River basin.The Yangtze River basin experiences a distinct wet season from about May to late September or early October.The corresponding dry season spans from late September or early October to spring.The summer monsoons contribute most of the wet season precipitation (Harvey et al., 2007).As seen in Fig. 8, generally speaking, the spatially averaged TWS standardized anomaly trends are not significant (< 95 % confidence level) and not monotonic (i.e. with a transition point) during the period 1979 to 2010.The only transition point during the 32 yr period of 1979 to 2010, at which point the TWS standardized anomalies began to decrease sharply, occurred in 2006.This trend reaches the 95 % confidence level in 2010.The transition point occurs one year earlier in the wet season, and two years later in the dry season.In the middle and lower Yangtze reaches, the transition point occurs around 2005 both in the wet and in the dry season.It is noted that there is a significant downward trend in 2009 and 2010 after four years of insignificant decrease.It was the first time this happened since the start of the study in 1979.In the upper Yangtze reaches, the TWS standardized anomalies experience mainly downward trends during the wet seasons of the past 3 decades and increasing trends during the dry seasons.In addition, transition points occur several times (1982, 1989, 1995, 2001, 2005,. . . ) and in the period of 1986-1988, the decrease is significant in the wet season in the upper Yangtze reaches.We also examined the transition points through MKS of the TWS standardized anomalies in the middle and lower Yangtze reaches during June-July, and the result is exactly the same as in the wet season, though the TWS standardized anomalies do differ from each other.
As seen in Figs. 8 and 9, the past 6 yr period (2005-2010) was the driest period in the Yangtze River basin (especially in the middle and lower reaches) since 1979.This result is quite consistent with the severe drought events documented for the basin by other research.Wei et al. (2012) documented that the Yangtze River basin suffered one of the driest rainy seasons during the 32 yr period of 1979 to 2010 in 2005, and Yan et al. (2007) noted that a widespread drought occurred over the southwestern part of the basin that same spring, and that it was the most serious drought since 1979 till 2007.Then the worst drought in more than a century struck southwest China and Sichuan in the summer of 2006, and Dai et al. (2008)    the 99% confidence level) in the middle and lower reaches, with a maximum of -3.93 mm year −1 .The upper reaches suffer milder decreases and even insignificant trends in some parts during the period 1998 to 2010.Between 1979 and 1997, it provides insignificant trends for most regions of the basin.This result indicates that the Yangtze River Basin is drying up, conclusion also reached by a new WWF study (http://www.asianscientist.com/topnews/yangtzeriver-basin-is-drying-up-wwf-china-2012/).
Monthly standardized anomalies of the TWS are calculated by the monthly TWS minus the corresponding monthly value of the annual cycle, and then divided by the standard deviation of the values of the same months within the period 1979 to 2010, in order to eliminate the influence of inter annual variability for intra annual analysis.Monthly standardized anomalies of the TWS can be expressed as follows: with: where M ij is the monthly TWS standardized anomaly in the ith month of the jth year.The subscripts i and j represent the ith month and jth year, respectively; T W S i is the TWS of the ith month averaged over all the years; σ i is the 560 standardized deviation of ith month TWS over all the years.
As discussed before, the ERA-Interim dataset has better accuracy, at least in the upper reaches of the Yangtze River, while GLDAS-Noah outputs show notable errors in 1996 and 1998.Therefore, only the ERA-Interim dataset is used in this 565 section.
The MKS test is applied to detect the transition points in the temporal behaviour of TWS standardized anomalies based on the annual mean, the wet season mean and the dry season mean, respectively.The definition of wet season and 570 dry season is based on the precipitation climatology of the Yangtze River Basin.The Yangtze River Basin experiences a distinct wet season from about May to late September or early October.The corresponding dry season spans from late September or early October to spring.The summer mon-575 soons contribute most of the wet season precipitation (Harvey et al., 2007).As seen in Fig. 8, generally speaking, the spatially averaged TWS standardized anomaly trends are not significant (< 95% confidence level) and not monotonic (i.e., with a transition point) during the period 1979 to 2010.The In the upper Yangtze reaches, the TWS standardized anomalies experience mainly downward trends during the wet seasons of the past 3 decades and increasing trends during the dry seasons.In addition, transition points occur several times (1982,1989,1995,2001,2005) and in 595 the period of 1986-1988, the decrease is significant in the wet season in the upper Yangtze reaches.We also examined the transition points through MKS of the TWS standardized anomalies in the middle and lower Yangtze reaches during June-July, and the result is exactly the same as in the wet sea-600 son, though the TWS standardized anomalies do differ from each other.
As seen in Fig. 8 and Fig. 9, the past 6-year period (2005-2010) was the driest period in the Yangtze River Basin (especially in the middle and lower reaches) since 1979.This 605 result is quite consistent with the severe drought events documented for the basin by other research.Wei et al. (2012) documented that the Yangtze River Basin suffered one of flood season.In 2007, the area around the Yangtze River suffered a severe drought again.In some places the water levels of the river dropped to their lowest levels since records began 142 yr ago.The drought was also severe in large areas of the normally wet south.Reservoirs and rivers shrunk and supplies of drinking water fell to alarmingly low levels.However, he extreme drought of 2009/2010 over southwestern Yangtze (including Yunnan, Sichuan and Guizhou) is the driest meteorological event with the lowest percentage rainfall anomaly and the longest rain-free period occurring during a winter season (October-February) in the past 50 years, and also the severest one with the lowest percentage rainfall anomaly since 1880, as documented by Yang et al. (2012) (http://factsanddetails.com/china.php?itemid= 1879&catid=10&subcatid=64).
Monthly standardized anomalies of precipitation from ERA-Interim, the GPCC, and PREC/L have been computed and compared to monthly standardized anomalies of TWS (not shown here) to examine possible correlation.The correlation between TWS anomalies from ERA-Interim and precipitation anomalies from ERA-Interim, the GPCC, and PREC/L concerning the Yangtze basin is reasonably high in the wet season (0.69, 0.53, and 0.49, respectively), but much lower, especially for GPCC and PREC/L, in the dry season (0.48, 0.21, and 0.25, respectively).From a regional perspective, the middle and lower Yangtze reaches exhibited greater agreement between TWS from ERA-Interim and precipitation from the three datasets, than the upper Yangtze reaches did.The notable negative TWS anomalies in the middle and lower Yangtze reaches are in clear agreement with the significant decrease in precipitation seen in the ERA-Interim data in the past 6 yr; the GPCC and PREC/L data exhibit more gentle negative precipitation anomalies during this period and do not show any special difference to the prior period (not shown here).The differences can also be seen clearly in Fig. 3, where there is a clear downward shift for ERA-Interim relative to both the GPCC and PREC/L data in the past 6 yr.This shift matches the general decline in values relative to the GPCC data for the past decade, which may be caused by a too low sea surface temperature (SST) in the ERA-Interim dataset, or fewer stations in the GPCC archive in recent years (Simmons et al., 2010).PREC/L uses fewer gauging stations since the 1990s as well, thus it is difficult to assess the recent huge and sudden drop in the ERA-Interim figures only by comparing them to the GPCC and PREC/L data.However, the dramatic precipitation decrease in the middle and lower Yangtze reaches had been examined by Zhu et al. (2011, Fig. 2) over the past decade (2000)(2001)(2002)(2003)(2004)(2005)(2006)(2007)(2008), and the rainfall anomalies based on a 160-station precipitation dataset of the last 58 yr  dropped sharply from positive to negative values around 2004.This pattern of precipitation is consistent with the dramatic decrease of precipitation seen in ERA-Interim data, changing from positive to negative values since 2004 in the middle and lower Yangtze reaches (Fig. 3), suggesting that the recent drop in precipitation is most likely the biggest contributor to the massive decline in TWS in the middle and lower Yangtze reaches.Ding et al. (2008) pointed out that the recent change in the summer rainfall pattern in the Yangtze River is strongly related to the variability of the East Asian summer monsoon (EASM) through its moisture transport and supply.Zhu et al. (2011) stated that the eastward recession of the western Pacific subtropical high (WPSH) and the significant changes in the global SST are the main causes of the rainfall deficit in the Yangtze River basin since the year 2000.Yan et al. (2007) and Liu et al. (2007) documented that the intensification and westward shift of the WPSH and the easterly anomaly over the northern Indian Ocean are two key causes of the 2005-spring drought over southwestern China.The northward shift of the WPSH and the negative snow cover anomaly over the Tibetan Plateau are important contributors to the 2006-summer drought (Zhu and Gao, 2007;Li et al., 2009)  TGD and the consistent drought in recent years, even though there has been no irrefutable evidence to prove that the TGR 720 is responsible for the extremely driest period that has occurred in the past several years, as the TGD has only been in operation for a short period.Apart from the TGR, numerous other reservoirs within the Yangtze catchment together reached 200 km 3 (Yang et al., 2005), more than five times 725 the storage capacity of the TGR.The impact of these reservoirs on the TWS should not be neglected.The Yangtze basin has witnessed remarkable changes in land use and cover induced by high population density and rapid but uneven economic growth (Long et al., 2007;Yin et al., 2010).These   TGD and the consistent drought in recent years, even though there has been no irrefutable evidence to prove that the TGR is responsible for the extremely driest period that has occurred in the past several years, as the TGD has only been in operation for a short period.Apart from the TGR, numerous other reservoirs within the Yangtze catchment together reached 200 km 3 (Yang et al., 2005), more than five times the storage capacity of the TGR.The impact of these reservoirs on the TWS should not be neglected.The Yangtze basin has witnessed remarkable changes in land use and cover induced by high population density and rapid but uneven economic growth (Long et al., 2007;Yin et al., 2010).These changes might alter the soil properties and soil-climate in-teractions, probably having great influence on the TWS and runoff distribution.It should be pointed out that the ERA-Interim TWS could contain significant uncertainties, as it relies heavily on satellite observations and modeling.Further  Since ERA-Interim TWS products do not include the impact of anthropogenic activity, such as the TGD and land cover change, in the model structure but rather in the assimilated observations, their effect on the regional climate in the Yangtze River basin is not obvious.However, some studies have tried to demonstrate the extent of human impact on the Yangtze River basin.For example, Dai et al. (2008) and Yang et al. (2010) show that the TGD reservoirs could have a direct impact on the intra-annual changes in the downstream Yangtze discharges, leading to a dumping of the seasonal variations in the Yangtze River discharge in the middle and lower reaches.Miller et al. (2005) and Wu et al. (2006) also documented that the land use change associated with the TGD would alter the regional pattern of precipitation, wind, and temperature.It could also impact the hydrologi-cal cycle of the river basin, and may lead to changes in the soil-climate interaction, which would probably alter the current dumping effect of soil wetness on the climate variability.As shown in Figs. 8 and 9, the consistent droughts in recent years and the operation of the TGD have occurred simultaneously.In 2003, the water level of the TGR reached 135 m.Coincidently, in 2004, the driest period of the past 32 yr began for the middle and lower Yangtze.Also, the whole basin suffered an abrupt change in 2006, when the TGR raised its water level from 135 to 156 m.This coincidence is very striking and may imply the possible connection between the TGD and the consistent droughts in recent years, even though there has been no irrefutable evidence to prove that the TGR is responsible for the extremely driest period that has occurred in the past several years, as the TGD has only been  (Yang et al., 2005), more than five times the storage capacity of the TGR.The impact of these reservoirs on the TWS should not be ignored.The Yangtze basin has witnessed remarkable changes in land use and cover induced by high population density and rapid but uneven economic growth (Long et al., 2007;Yin et al., 2010).These changes might alter the soil properties and soil-climate interactions, probably having great influence on the TWS and runoff distribution.It should be pointed out that the ERA-Interim TWS could contain significant uncertainties, as it relies heavily on satellite observations and modeling.Further investigation and analysis is needed to assess the significant impact of human activity on the TWS of the Yangtze River basin.

Conclusions
This study analyzes the spatial and temporal variations of the TWS in the Yangtze River basin during the period 1979 to 2010 based on ERA-Interim and GLDAS-Noah datasets.After comparison with field measurement data, ERA-Interim is found to perform better in the Yangtze River basin than GLDAS-Noah.Linear regression and the MK test have been used to detect trends and the significance of trends in TWS for each pixel in the whole basin.In addition, the MKS test has been used to detect transition points in the temporal series of the spatially averaged TWS annual standardized anomalies in the upper, middle/lower, and whole basin, respectively.We conclude that the TWS variations over the Yangtze River basin during the period 1979 to 2010 have the following characteristics: -Most of the Yangtze basin exhibits the highest positive TWS anomalies during June-July, mainly due to the Meiyu rain event.This intensified rainfall not only results in very high positive TWS anomalies, but probably also leads to high runoff anomalies, causing floods across the basin.
-The Yangtze River basin is drying up, especially since the year 1998.The TWS variation is strongly correlated to the precipitation variation derived from ERA-Interim, the GPCC and PREC/L data, especially during the wet season and in the middle and lower reaches, suggesting the TWS variation is mainly controlled by precipitation.
-In the middle and lower Yangtze reaches, the TWS behavior changed abruptly and started to decrease in 2004.Coincidently, the TGR started impoundment in 2003.
The year 2006 is detected as only transition point, initiating the TWS to suffer high negative anomalies, while the TGR raised its water level from 135 to 156 m that same year.Due to the short period of the TGD being in operation, it is hard to prove a link between the TGR and the consistent drought in recent years, but the coincidence suggests a possible connection.
-The past 6 yr period (2005)(2006)(2007)(2008)(2009)(2010) was the driest period in terms of TWS in the Yangtze River basin (especially in the middle and lower Yangtze reaches) since 1979.This is mainly the result of a dramatic decrease in precipitation, which may also be influenced by human activity.

Fig. 1 .
Fig. 1.Elevation map of the study area, the Yangtze River Basin.The green stars denote the locations of the main gauging stations.The position of the Three Gorges Dam is depicted by a red diamond.

Fig. 1 .
Fig. 1.Elevation map of the study area, the Yangtze River basin.The green stars denote the locations of the main gauging stations.The position of the Three Gorges Dam is depicted by a red diamond.
Spatial averaged time series of ERA-Interim estimated (red curve), GLDAS-Noah estimated (blue curve) and observed (black curve) unoff of the upper Yangtze reaches between January 1979 and December 2004.ig. 3. Spatially averaged time series of standardized anomalies of the annual precipitation in the middle and lower Yangtze reaches, based n PREC/L(blue curve), GPCC(red curve), ERA-Interim(black circle curve) and GLDAS-Noah (gray diamond curve) data from 1979 to 010.Firstly, we computed the accumulated monthly runoff rom ERA-Interim/GLDAS-Noah data at each pixel during he period 1979 to 2004.Secondly, we calculated the spatialean of the accumulated monthly runoff (mm) of all pixels ocated in the upper reaches of the Yangtze Basin.b) Discharge at the Yichang gauging station Firstly, we computed the accumulated monthly discharge m 3 ) from the daily discharge data (m 3 /s) of the Yichang tation.Secondly, we divided this figure by the area of the pper Yangtze reaches.The second step is supported by the act that the Yichang station forms the exit point of the upper eaches of the Yangtze River Basin.

Fig. 2 .
Fig. 2. Spatial averaged time series of ERA-Interim estimated (red curve), GLDAS-Noah estimated (blue curve) and observed (black curve) runoff of the upper Yangtze reaches between January 1979 and December 2004.

Fig. 2 .
Fig. 2. Spatial averaged time series of ERA-Interim estimated (red curve), GLDAS-Noah estimated (blue curve) and observed (black curve) runoff of the upper Yangtze reaches between January 1979 and December 2004.

Fig. 3 .
Fig. 3. Spatially averaged time series of standardized anomalies of the annual precipitation in the middle and lower Yangtze reaches, based on PREC/L(blue curve), GPCC(red curve), ERA-Interim(black circle curve) and GLDAS-Noah (gray diamond curve) data from 1979 to 2010.
Figure 2 shows that the ERA-Interim modeled runoff 325 fits the observed values better than the GLDAS-Noah modeled runoff does, for the period between 1979 and 2004.The coefficient of determination (R-squared) and the root mean square error (RMSE) between the modeled and observed values for ERA-Interim (R 2 E−O , RM SE E−O ) are 330 0.87 and 4.19mm, respectively, while for GLDAS-Noah (R 2 G−O , RM SE G−O ), they are 0.68 and 14.58mm, respectively.Note that the runoff is consistently underestimated by GLDAS-Noah, which is also confirmed by Zaitchik et al. (2010).GLDAS-Noah outputs show errors in 1996 and 1997.335 Apparently, ERA-Interim datasets show higher accuracy and reliability in the Yangtze River Basin.

Fig. 3 .
Fig. 3. Spatially averaged time series of standardized anomalies of the annual precipitation in the middle and lower Yangtze reaches, based on PREC/L (blue curve), GPCC (red curve), ERA-Interim (black circle curve) and GLDAS-Noah (gray diamond curve) data from 1979 to 2010.

Figure 2
Figure2shows that the ERA-Interim modeled runoff fits the observed values better than the GLDAS-Noah modeled runoff does, for the period between 1979 and 2004.The coefficient of determination (R-squared) and the root mean square error (RMSE) between the modeled and observed values for ERA-Interim (R 2 E-O , RMSE E-O ) are 0.87 and 4.19 mm, respectively, while for GLDAS-Noah (R 2 G-O , RMSE G-O ), they are 0.68 and 14.58 mm, respectively.Note that the runoff is consistently underestimated by GLDAS-Noah, which is also confirmed byZaitchik et al. (2010).GLDAS-Noah outputs show errors in 1996 and 1997.Apparently, ERA-Interim datasets show higher accuracy and reliability for the Yangtze River basin.
Fig. 4. TWS anomalies [cm] averaged for a seven years period (2004-2010) and obtained from ERA-Interim (red line) and GRACE RL0 (blue line) datasets for the Yangtze River Basin.

Fig. 5 .Fig. 5 .Fig. 6 .
Fig. 5. Spatial patterns of monthly averaged TWS annual standardized anomalies computed from ERA-Interim for the period January 1979 till December 2010 and from GLDAS-Noah for the periods January 1979 till December 1995 and January 1998 till December 2010.

Fig. 6 .
Fig. 6.Spatial patterns of monthly averaged TWSC annual standardized anomalies computed from ERA-Interim for the period January 1979 till December 2010 and from GLDAS-Noah for the periods January 1979 till December 1995 and January 1998 till December 2010.
580 only transition point during the 32-year period of 1979 to 2010, at which point the TWS standardized anomalies began to decrease sharply, occurred in 2006.This trend reaches the 95% confidence level in 2010.The transition point occurs one year earlier in the wet season, and two years later in the 585 dry season.In the middle and lower Yangtze reaches, the transition point occurs around 2005 both in the wet and in the dry season.It is noted that there is a significant downward trend in 2009 and 2010 after four years of insignificant decrease.It was the first time this happened since the start 590 of the study in 1979.

Fig. 7 .
Fig. 7. ERA-Interim estimated TWS annual trends between 1979 and 1997, and between 1998 and 2010, in millimeters per year.Gray grids cells represent insignificant trends; cells with an empty diamond indicate the trend surpasses the 95 % confidence level; cells with a filled diamond indicate the trend surpasses the 99 % confidence level; others indicate the trend surpasses the 90 % confidence level.

Fig. 8 .Fig. 9 .
Fig. 8.The forward (UF, red) and backward (UB,blue) Mann-Kendall statistic rank series for TWS standardized anomalies, per year (a, d, g), per wet season (b, e, h), and per dry season (c, f, i) during the period 1979 to 2010 for the upper reaches(a, b, c), the middle-lower reaches(d, e, f) and the whole (g, h, i) of the Yangtze River Basin(the horizontally dotted lines represent the critical values corresponding to the 95% confidence level)

730
changes might alter the soil properties and soil-climate in-teractions, probably having great influence on the TWS and runoff distribution.It should be pointed out that the ERA-Interim TWS could contain significant uncertainties, as it relies heavily on satellite observations and modeling.Further 735 investigation and analysis is needed to assess the significant impact of human activity on the TWS of the Yangtze River Basin.6 ConclusionsThis study analyzes the spatial and temporal variations of the 740 TWS in the Yangtze River Basin during the period 1979 to 2010 based on ERA-Interim and GLDAS-Noah datasets.Af-

Fig. 8 .
Fig. 8.The forward (UF, red) and backward (UB,blue) Mann-Kendall statistic rank series for TWS standardized anomalies, per year (a, d, g), per wet season (b, e, h), and per dry season (c, f, i) during the period 1979 to 2010 for the upper reaches (a, b, c), the middle/lower reaches (d, e, f) and the whole (g, h, i) of the Yangtze River basin (the horizontally dotted lines represent the critical values corresponding to the 95 % confidence level).

Fig. 8 .
Fig. 8.The forward (UF, red) and backward (UB,blue) Mann-Kendall statistic rank series for TWS standardized anomalies, per year (a, d, g), per wet season (b, e, h), and per dry season (c, f, i) during the period 1979 to 2010 for the upper reaches(a, b, c), the middle-lower reaches(d, e, f) and the whole (g, h, i) of the Yangtze River Basin(the horizontally dotted lines represent the critical values corresponding to the 95% confidence level)

735
investigation and analysis is needed to assess the significant impact of human activity on the TWS of the Yangtze River Basin.6 ConclusionsThis study analyzes the spatial and temporal variations of the 740 TWS in the Yangtze River Basin during the period 1979 to 2010 based on ERA-Interim and GLDAS-Noah datasets.Af-

Fig. 9 .
Fig. 9.The time series of the annual, wet season (May-October) and dry season averaged (November-April) TWS standardized anomalies for the upper Yangtze Reaches, the middle/lower Yangtze reaches and the whole Yangtze River Basin, respectively.

Hydrol. Earth Syst. Sci., 17, 1985-2000, 2013 www
showed that the middle and lower Yangtze River reaches suffered the lowest level of the past 50 yr during that .hydrol-earth-syst-sci.net/17/1985/2013/ 8 Y.Huang et al.: Analysis of Long-term Terrestrial Water Storage Variations in Yangtze River Basin . The extreme drought event of 2009/2010 over southwestern China is associated with the westward extension of the WPSH brought about by the Arabian Sea cyclonic anomaly and the El Nino Modoki event during 2009/2010.