Resumos

SEÇÃO VI- MANEJO E CONSERVAÇÃO DO SOLO E DA ÁGUA

Variabilidade espacial de atributos do solo para adoção do sistema de agricultura de precisão na cultura de cana-de-açúcar¹ 1 Parte da Tese de Mestrado do segundo autor, apresentado ao Curso de Pós-Graduação em Agronomia, Ciência do Solo, Departamento de Solos e Adubos, Faculdade de Ciências Agrárias e Veterinárias — UNESP/Campus de Jaboticabal (SP).

Assessment of spatial variability of soil attributes as a basis for the adoption of precision agriculture in sugarcane plantations

J. E. Corá^I; A. V. Araujo^II; G. T. Pereira^III; J. M. G. Beraldo^IV

^IProfessor do Departamento de Solos e Adubos da Faculdade de Ciências Agrárias e Veterinárias — UNESP. Via de Acesso Prof. Paulo Donato Castellane s/n, CEP 14884-900 Jaboticabal (SP). E-mail: cora@fcav.unesp.br

^IIProfessora da Universidade Estadual do Mato Grosso do Sul — UEMT. Rodovia MS 306 Km 6,4, CEP 79540-000 Cassilândia (MS). E-mail: anamari.v@uol.com.br

^IIIProfessor do Departamento de Ciências Exatas da Faculdade de Ciências Agrárias e Veterinárias, UNESP. E-mail: generta@fcav.unesp.br

^IVMestrando do Curso de Pós-Graduação em Produção Vegetal da Faculdade de Ciências Agrárias e Veterinárias, UNESP. Bolsista da CAPES. E-mail: jberaldo@terra.com.br

RESUMO

O objetivo do presente trabalho foi avaliar a variabilidade espacial de atributos de Latossolos sob cultivo de cana-de-açúcar na região de Jaboticabal (SP). Foram feitas amostragens do solo a intervalos regulares de 50 m, em uma área de 90 ha, nas profundidades de 0,00—0,20 e 0,60—0,80 m para determinação de pH, CTC, V % e teores de matéria orgânica, P, K⁺, Ca²⁺, Mg²⁺, H + Al e argila. Os dados foram submetidos às análises: estatística descritiva, geoestatística e interpolação por krigagem. Os alcances de dependência espacial para os atributos químicos do solo e teores de argila na camada de 0,60—0,80 m de profundidade foram menores, quando comparados àqueles referentes à camada de 0,00—0,20 m. Estes resultados demonstraram maior descontinuidade na distribuição espacial dos atributos do solo na camada de 0,60—0,80 m de profundidade dos Latossolos, indicando que essa classe de solos não apresentou homogeneidade de seus atributos como conceituadamente a ela foi atribuída. O manejo no solo alterou a dependência espacial dos atributos do solo na camada superficial de forma a diminuir a variabilidade espacial dos atributos químicos do solo em relação à camada mais profunda. A investigação da variabilidade espacial de atributos químicos e do teor de argila da camada superficial e subsuperficial dos solos proporcionou condições para a definição de zonas homogêneas de manejo, o que permite a adoção do sistema de agricultura de precisão.

Termos de indexação: manejo do solo, geoestatística, krigagem, argila do solo.

SUMMARY

The spatial variability of physical and chemical attributes of Latosols (Oxisols) cultivated with sugarcane was assessed. The study was carried out in a 90 ha area in the surroundings of Jaboticabal, São Paulo State, Brazil (21 º 15 ' S and 48 º 18 ' W). A 50 m regular grid was projected onto the experimental area. Soil samples were obtained at each 50 m grid point (421 points) from the 0 to 0.20 m and from the 0.60—0.80 m soil layers. Organic matter, pH, P, K⁺, Ca²⁺, Mg²⁺, H + Al, CEC, soil clay content, and V% values of each sample were determined. The data set was submitted to descriptive statistics, geostatistics, and kriging analyses. The range of the spatial dependence for the soil attributes, was lower in the 0.60—0.80 m than in the 0—0.20 m layer. Results showed a higher discontinuity in the spatial variability of soil attributes in the 0.60—0.80 m layer of Oxisols indicating that, in disagreement to the currently established concept, this soil class does not present homogeneity of soil attributes. Soil management altered the spatial dependence of soil chemical attributes in the surface, reducing the spatial variability structure in comparison to the deeper soil layer. The assessment of the spatial variability of soil chemical attributes and soil clay content in the top and the subsurface layers allowed the identification of homogeneous site-specific management zones which permit the adoption of precision agriculture for sugarcane.

