Next Article in Journal
Diverse Physiological and Physical Responses among Wild, Landrace and Elite Barley Varieties Point to Novel Breeding Opportunities
Previous Article in Journal
De Novo Transcriptome Sequencing in Kiwifruit (Actinidia chinensis var. deliciosa (A Chev) Liang et Ferguson) and Development of Tissue-Specific Transcriptomic Resources
Previous Article in Special Issue
A Soft Clustering Approach to Detect Socio-Ecological Landscape Boundaries Using Bayesian Networks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Conservation Strategy for Palm Groves: Optimal Chemical Control Model for Red Palm Weevil, Rhynchophorus ferrugineus

1
Center for Agribusiness Biotechnology Research, University of Almería, La Cañada de S. Urbano, s/n, ES-04120 Almería, Spain
2
Department of Mathematics, University of Almería, La Cañada de San Urbano, ES-04120 Almería, Spain
3
MTA-ELTE Theoretical Biology and Evolutionary Ecology Research Group, Department of Plant Systematics, Ecology and Theoretical Biology, Centre for Ecological Research, Institute of Evolution, L. Eötvös University, H-1117 Budapest, Hungary
4
Department of Mathematics and Modelling, Institute of Mathematics and Basic Science, Hungarian University of Agriculture and Life Sciences, H-2100 Gödöllő, Hungary
*
Author to whom correspondence should be addressed.
Submission received: 25 March 2021 / Revised: 30 April 2021 / Accepted: 4 May 2021 / Published: 7 May 2021

Abstract

:
Rhynchophorus ferrugineus (Olivier, 1790) is an invasive pest species that constitutes one of the most important problems around the Mediterranean region and has been responsible for the loss of over 100,000 palm trees with an estimated annual cost of EUR several hundred million since its introduction into Europe. Methodological approaches of conservation ecology, such as multidisciplinary modelling, also apply in the management of cultural landscapes concerning ornamental plants, such as palm trees of the area. In this paper, we propose a dynamic model for the control of the red palm weevil, contributing in this way to the sustainability of an existing cultural landscape. The primary data set collected is a sample from the density-time function of a two-cohort pest population. This data set suggests a bimodal analytic description. If, from this data set, we calculate a sample from the accumulated density-time function (the integral of the density-time function), it displays a double sigmoid function (with two inflections). A good candidate for the analytical description of the latter is the sum of two logistic functions. As for the dynamic description of the process, a two-dimensional system of differential equations can be obtained, where the solution’s second component provides the analytical description of the original density-time function for the two-cohort population. Since the two-cohort waves appear in all three cycle stages, this reasoning applies to the subpopulations of larvae, pupae and adults. The model fitting is always performed using the SimFit package. On this basis, a mathematical model is proposed, which is sufficiently versatile to be of help in the control of this pest species in other geographical areas.

1. Introduction

