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.     

Dispersal and the metapopulation paradigm in amphibian ecology and conservation: are all amphibian populations metapopulations?

Amphibians are frequently characterized as having limited dispersal abilities, strong site fidelity and spatially disjunct breeding habitat. As such, pond‐breeding species are often alleged to form metapopulations. Amphibian species worldwide appear to be suffering population level declines caused, at least in part, by the degradation and fragmentation of habitat and the intervening areas between habitat patches. If the simplification of amphibians occupying metapopulations is accurate, then a regionally based conservation strategy, informed by metapopulation theory, is a powerful tool to estimate the isolation and extinction risk of ponds or populations. However, to date no attempt to assess the class‐wide generalization of amphibian populations as metapopulations has been made. We reviewed the literature on amphibians as metapopulations (53 journal articles or theses) and amphibian dispersal (166 journal articles or theses for 53 anuran species and 37 salamander species) to evaluate whether the conditions for metapopulation structure had been tested, whether pond isolation was based only on the assumption of limited dispersal, and whether amphibian dispersal was uniformly limited. We found that in the majority of cases (74%) the assumptions of the metapopulation paradigm were not tested. Breeding patch isolation via limited dispersal and/or strong site fidelity was the most frequently implicated or tested metapopulation condition, however we found strong evidence that amphibian dispersal is not as uniformly limited as is often thought. The frequency distribution of maximum movements for anurans and salamanders was well described by an inverse power law. This relationship predicts that distances beneath 11–13 and 8–9 km, respectively, are in a range that they may receive one emigrating individual. Populations isolated by distances approaching this range are perhaps more likely to exhibit metapopulation structure than less isolated populations. Those studies that covered larger areas also tended to report longer maximum movement distances – a pattern with implications for the design of mark‐recapture studies. Caution should be exercised in the application of the metapopulation approach to amphibian population conservation. Some amphibian populations are structured as metapopulations – but not all.     

Consequences of metapopulation collapse: comparison of genetic attributes between two Allegheny woodrat metapopulations

Disruptions in metapopulation connectivity due to demographic pressure can leave genetically isolated subpopulations susceptible to genetic drift, accumulation of deleterious alleles, and inbreeding depression. Such a scenario may be playing out within Allegheny woodrat (Neotoma magister) metapopulations as a series of synergistic extrinsic pressures have contributed to the rangewide decline of the species over the last 40 years. Our goal was to elucidate the effects of demographic collapse on metapopulation function by using 11 microsatellites markers to quantify differences in patterns of connectivity and genetic diversity between a demographically stable metapopulation and one in severe demographic decline. The demographically diminished metapopulation had lower levels of genetic diversity than the stable metapopulation at all levels evaluated (metapopulation-, subpopulation-, and individual-scales). In contrast to patterns of connectivity observed within the stable metapopulation, peripheral subpopulations in the diminished metapopulation had become completely isolated and were drifting toward genetic fixation, likely as a result of the extirpation of stepping-stone subpopulations. The declining genetic parameters observed within these isolated peripheral subpopulations suggest that inbreeding depression may be contributing significantly to their demographic decline. Allegheny woodrats readily express the genetic consequences of metapopulation decline due to the low effective population sizes of subpopulations and the species’ limited dispersal capacity. Differences in genetic parameters observed between demographically stable and diminished Allegheny woodrat metapopulations emphasize the risks posed to metapopulation function and associated genetic processes introduced with demographic decline.     

Dispersal and gene flow in the rare, parasitic Large Blue butterfly Maculinea arion

