Bifurcation Analysis in a Coffee Berry-Borer-and-Ants Prey–Predator Model

: One of the most important agricultural activities worldwide, coffee cultivation, is severely affected by the Coffee Berry Borer (CBB), Hypothenemus hampei , considered the primary coffee pest. The CBB is a tiny beetle that diminishes the quantity and quality of coffee beans by penetrating them to feed on the endosperm and deposit its eggs, continuing its life cycle. One strategy to combat CBBs is using biological control agents, such as certain species of ants. Here, a mathematical model (consisting of a system of nonlinear ordinary differential equations) is formulated to describe the prey–predator interaction between CBBs and an unspecified species of ants. From this mathematical perspective, the model allows us to determine conditions for the existence and stability of extinction, persistence or co-existence equilibria. Transitions among those equilibrium states are investigated through the maximum per capita consumption rate of the predator as a bifurcation parameter, allowing us to determine the existence of transcritical and saddle-node bifurcations. Phase portraits of the system are presented for different values of bifurcation parameter, to illustrate stability outcomes and the occurrence of bifurcations. It is concluded that an increase in bifurcation parameters significantly reduces the CBB population, suggesting that ant predation is an effective control strategy, at least theoretically.


Introduction
Over the years, coffee has become a globally significant agricultural product and, consequently, a driver of economic growth in coffee-producing countries.According to statistics from 2020, approximately 125 million people worldwide depend on coffee production and/or trade for their livelihoods and quality of life.Importantly, coffee is a universally popular drink among consumers, with over US$50 billion in retail sales a year [1].Furthermore, in [2], it is stated that coffee is the most significant stimulant beverage on the planet.
Unfortunately, coffee cultivation faces a global threat: namely, the Coffee Berry Borer (CBB), Hypothenemus hampei (Ferrari, 1867) (Coleoptera: Curculionidae: Scolytidae), a beetle about the size of a pinhead that penetrates coffee beans to feed and reproduce within them, thereby reducing both the quantity and the quality of production.Its presence and impact on coffee crops worldwide designate it the primary coffee pest and pose a significant concern to the coffee industry [3][4][5].In fact, ref. [6] asserts that without proper plantation control measures, up to 100% of coffee beans can fall victim to CBB infestation.
The deterioration process begins when the adult female CBB reaches the coffee bean and pierces it at the navel to access the endosperm, which serves as both nourishment and a site for constructing galleries to complete its life cycle: egg, larva, pupa, and adult, with the larval stage being the most voracious consumer of the endosperm [4,7,8].Only adult female CBBs emerge from the bean to disperse and lay eggs in healthy beans, as males solely fulfill reproductive functions [4,8].The damage inflicted upon the coffee bean diminishes the quality of the product and, in many instances, results in complete bean loss due to premature shedding from the tree [9].
Due to the severity of the damage caused by the CBB in coffee plantations worldwide, its control represents a significant and ongoing challenge for the coffee industry.Addressing this challenge necessitates developing strategies to mitigate the effects of this beetle, with the goal of preventing or at least reducing coffee field infestations.Unfortunately, the CBB has proved to be a difficult insect to control, owing to its cryptic nature [6] and the protection provided by the coffee bean, which serves as the primary habitat for the CBB throughout most of its life cycle [8,9].
Among the various control strategies for the CBB, one approach is the application of insecticides, referred to as chemical control.However, according to [7,9,10] this method is not recommended as the sole means to mitigate the impact of this pest and its use may be considered erratic due to several drawbacks: (1) Insecticide sprays are only effective if applied timeously, specifically when the CBB is in flight or penetrating the coffee bean.These are the only moments when the CBB can come into contact with the chemical product, as once the pest reaches the endosperm, no insecticide offers satisfactory effectiveness.
(2) Insecticides are associated with environmental contamination issues that also affect beneficial fauna responsible for pest regulation within the coffee plantation.This can lead to ecological imbalances and the resurgence of secondary pests.(3) There is a risk of the development of insecticide resistance by the CBB, posing challenges for its control.Once the CBB reaches the endosperm, it can only be managed through timely coffee harvesting, considered a part of cultural control, or through the actions of other animals that may enter and attack the CBB within the bean, which constitutes a component of biological control [9].
In general, biological control or biocontrol practices focus on two primary approaches: classical biological control (i.e., the introduction of exotic enemies against exotic or native pests) and conservation biological control (i.e., the protection and enhancement of existing biological control agents within the agro-ecosystem) [7].In the specific case of the CBB, biocontrol can be carried out by promoting the beneficial fauna of coffee plantations or by introducing biological enemies of the pest, such as parasitoids, entomopathogens, and predators [9].
Predators make up the most numerous and diverse group of natural enemies of the CBB.This is because the role of a predator can be fulfilled by various organisms, such as beetles [11], lizards [8], and birds [12].However, a CBB predator that deserves special attention is the ant (Hymenoptera: Formicidae).Unlike the previously mentioned predators, ants naturally inhabit coffee plantations, and some of them can enter coffee beans, where they find the CBB in all its biological stages [9,13].
Ants play a significant role as natural regulators of CBBs and can be considered an essential component of an Integrated Borer Management (IBM) strategy under an environmentally friendly agro-ecological approach [10,[13][14][15].However, few studies have been conducted to quantify their potential as CBB predators and, in many cases, coffee growers are unaware of the benefits provided by ants, in terms of pest control services in coffee plantations [10].
Mathematical models also play a crucial role in understanding and designing effective control strategies to mitigate the spread and impact of the CBB in coffee plantations.In the literature, several studies have been reported that use mathematical models, in terms of differential equations to describe the population dynamics of the CBB, such as [6,16,17], which investigated the interaction of the CBB with coffee beans without the intervention of other species.However, what is of particular interest are those works that propose mathematical models in terms of differential equations but describe the interaction of the CBB with some form of biocontrol.Below, a brief overview of some of these studies is presented.
In [2], a deterministic mathematical model was formulated, consisting of a system of nonlinear ordinary differential equations, to describe the infestation dynamics of the CBB in interaction with coffee beans in plantations and with aware farmers.The proposed model was subsequently extended to an optimal control problem, where one of the control variables was the effort made to reduce the number of colonizing CBB females using a biological control agent-in this case, the entomopathogenic fungus Beauveria bassiana.
In [18], a system of nonlinear ordinary differential equations was proposed to model the infestation dynamics of the CBB in interaction with coffee beans in plantations.In this work, the CBB was divided into two populations: flying fertilized females looking for their host, and the infesting females, which corresponded to the females that lay eggs inside the beans.The formulated model was then extended to an optimal control problem, specifically, biological control based on the use of the entomopathogenic fungus Beauveria bassiana.
In [19], the interaction of the CBB with three types of parasitoids was modeled using systems of ordinary differential equations.Such parasitoids correspond to wasps, and are Cephalonomia stephanoderis, Prorops nasuta, and Phymasthichus coffea.Conclusions regarding the most effective parasitoid were drawn from numerical simulations.
In [5], the population dynamics of the CBB when preyed upon by ants were studied.In particular, modifications in dynamics were investigated based on variations in the death rate of adult CBBs due to factors other than predation-factors such as cultural and/or chemical control methods.
Considering the advantages of ants as biological control agents of the CBB, as previously discussed, this work presents a mathematical model based on nonlinear ordinary differential equations to describe the prey-predator interaction, where the CBB is consumed by an unspecified species of ants.The main objective of this research was to analyze the impact of the maximum per capita consumption rate of the predator on biocontrol efficiency.This rate corresponds to a parameter of the Holling Type II functional response incorporated into the model.Specifically, the study aimed to investigate the population dynamics of the CBB and the modifications in the system's dynamics as this rate is increased.In doing so, this study contributes academically to CBB management because, as noted by [20], studying the pest's population dynamics over time helps determine opportune moments for control.Moreover, this approach seeks to make a significant contribution to the limited literature regarding the predatory potential of ants on the CBB, as emphasized by [10].
This study, along with the article by [5], represents a novelty in modeling the interaction of the CBB with a beneficial predator for coffee growers, such as ants.However, it is important to highlight that the focus of this work is on analyzing the effects of predation by ants while excluding external factors.Although the specific species of predatory ants is not specified, this does not diminish the merit of the research.As [10] points out, the true predatory potential of ants on CBB populations under field conditions in coffee plantations has not been accurately quantified because most studies have been conducted in controlled laboratory conditions, and because the few field studies carried out have not specified the ant species involved.
The rest of the paper is organized as follows: In Section 2, we present the model that describes the prey-predator interaction between CBBs and ants of an unspecified species.In Section 3, the main properties of the model are proven.Some numerical simulations are shown in Section 4, and in Section 5 we explain the ecological meanings of the obtained analytical outcomes.

Biological Background
One of the key elements in the prey-predator relationship is the predator's functional response or consumption function, which describes how the predation rate (i.e., the rate at which predators consume their prey) varies depending on prey density [21].In particular, the Holling Type II functional response assumes that the predation rate is proportional to the encounter rate between prey and predators but also takes into account predator satiation as prey density increases.In other words, the Holling Type II functional response in a prey-predator model captures the following observed behavior in nature: predators cannot consume prey indefinitely at constant rates, but their efficiency decreases as prey density increases [22].One mathematical expression of the Holling Type II functional response is as follows: where the variable B and the involved parameters, according to [23], have the following meanings: • The variable B represents the prey density.In the context of the prey-predator model, prey density is the quantity of prey present in a specific area.

•
The parameter α represents the maximum per capita consumption rate of the predator, i.e., the maximum number of prey an individual predator can consume in each unit of time.

•
The parameter η represents the number of prey required to reach half of the maximum per capita consumption rate of the predator α.This parameter reflects the point at which predators begin to saturate and become less efficient at capturing additional prey.
In some works, it is referred to as the half-saturation constant, as in [2], for example.

Model Formulation
We consider a prey-predator interaction modeled by a system of nonlinear ordinary differential equations, where the role of the predator is played by ants, which have among their dietary habits the consumption of CBB.It is assumed that both populations are homogeneously distributed (position is irrelevant), and migration phenomena are not considered.Additionally, it is assumed that the described dynamics correspond to a finite period, where populations are not abruptly affected by factors such as coffee harvesting or crop change.
For the population of the CBB, denoted by the letter B, logistic growth is considered, with ω representing the intrinsic growth rate and q representing the carrying capacity.It is also assumed that the CBB population is reduced by the predatory effect of ants.According to works such as [4,5,7,10,13,15], different ant species feed on CBB adults and/or immatures (eggs, larvae, and pupae).In this case, it is assumed that the number of "consumptions" is proportional to the average number of CBB, regardless of their state, but we also take into account the satiation of predators as prey density increases.This situation is modeled using the Holling Type II functional response, which, according to [22], is generally exemplified by invertebrate predators, although [8] uses it in the same context of CBB being consumed by lizards of the genus Anolis.The parameters α and η are involved in this functional response, where the former corresponds to the maximum per capita consumption rate of the predator, while the latter is the half-saturation constant.Although η << q because, in this case, q is the carrying capacity of the CBB population, it is important to note that if η → q, a larger number of prey is needed to reach α/2 [23].Finally, the death of CBB due to factors other than predation is not considered, although it is important to remember that natural CBB mortality is implicit in the intrinsic growth rate ω.In accordance with the above, the population dynamics of the CBB are modeled by the following differential equation: For the ant population, denoted by the letter R, logistic growth is also assumed, with r representing the intrinsic growth rate and k representing the carrying capacity.Furthermore, an additional increase in the intrinsic growth of ants, attributed to predation, is assumed-a hypothesis incorporated in population dynamics works such as [24][25][26][27] in the context of Aedes aegypti and in [5] within the context of CBB.As is customary in prey-predator models, it should be noted that not all prey biomass is converted into predator biomass, so a biomass conversion rate denoted by ε is considered.The only inhibiting factor for the ant population's growth is assumed to be natural death, which is implicit in the intrinsic growth rate r, ignoring the effects of potential natural enemies or death due to insecticides.In accordance with the above, the population dynamics of the ants are modeled by the following differential equation: The system consisting of the two previous differential equations represents the model that, based on the considerations previously mentioned, describes the prey-predator dynamics between the CBB and an unspecified species of ants.Below, we present this system: with initial conditions B(0) = B 0 > 0 and R(0) = R 0 > 0, where all the parameters are positive and are described in Section 4, corresponding to the Numerical Results.

Main Results
The main results of model ( 1) are established in the form of propositions, which are stated and proved.In these propositions, the region of invariance, the biological sense of equilibrium points, and their stability are analyzed.In Proposition 1, the region of invariance Ω of model ( 1) is determined.It is an area of biological interest that is positively invariant with respect to the flow of system (1), meaning that any solution of system (1) that starts in Ω remains there for all t ≥ 0. Proposition 1.The set Ω provided below is a positively invariant region: (2) Proof.First, it is shown that the B-axis, the R-axis, and the horizontal straight line with equation R = k are invariant sets of system (1).Indeed, by setting B = 0 in the first differential equation and setting R = 0 and R = k in the second differential equation, we obtain This means that if non-negative initial conditions are considered (so that they make biological sense) and are located on any of these three sets (straight lines), the trajectories will remain within that set.Now, replacing B = q in the first differential equation results in This means that trajectories with non-negative initial conditions on the straight line B = q enter the region Ω. Figure 1 illustrates the results obtained in this proof.Note that in model ( 1) the absence of predation (α = 0) is a condition that guarantees the co-existence of species, as both populations acquire logistic dynamics with parameters ω > 0 and r > 0. However, the goal is to analyze the effect of the maximum per capita consumption rate of the predator α on the efficiency of biological control by ants on the CBB.To this end, we begin by determining the equilibrium points of system ( 1).An equilibrium point (equilibrium solution) is reached when the rates of change of all variables in the system are equal to zero.This means that at this point the variables are not increasing or decreasing, resulting in a constant and balanced population.To obtain the equilibrium points of model ( 1), we set the derivatives equal to zero, yielding the following: It is clear that equilibrium points E 1 , E 2 , and E 3 always have biological sense, but the same does not hold for E 4 and E 5 .Next, we show under which condition on the parameter α the radicand in (3) and ( 4) is non-negative.This is given by Therefore, the abscissas B4 (3) and B5 (4) of equilibrium points E 4 = ( B4 , k) and Now, the discussion of the biological sense of equilibria E 4 and E 5 focuses on the determining conditions that ensure the abscissas B4 (3) and B5 (4) are non-negative.This analysis is carried out in Proposition 2, varying the parameter α (the maximum per capita consumption rate of the predator), concerning the thresholds α * = ηω k and α * * = ω(q+η) 2 4kq , which satisfy the inequality α * < α * * .Indeed, since q >> η because q is the carrying capacity of the CBB population, then it is true that 0 < (q − η) 2 and, therefore, And by multiplying both sides by ω k , we obtain In summary, Proposition 2 establishes conditions in terms of the parameter α that determine the biological sense of the equilibrium points or their collision.
Proposition 2. If q > η, the biological sense of equilibria E 4 and E 5 in system (1) depends on the parameter α, the maximum per capita consumption rate of the predator, as follows: 1.
Next, we will state and prove Proposition 3, which analyzes the stability of equilibrium points E 1 , E 2 , and E 3 of system (1).Proposition 3. The stability of equilibrium points E 1 , E 2 , and E 3 of system (1) is as follows: 1.
Proof.The Jacobian matrix of the system (1) is given by The proof involves evaluating the Jacobian matrix (7) at each equilibrium point E i , i = 1, 2, 3, which is denoted as J(E i ), and calculating the eigenvalues in each case.For hyperbolic equilibrium points, stability conclusions are obtained based on the sign of the real part of the eigenvalues, following the methodology of the Hartman-Grobman theorem [28].For the nonhyperbolic equilibrium point, the Sotomayor theorem is used to demonstrate the occurrence of a transcritical bifurcation [28]: 1.
The eigenvalues of J(E 1 ) = J(0, 0) are λ 1 = ω and λ 2 = r.Since both eigenvalues are positive, we conclude that E 1 is an unstable equilibrium point of the source type.

3.
The matrix J(E 3 ) = J(0, k) has the following eigenvalues: As λ 2 < 0, the stability type of equilibrium point E 3 depends on the sign of λ 1 , which, in turn, depends on whether α is less than, equal to, or greater than α * .Thus, if α < α * then λ 1 > 0, and, therefore, E 3 is an unstable equilibrium point of the saddle type.And if α > α * , then λ 1 < 0 and, thus, the equilibrium point E 3 is a locally asymptotically stable sink.The case α = α * results in λ 1 = 0. Consequently, equilibrium point E 3 is nonhyperbolic, and the occurrence of the bifurcation is demonstrated by verifying that the conditions of Sotomayor's theorem for the case of transcritical bifurcation are satisfied [28].In this case, v = (1, 0) T and w = (1, 0) T are the eigenvectors associated with the zero eigenvalue of the matrices J(E 3 ) and J T (E 3 ), respectively.And the vector field associated with system (1) is from where On the other hand, since Finally, applying Then, Once conditions ( 9), ( 10) and ( 12) are verified, the proposition is proven.
From Proposition 2, it is known that the equilibrium point E 4 = ( B4 , k) makes biologi- cal sense if α ≤ α * * and that the equilibrium point E 5 = ( B5 , k) makes biological sense and does not collide with other equilibrium points if α * < α < α * * .Therefore, in Proposition 4, the stability of E 4 and E 5 is analyzed based on these same conditions.Proposition 4. The stability of equilibrium points E 4 and E 5 of system (1) depends on the value of α, the maximum per capita consumption rate of the predator, as follows: 1.
Proof.To simplify the writing of this proof, we denote ∆ as the radicand in ( 3) and ( 4), which is given by 1.

2.
When evaluating the Jacobian matrix ( 7) at E 5 , we obtain J(E 5 ) = J( B5 , k), the eigenvalues of which are given by 16) Substituting ( 4) and ( 13) into ( 16), we obtain Again, from ( 5) and ( 13), we know that if α < α * * ⇔ ∆ > 0 then the first factor in the numerator is positive.Next, we show that if α > α * then the second factor in the numerator is positive.
On both sides of this last inequality, we add the expression and then it is factored to obtain It has been shown that λ 1 > 0. Now, we substitute ( 4) and ( 13) into ( 17) to obtain Next, we show that if α > α * then the denominator is positive.To do this, we start with the following fact: allowing us to conclude that λ 2 < 0. Taking into account that λ 1 > 0 and λ 2 < 0, we conclude that equilibrium point E 5 is unstable, of the saddle type.
In Figure 2, the biological sense and stability results obtained in Propositions 2-4 are summarized: Bifurcation diagram summarizing the biological sense and stability results obtained in Propositions 2-4.The solid line and black circle indicate that the equilibrium point is locally and asymptotically stable, the dashed line and white circle indicate that the equilibrium point is unstable, and the green circle indicates that the equilibrium point is unstable, of the saddle-node type, as referred to in [28].

