Ionospheric correction using NTCM driven by GPS Klobuchar coefficients for GNSS applications

Global Navigation Satellite Systems (GNSS) require mitigation of ionospheric propagation errors because the ionospheric range errors might be larger than tens of meters at the zenith direction. Taking advantage of the frequency-dispersive property of ionospheric refractivity, the ionospheric range errors can be mitigated in dual-frequency applications to a great extent by a linear combination of carrier phases or pseudoranges. However, single-frequency GNSS operations require additional ionospheric information to apply signal delay or range error corrections. To aid single-frequency operations, the global positioning system (GPS) broadcasts 8 coefficients as part of the navigation message to drive the ionospheric correction algorithm (ICA) also known as Klobuchar model. We presented here an ionospheric correction algorithm called Neustrelitz TEC model (NTCM) which can be used as complementary to the GPS ICA. Our investigation shows that the NTCM can be driven by Klobuchar model parameters to achieve a significantly better performance than obtained by the mother ICA algorithm. Our research, using post-processed reference total electron content (TEC) data from more than one solar cycle, shows that on average the RMS modeled TEC errors are up to 40% less for the proposed NTCM model compared to the Klobuchar model during high solar activity period, and about 10% less during low solar activity period. Such an approach does not require major technology changes for GPS users rather requires only introducing the NTCM approach a complement to the existing ICA algorithm while maintaining the simplicity of ionospheric range error mitigation with an improved model performance.     

Brazilian active GNSS networks as systems for monitoring the ionosphere

This research shows the viability of using Global Navigation Satellite System (GNSS) stations from Brazilian active networks in monitoring the ionosphere. Various indexes of ionospheric irregularities and scintillation of GNSS signals, estimated in real-time and post-processed from GNSS data, are explored for this purpose. This way, an increase in the spatial resolution of ionospheric information is provided, allowing the generation of maps of scintillation and irregularities in observing the spatial and temporal behavior of the layer’s activity cycle, since the number of ionosondes, imagers, and radars is insufficient for monitoring the irregularities in Brazil. Experiments to evaluate the estimates of the indexes are performed for periods of high and low variability of electrons. Three Brazilian networks are used: the Brazilian Network for Continuous Monitoring (RBMC), the GNSS Active Network of Sao Paulo State (GNSS-SP), and CIGALA/CALIBRA. The results are compared with data from ionosondes and PolaRxS-PRO Septentrio receivers, proving compatible with moderate to high correlations. An analysis of the seasonal variation during the peak of solar cycle 24 is carried out. The maps allow identifying the displacement of ionospheric irregularities along the magnetic equator over Brazil, from northeast to southwest, starting at 7:00 pm and ending at 2:00 am local time. Real-time monitoring is carried out for the summer solstice in the southern hemisphere, and results are consistent with those from the post-processed mode. The indexes and maps can be applied to the analysis of GNSS positioning. Real-time ionospheric information can be used in important practical applications because the displacement monitoring of irregularities allows prior knowledge of whether there will be a deterioration of positioning accuracy in a certain region.     

适用于不同尺度区域的Klobuchar-like电离层模型

导航定位中运用最广泛的电离层修正模型是Klobuchar模型,但经典的Klobuchar模型不能满足日益增长的导航定位精度的需求,因此不同的精化模型被提出。本文利用GIMs分析了夜间电离层随地方时的变化和电离层电子总含量随纬度的变化情况,在对各种适用范围较广的模型精化方案进行归纳总结的基础上,提出了一种适用于不同尺度区域的Klobuchar-like模型,并利用不同太阳活动时期不同季节的GIMs建立了适用于单站、大区域和全球的Klobuchar-like模型、14参数Klobuchar模型和8参数Klobuchar模型。Klobuchar-like模型单站、区域、全球的修正率分别达到了92.96%、91.55%、72.67%,均高于14参数、8参数Klobuchar模型和GPS Klobuchar模型,表明了该模型的有效性与实用性。     

Global ionosphere map constructed by using total electron content from ground-based GNSS receiver and FORMOSAT-3/COSMIC GPS occultation experiment

