Amostragem e variabilidade espacial de atributos físicos de um latossolo vermelho em Mossoró, RN

Os solos, mesmo aqueles considerados homogêneos, exibem alguma variabilidade espacial e temporal em seus atributos físicos, químicos e biológicos. A definição dessa variabilidade é muito importante para um eficiente manejo do solo. Este trabalho objetivou estudar a variabilidade espacial de alguns atributos físicos de um Latossolo Vermelho, em Mossoró-RN, a fim de estabelecer o melhor procedimento na coleta de amostras de solo. Em uma área com dimensões de 100 x 100 m, marcada a cada 10 m em duas direções perpendiculares entre si, fez-se a coleta de amostras de solo na camada de 0,0-0,2 m de profundidade, totalizando 100 amostras. As análises compreenderam granulometria, densidade do solo, densidade de partículas, porosidade total e retenção de água nas tensões de 0,033 e 1,500 MPa. Os dados foram analisados por estatística descritiva e geoestatística. Concluiu-se que a variabilidade do solo, medida pelo coeficiente de variação, mostrou-se baixa para areia, densidade do solo, densidade de partículas e volume total de poros e média para silte, argila, capacidade de campo e ponto de murcha permanente; a maioria das variáveis analisadas segue distribuição normal de freqüências, excetuando-se areia, densidade do solo e volume total de poros; o mapeamento das variáveis estudadas sugere o estabelecimento de subáreas homogêneas, a fim de tornar o manejo dos recursos mais eficiente, especialmente o de água.

Métodos de ajuste e modelos de semivariograma aplicados ao estudo da variabilidade espacial de atributos físico-hídricos do solo

Avaliou-se a variabilidade e a estrutura de dependência espacial de atributos físico-hídricos do solo (volume total de poros - VTP, condutividade hidráulica saturada - k0, porosidade drenável - PD e conteúdo de água volumétrico na capacidade de campo - tetacc), na bacia hidrográfica do Ribeirão Marcela, representativa do domínio dos Latossolos na Região Alto Rio Grande, ajustando-se semivariogramas do tipo esférico e exponencial pelos métodos dos Quadrados Mínimos Ponderados (QMP) e Máxima Verossimilhança (MV). As amostras para caracterização físico-hídrica foram coletadas na camada de 0 a 0,15 m de profundidade, obedecendo a grids de 240 x 240 m e 60 x 60 m, totalizando 165 pontos amostrais. Os atributos VTP e tetacc apresentaram baixa variabilidade, com coeficiente de variação (CV) de 7,54 e 10,22 %, respectivamente, sendo a variabilidade do atributo PD média, com CV de 43,2 %, e alta para k0, com CV de 88,37 %. Todos os atributos apresentaram estrutura de dependência espacial, com forte grau para os atributos VTP, k0 e tetacc, com GD de 81,82, 75,31 e 82,61 %, respectivamente, e moderada para PD, com GD na ordem de 54,88 %. Os alcances para VTP, tetacc e PD foram da ordem de 1.000,0 m, enquanto o atributo k0 mostrou alcance de 77,46 m. Verificou-se desempenho similar entre os métodos empregados, sendo ambos indicados para ajuste de semivariogramas a atributos físico-hídricos do solo. Entre os modelos de semivariograma, baseando-se na autovalidação, o esférico seria indicado para os atributos VTP, k0 e tetacc, apresentando quadrado médio do erro na ordem de 0,00161; 0,2795 e 0,00155, respectivamente; e o exponencial, para o atributo PD, com quadrado médio do erro de 0,004070.

Variabilidade espacial de atributos químicos do solo em áreas sob sistema plantio convencional

