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

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

Challenges in Arctic Navigation and Geospatial Data : User Perspective and Solutions Roadmap

Navigation and location-based applications, including business such as transport, tourism, and mining, in Arctic areas face a variety of specific challenges. In fact, these challenges concern not only the Arctic Circle but certain other areas as well, such as the Gulf of Bothnia. This report provides a review on these challengs which concern a variety of technologies ranging from satellite navigation to telecommunications and mapping. In order to find out end-users' views on the significance of Arctic challenges, an online survey was conducted. The 77 respondents representing all Arctic countries, the majority being from Finland, highlighted the challenges in telecommunications as well as accuracy concerns for emerging applications dealing with precise navigation. This report provides a review of possible technologies for addressing the Arctic challenges, based on which a road map for solving them is developed. The road map also uses the results of expert working groups from the Challenges in Arctic Navigation workshop arranged in April 2018 in Olos, Muonio, Finland. This report was produced within the ARKKI project. It was funded by the Finnish Ministry of Foreign Affairs under the Baltic Sea, Barents and Arctic cooperation programme, and implemented by the Finnish Geospatial Research Institute in collaboration with the Finnish Ministry of Transport and Communications

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

GNSS transpolar earth reflectometry exploriNg system (G-TERN): mission concept

The global navigation satellite system (GNSS) Transpolar Earth Reflectometry exploriNg system (G-TERN) was proposed in response to ESA's Earth Explorer 9 revised call by a team of 33 multi-disciplinary scientists. The primary objective of the mission is to quantify at high spatio-temporal resolution crucial characteristics, processes and interactions between sea ice, and other Earth system components in order to advance the understanding and prediction of climate change and its impacts on the environment and society. The objective is articulated through three key questions. 1) In a rapidly changing Arctic regime and under the resilient Antarctic sea ice trend, how will highly dynamic forcings and couplings between the various components of the ocean, atmosphere, and cryosphere modify or influence the processes governing the characteristics of the sea ice cover (ice production, growth, deformation, and melt)? 2) What are the impacts of extreme events and feedback mechanisms on sea ice evolution? 3) What are the effects of the cryosphere behaviors, either rapidly changing or resiliently stable, on the global oceanic and atmospheric circulation and mid-latitude extreme events? To contribute answering these questions, G-TERN will measure key parameters of the sea ice, the oceans, and the atmosphere with frequent and dense coverage over polar areas, becoming a “dynamic mapper”of the ice conditions, the ice production, and the loss in multiple time and space scales, and surrounding environment. Over polar areas, the G-TERN will measure sea ice surface elevation (<;10 cm precision), roughness, and polarimetry aspects at 30-km resolution and 3-days full coverage. G-TERN will implement the interferometric GNSS reflectometry concept, from a single satellite in near-polar orbit with capability for 12 simultaneous observations. Unlike currently orbiting GNSS reflectometry missions, the G-TERN uses the full GNSS available bandwidth to improve its ranging measurements. The lifetime would be 2025-2030 or optimally 2025-2035, covering key stages of the transition toward a nearly ice-free Arctic Ocean in summer. This paper describes the mission objectives, it reviews its measurement techniques, summarizes the suggested implementation, and finally, it estimates the expected performance.Peer ReviewedPostprint (published version

The applications of satellites to communications, navigation and surveillance for aircraft operating over the contiguous United States. Volume 1 - Technical report

Satellite applications to aircraft communications, navigation, and surveillance over US including synthesized satellite network and aircraft equipment for air traffic contro

Altimetric system: Earth observing system. Volume 2h: Panel report

A rationale and recommendations for planning, implementing, and operating an altimetric system aboard the Earth observing system (Eos) spacecraft is provided. In keeping with the recommendations of the Eos Science and Mission Requirements Working Group, a complete altimetric system is defined that is capable of perpetuating the data set to be derived from TOPEX/Poseidon, enabling key scientific questions to be addressed. Since the scientific utility and technical maturity of spaceborne radar altimeters is well documented, the discussion is limited to highlighting those Eos-specific considerations that materially impact upon radar altimetric measurements

Detection of F-region electron density irregularities using incoherent-scatter radar

Thesis (M.S.) University of Alaska Fairbanks, 2014.Incoherent-scatter radar data from Poker Flat, Alaska has been used to determine size distributions of electron density structures in the evening time sector of the auroral zone. At high latitudes ionospheric plasma typically moves east-west with speeds of several hundred meters per second. Density irregularities that rapidly move through the radar beam are therefore observed as time-varying power fluctuations. The new phased array radar used for this study has been operated with several antenna directions with successive pulses transmitted in each direction. It is therefore possible to observe plasma Doppler velocities in multiple directions and determine the vector direction of the plasma motion. This near-simultaneous observation of the plasma velocity in conjunction with the electron density height profile data enable a new technique to determine the scale sizes of electron density fluctuations that move horizontally through the radar beam. The study focuses on the collision-less F-region ionosphere where the plasma drift is approximately constant with altitude. The experimental technique limits the range of scale sizes that may be studied to relatively large-scale sizes (i.e. greater than few tens of km). Results show that during magnetically disturbed conditions (Kp ≥ 4) when westward plasma velocities are relatively high (500-1000 m/s) the scale sizes of irregularities (often called plasma blobs) are in the range of 100-300 km and predominantly originate from the polar cap and are transported over long distances (~1000 km) due to the long chemical recombination times (30-90 minutes). Some irregularities are caused by local auroral particle precipitation and have been identified with associated electron temperature enhancements. For cases of low magnetic activity (Kp ≤ 1), when the radar is located in a region of low plasma velocities (100-500 m/s) well south of the auroral oval (essentially a mid-latitude type ionosphere), the density distribution is always biased strongly toward small-scale sizes (less than 50 km).Chapter 1. Introduction -- 1.1. Overview of the Ionosphere -- 1.2. Impact of Electron Density Irregularities on Navigation Systems and SAR Imaging -- 1.3. Historical Review and Objective of Thesis -- Chapter 2. Incoherent scatter radar -- 2.1. Introduction -- 2.2. Determination of Electron Density -- 2.3. Basic Experiment Design and Estimation of Line-of-sight Velocity -- 2.4. Determination of Ion Vector Velocities -- Chapter 3. Experimental approach -- 3.1. Introduction -- 3.2. AMISR System Description -- 3.3. Description of Experiments -- 3.4. Determination of Ionization Blob Scale Sizes -- Chapter 4. Observations of electron density irregularities using the poker flat incoherent scatter radar -- 4.1. Introduction -- 4.2. Results -- 4.2.1. August 28th, 2011 -- 4.2.2. August 31st, 2011 -- 4.2.3. March 9th, 2012 -- 4.2.4. March 25th, 2012 -- 4.2.5. October 9th, 2012 -- 4.2.6. October 12th, 2012 -- 4.2.7. October 19th, 2012 -- 4.2.8. November 4th, 2012 -- 4.3. Discussion -- Chapter 5. Conclusions and future work -- 5.1. Conclusions -- 5.2. Future Work

