Variabilidade espacial de atributos de solos coesos do leste maranhense

Spatial variability of the properties of cohesive soils from eastern Maranhão, Brazil

DIVISÃO 1 - SOLO NO ESPAÇO E NO TEMPO

COMISSÃO 1.3 - PEDOMETRIA

Variabilidade espacial de atributos de solos coesos do leste maranhense(¹)

Spatial variability of the properties of cohesive soils from eastern Maranhão, Brazil

José Maria do Amaral Resende^I; José Marques Júnior^II; Marcilio Vieira Martins Filho^II; Jussara Silva Dantas^I; Diego Silva Siqueira^III; Daniel De Bortoli Teixeira^IV

^IProfessor, Centro de Ciências Agrárias e Ambientais, Universidade Federal do Maranhão, Campus de Chapadinha. BR-222, km 04, Boa Vista. CEP 65500-000 Chapadinha (MA). E-mail: [e-mail], [e-mail]

^IIProfessor Adjunto, Departamento de Solos e Adubos, Faculdade de Ciências Agrárias e Veterinárias, Universidade Estadual Paulista. Via de acesso Prof. Paulo Donato Castellane, s/n. CEP 14884-900 Jaboticabal (SP). Grupo de Pesquisa CSME (www.csme.com.br). E-mail: [e-mail], [e-mail]

^IIIPós-doutorando em Agronomia, Ciência do Solo. E-mail: [e-mail]

^IVDoutorando do Programa de Pós-Graduação em Agronomia - Produção Vegetal, FCAV/UNESP. E-mail: [e-mail], [e-mail]

ABSTRACT

A demanda nacional e internacional por mapas de atributos do solo tem aumentado. Os objetivos deste trabalho foram caracterizar e associar a variabilidade da granulometria, os atributos químicos e a cor de solos coesos a diferentes formas de paisagem, em áreas agrícolas e de vegetação nativa, por meio de técnicas geoestatísticas. Para a instalação do experimento, foram escolhidas três áreas representadas por uma área cultivada com soja e, contígua a essa, uma com cobertura vegetal remanescente de três feições de cerrado, designadas de Cerradão (A), Cerrado (B) e Campo Cerrado (C). As áreas A e B estão localizadas em pedoforma côncava e a área C, na convexa. Em cada área, foram estabelecidas malhas de amostragem com 121 pontos; nas áreas com cobertura vegetal construíram-se transeções com cinco pontos espaçados em 20 m. Os maiores alcances médios, considerando ambas as profundidades avaliadas, foram encontrados para a área A, sendo de 115 m para granulometria, 157 m para atributos químicos e 168 m para a cor do solo. A área B apresentou alcances médios de 95, 64 e 160 m para granulometria, atributos químicos e cor do solo, respectivamente. A área C, por sua vez, exibiu alcances médios de 63, 65 e 58 m para granulometria, atributos químicos e cor do solo, respectivamente. O ambiente com maior variabilidade (área C) está relacionado com locais de ocorrência de vegetação do tipo Campo Cerrado e pedoforma convexa. Esse mesmo local evidenciou a menor capacidade de resposta ao manejo da cultura de soja, evidenciado pelos índices de fertilidade do solo para essa cultura, com destaque para os baixos valores de matéria orgânica (9,20 mg kg^-1) e CTC (29,60 mmol_c dm^-3). Assim, pode-se afirmar que para o compartimento geológico estudado, os ambientes de alta variabilidade sempre estão associados a áreas com menor resposta ao manejo do solo para a cultura da soja.

Termos de indexação: geoestatística, espectroscopia de reflectância difusa, cerrado.

SUMMARY

The national and international demand for maps of soil properties has increased. The aim of this study was to characterize cohesive soils and associate their variability of particle size, chemical properties, and color to different landforms in agricultural areas and native vegetation using geostatistical techniques. To set up the experiment, we selected three areas of soil consisting of an area planted to soybeans and, contiguous to it, another area with plant cover remaining from three types of cerrado (tropical savanna), designated Cerradão (A), Cerrado (B), and Campo Cerrado (C). The areas A and B are located on a concave landform, while area C is on a convex landform. Sampling grids with 121 points were established in each area, and in areas with native vegetation, transects were constructed with five points spaced at 20 m. The largest mean ranges, considering both depths evaluated, were found for area A, with 115 m for particle size, 157 m for chemical attributes, and 168 m for soil color. Area B had mean ranges of 95, 64, and 160 m for grain size, chemical attribute, and soil color, respectively. For its part, area C exhibited the mean ranges of 63, 65, and 58 m in grain size, chemical attributes, and soil color, respectively. The environment with the greatest variability (area C) is related to locations with the occurrence of Campo Cerrado vegetation and a convex landform. This same site had the lowest responsiveness to soybean crop management, shown by soil fertility indicators for this crop, especially the low values for organic matter (9.20 mg kg^-1) and CEC (29.60 mmol_c dm^-3). Thus, it can be stated that for the geological compartment studied, the environments of high variability will always be associated with areas with lower response to soil management for soybean.

Index terms: geostatistics, diffuse reflectance spectroscopy, savanna.

