Absolute Positioning with Single-Frequency GPS Receivers

The use of precise post-processed satellite orbits and satellite clock corrections in absolute positioning, using one GPS receiver only, has proven to be an accurate alternative to the more commonly used differential techniques for many applications in georeferencing. The absolute approach is capable of centimeter accuracy when using state-of-the-art, dual-frequency GPS receivers. When using observations from single-frequency receivers, however, the accuracy, especially in height, decreases. The obvious reason for this degradation in accuracy is the effect of unmodeled ionospheric delay. This paper discusses the availability of some empirical ionospheric models that are publicly available and quantifies their usefulness for absolute positioning using single-frequency GPS receivers. The Global Ionospheric Model supplied by International GPS Service (IGS) is the most accurate one and is recommended for absolute positioning using single-frequency GPS receivers. Using high-quality single-frequency observations, a horizontal epoch-to-epoch accuracy of better than 1 m and a vertical accuracy of approximately 1 m is demonstrated. ? 2002 Wiley Periodicals, Inc.     

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.     

Quality assessment of GPS rapid static positioning with weighted ionospheric parameters in generalized least squares

Precise GPS positioning requires the processing of carrier-phase observations and fixing integer ambiguities. With increasing distance between receivers, ambiguity fixing becomes more difficult because ionospheric and tropospheric effects do not cancel sufficiently in double differencing. A popular procedure in static positioning is to increase the length of the observing session and/or to apply atmospheric (ionospheric) models and corrections. We investigate the methodology for GPS rapid static positioning that requires just a few minutes of dual-frequency GPS observations for medium-length baselines. Ionospheric corrections are not required, but the ionospheric delays are treated as pseudo-observations having a priori values and respective weights. The tropospheric delays are reduced by using well-established troposphere models, and satellite orbital and clock errors are eliminated by using IGS rapid products. Several numerical tests based on actual GPS data are presented. It is shown that the proposed methodology is suitable for rapid static positioning within 50–70 km from the closest reference network station and that centimeter-level precision in positioning is feasible when using just 1 min of dual-frequency GPS data.     

Precision spacecraft navigation using a low-cost GPS receiver

Within the PROBA-2 microsatellite mission, a miniaturized single-frequency GPS receiver based on commercial-off-the-shelf (COTS) technology is employed for onboard navigation and timing. A rapid electronic fuse protects against destructive single-event latch-ups (SEL) and enables a quasi-continuous receiver operation despite the inherent sensitivity to space radiation. While limited to single-frequency C/A-code tracking with a narrow-band frontend, the receiver is able to provide precision navigation services through processing of raw GPS measurements on ground as well as a built-in real-time navigation system. In both cases, ionospheric path delays are eliminated through a combination of L1 pseudorange and carrier phase measurements, which also offers a factor-of-two noise reduction relative to code-only processing. By comparison with satellite laser ranging (SLR) measurements, a 0.3-m (3D rms) accuracy is demonstrated for the PROBA-2 reduced dynamic orbit determinations using post-processed GPS orbit and clock products. Furthermore, the experimental onboard navigation system is shown to provide real-time position information with a 3D rms accuracy of about 1?m, which notably outperforms the specification of the Standard Positioning Service (SPS). In view of their lower hardware complexity, mass budget and power requirements as well as the reduced interference susceptibility, legacy C/A-code receivers can thus provide an attractive alternative to dual-frequency receivers even for demanding navigation applications in low Earth orbit.     

An improvement to ionospheric delay correction for single-frequency GPS users – the APR-I scheme

 For the commonly used GPS wide-area augmentation systems (WAAS) with a grid ionospheric model, the efficient modelling of ionospheric delays in real time, for single-frequency GPS users, is still a crucial issue which needs further research. This is particularly necessary when differential ionospheric delay corrections cannot be broadcast, when users cannot receive them, or when there are ionospheric anomalies. Ionospheric delays have a severe effect on navigation performance of single-frequency receivers. A new scheme is proposed which can efficiently address the above problems. The robust recurrence technique is based on the efficient combination of single-frequency GPS observations by users and the high-precision differential ionospheric delay corrections from WAAS. Its effectiveness is verified with examples. Received: 24 December 1999 / Accepted 21 February 2001     

