Sorção de Cd e Pb nos solos da região do Médio Rio Paraíba - RJ, Brasil1 Parte da Dissertação de Mestrado da primeira autora apresentada à Universidade Federal do Rural do Rio de Janeiro/UFRRJ

Sorption of Cd and Pb in soils of the Middle Region of Paraíba River - RJ, Brazil

Amanda Guimarães de Mattos Nelson Moura Brasil do Amaral Sobrinho Erica Souto Abreu Lima Jair do Nascimento Guedes Fábio Freire de Sousa Sobre os autores

RESUMO

As concentrações de metais na solução do solo são reguladas, principalmente, pelos processos de sorção. Neste estudo foi avaliada a dinâmica de sorção de Cd e Pb em 19 amostras de solos da Região do Médio Rio Paraíba - RJ, Brasil, através da comparação dos parâmetros de diferentes modelos de isotermas obtidos pelo software ISOFT após ensaios de sorção pelo método da batelada. Os resultados mostraram que a isoterma Linear é mais adequada para representar a sorção de Cd e o modelo de Freundlich para Pb e que a adsorção de Cd é fortemente influenciada pelo pH e pelas concentrações de Fe e Ca, enquanto que a adsorção de Pb é mais influenciada pelo pH, e pelas concentrações de Mn e Ca.

Palavras-chave:
Metais Pesados; Modelos de sorção; Isotermas

ABSTRACT

Metal concentrations in soil solution are mainly regulated by sorption processes. In this study, Cd and Pb sorption dynamics have been assessed for 19 soil samples from the Middle Region of Paraíba River - RJ, Brazil, through parameters of different models of isotherms obtained by the ISOFT software following a sorption testing using the batch method. Results showed that the linear isotherm is more appropriate to describe Cd sorption and the Freundlich model for Pb. Cd sorption is strongly affected by pH and concentrations of Fe and Ca; while Pb sorption is mostly influenced by pH and the concentrations of Mn and Ca.

Key words:
Heavy metals; Models of sorption; Isotherms

INTRODUÇÃO

Quando é detectada a presença de altas concentrações de metais pesados no solo, seja ela de fonte natural ou antropogênica, preocupa-se com os problemas que poderão surgir como: a entrada na cadeia trófica, redução da produtividade agrícola devido aos efeitos fitotóxicos, acúmulo no solo, alteração da atividade microbiana e contaminação dos recursos hídricos (PIRES; MATTIAZZO, 2007PIRES, A. M. M; MATTIAZZO, M. E. Cinética de solubilização de metais pesados por ácidos orgânicos em solos tratados com lodo de esgoto. Revista Brasileira de Ciência de Solo, v. 31, p. 143-151, 2007.). Entretanto, é necessário conhecer a dinâmica dos metais no solo, que é governada por uma série de reações químicas, para poder predizer seus riscos.

Os impactos provocados pela presença destes elementos estão associados, principalmente, à sua biodisponibilidade e mobilidade, ou seja, quando se encontram retidos com baixa energia de ligação ao solo e/ou colóides ou livres na solução. A sorção é considerada o mais importante fenômeno que regula as concentrações de metais na solução do solo (CAMARGO; ALLEONI; CASAGRANDE, 2001CAMARGO, O. A.; ALLEONI, L. R. F.; CASAGRANDE, J. C. Reações dos micronutrientes e elementos tóxicos no solo. In: FERREIRA, M. E. et al. Micronutrientes e elementos tóxicos na agricultura. Jaboticabal: CNPQ:FAPESP:POTAFOS, 2001. cap. 5, p. 89-124.).

Em estudos de fitotoxicidade de Cd e Pb em diferentes gramíneas, em solução nutritiva, Fonseca (2012)FONSECA, F. G. Tolerância de gramíneas forrageiras ao cádmio e biodisponibilidade no solo. 2012. 63 f. Dissertação (Mestrado em Produção Vegetal) - Universidade Federal dos Vales do Jequitinhonha e Mucuri, Diamantina, MG, 2012. e Nascimento (2013)NASCIMENTO, S. S. Tolerância de gramíneas forrageiras ao chumbo e sua disponibilidade no solo. 2013. 63f. Dissertação (Mestrado em Produção Vegetal) - Universidade Federal dos Vales do Jequitinhonha e Mucuri, Diamantina, MG, 2013., determinaram as doses críticas de toxicidade (DCT), ou seja, concentração do elemento na solução que promoveria redução de 10% do crescimento relativo das espécies. Verificaram para Cd teores entre 0,03 - 0,22 mg L-1 e para Pb entre 4,2 - 15,5 mg L-1. Entretanto, Paiva, Carvalho e Siqueira (2000)PAIVA, H. P.; CARVALHO, J. G.; SIQUEIRA, J. O. Efeito de Cd, Ni, Pb e Zn sobre mudas de cedro (Cedrela fissilis Vell.) e ipê-roxo (Tabebuia impetiginosa (Mart.) Standley) em solução nutritiva. Revista Árvore, v. 24, p. 369-378, 2000., para cedro e ipê-roxo, encontraram valores de DCT nas faixas de 0,03 - 0,14 e 0,05 - 0,13 mg L-1 para Cd e Pb, respectivamente. Esses autores concluíram que a DCT dependia da espécie a ser implantada e da capacidade de sorção desses metais pelo solo. Visto que, mudanças nas condições do solo, como alterações do pH, podem causar a mobilização dos metais pesados da fase sólida para a fase líquida. Além disso, a sorção é também influenciada pelos demais atributos do solo, destacando-se o teor e tipo de argila, capacidade de troca catiônica, teor de ferro, manganês e matéria orgânica (PIRANGUELI et al, 2005PIRANGUELI, M. A. P. et al. Efeito do pH na adsorção e dessorção de cádmio em Latossolos brasileiros. Revista Brasileira de Ciência de Solo, v. 29, p. 523-532, 2005.).

