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.     

Use of global ionospheric maps for HF Doppler measurements interpretation

The HF Doppler technique, a method of measurement of Doppler frequency shift of ionospheric signal, is one of the well-known and widely used techniques of ionosphere research. It allows investigation of various disturbances in the ionosphere. There are different sources of disturbances in the ionosphere such as geomagnetic storms, solar flashes, meteorological effects and atmospheric waves. The HF Doppler technique allows us to find out the influence of earthquakes, explosions and other processes on the ionosphere, which occurs near the Earth. HF Doppler technique has high sensitivity to small frequency variations and high time resolution but interpretation of results is difficult. In this paper, we attempt to use GPS data for Doppler measurements interpretation. Modeling of Doppler frequency shift variations with use of TEC allows separation of ionosphere disturbances of medium scale.     

Accuracy and consistency of different global ionospheric maps released by IGS ionosphere associate analysis centers

Due to the differences of ionospheric modeling methods and selected tracking stations, the accuracy and consistency of Global Ionospheric Maps (GIMs) released by Ionosphere Associate Analysis Centers (IAACs) are different. In this study, we evaluate and analyze in detail the accuracy and consistency of GIMs final products provided by six IAACs from three different aspects. Firstly, the comparison of these GIMs shows that the mean bias (MEAN) is related to the modeling methods of various IAACs. The variation trend of the standard deviation (STD) is consistent with the solar activities, and accompanied by certain seasonal and annual periodic variations. The MEAN between IGS and each center is about −1.3 to 1.0 TECU, and the STD is about 1.4–2.5 TECU. Secondly, the validation with GPS TEC shows that the STD of CODE is the smallest at various latitudes, and the STD is about 0.7–4.5 TECU. Thirdly, The validation with the Jason2 VTEC shows that the STD between Jason2 and IAACs is about 4.4–5.2 TECU. In addition, the STD between Jason2 and six GIMs in the areas with more tracking stations is better than that of the regions with fewer tracking stations in different latitude regions. Regardless of whether the tracking stations are more or less, the MEAN and STD in high solar activity are larger than in low solar activity.     

Regional 4-D modeling of the ionospheric electron density from satellite data and IRI

Accurate knowledge of the electron density is the key point in correcting ionospheric delays of electromagnetic measurements and in studying ionosphere physics. During the last decade Global Navigation Satellite Systems (GNSS) have become a promising tool for monitoring ionospheric parameters such as the total electron content (TEC). In this contribution we present a four-dimensional (4-D) model of the electron density consisting of a given reference part, i.e., the International Reference Ionosphere (IRI), and an unknown correction term expanded in terms of multi-dimensional base functions. The corresponding series coefficients are calculable from the satellite measurements by applying parameter estimation procedures. Since satellite data are usually sampled between GPS satellites and ground stations, finer structures of the electron density are modelable just in regions with a sufficient number of ground stations. The proposed method is applied to simulated geometry-free GPS phase measurements. The procedure can be used, for example, to study the equatorial anomaly.     

Characterisation of residual ionospheric errors in bending angles using GNSS RO end-to-end simulations

Global Navigation Satellite System (GNSS) radio occultation (RO) is an innovative meteorological remote sensing technique for measuring atmospheric parameters such as refractivity, temperature, water vapour and pressure for the improvement of numerical weather prediction (NWP) and global climate monitoring (GCM). GNSS RO has many unique characteristics including global coverage, long-term stability of observations, as well as high accuracy and high vertical resolution of the derived atmospheric profiles. One of the main error sources in GNSS RO observations that significantly affect the accuracy of the derived atmospheric parameters in the stratosphere is the ionospheric error. In order to mitigate the effect of this error, the linear ionospheric correction approach for dual-frequency GNSS RO observations is commonly used. However, the residual ionospheric errors (RIEs) can be still significant, especially when large ionospheric disturbances occur and prevail such as during the periods of active space weather. In this study, the RIEs were investigated under different local time, propagation direction and solar activity conditions and their effects on RO bending angles are characterised using end-to-end simulations. A three-step simulation study was designed to investigate the characteristics of the RIEs through comparing the bending angles with and without the effects of the RIEs. This research forms an important step forward in improving the accuracy of the atmospheric profiles derived from the GNSS RO technique.     

