# Deformation of Earth caused by centrifugal force

The centrifugal force causes a distension at the equator. With a simple model we can estimate the dimension of this deformation ellipsoid. We start with a mass $$m$$ at the point $$P$$ on the surface of the deformed earth at the height $$h$$ over the non-deformed earth.

Assume all particles are in balanced state on the surface. Show that the height of the particle can be described as: $$h = \frac{R^2\omega^2}{6g}(3sin^2 \theta-2).$$

I received some tips how to approach this problem: $$U_{total} = mgh - \frac{m}{2}\omega^2x^2 + C$$ Determine h out of the potential energy, when the particle is in balance.

Use spherical coordinates.

Assume the volume of the deformed earth is the same as the volume of the non deformed earth.

Integrate over the spherical surface.(that's the thing I don't understandt)

Here's a picture of the problem:

first step is easy: $$mgh = \frac{m}{2}\omega^2x^2 + C$$ with $$x=r*sin(\theta)$$ $$h=\frac{\omega^2R^2sin^2(\theta)}{2g} + B$$ with $$B = -C/mg$$.

From this point on I don't understand how to keep processing. We have solved for $$h$$. Why do I need to integrate? And far more confusing, why integrating over the spherical surface?! What's the principle behind this?

• Note that this problem is neglecting an important effect: the deformed Earth's gravitational quadrupole moment, which accounts for an extra factor $5/2$, cf. this Phys.SE post. Sep 1, 2017 at 18:21
• Does anyone have an idea left? Sep 2, 2017 at 16:35

In your comments in your answer to @Floris, you are doing the integral wrong. The volume element is $dV=r^2\sin\theta\,\mathrm dr\,\mathrm d\theta\,\mathrm d\phi$ but because $r=R+h\approx R$, you can approximate (after taking the radial integral and azimuthal integrals) it as $2\pi R^2h\sin\theta\,\mathrm d\theta$. You were missing that factor of $\sin\theta$. Using $A=\frac{\omega^2R^2}{2g}$,
$$\int_0^\pi (A\sin^2\theta+B)\sin\theta\,\mathrm d\theta=\frac{4}{3}A+2B$$
The difference in volume from a sphere is therefore $4\pi R^2\left(\frac{\omega^2R^2}{3g}+B\right)$. Setting that equal to $0$ gets you a value of $-\frac{\omega^2R^2}{3g}$ for $B$, which is the result you're looking for.
• I spent way too long trying to figure out the value of B by integrating $\int_R^{R+h}r^2\,\mathrm dr$ without using any approximations. The value of that integral is $R^2h+Rh^2+\frac{1}{3}h^3$, but if we're already using $F=mg$, there's no point going beyond lowest order in $h$, giving that value of $R^2h$ I used in the answer. Sep 5, 2017 at 4:37