Ionospheric response to the 17 March and 22 June 2015 geomagnetic storms over Wuhan region using GNSS-based tomographic technique

By using the data of GNSS (Global Navigation Satellite System) observation from Crustal Movement Observation Network of China (CMONOC), ionospheric electron density (IED) distributions reconstructed by using computerized ionospheric tomography (CIT) technique are used to investigate the ionospheric storm effects over Wuhan region during 17 March and 22 June 2015 geomagnetic storm periods. F-region critical frequency (foF2) at Wuhan ionosonde station shows an obvious decrease during recovery phase of the St. Patrick’s Day geomagnetic storm. Moreover, tomographic results present that the decrease in electron density begins at 12:00 UT on 17 March during the storm main phase. Also, foF2 shows a long-lasting negative storm effect during the recovery phase of the 22 June 2015 geomagnetic storm. Electron density chromatography presents the evident decrease during the storm day in accordance with the ionosonde observation. These ionospheric negative storm effects are probably associated with changes of chemical composition, PPEF and DDEF from high latitudes.     

Tomographic reconstruction of the ionospheric electron density as a function of space and time

Electron density distribution is the major determining parameter of the ionosphere. Computerized Ionospheric Tomography (CIT) is a method to reconstruct ionospheric electron density image by computing Total Electron Content (TEC) values from the recorded Global Positioning Satellite System (GPS) signals. Due to the multi-scale variability of the ionosphere and inherent biases and errors in the computation of TEC, CIT constitutes an underdetermined ill-posed inverse problem. In this study, a novel Singular Value Decomposition (SVD) based CIT reconstruction technique is proposed for the imaging of electron density in both space (latitude, longitude, altitude) and time. The underlying model is obtained from International Reference Ionosphere (IRI) and the necessary measurements are obtained from earth based and satellite based GPS recordings. Based on the IRI-2007 model, a basis is formed by SVD for the required location and the time of interest. Selecting the first few basis vectors corresponding to the most significant singular values, the 3-D CIT is formulated as a weighted least squares estimation problem of the basis coefficients. By providing significant regularization to the tomographic inversion problem with limited projections, the proposed technique provides robust and reliable 3-D reconstructions of ionospheric electron density.     

Comparison of electron density measurements from CSES and Swarm satellites with GNSS ionospheric tomography data

Electron density measurements obtained from China Seismo‐Electromagnetic Satellite (CSES) and Swarm-B can play an increasingly important role in the study of ionosphere above F2 peak height. This study presented a comprehensive comparison of electron density products obtained from Langmuir probe mounted on CSES and Swarm-B with ionospheric tomography for a whole year period of 2019. CSES was fully compared with Swarm-B on a global scale, including both absolute and relative differences, and a new index called NFI was developed to better quantify the similarity between two latitudinal profiles of electron density. CSES and Swarm-B were then compared with tomography respectively in four regions, roughly located in America, Europe, Australia and China. Results indicated that CSES data are consistent with Swarm-B, as NFI values exceed 0.6 for most of the analyses. Tomography and Swarm-B were found to have a good agreement as their biases are less than 0.2 × 10⁵ el/cm³ in general. For the comparison between CSES and tomography, the bias increased to around 0.6 × 10⁵ el/cm³ but the standard deviation changed slightly, validating the underestimation of electron density by CSES. The spatiotemporal comparisons of CSES and Swarm-B with tomography showed that: 1) the differences in electron density were relatively low in middle latitudes and increased rapidly in the regions of equatorial ionization anomaly; 2) Swarm-B has a better consistent with tomography than CSES, but both are capable of detecting ionosphere anomalies such as midlatitude arcs; and 3) CSES and Swarm-B both can capture the seasonal changes of electron density, while their values are basically smaller than those from tomography in Spring and Summer months.     

Ionospheric path delay models for spaceborne GPS receivers flying in formation with large baselines

