A demanding challenge in atmospheric research is the night-time characterization of aerosols using passive techniques, that is, by extracting information from scattered light that has not been emitted by the observer. Satellite observations of artificial night-time lights have been used to retrieve some basic integral parameters, like the aerosol optical depth. However, a thorough analysis of the scattering processes allows one to obtain substantially more detailed information on aerosol properties. In this Letter we demonstrate a practicable approach for determining the aerosol particle size number distribution function in the air column, based on the measurement of the angular radiance distribution of the scattered light emitted by night-time lights of cities and towns, recorded from low Earth orbit. The method is self-calibrating and does not require the knowledge of the absolute city emissions. The input radiance data are readily available from several spaceborne platforms, like the VIIRS-DNB radiometer onboard the Suomi-NPP satellite.

The visual brightness of the night sky is not a single-valued function of its brightness in other photometric bands, because the transformations between photometric systems depend on the spectral power distribution of the skyglow. We analyze the transformation between the night sky brightness in the Johnson-Cousins V band (mV, measured in magnitudes per square arcsecond, mpsas) and its visual luminance (L, in SI units cd m^-2) for observers with photopic and scotopic adaptation, in terms of the spectral power distribution of the incident light. We calculate the zero-point luminances for a set of skyglow spectra recorded at different places in the world, including strongly light-polluted locations and sites with nearly pristine natural dark skies. The photopic skyglow luminance corresponding to mV=0.00 mpsas is found to vary between 1.11-1.34 x 10^5 cd m^-2 if mV is reported in the absolute (AB) magnitude scale, and between 1.18-1.43 x 10^5 cd m^-2 if a Vega scale for mV is used instead. The photopic luminance for mV=22.0 mpsas is correspondingly comprised between 176 and 213 microcd m^-2 (AB), or 187 and 227 microcd m^-2 (Vega). These constants tend to decrease for increasing correlated color temperatures (CCT). The photopic zero-point luminances are generally higher than the ones expected for blackbody radiation of comparable CCT. The scotopic-to-photopic luminance ratio (S/P) for our spectral dataset varies from 0.8 to 2.5. Under scotopic adaptation the dependence of the zero-point luminances with the CCT, and their values relative to blackbody radiation, are reversed with respect to photopic ones.

Determining the all-sky radiance distribution produced by artificial light sources is a computationally demanding task that generally requires an intensive calculation load. We develop in this work an analytic formulation that provides the all-sky radiance distribution produced by an artificial light source as an explicit and analytic function of the observation direction, depending on two single parameters that characterize the overall effects of the atmosphere. One of these parameters is related to the effective attenuation of the light beams, whereas the other accounts for the overall asymmetry of the combined scattering processes in molecules and aerosols. By means of this formulation a wide range of all-sky radiance distributions can be efficiently and accurately calculated in a short time. This substantial reduction in the number of required parameters, in comparison with other currently used approaches, is expected to facilitate the development of new applications in the field of light pollution research.

Light pollution poses a growing threat to optical astronomy, in addition to its detrimental impacts on the natural environment, the intangible heritage of humankind related to the contemplation of the starry sky and, potentially, on human health. The computation of maps showing the spatial distribution of several light pollution related functions (e.g. the anthropogenic zenithal night sky brightness, or the average brightness of the celestial hemisphere) is a key tool for light pollution monitoring and control, providing the scientific rationale for the adoption of informed decisions on public lighting and astronomical site preservation. The calculation of such maps from satellite radiance data for wide regions of the planet with sub-kilometric spatial resolution often implies a huge amount of basic pixel operations, requiring in many cases extremely large computation times. In this paper we show that, using adequate geographical projections, a wide set of light pollution map calculations can be reframed in terms of two-dimensional convolutions that can be easily evaluated using conventional fast Fourier-transform (FFT) algorithms, with typical computation times smaller than 10^-6 s per output pixel.

The light pollution levels experienced at any given site generally depend on a wide number of artificial light sources distributed throughout the surrounding territory. Since photons can travel long distances before being scattered by the atmosphere, any effective proposal for reducing local light pollution levels needs an accurate assessment of the relative weight of all intervening light sources, including those located tens or even hundreds of km away. In this paper we describe several ways of quantifying and visualizing these relative weights. Particular emphasis is made on the aggregate contribution of the municipalities, which are -- in many regions of the world -- the administrative bodies primarily responsible for the planning and maintenance of public outdoor lighting systems.

