GPS-assisted GLONASS orbit determination

 Using 1 week of data from a network of GPS/GLONASS dual-tracking receivers, 15-cm accurate GLONASS orbit determination is demonstrated with an approach that combines GPS and GLONASS data. GPS data are used to define the reference frame, synchronize receiver clocks and determine troposphere delay for the GLONASS tracking network. GLONASS tracking data are then processed separately, with the GPS-defined parameters held fixed, to determine the GLONASS orbit. The quality of the GLONASS orbit determination is currently limited by the size and distribution of the tracking network, and by the unavailability of a sufficiently refined solar pressure model. Temporal variations in the differential clock bias of the dual-tracking receivers are found to have secondary impact on the orbit determination accuracy. Received: 5 January 2000 / Accepted: 15 February 2001     

Analytical solar radiation pressure modelling for GLONASS using a pixel array

 A new method for calculating analytical solar radiation pressure models for GNSS spacecraft has been developed. The method simulates the flux of light from the Sun using a pixel array. The method can cope with a high level of complexity in the spacecraft structure and models effects due to reflected light. Models have been calculated and tested for the Russhar global navigation satellite system GLONASS IIv spacecraft. Results are presented using numerical integration of the force model and long-arc satellite laser ranging (SLR) analysis. The integrated trajectory differs from a precise orbit calculated using a network of global tracking stations by circa 2 m root mean square over a 160 000-km arc. The observed − computed residuals for the 400-day SLR arc are circa 28 mm. Received: 23 December 1999 / Accepted: 28 August 2000     

Laser-based validation of GLONASS orbits by short-arc technique

 The International GLONASS Experiment (IGEX-98) was carried out between 19 October 1998 and 19 April 1999. Among several objectives was the precise orbit determination of GPS and GLONASS satellites and its validation by laser ranging observations. Local laser-based orbit corrections (radial, tangential and normal components in a rotating orbital local reference frame) are computed using a geometrical short-arc technique. The order of magnitude of these corrections is at the level of few decimeters, depending on the considered components. The orbit corrections are analyzed as a function of several parameters (date, orbital plane, geographical area). The mean corrections are at the level of several centimeters. However, when averaging over the entire campaign and for all the satellites, no mean radial, tangential and normal orbit corrections are found. The origin of the observed corrections is considered (errors due to the geocentric gravitational constant, the non-gravitational forces, the thermal equilibrium of on-board equipment, the reference systems, the location and the signature of the retroreflector array, and the precision of the satellite laser ranges). Some features are also due to errors in the radio-tracking GLONASS orbits. Further investigations will be needed to better understand the origin of various biases. Received: 17 February 2000 / Accepted: 31 January 2001     

Combined processing of the IGS and the IGEX network

 Since the beginning of the International Global Navigation Satellite System (GLONASS) Experiment, IGEX, in October 1998, the Center for Orbit Determination in Europe (CODE) has acted as an analysis center providing precise GLONASS orbits on a regular basis. In CODE's IGEX routine analysis the Global Positioning System (GPS) orbits and Earth rotation parameters are introduced as known quantities into the GLONASS processing. A new approach is studied, where data from the IGEX network are combined with GPS observations from the International GPS Service (IGS) network and all parameters (GPS and GLONASS orbits, Earth rotation parameters, and site coordinates) are estimated in one processing step. The influence of different solar radiation pressure parameterizations on the GLONASS orbits is studied using different parameter subsets of the extended CODE orbit model. Parameterization with three constant terms in the three orthogonal directions, D, Y, and X (D = direction satellite–Sun, Y = direction of the satellite's solar panel axis), and two periodic terms in the X-direction, proves to be adequate for GLONASS satellites. As a result of the processing it is found that the solar radiation pressure effect for the GLONASS satellites is significantly different in the Y-direction from that for the GPS satellites, and an extensive analysis is carried out to investigate the effect in detail. SLR observations from the ILRS network are used as an independent check on the quality of the GLONASS orbital solutions. Both processing aspects, combining the two networks and changing the orbit parameterization, significantly improve the quality of the determined GLONASS orbits compared to the orbits stemming from CODE's IGEX routine processing. Received: 10 May 2000 / Accepted: 9 October 2000     

