Variabilidade espacial da textura de um latossolo vermelho eutroférrico sob cultivo de cana-de-açúca (página 2)

A área de estudo localiza-se no nordeste do Estado de São Paulo, no município de Guariba - SP. As coordenadas geográficas são 21º19' de latitude sul e 48º13' de longitude oeste, com altitude média de 640 m. O clima da região, segundo a classificação de Köeppen, é do tipo mesotérmico com inverno seco (Cwa), com precipitação média de 1.400 mm e chuvas concentradas no período de novembro a fevereiro. A vegetação natural era constituída por floresta tropical subcaducifólia e mata ciliar.

O relevo é predominantemente suave ondulado com declividades médias variando de 3 a 8%. A área experimental vem sendo cultivada há mais de 30 anos com cana-de-açúcar. O solo da área foi classificado como Latossolo Vermelho eutroférrico, textura muito argilosa (LVef) (EMBRAPA, 1999).

No terço inferior da encosta, por meio da análise detalhada das formas do relevo e da variação do gradiente do declive, caracterizaram-se dois compartimentos (I e II, Figura 1), sob os quais os solos foram amostrados nos pontos de cruzamento de uma malha, com intervalos regulares de 50 m, perfazendo um total de 206 pontos nas profundidades de 0-0,2 m e 0,6-0,8 m. Acrescentaram-se mais quatro submalhas com 25 pontos (refinamento), com intervalos regulares de 10 m (100 pontos), totalizando-se 306 pontos de amostragem em uma área de 42 ha. A composição granulométrica foi determinada pelo método da pipeta, utilizando solução de NaOH 0,1 N como dispersante químico e agitação mecânica em aparato de baixa rotação por 16 h, seguindo metodologia proposta pela EMBRAPA (1997), com modificações. Utilizaram-se 15 g de areia lavada com diâmetro entre 0,5 -1,0 mm para auxiliar na dispersão física das partículas. A fração argila foi separada por sedimentação, de acordo com a lei de Stokes, sendo a fração silte determinada por diferença. A fração areia, após secagem e pesagem, foi subdividida em cinco subfrações: areia muito grossa (AMG), areia grossa (AG), areia média (AM), areia fina (AF) e areia muito fina (AMF), em peneiras de 2-1 mm; 1-0,5 mm; 0,5-0,25 mm; 0,25-0,10 mm e 0,10-0,05 mm, respectivamente.

Os atributos granulométricos do solo foram analisados por meio da análise estatística descritiva e de técnicas geoestatísticas. A hipótese de normalidade dos dados foi testada pelo teste de Kolmogorov-Smirnov, SAS (1995). A dependência espacial foi analisada por meio de ajuste de semivariograma (VIEIRA et al., 1983 e ROBERTSON, 1998), com base nas pressuposições de estacionariedade da hipótese intrínseca, estimado por:

em que, N (h) é o número de pares experimentais de observações Z(x_i) e Z (x_i + h) separados por uma distância h. O semivariograma é representado pelo gráfico versus h. Do ajuste de um modelo matemático aos valores calculados de são estimados os coeficientes do modelo teórico para o semivariograma (o efeito pepita, C₀; variância estrutural, C₁; patamar, C₀ + C₁, e o alcance, a). Os modelos de semivariogramas considerados foram o esférico, o exponencial e o gaussiano (McBRATNEY & WEBSTER, 1986). A razão entre a variância estrutural e o patamar C₁/(C₁+C₀), expressa em porcentagem, permitiu a classificação do grau de dependência espacial, a qual, segundo ROBERTSON (1998), é considerada fraca se essa razão for < 25% , moderada quando essa razão estiver entre 25 e 75%, e forte se a razão for > 75%.

Esses modelos foram ajustados por meio do programa GS⁺ (ROBERTSON, 1998). Em caso de dúvida entre mais de um modelo para o mesmo semivariograma, considerou-se o melhor R² (coeficiente de determinação) e a menor SQR (soma de quadrados do resíduo). Para a elaboração dos mapas de distribuição espacial das variáveis, foi utilizado o programa Surfer (GOLDEN SOFTWARE, 1999).

