Simulação da Química da Atmosfera Poluída por Automóveis Movidos a Álcool

Simulation of the urban atmosphere chemistry polluted by alcohol fueled automobiles

troposphere; aldehyde pollution; numerical modeling

Simulação da Química da Atmosfera Poluída por Automóveis Movidos a Álcool

A. S. Pimentel e G. Arbilla

Departamento de Físico-Química - Instituto de Química - Universidade Federal do Rio de Janeiro - CT - Bloco A - Cidade Universitária - 21949-900 - Rio de Janeiro - RJ

Recebido em 7/5/96; aceito em 16/10/96

Simulation of the urban atmosphere chemistry polluted by alcohol fueled automobiles. Numerical modeling studies of urban atmospheres were performed with a photochemical box submodel which incorporates primary emission rates, temperature and solar irradiance rate coefficients. Actinic flux calculations were made for an appropriate latitude for Rio de Janeiro and initial concentrations and emission rates were taken from experimental data. The relative importance of individual reactions is discussed as well as the effect of enhanced aldehyde emissions from vehicles using ethanol and ethanol-containing fuel.

Keywords: troposphere; aldehyde pollution; numerical modeling.

INTRODUÇÃO

O estudo da atmosfera e dos problemas que acontecem nela tem despertado crescente interesse nos últimos anos em diversas áreas da ciência. A exploração da atmosfera e o desenvolvimento de modelos permitiu a melhor compreensão dos fenômenos envolvidos e levou à conclusão de que as atividades humanas podem produzir mudanças significativas na composição da mesma^1-3. Contudo, são muitos os pontos que ainda não estão bem esclarecidos, especialmente aqueles que são importantes só em alguns países, como as conseqüências do uso de combustíveis alternativos (por exemplo, gás natural e álcool).

A poluição do ar provocada por veículos automotores leves e pesados é um problema muito sério nos grandes centros urbanos brasileiros e estrangeiros^4-8. Sendo assim, um país como o Brasil, com uma frota tão ampla e diferenciada, emite uma variedade enorme de poluentes primários⁹ que, por sua vez, afetam grandemente a concentração de poluentes secundários.

Com o início da utilização do álcool como combustível, o Brasil tem sofrido até hoje problemas complexos de poluição ambiental, pois possui uma frota veicular bastante diversificada, onde existem carros antigos e novos, movidos por diferentes combustíveis, e veículos pesados.

Um dos principais objetivos dos estudos na área das ciên-cias atmosféricas é relacionar as emissões de poluentes primários em uma determinada região com as concentrações de poluentes secundários e com as mudanças das diversas espécies constituintes da troposfera. Para isto é necessário construir um modelo completo e comparar com os resultados experimentais. Este tipo de estudo, tipicamente, é realizado em duas etapas²: primeiramente é construído um "submodelo" puramente químico que, quando possível, é avaliado com medições em uma câmara de gases onde as reações são realizadas em condições controladas, independentemente das perturbações meteorológicas; e em uma segunda etapa, é construído um "submodelo" para reproduzir as condições meteorológicas do local. De uma forma geral, um modelo completo para a troposfera poluída deve incluir três aspectos:

· As emissões dos poluentes primários para cada horário do dia e para uma determinada região.

· As condições meteorológicas do local para cada estação do ano, como a temperatura, umidade relativa, ventos, difusão e movimentos de massas de ar que, junto ao tempo de vida de cada uma das espécies, irão determinar se o efeito de um determinado poluente será local, regional ou global.

· Os processos químicos.

Os processos que acontecem na atmosfera são complicados^10,11 e envolvem um grande número de espécies e interações simultâneas, tornando impossível o tratamento dos mecanismos envolvidos de forma analítica. É preciso, então, resolver de forma numérica o sistema de equações diferenciais acopladas resultante do modelo proposto^12-17.

