Modelagem da curva de retenção de água de Latossolos utilizando a Equação Duplo Van Genuchten

Os Latossolos da região do Cerrado brasileiro, em função da mineralogia da fração argila e estrutura granular, apresentam elevado volume de poros grandes e de poros extremamente pequenos, não mostrando significativo volume de poros intermediários, o que faz com que neles a disponibilidade de água para as plantas seja baixa. Este trabalho objetivou caracterizar e modelar o comportamento da retenção de água em Latossolos oxídicos da região do Cerrado, pertencentes a diferentes classes texturais. Foram coletadas amostras do horizonte Bw de 10 Latossolos oxídicos sob vegetação nativa. A retenção de água nos potenciais matriciais de 1, 2, 4, 6, 8 e 10 kPa foi obtida em unidade de sucção nos potenciais de 33, 60, 100, 500 e 1.500 kPa no extrator de Richards, e a água retida sob potenciais de 1.500 a 300.000 kPa foi quantificada utilizando o psicrômetro de termopar WP4-T. O modelo duplo van Genuchten foi proposto para ajustar os dados experimentais de retenção de água por meio de procedimentos de ajuste de modelos não lineares do software R 2.10.1; também foi avaliada a relação entre as estimativas dos parâmetros do modelo, bem como a inclinação do ponto de inflexão com as propriedades texturais dos solos, aplicando o teste de correlação de Pearson. Os resultados mostraram o bom ajuste do modelo e alto poder de predição, sendo observada correlação do conteúdo de argila do solo com os parâmetros da equação (Usat, Upmp, Ures), assim como, com a inclinação do segundo ponto de inflexão (Itex). A textura dos Latossolos influenciou o comportamento das curvas de retenção de água. As curvas de retenção de água dos Latossolos em estudo apresentaram comportamento bimodal no intervalo de potencial matricial estudado.

Desenvolvimento de um modelo fractal para a estimativa da condutividade hidráulica de solos não saturados

Baseado nos conceitos da geometria fractal e nas leis de Laplace e de Poiseuille, foi criado um modelo geral para estimar a condutividade hidráulica de solos não saturados, utilizando a curva de retenção da água no solo, conforme representada por um modelo em potência. Considerando o fato de que este novo modelo da condutividade hidráulica introduz um parâmetro de interpolação ainda desconhecido, e que, por sua vez, depende das propriedades dos solos, a validação do modelo foi realizada, utilizando dois valores-limite fisicamente representativos. Para a aplicação do modelo, os parâmetros de forma da curva de retenção da água no solo foram escolhidos de maneira a se obter o modelo de van Genuchten. Com a finalidade de obter fórmulas algébricas da condutividade hidráulica, foram impostas relações entre seus parâmetros de forma. A comparação dos resultados obtidos com o modelo da condutividade e a curva experimental da condutividade dos dois solos, Latossolo Vermelho-Amarelo e Argissolo Amarelo, permitiu concluir que o modelo proposto é simples em sua utilização e é capaz de predizer satisfatoriamente a condutividade hidráulica dos solos não saturados.

Modelos de ajuste e métodos para a determinação da curva de retenção de água de um Latossolo-vermelho-amarelo

A caracterização da capacidade de retenção de água de um solo é fundamental para a descrição do fluxo de água através dele e para o adequado manejo da irrigação. São apresentadas comparações entre curvas de retenção de água do solo: ajustadas pelos modelos propostos por van Genuchten e por Hutson & Cass; obtidas pelo método do WP4 usando processo de umedecimento e de secagem; obtidas pelo método da centrífuga utilizando amostras deformadas e indeformadas; e obtidas pelo WP4 e centrífuga, usando um processo de secagem e amostra deformada. Amostras deformadas e indeformadas foram coletadas com trados específicos em um Latossolo Vermelho-Amarelo (LVA), textura argilosa. Foram determinadas as propriedades hídricas do solo necessárias à elaboração das curvas de retenção de água obtidas por análise de regressão. O modelo de van Genuchten possibilitou o melhor ajuste nas diversas situações estudadas. Considerando o processo de secagem, o teor de água útil obtido superou em 13 % o resultado do processo de umedecimento e evidenciou reduzido efeito de histerese. A amostra deformada apresentou-se com um teor de água útil superior 61,7 % ao valor obtido para a amostra indeformada. Entre os métodos estudados, verificou-se que aquele que usa o WP4 subestimou os dados obtidos pela centrífuga. Verificou-se que há diferenças entre amostra deformada e indeformada e os métodos utilizados na obtenção da curva de retenção.

