Spatiotemporal Dynamics of *Bemisia tabaci* MEAM1 (Hemiptera: Aleyrodidae) in Commercial Soybean Crops

*Luciana Barboza Silva, Raimundo Henrique Ferreira Rodrigues, Thiago Ferreira Rodrigues, Maria Carolina Farias e Silva, Edivania de Araujo Lima and José Wellington Batista Lopes*

## **Abstract**

Spatiotemporal dynamics studies of crop pests enable the determination of the colonization pattern and dispersion of these insects in the landscape. Geostatistics is an efficient tool for these studies: to determine the spatial distribution pattern of the pest in the crops and to make maps that represent this situation. The aim of this study was to determine the spatiotemporal distribution of *B. tabaci* adults and nymphs in commercial soybean crops from planting to harvest using geostatistics. Infestation by adults and nymphs of *B. tabaci* started between 30 and 50 days after the emergence of the plants. The maximum population density of ten adults per plant and two nymphs per leaf was registered between 90 and 101 days after plant emergence. The colonization of soybean plants by *B. tabaci* may be divided into three stages: beginning infestation (at the outermost parts of the crop), whole area colonization, and dispersion colonization (when the whole crop area is infested). The density of adult insects was positively correlated with rainfall and relative humidity. Wind speed positively affected the dispersion of adult whiteflies. The distribution pattern of *B. tabaci* in the soybean crop was aggregated.

**Keywords:** Aleyrodidae, dispersion, Geostatistics, *Glycine max*, Hemiptera, whitefly

## **1. Introduction**

*Bemisia tabaci* (Gennadius) Middle East -Asia Minor 1 - MEAM1 (Hemiptera: Aleyrodidae), commonly known as whitefly, is one of the most invasive and destructive pests for global agriculture, causing significant yield losses in commercial crops [1, 2]. This insect has high fecundity and fertility rates. It is polyphagous and occurs in areas of temperate, subtropical, and tropical climates [3–5]. The whitefly has a wide host range; thus, crops such as cotton, beans, squash, melon, and tomato may serve as alternate host plants for this species infesting soybean crops.

The whitefly causes several problems in Brazilian crops, and the insect may cause direct damage by sucking phloem sap. Indirect damage is linked to the excretion of honeydew, which serves as a substrate for the growth of opportunistic fungi (Capnodium sp.) [6–8]. *B. tabaci* may inhibit the gathering of carotenoids and chlorophyll affecting the photosynthetic rate of plants [9]. In soybean [*G. max* (L.) Merr., Fabaceae], *B. tabaci* is the vector for the cowpea mild mottle virus (CpMMV) [10].

In recent years, soybean has been highly attacked by high whitefly populations, mainly in the reproductive stages (R1 stage) Until very recently, this species was only seen as an occasional pest, however it is now considered a key pest in soybean. Although the importance of the whitefly *B. tabaci* MEAM1 has increased its incidence in soybean cultivars [11–13], there is no sampling method currently defined for this pest in this crop. The recommendations do not specify the sample numbers, size of area to be sampled, location of the different stages, or the most suitable forms of evaluation. In general, management is based on information obtained from other countries and what has been observed in other crops, such as cotton [14–16] and melon [17], in the state of Arizona, USA. In Brazil [18].

Understanding the spatiotemporal dynamics of pest insects in crops provides important information that can be incorporated into integrated pest management programs [19–22]. The study of the spatiotemporal dynamics of insects using geostatistics can provide important information about the pattern of colonization, aggregation, and dispersion of pests in the field [20, 22–25]. This statistical tool uses a method that characterizes spatial variation by comparing similarities between distant and proximal points.

This technique provides results that allow colonization maps to be produced, zoning the different densities and determining the pattern of spatial distribution of the insects in the field [20, 25]. A sequence of these maps during crop development may indicate areas that demand greater attention for the observation of pest infestation over time. In these spatiotemporal dynamic's studies, plants in the landscape that can serve as pest sources can be identified and located. This may be useful for farmers when planning the location of crops in the landscape or even for predicting where pest attacks will start within the crop area [24, 26, 27].