Na resolução de um modelo completo, tridimensional (3D), para a troposfera, os processos de transporte vertical e horizontal envolvem entre 5 e 30% do tempo computacional^13-17. O maior custo computacional está relacionado aos fenômenos químicos, pois o conjunto de reações elementares origina um sistema de equações diferenciais acopladas, geralmente muito rígido, que deverá ser resolvido no tempo e no espaço. Entretanto, os processos físicos de transporte não podem ser negligenciados, especialmente na descrição da troposfera, onde as emissões (naturais e antropogênicas) e a difusão contribuem significativamente para a velocidade dos processos globais. A solução de um modelo 3D pode ser feita de forma paralela¹³ ou de forma seqüencial¹⁷, resolvendo em primeiro lugar os submodelos meteorológico e de emissões e, posteriormente, os submodelos químicos (reações homogêneas, heterogêneas, formação e deposição de aerossóis). Este tipo de tratamento requer tempo computacional e espaço em memória consideráveis, especialmente no período de otimização das condições de contorno e de cada uma das partes que compõem o modelo. Geralmente a parte química do modelo precisa de uma resolução temporal e espacial muito menor. Porém, o submodelo químico pode ser resolvido e analisado isoladamente, no tempo, para determinadas condições espaciais, usando um método apropriado para equações diferenciais ordinárias em um computador pessoal. O procedimento equivale à solução de um ponto representativo da matriz de pontos que descrevem o espaço tridimensional. Sendo assim, torna-se importante realizar primeiro um estudo detalhado da química do processo e da importância relativa das reações a fim de obter um submodelo químico apropriado e o mais reduzido possível, o que implica menos tempo computacional. Nesse estudo os processos físicos podem ser descritos em forma simplificada através do chamado "modelo em caixa" que envolve, também, só equações diferenciais ordinárias. O submodelo otimizado é o ponto de partida para a implementação do modelo completo. Existem atualmente várias propostas para simular a química da atmosfera urbana, realizadas com o objetivo de prever os níveis de ozônio e compostos nitrogenados em diferentes condições iniciais, a fim de gerar dados e conclusões para estabelecer estratégias de controle ambiental^17-20. Esses modelos são conhecidos genericamente como modelos cinéticos empíricos (EKMA), e não pretendem explicar alguns dos processos químicos relevantes da atmosfera urbana de cidades brasileiras, como as reações de aldeídos e compostos relacionados.

Um outro aspecto importante a ser considerado na modelagem da atmosfera é a mudança do fluxo solar actínico, que varia com o horário do dia e mês do ano²¹. Além disso, este varia com a latitude, comprimento de onda e altitude, fazendo com que cada cidade, ou seja, cada problema tenha características bem próprias. A figura 1 mostra como o ângulo zenital solar varia com o horário do dia e com mês do ano e a figura 2 mostra como varia o fluxo solar com o ângulo zenital solar e com o comprimento de onda para uma latitude de 20^oN. A latitude de 20^oN foi escolhida pois há uma extrema dificuldade de se obter valores deste tipo para o hemisfério sul² e, ainda, porque a cidade do Rio de Janeiro possui latitude próxima de 20^oS (23^oS), o que equivale ao inverso da figura 1, ou seja, ao invés de ter ângulo zenital solar próximo a 0^o no mês de junho, têm-se esse ângulo no mês de janeiro. Os dados representados não consideram a presença de nuvens e aerossóis, que podem diminuir sensivelmente a quantidade de energia radiante que atinge o solo. A sua influência nos processos fotoquímicos não está bem determinada e se constitui numa fonte importante de erro.

Os comprimentos de onda representados na figura 2 estão compreendidos entre 290 e 360 nm para uma melhor visualização da mesma, mas no cálculo foram considerados comprimentos de onda entre 290 e 600 nm pois existem compostos que absorvem nesta faixa de comprimento de onda. As figuras 1 e 2 apresentam resultados teóricos para o hemisfério norte²², como já foi dito anteriormente. Entretanto, alguns cálculos preliminares realizados para o hemisfério sul, utilizando o programa JSPECTRA²³, demonstraram que a aproximação feita não afeta o resultado deste trabalho.

