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

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.


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][4][5][6][7][8][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][11][12][13][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 statedependent 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 singlespecies 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.

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;  1 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 lim  → ∞ () = ; −ℎ() 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.
If  >  1 /ℎ, one of the roots of (3) is positive, then  1 is unstable saddle, and positive equilibrium point  2 exists.
Proof.If  <  1 /ℎ,  1 is locally stable.In this case, there is no any interior equilibrium in system (1), so there cannot be any periodic solutions in  2 + since a periodic solution must contain at least one equilibrium.By the Poincaré-Bendixson theory,  1 is globally asymptotic stable.
If  >  1 /ℎ, then  2 exists, and  2 is locally asymptotically stable.Let (, ), (, ) be the right sides of (5), respectively, then we have Set Dulac function (, ) =  −1  −1 , then we have which implies that there is no limit cycle in  2 + , so the equilibrium  2 is global asymptotically stable.This completes the proof.
Remark 2. Since  <  1 /ℎ can be rewritten as  <  1  * , under this condition the growth rate of the species becomes  −  1  * < 0. 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  >  1 /ℎ, it means that the growth rate of the species − 1  * > 0, and the species will be permanent and tend to the unique globally asymptotically stable positive equilibrium.
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  <  1 /ℎ.
Remark 4. In the model (7), () is the concentration of pollution in the environment, so 0 ≤ () ≤ 1.To assure that  0 () remains less than one, it is necessary that In the following, we always assume that (8) holds true.
The following lemma is a corollary of Theorem 2.3 in [23].
Remark 6.In model (9), () tends to a unique positive periodic solution c(), which shows that periodic impulsive pollution treatment causes a periodic behavior of () in system (9).
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.
Remark 11.From Theorems 8 and 10, we obtain that the impulsive period  0 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  0 , 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  0 , 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  <  0 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.

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 () > 0 for all  ≥ 0. It is easy to prove that the solution () of ( 22) is positive if and only if the initial (0) > 0.
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 0 <  <  0 , then there exists a unique positive -periodic solution ( x(), c()) which is globally asymptotically stable, that is, for any solution ((), ()) of system (7), we have Proof.From Lemma 5, we know that () → c() as  → ∞, so we only prove if the condition of Theorem 14 holds, then in system (7) there exists a positive -periodic, and () → x(), 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 x().Since c() is a -periodic solution of system (7), so x() is a -periodic solution of the system (7).And () → c(), as  → ∞, we obtain that () is asymptotical c(), and from Lemma 13, for any solution () of system (7), we have () → x() 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  0 , 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 x() is a positive periodic solution of system (7), we have (0) = () 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  >  0 .

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: Figure 1: Time series of system (7).(a) Extinction solution of system (7) with parameters  = 0.5;  = 2;  1 = 0.7;  = 0.8; ℎ = 0.9;  = 4;  = 0.55; (0) = 0.6; (0) = 0.5.(b) -periodic solution of system (7) with parameters  = 0.5;  = 2;  1 = 0.7;  = 0.8; ℎ = 0.9;  = 2;  = 0.55; (0) = 0.6; (0) = 0.5.
where  ( = 1, 2, . ..) is periodic impulsive point series at which pollution treatment is taken,  0 is the environmental threshold (ET) which threatens the survival of the species, and we assume that the initial concentration of pollution is less than  0 .
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 ET = 0.8;  1 = 0.4;  2 = 0.55; 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.

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  0 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.