Precise ionosphere-free single-frequency GNSS positioning

Ionospheric delays can be efficiently eliminated from single-frequency data using a combination of carrier phases and code ranges. Unfortunately, GPS and GLONASS ranges are relatively noisy which can limit the use of the positioning method. Nevertheless, position standard deviations are in the range of 6–8 cm (horizontal) and 7–9 cm (3d) obtained from diurnal data batches from selected IGS reference stations can be further reduced to 2–3 cm (3d) for weekly smoothed averages. GPS data sets collected in Ghana (Africa) reveal a typical level of 10 cm of deviation that must be anticipated under average conditions. Looking at the future of GNSS, the European Galileo system will, in contrast to GPS, provide the broadband signal E5 that is by far less affected by multipath thus providing rather precise range measurements. Simulated processing runs featuring both high ionospheric and tropospheric delay variations show a 3d position precision of 4 cm even for a data batch as short as just 1 h, whereas GPS L1/Galileo E1 performance is close to 13 cm for the same data set.     

Galileo系统的特点分析

随着“全球定位系统(GPS)”在各个领域应用的不断深入，欧盟也认识到了拥有自身卫星导航系统的重要性，经过长时间的协调和争论，终于在2002年3月26日由欧盟各国交通部长正式签署协议开始建设欧洲自己的卫星无线电导航系统-Galileo系统。我国在Galileo系统谈判阶段就开始与欧盟接触，寻求未来在Galileo系统上与欧洲相关的技术协作和系统共享，目前已经达成初步协议，并于2003年9月18日在北京成立了“中欧卫星导航系统培训与合作中心”。对该系统的总体设计思路、系统特点以及相关技术进行了分析。     

顾及卫星姿态的多系统精密钟差产品综合

国际GNSS服务(IGS)提供的GPS综合产品被广泛应用于各种高精度科学研究中. 随着各国卫星导航系统的发展,亟需研究针对多系统全球卫星导航系统(GNSS)产品的综合策略. 由于卫星姿态与钟差相互耦合,综合钟差时额外考虑姿态改正将进一步提高综合产品精度,因此研究了一种顾及卫星姿态的GNSS钟差综合策略,改正姿态后GPS综合残差最大可减小80%. 对142个IGS测站进行精密单点定位(PPP)解算发现,综合产品比单个分析中心产品更加稳定,东(E)、北(N)、高(U)方向的动态定位精度最大可提升22.7%、16.7%和18.3%. 相对于未顾及姿态改正的综合产品,顾及姿态改正的综合产品的动态定位精度最大可提升65.3%.      

Single epoch estimation of the Galileo integrity chain sensor station clock offsets

The Galileo integrity chain depends on a number of key factors, one of which is contamination of the signal-in-space errors with residual errors other than imperfect modelling of satellite orbits and clocks. A potential consequence of this is that the user protection limit is driven not by the errors associated with the imperfect orbit and clock modelling, but by the distortions induced by noise and bias in the integrity chain. These distortions increase the minimum bias the integrity chain can guarantee to detect, which is reflected in the user protection limit. A contributor to this distortion is the inaccuracy associated with the estimation of the offset between the Galileo sensor station (GSS) receiver clocks and the Galileo system time (GST). This offset is termed the receiver clock synchronization error (CSE). This paper describes the research carried out to determine both the CSE and its associated error using GPS data as captured with the Galileo System Test Bed Version 1 (GSTB-V1). In the study we simulate open access to a time datum using IGS data. Two methods are compared for determining CSE and the corresponding uncertainty (noise) across a global network of tracking stations. The single-epoch single-station method is an ‘averaging’ technique that uses a single epoch of data, and is carried out at individual sensor stations, without recourse to the data from other stations. The global network solution method is also single epoch based, but uses the inversion of a linearised model of the global system to solve for the CSE simultaneously at all GSS along with a number of other parameters that would otherwise be absorbed into the CSE estimate in the averaging technique. To test the effectiveness of various configurations in the two methods the estimated synchronisation errors across the GSS network (comprising 25 stations) are compared to the same values as estimated by the International GPS Service (IGS) using a global tracking network of around 150 stations, as well as precise orbit and satellite clock models determined by a combination of global analysis centres. The results show that the averaging technique is vulnerable to unmodelled errors in the satellite clock offsets from system time, leading to receiver CSE errors in the region of 12 ns (3.7 m), this value being largely driven by the satellite CSE errors. The global network approach is capable of delivering CSE errors at the level of 1.5 ns (46 cm) depending on the number of parameters in the linearised model. The International GNSS Service (IGS) receiver clock estimates were used as a truth model for comparative assessment.     

