Field

What the numbers refer to

A latitude and a longitude are numbers about a model, and the model has a shape, a placement, a set of marked points and a date. Change any of the four and the ground under the number moves.
The same coordinate on four datums. One pair of numbers — 2.0° west, 54.5° north — read as a coordinate on OSGB36, ED50, NAD27 and on WGS84, and plotted where each reading puts the mark on the ground. The spread runs to 195 metres. The numbers are identical; only what they refer to differs.

What a coordinate refers to

A latitude and a longitude are not a place. They are a claim about a model, and the model has a shape, a position, a set of marked stones and a date — change any one of the four and the ground under the numbers moves by hundreds of metres.

Datum · rung 1
OSGB36 to WGS84, one parameter at a time. Each bar is how far the mark at 2.0° west, 54.5° north moves under one of the seven parameters with the other six set to zero. seven of the seven are not zero for OSGB36. The units hide the comparison: one arcsecond of rotation moves this point 30.8 metres and one part per million of scale moves it 6.36 metres, so OSGB36's scale term contributes 130 metres — more than one of its three translations.

The seven parameters, and what each one does

A datum transformation is published as three metres, three arcseconds and one part per million, in three different units — which hides the fact that all seven are the same size of thing. Converted to ground displacement, the scale term beats two of the translations.

Datum · rung 2
Fitting a sphere to 10° of latitude at 54.5°. The Earth's Gaussian radius of curvature across the band, with two candidate spheres drawn against it. The global mean radius is 6371.0 kilometres and sits 2203 parts per million away from the ground here — 2.2 metres in every kilometre measured. The best local radius is 6385.0 kilometres and has no bias at all by construction, leaving 325 parts per million of spread that no sphere can remove, because the curvature varies across the band and a sphere's does not. The gain is a factor of 6.9.

A datum is fitted to a region

Forty countries adopted forty ellipsoids, and the usual explanation is that measurement was poor. It was not — a figure fitted to one country beats the global one over that country by a factor of seven, and what it removes is a systematic bias of two metres in every kilometre.

Datum · rung 2
What 4 parts per million of network strain leaves behind. Each arrow is where the best-fitting seven-parameter transformation leaves a marker, over 36 markers laid out as a grid across OSGB36's ground. The RMS residual is 1.64 metres and the worst is 3.00 metres. Seven parameters span the constant and linear parts of a displacement field; this one is quadratic, so no choice of the seven can reach it. The arrows are drawn 27× life size.

Where a fit leaves residuals

Seven parameters can carry a rigid motion and a size exactly. A triangulation network is neither, so the best possible transformation between two datums leaves metres on the table — in a pattern, not as noise — and which seven parameters come out depends on where the markers were.

Datum · rung 3
Molodensky against the exact route, NAD27. The distance on the ground between where each shortcut puts the transformed point and where the exact Cartesian route puts it, at 0 metres of ellipsoidal height, on a logarithmic scale. The full formulae stay within 0.8 centimetres across every latitude drawn; the abridged form, which replaces two ellipsoid-difference coefficients with one combined term, is worst at 40 centimetres — a factor of 48. The abridged curve dips near 45°, where the combined coefficient happens to equal the pair it replaces.

Molodensky's shortcut

The exact way to shift a datum needs an iteration nobody wanted to run on a 1950s machine, so a direct formula was derived instead. It lands within centimetres — and the abridged version everybody quotes is a different formula, worse by a factor of fifty.

Datum · rung 3
What grid spacing a tolerance buys. The worst bilinear interpolation error across OSGB36's ground, against the spacing of the table it is interpolated from, both axes logarithmic. The smooth shift falls by exactly four for every halving, which is what second-order interpolation does. Adding a 30 centimetre ripple 2 degrees across — three parts in a thousand of the 99 metre shift it rides on — moves the spacing needed for 20 millimetres from 1° to 0.125°, which is 64 times as many nodes.

When a formula is not enough

National mapping agencies distribute datum shifts as tables of numbers on a grid rather than as parameters. What decides the spacing that table needs is not the size of the shift — a hundred-metre shift tabulates coarsely, and a thirty-centimetre ripple riding on it does not.

Datum · rung 4
Where the ground is going, six plates. The velocity of the ground under a coordinate, computed as the cross product of each plate's rotation vector with the position, and drawn along the band where that rotation is fastest rather than over an outline this site does not hold. The quickest arrow is 75 millimetres a year on the Pacific plate, which is 1.89 metres in twenty-five years. drawn in Robinson — the arrows' directions are the projection's as much as the plates'.

The epoch is part of the coordinate