An approach to GLONASS ambiguity resolution

 When processing global navigation satellite system (GLONASS) carrier phases, the standard double-differencing (DD) procedure cannot cancel receiver clock terms in the DD phase measurement equations due to the multiple frequencies of the carrier phases. Consequently, a receiver clock parameter has to be set up in the measurement equations in addition to baseline components and DD ambiguities. The resulting normal matrix unfortunately becomes singular. Methods to deal with this problem have been proposed in the literature. However, these methods rely on the use of pseudo-ranges. As pseudo-ranges are contaminated by multi-path and hardware delays, biases in these pseudo-ranges are significant, which may result in unreliable ambiguity resolution. A new approach is addressed that is not sensitive to the biases in the pseudo-ranges. The proposed approach includes such steps as converting the carrier phases to their distances to cancel the receiver clock errors, and searching for the most likely single-differenced (SD) ambiguity. Based on the results from the theoretical investigation, a practical procedure for GLONASS ambiguity resolution is presented. The initial experimental results demonstrate that the proposed approach is useable in cases of GLONASS and combined global positioning system (GPS) and GLONASS positioning. Received: 19 August 1998 / Accepted: 12 November 1999     

The impact of accelerometry on CHAMP orbit determination

 The contribution of the STAR accelerometer to the CHAMP orbit precision is evaluated and quantified by means of the following results: orbital fit to the satellite laser ranging (SLR) observations, GPS reduced-dynamic vs SLR dynamic orbit comparisons, and comparison of the measured to the modeled non-gravitational accelerations (atmospheric drag in particular). In each of the four test periods in 2001, five CHAMP arcs of 2 days' length were analyzed. The mean RMS-of-fit of the SLR observations of the orbits computed with STAR data or the non-gravitational force model were 11 and 24 cm, respectively. If the accelerometer calibration parameters are not known at least at the few percent level, the SLR orbit fit deteriorates. This was tested by applying a 10% error to the along-track scale factor of the accelerometer, which increased the SLR RMS-of-fit on average to 17 cm. Reference orbits were computed employing the reduced-dynamic technique with GPS tracking data. This technique yields the most accurate orbit positions thanks to the estimation of a large number of empirical accelerations, which compensate for dynamic modeling errors. Comparison of the SLR orbits, computed with STAR data or the non-gravitational force model, to the GPS-based orbits showed that the SLR orbits employing accelerometer observations are twice as accurate. Finally, comparison of measured to modeled accelerations showed that the level of geomagnetic activity is highly correlated with the atmospheric drag model error, and that the largest errors occur around the geomagnetic poles. Received: 7 May 2002 / Accepted: 18 November 2002 Correspondence to: S. Bruinsma Acknowledgments. The TIGCM results were obtained from the CEDAR database. This study was supported by the Centre National d'Etudes Spatiales (CNES). The referees are thanked for their helpful remarks and suggestions.     

First use of IGEX precise ephemerides for intercontinental GLONASS P-code time transfer

 Until recently, the Global Positioning System (GPS) was the only operational means of distributing time to an arbitrary number of users and of synchronizing clocks over large distances with a high degree of precision and accuracy. Over the last few years it has been shown that similar performance can be achieved using the Russian Global Navigation Satellite System (GLONASS). GLONASS time transfer between continents was initially hampered by the lack of post-processed precise ephemerides. Results from the International GLONASS Experiment (IGEX) campaign are now available, however, and this paper reports on the first use of IGEX precise ephemerides for GLONASS P-code intercontinental time links. The results of GLONASS P-code and GPS C/A-code time transfer are compared under similar conditions. Received: 31 January 2000 / Accepted: 10 July 2000     

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     

