Analysis of Mathematics and Sustainability in an Impulsive Eutrophication Controlling System

and Applied Analysis 3 (Moina, Cyclops, and Brachionus) (10 piece/L), and the filterfeeding fish (silver carp and bighead carp) (1 fish/Cubic meter) at time t, respectively, Δx(t) = x(t) − x(t), Δy(t) = y(t + ) − y(t), and Δz(t) = z(t) − z(t). 0 < l < 1 is used to describe the intervals of time between the pulsed use of controls. r i (i = 1, 2) are the intrinsic growth rate, a i (i = 1, 2, 3) and b i (i = 1, 2, 3) measure the efficience of the prey in evading a predator attack. e i (i = 1, 2, 3) denote the efficiency with which resources are converted to new consumers, k i (i = 1, 2) are carrying capacity in the absence of the filterfeeding fish, and m is the mortality rates for filter-feeding fish. The decision-maker of biological treatment technology must implement an impulsive control strategy that can better maintain biological treatment process, so we not only release a certain amount of the filter-feeding fishr z (p > 0 is the release amount of the filter-feeding fishr z) at t = nT to graze the algae and the dominant zooplankton, but also consider the algae and the dominant zooplankton losses due to other reasons as well as the filter-feeding fish harvest at t = (n + l − 1)T (δ i (i = 1, 2, 3) (0 < δ i < 1) represent the fraction of three species), where T is the period of the impulsive effect, n ∈ N, andN is the set of all nonnegative integers. Themain purpose of this paper is to enhance a theoretical and constructive framework by making the attempt to deal with eutrophication controlling problem for some drinking reservoirs using an impulsive eutrophication controlling system. Meanwhile, the theory discussions on the global asymptotical stability and permanent would help to design some efficient parameter constraint representations, which can in turn provide a theoretical basis for the numerical simulation. Further, in the context of population growth dynamics, how to implement impulsive control strategy to prevent and control the algal bloom is studied by simulating the dynamics of the impulsive controlling eutrophication system. 2. Mathematical Analysis Nowwewill pursue the investigation of the constraint expression of the release amount p and the average harvest δ 3 of the filter-feeding fish. Let R + = [0,∞), and R + = {X ∈ R 3 | X ≥ 0}. Denote f = (f 1 , f 2 , f 3 ) as the map defined by the right hand of the first, second, and third equations of system (1). Let V: R + × R 3 + → R + ; then V is said to belong to class V 0 if (1) V is continuous in (nT, (n + l − 1)T) ∪ ((n + l − 1)T, (n + 1)T] × R 3 + , and for each X ∈ R + , n ∈ N, lim (t,y)→ (nT + ,X) V(t, y) = V(nT + , X) exists; (2) V is locally Lipschitzian in X. Definition 1. Let V ∈ V 0 ; then for (t, x) ∈ (nT, (n + l − 1)T) ∪ ((n + l − 1)T, (n + 1)T] × R 3 + , the upper right derivative of V(t, X) with respect to the impulsive differential system (1) is defined as


