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

Monitoring land surface deformation using persistent scatterers interferometric synthetic aperture radar technique

Land subsidence is one of the major hazards occurring globally due to several reasons including natural and human activities. The effect of land subsidence depends on the extent and severity. The consequences of this hazard can be seen in many forms including damaged of infrastructures and loss of human lives. Although land subsidence is a global problem, but it is very common in urban and sub urban areas especially in rapidly developing countries. This problem needs to be monitored effectively. Several techniques such as land surveying, aerial photogrammetry and Global Positioning System (GPS) can be used to monitor or detect the subsidence effectively but these techniques are mostly expensive and time consuming especially for large area. In recent decades, Interferometric Synthetic Aperture Radar (InSAR) technique has been used widely for the monitoring of land subsidence successfully although this technique has several limitations due to temporal decorrelation, atmospheric effects and so on. However, the uncertainties related to InSAR technique have been reduced significantly with the recent Persistent Scatterers Interferometric Synthetic Aperture Radar (PSInSAR) technique which utilized a stack of interferograms generated from several radar images to estimate deformation by finding a bunch of stable points. This study investigates the surface deformation focusing on Kuala Lumpur, a rapidly growing city and Selangor using PSInSAR technique with a set of ALOS PALSAR images from 2007 to 2011. The research methodology consists of several steps of image processing that incudes i) generation of Differential Interferometric Synthetic Aperture Radar (DInSAR), ii) selection of Persistent Scatterers (PS) points, iii) removal of noise, iv) optimization of PS point selection, and v) generation of time series deformation map. However, special consideration was given to optimize the PS selection process using two master images. Results indicate a complete variation of mean line-of-sight (LOS) velocities over the study area. Stable areas (mean LOS=1.1 mm/year) were mostly found in the urban center of Kuala Lumpur, while medium rate of LOS (from 20 mm/year to 30 mm/year) was observed in the south west area in Kuala Langat and Sepang districts. The infrastructures in Kuala Lumpur are mostly stable except in Kuala Lumpur International Airport (KLIA) where a significant subsidence was detected (28.7 mm/year). Meanwhile, other parts of the study area such as Hulu Langat, Petaling Jaya and Klang districts show a very low and non-continuous movement (LOS < 20 mm/year), although comparatively higher subsidence rate (28 mm/year) was detected in the mining area. As conclusion, PSInSAR technique has a potential to monitor subsidence in urban and sub urban areas, but optimization of PS selection processing is necessary in order to reduce the noise and get better estimation accuracy

Evaluating Multi-Temporal DInSAR Measurements of Ground Surface Deformation Around the Rhenish Coalfields in Germany Using Sentinel-1 SAR Imagery

Ground surface deformation caused by land subsidence is a critical aspect when assessing safety and environmental impacts of open-pit mining operations. These mining-induced deformations can cause severe damage to local infrastructure and buildings located in the region. Monitoring mining operations using traditional techniques is laborious, costly and time consuming, with many locations being difficult to access. As a result, in-situ observations of surface movement and structural stability are often sparsely completed along areas of high risk or concern. Remote sensing observations facilitate a reliable, temporally and spatially continuous monitoring process of mining operations, regardless of physical accessibility, at large scales with minimal cost. In this study, two commonly used C-band Differential Interferometric Synthetic Aperture Radar (DInSAR) methods, namely Small Baseline Subset (SBAS) and Persistent Scatterer Interferometry (PSI), are applied to Sentinel-1 data to show consistency when observing land displacement in the Rhenish coalfields region in Germany. Deformation rates around the Rhenish coalfields are assessed for 2015 and 2016, independently, using both techniques and compared using GNSS station data of three nearby stations part of the EUREF network. The primary focus on this study is to determine the performance differences between the two techniques. The PSI deformation estimates showed significantly lower variation in measurements but had a severely limited spatial coverage compared to the SBAS estimates. Comparison of annual deformation rates measured by the two techniques showed a stronger agreement for 2016 (RMSE = 8.7 mm/year) than 2015 (RMSE = 26 mm/year). The discrepancy between annual comparisons is caused by overestimation of surface subsidence rates by the 2015 SBAS measurements. This overestimation is likely the result of atmospheric artefacts propagating into the phase unwrapping due to a lack of sufficient SAR scenes available in 2015. It was determined that the SBAS technique benefits more heavily from additional data and when data is limited, the PSI technique proves to be more reliable for measuring annual deformation rates. Despite this, both techniques observed similar overall deformation patterns across the landscape but reported different magnitude of deformation; commonly observed in other studies. The two methods showed stronger agreement over targets characterized as hard impermeable surfaces (e.g. built-up areas). Although only one of the three publicly available GNSS sites coincided with the InSAR study domain, when compared with base station GNSS sites, the DInSAR estimates are within the uncertainty range of the GNSS measurements, indicating good potential. This thesis illustrates the applicability of using satellite remote sensing observations to monitor mining-induced surface deformations to supplement traditional monitoring techniques