Numerical Results
To illustrate the biological sense conditions from Proposition 2 and the stability results discussed in Propositions 3 and 4, numerical simulations of model (1) were performed using the parameter values provided in Table 1.It is important to note that these parameter values are used for illustrative purposes and may not necessarily represent reality accurately [29].Additionally, while many of these values are sourced from the literature, they should not be considered standard values, as different research may use different figures.This can be exemplified by two studies published in 2021 that used the maximum number of CBBs per tree but reported significantly different values.
In [30], a mathematical model based on nonlinear ordinary differential equations describing the interaction between coffee beans and CBBs was analyzed.The population dynamics of the pest were partially governed by logistic growth, and the carrying capacity, defined as the number of CBBs per tree, was assumed to be 20,000 individuals.On the other hand, in [20], the population dynamics, dispersal, and colonization of CBBs in Colombia were studied and, through fieldwork, the number of CBBs per tree was estimated to be 2674.Despite the significant differences in the figures used in these two studies, this work assumes, for illustrative purposes, that the carrying capacity for the CBB population in model (1), denoted as q, is 2674 CBB per tree, as indicated in Table 1.
Now, let us focus on the carrying capacity of the ant population, which is challenging to precisely measure, since the majority of colony individuals remain in the nest, and those that forage do so in restricted time intervals [31].However, in [14], it was determined through observations in coffee plantations that the average number of ants per branch is five.On the other hand, in [32], fieldwork established that the number of branches per low-stature coffee tree (e.g., the Caturra variety) is 81.Based on this, it is estimated that the carrying capacity for the ant population, denoted as k, is 400 individuals per tree, obtained by approximating the multiplication of (five ants) × (81 branches), as shown in Table 1.

