Parametrização do método do perfil instantâneo para a determinação da condutividade hidráulica do solo em experimentos com evaporação

Parameterization of the instantaneous profile method to determine soil hydraulic conductivity in evaporation experiments

Resumos

A condutividade hidráulica do solo é uma propriedade cuja quantificação é essencial para qualquer estudo que envolva o movimento da água no solo. Os métodos para sua determinação baseados na drenagem interna, como o do perfil instantâneo - MPI, são os mais empregados, restringindo-se à determinação da condutividade hidráulica entre a condição de solo saturado e a umidade na capacidade de campo. Com a finalidade de obter a condutividade hidráulica para menores valores de umidade, propôs-se um método de emprego do MPI com evaporação, apresentando-se o procedimento para a determinação do plano de fluxo zero e estimativa da densidade de fluxo. Demonstrou-se a aplicação do método, obtendo-se resultados de condutividade hidráulica para a camada superficial de um perfil monitorado a partir de um experimento de campo realizado durante um período de 12 dias. A profundidade máxima do plano de fluxo zero estimado foi de aproximadamente 0,50 m. A faixa de umidade obtida variou de 0,25 a 0,21 m³ m-3, valores complementares aos obtidos em experimentos sem evaporação no mesmo local. Os resultados indicaram a viabilidade do método, útil especialmente para se estender à faixa de umidade obtida na camada superficial de tais experimentos.

plano de fluxo zero; drenagem; ascersão

Soil hydraulic conductivity is an essential property and its quantification is essential for any study that involves soil water movement. Methods of determination based on internal drainage, such as the instantaneous profile method (IPM) are among the most commonly used, but their determination is restricted to the soil water range between saturation and field capacity. In order to obtain hydraulic conductivity for lower water contents, we propose an IPM that takes into account evaporation as well as a calculation routines to determine the zero flux depth and to estimate flux densities. The application of the method is demonstrated by obtaining results for hydraulic conductivity of the surface layer of a monitored profile in a 12 -day field experiment. The estimated maximum zero flux depth was about 0.50 m. Soil water content varied from 0.25 to 0.21 m³ m-3, complementing the range obtained in experiments without evaporation at the same site. Results confirm the applicability of the method, which is especially useful to extend the range of soil water contents in the surface layer obtained in this kind of experiment.

zero flux plante; drainage; ascersion

NOTA

Parametrização do método do perfil instantâneo para a determinação da condutividade hidráulica do solo em experimentos com evaporação

Parameterization of the instantaneous profile method to determine soil hydraulic conductivity in evaporation experiments

Ana Lúcia Berretta Hurtado^I; Rogerio Cichota^II; Quirijn de Jong Van Lier^III

^IDoutoranda em Agronomia pelo Programa de Pós-Graduação em Irrigação e Drenagem da Universidade de São Paulo – ESALQ/USP, Caixa Postal 9, CEP 13418-900 Piracicaba (SP). E-mail: [e-mail]

^IIDoutorando em Ciência do Solo, Massey University, Palmerston North, Nova Zelândia. E-mail: [e-mail]

^IIIProfessor do Departamento de Ciências Exatas, ESALQ/USP. E-mail:[e-mail]

RESUMO

A condutividade hidráulica do solo é uma propriedade cuja quantificação é essencial para qualquer estudo que envolva o movimento da água no solo. Os métodos para sua determinação baseados na drenagem interna, como o do perfil instantâneo – MPI, são os mais empregados, restringindo-se à determinação da condutividade hidráulica entre a condição de solo saturado e a umidade na capacidade de campo. Com a finalidade de obter a condutividade hidráulica para menores valores de umidade, propôs-se um método de emprego do MPI com evaporação, apresentando-se o procedimento para a determinação do plano de fluxo zero e estimativa da densidade de fluxo. Demonstrou-se a aplicação do método, obtendo-se resultados de condutividade hidráulica para a camada superficial de um perfil monitorado a partir de um experimento de campo realizado durante um período de 12 dias. A profundidade máxima do plano de fluxo zero estimado foi de aproximadamente 0,50 m. A faixa de umidade obtida variou de 0,25 a 0,21 m³ m^-3, valores complementares aos obtidos em experimentos sem evaporação no mesmo local. Os resultados indicaram a viabilidade do método, útil especialmente para se estender à faixa de umidade obtida na camada superficial de tais experimentos.