O objetivo deste trabalho é modelar as reações envolvendo aldeídos em diferentes condições, determinando os intermediários, produtos finais e a importância relativa das diferentes rea-ções envolvidas. Visa, também, correlacionar os aldeídos com a presença de outras espécies na atmosfera tais como HOx, NO_x, O₃, nitrato de peroxiacetila (PAN), peróxidos nitrados e ácidos em geral.

METODOLOGIA

Método numérico

O método numérico e os programas de cálculo utilizados foram discutidos detalhadamente em trabalhos anteriores^{24, 25}. Por este motivo, faremos aqui só uma breve descrição da metodologia utilizada.

O modelo pode ser descrito matematicamente^{26, 27} através de um conjunto de equações diferenciais ordinárias, acopladas, para as quais são conhecidas as condições iniciais. O sistema de equações é do tipo :

onde Y_i representa a espécie química i, e k representa o vetor das constantes de velocidade. Assim, para um mecanismo que envolve N espécies, resultará um sistema de N equações diferenciais ordinárias (ODE) acopladas, que podem ser resolvidas usando um método numérico apropriado. Para um sistema rígido, os dois métodos mais utilizados são o método de Gear²⁸ e o método Runge-Kutta-4-semi-implícito²⁹.

Neste trabalho foram utilizados ambos os métodos. A integração das equações de velocidade foi realizada usando o programa KINETICS, desenvolvido neste laboratório. Este programa, escrito em FORTRAN 77, utiliza o método de Gear, implementado nas subrotinas LSODE e VODE³⁰ cedidas pelo Prof. Hindmarsh (Livermoore). O programa computacional foi extensivamente discutido na referência 24. Neste trabalho, o programa foi adaptado e modificado para a inclusão da dependência das constantes de velocidade com a temperatura e fluxo solar, e para a inclusão da emissão de poluentes em função da hora do dia.

A importância relativa de cada etapa elementar e a evolução temporal de cada espécie foram calculadas usando o pacote de domínio público KINAL³¹, baseado no método Runge-Kutta-4-semi-implícito e no método de derivadas diretas. Nesse caso, os fundamentos e metodologia do cálculo já foram discutidos com detalhe²⁵. Note que este tipo de programa não permite a consideração de fenômenos como emissões veiculares, mudança da temperatura ou da radiação solar.

DESCRIÇÃO DO MODELO

Com o intuito de realizar uma modelagem exploratória, inicialmente foram resolvidos modelos simples, um deles envolvendo só a fotólise primária e secundária do formaldeído, um modelo mais elaborado incluindo as reações com os oxidantes HO_X, um modelo mais abrangente das reações do formaldeído com os oxidantes HO_X e um modelo da fotodecomposição do acetaldeído envolvendo a fotólise do acetaldeído com suas principais reações. Os resultados permitiram fazer um estudo preliminar para compreender a química dos aldeídos, demonstrando que os sistemas de equações resultantes são muito rígidos e acoplados, e que, para explicar os fenômenos químicos na troposfera, é necessário considerar a interação entre as diversas espécies, especialmente as famílias de espécies HO_x, NO_x e O_x.

Foi construído um modelo considerando um conjunto abrangente de espécies e reações, e a contribuição de cada um deles ao processo total foi estimada. Estes resultados permitiram eliminar mais de vinte reações elementares e algumas espécies que não contribuiram significativamente, reduzindo o tempo computacional. O modelo final é mostrado no ^Apêndice, e inclui 35 espécies e 55 reações químicas, que constituem um modelo bastante completo capaz de reproduzir satisfatoriamente a química da troposfera urbana. Como indicado no ^Apêndice, esse modelo ainda pode ser reduzido sem perda apreciável na sua capacidade de predição.

O modelo escrito em FORTRAN 77 foi implementado com um algoritmo que simula a emissão de poluentes na troposfera e outro que simula a variação do fluxo solar e temperatura com o decorrer do dia.