New IGS Station and Satellite Clock Combination

Following the principles set forth in the Position Paper #3 at the 1998 Darmstadt Analysis Center (AC) Workshop on the new International GPS Service (IGS) International Terrestrial Reference Frame (ITRF) realization and discussions at the 1999 La Jolla AC workshop, a new clock combination program was developed. The program allows for the input of both SP3 and the new clock (RINEX) format (ftp://igsch.jpl.nasa.gov//igscb/data/format/rinex_clock.txt). The main motivation for this new development is the realization of the goals of the IGS/BIPM timing project. Besides this there is a genuine interest in station clocks and a need for a higher sampling rate of the IGS clocks (currently limited to 15 min due to the SP3 format). The inclusion of station clocks should also allow for a better alignment of the individual AC solutions and should enable the realization of a stable GPS time-scale. For each input AC clock solution the new clock combination solves and corrects for reference clock errors/instabilities as well as satellite/station biases, geocenter and station/satellite orbit errors. External station clock calibrations and/or constraints, such as those resulting from the IGS/BIPM timing pilot project, can be introduced via a subset of the fiducial timing station set, to facilitate a precise and consistent IGS UTC realization for both station and satellite combined clock solutions. Furthermore, the new clock combination process enforces strict strict conformity and consistency with the current and future IGS standards. The new clock combination maintains orbit/clock consistency at millimeter level, which is comparable to the best AC orbit/clock solutions. This is demonstrated by static GIPSY precise point positioning tests using GPS week 0995 data for stations in both Northern and Southern Hemispheres and similar tests with the Bernese software using more recent data from GPS week 1081. ? 2001 John Wiley & Sons, Inc.     

The international global navigation satellite systems service (IGS): development and achievements

Since 21 June 1992 the International GPS Service (IGS), renamed International GNSS Service in 2005, produces and makes available uninterrupted time series of its products, in particular GPS observations from the IGS Global Network, GPS orbits, Earth orientation parameters (components x and y of polar motion, length of day) with daily time resolution, satellite and receiver clock information for each day with different latencies and accuracies, and station coordinates and velocities in weekly batches for further analysis by the IERS (International Earth Rotation and Reference Systems Service). At a later stage the IGS started exploiting its network for atmosphere monitoring, in particular for ionosphere mapping, for troposphere monitoring, and time and frequency transfer. This is why new IGS products encompass ionosphere maps and tropospheric zenith delays. This development became even more important when more and more space-missions carrying space-borne GPS for various purposes were launched. This article offers an overview for the broader scientific community of the development of the IGS and of the spectrum of topics addressed today with IGS data and products.     

GPS and GLONASS Integration: Modeling and Ambiguity Resolution Issues

The integration of GPS with GLONASS may be considered a major milestone in satellite-based positioning, because it can dramatically improve the reliability and productivity of said positioning. However, unlike GPS, GLONASS satellites transmit signals at different frequencies, which result in significant complexity in terms of modeling and ambiguity resolution for integrated GPS and GLONASS positioning systems. In this paper, a variety of mathematical and stochastic modeling methodologies and ambiguity resolution strategies are analyzed, and some remaining research challenges are identified. The exercise, of developing mathematical models and processing methodologies for integrated systems based on more than one satellite system, is a valuable one as it identified crucial issues concerned with the combination of any two or more microwave positioning systems, be they satellite-based or terrestrial. Hence these are experiences that can be applied to future projects that might integrate GPS with Galileo, or GLONASS and Galileo, or all three. ? 2001 John Wiley & Sons, Inc.     

The activities of the Ionosphere Working Group of the International GPS Service (IGS)