Em estudos de sorção, os resultados experimentais são descritos quantitativamente pelas isotermas de adsorção, que compreendem equações matemáticas utilizadas para a descrição de resultados experimentais de sorção de solutos por sólidos. São representadas graficamente, demonstrando as relações de equilíbrio entre a quantidade do elemento sorvido versus a quantidade remanescente na solução final (SPOSITO, 2008SPOSITO, G. The chemistry of soils. 2.ed. New York: Oxford University Press, 2008. 342p.). Através das isotermas é possível determinar a capacidade de sorção do solo, servindo também para a compreensão de outros mecanismos de interação, como a competição por sítios ou a heterogeneidade da superfície sólida do solo (BRADL, 2004BRADL, H. B. Adsorption of heavy metal ions on soils and soils constituents. Journal of Colloid and Interface Science, v. 277, p. 1-18, 2004.). Contudo, existem diferentes modelos matemáticos para a construção dos gráficos das isotermas, onde um modelo pode ser mais adequado que o outro para representar a sorção de diferentes metais.

A Região de governo do Médio Paraíba, localizada ao sul do Estado do Rio de Janeiro, é cortada pela Rodovia Presidente Dutra, que interliga as principais capitais do país: Rio de Janeiro e São Paulo. Apresenta uma atividade industrial muito intensa, abrigando importantes indústrias como a Companhia Siderúrgica Nacional (CSN), em Volta Redonda. O Rio Paraíba do Sul atravessa esta Região e abastece a Região Metropolitana através de sua transposição para o Rio Guandu.

Portanto, o objetivo deste trabalho foi selecionar o modelo de isoterma que melhor se ajusta aos dados de sorção de Cd e Pb e identificar os fatores que mais influenciam a sorção destes elementos em solos representativos da região do Médio Paraíba do Estado do Rio de Janeiro.

MATERIAL E MÉTODOS

Foram utilizadas 19 amostras superficiais (0 a 40 cm) de diferentes classes de solos (Latossolos, Argissolos, Cambissolos e Neossolos) coletadas em áreas de baixa interferência antrópica da Região do Médio Paraíba, localizada no Estado do Rio de Janeiro, Brasil.

A caracterização química e física foi feita em amostras de terra fina seca ao ar (TFSA), destorroadas em almofariz de ágata e passadas em peneira de aço inox de 2 mm e determinados os teores de carbono orgânico, pH em H2O (1:2,5), cátions trocáveis, acidez potencial, capacidade de troca catiônica (CTC) e granulometria (Embrapa ,1997EMPRESA BRASILEIRA DE PESQUISA AGROPECUÁRIA. Manual de métodos de análises de solos. 2. ed. Rio de Janeiro, 1997.).

A avaliação da sorção de Cd e Pb foi realizada no Laboratório de Química e Poluição do Solo, do Departamento de Solos da Universidade Federal Rural do Rio de Janeiro. Para os ensaios de sorção foi utilizado o método “batch’ ou método da batelada utilizado por Harter e Naidu (2001)HARTER, R. D.; NAIDU, R. An assessment of environmental and solution parameter impact on trace-metal sorption by soils. Soil Science Society of America Journal, v. 65, p. 597-612, 2001..

Para o procedimento da batelada foi tomado 1 g de terra fina seca ao ar (TFSA) passada em peneira de 0,150 mm, em pH natural, onde foram adicionados 20 ml de solução contendo concentrações crescentes de cada elemento, cujas diferentes concentrações de Cd e Pb foram estabelecidas a partir de ensaio preliminar. Visando manter a mesma força iônica em todas as soluções, estas foram preparadas em solução de KCl 0,02 mol L-1 utilizando-se água ultra pura coletada de sistema Milli-Q. Para o preparo da solução de Cd, foi utilizado CdCl2 nas concentrações crescentes de 0; 1; 2; 3; 4 e 5 mg L-1. Para Pb, utilizou-se Pb(NO3)2 nas concentrações de 0; 20; 50; 80; 100 e 130 mg L-1. Todas as concentrações foram realizadas em triplicata para cada amostra e cada concentração. Todo o procedimento foi mantido em temperatura ao redor de 25 °C.

Após, as amostras de solo contendo as soluções foram acondicionadas em tubos de polietileno (tipo falcon) e colocadas para agitar por 24 h em agitador horizontal a 120 rpm. Após agitação, foram realizadas as leituras do pH e colocadas para centrifugar a 3.000 rpm por 15 min e, posteriormente, filtradas em papel de filtro e armazenadas em frascos de polietileno. As determinações das concentrações dos elementos nos extratos foram realizadas em aparelho de espectrofotômetro de absorção atômica Variam EspectrAA 55B. De posse dos resultados dos extratos, a quantidade sorvida foi calculada, para cada amostra e concentração, através da equação 1:

onde: qe = é a quantidade sorvida de adsorbato por unidade de massa de adsorvente (mg kg-1); V = é o volume adicionado ao frasco contendo a concentração de metal (ml); M = é a massa de solo (g); Co = representa a concentração de metal inicialmente colocada em contato com o solo (mg L-1) e Ce = a concentração de equilíbrio da solução (mg L-1) (HARTER; NAIDU, 2001HARTER, R. D.; NAIDU, R. An assessment of environmental and solution parameter impact on trace-metal sorption by soils. Soil Science Society of America Journal, v. 65, p. 597-612, 2001.).

Através da utilização do software “ISOFIT” (Isotherm Fitting Tool), os parâmetros dos diversos modelos matemáticos de isotermas foram calculados de forma automatizada.

