顾及非线性变化的地球参考框架建立与维持的思考

介绍了全球及区域参考框架(ITRF)的发展现状,探讨了国际地球参考框架(ITRF)的不足,给出了使用顾及非线性变化的速度模型建立和维持参考框架的思路,并分析了解决不同类型非线性变化的处理方法。在此基础上,探讨了区域参考框架建立与维持的方法。最后展望了全球大地观测系统背景下的未来ITRF发展趋势。     

WGS84与ITRF框架间的坐标转换分析

简要介绍了WGS84坐标系和ITRF框架,给出了不同ITRF框架间的坐标转换流程,并利用实例对WGS84与ITRF框架间的转换关系进行了验证分析.结果表明,ITRF2008与WGS84坐标基本一致,但由于ITRF框架的站速度对站坐标的影响与时间成正相关,当需要采用ITRF框架时,应选用最新的国际地球参考框架.     

论地球参考框架的维持

地球参考框架的维持,是指给出随时间变化的框架点坐标.一般有两种形式,即给出速度模喇或坐标时问序列.全球参考框架(主要指ITRF)的维持小断趋于完善,区域参考框架一般足指对ITRF的加密.由于区域参考框架的维持上要是通过GPS技术建市CORS网来实现的,因而对IGS实现和维持的参考框架作简要介绍.讨论Ⅸ域参考框架维持的两种技术途径,并对我围地球参考框架的维持提出儿点看法.     

地球参考框架联合解算方法

地球参考框架是地球坐标系统的实现,本文系统介绍了地球参考框架的确定方法,结合ITRF2005参考框架,详细介绍了地球参考框架联合解算的方法,并对地球参考框架的质量分析方法进行了探讨。     

ITRF2005的实现与改进

2006年10月,IERS发布了继ITRF2000以后最新版本的国际地球参考框架——ITRF2005。ITRF2005的实现是基于多种空间技术解时间序列的联合处理而得到的,这在ITRF的实现历史中还是第一次。此外,ITRF2005除了包含作为参考框架体现的站点坐标和速率之外,还包含一起参与联合处理的地球定向参数:极移,极移速率,日长,UT1的时间序列。主要介绍ITRF2005的实现的基本情况及其相对于以前版本的改进。     

ITRF2005:基于测站位置和地球自转参数的最新地球参考框架

与之前的国际地球参考框架(ITRF)将全球长期解作为输入数据进行组合不同,ITRF2005将测站坐标(卫星技术每星期的数据和VLBI每24小时的数据)和每天的地球自转参数(EOPs)作为输入数据。使用测站位置时间序列的优势在于可以监控测站的非线性运动和非连续性,并检验框架物理参数即原点和尺度的时变特性。ITRF2005原点定义为:相对于由SLR技术13年的观测数据所得的地球质心的平移和平移速度为零;尺度定义为:相对于由VLBI技术26年的观测数据所得的尺度及其变化率为零;ITRF2005的定向(2000.0历元)及其速率与ITRF2000中70个高质量的测站一致。ITRF2005原点(2000.0历元)及其速率相对于ITRF2000沿X,Y,Z轴在0.1,0.8,5.8mm和0.2,0.1,1.8mm/y的水平上一致,其分量的误差分别为0.3mm和0.3mm/y。两个参考框架原点间一致性差可能是因为SLR网的几何图形差。ITRF2005组合中包含了84个并置站,尺度的不一致性在2000.0历元为1ppb(赤道处为6.3mm),SLR和VLBI由各自时间序列堆栈得到的长期解之间尺度不一致性为0.08ppb/yr。这些不一致性可能是因为SLR和VLBI网形差、并置站质量不好、局部联系的不确定性、系统误差影响以及数据分析中模型改正的不一致性。ITRF历史上,ITRF2005第一次采用了紧组合的方式给出了与之相一致的EOP序列,包括由VLBI和卫星技术得到的极移和仅从VLBI得到的UT和日长数据。     

NNR-NUVEL-1A板块运动模型

NNR-NUVEL-1A模型作为国际地球参考框架(ITRF)的国际标准模型应该是一个足够精确和严密的整网无残余旋转模型(NNR)。本文中采用了新方法精确计算了全球板块运动的惯性张量Q以及板块的总角动量。结果表明附加NNR-NUVEL-1A约束的参考框架相对于地球岩石圈有0.012°/Ma的旋转,采用NNR-NUVEL-1A模型并没有完全实现NNR约束。本文最后提出了新的NNR-NUVEL-1B这一简化模型。     

ITRF的发展及其在建立和维持地区性大地坐标系中的作用