This article is based on a position paper presented at the IGS Network, Data and Analysis Center Workshop 2002 in Ottawa, Canada, 8–11 April 2002, and introduces the IGS Ionosphere Working Group (Iono_WG). Detailed information about the IGS in general can be found on the IGS Central Bureau Web page: http://igscb.jpl.nasa.gov. The Iono_WG commenced working in June 1998. The working group's main activity currently is the routine production of ionosphere Total Electron Content (TEC) maps with a 2-h time resolution and daily sets of GPS satellite and receiver hardware differential code bias (DCB) values. The TEC maps and DCB sets are derived from GPS dual-frequency tracking data recorded with the global IGS tracking network. In the medium- and long-term, the working group intends to refine algorithms for the mapping of ionospheric parameters from GPS measurements and to realize near–real–time availability of IGS ionosphere products. The paper will give an overview of the Iono_WG activities that include a summary of activities since its establishment, achievements and future plans. Electronic Publication     

基于STK的全球卫星导航定位系统DOP值仿真

通过STK卫星数据库建立了GPS卫星星座;根据GALILEO卫星设计参数,利用轨道向导建立GALILEO卫星星座。比较了武汉地区GPS、GALILEO、GPS和GALILEO组合的卫星可见数和GDOP值,以及全球GDOP值的分布。     

A simplified yaw-attitude model for eclipsing GPS satellites

A simplified yaw-attitude modeling, consistent with Bar-Sever (1996), has been implemented and tested in the NRCan PPP software. For Block IIR GPS satellite it is possible to model yaw-attitude control during eclipsing periods by using the constant hardware yaw rate of 0.20°/s. The Block IIR satellites maintain the nominal yaw attitude even during a shadow crossing (Y. E. Bar-Sever, private communication, 2007), except for the noon and shadow midnight turn maneuvers, both of which can be modeled and last up to 15 min. Thus, for Block IIR satellites it is possible to maintain continuous satellite clock estimation even during eclipsing periods. For the Block II/IIA satellites, it is possible to model satisfactorily the noon turns and also shadow crossing, thanks to the permanent positive yaw bias of 0.5°, implemented in November 1995. However, in order to model the Block II/IIA shadow crossings, satellite specific yaw rates should be used, either solved for or averaged yaw-rate solutions. These yaw rates as estimated by the Jet Propulsion Laboratory (JPL) can differ significantly from the nominal hardware values. The Block II/IIA post-shadow recovery periods, which last about 30 min, should be considered uncertain and cannot be properly modeled. Data from post-shadow recovery periods should, therefore, not be used in precise global GPS analyses (Bar-Sever 1996). For high-precision applications, it is essential that users implement a yaw-attitude model, which is consistent with the generation of the satellite clocks. Initial testing and analyses, based on the IGS and AC Final orbits and clocks have revealed that during eclipsing periods, significant inconsistencies in yaw-attitude modeling still exist amongst the IGS Analyses Centers, which contribute to the errors of the IGS Final clock combinations.     

Impact of Earth radiation pressure on GPS position estimates

GPS satellite orbits available from the International GNSS Service (IGS) show a consistent radial bias of up to several cm and a particular pattern in the Satellite Laser Ranging (SLR) residuals, which are suggested to be related to radiation pressure mismodeling. In addition, orbit-related frequencies were identified in geodetic time series such as apparent geocenter motion and station displacements derived from GPS tracking data. A potential solution to these discrepancies is the inclusion of Earth radiation pressure (visible and infrared) modeling in the orbit determination process. This is currently not yet considered by all analysis centers contributing to the IGS final orbits. The acceleration, accounting for Earth radiation and satellite models, is introduced in this paper in the computation of a global GPS network (around 200 IGS sites) adopting the analysis strategies from the Center for Orbit Determination in Europe (CODE). Two solutions covering 9 years (2000–2008) with and without Earth radiation pressure were computed and form the basis for this study. In previous studies, it has been shown that Earth radiation pressure has a non-negligible effect on the GPS orbits, mainly in the radial component. In this paper, the effect on the along-track and cross-track components is studied in more detail. Also in this paper, it is shown that Earth radiation pressure leads to a change in the estimates of GPS ground station positions, which is systematic over large regions of the Earth. This observed “deformation” of the Earth is towards North–South and with large scale patterns that repeat six times per GPS draconitic year (350 days), reaching a magnitude of up to 1 mm. The impact of Earth radiation pressure on the geocenter and length of day estimates was also investigated, but the effect is found to be less significant as compared to the orbits and position estimates.     

GBM快速轨道产品及非差模糊度固定对其精度的改进

