Ionospheric Phase Compensation for InSAR Measurements Based on the Faraday Rotation Inversion Method.

Mapping high spatial resolution ionospheric total electron content by integrating Time Series InSAR with International Reference Ionosphere model

The ionospheric propagation path delay is a major error source in synthetic aperture radar (SAR) interferograms and, therefore, has to be estimated and corrected. Various methods can be used to extract different kinds of information about the ionosphere from SAR images, with different accuracies. This paper presents a general technique, based on a Bayesian inverse problem, that combines various information sources in order to increase the estimation accuracy, and thus the correction. A physically realistic fractal modeling of the ionosphere turbulence and a data-based estimation of the model parameters allow the avoidance of arbitrary filtering windows and coefficients. To test the technique, the differential ionospheric phase screen was estimated by combining the split-spectrum method with the azimuth mutual shifts between interferometric pair images. This combination is convenient since it can benefit from the strengths of both sources: range and azimuth variations from the split-spectrum method and small-scale azimuth variations from more sensitive azimuth shifts. Therefore, the two methods can recover the long and short wavelength components of the ionospheric phase screen, respectively. A theoretical comparison between the Faraday rotation method and the split-spectrum method is also reported. For the use in the combination, precedence was then given to the split-spectrum method because of the comparable precision level, lower susceptibility to biases, and wider applicability. Finally, Advanced Land Observing Satellite Phased Array type L-band SAR L-band images are used to show how the combined result is more accurate than that obtained with the simple split-spectrum method.

The L-band synthetic aperture radar (SAR) interferometry (InSAR) technique has a lower accuracy due to ionospheric phase distortions. Recently, a multiple-aperture interferometry (MAI)-based ionospheric correction method has been proposed. Using four types of ionosphere-distorted interferograms, the performance of the correction method was evaluated and then analyzed the feasibility of the correction method. The test interferograms contained severe azimuth streaking, low-frequency ionospheric phase distortion and drastic phase change due to the ionosphere. The results showed that (i) the existence and magnitude of ionospheric phase distortions can be recognized from MAI interferograms and (ii) the MAI-based ionospheric correction method efficiently reduced severe azimuth streaking and low-frequency distortion but did not mitigate the drastic phase change perfectly. It is allowed (i) to determine whether a given SAR interferogram has an ionospheric distortion and (ii) to predict whether the ionospheric distortion can be corrected by using the MAI-based ionospheric correction method.

Global ionospheric mapping (GIM) is a new and emerging technique for determining global ionospheric TEC (total electron content) based on measurements from a worldwide network of Global Positioning System (GPS) receivers. In this study, GIM accuracy in specifying TEC is investigated by comparison with direct ionospheric measurements from the TOPEX altimeter. A climatological model (Bent model) is also used to compare with the TOPEX altimeter data. We find that the GIM technique has much better agreement with TOPEX in TEC measurements, compared with the predictions of the climatological model. The difference between GIM and TOPEX in TEC measurements is very small (less than 1.5 TEC units (TECU)) within a 1500‐km range from a reference GPS station. The RMS gradually increases with increasing distance from the station, while the Bent model shows a constant large RMS, unrelated to any station location. Within a 1000‐km distance of a GPS site (elevation angle > 25°), GIM has a good correlation (R > 0.93) to TOPEX with respect to TEC measurements. The slope of the linear fitting line to the data set from two TOPEX cycles is 44.5° (near the ideal 45°). In the northern hemispheric regions, ionospheric specification by GIM appears to be accurate to within 3‐10 TECU up to 2000+ km away from nearest GPS station (corresponding to ∼1° elevation angle cutoff). Beyond 2000 km, GIM accuracy, on average, is reduced to the Bent model levels. In the equatorial region, the Bent model predictions are systematically lower (∼5.0 TECU) than TOPEX values and often show a saturation at large TEC values. During ionospheric disturbed periods, GIM sometimes shows differences from TOPEX values due to transient variations of the ionosphere. Such problems may be improved by the continuous addition of new GPS stations in data‐sparse regions. Thus, over a GPS station's measurement realm (up to 2000 km in radius), GIM can produce generally accurate TEC values. Through a spatial and temporal extrapolation of GPS‐derived TEC measurements, the GIM technique provides a powerful tool for monitoring global ionospheric features in near real time.