Several light pollution indicators are commonly used to monitor the effects of the transition from outdoor lighting systems based on traditional gas-discharge lamps to solid-state light sources. In this work we analyze a subset of these indicators, including the artificial zenithal night sky brightness in the visual photopic and scotopic bands, the brightness in the specific photometric band of the widely used Sky Quality Meter (SQM), and the top-of-atmosphere radiance detected by the VIIRS-DNB radiometer onboard the satellite Suomi-NPP. Using a single-scattering approximation in a layered atmosphere we quantitatively show that, depending on the transition scenarios, these indicators may show different, even opposite behaviors. This is mainly due to the combined effects of the changes in the sources' spectra and angular radiation patterns, the wavelength-dependent atmospheric propagation processes and the differences in the detector spectral sensitivity bands. It is suggested that the possible presence of this differential behavior should be taken into account when evaluating light pollution indicator datasets for assessing the outcomes of public policy decisions regarding the upgrading of outdoor lighting systems.

Light pollution is a worldwide phenomenon whose consequences for the natural environment and the human health are being intensively studied nowadays. Most published studies address issues related to light pollution inland. Coastal waters, however, are spaces of high environmental interest, due to their biodiversity richness and their economical significance. The elevated population density in coastal regions is accompanied by correspondingly large emissions of artificial light at night, whose role as an environmental stressor is increasingly being recognized. (...) At the same time, the marine surface environment provides a stage free from obstacles for measuring the dependence of the skyglow on the distance to the light polluting sources, and validating (or rejecting) atmospheric light propagation models. In this work we present a proof-of-concept of a gimbal measurement system that can be used for zenithal skyglow measurements on board both small boats and large vessels under actual navigation conditions. We report the results obtained in the summer of 2016 along two measurement routes in the Mediterranean waters offshore Barcelona, travelling 9 and 31.7 km away from the coast. The atmospheric conditions in both routes were different from the ones assumed for the calculation of recently published models of the anthropogenic sky brightness. They were closer in the first route, whose results approach better the theoretical predictions. The results obtained in the second route, conducted under a clearer atmosphere, showed systematic differences that can be traced back to two expected phenomena, which are a consequence of the smaller aerosol content: the reduction of the anthropogenic sky glow at short distances from the sources, and the slower decay rate of brightness with distance, which gives rise to a relative excess of brightness at large distances from the coastline.

The formal link between magnitudes per square arcsecond and luminance is discussed in this paper. Directly related to the human visual system, luminance is defined in terms of the spectral radiance of the source, weighted by the CIE V($\lambda$) luminous efficiency function, and scaled by the 683 lm/W luminous efficacy constant. In consequence, any exact and spectrum-independent relationship between luminance and magnitudes per square arcsecond requires that the latter be measured precisely in the CIE V($\lambda$) band. The luminance value corresponding to m_VC=0 (zero-point of the CIE V($\lambda$) magnitude scale) depends on the reference source chosen for the definition of the magnitude system. Using absolute AB magnitudes, the zero-point luminance of the CIE V($\lambda$) photometric band is 10.96 x 10^4 cd/m^2.

A recurring question arises when trying to characterize, by means of measurements or theoretical calculations, the zenithal night sky brightness throughout a large territory: how many samples per square kilometer are needed? The optimum sampling distance should allow reconstructing, with sufficient accuracy, the continuous zenithal brightness map across the whole region, whilst at the same time avoiding unnecessary and redundant oversampling. This paper attempts to provide some tentative answers to this issue, using two complementary tools: the luminance structure function and the Nyquist-Shannon spatial sampling theorem. The anaysis of several regions of the world, based on the data from the New world atlas of artificial night sky brightness (Falchi et al 2016, Sci. Adv. 2, e1600377) suggests that, as a rule of thumb, about one measurement per square kilometer could be sufficient for determining the zenithal night sky brightness of artificial origin at any point in a region to within 0.1 magV/arcsec2 (in the root-mean-square sense) of its true value in the Johnson-Cousins V band. The exact reconstruction of the zenithal night sky brightness maps from samples taken at the Nyquist rate seems to be considerably more demanding.

Under stable atmospheric conditions, the zenithal brightness of the urban sky varies throughout the night following the time course of the anthropogenic emissions of light. Different types of artificial light sources (e.g. streetlights, residential, and vehicle lights) present specific time signatures, and this feature makes it possible to estimate the amount of sky brightness contributed by each one of them. Our approach is based on transforming the time representation of the zenithal sky brightness into a modal coefficients one, in terms of the time course signatures of the sources. The modal coefficients, and hence the absolute and relative contributions of each type of source, can be estimated from the measured brightness by means of linear least squares fits. A method for determining the time signatures is described, based on wide-field time-lapse photometry of the urban nightscape. Our preliminary results suggest that artificial light leaking out of the windows of residential buildings may account for a significant share of the time-varying part of the zenithal sky brightness, whilst the contribution of the vehicle lights seems to be significantly smaller.