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

Investigating SAR algorithm for spaceborne interferometric oil spill detection

The environmental damages and recovery of terrestrial ecosystems from oil spills can last decades. Oil spills have been responsible for loss of aquamarine lives, organisms, trees, vegetation, birds and wildlife. Although there are several methods through which oil spills can be detected, it can be argued that remote sensing via the use of spaceborne platforms provides enormous benefits. This paper will provide more efficient means and methods that can assist in improving oil spill responses. The objective of this research is to develop a signal processing algorithm that can be used for detecting oil spills using spaceborne SAR interferometry (InSAR) data. To this end, a pendulum formation of multistatic smallSAR carrying platforms in a near equatorial orbit is described. The characteristic parameters such as the effects of incidence angles on radar backscatter, which support the detection of oil spills, will be the main drivers for determining the relative positions of the small satellites in formation. The orbit design and baseline distances between each spaceborne SAR platform will also be discussed. Furthermore, results from previous analysis on coverage assessment and revisit time shall be highlighted. Finally, an evaluation of automatic algorithm techniques for oil spill detection in SAR images will be conducted and results presented. The framework for the automatic algorithm considered consists of three major steps. The segmentation stage, where techniques that suggest the use of thresholding for dark spot segmentation within the captured InSAR image scene is conducted. The feature extraction stage involves the geometry and shape of the segmented region where elongation of the oil slick is considered an important feature and a function of the width and the length of the oil slick. For the classification stage, where the major objective is to distinguish oil spills from look-alikes, a Mahalanobis classifier will be used to estimate the probability of the extracted features being oil spills. The validation process of the algorithm will be conducted by using NASA’s UAVSAR data obtained over the Gulf of coast oil spill and RADARSAT-1 dat

Development of the TanDEM-X Calibration Concept: Analysis of Systematic Errors

The TanDEM-X mission, result of the partnership between the German Aerospace Center (DLR) and Astrium GmbH, opens a new era in spaceborne radar remote sensing. The first bistatic satellite synthetic aperture radar mission is formed by flying the TanDEM-X and TerraSAR-X in a closely controlled helix formation. The primary mission goal is the derivation of a high-precision global digital elevation model (DEM) according to High-Resolution Terrain Information (HRTI) level 3 accuracy. The finite precision of the baseline knowledge and uncompensated radar instrument drifts introduce errors that may compromise the height accuracy requirements. By means of a DEM calibration, which uses absolute height references, and the information provided by adjacent interferogram overlaps, these height errors can be minimized. This paper summarizes the exhaustive studies of the nature of the residual-error sources that have been carried out during the development of the DEM calibration concept. Models for these errors are set up and simulations of the resulting DEM height error for different scenarios provide the basis for the development of a successful DEM calibration strategy for the TanDEM-X mission

Conceptual Study and Performance Analysis of Tandem Dual-Antenna Spaceborne SAR Interferometry

Multi-baseline synthetic aperture radar interferometry (MB-InSAR), capable of mapping 3D surface model with high precision, is able to overcome the ill-posed problem in the single-baseline InSAR by use of the baseline diversity. Single pass MB acquisition with the advantages of high coherence and simple phase components has a more practical capability in 3D reconstruction than conventional repeat-pass MB acquisition. Using an asymptotic 3D phase unwrapping (PU), it is possible to get a reliable 3D reconstruction using very sparse acquisitions but the interferograms should follow the optimal baseline design. However, current spaceborne SAR system doesn't satisfy this principle, inducing more difficulties in practical application. In this article, a new concept of Tandem Dual-Antenna SAR Interferometry (TDA-InSAR) system for single-pass reliable 3D surface mapping using the asymptotic 3D PU is proposed. Its optimal MB acquisition is analyzed to achieve both good relative height precision and flexible baseline design. Two indicators, i.e., expected relative height precision and successful phase unwrapping rate, are selected to optimize the system parameters and evaluate the performance of various baseline configurations. Additionally, simulation-based demonstrations are conducted to evaluate the performance in typical scenarios and investigate the impact of various error sources. The results indicate that the proposed TDA-InSAR is able to get the specified MB acquisition for the asymptotic 3D PU, which offers a feasible solution for single-pass 3D SAR imaging.Comment: 16 pages, 20 figure

Reliable estimation of orbit errors in spaceborne SAR interferometry. The network approach

An approach to improve orbital state vectors by orbit error estimates derived from residual phase patterns in synthetic aperture radar interferograms is presented. For individual interferograms, an error representation by two parameters is motivated: the baseline error in cross-range and the rate of change of the baseline error in range. For their estimation, two alternatives are proposed: a least squares approach that requires prior unwrapping and a less reliable gridsearch method handling the wrapped phase. In both cases, reliability is enhanced by mutual control of error estimates in an overdetermined network of linearly dependent interferometric combinations of images. Thus, systematic biases, e.g., due to unwrapping errors, can be detected and iteratively eliminated. Regularising the solution by a minimum-norm condition results in quasi-absolute orbit errors that refer to particular images. For the 31 images of a sample ENVISAT dataset, orbit corrections with a mutual consistency on the millimetre level have been inferred from 163 interferograms. The method itself qualifies by reliability and rigorous geometric modelling of the orbital error signal but does not consider interfering large scale deformation effects. However, a separation may be feasible in a combined processing with persistent scatterer approaches or by temporal filtering of the estimates

Effects of Orbit and Pointing Geometry of a Spaceborne Formation for Monostatic-Bistatic Radargrammetry on Terrain Elevation Measurement Accuracy

