Imaging and mapping the impact of clouds on skyglow with all-sky photometry

Artificial skyglow is constantly growing on a global scale, with potential ecological consequences ranging up to affecting biodiversity. To understand these consequences, worldwide mapping of skyglow for all weather conditions is urgently required. In particular, the amplification of skyglow by clouds needs to be studied, as clouds can extend the reach of skyglow into remote areas not affected by light pollution on clear nights. Here we use commercial digital single lens reflex cameras with fisheye lenses for all-sky photometry. We track the reach of skyglow from a peri-urban into a remote area on a clear and a partly cloudy night by performing transects from the Spanish town of Balaguer towards Montsec Astronomical Park. From one single all-sky image, we extract zenith luminance, horizontal and scalar illuminance. While zenith luminance reaches near-natural levels at 5 km distance from the town on the clear night, similar levels are only reached at 27 km on the partly cloudy night. Our results show the dramatic increase of the reach of skyglow even for moderate cloud coverage at this site. The powerful and easy-to-use method promises to be widely applicable for studies of ecological light pollution on a global scale also by non-specialists in photometry.

Artificial skyglow is the part of artificial light at night (ALAN) that is scattered or reflected within the atmosphere and directed back towards the Earths surface 1 . It is one form of (indirect) light pollution 2 . ALAN has grown to a global phenomenon with recent growth rates on the order of by 3-6% per year 3 , with further increase to be expected 4 by efficient solid state lighting technology 5 . Recent research on ALAN has found that direct and indirect light pollution can negatively affect flora, fauna and human well-being 6 .
Ecological light pollution 7 (ELP) is the overarching term describing the negative effect of ALAN on the environment. For an introduction, see the review 8 and book by Rich and Longcore 7 . Recent work on ELP include studies on marine turtles 9 , bats 10 , insects 11 , plants 12,13 , and microorganisms 14 . Furthermore, concerns that ELP can affect whole ecosystems and biodiversity have been raised [15][16][17] , especially for skyglow as it can cause long range effects 17 . Studies of skyglow in the context of ELP are very sparse, with the exemption of seminal work by Moore on zooplankton 18 . This is mainly owing to the fact that the status-quo of light at night in general and of ALAN (including skyglow) is not well known. Skyglow is highly dynamic with changing atmospheric conditions. Clouds can dramatically increase ALAN in urban areas [19][20][21] , but can also reduce ALAN in rural areas 22,23 , as shown by experimental work. Theoretical work is in development but sparse so far 24,25 . To understand skyglow and ELP in detail, truly interdisciplinary studies are necessary. As a basis for this, it is essential to first fully know the availability of all light at night in the environment 26 . This should include unpolluted sites as well as heavily polluted sites and all states in between, at best at all weather conditions. Skyglow can be determined from the ground as well as from space. Satellite monitoring 27,28 , allows for global coverage with reasonable spatial resolution and revisiting times. In combination with radiative transfer models such data can be used to infer the status of skyglow worldwide, as demonstrated earlier with the world atlas of artificial night sky brightness 29 , that recently underwent a major update 30 . Other spaceborn sources are (individual) images from astronauts (typically from cities) with improved spatial resolution 31 . However, with space-born and airborne methods only the direct upwelling part of ALAN can be measured, while the downwelling part, most relevant for ELP, cannot be accessed directly. Furthermore, atmospheric conditions and especially clouds prevent any satellite data acquisition, rendering it a clear-weather-only method.
Therefore, ground based measurements are still essential to evaluate skyglow and its impact on the environment. Inexpensive single channel devices such as the sky quality meter (SQM, Unihedron, Ontario, Canada) as handheld or permanently installed devices have allowed long-term monitoring with high temporal resolution 22,23 , the conduction of local surveys 32 , the establishment of local networks 33 and global comparison 34 . While the SQM-L (with a lens) has evolved as a standard device, alternatives like the STARS4ALL night sky brightness photometer have been developed 35 . However, with both approaches only the luminance (or radiance) at zenith is measured. The zenith luminance is not a good proxy for the total light available at a site, and it is not straightforward to convert it to horizontal or scalar illuminance, commonly used by biologists, as the night sky brightness is usually non-uniform 36,37 , Commercial luxmeters have the drawback that they are either not sensitive enough or, when calibrated for low light levels, are relatively expensive. The IYA Lightmeter, a cheap and rigid solar-cell based device, could be an affordable alternative but is discontinued and requires self calibration 38 . All single channel devices suffer from the fact, that the spatial distribution and the origin of light pollution cannot be resolved.
Information of the spatial distribution of the night sky brightness is relevant for ELP studies and can be obtained from imaging sensors as nicely outlined by Duriscoe 36 . Several approaches to measure light pollution with imaging sensors exist, like the all-sky transmission monitor (ASTMON) 39,40 , CCD cameras with sophisticated mechanics to built all-sky mosaics like the US National Park Camera 41 as well as other custom built camera systems 42 . However, while these solutions are without question very good options for detailed and precise monitoring of astronomical light pollution on clear nights, they are rather complicated in handling (e.g. by using multiple filters or requiring photometric calibration with celestial objects), relatively slow in data acquiring (e.g. about 40 minutes for an all-sky mosaic 41 ) or require extensive data processing and or calibration 42 .
The importance of imaging devices for light pollution measurements has been pointed out earlier 36,43,44 . However, it is not widely recognized in the community investigating ELP, with the exception of work on marine turtle hatchlings 45 . We assume that this is mainly owing to the fact that the threshold of using the aforementioned systems by biologists in the field is high due to their complexity. Recently, commercial digital single lens reflex (DSLR) cameras with fisheye lenses have been used for night sky studies 22,44,[46][47][48] . These cameras come with reasonable factory calibration and promise to allow data acquisition by non specialists as they do not require special training or complex knowledge about photometry or electronics. The huge public demand for imaging devices in consumer electronics like smartphones, has also resulted in dropped prizes for stand alone cameras, making professional level DSLR cameras cheaper than calibrated luxmeters.
Here we show that off-the-shelf commercial DSLR cameras are well suited for measuring the night sky brightness in the field under overcast as well as clear conditions. The simple and fast data acquisition allowed us to perform two rapid transects from the town of Balaguer towards Montsec Astronomical Park in Spain, covering more than 20 km in about 2 hours. We were able to map the reach of skyglow from the peri-urban area into the remote area. From the data, we can obtain illustrative luminance maps, extract the zenith luminance as well as horizontal and scalar illuminance (see methods), giving a comprehensive set of data that cannot be inferred easily from single channel sensors or satellite imagery, especially not under cloudy conditions. We find that with the presence of clouds, areas that reach near-natural zenith brightnesses under clear sky conditions experience the same light levels as peri-urban area under clear conditions. Such a shift of anthropogenic skyglow into pristine areas has not been investigated at this level of detail with imaging sensors. We judge that the method is widely applicable for field studies up to the global level as it relies on off-the-shelf equipment, making it an ideal choice for non-specialists in photometry.  Fig. 1 data from a night with partly cloudy sky (May 5-6) is shown. The upper row (a, b) shows data taken at 2 km distance, the second row (c, d) data taken at 8.5 km distance and the third row (e, f) data taken at 17.7 km distance to the town center of Balaguer. The lower row shows data taken at Astronomical Park of Montsec at 27 km distance to Balaguer. The full data set, consisting of 18 all-sky luminance maps, is shown in the Supplement. The measurement locations, topography and other artificial skyglow data from the region are described in the methods section.

