Comparison of standard TEC models with a Neural Network based TEC model using multistation GPS TEC around the northern crest of Equatorial Ionization Anomaly in the Indian longitude sector during the low and moderate solar activity levels of the 24th solar cycle

The highest Total Electron Content (TEC) values in the world normally occur at Equatorial Ionization Anomaly (EIA) region resulting in largest ionospheric range delay values observed for any potential Space Based Augmentation System (SBAS). Reliable forecasting of TEC is crucial for satellite based navigation systems. The day to day variability of the location of the anomaly peak and its intensity is very large. This imposes severe limitations on the applicability of commonly used ionospheric models to the low latitude regions. It is necessary to generate a mathematical ionospheric forecasting and mapping model for TEC based on physical ionospheric influencing parameters. A model, IRPE-TEC, has been developed based on real time low latitude total electron content data using GPS measurements from a number of stations situated around the northern crest of the EIA during 2007 through 2011 to predict the vertical TEC values during the low and moderate solar activity levels of the 24th solar cycle. This model is compared with standard ionospheric models like International Reference Ionosphere (IRI) and Parameterized Ionospheric Model (PIM) to establish its applicability in the equatorial region for accurate predictions.     

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

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

Wavelet modelling in support of IRI

The slant total electron content (STEC) of the ionosphere is defined as the integral of the electron density along the ray-path of the signal between the transmitter and the receiver. So-called geometry free GPS measurements provide information on the electron density, which is basically a four-dimensional function depending on spatial position and time. Since ground-based measurements are not very sensitive to the vertical structure within the atmosphere, the ionosphere is often represented by a spherical layer, where all electrons are concentrated. Then the STEC is transformed into the vertical total electron content (VTEC), which is a three-dimensional function depending on longitude, latitude and time.In our approach, we decompose an ionospheric function, i.e. the electron density or the VTEC, into a reference part computed from a given model like the International Reference Ionosphere (IRI) and an unknown correction term expanded in a multi-dimensional series in terms of localizing base functions. The corresponding series coefficients are calculable from GPS measurements applying parameter estimation procedures. Since the GPS receivers are located rather unbalanced, finer structures are modelable just in regions with a sufficient number of observation sites. Due to the localizing feature of B-spline functions we apply a tensor product spline expansion to model the correction term regionally. Furthermore, the multi-resolution representation derived from wavelet analysis allows monitoring the ionosphere at different resolutions levels. We demonstrate the advantages of this procedure by representing a simulated VTEC data set over South America.     

Klobuchar电离层模型误差分析及预测

电离层延迟误差是全球导航卫星系统（global navigation satellite system，GNSS）中的重要误差源之一。目前在电离层延迟改正模型中，应用最广泛的是Klobuchar参数模型，但是该模型的改正率仅能达到60%左右，无法满足日益增长的精度需求。将国际GNSS监测评估系统（international GNSS monitoring & assessment system，iGMAS）发布的高精度电离层格网数据作为对照，对Klobuchar电离层模型误差进行计算和分析，结果发现在中纬度区域误差存在明显的周期性特征。为进一步提高Klobuchar电离层模型在中纬度区域的改正精度，建立了基于粒子群优化反向传播（back propagation，BP）神经网络的Klobuchar电离层误差预测模型,并以2019年10月的采样数据为例进行误差预测。结果表明，用该模型对中纬度区域电离层延迟提供误差补偿，可将精度提高到90%左右。     

基于CIT技术的广域增强系统电离层延迟修正算法

提出了一种基于电离层层析成像(CIT)技术的广域增强系统(WAAS)电离层延迟 修正算法. 该算法利用模式基函数与截断奇异值分解正则化组合的方式, 实现WAAS单频用户的电离层延迟误差修正. 基于中国区域23个广域基准站 和10个用户站的仿真结果分析表明, 传统的网格算法和基于CIT技术的电离 层延迟修正算法的电离层延迟修正精度与太阳活动、昼夜变化及地磁纬度 之间均存在明显的相关性. 基于CIT技术的电离层延迟修正算法精度优于网格算法, CIT算法的平均误差与标准差相比网格算法均有较大幅度的下降.      

北斗三号电离层模型性能分析