Research on longitudinal discrepancies in local ionospheric variability, especially in equatorial and low-latitude regions, is a focal point of interest for the space weather modeling community. The ionosphere over these regions is influenced by complex electrodynamics, wind, and temperature dynamics that can seriously impact dynamic technological systems such as satellite tracking and positioning, satellite radio communication, and navigation control systems. Here, we researched the longitudinal variability in the ionospheric total electron content (TEC) by analyzing observed global positioning system (GPS)-derived TEC values along with those extracted from the most reliable global ionospheric maps (GIMs) and the International Reference Ionosphere (IRI-2016) model at selected stations in the vicinity of the magnetic equator along the American, African, and Asian longitude sectors. The period of study covered the descending (2016–2017) and deep solar minimum (2018–2019) years in the 24th solar cycle. Apart from the decreasing trend of the TEC from the descending to deep solar minimum period irrespective of season and longitude sector, the results showed a relatively higher magnitude of TEC in the African longitude than the other two longitude sectors. Despite evident overestimation and underestimations of TEC in both models, GIM predictions generally looked better in terms of observed variation patterns, especially in the African longitude. The study also highlights the seasonal and semiannual effects of longitudinal variations in TEC, manifesting in local time offsets and some peculiar anomalies, which seemed to be different from previously reported results, especially during the solar minimum years at the three longitude sectors. The insignificant effects of longitudinal variations on the equinoctial asymmetry are attributed to the diverse electron density distribution and ionospheric morphology at the three longitude sectors that will prompt further investigations in the future. The outcomes from this study may augment the past efforts of scientists to understand the seasonal effects of the longitudinal variations in TEC, thereby complementing the improvements of ionospheric representations in global ionosphere models and maps.

The authors tested the applicability of coherence information for optimizing SAR (synthetic aperture radar) interferogram. The purpose of the optimization in this study is to improve the smoothness of the topographic fringe in a SAR interferogram, which results in better fringe for unwrapping. Coherence information was computed as the complex correlation coefficient between two SAR complex data. To optimize the interferogram, the maximum coherence value was searched in 3/spl times/5 window around the initial registered position between two data and then the phase difference was obtained from the pixel pair that has the maximum coherence. This optimizing method was attempted for several SAR complex data pairs. The effect of the optimization was found in some pairs. It was found that the increase of the coherence clue to the optimization was one of the important factors for the improvement of the smoothness.

Aside from the ionospheric total electron content (TEC) information, root-mean-square (RMS) maps are also provided as the standard deviations of the corresponding TEC errors in global ionospheric maps (GIMs). As the RMS maps are commonly used as the accuracy indicator of GIMs to optimize the stochastic model of precise point positioning algorithms, it is of crucial importance to investigate the reliability of RMS maps involved in GIMs of different Ionospheric Associated Analysis Centers (IAACs) of the International GNSS Service (IGS), i.e., the integrity of GIMs. We indirectly analyzed the reliability of RMS maps by comparing the actual error of the differential STEC (dSTEC) with the RMS of the dSTEC derived from the RMS maps. With this method, the integrity of seven rapid IGS GIMs (UQRG, CORG, JPRG, WHRG, EHRG, EMRG, and IGRG) and six final GIMs (UPCG, CODG, JPLG, WHUG, ESAG and IGSG) was examined under the maximum and minimum solar activity conditions as well as the geomagnetic storm period. The results reveal that the reliability of the RMS maps is significantly different for the GIMs from different IAACs. Among these GIMs, the values in the RMS maps of UQRG are large, which can be used as ionospheric protection level, while the RMS values in EHRG and ESAG are significantly lower than the realistic RMS. The rapid and final GIMs from CODE, JPL and WHU provide quite reasonable RMS maps. The bounding performance of RMS maps can be influenced by the location of the stations, while the influence of solar activity and the geomagnetic storm is not obvious.

