利用数据同化技术实现区域电离层TEC重构

电离层参量的提取是开展电离层研究的基础,而数据同化技术则是获取电离层参量的一种重要手段。以NeQuick模型的输出作为背景场,Kalman滤波作为同化算法,利用数据同化技术实现区域电离层TEC重构,结果表明,数据同化方法重构的倾斜总电子含量（TEC）和垂直TEC与实测值较为一致。相比NeQuick模型及全球电离层地图（GIM）数据,数据同化方法重构得到的TEC的平均误差和标准差均有明显的降低,实测数据验证了数据同化技术在区域TEC重构中的精度和可靠性。     

电离层TEC的预测模型

电离层总电子含量（TEC）的精确预报对提高GNSS导航精度,保障无线电空间远程通讯具有重要作用。分析了IGS发布的电离层格网点总电子含量（TEC）的时间序列特点,基于时间序列分析理论,以AR模型对格网点TEC随机时间序列平稳化后建模和预报。实例分析表明,研究的预报技术和方法是可行的。     

Eye on the Ionosphere: The Correlation between Solar 10.7 cm Radio Flux and Ionospheric Range Delay

Patricia Doherty joins the regular contributors of this column to discuss the correlation between measurements of solar 10.7 cm radio flux and ionospheric range delay effects on GPS. Mrs. Doherty has extensive experience in the analysis of ionospheric range delays from worldwide systems and in the utilization and development of analytical and theoretical models of the Earth's ionosphere. Ionospheric range delay effects on GPS and other satellite ranging systems are directly proportional to the Total Electron Content (TEC) encountered along slant paths from a satellite to a ground location. TEC is a highly variable and complex parameer that is a function of geographic location, local time, season, geomagnetic activity, and solar activity. When insufficiently accounted for, ionospheric TEC can seriously limit the performance of satellite ranging applications. Since the ionosphere is a dispersive medium, dual-frequency Global Positoning System (GPS) users can make automatic corrections for ionospheric range delay by computing the apparent difference in the time delays between the two signals. Single-frequency GPS users must depend on alternate methods to account for the ionospheric range delay. Various models of the ionosphere have been used to provide estimates of ionospheric range delay. These models range from the GPS system's simple eight-coefficient algorithm designed to correct for approximately 50% rms of the TEC, to state-of-the-art models derived from physical first principles, which can correct for up to 70 to 80% rms of the TEC but at a much greater computational cost. In an effort to improve corrections for the day-to-day variability of the ionosphere, some attempts have been made to predict the TEC by using the daily values of solar 10.7 cm radio flux (F_10,7). The purpose of this article is to show that this type of prediction is not useful due to irregular, and sometimes very poor, correlation between daily values of TEC and F_10.7. Long-term measurements of solar radio flux, however, have been shown to be well correlated with monthly mean TEC, as well as with the critical frequency of the inonospheric F2 region (foF2), which is proportional to the electron density at the peak of the ionospheric F2 region. ? 2000 John Wiley & Sons, Inc.     

基于IGS的南北半球TEC非对称性研究

利用2000~2012年的IGS电离层TEC数据,研究了南北半球电离层TEC的非对称性。首先对南北半球TEC日平均值的不对称指数(AI)进行了研究,结果表明,AI的极大值的绝对值集中在0.4左右,极小值的绝对值集中在0.3左右,表明南北半球电离层TEC存在数值上的非对称性;通过傅里叶变换,发现AI存在1 a和1/3 a周期;每年北半球电离层TEC高于南半球电离层TEC的时间比例均大于50%,平均值为58%,表明了南北半球电离层TEC存在时间上的非对称性。然后进一步分析了南北半球TEC年平均值的差异,发现了北半球TEC年平均值高于南半球的规律,2000~2012年的南北半球TEC年平均值的不对称指数在0.036左右。最后从全球大气分布及运动规律的南北非对称性和热层年度变化赤道非对称性两个角度,对南北半球电离层TEC的非对称性产生的原因做了初步分析。     

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

