基于机器学习集成算法的电离层层析算法迭代初值精化

GNSS电离层层析技术在电离层探测中扮演着重要角色,然而,由于观测数据不足或分布不均,导致层析模型中固有的不适定性问题成为制约该技术推广应用的主要瓶颈,其主要表现在无信号射线像素的反演电子密度对初值比较依赖,而初值通常是由精度不高的经验模型给出,从而拉低了电离层层析的整体精度.针对该问题,本文提出了一种基于机器学习集成算法(XGBoost)的电离层层析新方法(XGB-CIT),即基于传统层析算法,利用多个连续时段中有信号射线像素的特征及其电子密度反演值对机器学习模型进行训练,然后预测后续时段中所有像素的电子密度,并以此作为电离层层析的迭代初值,实现对电离层层析算法迭代初值的精化以及层析精度和效率的提升.利用湖南省连续两天共46个时段的CORS观测数据进行层析反演,并以23个连续时段为滑动窗口构建机器学习模型进行初值预测,并在此基础上利用正演误差和电离层测高仪数据对XGB-CIT的精度和适用性进行检验.其结果表明：相较于IRI2016模型,XGBoost算法提供的迭代初值精度提高了68%,而基于该初值得到的XGB-CIT模型,其精度和效率也比传统层析方法有所提高,其中收敛速度提高了20%...     

电离层三维层析成像的自适应联合迭代重构算法

在电离层层析成像过程中,联合迭代重构算法是一种常用的反演算法.然而,该算法迭代收敛较慢,反演结果精度不高.为此,本文发展了一种自适应的联合迭代重构算法,该算法利用上一轮的电离层电子密度反演结果,自适应地调整松弛因子和加权参数.通过模拟数据和实测数据对该算法的反演结果进行了验证,并将得到的反演结果与电离层测高仪数据进行了比较,结果表明,该算法能够有效地反演电离层电子密度,且反演结果精度优于常用的联合迭代重构算法.     

电离层TEC卡尔曼滤波成像研究

随着太空探测技术的进步,对TEC(Total Electron Content,简称TEC)探测精度要求越来越高.本文利用COSMOS 2414卫星数据资料获得观测TEC,在电离层NeQuick模型下,得到电离层电子密度,并使用卡尔曼滤波算法反演电子密度,最后结合电离层测高仪数据对实验结果进行判定.结果发现利用卡尔曼滤波反演信标资料算法,可以获得可靠的二维电子密度场.     

电离层层析重构的一种新算法

自适应联合迭代重构算法是电离层层析成像中一种收敛速度较快的反演算法.然而,在电离层电子密度重构过程中,该方法对迭代初值的精度要求较高.针对上述问题,本文提出了一种约束自适应联合迭代重构算法,该方法通过附加合理的平滑约束,减弱了没有观测信息的格网对迭代初值精度的依赖性,有效地提高了反演结果精度.数值模拟实验和实测数据的反演结果证实了该算法的可行性和在重构精度上的优越性.     

电离层GPS掩星反演技术研究

GPS无线电掩星技术是崭新的、高效的地球大气层和电离层探测技术,但仍在发展和完善之中.本文详细推导了Abel积分和绝对TEC电离层反演方法,研究了如何解决Abel积分产生的上下限异常问题;用COSMIC发布的GPS原始数据进行了反演计算,将结果与地面电离层测高仪数据进行了比较,最后讨论了周跳对反演结果的影响问题.结果表明：（1）在较高轨道高度（约800 km）,Abel积分与绝对TEC方法的反演结果基本一致,都与电离层测高仪反演结果符合良好;在较低轨道高度（约500 km）,绝对TEC反演精度优于Abel积分反演精度;（2）绝对TEC反演的最大电子密度N_m较Abel积分法获得的结果更接近于电离层测高获得的峰值电子密度N_mF₂,绝对TEC反演法更加严密和有效;（3）周跳对绝对TEC反演结果的影响较Abel积分反演结果的影响更为敏感,但无论哪种方法,周跳对反演精度都造成严重损失.综合而言,绝对TEC反演法是更优的方法.     

