Mapping land subsidence over the eastern Beijing city using satellite radar interferometry*

Beijing City has suffered from groundwater-induced subsidence since the late 1930s and the over-exploration of groundwater could lead to subsidence as much as ?12.0 cm?yr^?1. Previous studies on the ground deformation at Beijing City mainly focused on the period before the year of 2014 when a mega-engineering project was launched to reduce water shortage in Beijing. To study the most recent ground deformation, 19 L-band ALOS-1 PALSAR images (June 2007–January 2011), 24 C-band Sentinel-1 SAR images (June 2015–November 2016) together with 9 ALOS-2 PALSAR acquisitions (September 2014–February 2017) were analysed in this work. Levelling measurements were exploited to verify the ALOS-1-based time series InSAR (TS-InSAR) result while Sentinel-1 and ALOS-2 result were cross-verified with each other. Furthermore, the whole study area was divided into four sub-zones, and the result indicated that the subsidence rates over five townships, Cuigezhuan, Jinzhan, Liyuan, Songzhuang and Yanjiao were accelerating and more attentions should be paid. On the contrary, the town centre of Douge Zhuang township experienced a decreasing trend between these two temporal-periods. Additionally, the time series measurements with respect to five selected measurement points and the profile line along the subsidence hot spots were analysed.     

Land subsidence under different land use in the eastern Beijing plain,China 2005-2013 revealed by InSAR timeseries analysis

Land subsidence has been occurring in Beijing since the 1970s. Five major land subsidence areas have been formed: Dongbalizhuang–Dajiaoting, Laiguangying, Changping Shahe–Ba Xianzhuang, Daxing Yufa–Lixian, and Shunyi–Ping Gezhuang. In this paper, we studied on land subsidence in Dongbalizhuang–Dajiaoting and Laiguangying using small baseline subset interferometry and interferometric point target methods of 47 ENVISAT ASAR and 29 RADARSAT-2 data. The results showed that the degree of land subsidence in these areas varied significantly. The mean land subsidence rate ranged from 143.43 to 8.2 mm/a and from 132.11 to 7.3 mm/a during 2005–2010 and 2011–2013, respectively. We correlated the observed settlement with the land use (agricultural, residential, and industrial). Displacement in the agricultural areas was greater than that in the other areas from 2005 to 2013. Moreover, we compared the observed deformation and the groundwater level in phreatic and confined aquifers. There was a strong correlation between ground subsidence and the groundwater level and the ground settlement increased with a decrease in the groundwater level and the maximum correlation coefficient can reach 0.525. Furthermore, subsidence appeared to be associated with compressible deposits, suggesting that for 90–210-m thick compressible deposits, ground settlement is more likely to occur as the thickness of the compressible layer increases.     

Land subsidence prediction in Beijing based on PS-InSAR technique and improved Grey-Markov model

Land subsidence induced by excessive groundwater withdrawal has caused serious social, geological, and environmental problems in Beijing. Rapid increases in population and economic development have aggravated the situation. Monitoring and prediction of ground settlement is important to mitigate these hazards. In this study, we combined persistent-scatterer interferometric synthetic aperture radar with Grey system theory to monitor and predict land subsidence in the Beijing plain. Land subsidence during 2003–2014 was determined based on 39 ENVISAT advanced synthetic aperture radar (ASAR) images and 27 RadarSat-2 images. Results were consistent with global positioning system, leveling measurements at the point level and TerraSAR-X subsidence maps at the regional level. The average deformation rate in the line-of-sight was from ?124 to 7 mm/year. To predict future subsidence, the time-series deformation was used to build a prediction model based on an improved Grey-Markov model (IGMM), which adapted the conventional GM(1,1) model by utilizing rolling mechanism and integrating a k-means clustering method in Markov-chain state interval partitioning. Evaluation of the IGMM at both point level and regional scale showed good accuracy (root-mean-square error <3 mm; R² = 0.94 and 0.91). Finally, land subsidence in 2015–2016 was predicted, and the maximum cumulative deformation will reach 1717 mm by the end of 2016. The promising results indicate that this method can be used as an alternative to the conventional numerical and empirical models for short-term prediction when there is lack of detailed geological or hydraulic information.     

Land subsidence lagging quantification in the main exploration aquifer layers in Beijing plain,China

