Impact of seasonal station motions on VLBI UT1 intensives results

UT1 estimates obtained from the very long baseline interferometry (VLBI) Intensives data depend on the station displacement model used during processing. In particular, because of seasonal variations, the instantaneous station position during the specific intensive session differs from the position predicted by the linear model generally used. This can cause systematic errors in UT1 Intensives results. In this paper, we first investigated the seasonal signal in the station displacements for the 5 VLBI antennas participating in UT1 Intensives observing programs, along with the 8 collocated GPS stations. It was found that a significant annual term is present in the time series for most stations, and its amplitude can reach 8 mm in the height component, and 2 mm in horizontal components. However, the annual signals found in the displacements of the collocated VLBI and GPS stations at some sites differ substantially in amplitude and phase. The semiannual harmonics are relatively small and unstable, and for most stations no prevailing signal was found in the corresponding frequency band. Then two UT1 Intensives series were computed with and without including the seasonal term found in the previous step in the station movement model. Comparison of these series has shown that neglecting the seasonal station position variations can cause a systematic error in UT1 estimates, which can exceed 1  $\upmu $ s, depending on the observing program.     

Annual deformation signals from homogeneously reprocessed VLBI and GPS height time series

The first part of this paper compares homogeneously reprocessed Very Long Baseline Interferometry (VLBI) and Global Positioning System (GPS) long-term height series from 1994 to 2007. The data analysis used fully adapted state-of-the-art models (like VMF1 and a priori zenith delays from ECMWF) for the GPS and VLBI processing. The series are compared in terms of long-term non-linear behaviour, harmonic and mean annual signals (not necessarily of harmonic nature). The similarity between both techniques is very good (especially the mean annual signals), which is assumed to be due to the adapted models and consistent reprocessing of both series. As two almost independent observing techniques see the same annually recurring signals at almost all co-located sites, we expect a good geophysical interpretability as integral vertical deformation. For the second part of this paper, the height time series of 161 suitable GPS sites (of the same solution as before) are used to determine a harmonic and a mean annual signal for each of them. Comparing the annual signals for this big dataset visually to GRACE-determined load deformations described in other publications, we find good agreement. This puts emphasis to the assumption that our height data have a lot of potential to be interpreted as geophysical signals. Out of these 161, 131 are grouped to 55 clusters, if at least two nearby (some thousand kilometres) sites show similar mean annual signals, which are thus confirmed to be real regional deformation, not local or technical artefacts. These 55 signals are presented on a “world map” of regional average mean annual height signals, as easy-to-handle tool to validate geophysical models. The data of these measured regional mean annual signals can be downloaded from a web-page for numerical analysis.     

Multi-technique comparison of troposphere zenith delays and gradients during CONT08

CONT08 was a 15 days campaign of continuous Very Long Baseline Interferometry (VLBI) sessions during the second half of August 2008 carried out by the International VLBI Service for Geodesy and Astrometry (IVS). In this study, VLBI estimates of troposphere zenith total delays (ZTD) and gradients during CONT08 were compared with those derived from observations with the Global Positioning System (GPS), Doppler Orbitography and Radiopositioning Integrated by Satellite (DORIS), and water vapor radiometers (WVR) co-located with the VLBI radio telescopes. Similar geophysical models were used for the analysis of the space geodetic data, whereas the parameterization for the least-squares adjustment of the space geodetic techniques was optimized for each technique. In addition to space geodetic techniques and WVR, ZTD and gradients from numerical weather models (NWM) were used from the European Centre for Medium-Range Weather Forecasts (ECMWF) (all sites), the Japan Meteorological Agency (JMA) and Cloud Resolving Storm Simulator (CReSS) (Tsukuba), and the High Resolution Limited Area Model (HIRLAM) (European sites). Biases, standard deviations, and correlation coefficients were computed between the troposphere estimates of the various techniques for all eleven CONT08 co-located sites. ZTD from space geodetic techniques generally agree at the sub-centimetre level during CONT08, and??as expected??the best agreement is found for intra-technique comparisons: between the Vienna VLBI Software and the combined IVS solutions as well as between the Center for Orbit Determination (CODE) solution and an IGS PPP time series; both intra-technique comparisons are with standard deviations of about 3?C6?mm. The best inter space geodetic technique agreement of ZTD during CONT08 is found between the combined IVS and the IGS solutions with a mean standard deviation of about 6?mm over all sites, whereas the agreement with numerical weather models is between 6 and 20?mm. The standard deviations are generally larger at low latitude sites because of higher humidity, and the latter is also the reason why the standard deviations are larger at northern hemisphere stations during CONT08 in comparison to CONT02 which was observed in October 2002. The assessment of the troposphere gradients from the different techniques is not as clear because of different time intervals, different estimation properties, or different observables. However, the best inter-technique agreement is found between the IVS combined gradients and the GPS solutions with standard deviations between 0.2 and 0.7?mm.     