Dispersal is crucial for gene flow and often determines the long‐term stability of meta‐populations, particularly in rare species with specialized life cycles. Such species are often foci of conservation efforts because they suffer disproportionally from degradation and fragmentation of their habitat. However, detailed knowledge of effective gene flow through dispersal is often missing, so that conservation strategies have to be based on mark–recapture observations that are suspected to be poor predictors of long‐distance dispersal. These constraints have been especially severe in the study of butterfly populations, where microsatellite markers have been difficult to develop. We used eight microsatellite markers to analyse genetic population structure of the Large Blue butterfly Maculinea arion in Sweden. During recent decades, this species has become an icon of insect conservation after massive decline throughout Europe and extinction in Britain followed by reintroduction of a seed population from the Swedish island of Öland. We find that populations are highly structured genetically, but that gene flow occurs over distances 15 times longer than the maximum distance recorded from mark–recapture studies, which can only be explained by maximum dispersal distances at least twice as large as previously accepted. However, we also find evidence that gaps between sites with suitable habitat exceeding ～20 km induce genetic erosion that can be detected from bottleneck analyses. Although further work is needed, our results suggest that M. arion can maintain fully functional metapopulations when they consist of optimal habitat patches that are no further apart than ～10 km.     

Sex-biased dispersal in a rare butterfly and the implications for its conservation

The survival of many species may be dependent on their ability to exist in human-altered landscapes within metapopulations; in turn, metapopulation persistence is dictated by the ability of individuals to move effectively among patches to promote recolonization. The Taylor’s checkerspot butterfly (Euphydryas editha taylori) is a species that does not naturally occur in fragmented landscapes, yet it is now restricted to a handful of small isolated prairie habitats. Current recovery plans aim to establish a stable metapopulation; however, to date little is known about the species’ ability to move across the landscape. In 2010 and 2011, we conducted marking, tracking and boundary surveys to explore the movement dynamics of adults within two sites in Oregon, USA. Over the survey period, we marked 136 male butterflies, tracked 174 individuals and observed the behavior of 1,576 individual butterflies at site boundaries. Our study revealed a significant sex-bias in the movement dynamics of the Taylor’s checkerspot in both suitable habitat and surrounding matrix. Males were highly motile, whereas females appeared sedentary, rarely moving from their natal site. The limited dispersal behavior of females indicates that populations cannot persist naturally in a metapopulation and thus are at high risk of extinction. Based on our findings, we recommend that managers take proactive measures to increase or enable dispersal (including translocation) to existing and/or restored sites.     

Habitat preferences of Maculinea arion and its Myrmica host ants: implications for habitat management in Italian Alps

Projections of climate-change scenarios indicate that many Maculinea arion populations will disappear from Europe over the next 50 years. Extinctions will be particularly concentrated around the southern limits of the species?? range, such as the Italian peninsula, except mountain populations, mainly in the Alps. M. arion is a social parasite of Myrmica ants and its survival depends on the presence and abundance of two resources, i.e. (1) a specific larval foodplant and (2) a suitable Myrmica host ant. We collected data on Thymus plants distribution, Myrmica ants abundance, turf height around Myrmica nests, distance between nests and the closest thyme plant (Thymus pulegioides), and vegetation structure at 14 patches of a M. arion fragmented population occurring in the Western Italian Alps (Val Ferret: Aosta). We found that patches with the highest abundance of M. arion had significantly higher foodplant abundance. Among 407 nests belonging to nine Myrmica species, we found M. arion larvae in one Myrmica lonae and in two Myrmica sulcinodis nests. The distribution of M. sulcinodis nests showed the best spatial overlap with Thymus plants. M. sulcinodis was also significantly more abundant on pastures showing high M. arion density, and therefore offered higher adoption probabilities to develop butterfly larvae. On M. arion habitat patches, the mean turf height around Myrmica nests varied from 10.5 to 22.3 cm and best matched the habitat requirements of T. pulegioides. Thus, we suggest that extensive grazing is needed to maintain the two most important resources for M. arion. Moreover, patches where M. arion was abundant included more numerous plant species characteristic of these mountain pastures, as well as higher plant diversity (Shannon index). As a consequence, M. arion represents an umbrella species and a good indicator for the conservation status of mountain grasslands.     

Maintenance of the sexes and persistence of a clonal organism in spatially complex metapopulations

