Stabilization through spatial pattern formation in metapopulations with long-range dispersal

Many studies of metapopulation models assume that spatially extended populations occupy a network of identical habitat patches, each coupled to its nearest neighbouring patches by density-independent dispersal. Much previous work has focused on the temporal stability of spatially homogeneous equilibrium states of the metapopulation, and one of the main predictions of such models is that the stability of equilibrium states in the local patches in the absence of migration determines the stability of spatially homogeneous equilibrium states of the whole metapopulation when migration is added. Here, we present classes of examples in which deviations from the usual assumptions lead to different predictions. In particular, heterogeneity in local habitat quality in combination with long-range dispersal can induce a stable equilibrium for the metapopulation dynamics, even when within-patch processes would produce very complex behaviour in each patch in the absence of migration. Thus, when spatially homogeneous equilibria become unstable, the system can often shift to a different, spatially inhomogeneous steady state. This new global equilibrium is characterized by a standing spatial wave of population abundances. Such standing spatial waves can also be observed in metapopulations consisting of identical habitat patches, i.e. without heterogeneity in patch quality, provided that dispersal is density dependent. Spatial pattern formation after destabilization of spatially homogeneous equilibrium states is well known in reaction–diffusion systems and has been observed in various ecological models. However, these models typically require the presence of at least two species, e.g. a predator and a prey. Our results imply that stabilization through spatial pattern formation can also occur in single-species models. However, the opposite effect of destabilization can also occur: if dispersal is short range, and if there is heterogeneity in patch quality, then the metapopulation dynamics can be chaotic despite the patches having stable equilibrium dynamics when isolated. We conclude that more general metapopulation models than those commonly studied are necessary to fully understand how spatial structure can affect spatial and temporal variation in population abundance.     

Do invasive species undergo metapopulation dynamics? A case study of the invasive Caspian gull,Larus cachinnans,in Poland

Aim The mechanisms of initial dispersal and habitat occupancy by invasive alien species are fundamental ecological problems. Most tests of metapopulation theory are performed on local population systems that are stable or in decline. In the current study we were interested in the usefulness of metapopulation theory to study patch occupancy, local colonization, extinction and the abundance of the invasive Caspian gull (Larus cachinnans) in its initial invasion stages. Location Waterbodies in Poland. Methods Characteristics of the habitat patches (waterbodies, 35 in total) occupied by breeding pairs of Caspian gulls and an equal sample of randomly selected unoccupied patches were compared with t‐tests. Based on presence–absence data from 1989 to 2006 we analysed factors affecting the probability of local colonization, extinction and the size of local populations using generalized linear models. Results Occupied habitat patches were significantly larger and less isolated (from other habitat patches and other local populations) and were located closer to rivers than empty patches. The proximity of local food resources (fish ponds, refuse dumps) positively affected the occurrence of breeding pairs. The probability of colonization was positively affected by patch area, and negatively by distances to fish ponds, nearest habitat patch, nearest breeding colony and to a river, and by higher forest cover around the patch boundaries. The probability of extinction was lower in patches with a higher number of breeding pairs and with a greater area of islets. The extinction probability increased with distances to other local populations, other habitat patches, fish ponds and to refuse dumps and with a higher cover of forest around the patch boundaries. The size of the local population decreased with distances to the nearest habitat patch, local population, river, fish pond and refuse dump. Local abundance was also positively affected by the area of islets in the patch. Main conclusions During the initial stages of the invasion of Caspian gulls in Poland the species underwent metapopulation‐like dynamics with frequent extinctions from colonized habitat patches. The results prove that metapopulation theory may be a useful conceptual framework for predicting which habitats are more vulnerable to invasion.     

Influence of Habitat Quality and Patch Size on Occupancy and Persistence in two Populations of the Apollo Butterfly (Parnassius apollo)

Recent studies on butterflies emphasize habitat characteristics together with metapopulation parameters (patch area and isolation) giving a more thorough understanding of processes influencing population persistence and patch occupancy, than either of them alone. We studied a coastal and an archipelago population of the Apollo butterfly (Parnassius apollo) in SW Finland. Larvae were surveyed for four years in both populations. Counting larvae on three consecutive days and temporarily removing them tested the survey accuracy. The removals showed four times higher larval abundance in the archipelago than on the coast. Survey methods were reliable, provided that empty patch status was not based on single visits only, if larval abundance was low. On the coast, large patches, and patches with high host-plant abundance were often occupied. In the archipelago, patches rich in host-plant were often occupied whereas patch area did not affect patch occupancy. In both populations, the probability of patches being occupied for three consecutive years increased with increasing host-plant abundance and patch area. Conservation of P. apollo depends on securing host-plant abundance on large enough patches in both study systems. In these systems, even crude habitat measures prove useful for understanding ecological processes behind observed patterns.     