Os valores de Ce e qe foram utilizados como dados de entrada no software ISOFIT para a obtenção dos parâmetros das isotermas e critérios de seleção de seis modelos: BET, Langmuir, Freundlich, Freundlich-Langmuir, Linear e Toth. As equações referentes a cada modelo são apresentadas na Tabela 1.

Tabela 1
Equações referentes a cada modelo de isoterma de adsorção

Para a comparação dos modelos de isotermas foram utilizados como critérios o coeficiente de determinação (R2), Critério de Informação de Akaike corrigido (AICc) e peso de Akaike (AICw) (MATTOT; RABIDEAU, 2008MATTOT, L. S.; RABIDEAU, A. J. A program for fitting sorption isotherms to experimental data. Environmental Modelling & Software, v. 23, p. 670-676, 2008.). O peso de evidência ou peso de Akaike (AICw), ao contrário do R2 e do AICc, não é fornecido pelo ISOFIT, sendo calculado através da equação 2, proposta por Burnham e Anderson (2004)BURNHAM, K. P.; ANDERSON, D. R. Multimodel inference: understanding AIC and BIC in model selection. Sociological Methods and Research, v. 33, n. 2, p. 261-304, 2004.:

onde: AICcj = é o valor de AICc obtido pelo modelo; AICcmin = é o menor AICc dentre os modelos candidatos.

Para verificar a interação entre as variáveis e suas influências nos parâmetros das isotermas de sorção, foi realizada Análise de Componentes Principais e correlação de Pearson a 5% de probabilidade, calculados pelo Programa R Development Core Team (2013)PIRES, A. M. M; MATTIAZZO, M. E. Cinética de solubilização de metais pesados por ácidos orgânicos em solos tratados com lodo de esgoto. Revista Brasileira de Ciência de Solo, v. 31, p. 143-151, 2007..

RESULTADOS E DISCUSSÃO

Na Tabela 2, estão apresentados os atributos das 19 amostras de solo selecionadas. Os valores de pH indicam que a maioria das amostras, cerca de 80%, apresentam-se de fortemente ácidas (pH 4,4 a 5,3) a moderada acidez (pH 5,4 a 6,5) (RAIJ et al, 1996RAIJ, B. van. et al. Recomendações de adubação e calagem para o Estado de São Paulo. 2. ed. Campinas: Instituto Agronômico de Campinas, 1996. 285 p. (Boletim Técnico, 100).), porém as amostras S4 e S5 apresentam alcalinidade moderada, com pH em água em torno de 7,5. Em 63% das amostras a CTC é considerada média (valores entre 4,3 e 8,6 cmolc dm-3) (ALVAREZ VENEGAS; RIBEIRO, 1999ALVAREZ VENEGAS, V. H.; RIBEIRO, A. C. Calagem. In: RIBEIRO, A. C.; GUIMARÃES, P. T. G.; ALVAREZ VENEGAS. V. H. Recomendação para o uso de corretivos e fertilizantes em Minas Gerais: 5ª aproximação. Viçosa, MG: Imprensa Universitária, 1999. p. 43-60.). Os teores de ferro (Fe) e manganês (Mn) apresentaram grande variação, com valor mínimo de ferro de 12,11 g kg-1 e valor máximo de 45,53 g kg-1 e teores de manganês variando de 25,9 mg kg-1 a 760,2 mg kg-1. Os teores de argila também tiveram grande variação (36 a 503 g kg-1), sendo que nos solos S2, S5 e S7 o teor de argila foi superior a 400 g kg-1. Estes resultados irão representar uma capacidade adsortiva diferenciada entre os solos mais discrepantes, haja vista que os óxidos de ferro e manganês e os colóides da argila são importantes controladores da sorção de metais (ALVAREZ VENEGAS; RIBEIRO, 1999ALVAREZ VENEGAS, V. H.; RIBEIRO, A. C. Calagem. In: RIBEIRO, A. C.; GUIMARÃES, P. T. G.; ALVAREZ VENEGAS. V. H. Recomendação para o uso de corretivos e fertilizantes em Minas Gerais: 5ª aproximação. Viçosa, MG: Imprensa Universitária, 1999. p. 43-60.). Os teores de carbono orgânico (CO) apresentaram pouca variação, em que cerca de 63% das amostras apresentaram valores entre 1 e 2%, com valor mínimo de 0,35% e valor máximo de 3,29%, considerados, em geral, teores baixos (RAIJ et al., 1996RAIJ, B. van. et al. Recomendações de adubação e calagem para o Estado de São Paulo. 2. ed. Campinas: Instituto Agronômico de Campinas, 1996. 285 p. (Boletim Técnico, 100).). Os teores pseudototais de Cd em todas as amostras ficaram abaixo do limite de detecção do método, não sendo possível calcular os valores médio, máximo e mínimo. Para Pb, o valor médio obtido foi de 14,1 mg kg-1, com mínimo de 7,4 mg kg-1 e máximo de 24,9 mg kg-1.

Tabela 2
Atributos dos 19 solos, média, mediana, valores mínimos, máximos e desvio padrão

A comparação dos modelos de isotermas através do coeficiente de determinação (R2), Critério de Informação de Akaike corrigido (AICc) e peso de Akaike (AICw) estão apresentados nas Tabelas 3 e 4. Os valores de R2 para Cd e Pb não variaram entre os modelos para cada amostra. Dessa forma, não é adequada a utilização deste parâmetro para a seleção de modelos de sorção para os elementos Cd e Pb.

Tabela 3
Valores de R2, AICc e AICw para seis modelos de isotermas de sorção para Cd nas amostras de solo selecionadas
Tabela 4
Valores de R2, AICc e AICw para seis modelos de isotermas de sorção para Pb nas amostras de solo selecionadas

