AUTHOR=Small Christopher TITLE=Spatiotemporal Characterization of VIIRS Night Light JOURNAL=Frontiers in Remote Sensing VOLUME=Volume 2 - 2021 YEAR=2021 URL=https://www.frontiersin.org/journals/remote-sensing/articles/10.3389/frsen.2021.775399 DOI=10.3389/frsen.2021.775399 ISSN=2673-6187 ABSTRACT=The Visible Infrared Imaging Radiometer Suite (VIIRS) Day Night Band (DNB) on Suomi NPP satellite provides almost a decade of daily observations of night light. The most common approach to night light change detection involves direct attribution of observed changes to the phenomenon of interest. Implicit is the assumption that other forms of actual and apparent change in the light source are negligible or non-existent. An alternative approach is to characterize the spatiotemporal variability prior to deductive attribution of causation so the attribution can be made in the context of the full range of spatial and temporal variation. The primary objective of this study is to characterize night light variability over a range of spatial and temporal scales to provide a context for interpretation of night light changes observed on both subannual and interannual time scales. This analysis is based on temporal moments, spatial correlation and Empirical Orthogonal Function (EOF) analysis. A key result of this study is the pervasive heteroskedasticity of VIIRS monthly mean night light. Specifically, the monotonic decrease of variability with increasing mean brightness. Anthropogenic night light is remarkably stable on subannual time scales while background luminance varies considerably. The overall variance partition derived from the eigenvalues of the spatiotemporal covariance matrix are 88%, 2% and 2% for spatial, seasonal and interannual variance (respectively) in the most diverse geographic region on Earth (Eurasia). Heteroskedasticity is a pervasive characteristic of the monthly composites and is present in all areas for all months of the year, suggesting that much, if not most, of the month-to-month variability may be related to luminance of otherwise stable sources subjected to multiple aspects of the imaging process varying in time. Given the skewed distribution of all night light arising from radial peripheral dimming of bright sources subject to atmospheric scattering, even aggregate metrics using thresholds must be interpreted in light of the fact that much larger numbers of more variable low luminance pixels may statistically overwhelm smaller numbers of stable higher luminance pixels and cause apparent changes related to the imaging process to be interpreted as actual changes in the light sources.