1 Why the two-body ellipse is not quite enough

GPSORB01–GPSORB09 developed the ideal Kepler problem from

(1)

Because the force is central,

(2)

so the orbital plane is fixed. The eccentricity vector is also constant, so periapsis is fixed within that plane.

Real Earth violates the spherical assumption. Its equatorial radius is larger than its polar radius, and the external gravity field is therefore not exactly proportional to 1∕r. The largest correction is the degree-two zonal harmonic J₂.

The immediate consequences are important:

1.

the acceleration is no longer exactly parallel to −r;

2.

the specific angular momentum vector can slowly change direction;

3.

the ascending node can regress or advance;

4.

the periapsis direction can precess within the orbital plane.

These effects are small during one revolution, but they accumulate over many orbits. That makes them especially important for a long-lived navigation Constellation.

2 External gravity as a multipole expansion

Outside Earth’s mass distribution the gravitational potential satisfies Laplace’s equation,

(3)

For an axisymmetric body, there is no longitude dependence. Separation of variables in spherical coordinates gives an exterior solution of the form

(4)

Here 𝜃 is geocentric colatitude, R_E is a reference equatorial radius, P_ℓ is a Legendre polynomial, and J_ℓ is a dimensionless zonal coefficient.

The monopole term ℓ = 0 is the familiar spherical potential −μ∕r. The dipole term ℓ = 1 vanishes when the origin is placed at Earth’s center of mass. The leading departure from sphericity is therefore the quadrupole term ℓ = 2.

Using geocentric latitude ϕ = π∕2 − 𝜃 gives

(5)

and

(6)

Keeping only J₂ in Equation (4),

(7)

Equivalently,

(8)

The second term is the J₂ perturbing potential per unit mass.

For Earth, J₂ is positive and approximately

(9)

The small number is deceptive: a part-per-thousand distortion of the dominant gravity field is easily large enough to rotate orbital orientation by many degrees over a year.

3 From potential to the J2 acceleration

The acceleration follows from

(10)

Writing

(11)

Equation (8) can be differentiated in Cartesian coordinates. The total acceleration through degree two is

(12)

(13)

(14)

The correction is not generally parallel to r. Therefore

(15)

is generally nonzero. This is the physical origin of orbital-plane precession.

4 Geometry inside the orbital plane

Let

(16)

be the argument of latitude, where ω is the argument of periapsis and ν is true anomaly. For an orbit inclined by i to the equator, the geocentric latitude satisfies

(17)

This relation is the key bridge between the Earth-fixed symmetry of J₂ and the satellite’s orbital elements.

Substituting Equation (17) into the perturbing potential gives

(18)

It is convenient in the Lagrange planetary equations to define the disturbing function as the negative perturbing potential,

(19)

Hence

(20)

Equation (20) contains both short-period terms and terms that survive orbital averaging. The secular changes come from the orbit-averaged disturbing function.

5 Averaging the J2 disturbing function

For a Keplerian reference ellipse,

(21)

The relation between mean anomaly and true anomaly is

(22)

An average over one orbital period is therefore

(23)

First evaluate ⟨r⁻³⟩. Using the equations above,

(24)

Integrating from 0 to 2π removes the cosine term, giving

(25)

Next,

(26)

because the average of sin ²(ω + ν) is 1∕2, while the additional term proportional to e cos ν integrates to zero over a full revolution.

Applying Equations (25) and (26) to Equation (20),

(27)

This single averaged scalar function contains the first-order secular J₂ rates.

6 Lagrange planetary equation for the ascending node

For an averaged disturbing function that is independent of the node itself, the Lagrange planetary equation for Ω is

(28)

Differentiate Equation (27):

(29)

Substitute into Equation (28):

(30)

Using

(31)

and

(32)

we obtain the standard secular nodal rate

(33)

For a prograde orbit with 0 < i < 90^∘, cos i > 0, so

(34)

The ascending node regresses.

For a polar orbit, i = 90^∘ and the first-order secular node rate vanishes. For a retrograde orbit, cos i < 0 and the sign reverses.

7 Lagrange planetary equation for periapsis

The argument of periapsis is affected both by the eccentricity dependence of the disturbing function and by the moving node. With the present disturbing-function convention,

(35)

Differentiate Equation (27) with respect to eccentricity. Since

(36)

we obtain

(37)

The first term in Equation (35) becomes

(38)

The node-coupling term is

(39)

Adding the two pieces gives

(40)

This is apsidal precession: the eccentricity vector rotates inside the orbital plane.

8 The critical inclination

Equation (40) vanishes when

(41)

Therefore

(42)

and

(43)

