Resumo: A Amazônia tem despertado a atenção dos produtores no plantio de açaizeiros, muitos deles, buscam inovar técnicas de cultivo que possam auxiliar nas tomadas de decisões na agricultura, visando desenvolver sistemas de cultivo apropriados para aumentar a produtividade. Objetivo deste trabalho foi avaliar a variabilidade espacial dos atributos químicos do solo e sua interação com o cultivo do açaí (Euterpe oleracea Mart.), localizada na fazenda Fonte Nova no Município do Acará-PA. Os dados foram submetidos a estatísticas descritivas e teste de normalidade de Shapiro-Wilk. Para os elementos químicos do solo avaliou-se ainda a dependência espacial, geração modelos dos semivariogramas e mapas de krigagem. Por fim foi realizada uma matriz de correlação entre os atributos estudados. Os resultados mostraram que a maior parte das variáveis químicas do solo apresentou dependência espacial, indicando que o número de amostras coletadas e os pontos plotados foram suficientes para mapear área do plantio, por meio das técnicas geoestatísticas. Os modelos teóricos predominantes no estudo foram esférico, gaussiano e o exponencial. As variáveis apresentaram moderada correlação. Percebeu-se que os mapeamentos dos atributos do solo e do plantio foram importantes para identificação das áreas com diferentes níveis de fertilidades e as regiões com maiores potenciais produtivos.
Palavras-chave: Geoestatística, correlação espacial, agricultura de precisão.
Abstract: The Amazon has awakened the attention of producers in the planting of açaí trees, many of them seeking to innovate cultivation techniques that can assist in decision-making in agriculture, aiming to develop appropriate cultivation systems to increase productivity. The objective of this work was to evaluate the spatial variability of soil chemical attributes and its interaction with the cultivation of açaí (Euterpe oleracea Mart.). Located on the Fonte Nova farm in the Municipality of Acará-PA. The data were submitted to descriptive statistics and normality test by Shapiro-Wilk. For soil chemical elements, spatial dependence, generation of semivariograms and kriging maps were also evaluated. Finally, a correlation matrix was performed between the studied attributes. The results showed that most of the chemical variables of the soil presented spatial dependence, indicating that the number of samples collected and the plotted points were sufficient to map the planting area, using geostatistics techniques. The theoretical models predominant in the study were spherical, Gaussian and exponential. The variables showed a moderate correlation. It was noticed that the mapping of soil and planting attributes were important to identify areas with different levels of fertility and regions with greater productive potential.
Keywords: Geostatistics, spatial correlation, precision agriculture.
Artigo
Dispersão Espacial de Atributos Químicos do Solo de um Açaizeiro na Região Amazônica
Spatial Dispersion of Chemical Attributes of the Soil of an Açaí Tree in the Amazon Region
Recepção: 16 Janeiro 2021
Aprovação: 15 Dezembro 2021
A floresta Amazônica é caracterizada pela sua alta diversidade biológica. Dentro dessa diversidade, destacam-se as plantas da família Arecaceae denominadas popularmente de “palmeiras”, no qual podem se extrair diversos produtos, tanto do fruto como em gêneros derivados da planta. O açaí é um fruto típico e popular da região amazônica, ele é comercializado e consumido pela população local nas regiões produtoras do Pará, Maranhão, Amapá, Acre e Rondônia (Menezes, Torres & Srur 2008). Com o crescimento do mercado, o açaí passou a ser consumido também nas grandes capitais brasileiras, assim como a exportação que passou a ser feita para Estados Unidos, Japão, China e alguns países da Europa (Silva 2002). Essa demanda nacional e internacional por suco de açaí tem levado a uma grande expansão do cultivo em áreas irrigadas (Lindolfo et al. 2020).
Os produtores da mesorregião do Nordeste Paraense, polo de produção de açaizeiros em terra firme, vêm desenvolvendo técnicas inovadoras de cultivo, visando aumentar a produtividade, tanto na safra como na entressafra (Homma et al. 2014). Atualmente para maior precisão na tomada de decisão na agricultura, é importante que os processos de manejo e os sistemas de produção sejam acompanhados de inovações tecnológicas, como a agricultura de precisão.
Os fundamentos da agricultura de precisão se baseiam na aplicação de tecnologias de maneira localizada de acordo com a variabilidade espacial (Grego, Oliveira & Vieira 2014). Diretamente relacionada à agricultura de precisão, está a geoestatística que possibilita compreender a variabilidade espacial, promover a interpretação dos resultados com base no fundamento da variabilidade natural dos atributos analisados e consegue representar os mais diversos tipos de variações (Souza, Marques & Pereira 2010).
Dentre os atributos que podem ser avaliados neste tipo de estudo está o solo e planta, que fornecem uma série de informações que podem ser utilizadas em conjunto, visando o aumento da produtividade, otimização no uso dos recursos, fertilizantes e redução no impacto ambiental (Dalchiavon et al. 2017).
O conhecimento da variação de atributos químicos e sua relação com a planta é importante para o manejo e planejamento de esquemas de amostragem e gerenciamento de práticas agrícolas (Li et al. 2012). As interações existentes no sistema solo são fundamentais para sua compreensão e decisivas para o desenvolvimento de culturas agrícolas. Desta forma, a avaliação da variabilidade espacial dos atributos tem se tornado importante ferramenta, na determinação de estratégias de manejo do solo que buscam aumentar a produtividade agrícola (Minasny & McBratney 2007).
O estudo do comportamento da variabilidade espacial dos atributos químicos do solo com plantio, por meio da técnica de geoestatística, pode otimizar a aplicação localizada de corretivos e fertilizantes. Nesse sentido, o objetivo deste trabalho foi avaliar a variabilidade espacial dos atributos químicos do solo e seus efeitos na área plantio de açaizeiro (Euterpe oleracea Mart.) na Fazenda Fonte Nova (nordeste do Pará) visando o controle do sistema de produção do açaizeiro.
A área de desenvolvimento desse trabalho foi no município de Acará que pertence à mesorregião do nordeste Paraense e a microrregião de Tomé-Açú, na Fazenda Fonte Nova, próximo à rodovia PA 256, km 45. O local possui cultivo de açaí (Euterpe oleiracea Mart.) e é dividido em 4 talhões, totalizando de 16,4ha entre as coordenadas geográficas: 02º 27' 00'' e 02º 27' 28,8'' de latitude Sul, 48º 34' 26,4'' e 48º 34' 40,8'' de longitude Oeste.
A área de estudo foi destinada ao pasto em 1979, mas a partir de 2008 deu-se o início a preparação da terra para a implantação das mudas, com espaçamento 5 x 5 m. A primeira produção foi em 2012 com 5.665,98 kg/ha, e para os anos de 2013, 2014 e 2015 a produção foi de 8.909,55 kg/ha, 8.242,38 kg/ha e 10.033,78 kg/ha, respectivamente. No decorrer de cada ano, foram efetuadas três aplicações de adubação de produção e o processo de irrigação é do tipo microaspersão, com 120 litros de água/touceira/dia, ficando ligada por um período de 2 horas em cada talhão.
O solo presente na área é do tipo Latossolo amarelo, textura argilosa e média, clima Afi, correspondente ao tropical úmido, com temperatura mensal superior a 18ºC, precipitação acima de 60 mm, em relação às aferições mensais são equivalentes a 2.837 mm e amplitude térmica anual em torno de 5ºC, segundo a classificação de Köppen, (IDESP 2014).
Para a coleta dos dados a área foi dividida em uma malha regular de 60 x 60 m determinando assim 38 pontos amostrados sistematicamente, sendo que em cada ponto foi implantada uma parcela medindo 30 x 30 m, totalizando 38 parcelas, ilustrada na Figura 1. No centro de cada parcela foi coletado amostragem dos atributos químicos do solo e a sua posição espacial, utilizando-se um receptor de Sistema de Posicionamento Global (GPS), configurado no sistema de coordenadas métricas UTM- SIRGAS 2000. As amostras do solo foram coletadas um mês após a adubação de produção, pois é o tempo necessário para obter um resultado mais confiável das condições químicas do solo (Borges & Accioly 2007; Brasil, Cravo & Veloso 2020). Para realização da coleta utilizou-se o trado do tipo sonda para obtenção cinco amostras simples de solo na camada de 0-20 cm de profundidade que posteriormente constituíram uma amostra composta para cada parcela. O processo de coleta foi conforme recomendação da EMBRAPA (2009) e enviada para análises no laboratório de solos da Universidade Federal Rural da Amazônia.