Joint estimation of vertical total electron content (VTEC) and satellite differential code biases (SDCBs) using low-cost receivers

Vertical total electron content (VTEC) parameters estimated using global navigation satellite system (GNSS) data are of great interest for ionosphere sensing. Satellite differential code biases (SDCBs) account for one source of error which, if left uncorrected, can deteriorate performance of positioning, timing and other applications. The customary approach to estimate VTEC along with SDCBs from dual-frequency GNSS data, hereinafter referred to as DF approach, consists of two sequential steps. The first step seeks to retrieve ionospheric observables through the carrier-to-code leveling technique. This observable, related to the slant total electron content (STEC) along the satellite–receiver line-of-sight, is biased also by the SDCBs and the receiver differential code biases (RDCBs). By means of thin-layer ionospheric model, in the second step one is able to isolate the VTEC, the SDCBs and the RDCBs from the ionospheric observables. In this work, we present a single-frequency (SF) approach, enabling the joint estimation of VTEC and SDCBs using low-cost receivers; this approach is also based on two steps and it differs from the DF approach only in the first step, where we turn to the precise point positioning technique to retrieve from the single-frequency GNSS data the ionospheric observables, interpreted as the combination of the STEC, the SDCBs and the biased receiver clocks at the pivot epoch. Our numerical analyses clarify how SF approach performs when being applied to GPS L1 data collected by a single receiver under both calm and disturbed ionospheric conditions. The daily time series of zenith VTEC estimates has an accuracy ranging from a few tenths of a TEC unit (TECU) to approximately 2 TECU. For 73–96% of GPS satellites in view, the daily estimates of SDCBs do not deviate, in absolute value, more than 1 ns from their ground truth values published by the Centre for Orbit Determination in Europe.     

Anomalies in broadcast ionospheric coefficients recorded by GPS receivers over the past two solar cycles (1992–2013)

The anomaly phenomenon of broadcast ionospheric model coefficients of the Global Positioning System (GPS) is revealed after analyzing the navigation file data collected from all the IGS (International GNSS Service) stations worldwide over a 22-year period (1992–2013). GPS broadcast ionospheric coefficients widely used by many single-frequency users to correct the ionosphere errors for numerous GPS applications are usually believed to have only one set/version per day. However, it is found that GPS receivers from the IGS network can report as many as eight sets/versions of ionospheric coefficients in a day. In order to investigate the possible factors for such an anomalous phenomenon, the relationship between the number of coefficient sets and solar cycle, the receiver geographic locations, and receiver types/models are analyzed in detail. The results indicate that most of the coefficients show an annual variation. During the active solar cycle period from mid-1999 to mid-2001, all of the coefficients extracted from IGS navigation files behaved anomalously. Our analysis shows that the anomaly is also associated with GPS receiver types/models. Some types/models of GPS receivers report one set/version of ionospheric coefficients daily, while others report multiple sets. Our analysis also suggests that the ionospheric coefficient anomaly is not necessarily related to ionospheric scintillations. No correlation between the anomaly and geographic location of GPS receivers has been found in the analysis. Using the ionospheric coefficient data collected from 1998 to 2013, the impact of ionospheric coefficient anomaly on vertical total electron content (VTEC) calculation using the Klobuchar model has been evaluated with respect to the Global Ionospheric Maps generated by the Center for Orbit Determination in Europe. With different sets of coefficients recorded on the same day, the resulting VTEC values are dramatically different. For instance on June 1, 2000, the largest VTEC at one of our test stations can be as large as 153.3 TECu (total electron content unit) using one set of coefficients, which is 16.36 times larger than the smallest VTEC of 9.37 TECu computed from using another set of coefficients.     

Ionospheric error contribution to GNSS single-frequency navigation at the 2014 solar maximum