The red palm weevil, Rhynchophorus ferrugineus, hereinafter referred to as RPW (Olivier, 1790) (Col.: Dryophthoridae), is a species native to Asia and Polynesia, and is characterized as an invasive pest of great severity and economic importance worldwide [1,2]. This species has expanded from its area of origin to different continents. It has been recorded in 28 countries in Asia, six in Africa, one in North America, two in Central America and the Caribbean, 14 in Europe and five in Oceania, causing serious problems in coconut trees, both in crop cultivation and ornamental use [1,3].
In Spain, RPW has spread widely over the Mediterranean coast, the interior of the peninsula and the islands, causing significant ecological and economic damage in areas such as the “Palmeral de Elche” (in the Valencian Community), declared a World Heritage Site by UNESCO, and in the Phoenix canariensis forest on the Canary Islands, which is a natural heritage site and represents the most relevant source of genetic variation for this palm tree [1,4,5].
Many aspects of RPW biology and ecology have been widely studied and are reported in several publications (e.g., [1,6,7,8,9]). However, the cold survival strategy adopted by this species, native to warm geographic areas, remains poorly known [10].
RPW infestations threaten agriculture (e.g., date plantations), affect the aesthetics of infested palm landscapes and uncontrolled urban populations (e.g., UNESCO World Heritage palm groves) and threaten natural areas (e.g., palm oases and endemic palm) [1,8,9,11,12].
In Spain during the 2005–2009 period, between 47,000 and 50,000 infested palm trees were destroyed; this involved containment/eradication costs of EUR 44.5 million, particularly in the Valencian Community where professional and family managed nurseries were affected [11,13,14]. In this regard, the FCEC [11] considered that infested (but not destroyed) palm trees might result in losses of between EUR 96 and 288 million.
Insecticide applications are, at present, the most effective method for protecting palms from attack by palm weevils [8]. Larval control, mainly by injection of insecticide into the trunk, carried out during the early stages of infestation, sometimes accompanied by application of spray for adult control, is currently the most employed (e.g., [3,8,15,16,17,18,19]). This is usually accompanied by monitoring and mass trapping of adults using semiochemicals (e.g., [3,9,20]). It should be noted that other control methods have been studied and developed, e.g., biological control (with entomophages and entomopathogens) and genetic methods (sterile backcrosses, sterile insect technique and chemosterilization), including two or more of the above control measures in integrated pest management programmes (e.g., [3,8,20]). The severity of R. ferrugineus and the problems associated with its control, including the difficulty of detecting infested palm trees early, the risk of mass captures in non-infested areas and the loss of control strategy effectiveness [9], motivate the search for new tools that increase the bioecological knowledge for dealing with RPW populations. In this regard, mathematical models can offer a holistic approach to understanding the biological behaviour of species, and they are an instrument for predicting and/or simulating complex systems that include the population dynamics of the pest and its control strategies [21,22,23].
In the plant protection field, mathematical models have focused on using sigmoid functions, such as Verhulst–Pearl or logistics, Gompertz, Weibull, Richards and monomolecular, as well as those models that combine them [24,25,26,27]. These functions have also been integrated to optimize the use of the resources and time spent on combating pests, through the formulation of optimal control problems aimed at keeping a pest population below the economic damage level as well as reducing the need to implement control strategies [21,22,23,28,29].
The existing literature on mathematical models related to the ecological aspects of R. ferrugineus is scarce, with only the work of Rossini et al. [30], who proposed a linear model for RPW development as a function of temperature. However, Hansen et al. [31] adapted a logistic model for another Curculionidae species, Dendroctonus rufipennis (Kirby), based on temperature as a diapause prediction factor of the insect’s fourth larval stage. It should be noted that the larval diapause in Curculionidae has also been recorded for Curculio elephas Gyllenhal, C. sikkimensis (Heller) and Smicronyx fulvus LeConte, where the larvae were subjected to temperatures below 15 °C [32,33,34,35].
In the R. ferrugineus larvae, Martín and Cabello [1] recorded a slowdown in the development of their last stages when the temperature dropped and remained in the 10–15 °C range for a period of more than 80 days, arguing that this phenomenon was due to a larval diapause, because the delay in the larval development and metabolic activity was related to changes in temperature [36,37].
The importance of RPW within palm production systems and/or ecosystems, and the need to deepen our understanding of the population behaviour of its stages, is fundamental to the development of mathematical models that can describe the population dynamics of the larvae, pupae and adults of R. ferrugineus, and to propose an optimal solution to control their larvae. This was the objective of this work.
The resulting modelling–methodological development also contributes to the sustainable management of the involved cultural landscape of palm groves and ornamental palm trees.
The objectives of this work were to answer two initial issues:
(1)
Is it possible to fix a mathematical model for the biology of RPW in conditions of cold winter climates?
(2)
Once the above model has been fixed, is it possible to use an optimal control model to simulate the optimum time for application of chemical larval control?

2. Materials and Methods

2.1. Biological Data Used

The trial was carried out in the “Cortijo del Olivar”, Municipal District of Almuñecar, Granada, Spain, 50 m above sea level (13.748219 N, 3.690619 W). It was in a room, communicated with the outside, located inside a country building; it was carried out between July and September of the following year. The trial started with a total of 70 RPW eggs, less than 24 h old, which were obtained from the laboratory rearing populations (Entomology Laboratory, University of Almeria, Almeria, Spain). The food substrate for larvae and adults was sugar cane (cultivar: “Algarrobeña”; traditionally grown in southern Spain, Almuñecar). Specimen handling was carried out following the methodology described in [38,39]. After development, 13 pairs were used for the evaluation of adult parameters. Data recorded were number of days for egg, larval and pupal stages, as well as mortality at these stages; adult longevity and fecundity were also recorded. Temperature was recorded daily by thermograph.

2.2. Mathematical Model