Termos de indexação: plano de fluxo zero, drenagem, ascersão.

SUMMARY

Soil hydraulic conductivity is an essential property and its quantification is essential for any study that involves soil water movement. Methods of determination based on internal drainage, such as the instantaneous profile method (IPM) are among the most commonly used, but their determination is restricted to the soil water range between saturation and field capacity. In order to obtain hydraulic conductivity for lower water contents, we propose an IPM that takes into account evaporation as well as a calculation routines to determine the zero flux depth and to estimate flux densities. The application of the method is demonstrated by obtaining results for hydraulic conductivity of the surface layer of a monitored profile in a 12 -day field experiment. The estimated maximum zero flux depth was about 0.50 m. Soil water content varied from 0.25 to 0.21 m³ m^-3, complementing the range obtained in experiments without evaporation at the same site. Results confirm the applicability of the method, which is especially useful to extend the range of soil water contents in the surface layer obtained in this kind of experiment.

Index terms: zero flux plante, drainage, ascersion.

INTRODUÇÃO

As propriedades que determinam a densidade do fluxo de água no solo são a condutividade hidráulica, que mede a habilidade do solo em conduzir água, e as características de retenção de água no solo, que expressam a sua habilidade em armazenar água (Klute & Dirksen, 1986). Assim, o conhecimento da condutividade hidráulica do solo é essencial para qualquer estudo que envolva o movimento da água no solo, seja para estudar a própria dinâmica da água, seja para estudar o transporte de elementos químicos, nutrientes e defensivos agrícolas, bem como seus impactos potenciais ao ambiente.

A condutividade hidráulica K tem seu valor máximo no solo saturado e é altamente dependente da umidade q. O valor de K decresce acentuadamente com a diminuição de q e, por esse motivo, em geral, expressa-se K em função de q(K(q)). Alternativa-mente, pode-se também escrever K em função do potencial matricial y (K(y)).

Além de variar com a umidade, detecta-se alta variabilidade espacial da condutividade hidráulica no solo, tanto em determinações no campo como no laboratório (Warrick & Nielsen, 1980; Jong van Lier & Libardi, 1999). Os métodos de campo empregados com mais sucesso para sua medição são fundamentados em experimentos de drenagem interna, como o método do perfil instantâneo – MPI – (Hillel et al., 1972), que se baseia na análise dos componentes da equação de Richards em perfis de potenciais e fluxos transientes durante a drenagem vertical que ocorre após a saturação do solo por uma chuva ou irrigação (Paige & Hillel, 1993).

O MPI apresenta como vantagem principal a possibilidade de se medir diretamente no campo as propriedades hidráulicas, sem a necessidade de aguardar o estabelecimento de um equilíbrio dinâmico ("steady state"). Vários autores, comparando resultados do MPI com os de outros métodos, inclusive de laboratório, concluíram que o MPI in situ resulta em determinações mais consistentes, além de sua relativa facilidade de realização (Prevedello et al., 1981; Paige & Hillel, 1993; Marion et al., 1994), chegando a considerá-lo o método padrão.

Entretanto, o método apresenta algumas limitações. Baker et al. (1974) e Paige & Hillel (1993) ressaltam que a drenagem interna pode ser influenciada pela presença da água subterrânea quando esta se encontra próxima à superfície; em perfis de solos heterogêneos, não há uma redistribuição de água uniforme, pois a presença de camadas menos permeáveis impede o fluxo vertical de água, favorecendo o fluxo lateral; e solos com inclinação não permitem uma inundação e uma drenagem vertical homogênea. Contudo, Shouse et al. (1992) aplicaram o método do perfil instantâneo em um perfil verticalmente heterogêneo, utilizando o escalonamento das propriedades hidráulicas, e obtiveram bons resultados.

