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     

GPS vector configuration design for monitoring deformation networks

 The performance of geodetic monitoring networks is heavily influenced by the configuration of the measured GPS vectors. As an effective design of the GPS measurements will decrease GPS campaign costs and increase the accuracy and reliability of the entire network, the identification of the preferred GPS vectors for measurement has been highlighted as a core problem in the process of deformation monitoring. An algorithm based on a sensitivity analysis of the network, as dependent upon a postulated velocity field, is suggested for the selection of the optimal GPS vectors. Relevant mathematical and statistical concepts are presented as the basis for an improved method of vector configuration design. A sensitivity analysis of the geodetic geodynamic network in the north of Israel is presented, where the method is examined against two deformation models, the Simple Transform Fault and the Locked Fault. The proposed method is suggested as a means for the improvement of the design of monitoring networks, a common practice worldwide. Received: 30 July 2001 / Accepted: 3 June 2002 Acknowledgments. It is my pleasant duty to thank the Survey of Israel and Dr. E. Ostrovsky for providing the variance–covariance matrix of the G1 network in northern Israel. I would like to thank the reviewers of this paper for their constructive and helpful remarks.     

An improvement to ionospheric delay correction for single-frequency GPS users – the APR-I scheme

 For the commonly used GPS wide-area augmentation systems (WAAS) with a grid ionospheric model, the efficient modelling of ionospheric delays in real time, for single-frequency GPS users, is still a crucial issue which needs further research. This is particularly necessary when differential ionospheric delay corrections cannot be broadcast, when users cannot receive them, or when there are ionospheric anomalies. Ionospheric delays have a severe effect on navigation performance of single-frequency receivers. A new scheme is proposed which can efficiently address the above problems. The robust recurrence technique is based on the efficient combination of single-frequency GPS observations by users and the high-precision differential ionospheric delay corrections from WAAS. Its effectiveness is verified with examples. Received: 24 December 1999 / Accepted 21 February 2001     

Spherical harmonic analysis and synthesis for global multiresolution applications

 On the Earth and in its neighborhood, spherical harmonic analysis and synthesis are standard mathematical procedures for scalar, vector and tensor fields. However, with the advent of multiresolution applications, additional considerations about convolution filtering with decimation and dilation are required. As global applications often imply discrete observations on regular grids, computational challenges arise and conflicting claims about spherical harmonic transforms have recently appeared in the literature. Following an overview of general multiresolution analysis and synthesis, spherical harmonic transforms are discussed for discrete global computations. For the necessary multi-rate filtering operations, spherical convolutions along with decimations and dilations are discussed, with practical examples of applications. Concluding remarks are then included for general applications, with some discussion of the computational complexity involved and the ongoing investigations in research centers. Received: 13 November 2000 / Accepted: 12 June 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     

Auto-covariance estimation of variable samples (ACEVS) and its application for monitoring random ionospheric disturbances using GPS

 Ionospheric variation may be considered as a stationary time series under quiet conditions. However, the disturbance of a stationary random process from stationarity results in the bias of corresponding samples from the stationary observations, and in the change of statistical model parameters of the process. From a general mathematical aspect, a new method is presented for monitoring ionospheric variations, based on the characteristic of time-series observation of GPS, and an investigation of the statistical properties of the estimated auto-covariance of the random ionospheric delay when changing the number of samples in the time series is carried out. A preliminary scheme for monitoring ionospheric delays is proposed. Received: 18 August 2000 / Accepted: 12 April 2001     

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     

Simplified formulae for the BIQUE estimation of variance components in disjunctive observation groups

 General rigorous and simplified formulae are reported for the best invariant quadratic unbiased estimates of the variance–covariance components, which can be applied to all least-squares adjustments with the general linear stochastic model. Simplified procedures are given for two cases frequently recurring in geodetic applications: uncorrelated groups of correlated or uncorrelated observations, with more than one variance component in each group. Received: 19 November 1998 / Accepted: 21 March 2000     

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     

Improved convergence rates for the truncation error in gravimetric geoid determination

 When Stokes's integral is used over a spherical cap to compute a gravimetric estimate of the geoid, a truncation error results due to the neglect of gravity data over the remainder of the Earth. Associated with the truncation error is an error kernel defined over these two complementary regions. An important observation is that the rate of decay of the coefficients of the series expansion for the truncation error in terms of Legendre polynomials is determined by the smoothness properties of the error kernel. Previously published deterministic modifications of Stokes's integration kernel involve either a discontinuity in the error kernel or its first derivative at the spherical cap radius. These kernels are generalised and extended by constructing error kernels whose derivatives at the spherical cap radius are continuous up to an arbitrary order. This construction is achieved by smoothly continuing the error kernel function into the spherical cap using a suitable degree polynomial. Accordingly, an improved rate of convergence of the spectral series representation of the truncation error is obtained. Received: 21 April 1998 / Accepted: 4 October 1999     