Land subsidence is rapidly developing across the Beijing Plain, China. Long-term intense overexploitation of groundwater is the main reason for land subsidence in Beijing. In this study, an optimized Small Baseline Subset (SBAS) interferometry method was developed to process 46 RADATSAT-2 images from 2011 to 2015 to investigate the spatial and temporal dynamics of land subsidence in the Beijing Plain. The lag time between land subsidence and groundwater exploitation was first analyzed by the Continuous Wavelet Transform (CWT) and Cross Wavelet Transform (XWT) methods Our study found that the maximum subsidence rate reached 141 mm per year. The analysis of the areas and volumes of the annual subsidence rates indicated that the overall deformation trend slowed down from 2011 to 2015. Our results indicate that the subsidence center is always located in the southeast of Chaoyang District from 2011 to 2015. The lag time between the observed subsidence and the groundwater level drops in the main exploration aquifer layers was 0.57–1.76 months. This information is helpful to reveal the mechanism of land subsidence and build hydrogeological model.     

InSAR视角下准噶尔盆地农业灌溉区地表形变特征和演化规律

抽取地下水进行农业灌溉是导致地下水位快速下降的重要因素，而长期过度开采地下水往往会引发地面沉降灾害，这种现象在干旱和半干旱地区非常普遍。为了研究农业灌溉超采引发的地表形变特征和演化规律，本文以准噶尔盆地南缘、天山北麓地带为研究区域，利用SBAS-InSAR技术对2003年—2009年覆盖呼图壁县的ENVISAT/ASAR升、降轨数据进行处理，获取了该地区的地表形变场，并结合研究区的农业灌溉方式、水资源补给和季节变化等资料对地面沉降的时空变化特征进行分析，为水资源和农业可持续发展提供参考意义。实验表明，研究区内主要有两个沉降幅度较大的漏斗，且都位于农田区域。2007年以前，研究区地表没有显著形变，之后发生了较大量级的沉降。采用传统灌溉方式和时针式灌溉系统的农业区平均沉降速率最高分别可达50 mm/a和30 mm/a，前者在时间上呈线性变化，而后者具有显著的周期性变化特征。在冬季时，采用时针式灌溉系统的地区地面抬升量可达40 mm，远大于传统灌溉方式的农田区域，而夏季地面沉降速率可达200 mm/a。对研究区农业灌溉活动进行分析后发现，农业灌溉造成的地下水超采是该地区地面沉降的主要影响因素，其形变机制与季节变化具有较高的相关性，在灌溉活动休止期内地表形变取决于地下水的补给量。研究区内的形变特征和影响因素分析将为地下水资源的充分利用和农业的可持续发展提供有效的信息。     

浅谈基于InSAR的南票煤矿开采沉陷区沉降量外业检核

经过时间序列InSAR技术处理后,得到监测时间段内南票区的地面形变平均速率、累计形变量、时间序列形变量,确定沉降区域。本文通过南票煤矿开采沉陷区基岩点、基本水准点、监测点布设以及这些点测量工作,提出山区采用水准测量和GPS测量联合方法,实现沉陷区沉降量核查。     

Land subsidence in major cities of Central Mexico: Interpreting InSAR-derived land subsidence mapping with hydrogeological data

Significant structural damages to urban infrastructures caused by compaction of over-exploited aquifers are an important problem in Central Mexico. While the case of Mexico City has been well-documented, insight into land subsidence problems in other cities of Central Mexico is still limited. Among the cities concerned, we present and discuss the cases of five of them, located within the Trans-Mexican Volcanic Belt (TMVB): Toluca, Celaya, Aguascalientes, Morelia, and Queretaro. Applying the SBAS-InSAR method to C-Band RADARSAT-2 data, five high resolution ground motion time-series were produced to monitor the spatio-temporal variations of displacements and fracturing from 2012 to 2014. The study presents recent changes of land subsidence rates along with concordant geological and water data. It aims to provide suggestions to mitigate future damages to infrastructure and to assist in groundwater resources management.Aguascalientes, Celaya, Morelia and Queretaro (respectively in order of decreasing subsidence rates) are typical cases of fault-limited land subsidence of Central Mexico. It occurs as a result of groundwater over-exploitation in lacustrine and alluvial deposits covering highly variable bedrock topography, typical of horst-graben geological settings. Aguascalientes and Toluca show high rates of land subsidence (up to 10 cm/yr), while Celaya and Morelia show lower rates (from 2 to 5 cm/yr). Comparing these results with previous studies, it is inferred that the spatial patterns of land subsidence have changed in the city of Toluca. This change appears to be mainly controlled by the spatial heterogeneity of compressible sediments since no noticeable change occurred in groundwater extraction and related drawdown rates. While land subsidence of up to 8 cm/yr has been reported in the Queretaro Valley before 2011, rates inferior to 1 cm/yr are measured in 2013–2014. The subsidence has been almost entirely mitigated by major changes in the water management practices of the city, i.e., a 122 km long pipeline bringing surface water from an adjacent state allowed to cease pumping in half of the wells.     

