InSAR bias and uncertainty due to the systematic and stochastic tropospheric delay

We quantify the bias and uncertainty of interferometric synthetic aperture radar (InSAR) displacement time series and their derivatives, the displacement velocities, by analyzing the systematic and stochastic components of the temporal variation of the tropospheric delay. The biases due to the systematic seasonal delay depend on the SAR acquisition times, whereas the uncertainties depend on the standard deviation of the random delay, the number of acquisitions, the total time span covered, and the covariance of the time series of the stochastic delay between a pixel and the reference. We study the contribution of the wet delay to the InSAR observations along the western India plate boundary using (i) Moderate Resolution Imaging Spectroradiometer precipitable water vapor, (ii) stratified tropospheric delay estimated from the ERA-I global atmospheric model, and (iii) seven Envisat InSAR swaths. Our analysis indicates that the amplitudes of the annual delay vary by up to ~10 cm in this region equivalent to a maximum displacement bias of ~24 cm in InSAR line of sight direction between two epochs (assuming Envisat IS6 beam mode). The stratified tropospheric delay correction mitigates this bias and reduces the scatter due to the stochastic delay. For ~7 years of Envisat acquisitions along the western India plate boundary, the uncertainty of the InSAR velocity field due to the residual stochastic wet delay after stratified tropospheric delay correction using the ERA-I model is in the order of ~2 mm/yr over 100 km and ~4 mm/yr over 400 km. We discuss the implication of the derived uncertainties on the full variance-covariance matrix of the InSAR data

Ionospheric correction of interferometric SAR data with application to the cryospheric sciences

Thesis (Ph.D.) University of Alaska Fairbanks, 2018The ionosphere has been identified as an important error source for spaceborne Synthetic Aperture Radar (SAR) data and SAR Interferometry (InSAR), especially for low frequency SAR missions, operating, e.g., at L-band or P-band. Developing effective algorithms for the correction of ionospheric effects is still a developing and active topic of remote sensing research. The focus of this thesis is to develop robust and accurate techniques for ionospheric correction of SAR and InSAR data and evaluate the benefit of these techniques for cryospheric research fields such as glacier ice velocity tracking and permafrost deformation monitoring. As both topics are mostly concerned with high latitude areas where the ionosphere is often active and characterized by turbulence, ionospheric correction is particularly relevant for these applications. After an introduction to the research topic in Chapter 1, Chapter 2 will discuss open issues in ionospheric correction including processing issues related to baseline-induced spectrum shifts. The effect of large baseline on split spectrum InSAR technique has been thoroughly evaluated and effective solutions for compensating this effect are proposed. In addition, a multiple sub-band approach is proposed for increasing the algorithm robustness and accuracy. Selected case studies are shown with the purpose of demonstrating the performance of the developed algorithm. In Chapter 3, the developed ionospheric correction technology is applied to optimize InSAR-based ice velocity measurements over the big ice sheets in Greenland and the Antarctic. Selected case studies are presented to demonstrate and validate the effectiveness of the proposed correction algorithms for ice velocity applications. It is shown that the ionosphere signal can be larger than the actual glacier motion signal in the interior of Greenland and Antarctic, emphasizing the necessity for operational ionospheric correction. The case studies also show that the accuracy of ice velocity estimates was significantly improved once the developed ionospheric correction techniques were integrated into the data processing flow. We demonstrate that the proposed ionosphere correction outperforms the traditionally-used approaches such as the averaging of multi-temporal data and the removal of obviously affected data sets. For instance, it is shown that about one hundred multi-temporal ice velocity estimates would need to be averaged to achieve the estimation accuracy of a single ionosphere-corrected measurement. In Chapter 4, we evaluate the necessity and benefit of ionospheric-correction for L-band InSAR-based permafrost research. In permafrost zones, InSAR-based surface deformation measurements are used together with geophysical models to estimate permafrost parameters such as active layer thickness, soil ice content, and permafrost degradation. Accurate error correction is needed to avoid biases in the estimated parameters and their co-variance properties. Through statistical analyses of a large number of L-band InSAR data sets over Alaska, we show that ionospheric signal distortions, at different levels of magnitude, are present in almost every InSAR dataset acquired in permafrost-affected regions. We analyze the ionospheric correction performance that can be achieved in permafrost zones by statistically analyzing correction results for large number of InSAR data. We also investigate the impact of ionospheric correction on the performance of the two main InSAR approaches that are used in permafrost zones: (1) we show the importance of ionospheric correction for permafrost deformation estimation from discrete InSAR observations; (2) we demonstrate that ionospheric correction leads to significant improvements in the accuracy of time-series InSAR-based permafrost products. Chapter 5 summarizes the work conducted in this dissertation and proposes next steps in this field of research

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

Generation of Earth’s Surface Three-Dimensional (3-D) Displacement Time-Series by Multiple-Platform SAR Data

In this chapter, the recent advancements of differential synthetic aperture radar interferometry (DInSAR) technique are presented, with the focus on the DInSAR-based approaches leading to the generation of three-dimensional time-series of Earth’s surface deformation, based on the combination of multi-platform line-of-sight (LOS)-projected time-series of deformation. Use of pixel-offset (PO) measurements for the retrieval of North-South deformation components, which are difficult to be extracted from DInSAR data, only, is also discussed. A review of the principal techniques based on the exploitation of amplitude and phase signatures of sequences of SAR images will be first provided, by emphasizing the limitations and strength of each single approach. Then, the interest will be concentrated on the recently proposed multi-track InSAR combination algorithm, referred as minimum acceleration InSAR combination (MinA) approach. The algorithm assumes the availability of two (or more) sets of SAR images acquired from complementary tracks. SAR data are pre-processed through one of currently available multi-temporal DInSAR toolboxes, and the LOS-projected surface deformation time-series are computed. An under-determined system of linear equations is then solved, based on imposing that the 3-D displacement time-series have minimum acceleration (MA). The presented results demonstrate the validity of the MinA algorithm

