A first look at the effects of ionospheric signal bending on a globally processed GPS network

This study provides a first attempt at quantifying potential signal bending effects on the GPS reference frame, coordinates and zenith tropospheric delays (ZTDs). To do this, we homogeneously reanalysed data from a global network of GPS sites spanning 14 years (1995.0–2009.0). Satellite, Earth orientation, tropospheric and ground station coordinate parameters were all estimated. We tested the effect of geometric bending and dTEC bending corrections, which were modelled at the observation level based, in part, on parameters from the International Reference Ionosphere 2007 model. Combined, the two bending corrections appear to have a minimal effect on site coordinates and ZTDs except for low latitude sites. Considering five days (DOY 301–305, 28 October–1 November 2001) near ionospheric maximum in detail, they affect mean ZTDs by up to ~1.7 mm at low latitudes, reducing to negligible levels at high latitudes. Examining the effect on coordinates in terms of power-spectra revealed the difference to be almost entirely white noise, with noise amplitude ranging from 0.3 mm (high latitudes) to 2.4 mm (low latitudes). The limited effect on station coordinates is probably due to the similarity in the elevation dependence of the bending term with that of tropospheric mapping functions. The smoothed z-translation from the GPS reference frame to ITRF2005 changes by less than 2 mm, though the effect combines positively with that from the second order ionospheric refractive index term. We conclude that, at the present time, and for most practical purposes, the geometric and dTEC bending corrections are probably negligible at current GPS/reference frame precisions.     

Combination of GNSS and SLR observations using satellite co-locations

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.     

Present-day crustal deformation in China relative to ITRF97 kinematic plate model

 A global plate motion model is established based on the ITRF97 velocity fields and geological model NUVEL1. Sub-plate models are estimated by using the velocity fields derived from 45 global positioning system (GPS) sites under the ITRF97 reference frame in China. Comparisons between space geodesy and geological models are given. It is found that the Euler vector of the AFRC–EURA pair has an obvious discrepancy between space geodetic and geological models. The motion patterns of tectonic blocks in China predicted by GPS are consistent with those of geological data on the whole. Received: 9 November 2000 / Accepted: 17 September 2001     

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.     

GPS network monitors the Western Alps' deformation over a five-year period: 1993–1998

  The Western Alps are among the best studied collisional belts with both detailed structural mapping and also crustal geophysical investigations such as the ECORS and EGT seismic profile. By contrast, the present-day kinematics of the belt is still largely unknown due to small relative motions and the insufficient accuracy of the triangulation data. As a consequence, several tectonic problems still remain to be solved, such as the amount of N–S convergence in the Occidental Alps, the repartition of the deformation between the Alpine tectonic units, and the relation between deformation and rotation across the Alpine arc. In order to address these problems, the GPS ALPES group, made up of French, Swiss and Italian research organizations, has achieved the first large-scale GPS surveys of the Western Alps. More than 60 sites were surveyed in 1993 and 1998 with a minimum observation of 3 days at each site. GPS data processing has been done by three independent teams using different software. The different solutions have horizontal repeatabilities (N–E) of 4–7 mm in 1993 and 2–3 mm in 1998 and compare at the 3–5-mm level in position and 2-mm/yr level in velocity. A comparison of 1993 and 1998 coordinates shows that residual velocities of the GPS marks are generally smaller than 2 mm/yr, precluding a detailed tectonic interpretation of the differential motions. However, these data seem to suggest that the N–S compression of the Western Alps is quite mild (less than 2 mm/yr) compared to the global convergence between the African and Eurasian plate (6 mm/yr). This implies that the shortening must be accomodated elsewhere by the deformation of the Maghrebids and/or by rotations of Mediterranean microplates. Also, E–W velocity components analysis supports the idea that E–W extension exists, as already suggested by recent structural and seismotectonic data interpretation. Received: 27 November 2000 / Accepted: 17 September 2001     