Um dos procedimentos experimentais comuns do MPI inclui a cobertura da parcela experimental com lona plástica, o que foi sugerido primeiramente por Ogata & Richards (1957), evitando, assim, a evaporação ou a reaplicação de água, facilitando a análise dos dados. Posteriormente, Watson (1966) melhorou essa técnica, tornando o método mais preciso. Dessa forma, a única perda de água é por drenagem, e não se forma uma zona de fluxo ascendente no perfil. Alguns roteiros para processar manualmente os dados desses experimentos estão disponíveis em Libardi (2000). Também podem ser encontradas rotinas computacionais, como o software Kteta (Jong van Lier & Libardi, 1999).

Os experimentos com cobertura de lona plástica somente permitem medir a condutividade hidráulica entre a condição saturada e a capacidade de campo, tornando-se, após essa faixa de umidade, o movimento muito lento para ser detectado (Bouma, 1983; Paige & Hillel, 1993), havendo necessidade de se realizar um experimento que permita a evaporação de água do perfil para obter a condutividade hidráulica para valores mais baixos de umidade. Como no experimento com evaporação a área não é coberta, o monitoramento das propriedades hidráulicas pode ser feito até que haja a ocorrência de uma chuva, o que pode limitar o período de coleta dos dados.

Além das vantagens inerentes ao MPI, o procedimento com evaporação permite ampliar a faixa de valores de umidade em camadas próximas à superfície, ou seja, desde condições mais úmidas até condições mais secas, em um período de tempo bem menor que o necessário, quando se depende apenas da drenagem.

Queiroz et al. (1999), num experimento com evaporação, utilizaram o MPI e obtiveram resultados de K para a profundidade de 0,25 m em uma pequena faixa de umidade (0,23 a 0,25 m³ m^-3). Esses autores concluíram ser necessário um tempo maior de observações, para que haja mais evaporação de água no perfil e, conseqüentemente, obter valores de condutividade em umidades menores.

Quando se utiliza o MPI com evaporação, ocorrem dois fluxos um ascendente na parte superior do perfil, que se encontra sob a ação da evaporação, e um fluxo descendente na parte mais profunda, que compõe a drenagem. A profundidade limítrofe entre o fluxo ascendente e o descendente é denominada profundidade do plano de fluxo zero (z_PFZ), na qual não há movimento de água.

Além das dificuldades para determinar simultaneamente o perfil de potenciais e a variação temporal da armazenagem da água no solo, comuns em experimentos de perfil instantâneo e descritos em diversos trabalhos (Baker et al., 1974; Bouma, 1983; Paige & Hillel, 1993; Jong van Lier, 2001), quando se aplica o método com evaporação acrescenta-se a necessidade da determinação de z_PFZ em função do tempo.

Uma vez que não se encontra na literatura um método para a análise de resultados de experimentos de perfil instantâneo com evaporação, o presente trabalho teve por finalidade propô-lo e apresentar um exemplo de sua aplicação.

MÉTODOS

Fundamentos teóricos

O MPI baseia-se na determinação da condutividade hidráulica do solo pela quantificação dos componentes da equação de Richards, em condições transientes:

em que q (m³ m^-3) é a umidade volumétrica, t (h) é o tempo, K(q, y_m) (m h^-1) é a condutividade hidráulica em função da umidade ou do potencial matricial, y_t (m) é o potencial total da água e z (m) é a coordenada vertical de posição.

Considerando um solo no qual há uma zona de fluxo ascendente na parte superior do perfil e outra de fluxo descendente na parte inferior, separadas por um plano de fluxo zero (PFZ), e integrando-se a equação de Richards (Equação 1) com relação à profundidade, entre a profundidade do plano de fluxo zero (z_PFZ, m) e outra profundidade z (m), com as condições de contorno:

obtém-se:

ou, rearranjando

em que h_zPFZ-z (m³ m^-2) é a armazenagem de água no solo entre a profundidade do plano de fluxo zero (z = z_PFZ) e a profundidade z.

