Scinovex

Research article · 2022 · Full text

International Journal of Statistics and Applied Mathematics

Great circle theorem and the application of the spherical cosine rule to estimate distances on a globe

Maharshi Roy

Volume 7(3) · Pages 136–142 · DOI: 10.22271/maths.2022.v7.i3b.833 ↗

Research areasHistorical Geography and CartographyMathematics and ApplicationsData Management and Algorithms

↓ Read full-text PDF3 citations

Abstract

Author-written

This research paper explores my interest in the Great Circle Theorem and Non-Euclidean Geometry to apply its theory in Earth Sciences and explore navigating flight distances across the globe. I use three different equations, based on the situation of the two points on the map and their relation with each other, to obtain approximate great circle distances between the two points. The values acquired from the final equation, which is derived from the Spherical Law of Cosines, is compared with the actual great circle distances between the two points on the globe. A percentage error is calculated and the average is found, after which a correlation is drawn between the percentage error and the magnitude of the value of the great circle distance.

Summary

Source-checked

Plain-language reading of the full paper, prepared with AI and checked against the text. Each point links to the authors' own words.

  1. Why it matters

    • Because the Earth is a sphere, the shortest distance between two points on its surface is an arc rather than a straight line. View evidence ↗

      Why it matters

      Because the Earth is a sphere, the shortest distance between two points on its surface is an arc rather than a straight line.

      In the paper's own words

      “I realised that this was due to the fact that the Earth is a sphere and therefore the closest distance between two points on a curved surface would result in an arch rather than a straight line [1].”

      Page 1

      Found in: Introduction

      Read it in the full textOpen PDF at page 1 ↗
    • Aircraft utilize the shortest possible three-dimensional route, referred to as a geodesic or great circle route. View evidence ↗

      Why it matters

      Aircraft utilize the shortest possible three-dimensional route, referred to as a geodesic or great circle route.

      In the paper's own words

      “planes would travel through the shortest three-dimensional route possible, known as the “geodesic” or “great circle route”.”

      Page 1

      Found in: Introduction

      Read it in the full textOpen PDF at page 1 ↗
    • Projecting a 3D spherical globe onto a 2D map, such as a Mercator projection, causes polar regions to appear disproportionately larger. View evidence ↗

      Why it matters

      Projecting a 3D spherical globe onto a 2D map, such as a Mercator projection, causes polar regions to appear disproportionately larger.

      In the paper's own words

      “due to the fact that the globe is a spherical 3-dimensional object, when projecting the earth onto a 2-D map, areas to the northmost and southmost of the globe look disproportionately larger than they actually are.”

      Page 3

      Found in: Introduction

      Read it in the full textOpen PDF at page 3 ↗
  2. The question

    • The study aims to apply the Great Circle Theorem and Non-Euclidean Geometry to estimate and explore flight distances across the globe. View evidence ↗

      The question

      The study aims to apply the Great Circle Theorem and Non-Euclidean Geometry to estimate and explore flight distances across the globe.

      In the paper's own words

      “This research paper explores my interest in the Great Circle Theorem and Non-Euclidean Geometry to apply its theory in Earth Sciences and explore navigating flight distances across the globe.”

      Page 1

      Found in: Introduction

      Read it in the full textOpen PDF at page 1 ↗
  3. What they did

    • The author uses three different equations to calculate approximate distances, depending on the relative positions of the two points on the map. View evidence ↗

      What they did

      The author uses three different equations to calculate approximate distances, depending on the relative positions of the two points on the map.

      In the paper's own words

      “I use three different equations, based on the situation of the two points on the map and their relation with each other, to obtain approximate great circle distances between the two points.”

      Page 1

      Found in: Another major difference between Euclidean Geometry

      Read it in the full textOpen PDF at page 1 ↗
    • The Spherical Law of Cosines is used for distance calculations due to its relative simplicity and accuracy compared to other methods like the haversine or Vincenty formulas. View evidence ↗

      What they did

      The Spherical Law of Cosines is used for distance calculations due to its relative simplicity and accuracy compared to other methods like the haversine or Vincenty formulas.

      In the paper's own words

      “Most geoscientists prefer to use the “Spherical Law of Cosines” equation for its relative simplicity and accuracy.”

      Page 6

      Found in: Inserting the values for r and θ:

      Read it in the full textOpen PDF at page 6 ↗
    • The author calculated the distance between 10 pairs of random points on the globe to determine the percentage error of the Spherical Law of Cosines. View evidence ↗

      What they did

      The author calculated the distance between 10 pairs of random points on the globe to determine the percentage error of the Spherical Law of Cosines.

      In the paper's own words

      “I will be calculating the distance between random points on the globe and calculate the percentage error. I will do this for 10 great circle distances”

      Page 6

      Found in: Inserting the values for r and θ:

      Read it in the full textOpen PDF at page 6 ↗
  4. What they found

    • For the route between Madrid and JFK airports, the calculated distance was 15 km different from the actual distance of approximately 5775 km, resulting in a 0.26% error. View evidence ↗

      What they found

      For the route between Madrid and JFK airports, the calculated distance was 15 km different from the actual distance of approximately 5775 km, resulting in a 0.26% error.

      In the paper's own words

      “I found out the actual great circle distance between Madrid Airport and JFK Airport to be approx. 5775 km. My calculations were just 15km off with approximately 0.26% error.”

      Page 6

      Found in: Inserting the values for r and θ:

      Read it in the full textOpen PDF at page 6 ↗
    • The average percentage error across the calculated distances was 0.96%. View evidence ↗

      What they found

      The average percentage error across the calculated distances was 0.96%.

      In the paper's own words

      “Through my calculations, I found out that the average percentage error was 0.96%.”

      Page 7

      Found in: Inserting the values for r and θ:

      Read it in the full textOpen PDF at page 7 ↗
    • A trend was observed where the percentage error decreased as the distance being calculated increased. View evidence ↗

      What they found

      A trend was observed where the percentage error decreased as the distance being calculated increased.

      In the paper's own words

      “It was noticed that, as the distance calculated increased, the percentage error fell.”

      Page 7

      Found in: Inserting the values for r and θ:

      Read it in the full textOpen PDF at page 7 ↗
    • The error was higher (2.50%) for shorter distances, such as Delhi to Jakarta, and lower (0.04%) for longer distances, such as Doha to New York. View evidence ↗

      What they found

      The error was higher (2.50%) for shorter distances, such as Delhi to Jakarta, and lower (0.04%) for longer distances, such as Doha to New York.

      In the paper's own words

      “Delhi and Jakarta are close by, so the percentage error was greater, at 2.50%, whilst Doha and New York are quite far away, therefore resulting in a smaller percentage error of 0.04%.”

      Page 7

      Found in: Inserting the values for r and θ:

      Read it in the full textOpen PDF at page 7 ↗
  5. What it means

    • The author concludes that the Spherical Law of Cosines is more accurate when calculating larger distances across the globe. View evidence ↗

      What it means

      The author concludes that the Spherical Law of Cosines is more accurate when calculating larger distances across the globe.

      In the paper's own words

      “Through my research, I concluded that the Spherical Law of Cosines achieves a higher level of accuracy while being used to calculate larger distances across the Globe.”

      Page 7

      Found in: Conclusion

      Read it in the full textOpen PDF at page 7 ↗

