Abstract

This paper investigates the dynamics of a seasonally forced phytoplankton-zooplankton model with impulsive biological control. It shows that the periodic eradicated solution is unstable. Further, the condition for permanence of the system is established by relations between the model parameters and the intensity of the impulses. The numerical analysis is performed to study the effect of seasonality and impulsive perturbations on plankton dynamics. The numerical results imply that the seasonal forcing can trigger more periodic mode and the impulsive period for control of the size of phytoplankton is more practicable to the system than the impulsive release of zooplankton. These conclusions provide a better understanding of controlling harmful algae blooms.

1. Introduction

Algae are very diverse and found almost everywhere on the planet. They play an important role in marine, freshwater, and some terrestrial ecosystems. On one hand, phytoplankton are critical to sustain most aquatic food chains and produce half of the world’s oxygen in the process of photosynthesis [1]. On the other hand, phytoplankton populations can grow explosively and lead to severe oxygen depletion in the relevant waters when harmful algal blooms occur. As a consequence, human activities would be limited and economy would suffer. Therefore, in order to prevent and control harmful algae blooms a better understanding of the mechanisms that trigger the occurrence or explain the absence of a phytoplankton bloom is of considerable significance. Truscott and Brindley [2] were the first to offer a model considering the prey-predator system of phytoplankton and zooplankton as a nonlinear excitable system to explain the dynamics of harmful algal blooms. In the framework, the evolution of phytoplankton and zooplankton populations is formulated by where and represent population densities of phytoplankton and zooplankton, respectively. The first term in the right side of the first equation in (1) is the logistic growth function. The term is called Hollings type III grazing function [3]. is the maximum growth rate of phytoplankton when is small, is the environmental carry capacity of phytoplankton, reflects the maximum specific predation rate, governs how quickly that maximum is attained as prey densities increase, denotes the ratio of biomass consumed to biomass of new herbivores produced, and measures the zooplankton death rate. If we quote the typical parameter values [2, 4] for system (1), then the positive equilibrium of (1) is asymptotically stable. Assume that the initial phytoplankton concentration is fixed in the stationary value. Then, a phytoplankton bloom that is followed by a delayed zooplankton bloom is triggered by suppressing the initial zooplankton concentration sufficiently far below the stationary value [4].

As a classical predator-prey model, it has been extensively discussed by many researchers [59]. In fact in spite of noticing the remarkable impact of seasonal forcing on the phytoplankton birth-rate by Truscott and Brindley, there are few models explicitly taken into account. Even though Gao et al. [10] had considered the effect of seasonality and periodicity on the growth rate of phytoplankton in their model, the growth rate taken as a sinusoidal forcing function of time could not reflect sufficiently the explosive growth of phytoplankton. Freund et al. [4] filled the gap, presenting and discussing simulation results of the seasonally forced Truscott-Brindley model with an exponential function depicting the explosive growth. Later, Luo [11] developed the seasonally forced Truscott-Brindley model by including the growth rate and the intrinsic carrying capacity of phytoplankton changing with respect to time and nutrient concentration. For simplicity we only consider the maximum growth rate in the model (1) as periodically varying function of time due to seasonal variation. And, we adopt the same relation between the phytoplankton growth rate and seasonally varying temperature as Freund et al. [4] have suggested. The effect of changing temperature on growth rate of phytoplankton is where is the average value of the intrinsic growth rate of phytoplankton, is assumed to be a constant which asserts that a change of the temperature by will multiply the rate at mean temperature, denotes time (days), represents temperature on time which is adapted from a fit by using average temperature data, is the average temperature in one cycle, denotes the amplitude of temperature, is the angular velocity, and is the initial phase. For the sake of simplification, the intraspecific competition of phytoplankton is not affected by seasonally varying temperature.

The use of impulsive control for ecological systems is proved to be one of the most effective methods and has received much attention from both ecologists and applied mathematicians [1216]. However, almost all the work on the Truscott-Brindley models neglects the impulsive biological control of phytoplankton. Thus, we periodically release zooplankton in laboratories at a constant to reduce the population level of phytoplankton by grazing. Keeping these aspects in view, we establish the following model:where , , is the period of the impulsive effect, denotes the concentration change of zooplankton by releasing which is determined by the maximum amount of zooplankton produced by laboratories and , and is the set of all nonnegative integers. In system (3), all parameters are supposed to be positive constants.

This paper is organized as follows. In Section 2, we study the dynamics of system (3) without impulsive effect. In Section 3, we mainly focus on system (3) and obtain the stability of phytoplankton-eradication periodic solution and the condition of permanence of (3). In Section 4, the numerical analysis is performed to investigate the dynamics of (3). Finally, we close with a discussion in Section 5.

2. Dynamical Properties of (3) without Impulsive Effect

In order to analyze to the dynamical behavior of (3) without impulsive effect, we first consider a periodically logistic equation:where and are periodically continuous functions defined on with the common period . According to the results of [11, 17], we obtain the following conclusions.

Lemma 1 (see [17]). If for and , then (4) has a unique nonnegative -periodic solution which is globally asymptotically stable. That is, , for each positive solution of (4). Moreover, if , then for and if , then .

Consider the extended logistic model: Since , , and are positive constants and is -periodic, by means of Lemma 1, there are the following results.

Theorem 2 (see [17]). System (5) has a unique positive -periodic solution which is globally asymptotically stable; that is, there is a positive -periodic function such that for each solution of (5) with positive initial value , one has .
Consider the following system:where , , , , , , and are positive constants and is a periodic continuous function on with the common period .

Definition 3. System (6) is said to be permanent if there exist constants , , satisfying and when , where is any solution of system (6) with the initial values and .

Theorem 4 (see [11]). System (6) is permanent provided thatwhere is the uniquely periodic solution of (5) established by Theorem 2.

Theorem 5 (see [11]). Assume that then population of zooplankton tends to extinction; more precisely we havefor each solution of (6), where is the uniquely periodic solution of (5) established by Theorem 2.

3. Extinction and Permanence of (3)

Definition 6. System (3) is regarded as permanent if there exist constants , , satisfying and when , where is any solution of system (3) with the initial values and .

Theorem 7. Each solution of system (3) is ultimately bounded. That is, there exists a constant , such that and for each solution of system (3) with all large enough.

Proof. Suppose that is a solution of system (3). Let . Computing the upper right derivative of for ,Clearly, there exists satisfying the following:where .
For , .
Then, if , we can getHence, . is ultimately bounded by a constant and there exists a constant such that and for each solution of system (3) with all large enough. This completes the proof.

Theorem 8. System (3) admits a positive periodic solution that corresponds to phytoplankton eradication:where . Furthermore, .

Proof. If the phytoplankton population is absent, that is, , system (3) reduces to Then it yieldsThen, , .
Imposing , we obtain the periodicity condition , from which (13) is clear.
Note that the solution of (14) with initial value is gotten by Hence, .

Theorem 9. Let be a solution of (3). Then is unstable.

Proof. To investigate the local stability of periodic solution , we will use the method of small amplitude perturbations. To this purpose, define where and represent the small amplitude perturbations.
Therefore, the linearization of system (2) becomesThen, it results inwhere satisfiesand , the identity matrix.
The linearization of the third and the fourth equations of (3) becomesThe stability of periodic solution is determined by the eigenvalues of which areThus, according to Floquet theory [16], is locally stable if and . With the assumption of (3), we know that . So, . This establishes the theorem.

Theorem 10. System (3) is permanent if , where .

Proof. Suppose that is a solution of (3) with . From Theorem 7, assume that and for . From Theorem 8, we have for all large enough and some . Let , so for large enough. Therefore, we will next find such that for large enough. We will do it in the following two steps.
Step  1. Put . Since , we can select and small enough such thatWe will show that there exists such that . Otherwise, using the above hypothesis, we get . By Lemmas  2.2 in [18], we have and , where is the solution of the following equation:and , . Thus, there exists such that for .
Integrating (26) on , it getsThen, , , which is a contradiction to the boundedness of . Hence, there exists such that .
Step  2. If for all , then our goal is accomplished. If not, we set . Then, for . Since is continuous, we have and . By Step  1, there also exists such that . Put . Then, for and . Let this process continue by using Step  1. If the process is stopped in finitely many times, the theorem follows. Otherwise, define an interval sequence with such that , and , . Put . If , we can obtain a subsequence satisfying , . As shown in Step  1, it will lead to a contradiction because of the boundedness of . Thus, is true. Note that , where . Put and . It follows for sufficiently large. The proof is completed.

4. Numerical Analysis

In this section, we will study the influence of seasonal forcing and impulsive perturbations. We quote typical parameter values [4]: , μg/L, , μg/L, , , , , , and . Based on Theorems 2 and 4, we see that the sufficiently small value of leads to be permanent for (6). For enough large value of , in contrast, (6) possesses a periodic solution that corresponds to zooplankton eradication. Figure 1(b) illustrates our conclusion. By comparison with Figure 1, a seasonal forcing of the phytoplankton growth rate can trigger a bloom mode. This claim is pointed out to be important for considering the seasonal variation in study of the Truscott and Brindley model.

First, we analyze the effect of on the long-run dynamics of system (3). For fixed parameter values and days, the bifurcation diagrams with respect to have been plotted in Figure 2(b). It needs to highlight that with the change of the parameter value system (3) undergoes period-3 attractors and later enters into chaotic behaviors. We know that system (3) without impulsive release of zooplankton presents the harmful algae bloom mode (Figure 1(b)). However, it can be seen in Figure 2(b) that if the impulsive biological control is considered, then the nonbloom mode will occur. That is to say that the occurrence of phytoplankton blooms is successfully prevented. Besides, for the purpose of studying the influence of seasonal forcing of the phytoplankton growth rate on system (3), the bifurcation diagrams of system (3) without considering seasonality are shown in Figure 2(a). As we can see in Figure 2(a), system (3) without seasonal fluctuations in the growth rate of phytoplankton exhibits relatively simple dynamical behaviors: period-1 attractor. It can be observed in Figure 2 that adding seasonal forcing reduces the maximum density of phytoplankton and enhances the zooplankton population size. What is more is that adding seasonal forcing produces chaos. It is noted in Figure 2(b) that although system (3) shows chaotic dynamics, the phytoplankton population does not trigger bloom mode. Second, we will investigate the influence of impulsive period on the population dynamics of system (3). Thus, we set and μg/L. It follows that system (3) shown in Figure 3(a) exhibits the only dynamical behavior: period-1 attractor when we do not consider seasonality. It is worth to be stressed that the max density of phytoplankton is monotonically increasing with increase of , and the max density of zooplankton is monotonically decreasing. In Figure 3(a) we need to emphasize that changing of system (3) without seasonality does not exhibit the phytoplankton bloom. Conversely, if the seasonal factor is included, Figure 3(b) illustrates relatively complex dynamical behavior for system (3). Figures 4 and 5 display the detailed cases. It is worth pointing out that in Figure 4 when the time of releasing presents a 5-day delay, the periodic solution inflates to until . If days, the phytoplankton population of system (3) exhibits bloom dynamics (Figure 5). By comparison with Figure 4, it concludes that a seasonal forcing augments fluctuation of population. Further analysis shows that is more practicable to system (3) as a control parameter than . Thus, we suggest that in order to control the size of phytoplankton the impulsive period should be selected to be feasible.

5. Discussion

We have studied the dynamics of a seasonally forced phytoplankton-zooplankton model with impulsive biological control. By control of the releasing period of zooplankton, the size of phytoplankton could be reduced. However, the amount of is determined by laboratories and this kind of production is limited. Therefore, system (3) cannot be used for larger areas of water. We suggest that system (3) should be suitable for smaller areas of water such as urban lakes. In this paper we consider the dynamic behaviors of the system when periodic forcing is solely superimposed on the growth rate of the phytoplankton. In fact, there are a number of ways to apply periodic forcing in an ecological model [19]. Thus, in the future, applying periodic forcing to the Truscott and Brindley model with impulsive biological control will be considered. From the above, it is known that releasing zooplankton to control the size of phytoplankton is limited. Therefore, we need study other effective ways to reduce harmful algae blooms such as introducing competitive population for phytoplankton and harvesting for phytoplankton.

Competing Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (no. 11461058) and by the Scientific Research Fund of Sichuan Provincial Education Department (nos. 14ZA0296 and 15ZB0333).