GPS/BDS卫星姿态异常对PPP相位缠绕的影响及其改正模型

在精密单点定位中,相位缠绕是一项不可忽略的误差。相位缠绕的计算严格依赖于卫星姿态的确立,不同的卫星类型产生不同的异常。本文给出了卫星在正常情况下的姿态模型和在异常情况下的姿态改正模型。使用真实数据测试以验证本文所提出模型的正确性。观察滤波收敛后出现异常情况的卫星观测值的残差,结果表明:在异常时期残差最大可能超过20 cm,然而使用本文的改正模型,残差可降低到5 cm以下。使用不同分析中心的精密轨道和钟差产品,效果存在微小差异。II/IIA卫星通过地影区域的时间最长可达1 h,此期间卫星姿态完全受航向角偏差(II/IIA为+0.5°)控制,出了地影区域后30 min,姿态难以模型化,因此这30 min的观测数据不建议采用。     

Improved relativistic transformations in GPS

For GPS satellite clocks, a nominal (hardware) frequency offset and a conventional periodic relativistic correction derived as a dot product of the satellite position and velocity vectors, are used to compensate the relativistic effects. The conventional hardware clock rate offset of 38,575.008 ns/day corresponds to a nominal orbit semi-major axis of about 26,561,400 m. For some of the GPS satellites, the departures from the nominal semi-major axis can cause an apparent clock rate up to 10 ns/day. GPS orbit perturbations, together with the earth gravity field oblateness, which is largely responsible for the orbit perturbations, cause the standard GPS relativistic transformations to depart from the rigorous relativity transformation by up to 0.2 ns/day. In addition, the conventional periodic relativistic correction exhibits periodic errors with amplitudes of about 0.1 and 0.2 ns, with periods of about 6 h and 14 days, respectively. Using an analytical integration of the gravity oblateness term (J₂), a simple analytical approximation was derived for the apparent clock rate and the 6-h periodic errors of the standard GPS gravity correction. For daily linear representations of GPS satellite clocks, the improved relativistic formula was found to agree with the precise numerical integration of the GPS relativistic effects within about 0.015 ns. For most of the Block IIR satellites, the 6-h periodical errors of the GPS conventional relativistic correction are already detectable in the recent IGS final clock combinations.     

GPS卫星钟的特性与预报研究

经实验分析发现:GPS卫星钟差的预报精度与卫星的种类密切相关,最近发射的BLOCK ⅡR和BLOCKIIR-M类卫星比以往的BLOCK ⅡA类卫星要更加稳定,其卫星钟差的预报精度明显较高。直接利用IGS超快速产品和线性模型预报后6小时的卫星钟差,精度在0.5纳秒水平;但一些BLOCK ⅡA类卫星是不稳定的,通过对其预报残差的分析发现:同一颗卫星每天在相同时段用相同的模型去预报其卫星钟差,预报所得的残差呈周期性变化,并且这种周期性变化并不完全重合,还具有一定的随机性。依据这一特性本文构建了一个新的预报模型来实时预报GPS卫星钟差。该模型不仅能预报卫星钟差的总体变化趋势,还能预报残差的周期性变化以及随机项的变化,因此精度更高。预报结果均与IGS发布的最终产品相比,实验显示利用该方法实时预报GPS卫星钟差,预报精度可达0.5纳秒水平。     

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.     

An investigation of GNSS atomic clock behavior at short time intervals

A technique for obtaining clock measurements from individual GNSS satellites at short time intervals is presented. The methodology developed in this study allows for accurate satellite clock stability analysis without an ultra-stable clock at the ground receiver. Variations in the carrier phase caused by the satellite clock are isolated using a combination of common GNSS carrier-phase processing techniques. Furthermore, the white phase variations caused by the thermal noise of the collection and processing equipment are statistically modeled and removed, allowing for analysis of clock performance at subsecond intervals. Allan deviation analyses of signals collected from GPS and GLONASS satellites reveal distinct intervals of clock noise for timescales less than 100 s. The clock data collected from GPS Block IIA, IIR, IIR-M, and GLONASS satellites reveal similar stability performance at time periods greater than 20 s. The GLONASS clock stability in the 0.6–10 s range, however, is significantly worse than GPS. Applications that rely on ultra-stable clock behavior from the GLONASS satellites at these timescales may therefore require high-rate corrections to estimate and remove oscillator-based errors in the carrier phase.     

