Abstract

A stochastic prey-predator system in a polluted environment with Beddington-DeAngelis functional response is proposed and analyzed. Firstly, for the system with white noise perturbation, by analyzing the limit system, the existence of boundary periodic solutions and positive periodic solutions is proved and the sufficient conditions for the existence of boundary periodic solutions and positive periodic solutions are derived. And then for the stochastic system, by introducing Markov regime switching, the sufficient conditions for extinction or persistence of such system are obtained. Furthermore, we proved that the system is ergodic and has a stationary distribution when the concentration of toxicant is a positive constant. Finally, two examples with numerical simulations are carried out in order to illustrate the theoretical results.

1. Introduction and Model Formulation

The Lotka-Volterra model [13] is a classical model in the study of biological mathematics, and the continuous Lotka-Volterra model which is modeled by ordinary differential equations and delay differential equations is widely used to characterize the dynamics of biological systems [413]. The functional response functions are important in the population ecological models [14]. In general, functional responses fall into two categories: one depends only on the density of the prey, such as Holling I–III [1517]; the other depends on the density of both the prey and the predator, such as Beddington-DeAngelis type [18, 19]. Compared with the Holling II functional response, the Beddington-DeAngelis type functional response, , has an additional term in the denominator modeling mutual interference among predators. In other words, this type of functional response is affected by both predator and prey. Some biologists believe that if the predators compete with each other to obtain food, functional response should depend on the density of both the prey and the predator. Arditi et al. [20] and Jost et al. [21, 22] used the actual observation data to verify this point. In particular, having collected observation data from 19 predator-prey communities, Skalski and Gilliam [23] found that predator-dependent functional responses were in agreement with the observation data, and in many instances, the Beddington-DeAngelis type looked better than the others. The Beddington-DeAngelis functional response has been widely used in the modeling of ecosystems in which there is mutual interference among predators [24, 25]. In [19], DeAngelis et al. have extensively investigated the dynamical properties of the following prey-predator system:where and represent the density of the prey and the predator, respectively. is the intrinsic growth rate of the prey, , , and are the consumption rate, the saturation constant, and the saturation constant for an alternative prey, respectively. is the conversion rate of nutrients into the reproduction for the predator. The parameters and are the nonpredatory loss rate and the interspecific competition rate. We refer the reader to [19] for more details.

In many ecosystems, predators tend to be omnivorous, they have wide variety of food sources. For example, the giant panda is omnivorous animal, since it can eat both meat and plant such as bamboos. In the lake ecosystem, some fishes not only prey on aquatic invertebrates, but also feed on algae and other aquatic plants. Polis and Strong in [26] and McCann and Hastings in [27] studied omnivorous nature of animals in the food chain in 1996 and 1997, respectively. Based on the above literature, we established a kind of omnivorous model as follows:where represents the growth rate of due to omnivorous nature and denote the density-dependent coefficient of the prey and the predator, respectively. , , , , and are the consumption rate, the saturation constant, the predator interference, the saturation constant for an alternative prey, and the conversion rate, respectively. All parameters are positive in system (2).

It is well known that the biological population is inevitably affected by environment perturbation while the stochastic population model is more in line with the actual situation. Recently, various models based on stochastic differential equations (SDEs) have extensively been paid the attention of the researchers (see, e.g., [2837]). Parameter perturbation induced by white noise is an important and common form to describe the effect of stochasticity (see, e.g., [3748]). In this paper, we consider the white noise perturbation for the intrinsic growth rates of the prey and predator; that is, and , where , are mutually independent Brownian motions and , denote the intensities of the white noise. On the other hand, it can be seen from the recent literature that the environmental pollution has an important effect on the population systems [4960]. In 1983, Hallam et al. [61, 62] studied the influence of environmental pollution on the population and established a relationship model between environmental toxins and population. Subsequently, Hallam et al. [63, 64] studied the persistence and extinction of population in polluted environment. The mathematical model established by Hallam et al. considered only the toxins in the organism to cause a decrease in the birth rate or an individual death, which is reasonable in the case of lower concentration of the toxicant in the environment. When pollution is serious, the emission of pollutants may directly lead to the death of the species; see [6569]. The authors in [68] added the environmental toxic term directly to the model; this is reasonable in the heavily polluted environment. For example, in a lake ecosystem, the discharge of large amounts of industrial waste water may directly lead to the death of fish, aquatic invertebrates, and so on. Therefore, we assume that the emission of pollutants to the environment is impulsive and directly affects the survival of the species in such an environment, so we get the following system:where , are positive, nonconstant, and continuous functions of period , stands for the concentration of the toxicant in the environment, denotes the loss rate of toxicant at time , is the impulsive input period and is the impulsive input amount, and and represent the dose-response of the prey and predator to the environmental toxicant, respectively.

