Generic interferometric synthetic aperture radar atmospheric correction model and its application to co- and post-seismic motions

PhD ThesisThe tremendous development of Interferometric Synthetic Aperture Radar (InSAR) missions in recent years facilitates the study of smaller amplitude ground deformation over greater spatial scales using longer time series. However, this poses more challenges for correcting atmospheric effects due to the spatial-temporal variability of atmospheric delays. Previous attempts have used observations from Global Positioning System (GPS) and Numerical Weather Models (NWMs) to separate the atmospheric delays, but they are limited by (i) the availability (and distribution) of GPS stations; (ii) the time difference between NWM and radar observations; and (iii) the difficulties in quantifying their performance. To overcome the abovementioned limitations, we have developed the Iterative Tropospheric Decomposition (ITD) model to reduce the coupling effects of the troposphere turbulence and stratification and hence achieve similar performances over flat and mountainous terrains. Highresolution European Centre for Medium-Range Weather Forecasts (ECMWF) and GPS-derived tropospheric delays were properly integrated by investigating the GPS network geometry and topography variations. These led to a generic atmospheric correction model with a range of notable features: (i) global coverage, (ii) all-weather, all-time usability, (iii) available with a maximum of two-day latency, and (iv) indicators available to assess the model’s performance and feasibility. The generic atmospheric correction model enables the investigation of the small magnitude coseismic deformation of the 2017 Mw-6.4 Nyingchi earthquake from InSAR observations in spite of substantial atmospheric contamination. It can also minimize the temporal correlations of InSAR atmospheric delays so that reliable velocity maps over large spatial extents can be achieved. Its application to the post-seismic motion following the 2016 Kaikoura earthquake shows a success to recover the time-dependent afterslip distribution, which in turn evidences the deep inactive subduction slip mechanism. This procedure can be used to map surface deformation in other scenarios including volcanic eruptions, tectonic rifting, cracking, and city subsidence.This work was supported by a Chinese Scholarship Council studentship. Part of this work was also supported by the UK NERC through the Centre for the Observation and Modelling of Earthquakes, Volcanoes and Tectonics (COMET)

Applications of SAR Interferometry in Earth and Environmental Science Research

This paper provides a review of the progress in regard to the InSAR remote sensing technique and its applications in earth and environmental sciences, especially in the past decade. Basic principles, factors, limits, InSAR sensors, available software packages for the generation of InSAR interferograms were summarized to support future applications. Emphasis was placed on the applications of InSAR in seismology, volcanology, land subsidence/uplift, landslide, glaciology, hydrology, and forestry sciences. It ends with a discussion of future research directions

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

Fusion of remotely sensed displacement measurements: current status and challenges

International audienc

A Spatially Varying Scaling Method for InSAR Tropospheric Corrections Using a High‐Resolution Weather Model

Variation in tropospheric delay is a major limiting factor on the accuracy of interferometric synthetic aperture radar (InSAR) measurements. This is particularly the case when deformation and topography are correlated. To address limitations of previous InSAR tropospheric correction methods, here we present a new approach that combines the use of both external weather model data and the interferometric phase. We assume that vertical refractivity profiles calculated from a high‐resolution weather model data can generally describe the form of the relationship between tropospheric delay and height but that the magnitude can be incorrect. We estimate a magnitude correction by scaling the original delays to best match the interferometric phase. We validated our new method using simulated data and demonstrate that both coseismic and interseismic signals can be separated from strong tropospheric delays. We also applied our algorithm to the central portion of the Altyn Tagh Fault in northern Tibet, where deformation correlates strongly with topographic relief of 6,000 m, and show that the derived velocity field is more internally consistent and agrees better with independent Global Positioning System measurements

Digital image correlation (DIC) analysis of the 3 December 2013 Montescaglioso landslide (Basilicata, Southern Italy). Results from a multi-dataset investigation

Image correlation remote sensing monitoring techniques are becoming key tools for providing effective qualitative and quantitative information suitable for natural hazard assessments, specifically for landslide investigation and monitoring. In recent years, these techniques have been successfully integrated and shown to be complementary and competitive with more standard remote sensing techniques, such as satellite or terrestrial Synthetic Aperture Radar interferometry. The objective of this article is to apply the proposed in-depth calibration and validation analysis, referred to as the Digital Image Correlation technique, to measure landslide displacement. The availability of a multi-dataset for the 3 December 2013 Montescaglioso landslide, characterized by different types of imagery, such as LANDSAT 8 OLI (Operational Land Imager) and TIRS (Thermal Infrared Sensor), high-resolution airborne optical orthophotos, Digital Terrain Models and COSMO-SkyMed Synthetic Aperture Radar, allows for the retrieval of the actual landslide displacement field at values ranging from a few meters (2–3 m in the north-eastern sector of the landslide) to 20–21 m (local peaks on the central body of the landslide). Furthermore, comprehensive sensitivity analyses and statistics-based processing approaches are used to identify the role of the background noise that affects the whole dataset. This noise has a directly proportional relationship to the different geometric and temporal resolutions of the processed imagery. Moreover, the accuracy of the environmental-instrumental background noise evaluation allowed the actual displacement measurements to be correctly calibrated and validated, thereby leading to a better definition of the threshold values of the maximum Digital Image Correlation sub-pixel accuracy and reliability (ranging from 1/10 to 8/10 pixel) for each processed dataset

