Combining differential SAR interferometry and the probability integral method for three-dimensional deformation monitoring of mining areas

With the exploitation of coal resources, ground surface subsidence continues to occur in mining areas, destroying the ecological environment and significantly affecting the daily productivity and life of humans. The differential synthetic aperture radar interferometry (D-InSAR) technique is widely used to monitor ground surface deformation because of its unique advantages such as high accuracy and wide coverage. However, conventional D-InSAR technology provides only one-dimensional (1D) displacement monitoring along the radar line of sight (LOS). This article proposes a method based on an analysis of the mining subsidence law for true three-dimensional (3D) displacement monitoring by combining D-InSAR and a subsidence prediction model based on the probability integral method. In this approach, 1D displacement, obtained using D-InSAR, is then combined with the prediction model to obtain the 3D displacement of ground surface target points. Here, 3D displacement curves were obtained for the Fengfeng mining area (China) using RadarSat-2 images obtained on 9 January and 2 February 2011. True ground surface displacement was measured simultaneously by levelling when the 152under31 s working face was being exploited in Jiulong mine. Vertical displacement and inclined deformation calculated using the proposed method were compared with levelling survey data and the results showed average differences of 3.2 mm and 0.1 mm m^?1, respectively; the calculated maximum displacement in the east–west and south–north directions were 106 and 73 mm, respectively. The spatial distribution of the displacements was in accordance with the mining subsidence law. Thus, the new method can retrieve highly accurate 3D displacements caused by mining subsidence.     

Two‐scale surface deformation analysis using the SBAS‐DInSAR technique: a case study of the city of Rome,Italy

We have exploited the capability of the differential synthetic aperture radar (SAR) interferometry (DInSAR) technique, referred to as Small BAseline Subset (SBAS) approach, to analyse surface deformation at two distinct spatial scales: a low resolution, large scale, and a fine resolution, local scale. At the large scale, the technique investigates DInSAR data with a ground resolution of the order of 100 m×100 m and leads to generate mean deformation velocity maps and associated time series for areas extending to some thousands of square kilometres. At the local scale, the technique exploits the SAR images at full spatial resolution (typically of the order of 5 m×20 m), detecting and analysing localized deformation phenomena. The study is focused on the city of Rome, Italy, and we used the ERS‐1/2 satellite radar data relevant to the 1995–2000 time period. The presented results demonstrate the capability of the SBAS approach to retrieve, from the low‐resolution DInSAR data, large‐scale deformation information leading to identify several sites affected by significant displacements. Our analysis permitted us to conclude that a major contribution to the detected displacements is due to the consolidation of the alluvial soils present in the area, mostly enforced by the buildings' overload. Furthermore, in a selected area, a detailed analysis was carried out by exploiting the full resolution DInSAR data. In this case we investigated deformation phenomena at the scale of single buildings. As key result we showed that differential displacements of few mm a^?1, affecting single man‐made structures or building complexes, could be detected, thus allowing to identify sites that may potentially be involved in critical situations.     

Integration of TerraSAR-X and PALSAR PSI for detecting ground deformation

In connection with the detection of various spatial- and temporal-scale ground settlements, an integrated persistent scatterer interferometry (PSI) approach is discussed using multi-source, multi-temporal, and multi-resolution synthetic aperture radar (SAR) data. Based on the comprehensive analysis of characteristics of available radar sensors, two remote-sensing SAR data sets were selected: 1 m resolution X-band TerraSAR-X and 10 m resolution L-band Advanced Land Observing Satellite (ALOS) phased array L-band SAR. ‘Tianjin Binhai New Area’ has become one of the most important economic centres in China, and one of its fast-developing urban areas, Tanggu, was selected as the study area. PSI processing was conducted on both data sets. Substantial validation was performed for PSI results from both data sources using levelling measurement. The overall good agreement confirmed the ground deformation maps derived from both data sets. Integration of PSI results appears to be a potentially significant contribution to solving the problems related to common spatial and temporal gaps when using single-type data sets. Application of both data sets revealed the capability of integrated PSIs to measure ground deformation with strong temporal and spatial variation, thereby improving the interpretation of ground deformation characteristics which increases the confidence of hazard assessment and provides some insight into complex underlying mechanisms.     