A variabilidade espacial de atributos químicos de solo foi avaliada em uma lavoura comercial cultivada sob sistema plantio convencional, em Uberlândia (MG), no ano de 2004. O objetivo foi avaliar a distribuição e a dependência espacial dos atributos químicos do solo em uma lavoura sob sistema plantio convencional. A grade de amostragem foi estabelecida na Fazenda Santa Rosa, sendo o solo classificado como Latossolo Vermelho textura muito argilosa (680 g kg-1 de argila). Coletaram-se dados do solo, dispostos segundo uma malha de 121 pontos amostrais, espaçados de 50 m, analisados por meio da geoestatística, na profundidade de 0-0,2 m. Foram determinados o pH em água; P e K disponíveis; Ca, Mg e Al trocáveis, H + Al, B, Cu, Fe, Mn, Zn, S e MO. Calcularam-se a saturação por bases (V %), CTC total (T) e soma de bases (SB). Os dados foram avaliados por estatística descritiva e pela análise de dependência espacial, com base no ajuste de semivariogramas. A maioria dos atributos apresentou coeficiente de variação (CV) alto, sendo o maior encontrado para P, 73,51 %, e o menor, para pH em água, 5,96 %. A grande maioria dos atributos avaliados mostrou dependência espacial, a qual foi classificada como moderada e forte. A maioria dos dados se ajustou ao semivariograma de modelo gaussiano.

Variabilidade espacial de propriedades físicas e químicas de um Cambissolo sob dois sistemas de manejo de cana-de-açúcar

A prática da queima da palha da cana-de-açúcar antes da colheita tem sido objeto de intensas discussões tanto quanto a possíveis alterações de propriedades dos solos onde é cultivada como com respeito à emissão de gases nocivos à atmosfera. O objetivo deste estudo foi avaliar variações em propriedades físicas e químicas de duas áreas de Cambissolos cultivadas com cana-de-açúcar: uma submetida a cultivo contínuo com aração e gradagem e outra onde a lavoura se encontrava na qüinquagésima soca sem que nenhuma forma mecânica de mobilização do solo tivesse sido aplicada. As áreas em estudo estão localizadas na região norte do Estado do Rio de Janeiro, no município de Campos dos Goytacazes. Foram coletadas amostras de solo em uma área total de 4.422 m², num reticulado de sete linhas e sete colunas, com espaçamento constante de 9,5 m entre linhas e colunas. Dessa forma, formaram-se 49 células, em cujo centro foram localizados os pontos de amostragem. Foram utilizados semivariogramas para avaliar a estrutura da variância e caracterizar a dependência espacial das propriedades do solo em cada área. Verificaram-se diferenças quanto a granulometria do material de origem dos dois solos, uma vez que os semivariogramas apontaram correlações distintas para cada área. Ainda, o sistema de manejo tendeu a afetar as características relacionadas à porosidade do solo, principalmente na faixa de 0 a 0,20 m de profundidade. Quanto aos elementos Ca, Mg e K, os resultados da área com cultivo contínuo apresentaram menor correlação espacial. Já o tratamento sem cultivo, mesmo não recebendo adubação, apresentou níveis nutricionais próximos ou superiores aos do tratamento cultivo contínuo. Conclui-se que o comportamento espacial das frações granulométricas, diferentes nas duas áreas, aliado à diferença de médias desses atributos, indicou que os solos parecem ter material de origem de granulometrias diferentes; que o cultivo contínuo aumentou a continuidade da variabilidade espacial de atributos relacionados à porosidade e alterou também a correlação espacial para os elementos Ca, Mg e K.

Variabilidade espacial de atributos químicos do solo sob diferentes usos e manejos

O uso e manejo do solo e da cultura são importantes condicionadores da variabilidade de atributos do solo. O trabalho foi desenvolvido em Selvíria (MS), com o objetivo de avaliar a variabilidade espacial do pH, potássio (K), cálcio (Ca), magnésio (Mg) e saturação por bases (V) em Latossolo Vermelho sob diferentes usos e manejos. Os solos foram amostrados em malha, com intervalos regulares de 2 m, perfazendo o total de 64 pontos, nas profundidades de 0,0-0,1 e 0,1-0,2 m, nas seguintes áreas: vegetação natural (Cerrado), plantio direto, plantio convencional e pastagem. As maiores variabilidades, medidas por meio do coeficiente de variação, foram observadas para K, Mg e Ca; o pH apresentou o menor coeficiente de variação nos diferentes usos e manejo do solo, e o atributo V, coeficiente de variação médio. Os sistemas preparo convencional e pastagem apresentaram os menores alcances quando comparado aos sistemas Cerrado e plantio direto.