GNSS星载原子钟性能评估

星载原子钟作为导航卫星上维持时间尺度的关键载荷,其性能会对用户进行导航、定位与授时的精度带来影响。介绍了原子钟评估常用的三个指标（频率准确度、飘移率和稳定度）的定义及计算方法,利用事后卫星精密钟差数据,开展了全球卫星导航系统（global navigation satellite system,GNSS）星载原子钟性能评估,分析了GNSS星载原子钟特性。结果表明,GPS（global position system）BLOCKIIF星载铷钟与Galileo星载氢钟综合性能最优;北斗系统中地球轨道卫星与倾斜同步轨道卫星星载原子钟天稳定度达到2~4×10-14量级,与BLOCK IIR卫星精度相当;频率准确度达到1~4×10-11量级;频率漂移率达到10-14量级。     

Relativistic Time Transformations in GPS

Since Selective Availability was permanently switched off on 7 May 2000, most of the GPS satellite clocks have been well behaved. During a 24-h period precise satellite clock solutions, corrected for GPS conventional relativistic corrections, follow straight lines within a few nanoseconds. The linear clock fit RMS for the best satellite clocks are well below the 1-ns level, which is consistent with the nominal stability of the GPS frequency standards. Typically, the GPS satellite clocks show an Allan variance at or below one part in 10¹¹/100 s for the Cesium frequency standards and a few parts in 10¹²/100 s for the Rubidium frequency standards. These results correspond to clock RMSs for 15-min sampling at or below 3 and 0.3 ns, respectively. This already confirms experimentally that the conventional periodic relativity correction of the GPS system, also adopted for all the IGS clock solution products, is precise and correct to 0.6 ns or better. To establish the precision limits of the GPS conventional relativity treatment, the relativistic time transformations of GPS satellite frequency and clocks are critically reviewed, taking into account all the contributions larger than the 10⁻¹⁸ (or 0.001 ns). The conventional GPS relativity treatment was found to be accurate, i. e., correctly modeling the actual relativistic frequency (clock rate) effects of GPS satellites at about the 10⁻¹⁴ level. However, it is also affected by small periodic errors of the same magnitude. The integration of these small periodic frequency relativistic errors gives the approximation errors of the conventional periodic relativistic clock correction with amplitudes of about 0.1 ns and a predominant period equal to a half of the orbital period (∼ 6 h). These approximation errors of the conventional GPS relativistic clock correction are at about the same level as the current precision of the IGS clock solutions. ? 2002 Wiley Periodicals, Inc.     

偏航姿态对GPS和GLONASS精密轨道和钟差的影响

利用全球约110个国际GNSS服务（International GNSS Service,IGS）测站2013年全年观测数据,分析和研究了GPS和全球卫星导航系统（global navigation satellite system,GLONASS）卫星偏航姿态对其精密轨道和钟差的影响。结果表明,偏航姿态对不同型号GPS卫星轨道和钟差的影响程度不同,当采用偏航姿态改正后地影期的BLOCK ⅡA型卫星轨道改善可达17 mm,BLOCK ⅡF为近5 mm,而BLOCK ⅡR几乎不受影响。由于偏航姿态对GLONASS-M卫星定轨精度影响较大,因此,当改正偏航姿态后所有GLONASS卫星相对于IGS最终轨道平均一维差异提高10 mm,相对于德国地学中心（German Research Center for Geosciences,GFZ）最终钟差平均标准差提升0.034 ns。     

