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

Insar Role in the Study of Earth's Surface and Synergic Use with Other Geodetic Data: the 2014 South Napa Earthquake

This work focuses on the role of SAR Interferometry (InSAR) in the study of many phenomena characterizing the Earth's surface. We propose an advanced integration method in order to merge the InSAR data with other geodetic data, i.e. Multiple Aperture Interferometry (MAI), Pixel Offset Tracking (POT) and Global Positioning System (GPS). We apply the method to constrain the full 3D displacement field produced by the Mw 6.1 2014 South Napa Valley earthquake and then we used the results from the integration to perform the source modeling. The first Chapter is meant to introduce the topic of the progressive use of Remote Sensing geodetic data to support the activities of monitoring and hazard mitigation related to natural phenomena. Chapter 2 shows the application of the InSAR technique to reconstruct and model surface displacement fields induced by several phenomena. In Chapter 3, the 3D coseismic displacement map due to the 2014 Mw 6.1 South Napa earthquake, close the San Andreas Fault system (California), is estimated by using a method to merge InSAR and GPS data. InSAR data are provided by the latest satellite of the European Space Agency (ESA), i.e. Sentinel-1, whereas the GPS data were obtained from the BARD network and several online archives. In Chapter 4 we propose an improved algorithm for the data integration and test it on the Napa earthquake. Geodetic data from MAI and POT are added in the processing chain and the GPS data interpolation is modified according to the specific phenomenon. Futhermore, the source modeling is performed by inversion of the obtained 3D displacement component. The best fit is obtained by simulating a fracture in the fault segment in agreement with previous works. Finally, in the last chapter we discuss about the advantages and disadvantages of the data integration and the future perspectives

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

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

Detection and Time Series Analysis of Natural Hazards Through Utilization of Optical and Radar Satellite Remote Sensing Techniques

Advancements in satellite remote sensing techniques have revolutionized Earth science observation and are powerful tools used to detect, monitor, and analyze natural hazards. High resolution Digital Surface Models (DSMs) depict Earth&rsquo;s topographic surface including the natural environment, like trees, and artificial features such as buildings, at meter-submeter spatial resolution. When DSMs are acquired at various times overlapping the same region, they can be differenced to calculate areas or volumes of ground elevation change. First, I use high resolution DSMs to characterize the 2015 Taan Fiord tsunamigenic landslide in Alaska, which generated a displacement wave with a 193-m runup. I create six, 2-meter posting DSMs using DigitalGlobe/Maxar satellite imagery acquired near-annually between 2012-2019, and the Surface Extraction with TIN-based Search-space Minimization (SETSM) high-performance computing algorithm. By relatively aligning each acquisition and differencing them through time, I find that the landslide mobilized roughly 77. 0 &plusmn; 0.9 Mm3 of material, of which approximately 56.3 Mm2 were deposited in the fjord waters. Furthermore, I quantify an additional 27.2 &plusmn;3.8 Mm3 of material scoured from fjord-adjacent hillslopes and deposited in the fjord waters, providing new constraints on the subaqueous deposition. This is the first time that DSMs have been used to estimate the volume of scour caused by a tsunami and the subsequent changes in extent and volume with time. I also identify precursory motion prior to the 2015 landslide, characterize several smaller-scale landslides in the larger Taan Fiord region, delineate terminus positions and associated ice dynamics of Tyndall Glacier, and detail seasonal changes in vegetation growth and snow melt/accumulation Another remote sensing technique, Differential Interferometric Synthetic Aperture Radar (DInSAR), quantifies line-of-sight (LOS) surface deformation with cm-mm precision. Sentinel-1A/B satellites acquire images over regions 100km2 wide with an orbital repeat cycle of 6-12 days, collecting large datasets with dense spatio-temporal coverage. Global Navigation Satellite System (GNSS) and DInSAR may be combined to improve the accuracy of deformation results. Next, I integrate DInSAR with GNSS time series to create a fused dataset with enhanced accuracy of 3D ground motions from volcanic eruptions in Hawaii from November 2015 to April 2021. I present a comparison of the raw datasets against the fused time series and give a detailed account of observed ground deformation leading to the May 2018 and December 2020 volcanic eruptions. A DSM is required as input for SAR processing and plays an important role in isolating the deformation signal of interest by removing topographic noise contained in the data. In my final chapter, I compare Sentinel-1A/B DInSAR time series results over Hawaii using multiple elevation products sampled at various spatial resolutions to determine the optimal topographic correction. I construct WorldView imagery at half-meter spatial resolution and test new methods to correct image distortions, co-register vertical accuracy with ICESat-2 geolocated photons, and mosaic images together to achieve large spatial coverage. I find that high vertical accuracy, high spatial resolution, and the use of updated topography have a significant impact on time series results. During periods of stable ground conditions, WorldView DSMs can improve spatial and temporal coherence between SAR images by 10.3-cm.</p

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

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

Correcting atmospheric effects on InSAR with MERIS water vapour data and elevation-dependent interpolation model

