## Nonlinear Dynamics in Applied Sciences Systems: Advances and Perspectives

View this Special IssueResearch Article | Open Access

# Analysis of Stochastic Delay Predator-Prey System with Impulsive Toxicant Input in Polluted Environments

**Academic Editor:**Massimiliano Ferrara

#### Abstract

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: where all the parameters are positive constants and, and : the size of prey population and the predator population, respectively; : the intrinsic growth rate of the th population without toxicant; : the th population response to the pollutant present in the organism; : the concentration of toxicant in the th organism; : the concentration of toxicant in the environment; : the organism’s net uptake of toxicant from the environment; : the egestion and depuration rates of the toxicant in the th organism; : the toxicant loss from the environment itself by volatilization and so on; : the period of the impulsive effect about the exogenous input of toxicant; : the toxicant input amount at every time.

Yang et al. [15] showed that in the following Lemma holds.

Lemma 1. *For system (1), define
*(a)*If , then , .*(b)*If and , then and goes to extinction.*(c)*If , then ,.*

Some interesting and important problems arise naturally. 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? 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? 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 (see, e.g., [27–39]), where is a white noise and is the intensity of the noise. Moreover, taking time delays into account, we obtain the following model: with initial condition where , , is continuous on . Our main result is the following theorem.

Theorem 2. *For system (3), define
*(i)*If , then both and go to extinction almost surely (a.s.); that is, .*(ii)*If and , then goes to extinction and is stable in time average a.s.; that is,
*(iii)*If , then both and are stable in time average a.s.
*

*Remark 3. *By comparing Lemma 1 with our Theorem 2, we can see that on the one hand, if and , then , , , 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:

Lemma 4. *For any given initial value , there is a unique global positive solution to the first two equations of system (3) a.s.*

*Proof. *The proof is similar to Hung [29] by defining
and hence is omitted.

To begin with, let us consider the following subsystem of (3):

Lemma 5 (see [13, 15]). *System (10) has a unique positive -periodic solution , and for each solution of (10), , , and as . Moreover, and for all if and , , where
**
for and . In addition,
*

Lemma 6 (see [34]). * Suppose that . * *(I) If there exist and positive constants , such that
* *for , where are independent standard Brownian motions and are constants, , then one has the following: if , then a.s.; if , then .* *(II) If there exist positive constants , and such that
* *for , then a.s.*

Now, let us consider the following auxiliary system: with initial value .

Lemma 7. * If , then the solution of system (15) obeys
*

*Proof. *By Lemma 5,
Then, for all, there exists such that
An application of Itô’s formula to (15) yields
That is to say, we have shown that
When (18) is used in (20), we can see that for ,
Let be sufficiently small such that . Making use of (I) and (II) in Lemma 6 to (22) and (23), respectively, we have
It then follows from the arbitrariness of that
Substituting (17) and (25) into (20) and noting that , one can derive that
Employing (20) and (21) in the expression yields
In view of (25), we get
By (17), (26), (27), and (28), for all, there exists such that, for ,
If , then we can choose sufficiently small such that . Then, by (29) and (I) in Lemma 6, we obtain . If , then we can choose sufficiently small such that . An application of (I) and (II) in Lemma 6 to (29) and (30), respectively, makes one observe that
Therefore, using the arbitrariness of results in
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

It follows from (17) and (33) that
for sufficiently large . Since , then we can choose sufficiently small such that . Then, by (I) in Lemma 6,
When (36) is used in (34), one can see that
for sufficiently large , where obeys . In view of Lemma 6 again, ,.

By the stochastic comparison theorem [40], one can observe that
Note that and ; it then follows from Lemma 7 that ,. Making use of (38) gives ,. Thus, for all, there exists such that, for ,
Substituting the above inequalities into (33) and then using (18), we obtain
Let be sufficiently small such that , and then, applying (I) and (II) in Lemma 6 to (40) and (41), respectively, one can see that
An application of the arbitrariness of gives

