Fluctuations and correlations in an individual-based model of biological coevolution

Tangent map intermittency as an approximate analysis of intermittency in a high dimensional fully stochastic dynamical system: The Tangled Nature model

It is well known that low-dimensional nonlinear deterministic maps close to a tangent bifurcation exhibit intermittency and this circumstance has been exploited, e.g., by Procaccia and Schuster [Phys. Rev. A 28, 1210 (1983)], to develop a general theory of 1/f spectra. This suggests it is interesting to study the extent to which the behavior of a high-dimensional stochastic system can be described by such tangent maps. The Tangled Nature (TaNa) Model of evolutionary ecology is an ideal candidate for such a study, a significant model as it is capable of reproducing a broad range of the phenomenology of macroevolution and ecosystems. The TaNa model exhibits strong intermittency reminiscent of punctuated equilibrium and, like the fossil record of mass extinction, the intermittency in the model is found to be non-stationary, a feature typical of many complex systems. We derive a mean-field version for the evolution of the likelihood function controlling the reproduction of species and find a local map close to tangency. This mean-field map, by our own local approximation, is able to describe qualitatively only one episode of the intermittent dynamics of the full TaNa model. To complement this result, we construct a complete nonlinear dynamical system model consisting of successive tangent bifurcations that generates time evolution patterns resembling those of the full TaNa model in macroscopic scales. The switch from one tangent bifurcation to the next in the sequences produced in this model is stochastic in nature, based on criteria obtained from the local mean-field approximation, and capable of imitating the changing set of types of species and total population in the TaNa model. The model combines full deterministic dynamics with instantaneous parameter random jumps at stochastically drawn times. In spite of the limitations of our approach, which entails a drastic collapse of degrees of freedom, the description of a high-dimensional model system in terms of a low-dimensional one appears to be illuminating.

Random walk in genome space: a key ingredient of intermittent dynamics of community assembly on evolutionary time scales

Community assembly is studied using individual-based multispecies models. The models have stochastic population dynamics with mutation, migration, and extinction of species. Mutants appear as a result of mutation of the resident species, while migrants have no correlation with the resident species. It is found that the dynamics of community assembly with mutations are quite different from the case with migrations. In contrast to mutation models, which show intermittent dynamics of quasi-steady states interrupted by sudden reorganizations of the community, migration models show smooth and gradual renewal of the community. As a consequence, instead of the 1/f diversity fluctuations found for the mutation models, 1/f(2), random-walk like fluctuations are observed for the migration models. In addition, a characteristic species-lifetime distribution is found: a power law that is cut off by a "skewed" distribution in the long-lifetime regime. The latter has a longer tail than a simple exponential function, which indicates an age-dependent species-mortality function. Since this characteristic profile has been observed, both in fossil data and in several other mathematical models, we conclude that it is a universal feature of macroevolution.

Effects of demographic stochasticity on biological community assembly on evolutionary time scales

We study the effects of demographic stochasticity on the long-term dynamics of biological coevolution models of community assembly. The noise is induced in order to check the validity of deterministic population dynamics. While mutualistic communities show little dependence on the stochastic population fluctuations, predator-prey models show strong dependence on the stochasticity, indicating the relevance of the finiteness of the populations. For a predator-prey model, the noise causes drastic decreases in diversity and total population size. The communities that emerge under influence of the noise consist of species strongly coupled with each other and have stronger linear stability around the fixed-point populations than the corresponding noiseless model. The dynamics on evolutionary time scales for the predator-prey model are also altered by the noise. Approximate 1/f fluctuations are observed with noise, while 1/f2 fluctuations are found for the model without demographic noise.

Self-optimization, community stability, and fluctuations in two individual-based models of biological coevolution

We compare and contrast the long-time dynamical properties of two individual-based models of biological coevolution. Selection occurs via multispecies, stochastic population dynamics with reproduction probabilities that depend nonlinearly on the population densities of all species resident in the community. New species are introduced through mutation. Both models are amenable to exact linear stability analysis, and we compare the analytic results with large-scale kinetic Monte Carlo simulations, obtaining the population size as a function of an average interspecies interaction strength. Over time, the models self-optimize through mutation and selection to approximately maximize a community potential function, subject only to constraints internal to the particular model. If the interspecies interactions are randomly distributed on an interval including positive values, the system evolves toward self-sustaining, mutualistic communities. In contrast, for the predator-prey case the matrix of interactions is antisymmetric, and a nonzero population size must be sustained by an external resource. Time series of the diversity and population size for both models show approximate 1/f noise and power-law distributions for the lifetimes of communities and species. For the mutualistic model, these two lifetime distributions have the same exponent, while their exponents are different for the predator-prey model. The difference is probably due to greater resilience toward mass extinctions in the food-web like communities produced by the predator-prey model.

Individual-based predator-prey model for biological coevolution: fluctuations, stability, and community structure

We study an individual-based predator-prey model of biological coevolution, using linear stability analysis and large-scale kinetic Monte Carlo simulations. The model exhibits approximate 1/f noise in diversity and population-size fluctuations, and it generates a sequence of quasisteady communities in the form of simple food webs. These communities are quite resilient toward the loss of one or a few species, which is reflected in different power-law exponents for the durations of communities and the lifetimes of species. The exponent for the former is near -1 , while the latter is close to -2 . Statistical characteristics of the evolving communities, including degree (predator and prey) distributions and proportions of basal, intermediate, and top species, compare reasonably with data for real food webs.

The species–area relationship and evolution

Models relating to the species-area curve usually assume the existence of species, and are concerned mainly with ecological timescales. We examine an individual-based model of co-evolution on a spatial lattice based on the tangled nature model in which species are emergent structures, and show that reproduction, mutation and dispersion by diffusion, with interaction via genotype space, produces power-law species-area relations as observed in ecological measurements at medium scales. We find that long-lasting co-evolutionary habitats form, allowing high diversity levels in a spatially homogenous system.

Network properties, species abundance and evolution in a model of evolutionary ecology

We study the evolution of the network properties of a populated network embedded in a genotype space characterized by either a low or a high number of potential links, with particular emphasis on the connectivity and clustering. Evolution produces two distinct types of network. When a specific genotype is only able to influence a few other genotypes, the ecosystem consists of separate non-interacting clusters (i.e. isolated compartments) in genotype space. When different types may influence a large number of other sites, the network becomes one large interconnected cluster. The distribution of interaction strengths--but not the number of connections--changes significantly with time. We find that the species abundance is only realistic for a high level of species connectivity. This suggests that real ecosystems form one interconnected whole in which selection leads to stronger interactions between the different types. Analogies with niche and neutral theory and assembly models are also considered.

Stochastic models in population biology and their deterministic analogs

We introduce a class of stochastic population models based on "patch dynamics." The size of the patch may be varied, and this allows one to quantify the departures of these stochastic models from various mean-field theories, which are generally valid as the patch size becomes very large. These models may be used to formulate a broad range of biological processes in both spatial and nonspatial contexts. Here, we concentrate on two-species competition. We present both a mathematical analysis of the patch model, in which we derive the precise form of the competition mean-field equations (and their first-order corrections in the nonspatial case), and simulation results. These mean-field equations differ, in some important ways, from those which are normally written down on phenomenological grounds. Our general conclusion is that mean-field theory is more robust for spatial models than for a single isolated patch. This is due to the dilution of stochastic effects in a spatial setting resulting from repeated rescue events mediated by interpatch diffusion. However, discrete effects due to modest patch sizes lead to striking deviations from mean-field theory even in a spatial setting.