The Japanese Quasi-Zenith Satellite System (QZSS) is a regional navigation satellite system covering the entire Asia-Oceania region. Except for the standard satellite navigation signals, QZSS satellites also broadcast L6E augmentation signals with real time GNSS precise orbit every 30 s and clock messages every 1 s, which is very important and necessary for Real-Time precise point positioning (RTPPP) applications. In this paper, the MADOCA real-time services derived from L6E augmentation signals were evaluated for both accuracy and availability compared with IGS final products. To avoid the datum difference of GPS orbit between MADOCA real-time and IGS final products, the 7-parameters Helmert transformation was firstly used in this paper, and then the orbit was evaluated on radial, along, and cross-track directions. On the clock evaluation, the mean satellites clock errors were taken as reference clock error, and then the standard deviation (STD) was calculated for each satellite. Furthermore, the signal in space range errors (SISRE) were also summarized to evaluate the ranging-measurement accuracy. Seven-day evaluation results show that satellite orbit, clock errors, and the final SISRE errors range as being 1.8–3.9 cm, 0.04–0.15 ns (1.2–4.5 cm), and 5–10 cm, respectively. For the one-year long-term evaluation, daily SISRE errors in 2018 show consistent performance with that of seven days. Furthermore, the open source software RTKLIB was used to evaluate the kinematic PPP performance based on the MADOCA real-time products, and it shows that the daily positioning accuracy of the 20 globally distributed IGS stations can reach 4.9, 4.2, 11.7, and 12.1 cm in the east, north, up, and 3D directions, respectively. Hence, it is concluded that the current MADOCA real-time ephemeris products can provide orbit and clock products with SISRE on centimeters level with high interval, which could meet the demands of the RTPPP solution and serve real-time users who can access the MADOCA real-time products via L6E signal or internet.

Ionospheric delay is one of the most problematic errors in satellite-based positioning data processing. The Global Ionospheric Map (GIM), which is publicly available daily in various analysis centers, is thus vitally important for positioning users. There are variations in the accuracy and consistency of GIMs issued by Ionosphere Associate Analysis Centers (IAACs) due to the differences in ionospheric modeling methods and selected tracking stations. In this study on the International GNSS Service’s (IGS) final GIM, the ionospheric total electron content (TEC) (from 243 global navigation satellite system (GNSS) monitoring stations around the world) and the ionospheric TEC (from the Jason-3 altimeter satellite) are selected as reference. By using these three references, we evaluate the performance and consistency of final GIM products from seven IGS IAACs, including the Chinese Academy of Sciences (CAS), the Center for Orbit Determination in Europe (CODE), Natural Resources Canada (EMR), the European Space Agency (ESA), the Jet Propulsion Laboratory (JPL), Universitat Politècnica de Catalunya (UPC), and Wuhan University (WHU) in the mid-solar activity year (2022) and the low-solar activity year (2020). Firstly, the comparison with the IGS final GIM shows that the consistency of each GIMs is basically the same, with the mean value ranging from −0.3 TECu (total electron content unit) to 1.4 TECu. Secondly, the validation with Jason-3 altimeter satellite shows that the accuracy of several GIMs is almost the same, except for the JPL with the worst accuracy and an overall mean deviation (BIAS) between 2 and 6 TECu. Thirdly, the comparison with VTEC extracted from GNSS monitor stations shows that the CAS has the best accuracy in different latitude bands with a root mean square (RMS) of about 2.2–4.7 TECu. In addition, it is found that the accuracy in areas with more stations for ionospheric modelling is better than those with less stations in different latitude bands; meanwhile, the accuracy is closely related to the modeling methods of different GIMs.

As global navigation satellite system (GNSS)stations are sparsely distributed in oceanic area, oceanic areas usually have lower precision than continental areas on a global ionosphere maps (GIM). On the other hand, space-borne observations like satellite altimetry (SA) and ionospheric radio occultation (IRO) have substantial dual-frequency observations in oceanic areas, which could be used for total electron content (TEC) retrieval. In this paper, the Jason-2 SA and Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC) IRO products were used to assess the precision of IGS GIM products. Both the systematic biases and scaling factors between the international GNSS service (IGS) GIM TEC and space-borne TEC were calculated, and the statistical results show that the biases and the scaling factors obviously vary under different temporal-spatial conditions. This analysis shows that these differences are variable with diurnal and latitude factors, that is, the differences in biases during the day time are higher than those during the night time, and larger biases are experienced at lower latitude areas than at high latitude areas. The results also show that in the southern hemisphere middle-high latitude area and some other central oceanic areas, the space-borne TEC values are even higher than GIM TEC values. As the precision of space-borne TEC should be evenly distributed around different areas on Earth, it can be explain that the TEC in these areas is undervalued by the current GIM model, and the space-borne SA and IRO techniques could be used as complementary observations to improve the accuracy and reliability of TEC values in these areas.

