SciELO - Scientific Electronic Library Online

 
vol.22 número2Dispersión de adultos de Piezodorus guildinii (Hemiptera: Pentatomidae) entre cultivos de soja y de alfalfaEvaluación del riego por melgas en los suelos del sur del Uruguay. Uso del modelo WinSRFR, resultados preliminares índice de autoresíndice de materiabúsqueda de artículos
Home Pagelista alfabética de revistas  

Servicios Personalizados

Revista

Articulo

Links relacionados

Compartir


Agrociencia (Uruguay)

versión impresa ISSN 1510-0839versión On-line ISSN 2301-1548

Agrociencia Uruguay vol.22 no.2 Montevideo nov. 2018

https://doi.org/10.31285/agro.22.2.10 

Articles

Analysis of the Spatial Distribution of Cydia pomonella (L.) in Southern Uruguay Using Geostatistical Tools

Análisis de la distribución espacial de Cydia pomonella (L.) en el sur de Uruguay utilizando herramientas geoestadísticas

MV Calvo1 

F Duarte1  2 

A Borges3 

IB Scatoni1 

1Universidad de la República, Facultad de Agronomía, Departamento de Protección Vegetal. Avenida Garzón 780, 12900 Montevideo, Uruguay. Correo electrónico: vcalvo@fagro.edu.uy

2Ministerio de Ganadería, Agricultura y Pesca, Dirección General de Servicios Agrícolas, División Protección Agrícola. Avenida Millán 4703, 12900 Montevideo, Uruguay

3Universidad de la República, Facultad de Agronomía, Departamento de Biometría, Estadística y Computación. Avenida Garzón 780, 12900 Montevideo, Uruguay


Summary:

Geostatistics can support pest management by analyzing the spatial distribution of a pest at a local or regional level. Therefore, specific site pest control can be performed based on accurate spatio-temporal information. Cydia pomonella (L.) (Lepidoptera: Tortricidae), the most serious apple and pear pest, requires several insecticide applications to reduce injury to acceptable levels. The aim of this research was to study the temporal and spatial variation in distribution and abundance of C. pomonella in the main fruit-growing area of Uruguay. From 2007 to 2009, approximately 120 pheromone traps per year were installed in apple and pear orchards distributed over 50,000 ha, and georeferenced. Male adult captures were registered weekly from October to March. Spatial analysis of the captures was done, obtaining semivariograms for the accumulated captures analyzed by generation and by growing season. Maps were constructed based on significant models obtained, where the population level of the pest was estimated using ordinary kriging. The correlation range estimated was from 1720 to 2690 m. Hot spots of high population level and some areas with comparatively low populations were constant over the 2-year period.

Keywords: codling moth; Lepidoptera; Tortricidae; semivariogram; pest maps

Resumen:

La geoestadística puede ser utilizada como una herramienta para el manejo de plagas mediante el análisis de la distribución espacial de una plaga tanto a nivel local como regional. Por lo tanto, el control de plagas en un sitio específico se puede lograr en base a información espacio temporal. Cydia pomonella (L.) (Lepidoptera: Tortricidae), la plaga de manzanas y peras más importante, requiere de varias aplicaciones de insecticidas para reducir los daños a niveles aceptables. El objetivo de esta investigación fue estudiar la variación temporal y espacial de la distribución y abundancia de C. pomonella en la principal área frutícola de Uruguay. Desde 2007 a 2009, se instalaron y georeferenciaron aproximadamente 120 trampas de feromonas por año en montes de manzanos y perales distribuidos en más de 50.000 ha. Las capturas de adultos machos se registraron semanalmente de octubre a marzo. Se realizó un análisis espacial de las capturas, obteniéndose semivariogramas para las capturas acumuladas analizadas por generación y por temporada. Los mapas se construyeron sobre la base de los modelos significativos obtenidos, donde el nivel de población de la plaga se estimó utilizando un kriging ordinario. El rango de correlación estimado fue de 1720 a 2690 m. Las zonas de mayores densidades poblacionales y algunas áreas con poblaciones comparativamente bajas fueron constantes durante el período de dos años de estudio.

Palabras clave: Carpocapsa; Lepidoptera; Tortricidae; semivariograma; mapa de plagas

Introduction