For single-frequency users of the global satellite navigation system (GNSS), one of the main error contributors is the ionospheric delay, which impacts the received signals. As is well-known, GPS and Galileo transmit global models to correct the ionospheric delay, while the international GNSS service (IGS) computes precise post-process global ionospheric maps (GIM) that are considered reference ionospheres. Moreover, accurate ionospheric maps have been recently introduced, which allow for the fast convergence of the real-time precise point position (PPP) globally. Therefore, testing of the ionospheric models is a key issue for code-based single-frequency users, which constitute the main user segment. Therefore, the testing proposed in this paper is straightforward and uses the PPP modeling applied to single- and dual-frequency code observations worldwide for 2014. The usage of PPP modeling allows us to quantify—for dual-frequency users—the degradation of the navigation solutions caused by noise and multipath with respect to the different ionospheric modeling solutions, and allows us, in turn, to obtain an independent assessment of the ionospheric models. Compared to the dual-frequency solutions, the GPS and Galileo ionospheric models present worse global performance, with horizontal root mean square (RMS) differences of 1.04 and 0.49 m and vertical RMS differences of 0.83 and 0.40 m, respectively. While very precise global ionospheric models can improve the dual-frequency solution globally, resulting in a horizontal RMS difference of 0.60 m and a vertical RMS difference of 0.74 m, they exhibit a strong dependence on the geographical location and ionospheric activity.     

Calibration errors on experimental slant total electron content (TEC) determined with GPS

The Global Positioning System (GPS) has become a powerful tool for ionospheric studies. In addition, ionospheric corrections are necessary for the augmentation systems required for Global Navigation Satellite Systems (GNSS) use. Dual-frequency carrier-phase and code-delay GPS observations are combined to obtain ionospheric observables related to the slant total electron content (sTEC) along the satellite-receiver line-of-sight (LoS). This observable is affected by inter-frequency biases [IFB; often called differential code biases (DCB)] due to the transmitting and the receiving hardware. These biases must be estimated and eliminated from the data in order to calibrate the experimental sTEC obtained from GPS observations. Based on the analysis of single differences of the ionospheric observations obtained from pairs of co-located dual-frequency GPS receivers, this research addresses two major issues: (1) assessing the errors translated from the code-delay to the carrier-phase ionospheric observable by the so-called levelling process, applied to reduce carrier-phase ambiguities from the data; and (2) assessing the short-term stability of receiver IFB. The conclusions achieved are: (1) the levelled carrier-phase ionospheric observable is affected by a systematic error, produced by code-delay multi-path through the levelling procedure; and (2) receiver IFB may experience significant changes during 1 day. The magnitude of both effects depends on the receiver/antenna configuration. Levelling errors found in this research vary from 1.4 total electron content units (TECU) to 5.3 TECU. In addition, intra-day vaiations of code-delay receiver IFB ranging from 1.4 to 8.8 TECU were detected.     

Real-time kinematic positioning of LEO satellites using a single-frequency GPS receiver

Due to their low cost and low power consumption, single-frequency GPS receivers are considered suitable for low-cost space applications such as small satellite missions. Recently, requirements have emerged for real-time accurate orbit determination at sub-meter level in order to carry out onboard geocoding of high-resolution imagery, open-loop operation of altimeters and radio occultation. This study proposes an improved real-time kinematic positioning method for LEO satellites using single-frequency receivers. The C/A code and L1 phase are combined to eliminate ionospheric effects. The epoch-differenced carrier phase measurements are utilized to acquire receiver position changes which are further used to smooth the absolute positions. A kinematic Kalman filter is developed to implement kinematic orbit determination. Actual flight data from China’s small satellite SJ-9A are used to test the navigation performance. Results show that the proposed method outperforms traditional kinematic positioning method in terms of accuracy. A 3D position accuracy of 0.72 and 0.79 m has been achieved using the predicted portion of IGS ultra-rapid products and broadcast ephemerides, respectively.     

Ionospheric correction using NTCM driven by GPS Klobuchar coefficients for GNSS applications