Clonal organisms persist at a range of population sex ratios, from equal numbers of males and females to single-sex systems. When intersexual competition is strong enough to drive one sex locally extinct, the maintenance of the sexes is facilitated by the semi-independent dynamics of populations within a metapopulation. These semi-independent dynamics are influenced by dispersal and recolonization rates, which are affected by the spatial arrangement of populations. To establish the quantitative relationship between spatially complex metapopulations and the maintenance of the sexes, we used a mathematical model of the liverwort Marchantia inflexa. This clonal organism is found in discrete patches on rocks and along the banks of streams, which form single-sex and two-sex metapopulations. In this system, asexual propagules mainly disperse short distances. Long-distance between-patch dispersal and recolonization mainly occurs via sexual propagules, which require both sexes to be present. Dispersal of these two types of propagules could interact with the spatial arrangement of populations to affect the maintenance of the sexes. With our mathematical model, we found that at intermediate distances between populations, metapopulations maintained both sexes, and the spatial arrangement of populations changed the threshold at which one sex was lost. On the other hand, when populations were close to one another, one sex was lost and the single-sex metapopulation persisted through dispersal of asexual propagules. When populations were far apart, one sex was lost, and the metapopulation either went extinct due to lack of recolonization by asexual propagules or persisted because clumped populations facilitated recolonization. These idealized spatial arrangements help clarify the effects of the spatial arrangement on the maintenance of the sexes and the persistence of metapopulations of clonal organisms, which can help explain geographic parthenogenesis and the distribution of asexual populations, the persistence of asexual species, and inform the conservation of clonal organisms.     

Evolution of body condition‐dependent dispersal in metapopulations

Body condition‐dependent dispersal strategies are common in nature. Although it is obvious that environmental constraints may induce a positive relationship between body condition and dispersal, it is not clear whether positive body conditional dispersal strategies may evolve as a strategy in metapopulations. We have developed an individual‐based simulation model to investigate how body condition–dispersal reaction norms evolve in metapopulations that are characterized by different levels of environmental stochasticity and dispersal mortality. In the model, body condition is related to fecundity and determined either by environmental conditions during juvenile development (adult dispersal) or by those experienced by the mother (natal dispersal). Evolutionarily stable reaction norms strongly depend on metapopulation conditions: positive body condition dependency of dispersal evolved in metapopulation conditions with low levels of dispersal mortality and high levels of environmental stochasticity. Negative body condition‐dependent dispersal evolved in metapopulations with high dispersal mortality and low environmental stochasticity. The latter strategy is responsible for higher dispersal rates under kin competition when dispersal decisions are based on body condition reached at the adult life stage. The evolution of both positive and negative body condition‐dependent dispersal strategies is consequently likely in metapopulations and depends on the prevalent environmental conditions.     

How Levins’ dynamics emerges from a Ricker metapopulation model

Understanding the dynamics of metapopulations close to extinction is of vital importance for management. Levins-like models, in which local patches are treated as either occupied or empty, have been used extensively to explore the extinction dynamics of metapopulations, but they ignore the important role of local population dynamics. In this paper, we consider a stochastic metapopulation model where local populations follow a stochastic, density-dependent dynamics (the Ricker model), and use this framework to investigate the behaviour of the metapopulation on the brink of extinction. We determine under which circumstances the metapopulation follows a time evolution consistent with Levins’ dynamics. We derive analytical expressions for the colonisation and extinction rates (c and e) in Levins-type models in terms of reproduction, survival and dispersal parameters of the local populations, providing an avenue to parameterising Levins-like models from the type of information on local demography that is available for a number of species. To facilitate applying our results, we provide a numerical algorithm for computing c and e.     

Climate-dependent dispersal rates in metapopulations of burnet moths

Predicted climate change implies warmer weather and a higher frequency of extreme weather situations. The consequences of the warm July in 2003 was investigated in contrast to the cold July in 2004 in southern Sweden with focus on the dispersal rate of two species of burnet moths (Zygaenidae). During an extensive mark-release-recapture program in metapopulations of burnet moths substantial differences in inter habitat patch dispersal rates were observed. For two species of Zygaena it was 2.14 and 2.18 times higher during the warm year. Measured patch dispersal rate significantly declined towards the edges of the large study-area, suggesting that individuals disperse outside the study-area. No significant differences in dispersal rates were found between the two species studied. The dispersal rates were similar for both females and males, and no density dependent effects on dispersal or effects of amount of nectar rich flowers were found. Dispersal events appeared as distinct movements, often reaching several kilometres. We conclude that dispersal depend on climatic conditions. Exceptionally warm years may be especially important for survival of certain species in metapopulation systems and this may be crucial for the survival of populations during scenarios of climate change.     

