Performance of GPS slant total electron content and IRI-Plas-STEC for days with ionospheric disturbance

Total Electron Content （TEC） is an important observable parameter of the ionosphere which forms the main source of error for space based navigation and positioning systems. Since the deployment of Global Navigation Satellite Systems （GNSS）, cost-effective esti- mation of TEC between the earth based receiver and Global Positioning System （GPS） sat- ellites became the major means of investigation of local and regional disturbance for earthquake precursor and augmentation system studies. International Reference Iono- sphere （IRI） extended to plasmasphere （IRI-Plas） is the most developed ionospheric and plasmaspheric climatic model that provides hourly, monthly median of electron density distribution globally. Recently, IONOLAB group {www.ionolab.org） has presented a new online space weather service that can compute slant TEC （STEC） on a desired ray path for a given date and time using IRI-Plas model （IRI-Plas-STEC）. In this study, the performance of the model based STEC is compared with GPS-STEC computed according to the estimation method developed by the IONOLAB group and includes the receiver bias as IONOLAB-BIAS （IONOLAB-STEC）. Using Symmetric Kullback-Leibler Distance （SKLD）, Cross Correlation （CC） coefficient and the metric norm （L2N） to compare IRI-Plas-STEC and IONOLAB-STEC for the month of October 2011 over the Turkish National Permanent GPS Network （TNPGN- Active）, it has been observed that SKLD provides a good indicator of disturbance for both earthquakes and geomagnetic storms.

Longitudinal dependence of the forecast accuracy of the ionospheric total electron content in the equatorial zone

The longitudinal dependence of the behavior of ionospheric parameters has been the subject of a number of works where significant variations are discovered.This also applies to the prediction of the ionospheric total electron content(TEC),which neural network methods have recently been widely used.However,the results are mainly presented for a limited set of meridians.This paper examines the longitudinal dependence of the TEC forecast accuracy in the equatorial zone.In this case,the methods are used that provided the best accuracy on three meridians:European(30°E),Southeastern(110°E)and American(75°W).Results for the stations considered are analyzed as a function of longitude using the Jet Propulsion Laboratory Global Ionosphere Map(JPL GIM)for 2015.These results are for 2 h ahead and 24 h ahead forecast.It was found that in this case,based on the metric values,three groups of architectures can be distinguished.The first group included long short-term memory(LSTM),gated recurrent unit(GRU),and temporal convolutional networks(TCN)models as a part of unidirectional deep learning models;the second group is based on the recurrent models from the first group,which were supplemented with a bidirectional algorithm,increasing the TEC forecasting accuracy by 2-3 times.The third group,which includes the bidirectional TCN architecture(BiTCN),provided the highest accuracy.For this architecture,according to data obtained for 9 equatorial stations,practical independence of the TEC prediction accuracy from longitude was observed under the following metrics(Mean Absolute Error MAE,Root Mean Square Error RMSE,Mean Absolute Percentage Error MAPE):MAE(2 h)is 0.2 TECU approximately;MAE(24 h)is 0.4 TECU approximately;RMSE(2 h)is less than 0.5 TECU except Niue station(RMSE(2 h)is 1 TECU approximately);RMSE(24 h)is in the range of 1.0-1.7 TECU;MAPE(2 h)<1%except Darwin station,MAPE(24 h)<2%.This result was confirmed by data from additional 5 stations that formed latitudinal chains in the equatorial part of the three meridians.The complete correspondence of the observational and predicted TEC values is illustrated using several stations for disturbed conditions on December 19-22,2015,which included the strongest magnetic storm in the second half of the year(min Dst=-155 nT).

Ionospheric total electron content disturbance associated with May 12, 2008,Wenchuan earthquake

Possible ionospheric disturbances relating to the May 12, 2008, MsS.0 Wenchuan earthquake were identified by Global Positioning System （GPS）-derived total electron content （TEC）, ion- osonde observations, the global ionospheric map （GIM）, and electron density profiles detected by the Constellation Observation System for Meteorology Ionosphere and Climate （COSMIC）. We applied a statistical test to detect anomalous TEC signals and found that a unique enhancement in TEC, recorded at 16 GPS stations, appeared on May 9, 2008. The critical fre- quency at F2 peak （foF2）, observed by the Chinese ionosondes, and maximal plasma frequency, derived from COSMIC data, revealed a characteristic similar to GPS TEC variations. The GIM showed that the anomalous variations of May 9 were located southeast of the epicenter. Using GPS data from 13 stations near the epicenter, we analyzed the TEC variations of satellite orbit traces during 04：00-11：00 UT. We found that TEC decreased to the east and increased to the southeast of the epicenter during this period. Results showed that the abnormal disturbance on May 9 was probably an ionosphenc precursor of the Wenchuan earthquake of May 12, 2008.

Manifestations of the earthquake preparations in the ionosphere total electron content variations

Recent ionospheric observations report anomalous total electron content (TEC) deviations prior strong earthquakes. We discuss common fetures of the pre-earthquake TEC disturbances on the basis of statistics covering 50 strong seismic events during 2005-2006. The F2-layer ionospheric plasma drift under action of the electric fields of seismic origin is proposed as the main reason of producing TEC anomalies. The origin of such electric fields is discussed in terms of the lithosphere-atmosphere-ionosphere coupling system. This theory is supported by numerical simulations using global Upper Atmosphere Model (UAM). UAM calculations show that the vertical electric current with the density of about 20 - 40 nA/m2 flowing between the Earth and ionosphere over an area of about 200 by 2000 km is required to produce the TEC disturbances with the amplitude of about 30% - 50% relatively to the non-disturbed conditions. Ionosphere responses on the variations of the latitudinal position, direction and configuration of the vertical electric currents have been investigated. We show that not only the vertical component of the ionospheric plasma drift but also horizontal components play an important role in producing pre-earthquake TEC disturbances.

Ionospheric vertical total electron content prediction model in low-latitude regions based on long short-term memory neural network

Ionosphere delay is one of the main sources of noise affecting global navigation satellite systems, operation of radio detection and ranging systems and very-long-baseline-interferometry. One of the most important and common methods to reduce this phase delay is to establish accurate nowcasting and forecasting ionospheric total electron content models. For forecasting models, compared to mid-to-high latitudes, at low latitudes, an active ionosphere leads to extreme differences between long-term prediction models and the actual state of the ionosphere. To solve the problem of low accuracy for long-term prediction models at low latitudes, this article provides a low-latitude, long-term ionospheric prediction model based on a multi-input-multi-output, long-short-term memory neural network. To verify the feasibility of the model, we first made predictions of the vertical total electron content data 24 and 48 hours in advance for each day of July 2020 and then compared both the predictions corresponding to a given day, for all days. Furthermore, in the model modification part, we selected historical data from June 2020 for the validation set, determined a large offset from the results that were predicted to be active, and used the ratio of the mean absolute error of the detected results to that of the predicted results as a correction coefficient to modify our multi-input-multi-output long short-term memory model. The average root mean square error of the 24-hour-advance predictions of our modified model was 4.4 TECU, which was lower and better than5.1 TECU of the multi-input-multi-output, long short-term memory model and 5.9 TECU of the IRI-2016 model.

Questions on Optimization of Measurement of Total Electron Content in Ionosphere with GPS

The analysis of existing method for calculation of total content of electrons (TEC) in ionosphere using GPS occultation method does show that due to different values of signal/noise ration in GPS signals ?and , the new method of optimum measurements of relevant frequency components of TEC measured by phase and code methods should be developed. The optimum quantity of measurements of the above-mentioned frequency components is determined taking into account the limitation imposed on general number of necessary measurements.

Total Electron Content Diurnal and Seasonal Variations and Response to Solar Events at Koudougou Station in Burkina Faso

In this paper, we studied the seasonal behavior of the total electron content (TEC) during a part of solar cycle 24 ascending, maximum and decreasing phases at Koudougou station (Latitude: 12°15'09"N Longitude: 2°21'45"W). Response of TEC to solar recurrent events is presented. The highest values of the TEC in 2014, 2015 and 2016 were recorded on March and October, while in 2013 they were recorded on April and November, corresponding to equinox months. This observation shows that TEC values at the equinoxes are higher than those of solstices. Moreover, the monthly TEC varies in phase with the sunspots number showing a linear dependence of the TEC on solar activity. The ionospheric electron contents are generally very low both before noon and during the night, but quite high at noon and after noon. This pattern of TEC variation is due to the fluctuation of incident solar radiation on the Earth’s equatorial ionosphere. During quiet periods, the number of free electrons generated is lower than that generated during recurrent periods, which shows a positive contribution of recurrent activity to the level of the TEC. Investigations have also highlighted a winter anomaly and equinoctial asymmetry in TEC behavior at Koudougou station.

Total Electron Content during Recurrent and Quiet Geomagnetic Periods at the Koudougou Station in Burkina Faso

In this work, the comparative study of total electron content (TEC) between recurrent and quiet geomagnetic periods of solar cycle 24 at Koudougou station with geographical coordinates 12°15'N;- 2°20'E was addressed. This study aims to analyze how geomagnetic variations influence the behavior of TEC in this specific region. The geomagnetic indices Kp and Dst were used to select quiet and recurrent days. Statistical analysis was used to interpret the graphs. The results show that the mean diurnal TEC has a minimum before dawn (around 0500 UT) and reaches a maximum value around 1400 UT, progressively decreasing after sunset. In comparison, the average diurnal TEC on recurrent days is slightly higher than on quiet days, with an average difference of 7 TECU. This difference increases with the level of geomagnetic disturbance, reaching 21 TECU during a moderate storm. The study also reveals significant monthly variations, with March and October showing the highest TEC values for quiet and recurrent days, respectively. Equinox months show the highest mean values, while solstice months show the lowest. Signatures of semi-annual, winter and equatorial ionization anomalies were observed. When analyzing annual variations, it was found that the TEC variation depends significantly on F10.7 solar flux, explaining up to 98% during recurrent geomagnetic activity and 92% during quiet geomagnetic activity.

电离层TEC组合预报方法研究

针对电离层总电子含量(TEC)的时间序列存在非线性、非平稳性以及高噪声等特征的情况,本文建立奇异谱分析(SSA)与Holt-Winters模型的电离层TEC组合预报模型。使用本文提出模型与对比模型分别对地磁平静期与地磁活跃期的电离层TEC进行预报与分析,结果表明,在地磁平静期,本文提出组合模型的电离层TEC预报精度更高,预报残差在1 TECu以内;在地磁活跃期,本文提出组合模型预报精度较单一Holt-Winters加法模型有显著提升,预报残差值同样在1 TECu以内。通过对比实验与分析得出,本文提出的组合预报模型可有效提升电离层TEC预报精度。

基于注意力机制LSTM的电离层TEC预测

电离层总电子含量(Total Electron Content,TEC)的监测与预报是空间环境研究的重要内容,对卫星通讯和导航定位等有重要意义.TEC值影响因素较多,很难确定精确物理模型来对其进行预测.本文设计了基于注意力机制的LSTM模型(Att-LSTM),采用过去24小时TEC观测数据对未来TEC进行预测.选择北半球东经100°上,每2.5°纬度选择一个位置,共计36个位置来验证本文提出模型的性能,并与主流的深度学习模型如DNN、RNN、LSTM进行对比实验.取得了如下成果:(1)在选定的36个地区未来2小时单点预测上,基于本文的Att-LSTM模型的TEC预测性能明显优于其他对比模型;(2)讨论了纬度对Att-LSTM预测未来2小时TEC值时性能的影响,发现在北纬0°到60°之间,Att-LSTM预测性能随着纬度的升高而略有降低,在北纬62.5°~87.5°之间,模型预测性能出现扰动,预测效果略差;(3)讨论了磁暴期和磁静期模型的预测性能,发现无论是磁暴期还是磁静期,本文模型预测性能均较好;(4)还讨论了对未来多时点预测效果,实验结果表明,本文所提出的模型对未来2、4个小时的预测拟合度R-Square均超过0.95,预测结果比较可靠,对未来6、8、10个小时预测拟合度最高为0.7934,预测拟合度R-Square下降迅速,预测结果不可靠.

基于GPS TEC的2022年1月15日汤加火山喷发激起的电离层行扰分析

2022年1月15日南太平洋汤加海底火山发生剧烈喷发,是近30年来最大规模的火山爆发,产生的强烈大气波动为开展火山喷发电离层扰动研究提供了难得的机会。本文利用GPS数据探测火山附近、新西兰、澳大利亚和中国地区的电离层扰动,从波形、频率、传播速度和时空分布等角度分析汤加火山喷发电离层行扰(TIDs)的特征,并利用电离层测高站、海平面监测站和大气压监测站的观测数据进一步验证TIDs的传播特征。研究结果发现,汤加火山喷发在其附近区域、新西兰、澳大利亚和中国地区引起了3类TIDs。在火山附近东、西、南、北方向上均探测到第一类TIDs, TIDs的传播速度为617~972 m/s,该类TIDs极有可能由火山喷发产生的声波引起。汤加火山喷发仅在火山附近东、西方向引起第二类TIDs,其传播速度分别为472 m/s和418 m/s,可能由声波传播过程中衍生的声重力波或者混合波引起,形成机理有待进一步研究。汤加火山喷发在新西兰、澳大利亚和中国地区引发了第三类TIDs,其传播速度为328~352 m/s,该类TIDs与Lamb波密切相关。

Analysis of Anomalous Enhancement in TEC and Electron Density in the China Region Prior to the 17 March 2015 Geomagnetic Storm Based on Ground and Space Observations

Total Electron Content(TEC)and electron density enhancement were observed on the day before 17 March 2015 great storm in the China Region.Observations from ground-and space-based instruments are used to investigate the temporal and spatial evolution of the pre-storm enhancement.TEC enhancement was observed from 24°N to 30°N after 10:00 UT at 105°E,110°E and 115°E longitudes on March 16.The maximum magnitude of TEC enhancement was more than 10 TECU and the maximal relative TEC enhancement exceeded 30%.Compared with geomagnetic quiet days,the electron density of Equatorial Ionization Anomaly(EIA)northern peak from Swarm A/C satellites on March 16 was larger and at higher latitudes.NmF2 enhanced during 11:30—21:00 UT at Shaoyang Station and increased by 200%at~16:00 UT.However,TEC and electron density enhancement were not accompanied by a significant change of hmF2.Most research has excluded some potential mechanisms as the main driving factors for storm-time density enhancements by establishing observational constraints.In this paper,we observed pre-storm enhancement in electron density at different altitudes and Equatorial Electrojet(EEJ)strength results derived from ground magnetometers observations suggest an enhanced eastward electric field from the E region probably played a significant role in this event.

基于神经网络的北京地区电离层TEC经验模型研究

利用2009—2019年间的全球电离层图(Global Ionospheric Map,GIM)产品,提出了一种基于神经网络技术的北京地区电离层总电子含量(Total Electron Content,TEC)经验模型建立方法。模型精度验证结果表明,本文提出的方法能够有效提高电离层TEC经验模型精度,相对Klobuchar模型和BDGIM模型精度分别提高了62%和21%。BJFS站的单频定位验证结果表明,本文建立的电离层TEC经验模型能够帮助单频用户有效提高定位精度,相对Klobuchar模型和BDGIM模型的三维定位精度分别提高了32%和7.5%。

基于SSA-PSO-LSTM模型的电离层TEC预报

受多种因素影响,电离层电子总含量(TEC)时间序列具有非线性、非平稳性特征,为提升长短期记忆(LSTM)神经网络模型在电离层TEC预报中的精度,本文在该神经网络模型的基础上,引入奇异谱分析(SSA)与粒子群优化(PSO)算法,构建了新的SSA-PSO-LSTM模型。一方面,利用了SSA对TEC时间序列进行数据预处理;另一方面,利用粒子群优化算法改进LSTM神经网络模型参数。选用欧洲地球参考框架(EUREF)提供的格网点电离层TEC时间序列数据进行实验,实验结果表明,在磁平静期与磁暴期,该组合模型的TEC预报结果均方根误差分别为0.28个总电子含量单位(TECu)、0.83个TECu,平均相对精度分别为96.35%、91.33%,均优于对比模型,验证了本文提出的组合预报模型的有效性与优越性。平均相对精度分别为96.35%、91.33%,均优于对比模型,验证了本文提出的组合预报模型的有效性与优越性。

基于滑窗初始点dSTEC定权的全球实时电离层综合

自2017年起,国际GNSS服务组织(IGS)部分电离层分析中心先后提供全球实时电离层校正服务。针对实际应用中单一分析中心实时电离层产品存在的可靠性差问题,本文提出了基于滑窗初始点差分斜向总电子含量(dSTEC)定权的全球实时电离层综合方法。以中国科学院(CAS)、法国空间研究中心(CNES)、西班牙加泰罗尼亚理工大学(UPC)及武汉大学(WHU)实时产品为基础,实现了实时全球电离层地图(RT-GIM)产品的例行综合。从电离层延迟误差修正、单频标准单点定位(SF-SPP)及单频精密单点定位(SF-PPP) 3个方面,分析2022年2月15日至2022年3月15日综合RT-GIM产品的应用性能。以IGS事后GIM为参考,本文综合RT-GIM与UPC综合RT-GIM对应的RMS为分别为3.30和3.20 TECU,二者精度相当;以定位残差95%分位数为统计量,与IGS事后GIM产品相比,综合RT-GIM对应的SF-SPP及SF-PPP精度分别低7.7%和4.9%;与北斗三号广播电离层模型BDGIM相比,综合RT-GIM对应的SF-SPP及SF-PPP精度分别提升了15.9%和9.5%。自2022年起,CAS综合RT-GIM产品已例行提交至IGS。

中国地区电离层TEC暴扰动研究

电离层总电子含量(TEC)是空间天气研究和监测预报的重要参量.本文引入了电离层TEC扰动指数DI,对青岛等6个台站的DI数据进行分析,选取DI>0.35(DI≤-0.30)作为正(负)相电离层TEC扰动的强度标准,并以连续6h及以上的DI满足该值来判定电离层TEC暴扰动事件.对电离层TEC暴扰动事件的统计分析表明,在地方时日落后至子夜前为发生高峰时段,正(负)相暴扰动事件平均持续时间约为10.9h(10.5h),正相暴发生率以冬季为多,夏季为少,而负相暴则以夏季略高.发现位于赤道异常驼峰区的广州站和位于高中纬度的海拉尔站比典型中纬地区的北京站电离层TEC暴扰动更易发生,且低纬地区以正相暴扰动为主.分析表明,约有70%的电离层TEC暴扰动伴随着有地磁扰动,但是电离层TEC暴扰动并不完全由地磁扰动所引起,强烈气象活动等局地环境因素也可能对电离层TEC暴扰动有着重要影响.

汶川Ms8.0地震前电离层TEC异常分析

基于中国地壳运动观测网络的GPS观测数据,解算得到了汶川Ms8.0地震前中国区域上空电离层TEC分布图。用两种方法的研究显示,在相同时间和空间电离层TEC存在显著的异常扰动,即震前1～6天(除5月9日电离层TEC出现异常增加外)时间内多次出现显著的TEC异常减少,而且电离层TEC的异常驼峰有向磁赤道漂移的趋势,由于距离发震时刻较近,排除其他因素后认为该异常可能源于汶川Ms8.0地震引起的电离层效应。

基于GPS技术实时监测2009年7月22日日全食期间长三角区域电离层TEC变化

2009年7月22日上午发生的日全食是21世纪持续时间最长的日全食,其全食带覆盖了中国中部的长江流域,为研究日全食对电离层的影响提供了一次难得的机会.为此本文通过卡尔曼滤波算法实现了实时求解TEC和GPS系统硬件延迟,为实时监测日全食期间电离层变化提供了绝对的电离层TEC.采用上海和浙江区域内GPS网的观测数据,建立了实时区域电离层延迟模型,进而计算出了实时的VTEC和TEC变化率.同时考虑太阳和地磁活动参数,综合上述方法详细分析和讨论了长三角区域在此次日全食期间的TEC变化的电离层异常现象.

震前电离层TEC异常探测方法研究

考虑到震前电离层TEC异常扰动探测方法(主要包括平均数法、中位数法、四分位距法等)的局限性,提出采用电离层TEC拟合后的剩余残差来识别TEC异常,即先对电离层TEC时间序列进行拟合,然后利用滑动时窗法进行异常识别。基于中国地壳运动观测网络GNSS观测资料解算的电离层TEC数据,采用上述方法对汶川8.0级大地震进行分析,结果显示该方法能够合理有效地探测到震前电离层TEC的异常扰动。

基于CODE GIM的震前电离层TEC异常分析

利用CODE GIM数据分析了2008年5月12日汶川地震和2010年4月14日青海玉树地震震前电离层TEC变化。以震前10天TEC中位数为基准,1.5倍标准差为误差限值,并考虑相应时期的太阳和地磁活动水平(Dst指数和Kp指数),对震前电离层异常的分布进行了分析。结果表明,两次地震震前几天内都有TEC异常出现。