An annual variation analysis of the ionospheric spatial gradient over a regional area for GNSS applications

An ionospheric spatial gradient represents the ionosphere delay difference between different locations, and its variation over a specific area is important for implementing differential GNSS systems. An estimation method for the ionospheric spatial gradient over a small regional area is proposed. A plate map model is implemented for the direct estimation of the gradients. Nine years of GPS data were processed to figure out the annual variation of the mean gradient at the mid-geomagnetic latitude of 30° N. Gradients along the north–south direction have a mean of 0.65 mm/km and follow solar-cycle variations.     

Evaluation of COMPASS ionospheric model in GNSS positioning

As important products of GNSS navigation message, ionospheric delay model parameters are broadcasted for single-frequency users to improve their positioning accuracy. GPS provides daily Klobuchar ionospheric model parameters based on geomagnetic reference frame, while the regional satellite navigation system of China’s COMPASS broadcasts an eight-parameter ionospheric model, COMPASS Ionospheric Model(CIM), which was generated by processing data from continuous monitoring stations, with updating the parameters every 2 h. To evaluate its performance, CIM predictions are compared to ionospheric delay measurements, along with GPS positioning accuracy comparisons. Real observed data analysis indicates that CIM provides higher correction precision in middle-latitude regions, but relatively lower correction precision for low-latitude regions where the ionosphere has much higher variability. CIM errors for some users show a common bias for in-coming COMPASS signals from different satellites, and hence ionospheric model errors are somehow translated into the receivers’ clock error estimation. In addition, the CIM from the China regional monitoring network are further evaluated for global ionospheric corrections. Results show that in the Northern Hemisphere areas including Asia, Europe and North America, the three-dimensional positioning accuracy using the CIM for ionospheric delay corrections is improved by 7.8%–35.3% when compared to GPS single-frequency positioning ionospheric delay corrections using the Klobuchar model. However, the positioning accuracy in the Southern Hemisphere is degraded due apparently to the lack of monitoring stations there.     

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.     

Estimation and analysis of GPS receiver differential code biases using KGN in Korean Peninsula

The total electron content (TEC) estimation by the Global Positioning System (GPS) can be seriously affected by the differential code biases (DCB), referred to as inter-frequency biases (IFB), of the satellite and receiver so that an accuracy of GPS–TEC value is dependent on the error of DCBs estimation. In this paper, we proposed the singular value decomposition (SVD) method to estimate the DCB of GPS satellites and receivers using the Korean GPS network (KGN) in South Korea. The receiver DCBs of about 49 GPS reference stations in KGN were determined for the accurate estimation of the regional ionospheric TEC. They obtained from the daily solution have large biases ranging from +5 to +27 ns for geomagnetic quiet days. The receiver DCB of SUWN reference station was compared with the estimates of IGS and JPL global ionosphere map (GIM). The results have shown comparatively good agreement at the level within 0.2 ns. After correction of receiver DCBs and knowing the satellite DCBs, the comparison between the behavior of the estimated TEC and that of GIMs was performed for consecutive three days. We showed that there is a good agreement between KASI model and GIMs.     

用GPS 观测研究电离层TEC 水平梯度

双频GPS 用户能自动修正电离层总电子含量(TEC) 引起的延时误差, 但是对于电离层中的不规则体造成的信号闪烁而引起的误差则不能消除. 即使是差分GPS 系统, 电离层误差仍然是其主要的误差源, 其中电离层TEC 梯度将会影响到系统的定位精度和性能. 本文用GPS 方法研究了电离层TEC 的水平梯度问题, 用处于赤道异常区NTUS 台站的GPS 观测数据作了具体计算. 结果表明, 在日落以后到子夜前后电离层垂直TEC 出现了大的涨落, 电离层中的不规则体导致L 波段信号强的闪烁, 同时还伴随着大而快速变化的电离层~TEC 水平梯度. 对比发现, ROTI指数、电离层TEC 水平梯度和电离层垂直TEC 三者之间有很好的对应关系, 它们的变化特征均由电离层中的不规则体引起. 我们认为研究电离层闪烁, 特别是在缺乏S4指数时, 电离层TEC 梯度也可以作为一个重要的可选参数.      

基于CIT技术的广域增强系统电离层延迟修正算法

提出了一种基于电离层层析成像(CIT)技术的广域增强系统(WAAS)电离层延迟 修正算法. 该算法利用模式基函数与截断奇异值分解正则化组合的方式, 实现WAAS单频用户的电离层延迟误差修正. 基于中国区域23个广域基准站 和10个用户站的仿真结果分析表明, 传统的网格算法和基于CIT技术的电离 层延迟修正算法的电离层延迟修正精度与太阳活动、昼夜变化及地磁纬度 之间均存在明显的相关性. 基于CIT技术的电离层延迟修正算法精度优于网格算法, CIT算法的平均误差与标准差相比网格算法均有较大幅度的下降.      

