Forma da paisagem como critério para otimização amostral de latossolos sob cultivo de cana-de-açúcar (página 2)

O estudo das relações solo-geomorfologia tem sido realizado por meio de vários modelos de paisagem (Troeh, 1965; Ruhe et al., 1967; Daniels et al., 1971). A associação das formas da paisagem (pedoforma) à variabilidade espacial dos atributos dos solos tem contribuído na identificação e mapeamento de áreas mais homogêneas, com limites mais precisos entre elas (Souza et al., 2003), o que permite que técnicas agronômicas possam ser transferidas com facilidade e economia para ambientes semelhantes.

O comportamento diferenciado dos atributos dos solos em diferentes segmentos da paisagem pode ser entendido pela caracterização da variabilidade espacial dos solos por meio de técnicas de geoestatística, identificando locais específicos na paisagem. Porém, o custo da coleta e análise das amostras a fim de caracterizar a variabilidade espacial é um fator limitante em estudos de reconhecimento e mapeamento de locais específicos de manejo. Assim, técnicas com o objetivo de otimizar o esquema de amostragem e baratear custos são importantes (Webster & Oliver, 1990).

As pedoformas, em áreas de mesma classe de solo e mesmo histórico de manejo durante longos anos, influenciam a variabilidade dos atributos dos solos de maneira sistemática. Esta premissa pode ser aplicada às técnicas de simulação de resultados (Van Groenigen et al., 1999), possibilitando otimização do número de amostras e transferência de informação para ambientes similares.

Uma das grandes dificuldades encontradas nos estudos de variabilidade espacial do solo tem sido a determinação do espaçamento ideal de amostragem, o que, em algumas situações, inviabiliza a adoção das técnicas da agricultura de precisão. Assim, o estudo dos aspectos da amostragem do solo, a fim de subsidiar definições sobre a utilização e recomendação destas técnicas em ambientes distintos, tem sido preocupação constante dos pesquisadores (Webster & Oliver, 1990; Van Groenigen et al., 1999; Lark, 2000).

Em áreas côncavas, independentemente do histórico de manejo, há mais variabilidade de atributos químicos que em áreas lineares (Souza et al., 2003). Isto implica em maior número de amostras de solo a serem coletadas nestas áreas. Então, na avaliação da variabilidade espacial, o espaçamento com intervalos menores é necessário em áreas côncavas e uma das maneiras de avaliar esta premissa é pela otimização de amostragem por meio do programa Sanos (Van Groenigen et al., 1999).

O objetivo deste trabalho foi utilizar a forma da paisagem como critério auxiliar na otimização do esquema amostral na avaliação dos atributos químicos de latossolos, em áreas sob cultivo de cana-de-açúcar.

A área do presente estudo localiza-se na fazenda Santa Isabel, Município de Jaboticabal, SP, inserida na província geomórfica das Cuestas Basálticas, limite com o Planalto Ocidental Paulista. A região apresenta relevo suave ondulado com altitude média de 600 m. As coordenadas geográficas da área de estudo são: 21º17' a 21º18' S e 48º8' a 48º10' W. O clima foi classificado, segundo Köppen, como do tipo Aw, definido como tropical de verão chuvoso e inverno seco, com temperatura variando de 19,1ºC no mês mais frio a 24,8ºC no mês mais quente.

O material geológico na área está relacionado aos basaltos do grupo São Bento, formação Serra Geral, e material retrabalhado do basalto. O solo predominante na área foi classificado como Latossolo Vermelho eutroférrico (LVef) (Embrapa, 1999).

A região foi percorrida para identificação dos modelos da paisagem. Desse modo, escolheu-se uma área contínua, sob cultivo de cana-de-açúcar há mais de 30 anos, caracterizada por possuir duas pedoformas, segundo classificação de Troeh (1965): côncava, que ocorre nas posições mais elevadas; e linear, constituída pelos segmentos escarpa, meia encosta e encosta inferior, que grada para o aluvião do Rio Mogi-Guaçu (Figura 1).

As pedoformas estudadas são apresentadas na Figura 1, através do perfil altimétrico. Assim, na pedoforma côncava, o perfil de curvatura é convergente para o centro e caracteriza-se como um sistema aberto cuja altitude varia entre 620 e 640 m. A pedoforma linear possui perfil retilíneo, caracterizando uma encosta com declividade suave nos segmentos da meia encosta e mais declivosa no segmento da escarpa e encosta inferior. Esta pedoforma ocorre em seqüência à pedoforma côncava, seguindo a encosta até a área coluvial.

Os atributos químicos em ambas pedoformas foram estimados a partir da coleta de amostras de solos em uma malha regular nas profundidades de 0,0–0,2 m e 0,6–0,8 m, num total de 421 pontos em uma área de 300x3.000 m, com espaços de 50x50 m, abrangendo 94 hectares. As amostras foram secadas ao ar, destorroadas e passadas em peneira com malha de 2 mm de abertura. Nelas foram determinados pH (CaCl₂), Ca, Mg, K e o P disponível, pelo método da resina trocadora de íons, e o teor de matéria orgânica conforme Raij et al. (2001). Com base nos resultados das análises químicas, calcularam-se a CTC e saturação por bases (V%). Foram realizadas as análises estatísticas descritivas dos dados obtidos – média, mediana, coeficiente de variação, assimetria e curtose e teste de normalidade –, utilizando-se o programa estatístico SAS (SAS Institute, 1995).

