Evolution of cross-diffusion and self-diffusion

This article is concerned with the evolution of certain types of density-dependent dispersal strategy in the context of two competing species with identical population dynamics and same random dispersal rates. Such density-dependent movement, often referred to as cross-diffusion and self-diffusion, assumes that the movement rate of each species depends on the density of both species and that the transition probability from one place to its neighbourhood depends solely on the arrival spot (independent of the departure spot). Our results suggest that for a one-dimensional homogeneous habitat, if the gradients of two cross- and self-diffusion coefficients have the same direction, the species with the smaller gradient will win, i.e. the dispersal strategy with the smaller gradient of cross- and self-diffusion coefficient will evolve. In particular, it suggests that the species with constant cross- and self-diffusion coefficients may have competitive advantage over species with non-constant cross- and self-diffusion coefficients. However, if the two gradients have opposite directions, neither of the two dispersal strategies wins as these two species can coexist.     

Evolution of conditional dispersal: evolutionarily stable strategies in spatial models

We consider a two-species competition model in which the species have the same population dynamics but different dispersal strategies. Both species disperse by a combination of random diffusion and advection along environmental gradients, with the same random dispersal rates but different advection coefficients. Regarding these advection coefficients as movement strategies of the species, we investigate their course of evolution. By applying invasion analysis we find that if the spatial environmental variation is less than a critical value, there is a unique evolutionarily singular strategy, which is also evolutionarily stable. If the spatial environmental variation exceeds the critical value, there can be three or more evolutionarily singular strategies, one of which is not evolutionarily stable. Our results suggest that the evolution of conditional dispersal of organisms depends upon the spatial heterogeneity of the environment in a subtle way.     

The Evolution of Conditional Dispersal Strategies in Spatially Heterogeneous Habitats

To understand the evolution of dispersal, we study a Lotka–Volterra reaction–diffusion–advection model for two competing species in a heterogeneous environment. The two species are assumed to be identical except for their dispersal strategies: both species disperse by random diffusion and advection along environmental gradients, but with slightly different random dispersal or advection rates. Two new phenomena are found for one-dimensional habitats and monotone intrinsic growth rates: (i) If both species disperse only by random diffusion, i.e., no advection, it was well known that the slower diffuser always wins. We show that if both species have the same advection rate which is suitably large, the faster dispersal will evolve; (ii) If both species have the same random dispersal rate, it was known that the species with a little advection along the resource gradient always wins, provided that the other species is a pure random disperser and the habitat is convex. We show that if both species have the same random dispersal rate and both also have suitably large advection rates, the species with a little smaller advection rate always wins. Implications of these results for the habitat choices of species will be discussed. Some future directions and open problems will be addressed.     

Evolutionarily Stable and Convergent Stable Strategies in Reaction–Diffusion Models for Conditional Dispersal

We consider a mathematical model of two competing species for the evolution of conditional dispersal in a spatially varying, but temporally constant environment. Two species are different only in their dispersal strategies, which are a combination of random dispersal and biased movement upward along the resource gradient. In the absence of biased movement or advection, Hastings showed that the mutant can invade when rare if and only if it has smaller random dispersal rate than the resident. When there is a small amount of biased movement or advection, we show that there is a positive random dispersal rate that is both locally evolutionarily stable and convergent stable. Our analysis of the model suggests that a balanced combination of random and biased movement might be a better habitat selection strategy for populations.     

The evolution of density-dependent dispersal

Despite a large body of empirical evidence suggesting that the dispersal rates of many species depend on population density, most metapopulation models assume a density-independent rate of dispersal. Similarly, studies investigating the evolution of dispersal have concentrated almost exclusively on density-independent rates of dispersal. We develop a model that allows density-dependent dispersal strategies to evolve. Our results demonstrate that a density-dependent dispersal strategy almost always evolves and that the form of the relationship depends on reproductive rate, type of competition, size of subpopulation equilibrium densities and cost of dispersal. We suggest that future metapopulation models should account for density-dependent dispersal     

Evolution of conditional dispersal: a reaction–diffusion–advection model

