A mathematical model for fall armyworm management on maize biomass

*Correspondence: daudis@nm-aist.ac.tz 1School of Computational and Communication Science Engineering, The Nelson Mandela African Institute of Science and Technology (NM-AIST), P.O. Box 447, Arusha, Tanzania 2Department of Mathematics, Humanities and Social Science (MHSS), National Institute of Transport (NIT), Tanzania, P.O. Box 705, Dar-es-Salaam, Tanzania Full list of author information is available at the end of the article Abstract


Introduction
Fall armyworm (FAW), Spodoptera frugiperda, a highly destructive and fast spreading agricultural pest native to North and South America, poses a real threat to global food security. Spodoptera frugiperda remains an important pest of members of family Poaceace including major food crops such as corn, sorghum, rice wheat, maize, and diverse pasture [1]. According to FAO [2], food security is defined as a "situation that exists when all people, at all times, have physical, social, and economic access to sufficient, safe, and nutritious food that meets their dietary needs and food preferences for an active and healthy life" [2]. In recent years, the FAW has spread globally and emerged in countries where it had rarely or never before been presented, posing a real threat to global food security [3,4]. Prior studies suggest that FAW pest is native to and widely distributed in the tropical and subtropical regions of America [1], and its invasion into Africa was reported for the first time in January 2016 [1]. Since then, it has become an epidemic pest in and beyond several African countries [1,[3][4][5]. FAW is regarded to be a major pest of maize biomass and other crops, such as rice, millet, sorghum, and cotton [5]. Due to the importance of maize biomass in many countries, there is a need to explore the implications of FAW outbreak in a field of maize biomass to increase the harvest of this crop.
Mathematical modeling has become a tool used to explore many real-world phenomena. Ordinary differential equations (ODEs) and partial differential equations (PDEs) with and without memory effects are some of the tools that have been commonly used to formulate equation(s) that mirror the real-world problem(s) [6][7][8][9][10][11][12]. In recent years a number of mathematical models were developed to explore plant-pest interaction [11][12][13][14][15][16][17][18][19]. Tang et al. [17] proposed impulsive differential equation models or hybrid dynamical system to model the introduction of a periodic IPM strategy, which includes periodic spraying of pesticide and release of natural enemies at critical time [14][15][16][17]. On the other hand, Tang et al. [17] developed an impulsive pest-natural enemy model, in which pulsing actions such as spraying pesticide and releasing natural enemies were considered with the assumption that the pesticide kills a pest instantly, whereas Chowdhury [18] formulated and extensively investigated continuous and discrete predator-prey models concerning IPM strategy. Discrete host parasitoid models were also proposed in circumstances where the timing of pesticide application leads to the death of parasitoid, and four different cases involving the timing of applications were investigated [11,18].
The aforementioned studies and several other cited therein have certainly produced many useful results and improved the existing knowledge on plant-pest interaction. However, one of the limitations of these studies is that their models were based on integer-order ordinary differential equations. Recent studies suggest that models that use integer-order differential equations do not adequately capture memory effects and hereditary properties, which are inherent in many real-world problems [20]. As such, in recent years, fractional calculus has become an intriguing field. Several researchers have shown that models that utilize fractional calculus are more likely to replicate real-world problems compared those that use integer-order differential equations since fractional-order differential equations are able to capture memory effects [20,21].
Therefore the present work aims to utilize fractional calculus to explore the implications of FAW infestation in a field planted with initial number of maize biomass at time t = 0 and obtaining maximum harvest of the biomass at the end of the season. The proposed model incorporates all the relevant biological information. In particular, the FAW population has been subdivided into egg population, larvae population, pupae population, and the adult population, also known as moth. Although the FAW has six larval instar stages, we have considered this as a single group to reduce complexity of the model. The proposed model also incorporates the use of nonbiological control methods such as use of pesticides and commonly known traditional methods such as hand picking of caterpillars. The role of biological control on FAW dynamics was comprehensively explored in [12], and hence we did not consider this aspect. At larval stage cannibalism is known to occur in FAW dynamics, and we also incorporate this aspect. We support qualitative and quantitative analytical results obtained in this study by numerical illustration.
We organize the paper as follows. In Sect. 2, we give some necessary definitions and some known properties of fractional calculus. In Sect. 3, we propose a fractional-order model for fall armyworm and maize biomass interaction. We investigate the local and global stability of the model equilibria. To support analytical findings, we carry out numerical simulation and present their results in Sect. 4. In the last section, we present the conclusions of this paper.

Preliminaries on the Caputo fractional calculus
We begin by introducing the definition of Caputo fractional derivative and state related theorems (see [22][23][24]), which we will utilize to derive important results in this work.
where is the gamma function. 22,24]) The Riemann-Liouville fractional integral of arbitrary real order q > 0 of a function f (t) is defined by the integral and J 0 f (t) = f (t).
are continuous on [t 0 , ∞) and have the exponential order and that c Lemma 2.1 ([26]) Let x(·) be a continuous and differentiable function with x(t) ∈ R + . Then, for any time instant t ≥ t 0 , we have