Higher order ionospheric propagation effects on GPS radio occultation signals

With the increasing number of remote sensing satellites using the GPS radio occultation technique for atmospheric sounding, the estimation of higher order ionospheric effects and their mitigation have become relevant and important. Due to long ionospheric limb paths, GPS signals are strongly affected by ionospheric refraction during radio occultation. Standard dual-frequency GPS measurements may be used to estimate the first order term of the refractive index. However, non-linear terms such as the second and third order ionospheric terms and ray path bending effects are not considered in occultation measurements so far. Analysing selected CHAMP–GPS occultation events different higher order ionospheric terms are estimated and their effects on dual-frequency range estimation and total electron content (TEC) estimation are discussed. We have found that the separation between the GPS L1 and L2 ray paths exceeds the kilometer level during occultation for a vertical TEC level of more than 160 TEC units. Corresponding errors in the GPS dual-frequency range estimation and TEC estimation are found to exceed the meter and 10 TEC units level, respectively.     

The northern and southern mid-latitude ionospheric trough using global IGS vTEC maps

The aim of this work is the analysis of the mid-latitude ionospheric trough (MIT) using the Global Ionospheric Maps from IGS (GIMs) during the solar minimum, year 2008. This study was performed for different local times, 22, 00, 02 and 04 LT on the Northern and Southern hemisphere simultaneously. In the two hemispheres the MIT show asymmetric pattern. The high-latitude troughs are clearly distinguished in autumn and winter. Another feature is the longitudinal development towards the west of the geomagnetic pole covering a wider area in the Northern Hemisphere. Five empirical reference models were tested and compared with the MIT minimum position obtained from GIMs at different local times for both hemisphere. The results show a better agreement with the observations for the Northern Hemisphere specially with the Köehnlein & Raitt model. Fluctuations of 9 days and 27 days of the MIT minimum position are found, which could be related with the solar wind oscillations, especially for 00 and 02 LT in both hemisphere, suggesting a link between them.     

Study of GNSS-measured ionospheric total electron content variations generated by powerful HF-heating

The purpose of this work is to report the experimental evidences for the influence of perturbations in the electron density in the dayside mid-latitude ionosphere, that are caused by high-frequency heating of the F2 layer, on the GNSS signals. The experiments were carried out at the Sura heater (Radio Physical Research Institute, N. Novgorod). During the sessions of ionospheric heating with different time modulations of the radiated power the rays linking the navigational satellites with the ground receiver intersected the heated region. Variations in the total electron content (TEC) were studied; these variations are proportional to the reduced phases of navigational signals. It is shown that with the square-wave modulation of the radiated power (with periods of 1, 6, 10 and 15 min), perturbations with periods of the main modulation of heating and its harmonics appear in the spectrum of TEC variations. Examples are presented of identification of the heating-induced variations in TEC, including determination of the amplitudes and time characteristics of these variations.     

Dynamics of global electron content in 1998–2005 derived from global GPS data and IRI modeling

We analyzed the dynamics of global electron content (GEC) for the period 1998–2005 and compared the estimated GEC with variations of the 10.7-cm solar radio emission and with and with GEC values obtained with IRI-2001. We found a strong resemblance between the curves’ shapes for the experimental and modeled GEC: strong semiannual variations are discernible in these series and both curves tend to increase the absolute GEC value during the period of maximum of solar activity. However, there are some significant distinctions, such as absence of 27-day fluctuations in the series of GEC computed by the IRI-2001. On the contrary, observational GEC reflects well dynamics of solar activity: 27-day variations of GEC are very similar to the ones of the index F10.7, but GEC undergoes a lagging of about of 30–60 h as compared to value of the F10.7 index. The relative amplitude of 27-day variations decreases from 8% at the rising and falling solar activity to 2% at the period of its maximum.     