Seleção de modelos de variabilidade espacial para elaboração de mapas temáticos de atributos físicos do solo e produtividade da soja

Pesquisas sobre a variabilidade espacial dos atributos do solo que influenciam a produtividade são de uma grande importância para o desenvolvimento de novas técnicas que beneficiam a agricultura. A variabilidade desses atributos pode ser avaliada por técnicas de geoestatística e auxiliar no mapeamento e manejo do solo. Este trabalho teve por objetivo avaliar a qualidade do ajuste dos modelos teóricos espaciais segundo o Critério de Informação de Akaike, de Filliben, de Validação Cruzada e o valor máximo do logaritmo da função verossimilhança, de dados da umidade do solo, da densidade do solo e da resistência do solo à penetração, nas camadas de 0 a 0,1, 0,1 a 0,2 e 0,2 a 0,3 m, e de produtividade da soja do ano agrícola 2004-2005. Os parâmetros dos modelos de variabilidade espacial foram estimados por meio dos métodos de mínimos quadrados ordinários, mínimos quadrados ponderados e máxima verossimilhança. A pesquisa foi desenvolvida em uma área de 57 ha de um Latossolo Vermelho distroférrico, utilizando-se uma malha de 75 x 75 m georreferenciada. Concluiu-se que, dos métodos de avaliação de ajustes estudados, o da Validação Cruzada foi o mais adequado para escolha do melhor ajuste do modelo de variabilidade espacial; conseqüentemente têm-se mapas temáticos mais acurados.

Variabilidade espacial dos parâmetros hidrodinâmicos de duas parcelas agrícolas no estado da Paraíba

O conhecimento das curvas de retenção e da condutividade hidráulica da água do solo é indispensável para a modelagem do transporte de água e de poluentes em sistemas agrícolas. Essas propriedades são afetadas por numerosas fontes de variabilidade, geralmente associadas a fatores espacial e temporal, e relacionadas com o manejo agrícola. O objetivo desse trabalho foi analisar a variabilidade espacial dos parâmetros das curvas de retenção da água (van Genuchten) e da condutividade hidráulica (Brooks & Corey) da água do solo, obtidos com o método "Beerkan", que se baseia em ensaios de infiltração simplificados e nas informações texturais do solo. Esses parâmetros são de forma ou texturais (parâmetros de ajuste que dão forma às curvas) e de normalização ou estruturais (umidade volumétrica saturada, condutividade hidráulica saturada do solo e potencial de entrada de ar). Os ensaios de campo foram conduzidos em uma malha regular de 25 x 25 m num Neossolo Flúvico e num Latossolo Amarelo cultivados com feijão (Vigna unguiculata (L.) Walp.). Para ambos os solos, os valores dos parâmetros de forma, da umidade volumétrica saturada, dos logaritmos da condutividade hidráulica saturada e do potencial de entrada de ar apresentaram dependência espacial. Os valores da umidade volumétrica saturada e dos logaritmos do potencial de entrada de ar foram mais sensíveis à variação a pouca distância no Neossolo Flúvico do que do Latossolo Amarelo, com semivariogramas apresentando estruturas de variação com fraca dependência espacial.

Simulação da variabilidade espacial da erosão hídrica em uma sub-bacia hidrográfica de Latossolos no sul de Minas Gerais

