Refinement of global ionospheric coefficients for GNSS applications: Methodology and results

We developed the methodology for the optimal estimation of global ionospheric coefficients of the current Global Navigation Satellite Systems (GNSSs), including the eight- and ten-parameter Klobuchar-like as well as NeQuick models. The ionospheric coefficients of those correction models are calculated from two sets of globally distributed tracking stations of the International GNSS Services (IGS). Performance of the re-estimated Klobuchar-like and NeQuick coefficients are validated during 2002–2014 over the continental and oceanic areas, respectively. Over the continental areas, GPS TECs derived from 40 ground GPS receivers are selected as reference. The eight-, ten-parameter Klobuchar-like and NeQuick models can mitigate the ionospheric delay by 65.8, 67.3 and 75.0%, respectively. Over the global oceans, the independent TECs derived from Jason-1&2 altimeters are used as reference. The re-estimated ionospheric correction models can mitigate 56.1–66.7% of the delay errors. Compared to the original GPS Ionospheric Correction Algorithm (ICA), performance of those eight-, ten-parameter Klobuchar-like and NeQuick models has improved 3.4, 5.9 and 13.4% during the whole test period, respectively. The methodology developed here takes the advantage of high-quality ionospheric TECs derived from the global network of GNSS receivers. The re-estimated ionospheric coefficients can be used as precise ionospheric products to monitor and assess GNSS broadcast ionospheric parameters and to improve the performance of various single-frequency GNSS applications.     

The multi-source data fusion global ionospheric modeling software—IonoGim

We introduce a new global ionospheric modeling software—IonoGim, using ground-based GNSS data, the altimetry satellite and LEO (Low Earth Orbit) occultation data to establish the global ionospheric model. The software is programmed by C++ with fast computing speed and highly automatic degree, it is especially suitable for automatic ionosphere modeling. The global ionospheric model and DCBs obtained from IonoGim were compared with the CODE (Center for Orbit Determination in Europe) to verify its accuracy and reliability. The results show that IonoGim and CODE have good agreement with small difference, indicating that IonoGim owns high accuracy and reliability, and can be fully applicable for high-precision ionospheric research. In addition, through comparison between only using ground-based GNSS observations and multi-source data model, it can be demonstrated that the space-based ionospheric data effectively improve the model precision in marine areas where the ground-based GNSS tracking station lacks.     

Ionospheric climate index as a driving parameter for the NeQuick model

Ionospheric delay is one of the significant error sources for global navigation satellite system (GNSS) positioning. GNSSs broadcast the coefficients of the ionospheric model to correct ionospheric delay for single-frequency users. A modified three-dimensional model (NeQuick G) based on the NeQuick climatological model is adopted for Galileo users. The NeQuick G model uses the effective ionization level (Az) instead of the sunspot number as the driving parameter. In this study, we introduce the ionospheric climate index (ICI) as a new driving parameter for the NeQuick model. In comparison, the ICI-driven NeQuick model has a better performance than the Az-driven NeQuick G model at both low and high latitudes. In addition, only one GNSS station at low latitudes is required to calculate the ICI, which would save maintenance costs and improve the efficiency of updating the broadcast coefficients. This model has potential application value for future upgrades of Galileo’s ionospheric broadcast model.     

Ionosphere modelling for Galileo single frequency users: Illustration of the combination of the NeQuick model and GNSS data ingestion

The ionospheric effect remains one of the main factors limiting the accuracy of Global Navigation Satellite Systems (GNSS) including Galileo. For single frequency users, this contribution to the error budget will be mitigated by an algorithm based on the NeQuick global ionospheric model. This quick-run empirical model provides flexible solutions for combining ionospheric information obtained from various sources, from GNSS to ionosondes and topside sounders. Hence it constitutes an interesting simulation tool not only serving Galileo needs for mitigation of the ionospheric effect but also widening the use of new data.     

GNSS电离层模型改正精度评估与分析

电离层时延误差是导航定位信号在空间传播路径上的主要误差源之一,因此全面了解GNSS电离层模型的改正精度具有一定现实意义.根据GPS,BDS和Galileo系统所采用的电离层修正模型,利用2014年电离层校正参数,以高精度全球电离层图为基准,评估分析了三大系统电离层时延的改正精度.结果表明:目前GNSS使用的几种电离层修正模型的改正率在65~75%左右;Galileo系统使用的第二版NeQuick模型与第一版NeQuick模型相比在修正精度上并无显著提高;GPS使用的Klobuchar 8参数模型在北半球25°-45°N的中纬度地区精度很高,但是在全球其他区域精度较低,分布性较差,而NeQuick模型全球改正率分布则较为平均且平滑.      