利用GPS和GRACE研究澳大利亚地壳垂向季节性变化

为探究重力场恢复与气候实验(gravity recovery and climate experiment,GRACE)卫星与全球定位系统(global positioning system,GPS)两种独立技术获取的因陆地水储量变化引起的地壳垂向季节性位移的一致性,选取澳大利亚27个GPS站点5～10 a的高程时间序...     

Atmospheric loading corrections at the observation level in VLBI analysis

This paper investigates whether in very long baseline interferometry (VLBI) analysis atmospheric loading corrections should be applied a priori at the observation level or whether it is sufficient to correct for atmospheric loading effects a posteriori by adding constant values per session to the estimated station coordinates. Simulated observations at single stations corresponding to the precise point positioning approach of global navigation satellite systems show that the atmospheric loading effect can be fully recovered by a posteriori corrections, i.e., the height differences between both approaches stay well below 1 mm. However, real global VLBI network solutions with sessions from 1984 to 2008 reveal that the effect of neglected atmospheric loading corrections at the stations is distributed to the other stations in the network, thus resulting in station height differences between solutions with observation level and with a posteriori corrections which can be as large as 10 mm and a ‘damping’ effect of the corrections. As soon as the terrestrial reference frame and the corresponding coordinate time series are determined, it would be conceptually wrong to apply atmospheric loading corrections at the VLBI stations. We recommend the rigorous application of atmospheric loading corrections at the observation level to all stations of a VLBI network because the seven parameters for translation, rotation, and in particular the network-scale of VLBI networks are significantly affected.     

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.     

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

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级大地坐标框架的需求。     

空间技术测定地壳运动的相互比较

基线长度变化率反映了测站间的地壳运动情况。采用国际空间数据分析中心提供的GPS和VLBI实测站速度数据，解算了并置站间的基线长度变化率。并置站间不同技术测定的地壳运动的相互比较．可以发现不同技术测定地壳运动的系统差。研究结果表明，VLBI基线长度变化率比GPS的相应值系统性地偏大约2％，GPS基线长度变化率比VLBI的相应值快一个常值(约0．6mm／a)。     

Multi-technique comparisons of 10?years of wet delay estimates on the west coast of Sweden

We present comparisons of 10-year-long time series of the atmospheric zenith wet delay (ZWD), estimated using the global positioning system (GPS), geodetic very long baseline interferometry (VLBI), a water vapour radiometer (WVR), radiosonde (RS) observations, and the reanalysis product of the European Centre for Medium-Range Weather Forecasts (ECMWF). To compare the data sets with each other, a Gaussian filter is applied. The results from 10 GPS–RS comparisons using sites in Sweden and Finland show that the full width at half maximum at which the standard deviation (SD) is a minimum increases with the distance between each pair. Comparisons between three co-located techniques (GPS, VLBI, and WVR) result in mean values of the ZWD differences at a level of a few millimetres and SD of less than 7?mm. The best agreement is seen in the GPS–VLBI comparison with a mean difference of ?3.4?mm and an SD of 5.1?mm over the 10-year period. With respect to the ZWD derived from other techniques, a positive bias of up to ~7?mm is obtained for the ECMWF reanalysis product. Performing the comparisons on a monthly basis, we find that the SD including RS or ECMWF varies with the season, between 3 and 15?mm. The monthly SD between GPS and WVR does not have a seasonal signature and varies from 3 to 7?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.     

地表质量负载对格陵兰岛区域GPS站点坐标时序特征的影响

