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

Extrapolation of Airborne Polarimetric and Interferometric SAR Data for Validation of Bio-Geo-Retrieval Algorithms for Future Spaceborne SAR Missions

Spaceborne SAR system concepts and mission design is often based on algorithms developed and the experience gathered from airborne SAR experiments and associated dedicated campaigns. However, airborne SAR systems have better performance parameters than their future space-borne counterparts as their design is not impacted by mass, power, and storage constraints. This paper describes a methodology to extrapolate spaceborne quality SAR image products from long wavelength airborne polarimetric SAR data which were acquired especially for the development and validation of bio/geo-retrieval algorithms in forested regions. For this purpose not only system (sensor) related parameters are altered, but also those relating to the propagation path (ionosphere) and to temporal decorrelation

Toward Operational Compensation of Ionospheric Effects in SAR Interferograms: The Split-Spectrum Method

The differential ionospheric path delay is a major error source in L-band interferograms. It is superimposed to topography and ground deformation signals, hindering the measurement of geophysical processes. In this paper, we proceed toward the realization of an operational processor to compensate the ionospheric effects in interferograms. The processor should be robust and accurate to meet the scientific requirements for the measurement of geophysical processes, and it should be applicable on a global scale. An implementation of the split-spectrum method, which will be one element of the processor, is presented in detail, and its performance is analyzed. The method is based on the dispersive nature of the ionosphere and separates the ionospheric component of the interferometric phase from the nondispersive component related to topography, ground motion, and tropospheric path delay. We tested the method using various Advanced Land Observing Satellite Phased-Array type L-band synthetic aperture radar interferometric pairs with different characteristics: high to low coherence, moving and nonmoving terrains, with and without topography, and different ionosphere states. Ionospheric errors of almost 1 m have been corrected to a centimeter or a millimeter level. The results show how the method is able to systematically compensate the ionospheric phase in interferograms, with the expected accuracy, and can therefore be a valid element of the operational processor

Imaging ionospheric inhomogeneities using spaceborne synthetic aperture radar

We present a technique and results of 2-D imaging of Faraday rotation and total electron content using spaceborne L band polarimetric synthetic aperture radar (PolSAR). The results are obtained by processing PolSAR data collected using the Phased Array type L-band Synthetic Aperture Radar (PALSAR) on board the Advanced Land Observation Satellite. Distinguished ionospheric inhomogeneities are captured in 2-D images from space with relatively high resolutions of hundreds of meters to a couple of kilometers in auroral-, middle-, and low-latitude regions. The observed phenomena include aurora-associated ionospheric enhancement arcs, the middle-latitude trough, traveling ionospheric disturbances, and plasma bubbles, as well as ionospheric irregularities. These demonstrate a new capability of spaceborne synthetic aperture radar that will not only provide measurements to correction of ionospheric effects in Earth science imagery but also significantly benefit ionospheric studies

Estimation of Forest Biomass and Faraday Rotation using Ultra High Frequency Synthetic Aperture Radar

Synthetic Aperture Radar (SAR) data in the Ultra High Frequency (UHF; 300 MHz – 3 GHz)) band have been shown to be strongly dependent of forest biomass, which is a poorly estimated variable in the global carbon cycle. In this thesis UHF-band SAR data from the fairly flat hemiboreal test site Remningstorp in southern Sweden were analysed. The data were collected on several occasions with different moisture conditions during the spring of 2007. Regression models for biomass estimation on stand level (0.5-9 ha) were developed for each date on which SAR data were acquired. For L-band (centre frequency 1.3 GHz) the best estimation model was based on HV-polarized backscatter, giving a root mean squared error (rmse) between 31% and 46% of the mean biomass. For P-band (centre frequency 340 MHz), regression models including HH, HV or HH and HV backscatter gave an rmse between 18% and 27%. Little or no saturation effects were observed up to 290 t/ha for P-band. A model based on physical-optics has been developed and was used to predict HH-polarized SAR data with frequencies from 20 MHz to 500 MHz from a set of vertical trunks standing on an undulating ground surface. The model shows that ground topography is a critical issue in SAR imaging for these frequencies. A regression model for biomass estimation which includes a correction for ground slope was developed using multi-polarized P-band SAR data from Remningstorp as well as from the boreal test site Krycklan in northern Sweden. The latter test site has pronounced topographic variability. It was shown that the model was able to partly compensate for moisture variability, and that the model gave an rmse of 22-33% when trained using data from Krycklan and evaluated using data from Remningstorp. Regression modelling based on P-band backscatter was also used to estimate biomass change using data acquired in Remningstorp during the spring 2007 and during the fall 2010. The results show that biomass change can be measured with an rmse of about 15% or 20 tons/ha. This suggests that not only deforestation, but also forest growth and degradation (e.g. thinning) can be measured using P-band SAR data. The thesis also includes result on Faraday rotation, which is an ionospheric effect which can have a significant impact on spaceborne UHF-band SAR images. Faraday rotation angles are estimated in spaceborne L-band SAR data. Estimates based on distributed targets and calibration targets with high signal to clutter ratios are found to be in very good agreement. Moreover, a strong correlation with independent measurements of Total Electron Content is found, further validating the estimates

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

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

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