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

Use of the S-index to evaluate soil physical quality

NOTA

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

Use of the S-index to evaluate soil physical quality

Celsemy Eleutério Maia

Engenheiro-Agrônomo, DS em Recursos Naturais. Professor do Departamento de Ciências Ambientais e Tecnológicas da Universidade Federal Rural do Semiárido - Ufersa. BR 110, km 47, Costa e Silva, CEP 59625-900 Mossoró (RN). Bolsista do CNPq. E-mail: [e-mail]

RESUMO

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).

Termos de indexação: curva de retenção de água, qualidade do solo, modelagem.

SUMMARY

Nowadays the concept of soil physical quality is well-established, based mainly on the indices that quantify soil quality. The S index is a new indicator of soil physical quality related to important soil characteristics. The purpose of this study was to evaluate the method of obtaining the S index and the derivative of the function of the of water retention curve, using the transformation ln(h) or not, and assess the problem of units and the restriction m = 1 - 1/n for the index calculation. It was concluded that the S index value is higher when the soil moisture used to adjust the water retention curve of van Genuchten (1980) is expressed as volumetric than as gravimetric moisture. For the same unit of soil moisture, no difference was observed in the S index for any tension unit. The S index value differs according to the use or not use of the restriction m = 1 -1/n to estimate the parameters m and n of the van Genuchten retention curve.

Index terms: water retention curve, soil quality, modeling.

INTRODUCTION

Para avaliar a qualidade física dos solos, Dexter (2004) propôs um novo método, denominado de índice S, baseado no ponto de inflexão do modelo de curva de retenção da água proposta por van Genuchten (1980). Segundo Dexter (2004), maiores valores de S (maior inclinação) indicam melhor distribuição do tamanho de poros, condizente com condições estruturais que estabelecem um adequado funcionamento físico do solo. Assim, fatores relacionados com o solo (textura) e com o manejo (matéria orgânica e estrutura), que influenciam a distribuição do tamanho de poros, refletem em mudanças nos valores de S e, portanto, na qualidade física do solo. Este autor propõe valores de S > 0,035 e S < 0,035 para áreas com condições favoráveis e restritivas para o crescimento das raízes das plantas, respectivamente; valores de S < 0,020 indicam condições físicas do solo muito pobres e com alta restrição ao crescimento das raízes das plantas. Assim, o índice S apresenta sensibilidade para identificar a degradação da qualidade física do solo em diferentes sistemas de uso do solo. O maior valor S na floresta, por exemplo, indicou a melhor configuração de poros no solo e, portanto, menor restrição física para o crescimento de raízes das plantas, seja por aeração, por restrição mecânica ou por características de retenção de água, quando comparado com o de áreas cultivadas (Streck et al., 2008). Entretanto, Silva et al. (2010) afirmam que os valores do índice S proposto por Dexter (2004) são provisórios para classificar a qualidade física dos solos. Afirmam ainda que estudos realizados pela Universidade Estadual de Maringá em 2006 comprovaram a eficiência do índice S na avaliação da qualidade física do solo sob diferentes sistemas de uso e manejo, porém alertaram que o índice ainda carece de uma quantidade considerável de pesquisas que comprovem a sua eficácia como indicador da qualidade física do solo.

Segundo Dexter (2004), a determinação do valor de S pode ser feita diretamente a partir da curva de retenção de água, com a plotagem da umidade do solo em função do potencial matricial. Entretanto, é mais conveniente ajustar a curva a uma função matemática e, pela derivada, calcular a inclinação no ponto de inflexão com base nos parâmetros da função; isso porque a curva de retenção de água proposta por van Genuchten (1980) apresenta um ponto crítico, observado pela forma gráfica dela quando os valores do módulo do potencial matricial (h) são transformados usando ln(h). Contudo, igualando a segunda derivada da equação 1 a zero, verifica-se que, mesmo sem a transformação ln(h), o modelo proposto por van Genuchten (1980) apresenta ponto de inflexão. Nesse sentido, sem a transformação ln(h), observa-se que a derivada primeira da equação 1 é dada pela equação 2, diferindo da obtida por Dexter (2004) usando ln(h) - equação 3.

em que: θ = umidade do solo (m³ m^-3 ou kg kg^-1); θ_s = umidade de saturação (m³ m^-3 ou kg kg^-1); θ_r = umidade residual do solo (m³ m^-3 ou kg kg^-1); h = módulo do potencial matricial (MPa, kPa, hPa, cm H₂O); α, n e m = parâmetros do modelo, com a em 1/tensão (MPa, kPa, hPa, cm H₂O); e n e m são adimensionais.

