The observed sky

The shadow that climbs the zenith

After sunset the sky overhead is lit only above the Earth's own shadow, and the shadow climbs as the square of the Sun's depression. Scattering computed once from that geometry predicts a sky that dims faster and faster and is as dark as night by nine degrees. The real sky takes eighteen, and the difference is light that has been scattered more than once.

Assumes Twilight and Refraction.

The three definitions of twilight are angles — six, twelve and eighteen degrees of solar depression — and what they stand for is a brightness that the angles do not compute. The sky after sunset is lit by sunlight scattered in the air overhead, and how bright it is at a given depression is a question about where that sunlight can still reach. The angle is the proxy; the scattering is the thing.

The scattering can be computed from one piece of geometry and one assumption, and the result is instructive twice over. The geometry predicts a sky that darkens faster and faster, lit by a thin layer of air that climbs as the Sun sinks. The assumption — that each photon is scattered exactly once — predicts that the sky is as dark as a moonless night by nine degrees of depression. The real sky is not dark until eighteen, and the gap between those two numbers turns out to be the most important thing the calculation measures.

The zenith sky after sunset, as single scattering predicts it. The brightness of the zenith sky at 550 nm, in magnitudes per square arcsecond with brighter upward, against the Sun's depression below the horizon, computed by single scattering in a spherical atmosphere: US Standard Atmosphere densities, Rayleigh scattering, a 300 Dobson-unit ozone layer, and every photon scattered exactly once. The Earth's shadow climbs the zenith as R(sec d − 1) — 8.7 km at 3°, 35.1 at 6°, 142 at 12° — through air that thins by a factor of e every eight kilometres or so, so the brightness does not fall steadily: it falls 1.4 magnitudes a degree between 3° and 6°, and 2.7 a degree between 7° and 10°. The model gives 14.3 at the end of civil twilight and 27.0 at the end of nautical. The dashed line is the natural night sky, 21.9 magnitudes per square arcsecond. Single scattering reaches it at 9.25° of depression — 8.8 degrees before the 18° at which astronomical twilight is observed to end. The difference is not an error in the arithmetic; it is the light this model leaves out, scattered more than once.
Fig. 1 The brightness of the zenith sky at 550 nm, in magnitudes per square arcsecond with brighter upward, against the Sun’s depression, computed by single scattering in a spherical atmosphere with standard densities, Rayleigh scattering and a 300 Dobson-unit ozone layer. The Earth’s shadow climbs the zenith as R(sec d − 1): 8.7 km at 3°, 35.1 at 6°, 142 at 12°. The brightness falls 1.4 magnitudes a degree between 3° and 6° and 2.7 a degree between 7° and 10°, reaching 14.3 at 6° and 27.0 at 12°. It reaches the natural night sky, 21.9, at 9.25° — 8.8 degrees before twilight is observed to end.

The shadow over the zenith

After sunset the Earth stands between the observer’s sky and the Sun, and it casts a shadow into its own atmosphere. Directly overhead that shadow has a floor that is easy to compute. A point at height zz above the observer is in sunlight if the straight line from it towards the Sun clears the Earth. With the Sun a depression dd below the horizon, that line passes closest to the Earth’s centre at a distance (R+z)cosd(R+z)\cos d, so the point is lit only if that exceeds RR — only above

zs=RcosdR    Rd22.z_s = \frac{R}{\cos d} - R \;\approx\; \frac{R\,d^2}{2}.

Everything in the zenith column below zsz_s is in the shadow and contributes nothing. Everything above it is sunlit and scatters some of that light down.

The quadratic is the whole story in miniature. At three degrees of depression the shadow over the zenith is 8.7 kilometres up, below the height of the highest mountains’ summits; at six degrees, 35.1 kilometres, above nearly all the air; at twelve, 142 kilometres, where the atmosphere is a few billionths of its density at the ground. Doubling the depression quadruples the shadow’s height, and the air thins exponentially with height, so every degree of depression switches off more of the lit air than the degree before.

