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

Modeling the water retention curve in Oxisols using the Double Van Genuchten Equation

DIVISÃO 2 - PROCESSOS E PROPRIEDADES DO SOLO

2.2 - FÍSICA DO SOLO

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

Modeling the water retention curve in Oxisols using the Double Van Genuchten Equation

Carla Eloize Carducci^I; Geraldo César de Oliveira^II; Eduardo da Costa Severiano^III; Walmes Marques Zeviani^IV

^IDoutoranda em Ciência do Solo, Universidade Federal de Lavras – UFLA. Caixa Postal 37, CEP 37200-000 Lavras (MG). Bolsista do CNPq. E-mail: [e-mail]

^IIProfessor Associado da UFLA. Bolsista CNPq, E-mail: [e-mail]

^IIIProfessor D-I, Instituto Federal de Educação, Ciência e Tecnologia Goiano, Campus Rio Verde (IFGoiano/Campus Rio Verde), CEP 75901-970 Rio Verde (GO). E-mail: [e-mail]

^IVProfessor Assistente I, Departamento de Estatística Campus III - Centro Politécnico, Universidade Federal do Paraná – UFPR. Caixa Postal 19081, CEP 81531-990 Curitiba (PR). E-mail: [e-mail]

RESUMO

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 (U_sat, U_pmp, U_res), assim como, com a inclinação do segundo ponto de inflexão (I_tex). 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.

Termos de indexação: água no solo, Latossolo, modelo matemático.

SUMMARY

The clay fraction mineralogy and the granular structure of Oxisols in the Cerrado (Brazilian savanna) are the reason for the great volume of large pores as well as of extremely small pores, but no significant medium pore volume, resulting in low available water for plants. This study aimed to characterize and model the behavior of water retention in Oxisols with different texture classes of the Cerrado region. Samples were collected from the Bw horizon of 10 Oxisols under native vegetation. Water retention was measured at matric potentials of 1, 2, 4, 6, 8, 10 kPa using a Richards plate, and the water retained in potentials of 1.500 and 300.000 kPa determined by a thermocouple psychrometer model WP4-T Dewpoint Potentia Meter. The double van Genuchten model was proposed to fit the experimental water retention data by non-linear fitting using software R 2.10.1. The relationship between the estimated model parameters and the inflection point with soil texture properties was evaluated also by the Pearson correlation test. Results indicated good fitting of the model and great prediction power and correlation between the soil clay content and equation parameters (U_sat, U_pmp, U_res), as well as with the slope of the second inflection point (I_tex). Oxisol texture influenced the behavior of water retention curves. The water retention curves of Oxisols were bimodal in the studied range of water potential.

Index terms: water content, Oxisol, mathematical model.

INTRODUÇÃO

Os solos presentes na região do Cerrado brasileiro, com destaque aos da ordem Latossolos, são geralmente muito intemperizados e apresentam baixa fertilidade natural, contrastando com a elevada estabilidade dos agregados, em função da mineralogia da fração argila e da estrutura granular. São solos que apresentam elevada porosidade drenável (> 145 µm) e elevado volume de poros, com diâmetro extremamente pequeno (< 2,9 µm) (Oliveira et al., 2004), não apresentando elevado volume de poros intermediários. Os poros pequenos fazem com que grande quantidade de água possa permanecer retida no solo com energia muito alta, o que a torna indisponível às plantas (classificada como água higroscópica), particularmente nos mais argilosos (Giarola et al., 2002).

A curva de retenção de água no solo (CRA) tem sido utilizada para descrever a dinâmica da água no solo (van Genuchten, 1980; Dexter & Bird, 2001). Essa curva representa graficamente a relação entre a energia de retenção de água (potencial matricial, em escala logarítmica) e o conteúdo de água correspondente, o qual é dependente das características intrínsecas de cada solo, resultando da ação conjunta dos atributos do solo, como textura, estrutura, mineralogia e matéria orgânica (Gupta & Larson, 1979; Beutler et al., 2002).

Em termos fitotécnicos, a faixa de água de interesse no solo corresponde ao intervalo compreendido entre a capacidade de campo e o ponto de murcha permanente, convencionalmente denominado de água disponível para as plantas (Oliveira et al., 2004). Por esse motivo, a capacidade de retenção de água no solo tem sido extensamente estudada até o potencial matricial de 1.500 kPa, potencial em que o conteúdo de água no solo é considerado limítrofe, em termos de deficiência hídrica, para o desenvolvimento da maioria das plantas cultivadas (Richards & Weaver, 1943).