Apparent clock variations of the Block IIF-1 (SVN62) GPS satellite

The Block IIF satellites feature a new generation of high-quality rubidium clocks for time and frequency keeping and are the first GPS satellites transmitting operational navigation signals on three distinct frequencies. We investigate apparent clock offset variations for the Block IIF-1 (SVN62) spacecraft that have been identified in L1/L2 clock solutions as well as the L1/L5-minus-L1/L2 clock difference. With peak-to-peak amplitudes of 10?C40?cm, these variations are of relevance for future precision point positioning applications and ionospheric analyses. A proper characterization and understanding is required to fully benefit from the quality of the new signals and clocks. The analysis covers a period of 8?months following the routine payload activation and is based on GPS orbit and clock products generated by the CODE analysis center of the International GNSS Service (IGS) as well as triple-frequency observations collected with the CONGO network. Based on a harmonic analysis, empirical models are presented that describe the sub-daily variation of the clock offset and the inter-frequency clock difference. These contribute to a better clock predictability at timescales of several hours and enable a consistent use of L1/L2 clock products in L1/L5-based positioning.     

Yaw attitude modeling for BeiDou I06 and BeiDou-3 satellites

The estimated yaw angles of the BeiDou I06 satellite demonstrated that the satellite experienced midnight- or noon-turn maneuvers when the sun elevation angle above the orbital plane (β angle) was in the range of [??3°, +?3°] and the orbital angle was in the range of approximately [??6°, 6°] or [174°, 186°]. The behavior of yaw attitude maneuvers in the vicinity of the midnight and noon points was identical. An alternative yaw attitude model similar to that used for the Galileo Full-Operation-Capacity (FOC) satellites was developed on the basis of the estimated BeiDou I06 yaw angles with an accuracy of approximately 3.4° to reproduce the yaw attitude behaviors. However, a discrepancy in the form of a reversal in yaw direction during the midnight-turn maneuver was observed for BeiDou I06 when the β angle was extremely small (<?0.1°). The derived yaw attitude model was proved to model the yaw attitude of the BeiDou-3 experimental satellites, and reduces the observation residuals in the vicinity of the midnight and noon points to normal levels, and facilitates continuous satellite clock estimation during eclipse periods. Compared to the yaw attitude model developed by the European Space Operations Centre (ESOC), a similar performance has been achieved with maximum yaw differences up to 9.2° when the β angle is close to 0°. The average agreement between the models is about 1°. However, the ESOC model was developed based on a patented eclipsing model, the developed model in this study is open access.     

An evaluation of solar radiation pressure strategies for the GPS constellation

The subtle effects of different Global Positioning System (GPS) satellite force models are becoming apparent now that mature processing strategies are reaching new levels of accuracy and precision. For this paper, we tested several approaches to solar radiation pressure (SRP) modeling that are commonly used by International GNSS Service (IGS) analysis centers. These include the GPS Solar Pressure Model (GSPM; Bar-Sever and Kuang in The Interplanetary Network Progress Report 42-160, 2005) and variants of the so-called DYB model (Springer et al. in Adv Space Res 23:673–676, 1999). Our results show that currently observed differences between GPS orbit solutions from the various IGS analysis centers are in large part explained by differences between their respective approaches to modeling SRP. DYB-based strategies typically generate orbit solutions that have the smallest differences with respect to the IGS final combined solution, largely because the DYB approach is most commonly used by the contributing analysis centers. However, various internal and external metrics, including ambiguity resolution statistics and satellite laser ranging observations, support continued use of the GSPM-based approach for precise orbit determination of the GPS constellation, at least when using the GIPSY-OASIS software.     

基于IGS精密星历的完备性监测指标SISA研究