GNSS是实时定位导航最重要的方法,精密卫星轨道钟差产品是GNSS高精度服务的前提。国际GNSS服务中心(IGS)及其分析中心长期致力于GNSS数据处理的研究及高精度轨道和钟差产品的提供。GFZ作为分析中心之一,提供GBM多系统快速产品。本文基于2015—2021年GBM提供的精密轨道产品,阐述了数据处理策略,分析了轨道的精度,介绍了非差模糊度固定的原理和对精密定轨的影响。结果表明:GBM快速产品中的GPS轨道精度与IGS后处理精密轨道相比的精度约为11~13 mm,轨道6 h预报精度约为6 cm;GLONASS预报精度约为12 cm,Galileo在该时期的精度均值为10 cm,但是在2016年底以后精度提升到5 cm左右;北斗系统的中轨卫星(medium earth orbit,MEO)在2020年以后预报精度约为10 cm;北斗的静止轨道卫星(geostationary earth orbit,GEO)卫星和QZSS卫星的预报精度在米级;卫星激光测距检核表明,Galileo、GLONASS、BDS-3 MEO卫星轨道精度分别为23、41、47 mm;此外,采用150 d观测值的试验结果表明,采用非差模糊度固定能显著改善MEO卫星轨道精度,对GPS、GLONASS、Galileo、BDS-2和BDS-3的MEO卫星的6 h时预报精度改善率分别为9%~15%、15%~18%、11%~13%、6%~17%和14%~25%。     

GNSS receiver autonomous integrity monitoring (RAIM) performance analysis

The availability of GPS signals is a major concern for many existing and potential applications. Fortunately, with the development of Galileo by the European Commission (EC) and European Space Agency (ESA) and new funding for the restoration of the Russian GLONASS announced by the Russian Federation (Revnivykh et al., in European Navigation Conference 2005, Munich, Germany, July 19–22), the future for satellite-based positioning and navigation applications is extremely promising. With the complete cooperation from all these global navigation satellite systems (GNSS), greater levels of satellite visibility and therefore integrity can be expected. In this paper, a receiver autonomous integrity monitoring (RAIM) scheme along with reliability and separability measures are used to assess integrity performance levels of standalone GPS and integrated GPS/GLONASS, GPS/Galileo and GPS/GLONASS/Galileo systems where the clock offsets for each of the additional systems are estimated. It is shown, herein, that a minimum of three satellites must be visible in an additional system in order to provide a full integrity contribution when the system’s clock offset is to be estimated within the adjustment. A comparison of the integrity results obtained via system clock offsets estimated in the adjustment versus the case where the offsets are known and the measurements are corrected prior to the adjustment is also made for a high elevation mask scenario. Global simulation results for combined GPS/GLONASS/Galileo show that, theoretically, for the time of simulation and for any point on the globe, an outlier of 20 m can be detected with 80% probability at the 0.5% significance level and then separated from any other measurement with 90% probability. Corresponding values for the GPS only and combined GPS/GLONASS and GPS/Galileo systems, respectively, are approximately 435, 110 and 28 m, respectively, for the maximum MSBs and 312, 50 and 26 m, respectively, for the maximum MDBs. Temporal 24 h simulations for the GPS/GLONASS/Galileo scenario delivered agreeable results with the global snapshots for a 15° elevation mask. For the case where system clock offsets are estimated within the adjustment, it was shown that only the reliability measure was available for 100% of the time, with horizontal external reliability values of no more than about 12 m when a 30° masking angle was used. By assuming the clock offsets were determined and corrected for prior to the adjustment, the separability measure was markedly improved and was also available 100% of the time.     

精密卫星钟差加密方法及其对星载GPS低轨卫星定轨精度影响

系统分析、比较了几种精密卫星钟差加密方法，研究了利用全球分布的IGS永久跟踪站的GPS观测数据估计高采样率卫星钟差参数的原理与方法，并将各种卫星钟差加密方法得到的结果与IGS数据分析中心估计的卫星钟差结果相比较。最后将不同加密方法得出的精密卫星钟差结果用于基于星载GPS双频非差观测值的CHAMP低轨卫星的定轨，并将不同方法得到的定轨精度进行比较。结果表明，利用地面跟踪站的GPS观测数据，可高精度、高密度地估计GPS卫星钟差，估计精度可达0．1～0．5ns。经地面GPS跟踪站数据估计的GPS卫星钟差，应用于基于PPP方法的低轨卫星定轨，其定轨精度在10cm以内。     