Comparison between ionospheric peak parameters retrieved from COSMIC measurement and ionosonde observation over Sanya

In this paper we compared the ionospheric peak parameters (peak electron density of the F region, NmF2, and peak height of the F region, hmF2) retrieved from the FORMASAT-3/COSMIC (COSMIC for short) satellite measurement with those from ionosonde observation at Sanya (18.3°N, 109.6°E) during the period of 2008–2013. Although COSMIC NmF2 (hmF2) tends to be lower (higher) than ionosonde NmF2 (hmF2), the results show that the ionospheric peak parameters retrieved from COSMIC measurement generally agree well with ionosonde observation. For NmF2 the correlation between the COSMIC measurement and the ionosonde observation is higher than 0.89, and for hmF2 the correlation is higher than 0.80. The correlation of the ionospheric peak parameters decreases when solar activity increases. The performance of COSMIC measurement is acceptable under geomagnetic disturbed condition. The correlation of NmF2 between COSMIC and ionosonde measurements is higher (lower) during the nighttime (daytime), while the correlation of hmF2 is lower (higher) during the nighttime (daytime).     

Examining the use of the NeQuick bottomside and topside parameterizations at high latitudes

An examination of the high latitude performance of the bottomside and topside F-layer parameterizations of the NeQuick electron density model is presented using measurements from high latitude ionosonde and Incoherent Scatter Radar (ISR) facilities.For the bottomside, we present a comparison between modeled and measured B2Bot thickness parameter. In this comparison, it is seen that the use of the NeQuick parameterization at high latitudes results in significantly underestimated bottomside thicknesses, regularly exceeding 50%. We show that these errors can be attributed to two main issues in the NeQuick parameterization:(1) through the relationship relating foF2 and M3000F2 to the maximum derivative of F2 electron density, which is used to derive the bottomside thickness, and (2) through a fundamental inability of a constant thickness parameter, semi-Epstein shape function to fit the curvature of the high latitude F-region electron density profile.For the topside, a comparison is undertaken between the NeQuick topside thickness parameterization, using measured and CCIR-modeled ionospheric parameters, and that derived from fitting the NeQuick topside function to Incoherent Scatter Radar-measured topside electron density profiles. Through this comparison, we show that using CCIR-derived foF2 and M3000F2, used in both the NeQuick and IRI, results in significantly underestimated topside thickness during summer periods, overestimated thickness during winter periods, and an overall tendency to underestimate diurnal, seasonal, and solar cycle variability. These issues see no improvement through the use of measured foF2 and M(3000)F2 values. Such measured parameters result in a tendency for the parametrization to produce a declining trend in topside thickness with increasing solar activity, to produce damped seasonal variations, and to produce significantly overestimated topside thickness during winter periods.     

Ionospheric response to the 17 March and 22 June 2015 geomagnetic storms over Wuhan region using GNSS-based tomographic technique

By using the data of GNSS (Global Navigation Satellite System) observation from Crustal Movement Observation Network of China (CMONOC), ionospheric electron density (IED) distributions reconstructed by using computerized ionospheric tomography (CIT) technique are used to investigate the ionospheric storm effects over Wuhan region during 17 March and 22 June 2015 geomagnetic storm periods. F-region critical frequency (foF2) at Wuhan ionosonde station shows an obvious decrease during recovery phase of the St. Patrick’s Day geomagnetic storm. Moreover, tomographic results present that the decrease in electron density begins at 12:00 UT on 17 March during the storm main phase. Also, foF2 shows a long-lasting negative storm effect during the recovery phase of the 22 June 2015 geomagnetic storm. Electron density chromatography presents the evident decrease during the storm day in accordance with the ionosonde observation. These ionospheric negative storm effects are probably associated with changes of chemical composition, PPEF and DDEF from high latitudes.     

A new method for improving the performance of an ionospheric model developed by multi-instrument measurements based on artificial neural network