Clearly, implies , and then, by Lemma 7,
Thus, similar to the proof of (28), we get
Therefore, by (26), (28), and (38), we can observe that
Employing (33) and (34) in the expression yields
When (18), (46) and (47), are used in (48), one can obtain
for sufficiently large , where obeys . It then follows from in Lemma 6 that
By virtue of the arbitrariness of , we can see that
Consequently, for every , there is such that
Substituting the above inequality into (33) and then using (18) and (47), one can see that
for sufficiently large . Since , and then, by Lemma 6 and the arbitrariness of , one can observe that
When this inequality, (18) and (47), are used in (34), we can see that
for sufficiently large . Then, it follows from Lemma 6 and the arbitrariness of that
Substituting the above inequality and (18) into (33), we get
for sufficiently large . By in Lemma 6 and the arbitrariness of again, we obtain
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 , , , , , , , , , , , , , and . Then, By in Lemma 1, the solution of model (1) obeys 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 is different. In Figure 1(a), we choose . Therefore, Then, by in Theorem 2, both and are extinctive. Figure 1(a) confirms these. In Figure 1(b), we choose . That is to say and . It then follows from in Theorem 2 that is extinctive and is stable in time average: See Figure 1(b). In Figure 1(c), we choose . Then, . In view of in Theorem 2, we can obtain that both and are stable in time average: Figure 1(c) confirms these.

**(a)**

**(b)**

**(c)**

In Figure 2, we choose , , , , , , , , , , , , , , and . The only difference between conditions of Figures 1(c) and 2 is that the value of is different. In Figure 2, we choose . Then, . It follows from in Theorem 2 that both and 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.

#### 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. (A)Firstly, time delay is harmless for stability in time average and extinction of the stochastic system (3).(B)The white noise and can change the properties of the system greatly.(C)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 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).

#### References

