- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Annual Issues ·
- Article Processing Charges ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Recently Accepted Articles ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents

Discrete Dynamics in Nature and Society

Volume 2013 (2013), Article ID 412409, 8 pages

http://dx.doi.org/10.1155/2013/412409

## Modeling the Dynamics of a Single-Species Model with Pollution Treatment in a Polluted Environment

^{1}College of Mathematics and Information Science, Anshan Normal University, Anshan, Liaoning 114007, China^{2}Department of Mathematics, Liaoning Normal University, Dalian, Liaoning 116029, China^{3}Department of Mathematics, Beihang University, Beijing 100083, China

Received 19 November 2012; Accepted 25 December 2012

Academic Editor: Zhen Jin

Copyright © 2013 Bing Liu et al. 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.

#### Abstract

Without any question, environmental pollution is the main cause for the species extinction in recent times. In this paper, based on impulsive differential equation, the dynamics of a single-species model with impulsive pollution treatment at fixed time in a polluted environment is considered, in which we assume that the species is directly affected by the pollutants. Sufficient conditions for permanence and extinction of the species are given. The results show that the species is permanent when the impulsive period is less than some critical value, otherwise the species will be extinct. Although shortening the impulsive period can protect the species from extinction, it is expensive. To see how pollution treatment applications could be economical, we also establish a hybrid impulsive model involving periodic pollution treatment at fixed time with state-dependent pollution treatment applied when the pollution concentration reaches the given Environment Threshold (ET). It indicates that the hybrid method is the most effective method to protect the species from extinction. Numerical simulations confirm our theoretical results.

#### 1. Introduction

With the globalization of economy and trade, environmental pollution increasingly becomes one of the most serious problems faced by all the countries in the world. For instance, the greenhouse effect, acid rain, and ozone depletion are the environmental effects of air pollution. Another problem that the whole world is concerned about is how to protect endangered species. Of the ecosystem types, the largest species extinction happened to tropical rainforests, including many species not investigated and named by people. Tropical forests are home to more than one-half of the Earth’s known species. The scientists have predicted that in the following 30 years the tropical rainforests abundant in species will likely be destroyed by the contemporaries and that 5–10 species in the tropical rainforests will likely tend to extinction. Environmental pollution is a major contributor to species extinction or near extinction and the further reduction of biodiversity. It affects the distribution of species in nature, reduces habitat animals live in, and causes mass migration of animal species and also has a negative effect on growth and reproduction of organisms, and when pollution is serious, organisms will change markedly in morphological specificity, their number, and the like. To preserve the ecological environment, we must curb environmental pollution.

After the pioneering theoretical work [1, 2], some papers [3–9] have been researched on the effect of the continuous input pollution on the survival of a biological population. Considering the pollutants are emitted in regular pulses, Liu et al. have been focusing on the effect of impulsive pollution input on a single species [10–14] and two-species Lotka-Volterra competition system [15] and have obtained the conditions of the persistence and extinction for the population. Reference [16] investigated a logistic offshore fishery system with impulsive pollutant input in inshore areas and impulsive diffusion at different fixed time. They obtained the sufficient conditions of the existence of the positive periodic solution and the global asymptotic stability of both the trivial periodic solution and the positive periodic solution and gave the corresponding optimal harvesting strategy.

Based on the biological background mentioned above and with environmental pollution worsening, to protect endangered species, in this paper a single-species model with impulsive pollution treatment at fixed time in a polluted environment is established and studied. So far, no research has been done on how to use the mathematical model to study the effect of pollution control on population size. Besides, according to practical problems, we will consider constructing a hybrid impulsive single-species model combining periodic pollution treatment at fixed time with state-dependent pollution treatment applied when the pollution concentration reaches the given Environment Threshold and gives specific biological explanations based on numerical simulations, and our conclusion will provide the theoretical decision-making basis for practical problems.

