Abstract

The Gravity Recovery and Climate Experiment (GRACE) satellite mission provides measurements of Earth’s static and time-variable gravity fields with monthly resolution. In this study, changes of water storage in northwestern China were determined by GRACE monthly gravity field data obtained from 2003 to 2010. Comparisons of water storage change (WSC) simulated by a four-dimensional assimilation model (Noah) and observed by GRACE revealed similar patterns of change and a correlation coefficient of 0.71 . Trend analysis indicated significant changes in the spatiotemporal variation of WSC in northwestern China during the 8-year study period, which were stronger in the east than in the west and more pronounced in the south than in the north. The most pronounced increase in water storage occurred in Gansu and Qinghai provinces, but, overall, water storage increased by 0.61 mm/a over northwestern China during the study period. Clear seasonal variations of WSC and precipitation were found, because glacial meltwater and precipitation are the main sources of water in the hydrosphere; meanwhile, the distributions of glaciers and permafrost also affect the spatial distribution of WSC.

1. Introduction

The cryosphere is an important component of the global water cycle, and its sensitivity to climate change has been a subject of much concern. In particular, since the 1980s, the global cryosphere has undergone significant changes due to global warming, which have important consequences with regard to the global water cycle [1]. It has been in the principal cryospheric region of northwestern China, including Xinjiang (XJ), Gansu (GS), Ningxia (NX), and Qinghai (QH) provinces and Tibet, where these changes in the water cycle have been most pronounced [1]. Therefore, the study of the spatiotemporal variation of water storage change (WSC) in these areas has importance regarding the understanding of the water cycle, climate change, agricultural production, and natural disasters.

In northwestern China water sources are mainly melting snow and ice and groundwater recharge. For example, mountain runoff is the only source that recharges surface and subsurface runoff in the basin of the Xinjiang Plain [2]. The socioeconomic development of northwestern China has led to increased water demand, overexploitation of water resources, and consequent water shortages that have become a major restriction to the ongoing development of the region [3, 4]. WSC is an important measure that reflects the changes of regional water resources, which have become one of the most important issues regarding socioeconomic development in China.

It is difficult to estimate WSC using traditional long-term water balance calculations and estimations of the water resources of a watershed/region, and, therefore, it is often considered as zero on the annual scale. However, ignoring WSC could have significant impact on the estimated elements of the water balance under climate change scenarios. In estimating the water balance in the source region of the Yellow River, Xu et al. [5] showed that average evaporation calculated by the Sib2 model is 978.2 mm, but when considering WSC, the average evaporation is 507.4 mm. These results suggest that an assessment of WSC is crucial for understanding the changes of the water balance in northwestern China.

WSC can be derived by calculating the difference between precipitation, evapotranspiration, surface runoff, and soil exchange of moisture and groundwater in the water cycle [6]. However, it is difficult to obtain an accurate estimation of the WSC of a region/watershed because of the lack of observational data and the high cost of the measurement methods. Some researchers have used satellite-derived gravity data to monitor long-term WSC and they have made considerable advances [79]. The advantage of satellite-based methods is their ability to reflect uniform region/basin-scale observations of water storage variations. Traditional remote sensing satellites can be used to detect the soil moisture of the land surface to depths of only tens of centimeters, and the spatial distribution of field stations for verification purposes is poor. However, to a certain extent, the Gravity Recovery and Climate Experiment (GRACE) satellite can compensate for such disadvantages, and it offers a new opportunity for the quantitative study of regional land WSC.

