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.     

Assessment of the possible contribution of space ties on-board GNSS satellites to the terrestrial reference frame

The realization of the international terrestrial reference frame (ITRF) is currently based on the data provided by four space geodetic techniques. The accuracy of the different technique-dependent materializations of the frame physical parameters (origin and scale) varies according to the nature of the relevant observables and to the impact of technique-specific errors. A reliable computation of the ITRF requires combining the different inputs, so that the strengths of each technique can compensate for the weaknesses of the others. This combination, however, can only be performed providing some additional information which allows tying together the independent technique networks. At present, the links used for that purpose are topometric surveys (local/terrestrial ties) available at ITRF sites hosting instruments of different techniques. In principle, a possible alternative could be offered by spacecrafts accommodating the positioning payloads of multiple geodetic techniques realizing their co-location in orbit (space ties). In this paper, the GNSS–SLR space ties on-board GPS and GLONASS satellites are thoroughly examined in the framework of global reference frame computations. The investigation focuses on the quality of the realized physical frame parameters. According to the achieved results, the space ties on-board GNSS satellites cannot, at present, substitute terrestrial ties in the computation of the ITRF. The study is completed by a series of synthetic simulations investigating the impact that substantial improvements in the volume and quality of SLR observations to GNSS satellites would have on the precision of the GNSS frame parameters.     

New global positioning system reference station in Brazil

Co-located very long baseline interferometry (VLBI) and global positioning system (GPS) reference stations were installed near Fortaleza, Brazil, in 1993. Both have been important in the realization and maintenance of the International Terrestrial Reference Frame. A new-generation GPS system was installed in 2005 to replace the original station. Experience gained in the prior 12 years was used to improve the design of the GPS antenna mount. Preliminary indications are greatly improved data quality from the new station. Simultaneous observations from the nearly half-year of overlapping operation have been used to determine the local tie between the new and old GPS reference points to about 1 mm accuracy. This can be used to update the 1993 survey tie between the original GPS and the VLBI points, although there are questions about the accuracy of that measurement based on a comparison with space geodetic data. A test of removing the conical radome over the old GPS antenna indicates that it has biased the station height by about 16 mm downward, which probably accounts for most of the previous survey discrepancy.     

Effect of the selection of reference radio sources on geodetic estimates from VLBI observations

This paper evaluates the effect of the accuracy of reference radio sources on the daily estimates of station positions, nutation angle offsets, and the estimated site coordinates determined by very long baseline interferometry (VLBI), which are used for the realization of the international terrestrial reference frame (ITRF). Five global VLBI solutions, based on VLBI data collected between 1979 and 2006, are compared. The reference solution comprises all observed radio sources, which are treated as global parameters. Four other solutions, comprising different sub-sets of radio sources, were computed. The daily station positions for all VLBI sites and the corrections to the nutation offset angles were estimated for these five solutions. The solution statistics are mainly affected by the positional instabilities of reference radio sources, whereas the instabilities of geodetic and astrometric time-series are caused by an insufficient number of observed reference radio sources. A mean offset of the three positional components (Up, North, East) between any two solutions was calculated for each VLBI site. From a comparison of the geodetic results, no significant discrepancies between the respective geodetic solutions for all VLBI sites in the Northern Hemisphere were found. In contrast, the Southern Hemisphere sites were more sensitive to the selected set of reference radio sources. The largest estimated mean offset of the vertical component between two solutions for the Australian VLBI site at Hobart was 4 mm. In the worst case (if a weak VLBI network observed a limited number of reference radio sources) the daily offsets of the estimated height component at Hobart exceeded 100 mm. The exclusion of the extended radio sources from the list of reference sources improved the solution statistics and made the geodetic and astrometric time-series more consistent. The problem with the large Hobart height component offset is magnified by a comparatively small number of observations due to the low slewing rate of the VLBI dish (1°/ s). Unless a minimum of 200 scans are performed per 24-h VLBI experiment, the daily vertical positions at Hobart do not achieve 10 mm accuracy. Improving the slew rate at Hobart and/or having an increased number of new sites in the Southern Hemisphere is essential for further improvement of geodetic VLBI results for Southern Hemisphere sites.     

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.     

