Motional instabilities in prey-predator systems

A Comparison of the "Reduced Losses" and "Increased Production" Models for Mussel Bed Dynamics

Self-organised regular pattern formation is one of the foremost examples of the development of complexity in ecosystems. Despite the wide array of mechanistic models that have been proposed to understand pattern formation, there is limited general understanding of the feedback processes causing pattern formation in ecosystems, and how these affect ecosystem patterning and functioning. Here we propose a generalised model for pattern formation that integrates two types of within-patch feedback: amplification of growth and reduction of losses. Both of these mechanisms have been proposed as causing pattern formation in mussel beds in intertidal regions, where dense clusters of mussels form, separated by regions of bare sediment. We investigate how a relative change from one feedback to the other affects the stability of uniform steady states and the existence of spatial patterns. We conclude that there are important differences between the patterns generated by the two mechanisms, concerning both biomass distribution in the patterns and the resilience of the ecosystems to disturbances.

Dynamical study of infochemical influences on tropic interaction of diffusive plankton system

We propose a model for tropic interaction among the infochemical-producing phytoplankton and non-info chemical-producing phytoplankton and microzooplankton. Volatile information-conveying chemicals (infochemicals) released by phytoplankton play an important role in the food webs of marine ecosystems. Microzooplankton is an ecologically important grazer of phytoplankton for coexistence of a large number of phytoplankton species. Here, we discuss how information transferred by dimethyl sulfide shapes the interaction of phytoplankton. Phytoplankton deterrents may lead to propagation of IPP bloom. The interaction between IPP and microzooplankton follows the Beddington–DeAngelis-type functional response. Analytically, we discuss boundedness, stability and Turing instability of the model system. We perform numerical simulation for temporal (ODE model) as well as a spatial model system. Our numerical investigation shows that microzooplankton grazing refuse of IPP leads to oscillatory dynamics. Increasing diffusion coefficient of microzooplankton shows Turing instability. Time evolution also plays an important role in the stability of system dynamics. The results obtained in this paper are useful to understand the dominance of algal bloom in coastal and estuarine ecosystem.

Analysing diffusion and flow-driven instability using semidefinite programming

Diffusion and flow-driven instability, or transport-driven instability, is one of the central mechanisms to generate inhomogeneous gradient of concentrations in spatially distributed chemical systems. However, verifying the transport-driven instability of reaction-diffusion-advection systems requires checking the Jacobian eigenvalues of infinitely many Fourier modes, which is computationally intractable. To overcome this limitation, this paper proposes mathematical optimization algorithms that determine the stability/instability of reaction-diffusion-advection systems by finite steps of algebraic calculations. Specifically, the stability/instability analysis of Fourier modes is formulated as a sum-of-squares optimization program, which is a class of convex optimization whose solvers are widely available as software packages. The optimization program is further extended for facile computation of the destabilizing spatial modes. This extension allows for predicting and designing the shape of the concentration gradient without simulating the governing equations. The streamlined analysis process of self-organized pattern formation is demonstrated with a simple illustrative reaction model with diffusion and advection.

Influence of fast advective flows on pattern formation of Dictyostelium discoideum

We report experimental and numerical results on pattern formation of self-organizing Dictyostelium discoideum cells in a microfluidic setup under a constant buffer flow. The external flow advects the signaling molecule cyclic adenosine monophosphate (cAMP) downstream, while the chemotactic cells attached to the solid substrate are not transported with the flow. At high flow velocities, elongated cAMP waves are formed that cover the whole length of the channel and propagate both parallel and perpendicular to the flow direction. While the wave period and transverse propagation velocity are constant, parallel wave velocity and the wave width increase linearly with the imposed flow. We also observe that the acquired wave shape is highly dependent on the wave generation site and the strength of the imposed flow. We compared the wave shape and velocity with numerical simulations performed using a reaction-diffusion model and found excellent agreement. These results are expected to play an important role in understanding the process of pattern formation and aggregation of D. discoideum that may experience fluid flows in its natural habitat.

Turing patterns and long-time behavior in a three-species food-chain model

We consider a spatially explicit three-species food chain model, describing generalist top predator-specialist middle predator-prey dynamics. We investigate the long-time dynamics of the model and show the existence of a finite dimensional global attractor in the product space, L(2)(Ω). We perform linear stability analysis and show that the model exhibits the phenomenon of Turing instability, as well as diffusion induced chaos. Various Turing patterns such as stripe patterns, mesh patterns, spot patterns, labyrinth patterns and weaving patterns are obtained, via numerical simulations in 1d as well as in 2d. The Turing and non-Turing space, in terms of model parameters, is also explored. Finally, we use methods from nonlinear time series analysis to reconstruct a low dimensional chaotic attractor of the model, and estimate its fractal dimension. This provides a lower bound, for the fractal dimension of the attractor, of the spatially explicit model.

A mathematical biologist's guide to absolute and convective instability

Mathematical models have been highly successful at reproducing the complex spatiotemporal phenomena seen in many biological systems. However, the ability to numerically simulate such phenomena currently far outstrips detailed mathematical understanding. This paper reviews the theory of absolute and convective instability, which has the potential to redress this inbalance in some cases. In spatiotemporal systems, unstable steady states subdivide into two categories. Those that are absolutely unstable are not relevant in applications except as generators of spatial or spatiotemporal patterns, but convectively unstable steady states can occur as persistent features of solutions. The authors explain the concepts of absolute and convective instability, and also the related concepts of remnant and transient instability. They give examples of their use in explaining qualitative transitions in solution behaviour. They then describe how to distinguish different types of instability, focussing on the relatively new approach of the absolute spectrum. They also discuss the use of the theory for making quantitative predictions on how spatiotemporal solutions change with model parameters. The discussion is illustrated throughout by numerical simulations of a model for river-based predator-prey systems.

Spatial pattern formation in external noise: theory and simulation

Spatial pattern formation in fluctuating media is researched analytically from the point of view of the order parameters concept. A reaction-diffusion system with external noise is considered as a model of such media. Stochastic equations for unstable mode amplitudes (order parameters), the dispersion equation for averaged amplitudes of unstable modes, and the Fokker-Planck equation for the order parameters are obtained. The theory developed makes it possible to analyze different noise-induced effects including the variation of boundaries of ordering and disordering phase transitions depending on the parameters of external noise.

Propagation of Turing patterns in a plankton model

The paper is devoted to a reaction-diffusion system of equations describing phytoplankton and zooplankton distributions. Linear stability analysis of the model is carried out. Turing and Hopf stability boundaries are found. Emergence of two-dimensional spatial structures is illustrated by numerical simulations. Travelling waves between various stationary solutions are investigated. Transitions between homogeneous in space stationary solutions and Turing structures are studied.

Pattern formation in prey-taxis systems

In this paper, we consider spatial predator-prey models with diffusion and prey-taxis. We investigate necessary conditions for pattern formation using a variety of non-linear functional responses, linear and non-linear predator death terms, linear and non-linear prey-taxis sensitivities, and logistic growth or growth with an Allee effect for the prey. We identify combinations of the above non-linearities that lead to spatial pattern formation and we give numerical examples. It turns out that prey-taxis stabilizes the system and for large prey-taxis sensitivity we do not observe pattern formation. We also study and find necessary conditions for global stability for a type I functional response, logistic growth for the prey, non-linear predator death terms, and non-linear prey-taxis sensitivity.

The role of prey taxis in biological control: a spatial theoretical model

We study a reaction-diffusion-advection model for the dynamics of populations under biological control. A control agent is assumed to be a predator species that has the ability to perceive the heterogeneity of pest distribution. The advection term represents the predator density movement according to a basic prey taxis assumption: acceleration of predators is proportional to the prey density gradient. The prey population reproduces logistically, and the local population interactions follow the Holling Type II trophic function. On the scale of the population, our spatially explicit approach subdivides the predation process into random movement represented by diffusion, directed movement described by prey taxis, local prey encounters, and consumption modeled by the trophic function. Thus, our model allows studying the effects of large-scale predator spatial activity on population dynamics. We show under which conditions spatial patterns are generated by prey taxis and how this affects the predator ability to maintain the pest population below some economic threshold. In particular, intermediate taxis activity can stabilize predator-pest populations at a very low level of pest density, ensuring successful biological control. However, very intensive prey taxis destroys the stability, leading to chaotic dynamics with pronounced outbreaks of pest density.

Directed movement of predators and the emergence of density-dependence in predator-prey models

We consider a bitrophic spatially distributed community consisting of prey and actively moving predators. The model is based on the assumption that the spatial and temporal variations of the predators' velocity are determined by the prey gradient. Locally, the populations follow the simple Lotka-Volterra interaction. We also assume predator reproduction and mortality to be negligible in comparison with the time scale of migration. The model demonstrates heterogeneous oscillating distributions of both species, which occur because of the active movements of predators. One consequence of this heterogeneity is increased viability of the prey population, compared to the equivalent homogeneous model, and increased consumption. Further numerical analysis shows that, on the spatially aggregated scale, the average predator density adversely affects the individual consumption, leading to a nonlinear predator-dependent trophic function, completely different from the Lotka-Volterra rule assumed at the local scale.

Wave of chaos: new mechanism of pattern formation in spatio-temporal population dynamics

The dynamics of a simple prey-predator system is described by a system of two reaction- diffusion equations with biologically reasonable non-linearities (logistic growth of the prey, Holling type II functional response of the predator). We show that, when the local kinetics of the system is oscillatory, for a wide class of initial conditions the evolution of the system leads to the formation of a non-stationary irregular pattern corresponding to spatio-temporal chaos. The chaotic pattern first appears inside a sub-domain of the system. This sub-domain then steadily grows with time and, finally, the chaotic pattern invades the whole space, displacing the regular pattern.