A latitude and a longitude with no date on them are incomplete, because the ground they refer to is moving at between ten and seventy millimetres a year. Australia crosses an ordinary survey tolerance in under two years and moved 1.8 metres between two versions of its own datum.

Datum · rung 4
Three surfaces, and the two heights between them. The ellipsoid, the geoid and the ground, with the two heights a coordinate can carry. Ellipsoidal height h is what a satellite fix returns and is measured from a surface defined by four constants. Orthometric height H is what a level and a staff measure and is referred to the geoid — the equipotential surface that best fits mean sea level. They differ by the separation N, drawn here as 45 metres because that is a stated input rather than a computed one: a geoid model is a data product with a truncation degree in it, and this site computes rather than downloads. The arithmetic h = H + N is exact whatever N is.

Height above what?

A satellite reports a height above a mathematical surface. A level and a staff report a height above the surface water settles on. The two disagree by tens of metres, both are correct, and only one of them decides which way a pipe drains.

Height · rung 1
Normal gravity, derived from four constants. Gravity on the surface of the level ellipsoid, by Somigliana's closed form, for WGS84. Nothing here is measured: a, f, GM and ω go in and the whole curve comes out, rising 5186 milligal — 0.53 per cent — from equator to pole. The two open marks are the published values of equatorial and polar gravity for WGS84, which the derivation reproduces to ten significant figures rather than borrowing.

The ellipsoid is a level surface

WGS84 publishes two dozen constants and defines four of them. The other twenty are consequences — polar gravity, the potential of the ellipsoid, the coefficient that dominates the Earth's gravity field — and every one comes back here from a, f, GM and ω to the last digit published.

Height · rung 2
Why a levelled height is not a distance. The correction between raw levelling and orthometric height, for lines at 200, 500, 1000, 2000 metres above the geoid running north from 50°. It is not instrument error: level surfaces converge towards the pole, so a run that stays on one of them gains height relative to another. A line 2000 metres up reaches 647 millimetres over 400 kilometres. The dashed line is 10 millimetres, which is about what a first-order levelling network closes to over that distance — so this is not a refinement, it is the larger of the two numbers.

A levelled height is not a distance

Level surfaces converge towards the poles by five metres in a thousand, so a chain of perfectly executed levelling observations does not sum to a height difference. The correction over four hundred kilometres of northing is larger than the network's own closure.

Height · rung 3
A deflection of 10 arcseconds, and what it hides. The ellipsoid normal and the plumb line at one point, with the geoid tilted against the ellipsoid by 10 arcseconds — drawn 3000× steeper than life, because at true scale the two lines are indistinguishable. The relation is exact and linear: an arcsecond of deflection is the geoid rising 4.85 millimetres in a kilometre, so 10 arcseconds is 48.5 millimetres per kilometre. A star sight measures the plumb line's direction, so astronomic latitude differs from geodetic by exactly this angle — 309 metres of ground, at a point where the coordinate itself is correct.

The plumb line is not the normal

A latitude measured from the stars and a latitude that means a position on the ellipsoid are different angles, because a plumb bob hangs along gravity and gravity is not perpendicular to a mathematical surface. Ten arcseconds of difference is 309 metres of ground.

Height · rung 3
The two reductions, at a grid factor of 1.0004. A measurement made on the ground has to be brought to the ellipsoid and then to the grid, and the two corrections have opposite signs. The elevation factor is R/(R+h) and always shrinks; the grid factor here is 1.0004, which stretches. Their product is the only number a surveyor can use. They cancel exactly at 2551 metres. At 1500 metres the combined factor is 165 parts per million, which is 1.65 metres on a 10 kilometre baseline.

The ground is not the grid

A tape measure on a hillside has to be brought down to the ellipsoid and then out onto the map, and the two corrections have opposite signs. On a grid whose scale factor exceeds one there is exactly one elevation where they cancel — 2,551 metres, for a factor of 1.0004.

Height · rung 3
The same coordinate on four datums. One pair of numbers — 2.0° west, 54.5° north — read as a coordinate on OSGB36, ED50, NAD27 and on WGS84, and plotted where each reading puts the mark on the ground. The spread runs to 195 metres, and the height changes by 51, 44, 47 metres, which is the component nobody quotes.

The third coordinate moves too

A datum shift is quoted as a horizontal displacement because horizontal is what people look at. The transformation acts on a three-dimensional point, and its vertical component is between a quarter and a half of the horizontal one — 51 metres, for a British coordinate.

