Article 72EW6 The center of the earth is not straight down

The center of the earth is not straight down

by
John
from John D. Cook on (#72EW6)

If the earth were a perfect sphere, down" would be the direction to the center of the earth, wherever you stand. But because our planet is a bit flattened at the poles, a line perpendicular to the surface and a line to the center of the earth are not the same. They're nearly the same because the earth is nearly a sphere, but not exactly, unless you're at the equator or at one of the poles. Sometimes the difference matters and sometimes it does not.

From a given point on the earth's surface, draw two lines: one straight down (i.e. perpendicular to the surface) and one straight to the center of the earth. The angle that the former makes with the equatorial plane is geographic latitude. The angle that the latter makes with the equatorial plane is geocentric latitude.

For illustration we will draw an ellipse that is far more eccentric than a polar cross-section of the earth.

geographic_geocentric_ellipse.png

At first it may not be clear why geographic latitude is defined the way it is; geocentric latitude is conceptually simpler. But geographic latitude is easier to measure: a plumb bob will show you which direction is straight down.

There may be some slight variation between the direction of a plumb bob and a perpendicular to the earth's surface due to variations in surface gravity. However, the deviations due to gravity are a couple orders of magnitude smaller than the differences between geographic and geocentric latitude.

Conversion formulas

The conversion between the two latitudes is as follows.

geographic_geocentric.svg

Heree is eccentricity. The equations above work for any ellipsoid, but for earth in particular e^2 = 0.00669438.

The function atan2(y, x) returns an angle in the same quadrant as the point (x, y) whose tangent isy/x. [1]

As a quick sanity check on the equations, note that when eccentricity e is zero, i.e. in the case of a circle, = . Also, if = 0 then = for all eccentricity values.

Next we give a proof of the equations above.

Proof

We can parameterize an ellipse with semi-major axisa and semi-minor axisb by

geographic_geocentric1.svg

The slope at a point (x(t),y(t)) is the ratio

geographic_geocentric2.svg

and so the slope of a line perpendicular to the tangent, i.e tan , is

geographic_geocentric3.svg

Now

geographic_geocentric4.svg

and so

geographic_geocentric5.svg

where e^2 = 1 - b^2/a^2 is the eccentricity of the ellipse. Therefore

geographic_geocentric6.svg

and the equations at the top of the post follow.

Difference

For the earth's shape, e^2 = 0.006694 per WGS84. For small eccentricities, the difference between geographic and geocentric latitude is approximately symmetric around 45.

geographic_diff1.png

But for larger values of eccentricity the asymmetry becomes more pronounced.

geographic_diff2.png

Related posts

[1] There are a couple complications with programming language implementations of atan2. Some call the function arctan2 and some reverse the order of the arguments. More on that here.

The post The center of the earth is not straight down first appeared on John D. Cook.
External Content
Source RSS or Atom Feed
Feed Location http://feeds.feedburner.com/TheEndeavour?format=xml
Feed Title John D. Cook
Feed Link https://www.johndcook.com/blog
Reply 0 comments