Near real-time modeling of global ionospheric vertical total electron content using hourly IGS data

The International GNSS Service (IGS) has been providing reliable Global Ionospheric Maps (GIMs) since 1998. The Ionosphere Associate Analysis Centers (IAACs) model the global ionospheric Total Electron Content (TEC) and generate the daily GIM products within the context of the IGS. However, the rapid and final daily GIM products have a latency of at least one day and one week or so, respectively. This limits the value of GIM products in real-time GNSS applications. We propose and develop an approach for near real-time modeling of global ionospheric TEC by using the hourly IGS data. We perform an experiment in a real operating environment to generate near real-time GIM (named BUHG) products for more than two years. Final daily GIM products, Precise Point Positioning (PPP) based VTEC resources, and JASON-3 Vertical TEC (VTEC) measurements are collected for testing the performance of BUHG. The results show that the performance of BUHG is very close to that of the daily GIM products. Also, there is good agreement between BUHG and PPP-derived VTEC as well as with JASON-3 VTEC. It is possible that BUHG would be further improved with an increase in available hourly GNSS data.     

基于陆态网的区域电离层TEC空间变化特性分析

为研究中国陆态网区域电离层TEC在空间小尺度、高分辨率情况下的变化特性及适用精度范围,利用陆态网260个GNSS连续运行观测站数据,解算并生成2016-2017年731天陆态网区域电离层RIM格网,并进行精度验证.在同一RIM格网中,分别在经度和纬度方向上对间隔不同经纬度的TEC格网点作差分析.结果表明:陆态网区域内经度方向上TEC最大变化率和平均变化率分别为0.30TECU·（°）-1和0.11TECU·（°）-1;经度间隔1°时,TEC差值小于2TECU,且随着经度间隔的增大,其TEC差值也随之增大,并表现出一定的半年和周年变化规律;纬度方向上TEC最大变化率和平均变化率分别为1.7TECU·（°）-1和0.46TECU·（°）-1;陆态网区域内电离层TEC随纬度减小而增大,纬度间隔1°时,99.4%的TEC差值小于4TECU,且随着纬度间隔的增大,其TEC差值也随之增大,并表现出一定的半年和周年变化规律;间隔相同情况下,纬度方向上TEC的变化比经度方向上大.      

Ionospheric response to solar and magnetospheric protons during January 15–22, 2005: EAGLE whole atmosphere model results

We present an analysis of the ionosphere and thermosphere response to Solar Proton Events (SPE) and magnetospheric proton precipitation in January 2005, which was carried out using the model of the entire atmosphere EAGLE. The ionization rates for the considered period were acquired from the AIMOS (Atmospheric Ionization Module Osnabrück) dataset. For numerical experiments, we applied only the proton-induced ionization rates of that period, while all the other model input parameters, including the electron precipitations, corresponded to the quiet conditions. In January 2005, two major solar proton events with different energy spectra and proton fluxes occurred on January 17 and January 20. Since two geomagnetic storms and several sub-storms took place during the considered period, not only solar protons but also less energetic magnetospheric protons contributed to the calculated ionization rates. Despite the relative transparency of the thermosphere for high-energy protons, an ionospheric response to the SPE and proton precipitation from the magnetotail was obtained in numerical experiments. In the ionospheric E layer, the maximum increase in the electron concentration is localized at high latitudes, and at heights of the ionospheric F2 layer, the positive perturbations were formed in the near-equatorial region. An analysis of the model-derived results showed that changes in the ionospheric F2 layer were caused by a change in the neutral composition of the thermosphere. We found that in the recovery phase after both solar proton events and the enhancement of magnetospheric proton precipitations associated with geomagnetic disturbances, the TEC and electron density in the F region and in topside ionosphere/plasmasphere increase at low- and mid-latitudes due to an enhancement of atomic oxygen concentration. Our results demonstrate an important role of magnetospheric protons in the formation of negative F-region ionospheric storms. According to our results, the topside ionosphere/plasmasphere and bottom-side ionosphere can react to solar and magnetospheric protons both with the same sign of disturbances or in different way. The same statement is true for TEC and foF2 disturbances. Different disturbances of foF2 and TEC at high and low latitudes can be explained by topside electron temperature disturbances.     

The post sunset equatorial F- region zonal drift variability and its linkage with equatorial spread F onset and duration over Indian longitudes

The seasonal and solar activity variation of the post sunset F- region zonal plasma drift, at the magnetic equatorial region over Indian longitudes is analyzed using the Republic of China Satellite-1 data from January 2000 to April 2004. The post sunset F- region zonal drifts are observed to be higher in the years of high solar activity in comparison with low solar activity, while seasonally the drifts are minimum in summer with much higher values in other seasons. The seasonal and solar activity variations of zonal plasma drift are attributed to the corresponding variations in the neutral winds. The dependences of the F region peak vertical drift on the zonal plasma drift at 18.5 IST (Indian Standard Time) and the time difference of the conjugate points sunset times, are quantitatively analyzed. Further an integrated parameter (incorporating the above mentioned two independent factors), which is able to predict the peak vertical drift and growth rate of Rayleigh Taylor instability is proposed. The other major outcome of the study is the successful prediction of the Equatorial Spread F (ESF) onset time and duration using the new integrated parameter at 18.5 IST. ESF irregularities and associated scintillations adversely affect communication and navigation systems. Hence, the present methodology for the prediction of the characteristics of these nocturnal irregularities becomes relevant.     

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

A prediction model of short-term ionospheric foF2 based on AdaBoost

