Skip to main content

Dispersal and adaptation strategies of the high mountain butterfly Boloria pales in the Romanian Carpathians



Habitat quality is one main trigger for the persistence of butterflies. The effects of the influencing biotic and abiotic factors may be enhanced by the challenging conditions in high-alpine environments. To better our knowledge in this field, we performed a mark-release-recapture study with Boloria pales in the Southern Carpathians.


We analysed population structure, movement and foraging behaviour to investigate special adaptations to the alpine environment and to reveal differences between sexes. We compared these aspects in one sector with and one sector without grazing to address the effects of grazing intensity on habitat quality.


We observed “soft” protandry, in which only a small number of males appeared before females, and an extended emergence of individuals over the observed flight period, dividing the population’s age structure into three phases; both observations are considered adaptations to high mountain environments. Although both sexes were mostly sedentary, movement differences between them were obvious. Males flew larger distances than females and were more flight-active. This might explain the dimorphism in foraging behaviour: males preferred nectar sources of Asteraceae, females Caprifoliaceae. Transition from the grazed to the ungrazed sector was only observed for males and not for females, but the population density was higher and the flight distances of the individuals were significantly longer on the grazed sector compared with the ungrazed one.


Soft protandry, an extended emergence of the individuals and an adapted behavioural dimorphism between sexes render to represent a good adaptation of B. pales to the harsh environmental conditions of high mountain ecosystems. However, land-use intensity apparently has severe influence on population densities and movement behaviour. To protect B. pales and other high-alpine species from the negative consequences of overgrazing, areas without or just light grazing are needed.


The quality of habitats is of great importance for the adaptation of populations and hence is one major driver for the evolution of species [1]. Moreover, persistence of a species in a particular habitat is primarily depending from the ability to adapt to the respective living conditions [2]. However, comparatively little is known about adaptations to environmentally harsh environments. This lack of studies partly results from the complicated study conditions under such hostile environmental circumstances. Consequently, the main focus in autecological research has been on species whose ecological demands and specialisations are linked to conditions in temperate zones. High mountain systems are one important ecosystem type with extreme environmental conditions. These areas exhibit harsh environmental conditions, such as cold snaps that may occur at any time of year [3,4,5]. Furthermore, climate change and direct human activities aggravate the negative effects of these challenging environments on species’ living conditions. To survive, high mountain species consequently have to be adapted to these unpredictable environmental circumstances [6, 7]. Consequently, habitats like high-alpine areas offer excellent conditions to study more general biogeographical and evolutionary questions.

Specific adaptions of butterflies to harsh climatic conditions have been unravelled in several studies [8,9,10,11]. However, little is known about local adaptions and persistence of metapopulations of high-alpine Lepidoptera, especially for species of high-alpine systems other than the European Alps and the North American Rocky Mountains [12,13,14]. Additionally, the effects of land-use have recently become of increasing interest in studies concerning the impact assessment of land-use changes on the ecology of high-alpine regions and species [15,16,17]. Thus, Tappeiner and colleagues [18] predicted negative consequences resulting from land-use intensification, especially for species of formerly extensively used alpine habitats. The negative effects of changes in land-use and regional climate have the potential to particularly impact species with limited dispersal, with habitat changes in particular reducing the persistence probabilities of metapopulations [19].

Therefore, we conducted a field study in the Southern Carpathians to investigate general and local adaptations of the shepherd’s fritillary Boloria pales, a butterfly species with a disjunct distribution in the alpine belt of several European high mountain systems [20,21,22]. In general, Romania’s butterfly diversity is linked to the relatively large number of natural and semi-natural habitat types, including a variety of different grassland types, which are a heritage of the Romanian traditional way of living [23]. These grasslands are maintained by extensive grazing, manual mowing and other traditional agricultural techniques, often without the application of agrochemicals [24]. Changes in these traditional agricultural practices have already been identified and predicted as one major threat to Romania’s butterfly diversity [23, 24]. At high altitudes in the Carpathians, which are mostly dominated by natural grasslands, the main threat is the constantly increasing number of livestock, resulting in a continuously increasing rate of overgrazing.

In Romania, B. pales is a rare species restricted to some few locations of natural high altitude grasslands in the southern Carpathians [25]. Even at its possibly most important stronghold in the Bucegi Mountains, changes in land-use are apparent, especially changing grazing regimes [26, 27]. We therefore aim to reveal the demographic structure and senescence within one metapopulation of B. pales over one entire flight period. We intend to understand whether the species’ regional adaptation to these high mountain habitats is sufficient to maintain a protandrous population structure with all its advantages [28,29,30] or whether a close synchronization of the emergence of both sexes is necessary for surviving the high mountain challenge, as for example in the case of Euphydryas aurinia glaciegenita [31, 32]. We further analyse whether emergence is restricted to a short time window as in many univoltine lowland butterfly species [33,34,35,36,37,38], or is extended over a longer period to spread the risk of encountering highly unfavourable conditions caused by the unpredictable weather in high mountain ecosystems. The latter response is known for populations of this [39] and other species [40] assayed in the Alps.

While male butterflies are often assumed to be the more flight-active sex, for example due to patrolling for females [35, 37, 41,42,43], we test for differences in dispersal behaviour and movement patterns between sexes, by analysing the distances moved and the mean recruitment of flying individuals in the population for both sexes. It is already known that some male and female butterflies prefer nectar containing different ingredients [44, 45]. For example, males of Polyommatus bellargus prefer high amounts of sucrose and total sugar, while females opt for glucose and amino-acids [46]. Therefore, we test whether each sex exhibits preferences for different nectar sources and which flower traits may affect their nectar foraging strategy. Finally, we analyse whether the individuals are specialised on specific flowers or are mostly unspecific, as postulated for B. pales as a species [47].

Independent from the sex-dependent needs, habitat quality in general has an important impact on the behaviour and movement patterns of animals [48]. However, habitat quality for many butterfly species is strongly influenced by grazing intensity and general management. As strong changes in the use of high-altitude grasslands are apparent in many high mountain areas, we examine the effects of grazing intensity on behavioural traits and movement patterns of B. pales by studying patches within a metapopulation with different degrees of disturbance by sheep and cow pasturing. In this context, we also analyse the differential response of both sexes to habitat quality in combination with population density.

To tackle these general aspects, we examined the population ecology of B. pales in the Southern Carpathians from 02 July to 02 September 2014, thus covering most of the flight period, using mark-release-recapture. We selected one of the largest metapopulations known in Romania, in the Bucegi Mountains, a plateau-like region traditionally used as high-altitude sheep pasture, where grazing pressure has strongly increased during the past few decades. These data in combination with results obtained from the eastern Alps [39] allow the analysis of the adaptation of this species to high mountain conditions, the assessment of the different behavioural traits of both sexes and a better understanding of their needs, as well as deepening our understanding of the consequences of the increasing human impact in some parts of the European high mountain regions. We conclude with some general conclusions for nature conservation in high mountain ecosystems.

Materials and methods

Study species

The alpine butterfly Boloria pales (Denis & Schiffermüller, 1775) belongs to the family Nymphalidae. Except for the northern European mountains (Scandes, Scottish Highlands), the species occurs in most high mountain systems of Europe from 1500 to 2700 (maximum 3100) m asl.; in the Pyrenees, populations are known as low as 1100 m asl. [20,21,22, 49]. The species lives at flower-rich high-altitude grasslands, where adults feed at blossoms of the respectively dominant flowering plants [20]. The larvae feed on Viola, Plantago or Polygonum species [21, 49,50,51]. B. pales usually occurs in one generation per year [20, 21, 49, 50] and the flight period is comparatively long, often more than eight weeks; mostly from late June or early July to early September [21, 49, 51]. The larvae hibernate as caterpillars (of the 1st instar) and grow after the snow has melted in springtime [20]. Due to the often high abundance at their flight places, the species is a suitable model organism for analysing the population dynamics of high mountain species.

Study area

The mark-release-recapture (MRR) study area was located in the Bucegi Natural Park / Parcul Natural Bucegi (N 45°21′, E 25°30′) near the hiking trail Cota 2000 in the Southern Carpathians (Sinaia, Prahova county, Romania). The area studied had no natural boundaries and a mostly homogeneous vegetation cover. It was divided into two sectors because of the different grazing intensity (distance between sectors ca. 150 m; altitudinal difference: ca. 30 m). The lower sector of 14.2 ha reached from 1660 to 1980 m asl., the upper sector of 2.2 ha reached from 2010 to 2090 m asl. Both sectors were characterized by alpine and subalpine grasslands (partly dominated by Carex curvula) with generally little inclination, but interspersed with steep, partly rocky slopes and some few cliff edges (Fig. 1). At the lower sector, grazing by a herd of cows (about 80 individuals) and two herds of sheep (about 200 and 350 individuals) was observed each sampling day, whereas the upper sector was not used as pasture, probably because of the proximity to the Cota 2000 trail and other nearby hiking trails. Consequently, although similar in their floristic composition, the abundance of flower heads was considerably higher in the upper sector if compared with the lower areas. In the following, we call the lower sector grazed sector and the upper sector ungrazed sector.

Fig. 1

Map showing the study area in the Bucegi Natural Park with sector boundaries and habitat types

Mark-release-recapture (MRR)

