Acessibilidade / Reportar erro

Distribuição e incerteza da acidez de um latossolo vermelho-amarelo húmico sob cultivo de café

Distribution and uncertainty of acidity of an oxisol under coffee cultivation

Resumos

Com o desenvolvimento da agricultura, a utilização de novas alternativas na avaliação das propriedades que influenciam o rendimento das plantas tornou-se indispensável para melhor manejar o sistema agrícola. O objetivo deste trabalho foi utilizar uma ponderação por meio da modelagem fuzzy para estudar, com base em atributos químicos do solo, a propagação de incertezas da acidez de um Latossolo Vermelho-Amarelo húmico cultivado com café arábica. As amostragens de solo foram realizadas na profundidade de 0-20 cm, em uma malha totalizando 50 pontos. Os atributos avaliados foram: pH em H2O, acidez trocável (Al3+), acidez potencial (H + Al) e saturação por Al (m). Os dados foram analisados pela estatística descritiva e pela geoestatística. Utilizou-se um sistema de classificação fuzzy e os atributos descritos para inferir sobre a acidez do solo. A lógica fuzzy, pelo algoritmo utilizado, teve um bom desempenho na caracterização e no mapeamento das incertezas da acidez do solo em questão.

geoestatística; krigagem; lógica nebulosa; manejo do solo


With the development of agriculture, the use of new tools for the evaluation of the properties related to plant yields have become essential to improve the management of agricultural systems. The purpose of this study was to use a fuzzy model to investigate, based on soil chemical properties, the propagation of uncertainty of the acidity of a humic Red-Yellow Latosol (Oxisol) under Arabic coffee plantation. The soil was sampled at a depth of 0 - 0.2 m, in a grid, totaling 50 sampling points. The following attributes were determined: pH in H2O, exchangeable acidity (Al3+), potential acidity (H + Al) and Al saturation (m). The data were analyzed by descriptive statistics and geostatistics. A fuzzy classification system was used based on the attributes described to infer on the soil acidity. Fuzzy logic, with the algorithm used, was efficient to characterize and map the uncertainties of acidity of the soil under study.

geostatistics; kriging; hazy logic; soil management


NOTA

Distribuição e incerteza da acidez de um latossolo vermelho-amarelo húmico sob cultivo de café

Distribution and uncertainty of acidity of an oxisol under coffee cultivation

Samuel de Assis SilvaI; Julião Soares de Souza LimaII; Gustavo Soares de SouzaI; Alexandre Candido XavierIII

IMestre em Produção Vegetal, Universidade Federal do Espírito Santo - UFES. Caixa Postal 16, CEP 29500-000 Alegre (ES). E-mail: samuel-assis@hotmail.com

IIProfessor Associado, UFES. E-mail: juliaosslima@cca.ufes.br

IIIProfessor Adjunto, UFES. E-mail: xavier@cca.ufes.br

RESUMO

Com o desenvolvimento da agricultura, a utilização de novas alternativas na avaliação das propriedades que influenciam o rendimento das plantas tornou-se indispensável para melhor manejar o sistema agrícola. O objetivo deste trabalho foi utilizar uma ponderação por meio da modelagem fuzzy para estudar, com base em atributos químicos do solo, a propagação de incertezas da acidez de um Latossolo Vermelho-Amarelo húmico cultivado com café arábica. As amostragens de solo foram realizadas na profundidade de 0-20 cm, em uma malha totalizando 50 pontos. Os atributos avaliados foram: pH em H2O, acidez trocável (Al3+), acidez potencial (H + Al) e saturação por Al (m). Os dados foram analisados pela estatística descritiva e pela geoestatística. Utilizou-se um sistema de classificação fuzzy e os atributos descritos para inferir sobre a acidez do solo. A lógica fuzzy, pelo algoritmo utilizado, teve um bom desempenho na caracterização e no mapeamento das incertezas da acidez do solo em questão.

Termos de indexação: geoestatística; krigagem; lógica nebulosa; manejo do solo.

SUMMARY