Tectonic plate motion and co-seismic steps surveyed by DORIS field beacons: the New Hebrides experiment

 The New Hebrides experiment consisted of setting up a pair of DORIS beacons in remote tropical islands in the southwestern Pacific, between 1993 and 1997. Because of orbitography requirements on TOPEX/Poséidon, the beacons were only transmitting to SPOT satellites. Root-mean-square (RMS) scatters at the centimeter level on the latitude and vertical components were achieved, but 2-cm RMS scatters affected the longitude component. Nevertheless, results of relative velocity (123 mm/year N250°) are very consistent with those obtained using the global positioning system (GPS) (126 mm/yr N246°). The co-seismic step (12 mm N60°) related to the Walpole event (M _W = 7.7) is consistent with that derived from GPS (10 mm N30°) or from the centroid moment tensor (CMT) of the quake (12 mm N000°). Received: 19 November 1999 / Accepted: 17 May 2000     

Estimates of ocean tide loading displacements and its impact on position time series in Hong Kong using a dense continuous GPS network

Three-dimensional ocean tide loading (OTL) displacements of eight diurnal and semidiurnal constituents at 12 sites in Hong Kong were estimated using 3–7 years of continuous global positioning system (GPS) observations. OTL displacements were estimated using the precise point positioning (PPP) technique on a daily basis and then combined. The OTL displacements obtained by GPS were compared with predictions using seven recent global ocean tide models. The effect of OTL displacements on GPS position time series was also investigated. The study shows that the GPS-derived OTL displacements (excluding K1 and K2 constituents) agree best with those predicted by the GOT4.7 and NAO99b models. The GPS/model agreement is generally at the sub-millimeter level, except for S2, K1, and K2 constituents with relatively large errors. After systematic biases between the GPS and model values are removed, the misfits of all sites for M2, S2, N2, O1, P1, and Q1 are less than 0.5 and 1.0 mm in the horizontal and vertical components, respectively, while larger misfits (within 2.5 mm) are observed for K1 and K2. Integer ambiguity fixing slightly improves the east component of OTL displacement estimates. The study also finds that GPS-derived OTL corrections, instead of model predicts, can be used in daily data processing with the exception of K1 and K2. Including K2 corrections, a secular vertical rate of up to 1 mm/year in position time series can be induced, which needs to be confirmed by further studies.     

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

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.     

Single receiver phase ambiguity resolution with GPS data

Global positioning system (GPS) data processing algorithms typically improve positioning solution accuracy by fixing double-differenced phase bias ambiguities to integer values. These “double-difference ambiguity resolution” methods usually invoke linear combinations of GPS carrier phase bias estimates from pairs of transmitters and pairs of receivers, and traditionally require simultaneous measurements from at least two receivers. However, many GPS users point position a single local receiver, based on publicly available solutions for GPS orbits and clocks. These users cannot form double differences. We present an ambiguity resolution algorithm that improves solution accuracy for single receiver point-positioning users. The algorithm processes dual- frequency GPS data from a single receiver together with wide-lane and phase bias estimates from the global network of GPS receivers that were used to generate the orbit and clock solutions for the GPS satellites. We constrain (rather than fix) linear combinations of local phase biases to improve compatibility with global phase bias estimates. For this precise point positioning, no other receiver data are required. When tested, our algorithm significantly improved repeatability of daily estimates of ground receiver positions, most notably in the east component by approximately 30% with respect to the nominal case wherein the carrier biases are estimated as real values. In this “static” test for terrestrial receiver positions, we achieved daily repeatability of 1.9, 2.1 and 6.0 mm in the east, north and vertical (ENV) components, respectively. For kinematic solutions, ENV repeatability is 7.7, 8.4, and 11.7 mm, respectively, representing improvements of 22, 8, and 14% with respect to the nominal. Results from precise orbit determination of the twin GRACE satellites demonstrated that the inter-satellite baseline accuracy improved by a factor of three, from 6 to 2 mm up to a long-term bias. Jason-2/Ocean Surface Topography Mission precise orbit determination tests results implied radial orbit accuracy significantly below the 10 mm level. Stability of time transfer, in low-Earth orbit, improved from 40 to 7 ps. We produced these results by applying this algorithm within the Jet Propulsion Laboratory’s (JPL’s) GIPSY/OASIS software package and using JPL’s orbit and clock products for the GPS constellation. These products now include a record of the wide-lane and phase bias estimates from the underlying global network of GPS stations. This implies that all GIPSY–OASIS positioning users can now benefit from this capability to perform single-receiver ambiguity resolution.     

