This manuscript describes a multispectral optical sensor that effectively detected damage to early season cotton artificially infested with varying densities of the two-spotted spider mite populations.
The objective of this study was to evaluate a ground-based multispectral optical sensor as a remote sensing tool to assess foliar damage caused by the two-spotted spider mite (TSSM), Tetranychus urticae Koch, on greenhouse grown cotton. TSSM is a polyphagous pest which occurs on a variety of field and horticultural crops. It often becomes an early season pest of cotton in damaging proportions as opposed to being a late season innocuous pest in the mid-southern United States. Evaluation of acaricides is important for maintaining the efficacy of and preventing resistance to the currently available arsenal of chemicals and newly developed control agents. Enumeration of spider mites for efficacy evaluations is laborious and time consuming. Therefore, subjective visual damage rating is commonly used to assess density of spider mites. The NDVI (Normalized Difference Vegetation Index) is the most widely used statistic to describe the spectral reflectance characteristics of vegetation canopy to assess plant stress and health consequent to spider mite infestations. Results demonstrated that a multispectral optical sensor is an effective tool in distinguishing varying levels of infestation caused by T. urticae on early season cotton. This remote sensing technique may be used in lieu of a visual rating to evaluate insecticide treatments.
Two-spotted spider mite, Tetranychus urticae (Koch) is a polyphagous and cosmopolitan pest of many field and horticultural plants1,2. It lives inside webbings in colonies on the bottom surface of the plant3,4. It has evolved from being a late season to an early season pest in the mid-southern United States over the last decade5. TSSM was the 5th most damaging pest of cotton and caused an estimated loss of 57,441 bales of cotton and 0.167% reduction in yield in the United States in 20116,7. Its short life cycle, high fecundity and haploid-diploid sex determination combined with the ability to digest and detoxify xenobiotics have exacerbated the development of resistance to pesticides8. Presently, acaricides remain as the only dependable control mechanism for the suppression of T. urticae. Therefore, field entomologists constantly evaluate currently available and newly developed acaricides for efficacy.
The estimation of damage by the spider mites is usually conducted by scoring the damage on a subjective scale because of the difficulty encountered in manually counting the mites. Some conducted binomial sampling, where only the proportion of infested leaves was scored rather than the number of spider mites per leaf9. A leaf reddening index scale, which varied from stippling and reddening to extensive reddening of the vegetation, was used as a criterion in estimating damage. The spatial distribution pattern of T. urticae on cotton conformed to a clumped distribution pattern9. Mites are distributed on cotton foliage from sparse to heavily clustered and remain so under field condition. Such distribution pattern coupled with its small size, mobility and prolific reproduction makes enumeration of TSSM difficult. Reliable alternative techniques are needed for the assessment of mite density in order to quantitatively evaluate efficacy of acaricides against TSSM.
The objective of this study was to separate cotton plants damaged by varying densities of TSSM by using an optical multispectral sensor. Our intent was to determine if the ground-based optical sensor could classify and separate healthy cotton plants from those damaged by the spider mites.
1. Establish TSSM Colonies on Pinto Beans
Figure 1: Planting pinto beans in plastic trays. Pinto beans seeds were planted in plastic trays (56 x 28 x 5 cm3) in the greenhouse and were watered daily. Please click here to view a larger version of this figure.
Figure 2: Pinto beans with trifoliate leaves. The first true leaf that is formed after the cotyledons emerge from the soil is the simple or unifoliate leaf. Subsequent leaves are the trifoliate leaves with denticulate tips. Please click here to view a larger version of this figure.
2. Transfer TSSM to Cotton Plants
Figure 3. Cotton plants with 4-5 true leaf stage. The cotyledons emerge from the soil as leaf-like structures oriented opposite to each other on the seedling stem. The apical meristem emerges through the cotyledons and form the first true leaves. Please click here to view a larger version of this figure.
Figure 4. TSSM masses hanging on a trifoliate bean leaf. TSSM lives in colonies and when populations reach high density, they form masses or boll-like structures and congregate at the leaf tips for dispersal. Please click here to view a larger version of this figure.
Figure 5. Cutting pinto bean leaf tips containing a TSSM mass with scissors. Trifoliate bean leaf tips containing TSSM masses were removed with scissors for infesting cotton plants. Please click here to view a larger version of this figure.
Figure 6. TSSM masses on bean leaf tips placed inside the pan. When sufficient number of trifoliate bean leaf tips with TSSM were found in the test plants, they were removed and placed inside the pan. These samples were used to infest treatment categories: light, medium and heavy which received 3, 20 and 40 masses of TSSM, respectively. Please click here to view a larger version of this figure.
Figure 7. Turning pan upside down. The pans containing trifoliate bean leaf tips with TSSM were turned upside down on cotton canopy to infest the test plants. Please click here to view a larger version of this figure.
3. Scan TSSM Infested Cotton Plants with the Multispectral Optical Sensor
Figure 8. A multispectral optical sensor used to quantitatively measure the health of cotton plants infested with varying density levels of TSSM. The wheeled push cart with the test plants were slowly traversed under the sensor head to obtain spectral reflectance values. a is the sensor head; b is pocket PC; c is battery compartment and input/output ports; d is RS-232 serial data cable and e is the craft paper to provide uniform background. Please click here to view a larger version of this figure.
4. Data Analyses
Figure 9. Percent reduction or change in NDVI relative to days after treatment. JMP software was used to graphically illustrate the functional relationship between percent change in NDVI relative to days of sampling (DAT). Please click here to view a larger version of this figure.
The sensor emits red and infrared light and this light in turn gets reflected back from the plant canopy. The reflected light serves as a quantitative measurement of the visible and the near-infrared bands of the electromagnetic spectrum and is numerically recorded as NDVI (Normalized Difference Vegetation Index) reading. The NDVI values vary from 0 to 0.99. The higher the NDVI reading the healthier the plant canopy. Healthy vegetation absorbs visible light and reflects the near-infrared light and unhealthy vegetation reflects more visible light and less near-infrared light. NDVI serves as a surrogate to photosynthetic activity and this spectral property is highly correlated to photosynthetically active radiation15,16,17,18. Very low values of NDVI (0.1 and below) corresponds to barren areas of rock or sand, while moderate values (0.2 to 0.3) represent shrub and grassland and high values (0.6 to 0.8) indicate temperate and tropical vegetation. The reflectance measurements were obtained under daytime lighting in the red and near-infrared regions of the spectra.
NDVI was calculated from the following equation: NDVI = (NIR – RED)/(NIR + RED), where RED and NIR are the spectral reflectance values (0-255) in the red and near-infrared spectrums at 660 and 770 nm, respectively. The sensor records reflectance values every 100 msec. The maximum NDVI readings culled from an array of NDVI values taken each time by the sensor were used in the analyses to minimize background reflectance and to provide consistently repeatable numerical values.
The analysis of variance of the data revealed that significant differences in percent reduction in NDVI were observed between lightly, medially and heavily infested cotton plants compared to the untreated control (F = 436.4; P <0.0001; df = 3, 32). The NDVI values which described plant vigor during the test period varied significantly between days of observation (F = 1398.2; P <0.0001; df = 8, 256). Also, the percent reduction in NDVI values were inversely correlated to the treatments during the observation period (DAT), but seldom tended to deviate from this pattern and interacting significantly with DAT (F = 201.5; P <0.0001; df = 24, 256). The MANOVA test criteria for no DAT effect was significant as well (Wilk's λ = 0.00913; F = 339.0; P <0.0001; df = 8, 25). Similarly, the interaction between DAT and treatment was significant (Wilk's λ = 0.00101; F = 29.8; P <0.0001; df = 24, 73).
Figure 9 shows percent change in plant vigor as demonstrated by NDVI values during the observation period. A positive percent change in NDVI value indicates healthy growing plants, while a negative value indicates that the vigor of the plant has decreased since the first measurement (i.e. Day 0) was made. The non-infested control plants showed increased vegetative growth throughout the course of the study, while the TSSM-infested plants showed degradation in health over time. Mean separation of the treatments shown in Table 1 reveals that no definable difference in percent reduction in NDVI between treatment categories (light, medium and heavy) and the control was observed until Day 5 when infestation classes significantly deviated from the control and remained overwhelmingly so thereafter. This data demonstrates that the optical sensor can be effectively used in lieu of labor-intensive manual sampling to assess treatment efficacy against acaricides on cotton.
Days After Treatment (DAT) | |||||||||
Infestation Category | 1 | 5 | 6 | 7 | 9 | 10 | 12 | 13 | 14 |
Control | 1.18 ±0.33a | 2.70 ±0.40a | 4.0 ±0.36a | 3.94 ±0.37a | 3.68 ±0.53a | 2.57 ±0.42a | 2.96 ±0.47a | 3.48 ±0.38a | 3.08 ±0.22a |
Light | -0.13 ±0.13b | -0.71 ±0.29b | -0.65 ±0.28b | -2.02 ±0.47b | -5.68 ±0.72b | -11.17 ±0.94b | -15.73 ±1.76b | -19.54 ±1.68b | -24.9 ±1.90b |
Medium | -1.83 ±0.42c | -7.06 ±0.63c | -9.61 ±0.53c | -10.39 ±0.57c | -17.06 ±0.80c | -26.92 ±0.72c | -33.84 ±0.96c | -37.05 ±1.14c | -41.74 ±0.73c |
Heavy | -0.97 ±0.58bc | -11.76 ±0.29d | -13.83 ±0.86d | -15.20 ±0.63d | -25.0 ±1.0d | -34.63 ±0.54d | -39.07 ±0.94d | -42.68 ±0.62d | -46.71 ±0.63 |
Table 1: Percent reduction in Max NDVI after cotton plants were infested with varying number of clusters or masses of TSSM. Cotton plants grown in plastic trays in the greenhouse were infested with three categories of spider mite density. Category Light received 3 masses or clusters of TSSM per tray, Category Medium received 20 masses per tray and Category Heavy received 40 masses per tray. Means were separated from control according to Duncan's Multiple Range Test (P = 0.05). Means followed by the same lower case letter were not significantly different at the 5% level of probability.
Conventionally, insecticide efficacy tests conducted in the field include several treatments of the chemical applied at different rates and compared to an untreated check. Acaricides with varying toxicity profiles against instar and adult stages of TSSM are assessed to determine whether the damage caused by them could be reduced by the chemical treatment. TSSM samples are collected and brought to the laboratory where they are examined under the microscope and the varying stages of TSSM are counted and recorded. It is critically important to take adequate samples of the host plant to determine the damage in each treatment and differentiate them with statistically acceptable precision. The number of samples that are required to differentiate the treatments from each other depends upon the distribution profile of the organism. Highly patchy distribution of the TSSM leads to a considerable amount of variation between sample areas, and many plants must be sampled in order to ensure reproducibility of the population estimates. However, budget, manpower, time and statistical precision are important factors that impinge upon sampling techniques. It behooves the researcher to allocate the available resources optimally to conduct sampling with the least cost but with most precision.
Entomologists in lieu of counting TSSM stages visually score the damage based upon a scale of no damage to varying levels of damage. For instance, some researchers proposed binomial sampling, where only the proportion of infested leaves was scored rather than the number of spider mites per leaf9,19. Others estimated damage by TSSM on cotton based on a leaf reddening index scale, which varied from stippling and reddening to extensive reddening of vegetation canopy19. These methods are arbitrary, anecdotal and biased by individual perceptions of the degree of damage. A more robust and quantitative assessment of the damage caused by TSSM is required to evaluate and separate the treatments with statistical precision.
The ground-based multispectral optical sensor appears to be an improved sampling tool to quantitatively determine the damage inflicted by TSSM and to separate the treatments more precisely than the visual damage scoring system adopted by many researchers. However, researchers have reported that the data intensive hyperspectral remote sensing provided numerous spectral signatures to identify and detect crop stresses and canopy characteristics compared to multispectral remote sensing which is less data intensive with two wave lengths20,21. Using a hyperspectral spectrometer, Reisig and Godfrey found that the NIR reflectance wavelength ≈850 nm as an informative spectrum in distinguishing arthropod-infested from uninfested cotton22. In this study, we have shown that the multispectral reflectance values (NIR value being ≈770 nm) with just two spectral bands were able to identify and characterize cotton plants infested with varying densities of TSSM. Also, we have reported earlier that the multispectral optical sensor not only effectively separated cotton plants infested with widely varying density categories of TSSM, but also showed that spiromesifen was more effective than abamectin in controlling TSSM in early season cotton in the greenhouse at one-half rate of the lowest label rate23.
The multispectral optical sensor could be mounted on a mobile research platform and the reflectance values could be obtained from the treated areas of the vegetation canopies without human subjectivity. Acaricidal efficacy data can thus be obtained without much human labor. The NDVI data can be easily uploaded to the computer and analyzed using commercial statistical software. A GPS receiver also can be mounted onto the mobile platform to collect the GPS coordinates in order to generate a damage map of the field. Using multispectral reflectance signatures from the plant canopy, the multispectral optical sensor provides a rapid and cost effective means of identifying and quantifying plant stress. Moreover, a much larger area of the field can be covered in less time with higher spatial resolution of plant canopy compared to conventional field scouting. It is important to remember that the threshold damage level for TSSM on cotton is variable from region to region in the United States. For example, TSSM damage would be higher in an arid environment such as in California compared to the Midsouth region where rainfall and high humidity often prevail24. Therefore, yield loss due to TSSM damage will be variable and so is the damage threshold. However, reports from Mississippi, Arkansas and Tennessee reveal that the threshold level for TSSM on cotton appears to be when 30 to 50% of the plants are infested and populations are actively increasing5. Furthermore, field grown plants are exposed to multiple stresses including water stresses and feeding of arthropod herbivores and the interactions between these activities can substantially reduce plant productivity and are likely to influence damage threshold. The TSSM can reduce stomatal conductance, photosynthesis and transpiration rate in cotton25. Plants grown in the greenhouse are influenced by UV light radiation and it does significantly influence stomatal function, photosynthesis and canopy morphology26,27 and could probably have an additive effect on plant stress. However, TSSM is capable of avoiding UV radiation by having access to habitat protected from solar radiation on the lower surface of the plant canopy28,29,30, where it resides.
The height of the optical sensor above the target canopy and the orientation of the sensor relative to the target are important factors that affected significantly the reflectance values obtained by the multispectral optical sensor31. For instance, when the mobile research platform traverses through a row crop field such as when cotton canopy is open, the sensor is likely to produce different results depending upon the orientation of the sensor, either parallel to or perpendicular to the row. It is also likely that soil and other background materials may greatly influence the sensor readings, especially when the sensor is oriented perpendicular to the row. In order to obtain maximum response from the sensor, the sensor head should be oriented in-line with and directly above the rows. Although orienting the light beam perpendicular to the rows is more likely to pick up background soil reflectance, this may be acceptable, however, when cotton canopy is closed with lush vegetation. In addition, operators should follow manufacturer's recommendations on an operating height range of 81-122 cm and orienting the sensor head in-line with the target to obtain maximum signal response. It is important to charge the sensor battery before use or it should be kept plugged in to avoid troubleshooting. Low battery level is likely to produce erroneous readings.
The authors have nothing to disclose.
We appreciate the assistance of Chris Parker who daily scanned the plants and Curtis Hubbard who maintained the plants in the greenhouse.
GreenSeeker | Trimble Ag. Division | Model 505 | Red NDVI sensor |
Westminster, CO | |||
Pinto beans | Producer's Co-op., Bryan, TX | Not applicable | Free choice item |
Deltapine cotton seeds | Brazos Bottom Crop Care, Caldwell, TX77836 | Not applicable | 436 RR; NonBt & RoundUp |
Ready | |||
Plastic trays | BWI, Schulenberg, TX | FG1020NL7 | 56 x 28 cm |
Label sticks | Gempler's, Janesville, WI 53547 | Item # 151276 | Durable spike-style pot markers |
4-wheel Garden push cart | Farm Tek, Dyersville, IA 52040 | Item # 108676 | 61 x 122 cm |