We performed a MRR study from 02 July to 02 September 2014, i.e. almost the entire flight period of B. pales. Weather conditions during our time of study were not unusual for the Bucegi Mountains. However, days with rainfall were frequent (21 of 63 days) with only two longer periods without precipitation, i.e. 31 July to 7 August and 25 August to 2 September. On the other hand, rainy periods were never longer than three continuous days. The highest temperature was reached at 14 August with about 20 °C, whereas the coldest temperature was recorded in the night 28 to 29 August with about − 1 °C.

On each day with suitable weather conditions (weak or moderate wind; no rain or snow), we sampled the entire study area, netting all available individuals. On every sampling day, we started at a different point in the area, to ensure that sufficient information was collected from all parts of the study area for every time of the day. The grazed sector was examined over the total study period, whereas the ungrazed sector was only studied from 13 August to 02 September. The great initial expenditure of time to examine the grazed sector did not permit the analysis of the additional ungrazed sector before.

Every captured butterfly was marked with a durable fineliner using an individual code consisting of a letter for the sampling day and a running number, always starting with one every day. Prior to release, we recorded the following information in addition to the individual code: sex, GPS position of capture point (with a maximal deviation of three metres), capture time, wing wear (scored from 1 for fresh to 4 for heavily damaged [42, 52, 53]), behaviour (i.e. flying, resting, feeding, interaction) and, if applicable, the nectar plant. For all recaptures, made at least one day after the last capture event to avoid capture-release trauma [54, 55], we recorded identical information. With the information on wing conditions, we classified the age of the butterfly individuals [54, 56, 57]. Based on all captured individuals and all capture events, we calculated the age structure along the observed flight period (values based on less than five individuals per sex and day were excluded). An analysis correlating the time elapsed from capture to the next recapture with the average wing decay over this time period allowed the assessment of the sex-dependent wing decay rate per time unit. With the four behavioural categories, we performed χ2 homogeneity tests to detect differences in the behaviour between sexes (to avoid the bias of individual preferences, the behaviour of each individual was included only once at its first capture event).

