Next Article in Journal
Resistance Inducers for the Protection of Pedunculate Oak (Quercus robur L.) Seedlings against Powdery Mildew Erysiphe alphitoides
Next Article in Special Issue
Diagnosis and Mobile Application of Apple Leaf Disease Degree Based on a Small-Sample Dataset
Previous Article in Journal
Terpenoid Transport in Plants: How Far from the Final Picture?
Previous Article in Special Issue
MIX-NET: Deep Learning-Based Point Cloud Processing Method for Segmentation and Occlusion Leaf Restoration of Seedlings
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Big Data and Machine Learning to Improve European Grapevine Moth (Lobesia botrana) Predictions

by
Joaquín Balduque-Gil
1,*,
Francisco J. Lacueva-Pérez
2,
Gorka Labata-Lezaun
2,
Rafael del-Hoyo-Alonso
2,
Sergio Ilarri
3,
Eva Sánchez-Hernández
4,*,
Pablo Martín-Ramos
4 and
Juan J. Barriuso-Vargas
1
1
Department of Agricultural Sciences and Natural Environment, AgriFood Institute of Aragon (IA2), University of Zaragoza, Avenida Miguel Servet 177, 50013 Zaragoza, Spain
2
Department of Big Data and Cognitive Systems, Instituto Tecnológico de Aragón, ITAINNOVA, María de Luna 7-8, 50018 Zaragoza, Spain
3
Departamento de Informática e Ingeniería de Sistemas, Instituto de Investigación en Ingeniería de Aragón (I3A), Universidad de Zaragoza, María de Luna 1, 50018 Zaragoza, Spain
4
Department of Agricultural and Forestry Engineering, ETSIIAA, University of Valladolid, Avenida de Madrid 44, 34004 Palencia, Spain
*
Authors to whom correspondence should be addressed.
Plants 2023, 12(3), 633; https://doi.org/10.3390/plants12030633
Submission received: 22 December 2022 / Revised: 24 January 2023 / Accepted: 29 January 2023 / Published: 1 February 2023
(This article belongs to the Collection Application of AI in Plants)

Abstract

:
Machine Learning (ML) techniques can be used to convert Big Data into valuable information for agri-environmental applications, such as predictive pest modeling. Lobesia botrana (Denis & Schiffermüller) 1775 (Lepidoptera: Tortricidae) is one of the main pests of grapevine, causing high productivity losses in some vineyards worldwide. This work focuses on the optimization of the Touzeau model, a classical correlation model between temperature and L. botrana development using data-driven models. Data collected from field observations were combined with 30 GB of registered weather data updated every 30 min to train the ML models and make predictions on this pest’s flights, as well as to assess the accuracy of both Touzeau and ML models. The results obtained highlight a much higher F1 score of the ML models in comparison with the Touzeau model. The best-performing model was an artificial neural network of four layers, which considered several variables together and not only the temperature, taking advantage of the ability of ML models to find relationships in nonlinear systems. Despite the room for improvement of artificial intelligence-based models, the process and results presented herein highlight the benefits of ML applied to agricultural pest management strategies.

1. Introduction

Agriculture faces the challenge of increasing its productivity while reducing its environmental impact [1,2]. Nonetheless, agricultural systems are affected by different problems, such as environmental conditions, soil characteristics, water availability, and pest infestation. The latter has traditionally been managed through the knowledge and experience of the farmers themselves [3,4].
Today, agriculture is undergoing a major transformation in the collection and use of large amounts of data to make smarter decisions [5]. The applicability of new trends such as the Internet of Things (IoT) and Information and Communication Technologies (ICT) to agricultural activities enables farmers to adopt an approach based on the analysis of data from their farms to improve decision-making efficiency and mitigate pest risks [3,6]. Recent developments regarding open access to data, coupled with the unprecedented growth in the volume of data, referred to as Big Data (BD), have led to a shift in focus toward methods aimed at effectively managing such data for use in agri-environmental research [7]. In this sense, both traditional statistical techniques and Machine Learning (ML) techniques can be used to convert collected BD into timely and valuable information [8].
The grapevine moth, Lobesia botrana (Denis & Schiffermüller) 1975 (Lepidoptera: Tortricidae) is one of the main pests of grapevines [9,10,11,12], responsible for high productivity losses in some vineyards worldwide. Figure 1 shows how this pest can be identified in the field at different stages of development. L. botrana presents a facultative diapause [13,14,15,16], which results in a variable number of flights and generations per year, depending on the temperature and photoperiod [15,17]. L. botrana is multivoltine; in Mediterranean latitudes, it has three generations [15].
L. botrana causes different types of damage (Figure 2). On the one hand, direct damage caused by the pest occurs during the first generation, destroying a certain number of inflorescences, which generally does not affect crop yield [18]. In the second and third generations, the larvae directly attack the berries, causing more significant harvest losses [19]. Moreover, regarding indirect damage, in the second generation, the attack on berries favors the presence of Botrytis cinerea Pers.: Fr (Botryotinia fuckeliana (de Bary) Whetz.) [10,13,16].
In terms of crop protection, it has been experimentally demonstrated that integrated pest management (IPM) [20] is more effective than classical methods based solely on biological or chemical control alone [21]. IPM is defined as a long-term control strategy that combines biological, cultural, and chemical methods to reduce pathogenic populations to tolerable levels so that pests do not reach an economic threshold of damage [22]. In this line, as occurs with other pests, L. botrana control strategies are based on monitoring populations and applying control measures (such as the sex pheromone-based mating disruption technique, mentioned later in this section, or phytosanitary treatments, when necessary) to avoid reaching the aforementioned level of economic damage [23].
It is worth pointing out that this type of monitoring-based population tracking has sensitive aspects, such as the significant amount of time needed to carry it out, as well as a high level of technical expertise [4]. In addition to this, an insufficiently high monitoring frequency may represent a significant loss of information between field inspections [24]. In this sense, the potential of predictive models of future pest behavior makes modeling a powerful resource to complement monitoring, as well as being particularly useful for decision support systems (DSS), which help farmers to make decisions on crop protection measures [25,26]. According to Nansen et al. [27], a large part of the production systems are expected to face an increase in pest pressure, and therefore in potential phytosanitary applications. Because of this, the assessment of the sustainability of viticulture under future conditions, both from an environmental and economic point of view, can be considered of utmost importance for winegrowing operations [28]. Consequently, the importance of the use of DSS lies in their viability to increase the sustainability of agricultural production, which can lead to a reduction of the effects of agrochemicals on the environment and their economic cost at the farm level [29,30].
Considering that agricultural crops are managed biological systems, the utilization of applied mathematics developed for biological systems can be of great use in pest control through modeling and simulation, as Plant et al. [31] explained at the time. Nowadays, the role that modeling can play in plant health continues to arouse interest among entomologists and environmental scientists [32]. Bearing in mind that the objective of an optimal application of phytosanitary products is not only to control a pest but to predict and mitigate its appearance [2], predictive models are considered essential for efficient crop protection methods [33,34,35], because of their usefulness in deciding if and when treatments are necessary.
Since insects are ectothermic organisms, it has been shown that there is a direct relationship between their life cycle and abiotic factors, such as climatic factors (viz. temperature, relative humidity, wind direction, rainfall, photoperiod, and insolation) [24]. In the case of L. botrana, temperature is the variable with the greatest influence on its development [9,12,13,14,15,17,32,36,37,38,39,40,41]. Hence, the Touzeau predictive model for the development of L. botrana [42], developed in the Toulouse region (France), is based on the calculation of thermal integrals through temperature accumulation above a set threshold value. Flight monitoring of this pest, which is carried out via pheromone trapping of adults, provides potentially useful data to verify Touzeau’s predictions [38].
In recent years, the sex-pheromone-based mating disruption technique has been successfully implemented for L. botrana control [43,44], with a significant reduction in the use of chemical insecticides and higher protection of non-target insects [45,46]. This method is based on the use of pheromones emitted by insects, which they use as chemical messengers. The best known and most widely used are the sex pheromones emitted, in this case, by the females of Lepidoptera to attract males and facilitate mating. Thus, the strategy is to saturate the air with the female’s sex pheromone by arranging diffusers, disrupting communication between the female and the male [47]. As a consequence, mating is hindered or delayed, leading to a drastic and gradual reduction in the oviposition of fertile eggs in later generations, causing a noticeable decline in the pest population [46,48]. This technique began to be implemented in 2012 in the province of Zaragoza (Aragón, Spain), the study area on which this work focuses. At the present time, in Aragon, L. botrana is controlled almost exclusively using mating disruption (although, if necessary, phytosanitary treatments are carried out on an ad hoc basis).
The research presented herein focuses on the analysis of a correlation model between weather conditions and L. botrana development, analogous to the aforementioned Touzeau predictive model. This analysis is centered on three protected designations of origin (PDO), viz. Cariñena, Campo de Borja, and Calatayud, located in the province of Zaragoza. Given that the fit of the Touzeau model is not sufficiently accurate in this area, its improvement would entail two main advantages: (i) A contribution to the preventive control of L. botrana in case the established sex pheromone-based mating disruption system becomes ineffective (at present, the mating disruption method is still effective, but in other pests, it has been observed that males learn to differentiate between the pheromone of the diffusers and that of the females, which is much more complex); and (ii) the possibility to tackle the adaptation of this pest to the modification of meteorological conditions in the area due to climate change, which has already begun to be observed by several authors [9,15,16,40,49,50]. From both points of view, the optimization of the model would be very useful to schedule an efficient treatment calendar. It is therefore proposed to improve the current predictive model through the use of Artificial Intelligence (AI)-based systems, particularly using ML models based on neural networks, whose highly hierarchical structure and great learning capacity allow for particularly good classifications and predictions [1].

