# Pattern Formation in a Semi-Ratio-Dependent Predator-Prey System with Diffusion

Authors: Hunki Baek; Dong Ick Jung; Zhi-wen Wang
Journal: Discrete Dynamics in Nature and Society
 (2013)

---

## Abstract

We investigate spatiotemporal dynamics of a semi-ratio-dependent predator-prey system with reaction-diffusion and zero-flux boundary. We obtain the conditions for Hopf, Turing, and wave bifurcations of the system in a spatial domain by making use of the linear stability analysis and the bifurcation analysis. In addition, for an initial condition which is a small amplitude random perturbation around the steady state, we classify spatial pattern formations of the system by using numerical simulations. The results of numerical simulations unveil that there are various spatiotemporal patterns including typical Turing patterns such as spotted, spot-stripelike mixtures and stripelike patterns thanks to the Turing instability, that an oscillatory wave pattern can be emerged due to the Hopf and wave instability, and that cooperations of Turing and Hopf instabilities can cause occurrence of spiral patterns instead of typical Turing patterns. Finally, we discuss spatiotemporal dynamics of the system for several different asymmetric initial conditions via numerical simulations.


---

## Body

## 1. Introduction

In recent years, pattern formations in nonlinear complex systems have been one of the central problems of the natural, social, technological sciences and ecological systems [1–7]. Particularly, many researchers have studied the prey-predator system with reaction-diffusion. For example, in [8], Garvie and Trenchea presented the analysis of  reaction-diffusion predator-prey systems with the Holling type II functional response and provided an L∞ a priori estimate. Wang et al. in [9] investigated the spatial pattern formation of a predator-prey system with prey-dependent functional response of Ivlev type and reaction-diffusion. Also Zhang et al. in [10] studied a linear stability and bifurcation analysis including the Hopf and Turing bifurcation for a spatial Holling-type IV predator-prey model. In addition, Camara and Aziz-Alaoui, the authors of [11], considered a predator-prey system with a modified Leslie-Gower functional response modeled by a reaction-diffusion equation and derived the conditions for Hopf and Turing bifurcation in the spatial domain.In this context, in this paper, we will focus on the following a semi-ratio-dependent predator-prey system with reaction-diffusion:(1)∂N∂T=rN(1-NK)-aNN+bPP+D1∇2NinΩ,∂P∂T=eP(1-hPN)+D2∇2PinΩ,
where N,  P stand for the population densities of prey and predator with time T, respectively. The parameters r,  K,  a, and b are all positive constants standing for the prey intrinsic growth rate, the carrying capacity of prey, the capturing rate, and the half-saturation constant, respectively. The predator grows logically with the growth rate e and the carrying capacity N/h. In fact, the parameter h is a measure of the food quality that the prey provides for conversion into the predator. And D1 and D2 are diffusion coefficients. Here ∇2=∂2/∂x2+∂2/∂y2 is the usual Laplacian in two-dimensional and x and y stand for the space.In fact, the authors in [12–14] have explored a similar model to system (1) without reaction-diffusion, which is a kind of a temporal predator-prey system, and figured out dynamic properties of the temporal system such as stabilities of equilibria and permanence of the system. However, as mentioned above, it is important to investigate spatiotemporal dynamical behaviors of a diffusive predator-prey system like system (1). Thus, in this paper, we will concentrate our concerns on studying the bifurcation analysis and the spatiotemporal pattern formation analysis of system (1).Throughout this paper, we assume that no external input is imposed from the outside. Hence the boundary conditions are taken as∂N/∂n=∂P/∂n=0 on ∂Ω, Ω is the spatial domain in ℝ2, and n is the outward unit normal vector of the boundary ∂Ω. The main reason for choosing such boundary conditions is that we are interested in the self-organization of pattern. Also system (3) needs to be analyzed with the initial populations N(x,y,0)>0 and P(x,y,0)>0.In order to minimize the number of parameters in system (1), we set
(2)N=KN-,P=rKaP-,T=T-r.
For simplicity we will omit the bar notation in N-, P-, T- in the rest of the paper. Thus we can have the following system containing dimensionless quantities:
(3)∂N∂T=N(1-N)-NPN+αP+d1∇2N,∂P∂T=δP(1-βPN)+d2∇2P,
where rb/a=α,  e/r=δ,  rh/a=β,  D1/r=d1,  and D2/r=d2.For convenience, we setN(1-N)-NP/(N+αP)≡F(N,P) and δP(1-β(P/N))≡G(N,P).The main object of this paper is to look into the spatial dynamic behaviors of system (3). For this, in Section 2, we investigate bifurcation phenomena, Hopf, Turing, and wave bifurcations. In particular, we find out sufficient conditions for the Hopf bifurcation, the Turing instability, and the wave bifurcation of system (3). In Section 3, we give numerical simulations of system (3) with  random perturbation initial conditions of the stationary solution of the spatially homogeneous system to investigate various spatiotemporal pattern formations such as spotted, stripelike, and spiral patterns. Finally, we discuss dynamical behaviors of system (3) with several different initial conditions in Section 4.


