Next Article in Journal
Theoretical Validation of New Two-Dimensional One-Variable-Power Copulas
Next Article in Special Issue
Solutions of (2+1)-D & (3+1)-D Burgers Equations by New Laplace Variational Iteration Technique
Previous Article in Journal
Stability Switching in Lotka-Volterra and Ricker-Type Predator-Prey Systems with Arbitrary Step Size
Previous Article in Special Issue
Application of the Optimal Homotopy Asymptotic Approach for Solving Two-Point Fuzzy Ordinary Differential Equations of Fractional Order Arising in Physics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamics of a Double-Impulsive Control Model of Integrated Pest Management Using Perturbation Methods and Floquet Theory

by
Fahad Al Basir
1,
Jahangir Chowdhury
2 and
Delfim F. M. Torres
3,*
1
Department of Mathematics, Asansol Girls’ College, Asansol 713304, India
2
Department of Applied Science, RCC Institute of Information Technology, Canal South Road, Kolkata 700015, India
3
Center for Research and Development in Mathematics and Applications (CIDMA), Department of Mathematics, University of Aveiro, 3810-193 Aveiro, Portugal
*
Author to whom correspondence should be addressed.
Axioms 2023, 12(4), 391; https://doi.org/10.3390/axioms12040391
Submission received: 20 January 2023 / Revised: 31 March 2023 / Accepted: 15 April 2023 / Published: 18 April 2023
(This article belongs to the Special Issue Mathematical Methods in the Applied Sciences)

Abstract

:
We formulate an integrated pest management model to control natural pests of the crop through the periodic application of biopesticide and chemical pesticides. In a theoretical analysis of the system pest eradication, a periodic solution is found and established. All the system variables are proved to be bounded. Our main goal is then to ensure that pesticides are optimized, in terms of pesticide concentration and pesticide application frequency, and that the optimum combination of pesticides is found to provide the most benefit to the crop. By using Floquet theory and the small amplitude perturbation method, we prove that the pest eradication periodic solution is locally and globally stable. The acquired results establish a threshold time limit for the impulsive release of various controls as well as some valid theoretical conclusions for effective pest management. Furthermore, after a numerical comparison, we conclude that integrated pest management is more effective than single biological or chemical controls. Finally, we illustrate the analytical results through numerical simulations.

1. Introduction

In today’s farming systems, a variety of approaches are used for pest control. Maintaining high output while guaranteeing sustainability is crucial for the entire agriculture sector. Since the beginning of human civilization, insect and pest control has been one of the most significant difficulties in the agricultural sector [1,2]. Every day, people come up with fresh ideas for equipment and tactics to use in their fight against pests. As a result of human efforts to manage pests, our natural ecology and nature are on the verge of extinction [3].
Chemical controls are less expensive to implement, yet they result in significant environmental damage [4,5,6]. On the other hand, biological controls are more costly to implement but have less environmental impact [7]. However, frequently, the use of a single control method is not beneficial to control pest resistance and preserve environmental quality [8,9]. In order to reduce insect populations below economic levels, integrated pest management (IPM), a safer and more effective method, was developed. IPM is used for a variety of agronomic crops and is now widely used as an economical and environment-friendly pest control method in several nations [8,9,10]. When the ecological cost of management is added to the economic price of controls, a combination of chemical and biological controls yields a superior result when they are used with proper rate and with tolerable intervals. Modeling of this phenomenon leads naturally to the use of impulsive differential equations [11,12].
Many researchers have designed mathematical models for pest management through control strategies, some of which promote chemical agents [4,7], some advocate the use of biological agents to impose a total solution of pest and disease [1,2,10,13,14,15], and some researchers use both the chemical and biopesticides in their mathematical models [5,16,17,18]. Mathematical-model-based works using impulsive differential equations are also available in the literature, as already mentioned [19,20,21,22,23,24,25,26]. Recently, Li, Huang, and Liu proposed a pest management model to simulate the application of pesticides and build a pesticide function with residual and delayed effects of pesticides, proposing pest management with pollutant emission [27]. Liu et al. constructed a mathematical model for pest control in which susceptible and infected pests are separated from the pest population and only susceptible pests are harmful to crops [28]. They weighed the two approaches of spraying pesticides and releasing diseased pests and natural enemies to control vulnerable pests when completing their task. Alzabut established a mathematical model based on the sense of biological survey in the field of agriculture, and introduced various control methods to determine how to protect the crops from destructive pests [29]. In [24], an integrated pest management model using impulsive differential equations was proposed and analyzed for Jatropha curcas using the release of infective pests and spraying of chemical pesticides. The existence and stability of susceptible pest-eradication solutions were analyzed using Floquet theory and the small amplitude perturbation method. To the best of our knowledge, all available articles deal with single-impulse differential equation models where the stability analysis of the periodic pest extinction solution is obtained by the Floquet theory, the method of small amplitude perturbation, and the comparison theorem. However, none of the prior research available in the literature employs the concentration of chemical pesticides as a system variable as we do here. Moreover, in our study, we spray biological and chemical pesticides at two different time intervals, simultaneously varying the time period.
Pest control models using a single impulse are available (see, e.g., [30]), but using double-impulsive controls is rare [31,32]. The authors of [31] took a predator population along with biopesticides in an impulsive periodic way for the control of crop pests. The authors of [32] proposed a predator–prey model with disease in the prey and investigated it for the purpose of integrated pest management. The permanence of the system and global stability of the susceptible pest-eradication periodic solution were shown by means of the released amounts of infective prey and predator. In contrast, here, impulses on both chemical and biopesticides were assumed in the formulation of the mathematical model for crop pest management.
Our use of the concentration profile of the chemical pesticide as a model variable is a novel approach. We demonstrated the dynamics using both the chemical and biological pesticides in the system in an impulsive way, which is, to the best of our knowledge, a novel concept in crop pest control. The proposed double-impulsive system was analyzed with proper analytical methods, namely, using Floquet theory and the perturbation method.
Floquet theory is a powerful mathematical tool for analyzing periodic systems, and it can be extended to impulsive models with periodic impulses. In impulsive models, the system’s behavior is characterized by a sequence of discrete impulses applied at regular intervals. These impulses may arise in many practical scenarios, including electrical circuits, control systems, and biological systems. Floquet theory provides a robust framework for analyzing and designing control strategies for impulsive models with periodic impulses. One can analyze the stability of the impulsive system by examining the eigenvalues of the Floquet matrix [33]. In our analysis, we utilize small amplitude perturbation techniques and Floquet theory and obtain some valid theoretical results for successful management of pests. Moreover, we also establish the threshold time limit for the impulsive release of control agents. Our approach for using the Floquet theory is novel. Additionally, we examined the dynamics of the system for biological and chemical pesticides used as a sole control measure.
The paper is organized as follows. In Section 2, we derive the model by using impulsive differential equations for capturing the IPM system dynamics, taking plant, pest, virus, and chemical pesticide as model variables. The mathematical analysis of the model is then discussed in Section 3, which contains three subsections. In Section 3.1, we determine susceptible pest-eradication periodic solutions and check the feasibility–boundedness of the system variables discussed in Section 3.2. The local and global stability conditions around the susceptible pest-eradication periodic solutions are explored in Section 3.3. In Section 4, we exhibit our mathematical results through numerical simulations. Finally, in Section 5, we provide a discussion on the three types of control strategies: spraying chemical pesticide only, impulsively incorporating of infected pest only, and integrated control with a fixed and a variable impulse period, to make the final conclusion.

