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.     

Impact of loading displacements on SLR-derived parameters and on the consistency between GNSS and SLR results

Displacements of the Earth’s surface caused by tidal and non-tidal loading forces are relevant in high-precision space geodesy. Some of the corrections are recommended by the international scientific community to be applied at the observation level, e.g., ocean tidal loading (OTL) and atmospheric tidal loading (ATL). Non-tidal displacement corrections are in general recommended not to be applied in the products of the International Earth Rotation and Reference Systems Service, in particular atmospheric non-tidal loading (ANTL), oceanic and hydrological non-tidal corrections. We assess and compare the impact of OTL, ATL and ANTL on SLR-derived parameters by reprocessing 12 years of SLR data considering and ignoring individual corrections. We show that loading displacements have an influence not only on station long-term stability, but also on geocenter coordinates, Earth Rotation Parameters, and satellite orbits. Applying the loading corrections reduces the amplitudes of annual signals in the time series of geocenter and station coordinates. The general improvement of the SLR station 3D coordinate repeatability when applying OTL, ATL and ANTL corrections are 19.5 %, 0.2 % and 3.3 % respectively, w.r.t. the solutions without loading corrections. ANTL corrections play a crucial role in the combination of optical (SLR) and microwave (GNSS, VLBI, DORIS) space geodetic observation techniques, because of the so-called Blue-Sky effect: SLR measurements can be carried out only under cloudless sky conditions—typically during high air pressure conditions, when the Earth’s crust is deformed, whereas microwave observations are weather-independent. Thus, applying the loading corrections at the observation level improves SLR-derived products as well as the consistency with microwave-based results. We assess the Blue-Sky effect on SLR stations and the consistency improvement between GNSS and SLR solutions when ANTL corrections are included. The omission of ANTL corrections may lead to inconsistencies between SLR and GNSS solutions of up to 2.5 mm for inland stations. As a result, the estimated GNSS–SLR coordinate differences correspond better to the local ties at the co-located stations when applying ANTL corrections.     

Consistent realization of Celestial and Terrestrial Reference Frames

The Celestial Reference System (CRS) is currently realized only by Very Long Baseline Interferometry (VLBI) because it is the space geodetic technique that enables observations in that frame. In contrast, the Terrestrial Reference System (TRS) is realized by means of the combination of four space geodetic techniques: Global Navigation Satellite System (GNSS), VLBI, Satellite Laser Ranging (SLR), and Doppler Orbitography and Radiopositioning Integrated by Satellite. The Earth orientation parameters (EOP) are the link between the two types of systems, CRS and TRS. The EOP series of the International Earth Rotation and Reference Systems Service were combined of specifically selected series from various analysis centers. Other EOP series were generated by a simultaneous estimation together with the TRF while the CRF was fixed. Those computation approaches entail inherent inconsistencies between TRF, EOP, and CRF, also because the input data sets are different. A combined normal equation (NEQ) system, which consists of all the parameters, i.e., TRF, EOP, and CRF, would overcome such an inconsistency. In this paper, we simultaneously estimate TRF, EOP, and CRF from an inter-technique combined NEQ using the latest GNSS, VLBI, and SLR data (2005–2015). The results show that the selection of local ties is most critical to the TRF. The combination of pole coordinates is beneficial for the CRF, whereas the combination of \(\varDelta \hbox {UT1}\) results in clear rotations of the estimated CRF. However, the standard deviations of the EOP and the CRF improve by the inter-technique combination which indicates the benefits of a common estimation of all parameters. It became evident that the common determination of TRF, EOP, and CRF systematically influences future ICRF computations at the level of several \(\upmu \)as. Moreover, the CRF is influenced by up to \(50~\upmu \)as if the station coordinates and EOP are dominated by the satellite techniques.     

微型激光测距与甚长基线干涉测量的天线旋转中心的直接解法

精确确定SLR (Satellite Laser Rainging)和VLBI (Very Long Baseline Interferometry) 天线的旋转中心是并址站归心基线测量的关键问题。本文利用两类约束条件建立旋转中心与观测标志之间的直接关系,第1类约束是SLR或VLBI上的观测标志绕其旋转轴旋转形成1个由平面和球面相割得到的平面圆,第2类约束是SLR或VLBI的旋转中心与其垂直轴旋转圆心在同一铅垂线上,与其水平轴旋转圆心在同一水平面上。根据这两类约束条件建立相应的条件方程,利用标志点观测值直接解算旋转中心的坐标及其协方差阵。利用我国2个GNSS (Global Navigation Satellite System) 与SLR或VLBI并址站的实测数据,求解了基于本文直接解法的空间归心基线。结果表明,与已有分步解的差值优于1mm。     

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.     

