We present a theoretical analysis of processes of pattern formation that involves organisms distribution and their interaction of spatially distributed population with cross-diffusion in a Harrison-type predator-prey model. We analyze the global behaviour of the model by establishing a Lyapunov function. We carry out the analytical study in detail and find out the certain conditions for Turing’s instability induced by cross-diffusion. And the numerical results reveal that, on increasing the value of the half capturing saturation constant, the sequences “spots → spot-stripe mixtures → stripes → hole-stripe mixtures → holes” are observed. The results show that the model dynamics exhibits complex pattern replication controlled by the cross-diffusion.

1. Introduction

Understanding of spatial and temporal behaviors of interacting species in ecological systems is one of the central scientific problems in population ecology [115], since the pioneering work of Turing [16]. Throughout the history of theoretical ecology, reaction-diffusion equations have been intensively used to describe spatiotemporal dynamics [15, 17].

In recent years, the effect of cross-diffusion in reaction-diffusion systems has received much attention by both ecologists and mathematicians, for example, see [1825] and the references therein. Kerner [18] was the first to examine that cross-diffusion can induce pattern forming instability in an ecological situation. Cross-diffusion expresses the population fluxes of one species due to the presence of the other species. And Gurtin [19] developed some mathematical models for population dynamics with the inclusion of cross-diffusion as well as self-diffusion and showed that the effect of cross-diffusion may give rise to the segregation of two species.

In this paper, we are attempting to study the effect of cross-diffusion in a predator-prey model with Harrison-type functional response [26]. The model can be written as where and represent population density of prey and predator at time , respectively. is the intrinsic growth rate of prey, is the prey carrying capacity, is the capture rate, the half capturing saturation constant, is the death rate of predator, and is conversion rate. and are the self-diffusion coefficients of and , respectively, and are the cross-diffusion coefficients of and , respectively. We always assume that , and . The value of the cross-diffusion coefficient may be positive, negative, or zero. Positive cross-diffusion coefficient denotes, that one species tends to move in the direction of lower concentration of another species, while negative cross-diffusion expresses the population fluxes of one species in the direction of higher concentration of the other species [27]. is the usual Laplacian operator in 2-dimensional space. is a bounded domain with smooth boundary . The initial data and are continuous functions on .

We make a change of variables: For the sake of convenience, we still use variables , instead of , . Thus, considering zero-flux boundary conditions, model (1) is converted into where the new parameters are is the diffusion matrix, , , and . is the outward unit normal vector on and the zero-flux boundary conditions mean that model (3) is self-contained and has no population flux across the boundary [28, 29].

In particular, when , that is, the cross-diffusion coefficients are equal 0, we can obtain the following model: We call model (5) as self-diffusion model, while we call model (3) cross-diffusion model.

The corresponding kinetic equation to models (3) and (5) is:

In recent years, there has been considerable interest to investigate the stability behavior of a predator-prey system by taking into account the effect of self- as well as cross-diffusion [3, 6, 810, 1215]. But in the studies on the spatiotemporal dynamics of predator-prey system with functional response, little attention has been paid to study on the effect of cross-diffusion.

Mathematically speaking, an equilibrium in Turing’s instability (diffusion-driven instability) means that it is an asymptotically stable equilibrium of model (6) but is unstable with respect to the solutions of reaction-diffusion model (3) or (5). Especially, if is also stable with respect to the solutions of the self-diffusion model (3), that is, in the cross-diffusion model (5), then there is nonexistence of Turing’s instability in this situation.

And there comes a question: if there is nonexistence of Turing’s instability in the case of self-diffusion (i.e., ), does model (5) exhibit Turing’s instability induced by cross-diffusion?