基于等离子体GCPM模型对电离层薄壳模型高度的仿真研究

在基于GPS数据提取电离层总电子含量(TEC)的过程中,电离层薄壳高度的选择对解算电离层垂直TEC的精度有很大的影响.但由于不可能获得一个真实的从电离层D层到GPS卫星高度的电子密度剖面,关于电离层薄壳高度的选择一直是基于GPS数据解算电离层TEC方法中关注的一个问题.本文利用等离子体GCPM模型,对太阳活动高年(2002)和太阳活动低年(2008)情况下电离层有效薄壳高度的选择进行了仿真计算.结果表明,最佳的薄壳高度在2002年为560 km,而在2008年为695 km.通过对全球八个具有代表性地点的仿真计算,揭示了有效薄壳高度更复杂的变化特点.在白天,最佳薄壳的高度变化不大(500 km至750 km);但在夜晚,最佳薄壳高度变化范围很大,甚至可以超过2000 km.此外,本文还对不同卫星仰角的情况下斜向TEC转换为垂直TEC的误差进行了分析,结果表明:随着卫星仰角的增加,薄壳模型带来的转换误差基本上是单调减少的.因而,在实际应用中,尽可能地采用大仰角的卫星数据有助于提高解算的电离层垂直TEC的精度.最后,对全球不同地点的电离层TEC的仿真研究表明,在电子密度水平梯度较大的地区,应用电离层薄壳模型时会导致电子密度较高处的TEC被高估,而电子密度较低处的TEC被低估,在分析基于GPS数据提取的电离层TEC空间变化时要认识到这一点.     

一种基于LEO卫星信标的电离层层析成像新算法

LEO卫星信标是电离层监测的重要手段之一.利用电离层层析成像算法,LEO卫星信标能够实现区域电离层电子密度的快速重构.针对LEO卫星信标的特点,本文提出了一种函数基模型与像素基模型组合的电离层层析成像新算法.选择差分相对电离层总电子含量作为输入数据源,先通过函数基模型法获取电离层电子密度初始分布,再利用像素基模型法对初始分布进行二次迭代重构,该方法可有效降低电离层层析成像对背景电离层模型的依赖,同时能够实现电离层小尺度扰动结构的有效反演.利用数值仿真方法及低纬度电离层层析成像网的实测数据的反演结果验证了本文提出的新算法的可行性和可靠性.     

顶部电离层和等离子体层电子密度分布——基于GRACE星载GPS信标测量的CT反演

本文利用两颗跟飞的GRACE卫星载GPS信标测量数据和基于差分相对TEC的层析算法,实现了全球范围的顶部电离层和等离子体层(450~5000 km) 层析成像.反演结果表明,利用低轨道卫星载GPS信标测量数据可以有效地重建顶部电离层和等离子体层的全球二维分布图像.对不同地磁活动条件下的天基层析反演结果表明,等离子体层电子密度随纬度的分布是不均匀的;在低纬赤道带,从顶部电离层向上延伸直到等离子体层,以及等离子体层中局地的电离增强云团,经常出现近似垂直于磁力线的电子密度柱状增强结构.     

电离层GPS掩星观测改正TEC反演方法

电离层掩星观测中,当低轨卫星（LEO）轨道高度较低时,轨道以上的电离层电子总含量（TEC）对掩星反演的影响不能忽略.目前,一般采用指数函数等外推方法来处理该问题,对反演结果可能引起较大误差.为提高电离层掩星反演精度,本文研究利用LEO处于非掩星一侧GPS观测数据的改正TEC新反演方法.用三维射线追踪程序计算出电离层掩星观测模拟数据,分别应用改正TEC方法和外推方法进行反演,将反演结果与所用模式值进行比较.结果表明：对于轨道高度约800km的GPS/MET掩星模拟数据,外推方法和改正TEC方法反演结果都与模式值基本一致;对于轨道高度约400km的CHAMP掩星模拟数据,外推方法误差较大,改正TEC方法反演结果与模式值相符得较好.将改正TEC方法应用于GPS/MET实测数据的反演,取得了合理的结果.这些说明,改正TEC算法是一种有效的电离层掩星反演方法,尤其是对于轨道较低的LEO的电离层掩星观测反演特别有用.     