Evaluation of COMPASS ionospheric model in GNSS positioning

As important products of GNSS navigation message, ionospheric delay model parameters are broadcasted for single-frequency users to improve their positioning accuracy. GPS provides daily Klobuchar ionospheric model parameters based on geomagnetic reference frame, while the regional satellite navigation system of China’s COMPASS broadcasts an eight-parameter ionospheric model, COMPASS Ionospheric Model(CIM), which was generated by processing data from continuous monitoring stations, with updating the parameters every 2 h. To evaluate its performance, CIM predictions are compared to ionospheric delay measurements, along with GPS positioning accuracy comparisons. Real observed data analysis indicates that CIM provides higher correction precision in middle-latitude regions, but relatively lower correction precision for low-latitude regions where the ionosphere has much higher variability. CIM errors for some users show a common bias for in-coming COMPASS signals from different satellites, and hence ionospheric model errors are somehow translated into the receivers’ clock error estimation. In addition, the CIM from the China regional monitoring network are further evaluated for global ionospheric corrections. Results show that in the Northern Hemisphere areas including Asia, Europe and North America, the three-dimensional positioning accuracy using the CIM for ionospheric delay corrections is improved by 7.8%–35.3% when compared to GPS single-frequency positioning ionospheric delay corrections using the Klobuchar model. However, the positioning accuracy in the Southern Hemisphere is degraded due apparently to the lack of monitoring stations there.     

Large-scale traveling ionospheric disturbances of auroral origin according to the data of the GPS network and ionosondes

The intensity of large-scale traveling ionospheric disturbances (LS TIDs), registered using measurements of total electron content (TEC) during the magnetic storms on October 29–31, 2003, and on November 7–11, 2004, had been compared with that of local electron density disturbances. The data of TEC measurements at ground-based GPS receivers located near the ionospheric stations and the corresponding values of the critical frequency of the ionospheric F region f_oF2 were used for this purpose. The variations of TEC and f_oF2 were similar for all events mentioned above. The previous assumption that the ionospheric region with vertical extension from 150 to 200 km located near the F-layer maximum mainly contributes to the TEC variations was confirmed for the cases when the electron density disturbance at the F region maximum was not more than 50%. However, this region probably becomes vertically more extended when the electron density disturbance in the ionospheric F region is about 85%.     

星载单频GPS数据的电离层延迟改正方法分析

利用C/A码单点定位对LEO(Low Earth Orbit)卫星上的电离层延迟改正方法——"电离层比例因子法"进行了分析研究.计算的CHAMP卫星的轨道结果表明:采用电子密度峰值高度(hmF2,F2 region maximum electron density height)平均值和瞬时值计算的电离层比例因子α变化范围分别为0.3~0.4和0.2~0.65之间,两者最大差异可达0.3,相比较而言,hmF2瞬时值的结果更加合理,并且相应的大地高H方向的系统偏差要降低0.05~0.3m左右;与双频无电离层组合的普通单点定位结果相比表明该方法能较好地消除电离层一阶项所引入的H方向上的系统偏差;该方法适用的LEO卫星轨道高度范围大致在200~ 600km之间,当轨道高度超过700km时,该方法并不适用.     

The impact on tsunami detection from space using GNSS-reflectometry when combining GPS with GLONASS and Galileo

