Orbital Period & Gravity Field Calculator: Kepler, g(r), v
Calculate orbital periods using Kepler's third law and evaluate gravity field quantities. Compare satellites around Earth, Moon, Mars, Jupiter, Sun, or custom celestial bodies.
Kepler III Derivation from Energy Conservation
The ISS completes an orbit in 92 minutes while GPS satellites take 12 hours—yet both follow the same simple law discovered by Kepler four centuries ago. The orbital period T = 2π√(a³/μ) relates time to geometry through a single equation. A common calculation mistake: using altitude instead of semi-major axis. Altitude is height above the surface; semi-major axis is distance from the center. Miss that distinction, and your period is off by 50% or more.
Deriving Kepler's Third Law
For a circular orbit, centripetal acceleration equals gravitational acceleration:
v²/r = GM/r² → v = √(μ/r)
Period T = 2πr/v = 2πr/√(μ/r) = 2π√(r³/μ)
For ellipses: replace r with semi-major axis a
T = 2π√(a³/μ)
This holds for any bound orbit—circular or elliptical. The remarkable fact: period depends only on the semi-major axis a and the gravitational parameter μ = GM. Eccentricity doesn't appear at all.
| Variable | Symbol | Units | Example (Earth) |
|---|---|---|---|
| Orbital period | T | seconds | 5,550 s (LEO) |
| Semi-major axis | a | meters or km | 6,771 km (LEO) |
| Gravitational parameter | μ | km³/s² | 398,600 km³/s² |
Units Warning: μ is often given in km³/s² while altitudes are in km. Make sure everything is consistent—either all km or all m. Mixing them gives wrong answers by factors of 10⁹.
Period vs Eccentricity: The Semi-Major Axis Rule
Here's what trips up students: a highly elliptical orbit and a circular orbit with the same semi-major axis have exactly the same period. Eccentricity changes the shape but not the time. The math proves this—eccentricity doesn't appear in T = 2π√(a³/μ) at all.
| Orbit Type | Eccentricity | a (km) | Periapsis | Apoapsis | Period |
|---|---|---|---|---|---|
| Circular | 0.00 | 10,000 | 10,000 km | 10,000 km | 2.97 hr |
| Slightly elliptical | 0.20 | 10,000 | 8,000 km | 12,000 km | 2.97 hr |
| Highly elliptical | 0.60 | 10,000 | 4,000 km | 16,000 km | 2.97 hr |
| Very eccentric | 0.90 | 10,000 | 1,000 km | 19,000 km | 2.97 hr |
What Eccentricity Does Affect
Eccentricity determines the speed variation. At periapsis (closest point), velocity is highest; at apoapsis (farthest), it's lowest. A satellite in an e = 0.9 orbit spends most of its time crawling near apoapsis, then whips around periapsis. But the total time for one complete orbit? Same as a circular orbit with equal semi-major axis.
Practical Example: Molniya orbits (used for Russian communications) are highly elliptical (e ≈ 0.74) with 12-hour periods—same as GPS circular orbits. Both have a ≈ 26,560 km.
Scenario Table: LEO to GEO to Lunar
Compare orbits from LEO to the Moon. The period scales as a^(3/2)—double the semi-major axis and the period increases by 2^1.5 ≈ 2.83×. The relationship is dramatically nonlinear.
| Orbit | Altitude (km) | a (km) | Period | vs LEO |
|---|---|---|---|---|
| LEO (ISS) | 400 | 6,771 | 92.5 min | Baseline |
| Medium (Starlink) | 550 | 6,921 | 95.6 min | +3.4% |
| GPS | 20,200 | 26,571 | 12 hr | +678% |
| GEO | 35,786 | 42,164 | 23.93 hr | +1,451% |
| Moon's orbit | 384,400 | 384,400 | 27.3 days | +42,400% |
The a^(3/2) Scaling
GEO altitude is about 6× LEO altitude, but the semi-major axis ratio is only ~6.2× (because you measure from Earth's center). The period ratio? 23.93 hr / 1.54 hr ≈ 15.5× . That matches (6.2)^1.5 ≈ 15.4. The a^(3/2) rule holds precisely.
Why GEO is Special: A 23.93-hour period matches Earth's sidereal rotation rate. The satellite appears stationary over one spot—ideal for TV, weather, and communication satellites that need a fixed ground footprint.
Mean Motion and Angular Velocity
Mean motion n = 2π/T is the average angular velocity—how fast the satellite sweeps through its orbit in radians per second. For circular orbits, this is the actual angular velocity at all times. For elliptical orbits, actual angular velocity varies, but mean motion gives the average rate.
| Orbit | Period | Mean Motion (rad/s) | Revs/day |
|---|---|---|---|
| LEO (400 km) | 92.5 min | 1.13×10⁻³ | 15.5 |
| GPS | 12 hr | 1.45×10⁻⁴ | 2.0 |
| GEO | 23.93 hr | 7.29×10⁻⁵ | 1.0 |
| Moon | 27.3 days | 2.66×10⁻⁶ | 0.037 |
Why Mean Motion Matters
Two-Line Elements (TLEs) used to track satellites list mean motion in revolutions per day, not period in hours. ISS has n ≈ 15.5 rev/day. GPS satellites have n = 2.0 rev/day (exactly, by design—two orbits per sidereal day means the ground track repeats). Converting: n (rad/s) = 2π / T (seconds).
Elliptical Orbits: In a Molniya orbit, the satellite crawls at 0.1 km/s near apoapsis but races at 10 km/s through periapsis. Mean motion averages this: n = 2π/(12 hours) regardless of how the velocity varies.
Worked Comparison: ISS vs GPS vs Moon
Let's calculate orbital periods for three objects orbiting Earth and see exactly how Kepler's third law plays out with real numbers.
The Setup
Using T = 2π√(a³/μ) with Earth's μ = 398,600 km³/s²:
ISS (LEO)
- h = 400 km
- a = 6371 + 400 = 6771 km
- a³/μ = (6771)³/398600
- T = 2π√(778.5) = 5545 s = 92.4 min
GPS
- h = 20,200 km
- a = 6371 + 20200 = 26571 km
- a³/μ = (26571)³/398600
- T = 2π√(4.71×10⁷) = 43077 s = 12.0 hr
Moon
- a = 384,400 km
- a³/μ = (384400)³/398600
- T = 2π√(1.43×10¹¹)
- T = 2.36×10⁶ s = 27.3 days
| Metric | ISS | GPS | Moon |
|---|---|---|---|
| Semi-major axis | 6,771 km | 26,571 km | 384,400 km |
| Period | 92.4 min | 12.0 hr | 27.3 days |
| a ratio vs ISS | 1.0× | 3.9× | 56.8× |
| T ratio vs ISS | 1.0× | 7.8× | 425× |
Check: GPS has a 3.9× larger semi-major axis, so T ratio should be 3.9^1.5 = 7.7×. Moon has 56.8× larger a, so T ratio should be 56.8^1.5 = 428×. Both match our calculations within rounding. Kepler's law works.
Altitude-Period Trade-Off Curves
Higher altitude means longer period—but the relationship isn't linear. Small altitude changes at LEO barely affect period, while the same delta-h at GEO has larger effects (in absolute terms, though smaller percentages).
| Altitude Range | Period Change | dT/dh | Sensitivity |
|---|---|---|---|
| 200 → 400 km | 88.5 → 92.4 min | ~1.2 s per km | Low (deep in gravity well) |
| 400 → 600 km | 92.4 → 96.6 min | ~1.3 s per km | Low |
| 20,000 → 20,200 km | 11.9 → 12.0 hr | ~2.7 s per km | Medium |
| 35,600 → 35,800 km | 23.8 → 24.0 hr | ~3.6 s per km | High |
Why This Matters for GEO Stations
GEO satellites must maintain period within seconds of the sidereal day (86,164 s) to stay "stationary." A 100 km altitude error at GEO changes the period by ~6 minutes—enough to drift the satellite across the sky within weeks. Station-keeping maneuvers correct for this, but initial insertion must be precise.
Design Tip: To find the altitude for a desired period, invert Kepler's law: a = (μ × (T/2π)²)^(1/3), then h = a − R. For T = 24 hours around Earth: a = 42,164 km, h = 35,793 km.
J2 Oblateness and Real-World Precession
The two-body Kepler formula T = 2π√(a³/μ) assumes a perfectly spherical central body. Earth isn't—it bulges at the equator (J2 oblateness). This causes orbits to precess, and actual periods differ slightly from the Kepler prediction.
| Effect | Cause | LEO Impact | GEO Impact |
|---|---|---|---|
| Nodal regression | J2 oblateness | ~5-7°/day | ~0.01°/day |
| Apsidal precession | J2 oblateness | ~3-4°/day | ~0.005°/day |
| Period deviation | J2 oblateness | ~1-2 s | ~0.01 s |
| Orbital decay | Atmospheric drag | Significant (<600 km) | Negligible |
Sun-Synchronous Orbits: Using J2
Clever mission designers exploit J2 precession. A sun-synchronous orbit has inclination and altitude chosen so nodal regression exactly matches Earth's motion around the Sun (~1°/day). The satellite crosses the equator at the same local time every orbit—ideal for Earth observation with consistent lighting.
Educational Note: This calculator ignores J2 effects. For LEO satellites where J2 matters significantly, professional software (STK, GMAT) includes full perturbation models. The Kepler period is typically accurate to ~0.1% for first-order estimates.
IAU Standard Gravitational Parameters
Orbital mechanics calculations use the gravitational parameter μ = GM rather than G and M separately. Why? μ is measured directly from orbital observations with far higher precision than either G or M individually.
| Body | μ (km³/s²) | Uncertainty | Source |
|---|---|---|---|
| Sun | 1.32712440042×10¹¹ | ±1×10¹ | IAU 2015 |
| Earth | 398,600.4418 | ±0.0008 | IAU 2015 |
| Moon | 4,902.8001 | ±0.0003 | DE440 |
| Mars | 42,828.375 | ±0.002 | DE440 |
| Jupiter | 1.26687×10⁸ | ±10 | DE440 |
Why μ Precision Matters
Earth's μ is known to 10 significant figures. The gravitational constant G is only known to about 5 figures. Using G × M to compute μ introduces unnecessary error. For GPS satellites maintaining centimeter-level positioning, this precision matters—their orbits are predicted using the best available μ, not computed G × M.
Best Practice: Always use tabulated μ values from JPL's DE440/441 ephemerides or IAU resolutions. Don't compute μ = G × M unless you have no other option. The precision difference can matter for navigation applications.
Sources & References
- Bate, R. R., Mueller, D. D., & White, J. E. (2020). Fundamentals of Astrodynamics (2nd ed.). Dover Publications. — The classic textbook for orbital mechanics, Kepler's laws, and perturbation theory.
- Curtis, H. D. (2019). Orbital Mechanics for Engineering Students (4th ed.). Butterworth-Heinemann. — Comprehensive coverage of orbital period calculations with J2 perturbation analysis.
- IAU Resolution B3 (2015) — iau.org — Standard values for gravitational parameters of major solar system bodies.
- NASA JPL Horizons & DE440 — ssd.jpl.nasa.gov — High-precision ephemeris data for solar system objects.
- Vallado, D. A. (2013). Fundamentals of Astrodynamics and Applications (4th ed.). Microcosm Press. — Industry reference for satellite orbital mechanics including perturbations.
Formulas implemented using IAU 2015 standard values. For real mission planning requiring J2, drag, and multi-body effects, use professional software like STK or GMAT.
Fixing Kepler III and Orbital Period Mistakes
Real questions from students confused about altitude vs radius, circular vs elliptical periods, and real-world perturbations.
I plugged in 400 km altitude for ISS but got a period of 0.5 minutes instead of 92 minutes—what did I do wrong?
You used altitude instead of semi-major axis. The formula T = 2π√(a³/μ) requires the distance from Earth's center, not from the surface. For ISS at 400 km altitude, add Earth's radius: a = 6371 + 400 = 6771 km = 6.771×10⁶ m. With μ = 3.986×10¹⁴ m³/s², you get T = 5553 s ≈ 92.6 minutes. Using 400 km directly gives T = 2π√((4×10⁵)³/(3.986×10¹⁴)) ≈ 31.7 s—off by a factor of 175. Always convert altitude to orbital radius first.
My Hohmann transfer calculation gives the same period for a circular and elliptical orbit—that can't be right, can it?
It is right! Kepler's third law depends only on semi-major axis, not eccentricity. A circular orbit at a = 10,000 km and an ellipse with periapsis 7000 km, apoapsis 13,000 km both have a = 10,000 km and identical periods. The ellipse moves faster near periapsis and slower near apoapsis, but these effects cancel exactly. This is why transfer orbits work: the ellipse connecting LEO to GEO has a semi-major axis of (6671 + 42164)/2 ≈ 24,418 km, and that alone determines its half-period of ~5.3 hours.
I calculated GEO altitude as 35,786 km but my professor marked me wrong for using 42,164 km—which is correct?
Both are correct for different quantities. The semi-major axis is a = 42,164 km (from Earth's center). The altitude is h = a − R_Earth = 42,164 − 6,378 = 35,786 km (above the equator). The formula T = 2π√(a³/μ) needs 42,164 km. If your professor asked for altitude and you wrote 42,164 km, that's the error. Always label your answer clearly: 'altitude h = 35,786 km' or 'orbital radius a = 42,164 km.' The distinction matters for mission planning—launch guidance needs altitude, orbit mechanics need radius.
I'm using G = 6.674×10⁻¹¹ and M_Earth = 5.972×10²⁴ but my period is off by 0.1%—is my mass value wrong?
Your method is wrong, not your values. G is only known to ~5 significant figures, and M_Earth has similar uncertainty. Multiplying them compounds errors. Use μ = GM directly: Earth's μ = 3.986004418×10¹⁴ m³/s² is known to 10 significant figures from satellite tracking. JPL publishes these values for mission-critical calculations. For homework, μ ≈ 3.986×10¹⁴ m³/s² is fine. For precision work, never compute G×M separately—always look up μ from authoritative sources like IAU or JPL HORIZONS.
The ISS orbit is listed as 51.6° inclination—why doesn't that affect the period calculation?
Inclination affects where the orbit goes (ground track), not how long it takes. The period T = 2π√(a³/μ) is identical for equatorial (0°), polar (90°), and any inclined orbit with the same semi-major axis. However, inclination does matter for real ISS tracking because Earth's equatorial bulge (J2 effect) causes the orbit to precess differently at different inclinations. A 51.6° orbit precesses westward at about 5°/day, which slightly affects the apparent period relative to ground stations—but the inertial period remains ~92.6 minutes.
I need to find when a satellite passes overhead, but the period formula just gives me one number—how do I get a schedule?
The period tells you how often, not when. For overhead passes, you need the full orbital elements (epoch, semi-major axis, eccentricity, inclination, RAAN, argument of periapsis, true anomaly) plus your ground station coordinates. From there, propagate the orbit forward in time and check when the satellite's ground track comes within your visibility cone. Software like GMAT, STK, or even N2YO.com does this automatically. The period is just one input—you also need the initial position (epoch) and orientation (the other five elements).
My textbook says GPS satellites orbit twice per day, but 12-hour period gives me ~20,200 km altitude—Google says 20,180 km. Why the 20 km discrepancy?
GPS actually uses a half-sidereal-day period (~11 hours 58 minutes), not exactly 12 hours. A sidereal day is 23h 56m 4s, so half is 43,082 seconds. This gives a = (μT²/(4π²))^(1/3) = 26,560 km from Earth's center, or altitude = 26,560 − 6,378 = 20,182 km. The ~20 km variations between sources come from using different Earth radius values (equatorial vs mean) and whether they account for J2 secular effects. For homework, 20,200 km is acceptable; for navigation, millimeter-level ephemerides matter.
I derived v = √(μ/a) for circular orbit velocity but my friend got v = 2πa/T—are these different formulas?
They're identical! Substitute T = 2π√(a³/μ) into v = 2πa/T: v = 2πa/(2π√(a³/μ)) = a/√(a³/μ) = √(μa/a³) = √(μ/a). Both give the same velocity. The vis-viva equation generalizes this to elliptical orbits: v = √(μ(2/r − 1/a)), which reduces to √(μ/a) when r = a (circular orbit). Use whichever form has the variables you know—if you have period, use 2πa/T; if you have μ and radius, use √(μ/r).
My simulation shows the Moon's period as 27.3 days but my formula gives 27.5 days—what accounts for the 0.2-day error?
Several factors: (1) The Moon's orbit is elliptical (e ≈ 0.055), so the semi-major axis you used might be the current distance, not the mean. Use a = 384,400 km for the semi-major axis. (2) Earth-Moon μ should be the sum: μ_system = μ_Earth + μ_Moon ≈ 4.035×10¹⁴ m³/s², not just Earth's μ—though this correction is only ~1.2%. (3) The 27.3-day 'sidereal month' is measured relative to stars; if your simulation uses synodic months (~29.5 days), that's a different quantity. With a = 384,400 km and μ = 4.035×10¹⁴ m³/s², you get T ≈ 27.28 days.
Why do low Earth orbit satellites need reboosting every few months if orbital mechanics is supposed to be predictable?
Atmospheric drag. Even at 400 km, there's enough residual atmosphere (~10⁻¹² kg/m³) to slow satellites by meters per second over weeks. This lowers the orbit, which increases drag further (exponential decay). ISS loses about 2 km/month in altitude and needs periodic reboosts. The Kepler/two-body formula T = 2π√(a³/μ) assumes vacuum—it's exact for the Moon but an approximation for LEO. Below ~300 km, decay is rapid (days to weeks); above ~1000 km, orbits last centuries. Drag depends on solar activity, which heats and expands the atmosphere unpredictably.