Complex eco-evolutionary dynamics induced by the coevolution of predator–prey movement strategies

The coevolution of predators and prey has been the subject of much empirical and theoretical research that produced intriguing insights into the interplay of ecology and evolution. To allow for mathematical analysis, models of predator–prey coevolution are often coarse-grained, focussing on population-level processes and largely neglecting individual-level behaviour. As selection is acting on individual-level properties, we here present a more mechanistic approach: an individual-based simulation model for the coevolution of predators and prey on a fine-grained resource landscape, where features relevant for ecology (like changes in local densities) and evolution (like differences in survival and reproduction) emerge naturally from interactions between individuals. Our focus is on predator–prey movement behaviour, and we present a new method for implementing evolving movement strategies in an efficient and intuitively appealing manner. Throughout their lifetime, predators and prey make repeated movement decisions on the basis of their movement strategies. Over the generations, the movement strategies evolve, as individuals that successfully survive and reproduce leave their strategy to more descendants. We show that the movement strategies in our model evolve rapidly, thereby inducing characteristic spatial patterns like spiral waves and static spots. Transitions between these patterns occur frequently, induced by antagonistic coevolution rather than by external events. Regularly, evolution leads to the emergence and stable coexistence of qualitatively different movement strategies within the same population. Although the strategy space of our model is continuous, we often observe the evolution of discrete movement types. We argue that rapid evolution, coexistent movement types, and phase shifts between different ecological regimes are not a peculiarity of our model but a result of more realistic assumptions on eco-evolutionary feedbacks and the number of evolutionary degrees of freedom

Computational Study of the Structure and Thermodynamic Properties of Ammonium Chloride Clusters Using a Parallel J-Walking Approach

The thermodynamic and structural properties of (NH$_4$Cl)$_n$ clusters, n=3-10 are studied. Using the method of simulated annealing, the geometries of several isomers for each cluster size are examined. Jump-walking Monte Carlo simulations are then used to compute the constant-volume heat capacity for each cluster size over a wide temperature range. To carry out these simulations a new parallel algorithm is developed using the Parallel Virtual Machine (PVM) software package. Features of the cluster potential energy surfaces, such as energy differences among isomers and rotational barriers of the ammonium ions, are found to play important roles in determining the shape of the heat capacity curves.Comment: Journal of Chemical Physics, accepted for publicatio

Self-organization of collective escape in pigeon flocks

Bird flocks under predation demonstrate complex patterns of collective escape. These patterns may emerge by self-organization from local interactions among group-members. Computational models have been shown to be valuable for identifying what behavioral rules may govern such interactions among individuals during collective motion. However, our knowledge of such rules for collective escape is limited by the lack of quantitative data on bird flocks under predation in the field. In the present study, we analyze the first GPS trajectories of pigeons in airborne flocks attacked by a robotic falcon in order to build a species-specific model of collective escape. We use our model to examine a recently identified distance-dependent pattern of collective behavior: the closer the prey is to the predator, the higher the frequency with which flock members turn away from it. We first extract from the empirical data of pigeon flocks the characteristics of their shape and internal structure (bearing angle and distance to nearest neighbors). Combining these with information on their coordination from the literature, we build an agent-based model adjusted to pigeons’ collective escape. We show that the pattern of turning away from the predator with increased frequency when the predator is closer arises without prey prioritizing escape when the predator is near. Instead, it emerges through self-organization from a behavioral rule to avoid the predator independently of their distance to it. During this self-organization process, we show how flock members increase their consensus over which direction to escape and turn collectively as the predator gets closer. Our results suggest that coordination among flock members, combined with simple escape rules, reduces the cognitive costs of tracking the predator while flocking. Such escape rules that are independent of the distance to the predator can now be investigated in other species. Our study showcases the important role of computational models in the interpretation of empirical findings of collective behavior

Emergence of splits and collective turns in pigeon flocks under predation

Complex patterns of collective behaviour may emerge through self-organization, from local interactions among individuals in a group. To understand what behavioural rules underlie these patterns, computational models are often necessary. These rules have not yet been systematically studied for bird flocks under predation. Here, we study airborne flocks of homing pigeons attacked by a robotic falcon, combining empirical data with a species-specific computational model of collective escape. By analysing GPS trajectories of flocking individuals, we identify two new patterns of collective escape: early splits and collective turns, occurring even at large distances from the predator. To examine their formation, we extend an agent-based model of pigeons with a ‘discrete’ escape manoeuvre by a single initiator, namely a sudden turn interrupting the continuous coordinated motion of the group. Both splits and collective turns emerge from this rule. Their relative frequency depends on the angular velocity and position of the initiator in the flock: sharp turns by individuals at the periphery lead to more splits than collective turns. We confirm this association in the empirical data. Our study highlights the importance of discrete and uncoordinated manoeuvres in the collective escape of bird flocks and advocates the systematic study of their patterns across species