GPS relative navigation filters could benefit notably from an accurate modeling of the ionospheric delays, especially over large baselines (>100 km) where double difference delays can be higher than several carrier wavelengths. This paper analyzes the capability of ionospheric path delay models proposed for spaceborne GPS receivers in predicting both zero-difference and double difference ionospheric delays. We specifically refer to relatively simple ionospheric models, which are suitable for real-time filtering schemes. Specifically, two ionospheric delay models are evaluated, one assuming an isotropic electron density and the other considering the effect on the electron density of the Sun aspect angle. The prediction capability of these models is investigated by comparing predicted ionospheric delays with measured ones on real flight data from the Gravity Recovery and Climate Experiment mission, in which two satellites fly separated of more than 200 km. Results demonstrate that both models exhibit a correlation in the excess of 80% between predicted and measured double-difference ionospheric delays. Despite its higher simplicity, the isotropic model performs better than the model including the Sun effect, being able to predict double differenced delays with accuracy smaller than the carrier wavelength in most cases. The model is thus fit for supporting integer ambiguity fixing in real-time filters for relative navigation over large baselines. Concerning zero-difference ionospheric delays, results demonstrate that delays predicted by the isotropic model are highly correlated (around 90%) with those estimated using GPS measurements. However, the difference between predicted and measured delays has a root mean square error in the excess of 30 cm. Thus, the zero-difference ionospheric delays model is not likely to be an alternative to methods exploiting carrier-phase observables for cancelling out the ionosphere contribution in single-frequency absolute navigation filters.     

A new method for improving the performance of an ionospheric model developed by multi-instrument measurements based on artificial neural network

There are remarkable ionospheric discrepancies between space-borne (COSMIC) measurements and ground-based (ionosonde) observations, the discrepancies could decrease the accuracies of the ionospheric model developed by multi-source data seriously. To reduce the discrepancies between two observational systems, the peak frequency (foF2) and peak height (hmF2) derived from the COSMIC and ionosonde data are used to develop the ionospheric models by an artificial neural network (ANN) method, respectively. The averaged root-mean-square errors (RMSEs) of COSPF (COSMIC peak frequency model), COSPH (COSMIC peak height model), IONOPF (Ionosonde peak frequency model) and IONOPH (Ionosonde peak height model) are 0.58 MHz, 19.59 km, 0.92 MHz and 23.40 km, respectively. The results indicate that the discrepancies between these models are dependent on universal time, geographic latitude and seasons. The peak frequencies measured by COSMIC are generally larger than ionosonde’s observations in the nighttime or middle-latitudes with the amplitude of lower than 25%, while the averaged peak height derived from COSMIC is smaller than ionosonde’s data in the polar regions. The differences between ANN-based maps and references show that the discrepancies between two ionospheric detecting techniques are proportional to the intensity of solar radiation. Besides, a new method based on the ANN technique is proposed to reduce the discrepancies for improving ionospheric models developed by multiple measurements, the results indicate that the RMSEs of ANN models optimized by the method are 14–25% lower than the models without the application of the method. Furthermore, the ionospheric model built by the multiple measurements with the application of the method is more powerful in capturing the ionospheric dynamic physics features, such as equatorial ionization, Weddell Sea, mid-latitude summer nighttime and winter anomalies. In conclusion, the new method is significant in improving the accuracy and physical characteristics of an ionospheric model based on multi-source observations.     

磁暴期间中纬度电离层剖面结构变化的数值模拟

利用电离层理论模型模拟了磁暴期间热层大气温度、成分、中性风和电场扰动对电离层电子密度剖面结构，特别是峰值密度和峰值高度变化的影响，结果表明，热层大气温度变化所引起光化反应系数的改变对电离层剖面结构影响不大；热层大气成分特别是Ｎ２／Ｏ的变化能有效地引起密度剖面变化，Ｎ２增加足以使峰值密度产生所观测到的负相暴；由中性风和电场引起等离子体漂移是峰值高度ｈｍＦ２变化的主要原因，但对电子密度的影响不足以抵消     

Proposal for ionospheric tomography technique using a new HF radio source

Measurements of the line integral of the electron density along satellite-to-ground ray paths (i.e. TEC) using differential phase or Doppler of two coherent VHF/UHF signals transmitted from NNSS or GPS satellite networks have been used in the ionospheric tomography for mapping large-scale ionospheric images over region of interest. In this paper, we present theoretical studies of using a new signal source HF frequency in tomographic imaging. The initial phase problem inherent in the phase measurement can be eliminated by measuring Faraday rotations. Relative rotation on two adjacent HF frequencies is used to solve the ambiguity problem. A second-order approximation of the Faraday rotation incorporated with ray-tracing technique improves the reconstruction degradation due to ionospheric refractions. CASSIOPE is a multi-purpose small satellite that receives HF signals from ground radar facilities and it is scheduled for launch in early 2008. Simulations have demonstrated the potential applications of the ionospheric tomography in CASSIOPE/e-POP satellite experiment.     

