Resource use in two contrasting habitat types raises different challenges for the conservation of the dryad butterfly Minois dryas

The suitability of any location for a given species is determined by the available resources. However, there are many species that occur in more than one habitat type and their successful conservation may be particularly difficult. The dryad Minois dryas, a locally endangered butterfly, occurs in two contrasting habitats—xerothemic and wet grasslands. We investigated the influence of various habitat characteristics, such as vegetation height, grass cover, proximity of shrubs, plant species composition, Ellenberg indices of trophic and microclimatic conditions, on the microhabitat selection by the species. The nectaring of randomly selected butterflies was observed and habitat characteristics were compared at random points within the meadow and at the butterfly’s nectaring and resting places. The butterflies generally preferred to stay close to shrubs and avoided invasive goldenrods. Thermal conditions and the availability of nectar plants were the factors limiting the dryad’s use of wet grassland. In xerothermic habitats grass cover affected the distribution of butterflies. Concerning the availability of larval host plants, wet meadows proved potentially more favourable, whereas nectar resources for adults were more abundant in xerothermic grasslands. Based on our findings, conservation strategies for this butterfly must differ in the two habitats. Rotational mowing in xerothermic grasslands and the removal of invasive goldenrods in wet grasslands are the recommended actions. At a larger spatial scale, a habitat mosaic composed of xerothermic and wet grasslands in close proximity would seem to be the most suitable areas for the conservation of the dryad.


Introduction
Habitat availability has become a key problem in conservation biology (WallisDeVries et al. 2002;Collinge et al. 2003;Dennis et al. 2003;Timonen et al. 2011), because the worldwide problem of habitat fragmentation has led to a loss of continuous habitat, reduction in habitat patch sizes and an increase of patch isolation (Andrén 1994;Fahrig 2003). However, habitat delimitation is a complex task (Hall et al. 1997;Dennis et al. 2003;Vanreusel and van Dyck 2007;Dover and Settele 2009). Recently the concept of resource-based habitats has been introduced (Dennis et al. 2003;Dennis 2010;Turlure et al. 2009), which offers an alternative to the traditional classification of landscape into patches and matrix (Vanreusel and van Dyck 2007). The concept assumes that the landscape is a continuum of overlapping resource distributions. Resource availability may affect patch quality and structure, as well as interresource distances . The resource-based approach may help to recognize and explain the spatial distribution patterns of animals, and to capture the functional and spatial interaction of animals with their environment (Dennis et al. 2003).
The concept of resource-based habitats is particularly important for the conservation of many insects, such as butterflies or bees. For these insects, habitats are determined not only by the presence of foodplants and nectar plants, but also by the microclimate and vegetation structure, which provides crucial components such as resting, roosting, mate-location and oviposition places, as well as shelter from wind and refuge or escape from predators (Sparks and Parish 1995;Shreeve et al. 2001;Konvicka et al. 2002;Hardy and Dennis 2007;Slamova et al. 2011). Thus, for the successful conservation of an endangered species, knowledge of the crucial habitat components as well as the species biology is essential Turlure et al. 2009). However, information about habitat requirements of many endangered invertebrates is still missing and, therefore, the management strategies for their habitats remain unknown (WallisDeVries et al. 2002).
Furthermore, there are many animal species that occur in more than one habitat type (e.g. Bell et al. 2001;Bonte et al. 2002;Mildenstein et al. 2005;Skórka et al. 2006;Hahn et al. 2011). These habitat types are sometimes so distinctive in terms of their resources, vegetation and structure that the behavior and biology of a species may differ substantially depending upon the type (Boudjemadi et al. 1999;Potts et al. 2001;Allen and Thompson 2010). Species conservation, especially for endangered species using differing habitat types, is particularly challenging and requires a differentiated application of methods (Mayberry and Elle 2010). One such species is the dryad Minois dryas (Scopoli, 1763), a grassland butterfly of the family Nymphalidae.
The dryad occurs in two distinct habitat types, namely (1) wet grasslands (ordo Molinietalia) and (2) xerothermic grasslands and steppe-like habitats, dominated by the communities of the class Festuco-Brometea (Krzywicki 1982;Ebert and Rennwald 1991;Dąbrowski 1999;Beneš et al. 2002). Both habitat types are themselves endangered. Cessation of extensive management has led to the expansion of shrubs on both xerothermic and wet grassland, and specifically in our study area also to the invasion of alien goldenrods Solidago spp., causing substantial habitat deterioration Rosin et al. 2012).
The aim of this study was to answer the following questions: (1) are the resources used by the dryad in the two habitat types the same? (2) what are the key factors deciding the choice of microhabitats by the butterflies in the two habitats? and (3) what are the management implications for the conservation of the dryad butterfly?