Impact of Earth radiation pressure on GPS position estimates

GPS satellite orbits available from the International GNSS Service (IGS) show a consistent radial bias of up to several cm and a particular pattern in the Satellite Laser Ranging (SLR) residuals, which are suggested to be related to radiation pressure mismodeling. In addition, orbit-related frequencies were identified in geodetic time series such as apparent geocenter motion and station displacements derived from GPS tracking data. A potential solution to these discrepancies is the inclusion of Earth radiation pressure (visible and infrared) modeling in the orbit determination process. This is currently not yet considered by all analysis centers contributing to the IGS final orbits. The acceleration, accounting for Earth radiation and satellite models, is introduced in this paper in the computation of a global GPS network (around 200 IGS sites) adopting the analysis strategies from the Center for Orbit Determination in Europe (CODE). Two solutions covering 9 years (2000–2008) with and without Earth radiation pressure were computed and form the basis for this study. In previous studies, it has been shown that Earth radiation pressure has a non-negligible effect on the GPS orbits, mainly in the radial component. In this paper, the effect on the along-track and cross-track components is studied in more detail. Also in this paper, it is shown that Earth radiation pressure leads to a change in the estimates of GPS ground station positions, which is systematic over large regions of the Earth. This observed “deformation” of the Earth is towards North–South and with large scale patterns that repeat six times per GPS draconitic year (350 days), reaching a magnitude of up to 1 mm. The impact of Earth radiation pressure on the geocenter and length of day estimates was also investigated, but the effect is found to be less significant as compared to the orbits and position estimates.     

An assessment of Bernese GPS software precise point positioning using IGS final products for global site velocities

We assess the use of precise point positioning (PPP) within the Bernese GPS software (BSW) Version 5.0 over the period from 2000 to 2005. In our strategy, we compute a set of daily PPP solutions for international GNSS service (IGS) reference frame (IGb00) sites by fixing IGS final satellite orbits and clock products, followed by a Helmert transformation of these solutions into ITRF2000, forming a set of continuous position time series over the entire time span. We assess BSW PPP by comparing our set of transformation parameters with those produced by the IGS analysis centre coordinator (ACC) and our position time series with those of the Jet Propulsion Laboratory (JPL) and the Scripps Orbit and Permanent Array Centre at the Scripps Institute of Oceanography (SIO). The distributions of the north (N), east (E) and up (U) daily position differences are characterized by means and SD of +2.2 ± 4.8, −0.6 ± 7.9 and +4.8 ± 17.3 mm with respect to JPL, and of +0.1 ± 4.4, −0.1 ± 7.4 and −0.1 ± 11.8 mm with respect to SIO. Similarly, we find sub-millimetre mean velocity differences and SD for the N, E and U components of 0.9, 1.5 and 2.2 mm/year with JPL, and of 1.2, 1.6 and 2.3 mm/year with SIO. A noise analysis using maximum likelihood estimation (MLE) shows that when estimating global site velocities from our position time series, the series need to be on average up to 1.3 times longer than those of JPL and SIO, before an uncertainty of less than 0.5 mm/year is obtained.     

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.     

Rigorous transformation of variance–covariance matrices of GPS-derived coordinates and velocities

