Two-patch dispersal-linked compensatory-overcompensatory spatially discrete population models

We study the role of asynchronous and synchronous dispersals on discrete-time two-patch dispersal-linked population models, where the pre-dispersal local patch dynamics are of mixed compensatory and overcompensatory types. Single-species dispersal-linked models behave as single-species single-patch models whenever all pre-dispersal local patch dynamics are compensatory and dispersal is synchronous. However, the dynamics of the corresponding two-patch population model connected by asynchronous dispersal depends on the dispersal rates. The species goes extinct on at least one patch when the asynchronous dispersal rates are high, while it persists when the rates are low. We use numerical simulations to show that in both synchronous and asynchronous mixed compensatory and overcompensatory systems, symmetric and asymmetric dispersals can control and impede the onset of cyclic population oscillations via period-doubling reversal bifurcations. Also, we show that in mixed systems both asynchronous and synchronous dispersals are capable of altering the pre-dispersal local patch dynamics from overcompensatory to compensatory dynamics. Dispersal-linked population models with 'unstructured' overcompensatory pre-dispersal local dynamics connected by synchronous dispersal can generate multiple attractors with fractal basin boundaries. However, mixed compensatory and overcompensatory systems appear to exhibit single attractors and not coexisting (multiple) attractors.     

Interplay between local dynamics and dispersal in discrete-time metapopulation models

The effects of synchronous dispersal on discrete-time metapopulation dynamics with local (patch) dynamics of the same (compensatory or overcompensatory) or mixed (compensatory and overcompensatory) types are explored. Single-species metapopulation models behave as single-species single-patch models, whenever all local patches are governed by compensatory dynamics. Dispersal gives rise to multiple attractors with complex basin structures, whenever some local patches are under overcompensatory dynamics. In mixed systems, dispersal is capable of altering the local dynamics from compensatory to overcompensatory dynamics and vice versa. Examples are provided of metapopulation models supporting multiple attractors with intermingled basins of attraction.     

Dispersal effects on a discrete two-patch model for plant-insect interactions

A two-patch discrete time plant-insect model coupled through insect dispersal is studied. The model is based on three different phases: Plant growth is followed by the dispersal of insects followed by insect attacks. Our objective is to understand how different intensities of dispersal impact both local and global population dynamics of the two-patch model. Special attention is paid to two situations: When the single-patch model (i.e., in the absence of dispersal) is permanent and when the single-patch model exhibits Allee-like effects. The existence and stability of synchronous and asynchronous dynamics between two patches is explored. If the single-patch system is permanent, the permanence of the system in two patches is destroyed by extremely large dispersals and large attacking rates of insects, thus creating multiple attractors. If the single-patch model exhibits Allee-like effects, analytical and numerical results indicate that small intensity of dispersals can generate source-sink dynamics between two patches, while intermediate intensity of dispersals promote the extinction of insects in both patches for certain parameter ranges. Our study suggests a possible biology control strategy to stop the invasion of a pest by controlling its migration between patches.     

Density-dependent migration and synchronism in metapopulations

A spatially explicit metapopulation model with density-dependent dispersal is proposed in order to study the stability of synchronous dynamics. A stability criterion is obtained based on the computation of the transversal Liapunov number of attractors on the synchronous invariant manifold. We examine in detail a special case of density-dependent dispersal rule where migration does not occur if the patch density is below a certain critical density, while the fraction of individuals that migrate to other patches is kept constant if the patch density is above the threshold level. Comparisons with density-independent migration models indicate that this simple density-dependent dispersal mechanism reduces the stability of synchronous dynamics. We were able to quantify exactly this loss of stability through the frequency that synchronous trajectories are above the critical density.     

The dynamics of two diffusively coupled predator-prey populations

I analyze the dynamics of predator and prey populations living in two patches. Within a patch the prey grow logistically and the predators have a Holling type II functional response. The two patches are coupled through predator migration. The system can be interpreted as a simple predator-prey metapopulation or as a spatially explicit predator-prey system. Asynchronous local dynamics are presumed by metapopulation theory. The main question I address is when synchronous and when asynchronous dynamics arise. Contrary to biological intuition, for very small migration rates the oscillations always synchronize. For intermediate migration rates the synchronous oscillations are unstable and I found periodic, quasi-periodic, and intermittently chaotic attractors with asynchronous dynamics. For large predator migration rates, attractors in the form of equilibria or limit cycles exist in which one of the patches contains no prey. The dynamical behavior of the system is described using bifurcation diagrams. The model shows that spatial predator-prey populations can be regulated through the interplay of local dynamics and migration.     

