Predator and prey space use: dragonflies and tadpoles in an interactive game

Predator and prey spatial distributions have important population and community level consequences. However, little is known either theoretically or empirically about behavioral mechanisms that underlie the spatial patterns that emerge when predators and prey freely interact. We examined the joint space use and behavioral rules governing movement of freely interacting groups of odonate (dragonfly) predators and two size classes of anuran (tadpole) prey in arenas containing two patches with different levels of the prey's resource. Predator and prey movement and space use was quantified both when they were apart and together. When apart from predators, large tadpoles strongly preferred the high resource patch. When apart from prey, dragonflies weakly preferred the high resource patch. When together, large prey shifted to a uniform distribution, while predators strongly preferred the high resource patch. These patterns qualitatively fit the predictions of several three trophic level, ideal free distribution models. In contrast, the space use of small prey and predators did not deviate from uniform. Three measures of joint space use (spatial correlations, overlap, and co-occurrence) concurred in suggesting that prey avoidance of predators was more important than predator attraction to prey in determining overall spatial patterns. To gain additional insight into behavioral mechanisms, we used a model selection approach to identify behavioral movement rules that can potentially explain the observed, emergent patterns of space use. Prey were more likely to leave patches with more predators and more conspecific competitors; resources had relatively weak effects on prey movements. In contrast, predators were more likely to leave patches with low resources (that they do not consume) and more competing predators; prey had relatively little effect on predator movements. These results highlight the importance of investigating freely interacting predators and prey, the potential for simple game theory models to predict joint spatial distributions, and the utility of using model choice methods to identify potential key factors that govern movement.     

Indirect effects of prey swamping: differential seed predation during a bamboo masting event

Resource pulses often involve extraordinary increases in prey availability that "swamp" consumers and reverberate through indirect interactions affecting other community members. We developed a model that predicts predator-mediated indirect effects induced by an epidemic prey on co-occurring prey types differing in relative profitability/preference and validated our model by examining current-season and delayed effects of a bamboo mass seeding event on seed survival of canopy tree species in mixed Patagonian forests. The model shows that predator foraging behavior, prey profitability, and the scale of prey swamping influence the character and strength of short-term indirect effects on various alternative prey. When in large prey-swamped patches, nonselective predators decrease predation on all prey types. Selective predators, instead, only benefit prey of similar quality to the swamping species, while very low or high preference prey remain unaffected. Negative indirect effects (apparent competition) may override such positive effects (apparent mutualism), especially for highly preferred prey, when prey-swamped patches are small enough to allow predator aggregation and/or predators show a reproductive numerical response to elevated food supply. Seed predation patterns during bamboo (Chusquea culeou) masting were consistent with predicted short-term indirect effects mediated by a selective predator foraging in large prey-swamped patches. Bamboo seeds and similarly-sized Austrocedrus chilensis (ciprés) and Nothofagus obliqua (roble) seeds suffered lower predation in bamboo flowered than nonflowered patches. Predation rates on the small-seeded Nothofagus dombeyi (coihue) and the large-seeded Nothofagus alpina (rauli) were independent of bamboo flowering. Indirect positive effects were transient; three months after bamboo seeding, granivores preyed heavily upon all seed types, irrespective of patch flowering condition. Moreover, one year after bamboo seeding, predation rates on the most preferred seed (rauli) was higher in flowered than in nonflowered patches. Despite rapid predator numerical responses, short-term positive effects can still influence community recruitment dynamics because surviving seeds may find refuge beneath the litter produced by bamboo dieback. Together, our theoretical analysis and experiments indicate that indirect effects experienced by alternative prey during and after prey-swamping episodes need not be universal but can change across a prey quality spectrum, and they critically depend on predator-foraging rules and the spatial scale of swamping.     

Fine-scale habitat modeling of a top marine predator: do prey data improve predictive capacity?

