Simulação do deslocamento de potássio em colunas verticais de solo não-saturado

Potassium displacement simulation in vertical columns of unsaturated soil

ARTIGOS CIENTÍFICOS

ENGENHARIA DE ÁGUA E SOLO

Simulação do deslocamento de potássio em colunas verticais de solo não-saturado

Potassium displacement simulation in vertical columns of unsaturated soil

Jarbas H. Miranda^I; Sergio N. Duarte^II; Paulo L. Libardi^III; Marcos V. Folegatti^IV

^IEngº. Agrônomo, Prof. Dr., Departamento de Ciências Exatas, Escola Superior de Agricultura "Luiz de Queiroz" - ESALQ/USP, Av. Pádua Dias, 11, Piracicaba - SP, Fone: (0XX19) 3429.4283 ramal: 210, [e-mail]

^IIProf. Dr. Departamento de Engenharia Rural, ESALQ/USP, [e-mail]

^IIIProf. Titular, Departamento de Ciências Exatas, ESALQ/USP, [e-mail]

^IVProf. Associado, Departamento de Engenharia Rural, ESALQ/USP, [e-mail]

RESUMO

O estudo do transporte de água e potássio em solo não-saturado é importante, tanto do ponto de vista do ambiente quanto do econômico. Assim sendo, o uso da modelagem computacional é importante, pois permite de maneira precisa e rápida o monitoramento do deslocamento de solutos, importante na prevenção de impactos ao ambiente. No presente trabalho, teve-se o objetivo de avaliar a simulação do deslocamento do íon potássio em colunas de solo não-saturado, utilizando o modelo MIDI, bem como apresentar a determinação dos parâmetros de transporte do íon potássio em um Latossolo Vermelho-Amarelo, fase arenosa. Concluiu-se que o modelo foi capaz de simular de maneira satisfatória o perfil de umidade e o deslocamento do íon potássio.

Palavras-chave: dinâmica de solutos, modelagem, curva de eluição.

ABSTRACT

Water and solute transport studies in unsaturated soil are important for both economical and environmental points of view and, in this sense, it should be emphasized the increase of agricultural use of urban and industrial residues, to the water resources and fertilizers saving. Thus, the computational modeling use is important, because it allows the monitoring of solute displacement, necessary to the environmental impacts prevention in a precise and fast way. The main objective of the present work is to simulate the displacement of potassium ion in unsaturated soil columns, using the MIDI model, as well as to present transport parameters determination of the potassium ion in a Red Yellowish Latossol, sandy phase. The obtained results allowed concluding that the model was capable to adequately simulate the potassium ion displacement.

Keywords: miscible displacement, modeling, breakthrough curve.

INTRODUÇÃO

O transporte simultâneo de água e de solutos em condições de fluxo não-saturado é objeto de interesse, variado e importante, pois a eficiência do uso da água e de produtos químicos nas produções agrícolas é por demais significativa, seja do ponto de vista do ambiente, seja econômico. Além disso, a possibilidade do uso de resíduos urbanos ou industriais, como forma de corrigir a deficiência nutricional de solos, tem sido considerada bastante promissora (MATOS et al., 1999).

O potássio, segundo MALAVOLTA (1980), é um macronutriente que tem caminhamento no solo, predominantemente pelo processo de difusão. Entretanto, alta parcela desse nutriente pode movimentar-se no solo por fluxo de massa. No Brasil, a aplicação de potássio às culturas é feita quase que totalmente no plantio e/ou cobertura, em aplicação direta ao solo ou via fertirrigação. Algumas vezes, a aplicação é complementada com adubações foliares (ZANINI, 1991).

ZANINI (1991) verificou que, de maneira geral, os locais de maiores concentrações de K⁺ coincidiram com os locais de maiores valores de umidade, evidenciando seu caminhamento por fluxo de massa, concordando com MALAVOLTA (1980), isto é, a distribuição de potássio no solo correlacionou-se com a distribuição de água no solo, indicando que se pode ter elevado controle da localização desse íon no solo em função da fertirrigação e da irrigação.

