← TTastro Design & Theory 中文
SKY QUALITY METER · Field Measurements × Light-Pollution Map

Checking 16,074 Night-Sky Brightness Measurements
Against the Light-Pollution Map

Every SQM reading in the Unihedron global dark-sky database (2005–2026), checked point by point against the David Lorenz 2024 theoretical light-pollution map. Do they agree? And does the year of the data matter?

— readings — sites — in the fair-comparison subset 2005–2026 spanning 22 years
— mag
Median residual (measured − theoretical)
near-perfect agreement
—%
Share of residuals within ±0.3 mag
—
Slope of dark-site residual vs. year
(mag/yr, ≈ 0)
—
Multi-year tracked sites
sky brightness essentially stable
01

Where the data comes from

On one side, night-sky background brightness measured in the field by amateur astronomers worldwide with SQM instruments; on the other, the "theoretical value" returned by a light-pollution lookup service. Only by understanding where each comes from can you compare them correctly.

① Measured · Unihedron Dark-Sky Database

SQM readings submitted by users worldwide, in mag/arcsec² (higher means a darker sky). Each record carries fields such as time, latitude/longitude, weather, moon phase and instrument model, going back as far as 2005.

② Theoretical · 2024 Light-Pollution Map

The theoretical value decoded from the public tiles of the David Lorenz 2024 light-pollution atlas. It represents the background brightness under clear, moonless, zenith conditions.

The online endpoint POST /light_pollution requires authenticated login (tkid + token) and cannot be called anonymously. But in essence it merely reads public tiles and runs a deterministic calculation — this report fully reproduces that algorithm, matching the online endpoint bit for bit at the same coordinates. Theoretical values were obtained for 16,040 / 16,074 readings.
02

How well do they agree

The theoretical map describes an ideal night sky, so a fair comparison is needed: keep only the — readings taken zenith-up, under clear skies, with the moon below the horizon and a sensible value. Each point below is one measurement — theoretical on the x-axis, measured on the y-axis; the closer to the diagonal, the better the agreement.

Measured vs. Theoretical · Calibration Scatter
Grey points are individual readings; the gold line is the measured mean binned by theoretical value; the dashed line is the ideal "theoretical = measured" diagonal.
Binned measured mean Theoretical = measured Individual readings (—)
Residual distribution
Measured − theoretical, in mag/arcsec²
Comparison summary

Among the readings in the comparison, the measured mean is —, the theoretical mean —, and the median residual just — mag.

—
|residual| ≤ 0.1
—
|residual| ≤ 0.3
—
|residual| ≤ 0.5

For a global model map this is remarkably close agreement — it confirms the reliability of the theoretical map, and turns the "residual" into a magnifying glass for other factors.

03

Does the year of the data matter

The earliest data is from 2005, while the theoretical map is the 2024 edition. If a location brightened over these years, its older readings should read darker than the 2024 model (positive residual). Looking at residuals year by year — is the trend really there?

Residual by year (measured − 2024 theoretical)
Grey point size ∝ sample count; the orange line is all comparison points, the purple line dark-site points only (theoretical ≥ 20.5).
All comparison points Dark-site points (city auto-stations removed) Zero-residual baseline
For the full sample, the residual-vs-year correlation is only — with a slope of — mag/yr; once the recent city auto-station contamination is removed and only dark-site points remain, the correlation drops to — and the slope to — mag/yr — essentially zero. The year itself produces no systematic bias.
04

Why the "raw trend" is misleading

The real trap is not the passage of time, but that the sites sampled in different years are simply not the same set. Look at the per-year "mean theoretical brightness" and sample count — the last three years have been flooded by a handful of bright city auto-stations.

Mean theoretical SQM and sample count by year
Bars: comparison sample count that year; line: mean theoretical brightness of those sites that year (lower means brighter).
Sample count Mean theoretical SQM

2005–2023 was mostly amateur observation at dark sites (theoretical ≈ 20.5–21); from 2024, hourly-reporting SQM-L units like the one in suburban Cincinnati pull the sample toward the city (theoretical ≈ 18.7). Without controlling for site, this "change in sampling" gets misread as "the sky changing with the year".

