Celestial_Navigation_Tutorial
Celestial_Navigation_Tutorial
Position Lines 3
Corrections 10
Index Error
Dip
Refraction
Temperature and Pressure Corrections to Refraction
Semi Diameter
Augmentation of the Moon’s Semi-Diameter
Parallax
Reduction of the Moon’s Horizontal Parallax
Examples
Using the altitude of a celestial body is similar to using the altitude of a lighthouse or
similar object of known height, to obtain a distance.
One object or body provides a distance but the observer can be anywhere on a circle
of that radius away from the object. At least two distances/ circles are necessary for a
position. (Three avoids ambiguity.)
In practice, only that part of the circle near an assumed position would be drawn.
Using a Sextant for Celestial Navigation
After a few corrections, a sextant gives the true distance of a body if measured on an
imaginary sphere surrounding the earth.
Using a Nautical Almanac to find the position of the body, the body’s position could
be plotted on an appropriate chart and then a circle of the correct radius drawn around
it. In practice the circles are usually thousands of miles in radius therefore distances
are calculated and compared with an estimate.
Working on this sphere, the distance becomes [90° - Altitude.] The point on the
sphere corresponding to the Observer is known as his Zenith.
Position Lines
Because of the immense radius, the short length of interest can be considered a
straight line.
Comparing the observed distance to the body and the calculated distance from an
assumed position provides the distance towards or away from the body. The bearing
is found by calculation and then it can be a simple matter of marking the assumed
position on a chart, drawing a line in the direction of the body and marking off the
difference between the observed and calculated distances.
The observed distance is known as the True Zenith Distance (TZD.) The value
based on the assumed position is the Calculated Zenith Distance (CZD.) The
difference between the two is known as the Intercept.
The closest point on this circle is known as the Intercept Terminal Position (ITP)
and the line representing the circle at that position is called a Position Line.
Additional sights provide additional position lines that intersect to provide a Fix.
A Running Fix
A vessel is usually moving between sights therefore they are combined "on the run."
The position line from a first sight must be moved to allow it to be combined with
another position line for a different time.
After a second sight has been calculated, its position line can be plotted and combined
with the first to provide a fix.
The shorter the “run” between the two observations, the more accurate will be the
final result. This is why star sights are usually better than sun sights. The opportunity
to obtain more position lines of different directions is an additional benefit.
Notes on Running Fixes
Under normal conditions, one would expect an error of +/- 0’.3 in the Position Lines.
(This error is mainly due to the time recorded under practical conditions.) Land
Surveyors achieve accuracy comparable to a GPS using more sophisticated
instruments but the same calculations/ method.
Final accuracy is obviously improved by taking more observations. Six star sights will
typically provide a fix within 0’.2 of the true position.
Most people adopt some shortcuts in the interest of speed. These have a cost in terms
of accuracy. The Sun's Total Correction Tables assume that the Sun's semi-diameter is
either 15'.9 or 16'.2. A Sun Sight in April (SD = 16'.0) is immediately in error by 0'.2.
Tables are rounded to the nearest 0'.1 which could introduce a cumulative error of
0'05 for every item.
With Star sights, the short interval between the first and last sight means that many
people use a single position for all the sights and plot the results without allowing for
the vessel's movement. The error is larger than above, but quite acceptable in mid-
ocean.
The method used until the 1980s was the Haversine Formula and Log Tables. A few
commercial navigators used Sight Reduction Tables but most preferred the longer
method in the interests of accuracy and flexibility.
Correct the Sextant Altitude to find the true distance of the body
Calculate the bearing and distance from an assumed position
Using the difference in distances to obtain a Position Line
Altitudes of Stars do not need a Semi-Diameter correction while the Moon needs
more corrections. See examples at the end of the next section.
Positions for the observer and position lines can be plotted on a chart or calculated.
The section on “Sailings” deals with mathematical calculations.
The other terms in the following formulae are derived from a Nautical Almanac. (See
Nautical Almanac Information.)
The formulae for calculating the distance of the body and its altitude are
ABC tables are very easy to use and more than adequate for the bearing of a celestial
body.
These tables avoid the need to use a calculator or Log tables but are based on the
previous formulae.
A = Tan(Lat) / Tan(LHA)
B = Tan(Dec) / Sin(LHA)
C = Difference A ~ B = 1/ [Tan(Azimuth) x Cos(Lat) ]
The C Table gives a bearing of 22.0. The sign of C means that this bearing is south.
It is west because the LHA is less than 180.
The effect of rounding ABC Tables’ values is negligible (+/- 01.) This is not true of
the older Sight Reduction Tables where the calculated altitude is rounded to the
nearest minute. Furthermore the need to use a plotting sheet with a rounded, estimated
position provides considerable scope for inaccuracy. (Sight Reduction Tables were
known as the Air Navigation Tables until 2003.)
The author’s preferred manual method is to use a calculator for the Zenith Distance
and ABC tables for Azimuths. Without a calculator he would still use the Cosine
formula but with log tables.
Obtaining a Position Line
The difference between the True (TZD) and Calculated (CZD) Zenith Distances is the
Intercept.
If the TZD is less than the CZD then the assumed position must be moved in the
direction of the body by the amount of the Intercept. This gives a position of the
correct distance from the body. It is known as the Intercept Terminal Position or ITP.
As the radius of the circle is normally very large, it is considered to be a straight line
near this point. A line at 90 to the direction of the body is the Position Line.
A single Position Line must be combined with other observations for a fix. This can
be achieved using a plotting sheet and then transferring the ITP by the distance to the
next sight and redrawing the Transferred Position Line in the same direction as the
original.
For Sun sights, it is more usual to calculate the ITP of a morning sight and then
calculate the transferred position for the Sun's Meridian Passage (Noon.) The
difference between calculated and observed latitudes provides a longitude using
“Plane Sailing.” With a little practice, this will be found to be a faster, not to mention
more accurate method.
For Star Sights, many people use a single position and then plot the Position Lines
without allowing for the vessel's movement. This may appear a sloppy practice but a
few miles error mid-ocean is usually irrelevant. Even if the position at sunset was
perfect, there is no guarantee that the position an hour later is within a mile. Indeed
even if the position agrees perfectly with a GPS position, there is no guarantee that an
intervening military operation has not thrown the GPS position out let alone a fault in
the equipment/ aerial. “I am about here,” is a far safer assumption than “My
wheelhouse is/ was within 10m of this position.”
Corrections to a Sextant Altitude
Index Error
This error can be found using the horizon. The sextant’s altitude is set to zero and
then the two images of the horizon are aligned. The Index Error can then be read off.
If the sextant altitude reads high, the correction is subtractive and termed “On the
Arc.” “Off the Arc” is the opposite.
After Index Error has been applied, the Sextant Altitude it is referred to as the
Observed Altitude.
Refraction
The deflection of light as it enters/ passes through the atmosphere is known as
Refraction.
Refraction is stable and therefore predictable above about 15°, below that one needs
to consider the characteristics of the atmospheric layers through which the light
passes at that time. (Taking the altitude of bodies at less than 15° is usually avoided
for this reason.)
Refraction tables make assumptions on the layers for low altitudes and should be
treated with caution. +/- 2 is not uncommon at an altitude of 2.
Refraction is subtracted from the Apparent Altitude to obtain the True Altitude.
Temperature and Pressure Correction for Refraction
The correction for Refraction assumes a temperature of 10° C and pressure of
1010mb. This may be modified for actual temperature and pressure. A temperature
difference of 10° C will alter Refraction by 3% and a 10mb pressure difference will
change Refraction by 1%. (0’05 and 0’.02 for an altitude of 30°)
Semi-Diameter
When measuring the altitudes of the Sun, Moon, Venus and Mars, it is usual to use
either the top (Upper Limb) or bottom (Lower Limb) of the body. This offset must
then be removed before comparison with the calculated value.
The angular diameter of a body depends on its distance from the Earth. Thus for the
Sun the Semi-Diameter varies between 16’.3 in January, when the Sun is closest and
15’.7 in June when it is furthest away.
For a lower limb observation, the Semi-Diameter should be added to the True altitude.
Horizontal Parallax is used in the formula as the lunar distance is not provided in a
Nautical Almanac.
This correction is typically 0’.15 and should be added to the Moon’s Semi-Diameter
before applying the Semi-Diameter to the True Altitude.
Parallax in Altitude
The Parallax correction allows for the difference in the altitude measured on the
Earth’s surface versus the altitude that would be measured at the centre of the Earth.
The effect of parallax reduces with altitude. It is greatest when the body is at the
horizon (Horizontal Parallax) and declines to zero when the body is overhead.
The effect is also proportional to the distance of the body. Thus the Horizontal
Parallax for the Moon is about 1° but only 0’.15 for the Sun. This correction must be
included for the Moon but is usually ignored for the Sun. It can be significant for
Venus and Mars, depending on their distance, but is always insignificant for Jupiter
and Saturn. (< 0’.05)
Add or Subtract
For a Star
The GHA for the Sun, Moon and Planets is tabulated in a Nautical Almanac for each
hour.
An increment is applied to allow for the minutes and seconds. These assume that the
GHAs change by a uniform amount per hour. This assumption is corrected using a
“v” correction.
To find the Increment, the Increments and Corrections pages are entered with the
minutes and seconds of the sight. For a Sun sight taken at exactly 30 minutes past the
hour, the increment value would be 7° 30’.0,
On the left hand side of the Increments page for the appropriate minutes, there is the
correction to apply based on the tabulated “v” value for the body. The “v” corrections
are considered to be linear but are actually tabulated for halfway through the minute.
This is why a “v” of 12.0 produces 6.1 for 30 minutes. The “v” adjustment is always
positive for the Moon but can be negative for the planets.
The “v” correction for the planets is found at the bottom of the daily pages. The rapid
changes in the motion of the Moon mean that a “v” value is tabulated for each hour.
No “v” correction is supplied for the Sun. Instead the GHA values for each hour in the
tables are massaged.
The GHA for stars are treated differently. Here the SHA (Sidereal Hour Angle) or
angular distance west of Aries is tabulated for each three day period. The GHA of
Aries is then tabulated for each hour and the Increment value for minutes and seconds
is found in exactly the same way as for the Sun, Moon and planets. No “v” value is
given for simplicity. The GHA for a star is the sum of the GHA for Aries and the
SHA.
Declination (Dec)
The Declination for most bodies is tabulated for each hour. Due to the very small
movement of stars, it is only provided once for each three day period.
The hourly change of Declination is usually small therefore the adjustment is found
using the “d” value at the bottom of the page for the Sun and planets. (It is negligible
for stars.) The Moon however moves rapidly which means that its’ “d” value is
provided for each hour.
The “d” correction is found in the same manner as for “v” by going to the appropriate
Increments and Corrections page to obtain the appropriate correction from the right
hand “v and d” section.
The direction in which to apply the “d” correction (North or South) is determined by
examining the next hourly value.
It goes on to say; “In practice it may be expected that only one third of the values of
GHA and Dec taken out will have errors larger than 0’.05 and less than one tenth will
have errors larger than 0’.1.”
The superficial attraction of using data from a source such as the Astronomical
Almanac in the interests of accuracy is illusory. RA and Declination may be quoted in
arc-seconds to several decimal places but the bodies, particularly the Moon, do not
necessarily move in a linear manner during the day. Similarly the Hour Angle of
Aries does not change linearly. Irregular motions are included in Nautical Almanac
data.
Example of Calculations for the Moon ~ 00:30:05 GMT, 19th January 2011
The Declination value of 40° 15’.8 taken from the day’s pages is not adjusted.
Sailings
Tabulated values are found in Nautical Tables as “Traverse Tables.” However the
simplicity of the formulae are ideal for calculators.
Example
Initial Position 45° N 30° W
Course 045° T
Distance 100’
For a Noon calculation, dLat is the difference between the calculated and observed
latitudes. The "Course" becomes the direction of the Position Line.
Mercator Sailing.
Mercator Sailing allow for the oblate shape of the Earth. (A squashed sphere.)
DMP stands for Difference in Meridian Parts. Meridian Parts are the distance in
nautical miles from the equator to the required latitude. These are tabulated in
Nautical Tables.
Example
Initial Position 45° N 30° W
Course 045° T
Distance 100’
Meridian Parts
for 45° 3013.38
for 46° 10’.7 3114.08
Difference in Meridian Parts 100.7
The drawback with Mercator Sailing is the need to refer to tables. If you use a
programmable calculator then this is the formula to calculate Meridian Parts. (The
infinite series of terms in Bowditch et al is simply an expansion of this.)
Note that most nautical tables are based on the Clarke 1880 spheroid that uses a
different compression to WGS84. The difference is small but noticeable (0.27 at 45°.)
For Clarke 1880 use e = 0.08248.
Celestial Navigation Calculations
An imaginary sphere surrounding the Earth is used for calculations. This is known as
the Celestial Sphere. Declination corresponds to Latitude and Hour Angles to
Longitudes.
Solving a sight uses spherical trigonometry. The triangle is known as the PZX
triangle.
The reason that the calculations do not allow for the shape of the Earth is because the
calculations are performed on the Celestial Sphere. Gravity ensures that the horizon
(and thus Altitude,) corresponds with the equivalents on the celestial Sphere.
Q1 to Q2 is the Equator and the Vertical line from P to G is the Greenwich Meridian
(0° GHA and Longitude.)
The compliment of an angle is 90° - the angle. This is also true of the sides in a
spherical triangle.
Compliments enable formulae to be simplified because the Sine of an angle equals the
Co-Sine of the compliment of that angle. This also applies to Tangents and
Cotangents, Secants and Cosecants.
SPHERICAL FORMULAE
For an Angle
Tan(C) = Sin(A)/ [Sin(b)/ Tan(c) – Cos(b) x Cos(A)]
Or
If angle A is 90° then a diagram is drawn with A above the circle and the sectors
filled with the adjoining sides.
Notice that the three sectors of the lower half are marked “co-” In other words the
compliment of these angles is used.
And
If co-Dec is the mid-part then the two opposites are Azimuth and co-Lat.
The Amplitude of a body is measured from East or West rather than North. In other
words Amplitude = 90° - Azimuth = co-Azimuth.
This gives the formula that many readers will be familiar with of
Sin(Amplitude) = Sin(Dec) x Sec(Lat)