Kepler Orbital Elements — a, e, i, Ω, ω, ν Explained

Why Exactly Six Elements?

Newton’s law of gravitation is a second-order ordinary differential equation in three dimensions. Its general solution contains 2 × 3 = 6 constants of integration. Those six constants can be chosen in many ways; the Keplerian parameterisation chooses them to be geometrically meaningful: three describe the shape and size of the orbit, and three describe its orientation in space, plus the body’s position along it.

Equivalently, the instantaneous state of a body in orbit requires three position coordinates (r = [x, y, z]) and three velocity components (v = [vx, vy, vz]) at a specific epoch. These six numbers are called the state vector. Keplerian elements are simply a different, more intuitive basis for the same six-dimensional space.

6 Constants of integration for the 2-body gravitational ODE

3 Elements describing orbit shape & size (a, e, period)

3 Elements orienting the orbit in inertial space (i, Ω, ω)

1 Element giving position along the orbit at a given epoch (ν or M or t = t₀)

Two-body assumption: Keplerian elements are exact only in a pure two-body system (one point mass orbiting another). Real orbits deviate because of Earth’s oblateness, third-body gravity, atmospheric drag, and radiation pressure. The elements then become an osculating orbit — the best-fitting Keplerian ellipse at each instant, changing slowly over time.

The Six Keplerian Elements

The standard osculating orbital elements are defined below. The reference plane and direction vary with application: for Earth satellites the reference plane is the equatorial plane and the reference direction is the vernal equinox (γ); for heliocentric orbits both are referred to the ecliptic/J2000 frame.

a

Semi-major axis

Half the length of the longest diameter of the ellipse. Directly linked to the orbital period T and energy E via T² ∝ a³ (Kepler’s third law) and E = −GM/(2a).

Range: 0 < a < ∞ (km or AU)

e

Eccentricity

Shape of the conic section: e = 0 circle → 0 < e < 1 ellipse → e = 1 parabola (escape trajectory) → e > 1 hyperbola. The ratio c/a where c is the focus distance.

Range: 0 ≤ e < 1 (closed orbits)

i

Inclination

Angle between the orbital plane and the equatorial (or ecliptic) plane. i = 0 is prograde equatorial; i = 90° is polar; i > 90° is retrograde (orbits opposite Earth’s rotation).

Range: 0° ≤ i ≤ 180°

Ω

RAAN — Right Ascension of the Ascending Node

Angle from the vernal equinox γ to the ascending node (where the orbit crosses the equatorial plane going north), measured in the equatorial plane. Defines the rotation of the orbital plane about the polar axis.

Range: 0° ≤ Ω < 360°

ω

Argument of Periapsis

Angle from the ascending node to the periapsis (closest approach point), measured in the orbital plane. Orients the ellipse within the orbital plane. For equatorial orbits (i = 0), Ω + ω → longitude of periapsis ˜ω.

Range: 0° ≤ ω < 360°

ν

True Anomaly

Angle from periapsis to the satellite’s current position, measured at the focus. This is the only time-dependent element. At periapsis ν = 0; at apoapsis ν = 180°. Linked to time via Kepler’s equation.

Range: 0° ≤ ν < 360°

Periapsis and apoapsis distances

Once a and e are known, the closest and farthest points follow immediately:

rₚ = a (1 − e) ← periapsis (perigee / perihelion) rₛ = a (1 + e) ← apoapsis (apogee / aphelion) period T = 2π √(a³ / μ) ← Kepler III (μ = GM)

Naming conventions: The periapsis and apoapsis get different prefixes depending on the primary body — perigee/apogee (Earth), perihelion/aphelion (Sun), perijove/apojove (Jupiter), periapsis/apoapsis (generic).

Kepler’s Equation and Orbit Propagation

To find where a body is at a specific time t, we use a chain of three anomalies:

Mean anomaly M: grows linearly in time. M = n(t − t₀) where n = 2π/T is the mean motion and t₀ is the epoch of periapsis passage.
Eccentric anomaly E: a geometric angle defined on the auxiliary circle circumscribed about the ellipse.
True anomaly ν: the actual angle at the focus from periapsis to the body.

Solving Kepler’s equation

The link between mean and eccentric anomaly is Kepler’s transcendental equation:

M = E − e sin(E) ← Kepler’s equation

This cannot be inverted analytically. The standard numerical solution uses Newton–Raphson iteration starting from E₀ = M:

Eₙ₊₁ = Eₙ − (Eₙ − e sin Eₙ − M) / (1 − e cos Eₙ)

Convergence is typically achieved in 5–10 iterations for e < 0.9, and with Laguerre’s method for high-eccentricity orbits (comets, e → 1).

From eccentric to true anomaly

tan(ν/2) = √((1+e)/(1−e)) ⋅ tan(E/2) or equivalently: cos ν = (cos E − e) / (1 − e cos E) sin ν = √(1−e²) sin E / (1 − e cos E)

The orbital radius at any point follows from the conic section equation (vis-viva in angular form):