Visual summary

1 visuals

Illustrations explain the paper; they do not show study data.

AI-generated illustration
AI-generated conceptual illustration (a terrestrial globe). It explains the idea and does not show study data.↗ Enlarge
Great Circle vs Straight Line
AI-generated conceptual illustration (a terrestrial globe). It explains the idea and does not show study data.

AI-generated conceptual illustration (a terrestrial globe). It explains the idea and does not show study data.

Open full-size image ↗

Great Circle vs Straight Line

AI-generated conceptual illustration (a terrestrial globe). It explains the idea and does not show study data.

View source evidence ↗

Figures & tables

From the original paper
Flight path between Madrid and New York on 3D diagramFigure 1Page 2Flight path on a 2-D map between Madrid and New YorkFigure 2Page 2Great Circle Routes stretched to infinityFigure 4Page 3A line segmentFigure 3Page 3Longitude and Latitude LinesFigure 5Page 4Angles of a Triangle on a SphereFigure 6Page 4Percentage Error CalculationsTable 1Page 7
FULL TEXT Original article · reflowed for reading

Introduction

When I was a child, my father used to work for Etihad Airways. We as a family used to go on many trips and therefore, oftentimes I would find myself in a seat on a plane. On the plane, I would see both two-dimensional and three-dimensional maps of the globe and a small figure of a flight travelling from one point to another in an arc across the globe. Along with the displayed maps, other various information was also given, such as a compass, airspeed, ground speed, etc. One day when travelling from Madrid to New York, I noticed that, on the 2D map, the flight was travelling in a curved path, but on the 3D map, the flight was travelling in a straight line between the two points of Madrid and New York. I realised that this was due to the fact that the Earth is a sphere and therefore the closest distance between two points on a curved surface would result in an arch rather than a straight line [1] . This greatly intrigued me, and given my inquisitive nature as well as my interest in Earth Sciences and Mathematics, I set out to understand this phenomenon. I asked my father numerous questions and I found out that planes would travel through the shortest threedimensional route possible, known as the "geodesic" or "great circle route". I learned that when converting a flight path from a 3-D plane to a 2-D plane, it may cause the line to look more like an arc on the flat surface. So even though in actuality, the flight may be moving in a straight line, on a 2-D surface, this would result in an arc, and as I hypothesised earlier, this is due to the fact that the Earth is really a sphere, rather than a flat surface. In Euclidean geometry, this phenomenon is known as the Great Circle Theorem [2] .