在卫星导航定位系统中,空间信号精度SISA是反映完备性参数的重要指标之一。从SISE与SISA的基本概念出发,描述了SISA的计算算法,根据卫星星钟的不同,选取了不同类型的GPS卫星,并利用IGS精密星历,预报精密星历进行了SISE和SISA的计算,结果表明:IGR星历计算的SISA指标满足系统要求,SISA主要受卫星钟差的精度影响,不同原子钟类型的卫星计算出来的SISA值变化不大。该成果对于卫星导航的完备性监测数据处理具有参考意义。     

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.     

Estimation of elevation-dependent satellite antenna phase center variations of GPS satellites

A method for the estimation of the phase center variations of GPS satellite antennas using global GPS data is presented. First estimations have shown an encouraging repeatability from day to day and between satellites of the same block. Thus, two different satellite antenna patterns for Block II/IIA and for Block IIR with a range of about 4 cm and an accuracy of less than 1 mm could be found. The present approach allows the creation of a consistent set of receiver and satellite antenna patterns and phase center offsets. Thereby, it is possible to switch from relative to absolute phase center variations without a scale problem in global networks. This changeover has an influence on troposphere parameters, reduces systematic effects due to uncorrect antenna modeling and should diminish the elevation dependence of GPS results. AcknowledgmentsThe authors thank Prof. G. Seeber (University of Hannover) and Dr. G. Wübbena (Geo++ GmbH) and their groups for their kindness in making available the absolute field calibration results derived from robot measurements.     

姿态模式切换期间QZSS卫星轨道及其钟差产品特性分析

导航卫星姿态控制模式切换对精密定轨解算得到的轨道和钟差均有较大影响。本文首先从理论上分析了卫星偏航姿态及其对精密定轨的影响,然后分别以卫星激光测距检核和钟差多项式拟合的方法对IGS MGEX分析中心的QZSS卫星轨道和钟差产品精度进行评价,最后以谱分析方法和改进阿伦方差揭示了卫星钟差的周期特性。基于2014年全年的QZSS卫星轨道和钟差产品的研究表明,一年内有两次长约20 d的地影季,太阳角呈现半年周期的波动;QZSS卫星在低太阳角时有零偏保护,其卫星轨道和钟差精度都与太阳角有显著相关性;卫星钟差具有与轨道周期相近的周期项,且周期项振幅与太阳角的大小也具有相关性,表明现有的定轨策略存在不足。考虑到QZSS与目前北斗星座中IGSO和MEO卫星姿态控制模式的相似性,该结论对于研究我国BDS姿态切换期间的精密定轨有一定参考价值。     

Short-term analysis of GNSS clocks

A characterization of the short-term stability of the atomic frequency standards onboard GNSS satellites is presented. Clock performance is evaluated using two different methods. The first method derives the temporal variation of the satellite’s clock from a polynomial fit through 1-way carrier-phase measurements from a receiver directly connected to a high-precision atomic frequency standard. Alternatively, three-way measurements using inter-station single differences of a second satellite from a neighboring station are used if the receiver’s clock stability at the station tracking the satellite of interest is not sufficient. The second method is a Kalman-filter-based clock estimation based on dual-frequency pseudorange and carrier-phase measurements from a small global or regional tracking network. Both methods are introduced and their respective advantages and disadvantages are discussed. The analysis section presents a characterization of GPS, GLONASS, GIOVE, Galileo IOV, QZSS, and COMPASS clocks based on these two methods. Special focus has been set on the frequency standards of new generation satellites like GPS Block IIF, QZSS, and IOV as well as the Chinese COMPASS/BeiDou-2 system. The analysis shows results for the Allan deviation covering averaging intervals from 1 to 1,000 s, which is of special interest for real-time PPP and other high-rate applications like processing of radio-occultation measurements. The clock interpolation errors for different sampling rates are evaluated for different types of clocks and their effect on PPP is discussed.     

中国GPS增强系统性能分析及技术发展

