Showing papers in "Journal of Geodesy in 2011"

ITRF2008: an improved solution of the international terrestrial reference frame

Abstract: ITRF2008 is a refined version of the International Terrestrial Reference Frame based on reprocessed solutions of the four space geodetic techniques: VLBI, SLR, GPS and DORIS, spanning 29, 26, 12.5 and 16 years of observations, respectively. The input data used in its elaboration are time series (weekly from satellite techniques and 24-h session-wise from VLBI) of station positions and daily Earth Orientation Parameters (EOPs). The ITRF2008 origin is defined in such a way that it has zero translations and translation rates with respect to the mean Earth center of mass, averaged by the SLR time series. Its scale is defined by nullifying the scale factor and its rate with respect to the mean of VLBI and SLR long-term solutions as obtained by stacking their respective time series. The scale agreement between these two technique solutions is estimated to be 1.05 ± 0.13 ppb at epoch 2005.0 and 0.049 ± 0.010 ppb/yr. The ITRF2008 orientation (at epoch 2005.0) and its rate are aligned to the ITRF2005 using 179 stations of high geodetic quality. An estimate of the origin components from ITRF2008 to ITRF2005 (both origins are defined by SLR) indicates differences at epoch 2005.0, namely: −0.5, −0.9 and −4.7 mm along X, Y and Z-axis, respectively. The translation rate differences between the two frames are zero for Y and Z, while we observe an X-translation rate of 0.3 mm/yr. The estimated formal errors of these parameters are 0.2 mm and 0.2 mm/yr, respectively. The high level of origin agreement between ITRF2008 and ITRF2005 is an indication of an imprecise ITRF2000 origin that exhibits a Z-translation drift of 1.8 mm/yr with respect to ITRF2005. An evaluation of the ITRF2008 origin accuracy based on the level of its agreement with ITRF2005 is believed to be at the level of 1 cm over the time-span of the SLR observations. Considering the level of scale consistency between VLBI and SLR, the ITRF2008 scale accuracy is evaluated to be at the level of 1.2 ppb (8 mm at the equator) over the common time-span of the observations of both techniques. Although the performance of the ITRF2008 is demonstrated to be higher than ITRF2005, future ITRF improvement resides in improving the consistency between local ties in co-location sites and space geodesy estimates.

First GOCE gravity field models derived by three different approaches

Abstract: Three gravity field models, parameterized in terms of spherical harmonic coefficients, have been computed from 71 days of GOCE (Gravity field and steady-state Ocean Circulation Explorer) orbit and gradiometer data by applying independent gravity field processing methods. These gravity models are one major output of the European Space Agency (ESA) project GOCE High-level Processing Facility (HPF). The processing philosophies and architectures of these three complementary methods are presented and discussed, emphasizing the specific features of the three approaches. The resulting GOCE gravity field models, representing the first models containing the novel measurement type of gravity gradiometry ever computed, are analysed and assessed in detail. Together with the coefficient estimates, full variance-covariance matrices provide error information about the coefficient solutions. A comparison with state-of-the-art GRACE and combined gravity field models reveals the additional contribution of GOCE based on only 71 days of data. Compared with combined gravity field models, large deviations appear in regions where the terrestrial gravity data are known to be of low accuracy. The GOCE performance, assessed against the GRACE-only model ITG-Grace2010s, becomes superior at degree 150, and beyond. GOCE provides significant additional information of the global Earth gravity field, with an accuracy of the 2-month GOCE gravity field models of 10 cm in terms of geoid heights, and 3 mGal in terms of gravity anomalies, globally at a resolution of 100 km (degree/order 200).

The international reference ionosphere today and in the future

Abstract: The international reference ionosphere (IRI) is the internationally recognized and recommended standard for the specification of plasma parameters in Earth’s ionosphere. It describes monthly averages of electron density, electron temperature, ion temperature, ion composition, and several additional parameters in the altitude range from 60 to 1,500 km. A joint working group of the Committee on Space Research (COSPAR) and the International Union of Radio Science (URSI) is in charge of developing and improving the IRI model. As requested by COSPAR and URSI, IRI is an empirical model being based on most of the available and reliable data sources for the ionospheric plasma. The paper describes the latest version of the model and reviews efforts towards future improvements, including the development of new global models for the F2 peak density and height, and a new approach to describe the electron density in the topside and plasmasphere. Our emphasis will be on the electron density because it is the IRI parameter most relevant to geodetic techniques and studies. Annual IRI meetings are the main venue for the discussion of IRI activities, future improvements, and additions to the model. A new special IRI task force activity is focusing on the development of a real-time IRI (RT-IRI) by combining data assimilation techniques with the IRI model. A first RT-IRI task force meeting was held in 2009 in Colorado Springs. We will review the outcome of this meeting and the plans for the future. The IRI homepage is at http://www.IRI.gsfc.nasa.gov.

GOCE gravitational gradiometry