De acordo com Burnham e Anderson (2004)BURNHAM, K. P.; ANDERSON, D. R. Multimodel inference: understanding AIC and BIC in model selection. Sociological Methods and Research, v. 33, n. 2, p. 261-304, 2004. a utilização do Critério de Informação de Akaike corrigido (AICc) é recomendada quando a relação n/p (número de amostras/parâmetros do modelo) for pequena, ou seja, menor que 40, como no caso deste estudo. Através desse critério, o modelo que apresentou o menor valor dentre os demais, assumiu-se como o mais adequado para representar a sorção dos metais no solo, visto que esta distância é uma medida de aproximação entre o modelo verdadeiro (hipotético) e o modelo calculado, ou seja, o menor valor de AICc representa maior aproximação do real (BURNHAM; ANDERSON, 2004BURNHAM, K. P.; ANDERSON, D. R. Multimodel inference: understanding AIC and BIC in model selection. Sociological Methods and Research, v. 33, n. 2, p. 261-304, 2004.). O peso de evidência ou peso de Akaike (AICw) é uma outra alternativa para expressar a distância entre os modelos e ratificar a seleção do modelo mais adequado. A soma dos pesos equivale a 1, o valor que mais se aproximar de 1 representa o melhor modelo de ajuste.

Na sorção de Cd (Tabela 3), observou-se que 84% das amostras tiveram os menores valores de AICc obtidos para o modelo Linear, também denominado como modelo do coeficiente de distribuição (Kd). Entretanto, alguns outros modelos também se mostraram adequados em algumas amostras, como os modelos de Langmuir e Freundlich.

Os modelos que apresentaram valores de AICw mais próximos de 1, foram os mesmos que apresentaram os menores valores de AICc. De acordo com estes valores, optou-se pela utilização do modelo Linear para a descrição da sorção de Cd por ter sido o de melhor ajuste na maioria das amostras. Este modelo também se ajustou bem aos solos utilizados no estudo da sorção de Cd realizado por Oliveira et al. (2010a)OLIVEIRA, L. F. C. et al. Adsorção e deslocamento do íon cádmio em solos do cerrado. Revista Brasileira de Engenharia Agrícola e Ambiental, v. 14, n. 8, p. 848-855, 2010a.. Além disso, sua utilização justifica-se por ser matematicamente mais simples e facilmente medido e ainda ser recomendado pelos órgãos ambientais na avaliação da previsão de risco (ALLEONI; MELLO; ROCHA, 2009ALLEONI, L. R. F.; MELLO, J. W. V.; ROCHA, W. S. D. Eletroquímica, adsorção e troca iônica no solo. In: ALLEONI, L. R. F.; MELLO, V. F. Química e mineralogia do solo. Parte II - Aplicações. Viçosa: Sociedade Brasileira de Ciência do Solo, 2009. cap. 2, p. 69-130.).

Ao contrário do observado para Cd, o modelo Linear não foi o mais adequado para descrever a sorção de Pb. De acordo com os valores de AICc, o modelo de Freundlich mostrou-se mais adequado para 87,5% das amostras avaliadas (Tabela 4). Contudo, em algumas amostras, os modelos que melhor se ajustaram aos dados foram, além de Freundlich, os de BET e Langmuir e foram considerados igualmente adequados por apresentarem AICw com valores mais próximos de 1.

Para o estudo da sorção de Pb nas amostras selecionadas, a utilização do modelo empírico de Freundlich ou de Langmuir irão descrever adequadamente a dinâmica da sorção. Segundo Camargo, Alleoni e Casagrande (2001)CAMARGO, O. A.; ALLEONI, L. R. F.; CASAGRANDE, J. C. Reações dos micronutrientes e elementos tóxicos no solo. In: FERREIRA, M. E. et al. Micronutrientes e elementos tóxicos na agricultura. Jaboticabal: CNPQ:FAPESP:POTAFOS, 2001. cap. 5, p. 89-124., na maior parte dos artigos científicos, estes modelos são os mais utilizados para os estudos que envolvem isotermas de sorção e, além de serem simples, descrevem os resultados experimentais em ampla faixa de concentração (SPOSITO, 2008SPOSITO, G. The chemistry of soils. 2.ed. New York: Oxford University Press, 2008. 342p.). Portanto, escolheu-se o modelo de Freundlich por melhor descrever a sorção de Pb na maior parte das amostras deste estudo.

Os diferentes modelos de isotermas testados na sorção de Cd e Pb mostraram, claramente, que é necessário selecionar um modelo adequado para cada elemento em estudo, pois verifica-se uma dinâmica de sorção diferenciada entre eles. Estes resultados demonstram ainda que, apesar de o modelo Linear ou modelo do coeficiente de distribuição (Kd), ser recomendado pelos órgãos ambientais para o estudo de sorção de metais na análise de risco (CETESB, 2001COMPANHIA DE TECNOLOGIA DE SANEAMENTO AMBIENTAL. Relatório de estabelecimento de valores orientadores para solos e águas subterrâneas no Estado de São Paulo. São Paulo: CETESB, 2001. 73 p. (Relatório Técnico).), para Pb este parâmetro não descreveu adequadamente a sorção.

O parâmetro Kd do modelo Linear selecionado para Cd e parâmetros Kf e 1/n do modelo de Freundlich selecionado para Pb, foram obtidos pelo software ISOFIT e estão apresentados na Tabela 5. Os coeficientes Kd e Kf estão associados à capacidade de sorção do solo, onde os maiores valores indicam maior sorção do metal. Portanto, observa-se uma grande variação destes parâmetros entre as amostras de solos, com valores de Kd variando de 2,33 a 3.580 mg g-1 e Kf de 185,7 a 2.685,9 mg g-1 para Pb. O parâmetro 1/n é adimensional e demonstra se a sorção é favorável (valores entre 0 e 1) ou desfavorável, relacionando-se à distribuição dos sítios energéticos nos coloides do solo (SPOSITO, 2008SPOSITO, G. The chemistry of soils. 2.ed. New York: Oxford University Press, 2008. 342p.), sendo a sorção favorável em todas as amostras.