As expressões para os coeficientes de velocidade específica de cada reação foram extraídas da bibliografia^2,32 e corrigidas, quando necessário, para a temperatura ambiente. No caso das reações fotoquímicas, o fluxo solar actínico foi calculado a partir do ângulo zenital solar de cada mês do ano e horário do dia a uma latitude de 20^oS (próxima ao Rio de Janeiro). As concentrações iniciais e os dados de emissão de cada espécie estável foram obtidos de dados experimentais coletados da literatura^2,6-9.

A determinação da velocidade de emissão de poluentes, requer um inventário de emissões para a região que está sendo modelada¹². Este tipo de informação não está disponível para as cidades brasileiras, por isso foram utilizados dados estimados por Tanner et. al.⁴, que são consistentes com outros dados disponíveis⁹.

RESULTADOS E DISCUSSÃO

Inicialmente os perfis de concentrações de todas as espécies em função do tempo foram calculados utilizando as condições iniciais descritas na tabela 1, valores médios de temperatura (20^oC e 35^oC para o inverno e verão, respectivamente), e valores de fluxo solar correspondentes às 12 horas. Estes cálculos não consideram emissões veiculares, e as condições iniciais são valores médios para a atmosfera do Rio de Janeiro, segundo os dados experimentais disponíveis^6-8. As equações diferenciais foram resolvidas utilizando os programas KINETICS (LSODE e VODE) e KINAL, e os resultados são mostrados nas figuras 3-5. A concordância entre as curvas obtidas com ambos os programas é bastante boa, o que demonstra que os resultados destes métodos numéricos convergem para valores bem próximos.

Posteriormente, as concentrações de todas as espécies foram calculados em função do horário do dia, das 6 às 15 horas, usando o programa KINETICS (Figuras 6-9). A variação diurna da temperatura utilizada para o verão foi de 22^oC a 40^oC e para o inverno foi de 19^oC a 27^oC. Já a variação do ângulo zenital solar para o verão foi de 2^o a 83^o e para o inverno foi de 43^o a 85^o. O modelo foi resolvido com as condições iniciais e emissões veiculares mostradas na tabela 1, que correspondem a níveis de poluição urbana da Cidade do Rio de Janeiro (Vila Izabel) obtidos da literatura^4,5,7.

Foi calculada também a contribuição de cada uma das reações na velocidade total da mudança de concentração de cada uma das espécies. Estes cálculos foram feitos para uma temperatura de 20^oC no inverno e 35^oC no verão, com as condições iniciais da tabela 1 e coeficientes fotoquímicos que correspondem aos valores de 12 horas de um dia típico de verão e inverno, usando o programa KINAL.

A principal via de consumo do formaldeído na atmosfera natural é a fotodecomposição. O tempo de vida desta espécie com relação a fotólise é de aproximadamente 6,3 horas no verão e 8,1 horas no inverno num dia claro. Em ambientes lim-pos, esta constitui a única via importante de consumo do formaldeído. Entretanto, nas condições iniciais de modelagem (atmosfera levemente poluída ou com poluição moderada), 44% da decomposição total é causada pelas reações de HO. Note que em ambientes com concentrações de acetaldeído relativamente altas, o formaldeído é um poluente secundário, formado eficientemente a partir de CH₃O. Já no caso do acetaldeído, a sua principal via de reação é com o HO para formar CH₃CO₃, precursor do PAN. O tempo de vida do acetaldeído com relação a fotólise é de aproximadamente 3,3 dias no verão e 5 dias no inverno, sendo este processo muito menos importante em dias nublados ou com pouca luminosidade. A figura 7 ilustra a relação [acetaldeído]/[formaldeído] fornecida pelo modelo e os valores superiores a 1 são típicos de ambientes urbanos com altas emissões de aldeídos.

