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

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

Automated daily process for global ionospheric total electron content maps and satellite ocean altimeter ionospheric calibration based on Global Positioning System data

The accuracy of single-frequency ocean altimeters benefits from calibration of the total electron content (TEC) of the ionosphere below the satellite. Data from a global network of Global Positioning System (GPS) receivers provides timely, continuous, and globally well-distributed measurements of ionospheric electron content. For several months we have been running a daily automatic Global Ionospheric Map process which inputs global GPS data and climatological ionosphere data into a Kalman filter, and produces global ionospheric TEC maps and ocean altimeter calibration data within 24 h of the end-of-day. Other groups have successfully applied this output to altimeter data from the GFO satellite and in orbit determination for the TOPEX/Poseidon satellite. Daily comparison of the global TEC maps with independent TEC data from the TOPEX altimeter is performed as a check on the calibration whenever the TOPEX data are available. Comparisons of the global TEC maps against TOPEX data will be discussed. Accuracy is best at mid-to-high absolute latitudes (∣latitude∣>30°) due to the better geographic distribution of GPS receivers and the relative simplicity of the ionosphere. Our highly data-driven technique is relatively less accurate at low latitudes and especially during ionospheric storm periods, due to the relative scarcity of GPS receivers and the structure and volatility of the ionosphere. However, it is still significantly more accurate than climatological models.     

Prophet时序预测模型在电离层TEC异常探测中的应用

提出一种基于Facebook 开源的Prophet预测模型进行电离层TEC异常识别的新方法。 首先, 对比分析了该方法与传统时间序列预测方法(ARIMA模型等)预测电离层TEC建模背景值的精度, 以及与经典电离层TEC异常识别方法(滑动四分位法)提取前面对应一致的电离层TEC背景值的精度。 结果表明, Prophet预测模型预测建模背景值的精度要明显优于其他方法, 且预测的建模精度比ARIMA模型等方法高2.55倍左右, 比滑动四分位法高10.74倍左右。 同时, 在最佳预测建模区间时, 其精度值大小比较依次为RMSE_IQR=10.5841>RMSE_ARIMA=3.2780>RMSE_Prophet=0.8469, 说明传统探测法预测建模背景值时具有较大的不足。 随后, 以2017年8月8日九寨沟7.0级地震为例, 利用该方法分析了电离层TEC异常扰动情况, 并对比验证了该方法的有效性和准确性。 实验结果表明: 在震前第10 d和第2 d电离层TEC发生较为明显的负异常, 第7 d电离层TEC发生较为明显的正异常。 对比实验表明, Prophet预测模型的有效性和准确性明显优于滑动四分位法。     

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.     

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

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

汶川大地震前后电离层电子含量的扰动分析

为了分析地震前后电离层电子含量的变化,本文选取四川汶川发生8.0级地震前后10天电离层中的电子含量作为研究对象,使用IGS网站提供的VTEC数据,采取统计学分析方法,结合太阳和地磁活动信息,对此次地震前后电离层VTEC资料进行了详细的统计分析和讨论,数据分析结果表明大地震前电离层电子含量会出现明显扰动,电离层中电子含量的异常变化有可能是地震发生的前兆信息.     

中国电离层TEC现报系统

作为最重要的电离层参量之一,电离层电子浓度总含量(TEC)可以通过当前广泛使用的全球定位系统(GPS)的信标进行观测.我们在我国北起漠河、经北京和武汉、南到三亚四个观测站建立了GPS接收站,经单站数据处理后将原始的单站GPS TEC观测数据上载到北京数据处理中心;采用我们发展的经验基函数模式算法,用实测数据估算格点TEC并提供给用户,同时生成覆盖中国疆域的TEC地图并在因特网上实时发布.这一电离层TEC现报系统是我国首个类似的技术系统,在观测站布局和TEC地图算法上有所创新.该系统已用于实时监测我国电离层环境,并可为我国卫星定位导航和测控等技术系统的电波修正提供实测电离层数据.     

Prompt derivation of TEC from GEONET data for space weather monitoring in Japan

Continuous monitoring of ionospheric conditions is essential to monitoring and forecasting space weather. The worldwide use of global navigation satellite systems like the Gobal Positioning System (GPS) makes it possible to continuously monitor the total electron content (TEC) of the ionosphere and plasmasphere up to a height of about 20,000 km. We have developed a system for deriving the TEC from GEONET data rapidly and we use the TEC distribution over Japan in the daily operations of the Space Weather Forecast Center at NICT (RWC Tokyo of ISES). Using instrumental biases from a few days before enables us to drastically shorten the processing time for deriving TEC. The latest TEC values (with a delay of about 1 h) are obtained every 3 h, and most of the values are within 2 TEC units of the actual TEC. We have found our system for deriving TEC rapidly to be useful for continuously monitoring the progress of ionospheric storms under any ionospheric conditions, even those under which the usual ionosonde observations are unable to obtain F-region profiles.     

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

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

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.     

