Gravitational potential energy of the earth - A spherical harmonic approach

A spherical harmonic equation for the gravitational potential energy of the earth is derived for an arbitrary density distribution by conceptually bringing in mass-elements from infinity and building up the earth shell upon spherical shell. The zeroth degree term in the spherical harmonic expansion agrees with the usual expression for the energy of a radial density distribution. The second degree terms give a maximum nonhydrostatic energy in the crust and mantle of -2.77 x 10 to the 29th ergs, an order of magnitude below McKenzie's (1966) estimate. McKenzie's result stems from mathematical error. Our figure is almost identical with Kaula's (1963) estimate of the minimum shear strain energy in the mantle, a not unexpected result on the basis of the virial theorem. If the earth is assumed to be a homogeneous viscous oblate spheroid relaxing to an equilibrium shape, then a lower limit to the mantle viscosity of 1.3 x 10 to the 20th P is found by assuming that the total geothermal flux is due to viscous dissipation of energy. This number is almost six orders of magnitude below MacDonald's (1966) estimate of the viscosity and removes his objection to convection. If the nonequilibrium figure is dynamically maintained by the earth acting as a heat engine at 1% efficiency, then the viscosity is 10 to the 22nd P, a number preferred by Cathles (1975) and Peltier and Andrew (1976) as the viscosity of the mantle.