Cydia pomonella (L.) (Lepidoptera: Tortricidade) known as codling moth, is the main pest of apple, Malus domestica (Borkhausen), and pear, Pyrus communis L., but also attacks quince, Cydonia oblonga L., and walnut, Juglans regia L. The larva penetrates the fruit to directly feed on the seeds. In Uruguay, C. pomonella develops three generations by year, diapauses in the fifth instar larva occurring in the last generation (Núñez, 2006). Southern Uruguay has a heterogeneous agricultural orchard system which is composed of suitable host patches that vary in size, shape, species composition and management strategies. Apple and pear are the most important species present. Most control methods applied in these orchards are aimed at controlling apple scab and codling moth (Núñez & Scatoni, 2013). The systematic use of neurotoxic insecticides or their application based on forecasting systems which depend on the temporal variability of insect populations, fail to suppress the pest to acceptable damage levels (Núñez & Scatoni, 2013). Therefore, it is important to develop other management strategies.

Areawide pest management is the most effective approach to manage pests, reducing the use of insecticides and allowing alternative pest-control strategies that have a less negative impact on the environment and on farm growers. It involves synchronized and coordinated control measures over wide areas, instead of at farm levels (Faust et al., 2008). At the landscape scale, the fragmentation of farmland has resulted in a scattered resource distribution, so it is necessary to determine the final spatial pattern of a pest (Sciarreta & Trematerra, 2014). This is reflected in the spatial dynamics and dispersal of pest species from population sources toward population sinks. Understanding these patchy distributions can greatly influence the success of an area wide management programs (Caprio, Fave & Hankins, 2004). Therefore, to develop efficient area wide pest management strategies it is essential to know the insect population’s dynamic and spatial distribution, so control measures can be site-specific where populations are densest (Emmen, 2004). Analysis of the spatial distributions of insect populations is important for developing accurate and efficient sampling programs leading to more general pest monitoring programs (Southwood, 1978; Taylor, Woiwod & Perry, 1978; Taylor, 1984). Moreover, it helps to predict population density in unsampled locations (Liebhold, Rossi & Kemp, 1993), and thus to improve pest management strategies (Hughes & McKinlay, 1988), and to understand ecological relationships with different abiotic and biotic factors (Iwao, 1970; Hassel & May 1974). The development of two technologies has opened up new avenues for analyzing spatial patterns in insect populations: geostatistics and geographical information systems (GIS). These technologies constitute the basic tools of an approach to precision farming.

Geostatistics quantifies and models spatial and temporal correlations and enable the analysis of large data series in a very short time and thus carry out large-scale studies. This technique was first developed in the field of geology and has since been applied to agriculture (Emmen, 2004), and in a number of cases it has been employed to study the spatial distribution of insects (Midgarden, Youngman & Fleischer, 1993; Ribes-Dasi, Avilla & Bascuñana, 1998; Farias et al., 2004; Boiteau, 2005; Ramirez-Davila et al., 2005; Castillo et al., 2006; Moral García et al., 2004; Trematerra, Gentile & Sciarretta., 2004; Sciarretta & Trematerra, 2006, 2014; Basoalto et al., 2010; Calvo et al., 2011; Comas et al., 2012; Duarte, 2012; Duarte et al., 2015a, 2015b). Spatial analysis is based on the theory of regionalized variables put forward by Matheron (1970), whereby the value of a variable at a given point is related to the value of that variable at another point nearby (Rossi et al., 1992). On a small scale, biological variables usually have a chaotic structure, but on a larger scale, they have an identifiable spatial variability structure, and this can be modeled by the semivariogram using geostatistic techniques. The semivariogram is built up by graphically representing the semivariance of groups of pairs of data separated by similar distance intervals, in function of distance (Rossi et al., 1992). The semivariogram has three parameters. The nugget effect (Co) represents discontinuity at the origin, the sill (Co+C) is the maximum variability value, and the range is the limit of spatial correlation (Cuador, 2004). This information can be used to estimate the population by kriging geostatistical procedure. This estimation is based on the continuous variogram model. To predict a value in a not sampled point, it uses the existing knowledge of the data in its neighbourhood by taking account of the way that a property varies in space through the variogram model (Armstrong & Carignan, 1997; Journel & Huijbregts, 1978; Hohn, 1988; Isaaks & Srivastava, 1989; Rossi et al., 1992; Webster & Oliver, 2007). This makes possible to draw up population distribution maps that show the foci of greatest abundance, and this information can be used for defining pest-control strategies.

