Modeling thermal deformation of VLBI antennas with a new temperature model

Temperature variations at very long baseline interferometry (VLBI) sites cause thermal deformations of the VLBI antennas and corresponding displacements of the VLBI reference points. The thermal deformation effects typically contain seasonal and daily signatures. The amplitudes of the annual vertical motion of the antenna reference point can reach several millimeters, depending on the design of the antenna structure, on the material, and on the environmental effects such as global station position, station height and climatology effects. Simple methods to correct this effect use the difference of the environmental temperature with respect to a defined reference temperature, the antenna dimensions, the elevation of the antenna, the material of antenna structure. Applying these simple models for thermal deformation in the VLBI data analysis improves the baseline length repeatability by 3.5%. A comparison of these simple models with local thermal deformation measurements at the antennas in Onsala and Wettzell show that the local measurements and the modeled corrections agree well when the temperature of the antenna structure is used, but agree less good when the surrounding air temperatures are used. To overcome this problem we present a method to model temperature penetration into the antenna structures, that allows to model thermal deformation effects that agree with the observed vertical deformation of the Onsala and Wettzell radio telescopes with a root mean square deviation of 0.07 and 0.13 mm, respectively. Possible implementations in the VLBI analysis are presented, and the definition of an adequate reference temperature is discussed.     

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.     

Simulations of VLBI observations of a geodetic satellite providing co-location in space

We performed Monte Carlo simulations of very-long-baseline interferometry (VLBI) observations of Earth-orbiting satellites incorporating co-located space-geodetic instruments in order to study how well the VLBI frame and the spacecraft frame can be tied using such measurements. We simulated observations of spacecraft by VLBI observations, time-of-flight (TOF) measurements using a time-encoded signal in the spacecraft transmission, similar in concept to precise point positioning, and differential VLBI (D-VLBI) observations using angularly nearby quasar calibrators to compare their relative performance. We used the proposed European Geodetic Reference Antenna in Space (E-GRASP) mission as an initial test case for our software. We found that the standard VLBI technique is limited, in part, by the present lack of knowledge of the absolute offset of VLBI time to Coordinated Universal Time at the level of microseconds. TOF measurements are better able to overcome this problem and provide frame ties with uncertainties in translation and scale nearly a factor of three smaller than those yielded from VLBI measurements. If the absolute time offset issue can be resolved by external means, the VLBI results can be significantly improved and can come close to providing 1 mm accuracy in the frame tie parameters. D-VLBI observations with optimum performance assumptions provide roughly a factor of two higher uncertainties for the E-GRASP orbit. We additionally simulated how station and spacecraft position offsets affect the frame tie performance.     

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.     

Diurnal atmospheric forcing and temporal variations of the nutation amplitudes

A 29-year time-series of four-times-daily atmospheric effective angular momentum (EAM) estimates is used to study the atmospheric influence on nutation. The most important atmospheric contributions are found for the prograde annual (77 μas), retrograde annual (53 as), prograde semiannual (45 as), and for the constant offset of the pole (δψsinɛ₀=−86 as, δɛ=77 as). Among them only the prograde semiannual component is driven mostly by the wind term of the EAM function, while in all other cases the pressure term is dominant. These are nonnegligible quantities which should be taken into account in the new theory of nutation. Comparison with the VLBI corrections to the IAU 1980 nutation model taking into account the ocean tide contribution yields good agreement for the prograde annual and semiannual nutations. We also investigated time variability of the atmospheric contribution to the nutation amplitudes by performing the sliding-window least-squares analysis of both the atmospheric excitation and VLBI nutation data. Almost all detected variations of atmospheric origin can be attributed to the pressure term, the biggest being the in-phase annual prograde component (about 30 as) and the retrograde one (as much as 100200 as). These variations, if physical, limit the precision of classical modeling of nutation to the level of 0.1 mas. Comparison with the VLBI data shows significant correlation for the retrograde annual nutation after 1989, while for the prograde annual term there is a high correlation in shape but the size of the atmospherically driven variations is about three times less than deduced from the VLBI data. This discrepancy in size can be attributed either to inaccuracy of the theoretical transfer function or the frequency-dependent ocean response to the pressure variations. Our comparison also yields a considerably better agreement with the VLBI nutation data when using the EAM function without the IB correction for ocean response, which indicates that this correction is not adequate for nearly diurnal variations. Received: 10 September 1997 / Accepted: 5 March 1998     

