Earth Oblateness Perturbation

Author: Tianjiang Shuo
Website: https://cislunarspace.cn

Definition

Earth oblateness perturbation arises from the non-spherical mass distribution of the Earth -- flattened at the poles and bulging at the equator -- causing gravitational deviation from the homogeneous sphere assumption. In near-Earth space, the second-degree zonal harmonic ( $J_2$ term) is the dominant perturbation, with the corresponding perturbing gravitational potential:

\Delta U(r,\varphi) = -\frac{1}{2}\frac{\mu_E}{r}J_2\left(\frac{R_E}{r}\right)^2(3\sin^2\varphi - 1)

where $J_2 = 1.08263 \times 10^{-3}$ is the Earth oblateness coefficient, $R_E$ is the Earth's radius, and $\varphi$ is the geocentric latitude.

Core Elements

Long-Term Effects of J2 on Orbital Elements

Considering only the $J_2$ long-period terms, the rates of change of each orbital element are:

\left\{\begin{array}{l} \dot{a} = 0, \quad \dot{e} = 0, \quad \dot{i} = 0 \\ \dot{\Omega} = -\frac{3J_2R_E^2}{2p^2}n\cos i \\ \dot{\omega} = \frac{3J_2R_E^2}{2p^2}n(2 - \frac{5}{2}\sin^2 i) \\ \dot{M} = n + \frac{3J_2R_E^2}{2p^2}n(1 - \frac{3}{2}\sin^2 i)\sqrt{1-e^2} \end{array}\right.

Ascending Node Right Ascension Drift (Orbital Precession)

Condition	Drift Direction	Physical Mechanism
Prograde orbit ( $i < 90°$ )	$\dot{\Omega} < 0$ , ascending node regresses westward	Gravitational torque from equatorial bulge causes angular momentum precession
Retrograde orbit ( $i > 90°$ )	$\dot{\Omega} > 0$ , ascending node progresses eastward	Same mechanism, opposite direction
Polar orbit ( $i = 90°$ )	$\dot{\Omega} = 0$	$\cos i = 0$ , unaffected by oblateness long-period terms

The drift rate depends on orbital altitude and eccentricity: smaller semi-major axes and eccentricities closer to 1 result in faster drift.

Argument of Perigee Drift

The rate of change of the argument of perigee $\dot{\omega}$ is constant and depends on $a$ , $e$ , and $i$ . When $\sin^2 i = 4/5$ (i.e., $i = 63.4°$ or $116.6°$ ), $\dot{\omega} = 0$ and the apse line does not rotate. This inclination is called the critical inclination.

Perturbation Decomposition

Earth oblateness perturbation can be decomposed into:

Long-period terms $R_C$ : Orbit-averaged values that cause monotonic changes in orbital elements
Short-period terms $R_S$ : Terms whose orbit-average over one period is zero, producing no secular accumulation

Application Value

Earth oblateness perturbation is the dominant perturbation source for near-Earth orbits (on the order of $10^{-3}$ ). Its long-period effects form the physical basis for designing sun-synchronous orbits and frozen orbits. The ascending node drift pattern is of great significance for constellation design, orbit maintenance, and space situational awareness.

References

郑伟, 安雪滢, 周祥, 何睿智. 空天飞行力学[M]. 国防科技大学, 2026.
刘林. 航天器轨道理论[M]. 国防工业出版社.