In this paper, the AdaBoost-BP algorithm is used to construct a new model to predict the critical frequency of the ionospheric F2-layer (foF2) one hour ahead. Different indices were used to characterize ionospheric diurnal and seasonal variations and their dependence on solar and geomagnetic activity. These indices, together with the current observed foF2 value, were input into the prediction model and the foF2 value at one hour ahead was output. We analyzed twenty-two years’ foF2 data from nine ionosonde stations in the East-Asian sector in this work. The first eleven years’ data were used as a training dataset and the second eleven years’ data were used as a testing dataset. The results show that the performance of AdaBoost-BP is better than those of BP Neural Network (BPNN), Support Vector Regression (SVR) and the IRI model. For example, the AdaBoost-BP prediction absolute error of foF2 at Irkutsk station (a middle latitude station) is 0.32 MHz, which is better than 0.34 MHz from BPNN, 0.35 MHz from SVR and also significantly outperforms the IRI model whose absolute error is 0.64 MHz. Meanwhile, AdaBoost-BP prediction absolute error at Taipei station from the low latitude is 0.78 MHz, which is better than 0.81 MHz from BPNN, 0.81 MHz from SVR and 1.37 MHz from the IRI model. Finally, the variety characteristics of the AdaBoost-BP prediction error along with seasonal variation, solar activity and latitude variation were also discussed in the paper.     

Detection of the ionospheric disturbances on GPS-TEC using Differential Rate Of TEC (DROT) algorithm

The solar, geomagnetic, gravitational and seismic activities can cause spatial and temporal (hourly, diurnal, seasonal and annual) variabilities of the ionosphere. Main observable ionospheric parameters such as Total Electron Content (TEC) can be used to quantify these. TEC is the total number of electrons on a ray path crossing the atmosphere. The network of world-wide Global Positioning System (GPS) receivers provide a cost-effective solution in estimating TEC over a significant proportion of global land mass. This study is focused on the analysis of the variations of ionosphere over a midlatitude region using GPS-TEC estimates for three Sun Spot Numbers (SSN) periods. The investigation is based on a fast and automatic variability detection algorithm, Differential Rate Of TEC (DROT). The algorithm is tested using literature data on disturbances generated by a geomagnetic activity, a Solar Flare, a Medium Scale Travelling Ionospheric Disturbance (MSTID), a Large Scale TID (LSTID) and an earthquake. Very good agreement with the results in the literature is found. DROT is applied to IONOLAB-TEC estimates from nine Turkish National Permanent GPS Network (TNPGN Active) stations over Turkey to detect the any wave-like oscillations, sudden disturbances and other irregularities during December, March, June, September months for 2010, 2011, 2012 years. It is observed that DROT algorithm is capable of detecting both small and large scale variability due to climatic, gravitational, geomagnetic and solar activities in all layers of ionosphere. The highest DROT values are observed in 2010 during winter months. In higher solar activity years of 2011 and 2012, DROT is able to indicate both seasonal variability and severe changes in ionosphere due to increased number of geomagnetic storms and local seismic activities.     

On the estimation of regional covariance functions of TEC variations over Canada

Spatial and temporal variations of Total Electron Content (TEC) can affect GNSS high accuracy positioning. Enhanced estimation of ionospheric variations and their de-correlation can benefit differential and point positioning rapid solutions. Global and regional TEC maps can provide the overall state of ionopsheric variations in space and time domains within their accuracy limits. In this paper, these maps are exploited to retrieve ionospheric variations by means of variograms and their associated covariance functions of TEC residuals over Canadian region during geomagnetically quiet and disturbed conditions. A number of theoretical variogram functions are reviewed for modeling covariance of TEC residuals. The variogram modeling of residuals during a strong geomagnetic storm revealed variances of one order of magnitude larger compared to a rather quiet condition. Variogram models are also used in regional and local kriging interpolation experiments and their performances are evaluated. Global maps of TEC RMS by International GNSS Service and two of its analysis centres are also compared over the Canadian region during a two-year period. Realistic representation of regional variances using estimated variograms when compared to global ionospheric RMS maps are also presented.     

Low latitude ionospheric TEC responses to dynamical complexity quantifiers during transient events over Nigeria

In this study, the values of chaoticity and dynamical complexity parameters for some selected storm periods in the year 2011 and 2012 have been computed. This was done using detrended TEC data sets measured from Birnin-Kebbi, Torro and Enugu global positioning system (GPS) receiver stations in Nigeria. It was observed that the significance of difference (SD) values were mostly greater than 1.96 but surprisingly lower than 1.96 in September 29, 2011. The values of the computed SD were also found to be reduced in most cases just after the geomagnetic storm with immediate recovery a day after the main phase of the storm while the values of Lyapunov exponent and Tsallis entropy remains reduced due to the influence of geomagnetic storms. It was also observed that the value of Lyapunov exponent and Tsallis entropy reveals similar variation pattern during storm period in most cases. Also recorded surprisingly were lower values of these dynamical quantifiers during the solar flare event of August 8th and 9th of the year 2011. The possible mechanisms responsible for these observations were further discussed in this work. However, our observations show that the ionospheric effects of some other possible transient events other than geomagnetic storms can also be revealed by the variation of chaoticity and dynamical complexity.     

基于Hatch滤波法的GPS伪距差分定位分析

利用Hatch滤波法对导航数据进行伪距差分定位研究,对单频接收机采用Hatch滤波法时电离层的影响做进一步研究,分析电离层风暴对定位精度的影响。为了提高定位精度,采用Hatch滤波法辅助的差分定位技术,对解算过程进行详细分析,并通过MATLAB软件对数据进行处理仿真,验证利用差分Hatch滤波法提高定位精度的可行性。