Theory for designing nature reserves for single species

We examine the question of the optimal number of reserves that should be established to maximize the persistence of a species. We assume that the mean time to extinction of a single population increases as a power of the habitat area, that there is a certain amount of habitat to be reserved, and that the aim is to determine how this habitat is most efficiently divided. The optimal configuration depends on whether the management objective is to maximize the mean time to extinction or minimize the risk of extinction. When maximizing the mean time to extinction, the optimal number of independent reserves does not depend on the amount of available habitat for the reserve system. In contrast, the risk of extinction is minimized when individual reserves are equal to the optimal patch size, making the optimal number of reserves linearly proportional to the amount of available habitat. A model that includes dispersal and correlation in the incidence of extinction demonstrates the importance of considering the relative rate at which these two factors decrease with distance between reserves. A small number of reserves is optimal when the mean time to extinction increases rapidly with habitat area or when risks of extinction are high.     

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.     

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.     

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.     

Matrix Matters:?Differences of Grand Skink Metapopulation Parameters in Native Tussock Grasslands and Exotic Pasture Grasslands

Modelling metapopulation dynamics is a potentially very powerful tool for conservation biologists. In recent years, scientists have broadened the range of variables incorporated into metapopulation modelling from using almost exclusively habitat patch size and isolation, to the inclusion of attributes of the matrix and habitat patch quality. We investigated the influence of habitat patch and matrix characteristics on the metapopulation parameters of a highly endangered lizard species, the New Zealand endemic grand skink (Oligosoma grande) taking into account incomplete detectability. The predictive ability of the developed zxmetapopulation model was assessed through cross-validation of the data and with an independent data-set. Grand skinks occur on scattered rock-outcrops surrounded by indigenous tussock (bunch) and pasture grasslands therefore implying a metapopulation structure. We found that the type of matrix surrounding the habitat patch was equally as important as the size of habitat patch for estimating occupancy, colonisation and extinction probabilities. Additionally, the type of matrix was more important than the physical distance between habitat patches for colonisation probabilities. Detection probability differed between habitat patches in the two matrix types and between habitat patches with different attributes such as habitat patch composition and abundance of vegetation on the outcrop. The developed metapopulation models can now be used for management decisions on area protection, monitoring, and the selection of translocation sites for the grand skink. Our study showed that it is important to incorporate not only habitat patch size and distance between habitat patches, but also those matrix type and habitat patch attributes which are vital in the ecology of the target species.     

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.     

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.     

Metapopulation persistence in fragmented landscapes: significant interactions between genetic and demographic processes

We formulated a mathematical model in order to study the joint influence of demographic and genetic processes on metapopulation viability. Moreover, we explored the influence of habitat structure, matrix quality and disturbance on the interplay of these processes. We showed that the conditions that allow metapopulation persistence under the synergistic action of genetic and demographic processes depart significantly from predictions based on a mere superposition of the effects of each process separately. Moreover, an optimal dispersal rate exists that maximizes the range of survival rates of dispersers under which metapopulation persists and at the same time allows the largest sustainable patch removal and patch‐size reduction. The relative impact of patch removal and patch‐size reduction depends both on matrix quality and the dispersal strategy of the species: metapopulation persistence is more affected by patch‐size reduction (patch removal) for low (high)‐dispersing species, in presence of a low (high) quality matrix. Avoidance of inbreeding, through increased dispersal when the rate of inbreeding in a population is large, has positive effects on low‐dispersing species, but impairs the persistence of high‐dispersing species. Finally, size heterogeneity between patches largely influences metapopulation dynamics; the presence of large patches, even at the expense of other patches being smaller, can have positive effects on persistence in particular for species of low dispersing ability.     

On optimal size and number of reserves for metapopulation persistence

Habitat fragmentation is generally considered to be detrimental to the persistence of natural populations. In nature management, one therefore tends to prefer few large nature reserves over many small nature reserves having equal total area. This paper examines whether this preference is warranted in a metapopulation framework with circular reserves (patches) by formulating the dependence of metapopulation persistence on the size and number of reserves, both of which depend on reserve radius if the total area is kept constant. Two measures of metapopulation persistence are used: R(0), the number of patches colonized by an occupied patch during its lifetime as an occupied patch, and T(e), the expected time to extinction. These two measures are functions of the extinction and colonization rates of the metapopulation. Several mechanisms for the extinction and colonization processes are formulated from which the dependence of these rates on reserve radius is calculated. It turns out that T(e)generally increases with reserve radius for all mechanisms, which supports the preference of few large reserves. However, R(0)supports this preference only in the case of some special, rather unrealistic, mechanisms. In many other, more realistic, cases an intermediate reserve size exists for which metapopulation persistence measured by R(0)is optimal.     

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.     