GPS、SLR、VLBI联合处理站坐标

采用了GPS、VLBI、SLR三种技术的并置站坐标,计算了三种技术实现的参考框架的转换参数,联合处理得到并置站的坐标并与IERS公布的坐标进行了比较。     

相似变换模型在参考框架基准实现中的应用

详细探讨了相似变换在地球参考框架的基准定义实施中的相关模型及方法。利用四种空间大地测量技术（GPS、SLR、DORIS、VLBI）在2000．0历元的站坐标组文件,对ITRF2005的内符精度进行了实验评价。     

Precise station positions from VLBI observations to satellites: a simulation study

Very long baseline interferometry (VLBI) tracking of satellites is a topic of increasing interest for the establishment of space ties. This shall strengthen the connection of the various space geodetic techniques that contribute to the International Terrestrial Reference Frame. The concept of observing near-Earth satellites demands research on possible observing strategies. In this paper, we introduce this concept and discuss its possible benefits for improving future realizations of the International Terrestrial Reference System. Using simulated observations, we develop possible observing strategies that allow the determination of radio telescope positions in the satellite system on Earth with accuracies of a few millimeters up to 1–2 cm for weekly station coordinates. This is shown for satellites with orbital heights between 2,000 and 6,000 km, observed by dense regional as well as by global VLBI-networks. The number of observations, as mainly determined by the satellite orbit and the observation interval, is identified as the most critical parameter that affects the expected accuracies. For observations of global positioning system satellites, we propose the combination with classical VLBI to radio sources or a multi-satellite strategy. Both approaches allow station position repeatabilities of a few millimeters for weekly solutions.     

顾及测量方式的地球自转参数差异性分析

随着甚长基线干涉测量(VLBI)、卫星激光测距(SLR)、激光测月(LLR)、全球卫星导航系统(GNSS)、多里斯系统(DORIS)等多种空间大地测量手段的使用,地球自转参数(ERP)的测量精度不断提高,为航天器导航、深空探测等诸多领域提供了高精度的国际天球参考系(ICRS)和国际地表参考系统(ITRS)之间的转换参数. 以国际地球自转与参考系服务发布的C04序列为基础序列,选取500天ERP序列,分析不同测量手段得到的ERP数据的误差分布情况,为研究利用不同数据之间的一致性进行精度检核的可行性及精度水平提供数据基础,同时也为ERP预报提供更多的数据选择.      

The IGS-combined station coordinates, earth rotation parameters and apparent geocenter

The International GNSS Service (IGS) routinely generates a number of weekly, daily and sub-daily products. Station coordinates and velocities, earth rotation parameters (ERPs) and apparent geocenter are among these products generated weekly by the IGS Reference Frame Coordinator. They have been determined since 1999 by combining independent estimates from at least seven IGS Analysis Centers (ACs). Two Global Network Associate Analysis Centers (GNAACs) also provide independent combinations using the same AC weekly solutions and they are currently used to quality control the IGS combination. The combined solutions are aligned to an IGS realization (IGS05) of the ITRF2005 using a carefully selected set of the IGS Reference Frame (RF) stations (nominally 132). During the combination process, the contributing solutions are compared and outliers are removed to ensure a high level of consistency of the estimated parameters. The ACs and the weekly combined solution are consistent at the 1–2 and 3–4 mm levels for the horizontal and vertical components. Similarly, the excess Length of Day (LOD), the pole positions and pole rates are consistent at the 10μs, 0.03–0.05 mas and 0.10–0.20 mas/day levels, respectively. The consistency of the apparent geocenter estimate is about 5 mm in the X and Y components and 10 mm in the Z component. Comparison of the IGS-combined ERP estimates with the IERS Bulletin A suggests a small bias of the order of ?0.04 mas and + 0.05 mas (both ±0.05 mas) in the x and y components.     

Combined Earth orientation parameters based on homogeneous and continuous VLBI and GPS data