Symmetry and asymmetry of ionospheric weather at magnetic conjugate points for two midlatitude observatories

Variations of the ionospheric weather W-index for two midlatitude observatories, namely, Grahamstown and Hermanus, and their conjugate counterpart locations in Africa are studied for a period from October 2010 to December 2011. The observatories are located in the longitude sector, which has consistent magnetic equator and geographic equator so that geomagnetic latitudes of the line of force are very close to the corresponding geographic latitudes providing opportunity to ignore the impact of the difference of the gravitational field and the geomagnetic field at the conjugate points on the ionosphere structure and dynamics. The ionosondes of Grahamstown and Hermanus provide data of the critical frequency (foF2), and Global Ionospheric Maps (GIM) provide the total electron content (TECgps) along the magnetic field line up to the conjugate point in the opposite hemisphere. The global model of the ionosphere, International Reference Ionosphere, extended to the plasmasphere altitude of 20,200 km (IRI-Plas) is used to deliver the F2 layer peak parameters from TECgps at the magnetic conjugate area. The evidence is obtained that the electron gas heated by day and cooled by night at the summer hemisphere as compared with the opposite features in the conjugate winter hemisphere testifies on a reversal of plasma fluxes along the magnetic field line by the solar terminator. The ionospheric weather W-index is derived from NmF2 (related with foF2) and TECgps data. It is found that symmetry of W-index behavior in the magnetic conjugate hemispheres is dominant for the equinoxes when plasma movement along the magnetic line of force is imposed on symmetrical background electron density and electron content. Asymmetry of the ionospheric storm effects is observed for solstices when the plasma diffuse down more slowly into the colder winter hemisphere than into the warmer summer hemisphere inducing either plasma increase (positive phase) or decrease (negative phase of W-index) in the ionospheric and plasmaspheric plasma density.     

化学物质释放人工改变电离层

考虑中性气体在电离层高度的扩散过程和相应的电离层离子化学过程，研究了利用主动化学物质释放来改变电离层的方法，理论计算了H2O和SF6两种气体释放后电离层随时间的响应过程，结果表明，在电离层高度上气体的扩散过程非常迅速，电离层F区的电子密度有很大程度的减少，而扩散慢且化学反应快的气体对电离层的影响更大，就更加有利于电离层洞的形成。     

Computerized ionospheric tomography with the IRI model

Computerized ionospheric tomography (CIT) is a method to estimate ionospheric electron density distribution by using the global positioning system (GPS) signals recorded by the GPS receivers. Ionospheric electron density is a function of latitude, longitude, height and time. A general approach in CIT is to represent the ionosphere as a linear combination of basis functions. In this study, the model of the ionosphere is obtained from the IRI in latitude and height only. The goal is to determine the best representing basis function from the set of Squeezed Legendre polynomials, truncated Legendre polynomials, Haar Wavelets and singular value decomposition (SVD). The reconstruction algorithms used in this study can be listed as total least squares (TLS), regularized least squares, algebraic reconstruction technique (ART) and a hybrid algorithm where the reconstruction from the TLS algorithm is used as the initial estimate for the ART. The error performance of the reconstruction algorithms are compared with respect to the electron density generated by the IRI-2001 model. In the investigated scenario, the measurements are obtained from the IRI-2001 as the line integral of the electron density profiles, imitating the total electron content estimated from GPS measurements. It has been observed that the minimum error between the reconstructed and model ionospheres depends on both the reconstruction algorithm and the basis functions where the best results have been obtained for the basis functions from the model itself through SVD.     

Solar activity variations of ionosonde measurements and modeling results