根据北斗卫星导航系统三步走战略，北斗三号系统空中星座在2020年左右完成部署，信号的覆盖范围从亚太大部分地区扩展至全球.实际上自2018年底北斗三号系统已经可以初步提供全球服务.相应地，北斗系统的单频电离层模型从只适用于亚太大部分地区的BDSKlob模型发展到适用于全球范围的BDGIM模型.为分析北斗三号电离层模型的性能，以IGS GIM最终数据为参考，比较了2019年3月1-31日期间北斗BDGIM模型、BDSKlob模型以及GPSKlob模型在不同时段不同地区的改正效果.结果表明:北斗三号BDGIM模型在中国区域的改正率达75%以上，优于BDSKlob模型59%的改正率；在全球范围的改正率达60%，优于GPSKlob模型49%的改正率.除南半球高纬度地区外，BDGIM模型的电离层误差值均小于GPSKlob模型，其在垂直方向的电离层均方根误差值优于6TECU.      

Epochwise prediction of GPS single differenced ionospheric delays of formation flying spacecraft

Both single and dual frequency GPS relative navigation filters may benefit from proper predictions of single differenced ionospheric delays. In this article, the single differenced ionospheric delays of GPS observations are predicted for the GRACE formation during the switch manoeuvre.Two prediction methods are considered. The first is based on a Taylor expansion to first order of a mapping function that maps slant total electron content measurements to vertical total electron content estimates. The second method fits a shape profile through undifferenced ionospheric data available. It then raytraces through this profile to estimate the difference in total electron content along the path of the GPS signals.Continuously changing ionospheric conditions hamper the assessment of the quality of the predictions. Comparison of both methods shows that the raytracing method performs better. The difference of predictions and measurements generally shows a smaller RMS than the measurements alone. However, both methods suffer from a number of systematically unpredicted observations, which arise from small unaccounted differential variations in electron densities along the path of the GPS signals. These prediction methods perform better when spacecraft separation is small. Baselines considered here range from tens of kilometres down to several hundred metres. When smallest spacecraft separation occurs (0.4 km), the single differenced ionospheric delays exhibit RMS values of 0.0089 m. The first method shows a difference between measurements and predictions with an RMS of 0.0081 m. For the second method the difference RMS is found to be 0.0067 m.     

Assessment study of ionosphere correction model using single- and multi-shell algorithms approach over sub-Saharan African region

It is a known fact that ionosphere is the largest and the least predictable among the sources of error limiting the reliability and accuracy of Global Navigation Satellite Systems (GNSS) and its regional augmentation systems like Satellite Based Augmentation System (SBAS) in a safety-of-life application. The situation becomes worse in the Equatorial Ionization Anomaly (EIA) region, where the daytime ionization distribution is modified by the fountain effect that develops a crest of electron density at around ±15° to ±20° of the magnetic equator and a trough at the magnetic equator during the local noon hours. Related to this phenomenon is the appearance of ionosphere irregularities and plasma bubbles after local sunset. These may degrade further the quality of service obtained from the GNSS/SBAS system of the said periods. Considering the present operational augmentation systems, the accuracy and integrity of the ionosphere corrections estimate decreases as the level of disturbances increases. In order to provide a correct ionosphere correction to the user of GNSS operating in African EIA region and meet the integrity requirements, a certified ionosphere correction model that accurately characterizes EIA gradient with the full capacity to over-bound the residual error will be needed. An irregularities detector and a decorrelation adaptor are essential in an algorithm usable for African sub-Saharan SBAS operation. The algorithm should be able to cater to the equatorial plasma vertical drifts, diurnal and seasonal variability of the ionosphere electron density and also should take into account the large spatial and temporal gradients in the region. This study presents the assessment of the ionosphere threat model with single and multi-layer algorithm, using modified planar fit and Kriging approaches.     

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.     

Computerized ionospheric tomography with the IRI model

Computerized ionospheric tomography (CIT) is a method to estimate ionospheric electron density distribution by using the global positioning system (GPS) signals recorded by the GPS receivers. Ionospheric electron density is a function of latitude, longitude, height and time. A general approach in CIT is to represent the ionosphere as a linear combination of basis functions. In this study, the model of the ionosphere is obtained from the IRI in latitude and height only. The goal is to determine the best representing basis function from the set of Squeezed Legendre polynomials, truncated Legendre polynomials, Haar Wavelets and singular value decomposition (SVD). The reconstruction algorithms used in this study can be listed as total least squares (TLS), regularized least squares, algebraic reconstruction technique (ART) and a hybrid algorithm where the reconstruction from the TLS algorithm is used as the initial estimate for the ART. The error performance of the reconstruction algorithms are compared with respect to the electron density generated by the IRI-2001 model. In the investigated scenario, the measurements are obtained from the IRI-2001 as the line integral of the electron density profiles, imitating the total electron content estimated from GPS measurements. It has been observed that the minimum error between the reconstructed and model ionospheres depends on both the reconstruction algorithm and the basis functions where the best results have been obtained for the basis functions from the model itself through SVD.     