Evolution of urban monitoring with radar interferometry in Madrid City: performance of ERS-1/ERS-2, ENVISAT,COSMO-SkyMed,and Sentinel-1 products

During the last decades, synthetic aperture radar (SAR) image exploitation has matured with the launch of different satellite missions and the development of different techniques, which allow exploiting the capabilities of the radar images. Among these techniques, persistent scatterer interferometry (PSI) has proven to be a powerful tool to derive terrain deformations over urban areas. It is based on the use of a large number of images over wide areas in order to obtain terrain displacements time series. The imagery from the different SAR missions has led to an archive with data that covers up to 30 years in the past. Moreover, different methods and algorithms have been proposed in order to perform this complex task. In this line, this work aims at identifying if data from different missions and processed by different techniques can be combined in order to study the evolution of urban monitoring. Three different PSI techniques are used in order to process data from four SAR missions: European Remote Sensing (ERS)-1/2, Environmental Satellite, COSMO-SkyMed, and the recent Sentinel-1 A/B. The rapidly evolving urban area of Madrid, where numerous undergrounding works have been carried out in the last decade, has been chosen as the testing environment. The density of persistent scatterers, the deformation accuracy validated with GPS displacements and deformation trends are used as the key performance items for the assessment.     

Land deformation monitoring using coherent target-neighbourhood networking method combined with polarimetric information: a case study of Shanghai,China

In this article, an advanced approach for land deformation monitoring using synthetic aperture radar (SAR) interferometry combined with polarimetric information is presented. The linear and nonlinear components of the deformation, the error of the digital elevation model (DEM) and the atmospheric artefacts can be achieved by a coherent target (CT)-neighbourhood networking approach. In order to detect recent land deformation in Shanghai, China, 12 ENVISAT advanced synthetic aperture radar (ASAR) alternating polarization images acquired from January 2006 to August 2008 are employed for deformation analysis. Over a 2.5-year period, two deformation velocity fields from HH and VV modes over Shanghai are derived using the CT-neighbourhood networking SAR interferometry (InSAR), then integrated into a final deformation map by a fusion scheme. It is found that the annual subsidence rates in the study area range from??20 to 10 mm year^?1 and the average subsidence rate in the downtown area reaches??7.5 mm year^?1, which is consistent with the local government statistics published in 2007.     

Landslide monitoring by corner reflectors differential interferometry SAR

Although the differential interferometric synthetic aperture radar (DInSAR) technique has the potential for monitoring ground deformation with millimetric accuracy, it still suffers from temporal and geometric decorrelation. In this paper, the corner reflectors differential interferometry synthetic aperture radar (CRDInSAR) technique was used to overcome the limitations of conventional DInSAR. We studied the basic principles of CRDInSAR, discussed the calculation of the flat earth and topographic phases based on the geometry of satellite and corner reflectors, presented the phase unwrapping approach for the sparse grid of corner reflectors, then investigated the construction of, and solution to, the unwrapped phase system equation. Subsequently, we applied CRDInSAR to monitor the displacement of the Shuping landslide in Hubei Province, China. In this study, we processed five SAR images on a descending pass acquired by the Environmental Satellite (ENVISAT) Advanced Synthetic Aperture Radar (ASAR) sensor from September 2005 to March 2006, and compared the achieved results of CRDInSAR with Global Positioning System (GPS) measurements at the same time of the SAR data acquisition assumed as reference. The result indicates a good agreement between the measurements provided by the two different techniques, which shows that CRDInSAR allows monitoring of slow landslide deformation in low coherence areas and provides accurate results.     

A method of monitoring three-dimensional ground displacement in mining areas by integrating multiple InSAR methods

