Emissões atmosféricas e implicações potenciais sobre a biota terrestre devido as atividades antrópicas na Baía do Almirantado/ Ilha Rei George Antártica

Neste trabalho foi realizado um estudo sobre o impacto sobre a atmosfera e a biota terrestre devido às emissões antrópicas na Baia do Almirantado/Ilha Rei George Antártica. Foram monitoradas as emissões dos compostos orgânicos voláteis e semi-voláteis, tanto nas fontes emissoras como no entorno da Estação Antártica Comandante Ferraz e estimadas as emissões dos navios, dos geradores a diesel e da incineração de lixo. Na avaliação do entorno, coletaram-se amostras de ar, neve e penas de aves. Com os resultados das emissões, do estudo topográfico e da meteorologia, realizou-se uma modelagem de plumas gaussiana para avaliar os impactos. Quatro cenários foram avaliados: dois com a presença dos navios NApOc Ary Rongel e Maximiano apresentaram concentrações máximas de até 356 g m-3 de COV e 18 g m-3 de material particulado, enquanto os demais, sem a presença dos navios, apenas considerando as estações de pesquisa EACF e Arctowski, apresentaram concentrações máximas de até 2,5 g m-3 de COV e 1,3 g m-3 de material particulado. Amostras de COV coletadas foram compatíveis com o cenário mais crítico. O estudo de correlação para carbonilas e HPA atmosférico e Carbono elementar e HPA, depositados em neve, apontaram a EACF como a principal fonte de emissão. As concentrações de levoglucosano detectadas a aproximadamente 2 km da EACF apontaram para a prática de incineração de lixo da EACF. Todas as áreas de interesse biológicos, anteriormente mapeados, dentro da AAEG, são vulneráveis às emissões antrópicas, como sugeriu o modelo de dispersão e a sobreposição dos resultados encontrados.

Mathematical models for the prediction of temperature distributions resulting from the discharge of heated water into large bodies of water

Mathematical models for heated water outfalls were developed for three flow regions. Near the source, the subsurface discharge into a stratified ambient water issuing from a row of buoyant jets was solved with the jet interference included in the analysis. The analysis of the flow zone close to and at intermediate distances from a surface buoyant jet was developed for the two-dimensional and axisymmetric cases. Far away from the source, a passive dispersion model was solved for a two dimensional situation taking into consideration the effects of shear current and vertical changes in diffusivity. A significant result from the surface buoyant jet analysis is the ability to predict the onset and location of an internal hydraulic jump. Prediction can be made simply from the knowledge of the source Froude number and a dimensionless surface exchange coefficient. Parametric computer programs of the above models are also developed as a part of this study. This report was submitted in fulfillment of Contract No. 14-12-570 under the sponsorship of the Federal Water Quality Administration.

Avaliação do impacto radiológico atmosférico de uma unidade de mineração e beneficiamento de urânio

Este trabalho teve como objetivo avaliar o impacto radiológico atmosférico da Unidade de Concentrado de Urânio URA, Caetité, BA, através da modelagem da dispersão de radionuclídeos e a estimativa da dose efetiva anual (em mSv.ano-1). Para tal, utilizou-se o programa MILDOS-AREA que foi desenvolvido pelo Argonne National Laboratory (ANL) em conjunto com a U.S. Nuclear Regulatory Commission (USNRC), para avaliar impacto radiológico ambiental atmosférico nas instalações de mineração e beneficiamento de urânio. O incremento de dose efetiva anual para três grupos críticos hipotéticos e oito grupos populacionais reais foi estimado com base na medida de fluxos de radônio e na estimativa das concentrações de radionuclídeos em particulados no ar dos principais termos fontes da URA (cava da mina, depósito de estéril e britador). Paralelamente, as medidas de concentração de radônio e taxa de kerma no ar, reportadas nos relatórios dos programas de monitoração ambiental pré-operacional (PMAPO) e operacional (PMAO) da URA, foram avaliadas. Os valores de dose efetiva anual estimados para os grupos críticos hipotéticos variaram de 1,78E-02 a 2,10E-02 mSv.ano-1, enquanto que para os grupos populacionais, variaram de 7,49E-05 a 1,56E-02 mSv.ano-1. A maior contribuição para o incremento da dose foi devida a inalação do radônio, sendo responsável por quase a totalidade da dose efetiva anual estimada. A média da concentração de atividade de radônio no entorno da URA foi 137,21 Bq m-3 e não sendo observada diferenças significativas entre as concentrações de radônio reportadas nos programas de monitoramento ambiental pré-operacional (valores de background) e operacional. Os valores médios de taxa de kerma no ar no entorno da URA foram de 0,136 μGy h-1. No entanto, em todos os pontos de monitoramento, os valores reportados no programa operacional foram inferiores aos valores reportados no programa pré-operacional (background), o que sugere problemas de medidas ou de coleta de dados durante a realização deste programa. O operador da URA utilizou para avaliação de impacto radiológico atmosférico, resultados apresentados em seus relatórios finais de análise de segurança (RFAS), um modelo próprio de simulação de dispersão, denominado Impacto Ambiental Radiológico (IAR7). Uma comparação entre o MILDOS-AREA e o IAR7, utilizando os mesmos parâmetros de entrada reportados no RFAS sugere que o IAR7 subestimou as concentrações de radônio no ar para os grupos críticos hipotéticos. Os resultados de simulação com o MILDOS-AREA mostram que as doses efetivas estimadas para os grupos críticos hipotéticos são inferiores a 0,3 mSv.ano-1 que é a restrição de dose estabelecida pela Comissão Nacional de Energia Nuclear. Recomenda-se que o código MILDOS-AREA seja utilizado no Brasil, para fins de licenciamento e controle, tendo em vista que o mesmo é um código validado e já utilizado em outros países para avaliar impacto radiológico ambiental atmosférico em instalações de mineração e beneficiamento de urânio

