Ideal free distributions, evolutionary games, and population dynamics in multiple-species environments

In this article, we develop population game theory, a theory that combines the dynamics of animal behavior with population dynamics. In particular, we study interaction and distribution of two species in a two-patch environment assuming that individuals behave adaptively (i.e., they maximize Darwinian fitness). Either the two species are competing for resources or they are in a predator-prey relationship. Using some recent advances in evolutionary game theory, we extend the classical ideal free distribution (IFD) concept for single species to two interacting species. We study population dynamical consequences of two-species IFD by comparing two systems: one where individuals cannot migrate between habitats and one where migration is possible. For single species, predator-prey interactions, and competing species, we show that these two types of behavior lead to the same population equilibria and corresponding species spatial distributions, provided interspecific competition is patch independent. However, if differences between patches are such that competition is patch dependent, then our predictions strongly depend on whether animals can migrate or not. In particular, we show that when species are settled at their equilibrium population densities in both habitats in the environment where migration between habitats is blocked, then the corresponding species spatial distribution need not be an IFD. Thus, when species are given the opportunity to migrate, they will redistribute to reach an IFD (e.g., under which the two species can completely segregate), and this redistribution will also influence species population equilibrial densities. Alternatively, we also show that when two species are distributed according to the IFD, the corresponding population equilibrium can be unstable.     

Putting competition strategies into ideal free distribution models: habitat selection as a tug of war

When resources are patchily distributed in an environment, behavioral ecologists frequently turn to ideal free distribution (IFD) models to predict the spatial distribution of organisms. In these models, predictions about distributions depend upon two key factors: the quality of habitat patches and the nature of competition between consumers. Surprisingly, however, no IFD models have explored the possibility that consumers modulate their competitive efforts in an evolutionarily stable manner. Instead, previous models assume that resource acquisition ability and competition are fixed within species or within phenotypes. We explored the consequences of adaptive modulation of competitive effort by incorporating tug-of-war theory into payoff equations from the two main classes of IFD models (continuous input (CI) and interference). In the models we develop, individuals can increase their share of the resources available in a patch, but do so at the costs of increased resource expenditures and increased negative interactions with conspecifics. We show how such models can provide new hypotheses to explain what are thought to be deviations from IFDs (e.g., the frequent observation of fewer animals than predicted in "good" patches of habitat). We also detail straightforward predictions made uniquely by the models we develop, and we outline experimental tests that will distinguish among alternatives.     

Foraging on spatially distributed resources with sub‐optimal movement,imperfect information,and travelling costs: departures from the ideal free distribution

Ideal free distribution (IFD) theory offers an important baseline for predicting the distribution of foragers across resource patches. Yet it is well known that IFD theory relies on several over‐simplifying assumptions that are unlikely to be met in reality. Here we relax three of the most critical assumptions: (1) optimal foraging moves among patches, (2) omniscience about the utility of resource patches, and (3) cost‐free travelling between patches. Based on these generalizations, we investigate the distributions of a constant number of foragers in models with explicit resource dynamics of logistic type. We find that, first, when foragers do not always move to the patch offering maximum intake rate (optimal foraging), but instead move probabilistically according to differences in resource intake rates between patches (sub‐optimal foraging), the distribution of foragers becomes less skewed than the IFD, so that high‐quality patches attract fewer foragers. Second, this homogenization is strengthened when foragers have less than perfect knowledge about the utility of resource patches. Third, and perhaps most surprisingly, the introduction of travelling costs causes departures in the opposite direction: the distribution of sub‐optimal foragers approaches the IFD as travelling costs increase. We demonstrate that these three findings are robust when considering patches that differ in the resource's carrying capacity or intrinsic growth rate, and when considering simple two‐patch and more complex multiple‐patch models. By overcoming three major over‐simplifications of IFD theory, our analyses contribute to the systematic investigation of ecological factors influencing the spatial distribution of foragers, and thus help in deriving new hypotheses that are testable in empirical systems. A confluence of theoretical and empirical studies that go beyond classical IFD theory is essential for improving insights into how animal distributions across resource patches are determined in nature.     

The ideal free distribution: a review and synthesis of the game-theoretic perspective