Medium-scale traveling ionospheric disturbances (MSTID) modeling using a dense German GPS network

The MSTIDs are wave-like perturbations of the ionospheric plasma, which cause the most common ionospheric disturbances in mid-latitude regions. Generally the MSTIDs have velocities of several hundred meters per second and wavelengths of several hundred kilometers. The wave-like effect of the MSTID is one of the main obstacles for accurate interpolation of ionospheric corrections in a medium-scale reference GPS network. In this paper we show a new method of detecting and modeling MSTIDs using dense German GPS network. The between-epoch single difference ionospheric delays from a medium scale dense GPS network are used to estimate the parameter of the MSTID e.g. amplitude, wavelength and velocity. The efficiency of the approach is tested with data from about 320 GPS stations in and near Germany. A MSTID wave moving from east to west across Germany was observed at September 27 in 2009. Its wavelength is about 302 km, with a period of ∼7 min and velocity of about 700 m/s.     

Synergistic use of multitemporal RAMP,ICESat and GPS to construct an accurate DEM of the Larsemann Hills region,Antarctica

An enhanced digital elevation model (DEM) of the Larsemann Hills region, east Antarctica, is constructed synergistically by using highly accurate ground-based GPS measurements, satellite-derived laser altimetry (GLAS/ICESat) and Radarsat Antarctic Mapping Project (RAMPv2) DEM-based point elevation dataset. Our DEM has a vertical accuracy of about 1.5 times better than RAMPv2 DEM and seven times better than GLAS/ICESAT-based DEM. The accuracy is improved by validating the RAMPv2 DEM elevation by supplementing with GLAS/ICESat and DGPS survey data, when compared to that of DEM constructed by using GLAS/ICESat or RAMPv2 alone. With the use of accurate GPS data as ground control points reference elevations, the DEM extracted is much more accurate with least mean RMSE of 34.5 m than that constructed by using a combination of GLAS/ICESat and RAMPv2 as true reference. The newly constructed DEM 7 achieves highest accuracy with the least average elevation difference of 0.27 m calculated using 46 ground reference points. Available DEMs of Antarctic region generated by using radar altimetry and the Antarctic Digital Database indicate elevation variations in the range of 50–100 m, which necessitates the generation of local DEM and its validation by using ground truth. This is our first attempt of fusing multi-temporal, multi-sensor and multi-source elevation data to generate a DEM of any part of Antarctica, in order to address the ice elevation change to infer the ice mass balance. Our approach focuses on the strengths of each elevation data source to produce an accurate DEM.     

The use of ionospheric tomography and elevation masks to reduce the overall error in single-frequency GPS timing applications

Signals from Global Positioning System (GPS) satellites at the horizon or at low elevations are often excluded from a GPS solution because they experience considerable ionospheric delays and multipath effects. Their exclusion can degrade the overall satellite geometry for the calculations, resulting in greater errors; an effect known as the Dilution of Precision (DOP). In contrast, signals from high elevation satellites experience less ionospheric delays and multipath effects. The aim is to find a balance in the choice of elevation mask, to reduce the propagation delays and multipath whilst maintaining good satellite geometry, and to use tomography to correct for the ionosphere and thus improve single-frequency GPS timing accuracy. GPS data, collected from a global network of dual-frequency GPS receivers, have been used to produce four GPS timing solutions, each with a different ionospheric compensation technique. One solution uses a 4D tomographic algorithm, Multi-Instrument Data Analysis System (MIDAS), to compensate for the ionospheric delay. Maps of ionospheric electron density are produced and used to correct the single-frequency pseudorange observations. This method is compared to a dual-frequency solution and two other single-frequency solutions: one does not include any ionospheric compensation and the other uses the broadcast Klobuchar model. Data from the solar maximum year 2002 and October 2003 have been investigated to display results when the ionospheric delays are large and variable. The study focuses on Europe and results are produced for the chosen test site, VILL (Villafranca, Spain). The effects of excluding all of the GPS satellites below various elevation masks, ranging from 5° to 40°, on timing solutions for fixed (static) and mobile (moving) situations are presented. The greatest timing accuracies when using the fixed GPS receiver technique are obtained by using a 40° mask, rather than a 5° mask. The mobile GPS timing solutions are most accurate when satellites at lower elevations continue to be included: using a mask between 10° and 20°. MIDAS offers the most accurate and least variable single-frequency timing solution and accuracies to within 10 ns are achieved for fixed GPS receiver situations. Future improvements are anticipated by combining both GPS and Galileo data towards computing a timing solution.     