GPS地面台网和掩星观测结合的时变三维电离层层析

本文给出GPS地面台网和掩星观测结合的时变三维电离层层析的原理、算法和基于实测数据的反演结果.反演结果的比较表明,联合地基GPS与掩星观测数据进行重建,电子密度整体图像的重建质量特别是其垂直结构的重建质量得到了明显改善.在平静日和磁暴期间两种条件下利用实测数据的重建结果表明,GPS地面台网和掩星观测结合的电离层层析可以获得电离层电子密度在高度-纬度-经度-时间四维空间中的变化.重建结果清晰地显示了磁暴期间电离层负相暴效应,表明结合GPS地面台网和掩星观测的时变三维电离层层析可以有效地监测扰动条件下的大尺度电离层结构.     

Tomographic imaging of the ionosphere using the GPS/MET and NNSS data

The earlier experiments of ionospheric tomography were conducted by receiving satellite signals from ground-based stations and then reconstructing electron density distribution from measures of the total electron content (TEC). In June 1994, National Central University built up the low-latitude ionospheric tomography network (LITN) including six ground stations spanning a range of 16.7° (from 14.6°N to 31.3°N) in latitude within 1° of 121°E longitude to receive the naval navigation satellite system (NNSS) signals (150 and 400 MHz). In the study of tomographic imaging of the ionosphere, TEC data from a network of ground-based stations can provide detailed information on the horizontal structure, but are of restricted utility in sensing vertical structure. However, an occultation observation mission termed the global positioning system/meteorology (GPS/MET) program used a low Earth orbiting (LEO) satellite (the MicroLab-1) to receive multi-channel GPS carrier phase signals (1.5 and 1.2 GHz) and demonstrate active limb sounding of the Earth's atmosphere and ionosphere. In this paper, we have implemented the multiplicative algebraic reconstruction technique (MART) to reconstruct and compare two-dimensional ionospheric structures from measured TECs through the receptions of the GPS signals, the NNSS signals, and/or both of the systems. We have also concluded the profiles retrieved from tomographic reconstruction showing much reasonable electron density results than the original vertical profiles retrieved by the Abel transformation and being in more agreement in peak electron density to nearby ionosonde measurements.     

Ionospheric imaging at mid-latitudes using both GPS and ionosondes

Measurements from the Global Positioning System (GPS) satellites provide a valuable source of information about the ionosphere in the form of ray-path integrations of electron density. Total electron content (TEC) through the ionosphere can be estimated for specific satellite-to-ground paths using the two GPS frequencies and knowledge of the dispersive properties of the ionosphere. One approach used is the ionospheric imaging tool Multi Instrument Data Analysis System (MIDAS), which uses differential phase data from a number of GPS satellites and receivers to create an ionospheric movie of electron density. This paper addresses the accuracy with which MIDAS images the electron density at the F-layer peak. Firstly, the image accuracy is tested using a simulation of the imaging technique, representative of 1 year of data. Experimental GPS phase data are then used to image the electron density during a period of disturbed geomagnetic activity during April 2002. The images are compared to independent measurements from three ionosondes located across Europe and confirm the underestimate in peak electron density that was found in the simulation. Regardless of the peak density errors the vertical TEC in the images remains accurate. The accuracy of the imaged peak electron density is shown to improve across the image when measurements from ionosondes are included in the inversion process.     

Dependence of the ionospheric response on the solar flare parameters based on the theoretical modeling and GPS data

The measurements of an increase in the total electron content (TEC) of the ionosphere during solar flares, obtained based on the GPS data, indicated that up to 30% of TEC increments corresponded to the ionospheric regions above 300 km altitude in some cases, and TEC increased mainly below altitudes of 300 km in other cases. The theoretical model of the ionosphere and plasmasphere was used to study the obtained effects. The altitude-time variations in the charged particle density in the ionospheric region from 100 to 1000 km were used depending on the solar flare spectrum. An analysis of the modeling results indicated that an intensification of the flare UV emission in the 55–65 and 85–95 nm spectral ranges results in a pronounced increase in the electron density in the topside ionosphere (above 300 km). The experimental dependences of the ionospheric TEC response amplitude on the localization and peak power of flares on the Sun in the X-ray range, obtained based on the GPS data, are also presented in the work.     

