Introduction

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 surface1. It is one form of (indirect) light pollution2. ALAN has grown to a global phenomenon with recent growth rates on the order of by 3–6% per year3, with further increase to be expected4 by efficient solid state lighting technology5. Recent research on ALAN has found that direct and indirect light pollution can negatively affect flora, fauna and human well-being6.

Ecological light pollution7 (ELP) is the overarching term describing the negative effect of ALAN on the environment. For an introduction, see the review8 and book by Rich and Longcore7. Recent work on ELP include studies on marine turtles9, bats10, insects11, plants12, 13, and microorganisms14. Furthermore, concerns that ELP can affect whole ecosystems and biodiversity have been raised15,16,17, especially for skyglow as it can cause long range effects17. Studies of skyglow in the context of ELP are very sparse, with the exemption of seminal work by Moore on zooplankton18.

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 areas19,20,21, but can also reduce ALAN in rural areas22, 23, as shown by experimental work. Theoretical work is in development but sparse so far24, 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 environment26. 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 monitoring27, 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 brightness29, that recently underwent a major update30. Other spaceborn sources are (individual) images from astronauts (typically from cities) with improved spatial resolution31. 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 resolution22, 23, the conduction of local surveys32, the establishment of local networks33 and global comparison34. While the SQM-L (with a lens) has evolved as a standard device, alternatives like the STARS4ALL night sky brightness photometer have been developed35. 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-uniform36, 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 calibration38. 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 Duriscoe36. 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 Camera41 as well as other custom built camera systems42. 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 mosaic41) or require extensive data processing and or calibration42.

The importance of imaging devices for light pollution measurements has been pointed out earlier36, 43, 44. However, it is not widely recognized in the community investigating ELP, with the exception of work on marine turtle hatchlings45. 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 studies22, 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.

Results

Transect near Montsec Astronomical Park and the Spanish town of Balaguer

Figure 1 shows 8 (of 18 total) luminance maps of the night sky obtained during two transects from the town of Balaguer towards Montsec Astronomical Park in May 2016. The images are shown in Hammer–Aitoff equal area projection and were calculated from all-sky images obtained with a commercial DSLR camera. In the left column (a, c, e, g) of Fig. 1 data from a night with clear sky (May 3–4) is shown. In the right column (b, d, f, h) of 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.

Figure 1
figure 1

Luminance maps of the night sky in Hammer–Aitoff equal area projection calculated from all-sky images obtained. along the transects from the town of Balaguer towards Montsec Astronomical Park in May 2016. The left column (a,c,e,g) shows data obtained on the night of May 3rd/May 4th with a clear sky and the right column (b,d,f,h) shows data obtained atthe night of May 5th/May 6th with partly cloudy sky. The upper row (a,b) was acquired at 2 km distance, the second row (c,d) at 8.5 km distance, the third row (e,f) at 17.7 km distance to the city center of Balaguer. The lower row was obtained at the Astronomical Park of Montsec at 27 km distance to Balaguer. The full data set is shown in the Appendix.

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 Fig. 2, L v,zen is plotted as a function of the distance from the town of Balaguer on a linear a) and on a double logarithmic scale b). The black open circle and solid line show the transect data for the clear night and the red open circles and solid line show the data from the transect for the partly cloudy night. The diamonds in both plots show data taken at the Astronomical park site at 27 km distance to Balaguer (lower row in 1) with black diamonds for clear sky and red diamonds for partly cloudy sky. For the clear sky an extra data point closer to the town center was taken, which is indicated by the black dashed line. Additionally, the zenith luminance values measured with the Digilum luminance meter (Optronik GmbH, Berlin, Germany) are shown for the partly cloudy night (blue crosses and solid line). The individual values from the whole data set are summarized 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/m2)30, please note that there are other values in use36) and the green dashed line indicates a value twice as much (0.50 mcd/m2).

Figure 2
figure 2

Luminance at the zenith as a function of the distance from the city center of Balaguer measured at two differentnights taken at the transects with clear sky (black open circle and solid line) and partly cloudy sky (red open circles and solidline) in a linear scale (a) and a double logarithmic scale (b). Additionally data from the Astronomical park at 27.2 km distance isshown as black (clear sky) and red (partly cloudy sky) diamonds in both plots, respectively. For the clear sky case, an extra datapoint closer to the town center (600 m) was taken, which is indicated by the black dashed line. For the partly cloudy night, comparative data obtained with the Digilum luminance meter is shown in blue. The dashed magenta line indicates the value ofa “typical” clear sky (0.250 mcd/m2) and the green dashed line a value twice as much (0.50 mcd/m2).

Table 1 List of the zenith luminance values at the different locations and the two nights (see methods for details).