Research on total electron content (TEC) empirical models is one of the important topics in the field of space weather services. Global TEC empirical models based on Global Ionospheric Maps (GIMs) TEC data released by the International GNSS Service (IGS) have developed rapidly in recent years. However, the accuracy of such global empirical models has a crucial restriction arising from the non-uniform accuracy of IGS TEC data in the global scope. Specifically, IGS TEC data accuracy is higher on land and lower over the ocean due to the lack of stations in the latter. Using uneven precision GIMs TEC data as a whole for model fitting is unreasonable. Aiming at the limitation of global ionospheric TEC modelling, this paper proposes a new global ionospheric TEC empirical model named the TECM-GRID model. The model consists of 5183 sections, corresponding to 5183 grid points (longitude 5°, latitude 2.5°) of GIM. Two kinds of single point empirical TEC models, SSM-T1 and SSM-T2, are used for TECM-GRID. According to the locations of grid points, the SSM-T2 model is selected as the sub-model in the Mid-Latitude Summer Night Anomaly (MSNA) region, and SSM-T1 is selected as the sub-model in other regions. The fitting ability of the TECM-GRID model for modelling data was tested in accordance with root mean square (RMS) and relative RMS values. Then, the TECM-GRID model was validated and compared with the NTCM-GL model and Center for Orbit Determination in Europe (CODE) GIMs at time points other than modelling time. Results show that TECM-GRID can effectively describe the Equatorial Ionization Anomaly (EIA) and the MSNA phenomena of the ionosphere, which puts it in good agreement with CODE GIMs and means that it has better prediction ability than the NTCM-GL model.

Interferometric synthetic aperture radar (InSAR) products may be significantly distorted by microwave signals traveling through the ionosphere, especially with long wavelengths. The split-spectrum method (SSM) is used to separate the ionospheric and the nondispersive phase terms with lower and higher spectral sub-band interferogram images. However, the ionospheric path delay phase is very delicate to the synthetic aperture radar (SAR) parameters including orbit vectors, slant range, and target height. In this paper, we get the impact of SAR parameter errors on the ionospheric phase by two steps. The first step is getting the derivates of geolocation with reference to SAR parameters based on the range-Doppler (RD) imaging model and the second step is calculating the derivates of the ionospheric phase delay with respect to geometric positioning. Through the numerical simulation, we demonstrate that the deviation of ionospheric phase has a linear relationship with SAR parameter errors. The experimental results show that the estimation of SAR parameters should be accurate enough since the parameter errors significantly affect the performance of ionospheric correction. The root mean square error (RMSE) between the corrected differential interferometric SAR (DInSAR) phase with SAR parameter errors and the corrected DInSAR phase without parameter errors varies from centimeter to decimeter level with the L-band data acquired by the Advanced Land Observing Satellite (ALOS) Phased Array type L-band SAR (PALSAR) over Antofagasta, Chile. Furthermore, the effectiveness of SSM can be improved when SAR parameters are accurately estimated.

Obtaining ionospheric Total Electron Content (TEC) with high spatial and temporal resolution is crucial for satellite navigation and space weather monitoring. Global Ionospheric Maps (GIMs) are key sources of TEC data. However, the temporal resolution of GIMs from ionosphere associate analysis centers does not satisfy the requirements for space weather monitoring. To address this issue, the Massachusetts Institute of Technology (MIT) Haystack Observatory has provided high‐time‐resolution TEC maps since 1998, using data from ∼5,000 Global Navigation Satellite System (GNSS) tracking stations. However, our research reveals a significant systematic bias in MIT's TEC maps. Unlike the biases between GIMs, it exhibits an unexpected pattern of variation and is also present in comparisons with radio occultation TEC. Further analysis indicates that the systematic bias in MIT TEC decreases with increasing geomagnetic latitude and is correlated with solar activity. In low‐latitude regions, the bias can reach up to 30 TECU, while in high‐latitude regions, it can drop to nearly 0 TECU. During diurnal and seasonal variations, the systematic bias in MIT TEC follows TEC patterns but appears unaffected by small‐scale ionospheric structures. Finally, we analyze the factors potentially contributing to the systematic bias and find that Differential Code Biases (DCBs) provided by MIT exhibits errors of a magnitude comparable to the TEC bias, with a strong correlation between the two. Through further investigation, we confirm that the systematic bias in MIT TEC primarily stems from its DCB estimation strategy and choice of mapping function, and this bias persists both before and after MIT updated its DCB estimation approach.

