Improved Understanding of Groundwater Storage Changes under the Influence of River Basin Governance in Northwestern China Using GRACE Data

Groundwater is crucial for economic development in arid and semiarid areas. The Shiyang River Basin (SRB) has the most prominent water use issues in northwestern China, and overexploited groundwater resources have led to continuous groundwater-level decline. The key governance planning project of the SRB was issued in 2007. This paper synthetically combines remote-sensing data from Gravity Recovery and Climate Experiment (GRACE) data and precipitation, actual evapotranspiration, land use, and in situ groundwater-level data to evaluate groundwater storage variations on a regional scale. Terrestrial water storage anomalies (TWSA) and groundwater storage anomalies (GWSA), in addition to their influencing factors in the SRB since the implementation of the key governance project, are analyzed in order to evaluate the effect of governance. The results show that GRACE-derived GWS variations are consistent with in situ observation data in the basin, with a correlation coefficient of 0.68. The GWS in the SRB had a slow downward trend from 2003 to 2016, and this increased by 0.38 billion m³/year after 2018. As the meteorological data did not change significantly, the changes in water storage are mainly caused by human activities, which are estimated by using the principle of water balance. The decline in GWS in the middle and lower reaches of the SRB has been curbed since 2009 and has gradually rebounded since 2014. GWS decreased by 2.2 mm EWH (equivalent water height) from 2011 to 2016, which was 91% lower than that from 2007 to 2010. The cropland area in the middle and lower reaches of the SRB also stopped increasing after 2011 and gradually decreased after 2014, while the area of natural vegetation gradually increased, indicating that the groundwater level and associated ecology significantly recovered after the implementation of the project.


Introduction
Aquifers are the world's third largest reservoir of water resources after oceans and the cryosphere, and the largest liquid freshwater reservoir. Due to its safety and stability, as compared to surface water, groundwater accounts for approximately one-third of the world's water consumption [1]. It is estimated that groundwater supplies 50% of the world's domestic water, 40% of industrial water and 20% of irrigation water [2]. In recent decades, groundwater levels have plummeted in many parts of the world as a result of extensive groundwater exploitation, with groundwater levels falling to 200 m in some places [3,4]. Due to the extremely slow process of groundwater renewal, overexploited groundwater resources in these areas will not fully recover for a long period of time. Depletion of groundwater resources not only limits the sustainable development of the regional economy, but also increases energy consumption because of the need for more The limited observation data cannot fully provide information on groundwater-level changes before and after the key governance planning project of the SRB. The objective of this paper is to synthetically combine data, including meteorological, hydrological, groundwater-level and GRACE data to evaluate groundwater storage variations at regional scales. Terrestrial water storage anomalies (TWSA) and groundwater storage anomalies (GWSA), in addition to their influencing factors in the SRB since the implementation of the key governance project, are analyzed in order to evaluate the effect of governance. This study will provide a scientific reference for the rational development and utilization of water resources in the SRB and the improvement of ecological protection.