Underground mining often causes large-gradient vertical and horizontal displacement in ground surface, resulting in the losses of coherence and difficulties in phase unwrapping. This article presents an approach to determine three-dimensional ground displacements in mining areas with the large gradient or phase decorrelation by integrating multiple interferometric synthetic aperture radar (InSAR) methods. The core of the proposed method is that the offset-tracking method is employed to solve for the displacement with the large gradient or phase decorrelation. First, the displacements in the radar line-of-sight directions are obtained from two interferometric pairs with different viewing geometries by integrating the measurements of differential InSAR and offset tracking. Then, the displacements in the azimuthal directions are obtained from two interferometric pairs with different viewing geometries by integrating the measurements of multiple aperture interferometry and offset tracking. Finally, the three-dimensional ground displacement fields are inferred from these four independent, one-dimensional displacements using the least squares method and Helmert variance component estimation. We apply this method to obtain the three-dimensional ground displacement field in the Dongtan mine region. We compare the results with those of levelling and global positioning system surveys, and the root mean square errors of the results were 24 and 43 mm in the vertical direction and horizontal directions, respectively. The experimental results indicate that the proposed method can be used to estimate three-dimensional ground displacement fields in mining areas with large-gradient displacement and phase decorrelation.     

Post-emplacement lava subsidence and the accuracy of ERS InSAR digital elevation models of volcanoes

Repeat-pass synthetic aperture radar interferometry (InSAR) using data acquired by the ERS platforms is an attractive method for acquiring topographic data of volcanoes. Caution is advised, however, when using this technique in regions covered by young, thick lava flows. In this study, the magnitude of post-emplacement subsidence associated with the 1991-93 lava flow at Mount Etna, Sicily, was measured using differential radar interferometric techniques, and it was found that the rates of subsidence are large enough to contribute a significant component to the measured phase shift, even in ERS data acquired on consecutive orbits. It demonstrates the detrimental effect that such phase shifts have on the accuracy of digital elevation models derived by repeat-pass radar interferometry.     

An approach for accurately retrieving the vertical deformation component from two-track InSAR measurements

As is well known, both conventional differential synthetic aperture radar interferometry (D-InSAR) and multi-temporal synthetic aperture radar interferometry (MT-InSAR) have a common limitation that they only can measure the deformation component along the radar line of sight (LOS) direction. However, in the majority of disaster investigations, there is more interest in the vertical deformation component than that of the horizontal direction, for example, in measuring ground subsidence in urban areas, or ground subsidence due to underground mining. To estimate the vertical deformation component accurately, it is in theory necessary to create at least three independent equations, and solve the vertical, the North–South, and the East–West deformation components by exploiting at least three-track InSAR LOS measurements or combining at least two-track InSAR LOS measurements with azimuth measurements. However, these methods are greatly limited and sometimes not even practical because there is typically little chance of obtaining a three-track SAR data set covering the same area in the same time span, and the accuracy of the azimuth measurements is far lower than that of InSAR LOS measurements. In this article, we found that it is possible to solve for the vertical deformation component from two-track InSAR LOS measurements in some circumstances. Then, an approach for accurately retrieving the vertical deformation component from two-track InSAR LOS measurements is proposed, and the analytical expression is presented. The approach is illustrated through an investigation of the ground subsidence in an area of Beijing, China. Unlike previous methods, this approach can accurately retrieve the vertical deformation component from two-track InSAR LOS measurements, and provide more reliable results for improving the interpretation of ground subsidence phenomena.     

Monitoring and analysis of mining 3D time-series deformation based on multi-track SAR data

Accurate monitoring of the developing process of a surface subsidence basin is the basis of building damage assessment and surface deformation prediction. In this paper, the Synthetic Aperture Radar (SAR) data of three different imaging geometries, TerraSAR, Radarsat-2, and Sentinel-1A, were exploited. Firstly, two-dimensional (2D) time-series deformation of the surface subsidence basin caused by 15,235 working face mining was obtained based on Multidimensional Small Baseline Subset (MSBAS) technology from 19 December 2015 to 5 March 2016. By comparing vertical deformation with levelling data, it is shown that the root-mean-square error of vertical deformation is 3.2 mm and the standard deviation is 1.9 mm when the ascending-descending track SAR data is available. Otherwise, the root-mean-square error of vertical deformation is 18.1 mm and the standard deviation is 11.6 mm. Because of the low precision of the north–south horizontal movement monitored by the SAR sensor, the vertical deformation acquired by MSBAS technology and the rules of the mining subsidence (horizontal movement is proportional to tilt) were combined to obtain the north–south horizontal movement which was proven to be reliable by comparing the 2D time-series deformation obtained by MSBAS technology. Then, the deformation of the railway in the surface subsidence basin was analysed based on the three-dimensional (3D) time-series deformation. The results show that the subsidence, tilt, and horizontal movement strongly influence the railway in the monitoring period, but will not affect the normal traffic. This experiment lays a technical foundation for preventing the occurrence of mining disasters and verifies the ability to monitor the deformation of buildings and structures by interferometry synthetic aperture radar technology.     

