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.     

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.     

WAAS系统中电离层折射校正的新方法及计算结果

电离层介质的色散性是影响电磁波信号进行卫星导航定位精度的重要因素之一.配合北斗二代分系绩研制任务,提出了一种新的电离层折射校正算法,并利用2000年7月1日到3日的双频GPS观测数据对6个用户站进行试算,进一步将试算所得均方根误差和电离层网格算法得到的误差进行比较.结果表明,对于中纬区域的用户站,估算的TEC误差约为0.5 m左右;而低纬用户误差相对增大,为1 m左右.文中给出的算法与电离层网格模型所提供的精度相差不多,在未来中国自主的卫星增强系统中采用新方法进行电离层进行修正是可行的及有效的.      

Model based Computerized Ionospheric Tomography in space and time

Reconstruction of the ionospheric electron density distribution in space and time not only provide basis for better understanding the physical nature of the ionosphere, but also provide improvements in various applications including HF communication. Recently developed IONOLAB-CIT technique provides physically admissible 3D model of the ionosphere by using both Slant Total Electron Content (STEC) measurements obtained from a GPS satellite - receiver network and IRI-Plas model. IONOLAB-CIT technique optimizes IRI-Plas model parameters in the region of interest such that the synthetic STEC computations obtained from the IRI-Plas model are in accordance with the actual STEC measurements. In this work, the IONOLAB-CIT technique is extended to provide reconstructions both in space and time. This extension exploits the temporal continuity of the ionosphere to provide more reliable reconstructions with a reduced computational load. The proposed 4D-IONOLAB-CIT technique is validated on real measurement data obtained from TNPGN-Active GPS receiver network in Turkey.     

2004年11月一次磁暴期间全球电离层TEC扰动分析

利用全球分布的GPS原始观测数据提取的电离层总电子含量(TEC)分析了2004年11月6日至12日期间全球电离层暴的形态特点与发展过程.结果表明,11月8日磁暴主相期间电离层暴以大范围的强烈正暴为主,在11月10日的恢复相,Dst又一次降到最低值前后期间,电离层再次受到很强的扰动,大范围的正暴和负暴交替出现.这次磁暴期间夏季半球的负暴更加强烈,反映出负暴偏向于在夏季半球发生的季节变化特点.另外,磁暴期间,夜晚TEC值普遍比磁暴前的平静期要低,具体是什么机制导致还需要进一步收集数据和分析.      

地磁暴期间中国中低纬电离层的CIT研究

利用电离层层析成像技术(Computerized Ionospheric Tomography, CIT)处理115°E子午圈附近6个台站的GPS观测数据, 分析了2004年11月地磁暴期间中国中低纬电离层的响应情况. 结果表明, 电离层呈正相扰动, 且不同高度上的响应不同, 800 km以下电子密度有不同程度的增加, 且在峰值高度附近增幅最大, 800 km以上地磁暴的影响并不显著; 伴随地磁能量的注入, 赤道异常峰极向扩展; 随磁扰强度的降低, 电子密度也逐渐恢复至平静水平. 这些结果与以往的理论和观测结果一致, 初步估计扰动是由热层暴环流引起的, 并受到赤道异常峰移动的影响.      

A real-time ionospheric model based on GNSS Precise Point Positioning

This paper proposes a method of real-time monitoring and modeling the ionospheric Total Electron Content (TEC) by Precise Point Positioning (PPP). Firstly, the ionospheric TEC and receiver’s Differential Code Biases (DCB) are estimated with the undifferenced raw observation in real-time, then the ionospheric TEC model is established based on the Single Layer Model (SLM) assumption and the recovered ionospheric TEC. In this study, phase observations with high precision are directly used instead of phase smoothed code observations. In addition, the DCB estimation is separated from the establishment of the ionospheric model which will limit the impacts of the SLM assumption impacts. The ionospheric model is established at every epoch for real time application. The method is validated with three different GNSS networks on a local, regional, and global basis. The results show that the method is feasible and effective, the real-time ionosphere and DCB results are very consistent with the IGS final products, with a bias of 1–2 TECU and 0.4 ns respectively.     

Ionosphere tomography based on spaceborne SAR