Index terms: soil management, geostatistics, kriging, soil clay.

INTRODUÇÃO

O manejo regionalizado do solo e da cultura é parte integrante de um sistema de Agricultura de Precisão, o qual envolve conceitos de uso de informações sobre a variabilidade de propriedades locais e climáticas de uma área, visando ao aumento da produtividade, otimização no uso dos recursos e redução do impacto da agricultura ao meio ambiente. Os processos e os atributos do solo que determinam o desempenho e a produção das culturas, bem como o impacto da agricultura ao meio ambiente, variam no espaço e no tempo. Por essa razão, o conhecimento da variabilidade espacial e temporal dos fatores de produção da cultura é o primeiro passo para adoção, com êxito, do sistema de agricultura de precisão. Este sistema tem sido adotado por diferentes grupos, incluindo fabricantes de equipamentos, fornecedores de insumos, companhias de sementes, consultores agronômicos, cientistas e produtores, como um caminho para melhorar o retorno econômico da atividade agrícola (Runge & Hons, 1999).

A variabilidade do solo é conseqüência de complexas interações dos fatores e processos de sua formação. Além dos fatores e processos, práticas de manejo do solo e da cultura são causas adicionais de variabilidade (Corá, 1997). Áreas pedologicamente idênticas podem apresentar variabilidade distinta em atributos, quando submetidas às diferentes práticas de manejo. Da mesma forma, áreas pedologicamente diferentes, quando submetidas ao mesmo manejo, podem apresentar-se semelhantes em seus atributos. O manejo pode alterar atributos químicos, físicos, mineralógicos e biológicos, com impacto principalmente nas camadas superficiais do solo (Blevins et al., 1983; Seta et al., 1993; Corá, 1997).

O conhecimento da distribuição espacial dos atributos do solo em determinada área é importante para o refinamento das práticas de manejo e avaliação dos efeitos da agricultura sobre a qualidade ambiental (Cambardella et al., 1994), assim como é importante para a definição da intensidade de amostragem do solo para sua caracterização, possibilitando, dessa maneira, reduzir o erro-padrão da média, maximizando a eficiência da amostragem, e reduzir os custos e mão-de-obra do trabalho.

A variabilidade dos atributos do solo é freqüentemente descrita pela estatística clássica, a qual assume que a média de determinada variável em uma unidade amostral é a esperança do valor da variável em toda a unidade, com a estimativa do erro expressa pela variância. Esta abordagem assume que a variabilidade sobre a média é devida ao acaso e não contém referência sobre a distribuição no espaço de diferenças dentro de unidades amostrais. Entretanto, estudos têm mostrado que atributos do solo, geralmente, não ocorrem de maneira aleatória no espaço, mas segundo um arranjo estrutural com uma dimensão característica, que é seu domínio, e que corresponde à distância dentro da qual há interdependência dos valores medidos. Ou seja, os atributos do solo podem apresentar-se espacialmente correlacionados em determinada área, podendo o conhecimento desta dependência espacial contribuir para o entendimento sobre a influência dos atributos e do manejo do solo na produtividade das culturas. Todavia, os métodos tradicionais de classificação do solo e análises estatísticas não consideram esta questão (Trangmar et al., 1985).

Portanto, conhecer a variabilidade espacial de atributos do solo que controlam a produtividade das culturas e os riscos de contaminação do ambiente e investigar as causas dessa variabilidade são fatores importantes em um sistema de produção que vise sustentabilidade por meio do manejo regionalizado de insumos e práticas agrícolas, como é o sistema de agricultura de precisão.

O objetivo deste trabalho foi caracterizar a variabilidade espacial de atributos químicos e teor de argila dos horizontes: superficial e subsuperficial de Latossolos em áreas de intenso cultivo com cana-de-açúcar.