Considering the importance of studying the spatiotemporal dynamics of insect pests and to the best of our knowledge, there are few published studies on *B. tabaci* in soybean. Thus, the aim of this study was to determine the spatiotemporal distribution of adults and nymphs of *B. tabaci* in soybean crops from planting to harvest by using geostatistics and verifying abiotic factors (temperature, relative humidity, wind, and rainfall) that affect the dispersion of these insects.

## **2. Materials and methods**

The experiment was carried out at Celeiro Seeds Farm, and the farm's total area was estimated at 16,000 hectares, where grains and soybean seeds are produced. The farm is located in the municipality of Monte Alegre do Piauí and Serra do Quilombo in the Brazilian state of Piaui (09°21′12"S; 45°07′42" W, 642 m) (**Figure 1A**). The region's soil is characterized as a yellow Latosol (Oxisol)28, and the climate is tropical with a dry season in winter, Aw, according to the Köppen-Geiger classification 29. Climatic data such as the mean air temperature (°C), relative humidity (%), wind speed (m/s) and wind direction were obtained from NIMET (National Institute

*Spatiotemporal Dynamics of* Bemisia tabaci *MEAM1 (Hemiptera: Aleyrodidae)… DOI: http://dx.doi.org/10.5772/intechopen.105460*

#### **Figure 1.**

*(A) Brazil's map - Monte Alegre do Piauí, Piauí. (B) sample grid of the evaluated points and geographic location. (C) sampling unit (3 plants equidistant from the central point, total of 12 plants per unit).*

of Meteorology), and rainfall data were obtained with a rainfall meter during the period of execution of the experiment (**Figure 2**). This work was carried out for 2 seasons (2017/2018 and 2018/2019). During the first season, the cultivars Brasmax BonoIPRO® and BRS 9180IPRO® were used; in the second season, the cultivars were Brasmax Bonus IPRO® and Brasmax Extreme IPRO®. The characteristics of the cultivars are shown in **Table 1**.

The area of 1.6 ha (16,000 m<sup>2</sup> ), with 100 m length and 160 m width, was divided into 160 sampling units, with 10 m<sup>2</sup> . In the center of each sample unit, sets of four ordered points, equally spaced 2 m from a central point, were evaluated to improve the estimate of the nugget effect. Each sampling point was georeferenced with the aid of a GPS device model GPSMAP 60CS®x (Garmin). The distance from each field to the nearest native forest is described in **Table 1**.

**Figure 2.**

*Rainfall, temperature and average relative humidity of the region. Celeiro seed farm in Serra do Quilombo, Monte Alegre do Piauí. National Institute of meteorology (INMET) data. \* sowing month of cultivars.*

#### **2.1 Evaluation of the** *B. tabaci* **density**

The density of insects was evaluated in the sampling grid of 160 points previously georeferenced (**Figure 1B**). At each georeferenced point, 12 plants were evaluated; four subsamples of three plants were taken (sampling unit). The subsamples were 2.0 m equidistant from the central point (**Figure 1C**). The evaluated plants were

*Spatiotemporal Dynamics of* Bemisia tabaci *MEAM1 (Hemiptera: Aleyrodidae)… DOI: http://dx.doi.org/10.5772/intechopen.105460*


#### **Table 1.**

*Phenological cycle, fertilization, plant population, planting date, distance from cultivars to the nearest native forest, and evaluation days by cultivar in the two seasons studied.*

positioned along a regular grid pattern throughout the crop cycle to obtain systematic sampling points and avoid directional trends.