The organization of this paper is as follows. In the next section, a continuous model is proposed and its corresponding dynamic properties is obtained, then we establish a single-species model with impulsive pollution treatment at fixed time in a polluted environment. In Section 3, we obtain the conditions for permanence and extinction of such a system. The conditions of existence and stability of positive periodic solution of model are given in Section 4. In Section 5, we propose a hybrid impulsive single-species model with Environment Threshold and give biological explanations by numerical simulations. In the last section, we discuss our results and suggest future work.

#### 2. Model Formulation

Firstly, we assume that the emission of pollutants to the environment is continuous and directly affects the survival of the species in such an environment. Let us consider the following single species model with the pollution effect: where is the density of the species at time ; is the concentration of pollution in the environment at time ; is the carrying capacity of environment; is the intrinsic growth rate of the species in the absence of pollutants in the environment; is the dose-response parameter of species to the pollution concentration in the environment; represents the exogenous rate of the pollutants input into the environment and is bounded, for the convenience, here we assume ; represents the totality of losses from the environment including biological transformation, chemical hydrolysis, volatilization, microbial degradation, and photosynthetic degradation.

In this paper, we assume that the capacity of the environment is so large that the change of the concentration of pollution in the environment that comes from uptake and egestion by the species can be ignored.

Clearly in system (1), and . System (1) always has a species extinction equilibrium ; if , there exists a unique positive equilibrium .

Let be any equilibrium of system (1), then the characteristic equation of its equilibrium is For the equilibrium , (2) is reduced to If , both roots of (3) have negative real part, so is locally asymptotically stable.

If , one of the roots of (3) is positive, then is unstable saddle, and positive equilibrium point exists.

For this equilibrium , (2) is reduced to If , both roots of (4) have negative real part, so is locally asymptotically stable.

Theorem 1. *If , the species extinction equilibrium is globally asymptotic stable; whereas if , the positive equilibrium is globally asymptotic stable.*

* Proof. *If , is locally stable. In this case, there is no any interior equilibrium in system (1), so there cannot be any periodic solutions in since a periodic solution must contain at least one equilibrium. By the Poincaré-Bendixson theory, is globally asymptotic stable.

If , then exists, and is locally asymptotically stable. Let , be the right sides of (5), respectively, then we have
Set Dulac function , then we have
which implies that there is no limit cycle in , so the equilibrium is global asymptotically stable. This completes the proof.

*Remark 2. *Since can be rewritten as , under this condition the growth rate of the species becomes . So the species will tend to be extinct; that is, the concentration of pollution in the environment is violence against the survival of the species. If , it means that the growth rate of the species , and the species will be permanent and tend to the unique globally asymptotically stable positive equilibrium.

From the above theorem, we know if , is globally asymptotically stable, that is, the species will tend to be extinct. In order to protect the population from extinction, we must take steps to control environmental pollution. Now based on the impulsive differential equation [10–22], we formulate the following single species model with impulsive pollution treatment at fixed time: where , ; , , , are positive constants, and ; is the period of impulsive effect and is the decrease rate of pollution concentration with the effects of pollution treatment. The other meanings of parameters are the same as those of model (1).

*Remark 3. * In this paper, we consider how to protect the species from extinction when the pollution concentration is large in the environment, so we assume that .

*Remark 4. *In the model (7), is the concentration of pollution in the environment, so . To assure that remains less than one, it is necessary that
In the following, we always assume that (8) holds true.

#### 3. Extinction and Permanence

In this section, we investigate the extinction and permanence of system (7).

First, we consider the following subsystem of system (7): where .

The following lemma is a corollary of Theorem 2.3 in [23].

Lemma 5 (see [23]). *System (9) has a unique globally asymptotically stable positive -periodic solution , where
*

Therefore, system (7) has a species extinction periodic solution .

*Remark 6. *In model (9), tends to a unique positive -periodic solution , which shows that periodic impulsive pollution treatment causes a periodic behavior of in system (9).

