The Balkan Peninsula has rich biodiversity with a large number of endemic species (Griffiths and Kryštufek 2004; Petrova and Vladimirov 2010; Tomović et al. 2014; Vuksanović et al. 2016). That is why a part of its territory in borders of Mediterranean has been recognized as a World Biodiversity Hotspot (Mittermeier et al. 1999, 2005; Myers et al. 2000). There are 6 biogeographic regions located within the Balkans: Alpine, Black Sea, Continental, Mediterranean, Pannonian, and Steppic (European Commission 2019; European Environment Agency 2012). Despite nature conservation efforts and development of nature conservation networks like the Important Bird Area Network (Puzovič 2009; Velevski et al. 2010), the Important Plant Area Network (Angelova et al. 2012; Melovski et al. 2012; Nikolić 2009; Petrović 2009; Stevanović and Šinžar-Sekulić 2009), the NATURA 2000 Network (Gussev and Tzonev 2015; The Ministry of Environment and Water of Bulgaria 2014), and Emerald Network (Convention on the Conservation of European Wildlife and Natural Habitats 2008; Group of Experts for the setting up of the Emerald Network of Areas of Special Conservation Interest 2003, 2006; Prokić 2008) in countries of the region, anthropogenic influence to natural and semi natural ecosystems has been increasing for a long time (Griffiths and Kryštufek 2004; Kulakowski et al. 2017; Longman et al. 2018; Marinova et al. 2012; The European Environment Agency 2011). Moreover, new types of disturbance and pollution arise, one of the more recent pressures is artificial light at night (ALAN) that has already been noted for protected areas of autonomous province of Vojvodina in North Serbia (Bjelajac and Đerčan 2019).

Nowadays there is no exact data about the level of ALAN from which the impact on biodiversity starts (Falchi et al. 2016), but a lot of evidence about ALAN influence on reproduction, navigation, foraging, habitat selection, communication, trophic, and social interactions of the biota have already been discussed (Bennie et al. 2016; Dominoni et al. 2016; Gaston and Bennie 2014; Hölker et al. 2010b; Longcore and Rich 2004; Navara and Nelson 2007; Rich and Longcore 2006). For example, ALAN has been documented as a contributing factor to the global decline in insect populations (Grubisic et al. 2018; Hallmann et al. 2017; Leather 2018; Macgregor et al. 2015). Furthermore, ecological light pollution not only can cause cascading effects in ecosystems, restructuring ecological communities by modifying the interactions between species, and impacting pollination and seed dispersal (Bennie et al. 2015a), but also being a driver of evolution across urban-rural landscapes (Hopkins et al. 2018). Unfortunately, understanding of benefits which human society has from ALAN (Boyce 2019; Gaston et al. 2015b), as well as its dynamics (Cinzano et al. 2001; Falchi et al. 2016, 2019; Kyba et al. 2017) does not allow considering this disturbance factor as temporary or insignificant for species, especially nocturnal ones.

In this case, strict protected areas have a special value for biodiversity conservation because the impact of different anthropogenic factors is minimized there. It must concern ALAN too; however, it is not taken into account by either scientists or authority in the majority of countries. There are only few countries like Italy, which enforce laws against light pollution (Falchi 2018). Though the investigations of the level of ecological light pollution within protected areas of some countries and regions show its widespread incursion (Bjelajac and Đerčan 2019; Gaston et al. 2015a; Guetté et al. 2018; Jiang et al. 2017; Peregrym et al. 2018). Therefore, in view of the foregoing, it has been decided to estimate the level of ecological light pollution and its danger for strict protected areas (Nature Reserves, National and Nature Parks) of the Balkans as unique places for biodiversity conservation.

1 Material and Methods

The research covers strict protected areas of three Balkan countries: the Republic of Serbia, the Republic of Bulgaria, and Montenegro. The Republic of Serbia is considered together with Kosovo region because the Republic of Kosovo is a partially recognized state and disputed territory still, though the nature management of Šar planina National Park which is located in this region is carried out by Kosovo government at present. These three countries have been selected because of their different political and economic status that is important in the context of the dependence of the level of ecological light pollution from economic development (Bennie et al. 2015b) and proceeding from data availability. GIS layers showing the borders of strict protected areas for Serbia were computed via site Geosrbija (https://a3.geosrbija.rs) on the basis of Map of Protected areas of Serbia, available through http://serbia.gdi.net/zzps/. The data for strict protected areas of Bulgaria was uploaded from the website of the Executive Environment Agency, the Ministry of Environment and Water of the Republic of Bulgaria (http://eea.government.bg/zpo/bg/index_download.jsp). GIS layers with borders of Montenegrin National Parks were received by personal communication with these organizations; however, it is available schematically via https://nparkovi.me/mapa-parkova/

The objects of our study were 13 National Parks, 11 Nature Parks and 55 Reserves: 5 National Parks (Đerdap, Fruška Gora, Kopaonik, Šar Planina, Tara) and 29 Nature Reserves (Bagremara, Carska bara, Deliblatska peščara, Goč-Gvozdac, Gornje Podunavlje, Jelašnička klisura, Jerma, Karađorđevo, Klisura reke Mileševke, Klisura reke Trešnjice, Koviljsko-Petrovaradniski Rit, Kraljevac, Kukavica, Ludaško jezero, Mala jasenova glava, Obedska bara, Okanj bara, Paljevine, Pašnjaci velike droplje, Peštersko polje, Ritovi Donjeg Potisja, Selevenjske pustare, Slano Kopovo, Suva planina, Tesne jaruge, Titelski breg, Uvac, Venerina padina, Zasavica) in Serbia; 3 National Parks (Central Balkan, Pirin, Rila), 11 Nature Parks (Belasitsa, Bulgarka, Golden Sands, Persina, Rila Monastery, Rusenski Lom, Shumen Plateau, Sinite Kamani, Strandzha, Vitosha, Vrachanski Balkan), and 26 Nature Reserves (Ali Botush, Beglika, Beli Lom, Byala Krava, Chervanata Stena, Chuprene, Dupkata, Elenova Gora, Gorna Topchiya, Gornata Koria, Kaliakra, Kamenshchitza, Kastrakliy, Kamchia, Kazanite, Kupena, Leshnitsa, Mantaritza, Orelyak, Orlitzata, Ropotamo, Sokolata, Soskovcheto, Tisata, Tserna Reka, Velchi Dol) in Bulgaria; 5 National Parks (Biogradska gora, Durmitor, Lovćen, Prokletije, Skadarsko jezero) in Montenegro (Fig. 1). It is worth noting that there are 55 Nature Reserves in Bulgaria, but 29 of them are located within borders of National or Nature Parks, ; therefore, these territories have not been considered separately to avoid duplicating the results. Calculated data for these 29 areas are included in data for National and Nature Parks in a whole.

Fig. 1
figure 1

Studied strict protected areas in Bulgaria, Montenegro, and Serbia

According to accepted biogeographic zoning (European Commission 2019; European Environment Agency 2012), they are located in the next way: 7 National Parks, 4 Nature Parks, and 15 Nature Reserves in Alpine region; 2 Nature Parks and 3 Nature Reserves in Black Sea region; 3 National Parks, 5 Nature Parks, and 21 Nature Reserves in Continental region; 2 National Parks in Mediterranean region; 1 National Park and 16 Nature Reserves in Pannonian region (Fig. 1).

The study has been carried out using available tools from Google Earth Pro (version 7.3.2.5776; https://www.google.com/earth/). We used the New World Atlas of Artificial Sky Brightness in the form of a kmz (Keyhole Markup language Zipped) layer which was created by Falchi et al. (2016) and is available through its 3D Globe version (https://cires.colorado.edu/Artificial-light). We overlaid the GIS layers of the borders of the protected areas with the artificial sky brightness layer and counted the number of squares of each index of level of artificial sky brightness according to the legend of the atlas (Falchi et al. 2016).

2 Results and Discussion

The results of area calculations of different levels of artificial sky brightness within protected territories are presented in Table 1 for Serbia, in Table 2 for Bulgaria, and in Table 3 for Montenegro. Total results for different types of protected areas of three countries have been summarized in Table 4, as well as total results for different biogeographic regions are given in Table 5. To quantify an error within the calculations, we have added two columns to each table, one column with the calculated area, and the other one with the official area (according to information from the Institute for Nature Conservation of Serbia (Brusin 2018), the Ministry of Environment and Water of Bulgaria (https://www.moew.government.bg/bg/priroda/zastiteni-teritorii/obsta-informaciya-za-zastitenite-teritorii/), and the Public Enterprise for National Parks of Montenegro (http://nparkovi.me/dokumenti/) for every protected area). The highlighted discrepancy is generally not more than 3–5% for studied territories in Bulgaria and Montenegro, except a few mountain ones like Rila Monastery Nature Park (7.83%). The common situation with results of ecological light pollution levels calculations for Serbian protected areas is the same, but there are several ones (Paljevine, Tesne jaruge, Jelašnička klisura and some others) with very small areas (less than 2.64 km2), where the measurement errors are significantly higher because of rounding numbers, the possible inaccuracies in borders of protected areas, as well as the Google Earth measuring capabilities. Besides, calculated data for two nature reserves (Bagremara and Selevenjske pustare) is respectively 3.5 and 2.75 times bigger than official areas of these objects. That is connected with zoning of the territories. Official data include only areas of their strict protected zone, but the available GIS layer with borders of these protected areas does not contain information about zones. Therefore, we were not able to count the level of ecological light pollution only for strict protected zones there, and we have done it for whole territories of these nature reserves. Nevertheless, such big differentiation, in this case, does not have a significant impact on the total result, because total areas of Bagremara and Selevenjske pustare Nature Reserves are relatively small. In light of the above stated, the differentiation between total calculated area and total official area of all strict protected areas of the Republic of Serbia is 0.60%, for the Republic of Bulgaria, 0.86% and for Montenegro, 0.20%. If it is considered for different types of protected areas, then the differentiation between total calculated area and total official area for National Parks is 0.01%, for Nature Parks, 0.25%, and for Nature Reserves, 0.27%. The similar situation is in biogeographic regions: Alpine region, 1.15%; Black Sea region, 0.83%; Continental region, 0.74%; Mediterranean region, 0.12%; Pannonian region, 1.50%.

Table 1 Areas with different levels of artificial sky brightness in the National Parks and Nature Reserves of Serbia
Table 2 Areas with different levels of artificial sky brightness in the Nature Reserves, National, and Nature Parks of Bulgaria
Table 3 Areas with different levels of artificial sky brightness in the National Parks of Montenegro
Table 4 Total results of area calculation of different levels of artificial sky brightness for strict protected territories of Serbia, Bulgaria, and Montenegro
Table 5 Total results of area calculation of different levels of artificial sky brightness for strict protected territories according to biogeographic regions within Serbia, Bulgaria, and Montenegro

According to Falchi et al. (2016), a night sky can be considered “pristine”, if its level of artificial brightness is no more than 1% or < 1.74 μcd/m2. There are no such areas as in studied countries and as within the Balkans as a whole that we can set from the 3D Globe version of the New World Atlas of Artificial Sky Brightness even visually. This fact alone allows saying about the widespread incursion of ALAN within strict protected areas in all Balkan countries. However, let us consider it below on separated examples using obtained data (Tables 1, 2, 3, 4, 5, Figs. 2, 3).

Fig. 2
figure 2

Distribution of strict protected area types in respect to different levels of artificial sky brightness in Serbia, Bulgaria, and Montenegro

Fig. 3
figure 3

Distribution of territories with different levels of artificial sky brightness in respect to biogeographical regions within Serbian, Bulgarian, and Montenegrin strict protected areas

As seen from Tables 1, 2, and 3, the lowest level of artificial sky brightness for strict protected areas of Serbia, Bulgaria, and Montenegro is > 6.96–13.9 μcd/m2. There are 1.85% from the total area of national parks and nature reserves in Serbia with such level of ALAN, 18.83% from the total area of nature reserves, national and nature parks in Bulgaria, and 22.64% from the total area of Montenegrin National Parks. The cleanest protected areas in the context of ALAN pollution are Strandzha Nature Park, Nature Reserves “Dupkata,” “Beglika,” “Gornata Koria,” and “Kazanite” in Bulgaria, as well as National Park “Durmitor” in Montenegro. Besides, Đerdap National Park in Serbia has a plot with a low level of light pollution; however, there are also strong polluted ones at the same time. It is the last level when a night sky can be considered relatively clean, because a night sky with artificial brightness in > 13.9–27.8 μcd/m2 is polluted on an astronomical point of view (Falchi et al. 2016). There are 1200.64 km2 or 13.74% of strict protected areas of Serbia, Bulgaria, and Montenegro (Table 4, Fig. 2) with artificial sky brightness in > 6.96–13.9 μcd/m2. On one hand, it is a small part of total territory, but, on the other hand, that is a unique situation for Continental Europe if Eastern Europe (Peregrym et al. 2018, 2019) and the Scandinavian Peninsula are not taken into account, because the similar territories are practically absent there (Bennie et al. 2015a, b, c; Bennie et al. 2015b, c; Falchi et al. 2016, 2019; Gaston et al. 2015a).

Nevertheless, strict protected areas of Serbia and Montenegro are very strongly polluted by ALAN generally (Fig. 2). So, there are 60.25% of strict protected areas of Serbia with artificial sky brightness more than 55.7 μcd/m2 and the maximum value of 1780 μcd/m2. It is 47.51% for Montenegro, but with the maximum value of 890 μcd/m2. The current situation in Bulgaria is different. There are only 19.1% of Bulgarian strict protected areas with artificial sky brightness more than 55.7 μcd/m2 and the maximum value of 1780 μcd/m2. Moreover, though there is almost 40% which are polluted on an astronomical point of view, potentially these can be cleared using a right policy in combating among ecological light pollution (Dick 2014, 2018; Hölker et al. 2010a).

Considering a potential “cleaning up” in strict protected areas from ALAN separately except Bulgarian National Parks, Rila Monestry Nature Park, and some Nature Reserves, the most perspective territories are National Parks “Prokletije” and “Biogradska gora” in Montenegro, as well as Tara National Park in Serbia. That is because they are located far enough from big settlements with developed infrastructure. Unfortunately, at present, it is complicated to wait for significant success in decreasing of ecological light pollution in National Parks like “Lovćen” and “Skadarsko jezero” because they are surrounded by cities and villages oriented to touristic business with active night life.

The result of data analysis from Table 4 shows that as a percentage the most ALAN polluted one are Nature Reserves among different types of strict protected areas of Serbia, Bulgaria, and Montenegro. There are 67.4% of them with artificial sky brightness more than 55.7 μcd/m2 and the maximum value to 1780 μcd/m2. At the same time, it is 28.5% and 29.94% for national and nature parks, respectively. However, based on areas, the most polluted territories are national parks because there are 1268.18 km2 of their total area with artificial sky brightness more than 55.7 μcd/m2 and the maximum value to 1780 μcd/m2. This indicator is not significantly lower for nature reserves and nature parks—1047.51 km2 and 824.04 km2, respectively.

As seen from Table 5 and Fig. 3, the representation of Serbian, Bulgarian, and Montenegrin strict protected areas in biogeographic regions is different; therefore, it is not correct to compare them. However, it is clear that the most polluted strict protected areas by ALAN are located in the Mediterranean and Pannonian regions. Here is the lowest level in > 27.8–55.7 μcd/m2 for Pannonian region and in > 55.7–111 μcd/m2 for Mediterranean one. It is easy to explain because all seashore regions are strongly polluted by ALAN due to their attraction for tourists and transport infrastructure (Bennie et al. 2015c). Pannonian region within Serbia has a large number of population as well as some industry districts which provide the high level of ecological light pollution there. It was logical to expect the similar data for Black Sea region and the Mediterranean region which are quite popular among tourists. However, within Strandzha Nature Park, which is the biggest strict protected area and the cleanest one from ALAN in the region, there is the lowest level of ecological light pollution among all biogeographic regions represented in studied countries of the Balkans. It is a very positive moment because Strandzha is a unique place which considered to be a refugia of thermophilic plant species in the last glacial time (Patronov 2005). Also a relatively low level of ecological light pollution is in Alpine region that is connected with poorly developed infrastructure in high mountains. Nevertheless, there is 50.71% of strict protected areas polluted by ALAN on an astronomical point of view, and 42.82% has higher level of this pollution. The current situation is significantly worse in Continental region. Here is only 2.64% of strict protected areas which are relatively clean, with artificial brightness in range > 6.96–13.9 μcd/m2. Though there is a big potential for decreasing the level of ecological light pollution in this region, since 30.51% or 608.90 km2 has artificial brightness in range > 13.9–27.8 μcd/m2.

3 Conclusions

Thus, there is no doubt that ALAN significantly influences biota and its habitats in the Balkans, including strict protected areas. It is a new threat for nature conservation not just in Northern Serbia as mentioned earlier (Bjelajac and Đerčan 2019), but in the whole region. At least, this is confirmed by our results of the investigation for Serbia, Bulgaria, and Montenegro. The current research shows only a static situation for strict protected areas of Balkan countries because the New World Atlas of Artificial Sky Brightness (Falchi et al. 2016) does not allow making any dynamical estimation. However, taking into account the pace of ecological light pollution in the world increasing by a degree of ~ 2% per year (Kyba et al. 2017), and the uniqueness of the Balkan Peninsula, first of all, it is very important to save dark skies in this region wherever possible, as well as to start combating for decreasing of ecological light pollution level where it is more perspective. Such activity has to be carried out within strict protected areas firstly, and has to be directed to both practical and educational aims. More difficult goal is to expand dark territories around these strict protected areas. Certainly this will be very hard to realize because of fast infrastructure development in these countries which is apparently often accompanied by increasing area of polluted territories, as well as the level of light pollution. However, the ideal aim is to build a regional network of dark areas which will join territories without light pollution or with its minimal level. It is needed for removal of any biogeographical barriers which appeared due to development of settlement light systems, since they influence on gene drift in populations significantly (Hopkins et al. 2018).

We think that changes in approaches to street light and lighting of buildings within protected areas and their surroundings, according to accepted recommendations (Dick 2014, 2018; Kolláth et al. 2016), have to be among the first practical steps. The next stage is creation of special buffer zones around strict protected areas for decreasing the ALAN impact from light sources of settlements and their infrastructure as well as starting taking into account ALAN effects during development of conservation management plans for protected areas and creation of new ones.

However, these actions must be simultaneous with designing and implementation of educational strategies, because nowadays even scientific popular information about the ALAN problem are limited. Also, it is very important to show to local communities the benefits from saving the dark sky and its impact on well-being, and as a result, conservation activities should be more effective. It is possible to do it by international collaboration, as an example in the framework of the International Dark Sky Places conservation program (http://darksky.org/idsp/) which has been initiated by the International Dark-Sky Association since 2001 (Barentine 2016). Fulfilling the requirements for International Dark Sky Places should provide benefits for both biodiversity conservation and tourism within protected areas. Such a positive experience has been documented in many countries, for example, in adjacent Hungary (Gyarmathy and Kolláth 2017).

Finally, continuous gathering of scientific facts about mechanisms of ALAN influence on biodiversity and protected areas is very important. Balkan countries could be a profitable place for it because of their unique biodiversity and habitat richness on a relatively small area.