Somigliana–Pizzetti gravity: the international gravity formula accurate to the sub-nanoGal level

 The Somigliana–Pizzetti gravity field (the International gravity formula), namely the gravity field of the level ellipsoid (the International Reference Ellipsoid), is derived to the sub-nanoGal accuracy level in order to fulfil the demands of modern gravimetry (absolute gravimeters, super conducting gravimeters, atomic gravimeters). Equations (53), (54) and (59) summarise Somigliana–Pizzetti gravity Γ(φ,u) as a function of Jacobi spheroidal latitude φ and height u to the order ?(10⁻¹⁰ Gal), and Γ(B,H) as a function of Gauss (surface normal) ellipsoidal latitude B and height H to the order ?(10⁻¹⁰ Gal) as determined by GPS (`global problem solver'). Within the test area of the state of Baden-Württemberg, Somigliana–Pizzetti gravity disturbances of an average of 25.452 mGal were produced. Computer programs for an operational application of the new international gravity formula with (L,B,H) or (λ,φ,u) coordinate inputs to a sub-nanoGal level of accuracy are available on the Internet. Received: 23 June 2000 / Accepted: 2 January 2001     

Prediction of surface horizontal displacements, and gravity and tilt changes caused by filling the Three Gorges Reservoir

 Horizontal displacements, and gravity and tilt changes induced by filling the Three Gorges Reservoir are modeled using elastic loading Green functions. When the water surface reaches its highest level, the effects become maximum on the reservoir banks. The longitudinal and latitudinal components of the horizontal displacements reach −8.2 and 7.7 mm respectively, gravity is increased by up to 3.4 mGal, and the prime vertical and meridian components of the tilt changes are −7.8 and −17.5 arcseconds respectively. Accordingly, the filling of the reservoir will influence values observed from global positioning system (GPS), gravimetry and tilt measurements in the area. The results given can be used to provide important corrections for extracting earthquake-related signals from observed data. Received: 19 January 2001 / Accepted: 3 September 2001     

Local geoid determination combining gravity disturbances and GPS/levelling: a case study in the Lake Nasser area, Aswan, Egypt

 The use of GPS for height control in an area with existing levelling data requires the determination of a local geoid and the bias between the local levelling datum and the one implicitly defined when computing the local geoid. If only scarse gravity data are available, the heights of new data may be collected rapidly by determining the ellipsoidal height by GPS and not using orthometric heights. Hence the geoid determination has to be based on gravity disturbances contingently combined with gravity anomalies. Furthermore, existing GPS/levelling data may also be used in the geoid determination if a suitable general gravity field modelling method (such as least-squares collocation, LSC) is applied. A comparison has been made in the Aswan Dam area between geoids determined using fast Fourier transform (FFT) with gravity disturbances exclusively and LSC using only the gravity disturbances and the disturbances combined with GPS/levelling data. The EGM96 spherical harmonic model was in all cases used in a remove–restore mode. A total of 198 gravity disturbances spaced approximately 3 km apart were used, as well as 35 GPS/levelling points in the vicinity and on the Aswan Dam. No data on the Nasser Lake were available. This gave difficulties when using FFT, which requires the use of gridded data. When using exclusively the gravity disturbances, the agreement between the GPS/levelling data were 0.71 ± 0.17 m for FFT and 0.63 ± 0.15 for LSC. When combining gravity disturbances and GPS/levelling, the LSC error estimate was ±0.10 m. In the latter case two bias parameters had to be introduced to account for a possible levelling datum difference between the levelling on the dam and that on the adjacent roads. Received: 14 August 2000 / Accepted: 28 February 2001     

The AUSGeoid98 geoid model of Australia: data treatment, computations and comparisons with GPS-levelling data

 The AUSGeoid98 gravimetric geoid model of Australia has been computed using data from the EGM96 global geopotential model, the 1996 release of the Australian gravity database, a nationwide digital elevation model, and satellite altimeter-derived marine gravity anomalies. The geoid heights are on a 2 by 2 arc-minute grid with respect to the GRS80 ellipsoid, and residual geoid heights were computed using the 1-D fast Fourier transform technique. This has been adapted to include a deterministically modified kernel over a spherical cap of limited spatial extent in the generalised Stokes scheme. Comparisons of AUSGeoid98 with GPS and Australian Height Datum (AHD) heights across the continent give an RMS agreement of ±0.364 m, although this apparently large value is attributed partly to distortions in the AHD. Received: 10 March 2000 / Accepted: 21 February 2001     

The ERS-2 altimetric bias and gravity field enhancement using dual crossovers between ERS and TOPEX/Poseidon

 Dual satellite crossovers (DXO) between the two European Remote Sensing satellites ERS-1 and ERS-2 and TOPEX/Poseidon are used to (1) refine the Earth's gravity field and (2) extend the study of the ERS-2 altimetric range stability to cover the first four years of its operation. The enhanced gravity field model, AGM-98, is validated by several methodologies and will be shown to provide, in particular, low geographically correlated orbital error for ERS-2. For the ERS-2 altimetric range study, TOPEX/Poseidon is first calibrated through comparison against in situ tide gauge data. A time series of the ERS-2 altimeter bias has been recovered along with other geophysical correction terms using tables for bias jumps in the range measurements at the single point target response (SPTR) events. On utilising the original version of the SPTR tables the overall bias drift is seen to be 2.6±1.0 mm/yr with an RMS of fit of 12.2 mm but with discontinuities at the centimetre level at the SPTR events. On utilising the recently released revised tables, SPTR2000, the drift is better defined at 2.4±0.6 mm/yr with the RMS of fit reduced to 3.7 mm. Investigations identify the sea-state bias as a source of error with corrections affecting the overall drift by close to 1.2 mm/yr. Received: 25 May 2000 / Accepted: 24 January 2001     

The gravitational potential and its derivatives for the prism

 As a simple building block, the right rectangular parallelepiped (prism) has an important role mostly in local gravity field modelling studies when the so called flat-Earth approximation is sufficient. Its primary (methodological) advantage follows from the simplicity of the rigorous and consistent analytical forms describing the different gravitation-related quantities. The analytical forms provide numerical values for these quantities which satisfy the functional connections existing between these quantities at the level of numerical precision applied. Closed expressions for the gravitational potential of the prism and its derivatives (up to the third order) are listed for easy reference. Received: 18 August 1999 / Accepted: 15 June 2000     

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     

The multiresolution character of collocation

 An interesting theoretical connection between the statistical (non-stochastic) collocation principle and the multiresolution/wavelet framework of signal approximation is presented. The rapid developments in multiresolution analysis theory over the past few years have provided very useful (theoretical and practical) tools for approximation and spectral studies of irregularly varying signals, thus opening new possibilities for `non-stationary' gravity field modeling. It is demonstrated that the classic multiresolution formalism according to Mallat's pioneering work lies at the very core of some of the general approximation principles traditionally used in physical geodesy problems. In particular, it is shown that the use of a spatio-statistical (non-probabilistic) minimum mean-square-error criterion for optimal linear estimation of deterministic signals, in conjunction with regularly gridded data, always gives rise to a generalized multiresolution analysis in the Hilbert space L <sup>2</sup>(R), under some mild constraints on the spatial covariance function and the power spectrum of the unknown field under consideration. Using the theory and the actual approximation algorithms associated with statistical collocation, a new constructive framework for building generalized multiresolution analyses in L <sup>2</sup>(R) is presented, without the need for the usual dyadic restriction that exists in classic wavelet theory. The multiresolution and `non-stationary' aspects of the statistical collocation approximation procedure are also discussed, and finally some conclusions and recommendations for future work are given. Received: 26 January 1999 / Accepted: 16 August 1999     

Fast spherical collocation: theory and examples

 It has long been known that a spherical harmonic analysis of gridded (and noisy) data on a sphere (with uniform error for a fixed latitude) gives rise to simple systems of equations. This idea has been generalized for the method of least-squares collocation, when using an isotropic covariance function or reproducing kernel. The data only need to be at the same altitude and of the same kind for each latitude. This permits, for example, the combination of gravity data at the surface of the Earth and data at satellite altitude, when the orbit is circular. Suppose that data are associated with the points of a grid with N values in latitude and M values in longitude. The latitudes do not need to be spaced uniformly. Also suppose that it is required to determine the spherical harmonic coefficients to a maximal degree and order K. Then the method will require that we solve K systems of equations each having a symmetric positive definite matrix of only N × N. Results of simulation studies using the method are described. Received: 18 October 2001 / Accepted: 4 October 2002 Correspondence to: F. Sansò     

The northern European geoid: a case study on long-wavelength geoid errors

 The long-wavelength geoid errors on large-scale geoid solutions, and the use of modified kernels to mitigate these effects, are studied. The geoid around the Nordic area, from Greenland to the Ural mountains, is considered. The effect of including additional gravity data around the Nordic/Baltic land area, originating from both marine, satellite and ground-based measurements, is studied. It is found that additional data appear to increase the noise level in computations, indicating the presence of systematic errors. Therefore, the Wong–Gore modification to the Stokes kernel is applied. This method of removing lower-order terms in the Stokes kernel appears to improve the geoid. The best fit to the global positioning system (GPS) leveling points is obtained with a degree of modification of approximately 30. In addition to the study of modification errors, the results of different methods of combining satellite altimetry gravity and other gravimetry are presented. They all gave comparable results, at the 6-cm level, when evaluated for the Nordic GPS networks. One dimensional (1-D) and 2-D fast Fourier transform (FFT) methods are also compared. It is shown that even though methods differ by up to 6 cm, the fit to the GPS is essentially the same. A surprising conclusion is that the addition of more data does not always produce a better geoid, illustrating the danger of systematic errors in data. Received: 4 July 2001 / Accepted: 21 February 2002     

The success rate and precision of GPS ambiguities

 An application of a theorem on the optimality of integer least-squares (LS) is described. This theorem states that the integer LS estimator maximizes the ambiguity success rate within the class of admissible integer estimators. This theorem is used to show how the probability of correct integer estimation depends on changes in the second moment of the ambiguity `float' solution. The distribution of the `float' solution is considered to be a member of the broad family of elliptically contoured distributions. Eigenvalue-based bounds for the ambiguity success rate are obtained. Received: 11 January 1999 / Accepted: 2 November 1999     

GPS measurements of ocean loading and its impact on zenith tropospheric delay estimates: a case study in Brittany, France

 The results from a global positioning system (GPS) experiment carried out in Brittany, France, in October 1999, aimed at measuring crustal displacements caused by ocean loading and quantifying their effects on GPS-derived tropospheric delay estimates, are presented. The loading effect in the vertical and horizontal position time series is identified, however with significant disagreement in amplitude compared to ocean loading model predictions. It is shown that these amplitude misfits result from spatial tropospheric heterogeneities not accounted for in the data processing. The effect of ocean loading on GPS-derived zenith total delay (ZTD) estimates is investigated and a scaling factor of 4.4 between ZTD and station height for a 10° elevation cut-off angle is found (i.e. a 4.4-cm station height error would map into a 1-cm ZTD error). Consequently, unmodeled ocean loading effects map into significant errors in ZTD estimates and ocean loading modeling must be properly implemented when estimating ZTD parameters from GPS data for meteorological applications. Ocean loading effects must be known with an accuracy of better than 3 cm in order to meet the accuracy requirements of meteorological and climatological applications of GPS-derived precipitable water vapor. Received: 16 July 2001 / Accepted: 25 April 2002 Acknowledgments. The authors are grateful to H.G. Scherneck for fruitful discussions and for his help with the ocean loading calculations. They thank H. Vedel for making the HIRLAM data available; D. Jerett for helpful discussions; and the city of Rostrenen, the Laboratoire d'Océanographie of Concarneau, and the Institut de Protection et de S?reté Nucléaire (BERSSIN) for their support during the GPS measurement campaign. Reviews by C.K. Shum and two anonymous referees significantly improved this paper. This work was carried out in the framework of the MAGIC project (http://www.acri.fr/magic), funded by the European Commission, Environment and Climate Program (EC Contract ENV4-CT98–0745). Correspondence to: E. Calais, Department of Earth and Atmospheric Sciences, Purdue University, West Lafayette, IN 47907-1397, USA. e-mail: ecalais@purdue.edu Tel. : +1-765-496-2915; Fax:+1-765-496-1210