Proposição de valores de referência para a concentração natural de metais pesados em solos brasileiros

Proposition of reference values for natural concentration of heavy metals in Brazilian soils

GESTÃO E CONTROLE AMBIENTAL

Proposição de valores de referência para a concentração natural de metais pesados em solos brasileiros

Proposition of reference values for natural concentration of heavy metals in Brazilian soils

Francisco de S. Fadigas^I; Nelson M. B. do Amaral Sobrinho^II; Nelson Mazur^II; Lúcia H. C. dos Anjos^II; Alessandra A. Freixo^III

^IDepartamento de Química Agrícola e Solos - UFBA, CEP.44 380-000, Cruz das Almas, BA. Fone: (75) 26211220. E-mail: [e-mail]

^IIDepartamento de Solos - UFRRJ -, BR 465, Km 7, CEP.23850 000, Seropédica, RJ. Fone (21) 37873755. E-mail: [e-mail]

^IIIDepartamento de Botânica - UEFS, Av. Universitária, s/n - Km 03 da BR 116, CEP: 44031-460 Feira de Santana, BA. Fone (75) 2248200. E-mail: [e-mail]

RESUMO

Para avaliar a extensão da poluição de uma área, é comum se comparar os teores totais de metais pesados encontrados num dado solo com aqueles defrontados em condições naturais (não poluídos) ou com valores de referência (padrões). Propuseram-se através deste estudo: 1) valores de referência (VR) e limites de tolerância (LT) para os mais representativos solos brasileiros, e 2) um modelo para obtenção das concentrações naturais de Cd, Co, Cr, Cu, Ni, Pb e Zn em um solo, a partir dos teores de silte, argila, Mn, Fe e CTC. Reuniram-se 256 amostras de solo por similaridade, em sete grupos e se calcularam as concentrações correspondentes ao quartil superior dos dados. Essas concentrações são propostas como valores de referência (VR) em solos brasileiros; além disso, obtiveram-se os limites de tolerância (LT) para cada grupo a partir da expressão antilog (m + 2s), cujos m (média) e s (desvio padrão) dos dados foram transformados em log₁₀. As funções de classificação da análise discriminante mostraram se apropriadas para alocar novas amostras nos grupos estabelecidos. Assim, é possível avaliar solos sob atividade antrópica e saber se há motivo para suspeitar de que a área esteja poluída.

Palavras-chave: poluição do solo, teores naturais, níveis de metais, metais em solos

ABSTRACT

To evaluate the extension of soil pollution of a site it is common to compare the contents found in a specific soil with those in natural conditions (unpolluted) or with reference values (patterns). The goal of this study was to propose: 1) reference values (VR) and tolerance limits for the most representative Brazilian soils, and 2) a model for obtaining the levels of Cd, Co, Cr, Cu, Ni, Pb and Zn in a soil under natural conditions from the silt, clay, Mn, Fe contents and CEC value. A set of 256 soil samples was formed by similarity in seven groups and the concentrations corresponding to the upper quartile of the data were calculated. These concentrations are proposed as the reference values (VR) in Brazilian soils. Besides, tolerance limits admitted for each group were taken as antilog (m + 2s), where m (average) and s (standard deviation) were transformed to log₁₀. The classification functions of the discriminant analysis proved to be appropriate to allocate new samples in the established groups. Thus, it is possible to evaluate soils under anthropogenic activity and to know if there is a reason for suspecting of pollution of the site.

Key words: soil pollution, baselines levels, metal levels, metals in soil

INTRODUÇÃO

Para se avaliar a extensão da poluição de uma área, é comum comparar os teores totais de metais pesados encontrados em um solo com aqueles defrontados em condições naturais (não poluídos) ou com valores de referência (padrões); todavia, no Brasil ainda não foi elaborado um sistema contendo padrões de referência de metais pesados no solo, para avaliação de áreas em relação à contaminação. Diferentes metodologias têm sido utilizadas para se estabelecer os valores de referência (VR) de metais pesados em solos, em que as mais comuns, são: o uso de valores considerados normais, citados na literatura; o estabelecimento de faixas de referência obtidas a partir de amostras controle, tomadas em áreas sem atividade antrópica, o mais longe possível de fontes de poluição (Casarini, 2000; Crock et al., 1992) e pela separação entre valores normais e anômalos, obtidos de uma coleção de dados que inclui solos contaminados e não contaminados (Tobias et al., 1997a; Wang, 1994).