To study evolution of conditional dispersal, a Lotka-Volterra reaction-diffusion-advection model for two competing species in a heterogeneous environment is proposed and investigated. The two species are assumed to be identical except their dispersal strategies: both species disperse by random diffusion and advection along environmental gradients, but one species has stronger biased movement (i.e., advection along the environmental gradients) than the other one. It is shown that at least two scenarios can occur: if only one species has a strong tendency to move upward the environmental gradients, the two species can coexist since one species mainly pursues resources at places of locally most favorable environments while the other relies on resources from other parts of the habitat; if both species have such strong biased movements, it can lead to overcrowding of the whole population at places of locally most favorable environments, which causes the extinction of the species with stronger biased movement. These results provide a new mechanism for the coexistence of competing species, and they also imply that selection is against excessive advection along environmental gradients, and an intermediate biased movement rate may evolve.     

Interspecific Competition, Environmental Gradients, Gene Flow, and the Coevolution of Species' Borders

Darwin viewed species range limits as chiefly determined by an interplay between the abiotic environment and interspecific interactions. Haldane argued that species' ranges could be set intraspecifically when gene flow from a species' populous center overwhelms local adaptation at the periphery. Recently, Kirkpatrick and Barton have modeled Haldane's process with a quantitative genetic model that combines density-dependent local population growth with dispersal and gene flow across a linear environmental gradient in optimum phenotype. To address Darwin's ideas, we have extended the Kirkpatrick and Barton model to include interspecific competition and the frequency-dependent selection that it generates, as well as stabilizing selection on a quantitative character. Our model includes local population growth, movements over space, natural selection, and gene flow. It simultaneously addresses the evolution of character displacement and species borders. It reproduces the Kirkpatrick and Barton single-species result that limited ranges can be produced with sufficiently steep environmental gradients and strong dispersal. Further, in the absence of environmental gradients or barriers to dispersal, interspecific competition will not limit species ranges at evolutionary equilibrium. However, interspecific competition can interact with environmental gradients and gene flow to generate limited ranges with much less extreme gradient and dispersal parameters than in the single-species case. Species display character displacement in sympatry, yet the reduction in competition that results from this displacement does not necessarily allow the two species to become sympatric everywhere. When species meet, competition reduces population densities in the region of overlap, which, in turn, intensifies the asymmetry in gene flow from center to margin. This reduces the ability of each species to adapt to local physical conditions at their range limits. If environmental gradients are monotonic but not linear, the transition zone between species at coevolutionary equilibrium occurs where the environmental gradient is steepest. If productivity gradients are also introduced into the model, then patterns similar to Rapoport's rule emerge. Interacting species respond to climate change, as it affects the optimal phenotype over space, by a combination of range shifts and local evolution in mean phenotype, while solitary species respond solely by range shifts. Finally, we compare empirical estimates for intrinsic growth rates and diffusion coefficients for several species to those needed by the single-species model to produce a stable limited range. These empirical values are generally insufficient to produce limited ranges in the model suggesting a role for interspecific interactions.     

Evolution of Complex Density-Dependent Dispersal Strategies

The question of how dispersal behavior is adaptive and how it responds to changes in selection pressure is more relevant than ever, as anthropogenic habitat alteration and climate change accelerate around the world. In metapopulation models where local populations are large, and thus local population size is measured in densities, density-dependent dispersal is expected to evolve to a single-threshold strategy, in which individuals stay in patches with local population density smaller than a threshold value and move immediately away from patches with local population density larger than the threshold. Fragmentation tends to convert continuous populations into metapopulations and also to decrease local population sizes. Therefore we analyze a metapopulation model, where each patch can support only a relatively small local population and thus experience demographic stochasticity. We investigated the evolution of density-dependent dispersal, emigration and immigration, in two scenarios: adult and natal dispersal. We show that density-dependent emigration can also evolve to a nonmonotone, “triple-threshold” strategy. This interesting phenomenon results from an interplay between the direct and indirect benefits of dispersal and the costs of dispersal. We also found that, compared to juveniles, dispersing adults may benefit more from density-dependent vs. density-independent dispersal strategies.     

Density-dependent dispersal complicates spatial synchrony in tri-trophic food chains

Spatial synchrony can increase extinction risk and undermines metapopulation persistence. Both dispersal and biotic interactions can strongly affect spatial synchrony. Here, we explore the spatial synchrony of a tri-trophic food chain in two patches connected by density-dependent dispersal, namely the strategies of prey evasion (PE) and predator pursuit (PP). The dynamics of the food chain are depicted by both the Hastings–Powell model and the chemostat model, with synchrony measured by the Pearson correlation coefficient. We use the density-independent dispersal in the system as a baseline for comparison. Results show that the density-independent dispersal of a species in the system can promote its dynamic synchrony. Dispersal of intermediate species in the tri-trophic food chain is the strongest synchronizer. In contrast, the density-dependent PP and PE of intermediate species can desynchronize the system. Highly synchronized dynamics emerged when the basal species has a strong PE strategy or when the top species has a moderate PP strategy. Our results reveal the complex relationship between density-dependent dispersal and spatial synchrony in tri-trophic systems.     

