Correction of Ionosphere for InSAR by the Combination of Differential TEC Estimators

Low frequency spaceborne SAR configurations are favoured for global forest mapping applications and D-InSAR applications over natural terrain. Several missions have been scheduled to be launched / or proposed to be implemented in the next years: JAXA’s ALOS-II (L-band), NASA’s Destyni (L-band), DLR’s Tandem-L (L-band) and ESA’s BIOMASS (P-band) are some of them. A common challenge for all these missions is to control / compensate the disturbances induced by the ionosphere. At these lower frequencies the ionosphere effects several components of the SAR measurements performed: It delays the group velocity of the transmitting / receiving pulses, advances their phase(s) and rotates their polarisation state. Accordingly, it distorts not only intensity but also polarimetric, interferometric and polarimetric interferometric observation spaces. The total electron content (TEC) is the most decisive parameter in the characterisation of the ionosphere. It is defined as the integrated electron number density per unit volume along the direction of propagation. Most of the free electrons are distributed within a relatively narrow altitude range allowing modelling the ionosphere as a thin layer at a fixed altitude. In this case the ionosphere can be characterised by a 2-D scalar field of TEC [1], [2]. Depending now on the SAR configuration and its observation space different correction approaches are possible leading to a wide range of calibration algorithms. In this paper we propose a concept towards the generalisation of ionospheric calibration methodology by integrating a number of individual approaches / algorithms. In this sense, a novel generic correction schema based on a combined (and improved) estimation of the 2-D TEC field (or the associated differential TEC field in the interferometric case) from a set of individual data based TEC and/or TEC gradient estimates is introduced and discussed. As a special case a combined 2-D (differential) TEC field estimator based on (differential) TEC estimated from Faraday rotation measurements and (differential) TEC gradients obtained from the estimation of azimuth/range (differential) shifts is presented. Both observations are independent, allowing establishing an inverse problem for the (differential) TEC estimation. Geophysical knowledge as the anisotropic nature of the TEC distribution can be incorporated as a priori information in the “combined” (differential) TEC estimator. The performance of the proposed approach is tested using ALOS quad-pol interferometric data sets over several test sites in Alaska. The achieved estimates are characterised by a significantly improved performance: While the FR based estimator suffers from the random granular deviation pattern of TEC after conversion, the proposed combined estimator effectively is free of such artefacts. Emphasis is given in the role of polarisation in the TEC estimation procedure [3] and on the calibration of Pol-InSAR data. References [1] Franz J. Meyer and Jeremy Nicoll, “Prediction, detection, and correction of Faraday rotation in full-polarimetric L-band SAR data”, IEEE Trans. Geosci. And Remote Sensing, 46(10), Oct., 3076-3086, 2008 [2] Xiaoqing Pi, Anthony Freeman, Bruce Champman, Paul Rosen, and Zhenhong Li, “Imaging ionospheric inhomogeneities using spaceborne synthetic aperture radar”, Jour. of Geophysical Research, 116, A04303, 2011 [3] Jun Su Kim, Konstantinos Papathanassiou, Shaun Quegan and Neil Rogers, “Estimation and correction of scintillation effects on spaceborne P-band SAR images”, in Proceedings of IGARSS2012, 23-27. Jul., 201

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

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

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

SAR Observation of Ionosphere Using Range/Azimuth Sub-Bands

In SAR applications, the ionosphere is normally considered as a disturbance that has to be removed or compensated for. In this paper the ionosphere is the objective of observation. This is motivated by the fact that nowadays polarimetric SAR systems can provide high resolution ionospheric maps that are not possible using conventional ionospheric mapping tools. This paper investigates a set of new ionospheric parameters that can be observed and mapped by SAR by exploring range and azimuth sub-bands. The range sub-bands are used to estimate the ionosphere independently of polarimetry. The azimuth sub-bands are related to the three-dimensional geometry of the ionosphere, as well as, to its dynamic component. The potential to estimate these parameters from SAR sub-bands is discussed and first results using ALOS PALSAR data are presented

The European Space Agency BIOMASS mission: Measuring forest above-ground biomass from space

The primary objective of the European Space Agency's 7th Earth Explorer mission, BIOMASS, is to determine the worldwide distribution of forest above-ground biomass (AGB) in order to reduce the major uncertainties in calculations of carbon stocks and fluxes associated with the terrestrial biosphere, including carbon fluxes associated with Land Use Change, forest degradation and forest regrowth. To meet this objective it will carry, for the first time in space, a fully polarimetric P-band synthetic aperture radar (SAR). Three main products will be provided: global maps of both AGB and forest height, with a spatial resolution of 200 m, and maps of severe forest disturbance at 50 m resolution (where “global” is to be understood as subject to Space Object tracking radar restrictions). After launch in 2022, there will be a 3-month commissioning phase, followed by a 14-month phase during which there will be global coverage by SAR tomography. In the succeeding interferometric phase, global polarimetric interferometry Pol-InSAR coverage will be achieved every 7 months up to the end of the 5-year mission. Both Pol-InSAR and TomoSAR will be used to eliminate scattering from the ground (both direct and double bounce backscatter) in forests. In dense tropical forests AGB can then be estimated from the remaining volume scattering using non-linear inversion of a backscattering model. Airborne campaigns in the tropics also indicate that AGB is highly correlated with the backscatter from around 30 m above the ground, as measured by tomography. In contrast, double bounce scattering appears to carry important information about the AGB of boreal forests, so ground cancellation may not be appropriate and the best approach for such forests remains to be finalized. Several methods to exploit these new data in carbon cycle calculations have already been demonstrated. In addition, major mutual gains will be made by combining BIOMASS data with data from other missions that will measure forest biomass, structure, height and change, including the NASA Global Ecosystem Dynamics Investigation lidar deployed on the International Space Station after its launch in December 2018, and the NASA-ISRO NISAR L- and S-band SAR, due for launch in 2022. More generally, space-based measurements of biomass are a core component of a carbon cycle observation and modelling strategy developed by the Group on Earth Observations. Secondary objectives of the mission include imaging of sub-surface geological structures in arid environments, generation of a true Digital Terrain Model without biases caused by forest cover, and measurement of glacier and icesheet velocities. In addition, the operations needed for ionospheric correction of the data will allow very sensitive estimates of ionospheric Total Electron Content and its changes along the dawn-dusk orbit of the mission

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