Número mínimo de tensões para determinação da curva característica de retenção de água de um latossolo vermelho eutrófico sob sistema de semeadura direta

Para o conhecimento do comportamento físico-hídrico do solo, é fundamental a determinação da curva de retenção de água (CRA). A aquisição de dados para obtenção da CRA envolve processos demorados e custos elevados. A hipótese deste estudo foi a de que é possível determinar a CRA com menor número de tensões, reduzindo seu tempo de obtenção, sem, contudo, comprometer a acurácia dos resultados. Assim, o objetivo deste trabalho foi determinar quais tensões podem ser utilizadas para determinação da CRA de Latossolo Vermelho eutrófico de maneira acurada e em menor tempo. Foram determinados os conteúdos de água retidos em amostras de um Latossolo Vermelho eutrófico sob sistema de semeadura direta (SSD). As tensões utilizadas na determinação dos conteúdos de água foram: 0, 10, 20, 40, 60, 80, 100, 200, 300, 500, 700, 1.000, 3.000, 5.000 e 15.000 hPa. Os dados foram combinados e determinaram-se 450 CRAs ajustadas pelo modelo de van Genuchten. Os parâmetros α, m, n e Ug res, gerados pelo modelo, foram submetidos à análise de variância (teste F) e as médias comparadas pelo teste de Scott-Knott (p = 0,05). A hipótese deste estudo foi confirmada, ou seja, é possível determinar a CRA com menor número de tensões, reduzindo seu tempo de obtenção em até cinco vezes, sem, contudo, comprometer a acurácia dos resultados. A tensão de 15.000 hPa deve estar contida na combinação de tensões a serem utilizadas para determinação da curva característica de retenção de água, quando se emprega o modelo matemático de van Genuchten para ajuste dos dados. As curvas características de retenção de água do solo, definidas com as combinações de tensões (0, 60, 700 e 15.000 hPa), (0, 80, 700 e 15.000 hPa), (0, 100, 1.000 e 15.000 hPa), (0, 20, 60, 100, 700 e 15.000 hPa), (0, 80, 300, 700 e 15.000 hPa) e (0, 100, 300, 1.000 e 15.000 hPa), sendo os dados ajustados ao modelo de van Genuchten, proporcionaram estimativas dos conteúdos de água, retidos no Latossolo Vermelho eutrófico sob sistema de semeadura direta, com a mesma acurácia, quando comparada à curva de referência partindo-se de um conjunto de 15 valores de tensões.

Curva de retenção de água no solo pelo método do papel-filtro

O método mais tradicional para determinação da curva de retenção de água no solo emprega a câmara de pressão de Richards. Por questões práticas, a busca por métodos alternativos ao de Richards, financeiramente mais acessíveis e menos morosos, faz-se interessante. Diante disso, desenvolveu-se o presente trabalho com o objetivo de avaliar o uso do método do papel-filtro. Para isso, realizou-se inicialmente coleta, caracterização física e preparo de amostras indeformadas de um Latossolo Vermelho distroférrico típico e ensaio em câmara de Richards, utilizando pressões de 0, 10, 30, 60, 100, 300, 500, 1000 e 1500 kPa. No ensaio com o papel-filtro, mediu-se o potencial matricial da água em amostras cujas umidades foram previamente estabelecidas, utilizando-se curva de calibração adequada. Os ensaios resultaram em pontos de pressão versus umidade, que foram ajustados pelo modelo de van Genuchten, utilizando o programa RETC. Realizou-se uma análise comparativa de valores de umidade volumétrica estimados pelo modelo ajustado nos ensaios com o papel-filtro com a curva de retenção ajustada obtida pela câmara de Richards. Por meio dessa comparação, verificou-se a aplicabilidade do método do papel-filtro para determinação da curva de retenção de água no solo agrícola utilizado.

