Total electron content measurements in ionospheric physics

With the advent of modern global networks of dual-frequency Global Positioning System (GPS), total electron content (TEC) measurements along slant paths connecting GPS receivers and satellites at 22,000 km have become the largest data set available to ionospheric scientists. The TEC can be calculated from the time and phase delay in the GPS signal using the GPS Toolkit, but an unknown bias will remain. In addition, UHF/VHF radio beacons on board low-Earth-orbiting satellites can also be used to measure the electron content. However, the TEC measurements are obtained by integrating TEC differences between slant paths, but also contain biases. It is often necessary to use data assimilative algorithms like the Ionospheric Data Assimilation Three-Dimensional (IDA3D), and to treat both GPS- and LEO-beacon TEC measurements as relative data in order to conduct ionospheric studies.     

基于双频GPS观测的电离层TEC与硬件延迟反演方法

利用全球定位系统（Global Positioning System,GPS）的双频观测数据反演得到电离层的总电子含量（Total Electron Content,TEC）,使得广域甚至全球范围高时空分辨率的电离层观测研究成为可能,但由于GPS卫星和接收机对信号的硬件延迟可导致TEC测量系统偏差,因此,需要探索反演TEC并估测GPS卫星与接收机硬件延迟的有效算法.本文根据电离层电波传播理论,阐述了基于双频GPS观测提取电离层TEC的方法,给出TEC与硬件延迟的基本关系.综合研究了TEC与硬件延迟的反演方法,进行分析与归纳分类,在此基础上提出了有待深入研究的问题.      

月球电离层掩星探测研究

日本SELENE/KAGUYA探测任务提供了研究月球电离层的机会。采用无线电掩星探测技术和趋势外推算法,消除地球电离层和行星际等离子体的干扰影响,残余的信号相位信息的变化反映了月球电离层的信息,估算出月球周围附近近似对称分布的稀薄电离层中电子总含量(TEC)约为每立方米10-14个。     

Demonstration of the use of the Doppler Orbitography and Radio positioning Integrated by Satellite (DORIS) measurements to validate GPS ionospheric imaging

There is a lack of independent ionospheric data that can be used to validate GPS imaging results at mid latitudes over severe storm times. Doppler Orbitography and Radio positioning Integrated by Satellite (DORIS), a global network of dual-frequency ground to satellite observations, provides this missing data and here is employed as verification to show the accuracy of the ionospheric GPS images in terms of the total electron content (TEC). In this paper, the large-scale ionospheric structures that appeared during the strong geomagnetic storm of 20 November 2003 are reconstructed with a GPS tomographic algorithm, known as MIDAS, and validated with DORIS TEC measurements. The main trough shown in an extreme equatorward position in the ionospheric imaging over mainland Europe is confirmed by DORIS satellite measurements. Throughout the disturbed day, the variations of relative slant TECs between DORIS data and MIDAS results agree quite well, with the average of the mean differences about 2 TECu. We conclude that as a valuable supplement to GPS data, DORIS ionospheric measurements can be used to analyse TEC variations with a relatively high resolution, ∼10 s in time and tens of kilometres in space. This will be very helpful for identification of some highly dynamic structures in the ionosphere found at mid-latitudes, such as the main trough, TID (Travelling Ionospheric Disturbances) and SED (Storm Enhanced Density), and could be used as a valuable auxiliary data source in ionospheric imaging.     

Regional 4-D modeling of the ionospheric electron density

The knowledge of the electron density is the key point in correcting ionospheric delays of electromagnetic measurements and in studying the ionosphere. During the last decade GNSS, in particular GPS, has become a promising tool for monitoring the total electron content (TEC), i.e., the integral of the electron density along the ray-path between the transmitting satellite and the receiver. Hence, geometry-free GNSS measurements provide informations on the electron density, which is basically a four-dimensional function depending on spatial position and time. In addition, these GNSS measurements can be combined with other available data including nadir, over-ocean TEC observations from dual-frequency radar altimetry (T/P, JASON, ENVISAT), and TECs from GPS-LEO occultation systems (e.g., FORMOSAT-3/COSMIC, CHAMP) with heterogeneous sampling and accuracy.     

Ionospheric scintillations at Guilin detected by GPS ground-based and radio occultation observations

The occurrence of ionospheric scintillations with S4 ? 0.2 was studied using GPS measurements at Guilin, China (25.29°N, 110.33°E; geomagnetic: 15.04°N, 181.98°E), a station located near the northern crest of the equatorial anomaly. The results are presented for data collected from January 2009 to March 2010. The results show that nighttime amplitude scintillations only took place in February and March of the considered years, while daytime amplitude scintillations occurred in August and December of 2009. Nighttime amplitude scintillations, observed in the south of Guilin, always occurred with phase scintillations, TEC (Total Electron Content) depletions, and ROT (Rate Of change of TEC) fluctuations. However, TEC depletions and ROT fluctuations were weak during daytime amplitude scintillations, and daytime amplitude scintillations always took place simultaneously for most of the GPS satellites which appeared over Guilin in different azimuth directions. Ground-based GPS scintillation/TEC observations recorded at Guilin and signal-to-noise-ratio (SNR) measurements obtained from GPS-COSMIC radio occultation indicate that nighttime and daytime scintillations are very likely caused by ionospheric F region irregularities and sporadic E, respectively. Moreover, strong daytime amplitude scintillations may be associated with the plasma density enhancements in ionospheric E region caused by the Perseid and Geminid meteor shower activities.     

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.     