05

How a single location changes over the years

The cleanest approach is to stay fixed on a single point. For the — sites that span at least 4 different years, each is fitted for the slope of "measured SQM ~ year": positive = darkening, negative = brightening.

Distribution of per-site annual-change slopes
In mag/yr; near 0 = stable, right = darkening, left = brightening.
The trend is weak, and its direction is mixed

Mean slope —, median — mag/yr — ground measurements do not show the clear "year-on-year brightening" seen from satellites.

—
Darkening
—
Stable
—
Brightening

Sites tracked for the most years

SiteYear rangeYearsSamples Measured meanTheoretical (2024)Mean residualSlope (mag/yr)Trend
06

Where these sites are located

Plotting the — sites on a world map, coloured by theoretical brightness — bright spots (cities) glow red, dark ones (wilderness) turn blue-black, faintly tracing the map of human lights.

Global site distribution · coloured by 2024 theoretical SQM
Equirectangular projection; point size ∝ number of readings at that point.
07

The dark end: zenith SQM doesn't have the final say

The sections above proved the zenith theoretical map is accurate — but it only answers "how dark is the zenith". At the darkest places (theoretical SQM ≳ 21.7, Bortle 1–3), zenith brightness saturates: the zenith readings of classes 1, 2 and 3 are almost identical, and the zenith alone cannot tell them apart. What truly sets the observing class is space — how large the surrounding dark area is, whether light domes sit on the horizon, and whether mountains block those domes below the line of sight.

The zenith saturates, but space can separate them

Median measured zenith · by Bortle
Y-axis SQM 18–22; classes 1/2/3 sit at almost the same height — the zenith has saturated.
Distance to nearest light dome · by Bortle
The same three classes, cleanly separated by distance (km).
The median measured zeniths of classes 1/2/3 (21.72 / 21.60 / 21.58) almost coincide — zenith SQM cannot separate them; but the "distance to the nearest light dome" (105 / 77 / 44 km) separates them cleanly. So classifying the dark end must shift from the zenith to space. By class 4 (where the zenith has clearly brightened), the zenith takes the lead again.

Ground light domes × mountain shielding

The "artificial glow" at the dark end comes from ground light domes near the horizon: light from the surrounding towns is summed by the Walker/Garstang propagation law Σ A(s)/d^2.5 (radius ~120 km), capturing the whole horizon rather than just the nearest city. A further terrain-gating layer is added: elevation is sampled along each azimuth to compute the horizon angle, and the first source beyond >2° is blocked below the horizon by the mountains and excluded from the glow. For the same zenith darkness, a valley backed by mountains can be far darker than open wilderness.

Terrain shielding pushes ground light domes darker
Each point's shift from "raw glow" → "terrain glow"; further right is darker, hover for values.
Raw glow (without terrain) Terrain glow (after mountain shielding)

The more negative the glow index, the darker. Tivoli and Wiruna are open on all sides, so terrain does not change their glow; but at Crestone the 5–16° Sangre de Cristo ridge to the east blocks the Front Range conurbation 100 km away → trimming 56%, actually overtaking the open Wiruna; Kendrick Park is shielded 52% by the 3850 m San Francisco Peaks → rescued from a "3-class dark-pocket drop" back to class 2. This is exactly the layer the zenith SQM cannot see at all, yet which decides the observing class.

View raw data table
BortleEnvironmentMedian measured zenithNearest light domeDark-area share
1Pristine dark sky21.72105 km100%
2Typical dark sky21.6077 km96%
3Rural dark21.5844 km64%
4Rural/suburban21.2918 km18%
Validation pointZenithRaw glowTerrain glowBlocked by mountains
Tivoli · Namib Desert22.00−3.70−3.700%
Wiruna · Australian outback21.96−2.42−2.420%
Kendrick Park · near Flagstaff21.91−2.07−2.3952%
Crestone · backed by a ridge21.85−2.25−2.6156%

Further reading · The full recalibration table and the glow and terrain-gating formulas are in "Dynamic Light-Pollution Rating · Design & Technical Document"; the field analysis of why the dark-end zenith reading fails is in "The Dark-End Divergence".