Long-Term Continuously Updated Deformation Time Series From Multisensor InSAR in Xi'an, China From 2007 to 2021

According to historical records, land subsidence has been occurring in Xi'an, China, since the 1960s, characterized by complex land subsidence patterns. This subsidence has the potential to cause serious societal and economic problems during the process of urbanization. Long-term, large-scale monitoring and dynamic high precision tracking of the evolution of surface deformation associated with geohazards is a prerequisite for effective prevention or advance warning of geological disasters. Synthetic aperture radar interferometry (InSAR), a satellite remote sensing technology, can facilitate such monitoring. Both currently operating and planned SAR satellites would provide extensive amounts of SAR data. In this article, we describe an approach for continuously updating long-term multisensor InSAR deformation time series using robust sequential least squares. It is successfully applied to near-real-time monitoring of the long-term evolution of surface deformation in Xi'an, China, from January 3, 2007 to February 12, 2021, using four SAR satellites: ALOS/PALSAR-1, TerraSAR-X, ALOS/PALSAR-2, and Sentinel-1A. In order to analyze deformation evolution, temporal independent component analysis was used to interpret deformation patterns. We found that land subsidence in Xi'an has slowed and even halted in some areas. However, large areas are uplifting, which presents a potential for geohazards. We conclude that the proposed approach using continuously updated deformation time series from multisensor InSAR can provide near-real-time deformation measurements, which are necessary for an early warning system

Karst landslides detection and monitoring with multiple SAR data and multi-dimensional SBAS technique in Shuicheng, Guizhou, China

Shuicheng District is a karst mountain area, located in Guizhou Province, China. Its fragile stratum and frequent underground mining activities makes it prone to landslides. Owning to its wide coverage and frequent revisit, the InSAR technology has advantages in potential landslide identification and deformation monitor. However, affected by dense vegetation and atmospheric delay, it is much difficult to get sufficient effective targets to derive the deformation in this area. Besides, deformation derived from single orbit SAR data can result in the missing identification of some potential landslides and the misinterpreting of the real kinematics process of landslides. In this study, the multi-source SAR data, atmospheric error correction by quadratic tree image segmentation method, and phase-stacking method were selected to derive the surface deformation of this area. Besides, DS-InSAR and MSBAS method were combined to derive the deformation of Pingdi landslide. First, the potential landslides in this area were identified, surface deformation result, optical remote sensing images and geomorphological features were jointly considered. Then, the landslide distribution characteristics was analyzed in terms of slope, elevation and stratum. After that, the deformation along the LOS direction was acquired using the DS-InSAR method. The MSBAS method was used to retrieve the two-dimensional deformation of Pingdi landslide. Finally, the comprehensive analysis of triggering factors and failure process were conducted according to the spatial-temporal deformation characteristics and field investigation. The results indicated that landslides in Shuicheng district were mostly located in the junction of T1 and P3 stratum and mining related. Mining activity was the main cause of the Pingdi landslide deformation, the precipitation was the driving factor of the landslide instability. The research provides an insight into the explore the unstable slope distribution characteristic and the failure process of the landslides

Displacement monitoring in airport runways by persistent scatterers SAR interferometry