Plasmaspheric electron content derived from GPS TEC and FORMOSAT-3/COSMIC measurements: Solar minimum condition

The plasmaspheric electron content (PEC) was estimated by comparison of GPS TEC observations and FORMOSAT-3/COSMIC radio occultation measurements at the extended solar minimum of cycle 23/24. Results are retrieved for different seasons (equinoxes and solstices) of the year 2009. COSMIC-derived electron density profiles were integrated up to the height of 700 km in order to retrieve estimates of ionospheric electron content (IEC). Global maps of monthly median values of COSMIC IEC were constructed by use of spherical harmonics expansion. The comparison between two independent measurements was performed by analysis of the global distribution of electron content estimates, as well as by selection specific points corresponded to mid-latitudes of Northern America, Europe, Asia and the Southern Hemisphere. The analysis found that both kinds of observations show rather similar diurnal behavior during all seasons, certainly with GPS TEC estimates larger than corresponded COSMIC IEC values. It was shown that during daytime both GPS TEC and COSMIC IEC values were generally lower at winter than in summer solstice practically over all specific points. The estimates of PEC (h > 700 km) were obtained as a difference between GPS TEC and COSMIC IEC values. Results of comparative study revealed that for mid-latitudinal points PEC estimates varied weakly with the time of a day and reached the value of several TECU for the condition of solar minimum. Percentage contribution of PEC to GPS TEC indicated the clear dependence from the time with maximal values (more than 50–60%) during night-time and lesser values (25–45%) during day-time.     

A global survey of COSMIC ionospheric peak electron density and its height: A comparison with ground-based ionosonde measurements

With a network of ground-based ionosondes distributed around the world, the ionospheric peak electron density and its height measured by FORMOSAT-3/COSMIC satellites in terms of GPS radio occultation technique are extensively examined in this article. It is found that, in spite of the latitude, the mean values of the peak electron density measured by COSMIC satellites are systematically smaller than those observed by ground-based ionosondes. The discrepancy between them is dependent on the latitude, namely, it is small in low and mid-latitudes and large in high-latitude region. Moreover, statistical analysis shows that the slopes of the regression line that is best fitted to the scatter diagram of occultation-retrieved peak electron density (ordinate axis) versus ionosonde-observed peak density (abscissa axis) are universally less than one. This feature is believed to be the result of path average effect of non-uniform distribution of the electron density along the GSP ray during the occultation. A comparison between COSMIC-measured peak height and ionosonde-derived peak height hmF2 indicates that the former is systematically higher than the latter. The difference in the two can be as large as 20% or more in equatorial and low-latitude regions. This result implies that the peak height hmF2 derived from the virtual height through true height analysis based on Titheridge method seems to underestimate the true peak height. The correlation between COSMIC and ionosonde peak electron densities is analyzed and the result reveals that correlation coefficient seems to be dependent on the fluctuation of the occultation-retrieved electron density profile. The correlation will be higher (lower) for the electron density profiles with smaller (larger) fluctuations. This feature suggests that the inhomogeneous distribution of the electron density along the GPS ray path during the occultation plays an important role affecting the correlation between COSMIC and ionosonde measurements.     

Ionospheric responses to the 21 August 2017 great American solar eclipse – A multi-instrument study

Herein, we report on the ionospheric responses to a total solar eclipse that occurred on 21 August 2017 over the US region. Ground-based GPS total electron content (TEC) data along with ground-based measurements (Millstone Hill Observatory (MHO) and digital ionosondes) and space-based measurements (COSMIC radio occultation (RO) technique) allowed us to identify eclipse-associated ionospheric responses. TEC data at ~20°, ~30°, and ~40°N latitudes from the west to east longitudes show not only considerable depression but also wave-like characteristics in TEC both in the path of totality and away from it, exclusively on the day of eclipse. Interestingly, the observed depressions are associated with lesser (higher) magnitudes at stations over which the solar obscuration percentage was meager (significant), a clear indication of bow-wave-like features. The MHO observes a 30% reduction in F₂-layer electron densities between 180 and 220 km on eclipse day. Ionosonde-scaled parameters over Boulder (40.4°N, 100°E) and Austin (30.4°N, 94.4°E) show a significant decrease in critical frequencies while an altitude elevation is seen in the virtual heights of the F-layer only during the eclipse day and that decreases are associated with wave-like signatures, which could be attributed to eclipse-generated waves. The estimated vertical electron density profile from the COSMIC RO-based technique shows a maximum depletion of 40%. Relatively intense and moderate depths of TEC depression, considerable reductions in the F₂-layer electron densities measured by the MHO and COSMIC RO-measured densities at the F₂-layer peak, and elevations in virtual heights and reduction in the critical frequencies measured by ionosondes during the eclipse day could be due to the eclipse-induced dynamical effects such as gravity waves (GWs) and their associated electro-dynamical effects (modification of ionospheric electric fields due to GWs).     

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.     

