Análise da determinação da condutividade hidráulica do solo pelo método do perfil instantâneo

An analysis of soil hydraulic conductivity determination by means of the instantaneous profile method

DIVISÃO 2 - PROCESSOS E PROPRIEDADES DO SOLO

COMISSÃO 2.2 - FÍSICA DO SOLO

Análise da determinação da condutividade hidráulica do solo pelo método do perfil instantâneo¹

An analysis of soil hydraulic conductivity determination by means of the instantaneous profile method

Adriano Dicesar Martins de Araújo Gonçalves^I; Paulo Leonel Libardi^II

^IPós-doutorando em Engenharia Florestal, Departamento de Solos, Universidade Federal de Santa Maria - UFSM. Av. Roraima, 1000, Prédio 42. CEP 97105-900 Santa Maria (RS). E-mail: [e-mail]

^IIProfessor Titular, Departamento de Engenharia de Biossistemas, ESALQ. Av. Pádua Dias, 11. CEP 13418-900 Piracicaba (SP). E-mail: [e-mail]

RESUMO

A condutividade hidráulica do solo K é uma propriedade que expressa a facilidade com que a água nele se movimenta; é importante para o manejo do solo, para a produção das culturas e para a preservação do solo e do ambiente. A determinação da condutividade hidráulica pode ser feita por métodos de laboratório e de campo; o primeiro, com maior controle das condições experimentais e, o segundo, com menor grau de perturbação no solo. Dentre os métodos de campo, o mais utilizado é o do perfil instantâneo, tanto o original como o modificado, em que se obtém teoricamente uma relação entre o conteúdo de água θ e o tempo de redistribuição t, supondo gradiente de potencial total igual a um. No intuito de realizar uma análise mais profunda desse método, foi desenvolvido um estudo com o objetivo de verificar a possibilidade de se obter uma relação do conteúdo de água no solo, em função do tempo de redistribuição sem a suposição da existência de gradiente de potencial total unitário. O estudo foi desenvolvido com dados de quatro solos e a validação da relação proposta foi realizada por comparação da função K(θ) obtida a partir dela com a alcançada tradicionalmente. Com base na análise dos resultados, pode-se concluir que a relação proposta entre θ e t para o método do perfil instantâneo apresentou-se válida e que o gradiente de potencial total pode influir tanto o parâmetro γ como o parâmetro K₀ da equação K = K₀ exp [γ (θ-θ₀)], ajustada aos dados experimentais. No caso específico dos solos avaliados neste trabalho, o gradiente de potencial teve baixa influência no parâmetro γ, comparada à sua influência no parâmetro K₀, levando à possibilidade de uso de um gradiente de potencial total médio para obtenção da função K(θ). Além disso, observou-se que o procedimento para o cálculo da densidade de fluxo a partir da curva da armazenagem em função do tempo de redistribuição da água é muito mais simples do que o procedimento original utilizado no método do perfil instantâneo.

Termo de indexação: armazenagem de água, drenagem interna, conteúdo de água no solo, tempo de redistribuição, gradiente de potencial total.

ABSTRACT

The soil hydraulic conductivity K is a property that expresses the ability of the soil in conducting water. It is important for soil management, crops yield and the preservation of soil and environment. The soil hydraulic conductivity can be determined by means of laboratory and field methods, the first better experimentally controlled and the second with less degree of soil disturbing. Among the field methods, the instantaneous profile one is the most used, both the original and the modified in which a relation between the soil water content θ and the soil water redistribution time t is obtained assuming soil water total potential gradient equal to one. In order to analyse in more detail the method, a study was carried out with the objective of verifying the possibility to obtain the relation of θ as a function of t without the assumption of the existence of total potential gradient equal to one. The study was developed with data of four soils and the validation of the proposed relation was made by comparing the K(θ) function obtained from it and the K(θ) function obtained in the traditional way. From the results, it could be conclude that the proposed relation between θ and t for the instantaneous method showed to be valid and that the total potential gradient may affect both γ and K₀ parameters of the K = K₀ exp [γ(θ - θ₀)] equation fitted to the experimental data. For the specific case of the soils evaluated in this work, the potential gradient had low influence on the γ parameter as compared with its influence on the K₀ parameter, leading to the possibility of utilization of a mean potential gradient to obtain the K(θ) function. Besides, it was observed that the procedure to calculate the flux density from the curve of soil water storage as a function of redistribution time is much simpler than the original that has been used for this purpose in the instantaneous profile method.

