Statistics of GPS ionospheric scintillation and irregularities over polar regions at solar minimum

A statistical study of the occurrence characteristic of GPS ionospheric scintillation and irregularity in the polar latitude is presented. These measurements were made at Ny-Alesund, Svalbard [78.9°N, 11.9°E; 75.8°N corrected geomagnetic latitude (CGMLat)] and Larsemann Hills, East Antarctica (69.4°S, 76.4°E; 74.6°S CGMLat) during 2007–2008. It is found that the GPS phase scintillation and irregularity activity mainly takes place in the months 10, 11 and 12 at Ny-Alesund, and in the months 5, 6 at Larsemann Hills. The seasonal pattern of phase scintillation with respect to the station indicates that the GPS phase scintillation occurrence is a local winter phenomenon, which shows consistent results with past studies of 250 MHz satellite beacon measurements. The occurrence rates of GPS amplitude scintillation at the two stations are below 1%. A comparison with the interplanetary magnetic field (IMF) B _y and B _z components shows that the phase scintillation occurrence level is higher during the period from later afternoon to sunset (16–19 h) at Ny-Alesund, and from sunset to pre-midnight (18–23 h) at Larsemann Hills for negative IMF components. The findings seem to indicate that the dependence of scintillation and irregularity occurrence on geomagnetic activity appears to be associated with the magnetic local time (MLT).     

Contribution of ionospheric irregularities to the error of dual-frequency GNSS positioning

This paper investigates the third-order residual range error in the dual-frequency correction of ionospheric effects on satellite navigation. We solve the two-point trajectory problem using the perturbation method to derive second-approximation formulas for the phase path of the wave propagating through an inhomogeneous ionosphere. It is shown that these formulas are consistent with the results derived from applying perturbation theory directly to the eikonal equation. The resulting expression for the phase path is used in calculating the residual range error of dual-frequency global positioning system (GPS) observations, in view of second- and third-order terms. The third-order correction includes not only the quadratic correction of the refractive index but also the correction for ray bending in an inhomogeneous ionosphere. Our calculations took into consideration that the ionosphere has regular large-scale irregularities, as well as smaller-scale random irregularities. Numerical examples show that geomagnetic field effects, which constitute a second-order correction, typically exceed the effects of the quadratic correction and the regular ionospheric inhomogeneity. The contribution from random irregularities can compare with or exceed that made by the second-order correction. Therefore, random ionospheric irregularities can make a significant (sometimes dominant) contribution to the residual range error.     

GPS scintillation over the European Arctic during the November 2004 storms

Small-scale irregularities in the background electron density of the ionosphere can cause rapid fluctuations in the amplitude and phase of radio signals passing through it. These rapid fluctuations are known as scintillation and can cause a Global Positioning System (GPS) receiver to lose lock on a signal. This could compromise the integrity of a safety of life system based on GPS, operating in auroral regions. In this paper, the relationship between the loss of lock on GPS signals and ionospheric scintillation in auroral regions is explored. The period from 8 to 14 November 2004 is selected for this study, as it includes both geomagnetically quiet and disturbed conditions. Phase and amplitude scintillation are measured by GPS receivers located at three sites in Northern Scandinavia, and correlated with losses of signal lock in receivers at varying distances from the scintillation receivers. Local multi-path effects are screened out by rejection of low-elevation data from the analysis. The results indicate that losses of lock are more closely related to rapid fluctuations in the phase rather than the amplitude of the received signal. This supports the idea, suggested by Humphreys et al. (2005) (performance of GPS carrier tracking loops during ionospheric scintillations. Proceedings Internationsl Ionospheric Effects Symposium 3–5 May 2005), that a wide loop bandwidth may be preferred for receivers operating at auroral latitudes. Evidence from the Imaging Riometer for Ionospheric Studies (IRIS) appears to suggest that, for this particular storm, precipitation of particles in the D/E regions may be the mechanism that drives the rapid phase fluctuations in the signal.

GPS differential code biases considering the second-order ionospheric term

GPS Solutions - Differential code biases (DCBs) of the global positioning system (GPS) are generally estimated together with total electron content (TEC) along the signal transmission path through...     

Different mechanisms of the ionospheric influence on GPS occultation signals