Two models of ionosphere tomography based on spaceborne SAR (Synthetic Aperture Radar) are proposed. For HF-SAR the signal with sweeping frequency lower than the characteristic frequency of ionosphere will be scatted during the ionosphere propagation and completely reflected at a corresponding height. The ionospheric electron density isolines looked as series of random surfaces can be reconstructed from the HF-SAR echoes by using the inverse scattering technique for layered rough surfaces and the method of moment (MoM). The numerical simulation show that due to the MoM can provide a full wave solution, the ionosphere tomography with high resolution can be obtained as long as enough sampling data of HF-SAR echoes are used. For VHF/UHF/P/L-band SAR the TEC (Total Electron Content) can be obtained from the SAR echoes scattered by some strong point targets (such as the calibrators, etc.) appeared in the SAR imaged ground region, and the ionosphere tomography can be performed by CT technique.     

Multi-station observation of ionospheric irregularities over South Africa during strong geomagnetic storms

This paper presents results pertaining to the response of the mid-latitude ionosphere to strong geomagnetic storms that occurred from 31 March to 02 April 2001 and 07–09 September 2002. The results are based on (i) Global Positioning Systems (GPSs) derived total electron content (TEC) variations accompanying the storm, (ii) ionosonde measurements of the ionospheric electrodynamic response towards the storms and (iii) effect of storm induced travelling ionospheric disturbances (TIDs) on GPS derived TEC. Ionospheric data comprising of ionospheric TEC obtained from GPS measurements, ionograms, solar wind data obtained from Advanced Composition Explorer (ACE) and magnetic data from ground based magnetometers were used in this study. Storm induced features in vertical TEC (VTEC) have been obtained and compared with the mean VTEC of quiet days. The response of the mid-latitude ionosphere during the two storm periods examined may be characterised in terms of increased or decreased level of VTEC, wave-like structures in VTEC perturbation and sudden enhancement in hmF2 and h′F. The study reveals both positive and negative ionospheric storm effects on the ionosphere over South Africa during the two strong storm conditions. These ionospheric features have been mainly attributed to the travelling ionospheric disturbances (TIDs) as the driving mechanism for the irregularities causing the perturbations observed. TEC perturbations due to the irregularities encountered by the satellites were observed on satellites with pseudo random numbers (PRNs) 15, 17, 18 and 23 between 17:00 and 23:00 UT on 07 September 2002.     

Evaluation of the IRI-2016 and NeQuick electron content specification by COSMIC GPS radio occultation,ground-based GPS and Jason-2 joint altimeter/GPS observations

We examined performance of two empirical profile-based ionospheric models, namely IRI-2016 and NeQuick-2, in electron content (EC) and total electron content (TEC) representation for different seasons and levels of solar activity. We derived and analyzed EC estimates in several representative altitudinal intervals for the ionosphere and the plasmasphere from the COSMIC GPS radio occultation, ground-based GPS and Jason-2 joint altimeter/GPS observations. It allows us to estimate a quantitative impact of the ionospheric electron density profiles formulation in several altitudinal intervals and to examine the source of the model-data discrepancies of the EC specification from the bottom-side ionosphere towards the GPS orbit altitudes. The most pronounced model-data differences were found at the low latitude region as related to the equatorial ionization anomaly appearance. Both the IRI-2016 and NeQuick-2 models tend to overestimate the daytime ionospheric EC and TEC at low latitudes during all seasons of low solar activity. On the contrary, during high solar activity the model results underestimated the EC/TEC observations at low latitudes. We found that both models underestimated the EC for the topside ionosphere and plasmasphere regions for all levels of solar activity. For low solar activity, the underestimated EC from the topside ionosphere and plasmasphere can compensate the overestimation of the ionospheric EC and, consequently, can slightly decrease the resulted model overestimation of the ground-based TEC. For high solar activity, the underestimated EC from the topside ionosphere and plasmasphere leads to a strengthening of the model underestimation of the ground-based TEC values. We demonstrated that the major source of the model-data discrepancies in the EC/TEC domain comes from the topside ionosphere/plasmasphere system.     

Improvement of the real time total electron content based on the International Reference Ionosphere model

One of the methods to adapt the International Reference Ionosphere model to real time conditions is to use instantaneous values of the critical frequency of the ionosphere foF2. It is shown that there can be large discrepancies between model and experimental values of the total electron content TEC after this adaptation. Improvement can be provided by the use of an empirical model of the ionospheric slab thickness τ. This conclusion is based on analysis of contributions of various correction factors connected with foF2 and τ into discrepancies between model and experimental values of TEC.     

2000年7月和2003年10月大磁暴期间东亚地区中低纬电离层的GPS TEC的响应研究