Index terms: soil water storage, soil water redistribution time, total potential gradient, soil water content.

INTRODUÇÃO

A condutividade hidráulica do solo é uma propriedade que expressa a facilidade com que a água nele se movimenta, sendo de extrema importância ao uso agrícola e, consequentemente, à produção das culturas e à preservação do solo e do ambiente.

A primeira equação utilizada para quantificar o movimento da água no solo foi introduzida por Henry Darcy, em 1856, o qual trabalhou com colunas de areia saturada com água. Essa equação, conhecida como equação de Darcy, estabelece que a quantidade de água que passa por unidade de tempo e de área pelo meio poroso saturado é proporcional ao gradiente de potencial total da água nesse meio. A constante de proporcionalidade foi denominada por Darcy de condutividade hidráulica, hoje conhecida como condutividade hidráulica do solo saturado K₀.

A condutividade hidráulica do solo saturado descreve a funcionalidade de seu sistema poroso, englobando propriedades relacionadas com a sua porosidade, como quantidade, tamanho, morfologia, continuidade e orientação dos poros. Pelo fato de a condutividade hidráulica na condição saturada depender em grande parte da forma e continuidade do sistema poroso, essa varia fortemente de um local a outro, apresentando valores extremos em determinados locais e podendo diferir também nas distintas orientações do solo. A condutividade hidráulica na condição saturada é mais dependente, portanto, da estrutura do que da textura do solo e, ao elevar o grau de agregação de um solo, o valor dela aumenta. O efeito da estrutura e em especial dos macroporos (diâmetro maior do que 50 ∝m) possibilita a um solo argiloso exibir valores de K₀ similares ao arenoso. Para Libardi (1999), todas essas propriedades do espaço poroso que influenciam a condutividade do solo podem ser reunidas no termo único "geometria porosa dos solos".

Segundo Silva & Kato (1997), a macroporosidade é fator de extrema importância na condutividade hidráulica do solo saturado e sua redução provoca diminuição nos valores de K₀. Esses autores observaram que pequena diferença de macroporosidade entre áreas de manejo convencional e de plantio direto foi acompanhada por grande diferença nos valores de K₀, sendo constatado o valor de 0,0035 cm s^-1, na primeira, e o de 0,0231 cm s^-1, na segunda, evidenciando a importância da macroporosidade na condutividade hidráulica, ou seja, o efeito do manejo na estrutura do solo.

Buckingham (1907), com base na teoria dos fluxos de calor e elétrico, teorizou uma equação para o fluxo da água em meios porosos não saturados, que ficou muito semelhante à equação de Darcy, introduzindo uma relação funcional entre condutividade hidráulica e conteúdo de água no solo e, igualmente, entre o potencial mátrico e conteúdo de água no solo. Tal equação é hoje denominada equação de Darcy-Buckingham.

Richards (1931) combinou a equação de Darcy-Buckinghan com a equação da continuidade, obtendo a equação diferencial geral que descreve o movimento da água no solo, conhecida como equação de Richards. Esse autor foi o primeiro a apresentar um método de determinação da condutividade hidráulica do solo, em condições não saturadas.

Dentre os vários fatores que influenciam a condutividade hidráulica do solo não saturado, destaca-se o conteúdo de água no solo (Libardi & Melo Filho, 2006). A relação entre condutividade e o conteúdo de água no solo apresenta elevada sensibilidade, de tal modo que uma variação de 1 a 2 % no conteúdo de água pode influenciar a condutividade em valores superiores a 170 % (Falleiros et al., 1998).

Marques et al. (2002), relacionando as propriedades hidráulicas com os horizontes pedológicos de dois Latossolos, não constataram a coincidência entre os valores de condutividade hidráulica e os horizontes com homogeneidade textural e estrutural, concluindo que a dinâmica da água nem sempre se altera quando há mudança de horizonte pedológico. Já Carvalho (2002) constatou haver relação entre os horizontes pedológicos e a condutividade hidráulica em um Latossolo.

Os maiores valores de K₀ são encontrados associados aos maiores valores de porosidade. Porém, Bouma (1982) relatou que pequenos poros podem conduzir mais quando são contínuos, enquanto os maiores em uma dada seção podem não contribuir para o fluxo quando apresentam descontinuidade no perfil do solo. Dessa forma, a classificação por tamanho não necessariamente reflete o importante padrão de continuidade dos poros do solo.