Subsidence in the Kathmandu Basin,before and after the 2015 Mw 7.8 Gorkha Earthquake,Nepal Revealed from Small Baseline Subset-DInSAR Analysis

Land subsidence in densely urbanized areas is a global problem that is primarily caused by excessive groundwater withdrawal. The Kathmandu Basin is one such area where subsidence due to groundwater depletion has been a major problem in recent years. Moreover, on 25 April 2015, this basin experienced large crustal movements caused by the Gorkha earthquake (Mw 7.8). Consequently, the effects of earthquake-induced deformation could affect the temporal and spatial nature of anthropogenic subsidence in the basin. However, this effect has not yet been fully studied. In this paper, we applied the SBAS-DInSAR technique to estimate the spatiotemporal displacement of land subsidence in the Kathmandu Basin before and after the Gorkha earthquake, using 16 ALOS-1 Phased Array L-band Synthetic Aperture Radar (PALSAR) images during the pre-seismic period and 26 Sentinel-1 A/B SAR images during the pre- and post-seismic periods. The results showed that the mean subsidence rate in the central part of the basin was about ?8.2 cm/year before the earthquake. The spatial extents of the subsiding areas were well-correlated with the spatial distributions of the compressible clay layers in the basin. We infer from time-series InSAR analysis that subsidence in the Kathmandu basin could be associated with fluvio-lacustrine (clay) deposits and local hydrogeological conditions. However, after the mainshock, the subsidence rate significantly increased to ?15 and ?12 cm/year during early post-seismic (108 days) and post-seismic (2015–2016) period, respectively. Based on a spatial analysis of the subsidence rate map, the entire basin uplifted during the co-seismic period has started to subside and become stable during the early-post-seismic period. This is because of the elastic rebound of co-seismic deformation. However, interestingly, the localized areas show increased subsidence rates during both the early-post- and post-seismic periods. Therefore, we believe that the large co-seismic deformation experienced in this basin might induce the local subsidence to increase in rate, caused by oscillations of the water table level in the clay layer.     

PALSAR和ASARPSI显著地表沉降探测与分析

永久散射体雷达干涉(PSI)已被广泛应用于城市缓慢沉降监测,但针对显著性地表沉降监测及PSI适用性的探讨较少。本文选取天津市与廊坊市城郊结合部的典型工业化城镇为研究区,分别以2007至2010年间的23幅L波段(23.6 cm)ALOS PALSAR影像和2007至2009年间的23幅C波段(5.6 cm)ENVISAT ASAR影像为数据源,使用基于相位空间相关性分析的PSI方法进行时空沉降提取,并对两个平台的结果进行对比分析和交叉验证。在此基础上,结合长短波SAR对沉降敏感性不同的特征探讨该PSI方法应用于显著沉降解算的可靠性,并结合地质条件及地下水开采信息对实验区沉降的时空分布进行深入分析和解释。研究表明,两平台的沉降结果十分吻合,在沉降漏斗区二者均方根误差(RMSE)为6.5 mm/a;实验区地表沉降呈现出显著的非均匀性,最大沉降速率超过200 mm/a,造成严重沉降的主要原因为地下水的大量开采;证实了PSI方法在显著地表沉降监测方面的有效性和可靠性。     

利用GRACE数据进行京津冀地区地面沉降与地下水耦合分析

地下水的过量开采已导致京津冀地区出现严重的地面沉降,为了分析地下水与地面沉降之间的耦合关系,首先利用GRACE-FO数据与GLDAS数据反演出京津冀地区2016—2019年地下水变化时序;然后利用MCTSB-InSAR技术反演出该地区同时段的沉降变化时序。通过试验分别获取地下水与地面沉降的差分变化序列和变化趋势线,并引入非弹性存储系数分析地下水对地面沉降影响力的变化规律。结果表明:①当地下水迅速升高或降低时,地面沉降速率减小或增大;当地下水持续升高时,地面沉降接近停止。②沉降越严重的区域,地下水与地面沉降变化趋势的相关性越强。③由地下水变化引起沉降的能力,随沉降等级升高而变大,且该能力在不同沉降等级中随时间变化的趋势也不同。     

Temporal and Spatial Evolution of Land Subsidence Induced by Groundwater Exploitation and Construction in the Eastern Chaoyang District,Beijing, China