Model formulation
A fractional-order model we introduce consists of two populations, maize biomass and the FAW population, where one of the populations is a stage-structured giving a total of five populations. Meanwhile, the FAW population is divided into four classes, which represent the FAW life cycle and are the egg stage E(t), Larvae L(t), pupal P(t), and the adult stage (Moth) A(t). Although the FAW typically has six larval instars, to reduce complexity of the model in a biological sensible way, all larval instars are represented by class L(t). The life cycle of FAW starts when eggs are laid in masses on maize biomass, mostly underside of these biomass. The following equation describes the FAW egg population dynamics: where b represents the egg laying rate for an adult female FAW, that is, an average number of eggs each female adult FAW will lay per day, K E represents the egg carrying capacity, that is, the availability of space to lay eggs, W is the proportion of adult FAW that are females, α E is the egg hatching rate, and u E accounts for the effects of intervention strategies a farmer will implement once they observe eggs laid on the maize biomass, μ E is the egg mortality rate. FAW larvae generally emerge simultaneously three to five days following oviposition. The following equation summarizes the population dynamics of the larvae stage: where the transition rate from the egg stage to larvae is α E . Prior studies [3,27] suggest that whenever the food is limited, the older larvae of FAW exhibit a cannibalistic behavior on the smaller larvae. Hence to account for this aspect, we assume that the death rate due to lack of food is proportional to the smaller larvae α E E and to the coefficient L/K L that represent the availability of food for each larvae. Therefore K L models the availability of food and space for the larvae population, μ L represents the natural mortality rate of the larvae, and 1/α L is the average duration of the larvae stage, which is estimated to the of 14-30 days [3,28,29]. In particular, it is estimated that this duration is shorter, around 14 days during warm summer months and longer, and around 30 days during cooler weather [28,29]. The parameter u L models the role of intervention strategies implemented by the farmer, which may be use of pesticides or handpicking of the larvae. The term θ q LM represents the interaction of larvae and maize biomass, which results in conversion of maize biomass into larvae biomass. Hence we can write θ q = e q β q , where the parameter e represents the efficiency with which caterpillar (FAW larvae) convert consumed maize biomass.
Pupation of the FAW normally occurs in the soil at a depth of 2-8 cm [29]. Here the larva constructs a loose cocoon, oval in shape and 20-30 mm in length, through tying together particles of soil with silk [3]. In areas where the soil is too hard, larvae web together leaf debris and other material to form a cocoon on the soil surface [3]. The following equation represents the dynamics of pupae stage: where μ P is the natural mortality rate, 1/α P represents the duration of the pupae stage, which is approximately 8-9 days during the summer; however, during winter it may reach 20-30 days [3,28]. It is worth noting that the pupal stage of FAW does not enter a diapause period to withstand protracted periods winter or summer seasons in the absence of the host plant biomass [3,30]. The effects of intervention strategies on reducing the population of the pupae is modeled by u P . Adult FAW are 20-25-mm long with a wingspan of approximately 30-40 mm. Female and male adult FAW have different color pattern on their forewing. Adult female FAW are responsible for laying eggs on the surface of maize biomass, and this process usually starts after a preoviposition period of 3-4 days and continues until they become 3-week old. The following equation summarizes the population dynamics of the adult FAW: where α P accounts for the proportion of FAW pupa population that successfully progresses to the adult stage, u A denotes the effects of intervention strategies, and 1/μ A is the life span, which is estimated to average about 10 days, with a range of about 7-21 days [3]. It is worth noting that the duration of the life cycle for FAW lasts for about 30 days at 28 • C and may take longer, 60-90 days, when the weather is cooler [28]. In addition, under favorable conditions, the FAW larvae have a potential to feed and breed on maize biomass year-round [28]. Plant biomass (plant seeds) planted at time t = 0 emerges in a period of 0 to 7 days. We assume that planting of maize seed per hectare at the beginning of the season is done in a day. In this regard, we let M(t) represent the population density of maize biomass per hectare Therefore the population dynamics of the maize biomass is governed by the equation where r is the growth rate of maize biomass, K M is the maximum biomass of maize plants, and β is rate at which larvae (FAW larva) attacks the biomass of the maize plants.
Our assumptions on the dynamics of fall armyworm in a maize biomass population density are demonstrated in Fig. 1, and equations are presented in the system