Hillel (1971) salientou que a diferença mais importante entre o escoamento em meio saturado e não saturado está na condutividade hidráulica. Quando o solo se encontra saturado, todos os poros estão preenchidos e conduzindo água, sendo sua condutividade máxima nessa condição. Quando o solo se torna não saturado, alguns poros ficam preenchidos por ar e a condutividade hidráulica decresce rapidamente. Tanto para meios porosos saturados como para os não saturados, a maior dificuldade na quantificação do fluxo da água reside na determinação da sua condutividade hidráulica (Pauletto et al., 1988). A condutividade hidráulica envolve o conhecimento da permeabilidade intrínseca, relacionada ao volume total e à distribuição do tamanho de poros, da tortuosidade e da viscosidade dinâmica do fluido (Cintra, 1997).

Reichardt (1988), analisando o conceito de capacidade de campo, observou que o principal fator de redução do fluxo de drenagem é a condutividade hidráulica do solo e não o gradiente de potencial da água. O principal problema ao basearem-se os critérios para definir a capacidade de campo em condutividade hidráulica está na dificuldade de sua estimativa.

A aplicação da teoria do fluxo de água em solo não saturado no campo e em laboratório requer, portanto, o conhecimento da condutividade hidráulica e das características de retenção de água do solo (Klute, 1972).

Silva & Libardi (2000) comentaram que perturbações na superfície do solo podem influenciar nos valores da condutividade hidráulica e proporcionar resultados que não condizem com o solo estudado. Dessa forma, é desejável que os métodos e as técnicas utilizadas para quantificar essas propriedades não provoquem perturbações adicionais no próprio solo e forneçam valores representativos e de precisão adequada.

Bacchi et al. (1991) comentaram que, após o trabalho de Richards (1931), vários métodos foram desenvolvidos para uso em laboratório e campo, para determinar a condutividade hidráulica em solo não saturado e que a maioria deles se baseia em soluções da equação de Richards.

No método do perfil instantâneo para determinação da condutividade hidráulica K, em função do conteúdo de água θ ou função K(θ) (Watson, 1966), procura-se uma solução para a equação de Richards que é a combinação da equação da continuidade com a equação de Darcy-Buckingham, na direção vertical, durante o processo de redistribuição da água, após a saturação de um perfil de solo e prevenção do fluxo de água através de sua superfície. Dentre as vantagens desse método, pode-se citar a obtenção de medidas diretas no campo ou outro local de estudo, apresentando valores mais realísticos de condutividade hidráulica, quando comparados com os de laboratório.

Para perfis de solo heterogêneo, o método do perfil instantâneo pode apresentar dificuldades como a presença de camadas menos permeáveis que podem interferir no fluxo de água, ocasionando fluxo lateral. Se a água subterrânea se situar muito próxima à superfície, essa pode influenciar os resultados.

Uma alternativa que pode ser utilizada no método do perfil instantâneo é considerar o gradiente de potencial total igual a um. Libardi et al. (1980) simplificaram o procedimento para determinar a função K(θ), com medidas apenas do conteúdo de água, presumindo uma relação exponencial entre a condutividade hidráulica e esse conteúdo, conseguindo uma equação que relaciona o conteúdo de água com o tempo (t) de redistribuição de água no perfil, e que, a partir dessa, obtém-se a função K(θ).

A relação entre θ e t alcançada por Libardi et al. (1980) tem sido uma função linear entre θ e lnt, mesmo em situações em que o gradiente de potencial total difere da unidade. Assim, no intuito de realizar uma análise mais profunda do método do perfil instantâneo, pretendeu-se neste estudo verificar a possibilidade de obter uma equação do conteúdo de água no solo em função do tempo de redistribuição de modo similar ao modelo de Libardi et al. (1980), mas sem a suposição de gradiente de potencial igual a um.

MATERIAL E MÉTODOS

Solos e dados utilizados

Para o estudo, foram extraídos dados de potencial mátrico e de conteúdo de água com o tempo para cada profundidade de interesse de quatro trabalhos da literatura, em que foi desenvolvido o método do perfil instantâneo em condições de campo, os quais são descritos a seguir.

Brito (2010) desenvolveu seu estudo em um Latossolo Vermelho-Amarelo com textura média, designado por LVA, e Prevedello (1987) aplicou o método num Nitossolo, designado por Ni (Quadro 1).

