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

Ground-based synthetic aperture radar (GBSAR) interferometry for deformation monitoring

Ph. D ThesisGround-based synthetic aperture radar (GBSAR), together with interferometry, represents a powerful tool for deformation monitoring. GBSAR has inherent flexibility, allowing data to be collected with adjustable temporal resolutions through either continuous or discontinuous mode. The goal of this research is to develop a framework to effectively utilise GBSAR for deformation monitoring in both modes, with the emphasis on accuracy, robustness, and real-time capability. To achieve this goal, advanced Interferometric SAR (InSAR) processing algorithms have been proposed to address existing issues in conventional interferometry for GBSAR deformation monitoring. The proposed interferometric algorithms include a new non-local method for the accurate estimation of coherence and interferometric phase, a new approach to selecting coherent pixels with the aim of maximising the density of selected pixels and optimizing the reliability of time series analysis, and a rigorous model for the correction of atmospheric and repositioning errors. On the basis of these algorithms, two complete interferometric processing chains have been developed: one for continuous and the other for discontinuous GBSAR deformation monitoring. The continuous chain is able to process infinite incoming images in real time and extract the evolution of surface movements through temporally coherent pixels. The discontinuous chain integrates additional automatic coregistration of images and correction of repositioning errors between different campaigns. Successful deformation monitoring applications have been completed, including three continuous (a dune, a bridge, and a coastal cliff) and one discontinuous (a hillside), which have demonstrated the feasibility and effectiveness of the presented algorithms and chains for high-accuracy GBSAR interferometric measurement. Significant deformation signals were detected from the three continuous applications and no deformation from the discontinuous. The achieved results are justified quantitatively via a defined precision indicator for the time series estimation and validated qualitatively via a priori knowledge of these observing sites.China Scholarship Council (CSC), Newcastle Universit

Minimizing the residual topography effect on interferograms to improve DInSAR results: estimating land subsidence in Port-Said City, Egypt

The accurate detection of land subsidence rates in urban areas is important to identify damage-prone areas and provide decision-makers with useful information. Meanwhile, no precise measurements of land subsidence have been undertaken within the coastal Port-Said City in Egypt to evaluate its hazard in relationship to sea-level rise. In order to address this shortcoming, this work introduces and evaluates a methodology that substantially improves small subsidence rate estimations in an urban setting. Eight ALOS/PALSAR-1 scenes were used to estimate the land subsidence rates in Port-Said City, using the Small BAse line Subset (SBAS) DInSAR technique. A stereo pair of ALOS/PRISM was used to generate an accurate DEM to minimize the residual topography effect on the generated interferograms. A total of 347 well distributed ground control points (GCP) were collected in Port-Said City using the leveling instrument to calibrate the generated DEM. Moreover, the eight PALSAR scenes were co-registered using 50 well-distributed GCPs and used to generate 22 interferogram pairs. These PALSAR interferograms were subsequently filtered and used together with the coherence data to calculate the phase unwrapping. The phase-unwrapped interferogram-pairs were then evaluated to discard four interferograms that were affected by phase jumps and phase ramps. Results confirmed that using an accurate DEM (ALOS/PRISM) was essential for accurately detecting small deformations. The vertical displacement rate during the investigated period (2007–2010) was estimated to be −28 mm. The results further indicate that the northern area of Port-Said City has been subjected to higher land subsidence rates compared to the southern area. Such land subsidence rates might induce significant environmental changes with respect to sea-level rise

InSAR Modeling of Geophysics Measurements

In the present work, the geometry and basic parameters of interferometric synthetic aperture radar (InSAR) geophysics system are addressed. Equations of pixel height and displacement evaluation are derived. Synthetic aperture radar (SAR) signal model based on linear frequency modulation (LFM) waveform and image reconstruction procedure are suggested. The concept of pseudo InSAR measurements, interferogram, and differential interferogram generation is considered. Interferogram and differential interferogram are generated based on a surface model and InSAR measurements. Results of numerical experiments are provided

Remote Monitoring of Civil Infrastructure Based on TomoSAR

Structural health monitoring and damage detection tools are extremely important topics nowadays with the civil infrastructure aging and deteriorating problems observed in urban areas. These tasks can be done by visual inspection and by using traditional in situ methods, such as leveling or using traditional mechanical and electrical sensors, but these approaches are costly, labor-intensive and cannot be performed with a high temporal frequency. In recent years, remote sensing has proved to be a very promising methodology in evaluating the health of a structure by assessing its deformation and thermal dilation. The satellite-based Synthetic Aperture Radar Tomography (TomoSAR) technique, based on the exploitation of a stack of multi-temporal SAR images, allows to remotely sense the movement and the thermal dilation of individual structures with a centimeter-to millimeter-level accuracy, thanks to new generation high-resolution satellite-borne sensors. In this paper, the effectiveness of a recently developed TomoSAR technique in assessing both possible deformations and the thermal dilation evolution of man-made structures is shown. The results obtained using X-band SAR data in two case studies, concerning two urban structures in the city of Naples (Italy), are presented

Post-failure evolution analysis of a rainfall-triggered landslide by multi-temporal interferometry SAR approaches integrated with geotechnical analysis

