Kepler's Laws of Planetary Motion

Kepler's laws of planetary motion are three fundamental kinematic principles of celestial mechanics formulated by German mathematician and astronomer Johannes Kepler between 1609 and 1619 that govern the orbital motion of planets, moons, spacecraft, binary stars, and exoplanets under a central inverse-square gravitational force. Derived empirically through years of geometric analysis of Danish nobleman Tycho Brahe's pre-telescopic, arcminute-precision positional observations of Mars (~1–2′ accuracy), the first two laws—the Law of Ellipses and the Law of Equal Areas—were published in Astronomia Nova (New Astronomy, 1609), while the third law—the Harmonic Law—was published a decade later in Harmonices Mundi (The Harmony of the World, 1619).

By replacing the two-millennia-old Aristotelian and Ptolemaic dogma of uniform circular motion (and the complex epicycles retained even in Nicolaus Copernicus's 1543 heliocentric model) with a physically grounded heliocentric ellipse where orbital speed varies continuously with heliocentric distance, Kepler transformed astronomy from a purely geometric bookkeeping exercise into physical astrodynamics. Sixty-eight years later, in his Philosophiae Naturalis Principia Mathematica (1687), Isaac Newton demonstrated that all three of Kepler's empirical laws are direct mathematical consequences of his three laws of motion combined with the universal law of gravitation (F = G M₁ M₂ / r²), generalizing Kepler's laws to arbitrary conic-section trajectories (circles, ellipses, parabolas, and hyperbolas) and providing the dynamical foundation for the vis-viva equation, interplanetary Hohmann transfer orbits, and modern gravitational mass measurements across the universe.

---

First Law: The Law of Ellipses and Conic Sections

Kepler's First Law states that the orbit of every planet is an ellipse with the Sun (or more precisely, the two-body system's center of mass, or barycenter) located at one of the two foci. An ellipse is the locus of points in a plane for which the sum of the distances to two fixed focal points (F₁ and F₂, separated by 2c) equals a constant 2a, where a is the semi-major axis (half the longest diameter) and b is the semi-minor axis (half the shortest diameter). The shape of the ellipse is parameterized by its dimensionless orbital eccentricity (e):

e = c / a = √(1 − b² / a²) (where 0 ≤ e < 1 for a bound ellipse)

In plane polar coordinates (r, θ) centered on the primary focus (F₁, occupied by mass M₁) with the true anomaly θ measured from the point of closest approach (periapsis or perihelion), the equation of a conic section derived from the Newtonian radial equation of motion (d²u/dθ² + u = G(M₁ + M₂) / h², where u = 1/r and h = r² dθ/dt is the specific angular momentum) is:

r(θ) = p / (1 + e cos θ) = a(1 − e²) / (1 + e cos θ)

where p = a(1 − e²) = b² / a = h² / [G(M₁ + M₂)] is the semi-latus rectum. Evaluating r(θ) at θ = 0 and θ = π gives the minimum and maximum orbital radii:

• Periapsis (Perihelion) Distance (θ = 0): r_p = a(1 − e)
• Apoapsis (Aphelion) Distance (θ = π): r_a = a(1 + e)
• Geometric Mean Distance: b = √(r_p r_a) = a √(1 − e²)

In Newtonian mechanics, bound ellipses (0 ≤ e < 1, total specific mechanical energy ε = −GM / 2a < 0, with a circle as the degenerate case e = 0) are joined by unbound open trajectories that equal or exceed the local escape velocity: parabolic escape trajectories (e = 1, ε = 0) and hyperbolic flyby trajectories (e > 1, ε > 0).

Second Law: The Law of Equal Areas and Angular Momentum Conservation

Kepler's Second Law states that a line segment (radius vector r) joining a planet and the Sun sweeps out equal areas during equal intervals of time. That is, the areal velocity (dA/dt) is constant throughout the orbit. Over an infinitesimal time dt, the radius vector r rotates through an angle dθ, sweeping out an infinitesimal triangular sector of area dA = ½ r (r dθ) = ½ r² dθ. Dividing by dt links the areal velocity directly to the specific orbital angular momentum h = |r × v| = r v_⊥ = r² (dθ/dt):

