First Estimation of Global Trends in Nocturnal Power Emissions Reveals Acceleration of Light Pollution

: The global spread of artificial light is eroding the natural night-time environment. The estimation of the pattern and rate of growth of light pollution on multi-decadal scales has nonetheless proven challenging. Here we show that the power of global satellite observable light emissions increased from 1992 to 2017 by at least 49%. We estimate the hidden impact of the transition to solid-state light-emitting diode (LED) technology, which increases emissions at visible wavelengths undetectable to existing satellite sensors, suggesting that the true increase in radiance in the visible spectrum may be as high as globally 270% and 400% on specific regions. These dynamics vary by region, but there is limited evidence that advances in lighting technology have led to decreased emissions.


Introduction
Satellite images of the Earth at night strikingly illustrate the extent of the distribution of artificial light sources associated with human settlement, industry and transport links, across its surface [1]. The quantity of artificial light detected from satellites includes a mixed signal of direct emissions from outdoor lighting that is not vertically shielded, as well as reflected light from the ground, buildings and vegetation, and light scattered within the atmosphere. This signal has increasingly been used as an important indicator of urbanisation, industrial activity [2,3] and economic development [4,5]. It has also been used as a measure of light pollution [6,7], which has long been a concern of astronomers [8,9] but has more recently been proposed as an emerging major threat both to public health and natural ecosystems. Indeed, light pollution has been implicated in increased human cancer risk [10,11], population declines in wild species [12,13], the restructuring of ecological communities [14,15], and the disruption of key ecosystem services such as pollination [16,17].
To determine how the global magnitude of the impacts of light pollution is changing on multi-decadal scales, as has been performed for many other anthropogenic environmental pressures, time series data on artificial light emissions are required. This proved challenging for a combination of reasons, including technical difficulties in distinguishing faint artificial light sources from natural ones, limitations of the available satellite sensors, and a paucity of ground-truthed sources. In consequence, studies to date [18,19] have been severely limited in their spatial and temporal coverage, and in the inferences that could be drawn. Here we provide a global estimation of changes in artificial night-time lighting over a 25-year period, using ground-truth data and a conservative estimation of the potential errors associated particularly with a lack of information on emissions in the blue part of the spectrum.
The Defense Meteorological Satellite Program Operational Line scan System (DMSP-OLS) forms the basis for the longest publicly available time series of artificial light at night from space, with uncalibrated composite images freely obtainable from a digital archive since 1992 [20]. Since December 2011, the Visible Infrared Imaging Radiation Suite (VIIRS) Day/Night Band (DNB) on board the NASA/NOAA Suomi National Polarorbiting Partnership (Suomi-NPP) satellite has provided a radiance calibrated product with global coverage. However, quantifying changes in artificial light emissions at night at the global scale from satellite imagery was hindered by several factors.
First, the lack of on-board calibration of the DMSP-OLS images means that irradiance values from different images from the satellite cannot be directly compared without complex cross-calibration. The sensor performance degrades with time and varies between images, years and satellites, calibrations are adjusted "on the fly" and, due to the low memory range and low dynamic range of the sensor, archived radiance values are represented in images as an uncalibrated digital number (DN) between 1 and 63, which cannot be compared between different composite images in a straightforward manner.
Second, in brightly lit urban and industrial areas the signal in DMSP-OLS data frequently reaches saturation at a DN of 63, making change detection within such areas impossible [21].
Third, the resolution of DMSP-OLS images is fairly low (2.7 km), the accuracy of georeferencing is often poor and images suffer from "blooming" effects [22], resulting in a blurred image around bright light sources, making the identification of fine scale features problematic [23].
Fourth, the single-channel spectral response of the DMSP-OLS sensor is close to zero at wavelengths of visible light below 500 nm. White phosphor-based LED lighting typically emits at a peak radiance of around 465 nm, and so as outdoor lighting increasingly switches from gas-discharge technologies such as high-pressure sodium to solid-state LED lighting, an undetermined, but increasing, proportion of emitted light is undetected ( Figure 1). This is particularly significant, because LED lighting has been widely championed as a technological revolution with increased energy efficiency compared with other white light sources, that enables the potential of outdoor lighting at reduced carbon dioxide emissions and financial costs, but recent analysis suggests that these reduced costs may encourage increased overall light emissions [19].
Remote Sens. 2021, 13, x FOR PEER REVIEW 3 of 11 [29] in which data are collected at three different fixed gain levels, with a preflight sensor calibration, at the new moon. This method has been used to produce radiance calibrated composite images for the years 1996, 1999, 2000, 2002, 2004, 2005 and 2010.
(a) (b) Here, we produced a global time series of satellite-detectable emitted power of outdoor lighting between 500 nm and 900 nm from 1992 to 2017 by integrating the uncalibrated and radiance-calibrated DMSP-OLS time series, VIIRS imagery and a data set of regional energy consumption by municipal street lighting [30], which is globally the principal source of outdoor artificial light at night. Uncalibrated DMSP-OLS data were fitted Some, but not all, of the issues associated with tracking global artificial night-time light emissions were resolved with the availability of the VIIRS-DNB imagery from 2012 onwards. VIIRS is radiance calibrated and has a higher dynamic range without saturation, allowing comparison between images and detection of change even within brightly lit areas [23]. It has a higher spatial resolution (740 m), and the archived cloud-free composite products suffer less from errors in georeferencing and "blooming" effects. However, the spectral response of the sensor is similarly low below wavelengths of 500 nm (Figure 1), leaving a proportion of light emitted from LEDs undetected. Recently, images taken from the International Space Station (ISS) [24] and satellites such as JL-13B [25] have provided higher resolution night-time imagery with three-band spectral resolution (red, green, and blue visual bands), potentially allowing changes in lighting technology to be detected directly [24]. However, at present the DMSP-OLS data remain the only source for assessing decadal-scale changes in radiance.
To overcome some of the limitations, methods have been proposed for inter-calibrating archived DMSP-OLS imagery for change detection [2,26,27], and the data set in uncalibrated and cross-calibrated forms has been used widely to investigate economic development and ecological impacts. Since 2016, attempts have been made to provide calibrated ground truthing for night-time light measurements [28]. A radiance calibrated version of the DMSP-OLS data, which does not suffer from saturation effects, was developed [29] in which data are collected at three different fixed gain levels, with a preflight sensor calibration, at the new moon. This method has been used to produce radiance calibrated composite images for the years 1996, 1999, 2000, 2002, 2004, 2005 and 2010. Here, we produced a global time series of satellite-detectable emitted power of outdoor lighting between 500 nm and 900 nm from 1992 to 2017 by integrating the uncalibrated and radiance-calibrated DMSP-OLS time series, VIIRS imagery and a data set of regional energy consumption by municipal street lighting [30], which is globally the principal source of outdoor artificial light at night. Uncalibrated DMSP-OLS data were fitted to calibrated DMSP-OLS data sets using non-linear regression methods, then all three satellite data sets were calibrated using linear methods to fit to the energy consumption data. Using assumptions about the rate of conversion to solid-state LED lighting, we then estimated the potential range of this "hidden" component of global light at night.