Introduction
The eutrophication of lakes and reservoirs has been widely and intensively studied for several decades, which is a degradation process originating from introduction of nutrients from agricultural runoff and untreated industrial and urban discharges [1]. The eutrophication has been identified as a major and serious water quality management issue [2], not only in the developing countries [3][4][5], but also in the developed countries [6][7][8]. Impairment of water quality due to eutrophication can lead to a series of problems and result in losses of ecological integrity, sustainability, and safe use of aquatic ecosystems [9]. It is characterized by frequently recurring algal blooms, zooplankton aggregation, and reduction in species diversity in water bodies at all trophic levels [2].
In 1932, Bertalanffy firstly proposed a method of using mathematical model to study the biological system [10].
He pointed out that this approach is a combination of coordination, order, and purpose, which in turn formed three basic ideas of studying the biological system, namely, trophiclevel analysis, system perspective, and dynamic view. The main aim of modeling population dynamics is to improve the understanding of the interactions between populations and their dependence on internal and external conditions [11]. Hence, the mathematical models of biological population dynamics have to account not only for the growth and interactions between populations but also for the forecast of the disasters which can be caused by the plankton like blue-green algae, Moina and Brachionus. Early attempts began with the logistic growth, exponential growth, and Cui growth (this growth can explain the relationship between population increment and limiting resources) for biological population. With the increasing effects of the eutrophication, zooplankton aggregation, and algal blooms, the history of mathematical modeling plankton or zooplankton dynamics and biological eutrophication removal process is already quite long and has been initiated by biological science. These approaches have been refined to more realistic descriptions of the development of biological natural populations, especially algae population. In recent decades, the main researches focused on how to control the eutrophication, how to predict the outbreak of algae, how to simulate the spread tendency of algae, and how to develop eutrophication control strategy based on actual data. Hence, more and more researches have given an attempt to combine the forces of ecologist and mathematical engineers in an evaluation of possible control strategies for the management of eutrophication of water bodies, which is becoming an urgent watermanagement problem of drinking water reservoirs [12]. Malmaeus and Håkanson have depicted the development of a lake eutrophication model; they indicated that all changes that provided significant differences and improved predictions have been adopted [13]. Chen et al. said that their mathematical description provides a different insight into how to use nitrogen and phosphorus more effectively in the eutrophication control and elimination of blue-green blooms [14].
Zeya reservoir of Wenzhou is located in subtropical regions. Because of eutrophication, the nuisance algal blooms and the serious zooplankton aggregation frequently come forth, which caused clogging and blocking of the filtration system and resulted in millions of people without drinking water. Eutrophication removal is achieved by two major processes, physicochemical and biological treatment techniques; specifically biological eutrophication removal from drinking water in the lake and reservoir is usually considered to accomplish optimal and economic eutrophication treatment [15]. In order to understand and improve the performance of biological eutrophication removal treatment, it is necessary to briefly look into the mechanisms of eutrophication removal. The removal process of eutrophication in the Zeya reservoir is complex and dynamic with many variables. Since successful management of eutrophication in Zeya reservoir is often dependent on proper identification of the most likely growth limiting amount for algae, special emphasis is given to the limiting amount concept. In order to apply the biological principle to control eutrophication in the Zeya reservoir, we give good biological control agents for algae. First, we bring up two species of filter-feeding fish, silver carp (Hypophthalmichthys molitrix), and bighead carp (Aristichthys nobilis) [16], which are thought to be suited to control algae and zooplankton biomass directly in freshwater reservoir, especially cyanobacteria. Second, because the decision-maker must give cost-benefit and maintenance considerations, fish should be carried out and controlled at an appropriate time. At the same time, in order to prevent the rapid growth of plankton, especially algae, and maintain a good circulating development of biological process, we must release a certain amount of filter-feeding fish at a fixed time. These processes are subject to short-term perturbations whose durations are too small compared to the duration of the process. Consequently, it is natural to assume that these perturbations act instantaneously in the form of impulses. Thus impulsive differential equations and differential equations involving impulsive effects appear as a natural mathematical description (model) of observed evolution phenomena of several real world problems [17][18][19][20][21]. These works have introduced the theory of impulsive differential equations and its applications and provided a theoretical basis of the application of impulsive differential equations. In recent years, Fan et al. [22] have proposed the dynamics of a class of mutualismcompetition-predator interaction models with Beddington-DeAngelis functional responses and impulsive perturbations; this work is important for solving the global stability and the globally exponential stability of system. Yu and Lu [23] mainly focused on the permanence and global attractivity of a discrete almost periodic ratio-dependent Leslie system with time delays and feedback controls; these works vigorously promote the application of feedback control. Shen et al. [24] mainly have studied the permanence of a single species system with distributed time delay and feedback controls; it is a very important fact that the feedback control is harmless to the permanence of species; these works are significant for the application of feedback control in biological control. Furthermore, the main aim of this investigation is to answer that what is the fundamental mathematical controlling mechanism, which can support the biological principle to control eutrophication in the Zeya reservoir and to provide a different insight into how to make the effects of environmental factors in connection with the ability of controlling eutrophication.
Speaking on eutrophication of Zeya reservoir, as the algae population grows to a certain biomass level, the dominant zooplankton species (Moina, Cyclops, and Brachionus are mainly grazing algae) are bound to sire and grow, so that these results can lead to damage or even crash of aquatic ecosystems. Furthermore, based on the above discussion, modeling is playing an increasing role in helping to identify causes and effective control measures through what if studies, and strategies are being identified for the management and control of eutrophication [25,26]; the paper considers an eutrophication mathematical model, which can be described by the following differential equations: where ( ), ( ), and ( ) are the biomass of the algae population (cyanobacteria) (10 6 cell/L), the dominant zooplankton Abstract and Applied Analysis 3 (Moina, Cyclops, and Brachionus) (10 2 piece/L), and the filterfeeding fish (silver carp and bighead carp) (1 fish/Cubic meter) at time , respectively, Δ ( ) = ( + ) − ( ), Δ ( ) = ( + ) − ( ), and Δ ( ) = ( + ) − ( ). 0 < < 1 is used to describe the intervals of time between the pulsed use of controls. ( = 1, 2) are the intrinsic growth rate, ( = 1, 2, 3) and ( = 1, 2, 3) measure the efficience of the prey in evading a predator attack. ( = 1, 2, 3) denote the efficiency with which resources are converted to new consumers, ( = 1, 2) are carrying capacity in the absence of the filterfeeding fish, and is the mortality rates for filter-feeding fish. The decision-maker of biological treatment technology must implement an impulsive control strategy that can better maintain biological treatment process, so we not only release a certain amount of the filter-feeding fishr ( > 0 is the release amount of the filter-feeding fishr ) at = to graze the algae and the dominant zooplankton, but also consider the algae and the dominant zooplankton losses due to other reasons as well as the filter-feeding fish harvest at = ( + − 1) ( ( = 1, 2, 3) (0 < < 1) represent the fraction of three species), where is the period of the impulsive effect, ∈ , and is the set of all nonnegative integers.
The main purpose of this paper is to enhance a theoretical and constructive framework by making the attempt to deal with eutrophication controlling problem for some drinking reservoirs using an impulsive eutrophication controlling system. Meanwhile, the theory discussions on the global asymptotical stability and permanent would help to design some efficient parameter constraint representations, which can in turn provide a theoretical basis for the numerical simulation. Further, in the context of population growth dynamics, how to implement impulsive control strategy to prevent and control the algal bloom is studied by simulating the dynamics of the impulsive controlling eutrophication system.
The following lemma is obvious.
We will use an important comparison theorem on impulsive differential equation.
Lemma 4 (Bainov and Simeonov, 1993 [17]). Suppose ∈ 0 . Assume that where : Finally, we give some basic properties about the following subsystems of system (1): 4 Abstract and Applied Analysis ∈ , and * (0 is a positive periodic solution of system (5). Since is the solution of system (5) with initial value 0 = (0 + ) ≥ 0, where ∈ (( − 1) , ( + − 1) ) ∪ (( + − 1) , ], ∈ , we get the following. Therefore, we obtain the complete expression for a periodic solution (0, 0, * ( )) of system (1). Now, we study the stability of the periodic solution, which explains the threshold expression of the release amount and the average harvest 3 of the filter-feeding fish under the condition of some species extinction. Theorem 7. Let ( ( ), ( ), ( )) be any solution of system (1); then (0, 0, * ( )) is said to be globally asymptotically stable if Proof. The local stability of periodic solution (0, 0, * ( )) may be determined by considering the behavior of small amplitude perturbation of the solution. Define We will put (10) into (1). The linearization of the system becomes Therefore, we have where Φ( ) satisfies Abstract and Applied Analysis 5 and Φ(0) = , the identity matrix, and ) .
Now we study the stability of the periodic solution, which explains the threshold expression of the release amount and the average harvest 3 of the filter-feeding fish under the condition of all species persistence. Thus we only need to find 2 and 8 such that ( ) > 2 and ( ) > 8 for large enough. Now we will find a 2 .
Using the theoretical analysis, we obtain the threshold expression of the release amount and the average harvest 3 of the filter-feeding fish under the condition of some species extinction and all species persistence, which in turn provides a theoretical basis for the numerical simulation.