A análise de importância das reações mostra claramente que o aumento considerável na concentração de PAN é devido às emissões veiculares de aldeídos. As concentrações calculadas correspondem a níveis de poluição moderada segundo os critérios de poluição aceitos por grupos de pesquisa². O principal processo de formação do PAN é a reação 47: CH₃CO₃ +NO₂ = CH₃CO₃NO₂ , que compete com a reação 49: CH₃CO₃ +NO que em atmosfera de oxigênio produz radicais CH₃O₂. O CH₃CO₃ é formado a partir da reação de acetaldeído com HO. Pela modelagem pode-se observar que a formação de PAN au-menta bastante com o aumento da umidade relativa do ar e a diminuição da temperatura (Figura 10) e levemente com o aumento da radiação solar para uma mesma temperatura. O aumento da umidade relativa e da radiação solar ocasionam um acréscimo na concentração de radicais HO. Alguns cálculos preliminares para valores de umidade relativa do ar entre 30% e 100%, demostraram que a concentração de PAN aumenta em torno de 40% no final do dia com o aumento da umidade nesse intervalo. Já o aumento da radiação solar, de valores típicos do inverno a valores típicos do verão, contribui muito pouco para um aumento na concentração de PAN. Maiores níveis de emissões de NO irão produzir uma diminuição na concentração de PAN (Figuras 11 e 12). Por outro lado, o PAN é termicamente instável e decompõe-se unimolecularmente num processo muito dependente da temperatura. O tempo de vida a 10^oC, t =1/k₄₈, é de 8,6 horas, entretanto a 30^oC é de somente 22 minutos. Assim o PAN atua como um composto de reserva de NO₂, ou seja, em temperaturas e/ou concentrações de NO baixas, a decomposição do PAN é inibida. Pelo contrário, com o aumento da temperatura ou altos níveis de NO no ar, os processos de decomposição e a reação competitiva para formar CH₃CO₂ crescem em importância e a quantidade de PAN resultante é menor. A natureza não linear do sistema torna difícil predizer o seu comportamento para um determinado conjunto de condições iniciais. Dependendo da temperatura, umidade, radiação solar e concentração dos NO_x, as concentrações de PAN mudarão, sendo maiores em geral para dias úmidos, frescos e com sol.

As concentrações de ácidos orgânicos, como ácido fórmico e acético, calculadas com o modelo completo, foram relativamente baixas. Porém este resultado não pode ser usado como uma indicação de um baixo índice de poluição porque o modelo não deu ênfase à descrição das reações de hidrocarbonetos e compostos orgânicos como cetona, podendo ser estas as principais vias de formação de ácidos orgânicos. Entretanto, a concentração de ácidos inorgânicos (HNO₃) apresentou-se muito grande (10-30 ppb), o que corresponde a altos níveis de poluição (Figura 8). Este resultado deve-se à lentidão do processo de decomposição fotoquímica do ácido nítrico. Por exemplo, com ângulo zenital solar de 0^o, o tempo de vida desta espécie, calculado considerando só a decomposição fotoquímica, é de aproximadamente 20 dias. Contudo os processos de deposição na superfície terrestre e dissolução em água, principal motivo da chuva ácida, são relativamente importantes e contribuem para diminuir o tempo de vida do ácido na troposfera. Estes processos não foram considerados neste modelo simplificado, sendo assim os valores calculados provavelmente mais altos que os valores reais.

A figura 9 mostra que a concentração de ozônio no verão, em todos os horários do dia, é sempre maior que no inverno. Isso se deve ao balanço das reações de formação e consumo da espécie. A formação de O₃ acontece principalmente pela reação 1: O(³P) + O₂ →O₃ e o consumo pela reação 6: NO + O₃ →NO₂ + O₂. A análise de importância das reações mostra que a relação entre as contribuições das reações 1 (formação) e 6 (decomposição) é ligeiramente maior no verão do que no inverno (Tabela 2 e Figura 9). Inclusive, existem evidências expe-rimentais², obtidas em alguns lugares da Europa e Estados Unidos, de que a concentração de ozônio é maior no verão do que no inverno.

CONCLUSÕES