Observing pulsars and fast transients with LOFAR

Low frequency radio waves, while challenging to observe,are a rich source of information about pulsars. The LOw Frequency ARray (LOFAR) is a new radio interferometer operating in the lowest 4 octaves of the ionospheric “radio window”: 10–240 MHz, that will greatly facilitate observing pulsars at low radio frequencies. Through the huge collecting area, long baselines, and flexible digital hardware, it is expected that LOFAR will revolutionize radio astronomy at the lowest frequencies visible from Earth.LOFAR is a next-generation radio telescope and a pathfinder to the Square Kilometre Array (SKA), in that it incorporates advanced multi-beaming techniques between thousands of individual elements. We discuss the motivation for low-frequency pulsar observations in general and the potential of LOFAR in addressing these science goals.We present LOFAR as it is designed to perform high-time-resolution observations of pulsars and other fast transients, and outline the various relevant observing modes and data reduction pipelines that are already or will soon be implemented to facilitate these observations. A number of results obtained from commissioning observations are presented to demonstrate the exciting potential of the telescope. This paper outlines the case for low frequency pulsar observations and is also intended to serve as a reference for upcoming pulsar/fast transient science papers with LOFAR

Adaptive, reliable, and accurate positioning model for location-based services

This thesis was submitted for the degree of Doctor of Philosophy and awarded by Brunel University.This thesis presents a new strategy in achieving highly reliable and accurate position solutions fulfilling the requirements of Location-Based Services (LBS) pedestrians’ applications. The new strategy is divided into two main parts. The first part integrates the available positioning technology within the surrounding LBS application context by introducing an adaptive LBS framework. The context can be described as a group of factors affecting the application behaviour; this includes environmental states, available resources and user preferences. The proposed adaptive framework consists of several stages, such as defining the contextual factors that have a direct effect on the positioning performance, identifying preliminary positioning performance requirements associated with different LBS application groups, and introducing an intelligent positioning services selection function. The second part of this work involves the design and development of a novel positioning model that is responsible for delivering highly reliable, accurate and precise position solutions to LBS users. This new model is based on the single frequency GPS Standard Positioning Service (SPS). Additionally, it is incorporated within the adaptive LBS framework while providing the position solutions, in which all identified contextual factors and application requirements are accounted. The positioning model operates over a client-server architecture including two main components, described as the Localisation Server (LS) and the Mobile Unit (MU). Hybrid functional approaches were developed at both components consisting of several processing procedures allowing the positioning model to operate in two position determination modes. Stand-alone mode is used if enough navigation information was available at the MU using its local positioning device (GPS/EGNOS receiver). Otherwise, server-based mode is utilised, in which the LS intervenes and starts providing the required position solutions. At the LS, multiple sources of GPS augmentation services were received using the Internet as the sole augmentation data transportation medium. The augmentation data was then processed and integrated for the purpose of guaranteeing the availability of valid and reliable information required for the provision of accurate and precise position solutions. Two main advanced position computation methods were developed at the LS, described as coordinate domain and raw domain. The positioning model was experimentally evaluated. According to the reported results, the LS through the developed position computation methods, was able to provide position samples with an accuracy of less than 2 meters, with high precision at 95% confidence level; this was achieved in urban, rural, and open space (clear satellite view) navigation environments. Additionally, the integrity of the position solutions was guaranteed in such environments during more than 90% of the navigation time, taking into consideration the identified integrity thresholds (Horizontal Alert Limits (HAL)=11 m). This positioning performance has outperformed the existing GPS/EGNOS service which was implemented at the MU in all scenarios and environments. In addition, utilising a simulation evaluation facility the developed positioning model performance was quantified with reference to a hybrid positioning service that will be offered by future Galileo Open Service (OS) along with GPS/EGNOS. Using the statistical t-test, it was concluded that there is no significant difference in terms of the position samples’ accuracy achieved from the developed positioning model and the hybrid system at a particular navigation environment described as rural area. The p-value was 0.08 and the level of significance used was 0.05. However, a significant difference in terms of the service integrity for the advantage of the hybrid system was experienced in all remaining scenarios and environments more especially the urban areas due to surrounding obstacles and conditions