Global Navigation Satellite Systems (GNSS) require mitigation of ionospheric propagation errors because the ionospheric range errors might be larger than tens of meters at the zenith direction. Taking advantage of the frequency-dispersive property of ionospheric refractivity, the ionospheric range errors can be mitigated in dual-frequency applications to a great extent by a linear combination of carrier phases or pseudoranges. However, single-frequency GNSS operations require additional ionospheric information to apply signal delay or range error corrections. To aid single-frequency operations, the global positioning system (GPS) broadcasts 8 coefficients as part of the navigation message to drive the ionospheric correction algorithm (ICA) also known as Klobuchar model. We presented here an ionospheric correction algorithm called Neustrelitz TEC model (NTCM) which can be used as complementary to the GPS ICA. Our investigation shows that the NTCM can be driven by Klobuchar model parameters to achieve a significantly better performance than obtained by the mother ICA algorithm. Our research, using post-processed reference total electron content (TEC) data from more than one solar cycle, shows that on average the RMS modeled TEC errors are up to 40% less for the proposed NTCM model compared to the Klobuchar model during high solar activity period, and about 10% less during low solar activity period. Such an approach does not require major technology changes for GPS users rather requires only introducing the NTCM approach a complement to the existing ICA algorithm while maintaining the simplicity of ionospheric range error mitigation with an improved model performance.     

RINEX_HO: second- and third-order ionospheric corrections for RINEX observation files

When GNSS receivers capable of collecting dual-frequency data are available, it is possible to eliminate the first-order ionospheric effect in the data processing through the ionosphere-free linear combination. However, the second- and third-order ionospheric effects still remain. The first-, second- and third-order ionospheric effects are directly proportional to the total electron content (TEC), although the second- and third-order effects are influenced, respectively, by the geomagnetic field and the maximum electron density. In recent years, the international scientific community has given more attention to these kinds of effects and some works have shown that for high precision GNSS positioning these effects have to be taken into consideration. We present a software tool called RINEX_HO that was developed to correct GPS observables for second- and third-order ionosphere effects. RINEX_HO requires as input a RINEX observation file, then computes the second- and third-order ionospheric effects, and applies the corrections to the original GPS observables, creating a corrected RINEX file. The mathematical models implemented to compute these effects are presented, as well as the transformations involving the earth’s magnetic field. The use of TEC from global ionospheric maps and TEC calculated from raw pseudorange measurements or pseudoranges smoothed by phase is also investigated.     

Performance evaluation of different ionospheric models in single-frequency code-based differential GPS positioning

Differential ionospheric slant delays are obtained from a quiet-time, three-dimensional ionospheric electron density model, called the TaiWan Ionosphere Model (TWIM), to be used in code-based differential GPS positioning. The code observations are acquired from nine continuously operating GPS stations around Taiwan whose baseline ranged from 19 to 340 km. Daily 24-hour epoch-per-epoch positioning obtained for 70 most geomagnetic quiet days (2008–2010) for each of the 72 baselines. The performance of TWIM has been compared with the standard operational Klobuchar model (KLB) used by typical single-frequency receivers and the IGS global ionospheric model (GIM). Generally, TWIM performed well in reducing the differential ionospheric delay especially for long baselines and different levels of low solar activity. It has a much better performance compared to the operational KLB model. TWIM also performed similarly with GIM, though GIM has the best performance overall. GIM has the best ionospheric gradient estimates among the three models whose differential ionospheric delay-to-horizontal error ratio is more than 0.25. This is followed closely by TWIM with about 0.20. KLB only has a ratio of <0.10. The similarity of the performance of TWIM and GIM demonstrates the feasibility of TWIM in correcting for differential ionospheric delays in the C/A code pseudorange that is caused by electron density gradients in the ionosphere. It can provide decimeter-to-centimeter level accuracy in differential GPS positioning for single-frequency receivers during geomagnetic quiet conditions across all seasons and different levels of low solar activities.     

Performance of the L2C civil GPS signal under various ionospheric scintillation effects

As GPS is modernizing, there are currently fourteen satellites transmitting L2C civil code and seven satellites transmitting L5 signal. While the GPS observables are subject to several sources of errors, the ionosphere is one of the largest error sources affecting GPS signals. Small irregularities in the electrons density along the GPS radio signal propagation path cause ionospheric scintillation that is characterized by rapid fluctuations in the signal amplitude and phase. The ionospheric scintillation effects are stronger in equatorial and high-latitude geomagnetic latitude regions and occur mainly due to equatorial anomaly and solar storms. Several researchers have analyzed the L2C signal quality since becoming available in December, 2005. We analyze the performance of L2C using GPS data from stations in the equatorial region of Brazil, which is subject of weak, moderate and strong ionospheric scintillation conditions. The GPS data were collected by Septentrio PolaRxS–PRO receivers as part of the CIGALA/CALIBRA network. The analysis was performed as a function of scintillations indexes S4 and Phi60, lock time (time interval in seconds that the carrier phase is tracked continuously without cycle slips), multipath RMS and position variation of precise point positioning solutions. The analysis shows that L2C code solutions are less affected by multipath effects than that of P2 when data are collected under weak ionospheric scintillation effects. In terms of analysis of positions, the kinematic PPP results using L2C instead P2 codes show accuracy improvements up to 33 % in periods of weak or strong ionospheric scintillation. When combining phase and code collected under weak scintillation effects, the results by applying L2C against P2 provide improvement in accuracy up to 59 %. However, for data under strong scintillation effects, the use of L2C for PPP with code and phase does not provide improvements in the positioning accuracy.     

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.     