A method for minimising of low frequency and unwrapping artefacts from interferometric calculations

On 7 September 1999 a magnitude M_W = 5.9 earthquake occurred in the Athens area of Greece producing a subsidence of 6–7 cm detected by radar interferometry. This study introduces a processing technique, which produces a clear deformation pattern of the earthquake, mostly released from artefacts due to orbital effects, unwrapping cycle slipping errors and atmospherics disturbances. A set of 17 ERS‐1 and ERS‐2 SAR images acquired between December 1997 and January 2001 has been used. The contribution of each artefact to the interferograms was calculated, both in the frequency and spatial domains.     

Monitoring land deformation in Changzhou city (China) with multi-band InSAR data sets from 2006 to 2012

Over exploitation of groundwater in Changzhou city, China can cause land deformation, which in turn proves detrimental to the urban infrastructure. In this study, multi-band synthetic aperture radar (SAR) data sets (C-band Envisat ASAR, L-band ALOS PALSAR, and X-band COSMO-SkyMed) acquired from 2006 to 2012 were analysed using the synthetic aperture radar (SAR) interferometry (InSAR) time-series method to investigate the relationship between spatial–temporal distribution of land deformation and groundwater exploitation. Annual deformation rate inferred from multi-band interferograms ranges from ?58 to 24 mm year^?1. Levelling-survey data were used to validate the multi-band InSAR measurements. The results showed that these two types of measurements were generally in agreement. Correlating groundwater-table and multi-band InSAR measurements at six groundwater-well stations showed that with the rise of the water table, the land rebounded. But in some areas with larger subsidence, continual subsidence was observed even though the water table rose after the prohibition of groundwater exploitation. This may have been caused by the hysteresis effect due to the consolidation of strata (especially for the creep deformation). Our study provides scientific evidence on the management of groundwater extraction and the assessment of land-subsidence hazards.     

Detection of mining related ground instabilities using the Permanent Scatterers technique—a case study in the east of France

Ground stability is a major concern for land use planning and both natural and anthropogenic risk assessment, especially in urbanized areas. Space‐borne differential radar interferometry provides a unique tool able to give a synoptic view of ground deformation with centimetric to millimetric vertical precision. Approaches for combining a wide range of radar images such as the permanent scatterers (PS) technique allow the estimation of the deformation history of single buildings. The PS approach has been exploited to investigate a test site particularly exposed to ground deformation hazards, namely the iron mining basin in Lorraine (France). In this Letter, a specific focus was set on the case of Roncourt, where precursor signs of a collapse affecting an area of ～300×300?m² have been identified.     

The present status of subsiding land vulnerable to roof collapse in the Jharia Coalfield,India, as obtained from shorter temporal baseline C-band DInSAR by smaller spatial subset unwrapped phase profiling