Empirical model of the TEC response to the geomagnetic activity over the North American region

The paper presents an empirical model of the total electron content (TEC) response to the geomagnetic activity described by the K_p-index. The model is built on the basis of TEC measurements covering the region of North America (50°W–150°W, 10°N–60°N) for the period of time between October 2004 and December 2009. By using a 2D (latitude-time) cross-correlation analysis it is found that the ionospheric response to the geomagnetic activity over the considered geographic region and at low solar activity revealed both positive and negative phases of response. The both phases of the ionospheric response have different duration and time delay with respect to the geomagnetic storm. It was found that these two parameters of the ionospheric response depend on the season and geographical latitude. The presence of two phases, positive and negative, of the ionospheric response imposed the implementation of two different time delay constants in order to properly describe the two different delayed reactions. The seasonal dependence of the TEC response to geomagnetic storms is characterized by predominantly positive response in winter with a short (usually ∼5–6 h) time delay as well as mainly negative response in summer with a long (larger than 15 h) time delay. While the TEC response in March and October is more close to the winter one the response in April and September is similar to the summer one.     

GPS derived TEC and foF2 variability at an equatorial station and the performance of IRI-model

The ionosphere induces a time delay in transionospheric radio signals such as the Global Positioning System (GPS) signal. The Total Electron Content (TEC) is a key parameter in the mitigation of ionospheric effects on transionospheric signals. The delay in GPS signal induced by the ionosphere is proportional to TEC along the path from the GPS satellite to a receiver. The diurnal monthly and seasonal variations of ionospheric electron content were studied during the year 2010, a year of extreme solar minimum (F10.7 = 81 solar flux unit), with data from the GPS receiver and the Digisonde Portable Sounder (DPS) collocated at Ilorin (Geog. Lat. 8.50°N, Long. 4.50°E, dip −7.9°). The diurnal monthly variation shows steady increases in TEC and F2-layer critical frequency (foF2) from pre-dawn minimum to afternoon maximum and then decreases after sunset. TEC show significant seasonal variation during the daytime between 0900 and 1900 UT (LT = UT + 1 h) with a maximum during the March equinox (about 35 TECU) and minimum during the June solstice (about 24 TECU). The GPS-TEC and foF2 values reveal a weak seasonal anomaly and equinoctial asymmetry during the daytime. The variations observed find their explanations in the amount of solar radiation and neutral gas composition. The measured TEC and foF2 values were compared with last two versions of the International Reference Ionosphere (IRI-2007 and IRI-2012) model predictions using the NeQuick and CCIR (International Radio Consultative Committee) options respectively in the model. In general, the two models give foF2 close to the experimental values, whereas significant discrepancies are found in the predictions of TEC from the models especially during the daytime. The error in height dependent thickness parameter, daytime underestimation of equatorial drift and contributions of electrons from altitudes above 2000 km have been suggested as the possible causes.     

Higher order ionospheric propagation effects on GPS radio occultation signals

With the increasing number of remote sensing satellites using the GPS radio occultation technique for atmospheric sounding, the estimation of higher order ionospheric effects and their mitigation have become relevant and important. Due to long ionospheric limb paths, GPS signals are strongly affected by ionospheric refraction during radio occultation. Standard dual-frequency GPS measurements may be used to estimate the first order term of the refractive index. However, non-linear terms such as the second and third order ionospheric terms and ray path bending effects are not considered in occultation measurements so far. Analysing selected CHAMP–GPS occultation events different higher order ionospheric terms are estimated and their effects on dual-frequency range estimation and total electron content (TEC) estimation are discussed. We have found that the separation between the GPS L1 and L2 ray paths exceeds the kilometer level during occultation for a vertical TEC level of more than 160 TEC units. Corresponding errors in the GPS dual-frequency range estimation and TEC estimation are found to exceed the meter and 10 TEC units level, respectively.     

Response of the ionospheric F-region in the Brazilian sector during the super geomagnetic storm in April 2000 observed by GPS