Parameter
However, the equilibrium points E 4 and E 5 do not have a fixed numerical representation because their abscissas, given in ( 3) and ( 4), respectively, depend on the bifurcation parameter α, whose value varies.Nevertheless, Table 2 shows some numerical representations of equilibria E 4 and E 5 for different values of the parameter α, providing a numerical summary of Proposition 2. The numerical simulations of model ( 1) consist of phase portraits that arise when varying the maximum per capita consumption rate of the predator α, with a particular focus on situations generated when this parameter crosses the bifurcation values α * and α * * .Figures 3-6, created using pplane8.m,an interactive computational tool in MATLAB software, depict the phase portraits of the system.These portraits exhibit an important generality: in all cases, equilibrium point E 1 is an unstable source and equilibrium point E 2 is an unstable saddle, as demonstrated in Proposition 3.
Figure 3 corresponds to the phase portrait of model (1) when α < α * .In this case, equilibrium point E 3 is unstable, of the saddle type, equilibrium point E 4 is locally asymptotically stable, and equilibrium point E 5 lacks biological sense because its abscissa is negative.In this scenario, the populations co-exist regardless of the initial conditions, and the highest number of CBBs is observed-specifically, approximately 2042 individuals, according to Table 2.This situation is associated with the lowest assumed value for the maximum per capita consumption rate of the predator α: namely, α = 0.06.In fact, as α → 0, E 4 → (q, k).For α = α * , the phase portrait is not shown, as it is similar to that presented in Figure 3.However, it should be noted that in this case, equilibrium point E 5 collides with E 3 , and the abscissa of equilibrium point E 4 is reduced.Specifically, for α = 0.077625, B4 corresponds to approximately 1774 CBB individuals, according to Table 2.
Figure 4 shows the phase portrait of model (1) when α * < α < α * * ; specifically, α = 0.095.In this case, equilibrium points E 3 and E 4 are locally asymptotically stable.However, the latter has a significantly reduced abscissa compared to the previous case; specifically, B4 = 1385 according to Table 2.The main novelty in this scenario is the biological sense of equilibrium point E 5 , which is also unstable, of the saddle type, in accordance with Propositions 2 and 4.   1) when α = α * * .In this case, equilibrium point E 3 remains locally asymptotically stable, and equilibrium point E 5 collides with E 4 , giving rise to a saddle-node bifurcation, as stated in Proposition 4. Equilibrium point E 4 becomes unstable, of the saddle-node type, a denomination given in [28], and now has a much-reduced abscissa: precisely, E 4 = (887, 400) for α = α * * = 0.1030022723.Finally, in Figure 6, the phase portrait of model (1) when α > α * * is shown, specifically for α = 0.11.This scenario is ideal for the coffee grower because it corresponds to the extinction of the CBB population.Indeed, equilibrium points E 4 and E 5 lose their biological sense and E 3 is locally asymptotically stable.This means that any initial condition within the region Ω given in (2) will lead to trajectories that asymptotically approach equilibrium point E 3 , except for initial conditions of the form (B(0), 0), i.e., those located on the B-axis where the ant population is zero.It is important to highlight that the trajectories in Figures 3-6 provide a visual representation of the system dynamics and that the initial conditions have been carefully chosen to ensure that the associated trajectories illustrate the local and asymptotic stability or instability of the different equilibrium points.Furthermore, the different trajectories highlight the system's sensitivity to initial conditions and how small variations in these can have a significant impact on its long-term behavior.In the particular case of Figures 4 and 5, where the dynamics of system (1) are illustrated when α * < α ≤ α * * , the initial concentrations of CBBs and ants play a crucial role in the evolution of the CBB population concerning its extinction.This aspect is intrinsically related to the configuration of the stable manifold and the basins of attraction defined by equilibrium points E 3 and E 4 .The stable manifold of E 5 in Figure 4 and the stable manifold of E 4 (E 4 = E 5 ) in Figure 5 define these basins of attraction.Depending on the initial position with respect to these stable manifolds, the system trajectories may converge towards one of the equilibrium points that have been identified as locally and asymptotically stable, determining the persistence or extinction of the CBB population.This analysis underscores the importance of considering initial conditions in the planning of control strategies for CBB in coffee plantations.