Reconstructing Regional Ionospheric Electron Density: A Combined Spherical Slepian Function and Empirical Orthogonal Function Approach

The computerized ionospheric tomography is a method for imaging the Earth’s ionosphere using a sounding technique and computing the slant total electron content (STEC) values from data of the global positioning system (GPS). The most common approach for ionospheric tomography is the voxel-based model, in which (1) the ionosphere is divided into voxels, (2) the STEC is then measured along (many) satellite signal paths, and finally (3) an inversion procedure is applied to reconstruct the electron density distribution of the ionosphere. In this study, a computationally efficient approach is introduced, which improves the inversion procedure of step 3. Our proposed method combines the empirical orthogonal function and the spherical Slepian base functions to describe the vertical and horizontal distribution of electron density, respectively. Thus, it can be applied on regional and global case studies. Numerical application is demonstrated using the ground-based GPS data over South America. Our results are validated against ionospheric tomography obtained from the constellation observing system for meteorology, ionosphere, and climate (COSMIC) observations and the global ionosphere map estimated by international centers, as well as by comparison with STEC derived from independent GPS stations. Using the proposed approach, we find that while using 30 GPS measurements in South America, one can achieve comparable accuracy with those from COSMIC data within the reported accuracy (1 × 10¹¹ el/cm³) of the product. Comparisons with real observations of two GPS stations indicate an absolute difference is less than 2 TECU (where 1 total electron content unit, TECU, is 10¹⁶ electrons/m²).     

基于北斗、GLONASS和GPS系统的中低纬电离层特性联合探测

基于GNSS(Global Navigation Satellite Systems)的发展,我们利用具有北斗、GLONASS和GPS三系统信号接收功能的接收机观测的数据,结合电离层总电子含量(Total Electron Content, TEC)的反演算法,提取出GNSS三系统观测的电离层TEC;同时,将GNSS三系统获取的TEC应用到电离层TEC地图、行进式扰动、不规则体结构和电离层的太阳耀斑响应等方面的研究中,这也是首次使用三种GNSS系统数据对电离层进行联合探测研究.研究结果表明,增加了北斗系统的GNSS三系统在研究中国地区电离层TEC地图、周日变化、逐日变化,行进式扰动以及电离层的实时监测等方面较单系统的GPS具有明显的优势.     

Relationships between GPS-signal propagation errors and EISCAT observations

When travelling through the ionosphere the signals of space-based radio navigation systems such as the Global Positioning System (GPS) are subject to modifications in amplitude, phase and polarization. In particular, phase changes due to refraction lead to propagation errors of up to 50 m for single-frequency GPS users. If both the LI and the L2 frequencies transmitted by the GPS satellites are measured, first-order range error contributions of the ionosphere can be determined and removed by difference methods. The ionospheric contribution is proportional to the total electron content (TEC) along the ray path between satellite and receiver. Using about ten European GPS receiving stations of the International GPS Service for Geodynamics (IGS), the TEC over Europe is estimated within the geographic ranges –20° 40°E and 32.5° ø 70°N in longitude and latitude, respectively. The derived TEC maps over Europe contribute to the study of horizontal coupling and transport processes during significant ionospheric events. Due to their comprehensive information about the high-latitude ionosphere, EISCAT observations may help to study the influence of ionospheric phenomena upon propagation errors in GPS navigation systems. Since there are still some accuracy limiting problems to be solved in TEC determination using GPS, data comparison of TEC with vertical electron density profiles derived from EISCAT observations is valuable to enhance the accuracy of propagation-error estimations. This is evident both for absolute TEC calibration as well as for the conversion of ray-path-related observations to vertical TEC. The combination of EISCAT data and GPS-derived TEC data enables a better understanding of large-scale ionospheric processes.