A simulação da distribuição espacial da erosão do solo consiste em uma ferramenta poderosa para o planejamento conservacionista em bacias hidrográficas, sendo uma importante aplicação da Equação Universal de Perdas de Solo (EUPS) associada a princípios de interpolação espacial, principalmente a geoestatística. Este trabalho objetivou simular a distribuição espacial da erosão hídrica numa sub-bacia hidrográfica da região do Alto Rio Grande (MG), aplicando ferramentas geoestatísticas para distribuição espacial e mapeamento. Diferentes cenários de uso do solo foram analisados. A erosão foi estimada com base na EUPS aplicada a células, considerando as unidades pedológicas e diferentes usos em cada célula. O fator topográfico (LS) das células foi obtido com base no Modelo Digital de Elevação da sub-bacia, identificando comprimento e direção principal do escoamento. A erosividade média anual da região é de 8.030 MJ mm ha-1 h-1 ano-1, e a erodibilidade dos solos foi extraída da literatura. A EUPS foi aplicada a cada célula, levando-se em conta a situação atual do solo, áreas degradadas plantadas com eucalipto, pastagem plantada e plantio convencional de milho, considerando a ocupação de toda a sub-bacia hidrográfica. Na situação atual, a sub-bacia apresenta taxas de erosão inferiores aos limites de tolerância para os respectivos solos, com exceção das áreas degradadas ocupadas por eucalipto e pastagem em Cambissolo. No entanto, em todas as situações analisadas, seu lado leste apresentou as maiores perdas de solo, especialmente para os cenários de eucalipto nas condições atuais e plantio convencional de milho em Cambissolo e Latossolo Vermelho-Amarelo, sendo necessária a aplicação de técnicas conservacionistas.

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

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.

Variabilidade espacial de atributos mineralógicos de um latossolo sob diferentes formas do relevo: I - mineralogia da fração argila

A distribuição dos minerais da fração argila depende de condições pedoambientais específicas. O estudo das relações entre atributos desses minerais e locais específicos de sua ocorrência na paisagem torna-se importante para o entendimento das relações entre a mineralogia e os demais atributos do solo. Este trabalho teve como objetivo avaliar os atributos cristalográficos, as razões Gt/(Gt+Hm) e Ct/(Ct+Gb) dos minerais da fração argila (hematita, goethita, caulinita e gibbsita) e o padrão de distribuição espacial destes em um Latossolo Vermelho eutroférrico sob cultivo de cana-de-açúcar há mais de 30 anos, com colheita no sistema cana crua. As amostras de solo foram coletadas nos pontos de cruzamento de uma malha, georreferenciados, com intervalos regulares de 10 m, e dimensão de 100 x 100 m, na profundidade de 0,2-0,4 m, totalizando 119 pontos, em uma área de 1 ha. As amostras de solos foram analisadas quimicamente. Posteriormente, a fração argila foi submetida à difração de raios X. Foram feitas as análises estatísticas e geoestatísticas dos atributos cristalográficos dos minerais analisados. Mapas de krigagem foram realizados para a visualização da variabilidade espacial dos atributos. Todos os atributos estudados apresentaram dependência espacial, e os atributos mineralógicos apresentaram-se relacionados com variações das formas do relevo. A gibbsita apresentou o maior diâmetro médio do cristal (DMC) entre os minerais da fração argila estudados, e a goethita foi o mineral que apresentou a maior variação do DMC e da largura à meia altura. Os maiores valores do DMC da goethita, hematita e gibbsita foram encontrados no compartimento II, caracterizado por uma área côncavo-convexa.

Variabilidade espacial de atributos mineralógicos de um latossolo sob diferentes formas do relevo: II - correlação espacial entre mineralogia e agregados

A caulinita (Ct), a gibbsita (Gb), a goethita (Gt) e a hematita (Hm) são os principais minerais da fração argila dos Latossolos brasileiros. A cristalinidade desses minerais interfere nos atributos físicos do solo, que são dependentes do relevo. Este trabalho teve como objetivo correlacionar espacialmente a estabilidade e o diâmetro médio geométrico dos agregados com os atributos cristalográficos dos minerais da fração argila de um Latossolo Vermelho eutroférrico sob cultivo de cana-de-açúcar com colheita no sistema cana crua. As amostras de solo foram coletadas nos pontos de cruzamento de uma malha, georreferenciados, com intervalos regulares de 10 m e dimensão de 100 x 100 m, na profundidade de 0,2- 0,4 m, totalizando 119 pontos, em uma área de 1 ha, onde as variações da forma do relevo foram mapeadas. Foram determinados o diâmetro médio geométrico e a percentagem de agregados retidos nas classes > 2 e < 1 mm das amostras, e a mineralogia da fração argila foi caracterizada por difratometria de raios X. Esses resultados foram submetidos à análise estatística e geoestatística. Buscando correlacionar a estabilidade e as diferentes classes de tamanho dos agregados com a cristalinidade dos minerais da fração argila do solo, foram realizados semivariogramas cruzados. Observou-se que goethitas de menor grau de cristalinidade se relacionam espacialmente com agregados de maior diâmetro médio geométrico e tamanho dos agregados e que os atributos estudados e suas relações são dependentes das formas do relevo.

