Discrete Dynamics in Nature and Society

Volume 2016 (2016), Article ID 2560195, 9 pages

http://dx.doi.org/10.1155/2016/2560195

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

Department of Mathematics, Sichuan Minzu College, Kangding 626001, China

Received 27 May 2016; Accepted 22 June 2016

Academic Editor: Darko Mitrovic

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.

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