Monitoring active open-pit mine stability in the Rhenish coalfields of Germany using a coherence-based SBAS method

With the recent progress in synthetic aperture radar (SAR) technology, especially the new generation of SAR satellites (Sentinel-1 and TerraSAR-X), our ability to assess slope stability in open-pit mines has significantly improved. The main objective of this work is to map ground displacement and slope instability over three open-pit mines, namely, Hambach, Garzweiler and Inden, in the Rhenish coalfields of Germany to provide long-term monitoring solutions for open-pit mining operations and their surroundings. Three SAR datasets, including Sentinel-1A data in ascending and descending orbits and TerraSAR-X data in a descending orbit, were processed by a modified small baseline subset (SBAS) algorithm, called coherence-based SBAS, to retrieve ground displacement related to the three open-pit mines and their surroundings. Despite the continuously changing topography over these active open-pit mines, the small perpendicular baselines of both Sentinel-1A and TerraSAR-X data were not affected by DEM errors and hence could yield accurate estimates of surface displacement. Significant land subsidence was observed over reclaimed areas, with rates exceeding 500 mm/yr, 380 mm/yr, and 310 mm/yr for the Hambach, Garzweiler and Inden mine, respectively. The compaction process of waste materials is the main contributor to land subsidence. Land uplift was found over the areas near the active working parts of the mines, which was probably due to excavation activities. Horizontal displacement retrieved from the combination of ascending and descending data was analysed, revealing an eastward movement with a maximum rate of ∼120 mm/yr on the western flank and a westward movement with a maximum rate of ∼ 60 mm/yr on the eastern flank of the pit. Former open-pit mines Fortuna-Garsdorf and Berghein in the eastern part of Rhenish coalfields, already reclaimed for agriculture, also show subsidence, at locations reaching 150 mm/yr. The interferometric results were compared, whenever possible, with groundwater information to analyse the possible reasons for ground deformation over the mines and their surroundings.     

Heterogeneous surface displacement pattern at the Hatchobaru geothermal field inferred from SAR interferometry time-series

We estimated surface displacements using persistent scatterer SAR interferometry (PS-InSAR) around the Hatchobaru geothermal field, Japan, from 18 ALOS/PALSAR images acquired from July 2007 to December 2010. Generally, geothermal fields, covered with natural targets such as rocky terrain and vegetation, have been one of the difficult targets for PS-InSAR analysis. However, we applied space adaptive filtering to increase the number of pixels for measuring surface displacement. The results of our analysis demonstrate ground subsidence with decaying velocity over the observation period around the geothermal field. The spatial pattern of ground subsidence includes sharp boundaries of subsidence that can be interpreted as fault traces. We demonstrated the usefulness of PS-InSAR analysis with the space adaptive filtering to estimate surface displacements with high spatial resolution and high spatial density around a geothermal field.     

时间序列InSAR技术辅助下的北京市高速公路网沉降监测应用

利用高分辨率TerraSAR-X影像,通过多主影像相干目标小基线干涉技术（MCTSB-InSAR）获取了北京市主城区2011-2015年的地表形变信息;通过建立缓冲区,提取并分析了北京市12条高速公路及一条快速公路的沉降现状;在高速公路沉降信息的基础上对机场第二高速进行了路面平整度分析。利用水准测量数据验证了监测结果,精度达4.6 mm/a,表明MCTSB-InSAR技术不仅可以快速监测出大区域高速公路交通网络的地表沉降,而且监测精度较高。     

高分辨率PSI上海市沉降探测及分析

本文在改进的永久散射体(PS)探测方法基础上,应用高分辨率永久散射体雷达干涉(PSI)探测上海市地表沉降,并对沉降原因进行了详细分析。实验选取2008年4月至2010年1月间,由德国卫星TerraSAR-X(TSX)所获取的18幅X波段(波长为3.1 cm)高分辨率SAR影像为数据源,进行PS探测、PSI建模、形变提取和分析。实验结果表明,改进的PS探测方法探测出的PS点是合理和可靠的,且高分辨率SAR影像对地面硬目标识别能力较强,显著提高了PS点的密度和覆盖范围。沉降探测结果显示,最大相对沉降速率达30 mm/yr,平均沉降速率为11.5 mm/yr。     

基于SBAS-InSAR的矿区地表沉降监测与分析