Study species and area
The dryad has a Euro-Siberian distribution range, extending from the northern part of the Iberian Peninsula, through Europe and Asia to Japan (Dąbrowski 1999;Kudrna 2002). In Europe, it occurs in a number of countries, though usually only locally (Warecki and Sielezniew 2008). Despite the fact that in the past the species was recorded in various places in Poland, at present its occurrence is restricted to the southern part of the country, namely the Kraków region and the Eastern Carpathians (Buszko and Masłowski 2008). It is protected under Polish law and is listed as a critically endangered species in the Polish Red Data Book of Animals (Głowaciński and Nowacki 2004). On a global scale, in terms of the IUCN Red List its conservation status is of less concern (Van Swaay et al. 2010).
The larval food plants include Poaceae grasses occurring in the aforementioned habitats, mainly Molinia caerulea as well as Festuca rubra, Bromus erectus, Calamagrostis epigeios, Lolium perenne, Arrhenatherum elatius and, exceptionally, Carex segdes (Settele et al. 1999;Głowaciński and Nowacki 2004). The dryad flight period is from the end of July to the beginning of September.
The study was carried out in an extensive grassland area in Kraków, southern Poland, located approximately 8 km south-west of the city centre ( Fig. 1, Appendix S1). The dryad occurs there in a mosaic of wet meadows in the Vistula valley and xerothermic calcareous grasslands located on hills. Two meadows of each habitat type were selected, namely, the Skołczanka Reserve (A) and Uroczysko Kowadza (B) in the case of xerothermic grasslands, as well as two wet meadows located on either side of a small wooded hill (Fig. 1, Appendix S1). All grasslands under study are located in the Bielańsko-Tyniecki Landscape Park and are part of the Natura 2000 site ''Dębnicko-Tyniecki Obszar Łąkowy'' (PLH 120065).

Field sampling
The field surveys were carried out in August 2009 and 2010, with the exception of the Uroczysko Kowadza meadow, where the fieldwork only took place in 2010. Observations of butterflies were carried out on sunny and windless days, while the measurements of the microhabitats at random points were taken regardless of the weather.
To obtain data on the selection of resting and nectaring places by the dryad, randomly-spotted individuals were followed and the places where they sat on plants were marked with a bamboo pole. In order to describe the general availability of resources and the vegetation composition, 31 and 54 random points were selected respectively at wet and xerothermic meadows through random generation of geographic coordinates. The vegetation within a 1.5-m radius was described for both the random points and the places utilized by the butterflies. Five measurements of the vegetation height (naturally arranged plants), one at the central point and the remaining four 1.5 m apart to the north, south, west and east, were taken with a measuring tape of 1-cm precision, and their mean was calculated. In addition, the distance to the nearest shrub was measured with a precision of 10 cm. Finally, all plant species within a 1.5-m radius were counted and their ground cover was recorded. The cover was measured at the following scale: 1 B 10 %; 2 = 10-20 %, …, 9 = 80-90 % and 10 = 90-100 %. We have also recorded the percentage of bare ground (existent only in xerotermic meadows), but it was not included in the analyses as it was strongly negatively correlated with vegetation height and generally very low and little variable, hardly ever exceeding 10 %.
In each habitat type approximately 40 butterflies were followed. For each individual we noted its sex, and classified its behavior on plants as either resting or nectaring (other types of behavior, e.g. oviposition or copulation, were rarely observed and consequently excluded from the analysis). In order to increase sample size of nectaring observations and to receive the knowledge of the range of plant utilized by the dryads, the behavior of randomlyspotted individuals was also recorded. In such cases the butterflies were followed until they settled on a plant.
When the individual started nectaring, the duration of the behavior was measured and plant species was noted.

