Spherical geometry

Spherical geometry studies lengths and angles measured along the surface of a sphere, rather than through its interior. It is a useful model whenever travel is constrained to a round surface. Its straightest paths are great circle arcs: they are geodesics, even though they appear curved in the surrounding Euclidean space. Two distinct great circles meet at an antipodal pair, so the parallel-line picture of Euclidean geometry does not apply.

For a sphere of radius $R>0$, write two positions as $Ru,Rv$, with $u,v$ unit vectors. Their shorter surface route has great-circle distancewhere the dot product gives the central angle in radians. The straight chord through the surrounding Euclidean space instead has length $2R\sin (d/(2R))$. For a quarter-circle separation these are $\pi R/2\approx 1.571R$ and $\sqrt{2}R\approx 1.414R$: specifying which distance we mean has a measurable effect. On the unit sphere, $R=1$ and great-circle distance equals the central angle.

Consider the spherical triangle with vertices $(R,0,0)$, $(0,R,0)$ and $(0,0,R)$, joined by the shorter great circle arcs. Each side has length $\pi R/2$. At each vertex the two outgoing tangent vectors point along the other coordinate axes and are orthogonal, so all three interior angles are $\pi /2$. This is the positive coordinate octant of the sphere, not a flat triangle drawn inside it.

For a convex minor spherical triangle with interior angles $\alpha ,\beta ,\gamma$ in radians, the spherical excess formula statesOur octant therefore has area $R^2(3\pi /2-\pi )=\pi R^2/2$. Independently, eight congruent octants partition the sphere, giving $4\pi R^2/8$ again. Thus angles alone determine area when $R$ is fixed, unlike Euclidean geometry, where similar triangles can have arbitrarily different areas. The Gauss-Bonnet theorem explains the connection: the sphere's constant Gaussian curvature $1/R^2$ times the area equals this excess. For spherical triangles whose area is small compared with $R^2$, the excess is small, recovering the familiar flat picture.

A nondegenerate convex regular spherical $n$-gon contained in an open hemisphere has equal interior angles $\alpha$ in the displayed interval. The spherical excess formula gives area $n\alpha -(n-2)\pi >0$. Conversely, put equally spaced vertices at colatitude $t\in (0,\pi /2)$ around a pole. Dividing each central triangle in half gives $\cot (\alpha /2)=\cos t\tan (\pi /n)$ by the spherical law of cosines. As $t$ increases, the angle ranges continuously from the Euclidean limiting angle to $\pi$.

On a sphere of radius $R$, a spherical circle of radius $r\in (0,\pi R)$ is the locus at spherical distance $r$ from its center. Its polar angle about that center is $r/R$, so parametrization by longitude has speed $R\sin (r/R)$. Its circumference is therefore $2\pi R\sin (r/R)$. For small $r$, this is $2\pi r-\pi r^3/(3R^2)+O(r^5)$, with the deficit determined by the sphere's positive Gaussian curvature.

A spherical lune is the region between two great semicircles joining antipodal points. On a unit sphere a lune of opening angle $\alpha$ has area $2\alpha$, obtained by integrating area over longitude width $\alpha$. Lune area counting yields the spherical excess formula.

For points on a sphere of radius $R$, write their position vectors relative to the centre as $Ru,Rv$ with $u,v$ unit vectors. The great-circle distance is the shorter great circle arc length, $R$ times its central angle. Hence it equals $R\arccos (u\cdot v)$, including zero for identical points and $\pi R$ for antipodal points. At antipodal points the shortest arc is not unique, but its length is.

Spherical trigonometry calculates the sides and interior angles of spherical triangles using trigonometric functions. Its sides must be shorter great circle arcs, which are geodesics on a round sphere. A curve at constant nonzero latitude does not qualify. The spherical law of sines and spherical law of cosines play the roles of the familiar flat-triangle rules, but account for the sphere's curvature.

Use uppercase $A,B,C$ both for the vertices and for their interior angles, and lowercase $a,b,c$ for the opposite angular sides. For example, $a$ is the central angle subtended by the arc from $B$ to $C$. On a sphere of radius $R$, its physical arc length is ${\ell }_a=Ra$ when $a$ is in radians; similarly ${\ell }_b=Rb$ and ${\ell }_c=Rc$. Only on the unit sphere do the angular and physical side lengths have the same numerical values. See J. B. Tatum, Celestial Mechanics, chapter 3, section 3.5: great-circle sides and angular units.