Furthermore, the prey-predator model may be perturbed by telegraph noise which is distinguished by factors such as rain falls and nutrition and can be represented by switching among two or more regimes of environment [40, 60, 7080]. For example, population growth rates in different seasons are not the same. The intraspecific competition coefficient varies according to the changes in nutrition and food resources. Generally, the switching between different regimes is memoryless and the waiting time for the next switch is exponentially distributed [81, 82]. Therefore, it can be described by a continuous-time Markov chain taking values in a finite state space . Taking into account the influences of white noise and telegraph noise, we propose the following stochastic differential system with impulsive toxicant input:For any , , , , , , , , , and are all positive constants. In model (4), the population is inevitably affected by severe stochastic interference such as drought; the parameter switches one state into another state and it will switch into the next regime until the next major environmental change.

The rest of this paper is organized as follows. In Section 2, we provide preliminaries which are used in the following sections. In Section 3, we show that system (3) admits a nontrivial positive -periodic solution by constructing Lyapunov function. In Section 4, we explore the sufficient conditions for extinction and permanence in mean of system (4). Finally, some examples with numerical simulations have been given to illustrate our theoretical results.

2. Preliminaries

Throughout this paper, let be a complete probability space with a filtration satisfying the usual conditions, is one-dimensional Brownian motion on this space, and is a right-continuous Markov chain and independent of the Brownian motion . The state space of this Markov chain is . Suppose that the generator matrix of is , where stands for the transition rate from state to and satisfies the following conditions: here, if , while , As a standing hypothesis, we assume that the Markov chain is irreducible, which means that system (4) can switch from one regime to another. Under this assumption, the Markov chain has a unique stationary distribution which is the solution of the system of linear equations subject to and for all . Hence, for any vector , we have that Let us consider the following stochastic differential equation with Markov conversion.where , , , and is a -dimensional Brownian motion defined on the underlying probability space. The matrix is called the diffusion matrix. Let be twice continuously differentiable and which is defined as follows be the diffusion operator about : Particularly, for one-dimensional stochastic systemthe following two lemmas can be given from referring to the articles [72, 77].

Lemma 1. System (10) has a unique continuous positive solution . When it exists, the solution is global and stochastically ultimately bounded.

Lemma 2. Suppose that ; then
(i) system (10) is stochastic permanent if and only if ;
(ii) system (10) is extinct if and only if ;
(iii) when , system (10) is ergodic and there exists a unique stationary distribution , such that

Next, we consider the following stochastic differential equation:

Lemma 3 (see [78]). Suppose that the coefficients of (12) are -periodic in and there exists a function which is -periodic in , and satisfies the following conditions:
(i) as .
(ii) outside some compact set.
Then there exists a solution for (12) which is a -periodic Markov process.

Furthermore, we introduce some results from [80, 83] in Lemmas 4 and 5, which will be used in next section.

Lemma 4 (see [80]). Let . Then
(i) if there are two positive constants and such that holds for all and constants , then (ii) if there are three positive constants , , and such that holds for any , then

Finally, we give some basic properties of the following subsystem of system (3),

Lemma 5 (see [83]). System (17) has a unique -periodic solution which is globally asymptotically stable. Here , , , and .

For convenience and simplicity, define , , and , where is an integrable function on . If is a bounded function on , define .

3. Existence of Periodic Solutions of System (3)

In this section, we devote our attention to the investigation of the existence of periodic solutions of system (3). From Lemma 5, we know that system (17) has a globally asymptotically stable periodic solution ; therefore, the limit system of (3) iswhere , , and are all positive and continuous functions of period .

Now, we discuss the existence of periodic solutions of system (18).

Define Then, we have the following theorem about periodic solutions of system (18).

Theorem 6. If and , there exists a prey extinction periodic solution of system (18).