During the last decade a methodology for the reconstruction of surface relief by Synthetic Aperture Radar (SAR) measurements – SAR interferometry – has become a standard. Different techniques developed before, such as stereo-radargrammetry, have been experienced from space only in very limiting geometries and time series, and, hence, branded as less accurate. However, novel formation flying configurations achievable by modern spacecraft allow fulfillment of SAR missions able to produce pairs of monostatic-bistatic images gathered simultaneously, with programmed looking angles. Hence it is possible to achieve large antenna separations, adequate for exploiting to the utmost the stereoscopic effect, and to make negligible time decorrelation, a strong liming factor for repeat-pass stereo-radargrammetric techniques. This paper reports on design of a monostatic-bistatic mission, in terms of orbit and pointing geometry, and taking into account present generation SAR and technology for accurate relative navigation. Performances of different methods for monostatic-bistatic stereo-radargrammetry are then evaluated, showing the possibility to determine the local surface relief with a metric accuracy over a wide range of Earth latitudes

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

Elevation and Deformation Extraction from TomoSAR

3D SAR tomography (TomoSAR) and 4D SAR differential tomography (Diff-TomoSAR) exploit multi-baseline SAR data stacks to provide an essential innovation of SAR Interferometry for many applications, sensing complex scenes with multiple scatterers mapped into the same SAR pixel cell. However, these are still influenced by DEM uncertainty, temporal decorrelation, orbital, tropospheric and ionospheric phase distortion and height blurring. In this thesis, these techniques are explored. As part of this exploration, the systematic procedures for DEM generation, DEM quality assessment, DEM quality improvement and DEM applications are first studied. Besides, this thesis focuses on the whole cycle of systematic methods for 3D & 4D TomoSAR imaging for height and deformation retrieval, from the problem formation phase, through the development of methods to testing on real SAR data. After DEM generation introduction from spaceborne bistatic InSAR (TanDEM-X) and airborne photogrammetry (Bluesky), a new DEM co-registration method with line feature validation (river network line, ridgeline, valley line, crater boundary feature and so on) is developed and demonstrated to assist the study of a wide area DEM data quality. This DEM co-registration method aligns two DEMs irrespective of the linear distortion model, which improves the quality of DEM vertical comparison accuracy significantly and is suitable and helpful for DEM quality assessment. A systematic TomoSAR algorithm and method have been established, tested, analysed and demonstrated for various applications (urban buildings, bridges, dams) to achieve better 3D & 4D tomographic SAR imaging results. These include applying Cosmo-Skymed X band single-polarisation data over the Zipingpu dam, Dujiangyan, Sichuan, China, to map topography; and using ALOS L band data in the San Francisco Bay region to map urban building and bridge. A new ionospheric correction method based on the tile method employing IGS TEC data, a split-spectrum and an ionospheric model via least squares are developed to correct ionospheric distortion to improve the accuracy of 3D & 4D tomographic SAR imaging. Meanwhile, a pixel by pixel orbit baseline estimation method is developed to address the research gaps of baseline estimation for 3D & 4D spaceborne SAR tomography imaging. Moreover, a SAR tomography imaging algorithm and a differential tomography four-dimensional SAR imaging algorithm based on compressive sensing, SAR interferometry phase (InSAR) calibration reference to DEM with DEM error correction, a new phase error calibration and compensation algorithm, based on PS, SVD, PGA, weighted least squares and minimum entropy, are developed to obtain accurate 3D & 4D tomographic SAR imaging results. The new baseline estimation method and consequent TomoSAR processing results showed that an accurate baseline estimation is essential to build up the TomoSAR model. After baseline estimation, phase calibration experiments (via FFT and Capon method) indicate that a phase calibration step is indispensable for TomoSAR imaging, which eventually influences the inversion results. A super-resolution reconstruction CS based study demonstrates X band data with the CS method does not fit for forest reconstruction but works for reconstruction of large civil engineering structures such as dams and urban buildings. Meanwhile, the L band data with FFT, Capon and the CS method are shown to work for the reconstruction of large manmade structures (such as bridges) and urban buildings

Novel closed-loop approaches for precise relative navigation of widely separated GPS receivers in LEO

This paper deals with the relative navigation of a formation of two spacecrafts separated by hundreds of kilometers based on processing dual-frequency differential carrier-phase GPS measurements. Specific requirements of the considered application are high relative positioning accuracy and real-time on board implementation. These can be conflicting requirements. Indeed, if on one hand high accuracy can be achieved by exploiting the integer nature of double-difference carrier-phase ambiguities, on the other hand the presence of large ephemeris errors and differential ionospheric delays makes the integer ambiguities determination challenging. Closed-loop schemes, which update the relative position estimates of a dynamic filter with feedback from integer ambiguities fixing algorithms, are customarily employed in these cases. This paper further elaborates such approaches, proposing novel closed loop techniques aimed at overcoming some of the limitations of traditional algorithms. They extend techniques developed for spaceborne long baseline relative positioning by making use of an on-the-fly ambiguity resolution technique especially developed for the applications of interest. Such techniques blend together ionospheric delay compensation techniques, nonlinear models of relative spacecraft dynamics, and partial integer validation techniques. The approaches are validated using flight data from the Gravity Recovery and Climate Experiment (GRACE) mission. Performance is compared to that of the traditional closed-loop scheme analyzing the capability of each scheme to maximize the percentage of correctly fixed integer ambiguities as well as the relative positioning accuracy. Results show that the proposed approach substantially improves performance of the traditional approaches. More specifically, centimeter-level root-mean square relative positioning is feasible for spacecraft separations of more than 260 km, and an integer ambiguity fixing performance as high as 98% is achieved in a 1-day long dataset. Results also show that approaches exploiting ionospheric delay models are more robust and precise of approaches relying on ionospheric-delay removal techniques. © 2013 IAA