Spatiotemporal dynamics of reaction–diffusion models of interacting populations

The present investigation deals with the necessary conditions for Turing instability with zero-flux boundary conditions that arise in a ratio-dependent predator–prey model involving the influence of logistic population growth in prey and intra-specific competition among predators described by a system of non-linear partial differential equations. The prime objective is to investigate the parametric space for which Turing spatial structure takes place and to perform extensive numerical simulation from both the mathematical and the biological points of view in order to examine the role of diffusion coefficients in Turing instability. Various spatiotemporal distributions of interacting species through Turing instability in two dimensional spatial domain are portrayed and analyzed at length in order to substantiate the applicability of the present model.     

Spike patterns in a reaction–diffusion ODE model with Turing instability

We explore a mechanism of pattern formation arising in processes described by a system of a single reaction–diffusion equation coupled with ordinary differential equations. Such systems of equations arise from the modeling of interactions between cellular processes and diffusing growth factors. We focus on the model of early carcinogenesis proposed by Marciniak‐Czochra and Kimmel, which is an example of a wider class of pattern formation models with an autocatalytic non‐diffusing component. We present a numerical study showing emergence of periodic and irregular spike patterns because of diffusion‐driven instability. To control the accuracy of simulations, we develop a numerical code on the basis of the finite‐element method and adaptive mesh grid. Simulations, supplemented by numerical analysis, indicate a novel pattern formation phenomenon on the basis of the emergence of nonstationary structures tending asymptotically to a sum of Dirac deltas. Copyright © 2013 John Wiley & Sons, Ltd.     

Turing pattern formation in a three species model with generalist predator and cross-diffusion

In a natural ecosystem, specialist predators feed almost exclusively on one species of prey. But generalist predators feed on many types of species. Consequently, their dynamics is not coupled to the dynamics of a specific prey population. However, the defense of prey formed by congregating made the predator tend to move in the direction of lower concentration of prey species. This is described by cross-diffusion in a generalist predator–prey model. First, the positive equilibrium solution is globally asymptotically stable for the ODE system and for the reaction–diffusion system without cross-diffusion, respectively, hence it does not belong to the classical Turing instability scheme. But it becomes linearly unstable only when cross-diffusion also plays a role. This implies that cross–diffusion can lead to the occurrence and disappearance of the instability. Our results exhibit some interesting combining effects of cross-diffusion, predations and intra-species interactions. Furthermore, we consider the existence and non-existence results concerning non-constant positive steady states (patterns) of the system. We demonstrate that cross-diffusion can create non-constant positive steady-state solutions.     

Characterization of spatial patterns produced by a Turing instability in coupled dynamical systems

We quantify the degree of spatial order of patterns at fixed time generated by lattices of coupled dynamical systems, using correlation-based and recurrence-based numerical diagnostics. These patterns are obtained through numerical integration of differential equations describing the interplay between activator and inhibitor species generating Turing patterns. We consider different types of coupling: linear (diffusive) interaction with nearest-neighbors, global (all-to-all) coupling and intermediate (nonlocal) coupling. Numerical simulations are performed in one and two spatial dimensions. The effects of noise are briefly discussed. We introduce a recurrence-based quantity (recurrence-rate matrix) to characterize two-dimensional spatial patterns.     

On nonlinear coupled diffusions in competition systems

A class of reaction-diffusion systems modeling plant growth with spatial competition in saturated media is presented. We show, in this context, that standard diffusion can not lead to pattern formation (Diffusion Driven Instability of Turing). Degenerated nonlinear coupled diffusions inducing free boundaries and exclusive spatial diffusions are proposed. Local and global existence results are proved for smooth approximations of the degenerated nonlinear diffusions systems which give rise to long-time pattern formations. Numerical simulations of a competition model with degenerate/non degenerate nonlinear coupled diffusions are performed and we carry out the effect of the these diffusions on pattern formation and on the change of basins of attraction.     

Turing instability for a ratio-dependent predator-prey model with diffusion

Ratio-dependent predator-prey models have been increasingly favored by field ecologists where predator-prey interactions have to be taken into account the process of predation search. In this paper we study the conditions of the existence and stability properties of the equilibrium solutions in a reaction-diffusion model in which predator mortality is neither a constant nor an unbounded function, but it is increasing with the predator abundance. We show that analytically at a certain critical value a diffusion driven (Turing type) instability occurs, i.e. the stationary solution stays stable with respect to the kinetic system (the system without diffusion). We also show that the stationary solution becomes unstable with respect to the system with diffusion and that Turing bifurcation takes place: a spatially non-homogenous (non-constant) solution (structure or pattern) arises. A numerical scheme that preserve the positivity of the numerical solutions and the boundedness of prey solution will be presented. Numerical examples are also included.     

Turing instability and Hopf bifurcation in a diffusive Leslie–Gower predator–prey model

In this paper, a reaction‐diffusion predator–prey system that incorporates the Holling‐type II and a modified Leslie‐Gower functional responses is considered. For ODE, the local stability of the positive equilibrium is investigated and the specific conditions are obtained. For partial differential equation, we consider the dissipation and persistence of solutions, the Turing instability of the equilibrium solutions, and the Hopf bifurcation. By calculating the normal form, we derive the formulae, which can determine the direction and the stability of Hopf bifurcation according to the original parameters of the system. We also use some numerical simulations to illustrate our theoretical results. Copyright © 2016 John Wiley & Sons, Ltd.     

Hopf bifurcation and Turing instability in spatial homogeneous and inhomogeneous predator-prey models

This paper is concerned with two-species spatial homogeneous and inhomogeneous predator-prey models with Beddington-DeAngelis functional response. For the spatial homogeneous model, the asymptotic behavior of the interior equilibrium and the existence of Hopf bifurcation of nonconstant periodic solutions surrounding the interior equilibrium are considered. Furthermore, the direction of Hopf bifurcation and the stability of bifurcated periodic solutions are investigated. For the model with no-flux boundary conditions, Turing instability of the interior equilibrium solution is studied. In particular, Turing instability region regarding the parameters is established. Finally, to verify our theoretical results, some numerical simulations are also included.     

Spatially discrete reaction–diffusion equations with discontinuous hysteresis

We address the question: Why may reaction–diffusion equations with hysteretic nonlinearities become ill-posed and how to amend this? To do so, we discretize the spatial variable and obtain a lattice dynamical system with a hysteretic nonlinearity. We analyze a new mechanism that leads to appearance of a spatio-temporal pattern called rattling: the solution exhibits a propagation phenomenon different from the classical traveling wave, while the hysteretic nonlinearity, loosely speaking, takes a different value at every second spatial point, independently of the grid size. Such a dynamics indicates how one should redefine hysteresis to make the continuous problem well-posed and how the solution will then behave. In the present paper, we develop main tools for the analysis of the spatially discrete model and apply them to a prototype case. In particular, we prove that the propagation velocity is of order $a{t}^{?1/2}$ as $t\to \infty$ and explicitly find the rate a.     

Bifurcation and Turing pattern formation in a diffusive ratio-dependent predator–prey model with predator harvesting

The ratio-dependent predator–prey model exhibits rich dynamics due to the singularity of the origin. Harvesting in a ratio-dependent predator–prey model is relatively an important research project from both ecological and mathematical points of view. In this paper, we study the temporal, spatial and spatiotemporal dynamics of a ratio-dependent predator–prey diffusive model where the predator population harvest at catch-per-unit-effort hypothesis. For the spatially homogeneous model, we derive conditions for determining the direction of Hopf bifurcation and the stability of the bifurcating periodic solution by the center manifold and the normal form theory. For the reaction–diffusion model, firstly it is shown that Turing (diffusion-driven) instability occurs, which induces spatial inhomogeneous patterns. Then it is demonstrated that the model exhibit Hopf bifurcation which produces temporal inhomogeneous patterns. Finally, the existence and non-existence of positive non-constant steady-state solutions are established. Moreover, numerical simulations are performed to visualize the complex dynamic behavior.     

Turing vegetation patterns in a generalized hyperbolic Klausmeier model

The formation of Turing vegetation patterns in flat arid environments is investigated in the framework of a generalized version of the hyperbolic Klausmeier model. The extensions here considered involve, on the one hand, the strength of the rate at which rainfall water enters into the soil and, on the other hand, the functional dependence of the inertial times on vegetation biomass and soil water. The study aims at elucidating how the inclusion of these generalized quantities affects the onset of bifurcation of supercritical Turing patterns as well as the transient dynamics observed from an uniformly vegetated state towards a patterned state. To achieve these goals, linear and multiple-scales weakly nonlinear stability analysis are addressed, this latter being inspected in both large and small spatial domains. Analytical results are then corroborated by numerical simulations, which also serve to describe more deeply the spatio-temporal evolution of the emerging patterns as well as to characterize the different timescales involved in vegetation dynamics.     

Direct discontinuous Galerkin method for nonlinear reaction–diffusion systems in pattern formation

Nonlinear reaction–diffusion systems are often employed in mathematical modeling for pattern formation. Most of the work to date has been concerned within one-dimensional or rectangular domains. However, it is recognised that in most applications multidimensional complex geometrical domains are typically more important. In this paper we solve reaction–diffusion systems by combining direct discontinuous Galerkin (DDG) finite element methods with implicit integration factor (IIF) time integration method, on triangular meshes. This allows us solve the nonlinear algebraic systems on an element-by-element bases with significant gains in computational time. Numerical solutions of two reaction–diffusion systems, the well-studied Schnakenberg model and chloride–iodide–malonic acid (CIMA) reactive model, are presented to demonstrate effects of various domain geometries on the resulting biological patterns. Our numerical results are in good agreement with other numerical and analytical results, and with experimental results.     

Group classification of systems of non-linear reaction-diffusion equations with general diffusion matrix. II. Generalized Turing systems

Group classification of systems of two coupled non-linear reaction-diffusion equation with a diagonal diffusion matrix is carried out. Symmetries of diffusion systems with singular diffusion matrix and additional first order derivative terms are described.     

Stability and pattern formation in a coupled arbitrary order of autocatalysis system

Spatiotemporal structures arising in two identical cells, each governed by arbitrary order autocatalator kinetics and coupled via the diffusive interchange of a reactant, are discussed. The stability of two homogeneous steady states is obtained with the use of linear stability analysis. By studying the linearized equations, it is found that two steady states, in the uncoupled and coupled system respectively, may give rise to the possibility of bifurcations to spatially nonuniform pattern forms. Further information about Turing bifurcation solutions close to these bifurcation points are obtained by weakly nonlinear theory. It is seen that the coupling leads to bifurcations not present in the uncoupled system which give rise to locally stable nonuniform pattern forms. Finally the stability of the equilibrium points of the amplitude equation is discussed by weakly nonlinear theory, with the bifurcation branches about small coupled system with 0<α?1$0<\alpha ?1$ and large coupling for α?1$\alpha ?1$.     

Spatiotemporal dynamics of a predator-prey model incorporating a prey refuge

In this paper, we investigate the spatiotemporal dynamics of a ratio-dependent predator-prey model with cross diffusion incorporating proportion of prey refuge. First we get the critical lines of Hopf and Turing bifurcations in a spatial domain by using mathematical theory. More specifically, the exact Turing region is given in a two parameter space. Also we perform a series of numerical simulations. The obtained results reveal that this system has rich dynamics, such as spotted, stripe and labyrinth patterns which show that it is useful to use the predator-prey model to reveal the spatial dynamics in the real world.     

Modeling and analysis of a marine plankton system with nutrient recycling and diffusion

This article describes a nutrient‐phytoplankton‐zooplankton system with nutrient recycling in the presence of toxicity. We have studied the dynamical behavior of the system with delayed nutrient recycling in the first part of the article. Uniform persistent of the system is examined. In the second part of the article, we have incorporated diffusion of the plankton population to the system and dynamical behavior of the system is analyzed with instantaneous nutrient recycling. The condition of the diffusion driven instability is obtained. The conditions for the occurrence of Hopf and Turing bifurcation critical line in a spatial domain are derived. Variation of the system with small periodicity of diffusive coefficient has been studied. Stability condition of the plankton system subject to the periodic diffusion coefficient of the zooplankton is derived. It is observed that nutrient‐phytoplankton‐zooplankton interactions are very complex and situation specific. Moreover, we have obtained different exciting results, ranging from stable situation to cyclic oscillatory behavior may occur under different favorable conditions, which may give some insights for predictive management. © 2014 Wiley Periodicals, Inc. Complexity 21: 229–241, 2015     

Turing pattern formation in a predator–prey system with cross diffusion

The paper explores the impacts of cross-diffusion on the formation of spatial patterns in a ratio-dependent predator–prey system with zero-flux boundary conditions. Our results show that under certain conditions, cross-diffusion can trigger the emergence of spatial patterns which is however impossible under the same conditions when cross-diffusion is absent. We give a rigorous proof that the model has at least one spatially heterogenous steady state by means of the Leray–Schauder degree theory. In addition, numerical simulations are performed to visualize the complex spatial patterns.