# Nodal precession

Nodal precession is the precession of the orbital plane of a satellite around the rotational axis of an astronomical body such as Earth. This precession is due to the non-spherical nature of a rotating body, which creates a non-uniform gravitational field. The following discussion relates to low Earth orbit of artificial satellites, which have no measurable effect on the motion of Earth. The nodal precession of more massive, natural satellites like the Moon is more complex.

Around a spherical body, an orbital plane would remain fixed in space around the gravitational primary body. However, most bodies rotate, which causes an equatorial bulge. This bulge creates a gravitational effect that causes orbits to precess around the rotational axis of the primary body.

The direction of precession is opposite the direction of revolution. For a typical prograde orbit around Earth (that is, in the direction of primary body's rotation), the longitude of the ascending node decreases, that is the node precesses westward. If the orbit is retrograde, this increases the longitude of the ascending node, that is the node precesses eastward. This nodal progression enables heliosynchronous orbits to maintain a nearly constant angle relative to the Sun.

## Description

Equatorial bulge torques a satellite orbit, leading to nodal precession

A non-rotating body of planetary scale or larger would be pulled by gravity into a sphere. Virtually all bodies rotate, however. The centrifugal force deforms the body so that it has an equatorial bulge. Because of the bulge the gravitational force on the satellite is not directly toward the center of the central body, but is offset toward the equator. Whichever hemisphere the satellite is in it is preferentially pulled slightly toward the equator. This creates a torque on the orbit. This torque does not reduce the inclination; rather, it causes a torque-induced gyroscopic precession, which causes the orbital nodes to drift with time.

## Equation

The rate of precession depends on the inclination of the orbital plane to the equatorial plane, as well as the orbital eccentricity.

For a satellite in a prograde orbit around Earth, the precession is westward (nodal regression), the node and satellite move in opposite directions.[1] A good approximation of the precession rate is

${\displaystyle \omega _{\mathrm {p} }=-{\frac {3}{2}}\cdot {\frac {{R_{\mathrm {E} }}^{2}}{\left(a\left(1-e^{2}\right)\right)^{2}}}J_{2}\omega \cos i}$

where

ωp is the precession rate (in rad/s)
RE is the body's equatorial radius (6378137 m for Earth)
a is the semi-major axis of the satellite's orbit
e is the eccentricity of the satellite's orbit
ω is the angular velocity of the satellite's motion (2π radians divided by its period in seconds)
i is its inclination (in degrees)
J2 is the body's second dynamic form factor (5C20 = 1.08262668×10−3 for Earth).

This last quantity is related to the oblateness as follows:

${\displaystyle J_{2}={\frac {2\varepsilon _{\mathrm {E} }}{3}}-{\frac {{R_{\mathrm {E} }}^{3}{\omega _{\mathrm {E} }}^{2}}{3GM_{\mathrm {E} }}}}$

where

εE is the central body's oblateness
RE is central body's equatorial radius (6378137 m for Earth)
ωE is the central body's rotation rate (7.292115×10−5 rad/s for Earth)
GME is the product of the universal constant of gravitation and the central body's mass (3.986004418×1014 m3/s2 for Earth).

The nodal progression of low Earth orbits is typically a few degrees per day to the west (negative). For a satellite in a circular (e = 0) 800 km altitude orbit at 56° inclination about Earth:

{\displaystyle {\begin{aligned}R_{\mathrm {E} }&=6.378\,137\times 10^{6}{\text{ m}}\\J_{2}&=1.082\,626\,68\times 10^{-3}\end{aligned}}}

The orbital period is 6052.4 s, so the angular velocity is 0.001038 rad/s. The precession is therefore

{\displaystyle {\begin{aligned}\omega _{\mathrm {p} }&=-{\frac {3}{2}}\cdot {\frac {6\,378\,137^{2}}{\left(7\,178\,137\left(1-0^{2}\right)\right)^{2}}}\cdot \left(1.082\,626\,68\times 10^{-3}\right)\cdot 0.001\,038\cdot \cos 56^{\circ }\\&=-7.44\times 10^{-7}{\text{ rad/s}}\end{aligned}}}

This is equivalent to −3.683° per day, so the orbit plane will make one complete turn (in inertial space) in 98 days.

The apparent motion of the sun is approximately +1° per day (360° per year / 365.2422 days per tropical year ≈ 0.9856473° per day), so apparent motion of the sun relative to the orbit plane is about 2.8° per day, resulting in a complete cycle in about 127 days. For retrograde orbits ω is negative, so the precession becomes positive. (Alternatively, ω can be thought of as positive but the inclination is greater than 90°, so the cosine of the inclination is negative.) In this case it is possible to make the precession approximately match the apparent motion of the sun, resulting in a heliosynchronous orbit.