Proof. From the first equation of system (18), it is easy to see Applying Itô’s formula and then integrating from to , we obtainwhere is local martingale. From strong law of large numbers for martingales (see [84]), we have It then follows from (21) by dividing on both sides and letting that namely, tends to zero exponentially almost surely.
Since , a.s., from the second equation of system (18), its limit system isAccording to Theorem   in [85], when and (24) has a unique positive -periodic solution .
Overall, when and , there exists a prey extinction periodic solution of system (18).
The proof of this theorem is completed.

In order to investigate the existence of a nontrivial positive -periodic solution for system (18), first of all, we assume following conditions hold., and and

Theorem 7. Suppose that , , and hold, then there exists a positive -periodic solution for system (18).

Proof. Obviously, the coefficients of system (18) are continuous bounded positive periodic functions in . Now, we show that conditions (i) and (ii) of Lemma 3 hold. Define a nonnegative -function where , , , , and is a function defined on satisfying and . Obviously, is a -periodic function on . Therefore, the function is -periodic in and satisfies where . Therefore, condition (i) of Lemma 3 holds. Next, we will prove that condition (ii) of Lemma 3 also holds.
Applying Itô’s formula, one has Therefore,Define a bounded closed set where is a sufficiently small number such thatand , are quantities to be determined in the rest of the proof.
Denote Note that . Now, we prove .
Case 1. If , from (29), it implies that where . In fact, from condition , one can get that is to say, .
Case 2. If , from (29) and (32), we can get where . Using condition , one can get .
Case 3. If , then where By (33), we have .
Case 4. If , then where By (34), we obtain
Thus, Therefore, the proof of Theorem 7 is completed.

4. Extinction and Persistence in Mean of System (4)

In this section, we investigate the long-term dynamic behaviors of the prey-predator system (4) with white noise and telegraph noise in a polluted environment and then discuss the extinction and average persistence of prey and predator. According to Lemma 5, the periodic solution of the toxicant input is globally asymptotically stable, so the limit system of (4) isIn order to obtain the threshold conditions of persistence and extinction of system (44), we assume that,,,,.

Theorem 8. Given initial value for system (44), then
(i) if is established, the prey population will be extinct,
(ii) if is established, the predator population will be extinct,
(iii) if and are established, both the prey and the predator will die out.

Proof. (i) By Itô’s formula, we get then, By the ergodic theory of the Markov chain and the strong law of large number, we have from , we know (ii) Similarly, from the second equation of system (44), we have then, further, from , we know (iii) By the condition , one can get a.s., so that the limit system of the second equation of system (44) is Through Lemma 2, if holds, we obtain This completes the proof of the theorem.

Remark 9. If remains unchanged and or increases so that or , then condition or is established. That is to say, if the intrinsic growth rate and the predation intensity are relatively fixation, the increase of white noise intensity or pollutant concentration will lead to the extinction of the biological population.

Next, we will discuss the persistence of system (44). Applying Itô’s formula to the first equation of system (44), one can get then, when is large enough, we derive where is a sufficiently small positive number. In view of and Lemma 4, we deduce That is, the prey population of system (44) will be persistence in mean under condition .

Furthermore, the persistent property of the predator species of system (44) can be investigated as follows.

From the first equation of system (44), we have Consider the following stochastic differential equations:Obviously, , . Using Lemma 2, if , system (61) is ergodic and there exists a unique stationary distribution , such that Applying Itô’s formula to (60) and then integrating from to , we get From the first equation of system (44), it yields that Due to , one can getFrom the second equation of system (44) we have By Lemma 4, if

holds, then

In summary, one gets the following.

Theorem 10. Given initial value for system (44), then
(i) if is established, the prey population will be persistent in mean,
(ii) if conditions and are satisfied, the predator population will be persistent in mean.

Remark 11. (i) It can be seen from , in the case where the intrinsic growth rate and the predation intensity are relatively constant, only by reducing the intensity of white noise or pollutant concentration, so that can be established to ensure the lasting survival of the prey population.
(ii) Obviously, ; if the prey population is persistent, the predator population is persistent as long as the white noise interference intensity or the toxin concentration is small enough, such that is established. As can be seen from condition , the omnivorous nature of contributes to its permanence.