Furthermore, we performed five plant inventories with an adjusted method following the Braun-Blanquet system [58] to assess plant species and their abundance [59]. We also estimated the cover of all plant species with the Braun-Blanquet scale, covering the entire study area and not just delimited areas. With these data, we calculated the proportion of every plant species in relation to the total plant cover during the study period. Assuming a non-random use of nectar sources by the butterflies, we compared the observed use of nectar sources with the one expected if flower selection was at random using χ2 tests. This test was performed separately for both sexes. With the Jacobs’ index of selection [60] (Jacobs interpretational classification of Jacobs` index for our study: 1 to 0.33 preference, 0.33 to − 0.33 neutrality, − 0.33 to − 1 avoidance) and Bailey’s confidence intervals at p values < 0.05 [61, 62], we tested for differences in preferred nectar sources at the plant families and plant genera level. All univariate statistics were calculated in SPSS 22.0 [63].

Population demography

In the program package MARK 8.0, we used the module POPAN 5.0 [64] to estimate daily population sizes based on the Jolly-Seber model (for open populations) [65]. We estimated population sizes separately for the sexes and sectors. POPAN 5.0 estimates three primary parameters: survival probability (phi), capture probability (p) and proportional recruitment (pent). These parameters may be constant (.), dependent on sex (g), respond to time in factorial (t), linear (T) or quadratic (T2) manner, or display additive (g + t, g + T, g + T2) or interactive (g x t, g x T, g x T2) effects [33, 66, 67]. The capture probability might also depend on the sampling effort (hours) [64]. Sampling effort was calculated as the product of time spent conducting field work (ranging from one to six hours per day) and the number of persons involved (ranging from two to four) (see Additional file 1: Table S1). With a second model in the program MARK (Multi-state Recaptures only [64]), we calculated the transition between both sectors (psi) in the period from 13 August to 02 September [36, 68,69,70,71]. Beside transition between sectors, calculations for survival (S) and capture probability (p) within both sectors were performed separately for both sexes with this module. As in the first module, these parameters may be constant (.), different among sectors (o), time (t) or linear (T) or quadratic (T2) time and can also occur as additive (o + t, o + T, o + T2) or interactive (o x t, o x T, o x T2) effects. The capture probability might also depend on the sector size (size) [64]. After separate Goodness-of-Fit tests for both modules (option: RELEASE), we analysed different combinations of the parameters mentioned above and selected separately for each module the model with the lowest value for the corrected Akaike Information Criterion (AICC) [72, 73] and the smallest numbers of parameters [74] as best supported. With the computed values for Net Birth rates and for rates of proportional recruitment (pent), we completed our classification of the age of the butterflies, derived from wing conditions. With the estimated survival probability (phi), we computed longevity using the formula of Cook [75]: Longevity = − 1 / loge (phi). Additionally, we performed Mann-Whitney U tests to examine significant differences in phi between sexes (both sectors pooled and each sector separately), sectors and days with fine (i.e. sun or loose cloudiness) against days with weather conditions unsuitable for butterfly activity (i.e. with a maximal day temperature below 5 °C or with rain or snow events). Further Mann-Whitney U tests were used to examine differences in the transition between both sectors (psi from the chosen models). This analysis was performed for both sexes pooled and separately. For tests between sectors, we only used data collected after 13 August, when both sectors were analysed simultaneously. We applied Bonferroni correction for all Mann-Whitney U tests because of the increasing risk of type I errors when making multiple statistical tests based on the same data set.

Table 1 Comparison of the best models of the POPAN 5.0 analyses for estimating the daily population sizes of B. pales in 2014 in Parcul Natural Bucegi, Romania: Akaike information criterion (AICC) and number of considered parameters

Mobility and movement patterns

The GPS data set of capture/recapture events obtained in this MRR study was used to analyse the movement behaviour of B. pales. With ArcGIS 10.2.1 (ESRI), we measured the direct distance travelled from the first capture event to the second capture event for all individuals. Based on these data, we tested for differences in mobility between sexes with a Mann-Whitney U test (with Bonferroni-correction). With these single distances, we calculated the inverse cumulative proportion of individuals moving in certain distance classes, with each class representing an interval of 20 m. We fitted these data against two different mathematical models, which are commonly applied to find the best prediction of long-distance movements [31, 33, 37, 42, 76, 77]: the negative exponential function (NEF) and the inverse power function (IPF). To exclude artefacts based on the 20 m interval size, we performed similar analyses with 30 m and 50 m intervals. All calculations were performed separately for both sexes. For the NEF, the relative proportion of individuals moving to distance D is.

INEF = ae-kD respective ln I = ln a - kD.

The parameter a represents a scaling constant while k is the dispersal constant describing the shape of the exponential curve. Under the IPF, the proportion I is expressed as.

IIPF = cD-n respective ln I = ln c - n (ln D),

where c is a scaling constant and n a variable describing the effect of the distance on dispersal [33, 76]. With the calculated coefficient of determination R2 and adjusted coefficient of determination Radj2 for the computed curves, we chose the best model and interval size to predict the proportion of individuals moving greater distances than those covered by our MRR study [31, 33, 42, 76].


During 36 field days, we marked 633 individuals (males: 401; females: 232) with a total of 1202 capture events on both sectors (Additional file 2). The recapture ratio was 51.4% for males and40.1% for females. In total, more than half of the capture events were at the ungrazed sector (622, i.e. 51.7%) during the 11 sampling days there. We obtained multiple recapture events (up to six recaptures for both sexes). The longest residence time between capture and last recapture was 24 days for males and 17 days for females, whereas the calculated longevity was 4.5 ± 0.5 SD days for males and 3.3 ± 0.6 SD days for females.

For comparing the estimated survival probabilities (phi), we used Mann-Whitney U tests because none of the data was normally distributed and all samples from the POPAN calculation were independent from each other. The comparison of the sexes pooled for both sectors revealed significant differences (critical p value with Bonferroni correction: 0.05 / 4 = 0.0125) between the daily survival rates of males and females (males: 0.80 ± 0.14 SD; females: 0.74 ± 0.17 SD; U test: p = 0.002). These significant differences were also found on the grazed sector (males: 0.80 ± 0.14 SD; females: 0.73 ± 0.16 SD; U test: p < 0.001), but not on the ungrazed sector (males: 0.81 ± 0.15 SD; females: 0.77 ± 0.19 SD; U test: p = 0.649). No significant differences between the sectors were found, neither for males (U test: p = 0.893) nor females (U test: p = 0.704). Significantly higher survival probabilities during days with suitable against days with harsh weather conditions were observed for both sectors and sexes (U tests: grazed sector: 0.85 ± 0.03 vs. 0.59 ± 0.17, p < 0.001; ungrazed sector: 0.86 ± 0.13 vs. 0.61 ± 0.11, p = 0.001; males: 0.88 ± 0.06 vs. 0.65 ± 0.16, p < 0.001; females: 0.83 ± 0.07 vs. 0.55 ± 0.15, p < 0.001; all values are means ± SD).

Demography – Grazed sector

The most plausible model estimated a total population size of 550 males (± 45 SE) and 242 females (± 25 SE), i.e. a population density of 39 males (± 3 SE) and 17 females (± 2 SE) per hectare (Table 1). The estimated gender ratio was 2.3 males per female. The phenology based on the POPAN model (Fig. 1a) was protandrous. Thus, the number of males started to increase from the beginning of our sampling, whereas females only noticeably increased some three to four weeks later. With the beginning of the females’ emergence, the phenologies of both sexes ran more or less in parallel. Both sexes increased slowly until late July, then rose more quickly to a peak at 20 August. Thereafter, the male population decreased, whereas the female population remained at a mostly constant level. At the very end of the study period, both sexes showed a slight increase. The sensitivity to bad weather conditions (rain or wind events) was moderate for both sexes at the beginning of the study. In August, when strong winds occurred in combination with strong rainfall events, the effects of weather were greater, but more obvious in males than females (Fig. 2a).

Fig. 2

Estimated population size of B. pales for every sampling day. a Grazed Sector, (b) Ungrazed Sector; Vertical bars represent days with harsh weather conditions (grey bars: no sampling on days with rain events, green bar: day after frost event at night); error bars represent the standard error of the calculated population size from the program MARK; note that the graphs have different units on the x-axis

Demography – Ungrazed sector

The most plausible model with the second lowest value for the AICC and the smallest numbers of parameters estimated 271 males (± 13 SE) and 270 females (± 31 SE), i.e. 123 males (± 6 SE) and 123 females (± 14 SE) per hectare (Table 1). The estimated gender ratio was practically equal (one male per female). Although the ungrazed sector was studied for a shorter period of time, the estimated population sizes of both sectors revealed a by far higher density on the ungrazed sector. The ungrazed sector’s phenology had a higher male population size until 20 August, when both sexes reached their peak at a nearly identical population size. After this peak, the female population stayed constant until 25 August and exceeded the male population (which declined slightly) for the rest of the study period. Thereafter, both sexes decreased in parallel and slightly increased again at the end of August (Fig. 2b). The sensitivity to bad weather conditions was similar to the results on the grazed sector.

Wing conditions

The deterioration rate of wing condition was assessed by comparing the condition of captured and first recaptured individuals. For both sectors combined and each sector separately, both sexes showed a strong correlation between wing decay and the time elapsed between capture events (Pearson correlations, all p < 0.001). The average time needed for the deterioration of wing condition by one category was shorter for males than for females (males: 13.5 days ±1.3 SE; females: 15.6 ± 2.4 SE).

Wing conditions of B. pales were good to very good for most of the individuals throughout the observed flight period. Even at the end of the study period (after mid-August), the wings of most captured individuals were in good condition (i.e. about 75% category 1). Seriously damaged wings (i.e. categories 3 and 4) were exceptions throughout, and freshly emerged individuals (verified by the excretion of a red liquid from the abdomen) of both sexes were always recorded, even at the end of the study period. Thus, old butterflies must have been constantly replaced by freshly emerged ones until mid-August or even later in the flight period.

At the beginning of the studied flight period, average wing condition deteriorated gradually until 03 August (Pearson Correlation: p = 0.004, R2 = 0.59, Radj2 = 0,54). Afterwards, this remained at the same average level until 25August (Pearson Correlation: p = 0.734, R2 = 0.01, Radj2 = 0,00), when the start of a second phase of average wing deterioration was observed (Pearson Correlation: p = 0.021, R2 = 0.69, Radj2 = 0,63) (Fig. 3). This tripartition was also supported by a fourth-degree polynomic function obtaining the best adjusted coefficient of determination (Radj2 = 0.60) and coefficient of determination (R2 = 0.66) (Additonal file 2: Table S2).

Fig. 3

Age structure calculated based on the wing conditions of B. pales. Days with less than five data were excluded (no analysed data prior to 10 July); black trend line for weighted mean value of the total population (fourth-degree polynomic function); coloured linear trend lines for the three phases of the aging process (yellow: 10 July–03 August, green: 03 August–25 August, orange: 25 August–02 September). Only data of wing decay from the first to the second capture event of each individual were used for the calculations

Hence, deterioration in wing conditions over the observed flight period can be divided into three different phases (Fig. 3), which are well supported by the development of the Net Birth rates and the rates of proportional recruitment (pent) calculated with (Fig. 4; third degree polynomic function obtaining the best adjusted coefficient of determination (Radj2) for both values); both rates show long initial increases with maxima around 15 August and rapid decrease thereafter:

  • Initial phase: Slight deterioration of average wing condition. Most individuals are still young and their mortality rate is comparatively low. Therefore, the replacement of the few older individuals by newly emerging individuals cannot compensate the overall ageing of the already existing individuals. (emergence > mortality).

  • Equilibrium phase: Dying older individuals are directly replaced by the large number of newly emerging individuals, which completely compensate the ageing process of the existing individuals. As a consequence, the population has a mostly constant age structure (emergence = mortality).

  • Ageing phase: Due to the small number of newly emerging individuals, the ageing process of the population is not compensated. The high mortality rate of older individuals is not sufficient to counteract this ageing process (mortality > > emergence).

Fig. 4

Net Birth rate and rate of proportional recruitment (pent) and of B. pales calculated with POPAN 5.0 in the program MARK. Red trend line for weighted mean value of pent (third-degree polynomic function); blue trend line for weighted mean value of Net Birth (third-degree polynomic function)

Equation, coefficient of determination (R2) and adjusted coefficient of determination (Radj2) for trend lines of the fourth-degree polynomic function:

y = 8E-07 × 4–0.14 × 3 + 8725.23 × 2 - 2E+ 08x + 3E+ 12; R2 = 0.66; Radj2 = 0.60.

Yellow: R2 = 0.59; Radj2 = 0,54; p = 0.004.

Green: R2 = 0.01; Radj2 = 0,00; p = 0.734.

Orange: R2 = 0.69; Radj2 = 0,63; p = 0.021.

Equations, adjusted coefficients of determination (Radj2) for trend lines of the third-degree polynomic function:

Net Birth rate: y = 0.0007 × 3 + 84.57 × 2 - 4E+06x + 5E+ 10; Radj2 = 0.29.

Rate of proportional recruitment (pent): y = −2E-06 × 3 + 0.2005 × 2–8534.2x + 1E+ 08; Radj2 = 0.35.

Behaviour and foraging

Significant behavioural difference existed between sexes. Thus, females fed significantly more often than males, whereas males flew significantly more often than females; a similarly portion of both sexes was resting when sampled (χ2-homogeneity test: χ2 = 16.60, df = 3, p = 0.001). A significant behavioural difference between sexes was also observed when only the grazed sector was considered (χ2 = 12.18, df = 3, p = 0.007), but not for the ungrazed sector alone (χ2 = 3.86, df = 3, p = 0.277). In general, the individuals are more often flying on the grazed than on the ungrazed sector, whereas the individuals on the ungrazed sector were more frequently resting (males) or feeding (females) (χ2 = 44.30, df = 3, p < 0.001). This difference was also significant for each sex separately (males: χ2 = 33.11, df = 3, p < 0.001; females: χ2 = 10.54, df = 3, p = 0.014) (Table 2).

Table 2 Percentage of individuals of B. pales in four different behavioural categories


Butterflies fed only on plants belonging to the families Asteraceae and Caprifoliaceae, with a clear preference of males for Asteraceae (84 of 93 visits) and of females for Caprifoliaceae (61 of 86 visits). The predominant nectar source was Hieracium aurantiacum for males (73 of 93 visits) and Scabiosa lucida for females (62 of 88 visits). This sex-related difference in foraging behaviour was significant (χ2-homogeneity test: plant families χ2 = 72.48, df = 3, p < 0.001; genera χ2 = 84.04, df = 4, p < 0.001).

Comparing the observed and expected (from plant social inventories and the assumption of random flower-use) nectar sources revealed significant differences for both sexes (males: χ2 = 297.20, df = 33, p < 0.001; females: χ2 = 1171.14, df = 33, p < 0.001). Males preferred the genera Hieracium (73 of 93 visits) in the Asteraceae, whereas the genera Scabiosa and Leontodon were used as expected (cf. Jacobs` index and Bailey’s confidence intervals in Table 3). For females, preference existed for the genus Scabiosa, whereas there were no preferences or even avoidance for all genera of the Asteraceae observed during our study (cf. Jacobs` index and Bailey’s confidence intervals in Table 3).

Table 3 Results for the selection of nectar sources of Boloria pales using Jacobs’ index of selection (Jacobs 1974; interpretational classification of Jacobs` index for our study: 1 to 0.33 preference, 0.33 to − 0.33 neutrality, − 0.33 to −1 avoidance) and Bailey’s confidence intervals at p value < 0.05 (Bailey, 1980; Cherry, 1996); rating: “+” preference nectar source, “=” neutrality, “-”avoidance

Mobility and movement patterns

Most of the B. pales individuals were sedentary. Especially on the ungrazed sector, 93% of the males and 69% of the females did not move more than 100 m. Although both sexes flew larger distances on the grazed sector, most of the females were also sedentary and did not fly more than 200 m (i.e. 76%). In contrast, males at the grazed sector were more flight-active, a larger percentage moved distances exceeding 200 m (i.e. 46%), and even movements of more than 700 m were frequently observed (Fig. 5a and b). The longest observed distance moved was 1096 m for males and 654 m for females. Both distances were observed starting at the grazed sector, whereby also large distances were observed when individuals started on the ungrazed sector (males: 810 m from the first to the second recapture event; females: 615 m from the first capture to the first recapture event). Hence, long movement distances were observed for both sexes on both sectors.

