Next Article in Journal
A Highly Sensitive 3D Resonator Sensor for Fluid Measurement
Previous Article in Journal
TactiGraph: An Asynchronous Graph Neural Network for Contact Angle Prediction Using Neuromorphic Vision-Based Tactile Sensing
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluating the Hydrological Components Contributions to Terrestrial Water Storage Changes in Inner Mongolia with Multiple Datasets

1
China Aero Geophysical Survey and Remote Sensing Center for Natural Resources, China Geological Survey, Beijing 100083, China
2
Key Laboratory of Aerial Geophysics and Remote Sensing Geology, Ministry of Natural Resources, Beijing 100083, China
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(14), 6452; https://doi.org/10.3390/s23146452
Submission received: 9 May 2023 / Revised: 3 July 2023 / Accepted: 10 July 2023 / Published: 17 July 2023
(This article belongs to the Section Remote Sensors)

Abstract

:
In this study, multiple remote sensing data were used to quantitatively evaluate the contributions of surface water, soil moisture and groundwater to terrestrial water storage (TWS) changes in five groundwater resources zones of Inner Mongolia (GW_I, GW_II, GW_III, GW_IV and GW_V), China. The results showed that TWS increased at the rate of 2.14 mm/a for GW_I, while it decreased at the rate of 4.62 mm/a, 5.89 mm/a, 2.79 mm/a and 2.62 mm/a for GW_II, GW_III, GW_IV and GW_V during 2003–2021. Inner Mongolia experienced a widespread soil moisture increase with the rate of 4.17 mm/a, 2.13 mm/a, 1.20 mm/a, 0.25 mm/a and 1.36 mm/a for the five regions, respectively. Significant decreases were detected for regional groundwater storage (GWS) with the rate of 2.21 mm/a, 6.76 mm/a, 6.87 mm/a, 3.01 mm/a, and 4.14 mm/a, respectively. Soil moisture was the major contributor to TWS changes in GW_I, which accounted 58% of the total TWS changes. Groundwater was the greatest contributor to TWS changes in other four regions, especially GWS changes, which accounted for 76% TWS changes in GW_IV. In addition, this study found that the role of surface water was notable for calculating regional GWS changes.

1. Introduction

Global warming accelerates the water cycle processes and directly impacts the spatial and temporal distribution of water resources [1,2,3,4,5]. In addition, the intensive human activities exacerbated the uncertainty of water resources changes [6,7]. Under the pressures of climate changes and anthropogenic activities, arid regions being faced with complex water resource problems, as ecosystems in such regions are usually fragile and limited by water availability [8,9,10,11]. The question of how water resources in arid regions change under the dual influence of climate changes and increasing human activities was a topic of great concern and debate [12,13].
Inner Mongolia Autonomous Region (Inner Mongolia), China, is a typical arid and semiarid region with regional water shortage due to scarce precipitation. Furthermore, water resources overexploitation caused by human activities aggravated the local water shortage. Water resources shortage is the most important factor that restricts the social development, ecological security and agricultural production in this region. Previous studies preliminarily explored the water resources statues [14,15,16,17,18] and conducted the water resources supply and demand balance analyses [19,20,21]. The large regional span of Inner Mongolia makes the significant spatial differences in hydrological conditions. However, few studies revealed the spatial-temporal variations of water resources in Inner Mongolia, due to the lack of in situ monitoring data.
The advances in satellite remote sensing provide a valuable tool for monitoring water resources changes [22,23,24,25], such as precipitation [26,27], evaporation [28,29], lakes water [30,31], soil moisture [32,33,34], groundwater [35,36] and total terrestrial water storage (TWS) [37,38]. Guo et al. [39] evaluated the temporal and spatial changes of TWS in Inner Mongolia during 2003–2021 based on Gravity Recovery and Climate Experiment (GRACE) and GRACE-Follow on (GFO), which revealed the major driving factors. However, challenges still remain in studying the water resource changes in Inner Mongolia.
One challenge is how different hydrological components contribute to the total TWS changes. TWS is defined as all forms of water stored above and underneath the surface of the Earth [40,41]. GRACE and GFO proved to be reliable for monitoring TWS change at large spatial scales [42,43,44,45]. As GRACE/GFO is not able to measure vertical profile of TWS changes, auxiliary data are used to isolate TWS changes into individual hydrological component [46], such as evapotranspiration [47,48], basin discharge [49,50], soil moisture [51,52] and groundwater storage (GWS) [53,54,55]. Machine learning techniques are employed to extract valuable insights from the large volumes of remote sensing multiply data that are generated for mapping terrestrial water changes [56] and monitoring areas with high hydrological changes [57] and might be used for evaluate component contributions as well. Revealing the contributions of different storage components to total TWS changes is crucial for understanding the water cycle processes and water resources management. However, a critical evaluation of the contributions from individual hydrological component to TWS changes is yet to be conducted in Inner Mongolia.
In addition, it is difficult to quantitatively reveal the responses of water resource changes to natural climate changes and anthropogenic activities. For a given region, the TWS changes observed by GRACE/GFO are the comprehensive responses to both climate changes and human activities. Therefore, TWS changes based on GRACE/GFO have certain limitations in distinguishing the natural and human-induced changes [12]. Supplementary data are required to partition the total TWS changes into the natural and human-induced changes. For example, the human-induced TWS changes could be estimated by the differences between total TWS changes based on GRACE/GFO and natural TWS changes which estimate by meteorological factors.
The large regional span of Inner Mongolia makes the significant spatial differences in hydrological conditions. Therefore, in this study, Inner Mongolia was first divided into different subregions according to hydrological conditions. Then, GRACE/GFO data were used to examine the spatial and temporal behavior of water storage in Inner Mongolia. Additional remote sensing data, including satellite images and land surface models, were compiled to study the changes in water cycle elements (precipitation, evapotranspiration, surface water and soil moisture). The objectives of this study were to provide the comprehensive view of the contribution of different hydrological components to total water resources changes in arid regions, and comprehensively decompose the water resources changes caused by climate changes and human activities with the help of multi remote sensing observations.

2. Data and Methods

2.1. Study Area

Inner Mongolia Autonomous Region (referred to as “Inner Mongolia”) is located in the northern frontier of China, stretching diagonally from northeast to southwest in a narrow and long shape. The linear distance is more than 2400 km in east–west from 97°12′ to 126°4′, and 1700 km in north–south from 37°24′ to 53°23′. The study area covers a total area of 1.183 million square kilometers, accounting for 12.3% of China’s land area. Inner Mongolia is adjacent to 8 provinces (Heilongjiang, Jilin, Liaoning, Hebei, Shanxi, Shaanxi, Ningxia and Gansu province), spans the Northeast, North and Northwest China, and shares a 4200 km border with Mongolia and Russia in the north (Figure 1).
Inner Mongolia has a vast territory and complicated landforms with the ecosystem structure composed by forest, grassland, desert and agro-pastoral zone. The unique ecosystem structure laid an important environmental foundation for the rich mineral resource and grain production in Inner Mongolia, forming a resource pattern of “forest in the east and mining in the west, farming in the south and grazing in the north”.
The scarce precipitation and the increasing demand for water resources due to the rapid economic development seriously affected the already fragile ecological environment system, which is mainly manifested as degradation of ecosystem and desertification in agricultural and pastoral zone, reduction in forest area and regional water shortage. The per capita water resource of Inner Mongolia is 2081 m3/person, which is in the middle range, and the average annual water modulus of Inner Mongolia is only 44,000 m3/km2, far lower than the national average level of 295,000 m3/km2.
Under the influence of natural geographical location and topography, there is a great disparity in water resources between the northeast and southwest [18]. There are five groundwater resources zones in Inner Mongolia (Figure 1), respectively, Songhua River Groundwater Zone (GW_I), Liao River Groundwater Zone (GW_II), Yellow River Groundwater Zone (GW_III), Hexi Corridor Groundwater Zone (GW_IV: extreme arid region) and Inner Mongolia High Plain Groundwater Zone (GW_V:).
GW_I and GW_II are in the temperate continental climate, with the annual precipitation of 274–619 mm and 248–514 mm in 2003–2021. The evaporation was 343–461 mm and 312–416 mm in 2003–2020 for GW_I and GW_II. The major land use is forest in GW_I, and cultivated land is concentrated in GW_II. Due to the large amount of water withdrawal, groundwater level decreases in GW_II, and water resources shortage is serious.
The annual precipitation is 155–407 mm in 2003–2021 for GW_III, and the annual evaporation was 238–359 mm in 2003–2020. Most of this area is grassland, the rest is cultivated land. Yellow River transit water and groundwater are the main water source for regional industrial and agricultural development and non-zonal grassland ecosystem.
The land use is bare land in GW_IV. GW_IV is dry and water-poor, with annual precipitation of 58–176 mm in 2003–2021 and annual evaporation of 78–119 mm in 2003–2020. There is a shortage in regional water resources, there are no perennial rivers and mainly groundwater is used.
The land use of GW_V is mainly the grassland, and the rest is cultivated land. The annual precipitation and evaporation in this region were 149–384 in 2003–2021 and 213–315 in 2003–2020, respectively.