dA / dt = ½ r² (dθ / dt) = ½ r v sin φ = L / (2 μ) = ½ h = constant

where μ = M₁ M₂ / (M₁ + M₂) is the two-body reduced mass, L = μ h is the total orbital angular momentum, and φ is the angle between the position vector r and the velocity vector v. Because the gravitational force F = −(G M₁ M₂ / r²) r̂ acts strictly along the radial line connecting the two bodies (central force), the net torque about the focus vanishes identically (τ = r × F = 0). Therefore, dL/dt = 0, proving that Kepler's Second Law holds for any central force law F(r) r̂, independent of whether the radial dependence is inverse-square.

At perihelion and aphelion (φ = 90°, sin φ = 1), conservation of specific angular momentum requires r_p v_p = r_a v_a = h. Combining this with conservation of specific mechanical energy via the vis-viva equation (v² = GM (2/r − 1/a)), where v is the orbital velocity determining the body's orbital translational kinetic energy per unit mass (½ v²), yields the maximum and minimum orbital speeds:

v_p = √[ (GM / a) (1 + e) / (1 − e) ] and v_a = √[ (GM / a) (1 − e) / (1 + e) ]

Consequently, v_p / v_a = (1 + e) / (1 − e): Earth (e = 0.0167) travels at 30.29 km s⁻¹ at perihelion in January versus 29.29 km s⁻¹ at aphelion in July, whereas Halley's Comet (e = 0.967) accelerates from 0.91 km s⁻¹ beyond Neptune's orbit at aphelion (35.1 AU) to 54.6 km s⁻¹ at perihelion (0.586 AU).

---

Third Law: Period–Semimajor Axis Relation and Mass Determination

Kepler's Third Law (the Harmonic Law) states that the square of the sidereal orbital period (T) of a planet is directly proportional to the cube of the semi-major axis (a) of its orbit: T² ∝ a³. To derive the constant of proportionality in Newtonian mechanics, one integrates the constant areal velocity dA/dt = h / 2 over one complete orbital period T, setting the total swept area equal to the geometric area of an ellipse (A = π a b = π a² √(1 − e²)):

A = (h / 2) T = π a² √(1 − e²)

Squaring both sides and substituting the semi-latus rectum relation h² = G (M₁ + M₂) a (1 − e²) cancels out the eccentricity factor (1 − e²) identically, yielding Newton's generalized form of Kepler's Third Law for the two-body problem:

T² = [ 4π² / G (M₁ + M₂) ] a³

Notice that for an elliptical orbit of a given semi-major axis a, the orbital period T depends only on a and the total system mass (M₁ + M₂), and is completely independent of the orbital eccentricity e. When a is expressed in astronomical units (1 AU = 1.495978707 × 10¹¹ m), T in Julian years (1 yr = 3.15576 × 10⁷ s), and masses in solar masses (M_☉), the gravitational constant factor 4π² / (G M_☉) equals 1 AU³ yr⁻², simplifying the relation to:

(M₁ + M₂) / M_☉ = (a / AU)³ / (T / yr)²

Because every planet in the Solar System has a mass M₂ ≪ M_☉ (even Jupiter has M_Jup / M_☉ = 9.546 × 10⁻⁴), the ratio a³ / T² equals 1.000 to within 0.1% for all planets, comets, and asteroids orbiting the Sun. More importantly, Newton's inclusion of the mass sum (M₁ + M₂) transforms Kepler's Third Law into the primary scale for weighing objects across the cosmos: by measuring the semi-major axis a and orbital period T of a natural satellite (such as Io orbiting Jupiter or Titan orbiting Saturn), a visual or eclipsing binary star system, a transiting exoplanet, or the star S2 orbiting Sagittarius A\* (a = 1,018 AU, T = 16.05 yr → M_BH = 4.15 × 10⁶ M_☉), astronomers directly determine the gravitating mass of the central body.

