The Azimuth Project
Insolation

Contents

Idea

Insolation is short for ‘incoming solar radiation’. There are various precise things this might mean. Sometimes it means the daily average power of the solar radiation hitting the top of the Earth’s atmosphere at a given latitude and a given time of the year. This quantity is also known as Q¯ day. For details, see:

Details

Scientists studying the effect of Milankovitch cycles on the Earth’s glacial cycles often want to know the quantity Q¯ day at the summer solstice, typically at the latitude 65°. (Presumably the latitude here is a bit arbitrary, but the summer solstice makes some of the math easier). As you can see, this quantity changes in a complicated way with the passage of time:

The red line is the actual value of Q¯ day. The greenish-brown line is what Q¯ day would be if the eccentricity of the Earth’s orbit were zero. This graph was taken from Wikicommons. It as prepared by someone with the monicker ‘Incredio’. Should you believe him? Or should you check his calculation yourself?

To check it, you’ll need this formula for Q¯ day:

Q¯ day=S oπR o 2R E 2[sin(ϕ)sin(δ)+cos(ϕ)cos(δ)cos(h)]\overline{Q}^{day} = \frac{S_o}{\pi}\frac{R_o^2}{R_E^2}\left[ \sin(\phi) \sin(\delta) + \cos(\phi) \cos(\delta) \cos(h) \right]

This is derived in the above Wikipedia article using some spherical trigonometry. Here:

  • R 0 is the mean distance from the Earth to the Sun.

  • R E is the current distance from the Earth to the Sun (since the Earth’s orbit is an ellipse, this changes during the course of a year).

  • S 0 is the solar constant, about 1367W/m 2.

  • ϕ is the latitude.

  • δ is the declination.

  • h is the hour angle. ±h o is the hour angle when the sun comes up or down, obtained by solving cos(h o)=tan(ϕ)tan(δ) when 1tan(ϕ)tan(δ)1. When tan(ϕ)tan(δ)>1, the sun never sets and we say h o=π. For example, this may happen above the Arctic circle in summer. When tan(ϕ)tan(δ)<1, the sun never rises and we can skip to the end of the calculation, since we instantly know Q¯ day=0. This may happen above the Arctic circle in winter. Warning: The minus sign in the definition of the hour angle for sunrise/set may be due to a different definition of δ than as given above.

The following is partially based on something which Wikipedia has flagged as: “This article may be confusing or unclear to readers.”

Let θ be the polar angle describing the Earth’s position in its yearly orbit around the sun. For convenience, set θ=0 at the vernal equinox. Then the declination is given by

δ=εsin(θ)\delta = \varepsilon \sin(\theta)

where ε is the obliquity, i.e. the angle of title of the Earth’s axis relative to its orbital plane, where obliquity zero means its axis is at right angles to this plane.

The conventional longitude of perihelion ϖ is defined relative to the vernal equinox, so for an elliptical orbit:

R E=R o1+ecos(θϖ)R_E=\frac{R_o}{1+e\cos(\theta-\varpi)}

where e is the eccentricity of the Earth’s orbit. In other words:

R oR E=1+ecos(θϖ)\frac{R_o}{R_E}={1+e\cos(\theta-\varpi)}

At the summer solstice, things simplify:

  • The declination δ is equal to the obliquity ϵ,

and

  • θ=π/2, so R oR E=1+esin(ϖ). The quantity esin(ϖ) is called the precession index.

The longitude of perihelion ϖ, the obliquity ε and the eccentricity e change slowly due to the Milankovitch cycles. They are somewhat complicated functions of time, but luckily this website will compute them for you:

So, with some work, you could check Incredio’s graph for yourself!

For what happens after the Sun’s radiation hits the top of the atmosphere, see Solar radiation.

Code

category: climate