Discrete Dynamics in Nature and Society

Volume 2013 (2013), Article ID 454209, 8 pages

http://dx.doi.org/10.1155/2013/454209

Research Article

Pattern Formation in a Diffusive Ratio-Dependent Holling-Tanner Predator-Prey Model with Smith Growth

Bo Yang

School of Electronic and Information Engineering, Gansu Lianhe University, Lanzhou 730000, China

Received 3 February 2013; Accepted 4 March 2013

Academic Editor: Yonghui Xia

Copyright © 2013 Bo Yang. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

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.

1. Introduction

The problem of pattern formation is, perhaps, the most challenging in modern ecology, biology, chemistry, and many other fields of science [1]. 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 [2].

Turing [3] 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 [4] first used reaction-diffusion system to explain pattern formation in ecological context based upon the seminal work by Turing [3]. 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 [24]. This is supported by numerous fields and laboratory experiments and observations [25, 26]. In [24], 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 [27] 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 [28] 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 [29], 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

Theorem 1. 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 1: Phase portraits of model (1). The parameters are taken as , , , , , , . is a saddle point. is locally asymptotically stable. The dashed curve is the -nullcline, and the dotted vertical line is the -nullcline.

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.

Figure 2: Phase portraits of model (1). The parameters are taken as , , , , , , . Model (1) enters into Hopf bifurcation around , and there is a limit cycle. The dashed curve is the -nullcline, and the dotted vertical line is the -nullcline.

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 [35] 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 [36] 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)).

(a)

(b)

(c)

(d)

Figure 3: Stripes pattern formation for model (4) by taking . Other parameters are fixed as (33). Times: (a) 0; (b) 50; (c) 250; (d) 2500.

In Figure 4, with , we show the stripes-spots mixtures pattern for model (4).

(a)

(b)

(c)

(d)

Figure 4: Stripes-spots mixtures pattern formation for model (4) by taking . Other parameters are fixed as (33). Times: (a) 0; (b) 50; (c) 250; (d) 2500.

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.

(a)

(b)

(c)

(d)

Figure 5: Spots pattern formation for model (4) by taking . Other parameters are fixed as (33). Times: (a) 0; (b) 50; (c) 250; (d) 2500.

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.

References

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.

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.

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.

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.

L. S. Luckinbill, “The effects of space and enrichment on a predator-prey system,” Ecology, vol. 55, no. 5, pp. 1142–1147, 1974.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

F. E. Smith, “Population dynamics in Daphnia Magna and a new model for population growth,” Ecology, vol. 44, pp. 651–663, 1963.

E. C. Pielou, An Introduction to Mathematical Ecology, John Wiley & Sons, New York, NY, USA, 1969.

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.

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.

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.

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.

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.

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.