Astronomy & Celestial Mechanics
the sky, computed
Concept

Oblateness — where it appears

The flattening of a rotating body, whose leading gravitational effect is a quadrupole term that turns the orientation of any orbit around it. The Earth's is a thousandth of its field and produces essentially all of the secular drift a satellite mission has to design around.

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

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 Why a shell pulls a point inside it not at all. A double cone from a point inside a uniform shell. In the narrow-cone limit the far patch is 2.80 times further away and 2.80 times wider — the same number, checked to a part in a thousand before this figure is drawn. Its mass is greater by the square of that ratio and its pull weaker by the same square, so the two cancel in every direction.

A sphere pulls exactly like a point, and the proof is a pair of cones

Every orbit ever computed treats the Sun as a dot. That is not an approximation — for a spherical body it is exact, and the reason is a cancellation between two patches of a shell.

gravitation · Shell theorem 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 three revolutions, on a turning Earth. The ground track of a circular orbit at 420 km and 51.64° inclination, over 3 revolutions, on an equirectangular graticule. The latitude is a sine wave bounded by ±51.64° exactly — sin φ = sin i sin u, so the inclination is the highest latitude the orbit ever passes over, and it is reached twice per revolution. Each successive pass is displaced west by (ω⊕ − Ω̇) × 92.90 min = 23.61°, of which 0.32° is the orbital plane's own regression and the rest is the planet turning underneath: the vehicle comes back to nearly the same place in inertial space and the place has moved. The period used is the nodal one, 92.899 min against the Keplerian 92.970: J₂ makes the two differ by 4.31 s, which is 0.018° of walk per revolution and 102° in a year — the difference between a repeat track and a track that used to repeat. The map is equirectangular and therefore wrong about area everywhere; what it is right about is longitude difference, which is the whole of what this figure measures.

The line under a satellite

A ground track is an orbit seen from a frame that is turning, so every pass lands west of the last one. The track closes only when two periods are commensurable — which turns "look at the same place every day" into a condition on the altitude.

spaceflight · Ground tracks Measured moment-of-inertia factors, from the Sun's 0.07 to the Moon's 0.3931. Eleven bodies whose interiors have never been sampled, arranged by the one interior quantity that has been measured for all of them. C/MR² is 2/5 for a uniform sphere and falls as mass is concentrated toward the centre, and the values here span from 0.07 to 0.3931. That spread is the content. The Moon at 0.3931 is barely differentiated — whatever iron core it has is a few per cent of its radius, which is why the Moon is the one large body in the inner solar system without a magnetic field of its own. Mercury at 0.346 is nearly as low as the Earth despite being an eighth of its mass, and for a body that small the only way to get there is an iron core filling most of the radius. The Sun at 0.07 is off the scale of anything a two-layer model describes; a star is not a planet with a bigger core but a body whose density falls by five orders of magnitude between centre and surface. The faint curves behind are the two-layer relation at a few density contrasts, drawn to show what kind of interior each value is consistent with — and the horizontal placement of each body on them is an illustration rather than a result, since one factor never fixes one core. Every number here was obtained by watching the body turn: a precession rate, a libration amplitude, a gravity field sampled on a flyby, or in the Sun's case the frequencies of its own oscillations.

Whether the heavy material sank

A moment of inertia is 0.4 of MR² for a uniform sphere and less for everything that has differentiated, and it is measured by watching a body wobble. Mercury's 0.346 says most of the planet is iron core; the Moon's 0.393 says almost none of it is.

gravitation · Moment of inertia Where each gravity harmonic gets its signal: J₂ from the bulk, J₁₀ from the outer 12 per cent. Why a spacecraft that never enters a planet can say something about its depth. Each zonal harmonic of the external field is an integral over the interior density weighted by r to the power of the degree plus two, and the curves here are those integrals accumulated outward: the fraction of each coefficient that has been contributed by the time the integration reaches a given fractional radius, for an interior of polytropic index 1. The weighting climbs steeply with degree, so the curves separate. Half of J₂ comes from inside 73 per cent of the radius, and half of J₁₀ from inside 88 per cent — the higher coefficients barely know the deep interior exists. That ordering is the whole basis of gravity science as a probe. A single coefficient is one number and constrains almost nothing; a series of them, each weighted differently, is a coarse depth profile, and it is how Jupiter's core turned out to be smeared over half the planet rather than sitting as a distinct sphere at the middle. Two limits are worth stating with it. The information falls off fast: by degree ten the kernel is concentrated in a shell so thin that measuring the coefficient says little about anything below it. And every curve here assumes north–south symmetry, under which the odd harmonics vanish identically — so a measured J₃ or J₅ is not a deeper probe of the same thing but a measurement of something else entirely, which at a giant planet is how fast and how deep the winds run.