Metre: Standard metric for measuring distances in geometry.

It is interesting to note that its initial definition is defined as "one ten-millionth of the distance from the North Pole to the Equator travelling along the meridian through Paris."

Great Circle Route: A circle on the surface of a sphere in a plane that passes through the centre of the sphere. The great circle is a perfect route for ships and planes because it represents the shortest distance between any two points on the sphere [3] . Small Circles: Any latitude meridian that is not passing through the largest circumference of the sphere (the equator) is called a small circle. These circles have radii smaller than that of the Earth, i.e less than approx 6400km. Any place situated in the same latitude is due exactly east or west of each other, since all small circles are parallel to each other.

Introduction to Spherical Geometry

After going through numerous research papers, I quickly came to understand that Euclidean Geometry does not apply to spherical objects. Therefore, the Spherical Geometry and the Great Circle Theorem came about. A great circle route is the straight line equivalent of spherical geometry. It represents the shortest distance between any two points on a sphere. Whilst Euclidean Geometry uses two dimensional planes to plot points and lines on, Spherical Geometry requires a sphere to plot lines, points etc. Spherical Geometry has numerous basic rules that govern the subject.

These rules are 1. A point in Spherical Geometry is defined as having no parts and lies on the surface of a sphere. 2. The shortest distance between two points on the surface of a sphere is known as "Great Circles" 3. Great Circles can be extended to circumnavigate the sphere entirely and are always finitely the length of the circumference of the sphere 4. In Spherical Geometry, the plane is the sphere itself 5. Angles are usually measured in radians and regular Euclidean Geometry laws don't apply to shapes lying on the surface of a sphere. (Examples mentioned below) 6. All right angles on the surface of a sphere are congruent 7. When lines intersect, they always do so at two points 8. Perpendicular lines can be formed by 2 intersecting great circles and form exactly 8 right angles 9. Parallel lines do not exist on a sphere, any two great circles will always intersect at two places on the surface of the sphere 10. Any great circle has the same radius of the sphere and will divide a sphere into two equal halves

As stated above, Spherical Geometry is an extreme tangent from the two-dimensional Euclidean Geometry that I was taught in my junior classes. Spherical geometry tends to reject many of the basic premises of Euclidean Geometry. Therefore, it is important to analyse and understand some of these major differences to gain a better understanding of the theory. Some major differences between Euclidean Geometry and Spherical Geometry are:

1. In Euclidean Geometry, the shortest path between any two points on a plane is known as a line segment. Line segments have a constant gradient, meaning that the rate of change of its slope is zero and the ratio between the vertical displacement and horizontal displacement remains constant. An example of a line segment between two points (A and B) can be seen below. 2. Although, on a 2-dimensional plane, a line segment is the shortest path between two points, on any curved surface, such as the Earth, Euclidean Geometry rules don't apply. If we were to pick any two points on the surface of a sphere and draw a straight line between those two points, we would end up with a great circle route. Euclidean Geometry also states that any finite line connecting two points on a plane can be extended to infinity in both directions and will never meet. However this is not the case in Spherical Geometry. If a line connecting two points on a sphere is stretched out, it would travel along the sphere and reach its original starting point in a loop. If stretched to infinity, the line will wrap around the sphere an infinite amount of times over. As seen in figure 4 above, the great circle joining points P to Q is stretched to infinity. As a result, the line travels around the sphere reaching point v and u on the equator of the sphere and then circumnavigating back to point P. As a result, a circle has been formed across the sphere, with a radius of the sphere itself and a circumference of the equator of the sphere.

Figure 4. Great Circle Routes stretched to infinity
Great Circle Routes stretched to infinity↗ Enlarge
Figure 4
Great Circle Routes stretched to infinity

Great Circle Routes stretched to infinity