Study Area
The SRB (Figure 1) is located east of the Hexi Corridor in China and north of the Qilian Mountains, between 36~42 degrees north latitude and 101~104 degrees east longitude. The SRB is surrounded by two deserts to the east, west and north, i.e., the Tenggeli and Badain Jaran deserts, which are typical desert oases, with a total area of 41,600 km 2 . The terrain is high in the south and low in the north, tilting from the southwest to the northeast, and land-use classes in the SRB mainly include croplands, barren land and grasslands, accounting for approximately 10%, 49%, and 41% of the area, respectively. The entire basin can be divided into four geomorphic units: the southern Qilian Mountains, central corridor plains, northern hilly areas and desert areas. In the southern Qilian Mountains, which are 2000~5000 m above sea level, the mountains generally trend from the northwest to the southeast. The central corridor plain area is divided into northern and southern basins by the Hanmu Mountain, Hongya Mountain and Alagu Mountain. The southern basin includes the Dajing, Wuwei and Yongchang sub-basins with elevations ranging from 1400 m to 2000 m. The northern basin includes the Minqin subbasin and Jinchuan-Changning sub-basin with elevations of 1300~1400 m. The lowest elevation is in Baitinghai at only 1020 m. The low hilly area in the north has an altitude of less than 2000 m. The SRB is located deep within the hinterland of the mainland and has a continental temperate arid climate. Climatic characteristics include strong solar radiation, sufficient sunshine, a large temperature difference, low precipitation, strong evaporation and dry air. Precipitation in the SRB is mainly concentrated in the Qilian Mountains. Precipitation in mountainous areas is abundant and relatively stable annually. The average annual precipitation is 300~600 mm, which not only provides good growth conditions for forests and grasslands in mountainous areas, but also reliably guarantees water resources in the plain areas. Precipitation in the plain areas is low, with an average annual value of 200 mm. The basin water system mainly originates in the Qilian Mountains. From east to west, there are seven main rivers and many small ditches and rivers, including the Gulang River, Huangyang River, Zamu River, Jinta River, Xiying River, Dongda River and Xida River, which flow from southwest to northeast. Precipitation in the mountainous areas and melting water from alpine ice and snow are the main water supplies for the rivers. The aquifer system in the study area is mainly composed of the Quaternary porous aquifer system, which varies from the unconfined sand layer in the upper part of SRB to being multi-layered in the middle and lower section of the SRB [24]. For many years, the dynamic characteristics of groundwater in the SRB have included the continuous decline in the groundwater level. From the data investigated by Management Bureau of Shiyang River Basin of Gansu, the last 20 years, from 1980 to 2000, show that the groundwater level in the South Wuwei Basin has decreased by 6~7 m, the rate of decline has been 0.31 m/year, and in the Minqin Basin, the decrease has shown an average of 10~12 m, with a decline of 0.57 m/year, and the maximum decrease being 15~16 m. At the beginning of the 21st century, the consumption rate of water resources in the basin reached 109%, and the development and utilization of water resources reached 172%, which is far more than the reasonable carrying capacity of the water resources in the basin. The superposition of resource-based and structural water shortages has led to the deterioration of the ecological environment of the basin, such as the decline in the groundwater level, disappearance of lakes, acceleration of land desertification and salinization, and the death of natural vegetation in large areas. If the Minqin oasis disappears, the western part of China will face the threat of desert erosion. Therefore, the ecological environment of the SRB has been highly valued by the state and is widely considered by all sectors of society. 0.57 m/year, and the maximum decrease being 15~16 m. At the beginning of the 21st century, the consumption rate of water resources in the basin reached 109%, and the development and utilization of water resources reached 172%, which is far more than the reasonable carrying capacity of the water resources in the basin. The superposition of resource-based and structural water shortages has led to the deterioration of the ecological environment of the basin, such as the decline in the groundwater level, disappearance of lakes, acceleration of land desertification and salinization, and the death of natural vegetation in large areas. If the Minqin oasis disappears, the western part of China will face the threat of desert erosion. Therefore, the ecological environment of the SRB has been highly valued by the state and is widely considered by all sectors of society.

Methods
The flowchart for this study is shown in Figure 2. First, different data sources will be collected, and the GRACE-GLDAS-derived GWS will be validated using in situ groundwater data. Then, the changes in GWS and TWS will be discussed in detail, and combined with remote-sensing and observation data, the causes of GWS changes will be explained. Finally, groundwater storage variation under the influence of human activities will be further analyzed using the water balance equation.

Methods
The flowchart for this study is shown in Figure 2. First, different data sources will be collected, and the GRACE-GLDAS-derived GWS will be validated using in situ groundwater data. Then, the changes in GWS and TWS will be discussed in detail, and combined with remote-sensing and observation data, the causes of GWS changes will be explained. Finally, groundwater storage variation under the influence of human activities will be further analyzed using the water balance equation. According to the principle of water balance, important factors, such as precipitation (P), actual evapotranspiration (ET), runoff (R) and human activity influence (Q), may con-  According to the principle of water balance, important factors, such as precipitation (P), actual evapotranspiration (ET), runoff (R) and human activity influence (Q), may contribute to the variation in water storage in the SRB with intense human activities. The terrestrial water storage change (TWSC) can be expressed as Equation (1) [25]: where ds/dt means the change in terrestrial water storage with the given time (t). Given that other items are known, the influence of human activities can be isolated as: GRACE-based TWSC is calculated as the backward difference of the TWSA [26] as:

