SHAKING: Adjusted spherical harmonics adding KrigING method for near real-time ionospheric modeling with multi-GNSS observations

Precise positioning based on Global Navigation Satellite System (GNSS) technique requires high accuracy ionospheric total electron content (TEC) correction models to account for the ionospheric path delay errors. We present an adjusted Spherical Harmonics Adding KrigING method (SHAKING) approach for regional ionospheric vertical TEC (VTEC) modeling in real time. In the proposed SHAKING method, the VTEC information over the sparse observation data area is extrapolated by the Adjusted Spherical Harmonic (ASH) function, and the boundary distortion in regional VTEC modeling is corrected by the stochastic VTEC estimated using Kriging interpolation. Using real-time GPS, GLONASS and BDS-2/3 data streams of the Crust Movement Observation Network of China (CMONOC), the SHAKING-based regional ionospheric VTEC maps are re-constructed over China and its boundary regions. Compared to GNSS VTECs derived from the independent stations, the quality of SHAKING solution improves by 13–31% and 6–33% with respect to the ASH-only solution during high and low geomagnetic periods, respectively. Compared to the inverse distance weighting (IDW) generated result, significant quality improved of SHAKING-based VTEC maps is also observed, especially over the edge areas with an improvement of 60–80%. Overall, the proposed SHAKING method exhibits notable advantage over the existing regional VTEC modeling techniques, which can be used for regional TEC modeling and associated high-precision positioning applications.     

基于NTCM-BC模型的全球卫星导航系统单频电离层延迟修正

选择NTCM-BC模型作为单频电离层延迟修正模型,通过非线性最小二乘拟合的方法,利用提前一天预测的电离层图（COPG文件）,计算得到NTCM-BC模型修正系数;利用Klobuchar模型和IGS发布的GIM数据对NTCM-BC模型进行比较和分析.对太阳活动高、中、低年实测数据的分析结果表明:全球平均水平上,NTCM-BC模型的电离层延迟修正性能明显优于Klobuchar模型,NTCM-BC模型的TEC平均误差和均方根误差比Klobuchar模型分别下降了41%和30%;模型的TEC计算误差与太阳活动剧烈程度成正相关,即太阳活动高年模型误差较大,太阳活动低年误差相对较低.相较于磁静日,磁扰日期间Klobuchar模型和NCTM模型的误差均有一定程度的增加.此外,模型的电离层修正误差同时存在明显的纬度、季节和地方时差异.      

一种基于GNSS数据的全球电离层不规则结构活动表征方法