The Earth’s shadow is not what makes the Moon’s phases, but it is what makes twilight, and the dark blue band that rises from the eastern horizon after sunset — the shadow seen edge-on against the air — is the same surface, viewed from the side rather than from underneath. The pinkish band that sits just above it, the Belt of Venus, is air still in sunlight on the far side of the observer, lit by rays that have crossed so much atmosphere on the way that little but the red end of the spectrum is left to scatter back.

The same formula decides what else overhead is still sunlit, and two of its consequences can be seen without an instrument. Noctilucent clouds, the thin ice clouds of the upper atmosphere, float at around 83 kilometres, and directly overhead the shadow does not reach that height until the Sun is about nine degrees down — which is why they are seen glowing against a darkening sky, well into nautical twilight, long after every weather cloud below them has gone grey. A satellite is higher again. The shadow does not reach the International Space Station’s 420 kilometres over the zenith until the Sun is about twenty degrees down, past the end of astronomical twilight, so a pass overhead on a summer evening can be sunlit and brilliant against a sky that is already fully dark.

How the brightness is computed

The number on the vertical axis is an absolute surface brightness, and getting it takes three ingredients and no adjustable constant.

The first is the air itself: its density at every height, taken from the standard atmosphere, and its scattering cross-section, which for molecules much smaller than the wavelength is Rayleigh’s, falling as the fourth power of the wavelength. The second is the path. Each lit point on the zenith line receives sunlight that has travelled past the Earth along a grazing line, and the attenuation of that sunlight is the optical depth integrated along the actual line — dense near its tangent point, thin at both ends. The light scattered from the point then travels straight down to the observer and is attenuated again by the air below. The third is the zero point. The ratio of the sky’s radiance to the Sun’s own is the scattering integral multiplied by the solid angle the Sun subtends, divided by the full sphere it scatters into, so the Sun’s surface brightness, −10.59 visual magnitudes per square arcsecond, turns the integral directly into a magnitude on the same scale as the night sky.

Adding the contributions of every lit point on the line gives the sky’s brightness, and the same sum, kept point by point, gives where that brightness comes from. Nothing in it is fitted to a twilight observation, which is what allows the result to be compared with one.

Why the dimming accelerates

If the air thinned by a factor of ee every scale height HH, and nothing attenuated the light on its way, the brightness of the zenith would be proportional to the air above the shadow, ezs/He^{-z_s/H}. Written in magnitudes that is 1.086zs/H1.086\,z_s/H, and since zsRd2/2z_s \approx Rd^2/2 the magnitude grows as the square of the depression. Its slope — magnitudes per degree — therefore grows in proportion to the depression itself.

The computed curve shows exactly that. Between three and six degrees it falls 1.4 magnitudes a degree; between seven and ten, 2.7. A constant rate of dimming would be a straight line on the figure; single scattering draws a curve that bends ever more steeply down.

That shape is the first place the calculation can be held against the sky, and it does not match. The conventional description of the real twilight sky is a dimming by about a factor of a hundred for every six degrees — five magnitudes in six degrees, close to a straight line and much less steep than the single-scattering curve at the depressions where astronomers care. Something is keeping the real sky brighter as the Sun goes further down, and doing so more and more as the directly lit layer climbs out of the air.

Where the light comes from

The same calculation can report not only the total but its source: how much each kilometre of the zenith column contributes.

Where along the zenith the twilight comes from. The contribution each kilometre of the zenith column makes to the sky's brightness, for the Sun at 2, 4, 6, 8° below the horizon, each curve scaled to its own peak. Below the Earth's shadow nothing is lit; just above it the sunlight has grazed the lowest, densest air and been heavily attenuated; far above it there is too little air to scatter. The light therefore comes from a layer a few kilometres thick above the shadow, and the layer climbs as the Sun sinks: at 2° the shadow is at 3.9 km and the brightest layer at 15.1 km; at 4° the shadow is at 15.6 km and the brightest layer at 27.1 km; at 6° the shadow is at 35.1 km and the brightest layer at 48.3 km; at 8° the shadow is at 62.6 km and the brightest layer at 75.2 km. Each degree of depression moves the layer up into air a few times thinner than the last, which is why the brightness falls faster and faster: the twilight sky is a column of air being switched off from the bottom.
Fig. 2 The share of the zenith’s brightness contributed by each kilometre of height, for the Sun at 2°, 4°, 6° and 8° below the horizon, each curve scaled to its own peak. The dashed lines are the shadow’s height. At 2° the shadow is at 3.9 km and the brightest layer at 15.1 km; at 4°, 15.6 and 27.1 km; at 6°, 35.1 and 48.3 km; at 8°, 62.6 and 75.2 km. The light comes from a layer a few kilometres thick sitting roughly a dozen kilometres above the shadow.