The Ideal Free Distribution (IFD), introduced by Fretwell and Lucas in [Fretwell, D.S., Lucas, H.L., 1970. On territorial behavior and other factors influencing habitat distribution in birds. Acta Biotheoretica 19, 16-32] to predict how a single species will distribute itself among several patches, is often cited as an example of an evolutionarily stable strategy (ESS). By defining the strategies and payoffs for habitat selection, this article puts the IFD concept in a more general game-theoretic setting of the “habitat selection game”. Within this game-theoretic framework, the article focuses on recent progress in the following directions: (1) studying evolutionarily stable dispersal rates and corresponding dispersal dynamics; (2) extending the concept when population numbers are not fixed but undergo population dynamics; (3) generalizing the IFD to multiple species.For a single species, the article briefly reviews existing results. It also develops a new perspective for Parker’s matching principle, showing that this can be viewed as the IFD of the habitat selection game that models consumer behavior in several resource patches and analyzing complications involved when the model includes resource dynamics as well. For two species, the article first demonstrates that the connection between IFD and ESS is now more delicate by pointing out pitfalls that arise when applying several existing game-theoretic approaches to these habitat selection games. However, by providing a new detailed analysis of dispersal dynamics for predator-prey or competitive interactions in two habitats, it also pinpoints one approach that shows much promise in this general setting, the so-called “two-species ESS”. The consequences of this concept are shown to be related to recent studies of population dynamics combined with individual dispersal and are explored for more species or more patches.     

State-dependent ideal free distributions

Summary The standard ideal free distribution (IFD) states how animals should distribute themselves at a stable competitive equilibrium. The equilibrium is stable because no animal can increase its fitness by changing its location. In applying the IFD to choice between patches of food, fitness has been identified with the net rate of energetic gain. In this paper we assess fitness in terms of survival during a non-reproductive period, where the animal may die as a result of starvation or predation. We find the IFD when there is a large population that can distribute itself between two patches of food. The IFD in this case is state-dependent, so that an animal's choice of patch depends on its energy reserves. Animals switch between patches as their reserves change and so the resulting IFD is a dynamic equilibrium. We look at two cases. In one there is no predation and the patches differ in their variability. In the other, patches differ in their predation risk. In contrast to previous IFDs, it is not necessarily true that anything is equalized over the two patches.     

Dynamical facilitation of the ideal free distribution in nonideal populations

The ideal free distribution (IFD) requires that individuals can accurately perceive density‐dependent habitat quality, while failure to discern quality differences below a given perception threshold results in distributions approaching spatial uniformity. Here, we investigate the role of population growth in restoring a nonideal population to the IFD. We place a simple model of discrete patch choice under limits to the resolution by which patch quality is perceived and include population growth driven by that underlying quality. Our model follows the population's distribution through both breeding and dispersal seasons when perception limits differ in their likely influence. We demonstrate that populations of perception limited movers can approximate an IFD provided sufficient population growth; however, the emergent IFD would be temporally inconstant and correspond to reproductive events. The time to emergence of the IFD during breeding is shorter under exponential growth than under logistic growth. The IFD during early colonization of a community persists longer when more patches are available to individuals. As the population matures and dispersal becomes increasingly random, there is an oscillation in the observance of IFD, with peaks most closely approximating the IFD occurring immediately after reproductive events, and higher reproductive rates producing distributions closer to the IFD.     

Integrating individual behavior and population ecology: the potential for habitat-dependent population regulation in a reef fish

We used the predictions of the ideal free and ideal despoticdistributions (IFD and IDD, respectively) as a basis to evaluatethe link between spatial heterogeneity, behavior, and populationdynamics in a Caribbean coral reef fish. Juvenile three-spotdamselfish (Stegastes planifrons) were more closely aggregatedin patch reef habitat than on continuous back reef. Agonisticinteractions were more frequent but feeding rates were lowerin the patch versus the continuous reef habitat. Growth rateswere lower in patch reef habitat than on the continuous reef,but mortality rates did not differ. A separate experiment usingstandard habitat units demonstrated that the patterns observedin natural habitat were the result of the spatial distributionof the habitat patches rather than resource differences between habitats. Our results do not follow the predictions of simpleIFD or IDD models. This deviation from IFD and IDD predictionsmay be the result of a number of factors, including lack ofperfect information about habitat patches, high movement costs,and higher encounter rates of dispersed patches. Our resultsdemonstrate that behavioral interactions are an integral partof population dynamics and that it is necessary to considerthe spatial organization of the habitat in both behavioraland ecological investigations.     

Physiological Ecology Meets the Ideal-free Distribution: Predicting the Distribution of Size-structured Fish Populations Across Temperature Gradients