INTRODUÇÃO

A expansão da fronteira agrícola, na Região Nordeste do Brasil, incorporou ao sistema produtivo solos de fertilidade marginal com alto grau de intemperização que apresentam atributos físicos que limitam o crescimento das raízes das plantas cultivadas e facilitam o desenvolvimento do processo erosivo. Esse é o caso dos solos coesos, desenvolvidos a partir de material originário da Formação Barreiras. Esses solos, mesmo quando pertencentes à mesma classe taxonômica, podem apresentar variações em seus atributos, tanto horizontal, ao longo da paisagem, quanto verticalmente, ao longo do perfil do solo. Essas variações podem ser resultantes dos fatores responsáveis pela formação do solo, pois tanto o material de origem como a sua intemperização não ocorrem homogeneamente (Buol et al., 1997).

A caracterização da variabilidade horizontal dos atributos do solo representa importante atividade de pesquisa na compreensão de processos e fatores inerentes ao solo (Parkin, 1993). Os resultados dessas pesquisas são importantes para nortear futuros projetos e subsidiar a tomada de decisão, para uso e ocupação do solo, em locais de variabilidade desconhecida, como os solos coesos da Formação Barreiras. Esses solos ocupam cerca de 200.000 km² no território nacional. A caracterização detalhada dessas variações pode indicar a ocorrência de processos capazes de influenciar na produtividade diferenciada das culturas; em consequência, possibilita a identificação de locais específicos de manejo do solo (Sanchez et al., 2012).

As técnicas de análises geoestatísticas destacam-se nesse cenário por permitir a identificação e o delineamento dessas relações espaciais (Vieira, 2000). Para a condução dessas análises, assume-se que determinada propriedade numa posição definida na área de estudo varia de acordo com a distância de separação entre as amostras vizinhas. Essa continuidade geográfica característica das variáveis é identificada ao observar que os valores das variáveis tornam-se mais similares à medida que a distância de separação entre elas diminui (Isaaks & Srivastava, 1989).

Na área de ciências agrárias, as análises geoestatísticas têm sido utilizadas com grande êxito na caracterização de atributos físicos, químicos e mineralógicos do solo (Camargo et al., 2008), adsorção de fósforo (Barbieri et al., 2009), estudo da emissão de CO₂ (Teixeira et al., 2012), aplicação de insumos em taxa variada (Barbieri et al., 2008), perda de solo e nutrientes por erosão (Sanchez et al., 2009), cultivo de citros (Siqueira et al., 2010), produção de café (Sanchez et al., 2012), entre outras. Recentemente, as análises geoestatísticas vêm sendo associadas a informações da paisagem, visando à compreensão mais precisa dos padrões de variabilidade dos atributos do solo e como esses podem ser utilizados para explicar a produção agrícola (Siqueira et al., 2010). Dentre essas informações, podem-se citar a utilização da topografia (Noorbakhsh et al., 2008); os modelos de paisagem, como superfícies geomórficas (Campos et al., 2007); as formas de relevo (Souza et al., 2009); e as pequenas variações do relevo (Camargo et al., 2008; Marques Jr. et al., 2014).

Nesse contexto, os objetivos deste trabalho foram caracterizar e associar a variabilidade de atributos da granulometria, bem como os químicos e os da cor de solos coesos a diferentes formas de paisagem, em àreas agrícolas e de vegetação nativa, por meio de técnicas geoestatísticas.

MATERIAL E MÉTODOS

Localização, caracterização da área e planejamento amostral

O estudo foi desenvolvido no município de Brejo, MA, na Fazenda Typuana, localizado na mesorregião leste do Maranhão, microrregião de Chapadinha, MA, nas coordenadas geográficas de 03^º 36' S e 42^º 52' W (Figura 1a). O clima, segundo a classificação climática de Köppen-Gerger, é do tipo Aw, clima tropical com estação seca de inverno. A estação chuvosa está concentrada entre os meses de novembro e maio, com precipitação pluvial média de 1.835 mm, temperatura média anual acima de 27 ^ºC e umidade relativa anual entre 73 e 79 % (Maranhão, 2002). A Formação Barreiras representa uma cobertura sedimentar terrígena continental, de idade pliocênica, depositada por sistemas fluviais entrelaçados, associados a leques aluviais. As faces de sistemas fluviais entrelaçados apresentam depósitos de granulometria variada com cascalhos e areias grossas e finas, de coloração creme amarelada, com intercalações de microclastos de argila síltica, indicativo de ambiente de sedimentação calmo, como de planície aluvial (Alheiros et al., 1988).

Para a instalação do experimento, foram selecionadas três áreas de solo denominadas A, B e C, com base na vegetação nativa. Cada compartimento era composto por dois ambientes, um ambiente de produção, cultivado com soja e outro, contíguo a esse, com cobertura vegetal remanescente de três feições de cerrado, designadas de Cerradão, Cerrado e Campo Cerrado.