Lemma 7. *Let , then has a unique positive root , and one has*(1)*if , ;*(2)*if , .*

* Proof. *Obviously .

Calculate the derivative of , we have
Obviously, ; Let , , since . So we have which has a unique positive root . If , we have , and then holds true for ; if , we have . Since , , so we have , which has a unique positive root , and . Thus if , we have ; if , we have . The proof is completed.

Now, we give sufficient conditions which drive to be extinct.

Theorem 8. *If , then the species extinction periodic solution is globally asymptotically stable.*

*Proof. *From Lemma 7, if , . We choose sufficiently small such that
From Lemma 5, we have as , that is, for being large enough,
holds. For the sake of simplicity, we assume that (13) holds true for all . From system (7), we have
thus
so we get and as . Therefore as since for . This completes the proof.

*Definition 9. *It is said that the species of system (7) is permanent if there exists a positive number and a finite time such that for all .

On the permanence of system (7), we have the following result.

Theorem 10. *The species of system (7) is permanent if .*

*Proof. *From Lemma 7, if, . We choose sufficiently small such that
From Lemma 5, we know for the above , there exists a such that
holds true for .

We claim that the inequality cannot hold for all . Otherwise, we have
Let , and . Integrating inequality (18) on , we have
so , which is a contradiction to the boundness of . Hence, there exists a such that . In the following, we will prove that for all , where
Otherwise, there exists a such that . Then there exists a such that and for . Let , , we have
this is a contradiction. Let , we can obtain for all . This completes the proof.

*Remark 11. *From Theorems 8 and 10, we obtain that the impulsive period is the threshold value which determines permanence or extinction of the species. It implies that if the impulsive period is less than the threshold value , then the pollution in the environment is treated frequently and the concentration of the pollution in the environment decreases rapidly, and eventually the species will be permanent. In the next section, we will see in this case the species will tend to a positive globally asymptotically stable -periodic solution; if the impulsive period is larger than the threshold value , then the pollution in the environment is not treated in time and will drive the species to be extinct. Therefore, in order to protect the endangered species, we should ensure that holds true, which can be carried by several strategies, such as increasing the treatment rate , increasing the frequency of pollution treatment, or controlling the emission of pollutants.

#### 4. Existence and Stability of Positive Periodic Solution of Model (7)

Now we give the condition of existence and stability of positive periodic solution of model (7). First of all, we consider the following periodic logistic equation where , are -periodic continuous on . A solution of (22) is said to be positive, if for all . It is easy to prove that the solution of (22) is positive if and only if the initial .

Lemma 12 (see [24]). *If and , then (22) has a unique positive -periodic solution which is globally asymptotically stable, that is; as for any positive solution of (22).*

We give a continuous function , where is a -periodic continuous function; if , we call is asymptotical . Consider the following unique solution of the Cauchy function: where is asymptotical , is a -periodic continuous function.

Lemma 13 (see [24]). *If and , for all , is a unique positive -periodic solution of (22), is the solution with initial value of (23), then for any initial values , one has as .*

In the following, we give sufficient conditions of existence and stability of positive periodic solution by using Lemmas 12 and 13.

Theorem 14. *For system (7), if , then there exists a unique positive -periodic solution which is globally asymptotically stable, that is, for any solution of system (7), we have
*

*Proof. *From Lemma 5, we know that as , so we only prove if the condition of Theorem 14 holds, then in system (7) there exists a positive -periodic, and , as . In the following, we use the equivalent system of (7) to prove this theorem.

Consider the following equation:
For
from Lemma 12, we know that in system (25) there exists a unique globally asymptotically positive -periodic solution, denoted by . Since is a -periodic solution of system (7), so is a -periodic solution of the system (7). And , as , we obtain that is asymptotical , and from Lemma 13, for any solution of system (7), we have as . The proof is completed.