Variabilidade espacial de atributos de solos em unidades de manejo em área piloto de produção integrada de coco

O objetivo deste trabalho foi, inicialmente, acessar as variações espaciais dos solos de parcelas experimentais de produção integrada do coco por meio do mapeamento detalhado dos solos e da construção de mapas de contorno de atributos do solo. Foram feitas observações, medições e coleta de amostras em microtrincheiras e perfis de solo. O teor de argila dos solos das microtrincheiras foi avaliado a campo pela sensação ao tato. Em todas as amostras foram realizadas análises físicas e químicas. Foram determinados granulometria (perfis), Al, Ca e Mg, acidez potencial, K e Na, pH, P disponível e C orgânico. Nas amostras dos perfis, foram determinados Fe, Al, Ti, P e Si pelo ataque sulfúrico. O Fe também foi determinado na fração argila, por extração com ditionito-citrato-bicarbonato de Na. Foi realizada difração de raios X da fração argila desferrificada. A construção de semivariogramas possibilitou a classificação e a comparação dos atributos do solo por meio do grau de dependência espacial destes. A elaboração dos mapas de contorno foi realizada pelo método interpolador de krigagem. A área mapeada corresponde a um topo amplo de tabuleiros costeiros, com uma depressão fechada no seu extremo sul. Do centro da depressão para fora dela seguem, em seqüência, Cambissolo Háplico textura argilosa, Latossolo Amarelo textura média/argilosa e Latossolo Amarelo textura média. Esses Latossolos são coesos em subsuperfície. Predominam atributos com elevado grau de heterogeneidade. Os valores médios de P, K e Ca + Mg estão influenciados por correções e adubações sistemáticas realizadas na área de plantio. Todas as variáveis analisadas apresentaram dependência espacial, expressa pelos modelos de semivariogramas. A análise da relação C0/(C0 + C1) revelou um grau de dependência espacial, de todas as variáveis, de moderado a forte. Os mapas de contorno detalharam diferenças que o mapa de solos já apontava para vários atributos, influenciados pela presença da depressão no extremo sul da área. Foram estabelecidas duas unidades de manejo para a área de estudo, as quais exigem práticas de manejo de água e solo diferenciadas.

Variabilidade espacial de atributos químicos de um Latossolo Vermelho-Amarelo húmico cultivado com café

O objetivo do presente estudo foi analisar a variabilidade espacial de atributos químicos de um Latossolo Vermelho-Amarelo húmico. O estudo foi realizado em uma área cultivada com Coffea arabica L. variedade Catuaí. As amostras de solo foram coletadas na profundidade de 0-20 cm, distribuídas em uma malha amostral, com 50 pontos. Os atributos químicos estudados foram: P, Na e S disponíveis, Ca, Mg e Al trocáveis, pH, H + Al, SB, t, T, V, m, MO, ISNa, P-remanescente e micronutrientes (Zn, Fe, Mn, Cu e B). Os dados foram analisados por meio da análise da estatística descritiva, seguida de uma análise de componentes principais com posterior análise de agrupamento. A geoestatística foi utilizada para verificar a existência e quantificar o grau de dependência espacial dos componentes principais. O método de análise multivariada com base nos componentes principais forneceu componentes interpretáveis, sendo o primeiro relacionado com a acidez e o segundo com o teor de matéria orgânica do solo. As técnicas de análise multivariada, em associação com a geoestatística, facilitaram a avaliação da variabilidade do solo; os componentes principais 1 e 2 apresentaram dependência espacial moderada, com maior continuidade espacial observada para o componente 1, o que permitiu melhor caracterização da acidez do solo. O solo em estudo apresentou, ao longo da encosta, disponibilidade reduzida de nutrientes, considerável acidez, menores teores de matéria orgânica na porção central da área e valores elevados de Al trocável na porção superior desta.