2. Results

A total of 10,000 models were trained. For all of them, the projection of the data for the last two weeks was considered input data. During training, these models are able to capture knowledge and are expected to find non-linear relationships with more than one variable. The variables identified as most relevant by each model were extracted from their descriptions. The best-performing trained model was an Artificial Neural Network (ANN) of four layers whose composition is shown in Table 1.
Regarding the aforementioned variables, the best-performing model considered the following variables as most relevant:
  • The Touzeau index.
  • The Chilling index.
  • The rainfall (in this case, calculated using daily data or accumulating the half-hourly samples, considering the starting date of dormancy as the start date to accumulation).
  • The longitude.
The results of the F1 metric for the classic Touzeau models and the ML models were obtained after comparing the day of the year (DOY) when the flight peaks were observed in the field with the DOY predicted by the models. Table 2 shows the comparison of observed flight peak DOY with the DOY predicted by the Touzeau and ML models for the 2008−2011 seasons. These results highlighted, firstly, that the Touzeau model did not provide a good fit in the study area, yielding a 0.03 F1 score, which is very low. On the other hand, the ML models resulted in much higher F1 scores, of up to 0.63.

3. Discussion

3.1. On the Model Performance

The F1 score reaches its best value at 1, while 0 is the worst score [51]. It should be clarified that there is no specific threshold value for considering an F1 score as ‘good’, so the model that produces the highest F1 score is generally considered the best. In this work, the Touzeau model obtained a 0.03 F1 score, which is very low compared with that obtained by the ANNs (of up to 0.63), so ML models can be regarded as more accurate.
This F1 score value of 0.63 may be assessed in a similar way as Sepúlveda et al. [52] did in their research: They concluded that the F1 score of 0.6 that they obtained could be interpreted as a low false negative rate and a similar false positive rate, with good compensation between errors, confirming a good F1 score metric.
Another interpretation would be based on the approach proposed by Nieto et al. [53], who used F1 to compare the results of two classification models developed for monitoring phenology. The best model under the F1 metric was 0.52 times better than the worst one. Although the scenario is not the same as the one for our study (given that their models were similarly developed and here we are comparing a classical model with an elaborated ANN), the same reasoning would suggest that the obtained ML model was 20 times better than the Touzeau model.

3.2. Advantages of ML Models and Most Relevant Weather Variables Identified by the Best-Performing ANN

As noted above, the Touzeau model takes into account the daily temperature record as the only parameter, i.e., it does not take into account the rainfall, relative humidity, or any other parameter collected by weather stations. In this sense, the Touzeau model considers all conditions optimal for the development of L. botrana. However, according to Stellwaag [54], the highest population levels of this pest, and therefore the most serious economic damage, are systematically recorded in areas with conditions considered optimal for this insect (an average annual temperature of around 20 °C and a relative humidity of around 70%). From this perspective, it can be explained that at times of unsuitable climatic conditions (low humidity, extreme temperatures), the species suffers a slowdown in its physiological processes and even adult mortality. The ML models, on the other hand, consider all the parameters recorded by the weather stations (Table S1). In addition, they can consider the 48 daily inputs corresponding to the data recorded every 30 min. In this way, ML models calculate the contribution of each of the 48 parts of the day to determine the daily contribution to cumulative indices and values. Furthermore, BD technologies allow for improving the quality of information used to train AI models, balancing skewed datasets (e.g., in our case, the number of samples with flights is lower than the number not having), thus reducing the risk of overfitting [55]. In addition, once built, the way the ML models are trained allows them to be easily updated with new data to adapt them to potential new scenarios.
Following on from the ability of ML models to consider all the weather parameters discussed above, the best-performing ANN considered temperature (Touzeau and Chilling indexes), rainfall, and longitude as the most relevant variables. As mentioned above, many authors have addressed the influence of temperature on the development of L. botrana. Others found that both temperature and rainfall were the variables that most influenced this pest distribution [12]. In this same line, Comșa et al. [50] found that warmer and drier seasons were more favorable. Zhan et al. [56] also explained that high temperature and low humidity provide optimal conditions for L. botrana, while rain together with low-temperature conditions seems to reduce the mating frequency and, subsequently, egg production.
In addition, previous works considered the relationship between low temperatures, on which the chilling index is based [57], and L. botrana. Andreadis et al. [14] studied the cold hardiness of this pest. Other research works pointed out that very low temperature, experienced by larval instars, tends to avert pupal diapause, with a very prolonged development of larvae with low temperatures [17]. Zhan et al. [56], on their side, highlighted that low temperatures had a great influence on the occurrence of this pest.
Regarding longitude, some authors analyzed the correlation between this parameter and other Lepidoptera. Khaghaninia et al. [58] performed their study with Cydia pomonella (Linnaeus) 1758 (Lepidoptera: Tortricidae), belonging to the same order and family as L. botrana, and their results showed a low correlation. On the other hand, Zhou et al. [59] found a strong relationship between this variable and another insect of this same order, viz. Parocneria orienta (Chao) 1978 (Lepidoptera: Erebidae); in particular, these authors noted that the hazard centroid shifted significantly with respect to longitude.

3.3. Applicability of the Developed ML Models

The developed ML model can help winegrowers of the three PDOs make better decisions on crop protection measures against L. botrana, relying on more accurate model predictions to carry out treatments only if and when they are needed. In this regard, once the model predicts a risk in terms of a pest infestation, and in terms of the field implementation of measures based on that prediction, the cost of taking or not taking measures should be assessed in relation to the risk raised.
Although, at present, L. botrana population density remains below an economic threshold of damage in those areas where the mating disruption method using sex pheromones has been implemented (which, as noted above, is still effective), the optimization of the predictions could allow winegrowers to improve the preventive control of L. botrana in case the mating disruption method becomes ineffective. Furthermore, it should help farmers to deal with the adaptation of this pest to the modification of meteorological conditions in the area due to climate change.

3.4. Room for Improvement of the AI-Based Models and Future Work

