Astronomy & Celestial Mechanics
the sky, computed
Generator

The doppler-tracking generator

An eight-hour pass, and a 351 m s⁻¹ sinusoid that is the whole of the angle
An eight-hour pass, and a 351 m s⁻¹ sinusoid that is the whole of the angle. Above: the range rate a two-way Doppler measurement returns over one pass from Goldstone, for a spacecraft receding at 14.6 km s⁻¹. Nothing here is an angle. The measurement is the fractional shift of a carrier the spacecraft coherently turned around and sent back, and its interpretation is that the distance is changing at some rate. Below: the same data with the spacecraft's own smooth signature removed. What is left is a sinusoid of exactly one cycle per day — the station's own motion, carried east at 379 metres a second by the rotation of the Earth, projected onto the line of sight. Its amplitude is that speed times cos δ and returns a declination of 22.0°; its zero crossing is the moment the spacecraft passed the meridian and returns the right ascension. The Earth's rotation is the interferometer. With Doppler good to 0.05 mm s⁻¹ at a 60-second cadence, 480 samples fit that amplitude to 0.003 mm s⁻¹ and the declination to 23 nanoradians — which is 4.7 milliarcseconds, from an instrument with no image plane and no angular resolution of any kind. What the picture cannot show is the part that makes this hard in practice: the spacecraft's own signature is not a straight line but a trajectory with unmodelled accelerations in it, and separating a slow non-gravitational force from a slow drift in the angles is the whole art of the fit.

Above: the range rate a two-way Doppler measurement returns over one pass from Goldstone, for a spacecraft receding at 14.6 km s⁻¹. Nothing here is an angle. The measurement is the fractional shift of a carrier the spacecraft coherently turned around and sent back, and its interpretation is that the distance is changing at some rate. Below: the same data with the spacecraft's own smooth signature removed. What is left is a sinusoid of exactly one cycle per day — the station's own motion, carried east at 379 metres a second by the rotation of the Earth, projected onto the line of sight. Its amplitude is that speed times cos δ and returns a declination of 22.0°; its zero crossing is the moment the spacecraft passed the meridian and returns the right ascension. The Earth's rotation is the interferometer. With Doppler good to 0.05 mm s⁻¹ at a 60-second cadence, 480 samples fit that amplitude to 0.003 mm s⁻¹ and the declination to 23 nanoradians — which is 4.7 milliarcseconds, from an instrument with no image plane and no angular resolution of any kind. What the picture cannot show is the part that makes this hard in practice: the spacecraft's own signature is not a straight line but a trajectory with unmodelled accelerations in it, and separating a slow non-gravitational force from a slow drift in the angles is the whole art of the fit.

7 essays call doppler-tracking. 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.

An eight-hour pass, and a 351 m s⁻¹ sinusoid that is the whole of the angle. Above: the range rate a two-way Doppler measurement returns over one pass from Goldstone, for a spacecraft receding at 14.6 km s⁻¹. Nothing here is an angle. The measurement is the fractional shift of a carrier the spacecraft coherently turned around and sent back, and its interpretation is that the distance is changing at some rate. Below: the same data with the spacecraft's own smooth signature removed. What is left is a sinusoid of exactly one cycle per day — the station's own motion, carried east at 379 metres a second by the rotation of the Earth, projected onto the line of sight. Its amplitude is that speed times cos δ and returns a declination of 22.0°; its zero crossing is the moment the spacecraft passed the meridian and returns the right ascension. The Earth's rotation is the interferometer. With Doppler good to 0.05 mm s⁻¹ at a 60-second cadence, 480 samples fit that amplitude to 0.003 mm s⁻¹ and the declination to 23 nanoradians — which is 4.7 milliarcseconds, from an instrument with no image plane and no angular resolution of any kind. What the picture cannot show is the part that makes this hard in practice: the spacecraft's own signature is not a straight line but a trajectory with unmodelled accelerations in it, and separating a slow non-gravitational force from a slow drift in the angles is the whole art of the fit. Spaceflight