Density-dependent dispersal and spatial population dynamics

The synchronization of the dynamics of spatially subdivided populations is of both fundamental and applied interest in population biology. Based on theoretical studies, dispersal movements have been inferred to be one of the most general causes of population synchrony, yet no empirical study has mapped distance-dependent estimates of movement rates on the actual pattern of synchrony in species that are known to exhibit population synchrony. Northern vole and lemming species are particularly well-known for their spatially synchronized population dynamics. Here, we use results from an experimental study to demonstrate that tundra vole dispersal movements did not act to synchronize population dynamics in fragmented habitats. In contrast to the constant dispersal rate assumed in earlier theoretical studies, the tundra vole, and many other species, exhibit negative density-dependent dispersal. Simulations of a simple mathematical model, parametrized on the basis of our experimental data, verify the empirical results, namely that the observed negative density-dependent dispersal did not have a significant synchronizing effect.     

Negative density-dependent dispersal emerges from the joint evolution of density- and body condition-dependent dispersal strategies

Empirical studies have documented both positive and negative density-dependent dispersal, yet most theoretical models predict positive density dependence as a mechanism to avoid competition. Several hypotheses have been proposed to explain the occurrence of negative density-dependent dispersal, but few of these have been formally modeled. Here, we developed an individual-based model of the evolution of density-dependent dispersal. This model is novel in that it considers the effects of density on dispersal directly, and indirectly through effects on individual condition. Body condition is determined mechanistically, by having juveniles compete for resources in their natal patch. We found that the evolved dispersal strategy was a steep, increasing function of both density and condition. Interestingly, although populations evolved a positive density-dependent dispersal strategy, the simulated metapopulations exhibited negative density-dependent dispersal. This occurred because of the negative relationship between density and body condition: high density sites produced low-condition individuals that lacked the resources required for dispersal. Our model, therefore, generates the novel hypothesis that observed negative density-dependent dispersal can occur when high density limits the ability of organisms to disperse. We suggest that future studies consider how phenotype is linked to the environment when investigating the evolution of dispersal.     

Coexistence of competitors in metacommunities due to spatial variation in resource growth rates; does R * predict the outcome of competition?

Simple mathematical models are used to investigate the coexistence of two consumers using a single limiting resource that is distributed over distinct patches, and that has unequal growth rates in the different patches. Relatively low movement rates or high demographic rates of an inefficient resource exploiter allow it to coexist at a stable equilibrium with a more efficient species whose ratio of movement to demographic rates is lower. The range of conditions allowing coexistence depends on the between‐patch heterogeneity in resource growth rates, but this range can be quite broad. The between‐patch movement of the more efficient consumer turns patches with high resource growth rates into sources, while low‐growth‐rate patches effectively become sinks. A less efficient species can coexist with or even exclude the more efficient species from the global environment if it is better able to bias its spatial distribution towards the source patches. This can be accomplished with density independent dispersal if the less efficient species has a lower ratio of per capita between‐patch movement rate to demographic rates. Conditions that maximize the range of efficiencies allowing coexistence of two species are: a relatively high level of heterogeneity in resource growth conditions; high dispersal (or low demographic rates) of the superior competitor; and low dispersal (or high demographic rates) of the inferior competitor. Global exclusion of the more efficient competitor requires that the inferior competitor have sufficient movement to also produce a source‐sink environment.     

Density‐dependent dispersal and the formation of range borders.

Knowledge about the mechanisms of range formation is crucial for scientifically based species conservation strategies in the face of ongoing global climate change. In recent years an increasing amount of studies have focused on the influences of density‐dependent dispersal on demographic and biogeographical patterns. However, it still remains unclear, to what extent and in what ways this strategy would affect the range formation of species. In order to fill this gap, we present a study using individual‐based simulations of a species with discrete generations living along a dispersal mortality gradient. We compare the evolution of range sizes for species following density‐dependent and density‐independent emigration. Furthermore we assess the influence of environmental stochasticity and Allee effects on range formation, as both processes are known to play an important role for dispersal evolution. We find that density‐dependent dispersal always results in much wider ranges than unconditional dispersal. Increasing environmental stochasticity, a predicted consequence of climate change, can remarkably expand the ranges of species living in such connectivity gradients if dispersal decisions are based on local population density. A strong Allee effect causes range contraction for both strategies, but the effect is considerably less dramatic under density‐dependent compared to density‐independent emigration. We strongly recommend accounting for these findings in future attempts to model species’ range shifts due to climate change.     