O estabelecimento de faixas de concentração de referência e de limites máximos considerados normais no solo, tem sido feito por meio de diferentes procedimentos estatísticos (Fadigas et al., 2002). A expressão L - antilog (m ± 2s), em que m é a média e s é o desvio padrão dos dados transformados em log₁₀, foi utilizada para calcular os limites superior e inferior dos valores de referência em solos do Alaska (Crock et al., 1992) e por Tobias et al. (1997a) para estimar faixas de concentração normal de alguns elementos traço no nordeste da Espanha. Valores de referência para o estado de São Paulo foram definidos com base no quartil superior (75%) dos resultados analíticos (Lemos, 2000) e em solos da China, as concentrações de metais pesados em argilito foram usadas como padrões gerais (Wang, 1994). Técnicas da estatística multivariada como análise de agrupamento e análise discriminante foram utilizadas em estudos geoquímicos (Moura, 1985; Oliveira et al., 1998) e de inter relações entre atributos do solo e produtividade de cultivos (Correia, 1993), as quais apresentam potencial para uso no estabelecimento e proposição de valores de referência de metais em solos.

Alguns atributos do solo se destacam por sua relação com o conteúdo de metais pesados. Diversos autores encontraram correlações positivas entre o teor de argila e ferro e os teores totais de metais (Pérez et al., 1997; Rovers et al., 1983; Santos Filho & Rocha, 1982; Souza et al., 1996; Valadares & Catani, 1975; Valadares, 1975). Também pode haver associação entre metais pesados e a fração silte, visto que nesta fração do solo podem ser encontradas partículas finas (argila) e óxidos, devido à ineficiência no processo de dispersão dos solos durante a análise granulométrica (Klamt & van Reeuwijk, 2000; Lima, 1984); por outro lado, nem sempre tem sido obtida correlação significativa com a matéria orgânica e o pH do solo (Furlani et al., 1977; Lavado et al., 1998; Qian et al., 1996).

Objetivou-se com este estudo, propor: 1) valores de referência (VR) e limites de tolerância (LT) para os mais representativos solos brasileiros e 2) um modelo para obtenção dos teores naturais de Cd, Co, Cr, Cu, Ni, Pb e Zn em um solo, a partir dos teores de silte, argila, Mn, Fe e a CTC.

MATERIAL E MÉTODOS

Para este estudo, utilizou-se um conjunto de 256 amostras de solo, correspondendo aos horizontes A e B (ou C) de 110 perfis, das principais classes de solos encontradas no Brasil, especialmente os Argissolos e Latossolos (Tabela 1). Os solos selecionados estavam sob mínima ou nenhuma atividade antrópica. As análises de metais pesados foram feitas por extração com água régia e as determinações em Plasma de Acoplamento Indutivo (ICP-AES).

Os resultados obtidos (médias das repetições) foram submetidos à análise de agrupamento (AC), utilizando-se os teores de silte, argila, Mn e Fe (pseudo-total) e CTC como variáveis, considerando-se que estas possuem boa relação com o teor dos metais em estudo (Fadigas, 2002). Estas variáveis apresentam distribuição assimétrica e foram transformadas em log₁₀, embora a ausência de distribuição normal não seja considerada crítica em análise multivariada, quando se trabalha com grande número de amostras ou repetições (Statsoft, 1999). Para a análise de agrupamento, os dados transformados foram padronizados (valor observado menos a média/desvio padrão), para evitar a interferência das diferentes unidades de medida nos cálculos (Broschat,1979; Parks, 1966).