Fig. 5

Percentage of recaptured individuals of B. pales in combination with their movement distances between capture and first recapture event. The distances are divided into 50 m intervals; (a) males, (b) females. Green bars represent the ungrazed sector, orange bars represent the grazed sector

Consequently, the mean distance moved by males (288 m ± 30 SD) was significantly longer than for females (138 m ± 22 SD) on the grazed sector (U test: p = 0.002). On the ungrazed sector, no significant difference between males (58 m ± 8 SD) and females (73 m ± 18 SD) was observed (U test: p = 0.914). Nevertheless, in contrast to the grazed sector, the value for females was higher than for males. Analysing the sexes separately, the mean distance moved was significantly longer on the grazed than on the ungrazed sector for males (U tests: p < 0.001), but not for females (U test: p = 0.045). The critical p value with Bonferroni correction was 0.05 / 2 = 0.025).

The highest stability indices (i.e. Radj2) of fitting the inverse cumulative proportion values of individuals moving certain distance classes to NEF and IPF were obtained for 50 m intervals, for both sexes and both sectors (exception: NEF for females on the grazed sector; males have similar values for all three intervals on the ungrazed sector). On the grazed sector, NEF had higher Radj2 values than IPF in all cases, whereas IPF always had higher Radj2 values on the ungrazed sector (Table 4). Based on both algorithms and 50 m intervals, we estimated the proportion of individuals moving distances of 1 km, 2 km, 3 km, 5 km and 10 km (Table 5).

Table 4 Adjusted Stability Index (Radj2) for IPF and NEF based on 20 m, 30 m and 50 m intervals calculated with movement distances of B. pales
Table 5 Percentage share of individuals of B. pales which would travel 1 km, 2 km, 3 km, 5 km or 10 km; calculated with the NEF and IPF based on 50 m intervals

For comparing the estimated transition probabilities (psi), we used Mann-Whitney U tests because none of the data was normally distributed and all samples from the calculations of the most plausible model of the Multi-state Recaptures (males second model, females first model; Table 6) were independent from each other. Significant differences (critical p value with Bonferroni correction: 0.05 / 4 = 0.0125) between sexes were observed in all cases (U test: transition ungrazed to grazed sector: males 0.060 ± 0.020 vs. females 0.012 ± 0.008, p < 0.001; transition grazed to ungrazed: males 0.100 ± 0.020 vs females 0.029 ± 0.027, p = 0.008; complete transition: males 0.079 ± 0.030 vs. females 0.020 ± 0.017, p < 0.001). The transition probabilities between sectors were not significantly different for females (U tests: p = 0.737) and for both sexes combined (U test: transition ungrazed to grazed sector 0.035 ± 0.021 vs. transition grazed to ungrazed sector 0.065 ± 0.060, p = 0.068). However, the transition probabilities between sectors was significantly different for males (U test: p < 0.001; all values are means ± SD) with more males taking the direction to the ungrazed sector than the reverse.

Table 6 Comparison of the best models of the Multi-state Recaptures; only analyses for estimating the daily population sizes of B. pales in 2014 in Parcul Natural Bucegi, Romania: Akaike information criterion (AICC) and number of considered parameters


Adaptations to high mountain ecosystems

The time-shifted appearance of males and females within a species [6, 54] is a common phenomenon in many Lepidoptera species, e.g. Limenitis camilla [78], Bupalus piniarius [79], Gonepteryx rhamni [80], Euphydryas desfontainii [81]. There are different observations about the consequences of protandry for butterfly population. On the one hand, disadvantages for butterflies are known, i.e. a temporal mismatch between males and females [82, 83], but most studies support the benefits for males [30] and females [29, 84]. Whether protandry is beneficial or not might strongly depend on the degree of adaptation to the local environment. Hence, its existence in a butterfly population is assumed to be an indicator for its adaptation to the prevailing (climatic) conditions and resources. For example, studies with different populations of Euphydryas aurinia demonstrated that well-adapted lowland populations of the species showed a protandrous phenology [33, 53, 84, 85], whereas a population at high altitudes in the eastern Alps did not [31]; this non-protandrous population most likely represented the upper edge of the species’ ecological niche.

In our study, protandry was observed in the surveyed Southern Carpathians B. pales population, but only a comparably small number of males emerged prior to the females (Fig. 1a), as also observed in the High Tauern National Park in the eastern Alps [39]. Nevertheless, a sufficient number of males was present at the beginning of the females’ emergence so that both sexes can benefit from the known advantages of protandry [28,29,30]. A quite similar phenology was also observed for Erebia nivalis in the eastern Alps [40], where this species inhabits extreme high-altitude habitats. Apparently, conditions in really extreme habitats may be so unpredictable that protandry must be softened to guarantee the survival of butterfly populations in more hostile years and not risking a complete temporal mismatch between both sexes.

This “soft” protandry in the Bucegi B. pales population was combined with a month-long phase of mostly constant age structure (Fig. 2), a phenomenon also observed in an eastern Alpine population of the species [39]. This phase of constant age structure was due to an extended emergence period (reflected in the Net Birth rates and the rate of proportional recruitment) resulting in a relatively low average wing damage over time if compared to other butterfly species with mostly simultaneous emergence and hence continuously ageing population structures over the entire observed flight period [34, 37, 54, 86,87,88]. This capability of constant rejuvenation reduces the species’ susceptibility to short-term bad weather events such as cold snaps [3], one main trigger of increased mortality rates in butterflies [89], which hence reduces the time available for successful reproduction [90, 91]. This form of risk spreading (i.e. bet-hedging [92,93,94]) against the unpredictable conditions in high mountain environments, combined with the observed “soft” protandry calls for a generally good adaptation of the phenology of the local populations of B. pales to the climatic conditions prevailing in the Southern Carpathians.

Sex-dependent behavioural differences

Comparing the general behavioural and movement traits between sexes, the distances moved by males were significantly larger than for females (Fig. 5 and Table 5), which consequently were observed flying much less, but visited flowers to obtain nectar much more often than males (Tab. 5). Similar behavioural differences between sexes are a common phenomenon in butterflies [35, 37, 41,42,43].

These differences in flight activity and movement patterns, in combination with the different roles of the sexes (males: finding females; females: successful egg production and oviposition), might also be responsible for the sexual dimorphism in nectar source selection observed in our study, most likely driven by differential nutritional needs. In the Bucegi Mountains (this study) and in the Alps [5], B. pales exhibited high specialisation on some few plant species, even down to the species level, but with a strong difference between sexes.

As already known, the selection of nectar sources by butterflies mainly depends on the offered amounts of hexoses, sucrose, fructose and amino-acids, with sugars being more important for flight activity and amino-acids more for egg production [44,45,46]. Nevertheless, the sex-dependent preferences and differential needs for nectar ingredients continue to be discussed controversially. This specialisation could become problematic in times of climate change, in particular in high mountain regions because it might lead to shifts in the timing of plant flowering and butterfly appearances [95], causing problems for butterflies to find sufficient nectar plants if their phonologies are disengaged from the flowering periods of their preferred nectar sources.

Mobility and the impact of habitat quality

Habitat quality is influenced by abiotic and biotic factors [96, 97]. Consequently, human activities have paramount importance today. For example, a non-adapted grazing intensity often produces pastures with a low diversity and abundance of flowering plants [98], a situation also observable in our study area in the Bucegi Mountains [26, 27]. Consequently, negative effects of grazing on population sizes and general diversity of invertebrate species have frequently been reported [99,100,101]. These negative grazing effects are mostly linked to the reduction of essential resources, such as foraging plants or places for oviposition [102]. In our study, the intensively pastured parts of the study area also had a remarkably lower population density of B. pales than the grazed sector. This difference is reflecting a higher habitat quality of the ungrazed sector. Nevertheless, this difference might not be solely the result of the different grazing regimes, but other factors like the slightly different altitudes of both sectors and hence some subtle differences in the climatic (and also other environmental) conditions might have additional influence.

However, habitat quality is not only of paramount importance for the persistence of butterfly species in an area [103,104,105], but also strongly influences the individuals’ behaviour. Thus, movement activities in general increase with a decrease in habitat quality [106, 107]. This also applies to high mountain butterfly species. For example, differences in the movement behaviour in high and low quality habitats were observed for two Erebia species in the Hrubý Jeseník Mountains. Individuals of E. epiphron tended to be more sedentary in areas with high habitat quality, and also individuals of E. sudetica had a stronger dispersal behaviour in habitats of lower quality [48]. Similar results are known for Parnassius smintheus, a species which showed restricted movements on large homogenous alpine grasslands without obvious barriers (i.e. optimum habitats), but considerably higher mobility on low quality patches [12, 48, 108, 109]. In this context, the differences between both sectors observed in our study with respect to population density (Fig. 2a and b) and movement behaviour (Fig. 5a and b) of B. pales are fully in line with comparable observations made on other mountain butterfly species. Especially the males on the grazed sector showed long movements and consequently a high potential for long-distance dispersal (Table 5) if compared to males of other alpine butterfly species [39].