The third leaflet from top to bottom (apical third) of each plant was evaluated by direct counting. Leaves were handled with care, and the nymphs and adults of *B. tabaci* were counted. A magnifying glass was used to count the number of nymphs. We evaluated these leaves in particular using a direct counting method because these are the ideal sample type and technique for assessing the density of *B. tabaci* nymphs and adults in soybean crops 14. The sampling schedule was standardized; the evaluations started at 7 am and ended at 11 am. The sampling schedules for each cultivar are described in **Table 1**.

#### **2.2 Analysis of the spatial distribution of** *B. tabaci* **in soybean**

The data on adult and nymph densities of *B. tabaci* were submitted to statistical analysis. Subsequently, principal component analysis was performed between the range and density of whitefly adults with climatic parameters (temperature, relative humidity, wind speed, and rainfall). These analyses were performed using R®

software. 30 All geostatistical procedures were performed using the Geostatistical Analyst Tool for ArcGIS 10.5 (ESRI). The procedure in this tool can be simplified as follows: map and examination of data; preprocess of data if necessary (transform, detrend, decluster); definition of spatial structure model; definition of search strategy; prediction of values at unsampled locations; quantification of uncertainty of the predictions; checking if the model produces reasonable results for predictions and uncertainties and, using the information in risk analysis and decision making.

Subsequently, geostatistical analysis was performed; the spatial patterns and interpolations were determined using the parameters of adjusted experimental semivariograms and ordinary kriging, respectively. The semivariograms were calculated from the primary collected data, allowing us to detect differences between pairs of sampled points in relation to the distances (this procedure was used to adjust the theoretical semivariogram). Once the semivariance increases, there is a spatial dependence relationship between the densities of *B. tabaci* and the sampled points. The nugget effect and the sill value were calculated for each of the adjusted models (spherical, exponential, and Gaussian).

The experimental semivariograms were adjusted to the theoretical models, and the selection was made based on the cross-validation parameters, in which the measured and estimated values were compared using the standard mean error (SME) and the root mean square standardized error (RMSSE).

The kriging indicator was applied to spatial distribution maps of the insect in the crops, with the objective of modeling the probability of unsampled locations exceeding the quality reference values (QRVs). The kriging indicator does not use normal distribution assumptions, because it transforms the original data into binary values and, consequently, into cumulative distribution functions from pre-established values, in this case the QRVs. 31, 32 The data were transformed into log-normal values to minimize distribution errors and meet the requirements of common kriging based on the rejection of the null hypothesis of the Kolmogorov-Smirnov test for normal distribution. Spatial variability was determined from isotropic and anisotropic semivariograms. Anisotropic calculations were performed in four directions (0, 45, 90 and 135°).

To obtain reliable estimates, the theoretical model needs to show SME values close to 0 and RMSSE close to 1.33 The Akaike information criterion (AIC) was used as the last selection criterion. The spatial dependence rate (SDR) was calculated according to the formula: [C0/(C0 + C1) × 100]. 34 A nugget effect less than or equal to 25% of the plateau was considered strong. The value was considered moderate when it was between 25 and 75% and weak when it was above 75%. Variables that showed SDR less than one unit were not considered.

Spatial distribution maps of *B. tabaci* adults and nymphs were prepared for each crop. In these maps, the predominant direction of the winds is shown.

### **3. Results**

According to the statistical analysis, the values of variance were within the range 0 to 7.52 and were predominantly positive, where the peak of the normality curve of the results was higher than the standard value, the distribution curves were classified with positive asymmetries, and the coefficient of variation values were predominantly greater than 30% (Online Resource 1). In the two seasons evaluated, infestation by adults and nymphs of *B. tabaci* generally started between 30 and 50 days after the

*Spatiotemporal Dynamics of* Bemisia tabaci *MEAM1 (Hemiptera: Aleyrodidae)… DOI: http://dx.doi.org/10.5772/intechopen.105460*

emergence of the plants. The maximum population density of 10 adults per plant and approximately two nymphs per leaflet were recorded 101 days after plant emergence (Online Resource 1 and **Figure 3**).

#### **Figure 3.**