The main purpose of this paper is to focus on the effect of cross-diffusion on the spatiotemporal dynamics of the reaction-diffusion predation model. The paper is organized as follows. In Section 2, we give some properties of the solutions of the model. In Section 3, we give the linearized stability analysis to show (i.e., no cross-diffusion), deduce the conditions of Turing’s instability induced by cross-diffusion, and illustrate the different Turing patterns by using the numerical simulations. Finally, in Section 4, some conclusions and discussions are given.

2. Dynamics Analysis

In this section, we present some preliminary results, including dissipativeness, boundedness, permanence of the solutions, and the equilibria stability analysis of the models.

2.1. Dissipativeness

Theorem 1. For any solution of model (6), Hence, model (6) is dissipative.

Proof. From the first equation of model (6), it can be easily shown that If , from the second equation of model (6), one has: A standard comparison argument shows that If , we have the following differential inequality: and the same argument above yields In either case, the second inequality of (7) holds.

2.2. Boundedness

Theorem 2. All the solutions of model (6) which initiate in are uniformly bounded within the region , where

Proof. Let us define the function: Calculating the time derivative of along the trajectories of model (6), we get Then, Using the theory of differential inequality, for all , we have Hence, we have Hence, all the solutions of model (6) that initiate in are confined in the region .

2.3. Permanence

Theorem 3. If and , then model (6) has the permanence property.

Proof. By the first equation of model (6), we have Since , by the famous comparison theorem, we have Hence, for large , .
As a result, for large , satisfies Since , by the famous comparison theorem, we can get The proof is complete.

2.4. Stability Analysis of the Equilibria

The nonspatial model (6) has three equilibria, which correspond to spatially homogeneous equilibria of model (3) and model (5), in the positive quadrant:(i) (total extinct) is a saddle point;(ii) (extinct of the predator, or prey only) is a saddle when , or stable node when ;(iii) (coexistence of prey and predator), where , and satisfies It is easy to verify that it has a unique positive equilibrium if .

The Jacobian matrix for the positive equilibrium is given by Obviously,

Therefore, we can obtain the following.

Theorem 4. Assume that the positive equilibrium exists, then is locally stable for model (6).

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

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

Proof. We adopt the Lyapunov function: where , .
Then, where
By some computational analysis, we obtain Considering the zero-flux boundary conditions, we have
and in the curly brackets can be expressed in the form and , respectively, where is negative definite if the symmetric matrices and are positive. It can be easily shown that the symmetric matrix is positive definite if
The symmetric matrix is positive definite if Since due to (a1), is true. It is easy to verify that .
Since then Hence, is strictly increasing in , with respect to , and Consequently, if (a2) holds, . As a result, .
Hence, is a Lyapunov function and the positive equilibrium of model (3) is globally asymptotically stable. This completes the proof.

Remark 6. When , Theorem 5 is true, too. That is, the positive equilibrium of the self-diffusion model (5) is globally asymptotically stable.

3. Turing’s Instability and Pattern Formation

3.1. Nonexistence of Turing’s Instability in the Self-Diffusion Model (5)

And in the presence of diffusion, we will introduce small perturbations , , where . To study the effect of self-diffusion on model (5), we consider the linearized form of system about as follows: where are defined as (24).

Following Malchow et al. [11], we can know any solution of model (38) can be expanded into a Fourier series so that where , and , . and , are the corresponding wavenumbers.

Substituting and into (38), we obtain where .

A general solution of (40) has the form , where the constants and are determined by the initial conditions and the exponents are the eigenvalues of the following matrix:

Correspondingly, , are the solution of the following characteristic equation: where

From (24), one can easily obtain Then, we can conclude that the equilibrium is also stable for self-diffusion model (5). That is to say, there is nonexistence of Turing’s instability in model (5).

3.2. Turing’s Instability in the Cross-Diffusion Model (3)

The linearized form of the cross-diffusion model (3) about is as follows:

The characteristic equation of the linearized model (3) is: where

Diffusive instability occurs when at least one of the following conditions is violated [2]:

It is evident that the condition is not violated when the requirement is met because we assume and . Hence, only violation of the condition will give rise to diffusion instability, that is, Turing’s instability. Then the condition for diffusive instability is given by otherwise for all since and .

For Turing’s instability, we must have for some . And we notice that achieves its minimum: at the critical value when As a consequence, if and hold, then is an unstable equilibrium with respect to model (3). In this case, has two positive roots and which satisfy where . Therefore, if we can find some such that , then .

Summarizing the above calculation, we obtain the following.

Theorem 7. Assume that the positive equilibrium exists. If the following conditions are true:(i), that is, (ii), that is,

then the positive equilibrium of model (3) is Turing unstable if    for some .

3.3. Pattern Formation

In this section, we perform extensive numerical simulations of the spatially extended model (3) 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 (3) is performed by using a finite difference approximation for the spatial derivatives and an explicit Euler method for the time integration [30, 31] with a time stepsize of and the space stepsize . The initial condition is always a small amplitude random perturbation around the positive equilibrium . After the initial period during which the perturbation spread, 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 1 shows five typical Turing’s patterns of prey in model (3) arising from random initial conditions for several values of the control parameter .

In Figure 1(a)-pattern, , consists of red (maximum density of ) hexagons on a blue (minimum density of ) background, that is, isolated zones with high population densities. In this paper, we call this pattern as “spots.”

In Figure 1(b), when increasing to , a few of stripes emerge, and the remainder of the spots pattern remains time independent. Pattern (b) is called -hexagon-stripe mixtures pattern.

While increasing to , model dynamics exhibits a transition from stripes-spots growth to stripes replication, that is, spots decay and the stripes pattern emerges (cf. Figure 1(c)).

In Figure 1(d), , on increasing of , a few of blue hexagons (i.e., holes, named by Von Hardenberg et al. [32], associated with low population densities) fill in the stripes, that is, the stripes-holes pattern emerges. Pattern (d) is called -hexagon-stripe mixtures pattern.

When increasing to , model dynamics exhibits a transition from stripe-holes growth to spots replication, that is, stripes decay and the holes pattern (-pattern) emerges (cf. Figure 1(e)).

From Figure 1, one can see that, on increasing the control parameter , the sequences “spots → spot-stripe mixtures → stripes → hole-stripe mixtures → holes” are observed. Ecologically speaking, spots pattern shows that the prey population are driven by predators to a high level in those regions, while holes pattern shows that the prey population are driven by predators to a very low level in those regions. The final result is the formation of patches of high prey density surrounded by areas of low prey densities [3].

4. Conclusions and Remarks

In this paper, we study the spatiotemporal dynamics of a Harrison predator-prey model with self- and cross-diffusions under the zero-flux boundary conditions. The value of this study lies in twofold. First, it gives the global stability of the positive equilibrium of the model by establishing a Lyapunov function. Second, it rigorously proves that the Turing instability can be induced by cross-diffusion, which shows that the model dynamics exhibits complex pattern replication controlled by the cross-diffusion.

The most important observation in this paper is that the cross-diffusion terms are necessary for the emergence of Turing’s instability and pattern formation in the model. More precisely, with the help of the numerical simulations, the sequences “spots → spot-stripe mixtures → stripes → hole-stripe mixtures → holes” can be observed.

On the other hand, population dynamics in the real world is inevitably affected by environmental noise which is an important component in an ecosystem. The deterministic models, such as model (3) or (5), assume that parameters in the systems are all deterministic irrespective environmental fluctuations. It is well known that the fact that due to environmental noise, the birth rate, carrying capacity, competition coefficient, and other parameters involved in the system exhibit random fluctuation to a greater or lesser extent [33]. We think that there may be exist other noise-controlled self-replicating patterns in models (3) and (5). This is desirable in future studies.

It is believed that our results related to cross-diffusion in predator-prey interactions model would certainly be of some help to theoretical mathematicians and ecologists who are engaged in performing experimental work.