Neste trabalho, apresentamos um submodelo puramente químico com o objetivo de correlacionar a emissão veicular de poluentes na atmosfera de centros urbanos com a formação de poluentes secundários, e as mudanças na composição do ar. O modelo simula também a influência da variação do fluxo solar e da temperatura nos processos químicos. Ainda que muito sim-ples, este modelo pode explicar muitos dos aspectos essenciais dos fenômenos que ocorrem na atmosfera, e pode ser considerado como uma indicação de que os processos fotoquímicos são muito importantes na descrição da atmosfera de centros urbanos como o Rio de Janeiro, diferentemente do que acontece em outras cidades, onde os fatores topográficos e meteorológicos são os que determinam a composição atmosférica³³.

A modelagem numérica da fotodecomposição dos aldeídos mostrou uma boa concordância com os dados experimentais^4-7 coletados em diferentes horas do dia. A concordância é muito boa para quase todas as espécies, com exceção do HNO₃. As concentrações de PAN e O₃ apresentam um máximo para o qual a localização e intensidade dependem muito das condições atmosféricas (temperatura, umidade e radiação solar) e emissões veiculares. Para as condições predominantes do Rio de Janeiro no inverno, a concentração máxima é em torno de 0,30 ppb para o PAN e 7 ppb para o O₃ nas primeiras horas da tarde. Estes resultados possuem uma boa concordância com os valores experimentais para a maioria dos dias monitorados^4,34. Contudo, o modelo mostra que alguns episódios de maior poluição poderão ser observados em dias úmidos e mais frios que o normal, ou com condições meteorológicas especiais, como acontecido recentemente na cidade do Rio de Janeiro³⁵. Note-se que o fluxo solar que atinge o solo depende de outros fatores além daqueles mencionados acima, como a presença de nuvens e aerossóis, que diminui sensivelmente a quantidade de energia radiante. Neste modelo simplificado, este efeito não foi considerado, o que leva a uma maior incerteza nas determinações quantitativas.

A relação CH₃CHO/HCHO é uma característica fundamental dos países que utilizam álcool como combustível. Em lugares muito poluídos, como túneis com grande fluxo de veículos, esta relação pode atingir valores muito altos, em torno de 4-5. Em lugares moderadamente poluídos observam-se valores em torno de 1,5. No caso de países que não utilizam álcool como combustível esta relação é sempre menor que 1, e o formaldeí-do é, certamente, o aldeído mais abundante, com concentrações que superam a concentração dos outros aldeídos em uma ordem de grandeza. Em atmosferas limpas, a principal fonte de formaldeído é a seqüência de reações 29, 33 e 34. Entretanto, nas condições desta modelagem, as principais fontes de formaldeído são as emissões veiculares e a seqüência de reações 44, 45, 49, 33 e 34. Em ambos os casos acontece a reação: CH₃O + O₂ →HCHO + HO₂. Em ambientes não poluídos o CH₃O é um produto secundário do CH₄ naturalmente presente na atmosfera. Entretanto, pelo estudo da tabela 2, pode-se deduzir que, nesta modelagem, o CH₃O provém majoritariamente do CH₃CHO.

A análise de importância revelou que os processos de remoção, como a fotólise da água, não são importantes. Nesta análise e no cálculo dos coeficientes de sensibilidade das reações observou-se que as reações que se encontram marcadas com asterisco (^Apêndice) são as menos importantes. O modelo reduzido, sem a consideração destes processos, fornece resultados igualmente satisfatórios para as principais espécies.

O etanol, principal componente do álcool combustível, provoca um aumento acentuado da emissão de aldeídos na troposfera. Por outro lado, dos resultados deste trabalho, pode-se concluir que, embora a emissão direta de poluentes por veículos movidos a álcool seja menor do que a feita por veículos movidos a gasolina, os aldeídos sofrem reações que modificam apreciavelmente as concentrações de HO_x e NO_x e aumentam as concentrações de oxidantes secundários e aerossóis. Certamente, a elaboração de um modelo mais completo torna-se necessária para obter resultados quantitativos de melhor qualidade, ou seja, provavelmente a concentração de ácidos orgânicos deve aumentar quando forem incluídas no modelo reações de hidrocarbonetos, radicais CH₃CO₂ e compostos oxigenados. Atualmente está sendo desenvolvido um novo modelo no qual a química do etanol e dos hidrocarbonetos voláteis é considerada de forma explícita. Também serão considerados os processos de deposição e difusão em forma mais detalhada.

