I didn’t manage to cover everything I intended last time, so I’m moving the stuff about the eccentricity of the Earth’s orbit to this week, and expanding it.

### Sunshine and the Earth’s orbit

I bet some of you are hungry for some math. As I mentioned, it takes some work to see how changes in the eccentricity of the Earth’s orbit affect the annual average of sunlight hitting the top of the Earth’s atmosphere. Luckily Greg Egan has done this work for us. While the result is surely not new, his approach makes nice use of the fact that both gravity and solar radiation obey an inverse-square law. That’s pretty cool.

Here is his calculation with some details filled in.

Let’s think of the Earth as moving around an ellipse with one focus at the origin. Its angular momentum is then

where is its mass, and are its polar coordinates, and is the angular component of its velocity:

So,

and

Since the brightness of a distant object goes like , the solar energy hitting the Earth per unit time is

for some constant It follows that the energy delivered per unit of angular progress around the orbit is

Thus, the total energy delivered in one period will be

So far we haven’t used the the fact that the Earth’s orbit is elliptical. Next we’ll do that. Our goal will be to show that depends only very slightly on the eccentricity of the Earth’s orbit. But we need to review a bit of geometry first.

### The geometry of ellipses

If the Earth is moving in an ellipse with one focus at the origin, its equation in polar coordinates is

where is the **eccentricity** and is the somewhat dirty-sounding **semi-latus rectum**. You can think of as a kind of average radius of the ellipse—more on that in a minute.

Let’s think of the origin in this coordinate system as the Sun—that’s close to true, though the Sun moves a little. Then the Earth gets closest to the Sun when is as big as possible. So, the Earth is closest to the Sun when , and then its distance is

Similarly, the Earth is farthest from the Sun happens when , and then its distance is

We call the **perihelion** and the **aphelion**.

The **semi-major axis** is half the distance between the opposite points on the Earth’s orbit that are farthest from each other. This is denoted These points occur at and , so the distance between these points is , and

So, the semi-major axis is the arithmetic mean of the perihelion and aphelion.

The **semi-minor axis** is half the distance between the opposite points on the Earth’s orbit that are closest to each other. This is denoted

**Puzzle 1.** Show that the semi-minor axis is the geometric mean of the perihelion and aphelion:

I said the semi-latus rectum is also a kind of average radius of the ellipse. Just to make that precise, try this:

**Puzzle 2.** Show that the semi-latus rectum is the harmonic mean of the perihelion and aphelion:

This puzzle is just for fun: the Greeks loved arithmetic, geometric and harmonic means, and the Greek mathematician Apollonius wrote a book on conic sections, so he must have known these facts and loved them. The conventional wisdom is that the Greeks never realized that the planets move in elliptical orbits. However, the wonderful movie *Agora* presents a great alternative history in which Hypatia figures it all out shortly before being killed! And the mathematician Sandro Graffi (who incidentally taught a course I took in college on the self-adjointness of quantum-mechanical Hamiltonians) has claimed:

Now an infrequently read work of Plutarch, several parts of the

Natural Historyof Plinius, of theNatural Questionsof Seneca, and of theArchitectureof Vitruvius, also infrequently read, especially by scientists, clearly show that the cultural elite of the early imperial age (first century A.D.) were fully aware of and convinced of a heliocentric dynamical theory of planetary motions based on the attractions of the planets toward the Sun by a force proportional to the inverse square of the distance between planet and Sun. The inverse square dependence on the distance comes from the assumption that the attraction is propagated along rays emanating from the surfaces of the bodies.

I have no idea if the controversial last part of this claim is true. But it’s fun to imagine!

More importantly for what’s to come, we can express the semi-minor axis in terms of the semi-major axis and the eccentricity. Since

we have

so the semi-minor axis is

while

so the semi-major axis is

and thus they are related by

Remember this!

### How total annual sunshine depends on eccentricity

We saw a nice formula for the total solar energy hitting the Earth in one year in terms of its angular momentum :

How can we relate the angular momentum to the shape of the Earth’s orbit? The Earth’s energy, kinetic plus potential, is constant throughout the year. The kinetic energy is

and the potential energy is

At the aphelion or perihelion the Earth isn’t moving in or out, just around, so by our earlier work