2.2. Datasets

An overview of the datasets can be found in Table 1. A detailed description of these data is provided below.
(1)
Gridded TWS anomalies (TWSA) data (relative to the 2004–2009 mean baseline in equivalent water height) were derived from the Gravity Recovery and Climate Experiment (GRACE) and GRACE Follow-On (GFO) mascon solution, provided by University of Texas Center for Space Research (CSR). More details about the CSR RL06 mascon (CSR-M) can be found in [58]. The spatial resolution was 1° × 1°, and temporal resolution was monthly. The GRACE was launched in March 2002 and collected over 15 years of time-variable gravity measurements prior to decommissioning in November 2017. GFO was launched in May 2018, and it obtained over four years of gravity observations to date. Cubic spline interpolation was used to fill in the missing data.
(2)
Soil moisture (SM) and snow water equivalent depth (SWE) were derived from NOAH in Global Land Data Assimilation System (GLDAS), which ingested satellite- and ground-based observational data by land surface models and to generate the land surface states and fluxes [59]. The data spans, temporal resolution and spatial resolution were the same with GRACE mascon solution.
(3)
The monthly gridded precipitation data was obtained from Global Precipitation Measurement (GPM) with spatial resolution of 0.1° × 0.1° [26]. These data were validated by cross-validation and error analysis with gauge-based precipitation. The monthly gridded evapotranspiration was provided by Global Land Evaporation Amsterdam Model (GLEAM) with a spatial resolution of 0.25° × 0.25° [29].
(4)
The lake areas in Inner Mongolia were interpreted by satellite images from Landsat for 2000–2018 and Sentinel-2 for 2019–2021. Modified Normalized Difference Water Index (MNDWI) combined with Normalized Difference Vegetation Index (NDVI) and Enhance Vegetation Index (EVI) were calculated based on above satellite images to acquire information on lake areas. The water area and volume information of Hulun Lake were obtained from Hydroweb.
(5)
Yellow River runoff data were obtained from Yellow River Water Resources Bulletin (Yellow River Conservancy Commission of the Ministry of Water Resources, 2003–2020).

2.3. Methods

2.3.1. The GWS Changes Based on GRACE/GFO

Under the assumption that the changes of water stored in rivers, lakes and reservoirs are negligible, GWS anomalies could be calculated as the following equations:
G W S A = T W S A S M A S W E A  
G W S C = G W S A ( t ) G W S A ( t 1 ) t ( t 1 )  
where GWSA is the GWS anomalies, GWSC is the GWS changes, TWSA is from GRACE/GFO, SMA and SWEA are the anomalies of SM and SWE from GLDAS. To unify the parameters, the SMA and SWEA calculated monthly data from the 2004–2009 mean value.

2.3.2. Contribution of Different Hydrological Components to TWS Changes

The component contribution ratio (CCR) proposed by Kim et al. [60] was used to reflect the contributions of different hydrological components to TWS changes. The ratio was calculated as:
C C R = M A D T V  
M A D = t N | S t S ¯ | N  
where MAD is the mean absolute deviation of the hydrological component, TV is the total variability and is calculated as summation of all components’ MAD. S t is the value of component S (SWE, SM and GWS) at time t and N is the number of months, S ¯ is the average of S t .

2.3.3. Isolation the Natural and Human Induced TWS Changes

In order to quantify the human-induced TWS changes, the total TWS changes from GRACE/GFO and natural TWS changes based on climate water balance were compared.
The total TWS changes based on GRACE/GFO are expressed as Equation (5):
T W S C = T W S A ( t ) T W S A ( t 1 ) t ( t 1 )  
where TWSA is the TWS anomalies from GRACE/GFO.
The natural TWS changes are derived by water balance approach without considering the human activities. Natural TWS changes are deduced by monthly precipitation, evapotranspiration and surface runoff as in Equation (6):
T W S C n a t = P r e E v a p R u n o f f  
where Pre is the precipitation from GPM, Evap is the evapotranspiration from GLEAM, and Runoff the surface runoff. As Inner Mongolia is located in an endorheic arid region, where the surface runoff is limited, the surface runoff was ignored in this method.
The TWS changes based on Equation (5) were combined with TWS changes based on Equation (6) to isolate the human-activities-induced TWS changes (TWSChum).
T W S C h u m = T W S C T W S C n a t  

2.3.4. Statistics Analysis

The correlation analysis can reveal the relationship between two sequences. In this study, the Pearson correlation coefficient (r) was used to quantitative two sequences.
r = i = 1 n ( X i X ¯ ) ( Y i Y ¯ ) i = 1 n ( X i X ¯ ) 2 i = 1 n ( Y i Y ¯ ) 2  
Root mean square error (RMSE) was used to evaluate the influences of surface water changes on local TWS\GWS changes.
R M S E = 1 n i = 1 n ( X i Y i ) 2

3. Results

3.1. The TWS Changes Based on GRACE/GFO

Figure 2 showed the TWSA variations in Inner Mongolia based on GRACE/GFO in the period of 2003–2021. The TWS varied obviously in different groundwater zones. The average TWSA was −25.43 mm, −33.97 mm, −15.83 mm and −13.95 mm for GW_II, GW_III, GW_IV and GW_V, respectively, indicating the deficit of land water storage. I average TWSA was 13.82 mm for GW_I, indicating the surplus land water storage.
The trends of TWS changes were also different for different groundwater zones (Figure 2). In GW_I, the TWS increased with the rate of 2.14 mm/a during 2003–2021, which significantly increased with the rate of 47.22 mm/a in 2018–2021 (Figure 2b). In GW_II, the TWS decreased with the rate of 4.62 mm/a during 2003–2021, decreased with the rate of 5.46 mm/a in 2003–2015 and increased with the rate of 28.69 mm/a in 2018–2021 (Figure 2c). The TWS decreased with the rate of 5.89 mm/a and 2.79 mm/a during 2003–2021 for GW_III and GW_IV (Figure 2d,e). In GW_V, the TWS decreased with the rate of 2.62 mm/a during 2003–2021, decreased with the rate of 5.80 mm/a in 2003–2011 and increased with the rate of 5.73 mm/a in 2018–2021 (Figure 2f).

3.2. The Surface Water Changes in Inner Mongolia

As the river water is limited in Inner Mongolia, this part mainly introduced the variation characteristics of lakes area based on satellite images and snow water based on GLDAS.

3.2.1. The Changes of Lakes Area Based on Satellite Images

To better describe the changes of different-sized lakes, the lakes were categorized into three classes: small (<1 km2), medium (1–10 km2) and large (>10 km2). GW_I had the richest lake water resources, with a total of 107 lakes, of which small, medium and large lakes accounted for 33%, 55% and 12%, respectively. Hulun Lake, the largest lake in Inner Mongolia, was in GW_I with an area of more than 1500 km2. Influenced by Hulun Lake, the area decreased at the rate of 34.67 km2/a in 2003–2012, and increased at the rate of 24.17 km2/a in 2013–2021 for large lakes in GW_I.
A total of 74 lakes were detected across GW_II in 2003–2021, of which small, medium and large lakes accounted for 35%, 61% and 4%, respectively. The lake area was relatively stable for lakes in GW_II. The average areas for small, medium and large lakes were 12.18 km2, 115.70 km2 and 41.58 km2, respectively.
A total of 82 lakes were detected across GW_III in 2003–2021, of which small, medium and large lakes accounted for 20%, 76% and 5%, respectively. The lake area was also stable for lakes in GW_III, with average areas of 11.84 km2, 165.50 km2 and 203.96 km2 for small, medium and large lakes, respectively.
The lake water resources were limited in GW_IV. Only 28 lakes were detected across GW_IV in 2003–2021, 46%, 46% and 7% for small, medium and large lakes, respectively. The average areas were 6.20 km2, 32.03 km2 and 53.36 km2 for small, medium and large lakes, respectively. The lake area increased at the rate of 1.75 km2/a in 2003–2021.
A total of 132 lakes were detected in GW_V in 2003–2021, with 38%, 55% and 8% for small, medium and large lakes, respectively. The average total area was 645.32 km2, with 34.46 km2, 188.89 km2 and 421.98 km2 for small, medium and large lakes, respectively. The area decreased at the rate of 2.67 km2/a for lakes in GW_V in 2003–2021 (Figure 3).