Abstract: GOCE is the first gravitational gradiometry satellite mission. Gravitational gradiometry is the measurement of the second derivatives of the gravitational potential. The nine derivatives form a 3 × 3 matrix, which in geodesy is referred to as Marussi tensor. From the basic properties of the gravitational field, it follows that the matrix is symmetric and trace free. The latter property corresponds to Laplace equation, which gives the theoretical foundation of its representation in terms of spherical harmonic or Fourier series. At the same time, it provides the most powerful quality check of the actual measured gradients. GOCE gradiometry is based on the principle of differential accelerometry. As the satellite carries out a rotational motion in space, the accelerometer differences contain angular effects that must be removed. The GOCE gradiometer provides the components V xx , V yy , V zz and V xz with high precision, while the components V xy and V yz are of low precision, all expressed in the gradiometer reference frame. The best performance is achieved inside the measurement band from 5 × 10–3 to 0.1 Hz. At lower frequencies, the noise increases with 1/f and is superimposed by cyclic distortions, which are modulated from the orbit and attitude motion into the gradient measurements. Global maps with the individual components show typical patterns related to topographic and tectonic features. The maps are separated into those for ascending and those for descending tracks as the components are expressed in the instrument frame. All results are derived from the measurements of the period from November to December 2009. While the components V xx and V yy reach a noise level of about $${10\;\rm{\frac{mE}{\sqrt{Hz}}}}$$ , that of V zz and V xz is about $${20\; \rm{\frac{mE}{\sqrt{Hz}}}}$$ . The cause of the latter’s higher noise is not yet understood. This is also the reason why the deviation from the Laplace condition is at the $${20 \;\rm{\frac{mE}{\sqrt{Hz}}}}$$ level instead of the originally planned $${11\;\rm{\frac{mE}{\sqrt{Hz}}}}$$ . Each additional measurement cycle will improve the accuracy and to a smaller extent also the resolution of the spherical harmonic coefficients derived from the measured gradients.

Mission design, operation and exploitation of the gravity field and steady-state ocean circulation explorer mission

Abstract: The European Space Agency’s Gravity field and steady-state ocean circulation explorer mission (GOCE) was launched on 17 March 2009. As the first of the Earth Explorer family of satellites within the Agency’s Living Planet Programme, it is aiming at a better understanding of the Earth system. The mission objective of GOCE is the determination of the Earth’s gravity field and geoid with high accuracy and maximum spatial resolution. The geoid, combined with the de facto mean ocean surface derived from twenty-odd years of satellite radar altimetry, yields the global dynamic ocean topography. It serves ocean circulation and ocean transport studies and sea level research. GOCE geoid heights allow the conversion of global positioning system (GPS) heights to high precision heights above sea level. Gravity anomalies and also gravity gradients from GOCE are used for gravity-to-density inversion and in particular for studies of the Earth’s lithosphere and upper mantle. GOCE is the first-ever satellite to carry a gravitational gradiometer, and in order to achieve its challenging mission objectives the satellite embarks a number of world-first technologies. In essence the spacecraft together with its sensors can be regarded as a spaceborne gravimeter. In this work, we describe the mission and the way it is operated and exploited in order to make available the best-possible measurements of the Earth gravity field. The main lessons learned from the first 19 months in orbit are also provided, in as far as they affect the quality of the science data products and therefore are of specific interest for GOCE data users.

The ionosphere: effects, GPS modeling and the benefits for space geodetic techniques

Abstract: The main goal of this paper is to provide a summary of our current knowledge of the ionosphere as it relates to space geodetic techniques, especially the most informative technology, global navigation satellite systems (GNSS), specifically the fully deployed and operational global positioning system (GPS). As such, the main relevant modeling points are discussed, and the corresponding results of ionospheric monitoring are related, which were mostly computed using GPS data and based on the direct experience of the authors. We address various phenomena such as horizontal and vertical ionospheric morphology in quiet conditions, traveling ionospheric disturbances, solar flares, ionospheric storms and scintillation. Finally, we also tackle the question of how improved knowledge of ionospheric conditions, especially in terms of an accurate understanding of the distribution of free electrons, can improve space geodetic techniques at different levels, such as higher-order ionospheric effects, precise GNSS navigation, single-antenna GNSS orientation and real-time GNSS meteorology.

Regional reference network augmented precise point positioning for instantaneous ambiguity resolution

Abstract: Integer ambiguity fixing can significantly shorten the initialization time and improve the accuracy of precise point positioning (PPP), but it still takes approximate 15 min of time to achieve reliable integer ambiguity solutions. In this contribution, we present a new strategy to augment PPP estimation with a regional reference network, so that instantaneous ambiguity fixing is achievable for users within the network coverage. In the proposed method, precise zero-differenced atmospheric delays are derived from the PPP fixed solution of the reference stations, which are disseminated to, and interpolated at user stations to correct for L1, L2 phase observations or their combinations. With the corrected observations, instantaneous ambiguity resolution can be carried out within the user PPP software, thus achieving the position solutions equivalent to the network real-time kinematic positioning (NRTK). The strategy is validated experimentally. The derived atmospheric delays and the interpolated corrections are investigated. The ambiguity fixing performance and the resulted position accuracy are assessed. The validation confirms that the new strategy can provide comparable service with NRTK. Therefore, with this new processing strategy, it is possible to integrate PPP and NRTK into a seamless positioning service, which can provide an accuracy of about 10 cm anywhere, and upgrade to a few centimeters within a regional network.

A new automated cycle slip detection and repair method for a single dual-frequency GPS receiver