A comparison of GPS, VLBI and model estimates of ocean tide loading displacements

In recent years, ocean tide loading displacements (OTLD) have been measured using the Global Positioning System (GPS) and Very Long Baseline Interferometry (VLBI). This study assesses the accuracy of GPS measurements of OTLD by comparison with VLBI measurements and estimates derived from numerical ocean tide models. A daily precise point positioning (PPP) analysis was carried out on ∼11 years of GPS data for each of 25 sites that have previous OTLD estimates based on data from co-located VLBI sites. Ambiguities were fixed to integer values where possible. The resulting daily estimates of OTLD, at eight principal diurnal and semi-diurnal tidal frequencies, were combined to give GPS measurements of OTLD at each site. The 3D GPS and VLBI measurements of OTLD were compared with estimates computed (by convolution with Green’s functions) from five modern ocean tide models (CSR4.0, FES2004, GOT00.2, NAO99b and TPXO6.2). The GPS/model agreement is shown to be similar to the VLBI/model agreement. In the important radial direction, the GPS/model misfit is shown to be smaller than the VLBI/model misfit for seven of the eight tidal constituents; the exception being the K2 constituent. Fixing of GPS carrier-phase ambiguities to integer values resulted in a marginal improvement to the GPS/model agreement. Statistically, it is shown there is no significance to the difference between the fit of the GPS and VLBI measurements of OTLD to modelled values. Equally, differences in fit of either the complete set of GPS or VLBI estimates to the five sets of model-derived values cannot be identified with statistical significance. It is thus concluded that, overall, we cannot distinguish between GPS and VLBI measurements of OTLD, and that at the global scale, present ocean tide models are accurate to within the current measurement noise of these techniques.     

Sensitivity of UT1 determined by single-baseline VLBI to atmospheric delay model, terrestrial and celestial reference frames

 At the present time, the daily VLBI observations on the Westford-Wettzell baseline is the only continually running VLBI project for studies of high-frequency Earth rotation variations. An analysis of this experiment with regard to the potential errors in the atmospheric delay model and in adopted celestial and terrestrial reference frames is presented in the paper. A new VLBI geometric delay model is applied and an algorithm for global adjustment for this specific single-baseline VLBI developed. The results over three years show discrepancies at the milliarcsecond level between the daily observations and the adopted atmospheric model as well as the combined celestial reference frame. A significant number of these discrepancies are removed by the global adjustment. Received: 19 August 1996; Accepted: 13 September 1996     

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.     

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.     

VLBI在探月卫星定位中的应用分析

中国实施的"嫦娥"探月工程中,探月卫星的定轨测控系统由我国现有的S频段航天测控网(USB)和甚长基线干涉测量(VLBI)系统组成。系统中,VLBI技术主要为绕月卫星定轨提供卫星的角位置。本文分析了在探月项目中,VLBI单点定位的必要性。探讨了VLBI技术用于探月卫星单点定位的基本原理及其实现方法。通过算例对模拟数据进行处理,检验了方法的正确性。对结果进行分析,得出一些结论。     

Multi-reference evaluation of uncertainty in earth orientation parameter measurements

Uncertainties in polar motion and length-of-day measurements are evaluated empirically using several data series from the space-geodetic techniques of the global positioning system (GPS), satellite laser ranging (SLR), and very long baseline interferometry (VLBI) during 1997–2002. In the evaluation procedure employed here, known as the three-corner hat (TCH) technique, the signal common to each series is eliminated by forming pair-wise differences between the series, thus requiring no assumed values for the “truth” signal. From the variances of the differenced series, the uncertainty of each series can be recovered when reasonable assumptions are made about the correlations between the series. In order to form the pair-wise differences, the series data must be given at the same epoch. All measurement data sets studied here were sampled at noon (UTC); except for the VLBI series, whose data are interpolated to noon and whose UT1 values are also numerically differentiated to obtain LOD. The numerical error introduced to the VLBI values by the interpolation and differentiation is shown to be comparable in magnitude to the values determined by the TCH method for the uncertainties of the VLBI series. The TCH estimates for the VLBI series are corrupted by such numerical errors mostly as a result of the relatively large data intervals. Of the remaining data sets studied here, it is found that the IGS Final combined series has the smallest polar motion and length-of-day uncertainties.     