组合GPS和VLBI数据建立板块运动模型

利用GPS和VLBI的组合数据,解算了北美、欧亚、太平洋等12板块之间的相对运动欧拉矢量,得到了实测的板块运动模型GVM1。与地学模型(NUVEL1A)的比较指出,GVM1大体上与地学模型一致;EURANOAM的极位置与NUVEL1A的相应极比较接近,旋转速率略微偏大;澳大利亚板块在最近几年内是稳定的;太平洋板块与其他板块对的极位置与地学模型较为接近,这表明多种技术的组合数据提高了板块运动模型建立的精确性和可靠性。     

Ionospheric calibration of single frequency VLBI and GPS observations using dual GPS data

Summary The ionospheric effect is one of the main sources of error in Very Long Baseline Interferometry (VLBI) and Global Positioning System (GPS) high precision geodesy. Although the use of two frequencies allows the estimation of this effect, in some cases dual observations are not possible due to the available equipment or the type of observation. This paper presents the ionospheric calibration of single frequency VLBI and GPS observations based on the ionospheric electron content estimated from dual frequency GPS data. The ionospheric delays obtained with this procedure and the VLBI baseline length results have been compared with those obtained with dual frequency data. For the European geodetic VLBI baselines, both solutions agree at the 3–5 parts in 10^–9 level. The noise introduced by the GPS-based calibration is in the order of 3 cm for the VLBI observables and of 10 cm for the GPS observables.     

Height bias and scale effect induced by antenna gravitational deformations in geodetic VLBI data analysis

The impact of signal path variations (SPVs) caused by antenna gravitational deformations on geodetic very long baseline interferometry (VLBI) results is evaluated for the first time. Elevation-dependent models of SPV for Medicina and Noto (Italy) telescopes were derived from a combination of terrestrial surveying methods to account for gravitational deformations. After applying these models in geodetic VLBI data analysis, estimates of the antenna reference point positions are shifted upward by 8.9 and 6.7 mm, respectively. The impact on other parameters is negligible. To simulate the impact of antenna gravitational deformations on the entire VLBI network, lacking measurements for other telescopes, we rescaled the SPV models of Medicina and Noto for other antennas according to their size. The effects of the simulations are changes in VLBI heights in the range [−3, 73] mm and a net scale increase of 0.3–0.8 ppb. The height bias is larger than random errors of VLBI position estimates, implying the possibility of significant scale distortions related to antenna gravitational deformations. This demonstrates the need to precisely measure gravitational deformations of other VLBI telescopes, to derive their precise SPV models and to apply them in routine geodetic data analysis.     

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.     

The International VLBI Service for Geodesy and Astrometry (IVS): current capabilities and future prospects

Very Long Baseline Interferometry (VLBI) plays a unique and fundamental role in the maintenance of the global (terrestrial and celestial) reference frames, which are required for precise positioning in many research areas such as the understanding and monitoring of global changes, and for space missions. The International VLBI Service for Geodesy and Astrometry (IVS) coordinates the global VLBI components and resources on an international basis. The service is tasked by the International Association of Geodesy (IAG) and International Astronomical Union (IAU) to provide products for the realization of the Celestial Reference Frame (CRF) through the positions of quasars, to deliver products for the maintenance of the terrestrial reference frame (TRF), such as station positions and their changes with time, and to generate products describing the rotation and orientation of the Earth. In particular, VLBI uniquely provides direct observations of nutation parameters and of the time difference UT1-UTC. This paper summarizes the evolution and current status of the IVS. It points out the activities to improve further on the product quality to meet future service requirements.     

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

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

On the consistency of the current conventional EOP series and the celestial and terrestrial reference frames