Potential metapopulation structure and the effects of habitat quality on population size of the endangered False Ringlet butterfly

The False Ringlet (Coenonympha oedippus) is a European butterfly species, endangered due to the severe loss and fragmentation of its habitat. In Hungary, two remaining populations of the butterfly occur in lowland Purple Moorgrass meadows. We studied a metapopulation occupying twelve habitat patches in Central Hungary. Our aim was to reveal what measures of habitat quality affect population size and density of this metapopulation, estimate dispersal parameters and describe phenology of subpopulations. Local population sizes and dispersal parameters were estimated from an extensive mark–release–recapture dataset, while habitat quality was characterized by groundwater level, cover of grass tussocks, bush cover, height of vegetation and grass litter at each habitat patch. The estimated size of the metapopulation was more than 3,000 individuals. We estimated a low dispersal capacity, especially for females, indicating a very low probability of (re)colonization. Butterfly abundance and density in local populations increased with higher grass litter, lower groundwater level and larger area covered by tussocks. We suppose that these environmental factors affect butterfly abundance by determining the microclimatic conditions for both larvae and adult butterflies. Our results suggest that the long-term preservation of the studied metapopulation needs the maintenance of high quality habitat patches by appropriate mowing regime and water regulation. Management also should facilitate dispersal to strengthen metapopulation structure with creating stepping-stones or gradually increase habitat quality in present matrix.     

Does colonization asymmetry matter in metapopulations?

Despite the considerable evidence showing that dispersal between habitat patches is often asymmetric, most of the metapopulation models assume symmetric dispersal. In this paper, we develop a Monte Carlo simulation model to quantify the effect of asymmetric dispersal on metapopulation persistence. Our results suggest that metapopulation extinctions are more likely when dispersal is asymmetric. Metapopulation viability in systems with symmetric dispersal mirrors results from a mean field approximation, where the system persists if the expected per patch colonization probability exceeds the expected per patch local extinction rate. For asymmetric cases, the mean field approximation underestimates the number of patches necessary for maintaining population persistence. If we use a model assuming symmetric dispersal when dispersal is actually asymmetric, the estimation of metapopulation persistence is wrong in more than 50% of the cases. Metapopulation viability depends on patch connectivity in symmetric systems, whereas in the asymmetric case the number of patches is more important. These results have important implications for managing spatially structured populations, when asymmetric dispersal may occur. Future metapopulation models should account for asymmetric dispersal, while empirical work is needed to quantify the patterns and the consequences of asymmetric dispersal in natural metapopulations.     

Spatial ecology of host–parasitoid interactions: a threatened butterfly and its specialised parasitoid

Habitat conservation for threatened temperate insect species is often guided by one of two paradigms: a metapopulation approach focusing on patch area, isolation and number; or a habitat approach focusing on maintaining high quality habitat for the focal species. Recent research has identified the additive and interacting importance of both approaches for maintaining populations of threatened butterflies. For specialised host-parasitoid interactions, understanding the consequences of habitat characteristics for the interacting species is important, because (1) specialised parasitoids are particularly vulnerable to the consequences of fragmentation, and (2) altered interaction frequencies resulting from changes to habitat management or the spatial configuration of habitat are likely to have consequences for host dynamics. The spatial ecology of Cotesia bignellii, a specialist parasitoid of the threatened butterfly Euphydryas aurinia, was investigated at two spatial scales: within habitat patches (at the scale of individual aggregations of larvae, or ‘webs’) and among habitat patches (the scale of local populations). Parasitism rates were investigated in relation to larval web size, vegetation sward height and host density. Within patches, the probability of a larval webs being parasitized increased significantly with increasing number of larvae in the web, and parasitism rates increased significantly with increasing web isolation. The proportion of webs parasitized was significantly and negatively correlated with cluster density. Among habitat patches the proportion of parasitized webs decreased as cluster density increased. Clusters with a high proportion of larval webs parasitized tended to have lower parasitism rates per larval web. These results support the call for relatively large and continuous habitat patches to maintain stable parasitoid and host populations. Conservation efforts directed towards maintenance of high host plant density could allow E. aurinia to reduce parasitism risk, while providing C. bignellii with sufficient larval webs to allow population persistence.     

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 for spatially extended predator–prey interactions