A escolha das variáveis foi baseada no estudo de correlações simples (Pearson) e no estudo de componentes principais, segundo Fadigas (2002). A análise foi processada em modo Q, ou seja, com o uso de um coeficiente de distância, para medir a similaridade entre as amostras, após a padronização dos dados. Para o agrupamento adotou-se a distância Euclidiana (distância geométrica entre dois objetos no espaço multidimensional) como medida de dissimilaridade. Embora o uso da distância Euclidiana possa causar distorção na estimativa da distância entre objetos quando são consideradas muitas variáveis (Cruz & Regazzi, 2001), sua escolha neste estudo se justifica pela impossibilidade de se utilizar outras medidas, como a distância Euclidiana média ou à distância de Mahalanobis para um conjunto de 256 amostras nos softwares em que foram realizadas as análises estatísticas (SAEG 5.0, Statistica versão 5.5 e SPSS 10.0). Os casos foram ligados entre si pelo método de Ward (Statsoft, 1999), no qual se procura diminuir a soma de quadrados (SQ) dos grupos formados em cada etapa e é considerado muito eficiente.

Após o estabelecimento dos grupos utilizou-se a análise discriminante (AD) com o objetivo de avaliar o grau de acerto na alocação das amostras nos sete grupos e obter as funções de classificação que possibilitassem o enquadramento de novas amostras. Na análise discriminante, os dados não foram padronizados, pois nela a diferença entre unidades de medida não interfere nos resultados. Para a classificação das amostras empregou-se a distância generalizada de Mahalanobis (Statsoft, 1999) e a probabilidade proporcional ao tamanho do grupo. As análises foram processadas através do programa Statistica versão 5.5 (Statsoft, 1999). Os limites para os teores metais (Cd, Co, Cr, Cu, Ni, Pb) considerados normais no solo em cada grupo, foram calculados adotando se o quartil superior (QS), ou seja, o valor correspondente a 75% da distribuição da freqüência acumulada (Casarini, 2000) e os limites máximos obtidos pela expressão LT antilog (m + 2s), em que LT - limite de tolerância, m - média e s - desvio padrão dos valores em log₁₀, respectivamente (Tobias et al., 1997a, 1997b).

RESULTADOS E DISCUSSÃO

Obtenção de grupos de amostras por semelhança

O estabelecimento de VR a partir de resultados analíticos de concentração de metais pesados, seja com base na média, intervalos de confiança ou quartil superior tem, como principal limitação, a diversidade dos solos encontrados no Brasil; assim, sua definição a partir das médias ou intervalos globais de concentração, não possuiria maior significado, uma vez que existe grande variabilidade desses teores entre os diferentes tipos de solo e, por conseqüência, seria obtida uma ampla faixa de variação. Da mesma forma o estabelecimento de um único limite máximo para concentrações normais baseado no quartil superior, não seria adequado para avaliação de solos cujas concentrações de metais fossem naturalmente baixas, indicando a necessidade de se utilizar uma metodologia que reúna solos com características semelhantes, especialmente aqueles atributos correlacionados com o teor de metais pesados. Para contornar este problema, as amostras de solo foram classificadas com base no seu grau de semelhança, em grupos homogêneos, levando se em conta a contribuição conjunta das variáveis silte, argila, Mn, Fe e CTC (Figura 1).

A interpretação do dendrograma considerando se um valor de 20 para o coeficiente de distância (Figura 1), permitiu que se identificasse sete grupos (G1 a G7) de amostras mais próximas entre si. As delimitações foram estabelecidas com base na análise visual do dendrograma, na qual se avaliaram os pontos de alta mudança de nível (Cruz & Regazzi, 2001; Moura, 1985). Analisando-se a composição dos grupos, constatou-se que todos se compunham tanto por amostras do horizonte A quanto de camadas intermediárias (A/B, B/A ou E) ou do B, em proporções variadas, indicando a coerência de se utilizar os atributos do solo que têm relação com o teor de metais pesados na formação dos grupos, sem considerar a qual horizonte do solo cada amostra pertence.