With the development of agriculture, the use of new tools for the evaluation of the properties related to plant yields have become essential to improve the management of agricultural systems. The purpose of this study was to use a fuzzy model to investigate, based on soil chemical properties, the propagation of uncertainty of the acidity of a humic Red-Yellow Latosol (Oxisol) under Arabic coffee plantation. The soil was sampled at a depth of 0 - 0.2 m, in a grid, totaling 50 sampling points. The following attributes were determined: pH in H2O, exchangeable acidity (Al3+), potential acidity (H + Al) and Al saturation (m). The data were analyzed by descriptive statistics and geostatistics. A fuzzy classification system was used based on the attributes described to infer on the soil acidity. Fuzzy logic, with the algorithm used, was efficient to characterize and map the uncertainties of acidity of the soil under study.

Index terms: geostatistics; kriging; hazy logic; soil management.

INTRODUÇÃO

Os solos podem ser naturalmente ácidos devido à pobreza em bases no material de origem ou a processos de formação que favorecem a remoção de elementos de reação básicos como K, Ca, Mg, Na, etc. A acidificação inicia-se ou acentua-se, devido à remoção de bases do solo e, também, à ação de cultivos e a adubações que, dependendo da forma de condução, podem levar a aumentos significativos dos níveis de acidez (Sousa & Lobato, 2004).

As culturas e mesmo as variedades ou cultivares variam muito na capacidade de tolerância à acidez e à disponibilidade de nutrientes. Dessa forma, as classes de fertilidade devem ser interpretadas, considerando as exigências específicas de cada cultivo agrícola (Martins, 2005).

Avanços na agropecuária têm mostrado a importância de novas avaliações das propriedades que alteram o rendimento das culturas, com o intuito de otimizar o aproveitamento de recursos e diminuir custos (Carvalho et al., 2002). Técnicas de análise que objetivem relacionar causas e efeitos e que considerem a natureza dos processos agrícolas, tendo em vista os modelos incertos associados à dinâmica desses processos, devem ser adotadas no gerenciamento da condução das culturas.

Nos últimos anos, observa-se o desenvolvimento de modelos quantitativos de predição espacial, geralmente derivados da geoestatística e de métodos como a lógica fuzzy (McBratney et al., 2000).

A lógica fuzzy objetiva modelar, de modo aproximado, o raciocínio humano, visando manipular informações em um ambiente de incerteza, fornecendo uma resposta aproximada para uma questão baseada em um conhecimento inexato, incompleto ou não totalmente confiável (Meirelles et al., 2007).

Essa técnica tem como apoio a teoria dos conjuntos fuzzy para sua representação e, utilizando expressões linguísticas na sua aplicação, definem-se conjuntos aos quais os valores são alocados com diferentes graus de pertinência entre falso e verdadeiro, ou seja, um elemento pode ou não pertencer a determinado conjunto fuzzy. Caso pertença, estabelece-se o grau de pertinência com que isso acontece (Bönisch et al., 2004).

A modelagem fuzzy tem sido amplamente utilizada por basear-se na caracterização de classes que não possuem, ou não podem definir limites rígidos entre si (Burrough & Mcdonnell, 1998), como acontece nas classes de fertilidade do solo, sendo indicada para lidar com ambiguidades, abstrações e ambivalências de modelos matemáticos complexos que representam limites difusos comuns em processos naturais (Bönisch et al., 2004).

O objetivo desse trabalho foi utilizar uma ponderação por meio da modelagem fuzzy para estudar, com base em atributos químicos do solo, a propagação de incertezas da acidez de um Latossolo Vermelho-Amarelo húmico cultivado com café arábica.

MATERIAL E MÉTODOS

O trabalho foi conduzido em uma área da fazenda Jaguaraí no município de Reduto, na região leste do Estado de Minas Gerais, a 20 ° 45 ' 45,4 '' de latitude S e 41 ° 32 ' 9,75 '' de longitude W. A área, com diferença de nível acentuada e 30 ° de declividade, vem sendo cultivada há cinco anos com Coffea arabica L., variedade Catuaí, no espaçamento de 2,0 x 0,60 m, em um Latossolo Vermelho-Amarelo húmico, com horizonte A bastante espesso e rico em matéria orgânica, conforme Embrapa (1999).

As amostragens de solo foram realizadas em uma malha irregular totalizando 50 pontos georrefe-renciados (Figura 1), com as amostras coletadas na projeção da copa de três plantas (constituindo o ponto de amostragem), na profundidade de 0-20 cm.


