Accessibility / Report Error

Anisotropia no estudo da variabilidade espacial de algumas variáveis químicas do solo

Anisotropy to analyze spatial variability of some spatially referenced soil chemical variables

Resumos

No estudo do mapeamento da fertilidade do solo, pelas técnicas de geoestatística, algumas características estruturais da variabilidade espacial devem ser consideradas, tais como continuidade espacial e ausência de anisotropia. Neste contexto, o presente trabalho apresenta uma análise da anisotropia no estudo da variabilidade espacial das variáveis químicas do solo: ferro (Fe), acidez potencial (H + Al), matéria orgânica (MO) e Mn, de um conjunto de dados de 128 parcelas sem manejo químico localizado, espacialmente referenciados, estudados entre 1998 e 2002, em um Latossolo Vermelho distroférrico, em Cascavel-PR. A identificação da anisotropia foi realizada por meio da construção de semivariogramas direcionais com modelos ajustados, e a correção da anisotropia realizou-se por meio de transformações lineares e de modelos combinados. Em seguida, utilizou-se um modelo ajustado ao semivariograma omnidirecional para construção de mapas temáticos de variabilidade das variáveis estudadas. Observou-se a existência de anisotropia geométrica para a variável H + Al. Já as variáveis MO, Mn e Fe mostraram a presença de anisotropia combinada, sendo corrigida inicialmente a anisotropia geométrica e, posteriormente, a anisotropia zonal. Por meio do estudo da anisotropia, foi possível eliminar as direções privilegiadas, melhorando assim o ajuste dos semivariogramas e produzindo mapas temáticos das variáveis químicas estudadas com maior acurácia.

geoestatística; dependência espacial; agricultura de precisão


Studies on soil fertility mapping based on geostatistics should consider some spatial variability characteristics such as spatial continuity and absence of anisotropy. The present study reports an anisotropic analysis of the spatial variability of the following soil chemical variables: iron (Fe), potential acidity (H + Al), organic matter (OM) and manganese (Mn) of a data set of 128 plots without local chemical management, spatially referenced, studied between 1998 and 2002 in a dystrophic Red Latosol (Oxisol), in Cascavel, Paraná, Brazil. Anisotropy was identified through the construction of directional semivariograms with adjusted models and corrected through linear transformations and by combined models. Thereafter a model was used adjusted to the omnidirectional semivariogram for the construction of thematic variability maps of the studied variables. The existence of geometric anisotropy was observed for the variable H + Al. The variables MO, Mn and Fe evidenced the presence of combined anisotropy. The geometric anisotropy was first corrected and then the zonal anisotropy was adjusted. Based on anisotropy it was possible to eliminate the biased directions, enhancing the semivariogram fitting and precision of the thematic maps of chemical variables.

geostatistics; spatial dependence; precision agriculture


SEÇÃO I - FÍSICA DO SOLO

Anisotropia no estudo da variabilidade espacial de algumas variáveis químicas do solo1 1 Trabalho extraído da Dissertação de Mestrado de Engenharia Agrícola da Universidade Estadual do Oeste do Paraná - UNIOESTE.

Anisotropy to analyze spatial variability of some spatially referenced soil chemical variables

Luciana Pagliosa Carvalho GuedesI; Miguel Angel Uribe-OpazoII; Jerry Adriany JohannIII; Eduardo Godoy de SouzaIV

IMestre em Engenharia Agrícola pela Universidade Estadual do Oeste do Paraná e Doutora em Estatística e Experimentação Agronômica pela Escola Superior de Agricultura "Luiz de Queiroz" - ESALQ. Grupo de Pesquisa GGEA. E-mail: lpcarva@esalq.usp.br

IIProfessor Associado do Centro de Ciências Exatas e Tecnológicas - CCET, Universidade Estadual do Oeste do Paraná - UNIOESTE, Grupo de Pesquisa GGEA. Rua Universitária 2069, Sala 65, CEP 85819-110 Cascavel (PR). E-mail: mopazo@unioeste.br