Terrestrial Water Storage Variations
There are two main processing approaches for GRACE data: Parameterizing the Earth's gravity field using global spherical harmonics (SH) basis functions and parameterizing the gravity field with regional mass concentration functions (mascons). Spherical harmonic and unconstrained mascon data are similar as they are both based on the same fundamental GRACE satellite data and models that remove atmospheric, oceanic, and tidal signals. A basic difference between SH and mascons is that SH solutions are global, whereas mascons can be applied at regional to global scales [27]. In this study, we used the monthly 0.5 • GRACE level-3 mascon (mass concentration) datasets from the RL06 time-variable gravity field model provided by the Jet Propulsion Laboratory (JPL) (ftp://podaac-ftp.jpl.nasa.gov/allData/tellus/L3/landmass/RL06 (accessed on 18 March 2021)). The C20 coefficients were replaced with the solutions from satellite laser ranging [28], and the degree-1 coefficients were estimated using the method from Swenson et al. [29]. A glacial isostatic adjustment correction was applied for the JPL mascon based on the ICE6G-D model from Peltier et al. [30], and the data are presented as anomalies relative to the time-mean baseline of the period from January 2004 to December 2009. The gridded products used in this paper covered a total period of 173 months, spanning from January 2003 to December 2019. Data from 31 months during the study period are not available. The missing data during the period from July 2017 to May 2018 were due to the discontinuation of GRACE and GRACE-FO satellite data, and the other 20 months were interpolated by simply averaging the values of the two months before and after the month with the missing data [31]. For the consistency of data resolution, the 0.5 • of GRACE data was averaged to 1 • .

Groundwater Storage Variations
Variations in TWS included a vertically integrated measure of water storage changes in groundwater (GW), soil moisture (SM), surface water (SW) and snow water (SWE) (Equation (4)). Therefore, GWS variations can be isolated from TWS components when other components are known. ∆GWS can be calculated by removing the storage anomaly of the aboveground part from ∆TWS in accordance with Equation (5): Since reliable and spatially continuous measurements of soil moisture were not currently available throughout our study area, the outputs from land-surface models were Remote Sens. 2021, 13, 2672 6 of 19 used to provide SM and SWE [32,33]. The surface water we used was plant canopy surface water from the land-surface models.

GLDAS Model Data
GLDAS was developed by NASA's Goddard Space Flight Center and the National Centers of Environmental Prediction (NCEP). The goal of GLDAS is to ingest satellite and ground-based observational data products using advanced land-surface modelling and data assimilation techniques to generate optimal fields of land-surface states and fluxes [34]. The GLDAS-2 datasets consist of 0.25 • × 0.25 • Noah model data, covering the period from 2000 to the present. The temporal resolution for the GLDAS products was three hours. Monthly products were generated through temporal averaging or a summation of the three-hourly products [35]. The GLDAS data were retrieved from http://disc.gsfc.nasa.gov/hydrology/data-holdings (accessed on 12 March 2021). Noah was selected as the data source in this paper, which simulates changes in water storage with less bias and uncertainty than other models in GLDAS [36,37]. The dataset contained soil moisture, snow water equivalent, runoff, plant canopy surface water and ET at a spatial resolution of 0.25 • and 1 • and a temporal resolution of three hours or one month. For consistency with the GRACE data, the average SM, SWE and plant canopy anomalies versus time were computed.

Precipitation Data
This study mainly collected two precipitation products: TRMM and CMA datasets. TRMM (Tropical Rainfall Measuring Mission, TRMM 3B43, 1998-2019) provided the temporal and spatial distribution and changes in precipitation in the range of 50 • south latitude to 50 • north latitude, with a spatial resolution of 0.25 • . The monthly precipitation dataset (0.5 • ) in China from January 1980 to December 2019 was obtained from the China Meteorological Administration (CMA), which integrated the daily precipitation observed by 2472 national stations (including the National Climate Observatory and the first-level and second-level meteorological observation stations) and adopted the optimal interpolation method based on the climate background field, in order to generate the grid product of daily precipitation in China.