选取ITRF2008框架下格陵兰岛区域12个GPS站2013年1月-2016年12月期间的日解坐标时间序列作为研究对象,并利用极大似然估计分析地表质量负载改正前后各站点的噪声特性、速度场及周期项振幅。结果表明：站点最优噪声模型主要为白噪声+幂律噪声与白噪声+闪烁噪声,地表质量负载形变修正GPS坐标时序后,明显增加U方向闪烁噪声的成分,平均降低其速度约0.36 mm/a,对水平方向影响较小;同时分别降低高程方向44.1%、14.2%的1 a项、0.5 a项振幅,相反,却增加了水平方向的周期项振幅。     

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.     

Detection of station coordinate discontinuities within the Italian GPS Fiducial Network

The Detection Identification Adaptation (DIA) procedure was applied to the coordinate time-series of some permanent GPS stations belonging to the Italian GPS Fiducial Network (IGFN), of the Italian Space Agency (ASI), to detect discontinuities and to reject outliers. The daily solutions of the stations of Cagliari, Genoa, Medicina, Noto, Turin, Perugia and Venice were computed for the period 1997.0–2003.0 using Bernese GPS software v.4.2. The data were interpolated using a model with a linear term and an annual periodic term. The parameters were estimated by least squares. The DIA procedure was organized to automatically detect discontinuities and outliers. Approximately, 70% of the discontinuities present in the coordinate time-series were identified and their magnitudes were estimated. The identified discontinuities are basically caused by equipment replacement and reference frame changes, but in a few cases the reason is still unknown. With regard to the outliers, roughly 6% of the data were rejected. These data were considered outliers on the base of the level of significance and of the power of the test adopted in this work. Except for the stations of Perugia and Venice, the estimated coordinates agree with ITRF2000 values at the 10 mm level, and the estimated velocities are within a few mm/year of the ITRF2000 values.     

Vertical velocities of European coastal sites derived from continuous GPS observations

Vertical velocities of 30 European permanent Global Positioning System (GPS) stations at or close to tide gauge sites are estimated from more than 3 years of continuous observations. The results of two different solution strategies are presented and compared. The first approach accumulates the daily free network normal equations, the second introduces all sets of daily ellipsoidal height estimates and their covariance matrix into a subsequent common least squares adjustment. In both solutions, mean station heights at a reference epoch, linear vertical velocities, height discontinuities and short period height offsets are estimated. The second approach solves in addition for periodic annual signals and for site-specific pressure loading coefficients. The vertical velocities range from +8 mm/year in the center of the Fennoscandian uplift area to –4 mm/year at a few subsiding locations. Apart from these extrema, most of the sites experience only very small vertical motions. The standard deviations from the second approach providing more realistic error estimates are well below 0.15 mm/year. Some specific data problems are discussed.     

Estimation and evaluation of real-time precipitable water vapor from GLONASS and GPS

The revitalized Russian GLONASS system provides new potential for real-time retrieval of zenith tropospheric delays (ZTD) and precipitable water vapor (PWV) in order to support time-critical meteorological applications such as nowcasting or severe weather event monitoring. In this study, we develop a method of real-time ZTD/PWV retrieval based on GLONASS and/or GPS observations. The performance of ZTD and PWV derived from GLONASS data using real-time precise point positioning (PPP) technique is carefully investigated and evaluated. The potential of combining GLONASS and GPS data for ZTD/PWV retrieving is assessed as well. The GLONASS and GPS observations of about half a year for 80 globally distributed stations from the IGS (International GNSS Service) network are processed. The results show that the real-time GLONASS ZTD series agree quite well with the GPS ZTD series in general: the RMS of ZTD differences is about 8 mm (about 1.2 mm in PWV). Furthermore, for an inter-technique validation, the real-time ZTD estimated from GLONASS-only, GPS-only, and the GPS/GLONASS combined solutions are compared with those derived from very long baseline interferometry (VLBI) at colocated GNSS/VLBI stations. The comparison shows that GLONASS can contribute to real-time meteorological applications, with almost the same accuracy as GPS. More accurate and reliable water vapor values, about 1.5–2.3 mm in PWV, can be achieved when GLONASS observations are combined with the GPS ones in the real-time PPP data processing. The comparison with radiosonde data further confirms the performance of GLONASS-derived real-time PWV and the benefit of adding GLONASS to stand-alone GPS processing.     

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.     

Geodetic accuracy of the Macrometer model V-1000