Traditional metapopulation theory classifies a metapopulation as a spatially homogeneous population that persists on neighboring habitat patches. The fate of each population on a habitat patch is a function of a balance between births and deaths via establishment of new populations through migration to neighboring patches. In this study, we expand upon traditional metapopulation models by incorporating spatial heterogeneity into a previously studied two-patch nonlinear ordinary differential equation metapopulation model, in which the growth of a general prey species is logistic and growth of a general predator species displays a Holling type II functional response. The model described in this work assumes that migration by generalist predator and prey populations between habitat patches occurs via a migratory corridor. Thus, persistence of species is a function of local population dynamics and migration between spatially heterogeneous habitat patches. Numerical results generated by our model demonstrate that population densities exhibit periodic plane-wave phenomena, which appear to be functions of differences in migration rates between generalist predator and prey populations. We compare results generated from our model to results generated by similar, but less ecologically realistic work, and to observed population dynamics in natural metapopulations.     

Metapopulation Dynamics of the Mistletoe and Its Host in Savanna Areas with Different Fire Occurrence

Mistletoes are aerial hemiparasitic plants which occupy patches of favorable habitat (host trees) surrounded by unfavorable habitat and may be possibly modeled as a metapopulation. A metapopulation is defined as a subdivided population that persists due to the balance between colonization and extinction in discrete habitat patches. Our aim was to evaluate the dynamics of the mistletoe Psittacanthus robustus and its host Vochysia thyrsoidea in three Brazilian savanna areas using a metapopulation approach. We also evaluated how the differences in terms of fire occurrence affected the dynamic of those populations (two areas burned during the study and one was fire protected). We monitored the populations at six-month intervals. P. robustus population structure and dynamics met the expected criteria for a metapopulation: i) the suitable habitats for the mistletoe occur in discrete patches; (ii) local populations went extinct during the study and (iii) colonization of previously non-occupied patches occurred. The ratio of occupied patches decreased in all areas with time. Local mistletoe populations went extinct due to two different causes: patch extinction in area with no fire and fire killing in the burned areas. In a burned area, the largest decrease of occupied patch ratios occurred due to a fire event that killed the parasites without, however, killing the host trees. The greatest mortality of V. thyrsoidea occurred in the area without fire. In this area, all the dead trees supported mistletoe individuals and no mortality was observed for parasite-free trees. Because P. robustus is a fire sensitive species and V. thyrsoidea is fire tolerant, P. robustus seems to increase host mortality, but its effect is lessened by periodic burning that reduces the parasite loads.     

Biased correlated random walk and foray loop: which movement hypothesis drives a butterfly metapopulation?

Animals in fragmented landscapes have a major challenge to move between high-quality habitat patches through lower-quality matrix. Two current mechanistic hypotheses that describe the movement used by animals outside of their preferred patches (e.g., high-quality habitat or home range) are the biased, correlated random walk (BCRW) and the foray loop (FL). There is also a variant of FL with directed movement (FL_dm). While these have been most extensively tested on butterflies, they have never been tested simultaneously with data across a whole metapopulation and over multiple generations, two key scales for population dynamics. Using the pattern-oriented approach, we compare support for these competing hypotheses with a spatially explicit individual-based simulation model on an 11-year dataset that follows 12 patches of the federally endangered Fender’s blue butterfly (Plebejus icarioides fenderi) in Oregon’s Willamette Valley. BCRW and medium-scale FL and FL_dm scenarios predicted the annual total metapopulation size for ≥9 of 12 patches as well as patch extinctions. The key difference, however, was that the FL scenarios predicted patch colonizations and persistence poorly, failing to adequately capture movement dynamics; BCRW and one FL_dm scenario predicted the observed patch colonization and persistence with reasonable probabilities. This one FL_dm scenario, however, had larger prediction intervals. BCRW, the biologically simplest and thus most parsimonious movement hypothesis, performed consistently well across all nine different tests, resulting in the highest quality metapopulation predictions for butterfly conservation.     

Drivers of local densities of endangered Lycaena helle butterflies in a fragmented landscape