The time series of hourly electron density profiles N(h) obtained at several mid-latitude stations in Europe have been used to obtain N(h) profiles on a monthly basis and to extract both the expected bottomside parameters and a proxy of the ionospheric variability as functions of time and height. With these data we present advances on a “Local Model” technique for the parameters B0 and B1, its applicability to other ionospheric stations, to other bottomside ionospheric parameters, and to modeling the time/height variability of the profile. The Local Model (LM) is an empirical model based on the experimental results of the solar activity dependence of the daily and seasonal behavior of the above parameters. The LM improves the IRI-2001 prediction of the B0 and B1 by factor of two at mid-latitudes. Moreover, the LM can be used to simulate other ionospheric parameters and to build mean N(h) profiles and the deviations from them. The modeling of both the average N(h) profiles and their deviations is an useful tool for ionospheric model users who want to know both the expected patterns and their deviations.     

Deformation of the ionosphere structure during the space weather events on October–November 2003

We present the spatial maps of the ionosphere–plasmasphere slab thickness τ (ratio of the vertical total electron content, TEC, to the F-region peak electron density, NmF2) during the intense ionospheric storms of October–November 2003. The model-assisted technology for estimate of the upper boundary of the ionosphere, hup, from the slab thickness components in the bottomside and topside ionosphere – eliminating the plasmasphere contribution of τ – is applied at latitudes 35° to 70°N and longitudes −10° to 40°E, from the data of 20 observatories of GPS-TEC and ionosonde networks, for selected days and hours of October and November 2003. The daily–hourly values of NmF2, hmF2 and TECgps are used as the constrained parameters for the International Reference Ionosphere extended to the plasmasphere, IRI-Plas, during the ionospheric quiet days, positive and negative storm phases for estimate of τ and hup. Good correlation has been found between the slab thickness and the upper boundary of the ionosphere for the intense ionospheric storms at October–November 2003. During the negative phase of the ionospheric storm, when the ionospheric plasma density is exhausted, the nighttime upper boundary of the ionosphere is greatly uplifted towards the magnetosphere tail, while the daytime upper boundary of the ionosphere is reduced below 500 km over the Earth.     

新乡上空电离层总电子含量的经验模式

基于电离层总电子含量与太阳活动性之间的线性相关关系,利用在新乡长期测量所得的电离层总电子含量的资料,得出了总电子含量的一种经验模式。      

Validation of the STORM model in IRI-2001 at a high latitude station

The International Reference Ionosphere IRI-2001 model contains geomagnetic activity dependence based on an empirical storm time ionospheric correction (STORM model). An extensive validation of the STORM model for the middle latitude region has been performed. In this paper the ability of the STORM model to predict foF2 values at high latitudes is analyzed. For this, ionosonde data obtained at Base Gral. San Martin (68.1°S, 293°E) are compared with those obtained by the IRI-2001 model with or without storm correction during four geomagnetic storms that occurred in 2000 (Rz₁₂ = 117) and 2001 (Rz₁₂ = 111). The results show that predicted values with the STORM model follow the behaviour of foF2 experimental data better than without the STORM model. The relative deviation between measured and predicted foF2 reaches values of up to 24% and 43% with and without the STORM model in IRI-2001, during the main phase of the storms. In order to explain increases of electron density that occurred prior to the storm onset and also decreases of electron density observed during the first part of the recovery of the storm, possible physical mechanisms are discussed.     

Trans-ionospheric GPS signal delay gradients observed over mid-latitude Europe during the geomagnetic storms of October–November 2003

Ionospheric disturbances are known to have adverse effects on the satellite-based communication and navigation. One particular type of ionospheric effects, observed during major geomagnetic storms and threatening the integrity performance of both ground-based and space-based GNSS augmentation systems, is the sharp increase/decrease in the ionospheric delay that propagates in horizontal direction, thus called for convenience ‘moving ionospheric wall’. This paper presents preliminary results from researching such anomalous ionospheric delay gradients at European middle latitudes during the storm events of 29 October 2003 and 20 November 2003. For the purpose, 30-s GPS data from the Belgian permanent network was used for calculating and analysing the slant ionospheric delay and total electron content values. It has been found that, during these two particular storm events, substantial gradients did occur in Europe although they were not so pronounced as in the American sector.     

The impact of spherical symmetry assumption on radio occultation data inversion in the ionosphere: An assessment study