Height · rung 4
Clairaut's theorem, and the term it drops. Clairaut's theorem says the flattening of a rotating body plus the flattening of the gravity on it equals five halves of the ratio of centrifugal to gravitational acceleration at the equator. Measured on WGS84: f = 3.3528e-3, f* = 5.3024e-3, and their sum is 8.65525e-3 against the theorem's 8.62447e-3. The residual is 3.08e-5, which is 2.74 times f² — second order, which is what a first-order theorem is entitled to be wrong by.

The flattening is not a free parameter

An ellipsoid is usually presented as two numbers somebody fitted. One of them is not free — Clairaut's theorem relates the shape of a rotating body to the gravity on it, and the relation holds on WGS84 with a residual of 3.1×10⁻⁵ — which is 2.74 times f², exactly what a first-order theorem is entitled to.

Ellipsoid · rung 5
How badly two arcs determine the flattening. The reciprocal flattening recovered by inverting two measured degree lengths, at 1.5° and 66.33°, against an error introduced into the equatorial one. The exact pair returns 298.26. Ten metres of error — 90 parts per million of a 110-kilometre arc — returns 301.5, and the relative error in the flattening is 118 times the relative error in the arc. At 1000 metres the inversion returns a negative flattening: an Earth longer through the poles than across the equator, which is the answer the Paris Observatory defended for a generation.

The figure of the Earth was measured

Two expeditions went to Lapland and Peru to measure the length of a degree of latitude, and the whole signal separating a flattened Earth from a round one is a kilometre in a hundred and eleven. Inverting two arcs amplifies their error by 118 — and a kilometre of error returns a lemon-shaped planet.

Ellipsoid · rung 5
Two parameter sets 100 metres apart, over the region they were fitted to. Each marker is drawn at a size proportional to how far the two transformations put it apart. The second set differs from the first by 100 metres of translation along the direction this network can least see, with the rotations and the scale re-fitted to absorb it — which is what a second agency's adjustment does when it chooses a different constraint. The worst disagreement anywhere in the region is 5.64 metres and the mean is 3.61. Applied at south-eastern Australia the same two sets differ by 193 metres, because the rotation that absorbed the translation here is a rotation of the whole Earth.

Two parameter sets, one transformation

Agencies publish seven-parameter datum transformations that differ by hundreds of metres in translation, and the usual reading is that one of them is better. Over the region either was fitted to they are the same transformation: a hundred metres of translation, re-absorbed by the rotations and the scale, moves a British coordinate by 5.6 metres and an Australian one by 193.

Datum · rung 5

A height that is not a length

Level surfaces converge polewards, so the surface a lake sits on is 5.28 metres lower at the pole than at the equator and a height system that reports lengths says a lake runs downhill. The fix reports a number that is constant on the surface and is not a distance from anything: a hundred-metre climb raises it by 99.73 metres at the equator and 100.26 at the pole.

Height · rung 4

A coordinate on another body

Mars publishes two latitudes for every point and they differ by up to 0.338°, which is twenty kilometres of ground — almost exactly the same distance as the Earth's own 0.192°, because Mars is smaller by nearly the same factor. The libraries that compute either have known both numbers since this collection's early essays and had never been pointed anywhere but here.

Bodies · rung 1

The same projection on a different body

A projection's distortion does not depend on the body's size at all — only on its flattening — so Web Mercator's 0.3848° of angular deformation is not a fact about Web Mercator. It is 2f radians, and the same mistake made on Mars measures 0.6765°, on Jupiter 7.68°, and on a body scaled ten times larger than Mars the identical number to twelve figures.

Bodies · rung 2

A body that is not an ellipsoid

Vesta's three axes are 286.3, 278.6 and 223.2 kilometres, all different, so it has no axis of revolution — and on such a body the planetographic latitude of a point at 45° planetocentric swings by 1.4° as one walks round it in longitude, and by 6.4° on Phobos. Latitude stops being a function of position.

Bodies · rung 3

A map of a body with three axes

Every projection on this site is written in a latitude, and a body with three unequal axes has none: the coordinate lines are not perpendicular, Lambert's equal-area formula spreads areas by 28 per cent on Vesta, and the equal-area map that does work has a different one on every meridian.

Bodies · rung 4

A body with no sea level

On Earth the zero of height is found rather than chosen — water settles onto the equipotential surface by itself. Nowhere else has one, and the difference is measurable: Clairaut's theorem predicts the Earth's flattening from its own gravity field to twelve metres at the pole and misses Mars's by 2.2 kilometres.

Bodies · rung 5

A deflection is the slope of a mass