Transect near Montsec Astronomical Park and the Spanish town of Balaguer.
From the clear sky luminance map obtained at 2 km distance shown in Fig. 1a), the light dome originating from the skyglow of the town Balaguer is clearly visible as a big hump on the right. The skyglow from the city of Lleida (24.8 km distance to this position) is apparent by a smaller hump in the center of the image (left of the Balaguer). Within the clear sky data set, these two sky domes diminish with progressing distance and merge into one broader hump confined to the horizon.
From Fig. 1, the impact of the clouds on the luminance of the whole sky becomes immediately apparent. For all distances, the partly cloudy night (right column) has a higher luminance than the corresponding clear night (left column). Furthermore, the impact of the sky dome from Lleida seems to have a higher impact for cloudy than for clear nights, as can be seen by a spot-like feature in the center of Fig. 1b and d). Furthermore, some of the skyglow that is masked by e.g. mountains on clear nights (as it is confined to the horizon) becomes apparent for partly cloudy conditions. This can be seen when examining the lowest row in Fig. 1 and comparing the luminance distribution on the left side of the images. Also, without any further data analysis, Fig. 1a) appears most similar to Fig. 1f), with the latter appearing slightly darker.
From the all-sky images, the zenith luminance, L v,zen , was calculated (see methods). In    Table 1. To guide the eye, the dashed magenta line indicates the value often referred to be the luminance of a "typical" clear sky (0.250 mcd/m 2 ) 30 , please note that there are other values in use 36 ) and the green dashed line indicates a value twice as much (0.50 mcd/m 2 ).
For the clear sky case, the L v zen , values drop from 1.16 ± 0.12 mcd/m 2 at 1 km distance to the center of Balaguer to values as low as 0.20 ± 0.02 mcd/m 2 for the Astronomical Park at 27 km distance. The value of 0.62 ± 0.06 measured at 2 km is already in the vicinity of 0.50 mcd/m 2 (green line), while the value at 8.4 km is only 25% above the 0.25 mcd/m 2 reference (magenta line).
With the presence of clouds, L v zen , values increase up to a factor of 6 at the individual sites (see Table 1 right hand column), with L v zen , dropping from 3.27 ± 0.33 mcd/m 2 at 1 km distance to 0.49 ± 0.05 mcd/m 2 at the Astronomical Park. Later in the night, with the sky clearing up, L v zen , reaches 0.33 ± 0.04 mcd/m 2 at the latter, most distant site. It is interesting to note, that all values obtained along the transect (stop 1-7) are higher than that obtained at 2 km distance for the clear night. Furthermore, the L v zen , value at 5 km distance is as high as measured inside of the town (600 m to center) for the clear night.
It is further possible to calculate horizontal illuminance, E v,hor , and scalar illuminance, E v,scal,hem , (here for the hemisphere) in the commonly used unit lux from the all-sky data, as the zenith luminance might not be representative for all light incident at a specific site (see methods for detail). Figure 3a) shows the horizontal illuminance for clear conditions (black) and partly cloudy conditions (red), with a hypothetical ideal natural-sky illuminance value of E v,hor,ideal ≈ 0.78 mlux 30 (again other values are reported in literature 37 , see methods and discussion) indicated by the horizontal dashed magenta line. The trend is the similar as for the zenith luminance but weaker. For example, the ratios between cloudy and clear nights are all lower than for the zenith luminance. Here, the value of E v,hor at 2.2 km distance for the clear night is in the same range as E v,hor at 13.5 km for the cloudy night. Figure 3b) shows the scalar illuminance as a function of distance from the center of Balaguer, again for clear conditions (black) and partly cloudy conditions (red) as well as with the ideal scalar illuminance, E v,scal,hem,ideal ≈ 1.56 mlux, indicated by the horizontal dashed magenta line. The overall trend is about the same as for the horizontal illuminance. The individual values for horizontal and scalar illuminance and the ratios between clear and cloudy night are summarized in Table 2.

