Area of a “rectangle” on a globe
This post will derive the area of a spherical region bounded by lines of latitude and longitude. Such a region corresponds to an actual rectangle on a Mercator projection map, with sides aligned with the coordinate axes, and is approximately a rectangle on a sphere if the rectangle is not too big.
What do we know up front?Before we get into detailed equations, we know that the area will be proportional to the difference in longitude. If we're looking that the area between two parallels, such as the equator and the Arctic Circle, the area between 10 and 20 longitude is the same as the area between 80 and 90 longitude, and twice the area between 72 and 77 longitude.
The difficulty is latitude. Say we look at squares on a map that are 1 of longitude wide and 1 of latitude tall. Those squares are on the map will correspond to more area on the globe for latitudes near the equator, and less area at high latitudes.
So the area bounded by longitudes 1 and 2 and latitudes 1 and 2 will depend on 1 and 2 individually, but only on the difference 1 - 2.
Spherical capsThe region on a sphere above a fixed line of latitude is called a spherical cap. The northern hemisphere, the region above the equator, would be a very large spherical cap. The region inside the Arctic Circle would be a smaller spherical cap.
Let R be the radius of the earth. Then the surface area above a latitude is
A = 2R^2(1 - sin ).
You could derive this using calculus by thinking of the spherical cap as a surface of revolution.
Spherical bandsGiven two latitudes 1 and 2 with latitudes 1 > 2, the area of a band between latitude 1 and latitude 2 is the area of the spherical cap above latitudes 2 minus the area of the spherical cap above latitudes 1. This gives
A = 2R^2(sin 1 - sin 2).
Area of latitude/longitude gridNow we can find the area of the region bounded by longitudes 1 and 2 and latitudes 1 and 2. The total area between latitudes 1 and 2 is given by the equation above. The subset of this area between longitudes 1 and 2 is proportional to 1 - 2. If longitude is measured in radians then
A = R^2 (sin 1 - sin 2) (1 - 2).
If longitude is measured in degrees, we have
A = R^2 (sin 1 - sin 2) (1 - 2)/180.
Related postsThe post Area of a rectangle" on a globe first appeared on John D. Cook.