In this paper, we identified recently subsiding areas in Jharia Coalfield, Jharkhand, India from the shorter temporal baseline Radarsat-2 C-band interferometric synthetic aperture radar (InSAR) data pairs of 2012. Although shorter wavelength C-band differential InSAR (DInSAR) is more sensitive to slow deformation and better suited for higher precision land subsidence measurement, the dynamic and adverse land cover in mining areas and resulting temporal decorrelation problem poses a serious problem for DInSAR observation in mining areas. We used smaller temporal baseline data pairs and adopted InSAR coherence-guided incremental filtering with smaller moving windows to highlight the deformation fringes over temporal decorrelation noise. We identified the deformation fringes and validated them based on ground information to prepare the land subsidence map of the coalfield in 2012. Several new, previously unreported subsidence areas were detected in the present study with a total subsiding area of 6.9 km². The recent incidence of roof collapse on 15 November 2014 at Angar Patra village in Katras region of the coalfield where 45 houses collapsed and 10 people were injured is situated in a highly subsiding vulnerable area as obtained from the present study. Due to spatial discontinuities of InSAR coherence, DInSAR phase unwrapping for the entire study area in one go did not appear feasible. To avoid this problem, we performed DInSAR processing in smaller spatial subsets and unwrapping of the subset interferograms by a ‘minimum cost flow’ algorithm. Subsequently, we plotted unwrapped phase profiles across the deformation fringes and retrieved the maximum deformation phase with respect to background phase and translated them into radar line of sight (LOS) displacement rates. For obtaining the average subsidence rates, we adopted InSAR coherence-weighted LOS displacement rates taking into account the contribution of each data pair as a function of DInSAR phase quality of the fringe areas. Ground-based subsidence measurements by precision levelling were conducted in four test sites that had been undergoing active underground mining during the observation period. We compared space-borne DInSAR-based subsidence rates obtained by the adopted technique with precision levelling measurements. Overall, the results are found to agree well. In the four test sites with gentle to flat topography, land subsidence occurs at slow to moderate rates due to compression of in-filled material (resulting from sand stowing in underground mining), without any evidence of roof collapse. In such cases, the horizontal displacement component is less significant, and overall surface displacement occurs essentially in the vertical direction. However, we assessed the nature of subtle horizontal strain to infer relative shrinkage or dilation of the land surface which could be additive or subtractive to vertical displacement in DInSAR-based LOS displacement.     

Monitoring of recent ground surface subsidence in the Cangzhou region by the use of the InSAR time-series technique with multi-orbit Sentinel-1 TOPS imagery

ABSTRACT

By collecting 39 scenes of descending images and 26 scenes of ascending images of Sentinel-1 synthetic aperture radar (SAR), we monitored the ground subsidence situation of Cangzhou in China during the period from March 2015 to February 2017 with the advanced synthetic aperture radar interferometry (InSAR) time-series technique and obtained the time-series subsidence rates of Cangzhou. We then selected the two sets of results of the monitoring obtained during the same period (from July 2015 to December 2016) to verify the results’ accuracy, considering three overlapping areas (Xinhua District, Botou County, and Dongguang County). This analysis clearly indicated that both types of results have good consistency, and the maximum subsidence occurred in Dongguang County. By further study of the central area of Dongguang and the related urban development, we found that the subsidence reached about 80.0 mm over the study period and there was a close relationship between the subsidence trend and the main direction of the city development. Moreover, by combining the two sets of results, we confirmed that there has been subsidence of the high-speed railway line in the whole of the Cangzhou area, among which the most obvious subsidence has occurred in Dongguang and Qing Counties. Finally, it was proved that the Sentinel-1 data can be used to monitor ground surface subsidence, and the data are especially effective in identifying persistent scatterer points along a linear feature. Therefore, this article could provide reliable data to assist with important decisions about urban development projects in the Cangzhou area in the next few years.     

Mapping hydrologically sensitive areas on the Boreal Plain: a multitemporal analysis of ERS synthetic aperture radar data

Characterizing the spatial and temporal dynamics of hydrologically sensitive areas (HSAs) is vital to the effective management of the boreal forest. HSAs are defined as saturated or inundated areas that, if disturbed, might result in a significant change in the movement of water, nutrients and biota within landscapes. This study presents a remote sensing technique that uses archived European Remote Sensing Satellite (ERS)‐1 and ERS‐2 synthetic aperture radar (SAR) images to monitor HSAs in the Willow River watershed (1030 km²) on the western Boreal Plain of Canada. ERS images were used to generate a probability of HSA occurrence map for a 10‐year period (1991–2000). This map revealed the complexity of HSAs on the western Boreal Plain, where some areas remained consistently dry or wet whereas others were dynamic, transitioning from dry to wet and vice versa. A probability map of HSA occurrence provides spatial and temporal information previously unavailable for this region that may expand our understanding of the hydrological behaviour of drainage basins and serve as a planning tool for land management decisions.     

