Research Article | Open Access

Jianglin Zhao, Yong Yan, "Dynamics of a Seasonally Forced Phytoplankton-Zooplankton Model with Impulsive Biological Control", *Discrete Dynamics in Nature and Society*, vol. 2016, Article ID 2560195, 9 pages, 2016. https://doi.org/10.1155/2016/2560195

# Dynamics of a Seasonally Forced Phytoplankton-Zooplankton Model with Impulsive Biological Control

**Academic Editor:**Darko Mitrovic

#### 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 [5–9]. 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 [12–16]. 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.

**(a)**

**(b)**

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.

**(a)**

**(b)**

**(a)**

**(b)**

**(a)**

**(b)**

**(c)**

**(d)**

**(e)**

**(f)**

#### 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).

#### References

- D. M. Karl, E. A. Law, P. Morris, P. L. B. Williams, and S. Emerson, “Metabolic balance of the open sea,”
*Nature*, vol. 426, no. 11, p. 32, 2003. View at: Google Scholar - J. E. Truscott and J. Brindley, “Ocean plankton populations as excitable media,”
*Bulletin of Mathematical Biology*, vol. 56, no. 5, pp. 981–998, 1994. View at: Publisher Site | Google Scholar - C. S. Holling, “The components of predation as revealed by a study of small-mammal predation of the European pine sawfly,”
*The Canadian Entomologist*, vol. 91, no. 5, pp. 293–320, 1959. View at: Publisher Site | Google Scholar - J. A. Freund, S. Mieruch, B. Scholze, K. Wiltshire, and U. Feudel, “Bloom dynamics in a seasonally forced phytoplankton-zooplankton model: Trigger mechanisms and timing effects,”
*Ecological Complexity*, vol. 3, no. 2, pp. 129–139, 2006. View at: Publisher Site | Google Scholar - T. K. Kar and H. Matsuda, “Global dynamics and controllability of a harvested prey-predator system with Holling type III functional response,”
*Nonlinear Analysis. Hybrid Systems*, vol. 1, no. 1, pp. 59–67, 2007. View at: Publisher Site | Google Scholar | MathSciNet - I. Siekmann and H. Malchow, “Local collapses in the Truscott-Brindley model,”
*Mathematical Modelling of Natural Phenomena*, vol. 3, no. 4, pp. 114–130, 2008. View at: Publisher Site | Google Scholar | MathSciNet - P. Olla, “Effect of demographic noise in a phytoplankton-zooplankton model of bloom dynamics,”
*Physical Review E*, vol. 87, no. 1, Article ID 012712, 2013. View at: Publisher Site | Google Scholar - I. Bashkirtseva and L. Ryashko, “Stochastic bifurcations and noise-induced Chaos in a dynamic prey-predator plankton system,”
*International Journal of Bifurcation and Chaos*, vol. 24, no. 9, Article ID 1450109, 7 pages, 2014. View at: Publisher Site | Google Scholar - R. K. Upadhyay, “A solution to the evolution-related Truscott-Brindley model for the generalized phytoplankton-zooplankton populations,” http://arxiv.org/abs/1212.5420. View at: Google Scholar
- M. Gao, H. Shi, and Z. Li, “Chaos in a seasonally and periodically forced phytoplankton-zooplankton system,”
*Nonlinear Analysis. Real World Applications*, vol. 10, no. 3, pp. 1643–1650, 2009. View at: Publisher Site | Google Scholar | MathSciNet - J. Luo, “Phytoplankton-zooplankton dynamics in periodic environments taking into account eutrophication,”
*Mathematical Biosciences*, vol. 245, no. 2, pp. 126–136, 2013. View at: Publisher Site | Google Scholar | MathSciNet - H. Yu, M. Zhao, Q. Wang, and R. P. Agarwal, “A focus on long-run sustainability of an impulsive switched eutrophication controlling system based upon the Zeya reservoir,”
*Journal of the Franklin Institute*, vol. 351, no. 1, pp. 487–499, 2014. View at: Publisher Site | Google Scholar - J. Yang and M. Zhao, “A mathematical model for the dynamics of a fish algae consumption model with impulsive control strategy,”
*Journal of Applied Mathematics*, vol. 2012, Article ID 452789, 17 pages, 2012. View at: Publisher Site | Google Scholar | MathSciNet - L. Wang, L. Chen, and J. J. Nieto, “The dynamics of an epidemic model for pest control with impulsive effect,”
*Nonlinear Analysis: Real World Applications*, vol. 11, no. 3, pp. 1374–1386, 2010. View at: Publisher Site | Google Scholar | MathSciNet - H. Baek, Y. Do, and Y. Saito, “Analysis of an impulsive predator-prey system with Monod-Haldane functional response and seasonal effects,”
*Mathematical Problems in Engineering*, vol. 2009, Article ID 543187, 16 pages, 2009. View at: Publisher Site | Google Scholar | MathSciNet - P. Georgescu and H. Zhang, “An impulsively controlled pest management model with n predator species and a common prey,”
*Bio Systems*, vol. 110, no. 3, pp. 162–170, 2012. View at: Publisher Site | Google Scholar - Z. Teng, “Uniform persistence of the periodic predator-prey Lotka-Volterra systems,”
*Applicable Analysis*, vol. 72, no. 3-4, pp. 339–352, 1999. View at: Publisher Site | Google Scholar | MathSciNet - H. Guo and X. Song, “An impulsive predator-prey system with modified Leslie-Gower and Holling type II schemes,”
*Chaos, Solitons and Fractals*, vol. 36, no. 5, pp. 1320–1331, 2008. View at: Publisher Site | Google Scholar | Zentralblatt MATH - S. Rinaldi, S. Muratori, and Y. Kuznetsov, “Multiple attractors, catastrophes and chaos in seasonally perturbed predator-prey communities,”
*Bulletin of Mathematical Biology*, vol. 55, no. 1, pp. 15–35, 1993. View at: Publisher Site | Google Scholar

#### Copyright

Copyright © 2016 Jianglin Zhao and Yong Yan. 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.