电离层周日变化对解算GPS硬件延迟稳定性的影响

针对电离层周日变化特征分析了其可能对SCORE方法估算的硬件延迟稳定性的影响. 利用BJFS以及XIAM台站的GPS观测数据, 解算了位于太阳活动高年(2001年)和太阳活动低年(2009年)的卫星硬件延迟并分析了估算的硬件延迟的稳定性. 研究发现, 电离层周日变化对估算的硬件延迟稳定性具有一定影响, 但是利用不同台站所得到的卫星硬件延迟稳定性在昼夜不同时间上的解算结果存在一定差异. 电离层周日变化对利用 BJFS台站数据解算的硬件延迟稳定性日夜差异较为明显, 在太阳活动高年利用XIAM 台站数据解算的硬件延迟日夜稳定性差异不很明显, 由于XIAM台站处于电离层赤道异常峰附近, 夜间电离层变化很大, 因此对比中纬度地区, 电离层周日变化对赤道异常峰附近地区硬件延迟稳定性解算结果的影响相对较小, 但在太阳活动低年, 其影响仍较为显著.      

利用GPS技术实时监测上海区域电离层变化

基于GPS技术实时监测电离层变化原理, 利用载波平滑伪距观测值建立区域电离层模型的方法, 计算了电离层延迟量和硬件延迟, 根据硬件延迟值相对稳定的特点, 采取一定时段求解出硬件延迟量, 对实时硬件延迟量进行预报, 进而实时分离GPS信号传播路径上的垂直总电子含量VTEC. 利用上海区域内的GPS网的观测数据, 建立实时上海区域电离层延迟模型, 监测上海区域的电离层变化. 数据分析结果表明, 这种方法的内符合精度优于3 TECU.      

夜间135.6 nm气辉反演的电离层f0F2与测高仪观测比较研究

在夜间电离层，气辉135.6 nm谱线主要由F层的O⁺和电子的辐射复合过程以及O⁺和O^–的中性复合过程激发，该谱线强度和电离层峰值电子密度Nm F₂存在很强的相关性。利用夜气辉135.6 nm辐射强度与F₂层峰值电子密度Nm F₂的平方成正比的物理模型，建立了在不同经纬度、地方时、季节和太阳活动下均适用的反演算法。通过DMSP卫星上搭载的紫外光谱成像仪(SSUSI)实际观测的135.6 nm气辉辐射强度来反演相应时空的电离层F₂层临界频率f₀F₂，并将其与地基测高仪探测结果做了综合对比。结果表明，在太阳活动高年(2013年)，相对误差小于等于20%的数据占比93.0%，平均相对误差约为7.08%；在太阳活动低年(2017年)，相对误差小于等于20%的数据占比80.8%，平均相对误差约为12.64%。最后，对该算法在太阳活动高低年的反演精度差异进行了分析。     

太阳活动对电离层TEC变化影响分析ormalsize

为研究太阳活动对电离层TEC变化的影响,从整体到局部分析了2000—2016年的太阳黑子数、太阳射电流量F_10.7指数日均值与电离层TEC的关系,并重点分析了2017年9月6日太阳爆发X9.3级特大耀斑前后15天太阳活动与电离层TEC变化的相关性.结果表明:由2000—2016年的数据整体看来,太阳黑子数、太阳F_10.7指数、TEC两两之间具有很强的整体相关性,但局部相关性强弱不均;此次耀斑爆发前后太阳黑子数、太阳F_10.7指数和TEC具有很强的正相关特性,太阳活动对TEC的影响时延约为2天;太阳活动对全球电离层TEC的影响不同步,从高纬至低纬约有1天的延迟,且对低纬度的影响远大于中高纬度.太阳活动是影响电离层TEC变化的主要原因,但局部也可能存在其他重要影响因素.      

玉树地震前的电离层异常现象分析