Para o sucesso da simulação e resolução das equações que permitem predizer o deslocamento de solutos no solo, é necessária uma determinação realista dos parâmetros de transporte que influenciam na relação solo-soluto. Os parâmetros mais importantes que devem ser determinados para esse fim, são a velocidade da água no poro, os coeficientes de difusão e dispersão e o fator de retardamento, sendo esse último definido por VALOCCHI (1984) como sendo a relação entre a velocidade do soluto reativo e a velocidade média da água no poro. O coeficiente de difusão-dispersão é um parâmetro físico que considera a difusão iônica e a dispersão mecânica na solução do solo (GENUCHTEN & WIERENGA, 1986).

Portanto, devido ao grande número de variáveis que envolvem o transporte de solutos no solo, a utilização de modelos de simulação torna-se desejável pela rapidez, precisão e por permitir que grande número de fatores e efeitos sejam contabilizados (MIRANDA, 2001). Vários modelos de lixiviação têm sido desenvolvidos, por exemplo, MOUSE (STEENHUIS et al., 1984), GLEAMS (KNISEL et al., 2000), HYDRUS (KOOL & GENUCHTEN,1991), MIDI (MIRANDA, 2001). O modelo MIDI é um modelo unidimensional que consiste na resolução de equações diferenciais parciais de segunda ordem, ou seja, a equação diferencial geral que rege o movimento de água no solo e a equação diferencial difusão-convecção do transporte de solutos.

O presente trabalho traz como objetivo avaliar a simulação do deslocamento do íon potássio em colunas de solo não-saturado, utilizando-se do modelo MIDI, bem como apresentar a determinação dos parâmetros de transporte do íon potássio em um Latossolo Vermelho-Amarelo, fase arenosa.

MATERIAL E MÉTODOS

Descrição do modelo

O modelo computacional MIDI foi desenvolvido visando à simulação da dinâmica de solutos no solo por meio de soluções numéricas de equações diferenciais que descrevem esse transporte. O modelo matemático implementado consiste na solução de duas equações diferenciais parciais de segunda ordem, ou seja, a equação do movimento de água no solo que descreve a variação da distribuição de umidade e a equação do transporte de solutos que permite estimar as mudanças na concentração de solutos no solo com o tempo.

As taxas de transferência de água entre células contíguas são calculadas pela equação de Darcy- Buckingham, que pode ser apresentada da seguinte forma:

em que,

q - fluxo de água entre as seções da coluna, L T^-1;

K(q) - condutividade hidráulica do solo no anel, função da umidade média () do solo no anel, LT^-1;

Dft - diferença do potencial total da água no solo, L, e

DZ - altura de cada anel da coluna, L.

O valor da condutividade hidráulica do solo não-saturado, ou seja, a função K(q), é obtido a partir da condutividade hidráulica do solo saturado (K_o) e dos parâmetros da curva de retenção de água no solo, de acordo com o modelo de MUALEM (1976) e GENUCHTEN (1980), apresentados na eq.(2).

em que,

K(q) - condutividade hidráulica do solo não-saturado, L T^-1;

Ko - condutividade hidráulica do solo saturado, L T^-1;

w - saturação efetiva, adimensional;

- parâmetro empírico, considerado igual a 0,5 para a maioria dos solos, e

m - parâmetro da curva de retenção da água no solo, adimensional.

O valor da saturação efetiva é obtido pela eq.(3).

em que,

q - umidade volumétrica do solo na condição inicial, L³ L^-3;

q_r - umidade volumétrica residual do solo, L³ L^-3, e

q_S - umidade volumétrica do solo no ponto de saturação, L³ L^-3.

Na superfície superior da primeira célula, em que está sendo aplicada a solução de potássio, o modelo considera que os valores da condutividade hidráulica e da umidade inicial correspondem ao ponto de saturação; em seguida, é determinado o potencial mátrico, fornecido pela eq.(4):

em que,

a - parâmetro com dimensão igual ao inverso da tensão, L^-1, e

n - parâmetro da curva de retenção de água no solo, adimensional.