The devastating Sumatra tsunami in 2004 demonstrated the need for a tsunami early warning system in the Indian Ocean. Such a system has been installed within the German-Indonesian Tsunami Early Warning System (GITEWS) project. Tsunamis are a global phenomenon and for global observations satellites are predestined. Within the GITEWS project a feasibility study on a future tsunami detection system from space has therefore been carried out. The Global Navigation Satellite System Reflectometry (GNSS-R) is an innovative way of using GNSS signals for remote sensing. It uses ocean reflected GNSS signals for sea surface altimetry. With a dedicated Low Earth Orbit (LEO) constellation of satellites equipped with GNSS-R receivers, densely spaced sea surface height measurements could be established to detect tsunamis. Some general considerations on the geometry between LEO and GNSS are made in this simulation study. It exemplary analyzes the detection performance of a GNSS-R constellation at 900 km altitude and 60° inclination angle when applied to the Sumatra tsunami as it occurred in 2004. GPS is assumed as signal source and the combination with GLONASS and Galileo signals is investigated. It can be demonstrated, that the combination of GPS and Galileo is advantageous for constellations with few satellites while the combination with GLONASS is preferable for constellations with many satellites. If all three GNSS are combined, the best detection performance can be expected for all scenarios considered. In this case an 18 satellite constellation will detect the Sumatra tsunami within 17 min with certainty, while it takes 53 min if only GPS is considered.     

Comparison of electron density measurements from CSES and Swarm satellites with GNSS ionospheric tomography data

Electron density measurements obtained from China Seismo‐Electromagnetic Satellite (CSES) and Swarm-B can play an increasingly important role in the study of ionosphere above F2 peak height. This study presented a comprehensive comparison of electron density products obtained from Langmuir probe mounted on CSES and Swarm-B with ionospheric tomography for a whole year period of 2019. CSES was fully compared with Swarm-B on a global scale, including both absolute and relative differences, and a new index called NFI was developed to better quantify the similarity between two latitudinal profiles of electron density. CSES and Swarm-B were then compared with tomography respectively in four regions, roughly located in America, Europe, Australia and China. Results indicated that CSES data are consistent with Swarm-B, as NFI values exceed 0.6 for most of the analyses. Tomography and Swarm-B were found to have a good agreement as their biases are less than 0.2 × 10⁵ el/cm³ in general. For the comparison between CSES and tomography, the bias increased to around 0.6 × 10⁵ el/cm³ but the standard deviation changed slightly, validating the underestimation of electron density by CSES. The spatiotemporal comparisons of CSES and Swarm-B with tomography showed that: 1) the differences in electron density were relatively low in middle latitudes and increased rapidly in the regions of equatorial ionization anomaly; 2) Swarm-B has a better consistent with tomography than CSES, but both are capable of detecting ionosphere anomalies such as midlatitude arcs; and 3) CSES and Swarm-B both can capture the seasonal changes of electron density, while their values are basically smaller than those from tomography in Spring and Summer months.     

Use of principal component analysis in the identification of the spatial pattern of an ionospheric total electron content anomalies after China’s May 12, 2008, M = 7.9 Wenchuan earthquake

This paper uses principal component analysis (PCA) to determine the spatial pattern of total electron content (TEC) anomalies in the ionosphere post the China’s Wenchuan Earthquake of 12 May, 2008 (UTC) (M_w = 7.9). PCA and image processing are applied to the global ionospheric map (GIM) with transforms conducted for the time period 08:00–10:00 UT on 12 May 2008. Results show that at a height of approximately 200 km the anomaly is widespread and less intense; however, it becomes more localized with height reaching maximum intensity and localization at an altitude of 300 km. The spatial distribution is remarkably similar to that reported for a TEC anomaly previously identified as a precursor anomaly on May 9, 2008 for the same time period. Potential causes of the results are discussed with emphasis given to vertical acoustic gravity waves based on the spatial pattern identified.     

Stochastic modelling considering ionospheric scintillation effects on GNSS relative and point positioning

Global Navigation Satellite Systems (GNSS), in particular the Global Positioning System (GPS), have been widely used for high accuracy geodetic positioning. The Least Squares functional models related to the GNSS observables have been more extensively studied than the corresponding stochastic models, given that the development of the latter is significantly more complex. As a result, a simplified stochastic model is often used in GNSS positioning, which assumes that all the GNSS observables are statistically independent and of the same quality, i.e. a similar variance is assigned indiscriminately to all of the measurements. However, the definition of the stochastic model may be approached from a more detailed perspective, considering specific effects affecting each observable individually, as for example the effects of ionospheric scintillation. These effects relate to phase and amplitude fluctuations in the satellites signals that occur due to diffraction on electron density irregularities in the ionosphere and are particularly relevant at equatorial and high latitude regions, especially during periods of high solar activity. As a consequence, degraded measurement quality and poorer positioning accuracy may result.     