Persistent Scatterers Interferometry (PSI) represents one of the most powerful techniques for Earth's surface deformation processes' monitoring, especially for long-term evolution phenomena. In this work, a dataset of 34 TerraSAR-X StripMap images (October 2013–October 2014) has been processed by two PSI techniques - Coherent Pixel Technique-Temporal Sublook Coherence (CPT-TSC) and Small Baseline Subset (SBAS) - in order to study the evolution of a slow-moving landslide which occurred on February 23, 2012 in the Papanice hamlet (Crotone municipality, southern Italy) and induced by a significant rainfall event (185 mm in three days). The mass movement caused structural damage (buildings' collapse), and destruction of utility lines (gas, water and electricity) and roads. The results showed analogous displacement rates (30–40 mm/yr along the Line of Sight – LOS-of the satellite) with respect to the pre-failure phase (2008–2010) analyzed in previous works. Both approaches allowed detect the landslide-affected area, however the higher density of targets identified by means of CPT-TSC enabled to analyze in detail the slope behavior in order to design possible mitigation interventions. For this aim, a slope stability analysis has been carried out, considering the comparison between groundwater oscillations and time-series of displacement. Hence, the crucial role of the interaction between rainfall and groundwater level has been inferred for the landslide triggering. In conclusion, we showed that the integration of geotechnical and remote sensing approaches can be seen as the best practice to support stakeholders to design remedial works.Peer ReviewedPostprint (author's final draft

Feasibility of simultaneous intracranial EEG-fMRI in humans: a safety study

In epilepsy patients who have electrodes implanted in their brains as part of their pre-surgical assessment, simultaneous intracranial EEG and fMRI (icEEG-fMRI) may provide important localising information and improve understanding of the underlying neuropathology. However, patient safety during icEEG-fMRI has not been addressed. Here the potential health hazards associated with icEEG-fMRI were evaluated theoretically and the main risks identified as: mechanical forces on electrodes from transient magnetic effects, tissue heating due to interaction with the pulsed RF fields and tissue stimulation due to interactions with the switched magnetic gradient fields. These potential hazards were examined experimentally in vitro on a Siemens 3 T Trio, 1.5 T Avanto and a GE 3 T Signa Excite scanner using a Brain Products MR compatible EEG system. No electrode flexion was observed. Temperature measurements demonstrated that heating well above guideline limits can occur. However heating could be kept within safe limits (< 1.0 °C) by using a head transmit RF coil, ensuring EEG cable placement to exit the RF coil along its central z-axis, using specific EEG cable lengths and limiting MRI sequence specific absorption rates (SARs). We found that the risk of tissue damage due to RF-induced heating is low provided implant and scanner specific SAR limits are observed with a safety margin used to account for uncertainties (e.g. in scanner-reported SAR). The observed scanner gradient switching induced current (0.08 mA) and charge density (0.2 μC/cm2) were well within safety limits (0.5 mA and 30 μC/cm2, respectively). Site-specific testing and a conservative approach to safety are required to avoid the risk of adverse events

Characterizing slope instability kinematics by integrating multi-sensor satellite remote sensing observations

Over the past few decades, the occurrence and intensity of geological hazards, such as landslides, have substantially risen due to various factors, including global climate change, seismic events, rapid urbanization and other anthropogenic activities. Landslide disasters pose a significant risk in both urban and rural areas, resulting in fatalities, infrastructure damages, and economic losses. Nevertheless, conventional ground-based monitoring techniques are often costly, time-consuming, and require considerable resources. Moreover, some landslide incidents occur in remote or hazardous locations, making ground-based observation and field investigation challenging or even impossible. Fortunately, the advancements in spaceborne remote sensing technology have led to the availability of large-scale and high-quality imagery, which can be utilized for various landslide-related applications, including identification, monitoring, analysis, and prediction. This efficient and cost-effective technology allows for remote monitoring and assessment of landslide risks and can significantly contribute to disaster management and mitigation efforts. Consequently, spaceborne remote sensing techniques have become vital for geohazard management in many countries, benefiting society by providing reliable downstream services. However, substantial effort is required to ensure that such benefits are provided. For establishing long-term data archives and reliable analyses, it is essential to maintain consistent and continued use of multi-sensor spaceborne remote sensing techniques. This will enable a more thorough understanding of the physical mechanisms responsible for slope instabilities, leading to better decision-making and development of effective mitigation strategies. Ultimately, this can reduce the impact of landslide hazards on the general public. The present dissertation contributes to this effort from the following perspectives: 1. To obtain a comprehensive understanding of spaceborne remote sensing techniques for landslide monitoring, we integrated multi-sensor methods to monitor the entire life cycle of landslide dynamics. We aimed to comprehend the landslide evolution under complex cascading events by utilizing various spaceborne remote sensing techniques, e.g., the precursory deformation before catastrophic failure, co-failure procedures, and post-failure evolution of slope instability. 2. To address the discrepancies between spaceborne optical and radar imagery, we present a methodology that models four-dimensional (4D) post-failure landslide kinematics using a decaying mathematical model. This approach enables us to represent the stress relaxation for the landslide body dynamics after failure. By employing this methodology, we can overcome the weaknesses of the individual sensor in spaceborne optical and radar imaging. 3. We assessed the effectiveness of a newly designed small dihedral corner reflector for landslide monitoring. The reflector is compatible with both ascending and descending satellite orbits, while it is also suitable for applications with both high-resolution and medium-resolution satellite imagery. Furthermore, although its echoes are not as strong as those of conventional reflectors, the cost of the newly designed reflectors is reduced, with more manageable installation and maintenance. To overcome this limitation, we propose a specific selection strategy based on a probability model to identify the reflectors in satellite images