Com os valores obtidos do potencial mátrico adicionados ao potencial gravitacional de cada ponto, considerados na metade de cada célula da coluna, o modelo determina o gradiente de potencial entre as células e, assim, estima o fluxo de água entre as mesmas [eq.(1)], considerando os valores de K(q) médios entre as células.

De posse dos valores de fluxo de água e estando no primeiro intervalo de tempo de simulação, o modelo determina a variação de umidade que ocorreu em cada ponto da coluna, cujo valor é obtido pela eq.(5):

em que,

Dq - variação de umidade nas células da coluna, L³ L^-3;

Dq - diferença entre o fluxo que entrou e saiu da célula, L T^-1;

A - área da seção transversal da coluna, L²;

DT - intervalo de tempo simulado, T, e

V - volume de cada anel da coluna, L³.

Essa variação de umidade é acrescida a cada célula da coluna, quando o modelo iniciar o segundo intervalo de tempo de processamento e assim por diante, até atingir o tempo total de simulação.

Uma vez determinados os valores de fluxo de água entre as células da coluna, inicia-se a determinação da concentração de potássio em cada célula. O primeiro valor a ser obtido é o fluxo de massa do potássio, resultado do produto entre o fluxo de água, entre os anéis da coluna, já obtido, e a concentração de potássio da solução no solo que preenche a coluna.

Conforme já visto, o processamento do modelo baseia-se na discretização da coluna em células. A diferença de concentração (DC) utilizada para a determinação do fluxo de potássio por difusão-dispersão é determinada dentro de cada simulação, no respectivo intervalo de tempo. Esse fluxo entre duas células contíguas da coluna é representado por J_dh, M L^-2 T^-1, e apresentado na eq.(6):

Na eq.(6), D é o coeficiente de dispersão hidrodinâmica longitudinal, L² T^-1, dado por:

em que,

D_m - coeficiente de difusão do soluto no solo na umidade q, L² T^-1, e

D_disp - coeficiente de dispersão do soluto para o fluxo q e a umidade q, L² T^-1.

O coeficiente Dm é obtido pela eq.(8) a seguir:

em que,

- umidade média entre células consecutivas, L³ L^-3;

q_S - umidade volumétrica do solo no ponto de saturação, L³ L^-3, e

Do - coeficiente de difusão molecular do soluto em água pura, L² T^-1.

O coeficiente D_disp é calculado por:

em que,

l - dispersividade do meio poroso, parâmetro empírico que varia com a textura do solo e a distância coberta pelo soluto, L.

O fluxo total de soluto entre as células da coluna é aditivo, ou seja, é o somatório do fluxo de massa e o fluxo de difusão-dispersão.

Para determinar a concentração do soluto em cada célula, é necessário determinar a quantidade de soluto que se movimenta, isto é, a quantidade de soluto que "entra" e "sai" em cada célula da coluna. Essas quantidades são calculadas pelo programa MIDI, realizando um balanço de massa.

Experimento em laboratório

O trabalho de validação do modelo foi conduzido no laboratório de Física do Solo do Departamento de Engenharia Rural - ESALQ/USP, conforme esquema composto pelos componentes apresentados nas Figuras 1 e 2. As soluções foram preparadas em concentrações de 500 kg m^-3 de potássio (KCl), sendo o solo classificado como Latossolo Vermelho-Amarelo, fase arenosa.

As análises físico-hídricas do solo apresentaram os seguintes resultados: 69,5% de areia; 12,0% de silte; 18,5% de argila; densidade do solo = 1,56 e densidade de partículas = 2,575; porosidade total do solo = 39,41% e condutividade hidráulica do solo saturado = 18,19 cm h^-1. Os parâmetros da curva de retenção, ajustada segundo o modelo de GENUCHTEN (1980), foram: q₀ = 0,0216 cm³ cm^-3; q_r = 0,162 cm³ cm^-3; q_s = 0,443 cm³ cm^-3; a = 0,0449 cm^-1; n = 3,6732, e m = 0,1119. A análise química foi determinada por amostras extraídas do material de solo do Latossolo Vermelho-Amarelo, fase arenosa (Tabela 1).