MATERIAL E MÉTODOS

A área deste estudo localiza-se na região de Jaboticabal/SP, cujas coordenadas são 21 º 15 ' de latitude sul e 48 º 18 ' de longitude oeste, com altitude média de 670 m. O relevo é predominantemente suave ondulado, com declividades médias variando de 3 a 8 %. A área experimental está sob cultivo de cana-de-açúcar há mais de trinta anos, onde os solos foram classificados como Latossolo Vermelho distrófico típico textura média e Latossolo Vermelho eutroférrico típico argiloso epidistrófico, com predomínio deste último, sendo considerados homogêneos do ponto de vista pedológico e de manejo.

Uma malha de amostragem com espaçamento regular de 50 m entre pontos foi estabelecida em uma subárea com dimensões de 300 x 3.000 m (90 ha), totalizando 421 pontos de amostragem. Em cada ponto, foram coletadas, com a utilização de um trado tipo holandês, cinco subamostras de material de solo para compor uma amostra composta representativa do ponto de amostragem, sendo uma no centro e as outras nos quatro pontos cardeais espaçadas 2 m do ponto central. As amostragens foram feitas nas profundidades de 0,00—0,20 e 0,60—0,80 m.

Foram determinados nas amostras de solo o pH, teores de matéria orgânica (MO), P, K⁺, Ca²⁺, Mg²⁺, H + Al e, posteriormente, calculadas a CTC e a saturação por bases (V %), segundo procedimentos descritos por Raij et al. (1987). O teor de argila foi determinado pelo método da pipeta (Gee & Bauder, 1986).

As análises estatísticas descritivas foram realizadas pelo programa SAS (1995). A hipótese de normalidade dos dados foi testada pelo teste W a 5 % (Shapiro & Wilk, 1965). Os dados que não apresentaram distribuição normal foram transformados em sua forma logarítmica (y = lnx) e novamente submetidos ao teste W. Após transformação, os dados apresentaram distribuição normal.

A análise de dependência espacial das variáveis foi efetuada pelo ajuste dos dados ao semivariograma experimental, de acordo com a teoria das variáveis regionalizadas (TVR), utilizando o programa GS+ (Robertson, 1998). Os cálculos de semivariância foram baseados numa distância máxima de 1.100 m, na profundidade de 0,00—0,20 m, e de 700 m, na profundidade de 0,60—0,80 m, com intervalos médios de 50 m.

A distância limite de dependência espacial, denominada alcance (A₀), indica que amostras localizadas a distâncias menores que o alcance são correlacionadas umas com as outras. O alcance depende do tamanho da área amostrada e da escala de observação, sendo tanto maior quanto maior for o intervalo entre medidas (Trangmar et al., 1985). Os modelos exponenciais não exibem um valor de patamar finito, mas, para propostas práticas, há um ponto limite no qual a semivariância pára de aumentar (Webster, 1985). Neste estudo, foi adotado o alcance prático, correspondente a 95 % do patamar para discussão dos semivariogramas ajustados ao modelo exponencial.

A seleção dos modelos foi realizada com base no melhor coeficiente de determinação (R²) e menor soma de quadrados do resíduo (SQR). Os modelos dos semivariogramas ajustados foram validados pela técnica de validação cruzada (Cressie, 1991). Na elaboração dos mapas de distribuição espacial das variáveis, foi utilizado o programa Surfer 6.01 (Golden Software, 1995), com base nos valores estimados por krigagem, realizada com o programa GS+ (Robertson, 1998)

RESULTADOS E DISCUSSÃO

Os valores da média e da mediana para todas as variáveis foram similares, com as medianas apresentando valores um pouco menores que a média, apesar da ocorrência de algumas distribuições assimétricas (Quadro 1). Segundo Cambardella et al. (1994), isto pode ser um indicativo de que as medidas de tendência central não são dominadas por valores atípicos na distribuição.