Abundance�Coccupancy relationships in metapopulations: examples of rock pool Daphnia

Intraspecific positive relationships between abundance and occupancy are observed for many species, suggesting that the same processes drive local and regional species dynamics. Two main groups of mechanisms explain this relationship: spatiotemporal variation in local population growth rates due to variation in habitat quality, or dispersal effects that increase occupancy of a species when locally abundant. Several studies show that spatiotemporal variation in population growth rates causes positive abundance?Coccupancy relationships, but few have shown dispersal effects. It is believed that such effects should be more evident for species whose dispersal is limited, e.g. metapopulations, but those studies are limited. This study investigates abundance?Coccupancy relationships in three Daphnia metapopulations in rock pools and the degree to which dispersal or habitat quality affect their local abundances and occurrence. Daphnia longispina and Daphnia magna showed positive abundance?Coccupancy relationships, but not Daphnia pulex. No single ecological factor could explain the abundance?Coccupancy relationships of any given species. Instead, dispersal processes and rock pool quality (mainly salinity and depth) seem to act together to shape the abundance?Coccupancy relationships. Such a conclusion is also supported by an immigration experiment in natural rock pools. This study suggests that although positive abundance?Coccupancy relationships may be commonly found for metapopulations, both dispersal processes and variation in habitat quality can be factors determining the abundance?Coccupancy relationship of metapopulations experiencing habitat heterogeneity.     

Can the assumption of a non‐random search improve our prediction of butterfly fluxes between resource patches?

Abstract. 1. Understanding dispersal patterns that enable small, spatially isolated populations to survive in fragmented landscapes has become an important issue in conservation biology and landscape management. However, for most of the species of interest it is not known whether dispersing individuals navigate or follow systematic search strategies, as opposed to moving randomly. 2. Recently it was shown that individuals of the butterfly species Maniola jurtina do not seek resources by means of random flight. If true, this may be problematic for existing metapopulation models, including those based on the evolution of dispersal rates in metapopulations. 3. The study tested to what extent the non‐random dispersal patterns described in the literature can explain M. jurtina fluxes in its natural habitat. 4. A model based on literature assumptions of M. jurtina movements is presented in the work reported here, and its predictions are compared with 2 years of capture–recapture data on its fluxes in two landscapes. 5. The model provides a good fit to the data and gives better predictions than the model based only on patch sizes and distances between patches. 6. Thus, if data are available about the actual landscape under consideration, the model should be preferred to simpler approaches; however, in general theoretical considerations the simple approach based on patch size and the degree of its isolation will retain its value.     

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.     

Founder events as determinants of within-island and among-island genetic structure of Daphnia metapopulations

The genetic structure of metapopulations offers insights into the genetic consequences of local extinction and recolonization. We studied allozyme variation in rock pool metapopulations of two species of waterfleas (Daphnia) with the aim to understand how these dynamics influence genetic differentiation. We screened 138 populations of D. magna and 65 populations of D. longispina from an area in the archipelago of southern Finland. The pools from which they were sampled are separated by distances between 1.5 and 4710 m and located on a total of 38 islands. The genetic population structure of the two species was strikingly similar, consistent with their similar metapopulation ecology. The mean F(PT) value (differentiation among pools with respect to the total metapopulation) was 0.55 and a hierarchical analysis showed that genetic differentiation was strong (>0.25) among pools within islands as well as among whole islands. Within islands, pairwise genetic differentiation increased with geographic distance, indicating isolation by distance due to spatially limited dispersal. Previous studies have shown strong founder events occurring during colonization in our metapopulation. We suggest that the genetic population structure in the studied metapopulations is largely explained by three consequences of these founder events: (i) strong drift during colonization, (ii) local inbreeding, which results in hybrid vigour and increased effective migration rates after subsequent immigration, and (iii) effects of selection through hitchhiking of neutral genes with linked loci under selection.     

