Abstract

The reaction-diffusion Holling-Tanner prey-predator model considering the Allee effect on predator, under zero-flux boundary conditions, is discussed. Some properties of the solutions, such as dissipation and persistence, are obtained. Local and global stability of the positive equilibrium and Turing instability are studied. With the help of the numerical simulations, the rich Turing patterns, including holes, stripes, and spots patterns, are obtained.

1. Introduction

The Holling-Tanner prey-predator model is an important and interesting predator-prey model in both biological and mathematical sense [14]. The reaction-diffusion Holling-Tanner prey-predator model takes the following form: where and represent population density of prey and predator at time , respectively. The parameters ,  ,  ,  ,  , and   are all positive. stands for the intrinsic growth rate of prey, is the prey carrying capacity, is the maximum predation rate, is the self-saturation prey density, is predator intrinsic growth rate, and is conversion rate of prey into predator biomass. and are the diffusion coefficients of and , respectively, and we always assume that ,  . is the usual Laplacian operator in 2-dimensional space. is a bounded domain with smooth boundary .

Set For the sake of convenience, we still use variables , instead of , . Thus, model (1) is converted into where the new parameters are

The dynamics of the reaction-diffusion Holling-Tanner prey-predator model has proven quite interesting and received intensive study by both ecologists and mathematicians in many articles, see, for example, [510] and the references therein. Peng and Wang [5, 6] analyzed the global stability of the unique positive constant steady state and established the results for the existence and nonexistence of positive nonconstant steady states; Wang et al. [7] studied positive steady-state solutions and investigated the appearance of sharp spatial patterns arising from the model. Shi and coworkers [8] studied the global attractor and persistence property, local and global asymptotic stability of the unique positive constant equilibrium, and the existence and nonexistence of nonconstant positive steady states; Li et al. [9] considered the Turing and Hopf bifurcations of the equilibrium solutions; Liu and Xue [10] found the model exhibits the spotted, black-eye, and labyrinthine patterns.

On the other hand, in population dynamics, Allee effect [11] is an ecological phenomenon caused by any mechanism leading to a positive relationship between a component of individual fitness and either the number or density of conspecific [1215]. In ecological studies, the understanding of the influence of Allee effect plays a central role since the Allee effect can greatly increase the likelihood of local and global extinctions [16, 17]. As a result both ecologists and mathematicians are interested in Allee effect in the predator-prey model, and much progress has been seen in the study of Allee effect, see [14, 1827], and many more investigations were done in recent years. But there have been few papers discussing the impact of the Allee effect on predator in the predator-prey models.

Based on the previous discussion, in the present paper we adopt the reaction-diffusion Holling-Tanner prey-predator model with Allee effect on predator.

If we assume that the predator population is subject to an Allee effect, taking into account zero-flux boundary conditions, model (3) can be rewritten as where is the term for the Allee effect, and can be defined as the Allee effect constant. The per capita growth rate of the predator is reduced from to due to the Allee effect [17, 28]. is the outward unit normal vector on , and the zero-flux boundary conditions mean that model (5) is self-contained and has no population flux across the boundary [29, 30]. The initial data and are continuous functions on .

The corresponding kinetic equation to model (5) is

The plan of the paper is as follows. Section 2 is dedicated to furnish some properties of the solutions concerning the mathematical model used. In Section 3, the local and global stability of the positive equilibrium of the model is considered. Section 4 is devoted to the diffusion-driven instability (Turing effect) and illustrates the different Turing patterns by using the numerical simulations. Finally, in Section 5, some conclusions and discussions are given.

2. Large Time Behavior of Solution to Model (5)

In this section, we give some properties of the solutions, and these results will be often used later.

2.1. Dissipation

Theorem 1. For any solution of model (5) with nonnegative initial conditions, then Hence, for any , the rectangle is a global attractor of model (5) in .

Proof. satisfies
Let be a solution of the ordinary differential equation, Then, . From the comparison principle, one can get ; hence,
As a result, for any , there exists a , such that for all and . Similarly, satisfies Thus, . From the arbitrariness of , we can get that This ends the proof.

2.2. Persistence