Precise transformation between the celestial reference frames (CRF) and terrestrial reference frames (TRF) is needed for many purposes in Earth and space sciences. According to the Global Geodetic Observing System (GGOS) recommendations, the accuracy of positions and stability of reference frames should reach 1 mm and 0.1 mm year\(^{-1}\), and thus, the Earth Orientation Parameters (EOP) should be estimated with similar accuracy. Different realizations of TRFs, based on the combination of solutions from four different space geodetic techniques, and CRFs, based on a single technique only (VLBI, Very Long Baseline Interferometry), might cause a slow degradation of the consistency among EOP, CRFs, and TRFs (e.g., because of differences in geometry, orientation and scale) and a misalignment of the current conventional EOP series, IERS 08 C04. We empirically assess the consistency among the conventional reference frames and EOP by analyzing the record of VLBI sessions since 1990 with varied settings to reflect the impact of changing frames or other processing strategies on the EOP estimates. Our tests show that the EOP estimates are insensitive to CRF changes, but sensitive to TRF variations and unmodeled geophysical signals at the GGOS level. The differences between the conventional IERS 08 C04 and other EOP series computed with distinct TRF settings exhibit biases and even non-negligible trends in the cases where no differential rotations should appear, e.g., a drift of about 20 \(\upmu \)as year\(^{-1 }\)in \(y_{\mathrm{pol }}\) when the VLBI-only frame VTRF2008 is used. Likewise, different strategies on station position modeling originate scatters larger than 150 \(\upmu \)as in the terrestrial pole coordinates.     

Radio source instability in VLBI analysis

The source position time-series for many of the frequently observed radio sources in the NASA geodetic very long baseline interferometry (VLBI) program show systematic linear and non-linear variations of as much as 0.5 mas (milli-arc-seconds) to 1.0 mas, due mainly to source structure changes. In standard terrestrial reference frame (TRF) geodetic solutions, it is a common practice to only estimate a global source position for each source over the entire history of VLBI observing sessions. If apparent source position variations are not modeled, they produce corresponding systematic variations in estimated Earth orientation parameters (EOPs) at the level of 0.02–0.04 mas in nutation and 0.01–0.02 mas in polar motion. We examine the stability of position time-series of the 107 radio sources in the current NASA geodetic source catalog since these sources have relatively dense observing histories from which it is possible to detect systematic variations. We consider different strategies for handling source instabilities where we (1) estimate the positions of unstable sources for each session they are observed, or (2) estimate spline parameters or rate parameters for sources chosen to fit the specific variation seen in the position-time series. We found that some strategies improve VLBI EOP accuracy by reducing the biases and weighted root mean square differences between measurements from independent VLBI networks operating simultaneously. We discuss the problem of identifying frequently observed unstable sources and how to identify new sources to replace these unstable sources in the NASA VLBI geodetic source catalog.     

GPS非差相位精密单点定位方法与实现

探讨了精密单点定位的基本原理、处理方法、所涉及的误差改正及数据处理中的一些关键技术;采用直接内插IGS卫星精密星历的方法代替利用IGS跟踪站进行轨道精化方法计算卫星轨道参数,对现有精密单点定位计算方法进行了简化,使之更具有实用性。最后利用自主研发的精密非差单点定位软件计算和分析了实测数据。计算结果表明,经过大约15 min的初始化后,非差相位单历元的定位结果精确度在X,Y,Z方向上均优于20 cm。     

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.     

Present-day crustal deformation in China relative to ITRF97 kinematic plate model

 A global plate motion model is established based on the ITRF97 velocity fields and geological model NUVEL1. Sub-plate models are estimated by using the velocity fields derived from 45 global positioning system (GPS) sites under the ITRF97 reference frame in China. Comparisons between space geodesy and geological models are given. It is found that the Euler vector of the AFRC–EURA pair has an obvious discrepancy between space geodetic and geological models. The motion patterns of tectonic blocks in China predicted by GPS are consistent with those of geological data on the whole. Received: 9 November 2000 / Accepted: 17 September 2001     

Implementation and testing of the gridded Vienna Mapping Function 1 (VMF1)