Open full-size image ↗
Original paper figureSource PDF · page 3 ↗

Another major difference between Euclidean Geometry

and Spherical Geometry is regarding the statement that "parallel lines never intersect". Although this statement stands true in Euclidean Geometry on a two-dimensional plane, it doesn't apply for a three dimensional sphere.

There are no parallel lines on a sphere. Any two great circles drawn on the surface of a sphere will always intersect at two places on the sphere. On Earth, although at the equator all latitude lines are parallel to each other at 90 degrees, they eventually converge at the poles of the planet. The same can be seen in figure 5 below. 4. The way that angles are measured in Spherical Geometry is also vastly different from that of Euclidean Geometry. For example, In Spherical Geometry, angles are calculated between 2 or more great circles, resulting in differences in theory when compared to Euclidean trigonometry. In Euclidean Geometry, we are taught that all angles in a triangle must sum up to 180 degrees. However, this does not stand true in the case of Spherical Geometry. On a spherical plane, the sum of all interior angles of a triangle will always exceed 180 degrees. Formula Derivation and Calculations Now that I have a considerable amount of knowledge and understanding of the key concepts regarding Spherical Geometry and the Great Circle Theorem, I may now apply them in calculating distances across a sphere. There are three ways in which this can be done, depending on the location of the two points between which distance must be calculated [4] .

Figure 5. Longitude and Latitude Lines
Longitude and Latitude Lines↗ Enlarge
Figure 5
Longitude and Latitude Lines

Longitude and Latitude Lines

Open full-size image ↗
Original paper figureSource PDF · page 4 ↗

The three situations are  Both points on the globe lie on the equator or same meridian of longitude  Both points lie on parallel latitudes  Points lie on random parts of the globe For the first situation, the calculation is straightforward. Since both points lie on the same meridian of longitude, or on the equator, all we have to do is calculate the angle between the two points on the globe in relation to each other. If we assume that the equator is the circumference of a circle where the radius is the same as that of the Earth, then we can easily find the distance between two points on the circumference of that circle by deriving its formula from the basics of circle geometry [5] . For this situation, let's take two points that lie on the equator. As we already know, the Equator is a great circle with radius 6400 km, we can easily calculate the distance between any two points on the equator using the same equation of the length of an arc. A city that lies on the equator is Libreville in Gabon. Singapore lies on the equator as well. To find the distance between the two cities along the equator, we first require their exact location on the globe. Libreville is 9.27•E and Singapore is 103.81•E First, we must calculate the difference in angle between the two points [6] .

Let "l" be the length of the arc: Paris lies at approx 48.67• N and 2.33• E. We already know that each meridian is a great circle of approx 6400km radius. The angle between the Equator and the Latitude of Paris is already given to be 48.67•. To find the angle between the Latitude of Paris and the North Pole, we must subtract the angle between the latitude of Paris and the Equator from 90•.

Using the equation for the length of an arc of a circle with a given radius, we can calculate the distance between the North Pole and Paris.

Let "l" be the length of the arc Inserting the values for r and θ:

We can use the same method to estimate the distance between Paris and the Equator Let us use another example for the same situation. Melbourne is at 37.82• S and 144.97• E. Let's try to calculate the distance between Melbourne and the North Pole, The South Pole and the Equator, assuming that we are travelling along the meridian line. The angle between the Equator and Melbourne is already given to be 37.82•. We already know that each meridian line is a great circle with the radius of Earth, approx 6400km.

Let "l" be the length of the arc:

Inserting the values for r and θ:

To find the angle between Melbourne and the North Pole, we must add 90•. The reason for adding 90• instead of subtracting is due to the fact that Melbourne is located in the Southern Hemisphere.

Let "l" be the length of the arc:

To find the angle between Melbourne and the South Pole, we must subtract 90•.

Let "l" be the length of the arc: Now, let's try to determine the distance between two points on the globe, which share the same latitude. Sydney City and Margaret River both lie on the 33.5•S parallel latitude. Sydney lies on 151.13•E and Margaret River lies on 115.04•E. To calculate the distance between these two points, we would have to calculate the length of the small circle travelling along the 33.5•S latitude. First we require to calculate the radius of the small circle connecting the two points:

The longitudinal angle difference between the two points is:

Let "l" be the length of the arc: Finally, let's find the distance between any two points on the globe, using only their coordinates. To do this, there are three different techniques, each of varying difficulties. The easiest method is by using the "Spherical Law of Cosines". This method is often used to calculate distances on a perfect sphere, therefore, it tends to be slightly inaccurate while calculating distances on the globe, since the Earth is slightly ellipsoid in shape. The second method, the most commonly used by geoscientists, is the "haversine" calculation method. It is slightly tedious and taxing to calculate using pen and paper, so scientists use a computer to insert the values and calculate the distance. It is usually precise and can be used to find distances between two points on a sphere fairly accurately, however, it has issues finding any two points at each end of a diameter, known as "antipodal points". The third and most complex method is finding the distance between two points on an ellipsoid using the Vincenty Formula. In this case, the sphere is an ellipsoid with equal major and minor axes. However, this is a very complex formula and is not frequently used to calculate distances on the globe. Most geoscientists prefer to use the "Spherical Law of Cosines" equation for its relative simplicity and accuracy.

The Spherical Law of Cosines is given by the equation: I calculated the percentage error using this formula:

To further my research and for my own practice, I will be calculating the distance between random points on the globe and calculate the percentage error. I will do this for 10 great circle distances and create a spreadsheet displaying my values and the percentage errors. Through my calculations, I found out that the average percentage error was 0.96%. I also noticed a trend with the percentage errors. It was noticed that, as the distance calculated increased, the percentage error fell. This trend can be best noticed by the difference in percentage error values of the great circle distance calculated between Delhi and Jakarta and the great circle distance calculated between Doha and New York. Delhi and Jakarta are close by, so the percentage error was greater, at 2.50%, whilst Doha and New York are quite far away, therefore resulting in a smaller percentage error of 0.04%.

Conclusion

Through my research, I learned a variety of new concepts. Initially, I started off by understanding the key concepts of Non-Euclidean Geometry, and its comparisons and contrasts between Euclidean Geometry. Next, I learned about the Great Circle Theorem and its applications in the real world. I then applied these concepts to calculate the distance between two points on Earth, using only location coordinates. I further applied my knowledge to find the percentage error that I incurred while using the Spherical Law of Cosines. Through my research, I concluded that the Spherical Law of Cosines achieves a higher level of accuracy while being used to calculate larger distances across the Globe.

Figures and tables from the paper

Figure 1. Flight path between Madrid and New York on 3D diagram
Flight path between Madrid and New York on 3D diagram↗ Enlarge
Figure 1
Flight path between Madrid and New York on 3D diagram

Flight path between Madrid and New York on 3D diagram

Open full-size image ↗
Original paper figureSource PDF · page 2 ↗
Figure 2. Flight path on a 2-D map between Madrid and New York
Flight path on a 2-D map between Madrid and New York↗ Enlarge
Figure 2
Flight path on a 2-D map between Madrid and New York

Flight path on a 2-D map between Madrid and New York

Open full-size image ↗
Original paper figureSource PDF · page 2 ↗
Figure 3. A line segment
A line segment↗ Enlarge
Figure 3
A line segment

A line segment

Open full-size image ↗
Original paper figureSource PDF · page 3 ↗
Figure 6. Angles of a Triangle on a Sphere
Angles of a Triangle on a Sphere↗ Enlarge
Figure 6
Angles of a Triangle on a Sphere

Angles of a Triangle on a Sphere

Open full-size image ↗
Original paper figureSource PDF · page 4 ↗
Table 1. Percentage Error Calculations
Percentage Error Calculations↗ Enlarge
Table 1
Percentage Error Calculations

Percentage Error Calculations

Open full-size image ↗
Original paper figureSource PDF · page 7 ↗

References

6 in the original paper
  1. Porcu E, Bevilacqua M, Genton MG. Spatio-temporal covariance and cross-covariance functions of the great circle distance on a sphere. Journal of the American Statistical Association. 2016;111(514):888-898.
  2. Gluck H, Warner FW. Great circle fibrations of the threesphere. Duke Mathematical Journal. 1983;50(1):107-132.
  3. Guinness J, Fuentes M. Isotropic covariance functions on spheres: Some properties and modeling considerations. Journal of Multivariate Analysis. 2016;143:143-52.
  4. Gneiting T. Strictly and non-strictly positive definite functions on spheres. Bernoulli. 2013;19(4):1327-1349.
  5. Buss SR, Fillmore JP. Spherical averages and applications to spherical splines and interpolation. ACM Transactions on Graphics (TOG). 2001;20(2):95-126.
  6. Helgason S. The totally-geodesic Radon transform on constant curvature spaces. Contemp. Math. 1990;113:141-149.

Cited by 3

From OpenAlex, OpenCitations and reference lists read by Scinovex