Finalmente pode-se notar que do ponto de vista numérico e computacional, o método numérico que mostrou melhor desempenho para a solução deste sistema foi o Método de Gear, nas sub-rotinas LSODE e VODE, pois demandam menos tempo computacional e menor número de passos de integração para atingir convergência. Contudo os resultados são perfeitamente concordantes para qualquer uma das três sub-rotinas utilizadas, o que indica que os resultados não dependem apreciavelmente dos métodos numéricos que foram utilizados.

AGRADECIMENTOS

Os autores agradecem o apoio financeiro do CNPq, CAPES e FAPERJ, e a colaboração dos Professores Antônio H. Miguel e Jailson B. de Andrade. Agradecem também aos Prof. A. C. Hindmarsh (Livermore) e T. Turányi (Leeds) por terem cedido o LSODE/VODE e KINAL, respectivamente, para uso acadêmico e de pesquisa.

REFERÊNCIAS E NOTAS

1. Wayne, R. P.; Chemistry of Atmospheres, Clarendon Press, Oxford 1991.

2. Finlayson-Pitts, B. J.; Pitts, J. N. Jr.; Atmospheric Chemistry: Fundamentals and Experimental Techniques, John Wiley & Sons, Chichester 1986.

3. Seinfeld, J. H.; Atmospheric Chemistry and Physics of Air Pollution, John Willey, Chichester 1986.

4. Tanner, R. L.; Miguel, A. H.; de Andrade, J. B.; Gaffney, J. S.; Streit, G. E.; Environ. Technol. 1988, 22, 1026.

5. Miguel, A. H.; Quím. Nova 1992, 15, 118.

6. Daisey, J. M.; Miguel, A. H.; de Andrade, J. B.; Pereira, P. A. P.; Tanner, R. L.; J. Air Pollut. Control Assoc. 1987, 37, 15.

7. Grosjean, D.; Miguel, A. H.; Tavares, T. M.; Atm. Environment 1990, 24B, 101.

8. de Andrade, J. B; de Andrade, M. V.; Pinheiro, H. L. C.; Regional Conference on Global Change, São Paulo 1995, 151.

9. Relatório do Laboratório de Emissões Veículares (LEV) do CENPES - PETROBRAS, 1995.

10. Leone, J. A.; Seinfeld, J. H.; Atmospheric Environment 1988, 19, 437.

11. de Aurelliano, J. V.; Acta Chimica Hungarica - Models in Chemistry, 1992, 129, 889.

12. Rethinking the Ozone Problem in Urban and Regional Air Pollution, editado pelo Comittee on Tropospheric Ozo-ne Formation and Measurement, National Academy Press, Washington, D. C. 1991

13. Dabdub, D.; Seinfeld, J. H.; Atmospheric Environment 1994, 28, 1679.

14. Jang, J. C.; Jeffries, H. E.; Byun, D.; Pleim, J. E.; Atmospheric Environment 1995, 29, 3085.

15. Jang, J. C.; Jeffries, H. E.; Tonnesen, S.; Atmospheric Environment 1995, 29, 3101.

16. Informações recentes sobre o estado da arte, novos modelos e bases de dados, podem ser obtidos na INTERNET. Alguns dos endereços mais importantes são:

a) http://airsite.unc.edu/, página da Universidade da Carolina do Norte - EUA e Universidade de Leeds - Reino Unido dedicada a informação sobre química atmosférica e poluição. Possui conexão com páginas da NASA, da EPA (Environment Protection Agency, EUA) e de outras universidades.

b) http://www.agu.org/, página da American Geophysical Union, sociedade científica com a participação de 115 países nas áreas de química atmosférica, oceanografia e geofísica, e editora das principais revistas da área.

c) http://www.igbp.kva.se/, página do IGBP (International Geosphere-Biosphere Programme).