Os grupos obtidos por semelhança entre as amostras também se compõem de solos de diferentes classes pedológicas, uma vez que solos de classes pedológicas diferentes podem apresentar alguns atributos químicos e de composição granulométrica em comum (EMBRAPA, 1999); assim, amostras de horizontes A, B ou transicionais de um mesmo solo ou de solos distintos, podem estar presentes em um mesmo grupo, desde que as concentrações de algumas das variáveis consideradas na análise de agrupamento, sejam próximas (Tabela 2). Alguns autores observaram que o teor de metais também varia entre solos pertencentes a uma mesma classe, ou entre solos de classes diferentes, em função das variações nas características dos solos (Amaral Sobrinho, 1993; Mazur, 1997; Ramalho et al., 2000; Valadares, 1975).

A análise das médias obtidas para as variáveis constantes na Tabela 2 permite que as variáveis que melhor caracterizam cada grupo, sejam identificadas, deste modo, pode se dizer que o grupo 1 (G1) é caracterizado pelas maiores concentrações de Mn e Fe, seguidos por um teor de argila próximo daqueles encontrados nos grupos 2 (G2) e 3 (G3). Este grupo (G1), é composto, predominantemente, por amostras de Latossolo Vermelho Distroférrico, Nitossolo Vermelho Distroférrico, Latossolo Bruno e Argissolo Vermelho Eutroférrico, formados sobre basalto ou materiais derivados (Oliveira, 1996), naturalmente ricos em metais pesados (Alloway, 1995; Ross, 1994; Oliveira et al., 1998; Valadares, 1975) fazem parte deste grupo, ainda, amostras de Latossolo Vermelho Distrófico, Argissolo Vermelho, Luvissolo e Cambissolo, originados de basalto, migmatito, biotita gnaisse e argilito da Formação Solimões. Os teores mais elevados de Mn e Fe neste grupo refletem a origem dos solos que os compõem.

O grupo 2 se diferencia dos demais, especialmente por seus teores de silte e Mn e pela CTC. Deste grupo participam diversos solos com CTC elevada, como Chernossolos, Luvissolos e solos eutróficos, como Planossolo Háplico Eutrófico, Plintossolo Argilúvico Eutrófico e Gleissolo Háplico Eutrófico; ressalta-se ser ele composto, ainda, por amostras de Latossolo Amarelo, Latossolo Vermelho Amarelo e Argissolo Vermelho.

Outro grupo que se destaca é o grupo 7 (G7) no qual são encontradas as menores concentrações médias das variáveis utilizadas para o agrupamento das amostras de solos, exceto a CTC. Em termos gerais, os grupos G6 e G7 são similares em composição, mas se diferenciam pelas concentrações de Argila e Fe; neles se encontram principalmente amostras de Latossolos e Argissolos Amarelos Distróficos, seguidas por pequeno número de Latossolos e Argissolos Vermelho Amarelos Distróficos e Planossolo Háplico Distrófico, originados predominantemente de sedimentos Terciários e Quaternários.

Os grupos G3, G4 e G5 apresentam características intermediárias entre G1 e G6/G7 e são constituídos por amostras de grande variedade de Latossolos e Argissolos e, em menor quantidade, por amostras de Plintossolo, Cambissolo e Nitossolo, com predominância acentuada do caráter distrófico.

O uso da análise de agrupamento possibilitou a separação de grupos mais ou menos uniformes e indica que, com base nas características dos solos presentes em cada grupo (silte, argila, CTC, Mn e Fe) e que possuem relação com a concentração dos metais pesados no solo, será possível obter-se as concentrações naturais de Cd, Co, Cr, Cu, Ni, Pb e Zn, para uma ampla variedade de solos brasileiros.

Valores propostos como referência

O agrupamento das amostras possibilitou se conhecer a distribuição das concentrações desses elementos em diferentes níveis, a partir do conhecimento dos teores de silte, argila, CTC, Mn e Fe, ou seja, de acordo com importantes atributos do solo e que estão associadas a processos geoquímicos e estudos pedológicos (Oliveira et al., 1992). Esta distinção pode ser melhor visualizada nos dados encontrados na Tabela 3, em que são mostrados os valores considerados normais em solos brasileiros, tomando-se por base o valor correspondente ao QS (75%) da distribuição de freqüências dos dados de cada grupo formado.