There are remarkable ionospheric discrepancies between space-borne (COSMIC) measurements and ground-based (ionosonde) observations, the discrepancies could decrease the accuracies of the ionospheric model developed by multi-source data seriously. To reduce the discrepancies between two observational systems, the peak frequency (foF2) and peak height (hmF2) derived from the COSMIC and ionosonde data are used to develop the ionospheric models by an artificial neural network (ANN) method, respectively. The averaged root-mean-square errors (RMSEs) of COSPF (COSMIC peak frequency model), COSPH (COSMIC peak height model), IONOPF (Ionosonde peak frequency model) and IONOPH (Ionosonde peak height model) are 0.58 MHz, 19.59 km, 0.92 MHz and 23.40 km, respectively. The results indicate that the discrepancies between these models are dependent on universal time, geographic latitude and seasons. The peak frequencies measured by COSMIC are generally larger than ionosonde’s observations in the nighttime or middle-latitudes with the amplitude of lower than 25%, while the averaged peak height derived from COSMIC is smaller than ionosonde’s data in the polar regions. The differences between ANN-based maps and references show that the discrepancies between two ionospheric detecting techniques are proportional to the intensity of solar radiation. Besides, a new method based on the ANN technique is proposed to reduce the discrepancies for improving ionospheric models developed by multiple measurements, the results indicate that the RMSEs of ANN models optimized by the method are 14–25% lower than the models without the application of the method. Furthermore, the ionospheric model built by the multiple measurements with the application of the method is more powerful in capturing the ionospheric dynamic physics features, such as equatorial ionization, Weddell Sea, mid-latitude summer nighttime and winter anomalies. In conclusion, the new method is significant in improving the accuracy and physical characteristics of an ionospheric model based on multi-source observations.     

Collocated ionosonde and dense GPS/GLONASS network measurements of midlatitude MSTIDs

To analyze midlatitude medium-scale travelling ionospheric disturbances (MSTIDs) over Kazan (55.5°N, 49°E), Russia, the sufficiently dense network of GNSS receivers (more than 150 ground-based stations) were used. For the first time, daytime MSTIDs in the form of their main signature (band structure) on high-resolution two-dimensional maps of the total electron content perturbation (TEC maps) are compared with ionosonde data with a high temporal resolution. For a pair of events, a relationship between southwestward TEC perturbations and evolution of F2 layer traces was established. So F2 peak frequency varied in antiphase to TEC perturbations. The ionograms show that during the movement of plasma depletion band (overhead ionosonde) the F2 peak frequency is the highest, and vice versa, for the plasma enhancement band, the F2 peak frequency is the lowest. One possible explanation may be a greater inclination of the radio beam from the vertical during the placement of a plasma enhancement band above the ionosonde, as evidenced by the absence of multiple reflections and the increased occurrence rate of additional cusp trace. Another possible explanation may be the redistribution of the electron content in the topside ionosphere with a small decrease in the F peak concentration of the layer with a small increase in TEC along the line-of-sight. Analysis of F2 peak frequency variation shows that observed peak-to-peak values of TEC perturbation equal to 0.4 and 1 TECU correspond to the values of ΔN/N equal to 13% and 28%. The need for further research is evident.     

Trans-ionospheric GPS signal delay gradients observed over mid-latitude Europe during the geomagnetic storms of October–November 2003

Ionospheric disturbances are known to have adverse effects on the satellite-based communication and navigation. One particular type of ionospheric effects, observed during major geomagnetic storms and threatening the integrity performance of both ground-based and space-based GNSS augmentation systems, is the sharp increase/decrease in the ionospheric delay that propagates in horizontal direction, thus called for convenience ‘moving ionospheric wall’. This paper presents preliminary results from researching such anomalous ionospheric delay gradients at European middle latitudes during the storm events of 29 October 2003 and 20 November 2003. For the purpose, 30-s GPS data from the Belgian permanent network was used for calculating and analysing the slant ionospheric delay and total electron content values. It has been found that, during these two particular storm events, substantial gradients did occur in Europe although they were not so pronounced as in the American sector.     

GNSS single frequency ionospheric range delay corrections: NeQuick data ingestion technique