Background
To produce a consistent time series of radiance for each country in the world, we undertook four steps in data processing. First, we used a robust regression approach to calibrate the uncalibrated DMSP-OLS images to the nearest calibrated image in time. Second, we assimilated provincial data on energy consumption from street lighting (as a proxy for expected light emissions) to smooth and refine the calibration of all images to produce a consistent time series. Third, we utilised the time difference in image acquisition between the DMSP-OLS and VIIRS satellite images to adjust this time series to produce two separate time series of radiance for each nation state, for early evening and midnight respectively. Fourth, we estimated the effect that changes in light technology towards blue-rich solid-state (LED) lighting may had on detected light emissions since 2012, given the poor detection of blue light by the satellites, and produce high and low estimates of the actual radiance given the detected radiance under three different scenarios of technological change. A summary of the procedure is shown in Figure S1.

Calibration of Uncalibrated DMSP-OLS Data
The annual, uncalibrated DMSP-OLS data have no absolute or relative intercalibration, and the response of each image to absolute radiance is non-linear. Previous studies have attempted to intercalibrate these images by assuming that certain regions, such as Sicily [31], are essentially invariant in radiance over time and have used conventional regression techniques to intercalibrate global images using such "invariant regions" [32]. However, it is likely that few locations have truly stable lights over several years [33], for example, the island of Sicily increased its energy consumption by street lighting by 20% over 20 years [30]. Lamp "lumen depreciation," a decay in brightness over time in most bulb technologies as bulbs age, may also limit the effectiveness of intercalibrations, even if the number of light sources and their bulb type remains constant. Other approaches used "robust" regression techniques, such as quantile regression through the median [26] or ridgeline regression [34], which do not assume that entire regions are invariant, but rather that changes within the calibration region are incremental and localised. Such "robust regression" approaches are less sensitive to localised changes than least-squared regression but assume that the brightness of most pixels remains approximately constant between two images, even if localised pixels increase or decrease substantially. More recently, the TERRA-VEGA project [28] attempted to ground-truth artificial light using a mobile source. However, such ground -based calibration sources have only been available since 2016. Recently, calibration of the VIIRS data using stars was performed, but as far as we know not yet implemented for the monthly data composites [35].
In order to calibrate the uncalibrated DMSP-OLS data to radiance values, we used a robust regression method on a calibration area encompassing Italy (bounded between 4 and 19 • E and 35.5 and 47 • N). We made no assumption that this area has remained unchanged in radiance over time, but rather that most changes in anthropogenic light are incremental and localised in space. This region was selected as it has historical data of energy use in municipal street lighting at the provincial level for use in subsequent calibration steps. We first matched each uncalibrated image with the closest (in time) radiance-calibrated image. Calibrated images exist for 1996,1999,2000,2002,2004,2005,2010 and 2011. Each value in the uncalibrated DMSP images is a digital number (DN) between 0 and 63, while the radiance-calibrated VIIRS and DMSP images are calibrated to nW cm −2 sr −1 [2]. For each digital number (DN) value of the uncalibrated data within the calibration region, we fitted a kernel density estimation using the density function in the statistical software, R [36]. We then took the maximum value of the kernel density (the mode of the fitted distribution) as the most likely absolute radiance value represented by each digital number. This step was based on the assumption that, while many absolute radiance values may have changed between the acquisition of each uncalibrated image and the nearest calibrated image, the modal value is robust to changes in the extremes of the distribution. To these modal values we fitted a 10th order polynomial regression, to produce a smoothed conversion from DN values to calibrated radiances. The uncalibrated DMSP images were intercalibrated in this way using the closest calibrated DMSP images following the uncalibrated image; for example uncalibrated images from 1992 to 1996 were calibrated using the 1996 calibrated image. Uncalibrated images from 2011 to 2014 were intercalibrated against the last available calibrated image in 2010. For many years there is more than one uncalibrated image from different satellites. Since the relationship is fitted to a single modal radiance value for each digital number, this method is robust against outliers, scatter, and overfitting, despite the large number of parameters in the polynomial fit.
While these steps help to linearise the relationship between the intercalibrated images and absolute radiance values across most of the dynamic range of the instruments, the intercalibrated images retain the problem of saturation at high radiance levels. This limits the ability of the data set to detect change at high radiance levels of, say, sports stadia or urban centres. As different images saturate at different light levels, we capped all data at 100 nW cm −2 sr −1 , the lowest saturation point for the images.