## 2. Bifurcation Analysis

In order to investigate pattern formations of system (3), first we must think over the nonspatial system of system (3). In fact, the nonspatial system for system (3) has three equilibria, which correspond to spatially homogeneous equilibria of system (3), in the positive quadrant as follows:(i)
(0,0) (total extinct): a saddle point,(ii)
(1,0) (extinct of the predator): a saddle point,(iii)
a nontrivial positive stationary state(N*,P*) (coexistence of prey and predator) if α+β>1, where
(4)N*=1-1α+β,P*=1βN*=1β(1-1α+β).From the biological point of view, in this paper, we mainly focus on the dynamics of nontrivial stationary state (N*,P*).Now, to perform a linear stability analysis for the nontrivial stationary state(N*,P*), we must linearize the dynamical system (3) around the spatially homogeneous fixed point (N*,P*) for small space- and time-dependent fluctuations and expand them in Fourier space. For this, let
(5)N(x→,t)~N*eλteik→·x→,P(x→,t)~P*eλteik→·x→,
where x→=(x,y) and k→·k→=k2; k and λ are the wavenumber vector and frequency, respectively. Then we can obtain the corresponding characteristic equation as follows:
(6)|Ak-λI|=0,
where Ak=A-k2D, I is the identity matrix and D=diag(d1,d2) is the diffusion matrix and A is given by
(7)A=(∂F∂N∂F∂P∂G∂N∂G∂P)(N*,P*)=(1-2N*-αP*2(N*+αP*)2-N*2(N*+αP*)2δβP*2N*2δ-2δβP*N*)≡(FNFPGNGP).
Now, (6) can be solved, yielding the so-called characteristic polynomial:
(8)λ2-tr(Ak)λ+det(Ak)=0,
where
(9)tr(Ak)=FN+GP-k2(d1+d2)=-(d1+d2)k2-1-δ+1α+β+β(α+β)2,det(Ak)=FNGP-FPGN-k2(FNd2+GPd1)+k4d1d2=d1d2k4+(d2+d1δ-d2α+β-d2β(α+β)2)k2+δ-δα+β.
The solutions of (8) yield the dispersion relation
(10)λk±=12(tr(Ak)±tr(Ak)2-4det(Ak)).
The reaction-diffusion systems have led to the characterization of three basic types of symmetry-breaking bifurcations—Hopf, Turing, and wave bifurcation, which are responsible for the emergence of spatiotemporal patterns [3, 8, 9, 15–30].


### 2.1. Hopf Bifurcation