A study of the performance of the NeQuick model and the Klobuchar model for GNSS single frequency range delay correction on a global scale was done using data for moderate solar activity. In this study NeQuick was used in the way intended for Galileo. This study is to assess the performance of the two models at each ionospheric geographic region during moderate solar activity as previously published studies were concentrated only on high solar activity. The results obtained showed that NeQuick outperformed Klobuchar for the whole year at the three geographical regions of the ionosphere. In terms of monthly root mean square of mismodeling, NeQuick outperformed Klobuchar by 15 TECU or more at low-latitudes, 5 TEC or more at mid-latitudes, and 1 TECU or more at high-latitudes.     

NeQuick G model based scale factor determination for using SBAS ionosphere corrections at low earth orbit

A space-based augmentation system (SBAS) provides real-time correction data for global navigation satellite system (GNSS) users near ground. In order to use the SBAS ionosphere correction for low Earth orbit (LEO) satellites, the correction should be scaled down for the LEO altitude. This scale factor varies with ionosphere distribution and it is hard to determine the value at LEO in real time. We propose a real-time scale factor determination method by using Galileo GNSS’s NeQuick G model. A LEO satellite GPS data and SBAS data received on ground were used to evaluate the performance of the NeQuick G derived variable scale factor. The NeQuick G derived scale factor shows a significant accuracy improvement over NeQuick G model or pre-determined constant scale factor. It improves a vertical positioning accuracy of the LEO satellite. The error mean reductions of the vertical positioning over NeQuick G and the constant scale factor are 31.5% and 11.7%, respectively.     

GPS derived TEC and foF2 variability at an equatorial station and the performance of IRI-model

The ionosphere induces a time delay in transionospheric radio signals such as the Global Positioning System (GPS) signal. The Total Electron Content (TEC) is a key parameter in the mitigation of ionospheric effects on transionospheric signals. The delay in GPS signal induced by the ionosphere is proportional to TEC along the path from the GPS satellite to a receiver. The diurnal monthly and seasonal variations of ionospheric electron content were studied during the year 2010, a year of extreme solar minimum (F10.7 = 81 solar flux unit), with data from the GPS receiver and the Digisonde Portable Sounder (DPS) collocated at Ilorin (Geog. Lat. 8.50°N, Long. 4.50°E, dip −7.9°). The diurnal monthly variation shows steady increases in TEC and F2-layer critical frequency (foF2) from pre-dawn minimum to afternoon maximum and then decreases after sunset. TEC show significant seasonal variation during the daytime between 0900 and 1900 UT (LT = UT + 1 h) with a maximum during the March equinox (about 35 TECU) and minimum during the June solstice (about 24 TECU). The GPS-TEC and foF2 values reveal a weak seasonal anomaly and equinoctial asymmetry during the daytime. The variations observed find their explanations in the amount of solar radiation and neutral gas composition. The measured TEC and foF2 values were compared with last two versions of the International Reference Ionosphere (IRI-2007 and IRI-2012) model predictions using the NeQuick and CCIR (International Radio Consultative Committee) options respectively in the model. In general, the two models give foF2 close to the experimental values, whereas significant discrepancies are found in the predictions of TEC from the models especially during the daytime. The error in height dependent thickness parameter, daytime underestimation of equatorial drift and contributions of electrons from altitudes above 2000 km have been suggested as the possible causes.     

Ionospheric parameters in the European sector during the magnetic storm of August 25–26, 2018

Variations of ionospheric parameters Total Electron Content (TEC) by GNSS, critical frequency (foF2) by vertical sounding and electron density (Ne) by low-altitude satellite were studied at high, mid and low latitudes of the European sector during the magnetic storm of August 25–26, 2018. During the main phase of the storm the ionospheric F2-layer was under the positive disturbance at mid and low latitudes. Then the transition from the positive to negative ΔfoF2 values occurred at all latitudes. The recovery phase was characterized by negative ionospheric disturbance at all latitudes. This is due to the decrease of thermospheric O/N2 ratio during the recovery phase of the storm. The intense Es layers screened the reflections from the F2-layer on August 26th at high and at low latitudes but at different times. Some blackouts occurred due to the high absorption level at high latitudes. In general, foF2 and TEC data were highly correlated. The major Ne changes were at the low latitudes. In general, Ne data confirmed the ionospheric dynamics revealed with foF2 and TEC.     

The ionospheric response to high-intensity long duration continuous AE activity (HILDCAA) event (13–15 April 2005) over mid-latitude African region