Estudo de parâmetros físicos, químicos e hídricos de um Latossolo Amarelo, na região Amazônica

Estudou-se os atributos físicos, hídricos e químicos de um Latossolo sob sistema de agroflorestal na região Amazônica. O solo foi descrito morfologicamente e amostras com estrutura deformada foram coletadas para caracterização física e química. A análise das propriedades hidráulicas foi realizada a partir de amostras indeformadas coletadas a cada camada de 0,10 m até 1,00 m de profundidade com 5 repetições para determinação da condutividade hidráulica do solo saturado (Ko) pelo método do permeâmetro de carga decrescente e determinação das curvas de retenção (fim(teta)), utilizando as tensões de 1, 2 e 4 kPa (funis de placa porosa) e 10, 30, 50, 100, 500 e 1500 kPa (câmaras de pressão com placa porosa). As demais análises foram realizadas até 1,50 m de profundidade. Os valores de Ko nas camadas apresentaram grande homogeneidade estando compreendido entre 2,22 a 3,20 cm dia-1, sendo relacionados com variação da densidade do solo e porosidade. A retenção hídrica demonstra que o solo tem a capacidade de reter elevada quantidade de água, mesmo quando submetido a elevados potenciais, indicando uma baixa disponibilidade para as plantas. As curvas de retenção não apresentaram um bom ajuste pela equação de van Genuchten , sugerindo uma distribuição de poros bimodal.

Complexação de íons de metais por matéria orgânica dissolvida: modelagem e aplicação em sistemas reais

Os íons de metais pesados são conhecidos tanto pela sua importância fisiológica e industrial, bem como pelo risco ambiental e à saúde humana. Para elucidar o comportamento dessas espécies nos corpos hídricos, os quais recebem grande parte da descarga de metais, seja de origem antrópica ou por fontes naturais, é necessário entender as interações que elas apresentam com o meio, principalmente a especiação química. Um dos mais importantes processos pelos quais passam as espécies metálicas em corpos aquáticos naturais é a interação com a matéria orgânica dissolvida (MOD), que pode ser por adsorção, reações de troca iônica ou por complexação. Neste trabalho foram realizados vários experimentos com o objetivo de descrever o comportamento da complexação de três importantes cátions Cu(II), Cd(II) e Pb(II) com a matéria orgânica (ácido húmico comercial), sob condições diversas de força iônica em meio tamponado. Os resultados foram avaliados de acordo com o modelo de van den Berg e Kramer para a complexação de metais. O modelo foi aplicado na determinação da capacidade de complexação dos íons em amostras reais, oriundas de três rios maranhenses que integram a Amazônia legal: Itapecuru, Bacanga e Pericumã. Nas águas dos rios utilizou-se o parâmetro carbono orgânico dissolvido (COD) para expressar a MOD. Os resultados confirmaram forte interação entre a MOD e íons de metais pesados e que o modelo de van den Berg e Kramer é satisfatório para se estimar a constante de complexação (K) e a concentração de sítios de complexação (Lt). Nas amostras simuladas em laboratório a ordem de complexação dos metais foi Cu(II) > Cd(II) > Pb(II) e a capacidade de complexação mostrou ser linear em função da concentração de ácido húmico comercial. Acredita-se que por ter menor raio iônico, o íon Cu(II) possui maior afinidade com os sítios de complexação. Nas amostras retiradas dos corpos aquáticos, observou-se que o rio com maior concentração de COD (rio Bacanga) apresentou maior capacidade de complexação; entretanto, a ordem foi Pb(II) > Cu(II) > Cd(II) provavelmente devido à presença de íons Cu(II) em maior quantidade nas águas dos rios.

