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.     

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.     

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.     

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.     

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.     

Calibration errors in determining slant Total Electron Content (TEC) from multi-GNSS data

The global navigation satellite system (GNSS) is presently a powerful tool for sensing the Earth's ionosphere. For this purpose, the ionospheric measurements (IMs), which are by definition slant total electron content biased by satellite and receiver differential code biases (DCBs), need to be first extracted from GNSS data and then used as inputs for further ionospheric representations such as tomography. By using the customary phase-to-code leveling procedure, this research comparatively evaluates the calibration errors on experimental IMs obtained from three GNSS, namely the US Global Positioning System (GPS), the Chinese BeiDou Navigation Satellite System (BDS), and the European Galileo. On the basis of ten days of dual-frequency, triple-GNSS observations collected from eight co-located ground receivers that independently form short-baselines and zero-baselines, the IMs are determined for each receiver for all tracked satellites and then for each satellite differenced for each baseline to evaluate their calibration errors. As first derived from the short-baseline analysis, the effects of calibration errors on IMs range, in total electron content units, from 1.58 to 2.16, 0.70 to 1.87, and 1.13 to 1.56 for GPS, Galileo, and BDS, respectively. Additionally, for short-baseline experiment, it is shown that the code multipath effect accounts for their main budget. Sidereal periodicity is found in single-differenced (SD) IMs for GPS and BDS geostationary satellites, and the correlation of SD IMs over two consecutive days achieves the maximum value when the time tag is around 4?min. Moreover, as byproducts of zero-baseline analysis, daily between-receiver DCBs for GPS are subject to more significant intra-day variations than those for BDS and Galileo.     

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.     

Performance assessment of real-time orbit determination for the Haiyang-2D using onboard BDS-3/GPS observations

The Global Navigation Satellite System (GNSS) receivers equipped on the Haiyang-2D (HY-2D) satellite is capable of tracking the signals of both the third generation of BeiDou satellite navigation System (BDS-3) and the Global Positioning System (GPS), which make it feasible to assess the performance of real-time orbit determination (RTOD) for the HY-2D using onboard GNSS observations. In this study, the achievable accuracy and convergence time of RTOD for the HY-2D using onboard BDS-3 and GPS observations are analyzed. Benefiting from the binary-offset-carrier (BOC) modulation, the BDS-3 C1X signal includes less noise than the GPS C1C signal, which has the same signal frequency and chipping rate. The root mean squares (RMS) of the noises of C1X and C1C code measurements are 0.579 m and 1.636 m, respectively. Thanks to a ten-times higher chipping rate, the code measurements of BDS-3 C5P, GPS C1W and C2W are less noisy. The RMS of code noises of BDS-3 C5P, GPS C1W, and C2W are 0.044 m, 0.386 m, and 0.272 m, respectively. For the HY-2D orbit, the three-dimensional (3D) and radial accuracies can reach 31.8 cm and 7.5 cm with only BDS-3 observations, around 50 % better than the corresponding accuracies with GPS. Better performance of the BDS-3 in RTOD for the HY-2D is attributed to the high quality of its broadcast ephemeris. When random parameters are used to absorb ephemeris errors, substantial improvement is seen in the accuracy of HY-2D orbit with either BDS-3 or GPS. The 3D RMS of HY-2D orbit errors with BDS-3 and GPS are enhanced to 23.1 cm and 33.6 cm, and the RMS of the radial components are improved to 6.1 cm and 13.3 cm, respectively. The convergence time is 41.6 and 75.5 min for the RTOD with BDS-3 and GPS, while it is reduced to 39.2 and 27.4 min after the broadcast ephemeris errors are absorbed by random parameters. Overall, the achievable accuracy of RTOD with BDS-3 reaches decimeter level, which is even better than that with GPS, making real-time navigation using onboard BDS-3 observations a feasible choice for future remote sensing missions.     

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.     

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.     

Multi-wavelength coordinated observations of ionospheric irregularity structures from an anomaly crest location during unusual solar minimum of the 24th cycle

The present paper reports coordinated ionospheric irregularity measurements at optical as well as GPS wavelengths. Optical measurements were obtained from Tiny Ionospheric Photometer (TIP) sensors installed onboard the Constellation Observing System for Meteorology, Ionosphere, and Climate (COSMIC) satellites. GPS radio signals were obtained from a dual frequency GPS receiver operational at Calcutta (22.58°N, 88.38°E geographic; geomagnetic dip: 32.96°; 13.00°N, 161.63°E geomagnetic) under the SCIntillation Network Decision Aid (SCINDA) program. Calcutta is located near the northern crest of Equatorial Ionization Anomaly (EIA) in the Indian longitude sector. The observations were conducted during the unusually low and prolonged solar minima period of 2008–2010. During this period, four cases of post-sunset GPS scintillation were observed from Calcutta. Among those cases, simultaneous fluctuations in GPS Carrier-to-Noise ratios (C/N_o) and measured radiances from TIP over a common ionospheric volume were observed only on February 2, 2008 and September 25, 2008. Fluctuations observed in measured radiances (maximum 0.95 Rayleigh) from TIP due to ionospheric irregularities were found to correspond well with C/N₀ fluctuations on the GPS links observed from Calcutta, such effects being noted even during late evening hours of 21:00–22:00 LT from locations around 40° magnetic dip. These measurements indicate the existence of electron density irregularities of scale sizes varying over several decades from 135.6 nm to 300–400 m well beyond the northern crest of the EIA in the Indian longitude sector during late evening hours even in the unusually low solar activity conditions.     

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

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

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.     

