Impact of the VLBA on reference frames and earth orientation studies

The geodetic VLBI community began using VLBA antennas in 1989 for geodesy and astrometry. We examine how usage of the VLBA has improved the celestial reference frame, the terrestrial reference frame, and Earth orientation parameters. Without the VLBA, ICRF2 would have had only 1011 sources instead of 3414. ICRF3 will contain at least 4121 sources, with approximately 70 % or more coming exclusively from VLBA astrometry and geodesy sessions. The terrestrial reference frame is also more stable and precise due to VLBA geodesy sessions. Approximately two dozen geodesy stations that have participated in VLBA sessions show average position formal errors that are \(\sim \)13–14 % better in the horizontal components and \(\sim \)5 % better in the vertical component than would be expected solely from the increased number of observations. Also the Earth orientation parameters obtained from the RDV sessions represent the most accurate EOP series of any of the long-term VLBI session types.     

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.     

The effects of frequency-dependent quasar variability on the celestial reference frame

We examine the relationship between source position stability and astrophysical properties of radio-loud quasars making up the International Celestial Reference Frame (ICRF2). Understanding this relationship is important for improving quasar selection and analysis strategies, and therefore reference frame stability. We construct flux density time series, known as light curves, for 95 of the most frequently observed ICRF2 quasars at both the 2.3 and 8.4 GHz geodetic very long baseline interferometry (VLBI) observing bands. Because the appearance of new quasar components corresponds to an increase in quasar flux density, these light curves alert us about potential changes in source structure before they appear in VLBI images. We test how source position stability depends on three astrophysical parameters: (1) flux density variability at X band; (2) time lag between flares in S and X bands; (3) spectral index root-mean-square (rms), defined as the variability in the ratio between S and X band flux densities. We find that the time lag between S and X band light curves provides a good indicator of position stability: sources with time lags $<$ 0.06 years are significantly more stable ( $>$ 20 % improvement in weighted rms) than sources with larger time lags. A similar improvement is obtained by observing sources with low $(<$ 0.12) spectral index variability. On the other hand, there is no strong dependence of source position stability on flux density variability in a single frequency band. These findings can be understood by interpreting the time lag between S and X band light curves as a measure of the size of the source structure. Monitoring of source flux density at multiple frequencies therefore appears to provide a useful probe of quasar structure on scales important to geodesy. The observed astrometric position of the brightest quasar component (the core) is known to depend on observing frequency. We show how multi-frequency flux density monitoring may allow the dependence on frequency of the relative core positions along the jet to be elucidated. Knowledge of the position–frequency relation has important implications for current and future geodetic VLBI programs, as well as the alignment between the radio and optical celestial reference frames.     

Generalized inner constraints for geodetic network densification problems

The estimated coordinates from a minimum-constrained (MC) network adjustment are generally influenced by two different error sources, that is the data noise from the available measurements and the so-called datum noise due to random errors in the fiducial coordinates that are used for the datum definition with regard to an external reference frame. Although the latter source does not affect the estimable characteristics of a MC solution, it still contributes a datum-related noise to the estimated positions which reflects the uncertainty of the coordinate system itself for the adjusted network. The aim of this paper is to develop a new type of MCs which minimizes both of the aforementioned effects in the final coordinates of an adjusted network. This particular problem has not been treated in the geodetic literature and the solution which is presented herein offers an elegant unification of the classic inner constraints into a more general framework for the datum choice problem of network optimization theory. Furthermore, the findings of our study provide a useful and rigorous tool for frame densification problems by means of an optimal MC adjustment in geodetic networks.     

A Global Terrestrial Reference Frame from simulated VLBI and SLR data in view of GGOS

In this study, we assess the impact of two combination strategies, namely local ties (LT) and global ties (GT), on the datum realization of Global Terrestrial Reference Frames in view of the Global Geodetic Observing System requiring 1 mm-accuracy. Simulated Very Long Baseline Interferometry (VLBI) and Satellite Laser Ranging (SLR) data over a 7 year time span was used. The LT results show that the geodetic datum can be best transferred if the precision of the LT is at least 1 mm. Investigating different numbers of LT, the lack of co-located sites on the southern hemisphere is evidenced by differences of 9 mm in translation and rotation compared to the solution using all available LT. For the GT, the combination applying all Earth rotation parameters (ERP), such as pole coordinates and UT1-UTC, indicates that the rotation around the Z axis cannot be adequately transferred from VLBI to SLR within the combination. Applying exclusively the pole coordinates as GT, we show that the datum can be transferred with mm-accuracy within the combination. Furthermore, adding artificial stations in Tahiti and Nigeria to the current VLBI network results in an improvement in station positions by 13 and 12%, respectively, and in ERP by 17 and 11%, respectively. Extending to every day VLBI observations leads to 65% better ERP estimates compared to usual twice-weekly VLBI observations.     

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.     

