Construction of the normal equations matrix for modeling of local gravitational field

: pp. 29 - 34
Received: December 15, 2013
Revised: March 24, 2014
Department of Geodesy, Institute of Geodesy, Lviv Polytechnic National University
Department of Higher geodesy and astronomy, Lviv polytechnic National University

We consider the method of constructing the local gravity field using nonorthogonal basic functions, which are solution of the Laplace equation in spherical cap or spherical segment. This approach involves using of associated Legendre functions of integer degree and noninteger order. These functions form two sets of functions. They are mutually orthogonal over the spherical cap in each set. Thus, for using both of these sets of functions it is traditionally used least squares method. However, for higher orders it is quite difficult to compute eigenvalues of these functions. In such case it is possible to project the initial data on the hemisphere and to use associated Legendre functions of integer degree and integer order. The properties of these functions are similar to properties of functions on the spherical cap. Traditionally, initial data is selected in the nodes of grid, especially for fast computations. There are many kinds of uniform grids, which allow to speed up the process of computation the unknown harmonic coefficients. Among these grids it is possible to allocate the geographical grid, Gauss grid and others. Thus, grid is developed to accommodate the initial data and is defined its basic properties in the segment of sphere and hemisphere . Using the properties of grid technique for computing the matrix of normal equations is obtained, which leads to a time reducing procedure. Also formulas for computations of unknown coefficients are obtained which allow to switch from the inversion of matrix with order α² to matrix with order α. The proposed algorithm for the constructing of the normal equations matrix and determination of harmonic coefficients of the local gravitational field leads to a time reducing procedure without degradation of accuracy.

  1. Dzhuman B.B. Pro pobudovu modeli lokalnoho hravitatsiinoho polia [On the constraction of local gravitational field model]. Heodynamika, 2013, no.1(14), pp.29–33.
  2. Lukianchenko Yu.O. Pobudova normalnykh rivnian dlia opratsiuvannia danykh misii GOCE [Constraction of normal equations mission for processing data GOCE]. Heodynamika, 2013, no.1(14), pp.34–37.
  3. De Santis A. Conventional spherical harmonic analysis for regional modeling of the geomagnetic feld. Geophys. Res. Lett., 1992, issue 19, pp.1065–1067.
  4. De Santis A., Torta J. Spherical cap harmonic analysis : a comment on its proper use for local gravity field representation // J. of Geod., 1997, 71, p. 526–532.
  5. Eicker A. Gravity field refinement by radial basis functions from In-situ satellite data, 2008, p.137.
  6. Haines G.V. Computer programs for spherical cap harmonic analysis of potential and general felds // Comput. Geosci., 1988, 14, p.413–447.
  7. Haines G.V. Spherical cap harmonic analysis. J. Geophys. Res., 1985, issue 90, pp.2583–2591.
  8. Hofmann-Wellenhof B., Moritz H. Physical Geodesy, Wien New York: Springer Science + Busines Media, 2005, p.403.
  9. Hwang C., Chen S. Fully normalized spherical cap harmonics: application to the analysis of sea-level data from TOPEX/POSEIDON and ERS-1 // Geophys. J. Int., 1997, 129, p. 450–460.
  10. Hwang J., Gravity and geoid model in South Korea and its vicinity by spherical cap harmonic analysis / Hwang J., Han H., Han S., Kim K., Kim J., Kang M., Kim C. // J. of Geodynamics, 2012, 53, p. 27–33.
  11. Jiancheng L., Spherical cap harmonic expansion for local gravity field representation / Jiancheng L., Dingbo C., Jiancheng N. // Manuscr. Geod., 1995. – 20, – p. 265–277.
  12. NGA, Arctic Gravity Project, The National Imagery and Mapping Agency, 2008. Available at:
  13. Thebault E., Modeling the lithospheric magnetic field over France by means of revised spherical cap harmonic analysis (R-SCHA) / Thebault E., Mandea M., Schott J. // J. Of Geophys. Research, 2006. – 111. – 13 p.