The new gridded Vienna Mapping Function (VMF1) was implemented and compared to the well-established site-dependent VMF1, directly and by using precise point positioning (PPP) with International GNSS Service (IGS) Final orbits/clocks for a 1.5-year GPS data set of 11 globally distributed IGS stations. The gridded VMF1 data can be interpolated for any location and for any time after 1994, whereas the site-dependent VMF1 data are only available at selected IGS stations and only after 2004. Both gridded and site-dependent VMF1 PPP solutions agree within 1 and 2 mm for the horizontal and vertical position components, respectively, provided that respective VMF1 hydrostatic zenith path delays (ZPD) are used for hydrostatic ZPD mapping to slant delays. The total ZPD of the gridded and site-dependent VMF1 data agree with PPP ZPD solutions with RMS of 1.5 and 1.8 cm, respectively. Such precise total ZPDs could provide useful initial a priori ZPD estimates for kinematic PPP and regional static GPS solutions. The hydrostatic ZPDs of the gridded VMF1 compare with the site-dependent VMF1 ZPDs with RMS of 0.3 cm, subject to some biases and discontinuities of up to 4 cm, which are likely due to different strategies used in the generation of the site-dependent VMF1 data. The precision of gridded hydrostatic ZPD should be sufficient for accurate a priori hydrostatic ZPD mapping in all precise GPS and very long baseline interferometry (VLBI) solutions. Conversely, precise and globally distributed geodetic solutions of total ZPDs, which need to be linked to VLBI to control biases and stability, should also provide a consistent and stable reference frame for long-term and state-of-the-art numerical weather modeling.     

利用北斗系统建立和维持国家大地坐标参考框架的方法研究

2000中国大地坐标系统（China Geodetic Coordinate System 2000,CGCS2000）的建立和维持主要依赖于GPS技术,不利于保障国家时空信息安全。中国北斗卫星导航系统（BeiDou navigation satellite system,BDS）提供亚太区域服务,可满足中国及周边地区高精度定位导航应用需求,对建立和维持国家大地坐标参考框架具有重要意义。研究利用已建成的北斗基准站网观测数据,实现基于BDS技术、并与国际地球参考框架（International Terrestrial Reference Frame,ITRF）一致的国家大地坐标参考框架,为今后国家级和全球性北斗坐标参考框架（BeiDou Terrestrial Reference Frame,BTRF）的建立和维持提供理论基础和方法支撑。初步计算结果表明,积累2 a以上的观测数据,利用单独BDS数据可以获得与GPS精度相当的水平速度场,精度约为2~3 mm/a。基于单独BDS数据,测站残差平面和高程的重复性分别可优于0.8 cm和1.7 cm。利用BDS数据已可监测到测站高程方向的季节性变化。此外,还对单独BDS与GPS数据计算的坐标可能存在的与经纬度相关的系统误差进行了分析。总体来说,目前的北斗系统可满足建立和维持中国cm级大地坐标框架的需求。     

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.     

Co-location of space geodetic techniques carried out at the Geodetic Observatory Wettzell using a closure in time and a multi-technique reference target

The quality of the links between the different space geodetic techniques (VLBI, SLR, GNSS and DORIS) is still one of the major limiting factors for the realization of a unique global terrestrial reference frame that is accurate enough to allow the monitoring of the Earth system, i.e., of processes like sea level change, postglacial rebound and silent earthquakes. According to the specifications of the global geodetic observing system of the International Association of Geodesy, such a reference frame should be accurate to 1 mm over decades, with rates of change stable at the level of 0.1 mm/year. The deficiencies arise from inaccurate or incomplete local ties at many fundamental sites as well as from systematic instrumental biases in the individual space geodetic techniques. Frequently repeated surveys, the continuous monitoring of antenna heights and the geometrical mount stability (Lösler et al. in J Geod 90:467–486, 2016.  https://doi.org/10.1007/s00190-016-0887-8) have not provided evidence for insufficient antenna stability. Therefore, we have investigated variations in the respective system delays caused by electronic circuits, which is not adequately captured by the calibration process, either because of subtle differences in the circuitry between geodetic measurement and calibration, high temporal variability or because of lacking resolving bandwidth. The measured system delay variations in the electric chain of both VLBI- and SLR systems reach the order of 100 ps, which is equivalent to 3 cm of path length. Most of this variability is usually removed by the calibrations but by far not all. This paper focuses on the development of new technologies and procedures for co-located geodetic instrumentation in order to identify and remove systematic measurement biases within and between the individual measurement techniques. A closed-loop optical time and frequency distribution system and a common inter-technique reference target provide the possibility to remove variable system delays. The main motivation for the newly established central reference target, locked to the station clock, is the combination of all space geodetic instruments at a single reference point at the observatory. On top of that it provides the unique capability to perform a closure measurement based on the observation of time.