针对雷达差分干涉测量技术（differential interferometry synthetic aperture radar,DInSAR）易受时空失相关、大气相位延迟等因素的影响,采用小基线集技术（small baseline subset,SBAS）对13景TerraSAR-X数据进行时序处理,估算并去除了残余数字高程模型误差、大气延迟误差和轨道误差,分析了研究区域2012-2013年沉降速率,发现2310和1301两工作面最大下沉速率分别为40 mm/a和50 mm/a;分析时序累计沉降值发现,2306、2308、2310这3个工作面在2012年11月15日之前地面沉降并不明显;在2310和1301工作面分别提取3个失相干现象表现较为缓慢的候选点进行时序分析,发现沉降值和时间成线性变化关系,且开采时间越早,其沉降特性越符合线性变化;将SBAS和DInSAR两种方法获得的累计沉降值进行差值分析,发现两种方法的差值在5 mm以内;在2310工作面走向方向和倾向方向选取若干观测点,并提取各观测点的时序沉降值对工作面时序沉降进行量化分析,实验表明SBAS-InSAR技术在矿区地表沉降监测与分析方面具有良好的应用前景。     

改进的PS探测方法及其应用

提出一种改进的永久散射体PS(Persistent Scatterer)探测方法,首先基于时序SAR影像的振幅信息,利用振幅和振幅离差指数双阈值法获取PS候选点,然后根据干涉相位空间相关性,对所有候选点进行相位稳定性分析,筛选出既满足雷达反射强度高、又符合散射稳定性要求的真实PS点,最后基于所选PS点提取形变信息。选取上海市局部区域作为实验区,以16幅TerraSAR-X(TSX)影像为数据源,进行PS探测和形变信息提取。结果表明,相比于已有的PS探测方法,本文方法能够更准确、可靠地识别出时序SAR影像中的PS点。形变结果显示,研究区域内最大沉降速率为-51mm/a,平均沉降速率为-15.8mm/a,以水准实测数据进行检核,验证了基于改进的PS探测方法进行地表形变信息提取结果的可靠性。     

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.     

PS-InSAR和SBAS-InSAR技术对昆明主城区地面沉降监测的对比分析

以昆明主城区为例,分别利用PS-InSAR和SBAS-InSAR技术对2014—2017年间29景升轨Sentinel-1A数据进行沉降监测,对比两种技术得到的沉降结果,进行剖面图分析与时序分析。结果表明,PS-InSAR和SBAS-InSAR技术监测结果具有一致性、相关性和可靠性。研究发现,昆明市沉降漏斗主要位于居民区、地铁、道路、高速公路以及滇池等区域,最大年沉降速率可达-39.580mm/a,累积沉降量达到85mm。研究表明,昆明主城区地面沉降主要由于近几年城市和轨道交通建设的飞跃发展,导致居民区和交通网络密集,地面载荷增加,地下隧道开挖与地下水开采等原因引起地面软土地层下沉。     

利用SBAS时序技术识别地热资源开采引起的地表形变

本文以地热资源丰富的雄安新区为研究对象,利用SBAS技术对2015年7月至2019年2月的42幅Sentinel-1区域卫星影像和2016年8月至2018年12月的10幅ALOS-2区域卫星影像进行处理,以此获取雄安新区的时间序列地表形变。研究结果表明,雄安新区整体地表沉降严重,且沉降区主要集中在雄县和容城县,其中雄县沉降最为严重,沉降中心累计沉降量超过330 mm,容城县沉降中心累计下沉量超过120 mm。结合地热井位置和地热开采历史,由时间序列分析表明,大量开采地热资源会引起地表沉降现象。     

不同手段和尺度下地面沉降时空变化特征分析——以深圳大空港区为例

本文利用不同手段和不同尺度的地面监测方法,分析城市填海区地面沉降的时空变化特征。通过深圳市大空港区的地面沉降监测工程实例分析,结合精密水准测量和遥感监测手段,实现了不同尺度地面沉降监测,获取该区域地面沉降时空变化特征。研究结果表明,精密水准测量方法与遥感监测方法均发现了会展中心周围道路产生的沉降,具有相同的变化特征。精密水准监测方法能较精确地获取局部地面沉降量,在石围路的建筑基坑周边产生最大为-303.6mm的地面沉降。结合不同手段和不同尺度的监测方法,可以较全面获取监测对象的形变特征。     

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.     

Monitoring dam structural health from space: Insights from novel InSAR techniques and multi-parametric modeling applied to the Pertusillo dam Basilicata,Italy