Positivity and boundedness of solutions to model (1)
Theorem 3.1 There exists a unique solution for the fractional-order model (1) in (0, ∞). Moreover, the solution is nonnegative for all t > 0 and remains in R 5 + .
Proof We begin by demonstrating that R 5 + = {(M, E, L, P, A) ∈ R 5 + : M ≥ 0, E ≥ 0, L ≥ 0, P ≥ 0, A ≥ 0} is positively invariant. For that, we have to demonstrate that on each hyperplane bounding the nonnegative orthant the vector field points to R 5 + . Let us consider the following cases.
Case 1. Assume that the exists t * > t 0 such that M(t * ) = 0 and M(t) < 0 for t ∈ (t * , t 1 ], where t 1 is sufficiently close to t * . If M(t * ) = 0, then From the last two equations of system (1) we can easily verify that From the discussion above we observe that each hyperplane bounding the nonnegative orthant, the vector field points to R 5 + , that is, all solutions of system (1), remains nonnegative for all t ≥ 0. (1) for t ≥ 0. Then the solution X (t) is bounded above, that is, X (t) ∈ , where denotes the feasible region given by

Theorem 3.2 Let X (t) = (E(t), L(t), P(t), A(t)) be a unique solution of the model
Proof Here we will now demonstrate that the solutions of system (1) are bounded for all t ≥ 0. For biological relevance, the least possible lower bound for each of the variables in system (1) is zero. Based on this, our discussion will be on determining the upper bound for these variables. Moreover, we can easily establish that the following conditions should Therefore it follows that lim sup t→∞ M(t) ≤ K M . Based on these bounds, we have Applying the Laplace transform leads to Grouping like terms, we get .
Applying the inverse Laplace transform leads to Applying the Laplace transform leads to Combining the like terms, we get .
Applying the inverse Laplace transform leads to By applying the Laplace transform it follows that Grouping similar terms, we have Utilizing inverse Laplace transform, we get is bounded from above. This completes the proof.

Model equilibria
By direct calculations we can observe that system (1) has four equilibrium points: (i) Trivial equilibrium (ii) First axial equilibrium point We can observe that this equilibrium point makes biological sense whenever Let Biologically, R 0 is a threshold quantity that accounts for the persistence of the FAW population, and thus when R 0 > 1, the population of FAW persists and will be an attack on maize plants leaves, and finally the population of maize plants is extinct. Hence we can precisely define R 0 as the average number of off-spring generated by an adult female FAW during its entire life span. Precisely, we can note that a proportion W of moth will each lay b eggs per day for an average duration of Caterpillars that emerge following oviposition have the probability (iv) Interior equilibrium point where

Local stability of equilibrium points
In this section, we study the local stability behavior of the four equilibrium points computed earlier by using the Jacobian matrix (a) Trivial equilibrium point. Evaluating the Jacobian matrix (4) about E 1 leads to The trivial equilibrium point is locally stable if all eigenvalues λ i (i = 1, 2, 3, 4) of J(E 1 ) satisfy the condition |arg(λ i )| > qπ 2 [31]. We can observe that one of the eigenvalues of (5) is r q > 0. The other equilibrium points are obtained from the characteristic equation The Routh-Hurwitz criteria for local asymptotic stability of the equilibrium point E 1 are As we can observe, all the coefficients of the characteristic polynomial (6) are positive whenever R 0 < 1, implying that condition H 1 holds for R 0 < 1. Since we have established that the trivial equilibrium point E 1 has another eigenvalue r q , which is always positive, we will not investigate the positivity of condition H 2 , and hence we conclude that E 1 is an unstable equilibrium point.
(b) First axial equilibrium point E 2 . Evaluating the Jacobian matrix (4) about E 2 leads to From (8) we can observe that one of the eigenvalues is -r q < 0, and the other eigenvalues are roots of the characteristic equation The Routh-Hurwitz criteria for local asymptotic stability of the equilibrium point E 2 are If conditions specified in (10) hold, then the equilibrium point E 2 is locally asymptotically stable. (c) Second axial equilibrium point E 3 . Evaluating the Jacobian matrix (4) about E 3 , we get with From (11) we can observe that -r q ( β q L 3 r q -1) is an eigenvalue, and other eigenvalues can be determined from the characteristic polynomial Ahmed et al. [31] presented some Routh-Hurwitz stability conditions for fractional-order systems. One well-known Routh-Hurwitz condition is that an equilibrium point is locally stable if all eigenvalues of the community matrix satisfy the condition |arg(λ i )| > qπ 2 . The Routh-Hurwitz criteria for the local asymptotic stability of the equilibrium point E 3 are Since the existence of the equilibrium point E 3 is based on R 0 > 1, (2), we conclude that the equilibrium point E 3 is locally asymptotically stable provided that conditions (12) hold and (i) r q < β q L 3 and (ii) R 0 > 1.
The Routh-Hurwitz criteria necessary and sufficient for local asymptotic stability of the equilibrium point E 4 are that the Hurwitz determinants H i are all positive [32]. For a fifth-degree polynomial, these criteria are Thus we have the following result. (14) hold; otherwise, it is unstable.