For the clear sky case, the \({L}_{v,zen}\) values drop from 1.16 ± 0.12 mcd/m2 at 1 km distance to the center of Balaguer to values as low as 0.20 ± 0.02 mcd/m2 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/m2 (green line), while the value at 8.4 km is only 25% above the 0.25 mcd/m2 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/m2 at 1 km distance to 0.49 ± 0.05 mcd/m2 at the Astronomical Park. Later in the night, with the sky clearing up, \({L}_{v,zen}\) reaches 0.33 ± 0.04 mcd/m2 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 literature37, 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.

Figure 3
figure 3

(a) Horizontal and (b) scalar (hemispheric) illuminance as a function of the distance from the city center of Balaguermeasured for two different nights with clear sky (black) and partly cloudy sky (red). Magenta lines show estimated ideal values (see methods).

Table 2 List of the illuminance values for the two nights at the different locations along the transect (see methods for details).

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/m2 is only 24% above the value of 0.250 mcd/m2, 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/m2. 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 tools22, 36, 43, 44, 46,47,48 have been almost exclusively used by astronomers, physicists, and engineers, with the notable exception of scientists working with marine turtles45. However, spatial information might be vital for understanding the effects of ELP, as some species possibly use spatial features for orientation49. 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 information44. Recent studies have determined the spectral irradiance with spatially single channel sensors of the sky26 and even under water50 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 setting51.

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 recently37.

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 ELP49. 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 scientists52.

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 report53.

Figure 4 shows the “World Atlas” predictions for (clear sky) artificial skyglow30 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.

Figure 4
figure 4

Location of the transect near the Spanish town of Balaguer towards Montsec Astronomical Park: Night skybrightness maps (taken from ref. 30) of (a) Spain, (b) Catalonia and d) the local area. (c) Shows a map of the local area.

Table 3 List of the locations along the transect from Balaguer to Port d’ Àger.

The topography of the area is dominated by the mountains of the Montsec region. Along the transect, the highway C-12 climbs from 273 m elevation (stop 1) to 908 m (stop 7). The dominant sources for artificial skyglow are the light domes of Balaguer (27.2 km from PAM-COU) and Lleida (ca. 46 km from PAM-COU), while the larger light dome of Barcelona (ca. 140 km from PAM-COU) is also apparent near the horizon. According to Falchi et al.30 the calculated clear sky lumnance in the center of Balaguer is 1.07–1.96 mcd/m2 and falls off towards PAM to 0.2–0.23 mcd/m2 (see Fig. 5). However, we want to point out that this global model does not include local characteristics such as the topography or vegetation cover of the landscape, yet.

Figure 5
figure 5

Night sky brightness map (taken from ref. 30) as seen from Parc Astronomic Montsec towards Balaguer.

Fortunately for astronomical observations, the PAM-COU (measurement location 8) is located in a valley surrounded by hills and mountains, usually shielding the artificial skyglow present at the horizon. For the majority of clear nights, the zenith night sky brightness is near-natural with values of 22 mag SQM /arcsec2. On cloudy nights the zenith night sky brightness can either decrease or increase, as shown by an in depth-study linking SQM values and cloud height at PAM-COU23.

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 vision55. 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 campaigns53, 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:

$${E}_{v,hor}={\int }_{0}^{2\pi }{\int }_{0}^{\pi /2}{L}_{v,sky}(\theta ,\varphi )sin\,\theta \,cos\,\theta \,d\theta \,d\varphi ,$$
(1)

with θ being the zenith angle, ϕ being the azimuth angle, and \({L}_{v,sky}(\theta ,\varphi )\), 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}_{v,sky}(\theta ,\varphi )=const\mathrm{.}={L}_{v,zen}\), we can treat \({L}_{v,zen}\) as a constant and do not have to integrate over it, while the integral in equation (1) simply becomes 2π from integration over the azimuth and 0.5 from integration over the zenith angle, giving

$${E}_{v,hor}\approx \pi \cdot {L}_{v,sky,zen}\mathrm{.}$$
(2)

With the commonly used assumption of a typical clear night sky zenith luminance of \({L}_{v,zen,ideal}\) = 0.25 mcd/m2, we can approximately infer a typical illuminance for this case: \({E}_{v,hor,ideal}\approx \pi \cdot 0.25{\mathrm{mcd}/m}^{2}\approx 0.78\) 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 illuminance36:

$${E}_{v,scal,hem}={\int }_{0}^{2\pi }{\int }_{0}^{\pi /2}{L}_{v,sky}(\theta ,\varphi )sin\,\theta \,d\theta \,d\varphi ,$$
(3)

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 report53. The device has been used in other intercomparison measurements before25, 57.