Deformations monitoring in airport runways and the surrounding areas is crucial, especially in case of low-bearing capacity subgrades, such as the clayey subgrade soils. An effective monitoring of the infrastructure asset allows to secure the highest necessary standards in terms of the operational and safety requirements. Amongst the emerging remote sensing techniques for transport infrastructures monitoring, the Persistent Scatterers Interferometry (PSI) technique has proven effective for the evaluation of the ground deformations. However, its use for certain demanding applications, such a as the assessment of millimetric differential deformations in airport runways, is still considered as an open issue for future developments. In this study, a time-series analysis of COSMO-SkyMed satellite images acquired from January 2015 to April 2019 is carried out by employing the PSI technique. The aim is to retrieve the mean deformation velocity and time series of the surface deformations occurring in airport runways. The technique is applied to Runway 3 at the “Leonardo da Vinci” International Airport in Rome, Italy. The proposed PSI technique is then validated by way of comparison with the deformation outcomes obtained on the runway by traditional topographic levelling over the same time span. The results of this study clearly demonstrate the efficiency and the accuracy of the applied PSI technique for the assessment of deformations in airport runways

An investigation of ongoing displacements of active faults in the Gobi desert using persistent scatterer interferometric synthetic aperture radar technique to support the permanent disposal of high-level waste in Beishan, China

This research demonstrated the application of PSInSAR method in identifying and characterising the micro-displacements along active faults in Beishan to support the selection of GDF host rock. This research first distinguishes and separates the tectonic induced and non-tectonic induced deformation within three study areas at Suanjingzi, Jiujing and Xinchang. Through the application of coherence change detection, it found the granite outcrop areas characterised by high coherence provide more robust results of tectonic activity. The Quaternary sediments covered areas which are characterised by low coherence usually show higher deformation rates due to the impacts of erosion and deposition. The tectonic induced displacements generally range from -0.4 to 0.4 mma-1 and are dominated by fault bound tectonic movements. As a part of wrench faut zone, Beishan is impacted by a NE-SW trended maximum in situ compressive stress field (σ1). To correlate the visible valleys, gullies, or cracks in Google Earth imagery with the SAR image deformation discontinuities, this study mapped and characterised more than 40 active faults in the three study areas, these include (1) the NE-SW trended sinistral strike-slip faults triggered by extension and (2) the NW-SE/W-E trended reverse faults triggered by maximum compression. The fault activity is characterised by subtle (minor) displacement rate value difference between the two sides of the fault plane. This research successfully improved the understanding of local structural geology and provided moderate guidance for the selection of HLW disposal sites in China. It was indicated that Xinchang has the highest tectonic stability, and this is then followed by Jiujing and Suanjingzi. This kind of displacement rate difference is possible due to the angle difference towards the Sanweishan Fault Zone. To trace and characterise the undiscovered active fault planes, the PSInSAR approach also benefits the prediction of earthquake by improving the positioning of the potential epicentres.Open Acces

Advanced exploitation of Sentinel-1 data for supporting landslide risk analysis