Cada atributo do solo foi submetido ao ajuste de modelos matemáticos na definição dos semivariogramas pelo programa GS⁺ (Robertson, 1998) e técnica de interpolação de dados (krigagem), pelo programa Surfer (Golden Software, 1999), para comparar o grau de variabilidade dos atributos químicos estudados. A análise do grau de dependência espacial dos atributos foi realizada segundo Cambardella et al. (1994). Nesta análise, são considerados de dependência espacial forte os semivariogramas que têm um efeito pepita <25% do patamar, moderada, entre 25% e 75%, e fraca, >75%.

A otimização do esquema amostral na estimação dos atributos químicos do solo foi obtida por meio do programa Sanos for Windows 0.1, o qual baseia-se em um algoritmo de simulação espacial (Van Groenigen et al., 1999) que utiliza vários critérios quantitativos de otimização, sendo o utilizado neste trabalho o da minimização da variância da krigagem ordinária e o semivariograma estimado.

Segundo Troeh (1965), a pedoforma côncava condiciona fluxo de água convergente, enquanto a pedoforma linear caracteriza-se por processos de escoamento superficial. Ressalta-se que na área em estudo, a pedoforma côncava possui acúmulo relativo de água, já que possui um sistema aberto no seu interior para uma das laterais (Figura 1).

A CTC na profundidade de 0,0–0,2 m e V% nas duas profundidades estudadas apresentaram distribuição normal (Tabela 1). Oliveira et al. (1999), estudando a variabilidade espacial de atributos químicos, observaram distribuição normal somente em relação à CTC efetiva dos onze atributos químicos do solo estudado. A profundidade de 0,0–0,2 m, por ser afetada pelo processo de mecanização, manejo das culturas e erosão, pode, provavelmente, afetar o comportamento dos atributos químicos do solo, em ambas pedoformas. Logo, o arrastamento das partículas do solo pela erosão pode ter contribuído no comportamento desuniforme nos teores dos atributos químicos analisados. No entanto em se tratando de dados obtidos na natureza, o ajuste de uma distribuição teórica é apenas aproximado (Cressie, 1991).

Os valores do critério de classificação para o coeficiente de variação (CV) proposto por Warrick & Nielsen (1980) para avaliar a dispersão dos dados revelaram-se baixos (CV<12%) para o pH nas duas profundidades estudadas (Tabela 1). Resultados semelhantes foram encontrados por Carvalho et al. (1998) e Salviano et al. (1998). As variáveis V% na profundidade de 0,6–0,8 m, e CTC e MO em ambas as profundidades apresentaram CV médio (12%<CV<24%), corroborando os resultados de Souza et al. (1997). Coeficientes de variação altos (CV>24 %) foram encontrados em relação ao P, K, Ca e Mg nas duas profundidades e V% na profundidade de 0,0–0,2 m. Souza et al. (1997), Salviano et al. (1998) e Carvalho et al. (2003) também encontraram CV altos quanto ao P, K, Ca e Mg, possivelmente por causa do manejo do solo e da aplicação de fertilizantes.

Os resultados da análise geoestatística mostraram que todos os atributos analisados apresentaram dependência espacial nas profundidades estudadas (Tabela 2, Figuras 2 e 3). O modelo esférico ajustou-se à semivariância estimada de todos os atributos químicos com exceção do K e da MO que se ajustaram ao modelo exponencial na profundidade de 0,0–0,2 m. Na profundidade de 0,6–0,8 m, o modelo exponencial ajustou-se à semivariância de todos os atributos químicos, com exceção do Mg que se ajustou ao modelo esférico. Todas as variáveis apresentaram coeficiente de determinação (R²) acima de 0,83, ou seja, 83% da variabilidade existente nos valores da semivariância estimada são explicadas pelos modelos ajustados. Os modelos esférico e exponencial ajustaram-se à semivariância dos atributos em estudo, estando de acordo com resultados que indicam que esses modelos se ajustam melhor aos atributos do solo (Trangmar et al., 1985; Salviano et al., 1998; Carvalho et al., 2003).

A análise da relação C₀/(C₀+C₁) dos atributos químicos do solo mostrou que todos os atributos apresentaram dependência espacial moderada, segundo a classificação proposta por Cambardella et al. (1994), com exceção da variável Ca, que apresentou forte dependência espacial (Tabela 2). O alcance representa a distância em que os pontos amostrais estão dependentes espacialmente entre si (Journel & Huijbregts, 1991), ou seja, os pontos localizados numa área de raio menor ou igual ao alcance são mais similares. Houve dois grupos de alcances nas profundidades de 0,0–0,2 m e 0,6–0,8 m, ou seja, as variáveis pH, CTC, V%, Mg, Ca e P apresentaram valores próximos de alcance (744–558, 747–660, 785–893, 808–556, 837–774 e 880–699 m, respectivamente), e as variáveis K e MO também apresentaram valores mais próximos entre si (264–183 e 426–141 m, respectivamente).

