Punctuated equilibria and 1/f noise in a biological coevolution model with individual-based dynamics

We present a study by linear stability analysis and large-scale Monte Carlo simulations of a simple model of biological coevolution. Selection is provided through a reproduction probability that contains quenched, random interspecies interactions, while genetic variation is provided through a low mutation rate. Both selection and mutation act on individual organisms. Consistent with some current theories of macroevolutionary dynamics, the model displays intermittent, statistically self-similar behavior with punctuated equilibria. The probability density for the lifetimes of ecological communities is well approximated by a power law with exponent near -2, and the corresponding power spectral densities show 1/f noise (flicker noise) over several decades. The long-lived communities (quasi-steady states) consist of a relatively small number of mutualistically interacting species, and they are surrounded by a ``protection zone'' of closely related genotypes that have a very low probability of invading the resident community. The extent of the protection zone affects the stability of the community in a way analogous to the height of the free-energy barrier surrounding a metastable state in a physical system. Measures of biological diversity are on average stationary with no discernible trends, even over our very long simulation runs of approximately 3.4x10^7 generations.Comment: 20 pages RevTex. Minor revisions consistent with published versio

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 fitness 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.Comment: 26 pages, 12 figures. Discussion of early-time dynamics added. J. Math. Biol., in pres