ITRF2014-GPS解数据质量分析

现搜集了最新的ITRF2014-GPS解数据,并利用移动回归模型与快速傅里叶变换进行了数据处理,将处理结果与ITRF2008-GPS解处理结果进行比较发现,ITRF2014-GPS解数据质量比ITRF2008有明显提高,其中主要非线性周期规律已经较好地被模型化改正,主要包含周年规律和半周年规律.     

ITRF 2008与ITRF 2014的差异对GNSS数据处理的影响

针对全球导航卫星系统(GNSS)数据处理过程中旧的ITRF 2008参考框架现势性不足及新的ITRF2014框架在数据的数量与质量、参数模型、测站的分布合理性上均有提高等状况,该文以陆态网的最近两年的观测数据为例,对比分析了ITRF2008和ITRF2014框架下各测站的坐标、基线长度、水平速度场的差异,以期为当前高精度GNSS数据处理提供参考。实验表明:两个框架下的成果经基准转换后,测站在X、Y、Z方向的差异均为毫米级;基线差异平均在1 mm以内;水平速度场差值的最大值为5.75(mm·a~(-1)),最小值为-4.88(mm·a~(-1)),平均值为-0.45(mm·a~(-1)),方向上差值的平均值为0.02rad。目前两个框架的差异对一般工程应用基本上可以忽略,但对地震监测的陆态网来说,则必须考虑。     

ITRF2008与CGCS2000坐标系的转换

由于当前精密星历所对应解算的ITRF框架坐标为ITRF2008参考框架,而在1∶10 000基础测绘生产项目要求提供CGCS2000坐标系成果,论述了ITRF2008到CGCS2000间的框架转换的方法及转换后精度分析,并重点分析了转换的关键性问题。     

ITRF2008框架下VLBI与GPS确定地心坐标的精度分析

为了解ITRF2008框架下VLBI和GPS两种空间技术确定地心坐标的真正实现精度,在并置站上对VLBI和GPS两种空间技术测定的地心坐标进行了比较,经过偏心改正和七参数转换之后,得到两种空间技术地心坐标不符值的加权中误差,其可以认为是这两种空间技术的真正实现精度,经比较分析这两种地心坐标三个坐标轴方向分量的外符精度都在10mm之内,说明VLBI和GPS确定的地心坐标精度已达到毫米级。     

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

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

坐标系统任意框架和历元转换研究

研究国际陆地参考框架(ITRF)坐标转换问题,分析坐标框架的转换模型,构建一种实用的历元转换方法.研究结果表明在相同历元下,不同框架之间的坐标差异较小,不同活动块体下ITRF2008和ITRF2014在2020.00历元时各分量坐标差小于3 mm.受不同区域地壳活动程度影响,相同框架下不同历元之间的坐标差异明显,基于建立的历元转换方法,实际测试精度可达1 cm左右,能够在工程应用中快速获得CGCS2000坐标.     

陆态网络基准网数据处理策略分析

从基准网测站分区和整网平差策略方面对基准网数据处理进行了详细分析,并进行了强约束法和相似变换法的比较实验。结果表明,强约束法对框架精度要求较高,约束点精度对结果有直接影响;而相似变换法保留了区域框架的本质,避免了网形结构的扭曲。由于CGCS2000框架是基于现势性较差的ITRF97的,结果不理想,因此数据处理过程中宜采用现势性较好的ITRF2008或ITRF2014站点联合解算。     

ITRF2014参考框架的实现与改进

ITRF2014通过对多种空间大地测量技术解的联合处理,在ITRF2014建立过程中首次对非线性运动建模,包括季节性变化的估计和震后形变(PSD)模型的应用。针对基准定义、输入数据和数据处理策略等方面介绍ITRF2014实现的基本情况,并对之前版本进行优化改进。     

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和日长数据。     

现今板块运动模型的发展及其比较

从地质数据的不断完善到空间大地测量数据的广泛应用,现今板块运动模型的发展有了质的飞跃.评述了其现状与发展趋势;利用ITRF2005的速度场建立了一个完全基于现代空间大地测量实测结果的现今全球板块运动模型ITRF2005VEL,并与地质模型NNR-NUVEL1A进行了比较.结果表明,从整体上看ITRF2005VEL与NNR-NUVEL1A模型有较好的一致性,且较以往的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.     

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.     

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.     

Defining a DORIS core network for Jason-1 precise orbit determination based on ITRF2000: methods and realization

In view of the future adoption of the new precise orbit determination (POD) standards for the TOPEX/Poseidon and Jason-1 satellites, we propose a method to evaluate terrestrial reference frames for POD. We applied this method to the ITRF2000 realization of the DORIS network using local geodetic ties, plate motion models, the recent DORIS IGN04D02 cumulative solution and DORIS weekly time-series of coordinates. We propose to adopt a selection of the ITRF2000 realization based on specific criteria that we define here, and to extend it with ground stations for which we propose new coordinates and velocities. Only 13 out of 131 stations were considered to be inappropriate for POD activities. The result is a robust and well-distributed DORIS core network of 118 stations (DPOD2000) suitable for POD during the 1993–2008 period considered here.     

IGS08: the IGS realization of ITRF2008

On April 17, 2011, the International GNSS Service (IGS) stopped using the IGS05 reference frame and adopted a new one, called IGS08, as the basis of its products. The latter was derived from the latest release of the International Terrestrial Reference Frame (ITRF2008). However, the simultaneous adoption of a new set of antenna phase center calibrations by the IGS required slight adaptations of ITRF2008 positions for 65 of the 232 IGS08 stations. The impact of the switch from IGS05 to IGS08 on GNSS station coordinates was twofold: in addition to a global transformation due to the frame change from ITRF2005 to ITRF2008, many station coordinates underwent small shifts due to antenna calibration updates, which need to be accounted for in any comparison or alignment of an IGS05-consistent solution to IGS08. Because the heterogeneous distribution of the IGS08 network makes it sub-optimal for the alignment of global frames, a smaller well-distributed sub-network was additionally designed and designated as the IGS08 core network. Only 2?months after their implementation, both the full IGS08 network and the IGS08 core network already strongly suffer from the loss of many reference stations. To avoid a future crisis situation, updates of IGS08 will certainly have to be considered before the next ITRF release.     

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.