Concept

Orbital elements — where it appears

The set of numbers, usually six, that fixes an orbit's size, shape, orientation and the body's position along it. Two of the six are undefined for a circular orbit and one for an uninclined one, which are coordinate singularities in perfectly well-behaved motion.

Named by 12 essays across 2 fields — each of them below, with the objects they name alongside it.

An orbit at eccentricity 0.6. An orbit of eccentricity 0.6. The primary sits at a focus, offset from the centre by 0.6 of the semi-major axis, and the closest and furthest points differ by a factor of 4.00.

The orbit is an ellipse, and the Sun is not in the middle of it

Kepler's first law is usually drawn wrong. The interesting content is not the ellipse — it is the focus, and the fact that one of the two is empty.

orbits · The ellipse

Equal areas in equal times, which is angular momentum in disguise

Kepler's second law is a statement about the area a radius line sweeps. It looks like an odd thing to have noticed, and it turns out to be a conservation law arriving eighty years early.

orbits · Angular momentum

A vector that does not change, drawn five times. The Laplace–Runge–Lenz vector A = v × L − GM r̂, constructed at five points of one orbit at e = 0.45 under a force ∝ r^−2, each drawn from the body rather than from the focus so that its length and direction can be compared point by point. Under the inverse square every one of them is the same vector: identical to machine precision in length and in direction. Its length is 0.4500, which is the eccentricity, and it points at pericentre — so the orbit's orientation is a conserved quantity and not a constant of integration, which is why the ellipse does not turn.

The third thing that is conserved

Energy fixes an orbit's size and angular momentum fixes its shape. Neither fixes which way it points — so a curve that closes needs a third conserved quantity, and only two force laws in the universe supply one.

orbits · Conic sections

The three anomalies at E = 1.15 rad. The auxiliary circle construction. The eccentric anomaly E is measured at the centre, the true anomaly ν at the focus, and the mean anomaly M is time expressed as an angle. Here E = 1.150, ν = 1.827 and M = 0.602 radians, related by M = E − e sin E.

The position that has no formula, and is computed anyway

Kepler's equation relates where a body is to when it is there. It cannot be solved in elementary functions, Kepler said so, and nobody has managed it since — which has stopped nothing.

orbits · The ellipse

An orbit at i = 42°, Ω = 35°, ω = 55°. The three orientation elements. The orbit is generated in its own plane and rotated by the standard sequence, so the inclination, the node line and the argument of periapsis are the angles that produced the drawn curve rather than labels applied to it.

Six numbers that fix an orbit for all time, and the sixth is the awkward one

Five of the orbital elements describe a curve that never changes. The sixth says where on it the body is, and it is the only one that has to keep being measured.

orbits · Orbital elements

Three rotations, applied in order. An orbit of eccentricity 0.45 carried into space by the three orientation angles, one panel per rotation: the argument of periapsis ω = 40°, then the inclination i = 42°, then the longitude of the ascending node Ω = 55°. The last panel applies the same three angles in the reverse order and arrives somewhere else, because rotations about different axes do not commute.

Three rotations that put an orbit in space, and they do not commute

An orbit's orientation takes three angles. Give the same three angles in a different order and the orbit ends up somewhere else — which is why the convention is part of the data.

orbits · Orbital elements

The osculating semi-major axis of a perturbed orbit. The semi-major axis a test particle would have if the perturber vanished, computed from its position and velocity at every step of an integration over 26 orbits of the perturber. It is not constant: a short-period ripple rides on a slow trend, and only the trend accumulates.

Elements that do not stay constant

Six numbers fix an orbit for all time, and the phrase is only true in a universe containing two bodies. Add a third and the six start moving — some of them wandering and returning, one or two of them drifting in one direction forever, and the difference between those two behaviours is the whole of celestial mechanics after Newton.

orbits · Perturbations

Ceres from five directions and no distance. Ceres seen five times over 41 days, from an Earth on a circular orbit, reduced in the plane. Each sighting gives a direction and no range, so the object is somewhere on its sight line; the five lines here span 1.37° of geocentric arc altogether, and Earth's own motion supplies the only baseline there is — 0.403 AU of its 0.691 AU of travel lies across the sight lines. A planar orbit is four numbers, so five angles over-determine it and one orbit comes out: a = 2.7658 AU, e = 0.0785. That is the answer and not an input — the sightings were generated at a = 2.7658 AU and e = 0.0785, and the solve, which sees only the directions and the dates, returns them to 3e-12. The two shaded sectors are what closes the determination: between the first and middle sightings the radius vector sweeps 0.2993 AU² in 21.0 days and between the middle and last 0.2853 AU² in 20.0 days, a ratio of 1.04918 against a time ratio of 1.04918. Slide all three crossings out along their sight lines together and that equality fails at once, so it fixes the distance by itself, with no propagation anywhere in the argument — and it gives a = 2.7658 AU over again. The two dashed curves are candidates that thread the same three sight lines at 85% and 108% of the recovered distance: a = 1.86 AU at e = 0.35, sweeping its areas in a ratio 4.5% wrong; and a = 5.61 AU at e = 0.49, sweeping its areas in a ratio 2.9% wrong. A few per cent in the distance is an orbit of another kind, which is the same fact the conditioning panel measures: one arcsecond of angle error moves a by 0.60% on this arc.