Density-dependent dispersal in integrodifference equations

Many species exhibit dispersal processes with positive density- dependence. We model this behavior using an integrodifference equation where the individual dispersal probability is a monotone increasing function of local density. We investigate how this dispersal probability affects the spreading speed of a single population and its ability to persist in fragmented habitats. We demonstrate that density-dependent dispersal probability can act as a mechanism for coexistence of otherwise non-coexisting competitors. We show that in time-varying habitats, an intermediate dispersal probability will evolve. Analytically, we find that the spreading speed for the integrodifference equation with density-dependent dispersal probability is not linearly determined. Furthermore, the next-generation operator is not compact and, in general, neither order-preserving nor monotonicity-preserving. We give two explicit examples of non-monotone, discontinuous traveling-wave profiles.      

Density-dependent dispersal and relative dispersal affect the stability of predator-prey metacommunities

Although density-dependent dispersal and relative dispersal (the difference in dispersal rates between species) have been documented in natural systems, their effects on the stability of metacommunities are poorly understood. Here we investigate the effects of intra- and interspecific density-dependent dispersal on the regional stability in a predator-prey metacommunity model. We show that, when the dynamics of the populations reach equilibrium, the stability of the metacommunity is not affected by density-dependent dispersal. However, the regional stability, measured as the regional variability or the persistence, can be modified by density-dependent dispersal when local populations fluctuate over time. Moreover these effects depend on the relative dispersal of the predator and the prey. Regional stability is modified through changes in spatial synchrony. Interspecific density-dependent dispersal always desynchronizses local dynamics, whereas intraspecific density-dependent dispersal may either synchronize or desynchronize it depending on dispersal rates. Moreover, intra- and interspecific density-dependent dispersal strengthen the top-down control of the prey by the predator at intermediate dispersal rates. As a consequence the regional stability of the metacommunity is increased at intermediate dispersal rates. Our results show that density-dependent dispersal and relative dispersal of species are keys to understanding the response of ecosystems to fragmentation.     

Structure and dynamics of an amphibian metacommunity in two regions

1. The concept of metacommunity is based on the hypothesis that species occurrence depends on species dynamics and interactions on local and regional scales via the movements of individuals between localities. Metacommunity approaches are currently being applied to pond breeding taxa such as amphibians. 2. Given that animal movement is also influenced by the physical quality of the matrix to be crossed to reach a breeding habitat and by the affinity of the species for specific terrestrial habitats, matrix characteristics may enhance or hinder dispersal success. These characteristics would, in turn, affect the composition of larval assemblages at local level and, consequently, determine metacommunity structure and dynamics. 3. Here we compared the structures and dynamics of two metacommunities with the same pool of anurans along similar freshwater gradients in two regions that are well differentiated in terms of their respective terrestrial matrix. 4. Abundance of tadpole species and species assemblage in the two regions were determined principally by local processes (at pond level); however, the structure and dynamics of the communities differed. In one region species abundance was explained in part by landscape factors and consequently showed lower co-occurrence and lower colonization rates (species sorting models) indicating that terrestrial habitat could restrict animal movements, whereas in the other region higher co-occurrence and higher colonization rates (mass effect models) indicated low dispersal limitations.     

Dispersal among habitats varying in fitness: reciprocating migration through ideal habitat selection

Current evolutionary models of dispersal set the ends of a continuum where the number of individuals emigrating from a habitat either equals the number of individuals immigrating (balanced dispersal) or where emigrants flow from a source habitat to a corresponding sink. Theories of habitat selection suggest a more sophisticated conditional strategy where individuals disperse from habitats where they have the greatest impact on fitness to habitats where their per capita impact is lower. Asymmetries between periods of population growth and decline result in a reciprocating dispersal strategy where the direction of migration is reversed as populations wax and wane. Thus, for example, if net migration of individuals flows from high- to low-density habitats during periods of population growth, net migration will flow in the opposite direction during population decline. Stochastic simulations and analytical models of reciprocating dispersal demonstrate that fitness, carrying capacity, stochastic dynamics, and interference from dominants interact to determine whether dispersal is balanced between habitats, or whether one habitat or the other acts as a net donor of dispersing individuals. While the pattern of dispersal may vary, each is consistent with an underlying strategy of density-dependent habitat selection.     