Based on results from previous work [60,61], an ANN was selected for training, but other more complex models could be trained to see if they would be able to obtain better results.
Although many models in grapevine entomology refer to a single species, the future of modeling in pest management should move towards global frameworks, covering more aspects of the problem (e.g., economics, crop yields, etc.), where the insect pest becomes part of the whole. For example, it may be necessary to choose the right timing for insecticide spraying against two different pests that overlap in the same season, and if both pests are predicted by a model, it would be easier to choose the best active substance, which would make the decision more environmentally friendly and cheaper [62].
Despite the outlined room for improvement of AI-based models, the process and results presented in this work highlight the benefits of ML applied to plant health strategies and its potential contribution to increasing the sustainability of agricultural activity, through a lower environmental impact as well as a lower economic cost of production (e.g., as a result of determining the right time to initiate insecticide application to control a pest or, preferably, using IPM tactics such as sexual confusion).
As part of the GRAPEVINE project (hiGh peRformAnce comPuting sErvices for preVentIon and coNtrol of pEsts in fruit crops), in which the present study is framed, further research is planned primarily along two lines of work: One will be validating how the ML models perform with data collected after the implementation of the mating disruption method, i.e., from 2012 to 2022, after the mating disruption method implementation. The second main issue to tackle in the future is related to testing the performance of the models in other geographical areas beyond the study area considered in this work.

4. Materials and Methods

This section introduces the field sites and the study period. How the L. botrana flight monitoring was performed is also presented, as well as the weather data used. The section is completed with a description of Touzeau and data-driven models.

4.1. Monitoring Field Sites and Study Period

The vineyards used as control fields are located in three PDOs (Cariñena, Campo de Borja, and Calatayud) in the province of Zaragoza (Aragón, Northeastern Spain). The study area, shown in Figure 3, is characterized by a continental climate, with the typical wind of the region—called ‘cierzo’—as the most outstanding characteristic, being a cold and dry north-western wind, with a Mediterranean summer influence. The monitoring sites are commercially cultivated vineyards belonging to wineries and cooperatives of the PDOs. The activity of these producers in terms of plant health is coordinated through the RedFAra network [63]. This network has a database repository with historical records of in situ monitoring and pest incidence. Thus, L. botrana flight data from the monitoring sites were available for the study period.
In the study area, the sex-pheromone-based mating disruption technique began to be implemented for L. botrana in 2012. Therefore, the study period selected for this work extends from 2004 to 2011 to avoid confusing and biased information on pest behavior. A total of 172 monitoring sites were used for this work during the entire study period (2004−2011).
Figure 4 (left) shows the location of the monitoring sites. Data from the cadaster of the Spanish registry [65] and the Aragon Open Data portal [66] provided the geographical description of the monitoring points, also allowing, through their location, to link each point to the nearest weather station as well as to information related to their geopositioning [60].

4.2. Lobesia botrana Flight Monitoring

Data from L. botrana flight activity at the monitoring sites used in the study belong to the repository of the RedFAra database [63]. Delta traps were placed one meter above the ground, one per hectare, with a gummed bottom along with a synthetic sex pheromone (E/Z-7,9-dodecadienil acetate) to monitor the flight activity of L. botrana males during the study period (Figure 5). Traps were placed during the whole vegetative period of the crop, and captures were checked every week. Pheromones were renewed every 6 weeks, according to the effective duration indicated by the supplier (OpenNatur SL, Lleida, Spain). Data on the captures over time provided information on the time of occurrence and duration of each generation at the monitoring sites.
The L. botrana life cycle begins when the first emerging adults from pupae that remained in diapause during the winter are detected. When sustained captures are encountered in pheromone delta traps, the first flight is considered to have taken place. This flight is estimated to end when the number of captures trends to zero. Further capture records should follow the same behavior, although the timing of these records will vary with the increase in average temperatures as spring and summer advance [67]. Data from these captures indicate the maximum flight moments of each generation, called ‘flight peaks’. The information provided by flight monitoring was analyzed together with the meteorological data records collected by weather stations.

4.3. Weather Data

The Agroclimatic Information System for Irrigation (SiAR) network, operated by the Spanish Ministry of Agriculture, has weather stations throughout Spain, with 49 of these stations located in Aragon [68]. Weather data stored in the SiAR database were used in this work to analyze the modeling of L. botrana. These data correspond to the SiAR stations listed in Table S2, whose locations are shown in Figure 4 (right). Modeling was performed by managing, between raw and processed data, 30 GB of data updated every 30 min.
The weather data from the stations were combined with data from field observations about flight monitoring. This was performed, on the one hand, to train the ML models—and then to make the predictions—and, on the other hand, to assess the accuracy of both the Touzeau and the ML models. The selection of the most representative weather station for a given vineyard was based on the following rules:
  • First, valid stations were selected for each year: A station was considered valid when, for the year in question, it provided data for at least 90% of the days of the year and the number of days that had less than 90% of the total samples per day (i.e., 48 when the sampling frequency was 30 min) was less than 90%.
  • Once the set of stations was defined for a year, the field observations from a vineyard in that year were combined with the climate data from the nearest station.
The weather stations provided data during the study period on the temperature, relative humidity, precipitation, wind speed and direction, and radiation, with daily and semi-hourly (every 30 min) records (Table S1).
When a data gap was detected (e.g., station Z22 did not provide data from 18:30 to 20:30 on 25th September 2021), the missing data were replaced. This was achieved using data provided by the European Centre for Medium-Range Weather Forecasts (ECMWF) ERA5 [69], accessed via the Open-Meteo Historical API [70]. ERA5 is the fifth-generation ECMWF atmospheric reanalysis of the global climate and is produced by the Copernicus Climate Change Service (C3S) at ECMWF. It provides hourly data readings on a 30 km grid, which had to be transformed to fit the semi-hourly (every 30 min) SiAR records. For that, all variables were copied into two records, 00 min and 30 min, except for the accumulated precipitation, with which a single register was divided by 2.
Regarding the recorded temperatures, and to analyze the correlation between these and the development of L. botrana, the daily maximum (TM) and minimum (Tm) were adjusted to effective temperatures [71]. This was due to the fact that, when working with cumulative indexes, negative temperatures could alter the real value of these indexes. Therefore, to avoid this potential disruption, any temperature below 0 °C was adjusted to the effective temperature of 0 °C.

4.4. Touzeau Model

Classical models of L. botrana, such as the Touzeau model, are based on temperature accumulation (in °C), which is considered to be the most influential environmental factor on the seasonal behavior of this pest. Degree-day accumulation methods using thermal integrals are the most commonly used methods to analyze the relationship between insect development and temperature [72]. These methods are based on the amount of heat energy accumulated daily above a minimum temperature threshold called the base temperature (TB). The Touzeau model uses the Touzeau index to calculate the accumulated thermal integral for L. botrana, based on the daily maximum (TM) and minimum (Tm) temperatures, and sets TB = 10 °C [42]. In the present work, in addition to TB, a temperature limit of 30 °C (above which the temperature accumulation stops) was taken into account, in agreement with Gabel et al. [73], Del Tío et al. [74], and Gallardo et al. [75]. Therefore, all temperatures with values above 30 °C used to calculate the Touzeau index were adjusted to 30 °C.
The Touzeau model index function is mathematically expressed by:
T o u = DOY = 1 n T M + T m 2 T B ,
where DOY refers to the day of the year or day of the Julian calendar. Thus, the temperature accumulation for the Touzeau index starts on January 1st (DOY = 1).
Table S3 shows the different developmental stages of L. botrana for the three generations, together with the cumulative daily temperatures needed to reach them, expressed according to the Touzeau model index (Tou). Each of the stages has its calorific requirements. The thermal integral accumulated from a starting point using the Tou is used to predict when a stage will be reached. The monitoring methods followed in this paper focus on flights of adult generations. As shown in Table S3, the accumulated Tou for L. botrana generations were 125 for 1st G, 500 for 2nd G, and 950 for 3rd G. In this way, to assess the accuracy of the Touzeau model, the DOY when the maximum population of each generation takes place in the field was identified according to flight peaks. Then, this DOY was compared with the DOY when the accumulated Tou reached the values indicated by the model for each generation.