Predators and prey assort themselves relative to each other, the availability of resources and refuges, and the temporal and spatial scale of their interaction. Predictive models of predator distributions often rely on these relationships by incorporating data on environmental variability and prey availability to determine predator habitat selection patterns. This approach to predictive modeling holds true in marine systems where observations of predators are logistically difficult, emphasizing the need for accurate models. In this paper, we ask whether including prey distribution data in fine-scale predictive models of bottlenose dolphin (Tursiops truncatus) habitat selection in Florida Bay, Florida, U.S.A., improves predictive capacity. Environmental characteristics are often used as predictor variables in habitat models of top marine predators with the assumption that they act as proxies of prey distribution. We examine the validity of this assumption by comparing the response of dolphin distribution and fish catch rates to the same environmental variables. Next, the predictive capacities of four models, with and without prey distribution data, are tested to determine whether dolphin habitat selection can be predicted without recourse to describing the distribution of their prey. The final analysis determines the accuracy of predictive maps of dolphin distribution produced by modeling areas of high fish catch based on significant environmental characteristics. We use spatial analysis and independent data sets to train and test the models. Our results indicate that, due to high habitat heterogeneity and the spatial variability of prey patches, fine-scale models of dolphin habitat selection in coastal habitats will be more successful if environmental variables are used as predictor variables of predator distributions rather than relying on prey data as explanatory variables. However, predictive modeling of prey distribution as the response variable based on environmental variability did produce high predictive performance of dolphin habitat selection, particularly foraging habitat.     

Linking edge effects and patch size effects: Importance of matrix nest predators

The edge effect is usually considered to be the proximate cause of area sensitivity in forest birds. We tested if birds nesting in large patches are less vulnerable to the edge effect using a simple model that assumes an increase in patch size reduces the probability of a matrix predator moving to the core areas of forest and that larger perimeter/area ratios result in a higher number of matrix predators per unit of area. The probability of a nest being successful decreased asymptotically with an increase in either the patch penetration distance of predators or predator density, but those effects were reduced when patch size was increased. Large patches have a lower probability of being affected by an Allee effect and they can function as sink habitats only if penetration distance and predator density are largely increased. However, the transition from an Allee effect to a sink condition occurs with a small increase in penetration distance and predator density. Since birds nesting in large patches are less vulnerable to an increase in matrix predator populations, persistence of bird populations may be possible by increasing the size of habitat patches that can act as source populations.     

How population dynamics shape the functional response in a one-predator-two-prey system

The type III functional response has historically been associated with switching predators; when there is a choice of prey the predator favors the more abundant prey type. Although this functional response has been found in experiments where both prey densities are manipulated, in real world studies the type II functional response is more commonly found. In modeling, the type III functional response is often used in systems where the second prey type is, implicitly, assumed to be constant. Here we define a functional response that takes into account both prey densities. This causes the functional response to show both type II and type III behavior, dependent on the interaction between the two prey densities. If we take into account population dynamics, we find a type II functional response in most cases, because predation regulates the relative prey densities. This explains why type III functional responses are found in experiments where both prey densities are manipulated, but type II functional responses occur when the feedback of population dynamics on the functional response is important. Furthermore, the results show that switching can have a stabilizing or destabilizing effect and can even lead to predator extinction.     

The influence of intraguild predation on prey suppression and prey release: a meta-analysis

Intraguild predation (IGP) occurs when one predator species consumes another predator species with whom it also competes for shared prey. One question of interest to ecologists is whether multiple predator species suppress prey populations more than a single predator species, and whether this result varies with the presence of IGP. We conducted a meta-analysis to examine this question, and others, regarding the effects of IGP on prey suppression. When predators can potentially consume one another (mutual IGP), prey suppression is greater in the presence of one predator species than in the presence of multiple predator species; however, this result was not found for assemblages with unidirectional or no IGP. With unidirectional IGP, intermediate predators were generally more effective than the top predator at suppressing the shared prey, in agreement with IGP theory. Adding a top predator to an assemblage generally caused prey to be released from predation, while adding an intermediate predator caused prey populations to be suppressed. However, the effects of adding a top or intermediate predator depended on the effectiveness of these predators when they were alone. Effects of IGP varied across different ecosystems (e.g., lentic, lotic, marine, terrestrial invertebrate, and terrestrial vertebrate), with the strongest patterns being driven by terrestrial invertebrates. Finally, although IGP theory is based on equilibrium conditions, data from short-term experiments can inform us about systems that are dominated by transient dynamics. Moreover, short-term experiments may be connected in some way to equilibrium models if the predator and prey densities used in experiments approximate the equilibrium densities in nature.     

A predator-prey model with satiation and intraspecific competition