Cho et al. (1977) realizaram o estudo deles em um Quartzipsamment (areia de duna), designado por AD, em que 10 % do material constituíram-se de partículas com diâmetro menor que 0,22 mm; 20 %, menor que 0,28; 60 %, menor que 0,33; e 100 %, menor que 0,90 mm. Com base nesses valores, o solo foi classificado como areia no triângulo textural, com a densidade das partículas em torno de 2,7 Mg cm^-3 e a do solo, em torno de 1,5 Mg cm^-3.

Klein (1998) desenvolveu o estudo num Latossolo Roxo ácrico (LRa), sob três diferentes tratamentos: mata-LRa_m, em cultura de sequeiro-LRa_s e irrigada-LRa_i (Quadro 2).

Metodologias com base no método do perfil instantâneo utilizadas neste trabalho

Foi calculada a condutividade hidráulica do solo não saturado por três metodologias fundamentadas nos procedimentos experimentais do perfil instantâneo, quais sejam os modelos de armazenagem da água no solo que consideram ou não o gradiente unitário, além da proposição de um novo modelo.

Metodologia da armazenagem considerando o gradiente unitário

Considera-se que a armazenagem na camada 0-L m de profundidade, h_L, varia linearmente com o logaritmo do tempo de redistribuição, lnt:

Derivando, assim, h_L em relação a t:

Considerando igual a 1 o gradiente de potencial total, com o que dh_L/dt=K, na profundidade L, transforma-se a equação 2 em

Considerando também que o conteúdo de água no solo θ na profundidade L também varia linearmente com lnt,

tem-se que

Substituindo a equação 5 na equação 3:

Chamando e , resulta em:

Com isso, calcula-se a' e b' da relação lnK em função de θ a partir das regressões 1 e 4.

Observa-se que a e A são sempre negativos.

Metodologia da armazenagem sem a consideração de gradiente unitário

Tendo em conta a consideração da dependência linear entre h_L e lnt (equação 1), quando não se considera unitário o gradiente de potencial total, obtém-se:

em que q é a densidade de fluxo de água na profundidade L, função do tempo de redistribuição t.

Dividindo q em cada tempo t pelo gradiente de potencial total no mesmo tempo t, obtém-se K para cada tempo e assumindo-se que o logaritmo neperiano desse K seja função de θ, tem-se que:

Modelo proposto

É comum na determinação da densidade de fluxo da água pela equação de Darcy-Buckingham, calcular-se o gradiente de potencial total em z por:

sendo z = (z₁ + z₂)/2.

Entretanto, assumindo que o potencial total φ_t varia com a coordenada vertical de posição z, segundo uma equação de segundo grau do tipo

a equação 10, tendo em vista a equação 11, pode ser escrita como:

ou:

Como (z₁ + z₂) = 2z, então:

que é a derivada de (Equação 11) com relação a z.

Portanto, quando se calcula o gradiente de potencial total por , está subentendida a consideração que φ_t varia com z conforme a equação 11, visto que a derivada de φ_t com relação a z dessa equação é igual a 2A_pz.

Nas equações 10, 11 e 12, o eixo de z foi escolhido positivo para cima com sua origem coincidente com a superfície do solo.

No desenvolvimento de Libardi et al. (1980), chega-se à equação 13:

em que os autores consideram que:

e . Na equação 14, K₀ e θ₀ são a condutividade hidráulica e o conteúdo de água no solo no tempo zero de redistribuição, respectivamente, e γ o coeficiente angular da equação linear de ln (K) em função de θ.

Na equação 13, L é a profundidade do solo e z a coordenada vertical de posição.

Sem a consideração de gradiente unitário e assumindo que a equação 14 seja válida, tem-se, substituindo as equações 12 e 14 na equação 13, que:

Como L = -z , então:

Separando as variáveis da equação 16 e integrando a equação resultante de θ₀ a θ e de 0 a t, obtém-se:

Na equação 17, t é o tempo de redistribuição da água e a_L é um parâmetro desenvolvido por Libardi et al. (1980) para correlação entre o conteúdo de água θ na profundidade L e o conteúdo de água médio da camada 0-L; com isso, à semelhança desses autores, tem-se por esta proposta, que:

Para tempos longos, as equações 17 e 18 reduzem-se, respectivamente, a:

e

sendo , ou seja, é o valor médio da função A_p(t) no intervalo de 0 a t.