简要介绍了 ITRF(国际地球参考框架 )的发展情况以及它与 IGS的关系 ,讨论了它在建立和维持地区性地球参考系的作用 ,最后展望了它在建立和维持我国地心参考系中的作用     

ITRF2005参考框架的综合探讨

对近年来精度高、应用较多的国际地球参考架ITRF2005做了简单概述,指出了ITRF2000与ITRF2005之间在解的生成、基准的定义和实现等方面的不同。此外,ITRF2005除了包含作为参考框架体现的站点坐标和速率之外,还包含一起参与联合处理的地球定向参数：极移、极移速率、日长、UT1的时间序列[1]。重点阐明ITRF2005的实现的基本情况及其相对于ITRF2000所作改进的理由和合理性。     

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

详细总结了ITRF等参考框架的原理,分析了其优缺点和特性并分析了不同地球参考框架的区别。本文在协议研究的基础上,结合我国实际,提出了建立我国新一代地球参考系统的若干建议。     

ITRF2008: an improved solution of the international terrestrial reference frame

ITRF2008 is a refined version of the International Terrestrial Reference Frame based on reprocessed solutions of the four space geodetic techniques: VLBI, SLR, GPS and DORIS, spanning 29, 26, 12.5 and 16?years of observations, respectively. The input data used in its elaboration are time series (weekly from satellite techniques and 24-h session-wise from VLBI) of station positions and daily Earth Orientation Parameters (EOPs). The ITRF2008 origin is defined in such a way that it has zero translations and translation rates with respect to the mean Earth center of mass, averaged by the SLR time series. Its scale is defined by nullifying the scale factor and its rate with respect to the mean of VLBI and SLR long-term solutions as obtained by stacking their respective time series. The scale agreement between these two technique solutions is estimated to be 1.05 ± 0.13 ppb at epoch 2005.0 and 0.049 ± 0.010?ppb/yr. The ITRF2008 orientation (at epoch 2005.0) and its rate are aligned to the ITRF2005 using 179 stations of high geodetic quality. An estimate of the origin components from ITRF2008 to ITRF2005 (both origins are defined by SLR) indicates differences at epoch 2005.0, namely: ?0.5, ?0.9 and ?4.7?mm along X, Y and Z-axis, respectively. The translation rate differences between the two frames are zero for Y and Z, while we observe an X-translation rate of 0.3?mm/yr. The estimated formal errors of these parameters are 0.2?mm and 0.2?mm/yr, respectively. The high level of origin agreement between ITRF2008 and ITRF2005 is an indication of an imprecise ITRF2000 origin that exhibits a Z-translation drift of 1.8?mm/yr with respect to ITRF2005. An evaluation of the ITRF2008 origin accuracy based on the level of its agreement with ITRF2005 is believed to be at the level of 1?cm over the time-span of the SLR observations. Considering the level of scale consistency between VLBI and SLR, the ITRF2008 scale accuracy is evaluated to be at the level of 1.2?ppb (8?mm at the equator) over the common time-span of the observations of both techniques. Although the performance of the ITRF2008 is demonstrated to be higher than ITRF2005, future ITRF improvement resides in improving the consistency between local ties in co-location sites and space geodesy estimates.     

基于我国全球连续监测与评估系统的ERP测定的精度分析

利用GPS观测资料解算地球自转参数,用全球均匀分布的22个IGS跟踪站（IGS05）的连续观测资料估计地球自转参数（ERP）,并与IERSC04（UTC0时）的结果相比较,二者相差很小,均在IERS的ERP估计精度范围之内。基于即将建成的COMPASS全球连续监测与评估系统跟踪站,选择其网的8个IGS跟踪站的资料进行了解算并进行了分析和比对。     

DORIS and the Determination of the Earth’s Polar Motion

DORIS (Détermination d’Orbite et Radiopositionnement Intégrés par Satellite) is a system used for precise orbit determination (POD) and ground-station positioning. It has been implemented on-board various satellites: the SPOT (Système pour l’Observation de la Terre) remote sensing satellites SPOT-2, SPOT-3, SPOT-4, SPOT-5, TOPEX/Poseidon and more recently on its successors Jason-1 and ENVISAT. DORIS is also a terrestrial positioning system that has found many applications in geophysics and geodesy; in particular, it contributes to the realization of the International Terrestrial Reference Frame, ITRF2000 and the forthcoming ITRF2005. Although not its primary objective, DORIS can bring information on Earth orientation monitoring, mainly polar motion and length of day (LOD) variations that complement other astrogeodetic techniques. In this paper, we have analyzed various recent polar motion solutions derived from independent analysis centers using different software packages and applying various analysis strategies. Comparisons of these solutions to the International Earth Rotation and Reference Systems Service (IERS) C04 solution are performed. Depending on the solutions, the accuracy of DORIS polar components are in the range of 0.5–1 mas corresponding to a few centimeters on the Earth’s surface. This is approximately ten times larger than results derived from GPS, which are typically 0.06 mas in both components. This does not allow DORIS results to be taken into account in the IERS–EOP combinations. A gain in the precision could come from technical improvements to the DORIS system, in addition to improvement of the orbit, tropospheric, ionospheric and Earth gravity field modeling.     

