Satellite SAR Interferometry for Earth’s Crust Deformation Monitoring and Geological Phenomena Analysis

Synthetic aperture radar interferometry (InSAR) and the related processing techniques provide a unique tool for the quantitative measurement of the Earth’s surface deformation associated with certain geophysical processes (such as volcanic eruptions, landslides and earthquakes), thus making possible long-term monitoring of surface deformation and analysis of relevant geodynamic phenomena. This chapter provides an application-oriented perspective on the spaceborne InSAR technology with emphasis on subsequent geophysical investigations. First, the fundamentals of radar interferometry and differential interferometry, as well as error sources, are briefly introduced. Emphasis is then placed on the realistic simulation of the underlying geophysics processes, thus offering an unfolded perspective on both analytical and numerical approaches for modeling deformation sources. Finally, various experimental investigations conducted by acquiring SAR multitemporal observations on areas subject to deformation processes of particular geological interest are presented and discussed

Co-seismic deformation during the M_w 7.3 Aqaba earthquake (1995) from ERS-SAR interferometry

The M_w 7.3 1995 Aqaba earthquake is the largest instrumental earthquake along the Dead Sea Fault. We complement previous seismological studies by analyzing co‐seismic ground displacement from differential interferometry computed from ERS images spanning 3 different areas. They are compared with a synthetic model derived from seismological study. Only far‐field deformation related to the main sub‐event could be revealed because the near‐field area lies within the gulf. The interferometric data imply a 56 km long and 10 km wide fault segment, connecting the Elat Deep to the Aragonese Deep, which strikes N195°E and dips 65° to the west, with 2.1 m left‐lateral slip and a 15.5° rake indicating a slight normal component. The geodetic moment compares well with the seismic momen

Surface deformation analysis in Northeast Italy by using PS-InSAR and GNSS data

In the present study, we exploited the potential of satellite-based geodetic data for detecting and measuring surface displacement in Northeast Italy. In this contest, we focused mainly on 1) the estimation of the interseismic deformation during the satellites’ observation period, 2) the detection and analysis of the main deformation patterns, and 3) the correlation of the signals to the active tectonic structures. Despite the low convergence rates (~ 1.5-3 mm/yr), Northeast Italy is an active tectonic area, as testified by the instrumental and historical seismicity. The Adria-Eurasia convergence is mainly accommodated by the thrusts and strike-slip faults of the Southeastern Alps and the External Dinarides, located in the northern and northeastern sectors of the study area. The Venetian-Friulian plain and the Adriatic coasts, affected by active subsidence, dominate the southern region. We used the Stanford Method for Persistent Scatterers (StaMPS) applied to Sentinel-1 SAR images acquired along the ascending and descending orbit tracks between 2015 and 2019. Based on a stack of single-master differential interferograms, we detected coherent and temporally stable pixels based on amplitude and phase noise analysis. After applying spatial-temporal filters and additional post-processing operations to refine the measurements, we used Adria-fixed GNSS velocities derived by permanent stations in the study area to calibrate the InSAR velocities. The outcome consists of Line-OF-Sight (LOS) mean ground velocity maps derived by displacement time series along the radar directions for each satellite track. The combination of the LOS datasets yields vertical and east-west velocity maps, which are mostly in agreement with GNSS data and previous geodetic studies. Based on our measurements, we observe a significant positive velocity gradient of 1 mm/yr across the westernmost sector of the Alpine system, suggesting an aseismic motion of the root of the Bassano-Valdobbiadene thrust. The positive vertical gradients (~1 and up to 2 mm/yr) across the Alpine-Dinaric systems in the central and eastern sectors and the eastward motion that increases northeastward (1-2 mm/yr) may be related to the active Alpine-Dinaric thrusts and strike-slip faults. We also suggest that the detected westward motion of the Friulian plain (around Udine) might be attributed to the presence of tectonic structures characterized by transcurrent-transpressive kinematics. Finally, we detect other signals, such as the significant subsidence (2-4 mm/yr) along the coasts and on the southern Venetian-Friulian plain, confirming the correlation between subsidence and the geological setting of the study area. In conclusion, our study confirms the potential of MT-InSAR and GNSS data for the estimation of the surface deformations in response to active tectonics, even in areas characterized by low deformation rates, such as Northeast Italy.In the present study, we exploited the potential of satellite-based geodetic data for detecting and measuring surface displacement in Northeast Italy. In this contest, we focused mainly on 1) the estimation of the interseismic deformation during the satellites’ observation period, 2) the detection and analysis of the main deformation patterns, and 3) the correlation of the signals to the active tectonic structures. Despite the low convergence rates (~ 1.5-3 mm/yr), Northeast Italy is an active tectonic area, as testified by the instrumental and historical seismicity. The Adria-Eurasia convergence is mainly accommodated by the thrusts and strike-slip faults of the Southeastern Alps and the External Dinarides, located in the northern and northeastern sectors of the study area. The Venetian-Friulian plain and the Adriatic coasts, affected by active subsidence, dominate the southern region. We used the Stanford Method for Persistent Scatterers (StaMPS) applied to Sentinel-1 SAR images acquired along the ascending and descending orbit tracks between 2015 and 2019. Based on a stack of single-master differential interferograms, we detected coherent and temporally stable pixels based on amplitude and phase noise analysis. After applying spatial-temporal filters and additional post-processing operations to refine the measurements, we used Adria-fixed GNSS velocities derived by permanent stations in the study area to calibrate the InSAR velocities. The outcome consists of Line-OF-Sight (LOS) mean ground velocity maps derived by displacement time series along the radar directions for each satellite track. The combination of the LOS datasets yields vertical and east-west velocity maps, which are mostly in agreement with GNSS data and previous geodetic studies. Based on our measurements, we observe a significant positive velocity gradient of 1 mm/yr across the westernmost sector of the Alpine system, suggesting an aseismic motion of the root of the Bassano-Valdobbiadene thrust. The positive vertical gradients (~1 and up to 2 mm/yr) across the Alpine-Dinaric systems in the central and eastern sectors and the eastward motion that increases northeastward (1-2 mm/yr) may be related to the active Alpine-Dinaric thrusts and strike-slip faults. We also suggest that the detected westward motion of the Friulian plain (around Udine) might be attributed to the presence of tectonic structures characterized by transcurrent-transpressive kinematics. Finally, we detect other signals, such as the significant subsidence (2-4 mm/yr) along the coasts and on the southern Venetian-Friulian plain, confirming the correlation between subsidence and the geological setting of the study area. In conclusion, our study confirms the potential of MT-InSAR and GNSS data for the estimation of the surface deformations in response to active tectonics, even in areas characterized by low deformation rates, such as Northeast Italy

