Astronomy & Celestial Mechanics
the sky, computed
Generator

The orbit-averages generator

Four averages of one distance, and the two of them that are the semi-major axis
Four averages of one distance, and the two of them that are the semi-major axis. The average distance of a body from its primary, against eccentricity and in units of the semi-major axis, computed four ways: averaged over time, over true anomaly, over eccentric anomaly, and as the harmonic mean in time. Every curve is a quadrature over the orbit — 2,048 panels uniform in eccentric anomaly, with Kepler's equation supplying the time weight — and not a closed form. Two of the four are exactly a at every eccentricity, which is why they are drawn as one line: the eccentric-anomaly average, because the mean of cos E over a turn is zero, and the harmonic mean in time, because the time weight cancels 1/r at every node before the sum begins. The other two are not: the time average is a(1 + e²/2), which rises to 1.4050 a at e = 0.9, and the true-anomaly average is a√(1−e²) — the semi-minor axis — which falls to 0.4359 a there. So a is the average distance in two senses out of four, and the ordering b ≤ a ≤ ⟨r⟩ₜ holds at every eccentricity with equality only on the circle. At Earth's e = 0.0167 the four agree to 0.014%, and at Mercury's e = 0.2056 the spread is 2.14%. The distinction is invisible for the planets and unavoidable for a comet, and it is the reason a quoted "mean distance" has to say which mean.

The average distance of a body from its primary, against eccentricity and in units of the semi-major axis, computed four ways: averaged over time, over true anomaly, over eccentric anomaly, and as the harmonic mean in time. Every curve is a quadrature over the orbit — 2,048 panels uniform in eccentric anomaly, with Kepler's equation supplying the time weight — and not a closed form. Two of the four are exactly a at every eccentricity, which is why they are drawn as one line: the eccentric-anomaly average, because the mean of cos E over a turn is zero, and the harmonic mean in time, because the time weight cancels 1/r at every node before the sum begins. The other two are not: the time average is a(1 + e²/2), which rises to 1.4050 a at e = 0.9, and the true-anomaly average is a√(1−e²) — the semi-minor axis — which falls to 0.4359 a there. So a is the average distance in two senses out of four, and the ordering b ≤ a ≤ ⟨r⟩ₜ holds at every eccentricity with equality only on the circle. At Earth's e = 0.0167 the four agree to 0.014%, and at Mercury's e = 0.2056 the spread is 2.14%. The distinction is invisible for the planets and unavoidable for a comet, and it is the reason a quoted "mean distance" has to say which mean.

2 essays call orbit-averages. The drawing above is what it returns with no arguments at all; every call below passes it something, because a placement that passes nothing draws whichever member of the family the generator happens to default to rather than the one its essay argues about.

Where it is called

Every figure listed here is the same construction drawn at different numbers, so a correction to one is a correction to all of them.

Four averages of one distance, and the two of them that are the semi-major axis. The average distance of a body from its primary, against eccentricity and in units of the semi-major axis, computed four ways: averaged over time, over true anomaly, over eccentric anomaly, and as the harmonic mean in time. Every curve is a quadrature over the orbit — 2,048 panels uniform in eccentric anomaly, with Kepler's equation supplying the time weight — and not a closed form. Two of the four are exactly a at every eccentricity, which is why they are drawn as one line: the eccentric-anomaly average, because the mean of cos E over a turn is zero, and the harmonic mean in time, because the time weight cancels 1/r at every node before the sum begins. The other two are not: the time average is a(1 + e²/2), which rises to 1.4050 a at e = 0.9, and the true-anomaly average is a√(1−e²) — the semi-minor axis — which falls to 0.4359 a there. So a is the average distance in two senses out of four, and the ordering b ≤ a ≤ ⟨r⟩ₜ holds at every eccentricity with equality only on the circle. At Earth's e = 0.0167 the four agree to 0.014%, and at Mercury's e = 0.2056 the spread is 2.14%. The distinction is invisible for the planets and unavoidable for a comet, and it is the reason a quoted "mean distance" has to say which mean. Orbits

The average depends on what is being averaged

Four ways of averaging one orbit's distance from its primary give four different numbers, and only two of them are the semi-major axis. Which two is not a matter of convention, and the same arithmetic decides how much sunlight a planet receives in a year.

 June sunlight at 65°N, against where perihelion sits. Daily-mean insolation at latitude 65 degrees north on the June solstice, at an obliquity of 23.44 degrees, against the longitude of perihelion measured from the March equinox — the angle that precesses right round in about twenty-one thousand years. Three eccentricities are drawn. The swing is ±10.0 per cent at e = 0.05 and ±1.0 per cent at e = 0.005, in proportion to the eccentricity, because the Sun–Earth distance on a fixed date carries e cos of the precession angle and that is first order. Over the same range of eccentricity the annual mean at this latitude moves by 0.12 per cent, because the annual mean carries 1/√(1−e²) and that is second order. The two together are the whole of the precession term in Milankovitch's theory: the eccentricity does almost nothing to how much sunlight the Earth receives and a great deal to when it arrives, and the ice sheets of the northern hemisphere respond to the summer they might melt in rather than to the year's total. It also explains why the precession signal disappears when the orbit is nearly circular: multiply a large angular swing by a vanishing eccentricity and there is nothing left, which is what the innermost curve here is. Orbits

An average that precession cannot move

Sunlight arrives as the inverse square of the distance and time passes as its square, so the two cancel exactly in a year's integral. The longitude of perihelion therefore changes the annual mean insolation at every latitude by precisely nothing — and changes June at 65°N by ten per cent.

The whole library · All essays