Coherency Matrix Decomposition-Based Polarimetric Persistent Scatterer Interferometry

© 2019 IEEE. Personal use of this material is permitted. Permission from IEEE must be obtained for all other uses, in any current or future media, including reprinting/republishing this material for advertising or promotional purposes,creating new collective works, for resale or redistribution to servers or lists, or reuse of any copyrighted component of this work in other works.The rationale of polarimetric optimization techniques is to enhance the phase quality of the interferograms by combining adequately the different polarization channels available to produce an improved one. Different approaches have been proposed for polarimetric persistent scatterer interferometry (PolPSI). They range from the simple and computationally efficient BEST, where, for each pixel, the polarimetric channel with the best response in terms of phase quality is selected, to those with high-computational burden like the equal scattering mechanism (ESM) and the suboptimum scattering mechanism (SOM). BEST is fast and simple, but it does not fully exploit the potentials of polarimetry. On the other side, ESM explores all the space of solutions and finds the optimal one but with a very high-computational burden. A new PolPSI algorithm, named coherency matrix decomposition-based PolPSI (CMD-PolPSI), is proposed to achieve a compromise between phase optimization and computational cost. Its core idea is utilizing the polarimetric synthetic aperture radar (PolSAR) coherency matrix decomposition to determine the optimal polarization channel for each pixel. Three different PolSAR image sets of both full- (Barcelona) and dual-polarization (Murcia and Mexico City) are used to evaluate the performance of CMD-PolPSI. The results show that CMD-PolPSI presents better optimization results than the BEST method by using either $D_{\mathrm{ A}}$ or temporal mean coherence as phase quality metrics. Compared with the ESM algorithm, CMD-PolPSI is 255 times faster but its performance is not optimal. The influence of the number of available polarization channels and pixel's resolutions on the CMD-PolPSI performance is also discussed.Peer ReviewedPostprint (author's final draft

Long-term monitoring of geodynamic surface deformation using SAR interferometry

Thesis (Ph.D.) University of Alaska Fairbanks, 2014Synthetic Aperture Radar Interferometry (InSAR) is a powerful tool to measure surface deformation and is well suited for surveying active volcanoes using historical and existing satellites. However, the value and applicability of InSAR for geodynamic monitoring problems is limited by the influence of temporal decorrelation and electromagnetic path delay variations in the atmosphere, both of which reduce the sensitivity and accuracy of the technique. The aim of this PhD thesis research is: how to optimize the quantity and quality of deformation signals extracted from InSAR stacks that contain only a low number of images in order to facilitate volcano monitoring and the study of their geophysical signatures. In particular, the focus is on methods of mitigating atmospheric artifacts in interferograms by combining time-series InSAR techniques and external atmospheric delay maps derived by Numerical Weather Prediction (NWP) models. In the first chapter of the thesis, the potential of the NWP Weather Research & Forecasting (WRF) model for InSAR data correction has been studied extensively. Forecasted atmospheric delays derived from operational High Resolution Rapid Refresh for the Alaska region (HRRRAK) products have been compared to radiosonding measurements in the first chapter. The result suggests that the HRRR-AK operational products are a good data source for correcting atmospheric delays in spaceborne geodetic radar observations, if the geophysical signal to be observed is larger than 20 mm. In the second chapter, an advanced method for integrating NWP products into the time series InSAR workflow is developed. The efficiency of the algorithm is tested via simulated data experiments, which demonstrate the method outperforms other more conventional methods. In Chapter 3, a geophysical case study is performed by applying the developed algorithm to the active volcanoes of Unimak Island Alaska (Westdahl, Fisher and Shishaldin) for long term volcano deformation monitoring. The volcano source location at Westdahl is determined to be approx. 7 km below sea level and approx. 3.5 km north of the Westdahl peak. This study demonstrates that Fisher caldera has had continuous subsidence over more than 10 years and there is no evident deformation signal around Shishaldin peak.Chapter 1. Performance of the High Resolution Atmospheric Model HRRR-AK for Correcting Geodetic Observations from Spaceborne Radars -- Chapter 2. Robust atmospheric filtering of InSAR data based on numerical weather prediction models -- Chapter 3. Subtle motion long term monitoring of Unimak Island from 2003 to 2010 by advanced time series SAR interferometry -- Chapter 4. Conclusion and future work