The Sentinel-1 mission for the improvement of the scientific understanding and the operational monitoring of the seismic cycle

We describe the state of the art of scientific research on the earthquake cycle based on the analysis of Synthetic Aperture Radar (SAR) data acquired from satellite platforms. We examine the achievements and the main limitations of present SAR systems for the measurement and analysis of crustal deformation, and envision the foreseeable advances that the Sentinel-1 data will generate in the fields of geophysics and tectonics. We also review the technological and scientific issues which have limited so far the operational use of satellite data in seismic hazard assessment and crisis management, and show the improvements expected from Sentinel-1 dat

Interseismic and Postseismic Shallow Creep of the North Qaidam Thrust Faults Detected with a Multitemporal InSAR Analysis

Understanding the mechanisms by which earthquake cycles produce folding and accommodate shortening is essential to quantify the seismic potential of active faults and integrate aseismic slip within our understanding of the physical mechanisms of the long-term deformation. However, measuring such small deformation signals in mountainous areas is challenging with current space-geodesy techniques, due to the low rates of motion relative to the amplitude of the noise. Here we successfully carry out a multitemporal Interferometric Synthetic Aperture Radar analysis over the North Qaidam fold-thrust system in NE Tibet, where eight Mw> 5.2 earthquakes occurred between 2003 and 2009. We report various cases of aseismic slip uplifting the thickened crust at short wavelengths. We provide a rare example of a steep, shallow, 13-km-long and 6-km-wide afterslip signal that coincides spatially with an anticline and that continues into 2011 in response to a Mw 6.3 event in 2003. We suggest that a buried seismic slip during the 2003 earthquake has triggered both plastic an-elastic folding and aseismic slip on the shallow thrusts. We produce a first-order two-dimensional model of the postseismic surface displacements due to the 2003 earthquake and highlight a segmented slip on three fault patches that steepen approaching the surface. This study emphasizes the fundamental role of shallow aseismic slip in the long-term and permanent deformation of thrusts and folds and the potential of Interferometric Synthetic Aperture Radar for detecting and characterizing the spatiotemporal behavior of aseismic slip over large mountainous regions

Seismic Source Quantitative Parameters Retrieval from InSAR Data and Neural Networks

The basic idea of this thesis is to exploit the capabilities of neural networks in a very new framework: the quantitative modelling of the seismic source and the interferogram inversion for retrieving its geometric parameters. The problem can be sum up as follows. When a moderateto- strong earthquake occurs we can apply SAR Interferometry (InSAR) technique to compute a differential interferogram. The latter is used to detect and measure the surface displacement field. The earthquake has been generated by an active, seismogenic, fault having its own specific geometry. Therefore each differential interferogram contains the information concerning the geometry of the seismic source the earthquake comes from; its shape and size, the number of fringes, the lobe orientation and number are the main features of the surface effects field. Two problems have been analysed in this work. The first is the identification of the seismic source mechanism. The second is a typical inversion exercise concerning the fault plane parameter. To perform both exercises of the seismic fault a huge number of synthetic interferograms has been computed. Each of them is generated by a different combination of such geometric parameters. As far as the retrieval of the geometric parameters is concerned an artificial neural network has been properly generated and trained to provide an inversion procedure to single out the geometric parameters of the fault. Five among these latter, Length, Width, Dip, Strike, Depth, have been simultaneously inverted. The result is in agreement with those results based on different approaches. Furthermore the method seems very promising and leads to improve the studies concerning the combined use of neural networks and InSAR technique

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

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