We present a new predator-prey model where, except for the prey growth, assumed to be logistic, we endeavor to give some behavioral justification to all elements of the predator-prey interaction. The functional response takes account of predator satiation and predator competition. It is supported by some experimental evidence. We distinguish two contributions to the numerical response: the positive part, proportional to the functional response, is the birth rate of predators; the negative part is the death rate due to hunger.Two outcomes are possible. If the prey are unable to grow fast enough to replace the amount killed by the predators, both species become extinct. In the opposite case, both populations stabilize at a constant population. At this equilibrium level, the prey are not abundant enough to satiate the predators.The predation rate that allows the highest predator population is one half of the ideal prey growth rate. A higher exploitation rate can allow higher populations only temporarily. Evolved predator behavior, reguges for the prey, or other mechanisms can explain this regulation.Two more population behaviors (cycles and predator extinction) can be obtained with a time-lag in one of the responses. This is shown in a separate paper.The model is structurally stable. It can thus withstand small environmental perturbations.     

Critical patch sizes for food-web modules

Because patch size and connectivity may strongly impact the assemblage of species that occur on a patch, the types of food-web interactions that occur among those species may also depend on spatial structure. Here, we identify whether food-web interactions among salt-marsh-inhabiting arthropods vary with patch size and connectivity, and how such changes in trophic structure might feed back to influence the spatial distribution of prey. In a multiyear survey, patch-restricted predators exhibited steeper occupancy-patch-size relationships than herbivores, and species' critical patch sizes were correlated with overall rarity. As a result, the presence of food-web modules depended strongly on patch size: large and well-connected patches supported complex food-web modules, but only the simplest modules involving the most abundant species were found on small patches. Habitat-generalist spiders dominated on small patches, and predation pressure from such species may contribute to the observed lower densities of mesopredators on small patches. Overall, patch size and connectivity influenced the types of modules present on a patch through differential loss of rare, patch-restricted predators, but predation by generalist predators may be a key mechanism influencing the spatial structure of certain prey species.     

Behavioral response races, predator-prey shell games, ecology of fear, and patch use of pumas and their ungulate prey

The predator-prey shell game predicts random movement of prey across the landscape, whereas the behavioral response race and landscape of fear models predict that there should be a negative relationship between the spatial distribution of a predator and its behaviorally active prey. Additionally, prey have imperfect information on the whereabouts of their predator, which the predator should incorporate in its patch use strategy. I used a one-predator-one-prey system, puma (Puma concolor)-mule deer (Odocoileus hemionus) to test the following predictions regarding predator-prey distribution and patch use by the predator. (1) Pumas will spend more time in high prey risk/low prey use habitat types, while deer will spend their time in low-risk habitats. Pumas should (2) select large forage patches more often, (3) remain in large patches longer, and (4) revisit individual large patches more often than individual smaller ones. I tested these predictions with an extensive telemetry data set collected over 16 years in a study area of patchy forested habitat. When active, pumas spent significantly less time in open areas of low intrinsic predation risk than did deer. Pumas used large patches more than expected, revisited individual large patches significantly more often than smaller ones, and stayed significantly longer in larger patches than in smaller ones. The results supported the prediction of a negative relationship in the spatial distribution of a predator and its prey and indicated that the predator is incorporating the prey's imperfect information about its presence. These results indicate a behavioral complexity on the landscape scale that can have far-reaching impacts on predator-prey interactions.     

Mean free-path length theory of predator–prey interactions: Application to juvenile salmon migration

Ecological theory traditionally describes predator–prey interactions in terms of a law of mass action in which the prey mortality rate depends on the density of predators and prey. This simplifying assumption makes population-based models more tractable but ignores potentially important behaviors that characterize predator–prey dynamics. Here, we expand traditional predator–prey models by incorporating directed and random movements of both predators and prey. The model is based on theory originally developed to predict collision rates of molecules. The temporal and spatial dimensions of predators–prey encounters are determined by defining movement rules and the predator's field of vision. These biologically meaningful parameters can accommodate a broad range of behaviors within an analytically tractable framework suitable for population-based models. We apply the model to prey (juvenile salmon) migrating through a field of predators (piscivores) and find that traditional predator–prey models were not adequate to describe observations. Model parameters estimated from the survival of juvenile chinook salmon migrating through the Snake River in the northwestern United States are similar to estimates derived from independent approaches and data. For this system, we conclude that survival depends more on travel distance than travel time or migration velocity.     

Harvesting creates ecological traps: consequences of invisible mortality risks in predator-prey metacommunities