O modelo digital de elevação foi utilizado para verificar a localização das áreas na paisagem (Figura 1b). Esse foi gerado com base nas informações da missão Shuttle Radar Topography Mission (SRTM), com resolução horizontal de 90 m e precisão vertical da ordem de 5 m. Para detectar possíveis valores ruidosos nas informações SRTM, foi aplicado o filtro de mediana, subtração do resultado do filtro pelo dado original, retirada dos valores com variação de 10 m decorrente da falta de informação e interpolação pelo método Topogrid. Também foram acrescentadas às informações SRTM, as informações de campo da altitude de 120 pontos das três áreas estudadas.

Na área A, apresenta-se vegetação remanescente de Cerradão, que está localizada em pedoforma côncava, relevo suave ondulado (Figura 1b). Esse ambiente de produção vem sendo cultivado há 10 anos com soja. Na área B, encontra-se a vegetação remanescente Cerrado, que está situada em pedoforma côncava, relevo suave ondulado. Esse ambiente de produção vem sendo cultivado há seis anos com soja. Já na área C, identifica-se a vegetação remanescente de Campo Cerrado, localizada em pedoforma convexa, relevo suave ondulado, sendo cultivada há quatro anos com soja.

Para a classificação dos solos, segundo o Sistema Brasileiro de Classificação dos Solos (Embrapa, 2006), foram abertas e descritas trincheiras em cada uma das áreas no ambiente de vegetação. O solo das três áreas, desenvolvido a partir de material da Formação Barreiras, foi classificado como Argissolo Amarelo distrocoeso típico.

Nos ambientes de produção das três áreas, é feita a correção da acidez do solo toda vez que a necessidade de calagem, calculada pelo método do V %, é maior que 1 t ha^-1 de calcário com PRNT de 70 %. A adubação de manutenção tem sido a mesma para as três áreas que anualmente tem recebido 16 kg ha^-1 de N, 120 kg ha^-1 de P₂O₅ e 60 kg ha^-1 de K₂O.

Nos ambientes de produção das áreas A e B, foram estabelecidas malhas de amostragem quadradas com espaçamento regular de 25 m, totalizando 121 pontos de amostragem em uma área de 6,2 ha (Figura 1c); já no da área C, estabeleceu-se malha de amostragem quadrada com espaçamento regular de 15 m, totalizando 121 pontos de amostragem em uma área de 2,2 ha (Figura 1d). Para fins de comparação da variabilidade dos atributos do solo em ambientes natural e de produção, foram construídas transeções de 100 m nas vegetações remanescentes. Foram coletados nas transeções cinco pontos espaçados em 20 m. As amostras das malhas e transeções foram coletadas nas profundidades de 0,0-0,2 e 0,2-0,4 m.

Análises laboratoriais

A análise granulométrica foi realizada utilizando uma solução de NaOH 0,1 mol L^-1 como dispersante químico e agitação com aparato de baixa rotação por 16 h. A argila foi determinada pelo método da pipeta, a areia foi separada por tamisação e o silte foi calculado por diferença (Embrapa, 1997). Para as análises de fertilidade do solo, utilizaram-se os métodos propostos por Raij et al. (2001). Com base nesses resultados, foram calculados os valores da soma de bases (SB), da capacidade de troca catiônica potencial (CTC) e da porcentagem de saturação por bases (V%). O fósforo remanescente foi determinado de acordo com Alvarez V. et al. (2000). Essas análises foram realizadas para todas as amostras coletadas nas áreas estudadas.

Avaliação da espectroscopia de reflectância difusa e determinação da cor

Para obtenção dos espectros de reflectância difusa, as avaliações foram feitas com o sensor de laboratório Lambda 950. Foi moído aproximadamente 0,5 g de TFSA em ágata até alcançar a coloração constante. O material de solo foi colocado em um porta-amostra com um espaço cilíndrico de 16 mm. Os valores de reflectância foram determinados em espectrofotômetro equipado com esfera integradora de 80 mm, a cada 1 nm, com um tempo de integração de 0,2 s, fazendo uma varredura no intervalo de 380 a 780 nm. Após se obterem os espectros de reflectância difusa das amostras de solo, foram determinados os valores de triestímulo XYZ definidos pela Comisión Internacional de L'Eclairage. A partir das coordenadas XYZ, foram deduzidos os valores Munsell de matiz, croma e valor, utilizando o programa Munsell Conversion versão 6.4, conforme Barrón et al. (2000) e Viscarra Rossel (2011). Essas determinações foram realizadas para todas as amostras coletadas nas áreas estudadas.

Análises dos dados

Foram calculados a média, o coeficiente de variação e o desvio-padrão. Para a avaliação da diferença dos valores médios entre as áreas, identificadas com base na vegetação nativa, foi aplicado o teste de Tukey a 5 %. Nas comparações das médias, foram considerados os seguintes fatores no esquema fatorial: áreas de solos com três níveis (A, B, C), ambiente com dois níveis, ambiente de produção (AP) e ambiente de vegetação (VN), o que resultaram seis tratamentos: AAP, AVN, BAP, BVN, CAP e CVN. Dessa maneira, para a aplicação do teste de Tukey, cada ponto coletado em cada área foi interpretado como sendo uma pseudorrepetição. O teste de Tukey leva em conta o desbalanceamento dessas pseudorrepetições.

