异质种群动态模型:破碎化景观动态模拟的新途径

景观破碎化导致物种以异质种群方式存活,使得基于异质种群动态模拟破碎化景观动态成为可能。异质种群动态模型的发展为景观动态模拟奠定了良好基础。根据空间处理方式的不同,异质种群模型可分为三大类,可不同程度地用于描述破碎化景观动态。(1)空间不确定异质种群模型,假定所有局域种群间均等互联,模型中不包含空间信息,仅能用于景观斑块动态描述;(2)空间确定异质种群模型,假设局域种群在二维空间上以规则格子形式排列,是一种准现实的空间处理方式,可用于景观动态的简单描述;(3)空间现实异质种群模型,包含了破碎化景观中局域种群的几何特征,可直接用于真实景观动态的模拟研究。空间现实的和基于个体的异质种群模型不但是未来异质种群模型发展的主流,也将成为未来破碎化景观动态研究的重要工具。为了更加准确完整地描述破碎化景观动态,不但应该综合运用已有的各种异质种群模型方法,更要引进新模型来刎画多物种、多变量、高维度、复杂连接的破碎化景观格局与过程。     

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.     

Metapopulation dynamics with quasi-local competition

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

Metapopulation dynamics for spatially extended predator–prey interactions

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

Migration, coherence and persistence in a fragmented landscape

The chance of local extinction is high during periods of small population size. Accordingly, a metapopulation made of local communities that support internal population cycling may face the threat of regional extinction if the local dynamics is coherent (synchronized). These systems achieve maximum sustainability at an intermediate level of migration that allows recolonization but prevents synchronization. Here we implement an individual-based simulation technique to examine the maximum persistence condition for a system of patch habitats connected by passive migration. The models discussed in this paper take into consideration realistic elements of metapopulations, such as migration cost, disordered spatial structure, frustration and environmental noise. It turns out that the state with maximum anti-correlation between neighboring patches is the most sustainable one, even in the presence of these complications. The results suggest, at least for small systems, a model independent conservation strategy: coherence between neighboring local communities has, in general, a negative impact, and population will benefit from intervention that increases anti-correlations.     

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 models for seasonally migratory animals

Metapopulation models are widely used to study species that occupy patchily distributed habitat, but are rarely applied to migratory species, because of the difficulty of identifying demographically independent subpopulations. Here, we extend metapopulation theory to describe the directed seasonal movement of migratory populations between two sets of habitat patches, breeding and non-breeding, with potentially different colonization and extinction rates between patch types. By extending the classic metapopulation model, we show that migratory metapopulations will persist if the product of the two colonization rates exceeds the product of extinction rates. Further, we develop a spatially realistic migratory metapopulation model and derive a landscape metric-the migratory metapopulation capacity-that determines persistence. This new extension to metapopulation theory introduces an important tool for the management and conservation of migratory species and may also be applicable to model the dynamics of two host-parasite systems.     

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.     

Stochastic models for mainland-island metapopulations in static and dynamic landscapes

This paper has three primary aims: to establish an effective means for modelling mainland-island metapopulations inhabiting a dynamic landscape; to investigate the effect of immigration and dynamic changes in habitat on metapopulation patch occupancy dynamics; and to illustrate the implications of our results for decision-making and population management. We first extend the mainland-island metapopulation model of Alonso and McKane [Bull. Math. Biol. 64:913–958, 2002] to incorporate a dynamic landscape. It is shown, for both the static and the dynamic landscape models, that a suitably scaled version of the process converges to a unique deterministic model as the size of the system becomes large. We also establish that, under quite general conditions, the density of occupied patches, and the densities of suitable and occupied patches, for the respective models, have approximate normal distributions. Our results not only provide us with estimates for the means and variances that are valid at all stages in the evolution of the population, but also provide a tool for fitting the models to real metapopulations. We discuss the effect of immigration and habitat dynamics on metapopulations, showing thatmainland-like patches heavily influence metapopulation persistence, and we argue for adopting measures to increase connectivity between this large patch and the other island-like patches. We illustrate our results with specific reference to examples of populations of butterfly and the grasshopperBryodema tuberculata.     

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

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

Extinction dynamics in mainland-island metapopulations: an N-patch stochastic model

A generalization of the well-known Levins’ model of metapopulations is studied. The generalization consists of (i) the introduction of immigration from a mainland, and (ii) assuming the dynamics is stochastic, rather than deterministic. A master equation, for the probability that n of the patches are occupied, is derived and the stationary probability P _s (n), together with the mean and higher moments in the stationary state, determined. The time-dependence of the probability distribution is also studied: through a Gaussian approximation for general n when the boundary at n = 0 has little effect, and by calculating P(0, t), the probability that no patches are occupied at time t, by using a linearization procedure. These analytic calculations are supplemented by carrying out numerical solutions of the master equation and simulations of the stochastic process. The various approaches are in very good agreement with each other. This allows us to use the forms for P _s 0) and P(0, t) in the linearization approximation as a basis for calculating the mean time for a metapopulation to become extinct. We give an analytical expression for the mean time to extinction derived within a mean field approach. We devise a simple method to apply our mean field approach even to complex patch networks in realistic model metapopulations. After studying two spatially extended versions of this nonspatial metapopulation model—a lattice metapopulation model and a spatially realistic model—we conclude that our analytical formula for the mean extinction time is generally applicable to those metapopulations which are really endangered, where extinction dynamics dominates over local colonization processes. The time evolution and, in particular, the scope of our analytical results, are studied by comparing these different models with the analytical approach for various values of the parameters: the rates of immigration from the mainland, the rates of colonization and extinction, and the number of patches making up the metapopulation.     

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.     