Models of two-patch predator-prey metacommunities are used to explore how the global predator population changes in response to additional mortality in one of the patches. This could describe the dynamics of a predator in an environment that includes a refuge area where that predator is protected and a spatially distinct ("risky") area where it is harvested. The predator's movement is based on its perceived fitness in the two patches, but the risk from the additional mortality is potentially undetectable; this often occurs when the mortality is from human harvesting or from a novel type of top predator. Increases in undetected mortality in the risky area can produce an abrupt collapse of either the refuge population or of the entire predator population when the mortality rate exceeds a threshold level. This is due to the attraction of the risky patch, which has abundant prey due to its high predator mortality. Extinction of the refuge predator population does not occur when the refuge patch has a higher maximum per capita predator growth rate than the exploited patch because the refuge is then more attractive when the predator is rare. The possibility of abrupt extinction of one or both patches from high densities in response to a small increase in harvest is often associated with alternative states. In such cases, large reductions in mortality may be needed to avoid extinction in a collapsing predator population, or to reestablish an extinct population. Our analysis provides a potential explanation for sudden collapses of harvested populations, and it argues for more consideration of adaptive movement in designing protected areas.     

A review of predator diet effects on prey defensive responses

Numerous studies have examined how predator diets influence prey responses to predation risk, but the role predator diet plays in modulating prey responses remains equivocal. We reviewed 405 predator–prey studies in 109 published articles that investigated changes in prey responses when predators consumed different prey items. In 54 % of reviewed studies, prey responses were influenced by predator diet. The value of responding based on a predator’s recent diet increased when predators specialized more strongly on particular prey species, which may create patterns in diet cue use among prey depending upon whether they are preyed upon by generalist or specialist predators. Further, prey can alleviate costs or accrue greater benefits using diet cues as secondary sources of information to fine tune responses to predators and to learn novel risk cues from exotic predators or alarm cues from sympatric prey species. However, the ability to draw broad conclusions regarding use of predator diet cues by prey was limited by a lack of research identifying molecular structures of the chemicals that mediate these interactions. Conclusions are also limited by a narrow research focus. Seventy percent of reviewed studies were performed in freshwater systems, with a limited range of model predator–prey systems, and 98 % of reviewed studies were performed in laboratory settings. Besides identifying the molecules prey use to detect predators, future studies should strive to manipulate different aspects of prey responses to predator diet across a broader range of predator–prey species, particularly in marine and terrestrial systems, and to expand studies into the field.     

Optimal harvesting of ecologically interdependent fish species

The optimal exploitation of a two-species predator-prey system is considered, using Lotka-Volterra-type equations. Due to the density-dependence of ecological efficiency, both species should be harvested simultaneously over a range of relative prices. Beyond the limits of this price range, either the prey species should be utilized indirectly by harvesting the predator, or the predator should be eliminated to maximize the prey yield. Neglecting harvesting costs, the simultaneous harvest of prey and predators requires that a unit of prey biomass increase in value by being “processed” by predators. Certain results from single-species fishery models are shown not to apply to multispecies models. These are as follows: (i) Optimal regulation of a free access fishery may call for subsidizing instead of taxing the harvest of predator species. (ii) Increasing the discount rate may, at “moderate” levels, imply that the optimal standing stock of biomass increases instead of decreasing. (iii) A rising price or a falling cost per unit fishing effort of a species may raise and not lower the optimal standing stock of that species.     

Patch use and vigilance by sympatric lemmings in predator and competitor-driven landscapes of fear

Prey living in risky environments can adopt a variety of behavioral tactics to reduce predation risk. In systems where predators regulate prey abundance, it is reasonable to assume that differential patterns of habitat use by prey species represent adaptive responses to spatial variation in predation. However, patterns of habitat use also reflect interspecific competition over habitat. Collared (Dicrostonyx groenlandicus) and brown (Lemmus trimucronatus) lemmings represent such a system and possess distinct upland tundra versus mesic meadow habitat preferences consistent with interspecific competition. Yet, we do not know whether this habitat preference might also reflect differences in predation risk or whether the two species differ in their behavioral tactics used to avoid predation. We performed experiments where we manipulated putative predation risk perceived by lemmings by increasing protective cover in upland and meadow habitats while we recorded lemming activity and behavior. Both lemming species preferentially used cover more than open patches, but Dicrostonyx was more vigilant than Lemmus. Both species also constrained their activity to protective patches in upland and meadow habitats, but during different periods of the day. Use of cover and vigilance were independent of habitat, suggesting that both species live in a fearsome but flattened landscape of fear at Walker Bay (Nunavut, Canada), and that their habitat preference is a consequence of competition rather than predation risk. Future studies aiming to map the contours of fear in multi-prey–predator systems should consider how predation and competition interact to modify prey species’ habitat preference, patch use, and vigilance.     