r = a (1 − e²) / (1 + e cos ν) ← orbit equation (focus at origin)

The vis-viva equation gives the speed at any radius: v² = μ(2/r − 1/a), independent of the eccentricity once a and r are known.

Converting to Position and Velocity (State Vector)

The full conversion from orbital elements to state vector ‘r, v’ in an inertial frame is a three-stage process:

Solve Kepler’s equation for eccentric anomaly E, then find true anomaly ν.
Find position and velocity in the perifocal frame (origin at focus, x-axis toward periapsis, z-axis normal to orbital plane):

r⃗ₘₑₙᴿ = [r cos ν, r sin ν, 0] ← r = a(1−e²)/(1+e cos ν) v⃗ₘₑₙᴿ = (μ/h) [−sin ν, e + cos ν, 0] ← h = √(μ a (1−e²))

Rotate to the inertial (IJK) frame using three successive rotations:

R = R₃(−Ω) ⋅ R₁(−i) ⋅ R₃(−ω) where R₃(θ) rotates about z, R₁(θ) about x r⃗ = R ⋅ r⃗ₘₑₙᴿ , v⃗ = R ⋅ v⃗ₘₑₙᴿ

The inverse problem — converting a state vector to orbital elements — requires computing the angular momentum h = r × v, the eccentricity vector e = v × h/μ − r̂, and then deriving each element by geometry. This is the routine performed by satellite tracking networks from radar measurements.

TLE — Two-Line Element Set: The most common format for distributing satellite orbital data. A TLE encodes a modified set of mean elements (using the SGP4/SDP4 propagation model which accounts for J2 and drag) as two 69-character lines. Over 27 000 objects are tracked this way by US Space Surveillance Network.

Real-World Orbit Examples

ISS a = 6 731 km · e ≈ 0.0001 · i = 51.6° · T ≈ 92 min

GPS a = 26 560 km · e ≈ 0 · i = 55° · T = 12 h (half-sidereal)

GEO a = 42 164 km · e ≈ 0 · i ≈ 0° · T = 24 h exactly

Moon a = 384 400 km · e = 0.055 · i = 5.1° · T = 27.3 days

Special orbit types

Sun-synchronous orbit (SSO): Inclination ~98° makes the J2-driven nodal precession exactly one revolution per year, so the orbit plane always faces the same Sun angle — ideal for remote sensing.
Molniya orbit: e ≈ 0.74, i ≈ 63.4°, T = 12 h. The 63.4° inclination freezes apsidal precession (ω ̇ = 0) so the apogee stays over the northern hemisphere for ≈8 hours per pass.
Tundra orbit: like Molniya but T = 24 h, giving a single high-latitude coverage window per day.
Frozen orbit: a combination of e and ω chosen so that J2 and J3 perturbations cancel, keeping the periapsis altitude constant. Mars Reconnaissance Orbiter uses a frozen orbit.

Halley’s Comet — an extreme ellipse

Halley’s Comet (1P/Halley) has a = 17.9 AU and e = 0.967, giving a perihelion of 0.59 AU (inside Venus’s orbit) and an aphelion of 35.1 AU (beyond Neptune). Its period is ~75 years, yet its periapsis velocity of ~54 km/s is much higher than Earth’s orbital speed of ~30 km/s — all because vis-viva demands ½mv² + Eₘₙₙ = constant.

Perturbations and J2

Real orbits deviate from pure Keplerian motion due to several effects. The most important for Earth satellites is J2, the oblateness term in Earth’s geopotential (Earth’s equatorial radius exceeds its polar radius by ~21 km).

J2 secular effects

Ω̇ = − (3/2) n J₂ (R⊕/p)² cos i ← nodal regression ω̇ = (3/4) n J₂ (R⊕/p)² (5 cos²i − 1) ← apsidal precession where p = a(1−e²), J₂ = 1.08263×10⁻³, R⊕ = 6371 km

Setting ω̇ = 0 gives the critical inclination 5 cos²i − 1 = 0, solved by i = 63.4° or 116.6°. This is why Soviet Molniya satellites all use 63.4°.

Other perturbations

Atmospheric drag: Reduces a and e secularly, eventually causing reentry. The ISS loses ~2 km altitude per month without reboost.
Third-body gravity: The Moon and Sun raise tides that perturb high-altitude orbits (GEO, HEO) on year-long timescales.
Solar radiation pressure (SRP): A 1 m² flat panel at 1 AU experiences ~4.6×10−6 N/m². Critical for large solar-sail spacecraft.
General relativistic precession: Mercury’s perihelion advances 43″/century beyond Newtonian predictions — one of the first confirmations of GR.

Osculating elements: When perturbations are present, we define orbital elements at each instant as those of the best-fitting Keplerian ellipse (tangent to the actual trajectory in position and velocity). These “osculating” elements change slowly over time and are the quantities tabulated in ephemerides.

Try the simulations: Orbital Mechanics lets you visualise precession and drag in real time; Asteroid Deflection shows how a small Δv changes all six elements simultaneously.