使用参照物的射电望远镜中心坐标测量方法

为满足新建射电望远镜在单站或多站联合进行的深空探测或射电天文观测对中心坐标精确测定的要求,提出一种利用已知精确中心坐标的望远镜作为参照物,测量地平式射电望远镜中心点坐标的方法和测量数据处理方法。这一方法对场地和设备的要求较低,能够得到毫米级或亚毫米级的位置精度。尤其适合对天线阵列的中心位置进行测量。对国家天文台密云观测站的40 m射电望远镜进行了中心坐标测量,位置均方根误差为2.312 mm,满足了后续的观测工作对其位置的需求。     

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.     

坐标基准维持与动态监测网数据处理

坐标基准的维持不仅要求有高精度的连续大地测量观测或定期复测,还要求采用合理的数据处理原则与方法。分析了不同平差方法对基准的影响,其中包括经典整体平差、序贯平差和自适应序贯平差。利用GPS监测网数据进行了计算与分析,结果表明,当个别点位发生异常移动或变形时,序贯平差会使坐标基准产生实际的扭曲,而自适应序贯平差则能有效控制异常点位对平差基准的影响,保证了坐标基准的可靠性和现势性。     

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.     

Contributions to reference systems from Lunar Laser Ranging using the IfE analysis model

Lunar Laser Ranging (LLR) provides various quantities related to reference frames like Earth orientation parameters, coordinates and velocities of ground stations in the Earth-fixed frame and selenocentric coordinates of the lunar retro-reflectors. This paper presents the recent results from LLR data analysis at the Institut für Erdmessung, Leibniz Universität Hannover, based on all LLR data up to the end of 2016. The estimates of long-periodic nutation coefficients with periods between 13.6 days and 18.6 years are obtained with an accuracy in the order of 0.05–0.7 milliarcseconds (mas). Estimations of the Earth rotation phase \(\Delta \)UT are accurate at the level of 0.032 ms if more than 14 normal points per night are included. The tie between the dynamical ephemeris frame to the kinematic celestial frame is estimated from pure LLR observations by two angles and their rates with an accuracy of 0.25 and 0.02 mas per year. The estimated station coordinates and velocities are compared to the ITRF2014 solution and the geometry of the retro-reflector network with the DE430 solution. The given accuracies represent 3 times formal errors of the parameter fit. The accuracy for \(\Delta \)UT is based on the standard deviation of the estimates with respect to the reference C04 solution.     

极坐标测量变形监测网稳定基准点组的检验

分析了全站仪极坐标测量变形监测网的特点,提出了利用坐标转换修正法进行稳定基准点组的检验方法,实现了变形监测网稳定基准点组的检验。     

Error study for the determination of the center of mass of the earth from Pageos observations

The least squares adjustment for the station coordinates of the worldwide satellite triangulation net gives, in addition to the station positions, the positions of the satellites observed in a simultaneous event from two or more stations. Through successive events short arcs can be fitted and the coordinates of the center of mass of the earth can be determined. It is shown that the center of mass will be obtained with an accuracy consistent with the accuracy of the station positions, if well distributed Pageos arcs of about one quarter of a revolution are taken and if three events per arc are given.     

利用激光干涉测距三维网的加权秩亏自由网平差

激光跟踪仪在激光干涉测距模式下,其测距精度要远远高于测角精度,利用极坐标测量的空间点坐标精度主要受测角误差的影响,利用激光跟踪仪的高精度激光干涉测距值构成空间三维激光干涉测边网,消除测角误差对空间点位的影响,采用基于重心基准的加权秩亏网平差模型进行整网平差。定向点坐标近似值采用极坐标方法确定,测站中心坐标近似值采用距离后方交会解算,通过附加约束矩阵精密求解测站点和定向点的三维坐标值。实际计算结果表明,该模型在12m测量范围内将激光跟踪仪的点位精度由87μm(标称值)提高到了27μm。     