2. Derivation of the Impulsive Control Model

The following assumptions are taken to formulate the desired model: the crop plant and susceptible pest populations are denoted by x and y, respectively, and we denote z as the infected pest population.
Due to the finite size of a crop field, which, however, may be large, we assume logistic growth for the biomass of the crop, with net growth rate r and carrying capacity k. Crops become affected by pests, thereby causing considerable crop reduction.
Let α be the contact rate between crop and susceptible pest; let v ( t ) be the biopesticide (virus); and s be the concentration of chemical pesticide. A virus infects the susceptible pest at a rate, λ . The chemical pesticide kills the susceptible and infected pests at the rates m 1 and m 2 , respectively. Parameters c 1 and c 2 are the conversion factors of susceptible and infected pests, respectively, due to consumption of crop; d and d + δ are the mortality rates of susceptible and infected pest, respectively; θ is the virus replication rate; and γ is the lysis rate of the virus. Finally, we introduce a periodic application of biopesticide and chemical pesticide with different time intervals.
Based on the above assumptions, the desired impulsive system for integrated pest management is given as
d x d t = r x 1 x k α x y ϕ α x z , t ( n τ 1 , n τ 2 ) , d y d t = c 1 α x y λ y v d y m 1 s y , t ( n τ 1 , n τ 2 ) , d z d t = c 2 ϕ α x z + λ y v ( d + δ ) z m 2 s z , t ( n τ 1 , n τ 2 ) , d v d t = θ ( d + δ ) z γ v , t ( n τ 1 , n τ 2 ) , d s d t = μ s , t ( n τ 1 , n τ 2 ) , v ( t + ) = v ( t ) + v i , t = n τ 1 , s ( t + ) = s ( t ) + s i , t = n τ 2 ,
where v i and s i are the strength of biopesticide and chemical pesticide application in the system at t = n τ 1 and t = n τ 2 , respectively; n = 0 , 1 , 2 , 3 , , where τ 1 and τ 2 are the time periods. Here, v ( t ) and s ( t ) are the strength of biopesticide and chemical pesticide before the periodic input, and v ( t + ) and s ( t + ) are the strength of biopesticide and chemical pesticide after the periodic input.
In the impulsive model (1), we assumed the concentration of chemical pesticide as a model population, which is realistic and a novel idea. We use two different impulse intervals for two control agents (biopesticide and chemical pesticide) that will be analyzed both analytically and numerically. We proceed by analyzing the dynamics of model (1) by discussing the existence of equilibria with their stability.

3. Dynamics of the Impulsive Model

In this section, we analyze the boundedness of the solutions of system (1), we find out its pest-eradication steady state, and we analyze the local and global stability. Finally, we discuss the permanence of the impulsive system.

3.1. Boundedness of the Model Variables