Na descrição do comportamento físico-hídrico de solos nessas condições, o modelo de ajuste proposto por van Genuchten (1980) tem sido adotado universalmente e permite relacionar, com alto poder de predição, a energia de retenção e disponibilidade hídrica (Dexter, 2004). Essa equação matemática é caracterizada por apresentar duas assíntotas, relacionadas com os conteúdos de água no solo correspondentes à saturação (U_sat) e ao conteúdo residual (U_res), e um ponto de inflexão entre os platôs, o qual é dependente dos atributos do solo, sendo a sua forma e inclinação reguladas por parâmetros empíricos de ajuste do modelo ("α", "n" e "m").

Apesar da sua extensa utilização em relação aos demais modelos de ajuste disponíveis (Dourado Neto et al., 2001), esse modelo não se ajusta adequadamente a solos com distribuição bimodal de poros, isto é, solos que apresentam duas classes de poros contrastantes, classificados em poros estruturais e texturais (Dexter & Richard, 2009). Em razão disso, têm sido propostas modelagens utilizando equações capazes de identificar esse comportamento, em que essas classes de poros são quantificadas por meio de dois pontos de máximo, obtidos a partir da derivação da curva de retenção de água (Dexter et al., 2008; Alfaro Soto et al., 2008), e, consequentemente, dois pontos de inflexão.

O modelo proposto por Dexter et al. (2008), a partir do ajuste realizado pela equação dupla exponencial, permite identificar a distribuição bimodal de poros em solos temperados no intervalo de potencial referente a umidade de saturação (U_sat) até umidade residual (U_res). Entretanto, Alfaro Soto et al. (2008) identificaram, por meio da previsão de retenção de água no solo, a distribuição bimodal de poros em solos tropicais em um intervalo de potencial matricial superior à determinação-padrão da curva de retenção de água.

Apesar de ser considerada indisponível às plantas (Klein & Libardi, 2002), estudos visando quantificar a água residual retida nos solos tornam-se fundamentais, uma vez que essa pode ser determinante na regulação de processos bioquímicos microbianos no solo (Moreira & Siqueira, 2006), na reindução da tolerância à dessecação de sementes germinadas e plântulas quando submetidas a elevado potencial matricial (Ψ_m > 1.500 kPa) (Rodrigues, 2010), além de poder atuar como agente lubrificante entre agregados, quando o solo é submetido a pressões externas durante as operações mecanizadas (Dias Júnior, 2000). Nesse contexto, Severiano (2010) salienta que o tráfego de máquinas para a realização de tratos culturais na agricultura brasileira, em particular em Latossolos oxídicos na região do Cerrado, ainda é intenso em períodos do ano de menor intensidade pluvial, quando a água no solo deve estar retida sob potenciais matriciais acima de 1.500 kPa. Mesmo assim, esses solos apresentam elevada suscetibilidade à compactação em razão da sua baixa capacidade de suporte de carga.

Na obtenção dos resultados analíticos de retenção de água no solo, o uso de aparelhos de placas porosas é considerado método-padrão (Klute, 1986; Embrapa, 1997). Esses aparelhos foram desenvolvidos para quantificar a retenção de água na faixa disponível para as plantas. Em contrapartida, alguns psicrômetros permitem determinar essa relação em uma ampla faixa de potencial (Scalon et al., 2002; Albuquerque et al., 2005), apesar da menor precisão, quando utilizado na determinação de baixos potenciais matriciais (Ψ_m < 100 kPa), comparado aos extratores convencionais (Decagon Device, 2000). Assim, sugere-se o uso conjunto de ambos os métodos em estudos de retenção de água no solo, quando englobar extensa amplitude de potencial mátrico.

Diante do exposto, o presente trabalho objetivou caracterizar e modelar a retenção de água em Latossolos oxídicos da região do Cerrado, pertencentes a diferentes classes texturais.

MATERIAL E MÉTODOS