Experimento de campo

O experimento de campo consiste normalmente no monitoramento de instrumentos de medição do potencial matricial da água no solo (tensiômetros, tensímetros) e, ou, da umidade do solo (TDR, sonda de nêutrons) em algumas profundidades ao longo do tempo, partindo de uma condição úmida imposta por uma chuva ou irrigação. No caso de haver apenas a medição de y_m ou q, e não ambas, a curva de retenção deve ser determinada para permitir o processamento dos dados.

Estimativa da posição do plano de fluxo zero

No PFZ, por definição, o fluxo é zero e, pela Lei de Darcy, . Diversos experimentos mostraram que uma equação polinomial de segundo grau com parâmetros de ajuste a (m), b e c (m^-1) se ajusta geralmente bem aos dados experimentais de y_t versus z, de forma que:

e, sendo assim:

Igualando a Equação (5) a zero, obtém-se que z_PFZ é dado por:

Como o PFZ pode mudar de posição ao longo do tempo – a tendência comum é a de se aprofundar com o tempo – um valor para z_PFZ será obtido pela Equação (6) para cada tempo de observação experimental.

Para permitir interpolações e, ou, para suavizar o conjunto de dados obtidos, os valores dos dados estimados de z_PFZ, para cada tempo (t) de observação, podem ser ajustados a uma equação de ajuste.

Cálculo da condutividade hidráulica

A densidade de fluxo de água (q_z da Equação 3) é calculada para cada profundidade e tempo de observação.

Considera-se que o instrumento de medição de umidade instalado na profundidade z_i representa uma faixa de profundidades entre um limite superior (z_a) e um limite inferior (z_b), de forma que e , em que z_i-1 e z_i+1 são as profundidades de instalação imediatamente acima e abaixo de z_i. Nessa faixa, para fins de cálculo da variação de armazenagem (¶h da Equação 3), considera-se que a umidade seja igual à umidade observada em z_i. A densidade de fluxo, o gradiente do potencial total e a condutividade hidráulica são estimados para as mesmas profundidades onde se encontram instalados os instrumentos de medição.

Assim, a densidade de fluxo de água na profundidade de instalação dos instrumentos (profundidade média de cada camada) pode ser calculada como expressado na Equação (2), utilizando-se a taxa de variação da umidade pelo tempo, entre a profundidade média da camada e o z_PFZ.

A condutividade hidráulica pode ser obtida por:

em que K é a condutividade hidráulica (m h^-1) e Ñy é o gradiente de potencial total na camada (m m^-1), calculado pela Equação 5. Considerou-se essa condutividade hidráulica representativa para o teor de água e o potencial matricial no tempo médio e na profundidade de observação.

EXEMPLO DE APLICAÇÃO E DISCUSSÃO

Um experimento de evaporação foi realizado entre os dias 9 e 26 de abril de 2002 (401 h), em uma área situada no campus da Escola Superior de Agricultura "Luiz de Queiroz" – ESALQ/USP em Piracicaba (SP) a 22 ° 42 ' 30 " de latitude sul e 47 ° 38 ' 00 " de longitude oeste, em um Latossolo Vermelho-Amarelo (Typic Hapludox). O resumo de alguns dados meteorológicos do período do ensaio, observados no Posto Meteorológico do Departamento de Ciências Exatas da ESALQ/USP, localizado a 1 km da área experimental, encontra-se na figura 1. As medidas foram iniciadas após uma chuva que deixou o solo em condição muito úmida, não se observando, durante o período de monitoramento, precipitação. No local de observação, encontravam-se instalados cinco tensiômetros com manômetro de mercúrio em configuração circular, nas profundidades médias de 0,075, 0,225, 0,375, 0,525 e 0,675 m, com a finalidade de medir o potencial matricial. A umidade foi estimada por meio da equação de van Genuchten, com os parâmetros das curvas de retenção obtidos por meio de amostras retiradas a menos de 1 m do local e em cada profundidade (q_s = 0,446 m³ m^-3, q_r = 0,209 m³ m^-3, a = 0,26 m^-1, n = 1,496 e m = 13,545) (dados publicados por Cichota, 2003).