3.2.2. The Changes of Snow Water Based on NOAH

The SWE from NOAH, GLDAS, was used to reveal the changes of snow water. There were significant differences in snow water among different regions. The average monthly\annual snow water equivalent depth (SWE) was 3.19\3.58 mm for Inner Mongolia. The average monthly\annual SWE was 6.67\7.38 mm, 1.06\1.34 mm, 0.35\0.36 mm, 0.13\0.15 mm and 0.91\1.07 mm in GW_I, GW_II, GW_III, GW_IV and GW_V, respectively (Figure 4a). The average monthly\annual snow water equivalent depth anomalies (SWEA) were 1.04\1.75 mm, 0.60\0.88 mm, 0.05\0.06 mm, 0.05\0.07 mm and 0.42\0.58 mm for GW_I, GW_II, GW_III, GW_IV and GW_V, respectively (Figure 4b).

3.3. The SM Changes Based on NOAH

The average SM was 354.55 mm for Inner Mongolia in 2003–2021. The order of average SM was GW_III > GW_I > GW_II > GW_IV > GW_IV, with 382.10 mm, 332.72 mm, 406.34 mm 329.98 mm and 319.89 mm, respectively, while the average soil moisture anomalies (SMSA) were surplus for Inner Mongolia with the order of GW_I > GW_II > GW_V > GW_III > GW_IV (Figure 5).
In GW_I, the SMA increased with the rate of 4.17 mm/a from 2003 to 2021, with two increase stages and two decrease stages (Figure 6b). The SMA increased with a rate of 15.01 mm/a and 38.49 mm/a in 2008–2013 and 2018–2021 and decreased with a rate of 24.91 mm/a and 20.06 mm/a in 2004–2007 and 2014–2017, respectively. In GW_II, the SMA increased with the rate of 2.13 mm/a during 2003–2021 and 17.50 mm/a in 2018–2021 (Figure 6c). In GW_III, the SMA increased with the rate of 1.2 mm/a during 2003–2021 (Figure 6d). In GW_IV, the SMA was stable and only increased with the rate of 0.25 mm/a from 2003 to 2021 (Figure 6e). In GW_V, the SMA increased with the rate of 1.36 mm/a during 2003–2021 and 12.60 mm/a in 2018–2021 (Figure 6f).

3.4. The GWS Changes Based on Water Balance

There were successful studies calculating the GWS change in arid regions based on GRACE\GFO data [55,61,62,63,64,65]. Like most studies, the GWS changes in Inner Mongolia were calculated. The GWS decreased with a rate of 2.21 mm/a, 6.76 mm/a, 6.87 mm/a, 3.01 mm/a and 4.14 mm/a for GW_I, GW_II, GW_III, GW_IV and GW_V, respectively, for 2003–2021 (Figure 7).

3.5. The Spatial Changes of Water Resources in Inner Mongolia

There were great spatial differences for water resources changes in Inner Mongolia. TWS changed at the rate of −9.59–11.39 mm/a, decreased in most regions and only increased in the northeast region. Soil moisture changed with small rate of −1.21–8.70 mm/a, increased in most of other regions and only decreased in the western area. Compared with TWS and soil moisture, the change rates of lake water and snow cover were small and could be ignored. Therefore, the GWS changed at the rate of −13.27–5.32 mm/a, almost decreased in whole Inner Mongolia and only increased in a small region of the northeast regions (Figure 8).

4. Discussion

4.1. Contribution of Hydrological Components to TWS Changes

Like most studies based on GRACE\GFO data, we decomposed the TWS changes into three key vertical water storage components: snow water, soil moisture and groundwater, and we ignored the changes in surface runoff and lake water due to the limited surface water resources in arid and semiarid regions. Understanding the contributions of above different hydrological components to the total TWS changes is crucial for investigating how the changes in individual component can potentially affect the availability and utilization of water resources.
Firstly, the influences of snow water, soil moisture and groundwater to TWS changes were qualitatively evaluated by correlation analysis. Table 2 showed the Pearson coefficients of TWS changes with snow water, soil moisture changes and groundwater storage changes. In GW_I, TWS changes showed high correlation with soil moisture changes (R2 = 0.91). In GW_II, TWS changes were positively correlated with both soil moisture changes and groundwater changes (R2 = 0.56 and R2 = 0.69, respectively). TWS changes showed high correlation with groundwater changes in GW_III, GW_IV and GW_V, and the coefficients were 0.66, 0.79 and 0.65, respectively.
Secondly, the component contribution ratio (CCR) was used to quantitatively calculate the contributions of snow water, soil moisture and groundwater to TWS changes in different groundwater regions. Figure 9 showed the MAD and CCR in different groundwater regions based on Equations (3) and (4). The CCR indicated that the soil moisture was the major contributor to TWS changes in GW_I, followed by groundwater. Soil moisture and groundwater accounted for 58% and 32% of total TWS changes in GW_I, while groundwater was the major contributor to TWS changes in other four groundwater regions. The groundwater contributed 76% of total TWS changes in GW_IV, and soil moisture only accounted 23%. Groundwater contributed 62% and 57% of the total TWS changes in GW_II and GW_III, respectively, followed by soil moisture of 36% and 41%. In addition, changes in snow cover in the northeast also contributed to 10% of the total regional TWS changes.
In order to reveal the spatial distribution of contribution from groundwater to TWS changes, an index ® was defined by the change rate of GWS changes with the rate of TWS changes. From Figure 10, there were five conditions for this index:
(1)
The ratio was negative in the northeast region of Inner Mongolia. In this area, the TWS and soil moisture both showed an increase trend, while GWS showed a decrease trend. Soil moisture controlled the total TWS changes, and groundwater contributed negatively to the TWS changes.
(2)
The ratio was 0–0.3 in a small part of GW_I. In this area, TWS, soil moisture and groundwater all increased. In this condition, soil moisture and groundwater both contributed positively to the changes of total water storage, but the soil moisture was also the major contributor.
(3)
The ratio was between 0.3 and 0.6 in the area near Hulun Lake and Hei River. In these areas, TWS, soil moisture and groundwater all decreased, and soil moisture and groundwater contributed to the TWS changes in similar weight.
(4)
The ratio was between 0.6 and 1 in the area near the Hei River and eastern of GW_IV. In these regions, TWS, soil moisture and groundwater all decreased, with groundwater decreasing at a faster rate. TWS changes were mainly controlled by groundwater.
(5)
The ratio was higher than 1 in the rest areas, where the soil moisture increased and groundwater and TWS decreased. In these areas, groundwater was the major contributor of TWS changes.

4.2. The Role of Surface Runoff and Lake Water on Local Water Resources

Surface water is a key hydrological component, an important factor influencing the water cycle and ecological environment [66]. Kim, Yeh, Oki and Kanae [60] evaluated the role of a river in total TWS changes over global basins and concluded that neglecting river storage may lead to a mismatch in the amplitude and phase of TWS seasonal variations. In addition, studies showed that most of the lakes in semiarid and arid northern China experienced an obviously reduction in past decades [67,68,69]. Tao et al. [70] revealed that Inner Mongolia experienced significant lake shrinkage during the past several decades due to the unsustainable mining boom and agricultural irrigation, and the total water area of the lakes decreased from 4160.2 km2 in the late 1980s to 2900.6 km2 in 2010 with a decrease of 30.3%.
Despite this, most of the previous studies of GRACE/GFO hydrology applications did not consider the roles of surface water, and the impacts of surface runoff and lakes on TWS\GWS changes in arid regions did not receive enough attention relative to soil moisture and groundwater. Therefore, we raised a question that the lake water or surface runoff may have significance to the local water resources changes. To test this idea, we discussed the influences of Yellow River runoff changes and Hulun lake storage changes on local TWS\GWS changes in the respective groundwater resources zone.

4.2.1. The Influence of Yellow River on TWS and GWS Changes in GW_III