The primary data set collected is a sample from the density-time function of a two-cohort pest population. This data set suggests a bimodal analytic description.
Originally, the single logistic function is:
K 1 + K x 0 1 e r t
where x0 is the initial population; r: the logistic growth rate or steepness of the curve; K: the curve’s saturation value; in function (4) below, i subscripts (=1, 2) indicate the first or second cohort of the population.
Function (1) is the solution of the logistic differential equation:
x = r x 1 x K
with the initial condition   x 0 = x 0 . Function (1) is sigmoid; its derivative, with appropriate constants a 0 ,   a 1 ,   a 2 > 0 , is
y 1 t = a 0 e a 1 t 1 + a 2 e a 1 t 2 ,
displaying a unimodal curve. The sample from the density-time function of a two-cohort pest population suggests a bimodal analytic description, and for this, the sum of two functions of type (2) is a promising candidate [25]:
y 2 t = a 0 e a 1 t 1 + a 2 e a 1 t 2     + a 3 e a 4 t 1 + a 5 e a 4 t 2 ,
where a 0 ,   a 1 ,   a 2 ,   a 3 ,   a 4 , a 5 > 0 are the fitting parameters.
If we calculate a sample from the accumulated density-time function (the integral of the density-time function) from the data set of the density-time function, it shows a double sigmoid function (with two inflections). From the above, we conclude that a good candidate for the analytic description of the latter is the sum of two logistic functions (1):
K 1 1 + K 1 x 01 1 e r 1 t + K 2 1 + K 2 x 02 1 e r 2 t
Here, K i ,   r i , x 0 i > 0 (I = 1,2) are the fitting parameters.
It is worth emphasizing that, if we consider a two-cohort pest population of Rhynchophorus, the above reasoning can be applied to all the developmental stages: the larvae, pupae and adults.

2.3. Dynamic Model of Optimal Chemical Control

Considering any given stage, from functions (2) and (3), we easily obtain the following system of differential equations, where y 1 is the density-time function of a single cohort, and y 2 is the density-time function of the two-cohort population:
y 1 = a 1 y 1 1 2 1 + a 2 e a 1 t
y 2 = a 1 y 1 1 2 1 + a 2 e a 1 t + a 4 y 2 y 1 1 2 1 + a 5 e a 4 t
Although the dynamic model (5)–(6) can be adapted to any stage, for the chemical control of Rhynchophorus, the system (5)–(6) is used for a two-cohort larval population [1]. An operational time interval is fixed [0,T], and the set of admissible controls U ε 0 , T consisting of piecewise constant functions over a fixed uniform division of interval [0,T] with 0 u t ε   is defined. Constants μ 1 , μ 2 > 0   express the efficiency of the chemical on the respective populations. Then, from the (5)–(6) dynamics, we obtain the following control system:
y 1 = a 1 y 1 1 2 1 + a 2 e a 1 t μ 1 u   y 1
y 2 = a 1 y 1 1 2 1 + a 2 e a 1 t μ 1 u y 1 + a 4 y 2 y 1 1 2 1 + a 5 e a 4 t μ 2 u y 2
Now, on the basis of this control system, we set up the following optimal control problem: suppose that the damage caused by larval density y 2 in unit time is α y 2 ; then, the total damage during the time interval [0,T] is 0 T α y 2 t d t . If this damage is remains below a boundary, H, for the palm tree to survive, and we wish to minimize the cost of the chemical control proportional to 0 T u t d t , we obtain the following constrained optimal control problem:
Ψ u = 0 T u t d t m i n , u U ε 0 , T Equations   5 6   hold 0 T α y 2 t d t H y 1 0 ,   y 1 0     are   given.

2.4. Software Used

The above-mentioned mathematical models were fitted to the data using the TableCurve 2d v 5.0 [40] and SimFit v 6.9.8 [41] software packages.

3. Results

The results found were used to carry out the adjustment to the proposed mathematical model that explains the population dynamic of the insect’s life-cycle stages (Section 3.1). Subsequently, with this model, a control variable was used that would allow us to set the time for the chemical control of the pest in order to obtain an optimal result.

3.1. Results of Model Fitting

The results show that a portion of the larval population of R. ferrugineus developed normally (relative maximum in tp1), while another portion of the population lengthened its development up to 370 days (Table 1 and Figure 1). In the pupal and adult stages, two relative maxima (tp1 and tp3) are shown with a distance between them of approximately 160 days.
We note that the R2 value indicates that the goodness of fit is better (higher) for the double logistic function fitted to the accumulated data (0.99951) than for the fitting dynamics (5)–(6) (0.90319). Similar results were obtained for the rest of the analysed data. This observation justifies the use of a sigmoid function fitted to the accumulated data, smoother than the original density data, which displayed more random oscillations.
Subsequently, double logistic curves were also fitted to the accumulated density data (see Table 2 and Figure 2). In the latter, we also indicate the time moments tk, where the minima and maxima of the densities are attained, according to the inflexions of the accumulated densities.

3.2. Solution of the Optimal Control Problem

For a palm tree to survive, our objective is to keep the effect of the larval population under a boundary, H, at minimum cost. Using the Table 1 fitting parameters, we solved the optimal control problem (9) with H = 2000 units of biomass, T = 450, μ 1 = 0.95 , μ 2 = 0.5 ,   α = 0.7 ,   ε = 0.1 ,   y 1 0 = 13 ,   y 2 ( 0 ) = 14.
We solved the optimal control problem applying the toolbox developed for MatLab by Banga et al. [42], the results of which are shown in Figure 3.