电离层薄层高度对电离层模型化的影响

利用IRI2012模型分析了电离层薄层高度的时空变化规律,提出了基于应用中STEC的电离层改正误差分析理论,分析了电离层薄层高度变化的相关影响.结果表明,电离层薄层高度变化对电离层穿刺点位置、投影映射函数值、电离层建模结果、电离层模型精化和电离层模型精度评估结果的影响较大.高度截止角为10°时,电离层薄层高度变化导致电离层穿刺点的经纬度差异最大可达3.2°,投影映射函数最高可引入约15.46%的误差,电离层建模结果差异和建模实用误差最高分别达9.71%,3.64%,采用不同薄层高度数据的电离层模型参数拟合和模型精化结果最大可引入约9.26%的误差,采用不同电离层薄层高度数据进行模型精度评定时最大可引入约9.62%的误差.根据这些研究结果可知:在实际应用中应采用电离层薄层高度模型,并选取较大的卫星高度截止角来减小薄层高度变化引入的误差;采用固定高度时,区域电离层建模采用与实际电离层薄层一致的固定高度;进行精度评估时,参考数据的电离层薄层高度与需要精度评估的电离层模型薄层高度相等.      

WAAS系统中电离层折射校正的新方法及计算结果

电离层介质的色散性是影响电磁波信号进行卫星导航定位精度的重要因素之一.配合北斗二代分系绩研制任务,提出了一种新的电离层折射校正算法,并利用2000年7月1日到3日的双频GPS观测数据对6个用户站进行试算,进一步将试算所得均方根误差和电离层网格算法得到的误差进行比较.结果表明,对于中纬区域的用户站,估算的TEC误差约为0.5 m左右;而低纬用户误差相对增大,为1 m左右.文中给出的算法与电离层网格模型所提供的精度相差不多,在未来中国自主的卫星增强系统中采用新方法进行电离层进行修正是可行的及有效的.      

中国区域电离层异常数据野值检测

针对中国区域连续运行参考站接收机野值会干扰星基增强系统（SBAS）电离层异常事件提取的问题，提出了一种基于电离层垂直延迟时间梯度的野值检测方法。首先，介绍了电离层延迟数据的提取方法；然后，论述了依据野值和电离层异常的不同时空相关特性进行野值检测的方法，并用单双频定位误差结果验证了野值检测的正确性；最后，对检测结果进行了分析。结果表明:该方法可有效区分野值和电离层异常；在中国大陆构造环境监测网络200多个参考站中，野值检测所剔除的参考站数目在3~10个，对电离层穿透点空间分布的影响在可接受范围内。      

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.     

NeQuick G model based scale factor determination for using SBAS ionosphere corrections at low earth orbit

A space-based augmentation system (SBAS) provides real-time correction data for global navigation satellite system (GNSS) users near ground. In order to use the SBAS ionosphere correction for low Earth orbit (LEO) satellites, the correction should be scaled down for the LEO altitude. This scale factor varies with ionosphere distribution and it is hard to determine the value at LEO in real time. We propose a real-time scale factor determination method by using Galileo GNSS’s NeQuick G model. A LEO satellite GPS data and SBAS data received on ground were used to evaluate the performance of the NeQuick G derived variable scale factor. The NeQuick G derived scale factor shows a significant accuracy improvement over NeQuick G model or pre-determined constant scale factor. It improves a vertical positioning accuracy of the LEO satellite. The error mean reductions of the vertical positioning over NeQuick G and the constant scale factor are 31.5% and 11.7%, respectively.     

Effects of X2-class solar flare events on ionospheric GPS-TEC and radio waves over Brazilian sector