Inter-relações entre as propriedades físicas e os coeficientes da curva de retenção de água de um latossolo sob diferentes sistemas de uso

O comportamento físico e hídrico do solo fundamenta as práticas de seu uso e manejo. O objetivo deste trabalho foi avaliar as inter-relações entre as propriedades físicas e os coeficientes da curva de retenção de água de um Latossolo Vermelho distrófico sob diferentes sistemas de uso. Em 2004, foram selecionadas três áreas contíguas utilizadas por mais de vinte anos com mata (nativa), pousio (Brachiaria decumbens) e cultivo (culturas anuais). Foram coletadas 25 amostras de solo com estrutura não deformada na camada de 0-0,15 m de profundidade, que foram utilizadas para determinação das curvas de retenção de água do solo. As curvas de retenção de água foram ajustadas pela equação de van Genuchten, obtendo-se os coeficientes qs, qr, a e n. Também foram mensuradas as variáveis densidade do solo, teor de C orgânico e teor de argila do solo (grupo I). Essas variáveis e os coeficientes das curvas de retenção do solo (grupo II) foram submetidos às análises de correlações canônicas. Verificou-se que a intensificação do uso do solo (mata, pousio e cultivo) resultou em maiores valores de densidade do solo e na redução dos teores de C orgânico e argila. O primeiro par das variáveis canônicas indicou dependência de qs em relação a densidade, e o segundo par distinguiu os sistemas de uso do solo verificado pela dependência de qr em relação aos teores de argila e de C orgânico. As modificações dos teores de C orgânico e da densidade do solo sugerem a degradação da qualidade física e hídrica do solo.

Variabilidade espacial dos parâmetros hidrodinâmicos de duas parcelas agrícolas no estado da Paraíba

O conhecimento das curvas de retenção e da condutividade hidráulica da água do solo é indispensável para a modelagem do transporte de água e de poluentes em sistemas agrícolas. Essas propriedades são afetadas por numerosas fontes de variabilidade, geralmente associadas a fatores espacial e temporal, e relacionadas com o manejo agrícola. O objetivo desse trabalho foi analisar a variabilidade espacial dos parâmetros das curvas de retenção da água (van Genuchten) e da condutividade hidráulica (Brooks & Corey) da água do solo, obtidos com o método "Beerkan", que se baseia em ensaios de infiltração simplificados e nas informações texturais do solo. Esses parâmetros são de forma ou texturais (parâmetros de ajuste que dão forma às curvas) e de normalização ou estruturais (umidade volumétrica saturada, condutividade hidráulica saturada do solo e potencial de entrada de ar). Os ensaios de campo foram conduzidos em uma malha regular de 25 x 25 m num Neossolo Flúvico e num Latossolo Amarelo cultivados com feijão (Vigna unguiculata (L.) Walp.). Para ambos os solos, os valores dos parâmetros de forma, da umidade volumétrica saturada, dos logaritmos da condutividade hidráulica saturada e do potencial de entrada de ar apresentaram dependência espacial. Os valores da umidade volumétrica saturada e dos logaritmos do potencial de entrada de ar foram mais sensíveis à variação a pouca distância no Neossolo Flúvico do que do Latossolo Amarelo, com semivariogramas apresentando estruturas de variação com fraca dependência espacial.

Funil de haines modificado: curvas de retenção de solos próximos à saturação

