Metapopulation dynamics of the bog fritillary butterfly: movements between habitat patches

Recent studies on metapopulation dynamics have emphasized the need for improved methods for quantifying individual movements between local populations and habitat patches. In this paper, we report on a 6-yr study in which a network of 12 habitat patches occupied by the bog fritillary, Proclossiana eunomia , was surveyed, with special focus on quantifying movements between the habitat patches. We applied the Virtual Migration model which has been designed to estimate survival and migration parameters in a metapopulation of several connected local populations. The model was parameterized using mark-release-recapture data collected during 6 yr. Generally, the estimated parameter values indicated a high level of movements, with roughly half of butterfly-days spent outside the natal patch. Mortality within patches was higher in males than in females. Females tended to be more mobile and spent more time outside their natal patch than males. Further analysis of the MRR data shows that in this protandrous species males tend to move very little between habitat patches before substantial numbers of females have emerged.     

Metapopulation dynamics of the bog fritillary butterfly: modelling the effect of habitat fragmentation

Population viability analysis (PVA) and metapopulation theory are valuable tools to model the dynamics of spatially structured populations. In this article we used a spatially realistic population dynamic model to simulate the trajectory of a Proclossiana eunomia metapopulation in a network of habitat patches located in the Belgian Ardenne. Sensitivity analysis was used to evaluate the relative influence of the different parameters on the model output. We simulated habitat loss by removing a percentage of the original habitat, proportionally in each habitat patch. Additionally, we evaluated isolation and fragmentation effects by removing and dividing habitat patches from the network, respectively. The model predicted a slow decline of the metapopulation size and occupancy. Extinction risks predicted by the model were highly sensitive to environmental stochasticity and carrying capacity. For a determined level of habitat destruction, the expected lifetime of the metapopulation was highly dependent on the spatial configuration of the landscape. Moreover, when the proportion of removed habitat is above 40% of the original habitat, the loss of whole patches invariably leads to the strongest reduction in metapopulation viability.     

Metapopulation dynamics of the bog fritillary butterfly: demographic processes in a patchy population

Metapopulation models predicting the persistence of species on the verge of extinction in fragmented landscapes have to include local population dynamics. In this paper, we report on a 9-yr study in which a patchy population of the bog fritillary was monitored each year, with special focus on quantifying demographic parameters. We parameterised various constrained linear models (CLM) with capture-mark-recapture (CMR) data collected in the field each year and we tested the effects of several variables on demographic parameters. Selecting the best CLM, we estimated the mean and the variance of (1) individual parameters (survival and catchability) and (2) population features (population size, sex-ratio and recruitment). Survival was not related to the age of the individual but decreased continuously during the flight period, which means that lifetime expectancy is conversely related to the delay between the moment of the emergence of a given individual as compared to the emergence of the first adult in the population. Survival was lower in males than in females, as a consequence of male mate-locating behaviour. Daily population size fits a parabolic distribution in both sexes, males appearing, reaching their abundance peak and dying before females. Yearly variation in the sex-ratio is the consequence of the between-sex difference in the recruitment curve. Total adult population size shows large variations between years. Analysis of the relation between growth rate and abundance indicates density dependence, probably related to parasitism of the larval stages. Density at the equilibrium is about 200 adults/ha.     

Metapopulation viability analysis of the bog fritillary butterfly using RAMAS/GIS

Many species living in man-shaped landscapes are restricted to small natural habitat patches and form metapopulations; predicting their future is a central issue in applied ecology. We examined the viability of the bog fritillary butterfly Proclossiana eunomia Esper, a specialist glacial relict species, in a highly fragmented landscape (<1% of suitable habitat in 10 km²), by way of population viability analysis. We used comprehensive data from a long-term study in which a patchy population was monitored during ten consecutive years to parameterise a spatially structured metapopulation model using commercially available platform RAMAS/GIS 3.0. Population growth rate was density-dependent and modulated by various climatic variables acting on different developmental stages of the butterfly. Density dependence was probably related to larval parasitism by a specific parasitoid. Population size was negatively affected by an increase in the mean temperature. Dispersal was modelled as the observed proportion of movements between patches, taking into account the probability of emigration out of a given patch. Our model provided results close to the picture of the system drawn from the field data and was considered as useful in making predictions about the metapopulation. Demographic parameters proved to have a far higher impact on metapopulation persistence than dispersal or correlation of local dynamics. Scenarios simulating both global warming and management of habitat patches by rustic herbivore grazing indicated a decrease in the viability of the metapopulation. Our results prompted the regional nature conservation agency to modify the planned management regime. We urge conservation biologists to use structured population models including local population dynamics for viability analysis targeted to such threatened metapopulations in highly fragmented landscapes.     