分析了玉树地震前地基电离层探测临界频率、GPS TEC和卫星探测原位等离子体参量等多个参数的扰动变化信息, 研究了不同高度异常变化的时空关联性. 分析发现, 在地震前一天的4月13日, 多个电离层参量出现同步扰动异常, 电离层临界频率f₀F₂异常相对滑动中值增大40%, 异常空间上存在从震中东南向西南漂移的特性; GPS TEC异常增强15TECU (1TECU=1016m-2)左右, 分布于震中南部经度15°范围内, 且有明显的磁共轭效应; DEMETER观测的原位氧离子密度N_i(O+) 4月13日为1-4月中最强的一天, 异常分布偏向赤道区, 但仅局限在30°-50°左右的经度范围内. 综合三个参量的异常特征发现, 无论是空间的局地性还是时间上的密切关联均反映这次电离层扰动可能与玉树地震孕育有关. 结合其他观测信息, 进一步探讨了这次地震孕育过程的地震电离层耦合机理.      

基于双频GPS观测的电离层TEC与硬件延迟反演方法

利用全球定位系统（Global Positioning System,GPS）的双频观测数据反演得到电离层的总电子含量（Total Electron Content,TEC）,使得广域甚至全球范围高时空分辨率的电离层观测研究成为可能,但由于GPS卫星和接收机对信号的硬件延迟可导致TEC测量系统偏差,因此,需要探索反演TEC并估测GPS卫星与接收机硬件延迟的有效算法.本文根据电离层电波传播理论,阐述了基于双频GPS观测提取电离层TEC的方法,给出TEC与硬件延迟的基本关系.综合研究了TEC与硬件延迟的反演方法,进行分析与归纳分类,在此基础上提出了有待深入研究的问题.      

基于最大似然估计的远紫外遥感反演电离层电子密度算法

原子氧135.6 nm夜气辉主要由氧离子O+与电子的辐射复合反应生成,一些星载远紫外遥感观测任务证实135.6 nm夜气辉可用于反演电离层电子密度。针对远紫外临边遥感观测反演电离层电子密度,分析了135.6 nm夜气辉辐射强度与电子密度之间的非线型前向模型,基于离散反演理论设计了从夜间135.6 nm临边观测数据反演电子密度高度分布的反演算法,算法应用最大似然估计通过迭代求解电离层参数的最佳拟合值。通过仿真计算了TIMED卫星上全球紫外成像仪GUVI观测的反演结果,验证了本反演算法的可行性。对GUVI的实际观测数据进行反演,获得了电子密度高度分布。通过与GUVI数据的电离层参数对比分析得出,本文建立的反演模型使N_mF₂被高估,同时使h_mF₂被低估。对于不同的太阳活动强度,N_mF₂和 h_mF₂的系统误差分别在10%和5%以内,能较精确地获得电离层参数。精确获得电离层电子密度信息对于提高空间天气预报及电离层模型的修正具有重要意义。      

多GNSS全球电离层建模特性及其精度检验

分别以GPS单系统和融合BDS,GPS,GLONASS三系统两种方案,采用载波相位平滑伪距观测值和球谐函数,构建了全球电离层延迟模型并进行了对比和分析.本文以GPS单系统和融合三系统两种方法反演了2014年1月每日电离层变化过程,解算得出了频间偏差的月综合产品,并对结果进行了对比和分析.事实上,三系统融合不仅增加了可观测的卫星数,而且改善了穿刺点的几何分布.分析结果表明,三系统融合反演全球电离层在精度上优于GPS单系统,均有5~10 TECU的提高.计算得到的频间偏差结果显示,GPS优于GIONASS,BDS稳定性则较次之.     

利用北斗观测数据实时监测中国区域电离层变化

北斗卫星导航信号采用三个频点工作,可以利用伪距双频组差方法解算电离层电子含量,为实时监视中国区域电离层变化提供新的技术手段.中国中低纬度处于电离层赤道异常变化区,在北纬20°±5°区域时常发生较大梯度的电离层变化.利用北斗实时多频伪距和相位观测数据,采用相位平滑伪距方法计算电离层穿刺点电子含量,分析通过北斗系统GEO卫星监测的电离层周日变化特性;采用多面函数方法拟合中国区域1°×1°分辨率的电离层延迟量,每5min绘制一幅中国区域电离层图,观测区域所有电离层穿刺点拟合残差RMS为2.778TECU;分析北斗系统实时监测中国区域电离层异常情况,当发生电离层异常变化时,相邻两天的VTEC（Vertical Total Electronic Content）峰值相差约60TECU.      