4. Discussion

Aplasticity is well understood when it takes the form of responding to predictive environmental signals received during the individual’s lifespan. However, the so-called mixture of intragenotypic strategies has been neglected, perhaps due to the general assumption that it must be non-adaptive [43]. In this regard, several species of insects, especially from the Coleoptera order, may present two types of larval diapause (simple or prolonged), leading to diverse cohorts that differ in their life cycles, especially in terms of duration. The causes of what Ushatinskyaya [44] called “superdiapause” are poorly studied and not understood in detail. In our case, it seems clear from the results that RPW belongs to this group of insects.
The results for the first species, R. ferrugineus, showed that a portion of the larval population develop normally, but another portion could lengthen their development up to 370 days (Figure 1); this is demonstrated by the two relative maxima found in the pupal and adult stages. The above statement might indicate that the diapause exists as a survival mechanism. This phenomenon was later verified in other trials [1]. The presence or lack of a diapause in this species is crucial for the early detection and control of the pest, as well as for border inspections, since detection is carried out by ultrasound—if the larvae are in diapause, they are not feeding, and, therefore, do not make a noise with their jaws that can be detected [45]. Similarly, the results can be used to develop more complex predictive models to track the dynamics and expansion of this pest species.
The results of the logistic model developed for the population dynamics of the RPW larvae show a delay of larval development that would relate to the diapause manifesting as a survival mechanism. Martín and Cabello [1] noted that the diapause of RPW larvae was stimulated by a decline in temperature, of values between 10 °C and 15 °C, during the larval growth stage. Additionally, they indicated that this phenomenon promoted the division of the larval population into two groups: the first formed before the temperature drop occurred, in which the larvae continued their normal development (such as that recorded in tp1; Figure 1a), while the second group delayed its development by 170 days (similar to that shown from tp2 onwards; Figure 1a), with the consequent delay in pupal and adult morphogenesis (as in the relative maxima, tp1 and tp3, in Figure 1b,c, respectively).
The larval development behaviour of R. ferrugineus may represents a survival strategy that offers the maximum development opportunity in a univoltine cycle, as one group of larvae is able to accelerate their development to promote adult formation during the spring, while the other can slow down its development, influenced by the low winter temperatures, allowing adults to form that emerge during the autumn. This behaviour was also recorded by Hansen et al. [31] for the curculionid D. rufipennis, whose larval diapause was mathematically modelled in a logistic function that considers the thermal threshold as a predicting factor [31]. Other Curculionidae may also have a larval diapause in response to the drop in temperature, as is the case with Curculio elephas Gyllenhal [32,33], C. sikkimensis (Heller) [34] and Smicronyx fulvus LeConte [35].
Whether or not the RPW larvae are in diapause is crucial for early detection and pest control because most RPW inspection systems include equipment capable of detecting ultrasound produced by the larvae during feeding [46,47]. Therefore, detecting no noise will not mean the absence of larvae inside the palm stems, since diapause larvae do not feed [45].
The mathematical models developed for the larval, pupal and adult stages of RPW constitute a valuable tool for bioecological and/or prediction studies at the population level of this pest in various Spanish zones where the ambient temperature varies over the year, and can even stay below 15 °C during the winter months. At this time, the larvae can slow their development as a survival mechanism, as recorded by Martín and Cabello [1] in the municipality of Almuñecar (Granada, Spain) from the month of November to March, in which the temperature, both of the air and inside the palm tree, remained between 10 °C and 15 °C, with fluctuations within the palm tree of around 10 °C for more than 80 days.
In this regard, mathematical models can be used in various parts of Spain where palm trees are grown, such as the Andalusian coastal zones of the Mediterranean and the coastline from the Ebro Delta to Gibraltar, the Valle bajo of the Guadalquivir and Guadiana, as well as in coastal areas of the Rias Bajas in Galicia, where the ambient temperature can range from 10 °C to 12.5 °C during the month of January [48], also having an impact on the internal temperature of the palm trees.
Moreover, the optimal control model for RPW larvae reveals that the application of a chemical control measure can significantly reduce the initial larval population peak over the first 50 days. The action of chemical control also produces a significant decline in the second larval maximum, which begins to form from about day 100 in the uncontrolled population.
In turn, the rise in the chemical control trajectory seen in Figure 3, which begins after day 100, corresponds to the need to increase the use of insecticide to keep the palm biomass threshold below the economic damage level, since the larval development continues in this period, having been delayed by the temperature decline experienced by the population without chemical control. Similarly, this increase in the control trajectory can be related to the formation of RPW pupae and/or adults, which require greater control.
The logistic models developed for the population dynamics of the R. ferrugineus stages, as well as the optimal control model proposed, are instruments that can be incorporated into integrated pest control programmes, with the aim of improving decision making and reducing the cost of managing the RPW population. Similarly, these tools can be included within more complex prediction models, allowing us to understand RPW dynamics and expansion.
In this regard, it should be noted that this mathematical model has several advantages:
(A)
The model is sufficiently simple and flexible to allow its extension to this pest species in other geographical areas.
(B)
It saves some of the costs of experimentation on the effectiveness of the chemical control of RPW, as simulations can be carried out before expensive field trials, providing valuable information in advance.
(C)
With the introduction of optimisation into the model, the model could be used to calculate the ideal timing of chemical pest control, depending on climatic conditions and/or stage of infestation, with potential savings in treatment costs.
In summary, the results from the logistic models, for the pest species as well as for the different aspects of their biology/ecology, suggest that they could be an excellent tool to use in integrated pest management (IPM) programmes. This may also corroborate the results found regarding the dynamics of crop diseases [24,25], as their use can be generalized to crop protection.