Coupling between cavity mode and heavy-hole exciton and light-hole exciton in semiconductor microcavity

The reflectivity spectra at different incident angles of semiconductor microcavity having heavy-hole exciton and light-hole exciton are calculated ly transfer matrix method using the linear dispersion model. Meanwhile we calculate the energy of three cavity polaritons at different incident angles formed by the coupling between cavity mode and the two exciton modes using the three harmonic oscillators coupling model, and the weights of cavity mode and the two exciton modes in the three cavity polaritons. The results indicate that there is obvious anticross between the high energy cavity polariton and the two low energy cavity polaritons with increasing incident angles, and the weights of three modes(cavity mode, heavy-hole exciton mode and light-hole exciton mode) in the three cavity polaritons increase or decrease.

Risk of Exposure to Airborne Ambrosia Pollen from Local and Distant Sources in Europe – an Example from Denmark

Background. Ambrosia artemisiifolia L. is a noxious invasive alien species in Europe. It is an important aeroallergen and millions of people are exposed to its pollen. Objective. The main aim of this study is to show that atmospheric concentrations of Ambrosia pollen recorded in Denmark can be derived from local or more distant sources. Methods. This was achieved by using a combination of pollen measurements, air mass trajectory calculations using the HYPLIT model and mapping all known Ambrosia locations in Denmark and relating them to land cover types. Results. The annual pollen index recorded in Copenhagen during a 15-year period varied from a few pollen grains to more than 100. Since 2005, small quantities of Ambrosia pollen has been observed in the air every year. We have demonstrated, through a combination of Lagrangian back-trajectory calculations and atmospheric pollen measurements, that pollen arrived in Denmark via long-distance transport from centres of Ambrosia infection, such as the Pannonian Plain and Ukraine. Combining observations with results from a local scale dispersion model show that it is possible that Ambrosia pollen could be derived from local sources identified within Denmark. Conclusions. The high allergenic capacity of Ambrosia pollen means that only small amounts of pollen are relevant for allergy sufferers, and just a few plants will be sufficient to produce enough pollen to affect pollen allergy sufferers within a short distance from the source. It is necessary to adopt control measures to restrict Ambrosia numbers. Recommendations for the removal of all Ambrosia plants can effectively reduce the amount of local pollen, as long as the population of Ambrosia plants is small.

Inventories and concentration profiles of 137Cs in undisturbed soils in the northeast of Buenos Aires Province, Argentina

Inventories and vertical distribution of (137)Cs were determined in La Plata region undisturbed soils, Argentina. A mean inventory value of 891 ± 220 Bq/m(2) was established, which is compatible with the values expected from atmospheric weapon tests fallout. The study was complemented with pH, organic carbon fraction, texture and mineralogical soil analyses. Putting together Southern Hemisphere (137)Cs inventory data, it is possible to correlate these data with the mean annual precipitations. The large differences in (137)Cs concentration profiles were attributed to soil properties, especially the clay content and the pH values. A convection-dispersion model with irreversible retention was used to fit the activity concentration profiles. The obtained effective diffusion coefficient and effective convection velocity parameters values were in the range from 0.2 cm(2)/y to 0.4 cm(2)/y and from 0.23 cm/y to 0.43 cm/y, respectively. These data are in agreement with values reported in literature. In general, with the growth of clay content in the soil, there was an increase in the transfer rate from free to bound state. Finally, the highest transfer rate from free to bound state was obtained for soil pH value equal to 8.