The availability of new constellations of synthetic aperture radar (SAR) sensors is leading to important advances in infrastructure monitoring. These constellations offer the advantage of reduced revisit times, providing low-latency data that enable analysis that can identify infrastructure instability and dynamic deformation processes.In this paper we use COSMO-SkyMed (CSK) and TerraSAR-X (TSX) data to monitor seasonal induced deformation at the Pertusillo dam (Basilicata, Italy) using multi-temporal SAR data analysis. We analyzed 198 images spanning 2010–2015 using a coherent and incoherent PS approach to merge COSMO-SkyMed adjacent tracks and TerraSAR-X acquisitions, respectively. We used hydrostatic-seasonal-temporal (HST) and hydrostatic-temperature-temporal (HTT) models to interpret the non-linear deformation at the dam wall using ground measurements together with SAR time-series analysis. Different look geometries allowed us to characterize the horizontal deformation field typically observed at dams. Within the limits of our models and the SAR acquisition sampling we found that most of the deformation at the Pertusillo dam can be explained by taking into account only thermal seasonal dilation and hydrostatic pressure. The different models show slightly different results when interpreting the aging term at the dam wall. The results highlight how short-revisit SAR satellites in combination with models widely used in the literature for interpreting pendulum and GPS data can be used for supporting structural health monitoring and provide valuable information to ground users directly involved in field measurements.     

黄土填方区地表沉降时序InSAR监测分析

“削山填沟造地”等岩土工程在湿陷性黄土沟壑地区屡见不鲜,掌握填方区沉降情况具有重要意义。本文收集了2017年11月—2020年12月获取的56景TerraSAR-X StripMap模式影像,利用时序InSAR技术监测了陕北某湿陷性黄土填方地基工程的沉降信息,并与2017年11月—2020年12月期间监测区3个水准点的沉降测量结果比对。结果表明,在填方区地表以沉降为主,在挖方区地表以抬升为主,研究区存在有1处较为明显的地表沉降情况,位于填挖边界线附近填方区内,形变速率范围为-40~-20 mm/a,最大形变速率达-49.9 mm/a,累计量为-151.6 mm,时序InSAR形变结果和实地水准结果吻合性较好,垂直方向形变速率中误差为1.8 mm/a,表明时序InSAR技术在湿陷性黄土填挖方区变形监测中具有较好的应用价值。     

Identifying natural and anthropogenically-induced geohazards from satellite ground motion and geospatial data: Stoke-on-Trent,UK

Determining the location and nature of hazardous ground motion resulting from natural and anthropogenic processes such as landslides, tectonic movement and mining is essential for hazard mitigation and sustainable resource use. Ground motion estimates from satellite ERS-1/2 persistent scatterer interferometry (PSI) were combined with geospatial data to identify areas of observed geohazards in Stoke-on-Trent, UK. This investigation was performed within the framework of the EC FP7-SPACE PanGeo project which aimed to provide free and open access to geohazard information for 52 urban areas across Europe. Geohazards identified within the city of Stoke-on-Trent and neighbouring rural areas are presented here alongside an examination of the PanGeo methodology.A total of 14 areas experiencing ground instability caused by natural and anthropogenic processes have been defined, covering 122.35 km². These are attributed to a range of geohazards, including landslides, ground dissolution, made ground and mining activities. The dominant geohazard (by area) is ground movement caused by post-mining groundwater recharge and mining-related subsidence (93.19% of total geohazard area), followed by landsliding (5.81%). Observed ground motions along the satellite line-of-sight reach maxima of +35.23 mm/yr and −22.57 mm/yr. A combination of uplift, subsidence and downslope movement is displayed.‘Construction sites’ and ‘continuous urban fabric’ (European Urban Atlas land use types) form the land uses most affected (by area) by ground motion and ‘discontinuous very low density urban fabric’ the least. Areas of ‘continuous urban fabric’ also show the highest average velocity towards the satellite (5.08 mm/yr) and the highest PS densities (1262.92 points/km²) along with one of the lowest standard deviations. Rural land uses tend to result in lower PS densities and higher standard deviations, a consequence of fewer suitable reflectors in these regions. PSI is also limited in its ability to identify especially rapid ground motion. As a consequence the supporting geospatial data proved especially useful for the identification of landslides and some areas of ground dissolution. The mapped areas of instability are also compared with modelled potential geohazards (the BGS GeoSure dataset).     

基于PALSAR数据的青藏高原冻土形变检测方法研究

季节性冻胀和融沉导致的地面形变是青藏高原冻土区建设施工与维护的主要问题。对冻融造成的形变进行有效监测是青藏铁 路建设与维护的前提。差分干涉测量技术是地表形变监测的重要手段之一,PALSAR（L波段的合成孔径雷达）数据在非城市区域具 有较高的相关性,适合青藏高原冻土区的地表形变监测。本文选用4景覆盖研究区域的PALSAR数据,研究利用该数据进行冻土形变 检测的方法,并对其检测结果进行了分析。结果表明,该方法与水准测量方法有较好的一致性。     

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.     