Os alcances de dependência espacial na profundidade de 0,6–0,8 m foram menores em relação à profundidade de 0,0-0,2 m com exceção da variável V%, indicando a maior descontinuidade na distribuição espacial dos atributos químicos do solo na camada mais profunda (Tabela 2). O manejo do solo, pelas aplicações de fertilizantes e revolvimento do mesmo, caracterizou uma maior continuidade na distribuição dos atributos químicos na camada mais superficial da área estudada.

A maior concentração de pontos indicados pelo programa Sanos 0.1, revela os locais de maior variabilidade espacial dos atributos avaliados (Tabela 3). Observa-se que na caracterização da maioria dos atributos químicos estudados, há necessidade de um maior número de pontos de amostragem na pedoforma côncava. Isso reforça a idéia de que, na pedoforma côncava, as condições pedogeomórficas favorecem maior variabilidade espacial relativa dos atributos químicos (Zebarth et al., 2002; Souza et al., 2003).

Ao se admitir a hipótese de possível adubação diferenciada nas pedoformas, o maior número de pontos amostrais recomendados para a maioria dos atributos químicos estudados pelo programa Sanos 0.1, na profundidade de 0,6-0,8 m, corrobora a informação de maior variabilidade nos solos de locais côncavos na paisagem (Tabela 3). Aspecto relevante desse programa é a plotagem de pontos em locais de real necessidade, já que os cálculos são realizados com base nas coordenadas geográficas. Souza et al. (2004), estudando a influência da pedoforma na variabilidade espacial de atributos físicos e hídricos, observaram que pequenas variações na pedoforma condicionam variabilidade diferenciada em relação aos atributos físicos e hídricos do solo, sendo que a variabilidade espacial foi maior na pedoforma côncava em relação à linear, influenciada provavelmente pelo fluxo de água.

Souza et al. (2003) também constataram maior variabilidade da erosão na pedoforma côncava. Segundo Resende et al. (1997), a pedoforma côncava possui uma instabilidade maior nas partes mais altas (francamente de remoção) e maior estabilidade nas áreas mais baixas (de acúmulo). Essas diferentes posições criam situações locais variadas. Nesse sentido, estas informações são válidas como critério auxiliar no planejamento de amostragem de solo, que, no caso da agricultura de precisão, representa o fator limitante nos estudos de reconhecimento e mapeamento de locais específicos de manejo. Kravchenko & Bullock (2000), estudando a correlação da produção de milho e feijão com variações da topografia, verificaram que 40% da variabilidade da produção destas culturas e dos atributos do solo foi explicada pelas curvaturas do terreno. Sobieraj et al. (2002) e Zebarth et al. (2002) encontraram relação de atributos químicos com diferentes posições na paisagem.

O número de amostras varia segundo o relevo. Os locais da paisagem sob a mesma classe de solo e manejo semelhante, influenciados pela pedoforma e declive, apresentam padrões de variabilidade espacial dos atributos químicos do solo diferenciados (Tabela 3, Figuras 2 e 3). Deste modo, deve-se considerar a paisagem como fator relevante em estratégias de amostragem do solo, como também a incorporação de um modelo preditivo de variabilidade espacial. Segundo Souza et al. (2003), as diferenças na distribuição espacial dos atributos do solo nas diferentes pedoformas estão associadas às variações no relevo (fluxo de água), sendo este o principal responsável pelas maiores variabilidades na pedoforma côncava.

Malo et al. (1974) relataram que muitos dos atributos do solo seguem padrões sistemáticos e variam de acordo com a posição do solo na paisagem. Desse modo, o conhecimento da variabilidade do solo subsidia decisões sobre diferentes sistemas de manejo, além de outras aplicações como em levantamentos de solos, conhecimento e inferências nas propriedades físicas, químicas e mineralógicas dos solos, entendimento dos processos geomórficos e pedogênicos e estimativas de densidade amostral (Wilding & Dress, 1983; Trangmar et al., 1985; Burrough, 1986).

Ao se utilizar mapas de otimização de amostragem, obtêm-se informações que possibilitam entender melhor o padrão de distribuição espacial e definir diferentes zonas de manejo. Estes mapas podem ser de grande utilidade no planejamento experimental, sendo úteis também como ferramenta aos programas de agricultura de precisão, uma vez que o número de amostras a serem coletadas representa custo elevado.

A otimização amostral pelo programa Sanos 0.1, em ambas profundidades, indica que a pedoforma côncava apresenta maior variabilidade espacial para os atributos químicos do solo do que a pedoforma linear.

Ao CNPq, pela bolsa concedida a Rafael Montanari; à Fundunesp, pelo financiamento da pesquisa; à Fazenda Santa Izabel, pela parceria no desenvolvimento da pesquisa.