Coseismic deformation observed with radar interferometry: Great earthquakes and atmospheric noise

Spatially dense maps of coseismic deformation derived from Interferometric Synthetic Aperture Radar (InSAR) datasets result in valuable constraints on earthquake processes. The recent increase in the quantity of observations of coseismic deformation facilitates the examination of signals in many tectonic environments associated with earthquakes of varying magnitude. Efforts to place robust constraints on the evolution of the crustal stress field following great earthquakes often rely on knowledge of the earthquake location, the fault geometry, and the distribution of slip along the fault plane. Well-characterized uncertainties and biases strengthen the quality of inferred earthquake source parameters, particularly when the associated ground displacement signals are near the detection limit. Well-preserved geomorphic records of earthquakes offer additional insight into the mechanical behavior of the shallow crust and the kinematics of plate boundary systems. Together, geodetic and geologic observations of crustal deformation offer insight into the processes that drive seismic cycle deformation over a range of timescales. In this thesis, I examine several challenges associated with the inversion of earthquake source parameters from SAR data. Variations in atmospheric humidity, temperature, and pressure at the timing of SAR acquisitions result in spatially correlated phase delays that are challenging to distinguish from signals of real ground deformation. I characterize the impact of atmospheric noise on inferred earthquake source parameters following elevation-dependent atmospheric corrections. I analyze the spatial and temporal variations in the statistics of atmospheric noise from both reanalysis weather models and InSAR data itself. Using statistics that reflect the spatial heterogeneity of atmospheric characteristics, I examine parameter errors for several synthetic cases of fault slip on a basin-bounding normal fault. I show a decrease in uncertainty in fault geometry and kinematics following the application of atmospheric corrections to an event spanned by real InSAR data, the 1992 M5.6 Little Skull Mountain, Nevada, earthquake. Finally, I discuss how the derived workflow could be applied to other tectonic problems, such as solving for interseismic strain accumulation rates in a subduction zone environment. I also study the evolution of the crustal stress field in the South American plate following two recent great earthquakes along the Nazca- South America subduction zone. I show that the 2010 Mw 8.8 Maule, Chile, earthquake very likely triggered several moderate magnitude earthquakes in the Andean volcanic arc and backarc. This suggests that great earthquakes modulate the crustal stress field outside of the immediate aftershock zone and that far-field faults may pose a heightened hazard following large subduction earthquakes. The 2014 Mw 8.1 Pisagua, Chile, earthquake reopened ancient surface cracks that have been preserved in the hyperarid forearc setting of northern Chile for thousands of earthquake cycles. The orientation of cracks reopened in this event reflects the static and likely dynamic stresses generated by the recent earthquake. Coseismic cracks serve as a reliable marker of permanent earthquake deformation and plate boundary behavior persistent over the million-year timescale. This work on great earthquakes suggests that InSAR observations can play a crucial role in furthering our understanding of the crustal mechanics that drive seismic cycle processes in subduction zones

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

Investigating the postseismic deformation of strike-slip earthquakes on the Tibetan Plateau

InSAR is a useful technique to detect large-scale surface deformation from space. To place constraints on the rheological structure of the lithosphere in the Tibetan Plateau, two strike-slip earthquakes have been investigated. One is the Mw 7.6 Manyi earthquake, which occurred in the north-central Tibetan Plateau. The other is the Mw 6.5 Jiuzhaigou earthquake, which happened that on the eastern part of the Tibetan Plateau. My InSAR data cover 12 years following the Manyi earthquake, much longer than previous researchers’ dataset. I test three viscoelastic models (Maxwell, Standard linear solids, and Burgers body) and one afterslip model. The viscoelastic models cannot match the observed temporal–spatial deformation patterns. The distributions of deformation in the viscoelastic models extend into the far field and the residuals tend to increase, which are inconsistent with the data. The afterslip model has the lowest misfit and explains the temporal and spatial pattern of the observed deformation with decent result. A combined model that considers the effects of both afterslip and viscoelastic relaxation has also been tested. In this combined model, the viscoelastic relaxation that occurs with an elastic layer of thickness of 30 km over a half-space place, produce an estimate for viscosity of 5 × 1019 Pa s for this area. Therefore, either the afterslip model or the combined model can be used to explain the 12 years postseismic deformation of Manyi earthquake. The long time series of the Manyi earthquake enable us to distinguish between afterslip and viscoelastic relaxation. The seismogenic fault of the Jiuzhaigou earthquake was previously unidentified and no surface rupture is found after the earthquake. I first determined the fault geometry and calculated coseismic slip model. The slip model indicates a left-lateral strike-slip pattern, which is consistent with focal mechanisms were determined by different agencies. There is no visible postseismic deformation signal of the fault, which means the surface deformation generated by fault creeping is smaller than the noise of our observation method over that period. Therefore, I try to find the lower bound of the viscosity for this area. My preferred minimum possible viscosity of the underlying half-space is ∼6 × 1017 Pa s. Together with previous geodetic studies, the viscosities obtained from central Tibet show at least one order of magnitude difference with the viscosities obtained from the eastern Tibet. The heterogeneity indicates the rheology has a relatively large spatial change through the whole Plateau. The viscoelastic model always been proposed to explain long-term postseismic deformation and afterslip is used to explain the short-term deformation or localised deformation. Sometimes, the viscoelastic deformation signal is invisible in the moderate earthquakes as the stress is not large enough to generate observable deformation