*Occurrence of* Bemisia tabaci *adults and nymphs in the Brasmax BonoIPRO®, BRS 9180IPRO® and Brasmax ExtremaIPRO® cultivars for seasons I, II, and l.*

The results that showed a strong or moderate degree of spatial dependence (SDR) were subjected to geostatistical models and chosen based on cross-validation parameters using the standard mean error (SME) and the root mean square standardized error (RMSSE). Of the 126 models that were processed, 33 were selected, of which 13 were exponential, 10 Gaussian, 10 spherical, and 36 pure nugget effects (**Table 1**). All selected models were isotropic (i.e., the spatial autocorrelation was the same in all directions).

Differences in the variogram parameters (nugget, sill, and range) were observed in the three cultivars in the two sampled seasons. A nugget effect was observed in the variograms, always at a low density of *B. tabaci* adults or nymphs (**Figures 4** and **5**). The SDR of the models ranged from 0 to 100. This showed a significant effect of the nugget effect on the interpolation for some sampled days. These proportions displayed that the spatial component accounted for 66% of the total spatial variance. From the selected models, 15.15% demonstrated a strong SDR (<0.25), 51.52% moderate SDR (between 0.25 and 0.75), and 33.33 weak SDR (**Table 2**).

The ranges of the models varied from 13.03 to 132.53 m, and the maximum range obtained for adults was 67 m at 104 DAE in the Brasmax ExtremaIPRO® cultivar and approximately 130 m for nymphs 101 DAE in the BRS 9180IPRO® cultivar (**Table 2**).

We observed that colonization of soybean plants by *B. tabaci* may be divided into three stages: the start of infestation, colonization of the area, and dispersion in the area. The adults of *B. tabaci* began to infest the experimental area from the outermost area; in season I, the insects came from the west, and in season II, the insects originated from the north (**Figure 4**).

Colonization by *B. tabaci* nymphs occurred as adults disperse in the area. Near the end of the culture cycle, adults, and nymphs of *B. tabaci* had already colonized the whole experimental area. With the colonization of the total area, adults began the process of migration to nearby areas (**Figures 4** and **5**).

The soybean cultivars were divided into three groups according to the density of the pest. In the first group (Brasmax Bonus IPRO®), the lowest densities of the pest (2 adults and 1 nymph per leaf) were observed. In the second group (Brasmax ExtremaIPRO®), lands with a higher density of nymphs (up to 10 nymphs per leaf) were observed. In the third group (BRS 9180IPRO®), plants exhibited the highest densities of adults (up to 15 adults per leaf) (**Figures 4** and **5**).

We observed two patterns of variation in the densities of *B. tabaci* throughout the development of the plants. The first occurred in Brasmax Bonus IPRO®, where the pest population varied little over time and remained at low density. The second pattern occurred in Brasmax ExtremaIPRO® and BRS 9180IPRO®, where the pest population increased over time and was distributed throughout the crop area (**Figures 4** and **5**).

The mean air temperature was high (29–33°C) during the cropping period. On the other hand, rainfall ranged from 10 to 127 mm per month during the cropping period (**Figure 2**). The winds during the experimental periods occurred predominantly in the northeast direction in 70% of cases with a maximum speed of 2.1 m/s. Performing the main component analysis, it was apparent that the two axes explained 74.79% of the data variation. There were positive correlations between the range and wind

speed, and this trend was observed also, in addition, between adult whitefly densities and rainfall and relative humidity (**Figure 6**).

#### **Figure 4.**

*Kriging maps for the spatial variability of* Bemisia tabaci *adults in the Brasmax BonoIPRO®, BRS 9180IPRO® and Brasmax ExtremaIPRO® cultivars for the 2017/2018 season I and 2018/2019 season II.*

**Figure 5.**

*Kriging maps referring to spatial variability of* Bemisia tabaci *nymphs in the Brasmax BonoIPRO®, BRS 9180IPRO® and Brasmax ExtremaIPRO® cultivars for the 2017/2018 season I and 2018/2019 season II.*