The variation in the gravity field detected by GRACE is mainly due to land WSC, ice masses (e.g., polar ice sheets, alpine glaciers ice cap), and other geophysical signals (e.g., postglacial rebound and deformation caused by earthquakes). Interference signals in the GRACE data are caused by the atmosphere and ocean mass variations, which can be removed using numerical simulation [10], and, therefore, following the removal of the geophysical signals, the WSC can be determined. Results show that the GRACE satellite can monitor the water equivalent change of the ground to about 0.9 cm and the precision can reach 1–1.5 cm [11, 12]. Some researchers have compared the results retrieved by GRACE with those calculated using global land surface process models. The two methods produce results that are similar in most areas, but there are also discrepancies related to difficulties in obtaining precise input data in some areas. GRACE data can overcome such problems and, thus, can provide a reference for the improvement of the accuracy of global land surface process models [1316]. Strassberg et al. [17] used GRACE data to estimate the WSC over about 450,000 km2 of the American High Plains region and found good correlation () between the WSC retrieved by GRACE and the measured soil moisture and groundwater data. In addition, WSC retrieved by GRACE has been used in cryospheric-related research on permafrost activities and WSC in the Arctic. The results showed that WSC increased at a rate of 1.15 km3/a on the coast of Alaska, whereas it showed a rate of decline of 7.44 km3/a in the Yukon basin. Changes of the permafrost active layer might be an important reason behind the variation of WSC, and the results of the study suggest that GRACE satellite data could be used to analyze the variation of WSC in permafrost regions [18]. Chinese research on this subject has concentrated mainly on large areas, such as the Chinese continental region, Xinjiang and the surrounding mountainous area, and northern, northeastern, and southern regions of China [4, 15, 16, 19]. Generally, these researches have involved regional quantitative descriptions of WSC, and they have established that different regions exhibit different variational features. Some studies have used GRACE satellite data to analyze the magnitude of WSC in the Three Gorges Reservoir area and the Yangtze River and Heihe River basins. The amplitude of annual variation of WSC in the Yangtze River basin was up to 3.4 cm equivalent water height with the maximum value appearing in spring and early autumn. The root-mean-square error (RMSE) of the variation of WSC of the Three Gorges Reservoir was 1.3 cm, for the amplitude of the annual variation of 0.1 cm. In addition, the WSC of the Heihe River basin increased yearly from 2002 to 2008 [2024]. These studies reveal that applications of GRACE data are lacking in depth and breadth, although they do indicate that GRACE satellite data could enhance the understanding of WSC, and provide a basis for better predictions of drought and changes in global water resources due to climate change. However, a comprehensive study of the WSC of the cryosphere and of its causes in northwestern China has not been reported previously. In this study, gravity information derived from GRACE satellite observations was used to retrieve monthly WSC, which was compared with results simulated by the Noah model to analyze WSC in northwestern China from 2003 to 2010. This work expands the spatiotemporal scale of the analysis of WSC in northwestern China, by considering the water resources of five provinces (Xinjiang, Gansu, Ningxia, Tibet, and Qinghai) and the Qinghai–Tibet Plateau and by discussing the effects on WSC of glaciers and the underlying permafrost.

2. Study Area

The study area incorporates Tibet and Gansu (GS), Qinghai (QH), Ningxia (NX), and Xinjiang (XJ) provinces with an area of about 2,808,000 km2 (Figure 1). The climate is mainly influenced by monsoon and westerlies, and precipitation declines from the east to the west and from the south to the north. Annual precipitation is about 400 mm in the east, which decreases to 200 mm in the west or 50 mm below. The main topographic features of the area are plateaus, basins, and mountains, including the Loess and Qinghai–Tibet Plateaus, the Tarim, Zhungeer, and Qaidam basins, and the Tianshang, Qianlian, Karakorum, and Kunlun Mountains. The Qinghai–Tibet Plateau, sometimes referred to as the world’s “third pole,” is the largest glacial resource in China. The total ice volume is 2798 km3, which accounts for 50% of the ice reserves. The ice covers about 32320 km2, which accounts for 54.4% of the entire area. There are 22468 glaciers within this region, and the permafrost area extends over about  km2 [25, 26]. There are 18311 glaciers in XJ, accounting for 46.8% of the total number of glaciers in China, covering an area of 24721.93 km2. The total of ice reserves are 2623.4 km3, and the water storage of the glaciers is  m3 [27].