More specifically, each of the individuals in our study, independently of its sex, had the option of either selecting a sector with lower habitat quality caused by intensive grazing but also lower competition for oviposition places and, most likely, nectar sources, or a sector with higher habitat quality but also higher general competition.

As the movement between both sectors was not clearly favouring one of the two possible directions, we assume that the habitat quality was still adequate for both sexes on both sectors. This was in clear contrast to observations in other butterfly species, such as Maculinea nausithous and M. teleius [36], that showed a clear migration in the direction of the more favourable habitats. Nevertheless, although both sectors were used by B. pales, the negative consequences of intensive grazing were obvious. Under the initially probably similar abiotic and biotic conditions in both sectors, similar population densities could be expected to have occurred on them (but see above). Consequently, the observed differences were most likely caused by the effects of grazing, thus highlighting the strong impact of non-adapted management on the pastures. Except for some small patches of suitable habitats, most parts of the grazed sector were avoided by butterflies, with important consequences on behaviour and dispersal (ES2).

Conservation implications

Until the end of the 1990s, land-use in the Southern Carpathians was mostly traditional, involving transhumance. The natural and semi-natural grasslands, supported by this type of use, are of great importance for nature and biodiversity conservation [23]. Structural changes in politics, economy and society have led to an intensification of land-use (e.g. by decoupling the sheep per hectare ratio [110]), causing severe conservation problems, especially for natural and semi-natural grasslands [111].

Hence, other studies in the Bucegi Mountains already have shown changes in the plant composition of this area, where a decline in the grass species Nardus stricta was observed in an 18-year study. This (formerly dominant) species was increasingly replaced not only by other forage grasses, such as Agrostis capillaris and Poa pratensis, but also flowering plants like Trifolium repens [27]. Furthermore, the intensive sheep pasturing in the Bucegi Mountains has reduced the number of endemic species as well as the biodiversity and abundance of plants and gastropods [26]. Beside these biotic consequences, negative effects on abiotic factors caused by intensive grazing have also been observed, such as increased water pollution and progressive soil erosion [110]. Such negative influences of intensive grazing on the diversity and abundance of insect groups are also known from other studies [100, 112]. This also applies to other high altitude butterfly species, like the south-eastern Alpine endemic Erebia calcaria, which is negatively affected by intensive grazing and fragmentation [99].

The land-use changes in the Bucegi Mountains have not only led to local resource scarcity and habitat fragmentation, but most likely have also caused ehavioural changes in the studied population of B. pales, evident in the different movement patterns within the grazed and the ungrazed sectors. In general, movement patterns are influenced by diverse factors, like the trade-off between costs and success of a movement, or the combined effects of internal and external information on individuals [113, 114]. As our study was performed at sectors with originally, most likely, similar biotic and abiotic conditions, we assume that also the observed differences result from the differences in grazing intensity. In combination with the previously mentioned studies, we assume that the current level of grazing intensity in the Bucegi Mountains negatively affects (insect) species and biodiversity in general.

To hinder further disappearance of these high mountain species, conservation measures are necessary. One possible instrument is the Carpathian Convention, a treaty to foster the sustainable development and the conservation of the Carpathian region, signed by seven Carpathian states in 2003 [115]. This convention suggested the implementation of traditional land-use in a sustainable manner, with a focus on „the protection of mountain ecosystems, […], the importance of biological diversity, and the specific conditions of mountains as less favoured areas (Carpathian Convention Art. 7, Paragraph 1). Although this article is broadly defined, it is apparent that sustainable development of the Carpathians is not possible without modifications in land-use intensity towards a grazing with a lower intensity. In the context of high mountain ecosystems, this entails a traditional low-intensity land-use, but not abandonment. Another possible instrument is the European Union Law, in particular the Common Agricultural Policy (CAP). This policy serves to support farms and to ensure the ration of agricultural goods (CAP, Art. 39). However, the actual Common Agricultural Policy 2014–2020 is not suitable to support an appropriate and sustainable agriculture as well as to act adequately against the trend of abandonment of farmland in mountain areas. On the one hand, the ecological demands of the system of benefits is not sufficient to achieve an adapted grazing in mountain areas (e.g. Greening and Cross Compliance; Regulation (EU) No 1307/2013). On the other hand, an insufficient amount of money is spent supporting rural development by the European Agricultural Fund for Rural Development (Regulation (EU) No 1305/2013) to prevent rural migration. This is aggravated by the lack of willingness of most EU member states, including Romania, to use the possibility of an obligatory payment for areas with natural constraints (Art. 46 Regulation (EU) No 1307 /2013). This could be an additional source of income for Romanian farmers in high mountain areas. Hence, it would be desirable to strengthen the situation of mountain farming in the next Common Agricultural Policy 2021–2028 to ensure a permeant and sustainable development of the Carpathian meadows. At the same time, further funding resources would help here to realize nature conservation through concrete management actions. In this context and from a European perspective, so-called LIFE projects would be useful supplements to the law principles; unfortunately, no LIFE project is currently in progress to save extensively used Carpathians meadows;

Several studies have shown that light grazing even increases the number and abundance of butterfly species, whereas total long-term abandonment with subsequent succession on the long run reduces butterfly diversity [116,117,118,119]. A fast and complete implementation of the Carpathian Convention as well as modifications of the Common Agricultural Policy in the next funding period would therefore be of great importance in counteracting the recently observed habitat degradation, fragmentation and diversity loss.


Specific adaptations for risk-spreading of B. pales to cope with the impact of unpredictable weather conditions of the Carpathian Mountains were revealed in this study. As in the Alps [39], we observed an extended emergence of individuals over the observed flight period, dividing the population’s age structure into three phases. In combination with the detected “soft” form of protandry, in which only a small number of males appeared before females, these adaptations enable the species to conquer short bad weather events, because only the part of the whole population being at the wing by that time might be killed. Comparing both sexes, movement differences between males and females were obvious; male individuals flew larger distances and had a higher flight-activity. This might explain the need of different nectar ingredients, which consequently leads to a dimorphism in foraging behaviour: males preferred nectar sources of Asteraceae, females Caprifoliaceae. Although we were not able to study the ungrazed sector as long as the grazed sector, a deterioration of the habitat quality caused by intensive grazing on the grazed sector in comparison to the ungrazed sector was likely. The permanent grazing might have negatively influenced population density, causing significantly larger flight distances of the individuals and different behaviour at both sectors. Comparably negative results of intensive grazing are already known for other insect groups [100, 112]. Hence, the ongoing land-use change and intensification of grazing in the Southern Carpathians, here in particular the Bucegi Mountains [26, 27], are a threat to the native insects in the whole mountain area. To counter this progress, improvements in the actual management and jurisdiction are needed, to protect these insect species from the consequences of overgrazing (e.g. by establishing areas without or just light grazing). In this context, this study about B. pales might provide essentials for further research and a better understanding of high alpine butterflies in the Carpathian Mountains.



Akaike Information Criterion




Common Agricultural Policy


Electronic Supplement


European Union


Inverse power function

m asl:

Metres above sea level




Negative exponential function



R2 :

Coefficient of determination

Radj 2 :

Adjusted coefficient of determination


Standard deviation