17. Chang, J. S.; Brost, R. A.; Isaken, S. A.; Madronich, S.; Middleton, P.; Stockwell, W. R.; Walcek, C. J.; J. Geophys. Res. 1987, 92, 14681.

18. Jin, S.; Demerjian, K.; Atmospheric Environment 1993, 27B, 371.

19. Gery, M. W.; Whitten, G. Z.; Killus, J. P.; Dodge, M. C.; J. Geophys. Res., 94, D10, 12925.

20. The Updated OZIPR and OZIPW (OZIPR with Windows Interface). Volumes 1 and 2. User's Guide. Gery, M. W.; US EPA Order n^o 2D0848 NASA; North Carolina, 1995.

21. A intensidade total de radiação fora da atmosfera terrestre é a quantidade total de radiação por unidade de área normal à direção de propagação da luz. Esta quantidade tem um valor médio de 1368 W m^-2. Do ponto de vista fotoquímico, é mais útil expressar o fluxo solar para um certo intervalo de comprimento de onda. Quando a radiação passa através da atmosfera, ela é modificada na sua intensidade e composição espectral pela absorção e dispersão devido a gases e partículas. Neste artigo, o termo fluxo solar actínico é usado para a quantidade de radiação que atinge a superfície terrestre por unidade de área e para intervalos de 5 nm de comprimento de onda. As unidades na figura 2 são fotons cm^-2 s^-1.

22. Demerjiam, K.; Schere, K. L.; Peterson, J. T.; Adv. Environ. Sci. Technol. 1980, 10, 369.

23. JSPECTRA: UNC Solar Radiation Models; Jeffries, H. E.; Relatório final para a US EPA, November 1991, North Carolina 1991.

24. Fonseca, A. de A.; Campinho, M. A.; Corrêa, S. M.; Arbilla, G.; Quím. Nova 1996, 19, 108.

25. Fonseca, A. de A.; Campinho, M. A.; Corrêa, S. M.; Arbilla, G.; Quím. Nova 1996, 19, 361.

26. Hirst, D. M.; A Computational Approach to Chemistry, Blackwell Scientific Publications, Oxford 1990.

27. Steinfeld, J. I.; Francisco, J. S.; Hase, W. L.; Chemical Kinetics and Dynamics, Prentice Hall, New Jersey 1989.

28. Gear, C. W.; Numerical Initial Value Problems in Ordinary Differential Equations, Prentice-Hall, Englewood Cliffs, New Jersey, USA 1971.

29. Kaps, P.; Rentrop, P.; Numer. Math. 1979, 33, 55.

30. Hindmarsh, A. C.; ACM Signum Newsletter 1980, 15, 10.

31. Turányi, T.; Computers Chem. 1990, 14, 253.

32. Atkinson, R.; Baulch, D. L.; Cox, R. A.; Hampson Jr., R. F.; Kerr, J. A.; Troe, J.; J. Phys. Chem. Ref. Data 1992, 21, 1125.

33. Stein, A. F.; Toselli, B.; Atmospheric Environment, no prelo.

34. Conforme relatado na Referência 4, de 2 a 5 de julho de 1985, foram observados máximos de 2-3 ppb de PAN. Contudo, na maior parte dos dias (entre 29 de junho e 8 de julho) os valores de PAN foram menores do que 0,5 ppb e muitas vezes não atingiram o limite de detecção do aparelho (0,07 ppb). As concentrações medidas para ozônio foram £ 0,5 ppb.

35. Durante o inverno de 1996 foram registrados dias com alta poluição, provavelmente devido a uma combinação de fenômenos climáticos e problemas urbanos (engarrafamentos e obras). Ver, por exemplo, os artigos: "Dias de Poluição Máxima" ; Jornal do Brasil, 19 de junho de 1996 e " O Inverno da Inversão"; Jornal do Brasil; 20 de junho de 1996.

Apêndice

Datas de Publicação

• Publicação nesta coleção
  15 Set 2008
• Data do Fascículo
  Jun 1997

Histórico

• Aceito
  16 Out 1996
• Recebido
  07 Maio 1996