介绍了震前电离层TEC异常探测原理的研究进展和主要的异常探测方法,详细介绍了时间序列法、Kalman滤波和小波变换在电离层异常探测中的原理和应用,总结了传统方法和新方法在电离层异常探测的应用发展情况,并分析了每种方法的优势与不足。为未来中国采用电离层异常探测开展地震预测工作提出了建议。     

太阳耀斑的GPS监测方法及实例分析

利用GPS伪距与载波相位联合数据处理的方法，分析了2000年7月14日太阳耀斑爆发期间，武汉、北京、乌鲁木齐GPS观测数据得到的电离层TEC，提出了利用多项式拟合计算由耀斑引起的电离层TEC增加量的方法。     

单站区域电离层TEC建模及精度分析

为了分析单站区域电离层总电子含量(total electron content,TEC)模型的适用范围和精度,基于2～15阶次球谐函数,分别建立了欧洲区域16个单站区域电离层TEC模型,生成了区域格网TEC,并与欧洲定轨中心(Center for Orbit Determination in Europe,CODE)、...     

Evaluation of NTCM-BC and a proposed modification for single-frequency positioning

Ionospheric delay is a dominant factor that affects the accuracy of single-frequency positioning. Thus, an empirical ionospheric model with high accuracy is very important for single-frequency users. This study proposes a modified empirical broadcast ionospheric model, called MNTCM-BC, based on the Neustrelitz Total Electron Content (TEC) broadcast model NTCM-BC. Nine daily ionospheric coefficients of these models are estimated using datasets of the previous day from 30 globally distributed Global Navigation Satellite System monitor stations, and the prediction performance of the MNTCM-BC is evaluated with the datasets of the current day from all 300 verification stations. The results show that the complex behavior of the ionosphere is well described by the MNTCM-BC, including the visibility of two ionization crests on both sides of the geomagnetic equator and the TEC variations that depend on the local time and geomagnetic latitude. In terms of the prediction accuracy, compared with the NTCM-BC, the main improvement in the MNTCM-BC is achieved in summer, whereas the accuracy is comparable in other seasons. Hence, the following analyses are focused on summer. In the low-solar activity year of 2009, the prediction accuracy of the MNTCM-BC is improved by 0.11 TECU compared with that of the NTCM-BC. As to the high-solar activity year of 2014, the corresponding improvement is 0.35 TECU. In addition, when the number of monitor stations is increased from 30 to 300, the prediction accuracy of two models can be slightly improved by 0.06 TECU in 2009 and 0.13 TECU in 2014, respectively, while reliability enhances. Furthermore, the average three-dimensional positioning accuracy of 160 globally distributed stations for single-frequency point positioning using the Klobuchar model, the NTCM-BC and the MNTCM-BC is 1.83, 1.21 and 1.20 m during quiet day and 3.15, 2.31 and 2.21 m during perturbed day, respectively. Relative to the Klobuchar model and the NTCM-BC, the average accuracy improvements in the MNTCM-BC are about 30 and 3%, respectively.     

Spherical cap harmonic model for mapping and predicting regional TEC

An approach to modeling the regional ionospheric total electron content (TEC) based on spherical cap harmonic analysis is presented. This approach not only provides a better regional TEC mapping accuracy, but also the capability for ionospheric model prediction based on spectrum analysis and least squares collocation. Unlike conventional approaches, which predict the immediate TEC with models using current observations, the spherical cap harmonic approach utilizes models using past observations to predict a model which will provide future TEC values. A significant advantage in comparison with conventional approaches is that the spherical cap harmonic approach can be used to predict the long-term TEC with reasonable accuracy. This study processes a set of GPS data with an observation time span of 1 year from two GPS networks in China. The TEC mapping accuracy of the spherical cap harmonic model is compared with the polynomial model and the global ionosphere model from IGS. The results show that the spherical cap harmonic model has a better TEC mapping accuracy with smoother residual distributions in both temporal and spatial domains. The TEC prediction with the spherical cap harmonic model has been investigated for both short- and long-term intervals. For the short-term interval, the prediction accuracies for the latencies of 1-day, 2-days, and 3-days are 2.5 TECU, 3.5 TECU, and 4.5 TECU, respectively. For the long-term interval, the prediction accuracy is 4.5 TECU for a 2-month latency.     

Ionospheric TEC predictions over a local area GPS reference network