Numerical Analysis
The effective methods used for eutrophication removal, algal bloom control process, and associated public perception are a crucial consideration in developing effective systems of eutrophication removal management for Zeya reservoir; hence the sustainability analysis of the impulsive eutrophication controlling system must be investigated in detail. In order to study how to choose the value of the parameters and 3 so that they can satisfy the actual monitoring data and enable Theorems 7 and 8 to be feasible, it is necessary to present a more in-depth analysis of the relationship among and 3 . Set 1 ( , 3 (1 − 2 )). On the basis of the actual monitoring data of the Zeya reservoir from year 2005 to year 2011 and laboratory controlling experiment data, the values of the parameters can be taken as follows: 1 = 0.75, 2 = 0.8, 1 = 20, 2 = 10, . It is obvious to find from Figure 1(a) that there exist some values of and 3 so that 1 ( , 3 ) < 0 and 2 ( , 3 ) < 0 or 1 ( , 3 ) > 0 and 2 ( , 3 ) > 0, which in turn confirms the existence of Theorems 7 and 8. From the actual stocking amount of the little filter-feeding fish in Zeya reservoir to control the algal bloom, we can choose biologically feasible parametric numerical value of as = 35. Further, it can be seen from Figure 1(b) that if the value of 3 is less than 0.4, there exist 1 ( , 3 ) > 0 and 2 ( , 3 ) > 0, but if the value of 3 is greater than or equal to 0.4, then the system (1) remains locally asymptotically stable. It is to be noted from Figure 2 that the time series solution of the filter-feeding fish population oscillates with stable limit cycle with increasing time. In contrast, both the algae population and the dominant zooplankton population rapidly decrease to their extinction with the increasing time. Meanwhile, these results further demonstrate that Theorems 7 and 8 are correct and feasible.
Considering the great uncertainty that still surrounds biological treatment mechanism, it is most prudent to interpret the simulation output on the basis of information and data obtained from observation and experiment. Simulation design offers tremendous flexibility to study the behavior dynamic of biological treatment process, so it is necessary to investigate that the release amount is how to affect the dynamics of biological treatment process with the parameter 3 = 0.4 and 0 < ≤ 35. It is worth pointing out from Figure 3 that the system (1) with 0 < < 35 exhibits relatively complex dynamical behaviors: chaos, period-halving bifurcation, and periodic attractor, but the system (1) with = 35 remains locally asymptotically stable. Furthermore, it should be stressed that the dominant zooplankton population rapidly decreases to their extinction when the value of is greater than 4.4268; this phenomenon shows that the dominant zooplankton population will advance toward extinction with the increase of . But it is regretful for this phenomenon because the dominant zooplankton population is one of the main food for the filter-feeding fish and can graze the algae population, and the idea of using biological eutrophication removal treatment to control the algal bloom not only can control the algal Abstract and Applied Analysis bloom, but also consider the sustainable development of eutrophication removal process with the maximizing economic benefits. Thus, controlling the desirable biomass of the dominant zooplankton population is a primary aim of operation and maintenance, and how to well carry out the eutrophication removal process is an extremely vital problem of Zeya reservoir. In real ecosystem, the effect of some critical factors can be depicted by the properties and the dynamic behaviors of species [27,28]. To more carefully investigate how to legitimately implement eutrophication removal treatment, Figure 4 depicts the long-run dynamics of system (1) for 0 ≤ 3 ≤ 1 with = 5. It is easy to find from Figure 4 that the system exhibits relatively complex dynamical behaviors: chaos, period-halving bifurcation, and periodic attractor, but it is worth pointing out that, when the value of 3 is beyond 0.6, the algae population will rapidly go to extinction and the dominant zooplankton population will go toward a reasonable state, and the biomass of the dominant zooplankton population is so small that it does not cause the aggregation of the dominant zooplankton. Meanwhile, the biomass of the filter-feeding fish is feasible to maintain the sustainable development of eutrophication removal process with the maximizing economic benefits because of = 5. This result suggests that we can choose the value of and 3 as 5 and 0.6 for legitimately implementing eutrophication removal treatment. For further analysis of the feasibility of this result, Figure 5 gives the dynamics of the system (1) with the increase of as well as = 5 and 3 = 0.6. It is gratifying to see from Figure 5 that the algae population quickly goes to the extinct status with the increase of , and the dominant zooplankton population is close to a reasonable number with the increase of . It is worth encouraging that the filter-feeding fish can arrive to a viable amount so that the eutrophication removal process can be developed with the maximizing economic benefits. Therefore, it should be stressed that these results can perfectly confirm that the impulsive eutrophication controlling system can be well applied in the process of controlling the algal bloom in the theoretical work of eutrophication removal and mathematical model.
Based on the foregoing analysis, it is successful for impulsive eutrophication controlling system to implement the eutrophication removal process. The impulsive eutrophication controlling system not only can promote all wanted species persistence, but also is convenient for the manager to maintain the normal development of biological treatment process. Moreover, simulation analysis provides an approximation of the real biological controlling system behaviors; hence it is very easy to obtain the best value of and 3 , Abstract and Applied Analysis which is in accordance with observations and in turn shows that the selection of the best value about some critical parameters in the view of population dynamics is feasible.