The ionospheric responses to High-Intensity Long Duration Continuous Auroral Electrojet Activity (HILDCAA) event which happened following the CIR-driven storm were studied over the southern hemisphere mid-latitude in the African sector. The 13–15 April 2005 event was analysed to understand some of the mechanisms responsible for the ionospheric changes during HILDCAA event. The ionosonde critical frequency of F2 layer (foF2) and Global Navigation Satellite System (GNSS) Total Electron Content (TEC) were used to analyse the ionospheric responses. The daytime increase in foF2 and TEC values were observed on 13 April 2005. The TEC and foF2 enhancement could be attributed to Large Scale Traveling Ionospheric Disturbances (LSTIDs), increase in thermospheric neutral composition changes, Prompt Penetration Electric Field (PPEF) and an expansion of Equatorial Ionization Anomaly (EIA) to the mid-latitude.     

Using TensorFlow-based Neural Network to estimate GNSS single frequency ionospheric delay (IONONet)

In the last 20?years, and in particular in the last decade, the availability of propagation data for GNSS has increased substantially. In this sense, the ionosphere has been sounded with a large number of receivers that provide an enormous amount of ionospheric data. Moreover, the maturity of the models has also been increased in the same period of time. As an example, IGS has ionospheric maps from GNSS data back to 1998, which would allow for the correlation of these data with other quantities relevant for the user and space weather (such as Solar Flux and Kp). These large datasets would account for almost half a billion points to be analyzed. With the advent and explosion of Big Data algorithms to analyze large databases and find correlations with different kinds of data, and the availability of open source code libraries (for example, the TensorFlow libraries from Google that are used in this paper), the possibility of merging these two worlds has been widely opened. In this paper, a proof of concept for a single frequency correction algorithm based in GNSS GIM vTEC and Fully Connected Neural Networks is provided. Different Neural Network architectures have been tested, including shallow (one hidden layer) and deep (up to five hidden layers) Neural Network models. The error in training data of such models ranges from 50% to 1% depending on the architecture used. Moreover, it is shown that by adjusting a Neural Network with data from 2005 to 2009 but tested with data from 2016 to 2017, Neural Network models could be suitable for the forecast of vTEC for single frequency users. The results indicate that this kind of model can be used in combination with the Galileo Signal-in-Space (SiS) NeQuick G parameters. This combination provides a broadcast model with equivalent performances to NeQuick G and better than GPS ICA for the years 2016 and 2017, showing a 3D position Root Mean Squared (RMS) error of approximately 2?m.     

Diagnostics of equatorial and low latitude ionosphere by TEC mapping over Brazil

The total electron content (TEC) in the equatorial and low-latitude ionosphere over Brazil was monitored in two dimensions by using 2011 data from the ground-based global navigation satellite system (GNSS) receiver network operated by the Brazilian Institute for Geography and Statistics. It was possible to monitor the spatial and temporal variations in TEC over Brazil continuously during both day and night with a temporal interval of 10 min and a spatial resolution of about 400 km. The daytime equatorial ionization anomaly (EIA) and post-sunset plasma enhancement (PS-EIA) were monitored over an area corresponding to a longitudinal extension of 4000 km in South America. Considerable day-to-day variation was observed in EIA and PS-EIA. A large latitudinal and longitudinal gradient of TEC indicated a significant ionospheric range error in application of the GNSS positioning system. Large-scale plasma bubbles after sunset were also mapped over a wide range. Depletions with longitudinally separated by more than 800 km were observed. They were extended by more than 2000 km along the magnetic field lines and drifted eastward. It is expected that 2-dimensional TEC mapping can serve as a useful tool for diagnosing ionospheric weather, such as temporal and spatial variation in the equatorial plasma trough and crest, and particularly for monitoring the dynamics of plasma bubbles.     

Local ionospheric electron density profile reconstruction in real time from simultaneous ground-based GNSS and ionosonde measurements

The purpose of the LIEDR (local ionospheric electron density profile reconstruction) system is to acquire and process data from simultaneous ground-based total electron content (TEC) and digital ionosonde measurements, and subsequently to deduce the vertical electron density distribution above the ionosonde’s location. LIEDR is primarily designed to operate in real time for service applications and, for research applications and further development of the system, in a post-processing mode. The system is suitable for use at sites where collocated TEC and digital ionosonde measurements are available. Developments, implementations, and some preliminary results are presented and discussed in view of possible applications.