Lixiviação de naftaleno em solos urbanos da região metropolitana do Recife, PE

Leaching transport of naphthalene in urban soils from the metropolitan region of Recife, PE

DIVISÃO 3 - USO E MANEJO DO SOLO

COMISSÃO 3.5 - POLUIÇÃO, REMEDIAÇÃO DO SOLO E RECUPERAÇÃO DE ÁREAS DEGRADADAS

Lixiviação de naftaleno em solos urbanos da região metropolitana do Recife, PE¹

Leaching transport of naphthalene in urban soils from the metropolitan region of Recife, PE

André Inácio do Carmo^I; Antônio Celso Dantas Antonino^II; Jean Manuel Fonseca Martins^III; Valdinete Lins da Silva^IV; Marie Christine Morel^V; Jean Paul Gaudet^VI

^IDoutor pelo Programa de Pós-Graduação em Tecnologias Energéticas e Nucleares, PROTEN-UFPE. E-mail: [e-mail]

^IIProfessor Associado, Departamento de Energia Nuclear, UFPE. Av. Prof. Luiz Freire 1000, Cid. Universitária. CEP 50740-540 Recife (PE). Bolsista do CNPq. E-mail: [e-mail]

^IIIPesquisador CNRS, Laboratoire d´étude des Transferts en Hydrologie et Environnement, LTHE-CNRS, Université Grenoble I, Cedex 09, França. E-mail: [e-mail]

^IVProfessora Titular, Departamento de Engenharia Química, UFPE. Bolsista do CNPq. E-mail: [e-mail]

^VProfessora Associada CNAM Paris, STI Department, 292 Rue Saint Martin 75141 Paris Cedex 03 e Université Grenoble 1, Laboratoire d'étude des Transferts en Hydrologie et Environnement, LTHE-INPG, Grenoble 1, Cedex 09, França. E-mail: [e-mail]

^VIPesquisador Université Grenoble, Laboratoire d´étude des Transferts en Hydrologie et Environnement, LTHE-CNRS, Université Grenoble I, Cedex 09, França. E-mail: [e-mail]

RESUMO

Os principais atributos edáficos relacionados com a lixiviação de naftaleno nos solos consistem na granulometria, no teor de matéria orgânica e na composição mineralógica das frações silte e argila. Este trabalho objetivou determinar os valores da dispersão e do retardamento do naftaleno em amostras de solos coletadas nas proximidades de tanques de armazenamento de combustíveis localizados em Recife, PE. Os experimentos de lixiviação foram realizados em colunas contendo amostras deformadas e foram modelados usando a equação de convecção-dispersão, que foi ajustada aos dados experimentais com o emprego do software CXTFIT 2.0. Os resultados indicaram que solos arenosos oferecem menor resistência à lixiviação de naftaleno e que a presença de goethita, juntamente com a matéria orgânica, é possivelmente o principal fator responsável pelo retardamento do transporte desse composto no conjunto de amostras estudadas.

Termos de indexação: parâmetros hidrodispersivos, curva de lixiviação, modelo CDE, fator de retardamento, contaminação ambiental.

ABSTRACT

The main soil properties related to the leaching of naphthalene in soils are particle size, organic matter content and the mineralogical composition of silt and clay fractions. The aim of this study was to evaluate naphtalene reactive transfer in soil columns determinated by the dispersion and retardation of naphthalene in soil samples collected from the surroundings of fuel storage tanks in Recife, Brazil. Leaching experiments were carried out in columns containing disturbed soils and were modeled using the convection-dispersion equation, which was adjusted to the experimental data with the CXTFIT software 2.0. The results indicated that sandy soils offer less resistance to leaching of naphthalene and that the presence of goethite, together with organic matter appeared to be the main factors responsible for the retardation of naphtalene transport in the studied samples.

Index terms: hydrodispersive parameters, leaching curves, model CDE, retardation factor, environmental contamination.

INTRODUÇÃO

A contaminação dos solos por derivados de petróleo tem sido sério problema em países em desenvolvimento. A escala de perigo imposta no ambiente natural por esses contaminantes depende da superfície da área contaminada, composição química dos hidrocarbonetos e profundidade em que ocorre a contaminação (Wolicka et al., 2009).