Abstract: This paper develops a new automated cycle slip detection and repair method that is based on only one single dual-frequency GPS receiver. This method jointly uses the ionospheric total electron contents (TEC) rate (TECR) and Melbourne–Wubbena wide lane (MWWL) linear combination to uniquely determine the cycle slip on both L1 and L2 frequencies. The cycle slips are inferred from the information of ionospheric physical TECR and MWWL ambiguity at the current epoch and that at the previous epoch. The principle of this method is that when there are cycle slips, the MWWL ambiguity will change and the ionospheric TECR will usually be significantly amplified, the part of artificial TECR (caused by cycle slips) being significantly larger than the normal physical TECR. The TECR is calculated based on the dual-frequency carrier phase measurements, and it is highly accurate. We calculate the ionospheric change information (including TECR and TEC acceleration) using the previous epochs (30 epochs in this study) and use the previous data to predict the TECR for the epoch needing cycle slip detection. If the discrepancy is larger than our defined threshold 0.15 TECU/s, cycle slips are regarded to exist at that epoch. The key rational of method is that during a short period (1.0 s in this study) the TECR of physical ionospheric phenomenon will not exceed the threshold. This new algorithm is tested with eight different datasets (including one spaceborne GPS dataset), and the results show that the method can detect and correctly repair almost any cycle slips even under very high level of ionospheric activities (with an average Kp index 7.6 on 31 March 2001). The only exception of a few detected but incorrectly repaired cycle slip is due to a sudden increased pseudorange error on a single satellite (PRN7) under very active ionosphere on 31 March 2001. This method requires dual-frequency carrier phase and pseudorange data from only one single GPS receiver. The other requirement is that the GPS data rate ideally is 1 Hz or higher in order to detect small cycle slips. It is suitable for many applications where one single receiver is used, e.g. real-time kinematic rover station and precise point positioning. An important feature of this method is that it performs cycle slip detection and repair on a satellite-by-satellite basis; thus, the cycle slip detection and repair for each satellite are completely independent and not affected by the data of other satellites.

A new global TEC model for estimating transionospheric radio wave propagation errors

Abstract: Space-based navigation and radar systems operating at single frequencies of <10 GHz require ionospheric corrections of the signal delay or range error. Because this ionospheric propagation error is proportional to the total electron content of the ionosphere along the ray path, a user friendly TEC model covering global scale and all levels of solar activity should be helpful in various applications. Since such a model is not available yet, we present an empirical model approach that allows determining global TEC very easily. Although the number of model coefficients and parameters is rather small, the model describes main ionospheric features with good quality. Presented is the empirical approach describing dependencies on local time, geographic/geomagnetic location and solar irradiance and activity. The non-linear approach needs only 12 coefficients and a few empirically fixed parameters for describing the broad spectrum of TEC variation at all levels of solar activity. The model approach is applied on high-quality global TEC data derived by the Center for Orbit Determination in Europe (CODE) at the University of Berne over more than half a solar cycle (1998–2007). The model fits to these input data with a negative bias of 0.3 TECU and a RMS deviation of 7.5 TECU. As other empirical models too, the proposed Global Neustrelitz TEC Model NTCM-GLis climatological, i.e. the model describes the average behaviour under quiet geomagnetic conditions. During severe space weather events the actual TEC data may deviate from the model values considerably by more than 100%. A preliminary comparison with independent data sets as TOPEX/Poseidon altimeter data reveals similar results for NeQuick and NTCM-GL with RMS deviations in the order of 5 and 11 TECU (1 TECU = 1016 electrons/m2) for low and high-solar activity conditions, respectively. The more extended data base of ionosphere information that accumulates in the coming years will help in further improving the set of coefficients of the model.

Validation of GOCE gravity field models by means of orbit residuals and geoid comparisons

Abstract: Three GOCE-based gravity field solutions have been computed by ESA’s high-level processing facility and were released to the user community. All models are accompanied by variance-covariance information resulting either from the least squares procedure or a Monte-Carlo approach. In order to obtain independent external quality parameters and to assess the current performance of these models, a set of independent tests based on satellite orbit determination and geoid comparisons is applied. Both test methods can be regarded as complementary because they either investigate the performance in the long wavelength spectral domain (orbit determination) or in the spatial domain (geoid comparisons). The test procedure was applied to the three GOCE gravity field solutions and to a number of selected pre-launch models for comparison. Orbit determination results suggest, that a pure GOCE gravity field model does not outperform the multi-year GRACE gravity field solutions. This was expected as GOCE is designed to improve the determination of the medium to high frequencies of the Earth gravity field (in the range of degree and order 50 to 200). Nevertheless, in case of an optimal combination of GOCE and GRACE data, orbit determination results should not deteriorate. So this validation procedure can also be used for testing the optimality of the approach adopted for producing combined GOCE and GRACE models. Results from geoid comparisons indicate that with the 2 months of GOCE data a significant improvement in the determination of the spherical harmonic spectrum of the global gravity field between degree 50 and 200 can be reached. Even though the ultimate mission goal has not yet been reached, especially due to the limited time span of used GOCE data (only 2 months), it was found that existing satellite-only gravity field models, which are based on 7 years of GRACE data, can already be enhanced in terms of spatial resolution. It is expected that with the accumulation of more GOCE data the gravity field model resolution and quality can be further enhanced, and the GOCE mission goal of 1–2 cm geoid accuracy with 100 km spatial resolution can be achieved.

An iterative solution of weighted total least-squares adjustment