IIIProfessor Assistente do CCET, Grupo de Pesquisa GGEA, da Universidade Estadual do Oeste do Paraná. E-mail: jerry.johann@hotmail.com

IVProfessor Associado do CCET, Grupo de Pesquisa GROSAP, da Universidade Estadual do Oeste do Paraná. E-mail: godoy@unioeste.br

RESUMO

No estudo do mapeamento da fertilidade do solo, pelas técnicas de geoestatística, algumas características estruturais da variabilidade espacial devem ser consideradas, tais como continuidade espacial e ausência de anisotropia. Neste contexto, o presente trabalho apresenta uma análise da anisotropia no estudo da variabilidade espacial das variáveis químicas do solo: ferro (Fe), acidez potencial (H + Al), matéria orgânica (MO) e Mn, de um conjunto de dados de 128 parcelas sem manejo químico localizado, espacialmente referenciados, estudados entre 1998 e 2002, em um Latossolo Vermelho distroférrico, em Cascavel-PR. A identificação da anisotropia foi realizada por meio da construção de semivariogramas direcionais com modelos ajustados, e a correção da anisotropia realizou-se por meio de transformações lineares e de modelos combinados. Em seguida, utilizou-se um modelo ajustado ao semivariograma omnidirecional para construção de mapas temáticos de variabilidade das variáveis estudadas. Observou-se a existência de anisotropia geométrica para a variável H + Al. Já as variáveis MO, Mn e Fe mostraram a presença de anisotropia combinada, sendo corrigida inicialmente a anisotropia geométrica e, posteriormente, a anisotropia zonal. Por meio do estudo da anisotropia, foi possível eliminar as direções privilegiadas, melhorando assim o ajuste dos semivariogramas e produzindo mapas temáticos das variáveis químicas estudadas com maior acurácia.

Termos de indexação: geoestatística, dependência espacial, agricultura de precisão.

SUMMARY

Studies on soil fertility mapping based on geostatistics should consider some spatial variability characteristics such as spatial continuity and absence of anisotropy. The present study reports an anisotropic analysis of the spatial variability of the following soil chemical variables: iron (Fe), potential acidity (H + Al), organic matter (OM) and manganese (Mn) of a data set of 128 plots without local chemical management, spatially referenced, studied between 1998 and 2002 in a dystrophic Red Latosol (Oxisol), in Cascavel, Paraná, Brazil. Anisotropy was identified through the construction of directional semivariograms with adjusted models and corrected through linear transformations and by combined models. Thereafter a model was used adjusted to the omnidirectional semivariogram for the construction of thematic variability maps of the studied variables. The existence of geometric anisotropy was observed for the variable H + Al. The variables MO, Mn and Fe evidenced the presence of combined anisotropy. The geometric anisotropy was first corrected and then the zonal anisotropy was adjusted. Based on anisotropy it was possible to eliminate the biased directions, enhancing the semivariogram fitting and precision of the thematic maps of chemical variables.

Index terms: geostatistics, spatial dependence, precision agriculture.

INTRODUÇÃO

Com o advento da mecanização e o aumento da produção agrícola mundial, o processo de produção agrícola passou a ser objeto de atenção e estudo de diversos pesquisadores, na busca de uma agricultura sustentável, capaz de maximizar a produção de forma racional, com diminuição do uso excessivo e irracional de insumos. Esse sistema de agricultura sustentável é chamado de Agricultura de Precisão.

Na Agricultura de Precisão a utilização dos insumos agrícolas é feita de forma diferenciada, com aplicações localizadas, possíveis apenas com um estudo da área agrícola por meio do mapeamento da variabilidade espacial das propriedades do solo, cujo objetivo é identificar e quantificar as necessidades locais quanto ao solo. A análise da variabilidade espacial é realizada por um conjunto de técnicas e conceitos, como a geoestatística, que considera como fator primordial a localização geográfica e a dependência espacial das variáveis em estudo, que aqui são denominadas variáveis regionalizadas. O semivariograma é o principal método geoestatístico para se detectar a variabilidade espacial. Trabalhos como os de Souza et al. (1999), Johann et al. (2002, 2004) e Mercante et al. (2003), entre outros, apresentam variáveis regionalizadas de propriedades do solo com estruturas de dependência espacial.

