Swarm系列卫星非差运动学厘米级精密定轨

采用2015年5月24日—30日的Swarm星载GPS双频观测数据,基于Melbourne-Wübbena(MW)和消电离层线性组合,在精密单点定位技术的基础上,采用批处理最小二乘估计法对不同轨道高度的Swarm系列卫星进行非差运动学精密定轨。利用星载GPS相位观测值残差、与欧空局发布的简化动力学轨道对比,以及SLR检核3种方法对Swarm系列卫星非差运动学定轨结果进行精度评估。结果表明:①Swarm系列卫星星载GPS相位观测值残差RMS为6~7 mm;②与欧空局发布的简化动力学轨道进行求差,径向、切向及法向轨道差值RMS为2~4 cm;③与欧空局发布的运动学轨道进行求差,径向、切向及法向轨道差值RMS为1~2 cm;④SLR检核结果表明Swarm-A/B/C卫星轨道精度为3~4 cm。因此,采用非差运动学定轨方法与本文提供的定轨策略进行Swarm系列卫星精密定轨是切实可行的,定轨精度为厘米级。     

HY2A卫星的GPS/DORIS/SLR数据精密定轨

采用HY2A卫星2013年2月的实测数据,研究了GPS、星载多谱勒无线电定轨定位系统(DORIS)及卫星激光测距(SLR)三种观测数据的单独和联合定轨问题。通过与法国CNES的精密轨道数据比较发现:分别采用GPS、DORIS和SLR数据进行单独定轨,GPS数据确定轨道的径向平均精度为1.3cm,三维位置约为6.2cm;DORIS定轨的径向平均精度为1.6cm,比GPS结果略差;SLR确定轨道的径向平均精度为2.3cm。用GPS、DORIS和SLR三种数据联合定轨,确定轨道的径向平均精度为1.2cm,三维位置约为6.5cm。与星载GPS定轨结果比较,三种观测数据的联合定轨在提高卫星轨道确定精度上不明显,但联合定轨有利于保持计算轨道精度相对稳定。用站星间高度角大于60°的SLR数据检验GPS/DORIS联合确定的轨道,两者在测距方向的均方差为2.5cm,可见基于HY2A的观测数据可以实现cm级的定轨需求。     

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%。     

Implementation and testing of the gridded Vienna Mapping Function 1 (VMF1)

The new gridded Vienna Mapping Function (VMF1) was implemented and compared to the well-established site-dependent VMF1, directly and by using precise point positioning (PPP) with International GNSS Service (IGS) Final orbits/clocks for a 1.5-year GPS data set of 11 globally distributed IGS stations. The gridded VMF1 data can be interpolated for any location and for any time after 1994, whereas the site-dependent VMF1 data are only available at selected IGS stations and only after 2004. Both gridded and site-dependent VMF1 PPP solutions agree within 1 and 2 mm for the horizontal and vertical position components, respectively, provided that respective VMF1 hydrostatic zenith path delays (ZPD) are used for hydrostatic ZPD mapping to slant delays. The total ZPD of the gridded and site-dependent VMF1 data agree with PPP ZPD solutions with RMS of 1.5 and 1.8 cm, respectively. Such precise total ZPDs could provide useful initial a priori ZPD estimates for kinematic PPP and regional static GPS solutions. The hydrostatic ZPDs of the gridded VMF1 compare with the site-dependent VMF1 ZPDs with RMS of 0.3 cm, subject to some biases and discontinuities of up to 4 cm, which are likely due to different strategies used in the generation of the site-dependent VMF1 data. The precision of gridded hydrostatic ZPD should be sufficient for accurate a priori hydrostatic ZPD mapping in all precise GPS and very long baseline interferometry (VLBI) solutions. Conversely, precise and globally distributed geodetic solutions of total ZPDs, which need to be linked to VLBI to control biases and stability, should also provide a consistent and stable reference frame for long-term and state-of-the-art numerical weather modeling.     

Apropos laser tracking to GPS satellites