Effects of rapidly changing ionospheric weather are critical in high accuracy positioning, navigation, and communication applications. A system used to construct the global total electron content (TEC) distribution for monitoring the ionospheric weather in near-real time is needed in the modern society. Here we build the TEC map named Taiwan Ionosphere Group for Education and Research (TIGER) Global Ionospheric Map (GIM) from observations of ground-based GNSS receivers and space-based FORMOSAT-3/COSMIC (F3/C) GPS radio occultation observations using the spherical harmonic expansion and Kalman filter update formula. The TIGER GIM (TGIM) will be published in near-real time of 4-h delay with a spatial resolution of 2.5° in latitude and 5° in longitude and a high temporal resolution of every 5 min. The F3/C TEC results in an improvement on the GIM of about 15.5%, especially over the ocean areas. The TGIM highly correlates with the GIMs published by other international organizations. Therefore, the routinely published TGIM in near-real time is not only for communication, positioning, and navigation applications but also for monitoring and scientific study of ionospheric weathers, such as magnetic storms and seismo-ionospheric anomalies.     

Medium-scale traveling ionospheric disturbances observed by GPS receiver network in Japan: a short review

Medium-scale traveling ionospheric disturbances (MSTID) are wave-like perturbations of the ionospheric plasma with wavelengths of several hundred kilometres and velocities of several hundred metres per second. MSTID is one of the most common ionospheric phenomena that generally induce the perturbations of ionospheric total electron content (TEC) by ∼10¹⁶ electron/m², which corresponds to ∼54 ns (16.2 cm) delay at GPS L1 signal. In the past decade, several new characteristics on MSTIDs have been revealed by the TEC observations using the dense GPS receiver network in Japan. In this paper, we provide a short review of these new observations and summarize the morphological characteristics of MSTIDs in Japan.     

北斗三号全球导航卫星系统全球广播电离层延迟修正模型(BDGIM)应用性能评估

2020年6月23日,我国北斗三号全球导航卫星系统正式完成星座全球组网。北斗三号全球导航卫星系统采用新一代全球广播电离层延迟修正模型(BDGIM),为用户提供电离层延迟改正服务。本文利用高精度全球电离层格网(GIM)以及实测BDS/GPS数据提供的电离层TEC作为参考,从延迟改正精度及北斗单频伪距单点定位应用、模型系数性能等方面,对北斗三号系统组网前后(2020年5月1日至2020年7月20日)BDGIM模型的改正精度等应用性能进行了分析与研究,并将其与美国GPS播发的Klobuchar模型和北斗二号卫星导航系统播发的BDS Klobuchar模型进行对比。研究表明,BDGIM模型在对北斗三号系统组网完成前后电离层延迟修正精度没有发生显著变化。上述时段内,以国际GNSS服务(IGS)发布的最终GIM产品为参考,BDGIM模型在中国区域、亚太地区和全球范围内的电离层修正百分比分别达到84.45%、74.74%和64.57%;以选取的全球83个GNSS检测站BDS、GPS双频数据实测电离层TEC为参考,BDGIM在中国区域、亚太地区和全球范围内的电离层修正百分比分别为73.12%、70.18%及68.06%;当BDGIM模型应用于北斗单频伪距单点定位时,在中国区域、亚太地区和全球范围内分别实现了2.22、2.66和2.96 m的三维定位精度。     

On the effects of the ionospheric disturbances on precise point positioning at equatorial latitudes

In precise point positioning (PPP), the ionospheric delay is corrected in a first-order approximation from GPS dual-frequency observations, which should eliminate almost completely the ionosphere as a source of error. However, sudden plasma density variations can adversely affect the GPS signal, degrading accuracy and reliability of positioning techniques. The occurrence of plasma density irregularities is frequent at equatorial latitudes and is reflected in large total electron content (TEC) variations. We study the relation between large changes in the rate of TEC (ROT) and positioning errors in single-epoch PPP. At equatorial latitudes and during post-sunset hours, the estimated altitudes contain errors of several meters for a single-epoch position determination, and latitude and longitude estimates are also degraded. These results have been corroborated by the online CSRS-PPP (NRCan) program. Moreover, abrupt changes in the satellite geometry have been discarded as possible cause of such errors, suggesting an apparent relation between the occurrence of large ROT and degraded position estimates.     

The ionospheric eclipse factor method (IEFM) and its application to determining the ionospheric delay for GPS