The layer has a floor and a ceiling for two different reasons. Just above the shadow the air is densest, but the sunlight reaching it has skimmed the Earth through the thickest part of the atmosphere on the way, and most of it has been scattered out or absorbed before it arrives. Well above the shadow the light arrives almost unattenuated, but there is little air left to scatter it. Between the two the product peaks, about a dozen kilometres up — a little more than one scale height above the shadow, which is where a grazing ray’s optical depth has fallen to about one.

So twilight is not the whole sky slowly dimming. It is a thin, bright sheet of air riding a dozen kilometres above a rising shadow, and the brightness at any moment is the density of the air that sheet happens to occupy.

The sheet keeps climbing after civil twilight, and past a hundred kilometres the air it moves through changes character.

Where along the zenith the twilight comes from. The contribution each kilometre of the zenith column makes to the sky's brightness, for the Sun at 9, 10, 11, 12° below the horizon, each curve scaled to its own peak. Below the Earth's shadow nothing is lit; just above it the sunlight has grazed the lowest, densest air and been heavily attenuated; far above it there is too little air to scatter. The light therefore comes from a layer a few kilometres thick above the shadow, and the layer climbs as the Sun sinks: at 9° the shadow is at 79.4 km and the brightest layer at 90.8 km; at 10° the shadow is at 98.3 km and the brightest layer at 129.0 km; at 11° the shadow is at 119.2 km and the brightest layer at 151.9 km; at 12° the shadow is at 142.3 km and the brightest layer at 175.6 km. Each degree of depression moves the layer up into air a few times thinner than the last, which is why the brightness falls faster and faster: the twilight sky is a column of air being switched off from the bottom.
Fig. 3 The same construction at 9°, 10°, 11° and 12° of depression. The shadow reaches 79.4, 98.3, 119.2 and 142.3 km, and the brightest layer 90.8, 129.0, 151.9 and 175.6 km. The gap between the two grows from about 11 km to about 33 km, because above a hundred kilometres the air’s density falls off far more slowly with height than it does lower down.

Above about a hundred kilometres the atmosphere is heated by the Sun’s ultraviolet, its temperature climbs to many hundreds of degrees, and its scale height stretches from eight kilometres to tens. The directly lit layer spreads out accordingly. That is the thermosphere, and it is the same thin gas whose density decides how fast a satellite’s orbit decays — and whose density a working model can get wrong by a factor of two when the Sun is active. The last single-scattered light of evening comes from the air that drags on satellites.

Nine degrees, not eighteen

The dashed line on the first figure is the brightness of the natural night sky, 21.9 magnitudes per square arcsecond at a dark site — a floor made of airglow, zodiacal light and unresolved stars that no amount of waiting lowers. Twilight is over, for any practical purpose, when the scattered sunlight falls below it.

Single scattering puts that moment at 9.25 degrees of depression. The sky is observed to reach its night-time floor at about eighteen, which is why eighteen is the definition of astronomical twilight. The calculation is missing half of twilight.

The missing half is light scattered more than once. After the zenith column has dropped into the Earth’s shadow, the air towards the western horizon — along the line to the Sun, at much lower altitude — is still in sunlight and still bright. That bright glow is itself a light source, and some of its light is scattered a second time by the shadowed air overhead and sent down to the observer. It is weaker than the direct sunlight was, by roughly the fraction of light the first scattering sent in the right direction, but it does not care about the shadow’s height over the zenith, because it arrives sideways. As the directly lit sheet climbs into air too thin to matter, the secondarily lit air below it takes over, and the dimming slows to the near-straight line the real sky follows.