A coluna utilizada foi confeccionada em PVC e dividida em dez anéis, todos com dimensões de 5 cm de diâmetro e 7 cm de altura. O preenchimento foi feito com volumes iguais na tentativa de reduzir efeitos de fluxos preferenciais e da existência de camadas compactadas. Para impedir que houvesse vazamento da solução entre os anéis para a parte externa da coluna, utilizou-se de uma fita adesiva para a união dos anéis (Figura 2).

A aplicação das soluções foi feita por meio de um dosador de soro, para evitar variações na vazão de entrada na coluna de solo. Primeiramente, procedeu-se à regulagem do volume a ser aplicado por unidade de tempo, o qual ficou em torno de 150 mL h^-1. Após essa regulagem do dosador, dava-se início à aplicação da solução.

Devido ao fato de a coluna ser de PVC, não era possível acompanhar a frente de molhamento, razão pela qual foi necessário fazer vários testes prévios para determinar o intervalo de tempo ideal, no qual a solução não atingia o final da coluna. O intervalo de tempo estipulado foi de 1h 45.

Após transcorrido o intervalo de tempo estipulado, desmontava-se a coluna e, à medida que isso era feito, coletava-se uma pequena amostra de cada anel para a determinação da umidade do solo na respectiva coordenada da coluna. O restante de solo de cada anel era utilizado para o preparo da pasta saturada, cujos valores de umidade ficavam próximos ao ponto de saturação. Após o período de 24 h, foram extraídas as soluções e conduzidas ao laboratório para a determinação das concentrações do íon potássio, realizada mediante a utilização de sensores de íons específicos.

Após as determinações da concentração do íon potássio no extrato de saturação, realizou-se, com três repetições, a correção dos valores obtidos para a umidade, antes do preparo da pasta de extração da solução [eq.(10)]:

em que,

C_{íon corrigida} - concentração do íon potássio corrigida em função da umidade, g cm^-3;

U_pasta - umidade com base em peso da pasta utilizada para a extração da solução do solo, g g^-1;

C_íon - concentração do íon (potássio) determinada na extração da solução, g cm^-3, e

U_anel - umidade com base em peso do solo em cada anel da coluna, após o término da aplicação da solução, g g^-1.

RESULTADOS E DISCUSSÃO

Curvas de eluição

A determinação dos parâmetros de transporte relativo ao solo utilizado foi feita mediante a utilização do modelo CXTFIT, desenvolvido pelo U.S. Salinity Laboratory - USDA - Riverside - CA, versão 2.1 (TORIDE, 1999).

Os valores obtidos dos parâmetros de transporte do potássio foram: velocidade da água nos poros (V) = 1,672 cm min^-1; coeficiente de dispersão (D) = 2,703 cm² min^-1; fator de retardamento (R) = 4,687; dispersividade (l) = 1,61662 cm; coeficiente de difusão (Do) a 25 ºC = 0,001188 cm² min^-1; concentração média inicial do soluto no solo () = 49,31 kg m^-3, e concentração da solução aplicada de cada soluto (Cs) = 500 kg m^-3. A representação gráfica dos valores experimentais, bem como dos ajustados pelo modelo CXTFIT, encontra-se na Figura 3.

Validação da estimativa da umidade realizada pelo modelo

Na Figura 4(A), apresentam-se os perfis de umidade medidos em três repetições e, na Figura 4(B), o perfil médio medido, além daquele simulado pelo modelo. Verifica-se que o modelo MIDI simulou, de maneira uniforme, o perfil de umidade resultante da aplicação da solução de potássio. As variações apresentadas são devidas, possivelmente, aos efeitos de preenchimento da coluna.

Validação da estimativa da concentração de potássio realizada pelo modelo

Em relação à simulação da concentração do íon potássio, pode-se observar maior variação entre as repetições, quando comparado com os dados de umidade. O solo que apresentava concentração inicial de 49,31 kg m^-3, ao término da aplicação da solução de potássio, revelou os valores médios que ficaram em torno de 300 kg m^-3, nos anéis superiores. Há tendência, portanto, de maiores concentrações de potássio nas camadas superiores, em torno de 14 cm, enquanto o molhamento se deu até cerca de 30 cm. Esse resultado, possivelmente, deve-se ao fato de o potássio em solução ter sido adsorvido, de forma que a frente de molhamento que caminhou para profundidades maiores, provavelmente, apresentava menor concentração desse soluto (Figura 5).