Smoothing and Recalibration of All Data Series Using Energy Consumption Data
Despite pixel-by-pixel intercalibration, it is clear that some annual images tend to over or underestimate radiant light when sampled across large areas, as evidenced by interannual variation around the temporal trend due to the saturation effects in the uncalibrated DMSP, as well as differences in the sensitivity of different images to detect low levels of light in rural and sparsely populated areas. To ensure that all three data sets (raw DMSP-OLS, radiance calibrated DMSP-OLS and VIIRS) are comparable, we used energy consumption data for municipal street lighting for 75 Italian regions (publicly available from [30]) to smooth the calibrated data. This process incorporates four steps: (i) Using the 2012 VIIRS image, we fitted a linear least-squares regression between the total light emitted (per province, as the sum of values from all pixels within that province) and the energy consumption of municipal street lighting in 2012. The slope of this regression yields an estimated conversion factor between provincial energy consumption and the expected radiance produced by the lights. We then converted the radiance units of nW cm −2 sr −1 into total radiant energy within the range of wavelengths detected by the satellite (in W), assuming an isotropic distribution of radiance (equal radiance in all directions). This conversion was undertaken for all satellite images. (ii) We then plotted the relationship between annual total light emitted (converted to detectable radiant energy in W) and estimated detectable radiant energy from the energy consumption data for each province (as a proxy of the real change slope). We calculated the mean residual from this relationship for each province. In this case a positive residual value for a province indicates a higher-than-expected energy consumption value for a given observed radiance; this may be due to, for example, effective shielding of street lighting or energy inefficient lighting within the province. Conversely a negative residual value may indicate emissions from private or commercial sources, low shielding of lighting or more efficient street lighting. (iii) We then fitted separate linear regressions for each image, between the total light emitted (converted to radiant energy in W) and estimated detectable radiant energy from the energy consumption data, adjusted by subtracting the province-specific residual values. The slope for each line was recalculated at this stage. (iv) We then subtracted the mean residual value for each province across images (assumed to represent consistent, localised differences in lighting efficiency and type) from the provincial data and recalculated the regression slopes for each image. The coefficients of this regression were then used to recalibrate the values to convert each image to radiant energy in W. Steps ii-iv were repeated iteratively until radiance values converged on stable values, and the converged values were used to calibrate global images.
Following the calibration of global images, the total power emitted per nation state and continent was extracted to reconstruct national, continental and global long-term trends in light emissions. To validate the approach, we compared the annual power predicted for Brazil, Croatia, Czech Republic, Denmark, France, Germany, Greece, Ireland, Norway, Portugal and Switzerland, to that calculated from national-level light consumption data using the above method ( Figure S2) and provinces with reliable data; Italy, Spain, Portugal and Greece (Figures S3 and S4). For the global extraction of the data sets, the Google Earth Engine was used [37] and the rest of the processing was performed using standard Python code using the libraries Numpy [38], Astropy [39], Pandas [40], Scipy [41], Matplotlib [42] and Selenium for the web scraping [43].