地球参考框架ITRF2014结果分析

ITRF2014是地球参考系的最新实现。该框架利用正弦函数估计负荷对台站位置的季节性效应,与ITRF2008相比,可以得到更稳定、精确的速度场;另外,ITRF2014引入了震后形变模型,可以更好地分析测站的非线性运动。本文通过分析发现:ITRF2014其原点相较于ITRF2008,其符合精度为3.5 mm;两种技术(VLBI和SLR)在2010.0历元确定的尺度因子不符值为1.18 ppb;同时,局部测量解与空间大地测量解解算的本地连接向量仍存在较大不符。     

VLBI terrestrial reference frame contributions to ITRF2008

In late 2008, the Product Center for the International Terrestrial Reference Frame (ITRF) of the International Earth Rotation and Reference Systems Service (IERS) issued a call for contributions to the next realization of the International Terrestrial Reference System, ITRF2008. The official contribution of the International VLBI Service for Geodesy and Astrometry (IVS) to ITRF2008 consists of session-wise datum-free normal equations of altogether 4,539 daily Very Long Baseline Interferometry (VLBI) sessions from 1979.7 to 2009.0 including data of 115 different VLBI sites. It is the result of a combination of individual series of session-wise datum-free normal equations provided by seven analysis centers (ACs) of the IVS. All series are completely reprocessed following homogeneous analysis options according to the IERS Conventions 2003 and IVS Analysis Conventions. Altogether, nine IVS ACs analyzed the full history of VLBI observations with four different software packages. Unfortunately, the contributions of two ACs, Institute of Applied Astronomy (IAA) and Geoscience Australia (AUS), had to be excluded from the combination process. This was mostly done because the IAA series exhibits a clear scale offset while the solution computed from normal equations contained in the AUS SINEX files yielded unreliable results. Based on the experience gathered since the combination efforts for ITRF2005, some discrepancies between the individual series were discovered and overcome. Thus, the consistency of the individual VLBI solutions has improved considerably. The agreement in terms of WRMS of the Terrestrial Reference Frame (TRF) horizontal components is 1 mm, of the height component 2 mm. Comparisons between ITRF2005 and the combined TRF solution for ITRF2008 yielded systematic height differences of up to 5 mm with a zonal signature. These differences can be related to a pole tide correction referenced to a zero mean pole used by four of five IVS ACs in the ITRF2005 contribution instead of a linear mean pole path as recommended in the IERS Conventions. Furthermore, these systematics are the reason for an offset in the scale of 0.4 ppb between the IVS’ contribution to ITRF2008 and ITRF2005. The Earth orientation parameters of seven series used as input for the IVS combined series are consistent to a huge amount with about 50 μas WRMS in polar motion and 3 μs in dUT1.     

The 2008 DGFI realization of the ITRS: DTRF2008

A new realization of the International Terrestrial System was computed at the ITRS Combination Centre at DGFI as a contribution to ITRF2008. The solution is labelled DTRF2008. In the same way as in the DGFI computation for ITRF2005 it is based on either normal equation systems or estimated parameters derived from VLBI, SLR, GPS and DORIS observations by weekly or session-wise processing. The parameter space of the ITRS realization comprises station positions and velocities and daily resolved Earth Orientation Parameters (EOP), whereby for the first time also nutation parameters are included. The advantage of starting from time series of input data is that the temporal behaviour of geophysical parameters can be investigated to decide whether the parameters can contribute to the datum realization of the ITRF. In the same way, a standardized analysis of station position time series can be performed to detect and remove discontinuities. The advantage of including EOP in the ITRS realization is twofold: (1) the combination of the coordinates of the terrestrial pole—estimated from all contributing techniques—links the technique networks in two components of the orientation, leading to an improvement of consistency of the Terrestrial Reference Frame (TRF) and (2) in their capacity as parameters common to all techniques, the terrestrial pole coordinates enhance the selection of local ties as they provide a measure for the consistency of the combined frame. The computation strategy of DGFI is based on the combination of normal equation systems while at the ITRS Combination Centre at IGN solutions are combined. The two independent ITRS realizations provide the possibility to assess the accuracy of ITRF by comparison of the two frames. The accuracy evaluation was done separately for the datum parameters (origin, orientation and scale) and the network geometry. The accuracy of the datum parameters, assessed from the comparison of DTRF2008 and ITRF2008, is between 2–5?mm and 0.1–0.8?mm/year depending on the technique. The network geometry (station positions and velocities) agrees within 3.2?mm and 1.0?mm/year. A comparison of DTRF2008 and ITRF2005 provides similar results for the datum parameters, but there are larger differences for the network geometry. The internal accuracy of DTRF2008—that means the level of conservation of datum information and network geometry within the combination—was derived from comparisons with the technique-only multi-year solutions. From this an internal accuracy of 0.32?mm for the VLBI up to 3.3?mm for the DORIS part of the network is found. The internal accuracy of velocities ranges from 0.05?mm/year for VLBI to 0.83?mm/year for DORIS. The internal consistency of DTRF2008 for orientation can be derived from the analysis of the terrestrial pole coordinates. It is estimated at 1.5–2.5?mm for the GPS, VLBI and SLR parts of the network. The consistency of these three and the DORIS network part is within 6.5?mm.     