Author name used in this publication: X. L. Ding2011-2012 > Academic research: refereed > Publication in refereed journalVersion of RecordPublishe

Modelling co- and post-seismic displacements revealed by InSAR, and their implications for fault behaviour

The ultimate goal of seismology is to estimate the timing, magnitude and potential spatial extent of future seismic events along pre-existing faults. Based on the rate-state friction law, several theoretical physical earthquake models have been proposed towards this goal. Tectonic loading rate and frictional properties of faults are required in these models. Modern geodetic observations, e.g. GPS and InSAR, have provided unprecedented near-field observations following large earthquakes. In theory, according to the frictional rate and state asperity earthquake model, velocity-weakening regions holding seismic motions on faults should be separated with velocity-strengthening regions within which faults slip only aseismically. However, early afterslip following the 2011 MW 9.1 Tohoku-Oki earthquake revealed from GPS measurements was largely overlaid on the historical rupture zones, which challenged the velocity weakening asperity model. Therefore, the performance of the laboratory based friction law in the natural events needs further investigation, and the factors that may affect the estimates of slip models through geodetic modelling should also be discussed systematically. In this thesis, several moderate-strong events were investigated in order to address this important issue. The best-fit co- and post-seismic slip models following the 2009 MW 6.3 Haixi, Qinghai thrust-slip earthquake determined by InSAR deformation time-series suggest that the maximum afterslip is concentrated in the same area as the coseismic slip model, which is similar to the patterns observed in the 2011 Japan earthquake. In this case, complex geometric asperity may play a vital role in the coseismic nucleation and postseismic faulting. The major early afterslip after the 2011 MW 7.1 Van mainshock, which was revealed by one COSMO-SkyMed postseismic interferogram, is found just above the coseismic slip pattern. In this event, a postseismic modelling that did not allow slip across the coseismic asperity was also tested, suggesting that the slip model without slip in the asperities can explain the postseismic observations as well as the afterslip model without constraints on slip in the asperities. In the 2011 MW 9.1 Tohoku-Oki earthquake, a joint inversion with the GRACE coseismic gravity changes and inland coseismic GPS observations was conducted to re-investigate the coseismic slip model of the mainshock. A comparison of slip models from these different datasets suggests that significant variations of slip models can be observed, particularly the locations of the maximum slips. The joint slip model shows that the maximum slip of ~42 m appears near the seafloor surface close to the Japan Trench. Meanwhile, the accumulative afterslip patterns (slip >2 m) determined in previous studies appear in spatial correlation with the Coulomb stress changes generated using the joint slip model. As a strike-slip faulting event, the 2011 MW 6.8 Yushu earthquake was also investigated through co- and post-seismic modelling with more SAR data than was used in previous study. Best slip models suggest that the major afterslip is concentrated in shallow parts of the faults and between the two major coseismic slip patterns, suggesting that the performance of the rate and state frictional asperity model is appropriate in this event. Other postseismic physical mechanisms, pore-elastic rebound and viscoelastic relaxation have also been examined, which cannot significantly affect the estimate of the shallow afterslip model in this study. It is believed that the shallow afterslip predominantly controlled the postseismic behaviour after the mainshock in this case. In comparison to another 21 earthquakes investigated using geodetic data from other studies, complementary spatial extents between co- and post-seismic slip models can be identified. The 2009 MW 6.3 Qinghai earthquake is an exceptional case, in which the faulting behaviours might be dominated by the fault structure (e.g. fault bending). In conclusion, the major contributions from this thesis include: 1) the friction law gives a first order fit in most of natural events examined in this thesis; 2) geometric asperities may play an important role in faulting during earthquake cycles; 3) significant uncertainties in co- and post-seismic slip models can appreciably bias the estimation of fault frictional properties; 4) new insights derived from each earthquake regarding their fault structures and complex faulting behaviours have been observed in this thesis; and (5) a novel package for geodetic earthquake modelling has been developed, which can handle multiple datasets including InSAR, GPS and land/space based gravity changes

Volcanic Processes Monitoring and Hazard Assessment Using Integration of Remote Sensing and Ground-Based Techniques

The monitoring of active volcanoes is a complex task based on multidisciplinary and integrated analyses that use ground, drones and satellite monitoring devices. Over time, and with the development of new technologies and increasing frequency of acquisition, the use of remote sensing to accomplish this important task has grown enormously. This is especially so with the use of drones and satellites for classifying eruptive events and detecting the opening of new vents, the spreading of lava flows on the surface or ash plumes in the atmosphere, the fallout of tephra on the ground, the intrusion of new magma within the volcano edifice, and the deformation preceding impending eruptions, and many other factors. The main challenge in using remote sensing techniques is to develop automated and reliable systems that may assist the decision maker in volcano monitoring, hazard assessment and risk reduction. The integration with ground-based techniques represents a valuable additional aspect that makes the proposed methods more robust and reinforces the results obtained. This collection of papers is focused on several active volcanoes, such as Stromboli, Etna, and Volcano in Italy; the Long Valley caldera and Kilauea volcano in the USA; and Cotopaxi in Ecuador