Nonlinear analysis of the three-dimensional datum transformation [conformal group ℂ7(3)]

 The weighted Procrustes algorithm is presented as a very effective tool for solving the three-dimensional datum transformation problem. In particular, the weighted Procrustes algorithm does not require any initial datum parameters for linearization or any iteration procedure. As a closed-form algorithm it only requires the values of Cartesian coordinates in both systems of reference. Where there is some prior information about the variance–covariance matrix of the two sets of Cartesian coordinates, also called pseudo-observations, the weighted Procrustes algorithm is able to incorporate such a quality property of the input data by means of a proper choice of weight matrix. Such a choice is based on a properly designed criterion matrix which is discussed in detail. Thanks to the weighted Procrustes algorithm, the problem of incorporating the stochasticity measures of both systems of coordinates involved in the seven parameter datum transformation problem [conformal group ℂ₇(3)] which is free of linearization and any iterative procedure can be considered to be solved. Illustrative examples are given. Received: 7 January 2002 / Accepted: 9 September 2002 Correspondence to: E. W. Grafarend     

Estimated SLR station position and network frame sensitivity to time-varying gravity

This paper evaluates the sensitivity of ITRF2008-based satellite laser ranging (SLR) station positions estimated weekly using LAGEOS-1/2 data from 1993 to 2012 to non-tidal time-varying gravity (TVG). Two primary methods for modeling TVG from degree-2 are employed. The operational approach applies an annual GRACE-derived field, and IERS recommended linear rates for five coefficients. The experimental approach uses low-order/degree $4\times 4$ coefficients estimated weekly from SLR and DORIS processing of up to 11 satellites (tvg4x4). This study shows that the LAGEOS-1/2 orbits and the weekly station solutions are sensitive to more detailed modeling of TVG than prescribed in the current IERS standards. Over 1993–2012 tvg4x4 improves SLR residuals by 18 % and shows 10 % RMS improvement in station stability. Tests suggest that the improved stability of the tvg4x4 POD solution frame may help clarify geophysical signals present in the estimated station position time series. The signals include linear and seasonal station motion, and motion of the TRF origin, particularly in Z. The effect on both POD and the station solutions becomes increasingly evident starting in 2006. Over 2008–2012, the tvg4x4 series improves SLR residuals by 29 %. Use of the GRGS RL02 $50\times 50$ series shows similar improvement in POD. Using tvg4x4, secular changes in the TRF origin Z component double over the last decade and although not conclusive, it is consistent with increased geocenter rate expected due to continental ice melt. The test results indicate that accurate modeling of TVG is necessary for improvement of station position estimation using SLR data.     

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.     

GPS形变监测网基线处理中系统误差的分析

分析了高精度GPS形变监测网基线处理中系统误差产生的原因，分类及其对基线处理的影响，并在此基础上提出消除和削弱这些系统误差影响的一些原则和算法。     

Enhancement of the accuracy of single-epoch GPS positioning for long baselines by local ionospheric modelling

Single-epoch relative GPS positioning has many advantages, especially for monitoring dynamic targets. In this technique, errors occurring in previous epochs cannot affect the position accuracy at the current epoch, but careful processing is required, and resolving carrier phase ambiguities is essential. Statistical ambiguity resolution functions have been used to determine the best values of these ambiguities. The function inputs include as a minimum the known base station position, the approximate roving antenna “seed” position, and the dual-frequency carrier phase measurements from both receivers. We investigate different solutions to find the ambiguity function inputs that achieve the highest ambiguity resolution success rate. First, we address the rover seed position. A regionally filtered undifferenced pseudorange coordinate solution proves better than a double-differenced one. Multipath errors approximately repeat themselves every sidereal day in the case of static or quasi-static antennas; applying a sidereal filter to the pseudorange-derived positions mitigates their effects. Second, we address the relative carrier phase measurements, which for medium to long baselines are significantly affected by ionospheric propagation errors imperfectly removed during differencing. In addition to the International GNSS Service ionospheric model, we generate a local pseudorange-based ionospheric correction. Applying this correction improves the quality of the phase measurements, leading to more successful ambiguity resolution. Temporally smoothing the correction by means of a Kalman filter further improves the phase measurements. For baselines in the range 60–120 km, the mean absolute deviation of single-epoch coordinates improves to 10–20 cm, from 30–50 cm in the default case.     

海洋空间信息基准技术进展与发展方向

简要分析了海洋空间信息基准的基本发展历程和当前的技术状态,梳理和评述了近年来关于验潮站垂直基准确定方法的改进与完善、验潮站控制下的垂直基准的实现形式、长期验潮站在维持垂直基准方面的作用、海洋潮汐模型构建、特征潮面模型构建与转换等有关海洋垂直基准及与其相关的海洋潮汐研究的主要进展,预测了海洋空间信息基准在海上高精度位置服务、垂直基准面及转换模型精化、空间基准应用等重点发展方向。