Characterization of periodic variations in the GPS satellite clocks

The clock products of the International Global Navigation Satellite Systems (GNSS) Service (IGS) are used to characterize the timing performance of the GPS satellites. Using 5-min and 30-s observational samples and focusing only on the sub-daily regime, approximate power-law stochastic processes are found. The Block IIA Rb and Cs clocks obey predominantly random walk phase (or white frequency) noise processes. The Rb clocks are up to nearly an order of magnitude more stable and show a flicker phase noise component over intervals shorter than about 100 s. Due to the onboard Time Keeping System in the newer Block IIR and IIR-M satellites, their Rb clocks behave in a more complex way: as an apparent random walk phase process up to about 100 s and then changing to flicker phase up to a few thousand seconds. Superposed on this random background, periodic signals have been detected in all clock types at four harmonic frequencies, n × (2.0029 ± 0.0005) cycles per day (24 h coordinated universal time or UTC), for n = 1, 2, 3, and 4. The equivalent fundamental period is 11.9826 ± 0.0030 h, which surprisingly differs from the reported mean GPS orbital period of 11.9659 ± 0.0007 h by 60 ± 11 s. We cannot account for this apparent discrepancy but note that a clear relationship between the periodic signals and the orbital dynamics is evidenced for some satellites by modulations of the spectral amplitudes with eclipse season. All four harmonics are much smaller for the IIR and IIR-M satellites than for the older blocks. Awareness of the periodic variations can be used to improve the clock modeling, including for interpolation of tabulated IGS products for higher-rate GPS positioning and for predictions in real-time applications. This is especially true for high-accuracy uses, but could also benefit the standard GPS operational products. The observed stochastic properties of each satellite clock type are used to estimate the growth of interpolation and prediction errors with time interval.     

Anomalous harmonics in the spectra of GPS position estimates

Prior studies of the power spectra of GPS position time series have found pervasive seasonal signals against a power-law background of flicker noise plus white noise. Dong et al. (2002) estimated that less than half the observed GPS seasonal power can be explained by redistributions of geophysical fluid mass loads. Much of the residual variation is probably caused by unidentified GPS technique errors and analysis artifacts. Among possible mechanisms, Penna and Stewart (2003) have shown how unmodeled analysis errors at tidal frequencies (near 12- and 24-hour periods) can be aliased to longer periods very efficiently. Signals near fortnightly, semiannual, and annual periods are expected to be most seriously affected. We have examined spectra for the 167 sites of the International GNSS (Global Navigation Satellite Systems) Service (IGS) network having more than 200 weekly measurements during 1996.0–2006.0. The non-linear residuals of the weekly IGS solutions that were included in ITRF2005, the latest version of the International Terrestrial Reference Frame (ITRF), have been used. To improve the detection of common-mode signals, the normalized spectra of all sites have been stacked, then boxcar smoothed for each local north (N), east (E), and height (H) component. The stacked, smoothed spectra are very similar for all three components. Peaks are evident at harmonics of about 1 cycle per year (cpy) up to at least 6 cpy, but the peaks are not all at strictly 1.0 cpy intervals. Based on the 6th harmonic of the N spectrum, which is among the sharpest and largest, and assuming a linear overtone model, then a common fundamental of 1.040 ± 0.008 cpy can explain all peaks well, together with the expected annual and semiannual signals. A flicker noise power-law continuum describes the background spectrum down to periods of a few months, after which the residuals become whiter. Similar sub-seasonal tones are not apparent in the residuals of available satellite laser ranging (SLR) and very long baseline interferometry (VLBI) sites, which are both an order of magnitude less numerous and dominated by white noise. There is weak evidence for a few isolated peaks near 1 cpy harmonics in the spectra of geophysical loadings, but these are much noisier than for GPS positions. Alternative explanations related to the GPS technique are suggested by the close coincidence of the period of the 1.040 cpy frequency, about 351.2 days, to the “GPS year”; i.e., the interval required for the constellation to repeat its inertial orientation with respect to the sun. This could indicate that the harmonics are a type of systematic error related to the satellite orbits. Mechanisms could involve orbit modeling defects or aliasing of site-dependent positioning biases modulated by the varying satellite geometry.     

