Host-parasitoid extinction and colonization in a fragmented prairie landscape

Few field studies of natural populations have examined the factors influencing local extinctions and colonization of empty habitat patches for a prey species and its predator. In this study, I carried out a census of planthopper (Prokelisia crocea; Hemiptera: Delphacidae) and egg parasitoid (Anagrus columbi; Hymenoptera: Mymaridae) incidence and densities in 147 host-plant patches (Spartina pectinata; Poaceae) over seven planthopper generations in a tall-grass prairie landscape. For both species, the likelihood of going extinct in a patch was related to a number of patch-specific variables: density, temporal variability in density, proportion of hosts parasitized (planthopper only), host-plant density, patch size, patch isolation, and composition of the surrounding matrix. Colonization likelihood was related only to the physical attributes of the patch. There was high patch turnover in this prairie landscape. On average, planthoppers went extinct in 23% of the patches and A. columbi went extinct in 51% of the patches in each generation. For the planthopper, extinction likelihood increased with a decrease in patch size and the proportion of the matrix composed of mudflat. Parasitism of eggs had no effect on the extinction likelihood of local P. crocea populations, suggesting that A. columbi may not play a major role in the patch dynamics of its host. The likelihood of extinction for A. columbi was dependent on factors that spanned three trophic levels. An increase in plant density, decrease in host density and decrease in parasitoid density all increased the likelihood of A. columbi extinction within a patch. The dependency on multiple trophic levels may explain the higher extinction risk for the parasitoid than its host. A. columbi extinction was also affected by the matrix habitat surrounding the patch—the effect was the opposite of that for P. crocea. Finally, vacant patches were colonized at rates of 53% and 34% per generation for the planthopper and parasitoid, respectively. For both species, colonization probabilities decreased with an increase in patch isolation. High host densities in a patch also favored high rates of colonization by A. columbi. I discuss how anthropogenic changes to the prairie landscape can affect the metapopulation dynamics and persistence time of this host-parasitoid interaction.     

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.     

Extinction-colonization dynamics and host-plant choice in butterfly metapopulations

Species living in highly fragmented landscapes often occur as metapopulations with frequent population turnover. Turnover rate is known to depend on ecological factors, such as population size and connectivity, but it may also be influenced by the phenotypic and genotypic composition of populations. The Glanville fritillary butterfly (Melitaea cinxia) in Finland uses two host-plant species that vary in their relative abundances among distinct habitat patches (dry meadows) in a large network of approximately 1,700 patches. We found no effect of host species use on local extinction. In contrast, population establishment was strongly influenced by the match between the host species composition of an empty habitat patch and the relative host use by larvae in previous years in the habitat patches that were well connected to the target patch. This "colonization effect" could be due to spatially variable plant acceptability or resistance or to spatially variable insect oviposition preference or larval performance. We show that spatial variation in adult oviposition preference occurs at the relevant spatial scale and that the other possible causes of the colonization effect can be discounted. We conclude that the colonization effect is generated by host preference influencing the movement patterns of ovipositing females. Migrant females with dissimilar host preferences have different perceptions of relative patch quality, which influences their likelihood of colonizing patches with particular host composition.     

Metapopulation dynamics of a beetle species confined to burned forest sites in a managed forest region

Despite increasing awareness of the theoretical importance of habitat dynamics on metapopulations, only a few empirical studies have been conducted. We aimed to increase our understanding of how patch size, dynamics and connectivity affect colonization–extinction dynamics and the occurrence patterns of a beetle (Stephanopachys linearis), which breeds only in burned trees, existing as dynamic habitat patches that have become rare in managed forest landscapes. We assessed species’ presence/absence twice in all known habitat patches (i.e. > 1 ha sites where forest fires had occurred during the previous 2–15 yr) in a 200 × 150 km region of central Sweden, dominated by managed boreal forest. Evaluated over six years, the colonization rate was 47% and the local extinction risk was 65%. Probability of colonization increased with patch size (number of suitable trees in a site) and connectivity to occupied patches within 30 km, and decreased with increasing time since fire. Local extinction risk decreased with habitat patch size but increased, unexpectedly, with connectivity. Occurrence increased with patch size and decreased with increasing time since fire. At a regional scale, S. linearis tracks the fire dynamics by colonising sites with burned trees and by becoming extinct at rates which make the species rare at sites where burnt trees are more than eight years old. In managed boreal forest landscapes, a large proportion of sites may be created by prescribed burning (in our study area: 82%), and consequently human decisions strongly affect the future amount of habitat for fire‐dependent species and its spatial distribution. Stephanopachys linearis uses burned sites more often if more trees are retained and, to some extent, if sites are concentrated in those parts of a region that already support high population densities of the species.     