Segundo Ortiz (2002), o estudo e a análise estrutural do fenômeno por meio da geoestatística também devem considerar que a estrutura da variabilidade espacial está relacionada com a estrutura natural do fenômeno. Assim, algumas características devem ser verificadas, como a continuidade espacial e a ausência de anisotropia.

Uma variável regionalizada é isotrópica quando o padrão de continuidade espacial é idêntico em todas as direções. A isotropia é de suma importância para a estimativa de locais não-amostrados. É utilizada na construção de mapas de variabilidade, pois, por meio da identificação de que a variável analisada é isotrópica, o estudo da dependência espacial pode ser feito por meio de um único semivariograma, chamado de omnidirecional. O ajuste de um modelo para esse semivariograma pode ser feito por meio dos modelos transitivos apresentados na literatura (Isaaks & Srivastava, 1989), que consideram apenas variáveis isotrópicas. Além disso, a área de influência da variável regionalizada determina a forma de janela de estimativa a ser utilizada na krigagem (Li & Lake, 1994).

Quando o fenômeno em estudo revela diferentes padrões de dependência espacial, ou seja, apresenta uma variabilidade que não é a mesma em todas as direções, o fenômeno em estudo é chamado de anisotrópico. Segundo Clark (1979), a anisotropia é encontrada por meio de semivariogramas, que, construídos em diferentes direções, revelam uma privilegiada direção com um distinto padrão de dependência espacial.

Os principais tipos de anisotropia estudados na geoestatística são: geométrica, zonal e combinada. A anisotropia geométrica é aquela em que existe uma direção com maior continuidade espacial, isto é, maior valor de alcance (a) no semivariograma experimental em determinada direção; a anisotropia zonal, quando existe uma direção com maior valor de patamar (C) nos semivariogramas experimentais em relação às demais direções; e a anisotropia combinada, quando houver determinadas direções com diferentes valores de alcance e patamar nos semivariogramas experimentais.

Nesse contexto, o objetivo deste trabalho foi analisar a anisotropia no estudo da variabilidade espacial e sua influência na estimação de locais não-amostrados de algumas variáveis químicas do solo de uma área experimental com Latossolo Vermelho distroférrico.

MATERIAL E MÉTODOS

Os dados utilizados neste trabalho foram provenientes de uma área experimental localizada no Centro de Pesquisa Eloy Gomes da Cooperativa Central Agropecuária de Desenvolvimento Tecnológico e Econômico Ltda (COODETEC), em Cascavel (PR), com solo classificado como Latossolo Vermelho distroférrico (Embrapa, 1999). A área experimental total era de 1,74 ha, dividida em uma malha regular de 256 parcelas de 7,20 x 7,20 m, com corredor de 2,4 m em uma das direções, perfazendo um total de 128 parcelas com manejo químico localizado (CML) e 128 parcelas sem manejo químico localizado (SML). O delineamento utilizado nessa área é denominado de amostragem desalinhada sistemática estratificada (Wollenhaupt & Wolikowski, 1994; Souza et al., 1999).

Para coleta dos elementos amostrais do solo em cada parcela, a partir de um ponto de referência, foram tomadas quatro subamostras na camada de 0,0- 0,20 m de profundidade, dentro de um raio de 0,25 m, as quais, posteriormente, foram misturadas para compor uma amostra representativa de cada parcela. Essas amostras foram encaminhadas ao laboratório de análise de solo da COODETEC, para realização de análise química de rotina.