In system (4), if the concentration of the toxicant in the environment remains unchanged, that is, is a positive constant, then the system can be converted into the following system:

Lemma 12 (see [77]). System (7) is ergodic and positive recurrent if the following conditions are satisfied:
(i) For , .
(ii) For each ,for all , with some constant for all .
(iii) There exists a bounded open set with a smooth boundary satisfying that, for each , there is a twice continuously differentiable nonnegative function and that for some , , for any
Moreover, the Markov process has a unique ergodic stationary distribution . Hence, for any Borel measurable function , if , then

Recently, the ergodicity and stationary distribution have been explored by many authors. In the following, we give sufficient conditions for the existence of stationary distribution of system (68) and prove the following theorem by showing that system (68) satisfies the three conditions in Lemma 12.

Theorem 13. Assume that for , and, hold; then the stochastic process of system (68) is ergodic and has a unique stationary distribution in .

Proof. By the assumption for in Theorem 13, condition (i) in Lemma 12 is satisfied. Let ; then Define a bounded open subset as follows: where is a constant. Let ; here and . For , we have for all . Thus condition (ii) in Lemma 12 holds. Therefore, it remains for us to verify condition (iii) in Lemma 12.
Define a -function on , where , , is a sufficiently small positive constant satisfying , and . , are quantities to be determined below.
An application of the operator to yieldsDefine the vectors and with , . As the generator matrix is irreducible, for each and , there exists and , respectively, which is a solution of the Poisson system [78] Therefore we haveCombining (75), (77), and (78), we obtain Similarly, for , we calculate From conditions and , we can choose sufficiently small such that So can be estimated as follows: It is easy to see that Consequently, we derive that, for a sufficiently small , Using Lemma 12, we obtain the conclusion that is ergodic and positive recurrent; that is, system (68) is positive recurrent and has a unique stationary distribution.
This completes the proof of Theorem 13.

5. Conclusions and Numerical Simulations

In this article, we discussed the dynamics of stochastic prey-predator models with Beddington-DeAngelis functional response in polluted environment.

Firstly, for system (3), there are the following properties:

If and , the limit system of (3) has a prey extinction periodic solution .

If conditions , , and are established, the limit system of (3) has a positive periodic solution.

Secondly, system (4) possesses the following properties:

If , the prey population will be extinct.

If , the predator population will be extinct.

If and , then prey population and predator population will die out.

If , the prey population will be persistent in mean.

If and , the predator population will be persistent in mean.

To verify the correctness of the theoretical analysis, numerical simulations are employed in the following examples.

Assume that the Markov chain take values in with the generator By the linear equations , we can see which is the stationary distribution of . Furthermore, in the following examples, we suppose and , consistently.

5.1. The Existence of Periodic Solutions of System (3)

Example 14. Assume , , , , , , , , , , and .

Case 1. We choose the density of white noise as the following: , , and let .
Note that , The conditions of Theorem 6 hold, so there exists a boundary periodic solution of system (3) (see Figure 1).

Case 2. We change the density of the white noise to , , and . This gives , , , and ; choose according to and
From Theorem 7, we know that there exists a positive -periodic solution of system (3) (see Figure 2).

5.2. The Extinction and Persistence of System (4)

Example 15. Choose parameters , , , , , , , and , if , and , , , , , , , and , if .

Case 1. Let , , and ; we note , , and The conditions of Theorem 8 are satisfied, so the prey and predator are both extinct (see Figure 3).
Next we only change the density of the white noise to and and keep . Simple calculation shows that , , and . The conditions of Theorem 10 are satisfied, so the prey and predator are persistent (see Figure 4).
It is easy to see from Figures 3 and 4 that the increase of the intensity of white noise can result in the extinction of prey and predator.

Case 2. Let , , and , which gives , , and . and are extinct (see Figure 5).
Next we only change the amount of toxicant to and keep and , We note , , and . Thus and are persistent in mean (see Figure 6).
Figures 5 and 6 show that the increase of the amount of toxicant can also result in the extinction of the prey and predator.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work is supported by Shandong Provincial Natural Science Foundation of China (no. ZR2015AQ001), the National Natural Science Foundation of China (no. 11371230), and Research Funds for Joint Innovative Center for Safe and Effective Mining Technology and Equipment of Coal Resources by Shandong Province and SDUST (2014TDJH102).