Discussion
One of the big advantages of this method is that it makes it possible to compare different physical quantities acquired simultaneously in a single image. The photometric data from the all-sky images allowed us to extract zenith luminance, horizontal illuminance and (hemispheric) scalar illuminance simultaneously. Along the transects, we could investigate the impact of clouds on the amplification of skyglow at individual sites, and map the increase in the reach of skyglow. The most dramatic change in individual site brightening and extension of the reach of skyglow was observed for the zenith luminance levels, while for the illuminance measurements this shift was lowered. Zenith luminance observed on clear nights can underestimate the degree to which a site is exposed to skyglow under all weather conditions. While the zenith appears near-natural, there can be significant skyglow at the horizon. With the presence of clouds, skyglow can "creep" into the field of view of the measurement device. This is apparent for example for the brightness values obtained at 8.4 km distance to Balaguer in 1c and d). At this site, the zenith luminance for the clear night L v,zen = 0.31 ± 0.03 mcd/m 2 is only 24% above the value of 0.250 mcd/ m 2 , appearing as almost natural sky. The horizontal illuminance at the clear night E v,hor = 1.9 ± 0.2 mlux is already 2.4 times as much as the reference value of 0.78 mlux, while the (hemispherical) scalar illuminance E v,scal,hem is more 3.1 times the dark sky value of 1.56 mlux. The amplification factor for clouds is 4.9 for the zenith luminance, resulting in a L v,zen value 6.1 times brighter compared to the natural sky value of 0.250 mcd/m 2 . For the horizontal illuminance, the amplification factor is 3.7, resulting in an 8.9 times higher value of E v,hor with respect to the reference case. For E v,scal,hem the amplification factor is 3.2 resulting in 9.9 higher value compared to the ideal clear sky case.
In the context of ELP, there are hardly any animals that just sense the zenith brightness alone. Thus illuminances are a better proxy for ELP. Therefore, the amount of light available for ecological effects is underestimated by a factor of 2 and 2.5 by SQMs (or similar devices) for the clear sky case, when inferring illuminance levels from zenith brightness compared to measured illuminances. On the other hand, the brightening effects of clouds can be overestimated by single channel devices.
For all three photometric quantities, the reach of skyglow was enhanced with the presence of clouds. An extreme shift was observed in particular for the zenith luminance L v,zen , with higher luminance observed at 22.2 km from the town for the partly cloudy night than were seen at 2.2 km for the clear night. For the horizontal illuminance, this "distance-effect" is slightly reduced. For the partly cloudy night, E v,hor was higher at 13.5 km than was observed at 2.2 km for the clear night. For E v,scal,hem , the size of the difference was further reduced.
The other big advantage of the method is its potential wide applicability. So far, all-sky imaging tools 22, 36, 43, 44, 46-48 have been almost exclusively used by astronomers, physicists, and engineers, with the notable exception of scientists working with marine turtles 45 . However, spatial information might be vital for understanding the effects of ELP, as some species possibly use spatial features for orientation 49 . We assume that wide application of imaging tools was so far hampered by the complexity of existing methods. Operation of the commercial DSLR camera is straightforward, and well suited for fast acquisition during fieldwork with a minimum of training necessary. The camera is mobile, yet rigid and reliable, and we have so far never experienced any equipment failure. In our case, reaching a vantage point in the hilly terrain from the road sometimes took longer than the data acquisition itself. The data analysis is not too complicated, and the software as well as the code used here is freely available. However, one drawback is that the light level values measured are not currently immediately available as they would be with some other instruments.
While this work focuses on spatially resolved luminance measurement of the whole hemisphere, the method can be extended to use the three color channels of the camera to obtain further spectral information 44 . Recent studies have determined the spectral irradiance with spatially single channel sensors of the sky 26 and even under water 50 showing the dramatic alteration of the night light climate. Other recent work includes the use of hyperspectral imaging systems to determine the (direct) light available in an urban setting 51 .