Analisando-se os dados da Tabela 3, é possível identificar, no G1, as maiores concentrações para quase todos os metais devido, principalmente, aos maiores teores de Mn e Fe, associados a uma quantidade expressiva de argila, o que, por sua vez, está relacionado com a origem basáltica de vários solos componentes deste grupo. O inverso ocorre no G7 em que a concentração de todos os elementos é inferior à dos demais grupos, exceto para o Pb. Os valores de Pb para este grupo devem ser vistos com reserva, pois além de não ter sido obtida leitura para 40% das amostras, a concentração média foi estimada com maior erro. As menores concentrações encontradas neste grupo são, possivelmente, decorrentes do tipo de material de origem associado aos solos deste grupo, predominantemente sedimentar.

Os grupos G6 e G7 são relativamente próximos em termos de concentrações de metais. As variações encontradas são, sem dúvida, devidas aos maiores teores de Fe e Argila no G6 (Tabela 2). Os demais grupos apresentam concentrações intermediárias entre o G1 e G7, havendo tendência geral para que as concentrações de todos os elementos diminuam a partir do G1 (Tabela 3).

Os valores médios do quartil superior (QSm) apresentados na Tabela 3, são próximos dos valores citados por Casarini (2000) como valores de referência de qualidade para metais em solos do estado de São Paulo, cujas concentrações em mg kg^-1, são: Cd (0,5), Co (12,5), Cu (35,1), Cr (40,2), Ni (13,2), Pb (17), Zn (59,9). Esses valores também foram calculados com base no QS e apenas o Zn apresenta uma concentração nitidamente superior à que foi obtida para a média do QS, neste estudo. As concentrações correspondentes ao QS podem ser consideradas valores normais para solos com características similares às dos sete grupos (Tabela 2) e são propostas para adoção como VR para concentrações naturais de solos brasileiros.

Além de servir para se obter os VR, os dados de concentração dos metais pesados também foram utilizados para o cálculo do limite máximo, a fim de que a concentração de um elemento fosse admitida como ainda esperada para determinado grupo (Tabela 4); em outras palavras, concentrações superiores a este limite seriam indicativas de que dada amostra não apresenta mais a concentração normal nem o valor máximo, através do qual seria possível se encontrar amostras semelhantes às do grupo em que se enquadra (teores de Silte, Argila, CTC, Mn e Fe). Esses limites foram chamados de limites de tolerância (LT) e correspondem ao valor obtido pela expressão antilog (m + 2s), em que m é média e s o desvio dos dados transformados em log₁₀. Abaixo deste valor estarão contidos cerca de 97,25% dos casos que podem ocorrer, considerando-se a variação dentro de cada grupo.

Como o objetivo deste trabalho é sugerir valores para caracterização das concentrações de metais em solos, que sejam úteis para interpretações do ponto de vista da poluição ambiental, somente o limite superior foi calculado. A expressão visando aos limites superior e inferior antilog (m ± 2s) foi utilizada por outros autores para o cálculo de concentrações normais de metais pesados nos solos (Crock et al., 1992; Reaves & Berrow, 1984; Tobias et al., 1997a)

Avaliação do grau de acerto na classificação das amostras por agrupamento

Uma maneira de se verificar a qualidade da distribuição das amostras resultante da formação dos grupos, é avaliar o grau de acerto na classificação com a AD. A AD é usada para interpretar grupos de objetos, definidos, a priori, pelos métodos de agrupamento e ordenação (Valentin, 2000). Esta técnica não consiste em formar grupos, já que eles são previamente conhecidos, mas interpretá los com base em variáveis que possam discriminar dois ou mais grupos que ocorram naturalmente (Statsoft, 1999). No presente estudo, como não havia um conhecimento prévio dos grupos que provavelmente se formariam, admitiu-se que os grupos obtidos por semelhança são aqueles que ocorreriam naturalmente, representando ambientes edáficos, comuns no Brasil; como resultado, obteve-se a matriz de classificação para as amostras em cada grupo e a percentagem de acertos em relação à classificação anteriormente obtida na análise de agrupamento (Tabela 5).