Large-scale disturbances of auroral origin during strong magnetic storms of October 29–31, 2003, and November 7–11, 2004, according to the data of the GPS network and ionosondes

The intensity of large-scale traveling ionospheric disturbances (LS TIDs), registered according to measurements of the total electron content (TEC) during the magnetic storms of October 29–31, 2003, and November 7–11, 2004, has been compared with that of local electron density disturbances. The data of TEC measurements at ground-based GPS receivers located near the ionospheric stations and the corresponding values of the critical frequency of the ionospheric F region (foF2) were used for this purpose. The variations in TEC and foF2 were similar for all events mentioned above. The previous assumption that the region of thickness 150–200 km in the vicinity of the ionospheric F region mainly contributes to TEC modulation was confirmed for the cases when the electron density disturbance at an F region maximum was not more than 50%. However, this region probably becomes more extensive in vertical when the electron density disturbance in the vicinity of the ionospheric F region is about 85%.     

基于COSMIC资料分析电离层F层不规则体结构

根据电离层不规则体的产生会导致周围电子浓度发生起伏变化的原理,利用2007年COSMIC掩星系统的TEC数据,通过平滑滤波得到TEC的扰动值ΔTEC的变化,利用其研究F层不规则体的时空变化特征.统计结果表明:扰动较大的掩星事件主要发生在磁纬±20°之间和高纬地区,春季和秋季带状分布较为明显,不同经度地区较强扰动的掩星事件的分布也有不同特征;较强ΔTEC的掩星事件主要发生在地方时午夜前和午夜后两个时段,发生的高度主要在250~400km范围内.这些结果与已知的F层不规则体的时空分布特征较为一致,说明利用TEC的扰动量来分析电离层F层不规则体结构是可行的.     

基于GNSS数据的朝鲜核爆电离层效应观测研究

本文分析了电离层对2009年、2013年和2016年朝鲜核试验中地下核爆(Underground Nuclear Explosion,UNE)的响应.利用南北半球IGS站的GNSS-TEC观测数据,发现了在3次核试验期间的磁共轭电离层扰动现象.观测结果表明,3次UNE所产生的电离层扰动分别以228m·s^-1、173m·s^-1和147m·s^-1的速度从核试验爆心地区径向传播.本文研究提出,UNE所产生的电离层TEC扰动是岩石圈-大气层-电离层耦合(Lithosphere-Atmosphere-Ionosphere Coupling,LAIC)电场渗透到电离层高度引发电动力学过程的结果.此外,LAIC电场可以沿磁力线映射到共轭半球,从而在共轭地区造成电离层TEC扰动.因此,在核试验地区附近以及其对应的磁共轭地区,UNE所产生的LAIC电场在电离层TEC扰动的形成中起着关键性的作用.     

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

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

利用陆态网络实时监测电离层TEC及在地震监测中的应用分析

中国大陆构造环境监测网络(陆态网络)是中国最大的GPS综合服务网络,该网络包括260个分布全国的GPS基准站。本文首先介绍了利用GPS研究电离层TEC的基本原理与方法,然后阐述了陆态网络电离层TEC监测系统的数据处理及产品发布的详细流程,最后分析了震前电离层异常现象与地震发生的耦合关系,并利用陆态网络电离层TEC监测系统的产品分析了四川省芦山7.0级地震的同震电离层扰动现象。     

Prediction of global positioning system total electron content using Neural Networks over South Africa

Global positioning system (GPS) networks have provided an opportunity to study the dynamics and continuous changes in the ionosphere by supplementing ionospheric studies carried out using various techniques including ionosondes, incoherent scatter radars and satellites. Total electron content (TEC) is one of the physical quantities that can be derived from GPS data, and provides an indication of ionospheric variability. This paper presents a feasibility study for the development of a Neural Network (NN) based model for the prediction of South African GPS derived TEC. Three South African locations were identified and used in the development of an input space and NN architecture for the model. The input space included the day number (seasonal variation), hour (diurnal variation), Sunspot Number (measure of the solar activity), and magnetic index (measure of the magnetic activity). An analysis was done by comparing predicted NN TEC with TEC values from the IRI-2001 version of the International Reference Ionosphere (IRI), validating GPS TEC with ionosonde TEC (ITEC) and assessing the performance of the NN model during equinoxes and solstices. For this feasibility model, GPS TEC was derived for a limited number of years using an algorithm still in the early phases of validation. However, results show that NNs predict GPS TEC more accurately than the IRI at South African GPS locations, but that more good quality GPS data is required before a truly representative empirical GPS TEC model can be released.     