Hopf bifurcation is an instability induced by the transformation of the stability of a focus. In fact, the space-independent Hopf bifurcation breaks temporal symmetry of a system and gives rise to oscillations that are uniform in space and periodic in time. Mathematically speaking, the Hopf bifurcation occurs when(11)Im(λk)≠0,Re(λk)=0atk=0.
Thus we can figure out that the critical value of Hopf bifurcation parameter δ is equal to
(12)δH=β(α+β)2+1α+β-1.
We assume that the value δH is greater than 0 for the existence of Hopf bifurcation. In other words, from now on, we will consider the parameters α and β which are satisfied with the following condition:
(13)α2+(2β-1)α+β(β-2)<0.
Thus the positive stationary solution (N*,P*) of system (3) is an unstable equilibrium solution if δ<δH and oscillates periodically with respect to time with the frequency
(14)ωH=Im(λk)=det(A)=δ-δα+β.
The corresponding wavelength is
(15)λH=2πωH=2πα+βδ(α+β-1).


### 2.2. Turing Instability

Turing instability (or called Turing bifurcation) is a phenomenon that causes certain reaction-diffusion system to lead to spontaneous stationary configuration. That is why Turing instability is often calleddiffusion-driven instability. The Turing instability is not dependent on the geometry of the system but only on the reaction rates and diffusion. It can occur only when the inhibitor (P) diffuses faster than the activator (N) [3, 8, 9, 20, 21, 28].In fact, the Turing instability sets in when at least one of the solutions of (8) crosses the imaginary axis. In other words, the spatially homogeneous steady state will become unstable due to heterogeneous perturbation when at least one solution of (8) is positive. For this reason, at least one out of the following two inequalities is violated to occur the Turing instability phenomenon:
(16)tr(Ak)=FN+GP-(d1+d2)k2<0,det(Ak)=d1d2k4-(d1GP+d2FN)k2+FNGP-FPGN>0.
Since perturbations of zero wavenumber are stable when diffusive instability sets in by definition (due to the stability for the nonspatial steady state), tr(A)=FN+GP<0 and det(A)=FNGP-FPGN>0 are satisfied. It is seen from these facts that the first condition in (16) always holds. Hence we have only one for the instability condition; that is, H(k2)≡det(Ak)<0. Elementary calculations yield that
(17)H(k2)=d1d2k4+(d1δ+d2-d2α+β-d2β(α+β)2)k2+δ(1-1α+β)<0.
Thus the minimum of H(k2) occurs at the critical wavenumber kT2, where
(18)kT2=d2FN+d1GP2d1d2=12d1d2(-d1δ-d2+d2α+β+d2β(α+β)2).
By substituting k2=kT2 into H(k2), we can get a sufficient condition for the Turing instability as follows:
(19)H(kT2)=Bδ2+12Cδ+D216B<0,
where
(20)B=-d14d2,C=1-1α+β+β(α+β)2,D=-1+1α+β+β(α+β)2.
It is from replacing the inequality in (19) by the equality that the critical value of bifurcation parameter δ for the Turing bifurcation equals
(21)δT=d2d1(β(α+β)2(1-1α+β)1-1α+β+β(α+β)2ZZZ+2β(α+β)2(1-1α+β)).
At the Turing threshold δT, the spatial symmetry of the system is broken and the patterns are stationary in time and oscillatory in space with the wavelength
(22)λT=2πkT.


### 2.3. Wave Bifurcation