Tabela 5
Parâmetros Kd, Kf e 1/n das isotermas de sorção selecionadas para os elementos Cd e Pb

Verifica-se que a amostra S5, apresentou o valor de Kd muito superior às demais, refletindo uma maior retenção de Cd. Para Pb, as amostras S5, S7 e S15 não apresentaram teores de Pb nos extratos pelo fato de terem sorvido todo o Pb adicionado. Estas amostras também tiveram alta sorção de Cd, evidenciando elevada afinidade dos metais pelos sítios de sorção. Dessa forma, pode-se inferir que os valores de Kf, determinados para Pb, nestas três amostras são superiores ao maior valor, observado na amostra S14 (2.685,9 mg g-1). No estudo da capacidade de retenção de Pb em solos brasileiros, realizado por Pirangueli et al. (2001)PIRANGUELI, M. A. P. et al. Teor total e capacidade máxima de adsorção de chumbo em latossolos brasileiros. Revista Brasileira de Ciência de Solo, v. 25, p. 279-288, 2001., observaram que os Latossolos apresentaram grande capacidade em reter Pb.

De maneira geral, o Pb ficou mais retido ao solo que o Cd; estes resultados corroboram com os obtidos por alguns autores, que observaram maior sorção de Pb quando comparado ao Cd (LINHARES et al., 2009LINHARES, L. A. et al. Adsorção de cádmio e chumbo em solos tropicais altamente intemperizados. Pesquisa Agropecuária Brasileira, v. 44, n. 03, p. 291-299, 2009.; OLIVEIRA et al., 2010bOLIVEIRA, L. F. C. et al. Isotermas de sorção de metais pesados em solos do cerrado de Goiás. Revista Brasileira de Engenharia Agrícola e Ambiental, v. 14, n. 7, p. 776-782, 2010b.; PIRANGUELI et al., 2007PIRANGUELI, M. A. P. et al. Comportamento sortivo, individual e competitivo, de metais pesados em latossolos com mineralogia contrastante. Revista Brasileira de Ciência de Solo, v. 31, p. 819-826, 2007.; PIRANGUELI et al. , 2009PIRANGUELI, M. A. P. et al. Sorção de cádmio e chumbo em Latossolo Vermelho Distrófico sob efeito de calcário e fosfato. Revista Brasileira de Ciências Agrárias, v. 4, n. 1, p. 42-47, 2009.). Isto pode ser explicado pela energia de ligação diferente que estes elementos apresentam. Na retenção de Cd ocorrem predominantemente forças eletrostáticas, enquanto que no Pb ocorrem principalmente interações específicas, formando complexos internos, que apresentam uma maior energia de ligação (SPOSITO, 2008SPOSITO, G. The chemistry of soils. 2.ed. New York: Oxford University Press, 2008. 342p.). Além disso, o Cd possui um menor poder competitivo, como verificado por Pirangueli et al. (2007)PIRANGUELI, M. A. P. et al. Comportamento sortivo, individual e competitivo, de metais pesados em latossolos com mineralogia contrastante. Revista Brasileira de Ciência de Solo, v. 31, p. 819-826, 2007.. Estas características fazem com que a presença do Cd no solo determine maior risco do ponto de vista ambiental, por apresentar uma interação menos estável, ou seja, com uma maior mobilidade e biodisponibilidade, potencializando as chances de transferência para a cadeia alimentar e a contaminação das águas superficiais e subsuperficiais.

Em relação aos níveis críticos de fitotoxicidade de Cd e Pb, observou-se que para Cd apenas a amostra S6, de maior Kd (maior adsorção), apresentou todas as concentrações de equilíbrio (CE) abaixo da faixa de DCT para gramíneas (FONSECA, 2012FONSECA, F. G. Tolerância de gramíneas forrageiras ao cádmio e biodisponibilidade no solo. 2012. 63 f. Dissertação (Mestrado em Produção Vegetal) - Universidade Federal dos Vales do Jequitinhonha e Mucuri, Diamantina, MG, 2012.), onde os valores variaram de 0,006 - 0,027 mg L-1. Nos demais solos os teores ficaram acima da DCT, em todas as doses aplicadas. Para Pb, as amostras S4, S9, S11 e S14 com maior adsorção, apresentaram valores de Ce abaixo do DCT para gramíneas (NASCIMENTO, 2013NASCIMENTO, S. S. Tolerância de gramíneas forrageiras ao chumbo e sua disponibilidade no solo. 2013. 63f. Dissertação (Mestrado em Produção Vegetal) - Universidade Federal dos Vales do Jequitinhonha e Mucuri, Diamantina, MG, 2013.) e acima da DCT para cedro e ipê-roxo (PAIVA; CARVALHO; SIQUEIRA, 2000PAIVA, H. P.; CARVALHO, J. G.; SIQUEIRA, J. O. Efeito de Cd, Ni, Pb e Zn sobre mudas de cedro (Cedrela fissilis Vell.) e ipê-roxo (Tabebuia impetiginosa (Mart.) Standley) em solução nutritiva. Revista Árvore, v. 24, p. 369-378, 2000.).