ET Data
The ET data used in this study included actual ET from the GLEAM, MODIS and GLDAS-Noah models. The Global Land Evaporation Amsterdam Model (GLEAM) is a set of algorithms dedicated to the estimation of terrestrial evaporation and root-zone soil moisture from satellite data. Since its development in 2011, the model has been continuously revised and updated [38]. To date, three kinds of datasets produced using the latest version of GLEAM v3.5 are currently available, and the three v3.5 datasets differ only in their forcing and spatial coverage (http://www.gleam.eu (accessed on 1 March 2021)). GLEAM v3.5a is a global dataset spanning the 39-year period, 1980~2020, at a 0.25 • spatial resolution. For more detailed information about GLEAM data, see the previous study by Miralles et al. [39].
The actual MODIS/MOD16 ET datasets were estimated according to the Penman-Monteith equation [40]. Mu et al. [41,42] modified the MODI6 actual ET algorithm several times in order to improve the accuracy of ET components, including night time ET, canopy transpiration, and soil heat flux, which has been successfully used in many regions worldwide. Previous studies have assessed the accuracy of MOD16 ET in several river basins in China [43,44]. The results showed that the fitting effect improved at the station scale, with an average correlation coefficient of 0.76. The monthly actual mean ET at a spatial resolution of 1 km was averaged for different spatial resolutions over the period 2003-2019.

Land-Cover Data from MODIS
The MODIS land-cover type product (MCD12Q1) provides a suite of science datasets (SDSs) that map global land coverage at a 500 m spatial resolution at annual time steps for six different land-cover legends, and was produced using supervised classification techniques, such as decision trees [45] and ensemble classification methods [46]. The product contains 13 science datasets, including five legacy classification schemes (IGBP, UMD, LAI, BGC, and PFT) and a new three-layer legend based on the Land Cover Classification System (LCCS) from the Food and Agriculture Organization [47], and the dataset we used to analyze land-use change in the study area is the IGBP classification scheme. More details can be found at https://doi.org/10.5067/MODIS/MCD12Q1.006 (accessed on 5 March 2021).

Mutation Test Methods
We used two mutation test methods to analyze the mutation characteristics of the GWSA trend in the SRB.

Mann-Kendall Nonparametric Mutation Test
The Mann-Kendall [48,49] mutation test [50] is a nonparametric statistical test method. Its advantage is that it is not only easy to calculate, but it can also clarify the starting point of the mutation and can identify the mutation period; it is commonly used for trend analysis and variation testing of hydrological and meteorological elements, such as rainfall, temperature and runoff. For time series x with n sample sizes, the rank column S k of the sequence needs to be constructed in the following manner: where S k is the cumulative number of values at moment i, which is greater than that at moment j. Under the assumption of a random and independent time series, the following statistics can be determined: where UF 1 = 0, E(S k ) and Var(S k ) are the mean and variance of S k , respectively. When x 1 , x 2 , and x n are independently and continuously distributed, they can be calculated by the following equations: where UF i is a standard normal distribution. For a given significance level, the critical value U a can be obtained by querying the normal distribution table. If |UF i | > U a (U 0.05 = 1.96, U 0.01 = 2.58), it indicates that the sequence has an obvious trend of increasing or decreasing. The time series x is generated into its corresponding inverse sequence The calculation process above is repeated to obtain the inverted sequence statistical value UB k (UB k = −UF k , k = n, n − 1, . . . , 1, UB 1 = 0). When the UF k and UB k curves intersect between the critical boundary, the intersection point can be used as the beginning of the mutation [51].
Moving T-Test Method The moving T-test assumes that the front of the sliding point is sequence 1 and the rear of the sliding point is sequence 2. The samples with capacities of n 1 and n 2 are extracted from the two sequences, and the assumption that the distribution functions of the two sequences are equal is tested: where x 1 and x 2 are the mean values of the two samples and T obeys a T(n 1 + n 2 − 2) distribution. For a given significance level a, the original assumption is rejected when |T| > t a/2 , the two samples have a significant difference and the sequence is mutated.

Correlation Analysis and Data
This study used correlation analysis to study the relationship between different variables. The results are expressed by the correlation coefficient r, which is used to describe the strength of the linear correlation between two variables. The formula is as follows: where X i and Y i are sample values and X and Y are the mean values of the sample values, respectively. The data used in this study are listed in Table 1.  Figure 3 shows monthly and yearly average changes in precipitation from the TRMM and CMA datasets in the SRB. The two datasets have a high consistency, with a correlation coefficient of 0.97 at the monthly scale and an RMSE (root mean square error) of 6.18 mm/month, and the TRMM data show a 10 mm underestimation in the rainy season at the monthly scale. The yearly precipitation in the SRB showed a fluctuation at a small magnitude of 50 mm from 2003 to 2015. However, precipitation showed an obvious upward trend after 2015, with an average annual precipitation of 276 mm. Precipitation in the rainy season accounts for 70% of the annual value. Due to the high consistency of the two sets of data and the higher resolution of the TRMM dataset than that of the CMA dataset, TRMM Remote Sens. 2021, 13, 2672 9 of 19 data was chosen for the subsequent analysis. To facilitate the comparative analysis with GRACE data, the resolution of precipitation data was upscaled to 1 • . coefficient of 0.97 at the monthly scale and an RMSE (root mean square error) of 6.18 mm/month, and the TRMM data show a 10 mm underestimation in the rainy season at the monthly scale. The yearly precipitation in the SRB showed a fluctuation at a small magnitude of 50 mm from 2003 to 2015. However, precipitation showed an obvious upward trend after 2015, with an average annual precipitation of 276 mm. Precipitation in the rainy season accounts for 70% of the annual value. Due to the high consistency of the two sets of data and the higher resolution of the TRMM dataset than that of the CMA dataset, TRMM data was chosen for the subsequent analysis. To facilitate the comparative analysis with GRACE data, the resolution of precipitation data was upscaled to 1°.         [52] in the Hexi Corridor, and the average declining rates were 1.58 mm/year, and 2.95 mm/year, respectively. The correlation coefficient between TWSA and GWSA was 0.81, indicating that the variation in TWS was well correlated with the variation in GWS in the SRB. Moreover, annual precipitation and actual evapotranspiration in the SRB increased to varying degrees, and the rate of annual actual evapotranspiration was slightly higher than that of the annual precipitation, with rates of 2.55 mm/year, and 1.18 mm/year, respectively. By analyzing the relationship between meteorological factors and water storage, it can be found that, in the years with rapid declines in water storage, such as 2009, 2014, 2016 and 2017, the actual ET was greater than that of the precipitation in the previous year. Due to drought, upstream inflow was reduced, and exploitation increased, which eventually led to a decline in GWS [53]. Using correlation analysis among TWSA, GWSA and the difference between the previous year's precipitation and ET both show a correlation, and the correlation coefficients are greater than 0.6, indicating that the variation in water storage in the SRB was affected by meteorological factors to some extent.   [52] in the Hexi Corridor, and the average declining rates were 1.58 mm/year, and 2.95 mm/year, respectively. The correlation coefficient between TWSA and GWSA was 0.81, indicating that the variation in TWS was well correlated with the variation in GWS in the SRB. Moreover, annual precipitation and actual evapotranspiration in the SRB increased to varying degrees, and the rate of annual actual evapotranspiration was slightly higher than that of the annual precipitation, with rates of 2.55 mm/year, and 1.18 mm/year, respectively. By analyzing the relationship between meteorological factors and water storage, it can be found that, in the years with rapid declines in water storage, such as 2009, 2014, 2016 and 2017, the actual ET was greater than that of the precipitation in the previous year. Due to drought, upstream inflow was reduced, and exploitation increased, which eventually led to a decline in GWS [53]. Using correlation analysis among TWSA, GWSA and the difference between the previous year's precipitation and ET both show a correlation, and the correlation coefficients are greater than 0.6, indicating that the variation in water storage in the SRB was affected by meteorological factors to some extent. rapid declines in water storage, such as 2009, 2014, 2016 and 2017, the actual ET was greater than that of the precipitation in the previous year. Due to drought, upstream inflow was reduced, and exploitation increased, which eventually led to a decline in GWS [53]. Using correlation analysis among TWSA, GWSA and the difference between the previous year's precipitation and ET both show a correlation, and the correlation coefficients are greater than 0.6, indicating that the variation in water storage in the SRB was affected by meteorological factors to some extent.  Surface runoff mainly occurs in the Qilian Mountain area in the upper reaches of the SRB. There are two agricultural irrigated areas (G4 and G5 grid blocks), located in the middle and downstream areas in the SRB, which do not have runoff production. Since the Hongyashan Reservoir was built in 1958, surface water from the Shiyang River into Minqin has been completely controlled by the reservoir, and all the water has been introduced into the area for irrigation.

Variations of TWSA Data and Their Components
GRACE-derived GWS changes represent regional groundwater storage change under the ground. However, observations of well data showed a groundwater level change at the point scale. Many researchers have compared GRACE-derived GWS data with in situ measured data [54][55][56], and the correlation coefficients ranged from 0.60 to 0.85. In this study, the daily groundwater level (GWL) in 71 observation wells in the SRB, from January 2007 to December 2018, was recorded, and 18 wells had complete monitoring data, which were mainly concentrated in the middle and lower reaches of the SRB (G4, G5). Therefore, the observation of well data in these two grids were averaged and compared to the GRACE data. Figure 7a shows that the change in the groundwater level in the G4 and G5 grid blocks was basically consistent with the trend in the GRACE data, and the correlation coefficient was 0.68 (p < 0.01). Therefore, we think GRACE-derived GWSA is reasonable and consistent with previous studies, which can reflect the dynamic changes of groundwater in SRB. In addition, by observing the change in the groundwater level in this area, it can be found that the average annual decline in the groundwater level in the middle and lower reaches of the SRB before 2012 was approximately 0.6 m, the change in groundwater level in 2012-2015 was almost stable, and the groundwater level increased significantly after 2015, in which the change was approximately 0.4 m/year. In addition, the groundwater level in the middle and lower reaches of the SRB also showed strong annual fluctuations before 2012. From July to August, the groundwater level reached its lowest level due to irrigation and pumping in spring. From autumn to spring, the groundwater level recovered in the second year due to the decrease in irrigation and pumping, in addition to possible lateral replenishment from rivers [57]. Annual changes in groundwater levels after 2012 were more stable than before. This may have been related to the closure of a large number of drilling wells, the completion of field water savings and cross-basin water transfer and other measures in the key governance project of the SRB, indicating the effect of strengthening groundwater management and implementing strict groundwater control and mining systems. According to the water resource bulletin data of the SRB in 2010, the amount of groundwater extraction in the SRB was 0.61 billion m 3 , which was 47.6% lower than the value of 1.

Abrupt Change in GWS Variations
From the trend of GWS, as of 2016, GWS in the middle and lower reaches of the SRB have been in a downward trend, and the change was −44 mm EWH, as compared to 2003, which is equivalent to a decrease of 8.8 × 10 8 m³ in GWS and an annual decrease of 0.6 × 10 8 m³. In addition, the rate of decline in GWS changed significantly before and after 2009, accounting for 1.2 × 10 8 m³/year, and 0.05 × 10 8 m³/year, respectively. The M-K mutation test was used to detect the GWS mutation point (Figure 8a). From 2004 to 2016, GWS has been significantly decreasing ( < −1.96), and there was a possible mutation point from 2007 to 2008. Then, the moving T-test method was used for mutation detection ( Figure  8b). The step size was set to two, the significance level was 0.05, and the critical value was 4.303. At this significance level, the maximum absolute value of T appeared in 2009 (p < 0.05); therefore, according to the moving T-test method, a possible mutation point in 2009 was detected. By combining the trends of GWS and the two mutation testing methods, we found that, after a series of management plans in the SRB, the decline in GWS was curbed close to 2009.   Figure 1. The groundwater level of each well fluctuated greatly before 2012 and then gradually stabilized, and each well only differed in the amplitude of the water level. Groundwater level fluctuations in monitoring wells near the Shiyang River (W2) in G4 were small, and the amplitude was approximately 1 m. Groundwater level fluctuations in monitoring wells far from the river (W1) were large, which may have been related to the hydraulic connection between the surface and groundwater. The wells near the river easily received river recharge and maintained a stable water level. The groundwater level in the downstream area (G5) had obvious fluctuations before 2012 and then became stable.

Abrupt Change in GWS Variations
From the trend of GWS, as of 2016, GWS in the middle and lower reaches of the SRB have been in a downward trend, and the change was −44 mm EWH, as compared to 2003, which is equivalent to a decrease of 8.8 × 10 8 m 3 in GWS and an annual decrease of 0.6 × 10 8 m 3 . In addition, the rate of decline in GWS changed significantly before and after 2009, accounting for 1.2 × 10 8 m 3 /year, and 0.05 × 10 8 m 3 /year, respectively. The M-K mutation test was used to detect the GWS mutation point (Figure 8a). From 2004 to 2016, GWS has been significantly decreasing (UF k < −1.96), and there was a possible mutation point from 2007 to 2008. Then, the moving T-test method was used for mutation detection (Figure 8b). The step size was set to two, the significance level was 0.05, and the critical value was 4.303. At this significance level, the maximum absolute value of T appeared in 2009 (p < 0.05); therefore, according to the moving T-test method, a possible mutation point in 2009 was detected. By combining the trends of GWS and the two mutation testing methods, we found that, after a series of management plans in the SRB, the decline in GWS was curbed close to 2009. 8b). The step size was set to two, the significance level was 0.05, and the critical value was 4.303. At this significance level, the maximum absolute value of T appeared in 2009 (p < 0.05); therefore, according to the moving T-test method, a possible mutation point in 2009 was detected. By combining the trends of GWS and the two mutation testing methods, we found that, after a series of management plans in the SRB, the decline in GWS was curbed close to 2009.

