Can we know the shape of the Earth?
orHow Al-Biruni Calculated the Circumference of the Earth Using a Mountain in the 11th century
As you might recall from my earlier posts, I've used the Sun at local Solar Noon (when the Sun is at its zenith) on the Equinox as a way to know the shape of the Earth. But all we need is some point that is high enough to be seen at some significant distance and Mountain peaks are much easier to experiment with than a moving Sun.
If the ground is a flat plane then we can simply triangulate using the mountain peak at different distances and get the SAME answer for the height using simple trigonometry:
For example, if the ground was flat, we could find h by simply measuring our angle α at two different distances, where Δd is the distance between the two observation points using the formula:
Δd * tan(α₀) * tan(α₁)
h = ----------------------
tan(α₁) - tan(α₀)
Can you prove that?
This seems complex but lets think about what the tangent function REALLY does. You put in an ANGLE and you get back the matching SLOPE. And SLOPE is simply RISE over RUN (rise/run). That's why we have the relationship: tan(θ) = opposite/adjacent -- once you understand that it is simply converting an angle into a slope I think things get so much simpler. All we are really doing here is a unit conversion.
And since it is rise/run we can expand this equation by replacing the tangent function with the equation for each of our slopes and we get:
(d₀-d₁) * (h/d₀) * (h/d₁)
-------------------------
h/d₁ - h/d₀
-------------------------
h/d₁ - h/d₀
Now what we want to do is prove (mathematically) that this equation would give us the value for h if we plugged in all the values. To do that we need to show that we can simplify the equation and everything will cancel out, leaving us the value for h at the end.
First, let's multiply out our numerator and move the denominator over to keep things easier to read (dividing by something is the same as multiplying by the inverse):
h*h*(d₀-d₁) 1
------------- * ------------
d₀ * d₁ h/d₁ - h/d₀
Next, we can join the denominator on the left by multiplying both fractions by the Least Common Denominator (which is [d₀*d₁], and simplify our new denominator):
Next, we can join the denominator on the left by multiplying both fractions by the Least Common Denominator (which is [d₀*d₁], and simplify our new denominator):
h*h*(d₀-d₁) d₀ * d₁
------------- * ----------
d₀ * d₁ h*(d₀-d₁)
Now we can see that [d₀ * d₁] cancels out leaving:
h*h*(d₀-d₁)
-------------
h*(d₀-d₁)
And we can cancel out the common term [h*(d₀-d₁)] from the numerator and denominator which leaves us with just:
h
[Or let Wolfram|Alpha do it]
So we have proven mathematically that our starting (slightly complex) equation will give us the value we need. Our problem is we don't know the value of h to find h/d₀ directly but since we do know what the value of h/d₀ is (because we have the angle and h/d₀=tan(α₀)) we can plug in the values and get the answer back out of the equation. So we don't even need the distance to the mountain in this case, just the distance between our two observations and the angle to the peak is sufficient to calculate the height AS LONG AS our ground plane is flat.
However, when you try this over longer distances, you find that it doesn't quite work. As you get further away from the mountain the mountain appears to shrink. The further away you get the more it shrinks. If you try it with a streetlamp you'll find that it works very well because the ground is flatter than your measurement error over short distances. But with mountains or tall buildings you begin to measure the curvature of the Earth appearing to make the object shorter and shorter with distance.
And the amount it shrinks is pretty consistent no matter which direction you go from the mountain so we know the curvature is pretty even in all directions and at all points around the Earth where we try to measure tall objects.
It is because of this, almost** equal shrinkage in all directions, at all points around the Earth, that we can conclude the Earth is a spheroid (and with better measurement accuracy we can tell it isn't EXACTLY equal in all directions but we can tell that it's a little bit oblate or squished making it a little bit wider around the equator than around the poles).
** On a sphere it would be equal, on our slightly oblate/squished ellipsoid Earth it varies slightly, mainly with latitude.
So what we find is that our ground is NOT flat because our angles keep changing with distance as the ground is rotating underneath us, throwing off our calculation:
This doesn't mean that triangulation is useless (but it fairly conclusively, confirms the shape of the Earth as being curved). But what we need to be able to do is take our measurements from close enough to the mountain and from some height so we can still see the same base point on the mountain so that our baseline is flat, based on a direct line-of-sight.
Ok, so we can use this to calculate the height of the mountain as long as we make sure each observation point measures the angle between the same spot at the base of the mountain and the peak AND we now know that the Earth is curved.
What Al-Biruni figured out in the 11th century was that if he could find the height of the mountain peak as above (using angles and distances) AND find the angle between a level on the peak and the horizon (shown as angle β below) then this would also be the angle made at the center of the circle, and from this he could estimate the Radius of the Earth, and with the Radius of the Earth it is easy to find the circumference (r·2π).
Here is the geometry that Al-Biruni had in mind, at point A at the peak of the mountain we measure the angle of the dip to the horizon (greatly exaggerated here, it will usually be a smaller angle):
sin(x) = opposite / hypotenuse
So what we find is that our ground is NOT flat because our angles keep changing with distance as the ground is rotating underneath us, throwing off our calculation:
This doesn't mean that triangulation is useless (but it fairly conclusively, confirms the shape of the Earth as being curved). But what we need to be able to do is take our measurements from close enough to the mountain and from some height so we can still see the same base point on the mountain so that our baseline is flat, based on a direct line-of-sight.
Ok, so we can use this to calculate the height of the mountain as long as we make sure each observation point measures the angle between the same spot at the base of the mountain and the peak AND we now know that the Earth is curved.
What Al-Biruni figured out in the 11th century was that if he could find the height of the mountain peak as above (using angles and distances) AND find the angle between a level on the peak and the horizon (shown as angle β below) then this would also be the angle made at the center of the circle, and from this he could estimate the Radius of the Earth, and with the Radius of the Earth it is easy to find the circumference (r·2π).
Here is the geometry that Al-Biruni had in mind, at point A at the peak of the mountain we measure the angle of the dip to the horizon (greatly exaggerated here, it will usually be a smaller angle):
Instead of a tangent, this time we're going to look at the sin function.
What we know about sin is that it is also a ratio, just like our tangent function, consider a right triangle:
What 'opposite' means here is the length of the side opposite the angle in question - so we have the following relationships:
sin(A) = a/c
sin(B) = b/c
sin(C) = c/c
From this it is easy to show that the Law of Sines holds:
a b c
------ = ----- = ------
sin(A) sin(B) sin(C)
By simply expanding our formula, for example:
b b c
--- = -*- = c
b/c 1 b
We also know that the sum of all the angles must equal 180° and that one of our angles is 90° so when we measure the angle β we also can know that α is simply 90°-β and we also know that sin(90°-β) = cos(β) [I'll leave this as an exercise for the reader].
Now armed with these identities we can write out the equations for this diagram of the Earth using our Law of Sines:
b c
----- = ------
sin(B) sin(C)
And we substitute in our known values, where:
b = R
sin(B) = sin(90°-β) = cos(β)
c = R+h
sin(C) = sin(90°) = 1
So that gives us a starting equation where we can solve for R.
R
------- = R + h
cos(β)
Isolate R to one side by subtracting R(cos(β)/cos(β)) from both sides, you will see why we added the (cos(β)/cos(β)) in the next step:
R R cos(β) R cos(β)
------- - -------- = R - h - --------
cos(β) cos(β) cos(β)
Which simplifies to:
1
R(------ - 1) = h
cos(β)
We can then multiply both sides by the inverse, or 1/(1/cos(β)-1) rewritten as cos(β)/(1-cos(β)) which gives us:
h cos(β)
R = --------
1-cos(β)
Which can also be written:
h
R = -------- = h/(sec(β)-1)
sec(β)-1
So now we have the portion of our triangle represented by R in terms of the height of our mountain and the angle between level at the top of the mountain and a distant point on the horizon.
Let's say we're up on a Mountain at 2000 meters. The dip angle to the horizon should be about 1.43546°. Plugging that into our formula:
2000/(sec(1.43546°)-1) = 6.371×10⁶
What if we're up at 33,709 feet?
Let's say we're up on a Mountain at 2000 meters. The dip angle to the horizon should be about 1.43546°. Plugging that into our formula:
2000/(sec(1.43546°)-1) = 6.371×10⁶
We would expect this angle to be about 3.25° - and since the iPhone Field of View is about 63.5° and that image is about 1077 pixels high we expect that to be about 55 pixels high (this is not exactly true because the lens is not perfectly linear, but it's a small section near the center so this is a good approximation). PLEASE NOTE: The angle in question here is the Elevation Angle - the "Horizon Angle" is how much the horizon is tilted clockwise/counterclockwise.
What do you think, does that look like about where the horizon is? It's hard to tell because of the clouds and atmospheric distortion in this image but hey, it's what I have (the horizon most definitely does not begin at the white section).
Next time you are going on a flight grab a Theodolite app and maybe you'll get a better view of the horizon and you can estimate the circumference of the Earth. Can you beat Al-Biruni's value which was correct to within ~200 miles?
Post your results in the comments section!
I have several flights coming up this month so I'll try to get some better data to share.