Algorithm for constructing the subsoil density distribution function considering its value on the surface
Received: September 11, 2023
Department of Cartography and Geospatial Modeling, Institute of Geodesy, Lviv Polytechnic National University
Department of Cartography and Geospatial Modeling, Institute of Geodesy, Lviv Polytechnic National University
Department of Geodesy, Lviv Polytechnic National University

The conventional approach to constructing a three-dimensional distribution of the Earth's masses involves using Stokes constants incrementally up to a certain order. However, this study proposes an algorithm that simultaneously considers all of these constants, which could potentially provide a more efficient method. The basis for this is a system of equations obtained by differentiating the Lagrange function, which takes into account the minimum deviation of the three-dimensional mass distribution of the planet's subsoil from one-dimensional referential one. An additional condition, apart from taking into account the Stokes constants, for an unambiguous solution to the problem is to specify the value of the function on the surface of the ellipsoidal planet. It is possible to simplify the calculation process by connecting the indices of summation values in a series of expansions to their one-dimensional analogues in the system of linear equations. The study presents a control example illustrating the application of the given algorithm. In its implementation, a simplified variant of setting the density on the surface of the ocean is taken. The preliminary results of calculations confirm the expediency of this approach and the need to expand such a technique with other conditions for unambiguously solving the inverse problem of potential theory. Objectives. To create and implement the algorithm that takes into account the density of the planet’s subsoil on its surface. Method. The mass distribution function of the planet's subsoil is represented by a decomposition into biorthogonal series, the coefficients of decomposition which are determined from a system of linear equations. The system of equations is obtained from the condition of minimizing the deviation function of the desired mass distribution from the initially determined two-dimensional density distribution (PREM reference model). Results. On the basis of the described algorithm, a three-dimensional model of the density distribution of subsoil masses in the middle of the Earth is obtained, which takes into account Stokes constants up to the eighth order inclusively and corresponds to the surface distribution of masses of the oceanic model of the Earth. Its concise interpretation is also presented.

  1. Bullen, K. (1962). Earth’s central density. Nature London, 973 p.
  2. Dziewonski, A. M, Hager, B. H, & O’Connell, R. J. (1977). Largescale heterogeneties in the lower manthle. J. Geophys. Res., 82 (2), 239–255.
  3. Dziewonski, A. M. (1984). Mapping the lower mantle: Determination of lateral heterogeneity in PVelocity up to degree and order 6. J. Geophys. Res. 89, 5929–5952.
  4. Furman V. (2018) Global models of seismic tomography in the study of structures Electronics and information technologies. Electronics and information technologies, 9, 48–62. (in Ukrainian).
  5. Fys, M., Brydun, A., & Yurkiv, M. (2019). Researching the influence of the mass distribution inhomogeneity of the ellipsoidal planet’s interior on its stokes constants. Geodynamics, 1 (26), 17–27.
  6. Fys, M., Yurkiv, M., Brydun, A., & Lozynskyi, V. (2016). One option of constructing three-dimensional distribution of the mass and its derivatives for a spherical planet Earth. Geodynamics, 21 (2), 36–44.
  7. Fys, M. M., Brydun, A. M., & Yurkiv, M. I. (2018). Method for approximate construction of three-dimensional mass distribution function and gradient of an elipsoidal planet based on external gravitational field parameters. Geodynamics, 2 (25), 27–36.
  8. Fys, M. M., Brydun, A. M., Yurkiv, M. I., Sohor, A. R., & Hubar, Yu. P. (2021). The gradient construction approach analysis of the three-dimensional mass distribution function of the ellipsoidal planet. Geodynamics, 2 (31), 29–40.
  9. Fys M. M., Gubar Y. P. (1999) On one combined method of constructing a model distribution of subsoil masses. Collection of scientific works of the scientific and technical conference, "Modern achievements of geodesy, geodynamics and geodetic production", Lviv, No. 1. P. 41–42. (in Ukrainian).
  10. Laske G., Masters G., Ma Z. and Pasyanos M., (2013) Update on CRUST1.0 - A 1-degree Global Model of Earth's Crust, Geophys. Res. Abstracts, 15, Abstract EGU2013-2658.
  11. Marchenko O., Tretyak K., Kylchitskiy A., Golubinka Yu., Marchenko D. & Tretyak N. (2012) Investigation of the gravitational field, ocean topography and crustal movements in the Antarctic region. Lviv, Lviv Polytechnic Publishing House, 306. (In Ukrainian)
  12. Martinenc Z., Pec K. (1986) Three-Dimensional Density Distribution Generating the Observed Gravity Field of planets: Part II. The Moon. Proc. Int. Symp. Figure of the Earth, the Moon and other Planets. Czechoslovakia, Prague, 1, P. 153–163.
  13. Meshcheryakov G. A. (1989) On the inverse (geophysical) problem of the theory potential. Theses reported conference "New approaches to the solution differential equations", Drohobych, P. 113. (in Russian).
  14. Meshcheryakov G.A. (1991) Problems theories potential and generalized Earth. Moscow: Nauka, P. 216. (in Russian).
  15. Moritz G. (1973) Computatson ellipsoidal Mass distributions. Department of Geodetic Science, The Ohio State University, 206, p. 20.
  16. Moritz G. (1994) Figure of the Earth: Theoretical geodesy and internal structure Earth. Kyiv, 240 p. (in Russian).
  17. Panning, M., & Romanowicz, B. (2006). A three-dimensional radially anisotropic model of shear velocity in the whole mantle. Geophys. J. Int., 167, 361-379.
  18. Pavlis N.K., Holmes SA, Kenyon SC (2008) An Earth Gravitational Model that degree 2160: EGM2008. EGU General Assembly. Geophysical Reaseach Abstract, No. 10, P. 2.
  19. Shcherbakov A. M. (1978) Volumetric density distribution of the Moon. Astronomical Bulletin, ХІІ. №2, P. 88–95. (in Russian).
  20. Strakhov V. N. (1977) About one common form decisions reverse tasks gravimetry. Add. Academy of Sciences of the USSR, Vol. 235. No. 6, pp. 1281-1284. (in Russian).
  21. Tarakanov Yu.A., Cherevko T.I. (1978) Interpretatsiya krupneyshih gravitatsionnyih anomaliy Zemli. Izv. AN SSSR. Fizika Zemli. Vol. 4. pp. 25–42.
  22. Tserklevych А., Zayats О. & Shylo Ye. (2016) Approximation of the physical surface of the Earth by biaxial and triaxial ellipsoid. Geodynamics, 1 (20), 40–49. (in Ukrainian).
  23. Tserklevych, A. L. (2005). Method of building a model of the distribution of the Earth’s tectonosphere density, consistent with the gravitational field and other geological-geophysical information. Geofizicheskiy Zhurnal, 27(2), 309–314. (in Ukrainian).
  24. Woodhouse, J. H., & Dziewonski A. M. (1986). Mapping the upper mantle: Three-dimensional modeling of Earth structure by inversion of seismic wave forms. J. Geophys. Res., 86, 5953–5986.
  25. Zayats O. S. (2006) Modeling of the topography, structure and internal gravitational field of the moon: author's abstract. Thesis Ph.D. technical Sciences: 05.24.01; Acting National Lviv Polytechnic University. Lviv, 18 p. (in Ukrainian).