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
... ps. 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 from the conventional relationship for E through the Green's identity. The first component expresses some minimum amount of the work W and the second component represents a deviation from 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 E min W W Δ min W H 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×10 39 ergs of the energy E. Thus, we get a perfect agreement between = −2.5073×10 Gauss E 39 ergs, = E −2.4910×10 39 ergs derived from the piecewise Roche's density, = −2.4884×10 PREM E 39 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 20 C 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.