Os atributos avaliados foram: pH em H2O, acidez trocável (Al3+), acidez potencial (H + Al) e saturação por Al (m), e os valores encontrados foram analisados por meio das medidas de posição e dispersão da estatística descritiva e exploratória, bem como pelo teste Shapiro-Wilk's (p < 0,05) para testar a normalidade.

A análise da dependência espacial foi feita pela Geoestatística, com auxílio do software GS+ (Robertson, 1998), que realiza os cálculos das semivariâncias amostrais, cuja expressão pode ser encontrada em Vieira et al. (1983):

em que N(h) é o número de pares amostrais [z(xi); z(xi + h)] separados pelo vetor h, sendo z(xi) e z(xi + h) os valores numéricos observados do atributo analisado, para dois pontos xi e xi + h separados pelo vetor h. O semivariograma é representado pelo gráfico de g(h) função de h. Após o ajuste de um modelo matemático aos valores calculados de g(h), são definidos os parâmetros do modelo teórico para o semivariograma (efeito pepita, Co; alcance da dependência espacial, a; e patamar, C). Amostras separadas por distâncias menores que o alcance são espacialmente correlacionadas, ao passo que as separadas por distâncias maiores não o são.

Para análise do índice de dependência espacial (IDE), foram utilizados a relação C/(C0+C) e os intervalos propostos por Zimback (2001), que consideram as dependências espaciais fraca (IDE < 25 %); moderada (25 % < IDE < 75 %) e forte (IDE > 75 %).

Na determinação da dependência espacial, utilizou-se o exame de semivariogramas, por meio do programa GS+ (Robertson, 1998). Em caso de dúvida entre mais de um modelo para o mesmo semivariograma, utilizou-se o critério dos mínimos quadrados, optando-se, na seleção, pelos modelos com maior valor de R2 (coeficiente de determinação), menor SQR (soma de quadrado dos resíduos) e maior valor do coeficiente de correlação obtido pelo método de validação cruzada (Guimarães, 2000).

Visando testar o efeito de funções de álgebra de mapas na propagação de incertezas associadas às representações numéricas, os atributos numéricos pH, Al3+, H + Al e m foram integrados para espacialização da acidez. Nesta etapa, foi necessário, inicialmente, expressar cada atributo de acordo com as classes de acidez (Quadro 1).


As representações desses atributos foram submetidas a uma classificação contínua utilizando o mapeamento fuzzy. A função da associação escolhida para os valores de Al3+, H + Al e m foi a linear descrita por Bönisch et al. (2004). Para os valores de pH, utilizou-se a função exponencial apresentada por Meirelles et al. (2007).

- Conjunto de dados de Al3+, H + Al e m

MF A (Z) = 0 se z < q

MF A (Z) = (-1/ α)/(z-r) se q < z < r

MF A (Z) = 1 se z > r

- Conjunto de dados de pH

MF A (Z) = 1/(1 + d(z-p)2) se z < p

MF A (Z) = 1 se p < z < q

MF A (Z)= 1/(1+d(z-q)2) se z > q

em que MFA = função pertinência da variável Z;

α = r-q; d = q-p e p, q e r são valores dos limites de classe de dois conjuntos fuzzy pertencentes a um conjunto A.

A zona de transição fuzzy é dada pela inclinação das retas para o conjunto A com valores decrescentes (Al3+, H + Al e m) e com valores exponenciais (pH).

Os valores dos parâmetros p, q e r foram definidos em função dos valores dos limites dos intervalos dos atributos, que foram relacionados com a classe de acidez baixa (q e r) e com uma acidez elevada (p) (Quadro 2).


O efeito combinado dos atributos foi medido utilizando-se funções de álgebra de mapas por meio do operador fuzzy de soma ponderada. Foi atribuído peso igual para todas as representações, ou seja, considerou-se que os atributos contribuíram, de forma equitativa, na estimativa da acidez do solo. Dessa forma, a operação fuzzy se tornou uma média simples das representações.

RESULTADOS E DISCUSSÃO

Os resultados referentes à análise descritiva para o pH em H2O, Al3+ (cmolc dm-3), H + Al (cmolc dm-3) e m (%) do solo são apresentados no quadro 3. Tal análise, apesar de considerar os dados espacialmente independentes, tem a função de analisar seu comportamento geral e identificar influências que prejudicam a análise geoestatística.


Analisando o coeficiente de curtose, com exceção do H + Al, os atributos apresentaram distribuição platicúrtica, porém com achatamento reduzido da curva de distribuição normal, uma vez que esses valores encontram-se próximos de zero.

Considerando os coeficientes de variação (CV) em valores baixos (CV < 12 %); médios (12 % < CV < 60 %) e altos (CV > 60 %), segundo critérios de Warrick & Nielsen (1980), o pH apresenta baixa variação. Resultados semelhantes foram encontrados por Silva et al. (2003), Souza et al. (2004) e Zanão Júnior et al. (2007). O Al3+ e o H + Al apresentaram variação mediana, e a saturação por Al (m), variação elevada, justificado pela sua maior amplitude total, sugerindo influência direta dos extremos na sua variação. Souza (1992) afirma existir correlação diretamente proporcional entre a amplitude dos dados e o coeficiente de variação, o que foi constatado neste estudo.

Pelos valores médios dos atributos químicos do solo e com base na classificação apresentada por Alvarez V. et al. (1999), para o Estado de Minas Gerais, os níveis de pH e Al3+ são classificados como médios, enquanto os valores de m enquadram-se na classe de baixa saturação. A acidez potencial (H + Al) encontra-se classificada como alta.

Constata-se que todos os atributos apresentaram distribuição normal, analisada pelo teste Shapiro-Wilk's a 5 %, sendo esta distribuição simétrica, uma vez que os coeficientes de assimetria são próximos de zero, implicando valores de média e mediana coincidentes (Quadro 3). Segundo Isaaks & Srivastava (1989), mais importante que a normalidade dos dados é a ocorrência ou não do chamado efeito proporcional, ou seja, que a média e a variabilidade sejam constantes na área de estudo (Figura 4), em que se observa que os dados estão distribuídos homogeneamente em toda a área, não havendo regiões com concentração de valores altos ou baixos em zonas específicas da área estudada, ocorrendo assim a estacionaridade isotrópica - o que é essencial na aplicação da Geoestatística.


Devido à configuração da área (formato retangular com predomínio de pontos na direção norte-sul), assumiu-se a hipótese de isotropia dos dados para construção do semivariograma e posterior krigagem, uma vez que não é possível visualizar se há ou não anisotropia dos dados.

Analisando os semivariogramas (Figura 4), observa-se que todos os atributos apresentaram dependência espacial. Como os semivariogramas apresentaram patamares bem definidos, assumiu-se, nesse caso, estacionaridade intrínseca, uma vez que não se verificou tendência de variação dos atributos com as direções (Figura 3).


Ajustou-se aos valores de pH e H + Al o modelo esférico. Ao Al3+, o modelo que melhor se ajustou foi o gaussiano. Os alcances dos semivariogramas variaram de 12 m para o pH e de 56 m para o Al3+. O valor de alcance tem aplicação direta no plano de amostragem em agricultura de precisão, pois indica a distância máxima que uma variável está correlacionada espacialmente (Vieira, 2000), ou seja, determinações realizadas a distâncias maiores que o alcance têm distribuição aleatória e, por isso, são independentes entre si.

O índice de dependência espacial (IDE) mostra moderada dependência para Al3+, H + Al e m, e o pH apresenta forte dependência espacial, de acordo com classificação proposta por Zimback (2001).

Com exceção do pH, os valores do efeito pepita (C0) estão distantes de zero, o que reduz a exatidão de estimativas, por meio da krigagem, de valores em locais não medidos, uma vez que o C0 representa o componente da variabilidade espacial que não pode ser relacionado com uma causa específica (variabilidade ao acaso), ou seja, quanto menor a variação ao acaso, mais exata é a estimativa (Vieira, 2000).

A espacialização das possibilidades (incertezas) dos atributos, onde se observa que a quase totalidade da área apresenta valores de pertinência próximos de 1 para o pH, demonstrando baixa atividade de íons H+ em solução (Figura 5).


Graus de saturação de Al3+ e, ou, deficiência de Ca, frequentes em situações de pH baixos, restringem o crescimento do sistema radicular das plantas, tanto das camadas superficiais quanto subsuperficiais (Silva & Souza, 1998).

O Al3+ segue o comportamento dos atributos supra discutidos, com valores de pertinência próximos de 1, indicando sua precipitação na forma Al(OH)3, causando baixa concentração desse elemento na solução do solo, o que favorece o crescimento do sistema radicular das plantas, uma vez que o efeito primário da sua toxidez manifesta-se, geralmente, bloqueando o crescimento radicular (Epstein & Bloom, 2006).

A acidez potencial (H + Al), divergindo dos demais atributos, apresentou valores de pertinência próximos de zero. Como esta acidez é representada por H e Al na fase sólida do solo em equilíbrio com os íons da solução, sendo liberados quando há um aumento do pH, sugere-se redução da saturação por bases e elevado poder tamponante do solo, ou seja, maior resistência em ter o valor de seu pH alterado (Camargos, 2005). Essa característica é comum em solos ricos em matéria orgânica como os húmicos (Jacques, 2003).

A área em estudo apresenta graus de pertinência elevados para a acidez do solo, demonstrando que, de modo geral, a acidez na área é reduzida (Figura 6).


Os valores de pertinência indicam variação importante na área estudada e demonstram, de forma geral, boas possibilidades de crescimento da variedade Catuaí na área, uma vez que esta variedade é muito sensível à acidez, principalmente ao Al trocável (Rodrigues et al., 2006).

CONCLUSÕES

1. A lógica fuzzy, pelo algoritmo utilizado, teve um bom desempenho para caracterizar e mapear as incertezas da acidez do solo em questão.

2. A área em estudo apresentou graus de pertinência elevados para a acidez do solo, indicando acidez reduzida na área.

AGRADECIMENTOS

Á CAPES, pelo auxílio por meio da concessão de bolsa de mestrado.

LITERATURA CITADA

Recebido para publicação em setembro de 2008 e aprovado em abril de 2009.

  • ALVAREZ V., V.H.; NOVAIS, R.F.; BARROS, N.F.; CANTARUTTI, R.B. & LOPES, A.S. Interpretação dos resultados das análises de solos. In: RIBEIRO, A.C.; GUIMARÃES, P.T.G. & ALVAREZ V., V.H., eds. Recomendações para uso de corretivos e fertilizantes em Minas Gerais. 5Ş Aproximação. Viçosa, MG, Comissão de Fertilidade do Solo do Estado de Minas Gerais - CFSEMG, 1999. p.25-32.
  • BÖNISCH, S.; LOPES ASSAD, M.L.; CÂMARA, G. & MONTEIRO, A.M.V. Representação e propagação de incertezas em dados de solos, 1- Atributos categóricos. R. Bras. Ci. Solo, 28:21-32, 2004.
  • BURROUGH, P.A. & McDONNELL, R.A. Principles of geographic information systems. New York, Oxford University Press, 1998. 333p.
  • CAMARGOS, S.L. Acidez do solo e calagem (reação do solo). Material didático. Cuiabá, Universidade Federal do Mato Grosso, Faculdade de Agronomia e Medicina Veterinária, 2005. 26p.
  • CARVALHO, J.R.P.; SILVEIRA, P.M. & VIEIRA, S.R. Geoestatística na determinação da variabilidade espacial de características químicas do solo sob diferentes preparos. Pesq. Agropec. Bras., 37:1151-1159, 2002.
  • EMPRESA BRASILEIRA DE PESQUISA AGROPECUÁRIA - EMBRAPA. Centro Nacional de Pesquisa de Solos, Sistema brasileiro de classificação de solos. Rio de Janeiro, 1999. 412p.
  • EPSTEIN, E. & BLOOM, A.J. Nutrição mineral de plantas: Princípios e perspectivas. 2.ed. Londrina, Planta, 2006. 403p.
  • GUIMARÃES, E.C. Variabilidade espacial de atributos de uma Latossolo Vermelho-Escuro textura argilosa da região do cerrado, submetido ao plantio direto e ao plantio convencional. Campinas, Universidade Estadual de Campinas, 2000. 85p. (Tese de Doutorado)
  • ISAAKS, E.H. & SRIVASTAVA, R.M. An introdution to applied geostatistics. New York, Oxford University Press, 1989. 561p.
  • JACQUES, A.V.A. A queima das pastagens naturais - efeitos sobre o solo e a vegetação. Ci. Rural, 33:145-152, 2003.
  • MARTINS, C.E. Práticas agrícolas relacionadas à calagem do solo. Juiz de Fora, Embrapa - MG, 2005. 6p. (Comunicado Técnico)
  • McBRATNEY, A.B.; ODEH, I.O.A.; BISHOP, T.F.A.; DUNBAR, M.S. & SHATAR, T.M. An overview of pedometric techniques for use in soil survey. Geoderma, 97:293-327, 2000.
  • MEIRELLES, M.S.P.; MOREIRA, F.R. & CAMARA, G. Técnicas de inferência espacial. In: MEIRELLES, M.S.P.; CAMARA, G. & ALMEIDA, C.M., eds. Geomática: Modelos e aplicações ambientais. Brasilia, Embrapa, Informação Tecnológica, 2007. 593p.
  • ROBERTSON, G.P. GS+: Geostatistics for the environmental sciences - GS+ User´s guide. Plainwell, Gamma Desing Software, 1998. 152p.
  • RODRIGUES, L.A.; MARTINEZ, H.E.P.; NEVES, J.C.L.; NOVAIS, R.F. & MENDONÇA, S.M. Respostas nutricionais de cafeeiros catuaí e icatu a doses de calcário em subsuperfície. R. Bras. Ci. Solo, 30:985-995, 2006.
  • SILVA, C.R. & SOUZA, Z.M. Eficiência do uso de nutrientes em solos ácidos: Manejo de nutrientes e uso pelas plantas. Ilha Solteira, Universidade Estadual de São Paulo, 1998. (Caderno Didático)
  • SILVA, V.R.; REICHERT, J.M.; STORCK, L. & FEIJÓ, S. Variabilidade espacial das características químicas do solo e produtividade de milho em um Argissolo Vermelho-Amarelo distrófico arênico. R. Bras. Ci. Solo, 27:1013-1020, 2003.
  • SOUSA, D.M. & LOBATO, E. Cerrado: Correção do solo e adubação. 2.ed. Brasília, Embrapa Informação Tecnológica, 2004. 416p.
  • SOUZA, L.S. Variabilidade espacial do solo em sistemas de manejo. Porto Alegre, Universidade Federal do Rio Grande do Sul, 1992. 162p. (Tese de Doutorado)
  • SOUZA, Z.M.; MARQUES JÚNIOR, J.; PEREIRA, G.T. & MOREIRA, L.F. Variabilidade espacial do pH, Ca, Mg e V% do solo em diferentes formas do relevo sob cultivo de cana-de-açúcar. Ci. Rural, 34:32-39, 2004.
  • VIEIRA, S.R. Geoestatística em estudos de variabilidade espacial 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-53.
  • VIEIRA, S.R.; HATFIEL, J.L.; NIELSEN, D.R. & BIGGAR, J.W. Geoestatistical theory and application to variability of some agronomical properties. Hilgardia, 51:1-75, 1983.
  • WARRICK, A.W. & NIELSEN, D.R. Spatial variability of soil physical properties in the field. In: HILLEL, D., ed. Applications of soil physics. New York, Academic Press, 1980. p.319-44.
  • ZANÃO JÚNIOR, L.A.; LANA, R.M.Q. & GUIMARÃES, E.C. Variabilidade espacial do pH, teores de matéria orgânica e micronutrientes em profundidades de amostragem num Latossolo Vermelho sob semeadura direta. Ci. Rural, 37:PAGINAS, 2007.
  • ZIMBACK, C.R.L. Análise espacial de atributos químicos de solos para fins de mapeamento da fertilidade do solo. Botucatu, Universidade Estadual Paulista, 2001. 114p. (Tese de Livre-Docência)

Datas de Publicação

  • Publicação nesta coleção
    30 Out 2009
  • Data do Fascículo
    Ago 2009

Histórico

  • Aceito
    Abr 2009
  • Recebido
    Set 2008
Sociedade Brasileira de Ciência do Solo Secretaria Executiva , Caixa Postal 231, 36570-000 Viçosa MG Brasil, Tel.: (55 31) 3899 2471 - Viçosa - MG - Brazil
E-mail: sbcs@ufv.br