We describe a habitat selection model that predicts the distribution of size-structured groups of fish in a habitat where food availability and water temperature vary spatially. This model is formed by combining a physiological model of fish growth with the logic of ideal free distribution (IFD) theory. In this model we assume that individuals scramble compete for resources, that relative competitive abilities of fish vary with body size, and that individuals select patches that maximize their growth rate. This model overcomes limitations in currently existing physiological and IFD-based models of habitat selection. This is because existing physiological models do not take into account the fact that the amount of food consumed by a fish in a patch will depend on the number of competitors there (something that IFD theory addresses), while traditional IFD models do not take into account the fact that fish are likely to choose patches based on potential growth rate rather than gross food intake (something that physiological models address). Our model takes advantage of the complementary strengths of these two approaches to overcome these weaknesses. Reassuringly, our model reproduces the predictions of its two constituent models under the simple conditions where they apply. When there is no competition for resources it mimics the physiological model of habitat selection, and when there is competition but no temperature variation between patches it mimics either the simple IFD model or the IFD model for unequal competitors. However, when there are both competition and temperature differences between patches our model makes different predictions. It predicts that input-matching between the resource renewal rate and the number of fish (or competitive units) in a patch, the hallmark of IFD models, will be the exception rather than the rule. It also makes the novel prediction that temperature based size-segregation will be common, and that the strength and direction of this segregation will depend on per capita resource renewal rates and the manner in which competitive weight scales with body size. Size-segregation should become more pronounced as per capita resource abundance falls. A larger fish/cooler water pattern is predicted when competitive ability increases more slowly than maximum ration with body size, and a smaller fish/cooler water pattern is predicted when competitive ability increases more rapidly than maximum ration with body size.     

The effect of travel costs on the ideal free distribution in stickleback

The ideal free distribution (IFD) theory, which predicts that a population of individuals will match the distribution of a patchily distributed resource, is widely used in ecology to describe the spatial distribution of animals. While many studies have shown general support of its habitat matching prediction, others have described a systematic pattern of undermatching, where too many animals feed at patches with fewer resources, and too few animals feed in richer patches. These results have been attributed to deviations from several of the assumptions of the IFD. One possible variable, the cost of travelling between patches, has received little attention. Here, we investigated the impact on resource matching when travel costs were manipulated in a simple laboratory experiment involving two continuous input patches. This experiment allowed us to control for extraneous variables and decouple time costs from energetic costs of travel. Two experiments examined the impact of varying travel costs on movement rates between foraging patches and how these travel costs impact conformity to the IFD. Our data demonstrated that there was less movement between patches and greater discrepancies from the IFD predictions as the cost of travel increased.     

The Allee-type ideal free distribution

The ideal free distribution (IFD) in a two-patch environment where individual fitness is positively density dependent at low population densities is studied. The IFD is defined as an evolutionarily stable strategy of the habitat selection game. It is shown that for low and high population densities only one IFD exists, but for intermediate population densities there are up to three IFDs. Population and distributional dynamics described by the replicator dynamics are studied. It is shown that distributional stability (i.e., IFD) does not imply local stability of a population equilibrium. Thus distributional stability is not sufficient for population stability. Results of this article demonstrate that the Allee effect can strongly influence not only population dynamics, but also population distribution in space.     

The effect of travel loss on evolutionarily stable distributions of populations in space

A key assumption of the ideal free distribution (IFD) is that there are no costs in moving between habitat patches. However, because many populations exhibit more or less continuous population movement between patches and traveling cost is a frequent factor, it is important to determine the effects of costs on expected population movement patterns and spatial distributions. We consider a food chain (tritrophic or bitrophic) in which one species moves between patches, with energy cost or mortality risk in movement. In the two-patch case, assuming forced movement in one direction, an evolutionarily stable strategy requires bidirectional movement, even if costs during movement are high. In the N-patch case, assuming that at least one patch is linked bidirectionally to all other patches, optimal movement rates can lead to source-sink dynamics where patches with negative growth rates are maintained by other patches with positive growth rates. As well, dispersal between patches is not balanced (even in the two-patch case), leading to a deviation from the IFD. Our results indicate that cost-associated forced movement can have important consequences for spatial metapopulation dynamics. Relevance to marine reserve design and the study of stream communities subject to drift is discussed.     

Whether ideal free or not,predatory mites distribute so as to maximize reproduction

