地基单站GNSS的电离层VTEC高精度解算方法

利用IGS提供的双频GNSS观测数据，分析了Kalman方法解算电离层垂直总电子含量（Vertical Total Electron Content，VTEC）存在的问题，提出了Kriging-Kalman改进解算方法，并对两种方法解算的电离层VTEC进行分析和比较.结果表明:在低纬地区，当观测卫星数量发生改变时，Kalman方法解算的VTEC存在跳变异常，Kriging-Kalman方法解算的VTEC变化较为平稳，不存在跳变现象.对比分析耀斑期间两种方法解算VTEC的变化，发现Kalman方法解算的VTEC变化明显小于耀斑引起VTEC的增量；Kriging-Kalman方法解算结果与实际变化相一致.表明Kriging-Kalman方法计算精度更高，能够更精确计算耀斑等剧烈异常空间天气活动期间的VTEC及其变化，有利于电离层VTEC日常精确监测、研究和工程应用.      

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

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

GPS电离层反演方法研究及其在地震方面的应用

利用地基GPS数据计算了电离层单球壳模型穿刺点上的垂向总电子含量(VTEC), 根据VTEC和卫星及测站接收机的差分码间偏差DCB的不同时变特性采用了复弧法, 将VTEC作为局部变量(每30 min 一组, 可调), DCB作为一天的全局量进行解算. 在解算的过程中, 充分考虑VTEC的空间分布特性, 利用变异函数通过Kriging插值法建立电离层VTEC的二维格网模型, 并给出了卫星和接收机的差分码间偏差DCB. 通过与IGS结果的比较, 发现其结果可靠, 且时空分辨率和稳定性都有较大提高. 同时, 基于简化的三元样条插值基函数对电离层电子密度进行三维展开, 利用乘型代数重建技术MART算法构建了同批数据的四维层析成像结果, 获得了电离层电子密 度的四维分布. 其结果与CHAMP无线电掩星结果非常一致. 利用上面两种算法又分别对2008年5月长江三角洲地区地基GPS数据进行处理, 简要分析了该时段该地区上空电离层总电子含量和电子密度的变化情况及其对汶川地震的响应.      

太阳活动低年低纬地区VTEC 变化特性分析

利用福州台站(26.1°N, 119.3°E, 磁纬14.4°N)电离层闪烁与TEC监测仪2006-2010年的观测数据, 对该地区垂直总电子含量(VTEC)进行时间变化特性分析. 结果表明, 春秋冬三季的VTEC平均最高值出现在06:00UT, 夏季出现在08:00UT, 所有季节的平均最低值均出现在21:00UT; VTEC变化存在季节异常和弱冬季异常, 春秋季节高, 冬夏季节低, 夏季VTEC比冬季低且最大值出现时间延迟; VTEC在2006-2009年呈现下降的变化趋势, 2010年开始增强, 年际变化与太阳活动及地磁活动变化趋势具有较好的对应关系; VTEC变化与太阳活动存在很好的相关性, 相关系数达到0.5以上, 地磁活动则显示了弱相关的特性; F_10.7与VTEC的相关性随着每天Kp指数总值Σ_kp的增大而减小.      

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)纬度地区使用该模型可以更精确地描述该地区的电离层,帮助航空器实现更高精度的定位。     

一种用于测量高层大气风场的新型地基Fabry-Perot干涉仪

地表海拔大约250\,km高度处大气非常稀薄,目前被动光学观测是该层风场探测最有效手段.Fabry-Perot干涉仪(FPI)由于具有较高能量利用率及光谱分辨率等特点,是该层大气最有效的地基风场探测仪器之一.基于采用光路缩束系统及滤光片后置(标准具之后)方法研制的小型化FPI,2014年分别在河北廊坊(39.40°N,116.65°E)和山西岢岚(38.71°N,111.58°E)进行了地基观测试验,将观测数据与反演算法相结合得到高层大气风场数据,并将数据结果与美国大气研究中心两台FPI的风场数据进行了比较研究,在气辉整体辐射较弱的情况下得到岢岚站的风速平均反演偏差为11.8 m·s-1.      

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

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

星载SAR电离层探测及误差补偿

工作于低频波段的星载合成孔径雷达（SAR）信号会受到电离层的显著影响，因此在系统设计时必须考虑相应的补偿方法。基于此，首先开展了适用于低频大宽带SAR模式的电离层影响评估研究，建立了基于勒让德展开的五阶误差分析模型，可有效解决传统模型各阶次耦合问题。其次，针对背景电离层色散问题，开展了基于双频自聚焦算法的补偿研究，利用ALOS PALSAR数据进行了半物理仿真验证，电离层反演精度优于04TECU，可有效提升图像聚焦质量；针对法拉第旋转角误差，利用ALOS PALSAR回波散射矩阵信息开展了补偿研究，结果显示，相比官网提供的补偿参数，误差可进一步降低27%。上述基于回波数据本身的补偿研究，可避免第三方数据精度差、分辨率低、需地面接收机、传播路径不一致等问题，同时降低了载荷成本。最后，结合SAR高分辨率特性，基于回波的电离层反演可有效提高现有电离层探测能力，开展了PALSAR 垂测仪联合反演电子密度研究，其结果比仅垂测仪数据精度普遍提高了30%以上。研究成果可为未来低频星载SAR的系统设计提供技术支撑。     

基于多种观测手段的东亚低纬电离层不规则体事件分析