#### *Spatiotemporal Dynamics of* Bemisia tabaci *MEAM1 (Hemiptera: Aleyrodidae)… DOI: http://dx.doi.org/10.5772/intechopen.105460*


*bThreshold. c(C0/C0 + C1) \* 100) = Spatial Dependence Rate. dRoot Mean Square Standardized Error; and: Standard Mean Error.*

## **Table 2.**

*Models and parameters estimated by semivariogram for* Bemisia tabaci *adults and nymphs in the cultivars Brasmax BonoIPRO®, BRS 9180IPRO® and Brasmax ExtremaIPRO.*

*Spatiotemporal Dynamics of* Bemisia tabaci *MEAM1 (Hemiptera: Aleyrodidae)… DOI: http://dx.doi.org/10.5772/intechopen.105460*

**Figure 6.**

*Principal component analysis of the effects of rainfall and climatic elements on the range and abundance of*  Bemisia tabaci *in soybean crops. The first and second axes explained 74.79% of the variance.*

### **4. Discussion**

The data presented here provide information on the spatiotemporal dynamics of *B. tabaci* in commercial soybean crops. Under the conditions evaluated, the initial foci of colonization by adults occurred near the end of the vegetative stage of the crop, and it is common to observe the occurrence of nymphs and adults of *B. tabaci* at this stage of plant development . The peak infestation by adults and nymphs recorded in the present study was relatively low compared to other studies [22, 28].

The cultivar Brasmax BonoIPRO® had high attractiveness for adults and nymphs of *B. tabaci*. However, there was low infestation in the two evaluated seasons, which may indicate that cultivars of an early cycle with early planting may escape the whitefly attack [29], as mentioned above. *B. tabaci* infestation was more severe in late February and early March.

*B. tabaci* showed an aggregated distribution in the soybean crops, indicating that there were factors that influenced this outcome. The spatial distribution of pest insects in crops is a consequence of the colonization and dispersion of these organisms at these sites [24, 30]. Climate elements, terrain topography, characteristics of insect species, and landscapes are among the factors influencing the spatial distribution of pest insects in crops [11, 24, 31–34].

Some results obtained for *B. tabaci* adults and nymphs showed a weak degree of spatial dependence. This finding is characterized as a pure nugget effect (total absence of spatial dependence), probably due to the low infestation by adults and whitefly nymphs in each evaluation, making the presence of insects in geostatistical analyses

imperceptible [25, 35]. These results are a common finding in this type of study; other authors with experiments on the spatial distribution of *Euschistus heros* (Fabricius, 1798) (Hemiptera: Pentatomidae) in soybean have made similar observations [25, 36].

According to the semivariograms made, whitefly adults and nymphs had an isotropic distribution. This finding confirms the premise that *B. tabaci* disperses across the field in a diffusion pattern, allowing us to achieve an equal or similar number of adults throughout the sampling points regardless of direction, with the center of origin as the beginning of the evaluation [3]. Isotropic and aggregate pattern distributions were observed for adults and nymphs, and a reduction in range (dispersion radius) was observed as the density of adults per plant increased, while for nymphs, the opposite effect was verified. Other authors have observed the same pattern of whitefly distribution in soybean and other crops [3, 20, 22]. In experiments focused on the spatial distribution of *B. tabaci* in commercial watermelon crops, it was demonstrated that the dispersion of the insect was not influenced by a physical barrier or height gradient [20].

With the kriging maps, it was possible to verify that colonization of the area by *B. tabaci* proceeded in the western and northern regions. Near the experimental area, there was a permanent preservation area, native forest. Therefore, when favorable environmental conditions existed, the infestation of *B. tabaci* in soybean crops was to be expected. Migration from the native forest to the crop occurred when there was abundant food and shelter [26]. However, the northeast winds were predominant in both seasons; that is, the adults of *B. tabaci* flew against the wind searching for food. This fact is common since wind helps in the dispersion of volatiles in the area. Insects use volatiles as olfactory cues to find food, so insects can fly upwind, cross-wind and downwind when having perceived volatiles [37].

