A stochastic metapopulation model accounting for habitat dynamics

A stochastic metapopulation model accounting for habitat dynamics is presented. This is the stochastic SIS logistic model with the novel aspect that it incorporates varying carrying capacity. We present results of Kurtz and Barbour, that provide deterministic and diffusion approximations for a wide class of stochastic models, in a form that most easily allows their direct application to population models. These results are used to show that a suitably scaled version of the metapopulation model converges, uniformly in probability over finite time intervals, to a deterministic model previously studied in the ecological literature. Additionally, they allow us to establish a bivariate normal approximation to the quasi-stationary distribution of the process. This allows us to consider the effects of habitat dynamics on metapopulation modelling through a comparison with the stochastic SIS logistic model and provides an effective means for modelling metapopulations inhabiting dynamic landscapes.     

A spatially structured metapopulation model with patch dynamics

Metapopulation models that incorporate both spatial and temporal structure are studied in this paper. The existence and stability of equilibria are provided, and an extinction threshold condition is derived which depends on patch dynamics (patch destruction and creation) and metapopulation dynamics (patch colonization and extinction). These results refine threshold conditions given by previous metapopulation models. By comparing landscapes with different spatial heterogeneities with respect to weighted long-term patch occupancies, we conclude that the pattern of a landscape is of overwhelming importance in determining metapopulation persistence and patch occupancy. We show that the same conclusion holds when a rescue effect is considered. We also derive a stochastic differential equations (SDE) model of the It? type based on our deterministic model. Our simulations reveal good agreement between the deterministic model and the SDE model.     

A fully coupled, mechanistic model for infectious disease dynamics in a metapopulation: movement and epidemic duration

The drive to understand the invasion, spread and fade out of infectious disease in structured populations has produced a variety of mathematical models for pathogen dynamics in metapopulations. Very rarely are these models fully coupled, by which we mean that the spread of an infection within a subpopulation affects the transmission between subpopulations and vice versa. It is also rare that these models are accessible to biologists, in the sense that all parameters have a clear biological meaning and the biological assumptions are explained. Here we present an accessible model that is fully coupled without being an individual-based model. We use the model to show that the duration of an epidemic has a highly non-linear relationship with the movement rate between subpopulations, with a peak in epidemic duration appearing at small movement rates and a global maximum at large movement rates. Intuitively, the first peak is due to asynchrony in the dynamics of infection between subpopulations; we confirm this intuition and also show the peak coincides with successful invasion of the infection into most subpopulations. The global maximum at relatively large movement rates occurs because then the infectious agent perceives the metapopulation as if it is a single well-mixed population wherein the effective population size is greater than the critical community size.     

A unifying framework for metapopulation dynamics

Many biologically important processes, such as genetic differentiation, the spread of disease, and population stability, are affected by the (natural or enforced) subdivision of populations into networks of smaller, partly isolated, subunits. Such "metapopulations" can have extremely complex dynamics. We present a new general model that uses only two functions to capture, at the metapopulation scale, the main behavior of metapopulations. We show how complex, structured metapopulation models can be translated into our generalized framework. The metapopulation dynamics arising from some important biological processes are illustrated: the rescue effect, the Allee effect, and what we term the "antirescue effect." The antirescue effect captures instances where high migration rates are deleterious to population persistence, a phenomenon that has been largely ignored in metapopulation conservation theory. Management regimes that ignore a significant antirescue effect will be inadequate and may actually increase extinction risk. Further, consequences of territoriality and conspecific attraction on metapopulation-level dynamics are investigated. The new, simplified framework can incorporate knowledge from epidemiology, genetics, and population biology in a phenomenological way. It opens up new possibilities to identify and analyze the factors that are important for the evolution and persistence of the many spatially subdivided species.     

Shallow lake sediments provide evidence for metapopulation dynamics: a pilot study

The aim of this pilot study is to test the hypotheses that sediment cores can provide evidence for metapopulation dynamics and that these can be linked with site characteristics. We focus on temporal patterns of incidence and abundance of overwintering stages (statoblasts) produced by the freshwater bryozoan Cristatella mucedo, an organism characterised by a metapopulation ecology, in sediment cores retrieved from 18 UK lakes. Runs and goodness-of-fit tests provided evidence for population instability, periods of low abundance and absence, and of asynchrony—all signatures of metapopulation dynamics. Further hypothesis testing indicated that extinction risk is greater in more isolated sites and in sites of smaller size. Absence of statoblasts from the top sections of spatially separated, replicate cores provided independent evidence for extinction in one site. Our study demonstrates how the abundances of sedimentary-bound propagules may be analysed initially for metapopulation dynamics and subsequently how this may lead to working hypotheses regarding the drivers of such dynamics. The sediment archive represents a unique historical record whose potential for characterising metapopulation dynamics has previously been untapped but is broadly relevant for understanding the population biology of freshwater organisms.     

Synchronism in a metapopulation model

We consider a spatially explicit metapopulation model with interaction among the two nearest neighbors to relate, with a simple mathematical expression, chaos in the local, uncoupled, populations, the degree of interaction among patches, size of the metapopulation, and the stability of the synchronized attractor. Since synchronism is strongly correlated with extinction, our results can provide useful information on factors leading to population extinction.     

Competition between unit-restricted fungi: a metapopulation model

We aimed to provide a theoretical framework for dynamic studies of competition between fungi living on divided and ephemeral resources. We previously adapted the seminal Skellam's patch-occupancy model (Skellam, 1951) to describe the population dynamics of one species of unit-restricted fungus whose mycelial growth occurs within resource units and which colonizes new resource units by spore dispersal (Gourbiere et al., 1999). In this study, we extended this model to describe the competition between a pair of unit-restricted fungal species that interact with each other inside units by decreasing their spore production. Accordingly, we designed a discrete-time metapopulation model where all patches go extinct at each generation and species interact by lowering their propagule production in jointly occupied patches. We showed that the two species easily coexist although there is no trade-off between their competitive and colonization abilities. Furthermore, the outcome of the competition process can depend on a founder effect. Founder effect determines either which species is excluded or the relative densities of each species when they coexist. We investigated the implications of these results on the distribution and abundance of fungal species along environmental gradients. This work bridges the gap between the mycological theory of "Resource Units" and the metapopulation theory, showing the specificity of fungal exploitation competition. We suggest that unit-restricted fungal species are appropriate biological models to test the theoretical results of the metapopulation theory, such as the appearance of alternative stable equilibria.     

Stochastic dynamics in exotic and native blowflies: an analysis combining laboratory experiments and a two-patch metapopulation model

In this study we explored the stochastic population dynamics of three exotic blowfly species, Chrysomya albiceps, Chrysomya megacephala and Chrysomya putoria, and two native species, Cochliomyia macellaria and Lucilia eximia, by combining a density-dependent growth model with a two-patch metapopulation model. Stochastic fecundity, survival and migration were investigated by permitting random variations between predetermined demographic boundary values based on experimental data. Lucilia eximia and Chrysomya albiceps were the species most susceptible to the risk of local extinction. Cochliomyia macellaria, C. megacephala and C. putoria exhibited lower risks of extinction when compared to the other species. The simultaneous analysis of stochastic fecundity and survival revealed an increase in the extinction risk for all species. When stochastic fecundity, survival and migration were simulated together, the coupled populations were synchronized in the five species. These results are discussed, emphasizing biological invasion and interspecific interaction dynamics.     

Host ecotype generates evolutionary and epidemiological divergence across a pathogen metapopulation

The extent and speed at which pathogens adapt to host resistance varies considerably. This presents a challenge for predicting when—and where—pathogen evolution may occur. While gene flow and spatially heterogeneous environments are recognized to be critical for the evolutionary potential of pathogen populations, we lack an understanding of how the two jointly shape coevolutionary trajectories between hosts and pathogens. The rust pathogen Melampsora lini infects two ecotypes of its host plant Linum marginale that occur in close proximity yet in distinct populations and habitats. In this study, we found that within-population epidemics were different between the two habitats. We then tested for pathogen local adaptation at host population and ecotype level in a reciprocal inoculation study. Even after controlling for the effect of spatial structure on infection outcome, we found strong evidence of pathogen adaptation at the host ecotype level. Moreover, sequence analysis of two pathogen infectivity loci revealed strong genetic differentiation by host ecotype but not by distance. Hence, environmental variation can be a key determinant of pathogen population genetic structure and coevolutionary dynamics and can generate strong asymmetry in infection risks through space.     

Aggregate statistical measures and metapopulation dynamics

There are two main types of metapopulation models. Spatially implicit models are analytically tractable but neglect spatial heterogeneities. Spatially explicit models are more realistic but too complex. In this paper, I build a bridge between both approximations. I derive a new metapopulation model using a well-known technique in population genetics. Spatial heterogeneities are captured by an aggregate statistical measure of spatial correlation. When this correlation is zero, i.e., space is homogeneous, the model becomes the well-known Levins' model. As spatial correlation increases, equilibrium patch occupancy decreases from what would be expected under the spatially homogeneous assumption. I proceed by testing how well spatial complexities from a spatially explicit simulation can be encapsulated by such an aggregate statistical measure.     

Optimal patch use and metapopulation dynamics

Summary We compared the metapopulation dynamics of predator—prey systems with (1) adaptive global dispersal, (2) adaptive local dispersal, (3) fixed global dispersal and (4) fixed local dispersal by predators. Adaptive dispersal was modelled using the marginal value theorem, such that predators departed patches when the instantaneous rate of prey capture was less than the long-term rate of prey capture averaged over all patches, scaled to the movement time between patches. Adaptive dispersal tended to stabilize metapopulation dynamics in a similar manner to conventional fixed dispersal models, but the temporal dynamics of adaptive dispersal models were more unpredictable than the smooth oscillations of fixed dispersal models. Moreover, fixed and adaptive dispersal models responded differently to spatial variation in patch productivity and the degree of compartmentalization of the system. For both adaptive dispersal and fixed dispersal models, localized (stepping-stone) dispersal was more strongly stabilizing than global (island) dispersal. Variation among predators in the probability of dispersal in relation to local prey density had a strong stabilizing influence on both within-patch and metapopulation dynamics. These results suggest that adaptive space use strategies by predators could have important implications for the dynamics of spatially heterogeneous trophic systems.     

Bubonic plague: a metapopulation model of a zoonosis.

Bubonic plague (Yersinia pestis) is generally thought of as a historical disease; however, it is still responsible for around 1000-3000 deaths each year worldwide. This paper expands the analysis of a model for bubonic plague that encompasses the disease dynamics in rat, flea and human populations. Some key variables of the deterministic model, including the force of infection to humans, are shown to be robust to changes in the basic parameters, although variation in the flea searching efficiency, and the movement rates of rats and fleas will be considered throughout the paper. The stochastic behaviour of the corresponding metapopulation model is discussed, with attention focused on the dynamics of rats and the force of infection at the local spatial scale. Short-lived local epidemics in rats govern the invasion of the disease and produce an irregular pattern of human cases similar to those observed. However, the endemic behaviour in a few rat subpopulations allows the disease to persist for many years. This spatial stochastic model is also used to identify the criteria for the spread to human populations in terms of the rat density. Finally, the full stochastic model is reduced to the form of a probabilistic cellular automaton, which allows the analysis of a large number of replicated epidemics in large populations. This simplified model enables us to analyse the spatial properties of rat epidemics and the effects of movement rates, and also to test whether the emergent metapopulation behaviour is a property of the local dynamics rather than the precise details of the model.     

An epidemiological model for the dynamics of Chagas' disease.

A model for the transmission dynamics of Chagas' disease is presented. The structure of the model is similar to that of the Ross-Macdonald model for malaria but includes an extra infectious compartment (chronically ill individuals) which is characteristic of Chagas' disease. In Chagas' disease there are two-main forms of transmission, by blood transfusion and by vector biting. The former is more common in urban environments and the latter is characteristic of rural settings. The characteristic long chronic (frequently asymptomatic) stage of Chagas' disease is potentially a risk factor that could enhance disease transmission by blood transfusion. The model evaluates the relative importance of both transmission modes in populations of constant size. The main results indicate that there is a strong tendency of the disease to reach an asymptotically stable endemic equilibrium point. Also, the magnitude of the basic reproductive number is very sensitive to the length of the chronic stage of the disease and hence it follows that early detection of cases (reducing the length of this stage) is important for the eventual eradication of the disease.     

From individual behavior to metapopulation dynamics: unifying the patchy population and classic metapopulation models

Spatially structured populations in patchy habitats show much variation in migration rate, from patchy populations in which individuals move repeatedly among habitat patches to classic metapopulations with infrequent migration among discrete populations. To establish a common framework for population dynamics in patchy habitats, we describe an individual-based model (IBM) involving a diffusion approximation of correlated random walk of individual movements. As an example, we apply the model to the Glanville fritillary butterfly (Melitaea cinxia) inhabiting a highly fragmented landscape. We derive stochastic patch occupancy model (SPOM) approximations for the IBMs assuming pure demographic stochasticity, uncorrelated environmental stochasticity, or completely correlated environmental stochasticity in local dynamics. Using realistic parameter values for the Glanville fritillary, we show that the SPOMs mimic the behavior of the IBMs well. The SPOMs derived from IBMs have parameters that relate directly to the life history and behavior of individuals, which is an advantage for model interpretation and parameter estimation. The modeling approach that we describe here provides a unified framework for patchy populations with much movements among habitat patches and classic metapopulations with infrequent movements.     

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.     