Forecasting systems based on the temporal variability of pests have proved to be useful, but they are limited when it comes to optimizing phytosanitary protection systems and even more for reducing the negative environmental impacts of pest-control measures in large areas (Faust et al., 2008). In Lerida, Spain, a forecast service to manage Cydia pomonella and Pandemis heparana on deciduous fruit has been developed using geostatistical analysis, allowing the determination of population levels at different areas (Avilla & Ribes-Dasi, 2004; Ribes-Dasi, Avilla & Bascuñana, 1998).

The aim of this research was to study the temporal and spatial variation in distribution and abundance of C. pomonella in the main deciduous fruit-growing area in southern Uruguay using geostatistical tools and relate it to local weather conditions.

Materials and Methods

Study area

The spatial distribution studies were carried out from 2007 to 2009 in the main deciduous fruit-growing area in southern Uruguay, in approximately 50,000 ha, located between NW x = 555577 y = 6176372; SE x = 576335 y = 6149388 (UTM WGS 84 21H). About 90 % of the fruit production of the country is distributed within the study area. The fruit production area includes 3508 ha of apple, 2220 ha of peach and nectarine, 1003 ha of pear, 344 ha of plum, and 287 ha of quince. Most of the crops are grown on small farms, less than 10 ha, and all species are commonly grown together on the same farms (Figure 1) (DIEA, 2010).

Figure 1: Map of the study area showing the fruit production, main sites, routes and roads, and the trap positions, which coincide with the distribution of the pear and apple orchards. 

Data collection

The activity of C. pomonella adult males was monitored weekly using wing traps baited with sex pheromone (Isca Tecnologias Ltda, 8.10 - (E,E) - dodecadien-1-ol, Codelmone- 1,4g/kg (0.14 %). Latex - 998.6g/kg (99.86 %)). Approximately 120 traps were installed in apple or pear orchards set between 500 to 1000 meters apart (Figure 1, Table 1). Traps were georeferenced and placed in the tree top at 1.8-2.0 m above ground. Pheromone dispensers were replaced every month and sticky bottom when necessary. Trapped adults were removed and counted weekly from October to March from 2007 to 2009.

Table 1: Analysis of accumulated captures of codling moth in pheromone traps for two years of study by each generation and over the entire growing season. 

Geostatistical analysis

Geostatistical analysis was performed according to Duarte et al. (2015b). The captures of male moths in traps were summed over each generation and over the entire growing season. All generations were delimited by the thermal constant (600 DD) using as lower threshold temperature 10 °C, and starting estimates of degree-days from October 1st until the end of March when capture of male moths ceased (overwintering generation, generation 1 and generation 2) (Núñez & Scatoni, 2013). The calculation of degree days and the analysis of weather conditions was based on records from the weather station in Instituto Nacional de Investigación Agropecuaria (INIA Experimental Station) located at Rincón del Colorado, Canelones (x = 560346 y = 6163235). Exploratory data analysis was performed, and the assumption of normality was verified using the Shapiro-Wilk test. The data were normalized using square root transformations except for the 2007-2008 accumulated growing season that logarithmic transformations were used. Presence of trends of the regionalized variable was also checked using regression linear models in all cases before adjusting empirical semivariograms to the data. We found no evidence of the presence of a trend.

Empirical semivariograms were constructed to provide a graphical representation of the semivariance of the captures as a function of distance, with a maximum of 15 lag intervals. The maximum distance for the structural analysis was 13,700 m, with small variations depending on the growing season, equivalent to half the distance between the two farthest points (Díaz-Viera, 2002). For each empirical semivariogram, two theoretical models were fitted using the mixed procedure of SAS software (SAS Inst. Inc., Cary, NC): the default model based on the lowest residual sum squares (RSS) and higher R squared (R2), and the spherical model, considered by some authors to be more adequate for biological studies (Moral García, 2004; Ramirez-Dávila et al., 2005). In all cases, only isotropic semivariograms were adjusted. In order to determine the presence/absence of spatial correlation of captures in each generation and in the entire growing season, the likelihood ratio test was conducted. It was applied to compare the fit of each semivariogram model to the null model (independent errors). This test expresses how many times more likely the data are under one model than the other, and it allows computing a p-value to decide whether to reject the null model in favor of the alternative model (spatially correlated errors). Afterward, based on the pre-selected spatial correlation models, interpolations were performed with ordinary kriging, and cross validation was applied to select models after kriging, using GS+ software. In this process, each sampled point was individually removed and its value estimated through kriging, obtaining for each trap the observed and the estimated value. To compare performance of statistical models, the Akaike Information Criterion (AIC) was used as described in Littell et al. (1996). The AIC is calculated based on log-likelihood values, taking into account the number of model parameters. Lower AIC values indicate better performance of the statistical model. For each model, the degree of spatial dependence was measured through the ratio between the nugget and sill parameters Co/(Co+C), and the range value was determined. The Co/(Co+C) ratio indicates the degree of uncertainty when the interpolation takes place (Table 2). The quotient indicates which proportion of the variability the model explains, and which proportion of the variability is independent of distance. If the Co/(Co+C) ratio value is less than 0.25, the level of spatial dependence is high, if it is between 0.26 and 0.75, the level is moderate, and if it is 0.76 or above, the level is considered low (Cambardella et al., 1994). Finally, the data were re-transformed and C. pomonella distribution maps were drawn based on the kriging estimations calculated using the GS+ software.