Os valores da mediana e média das frações granulométricas analisadas foram similares (Tabela 1), indicando distribuições simétricas, o que pode ser confirmado pelos valores de assimetria próximos de zero, com exceção da areia grossa (AG), em ambas as profundidades, e areia fina, na profundidade de 0,6-0,8 m.

A variável argila apresentou distribuição normal nas duas profundidades estudadas e silte na profundidade de 0,6-0,8 m; as demais variáveis estudadas não apresentaram normalidade. VIEIRA (1997), ARAUJO (2002) e GONÇALVES & FOLEGATTI (2002) encontraram normalidade para a variável argila, trabalhando com amostragem em malha regular. Resultados semelhantes para a variável silte também foram observados por VIEIRA (1997), GONÇALVES et al. (2001) e SOUZA et al. (2003). A não-normalidade para as areias também foi verificada por POCAY (2000), ARAUJO (2002), EGUCHI et al. (2002) e SOUZA et al. (2003).

Como salientado por WARRICK & NIELSEN (1980) e WEBSTER (1985), em se tratando de dados obtidos na natureza, o ajuste a uma distribuição teórica é apenas aproximado. De acordo com CRESSIE (1991), a normalidade dos dados não é uma exigência da geoestatística, sendo conveniente apenas que a distribuição não apresente caudas muito alongadas, o que pode comprometer os resultados. Assim, a partir dos resultados da análise exploratória dos dados, é possível admitir que essas variáveis apresentam distribuições suficientemente simétricas para o estudo geoestatístico.

Os valores do coeficiente de variação (C.V.) referentes à granulometria indicaram variabilidade baixa para argila (4,77 - 3,05%); resultados semelhantes também foram encontrados por COELHO FILHO et al. (2001) e GONÇALVES et al. (2001); média para silte (16,64 - 12,71%), AT (14,98 - 12,25%), AM (22,09 - 16,22%), AF (14,92 - 10,52%) e AMF (21,23 - 19,78%) e alta para a AG (40,58 - 35,12%), concordando com os resultados obtidos por POCAY (2000) e SOUZA (2001), nas profundidades de 0-0,2 m e 0,6-0,8 m, respectivamente. O coeficiente de variação foi menor na camada de 0-0,2 m em relação à camada de 0,6-0,8 m para todas as variáveis estudadas. Os resultados da Tabela 1 revelaram que o solo é de textura muito argilosa quando avaliado pela média e valores de mínimo e máximo.

Argila e areia dos latossolos são atributos estáveis, ou seja, modificam-se pouco ao longo do tempo por ação do intemperismo, portanto apresentam baixos valores de C.V. No presente estudo, o atributo areia apresentou maiores valores de C.V., provavelmente indicando a influência da posição da paisagem, devido a processos erosionais ou deposicionais, associados ao fluxo de água, concordando com os estudos de POCAY (2000) e SOUZA (2001).

Para as variáveis argila, silte, areia grossa, areia fina e areia muito fina, na profundidade de 0-0,2 m e com exceção da argila na profundidade 0,6-0,8 m, estimaram-se os semivariogramas para os resíduos da tendência linear. A remoção dessa tendência foi necessária porque o semivariograma para os dados originais apresentou crescimento sem limites com a distância.

Os resultados da análise geoestatística (Tabela 2 e Figura 2) mostraram que todas as variáveis analisadas apresentaram dependência espacial nas duas profundidades estudadas. Na profundidade 0-0,2 m, os dados das variáveis argila, silte, AG e AM ajustaram-se ao modelo esférico, e as variáveis AT, AF e AMF ajustaram-se ao modelo exponencial, com R² acima de 0,94. Na profundidade 0,6-0,8 m, todas as variáveis foram ajustadas ao modelo exponencial, com exceção da argila e AF que se ajustaram ao modelo esférico, com R² acima de 0,90. O modelo que se ajustou com maior freqüência aos dados granulométricos foi o exponencial, juntamente com o esférico, sendo que BERTOLANI et al. (2000) afirmam que são esses os modelos que mais se ajustam aos dados de atributos físicos do solo.