O movimento da água através da matriz do solo é geralmente modelado pelas equações de Darcy, Darcy-Buckingam e Richards, cujo uso está baseado no conhecimento de algumas propriedades físicas do solo, como, por exemplo, a distribuição de poros e a retenção de água pelo solo. A retenção de água pelo solo é conhecida a partir da determinação de sua curva de retenção (CR) ou curva característica. O primeiro objetivo deste trabalho foi apresentar um aparato simples, desenvolvido pela modificação do funil de placa porosa (funil de Haines), para a investigação e o levantamento de CRs detalhadas, em amostras de solos em condições de umidade próximas à saturação e em amostras com potenciais mátricos (ψm) que vão desde 0 kPa a aproximadamente -12 kPa (120 cm de altura de coluna de água). Foram investigados agregados de um Latossolo Vermelho distrófico (LVd) de uma região do Estado do Paraná, cujos diâmetros médios e densidades variaram, respectivamente, entre 1,6 e 5,7 cm e 1,01 e 1,31 g cm-3, e amostras arenosas reconstituídas com areias de diâmetros médios entre 0,106 mm e 2,000 mm, com dimensões fractais de fracionamento (Df) entre 2,5 e 3,0. O segundo objetivo do trabalho foi inferir a distribuição de poros das amostras investigadas. Isso foi conseguido utilizando-se os parâmetros de ajuste da curva de van Genuchtenaos pontos das CRs obtidas para a determinação da Função Capacidade de Água (FCA). Pela análise dessas FCAs, observou-se que as amostras de agregados de solo apresentaram um sistema poroso de maior complexidade; e que a variação da granulometria do solo arenoso está diretamente relacionada às modificações de suas propriedades de retenção de água. Finalmente, a variação na densidade dos agregados investigados não alterou significativamente o comportamento das curvas de retenção na faixa de tensões estudadas.

Métodos de obtenção da capacidade de suporte de carga de um argissolo cultivado

O controle da capacidade de suporte à carga do solo pode minimizar os efeitos danosos da compactação. O objetivo deste estudo foi avaliar o uso de métodos de determinação da tensão de pré-consolidação num Argissolo Vermelho-Amarelo cultivado e submetido a diferentes condições. Os tratamentos empregados foram: semeadura direta (testemunha), semeadura direta com intensidade de tráfego de 24,79 Mg km ha-1 e 49,59 Mg km ha-1, cultivo mínimo há cinco meses e cultivo mínimo em solo compactado pelo tráfego de uma pá-carregadora (260 kPa de pressão ao solo e 24,67 Mg km ha-1). Os métodos usados para estimação da capacidade de suporte foram: (1) interseção da reta de compressão virgem (RCV) com o eixo x na deformação zero; (2) intercepto da RCV pela regressão com os primeiros dois pontos da curva; (3) intercepto da RCV pela regressão com os três pontos da curva; (4) resistência do solo a 2,5 % da deformação; (5) Casagrande, por meio do ajuste sigmoidal da equação de van Genuchten (1980); e (6) Pacheco & Silva (ABNT, 1990), também por esse meio. Os métodos de 2,5 % da deformação e a interseção da RCV com eixo x estimaram, no geral, valores abaixo da referência ao longo de todas as camadas analisadas. Já o método de regressão teve comportamento dependente da inclinação da RCV; em inclinações elevadas houve superestimação; e em inclinações baixas, ocorreu o contrário.

Índice S para avaliação da qualidade física de solos

O conceito da qualidade física do solo atualmente é bem estabelecido, principalmente pelos índices que avaliam a qualidade do solo. O índice S é um recente indicador da qualidade física do solo relacionado com características importantes do solo. O objetivo deste trabalho foi fazer uma avaliação do método de obtenção do índice S quanto à derivada da função de curva de retenção de água da curva usando ou não a transformação ln(h), bem como avaliar o problema de unidades e da restrição m = 1 - 1/n no cálculo do índice. Concluiu-se que o valor do índice S é superior quando a umidade do solo usada no ajuste da curva de retenção de água proposta por van Genuchten (1980) é expressa na forma volumétrica, quando comparada com a gravimétrica. Para a mesma unidade da umidade do solo, não se observa diferença nos valores do índice S para qualquer unidade de tensão. O valor do índice S difere de acordo com o uso ou não da restrição m = 1 -1/n para estimativa dos parâmetros m e n da curva de retenção proposta por van Genuchten (1980).