A área foi cultivada com soja e os dados são referentes às variáveis químicas do solo, a partir das amostras coletadas nos anos de 1998 e 2002. Especificamente, neste trabalho utilizou-se o conjunto de dados com as 128 parcelas sem manejo químico localizado. Diversas variáveis químicas do solo foram avaliadas e, entre elas selecionaram-se aquelas que apresentaram anisotropia, a saber: Fe e Mn, extraído com Mehlich-1 com amostras de 10 cm3 de solo para 100 cm3 de solução com decantação de 16 h; acidez potencial ou total (H + Al), obtida em solução-tampão de pH SMP de 7,5; e MO, determinada pelo método de Walkley Black.

Inicialmente, foi feito um estudo exploratório, com uma análise descritiva e gráfica dos dados distribuídos espacialmente. A análise gráfica foi composta inicialmente pela construção de gráficos espaciais de área classificados segundo os quartis (Post-plot), gráficos de dispersão das médias versus desvios-padrão de linhas e colunas, gráficos de médias móveis versus desvio-padrão e gráficos de dispersão dos valores da variável nas direções de 0 º, 45 º, 90 º e 135 º. Pelo desenvolvimento desses gráficos, foi possível obter informações sobre a existência ou não de tendência direcional, o efeito proporcional e a continuidade espacial das variáveis regionalizadas supramencionadas.

Dessa forma, os gráficos que apresentaram variação gradativa numa dada região ou direção indicaram possível existência de uma tendência direcional e, conseqüentemente, a falta de estacionariedade (Libardi et al., 1996).

Na análise geoestatística foi utilizado o semivariograma, que é um gráfico da função semivariância. O estimador da função semivariância proposto por Matheron (1963) para variáveis regionalizadas com distribuição normal de probabilidade é apresentado na equação 1:

em que n(h): número de pares em uma mesma distância h; e Z(si) e Z(si + h): valores da variável na posição si e si + h, com .

Para identificar a presença de anisotropia e avaliar a dependência espacial das variáveis regionalizadas em estudo, foram utilizados os semivariogramas nas direções convencionais de 0 º, 45 º, 90 º e 135 º, não sendo impedido trabalhar com outras direções (Isaaks & Srivastava, 1989).

A anisotropia geométrica foi corrigida por transformações lineares (Diggle & Ribeiro Junior, 2007), as quais são usadas na rotação e dilatação das coordenadas espaciais, utilizando a notação matricial apresentada na equação 2:

em que : matriz de rotação; e : matriz de dilatação.

Os parâmetros fator de anisotropia (Fa) e ângulo da direção de maior continuidade (a) foram calculados considerando . Para , a2 e a1 são, respectivamente, os alcances nas direções de maior e menor continuidade espacial.

Isaaks & Srivastava (1989) mostram que a anisotropia zonal pode ser modelada, por meio de um semivariograma direcional equivalente com a distância reduzida, considerando como patamar o maior valor de patamar apresentado entre os semivariogramas direcionais construídos. Adotou-se, dessa forma, o modelo de semivariograma corrigido descrito na equação 3:

em que , sendo w e a os valores de patamar e alcance do semivariograma na direção que apresentou a anisotropia zonal.

A anisotropia combinada consiste em uma combinação das anisotropias geométrica e zonal; assim, os modelos de semivariogramas direcionais apresentam diferentes alcances e patamares. Isaaks & Srivastava (1989) propuseram a equação 4 para o estudo da anisotropia combinada, em que a primeira etapa consiste em modelar a anisotropia geométrica e, a segunda, em modelar a anisotropia zonal.

em que:

Na equação anterior, w1 representa o patamar do semivariograma direcional que apresentou maior alcance, e w2, o que apresentou menor alcance, entre os semivariogramas direcionais. Já ax e ay representam os alcances nas direções x e y, respectivamente.

Corrigida a anisotropia, modelos geoestatísticos foram ajustados aos semivariogramas omnidirecionais e foram construídos mapas temáticos de variabilidade espacial para as variáveis regionalizadas na área em estudo, utilizando a técnica de interpolação por krigagem ordinária para pontos não-amostrados.