Population dynamics of epiphytic orchids in a metapopulation context

Background and Aims

Populations of many epiphytes show a patchy distribution where clusters of plants growing on individual trees are spatially separated and may thus function as metapopulations. Seed dispersal is necessary to (re)colonize unoccupied habitats, and to transfer seeds from high- to low-competition patches. Increasing dispersal distances, however, reduces local fecundity and the probability that seeds will find a safe site outside the original patch. Thus, there is a conflict between seed survival and colonization.

Methods

Populations of three epiphytic orchids were monitored over three years in a Mexican humid montane forest and analysed with spatially averaged and with spatially explicit matrix metapopulation models. In the latter, population dynamics at the scale of the subpopulations (epiphytes on individual host trees) are based on detailed stage-structured observations of transition probabilities and trees are connected by a dispersal function.

Key Results

Population growth rates differed among trees and years. While ignoring these differences, and averaging the population matrices over trees, yields negative population growth, metapopulation models predict stable or growing populations because the trees that support growing subpopulations determine the growth of the metapopulation. Stochastic models which account for the differences among years differed only marginally from deterministic models. Population growth rates were significantly lower, and extinctions of local patches more frequent in models where higher dispersal results in reduced local fecundity compared with hypothetical models where this is not the case. The difference between the two models increased with increasing mean dispersal distance. Though recolonization events increased with dispersal distance, this could not compensate the losses due to reduced local fecundity.

Conclusions

For epiphytes, metapopulation models are useful to capture processes beyond the level of the single host tree, but local processes are equally important to understand epiphyte population dynamics.     

Modelling epiphyte metapopulation dynamics in a dynamic forest landscape

We combine simulations with spatial statistics to estimate the parameters of a metapopulation model for the epiphytic lichen Lobaria pulmonaria specializing on aspen ( Populus tremula ) and goat willow ( Salix caprea ) in Fennoscandian boreal old-growth forests. We estimated the parameters of a forest landscape model (FIN-LANDIS) by repeatedly running simulations and selecting the set of parameters for tree ecology and fire regime that reproduced empirical host tree density and spatial patterns. Second, we tested which variables were important in epiphyte colonization and estimated the dispersal kernel. Third, we run a metapopulation model for the lichen across the estimated landscape scenarios and selected values for the remaining parameters that reproduced the empirical patterns of epiphyte occurrence. There was little variation in predicted dynamics, occupancy and spatial patterns between replicate metapopulation simulations. However, more data would be required for accurately estimating the parameters of FIN-LANDIS primarily because of the inherent stochasticity in large scale forest fires. Following the beginning of fire suppression in the study area 150 years ago, the model predicts that lichen occupancy first increases but subsequently declines. The lower occupancy in the past than at present is explained by high rate of tree destruction by fires, which increases local extinction rate in patch-tracking metapopulations. In the absence of fires, the occupancy increases because of lower extinction rate, but without forest fires or alternative means of host tree regeneration, the lichen is predicted to go ultimately extinct because of severely reduced density of aspen and goat willow.     

Extinction dynamics and the regional persistence of a tree frog metapopulation

The concept of a metapopulation acknowledges local extinctions as a natural part of the dynamics of a patchily distributed population. However, if extinctions are not balanced by recolonizations or if there is a high degree of spatial synchrony of local extinctions, this poses a threat to and will reduce the metapopulation persistence time. Here we show that, in a metapopulation network of 378 pond patches used by the tree frog (Hyla arborea), even though extinctions are frequent (mean extinction probability p(e) = 0.24) they pose no threat to the metapopulation as they are balanced by recolonizations (p(c) = 0.33). In any one year there was a pattern of large populations tending to persist while small populations became extinct. The total number of individuals belonging to populations that went extinct was small (< 5%) compared with those populations that persisted. A spatial autocorrelation analysis indicated no clustering of local extinctions. The tree frog metapopulation studied consisted of a set of larger, persistent populations mixed with smaller populations characterized by high turnover dynamics.