Pedotransfer functions to estimate water retention parameters of soils in northeastern Brazil

Pedotransfer functions (PTF) were developed to estimate the parameters (α, n, θr and θs) of the van Genuchten model (1980) to describe soil water retention curves. The data came from various sources, mainly from studies conducted by universities in Northeast Brazil, by the Brazilian Agricultural Research Corporation (Embrapa) and by a corporation for the development of the São Francisco and Parnaíba river basins (Codevasf), totaling 786 retention curves, which were divided into two data sets: 85 % for the development of PTFs, and 15 % for testing and validation, considered independent data. Aside from the development of general PTFs for all soils together, specific PTFs were developed for the soil classes Ultisols, Oxisols, Entisols, and Alfisols by multiple regression techniques, using a stepwise procedure (forward and backward) to select the best predictors. Two types of PTFs were developed: the first included all predictors (soil density, proportions of sand, silt, clay, and organic matter), and the second only the proportions of sand, silt and clay. The evaluation of adequacy of the PTFs was based on the correlation coefficient (R) and Willmott index (d). To evaluate the PTF for the moisture content at specific pressure heads, we used the root mean square error (RMSE). The PTF-predicted retention curve is relatively poor, except for the residual water content. The inclusion of organic matter as a PTF predictor improved the prediction of parameter a of van Genuchten. The performance of soil-class-specific PTFs was not better than of the general PTF. Except for the water content of saturated soil estimated by particle size distribution, the tested models for water content prediction at specific pressure heads proved satisfactory. Predictions of water content at pressure heads more negative than -0.6 m, using a PTF considering particle size distribution, are only slightly lower than those obtained by PTFs including bulk density and organic matter content.

Estimation of water percolation by different methods using TDR

Detailed knowledge on water percolation into the soil in irrigated areas is fundamental for solving problems of drainage, pollution and the recharge of underground aquifers. The aim of this study was to evaluate the percolation estimated by time-domain-reflectometry (TDR) in a drainage lysimeter. We used Darcy's law with K(θ) functions determined by field and laboratory methods and by the change in water storage in the soil profile at 16 points of moisture measurement at different time intervals. A sandy clay soil was saturated and covered with plastic sheet to prevent evaporation and an internal drainage trial in a drainage lysimeter was installed. The relationship between the observed and estimated percolation values was evaluated by linear regression analysis. The results suggest that percolation in the field or laboratory can be estimated based on continuous monitoring with TDR, and at short time intervals, of the variations in soil water storage. The precision and accuracy of this approach are similar to those of the lysimeter and it has advantages over the other evaluated methods, of which the most relevant are the possibility of estimating percolation in short time intervals and exemption from the predetermination of soil hydraulic properties such as water retention and hydraulic conductivity. The estimates obtained by the Darcy-Buckingham equation for percolation levels using function K(θ) predicted by the method of Hillel et al. (1972) provided compatible water percolation estimates with those obtained in the lysimeter at time intervals greater than 1 h. The methods of Libardi et al. (1980), Sisson et al. (1980) and van Genuchten (1980) underestimated water percolation.

Plant-available soil water capacity: estimation methods and implications

The plant-available water capacity of the soil is defined as the water content between field capacity and wilting point, and has wide practical application in planning the land use. In a representative profile of the Cerrado Oxisol, methods for estimating the wilting point were studied and compared, using a WP4-T psychrometer and Richards chamber for undisturbed and disturbed samples. In addition, the field capacity was estimated by the water content at 6, 10, 33 kPa and by the inflection point of the water retention curve, calculated by the van Genuchten and cubic polynomial models. We found that the field capacity moisture determined at the inflection point was higher than by the other methods, and that even at the inflection point the estimates differed, according to the model used. By the WP4-T psychrometer, the water content was significantly lower found the estimate of the permanent wilting point. We concluded that the estimation of the available water holding capacity is markedly influenced by the estimation methods, which has to be taken into consideration because of the practical importance of this parameter.