Kepler's Equation and Perturbations From Pure Two-Body Motion

Because radial distance r and angular speed dθ/dt vary nonlinearly with time along an eccentric orbit (e > 0), predicting a planet's position θ(t) at an arbitrary time t after perihelion passage t_0 cannot be written in closed algebraic form. Instead, Kepler introduced an auxiliary geometric angle, the eccentric anomaly (E), projected onto a circumscribed circle of radius a, and linked it to the uniform mean anomaly M = n (t − t_0) = (2π / T) (t − t_0) via the transcendental Kepler's equation:

M = E − e sin E where tan(θ / 2) = √[ (1 + e) / (1 − e) ] tan(E / 2)

In modern astrodynamics, Kepler's equation is solved iteratively to machine precision in a few steps using Newton–Raphson root finding (E_{k+1} = E_k − (E_k − e sin E_k − M) / (1 − e cos E_k)). Finally, while Kepler's laws are exact for an isolated two-body system of spherical point masses in Newtonian gravity, real multi-planet systems experience small secular and resonant perturbations: mutual gravitational tugs among planets cause apsidal and nodal precession and lock bodies into mean-motion orbital resonances (such as the 4:2:1 Laplace resonance of Jupiter's moons Io, Europa, and Ganymede), while general relativistic spacetime curvature adds an extra perihelion precession of Δω = 6π G M / [a c² (1 − e²)] per orbit (42.98 arcseconds per century for Mercury).

---

Why does a planet move faster at perihelion than at aphelion?

A planet moves faster at perihelion because of conservation of angular momentum (L = m r v_⊥ = constant), which is the physical foundation of Kepler's Second Law. Because the Sun's gravitational pull points directly along the line connecting the Sun and the planet, it exerts zero torque (τ = r × F = 0). As the planet moves closer to the Sun (smaller r), its tangential orbital speed v must increase inversely (r_p v_p = r_a v_a), converting gravitational potential energy into orbital translational kinetic energy.

Does the orbital eccentricity of a planet change its orbital period?

No. According to Kepler's Third Law (T² = [4π² / G(M₁ + M₂)] a³) and the vis-viva equation (ε = −GM / 2a), the orbital period T and total specific mechanical energy ε of a bound orbit depend only on the semi-major axis a, not on the eccentricity e. A planet in a circular orbit (e = 0) with a = 1 AU and a comet in a highly elongated ellipse (e = 0.99) with the same semi-major axis a = 1 AU both take 1.00 year to complete one orbit.

Why did Kepler need Mars in particular to discover that planetary orbits are ellipses?

Among the naked-eye planets known in antiquity, Venus (e = 0.0068) and Earth (e = 0.0167) have nearly circular orbits whose deviation from a circle (b/a = √(1 − e²) ≈ 0.9999) was smaller than pre-telescopic measurement uncertainties. Mars, by contrast, has a relatively large orbital eccentricity of e = 0.0934, which displaced its apparent sky position by 8 arcminutes from Kepler's best circular vicarious hypothesis—four times larger than Tycho Brahe's ~2-arcminute observational error, compelling Kepler to abandon circular motion in favor of an ellipse.

---

1. Kepler, J. (1609). Astronomia Nova ΑΙΤΙΟΛΟΓΗΤΟΣ, seu Physica Coelestis, tradita commentariis de motibus stellae Martis ex observationibus G. V. Tychonis Brahe. (Translated by W. H. Donahue, 1992, New Astronomy, Cambridge University Press). https://doi.org/10.1017/CBO9780511608360
2. Wilson, C. (1968). "Kepler's derivation of the elliptical path." Isis, 59(1), 4–25. https://doi.org/10.1086/350331
3. Murray, C. D., & Dermott, S. F. (1999). Solar System Dynamics. Cambridge University Press. https://doi.org/10.1017/CBO9781139174817
4. Bate, R. R., Mueller, D. D., White, J. E., & Saylor, W. W. (2020). Fundamentals of Astrodynamics (2nd ed.). Dover Publications. https://store.doverpublications.com/0486497046.html