The CONT02 campaign is of great interest for studies combining very long baseline interferometry (VLBI) with other space-geodetic techniques, because of the continuously available VLBI observations over 2 weeks in October 2002 from a homogeneous network. Especially, the combination with the Global Positioning System (GPS) offers a broad spectrum of common parameters. We combined station coordinates, Earth orientation parameters (EOPs) and troposphere parameters consistently in one solution using technique- specific datum-free normal equation systems. In this paper, we focus on the analyses concerning the EOPs, whereas the comparison and combination of the troposphere parameters and station coordinates are covered in a companion paper in Journal of Geodesy. In order to demonstrate the potential of the VLBI and GPS space-geodetic techniques, we chose a sub-daily resolution for polar motion (PM) and universal time (UT). A consequence of this solution set-up is the presence of a one-to-one correlation between the nutation angles and a retrograde diurnal signal in PM. The Bernese GPS Software used for the combination provides a constraining approach to handle this singularity. Simulation studies involving both nutation offsets and rates helped to get a deeper understanding of this singularity. With a rigorous combination of UT1–UTC and length of day (LOD) from VLBI and GPS, we showed that such a combination works very well and does not suffer from the systematic effects present in the GPS-derived LOD values. By means of wavelet analyses and the formal errors of the estimates, we explain this important result. The same holds for the combination of nutation offsets and rates. The local geodetic ties between GPS and VLBI antennas play an essential role within the inter-technique combination. Several studies already revealed non-negligible discrepancies between the terrestrial measurements and the space-geodetic solutions. We demonstrate to what extent these discrepancies propagate into the combined EOP solution.     

Results of the VLBI experiments conducted with Syowa Station, Antarctica

The first successful geodetic Very Long Baseline Interferometry (VLBI) observations to Antarctica were made on baselines from Syowa Station (Antarctica) to Tidbinbilla (Australia) and to Kashima (Japan) in January 1990. Regular geodetic experiments started in 1998 with the installation of a permanent VLBI terminal at Syowa Station. These observations are conducted at the standard geodetic VLBI frequencies of 2.3 and 8.4 GHz, S- and X-Bands. In the first year, the 11-m multipurpose antenna at Syowa Station observed together with the 26-m radio telescope of the University of Tasmania in Australia and the 26-m radio telescope of the Hartebeesthoek Radio Astronomy Observatory in South Africa. From 1999, the experiments were expanded to also include the O’Higgins Station in Antarctica, Fortaleza in Brazil and Kokee on Hawaii. From 1999 until the end of 2003, 25 observing sessions have been reduced and analyzed using the CALC/SOLVE geodetic VLBI data reduction package. The results show that the horizontal baseline of Syowa-Hobart is increasing at the rate of 57.0±1.9 mm/year. The baseline Syowa-Hartebeesthoek is also increasing, but at the lower rate of 9.8±1.9 mm/year. The VLBI result of 2.0±3.1 mm/year and the GPS result of −1.9±0.7 mm/year for the Syowa-O’Higgins horizontal baseline support the hypothesis of one rigid Antarctic plate without intra-plate deformation, which is consistent with the NNR-NUVEL-1A global plate motion model. The location of the Euler pole of the Antarctic plate by VLBI is estimated as 59.7°S and 62.6°E with a rotation rate of 0.190 deg/Myr, while that by GPS in our study is estimated as 60.6°S and 42.2°E with a rotation rate of 0.221 deg/Myr. These pole positions are slightly different to that implied by the NNR-NUVEL-1A model of 63.0°S and 64.2°E with a rotation rate of 0.238 deg/Myr. VLBI observations over a longer time span may resolve small discrepancy of current plate motion from the NNR-NUVEL-1A model. The consistency of the VLBI coordinates with the GPS coordinates at Syowa Station, after correction for the local tie vector components between the two reference markers, is also discussed.     

Geodetic VLBI with an artificial radio source on the Moon: a simulation study

We perform extensive simulations in order to assess the accuracy with which the position of a radio transmitter on the surface of the Moon can be determined by geodetic VLBI. We study how the quality and quantity of geodetic VLBI observations influence these position estimates and investigate how observations of such near-field objects affect classical geodetic parameters like VLBI station coordinates and Earth rotation parameters. Our studies are based on today’s global geodetic VLBI schedules as well as on those designed for the next-generation geodetic VLBI system. We use Monte Carlo simulations including realistic stochastic models of troposphere, station clocks, and observational noise. Our results indicate that it is possible to position a radio transmitter on the Moon using today’s geodetic VLBI with a two-dimensional horizontal accuracy of better than one meter. Moreover, we show that the next-generation geodetic VLBI has the potential to improve the two-dimensional accuracy to better than 5 cm. Thus, our results lay the base for novel observing concepts to improve both lunar research and geodetic VLBI.     

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

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

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.     

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.     

Evaluation of co-location ties relating the VLBI and GPS reference frames

