The estimation of the Earth’s gravitational potential energy *E *based on the given density distribution is considered. The global density model was selected as combined solution of the 3D continuous distribution and reference radial piecewise profile with basic density jumps as sampled for the PREM density. This model preserves the external gravitational potential from zero to second degree/order, the dynamical ellipticity, the planet’s flattening, and basic radial density-jumps. The rigorous error propagation of adopted density parameterization was derived to restrict a possible solution domain. Comparison of lateral density anomalies with estimated accuracy of density leads to values of the same order in uncertainties and density heterogeneities. As a result, radial-only density models were chosen for the computation of the potential energy *E*. *E*-estimates were based on the expression *E =* -(*W*min +Δ*W*) derived from the conventional relationship for *E *through the Green’s identity. The first component *W*min expresses some minimum amount of the work *W *and the second component Δ*W *represents a deviation from *W*min treated via Dirichlet’s integral on the internal potential. Relationships for the internal potential and *E*, including error propagation were developed for continuous and piecewise densities. Determination of *E *provides the inequality with two limits for *E*-values corresponding to different density models. The upper limit EH agrees with the homogeneous distribution. The minimum amount *E*Gauss corresponds to Gauss’ continuous radial density. All *E*-estimates were obtained for the spherical Earth since the ellipsoidal reduction gives two orders smaller quantity than the accuracy σ*E*= ±0.0025×1039 ergs of the energy *E*. Thus, we get a perfect agreement between *E*Gauss = −2.5073×1039 ergs, *E *= −2.4910×1039 ergs derived from the piecewise Roche’s density, *E*PREM = −2.4884×1039 ergs based on the PREM model, and *E-*values from simplest models separated into core and mantle only. Distributions of the internal potential and its first and second derivatives were derived for piecewise and continuous density models. Influence of the secular variation in the zonal coefficient *C*20 on global density changes is discussed using the adopted 3D continuous density model as restricted solution of the three-dimensional Cartesian moments problem inside the ellipsoid of revolution.

1. Idelson NI (1932). Theory of Potential with Applications to the Geophysics Problems. GTTI, Leningrag - Moscow. (in Russian).

2. Kurant R, Gilbert D (1951) Methods of the mathematical physics, 2nd edition GITTL, Moscow − Leningrad, (in Russian).

3. Marchenko O.N., Yarema N.P. (2003) Estimation of influence of uncertainties of principle fundamental constant on the accuracy of the density distribution, Proceedings of the international symposium "Modern Achievements of Geodetic Science and Industry", Lviv, April, 2003, pp.77–84 (in Ukrainian).

4. Mescheryakov GA. (1973) On estimation of some values characterizing the internal gravity field of the Earth. Geodesy, cartography and aerophotosurveying, Lvov, No. 17, pp. 34–40.

5. Mescheryakov GO. (1977) On the unique solution of the inverse problem of the potential theory. Reports of the Ukrainian Academy of Sciences. Kiev, Series A, No. 6, pp. 492–495

6. Mescheryakov GA. (1991) Problems of the potential theory and generalized Earth. Nauka, Moscow, 1991. 203 p. (in Russian).

7. Mescheryakov GA, Shopyak IN, Deyneka YuP. (1977) On the representation of a function inside the Earth's ellipsoid by means of the partial sum of a generalized Fourier series. Geodesy, cartography and aerophotosurveying, No 21, pp. 55–62, Lvov. (in Russian).

8. Mihlin SG (1968) Course of the mathematical physics, Nauka, Moscow, (in Russian).

9. Tichonov AN, Arsenin VY (1974) Methods of solution of ill-posed problem, Nauka, Moscow, (in Russian).

10. Bullen KE (1975) The Earth's Density. Chapman and Hall, London. 22. Marchenko AN (2009a) Current estimation of the Earth's mechanical and geometrical parameters. In M.G.Sideris (ed.), Observing our Changing Earth, Int. Assoc. of Geodesy Symp. 132. pp. 473–481, Springer Verlag, Berlin-Heidelberg 2009.

https://doi.org/10.1007/978-94-009-5700-8

11. Bursa M (1993) Distribution of gravitational potential energy within the Solar system, Earth, Moon and Planets, Vol. 62, pp. 149–159.

https://doi.org/10.1007/BF00572144