RESULTADOS E DISCUSSÃO

Análise descritiva geral

A estatística descritiva e os níveis descritivos (p-valores) obtidos pelo teste de normalidade de Shapiro-Wilk estão apresentados no quadro 1. Os dados obtidos indicam que as variáveis regionalizadas apresentaram baixa dispersão (DP), média homogeneidade na distribuição dos dados de Fe, MO e Mn (10 < CV < 20 %, Gomes, 2000) e alta heterogeneidade dos dados para H + Al (20 < CV < 30 %). Considerando um nível de 5 % de significância, todas as variáveis regionalizadas em estudo possuem distribuição normal de probabilidade.


Análise da tendência e estacionaridade

Na figura 1 são apresentados os gráficos espaciais da área em estudo, classificados segundo os quartis para as variáveis Fe, MO, Mn e H + Al. Nesses gráficos, pode-se observar que, no caso da variável Mn (Figura 1c), existe aumento gradativo no seu valor variável nas direções de 0 º, 45 º e 135 º, indicando assim a presença de tendência, principalmente na direção de 0 º. Nas variáveis MO (Figura 1b) e H + Al (Figura 1d), verifica-se variação gradativa nas direções de 45 º, 90 º e 135 º, o que indica a presença de tendência direcional destas variáveis, principalmente nas referidas direções. Não foi identificada tendência direcional na variável Fe (Figura 1a).



Como para atender à hipótese intrínseca não deve ser possível identificar padrões de tendência em qualquer direção, ajustou-se a essa tendência um modelo de regressão múltipla entre os valores da variável e as localizações geográficas. Por isso, a tendência foi eliminada e, assim, trabalhou-se com o resíduo resultante . A ausência de tendência nos resíduos pode ser identificada pelos gráficos espaciais de área, classificados segundo os quartis dos resíduos (Figura 2). Pode-se perceber nos resíduos uma distribuição aleatória dos pontos, evidenciando ausência de tendência direcional.


Pelos gráficos da relação entre médias e desviospadrão de cada variável em linhas e colunas e os gráficos de médias móveis versus desvio-padrão, não foi identificada relação entre a média e o desvio-padrão, verificando-se assim ausência de efeito proporcional. Trabalhos como os de Souza et al. (1999) e Mercante et al. (2003) efetuaram análises de variabilidade espacial cujos estudos preliminares foram feitos com variáveis físicas e químicas do solo.

Identificação da continuidade espacial

Para identificar, de forma exploratória, a existência da continuidade espacial de cada variável regionalizada, utilizou-se a técnica dos envelopes proposta por Diggle & Ribeiro Junior (2007), que consiste num recurso de simulação baseado em um teste gerado por simulação de Monte Carlo, de forma que, com base na hipótese de ausência de correlação espacial, podem-se permutar as posições dos valores amostrados, gerando um novo conjunto de dados, com o respectivo semivariograma. O envelope é construído pelos valores mínimos a máximos obtidos em todos os semivariogramas dos conjuntos de dados modificados, definindo assim uma região de independência espacial na variável em estudo. Na figura 3 é apresentado o semivariograma experimental omnidirecional para cada uma das variáveis em estudo, com os respectivos envelopes. Observa-se, em cada caso, que pelo menos um ponto do semivariograma apresenta uma semivariância fora dos limites de variações do envelope, indicando que os dados possuem alguma continuidade espacial e que a dependência espacial realmente existe nas variáveis em estudo.


Existência de anisotropia

Para identificar a existência da anisotropia nas referidas variáveis, foram construídos juntos semivariogramas direcionais (Figura 4), utilizando o estimador de Matheron definido na equação (1), variando a direção de 10 º em 10 º, com um ângulo de tolerância igual a 40 º, como sugerido por Isaaks & Srivastava (1989); esses valores foram escolhidos por meio da construção e análise de diversos semivariogramas, com diferentes ângulos e tolerâncias direcionais.