中国电离层TEC同化现报系统

数据同化是在基于物理机制的背景模型上,融合时空不规则分布的观测数据的一种现报方法.同化能够有效弥补数据的时空局限和模型的精度偏差,使二者相互匹配从而获得更加合理可信的模拟效果.本研究利用电离层数据同化方法,针对中国及周边区域（15°N-55°N,70°E-140°E）构建了电离层总电子含量（TEC）同化现报系统.系统使用国际参考电离层（IRI）作为背景场,利用中国科学院空间环境监测网和国际GNSS服务组织（IGS）的部分地基GNSS台站数据作为观测值,并采用三维变分与Gauss-Markov卡尔曼滤波相结合的算法进行背景场和观测值的数据同化,生成覆盖中国及周边区域的电离层TEC和GPS单频接收机延迟误差的格点化准实时现报地图,并在中国科学院空间环境预报中心（http：//sepc.ac.cn/TEC_chn.php）网上发布,每15 min进行更新.该系统是我国基于同化算法的电离层现报系统之一,已用于中国及周边区域的电离层环境实时监测,可为卫星导航、雷达成像、短波通信等科学研究和工程应用提供相对及时、准确、有效的电离层TEC和误差修正信息.     

电离层垂直TEC映射函数的实验观测与统计特性

利用GPS信标测量获得的电离层电子浓度总含量（TEC）是沿电波路径的斜向TEC.理论研究和实际应用中,常常需要通过映射函数将斜向TEC转换为垂直方向的TEC,这在当前主要采用对电子浓度分布模型的数值积分得到模型映射函数来实现.本文在考察现有不同模型映射函数的基础上,又提出了一种源于实际观测的实验映射函数的概念与估算方法.我们利用IGS的全球GPS观测站的斜向TEC和JPL提供的垂直TEC数据获得了2006年期间的实验映射函数,并对所得结果进行了初步统计分析.在卫星天顶角较小时,上述实验映射函数和模型映射函数之间相差甚微,均可很好描述垂直TEC与斜TEC之间关系;但卫星天顶角较大时,实验映射函数和常用的模型映射函数之间存在明显差异.本文认为,这种差异主要是因为现有模型映射函数中没有考虑到等离子体层的贡献.我们认为采用基于实验映射函数的模式,或者通过考虑等离子体层的贡献对现有模型映射函数进行改进,可以有效提高电离层TEC的估算精度.     

Dynamics of high-latitude patches and associated small-scale irregularities during the October and November 2003 storms

Observations from a network of specially equipped GPS scintillation receivers in Northern Europe are used to investigate the dynamics of ionospheric plasma during the storm events of 30 October and 20 November 2003. The total electron content (TEC) and scintillation data, combined with ionospheric tomography produced by the multi-instrument data analysis system (MIDAS), reveal strong enhancements and steep gradients in TEC during nighttime under a prevailing negative B_z component of the interplanetary magnetic field (IMF). Amplitude and phase scintillation maxima are often co-located with the TEC gradients at the edge of plasma patches, revealing the presence of small-scale irregularities and suggesting association with a tongue of ionization (TOI) convecting in an anti-sunward direction from the American sector across the polar cap. Similarities and differences between the ionospheric response to the two storms are investigated. The 30 October event reveals a quite complex scenario showing two phases of plasma dynamics: the former reflects the expected convection pattern for IMF B_z southward and the latter possibly indicates a sort of TEC plasma stagnation signature of the more complex convection patterns during several positive/negative excursions of IMF B_z.     

COSMIC低轨卫星GPS接收机差分码偏差估计

GPS接收机差分码偏差（Differential Code Bias,DCB）是利用COSMIC低轨卫星观测值反演电离层总电子含量TEC的一项重要误差源.本文将COSMIC卫星轨道高度以上的电离层作为一个单层,采用球谐函数来参数化电离层TEC值,并利用最小二乘法同时估算电离层球谐系数和DCB参数.运用这种方法对2012年12月份的所有COSMIC卫星GPS接收机DCB进行了解算,并与COSMIC数据分析与档案中心CDAAC提供的产品进行了比较.实验结果表明：在2012年12月期间,估计的接收机DCB与CDAAC结果符合的较好,二者DCB变化趋势相近,DCB差值的RMS值在2 TECU以内,且最大绝对差值小于3 TECU;此外,本文计算的接收机DCB估计误差主要分布在0.2~0.4 TECU之间,具有较高的内符合精度.