Relative importance of host plant patch geometry and habitat quality on the patterns of occupancy,extinction and density of the monophagous butterfly <Emphasis Type="Italic">Iolana iolas</Emphasis>

Habitat fragmentation is a major cause of species rarity and decline because it increases local population extinctions and reduces recolonisation rates of remnant patches. Although two major patch characteristics (area and connectivity) have been used to predict distribution patterns in fragmented landscapes, other factors can affect the occurrence of a species as well as the probability of it becoming extinct. In this paper, we study the spatial structure and dynamics of the butterfly Iolana iolas in a 75-patch network of its host plant (Colutea hispanica) to determine the relative importance of patch area, connectivity and habitat quality characteristics on occupancy, extinction and density over the period 2003–2006. Occupancy in 2003, incidence (proportion of years occupied) and probability of extinction were mostly affected by patch area. Smaller patches were less likely to be occupied because they had a higher probability of extinction, partly due to environmental stochasticity. The density of I. iolas was negatively related to patch area in all study years. Only in 2004 was the density of I. iolas positively influenced by fruit production per plant. Our results suggest that for I. iolas, and probably for other specialist butterflies with clearly delimited resource requirements, metapopulation dynamics can be satisfactorily predicted using only geometric variables because most habitat characteristics are subsumed in patch area. However, this hypothesis should be subject to further testing under diverse environmental conditions to evaluate the extent of its generalisation.     

Testing metapopulation concepts: effects of patch characteristics and neighborhood occupancy on the dynamics of an endangered lagomorph

Metapopulation ecology is a field that is richer in theory than in empirical results. Many existing empirical studies use an incidence function approach based on spatial patterns and key assumptions about extinction and colonization rates. Here we recast these assumptions as hypotheses to be tested using 18 years of historic detection survey data combined with four years of data from a new monitoring program for the Lower Keys marsh rabbit. We developed a new model to estimate probabilities of local extinction and colonization in the presence of nondetection, while accounting for estimated occupancy levels of neighboring patches. We used model selection to identify important drivers of population turnover and estimate the effective neighborhood size for this system. Several key relationships related to patch size and isolation that are often assumed in metapopulation models were supported: patch size was negatively related to the probability of extinction and positively related to colonization, and estimated occupancy of neighboring patches was positively related to colonization and negatively related to extinction probabilities. This latter relationship suggested the existence of rescue effects. In our study system, we inferred that coastal patches experienced higher probabilities of extinction and colonization than interior patches. Interior patches exhibited higher occupancy probabilities and may serve as refugia, permitting colonization of coastal patches following disturbances such as hurricanes and storm surges. Our modeling approach should be useful for incorporating neighbor occupancy into future metapopulation analyses and in dealing with other historic occupancy surveys that may not include the recommended levels of sampling replication.     

Spatially structured metapopulation models: global and local assessment of metapopulation capacity