Statistical analysis
In a preliminary analysis we compared the characteristics of the surveyed meadows, based on the parameters measured in random points. The multivariate ANOVA was applied with the two habitat types as the main groups and meadow ID nested within them (Appendix S2). The analyses revealed sharp differences between wet and xerotermic meadows, which is obvious taking into consideration different plant species composition in both habitats (cf. Appendix S3). In contrast, there were no differences between particular meadows of each type, apart from those in vegetation height and some of the Ellenberg indicator values (Ellenberg 1974) between the two xerotermic meadows. Nevertheless, since the above differences were relatively small as compared with those between habitats, in the subsequent analyses we tested the patterns of habitat selection by the dryad separately for each habitat type, but not for each meadow.
A detrended correspondence analysis DCA (Ter Braak and Prentice 1988) was performed in order to ascertain whether the butterflies occur more often in places with particular plant species composition. Plant species of both habitats were classified according to following groups: recorded nectar sources, potential nectar sources, grass species comprising potential larval host plant, Apiaceae, and other species, among which special attention was given to the invasive goldenrods (see Appendix S3). The analysis was conducted on the basis of data on the percentage covers of the plant species in the above groups. Plant species communities were compared for the random points, as well as for butterfly nectaring and resting places.
We also investigated microhabitat selection by the dryad, comparing the characteristics of butterfly observation points and random points within their habitats. The analyses were performed using multivariate ANOVA with random points and butterfly observations as the main groups, as well as with nectaring and resting behavior nested within the latter group. The tested microhabitat parameters included distance from shrubs, vegetation height, grass cover, and Ellenberg indicators of light, temperature, humidity and trophic conditions (Ellenberg 1974), determined on the basis of plant species composition using Zarzycki's scale, which is the most appropriate for Poland (Zarzycki et al. 2002). Distances from shrubs were ln-transformed to achieve normality.
In addition, we analysed the dryad preferences for different nectar sources. Nectar plants were classified into the following groups: composite plants (Asteraceae and Dipsacaceae), plants with radial symmetry flowers, plants with dorsiventral symmetry flowers and Apiaceae to analyze nectaring preferences of butterflies (see Appendix S3). The availability of the above groups was then calculated separately for each habitat type, using the entire cover of all the plants in a given group, pooled across all the random points within the habitat type. The proportions of nectaring observations in particular plant groups were compared with their availability, using Bailey's use-availability tests (Bailey 1980). The same procedure was applied to test the dryad's preferences towards invasive alien goldenrod and native nectar plants within the composite plants.
A general linear mixed model (GLMM) was applied to investigate the effects of habitat type, flower type, butterfly sex, as well as their interactions and random effect in the form of the meadow number on the ln-transformed nectaring duration. The Apiaceae plants were omitted from the analysis due to the very small number of nectaring cases in this group. Since the inclusion of random factors did not improve the fit of the model, the final one was reduced to a generalized linear model (GLM) comprising only significant factors.
The detrended correspondence analysis DCA was performed in Canoco for Windows 4.52 (Lepš and Š milauer 2003), while the GLM was done in JMP 9 (SAS Institute 2010). The remaining analyses as well as assumption testing for all the analyses were performed in Statistica 9.0 and MS Excel 2007.

Selection of microhabitats
The detrended correspondence analysis biplot in xerothermic habitat (Fig. 2) shows that plant species composition in the dryad resting sites did not differ from that of the random points. However, plant species communities in both the nectaring and resting sites were subsets of the plant communities at the random points, as indicated by the considerably smaller convex polygons. For wet meadows, plant species compositions in the dryad resting and nectaring sites were similar, but, in turn, they only partly overlapped with the plant species composition at the random points (Fig. 2). In wet meadows, several grass species with higher cover were observed, while in the xerothermic habitats there were fewer species of grass and their cover was relatively poor (see Appendices S1 and S3).
The analyses of microhabitat conditions based on Ellenberg's indicators showed that, in the wet meadows, there were no differences in either the quantities of light, humidity or the trophic conditions between the resting and nectaring places as well as in these sites versus random points (Table 1). However, in the case of thermal conditions, the butterflies selected places that were significantly warmer than random points (Table 1). On the other hand, in the xerothermic meadow, the thermal, light, humidity and trophic conditions had no effect on the occurrence of butterflies (Table 1).