Em média, 89% das amostras foram agrupados em conjuntos que, de fato, são distintos entre si, em função dos teores de silte, argila, CTC, Fe e Mn; entretanto, em média 11% poderiam estar compondo grupos diferentes daqueles em que foram inicialmente alocados. Embora a percentagem de acertos não tenha sido a ideal, pode-se considerar que a formação dos grupos possui consistência suficiente para a validação dos valores obtidos para a concentração normal dos metais pesados em cada grupo.

Além de possibilitar a avaliação do grau de acerto na distribuição das amostras entre os grupos, a análise discriminante também permite alocar-se novos casos nos grupos já formados. As funções de classificação podem ser usadas para determinar a que grupo uma nova amostra tem maior probabilidade de pertencer e, em geral, os novos casos são classificados no grupo para o qual a função de classificação correspondente resulta no maior valor (Statsoft, 1999). Este procedimento pode ser de grande utilidade na avaliação de áreas suspeitas de estarem poluídas e para o monitoramento de áreas onde se desenvolvem atividades agrícolas ou industriais. Amostras coletadas em áreas cuja concentração natural já não pode ser determinada, podem ser enquadradas em um dos grupos formados (Tabela 2) e, em seguida, pode ser obtido o VR no solo, para cada elemento (Tabela 3), como também o LT (Tabela 4). A comparação direta entre o valor de concentração obtido pela análise da nova amostra e o VR e LT, indicaria se há ou não motivos para se suspeitar que a área em questão contém concentrações do elemento além daquela que normalmente teria, em função das variáveis de solo consideradas para a alocação da amostra no grupo. Na Tabela 6 são apresentadas as funções de classificação para os sete grupos.

Teste para verificar a viabilidade de uso das funções de classificação

Na tentativa de se verificar se realmente o uso das funções é viável na prática, tomaram-se dez amostras de solo, que não foram utilizadas na formação dos grupos e cujos atributos considerados são mostrados na Tabela 7.

Os valores para cada função, resultantes da substituição das variáveis (silte, argila, CTC, Mn e Fe) pela concentração de cada uma nas dez amostras, foram comparados entre si para determinar em qual grupo estas deveriam ser alocadas, o que pode ser visto na Tabela 8. O maior valor obtido entre as funções, para cada amostra, indica de qual grupo ela está mais próxima (Statsoft, 1999). As concentrações de metais pesados determinadas por extração com água régia, das dez amostras, foram então comparadas com os VR's e LT's constantes nas Tabelas 3 e 4. Observou-se que a concentração determinada pela análise química foi 100% dos casos, inferior ao LT e entre 40 e 70% dos casos o valor determinado foi < VR; os resultados dão uma boa indicação do potencial de uso das funções de classificação para a avaliação do teor de metais em amostras de solo provenientes de áreas sujeitas a atividades antrópicas.

CONCLUSÕES

1. As concentrações de Cd, Co, Cr, Cu, Ni, Pb e Zn definidas em cada grupo poderiam ser utilizadas para a avaliação preliminar de áreas com suspeita de poluição ou com uso agrícola, cujos solos possuam características similares às deste estudo.

2. As funções de classificação mostraram ser uma ferramenta adequada para a alocação de novas amostras, nos grupos previamente estabelecidos.

LITERATURA CITADA

Protocolo 149.04 - 07/10/2004 - Aprovado em 02/01/2006