In this investigation, we present and discuss the effects of 6 X2-class solar flare events in the ionospheric F region over Brazilian sector that occurred during 2013 to 2015. For this investigation, we present the vertical total electron content (VTEC) observations from nearly 120 Global Positioning System (GPS) receivers all over the Brazilian sector for each event. Also, ionospheric sounding observations obtained in São José dos Campos (23.2°S, 45.9°W, dip latitude 17.6°S; hereafter referred to as SJC), under the southern crest of the equatorial ionospheric anomaly (EIA), Brazil, are presented. The observations show that the greatest TEC impact occurs with the EUV fluxes increases lasting for more than one hour and when the solar active region is located close to the solar disc center. We present a detailed study of the efficiency of the EUV flux with wavelengths ranging from 0.1 to 190?nm for the F region ionization. The largest increase of ΔTEC occurs below the magnetic equator line, covering mainly the central, northeast, southeast and south regions, which includes the equatorial ionospheric anomaly (EIA) region. The ionograms show partial or total fade out in the echoes traces observed causing blackouts of radio signals of up to 60?min, which can have serious consequences to technological systems of public and private agencies around Brazilian sector. This study can help to better understand the effects of solar flares in the ionospheric F region.     

基于BP神经网络技术的电离层VTEC融合

随着电离层探测技术的不断发展，电离层观测资料也越来越多，只使用单一的观测资料会出现电离层反演精度不高的问题。为了提高电离层的反演精度，使用BP神经网络技术将地基反演和国际参考电离层(international reference ionosphere，IRI)模型的垂直总电子含量（vertical total electron content，VTEC）数据进行有效融合。在温带地区\[35°(N)~45°(N)，60°(E)~80°(E)\]进行电离层反演试验，结果表明基于BP神经网络技术的电离层数据融合和地基反演获得的电离层VTEC精度都比较高，但是基于BP神经网络的电离层数据融合反演精度比地基反演更高，所以基于BP神经网络技术的数据融合能够提高电离层的反演精度。     

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

Improving regional ionospheric TEC mapping based on RBF interpolation

Radial basis function (RBF) interpolation with multi-quadric is developed to perform ionospheric total electron content (TEC) mapping for the Chinese region between 15°N ~ 40°N and 100°E ~ 125°E. TEC measurements from the Centre for Orbit Determination in Europe (CODE) covering the solar maximum year 2011 are used to investigate the performance of the proposed RBF interpolation method. The differences between the RBF interpolated TEC and the CODE TEC are within 0.5 TECU and the root mean square error (RMSE) is very small when 49 data points are used. The maximum difference is ~5 TECU and the error is less than 1 TECU with 25 samples. Our study suggests that a random distribution of measurement points gives smaller RMSEs than a homogenous distribution when the number of sample points is low. The study indicates that RBF interpolation offers a powerful and reliable tool for ionospheric TEC mapping.     

Comparison of observed ionospheric vertical TEC over the sea in Indian region with IRI-2016 model

The vertical ionospheric TEC values obtained from GAGAN grid based ionospheric delay correction values over the sea in the Indian equatorial region have been compared with the corresponding values derived from the International Reference Ionosphere model, IRI-2016. The objective of this work is to study the deviation of the vertical TEC derived from the IRI model from ground truths over the sea for different conditions. This will serve the basic intention of assessing the candidature of the IRI model as an alternative ionospheric correction model in navigation receivers in terms of accuracy. We have chosen different solar activity periods, seasons, geomagnetic conditions, locations etc. for our comparison and analysis. The TEC values by the IRI-2016 were compared with the actual measured values for the given conditions and errors were obtained. The measured vertical TEC values at the ionospheric grid points were derived from the GAGAN broadcast ionospheric delay data and used as reference. The IRI model with standard internal functions was used in estimating the TEC at the same ionospheric grid points. The errors in the model derived values are statistically analysed. Broadly, the results show that, for the Indian sector over the sea, the IRI model performs better on quiet days in off equatorial regions, particularly in the northern region. The overall performance degrades for other conditions with the model generally underestimating the true TEC values and most severely in the equatorial region. The performance is worst in this region for the disturbed days of the equinoctial period. The comparison study is also done with the TEC data measured directly by dual frequency GPS receivers. The results were found to be in general agreement with those obtained by comparing the model with GAGAN broadcast data as reference. This study will be useful in considering the IRI-2016 model for real time estimates of TEC as an alternative to the current parametric model in a satellite navigation receiver in absence of other options.