Evidence for positive density-dependent emigration in butterfly metapopulations

A positive effect of (meta)population density on emigration has been predicted by many theoretical models and confirmed empirically in various organisms. However, in butterflies, the most popular species for dispersal studies, the evidence for its existence has so far been equivocal, with negative relationships between density and emigration being reported more frequently. We analysed dispersal in sympatric metapopulations of two Maculinea butterflies, intensively surveyed with mark–release–recapture methods for 7 years. Dispersal parameters, derived using the virtual migration model, were assessed against butterfly densities, which fluctuated strongly over the study period. Emigration was positively correlated with density, and this effect was particularly strong at densities above carrying capacity, when emigration increased up to threefold in females and twofold in males compared with the normal levels. In turn, density had little impact on other dispersal parameters analysed. Our findings provide good evidence for positive density-dependence of emigration in butterflies. Emigrating at high densities is particularly beneficial for females, because it gives them a chance to lay part of their egg-load in less crowded patches, where offspring survival is higher due to lower intraspecific competition. Even though the rise in emigration becomes considerable at densities exceeding carrying capacity, i.e. relatively infrequently, it still has serious implications for many ecological phenomena, such as species range expansions, gene flow, and metapopulation persistence. Consequently, instead of treating emigration as a fixed trait, it is worth allowing for its density-dependence in applications such as population viability analyses, genetic models or metapopulation models.     

The influence of stage-dependent dispersal on the population dynamics of three amphipod species

In metapopulations, the maintenance of local populations can depend on source–sink dynamics, where populations with positive growth rate seed populations with negative growth rate. The pattern and probability of successful dispersal among habitats can therefore be crucial in determining whether local populations will become rare or increase in abundance. We present here data on the dispersal strategy and population dynamics of three marine amphipods living in pen shells (Atrina rigida) in the Gulf of Mexico. The three amphipod species in this study disperse at different life stages. Neomegamphopus hiatus and Melita nitida disperse as adults, while Bemlos unicornis disperses as juveniles. The two species that disperse as adults have the highest initial population sizes when a new shell becomes available, likely caused by the arriving females releasing their brood into these recently occupied shells. This dispersal pattern results in initially higher population growth, but fewer occupied shells, as noted by their clumped distribution. In contrast, the species that disperses as juveniles accumulates more slowly and more evenly across habitats, eventually dominating the other two in terms of numerical abundance. The metapopulation dynamics of the three species seems to be highly dependent on the life history stage involved in dispersal. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Metapopulation dynamics in a changing climate: Increasing spatial synchrony in weather conditions drives metapopulation synchrony of a butterfly inhabiting a fragmented landscape

Habitat fragmentation and climate change are both prominent manifestations of global change, but there is little knowledge on the specific mechanisms of how climate change may modify the effects of habitat fragmentation, for example, by altering dynamics of spatially structured populations. The long‐term viability of metapopulations is dependent on independent dynamics of local populations, because it mitigates fluctuations in the size of the metapopulation as a whole. Metapopulation viability will be compromised if climate change increases spatial synchrony in weather conditions associated with population growth rates. We studied a recently reported increase in metapopulation synchrony of the Glanville fritillary butterfly (Melitaea cinxia) in the Finnish archipelago, to see if it could be explained by an increase in synchrony of weather conditions. For this, we used 23 years of butterfly survey data together with monthly weather records for the same period. We first examined the associations between population growth rates within different regions of the metapopulation and weather conditions during different life‐history stages of the butterfly. We then examined the association between the trends in the synchrony of the weather conditions and the synchrony of the butterfly metapopulation dynamics. We found that precipitation from spring to late summer are associated with the M. cinxia per capita growth rate, with early summer conditions being most important. We further found that the increase in metapopulation synchrony is paralleled by an increase in the synchrony of weather conditions. Alternative explanations for spatial synchrony, such as increased dispersal or trophic interactions with a specialist parasitoid, did not show paralleled trends and are not supported. The climate driven increase in M. cinxia metapopulation synchrony suggests that climate change can increase extinction risk of spatially structured populations living in fragmented landscapes by altering their dynamics.     

