Introducción

O Sistema de Plantio Direto (SPD) propõe uso contínuo de cobertura do solo na lavoura, através da deposição de restos vegetais de culturas advindas do processo de rotação cultural (¹). A palha na superfície do solo oferece diversos benefícios, como a dissipação da energia do impacto da gota de chuva, aumento da infiltração de água, manutenção da umidade, melhoria da macro e micro fauna, redução da amplitude térmica, menores perdas por erosão e diminuição da incidência de plantas daninhas (²), bem como, amortecimento da pressão imputada no tráfego.

Anualmente o aporte de palha para o sistema de plantio direto deve ser de 10 a 12 t ha-1 (³), nesse sentido, uma das etapas fundamentais do manejo da palha ocorre no momento da colheita, uma vez que as máquinas empregadas na operação são responsáveis por fazer a distribuição dos restos culturais ao solo. É de fundamental importância que esses restos culturais sejam distribuídos de forma homogênea ao solo, buscando evitar locais de acúmulo de palha, haja visto que ocasionarão problemas de embuchamento na semeadora no momento da semeadura da próxima cultura, além de irregularidades na profundidade das sementes. Já na condição contrária, ausência de palha, o problema está na maior propensão ao desenvolvimento de plantas daninhas (⁴). Diferentes concentrações de palha também podem refletir no aparecimento de faixas de fertilidade diferentes, bem como, podem gerar redução do contato solo/semente, resultando em faixas com menor emergência de plântulas (⁵). Ademais, a maior concentração de palha gerada pela irregularidade de distribuição dos restos culturais realizada pela colhedora, pode resultar em menores percentuais de emergência da cultura subsequente (⁶).

A cultura de aveia preta (Avena strigosa) é muito utilizada no Sul do Brasil durante o inverno, para fins de cobertura de solo, por produzir alta quantidade de palha, ao redor de 4 t ha-1 (⁷). Sua utilização traz inúmeros benefícios para as culturas de verão, quando comparado a prática de pousio durante o período de inverno. Entretanto, a efetividade desta alternativa de cobertura de solo se dá, em grande parte, em virtude do estande de plantas e esse, por sua vez, está aliado a emergência uniforme das plântulas.

Neste sentido, o presente estudo objetivou avaliar o efeito espacial da irregularidade de distribuição de restos culturais da colheita de soja no estande final de aveia preta.

Material e métodos

O experimento foi conduzido via técnicas de agricultura de precisão, geoestatística, realizado em uma área agrícola de 2,6 hectares localizada no município de Sertão/RS no ano agrícola 2018/2019. Situado a 28°0’34”S e 52°25’4”O e altitude de 620 metros, o local possui clima subtropical úmido (cfa), com chuvas bem distribuídas e temperatura média anual de 18,3°C (⁸), o solo é classificado como Nitossolo Vermelho Distrófico (⁹). A área é manejada em SPD há mais de 30 anos, na qual são cultivadas grandes culturas, como soja, trigo e aveia preta.

Para implantação do experimento, um gride de 20 pontos foi determinado na área, sendo 10 pontos relativos ao centro da colhedora e os outros 10 à extremidade (figura 1). Os parâmetros avaliados foram a quantidade de palha e o estande final de plantas emergidas de aveia preta. Para o primeiro parâmetro, foi extraída a massa de palha presente em 1 m2 , esta foi pesada e estimada seu valor por hectare, parâmetro mensurado no dia que foi realizada a semeadura de aveia preta. Após a emergência total da aveia preta, foi determinada o estande de plantas por m2 nos mesmos pontos amostrais anteriores.

Figura 1 Gride da área experimental no município de Sertão, RS, Brasil. Pontos 1-6 e 12-16: centro da colhedora. Pontos 7-11 e 17-21: extremidade da colhedora. Fonte: Google Earth. 

A colheita da soja foi realizada com uma colhedora New Holland® 2017 modelo TC 5070 equipada com plataforma de 20 pés. Para delimitação da área, e georreferenciamento dos pontos, foi empregado um GNSS portátil Garmin® Etrex 20, sendo os dados tabelados em planilha eletrônica.

Os mapas de isolinhas foram construídos por meio de interpolação por krigagem ordinária, levando-se em conta os parâmetros do semivariograma determinados pela análise geoestatística, utilizando o software CR Campeiro 7®.

Afim de verificar diferenças estatísticas dos parâmetros, os dados foram submetidos à análise de variância (ANOVA), teste de normalidade e comparação de médias pelo teste de Tukey à 5% de significância, utilizando o software SISVAR®. A estatística exploratória dos dados (estatística descritiva) foi realizada com auxílio do software PAST®, sendo obtidas medidas de posição (mínimo, média e máximo) e de dispersão (desvio padrão, coeficientes de variação (CV%), assimetria e de curtose). Com base nos valores de assimetria (Cs) obtidos, visando caracterizar quanto e como a distribuição de frequências se afasta da simetria, foram classificados como: Cs = 0 distribuição simétrica, Cs > 0 distribuição assimétrica à direita e Cs < 0 distribuição assimétrica à esquerda. Os valores de curtose (Ck), que expressam a dispersão (achatamento) da distribuição em relação a um padrão, geralmente utilizada a curva normal, foram classificados em: Ck = 0 distribuição mesocúrtica, Ck < 0 distribuição platicúrtica e Ck > 0 distribuição leptocúrtica. Com base nos valores de CV(%) dos dados, a dispersão dos dados foi classificada em: baixa para CV 35% (¹⁰).

Foi realizado teste de correlação linear de Pearson (r) para avaliar o grau de correlação entre as variáveis analisadas. Este coeficiente expressa numericamente, tanto a força quanto o sentido da correlação. Seu valor varia entre 1 e -1 ou seja: -1 ≤ R ≤ 1. Se R = 1, as duas variáveis relacionadas têm correlação perfeitamente positivas; se R = -1, as duas variáveis relacionadas têm correlação perfeitamente negativa; se R = 0, não existe correlação entre as duas variáveis correlacionadas. Quanto mais próximo de 1 (+ ou -), maior a força de correlação (¹¹).

A análise geoestatística foi realizada por meio do ajuste de semivariogramas experimentais, através do software CR Campeiro 7®, determinando efeito pepita (C0 ), contribuição (C1 ), patamar (C0 +C1 ), alcance (a) e índice de dependencia espacial (IDE). O IDE foi calculado por meio das equações: IDE=x(C1 /(C0 +C1 ))*(a/(0,5MD))*100, em que x são constantes (0,375, 0,317 e 0,504 para os modelos esférico, exponencial e gaussiano, respectivamente) e MD é a máxima distancia entre pontos. Se a razão a/0,5MD for maior que 1, utiliza-se 1. Com base no IDE, classificou-se a dependência espacial para o modelo esférico: ≤ 7% (fraco), entre 7% e 15% (moderado) e > 15% (forte); para o modelo exponencial: ≤ 6% (fraco), entre 6% e 13% (moderado) e > 13% (forte); e para o modelo gaussiano: ≤ 9% (fraco), entre 9% e 20% (moderado) e > 20% (forte) (¹²).

Resultados e discussão

A quantidade de palha de soja foi maior na área correspondente ao centro da colhedora, diferindo significativamente em relação a extremidade (tabela 1). No estande final de plantas de aveia preta ocorreu o contrário, foi maior na área das extremidades, diferindo significativamente em relação ao centro, isto demonstra uma relação, em que conforme aumenta a quantidade de palha ocorre a redução do estande final de aveia preta.

Tabela 1 Quantidade de palha de soja (kg ha-1) e estande final de plantas de aveia preta (plantas m-2) em função da posição da colhedora. 

A quantidade de palha distribuída pela colhedora foi espalhada irregularmente, de maneira que se concentra na centro e na medida que se direciona para as extremidades vai reduzindo. Isso demostra uma desuniformidade na distribuição da palha dentro da largura da plataforma de corte, fato também é relatado por (³) (¹³) (¹⁴).

A quantidade média de palha de soja na área (tabela 2) foi de 5618,18 kg ha-1, variando de 1696,97 a 12000,00 kg ha-1, já estande médio foi de 238,04 plantas m-2, variando de 164,71 a 305,88 plantas m-2, essa variação é elevada, fato que confirma que há variação da distribuição pela máquina (tabela 1), sendo visualizado no coeficiente de variação deste, que está classificado como alta dispersão dos dados (heterocedasticidade), já no estande final de plantas de aveia preta há moderada dispersão (¹⁰).

Ambos os conjunto de dados não apresentaram distribuição normal (tabela 2), condição evidenciada nos coeficientes de assimetria e de curtose, sendo que valores distantes de zero podem ser indicativos da não normalidade (¹⁶). A variável quantidade de palha demonstrou coeficiente de assimetria deslocado à direita (Cs>0) enquanto a variável estande final de plantas de aveia preta deslocado a esquerda (Cs<0) (10). O coeficiente de curtose para ambas as variáveis demostrou distribuição platicúrtica a (Ck<0) (¹⁰), o que demostra que as curvas de distribuição dos dados da quantidade de palha e do estande final de plantas de aveia preta são mais achatadas do que a distribuição normal, apresentando variação destes atributos ao longo da área.

Tabela 2 Estatística descritiva da quantidade de palha de soja (kg ha-1) e do estande final de aveia preta (plantas m2 ). 

O efeito pepita demonstra a variabilidade não explicada pelo modelo geoestatístico, levando em consideração a distância de amostragem utilizada, assim a análise do efeito pepita é de grande importância, pois, quanto menor for o seu valor, maior será a continuidade do fenômeno e menor será a variância da estimativa (¹⁷). Os valores encontrados nesse estudo (tabela 3) indicam que o erro foi bastante elevado no estande de plantas, demonstrando a alta variabilidade dos dados dentro da variável.

Tabela 3 Análise geoestatística da quantidade de palha de soja e do estande final de plantas de aveia preta. 

Quanto menor for a proporção do efeito pepita para o patamar do variograma, maior será a continuidade do fenômeno e a variância da estimativa, e maior será a confiança que se pode ter na estimativa (¹⁸), sendo calculada pelo índice de dependencia espacial. A quantidade de palha apresentou dependência espacial classificada como forte, ou seja, tem variação na área, confirmando a discussão demonstrada anteriormente (tabela 1) na comparação entre as posições de coleta. Esse fato comprova que são necessárias alterações construtivas nas colhedoras agrícolas, haja visto que atualmente o que mais é focado é no aumento da capacidade de campo operacional. Avaliando a distribuição da palha de soja em diferentes colhedoras, (¹³) constataram que a distribuição transversal da palha foi desuniforme, independente do mecanismo de trilha e da largura da plataforma, sendo verificado acúmulo de palha no centro da linha de deslocamento e pouco nas extremidades. Na colheita de trigo, a distribuição também foi irregular, conforme relatado nos estudos de (²⁰) (²¹), sendo constatado também, maior acúmulo de palha no centro das máquinas. Isto demonstra que a irregularidade de distribuição de palha das colhedoras é um problema recorrente na agricultura.

No estande final da aveia preta, o IDE foi classificado como dependência fraca, fato que pode estar relacionado à rusticidade desta planta, pois consegue emergir mesmo em condições não tão favoráveis ao desenvolvimento pleno, haja vista, que é possível ver uma relação deste parâmetro com a quantidade de palha (tabela 1), de forma inversa, conforme um aumenta, o outro reduz.

Para a quantidade de palha o modelo teórico que possibilitou os melhores ajustes dos semivariogramas foi o gaussiano, já para o estande final de aveia preta foi o esférico. O modelo gaussiano apresenta um comportamento parabólico próximo da origem e representa fenômenos extremamente contínuos. Já o modelo esférico tem um comportamento linear próximo da origem que representa fenômenos contínuos, mas não diferenciáveis (¹⁵). Tais modelos expressam bem a variação dos dados, em que no primeiro parâmetro observou-se IDE forte, bem como, a sua variação nos dados de CV.

O alcance é um critério importante para estudos de geoestatística, significa a distância máxima que a variável em estudo está correlacionada espacialmente, ou seja, maiores valores de alcance contribuem para uma maior confiança nas estimativas feitas (¹⁹). Para os valores de alcance não houve diferença entre as variáveis analisadas, desta forma, a equidistância da grade amostral foi suficiente para descrever as propriedades avaliadas com exatidão, pois os valores de alcance encontrados para as diferentes propriedades são superiores a da grade. As variáveis analisadas apresentaram correlação forte, negativa e significativa (r = -0,887; p ˂ 0,05) (¹¹), logo, assim como já foi sinalizado nos dados anteriores, quanto maior a quantidade de palha de soja, menor é o estande de aveia preta (figura 2).

Figura 2 Correlação entre o estande final da aveia preta com a quantidade de palha de soja. 

Os mapas de isolinhas da quantidade de palha de soja (figura 3A) e estande final de plantas de aveia preta (figura 3B) demonstram a existência de variabilidade espacial ao longo da área. Os locais em que há maior presença de palha coincidem com aqueles em que há menor estande final de plantas de aveia preta, corroborando para a existência de correlação negativa entre as variáveis (figura 2).

Figura 3 Mapas interpolados por “krigragem” da variabilidade espacial da quantidade de palha de soja (kg ha-1) (A) e do estande final de plantas de aveia preta (plantas m-2) (B). 

A distribuição de palha ficou 51,70% da área concentrada entre 4560 a 7340 kg ha-1, zona mais central da área, local em que a emergência foi 36,98% concentrada no estande final de plantas médio, variando de 200,50 a 241,50 plantas m-2. O maior estande final de aveia preta, classe de 282,75 a 324,00 plantas m-2, ocorreu em apenas 7,92% da área, zona essa que teve a menor classe de quantidade de palha, 1780 a 4650 kg ha-1, demonstrando ação da correlação negativa (Figura 2). A maior parte do resíduo acumulado no centro da linha de deslocamento da colhedora corresponde ao palhiço (material oriundo das peneiras) (²¹), desta forma, (²⁰) destacam a importância de um sistema de distribuição de palhiço, visando a obtenção de uma distribuição de palha mais homogênea, visto que, atualmente, a maioria das colhedoras possuem apenas aletas difusoras para distribuir palha no picador, sem a presença de um mecanismo distribuidor de palha nas peneiras. Na colhedora avaliada, embora possuindo no bandejão divisões para distribuir melhor o volume de palha sobre o sistema de peneira, a oscilação do terreno acaba direcionando a palha de maneira desuniforme, resultando nesse problema.

Nos locais com acúmulo de palha, forma-se uma barreira física que pode atrasar e/ou impedir o surgimento de plântulas, comprometendo o estande de plantas. Conforme (²²), a condição heterotrófica das plântulas (sobrevivência às custas das reservas da semente) permite a elas sobrevir por um curto período de tempo na ausência de luz, no entanto, a com o passar do tempo, a medida em que se tornam autotróficas, é necessário haver a presença de luz para a manutenção das atividades fisiológicas. Portanto, a aveia preta germina, mas não consegue se desenvolver, devido a barreira física gerada pelo excesso de palha de soja que impede a entrada de luz até as plantas, levando-as à morte. Resultados similares foram encontrados por (⁶), que ao estudarem o efeito da palha de soja na emergência de trigo concluíram que o aumento da palha na superfície do solo reduz a emergência das plântulas. A emergência de soja e milho também são afetadas pelo excesso de palha, fato relacionado a maior dificuldade de corte dos mecanismos das semeadoras (²³).

A quantidade de palha não é importante só na emergência da aveia preta, a palha possui efeito no manejo de fertilidade do solo, uma vez que a decomposição libera nutrientes ao meio (²). Este processo, dentro de um sistema de agricultura de precisão, que preconiza pela homogeneidade dos atributos químicos do solo, é fundamental.

A desuniformidade de palha na área (figura 3A) irá repercutir em problemas com plantas daninhas, conforme relatam (²⁴), que a palha assume grande importância no manejo destas, uma vez que ela tem a capacidade de influenciar na dormência, germinação e mortalidade de determinadas espécies daninhas, provocando modificações nas comunidades infestantes. Desta forma, nos locais em que há menor presença de palha de soja, possivelmente haverá maior incidência de plantas daninhas competindo com a cultura da aveia preta.

Embora a manutenção da palha das culturas sobre a superfície do solo seja uma importante prática conservacionista de solo, quando mal manejada pode afetar a emergência das sementes de plantas cultivadas. Neste sentido, a palha precisa garantir a máxima proteção do solo e ao mesmo tempo causar o menor efeito possível sobre a emergência das culturas. A maneira mais prática de gerenciar os resíduos da colheita é com a colhedora. A palha, portanto, deve ser espalhada uniformemente por toda a largura da plataforma de corte.

Conclusões

A variabilidade da distribuição da palha na colheita da soja influencia no estande final da aveia preta, havendo forte correlação entre estes parâmetros e variabilidade espacial de ambos.