A stochastic delay predator-prey model in a polluted environment with impulsive
toxicant input is proposed and studied. The thresholds between stability
in time average and extinction of each population are obtained. Some
recent results are extended and improved greatly. Several simulation figures
are introduced to support the conclusions.
1. Introduction
Environmental pollution by industries, agriculture, and other human activities is one of the most important socio-ecological problems in the world today. Due to toxins in the environment, lots of species have gone extinct, and many are on the verge of extinction. Thus, controlling the environmental pollution and the conservation of biodiversity are the major focus areas of all the countries around the world. This motivates scholars to study the effects of toxins on populations and to find out a theoretical persistence-extinction threshold.
Recently, a lot of population models in a polluted environment have been proposed and investigated; here, we may mention, among many others, [1–23]. Particularly, Yang et al. [15] pointed out that in many cases toxicants should be emitted in regular pulses, for example, the use of pesticides and the pollution by heavy metals (see, e.g., [24]). Thus, they proposed the following two-species Lotka-Volterra predator-prey system in a polluted environment with impulsive toxicant input:
(1)dx1(t)dt=x1(t)[r10-r11C10(t)-a11x1(t)-a12x2(t)],dx2(t)dt=x2(t)[-r20-r21C20(t)+a21x1(t)-a22x2(t)],dC10(t)dt=k1Ce(t)-(g1+m1)C10(t),dC20(t)dt=k2Ce(t)-(g2+m2)C20(t),dCe(t)dt=-hCe(t),t≠nγ,n∈Z+,Δxi(t)=0,ΔCi0(t)=0,ΔCe(t)=b,t=nγ,n∈Z+,i=1,2,
where all the parameters are positive constants andΔf(t)=f(t+)-f(t), Z+={1,2,…};x1(t) and x2(t): the size of prey population and the predator population, respectively; ri0: the intrinsic growth rate of the ith population without toxicant; ri1: the ith population response to the pollutant present in the organism; Ci0(t): the concentration of toxicant in the ith organism; Ce(t): the concentration of toxicant in the environment; kCe(t): the organism’s net uptake of toxicant from the environment; gCi0(t)+mCi0(t): the egestion and depuration rates of the toxicant in the ith organism; hCe(t): the toxicant loss from the environment itself by volatilization and so on; γ: the period of the impulsive effect about the exogenous input of toxicant; b: the toxicant input amount at every time.
Yang et al. [15] showed that in the following Lemma holds.
Lemma 1.
For system (1), define
(2)Δ2=r10a21-r20a11,Δ¯2:=a21r11K1γ+a11r21K2γ,Ki=kibh(gi+mi).
If r10<r11K1/γ, then limt→+∞xi(t)=0, i=1,2.
If r10>r11K1/γ and Δ2<Δ¯2, then limsupt→+∞t-1∫0tx1(s)ds>0 and x2(t) goes to extinction.
If Δ2>Δ¯2, then limsupt→+∞t-1∫0txi(s)ds>0,i=1,2.
Some interesting and important problems arise naturally.
(Q1) In the real world, the growth of species depends on various environmental factors, such as temperature, humidity and parasites and so forth. Therefore population models should be stochastic rather than deterministic (May [25]). Thus, what happens if model (1) is subject to stochastic noises?
(Q2) In addition, time delays occur in almost every situation. Kuang [26] has pointed out that ignoring time delays means ignoring reality. Therefore, what happens if model (1) takes time delays into account?
(Q3) Can we improve the results given in Lemma 1?
The aim of this paper is to study the above problems. Suppose that stochastic noises mainly affect the growth rates, with ri0→ri0+αiB˙i(t) (see, e.g., [27–39]), where B˙i(t) is a white noise and αi2 is the intensity of the noise. Moreover, taking time delays into account, we obtain the following model:
(3)dx1(t)=x1(t)[r10-r11C10(t)-a11x1(t)-a12x2(t-τ1)]dt+α1x1(t)dB1(t),dx2(t)=x2(t)[-r20-r21C20(t)-a22x2(t)+a21x1(t-τ2)]dt+α2x2(t)dB2(t),dC10(t)dt=k1Ce(t)-(g1+m1)C10(t),dC20(t)dt=k2Ce(t)-(g2+m2)C20(t),dCe(t)dt=-hCe(t),t≠nγ,n∈Z+,Δxi(t)=0,ΔCi0(t)=0,ΔCe(t)=b,t=nγ,n∈Z+,i=1,2,
with initial condition
(4)xi(t)=ϕi(t)>0,t∈[-τ,0];ϕi(0)>0,i=1,2,
where τi≥0, τ=max{τ1,τ2}, ϕi(t) is continuous on [-τ,0]. Our main result is the following theorem.
Theorem 2.
For system (3), define
(5)θ1=r10-0.5α12,θ2=r20+0.5α22,Δ=a11a22+a12a21,Δ~1:=θ1a22+θ2a12,Δ~2:=θ1a21-θ2a11,Δ¯1:=a22r11K1γ-a12r21K2γ.
If θ1<r11K1/γ, then both x1 and x2 go to extinction almost surely (a.s.); that is, limt→+∞xi(t)=0a.s.,i=1,2.
If θ1>r11K1/γ and Δ~2<Δ¯2, then x2(t) goes to extinction and x1 is stable in time average a.s.; that is,
(6)limt→+∞t-1∫0tx1(s)ds=θ1-r11K1/γa11>0,a.s.
If Δ~2>Δ¯2, then both x1 and x2 are stable in time average a.s.
(7)limt→+∞t-1∫0tx1(s)ds=Δ~1-Δ¯1Δ>0,a.s.limt→+∞t-1∫0tx2(s)ds=Δ~2-Δ¯2Δ>0,a.s.
Remark 3.
By comparing Lemma 1 with our Theorem 2, we can see that on the one hand, if α1=α2=0 and τ1=τ2=0, then θi=ri0, Δ~i=Δi, i=1,2, and our stochastic delay system (3) becomes model (1); on the other hand, our results in Theorem 2 improve that in Lemma 1. Lemma 1 shows that the superior limit is positive, while Theorem 2 reveals that the limit exists and gives the explicit form of the limit. The contribution of this paper is therefore clear.
2. Proof
For the sake of simplicity, we introduce some notations:
(8)R+2={a=(a1,a2)∈R2∣ai>0,i=1,2},〈f(t)〉=t-1∫0tf(s)ds;〈f(t)〉*=limsupt→+∞t-1∫0tf(s)ds,〈f(t)〉*=liminft→+∞t-1∫0tf(s)ds.
Lemma 4.
For any given initial value ϕ(t)=(ϕ1(t),ϕ2(t))∈C([-τ,0],R+2), there is a unique global positive solution x(t)=(x1(t),x2(t))T to the first two equations of system (3) a.s.
Proof.
The proof is similar to Hung [29] by defining
(9)V(x)=[x1-1-lnx1]+a11a22a212[x2-1-lnx2]+0.5a11∫t-τ2tx12(s)ds
and hence is omitted.
To begin with, let us consider the following subsystem of (3):
(10)dC10(t)dt=k1Ce(t)-(g1+m1)C10(t),dC20(t)dt=k2Ce(t)-(g2+m2)C20(t),dCe(t)dt=-hCe(t),t≠nγ,n∈Z+,ΔC10(t)=0,ΔC20(t)=0,ΔCe(t)=b,t=nγ,n∈Z+,0≤C0(0)≤1,0≤Ce(0)≤1.
Lemma 5 (see [13, 15]).
System (10) has a unique positive γ-periodic solution (C~10(t),C~20(t),C~e(t))T, and for each solution (C10(t),C20(t),Ce(t))T of (10), C10(t)→C~10(t), C20(t)→C~20(t), and Ce(t)→C~e(t) as t→∞. Moreover, Ci0(t)>C~i0(t) and Ce(t)>C~e(t) for all t≥0 if Ci0(0)>C~i0(0) and Ce(0)>C~e(0), i=1,2, where
(11)C~10(t)=C~10(0)e-(g1+m1)(t-nγ)+k1b(e-(g1+m1)(t-nγ)-e-h(t-nγ))(h-g1-m1)(1-e-hγ),C~20(t)=C~20(0)e-(g2+m2)(t-nγ)+k2b(e-(g2+m2)(t-nγ)-e-h(t-nγ))(h-g2-m2)(1-e-hγ),C~e(t)=be-h(t-nγ)1-e-hγ,C~10(0)=k1b(e-(g1+m1)γ-e-hγ)(h-g1-m1)(1-e-(g1+m1)γ)(1-e-hγ),C~20(0)=k2b(e-(g2+m2)γ-e-hγ)(h-g2-m2)(1-e-(g2+m2)γ)(1-e-hγ),C~e(0)=b1-e-hγ
for t∈(nγ,(n+1)γ] and n∈Z+. In addition,
(12)limt→+∞t-1∫0tC~i0(s)ds=kibh(gi+mi)γ=Kiγ,i=1,2.
Lemma 6 (see [34]).
Suppose that x(t)∈C[Ω×[0,+∞),R+].
(I) If there exist σ and positive constants σ0,T such that
(13)lnx(t)≤σt-σ0∫0tx(s)ds+∑i=1nβiBi(t)
for t≥T, where Bi(t) are independent standard Brownian motions and βi are constants, 1≤i≤n, then one has the following: if σ≥0, then 〈x〉*≤σ/σ0 a.s.; if σ<0, then limt→+∞x(t)=0a.s.
(II) If there exist positive constants σ0,T and σ such that
(14)lnx(t)≥σt-σ0∫0tx(s)ds+∑i=1nβiBi(t)
for t≥T, then 〈x〉*≥σ/σ0 a.s.
Now, let us consider the following auxiliary system:
(15)dy1(t)=y1(t)[r10-r11C10(t)-a11y1(t)]dt+α1y1(t)dB1(t),dy2(t)=y2(t)[-r20-r21C20(t)+a21y1(t-τ2)-a22y2(t)]dt+α2y2(t)dB2(t),
with initial value ϕ(t)∈C([-τ,0],R+2).
Lemma 7.
If θ1=r10-0.5α12>r11K1/γ, then the solution y(t) of system (15) obeys
(16)limt→+∞〈y1(t)〉=θ1-r11K1/γa11;limt→+∞y2(t)=0a.s.,ifΔ~2<Δ¯2;limt→+∞〈y2(t)〉=Δ~2-Δ¯2a11a22a.s.,ifΔ~2>Δ¯2.
Proof.
By Lemma 5,
(17)limt→+∞t-1∫0tCi0(s)ds=limt→+∞t-1∫0tC~i0(s)ds=Kiγ,i=1,2.
Then, for allε>0, there exists T>0 such that
(18)Kiγ-ε≤〈Ci0(t)〉≤Kiγ+ε,t>T,i=1,2.
An application of Itô’s formula to (15) yields
(19)lny1(t)-lny1(0)=θ1t-r11∫0tC10(s)ds-a11∫0ty1(s)ds+α1B1(t),lny2(t)-lny2(0)=-θ2t-r21∫0tC20(s)ds+a21∫0ty1(s-τ2)ds-a22∫0ty2(s)ds+α2B2(t)=-θ2t-r21∫0tC20(s)ds+a21∫0ty1(s)ds-a21[∫t-τ2ty1(s)ds-∫-τ20y1(s)ds]-a22∫0ty2(s)ds+α2B2(t).
That is to say, we have shown that
(20)t-1lny1(t)y1(0)=θ1-r11〈C10(t)〉-a11〈y1(t)〉+t-1α1B1(t),(21)t-1lny2(t)y2(0)+t-1a21[∫t-τ2ty1(s)ds-∫-τ20y1(s)ds]=-θ2-r21〈C20(t)〉+a21〈y1(t)〉-a22〈y2(t)〉+t-1α2B2(t).
When (18) is used in (20), we can see that for t>T,
(22)t-1lny1(t)y1(0)≤θ1-r11K1γ+r11ε-a11〈y1(t)〉+t-1α1B1(t),(23)t-1lny1(t)y1(0)≥θ1-r11K1γ-r11ε-a11〈y1(t)〉+t-1α1B1(t).
Let ε be sufficiently small such that θ1-r11K1/γ-r11ε>0. Making use of (I) and (II) in Lemma 6 to (22) and (23), respectively, we have
(24)〈y1(t)〉*≤θ1-r11K1/γ+r11εa11,〈y1(t)〉*≥θ1-r11K1/γ-r11εa11.
It then follows from the arbitrariness of ε that
(25)limt→+∞〈y1(t)〉=θ1-r11K1/γa11.
Substituting (17) and (25) into (20) and noting that limt→+∞t-1B1(t)=0, one can derive that
(26)limt→+∞t-1lny1(t)=0,a.s.
Employing (20) and (21) in the expression a21ln(y1(t)/y1(0))+a11ln(y2(t)/y2(0)) yields
(27)a11t-1lny2(t)y2(0)+a21t-1lny1(t)y1(0)=Δ~2-r11a21〈C10(t)〉-r21a11〈C20(t)〉-a11a22〈y2(t)〉-t-1a11a21[∫t-τ2ty1(s)ds-∫-τ20y1(s)ds]+t-1[a21α1B1(t)+a11α2B2(t)].
In view of (25), we get
(28)limt→+∞t-1∫t-τ2ty1(s)ds=limt→+∞t-1(∫0ty1(s)ds-∫0t-τ2y1(s)ds)=0,a.s.
By (17), (26), (27), and (28), for allε>0, there exists T>0 such that, for t≥T,
(29)a11t-1lny2(t)y2(0)≤Δ~2-Δ¯2+ε-a11a22〈y2(t)〉+t-1[a21σ1B1(t)+a11σ2B2(t)],(30)a11t-1lny2(t)y2(0)≥Δ~2-Δ¯2-ε-a11a22〈y2(t)〉+t-1[a21σ1B1(t)+a11σ2B2(t)].
If Δ~2<Δ¯2, then we can choose ε sufficiently small such that Δ~2-Δ¯2+ε<0. Then, by (29) and (I) in Lemma 6, we obtain limt→+∞y2(t)=0a.s. If Δ~2>Δ¯2, then we can choose ε sufficiently small such that Δ~2-Δ¯2-ε>0. An application of (I) and (II) in Lemma 6 to (29) and (30), respectively, makes one observe that
(31)Δ~2-Δ¯2-εa11a22≤〈y2(t)〉*≤〈y2(t)〉*≤Δ~2-Δ¯2+εa11a22,a.s.
Therefore, using the arbitrariness of ε results in
(32)limt→+∞〈y2(t)〉=Δ~2-Δ¯2a11a22a.s.
This completes the proof.
We are now in the position to prove our main results.
Proof of Theorem 2.
Applying Itô’s formula to (3) leads to
(33)lnx1(t)-lnx1(0)=θ1t-r11∫0tC10(s)ds-a11∫0tx1(s)ds-a12∫0tx2(s-τ1)ds+α1B1(t)=θ1t-r11∫0tC10(s)ds-a12∫0tx2(s)ds+a12[∫t-τ1tx2(s)ds-∫-τ10x2(s)ds]-a11∫0tx1(s)ds+α1B1(t).(34)lnx2(t)-lnx2(0)=-θ2t-r21∫0tC20(s)ds+a21∫0tx1(s)ds-a21[∫t-τ2tx1(s)ds-∫-τ20x1(s)ds]-a22∫0tx2(s)ds+α2B2(t).
(i) It follows from (17) and (33) that
(35)t-1lnx1(t)-t-1lnx1(0)≤θ1-r11〈C10(t)〉-a11〈x1(t)〉+α1B1(t)t≤θ1-r11K1λ+ε-a11〈x1(t)〉+α1B1(t)t
for sufficiently large t. Since θ1-r11K1/λ<0, then we can choose ε sufficiently small such that θ1-r11K1/λ+ε<0. Then, by (I) in Lemma 6,
(36)limt→+∞x1(t)=0,a.s.
When (36) is used in (34), one can see that
(37)t-1lnx2(t)-lnx2(0)≤-θ2+ε-a22〈x2(t)〉+α2B2(t)t
for sufficiently large t, where ε>0 obeys -θ2+ε<0. In view of Lemma 6 again, limt→+∞x2(t)=0,a.s.
(ii) By the stochastic comparison theorem [40], one can observe that
(38)x1(t)≤y1(t),x2(t)≤y2(t).
Note that θ1>r11K1/γ and Δ~2<Δ¯2; it then follows from Lemma 7 that limt→+∞y2(t)=0,a.s. Making use of (38) gives limt→+∞x2(t)=0,a.s. Thus, for allε>0, there exists T>0 such that, for t≥T,
(39)ε2≤a12x2(t)≤ε2.
Substituting the above inequalities into (33) and then using (18), we obtain
(40)t-1lnx1(t)≤t-1lnx1(0)+θ1-r11〈C10(t)〉-a11〈x1(t)〉+ε2+α1B1(t)t≤θ1-r11K1γ+2ε-a11〈x1(t)〉+α1B1(t)t,(41)t-1lnx1(t)≥t-1lnx1(0)+θ1-r11〈C10(t)〉-a11〈x1(t)〉-ε2+α1B1(t)t≥θ1-r11K1γ-2ε-a11〈x1(t)〉+α1B1(t)t.
Let ε be sufficiently small such that θ1-r11K1/γ-ε>0, and then, applying (I) and (II) in Lemma 6 to (40) and (41), respectively, one can see that
(42)θ1-r11K1/γ-2εa11≤〈x1(t)〉*≤〈x1(t)〉*≤θ1-r11K1/γ+2εa11a.s.
An application of the arbitrariness of ε gives
(43)limt→+∞〈x1(t)〉=θ1-r11K1/γa11,a.s.
(iii) Clearly, Δ~2>Δ¯2 implies θ1>r11K1/γ, and then, by Lemma 7,
(44)limt→+∞〈y2(t)〉=Δ~2-Δ¯2a11a22.
Thus, similar to the proof of (28), we get
(45)limt→+∞t-1∫t-τ1ty2(s)ds=0a.s.
Therefore, by (26), (28), and (38), we can observe that
(46)limsupt→+∞t-1lnx1(t)≤limt→+∞t-1lny1(t)=0,(47)limt→+∞t-1∫t-τ2tx1(s)ds=0,limt→+∞t-1∫t-τ1tx2(s)ds=0,a.s.
Employing (33) and (34) in the expression a21ln(x1(t)/x1(0))+a11ln(x2(t)/x2(0)) yields
(48)t-1a21lnx1(t)x1(0)+t-1a11lnx2(t)x2(0)=a12a21t-1[∫t-τ1tx2(s)ds-∫-τ10x2(s)ds]-a11a21t-1[∫t-τ2tx1(s)ds-∫-τ20x1(s)ds]+Δ~2-a21r11〈C10(t)〉-a11r21〈C20(t)〉-Δ〈x2(t)〉+t-1a21α1B1(t)+t-1a11α2B2(t).
When (18), (46) and (47), are used in (48), one can obtain
(49)t-1a11lnx2(t)x2(0)≥Δ~2-Δ¯2-ε-Δ〈x2(t)〉+t-1a21α1B1(t)+t-1a11α2B2(t)
for sufficiently large t, where ε>0 obeys Δ~2-Δ¯2-ε>0. It then follows from (II) in Lemma 6 that
(50)〈x2(t)〉*≥Δ~2-Δ¯2-εΔ.
By virtue of the arbitrariness of ε, we can see that
(51)〈x2(t)〉*≥Δ~2-Δ¯2Δ.
Consequently, for every 0<ε<a12(Δ~2-Δ¯2)/Δ, there is T>0 such that
(52)a12〈x2(t)〉≥a12〈x2〉*-ε≥a12(Δ~2-Δ¯2)Δ-ε,t>T.
Substituting the above inequality into (33) and then using (18) and (47), one can see that
(53)t-1lnx1(t)x1(0)≤θ1-a12(Δ~2-Δ¯2)Δ+3ε-a11t-1∫0tx1(s)ds+t-1α1B1(t)=a11(Δ~1-Δ¯1)Δ+3ε-a11t-1∫0tx1(s)ds+t-1α1B1(t)
for sufficiently large t. Since Δ~1-Δ¯1>0, and then, by Lemma 6 and the arbitrariness of ε, one can observe that
(54)〈x1(t)〉*≤Δ~1-Δ¯1Δ.
When this inequality, (18) and (47), are used in (34), we can see that
(55)t-1lnx2(t)x2(0)≤-θ2+a21Δ~1-Δ¯1Δ+3ε-a22t-1∫0tx2(s)ds+t-1α2B2(t)=a22(Δ~2-Δ¯2)Δ+3ε-a22t-1∫0tx2(s)ds+t-1α2B2(t)
for sufficiently large t. Then, it follows from Lemma 6 and the arbitrariness of ε that
(56)〈x2(t)〉*≤Δ~2-Δ¯2Δ.
Substituting the above inequality and (18) into (33), we get
(57)t-1lnx1(t)x1(0)≥θ1-a12Δ~2-Δ¯2Δ-3ε-a11t-1∫0tx1(s)ds+t-1α1B1(t)=a11(Δ~1-Δ¯1)Δ-3ε-a11t-1∫0tx1(s)ds+t-1α1B1(t)
for sufficiently large t. By (II) in Lemma 6 and the arbitrariness of ε again, we obtain
(58)〈x1(t)〉*≥Δ~1-Δ¯1Δ.
Then, the required assertion follows from (51), (54), (56), and (58).
3. Numerical Simulations
Let us use the famous Milstein method (see, e.g., [41]) to illustrate the analytical results.
To begin with, we choose r10=0.85, r20=0.05, r11=r21=1, a11=0.4, a12=0.4, a21=0.3, a22=0.3, τ1=3,τ2=8, α22=0.1, ki=gi=mi=0.1, i=1,2, h=0.5, b=0.6, and γ=12. Then,
(59)Ki=kibh(gi+mi)=0.6,Δ2=r10a21-r20a11=0.235>Δ¯2=a21r11K1γ+a11r21K2γ=0.035.
By (c) in Lemma 1, the solution of model (1) obeys
(60)limsupt→+∞t-1∫0tx1(s)ds>0,limsupt→+∞t-1∫0tx2(s)ds>0.
However, when the white noises are taken into account, the properties of the system may be changed greatly. In Figure 1, we let the coefficients be same with the above. The only difference between conditions of Figures 1(a), 1(b), and 1(c) is that the value of α12 is different. In Figure 1(a), we choose α12/2=0.82. Therefore,
(61)θ1=r10-α122=0.03<r11K1γ=0.05.
Then, by (i) in Theorem 2, both x1 and x2 are extinctive. Figure 1(a) confirms these. In Figure 1(b), we choose α12/2=0.65. That is to say θ1=0.2>r11K1/γ=0.05 and Δ~2=0.02<Δ¯2=0.035. It then follows from (ii) in Theorem 2 that x2 is extinctive and x1 is stable in time average:
(62)limt→+∞〈x1(t)〉=θ1-r11K1/γa11=0.375.
See Figure 1(b). In Figure 1(c), we choose α12/2=0.2. Then, Δ~2=0.155>Δ¯2=0.035. In view of (iii) in Theorem 2, we can obtain that both x1 and x2 are stable in time average:
(63)limt→+∞〈x1(t)〉=Δ~1-Δ¯1Δ=0.240.24=1,limt→+∞〈x2(t)〉=Δ~2-Δ¯2Δ=0.120.24=0.5.
Figure 1(c) confirms these.
Solutions of system (3) for r10=0.85, r20=0.05, r11=r21=1, a11=0.4, a12=0.4, a21=0.3, a22=0.3, τ1=3, τ2=8, α22=0.1, ki=gi=mi=0.1, i=1,2, h=0.5, b=0.6, γ=12, x1(0)=0.9, x2(0)=0.5, C0(0)=Ce(0)=0.1, and step size Δt=0.001. (a) is with α12/2=0.82; (b) is with α12/2=0.65; (c) is with α12/2=0.2.
In Figure 2, we choose r10=0.85, r20=0.05, r11=r21=1, a11=0.4, a12=0.4, a21=0.3, a22=0.3, τ1=3, τ2=8, α12=0.4, α22=0.1, ki=gi=mi=0.1, i=1,2, h=0.5, and b=0.6. The only difference between conditions of Figures 1(c) and 2 is that the value of γ is different. In Figure 2, we choose γ=0.8. Then, θ1=0.65<r11K1/γ=0.75. It follows from (i) in Theorem 2 that both x1 and x2 are extinctive. Figure 2 confirms these. By comparing Figure 1(c) with Figure 2, one can see that the impulsive period γ plays a key role in determining the stability in time average and the extinction of the species.
Solutions of system (3) for r10=0.85, r20=0.05, r11=r21=1, a11=0.4, a12=0.4, a21=0.3, a22=0.3, τ1=3, τ2=8, α12=0.4, α22=0.1, ki=gi=mi=0.1, i=1,2, h=0.5, b=0.6, γ=0.8, x1(0)=0.9, x2(0)=0.5, C0(0)=Ce(0)=0.1, and step size Δt=0.001.
4. Conclusions and Future Directions
This paper is concerned with stochastic delay predator-prey model in a polluted environment with impulsive toxicant input. For each species, the threshold between stability in time average and extinction is established. Some recent results are improved and extended. Our Theorem 2 reveals some interesting and important results.
Firstly, time delay is harmless for stability in time average and extinction of the stochastic system (3).
The white noise α1dB1(t) and α2dB2(t) can change the properties of the system greatly.
The impulsive period γ plays an important role in determining the stability in time average and the extinction of the species.
Some interesting questions deserve further investigations. One may consider some more realistic but more complex systems, for example, stochastic delay model with Markov switching (see, e.g., [30, 32, 39]). It is also interesting to investigate what happens if aij is stochastic.
Acknowledgments
The author thanks the editor and reviewer for these valuable and important comments. This research is supported by NSFC of China (nos. 11301207, 11171081, 11301112 and 11171056), Natural Science Foundation of Jiangsu Province (No. BK20130411) and Natural Science Research Project of Ordinary Universities in Jiangsu Province (no. 13KJB110002).
HallamT. G.ClarkC. E.LassiterR. R.Effects of toxicants on populations: a qualitative approach—1. Equilibrium environmental exposure19831842913042-s2.0-0020571184HallamT. G.ClarkC. E.JordanG. S.Effects of toxicants on populations: a qualitative approach—2. First order kinetics198318125372-s2.0-0020959624HallamT. G.De LunaJ. T.Effects of toxicants on populations: a qualitative approach—3. Environmental and food chain pathways198410934114292-s2.0-0021217861MaZ. E.CuiG. R.WangW. D.Persistence and extinction of a population in a polluted environment19901011759710.1016/0025-5564(90)90103-6MR1067460ZBL0714.92027FreedmanH. I.ShuklaJ. B.Models for the effect of toxicant in single-species and predator-prey systems1991301153010.1007/BF00168004MR1130786ZBL0825.92125HuapingL.ZhienM.The threshold of survival for system of two species in a polluted environment1991301496110.1007/BF00168006MR1130788ZhienM.WengangZ.ZhixueL.The thresholds of survival for an n-food chain model in a polluted environment1997210244045810.1006/jmaa.1997.5387MR1453184BuonomoB.Di LiddoA.SguraI.A diffusive-convective model for the dynamics of population-toxicant interactions: some analytical and numerical results19991571-2376410.1016/S0025-5564(98)10076-7MR1686468JinxiaoP.ZhenJ.ZhienM.Thresholds of survival for an n-dimensional Volterra mutualistic system in a polluted environment2000252251953110.1006/jmaa.2000.6853MR1800196LiZ.ShuaiZ.WangK.Persistence and extinction of single population in a polluted environment20042004108, 1C5MR2108879ZBL1084.92032HeJ.WangK.The survival analysis for a population in a polluted environment20091031555157110.1016/j.nonrwa.2008.01.027MR2502965ZBL1160.92041SinhaS.MisraO. P.DharJ.Modelling a predator-prey system with infected prey in polluted environment20103471861187210.1016/j.apm.2009.10.003MR2594573ZBL1193.34071LiuB.ChenL.ZhangY.The effects of impulsive toxicant input on a population in a polluted environment20031132652742-s2.0-014205960210.1142/S0218339003000907LiuB.TengZ.ChenL.The effects of impulsive toxicant input on two-species Lotka-Volterra competition system200512208220MR2183246ZBL1097.34040YangX.JinZ.XueY.Weak average persistence and extinction of a predator-prey system in a polluted environment with impulsive toxicant input200731372673510.1016/j.chaos.2005.10.042MR2262307ZBL1133.92032TaoF.LiuB.Dynamic behaviors of a single-species population model with birth pulses in a polluted environment20083851663168410.1216/RMJ-2008-38-5-1663MR2457381ZBL1178.34051LiuB.ZhangL.Dynamics of a two-species Lotka-Volterra competition system in a polluted environment with pulse toxicant input2009214115516210.1016/j.amc.2009.03.065MR2541054ZBL1172.92036CaiS.A stage-structured single species model with pulse input in a polluted environment200957337538210.1007/s11071-008-9448-xMR2520217ZBL1176.92051JiaoJ.ChenL.Dynamical analysis of a chemostat model with delayed response in growth and pulse input in polluted environment200946250251310.1007/s10910-008-9474-4MR2583046ZBL1196.92041JiaoJ.LongW.ChenL.A single stage-structured population model with mature individuals in a polluted environment and pulse input of environmental toxin20091053073308110.1016/j.nonrwa.2008.10.007MR2523269ZBL1162.92330MengX.LiZ.NietoJ. J.Dynamic analysis of Michaelis-Menten chemostat-type competition models with time delay and pulse in a polluted environment201047112314410.1007/s10910-009-9536-2MR2576641ZBL1194.92075JiaoJ.YeK.ChenL.Dynamical analysis of a five-dimensioned chemostat model with impulsive diffusion and pulse input environmental toxicant2011441–3172710.1016/j.chaos.2010.11.001MR2782490ZBL1216.92062JiaoJ.CaiS.ChenL.Dynamics of the genic mutational rate on a population system with birth pulse and impulsive input toxins in polluted environment2012401-244545710.1007/s12190-012-0577-5MR2965342JohnstonE. L.KeoughM. J.Field assessment of effects of timing and frequency of copper pulses on settlement of sessile marine invertebrates20001375-6101710292-s2.0-003449623310.1007/s002270000420MayR. M.2001Princeton University PressKuangY.1993Boston, Mass, USAAcademic Pressxii+398MR1218880BeddingtonJ. R.MayR. M.Harvesting natural populations in a randomly fluctuating environment197719743024634652-s2.0-0000735971RudnickiR.PichórK.Influence of stochastic perturbation on prey-predator systems2007206110811910.1016/j.mbs.2006.03.006MR2311674ZBL1124.92055HungL.-C.Stochastic delay population systems20098891303132010.1080/00036810903277093MR2574331ZBL1216.60051ZhuC.YinG.On hybrid competitive Lotka-Volterra ecosystems20097112e1370e137910.1016/j.na.2009.01.166MR2671923ZBL1238.34059LiuM.WangK.Persistence and extinction in stochastic non-autonomous logistic systems2011375244345710.1016/j.jmaa.2010.09.058MR2735535ZBL1214.34045LiuM.WangK.Stationary distribution, ergodicity and extinction of a stochastic generalized logistic system2012251119801985LiuM.WangK.WuQ.Survival analysis of stochastic competitive models in a polluted environment and stochastic competitive exclusion principle20117391969201210.1007/s11538-010-9569-5MR2823174ZBL1225.92059LiuM.WangK.Persistence and extinction of a single-species population
system in a polluted environment with random perturbations and
impulsive toxicant input20124515411550LiuM.QiuH.WangK.A remark on a stochastic predator-prey system with time delays201326331832310.1016/j.aml.2012.08.015MR3002141ZBL06130500LiuM.WangK.Population dynamical behavior of Lotka-Volterra cooperative systems with random perturbations20133362495252210.3934/dcds.2013.33.2495MR3007696ZBL06157314LiuM.WangK.Dynamics of a Leslie-Gower Holling-type II predator-prey system with Lévy jumps20138520421310.1016/j.na.2013.02.018MR3040360LiuM.WangK.Analysis of a stochastic autonomous mutualism model2013402139240310.1016/j.jmaa.2012.11.043MR3023266ZBL06156131LiuM.WangK.Stochastic Lotka—Volterra systems with Lévy noise2014410750763IkedaN.WatanabeS.A comparison theorem for solutions of stochastic differential equations and its applications1977143619633MR0471082ZBL0376.60065KloedenP. E.ShardlowT.The Milstein scheme for stochastic delay differential equations without using anticipative calculus201230218120210.1080/07362994.2012.628907MR2891451ZBL1247.65005