利用武汉电离层观象台研制的GPS TEC的现报方法及现报系统,对东亚地区GPS台网的观测数据进行处理分析,特别对2000年7月14-18日和2003年10月28日至11月1日两次特大磁暴期间的数据进行了对比考察,文中分析了两次磁暴间的电离层响应,得到对应不同磁暴时段电离层TEC的不同变化情况,着重揭示了TEC赤道异常峰的压缩和移动以及赤道异常随时间的压缩—反弹—恢复的过程,并结合高纬电离层的部分响应机制进行了说明,结果显示,两次磁暴期的电离层响应表现出了各自不同的特点,从而反映出因季节变化引起的高纬电离层暴时能量注入的不同而造成的全球性电离层扰动的不同形态,由此看出,磁暴期间电离层TEC的变化直接与太阳扰动发生的时间及其对高纬电离层的耦合有关,若短时期内连续发生多次磁暴,则电离层反应更加复杂,不能简单地当做单一磁暴叠加处理。     

Regional 4-D modeling of the ionospheric electron density from satellite data and IRI

Accurate knowledge of the electron density is the key point in correcting ionospheric delays of electromagnetic measurements and in studying ionosphere physics. During the last decade Global Navigation Satellite Systems (GNSS) have become a promising tool for monitoring ionospheric parameters such as the total electron content (TEC). In this contribution we present a four-dimensional (4-D) model of the electron density consisting of a given reference part, i.e., the International Reference Ionosphere (IRI), and an unknown correction term expanded in terms of multi-dimensional base functions. The corresponding series coefficients are calculable from the satellite measurements by applying parameter estimation procedures. Since satellite data are usually sampled between GPS satellites and ground stations, finer structures of the electron density are modelable just in regions with a sufficient number of ground stations. The proposed method is applied to simulated geometry-free GPS phase measurements. The procedure can be used, for example, to study the equatorial anomaly.     

The impact of large solar events on the total electron content of the ionosphere at mid latitudes

Ionospheric disturbances associated with solar activity may occur via two basic mechanisms. The first is related to the direct impact on the ionosphere of EUV photons from a flare, and the second by prompt electric field penetration into the magnetosphere during geomagnetic storms. In this paper we examine the possibility that these two mechanisms may have an impact at mid latitudes by calculating the total electron content (TEC) from GPS stations in Mexico during several large X-ray flares. We have found that indeed large, complex flares, which are well located, may affect the mid latitude ionosphere. In fact, in the solar events of July 14, 2000 and April 2001 storms, ionospheric disturbances were observed to increase up to 138 and 150 TECu, respectively, due to the influence of EUV photons. Also, during the solar events of July 2000, April 2001, Halloween 2003, January 2005 and December 2006, there are large ionospheric disturbances (up to 393 TECu in the Halloween Storms), due to prompt penetration electric field, associated with CME producing geomagnetic storm.     

Response of the ionosphere to super geomagnetic storms: Observations and modeling

The relative importance of the main drivers of positive ionospheric storms at low-mid latitudes is studied using observations and modeling for the first time. In response to a rare super double geomagnetic storm during 07–11 November 2004, the low-mid latitude (17°–48°N geomag. lat.) ionosphere produced positive ionospheric storms in peak electron density (NmF2) in Japan longitudes (≈125°–145°E) on the day of main phase (MP1) onset (06:30 LT) and negative ionospheric storms in American longitudes (≈65°–120°W) on the following day of MP1 onset (13:00–16:00 LT). The relative effects of the main drivers of the positive ionospheric storms (penetrating daytime eastward electric field, and direct and indirect effects of equatorward neutral wind) are studied using the Sheffield University Plasmasphere Ionosphere Model (SUPIM). The model results show that the penetrating daytime (morning–noon) eastward electric field shifts the equatorial ionisation anomaly crests in NmF2 and TEC (total electron content) to higher than normal latitudes and reduces their values at latitudes at and within the anomaly crests while the direct effects of the equatorward wind (that reduce poleward plasma flow and raise the ionosphere to high altitudes of reduced chemical loss) combined with daytime production of ionisation increase NmF2 and TEC at latitudes poleward of the equatorial region; the later effects can be major causes of positive ionospheric storms at mid latitudes. The downwelling (indirect) effect of the wind increases NmF2 and TEC at low latitudes while its upwelling (indirect) effect reduces NmF2 and TEC at mid latitudes. The net effect of all main drivers is positive ionospheric storms at low-mid latitudes in Japan longitude, which qualitatively agrees with the observations.     