Precise estimation of residual tropospheric delays using a regional GPS network for real-time kinematic applications

 A new method called Trop_NetAdjust is described to predict in real time the residual tropospheric delays on the GPS carrier phase observables using the redundant measurements from a network of GPS reference stations. This method can not only enhance the effectiveness and reliability of real-time kinematic users within the network, but also provide a valid approach to tropospheric parameter variation forecasting. Trop_NetAdjust is theoretically based upon LS prediction criteria and enables the prediction of residual tropospheric delays remaining after a standard model has been applied to the raw GPS measurements. Two cases are analyzed, namely a first case when the delay is required for an existing satellite at a new point within the network and a second case when the delay is required for a new satellite. Field tests were conducted using data collected in a network of 11 reference stations covering a 400×600 km region in southern Norway. The results were analyzed in the measurement domain (ionospheric-free double-difference residuals) and showed improvements of 20 to 65% RMS errors using Trop_NetAdjust. The estimates of the Trop_NetAdjust prediction accuracy were also obtained using the covariance analysis method. The agreement was consistently better than 30% when compared with data from a real network. Received: 28 February 2000 / Accepted: 9 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     

A new method for computing the ellipsoidal correction for Stokes's formula

 This paper generalizes the Stokes formula from the spherical boundary surface to the ellipsoidal boundary surface. The resulting solution (ellipsoidal geoidal height), consisting of two parts, i.e. the spherical geoidal height N ₀ evaluated from Stokes's formula and the ellipsoidal correction N ₁, makes the relative geoidal height error decrease from O(e ²) to O(e ⁴), which can be neglected for most practical purposes. The ellipsoidal correction N ₁ is expressed as a sum of an integral about the spherical geoidal height N ₀ and a simple analytical function of N ₀ and the first three geopotential coefficients. The kernel function in the integral has the same degree of singularity at the origin as the original Stokes function. A brief comparison among this and other solutions shows that this solution is more effective than the solutions of Molodensky et al. and Moritz and, when the evaluation of the ellipsoidal correction N ₁ is done in an area where the spherical geoidal height N ₀ has already been evaluated, it is also more effective than the solution of Martinec and Grafarend. Received: 27 January 1999 / Accepted: 4 October 1999     

Short-arc analysis of intersatellite tracking data in a gravity mapping mission

 A technique for the analysis of low–low intersatellite range-rate data in a gravity mapping mission is explored. The technique is based on standard tracking data analysis for orbit determination but uses a spherical coordinate representation of the 12 epoch state parameters describing the baseline between the two satellites. This representation of the state parameters is exploited to allow the intersatellite range-rate analysis to benefit from information provided by other tracking data types without large simultaneous multiple-data-type solutions. The technique appears especially valuable for estimating gravity from short arcs (e.g. less than 15 minutes) of data. Gravity recovery simulations which use short arcs are compared with those using arcs a day in length. For a high-inclination orbit, the short-arc analysis recovers low-order gravity coefficients remarkably well, although higher-order terms, especially sectorial terms, are less accurate. Simulations suggest that either long or short arcs of the Gravity Recovery and Climate Experiment (GRACE) data are likely to improve parts of the geopotential spectrum by orders of magnitude. Received: 26 June 2001 / Accepted: 21 January 2002     

New solutions for the geodetic coordinate transformation

 The Cartesian-to-geodetic-coordinate transformation is approached from a new perspective. Existence and uniqueness of geodetic representation are presented, along with a clear geometric picture of the problem and the role of the ellipse evolute. A new solution is found with a Newton-method iteration in the reduced latitude; this solution is proved to work for all points in space. Care is given to error propagation when calculating the geodetic latitude and height. Received: 9 August 2001 / Accepted: 27 March 2002 Acknowledgments. The author would like to thank the Clifford W.␣Tompson scholarship fund, Dr. Brian DeFacio, the University of Missouri College of Arts &Sciences, and the United States Air Force. He also thanks a reviewer for suggesting and providing a prototype MATLAB code. A MATLAB program for the iterative sequence is presented at the end of the paper (Appendix A).     

