Research Article | Open Access

Chandrima Banerjee, Pritha Das, "Impulsive Control on Seasonally Perturbed General Holling Type Two-Prey One-Predator Model", *Discrete Dynamics in Nature and Society*, vol. 2016, Article ID 9268257, 20 pages, 2016. https://doi.org/10.1155/2016/9268257

# Impulsive Control on Seasonally Perturbed General Holling Type Two-Prey One-Predator Model

**Academic Editor:**Zhengqiu Zhang

#### Abstract

We investigate the dynamical behaviors of two-prey one-predator model with general Holling type functional responses. The effect of seasonal perturbation on the model has been discussed analytically as well as numerically. The periodic fluctuation is considered in prey growth rate and the predator mortality rate of the model. The impulsive effects involving biological and chemical control strategy, periodic releasing of natural enemies, and spraying pesticide at different fixed times are introduced in the model with seasonal perturbation. We derive the conditions of stability for impulsive system using Floquet theory, small amplitude perturbation skills. A local asymptotically stable prey (pest) eradicated periodic solution is obtained when the impulsive period is less than some critical value. Numerical simulations of the model with and without seasonal disturbances exhibit different dynamics. Also we simulate numerically the model involving seasonal perturbations without impulse and with impulse. Finally, concluding remarks are given.

#### 1. Introduction

Functional response plays an important role in dynamics of population ecology. There are several examples in nature where the predator population affects the number of their prey populations. When predators are faced with increasing local density of their prey, they often respond by changing their consumption rate. This relationship of an individual predatorâ€™s rate of food consumption with the prey density was first suggested by Holling [1] as functional response. Functional responses are generally classified into Holling types I, II, III, and IV. Type I is the simplest capture rate which increases in direct proportion to prey density until it abruptly saturates. Similarly in type II, the rate of capture increases with increasing prey density. In contrast to the linear increase of type I, type II approaches saturation gradually at higher prey density. This behavior is better expressed by an asymptotic relationship between per capita feeding rate and prey density which is of the functional form , where is the prey density. (per unit time) and (per unit prey) are positive constants describing the effects of capture rate and handling time on the feeding rate of the predator, respectively. Holling type II functional response is most common type of functional response among arthropod predators. Type III resembles type II in having an upper limit to prey consumption except at low prey density, and the predatorâ€™s response to prey is depressed. The mathematical expression of type III functional response is of the form , where is half-saturation constant. For example, at low prey density, predators switch to alternative prey species, if the focal prey is less accessible due to surplus refuges [1]. Holling types I, II, and III are all monotonically nondescending in the first quadrant [2]. But for some prey-predator systems, a nonmonotonic response (Holling IV) occurs because the growth of predator may be inhibited when the prey density reaches a high level [3]. Holling type IV functional response is written as This is also known as Monod-Haldane-type functional response [4], where as the denominator of the expression does not vanish for nonnegative . When , a simplified form is proposed by Sokol and Howell [5], and some researchers also called it Holling type IV [4, 6]. In population dynamics, this functional response describes the phenomenon of group defense which causes a decrease or prevention of predation due to the increased ability of the prey to better defend or disguise themselves when their numbers are large enough. Lone ox can be attacked successfully by wolves in regular basis. While small herds of musk ox grazing in group are attacked by wolves with less success rate. But successful attacks on larger herds are not observed in general [7] (see Tener, 1965).

Kooij and Zegeling [8] and Sugie et al. [9] considered the predatorâ€™s functional response term of the form , where is a positive constant. is the conversion rate of prey captured by predator. This particular type functional response with was introduced by Kazarinoff and van den Driessche [10]. This is known as general Holling type functional response.

A model for two-prey one-predator ecological system with general Holling type functional responses is discussed here. In this paper, apparent competition between the prey is treated. Only the predator hunting time for each prey depending on the selection and availability of its favorite prey is considered. Most of the models are considered with constant environment which is not ecologically realistic. There are many factors in the environment which are not constant but vary periodically with change in time such as concentration of atmospheric carbon dioxide, climate, and ecological succession. It is natural to identify that seasonality plays an important functional role in changing the behavior of individual population. In the last decades, many scientists [11â€“15] have studied the interactions between seasonality and internal biological rhythms of simple predator-prey systems. Without seasonal variations the two-species nonlinear autonomous dynamical system has either stable or periodic solution. Seasonality introduces complex dynamical behavior in the system such as the existence of multiple attractors, catastrophes, and chaos [15â€“19]. In this paper we investigate how two-prey one-predator model with general Holling type functional responses changes its dynamics with change in seasonal parameters like degree of seasonality. The periodic variations in the parameters may not be synchronous due to seasonality. Seasonality parameters with different phases reach their maximum influence at different times [19]. We consider seasonal variation in prey growth rate and predator mortality rate with different phase angles between them. We investigate the boundedness of the model and the condition of uniform persistence.