Adjusting National Data for Differing Times of Acquisition
DMSP-OLS and VIIRS images of the calibration region show good correspondence for the years in which data from both sensors are available ( Figure S5). However, VIIRS images are acquired between 00:00 and 02:00 solar time (hereafter "post-midnight"), while DMSP-OLS images are typically acquired around 21:30 solar time (hereafter "evening") [1]. Data obtained from each satellite may therefore not be comparable, as human activity (and hence light emissions) tends to vary throughout the night; furthermore, some nations, regions and cities switch off or dim sources such as street-lights later at night. Data from both satellites are available for 2012, so we used this year to adjust the calibrations for each country. To produce a post-midnight time series, we recalibrated the DMSP-OLS derived data to match the VIIRS data for 2012; to produce an evening time series we recalibrated the VIIRS data to match the DMSP-OLS data for 2012.

Adjusting Data for Increased Emission of Blue Light from LEDs
To estimate the effect of conversion from gas-discharge street lighting (such as high-or low-pressure sodium lamps) to solid-state lighting (LED) on our reconstructed time series, we adjusted the data post-2007 (prior to which LED lighting was not widely adopted) to estimate the radiance as perceived by the human eye, as opposed to detected by the satellite sensor. To estimate possible rates of conversion we used a post-2007 rate of replacement of gas-discharge street lighting equivalent to that recorded by the Spanish street lighting manufacturers association, Asociacion Espanola de Fabricantes de Iluminacion (ANFALUM) [44]. As a low range estimate, we assume that sodium lighting is replaced with LEDs at a correlated colour temperature (CCT) of 2000 K, with low levels of blue light, similar detectability to sodium lighting and broadly equivalent to continued use of sodium discharge lighting. As a mid-range estimate, replacement is assumed to be with LEDs with a CCT of 3000 K, the maximum recommended by the International Dark Sky Association (IDA) until 2020 and other organizations such as the European Joint Research Centre [45] or the Government of France [46]. As a high-end estimate we assume that gas-discharge lighting is progressively replaced by blue-rich LED lighting at a CCT of 4000 K. In all scenarios we assume that sodium lights are replaced after 24,000 h of continuous operation and are used for 4100 h per year. In order to estimate the emissions in the blue part of the spectrum, we used the Nikon Blue Band [24]. This is currently the only observation band available on a large scale through images acquired from the International Space Station [24]. To convert from the VIIRS-DMSP power to the Nikon Blue band, the 2000 K flux needs to be multiplied by 0.1, 3000 K by 0.25, and 4000 K by 0.5. These conversion factors were calculated using the procedures described in [24].