Molodensky potential telluroid based on a minimum-distance map. Case study: the quasi-geoid of East Germany in the World Geodetic Datum 2000

 A potential-type Molodensky telluroid based upon a minimum-distance mapping is derived. With respect to a reference potential of Somigliana–Pizzetti type which relates to the World Geodetic Datum 2000, it is shown that a point-wise minimum-distance mapping of the topographical surface of the Earth onto the telluroid surface, constrained to the gauge W(P)=u(p), leads to a system of four nonlinear normal equations. These normal equations are solved by a fast Newton–Raphson iteration. Received: 7 February 2000 / Accepted: 23 October 2001     

A bound for the Euclidean norm of the difference between the best linear unbiased estimator and a linear unbiased estimator

 A bound is established for the Euclidean norm of the difference between the best linear unbiased estimator and any linear unbiased estimator in the general linear model. The bound involves the spectral norm of the difference between the dispersion matrices of the two estimators, and the residual sum of squares, all evaluated at the assumed model, but is independent of the provenance of the observation vector at hand. The bound, a straightforward consequence of first principles in Gauss–Markov theory, generalizes previous results on the difference between the best linear unbiased estimator and the ordinary least-squares estimator. In a numerical example from repeated precise levelling, the bound is used to analyse the sensitivity of estimates of vertical motion to the choice of estimator. Received: 9 September 1999 / Accepted: 15 March 2002     

A synthetic Earth for use in geodesy

 A synthetic Earth and its gravity field that can be represented at different resolutions for testing and comparing existing and new methods used for global gravity-field determination are created. Both the boundary and boundary values of the gravity potential can be generated. The approach chosen also allows observables to be generated at aircraft flight height or at satellite altitude. The generation of the synthetic Earth shape (SES) and gravity-field quantities is based upon spherical harmonic expansions of the isostatically compensated equivalent rock topography and the EGM96 global geopotential model. Spherical harmonic models are developed for both the synthetic Earth topography (SET) and the synthetic Earth potential (SEP) up to degree and order 2160 corresponding to a 5′×5′ resolution. Various sets of SET, SES and SEP with boundary geometry and boundary values at different resolutions can be generated using low-pass filters applied to the expansions. The representation is achieved in point sets based upon refined triangulation of a octahedral geometry projected onto the chosen reference ellipsoid. The filter cut-offs relate to the sampling pattern in order to avoid aliasing effects. Examples of the SET and its gravity field are shown for a resolution with a Nyquist sampling rate of 8.27 degrees. Received: 6 August 1999 / Accepted: 26 April 2000     

Integer estimation in the presence of biases

 Carrier phase ambiguity resolution is the key to fast and high-precision GNSS (Global Navigation Satellite System) kinematic positioning. Critical in the application of ambiguity resolution is the quality of the computed integer ambiguities. Unsuccessful ambiguity resolution, when passed unnoticed, will too often lead to unacceptable errors in the positioning results. Very high success rates are therefore required for ambiguity resolution to be reliable. Biases which are unaccounted for will lower the success rate and thus increase the chance of unsuccessful ambiguity resolution. The performance of integer ambiguity estimation in the presence of such biases is studied. Particular attention is given to integer rounding, integer bootstrapping and integer least squares. Lower and upper bounds, as well as an exact and easy-to-compute formula for the bias-affected success rate, are presented. These results will enable the evaluation of the bias robustness of ambiguity resolution. Received: 28 September 2000 / Accepted: 29 March 2001     

Combination procedure for length-of-day time series according to the noise frequency behavior

 Length-of-day (LOD) estimates from seven Global Positioning System (GPS) and three satellite laser ranging (SLR) analysis centers were combined into an even-spaced time series for a 27-month period (1996–1998). This time series was compared to the multi-technique Earth-orientation-parameter (EOP) combined solution (C04) derived at the Central Bureau of the International Earth Rotation Service (IERS/CB). Due to inhomogeneities in the different series derived from the various techniques (time, length, quality, and spatial resolution), the concept of a combined solution is justified. The noise behavior in LOD for different techniques varies with frequency; the data series were divided into frequency windows after removing both biases and trends. Different weight factors were assigned in each window, discriminating by technique, and produced one-technique combined solutions. Finally, these one-technique combined solutions were combined to obtain the final multi-technique solution. The LOD combined time series obtained by the present method based on the frequency windows combined series (FWCS) is very close to the IERS C04 solution. It could be useful to generate a new LOD reference time series to be used in the study of high-frequency variations of Earth rotation. Received: 28 March 2000 / Accepted: 15 February 2001