The wave instability caused by the wave bifurcation plays an important part in pattern formations in many areas [7]. In fact, mathematically speaking, the wave bifurcation occurs when
(23)Im(λk)≠0,Re(λk)=0atk=kw≠0.
Thus, from elementary calculation, the critical value of wave bifurcation parameter δ can be obtained as
(24)δw=β2(α+β)2+12d12((d1+d2)Φ+2d1d2(1-1α+β)+βd22(α+β)2),
where
(25)Φ=β(α+β)2[(d12+d22)β(α+β)2β(α+β)2Z+2d1d2(2(1-1α+β)-β(α+β)2)].It is well known that, at the wave thresholdδw, both spatial and temporal symmetries are broken and the patterns are oscillatory in space and time with the wavelength
(26)λw=2πkw,
where
(27)kw2=1d1+d2(-1-δ+1α+β+β(α+β)2).If one takesα=0.4, β=5/6, and d1=1 in system (3), the bifurcation diagram, shown in Figure 1, can be obtained using the above linear stability analysis. Since Turing patterns occur when the inhibitor (P) diffuses faster than the activator (N), we take into account the range of the value d2 greater than d1=1. According to Figure 1, the Hopf bifurcation line and the Turing bifurcation line intersect at a point (4.0611,0.3587) and the Turing bifurcation line and the wave bifurcation line meet at a point (1,0.0932). In addition, these bifurcation lines split the parametric space into five distinct domains. In domain I, located above two bifurcation lines, the steady state is the only stable solution of system (3). Domain II is the region of the pure Turing instability and domain III is the region of Hopf and wave instabilities. Hopf and Turing instabilities take place in domain IV. On the other hand, all three instabilities occur in domain V. Furthermore, Figure 2 shows the dispersion relations corresponding to several values of two parameters d2,δ. The transitions of Hopf, Turing, and wave modes from stable to unstable for system (3) can be ascertained clearly via Figure 2. While Hopf is subcritical, Turing and wave bifurcations are supercritical.Figure 1
Bifurcation diagrams of system (3) with respect to d2 and δ when α=0.4, β=5/6, and d1=1. Hopf bifurcation line: δ=0.3587; Turing bifurcation line: δ=0.0932d2; wave bifurcation line: δ=(7/37)d2+(375/1369)d22-(1015/1369)375d22-232d2+375+(375/1369). Hopf-Turing bifurcation point: (4.0611,0.3587). Turing-wave bifurcation point: (1,0.0932).Figure 2
Dispersion relations: (I)(d2,δ)=(2,0.5); (II) (d2,δ)=(7,0.5); (III) (d2,δ)=(2,0.3); (IV1) (d2,δ)=(7,0.2); (IV2) (d2,δ)=(2,0.15); (V) (d2,δ)=(7,0.3). Solid line: the real parts of the eigenvalues; dashed line: imaginary parts of complex pair of eigenvalues.


## 2.1. Hopf Bifurcation

Hopf bifurcation is an instability induced by the transformation of the stability of a focus. In fact, the space-independent Hopf bifurcation breaks temporal symmetry of a system and gives rise to oscillations that are uniform in space and periodic in time. Mathematically speaking, the Hopf bifurcation occurs when(11)Im(λk)≠0,Re(λk)=0atk=0.
Thus we can figure out that the critical value of Hopf bifurcation parameter δ is equal to
(12)δH=β(α+β)2+1α+β-1.
We assume that the value δH is greater than 0 for the existence of Hopf bifurcation. In other words, from now on, we will consider the parameters α and β which are satisfied with the following condition:
(13)α2+(2β-1)α+β(β-2)<0.
Thus the positive stationary solution (N*,P*) of system (3) is an unstable equilibrium solution if δ<δH and oscillates periodically with respect to time with the frequency
(14)ωH=Im(λk)=det(A)=δ-δα+β.
The corresponding wavelength is
(15)λH=2πωH=2πα+βδ(α+β-1).


## 2.2. Turing Instability