Let V ( t ) = x ( t ) + y ( t ) + z ( t ) + v ( t ) + s ( t ) and
d V d t + m V r x + m x r x 2 k ( 1 c 1 ) α x y ( 1 c 2 ) ϕ α x z { ( d m } y { ( 1 θ ) ( d + δ ) m } z ( γ m ) v ( μ m ) s .
Now, let us define m = min { d , γ , μ , ( d + δ ) ( 1 θ ) } . As 0 < c 1 , c 2 , θ < 1 , then 1 c 1 > 0 , 1 c 2 > 0 , and 1 θ > 0 , so we can write that
d V ( t ) d t + m V ( t ) M 0 ,
where k ( m + r ) 2 4 r = M 0 . At t = n τ 1 , we have
V ( n τ 1 + ) V ( n τ 1 ) + v i .
By the comparison theorem, for t 0 we have
V ( t ) V ( 0 ) e m t + M 0 ( 1 e m t ) m + v i e m ( t τ 1 ) 1 e m τ 1 + v i e m τ 1 e m τ 1 1 M 0 m + v i e m τ 1 e m τ 1 1 as t .
When t = n τ 2 ,
V ( n τ 2 + ) V ( n τ 2 ) + v i
and from the comparison theorem it follows that
V ( t ) V ( 0 ) e m t + M 0 ( 1 e m t ) m + v i e m ( t τ 2 ) 1 e m τ 2 + v i e m τ 2 e m τ 2 1 M 0 m + v i e m τ 2 e m τ 2 1 as t .
Thus, V ( t ) is uniformly bounded and there exists a positive constant M > 0 such that x ( t ) M , y ( t ) M , z ( t ) M , v ( t ) M , and s ( t ) M for all t.
From the above discussion, we have the following theorem.
Theorem 1.
For the impulsive system (1), there exists a positive constant M such that x ( t ) M , y ( t ) M , z ( t ) M , v ( t ) M , and s ( t ) M for all t.
For non-negative solutions, the following lemma follows from [13].
Lemma 1.
Let X ( t ) be a solution of the impulsive system (1) with X ( 0 + ) 0 . Then X ( t ) 0 for all t > 0 .

3.2. Existence of the Pest-Free Periodic Orbit

Since both pests are assumed to be harmful for crops, we discuss stability at infected and susceptible pest-eradication solutions of the system when y = 0 and z = 0 , t ( n τ 1 , n τ 2 ) , and the linear forms of the fourth and fifth equation of (1) are
d v d t = γ v and d s d t = μ s ,
respectively. For an impulse control, we must have
v ( t + ) = v ( t ) + v i , for t = τ 1 , s ( t + ) = s ( t ) + s i , for t = τ 2 .
From (7) and (8) it is clear that v and s are independent of each other. Thus, the solution of Equation (7) can be given as follows:
v ( t ) = { v ( 0 + ) v ( 0 + ) } e γ t + v ( t ) , f o r t ( τ 1 , ( n + 1 ) τ 1 ] , s ( t ) = { s ( 0 + ) s ( 0 + ) } e μ t + s ( t ) , f o r ( τ 2 , ( n + 1 ) τ 2 ] ,
where v ( t ) and s ( t ) , the positive periodic solution of (7), are given by
v ( t ) = v i e γ ( t n τ 1 ) 1 e γ τ 1 , s ( t ) = s i e μ ( t n τ 2 ) 1 e μ τ 2 ,
with initial values
v ( 0 + ) = v i 1 e γ τ 1 , s ( 0 + ) = s i 1 e μ τ 2 .
If y ( t ) = 0 and z ( t ) = 0 , then the first equation of (1) is
d x d t = r x 1 x k ,
which is a logistic equation, and its solution is
x ( t ) = k x ( 0 ) x ( 0 ) + ( k x ( 0 ) ) e r t for t ( n τ 1 , n τ 2 ) .
Clearly, (13) has two equilibria, such as x = 0 and x = k . Therefore, (1) has two pest-eradication solutions, ( 0 , 0 , 0 , v , s ) and ( k , 0 , 0 , v , s ) . Obviously, at x = 0 the system (13) is impossible from the perspective of ecology. For this reason, in the following subsection we study the stability for the system (1) at E = ( k , 0 , 0 , v , s ) .

3.3. Stability of the Pest-Free Periodic Solution

We establish the following theorem for the stability of the pest-free periodic orbit.
Theorem 2.
System (1) is both locally and globally stable around the pest-free periodic solution E = ( k , 0 , 0 , v , s ) for the following:
(i)
Application of biopesticide and chemical pesticide with same time interval t = n τ , provided that
c 1 α d λ v i e γ ( t n τ ) 1 e γ τ m 1 s i e μ ( t n τ ) 1 e μ τ < 0 , c 2 ϕ k α ( d + δ ) m 2 s i e μ ( t n τ ) 1 e μ τ < 0 ;
(ii)
Application of biopesticide with time interval t = n τ 1 and chemical pesticide with time interval t = n τ 2 , i.e., for different time intervals, where τ 1 τ 2 , provided that
c 1 α d λ v i e γ ( t n τ 1 ) 1 e γ τ 1 < 0 , c 2 ϕ k α ( d + δ ) < 0 , c 1 α d m 1 s i e μ ( t n τ 2 ) 1 e μ τ 2 < 0 .
Proof. 
We need to prove the stability of the system in two cases:
(i) Application of chemical pesticide and biopesticide with same time interval;
(ii) Application of biopesticide and chemical pesticide with different time intervals.
(i) In this case, let t = n t 1 = n t 2 = n τ . We discuss the stability of the system through the small amplitude perturbation method at the periodic solution ( k , 0 , 0 , v , s ) . Let
x ( t ) = k + ϵ 1 ( t ) , y ( t ) = ϵ 2 ( t ) , z ( t ) = ϵ 3 ( t ) , v ( t ) = v ( t ) + ϵ 4 ( t ) , s ( t ) = s ( t ) + ϵ 5 ( t ) .
Here, ϵ 1 , ϵ 2 , ϵ 3 , ϵ 4 , and ϵ 5 denote small amplitude perturbations. Thus, the corresponding system of (1) at ( k , 0 , 0 , v , s ) is given by
d ϵ 1 d t = r { k + ϵ 1 ( t ) } 1 k + ϵ 1 ( t ) k α { k + ϵ 1 ( t ) } ϵ 2 ( t ) ϕ α { k + ϵ 1 ( t ) } ϵ 3 ( t ) , t n τ , d ϵ 2 d t = c 1 α { k + ϵ 1 ( t ) } ϵ 2 ( t ) λ ϵ 2 ( t ) { v ( t ) + ϵ 4 ( t ) } d ϵ 2 ( t ) m 1 { s ( t ) + ϵ 5 ( t ) } ϵ 2 ( t ) , t n τ , d ϵ 3 d t = c 2 ϕ α { k + ϵ 1 ( t ) } ϵ 3 ( t ) + λ ϵ 2 ( t ) { v ( t ) + ϵ 4 ( t ) } ( d + δ ) ϵ 3 ( t ) m 2 { s ( t ) + ϵ 5 ( t ) } ϵ 3 ( t ) , t n τ , d ϵ 4 d t = θ ( d + δ ) ϵ 3 ( t ) γ { v ( t ) + ϵ 4 ( t ) } , t n τ , d ϵ 5 d t = μ { s ( t ) + ϵ 5 ( t ) } , t n τ , { v ( t ) + ϵ 4 ( t ) } = v i , t = n τ , { s ( t ) + ϵ 5 ( t ) } = s i , t = n τ .
Now, the linear system corresponding to the system (17) is given as
d ϵ 1 d t = r ϵ 1 ( t ) α k ϵ 2 ( t ) ϕ α k ϵ 3 ( t ) , t n τ , d ϵ 2 d t = c 1 α k ϵ 2 ( t ) λ ϵ 2 ( t ) v ( t ) d ϵ 2 ( t ) m 1 s ( t ) ϵ 2 ( t ) , t n τ , d ϵ 3 d t = c 2 ϕ k α ϵ 3 ( t ) + λ ϵ 2 ( t ) v ( t ) ( d + δ ) ϵ 3 ( t ) m 2 s ( t ) ϵ 3 ( t ) , t n τ , d ϵ 4 d t = θ ( d + δ ) ϵ 3 ( t ) γ ϵ 4 ( t ) , t n τ , d ϵ 5 d t = μ ϵ 5 ( t ) , t n τ , ϵ 4 ( t ) = v i , t = n τ , ϵ 5 ( t ) = s i , t = n τ .
The fundamental matrix M ( t ) of (18) is obtained as
d M ( t ) d t = 1 m r α k ϕ α k 0 0 0 c 1 α k λ v ( t ) d m 1 s ( t ) 0 0 0 0 λ v ( t ) m 33 0 0 0 0 θ ( d + δ ) γ 0 0 0 0 0 μ
with initial condition M ( t ) = I 5 (the identity matrix) and m 33 = c 2 ϕ k α ( d + δ ) m 2 s ( t ) . Now, the fundamental solution matrix is given by
M ( t ) = 1 m exp ( r t ) M 1 ( t ) M 2 ( t ) 0 0 0 M 3 ( t ) 0 0 0 0 M 4 ( t ) M 5 0 0 0 0 M 5 ( t ) exp ( γ t ) 0 0 0 0 0 exp ( μ t ) .
Here, M 5 = exp 0 t { c 2 ϕ k α ( d + δ ) m 2 s ( t ) } d t ,
M 3 ( t ) = exp 0 τ { c 1 α k λ v ( t ) d m 1 s ( t ) } d t ,
where the other M i ( t ) s are not required for our further analysis. According to Floquet theory [33], the periodic solution E ( k , 0 , 0 , v v , s ) is asymptotically stable if the absolute values of the eigenvalues of M ( τ ) are less than one.
The eigenvalues of M ( t ) are
λ 1 = exp { r τ } , λ 2 = exp 0 τ { c 1 α k λ v ( t ) d m 1 s ( t ) } d t , λ 3 = exp 0 τ { c 2 ϕ k α ( d + δ ) m 2 s ( t ) } d t , λ 4 = exp { γ τ } , λ 5 = exp { μ τ } .
Clearly, 0 < λ 1 < 1 , 0 < λ 4 < 1 and 0 < λ 5 < 1 . Thus, when both pesticides are applied with the same time interval, then the system is locally stable around the periodic solution E = ( k , 0 , 0 , v v , s ) if 0 < λ 2 < 1 and 0 < λ 3 < 1 . From this, we obtain
exp 0 τ { c 1 α k λ v ( t ) d m 1 s ( t ) } d t < 1 , exp 0 τ { c 2 ϕ k α ( d + δ ) m 2 s ( t ) } d t < 1 .
From Equation (19), we can choose δ 1 > 0 such that
η 1 = exp n τ ( n + 1 ) τ { c 1 α λ ( v ( t ) δ 1 ) d m 1 ( s ( t ) δ 1 ) } d t < 1 , η 2 = exp n τ ( n + 1 ) τ { c 2 ϕ k α ( d + δ ) m 2 ( s ( t δ 1 ) ) } d t < 1 .
Since all state variables are positive,
d v d t γ v ; d s d t μ s .
Thus, according to the comparison theorem and Equations (10) and (11), for δ 1 > 0 there exists t 0 > 0 such that v ( t ) v δ 1 , s ( t ) s δ 1 for all t > t 0 .
From the second equation of system (1), it can be written that
y ˙ ( t ) y ( t ) { c 1 α λ ( v ( t ) δ 1 ) d m 1 ( s ( t ) δ 1 ) } , t n τ , y ( t + ) = y ( t ) , t = n τ .
Integrating (20) into [ n τ , ( n + 1 ) τ ] , it can be shown that
y { ( n + 1 ) τ } y ( n τ ) exp n τ ( n + 1 ) τ { c 1 α λ ( v ( t ) δ 1 ) d m 1 ( s ( t ) δ 1 ) } d t = y ( n τ ) η 1 .
Similarly,
y { n τ } y { ( n 1 ) τ } η 1 .
Hence, from (21) and (22),
y { ( n + 1 ) τ } y { ( n 1 ) τ } η 1 2 .
Proceeding in this way, we obtain
y { ( n + 1 ) τ } y ( τ ) η 1 n .
Since η 1 < 1 , one has η 1 n 0 whenever n . Hence, y { ( n + 1 ) τ } 0 as n . Now we take n τ < t ( n + 1 ) τ . Then, clearly, 0 < y ( t ) y ( n τ ) exp ( n τ ) . Thus, y ( t ) 0 as t . For
η 2 = exp n τ ( n + 1 ) τ { c 2 ϕ k α ( d + δ ) m 2 ( s ( t δ 1 ) ) } d t < 1 ,
we can similarly prove that z ( t ) 0 as t .
We now prove that v ( t ) v ( t ) as t . Since z ( t ) 0 as t , then for some 0 < δ 2 < γ θ ( d + δ ) there exists t 1 > 0 such that 0 < z ( t ) < δ 2 for all t > t 1 . Thus, for t > t 1 and from the fourth equation of system (1), we can write that
θ ( d + δ ) δ 2 γ v ( t ) v ˙ ( t ) θ ( d + δ ) δ 2 γ v ( t ) .
Let v 1 ( t ) and v 2 ( t ) be the solutions of
v ˙ 1 ( t ) = θ ( d + δ ) δ 2 γ v 1 ( t ) , t n τ , v 1 ( t + ) = v 1 ( t ) + v i , t = n τ ,
and
v ˙ 2 ( t ) = θ ( d + δ ) δ 2 γ v 2 ( t ) , t n τ , v 2 ( t + ) = v 2 ( t ) + v i , t = n τ ,
respectively. Then, the solution will be
v 1 ( t ) = v i e γ ( t n τ ) 1 e γ τ + θ ( d + δ ) δ 2 , v 2 ( t ) = v i e γ ( t n τ ) 1 e γ τ θ ( d + δ ) δ 2 .
From (25), it is clear that when δ 2 0 we have v 1 ( t ) v ( t ) and v 2 ( t ) v ( t ) . Hence, it follows from (24) that v ( t ) v ( t ) as t .
Similarly, we can choose 0 < δ 3 < μ and, in the same way, we can prove that s ( t ) s ( t ) as t .
Finally, we shall prove that x ( t ) k as t . We already proved that y ( t ) , z ( t ) 0 as t . Thus, for δ 3 > 0 , there exists t 3 > 0 such that y ( t ) , z ( t ) < δ 3 for all t > t 3 . Hence, from the first equation of system (1), we can write that
r x ( t ) r x 2 ( t ) k > x ˙ ( t ) { r δ 3 α ( 1 + ϕ ) } x ( t ) r x 2 ( t ) k ,
which implies that
k x 0 x 0 + ( k x 0 ) e r t x ( t ) k { r δ 3 α ( 1 + ϕ ) } x 0 r x 0 + [ k { r δ 3 α ( 1 + ϕ ) } r x 0 ] e { r δ 3 α ( 1 + ϕ ) } t .
Hence, for δ 3 0 , x ( t ) k as t . Thus, for application of biopesticide and chemical pesticide together with the same time interval t = n τ , we can say that system (1) is locally as well as globally stable if
c 1 α d λ v i e γ ( t n τ ) 1 e γ τ m 1 s i e μ ( t n τ ) 1 e μ τ < 0 , c 2 ϕ k α ( d + δ ) m 2 s i e μ ( t n τ ) 1 e μ τ < 0 .
Two subcases arise here, namely,
Subcase I. Application of biopesticide with time interval t = n τ 1 .
In this case, s i = 0 . Hence, system (1) is locally as well as globally stable around the periodic solution if
c 1 α k d λ v i e γ ( t n τ 1 ) 1 e γ τ 1 < 0 , c 2 ϕ k α ( d + δ ) < 0 .
Subcase II. Application of chemical pesticide with time interval t = n τ 2 .
In this case, v i = 0 , and hence system (1) is locally as well as globally stable around the periodic solution if
c 1 α k d m 1 s i e μ ( t n τ 2 ) 1 e μ τ 2 < 0 , c 2 ϕ k α ( d + δ ) m 2 s i e μ ( t n τ 2 ) 1 e μ τ 2 < 0 .
The proof is complete.    □

4. Numerical Simulations

Now we solve the impulsive system numerically and we graphically display the results found. We varied the crucial parameters within their feasible ranges to observe their impact on the impulsive model’s solution trajectories and equilibria. Precisely, we solved the impulsive system and plotted the results in figures using the ode45 MATLAB solver.
In Figure 1, the impulsive time interval for microbial biological pesticide release is 5 days, and releasing of biopesticide was considered in different rates: v i = 0 (i.e., without pesticides), v i = 6 , and v i = 12 . It is revealed that susceptible pest population decreases with an increase in the release rate of biopesticides.
In Figure 2, by taking different impulsive intervals, biopesticide is applied to the system. A better result is obtained for lower intervals (2 days) but, with a higher release of biopesticide, pests are present in the system. From Figure 1 and Figure 2, we can conclude that pest control using only biopesticides is very costly and a time-consuming process.
Recall that in our model we take τ 1 as the time period for biopesticide spraying (generally a virus particle) and τ 2 as the time period for chemical pesticide sprays. In Figure 3 we see the effect for the same time intervals, τ 1 = τ 2 = 5 days, whereas in Figure 4 we see the effect for different time intervals, τ 1 = 3 days and τ 2 = 2 days.
If both microbial biopesticides and chemical pesticides are released simultaneously, with an equal time interval of 5 days, then the extinction of both infected and susceptible pest populations is possible (see Figure 3).
Figure 4 illustrates the dynamics of the double impulse with different impulse intervals. Double impulses occur at the time which is the common multiple of the two intervals. For example, if we take τ 1 = 2 and τ 2 = 3 , then simultaneous impulses will occur at the times t = 6 , t = 12 , t = 18 , and so on. Figure 4 is the most important figure characterizing the impact of two different but simultaneous impulses on the total pest population with different time intervals. It is observed that for v i = 12 , s i = 0.15 , τ 1 = 3 , and τ 2 = 2 , the total pest population becomes extinct. In Figure 4d, the impact of double impulses occurs at t = 48 , t = 54 , t = 60 days, etc., which are common multiples of the two intervals τ 1 = 2 and τ 2 = 3 .
It is also numerically checked that when the rate of the impulse control is high, a comparatively lower interval can be taken for cost-effectiveness of the process.
Thus, the advantage of the impulsive control is that we can determine the proper rate and a suitable interval of giving controlling agents in the system.

5. Discussion and Conclusions

In the present research, we studied impulsive periodic applications of integrated pesticides, that is, simultaneous use of biopesticide and chemical pesticide in a pest management system. We proposed a two-impulse mathematical model using an impulsive differential equation to observe the impact of periodic application of the combined pesticides in impulsive modes.
In the previous models available in the literature, chemical and biological pesticides were used in the model in a continuous way. In contrast, here, we used them in an impulsive periodic way. Consequently, a two-impulse mathematical model was established. Moreover, in the proposed model, we took chemical pesticides concentration as the model population, which is a novel approach.
Stability theory (Floquet theory) and numerical calculations were used to examine the system dynamical behavior. We determined the conditions under which the impulsive system will be stable both locally and globally. For example, the local stability of a pest-free periodic orbit was established. The dynamics varied with the rate of both biopesticide recruitment and the chemical pesticide concentration.
Chemical pesticides minimize the oscillations in the system and make the system stable in a shorter time. Numerical and analytical analysis reveals that increasing frequency of pesticide application will require less administration of biopesticide and chemical pesticide, which is economically beneficial and environmentally safe.
Our research is directed to optimize and find the right combination of pesticides with maximum benefit to the crop plant. The numerical simulation also shows that control over the spraying of chemical pesticides is needed to control pests and minimize the cost of cultivation. On the other hand, chemical pesticides may have negative environmental implications due to their lingering effects; nonetheless, the best control approach provides the least amount of collateral damage to the environment.
In a nutshell, the promising feature of the system is the combined use of the pesticides in impulsive control methods that reduce the cost and negative effects on the environment. Using a combination of pesticides to deliver the pesticide can save the cost and reduce the side effects of chemical pesticides. Our obtained results will give a new perspective to farmers who implement this in a real-world setting.
In the future, one can extend this work to an optimal impulsive system for cost-effectiveness of the control process.

Author Contributions

Conceptualization, F.A.B. and J.C.; methodology, F.A.B., J.C. and D.F.M.T.; software, F.A.B. and J.C.; validation, F.A.B., J.C. and D.F.M.T.; formal analysis, F.A.B., J.C. and D.F.M.T.; investigation, F.A.B., J.C. and D.F.M.T.; writing—original draft preparation, F.A.B., J.C. and D.F.M.T.; writing—review and editing, F.A.B., J.C. and D.F.M.T.; visualization, F.A.B. and J.C. All authors have read and agreed to the published version of the manuscript.

Funding

Torres was funded by The Portuguese Foundation for Science and Technology (FCT) grant number UIDB/04106/2020.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable to this article.

Acknowledgments

The authors are very grateful to four anonymous reviewers for careful reading of the submitted manuscript and also for providing several important comments and suggestions.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Al Basir, F.; Banerjee, A.; Ray, S. Role of farming awareness in crop pest management—A mathematical model. J. Theor. Biol. 2019, 461, 59–67. [Google Scholar] [CrossRef]
  2. Ndolo, D.; Njuguna, E.; Adetunji, C.O.; Harbor, C.; Rowe, A.; Den Breeyen, A.; Sangeetha, J.; Singh, G.; Szewczyk, B.; Anjorin, T.S.; et al. Research and development of biopesticides: Challenges and prospects. Outlooks Pest Manag. 2019, 6, 267–276. [Google Scholar] [CrossRef]
  3. Chattopadhyay, P.; Banerjee, G.; Mukherjee, S. Recent trends of modern bacterial insecticides for pest control practice in integrated crop management system. 3 Biotech 2017, 7, 60. [Google Scholar] [CrossRef]
  4. Perring, T.M.; Gruenhagen, N.M.; Farrar, C.A. Management of plant viral diseases through chemical control of insect vectors. Annu. Rev. Entomol. 1999, 44, 457–481. [Google Scholar] [CrossRef] [PubMed]
  5. Ghosh, S.; Bhattacharya, D.K. Optimization in microbial pest control: An integrated approach. Appl. Math. Model. 2010, 34, 1382–1395. [Google Scholar] [CrossRef]
  6. Bhattacharyya, S.; Bhattacharya, D.K. Pest control through viral disease: Mathematical modeling and analysis. J. Theor. Biol. 2006, 238, 177–197. [Google Scholar] [CrossRef] [PubMed]
  7. Fest, C.; Schmidt, K. The Chemistry of Organophosphorus Pesticides; Springer: New York, NY, USA, 1973. [Google Scholar] [CrossRef]
  8. Lenteren, J.C.V. Integrated pest management in protected crops. In Integrated Pest Management; Dent, D., Ed.; Chapman and Hall: London, UK, 1995; pp. 311–320. Available online: https://cir.nii.ac.jp/crid/1573950400276916864 (accessed on 20 January 2023).
  9. Flint, M.L. Integrated Pest Management for Walnuts; University of California Statewide Integrated Pest Management Project; Division of Agriculture and Natural Resources, University of California: Davis, CA, USA, 1987; Volume 3270, p. 3641. [Google Scholar]
  10. Stern, V.M. The bioeconomics of pest control. Iowa State J. Res. 1975, 49, 467–472. Available online: https://core.ac.uk/download/pdf/224978255.pdf (accessed on 20 January 2023).
  11. Liu, B.; Teng, Z.; Chen, L. Analysis of a predator-prey model with Holling II functional response concerning impulsive control strategy. J. Comput. Appl. Math. 2006, 193, 347–362. [Google Scholar] [CrossRef]
  12. Páez Chávez, J.; Jungmann, D.; Siegmund, S. Modeling and analysis of integrated pest control strategies via impulsive differential equations. Int. J. Differ. Equ. 2017, 2017, 1820607. [Google Scholar] [CrossRef]
  13. Wang, L.; Chen, L.; Nieto, J.J. The dynamics of an epidemic model for pest control with impulsive effect. Nonlinear Anal. Real World Appl. 2010, 11, 1374–1386. [Google Scholar] [CrossRef]
  14. Al Basir, F.; Noor, M.H. A Model for Pest Control using Integrated Approach: Impact of latent and gestation Delays. Nonlinear Dyn. 2022, 108, 1805–1820. [Google Scholar] [CrossRef]
  15. Abraha, T.; Al Basir, F.; Obsu, L.L.; Torres, D.F.M. Pest control using farming awareness: Impact of time delays and optimal use of biopesticides. Chaos Solitons Fractals 2021, 146, 110869. [Google Scholar] [CrossRef]
  16. Chowdhury, J.; Basir, F.A.; Takeuchi, Y.; Ghosh, M.; Roy, P.K. A mathematical model for pest management in Jatropha curcas with integrated pesticides—An optimal control approach. Ecol. Complex. 2019, 37, 24–31. [Google Scholar] [CrossRef]
  17. Tang, S.; Chen, L. Modelling and analysis of integrated pest management strategy. Discret. Contin. Dyn. Syst. Ser. B 2004, 4, 759–768. [Google Scholar] [CrossRef]
  18. Paez Chavez, J.; Jungmann, D.; Siegmund, S. A comparative study of integrated pest management strategies based on impulsive control. J. Biol. Dyn. 2018, 12, 318–341. [Google Scholar] [CrossRef] [PubMed]
  19. Kalra, P.; Kaur, M. Stability analysis of an eco-epidemiological SIN model with impulsive control strategy for integrated pest management considering stage-structure in predator. Int. J. Math. Model. Numer. Optim. 2022, 12, 43–68. [Google Scholar] [CrossRef]
  20. Liu, J.; Hu, J.; Yuen, P. Extinction and permanence of the predator-prey system with general functional response and impulsive control. Appl. Math. Model. 2020, 88, 55–67. [Google Scholar] [CrossRef]
  21. Tian, B.; Li, J.; Wu, X.; Zhang, Y. Dynamic behaviour of a predator-prey system with impulsive control strategy. Int. J. Dyn. Syst. Differ. Equ. 2022, 12, 493–509. [Google Scholar] [CrossRef]
  22. Kumari, V.; Chauhan, S.; Dhar, J. Controlling pest by integrated pest management: A dynamical approach. Int. J. Math. Eng. Manag. Sci. 2020, 5, 769–786. [Google Scholar] [CrossRef]
  23. Pang, Y.; Wang, S.; Liu, S. Dynamics analysis of stage-structured wild and sterile mosquito interaction impulsive model. J. Biol. Dyn. 2022, 16, 464–479. [Google Scholar] [CrossRef]
  24. Chowdhury, J.; Al Basir, F.; Cao, X.; Roy, P.K. Integrated pest management for Jatropha Carcus plant: An impulsive control approach. Math. Methods Appl. Sci. 2021, in press. [CrossRef]
  25. Tian, Y.; Tang, S.; Cheke, R.A. Dynamic complexity of a predator-prey model for IPM with nonlinear impulsive control incorporating a regulatory factor for predator releases. Math. Model. Anal. 2019, 24, 134–154. [Google Scholar] [CrossRef]
  26. Jose, S.A.; Ramachandran, R.; Cao, J.; Alzabut, J.; Niezabitowski, M.; Balas, V.E. Stability analysis and comparative study on different eco-epidemiological models: Stage structure for prey and predator concerning impulsive control. Optim. Control. Appl. Methods 2022, 43, 842–866. [Google Scholar] [CrossRef]
  27. Li, J.; Huang, Q.D.; Liu, B. A pest control model with birth pulse and residual and delay effects of pesticides. Adv. Differ. Equ. 2019, 2019, 117. [Google Scholar] [CrossRef]
  28. Liu, B.; Kang, B.-L.; Tao, F.-M.; Hu, G. Modelling the effects of pest control with development of pesticide resistance. Acta Math. Appl. Sin. Engl. Ser. 2021, 37, 109–125. [Google Scholar] [CrossRef]
  29. Alzabut, J. An Integrated Eco-Epidemiological Plant Pest Natural Enemy Differential Equation Model with Various Impulsive Strategies. 2022. Available online: http://earsiv.ostimteknik.edu.tr:8081/xmlui/handle/123456789/207 (accessed on 20 January 2023).
  30. Zhang, H.; Xu, W.; Chen, L. A impulsive infective transmission SI model for pest control. Math. Methods Appl. Sci. 2007, 30, 1169–1184. [Google Scholar] [CrossRef]
  31. Al Basir, F.; Chowdhury, J.; Das, S.; Ray, S. Combined impact of predatory insects and bio-pesticide over pest population: Impulsive model-based study. Energy Ecol. Environ. 2022, 7, 173–185. [Google Scholar] [CrossRef]
  32. Shi, R.; Jiang, X.; Chen, L. A predator-prey model with disease in the prey and two impulses for integrated pest management. Appl. Math. Model. 2009, 33, 2248–2256. [Google Scholar] [CrossRef]
  33. Klausmeier, C.A. Floquet theory: A useful tool for understanding nonequilibrium dynamics. Theor. Ecol. 2008, 1, 153–161. [Google Scholar] [CrossRef]
Figure 1. Impact of biopesticide application in impulsive mode on system (1). Evolution of (a) crop; (b) uninfected pest; (c) infected pest; (d) biopesticide. The set of parameters are r = 0.1 , k = 1 , α = 0.2 , β = 0.003 , m 1 = 0.8 , m 2 = 0.6 , c 1 = 0.5 , c 2 = 0.8 , γ = 0.15 , δ = 0.2 , d = 0.05 , κ = 100 , s = 0.3 , and λ = 0.35 . Here, the time interval is τ 1 = 5 days and the rates of biopesticide release are v i = 0 (black line), v i = 6 (red line), and v i = 12 (blue line).
Figure 1. Impact of biopesticide application in impulsive mode on system (1). Evolution of (a) crop; (b) uninfected pest; (c) infected pest; (d) biopesticide. The set of parameters are r = 0.1 , k = 1 , α = 0.2 , β = 0.003 , m 1 = 0.8 , m 2 = 0.6 , c 1 = 0.5 , c 2 = 0.8 , γ = 0.15 , δ = 0.2 , d = 0.05 , κ = 100 , s = 0.3 , and λ = 0.35 . Here, the time interval is τ 1 = 5 days and the rates of biopesticide release are v i = 0 (black line), v i = 6 (red line), and v i = 12 (blue line).
Axioms 12 00391 g001
Figure 2. Impact of biopesticide on system (1) for different impulse intervals and rates. Evolution of (a) crop; (b) uninfected pest; (c) infected pest; (d) biopesticide. Red line indicates v i = 6 and τ 1 = 5 , green line indicates v i = 12 and τ 1 = 5 , and blue line indicates v i = 12 and τ 1 = 2 .
Figure 2. Impact of biopesticide on system (1) for different impulse intervals and rates. Evolution of (a) crop; (b) uninfected pest; (c) infected pest; (d) biopesticide. Red line indicates v i = 6 and τ 1 = 5 , green line indicates v i = 12 and τ 1 = 5 , and blue line indicates v i = 12 and τ 1 = 2 .
Axioms 12 00391 g002
Figure 3. Impact of both biopesticide and chemical pesticide on system (1) with the same impulse interval, τ 1 = τ 2 = 5 . Evolution of (a) crop; (b) uninfected pest; (c) infected pest; (d) biopesticide; (e) chemical pesticide. The rates of impulses are s i = 0.15 and v i = 6 for black dotted color; s i = 0.1 and v i = 6 for red dashed line; and s i = 0.05 and v i = 6 for blue solid line.
Figure 3. Impact of both biopesticide and chemical pesticide on system (1) with the same impulse interval, τ 1 = τ 2 = 5 . Evolution of (a) crop; (b) uninfected pest; (c) infected pest; (d) biopesticide; (e) chemical pesticide. The rates of impulses are s i = 0.15 and v i = 6 for black dotted color; s i = 0.1 and v i = 6 for red dashed line; and s i = 0.05 and v i = 6 for blue solid line.
Axioms 12 00391 g003
Figure 4. Impact of both biopesticide and chemical pesticide on system (1) with same impulse interval, τ 1 = 3 , τ 2 = 2 , and where the rates of impulses are s i = 0.15 and v i = 6 . Evolution of (a) crop; (b) uninfected pest; (c) infected pest; (d) pesticides (biopesticide in blue, chemical pesticide in red).
Figure 4. Impact of both biopesticide and chemical pesticide on system (1) with same impulse interval, τ 1 = 3 , τ 2 = 2 , and where the rates of impulses are s i = 0.15 and v i = 6 . Evolution of (a) crop; (b) uninfected pest; (c) infected pest; (d) pesticides (biopesticide in blue, chemical pesticide in red).
Axioms 12 00391 g004
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Al Basir, F.; Chowdhury, J.; Torres, D.F.M. Dynamics of a Double-Impulsive Control Model of Integrated Pest Management Using Perturbation Methods and Floquet Theory. Axioms 2023, 12, 391. https://doi.org/10.3390/axioms12040391

AMA Style

Al Basir F, Chowdhury J, Torres DFM. Dynamics of a Double-Impulsive Control Model of Integrated Pest Management Using Perturbation Methods and Floquet Theory. Axioms. 2023; 12(4):391. https://doi.org/10.3390/axioms12040391

Chicago/Turabian Style

Al Basir, Fahad, Jahangir Chowdhury, and Delfim F. M. Torres. 2023. "Dynamics of a Double-Impulsive Control Model of Integrated Pest Management Using Perturbation Methods and Floquet Theory" Axioms 12, no. 4: 391. https://doi.org/10.3390/axioms12040391

APA Style

Al Basir, F., Chowdhury, J., & Torres, D. F. M. (2023). Dynamics of a Double-Impulsive Control Model of Integrated Pest Management Using Perturbation Methods and Floquet Theory. Axioms, 12(4), 391. https://doi.org/10.3390/axioms12040391

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop