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: 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.     

卫星大地测量学的研究现状及发展趋势

从全球国际地球参考框架（International Terrestrial Reference Frame，ITRF）的建立、维护与发展，卫星测高、卫星重力等的发展及应用，全球卫星导航系统（Global Navigation Satellite System，GNSS）、卫星激光测距（Satellite Laser Ranging，SLR）、甚长基线干涉测量（very Long Baseline Interferometry，VLBI）、卫星多普勒定轨定位（Doppler Orbitography by Radiopositioning Integrated on Satellite，DORIS）的融合应用，海洋测绘和室内定位的发展等几个方面综述了大地测量学及卫星导航定位技术的最新进展，并提出中国2000国家大地坐标系与自主卫星导航系统的主要应用及发展目标。     

DORIS Contribution to ITRF2005

We examine the contribution of the Doppler Orbit determination and Radiopositioning Integrated by Satellite (DORIS) technique to the International Terrestrial Reference Frame (ITRF2005) by evaluating the quality of the submitted solutions as well as that of the frame parameters, especially the origin and the scale. Unlike the previous versions of the ITRF, ITRF2005 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 (EOPs), including full variance–covariance information. Analysis of the DORIS station positions’ time-series indicates an internal precision reaching 15 mm or better, at a weekly sampling. A cumulative solution using 12 years of weekly time-series was obtained and compared to a similar International GNSS Service (IGS) GPS solution (at 37 co-located sites) yielding a weighted root mean scatter (WRMS) of the order of 8 mm in position (at the epoch of minimum variance) and about 2.5 mm/year in velocity. The quality of this cumulative solution resulting from the combination of two individual DORIS solutions is better than any individual solution. A quality assessment of polar motion embedded in the contributed DORIS solutions is performed by comparison with the results of other space-geodetic techniques and in particular GPS. The inferred WRMS of polar motion varies significantly from one DORIS solution to another and is between 0.5 and 2 mas, depending on the strategy used and in particular estimating or not polar motion rate by the analysis centers. This particular aspect certainly needs more investigation by the DORIS Analysis Centers.     

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.     

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.     

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.     

Strategies to mitigate aliasing of loading signals while estimating GPS frame parameters

Although GNSS techniques are theoretically sensitive to the Earth center of mass, it is often preferable to remove intrinsic origin and scale information from the estimated station positions since they are known to be affected by systematic errors. This is usually done by estimating the parameters of a linearized similarity transformation which relates the quasi-instantaneous frames to a long-term frame such as the International Terrestrial Reference Frame (ITRF). It is well known that non-linear station motions can partially alias into these parameters. We discuss in this paper some procedures that may allow reducing these aliasing effects in the case of the GPS techniques. The options include the use of well-distributed sub-networks for the frame transformation estimation, the use of site loading corrections, a modification of the stochastic model by downweighting heights, or the joint estimation of the low degrees of the deformation field. We confirm that the standard approach consisting of estimating the transformation over the whole network is particularly harmful for the loading signals if the network is not well distributed. Downweighting the height component, using a uniform sub-network, or estimating the deformation field perform similarly in drastically reducing the amplitude of the aliasing effect. The application of these methods to reprocessed GPS terrestrial frames permits an assessment of the level of agreement between GPS and our loading model, which is found to be about 1.5 mm WRMS in height and 0.8 mm WRMS in the horizontal at the annual frequency. Aliased loading signals are not the main source of discrepancies between loading displacement models and GPS position time series.     

基于 EEMD 和 Husrt 指数的 GNSS 基准站的垂向速率估计

连续的全球卫星导航系统(GNSS)基准站的坐标时间序列中包含了复杂的噪声信号、非构造形变及其它因素的影响,尤其在垂直方向,对GNSS基准站在国际地球参考框架(ITRF)下的运动速率估计产生了较大的干扰。为进一步提高速率精度,文中采用整体模态分解(EEMD)方法对GNSS基准站的垂向观测时间序列进行分解,并根据各种信号的Hurst值进行分类及重构为噪声信号、季节性信号和长期趋势信号,采用最小二乘方法拟合长期趋势信号得到垂向速率。通过对中国大陆构造环境监测网络(CMONOC)的GNSS台站从2001—2013年近13a的垂向坐标时间序列的实例分析,采用基于EEMD和Husrt指数的最小二乘法能够准确地估计GNSS基准站的垂向速率。     

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

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

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.     