*Remark 15. *From Theorem 14, we can see that the condition of the existence and stability of positive periodic solution of model (7) is consistent with the permanent condition. Therefore, if the impulsive period is smaller than the critical value , through the effective pollution treatment, the concentration of pollution will not be sufficiently high to kill the species, and the species will be permanent and tends to a unique positive -periodic solution of system (7) (see Figure 1(a)). Otherwise, the species will tend to extinction (see Figure 1(b)).

Finally, we evaluate the average density of the species in one period when system (7) has a positive periodic solution.

Since is a positive periodic solution of system (7), we have which gives so which implies that the average number of species strictly depends on impulsive period and pollution emission rate. The population will be extinct if .

#### 5. A Hybrid Impulsive Model with Environment Threshold

Because of economic and technical reasons, although the factory takes measures to treat pollution at regular time before the pollutants are emitted to the environment, sometimes the concentration of the pollution in the environment is still large and threatens the survival of the species. In practice, the factory is required to make pollution treatment immediately when the concentration of pollution reaches the given Environment Threshold (ET). According to the above biological background, we establish the following hybrid impulsive model involving periodic pollution treatment at fixed time with state-dependent pollution treatment applied when the pollution concentration reaches the given Environment Threshold: where is periodic impulsive point series at which pollution treatment is taken, is the environmental threshold (ET) which threatens the survival of the species, and we assume that the initial concentration of pollution is less than .

From the biological point of view, experimental methods in combination with the model approaches presented in this paper are the most effective methods. We have shown that the species of system (7) with parameters in Figure 1(b) will tend to be extinct. If we only treat pollution at the fixed time to protect the species from extinction, we must shorten the impulsive period as shown in Figure 1(a). Sometimes it is difficult to put it into practice. Now we use the hybrid method to protect the species from extinction. Set ; ; ; and other parameters are the same as those in Figure 1(b). When the concentration of pollution reaches ET, we treat pollution immediately. Then from Figure 2, we can see that the solution of system (29) will tend to a globally stable periodic solution, which means the species will be permanent.

#### 6. Discussion

In this paper, we investigated a single-species model with impulsive pollution treatment at fixed time in a polluted environment. From Theorems 8 and 10, we can see that was the threshold condition under which the species became extinct or permanent. As long as was small enough and less than this threshold, the survival was predicted. Therefore, this threshold provided us with a way to choose the impulsive period and pollution treatment rate in order to protect the species from extinction. To see how pollution treatment applications could be economical, we also established a hybrid impulsive model involving periodic treatment at fixed time with state-dependent pollution treatment applied when the pollution concentration reached the given Environment Threshold. The results indicated that this hybrid method was the most effective method to protect the species from extinction.

Studying the effect of pollution treatment on the protection of the endangered species is at the very beginning. There is still much work to be done. For example, we have not studied the effect of pollution treatment on the single species with stage structure and time delay and two-species model in a polluted environment. We just leave these for the future work.

#### Acknowledgments

This work is supported by the National Natural Science Foundation of China (10971001, 111011021) and Science and Excellent Talents Support Project of Universities and Colleges in Liaoning China.