A new method for modeling the ionospheric delay using global positioning system (GPS) data is proposed, called the ionospheric eclipse factor method (IEFM). It is based on establishing a concept referred to as the ionospheric eclipse factor (IEF) λ of the ionospheric pierce point (IPP) and the IEF’s influence factor (IFF) . The IEF can be used to make a relatively precise distinction between ionospheric daytime and nighttime, whereas the IFF is advantageous for describing the IEF’s variations with day, month, season and year, associated with seasonal variations of total electron content (TEC) of the ionosphere. By combining λ and with the local time t of IPP, the IEFM has the ability to precisely distinguish between ionospheric daytime and nighttime, as well as efficiently combine them during different seasons or months over a year at the IPP. The IEFM-based ionospheric delay estimates are validated by combining an absolute positioning mode with several ionospheric delay correction models or algorithms, using GPS data at an international Global Navigation Satellite System (GNSS) service (IGS) station (WTZR). Our results indicate that the IEFM may further improve ionospheric delay modeling using GPS data.     

Ionospheric tomography using GNSS: multiplicative algebraic reconstruction technique applied to the area of Brazil

Experimental analysis was performed using multiplicative algebraic reconstruction technique (MART) to map the ionosphere over Brazil. Code and phase observations from the global navigation satellite system (GNSS) together with the international reference ionosphere (IRI) enabled the estimation of ionospheric profiles and total electron content (TEC) over the entire region. Twenty-four days of data collected from existing ground-based GNSS receivers during the recent solar maximum period were used to analyze the performance of the MART algorithm. The results were compared with four ionosondes. It was demonstrated that MART estimated the electron density peak with the same degree of accuracy as the IRI model in regions with appropriate geometrical coverage by GNSS receivers for tomographic reconstruction. In addition, the slant TEC, as estimated with MART, presented lower root-mean-square error than the TEC calculated by ionospheric maps available from the International GNSS Service (IGS). Furthermore, the daily variations of the ionosphere were better represented with the algebraic techniques, compared to the IRI model and IGS maps, enabling a correlation of the elevation of the ionosphere at higher altitudes with the equatorial ionization anomaly intensification. The tomographic representations also enabled the detection of high vertical gradients at the same instants in which ionospheric irregularities were evident.     

Application of SWACI products as ionospheric correction for single-point positioning: a comparative study