针对如何利用GNSS(Global Navigation Satellite System)数据进行电离层扰动监测的问题,提出了一种基于GNSS数据表征全球电离层扰动的方法.利用大约400个GNSS地面站点的观测数据,计算总电子含量(Total Electron Content,TEC)变化率的标准差——ROTI(Ra...     

地磁暴期间北半球高纬度地区电离层变化特征及对精密定位的影响

基于加拿大地区高纬度电离层观测网的电离层闪烁观测数据,分析了2018年8月26日地磁暴事件引发的北半球高纬度地区电离层总电子含量（TEC）异常变化、TEC变化率指数（ROTI）及电离层相位闪烁的变化特征.结果表明:加拿大地区最大异常值约6 TECU,磁暴引发全球电离层TEC异常峰值高达20 TECU;加拿大地区电离层相位闪烁发生率最大增至12.6%,而磁静日期间约为1%;强电离层闪烁期间,电离层相位闪烁指数与ROTI之间具有较强的一致性.对GPS双频精密单点定位（Precise Point Positioning,PPP）结果进行分析发现:无闪烁期间定位误差随测站纬度的增高呈现出增大趋势,但均方根误差小于0.4m;闪烁发生期间各测站的定位误差均显著增大,水平和垂直方向定位误差分别增至约0.9m及1.7m.      

南极GPS定位中的电离层延迟分析

基于南极地区国际GNSS服务组织(IGS,International GNSS Service)跟踪站的全球定位系统(GPS,Global Position System)双频实测数据,分析了南极地区电离层延迟的变化情况及其二阶项延迟对南极GPS定位结果的影响.结果表明:南极地区的总电子含量(TEC,Total Electron Content)日间波动频繁,其日间TEC最大值变化较中纬度地区剧烈;在南极地区夏季,电离层二阶项延迟对GPS定位结果的影响可达cm级.同时,由于欧洲定轨中心(CODE,The European Center for Orbit Determination)提供的全球电离层模型(GIM,Global Ionosphere Maps)在南极区域应用的局限性,通过选取南极地区6个IGS跟踪站作为基准站建立了区域电离层TEC模型(RIM,Regional Ionosphere Model).经实测数据计算证明,对于南极地区,RIM的定位精度在一定程度上优于全球电离层模型GIM.     

An improved approach to model regional ionosphere and accelerate convergence for precise point positioning

Given the severe effects of the ionosphere on global navigation satellite system (GNSS) signals, single-frequency (SF) precise point positioning (PPP) users can only achieve decimeter-level positioning results. Ionosphere-free combinations can eliminate the majority of ionospheric delay, but increase observation noise and slow down dual-frequency (DF) PPP convergence. In this paper, we develop a regional ionosphere modeling and rapid convergence approach to improve SF PPP (SFPPP) accuracy and accelerate DF PPP (DFPPP) convergence speed. Instead of area model, ionospheric delay is modeled for each satellite to be used as a priori correction. With the ionospheric, wide-lane uncalibrated phase delay (UPD) and residuals satellite DCBs product, the wide-lane observations for DF users change to be high-precision pseudorange observations. The validation of a continuously operating reference station (CORS) network was analyzed. The experimental results confirm that the approach considerably improves the accuracy of SFPPP. For DF users, convergence time is substantially reduced.     

45°(N)纬度带的Klobuchar like电离层延迟季节修正模型与评估

电离层延迟是全球卫星导航系统（GNSS）的主要误差源之一。对于装配GNSS单频接收机的航空器,选择简单有效的Klobuchar广播电离层模型来改正电离层延迟误差,其修正率为50%～60%。针对45°（N）纬度带,提出了更高电离层修正需求。考虑到季节因素对中高纬度地区电离层的显著影响,利用GIMs（Global Ionospheric Maps）分析了昼夜中TEC（Total Electron Content）的峰值和谷值随季节（年积日）的变化,建立了一种适用于45°（N）纬度带的Klobuchar like电离层模型。该模型不增加广播模型系数,新模型的夜间和VTEC高峰时电离层修正率分别达到了82%和80%,表明在穿刺点集中的45°(N)纬度地区使用该模型可以更精确地描述该地区的电离层,帮助航空器实现更高精度的定位。     

Estimation of QZSS differential code biases using QZSS/GPS combined observations from MGEX

As an important error source in Global Navigation Satellite System (GNSS) positioning and ionospheric modeling, the differential code biases (DCB) need to be estimated accurately, e.g., the regional Quasi-Zenith satellite system (QZSS). In this paper, the DCB of QZSS is estimated by adopting the global ionospheric modeling method based on QZSS/GPS combined observations from Multi-GNSS experiment (MGEX). The performance of QZSS satellite and receiver DCB is analyzed with observations from day of year (DOY) 275–364, 2018. Good agreement between our estimated QZSS satellite DCB and the products from DLR and CAS is obtained. The bias and root mean square (RMS) of DCB are mostly within ±0.3 ns. The day-to-day fluctuation of the DCB time series is less than 0.5 ns with about 96% of the cases for all satellites. However, the receiver DCB is a little less stable than satellite DCB, and their standard deviations (STDs) are within 1.9 ns. The result shows that the stability of the receiver DCBs is not significantly related to the types of receiver or antenna.     

Ionospheric anomalies observed over South Korea preceding the Great Tohoku earthquake of 2011

We investigated the ionospheric anomalies observed before the Tohoku earthquake, which occurred near the northeast coast of Honshu, Japan on 11 March, 2011. Based on data from a ground-based Global Positioning System (GPS) network on the Korean Peninsula, ionospheric anomalies were detected in the total electron content (TEC) during the daytime a few days before earthquake. Ionospheric TEC anomalies appeared on 5, 8 and 11 March. In particular, the ionospheric disturbances on 8 March evidenced a remarkable increase in TEC. The GPS TEC variation associated with the Tohoku earthquake was an increase of approximately 20 total electron content units (TECU), observed simultaneously in local and global TEC measurements. To investigate these pre-earthquake ionospheric anomalies, space weather conditions such as the solar activity index (F10.7) and geomagnetic activity indices (the Kp and Dst indices) were examined. We also created two-dimensional TEC maps to visual the spatial variations in the ionospheric anomalies preceding the earthquake.     

On the importance of total electron content enhancements during the extreme solar minimum

Observations of ionospheric vertical total electron content (vTEC) from European ground-based Global Navigation Satellite Systems (GNSS) receivers during the period January 2008–January 2010 are used to investigate, for the first time, vTEC sensitivity to weak geomagnetic disturbances under extreme solar minimum conditions. This study shows a significant number of events for the period in question, all of which exhibited some form of exceptionally large values of vTEC during small-magnitude geomagnetic disturbances. To illustrate our point on the importance of vTEC enhancements during the extreme solar minimum and its relevance for the current GNSS and future Galileo applications, we present in this paper the results associated with two significant events that both occurred in equinoctial months. The 10–12 October 2009 event of anomalous TEC enhancement at two distant mid-latitude locations HERS (0.3 E; 50.9 N) and NICO (33.4 E, 35.1 N) is discussed in the context of strong vTEC variations during the well established ionospheric storm on 11 October 2008. We conclude with a short summary of the new findings and their consequences on ionospheric monitoring and modelling for operational communication and navigation systems.     

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.     

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.     

基于空间统计方法的电离层折射修正技术

针对中国上空电离层所具有的特殊性和GPS观测站在中国西部分布相对稀疏的特点, 尝试探索中国卫星增强系统电离层时延信息修正技术, 为卫星导航定位以及遥感、遥测等空间应用工程的电波修正提供数据. 利用中国地壳形变监测网提供的双频GPS数据, 以空间统计方法为主要工具, 给出了普通Kriging电离层估计算法, 构建了平静期和磁暴期电离层理论变异模型, 详细分析了电离层折射修正的精度. 结果表明, 将空间统计方法应用于卫星增强系统中的电离层时延改正问题, 有利于提高增强系统的电离层折射修正精度, 特别是在观测样点相对较少的情况下, 有利于系统完整性的实现.      

Multi-station observation of ionospheric irregularities over South Africa during strong geomagnetic storms

This paper presents results pertaining to the response of the mid-latitude ionosphere to strong geomagnetic storms that occurred from 31 March to 02 April 2001 and 07–09 September 2002. The results are based on (i) Global Positioning Systems (GPSs) derived total electron content (TEC) variations accompanying the storm, (ii) ionosonde measurements of the ionospheric electrodynamic response towards the storms and (iii) effect of storm induced travelling ionospheric disturbances (TIDs) on GPS derived TEC. Ionospheric data comprising of ionospheric TEC obtained from GPS measurements, ionograms, solar wind data obtained from Advanced Composition Explorer (ACE) and magnetic data from ground based magnetometers were used in this study. Storm induced features in vertical TEC (VTEC) have been obtained and compared with the mean VTEC of quiet days. The response of the mid-latitude ionosphere during the two storm periods examined may be characterised in terms of increased or decreased level of VTEC, wave-like structures in VTEC perturbation and sudden enhancement in hmF2 and h′F. The study reveals both positive and negative ionospheric storm effects on the ionosphere over South Africa during the two strong storm conditions. These ionospheric features have been mainly attributed to the travelling ionospheric disturbances (TIDs) as the driving mechanism for the irregularities causing the perturbations observed. TEC perturbations due to the irregularities encountered by the satellites were observed on satellites with pseudo random numbers (PRNs) 15, 17, 18 and 23 between 17:00 and 23:00 UT on 07 September 2002.     

Mitigation of receiver biases in ionospheric observables from PPP with ambiguity resolution

PPP (Precise Point Positioning) is a GNSS (Global Navigation Satellite Systems) positioning method that requires SSR (State Space Representation) corrections in order to provide solutions with an accuracy of centimetric level. The so-called RT-PPP (Real-time PPP) is possible thanks to real-time precise SSR products, for orbits and clocks, provided by IGS (International GNSS Service) and its associate analysis centers such as CNES (Centre National d'Etudes Spatiales). CNES SSR products also enable RT-PPP with integer ambiguity resolution. In GNSS related literature, PPP with ambiguity resolution (PPP-AR) in real-time is often referred as PPP-RTK (PPP – Real Time Kinematic). PPP-WIZARD (PPP - With Integer and Zero-difference Ambiguity Resolution Demonstrator) is a software that is made available by CNES. This software is capable of performing PPP-RTK. It estimates slant ionospheric delays and other GNSS positioning parameters. Since ionospheric effects are spatially correlated by GNSS data from active networks, it is possible to model and provide ionospheric delays for any position in the network coverage area. The prior knowledge ionospheric delays can reduce positioning convergence for PPP-RTK users. Real-time ionospheric models could benefit from highly precise ionospheric delays estimated in PPP-AR. In this study, we demonstrate that ionospheric delays obtained throughout PPP-AR estimation are actu ally ionospheric observables. Ionospheric observables are biased by an order of few meters caused by the receiver hardware biases. These biases prohibit the use of PPP-WIZARD ionospheric delays to produce ionospheric models. Receiver biases correction is essential to provide ionospheric delays while using PPP-AR based ionospheric observables. In this contribution, a method was implemented to estimate and mitigate receiver hardware biases influence on slant ionospheric observables from PPP-AR. In order to assess the proposed approach, PPP-AR data from 12 GNSS stations were processed over a two-month period (March and April 2018). A comparison between IGS ionospheric products and PPP-AR based ionospheric observables corrected for receiver biases, resulted in a mean of differences of −39 cm and 51 cm standard deviation. The results are consistent with the accuracy of the IGS ionospheric products, 2–8 TECU, considering that 1 TECU is ~16 cm in L1. In another analysis, a comparison of ionospheric delays from 5 pairs of short baselines GNSS stations found an agreement of 0.001 m in mean differences with 22 cm standard deviation after receiver biases were corrected. Therefore, the proposed solution is promising and could produce high quality (1–2 TECU) slant ionospheric delays. This product can be used in a large variety of modeling approaches, since ionospheric delays after correction are unbiased. These results indicate that the proposed strategy is promising, and could benefit applications that require accuracy of 1–2 TECU (~16–32 cm in L1).     

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.     

Ionospheric TEC disturbances over China during the strong geomagnetic storm in September 2017

From September 7 to 8, 2017, a G4-level strong geomagnetic storm occurred, which seriously impacted on the Earth’s ionosphere. In this work, the global ionospheric maps released by Chinese Academy of Sciences are used to investigate the ionospheric responses over China and its adjacent regions during the strong storm. The prominent TEC enhancements, which mainly associated with the neutral wind and eastward prompt penetration electric field, are observed at equatorial ionization anomaly crests during the main phase of the storm on 8 September 2017. Compared with those on 8 September, the TEC enhancements move to lower-latitude regions during the recovery phase on 9 September. A moderate storm occurred well before the start of the strong storm causes similar middle-latitude TEC enhancements on 7 September. However, the weak TEC depletion is observed at middle and low latitude on 9–10 September, which could be associated with the prevailing westward disturbance electric field or storm-time neural composition changes. In addition, the storm-time RMS and STD values of the ionospheric TEC grids over China increase significantly due to the major geomagnetic storm. The maximum of the RMS reaches 12.0 TECU, while the maximum of the STD reaches 8.3 TECU at ~04UT on 8 September.     

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.     

The ionosphere,radio navigation,and global navigation satellite systems

This article is a review of Global Navigation Satellite Systems (GNSS) for space scientists who are interested in how GNSS signals and observables can be used to understand ionospheric dynamics and, conversely, how ionospheric dynamics affect the operational capabilities of GNSS receivers. The most common form of GNSS is the Global Positioning System (GPS); we will first review its operating principles and then present a discussion of errors, of which ionospheric propagation is the most significant. Methods and systems for mitigating errors will be introduced, along with a discussion of modernization plans for GPS and for entirely new systems such as Galileo. In the second half of this article the effects of the ionosphere on GPS signals will be examined in more detail, particularly ionospheric propagation, leading to a discussion of the relation of TEC to ranging errors. Next, the subject of scintillations will also be introduced and connected to the presence and scale sizes of irregularities. Scintillations will be examined as spatial and temporal structures. The method of measuring scintillation pattern drift and ionospheric velocity will be discussed. We conclude by examining ionospheric effects on GPS at midlatitudes.