Author Contributions

Conceptualization, T.C., Y.S.-R. and J.G.; Methodology and formal analysis: T.C., Y.S.-R., Z.V., M.G., I.L. and J.G.; Software, M.G. and I.L.; Visualization, I.L. and T.C.; Writing—original draft: J.G., Z.V., Y.S.-R., I.L., M.G. and T.C. All authors have read and agreed to the published version of the manuscript.

Funding

The research was funded by the Fundación Carolina, Spain.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

I.L. also acknowledges the support from CDTIME (University of Almería).

Conflicts of Interest

The authors declare no conflict of interest. The sponsors had no role in the design, execution, interpretation or writing of the study.

References

  1. Martin, M.M.; Cabello, T. Biología y Ecología del Curculiónido rojo de la Palmera, Rhynchophorus ferrugineus (Col.: Dryophthoridae); Universidad de Almeria: Almeria, Spain, 2005; pp. 1–202. [Google Scholar]
  2. Cabello, T. Rhynchophorus ferrugineus: Biología, dispersión y modo de acción. Phytoma 2012, 235, 36–38. [Google Scholar]
  3. CABI. Invasive species compendium: Rhynchophorus ferrugineus. Available online: https://www.cabi.org/isc/datasheet/47472 (accessed on 25 April 2021).
  4. Ferry, M.; Gómez, S. El picudo rojo de la palmera datilera: Gravedad de la plaga en España y necesidad de un cambio radical y urgente de estrategia en la lucha. Phytoma 2007, 184, 1–4. [Google Scholar]
  5. Ferry, M.; Gomez, S. El picudo rojo de la palmera datilera. Phytoma España 2007, 185, 42–46. [Google Scholar]
  6. Rochat, D.; Dembilio, O.; Jaques, J.A.; Suma, P.; Pergola, A.; Hamidi, R.; Kontodimas, D.; Soroker, V. Rhynchophorus ferrugineus: Taxonomy, distribution, biology, and life cicle. In Handbook of Major Palm Pests: Biology and Management; Soroker, V., Colazza, S., Eds.; John Wiley & Sons: Chichester, UK, 2017; pp. 69–104. [Google Scholar]
  7. Peri, E.; Rochat, D.; Belusic, G.; Ilic, M.; Soroker, V.; Barkan, S.; Guarino, S.; Bue, P.L.; Colazza, S. Rhynchophorus ferrugineus: Behavior, ecology and communication. In Handbook of Major Palm Pests: Biology and Management; Soroker, V., Colazza, S., Eds.; John Wiley & Sons: Chichester, UK, 2017; pp. 105–130. [Google Scholar]
  8. Milosavljevic, I.; El-Shafie, H.A.F.; Faleiro, J.R.; Hoddle, C.D.; Lewis, M.; Hoddle, M.S. Palmageddon: The wasting of ornamental palms by invasive plam weevils, Rhynchophorus spp. J. Pest Sci. 2019, 92, 143–156. [Google Scholar] [CrossRef]
  9. Giblin-Davis, R.M.; RomenoFaleiro, J.; Jacas, J.A.; Peña, J.E.; Vidyasagar, P.S.P.V. Biology and management of the red palm weevil, Rhynchophorus ferrugineus. In Potential Invasive Pests of Agricultural Crops; Peña, J.E., Ed.; CABI International: Walling-Ford, UK, 2013; pp. 1–34. [Google Scholar]
  10. Leon-Quinto, T.; Fimia, A.; Madrigal, R.; Serna, A. Morphological response of the red palm weevil, Rhynchophorus ferrugineus, to a transient low temperature analyzed by computer tomography and holographic microscopy. J. Therm. Biol. 2020, 94. [Google Scholar] [CrossRef]
  11. FCEC (Food Chain Evaluation Consortium). Quantification of Costs and Benefits of Amendments to the EU Plant Health Regime: Final Report. Available online: https://ec.europa.eu/food/sites/food/files/plant/docs/ph_biosec_rules_fcec_final_report_economic_study_plant_health_en.pdf (accessed on 9 November 2020).
  12. MacLeod, A.; Hussein, M. Economic and social impact of Rhynchophorus ferrugineus and Paysandisia archon on palms. In Handbook of Major Palm Pests: Biology and Management; Soroker, V., Colazza, S., Eds.; John Wiley & Sons: Chichester, UK, 2017; pp. 54–68. [Google Scholar]
  13. Cobos-Suárez, J.M. Situation of R. ferrugineus in Spain. In Proceedings of the Red Palm Weevil Control Strategy for Europe: International Conference, Valencia, Spain, 5–6 May 2010. [Google Scholar]
  14. Jacas, J.A.; Dembilio, O.; Llácer, E. Research activities focused on management of red palm weevil at the UJI-IVIA Associated Unit (Region of Valencia, Spain). EPPO Bull. 2011, 41, 122–127. [Google Scholar] [CrossRef]
  15. Rao, P.V.S.; Subramaniam, T.R.; Abraham, E.V. Control of the red palm weevil on coconut. J. Plant. Crop. 1973, 1, 26–27. [Google Scholar]
  16. El-Ezaby, F.A. Injection as a method to control the Red Indian date palm weevil Rhynchophorus ferrugineus. Arab J. Plant Prot. 1997, 15, 31–38. [Google Scholar]
  17. Cabello, T.; De La Peña, J.; Barranco, P. Laboratory test of imidacloprid and oxamyl against Rhynchophorus ferrugineus, new pest of palms in Spain. Ann. Appl. Biol. 1997, 130, 2–3. [Google Scholar]
  18. Barranco, P.; de la Pena, J.; Martin, M.M.; Cabello, T. Efficiency of chemical control of the new palm pest Rhynchophorus ferrugineus (Col.: Curculionidae). Bol. Sanid. Veg. Plagas 1998, 24, 301–306. [Google Scholar]
  19. Hernandez-Marante, D.; Folk, F.; Sanchez, A.; Fernandez-Escobar, R. Control del curculiónido ferruginoso de las palmeras (Rhynchophorus ferrugineus) mediante inyecciones al tronco y pulverización foliar. Bol. Sanid. Veg. Plagas 2003, 29, 563–573. [Google Scholar]
  20. Faleiro, J.R. A review of the issues and management of the red palm weevil Rhynchophorus ferrugineus (Col.: Rhynchophoridae) in coconut and date palm during the las one hundred years. Int. J. Trop. Insect Sci. 2006, 26, 135–154. [Google Scholar]
  21. Rafikov, M.; Balthazar, J.M. Optimal pest control problem in population dynamics. Comput. Appl. Math. 2005, 24, 65–81. [Google Scholar]
  22. Meng, X.; Song, Z.; Chen, L. A new mathematical model for optimal control strategies of integrated pest management. J. Biol. Syst. 2007, 15, 219–234. [Google Scholar] [CrossRef]
  23. Dabbs, K. Optimal Control in Discrete Pest Control Models. University of Tennessee Honors Thesis Projects. 2010. Available online: https://trace.tennessee.edu/utk_chanhonoproj/1375 (accessed on 6 May 2021).
  24. Amorim, L.; Bergamin-Filho, A.; Hau, B. Analysis of progress curves of sugarcane smut on different cultivars using functions of double sigmoid pattern. Phytopathology 1993, 83, 933–936. [Google Scholar] [CrossRef]
  25. Hau, B.; Amorim, L.; Bergamin-Filho, A. Mathematical functions to describe disease progress curves of double sigmoid pattern. Phytopathology 1993, 83, 928–932. [Google Scholar] [CrossRef]
  26. Carreño, R. 1996 Modelos logísticos. Aplicaciones a la Agronomía. Facultad de Ciencias Experimentales. Ph.D. Dissertation, Universidad de Almeria, Almeria, Spain, 1997; p. 168. [Google Scholar]
  27. Cabello, T.; Carreño, R. Métodos logísticos aplicados a la fenología de Noctuidos Plagas en el Sur de España (Lep.: Noctuidae). Bol. Sanid. Veg. Plagas 2002, 28, 226–319. [Google Scholar]
  28. Rafikov, M.; Angelelli, T. Optimization of biological pest control of sugarcane borer. In Proceedings of the 18th IEEE International Conference on Control Applications Part of 2009 IEEE Multi-Conference on Systems and Control, Saint Petersburg, Russia, 8–10 July 2009. [Google Scholar]
  29. Gallego, J.R.; López, I.; Gámez, M.; Cabello, T.; Varga, Z.; Garay, J. Simulation model applied to biological pest control by entomophagous species in commercial tomato greenhouses. Hung. Agric. Eng. 2013, 25, 67–70. [Google Scholar]
  30. Rossini, L.; Severini, M.; Contarini, M.; Speranza, S. A proposal of application in urban context to forecast pest insects’ life cycle. In Environmental and Territorial Modelling for Planning and Design; Leone, A., Gargiulo, C., Eds.; Federico II Open Access University Press: Napoli, Italy, 2018; pp. 169–178. [Google Scholar] [CrossRef]
  31. Hansen, E.M.; Bentz, B.J.; Powell, J.A.; Gray, D.R.; Vandygriff, J.C. Prepupal diapause and instar IV developmental rates of the spruce beetle, Dendroctonus rufipennis (Col.: Curculionidae). J. Insect Physiol. 2011, 57, 1347–1357. [Google Scholar] [CrossRef]
  32. Menu, F. Diapause development in the chestnut weevil Curculio elephas. Entomol. Exp. Appl. 1993, 69, 91–96. [Google Scholar] [CrossRef]
  33. Menu, F.; Debouzie, D. Larval development variation and adult emergence in the chestnut weevil Curculio elephans (Col., Curculionidae). J. Appl. Entomol. 1995, 119, 279–284. [Google Scholar] [CrossRef]
  34. Higaki, M. Effect of temperature on the termination of prolonged larval diapause in the chestnut weevil Curculio sikkimensis (Col.: Curculionidae). J. Insect Physiol. 2005, 51, 1352–1358. [Google Scholar] [CrossRef] [PubMed]
  35. Prasifka, J.R.; Rinehart, J.P.; Yocum, G.D. Nonconstant thermal regimes enhance overwintering success and accelerate diapause development for Smicronyx fulvus (Col.: Curculionidae). J. Econ. Entomol. 2015, 108, 1804–1809. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Danks, H.V. Insect Dormancy: An Ecological Perspective; Biological Survey of Canada: Monograph Series No. 1; National Museum of Natural Science: Ottawa, ON, Canada, 1987; pp. 1–439. [Google Scholar]
  37. Gordh, G.; Headrick, D.H. A Dictionary of Entomology; CABI Publishing: Wallingford, UK, 2001; pp. 1–1032. [Google Scholar]
  38. Rahalkar, G.W.; Harwalkar, M.R.; Rananavare, H.D. Development of red palm weevil, Rhynchophorus ferrugineus, on sugarcane. Indian J. Entomol. 1972, 34, 213–215. [Google Scholar]
  39. Rahalkar, G.W.; Harwalkar, M.R.; Rannavaer, H.D.; Tamhankar, A.J.; Shanthram, K. Rhynchophorus ferrugineus. In Handbook of Insect Rearing; Singh, P., Moore, R.F., Eds.; Elsevier: Amsterdam, The Netherlands, 1985; pp. 279–286. [Google Scholar]
  40. Jandel Scientific. Table Curve 2D-User´s Manual. Version 5.0; Jandel Scientific: San Rafael, CA, USA, 1994. [Google Scholar]
  41. Bardsley, W.G. SIMFIT: Simulation, Fitting, Statistics and Plotting. Reference Manual; University of Manchester: Manchester, UK, 2010. [Google Scholar]
  42. Banga, J.R.; Balsa-Canto, E.; Moles, C.G.; Alonso, A.A. Dynamic optimization of bioprocesses: Efficient and robust numerical strategies. J. Biotechnol. 2005, 117, 407–419. [Google Scholar] [CrossRef]
  43. Menu, F.; Debouzie, D. Coin-flipping plasticity and prolonged diapause in insects: Example of the chestnut weevil Curculio elephas (Col.: Curculionidae). Oecologia 1993, 93, 367–373. [Google Scholar] [CrossRef] [PubMed]
  44. Ushatinskaya, R.S. A critical review of the superdiapause in insects. Ann. Zool. 1984, 21, 3–30. [Google Scholar]
  45. Cabello, T. Population biology and dynamics of the Red Palm Weeevil, Rhynchophorus ferrugineus (Col.: Dryophothoridae) in Spain. In I Jornada Internacional sobre el Picudo Rojo de las Palmeras; Fundación Agrolimed, Generalitat Valenciana: Valencia, Spain, 2006; pp. 19–34. [Google Scholar]
  46. Potamitis, I.; Ganchev, T.; Kontodimas, D. On automatic bioacoustic detection of pests: The cases of Rhynchophorus ferrugineus and Sitophilus oryzae. J. Econ. Entomol. 2009, 102, 1681–1690. [Google Scholar] [CrossRef]
  47. Gutiérrez, A.; Ruiz, V.; Moltó, E.; Tapia, G.; Téllez, M.M. Development of a bioacoustic sensor for the early detection of red palm weevil (Rhynchophorus ferrugineus Oliver). Crop Prot. 2010, 29, 671–676. [Google Scholar] [CrossRef]
  48. Centro Nacional de Información Geográfica. Instituto Geográfico Nacional. Atlas Nacional de España. Available online: http://centrodedescargas.cnig.es/CentroDescargas/busquedaRedirigida.do?ruta=PUBLICACION_CNIG_DATOS_VARIOS/aneTematico/Espana_Temperatura-media-de-las-maximas-de-enero_1981–2010_mapa_14670_spa.pdf# (accessed on 11 November 2020).