A modelagem do variograma experimental, seguindo os princípios estabelecidos pela hipótese intrínseca (Isaaks & Srivastava, 1989), foi realizada visando a captação da variabilidade espacial da granulometria, os atributos químicos e a cor do solo. O variograma experimental foi determinado por meio do cálculo da variância em razão da distância de separação entre amostras (Equação 1):

em que é a variância experimental para uma distância de separação h, N(h) é o número de pares experimentais de observações Z(x_i) e Z(x_i+h) são separados pela distância h. A estacionariedade necessária ao uso da geoestatística foi avaliada por meio de análise de tendência, utilizando regressões lineares para os eixos X, Y e suas interações. Assim, pela diferença entre o valor medido e o valor da superfície do polinômio ajustada, obtiveram-se os valores do resíduo que foi utilizado para a construção do variograma. A escolha do melhor modelo ajustado aos variogramas baseou-se na soma de quadrado dos resíduos (SQR) e no coeficiente de determinação (R²).

Após a construção dos variogramas, foi utilizada a técnica da krigagem ordinária para a interpolação de valores em locais não amostrados. Essa técnica baseia-se em uma média móvel ponderada das amostras vizinhas (Equação 2); os pesos (λ_i) de cada vizinho são determinados utilizando o modelo de variograma ajustado, resultando em uma estimativa de variância mínima:

em que é o valor estimado da propriedade no ponto 0; N, o número de valores utilizados na predição; λ_i, a ponderação associada a cada valor; e z(x_i), o valor observado no ponto i.

RESULTADOS E DISCUSSÕES

Os resultados obtidos no ambiente de vegetação evidenciam que os solos das áreas A e B não apresentaram diferenças significativas entre si quanto à granulometria, porém diferem, na maioria das frações granulométricas do solo, do ambiente C, que apresenta maiores teores de argila e silte e menores teores de areias nas duas profundidades estudadas (Quadro 1).

Embora esses solos tenham a mesma denominação taxonômica (Argissolo Amarelo distrocoeso típico) e algumas semelhanças morfológicas, possuem diferentes padrões de variabilidade, tanto em ambiente de produção como em de vegetação, fato esse observado por meio dos diferentes valores de coeficiente de variação (CV). Comparando-se os valores dos CV das frações granulométricas do ambiente de vegetação, com os valores de CV do ambiente de produção, na profundidade de 0,0-0,2 m, pode-se afirmar que o solo que apresentou a maior variabilidade de seus atributos em razão das práticas agrícolas, consideradas um dos fatores responsáveis pela variabilidade granulométrica, foi o solo da área A, seguido dos das áreas B e C, indicando que o solo menos impactado pelo manejo foi o da área C. Essas mesmas observações são válidas para a profundidade de 0,2-0,4 m.

Os solos das três áreas estudadas, no ambiente de vegetação, não diferem entre si quanto aos atributos químicos, na profundidade de 0,0-0,2 m, sendo pobres em Ca, Mg, K e P, apresentando altos valores de acidez ativa determinada em CaCl₂, baixos valores de CTC potencial, soma de bases e porcentagem de saturação por bases, de acordo com Raij et al. (1997) (Quadro 2). Essas condições apresentam-se como desfavoráveis ao ótimo crescimento e desenvolvimento das plantas cultivadas. Tais valores são ainda mais baixos na profundidade de 0,2-0,4 m (Quadro 2).

Nesse ambiente, os solos apresentam baixos valores de adsorção de P (Alvarez V. et al., 2000), expresso na forma de P remanescente, não diferindo significativamente entre eles, o que pode ser creditado aos baixos teores de argila apresentados pelos solos estudados. Também não houve diferenças significativas entre os teores de matéria orgânica nos solos das três áreas de vegetação estudadas, que variaram de 12,8 a 18,8 g kg^-1, na profundidade de 0,0-0,2 m (Quadro 2).

No ambiente de produção, pode-se verificar que os solos responderam de forma diferenciada à adubação, sendo as maiores respostas encontradas no solo da área A; e as menores respostas, no da área C. Embora o solo da área A seja aquele manejado há mais tempo, isso por si só não justifica os menores aumentos nos atributos químicos no solo da área C, uma vez que vem recebendo doses maciças de corretivos e fertilizantes minerais. Em um solo com cultivo intenso de cana-de-açúcar por mais de 20 anos, Sanchez et al. (2009) também verificaram que os maiores valores médios dos atributos químicos ocorreram na pedoforma côncava. Esse comportamento pode ser atribuído às condições de drenagem e características convergentes da inclinação dos declives responsáveis pela acumulação dos nutrientes (Resende et al., 2007).

A cobertura vegetal também reflete a capacidade de resposta do solo à adubação. O solo que mais responde à adubação apresenta a vegetação denominada Cerradão, e o que menos responde, a vegetação denominada Campo Cerrado.

Os elementos constituintes de cor do solo, matiz, valor e croma, determinados pelo método da reflectância difusa, não apresentaram diferenças significativas entre os ambientes de vegetação e de produção, para solos de uma mesma área, evidenciando que a cor do solo é um atributo morfológico não influenciado pelas práticas de manejo (Quadro 3).