Minimizing predation risk in a landscape of multiple predators: effects on the spatial distribution of African ungulates

Studies that focus on single predator-prey interactions can be inadequate for understanding antipredator responses in multi-predator systems. Yet there is still a general lack of information about the strategies of prey to minimize predation risk from multiple predators at the landscape level. Here we examined the distribution of seven African ungulate species in the fenced Karongwe Game Reserve (KGR), South Africa, as a function of predation risk from all large carnivore species (lion, leopard, cheetah, African wild dog, and spotted hyena). Using observed kill data, we generated ungulate-specific predictions of relative predation risk and of riskiness of habitats. To determine how ungulates minimize predation risk at the landscape level, we explicitly tested five hypotheses consisting of strategies that reduce the probability of encountering predators, and the probability of being killed. All ungulate species avoided risky habitats, and most selected safer habitats, thus reducing their probability of being killed. To reduce the probability of encountering predators, most of the smaller prey species (impala, warthog, waterbuck, kudu) avoided the space use of all predators, while the larger species (wildebeest, zebra, giraffe) only avoided areas where lion and leopard space use were high. The strength of avoidance for the space use of predators generally did not correspond to the relative predation threat from those predators. Instead, ungulates used a simpler behavioral rule of avoiding the activity areas of sit-and-pursue predators (lion and leopard), but not those of cursorial predators (cheetah and African wild dog). In general, selection and avoidance of habitats was stronger than avoidance of the predator activity areas. We expect similar decision rules to drive the distribution pattern of ungulates in other African savannas and in other multi-predator systems, especially where predators differ in their hunting modes.     

The influence of size-specific indirect interactions in predator-prey systems

Nonlethal indirect interactions between predators often lead to nonadditive effects of predator number on prey survival and growth. Previous studies have focused on systems with at least two different predator species and one prey species. However, most predators undergo extreme ontological changes in phenotype such that interactions between different-sized cohorts of a predator and its prey could lead to nonadditive effects in systems with only two species. This may be important since different-sized individuals of the same species can differ more in their ecology than similar-sized individuals of different species. This study examined trait-mediated indirect effects in a two-species system including a cannibalistic predator with different-sized cohorts and its prey. I tested for these effects using larvae of two stream salamanders, Gyrinophilus porphyriticus (predator) and Eurycea cirrigera (prey), by altering the densities and combinations of predator size classes in experimental streams. Results showed that the presence of large individuals can significantly reduce the impact of density changes of smaller conspecifics on prey survival through nonlethal means. In the absence of large conspecifics, an increase in the relative frequency of small predators significantly increased predation rates, thereby reducing prey survival. However, with large conspecifics present, increasing the density of small predators did not decrease prey survival, resulting in a 14.3% lower prey mortality than predicted from the independent effects of both predator size classes. Small predators changed their microhabitat use in the presence of larger conspecifics. Prey individuals reduced activity in response to large predators but did not respond to small predators. Both predators reduced prey growth. These results demonstrate that the impact of a predator can be significantly altered by two different types of trait-mediated indirect effects in two-species systems: between different-sized cohorts and between different cohorts and prey. This study demonstrates that predictions based on simple numerical changes that assume independent effects of different size classes or ignore size structure can be strongly misleading. We need to account for the size structure within predator populations in order to predict how changes in predator abundance will affect predator-prey dynamics.     

Using stable isotopes to reveal shifts in prey consumption by generalist predators.

The effectiveness of generalist predators in biological control may be diminished if increased availability of alternative prey causes individual predators to decrease their consumption of crop pests. Farming practices that enhance densities of microbidetritivores in the detrital food web can lead to increased densities of generalist predators that feed on pest species. The ability to predict the net biocontrol impact of increased predator densities depends upon knowing the extent to which individual predators may shift to detrital prey and feed less on crop pests when prey of the detritus-based food web are more abundant. We addressed this question by comparing ratios of stable isotopes of carbon (delta13C) and nitrogen (delta15N) in generalist ground predators and two types of prey (crop pests and microbidetritivores) in replicated 8 x 8 m cucurbit gardens subjected to one of two treatments: a detrital subsidy or no addition of detritus (control). Small sheet-web spiders (Linyphiidae) and small wolf spiders (Lycosidae) had delta13C values similar to those of Collembola in both the detrital and control treatments, indicating that small spiders belong primarily to the detrital food web. In control plots the larger generalist predators had delta13C values similar to those of the major insect pests, consistent with their known effectiveness as biocontrol agents. Adding detritus may have caused delta13C of one species of large wolf spider to shift toward that of the microbi-detritivores, although evidence is equivocal. In contrast, another large wolf spider displayed no shift in delta13C in the detrital treatment. Thus, stable isotopes revealed which generalist predators will likely continue to feed on pest species in the presence of greater densities of alternative prey.     