One useful application of spherical trigonometry is finding great-circle distance from spherical coordinates. Write a unit vector as $u(\varphi ,\lambda )=(\cos \varphi \cos \lambda ,\cos \varphi \sin \lambda ,\sin \varphi )$, where $\varphi$ is latitude and $\lambda$ is the equatorial azimuth, or longitude, measured from a chosen reference direction. The dot product of two such unit vectors givesHere $\delta \in [0,\pi ]$ is the central angle and $d$ is the great-circle distance. The same identity follows from the spherical law of cosines by adding the north pole as the third vertex, with angular sides $\pi /2-{\varphi }_1$ and $\pi /2-{\varphi }_2$. This connects spherical geometry to position calculations in astronomy and to navigation on a spherical model of Earth. Choosing the surface model matters: Paula Neeley and Anthony Narkawicz, Map Projection Induced Variations in Locations of Polygon Geofence Edges, NASA/TM-2017-219675 (2017), sections 2-3 distinguishes spherical and ellipsoidal geodesics and gives this coordinate identity.

For a toy example, put two points at latitude $30^{\circ }$, with equatorial azimuths $-60^{\circ }$ and $60^{\circ }$. The great-circle distance and the arc length along constant latitude areThe great circle takes a shorter route by bending towards the pole in these spherical coordinates. To locate its highest latitude, the endpoint unit vectors are $p=(\sqrt{3}/4,-3/4,1/2)$ and $q=(\sqrt{3}/4,3/4,1/2)$. Their great circle lies in the plane through the origin with $z=2x/\sqrt{3}$. On its shorter arc, the maximum latitude occurs at $y=0$ and is $\arctan (2/\sqrt{3})\approx 49.106605^{\circ }$. Thus a geodesic can change latitude substantially even when its endpoints have equal latitude. These are lengths on an ideal sphere; the example illustrates the geometry rather than specifying an actual flight route.

Historically, spherical trigonometry was prominent in teaching as well as in astronomy. John Conway and Alex Ryba note its presence in early twentieth-century school textbooks in Remembering spherical trigonometry, The Mathematical Gazette 100(547), March 2016, pages 1-8, DOI 10.1017/mag.2016.3. Their introductory observation documents that teaching history; it does not identify a single inventor of the spherical law of sines or spherical law of cosines.

For a nondegenerate convex minor spherical triangle, with angular sides $a,b,c$ opposite interior angles $A,B,C$, the spherical law of sines statesThe sines of sides refer to central angles, so a physical side length ${\ell }_a$ on a sphere of radius $R$ must enter as $\sin ({\ell }_a/R)$. All six angles lie in $(0,\pi )$, and the denominators are positive. The reciprocal form, $\sin a/\sin A=\sin b/\sin B=\sin c/\sin C$, is equivalent. This is Tatum's sine formula, equation 3.5.1, in Celestial Mechanics, section 3.5: spherical sine and cosine formulas.

A short proof explains the symmetry of the spherical law of sines. Let $u,v,w$ be the vertex unit vectors at $A,B,C$. At $u$, the outgoing unit tangent vectors towards $v,w$ are $t_B=(v-\cos cu)/\sin c$ and $t_C=(w-\cos bu)/\sin b$. Their angle is $A$, so the magnitude of their cross product is $\sin A$. Since both tangent vectors are orthogonal to $u$, the scalar triple product obeysDividing by $\sin a\sin b\sin c$ proves the spherical law of sines. Its beauty is that the same three-dimensional scalar triple product can be calculated at any vertex, producing all three ratios at once.

For the spherical triangle worked out under the spherical law of cosines, $A=90^{\circ }$, $b=c=60^{\circ }$ and $\cos a=1/4$. Hence $\sin a=\sqrt{15}/4$, and the spherical law of sines givesThe inverse sine gives $B\approx 63.434949^{\circ }$, but the supplementary angle $116.565051^{\circ }$ has the same sine. The spherical law of sines alone therefore leaves two candidates for this angle. The spherical law of cosines, using the known three sides, gives $\cos B=1/\sqrt{5}>0$ and selects $63.434949^{\circ }$; likewise $C=B$. When only two sides and an opposite angle are given, two distinct valid spherical triangles can remain. A cosine relation resolves an inverse sine branch only when enough additional data are known; it does not make every ambiguous data set unique.