本文提出了一种GPS卫星导航增强系统技术性能分析方法,即利用IGS提供的高精度GPS星历、卫星钟差数据和电离层数据作为外符和检测标准,检验增强系统提供的卫星星历改正数据、卫星钟差改正数据和电离层改正数据的精度。通过实测数据分析表明,我国GPS卫星导航增强系统目前的服务性能与国外同类系统具有一定的差距,主要的技术薄弱环节在于GPS卫星的精密定轨与钟差解算技术。     

EVA: GPS-based extended velocity and acceleration determination

In this work, a new GPS carrier phase-based velocity and acceleration determination method is presented that extends the effective range of previous techniques. The method is named ‘EVA’, and may find applications in fields such as airborne gravimetry when rough terrain or water bodies make difficult or impractical to set up nearby GPS reference receivers. The EVA method is similar to methods such as Kennedy (Precise acceleration determination from carrier phase measurements. In: Proceedings of the 15th international technical meeting of the satellite division of the Institute of Navigation. ION GPS 2002, Portland pp 962–972, 2002b) since it uses L1 carrier phase observables for velocity and acceleration determination. However, it introduces a wide network of stations and it is independent of precise clock information because it estimates satellite clock drifts and drift rates ‘on-the-fly’, requiring only orbit data of sufficient quality. Moreover, with EVA the solution rate is only limited by data rate, and not by the available precise satellite clocks data rate. The results obtained are more robust for long baselines than the results obtained with the reference Kennedy method. An advantage of being independent of precise clock information is that, beside IGS Final products, also the Rapid, Ultra-Rapid (observed) and Ultra-Rapid (predicted) products may be used. Moreover, the EVA technique may also use the undifferenced ionosphere-free carrier phase combination (LC), overcoming baseline limitations in cases where ionosphere gradients may be an issue and very low biases are required. During the development of this work, some problems were found in the velocity estimation process of the Kennedy method. The sources of the problems were identified, and an improved version of the Kennedy method was used for this research work. An experiment was performed using a light aircraft flying over the Pyrenees, showing that both EVA and the improved Kennedy methods are able to cope with the dynamics of mountainous flight. A RTK-derived solution was also generated, and when comparing the three methods to a known zero-velocity reference the results yielded similar performance. The EVA and the improved-Kennedy methods outperformed the RTK solutions, and the EVA method provided the best results in this experiment. Finally, both the improved version of the Kennedy method and the EVA method were applied to a network in equatorial South America with baselines of more than 1,770 km, and during local noon. Under this tough scenario, the EVA method showed a clear advantage for all components of velocity and acceleration, yielding better and more robust results.     

GPS/GALILEO偏航姿态异常对动态PPP的影响及其改正模型

当GPS、GALILEO卫星运行至与太阳、地球近似共线时,卫星很难维持名义姿态,将出现一段时间的偏航姿态异常。本文基于不同分析中心所提供的精密轨道和钟差产品,在卫星偏航姿态异常时期,设计不同姿态改正策略,选取全球分布的7个MGEX站10 d实测数据,分析了GPS、GALILEO卫星的天线相位中心改正、相位缠绕改正对观测值残差及动态PPP定位结果的影响。研究表明,在卫星偏航姿态异常时期,采用名义偏航姿态对GPS、GALILEO卫星观测值残差的影响可分别达到8和11 cm,在此期间,GPS/GALILEO卫星采用模型偏航姿态,与采用名义偏航姿态相比,动态PPP的E、N、U 3个方向的定位精度可分别提高13.30%、15.77%和12.98%,相较于剔除卫星策略,采用模型偏航姿态的动态PPP定位精度在E、N、U方向可分别提高5.399%、4.430%、5.992%。     

几种精密卫星钟差加密方法的比较与分析

卫星钟差是影响GPS高精度单点定位的一个重要因素,通过分析几种加密GPS精密卫星钟差的方法,将加密结果与GFZ(German Research Centre for Geosciences)提供的数据进行比较,并通过精度分析,得出分段线性插值法是加密GPS精密卫星钟差一种较为可靠的方法。