Modeling reaction kinetics and mass transfer in ozonation in water solutions

This dissertation is based on four articles dealing with modeling of ozonation. The literature part of this considers some models for hydrodynamics in bubble column simulation. A literature review of methods for obtaining mass transfer coefficients is presented. The methods presented to obtain mass transfer are general models and can be applied to any gas-liquid system. Ozonation reaction models and methods for obtaining stoichiometric coefficients and reaction rate coefficients for ozonation reactions are discussed in the final section of the literature part. In the first article, ozone gas-liquid mass transfer into water in a bubble column was investigated for different pH values. A more general method for estimation of mass transfer and Henry’s coefficient was developed from the Beltrán method. The ozone volumetric mass transfer coefficient and the Henry’s coefficient were determined simultaneously by parameter estimation using a nonlinear optimization method. A minor dependence of the Henry’s law constant on pH was detected at the pH range 4 - 9. In the second article, a new method using the axial dispersion model for estimation of ozone self-decomposition kinetics in a semi-batch bubble column reactor was developed. The reaction rate coefficients for literature equations of ozone decomposition and the gas phase dispersion coefficient were estimated and compared with the literature data. The reaction order in the pH range 7-10 with respect to ozone 1.12 and 0.51 the hydroxyl ion were obtained, which is in good agreement with literature. The model parameters were determined by parameter estimation using a nonlinear optimization method. Sensitivity analysis was conducted using object function method to obtain information about the reliability and identifiability of the estimated parameters. In the third article, the reaction rate coefficients and the stoichiometric coefficients in the reaction of ozone with the model component p-nitrophenol were estimated at low pH of water using nonlinear optimization. A novel method for estimation of multireaction model parameters in ozonation was developed. In this method the concentration of unknown intermediate compounds is presented as a residual COD (chemical oxygen demand) calculated from the measured COD and the theoretical COD for the known species. The decomposition rate of p-nitrophenol on the pathway producing hydroquinone was found to be about two times faster than the p-nitrophenol decomposition rate on the pathway producing 4- nitrocatechol. In the fourth article, the reaction kinetics of p-nitrophenol ozonation was studied in a bubble column at pH 2. Using the new reaction kinetic model presented in the previous article, the reaction kinetic parameters, rate coefficients, and stoichiometric coefficients as well as the mass transfer coefficient were estimated with nonlinear estimation. The decomposition rate of pnitrophenol was found to be equal both on the pathway producing hydroquinone and on the path way producing 4-nitrocathecol. Comparison of the rate coefficients with the case at initial pH 5 indicates that the p-nitrophenol degradation producing 4- nitrocathecol is more selective towards molecular ozone than the reaction producing hydroquinone. The identifiability and reliability of the estimated parameters were analyzed with the Marcov chain Monte Carlo (MCMC) method. @All rights reserved. No part of the publication may be reproduced, stored in a retrieval system, or transmitted, in any form or by any means, electronic, mechanical, photocopying, recording, or otherwise, without the prior permission of the author.

Modélisation bayésienne avec des splines du comportement moyen d'un échantillon de courbes

Cette thèse porte sur l'analyse bayésienne de données fonctionnelles dans un contexte hydrologique. L'objectif principal est de modéliser des données d'écoulements d'eau d'une manière parcimonieuse tout en reproduisant adéquatement les caractéristiques statistiques de celles-ci. L'analyse de données fonctionnelles nous amène à considérer les séries chronologiques d'écoulements d'eau comme des fonctions à modéliser avec une méthode non paramétrique. Dans un premier temps, les fonctions sont rendues plus homogènes en les synchronisant. Ensuite, disposant d'un échantillon de courbes homogènes, nous procédons à la modélisation de leurs caractéristiques statistiques en faisant appel aux splines de régression bayésiennes dans un cadre probabiliste assez général. Plus spécifiquement, nous étudions une famille de distributions continues, qui inclut celles de la famille exponentielle, de laquelle les observations peuvent provenir. De plus, afin d'avoir un outil de modélisation non paramétrique flexible, nous traitons les noeuds intérieurs, qui définissent les éléments de la base des splines de régression, comme des quantités aléatoires. Nous utilisons alors le MCMC avec sauts réversibles afin d'explorer la distribution a posteriori des noeuds intérieurs. Afin de simplifier cette procédure dans notre contexte général de modélisation, nous considérons des approximations de la distribution marginale des observations, nommément une approximation basée sur le critère d'information de Schwarz et une autre qui fait appel à l'approximation de Laplace. En plus de modéliser la tendance centrale d'un échantillon de courbes, nous proposons aussi une méthodologie pour modéliser simultanément la tendance centrale et la dispersion de ces courbes, et ce dans notre cadre probabiliste général. Finalement, puisque nous étudions une diversité de distributions statistiques au niveau des observations, nous mettons de l'avant une approche afin de déterminer les distributions les plus adéquates pour un échantillon de courbes donné.