The impact of radio source structure on European geodetic VLBI measurements

European geodetic very long baseline interferometry (VLBI) sessions (also known as EUROPE sessions) have been carried out on a regular basis for the past 15 years to study relative crustal motions within Europe. These sessions are based on observations of extragalactic radio sources, which serve as distant fiducial marks to establish an accurate and stable celestial reference frame for long-term geodetic measurements. The radio sources, however, are not always point-like on milliarcsecond scales, as VLBI imaging has revealed. In this work, we quantify the magnitude of the expected effect of intrinsic source structure on geodetic bandwidth synthesis delay VLBI measurements for a subset of 14 sources regularly observed during the EUROPE sessions. These sources have been imaged at both X-band (8.4 GHz) and S-band (2.3 GHz) based on dedicated observations acquired with the European VLBI Network (EVN) in November 1996. The results of this calculation indicate that the reference source 0457+024 causes significant structural effects in measurements obtained on European VLBI baselines (about 10 picoseconds on average), whereas most of the other sources produce effects that are only occasionally larger than a few picoseconds. Applying the derived source structure models to the data of the EUROPE5-96 session carried out at the same epoch as the EVN experiment shows no noticeable changes in the estimated VLBI station locations.     

基于联邦卡尔曼滤波的月球车联合定位方法

针对月球车定位,设计和推导了一种基于甚长基线干涉测量（very long baseline interferometry,VLBI）和天文导航系统（celestial navigation system,CNS）相结合的月球车联合定位方法,并采用联邦卡尔曼滤波来实现位置信息的最优估计以增强系统的可靠性和容错能力。最后,通过嫦娥-3号（CE-3）实测数据的解算证明了该方法较使用最小二乘法进行联合解算和单独采用VLBI方法定位,可以获得更高的月球车定位精度。同时,也有效地保障了月球车定位的可靠性和稳定性。     

Conventions on thermal expansion modelling of radio telescopes for geodetic and astrometric VLBI

Thermal expansion of radio telescopes has long been recognized as an effect which cannot be neglected in geodetic and astrometric VLBI data analysis if millimeter accuracy is desired. In this article, the author documents the conventions which are being set by the International VLBI Service for Geodesy and Astrometry (IVS) for a consistent modelling of this effect in its routine product generation. For the largest telescopes, the annual cycle of thermal expansion may change the height of the VLBI reference point by as much as 20 mm. However, for telescopes which are used in present-day IVS operations, the variations rather range from 4 to 6 mm.     

Multi-technique comparison of tropospheric zenith delays derived during the CONT02 campaign

In October 2002, 15 continuous days of Very Long Baseline Interferometry (VLBI) data were observed in the Continuous VLBI 2002 (CONT02) campaign. All eight radio telescopes involved in CONT02 were co-located with at least one other space-geodetic technique, and three of them also with a Water Vapor Radiometer (WVR). The goal of this paper is to compare the tropospheric zenith delays observed during CONT02 by VLBI, Global Positioning System (GPS), Doppler Orbitography Radiopositioning Integrated by Satellite (DORIS) and WVR and to compare them also with operational pressure level data from the European Centre for Medium-Range Weather Forecasts (ECMWF). We show that the tropospheric zenith delays from VLBI and GPS are in good agreement at the 3–7 mm level. However, while only small biases can be found for most of the stations, at Kokee Park (Hawaii, USA) and Westford (Massachusetts, USA) the zenith delays derived by GPS are larger by more than 5 mm than those from VLBI. At three of the four DORIS stations, there is also a fairly good agreement with GPS and VLBI (about 10 mm), but at Kokee Park the agreement is only at about 30 mm standard deviation, probably due to the much older installation and type of DORIS equipment. This comparison also allows testing of different DORIS analysis strategies with respect to their real impact on the precision of the derived tropospheric parameters. Ground truth information about the zenith delays can also be obtained from the ECMWF numerical weather model and at three sites using WVR measurements, allowing for comparisons with results from the space-geodetic techniques. While there is a good agreement (with some problems mentioned above about DORIS) among the space-geodetic techniques, the comparison with WVR and ECMWF is at a lower accuracy level. The complete CONT02 data set is sufficient to derive a good estimate of the actual precision and accuracy of each geodetic technique for applications in meteorology.     