The influence of shrubs, vegetation height and grass cover on habitat use
In both habitats butterflies selected sites that were closer to shrubs than random points (Table 2; Fig. 3). In the xerothermic grasslands the dryad also appeared to prefer higher vegetation; the test statistic was close to the significance level (Table 2). In addition, although we have not detected a difference in vegetation height between random points and the dryad observation points, among the latter resting places were characterized by significantly higher grass cover than nectaring places (Table 2, Fig. 3). The ANOVA results for wet meadows revealed no significant differences in vegetation height and grass cover for nectaring, resting, and random points (Table 2).

Nectaring preferences
The Bailey's test results showed that, in the xerothermic habitat the dryad more frequently foraged on plants with dorsiventral symmetry flowers than on composite plants and those with radial symmetry flowers. In wet meadows, in turn, plants with radial symmetry flowers were avoided, and the dryad preferred composite plants. Thus, the butterflies foraged relatively more frequently on composite plants than those with dorsiventral symmetry flowers (Fig. 4). In both habitats, the dryad avoided Apiaceae (Fig. 4). Native species were selected while goldenrods were avoided (Fig. 5).
The dryad nectaring time was longer in wet meadow than in the xerothermic habitats (Fig. 6). The butterflies fed significantly longer on composite plants (Tukey post hoc, P = 0.001) and plants with dorsiventral symmetry flowers (Tukey post hoc, P = 0.002) than on those with radial symmetry flowers (Fig. 6). Females fed longer than males (Fig. 6). Nectaring time decreased with the progression of the season (GLM: F 2,420 = 16.905; P \ 0.001; slope = -0.048), (Fig. 6). The GLM showed no significant interactions between habitat type, flower type and sex.

Discussion
The results of this study demonstrate a differential use of resources and microhabitats by the dryad in wet meadows The dryad behavior, i.e. resting or nectaring, was nested within the latter group and xerothermic habitats; however, some similarities in the utilization of these two contrasting habitats have also been found. In wet meadows, the factor limiting the dryad's use of microhabitat seems to be the thermal conditions. In the xerotherimic habitat, on the other hand, none of the factors examined affected the distribution of butterflies during resting behavior; however, nectaring butterflies chose slightly drier places, as indicated by the Ellenberg indicators.
It can be concluded that the dryad is a stenothermic species and its distribution may be restricted by unsuitable thermal conditions within its potential habitats. Xerothermic grasslands are usually sunny and warm, as evidenced by the presence of such plant species as Seseli annuum and Vincetoxicum hirundinaria (Zarzycki et al. 2002); thus, in terms of thermal conditions, they are beneficial to the dryad. Nectaring places in the xerothermic habitat were generally open, with low vegetation. At these sites, limestone rocks often protruded above the ground. However, such places may be too dry or too hot for the dryad to carry out all its daily activities there and hence it chooses more humid places with higher vegetation and with higher grass cover for resting.
The results of the study have shown that shrubs are an important component of the dryad habitat in both types of meadow. The presence of shrubs in these meadows is a result of meadow abandonment and the resulting natural succession Rosin et al. 2012). Shrubs provide refuges for butterflies and safe places for mating as copulations were observed there (authors' unpublished data). Shrubs may also provide shelter for butterflies during windy and rainy conditions (Sparks and Parish 1995;Dover Fig. 3 Mean (with 95 % confidence intervals) distance to shrubs, vegetation height, and grass cover at random points (white bars) as well as in the dryad nectaring and resting places (grey bars) in both investigated habitat types: a, c, e xerotermic grasslands; b, d, f wet meadows Dover and Sparks 2000;Dennis and Sparks 2006).
The composition of plant species communities did not affect the distribution of the dryad in xerothermic habitats. Nectaring plant resources were widespread there and thus they apparently did not limit the occurrence of butterflies. The same refers to grasses, i.e. potential larval food plants, despite the fact that they were less numerous in xerothermic grasslands. In contrast, in the wet meadows, the butterflies were concentrated predominantly in places where nectar plants were abundant. In the wet meadows there were already less flowering plants in August than there were in the xerothermic grasslands and this may have affected the dryad's flight period. In both 2009 and 2010, the number of flying butterflies fell sharply in the wet meadows after 14th August, while they could be observed flying in substantial numbers till the end of that month in the xerothermic grasslands (authors' unpublished data). The dryad preferred composite flowers that offer spatially concentrated nectar resources. The longer nectaring time in the wet meadows may imply a lower availability of nectar, especially towards the end of the flight period. This longer nectaring time is in line with the predictions of the marginal value theorem (Charnov 1976;Parker and Stuart 1976;Parker 1992). If the resources, which, in our case are nectaring plants, are scarce, then individuals should opt for the maximal utilization of the resources available at one plant, rather than searching for a new one, because the energy losses during the search may exceed the potential profit to be obtained from new plants. The avoidance of flowers with radial symmetry in the wet meadows is also in agreement with the above explanation. These flowers may have less nectar, which is particularly true for the pollen flowers of Hypericum spp. (Kugler 1955).