These are the classical critical inclinations. At this first-order J₂ level, the argument of periapsis has no secular drift there.

GPS uses an inclination near 55^∘, not the critical inclination, so both nodal and apsidal secular effects are present.

9 What remains constant in the first-order secular J₂ model

The averaged J₂ model has an important structure. To first order in the secular approximation,

(44)

But

(45)

Thus the orbit keeps approximately the same size, shape, and inclination while its orientation changes.

This does not mean the instantaneous elements are perfectly constant. The unaveraged J₂ acceleration produces short-period variations in several elements. The statements above refer only to the leading first-order secular behavior.

There is also a useful symmetry interpretation. Because the J₂ potential is axisymmetric about Earth’s spin axis, the potential does not depend on inertial longitude. The component of angular momentum along the symmetry axis is therefore conserved. The full angular-momentum vector can precess while maintaining its projection on z.

10 Inclination dependence of the two secular rates

Both secular rates scale with

(46)

but their inclination dependence differs:

(47)

(48)

The node rate changes sign at 90^∘. The apsidal rate changes sign at the two critical inclinations.

11 Numerical example for a GPS-like orbit

Use representative values

(49)

with

(50)

and

(51)

The semilatus rectum is

(52)

The Keplerian mean motion is

(53)

corresponding to an orbital period of approximately

(54)

Substituting into Equation (33),

(55)

Over one year, the pure first-order J₂ contribution would accumulate to roughly

(56)

For periapsis,

(57)

or about

(58)

Across a two-hour interval, however, these changes are only about

(59)

(60)

The rates are small over one broadcast-ephemeris interval but large over the lifetime of a satellite.

12 Connection to GPS constellation behavior

The GPS space segment uses medium-Earth orbits with nominal inclination near 55^∘. The GPS Standard Positioning Service Performance Standard lists a 55 degree reference inclination, and official constellation descriptions place the satellites near 20,200 km altitude.

The J₂ nodal-regression estimate above is therefore physically representative of GPS. It explains why an inertial ascending node cannot be treated as a permanent constant over months or years.

It is important, however, to distinguish three different quantities:

1.

the ideal two-body constant Ω;

2.

the analytic first-order perturbation rate _J2;

3.

the broadcast ephemeris coefficient transmitted by GPS.

The third is not merely Equation (33) copied into the navigation message. The broadcast ephemeris parameters are fitted so that the user equations reproduce a high-fidelity propagated trajectory over a limited fit interval. The fitted rate therefore represents the behavior required by that broadcast model, including more physics than the single analytic J₂ term.

13 Why the legacy broadcast model contains a node rate but not an explicit periapsis rate

The legacy GPS ephemeris includes parameters such as

(61)

as well as six harmonic correction coefficients. The user algorithm explicitly propagates the node longitude with a rate and propagates inclination with IDOT.

The legacy model does not transmit a separate term. This does not imply that real periapsis is physically stationary. Rather, ω is part of a short-arc fitted ephemeris, and residual orbital behavior is represented by the fitted element set and harmonic corrections over the applicable interval.

This distinction matters for learning the physics: Equation (40) tells us why periapsis should precess in a real oblate gravity field, even though the legacy broadcast parameterization chooses not to expose that physical rate as a separate user parameter.

14 Why J2 also hints at the harmonic correction structure

Equation (20) contains

(62)

Using

(63)

we immediately see why degree-two gravity naturally produces second-harmonic structure involving

(64)

and, after the full perturbation solution is developed,

(65)

This is the physical clue behind the GPS correction form

(66)

with analogous radius and inclination corrections.

The broadcast coefficients should not be identified one-for-one with a single closed-form J₂ perturbation coefficient. They are fitted parameters. But the double-angle structure is exactly what one expects when the leading nonspherical gravity term is degree two.

15 Beyond J2

A realistic GPS orbit model contains additional perturbations,

(67)

Higher geopotential harmonics, third-body gravitation, solar radiation pressure, Earth radiation effects, and maneuver history all contribute to the precise trajectory.

The purpose of J₂ is not to replace that high-fidelity model. Its value is that it is the first perturbation for which the ideal Kepler picture visibly breaks in a systematic, derivable way. It explains two of the most important qualitative facts about real satellite orbits:

(68)

and

(69)

16 Summary

The main derivation chain is

(70)

followed by

(71)

The two key secular rates are

(72)

and

(73)

For a GPS-like orbit near 55^∘ inclination, the first is negative and the second is positive. These rates supply the physical bridge from the fixed Keplerian orientation of GPSORB05–GPSORB09 to the slowly changing orbital geometry represented by real GPS ephemerides.

References