Conclusions and Remarks
The theoretical ecologists and applied mathematicians are now working on controlling the algal bloom in ecosystem, but the optimal way to eradicate the algae population to control the algal bloom is yet to know, but in the context of population growth dynamics, how to implement impulsive control strategy to prevent and control the algal bloom with the maximizing economic benefits is worthy of studying.
In the paper, we take an approach to study how to better implement impulsive eutrophication controlling system to control the algal bloom and maintain the normal development of biological treatment process with maximizing economic benefits by using mathematical analysis and simulating the dynamics. On the basis of eutrophication ecology and differential equation, an impulsive eutrophication controlling system is studied analytically and numerically. A key advantage of the impulsive eutrophication controlling system is that it can be quite accurate to describe the interaction effect of some critical factors (fishermen catch and releasing small fry, etc.), which enables a systematic and logical procedure for fitting eutrophication mathematical system to real monitoring data and experiment data. mathematical theoretical works have been pursuing the investigation of two threshold functions 1 ( , 3 ) and 2 ( , 3 ) of two critical parameters and 3 , if 1 ( , 3 ) < 0 and 2 ( , 3 ) < 0, the system (1) remain locally asymptotically stable, and if 1 ( , 3 ) > 0 and 2 ( , 3 ) > 0, the system (1) can ensure all species persistence. These results in turn provide a theoretical basis for the numerical simulation. From the viewpoint of population dynamics, numerical simulation works have been pursuing the investigation of confirming the existence and feasibility of Theorems 7 and 8 and how to choose the best value of two critical parameters and 3 . Within this framework, it can be assumed that the sustainability of the impulsive eutrophication controlling system mainly depends on two critical parameters and 3 . Then the direct and indirect effects caused by critical parameters and 3 are investigated by means of bifurcation analysis. These results indicate that if = 5 and 3 = 0.6, the algae population will go to extinction, but the dominant zooplankton and the filterfeeding fish can survive in a feasible range, which can assure that the eutrophication removal process can be circulatingly developed with maximizing economic benefits. In addition, it is worthwhile to remark that mathematical theoretical results and numerical simulation results are complementary. It is 12 Abstract and Applied Analysis hopeful that these results have an important implication for eutrophication removal management; hence in the follow-up works, we will use the system to character the growth rule and the correlation between the algae population, the dominant zooplankton, and the filter-feeding fish.