3. Materials and Methods

3.1. Precipitation Data and Spatial Interpolation Method

Precipitation data from 2003 to 2010 were obtained from national meteorological stations within the study area and corrected for the influence of wind [28]. Furthermore, to improve the surface accuracy of precipitation within the study area, the data were interpolated using Kriging with consideration of elevation [29]. This method uses precipitation and site elevation data to perform a linear regression. The spatial was calculated using the formula for linear regression; the residuals (regressed precipitation minus measured precipitation) were interpolated using the Kriging spatial interpolation method, and then the spatial was added to the spatial residuals, which provides the spatial distribution of precipitation:where is the spatial precipitation, is the rasterized precipitation regressed by elevation, and is the spatial residual (units of precipitation: mm).

3.2. Change of Precipitation

Precipitation is the main source of land surface water and changes in precipitation affect evaporation, runoff, and WSC. In this study, the change of precipitation was calculated using the following formula:where is the change of precipitation, is the precipitation in month , is the precipitation in the preceding month, and is the month (from 1 to 12).

3.3. Retrieval Method for WSC Data

This study used Release-04 (R4) Level-3 grids provided by the GRACE Science Team centers. The global grids comprise 1-arc-degree water equivalent mass change complete to degree and order 60 [8, 12]. The GRACE R4 Level-3 (300 km Gauss-smoothing kernel) land and ocean monthly grids from 2003 to 2008 were combined to provide global coverage [20, 23]. Glacial isostatic adjustment or solid mass flow within Earth’s mantle, following the decay of the Pleistocene Laurentide and Cordilleran ice sheets in North America and Scandinavia, was applied to the GRACE grids [30].

3.4. Soil Moisture Simulated by Model

WSC represents a vertically integrated measure of variation of water storage that includes soil moisture, groundwater, surface water, snow and ice, and biomass. A previous study of WSC in the High Plains of the United States showed that soil moisture and groundwater are the primary contributors to terrestrial water storage variability and that the variations due to snow and ice, biomass, and surface water are relatively minor [31]. Hydrological models can be divided into two categories based on global and continental scales: land surface models (LSMs) and water budget models. A land surface process model describes the composition of WSC, which includes groundwater, surface water, water exchange, and storage.

In this study, in order to compare the results retrieved by GRACE, the Noah LSM was used to simulate soil moisture. The Noah model was developed by the Research Applications Laboratory of the National Center for Atmospheric Research (NCAR), and a four-dimensional assimilation system was used to calculate the simulated results. The unsaturated zone of the soil was divided into four layers, the thicknesses of which from the soil surface to the bottom were 0.1, 0.3, 0.6, and 1.0 m. Model input data included precipitation, solar radiation, surface air pressure, humidity, and surface wind speed. The shortwave and longwave radiation data were provided by the Air Force Weather Agency agrometeorological simulation system (AGMET), and all other data were obtained from national meteorological stations. Model output included four layers of soil moisture change with temporal resolution of 1 month and spatial resolution of . It is important to note that model simulations are ineffective in permafrost zones.

3.5. WSC Cycle

To study the extent and period of WSC within study area quantitatively, the least-squares method was used to calculate the cycle and amplitude of regional WSC:where is amplitude, is phase, and is the cycle. Amplitude represents the intensity of WSC and the cycle represents the length of WSC in the time series.

4. Results and Discussion

4.1. Results of WSC Retrieved by GRACE and Simulated by Noah

The results of WSC retrieved by GRACE and simulated by Noah for representative periods of summer (August 2007) and winter (December 2008) are shown in Figure 2. It can be seen that over a large area the two results have similarities. For example, the WSC exhibits a state of loss in the west of XJ and NX, whereas it shows a state of gain in the middle of GS and the Qilian Mountains. However, in some small parts, they exhibit differences. Overall, the WSC retrieved by GRACE demonstrates better continuity than that simulated by Noah.