A core weighed by something that never went in

The external gravity field of a planet is a series, and each term of it is an integral over the interior density weighted by a different power of radius. Measure enough terms and the series becomes a coarse depth profile — which is how Jupiter's core turned out to be smeared over half the planet rather than sitting at the middle.

gravitation · Oblateness A star 1.24 times wider than it is tall, and 19 per cent brighter pole-on. Left, the meridional section of a star rotating at ω = 0.9257 of its critical angular velocity, computed from the Roche potential rather than sketched: the equator sits at 1.245 polar radii, and at the critical rate that ratio is exactly 1.5 whatever the star is made of. The same rotation expressed as a fraction of the critical equatorial speed is 0.768, and the two conventions differ by the distortion itself — a figure that prints one under the other's name is wrong by an amount that looks like rounding. Effective gravity at the equator is 0.329 of its polar value, so von Zeipel's flux law makes the pole hotter than the equator by a factor 1.320 at the theoretical exponent 0.25 and 1.232 at the 0.188 that interferometric imaging actually fits. Right, the apparent bolometric brightness against viewing inclination, integrated over the visible gravity-darkened surface: pole-on the star is 1.19 times brighter than edge-on, and the apparent temperature falls with it. The consequence is that a rapid rotator's place on the Hertzsprung–Russell diagram is partly a statement about the observer's position, which no spectrum taken alone can undo.

A temperature that depends on where the observer stands

A star turning near its break-up rate is half again as wide as it is tall, and its equator is thousands of degrees cooler than its poles. Neither of those is a small correction to a spectrum — the effective temperature and the luminosity such a star appears to have are partly statements about which way its axis happens to point.

starlight · Gravity darkening A flattened shell that pulls nothing inside it. The strength of the gravitational field of a thin homoeoid — the shell between two similar concentric ellipsoids, here with axis ratio 0.5 — along its long axis (solid) and its short axis (dashed), in units of GM/a² against distance from the centre in units of the long semi-axis a, from 20,000 equal masses placed on it. The dotted curve is a point mass. Inside the shell the field is zero along both axes to within 0.0001 — the numerical residue of a finite sampling — however flattened the shell is: this is Newton's theorem for a homoeoid, the flattened generalisation of the spherical shell's empty interior. Outside, the shell does not pull like a point: at 1.2 a along the long axis the field is 1.003 against the point's 0.694, and at 0.60 a along the short axis, just outside the shell there, 0.901. The exterior field remembers the shape; the interior does not know there is one.

The shell can be flattened if it thins at the poles

A spherical shell exerts no force anywhere inside it, and a ring exerts a force that points outward. Between them lies a shape Newton found that keeps the sphere's empty interior however far it is flattened: the homoeoid, a shell bounded by two similar ellipsoids, thick at the equator and thin at the poles in exactly the right proportion. It is why the gravity inside any flattened body can be computed at all — and why the planes of stellar streams measure the shapes of dark haloes.

gravitation · Shell theorem The plane a moon's orbit turns about, between equator and ecliptic. The inclination of the Laplace plane — the plane about which a satellite's orbit precesses — to the Earth's equator, against distance from the Earth in Earth radii, for an obliquity of 23.44°. Close in, the Earth's equatorial bulge dominates and the plane is the equator; far out, the Sun dominates and it is the ecliptic, 23.44° away. The dashed curve uses today's bulge, which makes the transition at 8.4 Earth radii. The solid curve lets the bulge grow with the faster spin the Earth had when the Moon was closer — a 5.3-hour day at 5 Earth radii, from the system's angular momentum alone — which moves the transition out to about 14. The Moon today, at 60.3 Earth radii, has a Laplace plane within 0.0022° of the ecliptic, so its 5.15° inclination is entirely free: nothing forces it.

A tilt the Moon should not have

A moon's orbit precesses about a plane that sits on its planet's equator close in and on the planet's own orbit far out. A Moon born in the Earth's equator and carried outward by the tides should have followed that plane all the way to the ecliptic, because the angle it keeps to it is an adiabatic invariant. It is inclined by five degrees instead, and slow recession cannot have put them there.

orbits · Angular momentum A two-body arc flown in the real field. How far from its target a spacecraft arrives when it flies the velocity Lambert's two-body problem gives, through the Earth's real field with its equatorial bulge, against the angle of the transfer, from a 400 km orbit to a point on a 420 km orbit, on a logarithmic scale. In a pure two-body field every solution arrives within a metre of the target; exactly 180° is left out, where the two points and the Earth's centre are in line and the transfer's plane is not determined. In the real field the miss grows with the transfer: at 51.6° inclination, 4.2 km after 60° and 43.0 km after 300°. An equatorial orbit misses most, because the bulge's extra pull is strongest over the equator and acts on it all the way round; an inclined orbit feels it only part of the time but also has its plane turned.

An answer from two bodies, flown through a flattened Earth

Lambert's problem gives the exact velocity that carries a spacecraft between two points in a given time — exactly, for a spherical Earth. The real Earth bulges at the equator, by a thousandth of its central pull, and a spacecraft flying the exact two-body answer between two points in low orbit arrives tens of kilometres from where it was sent. The repair is to keep the two-body solver and lie to it about the target, and the lie converges by a factor of a hundred each time it is told.

orbits · Lambert's problem An iteration that slows toward the straight-through transfer. The contraction ratio of the corrected-target iteration — how much each pass shrinks the miss — against the transfer angle, for a transfer between orbits of 51.6° inclination at 400 and 420 km, with the Earth's bulge at its real strength and at 10 and 30 times it, on a logarithmic scale. The iteration converges where the ratio is below 1 (dashed). At the real strength it shrinks the miss 1857-fold per pass at 150°, 83-fold at 178° and 8-fold at 181°: always convergent, ever slower. At 10 times the strength it fails between 181° and 184°; at 30 times it fails between 182.5° and 192°. The slowing is geometric: near 180° the transfer's plane is set by a sideways shift of the aim point, so correcting the bulge's sideways miss tilts the whole arc.

A correction that tilts the arc it meant to aim

A transfer designed with Lambert's two-body solution and flown through the Earth's flattened field misses its target, and the standard repair is to aim at the target minus the miss and solve again. That converges fast — until the transfer angle approaches 180°. There the two end points and the Earth's centre are nearly in line, the arc's plane is decided by a sideways shift of the aim, and correcting the bulge's sideways miss rotates the whole arc. The iteration slows by a factor of a hundred, and with a larger perturbation it runs away. Correcting only within the plane fixes it, at the cost of a plane change paid at arrival.

orbits · Lambert's problem Every orbit turns by about the same angle, whatever its size. How far the Earth's equatorial bulge turns an ellipse with a 300-km perigee in one revolution, in degrees, against its semi-major axis, logarithmic: the line of apsides (solid) and the node (dashed), for inclinations of 28.5, 51.6, 98°. The turn per orbit depends on the semi-latus rectum, which for an elliptical orbit lies between one and two perigee radii, so it changes by less than a factor of four from a near-circular orbit to one reaching the Moon's distance: at 28.5°, 0.56° of apsides per orbit at 8,000 km and 0.20° at 300,000 km. Per day the large orbit turns far less, because it completes so few revolutions; per revolution, almost the same.

An orbit lost to a burn that did not fire

A weak stage leaving the Earth splits its departure into several burns at perigee, coasting round a staircase of ever larger ellipses between them. If one burn does not fire, the stage goes round once more. That costs a revolution of time — an hour and a half early in the staircase, more than a day at the top — and a revolution of turning by the Earth's bulge, which comes out at nearly the same angle for every step. The angle is cheap to undo at the slow apogee of a large orbit. The time is not, when the target is the Moon.

spaceflight · Orbital transfer

Named alongside it

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

J2InclinationNodal precessionPerturbationsShell theoremSun-synchronous orbitCorrected target iterationDifferential correctionGeoidGravity scienceHydrostatic equilibriumLambert's problem

All concepts