We model metapopulation dynamics in finite networks of discrete habitat patches with given areas and spatial locations. We define and analyze two simple and ecologically intuitive measures of the capacity of the habitat patch network to support a viable metapopulation. Metapopulation persistence capacity lambda(M) defines the threshold condition for long-term metapopulation persistence as lambda(M)>delta, where delta is defined by the extinction and colonization rate parameters of the focal species. Metapopulation invasion capacity lambda(I) sets the condition for successful invasion of an empty network from one small local population as lambda(I)>delta. The metapopulation capacities lambda(M) and lambda(I) are defined as the leading eigenvalue or a comparable quantity of an appropriate "landscape" matrix. Based on these definitions, we present a classification of a very general class of deterministic, continuous-time and discrete-time metapopulation models. Two specific models are analyzed in greater detail: a spatially realistic version of the continuous-time Levins model and the discrete-time incidence function model with propagule size-dependent colonization rate and a rescue effect. In both models we assume that the extinction rate increases with decreasing patch area and that the colonization rate increases with patch connectivity. In the spatially realistic Levins model, the two types of metapopulation capacities coincide, whereas the incidence function model possesses a strong Allee effect characterized by lambda(I)=0. For these two models, we show that the metapopulation capacities can be considered as simple sums of contributions from individual habitat patches, given by the elements of the leading eigenvector or comparable quantities. We may therefore assess the significance of particular habitat patches, including new patches that might be added to the network, for the metapopulation capacities of the network as a whole. We derive useful approximations for both the threshold conditions and the equilibrium states in the two models. The metapopulation capacities and the measures of the dynamic significance of particular patches can be calculated for real patch networks for applications in metapopulation ecology, landscape ecology, and conservation biology.     

Colonization and extinction dynamics of an annual plant metapopulation in an urban environment

The metapopulation framework considers that the spatiotemporal distribution of organisms results from a balance between the colonization and extinction of populations in a suitable and discrete habitat network. Recent spatially realistic metapopulation models have allowed patch dynamics to be investigated in natural populations but such models have rarely been applied to plants. Using a simple urban fragmented population system in which favourable habitat can be easily mapped, we studied patch dynamics in the annual plant Crepis sancta (Asteraceae). Using stochastic patch occupancy models (SPOMs) and multi‐year occupancy data we dissected extinction and colonization patterns in our system. Overall, our data were consistent with two distinct metapopulation scenarios. A metapopulation (sensu stricto) dynamic in which colonization occurs over a short distance and extinction is lowered by nearby occupied patches (rescue effect) was found in a set of patches close to the city centre, while a propagule rain model in which colonization occurs from a large external population was most consistent with data from other networks. Overall, the study highlights the importance of external seed sources in urban patch dynamics. Our analysis emphasizes the fact that plant distributions are governed not only by habitat properties but also by the intrinsic properties of colonization and dispersal of species. The metapopulation approach provides a valuable tool for understanding how colonization and extinction shape occupancy patterns in highly fragmented plant populations. Finally, this study points to the potential utility of more complex plant metapopulation models than traditionally used for analysing ecological and evolutionary processes in natural metapopulations.     

Habitat quality of source patches and connectivity in fragmented landscapes

Because spatial connectivity is critical to dispersal success and persistence of species in highly fragmented landscapes, the way that we envision and measure connectivity is consequential for biodiversity conservation. Connectivity metrics used for predictive modeling of spatial turnover and patch occupancy for metapopulations, such as with Incidence Function Models (IFM), incorporate distances to and sizes of possible source populations. Here, our focus is on whether habitat quality of source patches also is considered in these connectivity metrics. We propose that effective areas (weighted by habitat quality) of source patches should be better surrogates for population size and dispersal potential compared to unadjusted patch areas. Our review of a representative sample of the literature revealed that only 12.5% of studies incorporated habitat quality of source patches into IFM-type connectivity metrics. Quality of source patches generally was not taken into account in studies even if habitat quality of focal patches was included in analyses. We provide an empirical example for a metapopulation of a rare wetland species, the round-tailed muskrat (Neofiber alleni), demonstrating that a connectivity metric based on effective areas of source patches better predicts patch colonization and occupancy than a metric that used simple patch areas. The ongoing integration of landscape ecology and metapopulation dynamics could be hastened by incorporating habitat quality of source patches into spatial connectivity metrics applied to species conservation in fragmented landscapes.     

Is the population turnover of patchy‐distributed annuals determined by dormancy dynamics or dispersal processes?