Due to their specialised habitat requirements, butterflies are particularly vulnerable to habitat loss and fragmentation. Understanding the drivers of local abundances of species is essential for their effective conservation in fragmented landscapes. We investigated factors affecting population densities of an endangered European butterfly, the Violet Copper (Lycaena helle), occurring in a small metapopulation near the city of Kraków, southern Poland. The environmental parameters tested as predictors of the local densities of the species included both the variables associated with spatial structure of habitats such as patch sizes, their isolation and fragmentation as well as those potentially reflecting habitat quality. Patch area and vegetation height turned out to be the only factors significantly influencing L. helle densities, both having a positive effect. The positive impact of patch area is a bit surprising, since its relationship with population densities is typically negative in butterflies. In our study system it is likely to derive from source-sink dynamics as the smaller habitat patches are apparently too small to sustain viable local populations. In turn, the positive influence of vegetation height implies that the ongoing succession does not deteriorate the quality of the recently abandoned meadows yet, whereas higher turf may provide better sheltering places. The loss of almost half of L. helle habitat patches in the study area in recent years is alarming. However, its inclusion into the Natura 2000 system should help to conserve the species as long as this act is followed by proper management of its habitats.     

Population turnover and habitat dynamics in Notonecta (Hemiptera: Notonectidae) metapopulations

Simple metapopulation models assume that local populations occur in patches of uniform quality habitat separated by non-habitat. However field metapopulations tend to show considerable spatial and temporal variation in patch quality, and hence probability of occupancy. This may have implications for the adequacy of simple metapopulation models in describing and predicting regional population dynamics of natural systems. This study investigated the effects of habitat characteristics on landscape-scale occupancy dynamics of two species of backswimmer (Notonecta, Hemiptera: Notonectidae) in small freshwater ponds. The results demonstrated clear links between habitat, pond occupancy and population turnover, particularly local extinction. There were considerable changes in the habitat of individual ponds between years, but local changes were not spatially correlated and the frequency distribution of habitat conditions at the landscape level remained similar in different years. Stable occupancy levels of Notonecta species appears to result from a balance of the rates of creation and loss of suitable habitat due to spatially uncorrelated habitat change. Systems such as this, where turnover is driven by habitat dynamics, demonstrate the potential value of incorporating the dynamics of habitat change into metapopulation models. Such developments are likely to improve predictions of landscape-scale occupancy dynamics, whilst also allowing patch-level predictions of occupancy, based on local habitat conditions. Received: 18 August 1999 / Accepted: 3 December 1999     

Spatial Heterogeneity in Habitat Quality and Cross-Scale Interactions in Metapopulations

Abstract Integration of habitat heterogeneity into spatially realistic metapopulation approaches reveals the potential for key cross-scale interactions. Broad-scale environmental gradients and land-use practices can create autocorrelation of habitat quality of suitable patches at intermediate spatial scales. Patch occupancy then depends not only on habitat quality at the patch scale but also on feedbacks from surrounding neighborhoods of autocorrelated patches. Metapopulation dynamics emerge from how demographic and dispersal processes interact with relevant habitat heterogeneity. We provide an empirical example from a metapopulation of round-tailed muskrats (Neofiber alleni) in which habitat quality of suitable patches was spatially autocorrelated most strongly within 1,000 m, which was within the expected dispersal range of the species. After controlling for factors typically considered in metapopulation studies—patch size, local patch quality, patch connectivity—we use a cross-variogram analysis to demonstrate that patch occupancy by muskrats was correlated with habitat quality across scales ≤1,171 m. We also discuss general consequences of spatial heterogeneity of habitat quality for metapopulations related to potential cross-scale interactions. We focus on spatially correlated extinctions and metapopulation persistence, hierarchical scaling of source–sink dynamics, and dispersal decisions by individuals in relation to information constraints.     

Metapopulation dynamics with quasi-local competition