Ionospheric effects of the solar flares as deduced from global GPS network data

Results derived from analysing the ionospheric response to faint and bright solar flares are presented. The analysis used novel technology of a global detection of ionospheric effects from solar flares as developed by the authors (Afraimovich, 2000a; Afraimovich, 2000b), on the basis of phase measurements of the total electron content (TEC) in the ionosphere using an international GPS network. The essence of the method is that use is made of appropriate filtering and a coherent processing of variations in the TEC which is determined from GPS data, simultaneously for the entire set of visible GPS satellites at all stations used in the analysis. This technique is useful for identifying the ionospheric response to faint solar flares (of X-ray class C) when the variation amplitude of the TEC response to separate line-on-sight (LOS) is comparable to the level of background fluctuations. The dependence of the TEC variation response amplitude on the flare location on the Sun is investigated.     

Near real-time PPP-based monitoring of the ionosphere using dual-frequency GPS/BDS/Galileo data

Ionosphere delay is very important to GNSS observations, since it is one of the main error sources which have to be mitigated even eliminated in order to determine reliable and precise positions. The ionosphere is a dispersive medium to radio signal, so the value of the group delay or phase advance of GNSS radio signal depends on the signal frequency. Ground-based GNSS stations have been used for ionosphere monitoring and modeling for a long time. In this paper we will introduce a novel approach suitable for single-receiver operation based on the precise point positioning (PPP) technique. One of the main characteristic is that only carrier-phase observations are used to avoid particular effects of pseudorange observations. The technique consists of introducing ionosphere ambiguity parameters obtained from PPP filter into the geometry-free combination of observations to estimate ionospheric delays. Observational data from stations that are capable of tracking the GPS/BDS/GALILEO from the International GNSS Service (IGS) Multi-GNSS Experiments (MGEX) network are processed. For the purpose of performance validation, ionospheric delays series derived from the novel approach are compared with the global ionospheric map (GIM) from Ionospheric Associate Analysis Centers (IAACs). The results are encouraging and offer potential solutions to the near real-time ionosphere monitoring.     

Low-latitude ionospheric scintillations and total electron content obtained with the CITRIS instrument on STPSat1 using radio transmissions from DORIS ground beacons

The primary objective of the Scintillation and Tomography Receiver in Space (CITRIS) is to detect ionospheric irregularities from space at low latitude. For this purpose, the satellite receiver uses the UHF and S-Band transmissions of the ground network of Doppler Orbitography and Radiopositioning Integrated by Satellite (DORIS) beacons. CITRIS, developed at the Naval Research Laboratory, differs from the normal DORIS receiver by being able to capture and store the complex amplitude of the 401.25 and 2036.25 MHz transmissions at 200 Hz sample rate. Ground processing of the CITRIS data yields total electron content (TEC) and both phase and amplitude scintillations. With CITRIS flying on the US Space Test Program (STP) satellite STPSat1, 2 years of data were collected and processed to determine the fluctuations in ionospheric TEC and radio scintillations associated with equatorial irregularities. CITRIS flights over DORIS transmitters yield direct measurements of the horizontal plasma density fluctuations associated with equatorial plasma bubbles. Future flights of CITRIS can provide valuable complements to other satellite instruments such as GPS occultation receivers used to estimate vertical electron density profiles in the ionosphere.     

电离层TEC日落阶段涨落异常增强现象研究

重点报道了一种TEC信距测量的夜间抖动异常现象，分析了这种现象的观测特征，指出了夜间GPS伪距观测是的异常抖动，特别是日落后时段GPS伪距观测量的散开十分严重的乐是个别的偶然现象，在考察了周围环境之后初步排除了这种现象是由周围地物引起的可能性，而进一步与电离层不均匀体（Spread-F)的特征，特别是它发生率最大的地方时进行比较，认为这是一种电离层效应，由大尺度的电离层不规则结构造成的散射很可能是这种夜间散开的主要原因。本文仅就单站GPS观测资料认证了这种夜间强烈抖动的存在，分析了北京上空电离层中发生的不均匀体或各种扰动对GPS信号的可能影响及其观测特征，指出了夜间GPS伪距观测量的异常抖动，特别是日落后时段GPS伪距观测量的散开十分严重的现象，可以用于电离层不规则结构的研究。由于GPS设备比较简便，数据量大，精度高，适用于各种环境等特点，我们发展的单站数据分析方法可能会对研究电离层TEC夜间涨落异常增加现象，包括地域，时间上的统计特性比较有利。