Long-term groundwater storage changes and land subsidence development in the North China Plain (1971–2015)
August 2018, Volume 26, Issue 5, pp 1417–1427 | Cite as
Long-term groundwater storage changes and land subsidence development in the North China Plain (1971–2015)
AuthorsAuthors and affiliations
Huili GongYun PanLongqun ZhengXiaojuan LiLin ZhuChong ZhangZhiyong HuangZhiping LiHaigang WangChaofan Zhou
First Online: 16 April 2018
Received: 19 September 2017
Accepted: 30 March 2018
1 Shares 1.1k Downloads
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.
KeywordsOver-abstraction Subsidence Drought GRACE China
Published in the special issue “Groundwater sustainability in fast-developing China”
Variations à long terme du stockage des eaux souterraines et développement de l’affaissement des terrains dans la Plaine de Chine du Nord (1971–2015)
La Plaine de Chine du Nord (PCN) souffre d’une diminution en ressources en eaux souterraines (RES) et d’affaissement des terrains depuis longtemps. Cet article recueille des données sur les changements de stockage des eaux souterraines et les affaissements des terrains à partir de mesures in situ des niveaux piézométriques, de la littérature et d’observations satellitaires afin de fournir un aperçu de l’évolution du système aquifère au cours de la période 1971–2015 en mettant l’accent sur les variations subrégionales. On a constaté que les RES présentaient un taux de déclin prolongé de –17.8 ± 0.1 mm/an pour la période 1971–2015, avec une corrélation négative avec les volumes d’eau souterraine prélevés avant ~2000 et une corrélation positive après ~2000. Les corrélations statistiques entre le taux de subsidence et l’anomalie RES (ARES) , les prélèvements en eaux souterraines, et les précipitations annuelles indiquent que l’affaissement des terrains dans trois sous-régions (Beijing, Tianjin et Hebei) représente différentes variations temporelles dues à des facteurs variés. La sécheresse sans interruption a causé une intense diminution des RES (–76.1 ± 6.5 mm/an) et un affaissement des terraines à Beijing entre 1999 et 2012. Des corrélations négatives entre les prélèvements totaux en eaux souterraines et l’affaissement des terrains observées après les années 1980 indiquent qu’il est peut-être discutable de déduire la subsidence à partir des données régionales de prélèvements. Au lieu de cela, la ARES fournit généralement une corrélation fiable avec la subsidence. Cette étude met en évidence les variabilités spatio-temporelles de la diminution des RES et de l’affaissement des terrains dans la PNC sous contraintes naturelles et anthropiques, et l’importance des changements de RES pour la compréhension du développement de l’affaissement des terrains ou subsidence.
Cambios en el almacenamiento de agua subterránea y en el desarrollo de subsidencia del terreno a largo plazo en la Llanura Norte de China (1971–2015)
La Llanura Norte de China (NCP) ha estado sufriendo un agotamiento del almacenamiento de agua subterránea (GWS) y la subsidencia del terreno durante un largo período. Este trabajo recopila datos sobre los cambios de GWS y la subsidencia del terreno a partir de mediciones in situ del nivel del agua subterránea, bibliografía y observaciones satelitales para proporcionar una visión general de la evolución del sistema acuífero durante 1971–2015 con un enfoque en las variaciones subregionales. Se encontró que el GWS mostró una tasa decreciente prolongada de −17.8 ± 0.1 mm/año durante 1971–2015, con una correlación negativa con la extracción de aguas subterráneas antes del año ~2000 y una correlación positiva después del ~2000. Las correlaciones estadísticas entre la tasa de subsidencia y la anomalía de GWS (GWSA), la extracción de agua subterránea y la precipitación anual muestran que la subsidencia del terreno en tres subregiones (Beijing, Tianjin y Hebei) representa variaciones temporales diferentes debido a factores variables de las causas. La sequía continua causó el agotamiento intensivo de GS (−76.1 ± 6.5 mm/año) y la subsidencia del terreno en Beijing durante 1999–2012. Las correlaciones negativas entre la extracción total del agua subterránea y la subsidencia exhibida después de la década de 1980 indican que puede ser cuestionable inferir la subsidencia de los datos de extracción regional. En cambio, el GWSA generalmente proporciona una correlación confiable con la subsidencia. Este estudio destaca las variabilidades espacio-temporales del agotamiento de GWS y la subsidencia del terreno en el NCP bajo impactos naturales y antropogénicos, y la importancia de los cambios de GWS para entender el desarrollo de la subsidencia del terreno.
地下水消耗与地面沉降是困扰华北平原很长时间的问题。本文综合水井观测、文献、卫星遥感手段获取的地下水储量变化和地面沉降数据,对华北平原含水层系统的长期变化(1971–2015)进行了分析,重点关注其演化的时空特征。研究发现,华北平原1971–2015年间的地下水储量变化速率为–17.8 ± 0.1 mm/yr。总体上,地下水储量距平变化与地下水年开采量在2000年以前为负相关,在2000年以后为正相关。地面沉降速率与其影响因素(地下水储量距平变化、年开采量、年降水量)在北京、天津、河北三个子区上表现出不同的统计特征,反映了京津冀地面沉降演化机制的差异性。1999–2012年连续的干旱使得北京地下水迅速消耗(−76.1 ± 6.5 mm/yr),并伴随着快速发展的地面沉降。研究同时发现,1980年代后地面沉降量与地下水开采量之间变为负相关,而地下水储量距平变化却总体上始终与沉降量表现为正相关。此次研究再次凸显了,在气候变化与人类活动共同作用下,华北平原地下水消耗与地面沉降的时空差异性,以及地下水储量变化对于更好地认识地面沉降演化的重要性。
Mudanças no armazenamento das águas subterrâneas de longo termo e desenvolvimento de subsidência de terreno na Planície do Norte da China (1971–2015)
A Planície do Norte da China (PNC) tem sofrido o esgotamento de armazenamento de águas subterrâneas (AAS) e a subsidência de terreno por um longo período. Este artigo coleta dados sobre as mudanças no AAS e a subsidência de terreni a partir de medições, literatura e observações de satélites no nível da água subterrânea in situ, para fornecer uma visão geral da evolução do sistema aquífero durante 1971–2015 com foco nas variações sub-regionais. Verificou-se que o AAS mostrou uma taxa de declínio prolongada de −17.8 ± 0.1 mm/ano durante 1971–2015, com uma correlação negativa com a captação de águas subterrâneas antes do ano 2000 e uma correlação positiva após ~2000. As correlações estatísticas entre a taxa de subsidência e a anomalia no ASS (AAAS), a captação de águas subterrâneas e as precipitações anuais mostram que a subsidência de terreno em três sub-regiões (Pequim, Tianjin e Hebei) representa diferentes variações temporais devido a fatores forçantes variáveis. A seca contínua causou o esgotamento intensivo de GWS (−76.1 ± 6.5 mm/ano) e a subsidência de terreno em Pequim durante 1999–2012. As correlações negativas entre a captação total de águas subterrâneas e a subsidência de terreno exibida após a década de 1980 indicam que pode ser questionável inferir subsidência a partir de dados de abstração regionais. Em vez disso, a AAAS geralmente fornece uma correlação confiável com a subsidência. Este estudo destaca as variabilidades espaço-temporais da depleção de AAS e da subsidência da terra na PNC sob impactos naturais e antropogênicos e a importância das mudanças de AAS para entender o desenvolvimento de subsidência de terreno.
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 well-known aquifers worldwide (Chen et al. 2016a; Famiglietti et al. 2011; Rodell et al. 2009); however, groundwater over-abstraction 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.
Open image in new window
Location map of the study area and the in situ observation gauges, and spatial map of specific yield according to Zhang et al. (2009)
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 Capacity Model) from the Global Land Data 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 sites (Fig. 1) provided by China Meteorological 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.
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).
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.
Open image in new window
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
Fitted seasonal terms (annual and semi-annual amplitudes and phases) and long-term trends of the monthly GWSA time series (Fig. 2) derived from GRACE and in situ observations based on the least-squares harmonic analysis (Jin and Feng 2013; Yan et al. 2006)
Annual amplitude (mm)
Semi-annual amplitude (mm)
Annual phase (degree)
Semi-annual phase (degree)
−19.2 ± 1.7
−17.7 ± 1.1
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.
Open image in new window
The 1971–2013 yearly time series of precipitation, groundwater abstraction, land subsidence, and GWSA in three sub-regions of a Beijing, b Tianjin, and c Hebei. Land subsidence data of typical regions 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
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).
Open image in new window
Statistical relationships between the GWSA and a–c precipitation, d–f annual groundwater abstraction, and g–i cumulative groundwater abstraction for Beijing (left panel), Tianjin (middle panel) and Hebei (right panel). The shaded area indicates the 95% confidence interval
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 GW-level 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.
Open image in new window
The graduated color map delineates the spatial distribution of GRACE-derived GWSA trend (2003–2015) and the black contour lines show the accumulated land subsidence from 1971 to 2003 in NCP
Statistics of annual precipitation (P), groundwater (GW) abstraction, mean subsidence rate, and GWSA trend in the three sub-regions of Beijing, Tianjin, and Hebei during 1971–2013
Mean annual P (mm)
Mean annual GW abstraction (mm)
Mean subsidence rate (mm/yr)
GWSA trend (mm/yr)
−36.9 ± 5.1
3.9 ± 2.8
−76.1 ± 6.5
−34.1 ± 2.4
−3.1 ± 1.3
−20.2 ± 2.0
−7.0 ± 0.4
−14.7 ± 1.1
−19.8 ± 2.1
−20.7 ± 3.4
−15.2 ± 0.6
Open image in new window
Statistical relationships between annual subsidence and a–c precipitation, d–f groundwater abstraction, and g–i GWSA, for Beijing (left panel), Tianjin (middle panel) and Hebei (right panel). The shaded area indicates the 95% confidence interval
The different correlations represented by abstraction and GWSA, as well as their variations in space and time, 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 changes 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 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 land deformation in the Beijing plain area. The 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).
Open image in new window
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)
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 Follow-On mission as well as succeeding new SAR satellites can enhance the understanding of the regional aquifer system at multiple spatial and temporal scales.
This study is jointly funded by the Natural Science Foundation of China (41771456, 41771455) and Natural Science Foundation of Beijing (8152012).
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 CrossRefGoogle Scholar
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 CrossRefGoogle Scholar
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 CrossRefGoogle Scholar
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 CrossRefGoogle Scholar
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 CrossRefGoogle Scholar
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 CrossRefGoogle Scholar
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, BeijingGoogle Scholar
Famiglietti JS (2014) The global groundwater crisis. Nat Clim Chang 4:945–948. https://doi.org/10.1038/nclimate2425 CrossRefGoogle Scholar
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 CrossRefGoogle Scholar
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 ground-based measurements. Water Resour Res 49(4):2110–2118. https://doi.org/10.1002/wrcr.20192 CrossRefGoogle Scholar
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, BeijingGoogle Scholar
Guo H, Zhang Z, Cheng G, Li W, Li T, Jiao JJ (2015) Groundwater-derived land subsidence in the North China plain. Environ Earth Sci 74:1415–1427. https://doi.org/10.1007/s12665-015-4131-2 CrossRefGoogle Scholar
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 CrossRefGoogle Scholar
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 CrossRefGoogle Scholar
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 CrossRefGoogle Scholar
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–25Google Scholar
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 CrossRefGoogle Scholar
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–39Google Scholar
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 CrossRefGoogle Scholar
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 CrossRefGoogle Scholar
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 CrossRefGoogle Scholar
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 CrossRefGoogle Scholar
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 CrossRefGoogle Scholar
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 Resour Res 52:9412–9429. https://doi.org/10.1002/2016WR019494 CrossRefGoogle Scholar
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 CrossRefGoogle Scholar
Wang Y (2004) Study on land subsidence caused by overexploiting groundwater in Beijing (in Chinese). Site Invest Sci Technol 5:46–49Google Scholar
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 CrossRefGoogle Scholar
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–P114Google Scholar
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 CrossRefGoogle Scholar
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 CrossRefGoogle Scholar
Zhang Z, Fei Y, Chen Z et al (2009) Survey and evaluation of groundwater sustainable utilization in North China plain (in Chinese). Geological House, BeijingGoogle Scholar
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 CrossRefGoogle Scholar
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 CrossRefGoogle Scholar
© The Author(s) 2018
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://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.
Authors and Affiliations
Huili Gong12Yun Pan12Email authorLongqun Zheng12Xiaojuan Li12Lin Zhu12Chong Zhang34Zhiyong Huang5Zhiping Li6Haigang Wang7Chaofan Zhou121.Beijing Laboratory of Water Resources SecurityCapital Normal UniversityBeijingChina2.College of Resource Environment and TourismCapital Normal UniversityBeijingChina3.State Key Laboratory of Earth Surface Processes and Resource Ecology, Faculty of Geographical ScienceBeijing Normal UniversityBeijingChina4.Institute of Land Surface System and Sustainable Development, Faculty of Geographical ScienceBeijing Normal UniversityBeijingChina5.Department of Earth SciencesThe University of Hong KongHong KongChina6.Beijing Institute of Hydrogeology and Engineering GeologyBeijingChina7.China Geological and Environmental Monitoring InstituteBeijingChina