Long-term groundwater storage changes and land subsidence development in the North China Plain (1971–2015)

Hydrogeology Journal, Apr 2018

The North China Plain (NCP) has been suffering from groundwater storage (GWS) depletion and land subsidence for a long period. This paper collects data on GWS changes and land subsidence from in situ groundwater-level measurements, literature, and satellite observations to provide an overview of the evolution of the aquifer system during 1971–2015 with a focus on the sub-regional variations. It is found that the GWS showed a prolonged declining rate of −17.8 ± 0.1 mm/yr during 1971–2015, with a negative correlation to groundwater abstraction before year ~2000 and a positive correlation after ~2000. Statistical correlations between subsidence rate and the GWS anomaly (GWSA), groundwater abstraction, and annual precipitation show that the land subsidence in three sub-regions (Beijing, Tianjin, and Hebei) represents different temporal variations due to varying driver factors. Continuous drought caused intensive GWS depletion (−76.1 ± 6.5 mm/yr) and land subsidence in Beijing during 1999–2012. Negative correlations between total groundwater abstraction and land subsidence exhibited after the 1980s indicate that it may be questionable to infer subsidence from regional abstraction data. Instead, the GWSA generally provides a reliable correlation with subsidence. This study highlights the spatio-temporal variabilities of GWS depletion and land subsidence in the NCP under natural and anthropogenic impacts, and the importance of GWS changes for understanding land subsidence development.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:

https://link.springer.com/content/pdf/10.1007%2Fs10040-018-1768-4.pdf

Long-term groundwater storage changes and land subsidence development in the North China Plain (1971–2015)