The ionosphere widely affects Global Navigation Satellite Systems (GNSS) applications, inducing among others a delay in GNSS measurements. This delay is closely linked to the Total Electron Content (TEC) of the ionosphere, a major parameter which can hence be monitored using GNSS. To this extent, phase measurements are taken as a basis for their lower noise level. Levelling strategies have then to be defined for the phase measurements are obtained with an initial unknown number of cycles called ambiguity. The most common technique, referred to as carrier-to-code levelling, consists in using the differences between code and phase measurements and their average on a continuous set of epochs. This option, chosen at the Royal Meteorological Institute (RMI) of Belgium to compute TEC for Belgian GPS stations, requires code hardware delays estimation. Another has been proposed which takes benefit from Global Ionospheric Maps (GIMs) to compute a reference TEC used for ambiguity resolution. In order to understand the consequences of using one method or the other, we compare slant TEC data obtained from both techniques for a mid-latitude station (Brussels) during a high solar activity period (2002). We observed large differences (6.8 TECu on average) showing features apparently related to ionospheric and geomagnetic activity. We attribute these observations to a combination of effects originating in code delays estimation, multipath and noise as well as GIMs errors. We try to differentiate between these effects by focusing on several days and satellites. We concentrate for example on days presenting large TEC differences and geomagnetic disturbances simultaneously (or not) or on satellites displaying recurrent patterns on consecutive days. Finally we highlight the impact of the choice of GIMs involved in sTEC calibration. To this extent, we analyse vertical TEC statistics showing a general underestimation from RMI data. The highest bias (5.8 TECu) is obtained for the UPC GIMs used in the second levelling technique.

Temporal extrapolation of TEC using WNN during 2007–2018 and comparison against GIM, IRI2016 and Kriging

Global ionosphere models may not be accurate enough in regions where there are insufficient observations. Therefore, there is strong interest in regional ionosphere modeling using observations from dual-frequency Global Navigation Satellite System (GNSS) receivers in such regions. This paper presents a data assimilation method for hourly total electron content (TEC) modeling along with estimation of the receivers’ differential code biases (DCBs). For this purpose, first, initial values for slant TEC (STEC) and the receivers’ DCBs are estimated utilizing global ionosphere maps (GIMs). Second, STEC is expanded to a linear combination of spherical radial basis functions with time-dependent coefficients to reduce the number of unknowns, and a constrained least-squares adjustment is applied to resolve the corrected STEC and DCB values in the study area using observations from dual-frequency GNSS receivers. Third, the resolved STEC at ionospheric pierce points (IPPs) are converted to vertical TEC (VTEC) using a mapping function and assimilated into the background provided by spatial and temporal interpolation of GIM. Retrieved VTECs are weighted according to the satellites’ elevation angles, and an objective function is defined to determine the background’s covariance matrix. In this study, the Iranian permanent Global Positioning System (GPS) network is used to construct the regional hourly VTEC model and calculate hourly DCB variations. The proposed method is validated in two ways. First, some of the VTEC values are excluded from the dataset for use as test data and compared with modeling results. The root mean square (RMS) of the constructed regional assimilative model on test data is significantly better than GIM. Second, the RINEX file of Tehran station, whose observations did not play a role in modeling, is corrected using the regional model. The single-frequency positioning results with corrected observations show much better accuracy than GIMs.

Assessment of ALOS-2 PALSAR-2L-band and Sentinel-1 C-band SAR backscatter for discriminating between large-scale oil palm plantations and smallholdings on tropical peatlands

Investigation of the variability of total electron content (TEC) is one of the most important parameters of the observation and monitoring of space weather, which is the main cause of signal disturbance in space-based communication, positioning, and navigation systems. TEC is defined as the total number of electrons on a ray path. The Global Positioning System (GPS) provides a cost-effective solution for the estimation of TEC. Due to various physical and operational disturbances, TEC may have temporal and spatial domain gaps. Global ionospheric maps (GIMs) provide worldwide TEC with 1- to 2-h temporal resolution and $2.5^{\circ}\, \times 5^{\circ}$ spatial resolution in latitude and longitude, respectively. The GIM-TEC with the highest possible accuracy can be obtained 10 days after the recording of the signals. Therefore, a high-resolution and accurate interpolation of TEC is necessary to image and monitor the regional distribution of TEC in near-real time. In this study, a novel spatiotemporal interpolation algorithm with automatic gridding is developed for 2-D TEC imaging by data fusion of GPS-TEC and GIM-TEC. The algorithm automatically implements optimum spatial resolution and desired temporal resolution with universal kriging with linear trend for midlatitude regions and ordinary kriging for other regions. The theoretical semivariogram function is estimated from GPS network data using a Matern family, whose parameters are determined with a particle swarm optimization algorithm. The developed algorithm is applied to the Turkish National Permanent GPS Network (TNPGN-Active), a dense midlatitude GPS network. For the first time in the literature, high spatial resolution TEC maps are obtained between May 2009 and May 2012 with a 2.5-min temporal update period. These TEC maps will be used to investigate the spatiotemporal variability of the ionosphere over the diurnal and annual trend structure, including seasonal anomalies and geomagnetic and seismic disturbances over ionosphere.