Turing instability (or called Turing bifurcation) is a phenomenon that causes certain reaction-diffusion system to lead to spontaneous stationary configuration. That is why Turing instability is often calleddiffusion-driven instability. The Turing instability is not dependent on the geometry of the system but only on the reaction rates and diffusion. It can occur only when the inhibitor (P) diffuses faster than the activator (N) [3, 8, 9, 20, 21, 28].In fact, the Turing instability sets in when at least one of the solutions of (8) crosses the imaginary axis. In other words, the spatially homogeneous steady state will become unstable due to heterogeneous perturbation when at least one solution of (8) is positive. For this reason, at least one out of the following two inequalities is violated to occur the Turing instability phenomenon:
(16)tr(Ak)=FN+GP-(d1+d2)k2<0,det(Ak)=d1d2k4-(d1GP+d2FN)k2+FNGP-FPGN>0.
Since perturbations of zero wavenumber are stable when diffusive instability sets in by definition (due to the stability for the nonspatial steady state), tr(A)=FN+GP<0 and det(A)=FNGP-FPGN>0 are satisfied. It is seen from these facts that the first condition in (16) always holds. Hence we have only one for the instability condition; that is, H(k2)≡det(Ak)<0. Elementary calculations yield that
(17)H(k2)=d1d2k4+(d1δ+d2-d2α+β-d2β(α+β)2)k2+δ(1-1α+β)<0.
Thus the minimum of H(k2) occurs at the critical wavenumber kT2, where
(18)kT2=d2FN+d1GP2d1d2=12d1d2(-d1δ-d2+d2α+β+d2β(α+β)2).
By substituting k2=kT2 into H(k2), we can get a sufficient condition for the Turing instability as follows:
(19)H(kT2)=Bδ2+12Cδ+D216B<0,
where
(20)B=-d14d2,C=1-1α+β+β(α+β)2,D=-1+1α+β+β(α+β)2.
It is from replacing the inequality in (19) by the equality that the critical value of bifurcation parameter δ for the Turing bifurcation equals
(21)δT=d2d1(β(α+β)2(1-1α+β)1-1α+β+β(α+β)2ZZZ+2β(α+β)2(1-1α+β)).
At the Turing threshold δT, the spatial symmetry of the system is broken and the patterns are stationary in time and oscillatory in space with the wavelength
(22)λT=2πkT.


## 2.3. Wave Bifurcation

The wave instability caused by the wave bifurcation plays an important part in pattern formations in many areas [7]. In fact, mathematically speaking, the wave bifurcation occurs when
(23)Im(λk)≠0,Re(λk)=0atk=kw≠0.
Thus, from elementary calculation, the critical value of wave bifurcation parameter δ can be obtained as
(24)δw=β2(α+β)2+12d12((d1+d2)Φ+2d1d2(1-1α+β)+βd22(α+β)2),
where
(25)Φ=β(α+β)2[(d12+d22)β(α+β)2β(α+β)2Z+2d1d2(2(1-1α+β)-β(α+β)2)].It is well known that, at the wave thresholdδw, both spatial and temporal symmetries are broken and the patterns are oscillatory in space and time with the wavelength
(26)λw=2πkw,
where
(27)kw2=1d1+d2(-1-δ+1α+β+β(α+β)2).If one takesα=0.4, β=5/6, and d1=1 in system (3), the bifurcation diagram, shown in Figure 1, can be obtained using the above linear stability analysis. Since Turing patterns occur when the inhibitor (P) diffuses faster than the activator (N), we take into account the range of the value d2 greater than d1=1. According to Figure 1, the Hopf bifurcation line and the Turing bifurcation line intersect at a point (4.0611,0.3587) and the Turing bifurcation line and the wave bifurcation line meet at a point (1,0.0932). In addition, these bifurcation lines split the parametric space into five distinct domains. In domain I, located above two bifurcation lines, the steady state is the only stable solution of system (3). Domain II is the region of the pure Turing instability and domain III is the region of Hopf and wave instabilities. Hopf and Turing instabilities take place in domain IV. On the other hand, all three instabilities occur in domain V. Furthermore, Figure 2 shows the dispersion relations corresponding to several values of two parameters d2,δ. The transitions of Hopf, Turing, and wave modes from stable to unstable for system (3) can be ascertained clearly via Figure 2. While Hopf is subcritical, Turing and wave bifurcations are supercritical.Figure 1
Bifurcation diagrams of system (3) with respect to d2 and δ when α=0.4, β=5/6, and d1=1. Hopf bifurcation line: δ=0.3587; Turing bifurcation line: δ=0.0932d2; wave bifurcation line: δ=(7/37)d2+(375/1369)d22-(1015/1369)375d22-232d2+375+(375/1369). Hopf-Turing bifurcation point: (4.0611,0.3587). Turing-wave bifurcation point: (1,0.0932).Figure 2
Dispersion relations: (I)(d2,δ)=(2,0.5); (II) (d2,δ)=(7,0.5); (III) (d2,δ)=(2,0.3); (IV1) (d2,δ)=(7,0.2); (IV2) (d2,δ)=(2,0.15); (V) (d2,δ)=(7,0.3). Solid line: the real parts of the eigenvalues; dashed line: imaginary parts of complex pair of eigenvalues.