Abstract: Total least-squares (TLS) adjustment is used to estimate the parameters in the errors-in-variables (EIV) model. However, its exact solution is rather complicated, and the accuracies of estimated parameters are too difficult to analytically compute. Since the EIV model is essentially a non-linear model, it can be solved according to the theory of non-linear least-squares adjustment. In this contribution, we will propose an iterative method of weighted TLS (WTLS) adjustment to solve EIV model based on Newton–Gauss approach of non-linear weighted least-squares (WLS) adjustment. Then the WLS solution to linearly approximated EIV model is derived and its discrepancy is investigated by comparing with WTLS solution. In addition, a numerical method is developed to compute the unbiased variance component estimate and the covariance matrix of the WTLS estimates. Finally, the real and simulation experiments are implemented to demonstrate the performance and efficiency of the presented iterative method and its linearly approximated version as well as the numerical method. The results show that the proposed iterative method can obtain such good solution as WTLS solution of Schaffrin and Wieser (J Geod 82:415–421, 2008) and the presented numerical method can be reasonably applied to evaluate the accuracy of WTLS solution.

GPS-derived orbits for the GOCE satellite

Abstract: The first ESA (European Space Agency) Earth explorer core mission GOCE (Gravity field and steady-state Ocean Circulation Explorer) was launched on 17 March 2009 into a sun-synchronous dusk–dawn orbit with an exceptionally low initial altitude of about 280 km The onboard 12-channel dual-frequency GPS (Global Positioning System) receiver delivers 1 Hz data, which provides the basis for precise orbit determination (POD) for such a very low orbiting satellite As part of the European GOCE Gravity Consortium the Astronomical Institute of the University of Bern and the Department of Earth Observation and Space Systems are responsible for the orbit determination of the GOCE satellite within the GOCE High-level Processing Facility Both quick-look (rapid) and very precise orbit solutions are produced with typical latencies of 1 day and 2 weeks, respectively This article summarizes the special characteristics of the GOCE GPS data, presents POD results for about 2 months of data, and shows that both latency and accuracy requirements are met Satellite Laser Ranging validation shows that an accuracy of 4 and 7 cm is achieved for the reduced-dynamic and kinematic Rapid Science Orbit solutions, respectively The validation of the reduced-dynamic and kinematic Precise Science Orbit solutions is at a level of about 2 cm

4D GPS water vapor tomography: new parameterized approaches

Abstract: Water vapor is a key variable in numerical weather prediction, as it plays an important role in atmospheric processes. Nonetheless, the distribution of water vapor in the atmosphere is observed with a coarse resolution in time and space compared to the resolution of numerical weather models. GPS water vapor tomography is one of the promising methods to improve the resolution of water vapor measurements. This paper presents new parameterized approaches for the determination of water vapor distribution in the troposphere by GPS. We present the methods and give first results validating the approaches. The parameterization of voxels (volumetric pixels) by trilinear and spline functions in ellipsoidal coordinates are introduced in this study. The evolution in time of the refractivity field is modeled by a Kalman filter with a temporal resolution of 30 s, which corresponds to the available GPS-data rate. The algorithms are tested with simulated and with real data from more than 40 permanent GPS receiver stations in Switzerland and adjoining regions covering alpine areas. The investigations show the potential of the new parameterized approaches to yield superior results compared to the non parametric classical one. The accuracy of the tomographic result is quantified by the inter-quartile range (IQR), which is decreased by 10–20% with the new approaches. Further, parameterized voxel solutions have a substantially smaller maximal error than the non parameterized ones. Simulations show a limited ability to resolve vertical structures above the top station of the network with GPS tomography.

A global mean dynamic topography and ocean circulation estimation using a preliminary GOCE gravity model

Abstract: The Gravity and steady-state Ocean Circulation Explorer (GOCE) satellite mission measures Earth’s gravity field with an unprecedented accuracy at short spatial scales. In doing so, it promises to significantly advance our ability to determine the ocean’s general circulation. In this study, an initial gravity model from GOCE, based on just 2 months of data, is combined with the recent DTU10MSS mean sea surface to construct a global mean dynamic topography (MDT) model. The GOCE MDT clearly displays the gross features of the ocean’s steady-state circulation. More significantly, the improved gravity model provided by the GOCE mission has enhanced the resolution and sharpened the boundaries of those features compared with earlier satellite only solutions. Calculation of the geostrophic surface currents from the MDT reveals improvements for all of the ocean’s major current systems. In the North Atlantic, the Gulf Stream is stronger and more clearly defined, as are the Labrador and the Greenland currents. Furthermore, the finer scale features, such as eddies, meanders and branches of the Gulf Stream and North Atlantic Current system are visible. Similar improvements are seen also in the North Pacific Ocean, where the Kuroshio and its extension are well represented. In the Southern hemisphere, both the Agulhas and the Brazil-Malvinas Confluence current systems are well defined, and in the Southern ocean the Antarctic Circumpolar Current appears enhanced. The results of this preliminary analysis, using an initial GOCE gravity model, clearly demonstrate the potential of the GOCE mission. Already, at this early stage of the mission, the resolution of the MDT has been improved and the estimated surface current speeds have been increased compared with a GRACE satellite-only MDT. Future GOCE gravity models are expected to build further upon this early success.

Evaluation of the first GOCE static gravity field models using terrestrial gravity, vertical deflections and EGM2008 quasigeoid heights