Para o melhor entendimento de quais os atributos do solo, variáveis independentes ou explicativas, que mais influenciam nos parâmetros de sorção relacionados com a capacidade sortiva (Kd e Kf), foi realizada uma análise multivariada através de Componentes Principais (Figuras 1 e 2), visto que os atributos, em geral, exercem uma ação combinada que resultam em maior ou menor impacto ambiental devido ao controle da retenção de metais (PIRANGUELI et al., 2007PIRANGUELI, M. A. P. et al. Comportamento sortivo, individual e competitivo, de metais pesados em latossolos com mineralogia contrastante. Revista Brasileira de Ciência de Solo, v. 31, p. 819-826, 2007.).

Figura 1
Análise de Componentes Principais para o parâmetros Kd da isoterma Linear, selecionada para descrever a sorção de Cd nos solos da Região do Médio Rio Paraíba-RJ
Figura 2
Análise de Componentes Principais para os parâmetros Kf e 1/n da isoterma de Freundlich, selecionada para descrever a sorção de Pb nos solos da Região do Médio Rio Paraíba-RJ

Para verificar a contribuição individual destes atributos em relação aos parâmetros de sorção, foi realizada uma análise de correlação de Pearson (Tabela 5). Na análise de Componentes Principais, verifica-se que 52,1% da variância acumulada dos atributos do solo e do parâmetro Kd (Figura 1) foram explicadas nas duas primeiras componentes e no segundo gráfico (Figura 2), para os parâmetros Kfe 1/n, as duas primeiras componentes explicam 61,8% da variância acumulada.

Observa-se que o parâmetro Kd apresenta contribuições mais significativas dos teores de Fe e do pH para a sorção de Cd. Contudo, os teores de Mn, Ca e Mg, também expressaram efetiva contribuição na sorção deste elemento. Na análise de correlação de Pearson (Tabela 6), foi avaliada a relação linear do Kd com os atributos e verificou-se que este parâmetro apresenta correlação significativa e positiva com pH e Ca e correlação com o Fe de 38%, porém não significativa. Linhares et al. (2009)LINHARES, L. A. et al. Adsorção de cádmio e chumbo em solos tropicais altamente intemperizados. Pesquisa Agropecuária Brasileira, v. 44, n. 03, p. 291-299, 2009., também observaram correlação positiva do Cd com o pH e também com a CTC.

Tabela 6
Coeficientes de correlação linear de Pearson entre os parâmetros das isotermas selecionadas para Cd (Kd) e Pb (Kf e 1/n) e os atributos do solo mais significativos pela análise de componentes principais

Os resultados demonstram que estes atributos estão intimamente envolvidos na retenção desses metais. As amostras que mais sorveram Cd (S5, S15, S4 e S14) foram as que apresentaram os mais elevados valores de pH e/ou Fe, Mn, Ca e Mg, apresentados na caracterização do solo (Tabela 2).

A amostra S5 apresentou o maior valor de Kd (Tabela 5) representando a mais alta sorção de Cd, provavelmente o pH e o Fe foram os atributos que mais contribuíram para a maior retenção nesta amostra, visto que apresentou o pH mais elevado (7,6), dentre as amostras em estudo e o mais alto teor de Fe, de 45,53 g kg-1, sendo que o valor médio de Fe das amostras é de 27,71 g kg-1. Sipos et al. (2008)SIPOS, P. et al. Sorption of copper, zinc and lead on soil mineral phases. Chemosphere, v. 73, p. 461-469, 2008., verificaram um aumento da sorção e precipitação dos metais nas partículas minerais quando os solos apresentam condições alcalinas. Na amostra S15, além de apresentar pH 6,6, o teor de Ca foi o mais elevado (8,4 cmolcdm-3). Já a amostra S14, não apresentou pH elevado (5,7), entretanto o teor de Ca foi de 3,8 cmolcdm-3, acima da média das amostras e o teor de Mg teve o maior valor observado (2,3 cmolc dm-3). Contudo, aquelas que apresentaram valores baixos de Kd para Cd (S6, S12 e S19) tinham também baixos valores de pH e/ou Fe. A amostra S6, com o menor valor de Kd, apresenta o menor valor de pH (4,4) de todas as amostras avaliadas. Portanto, fica evidenciado que a retenção de Cd é controlada por mais de um atributo do solo.

De acordo com a Figura 2, os atributos que mais contribuem para a sorção de Pb, ou seja, no parâmetro Kf são: pH, e teores de Mn, Ca, Mg, K e P. Entretanto, a contribuição do Mn, que expressa o teor de óxidos de Mn, é bem maior que a contribuição de Fe. O Fe nas componentes 1 e 2 não expressou significativa contribuição para o parâmetro Kf. A análise de correlação de Pearson do Kf (Tabela 6) com os atributos foi muito semelhante ao resultado da análise multivariada, tendo correlação positiva e altamente significativa (<0,01%) com o Ca e também se correlacionou positivamente com o pH, e teores de P, Mg e Mn.

Como visto anteriormente, a maior parte das amostras que apresentaram as maiores retenções de Pb, coincidem com as que também apresentaram maior retenção de Cd. Porém, é mais evidente a contribuição do pH, Ca e dos óxidos de Mn, na sorção de Pb. Enquanto que para Cd, além do pH e Ca, os óxidos de Fe exercem maior influência. Dong et al. (2000)DONG, D. et al. Adsorption of Pb and Cd onto metal oxides and organic material in natural surface coatings as determined by selective extractions: new evidence for the importance of Mn and Fe oxides. Water Research, v. 34, p. 427-436, 2000., em estudo de sorção, observaram também que os óxidos de Fe são mais importantes no controle da sorção de Cd, comparado aos óxidos Mn. Entretanto, na sorção de Pb foi verificado que os óxidos de Mn são mais importantes.