Figure 1. Fitted curve of density-time function model (3) for larvae (a), pupae (b) and adults (c) of Rhynchophorus ferrugineus.
Figure 1. Fitted curve of density-time function model (3) for larvae (a), pupae (b) and adults (c) of Rhynchophorus ferrugineus.
Agronomy 11 00920 g001aAgronomy 11 00920 g001b
Figure 2. Mathematical models (double logistic functions) fitted to the accumulated data of the palm pest species, Rhynchophorus ferrugineus (a) larvae, (b) pupae and (c) adults stages, also displaying the corresponding bimodal density time functions.
Figure 2. Mathematical models (double logistic functions) fitted to the accumulated data of the palm pest species, Rhynchophorus ferrugineus (a) larvae, (b) pupae and (c) adults stages, also displaying the corresponding bimodal density time functions.
Agronomy 11 00920 g002aAgronomy 11 00920 g002b
Figure 3. (a) Trajectories for larva stage y_2 without control in model (5)–(6) and with optimal control. (b) Optimal control to be applied and solution of problem (7)–(8).
Figure 3. (a) Trajectories for larva stage y_2 without control in model (5)–(6) and with optimal control. (b) Optimal control to be applied and solution of problem (7)–(8).
Agronomy 11 00920 g003
Table 1. Values of fitting parameters of model (3) for the all stages of pest population of Rhynchophorus.
Table 1. Values of fitting parameters of model (3) for the all stages of pest population of Rhynchophorus.
StageFitting Parameters (*)Statistical Parmeters
a0a1a2a3a4a5d.f.R2 adj.p
Larvae673.7390.044103.566912117.260.0196780.0043380.97654<0.01
Pupa74980.40.05752792.26510500.90.02003115.495400.97744<0.01
Adult28372.70.04231485.7132187190.025462901.76400.95999<0.01
( * )   a 0 = K 1 X 01 1 K 1 r 1 ;    a 1 = r 1 ;    a 2 = K 1 X 01 1
a 3 = K 2 X 02 1 K 2 r 2 ;    a 4 = r 2 ;    a 5 = K 2 X 02 1
Table 2. Values of fitting parameters of double logistic functions for pest populations of Rhynchophorus.
Table 2. Values of fitting parameters of double logistic functions for pest populations of Rhynchophorus.
StageTerm (i)Fitting ParametersStatistical Parameters
Kix0irid.f.R2 adj.p
Larvae12951.43214.530.06771370.99764<0.01
21673.3957.61150.01661
Pupa11810.514.80240.05484400.99872<0.01
24621.5847.10670.01889
Adult11342.241.24920.04789430.99951<0.01
23031.581.324820.02461
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Solano-Rojas, Y.; Gámez, M.; López, I.; Garay, J.; Varga, Z.; Cabello, T. Conservation Strategy for Palm Groves: Optimal Chemical Control Model for Red Palm Weevil, Rhynchophorus ferrugineus. Agronomy 2021, 11, 920. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy11050920

AMA Style

Solano-Rojas Y, Gámez M, López I, Garay J, Varga Z, Cabello T. Conservation Strategy for Palm Groves: Optimal Chemical Control Model for Red Palm Weevil, Rhynchophorus ferrugineus. Agronomy. 2021; 11(5):920. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy11050920

Chicago/Turabian Style

Solano-Rojas, Y., M. Gámez, I. López, J. Garay, Z. Varga, and T. Cabello. 2021. "Conservation Strategy for Palm Groves: Optimal Chemical Control Model for Red Palm Weevil, Rhynchophorus ferrugineus" Agronomy 11, no. 5: 920. https://0-doi-org.brum.beds.ac.uk/10.3390/agronomy11050920

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