The evolution of social philopatry in female primates

The transition from solitary life to sociality is considered one of the major transitions in evolution. In primates, this transition is currently not well understood. Traditional verbal models appear insufficient to unravel the complex interplay of environmental and demographic factors involved in the evolution of primate sociality, and recent phylogenetic reconstructions have produced conflicting results. We therefore analyze a theoretical model for the evolution of female social philopatry that sheds new light on the question why most primates live in groups. In individual-based simulations, we study the evolution of dispersal strategies of both resident females and their offspring. The model reveals that social philopatry can evolve through kin selection, even if retention of offspring is costly in terms of within-group resource competition and provides no direct benefits. Our model supports the role of predator avoidance as a selective pressure for group-living in primates, but it also suggests that a second benefit of group-living, communal resource defense, might be required to trigger the evolution of sizable groups. Lastly, our model reveals that seemingly small differences in demographic parameters can have profound effects on primate social evolution

Some Causes of the Variable Shape of Flocks of Birds

Flocks of birds are highly variable in shape in all contexts (while travelling, avoiding predation, wheeling above the roost). Particularly amazing in this respect are the aerial displays of huge flocks of starlings (Sturnus vulgaris) above the sleeping site at dawn. The causes of this variability are hardly known, however. Here we hypothesise that variability of shape increases when there are larger local differences in movement behaviour in the flock. We investigate this hypothesis with the help of a model of the self-organisation of travelling groups, called StarDisplay, since such a model has also increased our understanding of what causes the oblong shape of schools of fish. The flocking patterns in the model prove to resemble those of real birds, in particular of starlings and rock doves. As to shape, we measure the relative proportions of the flock in several ways, which either depend on the direction of movement or do not. We confirm that flock shape is usually more variable when local differences in movement in the flock are larger. This happens when a) flock size is larger, b) interacting partners are fewer, c) the flock turnings are stronger, and d) individuals roll into the turn. In contrast to our expectations, when variability of speed in the flock is higher, flock shape and the positions of members in the flock are more static. We explain this and indicate the adaptive value of low variability of speed and spatial restriction of interaction and develop testable hypotheses

Fluctuation-Driven Flocking Movement in Three Dimensions and Scale-Free Correlation

Recent advances in the study of flocking behavior have permitted more sophisticated analyses than previously possible. The concepts of “topological distances” and “scale-free correlations” are important developments that have contributed to this improvement. These concepts require us to reconsider the notion of a neighborhood when applied to theoretical models. Previous work has assumed that individuals interact with neighbors within a certain radius (called the “metric distance”). However, other work has shown that, assuming topological interactions, starlings interact on average with the six or seven nearest neighbors within a flock. Accounting for this observation, we previously proposed a metric-topological interaction model in two dimensions. The goal of our model was to unite these two interaction components, the metric distance and the topological distance, into one rule. In our previous study, we demonstrated that the metric-topological interaction model could explain a real bird flocking phenomenon called scale-free correlation, which was first reported by Cavagna et al. In this study, we extended our model to three dimensions while also accounting for variations in speed. This three-dimensional metric-topological interaction model displayed scale-free correlation for velocity and orientation. Finally, we introduced an additional new feature of the model, namely, that a flock can store and release its fluctuations

Non-local kinetic and macroscopic models for self-organised animal aggregations

The last two decades have seen a surge in kinetic and macroscopic models derived to investigate the multi-scale aspects of self-organised biological aggregations. Because the individual-level details incorporated into the kinetic models (e.g., individual speeds and turning rates) make them somewhat difficult to investigate, one is interested in transforming these models into simpler macroscopic models, by using various scaling techniques that are imposed by the biological assumptions of the models. However, not many studies investigate how the dynamics of the initial models are preserved via these scalings. Here, we consider two scaling approaches (parabolic and grazing collision limits) that can be used to reduce a class of non-local 1D and 2D models for biological aggregations to simpler models existent in the literature. Then, we investigate how some of the spatio-temporal patterns exhibited by the original kinetic models are preserved via these scalings. To this end, we focus on the parabolic scaling for non-local 1D models and apply asymptotic preserving numerical methods, which allow us to analyse changes in the patterns as the scaling coefficient ϵ is varied from ϵ=1 (for 1D transport models) to ϵ=0 (for 1D parabolic models). We show that some patterns (describing stationary aggregations) are preserved in the limit ϵ→0, while other patterns (describing moving aggregations) are lost. To understand the loss of these patterns, we construct bifurcation diagrams