Standard error


  1. 1.

    Kocher TD. Adaptive evolution and explosive speciation: the cichlid fish model. Nat Rev Genet. 2004;5:288–98.

    CAS  Article  PubMed  Google Scholar 

  2. 2.

    Kawecki TJ, Ebert D. Conceptual issues in local adaptation. Ecol Lett. 2004;7:1225–41.

    Article  Google Scholar 

  3. 3.

    Altshuler DL, Dudley R. Adaptations to life at high elevation: an introduction to the symposium. Integr Comp Biol. 2006;46:3–4.

    Article  PubMed  Google Scholar 

  4. 4.

    Hodkinson ID. Terrestrial insects along elevation gradients: species and community responses to altitude. Biol Rev Camb Philos Soc. 2005;80:489–513.

    Article  Google Scholar 

  5. 5.

    Ehl S, Hostert K, Korsch J, Gros P, Schmitt T. Sexual dimorphism in the alpine butterflies Boloria pales and Boloria napaea: differences in movement and foraging behavior (Lepidoptera: Nymphalidae). Insect Sci. 2017.

  6. 6.

    Schtickzelle N, Le Boulenge E, Baguette M. Metapopulation dynamics of the bog fritillary butterfly: demographic processes in a patchy population. Oikos. 2002;97:349–60.

    Article  Google Scholar 

  7. 7.

    Lederhouse RC. Population structure, residency and weather related mortality in the black swallowtail butterfly, Papilio polyxenes. Oecologia. 1983;59:307–11.

    Article  PubMed  Google Scholar 

  8. 8.

    Duplouy A, Hanski I. Butterfly survival on an isolated island by improved grip. Biol Lett. 2013;9:20130020.

    Article  PubMed  PubMed Central  Google Scholar 

  9. 9.

    Diamond SE, Cayton H, Wepprich T, Jenkins CN, Dunn RR, Haddad NM, Ries L. Unexpected phenological responses of butterflies to the interaction of urbanization and geographic temperature. Ecology. 2014;95:2613–21.

    Article  Google Scholar 

  10. 10.

    Talavera G, Vila R. Discovery of mass migration and breeding of the painted lady butterfly Vanessa cardui in the sub-Sahara: the Europe-Africa migration revisited. Biol J Linn Soc. 2016.

  11. 11.

    Willmott KR, Robinson Willmott JC, Elias M, Jiggins CD. Maintaining mimicry diversity: optimal warning colour patterns differ among microhabitats in Amazonian clearwing butterflies. Proc Biol Sci. 2017.

  12. 12.

    Roland J, Keyghobadi N, Fownes S. Alpine Parnassius butterfly dispersal: effects of landscape and population size. Ecology. 2000;81:1642–53.;2.

    Article  Google Scholar 

  13. 13.

    Roland J, Matter SF. Pivotal effect of early-winter temperatures and snowfall on population growth of alpine Parnassius smintheus butterflies. Ecol Monogr. 2016;86:412–28.

    Article  Google Scholar 

  14. 14.

    Matter SF, Roland J. A priori prediction of an extreme crash in 2015 for a population network of the alpine butterfly, Parnassius smintheus. Ecosphere. 2015;6:art178.

    Article  Google Scholar 

  15. 15.

    Tasser E, Tappeiner U. Impact of land use changes on mountain vegetation. Appl Veg Sci. 2002;5:173–84.

    Article  Google Scholar 

  16. 16.

    Zöchling A, Fiedler K. Einfluss der Nutzungsintensität von Almen auf Tagfalter-Gemeinschaften (Lepidoptera: Papilionoidea) im Nationalpark Gesäuse (Österreich). Entomologica Austriaca. 2015;22:7–26.

    Google Scholar 

  17. 17.

    Chemini C, Rizzoli A. Land use change and biodiversity conservation in the Alps. Journal of Mountain Ecology. 2003;7.

  18. 18.

    Tappeiner U, Tasser E, Leitinger G, Tappeiner GL i d A. historische Entwicklung und zukünftige Szenarien. In: Psenner R, Lackner R, editors. Die Alpen im Jahr 2020. Innsbruck: innsbruck university press; 2006. p. 23–39.

    Google Scholar 

  19. 19.

    Hill JK, Thomas CD, Fox R, Telfer MG, Willis SG, Asher J, Huntley B. Responses of butterflies to twentieth century climate warming: implications for future ranges. Proc Biol Sci. 2002;269:2163–71.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  20. 20.

    Lafranchis T, Jutzeler D, Guillosson J-Y, Kan P, Kan B. La vie des papillons: Écologie, biologie et comportement des Rhopalocères de France. In: [Paris]: Diatheo; 2015.

    Google Scholar 

  21. 21.

    Tolman T, Lewington R. Schmetterlinge Europas und Nordwestafrikas: Alle Tagfalter; über 2000 Arten. 2nd ed. Stuttgart: Kosmos; 2012.

    Google Scholar 

  22. 22.

    Kudrna O, Pennerstorfer J, Lux K. Distribution atlas of European butterflies and skippers. Schwanfeld: Wissenschaftlicher Verlag Peks; 2015.

  23. 23.

    Schmitt T, Rákosy L. Changes of traditional agrarian landscapes and their conservation implications: a case study of butterflies in Romania. Divers Distrib. 2007;13:855–62.

    Article  Google Scholar 

  24. 24.

    Cuvelier S. Dinca V. new data regarding the butterflies (Lep.: Rhopalocera) of Romania, with additional comments (general distribution in Romania, habitat preferences, threats and protection) for ten localized Romanian species. Phegea. 2007;35:93–115.

    Google Scholar 

  25. 25.

    Rákosy L, Goia M, Kovács Z. Catalogul lepidopterelor României. Cluj-Napoca: Societatea lepidopterologică Romană; 2003.

  26. 26.

    Baur B, Cremene C, Groza G, Schileyko A, Baur A, Erhardt A. Intensified grazing affects endemic plant and gastropod diversity in alpine grasslands of the southern Carpathian mountains (Romania). Biologia. 2007.

  27. 27.

    Maruşca T, Dragomir N, Mocanu V, Blaj A, Tarjoc F, Dragomir C, Constantinescu S. Effect of some improvement works on the floristic composition of the vegetal cover in Nardus stricta grasslands. Scientific Papers: Animal Science and Biotechnologies. 2014;47.

  28. 28.

    Fagerström T, Wiklund C. Why do males emerge before females? Protandry as a mating strategy in male and female butterflies. Oecologia. 1982;52:164–6.

    Article  Google Scholar 

  29. 29.

    Larsen E, Calabrese JM, Rhainds M, Fagan WF. How protandry and protogyny affect female mating failure: a spatial population model. Entomol Exp Appl. 2013;146:130–40.

    Article  Google Scholar 

  30. 30.

    Nève G, Singer MC. Protandry and postandry in two related butterflies: conflicting evidence about sex-specific trade-offs between adult size and emergence time. Evol Ecol. 2008;22:701–9.

    Article  Google Scholar 

  31. 31.

    Junker M, Wagner S, Gros P, Schmitt T. Changing demography and dispersal behaviour: ecological adaptations in an alpine butterfly. Oecologia. 2010;164:971–80.

    Article  PubMed  Google Scholar 

  32. 32.

    Junker M, Zimmermann M, Ramos AA, Gros P, Konvička M, Nève G, et al. Three in one--multiple faunal elements within an endangered European butterfly species. PLoS One. 2015;10:e0142282.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  33. 33.

    Junker M, Schmitt T. Demography, dispersal and movement pattern of Euphydryas aurinia (Lepidoptera: Nymphalidae) at the Iberian Peninsula: an alarming example in an increasingly fragmented landscape? J Insect Conserv. 2010;14:237–46.

    Article  Google Scholar 

  34. 34.

    Knight A, Brower LP, Williams EH. Spring remigration of the monarch butterfly, Danaus plexippus (Lepidoptera: Nymphalidae) in north-Central Florida: estimating population parameters using mark-recapture. Biol J Linn Soc. 1999;68:531–56.

    Article  Google Scholar 

  35. 35.

    Slamova I, Klecka J, Konvička M. Diurnal behavior and habitat preferences of Erebia aethiops, an aberrant lowland species of a mountain butterfly clade. J Insect Behav. 2011;24:230–46.

    Article  Google Scholar 

  36. 36.

    Timus N, Czekes Z, Rákosy L, Nowicki P. Conservation implications of source-sink dynamics within populations of endangered Maculinea butterflies. J Insect Conserv. 2016.

  37. 37.

    Weyer J, Schmitt T. Knowing the way home: strong philopatry of a highly mobile insect species, Brenthis ino. J Insect Conserv. 2013;17:1197–208.

    Article  Google Scholar 

  38. 38.

    Szentirmai I, Mesterházy A, Varga I, Schubert Z, Sándor LC, Ábrahám L, Kőrösi Á. Habitat use and population biology of the Danube clouded yellow butterfly Colias myrmidone (Lepidoptera: Pieridae) in Romania. J Insect Conserv. 2014;18:417–25.

    Article  Google Scholar 

  39. 39.

    Ehl S, Ebertshäuser M, Gros P, Schmitt T. Population demography of alpine butterflies: Boloria pales and Boloria napaea (Lepidoptera: Nymphalidae) and their specific adaptations to high mountain environments. Acta Oecol. 2017;85:53–61.

    Article  Google Scholar 

  40. 40.

    Ehl S, Dalstein V, Tull F, Gros P, Schmitt T. Specialized or opportunistic-how does the high mountain endemic butterfly Erebia nivalis survive in its extreme habitats? Insect Sci. 2016.

  41. 41.

    Wiklund C, Åhrberg C. Host plants, nectar source plants, and habitat selection of males and females of Anthocharis cardamines (Lepidoptera). Oikos. 1978;31:169.

    Article  Google Scholar 

  42. 42.

    Zimmermann K, Fric Z, Filipova L, Konvicka M. Adult demography, dispersal and behaviour of Brenthis ino (Lepidoptera: Nymphalidae): how to be a successful wetland butterfly. Eur. J. Entomol. 2005;102:699–706.

    Article  Google Scholar 

  43. 43.

    Zimmermann K, Konvička M, Fric Z, Čihaková V. Demography of a common butterfly on humid grasslands: Argynnis aglaja (Lepidoptera: Nymphalidae) studied by mark-recapture. Pol J Ecol. 2009;57:715–27.

    Google Scholar 

  44. 44.

    Mevi-Schütz J, Erhardt A. Can Inachis io detect nectar amino acids at low concentrations? Physiol Entomol. 2002;27:256–60.

    Article  Google Scholar 

  45. 45.

    Rusterholz H-P, Erhardt A. Preferences for nectar sugars in the peacock butterfly, Inachis io. Ecol Entomol. 1997;22:220–4.

    Article  Google Scholar 

  46. 46.

    Rusterholz H-P, Erhardt A. Can nectar properties explain sex-specific flower preferences in the Adonis blue butterfly Lysandra bellargus? Ecol Entomol. 2000;25:81–90.

    Article  Google Scholar 

  47. 47.

    Neumayer J, Spaethe J. Blütenfarbe, Nektarmenge und Blütenbesuche von Schmetterlingen in einem alpinen Lebensraum in Zentraleuropa. entomologia. 2007;29:269–84.

    Article  Google Scholar 

  48. 48.

    Kuras T, Benes J, Fric Z, Konvicka M. Dispersal patterns of endemic alpine butterflies with contrasting population structures: Erebia epiphron and E. sudetica. Popul Ecol. 2003;45:115–23.

    Article  Google Scholar 

  49. 49.

    Tshikolovets VV. Butterflies of Europe & the Mediterranean Area Pardubice: Tshikolovets Publications; 2011.

  50. 50.

    Huemer P, Wieser C. Nationalpark Hohe Tauern. Schmetterlinge; [wissenschaftliche Schriften]. Innsbruck. In: Wien: Tyrolia; 2008.

    Google Scholar 

  51. 51.

    Schweizerischer Bund für Naturschutz. Tagfalter und ihre Lebensräume. Arten, Gefährdung, Schutz : Schweiz und angrenzende Gebiete. 4th ed. Basel: Schweizerischer Bund für Naturschutz; 1994.

    Google Scholar 

  52. 52.

    Thomas JA. The ecology and conservation of Lysandra bellargus (Lepidoptera: Lycaenidae) in Britain. The Journal of Applied Ecology. 1983;20:59.

    Article  Google Scholar 

  53. 53.

    Munguira ML, Martín J, García-Barros E, Viejo JL. Use of space and resources in a Mediterranean population of the butterfly Euphydryas aurinia. Acta Oecol. 1997;18:597–612.

    Article  Google Scholar 

  54. 54.

    Watt WB, Chew FS, Synder LRG, Watt AG, Rothschild DE. Population structure of pierid butterflies. Oecologia. 1977;27:1–22.

    Article  Google Scholar 

  55. 55.

    Gall LF. The effects of capturing and marking on subsequent activity in Boloria acrocnema (Lepidoptera: Nymphalidae), with a comparison of different numerical models that estimate population size. Biol Conserv. 1984;28:139–54.

    Article  Google Scholar 

  56. 56.

    Tabashnik BE. Population structure of pierid butterflies: III. Pest populations of Colias philodice eriphyle. Oecologia. 1980;47:175–83.

    Article  Google Scholar 

  57. 57.

    Iwasa Y, Odendaal FJ, Murphy DD, Ehrlich PR, Launer AE. Emergence patterns in male butterflies: a hypothesis and a test. Theor Popul Biol. 1983;23:363–79.

    Article  Google Scholar 

  58. 58.

    Pflanƶensoziologie B-BJ. Grundƶüge der Vegetationskunde. Vienna, s.l. Vienna: Springer; 1964.

    Google Scholar 

  59. 59.

    Mühlenberg M, Bogenrieder A, Behre GF. Freilandökologie. 3rd ed. Heidelberg: Quelle & Meyer; 1993.

    Google Scholar 

  60. 60.

    Jacobs J. Quantitative measurement of food selection: a modification of the forage ratio and Ivlev's electivity index. Oecologia. 1974;14:413–7.

    Article  PubMed  Google Scholar 

  61. 61.

    Cherry S. A comparison of confidence interval methods for habitat use-availability studies. J Wildl Manag. 1996;60:653.

    Article  Google Scholar 

  62. 62.

    Bailey BJR. Large sample simultaneous confidence intervals for the multinomial probabilities based on transformations of the cell frequencies. Technometrics. 1980;22:583–9.

    Article  Google Scholar 

  63. 63.

    IBM Cooperation. SPSS Statistics 22: SPSS Statistics 22 (IBM).

  64. 64.

    Cooch EG, White GC. Program MARK - a gentle introduction. 2016.

    Google Scholar 

  65. 65.

    Bellebaum J, Köppen U, Grajetzky B. Ermittlung von Überlebenswahrscheinlichkeiten aus Ringfunddaten. Vogelwarte. 2010;48:21–32.

    Google Scholar 

  66. 66.

    White GC, Burnham KP. Program MARK: survival estimation from populations of marked animals. Bird Study. 1999;46:120–39.

    Article  Google Scholar 

  67. 67.

    Kadlec T, Vrba P, Kepka P, Schmitt T, Konvicka M. Tracking the decline of the once-common butterfly: delayed oviposition, demography and population genetics in the hermit Chazara briseis. Anim Conserv. 2010;13:172–83.

    Article  Google Scholar 

  68. 68.

    Arnason AN. Parameter estimates from mark-recapture experiments on two populations subject to migration and death. Res Popul Ecol. 1972;13:97–113.

    Article  Google Scholar 

  69. 69.

    Arnason AN. The estimation of population size, migration rates and survival in a stratified population. Res Popul Ecol. 1973;15:1–8.

    Article  Google Scholar 

  70. 70.

    Brownie C, Hines JE, Nichols JD, Pollock KH, Hestbeck JB. Capture-recapture studies for multiple strata including non-Markovian transitions. Biometrics. 1993;49:1173.

    Article  Google Scholar 

  71. 71.

    Schwarz CJ, Schweigert JF, Arnason AN. Estimating migration rates using tag-recovery data. Biometrics. 1993;49:177.

    Article  Google Scholar 

  72. 72.

    Sugiura N. Further analysts of the data by akaike' s information criterion and the finite corrections. Communications in Statistics - Theory and Methods. 2007;7:13–26.

    Article  Google Scholar 

  73. 73.

    Hurvich CM, Tsai C-L. Regression and time series model selection in small samples. Biometrika. 1989;76:297–307.

    Article  Google Scholar 

  74. 74.

    Burnham KP, Anderson DR. Model selection and multimodel inference: a practical information-theoretic approach. 2nd ed. New York: Springer; 2002.

    Google Scholar 

  75. 75.

    Cook LM, Brower LP, Croze HJ. The accuracy of a population estimation from multiple recapture data. J Anim Ecology. 1967;36:57.

    Article  Google Scholar 

  76. 76.

    Baguette M. Long distance dispersal and landscape occupancy in a metapopulation of the cranberry fritillary butterfly. Ecography. 2003;26:153–60.

    Article  Google Scholar 

  77. 77.

    Fric Z, Konvička M. Dispersal kernels of butterflies: power-law functions are invariant to marking frequency. Basic and Applied Ecology. 2007;8:377–86.

    Article  Google Scholar 

  78. 78.

    Lederer G. Verhaltensweisen der Imagines und der Entwicklungsstadien von Limenitis camilla camilla L. (Lep. Nymphalidae). Z Tierpsychol. 1960;17:521–46.

    Article  Google Scholar 

  79. 79.

    Botterweg PF. Protandry in the pine Looper, Bupalus Piniarius (Lep., Geometridae); an explanatory model. Netherlands Journal of Zoology. 1981;32:169–93.

    Article  Google Scholar 

  80. 80.

    Wiklund C, Lindfors V, Forsberg J. Early male emergence and reproductive phenology of the adult overwintering butterfly Gonepteryx rhamni in Sweden. Oikos. 1996;75:227.

    Article  Google Scholar 

  81. 81.

    Pennekamp F, Garcia-Pereira P, Schmitt T. Habitat requirements and dispersal ability of the Spanish fritillary (Euphydryas desfontainii) in southern Portugal: evidence-based conservation suggestions for an endangered taxon. J Insect Conserv. 2014;18:497–508.

    Article  Google Scholar 

  82. 82.

    Bubová T, Kulma M, Vrabec V, Nowicki P. Adult longevity and its relationship with conservation status in European butterflies. J Insect Conserv. 2016;20:1021–32.

    Article  Google Scholar 

  83. 83.

    Calabrese JM, Ries L, Matter SF, Debinski DM, Auckland JN, Roland J, Fagan WF. Reproductive asynchrony in natural butterfly populations and its consequences for female matelessness. J Anim Ecol. 2008;77:746–56.

    Article  PubMed  Google Scholar 

  84. 84.

    Hula V, Konvicka M, Pavlicko A, Fric Z. Marsh fritillary (Euphydryas aurinia) in the Czech Republic: monitoring, metapopulation structure, and conservation of an endangered butterfly. Entomologica Fennica. 2004;15:231–41.

    Google Scholar 

  85. 85.

    Wang R, Wang Y, Chen J, Lei G, Xu R. Contrasting movement patterns in two species of chequerspot butterflies, Euphydryas aurinia and Melitaea phoebe, in the same patch network. Ecol Entomol. 2004;29:367–74.

    Article  Google Scholar 

  86. 86.

    De-Andrade RB, Freitas AVL. Population biology of two species of Heliconius (Nymphalidae: Heliconiinae) in a semi-deciduous Forest in southeastern Brazil. Journal of lepidopterists´. Society. 2005;59:223–8.

    Google Scholar 

  87. 87.

    Matsumoto K. Population dynamics of Luehdorfia japonica leech (Lepidoptera: Papilionidae): I. a preliminary study on the adult population. Res Popul Ecol. 1984;26:1–12.

    Article  Google Scholar 

  88. 88.

    Arnold RA. Ecological studies of six endangered butterflies (Lepidoptera, Lycaenidae): island biogeography, patch dynamics, and the design of habitat preserves. Berkeley: University of California Press; 1983.

    Google Scholar 

  89. 89.

    Nowicki P, Bonelli S, Barbero F, Balletto E. Relative importance of density-dependent regulation and environmental stochasticity for butterfly population dynamics. Oecologia. 2009;161:227–39.

    Article  PubMed  Google Scholar 

  90. 90.

    Dingle H, Mousseau TA, Scott SM. Altitudinal variation in life cycle syndromes of California populations of the grasshopper, Melanoplus sanguinipes (F.). Oecologia. 1990;84:199–206.

    Article  Google Scholar 

  91. 91.

    Roff D. Optimizing development time in a seasonal environment: the ‘ups and downs’ of clinal variation. Oecologia. 1980;45:202–8.

    Article  Google Scholar 

  92. 92.

    Garcia-Gonzalez F, Yasui Y, Evans JP. Mating portfolios: bet-hedging, sexual selection and female multiple mating. Proc Biol Sci. 2015;282:20141525.

    Article  PubMed  PubMed Central  Google Scholar 

  93. 93.

    Hopper KR. Risk-spreading and bet-hedging in insect population biology. Annu Rev Entomol. 1999;44:535–60.

    CAS  Article  PubMed  Google Scholar 

  94. 94.

    Olofsson H, Ripa J, Jonzén N. Bet-hedging as an evolutionary game: the trade-off between egg size and number. Proc Biol Sci. 2009;276:2963–9.

    Article  PubMed  PubMed Central  Google Scholar 

  95. 95.

    Kharouba HM, Vellend M. Flowering time of butterfly nectar food plants is more sensitive to temperature than the timing of butterfly adult flight. J Anim Ecol. 2015;84:1311–21.

    Article  PubMed  Google Scholar 

  96. 96.

    Weiss SB, Murphy DD, White RR. Sun, slope, and butterflies: topographic determinants of habitat quality for Euphydryas Editha. Ecology. 1988;69:1486–96.

    Article  Google Scholar 

  97. 97.

    Harborne AR, Mumby PJ, Kennedy EV, Ferrari R. Biotic and multi-scale abiotic controls of habitat quality: their effect on coral-reef fishes. Mar Ecol Prog Ser. 2011;437:201–14.

    Article  Google Scholar 

  98. 98.

    Niedrist G, Tasser E, Lüth C, Dalla Via J, Tappeiner U. Plant diversity declines with recent land use changes in European Alps. Plant Ecol. 2009;202:195–210.

    Article  Google Scholar 

  99. 99.

    De Groot M, Rebeusek F, Grobelnik V, Govedic M, Salamun A, Verovnik R. Distribution modelling as an approach to the conservation of a threatened alpine endemic butterfly (Lepidoptera Satyridae). Eur J Entomol. 2009;106:77–84.

    Article  Google Scholar 

  100. 100.

    Kruess A, Tscharntke T. Grazing intensity and the diversity of grasshoppers, butterflies, and trap-nesting bees and wasps. Conserv Biol. 2002;16:1570–80.

    Article  Google Scholar 

  101. 101.

    Bromham L, Cardillo M, Bennett AF, Elgar MA. Effects of stock grazing on the ground invertebrate fauna of woodland remnants. Austral Ecology. 1999;24:199–207.

    Article  Google Scholar 

  102. 102.

    Turlure C, Choutt J, van Dyck H, Baguette M, Schtickzelle N. Functional habitat area as a reliable proxy for population size: case study using two butterfly species of conservation concern. J Insect Conserv. 2010;14:379–88.

    Article  Google Scholar 

  103. 103.

    Thomas JA, Bourn NA, Clarke RT, Stewart KE, Simcox DJ, Pearman GS, et al. The quality and isolation of habitat patches both determine where butterflies persist in fragmented landscapes. Proc R Soc B. 2001;268:1791–6.

    CAS  Article  PubMed  Google Scholar 

  104. 104.

    van Halder I, Barnagaud J-Y, Jactel H, Barbaro L. Woodland habitat quality prevails over fragmentation for shaping butterfly diversity in deciduous forest remnants. For Ecol Manag. 2015;357:171–80.

    Article  Google Scholar 

  105. 105.

    Barbaro L, van Halder I. Linking bird, carabid beetle and butterfly life-history traits to habitat fragmentation in mosaic landscapes. Ecography. 2009;32:321–33.

    Article  Google Scholar 

  106. 106.

    Fernández P, Rodríguez A, Obregón R, de HS, Jordano D, Fernández-Haeger J. Fine scale movements of the butterfly Plebejus argus in a heterogeneous natural landscape as revealed by GPS tracking. J Insect Behav. 2016;29:80–98.

    Article  Google Scholar 

  107. 107.

    Kallioniemi E. Effects of morphology, habitat and weather on the movement behaviour of range-expanding butterfly species. Great Britain: University of East Anglia; 2013.

    Google Scholar 

  108. 108.

    Keyghobadi N, Roland J, Strobeck C. Influence of landscape on the population genetic structure of the alpine butterfly Parnassius smintheus (Papilionidae). Mol Ecol. 1999;8:1481–95.

    CAS  Article  PubMed  Google Scholar 

  109. 109.

    Schmitt T, Hewitt GM, Muller P. Disjunct distributions during glacial and interglacial periods in mountain butterflies: Erebia epiphron as an example. J Evol Biol. 2006;19:108–13.

    CAS  Article  PubMed  Google Scholar 

  110. 110.

    Urushibara-Yoshino K, Mori K. Degradation of Geoecological and hydrological conditions due to grazing in south Carpathian Mountains under the influence of changing social structure in Romania. Geogr Rev JPN, Chirigaku Hyoron. 2007;80:272–89.

    Article  Google Scholar 

  111. 111.

    Huband S, McCracken DI, Mertens A. Long and short-distance transhumant pastoralism in Romania: past and present drivers of change. Pastoralism: Research,Policy and. Practice. 2010;1.

  112. 112.

    Kruess A, Tscharntke T. Contrasting responses of plant and insect diversity to variation in grazing intensity. Biol Conserv. 2002;106:293–302.

    Article  Google Scholar 

  113. 113.

    Bonte D, van Dyck H, Bullock JM, Coulon A, Delgado M, Gibbs M, et al. Costs of dispersal. Biol Rev Camb Philos Soc. 2012;87:290–312.

    Article  PubMed  Google Scholar 

  114. 114.

    Clobert J, Le Galliard J-F, Cote J, Meylan S, Massot M. Informed dispersal, heterogeneity in animal dispersal syndromes and the dynamics of spatially structured populations. Ecol Lett. 2009;12:197–209.

    Article  PubMed  Google Scholar 

  115. 115.

    Secretariat of the Carpathian Convention. Carpathian Convention; 2003.

    Google Scholar 

  116. 116.

    Scohier A, Dumont B. How do sheep affect plant communities and arthropod populations in temperate grasslands? Animal. 2012;6:1129–38.

    CAS  Article  PubMed  Google Scholar 

  117. 117.

    Jugovic J, Crne M, Pecniker ZF. The impact of grazing, overgrowth and mowing on spring butterfly (Lepidoptera: Rhopalocera) assemblages on dry karst meadows and pastures. Natura Croatica. 2013;22:157–69.

    Google Scholar 

  118. 118.

    Jerrentrup JS, Klimek S, Marchiori E, Bittante G, Ramanzin M, Sturaro E, Marini L. Impact of dairy farming on butterfly diversity in alpine summer pastures. Agric Ecosyst Environ. 2016;232:38–45.

    Article  Google Scholar 

  119. 119.

    Erhardt A. Diurnal Lepidoptera: sensitive indicators of cultivated and abandoned grassland. The Journal of Applied Ecology. 1985;22:849.

    Article  Google Scholar 