Yellow River, the second longest river in China, flows through Inner Mongolia in GW_III. The runoff of Yellow River showed an increase trend from 2003 to 2021. The contribution of Yellow River runoff to the TWS in GW_III was calculated by the difference between the inflow and outflow of Yellow River. For consistency, the Yellow River runoff was processed as the GRACE/GFO data, i.e., anomalies were calculated relative to 2004–2009.
The average Yellow River runoff anomalies in 2003–2021 was 0.94 km3, which was obviously higher than that of the snow water in GW_III (0.11 km3) (Figure 11a). This indicated that the Yellow River runoff had higher weight compared with the snow water. However, the Yellow River runoff was not enough to change the TWS trend in GW_III. The Yellow River runoff anomalies increased at a rate of 0.15 km3/a, while the annual TWS in GW_III decreased at a rate of −0.91 km3/a.
GWS changes in GW_III were underestimated if the runoff of Yellow River was not considered (Figure 11b). The annual average GWSA in GW_III was −6.19 km3 ignoring the Yellow River and −7.13 km3 considering the Yellow River. The root mean square error (RMSE) was 1.45 km3 between GWSA ignoring the runoff and GWSA considering the runoff. In addition, the GWS in GW_III decreased at a rate of 1.25 mm/a considering the runoff, which was higher than that which ignored the runoff (1.1 mm/a). Therefore, the runoff of Yellow River should not be neglected when analyzing the GWS changes in GW_III.

4.2.2. The Influence of Hulun Lake on TWS and GWS Changes in GW_I

Hulun Lake, which is the largest lake in Inner Mongolia and the fifth largest lake in China, is located in GW_I. The water area and water level of Hulun Lake experienced obviously changed [71,72]. Based on the data from Hydroweb, the lake area decreased from more than 2000 km2 to around 1500 km2 in 1992–2012 and increased to more than 2100 km2 in 2022. The lake height decreased from around 545 m to 540 m in 1992–2012 and increased back to 545 m in 2022.
The lake water volume anomalies were calculated by abstracting the average value of 2004–2009. The average of lake water volume anomalies in 2003–2021 was 1.57 km3, less than that of snow water (2.48 km3) (Figure 11c). This was due to the rich snow resources in GW_I.
GWS changes were underestimated if the lake water volume was not considered (Figure 11d). The annual average GWSA was −7.85 km3 ignoring the lake water volume, and −9.43 km3 considering the lake water volume. The RMSE was 2.64 km3 between GWSA ignoring the lake water and GWSA considering the lake water. Therefore, Hulun Lake should not be neglected when analyzing the GWS changes in GW_I. The GWS considering the lake water in GW_I decreased at a higher rate (0.96 mm/a) than that ignoring the lake water (0.74 mm/a). Therefore, the runoff of Yellow River should not be neglected when analyzing the GWS changes in GW_III.
According to the above analyses, the surface water, which did not receive much attention in previous studies, is an important water storage component. Additionally, the influences of surface runoff and lake water to local GWS changes are not negligible, based on water balance theory.

4.3. Natural and Human Activities Induced TWS Changes

GRACE/GFO captures the total TWS changes caused by both natural and anthropogenic drivers. Natural TWS changes are estimated by climate water balance. Human-activities-induced TWS changes could be estimated by computing the difference between the total TWS changes and natural TWS changes. Figure 12 showed the comparison between monthly natural-induced, human-induced and integrated TWS changes.
As can be seen from Figure 12, a good agreement can be seen between total TWS changes and natural TWS changes in terms of the seasonal distribution for GW_I and GW_II. However, significant discrepancies were also apparent in terms of the magnitude. For example, the variation range of natural TWS in GW_I was −21.54–32.01 mm, which was greater than total TWS (−6.65–5.99 mm). This indicated that human activities significantly reduced the seasonality of TWS caused by meteorological factors, although TWS changes in GW_I tended to increase due to the natural precipitation replenishment.
For GW_III and GW_IV, there were obviously differences in the seasonality between total TWS changes and natural TWS changes. For example, the natural TWS showed the most surplus in September, but TWS changes based on GRACE/GFO was deficit in September. This indicated human activities changed the seasonality of TWS changes in these regions.
It was interesting that the natural TWS showed more serious losses than GRACE/GFO-based TWS in GW_V, −17.62 mm, and −1.01 mm, respectively. In spring (March –May), both natural TWS changes and total TWS changes were deficit, −27.47 mm, and −10.10 mm, respectively. In summer (June–August), both natural TWS changes and total TWS changes were surplus, 5.67 mm and 18.61 mm, respectively. The above phenomena indicated that human activities played a positive role in total TWS changes in this region in spring and summer. The average natural TWS change was 0.30 mm in autumn, while the total TWS change was −11.72 mm, indicating human activities played a negative role in total TWS changes. The average natural TWS changes and total TWS changes were similar in winter, which indicated that natural climate factors dominated the TWS of this region.
Figure 13 compared the total TWS changes with natural TWS changes and human-activities-induced TWS changes in annual scale. Table 3 showed the correlation coefficient (R2) and RMSE between total TWS changes with natural TWS changes and human-activities-induced TWS changes in annual scale. Based on the above results, the total TWS changes were dominated by climate factors for GW_I and GW_II. Human activities explained the major TWS changes in IV. For the GW_III and GW_V, the situations were more complicated, it was hard to disentangle the contribution from climate factors and human activities.
Some studies showed that the changes of land water storage changes based on climate water balance were highly consistent with that based on GRACE [73,74]. However, the situation was different in this study, and this was not surprising for the discernable differences. We speculated these differences may be induced by the large water consumption caused by human-associated activities would overestimate the natural TWS changes. Without enough surface water, Inner Mongolia has little choice but to increase the reliance on groundwater to replenish the water for agricultural irrigation and coal mining [39]. Additionally, water storage changes in Inner Mongolia during the study period followed the typical pattern of increased groundwater abstraction in arid and semiarid regions [75]. Therefore, the deterioration of groundwater storage is expected to continue in the following decades not only because of the change climate but also due to increasing exploitation.

5. Conclusions

This study revealed the temporal and spatial changes of total TWS and their individual components in Inner Mongolia, China, using multiple remote sensing data for the period of 2003–2021, evaluated the contribution of hydrological components to TWS changes and isolated the human-induced TWS from total TWS changes.
Soil moisture and groundwater are the major contributors of TWS changes in Inner Mongolia. To be specific, soil moisture primarily contributes to TWS changes in Songhua River Groundwater Zone (GW_I), soil moisture and groundwater both contribute to Liao River Groundwater Zone (GW_II), and groundwater mainly composes the TWS changes in Yellow River Groundwater Zone (GW_III), Hexi Corridor Groundwater Zone (GW_IV) and Inner Mongolia High Plain Groundwater Zone (GW_V).
The contributions of surface runoff and lake water to TWS changes are usually ignored in arid and semi-arid areas. However, our analyses showed the runoff and lake water cannot be considered as an invariant constant which cannot change the major result about the groundwater changes.
In relation to natural- and human-induced changes in the terrestrial water cycle, the natural climate condition controls the water resources changes in GW_I and GW_II, while human activities influenced the water resources changes in GW_IV. Human activities changed the seasonality of TWS changes in GW_III and GW_V.
Due to the lack of in situ data available for validation, it is crucial to recognize that the analyses and resulting conclusions were subject to inherent uncertainty. However, we believe that this combination of multiple remote sense data provides a very valuable alternative for understanding hydrologic changes occurring in data-scarce regions.

Author Contributions

Y.G. and N.X. processed the data and conducted the analysis and drafted the paper; F.G., J.B. and B.Y. worked collaboratively on structuring and reviewing the text, as well as discussing the results. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Remote Sensing Quantitative Survey and Monitoring of Water Cycle Elements and Natural Resources in River Basin (DD20230417), China Geological Survey.

Institutional Review Board Statement

Not applicable for studies not involving humans or animals.

Informed Consent Statement

Not applicable.

Data Availability Statement

The GRACE and GRACE Follow-On data were collected from http://www2.csr.utexas.edu/GRACE/RL06_mascons.html (accessed on 13 October 2022). The GLDAS data were collected from http://disc.sci.gsfc.nasa.gov/serVICes/grads-gds/gldas/, accessed on 13 October 2022). The gridded precipitation data were obtained from Global Precipitation Measurement (GPM) (https://gpm.com.hk/, accessed on 13 October 2022). The evapotranspiration data were from Global Land Evaporation Amsterdam Model (http://www.gleam.eu, accessed on 13 October 2022).