Dexter (2004), para derivar a equação 1 usando ln(h), fez uma transformação de variável, ln(h) = y → h = e^y; usando a regra da cadeia e sendo igual à equação 2 e , obtém-se a equação 3.

Os pontos críticos usando valores de h e ln(h) podem ser observados na figura 1, indicando que nos dois casos existe ponto de inflexão, porém que os valores desse ponto usando h são inferiores quando comparado com ln(h) para tensão em kPa. O módulo do potencial matricial onde se observa o ponto de inflexão é dado pelas equações 4 e 5, para dados não transformados ou usando ln(h), respectivamente, e assumindo a restrição m = 1 - 1/n, para n > 1.

Conceitualmente, a umidade do solo no potencial matricial igual a zero é a umidade de saturação (θ_s), o que não é válido quando utilizada a transformação ln(h). Além disso, dependendo da unidade da tensão, os valores transformados usando ln(h) serão negativos quando os valores de tensão estiverem entre 0 e 1. Nesse caso, seguindo o sistema internacional de unidades, usando a tensão em MPa, os valores variam de 0 até 1,5 MPa e todos os valores entre 0 e 1 apresentam ln(h) < 0; quando se usa kPa (os valores variam de 0 a 1.500 kPa), a curva desloca-se para o primeiro quadrante dos eixos cartesianos, mas uma parte dela ainda permanece com valores de ln(h) < 0, o que pode ser solucionado quando se usa a tensão em hPa (0 a 15.000 hPa) (Figura 2), como utilizado por Dexter (2004), porém sem a definição sobre a tensão zero equivaler à umidade de saturação do solo, devido ao problema da descontinuidade da função no ponto h=0.

Substituindo a equação 4 em 2 e a 5 em 3, têm-se as equações 6 e 7, sendo esta última proposta por Dexter (2004) para o cálculo do índice S. Nos dois casos, as equações calculam a declividade da reta tangente ao ponto de inflexão da curva de retenção de água, utilizando os parâmetros do ajuste da curva conforme modelo proposto por van Genuchten (1980).

Observando as equações 6 e 7, levando em consideração que n e m são adimensionais, a unidade da equação 6 é umidade de θ por unidade de h, ao passo que a da 7 é apenas umidade (m³ m^-3 ou kg kg^-1) e que, matematicamente, a diferença é o parâmetro α presente no numerador da equação 6 e ausente na 7. Sendo a unidade de α igual a h^-1, dividindo a equação 6 por α, a unidade passa a ser apenas a da umidade e a equação 6 se iguala a 7.

Nesse sentido, trabalhando com unidades diferentes, tanto de umidade quanto de tensão, o valor da derivada deve mudar. Mesmo assim, assumindo que a unidade de S é a da umidade do solo, de acordo com a equação 7, esta deveria ser expressa nos resultados, porém o que se observa é que o índice S é expresso como se fosse adimensional.

O objetivo deste trabalho foi avaliar o método do índice S quanto à derivada da função de retenção de água da curva proposta por van Genuchten (1980) 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.

MATERIAL E MÉTODOS

Os dados utilizados neste trabalho foram obtidos de Moreira & Silva (1987), localizado na Unidade de Execução e Pesquisa (UEP) da Empresa Pernambucana de Pesquisa Agropecuária (IPA), no município de Goiana, PE, com caracterização granulométrica de acordo com o quadro 1.

Com os dados das amostras indeformadas, foram ajustadas as curvas de retenção de água para cada camada do solo usando o modelo proposto por van Genuchten (1980), de acordo com a equação 1, assumindo ou não a restrição m = 1 - 1/n e usando a umidade do solo em m³ m^-3 ou kg kg^-1 e as tensões (h) em MPa, kPa, hPa e cmH₂O. O índice S foi calculado pela equação 7, e a declividade no ponto de inflexão sem a transformação ln(h), pela equação 6.

RESULTADOS E DISCUSSÃO

