Global navigation satellite systems (GNSS) provide a great data source about the ionosphere state. These data can be used for testing ionosphere models. We studied the performance of nine ionospheric models (Klobuchar, NeQuickG, BDGIM, GLONASS, IRI-2016, IRI-2012, IRI-Plas, NeQuick2, and GEMTEC) both in the total electron content (TEC) domain—i.e., how precise the models calculate TEC—and in the positioning error domain—i.e., how the models improve single frequency positioning. The whole data set covers 20 years (2000–2020) from 13 GNSS stations, but the main analysis involves data during 2014–2020 when calculations are available from all the models. We used single-frequency positioning without ionospheric correction and with correction via global ionospheric maps (IGSG) data as expected limits for errors. Improvements against noncorrected solution were as follows: GIM IGSG—22.0%, BDGIM—15.3%, NeQuick2—13.8%, GEMTEC, NeQuickG and IRI-2016—13.3%, Klobuchar—13.2%, IRI-2012—11.6%, IRI-Plas—8.0%, GLONASS—7.3%. TEC bias and mean absolute TEC errors for the models are as follows: GEMTEC—−0.3 and 2.4 TECU, BDGIM—−0.7 and 2.9 TECU, NeQuick2—−1.2 and 3.5 TECU, IRI-2012—−1.5 and 3.2 TECU, NeQuickG—−1.5 and 3.5 TECU, IRI-2016—−1.8 and 3.2 TECU, Klobuchar—1.2 and 4.9 TECU, GLONASS—−1.9 and 4.8 TECU, and IRI-Plas—3.1 and 4.2 TECU. While TEC and positioning domains differ, new-generation operational models (BDGIM and NeQuickG) could overperform or at least be at the same level as classical empirical models.

Recommendations of the International Reference Ionosphere (IRI) Workshop 2017 in Taoyuan City, Taiwan and International GNSS Service (IGS) Workshop 2018 in Wuhan, China included establishment of an ionosphere mapping service that would fuse measurements from two independent sensor networks: IGS permanent GNSS receivers providing the vertical total electron content (VTEC) measurements and ionosondes of the Global Ionosphere Radio Observatory (GIRO) that compute the bottomside vertical profiles of the ionospheric plasma density. Using available GAMBIT software at GIRO, we introduced new VTEC products to its data roster: previously unavailable global average (climate) maps of VTEC and slab thickness based on climatological capabilities of IRI. Incorporation of the VTEC and τ maps into the GAMBIT Explorer environment provided data analysts with nearly 10-year history of the reference average VTEC records and opened access to the GAMBIT toolkit for evaluation and validation of the τ computations. This result is the first step towards establishing an infrastructure and the data workflow to provide GAMBIT users with the low latency and consistent quality and usability of the ionospheric weather-climate specifications. Combination of IGS-provided VTEC and GIRO-provided peak density of F2 layer NmF2 allows ground-based evaluation of the equivalent slab thickness τ, a derived property of the near-Earth plasma that characterizes the skewness of its vertical profile up to the GNSS spacecraft altitudes.

At equatorial and high latitudes, the intense ionospheric irregularities and plasma density gradients can seriously affect the performances of radio communication and satellite-based navigation systems; that represents a challenging topic for the scientific and engineering communities and operational use of communication and navigation services. The GNSS-based ROTI (rate of TEC index) is one of the most widely used indices to monitor the occurrence and intensity of ionospheric irregularities. In this paper, we examined the long-term performance of the ROTI in terms of finding the climatological characteristics of TEC fluctuations. We considered the different scale temporal signatures and checked the general sensitivity to the solar and geomagnetic activity. We retrieved and analyzed long-term time-series of ROTI values for two chains of GNSS stations located in European and North-American regions. This analysis covers three full years of the 24th solar cycle, which represent different levels of solar activity and include periods of intense geomagnetic storms. The ionospheric irregularities’ geographical distribution, as derived from ROTI, shows a reasonable consistency to be found within the poleward/equatorward boundaries of the auroral oval specified by empirical models. During magnetic midnight and quiet-time conditions, the equatorward boundary of the ROTI-derived ionospheric irregularity zone was observed at 65–70${}^{°}$ of north magnetic latitude, while for local noon conditions this boundary was more poleward at 75–85 magnetic latitude. The ionospheric irregularities of low-to-moderate intensity were found to occur within the auroral oval at all levels of geomagnetic activity and seasons. At moderate and high levels of solar activity, the intensities of ionospheric irregularities are larger during local winter conditions than for the local summer and polar day conditions. We found that ROTI displays a selective latitudinal sensitivity to the auroral electrojet activity—the strongest dependence (correlation R > 0.6–0.8) was observed within a narrow latitudinal range of 55–70${}^{°}$N magnetic latitude, which corresponded to a band of the largest ROTI values within the auroral oval zone expanded equatorward during geomagnetic disturbances.