How Levins’ dynamics emerges from a Ricker metapopulation model

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

Synchronisation and stability in river metapopulation networks

Spatial structure in landscapes impacts population stability. Two linked components of stability have large consequences for persistence: first, statistical stability as the lack of temporal fluctuations; second, synchronisation as an aspect of dynamic stability, which erodes metapopulation rescue effects. Here, we determine the influence of river network structure on the stability of riverine metapopulations. We introduce an approach that converts river networks to metapopulation networks, and analytically show how fluctuation magnitude is influenced by interaction structure. We show that river metapopulation complexity (in terms of branching prevalence) has nonlinear dampening effects on population fluctuations, and can also buffer against synchronisation. We conclude by showing that river transects generally increase synchronisation, while the spatial scale of interaction has nonlinear effects on synchronised dynamics. Our results indicate that this dual stability – conferred by fluctuation and synchronisation dampening – emerges from interaction structure in rivers, and this may strongly influence the persistence of river metapopulations.     

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

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

A Stochastic Metapopulation Model with Variability in Patch Size and Position

Analytically tractable metapopulation models usually assume that every patch is identical, which limits their application to real metapopulations. We describe a new single species model of metapopulation dynamics that allows variation in patch size and position. The state of the metapopulation is defined by the presence or absence of the species in each patch. For a system of n patches, this gives 2ⁿ possible states. We show how to construct and analyse a matrix describing transitions between all possible states by first constructing separate extinction and colonisation matrices. We illustrate the model′s application to metapopulations by considering an example of malleefowl, Leipoa ocellata, in southern Australia, and calculate extinction probabilities and quasi-stationary distributions. We investigate the relative importance of modelling the particular arrangement of patches and the variation in patch sizes for this metapopulation and we use the model to examine the effects of further habitat loss on extinction probabilities.     

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.     

Structured models of metapopulation dynamics

I develop models of metapopulation dynamics that describe changes in the numbers of individuals within patches. These models are analogous to structured population models, with patches playing the role of individuals. Single species models which do not include the effect of immigration on local population dynamics of occupied patches typically lead to a unique equilibrium. The models can be used to study the distributions of numbers of individuals among patches, showing that both metapopulations with local outbreaks and metapopulations without outbreaks can occur in systems with no underlying environmental variability. Distributions of local population sizes (in occupied patches) can vary independently of the total population size, so both patterns of distributions of local population sizes are compatible with either rare or common species. Models which include the effect of immigration on local population dynamics can lead to two positive equilibria, one stable and one unstable, the latter representing a threshold between regional extinction and persistence.     

RESISTANCE AND VIRULENCE STRUCTURE IN TWO LINUM MARGINALE-MELAMPSORA LINI HOST-PATHOGEN METAPOPULATIONS WITH DIFFERENT MATING SYSTEMS

Different patterns of resistance to six pathotypes of Melampsora lini were detected in 11 populations of Linum marginale distributed across two metapopulations. The two metapopulations (mountains and plains of New South Wales, Australia) differed in the annual cycle of disease development, which barely overlapped, and in the growth cycle and mating system of the host. Host populations in the mountains metapopulation were highly inbred, whereas those on the plains showed appreciable levels of outcrossing. Within each metapopulation there was significant variation among component populations in (1) levels of host resistance to individual pathogen isolates; (2) mean levels of resistance to all six isolates; (3) the number of resistance phenotypes present and the evenness of their distribution within the population; and (4) the average number of pathogen lines to which individual hosts were resistant. A more limited comparison of pathogen populations from the two metapopulations (two from each) found greater similarities in the structure of populations and particular virulence frequencies within, rather than among, the two metapopulations. Differences in host outcrossing rates between the two metapopulations are reflected in marked differences in the overall level of resistance, its partitioning within and among populations, the number and distribution of resistance phenotypes in the two areas, and the level of polymorphism for specific virulence factors in the pathogen, with the plains metapopulation showing consistently higher values. However, these differences were not significant. In general, variation for all parameters was just as great among populations within a metapopulation as between the two metapopulations.     

Reconciling theoretical approaches to stochastic patch-occupancy metapopulation models

There appear to be two different kinds of theoretical results about stochastic patch-occupancy metapopulation models: those recently proposed by Gyllenberg and Silvestrov about metapopulations including a very stable patch, and those by Darroch and Seneta about more general metapopulations. Based on the spectral theory of linear operators, it is shown that the results by Gyllenberg and Silvestrov are a limiting case of those by Darroch and Seneta. Taking the examples proposed by Gyllenberg and Silvestrov as a case study, the application and relevance of these results are discussed, with a particular stress to their bearing on real metapopulations.