Orbital perturbations — J2 effect (oblateness), derivation of nodal precession: Question bank

This page is a stress-test for your understanding of the J2 nodal-precession result. No heavy arithmetic here — every item pokes at a concept: a sign, a limiting case, a "why this and not that". Read the prompt, say your answer out loud, THEN reveal.

Before you start, keep the two headline results in view so every reveal has something to point at:

Definition The symbols on this page (so it stands alone)

$Ω$ = right ascension of the ascending node — the compass-angle at which the tilted orbit plane crosses the equator going north. When $Ω$ drifts, the whole orbit plane swivels. This is nodal precession.
$ω$ = argument of perigee — the angle inside the orbit plane from the ascending node to the point of closest approach (perigee). When $ω$ drifts, the ellipse rotates within its own plane.
$i$ = inclination — the tilt of the orbit plane relative to the equator.
$e$ = eccentricity (0 = circle), $a$ = semi-major axis (orbit size).
$μ = G M_{E}$ = Earth's gravitational parameter (mass $\times$ big- $G$ ).
$R_{E}$ = Earth's equatorial radius; $J_{2} \approx 1.0826 \times 1 0^{- 3}$ = the oblateness coefficient (size of the equatorial bulge).
$n = μ / a^{3}$ = mean motion (radians of orbit swept per second); $p = a (1 - e^{2})$ = semi-latus rectum.
$u = ω + ν$ = argument of latitude (how far around the orbit the satellite is, measured from the node); $ν$ = true anomaly; $ϕ$ = geocentric latitude.

The three figures below give you the pictures these symbols name; every reveal points back at them.

Recall The derivation in one screen (so the reveals have Steps 1–4 to point at)

This page's answers cite "Step 3" and "Lagrange's node equation" — here is the whole chain the parent note builds, written out so you can see it.

Step 1 — the disturbing function. The bulge's extra potential energy per mass, above the point-mass baseline, is $R = \frac{μ J _{2} R _{E}^{2}}{2 r ^{3}} (3 sin^{2} ϕ - 1) .$ The bracket is $2 P_{2} (sin ϕ)$ , the degree-2 Legendre polynomial.
Step 2 — geometry. Spherical trig of Picture 2 gives $sin ϕ = sin i sin u$ , so $R = \frac{μ J _{2} R _{E}^{2}}{2 r ^{3}} (3 sin^{2} i sin^{2} u - 1) .$
Step 3 — average over one orbit. Using $⟨ sin^{2} u ⟩ = \frac{1}{2}$ and $⟨ 1/ r^{3} ⟩ = 1/ [a^{3} (1 - e^{2})^{3/2}]$ , the short-period wiggles cancel and only the secular part survives: $\overset{ˉ}{R} \propto (3 sin^{2} i - 2) .$
Step 4 — Lagrange's planetary equation for the node. $\frac{d Ω}{d t} = \frac{1}{n a ^{2} 1 - e ^{2} s i n i} \frac{\partial R ˉ}{\partial i} .$ Since $\overset{ˉ}{R} \propto (3 sin^{2} i - 2)$ , $\partial \overset{ˉ}{R} / \partial i \propto 6 sin i cos i$ ; the $sin i$ in the numerator cancels the $sin i$ in the denominator, leaving the boxed $\dot{Ω} \propto - cos i$ .

The bulge, the sideways pull, and the tilted geometry are the whole story. Look at all four pictures once before answering — the reveals lean on them.

Picture 1 — the bulge and the off-center pull that becomes a torque.

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

Picture 2 — the celestial triangle behind $sin ϕ = sin i sin u$ .

Picture 3 — the $sin i$ (from the derivative) and the $cos i$ (what survives the cancellation), plotted together so you see the division.

Picture 4 — why $J_{3}$ averages to zero: an odd-latitude bulge feels equal-and-opposite pulls above and below the equator.

Spot the error

"Since J2 is a real extra force, it must add energy and speed the satellite up over time."

The force is real but the secular averaging (Step 3) gives

\overset{a}{˙} = \overset{e}{˙} = 0

, so mean orbital energy

- μ /2 a

is unchanged. In Picture 1 the sideways pull rotates

L

(tips it), it does not push the satellite faster along its path — J2 redistributes into orientation drift, not speed.

"Polar orbits cross the bulge most steeply, so they precess the fastest."

The steep crossing is true but irrelevant to node motion: precession about the polar axis needs

L

pointing along that axis, which is the

cos i

factor in Step 4. On Picture 3,

i = 90°

is precisely the zero crossing — precession stops there.

"For an eccentric orbit just use

a

in place of

p

; they're basically the same."

Only for circular orbits. In general

p = a (1 - e^{2})

, and the hidden

(1 - e^{2})^{2}

\dot{Ω} \propto p^{- 2}

can raise precession noticeably; swapping in

a

silently drops that factor.

"Sun-synchronous orbits need a prograde tilt to make

\dot{Ω}

positive."

Backwards —

\dot{Ω} \propto - cos i

, so a positive

\dot{Ω}

needs

cos i < 0