Abstract: Recently, four global geopotential models (GGMs) were computed and released based on the first 2 months of data collected by the Gravity field and steady-state Ocean Circulation Explorer (GOCE) dedicated satellite gravity field mission Given that GOCE is a technologically complex mission and different processing strategies were applied to real space-collected GOCE data for the first time, evaluation of the new models is an important aspect As a first assessment strategy, we use terrestrial gravity data over Switzerland and Australia and astrogeodetic vertical deflections over Europe and Australia as ground-truth data sets for GOCE model evaluation We apply a spectral enhancement method (SEM) to the truncated GOCE GGMs to make their spectral content more comparable with the terrestrial data The SEM utilises the high-degree bands of EGM2008 and residual terrain model data as a data source to widely bridge the spectral gap between the satellite and terrestrial data Analysis of root mean square (RMS) errors is carried out as a function of (i) the GOCE GGM expansion degree and (ii) the four different GOCE GGMs The RMS curves are also compared against those from EGM2008 and GRACE-based GGMs As a second assessment strategy, we compare global grids of GOCE GGM and EGM2008 quasigeoid heights In connection with EGM2008 error estimates, this allows location of regions where GOCE is likely to deliver improved knowledge on the Earth’s gravity field Our ground truth data sets, together with the EGM2008 quasigeoid comparisons, signal clear improvements in the spectral band ~160–165 to ~180–185 in terms of spherical harmonic degrees for the GOCE-based GGMs, fairly independently of the individual GOCE model used The results from both assessments together provide strong evidence that the first 2 months of GOCE observations improve the knowledge of the Earth’s static gravity field at spatial scales between ~125 and ~110 km, particularly over parts of Asia, Africa, South America and Antarctica, in comparison with the pre-GOCE-era

Ionospheric electron density observed by FORMOSAT-3/COSMIC over the European region and validated by ionosonde data

Abstract: This research is motivated by the recent IGS Ionosphere Working Group recommendation issued at the IGS 2010 Workshop held in Newcastle, UK. This recommendation encourages studies on the evaluation of the application of COSMIC radio occultation profiles for additional IGS global ionosphere map (GIM) validation. This is because the reliability of GIMs is crucial to many geodetic applications. On the other hand, radio occultation using GPS signals has been proven to be a promising technique to retrieve accurate profiles of the ionospheric electron density with high vertical resolution on a global scale. However, systematic validation work is still needed before using this powerful technique for sounding the ionosphere on a routine basis. In this paper, we analyze the properties of the ionospheric electron density profiling retrieved from COSMIC radio occultation measurements. A comparison of radio occultation data with ground-based measurements indicates that COSMIC profiles are usually in good agreement with ionosonde profiles, both in the F2 layer peak electron density and the bottom side of the profiles. For this comparison, ionograms recorded by European ionospheric stations (DIAS network) in 2008 were used.

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

Abstract: 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–6 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.

Vertical deformations from homogeneously processed GRACE and global GPS long-term series

Abstract: Temporal variations in the geographic distribution of surface mass cause surface displacements. Surface displacements derived from GRACE gravity field coefficient time series also should be observed in GPS coordinate time series, if both time series are sufficiently free of systematic errors. A successful validation can be an important contribution to climate change research, as the biggest contributors to mass variability in the system Earth include the movement of oceanic, atmospheric, and continental water and ice. In our analysis, we find that if the signals are larger than their precision, both geodetic sensor systems see common signals for almost all the 115 stations surveyed. Almost 80% of the stations have their signal WRMS decreased, when we subtract monthly GRACE surface displacements from those observed by GPS data. Almost all other stations are on ocean islands or small peninsulas, where the physically expected loading signals are very small. For a fair comparison, the data (79 months from September 2002 to April 2009) had to be treated appropriately: the GPS data were completely reprocessed with state-of-the-art models. We used an objective cluster analysis to identify and eliminate stations, where local effects or technical artifacts dominated the signals. In addition, it was necessary for both sets of results to be expressed in equivalent reference frames, meaning that net translations between the GPS and GRACE data sets had to be treated adequately. These data sets are then compared and statistically analyzed: we determine the stability (precision) of GRACE-derived, monthly vertical deformation data to be ~1.2 mm, using the data from three GRACE processing centers. We statistically analyze the mean annual signals, computed from the GPS and GRACE series. There is a detailed discussion of the results for five overall representative stations, in order to help the reader to link the displayed criteria of similarity to real data. A series of tests were performed with the goal of explaining the remaining GPS–GRACE residuals.

Lunar gravity field determination using SELENE same-beam differential VLBI tracking data

Abstract: A lunar gravity field model up to degree and order 100 in spherical harmonics, named SGM 100i, has been determined from SELENE and historical tracking data, with an emphasis on using same-beam S-band differential VLBI data obtained in the SELENE mission between January 2008 and February 2009. Orbit consistency throughout the entire mission period of SELENE as determined from orbit overlaps for the two sub-satellites of SELENE involved in the VLBI tracking improved consistently from several hundreds of metres to several tens of metres by including differential VLBI data. Through orbits that are better determined, the gravity field model is also improved by including these data. Orbit determination performance for the new model shows improvements over earlier 100th degree and order models, especially for edge-on orbits over the deep far side. Lunar Prospector orbit determination shows an improvement of orbit consistency from I-day predictions for 2-day arcs of 6 m in a total sense, with most improvement in the along and cross-track directions. Data fit for the types and satellites involved is also improved. Formal errors for the lower degrees are smaller, and the new model also shows increased correlations with topography over the far side. The estimated value for the lunar GM for this model equals 4902.80080 +/- 0.0009 cu km/sq s (10 sigma). The lunar degree 2 potential Love number k2 was also estimated, and has a value of 0.0255 +/- 0.0016 (10 sigma as well).