Quando se comparam os solos das diversas áreas entre si, verifica-se que os das áreas A e B, que apresentam coloração mais amarelada, não diferem entre si, mas diferem do solo da área C, que exibe coloração com um matiz mais avermelhado. Os solos das áreas A e B estão situados em pedoformas côncavas do relevo, enquanto o solo da área C situa-se em forma convexa. A forma côncava favorece o armazenamento de água, tornando o ambiente redutor. Assim, o Fe³⁺, que confere a cor vermelha do solo, é reduzido a Fe²⁺, mais solúvel, que é retirado do sistema, ficando o solo com coloração mais amarelada (Resende et al., 2007).

O CV não tem sido considerado como bom indicador da variabilidade espacial dos atributos do solo, quando utilizado individualmente, pois em alguns casos pode mascarar a real interpretação da heterogeneidade espacial sobre os valores médios. Nesse caso, é recomendada a utilização da análise geoestatística (Sampaio et al., 2011).

Aos variogramas experimentais dos valores das granulometrias de todos os solos, foi ajustado o modelo esférico (Quadro 4). Esse modelo tem sido o mais utilizado para descrever o comportamento de atributos de plantas e solo por diversos autores (Cambardella et al., 1994; Barbieri et al., 2008; Gomes et al., 2008; Schaffrath et al., 2008; Sanchez et al., 2009; Souza et al., 2009; Silva Neto et al., 2011), e normalmente ajusta-se aos dados cuja distribuição apresenta mudanças abruptas.

Os solos apresentaram dependência espacial para as frações granulométricas nas duas profundidades estudadas, apresentando Grau de Dependência Espacial (GDE) de moderada a forte, segundo classificação de Cambardella et al. (1994), a qual permite classificar o GDE em dependência espacial forte, quando 25 %; dependência espacial moderada, quando 25 % < GDE 75 %; e dependência espacial fraca, quando GDE > 75 %. Na profundidade de 0,0-0,2 m, os valores médios do GDE foram de 11,65; 58,31; e 46,66 %, respectivamente nos solos das áreas A, B e C; e na profundidade de 0,2-0,4 m, de 15,66; 27,51 e 31,23 %, respectivamente nos solos das áreas A, B e C.

O alcance representa a distância em que os pontos amostrais estão correlacionados entre si, ou seja, os pontos localizados numa área de raio igual ao alcance são mais homogêneos entre si. Diferentes valores de alcance foram encontrados para as várias frações granulométricas, nos diferentes solos e nas profundidades estudadas. Na profundidade de 0,0-0,2 m, o solo da área A apresentou os maiores alcances em todas as frações granulométricas estudadas, quando comparado aos solos das demais áreas. Os valores médios foram de 117, 72 e 66 m, para os solos das áreas A, B e C, respectivamente. Na profundidade de 0,2-0,4 m, os alcances médios foram de 113, 119 e 59 m, para os solos das áreas A, B e C, respectivamente.

A análise dos valores de alcance (Quadro 4) e dos de CV (Quadro 1) permite melhor compreensão do comportamento da variabilidade dos atributos do solo. Assim, os resultados médios de alcance dos variogramas contradizem a variabilidade inferida unicamente por meio da estatística descritiva, demonstrando que a menor variabilidade ocorre no solo da área A e a maior, no da área C. Como o solo da área A é o manejado há mais tempo, bem como o da área C há menos tempo, fica assim demonstrada a contribuição do manejo para a homogeneização da camada superficial.

Todos os valores dos atributos químicos dos solos, das três áreas, quando houve dependência espacial, foram ajustados ao modelo esférico (Quadros 5 e 6). Ajuste semelhante também foi obtido por Silva Neto et al. (2011). A maioria dos atributos químicos dos solos estudados apresentou dependência espacial moderada a forte. Na profundidade de 0,0-0,2 m, à exceção de K, Mg e V% no solo da área B, que evidenciou Efeito Pepita Puro (EPP), todos os demais atributos químicos do solo apresentaram GDE moderada a forte, sendo o GDE moderado predominante na maioria dos atributos. Os valores médios de GDE foram de 39,14; 55,42; e 54,83 % para os solos das áreas A, B e C, respectivamente.

Para a profundidade de 0,20 a 0,40 m, foi observada ausência de dependência espacial, caracterizada pelo EPP, para os atributos: V% (área A); P, K, Ca, Mg e SB (área B); e Ca (área C). Sob essas condições, a avaliação desses atributos é feita exclusivamente pela estatística descritiva. Nessa profundidade, os valores médios do GDE foram de 36,68; 45,74; e 39,91 % para os solos das áreas A, B e C, respectivamente.

A forte dependência espacial dos atributos do solo é atribuída a fatores intrínsecos, enquanto a dependência espacial moderada, a fatores extrínsecos (Carvalho et al., 2003). Assim, Silva Neto et al. (2011) creditaram a forte dependência espacial encontrada para os atributos químicos do solo a qualquer um dos fatores de formação desse, principalmente ao material de origem e relevo, enquanto a dependência espacial moderada seria em razão da sua homogeneização, que as adubações e calagens proporcionam ao solo.

