Introdução

O entendimento do processo de compactação do solo e a busca por alternativas para mitigar seus efeitos têm motivado a realização de estudos em escala global, uma vez que a compactação é considerada uma das principais ameaças à manutenção da qualidade do solo em áreas agrícolas (^{Schjønning et al., 2009}; ^{Keller et al., 2013}). No cenário agrícola brasileiro, áreas manejadas em plantio direto ocupam mais de 30 milhões do total de 48,8 milhões de hectares cultivados com espécies temporárias (^{Federação Brasileira de Plantio Direto na Palha, 2014}). Nessas áreas, há registros de aumento da compactação do solo na camada localizada entre 0,07 e 0,15 m (^{Suzuki et al., 2008}; ^{Drescher et al., 2011}), decorrente de descuidos no manejo, como ausência de diversificação de culturas, excesso de calagem, uso de semeadoras desprovidas de hastes sulcadoras com deposição de adubo na superfície ou na camada superficial do solo, pastejo intensivo e aumento da frequência de tráfego de máquinas agrícolas (^{Denardin et al., 2011}).

Na região Sul do País, mais de 80% da área de produção de grãos é cultivada com soja ( Glycine max L.), milho ( Zea mays L.) e trigo ( Triticum aestivum L.) em plantio direto (^{Cooperativa dos Agricultores de Plantio Direto, 2012}). Nessas áreas, o aumento da compactação do solo tem provocado restrição ao crescimento radicular das plantas e consequente concentração de raízes na camada superficial do solo, o que pode resultar na diminuição da absorção de água e nutrientes, mesmo em curtos períodos de deficit hídrico, e na redução do crescimento da parte aérea e da produtividade das culturas (^{Klein et al., 2009}).

Como estratégia para descompactar ou retardar o aumento da compactação do solo, habitualmente utilizam-se máquinas para a escarificação (^{Jin et al., 2007}). Entretanto, o uso dessa prática deve ser limitado, pois requer elevado consumo energético, investimento em equipamentos e mudanças no sistema de cultivo (^{Chamen, 2015}), em que os custos aumentam proporcionalmente à profundidade da camada compactada, e os benefícios dependem da textura do solo e da profundidade e da frequência da operação. Além disso, os benefícios da escarificação do solo sobre o rendimento de culturas nem sempre se concretizam (^{Gubiani et al., 2013}; ^{Nunes et al., 2014}).

Outra limitação do uso da escarificação como estratégia de recuperação da estrutura do solo é a efemeridade dos benefícios sobre as propriedades do solo, os quais, muitas vezes, duram menos de 1 ano (^{Álvarez et al., 2009}). Contudo, diversos autores relataram que a duração do efeito da escarificação pode se estender por mais de 2 anos (^{Botta et al., 2006}; ^{Drescher et al., 2011}), de acordo com as características e as propriedades do solo, e com as práticas de manejo a serem utilizadas subsequentemente à escarificação.

Além das características do solo e das práticas de manejo, as divergências na duração dos benefícios da escarificação também se relacionam às propriedades do solo. Em trabalhos com variáveis indicadoras do estado estrutural do solo, como densidade e porosidade, observou-se duração inferior a 1 ano (^{Drescher et al., 2012}; ^{Silva et al., 2012}), mas, em avaliações das propriedades físico-hídricas, diretamente relacionadas a processos do solo, como condutividade hidráulica, verificou-se a manutenção dos efeitos da escarificação por aproximadamente 24 meses (^{Vieira & Klein, 2007}).

Assim, a escolha da propriedade do solo que será avaliada é ponto-chave para determinar a duração dos efeitos da escarificação e permitir a compreensão de processos diretamente associados à variável escolhida. Desse modo, os resultados obtidos podem ampliar ou diminuir o período no qual se afirma existir efeito da escarificação.

O objetivo do presente trabalho foi determinar o tempo de duração do efeito da descompactação do solo, por escarificação mecânica, por meio de indicadores físico-hídricos de Latossolo argiloso, manejado em sistema plantio direto.

Material e Métodos

O estudo foi realizado no campo experimental da Embrapa Trigo, localizado no município de Coxilha, RS. O clima, de acordo com a classificação de Köppen, é do tipo Cfa, subtropical úmido, e o solo é classificado como Latossolo Vermelho distrófico típico (^{Santos et al., 2013}), de textura argilosa (Tabela 1), profundo e bem drenado, localizado em relevo ondulado a suavemente ondulado.

Tabela 1: Composição granulométrica e densidade de partículas de Latossolo Vermelho distrófico típico da área experimental. 

O experimento foi realizado de 2009 a 2012, em área de produção de grãos manejada em sistema plantio direto (SPD) há 27 anos, ou seja, desde 1986. Considerou-se como tratamento o tempo transcorrido após a escarificação mecânica em SPD (E_tempo): E0, E6, E12, E18, E24, E30 e E36, em que os números representam os meses transcorridos após a escarificação, que foram comparados ao tratamento testemunha, constituído pelo SPD contínuo há 27 anos (SPD_cont). As escarificações foram realizadas até 0,25 m de profundidade, com escarificador equipado com cinco hastes, espaçadas em 0,30 m, e rolo destorroador, que dispensa a operação de gradagem subsequente. O experimento foi conduzido em delineamento de blocos ao acaso, com quatro repetições, o que totalizou 32 unidades experimentais de 14,4 m² (8,0x1,8 m).

No período compreendido entre as safras agrícolas de verão 2009/2010 e 2012/2013, a área experimental foi cultivada com a sequência de culturas: milho/trigo, soja/centeio ( Secale cereale L.), milho/trigo e soja. As culturas foram implantadas com semeadora equipada com discos duplos defasados, para atuar a 0,07 m de profundidade.

Em dezembro de 2012 e abril de 2013 (início e fim, respectivamente, do ciclo da cultura de soja, safra 2012/2013), a resistência mecânica do solo à penetração (RP) foi medida em campo, com uso de penetrômetro georreferenciado PNT-2000 (DLG Automação Industrial, Sertãozinho, SP), de acordo com a norma ASAE S313.3. A haste com ponta cônica de 30º e com área do cone de 129 mm² foi inserida no solo com velocidade de penetração de, aproximadamente, 2,0 m min^-1, e os dados de RP foram registrados a cada 0,01 m, até 0,40 m de profundidade. As medições da RP foram realizadas na linha de semeadura e a 0,05, 0,10 e 0,20 m à esquerda e à direita, na entrelinha de semeadura. No momento das medições da RP, foram coletadas amostras de solo, para determinação do conteúdo gravimétrico de água, a 0,03, 0,10, 0,15, 0,20, 0,30 e 0,40 m de profundidade, em três pontos de cada bloco experimental.

A taxa de infiltração de água no solo foi medida em março de 2013, com infiltrômetro de anéis concêntricos, com diâmetros de 0,20 e 0,40 m para anel interno e externo, respectivamente, ambos inseridos no solo até a profundidade de 0,15 m, seguindo método descrito em ^{Claessen (1997)}. Considerou-se como taxa de infiltração estável (TIE) a da última leitura, às 2 horas, em razão da pequena variação da taxa de infiltração, nesse tempo.

Amostras de solo com estrutura preservada foram coletadas com cilindros de aço inox de 0,04x0,055 m, em maio de 2013, para determinação da densidade e da porosidade do solo. As amostras foram coletadas na entrelinha de semeadura da soja, safra 2012/2013, em três profundidades, tendo-se utilizado, para a definição das camadas, a avaliação do perfil cultural, como em ^{Tavares Filho et al. (1999)}. De acordo com este método, consideraram-se, como camada superficial, a que tinha estrutura granular solta e apresentava maior concentração de raízes, aproximadamente de 0 a 0,07 m; como segunda camada, aquela que apresentava estrutura maciça, na qual se percebia reduzida concentração de raízes, aproximadamente de 0,07 a 0,15 m; e, como terceira camada, aquela que apresentava estrutura típica de Latossolo, com pequena presença de raízes, aproximadamente de 0,15 a 0,25 m.

Após a coleta, as amostras foram saturadas por capilaridade por 48 horas, pesadas e levadas à coluna de areia (^{Reinert & Reichert, 2006}); em seguida, foram submetidas à tensão de 6,0 kPa por 48 horas, pesadas e, posteriormente, ressaturadas e submetidas à análise de condutividade hidráulica, a qual foi avaliada em permeâmetro de carga constante, de acordo com Libardi (2012). A condutividade hidráulica (Ks) foi calculada pelo rearranjo da equação de Darcy: Ks = VH ∕At(H + h), em que: V (cm³) é o volume de água que se desloca ao longo da amostra H (cm) e passa pela área da seção transversal A (cm²) da amostra, no intervalo de tempo t (horas); e h é a altura da lâmina de água sobre a amostra (cm).

A porosidade total (Pt), a microporosidade, a macroporosidade e a densidade do solo (Ds) foram determinadas pelo método adotado por ^{Claessen (1997)}. A densidade relativa (Dr) foi calculada com a divisão da Ds pela densidade máxima do solo (Ds_máx). A Ds_máx foi estimada por meio da equação de ^{Marcolin & Klein (2011)}: Ds_máx = -0,0092Argila + 2,0138 (r² = 0,92), com argila em percentagem, tendo-se obtido valores de 1,52, 1,51 e 1,50 Mg m^-3, respectivamente, para as camadas de 0‒0,07, 0,07-0,15 e 0,15‒0,25 m.

As variáveis foram submetidas aos testes de Lilliefors para normalidade e de Cochran para homogeneidade de variâncias. Detectou-se distribuição não normal da Ks e da TIE, em que a transformação logarítmica foi capaz de normalizar apenas a TIE. A distribuição da Ks continuou não normal mesmo após a aplicação de outras transformações. Por isso, a Ks foi analisada pelo teste não paramétrico de Kruskal-Wallis, e os postos dos valores das variáveis, após ordenação, foram comparados pelo teste t, a 5% de probabilidade (^{McDonald, 2014}). As demais variáveis foram submetidas à análise de variância, a 5% de probabilidade. Quando o teste F foi significativo, os resultados foram comparados pelo teste de Dunnett, também a 5% de probabilidade. Todas as avaliações foram realizadas no programa estatístico Assistat, versão 7.7 beta (Universidade Federal de Campina Grande, Campina Grande, PB).

Resultados e Discussão

Alterações na Ds, na Dr, na Pt e na macroporosidade decorrentes da escarificação em SPD foram mantidas apenas na primeira safra agrícola - E0 (Tabela 2), na camada de 0,07 a 0,15 m; nas demais camadas avaliadas, não houve diferença entre os períodos de tempo após escarificação. ^{Silva et al. (2012)} também observaram efeitos da escarificação na Ds com tempo de duração inferior a 1 ano, em Latossolo Vermelho distrófico. A rápida reconsolidação da estrutura do solo, como indicam as propriedades Ds, Pt, macroporosidade e microporosidade, pode acarretar necessidade de repetição semestral de descompactação mecânica, o que se contrapõe aos fundamentos do SPD, uma vez que o uso frequente da escarificação pode romper agregados e acelerar a mineralização da matéria orgânica do solo, o que torna o solo mais suscetível à compactação subsequente em maior intensidade (^{Canarache et al., 2000}).

Tabela 2: Densidade do solo (Ds), densidade relativa (Dr), porosidade total (Pt), micro (Mi) e macroporosidade (Ma) de Latossolo Vermelho distrófico típico, manejado em sistema plantio direto contínuo há 27 anos (SPD_cont) e em SPD escarificado mecanicamente há 0 (E0), 6 (E6), 12 (E12), 18 (E18), 24 (E24), 30 (E30) e 36 (E36) meses⁽¹⁾. 

⁽¹⁾Médias seguidas de letras iguais, nas colunas, não diferem pelo teste de Dunnett, a 5% de probabilidade

^nsNão significativo

É importante ressaltar que a duração dos efeitos da escarificação pode ser intensificada ou retardada, dependentemente da intensidade da manifestação de fatores controladores da reconsolidação e da resiliência do solo, como precipitações, ciclos de umedecimento e secagem, e espécies cultivadas, particularmente quando ocorrem imediatamente após uma determinada escarificação. Contudo, esses fatores não podem ser controlados nesse tipo de experimento executado em campo.

A RP avaliada logo após a semeadura da soja (Figura 1) foi menor do que a avaliada no final do ciclo da cultura (Figura 2), embora o conteúdo de água nas camadas superficiais do solo tenha sido maior no final do ciclo (Figura 3). No presente trabalho, não houve variação do conteúdo de água no solo entre os tratamentos, como evidenciado pelos baixos valores de desvio-padrão obtidos para essa variável; portanto, discutiu-se apenas o efeito da densidade. Isso porque, no momento das medições de RP, no início e no fim do ciclo da soja, safra 2012/2013, houve pequena variação da umidade entre as unidades experimentais e as camadas avaliadas, em decorrência de chuvas ocorridas 2 dias antes do momento das medições.

Figura 1: Resistência do solo à penetração medida aos 20 dias após a semeadura de soja ( Glycine max), safra 2012/2013, em Latossolo Vermelho distrófico típico, manejado em sistema plantio direto contínuo há 27 anos (SPD_cont) e em SPD escarificado mecanicamente há 0 (E0), 6 (E6), 12 (E12), 18 (E18), 24 (E24), 30 (E30) e 36 (E36) meses. 

Figura 2: Resistência do solo à penetração medida ao final do ciclo de soja ( Glycine max), safra 2012/2013, em Latossolo Vermelho distrófico típico, manejado em sistema plantio direto contínuo há 27 anos (SPD_cont) e em SPD escarificado mecanicamente há 0 (E0), 6 (E6), 12 (E12), 18 (E18), 24 (E24), 30 (E30) e 36 (E36) meses. 

Figura 3: Conteúdo gravimétrico de água e desvio-padrão em Latossolo Vermelho distrófico típico, manejado em sistema plantio direto, no momento da avaliação da resistência do solo à penetração, no início, em 20/12/2012 (A), e no fim, em 18/4/2013 (B), do ciclo da cultura de soja ( Glycine max). 

Em todos os tratamentos, a RP dos primeiros 0,05 m foi inferior a 1,0 MPa, o que pode ter sido consequência da maior atuação dos ciclos de umedecimento e de secagem, bem como da atividade biológica, associada à recente ação dos discos da semeadora. Em profundidades superiores a 0,07 m, onde não houve atuação dos discos da semeadora, verificou-se incremento da RP, o que confirma a existência de uma camada mais compactada, localizada entre 0,07 e 0,15 m em áreas em plantio direto, o que também foi relatado por ^{Suzuki et al. (2008)} e ^{Drescher et al. (2011)}.

Os valores baixos de RP, obtidos logo após a semeadura da soja (Figura 1), foram observados nas áreas escarificadas há 0, 6, 12, 18 e 24 meses. Os resultados também foram similares na avaliação realizada no final do ciclo de soja (Figura 2), na qual a escarificação manteve o perfil com menor RP até 18 meses (E18), a partir do qual foi constatada RP superior a 2,0 MPa, valor frequentemente considerado crítico para o desenvolvimento radicular das culturas (^{Collares et al., 2006}; ^{Reichert et al., 2009}). Embora a duração de alterações na Ds e na Pt tenha se mantido por apenas uma safra agrícola, a duração de alterações na RP foi percebida por, no mínimo, 18 meses. Por isso, antes de realizar nova escarificação baseada apenas nos valores de Ds e Pt, é de suma importância verificar, também, os valores de RP, para avaliar se a RP indica compactação restritiva ao desenvolvimento radicular, uma vez que a resposta da planta dependerá da magnitude e da duração da resistência oferecida pelo solo (^{Bengough et al., 2006}). Assim, a avaliação da RP serve como parâmetro melhor que a Ds e a Pt para apoiar a decisão de repetir ou não a escarificação mecânica do solo.

Para a condutividade hidráulica do solo saturado, o efeito da escarificação foi mantido até 18 meses (E18), na camada de 0 a 0,07 m, e até 24 meses (E24) na camada de 0,07 a 0,15 m de profundidade (Tabela 3). Isso indica que, para uma mesma propriedade do solo, a duração da escarificação também pode ser influenciada pela condição inicial de compactação. Dessa forma, tanto para Ks quanto para Pt, microporosidade e Ds (Tabela 2), observou-se maior duração da escarificação na camada de 0,07 a 0,15 m, na qual o grau de compactação foi superior a 90% em quase todos os tratamentos, exceto em E0. Esse comportamento também foi evidenciado por ^{Drescher et al. (2011)} e ^{Nunes et al. (2014)}, que constataram manutenção do efeito da escarificação na densidade e na distribuição do tamanho de poros apenas na camada de 0,07 a 0,20 m, respectivamente, em Latossolo e Nitossolo, ambos de textura argilosa.

Tabela 3: Condutividade hidráulica do solo saturado (mm h^-1) em Latossolo Vermelho distrófico típico em diferentes camadas, manejado em sistema plantio direto contínuo há 27 anos (SPD_cont) e em SPD escarificado mecanicamente há 0 (E0), 6 (E6), 12 (E12), 18 (E18), 24 (E24), 30 (E30) e 36 (E36) meses⁽¹⁾. 

⁽¹⁾Médias seguidas de letras iguais, nas colunas, não diferem pelo teste t, a 5% de probabilidade

^nsNão significativo

⁽¹⁾SP, soma dos postos calculada no teste de Kruskal-Wallis.

Assim como a Ks, a taxa de infiltração de água no solo também detectou as modificações na estrutura do solo promovidas pela escarificação. A taxa de infiltração (Figura 4) decresceu a um valor praticamente estável após 2 horas de infiltração (TIE_2horas), o qual foi utilizado para a comparação dos tratamentos. Alterações na TIE_2horas foram percebidas até 24 meses após a escarificação (E24), apesar de terem sido mais expressivas nos primeiros 12 meses (E0 e E6). Portanto, não foram detectadas diferenças na TIE_2horas em períodos superiores há 24 meses (E30 e E36), em relação à testemunha (SPD_cont).

Figura 4: Taxa de infiltração de água em Latossolo Vermelho distrófico típico, manejado em sistema plantio direto contínuo há 27 anos (SPD_cont) e em SPD escarificado mecanicamente há 0 (E0), 6 (E6), 12 (E12), 18 (E18), 24 (E24), 30 (E30) e 36 (E36) meses, ao longo do tempo de medição (A) e ao final de 2 horas (TIE_2horas) (B). Letras iguais, sobre as colunas, não diferem pelo teste de Dunnett, a 5% de probabilidade. 

A duração dos efeitos da escarificação mecânica depende da propriedade do solo avaliada. As propriedades RP, Ks e TIE_2horas foram as mais sensíveis para detectar a duração das alterações na estrutura do solo provocadas pela escarificação. O uso da Ds e da Pt permite detectar alteração na ocupação do volume de solo por sólidos (Ds) ou vazios (Pt), para períodos de até uma safra agrícola após escarificação. Já o uso da Ks e da taxa de infiltração de água no solo permite detectar alteração no transporte de água no solo, para períodos mais prolongados.

Variáveis como Ks e taxa de infiltração estão diretamente relacionadas à geometria do espaço poroso do solo. De acordo com a Lei de Laplace, o raio dos poros altera a curvatura da superfície da água na interface ar-água, e a sucção em capilares é inversamente proporcional ao raio do poro, como estabelece a equação da capilaridade. Quanto ao fluxo de água, a Lei de Poiseuille determina que o volume de água que flui por unidade de tempo é proporcional à quarta potência do raio do tubo. Embora o sistema poroso do solo seja irregular e não possa ser descrito como um feixe de tubos paralelos, as equações da capilaridade e de Poiseuille são úteis para compreender a importância do raio aparente dos poros no transporte de água no solo. Desse modo, pequenas mudanças no raio dos poros por reconsolidação mudam pouco a Ds e a fração de classes de poros, mas podem mudar muito a taxa de infiltração (sucção capilar e fluxo saturado), a Ks e a TIE (fluxo saturado), o que confere maior sensibilidade a essas variáveis em resposta às alterações promovidas pelo manejo do solo.

Além disso, por haver maior ocupação do volume de solo por raízes após escarificação (^{Reichert et al., 2009}; ^{Nunes et al., 2015}), mais bioporos podem permanecer no solo mesmo que haja reconsolidação de sua estrutura, pois a ausência de revolvimento favorece a estabilização de poros. Os bioporos facilitam os fluxos, tanto de ar quanto de água no solo, o que confere maiores taxas de infiltração de água. Portanto, para conhecer a duração da descompactação mecânica do solo, em condições de campo, a avaliação de propriedades relacionadas ao processo de transporte de água mostra-se mais adequada, uma vez que apresenta o efeito da escarificação mecânica por mais tempo.

Conclusões