利用GPS观测数据研究高频地球自转参数

介绍了利用GPS解算地球自转参数的理论与方法,选取80个左右的IGS参考站以8h观测弧段进行GPS数据处理,分析解算了2h时间间隔地球自转参数序列值。对以上基于GPS观测解算的高频时间序列值进行小波降噪,并利用FFT频谱分析的结果获得了地球自转参数的周日项、半周日项。     

VLBI技术CONT观测的高频ERP解算

国际VLBI测天测地服务机构(IVS)已组织了多次VLBI连续加密观测(CONT),提供了高精度连续的原始观测数据,在地球自转参数(ERP)的连续高频解算中起到积极的作用,揭示了地球自转高频变化的观测资料和理论模型之间的差异,有助于进一步解析其激发机制改进模型.这里使用VLBI资料处理软件系统OCCAM处理了CONT02,CONT05和CONT08数据,并进行ERP高频解算及频谱分析.从各次CONT观测的残差频谱中发现较强周期信号,反映了地球自转的特性.特别是CONT08残差频谱中存在明显的周日项信患,揭示了北半球夏季月份大气激发对地球自转的作用.     

基于IGS网络资源的完备性监测研究

本文从国际知名的高精度数据处理软件BERNESE软件入手,选择欧洲三个IGS跟踪站2002年1年的原始观测数据,应用IGS网络资源包括GPS精密星历,GPS精密钟差,地球自转参数文件,进行精密单点定位计算,讨论卫星轨道对完备性监测的影响,并将IGS跟踪站视为监测站和用户站,实现对IGS监测站和用户站的完备性监测。     

VLBI和GPS观测联合解算地球自转参数和日长变化

利用GAMIT软件处理了2005年9月、2006年1月和2月每天17个IGS国际GNSS服务(International GNSS Service,IGS)站的GPS观测数据,获得了地球自转参数(Earth rotation parameter,ERP)和日长变化(ΔLOD),并与IGS综合解进行了对比。利用OCCAM 6.2软件处理了相同时间内的甚长基线干涉测量(very long baseline interferometry,VLBI)观测数据,将解算结果与国际VLBI服务(International VLBI Service,IVS)结果进行了对比。分别采用基于内符合精度和IERS 08C04序列的定权方法对VLBI解算结果与GPS解算结果进行了加权平均,获得了VLBI和GPS技术对ERP和ΔLOD的协议结果。研究结果表明,VLBI解算结果与采用的插值方法有较大的关系,基于IERS 08C04序列的加权平均方法达到了利用VLBI解算结果对GPS解算结果进行修正的目的。VLBI和GPS技术的联合弥补了VLBI观测数据密度不够和GPS解算结果不稳定的缺陷,使解算结果的稳定性和可靠性有所提高。     

联合GNSS和SLR观测对地球自转参数的解算与分析

选取不同数量的IGS站,分别利用GPS和GPS+GLONASS观测数据计算ERP参数,并将结果和IGS公布值进行比较,分析测站数量增加和加入GLONASS观测数据对解算ERP参数的影响。此外,还利用GLONASS卫星的全部卫星激光测距(SLR)数据进行ERP参数解算,并将SLR结果和GNSS解算结果联合计算ERP,结果表明,联合SLR可改善GNSS数据解算ERP参数及高频ERP参数的系统性误差影响和稳定性。     

多分析中心站坐标产品的综合方法研究

首先介绍了多分析中心产品融合处理的两种综合策略，然后基于解层面的综合策略，提出了站坐标和地球自转参数同时综合的方法。采用国际GNSS服务组织（International GNSS Service，IGS）9个分析中心1 a的数据进行试验，从站坐标、地球自转参数精度以及地心运动3个方面验证了该方法的正确性。结果表明，基于综合方法得到的综合解和IGS综合解处于同一精度水平。站坐标在平面和高程方向的一致性分别为0.5 mm和1.0 mm，极移和极移速率的一致性分别优于7.0×10-6"和40.0×10-6"/d，日长参数优于7.7×10-6 s/d。所提出的综合方法可用于全球连续监测评估系统（international GNSS monitoring and assessment system，iGMAS）的站坐标/地球自转参数产品综合。     

iGMAS全球电离层延迟模型及并行计算策略

