On the estimation of regional covariance functions of TEC variations over Canada

Spatial and temporal variations of Total Electron Content (TEC) can affect GNSS high accuracy positioning. Enhanced estimation of ionospheric variations and their de-correlation can benefit differential and point positioning rapid solutions. Global and regional TEC maps can provide the overall state of ionopsheric variations in space and time domains within their accuracy limits. In this paper, these maps are exploited to retrieve ionospheric variations by means of variograms and their associated covariance functions of TEC residuals over Canadian region during geomagnetically quiet and disturbed conditions. A number of theoretical variogram functions are reviewed for modeling covariance of TEC residuals. The variogram modeling of residuals during a strong geomagnetic storm revealed variances of one order of magnitude larger compared to a rather quiet condition. Variogram models are also used in regional and local kriging interpolation experiments and their performances are evaluated. Global maps of TEC RMS by International GNSS Service and two of its analysis centres are also compared over the Canadian region during a two-year period. Realistic representation of regional variances using estimated variograms when compared to global ionospheric RMS maps are also presented.     

A MARS-based method for estimating regional 2-D ionospheric VTEC and receiver differential code bias

The geometry-free linear combination of dual-frequency GNSS reference station ground observations are currently used to build the Vertical Total Electron Content (VTEC) model of the ionosphere. As it is known, besides ionospheric delays, there are differential code bias (DCB) of satellite (SDCB) and receiver (RDCB) in the geometry-free observation equation. The SDCB can be obtained using the International GNSS Service (IGS) analysis centers, but the RDCB for regional and local network receivers are not provided. Therefore, estimating the RDCB and VTEC model accurately and simultaneously is a critical factor investigated by researchers. This study uses Multivariate Adaptive Regression Splines (MARS) to estimate the VTEC approximate model and then substitutes this model in the observation equation to form the normal equation. The least squares method is used to solve the RDCB and VTEC model together. The research findings show that this method has good modeling effectiveness and the estimated RDCB has good reliability. The estimated VTEC model applied to GPS single-frequency precise point positioning has better positioning accuracy in comparison to the IGS global ionosphere map (GIM).     

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

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

南极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.     

一种用于电离层TEC监测的GNSS信号载波跟踪算法

全球卫星导航系统（GNSS）是电离层TEC监测中应用最普遍的手段. 目前方法通常是在传统导航用途的GNSS接收机输出的原始观测量基础上,经过数据后处理得到电离层TEC信息,其GNSS信号的跟踪处理算法依然采用GNSS导航接收机的算法. 针对GNSS系统用于电离层TEC监测的特殊性,提出一种称为GNSS双频信号和差联合跟踪的新算法,与传统方法相比,该算法直接跟踪电离层TEC的变化,可以提高电离层TEC跟踪的灵敏度和TEC的观测精度,改善电离层TEC监测性能.      

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

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

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.     

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).     

基于极大验后推估理论的全球电离层地图预报

提出了一种基于极大验后估计理论的全球电离层预报方法，基于中国科学院电离层分析中心（CAS）提供的快速全球电离层地图（GIM），实现了1天、2天和5天GIM的预报。以国际GNSS服务组织（IGS）最终GIM、Jason测高卫星提供的电离层观测信息及全球GNSS基准站实测电离层总电子含量（TEC）为基准，评估了2008－2017年CAS电离层预报GIM在全球大陆及海洋区域的精度，并与欧洲定轨中心（CODE）、欧洲空间局（ESA）和西班牙加泰罗尼亚理工大学（UPC）的预报GIM进行对比。在评估时段内，与IGS-GIM相比，CAS预报GIM精度为2.4~3.1 TECU；与测高卫星TEC相比，CAS预报GIM的精度为5.1~6.6 TECU；与全球基准站实测TEC相比，CAS预报GIM的电离层延迟修正精度优于80%。总体来看，CAS预报GIM与CODE预报GIM精度相当，显著优于ESA和UPC预报GIM。      

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.     

2009年1月平流层爆发性增温期间全球电离层响应的研究