Combining the orbits of the IGS Analysis Centers

Currently seven Analysis Centers of the International GPS Service for Geodynamics (IGS) are producing daily precise orbits and the corresponding Earth Orientation Parameters (EOP). These individual products are available at several IGS Data Centers (e.g. CDDIS, IGN, SIO, etc.). During 1993 no official IGS orbits were produced, but the routine orbit comparisons by IGS indicated that, after small orientation and scale alignments, the orbit consistency was approaching the 20 cm level (a coordinate RMS), and that some orbit combination should be possible and feasible. An IGS combined orbit could provide a precise and efficient extension of the IERS Terrestrial Reference Frame (ITRF). Another advantage of such a combined orbit would be reliability and precision.Two schemes of orbit combinations are considered here: (a) the first method consists of a weighted averaging process of the earth-fixed satellite positions as produced by the individual Centers; (b) the second method uses the individual IGS orbit files as pseudo-observations in an orbit determination process, where in addition to the initial conditions, different parameter sets may be estimated. Both orbit combination methods have been tested on the January 1993 orbit data sets (GPS weeks 680 and 681) with an impressive agreement at the 5 cm level (coordinate RMS). The quality of the combined orbits is checked by processing a set of continental baselines in two different regions of the globe using different processing softwares. Both types of combined orbits gave similar baseline repeatability of a few ppb in both regions which compared favorably to the best individual orbits in the region.     

地球参考框架ITRF2014结果分析

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

Anomalous harmonics in the spectra of GPS position estimates

Prior studies of the power spectra of GPS position time series have found pervasive seasonal signals against a power-law background of flicker noise plus white noise. Dong et al. (2002) estimated that less than half the observed GPS seasonal power can be explained by redistributions of geophysical fluid mass loads. Much of the residual variation is probably caused by unidentified GPS technique errors and analysis artifacts. Among possible mechanisms, Penna and Stewart (2003) have shown how unmodeled analysis errors at tidal frequencies (near 12- and 24-hour periods) can be aliased to longer periods very efficiently. Signals near fortnightly, semiannual, and annual periods are expected to be most seriously affected. We have examined spectra for the 167 sites of the International GNSS (Global Navigation Satellite Systems) Service (IGS) network having more than 200 weekly measurements during 1996.0–2006.0. The non-linear residuals of the weekly IGS solutions that were included in ITRF2005, the latest version of the International Terrestrial Reference Frame (ITRF), have been used. To improve the detection of common-mode signals, the normalized spectra of all sites have been stacked, then boxcar smoothed for each local north (N), east (E), and height (H) component. The stacked, smoothed spectra are very similar for all three components. Peaks are evident at harmonics of about 1 cycle per year (cpy) up to at least 6 cpy, but the peaks are not all at strictly 1.0 cpy intervals. Based on the 6th harmonic of the N spectrum, which is among the sharpest and largest, and assuming a linear overtone model, then a common fundamental of 1.040 ± 0.008 cpy can explain all peaks well, together with the expected annual and semiannual signals. A flicker noise power-law continuum describes the background spectrum down to periods of a few months, after which the residuals become whiter. Similar sub-seasonal tones are not apparent in the residuals of available satellite laser ranging (SLR) and very long baseline interferometry (VLBI) sites, which are both an order of magnitude less numerous and dominated by white noise. There is weak evidence for a few isolated peaks near 1 cpy harmonics in the spectra of geophysical loadings, but these are much noisier than for GPS positions. Alternative explanations related to the GPS technique are suggested by the close coincidence of the period of the 1.040 cpy frequency, about 351.2 days, to the “GPS year”; i.e., the interval required for the constellation to repeat its inertial orientation with respect to the sun. This could indicate that the harmonics are a type of systematic error related to the satellite orbits. Mechanisms could involve orbit modeling defects or aliasing of site-dependent positioning biases modulated by the varying satellite geometry.     