Observa-se que, com exceção dos valores de CTC e teores de K⁺, na camada de 0,00—0,20 m, e teores de argila, na camada de 0,60—0,80 m, os demais atributos não apresentaram resultados significativos (P < 0,05) no teste de normalidade Shapiro & Wilk (1965) (Quadro 1). A mesma situação foi observada por Cambardella et al. (1994), onde a log-transformação dos dados reduziu a assimetria, porém a distribuição ainda apresentou curtose, indicando maior falha na redução da curtose do que na assimetria. Segundo Schlotzhauer & Littell (1997), a ocorrência de valores de assimetria positivos nos atributos que seguem a distribuição lognormal é indicativo da maior freqüência de valores menores do que a média e poucos valores maiores do que ela, confirmando os resultados encontrados por Souza et al. (1998), segundo os quais, para áreas adubadas, isto pode ser o resultado da aplicação não-uniforme de adubos. Assim, ocorrem no campo alguns locais com menores teores de nutrientes e outros com maiores. Oliveira et al. (1999) encontraram distribuição normal para CTC e distribuição lognormal para K⁺ e Mg²⁺, não tendo o pH seguido nenhuma distribuição avaliada pelos autores.

De acordo com os limites estabelecidos para os níveis de fertilidade do solo, para a cultura da cana-de açúcar no estado de São Paulo (Raij et al., 1996), os valores médios dos atributos químicos na camada de 0,00—0,20 m apresentaram-se baixos para pH e V %; médios para P e MO e altos para Ca²⁺, Mg²⁺ e K⁺. Observou-se que, apesar de os teores das bases trocáveis (Ca²⁺, Mg²⁺ e K⁺) terem sido classificados como altos, o pH do solo e a V % mantiveram-se baixos, indicando que a maior parte do complexo de troca do solo estava sendo ocupado pelos cátions ácidos: H⁺ e Al³⁺.

Embora, em geral, a média dos valores dos atributos para a área estudada tenha permanecido nas faixas de fertilidade média ou alta, na profundidade de 0,00—0,20 m, os valores apresentaram grande amplitude. O P apresentou teores que variaram desde 7 até 62 mg dm^-3, isto é, verificou-se que, em determinados locais, os teores de P foram considerados baixos, em outros médios e em outros altos. O mesmo foi observado para os demais atributos. Essa grande amplitude de valores revelou os problemas que podem ocorrer quando se usa a média dos valores dos atributos como base para a tomada de decisão sobre a realização do manejo químico do solo, ou seja, em alguns locais, a dose recomendada de fertilizantes ou corretivos estará subdimensionada; em outros, será adequada e, em outros, poderá haver aplicação excessiva de determinado fertilizante ou corretivo. Isto resulta em prejuízos econômicos, tanto pela aplicação desnecessária, quanto pelo desequilíbrio entre as quantidades de nutrientes que serão disponibilizados para as plantas, além do efeito deletério que essa prática pode ocasionar sobre o meio ambiente.

Em geral, os valores médios dos atributos analisados na profundidade de 0,60—0,80 m foram menores. Os valores mais elevados dos teores de P, K⁺, e Mg²⁺ na camada superficial (0,00—0,20 m) podem ser atribuídos às frequentes aplicações de fertilizantes e calcário ao longo do tempo nesta camada. A MO, em maiores teores nesta camada, proporcionou maiores valores da CTC do solo. Já os valores de pH apresentaram-se mais baixos na camada superficial, provavelmente por causa da acidificação natural que ocorre no solo ao longo do tempo e pela provocada por adubações nitrogenadas. Observou-se maior amplitude de valores dos atributos do solo na camada superficial, indicando que o manejo do solo pode ser causa adicional da variabilidade, como constatado por Corá (1997).

Seguindo a classificação do CV, proposta por Wilding & Drees (1983), observou-se que os valores de pH e CTC e teores de MO e de argila na camada de 0,00—0,20 m apresentaram baixa variabilidade, ou seja, CV < 15 %. Freqüentemente, são encontrados na literatura estudos nos quais os teores de MO apresentam moderada variabilidade, isto é, 15 < CV < 35 % (Cambardella et al., 1994; Souza et al., 1997; Salviano et al., 1998).