Metapopulation structure and habitat quality in modelling dispersal in the butterfly Iolana iolas

Quantifying dispersal is fundamental to understanding the effects of fragmentation on populations. Although it has been shown that patch and matrix quality can affect dispersal patterns, standard metapopulation models are usually based on the two basic variables, patch area and connectivity. In 2004 we studied migration patterns among 18 habitat patches in central Spain for the butterfly Iolana iolas, using mark–release–recapture methods. We applied the virtual migration (VM) model and estimated the parameters of emigration, immigration and mortality separately for males and females. During parameter estimation and model simulations, we used original and modified patch areas accounting for habitat quality with three different indices. Two indices were based on adult and larval resources (flowers and fruits) and the other one on butterfly density. Based on unmodified areas, our results showed that both sexes were markedly different in their movements and mortality rates. Females emigrated more frequently from patches, but males that emigrated were estimated to move longer daily dispersal distances and suffer higher mortality than females during migration. Males were more likely to emigrate from small than from large patches, but patch area had no significant effect on female emigration. The effects of area on immigration rate and the within-patch mortality were similar in both sexes. Based on modified areas, the estimated parameter values and the model simulation results were similar to those estimated using the unmodified patch areas. One possible reason for the failure to significantly improve the parameter estimates of the VM model is the fact that resource quantity and butterfly population sizes were strongly correlated with patch area. Our results suggest that the standard VM modelling approach, based on patch area and connectivity, can provide a realistic picture of the movement patterns of I. iolas .     

Spatio‐temporal dynamics of density‐dependent dispersal during a population colonisation

Predicting population colonisations requires understanding how spatio‐temporal changes in density affect dispersal. Density can inform on fitness prospects, acting as a cue for either habitat quality, or competition over resources. However, when escaping competition, high local density should only increase emigration if lower‐density patches are available elsewhere. Few empirical studies on dispersal have considered the effects of density at the local and landscape scale simultaneously. To explore this, we analyze 5 years of individual‐based data from an experimental introduction of wild guppies Poecilia reticulata. Natal dispersal showed a decrease in local density dependence as density at the landscape level increased. Landscape density did not affect dispersal among adults, but local density‐dependent dispersal switched from negative (conspecific attraction) to positive (conspecific avoidance), as the colonisation progressed. This study demonstrates that densities at various scales interact to determine dispersal, and suggests that dispersal trade‐offs differ across life stages.     

Gene flow rise with habitat fragmentation in the bog fritillary butterfly (Lepidoptera: Nymphalidae)

Background  

The main components of the spatial genetic structure of the populations are neighbourhood size and isolation by distance. These may be inferred from the allele frequencies across a series of populations within a region. Here, the spatial population structure of Proclossiana eunomia was investigated in two mountainous areas of southern Europe (Asturias, Spain and Pyrenees, France) and in two areas of intermediate elevation (Morvan, France and Ardennes, Belgium).     

Natal habitat effects drive density—dependent scaling of dispersal decisions

The dispersal behavior of a species is critical for the stability and persistence of its populations across a landscape. How population density affects dispersal decisions is important for predicting these dynamics, as the form of density‐dependent dispersal influences the stability and persistence of populations. Natal habitat experience often has strong impacts on individual dispersal behavior as well, but its influence on density‐dependent dispersal behaviors remains unexplored. Here we address this conceptual gap in two experiments separately examining habitat selection and emigration from recently colonized patches for two species of flour beetle Tribolium sp. We found that interactions between the quality of habitat experienced during natal development and current habitat for dispersal capable adults can strongly affect the form of density dependence, including reversing the direction of nonlinearities (accelerating to decelerating), or even negating the influence of population density for individual dispersal decisions. Across heterogeneous landscapes, where individuals from different populations may experience different natal habitats, this altering of density‐dependent relationships is predicted by theory to fundamentally influence regional population dynamics. Our results indicate that species which occur across heterogeneous environments, such as during conservation reintroductions, or as invasive species spread, have much potential for natal experience to interact with density dependence and influence local and regional population dynamics.     

