Visual walkthrough — Orbital perturbations — J2 effect (oblateness), derivation of nodal precession

This page rebuilds the headline result of the parent note from nothing. No symbol is used before it has both a plain-word meaning and a picture. We are chasing one equation, and by the end you will see every piece of it:

$\dot{\Omega }=-\frac{3}{2}n{J}_2{\left(\frac{R_E}{p}\right)}^2\cos i$

Let us decode the destination in words before we walk toward it. The dot over $\Omega$ (a Greek capital "Omega") means "how fast $\Omega$ changes each second" — a rate. $\Omega$ itself is an angle that says which way the tilted orbit ring is pointed. So the whole equation is a single question: how fast does the orbit ring swivel? Everything below earns each symbol on the right-hand side, one step at a time.

---

Step 1 — What is an orbit, and what are we allowed to change?

WHAT. A satellite loops around Earth on an ellipse. To pin down that loop in space we need a few numbers. Four of them describe size and shape and tilt; two describe which way the loop points.

WHY. Before we can ask "how fast does it swivel," we must know what swivels. The bulge, it turns out, freezes the size, shape and tilt but rotates the two pointing angles. So we must first meet all six.

PICTURE. In the figure, the flat gray disk is Earth's equatorial plane — the plane of Earth's fat belly. The blue ellipse is the orbit, tilted up out of that plane.

The line where the tilted orbit slices the equatorial plane is the line of nodes. The point where the satellite crosses that line going northward is the ascending node — the orange dot. When $\Omega$ changes, that whole line of nodes rotates. That rotation is nodal precession. Everything from the Two-body problem & Keplerian orbits gives us these six numbers; here we ask what a bulge does to them.

---

Step 2 — Where does the satellite sit? Latitude from geometry

WHAT. We introduce $\varphi$, the geocentric latitude — the angle of the satellite above or below the equatorial plane, seen from Earth's center. We will show that $\sin \varphi =\sin i\sin u,u=\omega +\nu .$

WHY. The bulge only cares whether the satellite is above or below the equator — that is exactly what $\varphi$ measures. We need $\varphi$ written in terms of things we already have ($i$, and how far around the loop we are). We cannot use $\varphi$ in any formula until we know what it depends on.

PICTURE. The figure shows the spherical triangle sitting on the celestial sphere: one corner at the ascending node, one at the satellite, one at the point on the equator directly "below" the satellite.

Here $u=\omega +\nu$ is the argument of latitude — a single angle measuring "how far around the loop from the ascending node." Reading the figure:

$\underset{\text{height above equator}}{\underbrace{\sin \varphi }}=\underset{\text{how steep the ramp is}}{\underbrace{\sin i}}\cdot \underset{\text{how far up the ramp we’ve walked}}{\underbrace{\sin u}}$

This is pure spherical trigonometry — the same triangle you meet in the Two-body problem & Keplerian orbits.

---

Step 3 — The bulge adds one extra term to gravity