Dynamics of populations with individual variation in dispersal on bounded domains

Most classical models for the movement of organisms assume that all individuals have the same patterns and rates of movement (for example, diffusion with a fixed diffusion coefficient) but there is empirical evidence that movement rates and patterns may vary among different individuals. A simple way to capture variation in dispersal that has been suggested in the ecological literature is to allow individuals to switch between two distinct dispersal modes. We study models for populations whose members can switch between two different nonzero rates of diffusion and whose local population dynamics are subject to density dependence of logistic type. The resulting models are reaction–diffusion systems that can be cooperative at some population densities and competitive at others. We assume that the focal population inhabits a bounded region and study how its overall dynamics depend on the parameters describing switching rates and local population dynamics. (Traveling waves and spread rates have been studied for similar models in the context of biological invasions.) The analytic methods include ideas and results from reaction–diffusion theory, semi-dynamical systems, and bifurcation/continuation theory.     

Population synchrony and stability in environmentally forced metacommunities

A general prediction from simple metapopulation models is that spatially synchronized forcing can spatially synchronize population dynamics and destabilize metapopulations. In contrast, spatially asynchronous forcing is predicted to decrease population synchrony and promote temporal stability and population persistence, especially in the presence of dispersal. Only recently have studies begun to experimentally address these predictions. Moreover, few studies have experimentally examined how such processes operate in the context of competition communities. Stabilizing processes may continue to operate when placed within a metacommunity context with multiple competing consumers but only at low to intermediate levels of dispersal. High dispersal rates can reverse these predictions and lead to destabilization. We tested this under controlled conditions using an experimental aquatic system composed of three competing species of zooplankton. Metacommunities experienced different levels of dispersal and environmental forcing in the form of spatially synchronous or asynchronous pH perturbations. We found support that dispersal can have contrasting effects on population stability depending on the degree to which population dynamics were synchronized in space. Dispersal under synchronous forcing or no forcing had either neutral of positive effects on spatial population synchrony of all three zooplankton species. In these treatments, dispersal reduced population stability at the local and metapopulation levels for two of three species. In contrast, asynchronously varying environments reduced population synchrony relative to unforced systems, regardless of dispersal level. In these treatments, dispersal enhanced temporal stability and persistence of populations not by reducing population synchrony but by enhancing population minima and spatial averaging of abundances. High dispersal rates under asynchronous forcing reduced the abundance of one species, consistent with increasing regional competition and general metacommunity theory. However, no effects on its stability or persistence were observed. Our work highlights the context‐dependent effects of dispersal on population dynamics in varying environments.     

Warming induces synchrony and destabilizes experimental pond zooplankton metacommunities

The spatial insurance hypothesis predicts that intermediate rates of dispersal between patches in a metacommunity allow species to track favourable conditions, preserving diversity and stabilizing biomass at local and regional scales. However, theory is unclear as to whether dispersal will provide spatial insurance when environmental conditions are changing directionally. In particular, increased temperatures as a result of climate change are expected to cause synchronous growth or decline across species and communities, and this has the potential to erode the stabilizing compensatory dynamics facilitated by dispersal. Here we report on an experimental test of how dispersal affects the diversity and stability of metacommunities under warming using replicate two‐patch pond zooplankton metacommunities. Initial differences in local community composition and abiotic conditions were established by seeding each patch in the metacommunities with plankton and sediment from one of two natural ponds that differed in water chemistry and species composition. We exposed metacommunities to a 2°C increase in average ambient temperature, crossed with three rates of dispersal (none, intermediate, high). In ambient conditions, intermediate dispersal rates preserved diversity and stabilized metacommunities by promoting spatially asynchronous fluctuations in biomass, especially between local populations of the dominant genus, Ceriodaphnia. However, warming synchronized their populations so that these effects of dispersal were lost. Furthermore, because the stabilizing effect of dispersal was primarily due to asynchronous fluctuations between populations of a single genus, metacommunity biomass was stabilized, but dispersal did not stabilize local community biomass. Our results show that dispersal can preserve diversity and provide stability to metacommunities, but also show that this benefit can be eroded when warming is directional and synchronous across patches of a metacommunity, as is expected with climate warming.     

Resource matching and population dynamics in a two-patch system