GOCE gravitational gradients along the orbit

Abstract: GOCE is ESA’s gravity field mission and the first satellite ever that measures gravitational gradients in space, that is, the second spatial derivatives of the Earth’s gravitational potential. The goal is to determine the Earth’s mean gravitational field with unprecedented accuracy at spatial resolutions down to 100 km. GOCE carries a gravity gradiometer that allows deriving the gravitational gradients with very high precision to achieve this goal. There are two types of GOCE Level 2 gravitational gradients (GGs) along the orbit: the gravitational gradients in the gradiometer reference frame (GRF) and the gravitational gradients in the local north oriented frame (LNOF) derived from the GGs in the GRF by point-wise rotation. Because the V XX , V YY , V ZZ and V XZ are much more accurate than V XY and V YZ , and because the error of the accurate GGs increases for low frequencies, the rotation requires that part of the measured GG signal is replaced by model signal. However, the actual quality of the gradients in GRF and LNOF needs to be assessed. We analysed the outliers in the GGs, validated the GGs in the GRF using independent gravity field information and compared their assessed error with the requirements. In addition, we compared the GGs in the LNOF with state-of-the-art global gravity field models and determined the model contribution to the rotated GGs. We found that the percentage of detected outliers is below 0.1% for all GGs, and external gravity data confirm that the GG scale factors do not differ from one down to the 10?3 level. Furthermore, we found that the error of V XX and V YY is approximately at the level of the requirement on the gravitational gradient trace, whereas the V ZZ error is a factor of 2–3 above the requirement for higher frequencies. We show that the model contribution in the rotated GGs is 2–35% dependent on the gravitational gradient. Finally, we found that GOCE gravitational gradients and gradients derived from EIGEN-5C and EGM2008 are consistent over the oceans, but that over the continents the consistency may be less, especially in areas with poor terrestrial gravity data. All in all, our analyses show that the quality of the GOCE gravitational gradients is good and that with this type of data valuable new gravity field information is obtained.

Global ionosphere maps of VTEC from GNSS, satellite altimetry and FORMOSAT-3/COSMIC data

Abstract: For space geodetic techniques, operating in microwave band, ionosphere is a dispersive medium; thus signals traveling through this medium are in the first approximation affected proportional to inverse of the square of their frequencies. This effect allows gaining information about the parameters of the ionosphere in terms of Total Electron Content (TEC) or the electron density (Ne). TEC or electron density can then be expressed by means of spherical harmonic base functions to provide a Global Ionosphere Map (GIM). The classical input data for development of GIMs are obtained from dual-frequency observations carried out at Global Navigation Satellite Systems (GNSS) stations. However, GNSS stations are in-homogeneously distributed around the world, with large gaps particularly over the oceans; this fact reduces the precision of the GIM over these areas. On the other hand, dual-frequency satellite altimetry missions such as Jason-1 provide information about the ionosphere precisely above the oceans; and furthermore Low Earth Orbiting (LEO) satellites, such as Formosat-3/COSMIC (F/C) provide well-distributed information of ionosphere globally. This study investigates on global modeling of TEC through combining GNSS and satellite altimetry data with global TEC data derived from the occultation measurements of the F/C mission. The combined GIMs of vertical TEC (VTEC) show a maximum difference of 1.3–1.7 TEC units (TECU) with respect to the GNSS-only GIMs in the whole day. The root mean square error (RMS) maps of combined solution show a reduction of about 0.1 TECU in the whole day. This decrease of RMS can reach up to 0.5 TECU in areas where no or few GNSS observations are available, but high number of F/C measurement is carried out. This proves that the combined GIMs provide a more homogeneous global coverage and higher reliability than results of each single method. All comparisons and validations made within this study provide vital information regarding combination and integration of various observation techniques in the Global Geodetic Observing System of the International Association of Geodesy.

Monte Carlo simulations of the impact of troposphere, clock and measurement errors on the repeatability of VLBI positions

Abstract: Within the International VLBI Service for Geodesy and Astrometry (IVS) Monte Carlo simulations have been carried out to design the next generation VLBI system (“VLBI2010”). Simulated VLBI observables were generated taking into account the three most important stochastic error sources in VLBI, i.e. wet troposphere delay, station clock, and measurement error. Based on realistic physical properties of the troposphere and clocks we ran simulations to investigate the influence of the troposphere on VLBI analyses, and to gain information about the role of clock performance and measurement errors of the receiving system in the process of reaching VLBI2010’s goal of mm position accuracy on a global scale. Our simulations confirm that the wet troposphere delay is the most important of these three error sources. We did not observe significant improvement of geodetic parameters if the clocks were simulated with an Allan standard deviation better than 1 × 10−14 at 50 min and found the impact of measurement errors to be relatively small compared with the impact of the troposphere. Along with simulations to test different network sizes, scheduling strategies, and antenna slew rates these studies were used as a basis for the definition and specification of VLBI2010 antennas and recording system and might also be an example for other space geodetic techniques.

Combination of GNSS and SLR observations using satellite co-locations