4.5. Data-Driven Models Used in This Study

Data-driven models are capable of extracting valuable insights from huge amounts of data without the need to define a law or equation to explain a given phenomenon [76]. To do so, they apply BD, analytics technologies and AI algorithms for the automation of data capturing and processing. They can combine heterogeneous multi-sourced data for discovering, understanding, and evaluating new models, which represent relationships that are too complex to be found by mathematical ones [77]. In our context, these computational technologies make it possible to use climatic variables, such as the maximum temperature, minimum temperature, humidity, rainfall, or wind speed, which are known factors influencing pests, but whose influence is not yet well understood [55].
Fenu et al. [55] and Domingues et al. [77] performed a review of the AI, ML, and Deep Learning (DL) algorithms used by different authors for creating pest models. The most popular models are those that use only climate data [55]. Several algorithms with different computational complexity have been used, from Logistic Regression to Convolutional Neuronal Networks (CNN), passing through Support Vector Machines (SVM), Support Vector Regression (SVR), Multi-Linear Regression (MLR), Random Forest (RF), and ANN.
Usually, the predicted variable is a binary variable representing the presence or absence of the pest. The performance of the models is measured using accuracy, precision, recall, or F1 metrics. However, due to the great heterogeneity of experimental conditions (i.e., approaches, datasets, parameters, and performance metrics), it is difficult—if not unreliable—to make a systematic comparison of the performance of approaches presented in different papers [55].
For our study, and based on the results presented in previous work [60], we created a model using an ANN. The aim of the model obtained was to predict the probability that a given date, within a 2-week horizon, is the day of flight peak for each of the generations of L. botrana in a given season.
To train (and use) the models, the following transformations were performed on the SiAR and RedFAra datasets:
  • For each monitoring site and season pair, the day of the flight peak for each generation was labeled. This was the variable to be predicted.
  • Growing degree days [78], Touzeau [42], and Chilling [57] indexes were calculated based on both the proportional accumulation of half-hourly samples and daily observations.
  • The aforementioned indexes were calculated considering different accumulation start dates, viz. January 1st, February 1st, or the starting date of dormancy, calculated as the first day of Autumn when the temperature remains under 10 °C. This date was considered to try to increase the accuracy of grapevine phenology predictions [79].
  • Registered radiation values.
  • Wind direction and speed were transformed into daily indexes, calculated using the following rules: wind direction was classified into one of eight categories (N, NE, E, SE, S, SW, W, or NW), and the average speed per day was calculated for each of those categories.
  • The raw weather data were also provided to build the models.
  • The weather data for the 14 days before a given date were projected horizontally so that the model could make predictions over a 2-week horizon.
The models were trained using Scikit-learn [80], Keras [81], and Optuna [82] on Python [83]. As noted above, based on results from previous work [60], an ANN was chosen for training. The selection of the data for testing and validation was performed considering data from different monitoring sites over different years during the study period. The ANN training process used the hyperparameters shown in Table 3.

4.6. Model Performance Assessment

The F1 score performance metric was used to assess the accuracy of both the Touzeau and ML models. Although it is slightly more difficult to calculate for the Touzeau model results than other metrics such as R2, it has been used by several authors to compare the results of classification models [53,84,85]. F1 is a widely-used metric for assessing the performance of binary classification problems and for comparing the results obtained by different models. It is a good metric when the performance in terms of true positive results is evaluated. However, F1 does not depend on true negatives and is sensitive to the distribution of the labels [86]. In addition, it does enable the user to consider the cost of failing in the classification.
The F1 metric was computed based on the differences between the predicted and observed values; more specifically, the DOY when the flight peaks were observed in the field was compared with the DOY predicted by the models. This yields four possibilities: True positives, true negatives, false positives, and false negatives, which were used to obtain the F1 metric, which is the harmonic mean of the precision and recall. The F1 was therefore calculated with the following formulas:
Precision = TP TP + FP ,
Recall = TP TP + FN ,
F 1 = 2 × Precision   ×   Recall Precision + Recall ,
The sites selected to validate both models from all available monitoring sites are presented in Tables S4 and S5. These sites used for the validations were randomly selected, seeking to obtain a representative sample free of bias in the choice. Validation was carried out by comparing predicted and observed values. The data corresponding to the study period from the monitoring sites listed in Table S4 was used for validation of the Touzeau model. In the same way, the monitoring sites appearing in Table S5 were used for validation of the best-performing ANN model of 4 layers of neurons.

5. Conclusions

The Touzeau model for L. botrana, based only on the temperature, showed a poor prediction performance in the area of study (three PDOs in Aragón, Northeastern Spain), with an F1 score of 0.03. In comparison, artificial neural network-type ML models, trained using data from 172 monitoring sites in the study area and considering several meteorological parameters (viz. temperature, relative humidity, wind speed, wind direction, solar radiation, and accumulated precipitation), reached an F1 score of 0.63. The best-performing model, a four-layer ANN, considered the following variables as most relevant: The Touzeau index, the Chilling index, rainfall (calculated using daily data or accumulating the half-hourly samples considering the starting date of dormancy as the start date to accumulation), and the longitude. The improvement of classical predictive models by using the reported ML models can complement pest-monitoring methods in integrated pest management, making decision-support systems more accurate.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/plants12030633/s1, Table S1: Parameters recorded daily every 30 min by the weather stations used in the modeling; Table S2: List of SiAR weather stations used for modeling; Table S3: Stages of development and sum of cumulative temperatures required according to the Touzeau model; Table S4: Monitoring sites used for validation of the Touzeau model; Table S5: Monitoring sites used for validation of the ML models.

Author Contributions

Conceptualization, J.B.-G., F.J.L.-P., R.d.-H.-A., S.I. and J.J.B.-V.; methodology, J.B.-G., F.J.L.-P. and J.J.B.-V.; validation, R.d.-H.-A. and J.J.B.-V.; formal analysis, J.B.-G., F.J.L.-P. and J.J.B.-V.; investigation, J.B.-G., F.J.L.-P., E.S.-H., P.M.-R., G.L.-L., R.d.-H.-A., S.I. and J.J.B.-V.; resources, R.d.-H.-A. and J.J.B.-V.; writing—original draft preparation, J.B.-G., F.J.L.-P., E.S.-H., P.M.-R., G.L.-L., R.d.-H.-A., S.I. and J.J.B.-V.; writing—review and editing, J.B.-G., F.J.L.-P., E.S.-H. and P.M.-R.; visualization, J.B.-G. and E.S.-H.; supervision, J.J.B.-V.; project administration, R.d.-H.-A.; funding acquisition, R.d.-H.-A. and J.J.B.-V. All authors have read and agreed to the published version of the manuscript.

Funding

This research was co-financed by the European Union’s Connecting Europe Facility (CEF) Grant No. INEA/CEF/ICT/A2018/1837816 GRAPEVINE project (hiGh peRformAnce comPuting sErvices for preVentIon and coNtrol of pEsts in fruit crops).

Data Availability Statement

The data presented in this study are available upon request from the corresponding author. The data are not publicly available due to their relevance to an ongoing Ph.D. thesis.

Acknowledgments

The authors thank the support of the project PID2020-113037RB-I00/AEI/10.13039/501100011033 and the funds of the Government of Aragon to the research groups COSMOS and IODIDE.

Conflicts of Interest

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