Estudio del impacto ambiental del uso del bagazo como fuente de energía en centrales azucareros en Cuba. Estudio de caso "Melanio Hernández"

This work was carried out to study the environmental impacts produced by power generation from biomass (bagasse) in sugar mills in Cuba. For this purpose, with the collaboration of the Center for Energy and Industrial Processes (CEEPI), the University Center of Sancti Spiritus and using different research methods and techniques, conducted an environmental survey of the area. Gaseous emissions were characterized and the suspended solid particles, allowing knowing that these concentrations do not exceed the maximum emission limit set by the cubana Standard. In addition, the dispersion model applied DISPER allowed us to obtain information in crisis conditions of air emissions from sugar mill and distillery associated and concluded that emissions from the distillery are the pollutants that contribute more the atmosphere. The correlating emissions with respiratory diseases, (acute respiratory infections or subacute (IRA) and ASMA), concurrent Asthma Crisis (CAB) is the most affected. The calculated costs associated in these diseases, amounting to $ 119 599.23 per year. In order to minimize the negative alternatives are proposed to be implemented in industry and community.

Solute movement through intact columns of cryoturbated Upper Chalk

Cryoturbated Upper Chalk is a dichotomous porous medium wherein the intra-fragment porosity provides water storage and the inter-fragment porosity provides potential pathways for relatively rapid flow near saturation. Chloride tracer movement through 43 cm long and 45 cm diameter undisturbed chalk columns was studied at water application rates of 0.3, 1.0, and 1.5 cm h(-1). Microscale heterogeneity in effluent was recorded using a grid collection system consisting of 98 funnel-shaped cells each 3.5 cm in diameter. The total porosity of the columns was 0.47 +/- 0.02 m(3) m(-3), approximately 13% of pores were >15 mu m diameter, and the saturated hydraulic conductivity was 12.66 +/- 1.31 m day(-1). Although the column remained unsaturated during the leaching even at all application rates, proportionate flow through macropores increased as the application rate decreased. The number of dry cells (with 0 ml of effluent) increased as application rate decreased. Half of the leachate was collected from 15, 19 and 22 cells at 0.3, 1.0, 1.5 cm h(-1) application rates respectively. Similar breakthrough curves (BTCs) were obtained at all three application rates when plotted as a function of cumulative drainage, but they were distinctly different when plotted as a function of time. The BTCs indicate that the columns have similar drainage requirement irrespective of application rates, as the rise to the maxima (C/C-o) is almost similar. However, the time required to achieve that leaching requirement varies with application rates, and residence time was less in the case of a higher application rate. A two-region convection-dispersion model was used to describe the BTCs and fitted well (r(2) = 0.97-0-99). There was a linear relationship between dispersion coefficient and pore water velocity (correlation coefficient r = 0.95). The results demonstrate the microscale heterogeneity of hydrodynamic properties in the Upper Chalk. Copyright (C) 2007 John Wiley & Sons, Ltd.

Quantifying uncertainty in critical loads: (B) acidity mass balance critical loads on a sensitive site

This paper reports an uncertainty analysis of critical loads for acid deposition for a site in southern England, using the Steady State Mass Balance Model. The uncertainty bounds, distribution type and correlation structure for each of the 18 input parameters was considered explicitly, and overall uncertainty estimated by Monte Carlo methods. Estimates of deposition uncertainty were made from measured data and an atmospheric dispersion model, and hence the uncertainty in exceedance could also be calculated. The uncertainties of the calculated critical loads were generally much lower than those of the input parameters due to a "compensation of errors" mechanism - coefficients of variation ranged from 13% for CLmaxN to 37% for CL(A). With 1990 deposition, the probability that the critical load was exceeded was > 0.99; to reduce this probability to 0.50, a 63% reduction in deposition is required; to 0.05, an 82% reduction. With 1997 deposition, which was lower than that in 1990, exceedance probabilities declined and uncertainties in exceedance narrowed as deposition uncertainty had less effect. The parameters contributing most to the uncertainty in critical loads were weathering rates, base cation uptake rates, and choice of critical chemical value, indicating possible research priorities. However, the different critical load parameters were to some extent sensitive to different input parameters. The application of such probabilistic results to environmental regulation is discussed.