- T. G. Hallam, C. E. Clark, and R. R. Lassiter, “Effects of toxicants on populations: a qualitative approach—1. Equilibrium environmental exposure,”
*Ecological Modelling*, vol. 18, no. 4, pp. 291–304, 1983. View at: Google Scholar - T. G. Hallam, C. E. Clark, and G. S. Jordan, “Effects of toxicants on populations: a qualitative approach—2. First order kinetics,”
*Journal of Mathematical Biology*, vol. 18, no. 1, pp. 25–37, 1983. View at: Google Scholar - T. G. Hallam and J. T. De Luna, “Effects of toxicants on populations: a qualitative approach—3. Environmental and food chain pathways,”
*Journal of Theoretical Biology*, vol. 109, no. 3, pp. 411–429, 1984. View at: Google Scholar - Z. E. Ma, G. R. Cui, and W. D. Wang, “Persistence and extinction of a population in a polluted environment,”
*Mathematical Biosciences*, vol. 101, no. 1, pp. 75–97, 1990. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - H. I. Freedman and J. B. Shukla, “Models for the effect of toxicant in single-species and predator-prey systems,”
*Journal of Mathematical Biology*, vol. 30, no. 1, pp. 15–30, 1991. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - L. Huaping and M. Zhien, “The threshold of survival for system of two species in a polluted environment,”
*Journal of Mathematical Biology*, vol. 30, no. 1, pp. 49–61, 1991. View at: Publisher Site | Google Scholar | MathSciNet - M. Zhien, Z. Wengang, and L. Zhixue, “The thresholds of survival for an
*n*-food chain model in a polluted environment,”*Journal of Mathematical Analysis and Applications*, vol. 210, no. 2, pp. 440–458, 1997. View at: Publisher Site | Google Scholar | MathSciNet - B. Buonomo, A. Di Liddo, and I. Sgura, “A diffusive-convective model for the dynamics of population-toxicant interactions: some analytical and numerical results,”
*Mathematical Biosciences*, vol. 157, no. 1-2, pp. 37–64, 1999. View at: Publisher Site | Google Scholar | MathSciNet - P. Jinxiao, J. Zhen, and M. Zhien, “Thresholds of survival for an
*n*-dimensional Volterra mutualistic system in a polluted environment,”*Journal of Mathematical Analysis and Applications*, vol. 252, no. 2, pp. 519–531, 2000. View at: Publisher Site | Google Scholar | MathSciNet - Z. Li, Z. Shuai, and K. Wang, “Persistence and extinction of single population in a polluted environment,”
*Electronic Journal of Differential Equations*, vol. 2004, no. 108, 1C5, 2004. View at: Google Scholar | Zentralblatt MATH | MathSciNet - J. He and K. Wang, “The survival analysis for a population in a polluted environment,”
*Nonlinear Analysis. Real World Applications*, vol. 10, no. 3, pp. 1555–1571, 2009. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - S. Sinha, O. P. Misra, and J. Dhar, “Modelling a predator-prey system with infected prey in polluted environment,”
*Applied Mathematical Modelling*, vol. 34, no. 7, pp. 1861–1872, 2010. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - B. Liu, L. Chen, and Y. Zhang, “The effects of impulsive toxicant input on a population in a polluted environment,”
*Journal of Biological Systems*, vol. 11, no. 3, pp. 265–274, 2003. View at: Publisher Site | Google Scholar - B. Liu, Z. Teng, and L. Chen, “The effects of impulsive toxicant input on two-species Lotka-Volterra competition system,”
*International Journal of Information & Systems Sciences*, vol. 1, no. 2, pp. 208–220, 2005. View at: Google Scholar | Zentralblatt MATH | MathSciNet - X. Yang, Z. Jin, and Y. Xue, “Weak average persistence and extinction of a predator-prey system in a polluted environment with impulsive toxicant input,”
*Chaos, Solitons & Fractals*, vol. 31, no. 3, pp. 726–735, 2007. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - F. Tao and B. Liu, “Dynamic behaviors of a single-species population model with birth pulses in a polluted environment,”
*The Rocky Mountain Journal of Mathematics*, vol. 38, no. 5, pp. 1663–1684, 2008. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - B. Liu and L. Zhang, “Dynamics of a two-species Lotka-Volterra competition system in a polluted environment with pulse toxicant input,”
*Applied Mathematics and Computation*, vol. 214, no. 1, pp. 155–162, 2009. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - S. Cai, “A stage-structured single species model with pulse input in a polluted environment,”
*Nonlinear Dynamics*, vol. 57, no. 3, pp. 375–382, 2009. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - J. Jiao and L. Chen, “Dynamical analysis of a chemostat model with delayed response in growth and pulse input in polluted environment,”
*Journal of Mathematical Chemistry*, vol. 46, no. 2, pp. 502–513, 2009. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - J. Jiao, W. Long, and L. Chen, “A single stage-structured population model with mature individuals in a polluted environment and pulse input of environmental toxin,”
*Nonlinear Analysis. Real World Applications*, vol. 10, no. 5, pp. 3073–3081, 2009. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - X. Meng, Z. Li, and J. J. Nieto, “Dynamic analysis of Michaelis-Menten chemostat-type competition models with time delay and pulse in a polluted environment,”
*Journal of Mathematical Chemistry*, vol. 47, no. 1, pp. 123–144, 2010. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - J. Jiao, K. Ye, and L. Chen, “Dynamical analysis of a five-dimensioned chemostat model with impulsive diffusion and pulse input environmental toxicant,”
*Chaos, Solitons & Fractals*, vol. 44, no. 1–3, pp. 17–27, 2011. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - J. Jiao, S. Cai, and L. Chen, “Dynamics of the genic mutational rate on a population system with birth pulse and impulsive input toxins in polluted environment,”
*Journal of Applied Mathematics and Computing*, vol. 40, no. 1-2, pp. 445–457, 2012. View at: Publisher Site | Google Scholar | MathSciNet - E. L. Johnston and M. J. Keough, “Field assessment of effects of timing and frequency of copper pulses on settlement of sessile marine invertebrates,”
*Marine Biology*, vol. 137, no. 5-6, pp. 1017–1029, 2000. View at: Publisher Site | Google Scholar - R. M. May,
*Stability and Complexity in Model Ecosystems*, Princeton University Press, 2001. - Y. Kuang,
*Delay Differential Equations with Applications in Population Dynamics*, Academic Press, Boston, Mass, USA, 1993. View at: MathSciNet - J. R. Beddington and R. M. May, “Harvesting natural populations in a randomly fluctuating environment,”
*Science*, vol. 197, no. 4302, pp. 463–465, 1977. View at: Google Scholar - R. Rudnicki and K. Pichór, “Influence of stochastic perturbation on prey-predator systems,”
*Mathematical Biosciences*, vol. 206, no. 1, pp. 108–119, 2007. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - L.-C. Hung, “Stochastic delay population systems,”
*Applicable Analysis*, vol. 88, no. 9, pp. 1303–1320, 2009. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - C. Zhu and G. Yin, “On hybrid competitive Lotka-Volterra ecosystems,”
*Nonlinear Analysis. Theory, Methods & Applications*, vol. 71, no. 12, pp. e1370–e1379, 2009. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - M. Liu and K. Wang, “Persistence and extinction in stochastic non-autonomous logistic systems,”
*Journal of Mathematical Analysis and Applications*, vol. 375, no. 2, pp. 443–457, 2011. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - M. Liu and K. Wang, “Stationary distribution, ergodicity and extinction of a stochastic generalized logistic system,”
*Applied Mathematics Letters*, vol. 25, no. 11, pp. 1980–1985, 2012. View at: Google Scholar - M. Liu, K. Wang, and Q. Wu, “Survival analysis of stochastic competitive models in a polluted environment and stochastic competitive exclusion principle,”
*Bulletin of Mathematical Biology*, vol. 73, no. 9, pp. 1969–2012, 2011. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - M. Liu and K. Wang, “Persistence and extinction of a single-species population
system in a polluted environment with random perturbations and
impulsive toxicant input,”
*Chaos Solitons Fractals*, vol. 45, pp. 1541–1550, 2012. View at: Google Scholar - M. Liu, H. Qiu, and K. Wang, “A remark on a stochastic predator-prey system with time delays,”
*Applied Mathematics Letters*, vol. 26, no. 3, pp. 318–323, 2013. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - M. Liu and K. Wang, “Population dynamical behavior of Lotka-Volterra cooperative systems with random perturbations,”
*Discrete and Continuous Dynamical Systems A*, vol. 33, no. 6, pp. 2495–2522, 2013. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - M. Liu and K. Wang, “Dynamics of a Leslie-Gower Holling-type II predator-prey system with Lévy jumps,”
*Nonlinear Analysis. Theory, Methods & Applications*, vol. 85, pp. 204–213, 2013. View at: Publisher Site | Google Scholar | MathSciNet - M. Liu and K. Wang, “Analysis of a stochastic autonomous mutualism model,”
*Journal of Mathematical Analysis and Applications*, vol. 402, no. 1, pp. 392–403, 2013. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet - M. Liu and K. Wang, “Stochastic Lotka—Volterra systems with Lévy noise,”
*Journal of Mathematical Analysis and Applications*, vol. 410, pp. 750–763, 2014. View at: Google Scholar - N. Ikeda and S. Watanabe, “A comparison theorem for solutions of stochastic differential equations and its applications,”
*Osaka Journal of Mathematics*, vol. 14, no. 3, pp. 619–633, 1977. View at: Google Scholar | Zentralblatt MATH | MathSciNet - P. E. Kloeden and T. Shardlow, “The Milstein scheme for stochastic delay differential equations without using anticipative calculus,”
*Stochastic Analysis and Applications*, vol. 30, no. 2, pp. 181–202, 2012. View at: Publisher Site | Google Scholar | Zentralblatt MATH | MathSciNet

#### Copyright

Copyright © 2013 Meng Liu. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.