Which planet is closest to Earth, on average? I used to think it was Venus, followed by Mars. But this paper claims it is instead Mercury.
At first this seemed to make no sense. Just picture the orbits of the planets: they're a bunch of concentric circles (approximately). Venus' orbit completely encloses Mercury's. Every point of it is closer to the Earth's orbit than Mercury's orbit is. And indeed, that's how you get that Venus is the closest planet to Earth, followed by Mars, and then Mercury: take the difference between the radius of their orbits.
But that doesn't actually get you the average distance. If two planets happen to be lined up (at the same point in their orbit) then yes, the distance between them is the difference between their orbital radii. But if one of them is on the opposite side of the Sun as the other (at the opposite point in their orbit), then the distance between them is the sum of their radii, and Mercury is the closest planet to Earth!
So, to figure out what planet is closest to the Earth on average, you have to actually do the math.
Let's calculate the average distance between the Earth (circular orbit of radius ) and another planet (circular orbit of radius ). We'll suppose each planet has an independent uniform probability of being at each point in its orbit, but because of symmetry, we can take our coordinate system to have Earth in the x axis, and only worry about the varying position of the other planet relative to it.
The distance between the planets, given by the cosine rule, is . To find the average, we have to integrate this over between and , and divide by . That looks like a really unpleasant thing to integrate, but luckily our source paper tells us the answer, which turns out to be:
where is the "elliptic integral of the second kind". After looking up the formula for this elliptic integral on Wikipedia, I was able to wrangle our expression for the distance into the paper's formula:
Skippable Math
We start with the average distance by the cosine rule:
And we want to get to an expression involving the elliptic integral of the second kind, which is:
To turn that cosine into a sine squared, we substitute and use the identity , getting:
Now we have the sine squared, but the sign in front of it is positive (two minuses that cancel out), and we want it to be negative. So we substitute again, by , and use the fact that :
Rearranging a bit, we get:
After all the substitutions, the integral is between and . However, since is an even function, the parts of the integral from to and from to have the same value, so we can rewrite as:
Dividing and multiplying by , we get:
Where the integral finally looks like our elliptic integral, with equal to . Replacing it in the expression, it becomes:
Exactly the expression from the original paper!
According to the paper, this average distance strictly increases with the radius .[1] So, the lowest average distance is to the planet with the smallest orbit, i.e. Mercury. Problem solved!
...problem not really solved. While this does prove that Mercury is the closest planet to Earth, it doesn't actually help explain why. Is there a simple reason we can actually understand that explains why planets with smaller radii are closer to us?
Yes. Consider, instead of a random point on the inner planet's orbit, two points A and B, at the same angle above the vertical:
Between these two points, the average horizontal distance to the Earth is just the Earth's orbital radius, , and so doesn't depend on the other orbit's radius; and the vertical distance is the same for both points, . So increasing the radius doesn't change the average horizontal distance at all, and increases the average vertical distance; of course this means the average total distance is increased!
Of course, you may notice this was not completely valid reasoning, since horizontal distances and vertical distances don't add, they combine through Pythagoras' theorem. To turn this verbal insight into an actual proof, we need to write down the formula for the average between the two points A and B's distances, take the derivative with respect to , and see if it is positive. As it turns out, it is:
Skippable Math 2
The angle of the point B is , so its cosine is . The sum of the distances is then:
(The average is just half of this, and the factor of doesn't make a difference as to whether it's increasing with , so we're discarding it.)
The derivative of this expression with respect to is:
Which, replacing the expressions for and , is:
Since and are positive, this expression is greater than zero if and only if its denominator is positive, so we want:
Now, if in the expression for we replace with , its value becomes lesser or equal (with equality only when ):
And we can do the same for :
( is always positive because we're taking between and , so we don't need the absolute value here.)
Substituting in our equation, we get:
The term in the right side of this inequality is certainly greater than or equal to zero: if is less than zero, it's the sum of a term and its negative, which is 0; if it's more than zero, it's the sum of two positive terms, which is more than zero.
Since the inequality is strict when is not 0, our initial expression for the derivative may be zero when is zero, but is positive otherwise. This means that the sum of the distances grows with the radius , and so does their average.
The average distance over the entire circle is equivalent to an integral over averages like this (divided by ), with varying from to , so it also grows with the radius.
So the intuitive explanation does turn into a viable proof that Mercury really is, on average, the closest planet to Earth and to every other planet.
Actually, it doesn't even say that, it just says "the distance between two orbiting bodies is at a minimum when the inner orbit is at a minimum".