Tendo em vista a equação 4, que independe da suposição de gradiente de potencial total igual ou diferente da unidade, pode-se supor que:

com o que as equações 19 e 20 transformam-se, respectivamente, em:

e

Desse modo, o valor de γ da equação 14 é obtido a partir do coeficiente angular da regressão linear de θ₀-θ, em função de ln(t) (Equação 22), e o valor de K_0, a partir dos coeficientes linear e angular da regressão linear de , em função de ln(t) (Equação 23).

Comparação dos modelos

No intuito de validar o modelo proposto, foi tomado como referência o modelo da armazenagem sem a consideração de gradiente de potencial igual a um. Assim, deduz-se facilmente que, se para o modelo proposto e para o de referência , então:

ou

sendo, evidentemente, e .

Desse modo, se os modelos forem exatamente iguais, evidentemente K₀ = K'₀, γ = γ ' e, portanto,

Consequentemente, se num mesmo gráfico de lnK em função de lnK' forem traçadas as retas representativas das equações 25 e 26, quanto mais próximos entre si ou similares os modelos, mais próximas entre si serão as retas. Para verificar a proximidade dessas retas, calculou-se a área sob elas (integral da função) para os valores de lnK' correspondentes à faixa de conteúdo de água obtida no ensaio de determinação da função K (θ) do solo em estudo.

Para simplificar a notação, seja lnK' = x e lnK = y.

Com isso, a equação 25 transforma-se em:

e a área sob sua curva no intervalo de x_m a x_M é:

sendo x_m o valor mínimo obtido de lnK' e x_M seu valor máximo.

Resolvendo a integral da equação 28 obtém-se:

De modo semelhante, a equação 26 torna-se:

y = x

e a área sob sua curva (1:1) A_r em:

ou

A diferença relativa da integral A em relação à integral A_r:

será, portanto, substituindo as equações 29 e 31 na equação 32, igual a:

RESULTADOS E DISCUSSÃO

A partir dos dados de conteúdo de água à base de volume θ para os diversos tempos de redistribuição t e profundidades dos solos L, foram elaborados os gráficos de θ em cada L e da armazenagem de água até a L de interesse, h_L, em função de t.

Nesses gráficos, as linhas que se ajustaram aos pontos experimentais foram obtidas por meio de regressão linear, cujos coeficientes de determinação (R²) foram maiores do que 0,94 para todas as regressões, evidenciando que de fato descreveram adequadamente o fenômeno, conforme já esclareceram Libardi et al. (1980), os quais trabalharam com a hipótese de gradiente unitário para o método do perfil instantâneo. Na figura 1, é apresentado um exemplo dessas regressões com os dados de armazenagem plotados em função do tempo de redistribuição para a profundidade 0,745 m do Nitossolo.

Carvalho (2002) também obteve boa correlação entre o conteúdo de água médio do solo e o lnt em uma equação de terceiro grau e afirmou que se fizesse um ajuste linear também encontraria alto coeficiente de correlação.

Os resultados de R² das regressões lineares de θ e de h_L em função de lnt indicaram que o denominado método da armazenagem, com ou sem a consideração de gradiente unitário, para avaliar a função K(θ), é alternativa muito mais simplificada em comparação com o procedimento de Hillel et al. (1972), que ainda hoje é utilizado por muitos pesquisadores.

Conforme o desenvolvimento teórico do modelo proposto, além da relação linear entre θ e lnt e entre e lnt (Equações 22 e 23), é necessário também que o logaritmo neperiano de , parâmetro ligado ao gradiente de potencial total (Equação 12), também varie linearmente com lnt (equação 21). As relações em função de lnt apresentaram, para os solos estudados, coeficientes de determinação que variaram entre 0,81 e 0,98, em que se pode inferir que os resultados permitiram utilizar o modelo proposto. A figura 2 também é um exemplo de tais gráficos para o Nitossolo na profundidade 0,745 m. Já o quadro 3 apresenta os valores médios dos gradientes do potencial total, acompanhados dos respectivos desvios-padrão, medidos ao longo do tempo de redistribuição, cujos valores máximos foram 510 h para o Latossolo Vermelho-Amarelo, 720 h para os três tratamentos do Latossolo Roxo, 168 h para a areia de duna e 647 h para o Nitossolo.

