CONDUTIVIDADE HIDRÁULICA DE UM LATOSSOLO ROXO, NÃO SATURADO, SOB DIFERENTES SISTEMAS DE USO E MANEJO¹

 

UNSATURATED SOIL HYDRAULIC CONDUCTIVITY OF UN OXISOL, UNDER DIFFERENT USE AND MANAGEMENT SYSTEMS

 

Vilson Antonio Klein² Paulo Leonel Libardi³

 

 

RESUMO

O objetivo deste trabalho foi avaliar as alterações na condutividade hidráulica de um Latossolo Roxo ácrico (Typic Acrorthox) do município de Guaira, Região Norte do Estado de São Paulo, submetido ao plantio direto sem irrigação e ao plantio direto irrigado comparando com o solo de mata, condição natural. Determinou-se a condutividade hidráulica do solo (K), no campo, em função do conteúdo de água (q ), pelo método do perfil instantâneo, avaliando-se a redistribuição da água no perfil do solo, por meio de tensiômetros instalados de 0,1 em 0,1m até 1m de profundidade. Para determinar a função K(q ) também em laboratório, amostras de solo com estrutura preservada, de 1,15m de comprimento e 0,2m de diâmetro, foram coletadas nos três manejos. Esses monólitos de solo foram instrumentados com tensiômetros e submetidos a variados fluxos estacionários descendentes de água, para determinação da função, pela equação de Darcy-Buckingham. Os resultados demonstraram que a estrutura do solo foi afetada pelos manejos até a profundidade de 0,4m, com maior intensidade no solo sob plantio direto irrigado, provocando redução da condutividade hidráulica do solo não saturado. O movimento da água no solo não saturado também foi afetado pelo aumento da densidade, especificamente pela alteração na distribuição do diâmetro dos poros, demonstrando ainda que, com a metodologia adotada (campo e monólito) o solo, na condição natural, apresenta duas fases distintas com relação à condução de água, uma próxima à saturação e outra a partir da tensão de 40kPa, pelo que as equações de K(q ) não podem ser extrapoladas para fora da faixa do conteúdo de água em que elas foram obtidas.

Palavras-chave: manejo do solo, dinâmica da água no solo, plantio direto.

 

SUMMARY

The objective of this paper was to assess the changes of the soil hydraulic conductivity of an oxisol (Typic Acorthox) of the county of Guaira (SP), Brazil (20 ^o27’S, 48 ^o18’W, 518m) by comparing the managements of no-tillage-non-irrigated and no-tillage-irrigated with the natural condition (forestry). Soil hydraulic conductivity K as a function of soil water content q was determined in the field by means of the instantaneous profile method, by assessing the water redistribution in the soil profile, by using tensiometers installed at 0.10m soil depth intervals from 0.10m to 1.00m. The K(q ) function was also determined in the laboratory, using undisturbed soil samples of 1.15m length and 0.20m diameter; these samples or soil monoliths from the three treatments, were instrumented with tensiometers and submitted to several steady descendent water fluxes in the laboratory in order to determine the K(q ) functions by means of the Darcy-Buckingham equation. Results have shown that the soil structure was affected by managements till the soil depth of 0.40m, more intensively in the soil under no-tillage-irrigated, leading to a decrease in the unsaturated soil hydraulic conductivity. The unsaturated soil water movement was also affected by the increase of the soil bulk density, due to the change on the soil pore diameter distribution. With the adopted methodology (field and soil monolith), it could be perceived that, under the natural condition, the studied soil presented two distinct phases with relation to soil water conduction, one near the saturation and other from the soil water content at tensions equal and higher than 40kPa; this fact led to the conclusion that the K(q ) equations cannot be extrapolated to soil water content ranges out of those they were measured.

Key words: soil managament, soil water dynamics , no tillage.

 

 

INTRODUÇÃO

O fluxo da água em solo não saturado, na região do sistema radicular dos cultivos é muito pouco estudado. No entanto, o seu conhecimento é de fundamental importância no entendimento dos processos de infiltração, redistribuição e suprimento de água às culturas.

As propriedades físicas do solo, como a textura e a estrutura determinam o fluxo de água no solo. Nesse sentido, ELLIES et al. (1997) destacam a importância da funcionalidade do sistema poroso do solo, englobando propriedades tais como: quantidade, tamanho, morfologia, continuidade e orientação dos poros. Todas estas propriedades do espaço poroso, que influenciam a condutividade, podem ser reunidas no termo único "geometria porosa dos solos" (LIBARDI 1999).

A condutividade hidráulica K do solo não saturado, pode ser descrita como uma função K(q ), onde q = umidade volumétrica do solo, que traduz o quanto ele conduz de água em dada umidade. A contribuição que este fluxo pode representar para a região das raízes das plantas é, ainda, mais importante e foi estudado por JONG VAN LIER& LIBARDI (1997) que observaram que uma raiz individual de uma planta de arroz é capaz de retirar água do solo de uma distância de mais de 0,1m no intervalo de alguns dias. Esses resultados concordam com os de HULUGALLE & WILLATT (1983) que destacam que, em condições de baixa densidade de raízes, deve-se utilizar uma condutividade hidráulica crítica em vez de uma umidade crítica na avaliação do desenvolvimento das plantas.

Segundo SHARMA & UEHARA (1968), a condutividade hidráulica, para Latossolos, decresce rapidamente, quando pequenas tensões são aplicadas. Próximo à saturação, a condutividade foi de aproximadamente 150mm.h^-1, com 10kPa de tensão foi reduzida para 0,1mm.h^-1. OTHMER et al. (1991) atribuem isto à agregação desse solo, que proporciona uma distribuição bimodal do diâmetro dos seus poros os quais são, por isso, classificados em poros inter-agregados e poros intra-agregados. Os primeiros são rapidamente esvaziados fazendo decrescer rapidamente a condutividade, e os segundos, mais lentamente, proporcionando uma condutividade menor.

O comportamento da condutividade hidráulica, em solos com estrutura microagregada próximo a saturação é diferente de outros solos argilosos ( RADULOVICH et al., 1992 e CHEN et al., 1993). Nesses solos, na faixa de umidade de domínio dos macroporos (poros interagregados), a função K(q ) é governada pelo potencial gravitacional, enquanto na faixa de umidade abaixo da capacidade de campo (poros intra-agregados) o potencial mátrico passa a governar, valendo para esta condição a clássica teoria, baseada na equação de Richards.

Estudando as mudanças nas características físicas do solo sob diferentes manejos, BLEVINS et al. (1984) concluíram que a condutividade hidráulica do solo não saturado, em uma faixa de baixas tensões, foi maior para solos sob plantio direto do que em sistema de cultivo convencional. Em contraposição, ROTH et al. (1988) estudando os parâmetros físicos e hídricos do solo em diferentes manejos, constataram que em potenciais abaixo de 20kPa, o plantio direto apresentou menores valores de condutividade hidráulica do que manejo convencional, enquanto, acima deste potencial mátrico, apresentou valores superiores, o que segundo os autores é muito importante no contexto da determinação de água disponível às plantas. Esta aparente contradição pode ser explicada em função de alguns fatores, entre os quais a variação de manejos adotados dentro do plantio direto, uns com maior atividade biológica no solo, outros utilizando culturas com sistema radicular mais agressivo e profundo, e a própria variação das características físicas do solo.

O objetivo deste trabalho foi determinar a função K(q ) de um Latossolo roxo, sob três condições de uso e manejo: plantio direto sem irrigação, plantio direto irrigado e mata natural.

 

MATERIAL E MÉTODOS

A área experimental localiza-se no Município de Guaira, situado na Região Norte do Estado de São Paulo, a qual é muito expressiva em irrigação por aspersão, tendo atualmente em torno de 200 pivôs centrais. O solo da área estudada é o Latossolo Roxo ácrico, A moderado, textura argilosa, profundo, bem drenado em função da sua microagregação e ocupando 85% da superfície do município (TOGNON, 1991; IPT, 1995).

Para avaliar o efeito dos manejos de solo sobre as propriedades físicas e hídricas do solo, três sistemas de uso e manejo foram selecionados, situados em áreas adjacentes. A primeira, com vegetação de mata secundária, não sofreu quaisquer tipos de cultivo e manejo. Suas propriedades físicos e hídricas foram consideradas como referência para comparação com os outros sistemas.

O solo sob condições de sequeiro situa-se em uma área desmatada em 1959, que desde então vem sendo cultivada com duas culturas anuais, uma no período chuvoso e outra na safrinha. Essa área até 1991, foi cultivada no sistema convencional de preparo, com utilização de implementos de discos (arados e grades) e, a partir de então, passou a ser manejada de forma conservacionista, intercalando plantio direto e preparo reduzido.

A terceira área, com solo irrigado, foi desmatada em 1981, quando foi instalado o sistema de irrigação por pivô central, sendo cultivada desde então com culturas anuais em sistema de manejo conservacionista, com plantio direto intercalado com aração profunda, quando da instalação da cultura do tomate (aração profunda em 1990 e 1992).

Para determinar as alterações na estrutura do solo em profundidade, amostras com estrutura preservada foram coletadas, por meio de anéis volumétricos, os quais foram submetidos a diferentes tensões utilizando-se funis de placa porosa e câmaras de Richards. A porosidade foi classificada em função do diâmetro dos poros. Para esta classificação, chamou-se de "macroporos", poros com diâmetro maior que 0,05mm (que perdem a água em tensões menores que 6kPa), "microporos" aqueles com diâmetro entre 0,05 e 0,0002mm (que são esvaziados a tensões entre 6 e 1500kPa) e "criptoporos", poros com diâmetro menor que 0,0002mm (que perdem a água apenas para tensões maiores que 1500kPa).

Para determinação, a campo, da condutividade hidráulica do solo não saturado, nove baterias de 10 tensiômetros (três em cada tratamento) foram instaladas, nas profundidades de 0,1 em 0,1m até 1m, próximo ao local da abertura das trincheiras onde foram coletadas as amostras.

Utilizou-se o método do perfil instantâneo proposto por HILLEL et al. (1972). Para proceder à inundação das parcelas, utilizaram-se folhas de flandre e, unindo as extremidades por meio de parafusos e vedação de borracha, obteve-se um círculo com 5m de diâmetro, correspondendo a 19,60m² de área inundada. Abriu-se um sulco em forma circular, com 0,2m de profundidade, onde foi introduzida a folha de flandre. Após fechado o sulco, iniciou-se a inundação das parcelas, utilizando-se em torno de 10m³ de água. Quando as leituras nos tensiômetros acusavam a saturação do solo até a profundidade de 1,0m, cessava-se a inundação e, para evitar evaporação de água, cobriu-se a parcela com duas lonas de plástico, e com folhas de isopor junto aos tensiômetros, para evitar a variação da temperatura que provocaria a condensação à noite.

No tempo zero de redistribuição, iniciaram-se as leituras dos tensiômetros e após, na primeira hora da manhã, antes de o sol aquecer o plástico e a água nos tensiômetros, efetuavam-se as leituras do potencial mátrico durante um período de aproximadamente 30 dias.

Com os valores do potencial mátrico e das curvas de retenção, para cada camada de solo, gráficos do armazenamento da água no solo em função do tempo para cada profundidade e manejo do solo, foram ajustados aos dados à uma equação do tipo potencial.

onde: h_z = lâmina da água armazenada no solo (mm) até a profundidade z (m); t = tempo(h); a (mm h^-b) e b são os parâmetros empíricos.

Derivando-se a equação 1 com relação ao tempo, obtém-se a equação 2 para cálculo do fluxo (q) de água em cada profundidade e tempo.

Para o cálculo do gradiente de potencial total, obtido através dos tensiômetros, foram utilizados os valores de potencial total dos limites da camada de solo considerada (equação 3), razão pela qual não foi possível calcular a condutividade na primeira e última profundidades.

Calculou-se a condutividade hidráulica do solo não saturado, através da equação:

Para a obtenção de uma relação entre K e q procurou-se ajustar os valores de K(q ) (mm.h^-1) em função de q (m³.m^-3), para cada profundidade e manejo de solo (equação 5).

onde: K = condutividade hidráulica (mm.h^-1); q = conteúdo de água no solo (m³.m^-3); a e b = parâmetros empíricos.

Para determinar a condutividade hidráulica do solo não saturado em laboratório, amostras com estrutura preservada foram coletadas utilizando-se tubos de PVC (polivinil-cloreto) de 1,15m comprimento e 0,2m de diâmetro e com parede de 4,5mm. Utilizando-se uma estrutura de madeira, ancorada ao solo por meio de estacas de ferro e arames, que serviu de apoio a um macaco hidráulico com capacidade de 20kN, introduziu-se o tubo no solo.Conforme o tubo era cravado no solo, cavava-se em torno do mesmo, para facilitar a penetração e a retirada da amostra. Essa coluna foi instrumentada com aplicador de água e tensiômetros. Para determinar a condutividade hidráulica, promoviam-se diferentes fluxos descendentes de água através da coluna, o que era obtido com um sistema de gotejamento na superfície do solo. Esperava-se ocorrer fluxo estacionário, que era medido com uma proveta, em um dispositivo de nível constante, conectado ao aplicador de água da parte inferior da coluna.

 

RESULTADOS E DISCUSSÃO

Determinação a campo

Para o fluxo de água no solo (q), nas diferentes profundidades e tempos, calculado pela equação 2, observa-se uma maior depleção de água na camada superficial para o solo da mata no tempo inicial (até 48 horas) em comparação aos outros dois manejos (Tabela 1). Em termos de manejo da irrigação, isto pode ser considerado um fator negativo, pois nessa condição natural, em função da pouca capacidade de retenção, irrigações mais freqüentes e com menor lâmina deveriam ser aplicadas, apresentando esse solo argiloso comportamento semelhante a um solo arenoso.

Na tabela 2, estão apresentados os parâmetros empíricos do ajuste de K em função de q pela equação 5. Ajustou-se uma única equação aos dados experimentais para as três repetições por profundidade. Observa-se que o ajuste foi satisfatório (p<0,05), e o coeficiente de determinação foi superior a 0,79, indicando não haver para esta propriedade do solo uma variabilidade espacial grande, ao contrário do que constataram REICHARDT et al. (1993) para uma Terra Roxa estruturada em Piracicaba (SP), o que talvez possa ser atribuído a maior variabilidade espacial daquele solo e ao equipamento usado para determinar o conteúdo de água no solo, que foi a sonda de nêutrons.

Analisando-se os parâmetros de ajuste do K em função de q (Tabela 2) pode-se constatar uma menor condutividade para o solo irrigado, até a profundidade de 0,5m, coincidindo com a camada em que ocorreram as maiores alterações na porosidade total e na distribuição do tamanho dos poros em função da compactação do solo (Figura 1), concordando com ROTH et al. (1988) que também encontraram influência negativa sobre o K(q ) em função de manejos e cultivos sucessivos.

É importante destacar que a faixa de umidade (0,45 a 0,30m³.m^-3) na qual se determinou a condutividade hidráulica do solo não saturado no campo é relativamente estreita, devido à rápida drenagem da água no início da determinação e da elevada capacidade de retenção de água que este solo possui. No caso da mata, a maior parte da água drenou rapidamente, dificultando medidas iniciais em função do tempo de resposta dos tensiômetros, uma vez que a condutância hidráulica média das cápsulas porosas era de 1.10^-5cm².s^-1 demorando para atingir o equilíbrio. No sistema irrigado, em função da redução na porosidade do solo acarretada pelo aumento da densidade (Figura 1), a umidade estabilizou em valores bem mais altos drenando menor quantidade de água, permitindo efetuar as determinações em umidades próximas à saturação.

Para o tempo de 12 horas e profundidade de 0,2m em solo de mata, a K(0,391) era de 6,45mm.dia^-1, enquanto no solo de cultivo irrigado, a K(0,425) era de 3,97mm.dia^-1, obtidas a partir da equação 5. Isto evidencia que, no sistema de cultivo irrigado, mesmo mantendo a umidade do solo mais elevada, a condutividade hidráulica do solo não saturado é menor. A constatação de que a condutividade hidráulica do solo não saturado é afetada pelo aumento da densidade em função da diminuição da porosidade e alteração na distribuição dos poros é muito importante e deve ser considerada em termos de manejo da água no solo, conforme também enfatizam HORTON et al. (1994).

Determinação em laboratório

A determinação da K(q ) em amostras com estrutura preservada coletadas no campo se mostrou uma metodologia simples, rápida de ser executada, permitindo avaliações em uma faixa de umidade diferente daquela obtida a campo (0,60 – 0,45m³.m^-3), mais próxima à saturação do solo, condição importante para a boa compreensão do processo de infiltração da água no solo.

Na tabela 3, estão apresentados os parâmetros empíricos do ajuste de ln K versus q , permitindo obter a função K(q ). Destaca-se que só foi possível determinar a K(q ) em cinco profundidades (0,2; 0,3; 0,4; 0,5; 0,6), pois a primeira 0,1m foi utilizada para o cálculo do grad f t na profundidade 0,2m (0,1 - 0,3m), enquanto os tensiômetros instalados nas profundidades 0,8 e 0,9m não saíram da faixa quase saturada, devido ao nível constante de água estabelecido na parte inferior da coluna de solo.

Observou-se também, com essa metodologia, diferenças no comportamento da condutividade para os diferentes manejos e profundidades, com a mesma tendência observada na determinação a campo, tendo o sistema irrigado menores valores até a profundidade 0,5m.

Análise conjunta de K(q ) realizada no campo e em laboratório

Observando-se os parâmetros de ajuste a e b da equação 5 (Tabela 2 e 3) notam-se as diferenças entre os valores obtidos no campo e em laboratório (Figura 2).

Observou-se, neste trabalho, pela diferença nos valores do coeficiente angular das equações de ajuste (Tabela 2 e 3) a diferença no comportamento da condutividade hidráulica próxima à saturação, em relação à outra fase na qual o potencial mátrico governa o movimento. Estes dados corroboram com as afirmações de RADULOVICH et al., (1992) e CHEN et al., (1993) que descrevem este fenômeno. Observa-se que essa diferença no comportamento foi diminuindo, à medida que, em função do aumento da densidade do solo (Figura 1), houve redução da macroporosidade. Na figura 2, observa-se como este efeito é acentuado no solo da mata, onde existe uma uniformidade de macroporos em grande volume em todo o perfil (Figura 1), enquanto no irrigado na camada até 0,4m de profundidade, praticamente não se observa esse efeito.

A descontinuidade da função K(q ) não pode ser atribuída à utilização de diferentes metodologias, pois em várias profundidades observa-se que os valores de K obtidos em laboratório acompanham a tendência daqueles obtidos no campo.

Analisando o valor do conteúdo de água em que ocorre a mudança na inclinação das retas, em torno de 0,45m³.m^-3 para todos os manejos, relacionando com a curva de retenção da água no solo, para obter potencial mátrico equivalente, calculou-se o tamanho aproximado do poro desta descontinuidade. O diâmetro situou-se em torno de 0,1mm concordando com LUXMOORE (1981) que sugere ser este o menor diâmetro dos poros para serem classificados como macroporos, responsáveis portanto pelo fluxo rápido da água no solo.

Isto sugere que a metodologia de determinação da condutividade hidráulica do solo não saturado a campo e em laboratório poderá ser utilizada para definir o limite do diâmetro dos poros no qual a ação da capilaridade (potencial mátrico) começa a interferir no movimento da água no solo.

 

CONCLUSÕES

Alterações na estrutura do solo, com aumento da densidade, redução da porosidade total e alteração na distribuição do diâmetro dos, diminuem a condutividade hidráulica de um Latossolo roxo não saturado.

O procedimento experimental tradicional do método do perfil instantâneo a campo não permite avaliar corretamente a função K(q ) para conteúdos de água do solo próximos à saturação.

A determinação da função K(q ) na faixa de conteúdo de água próximo à saturação no laboratório, usando-se amostras grandes de solo com estrutura preservada, é uma técnica que pode ser utilizada com sucesso para complementar os dados do método do perfil instantâneo a campo.

 

REFERÊNCIAS BIBLIOGRÁFICAS

BLEVINS, R.L., SMITH, M.S., THOMAS, G.W. Changes in soil properties under no-tillage. In: PHILLIPS, R.E., PHILLIPS, S.H. (Ed.) No-tillage agriculture: principles and practices. New York : VNR, 1984. Cap.9. p.190-230.        [ Links ]

CHEN, C., THOMAS, D.M., GREEN, R.E., et al. Two-Domain estimation of hydraulic properties in macropore soils. Soil Science Society of America Journal, v.57, p.680-686, 1993.        [ Links ]

ELLIES, A., GREZ, R., RAMIREZ, C. La conductividad hidráulica en fase saturada como herramienta par el diagnostico de la estructura del suelo. Agro Sur, v.5, p.51-56, 1997.        [ Links ]

HILLEL, D., KRENTOS, V.D., STYLIANOU, Y. Procedure and test of an internal drainage method for measuring soil hydraulic characteristics in situ. Soil Science, v.114, p.395-400, 1972.        [ Links ]

HORTON, R., ANKENY, M.D., ALLMARAS, R.R. Effects of compaction on soil hydraulic properties. In: SOANE, B.D., VAN OUWERKERK, C. (Ed.). Soil compaction in crop production. Amsterdam : Elsevier, 1994. Cap.7, p.141-165.        [ Links ]

HULUGALLE, N.R., WILLATT, S.T. The role of soil resistance in determining water uptake by plant toot systems. Australian Journal of Soil Research, v.21, p.571-574, 1983.        [ Links ]

INSTITUTO DE PESQUISAS TECNOLÓGICAS DO ESTADO DE SÃO PAULO. Áreas de conflito de uso das terras do município de Guaira - SP. São Paulo, 1995. 60p. (IPT. Relatório Técnico, 33 042).        [ Links ]

JONG VAN LIER, Q. de, LIBARDI, P.L. Extraction of soil water by plants: development and validation of a model. Revista Brasileira de Ciência do Solo, v.212, p.535-542, 1997.        [ Links ]

LIBARDI, P.L. Dinâmica da água no solo. 2.ed. Piracicaba : O autor, 1999. 497p.        [ Links ]

LUXMOORE, R.J. Micro-,meso-, and microporosity of soil. Soil Science Society of America Journal, v.45, p.671-672, 1981.        [ Links ]

OTHMER, H., DIEKKRUGER, B., KUTILEK, M. Bimodal porosity and unsaturated hydraulic conductivity. Soil Science, v.52, p.139-150, 1991.        [ Links ]

RADULOVICH, R., SOLLINS, P., BAVEYE, P., et al. Bypass water flow through unsaturated microaggregated tropical soils. Soil Science Society of America Journal, v.56, p.721-726, 1992.        [ Links ]

REICHARDT, K., BACCHI, O.O.S., VILLAGRA, M.M., et al. Hydraulic variability in space and time in a dark red latosol on the tropics. Geoderma, v.60, p.159-168, 1993.        [ Links ]

ROTH, C.H., MEYER, B., FREDE, H.G., et al. Effect of mulch rates and tillage systems on infiltrability and other soil physical properties of an oxisol in Paraná, Brazil. Soil & Tillage Research, v.11, p.81-91, 1988.        [ Links ]

SHARMA, M.L., UEHARA, G. Influence of soil structure on water relations in low Humic Latosols: II. Water movement. Soil Science Society of America Proceedings, v.32, p.770-774, 1968.        [ Links ]

TOGNON, A.A. Propriedades físico-hídricas do Latossolo Roxo da região de Guaira-SP sob diferentes sistemas de cultivo. Piracicaba, 1991. 85p. Dissertação (Mestrado em Solos e Nutrição de Plantas) - Escola Superior de Agricultura "Luiz de Queiroz", Universidade de São Paulo,1991.        [ Links ]

 

 

¹ Extraído da Tese do Doutorado, apresentada pelo primeiro autor à Escola Superior de Agronomia Luis de Querioz (ESALQ), Universidade de São Paulo (USP), Piracicaba, SP.

² Engenheiro, Agrônomo. Doutor, Professor Faculdade de Agronomia e Medicina Veterinária, Universidade de Passo Fundo. CP. 611, 99001-970, Passo Fundo, RS. E-mail: vaklein@upf.tche.br Autor para correspondência.

³ Engenheiro Agrônomo. Professor Titular DCE, ESALQ, USP. Pesquisador do CNPq. E-mail: pllibard@ciagri.usp.br.

Recebido para publicação em 20.11.00. Aprovado em 30.11.02