Cost-Effective Suppression and Eradication of Invasive Predators

Abstract: Introduced predators can have pronounced effects on naïve prey species; thus, predator control is often essential for conservation of threatened native species. Complete eradication of the predator, although desirable, may be elusive in budget‐limited situations, whereas predator suppression is more feasible and may still achieve conservation goals. We used a stochastic predator–prey model based on a Lotka‐Volterra system to investigate the cost‐effectiveness of predator control to achieve prey conservation. We compared five control strategies: immediate eradication, removal of a constant number of predators (fixed‐number control), removal of a constant proportion of predators (fixed‐rate control), removal of predators that exceed a predetermined threshold (upper‐trigger harvest), and removal of predators whenever their population falls below a lower predetermined threshold (lower‐trigger harvest). We looked at the performance of these strategies when managers could always remove the full number of predators targeted by each strategy, subject to budget availability. Under this assumption immediate eradication reduced the threat to the prey population the most. We then examined the effect of reduced management success in meeting removal targets, assuming removal is more difficult at low predator densities. In this case there was a pronounced reduction in performance of the immediate eradication, fixed‐number, and lower‐trigger strategies. Although immediate eradication still yielded the highest expected minimum prey population size, upper‐trigger harvest yielded the lowest probability of prey extinction and the greatest return on investment (as measured by improvement in expected minimum population size per amount spent). Upper‐trigger harvest was relatively successful because it operated when predator density was highest, which is when predator removal targets can be more easily met and the effect of predators on the prey is most damaging. This suggests that controlling predators only when they are most abundant is the “best” strategy when financial resources are limited and eradication is unlikely.     

Effects of density-dependent dispersal behaviours on the speed and spatial patterns of range expansion in predator-prey metapopulations

Dispersal can strongly affect the spatiotemporal dynamics of a species (its spread, spatial distribution and persistence). We investigated how two dispersal behaviours, namely prey evasion (PE) and predator pursuit (PP), affect the dynamics of a predator-prey system. PE portrays the tendency of prey avoiding predators by dispersing into adjacent patches with fewer predators, while PP describes the tendency of predators to pursue the prey by moving into patches with more prey. Based on the Beddington predation model, a spatially explicit metapopulation model was built to incorporate PE and PP. Numerical simulations were run to investigate the effects of PE and PP on the rate of spread, spatial synchrony and the persistence of populations. Results show that both PE and PP can alter spatial synchrony although PP has a weaker desynchronising effect than PE. The predator-prey system without PE and PP expanded in circular waves. The effect of PE can push the prey to distribute in a circular ring front, whereas the effect of PP can change the circular waves to anisotropic expansion. Furthermore, weak PE and PP can accelerate the spread of prey while strong and disproportionate intensities slow down the range expansion. The effects of PE and PP further enhance the population size, break down the spatial synchrony and promote the persistence of populations.     

The effects of opportunistic and intentional predators on the herding behavior of prey

In this article, we study how predator behavior influences the aggregation of prey into herds. Game-theoretic models of herd formation are developed based on different survival probabilities of solitary prey and prey that join the herd and on the predator's preference of what type of prey to search for. For an intentional predator that will only pursue its preferred type of prey, a single herd with no solitaries cannot emerge unless the herd acts as a prey refuge. If neither prey choice provides a refuge, it is shown that an equilibrium always exists where there are both types of prey and the predator does not always search for the same type of prey (i.e., a mixed equilibrium exists). On the other hand, if the predator is opportunistic in that it sometimes shifts to pursue the type of prey that is observed first, there may be a single herd equilibrium that does not act as a prey refuge when there is a high level of opportunistic behavior. For low opportunistic levels, a mixed equilibrium is again the only outcome. The evolutionary stability of each equilibrium is tested to see if it predicts the eventual herding behavior of prey in its corresponding model. Our analysis confirms that both predator and prey preferences (for herd or solitary) have strong effects on why prey aggregate. In particular, in our models, only the opportunistic predator can maintain all prey in a single herd that is under predation risk.