Ajustando os dados ao modelo de retenção de água usando como variável dependente a umidade em m³ m^-3 e kg kg^-1, e como variável independente a tensão de água no solo (h), observa-se (Quadro 2) que os valores do índice S em módulo foram diferentes, sendo superiores quando a umidade foi expressa em m³ m^-3, e que estes não variaram com a unidade da tensão na mesma profundidade. Isso implica, para comparação dos índices S em diferentes estudos, as umidades devem ser expressas na mesma unidade (m³ m^-3 ou kg kg^-1). Dexter & Bird (2001), avaliando o ponto de inflexão na curva de retenção, trabalharam com as unidades de hPa e kg kg^-1, para tensão e umidade do solo, respectivamente - mesmas unidades utilizadas por Dexter (2004) para propor o índice S. Outros autores também calcularam o índice utilizando essas unidades, como Reynolds et al. (2009). A diferença do índice S calculado usando m³ m^-3 ou kg kg^-1 é matematicamente observada na equação 7. Considerando que, no ajuste dos dados para uma mesma unidade de tensão, os valores de α, n e m não se alteram, o mesmo não acontece com as umidades θ_s e θ_r, que são inferiores quando expressas em kg kg^-1, comparadas com a umidade em m³ m^-3, alterando o numerador da equação 7 e, consequentemente, o índice S (Quadro 2). Apesar de os valores de S em módulo calculados usando m³ m^-3 terem sido maiores quando comparados com aqueles em kg kg^-1, é importante salientar que isso se deve ao fato de a densidade do solo (Ds) ser superior a 1,0 Mg m^-3 para todas as camadas avaliadas. Dessa maneira, têm-se valores dos índices S superiores em m³ m^-3, quando comparados com os em kg kg^-1 , e isso ocorre se Ds > 1,0 Mg m^-3, e inferiores se Ds < 1,0 Mg m^-3.

O uso da unidade seria importante para efeito de comparação dos índices. Como matematicamente o índice S pode ser determinado dividindo dθ/dh_pi por α, sendo dθ/dh_pi dado em kg kg^-1 kPa^-1 e α em kPa^-1, isso demonstra também que a unidade do índice S é a unidade de umidade, no caso, kg kg^-1, que ao ser multiplicada pela densidade do solo resultaria na expressão do índice em m³ m^-3.

Comparando os valores do índice S com relação à mudança na unidade de tensão (h), observa-se (Quadro 2) que eles não variaram com as tensões avaliadas (MPa, kPa, hPa e cm H₂O), porém, matematicamente, era de se esperar uma mudança, como observado para os valores de dθ/dh_pi. Independentemente da unidade de tensão, os valores do índice S foram sempre negativos, mas, por conveniência, muitos autores adotam o valor absoluto de S, como recomendado por Dexter (2004).

Ainda com relação à unidade de h, avaliando os valores de dθ/dh_pi, estes foram maiores em valores absolutos quando utilizado o MPa, diminuindo o valor da inclinação no ponto de inflexão na sequência: MPa > kPa > hPa ≅ cm H₂O, tanto para umidade em m³ m^-3 ou kg kg^-1.

Outra forma de avaliar a transformação ln(h) foi ajustando os valores de umidade em função de ln(h) utilizando a curva de retenção proposta por van Genuchten (1980) e estimando os valores de θ_s, θ_r, α, n e m, para obtenção de dθ/dh_pi e do índice S, de acordo com as equações 6 e 7, respectivamente. Nesse sentido, observa-se no quadro 3 que os valores de dθ/dh_pi assemelham-se aos de S do quadro 2.

Outro ponto importante a ser observado no cálculo do índice S é a restrição proposta por Mualem (1976), que é indicada no ajuste dos dados ao modelo da curva de retenção proposta por van Genuchten (1980). Assim, no quadro 4, observa-se que para a umidade, tanto em m³ m^-3 quanto em kg kg^-1, o valor do índice S alterou quando comparado com os de S do quadro 1. Assim, é necessário indicar para o cálculo do índice S o uso da restrição m = 1 - 1/n, pois ela pode apresentar valores diferentes com ou sem uso dessa restrição no cálculo do índice S. Como exemplos, Streck et al. (2008) e Silva et al. (2008) não mencionam a restrição, mas comparam seus resultados com os de outros autores. O valor do índice S de 0,035, proposto por Dexter (2004) como ponto divisor para separar solos com problemas de degradação, foi calculado usando a umidade do solo em kg kg^-1 e a restrição m = 1 - 1/n, para solos com D_s > 1,0 Mg m^-3. Anteriormente ao trabalho de Dexter (2004), Dexter & Bird (2001) já haviam trabalhado com a umidade em kg kg^-1 e com a restrição m = 1 - 1/n para o cálculo do ponto de inflexão na curva de retenção.

Observando os valores do índice S com base na umidade em m³ m^-3 e kg kg^-1, para as camadas do solo analisadas, verifica-se para a mesma camada discrepância nos valores do índice S de acordo com as unidades dos dados de entrada e da restrição, sendo mais pronunciada a diferença dos valores de S na camada mais superficial (Quadro 5). Assim, recomenda-se que no cálculo do índice S para amostras indeformadas, a umidade do solo deve ser expressa em kg kg^-1 e usando a restrição m = 1 - 1/n, como sugerido por Dexter (2004) e usado por Reynolds et al. (2009). No quadro 6 são apresentadas as unidades da umidade, da tensão e o uso ou não da restrição m = 1 - 1/n no cálculo do índice S, chamando a atenção os trabalhos de Andrade & Stone (2009a,b), que utilizaram para o cálculo do índice S a umidade do solo em kg kg^-1 e m³ m^-3, respectivamente (Quadro 6); contudo, como demonstrado aqui, dependendo da densidade do solo, os valores do índice S diferem com a utilização de umidade com unidade diferente.