Neste trabalho, os valores de H + Al e V % apresentaram moderada variabilidade na profundidade de 0,00—0,20 m. Souza et al. (1997), trabalhando com Latossolo Amarelo distrófico de textura média, sob pomar de citros, encontraram resultados semelhantes para H + Al. Entretanto, Salviano et al. (1998) observaram elevada variabilidade (CV > 35 %) para os valores de H + Al e V % em um Argissolo Vermelho-Amarelo. No presente estudo, os teores de P, K⁺, Ca²⁺ e Mg²⁺ apresentaram elevada variabilidade, resultados semelhantes àqueles encontrados por Salviano et al. (1998). Chien et al. (1997), entretanto, obtiveram valores de CV elevados para os teores de Ca²⁺ e Mg²⁺ e extremamente altos para os de P (199 %), tendo esta alta variabilidade sido atribuída aos efeitos da aplicação de fertilizantes ao solo, segundo os autores. Oliveira et al. (1999), trabalhando com solo salino-sódico de origem aluvial, observaram, na camada de 0,30 m, menor CV para pH (11,06 %) e maior para Mg²⁺ (95,68 %), apresentando os demais atributos químicos (K⁺, Ca²⁺, Na⁺, CTC, percentagem de Na⁺ (PST) e condutividade elétrica (CE)) valores intermediários a esses.

Em geral, observaram-se menores valores de CV na camada de 0,60—0,80 m de profundidade, exceto para os teores de P e valores de CTC. Os teores de MO e de argila e valores de pH e V % apresentaram CV baixo; os teores de Mg²⁺ e valores de CTC apresentaram CV médio e os teores de P, K⁺ e Ca²⁺ apresentaram valores de CV altos. Segundo Wollenhaupt et al. (1997), mesmo que os valores do CV sejam moderados, este não é necessariamente um bom indicador da variabilidade espacial dos atributos do solo, haja vista a ocorrência de locais no campo com valores extremamente altos ou baixos.

Na seleção dos modelos dos semivariogramas, foram considerados os valores de R² e SQR (Quadro 2). Segundo Robertson (1998), a SQR é um parâmetro mais robusto do que o R² e propicia uma medida exata de quão bem o modelo se ajusta aos dados e o programa GS⁺ utiliza a combinação de parâmetros do semivariograma para minimizar a SQR para cada modelo.

Todos os atributos do solo analisados apresentaram dependência espacial nas duas profundidades. Na camada de 0,00—0,20 m, os teores de P, Ca²⁺, e Mg²⁺ e valores de pH e V % foram ajustados ao modelo esférico, com coeficiente de determinação do modelo (R²) acima de 0,95. Os dados referentes aos teores de MO, K⁺ e argila e valores de H + Al e CTC foram ajustados ao modelo exponencial, com R² acima de 0,90.

Na camada de 0,60—0,80 m, os dados se ajustaram, em sua maioria, ao modelo exponencial, com exceção dos teores de Mg²⁺ e de argila, que se ajustaram ao modelo esférico. Os dados referentes aos valores de H + Al e pH e teores de K⁺ ajustaram-se aos seus respectivos modelos com menores valores de R², ou seja, 0,592, 0,804 e 0,806, respectivamente, enquanto, para os demais atributos, os valores de R² foram superiores ou iguais a 0,97.

Os atributos do solo apresentaram diferentes alcances de dependência espacial, e, na camada superficial, os teores de P e Mg²⁺ apresentaram os maiores alcances, ou seja, 1.034 e 1.001 m, respectivamente, enquanto os teores de K⁺ e MO apresentaram os menores alcances, 102 e 126 m, respectivamente. A menor continuidade espacial dos teores de K⁺ pode ser explicada pela aplicação não-uniforme de vinhaça, efetuada na área por meio de irrigação por aspersão utilizando sistema autopropelido com canhão hidráulico, cuja uniformidade de aplicação de água é bastante baixa.