Human activities have a large impact on land subsidence in great metropolitan areas. In this study, the RADARSAT-2 observation data over a 4-year period (from November 2010 to September 2014) is used to investigate the trends of land subsidence in the eastern Chaoyang District in Beijing, China, and to analyze the impact of human activities on it. The observation results indicate that the temporal and spatial evolutions trend of land subsidence in this area is uneven, with deformation rates ranging from ??116.52 to 7.1 mm/year. There are two large subsidence areas located in the northern part of the study area, and we find that the groundwater exploitation and deformation rates are strongly linearly correlated. The effect of construction on land subsidence has also been investigated. Over time, land subsidence is generally increased in the study area, with an exception at the Longfor Changying Galleria area, which is under construction during the observation period. Based on the time-series analysis of permanent scattered points in the 200 m buffer area, construction activities are observed to cause both land subsidence and upheaval deformation of the surrounding soil. During the construction period, the displacement of the surrounding soil is disturbed. After completion, the initial displacement of surrounding soil is tended to involve the uplift of the surrounding soil, followed by a gradual subsidence with time.     

InSAR to support sustainable urbanization over compacting aquifers: The case of Toluca Valley,Mexico

This paper illustrates how InSAR alone can be used to delineate potential ground fractures related to aquifer system compaction. An InSAR-derived ground fracturing map of the Toluca Valley, Mexico, is produced and validated through a field campaign. The results are of great interest to support sustainable urbanization and show that InSAR processing of open-access Synthetic Aperture Radar (SAR) data from the Sentinel-1 satellites can lead to reliable and cost-effective products directly usable by cities to help decision-making.The Toluca Valley Aquifer (TVA) sustains the water needs of two million inhabitants living within the valley, a growing industry, an intensively irrigated agricultural area, and 38% of the water needs of the megalopolis of Mexico City, located 40 km east of the valley. Ensuring water sustainability, infrastructure integrity, along with supporting the important economic and demographic growth of the region, is a major challenge for water managers and urban developers. This paper presents a long-term analysis of ground fracturing by interpreting 13 years of InSAR-derived ground displacement measurements. Small Baseline Subset (SBAS) and Persistent Scatterer Interferometry (PSI) techniques are applied over three SAR datasets totalling 93 acquisitions from Envisat, Radarsat-2, and Sentinel-1A satellites and covering the period from 2003 to 2016.From 2003 to 2016, groundwater level declines of up to 1.6 m/yr, land subsidence up to 77 mm/yr, and major infrastructure damages are observed. Groundwater level data show highly variable seasonal responses according to their connectivity to recharge areas. However, the trend of groundwater levels consistently range from −0.5 to −1.5 m/yr regardless of the well location and depth. By analysing the horizontal gradients of vertical land subsidence, we provide a potential ground fracture map to assist in future urban development planning in the Toluca Valley.     

Land subsidence characteristics of Jakarta between 1997 and 2005, as estimated using GPS surveys

Jakarta is the capital city of Indonesia with a population of about 12 million people, inhabiting an area of about 625 km². It is well known that several areas in Jakarta are subsiding rapidly. There are four different types of land subsidence that can be expected to occur in the Jakarta basin, namely: subsidence due to groundwater extraction, subsidence induced by the load of constructions (i.e., settlement of high compressibility soil), subsidence caused by natural consolidation of alluvial soil and tectonic subsidence. In addition to the leveling method, Global Positioning System (GPS) survey methods have been used to study land subsidence in Jakarta. In this paper, we characterize subsidence in the Jakarta basin using eight episodic/campaign GPS surveys between 1997 and 2005. The estimated subsidence rates are 1–10 cm/year. The observed subsidence rates in several locations show a positive correlation with known abstraction volumes of groundwater extraction. These basin-wide series of GPS measurements show how this type of measurement can play an important role in multiple public policy decision making in this rapidly growing area.     

郑徐高铁开封段地面沉降分层监测方法

高铁沿线施工期间引起地面沉降的因素错综复杂、相互影响,本文提出通过对沉降区域地下水位、孔隙水压力等进行分层观测,根据分层观测得到地下水位变化、孔隙水压力变化等数据,并采用回归分析法构建了沉降区域地面沉降与诸影响因子的关联模型,通过对郑徐客运专线开封段进行的地面沉降监测量与同时采用二等精密水准测量得到的监测点沉降量进行对比分析,结果表明本方法的有效性和可靠性。     

From ERS-1/2 to Sentinel-1: two decades of subsidence monitored through A-DInSAR techniques in the Ravenna area (Italy)