Download references


We thank Andrew Liston for counter-checking our manuscript and improving our English. Additionally, we thank the Bucegi Natural Park for permission to perform our studies.


Stefan Ehl was a member of the interdisciplinary graduate school “Cooperation of Science and Jurisprudence in Improving Development and Use of Standards for Environmental Protection – Strategies for Risk Assessment and Management,” funded by the German Science Foundation (DFG, GRK 1319). The funder did not provide any input into this study, and do not require approval of the manuscript prior to submission and publication.

Availability of data and materials

The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request. All data generated or analysed during this study are included in this published article [and its supplementary information files].

Author information




SE, NB and TS conceived and designed the experiments. SE, NB and MW performed the experiments. SE, NB and TS analyzed the data. SE, NB, LR and TS wrote the manuscript; other authors provided editorial advice. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Stefan Ehl.

Ethics declarations

Ethics approval and consent to participate

Our study had permissions from the Bucegi Natural Park. No butterflies were harmed or killed for this study.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional files

Additional file 1:

Table S1. Daily sampling effort for ungrazed and grazed sector (in minutes) calculated with the daily time in field multiplied with number of workers. Table S2. Equation, coefficient of determination (R2) and adjusted coefficient of determination (Radj2) for trend lines of the different polynomic function. (DOCX 23 kb)

Additional file 2:

Distribution of capture events of B. pales in Parcul Natural Bucegi, Romania (2014). (TIFF 4916 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Ehl, S., Böhm, N., Wörner, M. et al. Dispersal and adaptation strategies of the high mountain butterfly Boloria pales in the Romanian Carpathians. Front Zool 16, 1 (2019).

Download citation


  • Mark-release-recapture
  • Boloria pales
  • Soft protandry
  • Grazing
  • Habitat quality
  • Dispersal behaviour
  • Conservation
  • Risk spreading
  • Nectar sources