The eight or nine degrees between where single scattering ends and where twilight actually ends is a measurement of the atmosphere’s multiple scattering, made with nothing but a clock and a photometer. It is also why the curve’s shape is the wrong test of the model and its end point is the right one: a model that includes only once-scattered light cannot be made to last eighteen degrees by adjusting any density, because at nine degrees there is almost no air left above the shadow for it to use.

What the angles mean in minutes

A depression angle becomes a time of evening through the latitude and the season, and the first essay on twilight computed that conversion directly.

Twilight at latitude 52°. Solar altitude through the second half of the day at 52°, on June solstice, equinox, December solstice, with the four thresholds that define twilight marked. Sunset is at −0.833° rather than at 0° because the Sun's own semi-diameter and the atmosphere's refraction together lift it by that much when it is geometrically already down. Civil twilight ends at −6°, when the brightest stars appear and outdoor work stops; nautical at −12°, when the horizon can no longer be seen against the sky and a sextant becomes useless; astronomical at −18°, when the sky stops contributing to a photometric measurement. Their durations at equinox here are 34 min, 40 min, 42 min — twilight is not a fixed length, it is the reciprocal of how steeply the Sun descends, and the Sun descends at an angle of about 90° − φ to the horizon. At this latitude the curve for June solstice never reaches −18° at all: astronomical twilight does not end, and there is no astronomically dark night from about 48.6° upward.
Fig. 4 Solar altitude through the second half of the day at 52°, on the June solstice, the equinox and the December solstice, with the four thresholds that define twilight marked. At the equinox civil twilight takes 34 minutes, nautical 40 more and astronomical 42 more. At this latitude the June solstice curve never reaches −18° at all.

At 52 degrees on the equinox the Sun reaches six degrees of depression 34 minutes after its apparent setting and twelve degrees 40 minutes after that, so single scattering’s nine and a quarter degrees arrive roughly fifty-five minutes after sunset. The sky is not dark for another hour. Everything after the first fifty-five minutes of an equinoctial evening at this latitude is multiply scattered light — the whole of nautical twilight past its midpoint and the whole of astronomical twilight.

What ozone does to the brightness

The ozone layer, concentrated around twenty-odd kilometres, absorbs light across the middle of the visible spectrum in a broad, weak band. For a vertical ray it removes a few per cent. For a grazing ray the story is different, because the sunlight that lights the sheet has crossed the ozone layer near its tangent point, along a path dozens of times longer than a vertical one.

The zenith sky after sunset, as single scattering predicts it. The brightness of the zenith sky at 550 nm, in magnitudes per square arcsecond with brighter upward, against the Sun's depression below the horizon, computed by single scattering in a spherical atmosphere: US Standard Atmosphere densities, Rayleigh scattering, no ozone, and every photon scattered exactly once. The Earth's shadow climbs the zenith as R(sec d − 1) — 8.7 km at 3°, 35.1 at 6°, 142 at 12° — through air that thins by a factor of e every eight kilometres or so, so the brightness does not fall steadily: it falls 1.4 magnitudes a degree between 3° and 6°, and 2.8 a degree between 7° and 10°. The model gives 13.3 at the end of civil twilight and 26.6 at the end of nautical. The dashed line is the natural night sky, 21.9 magnitudes per square arcsecond. Single scattering reaches it at 9.50° of depression — 8.5 degrees before the 18° at which astronomical twilight is observed to end. The difference is not an error in the arithmetic; it is the light this model leaves out, scattered more than once.
Fig. 5 The same zenith brightness at 550 nm with the ozone layer removed and everything else unchanged. The shadow is at the same heights and the dimming has the same shape, falling 1.4 magnitudes a degree between 3° and 6° and 2.8 between 7° and 10°. The sky is brighter: 13.3 at 6° against 14.3 with ozone, and it reaches the night sky at 9.50° instead of 9.25°.