A real-time ionospheric model based on GNSS Precise Point Positioning

This paper proposes a method of real-time monitoring and modeling the ionospheric Total Electron Content (TEC) by Precise Point Positioning (PPP). Firstly, the ionospheric TEC and receiver’s Differential Code Biases (DCB) are estimated with the undifferenced raw observation in real-time, then the ionospheric TEC model is established based on the Single Layer Model (SLM) assumption and the recovered ionospheric TEC. In this study, phase observations with high precision are directly used instead of phase smoothed code observations. In addition, the DCB estimation is separated from the establishment of the ionospheric model which will limit the impacts of the SLM assumption impacts. The ionospheric model is established at every epoch for real time application. The method is validated with three different GNSS networks on a local, regional, and global basis. The results show that the method is feasible and effective, the real-time ionosphere and DCB results are very consistent with the IGS final products, with a bias of 1–2 TECU and 0.4 ns respectively.     

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

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

Multi-GNSS contributions to differential code biases determination and regional ionospheric modeling in China

Due to the limited number and uneven distribution globally of Beidou Satellite System (BDS) stations, the contributions of BDS to global ionosphere modeling is still not significant. In order to give a more realistic evaluation of the ability for BDS in ionosphere monitoring and multi-GNSS contributions to the performance of Differential Code Biases (DCBs) determination and ionosphere modeling, we select 22 stations from Crustal Movement Observation Network of China (CMONOC) to assess the result of regional ionospheric model and DCBs estimates over China where the visible satellites and monitoring stations for BDS are comparable to those of GPS/GLONASS. Note that all the 22 stations can track the dual- and triple-frequency GPS, GLONASS, and BDS observations. In this study, seven solutions, i.e., GPS-only (G), GLONASS-only (R), BDS-only (C), GPS + BDS (GC), GPS + GLONASS (GR), GLONASS + BDS (RC), GPS + GLONASS + BDS (GRC), are used to test the regional ionosphere modeling over the experimental area. Moreover, the performances of them using single-frequency precise point positioning (SF-PPP) method are presented. The experimental results indicate that BDS has the same ionospheric monitoring capability as GPS and GLONASS. Meanwhile, multi-GNSS observations can significantly improve the accuracy of the regional ionospheric models compared with that of GPS-only or GLONASS-only or BDS-only, especially over the edge of the tested region which the accuracy of the model is improved by reducing the RMS of the maximum differences from 5–15 to 2–3 TECu. For satellite DCBs estimates of different systems, the accuracy of them can be improved significantly after combining different system observations, which is improved by reducing the STD of GPS satellite DCB from 0.243 to 0.213, 0.172, and 0.165 ns after adding R, C, and RC observations respectively, with an increment of about 12.3%, 29.4%, and 32.2%. The STD of GLONASS satellite DCB improved from 0.353 to 0.304, 0.271, and 0.243 ns after adding G, C, and GC observations, respectively. The STD of BDS satellite DCB reduced from 0.265 to 0.237, 0.237 and 0.229 ns with the addition of G, R and GR systems respectively, and increased by 10.6%, 10.4%, and 13.6%. From the experimental positioning result, it can be seen that the regional ionospheric models with multi-GNSS observations are better than that with a single satellite system model.     

Assessment of the NeQuick model at mid-latitudes using GNSS TEC and ionosonde data

The modelling of the total electron content (TEC) plays an important role in global navigation satellite systems (GNSS) accuracy, especially for single-frequency receivers, the most common ones constituting the mass market. For the latter and in the framework of Galileo, the NeQuick model has been chosen for correcting the ionospheric error contribution and will be integrated into a global algorithm providing the users with daily updated information.