Mexico City land subsidence in 2014-2015 with Sentinel-1 IW TOPS: results using the Intermittent SBAS (ISBAS) technique

Differential Interferometric Synthetic Aperture Radar (DInSAR) can be considered as an efficient and cost effective technique for monitoring land subsidence due to its large spatial coverage and high accuracy provided. The recent commissioning of the first Sentinel-1 satellite offers improved support to operational surveys using DInSAR due to regular observations from a wide-area product. In this paper we show the results of an intermittent small-baseline subset (ISBAS) time-series analysis of 18 Interferometric Wide swath (IW) products of a 39,000 km2 area of Mexico acquired between 3 October 2014 and 7 May 2015 using the Terrain Observation with Progressive Scans in azimuth (TOPS) imaging mode. The ISBAS processing was based upon the analysis of 143 small-baseline differential interferograms. After the debursting, merging and deramping steps necessary to process Sentinel-1 IW roducts, the method followed a standard approach to the DInSAR analysis. The Sentinel-1 ISBAS results confirm the magnitude and extent of the deformation that was observed in Mexico City, Chalco, Ciudad Nezahualcóyotl and Iztapalapa by other C-band and L-band DInSAR studies during the 1990s and 2000s. Subsidence velocities from the Sentinel-1 analysis are, in places, in excess of -24 cm/year along the satellite line-of-sight, equivalent to over ~-40 cm/year vertical rates. This paper demonstrates the potential of Sentinel-1 IW TOPS imagery to support wide-area DInSAR surveys over what is a very large and diverse area in terms of land cover and topography

Monitoring land subsidence of airport using InSAR time-series techniques with atmospheric and orbital error corrections

Land subsidence is one of the common geological hazards worldwide and mostly caused by human activities including the construction of massive infrastructures. Large infrastructure such as airport is susceptible to land subsidence due to several factors. Therefore, monitoring of the land subsidence at airport is crucial in order to prevent undesirable loss of property and life. Remote sensing technique, especially Interferometric Synthetic Aperture Radar (InSAR) has been successfully applied to measure the surface deformation over the past few decades although atmospheric artefact and orbital errors are still a concerning issue in this measurement technique. Multi-temporal InSAR, an extension of InSAR technique, uses large sets of SAR scenes to investigate the temporal evolution of surface deformation and mitigate errors found in a single interferogram. This study investigates the long-term land subsidence of the Kuala Lumpur International Airport (KLIA), Malaysia and Singapore Changi Airport (SCA), Singapore by using two multi-temporal InSAR techniques like Small Baseline Subset (SBAS) and Multiscale InSAR Time Series (MInTS). General InSAR processing was conducted to generate interferogram using ALOS PALSAR data from 2007 until 2011. Atmospheric and orbital corrections were carried out for all interferograms using weather model, namely European Centre for Medium Range Weather Forecasting (ECMWF) and Network De-Ramping technique respectively before estimating the time series land subsidence. The results show variation of subsidence with respect to corrections (atmospheric and orbital) as well as difference between multi-temporal InSAR techniques (SBAS and MInTS) used. After applying both corrections, a subsidence ranging from 2 to 17 mm/yr was found at all the selected areas at the KLIA. Meanwhile, for SCA, a subsidence of about less than 10 mm/yr was found. Furthermore, a comparison between two techniques (SBAS and MInTS) show a difference rate of subsidence of about less than 1 mm/yr for both study area. SBAS technique shows more linear result as compared to the MInTS technique which shows slightly scattering pattern but both techniques show a similar trend of surface deformation in both study sites. No drastic deformation was observed in these two study sites and slight deformation was detected which about less than 20mm/yr for both study areas probably occurred due to several reasons including conversion of the land use from agricultural land, land reclamation process and also poor construction. This study proved that InSAR time series surface deformation measurement techniques are useful as well as capable to monitor deformation of large infrastructure such as airport and as an alternative to costly conventional ground measurement for infrastructure monitoring

Time Series Analysis of Surface Deformation Associated With Fluid Injection and Induced Seismicity in Timpson, Texas Using DInSAR Methods