A local mechanism for strong ionospheric effects on radio occultation (RO) global positioning satellite system (GPS) signals is described. Peculiar zones centered at the critical points (the tangent points) in the ionosphere, where the gradient of the electron density is perpendicular to the RO ray trajectory, strongly influence the amplitude and phase of RO signals. It follows from the analytical model of local ionospheric effects that the positions of the critical points depend on the RO geometry and the structure of the ionospheric disturbances. Centers of strong ionospheric influence on RO signals can exist, for example, in the sporadic E-layers, which are inclined by 3–6° relative to the local horizontal direction. Also, intense F2 layer irregularities can contribute to the RO signal variations. A classification of the ionospheric influence on the GPS RO signals is introduced using the amplitude data, which indicates different mechanisms (local, diffraction, etc.) for radio waves propagation. The existence of regular mechanisms (e.g., local mechanism) indicates a potential for separating the regular and random parts in the ionospheric influence on the RO signals.     

GPS observations of the ionospheric F2-layer behavior during the 20th November 2003 geomagnetic storm over South Korea

The ionospheric F2-layer peak density (NmF2) and its height (hmF2) are of great influence on the shape of the ionospheric electron density profile Ne (h) and may be indicative of other physical processes within the ionosphere, especially those due to geomagnetic storms. Such parameters are often estimated using models such as the semiempirical international reference ionosphere (IRI) models or are measured using moderately priced to expensive instrumentation, such as ionosondes or incoherent scatter radars. Global positioning system (GPS) observations have become a powerful tool for mapping high-resolution ionospheric structures, which can be used to study the ionospheric response to geomagnetic storms. In this paper, we describe how 3-D ionospheric electron density profiles were produced from data of the dense permanent Korean GPS network using the tomography reconstruction technique. These profiles are verified by independent ionosonde data. The responses of GPS-derived parameters at the ionospheric F2-layer to the 20th November 2003 geomagnetic storm over South Korea are investigated. A fairly large increase in the electron density at the F2-layer peak (the NmF2) (positive storm) has been observed during this storm, which is accompanied by a significant uplift in the height of the F2 layer peak (the hmF2). This is confirmed by independent ionosonde observations. We suggest that the F2-layer peak height uplift and NmF2 increase are mainly associated with a strong eastward electric field, and are not associated with the increase of the O/N₂ ratio obtained from the GUVI instruments aboard the TIMED satellite. It is also inferred that the increase in NmF2 is not caused by the changes in neutral composition, but is related to other nonchemical effects, such as dynamical changes of vertical ion motions induced by winds and E × B drifts, tides and waves in the mesosphere/lower thermosphere region, which can be dynamically coupled upward to generate ionospheric perturbations and oscillations.     

Monitoring Seasonal Variations of Ionospheric TEC Using GPS Measurements

The regional ionospheric model is adopted to determine satellite-plus-receiver differential delay. The satellite-plus-receiver differential delay is estimated as constant values for each day. Dual-frequency GPS pseudo-ranges observables are used to compute vertical TEC (VTEC). All the monthly mean VTEC profiles are represented by graphs using GPS data of the Beijing IGS site between 2000 and 2004. The monthly averaged values and amplitudes of VTEC are also represented by graphs. The results indicate that the VTEC has seasonal dependency. The monthly averaged values and amplitudes of VTEC in 2000 are about 2 times larger than that in 2004. The maximum VTEC values are observed in March and April, while the minimum VTEC values are observed in December. The seasonal variations trend is found to be similar after polynomial fitting between 2000 and 2004.     

Determination of the ionospheric foF2 using a stand-alone GPS receiver

The critical frequency of ionospheric F2 layer (foF2) is a measure of the highest frequency of radio signal that may be reflected back by the F2 layer, and it is associated with ionospheric peak electron density in the F2 layer. Accurate long-term foF2 variations are usually derived from ionosonde observations. In this paper, we propose a new method to observe foF2 using a stand-alone global positioning system (GPS) receiver. The proposed method relies on the mathematical equation that relates foF2 to GPS observations. The equation is then implemented in the Kalman filter algorithm to estimate foF2 at every epoch of the observation (30-s rate). Unlike existing methods, the proposed method does not require any additional information from ionosonde observations and does not require any network of GPS receivers. It only requires as inputs the ionospheric scale height and the modeled plasmaspheric electron content, which practically can be derived from any existing ionospheric/plasmaspheric model. We applied the proposed method to estimate long-term variations of foF2 at three GPS stations located at the northern hemisphere (NICO, Cyprus), the southern hemisphere (STR1, Australia) and the south pole (SYOG, Antarctic). To assess the performance of the proposed method, we then compared the results against those derived by ionosonde observations and the International Reference Ionosphere (IRI) 2012 model. We found that, during the period of high solar activity (2011–2012), the values of absolute mean bias between foF2 derived by the proposed method and ionosonde observations are in the range of 0.2–0.5 MHz, while those during the period of low solar activity (2009–2010) are in the range of 0.05–0.15 MHz. Furthermore, the root-mean-square-error (RMSE) values during high and low solar activities are in the range of 0.8–0.9 MHz and of 0.6–0.7 MHz, respectively. We also noticed that the values of absolute mean bias and RMSE between foF2 derived by the proposed method and the IRI-2012 model are slightly larger than those between the proposed method and ionosonde observations. These results demonstrate that the proposed method can estimate foF2 with a comparable accuracy. Since the proposed method can estimate foF2 at every epoch of the observation, it therefore has promising applications for investigating various scales (from small to large) of foF2 irregularities.     

基于CORS数据反演江苏上空电离层电子密度

本文利用IRI2007模型和江苏CORS网数据,结合附加约束的同时迭代重构算法,反演了地磁平静日江苏上空电子密度的结构分布。结果表明,整个研究区域电子密度在不同纬度和不同高度上存在明显的日变特征,电子密度最大值出现在北京时13～15时,随着时间推移,峰值幅度衰减;白天出现明显的E层,夜晚消失;反演结果说明了电离层三维层析技术为监测电离层时空结构提供了一种强有力的实验支持。     

Understanding the role of the ionospheric delay in single-point single-epoch GPS coordinates

The ionospheric delay is the main source of error for single-point single-epoch (SPSE) GPS positioning when using single-frequency receivers. In contrast to the common slant approach, in this article we focus on its effect in final coordinates through the study of bias propagation in SPSE positioning: we first show an analytical resolution for the propagation problem with highly symmetric satellite configurations. To overcome some of the disadvantages of this first method, we use Santerre’s technique and, finally, present a new numerical methodology that allows us to generalize for a real geometry and obtain an average ionospheric positioning error over a given site. From the results obtained, four working hypotheses that relate the ionospheric shape above the receiver with final position errors are presented and tested. These four hypotheses, which agree with average ionospheric positioning error in 95% of the studied cases, can be related to the construction of the design matrix. Finally, these hypotheses have been used to address a situation where the ionospheric delay is corrected with an ionospheric model.     

On the accuracy of the GPS L2 observable for ionospheric monitoring

The introduction of the unencrypted global positioning system (GPS) L2 civil (L2C) signal has the potential to improve measurements made with the L2 frequency, an important observable in GPS-based ionospheric research and monitoring. Recent work has shown significant differences between the legacy L2P(Y) and L2C-derived total electron content rate of change index (ROTI). This difference is observed between L2P(Y) and L2C-derived ROTI with certain receiver models and between zero-baseline receiver pairs. We discuss the likely cause for these differences: L1-aided tracking used to track both the L2P(Y) and L2C signals. We also present L2C data that are confirmed to be from tracking independent of L1. Using the ionospheric-free linear combination, we show that the independently tracked carrier phase dynamics are significantly more accurate than the L1-aided observables. This result is confirmed by comparing the behavior of the L2C and L2P(Y) carrier phase observables upon a sudden antenna rotation.     

GPS反演上电离层剖面及实例

介绍了在利用GPS观测电离层TEC和反演上电离层剖面过程中，如何选择最适宜GPS星，如何计算卫星有关参数，如何计算卫星到测站间斜向TEC等问题，最后给出了一个剖面反演实例。     

The ionospheric eclipse factor method (IEFM) and its application to determining the ionospheric delay for GPS

A new method for modeling the ionospheric delay using global positioning system (GPS) data is proposed, called the ionospheric eclipse factor method (IEFM). It is based on establishing a concept referred to as the ionospheric eclipse factor (IEF) λ of the ionospheric pierce point (IPP) and the IEF’s influence factor (IFF) . The IEF can be used to make a relatively precise distinction between ionospheric daytime and nighttime, whereas the IFF is advantageous for describing the IEF’s variations with day, month, season and year, associated with seasonal variations of total electron content (TEC) of the ionosphere. By combining λ and with the local time t of IPP, the IEFM has the ability to precisely distinguish between ionospheric daytime and nighttime, as well as efficiently combine them during different seasons or months over a year at the IPP. The IEFM-based ionospheric delay estimates are validated by combining an absolute positioning mode with several ionospheric delay correction models or algorithms, using GPS data at an international Global Navigation Satellite System (GNSS) service (IGS) station (WTZR). Our results indicate that the IEFM may further improve ionospheric delay modeling using GPS data.     