Abstract: Satellite Laser Ranging (SLR) observations to Global Navigation Satellite System (GNSS) satellites may be used for several purposes. On one hand, the range measurement may be used as an independent validation for satellite orbits derived solely from GNSS microwave observations. On the other hand, both observation types may be analyzed together to generate a combined orbit. The latter procedure implies that one common set of orbit parameters is estimated from GNSS and SLR data. We performed such a combined processing of GNSS and SLR using the data of the year 2008. During this period, two GPS and four GLONASS satellites could be used as satellite co-locations. We focus on the general procedure for this type of combined processing and the impact on the terrestrial reference frame (including scale and geocenter), the GNSS satellite antenna offsets (SAO) and the SLR range biases. We show that the combination using only satellite co-locations as connection between GNSS and SLR is possible and allows the estimation of SLR station coordinates at the level of 1–2 cm. The SLR observations to GNSS satellites provide the scale allowing the estimation of GNSS SAO without relying on the scale of any a priori terrestrial reference frame. We show that the necessity to estimate SLR range biases does not prohibit the estimation of GNSS SAO. A good distribution of SLR observations allows a common estimation of the two parameter types. The estimated corrections for the GNSS SAO are 119 mm and −13 mm on average for the GPS and GLONASS satellites, respectively. The resulting SLR range biases suggest that it might be sufficient to estimate one parameter per station representing a range bias common to all GNSS satellites. The estimated biases are in the range of a few centimeters up to 5 cm. Scale differences of 0.9 ppb are seen between GNSS and SLR.

Combination of different space-geodetic observations for regional ionosphere modeling

Abstract: Most of the space-geodetic observation techniques can be used for modeling the distribution of free electrons in the Earth’s ionosphere. By combining different techniques one can take advantage of their different spatial and temporal distributions as well as their different observation characteristics and sensitivities concerning ionospheric parameter estimation. The present publication introduces a procedure for multi-dimensional ionospheric modeling. The model consists of a given reference part and an unknown correction part expanded in terms of B-spline functions. This approach is used to compute regional models of Vertical Total Electron Content (VTEC) based on the International Reference Ionosphere (IRI 2007) and GPS observations from terrestrial Global Navigation Satellite System (GNSS) reference stations, radio occultation data from Low Earth Orbiters (LEOs), dual-frequency radar altimetry measurements, and data obtained by Very Long Baseline Interferometry (VLBI). The approach overcomes deficiencies in the climatological IRI model and reaches the same level of accuracy than GNSS-based VTEC maps from IGS. In areas without GNSS observations (e.g., over the oceans) radio occultations and altimetry provide valuable measurements and further improve the VTEC maps. Moreover, the approach supplies information on the offsets between different observation techniques as well as on their different sensitivity for ionosphere modeling. Altogether, the present procedure helps to derive improved ionospheric corrections (e.g., for one-frequency radar altimeters) and at the same time it improves our knowledge on the Earth’s ionosphere.

The AUSGeoid09 model of the Australian Height Datum

Abstract: AUSGeoid09 is the new Australia-wide gravimetric quasigeoid model that has been a posteriori fitted to the Australian Height Datum (AHD) so as to provide a product that is practically useful for the more direct determination of AHD heights from Global Navigation Satellite Systems (GNSS). This approach is necessary because the AHD is predominantly a third-order vertical datum that contains a ~1 m north-south tilt and ~0.5 m regional distortions with respect to the quasigeoid, meaning that GNSS-gravimetric-quasigeoid and AHD heights are inconsistent. Because the AHD remains the official vertical datum in Australia, it is necessary to provide GNSS users with effective means of recovering AHD heights. The gravimetric component of the quasigeoid model was computed using a hybrid of the remove-compute-restore technique with a degree-40 deterministically modified kernel over a one-degree spherical cap, which is superior to the remove-compute-restore technique alone in Australia (with or without a cap). This is because the modified kernel and cap combine to filter long-wavelength errors from the terrestrial gravity anomalies. The zero-tide EGM2008 global gravitational model to degree 2,190 was used as the reference field. Other input data are ~1.4 million land gravity anomalies from Geoscience Australia, 1′ × 1′ DNSC2008GRA altimeter-derived gravity anomalies offshore, the 9′′ × 9′′ GEODATA-DEM9S Australian digital elevation model, and a readjustment of Australian National Levelling Network (ANLN) constrained to the CARS2006 mean dynamic ocean topography model. To determine the numerical integration parameters for the modified kernel, the gravimetric component of AUSGeoid09 was compared with 911 GNSS-observed ellipsoidal heights at benchmarks. The standard deviation of fit to the GNSS-AHD heights is ±222 mm, which dropped to ±134 mm for the readjusted GNSS-ANLN heights showing that careful consideration now needs to be given to the quality of the levelling data used to assess gravimetric quasigeoid models. The publicly released version of AUSGeoid09 also includes a geometric component that models the difference between the gravimetric quasigeoid and the zero surface of the AHD at 6,794 benchmarks. This a posteriori fitting used least-squares collocation (LSC) in cross-validation mode to determine a correlation length of 75 km for the analytical covariance function, whereas the noise was taken from the estimated standard deviation of the GNSS ellipsoidal heights. After this LSC surface fitting, the standard deviation of fit reduced to ±30 mm, one-third of which is attributable to the uncertainty in the GNSS ellipsoidal heights.

Transverse Mercator with an accuracy of a few nanometers

