REVIEW 4 major objections 5 minor 44 references
The nature of the diffuse light near cities detected in nighttime satellite imagery
T0 review · 4 major / 5 minor · reviewed 2026-08-14 · deepseek-v4-flash
Pith's one-line read The diffuse glow around cities in nighttime satellite imagery is real light scattered by the atmosphere, not an instrumental artifact.
desk verdict A convincing multi-sensor correlation that the diffuse glow around cities is real scattered light, but the causal identification needs radiative transfer support and more transparent fitting choices. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The mechanism is atmospheric scattering, specifically the up/down symmetry of Rayleigh and Mie scattering for near-horizontal light. The paper's central identity is that the same atmospheric volume that scatters city light downward to a ground observer also scatters a nearly equal amount upward to space, so satellite radiance in unlit pixels is an upward view of the same skyglow. This operates through logarithmic regressions between ground-based zenith sky brightness and space-based radiance, with residual analysis separating pixels containing only scattered light from pixels mixed with direct emissions.
What would settle it
Measure the diffuse glow in unlit pixels near a city with VIIRS while simultaneously measuring ground-level zenith sky brightness, under known aerosol conditions; if the two are uncorrelated once distance to the city is removed, or if the glow can be reproduced by a surface albedo model with no atmospheric scattering, the central claim collapses.
Extended reading notes
Core claim
The paper shows that the diffuse light around cities in nighttime satellite imagery is mostly skyglow propagated from urban areas, not detector blooming. Because Rayleigh and Mie scattering are nearly up/down symmetric for the near-horizontal light that dominates ground-level skyglow, light scattered upward toward space and light scattered downward to the ground share the same source. This is supported by correlations between ground sky-brightness data and radiance from DMSP, VIIRS Day/Night Band, and two composites of ISS astronaut photographs; the main residual component has a half width of about 0.35 mag per square arcsecond, roughly a 35 percent radiance difference. A second residual component is attributed to pixels that mix direct ground emissions with scattered light, while the diffuse signal itself is most visible in the darkest pixels.
Load-bearing premise
The argument assumes that the unlit pixels used for comparison contain only scattered skyglow, with no direct ground emissions, and that the atmosphere scatters near-horizontal city light about equally upward and downward.
Editorial extensions
If this is right
- A future space-based imaging radiometer could monitor changes in the diffuse artificial skyglow of cities from orbit.
- Dark pixels near cities in VIIRS and similar data can be used to build regional or global sky-brightness maps, complementing purely modeled skyglow products.
- Changes in street-lighting intensity or spectra that alter ground-observed skyglow should appear as corresponding changes in the diffuse glow around cities.
- Satellite radiance in the diffuse regime around cities should be treated as scattered light rather than surface emission, and subtracting a sky-brightness model can isolate direct ground emissions.
- The consistency of residuals across four independent imagers indicates that low-resolution composites mix direct and scattered light, so high-resolution data are preferable for skyglow monitoring.
Reading between the lines
- If the correlation holds globally, archived DMSP composites contain a decades-long record of skyglow changes that predate dedicated ground-based skyglow surveys.
- The up/down symmetry could be inverted: ground-based skyglow measurements may serve as a calibration for upward-directed light flux, connecting satellite radiance to light-pollution inventories.
- Pairing space-derived skyglow with simultaneous aerosol measurements could yield a nighttime aerosol optical thickness proxy, extending daytime remote sensing to dark hours.
- The identified pixel-mixing tail implies that low-resolution nighttime-light products systematically blur the boundary between direct emissions and skyglow, so skyglow-free products must either preserve high resolution or subtract a scattering model.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper investigates the diffuse glow surrounding cities in nighttime satellite imagery, often dismissed as instrumental 'blooming.' The authors compare ground-based zenith sky brightness measurements from the Universidad Complutense Madrid (UCM) SQM survey with satellite radiance from four very different imaging systems: DMSP (5 km), VIIRS DNB (750 m), and two ISS astronaut-photograph HDR composites (240 m and 54 m). They find strong correlations between ground sky brightness and satellite radiance in pixels away from directly lit areas, interpret the diffuse glow as upward-scattered artificial skyglow rather than a sensor artifact, and propose that future space radiometers could monitor artificial skyglow. They also compare VIIRS-derived sky brightness with the World Atlas of artificial sky brightness and discuss other contributors to diffuse light such as surface albedo, airglow, fog, and residual blooming.
Significance. If established, the paper's central claim would change how the diffuse glow in nighttime satellite imagery is interpreted and would open the door to space-based monitoring of artificial skyglow. The study has real strengths: it uses four independent sensors with very different point-spread functions, all showing a consistent correlation with independent ground-based SQM data, which is strong evidence that the glow is a real optical signal rather than a single-sensor artifact. The paper is clearly written, uses public data, and candidly acknowledges that not all diffuse light is skyglow. However, the identification of the glow specifically as atmospheric scattering of city light rests on an unquantified symmetry assumption and on post hoc pixel-selection choices, and the World Atlas comparison is partly circular because the same UCM survey was used in the Atlas calibration. These issues prevent the paper, in its current form, from fully supporting the strong claim in the abstract.
major comments (4)
- [Methods, 'Physical base'] The central inference of the paper depends on the assertion that 'because of the up/down symmetry of Rayleigh and Mie diffusion for light traveling horizontally, it is reasonable to suppose that the same amount of light that is scattered to the ground from the zenith is also scattered towards space from the same volume of atmosphere.' This is an assumption, not a derivation. The relationship between ground zenith sky brightness and space-borne upwelling radiance depends on the angular distribution of city light emission, the aerosol phase function (which is strongly forward-peaked for Mie scattering), the vertical profile of scatterers, and the observation geometry. The paper does not provide a radiative transfer calculation or a quantitative model that connects the two quantities for Madrid's source distribution. Without such a calculation, the observed correlation could also arise from a common distance-to-city gradient or from low-level direct emission in the nominally 'unlit' pixels, so the abstract's claim that the glow is 'a real detection of light scattered by the atmosphere' is not yet established by the evidence presented.
- [Results, Figs. 3-5 and Fig. 6] The linear fits that define the ground-space relationship are made after rejecting bright bins: the caption of Fig. 3 says the fit uses 'the darkest values and rejecting the bins brighter than 19 mag/arcsec2,' and Fig. 4 uses 18.8 mag/arcsec2. These thresholds are post hoc, and the double-Gaussian decomposition of the residuals (Fig. 6) attributes the second component to direct light contamination without an independent test. Because the regression itself defines the expected sky brightness, the residual analysis is partly circular. The authors should show that the correlation and the width of the main residual peak are robust to (a) varying the brightness threshold, (b) excluding the pixels closest to the city center, and (c) including distance to the city as a control variable. Without such checks, the selection of a 'skyglow-only' population is not justified.
- [Results, 'Comparison with the World Sky atlas' and Methods] The comparison with the World Atlas model (Fig. 13) is not an independent validation, because the manuscript itself states in the Methods that 'The UCM sky brightness survey was used for the calibration of the "World Atlas of artificial sky brightness".' The apparent 'perfect fit' between VIIRS-derived values and the Atlas therefore contains a circular component. A stronger test would use a sky-brightness model or ground survey that was not part of the Atlas calibration, or a first-principles radiative transfer model driven by the observed VIIRS radiance field.
- [Results, 'Relationship between sky brightness and diffuse light'] The paper describes the evidence as 'strong correlation' and 'similar dispersion' but does not report the coefficient of determination, fit parameters, uncertainties, or sample sizes for the regressions in Figs. 3-5. These numbers are necessary for the reader to assess the strength of the central correlation and to support the claim that the dispersion is similar across four sensors with very different PSFs. Visual inspection of 2D histograms is not sufficient for a quantitative claim of this kind.
minor comments (5)
- [Throughout] The manuscript refers to the diffuse glow as 'blooming' but later distinguishes residual instrument blooming (Fig. 15); the terminology should be clarified early so that 'blooming' is not used for two different phenomena.
- [Discussion] The statement 'The changes in fit are produced by the different mean composition of the atmosphere when the image was taken' is speculative, since no atmospheric composition data are presented for the imaging epochs.
- [Discussion] The sentence 'We can conclude that the radiance observed in the regime between 0.2-5 nW/cm2/sr is dominated by diffuse light' introduces a specific radiance range without explaining how the endpoints were determined.
- [Discussion] There is a typographical error in the sentence beginning 'Other components such as natural sky glow...': 'in other to extract' should be 'in order to extract'.
- [Methods and Fig. 13] The choice of the months 'closest to the four months second equinox (Aug to Nov)' and the averaging over 2012-2018 should be justified, as seasonal and interannual changes in VIIRS data are known to occur.
Circularity Check
Central ground–space correlation is independent; only the secondary World Atlas validation shares UCM calibration data.
-
other
[Results, 'Comparison with the World Sky atlas' (Figure 13); Methods, 'Physical base', final paragraph.]
"Another way to verify the accuracy of this VIIRS data is to compare it with the predictions of a Sky Brightness model31 based on the Madrid Sky brightness survey and other measurements."
The cited model (Falchi et al. 2016) was calibrated with the same UCM Madrid survey: the Methods state 'The UCM sky brightness survey was used for the calibration of the World Atlas of artificial sky brightness'. The VIIRS values compared in Figure 13 are themselves converted to sky brightness using the UCM-calibrated regression from Figures 3-4. Therefore the reported 'perfect fit' is partly a self-consistency check between two products sharing the same UCM calibration input, not a fully independent external validation. This step is not the load-bearing support for the central claim, which rests on the direct ground-space correlations in Figures 3 and 4; hence the circularity is minor and secondary.
full rationale
The main derivation chain is not circular: UCM ground-based zenith sky brightness measurements are independent of DMSP, VIIRS DNB, and ISS astronaut-photograph radiances, and the paper finds strong correlations across four instruments with different resolutions and spectral responses. The physical argument (up/down symmetry of scattering for near-horizontal light) is an explicitly stated assumption rather than an equation that forces the result. The residual 'double Gaussian' analysis is an interpretation of the scatter, not a fitted parameter renamed as a prediction. The one partially self-referential element is the World Atlas comparison, because both the model and the satellite-derived values share the UCM survey as calibration input; however, that comparison is presented as an additional verification and the central claim that the diffuse glow is real scattered light does not reduce to it. Accordingly, a low score is appropriate.
Assumptions & free parameters
free parameters (6)
- Per-sensor least squares slope and intercept for DMSP-to-SQM prediction
- Per-sensor least squares slope and intercept for VIIRS DNB-to-SQM prediction
- Per-sensor least squares slope and intercept for ISS040 HDR-to-SQM prediction
- Per-sensor least squares slope and intercept for ISS030 HDR-to-SQM prediction
- Brightness rejection thresholds for the fits =
19 mag/arcsec2 (DMSP), 18.8 mag/arcsec2 (ISS)
- Double Gaussian parameters for residual histograms =
HWHM of main component about 0.35 mag/arcsec2
assumptions (6)
- domain assumption Shallow-angle city light is the dominant source of ground skyglow, and Rayleigh and Mie scattering of horizontally traveling light have up/down symmetry, so the same scattered light reaches the ground and space.
- domain assumption Pixels classified as away from installed lighting contain only scattered skyglow and no significant direct light emission.
- ad hoc to paper Madrid's sky brightness, streetlight spectra, and atmospheric conditions are representative of the rest of the world.
- ad hoc to paper The Falchi et al. World Atlas model is a valid external benchmark for VIIRS-derived sky brightness.
- domain assumption Madrid street lighting power consumption did not change significantly between 2009 and 2014.
- domain assumption The relationship between SQM sky brightness and satellite radiance is log-linear over the fitted range.
Cite this review
Pith. "Pith review of The nature of the diffuse light near cities detected in nighttime satellite imagery." pith.science (2026). https://pith.science/paper/6IWA2VTO
@misc{pith2026190805482,
author = {Pith},
title = {Pith review of: The nature of the diffuse light near cities detected in nighttime satellite imagery},
year = {2026},
howpublished = {\url{https://pith.science/paper/6IWA2VTO}},
note = {Machine review of arXiv:1908.05482}
}
read the original abstract
Diffuse glow has been observed around brightly lit cities in nighttime satellite imagery since at least the first publication of large scale maps in the late 1990s. In the literature, this has often been assumed to be an error related to the sensor, and referred to as "blooming", presumably in relation to the effect that can occur when using a CCD to photograph a bright source. Here we show that the effect is not instrumental, but in fact represents a real detection of light scattered by the atmosphere. Data from the Universidad Complutense Madrid sky brightness survey are compared to nighttime imagery from multiple sensors with differing spatial resolutions, and found to be strongly correlated. These results suggest that it should be possible for a future space-based imaging radiometer to monitor changes in the diffuse artificial skyglow of cities.
Figures
Figures from the paper (8 more)
Reference graph
Works this paper leans on
-
[1]
Imhoff, M. L., Lawrence, W. T., Stutzer, D. C. & Elvidge, C. D. A technique for using composite DMSP/OLS “city lights” satellite data to map urban area. Remote Sensing of Environment 61, 361 –370 (1997)
work page 1997
-
[2]
Elvidge, C. D. et al. Radiance calibration of dmsp-ols low-light imaging data of human settlements. Remote Sensing of Environment 68, 77–88 (1999)
work page 1999
- [3]
-
[4]
Levin, N. The impact of seasonal changes on observed nighttime brightness from 2014 to 2015 monthly viirs dnb composites. Remote Sensing of Environment 193, 150–164 (2017)
work page 2017
-
[5]
Publications of the Lick Observatory of the University of California
Bumham, S. Publications of the Lick Observatory of the University of California. No. v. 1 in University of California publications (The University, 1887). URL http://books.google.es/books?id=goM2AQAAIAAJ
-
[6]
Rosebrugh, D. W. Sky-glow from large cities. J Roy Astron Soc Can 29, 79 (1935)
work page 1935
- [7]
-
[8]
Physical behaviour of anthropogenic light propagation into the nocturnal environment
Aube´, M. Physical behaviour of anthropogenic light propagation into the nocturnal environment. Philosophical Transactions of the Royal Society of London B: Biological Sciences 370, 20140117 (2015)
work page 2015
Show all 44 references
-
[9]
Walker, M. F. The california site survey. Publications of the Astronomical Society of the Pacific 672–698 (1970)
1970
-
[10]
Walker, M. F. Light pollution in california and arizona. Publications of the Astronomical Society of the Pacific 508–519 (1973)
1973
-
[11]
C., de Graeve, E
Bertiau, F. C., de Graeve, E. & Treanor, P. J. The artificial night-sky illumination in italy. Vatican Observatory Publications 1, 157–179 (1973)
1973
-
[12]
Hoag, A. A. City sky glow monitoring at Kitt Peak. Publ Astron Soc Pac 85, 503–507 (1973)
1973
-
[13]
Improved scattering formula for calculations of artificial night-sky illumination
Garstang, R. Improved scattering formula for calculations of artificial night-sky illumination. The Observatory 104, 196–197 (1984)
1984
-
[14]
Night-sky brightness at observatories and sites
Garstang, R. Night-sky brightness at observatories and sites. Publications of the Astronomical Society of the Pacific 306–329 (1989)
1989
-
[15]
Cleaver, O. P. Control of coastal lighting in anti-submarine warfare. Tech. Rep., Engineer Board Fort Belvoir VA (1943). URL https://apps.dtic.mil/dtic/tr/fulltext/u2/a954894.pdf. Report No. 746, GN 373. Access date: 2019 May 16
1943
-
[16]
Campaign of sky brightness and extinction measurements using a portable ccd camera
Falchi, F. Campaign of sky brightness and extinction measurements using a portable ccd camera. Monthly Notices of the Royal Astronomical Society 412, 33–48 (2011)
2011
-
[17]
Zamorano, J. et al. Night sky brightness and light pollution in Comunidad de Madrid (2013)
2013
-
[18]
Kyba, C. C. et al. Citizen science provides valuable data for monitoring global night sky luminance. Scientific Reports 3 (2013)
2013
-
[19]
D., Baugh, K
Elvidge, C. D., Baugh, K. E., Zhizhin, M. & Hsu, F.-C. Why viirs data are superior to dmsp for mapping nighttime lights. Proceedings of the Asia-Pacific Advanced Network 35, 62–69 (2013)
2013
-
[20]
Sa´nchez de Miguel, A. et al. Atlas of astronaut photos of earth at night. Astronomy & Geophysics 55, 4–36 (2014)
2014
-
[21]
Kyba, C. et al. High-resolution imagery of earth at night: New sources, opportunities and challenges. Remote Sensing 7 , 1–23 (2014)
2014
-
[22]
Variacio´n espacial, temporal y espectral de la contaminacio´n lum´ınica y sus fuentes: Metodolog´ıa y resultados
Sa´nchez de Miguel, A. Variacio´n espacial, temporal y espectral de la contaminacio´n lum´ınica y sus fuentes: Metodolog´ıa y resultados. Ph.D. thesis, Universidad Complutense de Madrid (2015)
2015
-
[23]
& Schwarz, R
Puschnig, J., Schwope, A., Posch, T. & Schwarz, R. The night sky brightness at potsdam-babelsberg including overcast and moonlit conditions. Journal of Quantitative Spectroscopy and Radiative Transfer 139, 76 –81 (2014)
2014
-
[24]
& Uttenthaler, S
Puschnig, J., Posch, T. & Uttenthaler, S. Night sky photometry and spectroscopy performed at the Vienna University Observatory. Journal of Quantitative Spectroscopy and Radiative Transfer 139, 64 –75 (2014)
2014
-
[25]
Sa´nchez de Miguel, A. et al. Sky quality meter measurements in a colour-changing world. Monthly Notices of the Royal Astronomical Society 467, 2966–2979 (2017). 8/19
2017
-
[26]
de Miguel, A. S. et al. Colour remote sensing of the impact of artificial light at night (i): The potential of the international space station and other dslr-based platforms. Remote Sensing of Environment 224, 92–103 (2019)
2019
-
[27]
Bara´, S. Characterizing the zenithal night sky brightness in large territories: how many samples per square kilometre are needed? Monthly Notices of the Royal Astronomical Society 473, 4164–4173 (2017)
2017
-
[28]
Gorelick, N. et al. Google earth engine: Planetary-scale geospatial analysis for everyone. Remote Sensing of Environment 202, 18–27 (2017)
2017
-
[29]
S., Ocan˜a, F
Zamorano, J., de Miguel, A. S., Ocan˜a, F. & Castan˜o, J. G. El uso de imagenes de sate´lite para combatir la contaminacio´n lum´ınica. Astronom´ıa 80–87 (2013)
2013
-
[30]
McHardy, T. et al. An improved method for retrieving nighttime aerosol optical thickness from the viirs day/night band (2015)
2015
-
[31]
Falchi, F. et al. The new world atlas of artificial night sky brightness. Science Advances 2 , e1600377 (2016)
2016
-
[32]
Roma´n, M. O. et al. Nasa’s black marble nighttime lights product suite. Remote Sensing of Environment 210, 113–143 (2018)
2018
-
[33]
& Falchi, F
Cinzano, P. & Falchi, F. The propagation of light pollution in the atmosphere. Monthly Notices of the Royal Astronomical Society 427, 3337–3357 (2012)
2012
-
[34]
Light-pollution model for cloudy and cloudless night skies with ground-based light sources
Kocifaj, M. Light-pollution model for cloudy and cloudless night skies with ground-based light sources. Applied optics 46, 3013–3022 (2007)
2007
-
[35]
& Royer, A
Aube´, M., O’Neill, N., Gigue`re, J.-D. & Royer, A. Using artificial sky glow to retrieve night time aerosol optical depth. In EGU General Assembly Conference Abstracts, vol. 11, 2591 (2009)
2009
-
[36]
& Solano Lamphar, H
Garc´ıa Gil, M., Francia Paya`s, P., San Mart´ı Pa´ramo, R. & Solano Lamphar, H. Contaminacio´n lum´ınica: Una visio´n desde el foco contaminante: el alumbrado artificial (Universitat Polite`cnica de Catalunya. Iniciativa Digital Polite`cnica, 2012)
2012
-
[37]
Unified model of radiance patterns under arbitrary sky conditions
Kocifaj, M. Unified model of radiance patterns under arbitrary sky conditions. Solar Energy 115, 40 –51 (2015)
2015
-
[38]
Night sky photometry with sky quality meter
Cinzano, P. Night sky photometry with sky quality meter. ISTIL Int. Rep 9 (2005)
2005
-
[39]
& Kyba, C
Schnitt, S., Ruhtz, T., Fischer, J., Ho¨ lker, F. & Kyba, C. Temperature stability of the sky quality meter. Sensors 13, 12166–12174 (2013)
2013
-
[40]
Ha¨nel, A. et al. Measuring night sky brightness: methods and challenges. Journal of Quantitative Spectroscopy and Radiative Transfer 205, 278–290 (2018)
2018
-
[41]
Mapa de brillo de fondo de cielo de la comunidad de madrid (2010)
Pila D´ıez, B. Mapa de brillo de fondo de cielo de la comunidad de madrid (2010). URL http://eprints.ucm.es/ 11364/
2010
-
[42]
& Pascual, S
Sa´nchez de Miguel, A., Zamorano, J., Go´mez Castan˜o, J. & Pascual, S. Evolution of the energy consumed by street lighting in spain estimated with dmsp-ols data. Journal of Quantitative Spectroscopy and Radiative Transfer 139, 109 – 117 (2014). URL http://www.sciencedirect.co...
2014
-
[43]
Zamorano, J. et al. Testing sky brightness models against radial dependency: A dense two dimensional survey around the city of madrid, spain. Journal of Quantitative Spectroscopy and Radiative Transfer 181, 52–66 (2016)
2016
-
[44]
Nightpod - nodding mechanism for the iss
Sabbatini, M. Nightpod - nodding mechanism for the iss. Tech. Rep. Experiment Record #9337, European Space Agency, Noordwijk, The Netherlands (2014). Acknowledgements This work was supported by the EMISSI@N project (NERC grant NE/P01156X/1), COST (European Cooperation in Scien...
2014
Reviewed August 14, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.