Temporal scales of density‐dependent habitat selection in a large grazing herbivore

Habitat selection is a density‐dependent process, but little is known regarding how this relationship may vary across different temporal scales. Over long time scales, grazing shapes the structure, diversity and functioning of terrestrial ecosystems, and grazing‐induced changes in forage production over time are likely to affect the level of density dependence in habitat selection. In this fully‐replicated, landscape‐scale experiment, we investigated how density‐dependent habitat selection by a large grazing herbivore, sheep Ovis aries, develops over the time scale of a decade. We also address an often‐neglected challenge in habitat selection studies; namely, whether there is variation in use within a particular habitat or vegetation type and why. We found clear evidence of density dependence in habitat selection, with a wider use of habitats at high density. Despite a change in the standing biomass of high‐productivity vegetation at high herbivore density over the years, with herb biomass declining and graminoid biomass increasing, there was no clear evidence that these grazing‐induced changes in habitat over the years were strong enough to affect the level of density‐dependent habitat selection. The difference in selection for high versus low‐productivity habitats remained similar, despite annual fluctuations in the strength of selection. We found strong variation in selection within each vegetation type, even when vegetation types were mapped at a fine‐resolution scale. Our study shows that despite the interactive effects of herbivores and habitats, they are not always sufficiently strong enough to affect the level of density‐dependent habitat selection.     

Evolution of density‐dependent movement during experimental range expansions

Range expansions and biological invasions are prime examples of transient processes that are likely impacted by rapid evolutionary changes. As a spatial process, range expansions are driven by dispersal and movement behaviour. Although it is widely accepted that dispersal and movement may be context‐dependent, for instance density‐dependent, and best represented by reaction norms, the evolution of density‐dependent movement during range expansions has received little experimental attention. We therefore tested current theory predicting the evolution of increased movement at low densities at range margins using highly replicated and controlled range expansion experiments across multiple genotypes of the protist model system Tetrahymena thermophila. Although rare, we found evolutionary changes during range expansions even in the absence of initial standing genetic variation. Range expansions led to the evolution of negatively density‐dependent movement at range margins. In addition, we report the evolution of increased intrastrain competitive ability and concurrently decreased population growth rates in range cores. Our findings highlight the importance of understanding movement and dispersal as evolving reaction norms and plastic life‐history traits of central relevance for range expansions, biological invasions and the dynamics of spatially structured systems in general.     

Ontogenetic variation in density‐dependent growth of brown trout through habitat competition

1. Density‐dependent growth has been widely reported in freshwater fishes, but the ontogenetic evolution of competition and its subsequent effects on growth through a life span remains unclear. 2. Patterns of competition can be described by integrating population abundance data with habitat‐modelling results. Weighted usable area (WUA; m² WUA ha^?1) curves are obtained for each flow value and are then coupled with demographic data to obtain the occupancy rates (trout m^?2 WUA, the density of a given age class related to its suitable habitat) of the WUA for every age class, year and site. 3. We examined a long‐term data series searching for temporal variation in the influence of habitat occupancy rate on the growth of brown trout Salmo trutta. We tested whether (i) mean cohort mass (mean mass of the cohort during the first 3 years of life) is affected by the occupancy rate experienced across a life span; and (ii) the occupancy rate experienced at different ages influenced mean body size. 4. We observed a consistent negative power relationship between average cohort mass and mean occupancy rate through a life span, indicating that stronger cohorts were related to a reduced growth, with likely consequences for individual fitness. 5. The effects of occupancy rate on size‐at‐age were mainly detected in the size attained at the second year of life, but they were because of the competition at different times. Thus, the level of competition varied through ontogeny, in some of the rivers affecting growth since the first year of life, whereas in most of the rivers the main effects on body size resulted from the competition during the second year of life. 6. Occupancy rate appears more appropriate than density for assessing the occurrence of habitat competition in freshwater fishes, since it encompasses the differences in quantity and quality of suitable habitat for each age class. 7. Our study highlights the importance of density‐dependent growth as a key process in the dynamics of brown trout populations, its temporal variation depending on the temporal changes of density and the variation of competition associated with the habitat capacity for each life stage.     