#### 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 Scopus - T. G. Hallam, C. E. Clark, and G. S. Jordan, “Effects of toxicants on populations: a qualitative approach. II. First order kinetics,”
*Journal of Mathematical Biology*, vol. 18, no. 1, pp. 25–37, 1983. View at Scopus - 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 · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - P. D. N. Srinivasu, “Control of environmental pollution to conserve a population,”
*Nonlinear Analysis. Real World Applications. An International Multidisciplinary Journal*, vol. 3, no. 3, pp. 397–411, 2002. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at 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 · View at Google Scholar · View at MathSciNet - D. Mukherjee, “Persistence and global stability of a population in a polluted environment with delay,”
*Journal of Biological Systems*, vol. 10, no. 3, pp. 225–232, 2002. View at Publisher · View at Google Scholar · View at Scopus - F. Wang and Z. Ma, “Persistence and periodic orbits for an SIS model in a polluted environment,”
*Computers & Mathematics with Applications*, vol. 47, no. 4-5, pp. 779–792, 2004. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Y. Xiao and L. Chen, “How do the spatial structure and time delay affect the persistence of a polluted species,”
*Applicable Analysis*, vol. 82, no. 3, pp. 253–267, 2003. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - A. K. Pal and G. P. Samanta, “A single species population in a polluted environment,”
*International Journal of Biomathematics*, vol. 3, no. 2, pp. 187–204, 2010. View at Publisher · View at Google Scholar · View at 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 · View at Google Scholar · View at Scopus - 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 · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - X.-Z. Meng, Q.-L. Zhao, and L.-S. Chen, “Global qualitative analysis of new monod type chemostat model with delayed growth response and pulsed input in polluted environment,”
*Applied Mathematics and Mechanics*, vol. 29, no. 1, pp. 75–87, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - B. Liu, L. Zhang, and Q. Zhang, “The effects of a single stage-structured population model with impulsive toxin input and time delays in a polluted environment,”
*Applicable Analysis*, vol. 88, no. 8, pp. 1143–1155, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Z. Zhao, L. Chen, and X. Song, “Extinction and permanence of chemostat model with pulsed input in a polluted environment,”
*Communications in Nonlinear Science and Numerical Simulation*, vol. 14, no. 4, pp. 1737–1745, 2009. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at 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 · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - B. Kang, B. Liu, and L. Xu, “Dynamics of an inshore-offshore fishery model with impulsive pollutant input in inshore area,”
*Nonlinear Dynamics*, vol. 67, no. 4, pp. 2353–2362, 2012. View at Publisher · View at Google Scholar · View at Scopus - B. Liu, Y. Duan, and S. Luan, “Dynamics of an SI epidemic model with external effects in a polluted environment,”
*Nonlinear Analysis. Real World Applications. An International Multidisciplinary Journal*, vol. 13, no. 1, pp. 27–38, 2012. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - S. Sun and L. Chen, “Permanence and complexity of the eco-epidemiological model with impulsive perturbation,”
*International Journal of Biomathematics*, vol. 1, no. 2, pp. 121–132, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Y. Xue, A. Kang, and Z. Jin, “The existence of positive periodic solutions of an eco-epidemic model with impulsive birth,”
*International Journal of Biomathematics*, vol. 1, no. 3, pp. 327–337, 2008. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - J. Jia and H. Cao, “Dynamic complexities of Holling type II functional response predator-prey system with digest delay and impulsive harvesting on the prey,”
*International Journal of Biomathematics*, vol. 2, no. 2, pp. 229–242, 2009. View at Publisher · View at Google Scholar · View at MathSciNet - S. Tang, G. Tang, and R. A. Cheke, “Optimum timing for integrated pest management: modelling rates of pesticide application and natural enemy releases,”
*Journal of Theoretical Biology*, vol. 264, no. 2, pp. 623–638, 2010. View at Publisher · View at Google Scholar · View at Scopus - W. Gao and S. Tang, “The effects of impulsive releasing methods of natural enemies on pest control and dynamical complexity,”
*Nonlinear Analysis. Hybrid Systems. An International Multidisciplinary Journal*, vol. 5, no. 3, pp. 540–553, 2011. View at Publisher · View at Google Scholar · View at Zentralblatt MATH · View at MathSciNet - Z. Hu and M. Han, “Periodic solutions and bifurcations of first-order periodic impulsive differential equations,”
*International Journal of Bifurcation and Chaos*, vol. 19, no. 8, pp. 2515–2530, 2009. View at Publisher · View at Google Scholar · View at Scopus - J. Lopez-Gomez, R. Ortega, and A. Tineo, “The periodic predator-prey Lotka-Volterra model,”
*Advances in Differential Equations*, vol. 1, no. 3, pp. 403–423, 1996. View at Zentralblatt MATH · View at MathSciNet