Earth Radius by Latitude (WGS 84)

This online calculator calculates Earth radius at given latitude using WGS 84 reference ellipsoid

The calculator below calculates Earth radius at given latitude. In fact, of course, it calculates the radius of WGS 84 reference ellipsoid at given latitude, and if you want some theory recap, you can find it below the calculator.


Earth Radius by Latitude (WGS 84)


Digits after the decimal point: 3
Radius (km)
Save the calculation to reuse next time or share with friends.

Earth radius

Since the Earth is flattened at the poles and bulges at the equator, geodesy represents the shape of the Earth with an oblate spheroid. The oblate spheroid, or oblate ellipsoid, is an ellipsoid of revolution obtained by rotating an ellipse about its shorter axis. It is the regular geometric shape that most nearly approximates the shape of the Earth. A spheroid describing the figure of the Earth or other celestial body is called a reference ellipsoid. The reference ellipsoid for Earth is called an Earth ellipsoid.

Earth's physical surface is irregular. It can be approximated by the geoid, which was an important concept for almost two hundred years in the history of geodesy and geophysics. According to Gauss, who first described it, it is the "mathematical figure of Earth", a smooth but highly irregular surface whose shape results from the uneven distribution of mass within and on the surface of Earth. The geoid surface is irregular, but considerably smoother than Earth's physical surface.

Because of their relative simplicity, reference ellipsoids are used as a preferred surface on which geodetic network computations are performed and point coordinates such as latitude, longitude, and elevation are defined. Currently the most common reference ellipsoid used, and that used in the context of the Global Positioning System, is the one defined by WGS 84.

An ellipsoid of revolution is uniquely defined by two quantities. Several conventions for expressing the two quantities are used in geodesy, but they are all equivalent to and convertible with each other:

  • Equatorial radius a (called semi-major axis), and polar radius b (called semi-minor axis);
  • a and eccentricity e;
  • a and flattening f.

WGS 84 defines ellipsoid parameters as:
Semi-major axis a = 6378137.0 meters
Semi-minor axis b = 6356752.3142 meters

Point on the ellipsoid surface can be defined by the parametric curve equation
(x,y)=(a  \, cos(t), b \, sin(t))

Radius can be found using the Pythagorean theorem
R(t)^2 = a^2 cos^2(t) + b^2 sin^2(t)

The problem, however, is that the angle t from the example above is the geocentric latitude, and coordinates provided in geodetic datum, like WGS 84, are geodetic. Geodetic latitude is determined by the angle between the equatorial plane and normal to the ellipsoid, whereas geocentric latitude is determined by the angle between the equatorial plane and line joining the point to the centre of the ellipsoid (see figure).

Geocentric vs geodetic latitude
Geocentric vs geodetic latitude

So, to find the radius, we need to relate geodetic latitude \alpha to geocentric latitude \beta.

Let's start from tangent of our curve, which can be obtained by differentiating curve equation.
(x,y)'=(\frac{dx}{dt}, \frac{dy}{dt}) = ( -a  \, sin(t), b \, cos(t))

This is the vector which points along the curve (along T-line on the figure).

We can rotate it clockwise 90 degrees (x,y) => (y, -x) and get normal vector (b \, cos(t), a  \, sin(t)) which points along N-line.

Parameter t is our \alpha. The slope of normal vector is also tangent of the angle \beta. Thus
tan(\beta)=\frac{a  \, sin(\alpha)}{b \, cos(\alpha)}=\frac{a}{b}tan(\alpha)

Using the relations between tangent and cosine
1+tan^2(\alpha)=\frac{1}{cos^2(\alpha)} => cos^2(\alpha)=\frac{1}{1+tan^2(\alpha)}
and between tangent and sine
1+cotan^2(\alpha)=\frac{1}{sin^2(\alpha)} => sin^2(\alpha)=\frac{1}{1+\frac{1}{tan^2(\alpha)}}=\frac{tan^2(\alpha)}{1+tan^2(\alpha)},
we can rewrite formula for radius as
R^2 = a^2 \frac{1}{1+tan^2(\alpha)} + b^2 \frac{1}{1+\frac{1}{tan^2(\alpha)}} = \frac{a^2}{1+tan^2(\alpha)} + b^2 \frac{tan^2(\alpha)}{1+tan^2(\alpha)}
and replace tangent of \alpha with tangent of \beta expression
R^2 = \frac{a^2}{1+\frac{b^2}{a^2}tan^2(\beta)} + b^2\frac{\frac{b^2}{a^2}tan^2(\beta)}{1+\frac{b^2}{a^2}tan^2(\beta)}=\frac{a^4+b^4tan^2(\beta)}{a^2+b^2tan^2(\beta)}

Then, we can simplify it a bit
R^2 =\frac{a^4+b^4tan^2(\beta)}{a^2+b^2tan^2(\beta)} = \frac{a^4cos^2(\beta)+b^4sin^2(\beta)}{a^2cos^2(\beta)+b^2sin^2(\beta)}= \frac{ (a^2cos(\beta))^2+(b^2sin(\beta))^2}{(a\,cos(\beta))^2+(b\,sin(\beta))^2}

And, finally, get the formula from wikipedia

R =\sqrt{\frac{ (a^2cos(\beta))^2+(b^2sin(\beta))^2}{(a\,cos(\beta))^2+(b\,sin(\beta))^2}}

The calculator above uses this formula.