Electron density retrieval from occulting GNSS signals using a gradient-aided inversion technique

In the coming years, opportunities for remote sensing of electron density in the Earth’s ionosphere will expand with the advent of Galileo, which will become part of the global navigation satellite system (GNSS). Methods for accurate electron density retrieval from radio occultation data continue to improve. We describe a new method of electron density retrieval using total electron content measurements obtained in low Earth orbit. This method can be applied to data from dual-frequency receivers tracking the GPS or Galileo transmitters. This simulation study demonstrates that the method significantly improves retrieval accuracy compared to the standard Abel inversion approach that assumes a spherically symmetric ionosphere. Our method incorporates horizontal gradient information available from global maps of Total Electron Content (TEC), which are available from the International GNSS Service (IGS) on a routine basis. The combination of ground and space measurements allows us to improve the accuracy of electron density profiles near the occultation tangent point in the E and F regions of the ionosphere.     

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

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

EUV-TEC proxy to describe ionospheric variability using satellite-borne solar EUV measurements: First results

Primary photoionisation of major ionospheric constituents is calculated from satellite-borne solar EUV measurements. Number densities of the background atmosphere are taken from the NRLMSISE-00 climatology. From the calculated ionisation rates, a proxy termed EUV-TEC, which is based on the global total ionisation is calculated, and describes the ionospheric response to solar EUV and its variability. The proxy is compared against the global mean ionospheric total electron content (TEC) derived from GPS data. Results show that the EUV-TEC proxy provides a better overall representation of global TEC than conventional solar indices like F10.7 do. The EUV-TEC proxy may be used for scientific research, and to describe the ionospheric effects on radio communication and navigation systems.     

一种用于电离层TEC监测的GNSS信号载波跟踪算法

全球卫星导航系统（GNSS）是电离层TEC监测中应用最普遍的手段. 目前方法通常是在传统导航用途的GNSS接收机输出的原始观测量基础上,经过数据后处理得到电离层TEC信息,其GNSS信号的跟踪处理算法依然采用GNSS导航接收机的算法. 针对GNSS系统用于电离层TEC监测的特殊性,提出一种称为GNSS双频信号和差联合跟踪的新算法,与传统方法相比,该算法直接跟踪电离层TEC的变化,可以提高电离层TEC跟踪的灵敏度和TEC的观测精度,改善电离层TEC监测性能.      

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.     

Radio occultation experiments with Comet Halley

Previous radio occultation investigations on cometary comae and tails have included refraction measurements and intensity scintillations of natural radio sources used to derive the density and structure of the cometary plasma. Significant improvements in the coverage and sensitivity of these measurements will be achieved during the present apparition of Comet Halley. The comet missions GIOTTO and VEGA will also feature passive radio science expeirments designed to measure comet-induced Doppler shifts of the dual-frequency spacecraft signals during Halley flyby. A brief survey of these radio occultation measurement techniques and their application in the specific case of Comet Halley are presented.     

南极GPS定位中的电离层延迟分析

基于南极地区国际GNSS服务组织(IGS,International GNSS Service)跟踪站的全球定位系统(GPS,Global Position System)双频实测数据,分析了南极地区电离层延迟的变化情况及其二阶项延迟对南极GPS定位结果的影响.结果表明:南极地区的总电子含量(TEC,Total Electron Content)日间波动频繁,其日间TEC最大值变化较中纬度地区剧烈;在南极地区夏季,电离层二阶项延迟对GPS定位结果的影响可达cm级.同时,由于欧洲定轨中心(CODE,The European Center for Orbit Determination)提供的全球电离层模型(GIM,Global Ionosphere Maps)在南极区域应用的局限性,通过选取南极地区6个IGS跟踪站作为基准站建立了区域电离层TEC模型(RIM,Regional Ionosphere Model).经实测数据计算证明,对于南极地区,RIM的定位精度在一定程度上优于全球电离层模型GIM.     

不同太阳活动条件下电离层形态对估算GPS系统硬件延迟的影响

利用两个中纬度台站GPS观测数据提取的GPS卫星硬件延迟,分析了不同太阳活动情况下估算的硬件延迟稳定性和统计特征,结合同期电离层观测数据,研究了电离层状态对硬件延迟估算结果的影响.研究结果表明,基于太阳活动高年(2001年)GPS观测数据估算的硬件延迟稳定性要低于太阳活动低年GPS观测数据的估算结果,利用2001年GPS数据估算的卫星硬件延迟年标准偏差(RMS)平均值约为1TECU,而2009年GPS数据估算的卫星硬件延迟年标准偏差平均值约为0.8TECU.通过对2001年和2009年北京地区电离层F₂层最大电子密度(N_mF₂)变化性分析,结合GPS硬件延迟估算方法对电离层时空变化条件的要求,认为硬件延迟稳定性与太阳活动强度的联系是由不同太阳活动条件下电离层变化的强度差异引起的.      

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.