A magnitude at six degrees is a factor of two and a half, removed by a layer whose total amount would be three millimetres thick if it were compressed to the pressure at the ground. The grazing geometry multiplies its effect, and it does so in a way that depends on colour, because ozone’s absorption band sits in the orange and red.

The zenith sky after sunset, as single scattering predicts it. The brightness of the zenith sky at 450 nm, in magnitudes per square arcsecond with brighter upward (on the scale of the Sun's visual surface brightness, so the zero point is exact only at 550 nm and the shape is exact everywhere), against the Sun's depression below the horizon, computed by single scattering in a spherical atmosphere: US Standard Atmosphere densities, Rayleigh scattering, a 300 Dobson-unit ozone layer, and every photon scattered exactly once. The Earth's shadow climbs the zenith as R(sec d − 1) — 8.7 km at 3°, 35.1 at 6°, 142 at 12° — through air that thins by a factor of e every eight kilometres or so, so the brightness does not fall steadily: it falls 1.4 magnitudes a degree between 3° and 6°, and 2.7 a degree between 7° and 10°. The model gives 13.4 at the end of civil twilight and 26.2 at the end of nautical. The dashed line is the natural night sky, 21.9 magnitudes per square arcsecond. Single scattering reaches it at 9.50° of depression — 8.5 degrees before the 18° at which astronomical twilight is observed to end. The difference is not an error in the arithmetic; it is the light this model leaves out, scattered more than once.
Fig. 6 The zenith brightness with ozone, computed at 450 nm, on the same magnitude scale. It is 13.4 at 6° and 26.2 at 12°, and reaches the night-sky level at 9.50°. At this wavelength ozone absorbs very little, and Rayleigh scattering is more than twice as strong as at 550 nm.
The zenith sky after sunset, as single scattering predicts it. The brightness of the zenith sky at 650 nm, in magnitudes per square arcsecond with brighter upward (on the scale of the Sun's visual surface brightness, so the zero point is exact only at 550 nm and the shape is exact everywhere), against the Sun's depression below the horizon, computed by single scattering in a spherical atmosphere: US Standard Atmosphere densities, Rayleigh scattering, a 300 Dobson-unit ozone layer, and every photon scattered exactly once. The Earth's shadow climbs the zenith as R(sec d − 1) — 8.7 km at 3°, 35.1 at 6°, 142 at 12° — through air that thins by a factor of e every eight kilometres or so, so the brightness does not fall steadily: it falls 1.5 magnitudes a degree between 3° and 6°, and 2.8 a degree between 7° and 10°. The model gives 14.1 at the end of civil twilight and 27.5 at the end of nautical. The dashed line is the natural night sky, 21.9 magnitudes per square arcsecond. Single scattering reaches it at 9.50° of depression — 8.5 degrees before the 18° at which astronomical twilight is observed to end. The difference is not an error in the arithmetic; it is the light this model leaves out, scattered more than once.
Fig. 7 The same construction at 650 nm, where ozone still absorbs and Rayleigh scattering is weaker. The zenith is 14.1 at 6° and 27.5 at 12°, and reaches the night-sky level at 9.50°. Against the 450 nm sky it is 0.7 magnitudes fainter at six degrees of depression.

The blue sky at six degrees is seven tenths of a magnitude brighter than the red. That is a colour, and it is not the colour Rayleigh scattering alone would produce: scattering favours blue by the fourth power of the frequency, but a long grazing path strips blue from the incoming sunlight by the same law, and the two very nearly cancel. What tips the balance is the ozone, and why the twilight zenith is blue at all is a question the colour of the sky answers better than its brightness.

A screen below the shadow

Real air is not clean all the way down. Cloud, haze and dust in the lowest dozen kilometres can be opaque to a grazing ray, and the effect is to lift the shadow: a ray must now clear the top of the opaque layer rather than the ground.

The zenith sky after sunset, as single scattering predicts it. The brightness of the zenith sky at 550 nm, in magnitudes per square arcsecond with brighter upward, against the Sun's depression below the horizon, computed by single scattering in a spherical atmosphere: US Standard Atmosphere densities, Rayleigh scattering, a 300 Dobson-unit ozone layer, an opaque layer below 12 km that no grazing ray can cross, and every photon scattered exactly once. The Earth's shadow climbs the zenith as (R + 12 km) sec d − R — 20.8 km at 3°, 47.2 at 6°, 155 at 12° — through air that thins by a factor of e every eight kilometres or so, so the brightness does not fall steadily: it falls 1.4 magnitudes a degree between 3° and 6°, and 2.7 a degree between 7° and 10°. The model gives 14.5 at the end of civil twilight and 27.0 at the end of nautical. The dashed line is the natural night sky, 21.9 magnitudes per square arcsecond. Single scattering reaches it at 9.00° of depression — 9.0 degrees before the 18° at which astronomical twilight is observed to end. The difference is not an error in the arithmetic; it is the light this model leaves out, scattered more than once.
Fig. 8 The zenith brightness at 550 nm with an opaque layer below 12 km that no grazing ray can cross, so that the shadow over the zenith is raised by the layer’s height. The sky at 6° is 14.5 against 14.3 without the screen, 27.0 at 12°, and it reaches the night-sky level at 9.00° instead of 9.25°.

A twelve-kilometre screen changes the sky at six degrees by two tenths of a magnitude and brings single scattering’s end forward by a quarter of a degree. That is much less than lifting the shadow by twelve kilometres would suggest, and the reason is in the profile figure: the rays that grazed below twelve kilometres were already so attenuated that the brightness hardly depended on them. The screen removes light that was mostly lost anyway. The same logic explains why the refraction that lifts the setting Sun matters so little to the twilight sky once the Sun is a few degrees down: the rays it bends are the ones that have been all but extinguished.

What the calculation leaves out

Multiple scattering, as its central omission. Every number above is once-scattered light, which is why it ends twilight at nine degrees. The eighteen-degree twilight of the real sky is the calculation’s own measure of what it lacks.

Aerosols. Particles comparable in size to the wavelength scatter strongly forward, weakly sideways and with little dependence on colour. They brighten the sky near the Sun, redden sunsets, and add a component whose height distribution varies by the hour. The model is clean air.

The Sun’s size and the view direction. The Sun is a disc half a degree across, so the shadow’s edge is a penumbra rather than a line, and the model looks only straight up. The sky near the western horizon, where most twilight light is, is a different and more complicated calculation.

And extinction on the way down is only the vertical part. The scattered light reaching the observer from the zenith crosses the whole atmosphere once; an observer looking lower crosses more, and every surface brightness here would be fainter in any other direction.

What single scattering is good for

The calculation is wrong about when twilight ends, and it is right about everything that makes it wrong. It says where the light comes from, at every depression, to within a kilometre or two: a sheet riding a dozen kilometres above a shadow that climbs as the square of the angle. It says why the directly lit sky must collapse by nine degrees, because by then the sheet is in the thermosphere. And by collapsing where the real sky does not, it isolates the one process — light scattered more than once — that carries the second half of every evening.

The same kind of calculation, applied to another piece of physics in the same geometry, answers a question the brightness curves raised in passing: why the evening sky overhead is blue when the path the light has taken should have made it anything but.

Still open: why the zenith stays blue

The grazing ray that lights the twilight sheet crosses the ozone layer almost tangentially, and ozone absorbs orange and red. Without it, the two effects of Rayleigh scattering — blue scattered in, blue stripped out on the way — nearly cancel, and the sky overhead would be a pale, colourless glow. With it, the zenith is blue. Past that is the same geometry on worlds whose scattering layers are taller, thinner or made of dust, where the angle over which the sky darkens is set by the height of the scatterers against the size of the planet.

What this makes readable

Essays that name this one as a prerequisite.

About the same objects

Not linked from either essay — found by the objects both name.

What links here

Essays that link to this one from their own argument.

The objects this essay names

Each one links to every other essay that touches it.

Depression angleEarths shadowMultiple scatteringOzoneRayleigh scatteringScale heightSingle-scatteringSky brightnessThermosphereTwilight