Global stability of equilibrium points
In this section, we study the global stability of the equilibrium points E 1 , E 2 , E 3 , and E 4 determined earlier.
(a) Trivial equilibrium point E 1 . Let us consider the Lyapunov function As we can observe, the Lyapunov functional U 1

(M, E, L, P, A) is defined, continuous, and positive definite for all M(t), E(t), L(t), P(t), and A(t).
It is evident that U 1 vanishes at E 1 . The fractional derivative of U(t) along the solutions of system (1) leads to Therefore we have the following theorem.

Evidently, the function U 2 (M, E, L, P, A) is defined, continuous, and positive definite for all M(t), E(t), L(t), P(t), and A(t).
Furthermore, U 2 vanishes at E 2 . Hence the fractional derivative of U 2 (t) along the solutions of the system satisfies .
. Therefore we have the following theorem.

Theorem 3.5
The trivial equilibrium point E 2 is globally asymptotically stable if the following conditions hold: (c) Global stability of equilibrium points E 3 and E 4 . We will use the following Lyapunov function to investigate the global stability of the equilibrium points E 3 and E 4 : where a 1 , a 2 , a 3 , and a 4 are positive constants to be determined. Let Recall that at this equilibrium, we have the following identities: Setting Let (x) = 1x + ln x for x > 0. It follows that (x) ≤ 0 with the equality if and only if x = 1. Using this relation, we have Similarly, we can write

-
Therefore c t 0 D q t U 3 (t) is negative definite if the following conditions hold: Therefore we have the following theorem.

Theorem 3.6
The equilibrium points E 3 and E 4 are globally asymptotically stable if the following conditions hold:

Model parameterization
In this section, we present the baseline values for the model parameters. Majority of the parameter values are taken from previously published studies, and a few not available in literature were estimated within plausible and reasonable ranges so as to draw reasonably realistic scenarios.
(i) Natural mortality rate of adult FAW μ A : The life span of female adult FAW is 15-21 days. It follows that the natural mortality rate of the moth is (ii) Egg laying rate b and life span of adult moth μ -1 A : During its entire life span of 15-21 days and adult female FAW's total egg production per female averages about 1500 with a maximum of over 2000 [28]. The average daily egg laying rate can be expressed as follows: eggs laid per day = eggs laid in a lifetime expected lifetime .
Westbrook et al. [33] estimated that a female adult moth with a life span of 18 days can oviposit about 125 eggs. Hence in our simulations, we set b = 125 eggs per day and μ A = 1/18 per day. (iii) Egg hatching rate α E and gender ratio W : Mathematically, the egg hatching rate is the inverse of average duration of the egg stage, that is, Depending on the climate, the duration of egg stage takes an average period of 2-3 days [3,28]. Westbrook et al. [33] estimated a gender ratio of 50:50 males/females. (iv) Average duration of the larval stage α -1 L : The duration of the larval stage is influenced by climate changes. During summer periods, the larval stage is about 14 days and 30 days during cool weather [3,28,29].
(v) Average duration of the pupal stage α -1 P : Similarly to the larval stage, the pupal stage also depends on the climate. It is about 8-9 days during summer but reaches 20-30 days during the winter [3,28].