A alta afinidade do Pb pelos óxidos de Mn também tem sido relatada por outros autores (BRADL, 2004BRADL, H. B. Adsorption of heavy metal ions on soils and soils constituents. Journal of Colloid and Interface Science, v. 277, p. 1-18, 2004.; MORENO et al., 2006MORENO, A. M. et al. Factors influencing lead sorption-desorption at variable added metal concentrations in Rhodoxeralfs. Chemosphere, v. 64, p. 758-763, 2006.). Sipos et al. (2008)SIPOS, P. et al. Sorption of copper, zinc and lead on soil mineral phases. Chemosphere, v. 73, p. 461-469, 2008., verificaram que as correlações de Fe e Mn com metais relacionam-se à formação de óxidos que apresentam uma alta capacidade de sorção.

De acordo com o estudo realizado por Mavropoulos et al. (2005)MAVROPOULOS, E. et al. Pb2+, Cu2+ and Cd2+ ions uptake by Brazilian phosphate rocks. Brazilian Chemical Society, v. 16, n. 1, p. 62-68, 2005., sobre a imobilização de Pb por rochas fosfatadas brasileiras, foi observado que minerais contendo Ca em sua composição podem imobilizar o Pb presente no solo. Desta forma, verificou-se a relação observada entre o Pb e Ca, onde as amostras que apresentaram teores mais elevados de Ca também adsorveram mais Pb. Os teores elevados de Ca observados em algumas amostras de solos podem estar relacionados com o material de origem, pois segundo o Levantamento Geoquímico realizado no Estado do Rio de Janeiro (CUNHA et al., 2000CUNHA, F. G. et al. Levantamento geoquímico do Estado do Rio de Janeiro. Rio de Janeiro: Companhia de Pesquisa de Recursos Minerais. Serviço Geológico do Brasil, 2000. 42 p.), as rochas gnáissicas do complexo Paraíba do Sul mostram-se enriquecidas em Ca.

O parâmetro 1/n apresentou posição oposta ao Kf como se observa no gráfico das componentes principais (Figura 2), tendo uma maior contribuição da argila. A correlação deste parâmetro foi positiva com a argila e negativa com o Ca e o Mn. Segundo Sposito (2008)SPOSITO, G. The chemistry of soils. 2.ed. New York: Oxford University Press, 2008. 342p., como o valor de n está relacionado às partículas coloidais do solo, provavelmente este foi o motivo deste parâmetro correlacionar-se positivamente com a argila.

Matos et al. (2001)MATOS, A. T. et al. Mobility of heavy metals as related to soil chemical and mineralogical characteristics of Brazilian soils. Environmental Pollution, v. 111, n. 3, p. 429-435, 2001. relatam que é comum encontrar resultados divergentes de correlação dos atributos do solo com a sorção de metais, devido às diferenças químicas, físicas e mineralógicas que os solos apresentam e a dinâmica de cada metal. Contudo, os resultados da correlação com o pH para Cd e Pb corroboram com os encontrados na literatura (ALVAREZ VENEGAS; RIBEIRO, 1999ALVAREZ VENEGAS, V. H.; RIBEIRO, A. C. Calagem. In: RIBEIRO, A. C.; GUIMARÃES, P. T. G.; ALVAREZ VENEGAS. V. H. Recomendação para o uso de corretivos e fertilizantes em Minas Gerais: 5ª aproximação. Viçosa, MG: Imprensa Universitária, 1999. p. 43-60.; PIRANGUELI et al, 2007PIRANGUELI, M. A. P. et al. Comportamento sortivo, individual e competitivo, de metais pesados em latossolos com mineralogia contrastante. Revista Brasileira de Ciência de Solo, v. 31, p. 819-826, 2007.).

CONCLUSÕES

1. A isoterma Linear apresentou-se mais adequada para avaliar a sorção de Cd e a isoterma de Freundlich para avaliar a sorção de Pb;

2. Os teores de Fe e o pH foram os atributos do solo que mais contribuíram para a sorção de Cd enquanto que a retenção de Pb é controlada principalmente pelo pH, e teores de Mn, Mg e Ca.

  • Parte da Dissertação de Mestrado da primeira autora apresentada à Universidade Federal do Rural do Rio de Janeiro/UFRRJ

REFERÊNCIAS

  • ALLEONI, L. R. F.; MELLO, J. W. V.; ROCHA, W. S. D. Eletroquímica, adsorção e troca iônica no solo. In: ALLEONI, L. R. F.; MELLO, V. F. Química e mineralogia do solo. Parte II - Aplicações. Viçosa: Sociedade Brasileira de Ciência do Solo, 2009. cap. 2, p. 69-130.
  • ALVAREZ VENEGAS, V. H.; RIBEIRO, A. C. Calagem. In: RIBEIRO, A. C.; GUIMARÃES, P. T. G.; ALVAREZ VENEGAS. V. H. Recomendação para o uso de corretivos e fertilizantes em Minas Gerais: 5ª aproximação. Viçosa, MG: Imprensa Universitária, 1999. p. 43-60.
  • BRADL, H. B. Adsorption of heavy metal ions on soils and soils constituents. Journal of Colloid and Interface Science, v. 277, p. 1-18, 2004.
  • BURNHAM, K. P.; ANDERSON, D. R. Multimodel inference: understanding AIC and BIC in model selection. Sociological Methods and Research, v. 33, n. 2, p. 261-304, 2004.
  • CAMARGO, O. A.; ALLEONI, L. R. F.; CASAGRANDE, J. C. Reações dos micronutrientes e elementos tóxicos no solo. In: FERREIRA, M. E. et al Micronutrientes e elementos tóxicos na agricultura Jaboticabal: CNPQ:FAPESP:POTAFOS, 2001. cap. 5, p. 89-124.
  • COMPANHIA DE TECNOLOGIA DE SANEAMENTO AMBIENTAL. Relatório de estabelecimento de valores orientadores para solos e águas subterrâneas no Estado de São Paulo São Paulo: CETESB, 2001. 73 p. (Relatório Técnico).
  • CUNHA, F. G. et al Levantamento geoquímico do Estado do Rio de Janeiro Rio de Janeiro: Companhia de Pesquisa de Recursos Minerais. Serviço Geológico do Brasil, 2000. 42 p.
  • DONG, D. et al Adsorption of Pb and Cd onto metal oxides and organic material in natural surface coatings as determined by selective extractions: new evidence for the importance of Mn and Fe oxides. Water Research, v. 34, p. 427-436, 2000.
  • EMPRESA BRASILEIRA DE PESQUISA AGROPECUÁRIA. Manual de métodos de análises de solos 2. ed. Rio de Janeiro, 1997.
  • FONSECA, F. G. Tolerância de gramíneas forrageiras ao cádmio e biodisponibilidade no solo 2012. 63 f. Dissertação (Mestrado em Produção Vegetal) - Universidade Federal dos Vales do Jequitinhonha e Mucuri, Diamantina, MG, 2012.
  • HARTER, R. D.; NAIDU, R. An assessment of environmental and solution parameter impact on trace-metal sorption by soils. Soil Science Society of America Journal, v. 65, p. 597-612, 2001.
  • LINHARES, L. A. et al Adsorção de cádmio e chumbo em solos tropicais altamente intemperizados. Pesquisa Agropecuária Brasileira, v. 44, n. 03, p. 291-299, 2009.
  • MATOS, A. T. et al Mobility of heavy metals as related to soil chemical and mineralogical characteristics of Brazilian soils. Environmental Pollution, v. 111, n. 3, p. 429-435, 2001.
  • MATTOT, L. S.; RABIDEAU, A. J. A program for fitting sorption isotherms to experimental data. Environmental Modelling & Software, v. 23, p. 670-676, 2008.
  • MAVROPOULOS, E. et al Pb2+, Cu2+ and Cd2+ ions uptake by Brazilian phosphate rocks. Brazilian Chemical Society, v. 16, n. 1, p. 62-68, 2005.
  • MORENO, A. M. et al Factors influencing lead sorption-desorption at variable added metal concentrations in Rhodoxeralfs. Chemosphere, v. 64, p. 758-763, 2006.
  • NASCIMENTO, S. S. Tolerância de gramíneas forrageiras ao chumbo e sua disponibilidade no solo 2013. 63f. Dissertação (Mestrado em Produção Vegetal) - Universidade Federal dos Vales do Jequitinhonha e Mucuri, Diamantina, MG, 2013.
  • OLIVEIRA, L. F. C. et al Adsorção e deslocamento do íon cádmio em solos do cerrado. Revista Brasileira de Engenharia Agrícola e Ambiental, v. 14, n. 8, p. 848-855, 2010a.
  • OLIVEIRA, L. F. C. et al Isotermas de sorção de metais pesados em solos do cerrado de Goiás. Revista Brasileira de Engenharia Agrícola e Ambiental, v. 14, n. 7, p. 776-782, 2010b.
  • PAIVA, H. P.; CARVALHO, J. G.; SIQUEIRA, J. O. Efeito de Cd, Ni, Pb e Zn sobre mudas de cedro (Cedrela fissilis Vell.) e ipê-roxo (Tabebuia impetiginosa (Mart.) Standley) em solução nutritiva. Revista Árvore, v. 24, p. 369-378, 2000.
  • PIRANGUELI, M. A. P. et al Efeito do pH na adsorção e dessorção de cádmio em Latossolos brasileiros. Revista Brasileira de Ciência de Solo, v. 29, p. 523-532, 2005.
  • PIRANGUELI, M. A. P. et al Comportamento sortivo, individual e competitivo, de metais pesados em latossolos com mineralogia contrastante. Revista Brasileira de Ciência de Solo, v. 31, p. 819-826, 2007.
  • PIRANGUELI, M. A. P. et al Teor total e capacidade máxima de adsorção de chumbo em latossolos brasileiros. Revista Brasileira de Ciência de Solo, v. 25, p. 279-288, 2001.
  • PIRANGUELI, M. A. P. et al Sorção de cádmio e chumbo em Latossolo Vermelho Distrófico sob efeito de calcário e fosfato. Revista Brasileira de Ciências Agrárias, v. 4, n. 1, p. 42-47, 2009.
  • PIRES, A. M. M; MATTIAZZO, M. E. Cinética de solubilização de metais pesados por ácidos orgânicos em solos tratados com lodo de esgoto. Revista Brasileira de Ciência de Solo, v. 31, p. 143-151, 2007.
  • R DEVELOPMENT CORE TEAM 3.0.1. R: a language and environment for statistical computing Vienna, Austria: R Foundation for Statistical Computing, 2013.
  • RAIJ, B. van. et al Recomendações de adubação e calagem para o Estado de São Paulo 2. ed. Campinas: Instituto Agronômico de Campinas, 1996. 285 p. (Boletim Técnico, 100).
  • SIPOS, P. et al Sorption of copper, zinc and lead on soil mineral phases. Chemosphere, v. 73, p. 461-469, 2008.
  • SPOSITO, G. The chemistry of soils 2.ed. New York: Oxford University Press, 2008. 342p.

Datas de Publicação

  • Publicação nesta coleção
    Mar 2016

Histórico

  • Recebido
    30 Set 2014
  • Aceito
    30 Set 2015
Universidade Federal do Ceará Av. Mister Hull, 2977 - Bloco 487, Campus do Pici, 60356-000 - Fortaleza - CE - Brasil, Tel.: (55 85) 3366-9702 / 3366-9732, Fax: (55 85) 3366-9417 - Fortaleza - CE - Brazil
E-mail: ccarev@ufc.br