Evaluation of the Water Balance before and after the Key Management Plan
From the change in GWS in each governance stage (Table 2), it can be seen that GWS still decreased by 24.3 mm EWH after the end of the first stage of governance (2007)(2008)(2009)(2010). With further strengthening of governance, GWS decreased by 2.2 mm EWH in 2011-2016, which was 91% lower than that in 2007-2010; and by 2019, GWS increased by 9.1 mm EWH, as compared to 2016. After the completion of the water diversion project from the Xiying River to downstream of Minqin, known as the Caiqi Special Water Diversion Canal, 1.33 × 10 8 m 3 and 1.35 × 10 8 m 3 of water was transferred to Minqin and Caiqi in 2010 and 2011, respectively, which greatly increased the inflow to the Minqin Basin. In addition, on the basis of data from 2010, water-saving reconstruction projects were further carried out in the Huangyang, Gulang, Donghe and Qinghe irrigation areas, which effectively improved the water utilization coefficient of irrigation channels and the utilization efficiency of water resources [58]. Table 2. The change in GWS in the middle and lower reaches of the SRB in each governance stage. The water volume is calculated by multiplying the area of the middle and lower reaches of the SRB (20,000 km 2 ). The monthly average GWSA in the middle and lower reaches of the SRB from 2003 to 2016 were calculated. As shown in Figure 9a, the loss of GWS in the middle and lower reaches of the SRB was the least in July and the most in October. Variations in GWS from January to June were relatively close, and from July to October, the decline in GWS intensified; after October, the sharp decline in GWS was gradually alleviated. In addition, the change in SM in the middle and lower reaches was significantly negatively correlated with the change in GWS (r = −0.9). Under the influence of the rainy season and agricultural irrigation, soil water reserves reached a high value from September to October, while groundwater reserves decreased due to lagging exploitation and irrigation recharge. The influence of human activities calculated by the water balance equation was compared to groundwater extraction from the water resource bulletin of the SRB (Figure 10). Since the human impact calculated by the water balance is expressed as the equivalent water height, the water volume can be converted by multiplying the area of the middle and lower reaches of the SRB (20,000 km 2 ). As shown in Figure 10, the influence of human activities had an upward trend before 2007 and then began to decline, which is consistent with the change in the extraction recorded in the water resource bulletin. After 2010, the extraction of groundwater in the SRB was basically stable at 8 × 10 8 m³/year, while the human influence obtained by the water balance equation still showed a downward trend after 2010, and the human influence estimated by the water balance was greater than the groundwater depletion, which may have been related to the implementation of watersaving renovation projects and industrial structural deployment in the irrigation area. Due to the large number of drilling wells in the study area, there may have been some wells that were not installed in metering facilities, which may have also led to the results deviation in the bulletin. In addition, the spatial resolution of the impacts from human activities obtained by the water balance method is low, and there may be large errors when applied to small-scale areas such as the SRB. In addition, the error of human influence using different remote-sensing data is about 20%. However, the overall trend still shows that the reduction in water storage, caused by human activities, was effectively reduced by implementing a series of governance activities in the SRB.  With little change in annual precipitation and evapotranspiration, the change in GWS from January to August was more stable than before, but from August to October, GWS still had an obvious downward trend. After 2011, the standard deviation of the monthly changes in GWS became 4.3 mm EWH, which was a little lower than the 7.4 mm EWH for the period from 2007 to 2010.