Mitigation of receiver biases in ionospheric observables from PPP with ambiguity resolution

PPP (Precise Point Positioning) is a GNSS (Global Navigation Satellite Systems) positioning method that requires SSR (State Space Representation) corrections in order to provide solutions with an accuracy of centimetric level. The so-called RT-PPP (Real-time PPP) is possible thanks to real-time precise SSR products, for orbits and clocks, provided by IGS (International GNSS Service) and its associate analysis centers such as CNES (Centre National d'Etudes Spatiales). CNES SSR products also enable RT-PPP with integer ambiguity resolution. In GNSS related literature, PPP with ambiguity resolution (PPP-AR) in real-time is often referred as PPP-RTK (PPP – Real Time Kinematic). PPP-WIZARD (PPP - With Integer and Zero-difference Ambiguity Resolution Demonstrator) is a software that is made available by CNES. This software is capable of performing PPP-RTK. It estimates slant ionospheric delays and other GNSS positioning parameters. Since ionospheric effects are spatially correlated by GNSS data from active networks, it is possible to model and provide ionospheric delays for any position in the network coverage area. The prior knowledge ionospheric delays can reduce positioning convergence for PPP-RTK users. Real-time ionospheric models could benefit from highly precise ionospheric delays estimated in PPP-AR. In this study, we demonstrate that ionospheric delays obtained throughout PPP-AR estimation are actu ally ionospheric observables. Ionospheric observables are biased by an order of few meters caused by the receiver hardware biases. These biases prohibit the use of PPP-WIZARD ionospheric delays to produce ionospheric models. Receiver biases correction is essential to provide ionospheric delays while using PPP-AR based ionospheric observables. In this contribution, a method was implemented to estimate and mitigate receiver hardware biases influence on slant ionospheric observables from PPP-AR. In order to assess the proposed approach, PPP-AR data from 12 GNSS stations were processed over a two-month period (March and April 2018). A comparison between IGS ionospheric products and PPP-AR based ionospheric observables corrected for receiver biases, resulted in a mean of differences of −39 cm and 51 cm standard deviation. The results are consistent with the accuracy of the IGS ionospheric products, 2–8 TECU, considering that 1 TECU is ~16 cm in L1. In another analysis, a comparison of ionospheric delays from 5 pairs of short baselines GNSS stations found an agreement of 0.001 m in mean differences with 22 cm standard deviation after receiver biases were corrected. Therefore, the proposed solution is promising and could produce high quality (1–2 TECU) slant ionospheric delays. This product can be used in a large variety of modeling approaches, since ionospheric delays after correction are unbiased. These results indicate that the proposed strategy is promising, and could benefit applications that require accuracy of 1–2 TECU (~16–32 cm in L1).     

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.     

Large-scale traveling ionospheric disturbances of auroral origin according to the data of the GPS network and ionosondes

The intensity of large-scale traveling ionospheric disturbances (LS TIDs), registered using measurements of total electron content (TEC) during the magnetic storms on October 29–31, 2003, and on November 7–11, 2004, had been compared with that of local electron density disturbances. The data of TEC measurements at ground-based GPS receivers located near the ionospheric stations and the corresponding values of the critical frequency of the ionospheric F region f_oF2 were used for this purpose. The variations of TEC and f_oF2 were similar for all events mentioned above. The previous assumption that the ionospheric region with vertical extension from 150 to 200 km located near the F-layer maximum mainly contributes to the TEC variations was confirmed for the cases when the electron density disturbance at the F region maximum was not more than 50%. However, this region probably becomes vertically more extended when the electron density disturbance in the ionospheric F region is about 85%.     

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.     

Estimation of QZSS differential code biases using QZSS/GPS combined observations from MGEX

As an important error source in Global Navigation Satellite System (GNSS) positioning and ionospheric modeling, the differential code biases (DCB) need to be estimated accurately, e.g., the regional Quasi-Zenith satellite system (QZSS). In this paper, the DCB of QZSS is estimated by adopting the global ionospheric modeling method based on QZSS/GPS combined observations from Multi-GNSS experiment (MGEX). The performance of QZSS satellite and receiver DCB is analyzed with observations from day of year (DOY) 275–364, 2018. Good agreement between our estimated QZSS satellite DCB and the products from DLR and CAS is obtained. The bias and root mean square (RMS) of DCB are mostly within ±0.3 ns. The day-to-day fluctuation of the DCB time series is less than 0.5 ns with about 96% of the cases for all satellites. However, the receiver DCB is a little less stable than satellite DCB, and their standard deviations (STDs) are within 1.9 ns. The result shows that the stability of the receiver DCBs is not significantly related to the types of receiver or antenna.     

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.     

Dual-frequency carrier smoothed code filtering with dynamical ionospheric delay modeling

In GNSS applications, carrier-smoothed-code is a widely used technique to combine code pseudo-range and carrier phase measurements. A dynamical ionospheric delay modeling method is proposed based on Kalman filter and least-squares theory. The level of the process noise is adaptively tuned along with the real-time KF state estimation, based on the online variance component estimation method. Meanwhile, the correlations of the time differenced carrier phase measurements are considered. This approach avoids overly optimistically evaluating the estimate and improves the transient accuracy of the estimates. A real GPS dataset is employed to check the performance of the proposed method under different conditions. The results show that the new algorithm can model the ionospheric delay variation well with different sampling intervals or even in ionospheric abnormal environment. The positioning accuracy can be confirmed, about 21%, 35% and 16% better are obtained in the N, E, and U direction than raw dataset.