3 Chapter
3 Chapter
3 Chapter
CHAPTER 3
GEODESY
Para
0301. Introduction and Scope of Chapter
0302. The Need for Geodesy
SECTION 1 - SPHEROIDS
SECTION 2 - DATUMS
0309. Datums
0310. Geodetic Datums
0311. Geodetic Datums Worldwide
0312. Datum Shifts
0313. Datum Transformation
3-1
A .. ... , ~n-t Q. \/o. rc i nn 1
BRd 45(1)
CHAPTER3
GEODESY
SECTION 1 - SPHEROIDS
GEOID
UNDULATION
'liEROfD (f:((fi
a. The Surfaces of the Earth. At least three distinct surfaces of the Earth can be
defined.
(1) The Physical Surface. The physical surface of the Earth, self-evidently,
is the land and water which physically makes up the Earth. This is not a perfet
sphere as it sits upon the oblate spheroid, and is then further deflected from the
mean sphere by the topography (mountains, bathymetry) and the fluid surfaces
of the oceans and ice caps. The physical surface of the Earth is thus both in
constant motion, and highly irregular in shape with local departures of several
kilometres from a pure sphere or even spheroid. The physical surface of the
Earth is thus not a suitable surface on which to base calculations, even with
present-day computing technology.
BRd 45(1)
(2) The Geoid. The geoid is a more regular shape, which is often considered
the 'true figure of the Earth'. It is the shape which would be defined by the
surface of the oceans under the influence of Earth's gravity and rotation alone, if
they extended in an unbroken line around the Earth (it assumes hypothetical
narrow canals through the continents rather than a 'water world'). Whilst it is
conventional to think of the geoid at sea level, it can be defined at any arbitrary
level of gravitational potential because what specifically defines the location of
the geoid is that it is an 'equipotential surface' meaning that all points on the
surface have the same effective potential (the sum of gravitational potential
energy and centrifu§al potential energy). Gravity is always perpendicular to the
geoid, which means that, due to the many density and gravitational anomalies
in the Earth, the geoid is still a highly irregular shape. The geoid is therefore still
not suitable to use as a computational surface. In fact, whilst it has long been
identified as a concept, it is only in recent decades with the advent of satellite
technology and sophisticated gravity sensing systems that the geoid has been
identified with any sort of confidence.
a Note. The term 'spheroid' is used in this book, as it clearly aligns with the
concept of the spheroid Earth rather than the spherical Earth. However, in
much of the literature relating to WECDIS and ECDIS systems it is
increasingly common to use the term 'reference ellipsoid' and it must be
empha§i_sed that the two terms have exactly the same meanina.
f = a-b
b
3-4
BRd 45(1)
<------ a
I= a -b
a
e =("':.Irr thus ez = az + bz
a2
A y
*
I a = Semi-major (Equatorial) Axis
I b = Semi-minor (Polar) Axis
b e = Eccentricity
I
X
s
,,,- j
c
e = ( a2 a~ b2 1/2
a2 -b2
e2
7
B
3-5
August 2018 Version 1
BRd 45(1)
b. Geocentric Latitude. The angle M0A which, by definition, must pass through
the centre of the spheroid 0 , describes the geocentric latitude of M and is
conventionally denoted by 9.
b2
tan 0 = ---:a-2--·tan <p
(1- f}2 tan.fJ
= (1- e) 2 tan0
d. Difference between Geodetic and Geocentric Latitude. There is no
difference between geodetic and geocentric latitudes at the Equator and the poles,
but elsewhere there is a non-zero difference between the two. This difference is at
its greatest when cp = 45°. When using WGS 84, the maximum difference is 11.54
minutes of arc.
3-6
August 2018 Version 1
BRd 45(1)
Geocentric Lat 9 K
A1I ~rY~ I
P1
b
tan/3' = a~tan({l
The difference between the geodetic latitude and parametric latitude is zero at the
Equator and at the poles, and has greatest value when <P= 45°. Using the WGS 84
spheroid, where f = 1/298.257223563, the greatest value is approximately 5.85
minutes of arc.
3-7
August 2018 Version 1
BRd 45(1)
W ( ±~-a 1H JE
Geodetic Lat 0,
Pl
a. Latitude Formula. The length of one minute of latitude may be found from the
formula po<P where p is the circular radius of curvature in the meridian and o(f) is a
small increase (measure in radians) in the geodetic latitude (f)_ It may be shown
that:
a(l - e) 2
p = (1 - e 2 sin2 q> ) 3 12
a(l - e)2
1; of latitude
(1 - e.r2 ~ini <p ) 312 sin 1'
When <P = zero:
3-8
August 2018 Version 1
BRd 45(1)
p
p'de
WI nl /~\ IE
Geodetic Lau,
P,
8
Note. Position C is the centre of a circle and arc pd<P is part of its
circumference. This on the spheroid (reference ellipsoid), the position of C
will move relative to O as arc E_d<P moves over the surface of the spheroid.
b. Longitude. 1' of longitude at the Equator= a sin 11• (minutes of arc). Thus, at
latitude <P:
a cos <p . ,
1' of longitude =-(1--e--;:;.2;-s-:-in---:;2;-<P:-:;);-;-1n.12 sin 1
0308. Geodesic
Just as a great circle gives the shortest distance between two points on the surface
of a sphere, a geodesic gives the shortest distance between two points on the spheroidal
Earth.
3-9
August 2018 Version 1
SECTION 2- DATUMS
0309. Datums
Charted data, both vertical and horizontal, must be referenced to a datum. Charted
depths and heights are referenced to a vertical datum, and this subject is covered in more
detail in Chapter 11 . Horizontal data is also reference to a datum; the horizontal or geodetic
datum refers to the spheroidal shape of the Earth used to establish the positional framework
for that chart.
-• "!!!l
l--'......,___.__.,..,,,,_..__.,,,-_,......, a,
.
S 11 9 \ \
----=..,, -,--o;-•.._...1
1°,
I.
When referenced to OSGB 36, position 'A' (S0° 21.S'N 004° 10.0'W) is 101 metres
from position 'B' (also so• 21.S'N 004° 10.0'W, when referenced to WGS 84 Datum)
3-10
August 2018 Version 1
BRd 45(1)
fi
consider the chosen spheroid as a rugby ball, and the local (or regional)
geodetic datum as a piece of orange peel placed over it where the position
and orientation fitted best. The precise observed position would then be a
pin stuck into the orange peel at the appropriate location. Elsewhere
however the oranqe oeel and ruqbv ball will be widelv seoarated.
d. The World Geodetic System (WGS). The World Geodetic System (WGS) is a
combined spheroid and datum. By design it is global in coverage, and is directly
related to the Earth's centre of mass; high accuracies in specific regions are
sacrificed in favour of a consistent accuracy across the entire planet.
(1) Concept. A true world geodetic datum/spheroid, fit for use with satellites
and satellite-derived positions, is defined by the system in which the satellite
orbital parameters are given. These in turn are dependent upon the precise co-
ordinates of the satellite tracking systems, the geopotential model of the Earth's
gravitational field, and a set of constants including GM (the Earth's mass M
times the gravitational constant G).
(2) WGS 84. From this satellite reference system data, a universal worldwide
spheroid and datum was established in 1984. The WGS 1984 (or almost
always just WGS 84) is the current system, and is NATO's preferred
spheroid/datum. It supersedes WGS 72, although some older model
navigational aids (such as legacy LORAN receivers in RN ships) may still
operate in WGS 72.
3-11
August 2018 Version 1
BRd 45(1)
e. Modern Regional Datums (Tied to GRS 80/WGS 84). In the modern era,
rather than rely upon networks of triangulation points ('trig' points, the familiar
concrete markers found on the tops of many hills in the UK), it is possible and
indeed common place to use satellite positioning systems with centimetric accuracy.
One of the consequences of the ubiquity and accuracy of this technology is that the
movement of the tectonic plates becomes both accurately detectable, and
geodetically significant. To account for this, regional geodetic datums tied to WGS
84 have been established at a specific time epoch, and then held fixed in relation to
the tectonic plate for the region . Whilst for many practical navigational purposes , the
modern regional datums may be considered to be equivalent to WGS 84, although
after three decades (or more) of divergence between the two then noticeable
differences are detectable even at mid-accuracy applications (precision farming, civil
engineering), so that the appropriate regional datum must be used (for example) for
terrestrial surveys. Provided the GNSS receiver and WECDIS are configured
correctly then this should be invisible to the operator, but may cause problems in
some circumstances.
a
Note. Strictly, many of the regional datums are based on the GRS (Geodetic
Reference System) 80 ellipsoid. WGS 84 was originally based on the GRS
80 ellipsoid when first defined, although since then the underlying parameters
of the two have drifted very slightly as WGS 84 has been refined. The
difference between the GRS 80 and WGS 84 spheroids amount to just over
0. 1mm on the semi-polar axis.
3-12
August 2018 Ver!':inn 1
BRd 45(1)
3-13
August 2018 Version 1
JRO 4!>\lJ
Datum Spheroid (a) Equatorial Flattening (f) & Eccentricity (e) Remarks
& (Type) Radius Reciprocal (1/f) & Squared (e2)
(b} Polar Radius f = (a-h)la e2 = (a2-b2)1a2
)SGB 36 Airy (a) 6,377,563.396m (f) 0.003340851 (e) 0.081673374 Some remaining
Regional) 1830 (b) 6,356,256.909m (1/f) 299.324964 (e 2) 0.006670540 UK charts.
)SSN 80 Airy (a) 6,377,563.396m (f) 0.003340851 (e) 0.081673374 Not for charts.
Regional) 1830 (b) 6,356,256.909m (1/f) 299.3249645 (e 2) 0.006670540 OSGB 36 related.
NAD27 Clarke (a) 6,378,206.4m (f) 0.003390075 (e) 0.082271854 Some remaining
(Regional) 1866 (b) 6,356,583.8m (1/f) 294.978698 (e 2 ) 0.006768658 USA charts.
Arc Clarke (a) 6,378,249m (f) 0.003407561 (e) 0.0824834 Some remaining
(Regional) 1880 (b) 6,356,515m (1/f) 293.465 (e 2 ) 0.006803511 S. African charts.
WGS72 WGS72 (a) 6,378,135m (f) 0.003352779 (e) 0.081818811 Being replaced by
(W'ldwide) (1972) (b) 6,356,750.520m (1/f) 298.26 (e 2 ) 0.006694318 WGS84.
NAD83 GRS80 (a) 6,378,137m (f) 0.00335281068 (e) 0.0818191910 Regional WGS 84
(Regional) (1980) (b) 6,356,752.3141m (1/f) 298.25722210 (e 2) 0.0066943800 for N. America
ETRS83 GRS80 (a) 6,378,137m (f) 0.00335281068 (e) 0.0818191910 Regional WGS 84
(Regional) (1980) (b) 6,356,752.3141m (1/f) 298.25722210 (e 2J 0.0066943800 for Europe
WGS84 WGS84 (a) 6,378,137m (f) 0.00335281066 (e) 0.0818191908 WGS 84 is NATO
(W'ldwide) (1984) (b) 6,356,752.3142m (1/f) 298.257223563 (e 2) 0.0066943800 preferred Datum.
0312.
ii Note. The defined parameters are 'a' and '1/f; all others are derived values.
Datum Shifts
a. Datum Shifts. Where two adjacent geodetic datums overlap, the same point
on the Earth's surface will have two distinct sets of geodetic co-ordinates, one in
relation to each geodetic datum. The difference between the two sets of co-
ordinates is called a 'datum shift'. The datum shift may not be consistent in
magnitude and direction (i.e. it is likely to very with geographic location).
., 4.
BRd 45(1)
0
0313.
Note. See also ADRS which provides further detail on this subject, including
the number of UKHO products currently based on each datum.
Datum Transformations
Several techniques exist to transform positional variety across a large area from one
geodetic datum to another (most commonly in the navigational context into WGS 84). These
are different from shifts, which apply to points or small areas, in that they operate across a
wide range of space and time. These are complicated operations which will not be
conducted manually, but can be carried out using appropriate software. It is rare that the
practical navigator will have to carry out such a transformation, but it is conceivable that the
NO might need to understand that such a transformation is required and request it from the
appropriate source.
b. Published Geographical Shifts. When the datum shift has been determined
by comparison of positions common to both datums, a simple block shift in
geographical coordinate values may be applied. As a result, many charts referred to
datums other than WGS 84 now show the datum shift between WGS 84 and the
horizontal datum in which the chart is published. These datum shifts are quoted to
an accuracy commensurate with the scale of that particular chart and should not
produce an error capable of being plotted on that chart at that scale.
3-15
August 2018 Version 1
BRd 45(1)
c. Datum Conversions. Full details may be found in NATO STANAG 2211 which
mandates that AGeoP-21 is the authoritative NATO document. STANAG 2211/
AGeoP-21 is available in the public domain via the NATO Standardisation Office
(NSO) website (https://nso.nato.inUnso/). To transform positions between different
geodetic datums, the geodetic datum coordinates of latitude and longitude are first
converted into three-dimensional Cartesian coordinates. Using the Molodensky
equations, translations along the X, Y and Z axes, rotations round these axes and a
scale factor are applied to the Cartesian coordinates. The modified Cartesian
coordinates are then converted into the second geodetic datum coordinates of
latitude and longitude.
3-16
August 2018 Version 1
BRd 45(1)
a. Principal Types of Earth Model. Naval systems concerned with what could be
termed 'navigational matters' (so WECDIS, command systems, fire control systems
etc)_ aU require an Earth model of some sort to be used in order to carry out
navigational calculations (whether of own ship movement, aircraft movement,
missile engagements etc). These fall into three broad categories; one could also
consider 'inertial space' as well although movements in inertial space would
inevitably have to then be converted into an Earth-referenced model in order to have
meaning for the navigator.
b. Choice of Earth Models. The choice of which Earth model to use is ultimately
driven by the purpose of the system, and hence the required accuracy, as well as
limits on the computational power available (although in practice this constraint
probably no longer exists today, but was a factor when legacy systems were
designed). As discussed at length above the spheroidal Earth model is the closest
approximation to the true Earth shape, and hence it is the most suitable surface on
which to base precise calculations for navigation, especially if they concern locations
at long ranges from the ship. Other Earth models are less accurate approximations
of the 'ground truth' and will give progressively less accurate results, especially
when dealing with long-range contacts beyond the horizon. The size of these errors
should be compared with those of the spheroidal Earth model.
3-17
August 2018 Version 1
BRd 45(1)
a. WGS 84. By default, the most commonly encountered spheroidal Earth model
in the Royal Navy will be WGS 84. This spheroid, based on a modified GRS 80
spheroid with an integrated gravity model, is a fully realised terrestrial reference
frame (TRF) and is the TRF on which GPS is based. In addition to being the basis
of GPS, WGS 84 is the NATO-preferred geodetic datum and spheroid for mapping
and charting products in NATO maritime areas, and it was adopted for all RN
applications in 1991. Furthermore, WGS 84 has been adopted by the IMO for
ECDIS, ECS, and AIS, and all electronic navigation charts (ENCs) and raster
navigation charts (RNCs) used with them.
Example. During a medium range gunnery exercise, a target position on the North East
corner of Garvie Island at 58° 37.09'N 004° 52.20'W (WGS 84) translates into a British
National Grid position of 233418E 973594N when using the correct Airy (1830) spheroid and
OSGB 36 datum. Fleet ARCS products are available with overprinted grids for certain
areas, but otherwise the transformation must be carried out mathematically. Some
transformation facilities are found within WECDIS (and may be available in other ECDIS
types), and are also available via the Fleet gunnery calculator. If the wrong spheroid/datum
combination is selected for the transformation then a significant error can occur; for example,
using the example of the target on Garvie Island, if the International Spheroid (1924) is used
instead of Airy (1830) the target's British National Grid position will be plotted 19 metres East
and 288 metres North of the correct position. Clearly, this sort of error is not acceptable
during gunnery, and so the utmost caution must be exercised when carrying out these sorts
of conversions (whether using WECDIS, the command system/gun system or the gunnery
calculator).
3-18
BRd 45(1)
Table 3-3. Errors Arising from a Spherical Earth Calculation compared with WGS 84
a
Note. In the above calculations and results, no account has been taken of
other sources of error. In particular, no account has been taken of the height
above (or below) the reference surface of either own ship or the contact.
This in itself will lead to further range errors, in addition to those directly due
to the use of approximate Earth mod~ls.
a. Flat Earth Assumptions. Flat Earth models assume a plane Earth, which then
enables the use of a two-dimensional system of Cartesian coordinates. These
Cartesian coordinates have axes arranged to lie along true (or theoretically
magnetic) North and due East from an arbitrarily chosen grid origin. This grid origin
will be chosen based on a number of factors, which may not always coincide with
own ship's position.
b. Errors Associated with Flat Earth Assumption. The errors introduced by the
flat Earth assumption are dependent upon the specific algorithms implemented in
the particular system. However, in general it will be deduced that the size of the
error will be a function of both the distance from the (arbitrarily selected) grid origin,
and the latitude of the grid origin. Thus errors in bearing and distance increase with
distance from the grid origin and with increase in the latitude of the grid origin. For
close ranges from own ship, the errors will either be negligible or insignificant for
most purposes.
3-19
Auam;t 2018 Version 1