References

  1. Kamilaris, A.; Prenafeta-Boldú, F.X. Deep learning in agriculture: A survey. Comput. Electron. Agric. 2018, 147, 70–90. [Google Scholar] [CrossRef]
  2. Zangina, U.; Buyamin, S.; Aman, M.N.; Abidin, M.S.Z.; Mahmud, M.S.A. A greedy approach to improve pesticide application for precision agriculture using model predictive control. Comput. Electron. Agric. 2021, 182, 105984. [Google Scholar] [CrossRef]
  3. Ip, R.H.L.; Ang, L.-M.; Seng, K.P.; Broster, J.C.; Pratley, J.E. Big data and machine learning for crop protection. Comput. Electron. Agric. 2018, 151, 376–383. [Google Scholar] [CrossRef]
  4. Gutiérrez, S.; Hernández, I.; Ceballos, S.; Barrio, I.; Díez-Navajas, A.M.; Tardaguila, J. Deep learning for the differentiation of downy mildew and spider mite in grapevine under field conditions. Comput. Electron. Agric. 2021, 182, 105991. [Google Scholar] [CrossRef]
  5. Pham, X.; Stack, M. How data analytics is transforming agriculture. Bus. Horiz. 2018, 61, 125–133. [Google Scholar] [CrossRef]
  6. Kamilaris, A.; Kartakoullis, A.; Prenafeta-Boldú, F.X. A review on the practice of big data analysis in agriculture. Comput. Electron. Agric. 2017, 143, 23–37. [Google Scholar] [CrossRef]
  7. Lokers, R.; Knapen, R.; Janssen, S.; van Randen, Y.; Jansen, J. Analysis of Big Data technologies for use in agro-environmental science. Environ. Model. Softw. 2016, 84, 494–504. [Google Scholar] [CrossRef]
  8. Lokhorst, C.; de Mol, R.M.; Kamphuis, C. Invited review: Big Data in precision dairy farming. Animal 2019, 13, 1519–1528. [Google Scholar] [CrossRef]
  9. Gutierrez, A.P.; Ponti, L.; Cooper, M.L.; Gilioli, G.; Baumgärtner, J.; Duso, C. Prospective analysis of the invasive potential of the European grapevine moth Lobesia botrana (Den. & Schiff.) in California. Agric. For. Entomol. 2012, 14, 225–238. [Google Scholar] [CrossRef]
  10. Aguilera Sammaritano, J.; Deymié, M.; Herrera, M.; Vazquez, F.; Cuthbertson, A.G.S.; López-Lastra, C.; Lechner, B. The entomopathogenic fungus, Metarhizium anisopliae for the European grapevine moth, Lobesia botrana Den. & Schiff. (Lepidoptera: Tortricidae) and its effect to the phytopathogenic fungus, Botrytis cinerea. Egypt. J. Biol. Pest Control 2018, 28, 83. [Google Scholar] [CrossRef]
  11. Lucchi, A.; Sambado, P.; Juan Royo, A.B.; Bagnoli, B.; Conte, G.; Benelli, G. Disrupting mating of Lobesia botrana using sex pheromone aerosol devices. Environ. Sci. Pollut. Res. 2018, 25, 22196–22204. [Google Scholar] [CrossRef]
  12. Rank, A.; Ramos, R.S.; da Silva, R.S.; Soares, J.R.S.; Picanço, M.C.; Fidelis, E.G. Risk of the introduction of Lobesia botrana in suitable areas for Vitis vinifera. J. Pest Sci. 2020, 93, 1167–1179. [Google Scholar] [CrossRef]
  13. Tzanakakis, Μ.Ε.; Savopoulou-Soultani, Μ.; Oustapassidis, C.S.; Verras, S.C.; Hatziemmanouel, H. Induction of dormancy in Lobesia botrana by long day and high temperature conditions. Entomol. Hell. 1988, 6. [Google Scholar] [CrossRef]
  14. Andreadis, S.S.; Milonas, P.G.; Savopoulou-Soultani, M. Cold hardiness of diapausing and non-diapausing pupae of the European grapevine moth, Lobesia botrana. Entomol. Exp. Appl. 2005, 117, 113–118. [Google Scholar] [CrossRef]
  15. Martín-Vertedor, D.; Ferrero-García, J.J.; Torres-Vila, L.M. Global warming affects phenology and voltinism of Lobesia botrana in Spain. Agric. For. Entomol. 2010, 12, 169–176. [Google Scholar] [CrossRef]
  16. Caffarra, A.; Rinaldi, M.; Eccel, E.; Rossi, V.; Pertot, I. Modelling the impact of climate change on the interaction between grapevine and its pests and pathogens: European grapevine moth and powdery mildew. Agric. Ecosyst. Environ. 2012, 148, 89–101. [Google Scholar] [CrossRef]
  17. Roditakis, N.E.; Karandinos, M.G. Effects of photoperiod and temperature on pupal diapause induction of grape berry moth Lobesia botrana. Physiol. Entomol. 2001, 26, 329–340. [Google Scholar] [CrossRef]
  18. Moschos, T. Yield loss quantification and assessment of economic injury level for the anthophagous generation of the European grapevine moth Lobesia botrana Den. et Schiff. (Lepidoptera: Tortricidae). Int. J. Pest Manag. 2005, 51, 81–89. [Google Scholar] [CrossRef]
  19. Pavan, F.; Bigot, G.; Cargnus, E.; Zandigiacomo, P. Influence of the carpophagous generations of the European grapevine moth Lobesia botrana on grape bunch rots. Phytoparasitica 2013, 42, 61–69. [Google Scholar] [CrossRef]
  20. Stern, V.M.; Smith, R.F.; van den Bosch, R.; Hagen, K.S. The integration of chemical and biological control of the spotted alfalfa aphid: The integrated control concept. Hilgardia 1959, 29, 81–101. [Google Scholar] [CrossRef] [Green Version]
  21. van Lenteren, J.C. Success in biological control of arthropods by augmentation of natural enemies. In Biological Control: Measures of Success; Gurr, G., Wratten, S., Eds.; Springer: Dordrecht, The Netherlands, 2000; pp. 77–103. [Google Scholar] [CrossRef]
  22. Tang, S.; Cheke, R.A. Models for integrated pest control and their biological implications. Math. Biosci. 2008, 215, 115–125. [Google Scholar] [CrossRef] [PubMed]
  23. Sáenz-de-Cabezón, E.; Hernández, L.J.; Rivas, M.T.; García-Ruiz, E.; Marco, V.; Pérez-Moreno, I.; Sáenz-de-Cabezón, F.J. A computer implementation of the partition of the unity procedure and its application to arthropod population dynamics. A case study on the European grape berry moth. Math. Comput. Simul. 2011, 82, 2–14. [Google Scholar] [CrossRef]
  24. Rossini, L.; Bruzzone, O.A.; Contarini, M.; Bufacchi, L.; Speranza, S. A physiologically based ODE model for an old pest: Modeling life cycle and population dynamics of Bactrocera oleae (Rossi). Agronomy 2022, 12, 2298. [Google Scholar] [CrossRef]
  25. Knight, J.D.; Mumford, J.D. Decision support systems in crop protection. Outlook Agric. 2016, 23, 281–285. [Google Scholar] [CrossRef]
  26. Zhai, Z.; Martínez, J.F.; Beltran, V.; Martínez, N.L. Decision support systems for agriculture 4.0: Survey and challenges. Comput. Electron. Agric. 2020, 170, 105256. [Google Scholar] [CrossRef]
  27. Nansen, C.; Ferguson, J.C.; Moore, J.; Groves, L.; Emery, R.; Garel, N.; Hewitt, A. Optimizing pesticide spray coverage using a novel web and smartphone tool, SnapCard. Agron. Sustain. Dev. 2015, 35, 1075–1085. [Google Scholar] [CrossRef]
  28. Santos, J.A.; Fraga, H.; Malheiro, A.C.; Moutinho-Pereira, J.; Dinis, L.-T.; Correia, C.; Moriondo, M.; Leolini, L.; Dibari, C.; Costafreda-Aumedes, S.; et al. A review of the potential climate change impacts and adaptation options for European viticulture. Appl. Sci. 2020, 10, 3092. [Google Scholar] [CrossRef]
  29. Fabre, F.; Plantegenest, M.; Yuen, J. Financial benefit of using crop protection decision rules over systematic spraying strategies. Phytopathology 2007, 97, 1484–1490. [Google Scholar] [CrossRef]
  30. Michels, M.; Bonke, V.; Musshoff, O. Understanding the adoption of smartphone apps in crop protection. Precis. Agric. 2020, 21, 1209–1226. [Google Scholar] [CrossRef]
  31. Plant, R.E.; Mangel, M. Modeling and simulation in agricultural pest management. SIAM Rev. 1987, 29, 235–261. [Google Scholar] [CrossRef]
  32. Rossini, L.; Severini, M.; Contarini, M.; Speranza, S. EntoSim, a ROOT-based simulator to forecast insects’ life cycle: Description and application in the case of Lobesia botrana. Crop Prot. 2020, 129, 105024. [Google Scholar] [CrossRef]
  33. Moravie, M.A.; Davison, A.C.; Pasquier, D.; Charmillot, P.J. Bayesian forecasting of grape moth emergence. Ecol. Model. 2006, 197, 478–489. [Google Scholar] [CrossRef]
  34. Rossi, V.; Giosuè, S.; Caffi, T. Modelling plant diseases for decision making in crop protection. In Precision Crop Protection—The Challenge and Use of Heterogeneity; Oerke, E.-C., Gerhards, R., Menz, G., Sikora, R.A., Eds.; Springer: Dordrecht, The Netherlands, 2010; pp. 241–258. [Google Scholar] [CrossRef]
  35. Olatinwo, R.; Hoogenboom, G. Weather-based Pest Forecasting for Efficient Crop Protection. In Integrated Pest Management; Abrol, D.P., Ed.; Academic Press: San Diego, CA, USA, 2014; pp. 59–78. [Google Scholar] [CrossRef]
  36. Briere, J.-F.; Pracros, P. Comparison of temperature-dependent growth models with the development of Lobesia botrana (Lepidoptera: Tortricidae). Environ. Entomol. 1998, 27, 94–101. [Google Scholar] [CrossRef]
  37. Milonas, P.G.; Savopoulou-Soultani, M.; Stavridis, D.G. Day-degree models for predicting the generation time and flight activity of local populations of Lobesia botrana (Den. & Schiff.) (Lep., Tortricidae) in Greece. J. Appl. Entomol. 2001, 125, 515–518. [Google Scholar] [CrossRef]
  38. Amo-Salas, M.; Ortega-López, V.; Harman, R.; Alonso-González, A. A new model for predicting the flight activity of Lobesia botrana (Lepidoptera: Tortricidae). Crop Prot. 2011, 30, 1586–1593. [Google Scholar] [CrossRef]
  39. Moshtaghi Maleki, F.; Iranipour, S.; Hejazi, M.J.; Saber, M. Temperature-dependent age-specific demography of grapevine moth (Lobesia botrana) (Lepidoptera: Tortricidae): Jackknife vs. bootstrap techniques. Arch. Phytopathol. Plant Prot. 2016, 49, 263–280. [Google Scholar] [CrossRef]
  40. Gutierrez, A.P.; Ponti, L.; Gilioli, G.; Baumgärtner, J. Climate warming effects on grape and grapevine moth (Lobesia botrana) in the Palearctic region. Agric. For. Entomol. 2018, 20, 255–271. [Google Scholar] [CrossRef] [Green Version]
  41. Aguirre-Zapata, E.; Morales, H.; Dagatti, C.V.; di Sciascio, F.; Amicarelli, A.N. Semi physical growth model of Lobesia botrana under laboratory conditions for Argentina’s Cuyo region. Ecol. Model. 2022, 464, 109803. [Google Scholar] [CrossRef]
  42. Touzeau, J. Modélisation de l’évolution de l’Eudémis de la Vigne pour la région Midi Pyrénées. Boll. Zool. Agrar. Bachic. Ser. II 1981, 16, 26–28. [Google Scholar]
  43. Lucchi, A.; Suma, P.; Ladurner, E.; Iodice, A.; Savino, F.; Ricciardi, R.; Cosci, F.; Marchesini, E.; Conte, G.; Benelli, G. Managing the vine mealybug, Planococcus ficus, through pheromone-mediated mating disruption. Environ. Sci. Pollut. Res. 2019, 26, 10708–10718. [Google Scholar] [CrossRef]
  44. Ricciardi, R.; Benelli, G.; Suma, P.; Cosci, F.; Di Giovanni, F.; Zeni, V.; Conte, G.; Marchesini, E.; Savino, F.; Ladurner, E.; et al. One device for two pests: A new double dispenser for mating disruption of Lobesia botrana and Planococcus ficus. Entomol. Gen. 2022, 42, 289–307. [Google Scholar] [CrossRef]
  45. Welter, S.C.; Pickel, C.; Millar, J.; Cave, F.; Van Steenwyk, R.A.; Dunley, J. Pheromone mating disruption offers selective management options for key pests. Calif. Agric. 2005, 59, 16–22. [Google Scholar] [CrossRef]
  46. Benelli; Lucchi; Thomson; Ioriatti. Sex pheromone aerosol devices for mating disruption: Challenges for a brighter future. Insects 2019, 10, 308. [Google Scholar] [CrossRef] [PubMed]
  47. Witzgall, P.; Kirsch, P.; Cork, A. Sex pheromones and their impact on pest management. J. Chem. Ecol. 2010, 36, 80–100. [Google Scholar] [CrossRef] [PubMed]
  48. Kong, W.N.; Li, J.; Fan, R.J.; Li, S.C.; Ma, R.Y. Sex-pheromone-mediated mating disruption technology for the oriental fruit moth, Grapholita molesta (Busck) (Lepidoptera: Tortricidae): Overview and prospects. Psyche A J. Entomol. 2014, 2014, 253924. [Google Scholar] [CrossRef]
  49. Iltis, C.; Moreau, J.; Pecharová, K.; Thiéry, D.; Louâpre, P. Reproductive performance of the European grapevine moth Lobesia botrana (Tortricidae) is adversely affected by warming scenario. J. Pest Sci. 2020, 93, 679–689. [Google Scholar] [CrossRef]
  50. Comșa, M.; Tomoiagă, L.L.; Muntean, M.-D.; Ivan, M.M.; Orian, S.M.; Popescu, D.M.; Chedea, V.S. The effects of climate change on the activity of the Lobesia botrana and Eupoecilia ambiguella moths on the grapevine cultivars from the Târnave vineyard. Sustainability 2022, 14, 14554. [Google Scholar] [CrossRef]
  51. Ting, K.M. Confusion Matrix. In Encyclopedia of Machine Learning and Data Mining; Sammut, C., Webb, G.I., Eds.; Springer: New York, NY, USA, 2017; p. 260. [Google Scholar] [CrossRef]
  52. Sepúlveda, J.; Velastin, S.A. F1 score assessment of Gaussian mixture background subtraction algorithms using the MuHAVi dataset. In In Proceedings of the 6th International Conference on Imaging for Crime Prevention and Detection (ICDP-15), London, UK, 15–17 July 2015. [Google Scholar]
  53. Nieto, L.; Schwalbert, R.; Prasad, P.V.V.; Olson, B.J.S.C.; Ciampitti, I.A. An integrated approach of field, weather, and satellite data for monitoring maize phenology. Sci. Rep. 2021, 11, 15711. [Google Scholar] [CrossRef]
  54. Stellwaag, F. Die Weinbauinsekten der Kulturländer: Lehr- und Handbuch; Springer: Berlin/Heidelberg, Germany, 1928; pp. 631–660. [Google Scholar]
  55. Fenu, G.; Malloci, F.M. Forecasting plant and crop disease: An explorative study on current algorithms. Big Data Cogn. Comput. 2021, 5, 2. [Google Scholar] [CrossRef]
  56. Zhan, T.; Xu, H. Nonlinear optimization of GM(1,1) model based on multi-parameter background value. In Computer and Computing Technologies in Agriculture, Proceedings of the 5th IFIP TC 5, SIG 5.1 International Conference, CCTA 2011, Beijing, China, 29–31 October 2011; Springer: Berlin/Heidelberg, Germany, 2011; pp. 15–19. [Google Scholar]
  57. Arlo Richardson, E.; Seeley, S.D.; Walker, D.R. A model for estimating the completion of rest for ‘Redhaven’ and ‘Elberta’ peach trees. HortScience 1974, 9, 331–332. [Google Scholar] [CrossRef]
  58. Khaghaninia, S.; Mohammadi, S.A.; Nejad, K.H.I.; Sarafrazi, A.M. Population variation of codling moth Cydia pomonella (Lep.; Tortricidae) based on molecular data from northwestern Iran. Turk. J. Zool. 2011, 35, 571–578. [Google Scholar] [CrossRef]
  59. Zhou, X.-Y.; Wang, J.-W.; Zhang, J.; Yang, H.; Yang, W.; Yang, C.-P.; Zhou, W.; Li, T.; Wang, R.-L. The relationships among longitude, latitude and elevation of occurrence of Parocneria orienta Chao (Lepidoptera: Erebidae) in China. J. Kans. Entomol. Soc. 2020, 92. [Google Scholar] [CrossRef]
  60. Lacueva-Pérez, F.J.; Ilarri, S.; Barriuso, J.J.; Balduque, J.; Labata, G.; del-Hoyo, R. Grapevine phenology prediction: A comparison of physical and machine learning models. In 24th International Conference on Big Data Analytics and Knowledge Discovery, Proceedings of the 24th International Conference, DaWaK 2022, Vienna, Austria, 22–24 August 2022; Springer: Cham, Switzerland, 2022; pp. 263–269. [Google Scholar] [CrossRef]
  61. Lacueva-Pérez, F.J.; Ilarri, S.; Labata Lezaun, G.; Barriuso Vargas, J.J.; del Hoyo Alonso, R. Multifactorial Evolutionary Prediction of Phenology and Pests: Can Machine Learning Help? In Proceedings of 16th International Conference on Web Information Systems and Technologies (WEBIST 2020), Online, 3–5 November 2020; pp. 75–82. [Google Scholar] [CrossRef]
  62. Lessio, F.; Alma, A. Models applied to grapevine pests: A review. Insects 2021, 12, 169. [Google Scholar] [CrossRef] [PubMed]
  63. RedFAra. RedFAra Homepage. Available online: http://web.redfara.es/ (accessed on 21 December 2022).
  64. QGIS Development Team. QGIS Geographic Information System. Open Source Geospatial Foundation Project. 2023. Available online: http://qgis.osgeo.org (accessed on 21 December 2022).
  65. Spanish Cadastral Registry. Spanish Cadastral Registry Electronic Homepage. Available online: http://www.catastro.minhap.es/webinspire/index.html (accessed on 21 December 2022).
  66. Government of Aragon. Aragon Open Data Homepage. Available online: https://opendata.aragon.es/ (accessed on 21 December 2022).
  67. Silva, S.; Burgardt, M.; Ohanian, G.; San Blas, G.; Valeria, B.; Buran, M.; Musri, L. Aplicación del modelo de Touzeau como herramienta para el control y manejo de la “polilla de la vid” Lobesia botrana (Denis & Schiffmüller) en Mendoza. In Proceedings of the Séptimo Encuentro de Investigadores y Docentes de Ingeniería, EnIDI 2013, Mendoza, Argentina, 12–14 September 2013; pp. 398–407. [Google Scholar]
  68. SiAR. Agro-Climatic Information System for Irrigation Homepage. Available online: https://eportal.mapa.gob.es//websiar/Inicio.aspx (accessed on 21 December 2022).
  69. Copernicus. Copernicus Climate Data Store. Available online: https://cds.climate.copernicus.eu/cdsapp#!/dataset/reanalysis-era5-single-levels?tab=overview (accessed on 21 December 2022).
  70. Open-Meteo. Open-Meteo. 60 Years of Historical Weather as Free API and Download. Available online: https://openmeteo.substack.com/p/60-years-of-historical-weather-as (accessed on 21 December 2022).
  71. Akyuz, A. Growing degree day calculation method comparison between two methods in the northern edge of the US Corn Belt. J. Appl. Serv. Climatol. 2015, 2015, 1–9. [Google Scholar] [CrossRef]
  72. Nietschke, B.S.; Magarey, R.D.; Borchert, D.M.; Calvin, D.D.; Jones, E. A developmental database to support insect phenology models. Crop Prot. 2007, 26, 1444–1448. [Google Scholar] [CrossRef]
  73. Gabel, B.; Mocko, V. A functional simulation of European vine moth Lobesia botrana Den. Et Schiff. (Lep., Tortricidae) population development. J. Appl. Entomol. 1986, 101, 121–127. [Google Scholar] [CrossRef]
  74. Del Tío, R.; Martínez, J.L.; Ocete, R.; Ocete, M.E. Study of the relationship between sex pheromone trap catches of Lobesia botrana (Den. & Schiff.) (Lep., Tortricidae) and the accumulation of degree-days in Sherry vineyards (SW of Spain). J. Appl. Entomol. 2008, 125, 9–14. [Google Scholar] [CrossRef]
  75. Gallardo, A.; Ocete, R.; López, M.A.; Maistrello, L.; Ortega, F.; Semedo, A.; Soria, F.J. Forecasting the flight activity of Lobesia botrana (Denis & Schiffermüller) (Lepidoptera, Tortricidae) in Southwestern Spain. J. Appl. Entomol. 2009, 133, 626–632. [Google Scholar] [CrossRef]
  76. Montáns, F.J.; Chinesta, F.; Gómez-Bombarelli, R.; Kutz, J.N. Data-driven modeling and learning in science and engineering. C. R. Mécanique 2019, 347, 845–855. [Google Scholar] [CrossRef]
  77. Domingues, T.; Brandão, T.; Ferreira, J.C. Machine learning for detection and prediction of crop diseases and pests: A comprehensive survey. Agriculture 2022, 12, 1350. [Google Scholar] [CrossRef]
  78. Bonhomme, R. Bases and limits to using ‘degree.day’ units. Eur. J. Agron. 2000, 13, 1–10. [Google Scholar] [CrossRef]
  79. Rodriguez-Galiano, V.F.; Sanchez-Castillo, M.; Dash, J.; Atkinson, P.M.; Ojeda-Zujar, J. Modelling interannual variation in the spring and autumn land surface phenology of the European forest. Biogeosciences 2016, 13, 3305–3317. [Google Scholar] [CrossRef]
  80. Scikit-Learn. Scikit-Learn Homepage. Available online: https://scikit-learn.org/stable/ (accessed on 21 December 2022).
  81. Keras. Keras Homepage. Available online: https://keras.io/ (accessed on 21 December 2022).
  82. Optuna. Optuna Homepage. Available online: https://optuna.org/ (accessed on 21 December 2022).
  83. Python. Python Homepage. Available online: https://www.python.org/ (accessed on 21 December 2022).
  84. Magalhães, I.A.L.; de Carvalho Júnior, O.A.; de Carvalho, O.L.F.; de Albuquerque, A.O.; Hermuche, P.M.; Merino, É.R.; Gomes, R.A.T.; Guimarães, R.F. Comparing machine and deep learning methods for the phenology-based classification of land cover types in the Amazon biome using Sentinel-1 time series. Remote Sens. 2022, 14, 4858. [Google Scholar] [CrossRef]
  85. Mendes, J.; Peres, E.; Neves dos Santos, F.; Silva, N.; Silva, R.; Sousa, J.J.; Cortez, I.; Morais, R. VineInspector: The Vineyard Assistant. Agriculture 2022, 12, 730. [Google Scholar] [CrossRef]
  86. Sitarz, M. Extending F1 metric, probabilistic approach. arXiv 2022, arXiv:2210.11997. [Google Scholar]