Variabilidade espacial dos teores de macronutrientes em latossolos sob sistema plantio direto

A eficácia dos planos de amostragem do solo pode ser aumentada se for conhecida e considerada a variabilidade espacial de seus atributos, e, para isso, devem ser estudados os fatores que a determinam. O objetivo deste trabalho foi avaliar a variabilidade espacial dos teores de macronutrientes P, K, Ca, Mg e S em dois Latossolos, um muito argiloso e outro de textura média, cultivados com manejo semelhante, há oito anos sob sistema plantio direto. A coleta das amostras do solo foi realizada em intervalos regulares de 50 m, totalizando 121 pontos, em duas camadas (0-10 e 10-20 cm). Os dados foram avaliados por estatística descritiva e geoestatística, com base no ajuste de semivariogramas. Verificou-se que a dependência espacial varia conforme o elemento, a textura do solo e a profundidade de coleta da amostra. Assim, o Latossolo de textura média, de maneira geral, apresentou maior variabilidade espacial para os teores dos nutrientes em relação ao de textura muito argilosa. A camada de 0-10 cm, nos dois Latossolos, proporcionou maior variabilidade espacial para os teores dos nutrientes avaliados. A análise de dependência espacial mostrou que, nos dois solos e nas duas camadas, a maioria dos nutrientes estudados apresentou moderada correlação espacial. Os modelos de semivariogramas ajustados foram o exponencial e o esférico, sendo o primeiro em maior quantidade. Maiores teores de todos os macronutrientes avaliados foram verificados nos primeiros 10 cm do solo. Os valores de alcance da dependência espacial foram menores no Latossolo de textura média, variando de 9 a 29 m; no de textura muito argilosa eles variaram de 31 a 399 m. Se for adotada a geoestatística no esquema de amostragem, o número de amostras a serem coletadas será menor no Latossolo de textura muito argilosa, devido aos maiores valores do alcance apresentados por todas as variáveis.

Variabilidade espacial dos teores foliares de nutrientes e da produtividade da soja em dois anos de cultivo em um latossolo vermelho

A análise dos teores de nutrientes de folhas de soja é um dos métodos mais eficientes para avaliar o estado nutricional da lavoura. O objetivo deste estudo foi caracterizar a variabilidade espacial dos teores foliares de nutrientes e da produtividade da cultura da soja num Latossolo Vermelho sob sistema semeadura direta durante dois anos. A área do experimento media 120 x 160 m, totalizando 1,92 ha. As amostras de folhas e os dados de produtividade da soja foram coletados em grade regular de 20 m, totalizando 63 pontos de amostragem nos anos de 1986 e 1988. O teor de nutrientes (N, P, K, Ca, Mg, S, B, Cu, Fe, Mn e Zn) foi determinado analisando-se a terceira folha com pecíolo a partir do ápice de cinco plantas, em locais circunvizinhos de cada um dos pontos de amostragem. As produtividades de soja foram quantificadas em subparcelas de 5 m², sendo posteriormente transformadas para kg ha-1. Os dados foram analisados pela estatística descritiva, a fim de verificar os parâmetros de tendência central e dispersão. A variabilidade espacial foi determinada pelo cálculo do semivariograma e construção de mapas de contorno por meio dos valores obtidos na interpolação por krigagem ordinária. Houve dependência espacial para os teores foliares de alguns nutrientes e para a produtividade de grãos de soja passível de ser mapeada em uma área com cerca de 2 ha, adubada de maneira homogênea. A dependência espacial não ocorreu de forma constante ao longo do tempo, o que deve ser considerado em estudos com cultivos sequenciais. A dependência espacial da produtividade de soja aumentou entre os anos estudados. Entre os nutrientes aplicados anualmente, via adubação, verificou-se a formação de um padrão de distribuição espacial, em 1986 e em 1988, especialmente para os teores de P.