Foram feitas 12 observações no tempo, cada uma resultando em um conjunto de dados (y_t, z), aos quais foi ajustada a Equação 4. Na figura 2, são apresentados esses conjuntos e a curva ajustada, utilizada com a finalidade de determinar não só a profundidade do plano de fluxo zero para cada tempo, mas também o gradiente de potencial total (Equação 5).

A equação polinomial utilizada apresentou um bom R² (maior que 0,8 em 10 dos 12 casos) para todos os tempos monitorados, mostrando que a equação representa bem o fenômeno estudado.

Observou-se que um PFZ determinado pelo máximo das curvas ajustadas apareceu a partir de aproximadamente dois dias (Figura 2c) do início do monitoramento. Pode-se notar, nos dois primeiros tempos, um valor máximo em 0,225 m, no entanto, esse máximo não foi evidenciado pelas curvas ajustadas. Optou-se pela estimativa do PFZ pelas curvas descritas pela Equação (4), em vez de utilizar os valores individuais, pois as curvas baseiam-se nos cinco pontos observados. Além disso, a discrepância diminuiu com a profundidade de observação.

Neste trabalho, a equação polinomial de segundo grau (Equação 8) com parâmetros de ajuste p (m), q (m h^-1) e r (m h^-2) foi utilizada para ajustar os dados de z_PFZ, para cada tempo (t) de observação:

Foram estimados valores para z_PFZ pela Equação (6) e estes, em função do tempo, foram ajustados à Equação (8), como mostra a figura 3 para o presente exemplo, evidenciando que o PFZ se aprofunda com o tempo e não ultrapassa a profundidade de aproximadamente 0,50 m, semelhantemente ao encontrado por Queiroz et al. (1999).

A partir do conhecimento do z_PFZ para qualquer tempo, o procedimento de cálculo para a obtenção da densidade de fluxo médio entre dois tempos de observação é ilustrado no esquema da figura 4. Tomando como exemplo a primeira camada do perfil de solo monitorado, o fluxo na profundidade de 0,075 m (q_z=0,075) é calculado por "uma das equações dessa figura".

Desse modo, a condutividade hidráulica do solo foi obtida por meio da Equação (7) para os tempos médios aos tempos de observações, encontrando-se resumidos os valores determinados para a primeira camada (0,075 m) do experimento com evaporação (Quadro 1).

O gradiente de potencial total variou desde valores em torno de 0, no início do experimento, a valores superiores a 10, no final, condizente ao esperado em condições de evaporação em equilíbrio dinâmico, onde a diminuição da umidade e da condutividade hidráulica é compensada pelo aumento de Ñy_t. Observou-se que apenas no primeiro tempo (t = 8,9 h) a drenagem predominou sobre a evaporação, fato indicado pelo gradiente de potencial total negativo.

Observando a figura 5, verifica-se que a faixa de umidade variou de 0,251 m³ m^-3, quando o solo encontrava-se úmido, até 0,209 m³ m^-3, no final do experimento. Acredita-se ser essa uma faixa considerável, uma vez que, em experimentos no mesmo local sem evaporação de água no perfil, a faixa observada para a mesma profundidade foi de 0,396, no início da drenagem, até 0,261 m³ m^-3, no final (Berretta-Hurtado, 2004). Dessa forma, a complementação dos dados com a faixa de umidade observada sob evaporação contribui para um aumento significativo da faixa de umidade avaliada.

Os valores de K versus q determinados para os mesmos tempos são representados na figura 5, conjuntamente com uma equação exponencial ajustada.

CONCLUSÕES

1. Equações polinomiais de segundo grau ajustaram-se bem aos conjuntos de dados de potencial total versus profundidade e da profundidade do plano de fluxo zero e o tempo.

2. Os resultados experimentais apresentados, analisados da forma proposta, indicaram a viabilidade do método, útil especialmente para se estender a faixa de umidade obtida em tais experimentos na camada superficial.