Kriging and semivariogram maps are fundamental tools that can assist in the identification and determination of the spatiotemporal dynamics of pests, in addition to providing information such as the dispersion pattern and range of insects. Based on this information, an adequate and efficient pest management monitoring program may be defined [20, 22, 25]. In addition, the elaboration of the maps allows the identification of the places with the highest incidence of whiteflies in crops, as well as the delimitation of the direction of colonization by adults. This allows the application of efficient control measures [20].

Geostatistics allowed us to verify the movements of *B. tabaci* adults and nymphs during the evaluations and soybean cycle. As previously explored, there are factors that can affect the distribution of pests in the field, features linked to insect characteristics such as population growth (reproduction, mortality) and dispersion (immigration, colonization, emigration) [3, 22, 25, 27, 36]. *B. tabaci* may fly to a height of up to 7 m and a distance of 7 km [38]. Climatic factors such as rainfall and the relative humidity of the air favored an increase in the density of whitefly adults, while wind speed favored the dispersion of adults, according to the literature [3, 20]. The development of monitoring and management techniques aimed at the precise control of *B. tabaci* depends on a better understanding of the flight behavior of the whitefly, especially short-range migration. The present study provides important information about the pattern of aggregation and distribution of *B. tabaci*, which may be used in future work to assist in the sampling method.

### **Acknowledgements**

We thank the National Council for Scientific and Technological Development — CNPq — Brazil for financial support. We also thank the Coordination for the

*Spatiotemporal Dynamics of* Bemisia tabaci *MEAM1 (Hemiptera: Aleyrodidae)… DOI: http://dx.doi.org/10.5772/intechopen.105460*

Improvement of Higher Education Personnel (CAPES), Brazil and the Piauí State Research Foundation (FAPEPI) for the scholarships and resources provided. To the Federal University of Piauí, Brazil, for providing logistical support. We also thank the Celeiro Seeds Farm for the supply of the area and inputs.

## **Author details**

Luciana Barboza Silva\*, Raimundo Henrique Ferreira Rodrigues, Thiago Ferreira Rodrigues, Maria Carolina Farias e Silva, Edivania de Araujo Lima and José Wellington Batista Lopes Graduate Program in Agricultural Sciences, Federal University of Piauí, Bom Jesus, Piauí, Brazil

\*Address all correspondence to: lubarbosabio@ufpi.edu.br

© 2022 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## **References**

[1] Lima CH, Sarmento RA, Pereira PS, et al. Feasible sampling plan for *Bemisia tabaci* control decision -making in watermelon fields. Pest Management Science. 2017;**73**:2345-2352

[2] Oliveira C, Auad A, Mendes S, et al. Economic impact of exotic insect pests in Brazilian agriculture. Journal of Applied Entomology. 2013;**137**:1-15

[3] Byrne DN, Rathman RJ, Orum TV, et al. Localized migration and dispersal by the sweet potato whitefly, *Bemisia tabaci*. Oecology. 1997;**105**:320-328

[4] Cruz P, Baldin E. Performance of *Bemisia tabaci* biotype B on soybean genotypes. Neotropical Entomology. 2017;**46**:210-215

[5] De Barro PJ, Liu SS, Boykin LM, et al. *Bemisia tabaci:* A statement of species status. Annual Review of Entomology. 2011;**56**:1-19

[6] Li JM, Ruan YM, Li FF, et al. Gene expression profiling of the whitefly (*Bemisia tabaci*) Middle East – Asia minor 1 feeding on healthy and tomato yellow leaf curl China virus -infected tobacco. Insect Science. 2011;**18**:11-22

[7] Stansly PAN, Naranjo SE, Brown JK, et al. *Bemisia*: Bionomics and Management of a Global Pest. Vol. 1. Netherland: Springer; 2010. p. 540