With the advances in the field of GPS positioning and the global densification of permanent GPS tracking stations, it is now possible to determine at the highest level of accuracy the transformation parameters connecting various international terrestrial reference frame (ITRF) realizations. As a by-product of these refinements, not only the seven usual parameters of the similarity transformations between frames are available, but also their rates, all given at some epoch t _k . This paper introduces rigorous matrix equations to estimate variance–covariance matrices for transformed coordinates at any epoch t based on a stochastic model that takes into consideration all a priori information of the parameters involved at epoch t _k , and the coordinates and velocities at the reference frame initial epoch t ₀. The results of this investigation suggest that in order to attain maximum accuracy, the agencies determining the 14-parameter transformations between reference frames should also publish their full variance–covariance matrix. Electronic Publication     

Geodetic measurement of horizontal strain across the Red River fault near Thac Ba, Vietnam, 1963–1994

Geodetic measurements from 1963 through 1994 are used to estimate horizontal strain rates across the Red River fault near Thac Ba, Vietnam. Whether or not this fault system is currently active is a subject of some debate. By combining: (1) triangulation from 1963, (2) triangulation in 1983, and (3) Global Positioning System (GPS) observations in 1994, horizontal shear strain rates are estimated without imposing any prior information on fixed stations. The estimated rates of shear strain in ten triangular subnetworks surrounding the fault trace are not significantly different from zero at 95% confidence. The maximum rate of dextral shear is less than 0.3 μrad/year in all but one of the triangles. The estimates help bound the slip rate in a simple elastic dislocation model for a locked, vertical strike-slip fault. By assuming a locking depth of 5–20 km, the most likely values for the deep slip rate are between 1 and 5 mm/year of right-lateral motion. These values delimit the 23% confidence interval. At 95% confidence, the slip rate estimate falls between 7 mm/year of left-lateral motion and 15 mm/year of right-lateral motion. Received: 18 November 1997 / Accepted: 28 January 1999     

New results in airborne vector gravimetry using strapdown INS/DGPS

A method for airborne vector gravimetry has been developed. The method is based on developing the error dynamics equations of the INS in the inertial frame where the INS system errors are estimated in a wave estimator using inertial GPS position as update. Then using the error-corrected INS acceleration and the GPS acceleration in the inertial frame, the gravity disturbance vector is extracted. In the paper, the focus is on the improvement of accuracy for the horizontal components of the airborne gravity vector. This is achieved by using a decoupled model in the wave estimator and decorrelating the gravity disturbance from the INS system errors through the estimation process. The results of this method on the real strapdown INS/DGPS data are promising. The internal accuracy of the horizontal components of the estimated gravity disturbance for repeated airborne lines is comparable with the accuracy of the down component and is about 4–8 mGal. Better accuracy (2–4 mGal) is achieved after applying a wave-number correlation filter (WCF) to the parallel lines of the estimated airborne gravity disturbances.     

Monument-antenna effects on GPS coordinate time series with application to vertical rates in Antarctica

We examine the electromagnetic coupling of a GPS antenna–monument pair in terms of its simulated affect on long GPS coordinate time series. We focus on the Earth and Polar Observing System (POLENET) monument design widely deployed in Antarctica and Greenland in projects interested particularly in vertical velocities. We base our tests on an absolute robot calibration that included the top ~0.15 m of the monument and use simulations to assess its effect on site coordinate time series at eight representative POLENET sites in Antarctica over the period 2000.0–2011.0. We show that the neglect of this calibration would introduce mean coordinate bias, and most importantly for velocity estimation, coordinate noise which is highly sensitive to observation geometry and hence site location and observation period. Considering only sub-periods longer than 2.5 years, we show vertical site velocities may be biased by up to ±0.4 mm/year, and biases up to 0.2 mm/year may persist for observation spans of 8 years. Changing between uniform and elevation-dependent observation weighting alters the time series but does not remove the velocity biases, nor does ambiguity fixing. The effect on the horizontal coordinates is negligible. The ambiguities fixed series spectra show noise between flicker and random walk with near-white noise at the highest frequencies, with mean spectral indices (frequencies <20 cycles per year) of approximately −1.3 (uniform weighting) and −1.4 (elevation-dependent weighting). While the results are likely highly monument specific, they highlight the importance of accounting for monument effects when analysing vertical coordinate time series and velocities for the highest precision and accuracy geophysical studies.     

