Phosphate deposits in south-western Iran are part of the South Tethyan Phosphogenic Province, a huge carbonate-dominated strata that extends to the Middle East. The Tethyan phosphorites of Iran are dated Eocene-Oligocene (Pabdeh Formation) and categorized as low-grade ore deposits on a global scale. Depositional conditions of the facies indicate that the Pabdeh Formation was deposited on a carbonate ramp setting as a distally steepened ramp. Under such an environment, turbidity currents transported phosphate particles from the back-shoal setting to the deeper middle and outer ramp of the ocean where they were suspended and deposited as shell-lag and phosphate lamination. Microfacies studies demonstrate that all the phosphatic ooids and phosphatized foraminifera, fish scales, bones and phosphatic intraclasts reworked from shallow parts of the Tethyan Ocean to deeper parts with the help of turbidity currents. Analysis and interpretation of the data reveal positive correlation between. REE+Y and P2O5 in all studied sections which attests to their strong coherence as a geochemical group. The shale normalized REE patterns of Mondun phosphorites are characterized by negative Ce anomalies. This anomaly indicate that the depositional environment was oxic and highly reworked, bioturbated with higher energy realm during phosphate deposition, conversely Nill section with Ce enrichment reflect conditions of relatively deeper water sedimentation. These geochemical findings are in accordance with microfacies studies which indicate shallow and high energy condition for Mondun section with negative cerium anomalies and a deep ramp setting for Nill and Siah sections which denote a positive cerium anomalies in REE patterns.
1.
Introduction
The relationship between predator and prey is an important research content in ecosystem, and many scholars have studied this interaction by differential equation models [1,2,3,4,5,6]. The direct relationship between predator and prey is the consumption of prey. One of the indirect relationship between predator and prey is the fear of predator. When the prey gets the predator's signal (chemical/vocal), they will increase their vigilance time and reduce their foraging [7], such as mule deer v.s. mountain lions [8], elk v.s. wolves [9]. Consider the fear effect, Panday et al. [10] proposed the following model
All parameters are positive. The biological interpretation of parameters is given in Table 1. They incorporated fear effect by modifying the prey intrinsic growth rate R0 as R01+Kv. By the scaling
model (1.1) is changed to (drop the bars)
Panday et al. [10] introduced time delay of perceiving predator signals to model (1.3), and mainly studied the boundedness, persistence, local and global behavior of the delayed model.
In the real world, in addition to the fear effect of the prey, the clever predator also has spatial-memory and cognition [11], which is often ignored in modeling the predator-prey interaction. For example, blue whales rely on memory for migration, which is presented by B. Abrahms et al. [12] and W. F. Fagan [13]. As another example, animals in polar regions usually determine their spatial movement by judging footprints, which record the history of species distribution and movement, involving time delay [14]. Obviously, highly developed animals can even remember the historical distribution or cluster of species in space. Much progress has been made in implicitly integrating spatial cognition or memory [14,15,16,17,18]. To incorporate the memory effect, Shi et al. proposed a single specie model with spatial memory by introducing a additional delayed diffusion term [14]. They supposed that in addition to the negative gradient of the density distribution function at the present time, there is a directed movement toward the negative gradient of the density distribution function at past time [14]. After this pioneering work, some recent works [19,20,21,22,23] about the population model with memory effect have emerged. In [23], Song et al. obtained a computing method for the normal forms of the Hopf bifurcations in the diffusive predator-prey model with memory effect, which is friendly to use.
Inspired by the above work, we suppose the predator has spatial-memory diffusion and the prey has fear effect, then modified the model (1.3) as follow
All parameters are positive. The biological description of parameters is given in Table 1. The term −d∇(v∇u(t−τ)) represents the memory-based diffusion effect of the predator. The Neumann boundary conditions is used. The aim of this paper is to study the effect of predator's memory-based diffusion and prey's fear on the model (1.4).
The rest of this paper is organized as follows. In Section 2, the stability of coexisting equilibrium and the existence of Hopf bifurcation are considered. In Section 3, the property of Hopf bifurcation is studied. In Section 3, some numerical simulations are given to analyze the effect of spatial-memory and fear effect. In Section 4, a short conclusion is given.
2.
Stability analysis
For simplicity, we choose Ω=(0,lπ). Denote N as positive integer set, and N0 as nonnegative integer set. It is easy to obtain (0,0) and (1,0) are two boundary equilibria of model (1.4). Next, we will give the existence of coexisting equilibrium.
Lemma 2.1. If a>(1+b)μ, model (1.4) has one unique coexisting equilibrium (u∗,v∗) where u∗=μa−bμ, v∗=12c(−1+√1+4c[a−(1+b)μ](a−bμ)2).
Proof. The coexisting equilibrium of (1.4) is a positive root of the following equations
From the second equation, we have u=μa−bμ. Substitute it into the first equation, we have
Then v=12c(−1±√1+4c(a−(1+b)μ)(a−bμ)2). Obviously, the conclusion holds.
In this paper, we mainly study the stability of coexisting equilibrium E∗(u∗,v∗). Linearize model (1.4) at E∗(u∗,v∗), we have
where
and a1=u∗(abv∗(1+bu∗)2−11+cv∗), a2=−u∗(1−u∗)(1+2cv∗)v∗(1+cv∗)2<0, b1=av∗(1+bu∗)2>0. The characteristic equations are
where
2.1. τ=0
When τ=0, the characteristic Eq (2.3) become
where Bn+Cn=−a2b1−(a2dv∗+a1d2)n2l2+d1d2n4l4. Make the following hypothesis
Theorem 2.1. For model (1.4) with τ=0, E∗(u∗,v∗) is locally asymptotically stable under the hypothesis (H1).
Proof. If (H1) holds, we can easily obtain that An>0 and Bn+Cn>0. Then the characteristic roots of (2.4) all have negative real parts. Then E∗(u∗,v∗) is locally asymptotically stable.
2.2. τ>0
In the following, we assume (H1) holds. Let iω (ω>0) be a solution of Eq (2.3), then we have
We can obtain cosωτ=ω2−BnCn, sinωτ=AnωCn>0 under hypothesis (H1). It leads to
Let p=ω2, then (2.5) becomes
and the roots of (2.6) are p±n=12[−(A2n−2Bn)±√(A2n−2Bn)2−4(B2n−C2n)]. By direct computation, we have
and Bn+Cn>0 under hypothesis (H1). Define z±=−(a2dv∗−a1d2)±√(a2dv∗−a1d2)2−4d1d2(−a2b1)2d1d2, d∗=a1d2a2v∗+2v∗√−b1d1d2a2, and M={n|n2l2∈(z−,z+),n∈N0}. Then we can obtain that
The existence of purely imaginary roots of Eq (2.3) can be divided into the following two cases.
Case1:a21+2a2b1>0. We can obtain A2n−2Bn>a21+2a2b1>0. For d>d∗ and n∈M, Eq (2.3) has a pair of purely imaginary roots ±iω+n at τj,+n for j∈N0 and n∈M. Otherwise, Eq (2.3) does not have characteristic roots with zero real parts.
Case2:a21+2a2b1<0. This case can be divided into the following two subcases.
∙ For d≤d∗ and n∈M1:={n|A2n−2Bn<0,(A2n−2Bn)2−4(B2n−C2n)>0,n∈N0}, Eq (2.3) has two pairs of purely imaginary roots ±iω±n at τj,±n for j∈N0 and n∈M1. Otherwise, Eq (2.3) does not have characteristic roots with zero real parts.
∙ For d>d∗ and n∈M2:={n|A2n−2Bn<0,(A2n−2Bn)2−4(B2n−C2n)>0,n∈N0,n∉M}, Eq (2.3) has two pairs of purely imaginary roots ±iω±n at τj,±n for j∈N0 and n∈M1. For d>d∗ and n∈M, Eq (2.3) has a pair of purely imaginary roots ±iω+n at τj,+n for j∈N0 and n∈M. Otherwise, Eq (2.3) does not have characteristic roots with zero real parts.
The ω±n and τj,±n are defined as follow
Define
We have the following lemma.
Lemma 2.2. Assume (H1) holds. ThenRe(dλdτ)|τ=τj,+n>0, Re(dλdτ)|τ=τj,−n<0 for τj,±n∈S and j∈N0.
Proof. By Eq (2.3), we have
Then
Therefore Re(dλdτ)|τ=τj,+n>0, Re(dλdτ)|τ=τj,−n<0.
Denote τ∗=min{τ0,±n|τ0,±n∈S}. We have the following theorem.
Theorem 2.2. Assume (H1) holds, then the following statements are true for model (1.4).
∙ E∗(u∗,v∗) is locally asymptotically stable for τ>0 when S=∅.
∙ E∗(u∗,v∗) is locally asymptotically stable for τ∈[0,τ∗) when S≠∅.
∙ E∗(u∗,v∗) is unstable for τ∈(τ∗,τ∗+ε) for some ε>0 when S≠∅.
∙ Hopf bifurcation occurs at(u∗,v∗) when τ=τj,+n (τ=τj,−n), j∈N0, τj,±n∈S.
Remark 2.1. In the Theorem 2, if E∗(u∗,v∗) is locally asymptotically stable, the densities of prey and predator will tend to the equilibrium state in the whole region when the initial densities of prey and predator is near E∗(u∗,v∗). When Hopf bifurcation occurs at(u∗,v∗), then the densities of prey and predator will produce periodic oscillation. Especially, spatially homogeneous periodic oscillations may occur when τ near the critical value τ=τj,+0 or τ=τj,−0, and spatially inhomogeneous periodic oscillations may occur when τ near the critical value τ=τj,+n or τ=τj,−n (n>0).
3.
Property of Hopf bifurcation
In this section, we use the algorithm in [23] to compute the normal form of Hopf bifurcation. We denote the critical value of Hopf bifurcation as ˜τ and the purely imaginary roots as ±iωn of Eq (2.3). Let ˉu(x,t)=u(x,τt)−u∗ and ˉv(x,t)=v(x,τt)−v∗. Drop the bar, the model (1.4) can be written as
Define the real-valued Sobolev space X={U=(u,v)T∈W2,2(0,lπ)2,(∂u∂x,∂v∂x)|x=0,lπ=0}, the inner product
and C=C([−1,0];X). Set τ=˜τ+ε, where ε is small perturbation. Then system (3.1) is rewritten as
where for φ=(φ,φ2)T∈C, d(ε)Δ, L(ε):C→X, F:C×R2→X. They are defined as
and
Denote L0(φ)=˜τAφ(0), and rewrite (3.2) as
where ˜F(φ,ε)=εAφ(0)+F(φ,ε)+Fd(φ,ε). The characteristic equation for the linearized equation dU(t)dt=d0Δ(Ut)+L0(Ut) is ˜Γn(λ)=det(˜Mn((λ))), where
The eigenvalue problem
has eigenvalues n2l2 and normalized eigenfunctions
Set β(j)n=zn(x)ej, j=1,2, where e1=(1,0)T and e2=(0,1)T. Define ηn(θ)∈BV([−1,0],R2), such that
C=C([−1,0],R2), C∗=C([0,1],R2∗), and
Let ⋀={i˜ω,−i˜ω}, the eigenspace P, and corresponding adjoint space P∗. Decompose C=P⊕Q, where Q={φ∈C:<ψ,φ>=0,∀ψ∈P∗}. Choose Φ(θ)=(ϕ(θ),ˉϕ(θ)), Ψ(θ)=col(ψT(s),ˉψT(s)), where
and
Then ϕ(θ) and ψ(s) are the bases of P and P∗, respectively, and such that <ϕ,ψ>=I2.
By direct computation, we have
where f(1)20=2abv∗(1+bu∗)3−21+cv∗, f(1)11=−a(1+bu∗)2+c(−1+2u∗)(1+cv∗)2, f(1)02=2c2(1−u∗)u∗(1+cv∗)3, f(1)30=−6ab2v∗(1+bu∗)4, f(1)21=2ab(1+bu∗)3+2c(1+cv∗)2, f(1)12=2c2(1−2u∗)(1+cv∗)3, f(1)03=−6c3(1−u∗)u∗(1+cv∗)4, f(2)20=−2abv∗(1+bu∗)3, f(2)11=a(1+bu∗)2, f(2)02=0, f(2)30=6ab2v∗(1+bu∗)4, f(2)21=−2ab(1+bu∗)3, f(2)12=0, f(2)03=0. We can compute the following parameters
and ˜Aj1j2=Aj1j2−2n2l2Adj1j2 for j1,j2=0,1,2,j1+j2=2. In addition, h0,20(θ)=1lπ(˜M0(2i˜ω))−1A20e2i˜ωθ, h0,11(θ)=1lπ(˜M0(0))−1A11, h2n,20(θ)=12lπ(˜M2n(2i˜ω))−1˜A20e2i˜ωθ, h2n,11(θ)=1lπ(˜M2n(0))−1˜A11.
Then we have
where b(1)2n=−n2l2, b(2)2n=−2n2l2, b(3)2n=−4n2l2. The normal form of the Hopf bifurcation is
where
By the coordinate transformation z1=ω1−iω2, z2=ω1+iω2, and ω1=ρcosξ, ω2=ρsinξ, the normal form (3.8) can be rewritten as
where K1=12Re(B1), K2=13!Re(B2).
By the work [23], we have the following theorem.
Theorem 3.1. If K1K2<0(>0), the Hopf bifurcation is supercritical (subcritical), and the bifurcating periodic orbits is stable(unstable) for K2<0(>0).
Remark 3.1. In the Theorem 3.1, when Hopf bifurcation is supercritical (subcritical), then the bifurcating periodic solutions exist for τ>˜τ (τ<˜τ), where ˜τ is some critical value τ=τj,±n. When the periodic solution is stable, the densities of prey and predator will produce periodic oscillation, and finally continue to oscillate.
4.
Numerical simulations
In this section, we give some numerical simulations to analyze the effect of spatial memory in predator and fear in prey on the model (1.4). Fix the following parameters
4.1. The effect of d
If we choose c=1, then model (1.4) has a unique coexisting equilibrium (u∗,v∗)≈(0.6667,0.0.6667), and a21+2a2b1≈0.0352>0, d∗≈0.6206. To study the effect of memory-based diffusion coefficient d on the model (1.4), we give the bifurcation diagram of model (1.4) with parameter d as in Figure 1. By the Theorem 2.2, we know that (u∗,v∗) is locally stable for τ≥0 when d<d∗. But when d>d∗, the inhomogeneous Hopf bifurcation curves exist. This means that increasing parameter d is not conducive to the stability of the equilibrium (u∗,v∗), and the densities of prey and predator will produce spatially inhomogeneous periodic oscillation.
4.2. The effect of τ
Choose d=0.7, we have M={2,3} and τ∗=τ0,+2≈17.4593<τ0,+3≈19.1380. When τ∈[0,τ∗), (u∗,v∗) is locally stable (Figure 2). By direct calculation, we can obtain K1≈0.0166, and K2≈−0.0699. Then, the Hopf bifurcation is supercritical and the bifurcating periodic solution is stable (Figure 3). At this time, the bifurcating periodic solution is spatially inhomogeneous and with mode-2.
Choose d=0.8, we have M={2,3} and τ∗=τ0,+3≈8.4754<τ0,+2≈9.9645. When τ∈[0,τ∗), (u∗,v∗) is locally stable (Figure 4). By direct calculation, we can obtain K1≈0.1236, K2≈−0.3994. Then, the Hopf bifurcation is supercritical and the bifurcating periodic solution is stable (Figure 5). At this time, the bifurcating periodic solution is spatially inhomogeneous with mode-3. When τ∗=τ0,+3<τ0,+2<τ=10, there is an unstably spatially inhomogeneous periodic solution with mode-2 which transitions to the stably spatially inhomogeneous periodic solution with mode-3 (Figure 6).
4.3. The effect of c
Next, we will study the effect of fear effect c on the model (1.4). Fix the parameters as (4.1), then model (1.4) has a unique coexisting equilibrium (u∗,v∗). And a21+2a2b1>0 when 0<c<2.6194. We give the figure of d∗ with parameter c as in Figure 7. Set parameter d=0.7 and d=0.8, we give the bifurcation diagrams of model (1.4) with parameter c as in Figure 8.
When d=0.7 and τ=20, increasing parameter c can destroy the stability of the coexisting equilibrium (u∗,v∗), and induce spatially inhomogeneous periodic solution (Figure 9). This means that increasing parameter c is not conducive to the stability of the coexisting equilibrium (u∗,v∗).
When d=0.8 and τ=9, increasing parameter c can destroy the stability of the coexisting equilibrium (u∗,v∗), and induce spatially inhomogeneous periodic solution initially (Figure 10). But when c is larger enough, increasing parameter c can rule out the spatially inhomogeneous periodic oscillation and stabilize the coexisting equilibrium (u∗,v∗) (Figure 10). This means that increasing parameter c is not conducive to the stability of the equilibrium (u∗,v∗), initially. But when c is large, increasing parameter c is conducive to the stability of the coexisting equilibrium (u∗,v∗).
5.
Conclusions
In this paper, we incorporate the memory effect in predator and fear effect in prey into a predator-prey model. By using time delay in the memory of predator as bifurcating parameter, we analyze the local stability of coexisting equilibrium, the existence of Hopf bifurcation, and the property of Hopf bifurcation by the method in [23]. Through the numerical simulations, we analyzed the effect of memory effect in predator and fear in prey on the model.
The spatial memory effect plays an important role in the dynamics of the predator-prey model. Through the numerical simulations, we observed that the memory-based diffusion coefficient d has destabilizing effect on the predator-prey model when it is larger than some critical value. In addition. when d crosses the critical value, time delay τ in the memory of predator can affect the stability of the equilibrium (u∗,v∗). In the numerical simulations, we observe that the first Hopf bifurcation curve is inhomogeneous bifurcation curve, and homogeneous Hopf bifurcation curve does not exist. This is different from the predator-prey model without the spatial memory effect. When τ crosses the critical value τ∗, the densities of prey and predator will produce spatially inhomogeneous periodic oscillation. When τ crosses the second critical value, the spatially inhomogeneous periodic oscillations with different modes exist, but the densities of prey and predator will converge to the spatially inhomogeneous periodic solution corresponding to the first bifurcation curve. This shows that the spatially memory effect in predator can destroy the stability of the coexisting equilibrium, and induce spatially inhomogeneous periodic oscillations.
In addition, the fear effect parameter c in prey can also affect the stability of the coexisting equilibrium (u∗,v∗). A small fear effect parameter c means a large birth rate 11+cv, then the large birth rate can support fluctuations. Increasing parameter c can destroy the stability of the coexisting equilibrium (u∗,v∗), and induce spatially inhomogeneous periodic solution. Hence, we observed the destabilizing effect on the the coexisting equilibrium (u∗,v∗). A large fear effect parameter c means a low birth rate, then the low birth rate can not support fluctuations. Increasing parameter c can rule out the spatially inhomogeneous periodic oscillation and stabilize the coexisting equilibrium (u∗,v∗). Hence, we observed the stabilizing effect on the the coexisting equilibrium (u∗,v∗).
Acknowledgements
This research is supported by the Fundamental Research Funds for the Central Universities (Grant No. 2572022BC01) and Postdoctoral Program of Heilongjiang Province (No. LBH-Q21060).
Conflict of interest
The authors have no relevant financial or non-financial interests to disclose.