[8] Stansly PAN, Steven E. *Bemisia:* Bionomics and Management of a Global Pest. Vol. 1. Dordrecht, Springer Netherlands: Imprint, Springer; 2010. p. 558. ISBN: 978-90-481-2459-6

[9] Masuda K, Kato MT, Saito T. Reduction in carotenoid and chlorophyll content induced by the sweet potato

whitefly, *Bemisia tabaci*. Scientia Horticulture. 2016;**200**:102-104

[10] Marubayashi JM, Yuki VA, Wutke EB. Transmission of cowpea mild mottle virus by the whitefly *Bemisia tabaci* biotype B to bean and soybean plants. Summa Phytopathologica. 2010;**36**:158-160

[11] Czepak C, Coelho ASG, Rezende JM, et al. *Bemisia taba* ci MEAM 1 population surveys in soybean cultivation. Entomology Experimentalis and Applicata. 2018;**166**:215-223

[12] Da Silva Oliveira CE, Carneiro DEF, Toscano LC, et al. Population dynamics of *Bemisia tabaci* biotype B (Gennadius, 1889) in transgenic soybean. Revista de Agricultura Neotropical. 2018;**5**:1-5

[13] Vieira S, Bueno A, Boff M, et al. Resistance of soybean genotypes to *Bemisia tabaci* (Genn.) biotype B (Hemiptera: Aleyrodidae). Neotropical Entomology. 2011;**40**:117-122

[14] Ellsworth PC, Martinez-Carrillo JL. IPM for *Bemisia tabaci:* A case study from North America. Crop Protection. 2001;**20**:853-869

[15] Ellsworth PC, Diehl JW, Dennehy TJ, et al. Sampling Sweet Potato Whiteflies in Cotton. Tucson, AZ, United States: College of Agriculture and Life Sciences, Cooperative Extension IPM Series, University of Arizona; 1995. p. 2

[16] Naranjo SE, Flint HM. Spatial distribution of adult *Bemisia tabaci* (Homoptera: Aleyrodidae) in cotton and development and validation of fixedprecision sampling plans for estimating population density. Environmental Entomology. 1995;**24**:261-270

*Spatiotemporal Dynamics of* Bemisia tabaci *MEAM1 (Hemiptera: Aleyrodidae)… DOI: http://dx.doi.org/10.5772/intechopen.105460*

[17] Palumbo J, Tonhasca AD Jr, Byrne D. Sampling Plans and Action Thresholds for Whiteflies on Spring Melons. Online. Coop. Ext. IPM Series 1. Publ. No. 194021. Tucson, Arizona: University of Arizona, College of Agriculture and Life Sciences, Cooperative Extension; 1994 Available from: http://cals.arizona.edu/ crops/vegetables/insects/wf/wfmelons. html

[18] Moura MFD, Picanço MC, Silva ÉD, et al. Sampling plan of the biotype B of *Bemisia tabaci* in the cucumber crop. Pesquisa Agropecuária Brasileira. 2003;**38**:1357-1363

[19] Galdino TVS, Ferreira DO, Santana Júnior PA, et al. The role of the beetle *Hypocryphalus mangiferae* (Coleoptera: Curculionidae) in the spatiotemporal dynamics of mango wilt. Journal of Economic Entomology. 2017;**110**:865-874

[20] Lima CH, Sarmento RA, Galdino TV, et al. Spatiotemporal dynamics of whitefly *Bemisia tabaci* (Hemiptera: Aleyrodidae) in commercial watermelon crops. Journal of Economic Entomology. 2018;**111**:1895-1903

[21] Naveen N, Chaubey R, Kumar D, et al. Insecticide resistance status in the whitefly, *Bemisia tabaci* genetic groups Asia-I, Asia-II-1 and Asia-II-7 on the Indian subcontinent. Scientific Reports. 2017;**7**:1-15