Long-term ERS/ENVISAT deformation time-series generation at full spatial resolution via the extended SBAS technique

We extend the small baseline subset (SBAS) differential synthetic aperture radar (SAR) interferometry (DInSAR) approach to allow the generation of deformation time-series by processing, at the full spatial resolution scale, long sequences of European Remote Sensing (ERS-1 and ERS-2) and Environmental Satellite (ENVISAT) SAR data acquired with the same illumination geometry. In particular, we avoid the generation of ERS/ENVISAT cross-interferograms, which are severely affected by noise phenomena due to the carrier frequency separation of the two SAR systems, and we focus on single-platform interferograms only (i.e. ERS/ERS and ENVISAT/ENVISAT interferograms) that are properly combined by applying the singular value decomposition (SVD)-based SBAS approach. Moreover, we exploit the Doppler centroid variations of the post-2000 acquisitions of the ERS-2 sensor and the carrier frequency difference between the ERS-1/2 and the ENVISAT systems, in order to maximize the number of investigated SAR pixels and to improve their geocoding. The presented results, achieved on two data sets relevant to the Napoli Bay area and to the Murge region, both located in southern Italy, confirm the effectiveness of the extended SBAS technique and demonstrate the relevance of deformation analysis carried out at the scale of single buildings or human-made structures with more than 15 years of ERS and ENVISAT acquisitions.     

Potential of diffuse scatterer interferometry for monitoring CO2 storage sites in European contexts (land cover types)

This article proposes to test the feasibility of long-term surface deformation monitoring based on synthetic aperture radar (SAR) interferometry on carbon dioxide (CO₂) storage sites with land cover representative of potential European injection sites (agricultural or forests with minimum built-up land cover). Because no operational injection site is currently active in Europe, a SAR data set (based on EnviSAT–ASAR spaceborne data) is simulated by combining SAR scenes acquired over a potential future European injection site with deformation measurements from SAR analysis carried out on the In-Salah (Algeria) CO₂ injection demonstrator site. The study shows that under such conditions, both persistent scatterer interferometry (PSI) and diffuse scatterer (DS) interferometry appear insufficient to provide a sufficiently dense measurement network to characterize surface deformation correctly. Alternative solutions, to be investigated in further studies, include the use of data archives with shorter acquisition time spans (e.g. Sentinel-1 data when available) or installation of corner reflectors. The cost of the latter mixed space/ground solution must be evaluated with respect to conventional ground-based measurement methods in the proposed context.     

Assessing the activity of a large landslide in southern Italy by ground-monitoring and SAR interferometric techniques

Landslides are recognized as one of the most damaging natural hazards in Italy. Campania region represents a complex geological setting, where mass movements of different types are widespread, and urban expansion can be increasingly seen by the presence of buildings and infrastructure in landslide-prone areas. In such a context, monitoring of unstable slopes represents a key activity in the process of landslide risk prevention and mitigation, in order to correctly establish a cause–effect correlation and to predict the possible reactivation phases that may result in high costs for the human society. This article focuses on the application of different methods of landslide analysis and monitoring, including those developed more recently and based on data acquired by satellites and processed by synthetic aperture radar (SAR) interferometric techniques. The study area is a small town, Calitri, known worldwide for the large landslide reactivated by the 23 November 1980 earthquake that destroyed a large sector of the historical centre. The site has been investigated by two ground-monitoring campaigns, the analysis of which allowed identification of the evolution of landslide activity over time. Furthermore, differential SAR interferometry (DInSAR), based upon two different approaches, allowed us to produce point-wise and wide area deformation maps after processing data sets of Earth Resource Satellite 1/2 (ERS-1/2) images, respectively acquired in 1992–2001 and 1992–1995. The results obtained from this analysis highlighted the potentiality of remote-sensing tools in landslide hazard assessment and led to development of a research project based on the installation of corner reflectors along unstable slopes and aimed at creating a field–Earth observation monitoring system.