We have compared the VLBI and GPS terrestrial reference frames, realized using 5 years of time-series observations of station positions and polar motion, with surveyed co-location tie vectors for 25 sites. The goal was to assess the overall quality of the ties and to determine whether a subset of co-location sites might be found with VLBI–GPS ties that are self-consistent within a few millimeters. Our procedure was designed to guard against internal distortion of the two space-geodetic networks and takes advantage of the reduction in tie information needed with the time-series combination method by using the very strong contribution due to co-location of the daily pole of rotation. The general quality of the available ties is somewhat discouraging in that most have residuals, compared to the space-geodetic frames, at the level of 1–2 cm. However, by a careful selection process, we have identified a subset of nine local VLBI–GPS ties that are consistent with each other and with space geodesy to better than 4 mm (RMS) in each component. While certainly promising, it is not possible to confidently assess the reliability of this particular subset without new information to verify the absolute accuracy of at least a few of the highest-quality ties. Particular care must be taken to demonstrate that possible systematic errors within the VLBI and GPS systems have been properly accounted for. A minimum of two (preferably three or four) ties must be measured with accuracies of 1 mm or better in each component, including any potential systematic effects. If this can be done, then the VLBI and GPS frames can be globally aligned to less than 1 mm in each Helmert component using our subset of nine ties. In any case, the X and Y rotations are better determined, to about 0.5 mm, due to the contribution of co-located polar motion.     

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.     

Adding source positions to the IVS combination—First results

The consistent estimation of terrestrial reference frames (TRF), celestial reference frames (CRF) and Earth orientation parameters (EOP) is still an open subject and offers a large field of investigations. Until now, source positions resulting from Very Long Baseline Interferometry (VLBI) observations are not routinely combined on the level of normal equations in the same way as it is a common process for station coordinates and EOPs. The combination of source positions based on VLBI observations is now integrated in the IVS combination process. We present the studies carried out to evaluate the benefit of the combination compared to individual solutions. On the level of source time series, improved statistics regarding weighted root mean square have been found for the combination in comparison with the individual contributions. In total, 67 stations and 907 sources (including 291 ICRF2 defining sources) are included in the consistently generated CRF and TRF covering 30 years of VLBI contributions. The rotation angles \(A_1\), \(A_2\) and \(A_3\) relative to ICRF2 are ?12.7, 51.7 and 1.8 \({\upmu }\) as, the drifts \(D_\alpha \) and \(D_\delta \) are ?67.2 and 19.1 \(\upmu \) as/rad and the bias \(B_\delta \) is 26.1 \(\upmu \) as. The comparison of the TRF solution with the IVS routinely combined quarterly TRF solution shows no significant impact on the TRF, when the CRF is estimated consistently with the TRF. The root mean square value of the post-fit station coordinate residuals is 0.9 cm.     

Is there utility in rigorous combinations of VLBI and GPS Earth orientation parameters?

Combinations of station coordinates and velocities from independent space-geodetic techniques have long been the standard method to realize robust global terrestrial reference frames (TRFs). In principle, the particular strengths of one observing method can compensate for weaknesses in others if the combination is properly constructed, suitable weights are found, and accurate co-location ties are available. More recently, the methodology has been extended to combine time-series of results at the normal equation level. This allows Earth orientation parameters (EOPs) to be included and aligned in a fully consistent way with the TRF. While the utility of such multi-technique combinations is generally recognized for the reference frame, the benefits for the EOPs are yet to be quantitatively assessed. In this contribution, which is a sequel to a recent paper on co-location ties (Ray and Altamimi in J Geod 79(4–5): 189–195, 2005), we have studied test combinations of very long baseline interferometry (VLBI) and Global Positioning System (GPS) time-series solutions to evaluate the effects on combined EOP measurements compared with geophysical excitations. One expects any effect to be small, considering that GPS dominates the polar motion estimates due to its relatively dense and uniform global network coverage, high precision, continuous daily sampling, and homogeneity, while VLBI alone observes UT1-UTC. Presently, although clearly desirable, we see no practical method to rigorously include the GPS estimates of length-of-day variations due to significant time-varying biases. Nevertheless, our results, which are the first of this type, indicate that more accurate polar motion from GPS contributes to improved UT1-UTC results from VLBI. The situation with combined polar motion is more complex. The VLBI data contribute directly only very slightly, if at all, with an impact that is probably affected by the weakness of the current VLBI networks (small size and sparseness) and the quality of local ties relating the VLBI and GPS frames. Instead, the VLBI polar motion information is used primarily in rotationally aligning the VLBI and GPS frames, thereby reducing the dependence on co-location tie information. Further research is needed to determine an optimal VLBI-GPS combination strategy that yields the highest quality EOP estimates. Improved local ties (including internal systematic effects within the techniques) will be critically important in such an effort.