Dos semivariogramas direcionais, pode-se observar nas variáveis resíduos de Mn (Figura 4c) e MO (Figura 4b), maiores valores de alcance nas direções de 40 º e 50 º e valores maiores de patamar nas direções de 130 º e 140 º, sugerindo a presença de anisotropia combinada, em que foi considerada direção de menor e maior continuidade espacial de 45 º e 135 º, respectivamente. Na variável resíduo de H + Al (Figura 4d) podem-se perceber maiores valores de alcance na direção de 0 º, sugerindo maior continuidade espacial nesta direção, isto é, há presença de anisotropia geométrica. Para a variável Fe (Figura 4a), observam-se valores maiores de alcance e patamar nas direções de 0 º e 90 º, sugerindo a presença de anisotropia combinada.

A correção da anisotropia e o ajuste de um modelo teórico para o semivariograma omnidirecional corrigido trazem vantagens econômicas, uma vez que permitem a obtenção de melhor acurácia dos mapas de variabilidade espacial, podendo-se planejar melhor a malha amostral para posteriores estudos da variabilidade espacial na mesma área.

Os parâmetros dos modelos ajustados aos semivariogramas (Quadro 2) foram: efeito pepita (C0), que representa a descontinuidade; alcance (a), que descreve a distância máxima nas quais as amostras são consideradas com dependência espacial; e patamar (C = C0 + C1), que representa o valor da semivariância quando a distância é igual ao valor do alcance. Neste quadro foram incluídos o coeficiente de efeito pepita (ε = C0/(C0 + C1)), que mede o grau de dependência espacial da variável na área em estudo, o fator de anisotropia (Fα) e o ângulo de maior continuidade espacial (α).


No quadro 2, para a variável Fe, as direções de 90 º e 135 º apresentam diferentes estruturas no modelo do semivariograma, com diferentes valores de alcance e patamar, em que a maior continuidade espacial (maior alcance) situa-se na direção de 135 º e o maior valor do patamar ocorreu na direção de 90 º, identificando-se assim a presença de anisotropia combinada. Para os resíduos das variáveis MO e Mn observam-se diferentes valores de alcance e patamar nas direções de 45 º e 135 º; a maior continuidade espacial ocorreu na direção 135 º; e o maior valor de semivariância, na direção de 45 º, identificando-se também para estas variáveis a presença de anisotropia combinada. Quanto à variável H + Al, houve diferença nos alcances das direções de 0 º e 90 º, em que a direção de maior continuidade espacial foi em 90 º, identificando-se, portanto, a anisotropia geométrica. Essas diferenças entre os semivariogramas que apresentaram anisotropia são também observadas nos gráficos (Figura 5), que se referem aos modelos ajustados e sobrepostos aos semivariogramas experimentais direcionais das variáveis Fe e resíduos de MO, Mn e H + Al.


Correção da anisotropia

Identificada a presença de anisotropia, esta foi corrigida mediante transformações lineares nas coordenadas espaciais, com o modelo proposto por Diggle & Ribeiro Junior (2007), para variáveis com anisotropia geométrica, e com o modelo aditivo de Isaaks & Srivastava (1989), para variáveis com anisotropia combinada. A correção, numa primeira etapa, é da anisotropia geométrica; em seguida, efetuase a correção da anisotropia combinada.

Assim, corrigindo as coordenadas espaciais das variáveis, construíram-se novos semivariogramas nas direções de maior e menor continuidade espacial (Quadro 3 e Figura 6).