The effects of impulsive perturbation in population ecology have been widely studied and discussed by number of researchers [20â€“26]. Impulsive perturbations bring sudden changes in the system. Ecological model with combined effect of impulsive and seasonal perturbation has been investigated by few researchers [27â€“30]. Pei et al. [20] have considered also two prey and one predator with group defense and Holling type II functional responses and investigated the effect of impulsive perturbation in pest control using pesticides. In this paper, we study general Holling type two prey (pests) and one predator with seasonal disturbances as well as proportional periodic impulsive poisoning (spraying pesticide) for all species and constant periodic releasing, or immigrating, for the predator at different fixed time. We study the local stability of the prey- (pest-) free periodic solution.

The rest of this paper is organized as follows. In Section 2, a detailed description about mathematical formulation of model with seasonal variation is given. Then the model is modified by incorporating impulsive perturbation in it. In Section 3, we discuss various important theorems and lemmas. In Section 4, the boundedness and uniform persistence of the solution of the system are proved. Expression for prey-free periodic solution of the system is obtained in Section 5. We derive the conditions for stability of the prey- (pest-) free periodic solution. In Section 6, results of numerical analysis are presented. In Section 7 the results are discussed thoroughly.

#### 2. Mathematical Model Formulation

##### 2.1. Model Formulation without Impulsive Effect

We consider the following deterministic model with general Holling type function responses and seasonally perturbed parameter is described by the following system of differential equations: where , , and are the biomass of the two prey and the predator at time , respectively. s are intrinsic birth rates of prey and is intrinsic mortality rate of predator, s are parameters representing half-saturation constants, and s are equal to the transformation rates of predator due to predation (rates of converting prey into predator). s are the per capita predation rates of the predator. All parameters are considered as positive constants . The parameters s represent the degree of seasonality where . , , and are the magnitudes of the perturbations in , , and , respectively [18]. is the angular frequency of the fluctuations caused by seasonality. Finally and , where and also , can be interpreted as the difference in phase angles between the seasonality of the above three parameters. Clearly, when which imply that there is in-phase synchronous variation in between the intrinsic growth rates of prey, while , the variation is antiphase synchronous. In this paper, the values of phase angles will be considered , and and .

##### 2.2. Model Formulation with Impulsive Effect

We develop system (1) by introducing a proportional periodic impulsive poisoning (spraying pesticide) for all species and a constant periodic releasing, or immigrating, for the predator at different fixed times as follows: where . Here , , and are considered as the biomass of the two prey (pest) and the predator (natural enemy) at time , respectively. () present the fraction of the prey and the predator which die due to the harvesting or pesticides and so forth. is the period of the impulsive immigration or stock of the predator. is the size of immigration or stock of the predator. In order to get some conditions guaranteeing the system is permanent, is the period of spaying pesticides (harvesting) and the impulsive immigration or stock of the predator, respectively. We present the fraction of the prey and the predator which die due to the harvesting or pesticides and so forth.

Incorporating (see [31]) in system (2), we get

#### 3. Preliminaries and Basic Lemmas

*Definition 1. *System (1) is uniformly persistent if, for every positive solution of system (1), there exist positive constants , , , , , , and such that , , and for .

Lemma 2 (see [32]). *If , , and , then for and one has .*

*Note 1. *From Lemma 2 we get

Let , . Let be the set of nonnegative integers and let be defined as the map by the right hand of the system of equations (1). Let ; then belongs to class if(1)is continuous in , and for every , and exist;(2) is locally Lipschitzian in .

*Definition 3. *Let ; then for and , the upper right derivative of with respect to the impulsive differential system (2) is defined as

The solution of system (2) is a piecewise continuous function , is continuous on , where and , and then and exist. Thus these properties guarantee the global existence of and uniqueness of solution of system (2) [33, 34].

*Definition 4. *System (2) is permanent if there exists a compact such that every solution of system (2) will enter and remain in the region .