CONCLUSÕES

1. Usando a tensão sem transformação logarítmica, verificou-se que a curva de retenção proposta por van Genuchten também apresenta ponto de inflexão semelhante ao observado quando usada a transformação ln(h), porém os valores dos índices S são diferentes, tanto para o uso ou não da restrição m = 1 - 1/n.

2. 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.

3. Para a mesma unidade da umidade do solo, não há diferença nos valores do índice S para qualquer unidade de tensão.

4. O uso ou não da restrição m = 1 - 1/n para estimar os parâmetros m e n da curva de retenção de água proposta por van Genuchten (1980) confere diferente valores de índice S.

5. Há necessidade de padronização no uso da umidade em kg kg^-1 e uso da restrição m = 1 - 1/n para comparação dos valores de referência sugeridos por Dexter (2004) para avaliar a qualidade física dos solos, usando o índice S.

LITERATURA CITADA

• AGUIAR, M.I. Qualidade física do solo em sistemas agroflorestais. Viçosa, MG, Universidade Federal de Viçosa, 2008. 79p. (Tese de Mestrado)
• ANDRADE, R.S. & STONE, L.F. Índice S como indicador da qualidade física de solos do Cerrado brasileiro. R. Bras. Eng. Agri. Amb., 13:382-388, 2009a.
• ANDRADE, R.S. & STONE, L.F. Uso do índice S na determinação da condutividade hidráulica não-saturada de solos do cerrado Brasileiro. R. Bras. Eng. Agrí. Amb., 13:376-381, 2009b.
• BEUTLER, A.N.; FREDDI, O.S.; LEONE, C.L. & CENTURION, J.F. Densidade do solo relativa e parâmetro "S" como indicadores da qualidade física para culturas anuais. R. Bio. Ci. Terra, 8:27-36, 2008.
• DEXTER, A.R. Soil physical quality. Part I. Theory, effects of soil texture, density, and organic matter, and effects on root growth. Geoderma, 20: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 Tillage Res., 57:203-212, 2001.
• FIDALSKI, J. & TORMENA, C.A. Homogeneidade da qualidade física do solo nas entrelinhas de um pomar de laranjeira com sistemas de manejo da vegetação permanente. R. Bras. Ci. Solo, 31:637-645, 2007.
• MATIAS, S.S.R.; BORBA, J.A.; TICELLI, M.; PANOSSO, A.R. & CÂMARA, F.T. Atributos físicos de um Latossolo Vermelho submetido a diferentes usos. R. Ci. Agron., 40:331-338, 2009.
• MOREIRA, J.A.A. & SILVA, C.J.C.G. Característica de retenção de água de um solo Podzólico Vermelho-Amarelo de Goiana, Pernambuco. Pesq. Agropec. Bras., 22:411-418, 1987.
• MUALEM, Y. A new model for predicting the hydraulic conductivity of unsaturated porous media. Water Res. Res., 12:513-522, 1976.
• REYNOLDS, W.D.; DRURY,C.F.; TAN, C.S.; FOX, C.A. & YANG, X.M. Use of indicators and pore volume-function characteristics to quantify soil physical quality. Geoderma, 152:252-263, 2009.
• SILVA, F.F.; FREDDI, O.S.; CENTURION, J.F.; ARATANI, R.G.; ANDRIOLI, F.F. & ANDRIOLI, I. Propriedades físicas de um Latossolo Vermelho cultivado no sistema plantio direto. Irriga, 13:191-204, 2008.
• SILVA, A.P.; TORMENA, C.A.; DIAS JÚNIOR, M.S.; IMHOFF, S. & KLEIN, V.A. Indicadores da qualidade física do solo. In: JONG van LIER, Q., ed. Física do solo. Viçosa, MG, Sociedade Brasileira de Ciência do Solo, 2010. p.541-281.
• 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.
• van GENUCHTEN, M.T. A closed-form equation for predicting the hydraulic conductivity of unsatured soil. Soil Sci. Soc. Am. J., 44:892-898, 1980.

1

Recebido para publicação em 11 de abril de 2011 e aprovado em 13 de setembro de 2011.

Datas de Publicação

• Publicação nesta coleção
  10 Fev 2012
• Data do Fascículo
  Dez 2011

Histórico

• Recebido
  11 Abr 2011
• Aceito
  13 Set 2011