Application of Differential and Polarimetric Synthetic Aperture Radar (SAR) Interferometry for Studying Natural Hazards

In the following work, I address the problem of coherence loss in standard Differential Interferometric SAR (DInSAR) processing, which can result in incomplete or poor quality deformation measurements in some areas. I incorporate polarimetric information with DInSAR in a technique called Polarimetric SAR Interferometry (PolInSAR) in order to acquire more accurate and detailed maps of surface deformation. In Chapter 2, I present a standard DInSAR study of the Ahar double earthquakes (Mw=6.4 and 6.2) which occurred in northwest Iran, August 11, 2012. The DInSAR coseismic deformation map was affected by decorrelation noise. Despite this, I employed an advanced inversion technique, in combination with a Coulomb stress analysis, to find the geometry and the slip distribution on the ruptured fault plane. The analysis shows that the two earthquakes most likely occurred on a single fault, not on conjugate fault planes. This further implies that the minor strike-slip faults play more significant role in accommodating convergence stress accumulation in the northwest part of Iran. Chapter 3 presents results from the application of PolInSAR coherence optimization on quad-pol RADARSAT-2 images. The optimized solution results in the identification of a larger number of reliable measurement points, which otherwise are not recognized by the standard DInSAR technique. I further assess the quality of the optimized interferometric phase, which demonstrates an increased phase quality with respect to those phases recovered by applying standard DInSAR alone. Chapter 4 discusses results from the application of PolInSAR coherence optimization from different geometries to the study of creep on the Hayward fault and landslide motions near Berkeley, CA. The results show that the deformation rates resolved by PolInSAR are in agreement with those of standard DInSAR. I also infer that there is potential motion on a secondary fault, northeast and parallel to the Hayward fault, which may be creeping with a lower velocity

Decomposing DInSAR Time-Series into 3-D in Combination with GPS in the Case of Low Strain Rates: An Application to the Hyblean Plateau, Sicily, Italy

Differential Interferometric SAR (DInSAR) time-series techniques can be used to derive surface displacement rates with accuracies of 1 mm/year, by measuring the one-dimensional distance change between a satellite and the surface over time. However, the slanted direction of the measurements complicates interpretation of the signal, especially in regions that are subject to multiple deformation processes. The Simultaneous and Integrated Strain Tensor Estimation from Geodetic and Satellite Deformation Measurements (SISTEM) algorithm enables decomposition into a three-dimensional velocity field through joint inversion with GNSS measurements, but has never been applied to interseismic deformation where strain rates are low. Here, we apply SISTEM for the first time to detect tectonic deformation on the Hyblean Foreland Plateau in South-East Sicily. In order to increase the signal-to-noise ratio of the DInSAR data beforehand, we reduce atmospheric InSAR noise using a weather model and combine it with a multi-directional spatial filtering technique. The resultant three-dimensional velocity field allows identification of anthropogenic, as well as tectonic deformation, with sub-centimeter accuracies in areas of sufficient GPS coverage. Our enhanced method allows for a more detailed view of ongoing deformation processes as compared to the single use of either GNSS or DInSAR only and thus is suited to improve assessments of regional seismic hazard

Rapidly accelerating subsidence in the Greater Vancouver region from two decades of ERS-ENVISAT-RADARSAT-2 DInSAR measurements

Rapidly accelerating ground subsidence in the south-western part of British Columbia, the third largest metropolitan area in Canada with over 2.3 million of inhabitants, is estimated using the Multidimensional Small Baseline Subset (MSBAS) advanced Differential Interferometric Synthetic Aperture Radar (DInSAR), an effective processing strategy for multi-mission, multi-temporal SAR data. The Synthetic Aperture Radar (SAR) data used in this study consists of seven independent data sets: one ascending and one descending ERS-1/2 and ENVISAT frames, together spanning July 1995–September 2010, and three RADARSAT-2 frames spanning February 2009–October 2012. During the July 1995–October 2012 period we observe fast ground subsidence with a maximum rate of about 2 cm/year in the Greater Vancouver region that includes the Fraser River Delta and the cities of Burnaby, Richmond, Surrey, New Westminster and Vancouver. The rapidly accelerating subsidence is observed beneath the Vancouver International Airport, SkyTrain terminal as well as several agricultural and industrial locales. These time series suggest that the subsidence rate at the studied regions does not decrease with time, as suggested in previous studies, but remains steady or increases. These results also demonstrate the importance of acquiring and appropriately estimating longer time series, as previous studies on the same Greater Vancouver area may have misinterpreted the long term ground deformation rate and direction and underestimated the potential hazard. The long term impact of this subsidence on urban infrastructure can be significant and needs to be investigated further.The work of PG was supported by Banting Postdoctoral Fellowship. The work of KFT was supported by an NSERC Discovery Grant.Peer reviewe

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

