Eficiência de malhas amostrais utilizadas na caracterização da variabilidade espacial de fósforo e potássio

Efficiency of sampling grids used in the characterization of phosphorus and potassium

ARTIGOS CIENTÍFICOS

CIÊNCIA DO SOLO

Eficiência de malhas amostrais utilizadas na caracterização da variabilidade espacial de fósforo e potássio

Efficiency of sampling grids used in the characterization of phosphorus and potassium

Maurício Roberto Cherubin^I,1; Antônio Luis Santi^II Mateus Tonini Eitelwein^I Diego Ricardo Menegol^II Clóvis Orlando Da Ros^II Osmar Henrique de Castro Pias^II Juliano Berghetti^II

^IEscola Superior de Agricultura “Luiz de Queiroz” (ESALQ), Universidade de São Paulo (USP), Avenida Pádua Dias, 11, 13418-900 Piracicaba, SP, Brasil. E-mail: [e-mail]

^IIDepartamento de Ciências Agronômicas e Ambientais, Centro de Educação Superior Norte do Rio Grande do Sul, Universidade Federal de Santa Maria (UFSM), campus de Frederico Westphalen, Frederico Westphalen, RS, Brasil

RESUMO

O objetivo do trabalho foi estudar a eficiência das malhas amostrais, utilizadas nas áreas manejadas com agricultura de precisão, para a caracterização da variabilidade espacial dos teores de fósforo (P) e potássio (K). O estudo foi conduzido em 30 áreas agrícolas, localizadas no Rio Grande do Sul, cujos solos foram classificados como Latossolos Vermelhos. As dimensões das malhas amostrais estudadas foram: 100x100m (10 áreas), 142x142m (10 áreas) e 173x173m (10 áreas). Foram analisados os teores de P e K na profundidade de 0,00-0,10m. Os dados foram submetidos à análise estatística descritiva e a análise geoestatística. De modo geral, as áreas estudadas apresentaram teores médios de P e K adequados ao desenvolvimento das culturas. Entretanto, observaram-se locais com baixos teores de P, justificando fertilizações em sítio-específico. Com base em parâmetros geoestatísticos, as dimensões de malhas amostrais utilizadas são eficientes em captar as diferentes escalas da variabilidade espacial de P e K, porém, recomenda-se o uso de malhas amostrais menores, que considerem a variabilidade espacial a curtas distâncias.

Palavras-chave: agricultura de precisão, amostragem de solo, fertilidade do solo, geoestatística, Latossolos.

ABSTRACT

The objective of the research was to evaluate the efficiency of the sampling grids, used in areas managed with precision agriculture, for the characterization of spatial variability levels of phosphorus (P) and potassium (K). The study was conducted in 30 agricultural areas, located in Rio Grande do Sul state, presenting soils classified as Oxisols. The sizes of the sampling grids studied were: 100x100m (10 areas), 142x142m (10 areas) and 173x173m (10 areas).Were analyzed the levels of P and K, in depth of sampling from 0.00-0.10m. The data were subjected to descriptive statistical analysis and geostatistical analysis. In, general, the areas showed average levels of P and K suitable for crops development. However, were observed subareas with low levels of P, justifying fertilization on site-specific. Based on geostatistical parameters, sizes of the sampling grids used are efficient to capture the different scales of spatial variability of P and K, however, it is indicated smaller sampling grids to consider the spatial variability over short distances.

Key words: precision agriculture, soil sampling, soil fertility, geostatistical, Oxisols.

INTRODUÇÃO

A adoção da agricultura de precisão (AP), na região Sul do Brasil, tem se intensificado nos últimos anos. Estima-se que, somente no Rio Grande do Sul, mais de dois milhões de hectares (ha) utilizem alguma ferramenta de AP (SANTI et al., 2009), especialmente aquelas voltadas ao manejo da fertilidade do solo.