In recent years, a rise in unconventional oil and gas production in North America has been linked to an increase in seismicity rate in these regions (Ellsworth, 2013). As fluid is pumped into deep formations, the state of stress within the subsurface changes, potentially reactivating pre-existing faults and/or causing subsidence or uplift of the surface. Therefore, hydraulic fracturing and/or fluid disposal injection can significantly increase the seismic hazard to communities and structures surrounding the injection sites (Barnhart et al., 2014). On 17th May 2012 an Mw4.8 earthquake occurred near Timpson, TX and has been linked with wastewater injection operations in the area (Shirzaei et al., 2016). This study aims to spatiotemporally relate, wastewater injection operations to seismicity near Timpson using differential interferometric synthetic aperture radar (DInSAR) analysis. Results are presented as a set of time series, produced using the Multidimensional Small Baseline Subset (MSBAS) InSAR technique, revealing two-dimensional surface deformation

Monitoring von Hangbewegungen mit InSAR Techniken im Gebiet Ciloto, Indonesien

In this doctoral thesis, the InSAR techniques are applied to detect the ground movement phenomenon and to assess the InSAR result geometrically in the Ciloto area, Indonesia. Mainly, one of those techniques, the SB-SDFP algorithm, overcomes the limitations of conventional InSAR in monitoring rural and agricultural areas and can observe extremely slow landslides. The InSAR strategy is positively known as a promising option to detect and quantify the kinematics of active landslides on a large areal scale. To minimize the bias of the InSAR displacement result, the correction of the tropospheric phase delay was carried out in a first step. This procedure is demonstrated in experiments both in the small study area in Ciloto and in a larger area. The latter is an area located in Northern Baja California, Mexico and is dominated by tectonic activity as well as groundwater-induced subsidence. A detailed investigation of the slope movement's behavior in the Ciloto district was conducted utilizing multi-temporal and multi-band SAR data from ERS1/2 (1996-1999), ALOS PALSAR (2007-2009) and Sentinel-1 (2014-2018) satellites. The region was successfully identified as a permanent active landslide prone area, especially in the vicinity of the Puncak Pass and Puncak Highway. The full 3D velocity field and the displacement time series were estimated using the inversion model. The velocity rate was classified from extremely slow to slow movement. To comprehend the landslide's behavior, a further examination of the relationship between InSAR results and physical characteristics of the area was carried out. For the long period of a slow-moving landslide, the relationship between precipitation and displacement trend shows a weak correlation. It is concluded that the extremely slow to slow deformation is not directly influenced by the rainfall intensity, yet it effectuates the subsurface and the groundwater flow. The run-off process with rainfall exceeding a soil's infiltration capacity was suspected as the main driver of the slow ground movement phenomenon. However, when analyzing rapid and extremely rapid landslide events at Puncak Pass, a significant increase in the correlation coefficient between precipitation and displacement rate could be observed.In dieser Doktorarbeit wird die Anwendung von erweiterten Verarbeitungsstrategien von InSAR Daten zur Erkennung und geometrischen Bewertung der Bodenbewegungen im Ciloto - Indonesien dargestellt. Dieser Ansatz überwindet die Beschränkungen konventioneller SAR-Interferometrie und ermöglicht sowohl ein kontinuierliches Monitoring dieses landwirtschaftich geprägten Gebietes als auch die Erfassung extrem langsamer Hangrutschungen. Um eine Verzerrung der InSAR Deformationsergebnisse zu minimieren, wurde zunächst eine Korrektur der troposphärischen Phase durchgeführt. Diese neuartige Strategie wird sowohl im Forschungsgebiet Ciloto als auch an einem größeren Gebiet demonstriert. Bei letzterem handelt es sich um einen Küstenstreifen im nördlichen Niederkalifornien, Mexiko, welcher durch hohe tektonische Aktivität und grundwasserinduzierte Landsetzungen charakterisiert ist. Die detaillierte Untersuchung des Verhaltens von Hangrutschungen im Ciloto erfolgte durch die Verarbeitung multi-temporaler SAR-Daten unter Nutzung verschiedener Frequenzbänder, darunter ESR1/2 (1996-1999), ALOS PALSAR (2007-2009) und Sentinel-1 (2014-2018) Daten. Die Region konnte erfolgreich als permanent aktives Hangrutschungsgebiet identifiziert werden, wobei der Puncak Pass und der Puncak Highway ein erhöhtes Gefahrenpotential aufweisen. Ein 3D- Geschwindig-keitsfeld der Deformation und die zugehörigen Zeitreihen wurden mit dem Inversionsmodell berechnet. Die Geschwindigkeitsrate wurde als langsam bis extrem langsam klassifiziert. Um das dynamische Verhalten der Hangrutschung zu verstehen wurde, in einer weiteren Untersuchung die Beziehung zwischen dem InSAR-Ergebnis und den physikalischen Begebenheiten im Forschungsgebiet analysiert. Es wird der Schluss gezogen, dass die langsame bis extrem langsame Verformung nicht direkt von der Niederschlagsintensität beeinflusst wird, diese sich aber auf den Untergrund und die Grundwasserströmung auswirkt. Es wird vermutet, dass der Oberflächenablauf, welcher die Infiltrationskapazität des Bodens übersteigt, ausschlaggebend für das Phänomen der langsamen Bodenbewegung ist. Für die schnellen und extrem schnellen Hangrutschungen jedoch konnte eine signifikante Erhöhung des Korrelationskoeffizienten zwischen Niederschlag und Verschiebungsrate bei Untersuchungen der Hangrutschung am Puncak-Pass nachgewiesen werden

