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.

     

Host–parasitoid dynamics in a fragmented landscape: Holly trees,holly leaf miners and their parasitoids

Many species inhabit fragmented landscapes, where units of resource have a patchy spatial distribution. While numerous studies have investigated how the incidence and dynamics of individual species are affected by the spatial configuration and landscape context of habitat patches, fewer studies have investigated the dynamics of multiple interacting resource and consumer species in patchy landscapes. We describe a model system for investigating host–parasitoid dynamics in a patchy landscape: a network of 166 holly trees, a specialised herbivore of holly (the leaf miner, Phytomyza ilicis (Curtis, 1948)), and its suite of parasitoids. We documented patch occupancy by P. ilicis, its density within patches, and levels of parasitism over a 6-year period, and manipulated patch occupancy by creating artificially vacant habitat patches. Essentially all patches were occupied by the herbivore in each year, suggesting that metapopulation dynamics are unlikely to occur in this system. The main determinants of densities for P. ilicis and its parasitoids were resource availability (patch size and host density, respectively). While P. ilicis is apparently not restricted by the spatial distribution of resources, densities of its parasitoids showed a weaker positive relationship with host density in more isolated patches. In patches where local extinctions were generated experimentally, P. ilicis densities and levels of parasitism recovered to pre-manipulation levels within a single generation. Furthermore, patch isolation did not significantly affect re-colonisation by hosts or parasitoids. Analysing the data at a variety of spatial scales indicates that the balance between local demography and dispersal may vary depending on the scale at which patches are defined. Taken together, our results suggest that the host and its parasitoids have dispersal abilities that exceed typical inter-patch distances. Patch dynamics are thus largely governed by dispersal rather than within-patch demography, although the role of demography is higher in larger patches.     

Heterogeneity in local density allows a positive evolutionary relationship between self‐fertilisation and dispersal