Stepping-stone models for the ecological dynamics of metapopulations are often used to address general questions about the effects of spatial structure on the nature and complexity of population fluctuations. Such models describe an ensemble of local and spatially isolated habitat patches that are connected through dispersal. Reproduction and hence the dynamics in a given local population depend on the density of that local population, and a fraction of every local population disperses to neighboring patches. In such models, interesting dynamic phenomena, e.g. the persistence of locally unstable predator-prey interactions, are only observed if the local dynamics in an isolated patch exhibit non-equilibrium behavior. Therefore, the scope of these models is limited. Here we extend these models by making the biologically plausible assumption that reproductive success in a given local habitat not only depends on the density of the local population living in that habitat, but also on the densities of neighboring local populations. This would occur if competition for resources occurs between neighboring populations, e.g. due to foraging in neighboring habitats. With this assumption of quasi-local competition the dynamics of the model change completely. The main difference is that even if the dynamics of the local populations have a stable equilibrium in isolation, the spatially uniform equilibrium in which all local populations are at their carrying capacity becomes unstable if the strength of quasi-local competition reaches a critical level, which can be calculated analytically. In this case the metapopulation reaches a new stable state, which is, however, not spatially uniform anymore and instead results in an irregular spatial pattern of local population abundance. For large metapopulations, a huge number of different, spatially non-uniform equilibrium states coexist as attractors of the metapopulation dynamics, so that the final state of the system depends critically on the initial conditions. The existence of a large number of attractors has important consequences when environmental noise is introduced into the model. Then the metapopulation performs a random walk in the space of all attractors. This leads to large and complicated population fluctuations whose power spectrum obeys a red-shifted power law. Our theory reiterates the potential importance of spatial structure for ecological processes and proposes new mechanisms for the emergence of non-uniform spatial patterns of abundance and for the persistence of complicated temporal population fluctuations.     

The effective size of a metapopulation living in a heterogeneous patch network

I analyze stochastic patch occupancy models (SPOMs), which record habitat patches as empty or occupied. A problem with SPOMs has been that if the spatial structure of a heterogeneous habitat patch network is taken into account, the computational effort needed to analyze a SPOM grows as a power of 2n, where n is the number of habitat patches. I propose a computationally feasible approximation method, which approximates the behavior of a heterogeneous SPOM by an "ideal" metapopulation inhabiting a network of identical and equally connected habitat patches. The transformation to the ideal metapopulation is based on weighting the individual patch occupancies by the dynamic values of the habitat patches, which may be calculated from the deterministic mean-field approximation of the original SPOM. Conceptually, the method resembles the calculation of the effective size of a population in the context of population genetics. I demonstrate how the method may be applied to SPOMs with flexible structural assumptions and with spatially correlated and temporally varying parameter values. I apply the method to a real habitat patch network inhabited by the Glanville fritillary butterfly, illustrating that the metapopulation dynamics of this species are essentially driven by temporal variability in the environmental conditions.     

Three ways of assessing metapopulation structure in the butterfly Plebejus argus

1. Three independent methods were used to investigate population structure in the butterfly Plebejus argus . First, migration and dispersal ability were measured by mark–release–recapture in seven adjacent habitat patches, and by release of butterflies in unoccupied habitat. Secondly, colonization of newly created habitat was observed over 7 years. Finally, genetic differentiation of local populations within a metapopulation was investigated. Sampled local populations included parts of the mark–release–recapture study area.
2. Plebejus argus is relatively sedentary: the maximum movement detected was 395 m, and only 2% of individuals moved further than 100 m between recaptures on different days. None the less, adjacent local populations in the mark–release–recapture study area were linked by occasional migration, with ≈ 1.4% of individuals moving between patches separated by 13–200 m.
3. Despite low mobility, observed colonizations occurred rapidly over distances of 1 km. Because P. argus occurs at high population densities, 1.4% migration can generate enough migrants to colonize newly suitable habitat quickly at this spatial scale.
4. Mark–release–recapture data were used to predict that there would be limited genetic differentiation through drift between local populations at this spatial scale. The prediction was supported by allele frequency data for the same local populations.
5. Genetic differentiation often indicates higher levels of migration than are revealed by the movements of marked individuals. This study shows that when experimental releases and extensive marking are undertaken in areas that are large relative to most movements, indirect measures of gene flow and direct measures of dispersal can concur.
6. Evidence from the three different approaches was complementary, indicating that P. argus occurs as metapopulations within the study area.     

荒漠破碎化生境中长爪沙鼠集合种群野外验证研究