The Application of ALOS/PALSAR InSAR to Measure Subsurface Penetration Depths in Deserts

Spaceborne Synthetic Aperture Radar (SAR) interferometry has been utilised to acquire high-resolution Digital Elevation Models (DEMs) with wide coverage, particularly for persistently cloud-covered regions where stereophotogrammetry is hard to apply. Since the discovery of sand buried drainage systems by the Shuttle Imaging Radar-A (SIR-A) L-band mission in 1982, radar images have been exploited to map subsurface features beneath a sandy cover of extremely low loss and low bulk humidity in some hyper-arid regions such as from the Japanese Earth Resources Satellite 1 (JERS-1) and Advanced Land Observing Satellite/Phased Array type L-band Synthetic Aperture Radar (ALOS/PALSAR). Therefore, we hypothesise that a Digital Elevation Model (DEM) derived by InSAR in hyper-arid regions is likely to represent a subsurface elevation model, especially for lower frequency radar systems, such as the L-band system (1.25 GHz). In this paper, we compare the surface appearance of radar images (L-band and C-band) with that of optical images to demonstrate their different abilities to show subsurface features. Moreover, we present an application of L-band InSAR to measure penetration depths in the eastern Sahara Desert. We demonstrate how the retrieved L-band InSAR DEM appears to be of a consistently 1–2 m lower elevation than the C-band Shuttle Radar Topography Mission (SRTM) DEM over sandy covered areas, which indicates the occurrence of penetration and confirms previous studies

SAR interferometry at Venus for topography and change detection

AbstractSince the Magellan radar mapping of Venus in the early 1990’s, techniques of synthetic aperture radar interferometry (InSAR) have become the standard approach to mapping topography and topographic change on Earth. Here we investigate a hypothetical radar mission to Venus that exploits these new methods. We focus on a single spacecraft repeat-pass InSAR mission and investigate the radar and mission parameters that would provide both high spatial resolution topography as well as the ability to detect subtle variations in the surface. Our preferred scenario is a longer-wavelength radar (S or L-band) placed in a near-circular orbit at 600km altitude. Using longer wavelengths minimizes the required radar bandwidth and thus the amount of data that will be transmitted back to earth; it relaxes orbital control and knowledge requirements. During the first mapping cycle a global topography map would be assembled from interferograms taken from adjacent orbits. This approach is viable due to the slow rotation rate of Venus, causing the interferometric baseline between adjacent orbits to vary from only 11km at the equator to zero at the inclination latitude. To overcome baseline decorrelation at lower latitudes, the center frequency of a repeated pass will be adjusted relative to the center frequency of its reference pass. During subsequent mapping cycles, small baseline SAR acquisitions will be used to search for surface decorrelation due to lava flows. While InSAR methods are used routinely on Earth, their application to Venus could be complicated by phase distortions caused by the thick Venus atmosphere

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

Constructing water vapor maps by fusing InSAR, GNSS and WRF data

The work aims at constructing maps of the total water vapor in the atmosphere by fusing InSAR, GNSS, and WRF data

Measuring velocities of a surge type glacier with SAR interferometry using ALOS-2 data

In recent years, in-situ measurements on Kongsvegen, a surge-type glacier located in the Kongsfjorden area, have showed an acceleration in the flow speeds of the glacier. This could indicate the onset of a surging event, which presents the opportunity to study the dynamics of a glacier surge using remote sensing techniques with in-situ data for reference. Synthetic aperture radar (SAR) is well suited for this, as it does not rely on the sun for illumination and is not obstructed by clouds. In addition, SAR can be used to measure displacement with high accuracy and resolution through the use of interferometric SAR (InSAR). This study investigates the acceleration of Kongsvegen using InSAR, MAI and offset tracking. Velocity measurements from the combination DInSAR - MAI are then compared to in-situ data as well as the offset tracking measurements. For image pairs where InSAR measurements are not possible due to phase decorrelation, offset tracking is attempted as a back-up. Data from 2015, 2018 and 2019 was available, and the evolution of flow speeds over time could therefore be evaluated. The image pairs from 2018-2019 were acquired with 14 days separation in time, while the 2015 image pairs were acquired with 28 and 42 days separation. Due to the longer separation in time, the 2015 image pairs decorrelated in time. In addition, a pair acquired in the summer of 2018 decorrelated as a result of surface melting on the glaciers. Therefore only 3 of the total 8 pairs available were suited for interferometric analysis. For the image pairs from 2018-2019, the InSAR measurements were in good agreement with the in-situ data, as they also indicated an acceleration of the flow speeds on Kongsvegen. The offset tracking results on these pairs overestimated the velocity magnitudes, but also showed an increase in time. Similar to the InSAR estimates, the offset tracking failed to produce reasonable results on the 2015 image pairs, likely because of the large temporal baseline and lack of surface features on Kongsvegen. Overall, InSAR could be used to measure flow speeds on Kongsvegen successfully, but more data with a short temporal baseline is needed for an in-depth analysis