The Azimuth Project



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\overline{Q}^{day}. For details, see:


Scientists studying the effect of Milankovitch cycles on the Earth’s glacial cycles often want to know the quantity Q¯ day\overline{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\overline{Q}^{day}. The greenish-brown line is what Q¯ day\overline{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\overline{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:

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

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

where ε\varepsilon 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 ϖ\varpi 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 ee 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 longitude of perihelion ϖ\varpi, the obliquity ε\varepsilon and the eccentricity ee 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.


category: climate