GPS relative positioning at a precision level of one part per billion

Six days of data from the GIG'91 experiment have been analyzed with a fiducial strategy. The results obtained with our orbital software GEOSAT, show daily horizontal and length repeatabilities at the level of 1 part in 10⁹ plus 2 mm for baseline lengths up to 4000 km. A direct comparison with results from the GIPSY software shows, with some exceptions, mean differences at the sub-cm level. After transformation to ITRF'90 the rms of the coordinate differences is 14.8 mm. Studies of orbital predictions and comparisons with external high precision orbits indicate a mean orbit precision and accuracy of around 35 cm in each cartesian coordinate. Correlations between the GEOSAT and GIPSY solutions indicate some common model deficiencies.     

Lithology-controlled subsidence and seasonal aquifer response in the Bandung basin,Indonesia, observed by synthetic aperture radar interferometry

Land subsidence in the Bandung basin, West Java, Indonesia, is characterized based on differential interferometric synthetic aperture radar (DInSAR) and interferometric point target analysis (IPTA). We generated interferograms from 21 ascending SAR images over the period 1 January 2007 to 3 March 2011. The estimated subsidence history shows that subsidence continuously increased reaching a cumulative 45 cm during this period, and the linear subsidence rate reached ∼12 cm/yr. This significant subsidence occurred in the industrial and densely populated residential regions of the Bandung basin where large amounts of groundwater are consumed. However, in several areas the subsidence patterns do not correlate with the distribution of groundwater production wells and mapped aquifer degradation. We conclude that groundwater production controls subsidence, but lithology is a counteracting factor for subsidence in the Bandung basin. Moreover, seasonal trends of nonlinear surface deformations are highly related with the variation of rainfall. They indicate that there is elastic expansion (rebound) of aquifer system response to seasonal-natural recharge during rainy season.     

武汉市地面沉降时空分异特征及地理探测机制

城市地区地面沉降造成地面高程损失,威胁各类设施的安全运行,影响地表径流和水文循环,监测地面沉降现状并揭示其形成机制,对于城市可持续发展具有重要意义。以2007—2011年的ALOS-PALSAR影像和2015—2019年的Radarsat-2影像为数据源,基于小基线集技术（small baseline subset interferometric synthetic aperture radar, SBAS-InSAR）获取武汉市两个监测阶段的地面沉降速率、沉降时间序列,并利用地理探测器揭示规划单元尺度地面沉降的主导驱动因子及驱动因子之间的交互作用机制。结果表明:(1)2007—2011?年和2015—2019年地面沉降平均速率分别为-3.53 mm/a和-1.48 mm/a。地面沉降较为显著的区域:2007—2011年,是汉口、沙湖沿岸及以北、南湖以西和白沙洲地区;2015—2019年,是汉口、沙湖北和白沙洲地区。(2)?局部性、阶段性、与自然条件及人类活动相关是两个时期武汉市地面沉降演变的3个特点。(3)水文地质条件作为必要条件,通过与地面荷载、地下空间开发、工程施工因素交互作用形成武汉市地面沉降时空格局,2007?—2011年的工程施工因素、2015—2019年的地面荷载因素与水文地质条件交互作用明显较强。     

Swarm系列卫星非差运动学厘米级精密定轨

采用2015年5月24日—30日的Swarm星载GPS双频观测数据,基于Melbourne-Wübbena(MW)和消电离层线性组合,在精密单点定位技术的基础上,采用批处理最小二乘估计法对不同轨道高度的Swarm系列卫星进行非差运动学精密定轨。利用星载GPS相位观测值残差、与欧空局发布的简化动力学轨道对比,以及SLR检核3种方法对Swarm系列卫星非差运动学定轨结果进行精度评估。结果表明:①Swarm系列卫星星载GPS相位观测值残差RMS为6~7 mm;②与欧空局发布的简化动力学轨道进行求差,径向、切向及法向轨道差值RMS为2~4 cm;③与欧空局发布的运动学轨道进行求差,径向、切向及法向轨道差值RMS为1~2 cm;④SLR检核结果表明Swarm-A/B/C卫星轨道精度为3~4 cm。因此,采用非差运动学定轨方法与本文提供的定轨策略进行Swarm系列卫星精密定轨是切实可行的,定轨精度为厘米级。