The response of the ionospheric F-region in the equatorial and low latitude regions in the Brazilian sector during the super geomagnetic storm on 06–07 April 2000 has been studied in the present investigation. The geomagnetic storm reached a minimum D_st of −288 nT at 0100 UT on 07 April. In this paper, we present vertical total electron content (VTEC) and phase fluctuations (in TECU/min) from GPS observations obtained at Imperatriz (5.5°S, 47.5°W; IMPZ), Brasília (15.9°S, 47.9°W; BRAZ), Presidente Prudente (22.12°S, 51.4°W; UEPP), and Porto Alegre (30.1°S, 51.1°W; POAL) during the period 05–08 April. Also, several GPS-based TEC maps are presented from the global GPS network, showing widespread and drastic TEC changes during the different phases of the geomagnetic storm. In addition, ion density measurements on-board the satellite Defense Meteorological Satellite Program (DMSP) F15 orbiting at an altitude of 840 km and the first Republic of China satellite (ROCSAT-1) orbiting at an altitude of 600 km are presented. The observations indicate that one of the orbits of the DMSP satellite is fairly close to the 4 GPS stations and both the DMSP F15 ion-density plots and the phase fluctuations from GPS observations show no ionospheric irregularities in the Brazilian sector before 2358 UT on the night of 06–07 April 2000. During the fast decrease of D_st on 06 April, there is a prompt penetration of electric field of magnetospheric origin resulting in decrease of VTEC at IMPZ, an equatorial station and large increase in VTEC at POAL, a low latitude station. This resulted in strong phase fluctuations on the night of 06–07 April, up to POAL. During the daytime on 07 April during the recovery phase, the VTEC observations show positive ionospheric storm at all the GPS stations, from IMPZ to POAL, and the effect increasing from IMPZ to POAL. This is possibly linked to the equatorward directed meridional wind. During the daytime on 08 April (the recovery phase continues), the VTEC observations show very small negative ionospheric storm at IMPZ but the positive ionospheric storm effect is observed from BRAZ to POAL possibly linked to enhancement of the equatorial ionospheric anomaly.     

Ionospheric correction using GPS Klobuchar coefficients with an empirical night-time delay model

We proposed an ionospheric correction approach called NKlob to mitigate the ionospheric delay errors. NKlob is a modification of the original GPS Ionospheric Correction Algorithm (ICA), which uses an empirical night-time model depending on the time, geomagnetic location and periodicities of the ionospheric behavior to replace the night-time constant delay in GPS ICA. Performance of NKlob was evaluated by the independent total electron contents (TECs) derived from Global Ionospheric Maps (GIMs) of the International GNSS Services (IGS) and Jason-2 altimetry satellite during 2013–2017. Compared to GIM TECs, NKlob corrects 51.5% of the ionospheric delay errors, which outperforms GPS ICA by 6.3%. Compared to Jason-2 TECs, NKlob mitigates the ionospheric errors by 58.1%, which is approximately 3.7% better than that of GPS ICA. NKlob shows significant improvement in low-latitude and equatorial regions with respect to GPS ICA, meanwhile exhibiting comparable performance at middle and high latitudes. Since NKlob only requires slight technical changes at the processing level of GPS receivers, we suppose that it can be easily implemented for better ionospheric delay corrections of real-time GPS single-frequency applications.     

Kalman滤波估算电离层延迟的一种优化方法

频间偏差(Inter Frequency Bias,IFB)通常会给电离层延迟的解算带来误差.目前从电离层延迟中消除频间偏差的方法是基于GPS双频观测数据建立垂直电离层模型,利用卡尔曼滤波实时估算电离层模型系数和频间偏差.然而滤波过程中的测量噪声协方差矩阵没有考虑系统观测量之间的相关性,这会导致滤波模型不准确,进而影响...     

Global characteristics of ionospheric electric fields and disturbances during the first hours of magnetic storms

The ionospheric plasma density can be significantly disturbed during magnetic storms. In the conventional scenario of ionospheric storms, the negative storm phases with plasma density decreases are caused by neutral composition changes, and the positive storm phases with plasma density increases are often related to atmospheric gravity waves. However, recent studies show that the global redistribution of the ionospheric plasma is dominated primarily by electric fields during the first hours of magnetic storms. In this paper, we present the measurements of ionospheric disturbances by the DMSP satellites and GPS network during the magnetic storm on 6 April 2000. The DMSP measurements include the F region ion velocity and density at the altitude of ∼840 km, and the GPS receiver network provides total electron content (TEC) measurements. The storm-time ionospheric disturbances show the following characteristics. The plasma density is deeply depleted in a latitudinal range of ∼20° over the equatorial region in the evening sector, and the depletions represent plasma bubbles. The ionospheric plasma density at middle latitudes (20°–40° magnetic latitudes) is significantly increased. The dayside TEC is increased simultaneously over a large latitudinal range. An enhanced TEC band forms in the afternoon sector, goes through the cusp region, and enters the polar cap. All the observed ionospheric disturbances occur within 1–5 h from the storm sudden commencement. The observations suggest that penetration electric fields play a major role in the rapid generation of equatorial plasma bubbles and the simultaneous increases of the dayside TEC within the first 2 h during the storm main phase. The ionospheric disturbances at later times may be caused by the combination of penetration electric fields and neutral wind dynamo process.     

