Search

Traditionally, sea level is observed at tide gauge stations, which usually also serve as height reference stations for national levelling networks and therefore define a height system of a country. Thus, sea level research across countries is closely linked to height system unification and needs to be regarded jointly. The project aims to make use of a new observation technique, namely SAR positioning, which can help to connect the GNSS basic network of a country to tide gauge stations and as such to link the sea level records of tide gauge stations to the geometric network. By knowing the geoid heights at the tide gauge stations in a global height reference frame with high precision, one can finally obtain absolute sea level heights of the tide gauge stations in a common reference system and can link them together. By this method, on the one hand national height systems can be connected and on the other hand the absolute sea level at the tide gauge stations can be determined. By analyzing time series of absolute sea level heights their changes can be determined in an absolute sense in a global reference frame and the impact of climate change on sea level can be quantified (e.g. by ice sheet and glacier melting, water inflow, global warming). The paper presents the main scientific questions to be addressed by the project, introduces the idea of using SAR transponders for this application and describes the observation network implemented for this feasibility study.

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.

As early as the 15th century Swedes noticed that rocks in their harbors were slowly rising out of the sea [Ekman, 1991]. These local observations were not sufficient to distinguish whether the rocks were rising or the sea level falling. Later, it was realized that Fennoscandia was still rebounding from the last Ice Age. This historical observation is still relevant today. How can you know whether a point on the Earth's surface is slowly moving up, down, or horizontally? One must relate local measurements to a stable and accurate reference frame, one whose scale is much larger than the problem at hand. We remain concerned with sea-level variations, but present-day studies recognize that change must be measured from a global point of view and with respect to a globally well-defined reference frame. Thus, the regional and national geodetic datums developed over the past 200 years are inappropriate for studying global-level problems.

Repeated GEODYSSEA-GPS measurements and additional GPS data from the IGS and APRGP97/98 constrain the motion of SE-Asia within a global reference frame. Sundaland i.e. Indochina as well as the western and central part of Indonesia, together with South-China, constitute an apparently stable tectonic block that is decoupled from Eurasia. In the ITRF97 this block moves to the east about an Euler pole approximating that of Eurasia, but with a velocity that is one third higher than the eastward directed movement of Eurasia. With respect to India and Australia the Sundaland-South China block is moving due south. This suggests that a) Sundaland-South China are moving coherently to the East along the boundaries studied and b) the eastward motion of India is compensated by the eastward motion of Sundaland-South China. The current rather homogeneous kinematic behaviour of Sundaland-South China that differs in rate (no more than 5 mm) rather than direction from the movement of N-China differs from the kinematics proposed by different modelling approaches throughout the literature. The data help to constrain locations and behaviour of the active fault zones and give information on the interplate and intraplate deformation in the area.

This paper presents the final geodetic results of the GEODYSSEA project. The GPS data from a 42 station network observed during two field campaigns (1994/1996) were analyzed by four groups using different software packages and analysis strategies. The precision of both campaign coordinate solutions was found to be 4–7 mm for the horizontal, and 1 cm for the vertical component. The campaign solutions were merged into one unique solution, which was accurately mapped into the ITRF-96 reference frame. The global accuracy of this solution with respect to ITRF-96 is ±1 cm, while the resolution of the relative horizontal velocities is estimated to be at the level of 2–3 mm/yr. This solution was used as the basis for all scientific interpretations, which are published in separate papers. The velocity estimates of a part of the network provided the first direct measurement of a relative motion of the Sundaland block with respect to Eurasian plate.

Geodetic reference systems are fundamental requisites for accurate and reliable geodetic results. Unambiguous reference systems are needed to refer the geodetic observations and estimated parameters to a unique global basis. Highly accurate, consistent and reliable realizations of the terrestrial reference systems are required for measuring and mapping the Earth’s surface and its variations in time. These terrestrial reference frames are the basis for many practical applications, such as national and regional geodetic networks, engineering, precise navigation, geo-information systems, etc. as well as for scientific investigations in the Earth’s system (e.g., tectonic plate motion, sea level change, seasonal and secular loading signals, atmosphere dynamics and Earth orientation excitation).

The GGOS-D terrestrial reference frame has been computed in a common adjustment of station positions and velocities together with the Earth orientation parameters and the quasar coordinates (celestial reference frame). The data were processed as datum-free normal equations from homogeneously generated VLBI, SLR and GPS observation time series using identical standards for the modelling and parameterization. A major focus was on the analysis of the station position time series, investigations regarding seasonal variations in station motions and on the combination methodology for the terrestrial reference frame computation.

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).

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.

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.