Com base no método da estimativa da armazenagem sem a consideração de gradiente unitário, considerado como referência, foram feitas as regressões lineares de lnK em função de θ, cujos coeficientes de determinação (R²) obtidos foram todos acima de 0,97. Esse comportamento exponencial da função K(θ) determinada sob condições de campo pelo método do perfil instantâneo parece ser universal (Hillel et al., 1972; Libardi et al., 1980).

Expressando as equações de ajuste das regressões lineares na forma K =K_oe^γ(θ-θo), obtiveram-se os parâmetros θ₀, γ e K₀ para os três modelos utilizados neste trabalho, cujos resultados são apresentados no quadro 4. Exceção à linearidade ocorreu na areia de duna, em que foram obtidas duas equações lineares de lnK em função de θ ; uma para os tempos iniciais de redistribuição (parte mais úmida) e outra para os tempos finais (parte mais seca) do processo. Evidentemente esse fato foi reflexo de semelhante comportamento das relações entre θ e lnt e de h_L e lnt nesse solo, o que também foi observado por Libardi & Reichardt (2001) para solos de textura média a arenosos.

Pode-se observar, nos solos estudados, que há tendência de a condutividade hidráulica do solo aumentar em profundidade (Quadro 4). Resultados semelhantes foram obtidos por Carvalho (2002), num Latossolo Roxo distrófico, e por Aragão Júnior et al. (1983), num Podzólico Vermelho-Amarelo. Já Cadima et al. (1980), em estudo num Latossolo Vermelho-Amarelo textura média, concluíram que a variabilidade do parâmetro condutividade hidráulica no sentido vertical foi muito pronunciada e salientaram que o uso de valores médios pode trazer erros consideráveis na estimativa de fluxo de água nesse solo. Porém, esses autores acreditaram que a variabilidade em profundidades maiores deve ser menor por causa da maior homogeneidade das camadas mais profundas.

Ainda enfatizando o problema da variabilidade espacial da condutividade hidráulica determinada sob condições de campo, Jong van Lier & Libardi (1999) constataram que as grandes diferenças ocorridas entre as repetições de determinações de K(θ) pelo método do perfil instantâneo devem-se, em menor grau, às diferenças entre γ encontrados e, em maior grau, aos valores de K₀ e θ₀, estimados pelos seus respectivos ajustes.

Calculando os desvios relativos entre os dois métodos (proposto e de referência) pela equação 33, obtiveram-se os valores apresentados na terceira coluna do quadro 5, para os solos estudados. Observa-se, pelos valores de δ, que, em algumas profundidades, o método proposto superestimou os valores de condutividade (valores positivos de δ) e, em outras, subestimou (valores negativos de δ). No entanto, nota-se que os desvios são baixos, com exceção das profundidades 0,745; 0,905 e 1,05 m do Nitossolo, em que os valores de δ foram 12; 34,4 e 9,3 %, respectivamente. Provavelmente, esses maiores valores de δ , notadamente na profundidade de 0,905 m, devem ter relação com o horizonte textural desse solo, cuja expressão máxima ocorreu nessa profundidade, o que pode ter influenciado a condução do ensaio do método no campo.

A quarta coluna do quadro 5 apresenta os valores do desvio δ do método da armazenagem, considerando-se o gradiente de potencial total igual à unidade, em relação ao método de referência, sem consideração de gradiente unitário. Observa-se que os valores de δ são altos, exceto nos casos em que o gradiente de potencial total foi próximo da unidade ao longo do tempo de redistribuição, como na areia de duna e na profundidade de 1,88 m do Nitossolo. A comparação da terceira e quarta colunas do quadro 5 comprova a validade do método proposto.

No Nitossolo, se for considerado x_m = 0 na equação 31, que equivale a um K' = 1 mm dia^-1, valor relativamente baixo, o δ à profundidade de 0,905 m diminuiria de 34,4 a 7,4 %.