Tesi en modalitat de compendi de publicacionsSatellite Synthetic Aperture Radar Interferometry (InSAR) and Persistent Scatterer Interferometry (PSI) are now consolidated tools for ground movement detection and monitoring. Sentinel-1 (S1) is the first satellite providing free data access and ensuring a regular acquisition worldwide, every 6 days, increasing its potential for long-term monitoring applications. Several regional and national ground motion services are already active, providing products based on S1 data. Soon in 2022 the first European Ground Motion Service (EGMS) will be available and freely provide a displacement map over the whole Europe, with annual updates. This implies a strong expansion of availability of PSI-based displacement maps and an easy access for anyone, with an increasing interest among a wider range of users, including public or governmental institutions, academia, industry, and citizens. The analysis and interpretation of this amount of data is difficult and time consuming, mostly for non-expert InSAR users. The objective of this work is developing methodologies to simplify the operational use of PSI displacement maps, generating derived products with a clear message, easy-to-interpret, and fast to read. We propose a method to be applied over regional scale PSI displacement maps, to fast detect the most significant Active Deformation Areas (ADAs). The ADA map is a first product that allows a fast focusing on the active areas, to prioritize further analysis and investigation. Starting from the ADAs, the potential phenomena are attributed to each area through a preliminary interpretation based on auxiliary data, to derive the Geohazard Activity Map. In this work, a methodology to include the ADA information in the Civil Protection Activities is proposed, with the main output called Vulnerable Elements Activity Maps (VEAM). An application of the VEAM is illustrated in the Canary Islands. Furthermore, the ADA map is used in the Valle d'Aosta Region (Northern of Italy) to generate vulnerability and potential loss maps. Finally, a methodology to derive potential damage maps of the exposed buildings, based on the spatial gradients of movement, is proposed, and applied in a coastal area of the Province of Granada (Spain). A pack of software tools has been developed based on the proposed methods to extract ADA and then classify them to generate a Geohazard Activity Map. The set of tools is called ADATools, it is open-access, easy to use and fast, improving the operational exploitation of PSI regional-scale displacement maps. All the methodologies have been developed in the frame of several European projects (Safety, U-Geohaz, MOMIT and RISKCOAST), and are aimed at supporting the multi-scale territorial management and risk analysis activities, with a specific focus on landslides.La interferometría satelital radar (InSAR) y la interferometría de dispersores persistentes (PSI) son herramientas consolidadas para la detección y el monitoreo de movimientos de la superficie de la Tierra. Sentinel-1 (S1) es el primer satélite que proporciona acceso gratuito a los datos y garantiza una adquisición regular en todo el mundo, cada 6 días, aumentando su potencial para aplicaciones de monitoreo a largo plazo. Varios Ground Motion Services regionales y nacionales ya están activos, proporcionando productos basados en datos S1. Pronto, en 2022, el primer servicio europeo (European Ground Motion Service - EGMS) estará disponible y facilitará libremente un mapa de movimientos de toda Europa, con actualizaciones anuales. Esto implica un aumento de la disponibilidad de mapas de movimientos basados en PSI y un fácil acceso para cualquier persona, con un interés creciente entre una amplia gama de usuarios, incluyendo instituciones públicas o gubernamentales, academias, industrias y ciudadanos. El análisis e interpretación de esta cantidad de datos es difícil y consume mucho tiempo, mayormente para usuarios no expertos en la técnica. El objetivo de este trabajo es desarrollar metodologías para simplificar el uso operativo de los mapas de desplazamiento PSI, generando productos derivados con un mensaje claro, fácil de interpretar, y rápido de leer. Se propone un método para detectar rápidamente las Áreas de Deformación Activas (ADAs) más significativas, a partir de mapas de desplazamiento PSI de escala regional. El mapa de las ADAs es un primer producto que permite un enfoque rápido en las áreas activas, útil para priorizar el análisis y las investigaciones adicionales. A partir de las ADAs, se propone una interpretación preliminar basada en datos auxiliares, que atribuye a cada área el fenómeno que está detrás del movimiento, generando el Geohazard Activity Map (GAM). Después, se propone una metodología para incluir la información de las ADAs en las actividades de protección civil, generando los Vulnerable Element Activity Maps (VEAM), a través de su aplicación en las Islas Canarias. Además, el mapa de las ADAs se utiliza en la región de Valle D'Aosta (norte de Italia) para generar mapas de vulnerabilidad y posibles pérdidas económicas. Finalmente, se propone una metodología para obtener mapas de daños potenciales de los edificios expuestos, basados en los gradientes espaciales de movimiento, y se aplica en un área costera de la provincia de Granada (España). A partir de los métodos propuestos para extraer y clasificar las ADAs, y de otros métodos de análisis existentes, se ha desarrollado un paquete de herramientas, los ADAtools, de acceso abierto, fáciles de usar y rápidas, que optimizan la explotación operativa de los mapas de desplazamiento de escala regional. Todas las metodologías se han desarrollado en el marco de varios proyectos europeos (Safety, U-Geohaz, MOMIT y RISKCOAST), y están dirigidos a apoyar las actividades de gestión territorial y análisis de riesgos, con un enfoque específico a los deslizamientos de tierra.Postprint (published version