O trabalho foi desenvolvido utilizando-se materiais de solos provenientes de 10 Latossolos sob vegetação nativa (Quadro 1), classificados conforme Embrapa (2006) e localizados no município de Rio Verde, GO (17 º 47 ' 01 " S, 50 º 57 ' 55 " W), microrregião do sudoeste goiano. O clima é classificado, de acordo com Köppen (1931), como tropical úmido (Aw), com inverno seco e verão chuvoso, temperatura média anual entre 20 e 25 ºC e média pluvial anual variando de 1.500 a 2.000 mm.

Para a realização deste estudo, foram abertas trincheiras com dimensões de 1 x 2 x 1 m em cada local e coletadas no horizonte B latossólico (Bw) de todos os perfis (0,80 a 1,00 de profundidade), quatro amostras de solo com estrutura preservada em cilindros de alumínio de 6,4 cm de diâmetro x 2,5 cm de altura, utilizando o amostrador de Uhland, totalizando 40 amostras. As amostras foram acondicionadas em filmes plásticos e parafinadas, com o objetivo de preservá-las até serem realizadas as análises físico-hídricas.

Em laboratório, as amostras com estrutura preservada foram preparadas de modo que o solo correspondente às porções excedentes das partes superior e inferior do volume de seus respectivos anéis fosse retirado, secas ao ar e passadas em peneiras de 2 mm, obtendo assim as amostras deformadas.

Realizou-se a análise granulométrica via agitação rápida (10.000 rpm) utilizando como dispersante físico e NaOH 1 mol L^-1 em contato com a amostra durante 24 h, como dispersante químico (Embrapa, 1997). Procedeu-se à quantificação das partículas de areia (2,00-0,05 mm), por tamisagem; de argila (< 0,002 mm), pelo método da pipeta; e de silte (0,05-0,002 mm), determinada pela diferença entre o teor total da amostra de solo e os teores de areia e argila.

Os conteúdos de SiO₂, Al₂O₃ e Fe₂O₃ foram determinados do extrato sulfúrico e utilizados no cálculo das relações moleculares Ki (SiO₂/Al₂O₃) e Kr [SiO₂/(Al₂O₃ + Fe₂O₃)] (Embrapa, 2006) para classificação do grau de intemperismo.

Para determinação das curvas de retenção de água no solo (CRA), as amostras com estrutura preservada foram inicialmente saturadas e submetidas às tensões matriciais de 1, 2, 4, 6, 8 e 10 kPa, utilizando as unidades de sucção funis de Buchner e de 33, 60, 100, 500, 1.500 kPa no extrator de Richards (Embrapa, 1997). Esses sistemas são dotados de placas cerâmicas porosas; o extrator é um equipamento de produção e controle de pressão, em que a extração de água ocorre por diferença de potencial entre a placa e a amostra de solo.

Com o intuito de obter os valores de conteúdo de água no solo retido em potenciais matriciais superiores aos determinados pelos extratores de placa porosa (Ψ_m > 1.500 kPa), utilizou-se o psicrômetro de termopar do modelo WP4-T Dewpoint Potential Meter (Decagon Device, 2000). Para isso, foram utilizadas amostras de solo deformadas, totalizando quatro repetições por solo. A amostra de solo foi inicialmente equilibrada ao conteúdo de água correspondente ao potencial de 1.500 kPa, utilizando-se então, a partir desta, subamostras de solo úmido contendo aproximadamente 2 g, que foram colocadas em recipiente próprio do aparelho e encaminhadas para dentro da câmara hermética de leitura. A determinação do potencial matricial foi obtida pelo método do ponto de orvalho, ou seja, no equilíbrio do conteúdo de água no interior da câmara com o potencial da água no solo (Scalon et al., 2002; Decagon Device, 2000).

Após a determinação de cada potencial matricial, as amostras foram sucessivamente submetidas ao secamento natural e, ou, artificial, utilizando temperatura ambiente e, ou, estufa de ventilação forçada, com temperatura controlada entre 55 e 60 ºC, gerando os pontos da curva de retenção de água no solo. Para obter melhor caracterização do segmento da curva referente aos elevados potenciais matriciais, tomou-se de 19 a 28 leituras entre os potenciais de 1.500 e 300.000 kPa (amplitude de leitura permitida pelo aparelho).

Os dados obtidos foram ajustados ao modelo sugerido com base no modelo original proposto por van Genuchten (1980):

com restrição de m=1-1/n (Mualem, 1976), considerados para ambos segmentos da curva, estrutural (m_est) e textural (m_tex). A umidade gravimétrica e o potencial matricial são representados por U e Ψ, respectivamente. Os parâmetros U_res, U_pmp e U_sat representam o platô assintótico inferior (Ψ → ∞) ou conteúdo de água residual assintótica; o platô intermediário, ou seja, valor do conteúdo de água que se observa ligeiramente constante ao redor do ponto de murcha permanente e o platô assintótico superior (Ψ → 0), que indica o conteúdo de água na saturação, respectivamente. Os parâmetros a e n estão associados à escala e forma da curva entre as assíntotas superior, intermediária e inferior; α_est e n_est (estrutural) correspondem ao primeiro segmento e α_tex e n_tex (textural) ao segundo segmento da curva. Esse modelo será denominado de duplo van Genuchten, o qual apresenta sete parâmetros a serem estimados e dois pontos de inflexão (Figura 1).

O modelo duplo van Genuchten proposto foi ajustado aos dados experimentais de retenção de água para cada local, nos quais foi utilizada a umidade gravimétrica média para cada potencial. Isso foi obtido por meio de procedimentos de ajuste de modelos não lineares do software R 2.10.1 (R Development Core Team, 2009; Ritz & Streibig, 2008). Considerou-se o algoritmo de Gauss-Newton com tolerância terminal de 10^-6 e número máximo de iterações igual a 700.

Foram obtidas as estimativas da inclinação nos dois pontos de inflexão (I_est e I_tex) por meio da derivada da curva ajustada conforme fórmulas propostas por Dexter & Bird (2001). Para o cálculo do parâmetro I_est usou-se n_est, U_sat e U_pmp, e para I_tex utilizou-se n_tex, U_pmp e U_res. Para isso, considerou-se o fato de que a inclinação da curva devida ao primeiro segmento para potenciais próximos do ponto de inflexão dado pelo segundo segmento é praticamente zero, e vice-versa (Figura 1).

Para avaliação das relações entre as estimativas dos parâmetros do modelo sugerido denominado duplo van Genuchten e estimativa da inclinação estrutural e textural (I_est; I_tex) para ambos os segmentos da curva de retenção de água, com as propriedades texturais do solo (conteúdo de argila), aplicou-se o teste de correlação de Pearson. Para todas as inferências estatísticas considerou-se nível nominal de significância de 5 %.

RESULTADOS E DISCUSSÃO

Verificou-se ampla variação da distribuição de partículas por tamanho nas amostras dos diferentes solos (Quadro 2), o que está relacionado à composição do material de origem. Latossolos originados de sedimentos argilosos ou de material proveniente de intemperismo do basalto tendem a ser mais argilosos que aqueles desenvolvidos de arenito. Essa relação torna-se evidente na composição da fração areia (> 0,05 mm), composta basicamente do mineral quartzo (Ferreira et al., 1999; Embrapa, 2006), que é extremamente resistente ao intemperismo e está presente no arenito e ausente no basalto (Curi & Franzmeier, 1987). Dessa forma, sugere-se que esses dois materiais parentais em diferentes proporções originaram os solos estudados.

As relações moleculares Ki e Kr são utilizadas na estimativa da mineralogia e classificação dos Latossolos (Embrapa, 2006). A relação Kr variou de 0,2 a 0,7 para os horizontes Bw estudados; por conseguinte, todos os Latossolos foram classificados como oxídicos (Quadro 2).

O processo de ajuste do modelo duplo van Genuchten convergiu para todos os solos. O modelo proposto mostrou alto poder explicativo, apresentando coeficiente de determinação (R²) acima de 99 % em todos os solos estudados. Conforme observado na figura 2, o modelo caracteriza bem a relação entre os conteúdos de água no solo dentro do intervalo de potencial matricial avaliado. Fica evidente a relação não linear entre as variáveis expressas por esse modelo, inclusive para potenciais acima de 1.500 kPa, caracterizada também por apresentar um ponto de inflexão nessas condições.

A equação duplo van Genuchten apresenta comportamento semelhante ao do modelo da equação dupla exponencial proposto por Dexter et al. (2008) para solos de regiões de clima temperado e ao escalonamento fractal de Alfaro Soto et al. (2008) para solos tropicais, embora considerassem intervalos de potenciais matriciais diferentes para o estabelecimento da curva de retenção de água.

Assim como para o modelo exponencial e o escalonamento fractal, a derivada do modelo duplo van Genuchten apresenta uma função densidade bimodal para o tamanho de poros do solo e caracteriza a porosidade deste em poros estruturais e texturais. Além do mais, o modelo proposto foi desenvolvido com base na fundamentação teórica da equação de van Genuchten (1980), o qual apresenta extensa aplicabilidade na avaliação do conteúdo de água no solo retido sob baixos potenciais (1 < Ψ_m < 1.500 kPa), o que faz dele potencial para estudos de retenção de água sob altos potenciais (Ψ_m > 1.500 kPa), haja vista também a disponibilidade de software para realização dessa modelagem matemática.

De acordo com Ferreira et al. (1999), em solos tropicais altamente intemperizados, os óxidos de ferro e de alumínio atuam como agentes cimentantes entre as partículas minerais, alterando o arranjamento dos seus constituintes. Micromorfologicamente, os Latossolos estudados apresentam a distribuição dos minerais em relação ao plasma, seguindo o padrão aglutinado, em agregados granulares muito pequenos (diâmetro < 300 µm), conferindo elevado índice de vazios, por sua vez, organizado em cavidades interconectadas (Vollant-Tuduri et al., 2005; Severiano, 2010). Consequentemente, a distribuição de poros por tamanho nesses solos é caracterizada por apresentar predominantemente duas classes de poros distintos: a primeira, relacionada àqueles muito grandes, que promovem rápida drenagem interna do solo; e a segunda, a poros muito pequenos, nos quais a água é retida com energia muito alta, caracterizando-a como água higroscópica (Oliveira et al., 2004).

O comportamento bimodal da curva de retenção de água no solo reflete essas classes de poros, em que o primeiro ponto de inflexão ocorre sob baixos potenciais matriciais (entre 1 e 3 kPa), ao passo que o segundo ponto ocorre sob potenciais muito altos (entre 10.000 e 20.000 kPa) (Figura 2; Quadro 3). Para solos de clima temperado, o comportamento bimodal tem sido verificado dentro do intervalo da curva-padrão de retenção de água no solo, isto é, na faixa de 1 a 1.500 kPa (Dexter & Richard, 2009), em razão da distribuição de poros uniforme, comparativamente aos solos de regiões de clima tropical.

Quando inseridos no processo produtivo, os poros de elevado diâmetro podem ser seriamente afetados logo nas primeiras passagens do maquinário agrícola (Secco et al., 2004), dada a elevada suscetibilidade à compactação desses solos (Ajayi et al., 2009; Severiano, 2010). Segundo Oliveira et al. (2007), a compressão do solo decorrente do manejo promove o incremento do potencial matricial, em função da redução da macroporosidade e do aumento na porosidade de retenção de água, pela transformação de parte desta em microporos. Por esse motivo, têm sido observados valores de potencial matricial, para o primeiro ponto de inflexão, superiores aos encontrados neste trabalho para solos pertencentes à mesma classe, quando submetidos ao uso intensivo (Mello et al., 2002).

As estimativas dos parâmetros foram significativas, com exceção do conteúdo residual de água (U_res) para os solos L₁₅₂, L₂₀₉ e L₂₆₃ (Quadro 3). Esse resultado indica que toda a água desses solos foi removida, considerando a amplitude de leitura de potencial matricial avaliada pelo aparelho WP4-T. Em contrapartida, os Latossolos mais argilosos, por apresentarem maior proporção de partículas finas, promovem incremento na energia de retenção de água, observando-se ainda pequena quantidade de água residual adsorvida na matriz do solo mesmo em potenciais de 300.000 kPa. Da mesma forma, ao se avaliar o valor do conteúdo de água no ponto de murcha permanente (U_pmp), observa-se o mesmo comportamento.

O que foi discutido no parágrafo anterior no que diz respeito ao U_pmp é de grande relevância, pois, embora este seja limítrofe para o metabolismo da maioria das plantas cultivadas, em termos de estresse hídrico (Tormena et al., 2007), ele não é indicativo do melhor momento para trafegabilidade de solos agrícolas. Nesse sentido, Severiano (2010), ao comprovar que os Latossolos mais argilosos em estudo apresentam baixa capacidade de suporte de carga, algumas operação de manejo mecanizado, poderá provocar algum grau de compactação adicional, mesmo em condições de solo considerado seco (Ψ_m > 1.500 kPa).

O conteúdo total de argila tem sido apontado como o principal atributo diretamente relacionado à retenção de água em solos altamente intemperizados, em razão da baixa atividade desses minerais e dos baixos conteúdos de C orgânico (Silva et al., 2008), comparativamente aos solos de clima temperado (Resende et al., 2007). Esse comportamento pode ser visualizado pela alta correlação, obtida mesmo com tamanho de amostra igual a 10, entre todos os parâmetros relacionados ao conteúdo de água no solo retido (U_sat, U_pmp e U_res) do modelo do duplo van Genuchten e o conteúdo de argila dos Latossolos (Figura 3). Esses resultados corroboram os de Beutler et al. (2002), que demonstraram maior retenção de água em solos com maiores conteúdo de argila. Por outro lado, não se obteve qualquer evidência de relação linear entre a argila e os parâmetros empíricos de ajuste do modelo (α_est, n_est, α_tex e n_tex).

A retenção de água no solo é influenciada pelos atributos intrínsecos do solo, dentre os quais se destacam a textura, a estrutura e o conteúdo de C orgânico (Giarola et al., 2002; Dexter et al., 2004; Severiano et al., 2010). Segundo Rawls et al. (1991), a retenção de água no solo sob baixos potenciais matriciais é decorrente do processo de capilaridade, sendo, portanto, fortemente alterada pelo arranjamento das partículas do solo, devido à presença de poros estruturais (poros inter-agregados). Por outro lado, com o secamento do solo, a composição granulométrica e mineralógica assume maior importância, o que se deve à maior superfície específica para a adsorção das moléculas de água (Gupta & Larson, 1979; Machado et al., 2008). Por esse motivo, neste trabalho, o valor da inclinação referente ao primeiro ponto de inflexão (baixos potenciais) é tratado como inclinação estrutural (I_est), ao passo que aquela referente à segunda inflexão (altos potenciais), como inclinação textural (I_tex), sendo os valores apresentados no quadro 3.

A análise de correlação de Pearson entre o valor da inclinação nos pontos de inflexão da curva de retenção de água e o conteúdo de argila está apresentada na figura 4. Segundo Dexter (2004), o valor da inclinação da primeira inflexão, definido como índice S, é considerado um indicador de qualidade estrutural, uma vez que a degradação física do solo pode ser detectada por meio de suas alterações.

Dessa forma, esse parâmetro sofre influência de todos os fatores que alteram a estrutura do solo, o que justifica a não significância da sua correlação (I_est) com o teor de argila (Figura 4), uma vez que, sob condições naturais, os Latossolos estudados apresentam elevada porosidade e friabilidade, em decorrência de sua estrutura microgranular, refletindo na sua ótima fertilidade física (Oliveira et al., 2004; Reatto et al., 2007; Severiano, 2010). Por esse motivo, sugere-se que outros atributos estruturais, não quantificados neste estudo, estejam influenciando a inclinação do primeiro ponto de inflexão, corroborando Streck et al. (2008), que não observaram relação do valor de S com o conteúdo de argila do solo.

Observa-se que todos os valores de I_est foram muito superiores ao considerado limítrofe por Dexter & Richard (2009) a uma condição estrutural muito boa (S > 0,05) (Quadro 3; Figura 4), o que é justificado pela ausência de atividade antrópica e pelas características intrínsecas do solo no horizonte estudado. Esses elevados valores refletem melhor configuração de poros no solo, o que favorece sobremaneira o desenvolvimento das culturas (Streck et al., 2008), e corroboram aqueles encontrados por Da Silva (2004) em Latossolo Vermelho-Amarelo sob condições naturais.

Em contrapartida, para I_tex, verifica-se relação direta com o conteúdo de argila (Figura 4), o que é reflexo do comportamento bimodal da distribuição de poros (Dexter & Bird, 2001; Dexter et al., 2008), evidenciando o desenvolvimento da porosidade textural dos Latossolos estudados, não sendo, portanto, considerado um indicador de qualidade estrutural do solo.

Assim, os Latossolos de textura média, em decorrência da sua menor superfície específica, apresentaram discreto segundo ponto inflexão (I_tex), ao passo que o aumento da proporção de partículas finas do solo promoveu, gradativamente, incremento na adsorção de água, em consequência da alta energia em que esta fica retida nos poros intramicroagregados, resultando na sua maior inclinação (Figuras 2 e 4).

CONCLUSÕES

1. O modelo duplo van Genuchten ajusta-se adequadamente aos dados experimentais com alto poder preditivo.

2. As curvas de retenção de água dos Latossolos em estudo apresentaram comportamento bimodal no intervalo de potencial matricial estudado.

3. Os parâmetros da equação relacionados ao conteúdo de água (U_sat, U_pmp, U_res) apresentaram forte dependência com o conteúdo de argila, assim como a inclinação no segundo ponto de inflexão (I_tex).

4. O incremento do conteúdo de argila dos Latossolos influenciou na retenção de água, aumentando a capilaridade e adsorção de água.

LITERATURA CITADA

Recebido para publicação em fevereiro de 2010 e aprovado em dezembro de 2010.

• AJAYI, A.E.; DIAS JUNIOR, M.S.; CURI, N.; ARAUJO JUNIOR, C.F.; SOUZA, T.T.T. & INDA JUNIOR, A.V. Strength attributes and compaction susceptibility of Brazilian Latosols. Soil Till. Res., 105, 1:122-127, 2009.
• ALBUQUERQUE, J.A.; MAFRA, A.L.; FONTOURA, S.M.V.; BAYER, C. & PASSOS, J.F.M. Avaliação de sistemas de preparo e calagem em um Latossolo Bruno alumínico. R. Bras. Ci. Solo, 29:963-975, 2005.
• ALFARO SOTO, M.A.; KIANG, C.H. & VILAR, O.M. Avaliação do escalonamento fractal de alguns solos brasileiros. R. Bras. Geoci., 38:253-262, 2008.
• BEUTLER, A.N.; CENTURION, J.F.; SOUZA, Z.M.; ANDRIOLI, I. & ROQUE, C.G. Retenção de água em dois tipos de Latossolos sob diferentes usos. R. Bras. Ci. Solo, 26:829-834, 2002.
• CURI, N. & FRANZMEIER, D. Effect of parent rocks on chemical and mineralogical properties of some Oxisols in Brazil. Soil Sci. Soc. Am. J., 51:153-158, 1987.
• DA SILVA, A.P. Interação manejo versus física. In: REUNIÃO BRASILEIRA DE MANEJO E CONSERVAÇÃO DO SOLO E DA ÁGUA, Santa Maria, 2004. Palestra. Santa Maria, 2004. CD-ROM.
• DECAGON DEVICES. Operator's manual version 1.3 WP4-T dewpointmeter. Pullman, 2000. 70p.
• DEXTER, A.R. Soil physical quality Part I. Theory, effects of soil texture, density, and organic matter, and effects on root growth. Geoderma, 120:201-214, 2004.
• DEXTER, A.R. & BIRD, N.R.A. Methods for predicting the optimum and the range of soil water contents for tillage based on the water retention curve. Soil Till. Res., 57:203-212, 2001.
• DEXTER, A.R.; CZY^–, E.A.; RICHARD, G. & RESZKOWSKA, A. A user-friendly water retention function that takes account of the textural and structural pore spaces in soil. Geoderma. 143:243-253, 2008.
• DEXTER, A.R. & RICHARD, G. Tillage of soils in relation to their bi-modal pore size distributions. Soil Till. Res., 103:113-118, 2009.
• DIAS JUNIOR, M.S. Compactação do solo. In: NOVAIS, R.F.; ALVAREZ V., V.H. & SCHAEFER, C.E. Topicos em ciência do solo. Viçosa, MG, Sociedade Brasileira de Ciência do Solo, 2000. v.1. p.55-94.
• DOURADO NETO, D.; NIELSEN, D.R.; HOPMANS, J.W.; REICHARDT, K. & BACCHI, O.O.S. Software to model soil water retention curves (SWRC, version 2.0). Sci Agric., 57:191-192, 2001.
• EMPRESA BRASILEIRA DE PESQUISA AGROPECUÁRIA - EMBRAPA. Centro Nacional de Pesquisa de Solos. Manual de métodos de análise de solo. 2.ed. Rio de Janeiro: Ministério da Agricultura e Abastecimento, 1997. 212p.
• EMPRESA BRASILEIRA DE PESQUISA AGROPECUÁRIA - EMBRAPA. Centro Nacional de Pesquisa de Solos. Sistema brasileiro de classificação de solos. 2.ed. Rio de Janeiro, Ministério da Agricultura e Abastecimento, 2006. 306p.
• FERREIRA, M.M.; FERNADES, B. & CURI, N. Influência da mineralogia da fração argila nas propriedades físicas de Latossolos da região sudeste do Brasil. R. Bras. Ci. Solo. 23, 3:515-524, 1999.
• GIAROLA, N.F.B.; SILVA, A.P. & IMHOFF, S. Relações entre propriedades físicas e características de solos da região sul do Brasil. R. Bras. Ci. Solo, 26:885-893, 2002.
• GUPTA, S.C. & LARSON, W.E. Estimating soil water retention characteristics from particle size distribution, organic matter percent, and bulk density. Water Res. Res., 15:1633-1635, 1979.
• KLEIN, V.A. & LIBARDI, P.L. Densidade e distribuição do diâmetro dos poros de um Latossolo Vermelho sob diferentes sistemas de uso e manejo. R. Bras. Ci. Solo, 26:857-867, 2002.
• KLUTE, A. Water retention: laboratory methods. In: KLUTE, A., ed. Methods of soil analysis: Physical and mineralogical methods. 2.ed. Madison, American Society of Agronomy, 1986. p.635-660.
• KÖPPEN, W.P. Grundriss der Klimakunde. Berlin, W. Gruyter, 1931. 388p.
• MACHADO, J.L.; TORMENA, C.A.; FIDALSKI, J. & SCAPIM, C.A. 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. R. Bras. Ci. Solo. 32:495-502, 2008.
• MELLO, C.R.; OLIVEIRA, G.C.; RESCK, D.V.S.; LIMA, J.M. & DIAS JÚNIOR, M.S. Estimativa da capacidade de campo baseada no ponto de inflexão da curva característica. Ci. Agrotec., 26, 4:836-841, 2002.
• MOREIRA, F.M.S. & SIQUEIRA, J.O. Microbiologia e bioquímica do solo. 2.ed. Lavras, Universidade Federal de Lavras, 2006. 729p.
• MUALEM, Y. A new model for predicting the hydraulic conductivity of unsaturated porous media. Water Res. Res., 12;513-522, 1976.
• OLIVEIRA, G.C.; SEVERIANO, E.C. & MELLO, C.R. Dinâmica da resistência à penetração de um Latossolo Vermelho da microrregião de Goiânia, GO. R. Bras. Eng. Agríc. Amb., 11:265-270, 2007.
• OLIVEIRA, G.C.; DIAS JUNIOR, M.S.; RESCK, D.V.S. & CURI, N. Caracterização química e físico-hídrica de um Latossolo Vermelho após vinte anos de manejo e cultivo do solo. R. Bras. Ci. Solo. 28:327-336, 2004.
• R DEVELOPMENT CORE TEAM (2009). R: A language and environment for statistical computing. R Foundation for Statistical Computing. Vienna, Austria. Disponível em:<http://www.R-project.org>. Acesso em 02 de fev. 2009.
• RAWLS, W.J.; GISH, T.J. & BRAKENSIEK, D.L. Estimating soil water retention from soil physical properties and characteristics. Adv. Soil Sci., 16:213-234, 1991.
• REATTO, A.; BRUAND, A.; MARTINS, E.S.; MULLER, F.; SILVA, E.M.; CARVALHO JUNIOR, O.A. & BROSSARD, M. Hydraulic properties of the diagnostic horizon of Latosols of a regional toposequencia across the Brazilian central platea. Geoderma, 139:251-59, 2007.
• RESENDE, M.; CURI, N.; REZENDE, S.B. & CORRÊA, G.F. Pedologia: Base para distinção de ambientes. 5.ed. Lavras, Universidade Federal de Lavras, 2007. 322p.
• RICHARDS, L.A. & WEAVER, L.R. Fitten-atmosphere percentage as related to the permanent wilting percentage. Soil Sci., 56:331-339, 1943.
• RITZ, C. & STREIBIG, J.C. Nonlinear regression with R. New York, Springer, 2008.
• RODRIGUES, A.C. Germinação e reindução da tolerância à dessecação em sementes de Bauhinia forficata LINK (FABACEAE). Lavras, Universidade Federal de Lavras, 2010. 64p. (Dissertação de Mestrado)
• SCALON, B.R.; ANDRASKI, B.J. & BILSKIE, J. Water potential: miscellaneous methods for measuring matric or water potential. In: DANE, J.H.; TOPP, C., eds. Methods of soil analysis: Physical methods. Madison, Soil Science Society of America, 2002. Part 4. p.643-670.
• SECCO, D.; REINERT, D.J.; REICHERT, J.M. & DA ROS, C.O. Produtividade e soja e propriedades físicas de um Latossolo submetido a sistemas de manejo e compactação. R. Bras. Ci. Solo, 28:797-804, 2004.
• SEVERIANO, E.C. Alterações estruturais de Latossolos representativos da região do cerrado e potencial de uso de solos cultivados com cana-de-açúcar. Lavras, MG, Universidade Federal de Lavras, 2010. 148p. (Tese de Doutorado).
• SILVA, A.P.; TORMENA, C.A.; FIDALSKI, J. & IMHOFF, S. Funções de pedotransferência para as curvas de retenção de água e de resistência do solo à penetração. R. Bras. Ci. Solo, 32:1-10, 2008.
• STRECK, C.A.; REINERT, D.J.; REICHERT, J.M. & HORN, R. Relações do parâmetro s para algumas propriedades físicas de solos do sul do Brasil. R. Bras. Ci. Solo, 32:2603-2612, 2008. (Número Especial)
• TORMENA, C.A.; ARAÚJO, M.A.; FIDALSKI, J. & COSTA, J.M. Variação temporal do intervalo hídrico ótimo de um Latossolo Vermelho distroférrico sob sistemas de plantio direto. R. Bras. Ci. Solo, 31:211-219, 2007.

*

Parte da Dissertação de Mestrado do primeiro autor apresentada ao programa de Pós-Graduação em Ciência do Solo da Universidade Federal de Lavras – UFLA.

Datas de Publicação

• Publicação nesta coleção
  25 Abr 2011
• Data do Fascículo
  Fev 2011

Histórico

• Aceito
  Dez 2010
• Recebido
  Fev 2010