A Study of Types of Sensors used in Remote Sensing

Of late, the science of Remote Sensing has been gaining a lot of interest and attention due to its wide variety of applications. Remotely sensed data can be used in various fields such as medicine, agriculture, engineering, weather forecasting, military tactics, disaster management etc. only to name a few. This article presents a study of the two categories of sensors namely optical and microwave which are used for remotely sensing the occurrence of disasters such as earthquakes, floods, landslides, avalanches, tropical cyclones and suspicious movements. The remotely sensed data acquired either through satellites or through ground based- synthetic aperture radar systems could be used to avert or mitigate a disaster or to perform a post-disaster analysis

Urban Deformation Monitoring using Persistent Scatterer Interferometry and SAR tomography

This book focuses on remote sensing for urban deformation monitoring. In particular, it highlights how deformation monitoring in urban areas can be carried out using Persistent Scatterer Interferometry (PSI) and Synthetic Aperture Radar (SAR) Tomography (TomoSAR). Several contributions show the capabilities of Interferometric SAR (InSAR) and PSI techniques for urban deformation monitoring. Some of them show the advantages of TomoSAR in un-mixing multiple scatterers for urban mapping and monitoring. This book is dedicated to the technical and scientific community interested in urban applications. It is useful for choosing the appropriate technique and gaining an assessment of the expected performance. The book will also be useful to researchers, as it provides information on the state-of-the-art and new trends in this fiel

Comparison of Small Baseline Interferometric SAR Processors for Estimating Ground Deformation

The small Baseline Synthetic Aperture Radar (SAR) Interferometry (SBI) technique has been widely and successfully applied in various ground deformation monitoring applications. Over the last decade, a variety of SBI algorithms have been developed based on the same fundamental concepts. Recently developed SBI toolboxes provide an open environment for researchers to apply different SBI methods for various purposes. However, there has been no thorough discussion that compares the particular characteristics of different SBI methods and their corresponding performance in ground deformation reconstruction. Thus, two SBI toolboxes that implement a total of four SBI algorithms were selected for comparison. This study discusses and summarizes the main differences, pros and cons of these four SBI implementations, which could help users to choose a suitable SBI method for their specific application. The study focuses on exploring the suitability of each SBI module under various data set conditions, including small/large number of interferograms, the presence or absence of larger time gaps, urban/vegetation ground coverage, and temporally regular/irregular ground displacement with multiple spatial scales. Within this paper we discuss the corresponding theoretical background of each SBI method. We present a performance analysis of these SBI modules based on two real data sets characterized by different environmental and surface deformation conditions. The study shows that all four SBI processors are capable of generating similar ground deformation results when the data set has sufficient temporal sampling and a stable ground backscatter mechanism like urban area. Strengths and limitations of different SBI processors were analyzed based on data set configuration and environmental conditions and are summarized in this paper to guide future users of SBI techniques