In species with fragmented distribution, regional turnover dynamics is given by the processes of local population extinction and patch (re)colonization by migrants spreading from neighboring occupied patches. In plants with dormant stages (e.g. seeds) and limited dispersal capacity, regional dynamics based on dispersal processes can be overridden by pseudo-turnover determined by signals inducing or breaking dormancy (e.g. due to changes in habitat quality) resulting in a low importance of habitat configuration and size.
In this study, I investigated the turnover dynamics of 5 annual plant species growing on ant mounds of Lasius flavus over three years. I analyzed whether the grassland-scale dynamics of these annuals is influenced by dispersal processes, or alternatively, by pseudo-turnover of soil seed populations. For that purpose I 1) searched for populations formed from soil seeds only, 2) compared the relative contribution of the soil seed bank and seed rain for population restoration after disappearance from the vegetation and 3) investigated whether colonization and extinction events are affected by patch isolation. I assumed if population turnover was rather a result of the soil seed bank dynamics then spatial effects would be hard to detect.
In spite of the presence of populations formed from soil seed and the relatively more important soil seed bank for potential population reestablishment, turnover dynamics followed the predictions of metapopulation theory. Population appearance was more probable in larger and less isolated patches. Probability of disappearance increased with decrease of population size that was negatively influenced by the patch size and its isolation. These findings indicate dispersal processes to be important in the turnover dynamics and only limited contribution of soil seed populations. Their small effectiveness is probably related to the low chance of recurrent disturbance on the mound surface.     

Generalizing Levins metapopulation model in explicit space: models of intermediate complexity

A recent study [Harding and McNamara, 2002. A unifying framework for metapopulation dynamics. Am. Nat. 160, 173-185] presented a unifying framework for the classic Levins metapopulation model by incorporating several realistic biological processes, such as the Allee effect, the Rescue effect and the Anti-rescue effect, via appropriate modifications of the two basic functions of colonization and extinction rates. Here we embed these model extensions on a spatially explicit framework. We consider population dynamics on a regular grid, each site of which represents a patch that is either occupied or empty, and with spatial coupling by neighborhood dispersal. While broad qualitative similarities exist between the spatially explicit models and their spatially implicit (mean-field) counterparts, there are also important differences that result from the details of local processes. Because of localized dispersal, spatial correlation develops among the dynamics of neighboring populations that decays with distance between patches. The extent of this correlation at equilibrium differs among the metapopulation types, depending on which processes prevail in the colonization and extinction dynamics. These differences among dynamical processes become manifest in the spatial pattern and distribution of “clusters” of occupied patches. Moreover, metapopulation dynamics along a smooth gradient of habitat availability show significant differences in the spatial pattern at the range limit. The relevance of these results to the dynamics of disease spread in metapopulations is discussed.     

Single-species metapopulation dynamics: concepts, models and observations

This paper outlines a conceptual and theoretical framework for single-species metapopulation dynamics based on the Levins model and its variants. The significance of the following factors to metapopulation dynamics are explored: evolutionary changes in colonization ability; habitat patch size and isolation; compensatory effects between colonization and extinction rates; the effect of immigration on local dynamics (the rescue effect); and heterogeneity among habitat patches. The rescue effect may lead to alternative stable equilibria in metapopulation dynamics. Heterogeneity among habitat patches may give rise to a bimodal equilibrium distribution of the fraction of patches occupied in an assemblage of species (the core-satellite distribution). A new model of incidence functions is described, which allows one to estimate species' colonization and extinction rates on islands colonized from mainland. Four distinct kinds of stochasticity affecting metapopulation dynamics are discussed with examples. The concluding section describes four possible scenarios of metapopulation extinction.     

Consequences of insular population structure: Distribution and extinction of spruce grouse populations

Summary The population structure of the spruce grouse (Canachites canadensis) was studied in the Adirondack Mountains of New York, U.S.A. Twenty-five isolated habitat patches exist and are occupied by spruce grouse, with 7 suitable but unoccupied patches existing at the periphery of the range. The regional distribution and abundance of spruce grouse is correlated with the amount of lowland coniferous forest habitat. Unoccupied patches were significantly smaller and significantly farther from occupied patches than were other occupied patches. For all patches, as distance from the nearest occupied patch increased, the percent of occupied patches decreased linearly. I incorporated birth and death rates for spruce grouse into the MacArthur-Wilson survivorship model which closely predicted the proportion of occupied patches for an average population density (2.8 spruce grouse/100ha). For the same demographic parameters, extinction times were calculated which indicate that the 15 habitat patches of a carrying capacity of 3 female spruce grouse (100 ha) would have an average extinction time of less than 6 years. This in part accounts for the high proportion of these patches which are unoccupied. Extinctions and recolonizations of patches were observed during the study. The patterns of patch occupancy can partially be predicted based on their size, spatial arrangement, and the demographic characteristics of the spruce grouse.     