Details of the Method, Limitations, and Validation
This method uses the energy consumption of the street lighting of provinces of Italy as a proxy for light emissions. We do not assume that all the lights visible are streetlights, and by utilising residuals from the linear relationship between energy consumption and detected radiance we explicitly allow for a component of detected light that does not come from street lighting (step (ii) in Section 2.3 above). Previous studies have shown that trends in street lighting energy consumption are often reasonable proxies of trends in detected radiance [2,33]. Nevertheless, the residuals fitted for each province allow for local variation in the magnitude of other lighting (for example domestic and industrial sources). The street lighting energy consumption was used as a reference for the trend in growth in lighting (the slope of the regression) but this trend was offset by a fitted residual (the intercept) that can be considered an estimate of a fixed proportion of lighting from other sources. Figures S6-S8 and on the data availability statement.

Additional validation and intermediate products can be found on
The proportion of night-time lighting that derives from streetlights is unclear and debated in the literature; it is likely to vary regionally between countries and cities ( [7] compile different values that range from 12% in Flagstaff, AZ, USA to nearly 100% in Spain).

Results and Discussion
Globally, our analysis shows that satellite detectable light has increased at an approximately constant rate from 1992 to 2017, and by 49% over the 25 years ( Figure 2). Utilising data of typical conversion and replacement rates to estimate the possible fraction of undetected light assuming different scenarios of uptake and replacement of gas-discharge lighting with LEDs, suggests however that there has been a marked recent acceleration to global emitted power (Figure 2). The projected blue light global total emission of 513 MW in 2017 would be as high as 806 MW if globally 40% of lighting were LEDs with a colour temperature of 3000 K, and 1294 MW at a colour temperature of 4000 K. Thus, not only did the "lighting revolution" associated with the transition to solid-state technology fail to result in the intended reduction in energy consumption and instead led to a rebound effect of increased use [29], but it actually resulted in a dramatic rise in that consumption compared with that which had previously been occurring. Such a recent global acceleration has also been seen in other emerging drivers of environmental change such as plastic waste production [47], and synthetic chemicals, including pesticides [48]. This raises the prospect that this change becomes increasingly challenging to redress, especially if previously unrealised synergies between the impacts of drivers emerge.
Remote Sens. 2021, 13, x FOR PEER REVIEW 7 of 11 in 2017 would be as high as 806 MW if globally 40% of lighting were LEDs with a colour temperature of 3000 K, and 1294 MW at a colour temperature of 4000 K. Thus, not only did the "lighting revolution" associated with the transition to solid-state technology fail to result in the intended reduction in energy consumption and instead led to a rebound effect of increased use [29], but it actually resulted in a dramatic rise in that consumption compared with that which had previously been occurring. Such a recent global acceleration has also been seen in other emerging drivers of environmental change such as plastic waste production [47], and synthetic chemicals, including pesticides [48]. This raises the prospect that this change becomes increasingly challenging to redress, especially if previously unrealised synergies between the impacts of drivers emerge. At the continental scale, persistent rates of growth in detected light are observed in Asia, South America, Oceania and Africa; while in Europe, detected light increased until around 2010 and the levelled off, and in North America it appears to be in decline ( Figure  3). This is consistent with other global analyses of lighting, showing rapid urbanisation and economic growth in developing countries, and some evidence of a slowing down or decrease in light intensity in some developed nations [49]. Some proportion of this apparent lack of growth may be attributed to a saturation effect, as few remaining areas are unlit; there may also be observable effects of the decline of extractive industries [26], and better lighting design, resulting in less light being emitted above the horizontal plane. However, any conclusions as to the apparent levelling off or decrease in lighting intensity At the continental scale, persistent rates of growth in detected light are observed in Asia, South America, Oceania and Africa; while in Europe, detected light increased until around 2010 and the levelled off, and in North America it appears to be in decline (Figure 3). This is consistent with other global analyses of lighting, showing rapid urbanisation and economic growth in developing countries, and some evidence of a slowing down or decrease in light intensity in some developed nations [49]. Some proportion of this apparent lack of growth may be attributed to a saturation effect, as few remaining areas are unlit; there may also be observable effects of the decline of extractive industries [26], and better lighting design, resulting in less light being emitted above the horizontal plane. However, any conclusions as to the apparent levelling off or decrease in lighting intensity in Europe and North America should be treated with caution as the shifting spectral distribution of light towards blue-rich LEDs is likely to mask the effects of increases [49].
Remote Sens. 2021, 13, x FOR PEER REVIEW in Europe and North America should be treated with caution as the shifting sp tribution of light towards blue-rich LEDs is likely to mask the effects of increase