The broadcast ionospheric model is mainly used to correct the ionospheric delay error for single-frequency users. Since the BeiDou global ionospheric delay correction model (BDGIM) is a novel broadcast ionospheric model for BDS-3, its performance was analyzed through single point positioning (SPP) in this study. Twenty-two stations simultaneously receiving B1C, B2a, B1I and B3I signals were selected from the International GNSS Service (IGS) and the International GNSS Monitoring and Assessment System (iGMAS) tracking networks for the SPP experiments. The differential code bias (DCB) parameters were used to correct the hardware delays in the signals of B1C and B2a. The results showed that the BDGIM performs the best in high-latitude areas, and can effectively improve the positioning accuracy compared with the Klobuchar model. The average 3D positioning accuracy of the four civil signals can reach 3.58 m in high-latitude areas. The positioning accuracies with the BDGIM in the northern hemisphere are better than those in the southern hemisphere, and the global average 3D positioning accuracy of the four civil signals is 4.60 m. The performance of the BDGIM also shows some seasonal differences. The BDGIM performs better than the Klobuchar model on the days of spring equinox and winter solstice, while the opposite is true on the days of summer solstice and autumn equinox. On the day of winter solstice, the average 3D accuracies with the BDGIM on the signals of B1C, B2a, B1I and B3I are 4.13 m, 5.32 m, 4.40 m and 4.49 m, respectively. Although the SPP accuracies are to some extent affected by the geomagnetic storm, the BDGIM generally performs better and are more resistant to the geomagnetic storm than the Klobuchar model.

Ionospheric delay is one of the most problematic errors in single-frequency (SF) global navigation satellite system (GNSS) data processing. Global/regional ionospheric maps (GIM/RIM) are thus vitally important for positioning users. Given the coexistence of multi-GNSS, the integration of quad-constellation observations is essential for improving the distribution of ionospheric penetration points (IPPs) and increasing redundant observations compared with the existing GIM products from the IGS analysis center. In this paper, quad-constellation (GPS/GLONASS/Galileo/BDS) observations are applied to set up the RIM over Australia with uncombined precise point positioning (UC-PPP) and a low-order spherical harmonic function. The generated RIMs are then introduced to ionosphere-corrected (IC) and ionosphere-weighted (IW) single-frequency PPP (IC-SFPPP and IW-SFPPP) to verify their performance in terms of positioning accuracy and convergence time. Taking the CODE GIM as a reference, the results show that the mean root mean square (RMS) of VTEC differences is 0.867 TECUs, and the quad-constellation RIM (referred as ‘RIM4′) can improve the RMS of RIMs compared to single-constellation mode at the edge of regional experiment area. The application of the RIM4 in the BDS IC-SFPPP results in a 18.38% improvement (from 100.47 cm to 82.00 cm) of 3D positioning RMS compared to the CODE-GIMs, whereas 35.36% enhancement (from 115.92 cm to 74.62 cm) of 3D positioning RMS is achievable during an active ionospheric period. Moreover, if the criterion of the convergence time is defined as when positioning errors in the horizontal and vertical directions are less than 0.3 m and 0.6 m for 20 consecutive epochs, the IW-SFPPP can significantly speed up the convergence time compared to the uncombined SFPPP; that is, the convergence time is reduced by 52.7% (from 37 min to 17.5 min), 37.2% (from 72.5 min to 45.5 min), and 37.1% (from 62.0 min to 39.0 min) in the north, east and up direction, respectively, at the 68% confidence level.

The Global Navigation Satellite System (GNSS) Observations can invert the characteristics of global ionospheric changes, which help understand the coupled physical processes in the atmosphere caused by earthquakes. On 6 February 2023, Turkey experienced two strong earthquakes: an Mw7.8 event at UTC 01:17:34 and an Mw7.5 event at UTC 10:24:50. The propagation direction of ionospheric disturbances aligns with the rupture zones, which is bidirectional northeast-southwest ruptures for the first event and east-west ruptures for the second one. According to the results of Turkey stations, it is found that the Co-seismic Ionospheric Disturbances (CIDs) of the first earthquake are primarily driven by acoustic waves, but there are also Rayleigh waves and mixed waves components. Due to the influence of the first earthquake and aftershocks, the CIDs of the second earthquake are dominated by acoustic wave and Rayleigh wave. Notably, the Mw7.5 earthquake resulted in stronger ionospheric disturbance amplitudes than the Mw7.8 event, while the Mw7.8 event affected a larger spatial disturbance area. Meanwhile, it is confirmed that significant CIDs will be detected when strike-slip earthquake magnitude is large enough.