Host-plant patch qualities and presence of a likely competitor species affect the distribution and abundance of a rare British moth, <Emphasis Type="Italic">Cucullia lychnitis</Emphasis>

Understanding the environmental factors that affect the distribution and abundance of rare and threatened species can help identify priority sites for conservation action such as habitat management. The moth Cucullia lychnitis is identified in UK legislation as a priority species for conservation. Its conspicuous larvae feed on a patchily distributed host plant Verbascum nigrum, but are entirely absent from some host plant patches within the species’ range. Environmental variability among patches was investigated to understand factors affecting (i) patch occupancy by C. lychnitis and (ii) abundance of C. lychnitis, where present. Occupancy of individual V. nigrum plants within occupied patches was also investigated. The likelihood of patch occupancy increased with patch size, and decreased with patch isolation. Abundance of C. lychnitis was negatively correlated with patch isolation and with abundance of the weevil Cionus nigritarsis, which exploits the same host plant and may be a competitor. Within occupied patches, larvae of C. lychnitis were significantly more likely to be found on taller plants with more flower spikes; such plants are typical of established rather than newly-created patches. These results may help to guide efforts to conserve C. lychnitis, highlighting the importance of habitat connectivity and indicating potential management actions to promote features positively associated with larval incidence and abundance.     

The importance of forest patch networks for the conservation of the Thorn-tailed Rayaditos in central Chile

Conservation of forest birds in fragmented landscapes requires not only determining the critical patch characteristics influencing local population persistence but also identifying patch networks providing connectivity and suitable habitat conditions necessary to ensure regional persistence. In this study, we assessed the importance of patch attributes, patch connectivity, and network components (i.e., groups of interconnected patches) in explaining the occupancy pattern of the Thorn-tailed Rayadito (Aphrastura spinicauda), a forest bird species of central Chile. Using a daily movement threshold distance, we identified a total of 16 network components of sclerophyllous forest within the study area. Among those components, patch area and vegetation structure-composition were important predictors of patch occupancy. However, the inclusion of patch connectivity and component size (i.e., the area of a network component) into the models greatly increases the models’ accuracy and parsimony. Using the best-fitted model, a total of 33 patches were predicted to be occupied by rayaditos within the study area, but such occupied patches were distributed in only six network components. These results suggest that persistence of rayaditos in central Chile requires the maintenance of large single patches and patch networks providing habitat and connectivity.     

Population structure of a monophagous moth in a patchy landscape

1. The population structure of a monophagous noctuid moth, Abrostola asclepiadis , living on a patchily distributed perennial herb, Vincetoxicum hirundinaria is described. The study took place over 5 years at a landscape scale (about 12 km²).
2. Patch occupancy rates and population densities were studied in relation to patch size, degree of patch isolation, level of sun exposure and distance from the coast. In addition, flight tests in the laboratory were performed to estimate the potential dispersal capacity of the moth.
3. Occupancy rates were high and the likelihood of extinction depended on patch size. Small patches were less likely to be occupied than were large patches (> 10 m²). Sun-exposed patches were occupied for a lower proportion of years than were shaded patches. No distance effects could be discerned at the spatial scale of study, presumably because the insect is a strong flier.
4. Population densities in occupied patches decreased with increasing patch size. Furthermore, insect densities tended to increase with distance from the coast. Density changes in patches were synchronized.
5. The studied insect population can be described as a 'patchy population' sensu Harrison (1991) with spatially correlated population dynamics. These dynamics are superimposed on a landscape gradient.     

Conspecific and heterospecific attraction in assessments of functional connectivity