Sensitivity analysis of the reproduction number
Analytical results of the model have shown that the basic reproduction number is an important threshold parameter for the persistence and extinction of FAW during any outbreak. Since the parameters of the proposed model have either been drawn from literature or estimated, there is need to investigate the influence of each parameter on the magnitude of the basic reproduction number R 0 so as to understand the uncertainty regarding their values. To infer on the relationship between the model parameters and individual parameters, we conduct sensitivity analysis as follows.
Definition 4.1 (See [34]) The normalized sensitivity index of R 0 , which depends differentiably on a parameter, say κ, is defined by The model parameters whose sensitivity index is positive will increase the size of R 0 whenever they are increased, whereas those with negative index decrease R 0 whenever they are increased. It follows from (15) that the normalized sensitivity of R 0 with regard to the model parameters that define it is given by As we can observe from (16), the model parameters b, W , and α j (j = E, L, P, A) increase the size of R 0 whenever they are increased, whereas the model parameters μ j and u j decrease the size of R 0 whenever they are increased. It is worth noting that an increase in either b or W by 10% may result in an increase in the magnitude of R 0 by 10%. However, an increase by 10% of α j increases the size of R 0 by a value less than 10%. In addition, note that u j has a negative effect on R 0 , implying that intervention strategy has an impact on extinction and persistence of FAW in the environment. Without loss of generality, we set u E = u L = u P = u A = 0.3 and computed the sensitivity index for each model parameter that defines R 0 . The results are presented in Table 1 and Fig. 2.
Numerical results in Table 1 suggests that pest control intervention strategies more effects on minimizing the FAW population in the field if such strategies target the adult FAW population. Simulation results in Fig. 3 demonstrate the effects of varying the intervention strategies on extinction and persistence of pests in the field. For simplicity, we set u = u E = u L = u P = u A , whereas the other parameter values are taken from Table 2. From the results we can note that any value of u > 0.5 leads to the extinction of the pest, and persistence of the pests occurs for u < 0.5.

Population level effects
In this section, we conducted additional simulations to numerically illustrate the dynamical behavior of system (1) and to validate the analytical results such as the stability of the equilibria. We used the fractional Adam-Bashforth-Moulton method [35] to conduct the  Table 2. The exact numerical indices are also presented in Table 1, whereas control parameters are set to u E = u L = u P = u A = 0.3. Overall, we can note that model parameters b and W are highly correlated with R 0 . In particular, an increase in either b or W by 10% may result in an increase in the magnitude of R 0 by 10% Figure 3 Effects of varying intervention strategies on the magnitude of R 0 . Here we set u = u E , u L = u P = u A . (a) Dynamics of R 0 on a wide ranges of values for u. (b) Zoomed graph, that is, an illustration of the dynamics of R 0 for a refined range of values for u. From the results we can note that any value of u > 0.5 leads to the extinction of the pest, and persistence of the pests occurs for u < 0.5 simulations, that is, for a differential equation the fractional variant of the one step Adam-Moulton method is given by where t i = ih with some fixed h, and To determine the error in this method, by assuming that t i = ih = iτ N with some N ∈ N, we have (see [35]) Simulating system (1), we assumed the following initial population levels: E(0) = 1000, L(0) = P(0) = 0, A(0) = 500, and M(0) = 15. Numerical results in Fig. 4 illustrate the dynamics of the pest and maize biomass whenever the reproduction ratio R 0 is less than unity. As we can note, if the moth cannot produce more than one off-spring, then within a period of 200 days, all the FAW populations (eggs, larvae, pupae, and moth) will become extinct, whereas the maize biomass will increase with time till it reaches the expected maximum biomass per plant (50 kg plant -1 ). We can also observe that the convergence of solutions to their respective limiting points in time depends on the fractional order q: as q approaches unity, the time taken by solutions to converge to the limiting point increases. From the simulation results shown in Fig. 5 we can observe that whenever each female moth reproduces more than one off-spring, that is, R 0 > 1, then the pest population will persist in the field till the final harvesting time t = 300 day. In addition, the final maize biomass per plant will be less than the expected 50 k,plant -1 . Precisely, maize biomass increases from the start and reaches a maximum of 50 kg plant -1 after approximately 100 days, and after that, it decreases gradually till it stabilizes at approximated 18 kg plant -1 . Figure 6 shows the solutions of model system (1) for an experiment set up with small population sizes for the pest, that is, E(0) = 100, L(0) = P(0) = 0, and A(0) = 50, together with a control rate of u E = u L = u P = u A = 0.45 day -1 , leading to R 0 = 1.3583. Furthermore, q was fixed to 0.8. As we can observe, the pest population increases rapidly within the first 100 days, and then it stabilizes. The maize biomass also increases during the first 50 days and attains a maximum approximately close to the expected value 50 kg plant -1 , and then the biomass decreases gradually for approximately 50 days before it becomes stable at approximately 50 kg plant -1 . Overall, the egg population will dominate all the pest populations.
Numerical results in Fig. 7 depicts the effects of a FAW outbreak with a large initial pest life cycle population, E(0) = 2000, L(0) = P(0) = 0, and A(0) = 15, combined with less effective control measures, that is, u E = u L = 0.45, and u P = u A = 0. More often, pesticides,  Table 2. In addition, we set u E = u L = u P = u A = 0.52 to obtain R 0 = 0.8630. We can note that whenever R 0 < 1, the pest population becomes extinct while maize biomass increases with time, reaching its expected maximum of 50 kg plant -1 which are known to effectively control FAW, are expensive such that farmers in some areas rely on traditional methods of controlling the pest such as hand picking of caterpillars, picking and destroying egg masses, spraying lime, salt, oil, and soap solution. Prior studies suggests that traditional methods are less effective and are likely to eliminate the egg and larvae population. Hence, in Fig. 7, we explore the effects of a FAW outbreak with a large initial pest life cycle population coupled with less effective control measures. As we can observe, an outbreak with a large pest population coupled with less effective control measure may result in the pest population increasing rapidly so that in less than 100 days, they will reach their respective maximum. After an initial increase, the maize biomass would gradually decrease to a level below its initial biomass. The results highlight the importance of effective control measure on increasing maize biomass whenever there is a FAW outbreak.  Table 2. In addition, we set u E = u L = u P = u A = 0.35 to obtain R 0 = 2.8931. As we can observe, the solutions suggest that the pest will persist in the field and the final maize biomass will be less than the expected 50 kg plant -1