Kalman filtering for GPS/magnetometer integrated navigation system

This paper investigated the data processing method for a GPS/IMU/magnetometer integrated system with Kalman filtering (KF). As a result of GPS/IMU/magnetometer land vehicle system, dead-reckoning of magnetometer and accelerometer integrated subsystem bridged very well the GPS signal outage due to the trees on the two sides of the road. Both differential GPS data processing method and the carrier-phase method with magnetometers’ outputs for predicting the car position, velocity, and acceleration (PVA) are presented. The results from DGPS with Kinematical Positioning (KINPOS) software shown that the averages of the north, east, and down direction standard deviation (short for “std”) are 0.014, 0.010, and 0.018 m, respectively. The std of velocities and accelerations derived by the position and velocity differentiation are 10, 7, 13 mm/s, 7, 5, 9 mm/s², respectively. This method for getting velocities and accelerations requires higher accurate position coordinates. But the position accuracy has frequently been degraded in this case when the car drove under the trees or other similar kinematical environments. That caused the larger velocity and acceleration errors. While the results from the carrier-phase method are std of the velocities = 2.1 mm/s, 1.3 mm/s, 3.7 mm/s in north, east, down, and std of the accelerations = 1.5 mm/s², 0.9 mm/s², 2.3 mm/s² for the static test period; as compared with KINPOS software results, std of the velocity difference between the carrier-phase method and the DGPS method = 7 mm/s, 6.9 mm/s, 9.7 mm/s in north, east, down direction, and std of acceleration difference = 5.0 mm/s², 4.5 mm/s², 7.5 mm/s² in north, east, down direction for the kinematical test period. Obviously, errors come from both the carrier-phase method and DGPS velocity and acceleration results derived directly by the position differentiation. In addition, better accuracy of positions than that before KF has been got by means of velocities and accelerations derived by the carrier-phase method after KF.     

电离层不规则结构对GPS性能的影响

电离层不规则结构的存在可引起无线电信号的幅度和相位发生随机起伏, 这 种电离层闪烁现象会影响全球定位系统(Global Positioning System, GPS)的 性能, 降低定位精度, 严重时导致信号失锁. 电离层不规则结构对GPS性能的 影响涉及电离层物理、接收机设计和表征卫星几何分布的精度衰减因子(Dilution of Positioning, DOP)等多方面因素. 本文通过对表征电离层不规则结构参数 的分析, 根据GPS接收机跟踪环路和闪烁信号模型, 综合研究了电离层闪烁对 GPS接收机载波跟踪环和码跟踪环跟踪误差的影响; 结合实际观测, 评述 了电离层不规则结构对单频和双频GPS接收机定位性能的影响, 在此基础上 提出了有待深入研究的问题及具体建议.      

Comparative analysis of TEC disturbances over tropical cyclone zones in the North–West Pacific Ocean

In this paper, we study ionospheric total electron content (TEC) disturbances associated with tropical cyclones (TCs). The study relies on the statistical analysis of six cyclones of different intensity which occurred in the North–West Pacific Ocean in September–November 2005. We have used TEC data from the international network of two-frequency ground-based GPS receivers and NCEP/NCAR meteorological archive. TEC variations of different period ranges (02–20 and 20–60 min) are shown to be more intense during TC peaks under quiet geomagnetic conditions. The highest TEC variation amplitudes are registered when the wind speed in the cyclone and the TC area are maximum. The intensification of TEC disturbances is more pronounced when several cyclones occur simultaneously. We have revealed that the ionospheric response to TC can be observed only after the cyclone has reached typhoon intensity. The ionospheric response is more pronounced at low satellite elevation angles.     

The calculation of ionospheric absorption with modern computers

New outcomes are proposed for ionospheric absorption starting from the Appleton–Hartree formula, in its complete form. The range of applicability is discussed for the approximate formulae, which are usually employed in the calculation of non-deviative absorption coefficient. These results were achieved by performing a more refined approximation that is valid under quasi-longitudinal (QL) propagation conditions. The more refined QL approximation and the usually employed non-deviative absorption are compared with that derived from a complete formulation. Their expressions, nothing complicated, can usefully be implemented in a software program running on modern computers. Moreover, the importance of considering Booker’s rule is highlighted. A radio link of ground range D = 1000 km was also simulated using ray tracing for a sample daytime ionosphere. Finally, some estimations of the integrated absorption for the radio link considered are provided for different frequencies.