Ideal free distribution (IFD) models predict that animals distribute themselves such that no individual can increase its fitness by moving to another patch. Many empirical tests assume that the interference among animals is independent of density and do not quantify the effects of density on fitness traits. Using two species of predatory mites, we measured oviposition as a function of conspecific density. Subsequently, we used these functions to calculate expected distributions on two connected patches. We performed an experimental test of the distributions of mites on two such connected patches, among which one had a food accessibility rate that was twice as high as on the other. For one of the two species, Iphiseius degenerans, the distribution matched the expected distribution. The distribution also coincided with the ratio of food accessibility. The other species, Neoseiulus cucumeris, distributed itself differently than expected. However, the oviposition rates of both species did not differ significantly from the expected oviposition rates based on experiments on single patches. This suggests that the oviposition rate of N. cucumeris was not negatively affected by the observed distribution, despite the fact that N. cucumeris did not match the predicted distributions. Thus, the distribution of one mite species, I. degenerans, was in agreement with IFD theory, whereas for the other mite species, N. cucumeris, unknown factors may have influenced the distribution of the mites. We conclude that density-dependent fitness traits provide essential information for explaining animal distributions.     

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

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

Maintaining social cohesion is a more important determinant of patch residence time than maximizing food intake rate in a group-living primate,Japanese macaque (Macaca fuscata)

Animals have been assumed to employ an optimal foraging strategy (e.g., rate-maximizing strategy). In patchy food environments, intake rate within patches is positively correlated with patch quality, and declines as patches are depleted through consumption. This causes patch-leaving and determines patch residence time. In group-foraging situations, patch residence times are also affected by patch sharing. Optimal patch models for groups predict that patch residence times decrease as the number of co-feeding animals increases because of accelerated patch depletion. However, group members often depart patches without patch depletion, and their patch residence time deviates from patch models. It has been pointed out that patch residence time is also influenced by maintaining social proximity with others among group-living animals. In this study, the effects of maintaining social cohesion and that of rate-maximizing strategy on patch residence time were examined in Japanese macaques (Macaca fuscata). I hypothesized that foragers give up patches to remain in the proximity of their troop members. On the other hand, foragers may stay for a relatively long period when they do not have to abandon patches to follow the troop. In this study, intake rate and foraging effort (i.e., movement) did not change during patch residency. Macaques maintained their intake rate with only a little foraging effort. Therefore, the patches were assumed to be undepleted during patch residency. Further, patch residence time was affected by patch-leaving to maintain social proximity, but not by the intake rate. Macaques tended to stay in patches for short periods when they needed to give up patches for social proximity, and remained for long periods when they did not need to leave to keep social proximity. Patch-leaving and patch residence time that prioritize the maintenance of social cohesion may be a behavioral pattern in group-living primates.     

Spatial and temporal variation in patch occupancy and population density in a model system of an arctic Collembola species assemblage

We employed an experimental model system to investigate the mechanisms underlying patterns of patch occupancy and population density in a high arctic assemblage of Collembola species inhabiting a sedge tussock landscape on Svalbard. The replicate model systems consisted of 5 cores of the tussocks (habitat patches) imbedded in a barren matrix. Four of the patches were open so that animals could migrate between them, while there was one closed patch per system to test the effect of migration on extinction rate. There were model systems of two types: one with long and one with short inter‐patch distances to test the effect of patch isolation on colonisation and extinction rates. Each of the four most common collembolan species at the field site were introduced to two open patches per system (source patches), with the other two functioning as colonisation patches for the species. The experiment was run in an ecotrone over three identical, simulated arctic summers separated by winters of 3 weeks. Six replicates of systems with short and long inter‐patch distances were sampled at the end of each summer. The species varied markedly in their performance in both open arenas and closed patches, indicating differential responses to patch humidity, consistent with their differential distribution along the moisture gradient in the field site. The extinction – colonisation dynamics differed markedly between species as predicted from our field studies. This could partly be ascribed to differential dispersal and colonisation ability, but also to different tolerance to spatially variable patch quality and/or tendency for aggregative behaviour. Three of the species exhibited dynamics that superficially resemble what could be expected from classical metapopulation dynamics. However, there was a striking discrepancy between what would be expected from the effect of migration on the extinction rate of isolated patches (in particular closed patches) and the observed rates. Thus, metapopulation processes, such as stochastic colonisation and extinction events due to demographic stochasticity, were relatively unimportant compared to other sources of spatial variability among which subtle differences in patch quality are probably most important. We discuss the value of combining field studies with model system experiments, in particular when habitat quality cannot easily be measured in the field. However, our field and laboratory studies also emphasise the need for a thorough knowledge of species‐specific life history traits for making biologically sound interpretations based on both observational and experimental data.     

Behavioral choices based on patch selection: a model using aggregation methods