12. Bursa M, Krivsky L, Hovorkova O (1996) Gravitational potential energy of the Sun, Studia geoph, et geod. 40 (1996), pp. 1–8. 23. Marchenko AN (2009b) The Earth's global density distribution and gravitational potential energy. In M.G.Sideris (ed.), Observing our Changing Earth, Int. Assoc. of Geodesy Symp. 132. pp. 483–491, Springer Verlag, BerlinHeidelberg 2009.

13. Dziewonski AM, Anderson DL (1981) Preliminary reference Earth model. Physics of the Earth and Planetary Interiors, Vol. 25, pp. 297–356.

https://doi.org/10.1016/0031-9201(81)90046-7

14. Gauss KF (1877). Allgemeine Lehrsätze in Beziehung auf die im verkehrten Verhaltnisse des Quadrats der Entfernung wirkenden Anziehungs und Abstossungskrafte. Resultate

https://doi.org/10.1007/978-3-642-49319-5_11

15. Grafarend E, Engels J, Varga P (2000) The temporal variation of the spherical and Cartesian multipoles of the gravity field: the generalization MacCullagh representation. Journal of Geodesy, Vol. 74, pp. 519–530.

https://doi.org/10.1007/s001900000114

16. Groten E (2004) Fundamental parameters and current (2004) best estimates of the parameters of common relevance to astronomy, geodesy and geodynamics, Journal of Geodesy, Vol. 77, pp. 724–731.

https://doi.org/10.1007/s00190-003-0373-y

17. Heiskanen WA, Moritz H (1967) Physical Geodesy, W.H. Freeman, San Francisco.

18. Holota P (1995) Two branches of the Newton potential and geoid, Proceed. of the International Symposium No 113 "Gravity and Geoid", Graz, Austria, 1994, pp. 205–214, Springer-Verlag Berlin Heidelberg.

19. Kellogg OD (1929) Foundations of potential theory, Springer, Berlin, 1929.

https://doi.org/10.1007/978-3-642-90850-7

20. Marchenko AN (1999) Earth's radial profiles based on Legendre-Laplace law. Geodynamics, 1(2) 1999, pp. 1–6.

21. Marchenko AN (2000) Earth's radial density profiles based on Gauss' and Roche's distributions. Bolletino di Geodesia e Scienze Affini, Anno LIX, No.3, pp. 201–220.

22. Marchenko AN (2009a) Current estimation of the Earth's mechanical and geometrical parameters. In M.G.Sideris (ed.), Observing our Changing Earth, Int. Assoc. of Geodesy Symp. 132. pp. 473–481, Springer Verlag, Berlin-Heidelberg 2009.

23. Marchenko AN (2009b) The Earth's global density distribution and gravitational potential energy. In M.G.Sideris (ed.), Observing our Changing Earth, Int. Assoc. of Geodesy Symp. 132. pp. 483–491, Springer Verlag, BerlinHeidelberg 2009.

24. Marchenko AN, Schwintzer P (2003) Estimation of the Earth's tensor of inertia from recent global gravity field solutions, Journal of Geodesy, Vol. 76, pp. 495–509.

https://doi.org/10.1007/s00190-002-0280-7

25. McCarthy D, Petit G (2004) IERS Conventions (2003), IERS Technical Note, No.32, Verlag des Bundesamts fur Kartographie und Geodasie, Frankfurt am Main, 2004.

26. Maxwell JK (1881) A Treatise on Electricity and Magnetism. 2nd Edition, Oxford, Vol.1.

27. Moritz H (1990) The Figure of the Earth. Theoretical Geodesy and Earth's Interior, Wichmann, Karlsruhe.

28. Rubincam DP (1979) Gravitational potential energy of the Earth: A spherical harmonic approach, Journal of Geoph. Res. Vol. 84, No. B11, 6219–6225.

https://doi.org/10.1029/JB084iB11p06219

29. Thomson W, Tait PG (1879) Treatise on Natural Philosophy, Vol.2, Cambridge University Press.

30. Wermer J (1981) Potential theory, Springer, Berlin Heidelberg New York.

31. Yoder, CF, Williams JG, Dickey JO, Schutz BE, Eanes R.J, Tapley BD (1983) Secular variation of Earth's gravitational harmonic J2 coefficient from Lageos and nontidal acceleration of Earth rotation. Nature, 303, pp. 757–762.

https://doi.org/10.1038/303757a0