Change of GWS During the Period Change Rate Equivalent Water Height
The influence of human activities calculated by the water balance equation was compared to groundwater extraction from the water resource bulletin of the SRB (Figure 10). Since the human impact calculated by the water balance is expressed as the equivalent water height, the water volume can be converted by multiplying the area of the middle and lower reaches of the SRB (20,000 km 2 ). As shown in Figure 10, the influence of human activities had an upward trend before 2007 and then began to decline, which is consistent with the change in the extraction recorded in the water resource bulletin. After 2010, the extraction of groundwater in the SRB was basically stable at 8 × 10 8 m 3 /year, while the human influence obtained by the water balance equation still showed a downward trend after 2010, and the human influence estimated by the water balance was greater than the groundwater depletion, which may have been related to the implementation of water-saving renovation projects and industrial structural deployment in the irrigation area. Due to the large number of drilling wells in the study area, there may have been some wells that were not installed in metering facilities, which may have also led to the results deviation in the bulletin. In addition, the spatial resolution of the impacts from human activities obtained by the water balance method is low, and there may be large errors when applied to small-scale areas such as the SRB. In addition, the error of human influence using different remote-sensing data is about 20%. However, the overall trend still shows that the reduction in water storage, caused by human activities, was effectively reduced by implementing a series of governance activities in the SRB. deviation in the bulletin. In addition, the spatial resolution of the impacts from human activities obtained by the water balance method is low, and there may be large errors when applied to small-scale areas such as the SRB. In addition, the error of human influence using different remote-sensing data is about 20%. However, the overall trend still shows that the reduction in water storage, caused by human activities, was effectively reduced by implementing a series of governance activities in the SRB.   After 2016, the cropland area decreased significantly, and the change rate was −106.6 km 2 /year. By analyzing the variation of GWS in the middle and lower reaches, it can be found that the variation of GWS had a strong negative correlation with the change in cropland area, and the correlation coefficient was −0.94 (p < 0.01), indicating that the variation in GWS in the middle and lower reaches was closely related to cropland cultivation; reducing cropland areas had a significant effect on restoring GWS. With the decrease in cropland area and the rise in GWS, the barren area in the middle and lower reaches was gradually reduced, and natural vegetation, such as grassland, was gradually restored, indicating that the ecosystems in the SRB were in a recovery stage [59]. After 2016, the cropland area decreased significantly, and the change rate was −106.6 km 2 /year. By analyzing the variation of GWS in the middle and lower reaches, it can be found that the variation of GWS had a strong negative correlation with the change in cropland area, and the correlation coefficient was −0.94 (p < 0.01), indicating that the variation in GWS in the middle and lower reaches was closely related to cropland cultivation; reducing cropland areas had a significant effect on restoring GWS. With the decrease in cropland area and the rise in GWS, the barren area in the middle and lower reaches was gradually reduced, and natural vegetation, such as grassland, was gradually restored, indicating that the ecosystems in the SRB were in a recovery stage [59].