国际参考电离层模型输出参数峰值高度性能

电离层峰值高度H_mF₂是描述电离层形态的重要参数之一,国际参考电离层模型IRI-2016中融入了大量电离层测高仪和无线电掩星探测数据,用以提升H_mF₂的预测精度.本文利用太阳活动低年（2007—2010年）气象、电离层和气候卫星联合观测系统COSMIC探测数据描述全球范围内COSMIC H_mF₂的三维形态变化,对比分析了IRI-2016与IRI-2012模型的预测结果,同时分析了IRI-2016模型输出H_mF₂的性能.结果表明,IRI模型在中高纬度地区的输出结果高于COSMIC反演结果,而赤道及低纬地区则大都偏低.与IRI-2012模型相比,IRI-2016模型的输出结果在夜间至凌晨时段呈现较为明显的纬向梯度变化且大部分区域输出值偏高,但在白天时段赤道附近区域的输出值大都偏低.上述结果为电离层IRI模型的完善提供了一定参考.      

2009年7月22日日全食期间电离层总电子含量变化研究

利用中国区域内五个GPS台站(一个台站处于日全食区域、四个台站处于日偏食区域)观测数据, 研究2009年7月22日日全食期间电离层总电子含量(TEC)的变化, 结果表明, 日全食期间, 电离层TEC值经历了下降和恢复的过程, 最小TEC相对于最大食偏的时间延迟约为1～10min; 台站测得最小TEC的星下点(IPP)越靠近日全食带TEC下降量越大, 在日食期间武汉站(114.35°E, 30.53°N) TEC相对于各参考日期的TEC, 其平均下降量最大, 达到4.58TECU.      

大耀斑期间向日面电离层总电子含量的响应个例分析

利用2001年4月15日1336UT耀斑爆发期间向日面GPS观测数据提取的总电子含量的时间变化曲线。分析了向日面电离层对这次耀斑的响应特点．结果表明，耀斑期间向日面电离层出现了总电子含量突增事件．最大总电子含量增加量约为2．6TECU，在0600LT和1800LT都观测到了总电子含量突增，世增加幅度仅为0．5-1TECU．在高纬地区，由于电离层闪烁，从TEC时间变化曲线提取不出来总电子含量增加值．从各卫星星下点处的TEC增加量和各星下点处的太阳天顶角的关系可以看到，TEC增加量与太阳天顶角有关，太阳天顶角越大，TEC增幅越小。另外，从总电子含量时间变化率曲线上还观测到了时间同步的小尺度扰动，通过与耀斑期间硬X射线辐射通量的比较，发现两者有明显的相关性，电离层中的这种扰动与耀斑期间的硬X射线或远紫外辐射有关．     

Plasmaspheric electron content derived from GPS TEC and FORMOSAT-3/COSMIC measurements: Solar minimum condition

The plasmaspheric electron content (PEC) was estimated by comparison of GPS TEC observations and FORMOSAT-3/COSMIC radio occultation measurements at the extended solar minimum of cycle 23/24. Results are retrieved for different seasons (equinoxes and solstices) of the year 2009. COSMIC-derived electron density profiles were integrated up to the height of 700 km in order to retrieve estimates of ionospheric electron content (IEC). Global maps of monthly median values of COSMIC IEC were constructed by use of spherical harmonics expansion. The comparison between two independent measurements was performed by analysis of the global distribution of electron content estimates, as well as by selection specific points corresponded to mid-latitudes of Northern America, Europe, Asia and the Southern Hemisphere. The analysis found that both kinds of observations show rather similar diurnal behavior during all seasons, certainly with GPS TEC estimates larger than corresponded COSMIC IEC values. It was shown that during daytime both GPS TEC and COSMIC IEC values were generally lower at winter than in summer solstice practically over all specific points. The estimates of PEC (h > 700 km) were obtained as a difference between GPS TEC and COSMIC IEC values. Results of comparative study revealed that for mid-latitudinal points PEC estimates varied weakly with the time of a day and reached the value of several TECU for the condition of solar minimum. Percentage contribution of PEC to GPS TEC indicated the clear dependence from the time with maximal values (more than 50–60%) during night-time and lesser values (25–45%) during day-time.     