Os atributos químicos do solo apresentaram diferentes alcances nas áreas e profundidades estudadas; entretanto, verificou-se que o solo da área A, com exceção dos atributos P-rem e V% na camada superficial, apresentou os maiores valores de alcance para todos os atributos, nas duas profundidades. O alcance médio para a camada superficial foi de 155, 66 e 82 m para os solos das áreas A, B e C, respectivamente.

Na camada de 0,2-0,4 m, as médias dos alcances foram de 161, 80 e 55 m para os solos das áreas A, B e C, respectivamente. Também os valores dos alcances evidenciaram maior homogeneidade do solo da área A, em relação aos demais. Observou-se que, tanto para os atributos da granulometria quanto para os químicos, nas duas profundidades estudadas, os valores médios dos alcances são sempre maiores nos solos das áreas A e B, os quais estão localizados em áreas côncavas, e menores no solo da área C, localizado em área convexa, mostrando a influência do relevo no comportamento desses atributos. Resultados semelhantes foram encontrados por Barbieri et al. (2008). Para Marques Jr & Lepsch (2000), o entendimento das causas da variabilidade do solo está intimamente relacionado ao conhecimento dos processos do solo que operam em locais específicos. Tais processos estão relacionados à água, sendo o relevo o principal controlador de sua intensidade e de seu fluxo em determinada área.

Segundo McBratney & Webster (1983), o conhecimento do alcance é importante na definição da ótima intensidade de amostragem, visando a reduzir o erro-padrão da média, além de aumentar a representatividade da amostra. Para Montanari et al. (2005), 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, devendo-se considerar a paisagem como fator relevante em estratégia de amostragem de solo.

Todos os constituintes da cor do solo apresentaram GDE, sendo considerado forte para o atributo matiz no solo da área A; fraco, para o da área B; e moderado, para os demais solos (Quadro 7).

Os maiores alcances para matiz, valor e croma foram observados nos solos das áreas A e B, localizados em pedoforma côncava e originalmente sob vegetação de Cerradão e Cerrado, respectivamente. Os menores alcances para todos os constituintes da cor foram observados no solo da área C, localizado em pedoforma convexa e originalmente sob vegetação de Campo Cerrado. O alcance médio dos constituintes da cor do solo foi de 168, 160 e 56 m para os solos das áreas A, B e C, respectivamente.

Os padrões espaciais dos principais atributos estudados obtidos por meio da krigagem ordinária são apresentados na figura 2. Esses atributos foram selecionados com base na sua relevância para a caracterização das diferentes áreas. Os padrões espaciais dos atributos da área A apresentaram-se mais homogêneos, em comparação aos das áreas B e C. Esse fato é observado em razão do maior espaçamento entre as classes de variabilidade da área A em relação às demais. Essa maior homogeneidade é reflexo dos maiores alcances obtidos nos variogramas dos atributos da área A, que, por sua vez, podem ser explicados pelo maior tempo de cultivo dessa área em relação às demais, associada à sua localização em pedoforma côncava. Segundo Souza et al. (2003), a forma côncava possui menor variabilidade em relação à forma linear, e essa em relação à convexa. O padrão espacial do matiz da área B apresentou a maior erraticidade entre os atributos estudados, caracterizada pela intensa irregularidade na divisão entre as classes de variabilidade. Esse fato resulta da alta contribuição do efeito pepita para o modelo de dependência espacial, caracterizado pelo maior valor de GDE (80,17 %), entre os atributos estudados. Segundo Souza et al. (2006), mesmo que pequenas, as variações da forma da paisagem influem significativamente na dinâmica do fluxo de água superficial e subsuperficial, influenciando na variabilidade dos atributos do solo.

A comparação entre os padrões espaciais dos atributos das diferentes áreas evidencia a maior semelhança entre os atributos da área A. Esses apresentam altas intensidades de coeficientes de correlação (p<0,01), variando de 0,74 (croma × matiz) a 0,95 (valor × MO), sendo essa última correlação negativa. Nas áreas B e C, os coeficientes de correlação (p<0,05) variaram de 0,04 (argila × Prem) a 0,85 (valor × croma) e de 0,03 (matiz × CTC) a 0,79 (valor × croma), respectivamente. Essa maior semelhança deve-se ao padrão mais homogêneo da área A em relação às demais, tornando as linhas de isovalores mais coincidentes.

CONCLUSÕES

1. O relevo, associado à vegetação nativa, pode ser utilizado no mapeamento de áreas mínimas de manejo. Os locais de ocorrência da vegetação denominada Campo Cerrado, em pedoforma convexa, apresentam maior variabilidade espacial e menor potencial de resposta ao manejo da cultura de soja.

2. Em solos coesos, onde a variação do relevo não é perceptível no campo, e há ausência de vegetação nativa, a cor do solo pode ser indicador dos diferentes padrões de variabilidade da granulometria e dos atributos químicos do solo, podendo auxiliar na definição de áreas mínimas de manejo para cultura da soja.

LITERATURA CITADA

Recebido para publicação em 16 de junho de 2013

Aprovado em 28 de abril de 2014