Os alcances de dependência espacial para os atributos do solo na profundidade de 0,60—0,80 m foram, com exceção dos teores de argila, bem menores, quando comparados àqueles referentes à profundidade de 0,00—0,20 m. O menor valor foi encontrado para os teores de K⁺ (28,7 m) e o maior para os teores de Mg²⁺ (681 m). Estes resultados demonstram maior descontinuidade na distribuição espacial dos atributos do solo na camada de 0,60—0,80 m de profundidade. Como o estudo foi realizado em área sob Latossolos, os quais apresentam como característica intrínseca a homogeneidade dos atributos do solo ao longo do perfil, esperava-se menor descontinuidade espacial dos atributos do solo na camada mais profunda, tanto pelo fato de os Latossolos se apresentarem homogêneos ao longo do perfil, como pelo fato de que a camada mais profunda está menos sujeita aos efeitos do manejo do solo realizado na superfície, mantendo, dessa maneira, suas características originais de homogeneidade. Contudo, o manejo no solo, realizado na área ao longo do tempo, tais como: aplicações de fertilizantes e calcário e preparo convencional com arações e gradagens, pode ter contribuído para maior homogeneização da camada superficial, tendo, como conseqüência, o aumento do alcance de dependência espacial dos atributos do solo, o que condicionou a maior continuidade na distribuição espacial dos mesmos nesta camada em relação à camada mais profunda. Berg & Oliveira (2000) observaram que o manejo dos solos com cana-de-açúcar provocou alterações até camadas mais profundas em relação à rotação soja/trigo. Segundo os autores, as camadas mais superficiais do solo sob a rotação soja/trigo apresentaram pouca variação como resultado do manejo muito intensivo.

De qualquer forma, é importante ressaltar, no presente estudo, a ocorrência de descontinuidade espacial dos atributos do solo nas camadas mais profundas dos Latossolos da área, a qual pode ser explicada baseando-se em duas hipóteses: a primeira é que os Latossolos não apresentam a homogeneidade vertical e horizontal dos atributos analisados ao longo do perfil como se imaginava anteriormente. Tal fato pôde ser constatado pelas novas técnicas de investigação da variabilidade espacial dos atributos do solo, como, por exemplo, o uso da geoestatística. A segunda é que o manejo do solo, realizado na superfície, estaria alterando a variabilidade espacial dos atributos químicos do solo ao longo do perfil.

Entretanto, ao analisar os dados (Quadros 1 e 2), observa-se que mesmo aqueles atributos do solo pouco influenciáveis pelo manejo, na profundidade de 0,60—0,80 m, como CTC e teor de argila, apresentaram amplitudes elevadas e valores de alcance de dependência espacial relativamente baixos, comprovando a baixa continuidade espacial desses atributos na camada mais profunda do solo. Da mesma maneira, analisando os demais atributos (Quadro 1), não é possível observar valores incompatíveis àqueles comuns em horizontes diagnósticos de subsuperfície de Latossolos, que poderiam indicar efeito do manejo realizado na superfície. Portanto, a segunda hipótese poderia ser descartada.

Segundo Trangmar et al. (1985), o alcance define o raio máximo para o qual amostras vizinhas são usadas para interpolação por técnicas de krigagem. Assim, baixos valores de alcance podem influir na qualidade das estimativas, uma vez que poucos pontos são usados para realização da interpolação. Portanto, as estimativas obtidas por interpolação por krigagem para os teores de MO e K⁺ e valores de H + Al, na profundidade de 0,00—0,20 m e, para os teores de MO e K⁺ e para os valores de pH e H + Al, na profundidade de 0,60—0,80 m, podem apresentar-se pouco confiáveis para fins práticos, em razão dos baixos valores de alcance de dependência espacial apresentados por estes atributos.

O efeito pepita (C₀) representa a variância não explicada ou ao acaso, freqüentemente causada por erros de medições ou variações dos atributos que não podem ser detectadas na escala de amostragem (Trangmar et al., 1985; Vieira, 2000). A relação em percentagem entre o efeito pepita e o patamar do semivariograma (C) pode indicar o grau de dependência espacial (Trangmar et al., 1985). Neste estudo, observou-se que, em geral, os atributos estudados apresentaram moderado grau de dependência espacial nas duas profundidades, conforme critérios estabelecidos por Cambardella et al. (1994). Entretanto, foram observados maiores valores da relação C₀/(C₀ + C) na profundidade de 0,00—0,20 m, indicando que o solo pode apresentar variabilidade em escalas menores do que a utilizada neste estudo e que amostragens mais intensivas podem revelar maior continuidade espacial dos atributos analisados.