The aim of this work is to study the influence of patch selection on the dynamics of a system describing the interactions between two populations, generically called 'population N' and 'population P'. Our model may be applied to prey-predator systems as well as to certain host-parasite or parasitoid systems. A situation in which population P affects the spatial distribution of population N is considered. We deal with a heterogeneous environment composed of two spatial patches: population P lives only in patch 1, while individuals belonging to population N migrate between patch 1 and patch 2, which may be a refuge. Therefore they are divided into two patch sub-populations and can migrate according to different migration laws. We make the assumption that the patch change is fast, whereas the growth and interaction processes are slower. We take advantage of the two time scales to perform aggregation methods in order to obtain a global model describing the time evolution of the total populations, at a slow time scale. At first, a migration law which is independent on population P density is considered. In this case the global model is equivalent to the local one, and under certain conditions, population P always gets extinct. Then, the same model, but in which individuals belonging to population N leave patch 1 proportionally to population P density, is studied. This particular behavioral choice leads to a dynamically richer global system, which favors stability and population coexistence. Finally, we study a third example corresponding to the addition of an aggregative behavior of population N on patch 1. This leads to a more complicated situation in which, according to initial conditions, the global system is described by two different aggregated models. Under certain conditions on parameters a stable limit cycle occurs, leading to periodic variations of the total population densities, as well as of the local densities on the spatial patches.     

Habitat assessment by parasitoids: consequences for population distribution

The ideal free distribution (IFD) is a stable distribution ofcompetitors among resource patches. For equally efficient competitors,equilibrium is reached when the per capita rate of intake equalizesacross patches. The seminal version of the IFD assumes omniscience,but populations may still converge toward the equilibrium providedthat competitors 1) accurately assess their environment by learningand 2) remain for an optimal (rate-maximizing) time on eachencountered patch. In the companion article (Tentelier C, DesouhantE, Fauvergue X. 2006. Habitat assessment by parasitoids: mechanismsfor patch time allocation. Behav Ecol. Forthcoming), it is shownthat the parasitoid wasp Lysiphlebus testaceipes adapts itsexploitation of aphid host colonies based on previous experience,in a manner consistent with these two conditions. We thereforepredicted that a randomly distributed population of initiallynaive wasps should converge toward the IFD. We tested this predictionby introducing 1300 L. testaceipes females into a 110-m² greenhousecontaining 40 host patches. Just after introduction, the parasitoidrate of gain was positively affected by host number and negativelyaffected by parasitoid number but, as predicted, these effectsvanished in the course of the experiment. Six hours after introduction,the expected rate of gain reached a constant. Surprisingly,this passage through equilibrium was not accompanied by a decreasein the coefficient of variation among gain rates or by a shiftfrom a random to an aggregated distribution of parasitoids.These results challenge our understanding of the link betweenindividual behavior and population distribution.     

The ideal free distribution of clonal plant's ramets among patches in a heterogeneous environment

The plastic response of clonal plant to different patch quality is not always the same and the degree is different too. So the result of this kind of foraging behaviour is different. In order to make clear whether the ramtes stay in favourable patches and get the quantitative relationship between the ramets distribution among patches and the available resource amount in heterogeneous environment, we develop a theoretical work under ideal free distribution (IFD) theory framework by neglecting some morphological plasticity of the spacer in this article. The results of our general model show that the ramet distribution should obey input matching rule at equilibrium. That means the ratio of ramet number in different patches should be equal to the ratio of available resource amount in these patches. We also use the simulation to predict the distribution pattern under history mattering. The results show that the initial ramets number has significant influence on the final distribution: over matching and under matching both can occur. More initial ramets in favourable patch result in over matching and more initial ramets in unfavourable patch result in under matching. The degree of the deviation from input matching rule is great when the difference of patches is small. These results prove that ideal free distribution theory works the same with animals. The ramets can stay in favourable patches sometimes in spite of the plasticity of the spacer, and the distribution depends on both patch quality and the history factors. But these results are true only when the functional response is type II.     

Distribution of a naturally fluctuating ungulate population among heterogeneous plant communities: ideal and free?

1. Herbivore distribution is often assumed to follow the ideal free distribution (IFD) model. This assumes that organisms are omniscient about forage quality and availability within the area available to them and are free to move, with negligible cost, throughout this environment. If this were the case we would expect that, at lowest densities, all animals would be found in the best habitat patches, with less desirable habitats being occupied stepwise as population density increases. We test this using data from a naturally fluctuating population of feral Soay sheep. 2. We show that, although the distribution of individuals is correlated positively with food quality, in line with patterns reported for hill sheep in Scotland, their distribution does not conform to the predictions of the IFD model. We argue that it is the dynamic nature of their food resource that causes this departure from the predictions of the IFD model and make the case that the IFD model, in its unmodified form, is inappropriate for use in modelling distribution among patches containing dynamic resources.