A comparison of the two results as time series (Figure 3) reveals that both exhibit seasonal change and demonstrate a fluctuating feature (correlation between them is 0.71). However, as can be seen, the range of the values retrieved by GRACE is much larger than that simulated by Noah. One reason for this is that GRACE data can reflect deeper soil water equivalent.

It should be pointed out that although the two results coincide, they have differences in scale, and there are two reasons for this. The first is that the expressions of the two results are different; Noah only reflects changes of soil humidity, whereas GRACE reflects soil, surface water, glaciers, and frozen soil moisture. The other reason is that the distribution of the stations is sparse, which limits the results simulated by the model.

4.2. Temporal and Spatial Patterns of WSC and Precipitation over Northwestern China

WSC was retrieved from GRACE satellite data from 2003 to 2010 over northwestern China, and the distributions of the average WSC and of the average change of precipitation, calculated using formulas (1) and (2), are shown in Figures 4(a) and 4(b), respectively. As can be seen, there is much less change in the south than in the north and much more change in the east than in the west. WSC is in a surplus state in QH and GS, but precipitation does not coincide with the distribution of WSC in this area, which means that the spatial distribution of precipitation is not the only factor driving the spatial distribution of WSC. Precipitation is the main source of WSC, but degradation of the permafrost and increased thickening of the active layer of seasonally frozen soil also increase the storage capacity of the area [5]. WSC is in a state of loss in most parts of XJ, especially in western regions where the largest loss value reaches −175 mm. The WSC shows an intense state of loss in the Himalayas where the value is about −291.1 mm. This is consistent with research that claims there has been glacial mass balance loss in this area and where the retreat has been faster in the middle and west of the region [32, 33]. However, WSC shows gains in the Karakorum region and the Altai Mountains, which could corroborate research that has shown that glaciers in these areas have exhibited mass stability or even expansion. As can be seen, the precipitation has increased in these areas and previous studies have shown that snowfall is less sensitive to warming [34]. The WSC shows serious losses in XJ, except in the north. There are 18311 glaciers in XJ with a total area of 24721.93 km2, but most of the glaciers in XJ are retreating [35, 36], which means that the water resource has decreased. Across the entire study area, WSC is between −291.1 and 248.3 mm. As shown in Figure 4(b), the range of variation of precipitation is −20.9 to 43.6 mm with the distribution showing the characteristic of a decrease from southeast to northwest. Precipitation shows a reduction in most of the region, although it shows an increase in northern XJ, eastern Tibet, and QH. The distributions of WSC and precipitation have spatial differences, but they do show localized consistency (such as southeastern QH).

It can be seen from Figure 5 that the amplitude of variation of WSC across the study area is −36.4 to 55.9 mm. This means that the maximum value of increased water is about 55.9 mm equivalent water height, and the minimum value is −36.4 mm equivalent water height. Monthly precipitation is about 1.1–70.5 mm. Generally, WSC shows a slight upward trend, but the difference between inflow and outflow is increasingly small. WSC and precipitation have obvious seasonal variation and consistent change processes. The study area is located in inland areas of the mid-latitude westerly belt, where sources of water vapor depend mainly on the zonal westerly circulation and the Indian Ocean monsoon [37, 38]. Rainfall is concentrated mainly from May to September, which is when the WSC generally increases and the peak value appears in this period. Precipitation is relatively small from December to March, but evaporation is large, and, therefore, the WSC is at a loss. The characteristics of variation of precipitation have a dominant effect on WSC in the time series. In northwestern China, the WSC increases slightly year by year from 2003 to 2010; the average rate of increase is about 0.61 mm/a, and the total increase in the quantity of water is about  m3.

4.3. Monthly Average WSC from 2003 to 2010