Modelling single nucleotide effects in phosphoglucose isomerase on dispersal in the Glanville fritillary butterfly: coupling of ecological and evolutionary dynamics

Dispersal comprises a complex life-history syndrome that influences the demographic dynamics of especially those species that live in fragmented landscapes, the structure of which may in turn be expected to impose selection on dispersal. We have constructed an individual-based evolutionary sexual model of dispersal for species occurring as metapopulations in habitat patch networks. The model assumes correlated random walk dispersal with edge-mediated behaviour (habitat selection) and spatially correlated stochastic local dynamics. The model is parametrized with extensive data for the Glanville fritillary butterfly. Based on empirical results for a single nucleotide polymorphism (SNP) in the phosphoglucose isomerase (Pgi) gene, we assume that dispersal rate in the landscape matrix, fecundity and survival are affected by a locus with two alleles, A and C, individuals with the C allele being more mobile. The model was successfully tested with two independent empirical datasets on spatial variation in Pgi allele frequency. First, at the level of local populations, the frequency of the C allele is the highest in newly established isolated populations and the lowest in old isolated populations. Second, at the level of sub-networks with dissimilar numbers and connectivities of patches, the frequency of C increases with decreasing network size and hence with decreasing average metapopulation size. The frequency of C is the highest in landscapes where local extinction risk is high and where there are abundant opportunities to establish new populations. Our results indicate that the strength of the coupling of the ecological and evolutionary dynamics depends on the spatial scale and is asymmetric, demographic dynamics having a greater immediate impact on genetic dynamics than vice versa.     

Condition‐dependent movement and dispersal in experimental metacommunities

Dispersal and the underlying movement behaviour are processes of pivotal importance for understanding and predicting metapopulation and metacommunity dynamics. Generally, dispersal decisions are condition‐dependent and rely on information in the broad sense, like the presence of conspecifics. However, studies on metacommunities that include interspecific interactions generally disregard condition‐dependence. Therefore, it remains unclear whether and how dispersal in metacommunities is condition‐dependent and whether rules derived from single‐species contexts can be scaled up to (meta)communities. Using experimental protist metacommunities, we show how dispersal and movement depend on and are adjusted by the strength of interspecific interactions. We found that the predicting movement and dispersal in metacommunities requires knowledge on behavioural responses to intra‐ and interspecific interaction strengths. Consequently, metacommunity dynamics inferred directly from single‐species metapopulations without taking interspecific interactions into account are likely flawed. Our work identifies the significance of condition‐dependence for understanding metacommunity dynamics, stability and the coexistence and distribution of species.     

Bayesian state-space modeling of metapopulation dynamics in the Glanville fritillary butterfly

The complexity of mathematical models of ecological dynamics varies greatly, and it is often difficult to judge what would be the optimal level of complexity in a particular case. Here we compare the parameter estimates, model fits, and predictive abilities of two models of metapopulation dynamics: a detailed individual-based model (IBM) and a population-based stochastic patch occupancy model (SPOM) derived from the IBM. The two models were fitted to a 17-year time series of data for the Glanville fritillary butterfly (Melitaea cinxia) inhabiting a network of 72 small meadows. The data consisted of biannual counts of larval groups (IBM) and the annual presence or absence of local populations (SPOM). The models were fitted using a Bayesian state-space approach with a hierarchical random effect structure to account for observational, demographic, and environmental stochasticities. The detection probability of larval groups (IBM) and the probability of false zeros of local populations (SPOM) in the observation models were simultaneously estimated from the time-series data and independent control data. Prior distributions for dispersal parameters were obtained from a separate analysis of mark–recapture data. Both models fitted the data about equally, but the results were more precise for the IBM than for the SPOM. The two models yielded similar estimates for a random effect parameter describing habitat quality in each patch, which were correlated with independent empirical measures of habitat quality. The modeling results showed that variation in habitat quality influenced patch occupancy more through the effects on movement behavior at patch edges than on carrying capacity, whereas the latter influenced the mean population size in occupied patches. The IBM and the SPOM explained 63% and 45%, respectively, of the observed variation in the fraction of occupied habitat area among 75 independent patch networks not used in parameter estimation. We conclude that, while carefully constructed, detailed models can have better predictive ability than simple models, this advantage comes with the cost of greatly increased data requirements and computational challenges. Our results illustrate how complex models can be helpful in facilitating the construction of effective simpler models.     

