Melon (Cucumis melo L.; Cucurbitaceae) is one of the 10 most-consumed fruits in the world. The whitefly Bemisia tabaci (Gennadius) (Hemiptera: Aleyrodidae) is one of the most important pests in melon crops worldwide. Conventional sampling plans are the starting point to establish decisionmaking systems for integrated pest management (IPM) programs. The purpose of this study was to determine a conventional sampling plan for B. tabaci in melon crops with plants at the vegetative, flowering, and fruiting stages. The best sampling units for B. tabaci were the 5th and 6th most apical leaves of the plant vine. The best sampling technique was direct counting of adult whiteflies. The most appropriate frequency distribution to describe B. tabaci densities in melon fields was the negative binomial. Whiteflies on melon fields with plants at different phenological stages showed a common aggregation parameter (Kcommon = 0.9134). The optimal number of samples from the sampling plan was 72 samples per field with a maximum error of 15% in population determination. The sampling plan determined by this study can be used by farmers because it is a low-cost (US$5.27 per sampling), fast (39 minutes per sampling) and feasible (15% maximum evaluation error). The same sampling plan can be used with melon plants at the vegetative, flowering, and fruiting stages.
Melon (Cucumis melo L.; Cucurbitaceae), often called muskmelon or honeydew melon, is one of the 10 most-consumed fruits in the world (3.6% of global fruit consumption). Global production is 29.6 million tons per year (FAO 2014). In Brazil, cultivation of C. melo varies from 55 to 75 days and during this period it can generate a profit of US$6,364 per ha, as well as 1 direct and 3 indirect jobs per ha of cultivation (Freitas et al. 2009).
One of the principal pests occurring in melon crops is the whitefly Bemisia tabaci (Gennadius) (Hemiptera: Aleyrodidae) (Picanço et al. 2002; Braga Sobrinho et al. 2011; Pessarakli 2016), which is a polyphagous species that occurs in all cultivated regions worldwide (Stansly & Naranjo 2010; CABI 2016). This pest has approximately 700 reported host plants, mainly from the families Asteraceae, Brassicaceae, Cucurbitaceae, Fabaceae, Malvaceae, and Solanaceae (Greathead 1986; Li et al. 2011). It causes damage to crops due to sap sucking, injection of toxins into the vascular system, and transmission of viruses, especially geminiviruses (Blackmer & Byrne 1999; Toscano et al. 2004; Varma et al. 2011). Additionally, opportunistic fungi develop on the sweet excretions (honeydew) of B. tabaci, forming a dark layer of fungi that covers the plants (sooty mold) resulting in a reduction of photosynthesis and depreciation of the fruit value (Gusmão et al. 2005; Stansly & Naranjo 2010).
Sampling plans are the starting point for decision-making systems in integrated pest management (IPM) programs because they are used to determine economic damage levels and to validate sequential sampling plans (Binns et al. 2000; Gusmão et al. 2006; Pedigo & Rice 2014; Rosado et al. 2014; Pereira et al. 2017). In conventional sampling plans, the number of samples is fixed, whereas in sequential plans this number is variable (Pereira et al. 2017; Pinto et al. 2017). Conventional sampling plans consist of a sample unit, sampling technique, and the number of samples (Gusmão et al. 2005; Moura et al. 2007; Bacci et al. 2008). Optimal sampling units and techniques must be accurate, fast, low-cost, and representative. Accurate sampling units and techniques have low variability. Thus, for units and techniques to be precise they must have a relative variance (RV) of less than 25% (Southwood 1978). On the other hand, representative units and sampling techniques are those that have relative densities (individuals per sample) that present the same variation as the absolute density (individuals per plant or vine). The number of samples in the sampling plans should be economically feasible. Therefore, sampling plans must be simple, fast, and cost effective in order to be viable for farmers (Gusmão et al. 2005; Moura et al. 2007; Bacci et al. 2008; Pinto et al. 2017).
Despite the importance of B. tabaci as a pest in melon crops, a sampling plan for use in IPM programs for this species has yet to be developed for melons. Therefore, the objective of this study was to determine a conventional sampling plan for B. tabaci in melon crops with plants during vegetative, flowering, and fruiting stages.
Materials and Methods
EXPERIMENTAL CONDITIONS
This study was carried out during 2015 and 2016 in 12 different commercial melon fields (variety ‘Valenciano’) in Formoso do Araguaia, Tocantins State, Brazil (11.9021055°S, 49.5616027°W, with an altitude of 240 masl, and a tropical climate with a dry winter). Each field consisted of approximately 20 ha. The plants were spaced by 1.5 × 2.0 m and established according to Braga Sobrinho et al. (2008) and Pessarakli (2016). The cultural treatments used were plant and fruit thinning, weed control, and fruit turning.
The research was divided into 3 parts. In the first part, the optimal sampling unit to assess B. tabaci was selected using the direct counting technique. In the second part, the optimal sampling technique (direct counting or leaf beating on a white tray) was selected. In the third part, the appropriate number of samples for the sampling plan was determined.
SELECTION OF THE OPTIMAL SAMPLING UNIT TO ASSESS BEMISIA TABACI
This part of the study was carried out in 3 melon fields at different stages in 2015. The plants were evaluated in the vegetative, flowering, and fruiting stages (Figs. 1A, B, C). In each field, 100 plants were selected at random and in each plant, the largest vine was evaluated. The leaves were numbered according to their position on the branch. Leaf number 1 was the most apical of the vine, number 2 the 2nd most apical leaf, and so on (Fig. 1D). On each leaf, the number of B. tabaci adults was evaluated by direct count. In this evaluation, the leaf was handled carefully (to avoid insect escape) and the number of B. tabaci adults present was counted (Fig. 1E). The adults were evaluated because this is the insect stage that causes most damage to plants (Gusmão et al. 2006) and they are easier to count due to their larger size (Chu et al. 2003).
The sampling unit was selected based on the criteria of leaf occurrence frequency in a vine, precision, and representativeness. In the first criterion, only leaves with an occurrence frequency (frequency with which a leaf in a certain position appears in a vine) above 80% were selected. This was done due to the long time taken to locate leaves that have a low occurrence frequency on the vines, which can increase sampling time (Rosado et al. 2014; Pinto et al. 2017). The occurrence frequency of each leaf on the vines of melon plants was calculated using formula (1):
where Fi = frequency of occurrence of leaf i (%) in the 100 vines evaluated, i = position of the leaf on the vine (left to right), Ni = number of times that the leaf i was present on the 100 evaluated vines, and Nr = total number of evaluated vines.
The precision criterion used to identify the sampling unit was the relative variances of B. tabaci densities on each leaf of the vine. Only relative variances of less than 25% were selected because samples and techniques with relative variance higher than 25% generate unfeasible sampling plans (Southwood 1978). The relative variances of B. tabaci densities on each leaf of the vine were calculated using formula (2) (Moura et al. 2003; Bacci et al. 2006):
where RV = relative variance (%), S(x) = standard error of the mean densities, and x = mean densities.
Under the representativeness criterion, the relative densities (adults per leaf-1) and the absolute densities (adults per vine-1) were calculated, and the sample units considered were those with relative densities that were significantly correlated (P < 0.05) with absolute densities (Rosado et al. 2014; Pinto et al. 2017). The final selection was made based on the sample unit being amongst those suitable for the evaluation of pests according to all criteria during the vegetative, flowering, and fruiting stages.
SELECTION OF THE OPTIMAL SAMPLING TECHNIQUE
This part of the study was performed in 3 melon fields in 2015. During the evaluations, the plants of each field were in the vegetative, flowering, or fruiting stages (Figs. 1A, B, C). Twenty plants were randomly selected from each field. Two vines were evaluated in each plant. The evaluated leaves were those previously selected (fifth and sixth leaves). With 1 of these vines, B. tabaci density was evaluated by the direct count technique (Fig. 1E). With the other vine of the plant, whitefly densities were evaluated by leaf beating onto a white tray (40 × 25 × 3 cm) and counting B. tabaci adults (Fig. 1F). These techniques were used because they are suitable for the evaluation of sucking insects such as B. tabaci (Moura et al. 2007; Pinto et al. 2017). The time required to execute each sampling technique also was evaluated.
The relative variance of B. tabaci densities obtained by each sampling technique was calculated. The calculation of relative variance was performed using formula (2). The data for B. tabaci densities and sampling times for each sampling technique were analyzed by analysis of variance at P < 0.05. The technique that provided a relative variance of less than 25% and a shorter sampling time was selected (Southwood 1978).
DETERMINATION OF THE NUMBER OF SAMPLES FROM THE SAMPLING PLAN
This stage of the study was divided into 2 parts. In the first part, the most appropriate frequency distribution was determined to describe the B. tabaci densities in the melon fields. In the second part, the maximum error to be used in pest sampling was determined.
DETERMINATION OF THE FREQUENC Y DISTRIBUTION THAT BEST DESCRIBES THE DENSITIES OF BEMISIA TABACI IN CROPS
This stage of the research was carried out in 12 melon fields during 2015 and 2016. During the evaluations, the plants of each field were in the vegetative, flowering, or fruiting stages (4 fields at each phenological stage). In each field, 300 plants were evaluated. These plants were located equidistantly to eliminate directional errors (Bacca et al. 2006; Rosado et al. 2014; Pereira et al. 2017; Pinto et al. 2017). The sampling techniques used were those previously selected (direct counting of B. tabaci adults on the fifth or sixth most apical leaves of the vine).
The observed and expected frequencies of B. tabaci densities were calculated according to the negative binomial, Poisson, and positive binomial distributions through Microsoft Excel spreadsheets. These frequencies were compared using the chi-square test (c2) (Young & Young 1998). The data for the B. tabaci densities in a field were adjusted to a certain frequency distribution when the difference between the observed and expected frequencies was not significant (P > 0.05) (Young & Young 1998). A given frequency distribution was considered descriptive of B. tabaci densities in crops when the whitefly distribution conformed to this frequency distribution in most fields (more than 70%) (Rosado et al. 2014; Pinto et al. 2017).
The values of the aggregation parameter for each field were calculated using formula (3) (Young & Young 1998; Moura et al. 2003):
where k = aggregation parameter in the field, x = average density of the pest in the field, and S 2 = variance of the pest densities.
The aggregation parameters (k) of each of the 12 fields were assessed with simple linear regression analysis according to Bliss and Owen (1958). In this regression, it is considered that the fields present a common aggregation parameter (kc) when their intercept is not significant and the slope is significant by the F test at P < 0.05 (Rosado et al. 2014; Pinto et al. 2017).
Determination of the Maximum Allowed Error for Calculating the Number of Samples
Initially, the sample numbers of the sampling plan for melon crops in the vegetative, flowering, and fruiting stages were calculated using formula (4) (Young & Young 1998; Gusmão et al. 2005; Moura et al. 2007):
where NA = number of sample units, C = permitted error, x = population mean, and kc = common parameter of aggregation of the negative binomial frequency distribution (0.9134) determined previously.
In these determinations, we used error values of 0.05, 0.10, 0.15, 0.20, and 0.25. These errors were used because they are the standard range of errors used in these calculations (Moura et al. 2007; Pinto et al. 2017). Given the number of samples in these situations, the sampling time for each of these plans was evaluated in 3 fields. The plants of each field were in the vegetative, flowering, or fruiting stages.
The sampling cost was calculated using the following formula (5):
where CS = cost of 1 sampling, CM = cost of sampling material (pencil, rubber, paper, and clipboard: US$0.07 per sample), TS = sampling time (h), and CL = labor costs (salary of a rural worker and non-wage labor costs: US$8.00 per h) (Moura et al. 2007; Rosado et al. 2014).
The selected sampling plan was the plan with the lowest level of error and with sampling time of less than 1 h for 3 fields in the vegetative, flowering, or fruiting stages. This was done because sampling plans with this duration are feasible (Gusmão et al. 2005). In addition, these criteria allow these sampling plans to be standardized because they can be used in any melon plantation irrespective of phenological stage or pest density (Pinto et al. 2017).
Results
SELECTION OF THE PLANT COMPONENTS TO BE USED FOR BEMISIA TABACI SAMPLING
In plants during the vegetative stage, only the 7 most apical leaves displayed frequencies of occurrence on the vines higher than 80%. From those 7 leaves, only the fifth and sixth leaves had whitefly density relative variances lower than 25%. From the third to seventh leaves, the relative densities (adults per leaf-1) presented significant correlations (P < 0.05) with absolute densities (adults per vine-1) (Table 1). Therefore, in the vegetative stage, the fifth and sixth leaves were suitable to sample B. tabaci adults because they had frequency of occurrence above 80%, relative variance lower than 25%, and significant correlations with absolute densities.
In the flowering plants, only the 10 most apical leaves presented frequencies of occurrence on the vines higher than 80%. From the third to tenth leaves, the relative variance of B. tabaci densities was lower than 25%. From the fourth to tenth leaves, the relative densities (adults per leaf-1) presented significant correlations (P < 0.05). Thus, the fourth to tenth leaves of flowering melon plants are sampling units that may be used to sample B. tabaci adults.
In fruiting plants, the 15 most apical leaves presented frequencies of occurrence on the vines higher than 80%. From the fifth to fifteenth leaves the relative variance of B. tabaci densities was lower than 25% and the relative densities (adults per leaf-1) presented significant correlations (P < 0.05) with absolute densities (adults per vine-1) (Table 1). Thus, the fifth to fifteenth leaves of fruiting melon plants are sampling units that may be used to sample B. tabaci adults.
Therefore, the sampling units that were suitable in vegetative, flowering, and fruiting stages were the fifth and sixth most apical leaves of the vine. These were the ideal sampling units to generate a standardized sampling plan to be used for all phenological stages.
SELECTION OF THE OPTIMAL SAMPLING TECHNIQUE
In plants at vegetative and fruiting stages, higher densities of B. tabaci adults were detected using the leaf beating on a white tray technique than with the direct counting technique. The opposite occurred when the melon plants were at the flowering stage (Fig. 2A).
In plants at the vegetative stage, the relative variance of B. tabaci densities was lower than 25% when using both the leaf beating on a white tray technique and the direct counting technique. In the flowering and fruiting plants, the relative variance of B. tabaci densities was lower than 25% only when the direct counting technique was used (Fig. 2B).
Sampling time using the direct counting technique was shorter (average of 7 min per stage-1) than with the leaf beating on a white tray technique (average of 14 min per stage-1) in melon plants for all 3 phenological stages (vegetative, flowering, and fruiting) (Fig. 2C).
Thus, the best technique to sample B. tabaci adults in melon plants in the 3 phenological stages (vegetative, flowering, and fruiting) was direct counting.
Table 1.
Selection of leaf position to be used in the sampling of Bemisia tabaci adults in melon plants during vegetative, flowering, and fruiting stages: leaf occurrence frequency on the vine, density (mean ± standard error), relative variance (RV), Pearson correlation coefficient (r) between relative densities (adults per leaf-1) and absolute density (adults per vine-1).
DETERMINATION OF THE NUMBER OF SAMPLES FOR BEMISIA TABACI SAMPLING IN MELON FIELDS
In 9 of the 12 fields evaluated (75%), the differences between observed and expected frequencies of B. tabaci densities according to the negative binomial distribution were not significant (P > 0.05). Therefore, in these fields, pest densities followed the negative binomial frequency distribution. Bemisia tabaci densities followed the Poisson frequency distribution in only 3 of the melon fields. On the other hand, B. tabaci densities that followed a positive binomial frequency distribution were not observed in any of the fields (Table 2). Therefore, the ideal formula to calculate samples for the sampling plan is the negative binomial frequency distribution.
The simple linear regression of the common aggregation parameter of the B. tabaci densities in 12 melon fields (Kcommon) as a function of the values of the aggregation parameter of each crop showed a significant curve (P < 0.05) and a non-significant intercept (P > 0.05) (Table 3). Therefore, the B. tabaci densities in 12 melon fields showed a common aggregation parameter (Kcommon).
The sampling time for sampling plans with errors of 5 and 10% was higher than 1 h. When the sampling error was 15, 20, and 25%, the sampling time was less than 1 h. Therefore, 15% was the smallest error that generated a sampling plan with sampling time less than 1 h. Using this error (15%), the number of samples from the sampling plan was 72 per crop with a duration of 39 min and cost of US$5.27 per sample (Table 4).
Discussion
The fifth and sixth most apical leaves of the plant were the optimal sample units for the evaluation of B. tabaci in melon crops because they allowed more precise sampling of the pest for all phenological stages. In this context, these units allowed faster sampling because they were present in 100% of plant vines. These sample units also enabled accurate sampling, because B. tabaci densities displayed a relative variance of less than 25%. According to Southwood (1978), these values generate feasible sampling plans. In addition, the fifth and sixth most apical leaves of the branch allowed representative samplings of pest attacks on the vine. This is shown by the fact that relative B. tabaci densities (adults per leaf-1) showed a positive and significant correlation with the absolute density of the pest (adults per vine-1). Moreover, these fast, precise, and representative samplings were observed for all phenological stages (vegetative, flowering, and fruiting) of melon plants. The preference of whitefly to attack younger leaves can be explained by the fact that these leaves present better nutritional content, besides having soft, fine cuticles, and a higher concentration of water available to the insect (Van Lenteren & Noldus 1990).
Table 2.
Chi-squared test (c2) results between the frequencies observed and expected by the negative binomial, Poisson, and positive binomial Bemisia tabaci density distributions (means ± standard errors) in 12 fields of melon.
The best sampling technique to sample B. tabaci adults in melon crops was direct counting due to the possibility of sampling the pest at all the phenological stages of plant growth. One of the characteristics that led to the choice of the direct counting technique was the fact that it enabled accurate sampling (relative variance lower than 25%) in all situations. In addition, direct counting was faster than the leaf beating on a white tray technique for all phenological stages of the plant. The fact that a technique has a shorter sampling time indicates that its use will generate more feasible sampling plans (Southwood 1978; Gusmão et al. 2005; Moura et al. 2007).
The frequency distribution that best described B. tabaci adult densities in melon crops was the negative binomial. This was probably due to variance of pest densities being higher than average for this characteristic (Tonhasca et al. 1994; Moura et al. 2003). Additionally, we observed that in melon crops with plants at vegetative, flowering, and fruiting stages, the variance in B. tabaci densities was up to 74.65 and 104.00 times higher than the averages, respectively.
In this context, some authors misinterpret the frequency distribution data, stating that it represents the spatial distribution of insects in fields (Maruyama et al. 2006; Costa et al. 2010; Souza et al. 2013). However, this is not the case, because there is no relationship between frequency distributions of insect densities and the spatial distribution of these organisms in fields (Young & Young 1990; Barrigossi et al. 2001). To evaluate spatial distribution of insect populations in fields, it is necessary to locate the geographical coordinates of each sampling point and to analyze the data using a suitable statistical tool such as geostatistics (Young & Young 1998; Rijal et al. 2014; Rosado et al. 2015). However, frequency distribution data is an appropriate tool to select the formula to use in the calculation of samples for sampling plans (Young & Young 1998; Gusmão et al. 2005; Bacci et al. 2008) and population estimates (Pedigo & Buntin 1993; Heersink et al. 2016).
Table 3.
Analysis of variance of Bemisia tabaci adult density, sampled in 12 melon fields, to assess the dispersion parameter of the common negative binomial distribution (Kcommon).
The existence of a common aggregation parameter (Kcommon) for B. tabaci densities in melon crops indicates that it is possible to determine a sampling plan for this pest that is applicable to all sampled crops (Bliss & Owen 1958). This is even more relevant because the fields sampled in this project were in different phenological stages (vegetative, flowering, or fruiting). Therefore, the sampling plan generated here has great potential to be adopted by farmers due to its capacity for being used throughout all melon crop stages.
For a B. tabaci sampling plan to be viable, it should present a sampling time of less than 1 h and require a small number of samples (Bacci et al. 2006; Rosado et al. 2014). The smallest error used in the calculation to sample this pest was 15%. Therefore, it is a feasible system. According to Southwood (1978), the permitted error for this type of study can be up to 25%. The number of samples obtained was 72 samples per crop with a collection time of 39 min and a cost of US$5.27 per sample. The advantage of these systems is that producers are able to obtain more reliable (smaller sampling error), representative, faster, and less expensive samples. According to Bacci et al. (2008), these sampling techniques can reduce pesticide use, cost, and application time, and preserve the ecosystem and human health. In addition to providing the farmer with environmentally friendly management options, IPM ensures greater profitability by reducing pesticide spraying (Picanço et al. 2004; Pedigo & Rice 2014).
In conclusion, B. tabaci sampling in melon crops should be performed by direct counting of pest adults on the fifth and sixth most apical leaves of the plant vine. The sampling plan is composed of 72 samples per crop. This sampling plan is feasible with a duration of 39 min, a cost of US$5.27 per sample, and can be used in plants during vegetative, flowering, and fruiting stages.
Table 4.
Number of samples, time involved, and cost of sampling Bemisia tabaci adults in melon plants as a function of the maximum error admitted in this process.
Acknowledgments
We thank the National Council for Scientific and Technological Development - CNPq - Brazil (Projects: 458946/2014-1 and 304178/2015-2) for financial support. We also thank the Coordination for the Improvement of Higher Education Personnel (CAPES), Brazil (Project: PROCAD-2013), and the Minas Gerais State Research Foundation (FAPEMIG) for the scholarships and resources provided. The Universidade Federal de Viçosa, Brazil, and the Federal University of Tocantins, Brazil, are thanked for providing logistical support.