The Macrometer Model V-1000 is a geodetic positioning instrument that uses the radio signals broadcast by the GPS satellites. The Macrometer is the only GPS user equipment commercially available that does not require any of the GPS codes. The Model V-1000 receives only the 19-cm wavelength signals from GPS although similar instruments, built for the U.S. Air Force Geophysics Laboratory, receive both the 19-and the 24-cm wavelengths. In this paper we summarize the results of two years of field testing of the V-1000. This instrument, observing four or five satellites fof a few hours, yields a point position accurate within several meters in each coordinate: latitude, longitude and ellipsoidal height. All three components of the relative position vector between a pair of points can be determined within 2 parts per million of the distance, given a similar schedule of observations. This accuracy has been obtained for intersite distances from one kilometer to several thousand kilometers. Macrometer is a trademark of Aero Service Division, Western Geophysical Company of American, 8100 Westpark Drive, Houston, Texas, 77063, U.S.A.     

On the long-term stability of terrestrial reference frame solutions based on Kalman filtering

The Global Geodetic Observing System requirement for the long-term stability of the International Terrestrial Reference Frame is 0.1 mm/year, motivated by rigorous sea level studies. Furthermore, high-quality station velocities are of great importance for the prediction of future station coordinates, which are fundamental for several geodetic applications. In this study, we investigate the performance of predictions from very long baseline interferometry (VLBI) terrestrial reference frames (TRFs) based on Kalman filtering. The predictions are computed by extrapolating the deterministic part of the coordinate model. As observational data, we used over 4000 VLBI sessions between 1980 and the middle of 2016. In order to study the predictions, we computed VLBI TRF solutions only from the data until the end of 2013. The period of 2014 until 2016.5 was used to validate the predictions of the TRF solutions against the measured VLBI station coordinates. To assess the quality, we computed average WRMS values from the coordinate differences as well as from estimated Helmert transformation parameters, in particular, the scale. We found that the results significantly depend on the level of process noise used in the filter. While larger values of process noise allow the TRF station coordinates to more closely follow the input data (decrease in WRMS of about 45%), the TRF predictions exhibit larger deviations from the VLBI station coordinates after 2014 (WRMS increase of about 15%). On the other hand, lower levels of process noise improve the predictions, making them more similar to those of solutions without process noise. Furthermore, our investigations show that additionally estimating annual signals in the coordinates does not significantly impact the results. Finally, we computed TRF solutions mimicking a potential real-time TRF and found significant improvements over the other investigated solutions, all of which rely on extrapolating the coordinate model for their predictions, with WRMS reductions of almost 50%.     

Local effects of redundant terrestrial and GPS-based tie vectors in ITRF-like combinations

Tie vectors (TVs) between co-located space geodetic instruments are essential for combining terrestrial reference frames (TRFs) realised using different techniques. They provide relative positioning between instrumental reference points (RPs) which are part of a global geodetic network such as the international terrestrial reference frame (ITRF). This paper gathers the set of very long baseline interferometry (VLBI)–global positioning system (GPS) local ties performed at the observatory of Medicina (Northern Italy) during the years 2001–2006 and discusses some important aspects related to the usage of co-location ties in the combinations of TRFs. Two measurement approaches of local survey are considered here: a GPS-based approach and a classical approach based on terrestrial observations (i.e. angles, distances and height differences). The behaviour of terrestrial local ties, which routinely join combinations of space geodetic solutions, is compared to that of GPS-based local ties. In particular, we have performed and analysed different combinations of satellite laser ranging (SLR), VLBI and GPS long term solutions in order to (i) evaluate the local effects of the insertion of the series of TVs computed at Medicina, (ii) investigate the consistency of GPS-based TVs with respect to space geodetic solutions, (iii) discuss the effects of an imprecise alignment of TVs from a local to a global reference frame. Results of ITRF-like combinations show that terrestrial TVs originate the smallest residuals in all the three components. In most cases, GPS-based TVs fit space geodetic solutions very well, especially in the horizontal components (N, E). On the contrary, the estimation of the VLBI RP Up component through GPS technique appears to be awkward, since the corresponding post fit residuals are considerably larger. Besides, combination tests including multi-temporal TVs display local effects of residual redistribution, when compared to those solutions where Medicina TVs are added one at a time. Finally, the combination of TRFs turns out to be sensitive to the orientation of the local tie into the global frame.     

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.