针对全球电离层延迟建模中传统串行处理方法效率低等问题,研究了基于全球分布的IGS跟踪站和iGMAS跟踪站观测数据实现全球电离层建模并行解算的基本方法、流程及策略。在Bernese软件基础上研制了一套iGMAS全球电离层延迟建模软件。为了验证并行解算方法的正确性和计算效率,利用全球200个左右IGS跟踪站和6个iGMAS跟踪站2014-08-20-2014-09-06共7周的观测数据,解算了快速电离层TEC格网。与IGS,CODE以及ESA最终电离层格网比较,结果表明:基于该方法解算的快速电离层TEC格网,与CODE,ESA以及IGS最终电离层TEC格网的互差,统计不同纬度带内偏差的均方根误差,全球范围内偏差的均方根误差均在1.5~2.5 TECu之间,南北半球高纬度地区在0.5~1.5 TECu之间,所有地区均优于5 TECu,整体精度与IGS,CODE以及ESA最终电离层TEC格网精度产品相当。     

ERP精度对“嫦娥一号”差分VLBI定位精度的影响

推导了差分VLBI（ΔDOR）用于月球探测器定位的数学模型,利用嫦娥一号实测数据对探测器的位置和地球自转参数（ERP）同时进行了解算,并分析了ERP先验精度对定位精度的影响。     

由CHAMP星载GPS相位双差数据解算地球引力场模型

利用7d的CHAMP星载GPS相位观测数据和48个IGS跟踪站的观测数据，构造星地双差相位观测量，进行GPS数据预处理；利用Cowell Ⅱ数值法进行轨道积分和分块Bayes最小二乘参数估计，解算了地球引力场位系数。该模型与EGM96相比（70阶次），大地水准面起伏差异最大为2．872m，差弄精度为0．522m，平均差异为-0．003m，这说明本文解算的地球重力场模型与EGM96没有系统性差异。     

测站数目及分布对地球自转参数解算精度的影响

为研究GPS数据解算地球自转参数(ERP)精度受测站数目及分布均衡性影响规律的问题,该文利用全球国际GNSS服务(IGS)站提供的GPS数据,设置不同测站数、不同测站分布均衡程度的解算策略,通过对比不同策略下ERP解算精度,来研究测站数目和均衡程度对ERP解算过程中的影响规律。结果表明,考虑到解算效率的情况下,测站数目选择40个时能达到最佳效果,此时极移在x方向的RMS值为0.223 081 mas,在y方向的RMS值为0.186 941 mas;对于测站分布均衡性,该文提出用观测网的网重心坐标转换为大地坐标作为评价指标,当网重心越接近地心,解算精度越高。研究成果表明在利用GPS数据解算ERP参数时,选择适当数目的测站以及分布均衡性好的解算策略可以提高解算效率及精度。     

基于VLBI与BDS联合解算ERP的算法研究

利用PANDA软件解算2016年第61~91天的MGEX(Multi-GNSS Experiment)服务站的北斗数据,获得地球自转参数(ERP)。利用VieVS2.2软件处理了同时段的甚长基线干涉测量(VLBI)数据。采用基于IERS 08C04序列的定权方法对BDS和VLBI的解算结果进行加权平均,得到综合的ERP值。结果表明,与IERS比较,极移在X方向差值的RMS为0.249 mas,Y方向差值的RMS为0.296 mas,UT1-UTC差值的RMS为0.053 ms.利用BDS与VLBI数据对ERP参数进行联合解算,弥补了BDS解算结果不稳定和VLBI观测不连续的缺陷,使解算结果的稳定性和可靠性均有所提高。      

利用空间大地测量数据探测地球膨胀效应

地球自转服务局(IERS)采用多种高精度的空间探测技术综合解算得到的国际地球参考框架(ITRF)是国际上公认的精度高、稳定性好的参考框架。为了研究地球的膨胀或收缩效应,本文采用ITRF2000的站坐标和速度,利用Delaunay算法生成的三角网逼近地球形体,计算出了地球的体积变化。     

利用2006—2015年VLBI数据进行地球定向参数解算与分析

目前正处在下一代甚长基线干涉测量（very long baseline interferometry，VLBI）系统的建设时期。利用维也纳VLBI与卫星软件（Vienna VLBI and satellite software，VieVS）解算了2006—2015年的VLBI数据，得到了10 a的地球定向参数（Earth orientation parameters，EOP）时间序列，并与国际地球自转服务机构的结果进行了对比。利用解算结果得到了10 a的日长变化时间序列，通过傅里叶分析得出了日长变化的短周期、半月周期、月周期、半年周期和周年周期，同时还分析得到了极移序列中的周年项和张德勒周期项以及章动改正序列中的自由核章动项。此次解算工作可为武汉大学卫星台站日后的VLBI数据解析积累一定的经验。     