The following lemma is evident.

Lemma 5. *Let be a solution of system (2) with ; then for all . Also , if .*

We will use a comparison result of impulsive differential inequalities. Suppose that satisfies the following hypotheses.(H) is continuous on and the limit exists, where and , and is finite for and .

Lemma 6 (see [33]). *Suppose and where satisfies (H) and are nondecreasing for all . Let be the maximal solution for the impulsive Cauchy problem defined on . Then, implies that , , where is any solution of (6).*

Lemma 7. *There exists a constant such that , , and for each solution of system (2) with all large enough.*

We use comparison theorem on impulsive differential equation: Using the lemma, we obtain the following result.

System (2) has a unique positive periodic solution which is globally asymptotically stable, where, where is set of natural numbers and is a positive periodic solution of system (7). We obtain that

Lemma 8. * for all positive periodic solution of system (2) with .*

Therefore, system (2) has a prey- (pest-) free periodic solution Next we find the expression of some critical parameters under the condition for all species persistence of system (1).

#### 4. Boundedness and Uniform Persistence

Theorem 9. *There exist position constants , and for every position solution of system (1) with all sufficiently large t.*

*Proof. *Assume that is an arbitrary positive solution of system (1); then from the first and second equations of system (1) we getWe have Similarly, We have Let Then, since , where According to the theorem on differential inequalities [35], we have and, for , Thus, all solutions of system (1) enter into the region This proves the theorem.

*Note 2. *Since is bounded above, all of are bounded above. Therefore there exist positive constants , and such that , , and for .

Theorem 10. *Ifâ€‰â€‰, and for every small constant , then system (1) is uniformly persistent.*

*Proof. *From Theorem 9, there exist positive constants , and such that for From the first equation of system (1), we get From Lemma 2, if , then there is positive constant such that for . From the second equation of (1) we obtain that It is known that if , there are two positive constants and such that for . From the third equation of system (1), we get . We know that if , there are two positive constants and such that for . Let ; then we have for . This completes the proof.

#### 5. Stability of a Periodic Solution with Prey Eradication

In this section, we examine the stability of the pest eradication periodic solution of system (2). We find stability analysis of system (2) for different values of , , , and . We have system (2) that represents two-prey and one-predator model, (i) with Holling type II functional responses for , (ii) with Holling type III functional responses for , (iii) with Holling type IV functional response for , , , , (iv) with mixed Holling type II-III functional responses for and , (v) with mixed Holling type IIâ€“IV functional responses for and , and (vi) with mixed Holling type III-IV functional responses for and . We find that Holling type II, Holling type IV, and mixed Holling types IIâ€“IV system have same Jacobian matrix for the solutions () and Holling type III, mixed Holling types II-III, and Holling type III-IV systems have different Jacobian matrices for the solutions ().

Theorem 11. *Let be any solution of system (3) with Holling type II, Holling type IV, and mixed Holing IIâ€“IV functional responses; then is locally asymptotically stable if for where , .*

*Proof. *The local stability of the periodic solution may be determined by considering the behavior of small amplitude perturbations of the solution. Define , , and .

From system (3), we have the following linear system: where satisfies and , where is the identity matrix. So the fundamental solution matrix isThe resetting impulsive conditions of system (28) become Note that all eigenvalues of the matrix are , , and .

So , if and only if where , .

According to Floquet theory (see [34]) of impulsive differential equation, the two-pest eradication solution is locally asymptotically stable. This completes the proof.

Theorem 12. *Let be any solution of (3) with Holling type III; then periodic solution is locally asymptotically stable if*

*Proof. *Let us take the transformation ; then the system becomes where satisfies Using Floquet theory (see [34]), we find that the following inequalities hold: Therefore, is locally asymptotically stable. The proof is completed.

Theorem 13. *The periodic solution of system (3) with mixed Holling type II-III functional responses is locally asymptotically stable ifwhere .*

*Proof. *Let us consider the transformation ; then the liberalization of the system is where satisfies Applying Floquet theory (see [34]), the eigenvalues are , , and .

By using Floquet theory, if and only if where , . Therefore is locally asymptotically stable if The proof is completed.

Theorem 14. *The periodic solution of system (3) with mixed Holling type III-IV functional responses is locally asymptotically stable if where .*

*Remark 15. *If we consider the condition (see [31]), then similar results can be obtained by the above procedure by replacing with ().

