Gravity of Earth

A non-rotating perfect sphere of uniform mass density, or whose density varies solely with distance from the centre (spherical symmetry), would produce a gravitational field of uniform magnitude at all points on its surface. The Earth is rotating and is also not spherically symmetric; rather, it is slightly flatter at the poles while bulging at the Equator: an oblate spheroid. There are consequently slight deviations in the magnitude of gravity across its surface.

Gravity on the Earth's surface varies by around 0.7%, from 9.7639 m/s² on the Nevado Huascarán mountain in Peru to 9.8337 m/s² at the surface of the Arctic Ocean.[6] In large cities, it ranges from 9.7806 m/s² [7] in Kuala Lumpur, Mexico City, and Singapore to 9.825 m/s² in Oslo and Helsinki.

Latitude

The differences of Earth's gravity around the Antarctic continent.

The surface of the Earth is rotating, so it is not an inertial frame of reference. At latitudes nearer the Equator, the outward centrifugal force produced by Earth's rotation is larger than at polar latitudes. This counteracts the Earth's gravity to a small degree – up to a maximum of 0.3% at the Equator – and reduces the apparent downward acceleration of falling objects.

The second major reason for the difference in gravity at different latitudes is that the Earth's equatorial bulge (itself also caused by centrifugal force from rotation) causes objects at the Equator to be further from the planet's center than objects at the poles. Because the force due to gravitational attraction between two bodies (the Earth and the object being weighed) varies inversely with the square of the distance between them, an object at the Equator experiences a weaker gravitational pull than an object on one of the poles.

In combination, the equatorial bulge and the effects of the surface centrifugal force due to rotation mean that sea-level gravity increases from about 9.780 m/s² at the Equator to about 9.832 m/s² at the poles, so an object will weigh approximately 0.5% more at the poles than at the Equator.[2][12]

Altitude

The graph shows the variation in gravity relative to the height of an object above the surface

Earth vs Mars vs Moon gravity at elevation

Gravity decreases with altitude as one rises above the Earth's surface because greater altitude means greater distance from the Earth's centre. All other things being equal, an increase in altitude from sea level to 9,000 metres (30,000 ft) causes a weight decrease of about 0.29%. (An additional factor affecting apparent weight is the decrease in air density at altitude, which lessens an object's buoyancy.[13] This would increase a person's apparent weight at an altitude of 9,000 metres by about 0.08%)

It is a common misconception that astronauts in orbit are weightless because they have flown high enough to escape the Earth's gravity. In fact, at an altitude of 400 kilometres (250 mi), equivalent to a typical orbit of the ISS, gravity is still nearly 90% as strong as at the Earth's surface. Weightlessness actually occurs because orbiting objects are in free-fall.[14]

The effect of ground elevation depends on the density of the ground (see Slab correction section). A person flying at 9,100 m (30,000 ft) above sea level over mountains will feel more gravity than someone at the same elevation but over the sea. However, a person standing on the Earth's surface feels less gravity when the elevation is higher.

The following formula approximates the Earth's gravity variation with altitude:

${\displaystyle g_{h}=g_{0}\left({\frac {R_{\mathrm {e} }}{R_{\mathrm {e} }+h}}\right)^{2}}$

Where

• g_h is the gravitational acceleration at height h above sea level.
• R_e is the Earth's mean radius.
• g₀ is the standard gravitational acceleration.

The formula treats the Earth as a perfect sphere with a radially symmetric distribution of mass; a more accurate mathematical treatment is discussed below.

Depth

Gravity at different internal layers of Earth (1 = continental crust, 2 = oceanic crust, 3 = upper mantle, 4 = lower mantle, 5+6 = core, A = crust-mantle boundary)

Earth's radial density distribution according to the Preliminary Reference Earth Model (PREM).[15]

Earth's gravity according to the Preliminary Reference Earth Model (PREM).[15] Two models for a spherically symmetric Earth are included for comparison. The dark green straight line is for a constant density equal to the Earth's average density. The light green curved line is for a density that decreases linearly from center to surface. The density at the center is the same as in the PREM, but the surface density is chosen so that the mass of the sphere equals the mass of the real Earth.

See also: Shell theorem

An approximate value for gravity at a distance r from the center of the Earth can be obtained by assuming that the Earth's density is spherically symmetric. The gravity depends only on the mass inside the sphere of radius r. All the contributions from outside cancel out as a consequence of the inverse-square law of gravitation. Another consequence is that the gravity is the same as if all the mass were concentrated at the center. Thus, the gravitational acceleration at this radius is[16]

${\displaystyle g(r)=-{\frac {GM(r)}{r^{2}}}.}$

where G is the gravitational constant and M(r) is the total mass enclosed within radius r. If the Earth had a constant density ρ, the mass would be M(r) = (4/3)πρr³ and the dependence of gravity on depth would be