We study resource matching – the relationship between resource supply and forager numbers – under conditions of fluctuating population dynamics in a two-patch system. For the inter-patch dispersal we apply the patch-departure rule following the principle of the ideal free distribution: leave the current patch of residence if local conditions are worse than conditions elsewhere on average. We show that such a dispersal rule synchronises cyclic and chaotic local population dynamics, but unlike many other dispersal rules, leaves the underlying population dynamics untouched. We also show that the IFD dispersal rule is not very sensitive to biased information and navigation failures during the dispersal phase. Even under such circumstances we observe a quick process of populations becoming synchronised, even when the population dynamics are chaotic. We conclude that an IFD patch-departure rule represents an ESS dispersal behaviour towards which the dispersal patterns should evolve.     

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.     

Patterns of dispersal and dynamics among habitat patches varying in quality

Both source-sink theory and extensions of optimal foraging theory ("balanced dispersal" theory) address dispersal and population dynamics in landscapes where habitat patches vary in quality. However, studying dispersal mechanisms empirically has proven difficult, and dispersal is rarely tied back to long-term spatial dynamics. We used a manipulable laboratory system consisting of bacteria and protozoa to investigate the ability of source-sink and optimal foraging theories to explain both dispersal and emergent spatial dynamics. Consistent with source-sink models and contrary to balanced dispersal models, there was a consistent net flux of protist individuals from high to low resource patches. However, unlike the simplest source-sink models, intermediate rates of dispersal led to highest abundances in low resource patches. Side experiments found strong density dependence in local population dynamics and differences in average protist body size in high and low resource patches. Parameterization and analysis of a two-patch model showed that high migration from high to low resource patches could have depressed population density in low resource patches, creating pseudosinks. The movement of individuals and biomass from sources to sinks (a form of ecosystem subsidy) resulted in the convergence of body size and population densities in sources and sinks. Our results indicate a need to carefully consider movement patterns and interaction with local dynamics in potential source-sink systems.     

EVOLUTION OF DISPERSAL RATES IN METAPOPULATION MODELS: BRANCHING AND CYCLIC DYNAMICS IN PHENOTYPE SPACE

We study the evolution of dispersal rates in a two patch metapopulation model. The local dynamics in each patch are given by difference equations, which, together with the rate of dispersal between the patches, determine the ecological dynamics of the metapopulation. We assume that phenotypes are given by their dispersal rate. The evolutionary dynamics in phenotype space are determined by invasion exponents, which describe whether a mutant can invade a given resident population. If the resident metapopulation is at a stable equilibrium, then selection on dispersal rates is neutral if the population sizes in the two patches are the same, while selection drives dispersal rates to zero if the local abundances are different. With non-equilibrium metapopulation dynamics, non-zero dispersal rates can be maintained by selection. In this case, and if the patches are ecologically identical, dispersal rates always evolve to values which induce synchronized metapopulation dynamics. If the patches are ecologically different, evolutionary branching into two coexisting dispersal phenotypes can be observed. Such branching can happen repeatedly, leading to polymorphisms with more than two phenotypes. If there is a cost to dispersal, evolutionary cycling in phenotype space can occur due to the dependence of selection pressures on the ecological attractor of the resident population, or because phenotypic branching alternates with the extinction of one of the branches. Our results extend those of Holt and McPeek (1996), and suggest that phenotypic branching is an important evolutionary process. This process may be relevant for sympatric speciation.     

An epidemic model in a patchy environment

An epidemic model is proposed to describe the dynamics of disease spread among patches due to population dispersal. We establish a threshold above which the disease is uniformly persistent and below which disease-free equilibrium is locally attractive, and globally attractive when both susceptible and infective individuals in each patch have the same dispersal rate. Two examples are given to illustrate that the population dispersal plays an important role for the disease spread. The first one shows that the population dispersal can intensify the disease spread if the reproduction number for one patch is large, and can reduce the disease spread if the reproduction numbers for all patches are suitable and the population dispersal rate is strong. The second example indicates that a population dispersal results in the spread of the disease in all patches, even though the disease can not spread in each isolated patch.     

Variation in dispersal mortality and dispersal propensity among individuals: the effects of age, sex and resource availability

1.  Dispersal of individuals between habitat patches depends on both the propensity to emigrate from a patch and the ability to survive inter-patch movement. Environmental factors and individual characteristics have been shown to influence dispersal rates but separating the effects of emigration and dispersal mortality on dispersal can often be difficult. In this study, we use a soil mite laboratory system to investigate factors affecting emigration and dispersal mortality.
2.  We tested the movement of different age groups in two-patch systems with different inter-patch distances. Differences in immigration among age groups were primarily driven by differences in emigration but dispersal mortality was greater for some groups. Immigration declined with increasing inter-patch distance, which was due to increasing dispersal mortality and decreasing emigration.
3.  In a second experiment, we compared the dispersal of recently matured males and females and tested the impact of food availability during the developmental period on their dispersal. Dispersal was found to be male biased but there was no significant sex bias in dispersal mortality. There was some evidence that food availability could affect emigration and dispersal mortality.
4.  These results demonstrate that both emigration and dispersal mortality can be affected by factors such as individual age and resource availability. Understanding these effects is likely to be important for predicting the fitness costs and population consequences of dispersal.     

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.     