Modeling and estimation of C1–P1 bias in GPS receivers

 Modern dual-frequency global positioning system (GPS) receivers are capable of providing direct measurements of both L1 C/A (C1) and P code (P1) without the use of the Y-codes under Anti-Spoofing. A discrepancy or bias between the C1 and P1 measurements from these receivers has however been of concern to operators and users of GPS reference networks. For the purpose of modeling and estimation, the nature and characteristics of the discrepancy must be investigated. The research results presented indicate that the discrepancy between the C1 and P1 measurements contains two different types of components: one is of constant type while another is time variant. A method has been developed for their modeling and estimation. The residual C1–P1 time series after a satellite-dependent bias removal agree at a few-centimeter level, indicating the effectiveness of the proposed model. This allows the C1–P1 discrepancy, both constant and non-constant components, to be removed from GPS reference network solutions. Numerical results are provided to support the analysis. Received: 12 October 1999 / Accepted: 11 May 2000     

Adaptive Kalman Filtering for INS/GPS

After reviewing the two main approaches of adaptive Kalman filtering, namely, innovation-based adaptive estimation (IAE) and multiple-model-based adaptive estimation (MMAE), the detailed development of an innovation-based adaptive Kalman filter for an integrated inertial navigation system/global positioning system (INS/GPS) is given. The developed adaptive Kalman filter is based on the maximum likelihood criterion for the proper choice of the filter weight and hence the filter gain factors. Results from two kinematic field tests in which the INS/GPS was compared to highly precise reference data are presented. Results show that the adaptive Kalman filter outperforms the conventional Kalman filter by tuning either the system noise variance–covariance (V–C) matrix `Q' or the update measurement noise V–C matrix `R' or both of them. Received: 14 September 1998 / Accepted: 21 December 1998     

Absolute seafloor vertical positioning using combined pressure gauge and kinematic GPS data

Knowledge of the position and motion of points on the seafloor can be critically important in both fundamental research (for example, global geodesy and plate tectonics) and for more practical applications such as seismic risk evaluation, off-shore construction and pipeline monitoring. In the Vanuatu subduction zone, for example, measuring deformation underwater could provide valuable information for modeling deformation and understanding the seismic cycle. We report a shallow water experiment in Vanuatu to measure the relative and absolute depth of seafloor points. The experiment differs from previous efforts mainly in that it uses the height of the sea surface determined by kinematic GPS, allowing us to locate the points in a global reference frame. The ITRF2005 ellipsoidal height of a seafloor benchmark was determined with a 1-sigma uncertainty of 0.7–2.1 cm. The estimated ellipsoidal height differs only by a few tenths of a centimeter between measurements made in 2004 and another set made in 2006. These results are encouraging and open new perspectives for vertical underwater deformation monitoring in shallow water areas. Sea-surface GPS measurements can also help to reduce the uncertainty in depth difference determination for relative measurements.     

Absolute phase center corrections of satellite and receiver antennas

Results of the estimation of azimuth-dependent phase center variations (PCVs) of GPS satellite antennas using global GPS data are presented. Significant variations of up to ±3–4 mm that are demonstrated show excellent repeatability over eight years. The application of the azimuthal PCVs besides the nadir-dependent ones will lead to a further reduction in systematic antenna effects. In addition, the paper focuses on the benefit of a possible transition from relative to absolute PCVs. Apart from systematic changes in the global station coordinates, one can expect the GPS results to be less dependent on the elevation cut-off angle. This, together with the significant reduction of tropospheric zenith delay biases between GPS and VLBI, stands for an important step toward more consistency between different space geodetic techniques.