A análise da relação C₁/(C₀ + C₁) dos atributos granulométricos do solo (Tabela 2) mostrou que as variáveis argila, silte e AMF apresentaram dependência espacial moderada, e as variáveis AT, AG, AM e AF apresentaram dependência espacial forte na profundidade de 0-0,2 m; na profundidade de 0,6-0,8 m, todas as variáveis em estudo tiveram dependência espacial moderada, segundo a classificação proposta por ROBERTSON (1998). As variáveis apresentaram diferentes alcances, sendo que, na camada superficial, as variáveis AG e AM apresentaram os maiores alcances, e silte e AF os menores (366,8; 287,7; 157,0 e 199,6 m, respectivamente). Na profundidade de 0,6-0,8 m, os alcances foram menores em relação à profundidade de 0-0,2 m, com máximo de 215,2 m para a variável argila e mínimo de 57,9 m para variável AT, indicando a maior descontinuidade na distribuição espacial das propriedades do solo na camada mais profunda. Isso está de acordo com um importante conceito em gênese do solo, que o maior tempo relativo de exposição dos horizontes superficiais condiciona sua maior taxa de intemperismo (BUOL et al., 1997). Esse fato pode ter contribuído para aumentar o alcance da dependência espacial, caracterizando maior continuidade na distribuição das variáveis na camada mais superficial da área estudada, concordando com estudo de SOUZA (2001) e ARAUJO (2002), trabalhando com malhas regulares espaçadas em 50 m.

Os parâmetros dos modelos de semivariogramas ajustados (Figura 2) foram utilizados para estimar valores em locais não amostrados por meio da krigagem (Figura 3). Os mapas espaciais das variáveis argila e areia total têm comportamento inverso, quando comparados entre si. O maior teor da variável argila está no compartimento I, relativamente menos declivoso, em relação ao compartimento II. A variável areia total apresenta o mesmo comportamento das variáveis AG, AM e AF com maiores teores no compartimento de maior declive da área, influenciado pelas maiores taxas de erosão nessa posição da paisagem (SOUZA et al., 2002).

Essa informação sugere retenção de água e teores de íons diferenciados para esses compartimentos, podendo interferir no comportamento da cultura. FRAISSE et al. (1999) e FRIDGES et al. (2000) afirmam que a utilização da altitude, declividade e formas da paisagem do solo para delinear zonas de manejo para atributos granulométricos e químicos, em diferentes tipos de solos, são de extrema importância. MACHADO (1994) afirma que a posição topográfica dos solos na paisagem, influindo na drenagem do solo, é um dos principais fatores condicionantes da variabilidade espacial dos solos.

A amplitude das variáveis estudadas enfoca que, numa mesma classe de solos, conceitualmente considerada altamente intemperizada (Latossolo Vermelho eutroférrico), se observa que ocorre variabilidade espacial dos atributos granulométricos. Considerando que os teores de argila estão relacionados a compartimentos com declive diferenciado e que esse atributo é estável no solo, sugere-se que essas pequenas variações do gradiente do declive possam contribuir para a identificação de limites de precisão entre locais específicos da paisagem.

A variabilidade do solo medida pelo coeficiente de variação registrou-se baixa para argila, média para silte, areia total, areia média, areia fina e areia muito fina, e alta para areia grossa, nas profundidades estudadas.

As variáveis estudadas apresentaram estrutura de dependência espacial. Os maiores alcances foram observados na profundidade de 0-0,2 m para todos os atributos estudados, refletindo os efeitos do maior tempo de exposição dos horizontes superficiais nas maiores taxas relativas de intemperismo.

Os latossolos sob cultivo de cana-de-açúcar, apesar de serem considerados homogêneos, apresentam variabilidade diferenciada para os atributos granulométricos.

A identificação de compartimentos da paisagem na área de estudo mostrou-se muito eficiente na compreensão da variabilidade espacial dos atributos estudados.