Observa-se similaridade nos semivariogramas ajustados, com valores de alcance e patamar muito próximos entre as referidas direções (Quadro 3), evidenciando assim a transformação de variável anisotrópica geométrica em isotrópica. Para as variáveis Fe e resíduos de MO e Mn, foi corrigida em uma primeira etapa a anisotropia geométrica e construídos semivariogramas nas direções de maior e menor continuidade espacial, constatando-se dessa forma uma similaridade em suas estruturas de variabilidade espacial, com valores similares de alcance (a) e também de patamar (C); assim, tanto a anisotropia geométrica quanto a zonal podem ser consideradas corrigidas. Esse fato provavelmente ocorreu, pois, além de a anisotropia zonal ser um caso particular da anisotropia geométrica, a diferença entre os patamares nos modelos direcionais não foi tão significativa; dessa forma, eliminando a anisotropia geométrica, eliminouse também parte da anisotropia zonal.

Corrigida a anisotropia, foi possível utilizar um único semivariograma que representa todas as direções, denominado semivariograma omnidirecional (Quadro 4). O modelo foi ajustado para cada semivariograma experimental das variáveis em estudo, utilizando-se o método dos mínimos quadrados ponderados (WLS1) (Isaaks & Srivastava, 1989). Por esses parâmetros, todas as variáveis apresentam uma moderada dependência espacial, pois seus valores de coeficiente de efeito pepita (ε) estão entre 25 e 75 % (Souza et al., 1999). Assim, chega-se à conclusão de que, para as variáveis Fe e resíduos de Mn, MO e H + Al, apenas 40,0, 22,23, 35,43 e 40,62 % da variação total dos dados, respectivamente, são explicados pela dependência espacial.


A partir dos parâmetros obtidos para o semivariograma omnidirecional, foi possível construir mapas que expressam a variabilidade das variáveis avaliadas na área em estudo, representando uma superfície contínua que caracteriza o comportamento das variáveis no campo (Figura 7). Observou-se que, para as variáveis MO (Figura 7b) e Mn (Figura 7c), existe maior continuidade das sub-regiões na direção de 135 º. Para Fe (Figura 7a) há maior continuidade na direção de 0 º e, para H + Al (Figura 7d), uma aparente continuidade espacial das sub-regiões na direção de 90 º.


 




Foram também construídos mapas de variabilidade segundo a classificação para interpretação de resultados de análise do solo (Quadro 5), de acordo com EMATER-PARANÁ (Oleynik et al., 1998).


Pela classificação adotada (Oleynik et al., 1998), as variáveis Fe (Figura 8a) e Mn (Figura 8c) apresentaram-se homogêneas em toda a área de estudo, com nível alto. A variável MO (Figura 8b) também se apresentou homogênea por toda a área em estudo, com um nível médio. Já a variável H + Al (Figura 8d) foi a única a apresentar duas classificações diferentes no mapa: nível médio na parte inferior e nível alto no restante do mapa.


 




CONCLUSÕES

1. Sugere-se aos pesquisadores a importância das verificações da existência da anisotropia e a utilização de sua correção, de forma que ela não influencie a construção de mapas temáticos.

2. No estudo da anisotropia, foi possível eliminar as direções privilegiadas que continham maiores valores do alcance, implicando a existência de áreas com maior continuidade espacial nessa direção, melhorando assim o ajuste dos semivariogramas e produzindo, conseqüentemente, uma predição espacial mais eficiente das variáveis em estudo em localizações não-amostradas. Isso possibilitou que os mapas temáticos que foram construídos proporcionassem uma noção visual da variabilidade das variáveis Fe, MO, Mn e H + Al.

3. Além de uma melhor acurácia dos mapas de variabilidade espacial gerados, a análise da anisotropia em variáveis com dependência espacial possibilita melhor planejamento das malhas amostrais para futuros estudos na mesma área experimental.

AGRADECIMENTOS

Ao CNPq e à Fundação Araucária, pelo apoio financeiro, e à COODETEC, pelo apoio técnico.

LITERATURA CITADA