Observando-se os dados de gradiente de potencial total, percebe-se que esses não são unitários, mas sua variação ao longo do tempo de redistribuição não foi alta para a maioria das profundidades. Assumindo por hipótese que o gradiente de potencial total seja o mesmo em todos os tempos (valor médio no tempo), verifica-se que basta dividir o valor de K₀ do modelo em que se assume gradiente unitário para obter a nova função K(θ), na forma da equação K =K_oe^γ(θ-θo). Fazendo essa divisão e calculando o desvio γ da nova função (gradiente médio), em relação ao modelo de referência, obtiveram-se os valores também demonstrados no quadro 5, quinta coluna, que foram baixos, à semelhança dos desvios do modelo alternativo aqui apresentado, com pequeno aumento nas profundidades do horizonte textural do Nitossolo. Esse fato evidenciou que o gradiente de potencial total pareceu não ter influência sobre o valor do γ da equação K =K_oe^γ(θ-θo). Entretanto, é importante observar que se considerando a equação 22, essa influência também existe teoricamente no parâmetro γ. Considerando-se os resultados nesses solos, portanto, a função K(θ) pode ser obtida simplesmente dividindo-se o valor de K₀ da equação K =K_oe^γ(θ-θo), obtido com o modelo que assume o gradiente de potencial total igual a um pelo gradiente de potencial total médio, ao longo do tempo de redistribuição da água.

São necessários procedimentos experimentais do método do perfil instantâneo, como tamanho da parcela e isolamento da parcela em profundidade, para se ter certeza de que não há movimento lateral de água durante os processos de saturação do perfil e redistribuição da água nele. Nesse sentido, recomenda-se que no método do perfil instantâneo a coluna de solo no campo seja delimitada em profundidade, com base na análise prévia do perfil no que diz respeito à ocorrência de horizontes ou às camadas que possam desviar o movimento da água de seu trajeto vertical.

CONCLUSÕES

1. A equação teórica que relaciona o conteúdo de água no solo com o tempo de redistribuição no método do perfil instantâneo sem a consideração de gradiente de potencial total unitário apresentou-se válida para essa estimativa. Por essa equação, o gradiente de potencial total pode influir tanto no valor do parâmetro γ como no do parâmetro K₀ da equação K = K₀ exp [γ(θ-θ₀)] da condutividade hidráulica (K) em função do conteúdo de água no solo (θ ). No entanto, para os solos estudados, a influência do gradiente sobre o valor do γ foi irrelevante, comparativamente à influência sobre o valor do K₀.

2. O procedimento para o cálculo da densidade de fluxo, a partir da curva da armazenagem em função do tempo de redistribuição da água no solo pelo método do perfil instantâneo, é muito mais simples que o proposto por Hillel et al. (1972) para esse fim.

LITERATURA CITADA

Recebido para publicação em 11 de maio de 2012

Aprovado em 24 de maio de 2013