#### 6. Numerical Simulation and Results

In this section we numerically simulate systems (1) and (2) to investigate the changing dynamics of the systems due to seasonal variation and impulsive perturbation using Matlab R2012a. Firstly, the effects of seasonality in the system will be studied without impulsive perturbation. Then the combined influences of seasonal disturbances and impulsive perturbation on system (2) are discussed. The biologically feasible parametric values are given below:The initial conditions are .

The degree of seasonality (), angular frequency of the fluctuations caused by seasonality, and difference in phase angle are varied to find the effect of seasonal disturbances on system (1). Also the impulsive parameter is gradually increased to observe the effect of periodic impulsive release of population in system (2).

##### 6.1. Effect of Seasonal Disturbances on the Model without Impulse

From the numerical simulation we find that Holling type III, type IV, mixed type II-III, and mixed type III-IV systems are initially stable but exhibit periodic oscillations in presence of seasonal variations, whereas Holling type II and mixed type IIâ€“IV system dynamics change from limit cycle to period two and quasiperiod oscillations, respectively.

Small and large amplitude fluctuations can be found depending upon the intensity of the disturbances (mainly degree of seasonality and difference in phase angles between the seasonal parameters). It is observed that biomass level of the species increases significantly in presence of seasonality. It is further observed that there is no possibility of extinction and all populations can coexist in seasonally varying environment. The population time series shows the asymptotic stable behavior for Holling type III system in absence of seasonality (see Figure 1(a)). The solution trajectory exhibits oscillatory behavior in presence of seasonal disturbances (see Figure 1(b)). The population time series shows periodic behavior for mixed Holling type IIâ€“IV system in absence of seasonal variations (see Figure 2(a)) while quasiperiodic behavior is found in presence of seasonal variations (see Figure 2(b)).

**(a)**

**(b)**

**(a)**

**(b)**

and show in-phase, outer phase, and antiphase synchronous oscillation (see Figure 3) for different values of . The effect of changing angular frequency in the dynamics of the system has been analyzed. Figure 4 depicts the local long term dynamics for , , when varies from 0 to 0.4. If increases from 0 to 0.09, the system has rich changes in dynamics and shows periodic windows, period-doubling, chaos, and a chaotic crisis (see Figure 4(a)). Figures 4(b) and 4(c) show same pattern when and , respectively. These figures show that the size of the chaotic attractor abruptly changes, exhibiting attractor crisis. Hence, it is observed that the seasonality has a strong effect on the dynamical behaviors of system (1). The bifurcation diagrams of system (1) are plotted involving , , and populations with respect to different seasonal parameters like and . The resulting bifurcation diagrams clearly show that system (1) has rich dynamics including periodic oscillation, period-doubling bifurcation, period-halving bifurcation, and chaos (see Figures 5 and 6).

**(a)**

**(b)**

**(c)**

**(a)**

**(b)**

**(c)**

**(a)**

**(b)**

**(a)**

**(b)**

**(c)**

##### 6.2. Combined Effect of Impulsive and Seasonal Perturbations

Here we examine how the period of release natural enemies affects the dynamical behavior of system. In general, it is found that when the values of are relatively high and is in low range , and (pests) population biomass grow whereas (natural enemy or predator) biomass significantly declines (see Figure 7). For moderate value of (), population (second pest) biomass attains lower values and population biomass does not increase significantly while population (first pest) increases relatively to a larger value (see Figure 8). But for excessive high values of (), and population biomass go to extinction. It is also observed that biomass density increases up to a level and then fluctuates within a certain range even if the value of is increased further (see Figure 9).

**(a)**

**(b)**

**(c)**

**(a)**

**(b)**

**(c)**

**(a)**

**(b)**

**(c)**

Mixed Holling type II-III system shows a slightly different result as compared to other Holling type systems. Initially three species are densely populated for (see Figure 10) but for increasing values of population goes into extinction, population biomass fluctuates, and biomass density increases (see Figure 11). For , population biomass shows an decreasing trend while population becomes extinct and population biomass increases (see Figure 12). The bifurcation diagrams with respect to illustrates quasiperiodic, periodic oscillation and stable behavior when varies from to . (see Figure 13).

**(a)**

**(b)**

**(c)**

**(a)**

**(b)**

**(c)**

**(a)**

**(b)**

**(c)**

**(a)**

**(b)**

**(c)**