Table 2: Accumulated rainfall per month (in mm) during the study period in southern Uruguay. 

A six level scale with different shades of grey was used; the lighter the shade, the lower the population level. Scale shades were defined starting a threshold of one capture week; the intensity of gray was darkener when added up one capture week.

Results

Trap captures

Table 1 shows the classic statistics for the two study seasons and the variations of populations in time. Capture averages per trap were similar in both growing seasons; the second generation had the highest captures, which is reflected in the high captures registered in the next overwintering generation.

Population fluctuation

The fluctuation of the average populations of C. pomonella from mid-October to the end of March for each year is shown in Figure 2, where peaks of flights of adults are indicated (Overwintering generation, G1 and G2). Captures of C. pomonella during the two seasons followed a common pattern; however, for each generation the number of the insect caught and the timing of each generation were different. Most captures occurred in early October for the overwintering generation, while captures of the next generations were more extended in time (Table 1 and Figure 2).

Figure 2: Weekly average captures of Cydia pomonella in pheromone traps during the two seasons 

The 2007-2008 overwintering generation presented the lower accumulated captures. The autumn and early spring of 2007 were characterized by abundant rainfall (Table 2). The temperatures of the months of October and November were not favorable either. The minimum average temperature was 11 °C, while the next generations developed in more favorable conditions (Figure 3) (Núñez, 2006).

Figure 3: Weekly average of maximum and minimum temperatures according to records of INIA Las Brujas 

Captures during the 2008-2009 overwintering generation were higher than the previous growing season, which reflects in the establishment of the generation 1 (Table 2). Temperatures from October to November (spring) were very favorable for the insect growth (Figure 3), the average minimum temperature of the period was 13 °C, three degrees above the development threshold of the species which allowed a greater thermal accumulation (DD). Therefore, the first and second generation flights occurred almost 10 days before than the previous season (Table 1 and Figure 2). Maximum temperatures above 25 °C are an indicator that the twilight temperatures are above the minimum requirements for copulation and oviposition (Table 3), during the last week of November the temperatures were exceptionally high, over 30 ºC (Table 4). These conditions remained during December and January (summer), and precipitation was scarce. From the beginning of November until the end of January the average catch in pheromone traps did not reach three adults/week/trap. After February, more favorable temperatures were registered that allowed the reestablishment of populations. The population density of the overwintering generation is the reflection of what happened at the end of the previous season (Figure 2).

Table 3: Number of days with maximum temperatures below 25 °C in southern Uruguay. 

Table 4: Number of days with maximum temperatures above 30 °C in southern Uruguay. 

Spatial distribution

The spatial analysis performed with the accumulated capture data resulted in semivariograms for each generation and complete growing season (Figure 4). A significant spatial autocorrelation was observed with the spherical model showing the best fit in all cases, with ranges between 1.7 and 2.7 km. Among the models calculated for each season and each generation, we selected those with better adjustment according to the AIC criterion, and better R2. For the 2007-2008 growing season, significant models were obtained for the accumulated generations and for generation 1 (R2 = 0.876 and 0.778 respectively). For the 2008-2009 growing season all generations presented significant models with R2 higher than 0.56 except the G2 that, although the spherical model was significant, the R2 was 0, therefore it was not considered in the analysis. These models presented levels of high or moderate spatial dependence (Table 5).

Figure 4: Semivariogram models for each generation and complete growing season of Cydia pomonella

Table 5: Cydia pomonella models, range, spatial dependence level (Co/(Co+C) and determination coefficient (R2) obtained for each generation and season. 

The spatial distribution maps for the accumulated captures of codling moth between October and March of each season were obtained from the estimated models and were constructed following the grey scale color convention described previously (Figure 5). We can observe in the 2007-2008 season the large darkest grey areas. These large foci were mainly associated with the captures of the second generation, leaving a significant population remnant for the following spring.

The 2008-2009 season began with favorable conditions for the development of C. pomonella, and the overwintering generation made the greatest contribution to the accumulated capture (Figure 2). The dark grey foci observed in the maps correspond to the large abundance of that generation (Figure 5). The extreme drought that occurred from 2008 until February of 2009 plus the high temperatures registered could explain that the following generations of codling moth were smaller.

Figure 5: Interpolation maps of each generation and complete growing season of Cydia pomonella obtained by kriging based on significant statistically models. 

Discussion

Multiple factors may influence the spatial distribution and the population fluctuation of codling moth populations. The fluctuation observed during the two observation periods showed important differences, which were associated with the climatic conditions of each year. For example, the low densities of the 2007-2008 overwintering generation could be explained because, as was mentioned by González (2003), rain affects both adult flights and neonatal larval mortality. In contrast, autumn, winter, and spring of 2008 were characterized by a marked deficit in rainfall, which favored the emergence of adults of the overwintering generation and the establishment of first generation larvae in fruits (Núñez, 2006). Probably the conditions observed during the last week of November through December (temperatures over 30 ºC and scarce precipitation) became limiting for flights of the first generation 2008-2009, which would explain the low abundance. Twilight temperatures above 28 °C also affect the flights of this insect (Núñez, 2006). As discussed before, although climatic conditions seem to be one of the factors that explain some of the variations in the abundance of this species, the population density of the overwintering generation is the reflection of what happened at the end of the previous season. When the population of the second generation is high (G2 2007-2008), during the next spring captures of the overwintering generation increases greatly (Ow 2008-2009), the opposite happens in the following season, the low populations of the second generation (G2 2008 -2009) determine low captures of the overwintering generation in the following spring (Ow 2009-2010, data not presented) (Calvo et al., 2011). The foci located in the same areas throughout the period may be due to particular factors of the pest and the heterogeneous agro-ecosystem. C. pomonella is a low dispersal species (50 m). When the fruit is available in the orchard, this pest overwinters in the same site that was attacked in the previous season (Núñez, 2006). If the populations are high, this condition will persist at the orchard until the next season, except that important natural mortality factors (such as climatic conditions) operate on diapausing larvae (González, 2003). External factors that control pest abundance are: sanitary management intensity carried out to control this pest, hosts (cultivated species and varieties), handling of damaged fruits, proximity to fruit packing areas, abandoned orchards, public lighting on roads and paths, etc. (Ribes-Dasi, Avilla & Bascuñana, 1998). In this sense, the 2007-2008 season was characterized by abundant rains at the time of the sprouting of the pomes fruit, affecting the flights of the overwintering generation of C. pomonella and giving ideal conditions for the development of apple and pear scab. Due to the severity of this disease on fruits, growers ceased the pesticide sprays for diseases and insects. This explains in most cases the large darkest grey foci that are observed in the map of 2007-2008 season.

All factors mentioned above influence the creation of foci that could act as a source of insects or work as corridors to disperse them to other areas (Trematerra, Gentile & Sciarretta, 2004; Basoalto et al., 2010). These authors found that C. pomonella distribution appears to be structured as metapopulations in heterogeneous agro-ecosystems. On the other hand, the high populations recorded in pheromone traps were related to high percentages of damages at harvest, and the conventional control methods adopted by individual farmers in small areas were not sufficient to mitigate these damages. Therefore, codling moth spatial distribution and management should be analyzed at larger scales, not at the individual farm level, as proposed in area-wide integrated pest management (Vreysen, Robinson & Hendrichs, 2007; Knight, 2008).

Geostatistic analysis has proved to be a useful tool for area-wide pest management for different insect populations, including Ceratitis capitata (Sciarreta & Trematerra, 2010), Grapholita molesta (Sciarreta & Trematerra, 2006; Duarte, 2012; Duarte et al., 2015a) and Xylella fastidiosa (Castillo et al., 2006). The spatial analysis of the distribution of C. pomonella done in the present study showed that its population is aggregated and also that significant spatial autocorrelation was found in here. These patterns were also observed in previous studies done in other areas such as Central Chile (Basoalto et al., 2010), Italy (Sciarreta & Trematerra, 2014) and Spain (Ribes-Dasi, Avilla & Bascuñana, 1998; Comas et al., 2012). We adjusted isotropic models for all the cases, as it was reported by Comas et al. (2012) that C. pomonella spatial patterns were in general isotropic. The best significant model adjusted in all cases was the spherical model, as it was also observed in previous studies (Ribes-Dasi, Avilla & Bascuñana, 1998; Basoalto et al., 2010). The range of spatial dependence estimated in this research varies from 1720 m to 2690 m. As it was suggested by Duarte et al. (2015a) the model with the lowest range (1720 m) should be considered to define the optimal distance for designing a pheromone trap sampling system, under our conditions. Ribes-Dasi, Avilla y Bascuñana (1998) obtained similar ranges (1500 m) in Spain, but Basoalto et al. (2010) and Trematerra, Gentile & Scierreta (2004) found ranges from 150 m to 300 m. These last authors explained that male codling moth may not disperse as far as other studies indicated, perhaps due to landscape features influencing the tendency to disperse.

One of the outcomes of geostatistical analysis is the creation of risk assessment maps (Sciarreta & Trematerra, 2014). As it was mentioned before, these maps allowed detection of all foci with a high population of C. pomonella in order to make a rational use of insecticide and predict where high populations will develop next season based on population levels at the end of the previous one. Maps also provide a useful outcome to be included in the forecasting systems provided to farmers. The information given would help farmers to decide where and when to apply any control measure. In this sense, the data obtained in this research offer important tools to modify the alert systems in Uruguay. This information supports the idea that geoestatistical analysis is a very promising tool in area-wide integrated pest management and is an important contribution for defining the control strategies to be used in area-wide pest management in deciduous fruits of southern Uruguay (Zoppolo et al., 2016).

Acknowledgements:

The authors give special thanks to Dr. Manuel Ribes-Dasi from the Universidad de Lleida, Spain and to Dr. Jorge Franco from the Universidad de la República, Uruguay for their invaluable contributions in this research topic. The authors would also like to thank the Fondo de Promoción de Tecnología Agropecuaria of the Instituto Nacional de Investigación Agropecuaria (FPTA, INIA), which financed this research line, and the staff of the Servicio de Pronósticos Fitosanitarios of the Dirección Gener-a de Servicios Agrícolas of the Ministerio de Ganadería, Agricultura y Pesca for their cooperation (DGSSAA, MGAP).

References:

Armstrong, M. & Carignan, J. (1997). Géostatistique linéaire, application au domaine minier. Paris: École de Mines de Paris. [ Links ]

Avilla, J. & Ribes-Dasi M. (2004). Aplicaciones de la geoestadística y de los sistemas de información geográfica en el estudio de las poblaciones de las plagas. Phytoma, 164, 22-24. [ Links ]

Basoalto, E., Miranda, M., Knight, A. L. & Fuentes-Contreras, E. (2010). Landscape analysis of adult codling moth (Lepidoptera: Tortricidae) distribution and dispersal within typical agroecosystems dominated by apple production in Central Chile. Environmental Entomology, 39(5), 1399-1408. [ Links ]

Boiteau, G. (2005). Within-field spatial structure of colorado potato beetle (Coleoptera: Chrysomelidae) populations in New Brunswick. Environmental Entomology , 34, 446-456. [ Links ]

Calvo, M. V., Duarte, F., Borges, A. & Scatoni, I. (2011). Caracterización espacial de los lepidópteros plaga de los frutales de pepita en la zona sur de Uruguay. Montevideo: INIA. (Serie FPTA, 31). [ Links ]

Cambardella, C. A., Moorman, T. B., Novak, J. M. , Parkin, T. B., Karlen, D. L., Turco, R. F. & Konopka, A. E. (1994). Field-scale variability of soil properties in central Iowa soils. Soil Science Society of America Journal, 58(5), 1501-1511. [ Links ]

Caprio, M. A., Fave, M. K. & Hankins, G. (2004). Evaluating the impact of refuge with on source-sink dynamic between transgenic and non-transgenic cotton. Journal of Insect Science, 4, 3. [ Links ]

Castillo, A., Espinoza, J. C., Valle Mora, J. & Infante, F. (2006). Dispersión del parasitoide africano Phymastichus coffea Lasalle (Hymenoptera: Eulophidae) en un nuevo agroecosistema. Folia Entomológica Mexicana, 45, 319-327. [ Links ]

Comas, C., Avilla, J., Sarasúa, M. J., Albajes, R. & Ribes-Dasi, M. (2012). Lack of anisotropic effects in the spatial distribution of Cydia pomonella pheromone trap catches in Catalonia, NE Spain. Crop Protection, 34, 88-95. [ Links ]

Cuador, J. Q. (2004). Elementos de geoestadística. Available from https://www.monografias.com/trabajos14/geoestadistica/geoestadistica2.shtmlLinks ]

Díaz-Viera, M. A. (2002). Geoestadística aplicada. Available from http://mmc2.geofisica.unam.mx/cursos/geoest/GeoEstadistica.pdf. [ Links ]

DIEA. (2010). Encuesta frutícola de hoja caduca: Zafra 2009/10. Montevideo: MGAP. (Serie encuestas, 294). Available from http://www.mgap.gub.uy/sites/default/files/encuesta_fruticola_de_hoja_caduca_zafra_2009_-_2010_no294.pdfLinks ]

Duarte, F. (2012). Caracterización espacio-temporal de Grapholita molesta (Lepidoptera: Tortricidae) mediante métodos geoestadísticos y sistemas de información geográfica (Magister Thesis). Facultad de Agronomía, Universidad de la República, Montevideo, Uruguay. Available from http://biblioteca.fagro.edu.uy/iah/tesisposgrado/textostesis/2012/0074dua.pdfLinks ]

Duarte, F., Calvo, M. V., Borges, A. & Scatoni, I. B. (2015a). Geostatistics and geographic information systems to study the spatial distribution of Grapholita molesta (Busck) (Lepidoptera: Tortricidae) in peach fields. Neotropical Entomology, 44, 319-327. [ Links ]

Duarte, F., Calvo, M. V., Borges, A. & Scatoni, I. B. (2015b) Geostatistics applied to the study of the spatial distribution of insects and its use in integrated pest management. Revista Agronómica del Noroeste Argentino, 35(2), 9-20. [ Links ]

Emmen, D. (2004). La agricultura de precisión: Una alternativa para optimizar los sistemas de producción. Investigación y Pensamiento Crítico, 2, 68-74. [ Links ]

Farias, O. S., Roberto, R., Lopes, J. & Perecin, D. (2004). Geostatistical characterization of the spatial distribution of Xylella fastidiosa sharpshooter vectors on citrus. Neotropical Entomology , 33, 13-20. [ Links ]

Faust, R. M., Koul, O., Cuperus, G. & Elliot, N. (2008). General introduction to area wide pest management. In Areawide pest management: Theory and implementation (pp. 1-14). Beltsville: CAB International. [ Links ]

González, R. (2003). Las polillas de la fruta en Chile. Santiago de Chile: Universidad de Chile. (Serie Ciencias Agronómicas, 9). [ Links ]

Hassell, M. P. & May, R. M. (1974). Aggregation of predators and insect parasites and its effect on stability. Journal of Animal Ecology, 43, 567-594. [ Links ]

Hohn, M. E. (1988). Geostatistics and petroleum geology. New York: Van Nostrand Reinhold. [ Links ]

Hughes, G. & McKinlay, R. G. (1988). Spatial heterogeneity in yield-pest relationships for crop loss assessment. Ecological Modelling, 41, 67-73. [ Links ]

Isaaks, E. H. & Srivastava, R. M. (1989). Applied geostatistics. New York: Oxford University Press. [ Links ]

Iwao, S. (1970). Analysis of spatial patterns in animal populations: Progress of research in Japan. Review of plant protection research, 3, 41-54. [ Links ]

Journel, A. G. & Huijbrebts, C. J. (1978). Mining geostatistics. London: Academic. [ Links ]

Knight, A. (2008). Codling moth areawide integrated pest management. In O. Koul, G. Cuperus & N. Elliot (Eds.). Areawide pest management (pp. 159-190). London: CABI. [ Links ]

Liebhold, A., Rossi, R. & Kemp, W. (1993). Geostatistic and geographic information systems in applied insect ecology. Annual Review of Entomology, 38, 303-327. [ Links ]

Littell, R. C., Milliken, G. A., Stroup, W. W. & Wolfinger, R. D. (1996). SAS system for mixed models. Cary, North Carolina: SAS Institute Inc. [ Links ]

Matheron, G. (1970). La théorie des variables regionalisées et ses applications. Paris: Ecole de Mines de Paris. (Les Cahiers du Centre de Morphologie Mathematique de Fontainebleau, 5). [ Links ]

Midgarden, D. G., Youngman, R. R. & Fleischer, S. J. (1993). Spatial analysis of counts of western corn rootworm (Coleoptera: Chrysomelidae) adults on yellow sticky traps in corn: Geostatistics and dispersion indices. Environmental Entomology , 22, 1124-1133. [ Links ]

Moral García, F. J. (2004). Aplicación de la geoestadística en las ciencias ambientales. Ecosistema, 13, 95-105. [ Links ]

Moral García, F. J., Rodríguez Bernabé, J. A., Arranz Romero, A., de La Cruz Blanco, J. I. & Honoris Guido, F. (2004). Técnicas geoestadísticas aplicadas al análisis de la distribución de capturas de Helicoverpa armígera (Hübner) (Lepidoptera: Noctuidae) mediante trampas con feromonas sexuales en una plantación de tomate. Boletín de sanidad vegetal. Plagas, 30, 733-744. [ Links ]

Núñez, S. (2006). Cydia pomonella. In C. Bentancourt & I. B. Scatoni (Eds.) Lepidópteros de importancia económica en el Uruguay, reconocimiento, biología y daños de las plagas agrícolas y forestales (pp. 189-209). Montevideo: Hemisferio Sur. [ Links ]

Núñez, S. & Scatoni, I. (2013). Tecnología disponible para el manejo de plagas en frutales de hoja caduca. Montevideo: INIA . (Serie técnica, 210). Available from http://www.inia.uy/Publicaciones/Documentos%20compartidos/St%20210_2013.pdfLinks ]

Ramirez-Dávila, J. F., Gonzalez-Andujar, J. L., Lopez, M. A. & Ocete, R. (2005). Modelización y mapeo de la distribución espacial de ninfas del mosquito verde Jacobiasca lybica (Bergevin & Zanon) (Hemiptera, Cicadellidae) en viñedo. Boletín de sanidad vegetal. Plagas , 31, 119-132. [ Links ]

Ribes-Dasi, M., Avilla, J. & Bascuñana, M. (1998). Estudio de la distribución espacial de Cydia pomonella (L.) y Pandemis heparana (Denis&Schiffermüller) en Torregossa (Lleida) mediante métodos geoestadísticos. Boletín de sanidad vegetal. Plagas , 24, 935-948. [ Links ]

Rossi, R., Mulla, D., Journel, J. & Franz, E. (1992). Geostatistical tools for modeling and interpreting ecological spatial dependence. Ecological Monographs, 2, 277-314. [ Links ]

Sciarretta, A. & Trematerra, P. (2006). Geostatistical characterization of the spatial distribution of Grapholita molesta and Anarsia lineatella males in an agricultural landscape. Journal of Applied Entomology, 130, 73-83. [ Links ]

Sciarretta, A. & Trematerra, P. (2010). Spatio-temporal distribution of Ceratitis capitata population in a heterogeneus landscape in Central Italy. Journal of Applied Entomology , 35, 241-251. [ Links ]

Sciarretta, A. & Trematerra, P. (2014). Geostatistical tools for the study of insect spatial distribution: practical implications in the integrated management of orchard and vineyard pests. Plant Protection Science, 50, 97-110. [ Links ]

Southwood, T. R. E. (1978). Ecological methods. New York: John Wiley & Sons. [ Links ]

Taylor, L. R. (1984). Assessing and interpreting the spatial distributions of insect populations. Annual Review of Entomology , 29, 321-357. [ Links ]

Taylor, L. R., Woiwod, I. P. & Perry, J. N. (1978). The density dependence of spatial behavior and the rarity of randomness. Journal of Animal Ecology , 47, 383-406. [ Links ]

Trematerra, P., Gentile, P. & Sciarretta, A. (2004). Spatial analysis of pheromone trap catches of codling moth (Cydia pomonella) in two heterogeneous agroecosystems using geostatistical techniques. Phytoparasitica, 32, 325-341. [ Links ]

Vreysen, M. J., Robinson, A. S. & Hendrichs, J. (2007). Area-wide control insect pests: From research to field implementation. Dordrecht: Springer. [ Links ]

Webster, R. & Oliver, M. A. (2007). Geostatistics for environmental scientists (2nd ed.). Chichester: John Wiley & Sons. [ Links ]

Zoppolo, R., Scatoni, I. B, Duarte, F., Mujica, M. V. & Gabard, Z. (2016). Areawide pest management in deciduous fruits of southern Uruguay. Acta Horticulturae, 1137, 153-159. [ Links ]

Received: July 03, 2018; Accepted: July 11, 2018

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License