Recebido para publicação em agosto de 2004 e aprovado em outubro de 2008.

  • CLARK, I. Practical geostatistics. London, Applied Science Publishers, 1979. 126p.
  • DIGGLE, P.J. & RIBEIRO JUNIOR, P.J. Model-based geostatistics. New York, Springer, 2007. 230p.
  • EMPRESA BRASILEIRA DE PESQUISA AGROPECUÁRIA - EMBRAPA. Centro Nacional de Pesquisa de Solos. Sistema brasileiro de classificação de solos. Brasília, 1999. 412p.
  • GOMES, F.P. Curso de estatística experimental. 14.ed. Piracicaba, Nobel, 2000. 477p.
  • ISAAKS, E.H. & SRIVASTAVA, R.M. An introduction to applied geostatistics. New York, Oxford University Press, 1989. 561p.
  • JOHANN, J.A.; URIBE-OPAZO, M.A.; SOUZA, E.G. & ROCHA, J.V. Modelagem da correlação espacial entre mapas de produtividade da soja e mapas de atributos do solo por meio de análise de componentes principais. R. Eng. Agríc., 22:383-390, 2002.
  • JOHANN, J.A.; URIBE-OPAZO, M.A.; SOUZA, E.G. & ROCHA, J.V. Variabilidade espacial dos atributos físicos do solo e da produtividade em uma área experimental de agricultura de precisão com e sem manejo localizado. R. Eng. Agric. Amb., 8:212-219, 2004.
  • LI, D. & LAKE, L.W. A moving window semivariance estimator. Water Res. Res., 30:1479-1489, 1994.
  • LIBARDI, P.L.; MANFRON, P.A.; MORAES, S.O. & TOUN, R.L. Variabilidade da umidade gravimétrica de um solo hidromórfico. R. Bras. Ci. Solo, 20:1-12, 1996.
  • MATHERON, G. Principles of geoestatistics. Econ. Geol., 58:1246-66, 1963.
  • MERCANTE, E.; URIBE-OPAZO, M.A. & SOUZA, E.G. Variabilidade espacial e temporal da resistência mecânica do solo à penetração em áreas com e sem manejo químico localizado. R. Bras. Ci. Solo, 27:1149-1153, 2003.
  • OLEYNIK, J.; BRAGAGNOLO, N.; BUBLITZ, U. & SILVA, J.C.C. Análise do Solo: Tabela para transformação de resultados analíticos e interpretação de resultados. 5.ed.Curitiba, Emater-Paraná, 1998. 64p.
  • ORTIZ, G.C. Aplicação de métodos geoestatísticos para identificar a magnitude e a estrutura da variabilidade espacial de variáveis físicas do solo. Piracicaba, Escola Superior de Agricultura Luiz de Queiroz, 2002. 91p. (Tese de Mestrado)
  • SOUZA, E.G.; JOHANN, J.A.; ROCHA, J.V.; RIBEIRO, S.R.A.; SILVA, M.S.; URIBE-OPAZO, M.A.; MOLIN, J.P.; OLIVEIRA, E.F. & NÓBREGA, L.H.P. Variabilidade espacial dos atributos do solo em um Latossolo Roxo distrófico na região de Cascavel-PR. R. Eng. Agríc., 8:8092, 1999.
  • WOLLENHAUPT, N.C. & WOLICOWSKI, R.P. Grid soil sampling. Better Crops Plant Food, 78:6-9, 1994.
  • 1
    Trabalho extraído da Dissertação de Mestrado de Engenharia Agrícola da Universidade Estadual do Oeste do Paraná - UNIOESTE.
  • Datas de Publicação

    • Publicação nesta coleção
      16 Fev 2009
    • Data do Fascículo
      Dez 2008

    Histórico

    • Aceito
      Out 2008
    • Recebido
      Ago 2004
    Sociedade Brasileira de Ciência do Solo Sociedade Brasileira de Ciência do Solo, Departamento de Solos - Edifício Silvio Brandão, s/n, Caixa Postal 231 - Campus da UFV, CEP 36570-900 - Viçosa-MG, Tel.: (31) 3612-4542 - Viçosa - MG - Brazil
    E-mail: sbcs@sbcs.org.br