Metapopulation dynamics in an aphid-parasitoid system

Metapopulation theory makes a number of predictions concerning the effects of dispersal on the persistence of predator-prey or host-parasitoid systems. While the stabilising effects of dispersal have been shown in a number of laboratory studies, evidence from field studies remains scarce due to a lack of suitable model systems. I describe a host-parasitoid system that shows a classical metapopulation structure with frequent extinctions and colonisations consisting of the aphidiid Lysiphlebus hirticornisand the aphid Metopeurum fuscoviride. Both the parasitoid and the aphid are specialists on their respective hosts. I followed the dynamics of host and parasitoid on individually marked tansy (Tanacetum vulgare) plants, the host of M. fuscoviride. Dynamics of host and parasitoid populations were characterized by frequent extinctions and colonisations. Mean longevity of aphid colonies was only 3.1 weeks. Parasitism by L. hirticorniswas a main cause of extinction for the aphid as rates of parasitism often reached 100%, in particular towards the end of the field season. Patchiness in this system occurs at two spatial scales. Aphid colonies form on single tansy ramets = shoots but movements of aphid individuals among ramets within a particular tansy genet are frequent. Because aphids can persist on a genet for a large numer of generations, it is argued that local populations form on genets rather than ramets. The number of host and parasitoid extinctions described in this study exceeds the number of extinctions usually observed in field studies of host-parasitoid metapopulations. It is suggested that aphid-parasitoid systems such as the one studied in this paper may be good models to test the predictions of metapopulation theory.     

Population structure, habitat features and genetic structure of managed red deer populations

Management of game ungulates alters population structure and habitat features, with potential effects on genetic structure. Here, we study 26 red deer (Cervus elaphus) populations in Spain. We used census data and habitat features as well as genetic information at 11 microsatellite markers from 717 individuals. We found that metapopulations presented a distribution associated with forest interruptions. Within metapopulations, fences did not have a significant effect on red deer genetic structure. The metapopulations we studied presented similar population structure, but they differed in habitat features and genetic structure. The metapopulation with higher resource availability showed a genetic structure pattern in which genetic relatedness between geographically close individuals was high while relatedness between geographically distant individuals was low. Contrarily, the metapopulation with lower resource availability presented a genetic structure pattern in which the genetic relatedness between individuals of different populations was independent of the geographic distance. We discuss the possible connection between resource availability and genetic structure. Finally, we did not find any population or environmental variable related to genetic differentiation within metapopulations.     

Source-Sink Colonization as a Possible Strategy of Insects Living in Temporary Habitats

Continuous colonization and re-colonization is critical for survival of insect species living in temporary habitats. When insect populations in temporary habitats are depleted, some species may escape extinction by surviving in permanent, but less suitable habitats, in which long-term population survival can be maintained only by immigration from other populations. Such situation has been repeatedly described in nature, but conditions when and how this occurs and how important this phenomenon is for insect metapopulation survival are still poorly known, mainly because it is difficult to study experimentally. Therefore, we used a simulation model to investigate, how environmental stochasticity, growth rate and the incidence of dispersal affect the positive effect of permanent but poor (“sink”) habitats on the likelihood of metapopulation persistence in a network of high quality but temporary (“source”) habitats. This model revealed that permanent habitats substantially increase the probability of metapopulation persistence of insect species with poor dispersal ability if the availability of temporary habitats is spatio-temporally synchronized. Addition of permanent habitats to a system sometimes enabled metapopulation persistence even in cases in which the metapopulation would otherwise go extinct, especially for species with high growth rates. For insect species with low growth rates the probability of a metapopulation persistence strongly depended on the proportions of “source” to “source” and “sink” to “source” dispersal rates.