2009年1月平流层爆发性增温(Stratospheric Sudden Warming, SSW)事件是有记录以来最强、持续时间最长的一次主增温事件(Major Warming Event, MWE), 期间太阳活动和地磁活动均处于较低的水平, 因此非常有利于研究电离层对平流层增温事件的响应情况. 本文利用COSMIC (Constellation Observing System for Meteorology, Ionosphere, and Climate)系统提供的掩星数据, 使用Kriging方法分别构建了此次SSW期间及平静期的电离层N_mF₂, h_mF₂和110～750km高度范围的垂直积分TEC (简称VTEC)地图, 同时从全球定位导航卫星系统服务组织(International GNSS Service, IGS)发布的全球电离层TEC地图(Global Ionospheric Maps, GIMs)中提取了日固坐标系(Sun-fixed, 地磁纬度和地方时)下全球TEC地图. 通过对比发现, SSW期间与平静期相比, 地磁纬度中低纬电离层参数存在早晨上升, 下午和夜间下降的现象. 利用OSTM/JASON-2卫星高度计观测值进行验证后的结果显示, 此前研究均未有提及的夜间时段电离层参数N_mF₂, h_mF₂和TEC (VTEC和IGS TEC)的下降现象的确存在.      

Analysis of the data processing strategies of spherical harmonic expansion model on global ionosphere mapping for moderate solar activity

Spherical harmonic (SH) expansion is widely used to model the global ionosphere map (GIM) of vertical total electron content (VTEC). According to the impact of different data processing methods of the SH expansion model on the VTEC maps, we specifically performed comprehensive analysis in terms of the data sampling rate, the time resolution, the spherical harmonic degree, and the relative constraint. One month of GPS data (January in 2016) from the International GNSS (Global Navigation Satellite System) Service (IGS) network in a moderate ionospheric activity period at the descending phase of Solar Cycle 24 was processed. To improve the computational efficiency of the daily GIM generation, the data sampling rate of 5?min was recommended allowing the GIM precision loss within 0.10 TECU (total electron content unit). The global VTEC map could be better represented in temporal and spatial domains with higher time resolution and higher spherical harmonic degree, especially at low latitude bands and in the southern hemisphere. The GIM precision improvement was about 10.91% for 1-h and about 15.15% for 0.5-h compared with the commonly used 2-h time resolution. The use of spherical harmonic degree 17 or 20 instead of 15 could improve the precision by 3.19% or 6.06%. We also found that an optimal relative constraint had to be found experimentally considering both the GIM precision and the GIM root mean square (RMS) map.     

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.     

Pre-earthquake ionospheric anomalies before three major earthquakes by GPS-TEC and GIM-TEC data during 2015–2017

The dual-frequency satellite-based measurements from Global Positioning System (GPS) may provide feasible ways of studying and potentially detecting of earthquake (EQ) related anomalies in the ionosphere. In this paper, GPS based Total Electron Content (TEC) data are studied for three major M?>?7.0 EQs in Nepal and Iran-Iraq border during 2015–2017 by implementing statistical procedures on temporal and spatial scale. Previous studies presented different time interval of pre-seismic ionospheric anomalies, however, this study showed that EQs ionospheric precursors may occur within 10?days. Furthermore, the ionospheric anomalies on the suspected day occurred during UT?=?08:00–12:00?h before the main shock. The Global Ionospheric Map TEC (GIM-TEC) data retrieved over the epicenter of M7.8 (Nepal EQ) showed a significant increase of 6 TECU on April 24, 2015 (one day before the main shock), which is recorded by the ground GPS station data of Islamabad (station lies within the EQ preparation zone). Furthermore, the spatial GIM-TEC result imply significant anomalies over the epicenter during the time interval between UT?=?08:00–12:00?h (LT?=?13:00–17:00). For M7.3 (Nepal EQ), the TEC anomalies were detected on May 10, 2015 (2?days before the event) in the temporal data. The spatial TEC data imply the huge clouds over the epicenter at about UT?=?08:00–12:00?h on May 10, 2015, that may be associated with this EQ in the quiet geomagnetic storm conditions. Similarly, temporal and spatial TEC showed anomaly on November 3, 2017, during UT?=?08:00–12:00 (9?days before the Iran-Iraq border EQ) after implementing the statistical method on it. Conversely, there exists a short-term but low magnitude TEC anomaly synchronized with a geomagnetic storm on November 7–8, 2017 (4 to 5?days prior to M7.3 Iran-Iraq border EQ). The diurnal and hourly GIM-TEC and VTEC data also imply the execution of ionospheric anomalies within 10?days prior to all events. All these positive anomalies in TEC may be due to the existence of a huge energy from the epicenter during the EQ preparation period.     