Application of BDS-GEO for studying TEC variability in equatorial ionosphere on different time scales

The Earth's ionosphere and especially its equatorial part is a highly dynamical medium. Geostationary satellites are known to be a powerful tool for ionospheric studies. Recent developments in BDS-GEO satellites allow such studies on the new level due to the best noise pattern in TEC estimations, which corresponds to those of GPS/GLONASS systems. Here we used BDS-GEO satellites to demonstrate their capability for studying equatorial ionosphere variability on different time scales. Analyzing data from the equatorial SIN1 IGS station we present seasonal variations in geostationary slant TEC for the periods of high (October 2013 - October 2014) and low (January 2017 - January 2018) solar activity, which show semi-annual periodicity with amplitudes about 10 TECU during solar maximum and about 5 TECU during the solar minimum. The 27-day variations are also prominent in geostationary slant TEC variations, which correlates quite well with the variations in solar extreme UV radiation. We found semi-annual pattern in small scale ionospheric disturbances evaluated based on geostationary ROTI index: maximal values correspond to spring and fall equinoxes and minimum values correspond to summer and winter solstices. The seasonal asymmetry in ROTI values was observed: spring equinox values were almost twice as higher than fall equinox ones. We also present results on the 2017 May 28–29 G3 geomagnetic storm, when ～30 TECU positive anomaly was recorded, minor and final major sudden stratospheric warmings in February and March 2016, with positive daytime TEC anomalies up to 15–20 TECU, as well as the 2017 September 6 X9.3 solar flare with 2 TECU/min TEC rate. Our results show the large potential of geostationary TEC estimations with BDS-GEO signals for continuous monitoring of space weather effects in low-latitude and equatorial ionosphere.     

GPS derived TEC and foF2 variability at an equatorial station and the performance of IRI-model

The ionosphere induces a time delay in transionospheric radio signals such as the Global Positioning System (GPS) signal. The Total Electron Content (TEC) is a key parameter in the mitigation of ionospheric effects on transionospheric signals. The delay in GPS signal induced by the ionosphere is proportional to TEC along the path from the GPS satellite to a receiver. The diurnal monthly and seasonal variations of ionospheric electron content were studied during the year 2010, a year of extreme solar minimum (F10.7 = 81 solar flux unit), with data from the GPS receiver and the Digisonde Portable Sounder (DPS) collocated at Ilorin (Geog. Lat. 8.50°N, Long. 4.50°E, dip −7.9°). The diurnal monthly variation shows steady increases in TEC and F2-layer critical frequency (foF2) from pre-dawn minimum to afternoon maximum and then decreases after sunset. TEC show significant seasonal variation during the daytime between 0900 and 1900 UT (LT = UT + 1 h) with a maximum during the March equinox (about 35 TECU) and minimum during the June solstice (about 24 TECU). The GPS-TEC and foF2 values reveal a weak seasonal anomaly and equinoctial asymmetry during the daytime. The variations observed find their explanations in the amount of solar radiation and neutral gas composition. The measured TEC and foF2 values were compared with last two versions of the International Reference Ionosphere (IRI-2007 and IRI-2012) model predictions using the NeQuick and CCIR (International Radio Consultative Committee) options respectively in the model. In general, the two models give foF2 close to the experimental values, whereas significant discrepancies are found in the predictions of TEC from the models especially during the daytime. The error in height dependent thickness parameter, daytime underestimation of equatorial drift and contributions of electrons from altitudes above 2000 km have been suggested as the possible causes.     

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.     