Metapopulation theory for fragmented landscapes

We review recent developments in spatially realistic metapopulation theory, which leads to quantitative models of the dynamics of species inhabiting highly fragmented landscapes. Our emphasis is in stochastic patch occupancy models, which describe the presence or absence of the focal species in habitat patches. We discuss a number of ecologically important quantities that can be derived from the full stochastic models and their deterministic approximations, with a particular aim of characterizing the respective roles of the structure of the landscape and the properties of the species. These quantities include the threshold condition for persistence, the contributions that individual habitat patches make to metapopulation dynamics and persistence, the time to metapopulation extinction, and the effective size of a metapopulation living in a heterogeneous patch network.     

The influence of patch demographics on metapopulations, with particular reference to successional landscapes

The classical view of metapopulations relates the regional abundance of a species to the balance between the extinction and colonization dynamics of identical local populations. Species in successional landscapes may represent the most appropriate examples of classical metapopulations. However, Levins‐type metapopulation models do not explicitly separate population loss due to successional habitat change from other causes of extinction. A further complication is that the chance of population loss due to successional habitat change may be related to the age of a patch. I developed simple patch occupancy models to include succession and included consideration of patch age structure to address two related questions: what are the implications of changes in patch demographic rates and when is a move to a structured patch occupancy model justified? Age‐related variation in patch demography could increase or decrease the equilibrium fraction of the available habitat occupied by a species when compared to the predictions of an unstructured model. Metapopulation persistence was enhanced when the age class of patches with the highest species occupancy suffered relatively low losses to habitat succession. Conversely, when the age class of patches with the highest species occupancy also had relatively high successional loss rates, extinction thresholds were higher that would be predicted by a simple unstructured model. Hence age‐related variation in patch successional rate introduces biases into the predictions of simple unstructured models. Such biases can be detected from field surveys of the fraction of occupied and unoccupied patches in each age class. Where a bias is demonstrated, unstructured models will not be adequate for making predictions about the effects of changing parameters on metapopulation size. Thinking in successional terms emphasizes how landscapes might be managed to enhance or reduce the patch occupancy by any particular metapopulation     

Information thresholds,habitat loss and population persistence in breeding birds

The combination of spatial structure and non‐linear population dynamics can promote the persistence of coupled populations, even when the average population growth rate of the patches seen in isolation would predict otherwise. This phenomenon has generally been conceptualized and investigated through the movement of individuals among patches that each holds many individuals, as in metapopulation models. However, population persistence can likewise increase as the result of individuals moving among sites (e.g. breeding territories) within in a single patch. Here I examine the latter: individuals making small‐scale informed decisions with respect to where to breed can promote population persistence in poor environments. Based on a simple algebraic model, I demonstrate information thresholds, and predict that greater information use is required for population persistence under lower spatial heterogeneity in habitat quality, all else equal. Second, I implement an individual‐based model to explore prior experience and prospecting on conspecific success within a more complex, and spatially heterogeneous environment. Uniquely, I jointly examine the effects of simulated habitat loss, spatial heterogeneity prior to habitat, and variation in information gathering on population persistence. I find that habitat loss accelerates population quasi‐extinction risk; however, information use reduces extinction probabilities in proportion to the level of information gathering. Per capita reproductive success declines with number of breeding sites, suggesting that information‐mediated Allee effects may contribute to extinction risk. In conclusion, my study suggests that populations in a changing world may be increasingly vulnerable to extinction where patch size and spatial heterogeneity constrain the effectiveness of information‐use strategies.     

Effects of patch characteristics and isolation on relative abundance of the scarce heath butterfly <Emphasis Type="Italic">Coenonympha hero</Emphasis> (Nymphalidae)

The scarce heath (Coenonympha hero) is an internationally threatened butterfly in Western Europe, where it occurs primarily on hay fields and abandoned arable land in a small-scale agricultural landscape of south-central Scandinavia. Due to afforestation, this habitat is becoming increasingly fragmented in Sweden, and it can be expected that the scarce heath will decline abruptly when threshold conditions for metapopulation persistence are no longer met. We used stepwise polychotomous logistic regression to compare habitat characteristics and isolation measures for patches that harbour large, small or no populations, respectively, in an area of south-western Sweden. We found that patch area, distance to the nearest large population and amount of Galium spp. explained a significant part of the variation in relative abundance among patches. Distance to nearest large population resulted in a better model to predict occupancy than both distance to the nearest inhabited patch and connectivity, which suggests that primarily large populations act as sources for small satellite populations. Today, sites of three of the eight larger populations in the study area have been planted with spruce or pine and will disappear within 20 years. We argue that the disappearance of these patches may very well lead to rapid extinction of the whole metapopulation system.     

