Abstract
New algorithms have been derived for computing terrain connections, all components of the attraction of the topography at the topographic surface and the gradients of these attractions. These algorithms utilize fast Fourier transforms, but, in contrast to methods currently in use, all divergences of the integrals are removed during the analysis. Sequential methods employing a smooth intermediate reference surface have been developed to avoid the very large transforms necessary when making computations at high resolution over a wide area.
A new method for the numerical solution of Molodensky's problem has been developed to mitigate the convergence difficulties that occur at short wavelengths with methods based on a Taylor series expansion. A trial field on a level surface is continued analytically to the topographic surface, and compared with that predicted from gravity observations. The difference is used to compute a correction to the trial field and the process iterated. Special techniques are employed to speed convergence and prevent oscillations.
Three different spectral methods for fitting a point-mass set to a gravity field given on a regular grid at constant elevation are described. Two of the methods differ in the way that the spectrum of the point-mass set, which extends to infinite wave number, is matched to that of the gravity field which is band-limited. The third method is essentially a space-domain technique in which Fourier methods are used to solve a set of simultaneous equations.
Similar content being viewed by others
References
R.N. BRACEWELL: The Fourier transformation and its applications. 2nd ed., McGraw-Hill, New York, 1978.
C.W. CAMPBELL: Geometric interpretations of the Discrete Fourier Transformation (DFT). NASA Technical paper 2332, 1984.
R. FORSBERG and C.C. TSCHERNING. The use of height data in gravity field approximation by collocation. Jour. Geophys. Res., 86, pp. 7843–7854, 1981.
R. FORSBERG. A study of terrain reductions, density anomalies and geophysical inversion methods in gravity field modeling. Report #355, Department of Geodetic Science and Surveying Ohio State University, 1984.
R. FORSBERG. Gravity field terrain effect computations by FFT. Bull. Geod., 59, pp. 342–360, 1985.
M.G. SIDERIS: A Fast Fourier Transform method for computing terrain corrections. Manuscripta Geodaetica, 10, pp. 66–73, 1985.
M.G. SIDERIS and K.P. SCHWARZ: Solving Molodensky's series by fast Fourier transform techniques. Bull. Geod., 60, pp. 51–63, 1986.
M.G. SIDERIS: On the application of spectral techniques to the gravimetric problem. Invited paper to the XIX IAG/IUGG General Assembly, Vancouver, B.C., Aug. 9–22, 1987.
M.G. SIDERIS and K.P. SCHWARZ: Advances in the numerical solution of the linear Molodensky problem. Bull. Geod., 62, pp. 59–69, 1988.
H. SÜNKEL: Cardinal interpolation. Report #312, Department of Geodetic Science and Surveying, Ohio State University, 1981.
A.A. VASSILIOU. The use of spectral methods for the spatial modeling of gravity data. Manuscripta Geodaetica, 10, pp. 235–244, 1985.
Author information
Authors and Affiliations
Rights and permissions
About this article
Cite this article
Harrison, J.C., Dickinson, M. Fourier transform methods in local gravity modeling. Bull. Geodesique 63, 149–166 (1989). https://doi.org/10.1007/BF02519148
Received:
Accepted:
Issue Date:
DOI: https://doi.org/10.1007/BF02519148