As análises químicas do solo foram realizadas no laboratório da Universidade Federal Rural da Amazônia e seguiram as normas do manual de métodos de análise de solos, (EMBRAPA 2009). Primeiramente as amostras foram submetidas ao processo, terra fina seca ao ar (TFSA) e posteriormente com o término do condicionamento, foi feito a extração das seguintes informações: pH em H2O (água), Potássio (K), Fósforo (P), Cálcio (Ca), Cálcio + Magnésio (Ca+Mg), Alumínio (Al), Acidez potencial (H+Al), Carbono orgânico (C) e Matéria Orgânica (M.O). A partir desses resultados foram calculados a soma de bases (SB), Saturação de bases (V%), Saturação por alumínio (m%), Capacidade de troca de cátions efetiva (CTCefetiva) e Capacidade de troca de cátions efetiva em pH7 (CTC pH7). Conforme as equações abaixo.
Para a amostragem dos atributos da planta foram coletadas na área de estudo, informações de 10 plantas por parcela, situadas o mais próximo do ponto amostral da malha. Em cada planta foi obtida os seguintes dados: altura, circunferência (CIR), número de cachos no estipe principal (NCEP) e quantidade de frutos maduros em ponto de colheita no estipe principal (QFMEP), apresentada na Figura 2. As medidas foram feitas de forma manual, através do uso de fita métrica, uma vara métrica e anotações das observações visuais em campo. Essa metodologia é similar aos estudos de Gama (2012), Aguiar (2016) e Mendonça & Moreira (2020).