Figur
Global emitted power (Nikon blue band [14]) detected by satellites from artificial light sources from 1992 to 2017. squares represent annual DMSP-OLS composite data, filled squares represent radiance calibrated DMSP-OLS data, circles represent VIIRS Day/Night band data. Plotted data points assume constant spectral composition of light emis The shaded areas represent the possible range of undetected light assuming a recent phased transition from high pre sodium lighting to LEDs of color temperature 3000 K (dark grey) or 4000 K (light grey).
At the country scale, marked patterns of a decrease and increase in detecte power of lighting are often related to countries in conflict or recovering from co ing this period ( Figure 2); Syria [1], Ukraine [50], Central African Republic, an [51] (see Figure S9) all show marked decreases in emitted power of lighting ov riod, while Bosnia and Herzegovina, Afghanistan, Angola, Mozambique and S show decreases. Declines in Nigeria and Uzbekistan are apparently due, in larg changes in light emissions in areas of oil extraction.
The recent growth, and relative uncertainty in the scale of global emissio wavelength blue light is significant as lower wavelengths of visible light are imp the disruption of circadian rhythms in humans and other animals and have be to be more effective than longer wavelengths at driving a range of ecological r from phototaxis in insects to disrupting reproduction in mammals [49]. This h the need for multispectral nighttime imagery to monitor the ongoing changes in tral composition of artificial light in the environment as technology changes, as w creative use of multiple data sources to reconstruct the changing spectral signatu turnal lightscapes in the past.  [14]) detected by satellites from artificial light sources from 1992 to 2017. Open squares represent annual DMSP-OLS composite data, filled squares represent radiance calibrated DMSP-OLS data, filled circles represent VIIRS Day/Night band data. Plotted data points assume constant spectral composition of light emissions. The shaded areas represent the possible range of undetected light assuming a recent phased transition from high pressure sodium lighting to LEDs of color temperature 3000 K (dark grey) or 4000 K (light grey).
At the country scale, marked patterns of a decrease and increase in detected emitted power of lighting are often related to countries in conflict or recovering from conflict during this period ( Figure 2); Syria [1], Ukraine [50], Central African Republic, and Yemen [51] (see Figure S9) all show marked decreases in emitted power of lighting over the period, while Bosnia and Herzegovina, Afghanistan, Angola, Mozambique and Somalia all show decreases. Declines in Nigeria and Uzbekistan are apparently due, in large part, to changes in light emissions in areas of oil extraction.
The recent growth, and relative uncertainty in the scale of global emissions of low wavelength blue light is significant as lower wavelengths of visible light are implicated in the disruption of circadian rhythms in humans and other animals and have been shown to be more effective than longer wavelengths at driving a range of ecological responses, from phototaxis in insects to disrupting reproduction in mammals [49]. This highlights the need for multispectral nighttime imagery to monitor the ongoing changes in the spectral composition of artificial light in the environment as technology changes, as well as the creative use of multiple data sources to reconstruct the changing spectral signature of nocturnal lightscapes in the past.