The distributions of monthly average WSC indicate differences during the years for which it was retrieved by GRACE (Figure 6). It can be seen that the WSC is in a state of loss in the Himalayan region from January to June and in December, with particularly serious losses in February and March. However, it is in a state of gain from July to November, although as the months proceed, the state of gain weakens. That pattern is due to the low precipitation and high evaporation during spring and winter and greater precipitation brought in summer by the Indian monsoon, which brings much more water vapor [39, 40]. The XJ region is influenced by the westerly circulation, which brings increased levels of water vapor and a state of gain in western parts from March to July. This state is particularly evident in the Karakorum and Kunlun Mountains; the state of gain moves west to middle. The WSC is in a state of gain from January to March in the Three Rivers source area, because spring snowmelt runoff begins, which results in a of state of loss from April to June. This area is influenced by the rainy season as well as by glacial meltwater supply during July, which means that the WSC is in a state of gain from July to February of the following year.

4.4. Annual Distribution of WSC from 2003 to 2010

We also calculated the annual distribution of WSC (Figure 7). It can be seen from the figure that the WSC is in a state of gain from March to July in XJ but that it is in a state of loss for the rest of the year. The maximum value (about 13.5 mm) appears in July and the minimum value (about −28.2 mm) is in October. The WSC of GS shows a state of loss in January, February, and December, whereas it is in a state of gain for the other months and the maximum value (about 21.6 mm) occurs in September. The WSC is in a state of loss in QH from April to June, principally because the area includes the Three Rivers source area and the runoff of spring snowmelt leads to supplementary losses. In the other months, WSC is in a state of gain and the maximum value (about 44.4 mm) occurs in September. The WSC is irregular in NX, showing a state of gain in March and April and from August to December, but is in a state of loss in the other months; the maximum and minimum values occur in September and December, respectively. The WSC is more consistent in Tibet and on the Qinghai–Tibet Plateau. It is in a state of gain from July to November but in a state of loss in the other months. The maximum value appears in September and the minimum in April.

4.5. Process of WSC in Northwestern China

We calculated the WSC and precipitation during 2003–2010 in XJ, GS, QH, NX, and Tibet and on the Qinghai–Tibet Plateau (Figure 8). It can be seen that regional WSC exhibits strong seasonal fluctuations, but changes in precipitation are not obvious. The WSC increases in GS and QH and on the Qinghai–Tibet Plateau from 2003 to 2010, at rates of 3.63, 9.68, and 2.49 mm/a, respectively. The WSC decreases in XJ, NX, and Tibet at rates of −1.49, −2.54, and −3.05 mm/a, respectively. The total amount of water increased by about  m3 in QH and it decreased by about  m3 in Tibet (Table 1).

The trends of WSC retrieved by GRACE are in line with the normal patterns and we analyzed the cycles using (3). The amplitude of the land surface water resource reflects the intensity of the hydrological cycle, and the cycle represents the duration of the cycle of WSC. As can be seen from Table 2, the smallest intensity of WSC is 3.4 mm in QH, which means that the WSC is in a steady state in this area. The largest amplitude is 11.5 mm in GS, indicating that the change of water resource in this area has greater complexity and uncertainty and that although the WSC of the area has increased overall, the distribution of WSC is unequal. The largest cycle of WSC is in NX, which is about 25 months, and the smallest is in XJ, which is about 18 months. This reflects the greater development of agriculture and industry in XJ, which accelerates the speed of the water cycle. Whether on temporal or spatial scales, the distribution of WSC in northwestern China is shown to be heterogeneous.

4.6. Influence of Glaciers and Permafrost on WSC