[22] Suekane R, Degrande PE, Melo EPD, et al. Spatial distribution of soybean plants infested with whitefly Bemisia tabaci (Gennadius, 1889) (Hemiptera: Aleyrodidae). Archives of the Biological Institute. 2018;**85**:1-6

[23] Felicio TNP, Costa TL, Sarmento RA, et al. Surrounding vegetation, climatic elements, and predators affect the spatial dynamics of *Bemisia tabaci* (Hemiptera: Aleyrodidae) in commercial melon

fields. Journal of Economic Entomology. 2019;**112**:2774-2781

[24] Martins JC, Picanço MC, Silva RS, et al. Assessing the spatial distribution of *Tuta absoluta* (Lepidoptera: Gelechiidae) eggs in open -field tomato cultivation through geostatistical analysis. Pest Management Science. 2018;**74**:30-36

[25] Pias OHDC, Cherubin MR, Berghetti J, et al. Sampling grids for monitoring *Euschistus heros* infestation in soybean. Bragantia. 2017;**76**:386-397

[26] Macfadyen S, Hopkinson J, Parry H, et al. Early-season movement dynamics of phytophagous pest and natural enemies across a native vegetation-crop ecotone. Agriculture, Ecosystems & Environment. 2015;**200**:110-118

[27] Sciarretta A, Trematerra P. Geostatistical tools for the study of insect spatial distribution: Practical implications in the integrated management of orchard and vineyard pests. Plant Protection Science. 2014;**50**:97-110

[28] Vieira SS, Bueno RCODF, Bueno ADF, et al. Different timing of whitefly control and soybean yield. Ciência Rural. 2013;**43**:247-253

[29] Ameren N, Amjad M, Murtaza B, et al. Biogeochemical behavior of nickel under different abiotic stresses: Toxicity and detoxification mechanisms in plants. Environmental Science and Pollution Research. 2019;**26**:10496-10514

[30] Rosado JF, Picanço MC, Sarmento RA, et al. Geostatistics as a tool to study mite dispersion in physic nut plantations. Bulletin of Entomological Research. 2015;**105**:381-389

[31] Carrière Y, Degain B, Liesner L, et al. Validation of a landscape-based model for whitefly spread of the cucurbit yellow stunting disorder virus to fall melons.

Journal of Economic Entomology. 2017;**110**:2002-2009

[32] Ludwig M, Schlinkert H, Meyhöfer R. Wind -modulated landscape effects on colonization of B Russels sprouts by insect pests and their syrphid antagonists. Agricultural and Forest Entomology. 2018;**20**:141-149

[33] Ruscha A, Valantin-Morisona M, Sarthou JP, et al. Effect of crop management and landscape context on insect pest populations and crop damage. Agriculture, Ecosystems & Environment. 2013;**166**:118-112

[34] Silva GA, Picanço MC, Bacci L, et al. Control failure likelihood and spatial dependence of insecticide resistance in the tomato pinworm, *Tuta absoluta*. Pest Management Science. 2011;**67**:913-920

[35] Reay-Jones FPF, Toews MD, Greene JK, et al. Spatial dynamics of *stink bugs* (Hemiptera: Pentatomidae) and associated boll injury in southeastern cotton fields. Environmental Entomology. 2010;**39**:956-969

[36] De Souza LA, Barbosa JC, Aguirre-Gil OJ, et al. Spatial distribution of nymphs and adults of *Euschistus heros* (Fabricius, 1794) (Heteroptera: Pentatomidae) in transgenic soybean cultivars of different maturing cycles. Semina: Agricultural Sciences. 2018;**39**:905-920

[37] Reinecke A, Hilker M. Plant semiochemicals – Perception and behavioral responses by insects. Annual Plant Reviews online. 2014;**47**:115-153

[38] Isaacs R, Byrne DN. Aerial distribution, flight behavior and egg load: Their inter -relationship during dispersal by the sweetpotato whitefly. Journal of Animal Ecology. 1998;**67**:741-750

## Section 2