## 3. Spatiotemporal Pattern Analysis via Numerical Simulations

In this section, we will investigate spatiotemporal pattern formations of the spatially extended system (3) in two-dimensional space via numerical examples. In fact, the nonuniform stationary states of system (3) that corresponds to spatial patterns cannot be found analytically. In other words, the analytical methods are not sufficient to fully understand the system, which is a reason why we need to use computer simulations. For this, first, take α=0.4, β=5/6, and d1=1 in system (3). Then the spatially homogeneous equilibrium (N*,P*)=(0.1892,0.2270) is always the only nontrivial positive stationary state for the nonspatial system of system (3) regardless of the values d2 and δ.In order to solve partial differential equations numerically, one has to discretize the space and time of the given problem. For this reason, the discrete domain for the continuous domainΩ in system (3) is assumed to be the 200×200 lattice sites and the spacing between the lattice points is assumed to be Δx=Δy=1. The time evaluation is also discrete; that is, the time goes in steps of Δt. And all our numerical simulations employ the zero-flux boundary condition.System (3) is solved numerically by using a finite difference scheme for the spatial derivatives and an explicit Euler method for the time integration with the time step Δt=0.01. These numerical schemes for the diffusion equation (3) give stable solutions as long as the equalities d1Δt/(Δx)2,d2Δt/(Δy)2<1/2 are satisfied (see [31]).It is well known that the spatiotemporal dynamics of a diffusion-reaction system depends on the choice of initial conditions [3, 9, 32, 33]. In this section, an initial condition is taken as a small amplitude random perturbation around the steady state (N*,P*) since it is very natural from the biological point of view. We stop the simulation when the numerical solutions either reach a stationary state or show oscillatory behaviors.In the numerical simulations, different types of spatiotemporal dynamics are observed and we have found that the distributions of predator and prey are always of the same type. Consequently, we can restrict our analysis of pattern formations to one distribution. In this section, we show the distribution of predator, for instance.Now, we will classify spatiotemporal patter formations of system (3) for each bifurcation domain in Figure 1 by taking advantage of numerical simulations.As mentioned in Section2, the domain I in Figure 1 is the region where the steady state is the only stable solution of system (3). Thus we will omit to mention dynamical behaviors of system (3) for this domain.First, we will investigate dynamical behaviors of system (3) with a random small perturbation initial condition of the stationary solution (N*,P*) by taking into account the parameters in domain II. As mentioned in Section 2, the domain II is the region of the pure Turing instability. If we take the parameter values d2=7 and δ=0.6, we can see from Figure 3 that the random initial distribution (see Figure 3(a)) leads to the formation of a regular macroscopic spotted pattern which prevails over the whole domain at last, and after that the dynamics of the system does not change their behaviors anymore.Snapshots of contour pictures of the time evolution of predator in system (3) when d2=7 and δH<δ=0.6<δT: (a) 0 iteration, (b) 30000 iterations, (c) 100000 iterations, and (d) 200000 iterations.
(a)
(b)
(c)
(d)From now on, through this section, we will not display the snapshots of the initial pattern as Figure3(a) since all following figures are obtained with a random small perturbation initial condition of the stationary solution as the previous case.If we take the parameter valuesd2=7 and δ=0.535, we can observe from Figure 4 that the steady state of spotted pattern and the stripelike pattern coexist. In addition, for the parameters d2=7 and δ=0.4, Figure 5 shows that stripelike spatial patterns are prevalent in the whole domain eventually. In fact, we can see that the stripelike patterns grow steadily with time and finally they overwhelm the space domain (see Figure 5(c)) after the spotted pattern and the stripelike pattern coexist (see Figures 5(a) and 5(b)).Snapshots of contour pictures of the time evolution of predator in system (3) when d2=7 and δH<δ=0.535<δT: (a) 30000 iterations, (b) 100000 iterations, and (c) 200000 iterations.
(a)
(b)
(c)Snapshots of contour pictures of the time evolution of predator in system (3) when d2=7 and δH<δ=0.4<δT: (a) 20000 iterations, (b) 40000 iterations, and (c) 200000 iterations.
(a)
(b)
(c)Next, in order to study dynamical behaviors of system (3) in domain V, set the parameters d2=7 and δ=0.34. We can know from Figure 1 that, in this case, all three instabilities occur at the same time. In this context, we display Figure 6 to show snapshots of predator spatial patterns at 50000, 100000, and 200000 iterations. Since Turing instabilities take place in domain V, we can see the stripelike spatial patterns as shown in Figure 5. However, an oscillatory pattern in time can be observed even though we cannot observe this phenomenon via snapshots of Figure 6. To make it clear, we exhibit the local phase portrait of system (3) for a fixed point (0.1892,0.2270) in Figure 7. As one sees, thanks to the Hopf instability, a stable cycle can be observed. Moreover, we can ascertain numerically that the frequency of the periodic oscillations in time ω=0.2685 and the corresponding wavelength λ=23.4. Also, it follows from (14) and (15) that the frequency of the periodic oscillations in time is ωH=0.2536 and the corresponding wavelength λH=24.7738.Snapshots of contour pictures of the time evolution of predator in system (3) when d2=7 and δw<δ=0.34<δH: (a) 50000 iterations, (b) 100000 iterations, and (c) 200000 iterations.
(a)
(b)
(c)Figure 7
Dynamical behavior of system (3) for a fixed point (0.1892,0.2270) when d2=7 and δ=0.34: (a) time series of prey, (b) time series of predator, and (c) the stable limit cycle of system (3).Thus Figures3, 4, 5, and 6 show typical Turing patterns of predator in system (3) for the fixed value d2=7 when the parameter δ varies. In fact, as the value δ decreases, the sequence spotted → spot-stripelike mixtures → stripelike pattern is observed.Now, consider the parametersd2=2 and δ=0.34 in domain III. It is from the linear stability analysis in Section 2 that Hopf and wave instabilities occur. As shown in Figure 8, a special pattern may not be formed such as spotted or stripelike patterns since Turing instabilities do not occur in this domain. However, in this case, an oscillatory wave pattern formation is observed. In order to deal with this phenomenon in more detail, we display the local phase portrait of system (3) for a fixed point (0.1892,0.2270) in Figure 9(a), which shows a stable limit cycle, and the population density of predator in Figure 9(b), which is periodic in time and the space-time plots in Figure 9(c), obtained by piling up the lines y=100 in each snapshot as time progresses. Thus, we can clearly figure out that Hopf and wave instabilities play an important role in forming a wave and periodic oscillatory pattern in space and time as mentioned in Section 2. In fact, numerical calculations yield that the frequency of the periodic oscillations in time ω=0.2474 and the corresponding wavelength λ=25.4. However, the theoretical values of the frequency and the corresponding wavelength are ωH=0.2536 and λH=24.7738, respectively. However, the theoretical values kw=0.0789 and λw=79.6768 can be obtained from (26) and (27). From these facts, we can surmise that, in this case, Hopf instabilities have an stronger effect on pattern formation in terms of time than wave instabilities.Snapshots of contour pictures of the time evolution of predator in system (3) when d2=2 and δT<δ=0.34<δH: (a) 1000 iterations, (b) 50000 iterations, and (c) 200000 iterations.
(a)
(b)
(c)Dynamical behavior of system (3) when d2=2 and δ=0.34: (a) phase portrait of system (3) for a fixed point (0.1892,0.2270), (b) time series of predator, and (c) space-time plots.
(a)
(b)
(c)On the contrary, domain IV is the region of Hopf and Turing instabilities according to Figure1. To discuss dynamical behaviors of system (3) in this domain, we think about two parameter sets as d2=2, δ=0.15 and d2=7, δ=0.2. As one sees, Figures 10 and 11 show one of the typical spatial patterns called a spiral pattern instead of patterns shown in Turing domains II and V even if Turing instabilities occur in this domain. In fact, in many recent researches [8, 9, 28, 33], spiral patterns have been reported. However, most of them used specific initial conditions which are not the same as the conditions used in this section to generate a spiral pattern. Thus a diffusive system starting with a small amplitude random perturbation around the steady state could have spiral patterns thanks to the cooperations of Turing and Hopf instabilities.Snapshots of contour pictures of the time evolution of predator in system (3) when d2=2 and δ=0.15<δw: (a) 5000 iterations, (b) 150000 iterations, and (c) 200000 iterations.
(a)
(b)
(c)Snapshots of contour pictures of the time evolution of predator in system (3) when d2=7 and δ=0.2<δw: (a) 5000 iterations, (b) 30000 iterations, and (c) 261000 iterations.
(a)
(b)
(c)


