Learn
Foundational articles on coordinates, projections, datums, GPS, and time zones. Every article cites primary sources; per-article accuracy expectations are documented in /accuracy.
Foundations
What Is Latitude and Longitude?
Latitude and longitude defined — angular coordinates on the WGS-84 ellipsoid, how they're measured, what datum modern GPS uses, and the misconceptions to avoid.
What Is Latitude?
Latitude explained: the angle north or south of the equator, geodetic vs geocentric, why one arcminute equals one nautical mile, and degree-by-degree distances.
What Is Longitude?
Longitude explained: the angle east or west of the prime meridian, the 102 m IERS shift from the Airy line, why meridians converge, and the link to time zones.
Coordinate Formats Explained
Six common coordinate formats compared — decimal degrees, DMS, DDM, UTM, MGRS, Plus Codes — with the Empire State Building in each and ISO 6709 ordering rules.
How to Read Coordinates
A practical guide to decoding a coordinate when you see one — identify the format (DD, DMS, DDM, UTM, MGRS, Plus Code), decode each component, recognise the precision, and avoid the common reading mistakes.
How to Write Coordinates
A practical guide to producing well-formed coordinates yourself — choosing the format, picking the right precision, formatting components correctly, and avoiding the typographic and convention errors that break downstream parsers.
Decimal Degrees vs DMS (and DDM)
The conversion math between the three angular coordinate formats — decimal degrees (DD), degrees-minutes-seconds (DMS), and degrees-decimal-minutes (DDM). Formulas with worked examples, precision trade-offs, and the conversion mistakes that produce off-by-one bugs.
Why Latitude Comes First (and When It Doesn’t)
ISO 6709 specifies latitude before longitude; GeoJSON, WKT, PostGIS, and KML specify longitude first. Where each convention came from, the 90° detection rule, and the defensive patterns that prevent order-swap bugs.
Positive and Negative Coordinates: Signs and Hemisphere Letters
The two equivalent ways to encode hemisphere in coordinates: signed numbers (positive north / east, negative south / west) and hemisphere letters (N, S, E, W). Conversion rules, when each is used, the sign-loss bugs, and the antimeridian discontinuity.
Precision vs. Accuracy in Coordinates
Precision is how finely a coordinate is recorded; accuracy is how close it is to the true position. Per NIST, the two are independent — and consumer-GPS coordinates routinely conflate them. The 5-decimal rule of thumb, three worked examples, and why dataset provenance matters more than decimal count.
Coordinate Systems
Coordinate Systems Overview: How Locations Become Numbers
The four families of coordinate system — geographic, projected, grid-based, code-based — per ISO 19111, with EPSG codes and the WGS-84 vs NAD83 datum split explained.
Geographic Coordinate Systems: Latitude and Longitude on an Ellipsoid
A geographic coordinate system locates points using angular latitude and longitude referenced to an ellipsoid plus a datum. WGS84, NAD83, and ETRS89 are the dominant examples. The article covers the ellipsoid model, axis-order conventions, units, ellipsoidal height, and when geographic is the right choice (almost always for exchange) vs. projected.
Projected Coordinate Systems: Earth Flattened onto a Plane
A projected coordinate system maps angular latitude and longitude onto a flat plane through a map projection, then records positions in linear units (metres or feet). UTM, Web Mercator, State Plane, the British National Grid, and Lambert-93 are the dominant examples. Every projection trades off shape, area, distance, and direction — and no projection preserves all four.
The UTM Coordinate System: 60 Zones, Transverse Mercator, Sub-Metre Accuracy
The UTM grid system explained — 60 zones, 6° wide, transverse Mercator math, the 0.9996 scale factor, Norway/Svalbard exceptions, and the UPS polar handover.
MGRS Explained: The NATO Military Grid Reference System
The Military Grid Reference System (MGRS) is the NATO grid built on UTM and UPS. It encodes a position as a 5-to-15-character alphanumeric string — Grid Zone Designator, 100 km square letters, and a variable-precision numerical pair — designed to be unambiguous over voice radio and printed media. The article covers the three-layer grid, the 100 km square lettering scheme, precision levels from 100 km down to 1 m, and the polar-cap UPS handover.
WGS 84 Explained: The Datum Every GPS Uses
WGS 84 explained — the GPS-broadcast datum on a 6,378,137 m semi-major-axis ellipsoid, the seven realizations from 1984 to G2139, and the relationship to ITRF.
NAD83 Explained
NAD 83 explained — US federal horizontal datum on GRS80, plate-fixed to North America, ~1-2 m offset from WGS-84 in CONUS, and replaced by NATRF2022 in 2025-2027.
WGS 84 vs NAD83
WGS 84 and NAD83 use virtually identical ellipsoids and were aligned within a metre in 1986. They have since diverged by 1-2 metres in the conterminous US because NAD83 is anchored to the North American Plate (which drifts ~2.5 cm/year relative to the global reference) while WGS 84 tracks the International Terrestrial Reference Frame. This support covers the divergence, when it matters, how to convert between them, and which datum is appropriate for which use case.
ETRS89 Explained
The European Terrestrial Reference System of 1989 is the European counterpart to NAD83 — a geodetic reference frame anchored to the stable part of the Eurasian Plate. Aligned with ITRF89 at epoch 1989.00, ETRS89 has since diverged from the global ITRF by about 75 cm due to Eurasian plate motion. It is the standard reference for the EU INSPIRE directive and for nearly all European national geodetic systems, and underpins the European-wide LAEA, LCC, and Transverse Mercator projections used for cross-border geospatial work.
The State Plane Coordinate System
The State Plane Coordinate System is a set of about 124 projected coordinate zones covering the US — one or more per state — designed so that scale distortion stays below 1 part in 10,000 within each zone. Each zone uses Lambert Conformal Conic, Transverse Mercator, or (for Alaska) Oblique Mercator depending on the state's shape. SPCS is the working coordinate system for nearly all US land surveys, civil engineering design, and state-level GIS. The 2022 SPCS modernisation is now rolling out alongside the broader NSRS Modernization.
The British National Grid
The British National Grid is the UK's national projected coordinate system, used by Ordnance Survey, HM Land Registry, and UK government data. It uses Transverse Mercator on the OSGB36 datum and the Airy 1830 ellipsoid — an older, regionally-fitted ellipsoid that diverges from WGS 84 by up to 120 m at the UK's corners. The grid uses two-letter 100 km squares (SP, SU, TQ, etc.) with numeric east-north references. Ordnance Survey provides the OSTN15 transformation grid for converting between OSGB36 and ETRS89/WGS 84.
Quadkey and Tile Coordinates
Tile coordinates are the discrete coordinate system used by every modern web map. Each tile is identified by a (z, x, y) triple — zoom level, column, row — based on a Web Mercator projection of Earth. Microsoft's Quadkey is an equivalent base-4 string encoding that captures the hierarchical pyramid: each (z, x, y) tile maps to a unique string of length z. This support covers the math, the pyramid structure, the Quadkey encoding, and why tile coordinates are a coordinate system in their own right.
Projections
What Is a Map Projection?
Map projections explained — Gauss's Theorema Egregium, the four classes (conformal, equal-area, equidistant, compromise), and which projection fits which task.
The Mercator Projection
Mercator projection explained — Gerardus Mercator's 1569 conformal cylindrical chart, the sec(φ) scale factor that inflates Greenland, and Web Mercator differences.
Web Mercator (EPSG:3857): The Projection of Every Web Map
Web Mercator is the spherical Mercator variant that Google Maps introduced in 2005, now used by every major slippy-map provider — Google, Mapbox, Apple, Bing, OpenStreetMap. EPSG:3857. The article covers the simplified spherical math applied to the WGS 84 ellipsoid's radius, the ±85.0511° cutoff that makes the projection square, the zoom-level tile arithmetic that drove the design, the small high-latitude differences from classical Mercator, and when Web Mercator is the right (or wrong) choice.
The Peters Projection
The Gall-Peters projection is an equal-area cylindrical projection — every region has the correct proportional area, at the cost of severe shape distortion. James Gall described it in 1855; Arno Peters reintroduced and popularised it in 1973, accompanied by a political argument that the Mercator projection was Eurocentric. The 1980s 'map wars' between Peters and the cartographic establishment ended with the recognition that Peters is one of many equal-area projections and that no projection can preserve every property simultaneously.
The Robinson Projection
The Robinson projection is a pseudocylindrical compromise projection — it preserves no single property but distorts shape, area, distance, and angle each by a moderate amount, producing a visually pleasing world map. Arthur H. Robinson designed it in 1963 by empirical visual trial rather than closed-form mathematics; the projection is defined by a lookup table. National Geographic adopted it as the standard world-map projection from 1988 to 1998, when it was replaced by Winkel Tripel.
The Winkel Tripel Projection
The Winkel Tripel projection is a compromise pseudoazimuthal world-map projection designed by Oswald Winkel in 1921 as the arithmetic mean of the equirectangular and Aitoff projections. 'Tripel' is German for the threefold minimisation of area, direction, and distance distortion. National Geographic adopted it in 1998, replacing Robinson. It distorts each property by a moderate amount, achieves better balance than Robinson in numerical distortion measures, and is defined by a closed-form formula rather than a lookup table.
Equal-Area Projections
Equal-area projections preserve area exactly: every region on the map has the same proportional area as on Earth. The trade-off is shape distortion. This support surveys the major families (cylindrical, pseudocylindrical, azimuthal, conic, interrupted, and the modern Equal Earth) and gives practical-choice guidance for thematic mapping. Any map showing area-derived quantities — population density, GDP per area, rainfall — must use an equal-area projection or systematically misrepresent the data.
Why Greenland Looks Huge on Maps
On a typical world map — the Mercator projection used in classrooms and on Google Maps — Greenland looks about the same size as Africa. In reality Africa is 14 times the size of Greenland. The reason is mathematical: Mercator's area scale at latitude φ is sec²(φ), which means areas at 70° N are inflated by a factor of nearly 10. This article explains the math accessibly, surveys the other countries Mercator visually distorts (Russia, Canada, Antarctica, Alaska), and points toward equal-area alternatives.
Conformal Projections
Conformal projections preserve angles locally — the angle between any two curves on Earth equals the angle between their images on the map. The trade-off is area distortion, sometimes severe. This support surveys the major conformal projections (Mercator, transverse Mercator, oblique Mercator, Lambert conformal conic, stereographic), explains why navigation and aeronautical charts depend on the property, and contrasts the family with equal-area projections.
Cylindrical vs Conic vs Azimuthal Projections
Map projections are classified by the developable surface they are conceptually built on: a cylinder (cylindrical projections), a cone (conic), or a plane (azimuthal). Each shape has a natural latitude band where it works best. Pseudocylindrical, pseudoconic, and pseudoazimuthal variants relax the strict definitions. Combined with the conformal/equal-area/compromise property taxonomy, the shape classification gives a complete 2D framework for understanding any projection.
Datums
What Is a Geodetic Datum?
Geodetic datums explained — ellipsoid + origin + orientation + gravity model, horizontal vs vertical, WGS-84 vs NAD83, and how datums are realized via GNSS.
The Geoid Explained
The geoid is the equipotential surface of Earth's gravity field that best approximates mean sea level — the irregular reference surface that defines orthometric heights. It rises by up to ~85 m near Iceland and dips to ~−106 m south of India, owing to uneven mass distribution inside Earth. Modern global geoid models (EGM2008, GEOID2022) are computed from satellite gravity missions (GRACE, GOCE) plus terrestrial gravity surveys.
Ellipsoid vs Geoid: Two Reference Surfaces, Two Jobs
The ellipsoid is a smooth mathematical surface defined by two parameters; the geoid is the irregular gravity-equipotential surface that approximates mean sea level. The article compares them side-by-side, derives the fundamental relationship h = H + N between ellipsoidal and orthometric heights, and explains why coordinate work needs both surfaces.
Why the Earth Is Not a Sphere
Earth is an oblate spheroid: slightly flattened at the poles and bulged at the equator. The equatorial radius (6,378.137 km) is about 21.4 km larger than the polar radius (6,356.752 km), a flattening of 1/298.257223563 that's the consequence of Earth's rotation. The article covers the historical Cassini-vs-Newton dispute that settled the geometry, the modern reference ellipsoids that quantify it, comparison with other rotating bodies, and why every coordinate reference system uses an ellipsoid rather than a sphere.
The Reference Ellipsoid
A reference ellipsoid is a mathematical primitive: a smooth surface defined by two parameters — the semi-major axis a and the flattening f. Every modern geodetic datum builds on one. The article covers the defining and derived parameters, the historical list of ellipsoids still in use (Clarke 1866 through GRS80 and WGS 84), how an ellipsoid relates to a datum, and how to pick one for a project.
Datum Transformations: Converting Between Reference Frames
A datum transformation converts coordinates from one datum to another — WGS 84 to NAD 83, ETRS89 to ITRF, NAD 27 to NAD 83. The article covers the three main methods (three-parameter translation, seven-parameter Helmert / Bursa-Wolf, and grid-based like NTv2 / GEOID22), when to use each, the NGS NCAT / NADCON / VERTCON US-authoritative tools, and how PROJ implements the EPSG-registered transformations globally.
Horizontal vs Vertical Datum
Coordinate work tracks horizontal and vertical datums independently: a horizontal datum (WGS 84, NAD 83, ETRS89) for latitude / longitude and a vertical datum (NAVD 88, EGM2008, Newlyn) for height. The article covers why the two are separated, how they pair in practice (NAD 83 + NAVD 88 is the US default), the differences between geometric and gravimetric vertical datums, and the conversion between ellipsoidal and orthometric heights.
NSRS Modernization: Replacing NAD 83 with NATRF2022
The US National Geodetic Survey is replacing NAD 83 (horizontal) and NAVD 88 (vertical) with new reference frames in 2024–2026: NATRF2022 for the North American plate, plus PATRF2022 / CATRF2022 / MATRF2022 for Pacific / Caribbean / Mariana territories, and NAPGD2022 as the new geopotential / vertical reference. The article covers why, what's changing, the timeline, and the implications for surveyors, infrastructure records, and software.
Distance
Great-Circle Distance
Great-circle distance — haversine on a sphere vs Vincenty/Karney on the WGS-84 ellipsoid, the 0.5% accuracy gap, and the JFK-LHR worked example (~5,541 km).
The Haversine Formula Explained
The haversine formula computes great-circle distance on a sphere using a numerically stable trigonometric identity. The article derives the formula from spherical geometry, explains why it's preferred over the spherical law of cosines, gives a worked example, traces the history from Inman's 1835 navigation tables to modern JavaScript implementations, and covers the practical considerations (mean Earth radius, edge cases, performance).
The Vincenty Formula Explained
Vincenty's 1975 formula computes geodesic distance on a biaxial ellipsoid like WGS 84, iterating to millimetre accuracy. The article covers the inverse and direct solutions, the convergence behaviour, the well-known antipodal failure mode (and how robust implementations fall back to haversine), Karney's 2013 algorithm as the modern successor, and a worked Empire State → Statue of Liberty example.
Why Flight Paths Look Curved on Flat Maps
The 'curved' flight paths on airline route maps aren't curved — they're the straightest path possible, distorted by the flat-map projection (usually Mercator). Great-circle paths look bowed toward the nearer pole on Mercator because the projection stretches high latitudes. The article explains the geometry with worked examples (LA → Tokyo, NY → Hong Kong) and a few side notes on jet streams and ATC constraints.
What Is a Rhumb Line?
A rhumb line (loxodrome) is a path on Earth that crosses every meridian at the same angle — a constant compass bearing. On a globe it spirals toward the pole; on a Mercator map it appears as a straight line, which is the property Gerardus Mercator engineered his 1569 projection to produce. The article covers the definition, the geometry, the rhumb-vs-great-circle distance comparison, the closed-form distance formula, and when rhumb lines are still preferred (marine navigation, short legs).
Initial and Final Bearing on a Great Circle
Along a great-circle path on a sphere, the compass bearing changes continuously: the initial bearing at the start point and the final bearing at the destination differ. The article covers the forward-azimuth formula, why the two bearings differ, the reverse-azimuth relationship, the magnetic-vs-true distinction, and a worked JFK → LHR example showing initial 51° and final 109° bearings.
Great-Circle vs Rhumb Line: Side-by-Side Comparison
Two ways to traverse Earth's surface: the great-circle path (shortest) and the rhumb line (constant bearing). The article compares them on every axis — distance, bearing behaviour, appearance on different projections, historical adoption, modern use. Worked distance tables show the rhumb-line penalty ranges from negligible at short distances to ~25–30 % on transoceanic paths between distant latitudes.
The Antipode of a Location
The antipode of a point is the geometrically opposite point on Earth — the place you would arrive if you tunneled straight through the centre. The formula is simple: negate the latitude and shift the longitude by 180°. The article covers the geometry, famous antipodal pairs (NYC ↔ off-Western-Australia Indian Ocean; London ↔ Antipodes Islands east of New Zealand; Madrid ↔ Wellington area; Shanghai ↔ near Buenos Aires), why most antipodes are ocean (Pacific bias), and why antipodes matter for distance algorithms.
GPS
How GPS Works
How GPS works — 31 satellites at 20,200 km altitude, atomic-clock trilateration, the four-unknown equation, relativistic corrections, and the ~5 m accuracy budget.
GPS Accuracy Explained
GPS accuracy by the numbers — ~4.9 m smartphone (95%), ~1-2 m with SBAS, ~1-2 cm with RTK, the full error budget (ionosphere, multipath, geometry), and what's tunable.
Why GPS Is Not Always Accurate
GPS sometimes places you in the wrong building, on the wrong side of a road, or several blocks off. The article covers the practical failure modes — urban-canyon multipath, dense canopy attenuation, indoor signal loss, ionospheric disturbances, atmospheric effects, jamming and spoofing — with environment-specific accuracy degradation and the user-side mitigations.
Assisted GPS (A-GPS) Explained
Assisted GPS (A-GPS) is what makes your smartphone fix its location in 1–3 seconds rather than 30–60. Instead of waiting to download satellite ephemeris from the satellites themselves, the phone fetches it from the network (over cellular or Wi-Fi). The article covers how A-GPS works, the 3GPP specifications that standardise it, the difference between control-plane and user-plane A-GPS, and what happens when the network is unavailable.
GPS vs GNSS: The Five Global Satellite Constellations
'GPS' is the colloquial term for satellite positioning, but it specifically refers to the US Department of Defense system. GNSS (Global Navigation Satellite System) is the generic term for all such systems — GPS, GLONASS, Galileo, BeiDou, plus regional QZSS and IRNSS. Modern multi-constellation receivers track all of them simultaneously for better accuracy and reliability.
Differential GPS (DGPS) Explained
Differential GPS uses a fixed reference station with known coordinates to compute real-time GPS error corrections that nearby roving receivers apply for improved accuracy. Originally deployed by the US Coast Guard in the 1990s for maritime navigation; now superseded by SBAS / WAAS for most users but still the foundational technique behind every modern GPS augmentation system.
RTK GPS: Centimetre Accuracy in Real Time
Real-Time Kinematic (RTK) GPS achieves 1–2 cm horizontal accuracy by using carrier-phase observations and a real-time correction stream from a fixed base station. The article covers how RTK works (carrier-phase ambiguity resolution), the base-rover geometry, RTK vs network RTK (VRS) vs PPP-RTK, the NTRIP distribution standard, the typical 30 km baseline limit, and the surveying / agriculture / autonomous-vehicle use cases that drive RTK demand.
GPS Jamming and Spoofing
GPS jamming (drowning out the signal with noise) and spoofing (transmitting fake satellite signals) are increasing threats to civilian navigation. The article covers how each works, the documented incidents (Black Sea maritime 2017+, Middle East aviation 2023+, US truck-driver jammers), why civilian GPS has no authentication, and the mitigations: M-code military signals, RAIM integrity monitoring, multi-GNSS, inertial backup, anti-jam antenna arrays.
GPS Multipath Error
Multipath is the dominant urban-environment GPS error: signals reflecting off buildings, water, glass, or other surfaces arrive at the receiver delayed, causing inflated pseudorange measurements. The article covers the geometric mechanism, the typical magnitudes (1–10 m in cities, sub-metre open sky), receiver-side mitigations (CRPA antennas, multipath estimator, carrier-phase smoothing), site-selection rules for survey work, and why multipath remains the hardest civilian GPS error to eliminate.
Dilution of Precision (DOP) in GPS
Dilution of Precision (DOP) quantifies how much the geometric arrangement of GPS satellites in the sky multiplies the user-equivalent range error (UERE) into final position error. HDOP, VDOP, PDOP, TDOP, GDOP — each measures a different component. Good geometry: DOP < 2; poor: DOP > 6. The article explains the math, the interpretation, modern multi-GNSS DOP improvements, and the practical rules of thumb for survey-grade work.
GLONASS Explained
GLONASS — the Russian Global'naya Navigatsionnaya Sputnikovaya Sistema — is the second-oldest operational GNSS, launched by the Soviet Union and revived by Russia in the 2000s. 24+ satellites at ~19,100 km in three orbital planes inclined 64.8° (better polar coverage than GPS's 55°). Originally FDMA across frequency channels per satellite; modernized GLONASS-K2 adds CDMA. Similar accuracy to GPS (~5 m unaugmented). Routine multi-GNSS component in every modern receiver.
The Galileo Satellite System
Galileo is the European Union's civilian-controlled GNSS. 30 satellites at ~23,222 km altitude in three orbital planes inclined 56°, broadcasting on E1 / E5a/E5b / E6 frequencies. Initial Operational Capability declared 2016; Full Operational Capability progressing through the 2020s. Five service tiers: Open Service (free, 1–4 m), Public Regulated Service (encrypted, government), High Accuracy Service (free, sub-20-cm globally), Search and Rescue, OS-NMA (authenticated, anti-spoofing).
The BeiDou Satellite System
BeiDou (BDS-3) is China's global GNSS, fully operational since July 2020. 35 satellites — uniquely combining 24 MEO satellites with 3 geostationary and 3 inclined geosynchronous orbit (IGSO) satellites. Civilian accuracy ~3–5 m globally, 2 m in Asia. Distinguishing features include a regional short-message service over GEO satellites and B1C signal interoperability with GPS L1. The newest operational global GNSS, with the most modern signal design.
GPS IIIF Explained
GPS IIIF (Follow-on) is the next generation of GPS satellites, in development by Lockheed Martin under a 2018 contract. 22 satellites planned; first launch expected ~2027-2028 after multiple delays. Builds on GPS III (currently being launched, 2018-present) with new features: Regional Military Protection (RMP) high-power M-code beams over specific theaters, search-and-rescue payload compatible with COSPAS-SARSAT, laser retroreflectors for Satellite Laser Ranging precise tracking, longer 15-year design life, improved atomic clocks. Civil signals unchanged from GPS III (L1 C/A, L1C, L2C, L5). Ground operational control via the long-delayed OCX program. The article covers the technology, the program history, and the strategic context.
GPS Time vs UTC
GPS Time is a continuous atomic timescale that started on January 6, 1980 at 00:00:00 UTC and has no leap seconds. As of 2026, GPS Time - UTC = 18 seconds (the number of leap seconds added since the epoch). GPS satellites broadcast GPS Time plus the GPS Time - UTC offset; receivers compute UTC by applying the offset. GPS Time has a 10-bit Week Number field that rolls over every 1024 weeks (~19.6 years); the first rollover was August 21, 1999, the second April 6, 2019, the third around November 2038. Different GNSS systems use different timescales — GPS Time, GLONASS Time (UTC-coupled), Galileo System Time, BeiDou Time — each with specific offsets to UTC. The article covers the relationships, the rollovers, the broadcast format, and cross-GNSS timescale comparisons.
Dual-Frequency GPS Explained
Dual-frequency GPS uses two satellite frequencies (typically L1 at 1575.42 MHz and L5 at 1176.45 MHz) to directly compute ionospheric delay rather than estimating it from a model. The result: significantly improved accuracy — sub-meter under open sky vs 3-5 meters for single-frequency. The Broadcom BCM47755 chipset in 2018 brought dual-frequency to smartphones (Xiaomi Mi 8 was first; Pixel 4 in 2019; iPhone 14 Pro in 2022). The L5 signal and its Galileo E5a, BeiDou B2a, and QZSS L5 equivalents are at the same 1176.45 MHz frequency, simplifying multi-GNSS dual-frequency receivers. The article covers the physics, the smartphone adoption timeline, accuracy improvements, and remaining limitations.
Geocoding
What Is Geocoding?
Geocoding explained — forward (address → coordinate), reverse (coordinate → address), the major providers, match-code confidence, and accuracy by reference data.
Forward vs Reverse Geocoding
Forward geocoding converts an address to coordinates; reverse geocoding converts coordinates to an address. Forward is harder (addresses are messy, ambiguous, parsing-intensive); reverse is more deterministic but has its own edge cases (border points, multiple addresses sharing a location, ocean coordinates). The article covers the asymmetry, the practical use cases for each direction, and the implementation considerations.
How Geocoding Works
Under the hood of a geocoding system: address parsing, standardisation, gazetteer indexing, candidate matching, street-segment interpolation, and confidence ranking. The article walks through the internal pipeline that converts a messy text query into a ranked list of coordinate candidates — and the algorithms (libpostal, Elasticsearch / spatial indices, TIGER address ranges) that make it work at scale.
Address Standardization
Address standardization is the cleaning step that turns messy human input into a canonical form a geocoder can look up. USPS Publication 28 defines US conventions: street-suffix abbreviations, directional prefixes, secondary-unit handling, ZIP+4 normalization. The Universal Postal Union defines international conventions. The article covers the standardization rules, the libpostal library, common pitfalls (over-aggressive normalization, multilingual edge cases), and why standardization is upstream of every geocoder.
Geocoding Accuracy Levels
Geocoders return results at tiered accuracy levels: rooftop (~5 m, building footprint known), parcel (~10–20 m, property lot known), street-segment interpolated (~20–50 m, address-range estimate), street-level (~50–100 m, street known), postal-code centroid (~1–5 km), and broader administrative tiers. The article maps the tiers to per-provider taxonomies (Mapbox confidence, Google location_type, HERE matchLevel), explains how to filter and weight by tier, and the data-quality issues that drive tier degradation.
Major Geocoding Services Compared
A comparison of the major civilian geocoding services in 2026: Google Maps Platform, Mapbox, HERE, OpenCage, Nominatim, LocationIQ, SmartyStreets, Esri / ArcGIS, Pelias, TomTom, AWS Location Service. Compared on coverage quality, accuracy tiers supported, pricing, ToS quirks, API quality, self-hostability. No single winner — each has a niche where it's the right answer.
Nominatim Explained
Nominatim is the open-source geocoder for OpenStreetMap data — both the forward (address-to-coordinate) and reverse (coordinate-to-address) directions. Built initially by Brian Quinion around 2010, maintained by Sarah Hoffmann at the OSMF since the late 2010s; the codebase migrated from C/PHP to Python in 2024. Architecture: PostgreSQL/PostGIS database with Nominatim-specific tables, REST API on top. The OSM Foundation runs a free public service at nominatim.openstreetmap.org with a strict 1 req/sec usage policy; production use requires self-hosting or using a commercial Nominatim-based service (Geoapify, LocationIQ, MapTiler). Self-hosting needs ~64 GB RAM and ~1.5 TB SSD for the global database. Licensed GPL-2.0; OSM data ODbL.
Mapbox vs Google vs HERE Geocoding
The three dominant commercial geocoders compared on price, coverage, accuracy, terms, and feature set. Google Maps Geocoding API: premium pricing (~$5+ per 1,000 requests), broadest coverage, best POI database, strict caching rules. Mapbox Geocoding API v6 (the Coordinately default): ~$0.50 per 1,000 above 100k free tier, no-retention ToS, strong autocomplete, confidence bands. HERE Geocoding: mid-tier pricing, especially strong in Europe, generous developer tier (250k/month free), batch endpoints, strong logistics/routing integration. Each has specific strengths: Google for global accuracy, Mapbox for cost-effective web apps, HERE for European and logistics use. Open-source Nominatim is the budget alternative.
Geocoding Rate Limits and Caching
Operational considerations for production geocoding. Rate limits vary by provider: Mapbox v6 ~600 req/min default; Google ~50 QPS for free, scaling with billing; HERE 250k/month free with no per-second limit; Nominatim public service 1 req/sec. Caching restrictions vary: Mapbox no-retention per ToS (covered in Coordinately arch §19.2); Google allows 30-day caching; HERE typically 30-day; OSM Nominatim encourages caching. Batch endpoints exist for HERE and Mapbox v6. Session-token billing (Google) charges per session rather than per request for autocomplete UX. Operational patterns: coarse-then-fine queries, type-ahead with full geocoding only on selection, multi-provider failover, circuit breakers, PII handling, monitoring.
History
A History of Latitude and Longitude
A 2,300-year history of latitude and longitude — Eratosthenes (240 BCE), Ptolemy (150 CE), Harrison's chronometer (1759), the 1884 Greenwich vote, and modern GPS.
Eratosthenes and the Measurement of Earth's Circumference
Around 240 BC, the Greek astronomer Eratosthenes of Cyrene — chief librarian of the Library of Alexandria — measured Earth's circumference by comparing the sun's angle at noon between two cities (Syene and Alexandria). His result of 250,000 stadia translates to ~39,400 km — within 1.6% of the modern value of 40,075 km. The article tells the full story, the modern reproductions, the stadion-length controversy, and what 'this was the first quantitative measurement of a planetary parameter' means for the history of science.
The Longitude Problem
From the 1500s to the 1700s, navigators could measure latitude at sea but not longitude. The lack of accurate longitude killed sailors, sank ships, and stranded fortunes. The British Parliament's 1714 Longitude Act offered £20,000 (~£3 million today) for a solution. John Harrison's marine chronometers (H1 through H4, 1735–1761) eventually solved it. The article traces the problem, the alternative methods (lunar distance, Galileo's moons), and why the chronometer approach won.
John Harrison and the Marine Chronometer
John Harrison (1693–1776) was a self-taught Yorkshire carpenter who solved the longitude problem with four successive marine chronometers — H1 (1735), H2 (1741), H3 (1759), and H4 (1761) — and a final pocket-watch design H5 (1772). His H4 chronometer kept time to ~5 seconds across an 81-day voyage from England to Jamaica in 1761–62. The article covers his life, the H1–H5 designs, Larcum Kendall's K1 (used by Captain Cook), Harrison's 50-year battle with the Board of Longitude, and the chronometers' legacy.
Celestial Navigation
Celestial navigation determines position at sea by measuring the angles of the sun, moon, planets, and 57 navigational stars above the horizon — combined with the exact time and printed tables (Nautical Almanac). The technique was the primary navigation method from roughly 1750 until GPS in the 1990s. Modern ships still carry sextants as the GPS backup of last resort. The article covers the method, the instruments (sextant + chronometer + almanac), the math (sight reduction), the 57 navigational stars, the Marcq St Hilaire intercept method, and modern relevance.
Dead Reckoning
Dead reckoning is the method of estimating current position by applying course, speed, and elapsed time to a known previous position. Used continuously from ancient Phoenician seafaring through modern GPS-failure backup, it's the most fundamental navigation technique — every modern position-fixing method builds on it. The article covers the math, the historical instruments (compass, chip log, traverse board), the error accumulation, and modern relevance (inertial navigation, GPS backup, autopilot).
The 1884 International Meridian Conference
In October 1884, delegates from 25 nations met in Washington DC and resolved that the meridian through the Royal Observatory at Greenwich would be the world's prime meridian for longitude and timekeeping. The conference passed seven resolutions across two weeks, with the key Greenwich vote going 22 in favor, 1 against (San Domingo), and 2 abstentions (France, Brazil). France did not officially adopt Greenwich until 1911 and kept its own time-reckoning practice until 1978. The article covers the pre-conference state, Sandford Fleming's railway-time advocacy, the resolutions, the vote, the French abstention, and modern relevance.
The Sextant Explained
The sextant is the optical instrument that measures the angle between a celestial body and the horizon, with arcminute precision, from a moving ship. Invented independently by John Hadley (1731) and Thomas Godfrey (1730), it replaced earlier cross-staffs, mariner's astrolabes, and Hadley quadrants. The article covers the optical principle (double reflection), the components (frame, index arm, mirrors, telescope, drum/vernier, shades), the five sight corrections (index, dip, refraction, parallax, semi-diameter), the predecessor instruments, modern manufacturers, and the technique of bringing a star down to the horizon.
A History of GPS
The Global Positioning System evolved from 1957 Sputnik-tracking research, through the U.S. Navy's Transit Doppler system (1964 operational), to the 1973 DoD-mandated Navstar GPS program. First Block I satellite launched 1978; 24-satellite constellation fully operational 1994; Selective Availability turned off May 1, 2000, improving civilian accuracy from ~100 m to ~10 m overnight. GPS III generation began launching in 2018; GPS IIIF is in development. The article traces the technical milestones, the political decisions, the parallel GNSS programs (GLONASS, Galileo, BeiDou), and the present state of the constellation.
Time and Time Zones
Time Zones Explained
Time zones explained — 38 distinct UTC offsets (not 24), India +05:30 and Nepal +05:45, the IANA tz database, and why offsets follow politics, not longitude lines.
UTC Explained
UTC (Coordinated Universal Time) is the world's primary time reference, defined by ITU-R Recommendation TF.460-6 and coordinated by the Bureau International des Poids et Mesures (BIPM). UTC combines two ingredients: International Atomic Time (TAI), computed from ~400 atomic clocks at ~80 laboratories in ~50 countries, and occasional leap seconds that keep UTC within 0.9 seconds of UT1 (Earth-rotation time). Since 1972, 27 leap seconds have been added; the most recent was December 31, 2016. The 2022 CGPM Resolution 4 mandated suspending leap-second additions by 2035. The article covers the BIPM coordination, the leap-second mechanism, UTC's relationship to TAI/UT1/GPS-time, and the planned 2035 transition.
GMT vs UTC
Greenwich Mean Time (GMT) and Coordinated Universal Time (UTC) look identical in practice — synchronized to within 0.9 seconds — but they're technically distinct. GMT is based on Earth rotation (essentially UT1 for civil use); UTC is an atomic timescale with leap-second corrections. GMT has been in use since 1675; UTC was introduced in 1972. GMT remains the UK's legal civil-time reference (per the Interpretation Act 1978); UTC is the world technical standard for science, telecoms, internet protocols, and GPS. The article covers the historical origin of GMT, the metrological origin of UTC, the technical difference, the legal status, and when the distinction actually matters.
Leap Seconds Explained
Leap seconds are one-second adjustments inserted into UTC to keep it within 0.9 seconds of UT1 (the astronomical timescale based on Earth's actual rotation). The IERS announces them via Bulletin C with about six months of forward notice; insertion happens at the end of UTC June 30 or December 31, with the unusual second 23:59:60 appearing in the timestream. Since 1972, 27 leap seconds have been added — the most recent on December 31, 2016. The 2012 leap second triggered a Linux-kernel bug that brought down Reddit, LinkedIn, Mozilla, Qantas reservations, and others. CGPM Resolution 4 (2022) mandates that leap-second additions be suspended by 2035.
The IANA Time Zone Database
The IANA Time Zone Database (the 'tz database' or 'Olson database') is the authoritative reference for civil time zones worldwide. Created in 1986 by Arthur David Olson at the U.S. National Institutes of Health, coordinated by Paul Eggert at UCLA since around 2005, and formally hosted by IANA since 2011 (after the Astrolabe lawsuit). It tracks ~600 named zones in Region/City format, releases several updates per year, and is the time-zone authority for every Unix-like operating system, the Java JVM, Python, Node.js, .NET, Go, Ruby, and most major databases. The article covers the history, naming conventions, source-file structure, release process, software integration, and the Windows-vs-IANA mapping problem.
Daylight Saving Time Explained
Daylight Saving Time (DST) is the practice of advancing civil clocks by one hour during warmer months to shift daylight from morning to evening. First nationally adopted by Germany and Austria-Hungary on April 30, 1916 as a wartime fuel-saving measure. About 70 countries observe DST as of 2026; roughly 120 don't. The transition asymmetry creates two software headaches: some local times don't exist (spring forward), some occur twice (fall back). The EU voted in 2018 to abolish mandatory DST, with implementation delayed; the U.S. Sunshine Protection Act passed the Senate in 2022 but stalled. Health impacts include documented spikes in heart attacks and traffic accidents around the spring transition.
The International Date Line
The International Date Line (IDL) is the conventional line at approximately ±180° meridian where the calendar date changes. It is not a formal treaty — it emerged from the 1884 Conference's ±180° convention combined with later sovereign decisions. The line zigzags around Pacific island nations rather than running along the exact 180° meridian. Kiribati moved the line eastward in 1995 (eastern Line Islands jumped from UTC-10 to UTC+14, the most extreme positive offset on Earth). Samoa jumped to the western side on December 30, 2011 — a date that did not exist in Samoa. The article covers the conventions, the historical zigzags, the travel rules, the Magellan first-circumnavigation date-loss incident, and the software implications.
ISO 8601 Date and Time Format
ISO 8601 is the international standard for unambiguous date and time representation. Published in 1988 and most recently updated as ISO 8601-1:2019 and ISO 8601-2:2019. The canonical extended format is YYYY-MM-DDThh:mm:ss with a UTC offset suffix (Z for UTC or ±hh:mm). The standard also defines basic compact formats, ordinal dates (YYYY-DDD), week dates (YYYY-Www-D), durations (P3Y6M4DT12H30M5S), intervals (start/end), and repeating intervals (Rn/start/end). RFC 3339 (2002) is a stricter IETF subset used by most internet protocols. The article covers the formats, the RFC 3339 differences, common encodings (XML, JSON, HTML5), the sortable property, and frequent pitfalls.
The Prime Meridian
The prime meridian explained — 0° longitude through Greenwich since 1884, the IERS Reference Meridian 102 m east of the Airy Transit, and the 8 countries crossed.
How Many Time Zones Are There?
The textbook answer is 24. The real answer is about 38-40 distinct UTC offsets currently in use, because of half-hour offsets (India, Iran, Newfoundland), quarter-hour offsets (Nepal, Chatham, Eucla), and the extended range from UTC-12 to UTC+14 — a total span of 26 hours. The IANA time-zone database lists about 600 named zones to handle historical and political complexity. This support covers the gap between the textbook 24 and the operational reality.
How Time Zones Were Created
Before railways, every town set its own clock by local solar time. Coordinating a train schedule across a continent was impossible without a common reference. This support traces how time zones came into existence — from 1840s British railway time, through Charles Dowd's 1869 US schoolteacher proposal, Sandford Fleming's 1879 international scheme, the 18 November 1883 'Day of Two Noons' when US railways unified, the 1884 International Meridian Conference, and the patchwork country-by-country adoption that followed.
Why Greenwich Is the Prime Meridian
The prime meridian could have run through Paris, Pulkovo, Washington, or any other longitude — there is no physical reason for 0° to pass through Greenwich. The 1884 International Meridian Conference chose Greenwich primarily because by then about 72% of world shipping tonnage already used Greenwich-referenced charts. This support covers the pre-1884 prime-meridian landscape, the competing candidates at the conference, the practical case that won the day, France's 27-year holdout, and the modern IERS Reference Meridian.
Time Zones Without Hour Offsets
About a dozen UTC offsets are not whole-hour multiples — they are 30 minutes, 45 minutes, or even 15 minutes off the simple hourly pattern. India's UTC+5:30 covers about 1.4 billion people, the largest non-hour zone by population. Iran, Afghanistan, Myanmar, Newfoundland, and several Australian zones use 30-minute offsets; Nepal, Eucla, and the Chatham Islands use 45-minute offsets. This support covers each non-hour zone, why each country chose it, and the software implications.
Countries with Multiple Time Zones
About 20 countries use more than one time zone. Russia uses 11 — the most of any country — spanning UTC+2 to UTC+12. France technically uses 12 counting overseas territories, but only one in mainland France. The US has 6, Canada has 6, Australia 3-5, Mexico 4, Brazil 4. China and India each cover huge east-west extents on a single zone. This support covers each major multi-zone country, why each chose its zoning, and the operational consequences.
Magnetic and Compass
Magnetic Declination Explained
Magnetic declination explained — the angle between compass north and true north, ranging −17° to +20° in CONUS, with the WMM 2025 model and the drifting magnetic pole.
Magnetic North vs True North
Practical reference for converting between magnetic, true, and grid north — the three norths that appear on every printed topographic map. True north is the direction along a meridian to the geographic North Pole. Magnetic north is where a compass needle points. Grid north is along the y-axis of the map's projection grid (UTM, BNG, etc.). The article covers the geometric relationships, the conversion arithmetic with worked examples, the 'east is least, west is best' mnemonic, the topographic-map declination diagram, and the operational practice (USGS, Ordnance Survey, marine charts, hiking compasses with adjustable bezel, smartphone APIs).
The World Magnetic Model
The World Magnetic Model (WMM) is the standard mathematical representation of Earth's main magnetic field, jointly produced by NOAA NCEI and the British Geological Survey since 1990. Updated every 5 years on a regular cycle (current WMM2025 released late 2024, valid 2025–2030); an out-of-cycle update was issued in February 2019 due to faster-than-expected magnetic-pole drift. The model uses spherical-harmonic expansion to degree 12, fitted to satellite measurements (Swarm since 2013, CHAMP 2000–2010, Ørsted 1999–2014) and ~150 ground observatories. The WMM is mandated by U.S. DoD spec MIL-STD-3034, used by NATO, the FAA, and every smartphone with a compass. The article covers the mathematics, the data sources, the production cycle, the 2019 update, and the comparison with IGRF.
Aviation Runway Numbering
Aviation runway numbers are based on the runway's magnetic heading, rounded to the nearest 10 degrees and divided by 10. A runway numbered '22' has a magnetic heading of approximately 220°; the opposite-direction designation is '04' (040° = 220° - 180°, rounded). Parallel runways get L/C/R suffixes. The convention is mandated by ICAO Annex 14 and was adopted globally in the mid-20th century when pilots relied on magnetic compasses. As magnetic declination drifts, runways are periodically renumbered to keep the magnetic-heading correspondence accurate — notable examples include Tampa International (2011), Fairbanks (2008), and several Alaskan airports. Iceland is the only ICAO member using true-north runway numbering. The article covers the standard, the process, the notable renumbering events, and the slow movement toward true-north designation.
The Earth as a Magnet
Earth's magnetic field is generated by a self-sustaining geodynamo in the liquid outer core — a ~2,260 km thick shell of molten iron-nickel, between 2,890 and 5,150 km depth, that convects under heat flow from the inner core and Earth's rotation. The field is approximately a dipole tilted ~11° from the rotation axis, with surface intensities of 25,000–65,000 nanoteslas. William Gilbert proposed Earth was a giant magnet in De Magnete (1600); modern dynamo theory developed in the 20th century. The South Atlantic Anomaly is a region of weakened field over South America, slowly growing. Pole reversals happen on geological timescales (last 780,000 years ago, average interval ~200,000 years). The article covers the geodynamo, the field structure, comparisons with other planets, the South Atlantic Anomaly, reversal history, and the magnetosphere's role in protecting life.
Magnetic Inclination and Intensity
Declination describes the horizontal direction of Earth's magnetic field, but the field is a three-dimensional vector with two other quantities: inclination (the angle below horizontal, also called magnetic dip) and intensity (the total field strength, typically 25,000–65,000 nanoteslas at the surface). Inclination is 0° at the magnetic equator and ±90° at the magnetic poles; intensity is weakest at the equator and strongest at the poles. Robert Norman invented the dip needle in 1581. The aclinic line (where inclination = 0) is the magnetic equator. Compass needles are zone-balanced because of latitude-varying inclination. The article covers the geometric setup, the dip needle, mapping conventions, compass balancing, and the unit history (gauss → nanotesla → SI).
Geomagnetic Storms and Navigation
Geomagnetic storms are temporary disturbances in Earth's magnetosphere caused by solar events — coronal mass ejections (CMEs) and high-speed solar wind streams. The Carrington Event of September 1859 was the most intense in recorded history, disabling global telegraph systems. The March 1989 Quebec blackout cut power to 6 million people for nine hours. NOAA's Space Weather Prediction Center grades storms on G1–G5 (geomagnetic), R1–R5 (radio blackout), and S1–S5 (solar radiation) scales. Modern storms cause GPS position errors of meters during major events, polar aviation route disruptions, and power-grid GIC (geomagnetically induced currents). The article covers the solar origin, the historical events, the operational scales, the modern impacts, and the monitoring infrastructure.
Magnetic Compasses: History and Types
The compass evolved over 2,000 years from Chinese lodestone fish around 200 BC to today's solid-state smartphone magnetometers. Petrus Peregrinus described the first pivoted dry compass in his 1269 Epistola de Magnete. Liquid-damped marine compasses dominated the 18th–20th centuries. The Anschütz-Kaempfe gyrocompass (1908) freed ships from magnetic dependence. Modern compass types include marine (gimbal-mounted, liquid-damped), aircraft (vertical card), hiking (Silva-style baseplate), military lensatic (M-1950), gyrocompass (mechanical, true north), fluxgate (electronic), and smartphone magnetoresistive sensors. Each has different precision, environmental requirements, and use cases.
Coordinate Format Deep Dives
Plus Codes Explained
Plus Codes (formally Open Location Code) is a Google-developed open-source location encoding that gives any place on Earth a short alphanumeric identifier. Created in 2014 and released under Apache 2.0 license. A full 10-character plus code (e.g., 8FVC9G8F+6X) identifies a ~14 m × 14 m square; the '+' separator always appears between the 8th and 9th characters. Codes can be shortened in known local context. The 20-character alphabet (23456789CFGHJMPQRVWX) excludes easily-confused letters. Coverage is global, including ocean. Google Maps has displayed Plus Codes since 2018, and the system has been adopted as an addressing standard in parts of India and several African countries.
what3words Explained
what3words is a proprietary location-encoding system that divides Earth's surface into ~57 trillion 3 m × 3 m squares and assigns each a unique three-word identifier from a curated ~40,000-word dictionary. Founded in London in 2013 by Chris Sheldrick. Adopted by UK Police, RNLI, Mongolia Post Office, Mercedes-Benz, and others. The system is licensed proprietary — algorithm and dictionary are not open. Independent security researchers (Andrew Tierney 2019, Aaron Toponce 2021) have documented similar-sounding-word ambiguities that can produce misleading addresses. The article covers the algorithm sketch, adoption, real-world use, documented criticisms, and comparison with the open Plus Codes alternative.
Geohash Explained
Geohash is a public-domain location encoding invented by Gustavo Niemeyer in February 2008. It interleaves latitude and longitude bits and base32-encodes the result, producing strings like 'dr5ru6gh' that identify rectangular cells of decreasing size as length grows. The 32-character alphabet (0-9 b-z excluding a/i/l/o) is chosen to avoid confusion. Key property: lexicographic prefix containment — a shorter geohash contains all longer geohashes that start with it. Used as the geo-index primitive in Elasticsearch, MongoDB, Redis, PostGIS, Solr, and many spatial databases. The article covers the algorithm, precision table, the meridian/antimeridian edge cases, and the database-integration patterns.
Elevation and Vertical Datums
Mean Sea Level Explained
Mean Sea Level (MSL) is the average height of the ocean surface over a defined period — typically a 19-year average to span the 18.6-year lunar nodal cycle. Tide gauges and satellite altimetry are the primary measurement methods. MSL is location-specific: the Pacific is ~80 cm higher than the Atlantic across the Panama Canal due to currents and density differences. The geoid approximates MSL but they differ by sea surface topography (up to ±1 m). Modern satellite altimetry (Jason and Sentinel-6) measures global MSL change at ~3.4 mm/year. The article covers measurement, the lunar nodal cycle, chart datums (MLLW, MLW, LAT), sea surface topography, sea level rise, and the Mount Everest height controversy.
Vertical Datums Explained
A vertical datum is a reference surface for measuring elevations. Two fundamental types: ellipsoidal (height above a reference ellipsoid, used by GPS) and orthometric/gravimetric (height above the geoid, used in engineering and water-flow contexts). The North American Vertical Datum of 1988 (NAVD88) is the current US standard, anchored at the Father Point tide gauge near Rimouski, Quebec, and realized through ~70,000 km of leveling. NAVD88 is being replaced by NAPGD2022, a purely gravimetric datum from the GRAV-D program. The European Vertical Reference Frame (EVRF, multiple realizations) is anchored at the Normaal Amsterdams Peil. EGM2008 and EGM2020 are the global gravity-field models. The article covers the datum types, major regional datums, geoid models, and the long-running effort toward a unified International Height Reference System.
Digital Elevation Models Explained
A Digital Elevation Model (DEM) is a gridded representation of terrain elevation — typically a raster file where each pixel stores the elevation at that location. Three common terms are distinguished: DEM (generic), DSM (Digital Surface Model, top-of-surface including trees and buildings), and DTM (Digital Terrain Model, bare ground). The major global DEMs are SRTM (NASA Shuttle mission, 2000, 30 m), ASTER GDEM (NASA/METI, 30 m, latest v3 2019), ALOS World 3D (JAXA, 30 m), Copernicus DEM (ESA, 30 m global / 10 m European), and TanDEM-X (DLR, 12 m commercial). National lidar programs (USGS 3DEP, UK EA) provide sub-meter resolution over their countries. The article covers terminology, the major datasets, resolution vs accuracy, common uses, and file formats.
Contour Lines and Intervals
Contour lines connect points of equal elevation on a topographic map, allowing 3D terrain shape to be read from a 2D image. Charles Hutton drew the first land contour map at Schiehallion in 1774; Dutch nautical charts had used isobaths (depth contours) since Pieter Bruinsz in 1584. Standard intervals vary by scale: USGS 1:24,000 topo maps use 10/20/40 ft depending on terrain; UK Ordnance Survey 1:25,000 Explorer maps use 5 m. Index contours (every fifth, drawn thicker with labels) help reading. V-shapes pointing upstream indicate valleys; concentric circles indicate peaks or pits; closely spaced lines indicate steep slopes. Modern digital contours are generated from DEMs using algorithms like marching squares.
Elevation vs Altitude vs Height
Three related terms used inconsistently across surveying, aviation, and everyday speech. Elevation conventionally refers to a terrestrial feature's vertical distance above mean sea level (orthometric height). Altitude is ambiguous — it has multiple specific meanings in aviation (indicated, true, absolute, pressure, density, geopotential), can also mean ellipsoidal height in geodesy, and means angle above horizon in astronomy. Height is the most generic term, ambiguous without context. The article distinguishes the meanings, explains why aviation has six altitude types, covers the geodetic ellipsoidal vs orthometric vs geopotential distinctions, and notes the linguistic origins.
Lidar Explained
Lidar (Light Detection and Ranging) is a laser-based ranging technology that produces 3D point clouds of the environment. Originally a portmanteau of light and radar in the 1960s, now widely used as a word. Three platform types: airborne (ALS, the most common for terrain mapping), terrestrial (TLS, mm-precision close-range), and spaceborne (ICESat-2 ATLAS, GEDI on ISS). Typical accuracy: ±5–15 cm vertical for airborne, mm-level for terrestrial. The ASPRS LAS format is the standard data container. Multiple returns (first, last, intermediate) enable separating canopy from bare ground. Applications include topographic mapping (USGS 3DEP), forestry, autonomous vehicles, building information modeling, and archaeology (revealing features under canopy).
Bathymetry Explained
Bathymetry is the measurement of water depth, especially of seafloors. The Greek roots bathys (deep) + metria (measurement) gave the modern term. Methods evolved from lead lines (ancient) to single-beam echo sounding (1920s) to modern multibeam sonar (1960s+) and satellite-derived bathymetry (2000s+). GEBCO (General Bathymetric Chart of the Oceans, ~450 m resolution global) is the standard product, with the Seabed 2030 program (Nippon Foundation–GEBCO partnership) targeting complete ocean-floor mapping by 2030. Coverage was ~25% complete at MBES resolution in 2024; growing through national surveys, crowdsourced bathymetry from ships of opportunity, and improved processing. The IHO S-44 standards define accuracy requirements (Special, Order 1a, 1b, 2).
Isostasy and Post-Glacial Rebound
Isostasy is the gravitational equilibrium between Earth's crust and the underlying mantle. Two classical theories: Airy isostasy (variable crustal thickness floating in denser mantle) and Pratt isostasy (variable crustal density at constant thickness); modern understanding combines both plus elastic flexural rigidity. Post-Glacial Rebound (PGR), also called Glacial Isostatic Adjustment (GIA), is the slow ongoing rise of land that was depressed by ice-age glaciers. Fennoscandia rises at up to 9 mm/year; the Hudson Bay region at up to 12 mm/year. The process started ~20,000 years ago and continues today. The article covers the physics, the major regions, the connection to vertical datum drift, and the importance for isolating climate-driven sea-level rise.
Sea Level Rise Explained
Global mean sea level has risen approximately 21-24 cm since 1880, with the rate accelerating from ~1.4 mm/year in the early 20th century to ~3.4 mm/year today (NASA Sentinel-6 measurement). Four major components drive it: thermal expansion (~40% of recent), mountain glacier melt (~20%), Greenland Ice Sheet melt (~15%), and Antarctic Ice Sheet melt (~5%); land water storage changes contribute small variable amounts. IPCC AR6 projects an additional 0.32-1.01 m by 2100 depending on emissions scenario. Local rates vary substantially due to vertical land motion (Jakarta +50 mm/year; Stockholm -3 mm/year). The article covers measurement, components, projections, regional variation, and adaptation.
Tide and Tidal Datum Explained
Tides are the periodic rise and fall of sea level caused primarily by the gravitational attraction of the Moon and Sun. The dominant constituent is M2 (lunar semidiurnal, 12.42-hour period); other major constituents include S2 (solar semidiurnal), O1 and K1 (diurnal), and N2. Tide types vary by location: semidiurnal (two highs/lows per day, US East Coast), diurnal (one of each, Gulf of Mexico), and mixed semidiurnal (two of each but unequal, US Pacific Coast). The Bay of Fundy has the world's largest tidal range at ~16 m. Tidal datums (MHHW, MHW, MTL, MSL, MLW, MLLW, LAT, HAT) are statistical references derived from 19-year averaging. Modern tide prediction uses harmonic analysis of 396 constituents.
Barometric Altimetry Explained
A barometric altimeter is a precise barometer calibrated to display altitude. It converts atmospheric pressure to height using the International Standard Atmosphere (ISA): 1013.25 hPa / 29.92 inHg at sea level with a 6.5°C/km lapse rate. Three altimeter setting conventions: QNH (gives elevation above MSL on the ground), QFE (gives 0 on the runway), QNE (the standard 1013.25/29.92 used above transition altitude for flight levels). Errors arise from non-standard temperature and pressure: 'from high to low, look out below'. Lucien Vidi invented the aneroid barometer in 1844; aviation use developed in the early 20th century. Modern aircraft typically have a primary barometric altimeter plus radar altimeter and GPS-based geometric altitude.
Web Mapping and Tile Systems
Slippy Map Tiles Explained
Slippy maps are interactive web maps built from pre-rendered tile images organized in a Z/X/Y pyramid. Z is the zoom level (0 = world, 18-19 typically the most detailed); X is the column (longitudinal); Y is the row (latitudinal). Each zoom level has 2^Z × 2^Z tiles of 256×256 pixels (or 512×512 for retina). The coordinate system is Web Mercator (EPSG:3857), introduced by Google Maps in 2005 and adopted by OpenStreetMap, Mapbox, MapTiler, and most others. Y axis convention differs: XYZ (top-left origin) vs TMS (bottom-left). Bing uses the related quadkey system. The article covers the addressing math, the tile pyramid, URL conventions, and the standard providers.
Vector Tiles Explained
Vector tiles contain encoded vector geometry (points, lines, polygons) rather than pre-rendered raster images. The dominant format is the Mapbox Vector Tile (MVT) specification — Protocol Buffers-encoded files at .pbf or .mvt extensions, typically 5-10× smaller than equivalent raster tiles. Vector tiles are rendered client-side by MapLibre GL JS or Mapbox GL JS using the open Mapbox Style Specification. Advantages: runtime restyling, interactive features, sharper at all zoom levels. Tradeoffs: more client-side compute, larger initial download per tile. PMTiles (Brandon Liu) provides a single-file archive alternative to tile-server architecture. The article covers the format, encoding, rendering, generation, and major use cases.
Tile Servers and Caching Explained
Tile servers deliver map tiles over HTTP. Architectures include pre-generated static tiles served from a CDN (cheap, immutable), on-demand rendering from a database (flexible, expensive at scale), and hybrid render-and-cache approaches. Major tile servers: TileServer GL, Tegola, Martin, pg_tileserv, GeoServer, Mapnik via mod_tile, MapTiler Server. Caching strategies span browser → CDN edge → CDN regional → origin tile server, with hit rates typically 95%+. Pricing varies from ~$0.50-2 per 1,000 requests commercial to free self-hosted with operational cost. The article covers architectures, caching tiers, cost models, and operational considerations including bot scraping and DDoS mitigation.
Mapbox vs MapLibre Explained
Mapbox GL JS was open-source under BSD-3 from its 2014 launch through version 1.13.x. In December 2020, Mapbox released v2.0 under a proprietary license, sparking a community fork. MapLibre GL JS launched December 11, 2020, forked from Mapbox GL JS 1.13. Backed by Apple, AWS, MapTiler, Stadia Maps, Microsoft, Esri and other industry contributors, coordinated through OSGeo. MapLibre is now the open-source standard for vector-tile rendering, with API compatibility for Mapbox GL JS 1.x styles. Mapbox v2+ adds proprietary features (globe rendering, certain advanced 3D) but at the cost of vendor lock-in and required Mapbox account. MapLibre 5 (2024) added globe rendering, closing much of the feature gap. Coordinately uses MapLibre.
OpenStreetMap Explained
OpenStreetMap (OSM) is a free, open, crowd-sourced map of the world founded by Steve Coast in 2004 and maintained by the OpenStreetMap Foundation. Over 10 million registered users; ~100,000 monthly active contributors; ~2 million with at least one edit. The data model has three primitives — nodes (points), ways (lines/polygons from node lists), and relations (groups of elements) — each tagged with key-value attributes. Licensed under ODbL (Open Database License) since 2012, requiring attribution and share-alike for derivative databases. The Planet file is ~150 GB compressed, updated weekly; regional extracts from Geofabrik. Downstream users include Apple Maps, Microsoft Bing, Amazon Location, Tesla, Snapchat, Foursquare, Pokemon Go, and most web-mapping services.
Map Projections for the Web
Web Mercator (EPSG:3857) dominates web mapping because of its alignment with the Z/X/Y tile pyramid, but it isn't always the right choice. For European thematic maps, ETRS89-LAEA-Europe (EPSG:3035) is the EU-recommended equal-area projection. For polar regions, EPSG:3995 (Arctic) and EPSG:3031 (Antarctic) polar stereographic projections are standard. For North American thematic maps, Albers Equal Area (USGS national maps) is common. Library support varies: Mapbox GL JS / MapLibre GL JS historically Mercator-only (globe view in v5+); OpenLayers has best-in-class projection support; Leaflet + Proj4Leaflet supports custom projections; D3.js programmatic projections. The globe-view alternative trades flat-map convenience for projection-free 3D rendering.
Data Standards and Specifications
ISO 19115 Metadata Explained
ISO 19115 is the international standard for geographic-information metadata. Original version published 2003; current ISO 19115-1:2014 (Part 1: Fundamentals); ISO 19115-2 covers imagery and gridded data; ISO 19115-3 specifies XML encoding (2016 replaced the older ISO 19139). Defines a comprehensive metadata schema covering title, abstract, topic categories, spatial/temporal extent, CRS, lineage, data quality, distribution, contact, use constraints. Profiles: INSPIRE (EU directive), NAP (North American Profile), MEDIN (UK marine). Encoded via XML (ISO 19139 or 19115-3) or JSON (DCAT-AP). Major catalogs: GeoNetwork, INSPIRE Geoportal, data.gov. Modern alternatives: STAC (less rich, more practical), DCAT (W3C broader), schema.org Dataset.
ISO 19111 Spatial Reference Systems Explained
ISO 19111 is the international standard defining the conceptual framework for spatial reference systems (SRS). Published 2003; major revisions 2007 and 2019. The 2019 revision added explicit support for dynamic datums (positions that move with plate tectonics) and datum ensembles. Key concepts: Coordinate Reference System (CRS) = Coordinate System + Datum; types include geodetic, projected, vertical, engineering, temporal, compound. Datum types: geodetic, vertical, engineering, temporal. Coordinate operations: conversion (same datum) vs transformation (different datums). The EPSG Registry (maintained by IOGP) is the standard implementation with ~6,000+ codes. WKT 2 (ISO 19162:2019) is the modern text encoding. PROJ 6 (2019) is the reference C library. The article covers the conceptual framework, datum types, dynamic datums, and the PROJ ecosystem.
GeoJSON (RFC 7946) Explained
GeoJSON is the dominant JSON-based format for encoding geographic data on the web. RFC 7946 (IETF, August 2016) standardized the format, replacing the 2008 community specification. Object types: Geometry (Point, LineString, Polygon, MultiPoint, MultiLineString, MultiPolygon, GeometryCollection), Feature (geometry + properties), FeatureCollection. Coordinates are [longitude, latitude] order — important difference from some other formats. RFC 7946 mandates WGS 84 (EPSG:4326) and removed the optional CRS member from the 2008 spec. Polygon orientation: exterior rings counter-clockwise (right-hand rule), interior rings clockwise. Bounding box (bbox) optional member. MIME type: application/geo+json. Used by Leaflet, MapLibre, Mapbox GL JS, OpenLayers, all major web-mapping libraries. Variants: GeoJSON-LD, TopoJSON (Bostock), GeoJSON-Sequence (RFC 8142), FlatGeobuf (binary).
Grid
The Equator
The equator explained — 0° latitude, 40,075 km circumference on WGS-84, the 21 km equatorial bulge, the 13 countries crossed, and the Coriolis-zero zone.
Parallels of Latitude
Parallels of latitude are the family of circles on Earth's surface parallel to the equator. The equator is the longest parallel and the only one that is also a great circle; every other parallel is a small circle whose radius shrinks toward zero at the poles. This support covers the geometry, the five named parallels (equator, two tropics, two polar circles), the variable length of one degree of latitude across the WGS 84 ellipsoid, and the role parallels play in projections, navigation, and climate.
Meridians of Longitude
Meridians of longitude are the family of half-circles from pole to pole on Earth's surface, each contained in a plane that includes the rotation axis. Unlike parallels of latitude, every meridian is a great-circle arc, and every meridian has the same length on the WGS 84 ellipsoid (about 20,004 km). This support covers the geometry, the dramatic cosine variation in 1° of longitude with latitude, how meridians define local solar noon, and the historical meridian-arc measurements that established Earth's shape.
The Tropic of Cancer
The Tropic of Cancer is the northernmost latitude where the Sun reaches zenith at solar noon, currently at +23°26'22" N. Its position equals Earth's axial obliquity exactly — the Sun reaches the Tropic at the June solstice and moves no further north. This support covers the geometry, the geography of the 16 countries the line crosses across three continents, the climate of the descending Hadley-cell branch that turns much of the Tropic into a desert belt, and why precession has shifted the June-solstice Sun out of the constellation that gave the line its name.
The Tropic of Capricorn
The Tropic of Capricorn is the southernmost latitude where the Sun reaches zenith at solar noon, currently at −23°26'22" S. The Sun reaches the Tropic at the December solstice. The line crosses 10 sovereign states across South America, southern Africa, Madagascar, and Australia, passing through three of the planet's largest deserts. Same axial-obliquity definition as the Tropic of Cancer; same precessional drift; same Hadley-cell descending branch creating subtropical highs.
The Arctic Circle
The Arctic Circle is the northernmost parallel where the Sun remains above the horizon for 24 continuous hours at the June solstice and below the horizon for 24 hours at the December solstice. Its latitude is 90° minus Earth's axial obliquity — currently +66°33′38″ N. The line crosses 8 sovereign states, defines the southern boundary of the Arctic in most geographic frameworks, and is the threshold for both the midnight sun and the polar night.
The Antarctic Circle
The Antarctic Circle is the southernmost parallel where the Sun remains above the horizon for 24 continuous hours at the December solstice and below the horizon for 24 hours at the June solstice. Currently at −66°33′38″ S. The only land it crosses is the Antarctic mainland; the Antarctic Peninsula reaches only to about 63°S. James Cook first crossed it on 17 January 1773. The Antarctic Treaty of 1959 governs the continent under a separate framework that does not use the line directly.
The North Pole
The North Pole is the northernmost point on Earth — the intersection of the rotation axis with the surface at latitude exactly +90°. Longitude is undefined there; every meridian converges. The geographic pole lies in the Arctic Ocean over 4,200 m of water, with no land at all. This support disambiguates the four distinct 'north poles' (geographic, magnetic, geomagnetic, pole of inaccessibility), covers polar motion and the IERS Reference Pole, and tracks the disputed history of reaching the pole.
The South Pole
The South Pole is the southernmost point on Earth — the intersection of the rotation axis with the surface, at latitude exactly −90°. Unlike the North Pole, it sits on land: on the East Antarctic Ice Sheet at about 2,835 m elevation, with ~2,700 m of ice beneath. The Amundsen-Scott South Pole Station has been continuously occupied since 1957. Four distinct south poles need to be kept straight: geographic, magnetic, geomagnetic, and the Pole of Inaccessibility.
The Northern and Southern Hemispheres
The two halves of Earth separated by the equator — geometrically symmetric but dramatically asymmetric in geography and demography. About 68% of Earth's land sits in the Northern Hemisphere, and about 88% of humans live there. The two halves run opposite seasons, see different stars, and host opposite-rotating cyclones. This support covers the geometry, the land-sea and population asymmetries, the climatic consequences, and the celestial differences.
The Eastern and Western Hemispheres
The two halves of Earth separated by the prime meridian and antimeridian. Unlike the equator-bisected hemispheres, this split is set by treaty — the 1884 International Meridian Conference chose Greenwich as the longitude origin, and the Eastern/Western Hemisphere convention inherits that choice. The Western Hemisphere contains the Americas; the Eastern contains Europe, Africa, Asia, and Australia. About 85% of Earth's people live in the Eastern Hemisphere.
Why Magnetic North Moves
The Magnetic North Pole has moved from northern Canada (where James Clark Ross found it in 1831) to the Arctic Ocean north of Russia (where the 2025 World Magnetic Model places it) — about 2,800 km in 194 years, with the drift accelerating from roughly 10 km/year in the early 20th century to ~55 km/year in the late 2010s. This support covers the geodynamo physics in Earth's liquid outer core that produces the field, the field's secular variation, the WMM 2019 emergency update prompted by the rapid drift, and the geomagnetic-reversal record.