Increasing GNSS RTK availability with a new single-epoch batch partial ambiguity resolution algorithm

GPS Single-epoch Real-Time Kinematic positioning is immune to cycle slips and can be immediately re-initialized after loss-of-lock, providing high availability. This technique requires reliable ambiguity resolution: incorrect ambiguities can cause position errors of several meters, and failed ambiguity resolution reduces availability. However, a bias or inaccuracy in a single phase observation can prevent successful resolution of the whole set of ambiguities. Partial ambiguity resolution allows a subset of ambiguities to be resolved with greater probability of success than the full set. A new algorithm for resolving a subset of ambiguities with validation from previous epochs is described. If normal ambiguity resolution fails, all ambiguity subsets are generated and ordered with the best subsets first. Each subset is then resolved in turn. Fixed subsets are validated against values from previous epochs; this validation procedure greatly reduces the proportion of epochs with incorrect ambiguities. An additional algorithm is described that uses the fixed ambiguities as precise ranges to resolve the remaining unfixed ambiguities. In order to test these new algorithms, GPS data were collected from static and ship-based GPS receivers around Harwich harbor and processed from reference stations at distances up to 111 km. In the static tests the distance over which a 90% ambiguity resolution success rate for dual-frequency data was achieved was increased from 15 to 76 km. However, in some cases, the processing time was too long for this algorithm to be practical without a time-based cut-off. There is also a risk of incorrect ambiguities being propagated, particularly for single-frequency processing. In a ship-based test, the distance over which sufficient availability to support harbor navigation was achieved using single-epoch dual-frequency RTK was increased from 1 to 66 km.     

Ionospheric delay corrections for single-frequency GPS receivers over Europe using tomographic mapping

The majority of navigation satellites receivers operate on a single frequency and experience a positioning error due to the ionospheric delay. This can be compensated for using a variety of approaches that are compared in this paper. The study focuses on the last solar maximum. A 4D tomographic imaging technique is used to map the ionospheric electron density over the European region during 2002 and 2003. The electron density maps are then used to calculate the excess propagation delay on the L1 frequency experienced by GPS receivers at selected locations across Europe. The excess delay is applied to correct the pseudo-range single frequency observations at each location and the improvements to the resulting positioning are calculated. The real-time tomographic technique is shown to give navigation solutions that are better than empirical modelling methods and approach the accuracy of the full dual-frequency solution. The improvements in positioning accuracy vary from day to day depending on ionospheric conditions but can be up to 25 m during mid-day during these solar maximum conditions at European mid-latitudes.

Performance comparison of semicodeless GPS receivers for LEO satellites

This report provides a detailed performance analysis of three semicodeless dual-frequency GPS receivers for use in low Earth orbit (LEO). The test set comprises the IGOR receiver, which represents a follow-on of the flight-proven BlackJack receiver, as well as two geodetic receivers (NovAtel OEM4-G2 and Septentrio PolaRx2), which are entirely based on commercial-off-the-shelf technology (COTS). All three receivers are considered for upcoming flight projects or experiments and have undergone at least a preliminary environmental qualification program. Using extensive signal simulator tests, the cold start signal acquisition, tracking sensitivity, differential code biases, raw measurement accuracy, and navigation accuracy of each receiver have been assessed. All tests are based on a common scenario that is representative of an actual space mission and provides a realistic simulation of the signal dynamics and quality on a scientific LEO satellite. Compared to the other receivers, the IGOR instrument exhibits a superior tracking sensitivity and is thus best suited for occultation measurements with low tangent point altitudes. The OEM4-G2 and PolaRx2 receivers are likewise shown to properly track dual-frequency GPS signals and normal signal levels and to provide accurate code and carrier phase measurements. Given their limited resource requirements, these receivers appear well suited for precise orbit determination applications and ionospheric sounding onboard of microsatellites with tight mission budgets.     