A eficiência do manejo da fertilidade do solo, realizado nas áreas de AP, é diretamente condicionada pela qualidade da amostragem de solo, tradicionalmente realizada de modo sistemático por meio de malhas regulares (KERRY et al., 2010; WEBSTER & LARK, 2012), a qual deve possibilitar captar as diferentes escalas de variabilidade espacial dos atributos do solo (BOTTEGA et al., 2013). Para tanto, a definição da dimensão da malha amostral, necessariamente, deve levar em consideração princípios geoestatísticos, como a dependência espacial dos atributos de solo estudados (WEBSTER & OLIVER, 2007; KERRY et al., 2010). Entretanto, o que se tem observado comercialmente nas áreas manejadas com AP no sul do Brasil é a utilização de malhas regulares, com amostras espaçadas por 100m (uma amostra a cada um ha) a 175m (uma amostra a cada 3ha), definidas basicamente por razões de ordem econômica e prática.

Diversos estudos têm demonstrado que a dependência espacial dos atributos é variável conforme o tipo de solo (MONTANARI et al., 2008), relevo (SOUZA et al., 2006; MONTANARI et al., 2008), atributo estudado (CHERUBIN et al., 2011; NANNI et al., 2011), distância entre amostras (KERRY et al., 2010; NANNI et al., 2011) e uso e manejo do solo (BOTTEGA et al., 2013). Portanto, as diferentes escalas de variação dos atributos do solo induzem grande dificuldade no desenvolvimento de um plano de amostragem, que utilize uma malha amostral com espaçamento único, quando vários atributos do solo estão envolvidos (MONTANARI et al., 2012).

Nesse contexto, objetivou-se estudar a eficiência das malhas amostrais, utilizadas em áreas de Latossolos Vermelhos manejadas com AP, na caracterização da variabilidade espacial do fósforo (P) e potássio (K).

MATERIAL E MÉTODOS

O estudo foi conduzido em 30 áreas agrícolas, localizadas na região Norte do Rio Grande do Sul, Brasil, abrangendo 2.467,82ha. A região estudada está localizada na província geomorfológica do Planalto, com relevo suave ondulado. Os solos predominantes foram classificados como Latossolos Vermelhos de classe textural argilosa, profundos, com alta acidez e baixa fertilidade natural (STRECK et al., 2008).

O manejo adotado nas áreas é o sistema plantio direto, com cultivos rotacionados de soja e milho nas safras de verão e aveia branca (grãos), aveia preta e nabo forrageiro (cobertura do solo) nas safras de inverno.

As áreas foram amostradas a partir das dimensões de malhas amostrais que comumente são utilizadas em nível de campo no sul do Brasil, sendo essas: 100x100m; 142x142m; e 173x173m, que correspondem a uma amostra a cada um, dois e três ha, respectivamente. Assim, cada malha amostral foi utilizada em 10 áreas, totalizando 1.461 amostras de solo.

As amostras foram compostas por 10 subamostras, coletadas em um raio de 10 m dos pontos georreferenciados, utilizando um trado de rosca acoplado a um quadriciclo, na profundidade de 0,00-0,10m. Esse procedimento de coleta de solo segue o que está sendo realizado atualmente nas áreas de AP da região estudada. Posteriormente, os teores de P (mg dm^-3) e K (mg dm^-3) foram extraídos por Mehlich-1, seguindo metodologias recomendadas pela CQFS-RS/SC (2004). Quanto aos parâmetros interpretativos de P e K, as áreas apresentaram teores médios de argila classificados na classe 2 (40-60%) e capacidade de troca de cátions potencial (CTC_pH7,0) na classe média (5,1-15,0_cmolcdm^-3).