• Alloway, B. J. The origin of heavy metals in soils. In: Alloway, B. J. (ed.) Heavy metals in soils. 2.ed. Glasgow: Blackie Academic & Professional, 1995. p.38 57.
• Amaral Sobrinho, N. M. B. Interação dos metais pesados de resíduos siderúrgicos com um podzólico Vermelho Amarelo, 1993. Viçosa UFV. 163p. Tese Doutorado
• Broschat, T. K. Principal component analysis in horticultural research. Hortscience, Alexandria, v.14, n.2, p.114 117, 1979.
• Casarini, D. C. P. Proposta de valores de referência de qualidade e intervenção para solos e águas subterrâneas no estado de São Paulo. In: Seminário Internacional sobre Qualidade de Solos e Águas Subterrâneas, 2., 2000, São Paulo. Anais. São Paulo: CETESB, 2000. 165p.
• Correia, J. R. Uso de técnicas multivariadas no estudo das inter relações de características do solo e a produtividade do eucalipto.. Viçosa, UFV, 1993, 110p. Dissertação Mestrado
• Crock, J. G.; Severson, R. C.; Gough, L. P. Determining baselines and variability of elements in plants and soils near the Kenai National Wildlife Refuge, Water, Air, and Soil Pollution, Netherlands, v.63, n.3/4, p.253-271, 1992.
• Cruz, C. D.; Regazzi, A. J. Modelos biométricos aplicados ao melhoramento genético. Viçosa: UFV, 2001. 390p.
• EMBRAPA - Empresa Brasileira de Pesquisa Agropecuária.Centro Nacional de Pesquisa de Solos. Sistema brasileiro de classificação de solos. Brasília: Embrapa Produção de Informação, 1999. 412p.
• Fadigas, F. S. Estimativa das concentrações naturais (pseudo total) de Cd, Co, Cr, Cu, Ni, Pb e Zn em solos brasileiros e proposição de valores de referência utilizando técnicas da estatística multivariada. Seropédica. UFRRJ, 2002. 116p. Tese Doutorado
• Fadigas, F. S.; Amaral Sobrinho, N. M. B.; Mazur, N.; Anjos, L. H. C.; Freixo, A. A. Concentrações naturais de metais pesados em algumas classes de solos brasileiros. Bragantia, Campinas, v.2, p.151-159, 2002.
• Furlani, P. R.; Bataglia, O. A.; Valadares, J. M. A. Cobalto em solos do Estado de São Paulo. Revista Brasileira de Ciência do Solo, Piracicaba, v.1, p.65-67, 1977.
• Klamt, E.; van Reeuwijk, L. P. Evaluation of morphological, physical and chemical characteristics of ferralsols and related soils. Revista Brasileira de Ciência do Solo, Viçosa, v.24, n.3, p.573 587, 2000.
• Lavado, R. S.; Rodriguez, M. B.;Sheiner, I. D.; Taboada, M. A.; Rubio, G.; Alvarez. R.; Alconada, M.; Zubillaga, M. S. Heavy metals in soils of Argentina: comparison between urban and agricultural soils. Communications in Soil Science and Plant Analysis, Athens, v.29, n.11 14, p.1913-1917, 1998.
• Lemos, M. M. G. Metodologia adotada para o estabelecimento dos valores de referência de qualidade para solos e águas subterrâneas no estado de São Paulo. In: CETESB. Prevenção e controle da poluição do solo e das águas subterrâneas. São Paulo: CETESB, 2000. p.68-77.
• Lima, P. C. Micromorfolgia de horizontes B de Latossolos do Sudeste e Sul do Brasil. In: EMBRAPA - Empresa Brasileira de Pesquisa Agropecuária. Investigações coligadas em Latossolos do Sudeste e Sul do Brasil. Rio de Janeiro: EMBRAPA-SNCLS, 1984. p.168-188.
• Mazur, n. Níquel, chumbo, zinco e cobre em solos que receberam composto de resíduos sólidos urbanos. Viçosa, UFV,1997 135p. Tese Doutorado
• Moura, C. V. A. Aplicação de tratamento estatístico multivariante em dados geoquímicos de solo no mapeamento geológico na província de Carajás (alvo 2 - corpo 4). Revista Brasileira de Geociências, São Paulo, v.15, n.3, p.241 248, 1985.
• Oliveira, J. B.; Jacomine, P. K. T.; Camargo, M. N. Classes gerais de solos do Brasil. 2.ed. Jaboticabal: FUNEP, 1992. 201p.
• Oliveira, T. S. Metais pesados como indicadores de materiais de origem de solos.. Viçosa, UFV, 1996, 128p. Tese Doutorado
• Oliveira, T. S.; Costa, L. M.; Cruz, C. D. Importância relativa dos metais pesados do solo na identificação e separação de materiais de origem. Revista Ceres, Viçosa, v.45, n.260, p.359 371, 1998.
• Parks, J. M. Cluster analysis applied to multivariate geologic problems. Journal of Geology, Chicago, v.74, n.5, p.703 715, 1966.
• Pérez, D. V.; Saldanha, M. F. C.; Meneguelli, N. A.; Moreira, J. C; Vaitsman, D. S. Geoquímica de alguns solos brasileiros. Pesquisa em andamento. Rio de Janeiro, n.4, p.114, 1997.
• Qian, J.; Wang, Z.; Shan, X.; Tu Q.; Wen B.; Chen, B. Evaluation of plant availability of soil trace metals by chemical fractionation and multiple regression analysis. Environmental Pollution, Amsterdam, v.91, n.3, p.309 315, 1996.
• Ramalho, J. F. P. G.; Amaral Sobrinho, N. M. B.; Velloso, A. C. X. Acúmulo de metais pesados em solos cultivados com cana-de-açucar pelo uso contínuo de adubação fosfatada e água de irrigação. Revista Brasileira de Ciência do Solo, Viçosa, v.23, n.4, p.971 979, 1999.
• Ramalho, J. F. P. G.; Amaral Sobrinho, N. M. B.; Velloso, A. C. X. Contaminação da microbacia de Caetés com metais pesados pelo uso de agroquímicos. Pesquisa Agropecuária Brasileira, Brasília, v.35, p.1289 1303, 2000.
• Reaves, G. A.; Berrow, M. L. Total lead concentration in Scottish soils. Geoderma, Amsterdam, v.32, p.1 8, 1984.
• Ross, S. M. Toxic metals in soil-plant systems. Chichester: John Willey & Sons, 1994. 469p.
• Rovers, H.; Camargo, O. A.; Valadares, J. M. A. S. Níquel total e solúvel em DTPA em solos do Estado de São Paulo. Revista Brasileira de Ciência do Solo, Piracicaba, v.7, n.3, p.217 220, 1983.
• Santos Filho, A.; Rocha, H. O. O cobre em solos do Estado do Paraná. Revista do Setor de Ciências Agrárias, Curitiba, v.4, n.1 2, p.23-26, 1982.
• Souza, M. L. P.; Andreoli, C. V.; Amaral, M. B., et al. Levantamento preliminar dos teores de metais pesados em alguns solos do Paraná. Revista Sanare, Curitiba, v.5, n.5, p.68-75, 1996.
• Statsoft, Inc. Statistica for Windows - computer program manual. Tulsa (OK): StatSoft, 1999.
• Tobias, F. J.; Bech, J.; Algarra, P. S. Establishment of the background levels of some trace elements in soils of NE Spain with probability plots. The Science of Total Environment, Amsterdam,v.206, p.255 265, 1997a.
• Tobias, F. J.; Bech, J.; Algarra, P. S. Statistical approach to discriminate background and anthropogenic input of trace elements in soils of Catalonia, Spain. Water, Air, and Soil Pollution, Netherlands, v.100, p.63-78, 1997b.
• Valadares, j. m. a. s. Cobre em solos do Estado de São Paulo: cobre total. Bragantia, Campinas, v.34, n.4, p.125-132, 1975.
• Valadares, j. m. a. s.; Catani, r. a. Zinco em solos do Estado de São Paulo: Zinco total. Bragantia, Campinas, v.34, p.134 139, 1975.
• Valentin, J. L. Ecologia numérica - uma introdução à análise multivariada de dados ecológicos. Rio de Janeiro: Interciência, 2000. 117p.
• Wang, H. K. Research methodology for determining background levels of trace elements in Chinese soils. Environmental Geochemistry and Health, 1(Special Issue), Amsterdam, p.397-413, 1994.

Datas de Publicação

• Publicação nesta coleção
  06 Fev 2007
• Data do Fascículo
  Set 2006

Histórico

• Recebido
  07 Out 2004
• Aceito
  02 Jan 2006