Os dados obtidos foram submetidos a uma análise exploratória pelo método da estatística descritiva, obtendo-se média, mediana, valores mínimos e máximos, desvio padrão, coeficiente de variação, curtose, assimetria e coeficiente de variação. A hipótese de normalidade foi realizada através do teste de Shapiro-Wilk a 5% de probabilidade. Para as variáveis que não apresentaram distribuição normal, procedeu-se a transformação logarítmica dos dados verificando se estes seguiam a distribuição log-normal. Para obtenção dos cálculos utilizou-se o software livre BioEstat, versão 5.0.
Os dados de atributos químicos do solo foram submetidos a análise geoestatística, no qual a dependência espacial foi analisada por meio de ajustes de semivariogramas, que realiza os cálculos das semivariâncias amostrais, expressa pela equação abaixo:
Em que N(h) representa o número de pares de valores determinados, Z(xi) e Z(xi + h), separados por uma distância estabelecida (h). Os valores de Z podem ser qualquer um dos parâmetros analisados, enquanto os valores de xi e xi + h são determinados de acordo com as posições dos pontos amostrais no campo.
Após o cálculo do semivariograma, os valores de semivariância y(h) e da distância (h) foram dispostos em gráficos de dispersão, tendo como valores de Y as semivariâncias, e de X, as distâncias. Ajustou-se o modelo que melhor representou a relação entre semivariância experimental e a distância h, e estabeleceu os parâmetros: efeito pepita (C0), patamar (C0 + C) e alcance (A).
Para a análise do grau de dependência espacial das variáveis foi utilizada a classificação de (Cambardella et al. 1994), a fim de mensurar se a dependência espacial é fraca, forte ou moderada, para isso basta verificar a dependência espacial que existe entre o efeito pepita e o patamar através do índice de dependência espacial (IDE), se o IDE ≤ 25% é considerado forte dependência espacial, para o intervalo de o 25%≤ IDE ≤75% significa que há moderada dependência e acima de 75% ocorre fraca dependência espacial.
Os mapas de distribuição espacial dos atributos químicos do solo foram elaborados por interpolação dos dados (krigagem ordinária). Em seguida, os dados desenvolvimentos do plantio, referentes as 38 parcelas, foram convertidos do formato tabulado (Excel) para pontos geográficos vetoriais, com intuito para executar uma análise espacial. Por fim, após a obtenção das superfícies dos dados, foi aplicada uma matriz de correlação, para verificar a existência correlação entres as variáveis em estudo. Para realização dessas etapas de processamento de dados foi utilizado o software ArcGis, versão 10.1.
Os resultados dos dados de desenvolvimento da produção (Tabela 1) indicam que o teste de Shapiro-Wilk e a estatística descritiva confirmam a distribuição de normalidade dos dados, no qual observa-se pelo coeficiente de variação a baixa dispersão e pelos valores da média e mediana de que quase todas as variáveis são próximas entre si, caracterizando distribuição tendendo a simetria. Com exceção da altura da planta, no qual observou-se pelos valores de assimetria e curtose que tinham tendência a ser uma curva mais à esquerda e aberta que normal.