Single layer ionosphere models are frequently used for ionospheric modeling and estimation using GPS measurements from a network of GPS reference stations. However, the accuracies of single layer models are inherently constrained by the assumption that the ionospheric electrons are concentrated in a thin shell located at an altitude of about 350 km above Earths surface. This assumption is only an approximation to the physical truth because the electrons are distributed in the entire ionosphere region approximately from 50 to 1,000 km. To provide instantaneous ionospheric corrections for the real-time GPS positioning applications, the ionospheric corrections need to be predicted in advance to eliminate the latency caused by the correction computation. This paper will investigate ionospheric total electron content (TEC) predictions using a multiple-layer tomographic method for ionospheric modeling over a local area GPS reference network. The data analysis focuses on the accuracy evaluation of short-term (5 min in this study) TEC predictions. The results have indicated that the obtainable TEC prediction accuracy is at a level of about 2.8 TECU in the zenith direction and 95% of the total electron content can be recovered using the proposed tomography-based ionosphere model.     

不同TEC分布模式对星载SAR立体定位精度的影响及校正

基于SAR成像模型和R-D定位模型,分析了电离层对星载SAR立体定位精度的影响,提出了相应的校正方法,并进行了仿真实验。实验表明对于特定波长的SAR信号,不同TEC分布模式组合对SAR立体像对定位精度的影响不同,校正效果也不同:当两幅影像对应的TEC分布均匀且大小相同时,电离层对SAR立体定位精度的影响小到可以忽略;当两幅影像由于拍摄季节不同导致TEC不同时,电离层影响很大,特别是距离向精度,其影响程度与两幅影像对应的TEC差值正相关,此时定向参数精化法校正效果很好;当电离层发生扰动时,电离层影响随着扰动程度的增大而增大,定向参数精化法有一定效果,但当电离层扰动很大时,校正效果明显减弱,还需要通过斜距直接校正法等其他方法予以校正。     

Application of the TaiWan Ionospheric Model to single-frequency ionospheric delay corrections for GPS positioning

The performance of a three-dimensional ionospheric electron density model derived from FormoSat3/COSMIC GPS Radio Occultation measurements, called the TaiWan Ionosphere Model (TWIM), in removing the ionospheric delays in single-frequency pseudorange observations is presented. Positioning results using TWIM have been compared with positioning results using other ionospheric models, such as the Klobuchar (KLOB) and the global ionospheric model (GIM). C/A code pseudoranges have been observed at three International GPS Service reference stations that are representative of mid-latitude (BOR1 and IRKJ) and low-latitude (TWTF) regions of the ionosphere. The observations took place during 27 geomagnetically quiet days from April 2010 to October 2011. We perform separate solutions using the TWIM, KLOB, GIM ionospheric models and carry out a solution applying no ionospheric correction at all. We compute the daily mean horizontal errors (DMEAN) and the daily RMS (DRMS) for these solutions with respect to the published reference station coordinates. It has demonstrated that TEC maps generate using the TWIM exhibit a detailed structure of the ionosphere, particularly at low-latitude region, whereas the Klobuchar and the GIM only provide the basic diurnal and geographic features of the ionosphere. Also, it is shown that even for lower satellite elevations, the TWIM provides better positioning than the Klobuchar and GIM models. Specifically, using TWIM, the difference of the uncorrected solution (no ionospheric correction), and the other solutions, relative to the uncorrected solution, is 45 % for the mean horizontal error (DMEAN) and 42 % for the horizontal root-mean-square error (DRMS). Using Klobuchar and GIM, the percent for DMEAN only reaches to about 12 % and 3 %, while the values for the DRMS are only 12 and 4 %, respectively. In the vertical direction, all models have a percentage of about 99 and 70 % for the mean vertical error (VMEAN) and vertical root-mean-square error (VRMS), respectively. These percentages show the greater impact of TWIM on the ionospheric correction compared to the other models. In at least 40 % of the observed days and across all stations, TWIM has the smallest DMEAN, VMEAN, DRMS, and VRMS daily values. These values reach 100 % at station TWTF. This shows the overall performance of TWIM is better than the Klobuchar and GIM.     