WHAT. A perfect sphere pulls with potential $-\mu /r$, where $\mu =G{M}_E$ (Earth's mass times the gravitational constant) and $r$ is the distance from center. The bulge adds a correction. The whole thing: $U(r,\varphi )=-\frac{\mu }{r}\left[1-J_2{\left(\frac{R_E}{r}\right)}^2\frac{3{\sin }^2\varphi -1}{2}\right].$

WHY. Gravity from a lumpy Earth is built by adding correction shells — the Spherical harmonics of gravity fields expansion. $J_2$ is the biggest correction: it encodes the equatorial bulge. We keep only this one because for low orbits it dominates all others.

PICTURE. The figure contrasts a round Earth (perfectly circular pull, gray) with the true squashed Earth (extra outward pull near the equator, orange arrows).

Term by term: $U=-\frac{\mu }{r}\left[1\underset{\begin{array}{c}\text{bulge size,}\\ \approx 1.08\times 10^{-3}\end{array}}{\underbrace{-J_2}}\underset{\begin{array}{c}\text{fades fast}\\ \text{with distance}\end{array}}{\underbrace{{\left(\frac{R_E}{r}\right)}^2}}\underset{\begin{array}{c}\text{latitude shape}\\ =P_2(\sin \varphi )\end{array}}{\underbrace{\frac{3{\sin }^2\varphi -1}{2}}}\right]$

• $R_E$ = Earth's equatorial radius. The ratio $(R_E/r{)}^2$ says the bulge matters most when you skim close.
• $\frac{1}{2}(3{\sin }^2\varphi -1)$ is the Legendre polynomial $P_2$ (see Legendre polynomials). It is positive near the poles and negative near the equator — the mathematical fingerprint of "fat belly, flat poles."

---

Step 4 — Feed geometry into the disturbing function

WHAT. Substitute Step 2's $\sin \varphi =\sin i\sin u$ into $R$: $R=\frac{\mu J_2{R}_E^2}{2r^3}\left(3{\sin }^{2}i{\sin }^{2}u-1\right).$

WHY. $R$ was written using latitude $\varphi$, but latitude is not one of our six orbit numbers — $i$ and $u$ are. To use the machinery of the Lagrange planetary equations we must express $R$ purely in orbital quantities. This substitution is that translation.

PICTURE. The figure plots $R$ over one loop (as $u$ runs $0\to 360°$): it oscillates up and down twice per orbit, because ${\sin }^{2}u$ has two humps.

$R=\frac{\mu J_2{R}_E^2}{2r^3}(3\underset{\text{tilt, fixed}}{\underbrace{{\sin }^{2}i}}\underset{\text{position, wiggling}}{\underbrace{{\sin }^{2}u}}-1)$

Notice the structure: one part depends on the fixed tilt $i$, the other on the rapidly changing position $u$. That split is the key to the next step.

---

Step 5 — Average over one loop: keep only what accumulates

WHAT. Replace the wiggling quantities by their averages over one full orbit. Two facts do the work: $⟨{\sin }^{2}u⟩=\frac{1}{2},⟨\frac{1}{r^3}⟩=\frac{1}{a^3(1-e^2{)}^{3/2}}.$ The result is the mean disturbing function $\bar{R}$, proportional to $(3{\sin }^{2}i-2)$.

WHY. In one loop the wiggles push the orbit one way, then pull it back — they cancel. What is left over, the non-cancelling constant, is what quietly builds up loop after loop. That leftover is the secular drift. We chase only the drift, so we average.

PICTURE. The figure shows the wiggly $R(u)$ curve of Step 4 with a flat dashed line through its mean — the wiggles above and below the line have equal area, so only the dashed level survives.

$⟨{\sin }^{2}u⟩=\frac{1}{2}\Rightarrow (3{\sin }^{2}i\cdot \frac{1}{2}-1)=\frac{1}{2}(3{\sin }^{2}i-2).$

• $⟨\cdot ⟩$ means "average over one orbit."
• The $1/r^3$ average brings in $p=a(1-e^2)$, the semi-latus rectum — the orbit's "effective radius" for these purposes.
• $n=\sqrt{\mu /a^3}$ is the mean motion: the average angular speed, i.e. $2\pi$ divided by the orbital period. It appears because averaging over time weights positions by how long the satellite lingers there.

$\bar{R}\propto \frac{n^2{J}_2{R}_E^2}{(1-e^2{)}^{3/2}}(3{\sin }^{2}i-2).$

---

Step 6 — The lever that turns tilt into swivel

WHAT. Lagrange's planetary equation for the node is $\frac{d\Omega }{dt}=\frac{1}{n{a}^2\sqrt{1-e^2}\sin i}\frac{\partial \bar{R}}{\partial i}.$

WHY. This is angular-momentum bookkeeping in disguise. The node swivels precisely when the averaged energy $\bar{R}$ would change if you tilted the orbit a little. That "change per tilt" is the partial derivative $\partial \bar{R}/\partial i$ — it answers "how sensitive is the bulge energy to inclination?" A derivative is the right tool because we want a rate of change with respect to a small nudge, which is exactly what a derivative measures.

PICTURE. The figure shows the orbit ring as a gyroscope: the bulge applies a torque (red arrow) perpendicular to both the spin axis and the equator, forcing the ring to precess sideways rather than tumble.

Now differentiate. Since $\bar{R}\propto (3{\sin }^{2}i-2)$: $\frac{\partial }{\partial i}(3{\sin }^{2}i-2)=6\sin i\cos i.$

Slot it back and watch the cancellation: $\dot{\Omega }\propto \frac{1}{\sin i}\cdot \underset{\text{from }{\partial }_i\bar{R}}{\underbrace{6\sin i\cos i}}=6\cos i.$

The $\sin i$ in the denominator cancels the $\sin i$ from the derivative, and only $\cos i$ survives. Collecting every constant carefully gives the boxed result:

$\boxed{\dot{\Omega }=-\frac{3}{2}n{J}_2{\left(\frac{R_E}{p}\right)}^2\cos i}$

$\dot{\Omega }=\underset{\text{sign + size}}{\underbrace{-\frac{3}{2}}}\underset{\text{orbit speed}}{\underbrace{n}}\underset{\text{bulge}}{\underbrace{J_2}}\underset{\text{how deep in bulge field}}{\underbrace{{\left(\frac{R_E}{p}\right)}^2}}\underset{\text{tilt factor}}{\underbrace{\cos i}}$

Every symbol on the right has now been seen: $n$ from the loop timing (Step 5), $J_2$ and $R_E/p$ from the bulge potential (Steps 3, 5), and the all-important $\cos i$ from the lever cancellation (this step).

---

Step 7 — Every case, no gaps: reading $\cos i$

WHAT. Walk through every possible tilt and check the sign and size of $\dot{\Omega }$.

WHY. The Contract: the reader must never meet a scenario we did not show. $\cos i$ changes sign, so the physics flips character — we must cover it all.

PICTURE. The figure plots $\dot{\Omega }$ versus inclination $i$ from $0°$ to $180°$, marking each regime.

Inclination $i$	$\cos i$	$\dot{\Omega }$	Meaning
$0°$ (equatorial, prograde)	$+1$	most negative	ring swivels westward fastest
between $0°$ and $90°$	positive	negative	westward, slowing as $i$ rises
$90°$ (polar)	$0$	zero	ring does not swivel at all
between $90°$ and $180°$	negative	positive	ring swivels eastward
$180°$ (equatorial, retrograde)	$-1$	most positive	eastward, fastest

---

Step 8 — The companion: perigee rotation and the critical angle

WHAT. The same machinery, but differentiating for $\omega$ instead of $\Omega$, gives $\dot{\omega }=\frac{3}{4}n{J}_2{\left(\frac{R_E}{p}\right)}^2\left(5{\cos }^{2}i-1\right).$

WHY. We include it because it hides a second special inclination, and because it shows the method is general — not a one-trick derivation.

PICTURE. The figure plots the tilt factor $(5{\cos }^{2}i-1)$; it crosses zero at $i=63.4°$ and $116.6°$.

$5{\cos }^{2}i-1=0\Rightarrow {\cos }^{2}i=\frac{1}{5}\Rightarrow i=63.4°\text{ or }116.6°.$

At this critical inclination the perigee stops rotating — it "freezes." Molniya orbits exploit exactly this so their apogee hangs over the northern hemisphere for hours. See Molniya & critical inclination orbits.

---

The one-picture summary

This single figure compresses the whole derivation: bulge → off-center pull (Step 3) → torque on the tilted ring (Step 6) → averaged over a loop (Step 5) → steady swivel of the node scaled by $\cos i$ (Steps 6–7).

Recall Feynman retelling of the whole walkthrough

Picture a slightly squashed Earth — fat around the middle. A satellite rides a tilted ring around it. Whenever the satellite drifts above or below that fat belly, the extra belly-mass tugs it a little sideways, back toward the equator (Step 3). That sideways tug is not straight toward Earth's center, so it twists the ring — a torque (Step 6). Over one full loop, most of the twisting cancels out; only a tiny steady leftover survives (Step 5). That leftover slowly turns the ring's pointing direction — the node angle $\Omega$ drifts. How fast? It depends on the orbit's speed $n$, the size of the bulge $J_2$, how close you skim ($R_E/p$), and above all the tilt through $\cos i$ (Step 6). Straight-over-the-poles ($i=90°$): no turn at all. Flat by the equator: fastest turn, westward. Lean past vertical ($i>90°$): the turn flips eastward — and if you tune that eastward turn to match Earth's yearly loop around the Sun, your orbit keeps the same face to the Sun forever (Step 7). The very same idea, applied to the ellipse's perigee instead of the ring, gives a magic tilt of $63.4°$ where the ellipse stops rotating (Step 8). One bulge, two beautiful design tricks.

Recall Quick self-test

Which orbit tilt gives zero nodal precession? ::: $i=90°$, the polar orbit, because $\cos 90°=0$. Why does the $\sin i$ from the derivative not appear in $\dot{\Omega }$? ::: It cancels the $\sin i$ in the denominator of Lagrange's equation, leaving pure $\cos i$. What makes retrograde ($i>90°$) orbits special for Sun-synchronicity? ::: $\cos i<0$ flips $\dot{\Omega }$ positive (eastward), which can be matched to Earth's yearly motion around the Sun.