Land subsidence due to underground resources exploitation is a well-known problem that affects many cities in the world, especially the ones located along the coastal areas where the combined effect of subsidence and sea level rise increases the flooding risk. In this study, 25 years of land subsidence affecting the Municipality of Ravenna (Italy) are monitored using Advanced Differential Interferometric Synthetic Aperture Radar (A-DInSAR) techniques. In particular, the exploitation of the new Sentinel-1A SAR data allowed us to extend the monitoring period till 2016, giving a better understanding of the temporal evolution of the phenomenon in the area. Two statistical approaches are applied to fully exploit the informative potential of the A-DInSAR results in a fast and systematic way. Thanks to the applied analyses, we described the behavior of the subsidence during the monitored period along with the relationship between the occurrence of the displacement and its main driving factors.     

高速铁路沉降评估中工程沉降和区域沉降的可区分性研究

某高速铁路北京至济南段地处华北平原,沿线经过的大部分地区都存在着不同程度的地面沉降。对穿越不同地面沉降速率区段的高速铁路工程,沉降观测得到的桥墩沉降量是工程沉降与区域沉降的综合体现。由于地面沉降速率在空间上是不均匀的,本文对区分工程沉降和区域沉降进行了研究,并提出了对沉降观测的改进措施。     

利用水准数据的天津市地面沉降时空特征分析

地面不均匀沉降可能对城市的发展与人民的安全造成危害，天津市的地面沉降情况尤为严重。针对该问题，本文收集天津市2005—2012年、2016—2017年水准观测数据，以固定水准点位的沉降量、沉降速率、沉降加速率为状态向量，构建卡尔曼滤波模型，对天津市历年的水准观测数据进行滤波；根据滤波后的结果，本文利用多项式加权内插的方法，以距离、沉降速率、沉降加速率信息确定权值大小，对地面沉降情况进行内插；并以中误差作为精度评定参数，比较多种内插方法的精度。通过对内插结果的试验分析发现，2005—2017年天津市地面累计平均沉降量为394.477 5 mm，最大沉降量为1 143.5 mm；主要沉降区域为北辰、大港、塘沽等地区，且随着时间的增长，这些区域呈现漏斗式下沉。试验证明本文结合卡尔曼滤波与多项式加权内插的方式能够较好地反映地面沉降的时空特征分布情况并对未来一段时间的沉降情况进行预测，对天津市的城市发展及建设有一定的参考意义。     

Spatiotemporal changes of urban impervious surface area and land surface temperature in Beijing from 1990 to 2014

This study examined changes in urban expansion and land surface temperature in Beijing between 1990 and 2014 using multitemporal TM, ETM+, and OLI images, and evaluated the relationship between percent impervious surface area (%ISA) and relative mean annual surface temperature (RMAST). From 1990 to 2001, both internal land transformation and outward expansion were observed. In the central urban area, the high-density urban areas decreased by almost 7 km², while the moderate- and high-density urban land areas increased by 250 and 90 km², respectively, outside of the third ring road. From 2001 to 2014, high-density urban areas between the fifth and sixth ring roads experienced the greatest increase by more than 210 km², and RMAST generally increased with %ISA. During 1990–2001 and 2001–2014, RMAST increased by more than 1.5 K between the south third and fifth ring roads, and %ISA increased by more than 50% outside of the fifth ring road. These trends in urban expansion and RMAST over the last two decades in Beijing can provide useful information for urban planning decisions.     

基于SRTM数据的二轨D-InSAR监测城市地表沉降研究

合成孔径雷达差分干涉测量技术是新近发展起来的用于监测大范围地表形变的新技术,具有精度高、监测范围广等特点。本文首先深入分析了D-InSAR技术的基本原理,然后通过实验成功获取了上海市中心城区相隔7个月的地面沉降结果,通过与常规方法获取的沉降图的对比,发现二者具有高度的一致性,从而证明了利用D-InSAR技术监测城市地面沉降具有良好的应用前景。     

时序InSAR监测京雄城际铁路河北段地面沉降

本文以京雄城际铁路河北段固安站至雄安站沿线作为研究区,利用2018—2020年共34景Sentinel-1B影像,基于小基线集雷达干涉测量技术(SBAS-InSAR)获取京雄城际铁路河北段沿线的地面沉降时空分布信息,结合空间自相关分析方法,揭示研究区地面沉降的空间分布格局,并对沉降原因进行初步分析。研究结果表明,京雄城际铁路河北段沿线地面沉降发展由北向南存在一定的差异。北部年均沉降速率小于10 mm/a,南部沉降较为严重,最大年均沉降速率达-105.6 mm/a,且沿线西部年均沉降速率高于东部区域。通过分析影响因素得知,地面沉降量与地下水埋深值存在相关性,地下水埋深高的地区地面沉降量较高。同时结合研究区土地利用变化结果发现,城市化建设所产生的静载荷对京雄城际铁路沿线的地面沉降产生一定的影响。