Analysis of Spatiotemporal Variation of Land Subsidence in Beijing Plain , China

Since the 1970s, land subsidence has been developing rapidly in the Beijing Plain, the systematic study of its evolution mechanism is of great significance to the sustainable development of the regional economy. First, based on ENVISAT ASAT and RADARSAT2 data, the land subsidence data in Beijing Plain were obtained using permanent interferometer technology. Second, based on the GIS platform and using fishing net tools, vector data of ground settlement with different resolutions were obtained. Through a series of tests, a scale of 960 metres was selected as the research unit, and the subsidence rate of the grid was obtained from 2004 to 2015. Finally, based on the Mann-Kendall mutation test method, a trend analysis of land subsidence changes in various grids was carried out. The results showed that single-year mutation mainly distributed in the middle and lower parts of the Yongding River alluvial fan and the Chaobai River alluvial fan, mainly occurring in 2015, 2005 and 2013, respectively. The upper and middle alluvial fan of the Chaobai River, the vicinity of the emergency water source and the edge velocity of the groundwater funnel have undergone several sudden changes. Combined with hydrogeology, basic geological conditions and the impact of the South-to-North Water transfer project, we analysed the causes of the mutations in the grid. The research results can provide a basis for the study and prevention of land subsidence in this area and help to further explore the trend characteristics of land subsidence in this area.


nce in the st
dy area.

Introduction

