Survival analysis of a stochastic delay single-species system in polluted environment with psychological effect and pulse toxicant input

*Correspondence: daiaga0921@126.com 1School of Data Science of TongRen University, 554300 Tongren, P.R. China Full list of author information is available at the end of the article Abstract We propose and study a stochastic delay single-species population system in polluted environment with psychological effect and pulse toxicant input. We establish sufficient conditions for the extinction, nonpersistence in the mean, weak persistence, and strong persistence of the single-species population and obtain the threshold value between extinction and weak persistence. Finally, we confirm the efficiency of the main results by numerical simulations.


Introduction
Along with fast development of agriculture and modern industry, a large number of toxic pollutants are discharged into the ecosystem, and therefore it is an undeniable fact that environmental pollution becomes increasingly serious, such as pollution of pollutants from burning agricultural plant straw, heavy metal pollution, water pollution caused by crop fertilization and pesticide application. As it is well known, the existence of various poisons are becoming a threat to the survival of unprotected populations, which has prompted many scholars to investigate the impact of toxins on the population and assess the risk of the population. An important tool to analyze the effects of toxins on population is establishing a mathematical model [1][2][3][4][5][6][7].
Wei and Chen [8] proposed a mathematical model for the first time to study the physiological effects of vertebrates on population in a polluted environment. The so-called "psychological" effect refers to: in the heavy pollution environment, because the organism with spine has a good sensory nervous system, which can transmit the information of the polluted environment to the area where the brain can explain the information, the effective contact between the organism and the live environment will be reduced, which plays the role of self-protection; for instance, fish can identify the information in the polluted environment through their own neurosensory system and make decisions to either escape from the polluted area or bear the environmental toxicant [8]. Afterwards, Lan and Wei [9] considered that the exogenous input of toxins is regular in some practical situations. They proposed the following single-species population model with psychological effect and impulsive toxicant in a pollution environment: 1+αc 2 e (t) ), c 0 (t) = (kc e (t) -(g + m + b)c 0 (t)), c e (t) = -hc e (t), t = nγ , n ∈ z + , where x(t) represents the density of the population, c 0 (t) and c e (t) represent the concentration of toxins in the organism and the concentration of toxins in the environment, respectively, r and r 0 stand for the net growth rate in nonpolluted environment and response intensity of biological growth to toxins, kc e (t) stands for the uptake of toxins in the environment, gc 0 (t) represents the emission rate of toxins, mc 0 (t) represents the purification rate of toxins because of metabolic process of organisms, bc 0 (t) is the loss due to giving birth at time t, -hc e (t) denotes the amount of reduction in the purification of toxins by the environment itself, μ and γ denote the toxicant input amount and the period of pulse input toxin, respectively, and In reality, population growth is more or less disturbed by environmental factors, such as temperature, humidity, and seasonal climate change, and almost all the observed data show that there are obvious random fluctuations in the growth process of organisms. Therefore, in some practical cases, ignoring the randomness of the system and using deterministic models to describe and predict the system behavior are not always satisfactory; especially, it is not suitable to use deterministic population model to study how to protect endangered species [10]. May [11] pointed out that due to the impact of environmental noises, the birth rate, death rate, carrying capacity, competition coefficients, and other parameters involved with the system exhibit random fluctuation to a greater or lesser extent [12][13][14][15][16][17]. On the other hand, the influence of delay has played an important role in the population dynamic [18][19][20][21][22].
In this paper, we suppose that the environmental noises affect all parameters of model (1) (see e.g. [23,24]). Considering the work of population psychological effect, we propose the following stochastic impulsive single-species population model with delay and psychological effect in polluted environment: 1+αc 2 e (t) dB 5 (t), dc 0 (t) = (kc e (t) -(g + m + b)c 0 (t)) dt, dc e (t) = -hc e (t) dt, t = nγ , n ∈ z + , where τ represents the time delay, B i (t) stands for a standard Brownian motion defined on a complete probability space (Ω, F, P) with filtration {F t } t∈R + , and σ 2 i (i = 1, 2, 3, 4, 5) represents the intensity of noise. Let ϕ(θ ) be a continuous function on [-τ , 0], and let the solution x(t) of system (2) satisfy the initial condition Remark 1.1 (see [7]) In model (2), since c 0 (t) and c e (t) denote the concentrations of toxicant, we must have 0 ≤ c 0 (t) ≤ 1 and 0 ≤ c e (t) ≤ 1 for all t ≥ 0. To this end, we need the following constraints: k ≤ g + m and b ≤ 1e -hγ .