IGS contribution to the ITRF

We examine the contribution of the International GNSS Service (IGS) to the International Terrestrial Reference Frame (ITRF) by evaluating the quality of the incorporated solutions as well as their major role in the ITRF formation. Starting with the ITRF2005, the ITRF is constructed with input data in the form of time series of station positions (weekly for satellite techniques and daily for VLBI) and daily Earth Orientation Parameters. Analysis of time series of station positions is a fundamental first step in the ITRF elaboration, allowing to assess not only the stations behavior, but also the frame parameters and in particular the physical ones, namely the origin and the scale. As it will be seen, given the poor number and distribution of SLR and VLBI co-location sites, the IGS GPS network plays a major role by connecting these two techniques together, given their relevance for the definition of the origin and the scale of the ITRF. Time series analysis of the IGS weekly combined and other individual Analysis Center solutions indicates an internal precision (or repeatability) <2 mm in the horizontal component and <5 mm in the vertical component. Analysis of three AC weekly solutions shows generally poor agreement in origin and scale, with some indication of better agreement when the IGS started to use the absolute model of antenna phase center variations after the GPS week 1400 (November 2006).     

利用卫星激光测距技术研究全球垂直地壳运动的参考基准

当前地震、地质灾害等现象的频繁发生,表明对于地壳运动尤其是垂直方向的研究是迫切而必需的,而研究垂直地壳运动的一个重要前提就是其参考基准的确立。针对目前国内外对于全球垂直地壳运动参考基准不统一的现状,分析了研究全球统一的垂直地壳运动参考基准的必要性,提出将ITRF2008框架下获取的SLR的垂向速度场作为全球垂直地壳运动参考基准的方案,论证了该方案的合理性,并利用GPS和VLBI的观测数据,求取了其速度场与SLR垂直基准之间的系统差,分析了产生不同大小且方向各异系统差的原因,实现了各速度场与全球垂直地壳运动参考基准的统一。由于该基准包含了较为丰富的地球物理信息,具有可发展、可提高以及直接相对于地球质心等特性,因此,我们推荐将ITRF2008框架下SLR技术获取的垂向速度场定义为全球垂直地壳运动的参考基准。     

Evaluation of the ITRF2008 GPS vertical velocities using satellite antenna z-offsets

We develop a method to evaluate the terrestrial reference frame (TRF) scale rate error using Global Positioning System (GPS) satellite antenna phase center offset (APCO) parameters and apply it to ITRF2008. We search for the TRF in which z-APCO parameters have the smallest drift. In order to provide realistic error bars for the z-APCO drifts, we pay attention to model periodic variations and auto-correlated noise processes in the z-APCO time series. We will show that the GPS scale rate with respect to a frame is, as a first approximation, proportional to the estimated mean z-APCO trend if that frame is used to constrain station positions. Thus, an ITRF2008 scale rate error between ?0.27 and ?0.06 mm/yr depending on the GPS analysis center can be estimated, which demonstrates the high quality of the newly constructed ITRF2008. We will also demonstrate that the traditional estimates of the GPS scale rate from 7-parameter similarity transformations are consistent with our newly derived GPS scale rates with respect to ITRF2008 within two sigmas. We find using International GNSS Service (IGS) products that the traditional approach is relevant for scale rate determination even if some of the z-APCO values supplied by the IGS were not simultaneously calibrated. As the scale rate is related to the accuracy of vertical velocities, our estimates supply a conservative evaluation that can be used for error budget computation.