An experimental examination of the effects of habitat quality on the dispersal and local abundance of the butterfly Parnassius smintheus

Abstract 1. Nectar flower abundance was manipulated through flower removal, and sex ratio was manipulated by moving individual butterflies within a series of nine alpine meadows. The movement and abundance of the butterfly Parnassius smintheus in the meadows were monitored using mark–release–recapture methods.
2. A total of 937 butterflies, 698 males and 239 females, was captured. There were 223 observed between-meadow movements. Fifty-two per cent of males and 35% of females moved among meadows.
3. The immigration of male butterflies was related positively to nectar flowers, host plant abundance, and female butterflies. Male emigration was not affected by any of the treatments. The number of males captured was related positively to nectar flowers and host plants but not affected by sex ratio. The number of resident male butterflies was greater in meadows containing flowers and was related positively to host plant abundance, but unaffected by sex ratio.
4. Flower removal, sex ratio, and abundance of Sedum had no significant effect on the abundance, movement, or residence time for female butterflies, in part due to small sample size.
5. The fact that males immigrate to higher quality meadows suggests that male butterflies are assessing meadow quality, either by sampling meadows or potentially from a distance using olfactory cues.     

Long distance dispersal and landscape occupancy in a metapopulation of the cranberry fritillary butterfly

Movements between habitat patches in a patchy population of the butterfly Boloria aquilonaris were monitored using capture-mark-recapture methods during three successive generations. For each data set, the inverse cumulative proportion of individuals moving 100 m distance classes was fitted to the negative exponential function and the inverse power function. In each case, the negative exponential function provided a better fit than the inverse power function. Two dispersal kernels were generated using both negative exponential and inverse power functions. These dispersal kernels were used to predict movements between 14 suitable sites in a landscape of 220 km². The negative exponential function generated a dispersal kernel predicting extremely low probabilities for movements exceeding 1 km. The inverse power function generated probabilities predicting that between site movements were possible, according to metapopulation size. CMR studies in the landscape revealed that long distance movements occurred at each generation, corresponding to predictions of the inverse power function dispersal kernel. A total of 26 movements between sites (up to 13 km) were detected, together with recolonisation of empty sites. The spatial scale of the metapopulation dynamics is larger than ever reported on butterflies and long distance movements clearly matter to the persistence of this species in a highly fragmented landscape.     

Heritability of dispersal rate and other life history traits in the Glanville fritillary butterfly

Knowing the variances and heritabilities (h(2)) of life history traits in populations living under natural conditions is necessary for a mechanistic understanding of respective evolutionary processes. I estimated heritabilities of several life history traits, including dispersal rate, body mass, age at first reproduction, egg mass, clutch size and lifetime reproductive success, in the Glanville fritillary butterfly (Melitaea cinxia) using parent-offspring regression. Experiments were conducted under field conditions in a large population cage (32 x 26 m). Heritability estimates ranged from zero to almost one and several were significantly different from zero. Body size for both sexes, female age at first reproduction and egg weight were all moderately to highly heritable, whereas heritabilities were low or non-existent in clutch size and lifetime egg production. Heritability estimates for dispersal rate varied between the sexes, so that dispersal was heritable from mother to her female offspring only. This finding is consistent with previous results showing that the F1 female but not male offspring of females that naturally established new populations in the field are significantly more dispersive than butterflies in old populations.