近年来,人类活动和自然干扰,导致内蒙古阿拉善荒漠区生境的破碎化,出现了长爪沙鼠在不同斑块间的不连续分布,每一斑块内可能存在一个局域种群,而集合种群建立的前提条件,是局域种群斑块状分布在离散的栖息地环境中。2002~2012年每年的4~10月,在阿拉善荒漠区禁牧、轮牧、过牧和开垦4种人为不同利用方式形成的生境斑块中,采用标志重捕法对长爪沙鼠（Meriones unguiculatus）种群进行定点监测。通过分析长爪沙鼠种群动态,计算各局域种群的灭绝风险,利用Spearman秩相关系数检验种群动态的空间同步性,同时以种群周转率对长爪沙鼠扩散能力进行评估,以检验阿拉善荒漠区长爪沙鼠种群空间结构是否具有经典集合种群的功能。结果表明：（1） 不同生境斑块可被长爪沙鼠局域种群占据,11年间捕获长爪沙鼠2~7次不等;（2） 长爪沙鼠所有局域种群均具有灭绝风险,在轮牧区和禁牧区灭绝率高达1.000 0,开垦区灭绝率最低,也达到0.333 4,而本研究期间最大局域种群（2008年过牧区,26只/hm²）,在2010年发生了局域灭绝;（3） 不同生境斑块间没有明显的空间隔离而阻碍局域种群的重新建立,长爪沙鼠扩散能力较强,绝大部分月份的种群周转率在50.0%以上,特别是周转率达到100.0%的月份较多;（4） 不同生境斑块间仅轮牧区和禁牧区中长爪沙鼠种群密度显著正相关（P<0.05）,而其他生境斑块间相关性均不显著（P >0.05）,长爪沙鼠局域种群整体显示出明显的非同步空间动态。阿拉善荒漠区长爪沙鼠种群满足作为经典集合种群物种区域续存的4个条件,具有作为研究小哺乳动物集合种群的潜在价值。     

Contrasting movement patterns in two species of chequerspot butterflies, Euphydryas aurinia and Melitaea phoebe, in the same patch network

Abstract.  1. Mark–release–recapture studies were conducted on two species of chequerspot butterfly, Euphydryas aurinia and Melitaea phoebe , in the same habitat patch network in Yanjiaping, a small basin in the Taihang Mountains, north-west of Beijing, China, in 2000.
2. Euphydryas aurinia tended to stay in the habitat patches and to move to neighbouring patches, whereas M. phoebe moved widely among patches in the entire network.
3. The parameters of the virtual migration model showed higher daily emigration propensity in M. phoebe and in E. aurinia males than in E. aurinia females, and significantly greater average daily movement distance in M. phoebe than in E. aurinia .
4. The results are consistent with the previous findings showing genetic structuring among local populations of E. aurinia but not among local populations of M. phoebe .
5. Based on the genetic and ecological results, it was concluded that E. aurinia has a classic metapopulation in the study area, whereas M. phoebe appears to have a source–sink metapopulation.
6. In 2000, when there was an overall increase in the abundance of the two species, the limited mobility of E. aurinia resulted in an increase in the average local population size, whereas the increase in the number of local populations in M. phoebe was due to its high mobility .     

The fragmented population structure of the Sardinian chalk hill blue butterfly (Lepidoptera, Lycaenidae)

The endemic Sardinian chalk hill blue butterfly, Polyommatus coridon gennargenti, is considered vulnerable to extinction because of its low genetic variation and restricted distribution. The species also has a fragmented distribution, which follows the patchy distribution pattern of its larval host-plant. A preliminary investigation of the population structure of P. coridon gennargenti was carried out on a small network of four local populations by means of capture–recapture methods. Estimated population sizes and movement rates among the four adjacent local populations suggest that this taxon has a metapopulation structure composed of loosely connected small local populations. Natural fragmentation, isolation, and traditional land use contribute to the vulnerability of P. coridon gennargenti to extinction. Low effective population sizes and restricted movement between habitat patches lead to inbreeding and an increased vulnerability to extinction of this island population.     

Long-term persistence of species and the SLOSS problem

The single large or several small (SLOSS) problem has been addressed in a large number of empirical and theoretical studies, but no coherent conclusion has yet been reached. Here I study the SLOSS problem in the context of metapopulation dynamics. I assume that there is a fixed total amount A(0) of habitat available, and I derive formulas for the optimal number n and area A of habitat patches, where n=A(0)/A. I consider optimality in two ways. First, I attempt to maximize the time to metapopulation extinction, which is a relevant measure for metapopulation viability for rare and threatened species. Second, I attempt to maximize the metapopulation capacity of the habitat patch network, which corresponds both with maximizing the distance to the deterministic extinction threshold and with maximizing the fraction of occupied patches. I show that in the typical case, a small number of large patches maximizes the metapopulation capacity, while an intermediate number of habitat patches maximizes the time to extinction. The main conclusion stemming from the analysis is that the optimal number of patches is largely affected by the relationship between habitat patch area and rates of immigration, emigration and local extinction. Here this relationship is summarized by a single factor zeta, termed the patch area scaling factor.