CONCLUSÕES

A simulação realizada pelo modelo MIDI foi capaz de prever, com bom ajuste, o perfil de umidade e o de deslocamento do íon potássio, incluindo o retardamento de movimentação desse íon em relação à frente de molhamento. Portanto, para que a simulação atinja padrões cada vez mais confiáveis, é necessário ressaltar a importância da determinação dos parâmetros de transporte, principalmente quando se trata da previsão do deslocamento de solutos no solo.

Recebido pelo Conselho Editorial em: 7-6-2004

Aprovado pelo Conselho Editorial em: 20-9-2005

• GENUCHTEN, M.T. Van. A closed-form equation for predicting the hydraulic conductivity of unsaturated soils. Soil Science Society of America Journal, Madison, v.44, n.3, p.892-8, 1980.
• GENUCHTEN, M.T. Van; WIERENGA, P.J. Solute dispersion coefficients and retardation factors. In: BLACK, C.A. (Ed.) Methods of soil analysis Madison: Soil Science Society of America, 1986. p.1025-54. (American Society of Agronomy, 9)
• KNISEL, W.G; DAVIS, F.M. GLEAMS - Groundwater loading effects of agricultural management systems. Version 3.0. User Manual. Tifton: U.S. Department of Agriculture, Agricultural Research Service, Southeast Watershed Research Laboratory, 2000. 194 p.
• KOOL, J.B.; GENUCHTEN, M.T. Van. HYDRUS one-dimensional variably saturated flow and transport model, including hysteresis and root water uptake. Riverside: U.S. Salinity Laboratory USDA-ARS, 1991. 101 p.
• MALAVOLTA, E. Potássio, magnésio e enxofre nos solos e culturas brasileiras 2.ed. Piracicaba: Instituto Internacional de Pesquisa da Potassa, 1980. 91 p. (Boletim Técnico, 4)
• MATOS, A.T.; COSTA, L.M.; FONTES, M.P.F.; MARTINEZ, M.A. Retardation factors and the dispersion-diffusion coefficients of Zn, Cd, Cu and Pb in soils from Viçosa-MG, Brazil. Transactions of the ASAE, St. Joseph, v.42, n.4, p.903-10, 1999.
• MIRANDA, J.H. Modelo para simulação da dinâmica de nitrato em colunas verticais de solo não saturado. 2001. 79 f. Tese (Doutorado em Irrigação e Drenagem) - Escola Superior de Agricultura "Luiz de Queiroz", Universidade de São Paulo, Piracicaba, 2001.
• MUALEM, Y. A new model for predicting the hydraulic conductivity of unsaturated porous media. Water Resources Research, New York, v.12, n.3, p.513-22, 1976.
• STEENHUIS, T.S.; HUGHES, H.B.F.; PACENKA, S.; GROSS, M. MOUSE User's Manual. New York: Northeast Regional Agricultural Engineering Service, Cornell University, 1984. 44 p.
• TORIDE, N.; LEIJ, F.J.; GENUCHTEN, M.T. Van. The CXTFIT Code for Estimating Transport Parameters from Laboratory or Field Tracer Experiments. Version 2.0. Riverside: U.S. Salinity Laboratory, 1999. 119 p. (Research Report, 137)
• VALOCCHI, A.J. Describing the transport of ion-exchanging contaminants using an effective Kd approach. Water Resources Research, St. Louis, v.20, n.4, p.499-503, 1984.
• ZANINI, J.R. Distribuição de água e do íon K⁺ no solo, aplicados por fertirrigação em gotejamento. II - Teores de K⁺ no bulbo molhado. ITEM - Irrigação e Tecnologia Moderna, Brasília, v.46, n.1, p.24-38, 1991.

Datas de Publicação

• Publicação nesta coleção
  08 Mar 2006
• Data do Fascículo
  Dez 2005

Histórico

• Aceito
  20 Set 2005
• Recebido
  07 Jun 2004