Acknowledgments

This research was supported by the Remote Sensing Quantitative Survey and Monitoring of Water Cycle Elements and Natural Resources in River Basin (DD20230417), China Geological Survey. The financial support is highly appreciated. We would like to express thanks for the GRACE/FO Mascon data supported by the CSR and JPL, GLDAS CLSM/NOAH/VIC V2.1 from the Goddard Earth Sciences Data and Information Services Center, NASA. We are also thankful for the insightful comments raised by anonymous reviewers and the associate editor, who helped to improve this paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Barnett, T.P.; Adam, J.C.; Lettenmaier, D.P. Potential impacts of a warming climate on water availability in snow-dominated regions. Nature 2005, 438, 303–309. [Google Scholar] [CrossRef] [PubMed]
  2. Scherler, D.; Bookhagen, B.; Strecker, M.R. Spatially variable response of Himalayan glaciers to climate change affected by debris cover. Nat. Geosci. 2011, 4, 156–159. [Google Scholar] [CrossRef]
  3. Andresen, C.S.; Straneo, F.; Ribergaard, M.H.; Bjørk, A.A.; Andersen, T.J.; Kuijpers, A.; Nørgaard-Pedersen, N.; Kjær, K.H.; Schjøth, F.; Weckström, K.; et al. Rapid response of Helheim Glacier in Greenland to climate variability over the past century. Nat. Geosci. 2011, 5, 37–41. [Google Scholar] [CrossRef] [Green Version]
  4. Kabel, K.; Moros, M.; Porsche, C.; Neumann, T.; Adolphi, F.; Andersen, T.J.; Siegel, H.; Gerth, M.; Leipe, T.; Jansen, E.; et al. Impact of climate change on the Baltic Sea ecosystem over the past 1000 years. Nat. Clim. Change 2012, 2, 871–874. [Google Scholar] [CrossRef]
  5. Taylor, R.G.; Scanlon, B.; Döll, P.; Rodell, M.; van Beek, R.; Wada, Y.; Longuevergne, L.; Leblanc, M.; Famiglietti, J.S.; Edmunds, M.; et al. Ground water and climate change. Nat. Clim. Change 2012, 3, 322–329. [Google Scholar] [CrossRef] [Green Version]
  6. Vorosmarty, C.J.; Green, P.; Salisbury, J.; Lammers, R.B. Global Water Resources Vulnerability from Climate Change and Population Growth. Science 2000, 289, 284–288. [Google Scholar] [CrossRef] [Green Version]
  7. Pritchard, H.D. Asia’s shrinking glaciers protect large populations from drought stress. Nature 2019, 569, 649–654. [Google Scholar] [CrossRef]
  8. Wang, Y.-J.; Qin, D.-H. Influence of climate change and human activity on water resources in arid region of Northwest China: An overview. Adv. Clim. Change Res. 2017, 8, 268–278. [Google Scholar] [CrossRef]
  9. Xu, C.; Chen, Y.; Yang, Y.; Hao, X.; Shen, Y. Hydrology and water resources variation and its response to regional climate change in Xinjiang. J. Geogr. Sci. 2010, 20, 599–612. [Google Scholar] [CrossRef]
  10. Chen, Y.; Deng, H.; Li, B.; Li, Z.; Xu, C. Abrupt change of temperature and precipitation extremes in the arid region of Northwest China. Quatern Int. 2014, 336, 35–43. [Google Scholar] [CrossRef]
  11. Li, Z.; Qi, F.; Li, Z.; Yuan, R.; Gui, J.; Lv, Y. Climate background, fact and hydrological effect of multiphase water transformation in cold regions of the Western China: A review. Earth-Sci. Rev. 2019, 190, 33–57. [Google Scholar] [CrossRef]
  12. Felfelani, F.; Wada, Y.; Longuevergne, L.; Pokhrel, Y.N. Natural and human-induced terrestrial water storage change: A global analysis using hydrological models and GRACE. J. Hydrol. 2017, 553, 105–118. [Google Scholar] [CrossRef] [Green Version]
  13. Zhang, Y.; He, B.; Guo, L.; Liu, J.; Xie, X. The relative contributions of precipitation, evapotranspiration, and runoff to terrestrial water storage changes across 168 river basins. J. Hydrol. 2019, 579, 124194. [Google Scholar] [CrossRef]
  14. Zhao, M. Water Resource Sustainable Development and Utilization Strategies of Inner Mongolia. J. Arid. Land Resour. Environ. 2000, 14, 32–36, (In Chinese Abstract). [Google Scholar]
  15. Jiang, Y.F. Analysis of Water Resource Statue under Climate Change in Inner Mongolia. Meteorol. Inn. Mong. 2018, 3, 22–24, (In Chinese Abstract). [Google Scholar]
  16. Meng, H.F. The present situation and countermeasures of water resources in Inner Mongolia. Inn. Mong. Water Conserv. 2006, 3, 68–70, (In Chinese Abstract). [Google Scholar]
  17. Yu, L.L.; Tang, S.N.; Chen, F.; Ding, Y.Y.; Yang, Y.; Liu, Y.Z. The situation and countermeasure analysis of water resources exploitation and utilization in Inner Mongolia Autonomous Region. Water Conserv. Plan. Des. 2019, 7, 16–20, (In Chinese Abstract). [Google Scholar]
  18. Yang, G.Y.; Wang, H.; Li, S.Y.; Duan, N. Water resources conditions of Inner Mongolia Autonomous Region and measures for water development and utilization. China Water Resour. 2017, 23, 22–25, (In Chinese Abstract). [Google Scholar]
  19. Yang, Y.Z.; Zhang, W.K.; Liu, D.W. Water and land balance and water resources carrying capacity in Inner M ongolia. Arid. Land Geogr. 2008, 31, 436–441, (In Chinese Abstract). [Google Scholar]
  20. Gao, H.M.; Hu, X.S.; Li, F.L. Countermeasures of water resource supply and demand in Inner Mongolia in the early 21st century. Inn. Mong. Water Conserv. 2000, 1, 26–31, (In Chinese Abstract). [Google Scholar]
  21. Wang, Y.; Li, X.; Dai, W.J.; Li, T. Supply and Demand Analysis of Water Resources in Inner Mogolia. Yellow River 2012, 34, 42–44, (In Chinese Abstract). [Google Scholar]
  22. Dou, X. A critical review of groundwater utilization and management in China’s inland water shortage areas. Water Policy 2016, 18, 1367–1383. [Google Scholar] [CrossRef]
  23. Pekel, J.F.; Cottam, A.; Gorelick, N.; Belward, A.S. High-resolution mapping of global surface water and its long-term changes. Nature 2016, 540, 418–422. [Google Scholar] [CrossRef] [PubMed]
  24. Ma, R.; Duan, H.; Hu, C.; Feng, X.; Li, A.; Ju, W.; Jiang, J.; Yang, G. A half-century of changes in China’s lakes: Global warming or human influence? Geophys. Res. Lett. 2010, 37, L24106. [Google Scholar] [CrossRef]
  25. Ma, R.; Yang, G.; Duan, H.; Jiang, J.; Wang, S.; Feng, X.; Li, A.; Kong, F.; Xue, B.; Wu, J.; et al. China’s lakes at present: Number, area and spatial distribution. Sci. China Earth Sci. 2011, 54, 283–289. [Google Scholar] [CrossRef]
  26. Huffman, G.J.; Stocker, E.F.; Bolvin, D.T.; Nelkin, E.J.; Tan, J. GPM IMERG Final Precipitation L3 Half Hourly 0.1 Degree × 0.1 Degree V06; Goddard Earth Sciences Data and Information Services Center (GES DISC): Greenbelt, MD, USA, 2019.
  27. Zuo, J.; Xu, J.; Chen, Y.; Wang, C. Downscaling Precipitation in the Data-Scarce Inland River Basin of Northwest China Based on Earth System Data Products. Atmosphere 2019, 10, 613. [Google Scholar] [CrossRef] [Green Version]
  28. Miralles, D.G.; Holmes, T.R.H.; De Jeu, R.A.M.; Gash, J.H.; Meesters, A.G.C.A.; Dolman, A.J. Global land-surface evaporation estimated from satellite-based observations. Hydrol. Earth Syst. Sci. 2011, 15, 453–469. [Google Scholar] [CrossRef] [Green Version]
  29. Martens, B.; Miralles, D.G.; Lievens, H.; van der Schalie, R.; de Jeu, R.A.M.; Fernández-Prieto, D.; Beck, H.E.; Dorigo, W.A.; Verhoest, N.E.C. GLEAM v3: Satellite-based land evaporation and root-zone soil moisture. Geosci. Model Dev. 2017, 10, 1903–1925. [Google Scholar] [CrossRef] [Green Version]
  30. Luo, S.; Song, C.; Ke, L.; Zhan, P.; Fan, C.; Liu, K.; Chen, T.; Wang, J.; Zhu, J. Satellite Laser Altimetry Reveals a Net Water Mass Gain in Global Lakes With Spatial Heterogeneity in the Early 21st Century. Geophys. Res. Lett. 2022, 49, GL096676. [Google Scholar] [CrossRef]
  31. Wan, W.; Long, D.; Hong, Y.; Ma, Y.; Yuan, Y.; Xiao, P.; Duan, H.; Han, Z.; Gu, X. A lake data set for the Tibetan Plateau from the 1960s, 2005, and 2014. Sci. Data 2016, 3, 160039. [Google Scholar] [CrossRef] [Green Version]
  32. Tangdamrongsub, N.; Han, S.-C.; Yeo, I.-Y.; Dong, J.; Steele-Dunne, S.C.; Willgoose, G.; Walker, J.P. Multivariate data assimilation of GRACE, SMOS, SMAP measurements for improved regional soil moisture and groundwater storage estimates. Adv. Water Resour. 2020, 135, 103477. [Google Scholar] [CrossRef]
  33. Girotto, M.; Reichle, R.H.; Rodell, M.; Liu, Q.; Mahanama, S.; De Lannoy, G.J.M. Multi-sensor assimilation of SMOS brightness temperature and GRACE terrestrial water storage observations for soil moisture and shallow groundwater estimation. Remote Sens. Environ. 2019, 227, 12–27. [Google Scholar] [CrossRef]
  34. Jung, H.C.; Getirana, A.; Arsenault, K.R.; Kumar, S.; Maigary, I. Improving surface soil moisture estimates in West Africa through GRACE data assimilation. J. Hydrol. 2019, 575, 192–201. [Google Scholar] [CrossRef]
  35. Frappart, F.; Ramillien, G. Monitoring Groundwater Storage Changes Using the Gravity Recovery and Climate Experiment (GRACE) Satellite Mission: A Review. Remote Sens. 2018, 10, 829. [Google Scholar] [CrossRef] [Green Version]
  36. Podgorski, J.; Berg, M. Global threat of arsenic in groundwater. Science 2020, 368, 845–850. [Google Scholar] [CrossRef]
  37. Chen, Z.; Jiang, W.; Wu, J.; Chen, K.; Deng, Y.; Jia, K.; Mo, X. Detection of the spatial patterns of water storage variation over China in recent 70 years. Sci. Rep. 2017, 7, 6423. [Google Scholar] [CrossRef] [Green Version]
  38. Wang, J.; Song, C.; Reager, J.T.; Yao, F.; Famiglietti, J.S.; Sheng, Y.; MacDonald, G.M.; Brun, F.; Schmied, H.M.; Marston, R.A.; et al. Recent global decline in endorheic basin water storages. Nat. Geosci. 2018, 11, 926–932. [Google Scholar] [CrossRef] [PubMed]
  39. Guo, Y.; Gan, F.; Yan, B.; Bai, J.; Xing, N.; Zhuo, Y. Evaluation of Terrestrial Water Storage Changes and Major Driving Factors Analysis in Inner Mongolia, China. Sensors 2022, 22, 9665. [Google Scholar] [CrossRef]
  40. Huntington, J.L.; Niswonger, R.G. Role of surface-water and groundwater interactions on projected summertime streamflow in snow dominated regions: An integrated modeling approach. Water Resour. Res. 2012, 48, W11524. [Google Scholar] [CrossRef] [Green Version]
  41. Tangdamrongsub, N.; Han, S.-C.; Jasinski, M.F.; Šprlák, M. Quantifying water storage change and land subsidence induced by reservoir impoundment using GRACE, Landsat, and GPS data. Remote Sens. Environ. 2019, 233, 111385. [Google Scholar] [CrossRef]
  42. Han, S.-C.; Shum, C.K.; Bevis, M.; Ji, C.; Kuo, C.-Y. Crustal Dilatation Observed by GRACE After the 2004 Sumatra-Andaman Earthquake. Science 2006, 313, 658–662. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Tapley, B.D.; Watkins, M.M.; Flechtner, F.; Reigber, C.; Bettadpur, S.; Rodell, M.; Sasgen, I.; Famiglietti, J.S.; Landerer, F.W.; Chambers, D.P.; et al. Contributions of GRACE to understanding climate change. Nat. Clim. Change 2019, 5, 358–369. [Google Scholar] [CrossRef]
  44. Tu, M.Z.; Liu, Z.F.; He, C.Y.; Ren, Q.; Lu, W.L. Research progress of groundwater storage changes monitoring in China based on GRACE satellite data. Adv. Earth Sci. 2020, 35, 643–656, (In Chinese Abstract). [Google Scholar]
  45. Chen, J.; Tapley, B.; Tamisiea, M.E.; Save, H.; Wilson, C.; Bettadpur, S.; Seo, K.W. Error Assessment of GRACE and GRACE Follow-On Mass Change. J. Geophys. Res. Solid Earth 2021, 126, JB022124. [Google Scholar] [CrossRef]
  46. Tapley, B.D.; Bettadpur, S.; Ries, J.C.; Thompson, P.F.; Watkins, M.M. GRACE Measurements of Mass Variability in the Earth System. Science 2004, 305, 503–505. [Google Scholar] [CrossRef] [Green Version]
  47. Pan, Y.; Zhang, C.; Gong, H.; Yeh, P.J.F.; Shen, Y.; Guo, Y.; Huang, Z.; Li, X. Detection of human-induced evapotranspiration using GRACE satellite observations in the Haihe River basin of China. Geophys. Res. Lett. 2017, 44, 190–199. [Google Scholar] [CrossRef]
  48. Zhong, Y.; Zhong, M.; Mao, Y.; Ji, B. Evaluation of Evapotranspiration for Exorheic Catchments of China during the GRACE Era: From a Water Balance Perspective. Remote Sens. 2020, 12, 511. [Google Scholar] [CrossRef] [Green Version]
  49. Chen, Y.; Fok, H.S.; Ma, Z.; Tenzer, R. Improved Remotely Sensed Total Basin Discharge and Its Seasonal Error Characterization in the Yangtze River Basin. Sensors 2019, 19, 3386. [Google Scholar] [CrossRef] [Green Version]
  50. Xie, J.; Xu, Y.P.; Gao, C.; Xuan, W.; Bai, Z. Total Basin Discharge From GRACE and Water Balance Method for the Yarlung Tsangpo River Basin, Southwestern China. J. Geophys. Res. Atmos. 2019, 124, 7617–7632. [Google Scholar] [CrossRef]
  51. Sadeghi, M.; Gao, L.; Ebtehaj, A.; Wigneron, J.-P.; Crow, W.T.; Reager, J.T.; Warrick, A.W. Retrieving global surface soil moisture from GRACE satellite gravity data. J. Hydrol. 2020, 584, 124717. [Google Scholar] [CrossRef]
  52. Tian, S.; Renzullo, L.J.; van Dijk, A.I.J.M.; Tregoning, P.; Walker, J.P. Global joint assimilation of GRACE and SMOS for improved estimation of root-zone soil moisture and vegetation response. Hydrol. Earth Syst. Sci. 2019, 23, 1067–1081. [Google Scholar] [CrossRef] [Green Version]
  53. Asoka, A.; Gleeson, T.; Wada, Y.; Mishra, V. Relative contribution of monsoon precipitation and pumping to changes in groundwater storage in India. Nat. Geosci. 2017, 10, 109–117. [Google Scholar] [CrossRef] [Green Version]
  54. Seo, J.Y.; Lee, S.-I. Spatio-Temporal Groundwater Drought Monitoring Using Multi-Satellite Data Based on an Artificial Neural Network. Water 2019, 11, 1953. [Google Scholar] [CrossRef] [Green Version]
  55. Ahmed, K.; Shahid, S.; Demirel, M.C.; Nawaz, N.; Khan, N. The changing characteristics of groundwater sustainability in Pakistan from 2002 to 2016. Hydrogeol. J. 2019, 27, 2485–2496. [Google Scholar] [CrossRef]
  56. Bartold, M.; Kluczek, M. A Machine Learning Approach for Mapping Chlorophyll Fluorescence at Inland Wetlands. Remote Sens. 2023, 15, 2392. [Google Scholar] [CrossRef]
  57. Shen, Y.; Zheng, W.; Yin, W.; Xu, A.; Zhu, H.; Yang, S.; Su, K. Inverted Algorithm of Terrestrial Water-Storage Anomalies Based on Machine Learning Combined with Load Model and Its Application in Southwest China. Remote Sens. 2021, 13, 3358. [Google Scholar] [CrossRef]
  58. Save, H.; Bettadpur, S.; Tapley, B.D. High-resolution CSR GRACE RL05 mascons. J. Geophys. Res. Solid Earth 2016, 121, 7547–7569. [Google Scholar] [CrossRef]
  59. Rodell, M.; Houser, P.R.; Jambor, U.; Gottschalck, J.; Mitchell, K.; Meng, C.J.; Arsenault, K.; Cosgrove, B.; Radakovich, J.; Bosilovich, M.; et al. The Global Land Data Assimilation System. Bull. Am. Meteorol. Soc. 2004, 85, 381–394. [Google Scholar] [CrossRef] [Green Version]
  60. Kim, H.; Yeh, P.J.F.; Oki, T.; Kanae, S. Role of rivers in the seasonal variations of terrestrial water storage over global basins. Geophys. Res. Lett. 2009, 36, L17402. [Google Scholar] [CrossRef] [Green Version]
  61. Hu, Z.; Zhou, Q.; Chen, X.; Chen, D.; Li, J.; Guo, M.; Yin, G.; Duan, Z. Groundwater Depletion Estimated from GRACE: A Challenge of Sustainable Development in an Arid Region of Central Asia. Remote Sens. 2019, 11, 1908. [Google Scholar] [CrossRef] [Green Version]
  62. Melati, M.D.; Fleischmann, A.S.; Fan, F.M.; Paiva, R.C.D.; Athayde, G.B. Estimates of groundwater depletion under extreme drought in the Brazilian semi-arid region using GRACE satellite data: Application for a small-scale aquifer. Hydrogeol. J. 2019, 27, 2789–2802. [Google Scholar] [CrossRef]
  63. Wang, S.; Liu, H.; Yu, Y.; Zhao, W.; Yang, Q.; Liu, J. Evaluation of groundwater sustainability in the arid Hexi Corridor of Northwestern China, using GRACE, GLDAS and measured groundwater data products. Sci. Total Environ. 2020, 705, 135829. [Google Scholar] [CrossRef] [PubMed]
  64. Guo, Y.; Gan, F.; Yan, B.; Bai, J.; Wang, F.; Jiang, R.; Xing, N.; Liu, Q. Evaluation of Groundwater Storage Depletion Using GRACE/GRACE Follow-On Data with Land Surface Models and Its Driving Factors in Haihe River Basin, China. Sustainability 2022, 14, 1108. [Google Scholar] [CrossRef]
  65. Rodell, M.; Velicogna, I.; Famiglietti, J.S. Satellite-based estimates of groundwater depletion in India. Nature 2009, 460, 999–1002. [Google Scholar] [CrossRef] [Green Version]
  66. Yamazaki, D.; Trigg, M.A. The dynamics of Earth’s surface water. Nature 2016, 540, 348–349. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  67. Liu, H.; Yin, Y.; Piao, S.; Zhao, F.; Engels, M.; Ciais, P. Disappearing lakes in semiarid Northern China: Drivers and environmental impact. Environ. Sci. Technol. 2013, 47, 12107–12114. [Google Scholar] [CrossRef]
  68. Zhang, G.; Yao, T.; Piao, S.; Bolch, T.; Xie, H.; Chen, D.; Gao, Y.; O’Reilly, C.M.; Shum, C.K.; Yang, K.; et al. Climate change drives extensive and drastically different alpine lake changes on Asia’s high plateaus during the past four decades. Geophys. Res. Lett. 2017, 44, 252–260. [Google Scholar] [CrossRef] [Green Version]
  69. Zhang, G.; Yao, T.; Chen, W.; Zheng, G.; Shum, C.K.; Yang, K.; Piao, S.; Sheng, Y.; Yi, S.; Li, J.; et al. Regional differences of lake evolution across China during 1960s–2015 and its natural and anthropogenic causes. Remote Sens. Environ. 2019, 221, 386–404. [Google Scholar] [CrossRef]
  70. Tao, S.; Fang, J.; Zhao, X.; Zhao, S.; Shen, H.; Hu, H.; Tang, Z.; Wang, Z.; Guo, Q. Rapid loss of lakes on the Mongolian Plateau. PNAS 2015, 112, 2281–2286. [Google Scholar] [CrossRef]
  71. Wu, Y.N.; Pang, Z.G.; Lu, J.X.; Qu, W. Dynamic Changes of Hulun Lake Wetland Area Through Remote Sensing Monitoring and Driving Forces Analysis. J. China Institue Water Resour. Hydropower Res. 2013, 11, 20–26, (In Chinese Abstract). [Google Scholar]
  72. Wang, P.F.; Guo, Y.Y.; Zhou, K.; Zheng, S.F.; Jiang, X.; Wang, S.H. Variation of Lake area of Hulun Lake during 1961—2018 and Its Response to Climate Changes. Res. Environ. Sci. 2021, 34, 792–800, (In Chinese Abstract). [Google Scholar]
  73. Idowu, D.; Zhou, W. Performance Evaluation of a Potential Component of an Early Flood Warning System—A Case Study of the 2012 Flood, Lower Niger River Basin, Nigeria. Remote Sens. 2019, 11, 1970. [Google Scholar] [CrossRef] [Green Version]
  74. Sun, W.; Jin, Y.; Yu, J.; Wang, G.; Xue, B.; Zhao, Y.; Fu, Y.; Shrestha, S. Integrating satellite observations and human water use data to estimate changes in key components of terrestrial water storage in a semi-arid region of North China. Remote Sens. Environ. 2020, 698, 134171. [Google Scholar] [CrossRef] [PubMed]
  75. Voss, K.A.; Famiglietti, J.S.; Lo, M.; Linage, C.; Rodell, M.; Swenson, S.C. Groundwater depletion in the Middle East from GRACE with implications for transboundary water management in the Tigris-Euphrates-Western Iran region. Water Resour. Res. 2013, 49, 904–914. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Location and the land use of Inner Mongolia.
