# Clairaut's theorem

File:Gnuplot ellipsoid.svg Clairaut's theorem, published in 1743 by Alexis Claude de Clairaut in his Théorie de la figure de la terre, tirée des principes de l'hydrostatique, synthesized physical and geodetic evidence that the Earth is an oblate rotational ellipsoid. It is a general mathematical law applying to spheroids of revolution. It was initially used to relate the gravity at any point on the Earth's surface to the position of that point, allowing the ellipticity of the Earth to be calculated from measurements of gravity at different latitudes.

## Formula

Clairaut's formula for the acceleration of gravity g on the surface of a spheroid at latitude φ, was: $g = G \left[ 1 + \left(\frac{5}{2} m - f\right) \sin^2 \phi \right] \ ,$

where G is the value of the acceleration of gravity at the equator, m the ratio of the centrifugal force to gravity at the equator, and f the flattening of a meridian section of the earth, defined as: $f = \frac {a-b}{a} \ ,$

(where a = semimajor axis, b=semiminor axis ).

Clairaut derived the formula under the assumption that the Earth was composed of concentric coaxial spheroidal layers of constant density. This work was subsequently pursued by Laplace, who relaxed the initial assumption that surfaces of equal density were spheroids. Stokes showed in 1849 that the theorem applied to any law of density so long as the external surface is a spheroid of equilibrium. A history of the subject, and more detailed equations for g can be found in Khan.

The above expression for g has been supplanted by the Somigliana equation: $g = G \left[ \frac{1+k\sin^2 \phi}{\sqrt{1-e^2 \sin^2 \phi }} \right] \ ,$

where, for the Earth, G =9.7803267714 ms−2; k =0.00193185138639 ; e2 =0.00669437999013.

## Clairaut's relation

See main article Clairaut's relation

A formal mathematical statement of Clairaut's theorem is:

Let γ be a geodesic on a surface of revolution S, let ϝ be the distance of a point of S from the axis of rotation, and let ψ be the angle between γ and the meridians of S. Then ϝ sin ψ is constant along γ. Conversely, if ϝ sin ψ is constant along some curve γ in the surface, and if no part of γ is part of some parallel of S, then γ is a geodesic.

Andrew Pressley: Elementary Differential Geometry, p. 183

Pressley (p. 185) explains this theorem as an expression of conservation of angular momentum about the axis of revolution when a particle slides along a geodesic under no forces other than those that keep it on the surface.

## Geodesy

The spheroidal shape of the Earth is the result of the interplay between gravity and centrifugal force caused by the Earth's rotation about its axis. In his Principia, Newton proposed the equilibrium shape of a homogeneous rotating Earth was a rotational ellipsoid with a flattening f given by 1/230. As a result gravity increases from the equator to the poles. By applying Clairaut's theorem, Laplace was able to deduce from 15 gravity values that f = 1/330. A modern estimate is 1/298.25642. See Figure of the Earth for more detail.

For a detailed account of the construction of the reference Earth model of geodesy, see Chatfield.