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–7 mm for the positions and 1.0–1.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.

Consistent and homogeneous long-time series of the space geodetic techniquesGlobal Positioning System (GPS), Satellite Laser Ranging (SLR), andVery Long Baseline Interferometry (VLBI) provide the basis for thecombination efforts of GGOS-D. For a consistent combination, thedefinition of common standards for parameterization and modeling isessential. These standards and the technique-specific processingoptions of all individual GPS, SLR, and VLBI solutions as well asthe combined SLR and VLBI solutions are discussed.

The GGOS-D global terrestrial reference frame has been computed from a combination of homogeneously processed VLBI, SLR and GPS observation time series. A major focus was on the analysis of station position time series, investigations regarding the handling of non-linear station motions, and the development of refined combination methods. The terrestrial reference frame (station positions and velocities) has been estimated simultaneously with the Earth orientation parameters and the celestial reference frame (quasar coordinates).

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.

In its function as an ITRS Combination Centre DGFI has computed a solution of the International Terrestrial Reference Frame 2005 (ITRF2005). It is based on the combination of epoch normal equations (weekly or session data sets, respectively) of station positions and Earth Orientation Parameters (EOPs) from the geodetic space techniques-specific normal equations, the inter-technique combination using local tie measurements at co-location sites, and the computation of the ITRF2005 solution.

ITRF2005 is the first terrestrial reference frame computed from weekly/ session-wise data sets of the geodetic space techniques GPS, VLBI, SLR and DORIS. This allows to detect time variable effects in station positions, such as discontinuities or seasonal variations. In the combination of the terrestrial reference frame these time variable effects need to be taken into account.As part of the reference frame computation with the ITRF2005 data sets at DGFI, we have computed multi-year solutions of each technique with the combination software DOGS-CS. We compared individual weekly data sets to the multi-year solution by aligning them with a 7 parameter similarity transformation and analyzed the resulting time series of station positions and transformation parameters. This is done with respect to discontinuities (caused e. g. by instrumental changes or earthquakes) or periodic signals such as annual variations. This information is then used to compute a second improved iteration of technique multi-year solutions, where these effects are taken into account.

A. NOTHNAGEL, M. ROTHACHER, D. ANGERMANN, T. ARTZ, S. BOCKMANN, W. BOSCH, H. DREWES, M. GERSTL, R. KELM, M. KRUGEL, D. KONIG, R. KONIG, B. MEISEL, H. MULLER, B. RICHTER, N. PANAFIDINA, S. RUDENKO, W. SCHWEGMANN, P. STEIGENBERGER, V. TESMER, D. THALLER 1 Institut fur Geodasie und Geoinformation der Universitat Bonn (IGGB) Nusallee 17, D-53115 Bonn, Germany e-mail: nothnagel@uni-bonn.de 2 GeoForschungsZentrum (GFZ) Telegrafenberg D-14473 Potsdam, Germany 3 Deutsches Geodatisches Forschungsinstitut (DGFI) Alfons-Goppel-Strase 11, D-80539 Munchen 4 Bundesamt fur Kartographie und Geodasie (BKG) Richard-Straus-Allee 11, D-60598 Frankfurt am Main, Germany

In its function as an ITRS Coinbination Center, DGFI has developed refined methods for the terrestrial reference frame computation. The advanced approach is based on the combination of epoch normal equations (weekly/daily data sets) containing station positions and Earth orientation parameters (EOP) obtained from different geodetic space techniques such as VLBI, SLR, GPS and DORIS. This refined approach allows to account for nonlinear effects (e.g., periodic signals and discontinuities) in station positions and to ensure consistency between TRF and EOP. The ITRF2004 submissions were used as input for a refined realization of the terrestrial reference frame. This paper presents the combination methodology and the current status of the ITRF2004 computations at DGFI.