Os teores de P e Ca²⁺ e os valores de CTC apresentaram forte grau de dependência espacial na camada mais profunda (Quadro 2). Segundo Cambardella et al. (1994), os atributos que apresentaram forte dependência espacial são mais influenciados por propriedades intrínsecas do solo, como textura e mineralogia, enquanto os que apresentaram fraca dependência são mais influenciados por fatores externos, tais como: aplicações de fertilizante, preparo e cultivo do solo, ou seja, pelo manejo do solo.

Os parâmetros dos modelos de semivariogramas ajustados foram utilizados para estimar valores em locais não amostrados por meio da krigagem. A magnitude do efeito pepita é importante na krigagem, pois, quanto maior for a diferença do efeito pepita em relação ao patamar do semivariograma, maior a continuidade do fenômeno, menor a variância da estimativa, ou maior a confiança que se pode ter na estimativa (Isaaks & Srivastava, 1989).

Com os valores estimados por krigagem, foram gerados mapas de distribuição espacial dos atributos analisados (Figuras 1, 2, 3 e 4), e os valores observados foram divididos em quatro classes para todos os atributos. Observou-se um arranjo de distribuição espacial dos atributos mais estruturados na profundidade de 0,00—0,20 m (Figuras 1 e 2). Em geral, foram observados maiores valores dos atributos analisados na parte esquerda dos mapas, correspondente às distâncias aproximadas entre 0 e 1.500 m na direção do eixo das abscissas, o que possibilita, de maneira menos rigorosa, dividir a área em duas partes, com diferentes níveis de fertilidade. Entretanto, de maneira mais rigorosa, dentro de cada uma dessas subáreas, podem-se definir várias zonas homogêneas de manejo específico, com o objetivo de, por exemplo, efetuar a aplicação de calcário e fertilizantes a taxas variáveis, proporcionando, desta forma, um manejo mais eficiente e econômico da cultura.

Na profundidade de 0,60—0,80 m (Figuras 3 e 4), observou-se uma distribuição espacial mais aleatória dos atributos analisados, indicando que os Latossolos, conceituadamente considerados homogêneos, apresentaram variabilidade espacial de seus atributos mesmo em horizontes pouco influenciados pelo manejo do solo. Contudo, ainda foi possível perceber, de maneira geral, maiores valores dos atributos na parte esquerda dos mapas, semelhantemente ao que se observou na camada superficial, confirmando a existência de diferentes zonas de manejo na área.

CONCLUSÕES

1. Os alcances de dependência espacial para os atributos químicos do solo na camada de 0,60—0,80 m de profundidade foram bem menores, quando comparados àqueles referentes à camada de 0,00—0,20 m. Estes resultados demonstraram maior descontinuidade na distribuição espacial dos atributos do solo na camada de 0,60—0,80 m de profundidade dos Latossolos, indicando que essa classe de solos não apresentou homogeneidade de seus atributos como conceituadamente a ela foi atribuída.

2. O manejo no solo realizado na área ao longo do tempo alterou a dependência espacial dos atributos do solo na camada superficial de forma a diminuir a variabilidade espacial dos atributos químicos do solo em relação à camada mais profunda.

3. A investigação da variabilidade espacial de atributos químicos e teor de argila da camada superficial e subsuperficial dos solos da área estudada proporcionou a visualização e definição de zonas homogêneas de manejo, o que permite a adoção do sistema de agricultura de precisão.

AGRADECIMENTOS

Os autores agradecem à Fundação de Amparo à Pesquisa do estado de São Paulo (FAPESP), pela bolsa de estudos à segunda autora; à Fundação para o Desenvolvimento da UNESP (Fundunesp), pelo auxílio financeiro à pesquisa; à Fazenda Santa Isabel, na pessoa do Engenheiro-Agrônomo Paulo de Araújo Rodrigues, pela concessão da área experimental e total apoio nos trabalhos de campo.

LITERATURA CITADA

Recebido para publicação em maio 2003 e aprovado em outubro 2004.

Datas de Publicação

Histórico