Dentre os hidrocarbonetos potencialmente perigosos para a saúde humana, destaca-se o grupo dos hidrocarbonetos aromáticos policíclicos (HPAs), que são compostos constituídos por dois ou mais anéis aromáticos, potencialmente cancerígenos e altamente resistentes à biodegradação, permanecendo, assim, por longos períodos no ambiente (Yan et al., 2009; Hollosi et al., 2011; Ruiz et al., 2011).

Os HPAs são introduzidos no ambiente por meio de várias rotas e estão presentes em vários ambientes, como atmosfera (Qi et al., 2001), água (Zhou & Maskaoui, 2003), solos (Mielke et al., 2001) e sedimentos (McCready et al., 2000). Esses compostos são derivados principalmente da combustão incompleta de combustíveis fósseis, da queima da vegetação e de outros materiais orgânicos.

Nos centros urbanos, grandes quantidades de HPAs solúveis em água, como o naftaleno, podem ser transportados em direção ao lençol freático, promovendo a contaminação e deterioração dos aquíferos subterrâneos (Smith et al., 2003). Ainda são poucas as pesquisas sobre o transporte desse contaminante em solos urbanos (Lee et al., 2002; Motelay-Massei et al., 2004; Schwarz et al., 2011), sendo este estudo de fundamental importância para o entendimento dos processos e dos mecanismos físicos envolvidos no transporte do naftaleno nesses solos, com o intuito de prevenir a contaminação dos aquíferos subterrâneos existentes.

O naftaleno é um hidrocarboneto aromático policíclico, presente na composição química dos derivados de petróleo, especialmente na gasolina e no diesel. É utilizado comercialmente como matéria-prima de produtos utilizados na fabricação de cloreto de polivinila (PVC) e inseticidas. As características químicas do naftaleno são: massa molecular de 128 g mol^-1, fórmula molecular C₁₀H₈, ponto de fusão 80,35 °C, ponto de ebulição 218 °C, densidade de 1,16 g mL^-1 e solubilidade em água de 31,7 mg L^-1.

Este trabalho objetivou avaliar o transporte do naftaleno e HPA de maior solubilidade em água, em três solos contrastantes, localizados próximos aos tanques de armazenamento de combustíveis situados em Recife, PE.

MATERIAL E MÉTODOS

Solos

As amostras de solos foram coletadas em dois sítios na cidade de Recife, PE, que, segundo a classificação de Köppen, apresenta clima tropical. A precipitação pluvial anual local de 1.500 mm é concentrada principalmente no outono e inverno e a temperatura média anual é de 25 ^oC (Andrade, 2007).