A position measured from a frequency

A spacecraft is unresolvable and unreachable, and everything known about where it is comes from two scalars — a round-trip light time and a Doppler shift. Neither is an angle. The orbit solution returns two angles anyway, because the antenna is bolted to a rotating planet.

 A nanosecond across the Earth is 35.69 nanoradians on the sky. The angular accuracy of a differenced-delay measurement against the length of the baseline it is measured on, both axes logarithmic, for three levels of delay precision. The relation is σ_θ = cσ_τ/B and nothing else, so every curve is a straight line of slope −1.00: the only two ways to measure an angle better are a better clock or a wider Earth, and only one of those is available. The three marked baselines are the ones that exist — the deep-space complexes in California, Spain and Australia, 8,400, 10,600, 11,700 kilometres apart. On the longest of them a delay good to 0.05 nanoseconds is 1.28 nanoradians, which at 0.52 astronomical units is 100 metres across the line of sight; a more typical 0.15-nanosecond measurement on the shortest baseline is 5.35 nanoradians. What makes any of this survivable is that the same pair of antennas observes a quasar a few degrees away immediately afterwards. The quasar is at infinity, its position is known better than the measurement, and subtracting its delay from the spacecraft's removes the clock offsets, the water vapour over each dish and the station coordinates in one step — so the number that comes out is not a delay at all but an angular separation from a fixed point in the sky. Spaceflight

An angle measured against a quasar

A tracking station measures how fast a spacecraft is receding, which is one number where three are wanted. The two missing angles come from the Earth's rotation, slowly, and near a planetary encounter there is no time for slowly — so the position is instead measured directly, as a difference of arrival times between two antennas, referred to a quasar a few degrees away.

 Two interiors, 2.59 mm/s apart, against a floor of 0.02. What a radio link measures when a spacecraft flies past a moon. The horizontal axis is time from closest approach in minutes and the vertical axis is the accumulated change in the line-of-sight velocity in millimetres per second, after the pull of the moon as a point mass has been fitted and removed. What is left is the part of the field that is not spherically symmetric, and the two curves are what two interiors predict for it. The upper one is a body whose tidal Love number is 0.616 — a shell floating on a global liquid layer, free to deform almost as a fluid would. The lower one is a body solid throughout, at 0.03. They differ by 2.59 millimetres per second in the accumulated deflection, against a floor of 0.02 for a coherent two-way X-band link integrated over tens of seconds: 130 standard deviations in a single pass. That ratio is the whole reason the measurement is possible, and it is worth stating what is being compared. The point-mass deflection itself is 837 metres per second — five orders of magnitude larger than the signature of interest — so the interior is not read off the Doppler curve but off the residual left after a model of everything larger has been subtracted: the moon's mass, the planet's, the spacecraft's own thrusting and outgassing, the plasma along the path, the station's motion, and relativity. Every one of those has to be right to a part in a hundred thousand before the last curve here means anything, which is why gravity science needs many passes and a global fit rather than one flyby and a subtraction. The published uncertainty on Titan's k₂ is eleven per cent rather than the fraction of a per cent this signal-to-noise would suggest, and the difference is entirely correlations with the other parameters in that fit — a reminder that a formal error on a curve is not the error on the number extracted from it. The straight-line path assumed here is exact only in the limit of a fast flyby; a slow one bends, and the bending is solved for rather than approximated. Spaceflight

An ocean found in a Doppler residual

A spacecraft flying past a moon is deflected by hundreds of metres a second, and the part of that deflection which says whether the moon has an ocean is two millimetres a second. Everything larger has to be modelled and removed first, exactly, and what is left over is an interior.

 Three passes, and the fourth is below the noise. How far the position is still wrong after each pass of the light-time iteration, for five targets, on a logarithmic scale of kilometres. The first pass uses the body's position now and is wrong by exactly the distance the body travels while its light is in transit — 6254 kilometres for Mars near opposition, which is about twenty-five arcseconds and larger than any residual in an ephemeris fit. Each further pass multiplies the error by the body's speed divided by the speed of light, so the lines are straight and their slopes are that ratio and nothing else. Three passes take every one of these below a metre, which is why the loop in an ephemeris code has a fixed trip count rather than a convergence test. The same argument run backwards is why the correction cannot be applied to the observation instead: the direction light arrived from is what the observation is, and the position it corresponds to is not known until the body's orbit is. Orbits