This paper investigates the time evolution of the terrestrial reference frame (TRF). We analysed time series of site positions and datum parameters obtained from VLBI, SLR, GPS and DORIS solutions with respect to non-linear motions (e.g. jumps, seasonal signals) and systematic differences. The time series of the translation and scale parameters were used to identify solution- and technique-related problems, and to investigate the contribution of the different techniques to realise the TRF origin and scale. The position time series reveal non-linear motions and jumps for many sites, which have to be considered for future TRF realisations.

Results from repeated measurements across a 90-sites global positioning system network in Central Asia provide direct evidence of current high rates of tectonic deformation far north of the India–Eurasia suture zone. The data help to quantify both the partitioning of deformation within the seismically active Tien Shan and Northern Pamirs, as well as the ongoing rotation of virtually undeformed blocks such as the Tarim with respect to stable Eurasia. The NNW to SSE shortening rates derived from the computations reach values of 23±3 mm yr −1 over the ranges studied, suggesting that regional deformation rates have increased considerably within the study area since the onset of Cenozoic shortening 20–25 Myr ago.

With the advent of global positioning system (GPS) technology and its applications to high precision geodesy, earth scientists have gained an unprecedented opportunity to study the kinematics and dynamics of present-day deformation processes along plate boundary zones. However, until now our knowledge of the deformation processes associated with the Andean subduction along the western coast of South America has been limited. Here we present first estimates of present-day crustal deformation rates in the central and southern sections of the Andes, between latitudes 22‡S and 42‡S. We find that the deformation in the central section of our study area is dominated by the interseismic phase of an earthquake deformation cycle, caused by 100% coupling of the thrust interface between the subducting Nazca and the overriding South America plates. The estimated depth of coupling is not uniform along strike: north of 30‡S it is V33 km deep, while south of 35‡S it reaches 50 km depth. In addition, we observed postseismic relaxation effects in the northern part of our network in the area of the 1995 Mw8.0 Antofagasta earthquake. South of 38‡S, we detected a similar deformation pattern as in the Antofagasta area which we attribute to postseismic relaxation effects of the 1960 Mw9.5 Valdivia earthquake. fl 2001 Elsevier Science B.V. All rights reserved.

GPS data from four sites in the Nazca Plate (Easter Island, Galapagos, Robinson Crusoe and San Felix Islands) and from five sites in the stable core of the South American Plate enabled us to estimate the Euler vector of the Nazca Plate with respect to South America. The observed velocities of Easter Island (6.6 cm/yr at 102.3< degrees >), Galapagos (5.1 cm/yr at 90.0< degrees >), Robinson Crusoe (6.6 cm/yr at 80.1< degrees >) and San Felix (6.0 cm/yr at 82.1< degrees >) are significantly slower than the global plate model NUVEL-1A velocities for those four sites. Our estimated Euler pole is located at 48.8< degrees >N, 91.7< degrees >W with a rate of 0.59< degrees >/m.y. The NUVEL-1A and earlier space-geodetic studies give rotation rates that are 20% faster.

In order to study both the interplate seismic loading cycle and the distribution of intraplate deformation of the Andes, a 215 site GPS network covering Chile and the western part of Argentina was selected, monumented and observed in 1993 and 1994. A dense part of the network in northern Chile and northwest Argentina, comprising some 70 sites, was re-observed after two years in October/November, 1995. The M w = 8.0 Antofagasta (North Chile) earthquake of 30th July, 1995 took place between the two observations. The city of Antofagasta shifted 80 cm westwards by this event and the displacement still reached 10 cm at locations 300 km from the trench. Three different deformation processes have been considered for modeling the measured displacements: (1) interseismic accumulation of elastic strain due to subduction coupling, (2) coseismic strain release during the Antofagasta earthquake and (3) crustal shortening in the Sub-Andes.

In the framework of the GEODYnamics of South and South-East Asia (GEODYSSEA) project a network of 42 stations was observed by GPS during two field campaigns in 1994 and 1996. BKG and GFZ realized a third observation of the complete network in 1998. The data was analyzed independently by two analysis centers and a precision of the coordinate solutions was found to be 4–7 mm for the horizontal, and 10 mm for the vertical component. Subsequently, these campaign solutions were merged into one unique solution, which was accurately mapped into the International Terrestrial Reference Frame of 1997 (ITRF97). The global accuracy of this solution with respect to ITRF97 is of the order of 10 mm, while the resolution of the relative horizontal velocities is estimated to be at the level of 2–3 mm/yr. Comparison of the new solution with the previously published GEODYSSEA solution indicates, that improved results for the motion of Sundaland, South China and the overall deformation in the area could be found. Improvements are due to the increased time span and due to a better connection to and the advances in the ITRF reference frame. The new set of coordinates and velocities is used as the basis for scientific interpretations, which have been and will be published in separate papers.