Dispersal-induced desynchronization: from metapopulations to metacommunities

The conceptualization of fragmented populations in terms of metapopulation theory has become standard over the last three decades. It is well known that increases in between‐patch migration rates cause more synchronous population fluctuations and that this coherence increases the risk of global metapopulation extinction. Because species’ migration rates and the probability of individuals surviving migration events depend on the effective distance between patches, the benefit of improving conservation corridors or the matrix between habitat patches has been questioned. As populations occur in the context of larger communities, moving from a metapopulation to a metacommunity model framework is a natural extension to address the generality of these conclusions. We show how considering a metacommunity can modify the conclusion that decreasing the effective distance between habitat patches (via improving matrix quality or other measures) necessarily increases the degree of metapopulation synchrony. We show that decreases in effective between‐patch distance may deter population synchrony because of the simultaneous effect this change has on the migration patterns of other species. These results indicate that species interactions need to be considered when the effect of conservation measures on population synchrony, and ultimately persistence, is addressed.     

Long-term dynamics in a metapopulation of the American pika

A 20-yr study of a metapopulation of the American pika revealed a regional decline in occupancy in one part of a large network of habitat patches. We analyze the possible causes of this decline using a spatially realistic metapopulation model, the incidence function model. The pika metapopulation is the best-known mammalian example of a classical metapopulation with significant population turnover, and it satisfies closely the assumptions of the incidence function model, which was parameterized with data on patch occupancy. The model-predicted incidences of patch occupancy are consistent with observed incidences, and the model predicts well the observed turnover rate between four metapopulation censuses. According to model predictions, the part of the metapopulation where the decline has been observed is relatively unstable and prone to large oscillations in patch occupancy, whereas the other part of the metapopulation is predicted to be persistent. These results demonstrate how extinction-colonization dynamics may produce spatially correlated patterns of patch occupancy without any spatially correlated processes in local dynamics or extinction rate. The unstable part of the metapopulation gives an empirical example of multiple quasi equilibria in metapopulation dynamics. Phenomena similar to those observed here may cause fluctuations in species' range limits.     

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.     

Metapopulation dynamics: Does it help to have more of the same?

With the interest in conservation biology shifting towards processes from patterns, and to populations from communities, the theory of metapopulation dynamics is replacing the equilibrium theory of island biogeography as the population ecology paradigm in conservation biology. The simplest models of metapopulation dynamics make predictions about the effects of habitat fragmentation - size and isolation of habitat patches - on metapopulation persistence. The simple models may be enriched by considerations of the effects of demographic and environmental stochasticity on the size and extinction probability of local populations. Environmental stochasticity affects populations at two levels: it makes local extinctions more probable, and it also decreases metapopulation persistence time by increasing the correlation of extinction events across populations. Some controversy has arisen over the significance of correlated extinctions, and how they may affect the optimal subdivision of metapopulations to maximize their persistence time.     

Habitat destruction,habitat restoration and eigenvector-eigenvalue relations

According to metapopulation theory, the capacity of a habitat patch network to support the persistence of a species is measured by the metapopulation capacity of the patch network. Mathematically, metapopulation capacity is given by the leading eigenvalue lambda(M) of an appropriately constructed non-negative n x n matrix M, where n is the number of habitat patches. Both habitat destruction (in the sense of destruction of entire patches) and habitat deterioration (in the sense of partial destruction of patches) lower the metapopulation capacity of the patch network. The effect of gradual habitat deterioration is given by the derivative of lambda(M) with respect to patch attributes and may be straightforwardly evaluated by sensitivity analysis. In contrast, destruction of entire patches leads to a rank modification of matrix M, the effect of which on lambda(M) may be derived from eigenvector-eigenvalue relations. Eigenvector-eigenvalue relations have previously been analyzed only for symmetric matrices, which restricts their use in biological applications. In this paper I generalize some of the previous results by deriving eigenvector-eigenvalue relations for general non-symmetric matrices. In addition to the exact eigenvector-eigenvalue relations, I also derive eigenvalue perturbation formulae for rank-one modifications. These results lead to simple and intuitive approximation formulae, which may be used e.g. to assess the contribution of particular habitat patches to the metapopulation capacity of the landscape. The mathematical results presented are not restricted to the metapopulation context, but they should find a number of useful applications in biology, engineering and other applied sciences, where the removal (or addition) of matrix rows and columns often corresponds in a natural manner to decreasing (or increasing) the degrees of freedom of the focal system.     

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.