In Global Navigation Satellite Systems (GNSS) using L-band frequencies, the ionosphere causes signal delays that correspond with link related range errors of up to 100 m. In a first order approximation the range error is proportional to the total electron content (TEC) of the ionosphere. Whereas this first order range error can be corrected in dual-frequency measurements by a linear combination of carrier phase- or code-ranges of both frequencies, single-frequency users need additional information to mitigate the ionospheric error. This information can be provided by TEC maps deduced from corresponding GNSS measurements or by ionospheric models. In this paper we discuss and compare different ionospheric correction methods for single-frequency users. The focus is on the comparison of the positioning quality using dual-frequency measurements, the Klobuchar model, the NeQuick model, the IGS TEC maps, the Neustrelitz TEC Model (NTCM-GL) and the reconstructed NTCM-GL TEC maps both provided via the ionosphere data service SWACI (http://swaciweb.dlr.de) in near real-time. For that purpose, data from different locations covering several days in 2011 and 2012 are investigated, including periods of quiet and disturbed ionospheric conditions. In applying the NTCM-GL based corrections instead of the Klobuchar model, positioning accuracy improvements up to several meters have been found for the European region in dependence on the ionospheric conditions. Further in mid- and low-latitudes the NTCM-GL model provides results comparable to NeQuick during the considered time periods. Moreover, in regions with a dense GNSS ground station network the reconstructed NTCM-GL TEC maps are partly at the same level as the final IGS TEC maps.     

A preliminary evaluation of the performance of multiple ionospheric models in low- and mid-latitude regions of China in 2010–2011

Ionospheric delay is a dominant error source in Global Navigation Satellite System (GNSS). Single-frequency GNSS applications require ionospheric correction of signal delay caused by the charged particles in the earth’s ionosphere. The Chinese Beidou system is developing its own ionospheric model for single-frequency users. The number of single-frequency GNSS users and applications is expected to grow fast in the next years in China. Thus, developing an appropriate ionospheric model is crucially important for the Chinese Beidou system and worldwide single-frequency Beidou users. We study the performance of five globally accessible ionospheric models Global Ionospheric Map (GIM), International Reference Ionosphere (IRI), Parameterized Ionospheric Model (PIM), Klobuchar and NeQuick in low- and mid-latitude regions of China under mid-solar activity condition. Generally, all ionospheric models can reproduce the trend of diurnal ionosphere variations. It is found that all the models have better performances in mid-latitude than in low-latitude regions. When all the models are compared to the observed total electron content (TEC) data derived from GIM model, the IRI model (2012 version) has the best agreement with GIM model and the NeQuick has the poorest agreement. The RMS errors of the IRI model using the GIM TEC as reference truth are about 3.0–10.0 TECU in low-latitude regions and 3.0–8.0 TECU in mid-latitude regions, as observed during a period of 1 year with medium level of solar activity. When all the ionospheric models are ingested into single-frequency precise point positioning (PPP) to correct the ionospheric delays in GPS observations, the PIM model performs the best in both low and mid-latitudes in China. In mid-latitude, the daily single-frequency PPP accuracy using PIM model is ~10 cm in horizontal and ~20 cm in up direction. At low-latitude regions, the PPP error using PIM model is 10–20 cm in north, 30–40 cm in east and ~60 cm in up component. The single-frequency PPP solutions indicate that NeQuick model has the lowest accuracy among all the models in both low- and mid-latitude regions of China. This study suggests that the PIM model may be considered for single-frequency GNSS users in China to achieve a good positioning accuracy in both low- and mid-latitude regions.     

海口地区扩展F和GPS电离层闪烁发生时间相关性分析

对2012年8月至2013年7月太阳活动高年海口地区扩展F和GPS L1频段电离层闪烁进行了发生时间相关性分析。结果表明,该地区的GPS电离层闪烁与扩展F具有较好的相关性,二者发生的相关系数为095.受太阳活动高年低纬电离层不均匀体发展演化特性和设备观测方式的影响,观测到的电离层闪烁起始时间稍早于扩展F。扩展F结束时间比电离层闪烁结束时间有所滞后,分析认为,在低纬电离层不均匀体的消亡阶段,1 km以下中小尺度的不均匀体首先消失。      

RINEX_HO: second- and third-order ionospheric corrections for RINEX observation files

When GNSS receivers capable of collecting dual-frequency data are available, it is possible to eliminate the first-order ionospheric effect in the data processing through the ionosphere-free linear combination. However, the second- and third-order ionospheric effects still remain. The first-, second- and third-order ionospheric effects are directly proportional to the total electron content (TEC), although the second- and third-order effects are influenced, respectively, by the geomagnetic field and the maximum electron density. In recent years, the international scientific community has given more attention to these kinds of effects and some works have shown that for high precision GNSS positioning these effects have to be taken into consideration. We present a software tool called RINEX_HO that was developed to correct GPS observables for second- and third-order ionosphere effects. RINEX_HO requires as input a RINEX observation file, then computes the second- and third-order ionospheric effects, and applies the corrections to the original GPS observables, creating a corrected RINEX file. The mathematical models implemented to compute these effects are presented, as well as the transformations involving the earth’s magnetic field. The use of TEC from global ionospheric maps and TEC calculated from raw pseudorange measurements or pseudoranges smoothed by phase is also investigated.     

Temporal and spatial variations of global ionospheric total electron content under various solar conditions

By utilizing the numerical technique of principal component analysis (PCA), this work analyses temporal and spatial variations of the ionosphere under various solar conditions during the period 1999–2013. Applying the PCA technique to the time series of the global ionospheric total electron content (TEC) maps provides an efficient method for analyzing the main ionospheric variability on a global scale that is able to decompose periodic variations (e.g., annual and semiannual oscillations) while retaining the asymmetry in the temporal and spatial domains (e.g., seasonal and equator anomalies). The TEC series of different local times are processed separately at two time scales: (1) the whole 15 years of the period of study and (2) the individual years. In contrast with previous studies, the analysis of the dataset of the 15 years shows that dawn (e.g., LT4–6) and late morning (LT10–12) are the more remarkable characteristic times for ionospheric variability. This study also reveals a cyclic trend of the variability with respect to local times. The first two modes, which contain 80–90% of the total variance, represent spatial distributions and temporal variations with respect to the different stages of the solar cycle and local times. Annual and semiannual variations are demodulated from the first two modes, and the results show that these variations evidently have distinct trends for daytime and nighttime. An exception is that, under active solar conditions, extremely strong solar irradiance during the daytime has a residual effect on the variability of the nighttime.     

Two-step method for the determination of the differential code biases of COMPASS satellites

The differential code bias (DCB) in satellites of the Global Navigation Satellite Systems (GNSS) should be precisely corrected when designing certain applications, such as ionospheric remote sensing, precise point positioning, and time transfer. In the case of COMPASS system, the data used for estimating DCB are currently only available from a very limited number of global monitoring stations. However, the current GPS/GLONASS satellite DCB estimation methods generally require a large amount of geographically well-distributed data for modeling the global ionospheric vertical total electron content (TEC) and are not particularly suitable for current COMPASS use. Moreover, some satellites with unstable DCB (i.e., relatively large scatter) may affect other satellite DCB estimates through the zero-mean reference that is currently imposed on all satellites. In order to overcome the inadequacy of data sources and to reduce the impact of unstable DCB, a new approach, designated IGGDCB, is developed for COMPASS satellite DCB determination. IGG stands for the Institute of Geodesy and Geophysics, which is located in Wuhan, China. In IGGDCB, the ionospheric vertical TEC of each individual station is independently modeled by a generalized triangular series function, and the satellite DCB reference is selected using an iterative DCB elimination process. By comparing GPS satellite DCB estimates calculated by the IGGDCB approach based on only a handful (e.g., seven) of tracking stations against that calculated by the currently existing methods based on hundreds of tracking stations, we are able to demonstrate that the accuracies of the IGGDCB-based DCB estimates perform at the level of about 0.13 and 0.10?ns during periods of high (2001) and low (2009) solar activity, respectively. The iterative method for DCB reference selection is verified by statistical tests that take into account the day-to-day scatter and the duration that the satellites have spent in orbit. The results show that the impact of satellites with unstable DCB can be considerably reduced using the IGGDCB method. It is also confirmed that IGGDCB is not only specifically valid for COMPASS but also for all other GNSS.     

不同全球电离层格网产品在中国区域的应用精度评估与分析

基于全球卫星导航系统(GNSS)跟踪站和Jason测高卫星获得的电离层总电子含量(TEC),系统评估了国内外5家国际GNSS服务(IGS)电离层分析中心以及全球连续监测评估系统(iGMAS)综合中心的预报、快速和最终全球电离层格网产品(GIM)在中国区域的精度和标准单点定位应用性能. 结果表明:不同类型GIM产品在中国区域的精度由高到低分别是最终、快速和预报GIM产品;在太阳活动水平较低时,不同GIM产品精度大致相当;在太阳活动水平较高时,西班牙加泰罗尼亚理工大学(UPC)和iGMAS的快速和最终GIM产品精度优于其他机构同类型产品.      

Eye on the Ionosphere: The Correlation between Solar 10.7 cm Radio Flux and Ionospheric Range Delay

Patricia Doherty joins the regular contributors of this column to discuss the correlation between measurements of solar 10.7 cm radio flux and ionospheric range delay effects on GPS. Mrs. Doherty has extensive experience in the analysis of ionospheric range delays from worldwide systems and in the utilization and development of analytical and theoretical models of the Earth's ionosphere. Ionospheric range delay effects on GPS and other satellite ranging systems are directly proportional to the Total Electron Content (TEC) encountered along slant paths from a satellite to a ground location. TEC is a highly variable and complex parameer that is a function of geographic location, local time, season, geomagnetic activity, and solar activity. When insufficiently accounted for, ionospheric TEC can seriously limit the performance of satellite ranging applications. Since the ionosphere is a dispersive medium, dual-frequency Global Positoning System (GPS) users can make automatic corrections for ionospheric range delay by computing the apparent difference in the time delays between the two signals. Single-frequency GPS users must depend on alternate methods to account for the ionospheric range delay. Various models of the ionosphere have been used to provide estimates of ionospheric range delay. These models range from the GPS system's simple eight-coefficient algorithm designed to correct for approximately 50% rms of the TEC, to state-of-the-art models derived from physical first principles, which can correct for up to 70 to 80% rms of the TEC but at a much greater computational cost. In an effort to improve corrections for the day-to-day variability of the ionosphere, some attempts have been made to predict the TEC by using the daily values of solar 10.7 cm radio flux (F_10,7). The purpose of this article is to show that this type of prediction is not useful due to irregular, and sometimes very poor, correlation between daily values of TEC and F_10.7. Long-term measurements of solar radio flux, however, have been shown to be well correlated with monthly mean TEC, as well as with the critical frequency of the inonospheric F2 region (foF2), which is proportional to the electron density at the peak of the ionospheric F2 region. ? 2000 John Wiley & Sons, Inc.     

The ionosphere: effects,GPS modeling and the benefits for space geodetic techniques

The main goal of this paper is to provide a summary of our current knowledge of the ionosphere as it relates to space geodetic techniques, especially the most informative technology, global navigation satellite systems (GNSS), specifically the fully deployed and operational global positioning system (GPS). As such, the main relevant modeling points are discussed, and the corresponding results of ionospheric monitoring are related, which were mostly computed using GPS data and based on the direct experience of the authors. We address various phenomena such as horizontal and vertical ionospheric morphology in quiet conditions, traveling ionospheric disturbances, solar flares, ionospheric storms and scintillation. Finally, we also tackle the question of how improved knowledge of ionospheric conditions, especially in terms of an accurate understanding of the distribution of free electrons, can improve space geodetic techniques at different levels, such as higher-order ionospheric effects, precise GNSS navigation, single-antenna GNSS orientation and real-time GNSS meteorology.     

Tridimensional reconstruction of the Co-Seismic Ionospheric Disturbance around the time of 2015 Nepal earthquake

The Co-Seismic Ionospheric Disturbance of the 2015 Nepal earthquake is analyzed in this paper. GNSS data are used to obtain the Satellite-Station TEC sequences. After removing the de-trended TEC variation, a clear ionospheric disturbance was observed 10 min after the earthquake, while the geomagnetic conditions, solar activity, and weather condition remained calm according to the Kp, Dst, F10.7 indices and meteorological records during the period of interest. Computerized ionosphere tomography (CIT) is then used to present the tridimensional ionosphere variation with a 10-min time resolution. The CIT results indicate that (1) the disturbance of the ionospheric electron density above the epicenter during the 2015 Nepal earthquake is confined at a relatively low altitude (approximately 150–300 km); (2) the ionospheric disturbances on the west side and east sides of the epicenter are precisely opposite. A newly established electric field penetration model of the lithosphere–atmosphere–ionosphere coupling is used to investigate the potential physical mechanism.     

Performance of the L2C civil GPS signal under various ionospheric scintillation effects

As GPS is modernizing, there are currently fourteen satellites transmitting L2C civil code and seven satellites transmitting L5 signal. While the GPS observables are subject to several sources of errors, the ionosphere is one of the largest error sources affecting GPS signals. Small irregularities in the electrons density along the GPS radio signal propagation path cause ionospheric scintillation that is characterized by rapid fluctuations in the signal amplitude and phase. The ionospheric scintillation effects are stronger in equatorial and high-latitude geomagnetic latitude regions and occur mainly due to equatorial anomaly and solar storms. Several researchers have analyzed the L2C signal quality since becoming available in December, 2005. We analyze the performance of L2C using GPS data from stations in the equatorial region of Brazil, which is subject of weak, moderate and strong ionospheric scintillation conditions. The GPS data were collected by Septentrio PolaRxS–PRO receivers as part of the CIGALA/CALIBRA network. The analysis was performed as a function of scintillations indexes S4 and Phi60, lock time (time interval in seconds that the carrier phase is tracked continuously without cycle slips), multipath RMS and position variation of precise point positioning solutions. The analysis shows that L2C code solutions are less affected by multipath effects than that of P2 when data are collected under weak ionospheric scintillation effects. In terms of analysis of positions, the kinematic PPP results using L2C instead P2 codes show accuracy improvements up to 33 % in periods of weak or strong ionospheric scintillation. When combining phase and code collected under weak scintillation effects, the results by applying L2C against P2 provide improvement in accuracy up to 59 %. However, for data under strong scintillation effects, the use of L2C for PPP with code and phase does not provide improvements in the positioning accuracy.