The ionospheric response to high-intensity long duration continuous AE activity (HILDCAA) event (13–15 April 2005) over mid-latitude African region

The ionospheric responses to High-Intensity Long Duration Continuous Auroral Electrojet Activity (HILDCAA) event which happened following the CIR-driven storm were studied over the southern hemisphere mid-latitude in the African sector. The 13–15 April 2005 event was analysed to understand some of the mechanisms responsible for the ionospheric changes during HILDCAA event. The ionosonde critical frequency of F2 layer (foF2) and Global Navigation Satellite System (GNSS) Total Electron Content (TEC) were used to analyse the ionospheric responses. The daytime increase in foF2 and TEC values were observed on 13 April 2005. The TEC and foF2 enhancement could be attributed to Large Scale Traveling Ionospheric Disturbances (LSTIDs), increase in thermospheric neutral composition changes, Prompt Penetration Electric Field (PPEF) and an expansion of Equatorial Ionization Anomaly (EIA) to the mid-latitude.     

磁暴期间中国中低纬电离层不规则体与扰动分析

2017年9月8日发生了一次强磁暴，Kp指数最大值达到8.利用区域电离层格网模型（Regional Ionosphere Map，RIM）和区域ROTI（Rate of TEC Index）地图，分析了磁暴期间中国及其周边地区电离层TEC扰动特征和低纬地区电离层不规则体的产生与发展情况，同时利用不同纬度IGS（International GNSS Service）测站BJFS（39.6°N，115.9°E），JFNG（30.5°N，114.5°E）和HKWS（22.4°N，114.3°E）的GPS双频观测值，获取各测站的ROTI和DROT（Standard Deviation of Differential ROT）指数变化趋势.结果表明:此次磁暴发生期间电离层扰动先以正相扰动为主，主要发生在中低纬区域，dTEC（differential TEC）最大值达到14.9TECU，随后电离层正相扰动逐渐衰减，在低纬区域发生电离层负相扰动，dTEC最小值达到-7.2TECU；在12:30UT-13:30UT时段，中国南部低纬地区发生明显的电离层不规则体事件；相比BJFS和JFNG两个测站，位于低纬的HKWS测站的ROTI和DROT指数变化更为剧烈，这表明电离层不规则体结构存在纬度差异.      

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.     

South American regional ionospheric maps computed by GESA: A pilot service in the framework of SIRGAS

SIRGAS (Geocentric Reference Frame for the Americas) is an international enterprise of the geodetic community that aims to realize the Terrestrial Reference Frame in the America’s countries. In order to fulfill this commitment, SIRGAS manages a network of continuously operational GNSS receivers totalling around one hundred sites in the Caribbean, Central, and South American region. Although the network was not planed for ionospheric studies, its potential to be used for such a purpose was recently recognized and SIRGAS started a pilot experiment devoted to establish a regular service for computing and releasing regional vertical TEC (vTEC) maps based on GNSS data. Since July, 2005, the GESA (Geodesia Espacial y Aeronomía) laboratory belonging to the Facultad de Ciencias Astronómicas y Geofísicas of the Universidad Nacional de La Plata computes hourly maps of vertical Total Electron Content (vTEC) in the framework of the SIRGAS pilot experiment. These maps exploit all the GNSS data available in the South American region and are computed with the LPIM (La Plata Ionospheric Model). LPIM implements a de-biasing procedure that improves data calibration in relation to other procedures commonly used for such purposes. After calibration, slant TEC measurements are converted to vertical and mapped using local-time and modip latitude. The use of modip latitude smoothed the spatial variability of vTEC, especially in the South American low latitude region and hence allows for a better vTEC interpolation. This contribution summarizes the results obtained by GESA in the framework of the SIRGAS pilot experiment.