Figure 1. Location and the land use of Inner Mongolia.
Sensors 23 06452 g001
Figure 2. The temporal changes of TWSA in Inner Mongolia and different groundwater zones in 2003–2021. (a) The temporal changes of TWSA in Inner Mongolia in 2003–2021, (b) The temporal changes of TWSA in GW_I in 2003–2021, (c) The temporal changes of TWSA in GW_II in 2003–2021, (d) The temporal changes of TWSA in GW_III in 2003–2021, (e) The temporal changes of TWSA in GW_IV in 2003–2021, (f) The temporal changes of TWSA in GW_V in 2003–2021.
Figure 2. The temporal changes of TWSA in Inner Mongolia and different groundwater zones in 2003–2021. (a) The temporal changes of TWSA in Inner Mongolia in 2003–2021, (b) The temporal changes of TWSA in GW_I in 2003–2021, (c) The temporal changes of TWSA in GW_II in 2003–2021, (d) The temporal changes of TWSA in GW_III in 2003–2021, (e) The temporal changes of TWSA in GW_IV in 2003–2021, (f) The temporal changes of TWSA in GW_V in 2003–2021.
Sensors 23 06452 g002
Figure 3. The number of lakes in different size (a), the temporal changes of lake area in small lakes (b), medium lakes (c) and large lakes (d) in 2003–2021.
Figure 3. The number of lakes in different size (a), the temporal changes of lake area in small lakes (b), medium lakes (c) and large lakes (d) in 2003–2021.
Sensors 23 06452 g003
Figure 4. The temporal changes of snow water equivalent depth (SEW) (a) and snow water equivalent depth anomalies (SEWA) (b) in Inner Mongolia and different groundwater zones in 2003–2021.
Figure 4. The temporal changes of snow water equivalent depth (SEW) (a) and snow water equivalent depth anomalies (SEWA) (b) in Inner Mongolia and different groundwater zones in 2003–2021.
Sensors 23 06452 g004
Figure 5. The box plot of soil moisture (SM) and soil moisture anomalies (SMA) in Inner Mongolia and different groundwater zones.
Figure 5. The box plot of soil moisture (SM) and soil moisture anomalies (SMA) in Inner Mongolia and different groundwater zones.
Sensors 23 06452 g005
Figure 6. The temporal changes of soil moisture (SMSA) in Inner Mongolia and different groundwater zones in 2003–2021. (a) The temporal changes of SMSA in Inner Mongolia in 2003–2021, (b) The temporal changes of SMSA in GW_I in 2003–2021, (c) The temporal changes of SMSA in GW_II in 2003–2021, (d) The temporal changes of SMSA in GW_III in 2003–2021, (e) The temporal changes of SMSA in GW_IV in 2003–2021, (f) The temporal changes of SMSA in GW_V in 2003–2021.
Figure 6. The temporal changes of soil moisture (SMSA) in Inner Mongolia and different groundwater zones in 2003–2021. (a) The temporal changes of SMSA in Inner Mongolia in 2003–2021, (b) The temporal changes of SMSA in GW_I in 2003–2021, (c) The temporal changes of SMSA in GW_II in 2003–2021, (d) The temporal changes of SMSA in GW_III in 2003–2021, (e) The temporal changes of SMSA in GW_IV in 2003–2021, (f) The temporal changes of SMSA in GW_V in 2003–2021.
Sensors 23 06452 g006
Figure 7. The temporal changes of GWSA in Inner Mongolia and different groundwater zones in 2003–2021. (a) The temporal changes of GWSA in Inner Mongolia in 2003–2021, (b) The temporal changes of GWSA in GW_I in 2003–2021, (c) The temporal changes of GWSA in GW_II in 2003–2021, (d) The temporal changes of GWSA in GW_III in 2003–2021, (e) The temporal changes of GWSA in GW_IV in 2003–2021, (f) The temporal changes of GWSA in GW_V in 2003–2021.
Figure 7. The temporal changes of GWSA in Inner Mongolia and different groundwater zones in 2003–2021. (a) The temporal changes of GWSA in Inner Mongolia in 2003–2021, (b) The temporal changes of GWSA in GW_I in 2003–2021, (c) The temporal changes of GWSA in GW_II in 2003–2021, (d) The temporal changes of GWSA in GW_III in 2003–2021, (e) The temporal changes of GWSA in GW_IV in 2003–2021, (f) The temporal changes of GWSA in GW_V in 2003–2021.
Sensors 23 06452 g007
Figure 8. The spatial distribution of the change rate of TWS, SM, GWS in Inner Mongolia.
Figure 8. The spatial distribution of the change rate of TWS, SM, GWS in Inner Mongolia.
Sensors 23 06452 g008aSensors 23 06452 g008b
Figure 9. The mean absolute deviation (MAD) (a) and component contribution ratio (CCR) (b) in different groundwater regions in Inner Mongolia.
Figure 9. The mean absolute deviation (MAD) (a) and component contribution ratio (CCR) (b) in different groundwater regions in Inner Mongolia.
Sensors 23 06452 g009
Figure 10. The spatial distribution of the ratio of GWS change rate with TWS change rate in Inner Mongolia.
Figure 10. The spatial distribution of the ratio of GWS change rate with TWS change rate in Inner Mongolia.
Sensors 23 06452 g010
Figure 11. The changes of Yellow River runoff and snow water (a), the GWSA ignoring Yellow River runoff and considering the runoff in GW_III (b), the changes of Hulun Lake water volume and snow water (c) and GWSA ignoring Hulun lake water and considering Hulun lake water in GW_I (d).
Figure 11. The changes of Yellow River runoff and snow water (a), the GWSA ignoring Yellow River runoff and considering the runoff in GW_III (b), the changes of Hulun Lake water volume and snow water (c) and GWSA ignoring Hulun lake water and considering Hulun lake water in GW_I (d).
Sensors 23 06452 g011
Figure 12. The monthly variations of natural-induced and human-induced TWS changes.
Figure 12. The monthly variations of natural-induced and human-induced TWS changes.
Sensors 23 06452 g012
Figure 13. The relationship between natural TWS changes (TWSC_nat) with total TWS changes from GRACE/GFO data (TWSC) (a) and human-induced TWS changes (TWSC_hum) with TWSC (b).
Figure 13. The relationship between natural TWS changes (TWSC_nat) with total TWS changes from GRACE/GFO data (TWSC) (a) and human-induced TWS changes (TWSC_hum) with TWSC (b).
Sensors 23 06452 g013
Table 1. Information of datasets used in this study.
Table 1. Information of datasets used in this study.
VariablesDataPeriodSpatial ResolutionTemporal ResolutionSource
TWS changeGRACE/GFO2003–20211° × 1°monthlyhttp://www2.csr.utexas.edu (accessed on 13 October 2022)
Soil moistureGLDAS2003–20211° × 1°monthlyhttps://disc.gsfc.nasa.gov/ (accessed on 13 October 2022)
Snow water
PrecipitationGPM2000–20210.1° × 0.1°30 minhttps://gpm.com.hk/ (accessed on 13 October 2022)
gauges2000–2020-monthlyChina Meteorological Data Service Center (accessed on 13 October 2022)
EvaporationGLEAM2003–20210.25° × 0.25°dailyhttp://www.gleam.eu (accessed on 13 October 2022)
Lake areaLandsat 82000–201830 mmonthlyhttps://glovis.usgs.gov/ (accessed on 13 October 2022)
Sentinel 22019–202130 mmonthlyhttps://scihub.copernicus.eu/dhus/#/home (accessed on 13 October 2022)
Hydroweb -dailyhttp://hydroweb.theia-land.fr (accessed on 13 October 2022)
Runoff dataYellow River -annualYellow River Conservancy Commission of the Ministry of Water Resources (accessed on 13 October 2022)
Table 2. The Pearson coefficient of TWS changes with its hydrological components.
Table 2. The Pearson coefficient of TWS changes with its hydrological components.
GW_ITWSSMSSWSGWS
TWS10.906 **0.3820.183
SMS 10.436−0.226
SWS 1−0.419
GWS 1
GW_IITWSSMSSWSGWS
TWS10.556 *0.4180.686 **
SMS 10.438−0.222
SWS 10.050
GWS 1
GW_IIITWSSMSSWSGWS
TWS10.399−0.3770.660 **
SMS 1−0.412−0.426
SWS 1−0.061
GWS 1
GW_IVTWSSMSSWSGWS
TWS10.2910.3310.788 **
SMS 1−0.046−0.352
SWS 10.266
GWS 1
GW_VTWSSMSSWSGWS
TWS10.3450.3330.649 **
SMS 10.229−0.483 *
SWS 10.027
GWS 1
Notes: ** indicates the confidence coefficient is 0.01, and * indicates the confidence coefficient is 0.5.
Table 3. The correlation coefficient (R) and RMSE between total TWS changes with natural TWS changes and human-activities-induced TWS changes in annual scale.
Table 3. The correlation coefficient (R) and RMSE between total TWS changes with natural TWS changes and human-activities-induced TWS changes in annual scale.
GW_IGW_IIGW_IIIGW_IVGW_V
NHNHNHNHNH
R0.830.450.91−0.500.280.59−0.740.880.76−0.53
RMSE5.056.823.395.083.983.762.852.183.323.83
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Guo, Y.; Xing, N.; Gan, F.; Yan, B.; Bai, J. Evaluating the Hydrological Components Contributions to Terrestrial Water Storage Changes in Inner Mongolia with Multiple Datasets. Sensors 2023, 23, 6452. https://doi.org/10.3390/s23146452

AMA Style

Guo Y, Xing N, Gan F, Yan B, Bai J. Evaluating the Hydrological Components Contributions to Terrestrial Water Storage Changes in Inner Mongolia with Multiple Datasets. Sensors. 2023; 23(14):6452. https://doi.org/10.3390/s23146452

Chicago/Turabian Style

Guo, Yi, Naichen Xing, Fuping Gan, Baikun Yan, and Juan Bai. 2023. "Evaluating the Hydrological Components Contributions to Terrestrial Water Storage Changes in Inner Mongolia with Multiple Datasets" Sensors 23, no. 14: 6452. https://doi.org/10.3390/s23146452

APA Style

Guo, Y., Xing, N., Gan, F., Yan, B., & Bai, J. (2023). Evaluating the Hydrological Components Contributions to Terrestrial Water Storage Changes in Inner Mongolia with Multiple Datasets. Sensors, 23(14), 6452. https://doi.org/10.3390/s23146452

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop