Figures
Abstract
Habitat demands and species mobility strongly determine the occurrence of species. Sedentary species with specific habitat requirements are assumed to occur more patchy than mobile habitat generalist species, and thus suffer stronger under habitat fragmentation and habitat deterioration. In this study we measured dispersal and habitat preference of three selected butterfly species using mark-release-recapture technique. We used data on species abundance to calculate Species Distribution Models based on high-resolution aerial photographs taken using RGB / NIR cameras mounted on a UAV. We found that microhabitats for species with specific habitat requirements occur spatially restricted. In contrast, suitable habitats are more interconnected and widespread for mobile habitat generalists. Our models indicate that even managed grassland sites have comparatively little habitat quality, while road verges provide high quality micro-habitats. In addition, dispersal was more restricted for specialist butterfly species, and higher for the two other butterfly species with less ecological specialisation. This study shows synergies arising when combining ecological data with high precision aerial pictures and Species Distribution Models, to identify micro-habitats for butterflies. This approach might be suitable to identify and conserve high quality habitats, and to improve nature conservation at the ground.
Citation: Habel JC, Teucher M, Rödder D (2018) Mark-release-recapture meets Species Distribution Models: Identifying micro-habitats of grassland butterflies in agricultural landscapes. PLoS ONE 13(11): e0207052. https://doi.org/10.1371/journal.pone.0207052
Editor: Petr Heneberg, Charles University, CZECH REPUBLIC
Received: June 21, 2018; Accepted: October 9, 2018; Published: November 28, 2018
Copyright: © 2018 Habel et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Data Availability: All data underlying the study are within the paper and its Supporting Information files.
Funding: The authors received no specific funding for this work.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Habitat configuration and habitat quality significantly influence the occurrence and persistence of species [1]. Today, many remaining habitats are small and geographically isolated from each other, as most semi-natural calcareous grasslands across agricultural landscapes of Central Europe [2]. Recent studies showed that continuing fragmentation of remaining habitats and further intensification of agriculture in between lead to barriers for many species [3]. Furthermore, habitat quality in remaining habitat patches becomes deteriorated due to the influx of pesticides, pollutants and nitrogen [4–5]. These factors drive species losses and the reduction of species´ abundances, as identified for grassland butterfly species, particularly on calcareous grasslands [6]. This situation becomes even more precarious due to the fact that most remaining habitats are small and geographically isolated and many butterfly species show sedentary dispersal behaviour, with limited migrations across the landscape matrix [7–8]. This may accelerate local extinctions with the result that many potential suitable habitats remain `un-occupied´ [9]. Species demanding specific habitat requirements and with restricted dispersal suffer in particular under habitat fragmentation and the deterioration of habitat quality [10].
In nature conservation, especially high quality habitats considered a priority for habitat management and the protection of sites. Remote sensing techniques and high resolution aerial photographs may help to identify such high value sites. In this study we measure the availability of suitable microhabitats for three butterfly species differing in habitat specialisation and dispersal behaviour [11–20]. To identify these microhabitats we assessed species abundances at points which were randomly set. In parallel, we performed mark-release-recapture techniques to quantify the dispersal behaviour of our targeted species. To identify the spatial composition availability of microhabitats we used high-resolution aerial pictures taken from a UAV equipped by RGB / NIR cameras. Subsequently, Species Distribution Models were run (i) based on species’ presence/absence, and (ii) species abundance of individuals of the respective species reflecting home range preferences. We will interpret our model-outcome against the background of species´ specific dispersal behaviour obtained from our mark-release-recapture data. Based on our data we test the following hypothesis: (i) Suitable habitats for specialist butterflies occur patchy and are mainly restricted to managed sites if compared with habitat generalist butterflies; (ii) Specialist species show more restricted dispersal behaviour if compared with generalist species.
Material and methods
Study species
Our selected butterfly species occur in high population densities across our study area. Melanargia galathea is a common habitat generalist found at flower-rich fresh to dry grassland sites, which uses a comparatively broad ecological niche and occurs in high densities in extensively used grasslands, along road verges and forest edges [11–14]. Caterpillars feed on different species of grasses [11–14]. The species can be classified as a mobile butterfly [13], see also previous MRR studies [15]). The two other butterflies, Erebia medusa and Coenonympha arcania are more specialised than M. galathea. Erebia medusa occurs at nutrient-poor semi-open habitats like extensively used meadows and young fallows, fens, marshy forest meadows and clear cuttings and needs wind protecting structures as forest edges or hedge rows [16]. Its larval food plants are grasses (Poaceae mainly) [12]. The butterfly species requires specific open sites within the swards for oviposition. Eggs are deposited mostly more than 10 cm above ground so that moving before the eclosion of the caterpillars (i.e. before mid- to end-June) will destroy the next generation [16]. Consequently, the species decreases strongly on the landscape level due to agricultural intensification and afforestation [17]. This butterfly shows a sedentary dispersal behaviour [18]. Similarly specialised and sedentary is Coenonympha arcania. This butterfly occurs restricted to structure-rich, nutrient-poor grasslands, especially near forests, clearings in flood plain woodlands, but also in light sunny forests with a lot of understory [13]. Its larvae feed on various grass species. Coenonympha arcania is strongly declining in some regions of Central Europe due to the vanishing of extensively used fresh grasslands and young heterogeneous fallow land (mosaic of woody-open structures) [17, 19–20].
Study area
As study area we selected the managed calcareous grassland site `Dietersheimer Brenne´ (45ha), surrounded by agricultural land, forest, roads and paths. Our study area is located south of the city Freising (48.29N, 11.68E). The `Dietersheimer Brenne´ is characterised by high intensity of solar radiation, high temperature peaks, and dry climatic conditions [21]. Mesobromion is the dominating vegetation formation which consists of a species rich flora and fauna, including many rare and protected taxa [2]. Current conservation management antagonise the reforestation of this grassland site.
Aerial pictures
Aerial pictures were taken with a UAV (DJI Phantom 3 Professional) equipped with a 3-axis gimbal, standard DJI RGB camera (1/2.3” CMOS sensor), and a modified DJI camera for acquisition of Near Infrared Imagery (3.97mm f/2.8 82° HFOV Red + NIR channel) [22]. The flight was conducted in August 2017 at a height of 60 m aboveground using an orthogonal camera attachment with 12MP resolution, 94° field of view (FOV), 20mm focus equivalent and 2s time lapse interval while flying with 5 m/s. Flight path was prepared using the software Litchi for DJI Phantom [23], with stop and turn curve mode settings. Missions were prepared with minimum 50% overlap between consecutive pictures and between neighbouring flight lines to minimize lens distortion effects on results. In total, four flight lines in one flight mission (15 min) were conducted with each flight line of 250 m length composed as rectangular shape. For NIR photographs we used the MAPIR Camera Reflectance Calibration Ground Target Package (V1) [24]. In total 1825 usable pictures were taken for RGB and NIR. Single aerial pictures were post-processed using AgiSoft Photoscan Professional software with high-quality dense cloud processing and mesh construction settings. The calculated Digital Surface Model (DSM) was used in Photoscan Professional to create minimal distorted orthomosaic, due to undulating relief in the study area. NIR aerial imagery was reflectance calibrated using MAPIR Camera Control Application (MCC) and processed with identical settings. Based on geocoded imagery and visual verification of accuracy processed imagery was exported as orthomosaic in geotif raster file format with 0.02 m cell size and resampled to 0.2 m cell size for subsequent use in species distribution modelling. The processed orthomosaic achieved absolute projection error below 3m for XY-errors, due to missing high precision ground references. Aerial pictures were decomposed into RGB + NIR channels, which were used as input variables for the habitat suitability models. In preparation to modelling approach a Digital Surface Model DSM was created, too, using 0.04m resolution and clipped and resampled to same extent as orthomosaic using bilinear interpolation. For both, occupied and unoccupied plant sites we approved all locations again in the field using the previously taken GPS coordinates which are now combined with the generated and georeferenced orthomosaic. As all point localities coincide with our field sites we used all data for further habitat suitability models.
Mark-release-recapture
Mark-release-recapture was conducted inside the managed grassland site `Dietersheimer Brenne´. Points for butterfly capturing were randomly set across the study site (with 65 points in 2014 and 20 points in 2015 and 2017). Point selection was done using “random points in layer bounds” in QGIS 2.18 [25]. At each of these points we captured the targeted butterfly species for 10 min during its flight season. Data were collected during the flight period of the respective study species. Capturing of butterflies took place between 10am and 4pm during days with sunny weather conditions and air temperatures above 20°C. Butterflies were netted in the field and marked with a running number with a water proof pen. For each individual we recorded date, time and sex. Individuals which were already marked were classified as recaptures. For all species we considered the number of individuals marked at each study point, as well as the proportion of migrating individuals (in general, and sex-wise). Furthermore we analysed mean, minimum and maximum distances across all individuals, and for both sexes, separately. Radius of action was calculated using minimum bounding geometries in QGIS 2.18 [25] covering sampling plots were marked individuals were captured within one season. All raw data are provided in S1 and S2 Appendix.
Species distribution models
For SDM development we pooled the observations at the randomly selected points for each species. As the datasets also allow to quantify relative abundances at the sampled plots we followed two different approaches: we (1) run standard SDMs using only one observation per plot if the species was generally present and (2) added each observation to the training dataset for SDM development, i.e. including multiple observations at the same site if present. The latter approach weights the modelled preference of the species in environmental space according to the measured abundance in the field, while the first approach assesses the relative suitability of a given site based on the general distribution of the species in environmental space across occupied sites throughout the study area irrespective of local abundances. The movement information of each species as described by our mark-release-recapture data was used to assess the spatial distribution and accessibility of suitable micro-habitats.
The potential micro-habitats of the study species were quantified based on the remote sensing variables derived from areal pictures taken by the UAV, i.e. the DSM, the red, green and blue band of the aerial photo as well as the normalized differentiated vegetation index (NDVI; calculated as NDVI = (NIR—red)/(NIR + red)). These variables are suitable to describe both the general surface cover and its spatial structure as well as to discriminate different vegetation types. All environmental variables had a pairwise inter-correlation of R2 < 0.75 and can hence be treated as predominately independent variables, as required by the SDM algorithm.
As presence-pseudoabsence algorithm, we chose Maxent 3.4.1. as machine learning technique [26–27]. Maxent follows a maximum entropy approach and derives the potential distribution of a species from environmental conditions at its presence locations and the available environmental space within the study region. Testing various combinations of environmental variables and features thereof, Maxent uses a random walk methodology to derive the best suited model to describe the distribution of species’ occurrences without making too complex assumptions. As environmental background to select pseudo-absences, we defined an area defined by a buffer of 50 m enclosing the sample plots. Restricting the background is necessary to account for dispersal limitations of the target species as pseudo-absences are ideally selected from potentially accessible sites. As the results of single Maxent SDMs due to the machine-learning approach and data splits for model evaluation may provide multiple alternative best fitting models, multiple SDMs were computed. Therefore, for model evaluation across our both approaches we computed 100 single SDMs, each splitting the species records in 80% used for model training and 20% used for model evaluation based on the area under the receiver operating characteristic curve (AUC; [28]). The averages across all single SDM projections using the Cloglog output format were used for further processing and rescaled based on the minimum training presence threshold. Potential extrapolation beyond the training range of the variables was quantified using multivariate environmental similarity surfaces (MESS; [29]), which allow the identification of those sites characterized by environmental conditions beyond the training conditions of the SDMs. These sites are accompanied by a high prediction uncertainty and were omitted in subsequent analyses.
Results
Mark-release-recapture
The average radius of action (means across all years, 2014, 2015 and 2017) varied strongly among the three species analysed (χ2 = 59.4, P < 0.05). We found strongest dispersal for M. galathea (number of recaptured individuals, n = 85, 70.3m (SE = 3.4) (ranging from 18-147m), intermediate dispersal for E. medusa (n = 101) 52.9m (SE = 3.5), ranging from 8–150m), and most restricted dispersal behaviour for C. arcania (n = 100, 34.5m (SE = 3.0), ranging from 5–138m).
Habitat suitability models / Species distribution models
From a statistical point of view both approaches revealed good results in terms of AUC scores (AUCunique, training: range = 0.815–0.858; AUCunique, test: range = 0.785–0.829; AUCabundance, training: range = 0.822–0.900; AUCabundance, test: range = 0.819–0.896; Table 1). Across all models, RGB Red had the highest variable contribution, followed by DSM in most cases (Table 1). RGB Green had in all models the lowest explanative power. The same general pattern of variable contributions was suggested when using a permutation procedure. Within the projection area of the SDMs, the MESS analyses revealed no sites which would require model extrapolation beyond the training range. Hence, MESS results are not shown.
Comparing the results obtained with single observations per grid cell most open areas within the study area are suitable for all three species. However, accounting for local abundances the model suggests a much more patchy occupancy of micro habitats reflecting conditions of homerange cores, wherein M. galathea may occupy a much larger area then E. medusa and C. arcania (Fig 1). The results suggested by the SDMs are in concert with the expectations on the habitat requirements and the availability of species´ specific microhabitat niche.
Warmer colours suggest higher environmental suitability. Circles indicate sampling locations where mark-release-recapture was conducted. Black areas were not covered by UAV flights.
Discussion
Ecological requirements and dispersal
Our models show that the degree of ecological specialization is reflected by the distribution of suitable habitats. Most suitable habitats are identified for Melanargia galathea, while suitable habitats occur more patchy for the two species Erebia medusa and Cynonympha arcania. Our data are congruent with previous studies on butterfly species supporting that increasing habitat specialization lead to the reduction of the availability of suitable habitats, and thus reduces species´ occurrences [30]. Studies showed that such specialized species are only hardly able to respond on environmental changes [31]. Butterflies on the need for specific micro-habitat structures, specific larval food plants and/or specific nectar sources suffer in particular under recent landscape homogenization [13, 32]. In addition, specialist species are characterized by restricted dispersal behavior [13]. This combination of patchy occurrences and restricted movement might aggravate potential negative effects arising from habitat fragmentation and the deterioration of habitat quality. Subsequently, such taxa are assumed to suffer more under demographic and environmental stochasticity than species equipped with higher ecological plasticity [33–35]. Restricted movement and the fact that most high-quality habitat patches are small and geographically isolated from each other causes a silent vanishing of many taxa across our landscape, and thus, many suitable habitats remain unoccupied [9, 36].
Our models identified a high habitat quality along roads and paths for all three model species. This might arise from occasional disturbances which may produce suitable habitat structures for butterflies. These habitats along linear pathways may also serve as suitable corridors for the migration of butterflies, and its exchange among suitable habitats [37–39]. We argue that, particularly habitat specialists with restricted dispersal behavior such as Erebia medusa and Cynonympha arcania, may profit from such surrogate habitats and corridors which may help to increase landscape permeability.
Habitat quality
Our study area consists of a mosaic of managed extensively used grassland (Mesobromion), forests, agricultural fields, roads and paths. Our SDMs show a high habitat quality for all three study species along roads and paths. In contrast, major parts of the managed grassland are of rather low quality, despite land management. This may have various reasons: Atmospheric nitrogen deposition may diminish habitat quality, particularly of nitrogen-limited ecosystems such as calcareous grasslands [40–41]. Previous studies showed that typical grassland butterflies suffer under current atmospheric nitrogen deposition which changes vegetation structures. Atmospheric nitrogen increase growth rates of vegetation and thus may lead to higher standing grasses and denser vegetation; This affects microclimatic conditions (increasing shadow and subsequent higher humidity and lower temperatures), and the outcompete of warm-loving typical calcareous grassland species. This either reduces micro-habitat quality or may completely lead to the vanishing of specific grassland species, with detrimental effects on the successful larval development of butterfly species [6]. This particularly affects butterfly species which rely on typical grassland plant species. This trend of deterioration of habitat quality is currently observable in many nature reserves of the EU, which are embedded in agricultural landscapes [40–41]. The fact that road verges and paths, characterized by occasional disturbances and thus a rather open vegetation might be another evidence, that atmospheric nitrogen causes the deterioration of (undisturbed) grassland habitats. Another potential reason may arise from management strategies itself: previous studies showed that specialized butterfly species may respond highly sensitive on marginal changes in habitat structures [42–44]. Thus, conservation management focusing on other taxa such as orchids may lead to the vanishing of endangered arthropod populations [45–46].
Conclusion
With our data we approve both hypothesis raised at the end of our introduction. The combination of detailed field observations (such as MRR) with high resolution aerial pictures allows the projection of suitable micro-habitats as well as the quantification of micro-habitats used during transitional movement or as home range cores. This approach underlines that classical field observations provide ecological and behavior (as habitat use and dispersal behavior) needed. Results obtained from SDMs can only become interpreted against the background of field observations, such as dispersal behavior. For example, suitable habitat patches might exist across the landscape, however, its (re)colonisation by a butterfly species might be predicted to be unlikely due to its restricted dispersal behavior, according to MRR data. Important synergies may arise from both approaches and allow the combination of two spatial scales (from microhabitat to landscape scale). Advantages of the UAV-acquired imagery compared to very high-resolution (VHR) imagery are the flexibility for the date of acquisition to ensure optimal differentiation between different habitat types. Furthermore, the simultaneously acquired data were necessary for correct SDM calculations. Although the mentioned missing high precision ground reference points reduce absolute positioning accuracy the usage of VHR imagery (i.e. Quickbird, Pleiades) could not provide the needed landscape details and habitat parameters.
Supporting information
S1 Appendix. Mark-release-recapture data of all three butterfly species studied.
Given are individuals captured and recaptured per plot and sex of each individual, as well as exact GPS coordinates of each plot.
https://doi.org/10.1371/journal.pone.0207052.s001
(XLSX)
S2 Appendix. Dispersal data assessed for each of the three butterfly species.
https://doi.org/10.1371/journal.pone.0207052.s002
(DOCX)
Acknowledgments
We thank Edyta Czubernat, Krishna Cholleti, Max Prietzel, Klara Forstner, Johannes Frost, Simone Kainz, Iris Dallmann, Kira Kniffka and Evi Weidner (Freising, Germany) for their help in the field (MRR), and Stefan Reuter (Freising, Germany) for taking aerial pictures with the UAV. We thank Thomas Schmitt and two anonymous referees for helpful comments on a first version of this article.
References
- 1. Dennis RLH, Eales HT (1997) Patch occupancy in Coenonympha tullia (Muller, 1764) (Lepidoptera: Satyrinae): habitat quality matters as much as patch size and isolation. J Insect Conserv 1: 167–176.
- 2. Habel JC, Rödder D, Lens L, Schmitt T (2013) The genetic signature of ecologically different grassland Lepidopterans. Biodiv Conserv 22: 2401–2411.
- 3. Batáry P, Gallé R, Riesch F, Fischer C, Dormann CF, Mußhoff O et al. (2017) The former Iron Curtain still drives biodiversity–profit trade-offs in German agriculture. Nature Ecol Evol 1: 1279–1284.
- 4. Stevens CJ, Dise NB, Mountford JO, Gowing DJ (2004) Impact of nitrogen deposition on the species richness of grasslands. Science 303: 1876–1879. pmid:15031507
- 5. Geiger F, Bengtsson J, Berendse F, Weisser WW, Emmerson M, Morales MBet al. (2010) Persistent negative effects of pesticides on biodiversity and biological control potential on European farmland. Basic Appl Ecol 11: 97–105.
- 6. Habel JC, Segerer A, Ulrich W, Torchyk O, Weisser WW, Schmitt T (2016a). Butterfly community shifts over two centuries. Conserv Biol 30: 754–762.
- 7. Van Dyck H, Baguette M (2005) Dispersal behaviour in fragmented landscapes: Routine or special movements? Basic Appl Ecol 6: 535–545.
- 8. Stevens VM, Turlure C, Baguette M (2010) A meta-analysis of dispersal in butterflies. Biol Rev 85: 625–642. pmid:20055815
- 9. Habel JC, Schmitt T (2018) Vanishing of the common species: Empty habitats and the role of genetic diversity. Biol Conserv 218: 211–216.
- 10. Hughes JB (2000) The scale of resource specialization and the distribution and abundance of lycaenid butterflies. Oecologia 123: 375–383. pmid:28308592
- 11.
Weidemann HJ (1995) Tagfalter: beobachten, bestimmen. 2. Auflage. Weltbild Verlag. Augsburg.
- 12.
Ebert G, Rennwald E (1991) Die Schmetterlinge Baden-Württembergs. Band 1,E. Ulmer, Stuttgart, 522 pp.
- 13.
Bink FA (1992) Ecologische atlas van de dagvlinders van Noordwest-Europa. Haarlem. Schuyt.
- 14.
Asher J, Warren M, Fox R, Haidring P, Jeffcoate G, Jeffcoate S (2001) The Millennium Atlas of Butterflies in Britain and Ireland. Oxford University Press, Oxford.
- 15. Habel JC, Junker M, Schmitt T (2010) High dispersal ability and low genetic differentiation in the widespread butterfly species Melanargia galathea. J Insect Conserv 14: 467–478.
- 16. Schmitt T (2002) The biology of Erebia medusa ([Denis & Schiffermüller],1775) in Central Europe (Lepidoptera). Acta Biol Debr 24: 113–129.
- 17. Reinhardt R, Bolz R (2012) Rote Liste und Gesamtartenliste der Tagfalter (Rhapalocera) (Lepidoptera: Papilionoidea et Hesperioidea) Deutschlands. Naturschutz und Biologische Vielfalt 70: 167–194.
- 18. Schmitt T, Varga Z, Seitz A (2000) Forests as dispersal barriers for Erebia medusa (Nymphalidae, Lepidoptera). Basic Appl Ecol 1: 53–59.
- 19. Schneider C (2003) The influence of spatial scale on quantifying insect dispersal: an analysis of butterfly data. Ecol Entomol 28: 252–256.
- 20. Schneider C, Dover J, Fry GLA (2003) Movement of two grassland butterflies in the same habitat network: the role of adult resources and size of the study area. Ecol Entomol
- 21.
Autengruber H, Denic M, Klein T (2007) Kalkmagerrasenrelikte der Dietersheimer Brenne. Zustand, Gefährdung, Pflege und Nutzungskonflikte. Großpraktikum am Institut für Landschaftspflege und Botanik (TU München), Freising.
- 22.
Peau Production, https://www.peauproductions.com/collections/dji-phantom-3-lenses/products/gp39728 (accessed 19.2.2018)
- 23.
Flylitchi, https://flylitchi.com/ (accessed 19.2.2018)
- 24.
Mapir, https://www.mapir.camera/products/mapir-camera-reflectance-calibration-ground-target-package-v1 (accessed 19.2.2018)
- 25. QGIS (2018) QGIS Geographic Information System. Open Source Geospatial Foundation. URL http://qgis.osgeo.org (accessed 19.2.2018)
- 26. Phillips SJ, Dudík M (2008) Modeling of species distributions with Maxent: new extensions and comprehensive evaluation. Ecography 31: 161–175.
- 27. Elith J, Phillips SJ, Hastie T, Dudik M, Chee YE, Yates CJ (2011) A statistical explanation of MaxEnt for ecologists. Div Distr 17: 43–57.
- 28. Swets JA (1988) Measuring the accuracy of diagnostic systems. Science 240:1285–1293. pmid:3287615
- 29. Elith J, Kearney M, Phillips SJ (2010) The art of modelling range-shifting species. Methods Ecol Evol 1: 330–342.
- 30. Louy D, Habel JC, Schmitt T, Assmann T, Meyer M, Müller P (2007). Strongly diverging population genetic patterns of three skipper species: the role of habitat fragmentation and dispersal ability. Conserv Genet 8: 671–681
- 31. Henle K, Davies KF, Kleyer M, Margules C, Settele J (2004) Predictors of species sensitivity to fragmentation. Biodiv Conserv 13: 207–251.
- 32. Stefanescu C, Carnicer J, Peñuelas J (2011) Determinants of species richness in generalist and specialist Mediterranean butterflies: the negative synergistic forces of climate and habitat change. Biol Conserv 34: 353–363.
- 33. Krauss J, Steffan-Dewenter I, Tscharntke T (2003) How does landscape context contribute to effects of habitat fragmentation on diversity and population density of butterflies? J Biogeogr 30: 889–900.
- 34. Polus E, Vandewoestijne S, Choutt J, Baguette M (2007) Tracking the effects of one century of habitat loss and fragmentation on calcareous grassland butterfly communities. Biodiv Conserv 16: 3423–3436.
- 35. Devictor V, Julliard R, Jiguet F (2008) Distribution of specialist and generalist species along spatial gradients of habitat disturbance and fragmentation. Oikos 117: 507–514.
- 36. León-Cortés J, Lennon JJ, Thomas CD (2003) Ecological dynamics of extinct species in empty habitat networks. 1. The role of habitat pattern and quantity, stochasticity and dispersal. Oikos 102: 449–464.
- 37. Saarinen K, Valtonen A, Jantunen J, Saarino S (2005) Butterflies and diurnal moths along road verges: Does road type affect diversity and abundance? Biol Conserv 123: 403–412.
- 38. Cousins SAO (2006) Plant species richness in midfield islets and road verges–The effect of landscape fragmentation. Biol Conserv 127: 500–509.
- 39. Valtonen A, Saarinen K, Jantunen J (2006) Effect of different mowing regimes on butterflies and diurnal moths on road verges. Anim Biodiv Conserv 29: 133–148.
- 40. Bobbink R, Hornung M, Roelofs JGM (1998) The effects of air-borne nitrogen pollutants on species diversity in natural and semi-natural European vegetation. J Ecol 86: 717–738.
- 41. Lethmate J (2005) Stickstoff-Regen: Ein globales Eutrophierungsexperiment. Biologie in unserer Zeit 35: 118–135.
- 42. Habel JC, Teucher M, Ulrich W, Bauer M, Rödder D (2016b) Drones for butterfly conservation: larval habitat assessment with an unmanned aerial vehicle. Landscape Ecol 31: 2385–2395.
- 43. Filz K, Engler JO, Stoffels J, Weitzel M, Schmitt T (2013) Missing the target? A critical view on butterfly conservation efforts on calcareous grasslands in south-western Germany. Biodiv Conserv 22: 2223–2241.
- 44. Schmitt T (1993) Biotopansprüche von Erebia medusa brigobanna Fruhstorfer, 1917 (Rundaugen-Mohrenfalter) im Nordsaarland. Atalanta 24: 33–56.
- 45. Slámová I, Klecka J, Konvicka M (2013) Woodland and grassland mosaic from a butterfly perspective: habitat preferences of Erebia aethiops (Lepidoptera: Satyridae) Insect Conservation and Diversity 6: 243–254.
- 46. Konvicka M, Benes J, Cizek O, Kopecek F, Konvicka O, Vitaz L (2008) How too much care kills species: grassland reserves, agrienvironmental schemes and extinction of Colias myrmidone (Lepidoptera: Pieridae) from its former stronghold. J Insect Conserv 12: 519–525.