Some modifications of Stokes' formula that account for truncation and potential coefficient errors

 Stokes' formula from 1849 is still the basis for the gravimetric determination of the geoid. The modification of the formula, originating with Molodensky, aims at reducing the truncation error outside a spherical cap of integration. This goal is still prevalent among various modifications. In contrast to these approaches, some least-squares types of modification that aim at reducing the truncation error, as well as the error stemming from the potential coefficients, are demonstrated. The least-squares estimators are provided in the two cases that (1) Stokes' kernel is a priori modified (e.g. according to Molodensky's approach) and (2) Stokes' kernel is optimally modified to minimize the global mean square error. Meissl-type modifications are also studied. In addition, the use of a higher than second-degree reference field versus the original (Pizzetti-type) reference field is discussed, and it is concluded that the former choice of reference field implies increased computer labour to achieve the same result as with the original reference field. Received: 14 December 1998 / Accepted: 4 October 1999     

A Meissl-modified Vaníček and Kleusberg kernel to reduce the truncation error in gravimetric geoid computations

A deterministic modification of Stokes's integration kernel is presented which reduces the truncation error when regional gravity data are used in conjunction with a global geopotential model to compute a gravimetric geoid. The modification makes use of a combination of two existing modifications from Vaníček and Kleusberg and Meissl. The former modification applies a root mean square minimisation to the upper bound of the truncation error, whilst the latter causes the Fourier series expansion of the truncation error to coverage to zero more rapidly by setting the kernel to zero at the truncation radius. Green's second identity is used to demonstrate that the truncation error converges to zero faster when a Meissl-type modification is made to the Vaníček and Kleusberg kernel. A special case of this modification is proposed by choosing the degree of modification and integration cap-size such that the Vaníček and Kleusberg kernel passes through zero at the truncation radius. Received: 14 October 1996 / Accepted: 20 October 1997     

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     

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     

A solution to the downward continuation effect on the geoid determined by Stokes' formula

 The analytical continuation of the surface gravity anomaly to sea level is a necessary correction in the application of Stokes' formula for geoid estimation. This process is frequently performed by the inversion of Poisson's integral formula for a sphere. Unfortunately, this integral equation corresponds to an improperly posed problem, and the solution is both numerically unstable, unless it is well smoothed, and tedious to compute. A solution that avoids the intermediate step of downward continuation of the gravity anomaly is presented. Instead the effect on the geoid as provided by Stokes' formula is studied directly. The practical solution is partly presented in terms of a truncated Taylor series and partly as a truncated series of spherical harmonics. Some simple numerical estimates show that the solution mostly meets the requests of a 1-cm geoid model, but the truncation error of the far zone must be studied more precisely for high altitudes of the computation point. In addition, it should be emphasized that the derived solution is more computer efficient than the detour by Poisson's integral. Received: 6 February 2002 / Accepted: 18 November 2002 Acknowledgements. Jonas ?gren carried out the numerical calculations and gave some critical and constructive remarks on a draft version of the paper. This support is cordially acknowledged. Also, the thorough work performed by one unknown reviewer is very much appreciated.     

Geoid determination with density hypotheses from isostatic models and geological information

 Geoid determination by Stokes's formula requires a complete knowledge of the topographical mass density distribution in order to perform gravity reductions to the geoid boundary. However, deeper masses are also of interest, in order to produce a smooth field of gravity anomalies which will improve results from interpolation procedures. Until now, in most cases a constant mass density has been considered, which is a very rough approximation of reality. The influence on the geoid height coming from different mass density hypotheses given by the isostatic models of Pratt/Hayford, Airy/Heiskanen and Vening Meinesz is studied. Apart from a constant mass density value, additional density information deduced from geological maps and thick sedimentary layers is considered. An overview of how mass density distributions act within Stokes's theory is given. The isostatic models are considered in spherical and planar approximation, as well as with constant and lateral variable mass density of the topographical and deeper masses. Numerical results in a test area in south-west Germany show that the differences in the geoid height due to different density hypotheses can reach a magnitude of more than 1 decimetre, which is not negligible in a precise geoid determination with centimetre accuracy. Received: 7 January 2002 / Accepted: 20 September 2002 M. Kuhn now at: Western Australian Centre for Geodesy, Curtin University of Technology, GPO Box U1987, Perth, WA 6845, Australia Acknowledgements. The author would gratefully thank Prof. Dr.-Ing. B. Heck, who was the supervisor of my PhD thesis, and the second examiner Prof. Dr.-Ing. K.H. Ilk, as well as all other colleagues for their support of this work. Particular thanks go to the Landesvermessungsamt Baden–Württemberg (Survey Department of Baden–Württemberg), Bureau Gravimetrique International (BGI, France) for providing the gravity data and the Geologisches Landesamt Baden–Württemberg (Geological Department of Baden–Württemberg) for providing data and maps of the sediment layers within the Rhine Valley. Grateful thanks goes to Prof. W.E. Featherstone and the reviewers Prof. S.D. Pagiatakis, Dr. U. Marti as well as an unknown reviewer for their helpful comments on this paper.     

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 temporal variation of the spherical and Cartesian multipoles of the gravity field: the generalized MacCullagh representation

 The Cartesian moments of the mass density of a gravitating body and the spherical harmonic coefficients of its gravitational field are related in a peculiar way. In particular, the products of inertia can be expressed by the spherical harmonic coefficients of the gravitational potential as was derived by MacCullagh for a rigid body. Here the MacCullagh formulae are extended to a deformable body which is restricted to radial symmetry in order to apply the Love–Shida hypothesis. The mass conservation law allows a representation of the incremental mass density by the respective excitation function. A representation of an arbitrary Cartesian monome is always possible by sums of solid spherical harmonics multiplied by powers of the radius. Introducing these representations into the definition of the Cartesian moments, an extension of the MacCullagh formulae is obtained. In particular, for excitation functions with a vanishing harmonic coefficient of degree zero, the (diagonal) incremental moments of inertia also can be represented by the excitation coefficients. Four types of excitation functions are considered, namely: (1) tidal excitation; (2) loading potential; (3) centrifugal potential; and (4) transverse surface stress. One application of the results could be model computation of the length-of-day variations and polar motion, which depend on the moments of inertia. Received: 27 July 1999 / Accepted: 24 May 2000     

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 spherical horizontal and spherical vertical boundary value problem – vertical deflections and geoidal undulations – the completed Meissl diagram

 In a comparison of the solution of the spherical horizontal and vertical boundary value problems of physical geodesy it is aimed to construct downward continuation operators for vertical deflections (surface gradient of the incremental gravitational potential) and for gravity disturbances (vertical derivative of the incremental gravitational potential) from points on the Earth's topographic surface or of the three-dimensional (3-D) Euclidean space nearby down to the international reference sphere (IRS). First the horizontal and vertical components of the gravity vector, namely spherical vertical deflections and spherical gravity disturbances, are set up. Second, the horizontal and vertical boundary value problem in spherical gravity and geometry space is considered. The incremental gravity vector is represented in terms of vector spherical harmonics. The solution of horizontal spherical boundary problem in terms of the horizontal vector-valued Green function converts vertical deflections given on the IRS to the incremental gravitational potential external in the 3-D Euclidean space. The horizontal Green functions specialized to evaluation and source points on the IRS coincide with the Stokes kernel for vertical deflections. Third, the vertical spherical boundary value problem is solved in terms of the vertical scalar-valued Green function. Fourth, the operators for upward continuation of vertical deflections given on the IRS to vertical deflections in its external 3-D Euclidean space are constructed. Fifth, the operators for upward continuation of incremental gravity given on the IRS to incremental gravity to the external 3-D Euclidean space are generated. Finally, Meissl-type diagrams for upward continuation and regularized downward continuation of horizontal and vertical gravity data, namely vertical deflection and incremental gravity, are produced. Received: 10 May 2000 / Accepted: 26 February 2001     

Meissel-Stokes核函数应用于区域大地水准面分析

为提高区域大地水准面计算精度,基于EGM2008地球重力场位系数模型分析Meissel-Stokes核函数、截断误差系数以及截断误差。选取实验区,采用移去-恢复法评价Meissel-Stokes核函数计算大地水准面的精度。结果表明:Meissel-Stokes核函数及其截断误差系数收敛速度快;截断误差小且稳定。在积分半径不易扩展的情况下,应用Meissel-Stokes核函数计算区域大地水准面,比标准Stokes计算大地水准面精度略高。     

Altimetry–gravimetry problems with free vertical datum

 The definition and connection of vertical datums in geodetic height networks is a fundamental problem in geodesy. Today, the standard approach to solve it is based on the joint processing of terrestrial and satellite geodetic data. It is generalized to cases where the coverage with terrestrial data may change from region to region, typically across coastlines. The principal difficulty is that such problems, so-called altimetry–gravimetry boundary-value problems (AGPs), do not admit analytical solutions such as Stokes' integral. A numerical solution strategy for the free-datum problem is presented. Analysis of AGPs in spherical and constant radius approximation shows that two of them are mathematically well-posed problems, while the classical AGP-I may be ill posed in special situations. Received: 2 December 1998 / Accepted: 30 November 1999     

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     

On long-period polar motion

 Five separate polar motion series are examined in order to understand what portion of their variations at periods exceeding several years represents true polar motion. The data since the development of space-geodetic techniques (by themselves insufficient for study of long-period motion), and a variety of historical astrometric data sets, allow the following tentative conclusions: retrograde long-period polar motion below about −0.2 cpy (cycles per year) in pre-space-geodetic data (pre-1976) is dominantly noise. For 1976–1992, there is poor agreement between space-geodetic and astrometric series over the range −0.2 to +0.2 cpy, demonstrating that classical astrometry lacked the precision to monitor polar motion in this frequency range. It is concluded that all the pre-1976 astrometric polar motion data are likely to be dominated by noise at periods exceeding about 10 years. The exception to this is possibly a linear trend found in some astrometric and space geodetic series. At frequencies above prograde +0.2 cpy (periods shorter than about 5 years), historical astrometric data may be of sufficient quality for comparisons with geophysical excitation time series. Even in the era of space geodesy, significant differences are found in long-period variations in published polar motion time series. Received: 27 March 2001 / Accepted: 15 October 2001     

Efficient gravity field recovery from GOCE gravity gradient observations

 An efficient algorithm is proposed for gravity field recovery from Gravity Field and Steady-State Ocean Circulation Explorer (GOCE) satellite gravity gradient observations. The mathematical model is formulated in the time domain, which allows the inclusion of realistic observational noise models. The algorithm combines the iterative solution of the normal equations, using a Richardson-type iteration scheme, with the fast computation of the right-hand side of the normal equations in each iteration step by a suitable approximation of the design matrix. The convergence of the iteration is investigated, error estimates are provided, and the unbiasedness of the method is proved. It is also shown that the method does not converge to the solution of the normal equations. The performance of the approach for white noise and coloured noise is demonstrated along a simulated GOCE orbit up to spherical harmonic degree and order 180. The results also indicate that the approximation error may be neglected. Received: 30 November 1999 / Accepted: 31 May 2000     

Analysis of the first year of Earth rotation parameters with a sub-daily resolution gained at the CODE processing center of the IGS

 The solutions of the CODE Analysis Center submitted to the IGS, the International Global Position System (GPS) Service for Geodynamics, are based on three days of observation of about 80–100 stations of the IGS network. The Earth rotation parameters (ERPs) are assumed to vary linearly over the three days with respect to an a priori model. Continuity at the day boundaries as well as the continuity of the first derivatives are enforced by constraints. Since early April 1995 CODE has calculated a new ERP series with an increased time resolution of 2 hours. Again continuity is enforced at the 2-hours-interval boundaries. The analysis method is described, particularly how to deal with retrograde diurnal terms in the ERP series which may not be estimated with satellite geodetic methods. The results obtained from the first year of data covered by the time series (time interval from 4 April 1995 to 30 June 1996) are also discussed. The series is relatively homogeneous in the sense of the used orbit model and the a priori model for the ERPs. The largest source of excitation at daily and sub-daily periods is likely to be the effect of the ocean tides. There is good agreement between the present results and Topex/Poseidon ocean tide models, as well as with models based on Very Long Baseline Interferometry (VLBI) and Satellite Laser Ranging (SLR) data. Non-oceanic periodic variations are also observed in the series. Their origin is most probably a consequence of the GPS solution strategy; other possible sources are the atmospheric tides. Received: 13 July 1999 / Accepted: 21 March 2000     

Gravity/magnetic potential of uneven shell topography

 A fast spherical harmonic approach enables the computation of gravitational or magnetic potential created by a non-uniform shell of material bounded by uneven topographies. The resulting field can be evaluated outside or inside the sphere, assuming that density of the shell varies with latitude, longitude, and radial distance. To simplify, the density (or magnetization) source inside the sphere is assumed to be the product of a surface function and a power series expansion of the radial distance. This formalism is applied to compute the gravity signal of a steady, dry atmosphere. It provides geoid/gravity maps at sea level as well as satellite altitude. Results of this application agree closely with those of earlier studies, where the atmosphere contribution to the Earth's gravity field was determined using more time-consuming methods. Received: 14 August 2000 / Accepted: 19 March 2001     

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     

Seasonal sea level change from TOPEX/Poseidon observation and thermal contribution

Seasonal steric sea-level change due to temperature variation in the mixing layer is assessed using space-measured sea-surface temperature data and historical in situ temperature measurements. The results are compared with TOPEX/Poseidon satellite altimeter measurement at different large spatial scales. It is indicated that thermal effect accounts for much of the observed seasonal variability, especially when averaging over zonal regions. Some regional seasonal patterns of sea-level anomalies in the tropical oceans are well represented by the thermal model prediction. Systematic differences are shown between TOPEX/Poseidon observation and thermal contribution at a 1–2 cm level. The potential causes for these differences are discussed, including water mass exchanges among the atmosphere, land, and oceans, and error sources in the steric result and geophysical corrections applied in TOPEX/Poseidon data. Received: 25 September 1998 / Accepted: 13 July 1999