LITERATURA CITADA

Recebido para publicação em outubro de 2003 e aprovado em janeiro de 2005.

• BAKER, F.G.; VENEMAN, P.L.M. & BOUMA, J. Limitations of the instantaneous profile method for field measurement of unsaturated hydraulic conductivity. Soil Sci. Soc. Am. Proc., 38:885-888, 1974.
• BERRETTA-HURTADO, A.L. Variabilidade da condutividade hidráulica do solo em valores fixos de umidade e de potencial matricial. Piracicaba, Escola Superior de Agricultura "Luiz de Queiroz", 2004. 94p. (Tese de Doutorado)
• BOUMA, J. Use of soil survey data to select measurement techniques for hydraulic conductivity. Agric. Water Maneg., 6:177-190, 1983.
• CICHOTA, R. Avaliação no campo de um TDR segmentado para determinação da umidade do solo. Piracicaba, Escola Superior de Agricultura "Luiz de Queiroz", 2003. 100p. (Tese de Mestrado)
• HILLEL, D.; KRENTOS, V.D. & STYLIANOU, Y. Procedure and test of an internal drainage method for measuring soil hydraulic characteristics in situ. Soil Sci., 114:395-400, 1972.
• JONG VAN LIER, Q. Análise de sensibilidade de um algoritmo computacional para processamento de dados de experimentos de perfil instantâneo. Piracicaba, Escola Superior de Agricultura "Luiz de Queiroz", 2001. 73p. (Tese de Livre Docência)
• JONG VAN LIER, Q. & LIBARDI, P. 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.
• KLUTE, A. & DIRKSEN, C. Hydraulic conductivity and diffusivity: Laboratory methods. In: KLUTE, A., ed. Methods of soil analysis. Madison, Soil Science Society of America, 1986. p.697-734.
• LIBARDI, P. Dinâmica da água no solo. 2.ed. Piracicaba, P.L. Libardi, 2000. 509p.
• MARION, J.M.; OR, D.; ROLSTON, D.E.; KAVVAS, M.L. & BIGGAR, J.W. Evaluation of methods for determining soil-water retentivity and unsaturated hydraulic conductivity. Soil Sci., 158:1-13, 1994.
• OGATA, G. & RICHARDS, L.A. Water content changes following irrigation of bare-field soil that is protected from evaporation. Soil Sci. Soc. Am. Proc., 21:355-356, 1957.
• PAIGE, G.B. & HILLEL, D. Comparison of three methods for assessing soil hydraulic properties. Soil Sci., 155:175-189, 1993.
• PREVEDELLO, B.M.S.; PREVEDELLO, C.L. & LIBARDI, P.L. Simplificação analítica do método do perfil instantâneo para obtenção da condutividade hidráulica não saturada em condições de campo. R. Bras. Ci. Solo, 5:93-97, 1981.
• QUEIROZ, M.M.F.; RUSTEBERG, B.; GHEYI, H.R. & GUERRA, H.O.C. Estudo da condutividade hidráulica de um solo aluvial em lisímetros. R. Bras. Eng. Agric. Amb., 3:161-166, 1999.
• SHOUSE, P.J.; SISSON, J.B.; ELLSWORTH, T.R. & JOBES, J.A. Estimating in situ unsaturated hydraulic properties of vertically heterogeneous soils. Soil Sci. Soc. Am. Proc., 56:1673-1679, 1992.
• WARRICK, A.W. & NIELSEN, D.R. Spatial variability of soil physical properties in the field. In: HILLEL, D., ed. Applications of soil physics. New York, Academic Press, 1980. p.319-344.
• WATSON, K.K. An instantaneous profile method for determining the hydraulic conductivity of unsaturated porous materials. Water Res. Res., 2:709-715, 1966.

Datas de Publicação

• Publicação nesta coleção
  30 Jun 2005
• Data do Fascículo
  Abr 2005

Histórico

• Aceito
  Jan 2005
• Recebido
  Out 2003