Definition 2 (see [31]). Model (5) is said to have the persistence property if for any nonnegative initial data , there exists a positive constant , such that the corresponding solution of model (5) satisfies
In the following, we will show that model (5) is persistent. From the viewpoint of biology, this implies that the two species of prey and predator will always coexist at any time and any location of the inhabit domain, no matter what their diffusion coefficients are, under certain conditions on parameters.

Theorem 3. If and , then model (5) has the persistence property.

Proof. The proof is based on comparison principles. From (12), for , it is clear that there exists a , such that for all and . Hence, is an upper solution of the following problem:
Let be the unique positive solution to the following problem:
Since , then . By comparison, it follows that . This implies that Hence, for and .
Similarly, by the second equation in model (5), we have that is an upper solution of problem
Let be the unique positive solution to the following problem: Since , there exists a , such that . Hence, we have for the arbitrariness of , and an application of the comparison principle gives The proof is complete.

3. Stability

In this section, we will devote consideration to the stability of the positive equilibrium for model (5).

Clearly, model (5) has a unique positive equilibrium , where and with .

For the sake of simplicity, we rewrite model (5) as the vectorial form where ,  , and

Let be the eigenvalues of the operator on with the zero-flux boundary conditions. And set   is an orthonormal basis of , and , then where .

The linearization of model (5) at the positive equilibrium can be expressed by where

From [32], it is known that if all the eigenvalues of the operator have negative real parts, then is asymptotically stable; if there is an eigenvalue with positive real part, then is unstable; if all the eigenvalues have nonpositive real parts while some eigenvalues have zero real part, then the stability of cannot be determined by the linearization.

For each , is invariant under the operator and is an eigenvalue of if and only if is an eigenvalue of the matrix for some .

So, the local stability of the positive equilibrium can be analyzed as follows.

Theorem 4. Assume that and the first eigenvalue subject to the zero-flux boundary conditions satisfies Then, the positive equilibrium is uniformly asymptotically stable.

Proof. The stability of the positive equilibrium is reduced to consider the characteristic equation with
In view of , it follows that Remark that for any , we have .
In view of the relation , one can calculate that Recall that , we conclude that for all .
Therefore, the eigenvalues of the matrix have negative real parts. It thus follows from the Routh-Hurwitz criterion that, for each , the two roots and of all have negative real parts.
In the following, we prove that there exists such that
Let , then Since as , it follows that
By the Routh-Hurwitz criterion, it follows that the two roots , of all have negative real parts. Thus, there exists a positive constant , such that . By continuity, we see that there exists such that the two roots , of satisfy ,  , for all . In turn, , for all . Let Then and (33) holds for .
Consequently, the spectrum of , which consists of eigenvalues, lies in . In the sense of [32], we obtain that the positive constant steady-state solution of model (5) is uniformly asymptotically stable. This ends the proof.

In the following, we shall prove that the positive equilibrium of model (5) is globally asymptotically stable.

Theorem 5. Suppose that , , and . The positive equilibrium of model (5) is globally asymptotically stable, if(a1),(a2), where

Proof. We adopt the Lyapunov function where , . It can be easily verified that the function is zero at the positive equilibrium and is positive for all other positive values of and .
Then, where
It is obvious that if is negative definite.
can be expressed in a quadratic form , where
is negative definite if the symmetric matrices is positive. It can be easily shown that the symmetric matrix is positive definite if the following conditions are true:(i),(ii).
Proof of (i). Applying Theorems 1 and 3, we get Therefore, if (a1) holds, then for all .
Proof of (ii). Consider Consequently, if (a2) holds, .
Hence, is globally asymptotically stable for model (5) following the well-known theorem of Lyapunov stability.

4. Diffusion-Driven Instability: Turing Effect

In this section, we will investigate Turing instability and bifurcation for our model problem. We will also study pattern formation of the predator-prey solutions.

4.1. Turing Instability

Mathematically speaking, an equilibrium is Turing instability (diffusion-driven instability) means that it is an asymptotically stable equilibrium of model (6) but is unstable with respect to the solutions of diffusion model (5). In this subsection, we mainly focus on the emergency of the Turing instability of the positive equilibrium .

Now, the conditions for the positive equilibrium to be stable for the ODE are given by Hence, (the matrix ) has an eigenvalue with a positive real part, then it must be a real value and the other eigenvalue must be a negative real one. A necessary condition for the Turing instability of model (5) is Otherwise, for all eigenvalues of the operator since . For the Turing instability, we must have for some . And we notice that achieves its minimum at the critical value when