‘Onion-peeling’ is a very common technique used to invert Radio Occultation (RO) data in the ionosphere. Because of the implicit assumption of spherical symmetry for the electron density (N(e)) distribution in the ionosphere, the standard Onion-peeling algorithm could give erroneous concentration values in the retrieved electron density profile. In particular, this happens when strong horizontal ionospheric electron density gradients are present, like for example in the Equatorial Ionization Anomaly (EIA) region during high solar activity periods. In this work, using simulated RO Total Electron Content (TEC) data computed by means of the NeQuick2 ionospheric electron density model and ideal RO geometries, we tried to formulate and evaluate an asymmetry level index for quasi-horizontal TEC observations. The asymmetry index is based on the electron density variation that a signal may experience along its path (satellite to satellite link) in a RO event and is strictly dependent on the occultation geometry (e.g. azimuth of the occultation plane). A very good correlation has been found between the asymmetry index and errors related to the inversion products, in particular those concerning the peak electron density NmF2 estimate and the Vertical TEC (VTEC) evaluation.     

Ionospheric anomalies observed over South Korea preceding the Great Tohoku earthquake of 2011

We investigated the ionospheric anomalies observed before the Tohoku earthquake, which occurred near the northeast coast of Honshu, Japan on 11 March, 2011. Based on data from a ground-based Global Positioning System (GPS) network on the Korean Peninsula, ionospheric anomalies were detected in the total electron content (TEC) during the daytime a few days before earthquake. Ionospheric TEC anomalies appeared on 5, 8 and 11 March. In particular, the ionospheric disturbances on 8 March evidenced a remarkable increase in TEC. The GPS TEC variation associated with the Tohoku earthquake was an increase of approximately 20 total electron content units (TECU), observed simultaneously in local and global TEC measurements. To investigate these pre-earthquake ionospheric anomalies, space weather conditions such as the solar activity index (F10.7) and geomagnetic activity indices (the Kp and Dst indices) were examined. We also created two-dimensional TEC maps to visual the spatial variations in the ionospheric anomalies preceding the earthquake.     

GNSS single frequency ionospheric range delay corrections: NeQuick data ingestion technique

A study of the performance of the NeQuick model and the Klobuchar model for GNSS single frequency range delay correction on a global scale was done using data for moderate solar activity. In this study NeQuick was used in the way intended for Galileo. This study is to assess the performance of the two models at each ionospheric geographic region during moderate solar activity as previously published studies were concentrated only on high solar activity. The results obtained showed that NeQuick outperformed Klobuchar for the whole year at the three geographical regions of the ionosphere. In terms of monthly root mean square of mismodeling, NeQuick outperformed Klobuchar by 15 TECU or more at low-latitudes, 5 TEC or more at mid-latitudes, and 1 TECU or more at high-latitudes.     

NeQuick2 validation using multi-satellite and ground data

Global Navigation Satellite System’s (GNSS) positioning calculation is prone to ionospheric errors. Single frequency GNSS users receive ionospheric corrections through broadcast ionospheric models. Therefore, the accuracy of ionospheric models must be validated based on various geographic and geomagnetic conditions. In this work, an attempt is made to validate NeQuick2 electron density (Ne) using multiple sources of space-based and ground-based data at the Arabian Peninsula and for low solar activity conditions. These sources include space-based data from Swarm, DMSP and COSMIC-2 satellite constellations and ground-based data from GNSS receiver and the ionosonde. The period of this study is 1 year from October 2019 to September 2020. Analysis shows that the agreement between NeQuick2 and experimental Ne close to the peak density height depends on seasons and time of the day with the largest errors found in Autumn and during the daytime. NeQuick2 generally overestimates Ne during the daytime. During the early morning and evening hours, Ne estimates seem to be fairly accurate with slight underestimation in Winter and Spring. Estimation of slab thickness by NeQuick2 is found to be close to the values calculated using collocated ionosonde and GNSS receiver.     

Galileo卫星导航系统广播Nequick模型及其参数拟合

提出了基于IGRF模型的Galileo广播Nequick模型及其参数拟合算法, 解决了Galileo信号仿真中地理场景映射与地磁坐标下的电离层延时修正参数拟合问题. 应用IGRF模型, 可计算出任意给定位置和时间点的地磁参数以及E层、 F₁层、F₂层的电子密度, 从而计算出Galileo电离层修正参数. 仿真结果表明, 该算法拟合的全球电离层延时与IGS提供的实际观测值基本一致, 仿真精度高于一般的经验电离层模型, 实现了Galileo卫星信号的电离层延时修正参数的精确仿真.