, i.e. a retrograde inclination like

i \approx 98°

(the right-hand shaded region of Picture 3). A prograde orbit gives

\dot{Ω} < 0

and can never sun-synchronize.

"The formula

sin ϕ = sin i sin u

only works for circular orbits."

It holds for any eccentricity — it is pure spherical trigonometry of the celestial triangle in Picture 2 (equator, orbit plane, satellite meridian), linking geocentric latitude

ϕ

to inclination

i

and argument of latitude

u = ω + ν

. The triangle knows nothing about the orbit's shape, only its angles.

"We can skip averaging over the orbit and just use the instantaneous

R = \frac{μ J _{2} R _{E}^{2}}{2 r ^{3}} (3 sin^{2} ϕ - 1)

Then you'd track meaningless short-period wiggles that cancel each revolution. Only the averaged

\overset{ˉ}{R} \propto (3 sin^{2} i - 2)

isolates the secular drift that accumulates orbit after orbit — that's the whole point of Step 3.

"The Legendre term

P_{2} (sin ϕ) = \frac{1}{2} (3 sin^{2} ϕ - 1)

is largest at the equator."

It's most negative at the equator (

ϕ = 0 \Rightarrow P_{2} = - \frac{1}{2}

) and largest (positive) at the poles (

ϕ = \pm 90° \Rightarrow P_{2} = + 1

); this is exactly the bracket sitting inside

R

in Step 1, and its sign flip across latitude is what produces the restoring pull toward the equator that Picture 1 draws.

It's what makes retrograde Sun-synchronous orbits possible — a wrong sign flips your designed inclination to the wrong side of Picture 3's zero crossing and the mission geometry fails; precession direction is a hard design constraint, not a cosmetic detail.

Why does the bulge exert a torque rather than just a stronger central pull?

Look at Picture 1: when the satellite is above or below the equatorial plane, the extra equatorial mass pulls it slightly back toward the equator (the pink arrow) — a component not aimed at Earth's center. An off-center force is exactly a torque, which rotates

L

Why does a top under gravity precess instead of falling over, and how is that analogous here?

Gravity torques the top's spin angular momentum sideways, so its axis circles rather than toppling; likewise the bulge torques the orbit's

L

sideways (Picture 1's arrow), so the plane swivels (precesses) instead of tumbling randomly.

Edge cases

What happens to

\dot{Ω}

exactly at

i = 90°

It is exactly zero:

cos 90° = 0

(Picture 3's axis crossing). A polar orbit's node line does not precess at all under J2 (higher-order terms aside).

What happens as inclination passes

90°

into retrograde (

i > 90°

cos i

turns negative, so

\dot{Ω} = - \frac{3}{2} n J_{2} (R_{E} / p)^{2} cos i

flips sign and the node drifts eastward — the right-hand shaded region of Picture 3; this regime is what Sun-synchronous designs exploit.

What is

\dot{Ω}

in the degenerate case

J_{2} = 0

(perfect sphere)?

Identically zero everywhere — no bulge, so

R = 0

in Step 1, no off-center pull, no torque, no precession, for any inclination or eccentricity.

At the critical inclination

i = 63.4°

, is the orbit fully "frozen"?

No — only the perigee is frozen (

\overset{ω}{˙} = 0

); the node still precesses because

cos 63.4° \approx 0.447 \neq = 0

, so

\dot{Ω} \neq = 0

e \to 1

(nearly parabolic,

a

fixed), what does the formula predict for

\dot{Ω}

It blows up, since

p = a (1 - e^{2}) \to 0

makes

(R_{E} / p)^{2} \to \infty

; physically the secular averaging of Step 3 breaks down for such extreme orbits, so the formula's validity, not just its value, is the real limit here.

In the limit of very large

a

(deep space), what happens to precession?

Both

n \to 0

(as

a^{- 3/2}

) and

(R_{E} / p)^{2} \to 0

, so

\dot{Ω} \to 0

; far from Earth the tiny bulge field is negligible and the orbit plane is effectively fixed.

At the equator (

i = 0°

, prograde circular), is precession maximal or minimal?

Maximal in magnitude:

cos 0° = 1

gives the largest

∣ \dot{Ω} ∣

(the leftmost, deepest point of Picture 3's

- cos i

), a westward drift — the "flat loop near the belly" case that feels the bulge asymmetry most.

Recall One-line self-check before you leave

Say why $\dot{Ω} \propto - cos i$ produces (a) max westward drift at $i = 0°$ , (b) zero at $i = 90°$ , (c) eastward drift for $i > 90°$ . If all three come instantly, you own the sign structure. Answer ::: (a) $cos 0 = 1 \Rightarrow$ large negative $\dot{Ω}$ (west); (b) $cos 90 = 0 \Rightarrow \dot{Ω} = 0$ ; (c) $cos i < 0 \Rightarrow \dot{Ω} > 0$ (east), enabling Sun-synchronous orbits.

Question bank — Orbital perturbations — J2 effect (oblateness), derivation of nodal precession

True or false — justify

Spot the error

Why questions

Edge cases