O sítio 1 localiza-se ao lado do posto de combustível da Universidade Federal de Pernambuco - UFPE (8^o 03' 20.42" S e 34^o 57' 19.98" O), e o sítio 2, ao lado de um posto de combustível da Avenida General Polidoro, em Recife, PE (8^o 02' 31.14" S e 34^o 56' 36.86" O), sendo encontrados nesse sítio duas camadas (A e B).

As amostras de solo foram coletadas na profundidade de 0-20 cm para os sítios 1 e 2 (camada A) e na profundidade de 20-80 cm para o sítio 2 (camada B), secas ao ar, destorroadas e passadas em peneira de 2 mm. Em seguida, procedeu-se à caracterização granulométrica dessas, empregando-se o método da pipeta (Loveland & Whalley, 1991). As análises químicas e eletroquímicas consistiram nas determinações quantitativas de pH em água e em KCl, condutividade elétrica do extrato de saturação, carbono oxidável, Ca²⁺, Mg²⁺, K⁺, Na⁺, PO₄^3- e H+Al, a partir dos quais foram calculadas: soma de bases, capacidade de troca de cátions (CTC potencial), percentual de saturação por bases e percentual de sódio trocável (Embrapa, 1997).

A composição mineralógica da fração argila foi avaliada por meio de um difratômetro de raio-X Shimadzu XRD 6000, operando a uma tensão de 40 kv, com corrente de 20 mA, radiação Kα, anodo de cobre (Cu) e com monocromador de grafite. As amostras foram preparadas e analisadas na forma de pó (não orientadas), empregando um intervalo de varredura de 5 a 70^o 2θ e velocidade de 1,5^o 2θ min^-1. As amostras de argila foram saturadas tanto com KCl a 25 ^oC e aquecidas a 550 ^oC (amostras K25 e K550) quanto com MgCl₂ (amostra Mg). As argilas foram também analisadas na forma de esfregaço sobre lâminas de vidro (amostras orientadas), com amplitude de varredura entre 3 e 35^o 2θ e velocidade de registro de 1,5^o 2θ min^-1 (Jackson et al., 1986).

Procedimento experimental

Os ensaios de transporte foram realizados seguindo o protocolo experimental adotado por Carmo et al. (2010). O dispositivo experimental foi composto por colunas de solo em vidro com 5 cm de comprimento e 2,5 cm de diâmetro interno; uma bomba peristáltica; um coletor de frações; um condutivímetro digital para leitura das medidas de condutividade elétrica das soluções de KBr no efluente; capilares de teflon flexíveis com 2 mm de diâmetro interno; e balanças digitais para determinação da massa de solução deslocadora aplicada nas colunas de solo.

O acondicionamento do solo na coluna foi feito em camadas de aproximadamente 1 cm, levemente compactadas até que atingisse a massa específica aparente encontrada no campo para cada solo. Após a montagem, as colunas foram saturadas, no sentido ascendente delas, com solução iônica de CaCl₂ na concentração de 5,5 g L^-1, próxima à da solução do solo, para que os coloides das amostras de solo não sofressem desestabilização, comprometendo a permeabilidade em razão da diminuição da força iônica (Johnson et al., 1996). O volume da solução de CaCl₂ utilizado para saturação das colunas, até a formação de uma fina lâmina de água no topo dessas, foi utilizado como volume de poros (V₀).

Os ensaios consistiram em deslocar, por meio de uma bomba peristáltica conectada à parte superior da coluna, certo volume de líquido V_0, que ocupou o espaço poroso contido em uma coluna de solo, por meio de uma solução contendo o soluto (traçador ou soluto interativo) de concentração C₀, a uma velocidade aparente média v, tendo os efluentes da solução coletados na base da coluna, por meio de um coletor de frações.

O soluto tende a se difundir ao mesmo tempo em que se infiltra a velocidades variáveis, por meio dos poros do solo, originando a formação de uma zona de mistura característica do estado de dispersão do soluto.

Monitorou-se a progressão do avanço do soluto, medindo-se a concentração C do efluente no curso do tempo. A evolução da razão C/C₀ em razão do número de volumes de poros do efluente coletado (V/V₀) forneceu a curva de eluição do soluto.

Nos ensaios com KBr para definição dos parâmetros hidrodinâmicos e hidrodispersivos dos três solos, foram utilizadas as vazões de 0,2; 0,45; e 0,7 cm³ min^-1, enquanto nos ensaios para determinação da reatividade do naftaleno nesses solos, utilizaram-se apenas a vazão de 0,45 cm³ min^-1 (correspondente a uma precipitação pluvial média de 55 mm h^-1 na Região Metropolitana do Recife, PE) e a concentração de 20 mg L^-1 de naftaleno, correspondente a aproximadamente 63 % da solubilidade máxima em água para evitar erros por precipitação desse composto.

As amostras do efluente foram coletadas em tubos de vidro para evitar perdas por adsorção e os coletores de fração foram colocados dentro de um freezer a 0 ^oC para baixar a pressão do vapor do naftaleno, limitando com isso perdas por volatilização.

A solução-estoque de naftaleno foi preparada a partir da pesagem de 100 mg de naftaleno sólido, com 99,8 % de pureza, obtido da Sigma Aldrich (França), em uma balança de precisão, sendo em seguida diluída em 100 mL de metanol, em virtude da sua baixa solubilidade em água, obtendo-se uma solução final de 1 g L^-1. Após a dissolução em metanol, 20 mL da solução de naftaleno a 1 g L^-1 foi transferida para um balão de 1.000 mL, obtendo-se uma solução-estoque com concentração final de 20 mg L^-1.

A quantificação do naftaleno foi realizada por cromatografia líquida de alta eficiência (CLAE), utilizando o método descrito por Martins & Mermoud (1998). Foi utilizado um cromatógrafo, modelo Shimadzu 2000, no Laboratoire d'étude des Transferts en Hydrologie et Environnement (LTHE) na cidade de Grenoble, França. As amostras do efluente foram centrifugadas a 5.000 rpm por 10 min. Em seguida, realizou-se a análise quantitativa do sobrenadante (10 µL), usando uma coluna C18 (Supelco, 25 cm x 4,6 mm, 5µm), operando com a seguinte condição cromatográfica: 90 % metanol e 10 % água, vazão de 1 mL min^-1 e detecção a 267 nm.

Os valores de concentração relativa (C/C₀) e de seus respectivos valores de número de volume de poros foram submetidos ao software CXTFIT 2.0 (Parker & van Genuchten, 1984) para a resolução numérica do modelo de convecção-dispersão, cuja equação diferencial parcial é dada por:

em que C é a concentração do soluto expressa em massa de soluto por volume de solução [M L^-3]; z, a coordenada espacial [L]; t, o tempo [T]; P, o número de Peclet [1]; e R, o fator de retardamento [1].

As condições de contorno para a equação (1) foram:

- Condição inicial:

- Condição de contorno inferior:

- Condição de contorno superior:

em que C_o é a concentração do soluto na solução deslocadora (M L^-3); t, o tempo (T); e t_o, o tempo de aplicação da solução deslocadora [T].

Caracterização dos solos

No quadro 1, encontram-se as caracterizações granulométrica e mineralógica dos três solos, bem como a classificação textural desses, de acordo com Embrapa, 1997.

De acordo com os resultados da análise granulométrica do solo do sítio 1 (Quadro 1), constatou-se que esse classifica-se como Franco-arenoso (denominado de Franco-arenoso 1).

Os solos localizados no sítio 2 (camadas A e B) apresentaram grande diferença na textura, sendo a amostra da camada A classificada como Franco-arenoso (denominado de Franco-arenoso 2) e a da camada B, como Franco-argilossiltoso.

Pela diferença entre o pH em KCl e em água (Quadro 2), verificou-se que a carga líquida em todas as amostras é negativa, condizente com a mineralogia dos solos Franco-arenoso 1 e Franco-argilossiltoso, visto que a maioria dos minerais encontrada nessas amostras é do tipo 2:1 expansivo, como a esmectita e a vermiculita.

É importante ressaltar que essas amostras com minerais 2:1 expansivos, apesar de possuírem carga permanente negativa, apresentaram baixa capacidade de troca de cátions, como observado no quadro 2. Já o solo Franco-arenoso 2, por apresentar minerais do tipo 1:1, evidenciou baixa CTC advinda principalmente de cargas negativas de caráter variável.

Observou-se que o solo Franco-arenoso 2 apresenta comportamentos químico, físico e mineralógico diferenciados em relação aos demais, provavelmente por ter sido importado de outras áreas/regiões para fins de aterro, o que é constatado principalmente pela sua mineralogia diferenciada, em que prevalece principalmente caulinita e óxidos de ferro (goethita). Além disso, esse apresenta elevado percentual de sódio trocável (PST), que o classifica como solódico (6-15 % de sódio), diferentemente do solo Franco-argilossiltoso, que evidenciou valores baixos de PST, classificando-o como não sódico, sendo esse provavelmente o solo original dessa região. Essa última afirmação pode ser feita também para o solo Franco-arenoso 1.

Em relação à condutividade elétrica (CE), verificou-se que todos os solos apresentaram CE mediana e, conforme os valores encontrados de RAS, pode-se constatar que esse fato não foi em razão da presença do Na, mas provavelmente de outros cátions como H⁺ e Al³⁺, uma vez que, de forma geral, todos os solos apresentaram baixo somatório de bases (S) e consequentemente baixa CTC potencial (T).

Em relação ao carbono oxidável total, verificou-se que o solo Franco-arenoso 1 apresentou os maiores teores, principalmente em razão da deposição e decomposição de resíduos orgânicos. Como muitas estruturas da matéria orgânica não apresentam cargas elétricas, essas estruturas não iônicas podem reagir fortemente com partes não carregadas de compostos orgânicos, como o naftaleno, por meio de interações do tipo van der Waals (Sposito, 1989).

RESULTADOS E DISCUSSÃO

Caracterização hidrodinâmica com KBr

Condições experimentais

Os valores das variáveis determinadas experimentalmente para os ensaios de deslocamento miscível do KBr são apresentados no quadro 3.

As curvas médias de lixiviação do KBr ajustadas pelo modelo CDE, em colunas com os solos Franco-arenoso 1, Franco-arenoso 2 e Franco-argilossiltoso, saturadas nas vazões de 0,2; 0,45; e 0,7 mL min^-1 e com a concentração de 1 g L^-1 de KBr são apresentadas na figura 1.

O ajuste do modelo CDE aos pontos da curva de eluição do KBr para os três solos nas vazões de 0,2; 0,45; e 0,7 ml min^-1 foi excelente (Figura 1).

Em todos os ensaios, além de ter ocorrido simetria no trecho ascendente e descendente de cada uma das curvas, também as curvas de eluição e as ajustadas pelo modelo CDE passaram pelo ponto (0,5 C/C₀; 1,0 V/V₀) (Figura 1), confirmando a boa qualidade do KBr como traçador químico, estando esses resultados de acordo com a literatura (Martins & Mermoud, 1998; Costa et al., 2006; Milfont et al., 2008; Carmo et al., 2010).

Propriedades hidrodinâmicas dos solos

Os valores médios das condições experimentais e dos parâmetros hidrodispersivos dos ensaios de deslocamento miscível com o KBr, na concentração de 1 g L^-1, para os solos Franco-arenoso 1, Franco-arenoso 2 e Franco-argilossiltoso, na vazão de 0,2; 0,45; e 0,7 mL min^-1, são apresentados no quadro 4.

Os valores médios do fator de retardamento, R, nos Franco-arenoso 1 e Franco-argilossiltoso, nas três vazões utilizadas e na concentração de 1 g L^-1, ficaram próximos à unidade, indicando que o KBr não sofreu interações nesses solos (adsorção ou exclusão) (Quadro 4). No Franco-arenoso 2, foi constatada ligeira adsorção, possivelmente pela presença de goethita na fração argila.

A dispersividade, λ, foi obtida considerando-se a relação linear entre o coeficiente de dispersão hidrodinâmico, D, e a velocidade média da água nos poros, v, ou seja, D= λ v. A partir dos valores de D ajustados e de v, determinou-se λ para os três solos (Quadro 4), uma vez que a dispersividade está diretamente relacionada com o diâmetro médio das partículas de solo (Milfont et al., 2006).

O solo Franco-arenoso 2 foi o mais dispersivo dos três, possivelmente em razão de diferenças no tamanho e arranjo de suas partículas, que tendem a induzir a formação de poros largos, favoráveis à importante distribuição de velocidades da água.

Em relação ao número de Péclet, verificou-se que nos três solos e nas três vazões estudadas, os valores de Pe foram maiores que 10 (Quadro 4), indicando que o processo predominante de transporte é do tipo convectivo (Novy Quadri, 1993).

Transporte reativo e caracterização hidrodinâmica do naftaleno

Os ajustes realizados pelo modelo CDE, utilizando os valores de D obtidos nos ensaios com KBr, aos pontos das curvas médias de eluição do naftaleno, nas colunas do solo Franco-arenoso 1, Franco-arenoso 2 e Franco-argilossiltoso, na vazão de 0,45 mL min^-1 e na concentração de 20 mg L^-1, são apresentados na figura 2.

O solo Franco-arenoso 1 apresentou menor tempo de retenção para o naftaleno (Figura 2), enquanto Franco-arenoso 2, o maior. O modelo CDE ajustou-se bem aos pontos das curvas médias de eluição do naftaleno dos três solos (Figura 2).

O aumento no tempo de retenção do naftaleno no solo Franco-arenoso 2 pode ser explicado pela presença de goethita. Conforme demonstrado na literatura (Tunega et al., 2009; Muller et al., 2007), tem sido observada uma forte tendência de sorção de compostos hidrofóbicos, como o naftaleno, com os óxidos de ferro, os quais contribuem de forma significativa junto com a matéria orgânica para a retenção do naftaleno no referido solo.

Os valores médios das condições experimentais e dos parâmetros hidrodispersivos dos ensaios de deslocamento miscível com naftaleno na concentração de 20 mg L^-1 para os solos Franco-arenoso 1, Franco-arenoso 2 e Franco-argilossiltoso, na vazão de 0,45 mL min^-1, são apresentados no quadro 5.

O solo Franco-arenoso 2 apresentou o maior valor do fator de retardamento (Quadro 5), sendo o mais reativo, enquanto o Franco-arenoso 1, o menor, estando esses resultados de acordo com literatura (Lee et al., 2002).

O solo Franco-arenoso 2, dentre os demais, foi o que evidenciou o menor risco de contaminação do lençol freático existente nessa região, frente à possível contaminação pelo naftaleno, enquanto o Franco-arenoso 1, o maior risco.

CONCLUSÕES

1. Os solos arenosos oferecem menor resistência à lixiviação do naftaleno.

2. O solo Franco-arenoso 1 evidenciou o maior risco de contaminação por lixiviação dentre os solos estudados.

3. O modelo CDE apresentou bom ajuste aos pontos das curvas médias de eluição do traçador (KBr) e do naftaleno, indicando que os processos de convecção-dispersão são os predominantes envolvidos na transferência reativa desse composto.

4. A presença de goethita pareceu ser a principal responsável pelo retardamento da lixiviação do naftaleno, no conjunto de amostras estudadas.

5. A presença da matéria orgânica exerceu grande efeito reativo no processo de transferência do naftaleno nos três solos.

AGRADECIMENTOS

À CAPES (Processo N^o 677/10 - CAPES/COFECUB), à FACEPE (Processo N^o 1262-3.07/08), ao CNPq (Processos N^o 475574/2009-5, Edital MCT/CNPq 14/2009 - Universal N^o 01/2001, 305727/2009-5, 371970/2011-3 e 374299/2011-0), pela concessão das bolsas e dos recursos necessários ao desenvolvimento deste trabalho.

LITERATURA CITADA

Recebido para publicação em 2 de agosto de 2012

Aprovado em 24 de junho de 2013

• ANDRADE, M.C.O. Pernambuco e o trópico. R. Inst. Est. Bras., 45:11-20, 2007.
• CARMO, A.I.; ANTONINO, A.C.D.; NETTO, A.M. & CORREA, M.M. Caracterização hidrodispersiva de dois solos da região irrigada do Vale do São Francisco. R. Bras. Eng. Agric. Amb., 14:698-704, 2010.
• COSTA, C.T.F.; ANTONINO, A.C.D. & NETTO, A.M. Ensaios de deslocamento de líquido miscível na determinação dos parâmetros hidrodispersivos de um solo aluvial. R. Bras. Rec. Hídrol., 11:111-122, 2006.
• EMPRESA BRASILEIRA DE PESQUISA AGROPECUÁRIA - EMBRAPA. Serviço Nacional de Levantamento e Conservação do Solo. Manual de métodos de análise de solo. 2.ed. Rio de Janeiro, 1997. 212p.
• JACKSON, M.L.; LIM, C.H. & ZELAZNY, L.W. Oxides, hydroxides, and aluminosilicates. In: KLUTE, A. ed. Methods of soil analysis: Physical and mineralogical methods. Madison, American Society of Agronomy, 1986. Part 1. 49p.
• JOHNSON, W.P.; MARTIN, M.J.; GROSS, M.J. & LOGAN, B.E. Facilitation of bacterial transport through porous media by changes in solution and surface properties. Colloids Surf., 107:263-271, 1996.
• LEE, J.; HUNDAL, L.S.; HORTON, R. & THOMPSON, M.L. Sorption and transport behaviour of naphthalene in an aggregated soil. J. Environ. Qual., 31:1716-1721, 2002.
• HOLLOSI, L. & WENZL, T. Development and optimisation of a dopant assisted liquid chromatographic-atmospheric pressure photo ionisation-tandem mass spectrometric method for the determination of 15 + 1 EU priority PAHs in edible oils. J. Chromatogr. A., 1218:23-31, 2011.
• LOVELAND, P.J. & WHALLEY, W.R. Particle size analysis. In: SMITH, K.A. & MULLIS, CH.E., eds. Soil analysis: Physical methods. New York, Marcel Dekker, 1991. p.271-328.
• MARTINS, J.M.F. & MERMOUD, A. Sorption and biodegradation of four nitroaromatic herbicides in mono and multi-solute saturated/unsaturated soil batch systems. J. Contam. Hydrol., 33:187-210, 1998.
• MCCREADY, S.; SLEE, D.J.; BIRCH, G.F. & TAYLOR, S.E. The distribution of polycyclic 12 aromatic hydrocarbons in surficial sediments of Sydney Harbour, Australia. Mar. Pollut. Bull., 40:999-1006, 2000.
• MIELKE, H.W.; WANGB, G.; GONZALES, C.R.; LEE, B.; QUACH, V.N. & MIELKEC, P.W. PAH and metal mixtures in New Orleans soils and sediments. Sci. Total Environ., 281:217-227, 2001.
• MOTELAY_MASSEI, A.; OLLIVON, D.; GARBAN, B.; TEIL, M.J.; BLANCHARD, M. & CHEVREUIL, M. Distribution and spatial trends of PAHs and PCBs in soils in the Seine River basin, France. Chemosphere, 55:555-565, 2004.
• MULLER, S.; TOTSCHE, K.U. & KÖGEL-KNABER, I. Eur. J. Soil Sci., 58:918-931, 2007.
• MILFONT, M.L.B.; ANTONINO, A.C.D.; MARTINS, J.M.F.; NETTO, A.M.; GOUVEIA, E.R. & CORREA, M.M. Transporte do paclobutrazol em colunas de solos. R. Bras. Ci. Solo, 32:2165-2175, 2008.
• MILFONT, M.L.B.; ANTONINO, A.C.D.; MARTINS, J.M.F.; NETTO, A.M. & CORREA, M.M. Caracterização hidrodispersiva de dois solos do Vale do Rio São Francisco. R. Bras. Ci. Agron., 1:81-87, 2006.
• NOVY QUADRI, M.G. Transferts de solutés dans les sols saturés et non saturés. Application au pentachlorophénol. Grenoble, Université Joseph Fourier, 1993. 204p. (Tese de Doutorado)
• PARKER, J.C. & van GENUCHTEN, M.T.H. Boundary conditions for displacement experiments through short laboratory soil columns. Soil Sci. Soc. Am. J., 48:703-708, 1984.
• QI, S.; YAN, J.; ZHANG, G.; FU, J.; SHENG, G.; WANG, Z.; TONG, S.M.; TANG, U.W. & MIN, Y. Distribution of polycyclic aromatic hydrocarbons aerosols and dustfall in Macao. Environ. Monit. Assess., 72:115-127, 2001.
• SMITH, J.A.; SHANNON, L.B. & SUSAN, E.B. Sorption and permeability of gasoline hydrocarbons in organobentonite porous media. J. Hazard. Mater., 96:91-97, 2003.
• SPOSITO, G. The chemistry of soils. New York, Oxford University Press, 1989. 277p.
• RUIZ, Y.; SUAREZ, P.; ALONSO, A.; LONGO, E.; VILLAVERDE, A. & JUAN, A.S. Environmental quality of mussel farms in the Vigo estuary: Pollution by PAHs, origin and effects on reproduction. Environ. Pollut., 159:250-265, 2011.
• SCHWARZ, K.; GOCHT, T. & GRATHWOHL, P. Transport of polycyclic aromatic hydrocarbons in highly vulnerable karst systems. Environ. Pollut., 159:133-139, 2011.
• TUNEGA, D.; GERZABEK, M.; HABERHAUER, G.; TOTSCHE, K. & LISCHKA, H.J. Model study on sorption of polycyclic aromatic hydrocarbons to goethite. J. Colloid Interf. Sci., 330:244-249, 2009.
• YAN, W.; CHI, J.; WANG, Z.; HUANG, W. & ZHANG, G. Spatial and temporal distribution of polycyclic aromatic hydrocarbons (PAHs) in sediments from Daya Bay, South China. Environ. Pollut., 157:1823-1830, 2009.
• ZHOU, J.L. & MASKAOUI, K. Distribution of polycyclic aromatic hydrocarbons in water and surface sediments from Daya Bay, China. Environ. Pollut., 121:269-281, 2003.
• WOLICKA, D.; SUSZEK, A.; BORKOWSKI, A. & BIELECKA, A. Application of aerobic microorganisms in bioremediation in situ of soil contaminated by petroleum products. Biores. Technol., 100:3221-3227, 2009.

1

Parte da Tese de Doutorado do primeiro autor apresentada ao Programa de Pós-Graduação em Tecnologias Energéticas e Nucleares da Universidade Federal de Pernambuco - UFPE, Bolsista da FACEPE.

Datas de Publicação

• Publicação nesta coleção
  17 Dez 2013
• Data do Fascículo
  Out 2013

Histórico

• Recebido
  02 Ago 2012
• Aceito
  24 Jun 2013