GGOS-D: homogeneous reprocessing and rigorous combination of space geodetic observations

In preparation of activities planned for the realization of the Global Geodetic Observing System (GGOS), a group of German scientists has carried out a study under the acronym GGOS-D which closely resembles the ideas behind the GGOS initiative. The objective of the GGOS-D project was the investigation of the methodological and information-technological realization of a global geodetic-geophysical observing system and especially the integration and combination of the space geodetic observations. In the course of this project, highly consistent time series of GPS, VLBI, and SLR results were generated based on common state-of-the-art standards for modeling and parameterization. These series were then combined to consistently and accurately compute a Terrestrial Reference Frame (TRF). This TRF was subsequently used as the basis to produce time series of station coordinates, Earth orientation, and troposphere parameters. In this publication, we present results of processing algorithms and strategies for the integration of the space-geodetic observations which had been developed in the project GGOS-D serving as a prototype or a small and limited version of the data handling and processing part of a global geodetic observing system. From a comparison of the GGOS-D terrestrial reference frame results and the ITRF2005, the accuracy of the datum parameters is about 5?C7?mm for the positions and 1.0?C1.5?mm/year for the rates. The residuals of the station positions are about 3?mm and between 0.5 and 1.0?mm/year for the station velocities. Applying the GGOS-D TRF, the offset of the polar motion time series from GPS and VLBI is reduced to 50 ??as (equivalent to 1.5?mm at the Earth??s surface). With respect to troposphere parameter time series, the offset of the estimates of total zenith delays from co-located VLBI and GPS observations for most stations in this study is smaller than 1.5?mm. The combined polar motion components show a significantly better WRMS agreement with the IERS 05C04 series (96.0/96.0???as) than VLBI (109.0/100.7???as) or GPS (98.0/99.5???as) alone. The time series of the estimated parameters have not yet been combined and exploited to the extent that would be possible. However, the results presented here demonstrate that the experiences made by the GGOS-D project are very valuable for similar developments on an international level as part of the GGOS development.     

Combination procedure for length-of-day time series according to the noise frequency behavior

 Length-of-day (LOD) estimates from seven Global Positioning System (GPS) and three satellite laser ranging (SLR) analysis centers were combined into an even-spaced time series for a 27-month period (1996–1998). This time series was compared to the multi-technique Earth-orientation-parameter (EOP) combined solution (C04) derived at the Central Bureau of the International Earth Rotation Service (IERS/CB). Due to inhomogeneities in the different series derived from the various techniques (time, length, quality, and spatial resolution), the concept of a combined solution is justified. The noise behavior in LOD for different techniques varies with frequency; the data series were divided into frequency windows after removing both biases and trends. Different weight factors were assigned in each window, discriminating by technique, and produced one-technique combined solutions. Finally, these one-technique combined solutions were combined to obtain the final multi-technique solution. The LOD combined time series obtained by the present method based on the frequency windows combined series (FWCS) is very close to the IERS C04 solution. It could be useful to generate a new LOD reference time series to be used in the study of high-frequency variations of Earth rotation. Received: 28 March 2000 / Accepted: 15 February 2001     

Two-step method for the determination of the differential code biases of COMPASS satellites

The differential code bias (DCB) in satellites of the Global Navigation Satellite Systems (GNSS) should be precisely corrected when designing certain applications, such as ionospheric remote sensing, precise point positioning, and time transfer. In the case of COMPASS system, the data used for estimating DCB are currently only available from a very limited number of global monitoring stations. However, the current GPS/GLONASS satellite DCB estimation methods generally require a large amount of geographically well-distributed data for modeling the global ionospheric vertical total electron content (TEC) and are not particularly suitable for current COMPASS use. Moreover, some satellites with unstable DCB (i.e., relatively large scatter) may affect other satellite DCB estimates through the zero-mean reference that is currently imposed on all satellites. In order to overcome the inadequacy of data sources and to reduce the impact of unstable DCB, a new approach, designated IGGDCB, is developed for COMPASS satellite DCB determination. IGG stands for the Institute of Geodesy and Geophysics, which is located in Wuhan, China. In IGGDCB, the ionospheric vertical TEC of each individual station is independently modeled by a generalized triangular series function, and the satellite DCB reference is selected using an iterative DCB elimination process. By comparing GPS satellite DCB estimates calculated by the IGGDCB approach based on only a handful (e.g., seven) of tracking stations against that calculated by the currently existing methods based on hundreds of tracking stations, we are able to demonstrate that the accuracies of the IGGDCB-based DCB estimates perform at the level of about 0.13 and 0.10?ns during periods of high (2001) and low (2009) solar activity, respectively. The iterative method for DCB reference selection is verified by statistical tests that take into account the day-to-day scatter and the duration that the satellites have spent in orbit. The results show that the impact of satellites with unstable DCB can be considerably reduced using the IGGDCB method. It is also confirmed that IGGDCB is not only specifically valid for COMPASS but also for all other GNSS.     