and the kinetic energy is

Equating the Earth’s energy at aphelion and perihelion, we thus get

and doing some algebra:

and solving for

But remember that the semi-major and semi-minor axis of the Earth’s orbit are given by

respectively! So, we have

This lets us rewrite our old formula for the energy in the form of sunshine that hits the Earth each year:

But we’ve also seen that

so we get the formula we’ve been seeking:

This tells us as a function of semi-major axis and eccentricity.

As we’ll see later, the semi-major axis is almost unchanged by small perturbations of the Earth’s orbit. The main thing that changes is the eccentricity . But if is small, is even smaller, so doesn’t change much when we change

We can make this more quantiative. Let’s work out how much the actual changes in the Earth’s orbit affect the amount of solar radiation it gets! As we’ll see, the semi-major axis is almost constant, so we can ignore that. Complicated calculations we can’t redo here show that the eccentricity varies between 0.005 and 0.058. We’ve seen the total energy the Earth gets each year from solar radiation is proportional to

When the eccentricity is at its lowest value, we get

When the eccentricity is at its highest value, we get

So, the solar power hitting the Earth each year changes by a factor of

In other words, it changes by merely 0.167%.

That’s very small And the effect on the Earth’s temperature would *naively* be even less!

Naively, we can treat the Earth as a **greybody**: an ideal object whose tendency to absorb or emit radiation is the same at all wavelengths and temperatures. Since the temperature of a greybody is proportional to the fourth root of the power it receives, a 0.167% change in solar energy received per year corresponds to a percentage change in temperature roughly one fourth as big. That’s a 0.042% change in temperature. If we imagine starting with an Earth like ours, with an average temperature of roughly 290 kelvin, that’s a change of just 0.12 kelvin!

The upshot seems to be this: *in a naive model without any amplifying effects, changes in the eccentricity of the Earth’s orbit would cause temperature changes of just 0.12 °C!*

This is much less than the roughly 5 °C change we see between glacial and interglacial periods. So, if changes in eccentricity are important in glacial cycles, we have some explaining to do. Possible explanations include season-dependent phenomena and climate feedback effects, like the ice albedo effect we’ve been discussing. Probably *both* are very important!

### Adiabatic invariance

Why does the semi-major axis of the Earth’s orbit remain almost unchanged under small perturbations? The reason is that it’s an ‘adiabatic invariant’. This is basically just a fancy way of saying it remains almost unchanged. But the point is, there’s a whole *theory* of adiabatic invariants… which supposedly explains the near-constancy of the semi-major axis.

According to Wikipedia:

The Earth’s eccentricity varies primarily due to interactions with the gravitational fields of Jupiter and Saturn. As the eccentricity of the orbit evolves, the semi-major axis of the orbital ellipse remains unchanged. From the perspective of the perturbation theory used in celestial mechanics to compute the evolution of the orbit, the semi-major axis is an adiabatic invariant. According to Kepler’s third law the period of the orbit is determined by the semi-major axis. It follows that the Earth’s orbital period, the length of a sidereal year, also remains unchanged as the orbit evolves. As the semi-minor axis is decreased with the eccentricity increase, the seasonal changes increase. But the mean solar irradiation for the planet changes only slightly for small eccentricity, due to Kepler’s second law.

Unfortunately, even though I understand a bit about the general theory of adiabatic invariants, I have not gotten around to convincing myself that the semi-major axis is such a thing, for the perturbations experienced by the Earth.

Here’s something easier: checking that the semi-major axis of the Earth’s orbit determines the period of the Earth’s orbit, say . To do this, first relate the angular momentum to the period by integrating the rate at which orbital area is swept out by the planet:

over one orbit. Since the area of an ellipse is , this gives us:

On the other hand, we’ve seen

Equating these two expressions for shows that the period is:

So, the period depends only on the semi-major axis, not the eccentricity. Conversely, we could solve this equation to see that the semi-major axis depends only on the period, not the eccentricity.

I’m treating and as constants here. If the mass of the Sun decreases, as it eventually will when it becomes a red giant and puffs out lots of gas, the semi-major axes of the Earth’s orbit will change. It will actually increase! This is one reason people are still arguing about just when the Earth will get swallowed up by the Sun:

• David Appell, The Sun will eventually engulf the Earth—maybe, *Scientific American*, 8 September 2008.

And, to show just how subtle these things are, if the mass of the Sun slowly changes, while the semi-major axis of the Earth’s orbit will change, the eccentricity will remain almost unchanged. Why? Because for *this* kind of process, it’s the *eccentricity* that’s an adiabatic invariant!

Indeed, I got all excited when I started reading a homework problem in Landau and Lifschitz’s book *Classical Mechanics*, which describes adiabatic invariants for the gravitational 2-body problem. But I was bummed out when they concluded that the eccentricity was an adiabatic invariant for gradual changes in . They didn’t discuss any problems for which the semi-major axis was an adiabatic invariant.

I’ll have to get back to this later sometime, probably with the help of a good book on celestial mechanics. If you’re curious about the concept of adiabatic invariant, start here:

• Adiabatic invariant, Wikipedia.

and then try this:

• Marko Robnik, Theory of adiabatic invariants, February 2004.

And if you know how to show the Earth’s semi-major axis is an adiabatic invariant, please tell me how!

A quick and dirty solution of puzzle 1. We need to maximize . For that purpose we use and solve the equation of the ellipse

It suffices to differentiate ‘below the square root’ and to solve the resulting linear equation for the maximum . Plugging this into the equation for yields the result.

Is your the same as my , i.e. the angle in polar coordinates centered at a focus of the ellipse? Or is it something else, like the angle in polar coordinates centered at the center of the ellipse? My elliptology is sort of rusty.

Yes, I meant your polar coordinates. As one can see, my elliptology is nonexistent.

It is very pretty how this discussion of the ellipse incorporates five commonly encountered special cases of the generalized mean of r1 and r2 (arithmetic, geometric, harmonic, min (r1 itself), and max (r2 itself)).

Is it too much to hope that the sixth commonly encountered special case of the generalized mean (namely, root mean square) can be naturally incorporated as well?

I don’t know! This book mentions the

root mean squareorquadratic meanof the perihelion and aphelion :and says that it equals

where is the semi-major axis and is the eccentricity, but it doesn’t give any geometrical interpretation of this quantity. It just does so for the arithmetic, geometric and harmonic means.

I would like to know who first discovered these facts about means and ellipses. If Apollonius didn’t know them, I’ll be very disappointed.

Thanks, very interesting. BTW an apparent typo in your original posting:

1) 2nd occurrence of “so the semi-major axis is”: change major to minor

I’ll check that out and fix it if needed. Thanks!

1) The semi-latus rectum can be geometrically interpreted as a line from a focal pt to the ellipse, drawn parallel to the minor axis. Five down, one to go:-)

2) I too couldn’t find (in the sense of “search”, not in the sense of “think hard” like Apollonius) a geometrical interpretation of the Quadratic mean of r1 and r2, though I suspect one exists.

3) A nice geometric interpretation of the six means – min, H, G, A, Q, & max – of two numbers, in the context of a circle rather than of a general ellipse, is at: http://en.wikipedia.org/w/index.php?title=File:MathematicalMeans.svg&page=1

4) Nothing major but the major you made minor should be major and the major you kept major should be minor:-)

By the way, in class I titled one section of my talk ‘Sunshine and eccentricity’, and mentioned that this was a good description of California.

And perhaps Florida in equal degree.

If you accept that the sum of the distances between any point on an ellipse and the two foci is a constant, Puzzle 1 is easy.

Let be half the distance between the foci, the perihelion distance and the aphelion distance. Then the sum of distances to the foci at perihelion is , the sum of distances to the foci at aphelion is , and equating them gives us and .

When the planet lies on the minor axis, which bisects the line between the foci, by symmetry and Pythagoras the semi-minor axis is given by:

The nicest proof that the sum of the distances to the foci is constant involves the Dandelin spheres, and you can see it illustrated

here.

Puzzle 2 is straightforward algebra from , , but there’s a typo in the formula for the harmonic mean, which should be:

The Milankovich cycles that most closely match this are two cycles in the eccentricity of the Earth’s orbit which have periods of 95 and 123 thousand years. But as we saw last time, these have very tiny effects on the average solar energy hitting the Earth year round […]