Theory and Applications of Periodic Solutions and Almost Periodic SolutionsView this Special Issue
Pattern Formation in a Diffusive Ratio-Dependent Holling-Tanner Predator-Prey Model with Smith Growth
The spatiotemporal dynamics of a diffusive ratio-dependent Holling-Tanner predator-prey model with Smith growth subject to zero-flux boundary condition are investigated analytically and numerically. The asymptotic stability of the positive equilibrium and the existence of Hopf bifurcation around the positive equilibrium are shown; the conditions of Turing instability are obtained. And with the help of numerical simulations, it is found that the model exhibits complex pattern replication: stripes, spots-stripes mixtures, and spots Turing patterns.
The problem of pattern formation is, perhaps, the most challenging in modern ecology, biology, chemistry, and many other fields of science . Patterns generated in abiotically homogeneous environments are particularly interesting because they require an explanation based on the individual behavior of organisms. They are commonly called “emergent patterns,” because they emerge from interactions in spatial scales that are much larger than the characteristic scale of individuals .
Turing  showed how the coupling of reaction and diffusion can induce instability and pattern formation. Turing’s revolutionary idea was that the passive diffusion could interact with chemical reaction in such a way that even if the reaction by itself has no symmetry-breaking capabilities, diffusion can destabilize the symmetry so that the system with diffusion can have them. Segel and Jackson  first used reaction-diffusion system to explain pattern formation in ecological context based upon the seminal work by Turing . Since then, a lot of studies have been devoted to spatiotemporal patterns which were produced by reaction-diffusion predator-prey, models with either a prey-dependent or a ratio-dependent predator functional response, for example, [1, 2, 5–20] and references cited therein.
Recently, there is a growing explicit biological and physiological evidence [21–23] that in many situations, especially, when the predator has to search for food (and therefore has to share or compete for food), a more suitable general predator-prey theory should be based on the so-called ratio-dependent function which can be roughly stated as that the per capital predator growth rate should be a function of the ratio of prey to predator abundance, and so would be the so-called predator functional responses . This is supported by numerous fields and laboratory experiments and observations [25, 26]. In , the authors investigated the effect of time delays on the stability of the model and discussed the local asymptotic stability and the Hopf bifurcation. Liang and Pan  have studied the local and global asymptotic stability of the coexisting equilibrium point and obtained the conditions for Poincare-Andronov-Hopf-bifurcating periodic solution. M. Banerjee and S. Banerjee  have studied the local asymptotic stability of the equilibrium point and obtained the conditions for the occurrence of Turing-Hopf instability for reaction-diffusion model. It is shown that prey and predator populations exhibit spatiotemporal patterns resulting from temporal oscillation of both the population and spatial instability.
Besides, in , Smith has shown that the logistic equation is not realistic for a food-limited population under the effects of environmental toxicants and established a new growth function—Smith growth function. And it has been proposed by several authors [29–34] for the dynamics of a population where the growth limitations are based on the proportion of available resources not utilized. However, pattern formation in the case of Holling-Tanner type predator-prey models with ratio-dependent functional response and Smith growth still remains an interesting area of research.
In this present work, we will focus on the ratio-dependent Holling-Tanner model with Smith growth for predator-prey interaction where random movement of both species is taken into account. The rest of the paper is organized as follows. In Section 2, we establish the ratio-dependent Holling-Tanner predator-prey model with Smith growth and study the local asymptotic stability of the positive equilibrium, existence of Hopf bifurcation around the positive equilibrium, and the conditions for the occurrence of Turing instability. In Section 3, we present and discuss the results of pattern formation via numerical simulation, which is followed by the last section, that is, conclusions and discussions.
2. The Model and the Linear Stability Analysis
2.1. The Model
In this paper, we rigorously consider the radio-dependent Holling-Tanner predator-prey model with Smith growth taking the form: where and stand for prey and predator population (density) at any instant of time . , , , , , are positive constants that stand for prey intrinsic growth rate, carrying capacity, capturing rate, half capturing saturation constant, predator intrinsic growth rate, conversion rate of prey into predators biomass, respectively. And is the replacement of mass in the population at . The model with Smith growth takes into account both environmental and food chain effects of toxicant stress.
From the standpoint of biology, we are interested only in the dynamics of model (1) in the closed first quadrant . Thus, we consider only the biologically meaningful initial conditions which are continuous functions due to their biological sense.
Straightforward computation shows that model (1) is continuous and Lipschitzian in if we redefine Hence, the solution of model (1) with positive initial conditions exists and is unique.
Also considering the spatial dispersal and environmental heterogeneity, in this paper we study the diffusive Holling-Tanner model obtained from the temporal model (1) by incorporating diffusion terms as follows: where the nonnegative constants and are the diffusion coefficients of and , respectively. , the usual Laplacian operator in two-dimensional space, is used to describe the Brownian random motion.
Model (4) is to be analyzed under the following nonzero initial conditions: and zero-flux boundary conditions: where is a bounded domain with a smooth boundary and is the outward unit normal vector on . The zero-flux boundary condition indicates that predator-prey system is self-contained with zero population flux across the boundary.
2.2. The Stability of the Nonspatial Model (1)
In this subsection, we restrict ourselves to the stability analysis of the nonspatial model (1). It is easy to verify that model (1) has a trivial equilibrium point . Simple computation shows that if , model (1) possess a unique positive equilibrium, denoted by , where
The Jacobian matrix at is Clearly, is a saddle point.
In the following, we will discuss the stability of the positive equilibrium of model (1). The Jacobian matrix at is given by where Then we can get
The positive equilibrium is locally asymptotically stable if and only if
The positive equilibrium is unstable if and only if
The model enters into a Hopf-bifurcation around at , where satisfies the equality
Proof. (i) If , then . Thus, the equilibrium point is locally asymptotically stable, similar to the proof of (ii).
(iii) A Hopf bifurcation occurs if and only if there exists a such that where is a root of the characteristic equation of : The condition gives . Thus for , both eigenvalues will be purely imaginary and there are no other eigenvalues with negative real part. Now we verify the transversality condition .
Substituting into the equation and separating real and imaginary parts we obtain Differentiating (17) both sides with respect to , we get where . Thus, we obtain which verify the transversality condition. Hence, the system undergoes a Hopf bifurcation at as passes through the value . This ends the proof.
In Figure 1, we show the phase portraits of (1) with , , , , , , and . The horizontal axis is the prey population , and the vertical axis is the predator population . The dashed curve is the -nullcline, and the dotted vertical line is the -nullcline. It is easy to see that the equilibrium is a saddle and is locally asymptotically stable.
Figure 2 illustrates a Hopf-bifurcation situation of the model around for . In this case, limit cycle arising through Hopf bifurcation is a stable limit cycle which attracts all trajectories starting from a point in the interior of first quadrant.
2.3. The Stability of the Spatial Model (4)
In this subsection, we will focus on the effect of diffusion on the model system about the positive equilibrium.
Now, we study the nonlinear evolution of a perturbation around . The corresponding linearized model (4) then takes the form where , are defined the same as (10).
We use to denote a column vector, and let
Let and Then forms a basis of the space of functions in that satisfy zero-flux boundary conditions. We look for a normal mode corresponding to model (21) as following form: where is a vector depending on . Plugging (24) into model (21) yields where . A nontrivial normal mode can be obtained by setting This leads to the following dispersion formula for : where
Mathematically speaking, a positive equilibrium of model (4) is Turing unstable, which means that it is an asymptotically stable steady-state solution of the model (1) without diffusion but is unstable with respect to the solutions of the model (4) with diffusion.
Therefore, the Turing instability sets in when at least one of the following conditions is violated: But it is evident that is not violated if . Hence only the violation of condition gives rise to diffusion instability. As a consequence, a necessary condition is Otherwise for all . For instability we must have for some , and we notice that achieves its minimum: at the critical value when
Summarizing the previous calculation, we conclude the following theorem.
Theorem 2. Assume that(A1),(A2) + + + + 2 + + − .
Then the positive equilibrium of model (4) is Turing unstable.
From Theorem 2, we can know that there is Turing instability in model (4) if conditions (A1) and (A2) hold. In this situation, the solutions to model (4) may be unstable and Turing patterns can emerge in the model.
3. Turing Pattern Formation
In this section, we perform extensive numerical simulations of the spatially extended model (4) in two-dimensional space, and the qualitative results are shown here. All our numerical simulations employ the zero-flux boundary conditions with a system size of . Other parameters are set as
The numerical integration of model (4) is performed by using a finite difference approximation for the spatial derivatives and an explicit Euler method for the time integration  with a time step size of . The initial condition is always a small amplitude random perturbation around the positive constant steady-state solution . After the initial period during which the perturbation spread, the model either goes to a time-dependent state or to an essentially steady-state solution (time independent).
We use the standard five-point approximation  for the 2D Laplacian with the zero-flux boundary conditions. More precisely, the concentrations at the moment at the mesh position are given by with the Laplacian defined by where , , and the space step size .
In the numerical simulations, different types of dynamics are observed and it is found that the distributions of predator and prey are always of the same type. Consequently, we can restrict our analysis of pattern formation to one distribution. In this section, we show the distribution of prey , for instance. We have taken some snapshots with red (blue) corresponding to the high (low) value of prey .
Now, we show the Turing patterns for the different values of the control parameter . Via numerical simulations, one can see that the model dynamics exhibits spatiotemporal complexity of pattern formation, including stripes, stripes-spots mixtures, and spots Turing patterns.
In Figure 3, with , starting with a homogeneous state (cf. Figure 3(a)), the random perturbations lead to the formation to stripes spots (cf. Figure 3(c)), and the latter random perturbations make these spots decay, ending with the time-independent stripes pattern (cf. Figure 3(d)).
In Figure 4, with , we show the stripes-spots mixtures pattern for model (4).
Figure 5 shows the time process of spots pattern formation of prey for . In this case, the pattern takes a long time to settle down, starting with a homogeneous state (cf. Figure 5(a)), and the random perturbations lead to the formation of stripes and spots (cf. Figure 5(b)), ending with spots only (cf. Figure 5(d))—the prey is isolated zones with high population density, and the remainder region is of low density.
From Figure 3 to Figure 5, we can see that, on increasing the control parameter from to , the pattern sequence “stripe → stripes-spots mixtures → spots” can be observed.
4. Conclusions and Remarks
In summary, in this paper, we have investigated the spatiotemporal dynamics of a diffusive predator-prey model where the interaction between prey and predator follows Holling-Tanner formulation with ratio-dependent functional response and Smith growth. The value of this study is threefold. First, it presents the conditions for the stability of the equilibrium and the existence of Hopf bifurcation for the nonspatial model. Second, it rigorously proves Turing instability by linear stability analysis for the spatial model. Third, it illustrates the Turing pattern formation via numerical simulations, which shows that the spatial model dynamics exhibits complex pattern replication.
By a series of numerical simulations, we find that the spatial model (4) has rich Turing pattern replications, such as stripes, stripes-spots mixtures, and spots patterns. In the viewpoint of population ecology, in the case of stripe pattern (cf. Figure 3), the prey is the isolated “stripes-like region” with high density, and the remainder stripes-like region is of low density. And in the case of spots pattern (cf. Figure 5), the prey is the isolated “cycle region” with high density, and the remainder region is of low density, which is larger than the “spots” region.
S. V. Petrovskii and H. Malchow, “A minimal model of pattern formation in a prey-predator system,” Mathematical and Computer Modelling, vol. 29, no. 8, pp. 49–63, 1999.View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
D. Alonso, F. Bartumeus, and J. Catalan, “Mutual interference between predators can give rise to turing spatial patterns,” Ecology, vol. 83, no. 1, pp. 28–34, 2002.View at: Google Scholar
A. M. Turing, “The chemical basis of morphogenesis,” Philosactions Transactions of the Royal Society of London. Series B, vol. 237, no. 641, pp. 37–72, 1952.View at: Google Scholar
L. A. Segel and J. L. Jackson, “Dissipative structure: an explanation and an ecological example,” Journal of Theoretical Biology, vol. 37, no. 3, pp. 545–559, 1972.View at: Google Scholar
L. S. Luckinbill, “The effects of space and enrichment on a predator-prey system,” Ecology, vol. 55, no. 5, pp. 1142–1147, 1974.View at: Google Scholar
S. A. Levin, “The problem of pattern and scale in ecology: the Robert H. MacArthur award lecture,” Ecology, vol. 73, no. 6, pp. 1943–1967, 1992.View at: Google Scholar
M. Pascual, “Diffusion-induced chaos in a spatial predator-prey system,” Proceedings of the Royal Society B, vol. 251, no. 1330, pp. 1–7, 1993.View at: Google Scholar
A. B. Medvinsky, S. V. Petrovskii, I. A. Tikhonova, H. Malchow, and B.-L. Li, “Spatiotemporal complexity of plankton and fish dynamics,” SIAM Review, vol. 44, no. 3, pp. 311–370, 2002.View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
J. D. Murray, Mathematical Biology. II: Spatial Models and Biomedical Applications, vol. 18 of Interdisciplinary Applied Mathematics, Springer, New York, NY, USA, 3rd edition, 2003.View at: MathSciNet
A. B. Peet, P. A. Deutsch, and E. Peacock-López, “Complex dynamics in a three-level trophic system with intraspecies interaction,” Journal of Theoretical Biology, vol. 232, no. 4, pp. 491–503, 2005.View at: Publisher Site | Google Scholar | MathSciNet
D. A. Griffith and P. R. Peres-Neto, “Spatial modeling in ecology: the flexibility of eigenfunction spatial analyses,” Ecology, vol. 87, no. 10, pp. 2603–2613, 2006.View at: Publisher Site | Google Scholar
I. G. Pearce, M. A. J. Chaplain, P. G. Schofield, A. R. A. Anderson, and S. F. Hubbard, “Modelling the spatio-temporal dynamics of multi-species host-parasitoid interactions: heterogeneous patterns and ecological implications,” Journal of Theoretical Biology, vol. 241, no. 4, pp. 876–886, 2006.View at: Publisher Site | Google Scholar | MathSciNet
W. Wang, Q.-X. Liu, and Z. Jin, “Spatiotemporal complexity of a ratio-dependent predator-prey system,” Physical Review E, vol. 75, no. 5, Article ID 051913, 9 pages, 2007.View at: Publisher Site | Google Scholar | MathSciNet
W. Wang, L. Zhang, H. Wang, and Z. Li, “Pattern formation of a predator-prey system with Ivlev-type functional response,” Ecological Modelling, vol. 221, no. 2, pp. 131–140, 2010.View at: Publisher Site | Google Scholar
W. Wang, Y. Lin, F. Rao, L. Zhang, and Y. Tan, “Pattern selection in a ratio-dependent predator-prey model,” Journal of Statistical Mechanics: Theory and Experiment, vol. 2010, Article ID P11036, 2010.View at: Google Scholar
W. Wang, Y. Lin, L. Zhang, F. Rao, and Y. Tan, “Complex patterns in a predator-prey model with self and cross-diffusion,” Communications in Nonlinear Science and Numerical Simulation, vol. 16, no. 4, pp. 2006–2015, 2011.View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
X. Guan, W. Wang, and Y. Cai, “Spatiotemporal dynamics of a Leslie-Gower predator-prey model incorporating a prey refuge,” Nonlinear Analysis: Real World Applications, vol. 12, no. 4, pp. 2385–2395, 2011.View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
Y. Zhu, Y. Cai, S. Yan, and W. Wang, “Dynamical analysis of a delayed reaction-diffusion predator-prey system,” Abstract and Applied Analysis, vol. 2012, Article ID 323186, 23 pages, 2012.View at: Google Scholar | Zentralblatt MATH | MathSciNet
Y. Cai, W. Wang, and J. Wang, “Dynamics of a diffusive predator-prey model with additive Allee effect,” International Journal of Biomathematics, vol. 5, no. 2, Article ID 1250023, 11 pages, 2012.View at: Publisher Site | Google Scholar | MathSciNet
W. Wang, Z. Guo, R. K. Upadhyay, and Y. Lin, “Pattern formation in a cross-diffusive Holling-Tanner model,” Discrete Dynamics in Nature and Society, vol. 2012, Article ID 828219, 12 pages, 2012.View at: Publisher Site | Google Scholar
R. Arditi, L. R. Ginzburg, and H. R. Akcakaya, “Variation in plankton densities among lakes: a case for ratio-dependent models,” The American Naturalist, vol. 138, no. 5, pp. 287–1296, 1991.View at: Google Scholar
R. Arditi and H. Saiah, “Empirical evidence of the role of heterogeneity in ratio-dependent consumption,” Ecology, vol. 73, no. 5, pp. 1544–1551, 1992.View at: Google Scholar
A. P. Gutierrez, “Physiological basis of ratio-dependent predator-prey theory: the metabolic pool model as a paradigm,” Ecology, vol. 73, no. 5, pp. 1552–1563, 1992.View at: Google Scholar
T. Saha and C. Chakrabarti, “Dynamical analysis of a delayed ratio-dependent Holling-Tanner predator-prey model,” Journal of Mathematical Analysis and Applications, vol. 358, no. 2, pp. 389–402, 2009.View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
R. Arditi and L. R. Ginzburg, “Coupling in predator-prey dynamics: ratio-dependence,” Journal of Theoretical Biology, vol. 139, no. 3, pp. 311–326, 1989.View at: Google Scholar
R. Arditi, N. Perrin, and H. Saiah, “Functional responses and heterogeneities: an experimental test with cladocerans,” Oikos, vol. 60, no. 1, pp. 69–75, 1991.View at: Google Scholar
Z. Liang and H. Pan, “Qualitative analysis of a ratio-dependent Holling-Tanner model,” Journal of Mathematical Analysis and Applications, vol. 334, no. 2, pp. 954–964, 2007.View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
M. Banerjee and S. Banerjee, “Turing instabilities and spatio-temporal chaos in ratio-dependent Holling-Tanner model,” Mathematical Biosciences, vol. 236, no. 1, pp. 64–76, 2012.View at: Publisher Site | Google Scholar | MathSciNet
F. E. Smith, “Population dynamics in Daphnia Magna and a new model for population growth,” Ecology, vol. 44, pp. 651–663, 1963.View at: Google Scholar
E. C. Pielou, An Introduction to Mathematical Ecology, John Wiley & Sons, New York, NY, USA, 1969.View at: MathSciNet
T. G. Hallam and J. T. Deluna, “Effects of toxicants on populations: a qualitative approach III,” Journal of Theoretical Biology, vol. 109, no. 3, pp. 411–429, 1984.View at: Google Scholar
K. Gopalsamy, M. R. S. Kulenović, and G. Ladas, “Environmental periodicity and time delays in a “food-limited” population model,” Journal of Mathematical Analysis and Applications, vol. 147, no. 2, pp. 545–555, 1990.View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
W. Feng and X. Lu, “On diffusive population models with toxicants and time delays,” Journal of Mathematical Analysis and Applications, vol. 233, no. 1, pp. 373–386, 1999.View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
M. Fan and K. Wang, “Periodicity in a “food-limited” population model with toxicants and time delays,” Acta Mathematicae Applicatae Sinica, vol. 18, no. 2, pp. 309–314, 2002.View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet
M. R. Garvie, “Finite-difference schemes for reaction-diffusion equations modeling predator-prey interactions in MATLAB,” Bulletin of Mathematical Biology, vol. 69, no. 3, pp. 931–956, 2007.View at: Publisher Site | Google Scholar | MathSciNet
A. Munteanu and R. V. Solé, “Pattern formation in noisy self-replicating spots,” International Journal of Bifurcation and Chaos, vol. 16, no. 12, pp. 3679–3685, 2006.View at: Publisher Site | Google Scholar