Conclusions
Despite decades of remote sensing of night-time lights, reconstructing recent global trends in light pollution remains challenging due to the lack of calibration and saturation of the early DMSP-OLS datasets, the lack of spectral information from panchromatic sensors, and the relatively poor sensitivity of existing satellite sensors to blue light. This analysis places bounds on the uncertainty in global trends in nocturnal power emissions by combining remote sensing data from different sources with power consumption data from municipal street lighting. While considerable uncertainty remains, it is clear that light emissions have increased by at least 49% over the 25-year period studied, and that this estimate is likely to be conservative; the transition to blue-rich solid-state LED lighting likely masks a much larger increase that is under-detected by existing panchromatic satellite data.
Supplementary Materials: The following are available online at https://www.mdpi.com/article/10 .3390/rs13163311/s1. Figure S1: Flow diagram explaining the procedure of intercalibration, Figure  S2: Predicted vs. measured power of street lighting for countries and years within the study period for which it was available: Brazil, Croatia, Czech Republic, Denmark, France, Germany, Greece, Ireland, Norway, Portugal and Switzerland. Measured power obtained from publicly available data from national and regional governments in each case, Figure S3: Example of test of the trends for Portugal (a), Italy (b), Greece (c) and Spain (d). In blue are the DMSP uncalibrated data, in green the DMSP calibrated and in red the VIIRS data. In black, the estimated power expected from the street lighting statistics. These plots are using the midnight reference, aka. VIIRS. In the case of Spain, [33] showed how the national values were not reliable but that some of the provinces were reliable. Data for Greece were like those for Portugal but on a smaller scale. We found that the statistics have a delay of 2-3 years because they correspond not to the actual year when the energy is used but when the payment is effectively performed (Private communication Victor Kouloumpis for Greece and Raul Lima for Portugal). Note that the economic crisis was visible in the night lights trend of Greece way sooner than for Spain, Italy or Portugal. In Greece the change in slope happened in 2006, while in Portugal, Spain and Italy it happened in 2008-9, All before LEDs were introduced, Figure  S4: Example of test of the trend of regions of Faro (Portugal) (a), Grosseto (Italy) (b), Thessaloniki (Greece) (c) and Cantabria (Spain) (d). In blue are the DMSP uncalibrated data, in green the DMSP calibrated and in red the VIIRS data. In black, the estimated power expected from the street lighting statistics. The four provinces have very different behaviour, but smooth trends that match the slope of the corresponding estimation from street lighting energy consumption. The regions are in some cases provinces in others NUTS, depending on the availability of the street lighting statistics and shape files. Spain, Italy and Greece correspond to NUTS3, and Portugal to provinces, Figure S5: Relationship between DMSP-OLS 2010-2011 and SNPP/VIIRS/DNB 2012 for Italian provinces, Figure S6: Ratio of DMSP-OLS, VIIRS-DNB, versus G/R ratio from the ISS. Calculated following the same procedure as [24]. Same legend as plots in [24]. For any given spectra, the difference is lower than 15%, with the exception of black body lamps, that can reach 20% maximum, Figure S7: Examples of not rescaled raw DMSP, calibrated DMSP and VIIRS for different locations in Italy. In blue are the DMSP uncalibrated data, in green the DMSP calibrated and in red the VIIRS data. In black, the estimated power expected from the street lighting statistics. The gap between calibrated DMSP 2010 and VIIRS 2012 can provide an estimation of how much of the light is emitted in the first part of the night and how much in the second part. As can be seen, for three of the provinces, this gap is smaller than 10%, with the exception of Grosseto (a). Grosseto (a) is also the least bright province is more susceptible to being affected by single installations. The gaps for Frosinone (b), Pisa (c) and Macerata (d) can easily by explained by the regular evolution of the street lighting, Figure S8: Examples of not rescaled raw DMSP, calibrated DMSP and VIIRS for different locations in continents. In blue are the DMSP uncalibrated data, in green the DMSP calibrated and in red the VIIRS data. The gap between calibrated DMSP 2010 and VIIRS 2012 can be explained because of the time of observation, first during the evening, second during midnight. This provides information on the amount of variable lights on each continent. Europe (a), North America (b), Asia (c) and Oceania (d) have big gaps, probably because of the turning off of private lighting, although for South America (e) and Africa (f) this gap is much smaller. DMSP 2010 in Europe is an outlier because of the large amount of data with snow in north and mid Europe used for that composite, Figure S9