Conclusions
To understand the influences on water resources of the implemented key governance project of the Shiyang River Basin in 2007, this article studied the changes in terrestrial water and groundwater storage from 2003 to 2019. GRACE-GLDAS data were used to quantitatively evaluate the changes in groundwater storage before and after the SRB key management project. Combined with precipitation, actual ET, land use, bulletin and in situ observation data, the natural and human factors affecting the changes in water storage were analyzed, and the situation of ecological restoration in the SRB was evaluated from the perspective of water resources. The conclusions are as follows: 1. The variation in GWS obtained by GRACE and GLDAS was compared to the data from groundwater observation wells, and the results from the two sources show consistency of the wells (r = 0.68), indicating that the GRACE was reliable for GWS monitoring in the SRB. GWS in the SRB had a slow downward trend from 2003 to 2016. After 2018, GWS increased by approximately 0.4 mm/mon, which was approximately

Conclusions
To understand the influences on water resources of the implemented key governance project of the Shiyang River Basin in 2007, this article studied the changes in terrestrial water and groundwater storage from 2003 to 2019. GRACE-GLDAS data were used to quantitatively evaluate the changes in groundwater storage before and after the SRB key management project. Combined with precipitation, actual ET, land use, bulletin and in situ observation data, the natural and human factors affecting the changes in water storage were analyzed, and the situation of ecological restoration in the SRB was evaluated from the perspective of water resources. The conclusions are as follows: 1.
The variation in GWS obtained by GRACE and GLDAS was compared to the data from groundwater observation wells, and the results from the two sources show consistency of the wells (r = 0.68), indicating that the GRACE was reliable for GWS monitoring in the SRB. GWS in the SRB had a slow downward trend from 2003 to 2016. After 2018, GWS increased by approximately 0.4 mm/mon, which was approximately 3.84 × 10 8 m 3 per year. Additionally, the GRACE-observed TWSA was fairly well-correlated with the GWSA (r = 0.81).