On the survival of populations in a heterogeneous and variable environment

Summary The survival time of small and isolated populations will often be relatively low, by which the survival of species living in such a way will depend on powers of dispersal sufficiently high to result in a rate of population foundings that about compensates the rate of population extinctions. The survival time of composite populations uninterruptedly inhabiting large and heterogeneous areas, highly depends on the extent to which the numbers fluctuate unequally in the different subpopulations. The importance of this spreading of the risk of extinction over differently fluctuating subpopulations is demonstrated by comparing over 19 years the fluctuation patterns of the composite populations of two carabid species, Pterostichus versicolor with unequally fluctuating subpopulations, and Calathus melanocephalus with subpopulations fluctuating in parallel, both uninterruptedly occupying the same large heath area. The conclusions from the field data are checked by simulating the fluctuation patterns of these populations, and thus directly estimating survival times. It thus appeared that the former species can be expected to survive more than ten times better than the latter (other things staying the same). These simulations could also be used to study the possible influence of various density restricting processes in populations already fluctuating according to some pattern. As could be expected, the survival time of a population, which shows a tendency towards an upward trend in numbers, will be favoured by some kind of density restriction, but the degree to which these restrictions are density-dependent appeared to be immaterial. Density reductions that are about adequate on the average need even not occur at high densities only, if only the chance of occurrence at very low densities is low. The density-level at which a population is generally fluctuating appeared to be less important for survival than the fluctuation pattern itself, except for very low density levels, of course. The different ways in which deterministic and stochastic processes may interact and thus determine the fluctuations of population numbers are discussed. It is concluded that some stochastic processes will operate everywhere and will thus necessarily result in density fluctuations; such an omnipresence is much less imperative, however, for density-dependent processes, by which population models should primarily be stochastic models. However, if density-dependent processes are added to model populations, that are already fluctuating stochastically the effects are taken up into the general, stochastic fluctuation pattern, without altering it fundamentally.Communication No. 228 of the Biological Station WijsterDedicated to Professor Michael Evenari     

Self-diffusion of polymers in cartilage as studied by pulsed field gradient NMR

Pulsed field gradient (PFG) nuclear magnetic resonance (NMR) was used to investigate the self-diffusion behaviour of polymers in cartilage. Polyethylene glycol and dextran with different molecular weights and in different concentrations were used as model compounds to mimic the diffusion behaviour of metabolites of cartilage. The polymer self-diffusion depends extremely on the observation time: The short-time self-diffusion coefficients (diffusion time Delta approximately 15 ms) are subjected to a rather non-specific obstruction effect that depends mainly on the molecular weights of the applied polymers as well as on the water content of the cartilage. The observed self-diffusion coefficients decrease with increasing molecular weights of the polymers and with a decreasing water content of the cartilage. In contrast, the long-time self-diffusion coefficients of the polymers in cartilage (diffusion time Delta approximately 600 ms) reflect the structural properties of the tissue. Measurements at different water contents, different molecular weights of the polymers and varying observation times suggest that primarily the collagenous network of cartilage but also the entanglements of the polymer chains themselves are responsible for the observed restricted diffusion. Additionally, anomalous restricted diffusion was shown to occur already in concentrated polymer solutions.     

A multivariate analysis of fine-scale species density in the plant communities of a saltwater lagoon – the importance of disturbance intensity

Several factors might influence an organism's tendency or willingness to leave a given patch. One such factor is conspecific density, which may affect the per capita emigration rate. Some previous field studies on butterflies have reported positively density-dependent dispersal (emigration increases with population density) whereas the opposite, negatively density-dependent dispersal, has been found in other species. We investigated the effect of conspecific density on both the tendency to cross a patch boundary and within-patch mobility in Melitaea cinxia , by experimentally manipulating density in large outdoor cages divided into two habitat patches, separated by a barrier of unsuitable habitat. In contrast to previous results for M. cinxia , we found that the butterflies moved away from a patch at higher rates in high conspecific density (positively density-dependent emigration). The within-patch mobility, measured as the distance travelled per time unit, was however unaffected by butterfly density. A possible explanation for the seeming discrepancy with previous results could be that we used higher butterfly densities. For species with fluctuating population dynamics, such as M. cinxia , dispersal activity both at low and at high local density will be important for population phenomena such as fluctuations in distributional range over good and bad years.