Discussion
Mathematical models are relatively sophisticated representations of reality, but they are not exact.Their utility lies in their role as substitutes for scientific experiments.In the field of population ecology, mathematical models play a fundamental role as they are crucial for developing theoretical foundations related to the factors that influence population size [34].This importance is evident in the predator-prey model formulated and studied in this work, as it establishes conditions in terms of the parameter α that determine the co-existence of the two considered populations: CBB and an unspecified species of predatory ants, or the extermination of the pest.
The parameter α corresponds to the maximum per capita consumption rate of the predator, and it is clear from Table 2 that increasing its value disadvantages the CBB population.Indeed, by increasing the value of α, the CBB population located at the abscissa of equilibrium point E 4 decreases.The above statement can be corroborated by Figures 3-6.Note that when the value of the parameter α is low, specifically α < α * , Figure 3 indicates that the populations co-exist, and the same situation occurs for α = α * , although this scenario is not shown.However, for α * < α ≤ α * * , Figures 4 and 5 show that increasing α starts to be useful.In these cases, some initial conditions lead to trajectories that asymptotically approach equilibrium point E 3 , highlighting the possibility of the CBB population going extinct.Fortunately, a sufficient increase in α, specifically α > α * * , guarantees the extinction of the CBB population, as observed in Figure 6, where the only locally asymptotically stable equilibrium point is E 3 .This means that any initial condition in Ω where the ant population is not zero will lead to trajectories that asymptotically approach equilibrium point E 3 .
Taking into account the previous discussion, it is concluded that an increase in the bifurcation parameter α significantly reduces the CBB population, suggesting that ant predation is an effective control strategy, at least theoretically.In other words, an effective strategy to mitigate the impact of CBBs in coffee plantations is to stimulate an increase in the ant's predation rate.This can be achieved by eliminating potential food sources in the crops, as assumed in the proposed model, where ants obtain nutrients from sources other than CBB.
An important aspect to highlight in the results of this work is that regardless of whether the populations co-existed or the CBB population went extinct, the ant population always tended to its carrying capacity, as occurs in reality.In other words, according to the stability results of Propositions 3 and 4, and by Figures 3-6, the equilibrium points that were locally and asymptotically stable, i.e., E 3 and E 4 , had an ordinate of k.This means that in any stability situation, the ant population would always tend to its carrying capacity, corresponding to a natural logistic growth behavior.This situation arose because the conversion of CBB biomass into ant biomass was modeled as an extra increase in intrinsic growth attributed to predation.Hence, from the obtained results, it can be suggested that the second differential equation of system (1) is a better choice in modeling than the usual form dR dt = rR 1 − R k + ε αB B+η R. In fact, if this latter differential equation is incorporated into model (1), the equilibria E 4 and E 5 no longer have the ordinate k.
As a future work, it is proposed to incorporate other types of control for CBBs, such as cultural control when the presence of ants with a predatory role is not sufficient to mitigate the effects of CBBs, i.e., when the maximum per capita consumption rate of the predator α satisfies α ≤ α * * .Cultural control aims to minimize the availability of food and shelter for the pest through the implementation of various manual practices, such as frequent and efficient harvesting of ripe beans in the plantation or the use of homemade traps with alcoholic attractants [9].Assuming that this capture is carried out based on the quantity of CBB in the crop, we can consider the application of cultural control if the amount of CBB in the coffee field exceeds or does not meet a certain threshold.In this way, the system dynamics can be modeled with a more comprehensive approach that takes into account both the effects of biological control and cultural control, which can be achieved with PieceWise Smooth (PWS) dynamical systems.

Figure 5
Figure 5 corresponds to the phase portrait of model (1) when α = α * * .In this case, equilibrium point E 3 remains locally asymptotically stable, and equilibrium point E 5 collides with E 4 , giving rise to a saddle-node bifurcation, as stated in Proposition 4. Equilibrium point E 4 becomes unstable, of the saddle-node type, a denomination given in[28], and now has a much-reduced abscissa: precisely, E 4 = (887, 400) for α = α * * = 0.1030022723.

Table 2 .
Numerical representation of equilibrium points E 4 and E 5 of system (1).