Figure 1. Different stages of development of Lobesia botrana: (a,b) Eggs; (ce) larvae; (f) pupae; (g) adults. Photo credit: Julio Prieto-Díaz.
Figure 1. Different stages of development of Lobesia botrana: (a,b) Eggs; (ce) larvae; (f) pupae; (g) adults. Photo credit: Julio Prieto-Díaz.
Plants 12 00633 g001
Figure 2. Types of damage caused by Lobesia botrana and associated plant diseases: (a) Inflorescence; (b,c) berries; (d,e) Botrytis cinerea. Photo credit: Julio Prieto-Díaz.
Figure 2. Types of damage caused by Lobesia botrana and associated plant diseases: (a) Inflorescence; (b,c) berries; (d,e) Botrytis cinerea. Photo credit: Julio Prieto-Díaz.
Plants 12 00633 g002
Figure 3. Map of the study area with the three protected designations of origin (Cariñena, in violet color; Campo de Borja, in green color; and Calatayud, in red color), in the province of Zaragoza, region of Aragón, Northeastern Spain. Map generated with QGIS [64]. Source: Government of Aragón-Geographical Institute of Aragón (https://idearagon.aragon.es/portal/, accessed 21 December 2022).
Figure 3. Map of the study area with the three protected designations of origin (Cariñena, in violet color; Campo de Borja, in green color; and Calatayud, in red color), in the province of Zaragoza, region of Aragón, Northeastern Spain. Map generated with QGIS [64]. Source: Government of Aragón-Geographical Institute of Aragón (https://idearagon.aragon.es/portal/, accessed 21 December 2022).
Plants 12 00633 g003
Figure 4. (Left) Map with the location of the monitoring sites; (right) map with the locations of the weather stations. Source: Government of Aragón-Geographical Institute of Aragón (https://idearagon.aragon.es/portal/, accessed 21 December 2022).
Figure 4. (Left) Map with the location of the monitoring sites; (right) map with the locations of the weather stations. Source: Government of Aragón-Geographical Institute of Aragón (https://idearagon.aragon.es/portal/, accessed 21 December 2022).
Plants 12 00633 g004
Figure 5. Delta traps placed in the field for flight monitoring: (a,b) Traps; (c) gummed bottom; (d) sex pheromone; (e) adults captured.
Figure 5. Delta traps placed in the field for flight monitoring: (a,b) Traps; (c) gummed bottom; (d) sex pheromone; (e) adults captured.
Plants 12 00633 g005
Table 1. Description of the best-performing artificial neural network (ANN).
Table 1. Description of the best-performing artificial neural network (ANN).
LayerNeuronsActivation Function
096linear
1336selu
264linear
3144selu
Table 2. Comparison of observed (Obs.) flight peak days of the year (DOY) with the DOY predicted by the Touzeau (Tou.) and machine learning (ML) models for the 2008−2011 seasons.
Table 2. Comparison of observed (Obs.) flight peak days of the year (DOY) with the DOY predicted by the Touzeau (Tou.) and machine learning (ML) models for the 2008−2011 seasons.
SeasonGenerationObservation/
Prediction
Site Code
50024A0030001350073A0350009450004A0340017550073A0070004350073A0100000950073A0270001450073A0340012850073A0410002150073A0480001650073A0580004250073A0840003150073A0890001150073A0950007350073A1000005150098A0300002750201A0020012450201A0220001250268A00100050
20081stObs.126126 153 153 126 153
Tou.111111 111 111 111 111 111
ML126 14993 149 153 126 153
2ndObs.184184157159 157157157159157157160159157 159 157
Tou.172172172172 172172172172172172172172172 172 172
ML186 157157 157157157159157157 159157 159 156
3rdObs.221221224216 224224224221224224224214 216 224
Tou.210210210210 210210210210210210210210210 210 210
ML220 223215 223223223221223223 214223 216 223
20091stObs.134125131128 131131131125131131131128 128 131
Tou.127127127127 127127127127127127127127 127 127
ML 130128 130130130124130130130128 128
2ndObs.175173152158 152152152159152152152152152 159 152
Tou.166166166166 166166166166166166166166166 166 166
ML 151157 151151151159151151151152 151
3rdObs.223212217215 215217217215217215220215218 223 224
Tou.201201201201 201201201201201201201201201 201 201
ML 217215 215217217214217215220215 223
20101stObs.139139152122 157157124 131124152127124 156127124
Tou.116116116116 116116116 116116116116116 116116116
ML123123152122 120120123120123123152123123 156123123
2ndObs.181181158163 158158183157158183158163183 163160158
Tou.173173173173 173173173173173173173173173 173173173
ML159175157162 157157175156157175157162175 162159157
3rdObs.224221218224 218224224213214214214224 224224224
Tou.207207207207 207207207207207207207207 207207207
ML223223218223 217223223212214214214223 223223223
20111stObs.123123131128129131128128128128128128128128129121128128
Tou.103103103103103103103103103103103103103103103103103103
ML123123 126126130126126126126126126126126126121126126
2ndObs.164164157184157158158158 182157182 157157 182
Tou.162162162162162162162162 162162162 162162 162
ML 152 157157157 176150175 150150 150176
3rdObs.216216224220223214220220220224220220223213223 223213
Tou.202202202202202202202202202202202202202202202 202202
ML216216 216222210220220220223220220 213222 222213
Information for each of the site codes is available in Tables S4 and S5.
Table 3. Hyperparameters used for training the ML model.
Table 3. Hyperparameters used for training the ML model.
HyperparameterValue
Number of layers3−10
Number of neurons in intermediate layers0, 32, 512
Number of neurons in last 2 layers16, 464
Activation functions“selu”,”linear”,”tanh”,”softmax”
Exist activation functionsigmoid
Learning rates10.0 × 10−3, 10.0 × 10−2, 10.0 × 10−1
Optimizers‘sgd’,’adam’,’rmsprop’
CallbackVal_loss, patient = 17.0, min_delta = 0.17
Epochs1000
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Balduque-Gil, J.; Lacueva-Pérez, F.J.; Labata-Lezaun, G.; del-Hoyo-Alonso, R.; Ilarri, S.; Sánchez-Hernández, E.; Martín-Ramos, P.; Barriuso-Vargas, J.J. Big Data and Machine Learning to Improve European Grapevine Moth (Lobesia botrana) Predictions. Plants 2023, 12, 633. https://doi.org/10.3390/plants12030633

AMA Style

Balduque-Gil J, Lacueva-Pérez FJ, Labata-Lezaun G, del-Hoyo-Alonso R, Ilarri S, Sánchez-Hernández E, Martín-Ramos P, Barriuso-Vargas JJ. Big Data and Machine Learning to Improve European Grapevine Moth (Lobesia botrana) Predictions. Plants. 2023; 12(3):633. https://doi.org/10.3390/plants12030633

Chicago/Turabian Style

Balduque-Gil, Joaquín, Francisco J. Lacueva-Pérez, Gorka Labata-Lezaun, Rafael del-Hoyo-Alonso, Sergio Ilarri, Eva Sánchez-Hernández, Pablo Martín-Ramos, and Juan J. Barriuso-Vargas. 2023. "Big Data and Machine Learning to Improve European Grapevine Moth (Lobesia botrana) Predictions" Plants 12, no. 3: 633. https://doi.org/10.3390/plants12030633

APA Style

Balduque-Gil, J., Lacueva-Pérez, F. J., Labata-Lezaun, G., del-Hoyo-Alonso, R., Ilarri, S., Sánchez-Hernández, E., Martín-Ramos, P., & Barriuso-Vargas, J. J. (2023). Big Data and Machine Learning to Improve European Grapevine Moth (Lobesia botrana) Predictions. Plants, 12(3), 633. https://doi.org/10.3390/plants12030633

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