利用海南台站（19.5°N,109.1°E,dip:13.6°N）和磁赤道区的多种地基和天基观测数据,对2011年11月20日观测到的电离层不规则体事件进行了分析.海南台站VHF雷达、电离层闪烁和数字测高仪的综合观测结果表明,当天日落附近发生了强的电离层不规则体事件,主要表现为雷达羽和强闪烁的形态.结合磁赤道区GPS和C/NOFS卫星观测结果进行分析可知,海南台站日落附近出现的雷达羽和强闪烁与南海磁赤道区产生的主等离子体泡存在明显联系.      

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

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

基于小波分解与重构方法研究电离层偶发E层

电离层不规则体对卫星导航、通信、雷达系统等有重要影响.通过数值模拟及与实测数据的对比,论证基于小波分解与重构方法实现利用掩星数据反演电离层不规则体的可行性.以电离层偶发E层为例,利用国际参考电离层(IRI)模拟背景电离层电子密度分布,利用掩星探测的水平电子密度总含量δht反演不规则体信息,并与模拟数据进行比较.对200...     

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

The multi-source data fusion global ionospheric modeling software—IonoGim

We introduce a new global ionospheric modeling software—IonoGim, using ground-based GNSS data, the altimetry satellite and LEO (Low Earth Orbit) occultation data to establish the global ionospheric model. The software is programmed by C++ with fast computing speed and highly automatic degree, it is especially suitable for automatic ionosphere modeling. The global ionospheric model and DCBs obtained from IonoGim were compared with the CODE (Center for Orbit Determination in Europe) to verify its accuracy and reliability. The results show that IonoGim and CODE have good agreement with small difference, indicating that IonoGim owns high accuracy and reliability, and can be fully applicable for high-precision ionospheric research. In addition, through comparison between only using ground-based GNSS observations and multi-source data model, it can be demonstrated that the space-based ionospheric data effectively improve the model precision in marine areas where the ground-based GNSS tracking station lacks.     

The spatio-spectral localization approach to modeling VTEC over the western part of the USA using GPS observations

The precise ionosphere modeling is crucial and remains a challenge for GPS positioning and navigation, as well as many other Earth observation systems. In this research, two approaches have been proposed to model the vertical total electron content (VTEC) of the ionosphere using spherical slepian function. For the two-dimensional case, VTEC has been modeled in a sun-fixed reference frame and the three-dimensional approach based on the system of the three-dimensional base functions has been defined as the tensor product of the spherical slepian function for the longitude and latitude in an Earth-fixed reference frame, and the polynomial B-spline function for time. Rather than the spherical harmonics, the spherical slepian functions can be employed to produce the locally and globally orthogonal bases to optimally represent the data in any arbitrary region up to a given degree. The spherical slepian functions have been applied to the real data obtained from the ground-based GPS observation across the western part of the USA.     

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.     

基于因果卷积与LSTM网络的电离层总电子含量预报

电离层总电子含量（TEC）不仅是分析电离层形态的关键参数之一,同时为导航及定位等空间应用系统消除电离层附加时延提供重要支撑。由于电离层TEC的时空变化特征,本文融合因果卷积和长短时记忆网络,以太阳活动指数F_10.7、地磁活动指数Dst和电离层TEC历史数据作为特征输入,构建深度学习模型,实现提前24 h预报电离层TEC。进一步利用2005－2013年连续9年的CODE TEC数据,全面评估了模型在北京站（40°N,115°E）、武汉站（30.53°N,114.36°E）和海口站（20.02°N,110.38°E）的预报性能。结果显示不同太阳活动条件下三个站的TEC值与真实测量值的相关系数都大于0.87,均方根误差大都集中在0～1 TECU以内,且模型预报精度与纬度、太阳、地磁活动程度、季节变化相关。与仅由长短时记忆网络构成的预报模型相比,本实验模型均方根误差降低了15%,为电离层TEC预报模型的实际应用提供了参考。      

Changes in total electron content (TEC) during the annular solar eclipse of 15 January 2010

The solar eclipse of 15 January 2010 was an annular eclipse of the Sun with a maximum magnitude of 0.96 at 1.62°N, 69.29°E. To study the effect of this solar eclipse on the ionosphere the GPS data recorded at three different Indian stations Varanasi (Geographic latitude 25°, 16′N, longitude 82°, 59′E), Hyderabad (Geographic latitude 17°, 20′N, longitude 78°, 30′E) and Bengaluru (Geographic latitude 12°, 58′N, longitude 77°, 33′E) have been used to retrieve ionospheric total electron content (TEC). The ionospheric response to this rare event has been studied in terms of GPS-derived TEC observed at all the three Indian stations. A significant reduction in TEC reflected by all PRNs at all the three stations has been observed. The magnitude of the reduction in VTEC compared to quiet mean VTEC depends on latitude as well as longitude. The amount of reduction observed from different satellites (PRN) is different and depends on the location of the satellite from the solar eclipse path.     

一种基于系数择优的低阶球谐电离层延迟改正模型

针对单频接收机的电离层延迟改正问题, 提出了一种基于系数择优的低阶球谐电离层延迟改正模型. 按照电离层延迟改正模型参数择优问题的描述, 明确参数优化的目标和约束条件, 根据参数选择可编码的特点, 提出了利用遗传算法进行参数择优的方法及步骤. 以欧洲定轨中心(CODE)提供的电离层数据作为参考标准, 对参数择优模型、 低阶球谐模型和Klobuchar模型模拟的区域电离层VTEC精度进行了比较分析. 结果表明, 较之相同系数个数的低阶球谐模型, 参数择优模型精度平均改进了1～2TECU, 而且比Klobuchar模型及低阶球谐模型能更好地反映电离层的周日变化及纬度变化特征.      

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.