Evaluating pseudorange multipath effects at stations in the National CORS Network

A preliminary study was conducted to evaluate the amount of pseudorange multipath at 390+ sites in the National Continuously Operating Reference Station (CORS) Network. The National CORS Network is a cooperative effort involving over 110 different agencies, universities, and private companies who seek to make GPS data from dual-frequency receivers located throughout the United States and its territories available to the general public. For CORS users, pseudorange multipath can seriously degrade the accuracy of any application that relies on precise measurements of the pseudorange observable over a short period of time, including differential pseudorange navigation, kinematic and rapid-static surveying, and ionospheric monitoring. The main objectives of this study were to identify the most affected and least affected sites in the network, to closely investigate problematic sites, and to compare various receiver/antenna combinations. Dual-frequency carrier phase and pseudorange measurements were used to estimate the amount of L1 and L2 pseudorange multipath at each site over a one-year period. Some of the most severely affected sites were maritime Differential GPS and Nationwide Differential GPS (DGPS/NDGPS) sites. Photographs obtained for these sites verified the presence of transmission towers and other reflectors in close proximity to the GPS antennas. Plotting the variations of the L1 and L2 pseudorange multipath with respect to azimuth and elevation further verified that even above a 60° elevation angle there was still as much as five meters of pseudorange multipath at some sites. The least affected sites were the state networks installed in Ohio and Michigan; these sites used excellent antenna mounts, choke ring antennas, and new receiver technology. A comparison of the 12 most commonly used receiver/antenna combinations in the CORS Network indicated that newer receivers such as the Ashtech UZ-12, Leica RS-500, and Trimble 5700 help to significantly mitigate pseudorange multipath, while the receivers/antennas at some DGPS/NDGPS sites, and the antennas formerly used at the Wide Area Augmentation System (WAAS) sites, are among those most affected by pseudorange multipath. The receiver/antenna comparison did not take into account the potential presence of reflectors at the sites (i.e., it is possible that a well-performing receiver/antenna combination could have been consistently placed at very poor site locations, and vice-versa).Product Disclaimer: Mention of a commercial company or product does not constitute an endorsement by the National Oceanic and Atmospheric Administration. Use for publicity or advertisement purposes of information from this paper concerning proprietary products or the comparison of such products is not authorized.     

A double-differenced cycle slip detection and repair method for GNSS CORS network

The difficulty to detect and repair cycle slip of carrier phase measurements is a key limit for continuously high accuracy of GNSS positioning and navigation services. We propose an automated cycle slip detection and repair method for data preprocessing of a CORS network. The method jointly uses double-differenced (DD) geometry-free (GF) combination and ionospheric-free observation corrected for the computed geometrical distance (IF-OMC) to estimate the cycle slips in dual-frequency observations. The DD GF combination, which is only affected by the ionospheric residual, can be used to detect cycle slips with high reliability except for special pairs such as (77, 60) on GPS L1/L2 frequencies. The detection principle of the IF-OMC observable is such that there is a large discontinuity related to the previous epoch when cycle slips occur at the present epoch. The disadvantages of these two combinations can be overcome employing the proposed detection method. The cycle slip pair (77, 60) has no effect on the GF combination, while a change of 14.65 m is derived from GPS L1/L2 observations using the IF-OMC algorithm. Using pre-determined station coordinates as precise values, we found that the accuracy of the DD IF-OMC combination was 18 mm for a 200-km CORS baseline. Therefore, cycle slips in dual-frequency observations can be correctly and uniquely determined using DD GF and IF-OMC equations. The proposed method was verified by adding simulated cycle slips in observations collected from the CORS network under a quiet ionosphere and shown to be effective. Moreover, the method was assessed with observations made during intense ionospheric activity, which generated extensive cycle slips. The results show that the algorithm can detect and repair all cycle slips apart from two exceptions relating to long data gaps.