${\displaystyle g(r)={\frac {4\pi }{3}}G\rho r.}$

The gravity g′ at depth d is given by g′ = g(1 − d/R) where g is acceleration due to gravity on the surface of the Earth, d is depth and R is the radius of the Earth. If the density decreased linearly with increasing radius from a density ρ₀ at the center to ρ₁ at the surface, then ρ(r) = ρ₀ − (ρ₀ − ρ₁) r / R, and the dependence would be

${\displaystyle g(r)={\frac {4\pi }{3}}G\rho _{0}r-\pi G\left(\rho _{0}-\rho _{1}\right){\frac {r^{2}}{R}}.}$

The actual depth dependence of density and gravity, inferred from seismic travel times (see Adams–Williamson equation), is shown in the graphs below.

Local topography and geology

See also: Physical geodesy

Local differences in topography (such as the presence of mountains), geology (such as the density of rocks in the vicinity), and deeper tectonic structure cause local and regional differences in the Earth's gravitational field, known as gravitational anomalies.[17] Some of these anomalies can be very extensive, resulting in bulges in sea level, and throwing pendulum clocks out of synchronisation.

The study of these anomalies forms the basis of gravitational geophysics. The fluctuations are measured with highly sensitive gravimeters, the effect of topography and other known factors is subtracted, and from the resulting data conclusions are drawn. Such techniques are now used by prospectors to find oil and mineral deposits. Denser rocks (often containing mineral ores) cause higher than normal local gravitational fields on the Earth's surface. Less dense sedimentary rocks cause the opposite.

A map of recent volcanic activity and ridge spreading. The areas where NASA GRACE measured gravity to be stronger than the theoretical gravity have a strong correlation with the positions of the volcanic activity and ridge spreading.

There is a strong correlation between the gravity derivation map of earth from NASA GRACE with positions of recent volcanic activity, ridge spreading and volcanos: these regions have a stronger gravitation than theoretical predictions.

A plumb bob determines the local vertical direction

Gravity acceleration is a vector quantity, with direction in addition to magnitude. In a spherically symmetric Earth, gravity would point directly towards the sphere's centre. As the Earth's figure is slightly flatter, there are consequently significant deviations in the direction of gravity: essentially the difference between geodetic latitude and geocentric latitude. Smaller deviations, called vertical deflection, are caused by local mass anomalies, such as mountains.

Tools exist for calculating the strength of gravity at various cities around the world.[18] The effect of latitude can be clearly seen with gravity in high-latitude cities: Anchorage (9.826 m/s²), Helsinki (9.825 m/s²), being about 0.5% greater than that in cities near the equator: Kuala Lumpur (9.776 m/s²). The effect of altitude can be seen in Mexico City (9.776 m/s²; altitude 2,240 metres (7,350 ft)), and by comparing Denver (9.798 m/s²; 1,616 metres (5,302 ft)) with Washington, D.C. (9.801 m/s²; 30 metres (98 ft)), both of which are near 39° N. Measured values can be obtained from Physical and Mathematical Tables by T.M. Yarwood and F. Castle, Macmillan, revised edition 1970.[19]

If the terrain is at sea level, we can estimate, for the Geodetic Reference System 1980, ${\displaystyle g\{\phi \}}$, the acceleration at latitude ${\displaystyle \phi }$:

${\displaystyle {\begin{aligned}g\{\phi \}&=9.780327\,\,\mathrm {m} \cdot \mathrm {s} ^{-2}\,\,\left(1+0.0053024\,\sin ^{2}\phi -0.0000058\,\sin ^{2}2\phi \right),\\&=9.780327\,\,\mathrm {m} \cdot \mathrm {s} ^{-2}\,\,\left(1+0.0052792\,\sin ^{2}\phi +0.0000232\,\sin ^{4}\phi \right),\\&=9.780327\,\,\mathrm {m} \cdot \mathrm {s} ^{-2}\,\,\left(1.0053024-0.0053256\,\cos ^{2}\phi +0.0000232\,\cos ^{4}\phi \right),\\&=9.780327\,\,\mathrm {m} \cdot \mathrm {s} ^{-2}\,\,\left(1.0026454-0.0026512\,\cos 2\phi +0.0000058\,\cos ^{2}2\phi \right)\end{aligned}}}$

This is the International Gravity Formula 1967, the 1967 Geodetic Reference System Formula, Helmert's equation or Clairaut's formula.[20]

An alternative formula for g as a function of latitude is the WGS (World Geodetic System) 84 Ellipsoidal Gravity Formula:[21]

${\displaystyle g\{\phi \}=\mathbb {G} _{e}\left[{\frac {1+k\sin ^{2}\phi }{\sqrt {1-e^{2}\sin ^{2}\phi }}}\right],\,\!}$

where,