Earth orientation parameters estimated from VLBI during the CONT11 campaign

In this paper we investigate the accuracy of the earth orientation parameters (EOP) estimated from the continuous VLBI campaign CONT11. We first estimated EOP with daily resolution and compared these to EOP estimated from GNSS data. We find that the WRMS differences are about 31  $\upmu $ as for polar motion and 7  $\upmu $ s for length of day. This is about the precision we could expect, based on Monte Carlo simulations and the results of the previous CONT campaigns. We also estimated EOP with hourly resolution to study the sub-diurnal variations. The results confirm the results of previous studies, showing that the current IERS model for high-frequency EOP variations does not explain all the sub-diurnal variations seen in the estimated time series. We then compared our results to various empirical high-frequency EOP models. However, we did not find that any of these gave any unambiguous improvement. Several simulations testing the impact of various aspects of, e.g. the observing network were also made. For example, we made simulations assuming that all CONT11 stations were equipped with fast VLBI2010 antennas. We found that the WRMS error decreased by about a factor five compared to the current VLBI system. Furthermore, the simulations showed that it is very important to have a homogenous global distribution of the stations for achieving the highest precision for the EOP.     

利用LS_AR模型对极移参数的中长期预报

为了满足深空探测器自主导航定位对极移参数中长期预报的需求,阐述了基于LS_AR模型的极移参数中长期预报和精度评定的原理,提出了4种改进方案对LS_AR模型的构建进行优化,并利用IERS提供的1990～1996年的极移参数的时间序列检验4种优化方案,得到了最优的预报模型,在400 d跨度上,其预报结果的平均绝对误差比未优化的模型小3 mas左右。     

COMPASS地球参考框架的建立和维持

探讨了COMPAAS地球参考框架所对应的地球参考系统的定义、地球参考框架的建立和维持方案,并选取我国现有的COMPASS跟踪站和IGS跟踪站数据进行了一系列的仿真验证。分析了跟踪站分布对地球参考框架的影响以及多分析中心组合对参考框架周解稳定性和可靠性的影响。试验结果表明,全球均匀分布的跟踪站在20个左右时就能基本满足建立地球参考框架的需求,为了保障地球参考框架的长期稳定性和精度,全球均匀分布的跟踪站应为30个左右;引入多个分析中心组合处理生成最终解不但可以探测出单个数据处理中心所引入的粗差,也可进一步增强周解的稳定性和可靠性。     

IGS Earth Rotation Parameters

Since its official start in January 1994, the International GPS Service (IGS) has been distributing, as part of its product combination, two distinct Earth rotation parameter (ERP) series: the IGS Rapid series and the IGS Final series. Initially, the IGS Rapid ERP values were interpolations of the International Earth Rotation Service (IERS) Bulletin A, whereas the IGS Final ERP series was based on the IERS Bulletin B. Since June 1996, the IGS has been generating its own Final ERP series consistent with the IGS combined orbit products and based on weighted means of individual IGS analysis center (AC) solutions. At first, only the polar motion (PM) coordinates and their rates were combined. Length of Day (LOD) and Universal Time (UT) solutions, also based on separate weighted mean combinations, followed in March 1997. Currently, the IGS Rapid and Final combinations are produced and made available within 17 hours and 11 days, respectively, after the last observation. Both IGS and the best AC series are consistent and precise at the 0.1-milliarcsecond (mas) level for PM and at about 30 μs for LOD. Biases in some AC solutions may exceed these consistency levels. Comparisons of both IGS ERP series with external standards, such as the IERS multitechnique Bulletins and atmospheric angular momentum series, confirm the estimated precisions. ? 1999 John Wiley & Sons, Inc.