. Laser tracking to GPS satellites (PRN5 and 6) provides an opportunity to compare GPS and laser systems directly and to combine data of both in a single solution. A few examples of this are given in this study. The most important results of the analysis are that (1) daily SLR station coordinate solutions could be generated with a few cm accuracy; (2) coordinates of nine stations were determined in a 2.3-year-long arc solution; (3) the contribution of laser data on the `SLR-GPS' combined orbit, resulting from the simultaneous processing of SLR and GPS data, is significant and (4) laser-only orbits have an accuracy of 10–20 cm, 1-day predictions of SLR orbits differ from IGS orbits by about 20–40 cm, 2-day predictions by 50–60 cm. Received: 1 October 1996 / Accepted: 14 February 1997     

SWARM卫星简化动力学厘米级精密定轨

联合星载GPS双频观测值与简化的动力学模型,在卫星运动方程中引入适当的伪随机脉冲参数,对SWARM卫星进行精密定轨。采用星载GPS相位观测值残差、重叠轨道以及与外部轨道对比等3种方法对SWARM卫星简化动力学定轨结果进行检核。结果表明:SWARM星载GPS相位观测值残差RMS为7~10mm;径向、切向以及法向6h重叠轨道差值RMS均在1cm左右,3个方向均无明显的系统误差。通过与欧空局(ESA)发布的精密轨道进行对比分析,径向轨道差值RMS为2~5cm,切向轨道差值RMS为2~5cm,法向轨道差值RMS为2~4cm,3D轨道差值RMS为4~7cm;SWARM-B定轨精度优于SWARM-A与SWARM-C。因此,采用简化动力学法与本文提供的定轨策略进行SWARM卫星精密定轨是切实可行的,定轨结果良好且稳定,定轨精度达到厘米级。     

基于星载GPS数据的Jason-3卫星简化动力学和运动学法精密定轨

利用Jason-3星载GPS观测数据,采用简化动力学方法和运动学方法对Jason-3卫星进行精密定轨研究.通过载波相位残差、重叠轨道对比、参考轨道对比和卫星激光测距(SLR)轨道检核四种方式评定轨道精度.计算相位残差均方根(RMS)值,简化动力学轨道的RMS值在0.7～0.8cm,运动学轨道的RMS值在0.50～0.5...     

Combined processing of the IGS and the IGEX network

 Since the beginning of the International Global Navigation Satellite System (GLONASS) Experiment, IGEX, in October 1998, the Center for Orbit Determination in Europe (CODE) has acted as an analysis center providing precise GLONASS orbits on a regular basis. In CODE's IGEX routine analysis the Global Positioning System (GPS) orbits and Earth rotation parameters are introduced as known quantities into the GLONASS processing. A new approach is studied, where data from the IGEX network are combined with GPS observations from the International GPS Service (IGS) network and all parameters (GPS and GLONASS orbits, Earth rotation parameters, and site coordinates) are estimated in one processing step. The influence of different solar radiation pressure parameterizations on the GLONASS orbits is studied using different parameter subsets of the extended CODE orbit model. Parameterization with three constant terms in the three orthogonal directions, D, Y, and X (D = direction satellite–Sun, Y = direction of the satellite's solar panel axis), and two periodic terms in the X-direction, proves to be adequate for GLONASS satellites. As a result of the processing it is found that the solar radiation pressure effect for the GLONASS satellites is significantly different in the Y-direction from that for the GPS satellites, and an extensive analysis is carried out to investigate the effect in detail. SLR observations from the ILRS network are used as an independent check on the quality of the GLONASS orbital solutions. Both processing aspects, combining the two networks and changing the orbit parameterization, significantly improve the quality of the determined GLONASS orbits compared to the orbits stemming from CODE's IGEX routine processing. Received: 10 May 2000 / Accepted: 9 October 2000     

基于Kalman滤波的大规模GNSS网参数估计方法

针对传统方法存在的缺陷,研究了利用Kalman滤波技术进行大规模GNSS网参数（主要包括测站位置参数、卫星轨道参数及极移参数）估计的理论方法与关键技术,并利用40个全球均匀分布的IGS站多天的观测数据对理论成果进行了验证。结果表明,本文估计得到的测站位置参数与IGS结果各分量较差的RMS值分别为0.85、1.1、1.21 cm,得到的卫星轨道参数外推1 h后与IGS最终星历各分量较差的RMS值分别为9.8、8.6、7.2 cm,得到的极移参数与IERS结果的较差基本在1 mas之内;该方法具有较高的估值精度,可有效地用于GNSS网各类参数的估计。     

Precise orbit determination of the Fengyun-3C satellite using onboard GPS and BDS observations

The GNSS Occultation Sounder instrument onboard the Chinese meteorological satellite Fengyun-3C (FY-3C) tracks both GPS and BDS signals for orbit determination. One month’s worth of the onboard dual-frequency GPS and BDS data during March 2015 from the FY-3C satellite is analyzed in this study. The onboard BDS and GPS measurement quality is evaluated in terms of data quantity as well as code multipath error. Severe multipath errors for BDS code ranges are observed especially for high elevations for BDS medium earth orbit satellites (MEOs). The code multipath errors are estimated as piecewise linear model in \(2{^{\circ }}\times 2{^{\circ }}\) grid and applied in precise orbit determination (POD) calculations. POD of FY-3C is firstly performed with GPS data, which shows orbit consistency of approximate 2.7 cm in 3D RMS (root mean square) by overlap comparisons; the estimated orbits are then used as reference orbits for evaluating the orbit precision of GPS and BDS combined POD as well as BDS-based POD. It is indicated that inclusion of BDS geosynchronous orbit satellites (GEOs) could degrade POD precision seriously. The precisions of orbit estimates by combined POD and BDS-based POD are 3.4 and 30.1 cm in 3D RMS when GEOs are involved, respectively. However, if BDS GEOs are excluded, the combined POD can reach similar precision with respect to GPS POD, showing orbit differences about 0.8 cm, while the orbit precision of BDS-based POD can be improved to 8.4 cm. These results indicate that the POD performance with onboard BDS data alone can reach precision better than 10 cm with only five BDS inclined geosynchronous satellite orbit satellites and three MEOs. As the GNOS receiver can only track six BDS satellites for orbit positioning at its maximum channel, it can be expected that the performance of POD with onboard BDS data can be further improved if more observations are generated without such restrictions.     

基于Helmert变换的GPS动力学轨道平滑(英文)

Orbit fitting is used in many GPS applications. For example, in Precise Point Positioning (PPP), GPS orbits (SP3 orbits) are normally retrieved either from IGS or from one of its Analysis Centers (ACs) with 15 minutes’ sampling, which is much bigger than the normal observation sampling. Therefore, algorithms should be derived to fit GPS orbits to the observation time. Many methods based on interpolation were developed. Using these methods the orbits fit well at the sampling points. However, these methods ig...     

Orbit fitting based on Helmert transformation

Orbit fitting is used in many GPS applications. For example, in Precise Point Positioning (PPP), GPS orbits (SP3 orbits) are normally retrieved either from IGS or from one of its Analysis Centers (ACs) with 15 minutes’ sampling, which is much bigger than the normal observation sampling. Therefore, algorithms should be derived to fit GPS orbits to the observation time. Many methods based on interpolation were developed. Using these methods the orbits fit well at the sampling points. However, these methods ignore the physical motion model of GPS satellites. Therefore, the trajectories may not fit the true orbits at the periods in between 2 sampling epochs. To solve this problem, we develop a dynamic approach, in which a model based on Helmert transformation is developed in GPS orbit fitting. In this orbit fitting approach, GPS orbits at sampling points are treated as pseudo-observations. Thereafter, Helmert transformation is built up between the pseudo-observations and dynamically integrated orbits at each epoch. A set of Helmert parameters together with corrections of GPS initial orbits are then modeled as unknown parameters. Results show that the final fit orbits have the same precision as the IGS final orbits.     

GNSS实时精密轨道快速计算方法及服务

多系统全球导航卫星系统（Global Navigation Satellite System,GNSS）精密轨道确定及其预报是实现高精度实时精密定位的前提。针对多GNSS系统超快速轨道解算时效性及轨道预报精度随时间下降的问题,提出一种基于分块递推最小二乘配置方法,该方法通过对动力学和几何学待估参数松弛、连接以及轨道状态参数转移递推,能够同时兼容事后及实时滤波定轨方法。该方法能够有效地提高多GNSS系统轨道解算效率,缩短实时轨道更新时间。基于全球实测数据验证了该方法的可靠性和有效性,轨道精度优于国际GNSS服务组织发布的GPS超快速轨道及德国地学研究中心发布的超快速轨道,实验结果表明,采用该方法,GPS/GLONASS/Galileo/BDS四系统120个地面测站精密定轨可以实现1 h更新,延迟30 min发布,统计GPS/GLONASS/Galileo/BDS实时轨道可用部分3D均方根分别为2.8 cm、8.5 cm、5.0 cm及11.5 cm（IGSO/MEO）。目前,1 h更新多GNSS系统轨道及实时产品服务系统已业务化发布,较之前发布的3 h更新及6 h更新轨道分别有20%~40%的精度提升。     

Homogeneous reprocessing of GPS,GLONASS and SLR observations

The International GNSS Service (IGS) provides operational products for the GPS and GLONASS constellation. Homogeneously processed time series of parameters from the IGS are only available for GPS. Reprocessed GLONASS series are provided only by individual Analysis Centers (i. e. CODE and ESA), making it difficult to fully include the GLONASS system into a rigorous GNSS analysis. In view of the increasing number of active GLONASS satellites and a steadily growing number of GPS+GLONASS-tracking stations available over the past few years, Technische Universität Dresden, Technische Universität München, Universität Bern and Eidgenössische Technische Hochschule Zürich performed a combined reprocessing of GPS and GLONASS observations. Also, SLR observations to GPS and GLONASS are included in this reprocessing effort. Here, we show only SLR results from a GNSS orbit validation. In total, 18 years of data (1994–2011) have been processed from altogether 340 GNSS and 70 SLR stations. The use of GLONASS observations in addition to GPS has no impact on the estimated linear terrestrial reference frame parameters. However, daily station positions show an RMS reduction of 0.3 mm on average for the height component when additional GLONASS observations can be used for the time series determination. Analyzing satellite orbit overlaps, the rigorous combination of GPS and GLONASS neither improves nor degrades the GPS orbit precision. For GLONASS, however, the quality of the microwave-derived GLONASS orbits improves due to the combination. These findings are confirmed using independent SLR observations for a GNSS orbit validation. In comparison to previous studies, mean SLR biases for satellites GPS-35 and GPS-36 could be reduced in magnitude from \(-35\) and \(-38\)  mm to \(-12\) and \(-13\)  mm, respectively. Our results show that remaining SLR biases depend on the satellite type and the use of coated or uncoated retro-reflectors. For Earth rotation parameters, the increasing number of GLONASS satellites and tracking stations over the past few years leads to differences between GPS-only and GPS+GLONASS combined solutions which are most pronounced in the pole rate estimates with maximum 0.2 mas/day in magnitude. At the same time, the difference between GLONASS-only and combined solutions decreases. Derived GNSS orbits are used to estimate combined GPS+GLONASS satellite clocks, with first results presented in this paper. Phase observation residuals from a precise point positioning are at the level of 2 mm and particularly reveal poorly modeled yaw maneuver periods.     

IGR轨道外推及其在对流层实时监测中的应用研究

在精密定轨和轨道外推理论的基础上，对比分析了由IGS快速轨道（IGS rapid orbit，IGR）外推1d的轨道与IGS最终轨道间的差异，用此外推轨道和最终轨道分别按长基线双差网解方案计算了BRUS、IENG、PTBB等欧洲一些IGS／TAI站的天顶湿延迟（zenith wet delay，ZWD），并对两种轨道下各站的ZWD进行了对比。轨道对比结果显示，除异常卫星外，“病态”卫星的IGR外推1d的轨道仍具有一定的高精度。而ZWD的对比结果则表明，当剔除异常卫星的影响并对“病态”卫星加以适当处理后，IGR外推1d的轨道用于对流层实时监测是可行的。     

远海实时GPS潮汐的实时精密单点定位观测

针对传统事后精密单点定位技术的时间延迟问题,该文基于IGS RTS实时数据流产品,开展了实时精密单点定位技术在远海实时GPS验潮中的应用研究.对RTS改正的实时精密卫星轨道和钟差进行了精度验证和分析,给出了RT-PPP的数据处理策略以及实时GPS验潮的基本流程;组织和实施了渤海湾船载GPS验潮试验,以压力式验潮仪数据为参考,对远距离实时GPS潮汐测量结果进行了精度分析.结果表明:①以IGS最终卫星轨道和钟差产品为参考,RTS实时精密卫星轨道在X、y、Z方向的精度(RMS)均优于3 cm,卫星钟差的精度优于0.15 ns;②采用傅里叶低通滤波方法,消除波浪对潮汐观测的影响,进一步提取潮位信息.在忽略船体姿态改正的情况下,实时精密单点定位验潮相对于压力式验潮仪结果的最大偏差优于20 cm,RMS达到7.5 cm.     

实时精密单点定位的远海实时GPS潮汐观测

针对IGS实时数据流产品,该文开展了实时精密单点定位技术在远海实时GPS验潮中的应用研究。对RTS改正的实时精密卫星轨道和钟差进行了精度验证和分析,给出了实时精密单点定位的数据处理策略以及实时GPS验潮的基本流程;组织和实施了渤海湾船载GPS验潮试验,以压力式验潮仪数据为参考,对远距离实时GPS潮汐测量结果进行了精度分析。结果表明:以IGS最终卫星轨道和钟差产品为参考,实时数据流产品实时精密卫星轨道在X、Y、Z方向的精度均优于3cm,卫星钟差的精度优于0.15ns;采用傅里叶低通滤波方法,消除波浪对潮汐观测的影响,进一步提取潮位信息。在忽略船体姿态改正的情况下,实时精密单点定位验潮相对于压力式验潮仪结果的最大偏差优于20cm,RMS达到7.5cm。     

多模全球导航卫星系统融合精密定轨

基于武汉大学自主研制的卫星导航系统综合处理软件(PANDA),利用全球实测的GPS/GLONASS、GPS/Gali-leo试验卫星(GIOVE)多模接收机数据进行GPS、GLONASS、GIOVE卫星的融合精密定轨理论与方法研究。通过与IGS提供的GPS与GLONASS卫星精密轨道比较、轨道重叠弧段互差以及SLR观测数据检核等多种方法对融合计算的精密轨道精度进行了评定。     

Joint estimation of GPS/BDS real-time clocks and initial results

We present the joint estimation model for Global Positioning System/BeiDou Navigation Satellite System (GPS/BDS) real-time clocks and present the initial satellite clock solutions determined from 106 stations of the international GNSS service multi-GNSS experiment and the BeiDou experimental tracking stations networks for 1 month in December, 2012. The model is shown to be efficient enough to have no practical computational limit for producing 1-Hz clock updates for real-time applications. The estimated clocks were assessed through the comparison with final clock products and the analysis of post-fit residuals. Using the estimated clocks and corresponding orbit products (GPS ultra-rapid-predicted and BDS final orbits), the root-mean-square (RMS) values of coordinate differences from ground truth values are around 1 and 2–3 cm for GPS-only and BDS-only daily mean static precise point positioning (PPP) solutions, respectively. Accuracy of GPS/BDS combined static PPP solutions falls in between that of GPS-only and BDS-only PPP results, with RMS values approximately 1–2 cm in all three components. For static sites, processed in the kinematic PPP mode, the daily RMS values are normally within 4 and 6 cm after convergence for GPS-only and BDS-only results, respectively. In contrast, the combined GPS/BDS kinematic PPP solutions show higher accuracy and shorter convergence time. Additionally, the BDS-only kinematic PPP solutions using clock products derived from the proposed joint estimation model were superior compared to those computed using the single-system estimation model.     

Zero-difference GPS ambiguity resolution at CNES–CLS IGS Analysis Center

CNES (Centre National d’Etudes Spatiales) and CLS (Collecte Localisation Satellites) became an International GNSS Service (IGS) Analysis Center (AC) the 20th of May 2010. Since 2009, we are using the integer ambiguity fixing at the zero-difference level strategy in our software package (GINS/Dynamo) as an alternative to classical differential approaches. This method played a key role among all the improvements in the GPS processing we made during this period. This paper provides to the users the theoretical background, the strategies and the models used to compute the products (GPS orbits and clocks, weekly station coordinate estimates and Earth orientation parameters) that are submitted weekly to the IGS. The practical realization of the two-step, ambiguity-fixing scheme (wide-lane and narrow-lane) is described in detail. The ambiguity fixing improved our orbit overlaps from 6 to 3?cm WRMS in the tangential and normal directions. Since 2008, our products have been also regularly compared to the IGS final solutions by the IGS Analysis Center Coordinator. The joint effects of ambiguity fixing and dynamical model changes (satellite solar radiation pressure and albedo force) improved the consistency with IGS orbits from 35 to 18?mm 3D-WRMS. Our innovative strategy also gives additional powerful properties to the GPS satellite phase clock solutions. Single receiver (zero-difference) ambiguity resolution becomes possible. An overview of the applications is given.