Dependência espacial da resistência do solo à penetração e do teor de água do solo sob cultivo contí (página 2)

A área de estudo localiza-se no nordeste do Estado de São Paulo, no município de Jaboticabal, SP. As coordenadas geográficas são 21°22' de latitude Sul e 48°15' de longitude Oeste, com altitude média de 600m. O clima da região, segundo a classificação de Köppen, é do tipo mesotérmico com inverno seco (Cwa). A precipitação pluviométrica anual média é de 1.400 mm com chuvas concentradas no período de novembro a fevereiro. O relevo é predominantemente suave ondulado com declividades médias variando de 6 a 11%. O solo é classificado como Latossolo Vermelho distrófico, textura argilosa (LVd), desenvolvido sobre um arenito do Grupo Bauru, Formação Adamantina.

A área estudada tem 40,5ha e vem sendo cultivada com cana-de-açúcar há mais de 40 anos e atualmente está submetida à colheita mecanizada (cana crua). Uma malha com intervalos regulares de 10m, perfazendo um total de 100 pontos georeferenciados em 1ha, foi disposta no topo da área e as amostras foram coletadas nas profundidades de 0,00-0,15m, 0,15-0,30m e 0,30-0,45m, em outubro de 2004.

O teor de água no solo foi obtido pelo método gravimétrico em amostras deformadas (EMBRAPA, 1997). Para determinação da resistência do solo à penetração, utilizou-se um penetrômetro de impacto modelo IAA/Planalsucar com ângulo de cone de 30°. A transformação da penetração da haste do aparelho no solo (cm/impacto) em resistência à penetração foi obtida pela fórmula dos "holandeses", segundo STOLF (1991), a seguir:

onde R é a resistência à penetração, kgf cm^-2 (kgf cm^-2 * 0,098 = MPa); M a massa do êmbolo, 4kg (Mg – 4kgf); m a massa do aparelho sem êmbolo, 3,2kg (mg – 3,2kgf); h a altura de queda do êmbolo, 40cm; x a penetração da haste do aparelho, cm/impacto, e A a área do cone, 1,29cm².

Realizaram- se determinações da resistência do solo à penetração em duas épocas. A primeira foi 24 horas após uma chuva de 38mm (Época 1) e a segunda, 72 horas após a mesma chuva (Época 2). Este procedimento permitiu simular diferentes teores de água no solo e relacioná-los com a variabilidade espacial da resistência do solo à penetração.

Os atributos resistência do solo à penetração e teor de água no solo foram analisados através da análise estatística descritiva e de técnicas de geoestatística. A hipótese de normalidade dos dados foi testada pelo teste de Kolmogorov-Smirnov, por meio do programa computacional SAS (SCHLOTZHAVER & LITTELL, 1997). A dependência espacial foi analisada por meio de ajustes de semivariogramas através do programa GS⁺, com base na pressuposição de estacionariedade da hipótese intrínseca (VIEIRA et al., 1983), a qual é estimada por:

onde 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 estimado é 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₀; patamar, C₀+C₁; e o alcance, a). Os modelos de semivariogramas foram o esférico, o exponencial, o linear e o gaussiano (McBRATNEY & WEBSTER, 1986), sendo ajustados pelo programa de computador GS⁺ (ROBERTSON, 1998). Para analisar o grau da dependência espacial dos atributos em estudo, utilizou-se a classificação de CAMBARDELLA et al. (1994), em que são considerados de dependência espacial forte os semivariogramas que têm um efeito pepita menor ou igual a 25% do patamar, moderada quando está entre 25 e 75% e fraca quando for maior que 75%.

No caso de mais de um modelo para o mesmo semivariograma, foi escolhido o de maior valor do coeficiente de correlação obtido pelo método de validação cruzada e o de melhor coeficiente de determinação (R²).

Os valores médios de resistência do solo à penetração na época 1 e na profundidade de 0,00-0,15m, foram classificados como moderados de acordo com (SOIL SURVEY STAFF, 1993) e nas profundidades de 0,15-0,30m e 0,30-0,45m como alto (Tabela 1). Na época 2, os valores foram classificados como alto em todas as profundidades estudadas. Houve aumento da resistência do solo à penetração da profundidade de 0,00-0,15m para a profundidade de 0,15-0,30m, passando de 1,99 para 2,93MPa na época 1 e 2,30 para 3,77MPa na época 2, com diminuição da resistência do solo à penetração na profundidade de 0,30-0,45m, o mesmo comportamento foi verificado nas duas épocas de amostragem estudadas.

Observa-se, na tabela 1, valores altos de resistência do solo à penetração nas duas épocas de amostragem. Embora os teores de água no solo estejam próximo da capacidade de campo (época 1),os valores de resistência à penetração são considerados impeditivos ao desenvolvimento radicular das culturas, resultados semelhantes foram encontrados em latossolo por TORMENA (1998). DEXTER (1987) afirma que a compactação do solo é mais prejudicial em solo seco, e em condições de maior teor de água no solo pode haver crescimento radicular em valores de resistência do solo à penetração superiores a 4,0MPa.