Lemma 2.2
The positive γ -periodic solution (c o (t),c e (t)) has the following properties: Proof By the lemma of Yang [26] and Lemma 2.1 of Lan [9] we get Next, we will prove the last two limits. From , Similarly, In view of the periodicity ofc 0 (t) andc e (t) we can observe that lim t→+∞ x(t) = 0 a.s., λ < 0.
(2) If there are positive constants λ, λ 0 , and T such that Lemma 2.4 (see [14]) Consider the stochastic differential equation where r, a, and σ are positive constants. If r -0.5σ 2 > 0, then Proof Considering the stochastic differential equation

Main results
The proof of the existence and uniqueness of the global positive solution of the equation is similar to that of Dai et al. [24] by defining the nonnegative function Hence we omit it.
Next, we will focus on proving that lim t→+∞ Using Itô's formula again, we get Integrating both sides of equality (6) from 0 to t, we get where M 1 (t) = Using the exponential martingale inequality, for all δ > 0, β > 0, T > 0, and θ > 1, we have Taking δ = e -k , β = θ e k ln k, and T = k, we have By the Borel-Cantelli lemma there are an event Ω and positive integers k 1 = k 1 (ω) such that P(Ω) = 1 and for all ω ∈ Ω and k > k 1 , we have Similarly, By (7) we get Because a > 0, we can easily obtain that there is a positive constant K > 0 such that ln x + rax ≤ K , and from (8) it follows that e t ln x(t) ≤ ln x(0) + K e t -1 + 5θ e k ln k, 0≤ t ≤ k.

Theorem 3.2 Let x(t) be the solution of system (2).
(i) If A < 0, then the population x(t) will die out almost surely.

(ii) If A = 0, then the population x(t) is nonpersistent in the mean almost surely. (iii) If A > 0, then the population x(t) is weakly persistent almost surely.
Proof Integrating from 0 to t both sides of equation (5), we get where The quadratic variations of N 2 (t), N 3 (t), N 4 (t), and N 5 (t) are For all t ∈ (nγ , (n + 1)γ ], n ∈ N + , we have Noting that 0 ≤ c 0 (t) ≤ 1 and 0 ≤ c e (t) ≤ 1, by Lemmas 2.1 and 2.2 we get Similarly, By the strong law of large numbers it follows that By the exponential martingale inequality, choosing δ = 1, β = 2 ln k, and T = k, have It follows from Borel-Cantelli lemma that there exists a positive constant k 1 such that, for k > k 1 , we have Substituting this inequality into (9), we obtain that For all t ∈ (nγ , (n + 1)γ ], we have ds.