卫星VLBI测轨资料处理中的电离层延迟改正

结合我国探月项目卫星VLBI测轨资料分析中的实际需求讨论了两个问题：一是在S、X波段时延测量精度均为1ns情况下，电离层延迟改正所能够达到的精度；二是在飞行器VLBI测轨过程中，不能确保S、X波段双频观测情况下获取电离层时延改正的可能途径，包括借助于相关电离层模型、利用常规VLB1历史观测资料积累、借助于局域GPS观测网和IGS网单站GPS测量以及借助于专门设计的单站GPS测量等。最后对电离层VLB1和GPS技术实测结果进行了比较和问题分析。     

Ambiguity resolution for triple-frequency geometry-free and ionosphere-free combination tested with real data

The recent GPS Block IIF satellites SVN62 and SVN63 and the Galileo satellites GIOVE-A, GIOVE-B, PFM and FM2 already send signals on more than two frequencies, and more GNSS satellites will provide tracking data on at least three frequencies in the near future. In this paper, a simplified general method for ambiguity resolution minimizing the noise level for the triple-frequency geometry-free (GF) and ionosphere-free (IF) linear combinations is presented, where differently scaled code noise on the three frequencies was introduced. For the third of three required linear combinations, the most demanding one in triple-frequency ambiguity resolution, we developed a general method using the ambiguity-corrected phase observations without any constraints to search for the optimal GF and IF linear combination. We analytically demonstrate that the noise level of this third linear combination only depends on the three frequencies. The investigation concerning this frequency-dependent noise factor was performed for GPS, Galileo and Compass frequency triplets. We verified the theoretical derivations with real triple-frequency GPS and Galileo data from the Multi-GNSS Experiment (M–GEX) of the International GNSS Service (IGS). The data of about 30 M–GEX stations around the world over 11 days from 29 April 2012 to 9 May 2012 were used for the test. For the third linear combinaton using Galileo E1, E5b and E5a, which is expected to have the worst performance among all the GNSS frequency triplets in our investigation, the formal errors of the estimated ambiguities are in most cases below 0.2 cycles after 400 observation epochs. If more GPS satellites sending signals on three frequencies or more stations tracking Galileo E6 signal are available in the future, an improvement by a factor of two to three can be expected.     

Quality assessment of GPS reprocessed terrestrial reference frame

The International GNSS Service (IGS) contributes to the construction of the International Terrestrial Reference Frame (ITRF) by submitting time series of station positions and Earth Rotation Parameters (ERP). For the first time, its submission to the ITRF2008 construction is based on a combination of entirely reprocessed GPS solutions delivered by 11 Analysis Centers (ACs). We analyze the IGS submission and four of the individual AC contributions in terms of the GNSS frame origin and scale, station position repeatability and time series seasonal variations. We show here that the GPS Terrestrial Reference Frame (TRF) origin is consistent with Satellite laser Ranging (SLR) at the centimeter level with a drift lower than 1 mm/year. Although the scale drift compared to Very Long baseline Interferometry (VLBI) and SLR mean scale is smaller than 0.4 mm/year, we think that it would be premature to use that information in the ITRF scale definition due to its strong dependence on the GPS satellite and ground antenna phase center variations. The new position time series also show a better repeatability compared to past IGS combined products and their annual variations are shown to be more consistent with loading models. The comparison of GPS station positions and velocities to those of VLBI via local ties in co-located sites demonstrates that the IGS reprocessed solution submitted to the ITRF2008 is more reliable and precise than any of the past submissions. However, we show that some of the remaining inconsistencies between GPS and VLBI positioning may be caused by uncalibrated GNSS radomes.     

GPS、Galileo及其组合系统导航定位的DOP值分析

利用Galileo系统的轨道参数模拟了Galileo系统的卫星位置，并计算出DOP值，比较了同一地点同时段内GPS和Galileo的能见卫星个数及DOP值。与单个系统相比，GPS／Galileo组合在任意地点、任意时刻的观测卫星个数增加．其DOP值明显变小；同时计算了组合系统在不同卫星高度角情况下的DOP值变化，发现在卫星高度角较高时仍可以接收到四颗以上的卫星．其DOP值也能达到导航定位的要求。