Conclusions
We have shown that DSLR cameras with fisheye lenses are well suited for all-sky field studies in the context of ELP under both clear and cloudy conditions. By performing two transects from a peri-urban area into a remote area, we observed a shift in the reach of skyglow under clear and partly cloudy conditions. The data support the hypothesis that clouds can dramatically extend the reach of skyglow, resulting in ALAN penetrating areas that appear pristine on clear nights, and filling the gaps between illuminated areas. At our study site, zenith luminance levels at 22 km from a small city on a cloudy night are similar to those at 2 km on clear night. This trend is also observed in the more ELP-relevant illuminance levels, with slightly lower amplification and extension. We further find evidence that in some scenarios, the zenith brightness is not a good proxy for the light available at a site, supporting concerns raised recently 37 .
The strength of all-sky imaging is that several photometric quantities can be extracted from one image, and that individual parts of the night sky can be investigated. This might be essential for studies on animal behavior in the context of ELP 49 . This method uses off-the-shelf equipment, and could be easily adopted by the ELP community. This would hopefully result in increased and global data on the state of skyglow under cloudy as well as clear conditions. Single site studies would be useful, as well as other transect data. The method may also lend itself to long-term monitoring, and possibly data acquisition by citizen scientists 52 .

Methods
Study site. Transect from Balaguer to Montsec Astronomical Park, and the status of artificial skyglow in the area. During the 2016 Stars4All/LoNNe intercomparison campaign held at Montsec Astronomical Park (Parc Astronómic Montsec, Centre d' Observació de l' Univers: PAM-COU), two transects were conducted on a clear night (May 3rd to May 4th) and a partly cloudy night (May 5th to May 6th). Details on the intercomparison campaign, including other activities and other equipment used, can be found in the campaign report 53 . Figure 4 shows the "World Atlas" predictions for (clear sky) artificial skyglow 30 for a) the Iberian Peninsula b) Catalonia and surroundings, and d) the area of the transect. In Fig. 4c) a local map is shown. The transects were started at the outskirts of Balaguer (41.791100N, 0.797494E, 273 m elevation) and were conducted along the highway C-12 leading North towards Àger and PAM-COU, respectively. The last stop on the route was at Port d' Àger: (41.979070N, 0.750630E, 908 m elevation). In between these two points, measurements at five further stops have been made for each of the two nights. The seven stops are labeled with numbers in the map. Additional indicators in Fig. 4c and d) show the center of Balaguer and the location of the PAM-COU, where additional measurements were made synchonously with the transect. Details on all measurement sites can be found in Table 3.
The topography of the area is dominated by the mountains of the Montsec region. Along the transect, the  Weather conditions. On 3 May 2017, the sky was clear during the day with some cirrus, still significant around sunset but disappearing by around nautical twilight. No clouds were picked up by ceilometer at PAM after 10:20 pm local time (8:20 pm UTC). There were no clouds at all by astronomical twilight, and this condition lasted through the entire night, including during the morning twilight. At the end of the night, the humidity level increased to 78% (PAM-COU Weather Station).
On May 5th, the sky became progressively cloudier throughout the day, with high clouds. The entire evening was characterized by high thin clouds. At PAM, stars were visible, and there was a moisture layer up to the top of the mountain which allowed scattered light from very remote cities to be seen. The scattering from upward-directed lights from the villages to the North of the mountain was conspicuous. The sky was covered the whole night with medium and high clouds. The ceilometer measured cloud heights around 6 km for the relevant measurement period. Early in the morning low clouds also appeared, according to ceilometer measurements. DSLR camera. The all-sky images were obtained with a commercial DSLR camera (Canon EOS 6D). This camera has a full-frame CMOS sensor with 20.2 Megapixel (5496 × 3670 pixels) and allows ISO settings from 100 to 25,600 and shutter speed ranges of 1/4,000 to 30 s without extra equipment. It has a built in GPS sensor. The camera was operated with a circular fisheye lens (Sigma EX DG with 8 mm focal length). The lens was focused to a bright celestial object in the live view of the camera. If celestial objects are masked by clouds, focusing can be done at an artificial light source at large distance (e.g. on the horizon).
To be both mobile in the field yet acquire still images, the camera was mounted on a tripod. To acquire all-sky images, the camera was aligned with the center of the lens oriented towards the zenith. For the light levels present during the two nights, ISO 1600 was used for images closer to the town and ISO 3200 for images further away from the town. The shutter speed was varied between 15 s and 30 s along the transect, while a long exposure of 120 s was used at the Montsec Astronomical Park location. For the cloudy night, several images were obtained and averaged.
The method applied here of using commercial DSLR cameras with fisheye lenses is simple as well as robust, although not as rigid as a commercial outdoor or even waterproof light measuring devices like luxmeters. The point and click ability offers great opportunities, especially if measurements have to be acquired quickly.  Furthermore, commercial DSLR cameras are now quite affordable (in the same price range as a calibrated luxmeter or PAR sensor). The spatially resolved data allow for the extraction of several parameters acquired at the same time (see below).
DiCaLUM software. Image processing was performed using DiCaLum (Version 1.01, by Zoltán Kolláth) (for a brief description see ref. 54). The code is programmed in GNU Octave and relies on the free software DCRaw. It can be obtained from the authors directly. The RAW images are read out, and the radiance or luminance can be obtained from the pixel values of the three different channels. In the present version, the software uses the green channel, which has significant overlap with the V(λ) curve referenced to human photopic vision 55 . The pixel values from the green channel are then used to map the luminance values for each individual pixel. Although highly standardized, DSLR cameras and lenses still require thorough calibration, especially for the vignetting of the fisheye lens. DiCaLum has several cameras and fisheye lenses in its library. Initial calibration was done in the laboratory, while some cameras were cross calibrated in several intercomparison campaigns 53,56 , For further reading on camera calibration we refer to refs 44, 47, where the method is described in more detail.
A common measure to quantify how much light is available in the environment is horizontal illuminance, the total luminous flux incident on a flat horizontal surface (here the Earth's surface), per unit area. When assuming the sky to be a hemisphere, the horizontal illuminance is defined as: , being the luminance of the sky. The SI unit is lux. The horizontal illuminance is sometimes also referred to be cosine corrected. Commonly, the horizontal illuminance is inferred from the zenith luminance, L v sky zen , , . When assuming a uniform sky luminance θ φ = .= L c onst L ( , ) mlux. As discussed in more detail in ref. 37, this might not be a valid approximation for most real scenarios, as the sky brightness is rather nonuniform (evident in our own data, see Fig. 1).
The horizontal illuminance may not be the best measure for all organisms affected by ecological light pollution, as for some species the angle of incidence may not be of primary importance. For this case, the total luminous flux incident on a small spherical surface gives the so-called scalar illuminance. For the sky and the all-sky images, it makes sense to define a hemispherical scalar illuminance 36 : v scal hem v sky , , 0 2 0 /2 , which does not incorporate the cosine correction. The DiCaLum code can produce luminance maps, extract zenith luminance, and horizontal and (hemispherical) scalar illuminance, by summing the solid angle weighted luminance values with or without additional cosine correction.
Digilum luminance meter. The Digilum (Optronik GmbH, Germany) is a custom luminance meter specifically tailored and calibrated for low light level detection at the zenith. It has a 5° aperture and a spectral response that is well matched to photopic vision, V(λ). The Digilum was mounted on the back of the car for the second transect during the partly cloudy sky and is constructed on a gimbal, guaranteeing straight pointing to the zenith. Averaged luminance values were obtained for the measurement duration of the corresponding DSLR imaging. A detailed description of dark current subtraction is given in the Stars4All intercomparison campaign report 53 . The device has been used in other intercomparison measurements before 25,57 .