Os dados foram submetidos à análise estatística descritiva no programa Statistical Analysis System – SAS 8.0 (SAS Inc, Cary, EUA). A análise da variabilidade espacial dos dados de P e K foi realizada por meio de semivariogramas, cujos ajustes foram realizados por modelos matemáticos teóricos, utilizando o programa computacional GS+ (Gamma Design Software, LLC, Plainwell, EUA). Os modelos dos semivariogramas foram ajustados com base na menor soma de quadrados do resíduo (SQR) e no melhor coeficiente de determinação (r²), e avaliados pela técnica de validação cruzada. Do ajuste de um modelo matemático aos dados, foram definidos os parâmetros do semivariograma: efeito pepita (C₀), contribuição (C₁), patamar (C₀+C₁) e alcance (A). O grau de dependência espacial (GDE) foi calculado por meio da relação C₀/(C₀+C₁)*100, classificando-o como forte, para GDE ≤25%; moderado, para GDE entre 25 e 75%, e fraco, para GDE >75% (CAMBARDELLA et al., 1994).

RESULTADOS E DISCUSSÃO

A análise estatística descritiva dos valores de P (Tabela 1) possibilitou verificar que, em geral, os teores médios encontram-se adequados, entretanto, a maioria das áreas ainda apresenta locais com P abaixo do teor crítico (<9mg dm^-3 para a classe 2 de argila), sendo considerados limitantes ao desenvolvimento das culturas pela CQFS-RS/SC (2004).

Os valores de P apresentaram grande amplitude, com teores máximos de 2,32 a 18,57 vezes maiores que os valores mínimos, indicando elevada potencialidade do manejo em sítio-específico preconizado na AP (AMADO et al., 2009). A variação dos teores é comprovada pelos altos valores de CV (>35% em 63,34% das áreas) e pela distribuição não normal dos dados, em 73,34% das áreas. Altas variações nos teores de P concordam com diversos estudos em Latossolos reportados na literatura (MONTANARI et al., 2008; AMADO et al., 2009; CHERUBIN et al., 2011; NANNI et al., 2011; BOTTEGA et al., 2013). De acordo com muitos dos autores acima citados, a elevada variação dos teores de P decorre das sucessivas aplicações de fertilizantes a taxa fixa na linha de semeadura, contribuindo para a manutenção e, ou, aumento da micro e macrovariabilidade espacial nas áreas.

Em relação ao K (Tabela 2), verificou-se que as áreas estudadas apresentaram teores acima do nível crítico (60mg dm^-3), classificados como altos (61 a 120mg dm^-3) e muito altos (>120mg dm^-3) (CQFS-RS/SC, 2004). Altos teores de K também foram verificados por AMADO et al. (2009) e CHERUBIN et al. (2011) em Latossolos no sul do Brasil, e estão atrelados ao intenso intemperismo dos minerais primários (micas, feldspatos e feldspatoides) e secundários (illita, vermiculita e argilominerais interestratificados) que contém este elemento (ERNANI et al., 2007), associado ao histórico de sucessivas adubações, utilizando formulações com elevadas quantidades de K adicionadas às culturas.

Os teores máximos de K foram de 1,85 a 6,74 vezes maiores que os valores mínimos. Essa menor variação, em relação à observada no P, é refletida pelos valores de CV, entre 15-35% em 93,33% das áreas, e pela distribuição de frequências normais, observadas em 56,67% das áreas estudadas. As menores variações do K, também verificadas por AMADO et al. (2009), CHERUBIN et al. (2011) e BOTTEGA et al. (2013), podem ser atribuídas à maior mobilidade do K no solo (ERNANI et al., 2007) em relação ao P, minimizando a microvariabilidade a partir da concentração de fertilizantes potássicos na linha. Além disso, a rápida liberação do K presente nos restos culturais, dispostos uniformemente na área, contribui na redução da variabilidade desse atributo no solo.

Efetuada por meio dos semivariogramas, a análise geoestatística permitiu detectar que as áreas estudadas apresentaram diferentes escalas de variabilidade espacial para o P (Tabela 3) e para o K (Tabela 4), mesmo quando amostradas com a mesma dimensão de malha. Essa constatação é comprovada pelas distâncias em que os atributos apresentaram continuidade espacial (alcances), que variaram de 273 a 3.110m para o P e de 315 a 3.280m para o K, com grau de dependência espacial classificado de moderado a forte.

Em todas as malhas estudadas, verificaram-se áreas em que os valores de P e K não apresentaram dependência espacial, caracterizando distribuições aleatórias (efeito pepita puro - EPP). O EPP pode ocorrer devido a erros de medidas, de amostragem ou microvariações não detectadas, considerando que o espaçamento de amostragem utilizado é maior que o necessário para detectar dependência espacial (CAMBARDELLA et al., 1994; WEBSTER & OLIVER, 2007). Nesse sentido, observou-se que há uma tendência de aumento na ocorrência de distribuições aleatórias, à medida que aumenta a dimensão da malha amostral, confirmando que as malhas amostrais utilizadas possuem pontos muito espaçados e/ou em número insuficientes para detectar a dependência que, se existir, será manifestada a distâncias menores que o menor espaçamento utilizado entre amostras.

A não detecção da dependência espacial, ou o comportamento aleatório do P e K, observado em 33,33%, das áreas estudadas, deve-se ao fato de que esses atributos são altamente afetados por práticas de manejo (adubações em linha, exportações via colheita, sistema plantio direto e cultura cultivada), condicionando a variações mais pronunciadas em curtas distâncias (NANNI et al., 2011). Resultados semelhantes foram obtidos por CORÁ & BERALDO (2006) e NANNI et al. (2011), demonstrando que a coleta de uma amostra por ha não foi suficiente para captar a variabilidade espacial de P e K, indicando a necessidade do uso de malhas mais adensadas para recomendar fertilização com maior acurácia. Os resultados obtidos neste estudo também concordam com a realidade verificada nas áreas comerciais que adotam a AP nos Estados Unidos, em que as amostragens de solo, comumente utilizando malhas de 100x100m (uma amostra por ha), não são eficientes na detecção das diferentes escalas de variação dos atributos do solo, em especial de P e K (KERRY et al., 2010).

A partir dos resultados, poder-se-ia inferir que as malhas amostrais utilizadas possibilitaram captar a variabilidade espacial do P e K para a maioria das áreas estudadas, entretanto, tal afirmação merece cautela, uma vez que os semivariogramas foram gerados utilizando um limitado número de pontos (n), o que pode comprometer a consistência dos resultados. De acordo com LANDIM (2006), para utilizar a geoestatística, indica-se que sejam adotados no mínimo 30 pares para cada distância estudada.

Portanto, de acordo com os resultados e argumentações apresentadas, afirma-se que as dimensões de malhas amostrais, comumente utilizadas nas áreas de Latossolos manejados com AP no sul do Brasil, de modo geral, possibilitaram captar as diferentes escalas da variabilidade espacial de P e K. Apesar disso, recomenda-se evitar a utilização de malhas amostrais que condicionem a obtenção de poucas amostras nas áreas, uma vez que limita a acurácia da caracterização e o manejo da variabilidade espacial dos atributos analisados.

Ressalta-se também que a caracterização e manejo da microvariabilidade espacial dos teores de P e K, observada em parte das áreas e condicionada pelas práticas de manejo adotadas, exigem a utilização de malhas amostrais menores (mais de uma amostra por ha) do que as utilizadas atualmente, associando, assim, menores distâncias entre amostras e maior número de amostras coletadas na área. Dessa forma, não é possível extrapolar regionalmente um modelo confiável de variabilidade espacial destes atributos. Entretanto, o estudo fornece uma noção geral da escala de variação de P e K, constituindo-se um importante referencial que poderá ser considerado no planejamento de futuras estratégias de amostragem de solo a serem adotadas nas áreas de AP da região estudada.

Os resultados demonstraram que a acurácia dos mapeamentos da variabilidade espacial de P e K no solo e, consequentemente, o sucesso de intervenções localizadas em sítios-específicos depende da utilização de malhas amostrais adensadas, que possibilitem coletar um número de amostras suficiente para captar as diferentes escalas de variabilidade que ocorrem na área. Dessa forma, verifica-se que um dos grandes desafios da AP está na busca por metodologias alternativas, que permitam obter mapas acurados de atributos do solo, gerados a partir de maior número de informações possíveis, em menor tempo e com baixo custo.

Nesse sentido, o uso de sensores baseados em técnicas, como a espectroscopia do infravermelho próximo (NIR) e visível (VIS), tem apresentado resultados promissores na estimativa de atributos do solo, como o P e o K (KODAIRA & SHIBUSAWA, 2013) e tem sido apontado com um dos principais direcionamentos da AP.

CONCLUSÃO

As dimensões de malhas amostrais utilizadas nas áreas manejados com AP, sob procedimentos geoestatísticos, são eficientes para captar as diferentes escalas da variabilidade espacial dos teores de P e K no solo. Entretanto, recomenda-se evitar a utilização de malhas amostrais que condicionem a obtenção de poucas amostras nas áreas e, consequentemente, a caracterização e manejo da variabilidade espacial dos atributos analisados com menor acurácia.

A variabilidade espacial expressa em curtas distâncias deverá ser considerada em futuros planos de amostragem nas áreas manejadas com AP, indicando a necessidade da utilização de malhas amostrais menores (mais uma amostra por ha) em relação às utilizadas atualmente.

AGRADECIMENTOS

Ao Programa de Pós-graduação em Agronomia: Agricultura em Ambiente da Universidade Federal de Santa Maria (UFSM). A Coordenação Aperfeiçoamento de Pessoal de Nível Superior (CAPES) pela concessão da bolsa ao primeiro autor.

Recebido 01.06.13

Aprovado 19.09.13

Devolvido pelo autor 04.12.13

CR-2013-0771.R1

• AMADO, T.J.C. et al. Atributos químicos e físicos de Latossolos e sua relação com os rendimentos de milho e feijão irrigados. Revista Brasileira de Ciência do Solo, v.33, n.4, p.831-843, 2009. Disponível em: <http://www.scielo.br/pdf/rbcs/v33n4/08.pdf>. Acesso em: 20 nov. 2012. doi:10.1590/S0100-06832009000400008.
• BOTTEGA, E.L. et al. Variabilidade espacial de atributos do solo em sistema de semeadura direta com rotação de culturas no cerrado brasileiro. Revista Ciência Agronômica, v.44, n.1, p.1-9, 2013. Disponível em: <http://www.scielo.br/pdf/rca/v44n1/a01v44n1.pdf>.Acesso em: 2 jan. 2013.doi: 10.1590/S1806-66902013000100001.
• CAMBARDELLA, C.A. et al. Fieldscale variability of soil properties in central Iowa soils. Soil Science Society of America Journal, v.8, n.6, p.1501-1511, 1994. Disponível em: <http://ddr.nal.usda.gov./bitstream/10113/18168/1/IND20518015.pdf>. Acesso em: 12 set. 2009.
• CHERUBIN, M.R. et al. Caracterização e estratégias de manejo da variabilidade espacial dos atributos químicos do solo utilizando a análise dos componentes principais. Enciclopédia Biosfera, v.7, p.196-210, 2011. Disponível em: <http://www.conhecer.org.br/enciclop/2011b/ ciencias%20agrarias/caracterizacao%20e%20estrategias.pdf>. Acesso em: 12 dez. 2011.
• COMISSÃO DE QUÍMICA E FERTILIDADE DO SOLO (CQFS-RS/SC). Manual de adubação e calagem para os estados do Rio Grande do Sul e Santa Catarina Porto Alegre: Sociedade Brasileira de Ciência do Solo - Núcleo Regional Sul, 2004. 400p.
• CORÁ, J.E.; BERALDO, J.M.G. Variabilidade espacial de atributos do solo antes e após calagem e fosfatagem em doses variadas na cultura de cana-de-açúcar. Engenharia Agrícola, v.26, n.2, p.374-387, 2006. Disponível em: <http://www.scielo.br/pdf/eagri/v26n2/05.pdf>. Acesso em: 10 set. 2012. doi:10.1590/S0100-69162006000200005.
• ERNANI, P.R. et al. Potássio. In: NOVAIS, R.F. et al. (Eds.). Fertilidade do solo Viçosa-MG: Sociedade Brasileira de Ciência do Solo, 2007. p.551-594.
• KERRY, R. et al. Sampling in precision agriculture. In: OLIVER, M.A. (Org.). Geostatistical applications for precision agriculture Heidelberg: Springer-Verlag, 2010. p.35-63.
• KODAIRA, M.; SHIBUSAWA, S. Using a mobile real-time soil visible-near infrared sensor for high resolution soil property mapping. Geoderma, v. 199, n.5, p.64-79, 2013. Disponível em: <http://www.sciencedirect.com/science/article/pii/S0016706112003412>. Acesso em: 27 nov. 2012. doi: http://dx.doi.org/10.1016/j.geoderma.2012.09.007.
• LANDIM, P.M.B. Sobre geoestatística e mapas. Terræ Didatica, v.2, n.1, p.19-33, 2006. Disponível em: <http://www.ige.unicamp.br/terraedidatica/v2/pdfv2/t_didatica_2006_ v02n01_p019-033_landim.pdf>. Acesso em: 9 dez. 2012.
• MONTANARI, R. et al. Variabilidade espacial de atributos químicos em Latossolo e Argissolos. Ciência Rural, v.38, n.5, p.1266-1272, 2008. Disponível em: <http://www.scielo.br/pdf/cr/v38n5/a10v38n5.pdf>. Acesso em: 11 nov. 2012. doi:10.1590/S0103-84782008000500010.
• MONTANARI, R. et al. The use of scaled semivariograms to plan soil sampling in sugarcane fields. Precision Agriculture, v.13, n.5, p.542-552, 2012. Disponível em: <http://link.springer.com/content/pdf/10.1007%2Fs11119-012-9265-6.pdf>. Acesso em: 10 out. 2012.doi:10.1007/s11119-012-9265-6.
• NANNI, M.R. et al. Optimum size in grid soil sampling for variable rate application in site-specific management. ScientiaAgrícola, v.68, n.3, p.386-392, 2011. Disponível em: <http://www.scielo.br/pdf/sa/v68n3/v68n3a17.pdf>. Acesso em: 18 dez. 2012. doi: 10.1590/S0103-90162011000300017.
• SANTI, A.L. et al. É chegada a hora da integração do conhecimento. Revista Plantio Direto, v.129, n.1, p.24-30, 2009. Disponível em: <http://www.plantiodireto.com.br/?body=cont_ int&id=907>. Acesso em: 18 nov. 2012.
• SOUZA, Z.M. de. et al. Otimização amostral de atributos de latossolos considerando aspectos solo-relevo. Ciência Rural, v.36, n.3 p.829-836, 2006. Disponível em: <http://www.scielo.br/pdf/cr/v36n3/a16v36n3.pdf>. Acesso em: 25 nov. 2011. doi: 10.1590/S0103-84782006000300016.
• STRECK, E.V. et al. Solos do Rio Grande do Sul 2.ed. Porto Alegre: EMATER/RS-ASCAR, 2008. 222p.
• WEBSTER, R.; LARK, M. Field sampling for environmental science and management London: Routledge. 2012. 200p.
• WEBSTER, R.; OLIVER, M.A. Geostatistics for environmental scientists 2.ed. Chichester: John Wiley& Sons, 2007. 330p.

1

Autor para correspondência

Datas de Publicação

• Publicação nesta coleção
  18 Mar 2014
• Data do Fascículo
  Mar 2014

Histórico

• Aceito
  19 Set 2013
• Recebido
  01 Jun 2013