Despite empirical evidence for a positive relationship between dispersal and self‐fertilization (selfing), theoretical work predicts that these traits should always be negatively correlated, and the Good Coloniser Syndrome of high dispersal and selfing (Cf. Baker's Law) should not evolve. Critically, previous work assumes that adult density is spatiotemporally homogeneous, so selfing results in identical offspring production for all patches, eliminating the benefit of dispersal for escaping from local resource competition. We investigate the joint evolution of dispersal and selfing in a demographically structured metapopulation model where local density is spatiotemporally heterogeneous due to extinction‐recolonization dynamics. Selfing alleviates outcrossing failure due to low local density (an Allee effect) while dispersal alleviates competition through dispersal of propagules from high‐ to low‐density patches. Because local density is spatiotemporally heterogeneous in our model, selfing does not eliminate heterogeneity in competition, so dispersal remains beneficial even under full selfing. Hence the Good Coloniser Syndrome is evolutionarily stable under a broad range of conditions, and both negative and positive relationships between dispersal and selfing are possible, depending on the environment. Our model thus accommodates positive empirical relationships between dispersal and selfing not predicted by previous theoretical work and provides additional explanations for negative relationships.     

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.     

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.     

On the evolution of specialization with a mechanistic underpinning in structured metapopulations

We analyze the evolution of specialization in resource utilization in a discrete-time metapopulation model using the adaptive dynamics approach. The local dynamics in the metapopulation are based on the Beverton-Holt model with mechanistic underpinnings. The consumer faces a trade-off in the abilities to consume two resources that are spatially heterogeneously distributed to patches that are prone to local catastrophes. We explore the factors favoring the spread of generalist or specialist strategies. Increasing fecundity or decreasing catastrophe probability favors the spread of the generalist strategy and increasing environmental heterogeneity enlarges the parameter domain where the evolutionary branching is possible. When there are no catastrophes, increasing emigration diminishes the parameter domain where the evolutionary branching may occur. Otherwise, the effect of emigration on evolutionary dynamics is non-monotonous: both small and large values of emigration probability favor the spread of the specialist strategies whereas the parameter domain where evolutionary branching may occur is largest when the emigration probability has intermediate values. We compare how different forms of spatial heterogeneity and different models of local growth affect the evolutionary dynamics. We show that even small changes in the resource dynamics may have outstanding evolutionary effects to the consumers.     

Modelling single nucleotide effects in phosphoglucose isomerase on dispersal in the Glanville fritillary butterfly: coupling of ecological and evolutionary dynamics

Dispersal comprises a complex life-history syndrome that influences the demographic dynamics of especially those species that live in fragmented landscapes, the structure of which may in turn be expected to impose selection on dispersal. We have constructed an individual-based evolutionary sexual model of dispersal for species occurring as metapopulations in habitat patch networks. The model assumes correlated random walk dispersal with edge-mediated behaviour (habitat selection) and spatially correlated stochastic local dynamics. The model is parametrized with extensive data for the Glanville fritillary butterfly. Based on empirical results for a single nucleotide polymorphism (SNP) in the phosphoglucose isomerase (Pgi) gene, we assume that dispersal rate in the landscape matrix, fecundity and survival are affected by a locus with two alleles, A and C, individuals with the C allele being more mobile. The model was successfully tested with two independent empirical datasets on spatial variation in Pgi allele frequency. First, at the level of local populations, the frequency of the C allele is the highest in newly established isolated populations and the lowest in old isolated populations. Second, at the level of sub-networks with dissimilar numbers and connectivities of patches, the frequency of C increases with decreasing network size and hence with decreasing average metapopulation size. The frequency of C is the highest in landscapes where local extinction risk is high and where there are abundant opportunities to establish new populations. Our results indicate that the strength of the coupling of the ecological and evolutionary dynamics depends on the spatial scale and is asymmetric, demographic dynamics having a greater immediate impact on genetic dynamics than vice versa.     

Regional and local scale metapopulation dynamics in the interaction between Callosobruchus maculatus and Anisopteromalus calandrae

Habitat structure increases the persistence of many extinction‐prone resource–consumer interactions. Metapopulation theory is one of the leading approaches currently used to explain why local, ephemeral populations persist at a regional scale. Central to the metapopulation concept is the amount of dispersal occurring between patches, too much or too little can result in regional extinction. In this study, the role of dispersal on the metapopulation dynamics of an over‐exploitative host–parasitoid interaction is assessed. In the absence of the parasitoid the highly vagile bruchid, Callosobruchus maculatus, can maintain a similar population size regardless of the permeability of the inter‐patch matrix and exhibits strong negative density‐dependence. After the introduction of the parasitoid the size of the bruchid population decreases with a corresponding increase in the occurrence of empty patches. In this case, limiting the dispersal of both species decouples the interaction to a greater extent and results in larger regional bruchid populations. Given the disparity between the dispersal rates of the two species, it is proposed that the more dispersive host benefits from the reduction in landscape permeability by increasing the opportunity to colonise empty patches and rescue extinction prone populations. Associated with the introduction of the parasitoid is a shift in the strength of density‐dependence as the population moves from bottom–up towards top–down regulation. The importance of local and regional scale measurements is apparent when the role of individual patches on regional dynamics is considered. By only taking regional dynamics into account the importance of dispersal regime on local dynamics is overlooked. Similarly, when local dynamics were examined, patches were found to have different influences on regional dynamics depending on dispersal regime and patch location.     

The effects of density-dependence and local dispersal in individual-based stochastic metapopulations

We describe simulation models for metapopulations of individual‐based random walk populations with local dispersal on a coupled map lattice. The models were used to assess the factors determining persistence time, in particular the incidence of density‐dependence required for long‐term persistence of a temporally and spatially stochastic metapopulation, the extent of persistence possible in the absence of density‐dependence, and the factors that affect this.
Metapopulation persistence depended on the overall rate of increase of the metapopulation mean. This was maximised by (in order of importance) high mean and variance in the local rate of increase, high dispersal rates (20% or more of individuals dispersing each generation), large lattice size, and large dispersal range (to at least 24 neighbouring subpopulations).
With density‐dependence, the emergent dynamics of the metapopulation mean following global perturbation (reduction in density) resembled those of the logistic growth model. However, the overall metapopulation rate of increase and equilibrium level bore no resemblance to those of the subpopulations: rate of increase was higher (negative mean local rates of increase may give positive overall growth), and equilibrium mean metapopulation density was well below the local carrying capacity. This highlights the need to sample populations at an appropriate scale when seeking to understand regulatory mechanisms.
Metapopulations with the strongest tendency to grow gave the highest equilibrium mean density, the highest incidence of density‐dependence, and the longest persistence time. However, long‐term persistence with low average density and very low incidence of density‐dependence was possible on a sufficiently large lattice. For example, with 40×40 subpopulations, mean metapopulation persistence time was around 10⁴ generations, with mean subpopulation size of 2% of the carrying capacity, and local density‐dependence acting just once every 2500 generations on average. Metapopulation processes may explain our inability to detect density‐dependence in many real populations, and may also play an important part in the persistence of rare species.     

The interactive effects of competition and predation risk on dispersal in an insect

Dispersal dynamics have significant consequences for ecological and evolutionary processes. Previous work has demonstrated that dispersal can be context-dependent. However, factors affecting dispersal are typically considered in isolation, despite the probability that individuals make dispersal decisions in response to multiple, possibly interacting factors. We examined whether two ecological factors, predation risk and intraspecific competition, have interactive effects on dispersal dynamics. We performed a factorial experiment in mesocosms using backswimmers (Notonecta undulata), flight-capable, semi-aquatic insects. Emigration rates increased with density, and increased with predation risk at intermediate densities; however, predation had minimal effects on emigration at high and low densities. Our results indicate that factorial experiments may be required to understand dispersal dynamics under realistic ecological conditions.     

Support for a metapopulation structure among mammals

• 1 The metapopulation metaphor is increasingly used to explain the spatial dynamics of animal populations. However, metapopulation structure is difficult to identify in long‐lived species that are widely distributed in stochastic environments, where they can resist extinctions. The literature on mammals may not provide supporting evidence for classic metapopulation dynamics, which call for the availability of discrete habitat patches, asynchrony in local population dynamics, evidence for extinction and colonization processes, and dispersal between local populations.
• 2 Empirical evidence for metapopulation structure among mammals may exist when applying more lenient criteria. To meet these criteria, mammals should live in landscapes as discrete local breeding populations, and their demography should be asynchronous.
• 3 We examined the literature for empirical evidence in support of the classical criteria set by Hanski (1999 ), and for the more lenient subset of criteria proposed by Elmhagen & Angerbjörn (2001 ). We suggest circumstances where metapopulation theory could be important in understanding population processes in mammals of different body sizes.
• 4 The patchy distribution of large (>100 kg) mammals and dispersal often motivate inferences in support of a metapopulation structure. Published studies seldom address the full suite of classical criteria. However, studies on small mammals are more likely to record classic metapopulation criteria than those on large mammals. The slow turnover rate that is typical for medium‐sized and large mammals apparently makes it difficult to identify a metapopulation structure during studies of short duration.
• 5 To identify a metapopulation structure, studies should combine the criteria set by Hanski (1999 ) and Elmhagen & Angerbjörn (2001 ). Mammals frequently live in fragmented landscapes, and processes involved in the maintenance of a metapopulation structure should be considered in conservation planning and management.

     

Density-Dependent Natal Dispersal Patterns in a Leopard Population Recovering from Over-Harvest

Natal dispersal enables population connectivity, gene flow and metapopulation dynamics. In polygynous mammals, dispersal is typically male-biased. Classically, the ‘mate competition’, ‘resource competition’ and ‘resident fitness’ hypotheses predict density-dependent dispersal patterns, while the ‘inbreeding avoidance’ hypothesis posits density-independent dispersal. In a leopard (Panthera pardus) population recovering from over-harvest, we investigated the effect of sex, population density and prey biomass, on age of natal dispersal, distance dispersed, probability of emigration and dispersal success. Over an 11-year period, we tracked 35 subadult leopards using VHF and GPS telemetry. Subadult leopards initiated dispersal at 13.6 ± 0.4 months. Age at commencement of dispersal was positively density-dependent. Although males (11.0 ± 2.5 km) generally dispersed further than females (2.7 ± 0.4 km), some males exhibited opportunistic philopatry when the population was below capacity. All 13 females were philopatric, while 12 of 22 males emigrated. Male dispersal distance and emigration probability followed a quadratic relationship with population density, whereas female dispersal distance was inversely density-dependent. Eight of 12 known-fate females and 5 of 12 known-fate male leopards were successful in settling. Dispersal success did not vary with population density, prey biomass, and for males, neither between dispersal strategies (philopatry vs. emigration). Females formed matrilineal kin clusters, supporting the resident fitness hypothesis. Conversely, mate competition appeared the main driver for male leopard dispersal. We demonstrate that dispersal patterns changed over time, i.e. as the leopard population density increased. We conclude that conservation interventions that facilitated local demographic recovery in the study area also restored dispersal patterns disrupted by unsustainable harvesting, and that this indirectly improved connectivity among leopard populations over a larger landscape.     

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.     

Prey (Moina macrocopa) population density drives emigration rate of its predator (Trichocorixa verticalis) in a rock-pool metacommunity

Dispersal connects spatially separated local food webs at a larger, metacommunity scale, and as a result, dispersal may both influence and be influenced by local food-web dynamics. Here, I focused on a rock-pool metacommunity and used a combination of observational, experimental, and theoretical approaches to explore the role of local prey (Moina macrocopa) density on the rate of emigration by their predator (Trichocorixa verticalis) and in turn, the effect of predator emigration on the per capita predation rate experienced by local prey populations. A lab feeding experiment quantified predation rates, demonstrating that indeed adult T. verticalis are voracious predators of M. macrocopa. M. macrocopa densities vary over five orders of magnitude across both space and time in rock pools, and a mesocosm experiment showed that this variation significantly influences T. verticalis emigration: predators emigrated more rapidly when prey were in lower densities. Finally, computer simulations demonstrated that this pattern of dispersal by T. verticalis has the potential to relieve local M. macrocopa populations from predation when the prey are at low densities, thereby reducing the likelihood that local M. macrocopa populations will be driven extinct by predation from T. verticalis.     

Evolutionarily Stable Dispersal Rates Do Not Always Increase with Local Extinction Rates

Earlier models on the evolution of dispersal have suggested that evolutionarily stable dispersal rates should increase with the frequency of local extinctions. Most metapopulation models assume site saturation (i.e., no local population dynamics), yet the majority of species distributed as metapopulations rarely attain carrying capacity in all occupied patches. In this article, we relax this assumption and examine the evolutionarily stable dispersal rate under nonsaturated but still competitive demographic conditions. Contrary to previous predictions, we show that increasing local extinction rates may allow decreasing dispersal rates to evolve.     

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.     

The evolution of an 'intelligent' dispersal strategy: biased, correlated random walks in patchy landscapes

Theoretical work exploring dispersal evolution focuses on the emigration rate of individuals and typically assumes that movement occurs either at random to any other patch or to one of the nearest‐neighbour patches. There is a lack of work exploring the process by which individuals move between patches, and how this process evolves. This is of concern because any organism that can exert control over dispersal direction can potentially evolve efficiencies in locating patches, and the process by which individuals find new patches will potentially have major effects on metapopulation dynamics and gene flow. Here, we take an initial step towards filling this knowledge gap. To do this we constructed a continuous space population model, in which individuals each carry heritable trait values that specify the characteristics of the biased correlated random walk they use to disperse from their natal patch. We explore how the evolution of the random walk depends upon the cost of dispersal, the density of patches in the landscape, and the emigration rate. The clearest result is that highly correlated walks always evolved (individuals tended to disperse in relatively straight lines from their natal patch), reflecting the efficiency of straight‐line movement. In our models, more costly dispersal resulted in walks with higher correlation between successive steps. However, the exact walk that evolved also depended upon the density of suitable habitat patches, with low density habitat evolving more biased walks (individuals which orient towards suitable habitat at quite large distances from that habitat). Thus, low density habitat will tend to develop individuals which disperse efficiently between adjacent habitat patches but which only rarely disperse to more distant patches; a result that has clear implications for metapopulation theory. Hence, an understanding of the movement behaviour of dispersing individuals is critical for robust long‐term predictions of population dynamics in fragmented landscapes.     

Finite metapopulation models with density-dependent migration and stochastic local dynamics

The effects of small density-dependent migration on the dynamics of a metapopulation are studied in a model with stochastic local dynamics. We use a diffusion approximation to study how changes in the migration rate and habitat occupancy affect the rates of local colonization and extinction. If the emigration rate increases or if the immigration rate decreases with local population size, a positive expected rate of change in habitat occupancy is found for a greater range of habitat occupancies than when the migration is density-independent. In contrast, the reverse patterns of density dependence in respective emigration and immigration reduce the range of habitat occupancies where the metapopulation will be viable. This occurs because density-dependent migration strongly influences both the establishment and rescue effects in the local dynamics of metapopulations.     

The role of density-dependent dispersal in source-sink dynamics

I investigate two aspects of source-sink theory that have hitherto received little attention: density-dependent dispersal and the cost of dispersal to sources. The cost arises because emigration reduces the per capita growth rate of sources, thus predisposing them to extinction. I show that source-sink persistence depends critically on the interplay between these two factors. When the emigration rate increases with abundance at an accelerating rate, dispersal costs to sources is the lowest and risk of source-sink extinction the least. When the emigration rate increases with abundance at a decelerating rate, dispersal costs to sources is the highest and the risk of source-sink extinction the greatest. Density-independent emigration has an intermediate effect. Thus, density-dependent dispersal per se does not increase or decrease source-sink persistence relative to density-independent dispersal. The exact mode of dispersal is crucial. A key point to appreciate is that these effects of dispersal on source-sink extinction arise from the temporal density-dependence that dispersal induces in the per capita growth rates of source and sink populations. Temporal density-dependence due to dispersal is beneficial at low abundances because it rescues sinks from extinction, and detrimental at high abundances because it drives otherwise viable sources to extinction. These results are robust to the nature of population dynamics in the sink, whether exponential or logistic. They provide a means of assessing the relative costs and benefits of preserving sink habitats given three biological parameters.     

The evolution of dispersal – the importance of information about population density and habitat characteristics

The evolution of mobility patterns and dispersal strategies depend on different population, habitat and life history characteristics. The ability to perceive and make use of information about the surrounding environment for dispersal decisions will also differ between organisms. To investigate the evolutionary consequences of such differences, we have used a simulation model with nearest-neighbour dispersal in a metapopulation to study how variation in the ability to obtain and make use of information about habitat quality and conspecific density affects the evolution of dispersal strategies. We found a rather strong influence of variation in information on the overall rate of dispersal in a metapopulation. The highest emigration rate evolved in organisms with no information about either density or habitat quality and the lowest rate was found in organisms with information about both the natal and the neighbouring patches. For organisms that can make use of information about conspecific density, positively density-dependent dispersal evolved in the majority of cases, with the strongest density dependence occurring when an individual only has information about density in the natal patch. However, we also identified situations, involving strong local population fluctuations and frequent local extinctions, where negatively density-dependent dispersal evolved.