The site has computed the exact relation between a geoid slope and a deflection of the vertical — one arcsecond is 4.85 millimetres in a kilometre — and computed no deflection anywhere, because a deflection is the gradient of a geoid and the geoid is a citation. State a mass instead, and every quantity is a closed form: a salt dome five kilometres across bends the plumb line by 2.21 arcseconds.

Height · rung 5

A chain of transformations does not close

Two datum transformations applied one after the other are not the sum of their fourteen published numbers. The gap is seven millimetres in Britain, doing the two in the other order moves the answer twenty, and the rotation matrix everybody prints is not a rotation.

Datum · rung 6

A conformal map of a body with three axes

This site built an equal-area map of a triaxial body and wrote down what it could not do: the conformal one, which needs an isothermal coordinate that a surface with no axis of revolution was said not to have. It has one, Jacobi found it in 1839 for a different reason, and it takes two one-dimensional integrals.

Bodies · rung 6

The line a height is measured along

Five rungs have argued about the surface a height is measured *from* and every one of them took the line it is measured *along* to be straight and known. It is neither: through a stated buried mass a plumb line arrives 47 millimetres from the point below the summit, and the height it gives depends on the density of rock nobody has seen — 342 millimetres of spread at 4,000 metres and 1.37 metres at 8,000.

Height · rung 6

A conformal map of a body that is not a quadric

Jacobi's ellipsoidal coordinates give a triaxial body a conformal map by two quadratures, and this collection wrote down what that argument uses: the surface has to be a quadric. A real body is not. Solving the discrete Cauchy–Riemann equations instead — one complex equation per triangle, two thousand triangles, conjugate gradients — gives a conformal map of a bumped body to a median of 1.10°, converging at first order in the mesh, with the areal factor spreading by 3.07 and refusing to converge at all.

Bodies · rung 7

The seven parameters have their own uncertainty

Nine essays on this ladder print a datum transformation as seven exact numbers. Every published set is the output of a least-squares fit and arrives with standard errors as much a part of the result as the parameters — and pushing those widths through to the ground gives an ellipse, not a number, that is 68 mm across at the equator and 43 mm at 70°.

Datum · rung 7

The map depends on where it was cut

The previous rung solved the discrete conformal equations on a triangulated body and reported an areal spread of 3.07, then recorded that the number might belong to the patch, the boundary and the two pinned vertices rather than to the surface. It belongs to the pins: hold a different pair and it runs from 2.68 to 7.44, while the typical angular deformation does not move at all.

Bodies · rung 8

How far the plumb line bends

The previous rung dropped a plumb line down a four-kilometre column of rock beside one buried mass and found it arrived 47 millimetres from the point below the summit. That is a number about that mass. Parameterising by the deflection of the vertical instead — the quantity surveyors actually measure — gives 12.54 mm per arcsecond, exactly linear across a fortyfold range.

Height · rung 7

The geoid model stops at a degree

Eleven essays treat the geoid as a surface that exists. Every geoid anybody uses is a series truncated at a degree, so every orthometric height derived from one carries an omission error nobody quotes with the height — eighteen centimetres at degree 360 — and the same truncation removes two thirds of the slope, which does not converge at all.

Height · rung 9

Cuts of the same size in different places

Where a body is cut decides how well it can be mapped, by a factor of ten — established with five windows of five different sizes, so *where* and *how much* were confounded and the factor could have been entirely about extent. Held to the same surface area to a quarter of a per cent, the answer survives at a factor of 1.32, and what predicts it is the curvature the window encloses.

Bodies · rung 9

A mountain is not a buried sphere

The plumb line's drift was measured over a compact buried body and grows as the 0.69 power of the column's height — an exponent that is a statement about how quickly a buried sphere's field weakens with distance rather than about mountains. Spread the same mass into a crustal root and the exponent climbs to 0.84, while the proportionality to the deflection survives exactly.

Height · rung 8

A long window and a square one

Rung nine held the patch's area and found that where the cut goes still changes the map by a third, with the curvature it encloses predicting the change at r = 0.969. It recorded that it had held area and curvature and not shape. Sweeping the shape at constant area separates two things that had looked like one: the map is curvature and the cost is shape.

Bodies · rung 10

What another common point buys

Rung three finds that a seven-parameter datum fit leaves a pattern rather than noise. Six per cent of the residual it reports is the transformation's own error and the other ninety-four is distortion no seven parameters can follow — so adding common points improves a term that was already small and cannot touch the one that is quoted.

Datum · rung 8

A longitude that drifts with the rotation rate

Ten essays here map bodies whose shape is the problem. A longitude is not about shape: it is a landmark plus an extrapolation over however many days have passed, and an error in the last published decimal of a rotation rate is a coordinate error that grows without bound in time.