Spatial variation in branch size promotes metapopulation persistence in dendritic river networks

1. Despite years of attention, the dynamics of species constrained to disperse within riverine networks are not well captured by existing metapopulation models, which often ignore local dynamics within branches.
2. We develop a modelling framework, based on traditional metapopulation theory, for patch occupancy dynamics subject to local colonisation–extinction dynamics within branches and regional dispersal between branches in size-structured, bifurcating riverine networks. Using this framework, we investigate whether and how spatial variation in branch size affects species persistence for dendritic systems with directional dispersal, including one-way (up- or downstream only) and two-way (both up- and downstream) dispersal.
3. Variation in branch size generally promotes species persistence more obviously at higher relative extinction rate, suggesting that previous studies ignoring differences in branch size in real riverine systems might overestimate species extinction risk.
4. Two-way dispersal is not always superior to one-way dispersal as a strategy for metapopulation persistence especially at high relative extinction rate. The type of dispersal that maximises species persistence is determined by the hierarchical level of the largest, and hence most influential, branch within the network. When considering the interactive effects of up- and downstream dispersal, we find that moderate upstream-biased dispersal maximises metapopulation viability, mediated by spatial branch arrangement.
5. Overall, these results suggest that both branch-size variation and species traits interact to determine species persistence, theoretically demonstrating the ecological significance of their interplay.

     

Expansion or extinction: deterministic and stochastic two-patch models with Allee effects

We investigate the impact of Allee effect and dispersal on the long-term evolution of a population in a patchy environment. Our main focus is on whether a population already established in one patch either successfully invades an adjacent empty patch or undergoes a global extinction. Our study is based on the combination of analytical and numerical results for both a deterministic two-patch model and a stochastic counterpart. The deterministic model has either two, three or four attractors. The existence of a regime with exactly three attractors only appears when patches have distinct Allee thresholds. In the presence of weak dispersal, the analysis of the deterministic model shows that a high-density and a low-density populations can coexist at equilibrium in nearby patches, whereas the analysis of the stochastic model indicates that this equilibrium is metastable, thus leading after a large random time to either a global expansion or a global extinction. Up to some critical dispersal, increasing the intensity of the interactions leads to an increase of both the basin of attraction of the global extinction and the basin of attraction of the global expansion. Above this threshold, for both the deterministic and the stochastic models, the patches tend to synchronize as the intensity of the dispersal increases. This results in either a global expansion or a global extinction. For the deterministic model, there are only two attractors, while the stochastic model no longer exhibits a metastable behavior. In the presence of strong dispersal, the limiting behavior is entirely determined by the value of the Allee thresholds as the global population size in the deterministic and the stochastic models evolves as dictated by their single-patch counterparts. For all values of the dispersal parameter, Allee effects promote global extinction in terms of an expansion of the basin of attraction of the extinction equilibrium for the deterministic model and an increase of the probability of extinction for the stochastic model.     

Genetic influences on social dominance: cow wars

Analyzing population dynamics in changing habitats is a prerequisite for population dynamics forecasting. The recent development of metapopulation modeling allows the estimation of dispersal kernels based on the colonization pattern but the accuracy of these estimates compared with direct estimates of the seed dispersal kernel has rarely been assessed. In this study, we used recent genetic methods based on parentage analysis (spatially explicit mating models) to estimate seed and pollen dispersal kernels as well as seed and pollen immigration in fragmented urban populations of the plant species Crepis sancta with contrasting patch dynamics. Using two independent networks, we documented substantial seed immigration and a highly restricted dispersal kernel. Moreover, immigration heterogeneity among networks was consistent with previously reported metapopulation dynamics, showing that colonization was mainly due to external colonization in the first network (propagule rain) and local colonization in the second network. We concluded that the differences in urban patch dynamics are mainly due to seed immigration heterogeneity, highlighting the importance of external population source in the spatio-temporal dynamics of plants in a fragmented landscape. The results show that indirect and direct methods were qualitatively consistent, providing a proper interpretation of indirect estimates. This study provides attempts to link genetic and demographic methods and show that patch occupancy models may provide simple methods for analyzing population dynamics in heterogeneous landscapes in the context of global change.     

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.