Principal component analysis of MSBAS DInSAR time series from Campi Flegrei, Italy

Because of its proximity to the city of Naples and with a population of nearly 1 million people within its caldera, Campi Flegrei is one of the highest risk volcanic areas in the world. Since the last major eruption in 1538, the caldera has undergone frequent episodes of ground subsidence and uplift accompanied by seismic activity that has been interpreted as the result of a stationary, deeper source below the caldera that feeds shallower eruptions. However, the location and depth of the deeper source is not well-characterized and its relationship to current activity is poorly understood. Recently, a significant increase in the uplift rate has occurred, resulting in almost 13 cm of uplift by 2013 (De Martino et al., 2014; Samsonov et al., 2014b; Di Vito et al., 2016). Here we apply a principal component decomposition to high resolution time series from the region produced by the advanced Multidimensional SBAS DInSAR technique in order to better delineate both the deeper source and the recent shallow activity. We analyzed both a period of substantial subsidence (1993–1999) and a second of significant uplift (2007–2013) and inverted the associated vertical surface displacement for the most likely source models. Results suggest that the underlying dynamics of the caldera changed in the late 1990s, from one in which the primary signal arises from a shallow deflating source above a deeper, expanding source to one dominated by a shallow inflating source. In general, the shallow source lies between 2700 and 3400 m below the caldera while the deeper source lies at 7600 m or more in depth. The combination of principal component analysis with high resolution MSBAS time series data allows for these new insights and confirms the applicability of both to areas at risk from dynamic natural hazards.The work of PG was supported by a Banting Postdoctoral Fellowship. This research was partially funded by the NERC project Centre for the Observation and Modelling of Earthquakes, Volcanoes and Tectonics (COMET, GA/13/M/031, http://comet.nerc.ac.uk). The work of KFT was supported by an NSERC Discovery Grant and CIRES, University of Colorado Boulder. This work also was supported by the AQUARISK Spanish project and the EU MED-SUV project.Peer reviewe

Spatial and temporal patterns of deformation at the Tendaho geothermal prospect, Ethiopia

Observations of ground deformation in East Africa have been fundamental for unveiling the tectonics of continental rifting, assessing the seismic and volcanic hazard to development, and identifying geothermal resources. Here we investigate the active natural and anthropogenic processes in the Tendaho Graben, Afar using Interferometric Synthetic Aperture Radar (InSAR) collected by the Envisat satellite in 2004–2010. We used the Poly-Interferometric Rate And time series Estimation (π-RATE) method to calculate displacement in satellite line-of-sight, and a least-square inversion to decompose the line-of-sight displacement into vertical and rift perpendicular components. We observe two zones of deformation: a 20 km wide circular region of subsidence located 10 km northeast of the town of Semera with a maximum displacement rate of ∼5 cm/yr; and elongated zone (50 km) of subsidence in the area of the geothermal prospect, maximum rate of ∼4 cm/yr. The temporal characteristics of subsidence varies between these zones, with an increase in subsidence rate observed in the circular region in August 2008. We used a Bayesian inversion to find the best fitting source models and compared this to locations of seismicity and other geophysical observations. The pattern of deformation is consistent with a combination of magmatic and geothermal processes, but there does not appear to be a direct link to a sequence of dyke intrusions during 2005–2010 at Manda Hararo graben ∼60 km away, but dynamic stress changes or deep crustal flow could account for the observations

Measuring large topographic change with InSAR: lava thicknesses, extrusion rate and subsidence rate at Santiaguito volcano, Guatemala

Lava flows can produce changes in topography on the order of 10s-100s of metres. A knowledge of the resulting volume change provides evidence about the dynamics of an eruption. We present a method to measure topographic changes from the differential InSAR phase delays caused by the height differences between the current topography and a Digital Elevation Model (DEM). This does not require a pre-event SAR image, so it does not rely on interferometric phase remaining coherent during eruption and emplacement. Synthetic tests predicts that we can estimate lava thickness of as little as �9 m, given a minimum of 5 interferograms with suitably large orbital baseine separations. In the case of continuous motion, such as lava flow subsidence, we invert interferometric phase simultaneously for topographic change and displacement. We demonstrate the method using data from Santiaguito volcano, Guatemala, and measure increases in lava thickness of up to 140 m between 2000 and 2009, largely associated with activity between 2000 and 2005. We find a mean extrusion rate of 0.43 +/- 0.06 m3/s, which lies within the error bounds of the longer term extrusion rate between 1922-2000. The thickest and youngest parts of the flow deposit were shown to be subsiding at an average rate of �-6 cm/yr. This is the first time that flow thickness and subsidence have been measured simultaneously. We expect this method to be suitable for measurment of landslides and other mass flow deposits as well as lava flows

Innovative Techniques for the Retrieval of Earth’s Surface and Atmosphere Geophysical Parameters: Spaceborne Infrared/Microwave Combined Analyses

With the advent of the first satellites for Earth Observation: Landsat-1 in July 1972 and ERS-1 in May 1991, the discipline of environmental remote sensing has become, over time, increasingly fundamental for the study of phenomena characterizing the planet Earth. The goal of environmental remote sensing is to perform detailed analyses and to monitor the temporal evolution of different physical phenomena, exploiting the mechanisms of interaction between the objects that are present in an observed scene and the electromagnetic radiation detected by sensors, placed at a distance from the scene, operating at different frequencies. The analyzed physical phenomena are those related to climate change, weather forecasts, global ocean circulation, greenhouse gas profiling, earthquakes, volcanic eruptions, soil subsidence, and the effects of rapid urbanization processes. Generally, remote sensing sensors are of two primary types: active and passive. Active sensors use their own source of electromagnetic radiation to illuminate and analyze an area of interest. An active sensor emits radiation in the direction of the area to be investigated and then detects and measures the radiation that is backscattered from the objects contained in that area. Passive sensors, on the other hand, detect natural electromagnetic radiation (e.g., from the Sun in the visible band and the Earth in the infrared and microwave bands) emitted or reflected by the object contained in the observed scene. The scientific community has dedicated many resources to developing techniques to estimate, study and analyze Earth’s geophysical parameters. These techniques differ for active and passive sensors because they depend strictly on the type of the measured physical quantity. In my P.h.D. work, inversion techniques for estimating Earth’s surface and atmosphere geophysical parameters will be addressed, emphasizing methods based on machine learning (ML). In particular, the study of cloud microphysics and the characterization of Earth’s surface changes phenomenon are the critical points of this work

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

Wenchuan Earthquake Deformation 3D Modelling based on ALOS/PALSAR Data

A devastating earthquake of magnitude Mw 7.9 occurred in Wenchuan area of Sichuan Province, China on 12th May 2008 and caused great casualties and economic damage. This study is aiming to investigate the faulting geometry and motion of the major seismic faults in Longmenshan fault thrust belt that caused this earthquake, based on the surface rupture displacement data measured using differential interferometric synthetic aperture radar (DInSAR) and SAR amplitude pixel-offset techniques. The cross-event Japanese ALOS PALSAR data have been used for this study. First, the methodology for recovering the missing data in the decoherence zone of the DInSAR line-of-sight (LOS) surface motion maps was developed. In the area along the seismic fault zone, the coherence between pre- and post-event SAR images is completely lost because of the earthquake induced violent and chaotic destruction on the land surface and as the result, no surface displacement can be measured using the DInSAR technique. An Adaptive Local Kriging (ALK) technique has then been developed to retrieve the interferometric fringe patterns in the decoherence zone. The novel ALK operating in a multi-step approach enables to retrieve and interpolate the values with high fidelity to the original dataset. Thus a map of continuous radar LOS displacement was generated. Then, the horizontal displacement motion maps in ground range and azimuth direction were derived from cross-event SAR amplitude image pairs using advanced sub-pixel offset technique, Phase Correlation based Image Analysis System (PCIAS). Though the ground range pixel-offset is proportional to the LOS displacement, the azimuth pixel-offset data provide extra information of the coseismic motion. Thus the horizontal displacement vector field can be obtained in order to constrain the faulting motions in key areas. Finally, with the constraints of the ALK refined DInSAR data and the horizontal displacement data together with the published seismic focal mechanism solutions, seismic reflection profiles and field observations, forward modelling was proceeded using the Poly3D software to decide the most likely faulting geometry based on the optimal matching between the simulated and the measured surface displacement. In the much disputed Beichuan – Pengguan area, the best fit is achieved only when the Pengguan fault is set as the main fault that intercept the Yingxiu-Beichuan fault at a depth of about 13 kilometres. This geometric relationship between the two faults and the distribution of slip is compatible with them being two adjacent splay faults on a propagating thrust