For a nondegenerate convex minor spherical triangle, the spherical law of cosines relates two angular sides $b,c$, their included interior angle $A$, and the opposite angular side $a$:Relabelling the vertices gives the corresponding formulas for $b,B$ and $c,C$. As with the spherical law of sines, physical side lengths must first be divided by the sphere's radius. The spherical law of cosines determines $a\in (0,\pi )$ uniquely from $b,c,A$ when the data describe such a spherical triangle. It also recovers each interior angle from three sides, for exampleUsing the cosine and the range $0<A<\pi$ selects a unique angle, avoiding the supplementary-angle ambiguity of an inverse sine calculation.

To derive the spherical law of cosines, let $u,v,w$ be the vertex unit vectors at $A,B,C$, so $u\cdot v=\cos c$, $u\cdot w=\cos b$ and $v\cdot w=\cos a$. The outgoing unit tangent vectors at $u$ are $t_B=(v-\cos cu)/\sin c$ and $t_C=(w-\cos bu)/\sin b$. Their dot product is both $\cos A$ and $(\cos a-\cos b\cos c)/(\sin b\sin c)$, proving the identity. This gives a geometric explanation: spherical trigonometry compares angles in the tangent plane using vectors in the surrounding Euclidean space.

For a worked example on the unit sphere, take the verticesAt $A$, the outgoing unit tangent vectors are $(1,0,0)$ and $(0,1,0)$, so the interior angle is a right angle. The vertex dot products give $b=c=\pi /3=60^{\circ }$. The spherical law of cosines then reduces to the right-angle identity $\cos a=\cos b\cos c$, givingTo find $B$, use the relabelled spherical law of cosines:The physical side lengths on a sphere of radius $R$ are ${\ell }_a\approx 1.318116R$ and ${\ell }_b={\ell }_c=\pi R/3\approx 1.047198R$. For example, on a toy sphere with $R=10$ metres they are approximately $13.18116$, $10.47198$ and $10.47198$ metres; the interior angles stay the same. Applying the planar Pythagorean theorem directly to the angular sides would instead give $\sqrt{60^2+60^2}\approx 84.852814^{\circ }$. That flat right triangle calculation does not give the great-circle distance between $B$ and $C$.

The familiar law of cosines reappears as a local approximation. With $a,b,c$ all small in radians, substitute the small-angle approximations $\sin x\approx x$ and $\cos x\approx 1-x^2/2$ into the spherical law of cosines and keep terms through second order:Multiplication by $R^2$ gives the planar law of cosines for the physical side lengths. For a fixed small surface triangle on a progressively larger sphere, the angular sides shrink and this approximation improves. Thus the curved and flat formulas fit together: the Pythagorean theorem is recovered locally when $A$ is a right angle.

A spherical triangle is bounded by three great circle arcs. For the usual convex minor spherical triangle, the vertices are neither antipodal nor on one great circle, and its angular sides and interior angles lie in $(0,\pi )$. An interior angle is measured between outgoing tangent vectors to the two sides at a vertex, rather than between straight chords inside the sphere. The spherical law of sines and spherical law of cosines determine its side-angle relations; the spherical excess formula determines its area.

Reflections in the sides of a spherical triangle of angles $\pi /p,\pi /q,\pi /r$ generate a finite group when its images tessellate the sphere. The triangle area is $\pi (1/p+1/q+1/r-1)$ by the spherical excess formula. Label its three vertex types and continue those labels across each reflected side; they are consistent around every vertex and, because the sphere is simply connected, globally. The generated group preserves these labels, so a chamber stabilizer fixes all three vertices and is trivial. Dividing sphere area by chamber area gives its order. Up to permutation the possibilities are $(2,2,n)$ for $n\ge 2$, and $(2,3,3),(2,3,4),(2,3,5)$, of orders $4n,24,48,120$.