自回归移动平均模型的电离层总电子含量短期预报

摘 要：本文在充分考虑乘积性季节模型的情况下,利用差分法对电离层总电子含量（Total Electron Content,TEC）样本序列进行平稳化处理后,采用时间序列分析中的求和自回归移动平均模型（简称ARIMA,Autoregressive Integrated Moving Average)对TEC值序列进行预报分析。以欧洲定轨道中心（CODE）提供的2008-2012年电离层TEC值为样本数据,分析了该方法在电离层平静期、活跃期预报高、中、低不同纬度电离层TEC值的精度以及TEC样本数据的长短对预报精度的影响等。实验结果表明：在电离层平静期和活跃期预报6天的平均相对精度可达83.3%和86.6%;而平均预报残差分别为0.18±1.9TECU和0.69±2.6TECU,其中预报残差小于3TECU分别达到90%和81%以上;而且两个时期都具有纬度越高相对精度越低而绝对精度越高的规律。此外,预报精度会随TEC样本序列长度增加而提高,但40天左右为其最佳样本长度,如超过此长度,其精度会逐渐降低;而相同样本数据的预报精度会随预报长度的增加而减小,初期并不明显,但超过30天其相对精度将随时间明显降低。     

The activities of the Ionosphere Working Group of the International GPS Service (IGS)

This article is based on a position paper presented at the IGS Network, Data and Analysis Center Workshop 2002 in Ottawa, Canada, 8–11 April 2002, and introduces the IGS Ionosphere Working Group (Iono_WG). Detailed information about the IGS in general can be found on the IGS Central Bureau Web page: http://igscb.jpl.nasa.gov. The Iono_WG commenced working in June 1998. The working group's main activity currently is the routine production of ionosphere Total Electron Content (TEC) maps with a 2-h time resolution and daily sets of GPS satellite and receiver hardware differential code bias (DCB) values. The TEC maps and DCB sets are derived from GPS dual-frequency tracking data recorded with the global IGS tracking network. In the medium- and long-term, the working group intends to refine algorithms for the mapping of ionospheric parameters from GPS measurements and to realize near–real–time availability of IGS ionosphere products. The paper will give an overview of the Iono_WG activities that include a summary of activities since its establishment, achievements and future plans. Electronic Publication     

局部均值分解和奇异值分解在GNSS站坐标时间序列信号降噪中的应用

为了有效地提取GNSS站坐标时间序列的有用信息，降低噪声干扰，本文提出一种局部均值分解和奇异值分解相结合的信号降噪方法，并利用5个测站的实测坐标时间序列对新方法进行了验证。首先通过局部均值分解将坐标时间序列分解成一系列PF分量和余项，然后利用连续均方误差方法确定高频分量与低频分量的分界点，保持低频分量不变，运用奇异值分解方法对高频分量进行降噪重构，最后将重构的高频分量与低频分量叠加得到最终的降噪坐标时间序列，并对降噪效果进行对比分析。结果表明，与单纯的奇异值分解方法相比，局部均值分解和奇异值分解相结合方法能够自适应地选择合适的奇异值个数进行信号重构，提高了降噪效果。     

利用BP神经网络改进电离层短期预报模型

电离层总电子含量TEC(Total Electron Content)是电离层的一个重要特征参数。对TEC的预报也已经成为电离层研究的一个热点。根据JS CORS中心提供的GPS观测数据,建立了区域实时多站多项式模型;并分别以模型计算得到的南京地区的电离层电子含量数据和苏州地区的电离层电子含量数据为样本,采用时间序列和BP神经网络融合模型进行了预报。结果表明,采用融合模型在短期预报中能够取得较好的效果,精度比时间序列模型提高20%左右。     

时间序列分析在电离层短期预报中的应用

根据IGS提供的2012年TEC数据,在得到TEC值残差序列的基础上,利用谱分析去掉周期项和Matlab工具箱去掉趋势项并采用平滑算法去噪后,对随机项进行时间序列分析。目前,电离层格网点的预报均建立在原始TEC上,阐述了采用差分后的TEC值进行预报的方法,以减少日变化周期项的影响。根据AR（p）模型预报的结果加上周期项和趋势项后,再加上前一天对应时段的TEC值与IGS发布的数据比较,结果表明,该方法利用短期IGS发布的TEC值进行电离层预报能取得较高的精度。     