Concluding remarks
In this study, we presented a Caputo fractional-order model for fall armyworm (Spodoptera frugiperda) infestation in a maize field. Fall armyworm (Spodoptera frugiperda) commonly known as FAW remain a major pest of maize. The pest has already been considered as one of the greatest threat to food security in Africa despite the fact that it was first detected in 2016. The pest is highly destructive and fast spreading. Moth are capable of flying up to 100 km in one night. Based on its destructiveness, there is need to gain understanding the dynamics of this pest and the maize plant whenever there is an outbreak. Through the use of mathematical models, it is possible to predict many real-world problems in various fields such as agriculture, economics, biology, engineering, and so on. In particular, mathematical models are capable of providing solutions to phenomena that are difficult to measure in the field. Here fractional derivatives model have been utilized to model the dynamics of FAW infestation in a maize field based on the fact that fractional calculus is naturally related to many adaptive systems with memory and hereditary properties, which widely exist in several fields such as biology, agriculture, medicine, physics, chemistry, and engineering [21]. Mathematical analysis of the proposed model reveals that there exists a threshold parameter, the basic reproduction number, which governs the persistence and extinction of FAW in the field. Biologically, this basic reproduction number represents the average number of newborns produced by one individual female moth during its life span. We have noted that if one female moth is not capable of producing more than one off-spring, then the pest population becomes extinct; otherwise, it persists. The basic reproduction number was qualitatively and quantitatively used to investigate the local and global stability of the model steady state. For two steady states, the model has the pest-free equilibrium and the pest persistence equilibrium. We have observed that both are locally and globally stable. In particular, the pest-free equilibrium point is both locally and globally stable whenever the basic reproduction number is less than unity. However, when the basic reproduction number is greater than unity, there exists a pest persistence equilibrium point, which is also both locally and globally stable.
We have also noted that the model parameters, the egg laying rate and proportion of female moth in the environment, have a strong positive influence on increasing the size of the basic reproduction number. Precisely, increasing in the proportion of female moth results and to explore optimal intervention levels essential to minimize persistence of the pest population. We also used numerical simulations to illustrate the impact of initial pest population level during an outbreak on maize growth in a field.
The proposed model is not exhaustive. In the future work, we will explore the effects of temperature and seasonal variation on the dynamics of FAW and its implications on maize growth. In the current study, we found that there is a need for better metadata in plant population studies to help explain calibration and validation of proposed models.
Although we did not manage to validate the proposed model with data, due to its unavailability, the proposed model and results will certainly improve the existing knowledge on FAW dynamics and its implications in maize crops.