• ARAGÃO JÚNIOR, T.C.A.; MAGALHÃES, C.A. & CASTRO, P.T. Determinação da condutividade hidráulica em um solo Podzólico Vermelho-Amarelo em condições de campo. Pesq. Agropec. Bras., 18:805-810, 1983.
• BACCHI, O.O.S.; CORRENTE, J.E. & REICHARDT, K. Análise comparativa de dois métodos simplificados de determinação da condutividade hidráulica de solos. R. Bras. Ci. Solo, 15:249-252, 1991.
• BOUMA, J. Measuring the conductivity of soil horizons with continuous macropores. Soil Sci., Soc. Am. J., 46:438-441, 1982.
• BRITO, A.S. Variabilidade espacial da condutividade hidráulica e da permeabilidade ao ar em função dos conteúdos de água e ar no solo. Piracicaba, Escola Superior de Agricultura "Luiz de Queiroz", Universidade de São Paulo, 2010. 143p. (Tese de Doutorado)
• BUCKINGHAM, E. Studies on the movement of soil moisture. Washington, USDA, 1907. 61p. (Soils Bulletin, 38)
• CADIMA, Z.A.; LIBARDI, P.L. & REICHARDT, K. Variabilidade espacial da condutividade hidráulica em um Latossolo Vermelho-Amarelo textura média, no campo. R. Bras. Ci. Solo, 4:63-66, 1980.
• CARVALHO, L.A. Condutividade hidráulica do solo no campo: simplificações do método do perfil instantâneo. Piracicaba, Escola Superior de Agricultura "Luiz de Queiroz", Universidade de São Paulo, 2002. 89p. (Dissertação de Mestrado)
• CHO, T.; NOMURA, Y.; YANO, T.; SKIKASSHO, S. & INOUE, M. The use of nêutron moisture meter in studies of soil water regimes - water management in a sand dune area. Tottory, Tottori University, Faculty of Agriculture, 1977. 48p
• CINTRA, F.L.D. Disponibilidade de água no solo para porta-enxertos de citros em ecossistema de tabuleiro costeiro. Piracicaba, Escola Superior de Agricultura "Luiz de Queiroz", Universidade de São Paulo, 1997. 90p. (Tese de Doutorado)
• FALLEIROS, M.C.; PORTEZAN, O.; OLIVEIRA, J.C.M.; BACCHI, O.O.S. & REICHARDT, K. Spatial and temporal variability of soil hydraulic conductivity in relation to soil water redistribution, using an exponential model. Soil Till. Res., 45:279-285, 1998.
• HILLEL, D. Soil water: Physical principles and processes. New York, Academic Press, 1971. 288p.
• HILLEL, D.; KRENTOS, V.K. & STILIANOV, Y. Procedure a test an internal drainage method for measuring soil hydraulic characteristics in situ. Soil Sci., 114:395-400, 1972.
• JONG van LIER, Q. & LIBARDI, P.L. Variabilidade dos parâmetros da equação que relaciona a condutividade hidráulica com a umidade do solo no método do perfil instantâneo. R. Bras. Ci. Solo, 23:1005-1014, 1999.
• LIBARDI, P.L. Determinação da condutância hidráulica e da pressão de borbulhamento de cápsulas porosas para tensiômetros. Piracicaba, Escola Superior de Agricultura Luiz de Queiroz, 1999. (Nota de aula)
• LIBARDI, P.L. & REICHARDT, K. Libardi's method refinement for soil hydraulic conductivity measurement. Aust. J. Soil Res., 39:851-860, 2001.
• LIBARDI, P.L. & MELO FILHO, J.F. Análise exploratória e variabilidade dos parâmetros da equação da condutividade hidráulica, em um experimento de perfil instantâneo. R. Bras. Ci. Solo, 30:197-206, 2006.
• LIBARDI, P.L.; REICHARDT, K.; NIELSEN, D.R. & BIGGAR, J.W. Simple field methods for estimating hydraulic conductivity. Soil Sci. Soc. Am. J., 44:3-7, 1980.
• KLEIN, V.A. Propriedades físico-hídrico-mecânicas de um Latossolo Roxo, sob diferentes sistemas de uso. Piracicaba, Escola Superior de Agricultura "Luiz de Queiroz", Universidade de São Paulo, 1998. 150p. (Tese de Doutorado)
• KLUTE, A. The determination of the hydraulic conductivity and diffusivity of unsaturated soils. Soil Sci., 113:264-276, 1972.
• MARQUES, J.D.; LIBARDI, P.L. & JONG van LIER, Q. Relação entre horizontes pedológicos e propriedades hidráulicas em dois Latossolos. R. Bras. Ci. Solo, 26:567-577, 2002.
• PAULETTO, E.A.; LIBARDI, P.L.; MANFRON, P.A. & MORAES, S.O. Determinação da condutividade hidráulica a partir da curva de retenção de água. R. Bras. Ci. Solo, 12:189-195, 1988.
• PREVEDELLO, C.L. Teoria do fluxo da água em solos não saturados: novos conceitos e aplicações Piracicaba, Escola Superior de Agricultura "Luiz de Queiroz", Universidade de São Paulo, 1987. 264 p. (Tese de Doutorado)
• REICHARDT, K. Capacidade de campo. R. Bras. Ci. Solo, 12:211-216, 1988.
• RICHARDS, L.A. Capillary conduction of liquids in porous medium. Physics, 1:381-333, 1931.
• SILVA, C.L. & KATO, E. Efeito do selamento superficial na condutividade hidráulica saturada da superfície de um solo sob cerrado. Pesq. Agropec. Bras., 32:213-220, 1997.
• SILVA, E.B. & LIBARDI, P.L. Propriedades hidráulicas de um solo não saturado utilizando-se o permeâmetro a disco. Sci. Agric., 57:525-529, 2000
• WATSON, K.K. An instantaneous profile method for determining the hydraulic conductivity of unsaturated porous materials. Water Resour. Res., 2:709-715, 1966.

1

Parte da Tese de Doutorado do primeiro autor apresentada ao Programa de Pós-Graduação em Física do Ambiente Agrícola, Escola Superior de Agricultura "Luiz de Queiroz" - ESALQ.

Datas de Publicação

• Publicação nesta coleção
  17 Dez 2013
• Data do Fascículo
  Out 2013

Histórico

• Recebido
  11 Maio 2012
• Aceito
  24 Maio 2013