Para as variáveis químicas do solo, observou-se que os valores de média e mediana estão próximos, evidenciando a simetria, confirmado pelos valores de assimetria e curtose próximos de zero. Com exceção das variáveis: potássio (K) e Fósforo (P), que apresentam uma distribuição mais acentuada. A hipótese da normalidade pode-se confirmada pelo teste de Shapiro-Wilk (Tabela 1). Apenas algumas variáveis do solo se ajustaram à distribuição de normalidade, estes apontam que: pH em H20, Acidez potencial (H+Al), Matéria Orgânica (M.O), Cálcio (Ca), CTC pH 7 e Saturação por alumínio (m%) apresentam distribuição normal. Entretanto, algumas variáveis não apresentaram distribuição normal, teve-se que proceder a transformação log-normal, os atributos químicos foram: Fósforo (P), Cálcio + Magnésio (Ca+Mg), Alumínio (Al), Soma de Bases (SB) e Saturação por Bases (V%). Apesar do procedimento de transformação, o elemento químico potássio (K) foi o único que não obteve normalidade. De maneira geral, as distribuições de frequência do tipo normal para os demais atributos químicos do solo ficaram, em concordância observada pelo autor (Matias et al. 2015).
Os coeficientes de variação dos elementos químicos do solo, indicam que somente o atributo químico pH em H20 apresentou baixa variação, caracterizando uma baixa variabilidade dos dados, conforme proposta feita por Gomes (1985), que classifica os valores como: baixos (menores que 10%), médios (entre 10 e 20%), representando pouca variabilidade para essas variáveis, entre 20 e 30% considerado alto, e acima de 30% muito alto, representando alta variabilidade. As demais variáveis químicas ficaram acima de 10%, indicando presença de médios ou altos valores discrepantes. Um alto coeficiente de variação (CV) é o primeiro indicador da heterogeneidade dos dados de uma variável, o que pode estar relacionado a efeitos de formação e manejo do solo (Silva et al. 2020).
Na Tabela 2, mostra-se a análise geoestatística dos atributos químicos do solo, no qual as variáveis: potássio (K), Alumínio (Al), Matéria Orgânica (M.O) e Saturação por alumínio (m%) não apresentaram dependência espacial. Para estes dados, tem-se uma distribuição completamente aleatória, não tendo sob o ponto de vista geoestatístico, uma função que descreva sua variação no espaço. Isso denota que o número de amostras utilizadas para o estudo desses parâmetros não foi suficiente para detectar dependência espacial e, nesses casos, uma maior densidade de amostras poderia ser considerada para uma reavaliação de sua dependência espacial (Silva et al. 2020).
As variáveis: P, Ca, Ca + Mg, H + Al, pH em H2O, CTC pH7, V% e SB apresentaram estrutura espacial, assim sendo, esses elementos foram ajustados aos modelos teóricos. Observou-se que Fósforo (P) se ajustou melhor ao modelo exponencial, enquanto o Ca+Mg e Soma de Bases (SB) mostram que estas variáveis têm uma melhor adaptação ao modelo gaussiano, já as variáveis Cálcio (Ca), Acidez potencial (H+Al), pH em H2O, CTC pH7 e Saturação por Bases (V%) se ajustam melhor ao modelo teórico ajustado tipo esférico, atingindo assim uma correlação espacial adequada. Os resultados estão em concordância com Matias et al. (2015), confirmando assim, que os modelos esférico, gaussiano e exponencial têm sido os mais frequentes para estudos variabilidade espacial de atributos químicos do solo.

