versão impressa ISSN 0103-9016
Sci. agric. (Piracicaba, Braz.) v.63 n.4 Piracicaba jul./ago. 2006
Distribuição espacial de Pseudaletia sequax Franclemont na cultura do triticale, em plantio direto
Célia Regina GregoI, *; Sidney Rosa VieiraI; André Luiz LourençãoII
IAPTA/IAC - Centro de Pesquisa e Desenvolvimento de Solos e Recursos Ambientais, C.P. 28 - 13020-902 - Campinas, SP - Brasil
IIAPTA/IAC - Centro de Fitossanidade
The knowledge of the spatial distribution of insect populations in crops allows the estimation of the population density in space, the classification of damage impact on yield as well as decisions to be made about control measures. To evaluate the spatial distribution of the wheat armyworm Pseudaletia sequax Franclemont, in relation to the development and yield of triticale (Triticum secale Wittmack) under no-till management, samplings were taken at 302 points in a square grid measuring 10 x 10 m, in Campinas, SP, Brazil. The caterpillar population was evaluated by counting their number per meter at the base of triticale plants. Grain yield (kg ha-1), cover of the above-ground part (kg ha-1), mean plant height (m), and ground cover percentage were evaluated after the triticale crop was mechanically harvested, with the measurement of the mass of straw left on soil surface. Geostatistics was used to analyze spatial variability, by means of semivariogram analysis, and interpolation of data by both ordinary and indicator kriging, considering the injury threshold level of ten caterpillars per square meter to construct isoline maps. Spatial dependence occurred for all variables, demonstrating an aggregate spatial distribution pattern of the caterpillar in triticale. Random sampling in this area would fail to reveal that grain yield was affected at the location with the greatest occurrence of the caterpillar. At that location, the economic injury level was exceeded, which would justify the adoption of control measures against the pest.
Key words: Triticum secale Wittmack, Insecta, wheat armyworm, geostatistics, indicator kriging
O conhecimento da distribuição espacial de populações de insetos de uma cultura permite estimar a densidade populacional no espaço, classificar seus danos causados à produção e tomar uma decisão sobre a medida de controle. Assim, este trabalho foi desenvolvido com o objetivo de avaliar a distribuição espacial da lagarta do trigo Pseudaletia sequax Franclemont, em relação ao desenvolvimento e produção do triticale (Triticum secale Wittmack) em plantio direto, Campinas, SP. Foram realizadas amostragens em 302 pontos num reticulado quadrado de 10 x 10 m. Foi avaliada a população de lagartas pela contagem do número de lagartas por metro na base das plantas de triticale. Avaliou-se ainda a produção de grãos (kg ha-1), a quantidade de cobertura da parte aérea (kg ha-1), a altura média (m) de plantas e a porcentagem de cobertura vegetal após a colheita mecanizada, permanecendo a palha na superfície do solo. Para analisar a variabilidade espacial, utilizou-se a geoestatística, através da análise de semivariogramas, interpolação dos dados por krigagem simples e indicatriz, considerando o nível de dano de dez lagartas por metro quadrado para construção de mapas de isolinhas. Ocorreu dependência espacial para todas as variáveis, evidenciando a distribuição espacial da lagarta de forma agregada na cultura do triticale. Amostragem ao acaso nesta área falharia em revelar que a produção de grãos foi afetada no local de maior ocorrência da lagarta. Nesse local, o nível de dano econômico foi ultrapassado, o que justificaria a adoção de medidas de controle da praga.
Palavras-chave: Triticum secale Wittmack, Insecta, lagarta-do-trigo, geoestatística, krigagem indicatriz
Insects only become pests when they reach population levels that cause economic damage to the crop yield, justifying, in such situations, the application of control measures that can reduce the number of insecticide applications and preserve the natural enemy populations. Their spatial distribution, defined as the position that individuals occupy in space relative to each another, at a given focus (Vieira, 2000), influences insect population dynamics. Differences in structure can be observed in relation to the sampled area, crop where they live (Knudsen et al., 1994), species (Nestel & Klein, 1995), and developmental stage (Farias et al., 2001; Pereira et al., 2004). Insect spatial distribution studies have been performed using geostatistical analysis (Midgarden et al., 1993; Nestel & Klein, 1995; Ellsbury et al., 1998; Pereira et al., 2004) and have shown that random distribution patterns rarely occur in nature. According to Jesus et al. (2002), the aggregate type of distribution is predominant, since these patterns are influenced by specific differences in each environment. According to Reis & Miranda Filho (2003), these evidences compromise the validity and effectiveness of analysis models that assume independent errors. These studies allow population density estimates to be obtained and help control measures to be adopted more safely.
Many studies focus on the spatial distribution of insects, but few show results with regard to the development and yield of the affected crops. Among the pests that damage the aerial part of triticale (Triticum secale Wittmack), according to Gallo et al. (2002), Pseudaletia sequax Franclemont is one of the most important, due to defoliation and because of its difficult control. According to Gassen (2005), the economic injury level is ten caterpillars per square meter. The objective of this work was to evaluate the spatial distribution of P. sequax caterpillars in relation to triticale development and yield under no-till management.
MATERIAL AND METHODS
The field experiment was carried out at the "Centro Experimental" of the "Instituto Agronômico de Campinas" (IAC), SP, Brazil, and covers an area of 3.42 ha, (22º53' S and 47º04' W, at a mean elevation of 600 m above sea level) presenting a 6.5% slope. According to Oliveira et al. (1989) the soil is a clayey-textured Rhodic Eutrudox with a moderate A-horizon. According to Setzer (1966) international climate classification, the Campinas region presents a transition between the climatic types Cwa and Cfa, i.e., a tropical climate with a dry winter and a wet summer. The no-till management, system was adopted on area since 1985, with summer and winter crops for the production of grain and ground cover.
Three hundred and two points (10 x 10 m) were marked on the field for sampling in square grid (Figure 1). The triticale cultivar IAC 2 was sown as a winter crop in April 2004 under the no-till system. A spacing of 20 cm between rows was used, at a density of 60-70 seeds per meter. Fertilization at planting consisted of 150 kg ha-1 of a 8-28-16 NPK formulation, following, soil analyses and recommendations of Raij et al. (1996). The crop was sprayed 32 days after seeding to control Sitobion avenae (Fabr.), with a 30-meter range cannon, using 2 mL L-1 monocrotophos, safety interval of 21 days, at a spray volume of 0.8 L ha-1. The P. sequax population was evaluated 69 days after sowing, by counting the number of caterpillars found on the ground, around the base of the stem of plants, along one meter of row in each of the 302 points, to obtain the number of caterpillars per m2. At this time, the triticale was at the grain-filling stage.
The crop was harvested in September 2004, on a 1.5 x 1.5 m plot (2.25 m2), and grain yield (kg ha-1), dry mass yield of the above-ground part (kg ha-1), and mean height (m) were obtained for 10 plants per point. Percentage of residue ground cover was obtained after harvesting, using a 1 m2 square divided by lines into one hundred 0.1 x 0.1 m sub-squares. The square was inserted at random into each of the 302 sampling points, and ground cover was estimated at the intersection points of lines that coincided with straw on the soil surface. Data were converted to percentage of ground cover on soil surface.
Data were initially analyzed by descriptive statistics, and mean, variance, standard deviation, coefficient of variation, maximum value, minimum value, skewness, and kurtosis were obtained using the STAT software, according to Vieira et al. (2002). This analysis was carried out to detected outlier values and the frequency distribution of the data. When a set of data approaches the normal distribution, their skewness and kurtosis values approach zero; in this case, the central tendency value (mean) is equal to the measure that occurs most often (mode). If the skewness and kurtosis coefficients are greater than zero, the frequency distribution of data will tend to be lognormal, which is characterized by having a great number of small values, plus some values so large that will cause differences between the mean, the median, and the mode. This happens when the variable has a very localized distribution, i.e., there are sites with high values but the vast majority has low values.
To verify the spatial variability of variables, data were analyzed by semivariograms g(h), described by Vieira (2000), based on the assumptions of stationarity of the intrinsic hypothesis:
where N(h) is the number of pairs of measured values Z(xi), Z(xi+h), separated by a vector h, which is the distance determined from the Z(xi) and Z(xi + h) coordinates. Calculation of equation (1) generates g(h) values corresponding to h distances, for the construction of the semivariogram. According to Vieira (2000), it is expected that measurements located near one another are more similar among themselves than those separated by great distances, i.e., where g(h) increases with the h distance until a maximum value at which it stabilizes, at a level corresponding to the limit distance of spatial dependence, which is the range. Measurements located at greater distances than the range have random distribution and are therefore independent among themselves; beyond this distance, classic statistics can be applied.
Mathematical models were fitted to the semivariograms, allowing the nature of spatial variation to be displayed for any distances within the measured field. The criteria and procedures for fitting the semivariogram models were made according to Vieira et al. (1983). Based on the model used to fit the data, the following semivariogram parameters were defined: a) nugget effect (C0), which is the g value when h = 0; b) range of the spatial dependence (a), which is the distance at which g(h) remains approximately constant, after increasing as h increases constant; c) threshold (C0+C1) which is the g(h) value beyond the range approaching the data variance, if it exists.
Spatial dependence can be expressed by the DD parameter (degree of spatial dependence), which is the proportion of structural variance (C1) in relation to the threshold (C0 + C1), and can be calculated by the equation:
According to Zimback (2001), this can be used to classify the spatial dependence into weak if DD< 25%; moderate for DD between 26% and 75%; and strong if DD> 75%.
When semivariance is dependent upon distance, i.e., when spatial autocorrelation between samples exists, values can be interpolated in the study field, without bias and with minimum variance, using the kriging method described by Vieira (2000). Kriging is a technique used in geostatistics to estimate property values for locations where the property has not been measured. For this tool to be used it is required that the spatial dependence defined by the semivariogram exists (Vieira, 2000). Therefore, kriging is a process by which spatially-distributed property values are estimated from interdependent adjacent values (Vieira et al., 1983). The estimated value for the property, Z*(X0), is given by the expression:
where N is the number of neighbors measured, Z(xi) are used to estimate the property and li are the weights applied to each Z(xi), which are selected so that the estimate obtained is unbiased. This implies assuming two conditions:
Values interpolated by kriging were used to construct isoline maps for the variables under study as a function of their coordinates, using the Surfer program (Golden Software, 1999). The analysis of maps obtained via the kriging-estimated values allows the identification of sites with higher or lower occurrence of caterpillars or triticale grain yield.
The spatial variability in the number of caterpillars was also analyzed by indicator kriging (Goovaerts, 1997). Since the economic injury level is 10 caterpillars per square meter, levels of zero and one were established for values smaller and greater than ten, respectively, and values were interpolated and a map was constructed for caterpillar injury level in triticale, indicating the locations where control would be economically viable.
RESULTS AND DISCUSSION
Mean, variance, standard deviation, coefficient of variation (C.V.), minimum value, maximum value, skewness, and kurtosis (Table 1) indicated high variability for the caterpillar population, with a coefficient of variation in the order of 100%. In terms of plant height, grain yield, and ground cover, low coefficients of variation were obtained, smaller than 30%, according to criteria adopted by Warrick & Nielsen (1980). Ground cover showed low coefficients of variation, 10.26%, after harvesting the triticale, which could be attributed to the mass of straw distribution on the soil surface imposed by the harvesting machine.
The triticale crop grew under a very favorable condition in terms of water availability (Figure 2) and edaphic conditions, the main reason for the high grain (2,660 kg ha-1) and dry mass (average 3,243 kg ha-1) mean yields obtained, even when mean coefficient of variation (34.72%) and the attack by P. sequax are taken into consideration (Table 1). The variables grain yield and dry mass yield showed normal distribution with skewness and kurtosis values near zero, except for plant height, which may indicate a lognormal distribution.
Semivariograms (Figure 3) and their fitted parameters (Table 2) indicate that the spatial dependence can be considered weak to moderate, with a high nugget effect, near the structural variance and with long range. The nugget effect indicates the variability at distances smaller than the sampling spacing (10 m) whose value include variability due to sampling errors and to point to point variation. Its value also reveals the amount of population aggregation in the spatial distribution with respect to sampling interval. The nugget effect values found are comparable to those presented by Ellsbury et al. (1998) for Diabrotica spp. in corn roots.
Height in triticale plants showed a smaller degree of spatial dependence (Table 2) and also a very low coefficient of variation (5.8%) according to Table 1. This indicates that plant height had little spatial variation, i.e., it was extremely uniform in terms of height, and therefore it was not affected by the caterpillar, which caused defoliation when the plants probably had already reached their maximum growth. The obtained moderate degree of spatial dependence obtained justifies the sampling intensity adopted and is in agreement with results of insect spatial distribution studies (Farias et al., 2001; 2004, Jesus et al., 2002; Reis & Miranda Filho, 2003), since occurrence of the caterpillar was of the aggregate form, not at random. It can also be observed that it was necessary to remove the linear tendency for a better characterization of the spatial variability to fit a model to the semivariogram of dry mass (Figure 3e). However, dry mass ranges on one direction can be perceived on the map as a result of that tendency (Figure 4).
Grain yield was visually reduced in places where there was greater caterpillar attack (Figure 4). In this case, an indicator kriging can map, with greater precision, regions of the field where chemical control should be necessary based on the injury level of 10 caterpillars per square meter. Notwithstanding, the attack was detrimental to grain yield in areas where the number of caterpillars exceeded the economic injury level.
Reported results can not be extrapolated to other fields and also for the description of the insect distribution in other conditions, as this was not the objective of the study. However, the ranges of the semivariograms for insect population (95 m as shown in figure 3a) and for the grain yield (80m as shown in figure 3d) indicate that the sampling interval of 10 m was adequate to characterize the insect population and the damage caused on grain yield. Therefore, the sampling distance used can adopted in other fields having similar conditions.
To CNPq for a research fellowship granted to the first author and for research productivity fellowships granted to the other authors, and to FAPESP for financial support to this project.
ELLSBURY, M.M.; WOODSON, W.D.; CLAY, S.A.; MALO, D.; SCHUMACHER, J.CLAY,D.E.; CARLSON, C.G. Geostatistical characterization of the spatial distribution of adult corn rootworm (Coleoptera: Chrysomelidae) emergence. Environmental Entomology, v.27, p.910-917, 1998. [ Links ]
FARIAS, P.R.S.; BARBOSA, J.C.; BUSOLI, A.C. Distribuição espacial da lagarta-do-cartucho, Spodoptera frugiperda (J.E. Smith) (Lepidoptera: Noctuidae), na cultura do milho. Neotropical Entomology, v.30, p.681-689, 2001. [ Links ]
FARIAS, P.R.S.; ROBERTO, S.R.; LOPES, J.R.S.; PERECIN, D. Geostatitical characterization of the spatial distribution of Xylella fastidiosa sharpshooter vectors on citrus. Neotropical Entomology, v.33, p.13-20, 2004. [ Links ]
GALLO, D.; NAKANO, O.; SILVEIRA NETO, S.; CARVALHO, R.P.L.; BATISTA, G.C.; BERTI FILHO, E.; PARRA, J.R.P.; ZUCCHI, R.A.; ALVES, S.B.; VENDRAMIN, J.D.; MARCHINI, L.C.; LOPES, J.R.S.; OMOTO, C. Entomologia agrícola. Piracicaba: FEALQ, 2002. 920p. [ Links ]
GOLDEN SOFTWARE. Surfer 7.0: Contouring and 3D surface mapping for scientist's engineers. User's Guide. New York, 1999. 619p. [ Links ]
GOOVAERTS, P. Geostatistics for natural resources evaluation. New York: Oxford University Press, 1997. 483p. (Applied Geostatistics Series). [ Links ]
JESUS, C.R.; REDAELLI, L.R.; ROMANOWSKI, H.P. Population dynamics and spatial distribution of Spartocera dentiventris (Berg) (Hemiptera: Coreidae) adults on Nicotiana tabacum L. (Solanaceae). Neotropical Entomology, v.31, p.541-549, 2002. [ Links ]
KNUDSEN, G.R.; SCHOTZKO, D.J.; KRAG, C.R. Fungal entomopathogen effect on numbers and spatial patterns of the Russian wheat aphid (Homoptera: Aphididae) on preferred and nonpreferred host plans. Envirommental Entomology, v.23, p.1558-1567, 1994. [ Links ]
MIDGARDEN, D.G.; YOUNGMAN, R.R.; FLEISCHER, S.J. Spatial analysis of counts of Western corn rootworm (Coleoptera: Chysomelidae) adults on yellow sticky traps in corn. Geostatistics and Dispersion Indices, v.22, p.1124-1133, 1993. [ Links ]
NESTEL, D.; KLEIN, M. Geostitistical analysis of leafhopper (Homoptera: Cicadellidae) colonization and spread in deciduous orchards. Environmental Entomology, v.24, p.1032-1039, 1995. [ Links ]
OLIVEIRA, M.N.; OLIVEIRA, J.B.; LEPSCH, I.F. Mapa por compilação de solos do estado de São Paulo. Rio de Janeiro: EMBRAPA, SNLCS; IAC, 1989. [ Links ]
PEREIRA, M.F.A.; BOIÇA JR., A L.; BARBOSA J.C. Distribuição espacial de Bemisia tabaci (Genn.) biótipo B (Hemiptera: Aleyrodidae) em feijoeiro (Phaseolus vulgaris L.) Neotropical Entomology, v.33, p.493-498, 2004. [ Links ]
RAIJ, B. van; CANTARELLA, H.; QUAGGIO, J.A.; FURLANI, A.M.C. Recomendações de adubação e calagem para o estado de São Paulo. 2.ed. Campinas: Instituto Agronômico; Fundação IAC, 1996. 285p. (Boletim Técnico, 100). [ Links ]
REIS, A J.S.; MIRANDA FILHO, J.B. Autocorrelação espacial na avaliação de compostos de milho para resistência da lagarta do cartucho (Spodoptera frugiperda). Pesquisa Agropecuária Tropical, v.33, p.65-72, 2003. [ Links ]
SETZER, J. Atlas climático e ecológico do Estado de São Paulo. São Paulo: Comissão Interestadual da Bacia Paraná-Uruguai CESP, 1966. 57p. [ Links ]
VIEIRA, S.R. Geoestatística em estudos de variabilidade espacial do solo. In: NOVAIS, R.F.; ALVAREZ, V.H.; SCHAEFER, G.R. (Ed.) Tópicos em ciência do solo. Viçosa: SBCS, 2000. v.1, p.1-54. [ Links ]
VIEIRA, S.R.; MILLETE, J.; TOPP, G.C.; REYNOLDS, W.D. Handbook for geostatistical analysis of variability in soil and climate data. In: ALVAREZ, V.V.H.; SCHAEFER, C.E.G.R.; BARROS, N.F.; MELLO, J.W.V.; COSTA, L.M. (Ed.) Tópicos em ciência do solo. Viçosa: SBCS, 2002. v.2, p.1-45. [ Links ]
VIEIRA, S.R.; HATFIELD, J.L.; NIELSEN, D.R.; BIGGAR, J.W. Geostatistical theory and application to variability of some agronomical properties. Hilgardia, v.51, p.1-75, 1983. [ Links ]
WARRICK, A.W.; NIELSEN, D.R. Spatial variability of soil physical properties in the field. In: HILLEL, D. (Ed). Applications of soil physics. New York: Academic Press, 1980. p.319-44. [ Links ]
ZIMBACK, C.R.L. Análise especial de atributos químicos de solo para o mapeamento da fertilidade do solo. Botucatu, UNESP, 2001. 114p. (Tese - Livre Docência). [ Links ]
Received December 20, 2005
Accepted June 18, 2006s