Land subsidence is a geological phenomenon that is caused by human engineering activities (e.g., groundwater exploitation) or natural factors, and, in this phenomenon, the ground elevation within a certain area decreases [1,2].The occurrence of land subsidence is nearly irreversible.When severe local surface subsidence occurs, it can induce a series of geological environmental disasters, such as foundation sinking, house cracking, and underground pipeline damage [3][4][5][6]; furthermore, subsidence is characterized as necessitating a long formation time, having a wide influential range, and being highly difficult to prevent and control [2].At present, more than 150 countries and regions worldwide have suffered from land subsidence, among which the most severe areas are central Mexico [7,8], central Iran [9,10], the Bandung basin of Indonesia [11], and northern Italy [12,13].In China, land subsidence has mainly occurred on the North China Plain and within the Yangtze River Delta, the Fenwei Basin, and the Pearl River Delta [6,[14][15][16][17][18][19].Among these areas, land subsidence on the North China Plain (the Beijing-Tianjin-Hebei region) continues to increase, the subsidence within the Yangtze River Delta region has been effectively controlled, and the subsidence within the Fenwei Basin region is still rapidly developing.

In 2000, Persistent Scatterers Interferometry Synthetic Aperture Radar (PS-InSAR), which was proposed by Ferretti, et al. [20], was used to effectively reduce the influence of phase incoherence and atmospheric delay in traditional differential interferometry SAR (D-InSAR).In 2002 and 2004, Berardino, et al. [21] and Lanari,et al. [22] proposed the small baseline set interferometry (SBAS-InSAR) te hnique, which is more suitable for the slow deformation of surface monitoring over the long term.Many researchers have used this method to carry out in-depth research.In recent years, related research has mainly focused on severe land subsidence areas.In 2014, Chaussard, Wdowinski, Cabral-Cano, and Amelung [1] used ALOS data to monitor land subsidence in central Mexico and confirmed that groundwater overdraft was the main cause of land subsidence in that region.In 2014, Qu, et al. [23] used ENVISAT ASAR, ALOS, and TerraSAR-X data to reveal the temporal and spatial variation characteristics of ground fractures and surface subsidence in Xi'an, China.In 2016, Amighpey and Arabi [9] used ENVISAT ASAR data to obtain surface deformation information of the Yazd-Ardakan Plain in central Iran and compared it with changes in water levels to estimate the skeleton water storage coefficient.In 2017, Shviro, et al. [24] used COSMO-SkyMed data to study the quantitative impacts of floods on land subsidence in the Dead Sea region.In 2018, Zhao, et al. [25] used ALOS data to reveal the surface deformations in the Linfen-Yuncheng (China) basin and successfully detected the active faults and ground fractures.

As the capital of China, Beijing resides on one of the most severe areas of land affected by subsidence worldwide.Beijing has been experiencing rapid urban development and populatio

growth sinc
the 1960s, and its demand for water resources continues to increase [19].Over the years, groundwater has accounted for two-thirds of the total water supply in Beijing [26].Moreover, Beijing is one of the few global metropolises whose main source of water supply is groundwater [27].Land subsidence has been increasing rapidly in the Beijing area, with the long-term overexploitation of groundwater, and the land subsidence funnels have tended to be interconnected [28].At the end of 2015, the maximum cumulative settlement of the Beijing Plain area was approximately 1.4 metres, and the area whose cumulative settlement exceeded one metre reached 4000 km 2 .Thus, land subsidence has become a major threat to Beijing's urban construction and residential safety [29].

Using InSAR technology, many researchers have conducted in-depth studies on land subsidence in Beijing.In 2009, Gong,et al. [30] used ENVISAT ASAR data from 2003 to 2006 to monitor land subsidence on the Beijing Plain via Persistent Scatterers Interferometry (PSI).The results showed that Quaternary faults have certain influence on the spatial pattern of ground subsidence, which often occurs in areas whose clay layer is thicker than 50 m.In 2015, Chen, et al. [31] used ENVISAT ASAR data to obtain land subsidence information regarding Beijing, and, by combining GIS spatial analysis technology, they analysed the spatial and temporal evolutionary characteristics of uneven settlements.In 2015, Zhu, et al. [19] combined ENVISAT ASAR data, LANDSAT data, and hydrogeological data to study the land subsidence on the northern plain of Beijing.The results showed that the maximum subsidence rate reached 52 mm/year and that the silty clay layer strongly contributed to land subsidence.In 2016, Chen, et al. [32] used TerraSAR-X and ENVISAT ASAR data to analyse the spatial and temporal distribution characteristics of land subsidence on the Beijing Plain and revealed the main inducing factors.Chen, et al. [33] subsequently used the RADARSAT-2 data to reveal the spatial and temporal variation characteristics of the subsidence of Beijing Metro Line 6.In 2017, Zhou, et al. [31] used TerraSAR-X data to monitor the surface deformation information concerning the eastern Beijing Plain from 2010 to 2015 and analysed the correlations between the different land use types and land subsidence.In addition, Zhou,et al. [34] used TerraSAR-X and ENVISAT ASAR data to determine the surface deformation characteristics of the Beijing Plain.The authors determined that the maximum expansion direction of land subsidence on the Beijing Plain was East and that the growth rate was 11.3 km 2 /year using the equal-sector analysis method.In 2018, using TerraSAR-X, ENVISAT ASAR, and urban building information data, Yang, et al. [35] analysed the relationships between land subsidence and building density at the regional, block, and building scales.In 2019, Chen, et al. [36] used RADARSAT-2 data, SBAS technology, and wavelet change methodology to analyse the lag time between ground subsidence and groundwater levels quantitatively on the Beijing Plain.Relevant studies have mainly focused on the spatial distribution characteristics and formation mechanisms of land subsidence in Beijing; however, few studies have investigated the variation characteristics of land subsidence in conjunction with time series.

At present, the problem of land subsidence mainly concerns its spatial distribution and formation factors, and temporal issues that are associated with mutations of land subsidence have never been reported.Studying the mutation time and its corresponding region can help to quantitatively reveal the temporal and spatial variation characteristics of land subsidence.The Mann-Kendall method can be used to quantitatively analyse the mutation information of time series of precipitation, temperature, water quality, and other factors [37][38][39].This method has never been applied to land subsidence.We can use this method to study the mutation information of land subsidence in conjunction with time series quantitatively to further grasp the spatial-temporal variation characteristics of land subsidence in a study area.Therefore, this study used ENVISAT ASAR and RADARSAT-2 to obtain the surface deformation information of the Beijing Plain area from 2004 to 2015; furthermore, the research area was divided into 5575 grids that were 960 m by 960 m via the Fishnet tool.A Mann-Kendall test was subsequently carried out for each grid via Python, and the land subsidence mutation information of the Beijing Plain area was obtained on the basis of the grids.The combination of hydrogeology and basic geological conditions ultimately resulted in an analysis of the causes of mutations.


Study Area

As the capital of China, Beijing is the centre of the national politics, economy, and culture, and it is an international metropolis with a population of nearly 20 million people.The area of Beijing is located on the north-western edge of the North China Plain and it occupies the region between 39.4 • N, 115.7 • E and 41.6 • N, 117.4 • E. In general, the topography is higher in the north-western region and lower in the sout

ogic features
are divided into three parts: the western mountainous area, the northern mountainous area, and the south-eastern plain (Figure 1).The total area of Beijing is 1,680,780 km 2 and the plain covers 38% of the total area.From 2000 to 2011, Beijing suffered a continual drought, with an average rainfall of only 459.11 mm.The Beijing Plain is formed by the combination of five major water systems (i.e., the Juma River system, Yongding River system, Wenyu River system, Chaobai River system, and Jiyunhe River system) [3].As shown in Figure 2, the Quaternary loose sediments in the study area are widely distributed.At the top of the alluvial fan, the Quaternary thickness is approximately 20~40 m; a single gravel layer or a thin viscous soil layer that is covered with gravel represents this thickness [19].In the lower part of the alluvial fan, the Quaternary loose sediment thickness gradually increases, the gradation increases, and the particles gradually become smaller.The lithology gradually transitions to a point where sand, gravel, and viscous soil intersect with each other and viscous soil dominates it [36].There are many fractur

in the study area, and these fractures tend to
ostly occur in the SE-NW and NE-SW directions.The Beijing Plain is formed by the combination of five major water systems (i.e., the Juma River system, Yongding River system, Wenyu River system, Chaobai River system, and Jiyunhe River system) [3].As shown in Figure 2, the Quaternary loose sediments in the study area are widely distributed.At the top of the alluvial fan, the Quaternary thickness is approximately 20~40 m; a single gravel layer or a thin viscous soil layer that is covered with gravel represents this thickness [19].In the lower part of the alluvial fan, the Quaternary loose sediment thickness gradually increases, the gradation increases, and the particles gradually become smaller.The lithology gradually transitions to a point where sand, gravel, and viscous soil intersect with each other and viscous soil dominates it [36].There are many fractures in the study area, and these fractures tend to mostly occur in the SE-NW and NE-SW directions.Land subsidence, as a surface response to the development, utilization, and evolution of underground space, has become a global and multidisciplinary complex geological environmental problem [29].Land subsidence first occurred in Beijing in 1935, and the main subsidence area was located in the area from Xidan to Dongdan.Since the 1970s, land subsidenc

on the Beijing Plain has entered a period o
rapid development.In general, the subsidence area can be divided into two major areas (the northern and the southern areas) and it has seven subsidence centres [27].The northern subsidence area is large and it includes six centres: Beiqijia in Changping; Xixiaoying in Haidian; JinZhan, Sanjianfang, Heizhuanghu in Chaoyang; and, urban Tongzhou.The southern subsidence area is small and it includes the Yufa subsidence centre in Daxing.By 2012, the cumulative settlement greater than 50 mm reached 4,323 km 2 , the maximum historical settlement rate was 159.6 mm/year and the maximum cumulative settlement (1955-2012) was 1414 mm (Baxianzhuang Settlement Center) [29].At of the end of 2015, the area of land subsidence in Beijing involved an area of 4,842 km 2 , with an average annual subsidence of 21.6 mm/year and an annual maximum value of nearly 141 mm/year [36].Land subsidence, as a surface response to the development, utilization, and evolution of underground space, has become a global and multidisciplinary complex geological environmental problem [29]

and subsidence first occurr
d in Beijing in 1935, and the main subsidence area was located in the area from Xidan to Dongdan.Since the 1970s, land subsidence on the Beijing Plain has entered a period of rapid development.In general, the subsidence area can be divided into two major areas (the northern and the southern areas) and it has seven subsidence centres [27].The northern subsidence area is large and it includes six centres: Beiqijia in Changping; Xixiaoying in Haidian; JinZhan, Sanjianfang, Heizhuanghu in Chaoyang; and, urban Tongzhou.The southern subsidence area is small and it includes the Yufa subsidence centr in Daxing.By 2012, the cumulative settlement greater than 50 mm r mm/year and the maximum cumulative settlement (1955-2012) was 1414 mm (Baxianzhuang Settlement Center) [29].At of the end of 2015, the area of land subsidence in Beijing involved an of 21.6 mm/year and an annual maximum value of nearly 141 mm/year [36].


Methodology


Data Sources

Remote Sens. 2019, 11, 1170 6 of 20


Methodology


Data Sources

In this study, ENVISAT ASAR data and RADARSAT-2 data were selected to obtain the sur December 2015.ENVISAT is a solar-synchronous polar orbiting satellite that was launched by the European Space Agency on 1 March 2002, and it has been out of contact since 8 April 2012; additionally, the largest sensor on it is synthetic aperture radar.The ENVISAT ASAR sensor has an orbital altitude of 799.8 km, its operating band is C-band with a wavelength of 5.6 cm, and its orbital repeat observation interval is 35 days.We employed 47 images of the ascending-o

ta from June 2003 to August 2010, and the images covered m
st of the Beijing Plain, except for the southernmost tip of Beijing.RADARSAT-2 is a radar satellite launched on 14 December 2007, by the Canadian Space Agency in collaboration with the MacDonald Dettwiler Associates (MDA).This satellite operates in the C-band with a wavelength of 5.6 cm and the orbital repeat observation interval is 24 days.We employed 48 RADARSAT-2 images from the descending-orbit data from November 2010 to October 2015 (Table 1).


Processing of SBAS-InSAR on the Basis of StaMPS

Small baseline interferometry (SBAS-InSAR), as originally proposed by Berardino and Lanari in 2002 and 2004, is a reliable method for studying long-term slow surface deformation [21,22].The principle of this method is that the SAR image in the same study area is formed into several baseline pairs on the basis of the sizes of the spatial baseline and time baseline.Thus, the small baseline data are combined into several sets.Therefore, within a set, the baseline of the interference pair is small; between the sets, the interference pair is relatively large.Hooper initially proposed the method adopted in this paper in 2008, and this method differs by the identification of single-view slowly decorrelated filtered phase (SDFP) pixels in a single-view image [37].In this study, on the basis of Delft object-oriented radar interferometric software (DORIS), we deduced that the time bas

ine threshold and
the space baseline threshold of the ASAR images were 500 days and 500 metres, respectively, and the interference processing generated 46 interference images.The external Digital Elevation Model (DEM) that was applied to remove the topographic phase and flatten any effects that originated from the Shuttle Radar Topography Mission (SRTM), with a spatial resolution of 90 m (ht p://dds.cr.usgs.gov/srtm/).


Quasi-PS InSAR Processing Based on SARPROZ

Unlike traditional PS-InSAR, which can only process stable Persistent Scatter (PS) points, Quasi-PS InSAR technology (QPS) is an extension of PS-InSAR.The QPS method employs a spatial filtering method to improve the signal-to-noise ratio of the high-point target and selects interference pairs with high coherence in the time series data set to extract the elevation and deformation rate information of the target to involve more points in modelling and extract high-precision deformation information [38].

In the present study, the QPS technical method in SARPROZ software was selected to process the 48 RADARSAT-2 images.The RADARSAT-2 data resolution is 30 m and the time span ranged from 22 November 2010 to 20 November 2015.The time baselines and spatial baselines of the RADARSAT-2 images during 2010-2015 were 500 days and 500 metres, respectively.A total of 266 interferometric image pairs were generated during processing.The external DEM data were also retrieved from the

RM3 data (spatial
resolution of 90 m).


Mann-Kendall Mutation Test

Also known as the non-distribution test, the Mann-Kendall test is a nonparametric statistical method.The advantage of this test is that it does not require samples to obey a certain distribution, and a few outliers do not affect the results [40].The Mann-Kendall method has been used by many researchers to analyse the trend hanges of various elements, such as precipitation, runoff, and temperature, in conjunction with time series [37][38][39].In this paper, this test method was used to evaluate the change trend of land subsidence on the Beijing Plain.The method is summarized, as follows:

For time series X with n samples, an ordered list is constructed:
S k = k i=0 r i c r i = 1, x i > x j 0, else j = 1, 2, . . . , i(1)
S k repre

nts the cumulative number of values greater at time i tha
at time j.Under the assumption of stochastic independence of time series, we define the following statistics:
UF K = S k − E(S k ) Var(S k ) k = 1, 2, . . . , n(2)
Among them, UF 1 = 0 and E(S k ) and Var(S k ) are the mean and variance of the cumulative S k , respectively.When X 1 , X 2 , ..., X n are independent of each other and they present the same continuous distribution, the following formula is obtained:
E(S k ) = k(k − 1) 4 Var(s k ) = k(k − 1)(2k + 5) 72(3)
where UF k is a standard normal distribution.A statistical sequence is calculated in the order of time series X (x 1 , x 2 , ..., x n ).The significance level is given as α, after which the normal distribution table can be queried.If |UF k | > Uα the sequence exhibits an obvious trend change.Similarly, the inverse order of UFk is calculated as UBk.The significance level that was selected in this study was α = 0.05 and the critical value was U 0.025 = ±1.96.


Results


Acquisition of Ground Deformation Time Series Information

On the basis of the surface deformation information that was monitored by SBAS-InSAR and QPS, a displacement rate distribution map was obtained for the Beijing Plain.As shown in Figure 3, the spatial distribution of land subsidence on the Beijing Plain is widespread, and the spatial difference is dramatic.The areas with severe settlement are mainly in the eastern region of Chaoyang, the north-western region of Tongzhou, the southern region of Changping, the north-western region of Shunyi, and the southern region of Daxing.From 2004 to 2010, the number of PS points that were obtained on the basis of the ENVISAT ASAR data was 266,581 and the average deformation rate ranged from −114.4 mm/year to +18.2 mm/year.Among the PS points, the area whose settlement rate was greater than 25 mm/year reached 1078.5 km 2 , which ac ounted for 17.2% of the total area of the Beijing Plain.From 2011 to 2015, 100,515 PS points were obtained on the basis of the RADARSAT-2 data, and the average deformation rate ranged from −133.6 mm/year to +18.4 mm/year.The area whose settlement rate was greater than 25 mm/year reached 1139 km 2 , which accounted for 17.8% of the area of the Beijing Plain.By the end of 2015, the area whose cumulative settlement across the Beijing Plain exceeded 300 mm reached 1271 km 2 , which accounted for approximately 20% of the tot

alidation
his work was validated in two steps.First, the second-class level deformation isoline of 2007 was selected and the results of SBAS-InSAR were preliminarily tested as a whole.Figure 4 shows that the settlement of the PS points that were obtained by SBAS-InSAR in 2007 were sufficiently consistent with that of the level measurements across the overall spatial distribution.

Furthermore, the data of 25 levelling benchmarks from 2003 to 2013 were selected for validation.The levelling benchmarks were taken as the original points, and all of the monitoring points within a radius of 150 m were extracted.The mean value of the extracted PS points was taken as the settlement estimation value and then verified according to the settlement of the levelling benchmarks.As shown in Figure 5, the correlation coefficient of the InSAR monitoring results and the level monitoring results from 2003 to 2010 was 0.95, the maximum absolute error was 17.9 mm/year, the minimum absolute error was 0.1 mm/year, and the standard deviation was 4 mm/year.From 2011 to 2013, the correlation coefficient of the InSAR monitoring results and the level monitoring results was 0.99, with a maximum absolute error of 17.3 mm/year, a minimum absolute error of 0.6 mm/year, and a standard deviation of 3 mm/year.Although the PS points fail to coincide with the corresponding levelling points, they still reflect the high accuracy of the InSAR results.


InSAR Validation

This work was validated in two steps.First, the second-class level deformation isoline of 2007 was selected and the results of SBAS-InSAR were preliminarily tested as a whole.Figure 4 shows that the settlement of the PS points that were obtained by SBAS-InSAR in 2007 were sufficiently consistent with that of the level measurements across the overall spatial d

tribution.

Furthermore, the data of 25
levelling benchmarks from 2003 to 2013 were selected for validation.The levelling benchmarks were taken as the original points, and all of the monitoring points within a radius of 150 m were extracted.The mean value of the extracted PS points was taken as the settlement estimation value and then verified according to the settlement of the levelling benchmarks.As shown in Figure 5, the correlation coefficient of the InSAR monitoring results and the level monitoring results from 2003 to 2010 was 0.95, the maximum absolute error was 17.9 mm/year, the minimum absolute error was 0.1 mm/year, and the standard deviation was 4 mm/year.From 2011 to 2013, the correlation coefficient of the InSAR monitoring results and the level monitoring results was 0.99, with a maximum absolute error of 17.3 mm/year, a minimum absolute error of 0.6 mm/year, and a standard deviation of 3 mm/year.Although the PS points fail to coincide with the corresponding levelling points, they still reflect the high accuracy of the InSAR results.


Acquisition of Mutation Information of the Beijing Plain

The spatial distribution characteristics of settlement information cannot be intuitively revealed, owing to the different densities of PS points in different areas and the large number of PS points.Fishnet analysis is the basic method of geospatial statistics; this method combines various data with natural geographic element data to express the spatial distribution characteristics for each unit.Therefore, we used the Fishnet tool to subdivide the research area into grids and study the spatial distribution characteristics of land subsidence on the Beijing Plain.The grid scale starts at 30 m (a regular grid), after which it exponentially increases up to 15,360 m.The minimum scale is 30 m, because the spatial resolution of the remote sensing image is 30 m.In 2017, Luo, et al. [41] divided the land subsidence on the Beijing Plain into 37 units based on the hydrogeological conditions.The maximum scale of this study was 15,360 m to maintain consistency with their data.We obtained the number of PS points at different scales via statistical methods to determine the optimal grid scale (Table 2).According to the information of the PS points from 2004 to 2015 that was acquired by InSAR technology, using the spatial join method, we subsequently calculated the corresponding PS points at different scales within the grid, and the settlement information of the corresponding grid in 2004-2015 was obtained.Finally, a grid of 960 m × 960 m was chosen, because it guaranteed that 88% of the grids had PS points from 2003-2015; moreover, the grid can be refined as m ch as possible to represent the research area.In total, 5575 grids were created.The Mann-Kendall mutation test was ultimately applied to each grid via Python, and the settlement mutation information in the Beijing Plain area was obtained (Figure 6).

In total, 4887 grids had settlement information from 2004 to 2015, of which 3792 grids passed the significance test (p = 0.05).Among them, 2744 grids had a single mutation and 1048 grids had multiple mutations.Among the single-mutation grids, 1344 grids only underwent mutation in 2015, 915 grids only underwent mutations in 2005, and 152 grids only underwent mutations in 2013; moreover, fewer single mutations occurred in other years (Table 3).In total, 4,887 grids had settlement information from 2004 to 2015, of which 3,792 grids passed the significance test (p=0.05).Among them, 2,744 grids had a single mutation and 1,048 grids had multiple mutations.Among the single-mutation grids, 1,344 grids only underwent mutation in 2015, 915 grids only underwent mutations in 2005, and 152 grids only underwent mutations in 2013; moreover, fewer single mutations occurred in other years (Table 3).Among the 1048 grids with multiple mutation years of the settlement rate, the grids in which mutations occurred in 2015 contained the most mutations; in 2015, there were 768 grids with mutations.There were a large number of grids (238 in total) in which mutations occurred both in 2013 and 2015.

In addition, there were 153 grid mutations from 2013 to 2015 and 63 grid mutations in both 2011 and 2015.For the grids in which mutations occurred in 2005, 127 grids underwent mutations.In addition, the number of grids in which multiple mutations occurred in other years was relatively small.


Discussion

According to the map of the overall trend of the settlement mutation rate (Figure 6), grids with single-year mutations were mostly distributed in the middle and lower parts of the Chaobai River alluvial-diluvial fan and the Yongding River alluvial-diluvial fan, and the grids with multiple-year mutations were mostly distributed at the top of the alluvial-diluvial fan.The groundwater level of the Quaternary system has dropping drastically may be a reason for this distribution, which results in a decrease in pore water pressure in the overburdened layer and an effective stress increase between particles.However, in the middle and upper parts of the alluvial fan, most of the Quaternary diving or shallow confined water has good permeability, and the groundwater level greatly fluctuates due to the influence of precipitation.Rainfall quantity and the rainfall intensity are uneven in different seasons.These differences may cause the groundwater level in the middle and upper parts of the alluvial fan to experience relatively large fluctuations; thus, the grids with mutations may be more variable.On the other hand, the Huairou emergency water source, the Machikou emergency water source, the Xishan emergency water source, and the Pinggu emergency water source are in the middle and upper parts of the alluvial fan.The groundwater level has greatly changed, being affected by the distribution of exploitation, resulting in greater fluctuations of land subsidence.Figure 6 shows that, at the edge of the Quaternary groundwater drop funnel, the grids with mutations in the land subsidence rate greatly vary.In the following sections, by considering the hydrogeological conditions, we will attempt to discuss the causes of grid changes in a single year and the grid changes in multiple years.


Analysis of Single-Year Mutation Grids

Among the 2744 land subsidence grids with a single mutation, 1344 grids occurred in 2015 (Figure 7).These grids are mainly distributed in the alluvial fans of the Chaobai River and Yongding River.The main reason is that the "South-to-North Water Transfer Middle Line" project was initiated beginning at the end of December 2014, and the water transfer year from 2014 to 2015 (the water transfer year refers to the annual water supply from November 1 to October 31) supplied 703 million cubic metres of water that was diverted to Beijing.Every day, eight water plants used more than 2.2 million cubic metres of the "southern water" (accounting for 70% of the total water supply in Beijing's urban areas).Moreover, the Beijing municipal government has carried out a large-scale replacement project for self-provided wells, and the overall amount of groundwater extraction has decreased by 250 million cubic metres.This phenomenon not only reduces the amount of groundwater exploitation, but it also changes the distribution of groundwater exploitation, which may be one of the main reasons for the sudden change in the grids in 2015.The Water Resources Bulletin pointed that the average precipitation in Beijing in 2015 was 583 millimetres, which was 33% more than the 439 millimetres of precipitation recorded in 2014.In terms of the temporal precipitation distribution, the cumulative precipitation during the flood season (i.e., from June to September) was 447 millimetres, which accounted for 77% of the annual precipitation, and this value was 8% less than the annual average of 488 millimetres.This result means that the excessive rainfall is relatively lower and the effective supply of precipitation is relatively greater.The data of the Beijing Geological Environment Monitoring Station showed that 150 million cubic metres of water were recharged to the Beijing Plain every year from 2015 to 2020.Groundwater systems are interlinked, although the depth of recharge is typically approximately 50 metres, which may also be one of the reasons for the sudden change in the grids in 2015.Finally, dynamic load is one of the causes of land subsidence in Beijing.Three subways were opened in 2015, according to the official Beijing subway website.In addition, Chen et al. reported that, because subway construction severely disturbs the strata, the deformation rate changes markedly during construction [33] and that the groundwater level will be artificially lowered during the subway construction period.However, six subways were newly built in Beijing during 2015.These events may have caused some grids to mutate during that year.

the effective supply of precipitation is relatively greater.The data of the Beijing Geological Environment Monitoring Station showed that 150 million cubic metres of water were recharged to the Beijing Plain every year from 2015 to 2020.Groundwater systems are interlinked, although the depth of recharge is typically approximately 50 metres, which may also be one of the reasons for the sudden change in the grids in 2015.Finally, dynamic load is one of the causes of land subsidence in Beijing.Three subways were opened in 2015, according to the official Beijing subway website.In addition, Chen et al. reported that, because subway construction severely disturbs the strata, the deformation rate changes markedly during construction [33] and that the groundwater level will be artificially lowered during the subway construction period.However, six subways were newly built in Beijing during 2015.These events may have caused some grids to mutate during that year.Figure 8 shows that single mutations occurred in 2005 in land subsidence, mostly in the upper and middle parts of the Yongding River alluvial fan, the Nankou alluvial fan, and the lower part of the Chaobai River alluvial fan.From the perspective of the district, most of these sites were located in south-eastern Haidian, south-eastern Changping, south-eastern Chaoyang, and north-eastern Tongzhou.With respect to the alluvial-diluvial fan of the Yongding River corresponding to the grid location, the reasons for the sudden change in 2005 may be as follows.First, the Xishan emergency water source area was built in 2005 in the area of groundwater exploitation.Quaternary groundwater is contained in the mining target layer of the emergency water source, which is located in the upper part of the alluvial fan of the Yongding River.This region is mainly composed of a single layer of sand and gravel with good permeability.A sudden increase in groundwater exploitation will lead to a sudden drop of the groundwater level, which may cause a sudden change in the corresponding grid land subsidence.In addition, according to the Water Resources Bulletin, the amount of precipitation of the Yongding River was the lowest, at only 386 mm.This result shows that the groundwater system of the Yongding River has little recharge, which may be one of the reasons for the sudden change in the land subsidence of the alluvial fan of the Yongding River in 2005.Figure 8 shows that single mutations occurred in 2005 in land subsidence, mostly in the upper and middle parts of the Yongding River alluvial fan, the Nankou alluvial fan, and the lower part of the Chaobai River alluvial fan.From the perspective of the district, most of these sites were located in south-eastern Haidian, south-eastern Changping, south-eastern Chaoyang, and north-eastern Tongzhou.With respect to the alluvial-diluvial fan of the Yongding River corresponding to the grid location, the reasons for the sudden change in 2005 may be as follows.First, the Xishan emergency water source area was built in 2005 in the area of groundwater exploitation.Quaternary groundwater is contained in the mining target layer of the emergency water source, which is located in the upper part of the alluvial fan of the Yongding River.This region is mainly composed of a single layer of sand and gravel with good permeability.A sudden increase in groundwater exploitation will lead to a sudden drop of the groundwater level, which may cause a sudden change in the corresponding grid land subsidence.In addition, according to the Water Resources Bulletin, the amount of precipitation of the Yongding River was the lowest, at only 386 mm.This result shows that the groundwater system of the Yongding River has little recharge, which may be one of the reasons for the sudden change in the land subsidence of the alluvial fan of the Yongding River in 2005.The reason for the sudden change in the land subsidence in the Nankou alluvial fan in 2005 is that, on the one hand, the Wenyu River groundwater system (where the Nankou alluvial fan is located) and the Chaobai River groundwater system belong to a subsystem of the same groundwater system, which has certain boundaries and hydraulic connections.However, the Huairou emergency water source, which is located at the top of the Chaobai River groundwater system, started operating

August 2003.After a period of operation, the water source will
destroy the balance of the groundwater system, and the nearby Wenyu River groundwater system will supply it, which may have led to the sudden change at the southern mouth diluvial fan area at the lower part in 2005.Second, in terms of groundwater recharge, from the perspective of precipitation, the depositional particles in the Nankou diluvial fan area were relatively larger and they had better permeability, and the groundwater level greatly fluctuated under the influence of precipitation.The average annual precipitation in Beijing was 468 millimetres in 2005, which was 13% less than that of the same p

iod in 2004.From the regional distribution of precipitation, th
precipitation in mountainous areas was greater than that in the plain areas, and the former mountainous areas in the north-western and south-western regions received increased amounts of precipitation.Part of the rainstorm centre was located in Nankou district.This result shows that, in 2005, the region not only supplied less water, but also supplied water less effectively.For a small portion of the grid in the lower part of the alluvialdiluvial fan of the Chaobai River, there may have been very little precipitation in that year.On the other hand, it is possible that water was released from the Huairou emergency water source in the upper part of the area (the water source was officially operated at the end of 2003) and from the Pinggu water source (the water source started operation in August 2004), and the exploitation of the emergency water sources will cause the rapid decline in nearby groundwater levels.Thus, the direction of the recharge-runoff-discharge of the groundwater system changes, which indirectly affects the groundwater level in the lower part of the alluvial fan.In summary, the reason for the mutation in land subsidence in 2005 may have been the increase in groundwater exploitation and the small amount of effective precipitation replenishment.The reason for the sudden change in the land subsidence in the Nankou alluvial fan in 2005 is that, on the one hand, the Wenyu River groundwater system (where the Nankou alluvial fan is located) and the Chaobai River groundwater system belong to a subsystem of the same groundwater system, which has certain boundaries and hydraulic connections.However, the Huairou emergency water source, which is located at the top of the Chaobai River groundwater system, started operating in August 2003.After a period of operation, the water source will destroy the balance of the groundwater system, and the nearby Wenyu River groundwater system will supply it, which may have led to the sudden change at the southern mouth diluvial fan area at the lower part in 2005.Second, in terms of groundwater recharge, from the perspective of precipitation, the depositional particles in the Nankou diluvial fan area were relatively larger and they had better permeability, and the groundwater level greatly fluctuated under the influence of precipitation.The average annual precipitation in Beijing was 468 millimetres in 2005, which was 13% less than that of the same period in 2004.From the regional distribution of precipitation, the precipitation in mountainous areas was greater than that in the plain areas, and the former mountainous areas in the north-western and south-western regions received increased amounts of precipitation.Part of the rainstorm centre was located in Nankou district.This result shows that, in 2005, the region not only supplied less water, but also supplied water less effectively.For a small portion of the grid in the lower part of the alluvial-diluvial fan of the Chaobai River, there may have been very little precipit tion in that year.On the other hand, it is possible that water was released from the Huairou emergency water source in the upper part of the area (the water source was officially operated at the end of 2003) and from the Pinggu water source (the water source started operation in August 2004), and the exploitation of the emergency water sources will cause the rapid decline in nearby groundwater levels.Thus, the direction of the recharge-runoff-discharge of the groundwater system changes, which indirectly affects the groundwater level in the lower part of the alluvial fan.In summary, the reason for the mutation in land subsidence in 2005 may have been the increase in groundwater exploitation and the small amount of effective precipitation replenishment.

The 152 grids where land subsidence changed abruptly in 2013 were mainly distributed in the Chaobai River alluvial fan (as shown in Figure 9).In the administrative area, they are mainly distributed in the Tongzhou district.The reason is that, in September 2012, the Beijing municipal government announced the selection of Tongzhou district as a new city deputy centre.After 2013, Tongzhou rapidly developed.Second, the Tongzhou area became dramatically "higher" and "larger".However, the increasing number of buildings and the growth in height constituted a single cause of Beijing's l nd subsidence.Third, after 2013, the density of both roads and subways in the Tongzhou area dramatically increased.In summary, the above three reasons may explain the sudden change in land subsidence in some areas of Tongzhou in 2013.The 152 grids where land subsidence changed abruptly in 2013 were mainly distributed in the Chaobai River alluvial fan (as shown in Figure 9).In the administrative area, they are mainly distributed in the Tongzhou district.The reason is that, in September 2012, the Beijing municipal government announced the selection of Tongzhou district as a new city deputy centre.Aft

2013, Tong
hou rapidly developed.Second, the Tongzhou area became dramatically "higher" and "larger".However, the increasing number of buildings and the growth in height constituted a single cause of Beijing's land subsidence.Third, after 2013, the density of both roads and subways in the Tongzhou area dramatically increased.In summary, the above three reasons may explain the sudden change in land subsidence in some areas of Tongzhou in 2013.


Analysis of Grids with Multiple Land Subsidence Mutation Years

The 1,048 grids with multiple mutation years were mainly distributed at the top of the alluvial fan, the boundary of the groundwater system, an the edge of the groundwater funnel.The largest number of grids containing mutations in 2015 was 768, followed by 152 grids that contained mutations in 2005.The grids that contained multiple mutations in years that included 2015 were mainly distributed at the top of the Chaobai River alluvial fan, near the Pinggu emergency water source, at the top of the Yongding River alluvial fan and at the bottom of the Chaobai River alluvial fan (as shown in Figure 10).


Analysis of Grids with Multiple Land Subsidence Mutation Years

The 1048 grids with multiple mutation years were mainly distributed at the top of the alluvial fan, the boundary of the groundwater system, and the edge of the groundwater funnel.The largest number of grids containing mutations in 2015 was 768, followed by 152 grids

at contained
mutations in 2005.The grids that contained multiple mutations in years that included 2015 were mainly distributed at the top of the Chaobai River alluvial fan, near the Pinggu emergency water source, at the top of the Yongding River alluvial fan and at the bottom of the Chaobai River alluvial fan (as shown in Figure 10).In terms of the administrative area, the sites were mainly distributed in southern Huairou, northern Shunyi, south-west Pinggu, and eastern Tongzhou.The reasons may be as follows.First, after the water diversion project of the south-to-north water transfer project at the end of 2014, the underground water was reduced in man places, such as the Huairou emergency source, Chaobai River greening water source well, and Pinggu emergency water source.The "southern water" was used to recharge groundwater in well-penetrated areas (such as in the Huairou emergency water source and Pinggu emergency water source area), and the relevant departments have completed the replacement of several self-p oviding wells, which not only reduces the amount of groundwater exploitation and optimizes the groundwater exploitation layout, but also effectively recharges and conserves the groundwater system.Second, the precipitation in this region was higher in 2015 than in the past 10 years, and the groundwater recharge was relatively sufficient.Third, the Beijing government has taken a series of measures to save water and strengthen the use of reclaimed water, which indirectly protects groundwater resources.In terms of the administrative area, the sites were mainly distributed in southern Huairou, northern Shunyi, south-west Pinggu, and eastern Tongzhou.The reasons may be as follows.First, after the water diversion project of the south-to-north water transfer project at the end of 2014, the underground water was reduced in many places, such as the Huairou emergency source, Chaobai River greening water source well, and Pinggu emergenc water source.The "southern water" was used to recharge groundwater in well-penetrated areas (such as in the Huairou emergency water source and Pinggu emergency water source area), and the relevant departments have completed the replacement of several self-providing wells, which not only reduces the amount of groundwater exploitation and optimizes the groundwater exploitation layout, but also effectively recharges and conserves the groundwater system.Second, the precipitation in this region was higher in 2015 than in the past 10 years, and the groundwater recharge was relatively sufficient.Third, the Beijing government has taken a series of measures to save water and strengthen the use of reclaimed water, which indirectly protects groundwater resources.

Among the grids ith multiple mutations (including one during 2015), the number with simultaneous mutations was the highest in 2013 and 2015, at 238 (as shown in Figure 11).The spatial distribution almost completely encompassed the upper and middle parts of the alluvial-diluvial fan of the Chaobai River.From the perspective of administrative regions, most of the sites were distributed in the south-western region of Huairou, the southern region of Miyun, the western region of Pinggu, and Tongzhou district.The reason for the sudden change in 2013 may be the long-term overexploitation of emergency water sources, which led to the imbalance of the groundwater system.Among the gr ds with multiple mutations (including one during 2015), the number with simultaneous mutations was the highest in 2013 and 2015, at 238 (as shown in Figure 11).The spatial distribution almost completely encompassed the upper and middle parts of the alluvial-diluvial fan of the Chaobai River.From the perspective of administrative regions, most of the sites were distributed in the south-western region of Huairou, the southern region of Miyun, the western region of Pinggu, and Tongzhou district.The reason for the sudden change in 2013 may be the long-term overexploitation of emergency water sources, which led to the imbalance of the groundwater system.

Subsequently, there were 153 grids with multiple mutations in 2013-2015, and these grids were mainly distributed in the Tongzhou area.The reason may be tha , after September 2012, Tongzhou was selected as the city deputy centre of Beijing.From 2013 to 2015, the Tongzhou area continuously expanded, the demand for water resources increased, and its dynamic and static load density increased.These reasons may explain the great change in land subsidence in Tongzhou during this period