Hydrogeology Journal https://doi.org/10.1007/s10040 Long-term groundwater storage changes and land subsidence development in the North China Plain (1971-2015) Huili Gong 0 1 2 3 4 5 6 Yun Pan 0 1 2 3 4 5 6 Longqun Zheng 0 1 2 3 4 5 6 Xiaojuan Li 0 1 2 3 4 5 6 Lin Zhu 0 1 2 3 4 5 6 Chong Zhang 0 1 2 3 4 5 6 Zhiyong Huang 0 1 2 3 4 5 6 Zhiping Li 0 1 2 3 4 5 6 Haigang Wang 0 1 2 3 4 5 6 Chaofan Zhou 0 1 2 3 4 5 6 Yun Pan 0 1 2 3 4 5 6 0 State Key Laboratory of Earth Surface Processes and Resource Ecology, Faculty of Geographical Science, Beijing Normal University , Beijing 100875 , China 1 College of Resource Environment and Tourism, Capital Normal University , Beijing 100048 , China 2 Beijing Laboratory of Water Resources Security, Capital Normal University , Beijing 100048 , China 3 China Geological and Environmental Monitoring Institute , Beijing 100081 , China 4 Beijing Institute of Hydrogeology and Engineering Geology , Beijing 100048 , China 5 Department of Earth Sciences, The University of Hong Kong , Hong Kong , China 6 Institute of Land Surface System and Sustainable Development, Faculty of Geographical Science, Beijing Normal University , Beijing 100875 , China The North China Plain (NCP) has been suffering from groundwater storage (GWS) depletion and land subsidence for a long period. This paper collects data on GWS changes and land subsidence from in situ groundwater-level measurements, literature, and satellite observations to provide an overview of the evolution of the aquifer system during 1971-2015 with a focus on the sub-regional variations. It is found that the GWS showed a prolonged declining rate of −17.8 ± 0.1 mm/yr during 1971-2015, with a negative correlation to groundwater abstraction before year ~2000 and a positive correlation after ~2000. Statistical correlations between subsidence rate and the GWS anomaly (GWSA), groundwater abstraction, and annual precipitation show that the land subsidence in three sub-regions (Beijing, Tianjin, and Hebei) represents different temporal variations due to varying driver factors. Continuous drought caused intensive GWS depletion (−76.1 ± 6.5 mm/yr) and land subsidence in Beijing during 1999-2012. Negative correlations between total groundwater abstraction and land subsidence exhibited after the 1980s indicate that it may be questionable to infer subsidence from regional abstraction data. Instead, the GWSA generally provides a reliable correlation with subsidence. This study highlights the spatio-temporal variabilities of GWS depletion and land subsidence in the NCP under natural and anthropogenic impacts, and the importance of GWS changes for understanding land subsidence development. Over-abstraction; Subsidence; Drought; GRACE; China Introduction Groundwater sustains social and environmental development in many arid and semi-arid regions; however, over-abstraction of groundwater, induced by drought and large water demand for agricultural and domestic use, has become a global issue (Famiglietti 2014; Richey et al. 2015) . The Gravity Recovery and Climate Experiment (GRACE) satellites show that many aquifers around the world are being depleted at amazing rates such as in northwest India (Long et al. 2016; Rodell et al. 2009) , Central Valley in California, USA (Famiglietti et al. 2011), and the North China Plain (NCP; Feng et al. 2013; Huang et al. 2015; Pan et al. 2017) . More information can be found from the review work by Chen et al. (2016a and Wada (2016) . If only focusing on the groundwater depletion rate, the NCP (Fig. 1) does not seem to be the worst among the wellknown aquifers worldwide (Chen et al. 2016a; Famiglietti et al. 2011; Rodell et al. 2009) ; however, groundwater overabstraction accompanied by the compression of aquifer systems has led to serious land subsidence in NCP (Guo et al. 2015; Ye et al. 2016) . The subsidence-affected area in the NCP is around 120,000 km2 (86% of the total area), and an area of 17,500 km2 has been subsiding at a rate of ≥20 mm/yr (Ye et al. 2016). The Interferometric Synthetic Aperture Radar (InSAR) measured a subsidence rate up to 15 cm/yr in Beijing during 2003–2011 (Chen et al.2016b; Zhang et al. 2016) . The accumulated displacement of the city of Cangzhou, Hebei province, one of the most serious land subsidence sites in NCP, is up to around 3 m. Land subsidence may cause increased flood risk, structural damage, and seawater intrusion, and result in significant socio-economic distress (Hwang et al. 2016). The highly developed megacities (e.g., Beijing) with lots of large constructions, high-speed railways, and dense population, are under the potential risk of land surface subsiding. One may expect amelioration regarding the situation of groundwater over-abstraction and land subsidence after the implementation of the South-to-North Water Diversion Project in December 2014, which was designed to transfer 9.6 km3 of water per year from the south (Danjiangkou reservoir) to the north (including Beijing, Tianjin, Hebei and Henan). The government plans to reduce groundwater use in Beijing, Tianjin and Hebei province to 65% of current levels by 2023. Such a dramatic change to groundwater consumption may recover the groundwater systems in NCP to some extent. Numerous studies have explored the groundwater storage (GWS) depletion and land subsidence in NCP; however, previous studies usually focused on a limited time span in this area. For instance, Huang et al. (2015) and Pan et al. (2017) quantified and analyzed the GWS depletion during 2003–2014, while a few other researchers studied the land subsidence during 1992–2014 (e.g. Chen et al. 2016b, 2017; Zhang et al. 2016; Zhu et al. 2015) . In order to provide more information for future groundwater resources management, more studies are still needed to investigate the GWS dynamics and the evolution of the aquifer system from a long-term perspective (Castellazzi et al. 2016a, b) . This study aims to explore the long-term changes (1971– 2015) of GWS and land subsidence in the NCP, as well as their relationship using in situ observations and satellite measurements from GRACE and InSAR data collected from literature. The spatial and temporal variations of land subsidence will be discussed to provide a historical overview of the groundwater system in NCP and its three sub-regions, i.e., Beijing, Tianjin and Hebei provinces. Data and methods In situ groundwater level observations Historical monthly in situ groundwater level data (Fig. 1) from 1971 to 2013 were collected, including the records of 44 wells in NCP during 1971–1983, obtained from Water Resources Yearbook of Haihe River Basin compiled by the Ministry of Water Resources of China; 5 wells in Beijing during 1971– 2013 provided by Beijing Institute of Hydrogeology and Engineering Geology; and 81 wells in NCP during 2005– 2013 from the Groundwater Level Almanac compiled by the China Institute of Geological Environment Monitoring. There was a data gap for the entire NCP during 1984–2002. In this study, the continuous groundwater level observations of the five wells in Beijing, accompanied with reported data about the groundwater level changing rate in NCP, were interpolated to obtain a regional averaged estimate of monthly groundwater storage anomaly (GWSA) during 1984–2002. The observed groundwater level changes were converted into GWS changes by multiplying the specific yield according to Zhang et al. (2009; Fig. 1b) as used by Huang et al. (2015) . A 20% uncertainty is assumed for the specific yield. The monthly GWSA was then obtained by subtracting mean GWS changes of the whole period (1971–2015). GRACE-derived groundwater storage anomalies Monthly GRACE level-3 Release-05 mascon (i.e., mass concentration) products (2003–2015) from the Center for Space Research (CSR) at the University of Texas and Jet Propulsion Laboratory (JPL), USA are utilized. GRACE mascon solutions represent total terrestrial water storage (TWS) anomalies. In arid region like NCP, the TWS values are the aggregation of surface-water storage in lakes and reservoirs, soil moisture and groundwater. In order to estimate GRACE-derived GWS anomalies, the observed reservoir water storage (from 22 primary reservoirs) and the averaged soil moisture storage (SMS) simulated by four land surface models (Noah, Mosaic, Community Land Model and Variable Infiltration C a p a c i t y M o d e l ) f r o m t h e G l o b a l L a n d D a t a Assimilation System (GLDAS; Rodell et al. 2004) were subtracted from the mascon-based TWS anomalies (TWSA). The modeled SMS data from GLDAS are validated against in situ observations from 97 monitoring si tes (Fig. 1 ) p ro vi de d b y C h in a Me te or ol og ic al Administration. The GRACE mascon solution has reduced leakage from land to ocean, and no additional leakage correction is required (Scanlon et al. 2016) . More details about the mascon solution can be found in Save et al. (2016) and Watkins et al. (2015) . The missing values in GRACE data are conservatively estimated through a cubic-spline interpolation. The uncertainties of GWSA are estimated using the propagated errors from GRACE TWSA provided by mascon solutions and the modeled SMS estimated using the standard deviation of the four model simulations. Land subsidence The regional averaged annual land subsidence rates in Beijing eastern plain (1971–2014), Cangzhou (1971–2004, 2007, 2009–2010, 2013–2014), and Tianjin urban area (1971– 2013) were obtained from previous studies (Wang 2004; Liu et al. 2016; Chen et al. 2010; Pan et al. 2004; Guo 2016) . These data were obtained through ground leveling conducted by the institutes of local geologic survey; additionally, the land subsidence in Beijing plain area was estimated from 29 RADARSAT-2 InSAR images from 2011 to 2015. They were processed using the small baseline approaches following Zhou et al. (2016) . Auxiliary data Groundwater abstraction and precipitation data were used for correlation analysis with the GRACE-derived GWSA. Annual groundwater abstraction data (1971–2015) were provided by Beijing Institute of Hydrogeology and Engineering Geology and collected from Water Resources Yearbook of Haihe River Basin. Monthly gridded (0.5° × 0.5°) precipitation data (version 2.0) are obtained from the China Metrological Administration (CMA). This dataset has been produced through a thin plate spline interpolation using the data from 2,472 rain gauge stations across China. It has been well validated over the study area with 1,684 rain gauge records (2010–2012) in a previous study by Pan et al. (2017) . Results and discussion Long-term GWS depletion exacerbated by droughts Figure 2 plots the monthly GWSA of NCP during the period 1971–2015. It is apparent that the GWS in NCP showed prolonged declining trend since the 1970s with a rate of −17.8 ± 0.1 mm/yr (Fig. 2a; Table 1). However, variabilities of GWS depletion associated with the changes of precipitation and abstraction can be found for the sub-periods (Fig. 2b,c). The GWSA slightly decreased (−6.2 ± 1.3 mm/yr) during 1971–1980, followed by a rapid declining (−20.6 ± 0.2 mm/yr) during 1981–2002 with less precipitation and more abstraction. In recent years (2003–2015), the GWSA represented a similar declining trend (−18.2 ± 0.2 mm/yr) with slightly increased mean precipitation and abstraction as compared with that during 1981–2002, despite the fact that the abstraction gradually reduced year by year. The long-term declining trend is a combined result from natural and anthropogenic factors, as drought usually lead to more intensive groundwater abstraction due to increased agricultural water demand (Pan et al. 2017) . For continuous dry years (e.g. 1997– 2002) , the groundwater depletion rate can be as large as −34.7 ± 0.9 mm/yr. During 1971–2015, Beijing experienced the largest GWS depletion rate (−34.1 ± 2.4 mm/yr) as compared with Tianjin (−7.0 ± 0.4 mm/yr) and Hebei (−15.2 ± 0.6 mm/yr; Fig. 3). An aggravated depletion (−76.1 ± 6.5 mm/yr) occurred in Beijing since 1999 due to the continuous drought (Chen et al. 2017) . In Hebei, the GWS depletion in recent years (2005–2013) is only half of the historical record (−14.7 ± 1.1 mm/yr) of 1971– 1990. The GWS in Tianjin experienced accelerated depletion in recent years (2005–2013) at the rate of −20.2 ± 2.0 mm/yr. Figure 4 presents temporal relationships between the GWSA, precipitation, and groundwater abstraction for three sub-regions, i.e. Beijing, Tianjin and Hebei. Poor relationships exist between GWSA and precipitation (Fig. 4a–c), while significant correlations can be found between the GWSA and abstraction (Fig. 4d–f). In the early depletion stage (GWSA >0), negative correlations indicate that increased abstraction can significantly deplete the aquifer. In contrast, positive relationships exist in the late depletion stage (GWSA <0), indicating that the declining trend cannot be quickly recovered through reducing abstraction. This is because the reduced abstraction is still unsustainable and it cannot be totally offset by groundwater recharge (Zhang et al. 2009) . The response of the GWSA to the changes of abstraction varies among the sub-regions. The GWSA of Beijing is more sensitive to precipitation and abstraction as compared to Tianjin and Hebei (Fig. 4g–i). It means that the GWS of Beijing can be significantly recovered by intensive precipitation (e.g. 1985–1997), or be depleted by continuous droughts (e.g. 1999–2012) . Table 1 illustrates the least-squares-fitted seasonal terms (i.e., annual and semi-annual amplitudes and phases) and long-term trends of GWSA. In situ GWlevel observations from 2005 to 2013 are used in comparison with GRACE-derived GWSA. The two GWSA estimates show a close depletion rate of −17.7 ± 1.1 mm/yr (GRACE) and − 19.2 ± 1.7 mm/yr (in situ) during 2005– 2013, with a correlation coefficient of 0.84 and root mean square error (RMSE) of 37.6 mm. Regarding the larger amplitude represented by water wells comparing to that from GRACE (Fig. 2; Table 1), the reasons are probably due to errors associated with specified yield and the GRACE data (Scanlon et al. 2012) . Joint use of the results from GRACE and in situ observations lead to caution about the uncertainties. Sub-regional land subsidence responding to natural and anthropogenic factors Continuous GWS depletion has led to a large area of land subsidence in the NCP, with the accumulated deformation up to ~3 m in Cangzhou (Fig. 5). The changes of the GWS and land subsidence of sub-periods are given in Table 2. Large subsidence rates (>50 mm/yr) mainly occurred in the late (2003–2015), middle (1980– 2000), and early (1971–1987) periods in Beijing, Hebei, and Tianjin, respectively. The land subsidence generally responds to the changes of groundwater abstraction and GWSA with a good correlation, while it shows no significant correlations with annual precipitation (Fig. 6). During the early stage (before 1980s), positive correlations demonstrate that the subsidence was mainly regulated by abstraction and, thus, increased pumping would lead to worse subsidence. Negative correlations existing in the late stage (after the 1980s) indicate that the land subsidence continued with reduced regional abstraction during this period. This has probably resulted from the changed hydrogeological conditions due to the drawdown of groundwater level. Budhu and Adiyaman (2013) found that subsidence continued when pumping is stopped if the groundwater level is below the elevation of the clay zones. In contrast, the GWSA represents consistent negative correlations with subsidence during the whole period, except for Hebei in the late stage because subsidence data of the city of Cangzhou used for analysis fails to represent regional land subsidence of Hebei province during this period. The different correlations represented by abstraction and GWSA, as well as their variations in space and time, Fig. 2 a Monthly time series of the GWSA in NCP from 1971 to 2015. The blue curves are the GWSA estimated from in situ groundwater-level observations, while the red curves are that from GRACE data. The subplot shows the blown-up GWSA time series from GRACE and in situ observations during 2005–2013. The shaded areas are the monthly GWSA errors. Annual rainfall (b) and annual groundwater abstraction data (c) are shown for reference. The straight lines are the linearly fitted trends (a) and long-term mean annual rainfall (b) and groundwater abstraction (c) during the corresponding sub-period resulted from the combined effects of natural and anthropogenic factors. Taking Beijing as an example, GWS depletion due to increasing abstraction caused a notable development of subsidence during 1970s, followed by stable GWS c ha nges and subsidence during 1980–1998; however, the 9-year continuous drought that started from 1999 again caused a rapid increase of subsidence, despite reduced abstraction. In this case, the abstraction failed to represent subsidence development with a positive correlation shown in the early stage. In contrast, the GWSA well 20.4 −19.2 ± 1.7 −17.7 ± 1.1 of Beijing eastern plain, Tianjin urban area and Cangzhou was used for Beijing, Tianjin and Hebei, respectively. There are missing records in parts of the years reflected the variations of subsidence with negative correlations during the whole period as it integrated the information of both discharge and recharge. It should be noted that the lack of long-term land subsidence data at large scale limits the understanding of the spatio-temporal variabilities of the aquifer system. The relationships derived from this study may vary as the scale changes; besides, the hysteresis of subsidence and lithological differences may also contribute to the spatio-temporal variations (Guo et al. 2015) . Seasonal variations revealed by satellite data Figure 7a plots the spatial pattern of mean land subsidence rate (2011–2015) from RADARSAT-2 InSAR data and Fig. 7b compares the mean annual cycle (2011–2015) of GWSA from in situ groundwater-level observations and InSAR-derived l a n d d e f o r m a t i o n i n t h e B e i j i n g p l a i n a r e a . T h e RADARSAT-2 InSAR data are of high spatial resolution for explaining the distribution and driving factors of land subsidence (Zhou et al. 2016) and provide much temporal information to reveal its seasonal variations. As shown in Fig 7b, the land surface deformation represents a similar seasonal cycle comparing to GWSA. Both GWS depletion and land surface deformation evidence accelerating rates from March to June, and recover in the second half year. It is believed that this seasonal cycle has resulted from intensive groundwater pumping during March–May when little precipitation is available for wheat growth, while precipitation during the rainy season (June–September) sustains the water demand for maize (Pan et al. 2017) . Fig. 7 a Spatial distribution of the mean subsidence rate in Beijing for the period 2011–2015. b Seasonal cycle of in situ GWSA and InSAR-derived land deformation during 2011 and 2015. The diagonal-patterned area (a) indicates the cone of depression, and its seasonal deformation is illustrated with the green-diamond curve (b) Conclusions This paper provides an overview of the long-term evolution of the GWS and land subsidence in the NCP, with a focus on the spatio-temporal variation of the aquifer systems responding to natural and anthropogenic factors. The GWS shows a prolonged declining rate of −17.8 ± 0.1 mm/yr during the whole study period 1971–2015, while more intensive depletion (−76.1 ± 6.5 mm/yr) can be found in Beijing during 1999–2012 due to a precipitation deficit. Large subsidence rates (>50 mm/yr) can generally be found in the late (2003– 2015), middle (1980–2000), and early (1971–1987) period for Beijing, Hebei, and Tianjin, respectively. The land subsidence showed varied correlations to GWSA and abstraction in space and time. Contrast correlations were found between subsidence and abstraction during early (positive) and late (negative) stages, while the GWSA represented a well indicator for subsidence due to its integrated representation of discharge and recharge. This study also demonstrates the potential and the limitations of remote sensing approaches (satellite gravimetry, InSAR) for understanding the spatial and temporal dynamics of the aquifer system. However, the data absence for some periods (mainly 1984–2002) and the lack of overall land surface deformation data at larger scale may lead to uncertainties associated with the correlation analysis. It is expected that the GRACE FollowOn mission as well as succeeding new SAR satellites can enhance the understanding of the regional aquifer system at multiple spatial and temporal scales. Funding Information This study is jointly funded by the Natural Science Foundation of China (41771456, 41771455) and Natural Science Foundation of Beijing (8152012). Open Access This article is distributed under the terms of the Creative C o m m o n s A t t r i b u t i o n 4 . 0 I n t e r n a t i o n a l L i c e n s e ( h t t p : / / creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Budhu M , Adiyaman I ( 2013 ) The influence of clay zones on land subsidence from groundwater pumping . Groundwater 51 ( 1 ): 51 - 57 . https://doi.org/10.1111/j.1745- 6584 . 2012 . 00943 .x Castellazzi P , Martel R , Galloway DL , Longuevergne L , Rivera A ( 2016a ) Assessing groundwater depletion and dynamics using GRACE and InSAR: potential and limitations . Groundwater 54 : 768 - 780 . https://doi.org/10.1111/gwat.12453 Castellazzi P , Martel R , Rivera A , Huang J , Pavlic G , Calderhead AI , Chaussard E , Garfias J , Salas J ( 2016b ) Groundwater depletion in central Mexico: use of GRACE and InSAR to support water resources management . Water Resour Res 52 : 5985 - 6003 . https:// doi.org/10.1002/2015WR018211 Chen B , Gong H , Li X , Lei K , Zhu L , Gao M , Zhou C ( 2017 ) Characterization and causes of land subsidence in Beijing, China . Int J Remote Sens 38 : 808 - 826 . https://doi.org/10.1080/01431161. 2016 .1259674 Chen J , Famiglietti JS , Scanlon BR , Rodell M ( 2016a ) Groundwater storage changes: present status from GRACE observations . Surv Geophys 37 : 397 - 417 . https://doi.org/10.1007/s10712-015-9332-4 Chen M , Tomás R , Li Z , Motagh M , Li T , Hu L , Gong H , Li X , Yu J , Gong X ( 2016b ) Imaging land subsidence induced by groundwater extraction in Beijing (China) using satellite radar interferometry . Remote Sens 8 : 468 . https://doi.org/10.3390/rs8060468 Chen Y , Chu L , Fan F , Zhang J , Lu Z , Shen R ( 2010 ) Present situation and genetic analysis of land subsidence in Hebei plain (in Chinese) . Hebei Geology Survey Annual Meeting 2010 , Beijing Famiglietti JS ( 2014 ) The global groundwater crisis . Nat Clim Chang 4 : 945 - 948 . https://doi.org/10.1038/nclimate2425 Famiglietti JS , Lo M , Ho SL , Bethune J , Anderson KJ , Syed TH , Swenson SC , de Linage CR , Rodell M ( 2011 ) Satellites measure recent rates of groundwater depletion in California's Central Valley . Geophys Res Lett 38 ( 3 ):L3403. https://doi.org/10.1029/ 2010GL046442 Feng W , Zhong M , Lemoine JM , Biancale R , Hsu HT , Xia J ( 2013 ) Evaluation of groundwater depletion in North China using the gravity recovery and climate experiment (GRACE) data and groundbased measurements . Water Resour Res 49 ( 4 ): 2110 - 2118 . https:// doi.org/10.1002/wrcr.20192 Guo S ( 2016 ) Study on the relationship among exploitation of deep confined groundwater, groundwater level and land subsidence in Tianjin plain (in Chinese) . MSc Thesis , China University of Geosciences, Beijing Guo H , Zhang Z, Cheng G , Li W , Li T , Jiao JJ ( 2015 ) Groundwaterderived land subsidence in the North China plain . Environ Earth Sci 74 : 1415 - 1427 . https://doi.org/10.1007/s12665-015-4131-2 Huang Z , Pan Y , Gong H , Yeh PJF , Li X , Zhou D , Zhao W ( 2015 ) Subregional-scale groundwater depletion detected by GRACE for both shallow and deep aquifers in North China plain . Geophys Res Lett 42 : 1791 - 1799 . https://doi.org/10.1002/2014GL062498 Hwang C , Yang Y , Kao R , Han J , Shum CK , Galloway DL , Sneed M , Hung W , Cheng Y , Li F ( 2016 ) Time-varying land subsidence detected by radar altimetry: California, Taiwan and North China . Sci Rep 6 : 28160 . https://doi.org/10.1038/srep28160 Jin S , Feng G ( 2013 ) Large-scale global groundwater variations from satellite gravimetry and hydrological models, 2002 - 2012 . Glob Planet Chang 106 : 20 - 30 . https://doi.org/10.1016/j.gloplacha. 2013 . 02 .008 Liu M , Kou W , Luo Y , Wang R , Tian F , Zhao C , Liu J ( 2016 ) Analysis of the relationship between groundwater exploitation and land subsidence in Beijing city (in Chinese) . Geology 11 : 21 - 25 Long D , Chen X , Scanlon BR , Wada Y , Hong Y , Singh VP , Chen Y , Wang C , Han Z , Yang W ( 2016 ) Have GRACE satellites overestimated groundwater depletion in the Northwest India aquifer? Sci Rep 6 : 24398 . https://doi.org/10.1038/srep24398 Pan Y , Pan J , Gong H , Zhao W ( 2004 ) Research on the relation between groundwater exploitation and subsidence in Tianjin urban area (in Chinese) . Earth Environ 32 : 36 - 39 Pan Y , Zhang C , Gong H , Yeh PJF , Shen Y , Guo Y , Huang Z , Li X ( 2017 ) Detection of human-induced evapotranspiration using GRACE satellite observations in the Haihe River basin of China . Geophys Res Lett 44 : 190 - 199 . https://doi.org/10.1002/2016GL071287 Richey AS , Thomas BF , Lo M , Reager JT , Famiglietti JS , Voss K , Swenson S , Rodell M ( 2015 ) Quantifying renewable groundwater stress with GRACE . Water Resour Res 51 : 5217 - 5238 . https://doi. org/10.1002/2015WR017349 Rodell M , Houser PR , Jambor UEA , Gottschalck J , Mitchell K , Meng CJ , Arsenault K , Cosgrove B , Radakovich J , Bosilovich M , Entin JK , Walker JP , Lohmann D , Toll D ( 2004 ) The global land data assimilation system . Bull Am Meteor Soc 85 ( 3 ): 381 - 394 . https://doi.org/ 10.1175/BAMS-85-3-381 Rodell M , Velicogna I , Famiglietti JS ( 2009 ) Satellite-based estimates of groundwater depletion in India . Nature 460 : 999 - 1002 . https://doi. org/10.1038/nature08238 Save H , Bettadpur S , Tapley BD ( 2016 ) High-resolution CSR GRACE RL05 mascons . J Geophys Res Solid Earth 121 : 7547 - 7569 . https:// doi.org/10.1002/2016JB013007 Scanlon BR , Longuevergne L , Long D ( 2012 ) Ground referencing GRACE satellite estimates of groundwater storage changes in the California Central Valley , USA. Water Resour Res 48 : 1TW045201T. https://doi.org/10.1029/2011WR011312 Scanlon BR , Zhang Z , Save H , Wiese DN , Landerer FW , Long D , Longuevergne L , Chen J ( 2016 ) Global evaluation of new GRACE mascon products for hydrologic applications . Water R e s o u r R e s 5 2 : 9 4 1 2 - 9 4 2 9 . h t t p s : / / d o i . o r g / 1 0 . 1 0 0 2 / 2016WR019494 Wada Y ( 2016 ) Modeling groundwater depletion at regional and global scales: present state and future prospects . Surv Geophys 37 : 419 - 451 . https://doi.org/10.1007/s10712-015-9347-x Wang Y ( 2004 ) Study on land subsidence caused by overexploiting groundwater in Beijing (in Chinese) . Site Invest Sci Technol 5 : 46 - 49 Watkins MM , Wiese DN , Yuan D , Boening C , Landerer FW ( 2015 ) Improved methods for observing Earth's time variable mass distribution with GRACE using spherical cap mascons . J Geophys Res Solid Earth 120 : 2648 - 2671 . https://doi.org/10.1002/2014JB011547 Yan HM , Zhong M , Zhu YZ ( 2006 ) Scientific computation library WHIGG F90 LIB and it's application in geosciences (in Chinese) . J Geodesy Geodynam 26 ( 1 ): P109 - P114 Ye S , Xue Y , Wu J , Yan X , Yu J ( 2016 ) Progression and mitigation of land subsidence in China . Hydrogeol J 24 : 685 - 693 . https://doi.org/10. 1007/s10040-015-1356-9 Zhang Y , Wu H , Kang Y , Zhu C ( 2016 ) Ground subsidence in the Beijing-Tianjin-Hebei region from 1992 to 2014 revealed by multiple SAR stacks . Remote Sens 8 : 675 . https://doi.org/10.3390/ rs8080675 Zhang Z , Fei Y , Chen Z et al ( 2009 ) Survey and evaluation of groundwater sustainable utilization in North China plain (in Chinese) . Geological House , Beijing Zhou C , Gong H , Chen B , Zhu F , Duan G , Gao M , Lu W ( 2016 ) Land subsidence under different land use in the eastern Beijing plain, China 20052013 revealed by InSAR timeseries analysis . GISci Remote Sens 53 : 671 - 688 . https://doi.org/10.1080/15481603. 2016 . 1227297 Zhu L , Gong H , Li X , Wang R , Chen B , Dai Z , Teatini P ( 2015 ) Land subsidence due to groundwater withdrawal in the northern Beijing plain , China. Eng Geol 193 : 243 - 255 . https://doi.org/10.1016/j. enggeo. 2015 . 04 .020


This is a preview of a remote PDF: https://link.springer.com/content/pdf/10.1007%2Fs10040-018-1768-4.pdf

Huili Gong, Yun Pan, Longqun Zheng, Xiaojuan Li, Lin Zhu, Chong Zhang, Zhiyong Huang, Zhiping Li, Haigang Wang, Chaofan Zhou. Long-term groundwater storage changes and land subsidence development in the North China Plain (1971–2015), Hydrogeology Journal, 2018, 1-11, DOI: 10.1007/s10040-018-1768-4