Functional connectivity is known to have an important, positive influence on species persistence. Measurements of functional connectivity traditionally focus on structural attributes of landscapes such as the distance and matrix type between habitat patches as well as on how species interact with those structural attributes. However, we propose that the social behavior of a species, through conspecific and heterospecific attraction, will also impact connectivity by changing how dispersers move with respect to each other and occupied patches. We analyzed functional connectivity patterns using circuit and graph theory for golden-headed lion tamarins (Leontopithecus chrysomelas) in Brazil under three scenarios. In the first scenario, we looked at connectivity without the effects of attraction under varying maximum dispersal distance and ecological movement cost thresholds. In the second scenario, we allowed dispersers to travel over more hostile matrix than they normally would to reach an occupied patch. In the final scenario, we allowed dispersers to move only to occupied patches. We found that, according to the first scenario, range-wide functional landscape connectivity for golden-headed lion tamarins is low at realistic maximum dispersal distance and movement cost thresholds. Incorporating the effects of conspecific or heterospecific attraction would increase functional connectivity, in the case of scenario two, or decrease functional connectivity, in the case of scenario three. Because conspecific/heterospecific attraction can have an impact on movement for some species, this factor should be incorporated in assessments of functional connectivity patterns for social species and others where patch occupancy is likely to influence the movements of dispersers.     

Spatial heterogeneity, source-sink dynamics, and the local coexistence of competing species

Patch occupancy theory predicts that a trade-off between competition and dispersal should lead to regional coexistence of competing species. Empirical investigations, however, find local coexistence of superior and inferior competitors, an outcome that cannot be explained within the patch occupancy framework because of the decoupling of local and spatial dynamics. We develop two-patch metapopulation models that explicitly consider the interaction between competition and dispersal. We show that a dispersal-competition trade-off can lead to local coexistence provided the inferior competitor is superior at colonizing empty patches as well as immigrating among occupied patches. Immigration from patches that the superior competitor cannot colonize rescues the inferior competitor from extinction in patches that both species colonize. Too much immigration, however, can be detrimental to coexistence. When competitive asymmetry between species is high, local coexistence is possible only if the dispersal rate of the inferior competitor occurs below a critical threshold. If competing species have comparable colonization abilities and the environment is otherwise spatially homogeneous, a superior ability to immigrate among occupied patches cannot prevent exclusion of the inferior competitor. If, however, biotic or abiotic factors create spatial heterogeneity in competitive rankings across the landscape, local coexistence can occur even in the absence of a dispersal-competition trade-off. In fact, coexistence requires that the dispersal rate of the overall inferior competitor not exceed a critical threshold. Explicit consideration of how dispersal modifies local competitive interactions shifts the focus from the patch occupancy approach with its emphasis on extinction-colonization dynamics to the realm of source-sink dynamics. The key to coexistence in this framework is spatial variance in fitness. Unlike in the patch occupancy framework, high rates of dispersal can undermine coexistence, and hence diversity, by reducing spatial variance in fitness.     

Implications of empirical data quality to metapopulation model parameter estimation and application

Parameter estimation is a critical step in the use of any metapopulation model for predictive purposes. Typically metapopulation studies assume that empirical data are of good quality and any errors are so insignificant that they can be ignored. However, three types of errors occur commonly in metapopulation data sets. First, patch areas can be mis-estimated. Second, unknown habitat patches may be located within or around the study area. Third, there may be false zeros in the data set, that is, some patches were observed to be empty while there truly was a population in the patch. This study investigates biases induced into metapopulation model parameter estimates by these three types of errors. It was found that mis-estimated areas influence the scaling of extinction risk with patch area; extinction probabilities for large patches become overestimated. Missing patches cause overestimation of migration distances and colonization ability of the species. False zeros can affect very strongly all model components, the extinction risk in large patches, intrinsic extinction rates in general, migration distances and colonization ability may become all overestimated. Biases in parameter estimates translate into corresponding biases in model predictions, which are serious particularly if metapopulation persistence becomes overestimated. This happens for example when the migration capability of the species is overestimated. Awareness of these biases helps in understanding seemingly anomalous parameter estimation results. There are also implications for field work: it may be reasonable to allocate effort so that serious errors in data are minimized. It is particularly important to avoid observing false zeros for large and/or isolated patches.     

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.