Bodies · rung 11

The rotation has two sign conventions

Eleven essays price a datum transformation's parameters, their fit, their residuals and their uncertainty. None asks what the numbers mean: three of the seven are published under two conventions whose rotations differ in sign, and reading one set with the other formula costs exactly twice the rotation — 55 metres on OSGB36 and 152 on DHDN.

Datum · rung 9

The radius of curvature is two numbers

Twelve essays work on the ellipsoid and every one of them takes a radius when it needs one. There are two at every point, they differ by 42.70 kilometres at the equator, and the single number every table prints is out by 5,583 parts per million on a line running north.

Ellipsoid · rung 9

The correlation is not the same in every direction

Every number in the previous rung came out of Σ cₙ Pₙ(cos ψ) — a covariance that depends on the angular distance and nothing else. Ground has grain: at a modest anisotropy the same 20 km baseline omits 165 millimetres along it and 241 across, and the isotropic answer understates the worse direction by 18.4 per cent.

Height · rung 10

The ground has an indicatrix too

Two hundred and forty-six essays hold the Earth still while the page is measured. The ground is moving at tens of millimetres a year, and the motion is a map with four partial derivatives — so the same construction that draws Tissot's ellipse draws one for the ground, and the fastest plate in the model returns exactly nothing.

Strain · rung 1

A velocity needs a frame and a strain rate does not

One place on a plate boundary moves at 5.3, 7.0, 10.2 or 54.6 millimetres a year, and towards the south-east or the north-west, depending on which rotation was subtracted first. Every one of those readings reports the same strain rate, to two parts in a hundred thousand million.

Strain · rung 2

A grid stops fitting the ground it was laid on

A hundred-kilometre baseline across a plate boundary changes by 815 millimetres in fifty years, and the same baseline on the fastest plate in the model changes by nothing at all in the same fifty. The interval before an ordinary first-order specification is broken is six years in one place and never in the other.

Strain · rung 3

The strain a map adds to the ground's

A ground carried rigidly at forty millimetres a year deforms nothing, and read off a Mercator sheet at sixty degrees north it reports 10.9 nanostrain a year. Along a profile through a real boundary the invented part is 0.0 per cent of the answer where the zone is loud and 884 per cent where it is quiet.

Strain · rung 4

A vertical rate needs a height system

Four rungs of this anchor measure the two-by-two horizontal tensor, because that is what a tangent chart returns. The larger signal in a subsiding basin is vertical — 126 nanoradians a year of tilt against 0.76 nanostrain a year of horizontal strain — and it is not a measurement at all until the surface it is measured against is named, because that surface is moving too.

Strain · rung 5

Every country's zero is a different surface

Ten rungs measure a height against a geoid and treat the geoid as one object. No national datum is on it: each is pinned to the mean sea level at one tide gauge, the sea surface stands up to two metres from the geoid, and the eight European zeros measured here spread over 462 millimetres — a step no levelling can remove and nobody's error.

Height · rung 11

The parameters are not independent

Rung seven gives the seven parameters their own uncertainty and stops at seven numbers. There are twenty-eight, and the twenty-one nobody publishes are not small: a translation and the rotation that mimics it correlate at 0.94, the normal matrix has a condition number of 4 × 10¹⁶, and propagating from the diagonal alone overstates the transformation's uncertainty by up to a factor of thirty-six.

Datum · rung 10

A ray from the centre hits the surface twice

Eleven rungs map bodies that are lumpy, triaxial and turning at a drifting rate, and every one assumes the surface is star-shaped about the centre — which is what makes a longitude and a latitude a coordinate at all. A contact binary is not: on a stated body with a neck a third of a lobe wide, 10.9 per cent of the sky has no single radius, and the shape model everybody publishes fills the neck in and adds 1.67 per cent of the volume.

Bodies · rung 12

On a body with a hole, north can be up everywhere

Twelve rungs vary the body's shape and none varies its topology, which is what every impossibility here actually rests on. A torus has a nowhere-zero tangent field, a total curvature of zero rather than 4π, and a conformal world map with no cut and no singular point — and it still cannot be flattened, for the one reason that survives.

Bodies · rung 13

The height a coordinate does not carry

Ten rungs treat a datum shift as a map from one pair of angles to another. It is not one: the transformation runs through Cartesian coordinates, so the answer depends on the height — by 7.7 to 121 millimetres per kilometre depending on the datum and the place, which puts Lhasa 432 millimetres from where a height of zero would have said.

Datum · rung 11

49 essays in this field, the first 16 of them shown with their opening figure.

All essays