Com base na metodologia proposta por (Cambardella et al. 1994), variaram de moderado a forte dependência espacial (Tabela 2), o maior alcance constatado foi para o Ca + Mg com 339m, indicando maior continuidade espacial e já o menor alcance foi para o fósforo com 171 m. Os valores de alcance influenciam diretamente na estimativa. Os atributos que apresentam maior alcance de dependência espacial tendem a se apresentar mais homogêneos espacialmente (Santos et al. 2016).
De todas as variáveis em estudo, o cálcio (Ca) representou melhor explicação para a variabilidade espacial, cujo melhor modelo ajustado foi o esférico, que predomina nos trabalhos em ciência do solo (Artur et al. 2014; Malvezi et al. 2019). Essa variável, apresentou um raio de 334 m (alcance), que é referente a uma área de 350.463 m², no espaçamento de 5 x 5 m, com uma aleatoriedade de 11%, significando uma forte correlação com o espaço, ou seja, 89% da variabilidade é explicada pela correlação espacial.
Observa-se nos mapas de krigagem (Figura 3) a distribuição dos atributos químicos do solo foi de grande variação na área, confirmando os resultados preliminares da análise descritiva (Tabela 1). Esses resultados são fundamentais, pois mostram as menores concentrações dos atributos: Cálcio + Magnésio, Fósforo (P), Cálcio (Ca) e Soma de Bases (SB) na borda inferior do plantio, mesmo após a adubação de produção, evidenciando assim a aplicação do nutriente em questão no decorrer da área de cultivo, a taxas variadas, visando à homogeneização da fertilidade do solo, minimizando as chamadas manchas de fertilidade do solo (Carneiro et al. 2017).
Essa relação do cultivo com os elementos químicos do solo pode-se observado na Figura 3, na qual mostra o crescimento e produtividade em diferentes escalas do plantio sob a distribuição dos nutrientes do solo presente na área em concentrações diferenciadas. Essa falta ou excesso dos atributos em determinadas regiões do plantio, pode interferir na produtividade total da área. Contudo, vale ressaltar que outros fatores, como: pragas, doenças, competição por nutrientes, os atributos físicos do solo, a luminosidade, tipo de trato de manejo adotado, número de perfilhos e própria biologia da planta influenciam no desenvolvimento e na produtividade (Jardim & Anderson 1987; Gama 2012).

Essa associação das variáveis em estudo se faz necessária para entender as razões dessas variações, cujo é reflexo de vários das interações, tais como: formação do solo, manejo do solo e da cultura. A Tabela 3 apresenta a matriz de correlação entre os atributos químicos do solo e do plantio em que observou-se, que no geral, os atributos físicos da planta (Altura, Circunferência, NCEP e QFMEP), apresentaram moderada a fraca correlação com os elementos químicos do solo, o que pode ser justiçado pelas manchas de baixa fertilidade do solo presente na área, apesar da adubação de produção. Além disso, o solo em questão caracterizou-se como ácido o que pode interferir na absorção de nutrientes.
Para o Cálcio (Ca), em específico, há correlação direta com a (SB) e o (V%), pelo fato do Ca fazer parte direta no cálculo da (SB) e no cálculo da (V%), uma vez que o Cálcio é o principal elemento em termos quantitativos. Observa-se, ainda, que a relação linear entre SB e V% é de 0,89, essa elevada correlação positiva entre os dois se deve ao valor de V % pois este é expresso pela relação entre SB e CTC. Outra observação feita foi o grau de associação entre o Cálcio (Ca) e o pH em H2O ser de (0,89), pelo fato do Ca ser um nutriente catiônico básico trocável, ou seja, a disponibilidade do elemento, depende de uma maior basicidade do solo.
Diante disso, conhecer a distribuição espacial do elemento em estudo pode servir como importante ferramenta para saber a localização de áreas com maior e menor disponibilidade de nutrientes. A partir dessas análises da variabilidade espacial das variáveis químicas do solo de acordo com a quantidade presente no solo pode-se efetuar um melhor manejo a área do plantio de açaizeiro, como adubações, processo de calagem e aplicação de fertilizante em diferente dosagem.

Os modelos de semivariogramas mais indicados para os atributos químicos foram: o esférico, gaussiano e o exponencial.
Para as condições de estudo, o Cálcio (Ca) indicou grau de dependência espacial forte e o melhor coeficiente de determinação (R²= 1), indicando assim que 89% da variabilidade é explicada pela correlação espacial.
Ao analisar os dados apurados através dos mapas, observou-se que a aplicação de adubos e corretivos, pode ser aplicada de forma localizada conforme e necessidade do solo, diminuindo assim, a perda de fertilizantes agrícolas.
Agradecemos aos proprietários da Fazenda Fonte Nova, aos técnicos do laboratório de análise química do solo da Universidade Rural da Amazônia - UFRA e a toda equipe da empresa EKOS, pelo apoio logístico, profissional e pessoal, o qual foi extremamente importante para a prática do trabalho.
E-mails:luceliasbarros@gmail.com; ewelynrocha@gmail.com; nazare.maciel@ufra.edu.br; vania.melo@ufra.edu.br; rejane_ig@unb.br; tati_almeida@unb.brAutor correspondente: Lucelia Souza de Barros; luceliasbarros@gmail.com