## 4. Conclusion and Discussion

Various spatiotemporal patterns generated by a semi-ratio-dependent predator-prey system with reaction-diffusion are studied theoretically and numerically. The spatial parametric domain is divided into five regions by three bifurcation lines, Hopf, Turing, and wave bifurcation lines, shown in Figure1. We show that in domains II and V, the sequence spotted → spot-stripelike mixtures → stripelike pattern is emerged as the growth rate of predator decreases since the Turing instability occurs in domain III, an oscillatory wave pattern emerges due to the Hopf and wave instabilities and in domain IV a typical spiral pattern is investigated because of Hopf and Turing instabilities. Moreover, we calculate numerically the values of the frequency and wavelength of the periodic oscillations in time and compare these values to the theoretical values obtained from (14), (15), (26), and (27).The authors in [3, 9, 32, 33] have investigated the patterns arising from special initial conditions. They indicated that spatiotemporal dynamics of a reaction-diffusion system depend on the choice of initial conditions. For discussing the evolution process of the spatial pattern formation of system (3) with some specific initial conditions, we employ two specific initial conditions which were considered in [3, 9, 32, 33].Firstly, we employ an initial condition that predators are introduced in a small, localized region of the circle domain while the population of prey is fixed:(28)N(x,y,0)=1.0,(29)P(x,y,0)={0.2,(x-100)2+(y-100)2<100,0,otherwise.
Secondly, we take an initial condition that describes prey patch placed into a domain with a constant-gradient predator distribution as follows:
(30)N(x,y,0)=N*-ϵ1(x-180)(x-720)-ϵ2(y-90)(y-210),P(x,y,0)=P*-ϵ3(x-450)-ϵ4(y-135),
where ϵ1=2×10-7, ϵ2=6×10-7, ϵ3=3×10-5, and ϵ4=6×10-5.Figures12 and 13 illustrate snapshots of contour pictures of the time evolution of the predator in system (3) with the initial conditions (29) and (30), respectively, when d2=7 and δ=0.6, for 0,100000, and 1000000 iterations.Snapshots of contour pictures of the time evolution of predator in system (3) with the initial condition (29) when d2=7 and δ=0.6: (a) 0 iteration, (b) 100000 iterations, and (c) 1000000 iterations.
(a)
(b)
(c)Snapshots of contour pictures of the time evolution of predator in system (3) with the initial condition (30) when d2=7 and δ=0.6: (a) 0 iteration, (b) 100000 iterations, and (c) 1000000 iterations.
(a)
(b)
(c)These figures beg the question that the spotted pattern could be a globally stable pattern even though we do not have any theoretical evidences. Similarly, we can figure out numerically that system (3) with the parameters d2=7 and δ=0.3, which lies on domain V, has the same phenomena even if we do not exhibit their snapshots. In this case, the system has a stripelike pattern as shown in Figure 5.



---

*Source: 657286-2013-04-04.xml*