不同GPS掩星电离层剖面产品相关性分析

将在一定时空限定范围内的不同低轨卫星COSMIC、GRACE、CHAMP、FY3C的电离层掩星电子密度剖面定义为一个掩星对来对比分析不同类型掩星电离层产品。结果表明：COSMIC掩星对之间的电子密度剖面整体轮廓符合得很好,电子密度剖面主要在250 km以下和500 km以上存在较大的偏差,250~500 km的电子密度整体偏差较小,统计得到的COSMIC掩星对的电子密度参量NmF2和hmF2的相关系数能分别达到0.99和0.97,具有高度相关性,不同COSMIC卫星之间没有明显的系统误差;COSMIC、GRACE、CHAMP和FY3C不同低轨卫星间的电子密度剖面略有差异,通过统计电子密度参量NmF2和hmF2之间的相关系数,COSMIC和CHAMP的相关系数分别为0.95和0.86,COSMIC和GRACE的相关系数分别为0.98和0.94,COSMIC和FY3C的相关系数分别为0.96和0.92,不同掩星类型之间的电子密度参量之间也具有高度相关性,验证了不同卫星任务GPS掩星电离层剖面的一致性。     

Comparison of 4D tomographic mapping versus thin-shell approximation for ionospheric delay corrections for single-frequency GPS receivers over North America

The majority of navigation satellite receivers operate on a single frequency. They compensate for the ionospheric delay using either an ionospheric model which typically only corrects for 50% of the delay or a thin-shell map of the ionosphere. A 4D tomographic imaging technique is used to map the free electron density over the full-height of the ionosphere above North America during autumn 2003. The navigation solutions computed using correction based upon the thin-shell and the full-height maps are compared in this paper. The maps are used to calculate the excess propagation delay on the L1 frequency experienced by GPS receivers at selected locations across North America. The excess delay is applied to correct the single-frequency pseudorange observations at each location, and the improvements to the resulting positioning are calculated. It is shown that the thin-shell and full-height maps perform almost as well as a dual-frequency carrier-smoothed benchmark and for most receivers better than the unfiltered dual-frequency benchmark. The full-height corrections perform well and are considerably better than thin-shell corrections under extreme storm conditions.     

Impact of the Halloween 2003 ionospheric storm on kinematic GPS positioning in Europe

Using dual-frequency data from 36 GPS stations from the EUREF Permanent Network (EPN), the influence of the October 30, 2003 Halloween geomagnetic storm on kinematic GPS positioning is investigated. The Halloween storm induced ionospheric disturbances above the northern part of Europe and Scandinavia. It is shown that kinematic position repeatabilities for this period are mainly affected for stations in northern Europe with outliers reaching 12 cm in the horizontal, and 26 cm in the vertical. These magnitudes are shown to be possibly due to the second-order ionospheric delays on GPS signals, not accounted for in the kinematic GPS positioning analysis performed. In parallel, we generate hourly TEC (Total Electron Content) maps on a 1° × 1° grid using the dense EPN network. These TEC maps do not use any interpolation but provide a high resolution in the time and space and therefore allow to better evidence small structures in the ionosphere than the classical 2-hourly 2.5° × 5° grid Global Ionospheric TEC Maps (GIM). Using the hourly 1° × 1° TEC maps, we reconstruct and refine exactly the zones of intense ionosphere activity during the storm, and we show the correlation between the ionospheric activity and assess the quality of GPS-based kinematic positioning performed in the European region.     

Further observations of GPS satellite oscillator anomalies mimicking ionospheric phase scintillation

An incident has previously been reported where the signal from the Navstar 43 Global Positioning System (GPS) satellite contained phase anomalies in such a way as to mimic ionospheric scintillation. We have observed another 25 events from the same satellite, plus events from three more satellites. Our data includes simultaneous observations from widely spaced receivers (up to 6,590 km apart), from different manufacturers, further ruling out the possibility of local effects. Two of the events involved a satellite (GPS IIF SV-2) broadcasting the L2C signal. This signal contained phase deviations matching those of the L1 signal, but with a 120/154 multiplicative factor. This rules out the possibility of a genuine ionospheric scintillation event, as it does not match the plasma dispersion relation. It does, however, agree with what can be expected from an anomaly in the satellite’s oscillator. While the previously reported event could be dismissed as a freak occurrence, it is now apparent that these events are a persistent phenomenon. They have the potential to corrupt geophysical research with false data and to generate false alarms in systems to forewarn of GPS outages due to scintillation.