The distributions of glaciers, permafrost, and talik (a layer of ground that remains unfrozen year-round) affect the characteristics of surface drainage and groundwater storage changes in northwestern China [5]. The dynamics of permafrost, talik, and the active layer exert significant influences on surface and groundwater hydrology and geomorphology in area A of Figure 9. Area A is a region of numerous and sometimes quite large thermokarst thaw ponds, lakes, and drained thaw-lake basins [39, 40]. Area B has a large number of glaciers that have retreated considerably [36]. Glacial meltwater flows into area B from the northern slopes of the Himalayas. However, between areas A and B, a considerable amount of predominantly continuous permafrost prevents movement of the surface water into the groundwater. Thus, the active layer increases and thickens, which leads to greater infiltration of surface water into the groundwater, which results in increased water storage in area A. The GRACE monthly water equivalent changes capture the changes of the base flow from groundwater and surface water in the permafrost area of the Qinghai–Tibet Plateau. In related research, GRACE-derived WSC has been found to increase in the watershed of the Lena and Yenisei rivers [41]. These watersheds encompass large areas of both continuous and discontinuous permafrost zones and talik.

Glaciers have exhibited mass stability or even expansion in area D of Figure 9 [4244]. The results retrieved by GRACE also reflect that the WSC has increased in this area. Research has indicated that the Karakorum seasonal cycle is dominated by nonmonsoonal winter precipitation, which has uniquely protected it from reductions in annual snowfall under climate warming during the early part of the 21st century [34]. Glacial mass loss has occurred across areas C and A in Figure 9 [35]. Since the early 2000s, the trend of precipitation has not changed much [45, 46], whereas the WSC exhibits large differences (Figure 4(a)). The WSC in the Tianshan Mountains part of area C has decreased, whereas it has increased in the Qilian Mountains part of area A. The total area of the Tianshan Mountains is 270,000 km2, and the area of glaciers is approximately 9225 km2, which accounts for 3.41% of the total area. The area of permafrost is about 63,000 km2, which accounts for 23% of the total area. The area of the Qilian Mountains is 198,000 km2 and the area of glaciers is approximately 1931 km2, which accounts for about 1% of the total area. The area of permafrost is about 100,000 km2, which accounts for 50% of the total area. Both regions rely mainly on the supplement of meltwater runoff from the mountains. However, the WSC is in a state of loss in area C, which shows that the glacial meltwater flows out this area but that less surface water infiltrates into the groundwater. Furthermore, area A has a greater distribution of permafrost, and thus the thickened active layer as well as more ponds might store more water which caused an increase of WSC (Figure 10).

5. Conclusions

(1)The WSC in the northwestern China was retrieved using GRACE data and the results were compared with the simulations of the Noah model. The total amount of WSC in northwestern China was calculated and analyzed and discussed based on the regional WSC of five provinces and the Qinghai–Tibet Plateau. The retrieved results agree well with the simulations by Noah; the correlation coefficient is 0.71 (). And retrieved results reflect the seasonal change, interannual variation, and spatial distribution of the total water resource within the study area.(2)The interannual variability of WSC was found to have decreased in XJ, NX, and Tibet but increased in GS and QH and on the Qinghai–Tibet Plateau. The study area encompasses the principal regions of glaciers and permafrost in China, which have an effect on the annual regional differences in WSC in XJ, GS, QH, NX, Tibet, and the Qinghai–Tibet Plateau. Seasonal changes of WSC were clear, but the seasonal differences and cycles and annual magnitudes of change reflect the imbalance of water resource distribution in northwestern China.(3)Recent warming leads to the degradation of permafrost and retreat of glaciers. Permafrost plays an important role in the water cycle, and the increase and thickening of the active layer can lead to greater infiltration of surface water into the groundwater, which can result in increased water storage and changes in the regional water balance.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This study was supported by the National Key Program for Developing Basic Sciences in China (no. 2013CBA01806), Chinese Academy of Sciences (KJZD-EW-G03-04), National Natural Science Foundation of China (41121001, 41225002, and 41471060), and the Foundation for Young Talents in Cold and Arid Regions Environmental and Engineering Research Institute of CAS (no. Y451191001). The authors would like to thank the editors and the anonymous reviewers for their crucial comments, which improved the quality of this paper.