基于深度学习递归神经网络的电离层总电子含量经验预报模型

利用行星际太阳风参数与太阳活动指数、地磁活动指数、电离层总电子含量格点化地图数据,首次基于一种能处理时间序列的深度学习递归神经网络（Recurrent Neural Network,RNN）,建立提前24h的单站电离层TEC预报模型.对北京站（40°N,115°E）的预测结果显示,RNN对扰动电离层的预测误差低于反向传播神经网络（Back Propagation Neural Network,BPNN）0.49～1.46TECU,将太阳风参数加入预报因子模型后对电离层正暴预测准确率的提升可达16.8%.RNN对2001和2015年31个强电离层暴预报的均方根误差比BPNN低0.2TECU,将太阳风参数加入RNN模型可使31个事件的平均预报误差降低0.36～0.47TECU.研究结果表明深度递归神经网络比BPNN更适用于电离层TEC的短期预报,且在预报因子中加入太阳风数据对电离层正暴的预报效果有明显改善.      

Estimation of GPS instrumental biases from small scale network

With 4 GPS receivers located in the equatorial anomaly region in southeast China, this paper proposes a grid-based algorithm to determine the GPS satellites and receivers biases, and at the same time to derive the total electron content (TEC) with time resolution of 15 min and spatial resolution of 1° by 3.5° in latitude and longitude. By assuming that the TEC is identical at any point within a given grid block and the biases do not vary within a day, the algorithm arranges unknown biases and TECs with slant path TEC from the 4 receivers’ observations into a set of equations. Then the instrumental biases and the TECs are determined by using the least squares fitting technique. The performance of the method is examined by applying it to the GPS receiver chain observations selected from 16 geomagnetically quiet days in four seasons of 2006. It is found that the fitting agrees with the data very well, with goodness of fit ranging from 0.452 TECU to 1.914 TECU. Having a mean of 0.9 ns, the standard deviations for most of the GPS satellite biases are less than 1.0 ns for the 16 days. The GPS receiver biases are more stable than that of the GPS satellites. The standard deviation in the 4 receiver bias is from 0.370 ns to 0.855 ns, with a mean of 0.5 ns. Moreover, the instrumental biases are highly correlated with those derived from CODE and JPL with independent methods. The typical precision of the derived TEC is 5 TECU by a conservative estimation. These results indicate that the proposed algorithm is valid and qualified for small scale GPS network.     

Comparison of observed TEC values with IRI-2007 TEC and IRI-2007 TEC with optional foF2 measurements predictions at an equatorial region,Chumphon, Thailand

In this research, as part of working towards improving the IRI over equatorial region, the total electron content (TEC) derived from GPS measurements and IRI-2007 TEC predictions at Chumphon station (10.72°N, 99.37°E), Thailand, during 2004–2006 is analyzed. The seasonal variation of the IRI-2007 TEC predictions is compared with the TEC from the IRI-2007 TEC model with the option of the actual F2 plasma frequency (foF2) measurements as well as the TEC from the GPS and International GNSS service (IGS). The Chumphon station is located at the equatorial region and the low latitude of 3.22°N. For a declining phase of the solar cycle (2004–2006), the study shows that the IRI-2007 TEC underestimates the IRI-2007 TEC with the foF2 observation at the nighttime by about 5 TECU. The maximum differences are about 15 TECU during daytime and 5 TECU during nighttime. The overestimation is more evident at daytime than at nighttime. When compared in terms of the root-mean square error (RMSE), we find that the highest RMSE between GPS TEC and IRI 2007 TEC is 14.840 TECU at 1230 LT in 2004 and the lowest average between them is 1.318 TECU at 0630 LT in 2006. The noon bite-out phenomena are clearly seen in the IRI-2007 TEC with and without optional foF2 measurements, but not on the GPS TEC and IGS TEC. The IRI TEC with optional foF2 measurements gives the lowest RMSE values between IRI TEC predicted and TEC measurement. However, the TEC measurements (GPS TEC and IGS TEC) are more correct to use at Chumphon station.