Conservation recommendations
Our findings allow a conclusion that a xerothermic habitat is better for the dryad in terms of nectar resources. In turn, a wet meadow habitat is likely to be better in terms of larval food-plant availability, owing to greater abundance of grasses. The conservation strategies for this endangered butterfly must therefore differ for these two habitats. In wet habitats, it would be profitable to increase the cover of favorable nectar sources by proper mowing regimes.
In xerothermic meadows, vegetation height should be diverse to provide both favorable resting and nectaring sites. Thus, it may be good practice to apply rotational grazing or mowing (Dolek and Geyer 1997;WallisDeVries et al. 2002;Cremene et al. 2005), since patches under a diverse mowing regime are preferred by butterflies and other insects (Cizek et al. 2012). A single fragment should be mown every 2-3 years. Mowing should be conducted either in the spring, at the end of May, or in the fall, in mid-September at the earliest (Johst et al. 2006;Grill et al. 2008;Wynhoff et al. 2011). In wet meadows, however, mowing should only take place in the late autumn, because they are also an important habitat for endangered Maculinea (= Phenagris) butterflies (Nowicki et al. 2005(Nowicki et al. , 2007. Mowing is necessary not only to sustain the nectar and host plants favorable to the butterflies, but also to prevent the invasion of some notorious alien plants, such as goldenrods Solidago sp. These alien invasive species are a serious threat for biodiversity, predominantly in wet meadows Powell et al. 2011;Vilà et al. 2011) and the dryad did not use this invasive plant as a nectar source. Alien plants can outcompete native plants such as, for example, grasses (Callaway and Aschehoug 2000) leading to deterioration of the habitat quality. Invasive plants may be actively removed, preferably in the early summer, before they flower, and in the autumn (Güsewell 2003).
In many grassland conservation programmes, shrubs are indicators of habitat deterioration and thus their removal is postulated (Erhardt 1985;Krauss et al. 2004;Lenda and Skórka 2010). In the dryad, however, this is not the case. Shrubs appeared to be a key component of the dryad habitat in both the wet and xerothermic grasslands. Therefore, shrubs should be maintained in these meadows, even when grazing and mowing is applied. However, they should be scarce and scattered within a patch, in order not to overgrow flowering plants.
On a larger spatial scale, the sustaining or creation of a mosaic of both wet and xerothermic meadows located close to one another may be a favorable conservation strategy for the dryad. Obviously, it can only be applied in areas where the topographic conditions allow such a mosaic. Butterflies are known to be able to move between different habitat types in order to seek the appropriate resources (Schwarzwälder et al. 1997). In the case of the dryad, the occurrence of the both wet and xerothermic meadows in close proximity may be also advantageous because local populations of the dryad could utilize both habitat types differentially, depending on annual weather patterns and management. An arrangement of both habitats such as this should increase the viability of the entire metapopulation by reducing the impacts of natural catastrophes such as fires on xerothermic grasslands and flooding in wet meadows. Finally, a mosaic of wet and xerothermic meadows may be the most suitable areas for the reintroduction of the dryad to its former localities.