Absolute phase center corrections of satellite and receiver antennas

Results of the estimation of azimuth-dependent phase center variations (PCVs) of GPS satellite antennas using global GPS data are presented. Significant variations of up to ±3–4 mm that are demonstrated show excellent repeatability over eight years. The application of the azimuthal PCVs besides the nadir-dependent ones will lead to a further reduction in systematic antenna effects. In addition, the paper focuses on the benefit of a possible transition from relative to absolute PCVs. Apart from systematic changes in the global station coordinates, one can expect the GPS results to be less dependent on the elevation cut-off angle. This, together with the significant reduction of tropospheric zenith delay biases between GPS and VLBI, stands for an important step toward more consistency between different space geodetic techniques.     

河外射电源地面VLBI观测的相对论时延模型

重新研究了河外射电源地面VLBI观测的相时论时间延迟模型，给出了一个严格解析的表达式。根据公式可以得到数据处理中常用的Zhu-Groten模型、Shapiro模型和IERS推荐模型，详细地给出了时间延迟理论模型的计算步骤，介绍了各种时间系统的转换和使用方法。     

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.     

Methodology for the combination of sub-daily Earth rotation from GPS and VLBI observations

A combination procedure of Earth orientation parameters from Global Positioning System (GPS) and Very Long Baseline Interferometry (VLBI) observations was developed on the basis of homogeneous normal equation systems. The emphasis and purpose of the combination was the determination of sub-daily polar motion (PM) and universal time (UT1) for a long time-span of 13 years. Time series with an hourly resolution and a model for tidal variations of PM and UT1-TAI (dUT1) were estimated. In both cases, 14-day nutation corrections were estimated simultaneously with the ERPs. Due to the combination procedure, it was warranted that the strengths of both techniques were preserved. At the same time, only a minimum of de-correlating or stabilizing constraints were necessary. Hereby, a PM time series was determined, whose precision is mainly dominated by GPS observations. However, this setup benefits from the fact that VLBI delivered nutation and dUT1 estimates at the same time. An even bigger enhancement can be seen for the dUT1 estimation, where the high-frequency variations are provided by GPS, while the long term trend is defined by VLBI. The estimated combined tidal PM and dUT1 model was predominantly determined from the GPS observations. Overall, the combined tidal model for the first time completely comprises the geometrical benefits of VLBI and GPS observations. In terms of root mean squared (RMS) differences, the tidal amplitudes agree with other empirical single-technique tidal models below 4 μas in PM and 0.25 μs in dUT1. The noise floor of the tidal ERP model was investigated in three ways resulting in about 1 μas for diurnal PM and 0.07 μs for diurnal dUT1 while the semi-diurnal components have a slightly better accuracy.     

VLBI观测双线极化条纹拟合方法

提出了一种可用于下一代VLBI观测系统(VGOS)的双线极化条纹拟合方法。现有的VLBI观测模式采用的是右圆极化(RCP),而VGOS系统采用的是双线极化。本文方法包括校正和组合条纹拟合两部分。校正部分选择一颗强源作为参考源,分别得到不同极化方式下的通道时延及相位校正数据,用于目标源的校正。组合条纹拟合部分将4种极化分量的可见度数据组合成伪Stokes分量,通过搜索差分星位角使伪Stokes分量的幅值达到最大,从而获得最终的时延观测量。与单极化条纹拟合相比,组合极化获得的条纹具有更高的信噪比(SNR)及更小的条纹相位弥散度。