O teste de normalidade indicou normalidade para teor de água no solo em todas as profundidades e nas duas épocas amostradas e, resistência do solo à penetração nas profundidades de 0,15-0,30m e 0,30-0,45m na época 1 (Tabela 1). Os valores da média e mediana, para todas as variáveis, estão próximos, mostrando distribuições simétricas, o que pode ser confirmado pelo coeficiente de assimetria próximo de zero.

A variável resistência do solo à penetração apresentou coeficiente de variação (CV) alto de acordo com a classificação proposta por WARRICK & NIELSEN (1980) em todas profundidades e épocas estudadas (Tabela 1), concordando com os resultados obtidos por SOUZA et al. (2001), UTSET & CID (2001) e ABREU et al. (2003). Para o teor de água no solo o CV foi baixo em ambas as profundidades para a época 2 e médio em todas as profundidades na época 1. O CV para resistência do solo à penetração na época 1 foi maior quando comparado com a época 2, isto deve-se provavelmente ao maior teor de água no solo na época 1, corroborando com os resultados de UTSET & CID (2001) que verificaram para condições de solo seco um CV de 25%, enquanto sob condição de solo úmido, o CV foi de 80%.

A resistência do solo à penetração na época 2 (Tabela 2 e Figura 1), apresentou dependência espacial em todas as profundidades estudadas expressa por meio dos modelos de semivariogramas. Na época 1 a resistência do solo à penetração apresentou distribuição aleatória (efeito pepita puro), provavelmente devido ao alto teor de água no solo. Na época 2, os dados da variável resistência do solo à penetração ajustaram-se ao modelo exponencial na profundidade 0,00-0,15m e esférico nas profundidades de 0,15-0,30m e 0,30-0,45m. Para o teor de água no solo avaliados nas duas épocas de amostragem (Tabela 2 e Figura 2), os dados ajustaram-se ao modelo esférico nas profundidades de 0,00-0,15m e 0,15-0,30m e exponencial na profundidade de 0,30-0,45m.

Na condição de menor teor de água no solo, os dados apresentaram dependência espacial (Tabela 2 e Figura 2). Resultados semelhantes com aqueles obtidos por UTSET & CID (2001) os quais avaliando a resistência do solo à penetração, observaram o modelo esférico ajustando-se aos dados sob condição de solo seco e um efeito pepita puro em condições de solo úmido. Os modelos esférico e exponencial que se ajustaram aos dados deste estudo concordam com os resultados de várias pesquisas que indicam estes modelos como os de maior ocorrência para atributos do solo (SOUZA et al., 1997; SALVIANO et al., 1998; SOUZA et al., 2004).

A variabilidade espacial para resistência do solo à penetração, na época 1, é muito grande, ou seja, o esquema amostral de 10m entre pontos não foi suficiente para captar a dependência espacial para esta variável (Figura 1). Portanto, recomenda-se uma malha mais adensada que a utilizada neste estudo para avaliar a variabilidade espacial da resistência do solo à penetração, quando o teor de água no solo estiver próximo ou acima da capacidade de campo.

O efeito pepita é um parâmetro importante do semivariograma e indica uma variabilidade não explicada considerando a distância de amostragem utilizada (McBRATNEY & WEBSTER, 1986). Esse parâmetro pode ser expresso como percentagem do patamar, e tem por objetivo facilitar a comparação do grau de dependência espacial das variáveis em estudo (TRANGMAR et al., 1985). A análise da relação C₀/(C₀+C₁) mostrou que as variáveis resistência do solo à penetração e teor de água no solo apresentaram dependência espacial forte na profundidade de 0,15-0,30m e moderada nas profundidades de 0,00-0,15m e 0,30-0,45m nas épocas amostradas (Figuras 1 e 2).

O alcance (a) indica a distância limite entre pontos correlacionados entre si. Pontos coletados com distância maiores que o alcance são independentes e, para sua análise, pode-se utilizar a estatística clássica (VIEIRA, 2000). As variáveis em estudo apresentaram alcances da ordem de 48,6m para resistência do solo à penetração na profundidade de 0,00-0,15m e época 2 e 20,4m para teor de água no solo na profundidade de 0,15-0,30m e época 1. Observaram-se duas classes de valores para o alcance, na profundidade de 0,00-0,15m o alcance variou de 30 a 50m e nas profundidades de 0,15-0,30m e 0,30-0,45m variou de 19 a 30m em ambas as épocas de estudo, respectivamente (Figuras 1 e 2). MERCANTE et al. (2003) estudaram a variabilidade espacial e temporal da resistência do solo à penetração em áreas com e sem manejo químico localizado em latossolo, e constataram valores de alcances variando de 60 a 80m para uma área de 1ha, sendo que os maiores alcances foram observados nas áreas com manejo químico localizado e resultados similares a este trabalho nas áreas sem manejo químico localizado.

O teor de água no solo afeta a variabilidade espacial da resistência do solo à penetração e no maior teor de água no solo, houve uma menor dependência espacial 24 horas após uma chuva de 38mm.

A malha de amostragem deve ser mais adensada em relação à utilizada neste estudo quando o teor de água no solo for superior ao observado 72 horas após a chuva, a fim de se avaliar a variabilidade espacial para a resistência do solo à penetração.

À Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP), pela concessão de bolsa de pesquisa a Souza.