The contribution of Very Long Baseline Interferometry to ITRF2005

The contribution of the International VLBI Service for Geodesy and Astrometry (IVS) to the ITRF2005 (International Terrestrial Reference Frame 2005) has been computed by the IVS Analysis Coordinator’s office at the Geodetic Institute of the University of Bonn, Germany. For this purpose the IVS Analysis Centres (ACs) provided datum-free normal equation matrices in Solution INdependent EXchange (SINEX) format for each 24 h observing session to be combined on a session-by-session basis by a stacking procedure. In this process, common sets of parameters, transformed to identical reference epochs and a prioris, and especially representative relative weights have been taken into account for each session. In order to assess the quality of the combined IVS files, Earth orientation parameters (EOPs) and scaling factors have been derived from the combined normal equation matrices. The agreement of the EOPs of the combined normal equation matrices with those of the individual ACs in terms of weighted root mean square (WRMS) is in the range of 50–60 μas for the two polar motion components and about 3 μs for UT1−UTC. External comparisons with International GNSS Serive (IGS) polar motion components is at the level of 130–170 μas and 21 μs/day for length of day (LOD). The scale of the terrestrial reference frame realized through the IVS SINEX files agrees with ITRF2000 at the level of 0.2 ppb.     

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.     

Transformation between SLR/VLBI and WGS-84 reference frames

In geodetic and geophysical applications of GPS, it is important to realize the ephemerides of the GPS satellites and the coordinates of station positions in a consistent reference system. At present, more than one reference system is being used by various GPS users depending on their specific applications. The WGS-84 and various reference frames based on satellite laser ranging (SLR), very long baseline interferometry (VLBI), or a combination of SLR and VLBI are the most commonly used in high precision geophysical applications. The WGS-84 is widely used in applications which rely on the GPS broadcast ephemeris. Station coordinates estimated in one system may have to be transformed to another for further use or for evaluation/comparison purposes. This paper presents a seven-parameter transformation between the WGS-84 and SLR/VLBI reference frames. The GPS double-differenced phase measurements for two consecutive weeks from a set of five Defense Mapping Agency (DMA) sites (defined in the WGS-84 frame) and from an augmented set of fifteen CIGNET sites (defined in the SLR/VLBI frame) were processed in a least squares estimation scheme to determine station coordinates, from which the transformation parameters were determined. A scale difference of about 0.2 ppm and an orientation difference in longitude of about 31 milliarcseconds were found to be the only parameters of significance between the adopted SLR/VLBI and the WGS-84 frames. Transformation between WGS-84 and the ITRF90 is also included, in which the scale difference is the same as before but the longitude rotation is about 16 mas.     

Surveying co-located space-geodetic instruments for ITRF computation

A new and comprehensive method is presented that can be used for estimating eccentricity vectors between global positioning system (GPS) antennas, doppler orbitography and radiopositioning integrated by satellites (DORIS) antennas, azimuth-elevation (AZ-EL) very long baseline interferometry (VLBI) telescopes, and satellite laser ranging (SLR) and lunar laser ranging (LLR) telescopes. The problem of reference point (RP) definition for these space-geodetic instruments is addressed and computed using terrestrial triangulation and electronic distance measurement (EDM) trilateration. The practical ground operations, the surveying approach and the terrestrial data processing are briefly illustrated, and the post-processing procedure is discussed. It is a geometrically based analytical approach that allows computation of RPs along with a rigorous statistical treatment of measurements. The tight connection between the geometrical model and the surveying procedure is emphasized. The computation of the eccentricity vector and the associated variance–covariance matrix between an AZ-EL VLBI telescope (with or without intersecting axes) and a GPS choke ring antenna is concentrated upon, since these are fundamental for computing the International Terrestrial Reference Frame (ITRF). An extension to RP computation and eccentricity vectors involving DORIS, SLR and LLR techniques is also presented. Numerical examples of the quality that can be reached using the authors approach are given. Working data sets were acquired in the years 2001 and 2002 at the radioastronomical observatory of Medicina (Italy), and have been used to estimate two VLBI-GPS eccentricity vectors and the corresponding SINEX files.