Five directions and no distance among them

An image of a moving point records an angle and throws the range away, so an orbit has to be assembled out of angles alone. How many angles are needed is not a detail of the method — it is the whole of what a determination is.

orbits · Orbit determination

One admissible root, 0.01% from the truth. Gauss's reduction of three directions to a distance, drawn as the two relations whose intersection it is. Three observations of Ceres on days 0, 20, 40 of an arc, generated from its elements and used only as sight directions — no range, no radial velocity. The rising curve is geometry: the heliocentric distance a candidate at geocentric distance ρ₂ would have, r₂² = ρ₂² + 2ρ₂(R₂·L̂₂) + R₂², which contains no dynamics at all. The falling curve is dynamics: ρ₂ = A + µB/r₂³, with A and B built from the three sight vectors, the three observer positions and the three times, and containing no orbit. Eliminating ρ₂ between them gives r₂⁸ + a r₂⁶ + b r₂³ + c = 0 — an eighth-degree equation, from a problem with exactly as many equations as unknowns. Here they cross once at a positive ρ₂, at r₂ = 2.5893 AU against the true 2.5890. The other 2 real roots are rejected not by fitting but by sign: the ρ₂ each implies is negative, and an object behind the observer was not the thing observed.

Three observations and no orbit at all

Three directions in space give six numbers for the six elements of an orbit, which sounds like a solved problem. The algebra that solves it is of the eighth degree, and for a near-Earth asteroid three perfect observations can be consistent with three different orbits.

orbits · Orbit determination

A separation of 10⁻¹⁰ becomes 0.0285. The distance between two copies of Burrau's problem, started 10⁻¹⁰ apart in one coordinate of one body and then stepped in lockstep — one loop, one step size, both states advanced by it — for 60 time units. The vertical axis is logarithmic, so the straight stretch is exponential growth, and its slope is the Lyapunov exponent: 0.3413 per time unit, fitted over the 604 samples lying between ten times the starting offset and a tenth of the system's own size. That is an e-folding every 2.93 time units, so the separation multiplies by ten every 6.75 and by 7.8·10⁸ across the whole run. It has not saturated within the drawn interval: the growth rate over the last tenth of the run is still 0.6002 per time unit, 176% of the fitted exponent, and the separation has reached only 0.90% of the system's own size. The exponent is what sets a prediction horizon, and this run is drawn short of it on purpose: the straight line is the measurement, and the flat part that follows is arithmetic about how far apart two bounded systems can get.

A prediction with an expiry date

The inner solar system's Lyapunov time is about five million years, so a centimetre of error becomes an orbit in a hundred million. The ephemeris dies while the system survives — because the elements stay bounded when the phase does not, and only one of those is what stability means.

gravitation · Chaos $A thousandth of the field, and all of the precession. Left: the Earth's figure against a sphere of the same equatorial radius, with the flattening drawn 28× its true value. The real difference between the equatorial and polar radii is 21.4 km on 6378 — 1 part in 298 — which at this size would be 0.5 pixels and invisible, so the drawing is a schematic and the number is here instead. Right: the two components of the J₂ perturbation at the surface, each as a fraction of the monopole μ/r², both differentiated from the potential rather than quoted. The radial one strengthens the inward pull by 1.62×10⁻³ over the equator, where the extra mass is, and weakens it by 3.25×10⁻³ over the poles, vanishing at ±35.26° where P₂ does. The transverse one is zero at the equator and at both poles and peaks at ±45°, at 1.62×10⁻³ — and that is the component that does the work. It pulls an inclined orbit back towards the equatorial plane, which is a torque about the line of nodes, and a torque applied to something already turning moves it sideways rather than back. Averaged over an orbit the pair leave a, e and i untouched and turn the whole plane instead, which is why a term a thousandth of the field is the largest single perturbation on almost every satellite ever flown.$

The Earth's shape, read off a satellite's node

The Earth is a thousandth of a part from being a sphere, and that thousandth turns every satellite's orbital plane. Vanguard 1 measured it in 1959 — and one retrograde inclination turns the plane at exactly the rate the Sun moves, which is a perturbation used as a design constraint rather than corrected for.

gravitation · Oblateness

An invariant that is only almost one

The Tisserand parameter survives a close encounter with Jupiter exactly, and comet families are separated by boundaries in it drawn to two decimal places. The exactness holds for a Jupiter on a circle. Jupiter's eccentricity is 0.0489, and integrating the same encounter twice shows what that costs.

orbits · Tisserand parameter

Named alongside it

The objects these essays reach for when they reach for this one.

Eccentricity Periapsis Anomaly Epoch Kepler's equation Argument of periapsis Ascending node Astrometry Central force Focus Kepler's second law Mean anomaly

All concepts