A new method for evaluating regional air quality forecasts

A Kriging interpolation method is combined with an object-based evaluation measure to assess the ability of the UK Met Office's dispersion and weather prediction models to predict the evolution of a plume of tracer as it was transported across Europe. The object-based evaluation method, SAL, considers aspects of the Structure, Amplitude and Location of the pollutant field. The SAL method is able to quantify errors in the predicted size and shape of the pollutant plume, through the structure component, the over- or under-prediction of the pollutant concentrations, through the amplitude component, and the position of the pollutant plume, through the location component. The quantitative results of the SAL evaluation are similar for both models and close to a subjective visual inspection of the predictions. A negative structure component for both models, throughout the entire 60 hour plume dispersion simulation, indicates that the modelled plumes are too small and/or too peaked compared to the observed plume at all times. The amplitude component for both models is strongly positive at the start of the simulation, indicating that surface concentrations are over-predicted by both models for the first 24 hours, but modelled concentrations are within a factor of 2 of the observations at later times. Finally, for both models, the location component is small for the first 48 hours after the start of the tracer release, indicating that the modelled plumes are situated close to the observed plume early on in the simulation, but this plume location error grows at later times. The SAL methodology has also been used to identify differences in the transport of pollution in the dispersion and weather prediction models. The convection scheme in the weather prediction model is found to transport more pollution vertically out of the boundary layer into the free troposphere than the dispersion model convection scheme resulting in lower pollutant concentrations near the surface and hence a better forecast for this case study.

Tracer concentration profiles measured in central London as part of the REPARTEE campaign

There have been relatively few tracer experiments carried out that have looked at vertical plume spread in urban areas. In this paper we present results from two tracer (cyclic perfluorocarbon) experiments carried out in 2006 and 2007 in central London centred on the BT Tower as part of the REPARTEE (Regent’s Park and Tower Environmental Experiment) campaign. The height of the tower gives a unique opportunity to study vertical dispersion profiles and transport times in central London. Vertical gradients are contrasted with the relevant Pasquill stability classes. Estimation of lateral advection and vertical mixing times are made and compared with previous measurements. Data are then compared with a simple operational dispersion model and contrasted with data taken in central London as part of the DAPPLE campaign. This correlates dosage with non-dimensionalised distance from source. Such analyses illustrate the feasibility of the use of these empirical correlations over these prescribed distances in central London.

A study of the arrival over the United Kingdom in April 2010 of the Eyjafjallajökull ash cloud using ground-based lidar and numerical simulations

We make a qualitative and quantitative comparison of numericalsimulations of the ashcloud generated by the eruption of Eyjafjallajökull in April2010 with ground-basedlidar measurements at Exeter and Cardington in southern England. The numericalsimulations are performed using the Met Office’s dispersion model, NAME (Numerical Atmospheric-dispersion Modelling Environment). The results show that NAME captures many of the features of the observed ashcloud. The comparison enables us to estimate the fraction of material which survives the near-source fallout processes and enters into the distal plume. A number of simulations are performed which show that both the structure of the ashcloudover southern England and the concentration of ash within it are particularly sensitive to the height of the eruption column (and the consequent estimated mass emission rate), to the shape of the vertical source profile and the level of prescribed ‘turbulent diffusion’ (representing the mixing by the unresolved eddies) in the free troposphere with less sensitivity to the timing of the start of the eruption and the sedimentation of particulates in the distal plume.

Quantifying volcanic ash dispersal and impact of the Campanian Ignimbrite super-eruption

1] We apply a novel computational approach to assess, for the first time, volcanic ash dispersal during the Campanian Ignimbrite (Italy) super-eruption providing insights into eruption dynamics and the impact of this gigantic event. The method uses a 3D time-dependent computational ash dispersion model, a set of wind fields, and more than 100 thickness measurements of the CI tephra deposit. Results reveal that the CI eruption dispersed 250–300 km3 of ash over ∼3.7 million km2. The injection of such a large quantity of ash (and volatiles) into the atmosphere would have caused a volcanic winter during the Heinrich Event 4, the coldest and driest climatic episode of the Last Glacial period. Fluorine-bearing leachate from the volcanic ash and acid rain would have further affected food sources and severely impacted Late Middle-Early Upper Paleolithic groups in Southern and Eastern Europe