• ALHEIROS, M.M.; LIMA FILHO, M.F.; MONTEIRO, F.A.J. & OLIVEIRA FILHO, J.S. Sistemas deposicionais na Formação Barreiras no Nordeste Oriental. In: CONGRESSO BRASILEIRO DE GEOLOGIA, 35., Belém, 1988. Anais... Belém, Sociedade Brasileira de Geologia, 1988. p.753-760.
• ALVAREZ V., V.H.; NOVAIS, R.F.; DIAS, L.E. & OLIVEIRA, J.A. Determinação e uso do fósforo remanescente. B. Inf. Soc. Bras. Ci. Solo, 25:27-32, 2000.
• BARBIERI, D.M.; MARQUES JR, J.; ALLEONI, L.R.F.; GARBUIO, F.J. & CAMARGO, L.A. Hillslope curvature, clay mineralogy, and phosphorus adsorption in an Alfisol cultivated with sugarcane. Sci. Agric., 66:819-826, 2009.
• BARBIERI, D.M.; MARQUES JÚNIOR, J.M. & PEREIRA, G.T. Variabilidade espacial de atributos químicos de um Latossolo para aplicação de insumos à taxa variável em diferentes formas de relevo. Eng. Agríc., 28:645-653, 2008.
• BARRÓN, V.; MELLO, J.W.V. & TORRENT, J. Caracterização de óxidos de ferro em solos por espectroscopia de reflectância difusa. In: NOVAIS, R.F.; ALVAREZ V., V.H. & SCHAEFER, C.E.G.R., eds. Tópicos em ciência do solo. Viçosa, Sociedade Brasileira de Ciência do solo, 2000. v.1, p.139-162.
• BUOL, S.W.; HOLE, F.D.; McCRACKEN, R.J. & SOUTHARD, R.J. Soil Genesis and Classification. 4.ed. Iowa, Wiley-Blackwell, 1997. 527p.
• CAMARGO, L.A.; MARQUES JÚNIOR, J.; PEREIRA, G.T. & HORVAT, R.A. Variabilidade espacial de atributos mineralógicos de um Latossolo sob diferentes formas de relevo. I - Mineralogia da fração argila. R. Bras. Ci. Solo, 32:2269-2277, 2008.
• CAMBARDELLA, C.A.; MOOMAN, T.B.; NOVAK, J.M.; PARKIN, T.B.; KARLEM, D.L.; TURVO, R F. & KANOPA, A.E. Field scale variability of soil properties in central Iowa soil. Soil Sci. Soc. Am. J., 47:1501-1511, 1994.
• CAMPOS, M.C.C.; MARQUES JÚNIOR, J.; PEREIRA, G.T.; MONTANARI, R. & CAMARGO, L.A. Relações solo-paisagem em uma litossequência arenito-basalto na região de Pereira Barreto - SP. R. Bras. Ci. Solo, 31:519-539, 2007.
• CARVALHO, M.P.; TAKEDA, E.Y. & FREDDI, O.S. Variabilidade espacial de atributos de um solo sob videira em Vitória Brasil (SP). R. Bras. Ci. Solo, 27:695-70, 2003.
• EMPRESA BRASILEIRA DE PESQUISA AGROPECUÁRIA - EMBRAPA. Centro Nacional de Pesquisa de Solos. Manual de métodos de análise de solo. 2.ed. Rio de Janeiro, 1997. 212p.
• EMPRESA BRASILEIRA DE PESQUISA AGROPECUÁRIA - EMBRAPA. Centro Nacional de Pesquisa de Solos. Sistema brasileiro de classificação de solos. 2 ed. Rio de Janeiro, 2006. 306p.
• GOMES, J.B.V.; BOLFE, E.L.; CURI, N.; FONTES, H.R.; BARRETO, A.C. & VIANA, R.D. Variabilidade espacial de atributos dos solos em unidades de manejo em área-piloto de produção integrada de coco. R. Bras. Ci. Solo, 32:2471-2482, 2008.
• MARANHÃO. Governo do Estado. Atlas do Maranhão. 2.ed. São Luís, Gerência de Planejamento e Desenvolvimento Econômico, 2002. 39p.
• ISAAKS, E.H. & SRIVASTAVA, R.M. Applied geostatistics. Nova York, Oxford University Press, 1989. 561p.
• MARQUES JR, J. & LEPSCH, I.F. Depósitos superficiais neocenozoicos em superfícies geomórficas e solos em Monte Alto - SP. Geociência, 19:90-106, 2000.
• MARQUES JR., J.; SIQUEIRA, D.S.; CAMARGO, L.A.; TEIXEIRA, D.D.B.; BARRÓN, V. & TORRENT, J. Magnetic susceptibility and diffuse reflectance spectroscopy to characterize the spatial variability of soil properties in a Brazilian Haplustalf. Geoderma, 219:63-71, 2014.
• McBRATNEY, A.B. & WEBSTER, R. How many observations are needed for regional estimation of soil properties. Soil Sci. Soc. Am. J., 135:177-183, 1983.
• MONTANARI, R.; MARQUES JÚNIOR, J.; PEREIRA, G.T. & SOUZA, Z.M. Forma da paisagem como critério para otimização amostral de Latossolos sob cultivo de cana-de-açúcar. Pesq. Agropec. Bras., 40:69-77, 2005.
• NOORBAKHSH, S.; SCHOENAU, J.; SI, B.; ZELEKE, T. & QIAN, P. Soil properties, yield, and landscape relationships in South-Central Saskatchewan, Canada. J. Plant Nutr., 31:539-556, 2008.
• RAIJ, B.van; ANDRADE, J.C.; CANTARELA, H. & QUAGGIO, J.A. Análise química para avaliação da fertilidade de solos tropicais. Campinas, Instituto Agronômico de Campinas, 2001. 285p.
• RAIJ, B.van; CANTARELA, 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, 1997. 285p. (Boletim Técnico, 100)
• RESENDE, M.; CURI, N.; REZENDE, S.B. & CORRÊA, G.F. Pedologia: Base para distinção de ambientes. 5.ed. Lavras, Universidade Federal de Lavras, 2007. 322p.
• SAMPAIO, M.S.; ALVES, M.C.; SILVA, F.M.; POZZA, E.A. & OLIVEIRA, M.S. Avaliação do comportamento da variabilidade espacial do fósforo remanescente no solo de lavoura cafeeira. In: SIMPÓSIO BRASILEIRO DE SENSORIAMENTO REMOTO, 15., Curitiba, 2011. Anais... Curitiba, Sociedade Brasileira de Sensoriamento Remoto, 2011. p.9203-9210.
• SANCHEZ, R.B.; MARQUES JÚNIOR., J.; PEREIRA, G.T.; SOUZA, Z.M. & MARTINS FILHO, M.V. Variabilidade espacial de atributos do solo e de fatores de erosão em diferentes pedoformas. Bragantia, 68:873-884, 2009.
• SANCHEZ, R.B.; MARQUES JUNIOR, J.; PEREIRA, G.T.; BARACAT NETO, J.; SIQUEIRA, D.S. & SOUZA, Z.M. Mapeamento das formas do relevo para estimativa de custos de fertilização em cana-de-açúcar. Eng. Agríc., 32:280-292, 2012.
• SCHAFFRATH, V.R.; TORMENA, C.A.; FIDALSKI, J. & GONÇALVES, A.C.A. Variabilidade e correlação espacial de propriedades físicas de solo sob plantio direto e preparo convencional. R. Bras. Ci. Solo, 32:1369-1377, 2008.
• SILVA NETO, S.P.; SANTOS, A.C.; LEITE, R.L.L.; DIM, V.P.; CRUZ, R.S.; PEDRICO, A. & NEVES NETO, D.N. Análise espacial de parâmetros da fertilidade do solo em região de ecótono sob diferentes usos e manejos. Semina: Ci. Agron., 32:541-552, 2011.
• SIQUEIRA, D.S.; MARQUES JR., J. & PEREIRA, G.T. The use of landforms to predict the variability of soil and orange attributes. Geoderma, 155:55-66, 2010.
• SOUZA, C.K.; MARQUES JÚNIOR, J.; MARTINS FILHO, M.V. & PEREIRA, G.T. Influência do relevo e erosão na variabilidade espacial de um Latossolo em Jaboticabal (SP). R. Bras. Ci. Solo, 27:1067-1074, 2003.
• SOUZA, Z.M.; MARQUES JUNIOR, J.; PEREIRA, G.T. & BARBIERI, D.M. Small relief shape variations influence spatial variability of soil chemical attributes. Sci. Agric., 63:161-168, 2006.
• SOUZA, Z.M.; MARQUES JÚNIOR. J. & PEREIRA, G.T. Spatial variability of the physical and mineralogical properties of the soil from the areas with variation in landscape shapes. Braz. Arch. Biol. Technol., 52:305-316, 2009.
• PARKIN, T.B. Spatial variability of microbial processes in soil - a review. J. Environ. Qual., 22:409-417, 1993.
• TEIXEIRA, D.D.B.; BICALHO, E.S.; PANOSSO, A.R.; PERILLO, L.I.; IAMAGUTI, J.L.; PEREIRA, G.T. & LA SCALA Jr., N. Uncertainties in the prediction of spatial variability of soil CO₂ emissions and related properties. R. Bras. Ci. Solo, 36:1466-1475, 2012.
• VIEIRA, S.R. Uso de geoestatística em estudos de variabilidade espacial de propriedades do solo. In: NOVAIS, R.F.; ALVAREZ V., V.H. & SCHAEFER, C.E.G.R., eds. Tópicos em ciência do solo. Viçosa, MG, Sociedade Brasileira de Ciência do Solo, 2000. v.1, p.1-54.
• VISCARRA ROSSEL, R.A. Fine-resolution multiscale mapping of clay minerals in Australian soils measured with near infrared spectra. J. Geophys. Res., 116:1-15, 2011.

(1

) Parte da Tese de Doutorado do primeiro autor.

Datas de Publicação

• Publicação nesta coleção
  02 Out 2014
• Data do Fascículo
  Ago 2014

Histórico

• Aceito
  28 Abr 2014
• Recebido
  16 Jun 2013