电离层对星载SAR矩形包络线性调频信号的影响

讨论了电离层对星载SAR矩形包络线性调频信号的影响 ,推导出由电离层引起的图像偏移和峰值二次相位误差的计算公式 ,基于这些公式对P ,L ,C 3个波段的影响作了比较分析。结果表明 :电离层对P波段的影响很大 ,对L波段有一定影响 ,对C波段的影响基本可以忽略。     

Variability of Altimetric Range Correction Parameters over Indian Tropical Region using JASON-1 & JASON-2 Radar Altimeters

Region-specific atmospheric range correction maps are generated over the Indian tropical region from Jason-1 & Jason-2 radar altimeters data. Seasonal and spatial variability of wet tropospheric correction (WTC), ionospheric correction (IC), dry tropospheric correction (DTC), and sea state bias (SSB) correction are analyzed over the Bay of Bengal and the Arabian Sea. Two year atmospheric range correction data from JASON-1 (2008) & JASON-2 altimeters (2009) has been analyzed where each Jason cycle is exactly 9.9156?days repeat. The monthly and yearly mean variation of the range correction parameters has been studied over the Indian continent. For precise study, four different regions were selected as the Region of Interest in the North & South of the Arabian Sea and Bay of Bengal. WTC, Significant Wave Height (SWH), Wind Speed (WS) and SSB show the higher values during monsoon months. The yearly mean WTC over Indian Tropical region was 26.22?cm in 2008 and 26.20?cm in 2009. SSB Correction values mainly depend on the SWH and wind speed. The yearly mean SSB correction over Indian Tropical region was 6.87?cm in 2008 and 7.02?cm in 2009. DTC values are less during monsoon season and it shows a high value in the month of January. The yearly mean DTC over Indian Tropical region was 230.42?cm in 2008 and 230.43?cm in 2009.The IC values mainly depend on frequency and total electron content (TEC) in the ionosphere which further depends on the solar activity. The yearly mean IC over Indian Tropical region was higher in 2008 (2.98?cm) in comparison to mean IC in 2009 (2.29?cm). This study is useful to understand the variability of atmospheric correction parameters especially over Indian continent.     

Global Ionosphere Maps of VTEC from GNSS,satellite altimetry,and formosat-3/COSMIC data

For space geodetic techniques, operating in microwave band, ionosphere is a dispersive medium; thus signals traveling through this medium are in the first approximation affected proportional to inverse of the square of their frequencies. This effect allows gaining information about the parameters of the ionosphere in terms of Total Electron Content (TEC) or the electron density (N _e ). TEC or electron density can then be expressed by means of spherical harmonic base functions to provide a Global Ionosphere Map (GIM). The classical input data for development of GIMs are obtained from dual-frequency observations carried out at Global Navigation Satellite Systems (GNSS) stations. However, GNSS stations are in-homogeneously distributed around the world, with large gaps particularly over the oceans; this fact reduces the precision of the GIM over these areas. On the other hand, dual-frequency satellite altimetry missions such as Jason-1 provide information about the ionosphere precisely above the oceans; and furthermore Low Earth Orbiting (LEO) satellites, such as Formosat-3/COSMIC (F/C) provide well-distributed information of ionosphere globally. This study investigates on global modeling of TEC through combining GNSS and satellite altimetry data with global TEC data derived from the occultation measurements of the F/C mission. The combined GIMs of vertical TEC (VTEC) show a maximum difference of 1.3–1.7 TEC units (TECU) with respect to the GNSS-only GIMs in the whole day. The root mean square error (RMS) maps of combined solution show a reduction of about 0.1 TECU in the whole day. This decrease of RMS can reach up to 0.5 TECU in areas where no or few GNSS observations are available, but high number of F/C measurement is carried out. This proves that the combined GIMs provide a more homogeneous global coverage and higher reliability than results of each single method. All comparisons and validations made within this study provide vital information regarding combination and integration of various observation techniques in the Global Geodetic Observing System of the International Association of Geodesy.