However, the inequality is necessary but not sufficient for the Turing instability in the bounded domain . The possible eigenvalues are discrete. In this case, has two positive roots and so if we can find some such that , then , and the positive equilibrium of model (5) is unstable.

Summarizing the previous analysis and calculations, we have the following results.

Theorem 6. Assume that the positive equilibrium exists. If the following conditions are true:(i) , that is, (ii) , that is, (iii) , that is, then the positive equilibrium of model (5) is Turing unstable if for some .

In Figure 1, we show the Turing bifurcation diagram for model (5) with parameters ,  ,  ,  , and   in parameters plane. The Turing bifurcation breaks spatial symmetry, leading to the formation of patterns that are stationary in time and oscillatory in space [33]. Below the Turing bifurcation curve, the solution of the model is unstable, and Turing instability emerges, that is, Turing patterns emerge. This domain is called the “Turing space.” We will focus on the Turing pattern formation in this domain.

4.2. Pattern Formation

In this section, we perform extensive numerical simulations of the spatially extended model (5) 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 ,  ,  ,  , and  .

The numerical integration of model (5) is performed by using a finite difference approximation for the spatial derivatives and an explicit Euler method for the time integration [34, 35] with a time stepsize of and the space stepsize of . The initial condition is always a small amplitude random perturbation around the positive equilibrium . After the initial period during which the perturbation spreads, either the model goes into a time dependent state, or to an essentially steady-state solution (time independent).

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 .

Figure 2 shows the evolution process of the holes pattern of prey for the parameters at ,   and iterations. In this case, one can see that for model (5), the random perturbations lead to the formation of stripes holes (cf. Figure 2(b)), and the later random perturbations make these stripes decay ending with the holes pattern (cf. Figure 2(d))—the prey are isolated zones with low population density.

Figure 3 shows the process of spatial pattern formation of prey for the parameters at ,  , and iterations. The random perturbations lead to the formation of stripe-holes patterns (cf. Figure 3(b)), and the later random perturbations make these holes decay, ending with a time-independent stripe pattern (cf. Figure 3(d))—the prey are interlaced stripes of high and low population densities.

Figure 4 shows the process of spot pattern formation of prey for at , , and iterations. There is a competition exhibited between stripes and spots. The pattern takes a long time to settle down, starting with a homogeneous state (cf. Figure 4(a)), and the random perturbations lead to the formation of stripes and spots (cf. Figure 4(b)), ending with spots only (cf. Figure 4(d))—the prey are isolated zones with high population density.

Ecologically speaking, spots pattern shows that the prey population is driven by predators to a high level in those regions, while holes pattern shows that the prey population is driven by predators to a very low level in those regions.

5. Conclusions and Remarks

In this paper, we have studied the dynamics of a reaction-diffusion Holling-Tanner prey-predator model where the predator population is subject to Allee effect under the zero-flux boundary conditions. The value of this study lies in threefolds. First, it investigates qualitative properties of solutions to this reaction-diffusion model. Second, it gives local and global stability of the positive equilibrium of the model. Third, it rigorously proves the Turing instability and illustrates three categories of Turing patterns close to the onset Turing bifurcation, which shows that the model dynamics exhibits complex pattern replication.

It is seen that if Allee effect constant is low, then the persistence of the model is guaranteed. It is interesting to notice that, from the result of Theorem 5, the condition for global stability of is independent of the diffusion coefficients (). So, it can be said that when the conditions on the parameters are satisfied, is stabilized under arbitrary spatially inhomogeneous perturbation. A very interesting observation can be made from the result of the numerical simulations. It indicates that the spatial model dynamics exhibits a diffusion-controlled formation growth not only to holes (cf. Figure 1) and stripes (cf. Figure 2) but also to spots replication (cf. Figure 3).

Comparing Figures 1 and 3, we can conclude that the model exhibits holes and spots Turing patterns due to Allee effect constant . It is believed that the observations made in this investigation related to Allee effect on predator population remind us of the importance of the Allee effect. Therefore, if the prey or the predator to be protected is subject to an Allee effect, the measures taken should take this into account.