Theorem 3.3 Suppose that σ 3 = σ 4 = 0, and let x(t) be the positive solution of model (2)
with initial value x(θ ) = ϕ(θ ) ∈ C([-τ , 0], R + ). Then: (I) If A < 0, then the population x goes to die out a.s., that is, lim t→+∞ x(t) = 0 a.s. (II) If A > 0, then the population x is strongly persistent in the mean a.s.; moreover, Proof If σ 3 = σ 4 = 0, then by (19) we have It follows from (12) and (16) that for all > 0, there exists a positive constant T 1 such that for t > T 1 , we have (I) If A < 0, then by Lemma 2.3 and (21) we have lim t→+∞ x(t) = 0 a.s. Now we will prove (II). Let us consider the following auxiliary equation: dt + σ 1 y(t) dB 1 (t) + σ 2 c 0 (t)y(t) dB 2 (t) + σ 5 y(t)c e (t) 1 + αc 2 e (t) where It follows from the stochastic comparison theorem [14] that x(t) ≤ y(t). Hence, for all > 0, there is a positive constant T 2 such that for all t > T 2 , we have By (24), (25), and Lemma 2.3 this implies that Due to the arbitrariness of , from (26) we get Since x(t) ≤ y(t), by (27) we easily see that It follows from (28) that, for all > 0, there exists a positive constant T 3 such that for t > T 3 , have For all > 0, there exists T = max{T 1 , T 3 } such that, for t > T, have By (29) and (30) (2) for the single-species population extinction and weak persistence. If A < 0, then the population will be extinct, and if A > 0, then the population is weakly persistent. Particularly, when σ 3 = σ 4 = 0, A is also the threshold of system (2) for the single-species population extinction and strong persistence; moreover, if A > 0, then the single-species population is stable in the mean.
Remark 3.5 Note that from the expression for A = r -r 0 kμ h(g+m+b)γλη -0.5σ 2 1 -0.5σ 2 2 ν -0.5σ 2 5 ξ we can find that the parameters μ and γ obviously affect the persistence and extinction of system (2), that is, we can control the persistence and extinction of population x by controlling the toxicant input amount μ and the period of pulse input toxicant γ . Remark 3.6 Theorems 3.2 and 3.3 show that the persistence population x of deterministic system (1) may be extinct when σ 1 , σ 2 , and σ 5 are large enough; however, τ , σ 3 , and σ 4 have no effect on the persistence and extinction for population x.

Numerical simulations and discussion
Next, we show the numerical simulation results to illustrate the accuracy of analytical results in the previous section by using the famous Milstein method [28]. We choose the parameters of system (2) as follows: To understand the effects of white noise and pulse toxicant input on population dynamics, we change the values of μ, γ , and σ i (i = 1, 2, 3, 4, 5).
(2) Next, to analyze the effect of the toxicant input amount each time μ on the persistence of the single species, we adopt u = 0.8, σ 1 = 0.3, σ 2 = 0.2, σ 3 = 0.2, σ 4 = 0.4, σ 5 = 0.4, and γ = 5. Simple calculation shows that A = -0.0555 < 0, and from Theorem 3.2 it follows that the population x will die out a.s. (see Fig. 3). From Figs. 1 and 3 we can see that the population x will die out when the environmental toxicant amount of each time μ is large enough. (3) On the other hand, we will focus on the influence of the intensity of white noises on the survival for the population x. We adopt u = 0.3 and γ = 5 and suppose that σ 1 = 0.8, σ 2 = 0.2, σ 3 = 0.2, σ 4 = 0.6, and σ 5 = 0.6. Simple calculation shows that A = -0.0823 < 0, so that the population x will die out (see Fig. 4). Suppose σ 1 = 0.3, σ 2 = 0.2, σ 3 = 0, σ 4 = 0, and σ 5 = 0.4. Then A = 0.1955 > 0, and by Theorem 3.3 we obtain that the population x is strongly persistent in the mean; moreover, lim t→+∞ t -1 t 0 x(s) ds = 0.2443, that is, the population is stable in the mean (see Fig. 5). From Figs. 1 and 5 we can see that the population x will die out when σ 1 , σ 2 , σ 5 are large enough, but σ 3 and σ 4 have no effect on the survival of the population x.

Conclusions
We studied a stochastic impulsive single-species population model with psychological effect and delay in pollution environment. We obtained the threshold A between weak persistence and extinction. Particularly, A is also a threshold of system (2) for the strong persistence and extinction when σ 3 = σ 4 = 0. These results have revealed that the environmental noise, the impulsive period, and the amount of toxicant input for each time have