• ${\displaystyle a,\,b}$ are the equatorial and polar semi-axes, respectively;
• ${\displaystyle e^{2}=1-(b/a)^{2}}$ is the spheroid's eccentricity, squared;
• ${\displaystyle \mathbb {G} _{e},\,\mathbb {G} _{p}\,}$ is the defined gravity at the equator and poles, respectively;
• ${\displaystyle k={\frac {b\,\mathbb {G} _{p}-a\,\mathbb {G} _{e}}{a\,\mathbb {G} _{e}}}}$ (formula constant);

then, where ${\displaystyle \mathbb {G} _{p}=9.8321849378\,\,\mathrm {m} \cdot \mathrm {s} ^{-2}}$,[21]

${\displaystyle g\{\phi \}=9.7803253359\,\,\mathrm {m} \cdot \mathrm {s} ^{-2}\left[{\frac {1+0.001931852652\,\sin ^{2}\phi }{\sqrt {1-0.0066943799901\,\sin ^{2}\phi }}}\right]}$.

where the semi-axes of the earth are:

${\displaystyle a=6378137.0\,\,{\mbox{m}}}$

${\displaystyle b=6356752.314245\,\,{\mbox{m}}}$

The difference between the WGS-84 formula and Helmert's equation is less than 0.68 μm·s⁻².

Further reductions are applied to obtain gravity anomalies (see: Gravity anomaly#Computation).

From the law of universal gravitation, the force on a body acted upon by Earth's gravitational force is given by

${\displaystyle F=G{\frac {m_{1}m_{2}}{r^{2}}}=\left(G{\frac {M_{\oplus }}{r^{2}}}\right)m}$

where r is the distance between the centre of the Earth and the body (see below), and here we take ${\displaystyle M_{\oplus }}$ to be the mass of the Earth and m to be the mass of the body.

Additionally, Newton's second law, F = ma, where m is mass and a is acceleration, here tells us that

${\displaystyle F=mg}$

Comparing the two formulas it is seen that:

${\displaystyle g=G{\frac {M_{\oplus }}{r^{2}}}}$

So, to find the acceleration due to gravity at sea level, substitute the values of the gravitational constant, G, the Earth's mass (in kilograms), m₁, and the Earth's radius (in metres), r, to obtain the value of g:[22]

${\displaystyle g=G{\frac {M_{\oplus }}{r^{2}}}=6.674\cdot 10^{-11}\mathrm {{m}^{3}\,{kg}^{-1}{s}^{-2}} \times {\frac {6\times 10^{24}\mathrm {kg} }{(6.4\times 10^{6}\mathrm {m} )^{2}}}=9+{\frac {795}{2^{10}}}\approx 9.77637\,\mathrm {{m}\,{s}^{-2}} }$

This formula only works because of the mathematical fact that the gravity of a uniform spherical body, as measured on or above its surface, is the same as if all its mass were concentrated at a point at its centre. This is what allows us to use the Earth's radius for r.

The value obtained agrees approximately with the measured value of g. The difference may be attributed to several factors, mentioned above under "Variations":

• The Earth is not homogeneous
• The Earth is not a perfect sphere, and an average value must be used for its radius
• This calculated value of g only includes true gravity. It does not include the reduction of constraint force that we perceive as a reduction of gravity due to the rotation of Earth, and some of gravity being counteracted by centrifugal force.

There are significant uncertainties in the values of r and m₁ as used in this calculation, and the value of G is also rather difficult to measure precisely.

If G, g and r are known then a reverse calculation will give an estimate of the mass of the Earth. This method was used by Henry Cavendish.

The measurement of Earth's gravity is called gravimetry.

Satellite measurements

This section is an excerpt from Gravimetry § Satellite gravimetry.[edit]

Gravity anomaly map from GRACE

Currently, the static and time-variable Earth's gravity field parameters are being determined using modern satellite missions, such as GOCE, CHAMP, Swarm, GRACE and GRACE-FO.[23][24] The lowest-degree parameters, including the Earth's oblateness and geocenter motion are best determined from Satellite laser ranging.[25]

Large-scale gravity anomalies can be detected from space, as a by-product of satellite gravity missions, e.g., GOCE. These satellite missions aim at the recovery of a detailed gravity field model of the Earth, typically presented in the form of a spherical-harmonic expansion of the Earth's gravitational potential, but alternative presentations, such as maps of geoid undulations or gravity anomalies, are also produced.

The Gravity Recovery and Climate Experiment (GRACE) consists of two satellites that can detect gravitational changes across the Earth. Also these changes can be presented as gravity anomaly temporal variations. The Gravity Recovery and Interior Laboratory (GRAIL) also consisted of two spacecraft orbiting the Moon, which orbited for three years before their deorbit in 2015.

• Altitude gravity calculator
• GRACE – Gravity Recovery and Climate Experiment
• GGMplus high resolution data (2013)
• Geoid 2011 model Potsdam Gravity Potato