Abstract: Implementations of two algorithms for the transverse Mercator projection are described; these achieve accuracies close to machine precision. One is based on the exact equations of Thompson and Lee and the other uses an extension of Kruger’s series for the mapping to higher order. The exact method provides an accuracy of 9 nm over the entire ellipsoid, while the errors in the series method are less than 5 nm within 3900 km of the central meridian. In each case, the meridian convergence and scale are also computed with similar accuracy. The speed of the series method is competitive with other less accurate algorithms and the exact method is about five times slower.

Global sea-level rise and its relation to the terrestrial reference frame

Abstract: We examined the sensitivity of estimates of global sea-level rise obtained from GPS-corrected long term tide gauge records to uncertainties in the International Terrestrial Reference Frame (ITRF) realization. A useful transfer function was established, linking potential errors in the reference frame datum (origin and scale) to resulting errors in the estimate of global sea level rise. Contrary to scale errors that are propagated by a factor of 100%, the impact of errors in the origin depends on the network geometry. The geometry of the network analyzed here resulted in an error propagation factor of 50% for the Z component of the origin, mainly due to the asymmetry in the distribution of the stations between hemispheres. This factor decreased from 50% to less than 10% as the geometry of the network improved using realistic potential stations that did not yet meet the selection criteria (e.g., record length, data availability). Conversely, we explored new constraints on the reference frame by considering forward calculations involving tide gauge records. A reference frame could be found in which the scatter of the regional sea-level rates was limited. The resulting reference frame drifted by 1.36 ± 0.22 mm/year from the ITRF2000 origin in the Z component and by −0.44 ± 0.22 mm/year from the ITRF2005 origin. A bound on the rate of global sea level rise of 1.2 to 1.6 mm/year was derived for the past century, depending on the origin of the adopted reference frame. The upper bound is slightly lower than previous estimates of 1.8 mm/year discussed in the IPCC fourth report.

Precise orbit determination of GIOVE-B based on the CONGO network

Abstract: GIOVE-B is one of two test satellites for the future European Global Navigation Satellite System Galileo. The Cooperative Network for GIOVE Observation (CONGO) is a global tracking network of GIOVE-capable receivers established by Deutsches Zentrum fur Luft- und Raumfahrt (DLR) and Bundesamt fur Kartographie und Geodasie (BKG). This network provides the basis for the precise orbit determination of the GIOVE-B satellite for the time period 29 June till 27 October 2009 with a modified version of the Bernese GPS Software. Different arc lengths and sets of orbit parameters were tested. These tests showed that the full set of nine radiation pressure parameters resulted in a better performance than the reduced set of five parameters. An internal precision of about one to two decimeters could be demonstrated for the central day of 5-day solutions. The orbit predictions have a precision of about 1 m for a prediction period of 24 h. External validations with Satellite Laser Ranging (SLR) show residuals on the level of 12 cm. The accuracy of the final orbits is expected to be on the few decimeter level.

Assessing Greenland ice mass loss by means of point-mass modeling: a viable methodology

Abstract: Greenland ice mass loss is one of the most serious phenomena of present-day global climate change. In this context, both the quantification of overall deglaciation rates and its spatial localization are highly significant. We have thoroughly investigated the technique of point-mass modeling in order to derive mass-balance patterns from GRACE (Gravity Recovery And Climate Experiment) gravimetry. The method infers mass variations on the Earth’s surface from gravitational signals at satellite altitude. In order to solve for point-mass changes, we applied least-squares adjustment. Due to downward continuation, numerical stabilization of the inversion process gains particular significance. We stabilized the ill-posed problem by Tikhonov regularization. Our simulation and real data experiments show that point-mass modeling provides both rational deglaciation rates and high-resolution spatial mass variation patterns.

Evaluation of the impact of atmospheric pressure loading modeling on GNSS data analysis

Abstract: In recent years, several studies have demonstrated the sensitivity of Global Navigation Satellite System (GNSS) station time series to displacements caused by atmospheric pressure loading (APL). Different methods to take the APL effect into account are used in these studies: applying the corrections from a geophysical model on weekly mean estimates of station coordinates, using observation-level corrections during data analysis, or solving for regression factors between the station displacement and the local pressure. The Center for Orbit Determination in Europe (CODE) is one of the global analysis centers of the International GNSS Service (IGS). The current quality of the IGS products urgently asks to consider this effect in the regular processing scheme. However, the resulting requirements for an APL model are demanding with respect to quality, latency, and—regarding the reprocessing activities—availability over a long time interval (at least from 1994 onward). The APL model of Petrov and Boy (J Geophys Res 109:B03405, 2004) is widely used within the VLBI community and is evaluated in this study with respect to these criteria. The reprocessing effort of CODE provides the basis for validating the APL model. The data set is used to solve for scaling factors for each station to evaluate the geophysical atmospheric non-tidal loading model. A consistent long-term validation of the model over 15 years, from 1994 to 2008, is thus possible. The time series of 15 years allows to study seasonal variations of the scaling factors using the dense GNSS tracking network of the IGS. By interpreting the scaling factors for the stations of the IGS network, the model by (2004) is shown to meet the expectations concerning the order of magnitude of the effect at individual stations within the uncertainty given by the GNSS data processing and within the limitations due to the model itself. The repeatability of station coordinates improves by 20% when applying the effect directly on the data analysis and by 10% when applying a post-processing correction to the resulting weekly coordinates compared with a solution without taking APL into account.