A theoretical model for optimal control of banana Moko ( Musa AAB Simmonds)

A population simulation model with non-linear ordinary differential equations is presented, which interprets the dynamics of the banana Moko, with prevention of the disease and population of susceptible and infected plants over time. A crop with a variable population of plants and a logistic growth of replanting is assumed, taking into account the maximum capacity of plants in the delimited study area. Also, with the help of farmers, the costs of implementing prevention strategies and elimination of infected plants were calculated per week in order to determine the optimal conditions that control the disease and reduce production costs. We found that the implementation of prevention strategies (f) plays an important role, but the parameter that most influences the threshold value is the elimination of infected plants g . However, to reduce production costs due to the high implementation of prevention strategies and to maintain the disease in a controlled state, both controls u 1 and u 2 should be implemented between 40% and 60%, obtaining with this percentage an approximate reduction of 51.37% in production costs per week, where in 23 weeks following the same conditions it is expected to have a healthy plantation without infected plants.


Introduction
The banana is a fruit of great economic importance and food sovereignty, because it is found in the shopping basket of people across different social strata and because of its nutritional content.However, its production is threatened by re-emerging diseases such as Moko, caused by the bacterium Ralstonia solanacearum race 2 philotype II (Fegan & Prior, 2006), which causes wilting and deterioration of the plant.Symptoms are usually visible when there is a great spread of bacteria and adjacent plants may have already been infected (Viljoen et al., 2016).Moko is a peculiar manifestation of bacterial wilt; it is a quarantine pest that, once inside the host, moves through the vascular bundles.Being a vascular disease, this bacterium not only affects the vegetative part but also the daughter hill, promoting the spread of the disease, which is accelerated by the high optimal minimum, optimal and maximum temperatures of 10°C, 35°C and 41°C, respectively, infecting triploid plantains, heliconia (Heliconia spp.) and other ornamental Musacea plants (Jiang et al., 2017).
The development of mathematical models has contributed through the use of a wide range of techniques to the study of epidemics and diseases, helping to answer biological questions and raising new questions related to the epidemiology and ecology of pathogens and the diseases they cause; in most cases, mathematical models lead directly to applications in the control of the disease (Jeger et al., 2018).In plants, some prediction models have evaluated the control of some diseases in plantain and bananas, such as banana wilt by Xanthomonas (BXW) using mathematical models, describing a deterministic SI-type epidemic model for control of BXW focusing on the integrated management of the disease through cultural control as in Nannyonga et al. (2015), who considered the optimal control strategies associated with the prevention of transmission by the use of contaminated tools.The researchers assumed a model with three modes of transmission: vertical (from the mother plant to its child), horizontal (indirect) from the vector to plant, and through contaminated agricultural tools (Nannyonga et al., 2015).Likewise, Nakakawa et al. (2016) presented a mathematical model for BXW propagated by an insect vector.The mathematical model they formulated takes into account inflorescence infection and vertical transmission from the mother corm to the daughter hills, but not tool-based transmission by humans (Nakakawa et al., 2016).In this context, a dynamic system is formulated based on ordinary two-dimensional differential equations that interprets the dynamics of incidence of banana Moko disease, including prevention and treatment.
In the research carried out by (Bautista-Montealegre et al., 2016), the state of the disease in 2016 is shown, these authors to contribute to the management of banana Moko disease in the department of Quindío-Colombia, evaluated the relationship between the incidence of the disease and variables related to physical and chemical properties of the soil, as well as the use of the soil and the altitudinal location in 269 farms analyzing soils and foliar tissues, as well as the symptoms of the disease to establish the effect of the variables on the probability of occurrence of the disease, finding a positive and significant correlation between the incidence of the disease, the hydraulic conductivity and the saturation of potassium in the soil; and negative and significant with the altitude, foliar copper concentration and presence of associated crops; Likewise, they argue that 10 of the 12 municipalities in the department have a high percentage of the disease, demonstrating the inadequate phytosanitary management that is still being carried out (Bautista-Montealegre et al., 2016).

The model
Mathematical models are a tool of a growing scientific branch and of a notorious and marked interdisciplinary nature, linking mainly biologists and mathematicians, but also researchers from other areas with the challenge of applying mathematical techniques to the study of biological processes (García-Macías & Ubertini, 2019).Previously, they were used mainly in epidemiology in SI, SIS, SIR, SIRS epidemic models among others; at present, mathematical models are, therefore, a useful tool in biology, agronomy, phytopathology, chemistry, environment and among many other areas, since they allow to make a representation of a biological system,

Amendments from Version 1
In this version we include the evaluators' suggestions; The explanation of the mathematical model that was not clear was included, the sensitivity and stability analysis was organized into independent items and the optimal control problem was included with the new graphs complementing the results.This improved version is longer but more complete.
Any further responses from the reviewers can be found at the end of the article REVISED the behavior of a certain disease etc., and with them facilitate the understanding of its dynamics in order to make predictions for future decisions on actions that facilitate control (Jeger et al., 2018).
The model presents the following assumptions: • In the plantation the total plant population is assumed to be positive.
• It is considered a plantation with a maximum capacity of banana plants k.
• It is considered that the disease of the banana Moko, follows a model Epidemic type SI (Susceptible-Infectious).
• Banana plants in asymptomatic state are not considered.1.
The model parameters are: γ, constant overseeding rate; k, load capacity (maximum capacity) of banana plants in the study region; and β, probability of transmission of infection.Preventive controls are: g, fraction of infected banana plants removed; and f, fraction of susceptible banana plants that receive prevention of contagion of the bacteria.The dynamic system that interprets the infectious process including prevention and elimination, is formed by the following two nonlinear differential equations: The region of eco-epidemiological sense is defined where the trajectories of the plant infection dynamics make sense, Similarly, in the second equation the derivative dy dt represents the variation of the population of infected plants with respect to time in weeks t, given by the inflow the incidence of banana plants minus the outflow the population of plants infected deleted gy.

Stability analysis
We start by finding the equilibrium populations, the constant solutions of the system, where the population variation of susceptible plants and variation of infected plants become zero, that is, 0; 0 dx dy We solve this non-linear algebraic system for x and y, determining the equilibrium point according to agricultural conditions.
From Equation (4) y = 0 o Substituting y = 0 in Equation (3), we obtain Of which x = 0 or x = k.Therefore, we have the point E 0 = (0,0), which does not make agricultural sense, since P > 0 and the equilibrium point E 1 = (k,0), disease free and in which the susceptible population equals the capacity maximum.
Since P = x + y, the total plant population in equilibrium is, ˆˆˆ.
Therefore, ˆ( 1).P k ρ = − ξ 0 , indicates the average number of infected plants that an infected plant produces during the infectious period (before being killed) in the population of susceptible plants and is considered the threshold of infected plants.We can consider this threshold as a function that depends on f and g, To determine the stability of each equilibrium point (E), we apply the Hartman-Grobman theorem (Perko, 2008), linearizing the system of non-linear Equation (1) -Equation ( 2), obtaining the linearization matrix (Jacobian matrix) of the form: With the following partial derivative elements, These elements of the matrix J (E) are the coefficients of the linear system Where, U = (u, v) t (transposed vector).
We conclude that the free equilibrium point of Moko disease is locally and asymptomatically stable.
For case 3 ˆ( , ), = E x y in matrix ( 7) we obtain the trace and the determinant of J(E 3 ), respectively, We conclude that the equilibrium point with susceptible plants and infected plants is locally and asymptomatically stable if the threshold inequalities (ξ 0 > 1 and ρ < 1) and the inequalities are met, These analytical results are shown in the phase planes of Figure 2, made with Maple 18 software (free trial available; SageMath is an openly available alternative), for different scenarios varying initial conditions.

Sensitivity analysis
The local sensitivity analysis was performed that is a measure of the relative change in a variable when its parameters change (Chitnis et al., 2008;Hamby, 1994;Rodrigues et al., 2013).That is, With respect to the values of the threshold indices of the disease and the Figure 3., the following observations are made: • For f = 0.5 it is true that • When the percentage of plants receiving prevention increases, the disease threshold decreases, that is, they are inversely proportional.This behavior is good in managing the disease.
• For values of f,g < 0.5, the disease threshold index is lower for the elimination of infected plants (g) and in the case that f,g > 0.5, the disease threshold index is lower in the case of prevention of susceptible plants.
• The disease threshold increases proportionally with respect to the transmission probability.
It is concluded that mathematical simulation models are a useful tool for research in banana Moko disease.
With them it was determined that the elimination of banana plants infected with the disease plays an essential role in the good agronomic management of the crop.

Optimal control problem
An objective, quadratic and cost functional linked to a system is presented of nonlinear ordinary differential equations, which interprets the dynamics of banana Moko (Figure 4), including a constant rate of disease prevention in the population of susceptible plants over time.It assumes a variable population of plants and a logistic growth of replanting having taken into account the maximum capacity of plants in the study region (Cherruault & Gallego, 1985;Louadj et al., 2018).The variables and parameters of the optimal control problem are described in Table 1.
The functional objective of direct and indirect costs is proposed:  Linked to the system of differential equations: It is about finding optimal control 1 2 ( ( ), ( )) u t u t such that: Where, is the space of admissible controls and L 2 is the space of integrable functions, and T is the control terminal time.

Optimal control problem analysis
The Hamiltonian function or (Pontryagin function) is of the form: x u λ where x = (x, y) is the vector of state variables, u = (u 1 ; u 2 ) the vector of controls, λ = (λ 1 , λ 2 ) the vector of attached or conjugate variables and L is the Lagrangian.That is to say, Applying the first order condition ∂ the optimal control is obtained: The conjugate system (or adjoint system) has the form: Where, (1 ) ( ) (1 ) ( ) (1 ) ( ) The contour problem is formed by the system of state variables of the dynamics of the Moko with their respective initial conditions, the conjugated system and the terminal conditions and the optimal controls:

Results and
Numerical analysis of the contour problem: with this analysis we can observe the decrease of the infected plants varying different conditions of the controls u 1 and u 2 , showing that when u 1 and u 2 are not implemented as shown in the figure at week 120, the population of infected plants tends to increase (Figure 5), and that the implementation of control u 1 in approximately 45% (Figure 6) and the implementation of control u 2 in approximately 60% (Figure 7), produces a complete decrease in infected plants at week 23.With this we can determine that if both controls are implemented in 40% and 60% respectively, the banana Moko disease in a crop with good agronomic management tends to disappear, likewise it contributes to reducing the production costs associated with the loss of plants due to infection and the costs of implementing prevention strategies by 51.13% weekly, which is equivalent to 60,756 Colombian pesos.
We conclude that in order to reduce production costs and maintain the disease in a controlled state, the recommended prevention strategies should be implemented, and with greater relevance the detection and rapid elimination of infected plants.

Reference Source
The authors made all the suggested changes.In my opinion, the article is a great scientific contribution to understand the transmission dynamics of banana moko and also establishes control measures to guide farmers in the measures to be taken.
The following paragraph lacks clarity.Please rewrite it.Some prediction models that calculate the propagation threshold R0 have evaluated the control of some diseases in plantain and bananas, such as banana wilt by Xanthomonas (BXW) using mathematical models, describing a deterministic SI-type epidemic model for control of BXW focusing on the integrated management of the disease through cultural control as in Nannyonga et al. (2015), who considered the optimal control strategies associated with the prevention of transmission by the use of contaminated tools.The researchers assumed a model with three modes of transmission: vertical (from the mother plant to its child), horizontal (indirect) from the vector to plant, and through contaminated agricultural tools (Nannyonga et al., 2015).

3.
The following sentences should be included in the section Model Local sensitivity is a measure of the relative change in a variable when its parameters change (Chitnis et al., 2008;Hamby, 1994;Rodrigues et al., 2013).That is….

4.
What do authors mean when they talk about agronomic sense? 5.
The results need to be expanded, including a better description.In the case of figure two, the authors do not present any description.

6.
Conclusions generally present the main findings.Please rearrange it to meet the highquality standards of the journal.

7.
The presentation of the research perspective seems messy and needs to be improved.The authors present a research perspective, but it is not considered in the paper's main body.The author should consider it as part of the analysis.The optimal control problem should be in the Section Model, describe the development in the results and the implications or expected results with this analysis in conclusion.

8.
The list of references requires an update to include recent studies.9.
Authors are encouraged to look for other models about diseases in banana crops and discuss their results to see advantages and differences.10.
Overall, the paper looks interesting, but it needs a major revision before any consideration.Reviewer Expertise: Public health, environment, sustainability, environmental management, econometrics I confirm that I have read this submission and believe that I have an appropriate level of expertise to confirm that it is of an acceptable scientific standard, however I have significant reservations, as outlined above.
The benefits of publishing with F1000Research: Your article is published within days, with no editorial bias • You can publish traditional articles, null/negative results, case reports, data notes and more • The peer review process is transparent and collaborative • Your article is indexed in PubMed after passing peer review • Dedicated customer support at every stage • For pre-submission enquiries, contact research@f1000.com

Figure 2 .
Figure 2. Local stability of the susceptible plant population and infected plant population corresponding to ξ 0 = 7 y ξ 0 = 0.79.

Figure 3 .
Figure 3.In the graph the line (••••••) corresponds to the index 0 I g ξ , the line (---) at the index 0 I β ξ and the line ( ______ ) at the index 0 I f ξ , p indicates each parameter g,β and f.

Figure 4 .
Figure 4. Dynamics of banana Moko with variable controls.

Figure 5 .
Figure 5. Behavior of infected plants of banana Moko disease in the time t.

Figure 6 .
Figure 6.The trajectories of the control u 1 (t).
Louadj K, Spiteri P, Demim F, et al.: Open Archive Toulouse Archive Ouverte OATAO is an open access repository that collects the work of Toulouse researchers and makes it freely available over the web where possible.This is an author ' s version published in : Official.2018.Nakakawa J, Mugisha JYT, Shaw MW, et al.: A mathematical model for the dynamics of banana Xanthomonas wilt with vertical transmission and inflorescence infection.J Biol Syst.2016; 24(1): 147-165.Publisher Full Text Rodrigues HS, Monteiro MTT, Torres DFM: Sensitivity Analysis in a Dengue Epidemiological Model.Hindawi Publ Corp Conference Papers in Science.2013.Publisher Full Text Viljoen A, Ssali C, Ndayihanzamaso G, et al.: Banana Pests and Diseases Field Guide for Disease Diagnostics and Data Collection Improvement of banana for smallholder farmers in the Great Lakes Region of Africa.2016.
Is the work clearly and accurately presented and does it cite the current literature?PartlyIs the study design appropriate and is the work technically sound?YesAre sufficient details of methods and analysis provided to allow replication by others?YesIf applicable, is the statistical analysis and its interpretation appropriate?Not applicableAre all the source data underlying the results available to ensure full reproducibility?Yes Are the conclusions drawn adequately supported by the results?Partly Competing Interests: No competing interests were disclosed.