Where a planet is and where it is seen

An ephemeris is fitted to observations, and no observation is of a position. It is of a direction light arrived from, at a time that is not the time the light left, bent by a Sun that is nowhere near the line of sight. Three corrections stand between the two, and all three are larger than the residuals.

 Four known pieces of hardware, and the anomaly is the sum of them. The reported anomalous acceleration of a deep-space probe, in units of 10⁻¹⁰ metres per second squared, built up from the heat the spacecraft was known to be radiating. The generators put out about two and a half kilowatts of waste heat and sit on booms beside a large dish that reflects a share of it backwards, which is 62 per cent of the total on its own; the instrument compartment radiates through louvres on one face; the radio transmitter beams eight watts at the Earth, which is a torch pointing the wrong way. Sunlight is negligible this far out and is drawn to show that it is. The four sum to 8.65 against a measured 8.74 ± 1.33, and the agreement is the answer. What makes the episode worth keeping is that none of these numbers was discovered later: every one was in the spacecraft's own thermal documentation from before launch, and the model that produced the anomaly was a model of a point mass. Orbits

An acceleration that was the spacecraft's own heat

Two probes leaving the solar system were tracked for thirty years and both drifted from their predicted paths by a tenth of a nanometre per second squared. The residual was real, it was constant, and it was the same on both. It was also the waste heat of the reactors that powered them, radiating slightly more one way than the other.

 TCB has gained 23.5 seconds on TT since 1977. How far four of the solar system's time scales have drifted apart, against years since they were set equal. TT is the time a clock on the Earth's geoid keeps, and it is flat here by definition. TCG is the time a clock at rest just outside the Earth's gravitational well would keep, and it runs faster by seven parts in ten thousand million — 0.0220 seconds a year. TCB is the time a clock at rest outside the Sun's well would keep, and it runs faster by fifteen parts in a thousand million, or 0.4893 seconds a year. TDB is TCB with that rate divided out so that it stays within milliseconds of TT, which is what makes it usable as an ephemeris argument and what makes it a coordinate rather than a proper time. The differences are constants and they are not corrections: an ephemeris tabulated against one of these and evaluated against another is wrong by the whole of this plot. The observed sky

A second that depends on where the clock is

A clock in orbit runs fast, a clock at the bottom of a gravity well runs slow, and a clock at rest far from the Sun runs faster than anything on Earth by fifteen parts in a thousand million. Astronomy therefore has four different seconds, two of them longer than the others, and an ephemeris has to declare which one it is tabulated against.

 Four media, and only two of them care about the frequency. The four propagation delays between a tracking station and a spacecraft, in metres of round-trip range, against the angle between the target and the Sun. The troposphere contributes a few metres and depends only on the elevation; the ionosphere is smaller at X band and scales as the inverse square of the frequency; the solar plasma rises steeply towards conjunction and scales the same way; and the Sun's gravitational delay rises as a logarithm and has no frequency dependence at all. That last distinction is the whole of the calibration strategy: transmitting and receiving at two widely separated frequencies measures the two plasma terms and removes them, leaving the troposphere to a weather model and the relativistic term to theory. Near conjunction the plasma exceeds everything else by an order of magnitude, which is why a spacecraft passing behind the Sun is both untrackable and the best available laboratory for measuring the relativistic term. Spaceflight

Four media between the antenna and the spacecraft

A deep-space range measurement is a round-trip time, and the signal spends that time passing through a troposphere, an ionosphere, the solar wind and a curved region of spacetime. All four delay it, none of them is the orbit, and the whole of the navigation depends on removing them.

The whole library · All essays