2.
Different datasets of precipitation and actual ET data in the SRB were used to analyze their relationship with the GWSA. There was a correlation between the variation in GWS in the current year and the difference between precipitation and actual ET change in the previous year (r = 0.61 By using the water balance method, the change in GWS affected by human activities was obtained, and the influencing factors were discussed. The change in GWS affected by human activities rose before 2007 and then decreased, which was consistent with the change in groundwater depletion in the SRB. Land-use change in the middle and lower reaches of the SRB also showed that the increasing rate of the crop area gradually slowed down after 2011 and significantly decreased after 2016. Natural vegetation areas gradually increased, and GWS also gradually rose, which indicated that, through a series of governance activities in the SRB, the reduction in water storage caused by human activities was effectively reduced, and the ecological environment of the SRB gradually recovered.
GRACE/GRACE-FO data could provide accurate estimation of regional groundwater storage changes to support groundwater resource management. In desert and remote areas, GRACE is the only means to evaluate groundwater reserve change. However, due to the coarse resolution of GRACE data, the applicability of GRACE-derived GWSA data in smallscale areas should be explored. In addition, the uncertainties in multisource remote-sensing data used in the principle of water balance may also affect the accuracy of human impacts. However, in the absence of groundwater extraction data, multi-source remote-sensing data and water balance methods can provide feasible solutions for assessing the impact of human activities on groundwater storages.
Author Contributions: Conceptualization, X.L. and L.H.; data curation, X.L., K.S. and Z.Y.; formal analysis, X.L. and K.S.; methodology, X.L. and L.H.; visualization, Z.Y. and J.S.; writing-original draft, X.L.; writing-review and editing, X.L., L.H. and W.Y. All authors will be informed about each step of manuscript processing, including submission, revision, revision reminder, etc. via emails from our system or assigned Assistant Editor. All authors have read and agreed to the published version of the manuscript.