Sibilância em crianças e adolescentes vizinhos a uma indústria petroquímica no Rio Grande do Norte, Brasil

OBJETIVO: Verificar a associação entre relato de sibilância em crianças e adolescentes e o local de residência em relação à dispersão dos poluentes atmosféricos emitidos pelo Pólo Petroquímico (PPQ) de Guamaré (RN). MÉTODOS: Estudo transversal de relato de sibilância em crianças e adolescentes de 0 a 14 anos de idade, residentes no entorno do PPQ de Guamaré, em 2006. Foi utilizado o questionário padronizado do International Study of Asthma and Allergies in Childhood, acrescido de questões relativas ao tabagismo, renda, moradia e escolaridade. Concentrações diárias de PM10, PM2,5, carbono grafítico, SO2, NO2, O3, benzeno, tolueno e xilenos foram medidas em uma estação de monitoramento fixa. As comunidades residentes na área de influência das emissões do PPQ foram classificadas, segundo a direção preferencial dos ventos, em expostas e de referência. RESULTADOS: Participaram do estudo 209 crianças e adolescentes. As concentrações médias diárias dos poluentes monitorados mantiveram-se abaixo dos limites estabelecidos nos padrões de qualidade do ar. A prevalência de sibilos nos últimos 12 meses foi de 27,3%. Associações estatisticamente significantes com sibilos nos últimos 12 meses foram verificadas mesmo após ajustamentos para comunidades expostas [razão de chances (odds ratio, ORajust) = 2,01; intervalo de confiança de 95% (IC95%) 1,01-4,01], gênero masculino (ORajust = 2,50; IC95% 1,21-5,18) e idade de 0 a 6 anos (ORajust = 5,00; IC95% 2,41-10,39). CONCLUSÃO: Mesmo em baixas concentrações de poluentes atmosféricos, a ocorrência de sintomas respiratórios em crianças e adolescentes nas comunidades no entorno de um PPQ esteve associada a residência na direção preferencial dos ventos, mostrando-se mais vulnerável o grupo de pré-escolares do gênero masculino.

Caracterização química e toxicológica da água produzida descartada em plataformas de óleo e gás na costa brasileira e seu comportamento dispersivo no mar

The production of oil and gas is usually accompanied by the production of water, also known as produced water. Studies were conducted in platforms that discharge produced water in the Atlantic Ocean due to oil and gas production by Petrobras from 1996 to 2006 in the following basins: Santos (Brazilian south region), Campos (Brazilian southeast region) and Ceara (Brazilian northeast region). This study encompasses chemical composition, toxicological effects, discharge volumes, and produced water behavior after releasing in the ocean, including dispersion plumes modeling and monitoring data of the marine environment. The concentration medians for a sampling of 50 samples were: ammonia (70 mg L-1), boron (1.3 mg L1), iron (7.4 mg L-1), BTEX (4.6 mg L-1), PAH (0.53 mg L-1), TPH (28 mg L-1); phenols (1.3 mg L-1) and radioisotopes (0.15 Bq L-1 for 226Ra and 0.09 Bq L-1 for 228Ra). The concentrations of the organic and inorganic parameters observed for the Brazilian platforms were similar to the international reference data for the produced water in the North Sea and in other regions of the world. It was found significant differences in concentrations of the following parameters: BTEX (p<0.0001), phenols (p=0.0212), boron (p<0.0001), iron (p<0.0001) and toxicological response in sea urchin Lytechinus variegatus (p<0.0001) when considering two distinguished groups, platforms from southeast and northeast Region (PCR-1). Significant differences were not observed among the other parameters. In platforms with large gas production, the monoaromatic concentrations (BTEX from 15.8 to 21.6 mg L-1) and phenols (from 2 to 83 mg L-1) were higher than in oil plataforms (median concentrations of BTEX were 4.6 mg L-1 for n=53, and of phenols were 1.3 mg L-1 for n=46). It was also conducted a study about the influence of dispersion plumes of produced water in the vicinity of six platforms of oil and gas production (P-26, PPG-1, PCR-1, P-32, SS-06), and in a hypothetical critical scenario using the chemical characteristics of each effluent. Through this study, using CORMIX and CHEMMAP models for dispersion plumes simulation of the produced water discharges, it was possible to obtain the dilution dimension in the ocean after those discharges. The dispersion plumes of the produced water modelling in field vicinity showed dilutions of 700 to 900 times for the first 30-40 meters from the platform PCR-1 discharge point; 100 times for the platform P-32, with 30 meters of distance; 150 times for the platform P-26, with 40 meters of distance; 100 times for the platform PPG-1, with 130 meters of distance; 280 to 350 times for the platform SS-06, with 130 meters of distance, 100 times for the hypothetical critical scenario, with the 130 meters of distance. The dilutions continue in the far field, and with the results of the simulations, it was possible to verify that all the parameters presented concentrations bellow the maximum values established by Brazilian legislation for seawater (CONAMA 357/05 - Class 1), before the 500 meters distance of the discharge point. These results were in agreement with the field measurements. Although, in general results for the Brazilian produced water presented toxicological effects for marine organisms, it was verified that dilutions of 100 times were sufficient for not causing toxicological responses. Field monitoring data of the seawater around the Pargo, Pampo and PCR-1 platforms did not demonstrate toxicity in the seawater close to these platforms. The results of environmental monitoring in seawater and sediments proved that alterations were not detected for environmental quality in areas under direct influence of the oil production activities in the Campos and Ceara Basin, as according to results obtained in the dispersion plume modelling for the produced water discharge

Caracterização química e toxicológica da água produzida descartada em plataformas de óleo e gás na costa brasileira e seu comportamento dispersivo no mar

The production of oil and gas is usually accompanied by the production of water, also known as produced water. Studies were conducted in platforms that discharge produced water in the Atlantic Ocean due to oil and gas production by Petrobras from 1996 to 2006 in the following basins: Santos (Brazilian south region), Campos (Brazilian southeast region) and Ceara (Brazilian northeast region). This study encompasses chemical composition, toxicological effects, discharge volumes, and produced water behavior after releasing in the ocean, including dispersion plumes modeling and monitoring data of the marine environment. The concentration medians for a sampling of 50 samples were: ammonia (70 mg L-1), boron (1.3 mg L1), iron (7.4 mg L-1), BTEX (4.6 mg L-1), PAH (0.53 mg L-1), TPH (28 mg L-1); phenols (1.3 mg L-1) and radioisotopes (0.15 Bq L-1 for 226Ra and 0.09 Bq L-1 for 228Ra). The concentrations of the organic and inorganic parameters observed for the Brazilian platforms were similar to the international reference data for the produced water in the North Sea and in other regions of the world. It was found significant differences in concentrations of the following parameters: BTEX (p<0.0001), phenols (p=0.0212), boron (p<0.0001), iron (p<0.0001) and toxicological response in sea urchin Lytechinus variegatus (p<0.0001) when considering two distinguished groups, platforms from southeast and northeast Region (PCR-1). Significant differences were not observed among the other parameters. In platforms with large gas production, the monoaromatic concentrations (BTEX from 15.8 to 21.6 mg L-1) and phenols (from 2 to 83 mg L-1) were higher than in oil plataforms (median concentrations of BTEX were 4.6 mg L-1 for n=53, and of phenols were 1.3 mg L-1 for n=46). It was also conducted a study about the influence of dispersion plumes of produced water in the vicinity of six platforms of oil and gas production (P-26, PPG-1, PCR-1, P-32, SS-06), and in a hypothetical critical scenario using the chemical characteristics of each effluent. Through this study, using CORMIX and CHEMMAP models for dispersion plumes simulation of the produced water discharges, it was possible to obtain the dilution dimension in the ocean after those discharges. The dispersion plumes of the produced water modelling in field vicinity showed dilutions of 700 to 900 times for the first 30-40 meters from the platform PCR-1 discharge point; 100 times for the platform P-32, with 30 meters of distance; 150 times for the platform P-26, with 40 meters of distance; 100 times for the platform PPG-1, with 130 meters of distance; 280 to 350 times for the platform SS-06, with 130 meters of distance, 100 times for the hypothetical critical scenario, with the 130 meters of distance. The dilutions continue in the far field, and with the results of the simulations, it was possible to verify that all the parameters presented concentrations bellow the maximum values established by Brazilian legislation for seawater (CONAMA 357/05 - Class 1), before the 500 meters distance of the discharge point. These results were in agreement with the field measurements. Although, in general results for the Brazilian produced water presented toxicological effects for marine organisms, it was verified that dilutions of 100 times were sufficient for not causing toxicological responses. Field monitoring data of the seawater around the Pargo, Pampo and PCR-1 platforms did not demonstrate toxicity in the seawater close to these platforms. The results of environmental monitoring in seawater and sediments proved that alterations were not detected for environmental quality in areas under direct influence of the oil production activities in the Campos and Ceara Basin, as according to results obtained in the dispersion plume modelling for the produced water discharge

Avaliação da dispersão de efluente líquido de refinaria de petróleo

Neste trabalho são apresentadas duas técnicas para a avaliação da dispersão, num corpo d'água receptor, do efluente líquido de uma refinaria de petróleo. Assim, a pluma de dispersão foi caracterizada por medidas em campo de condutividade elétrica e por simulação computacional (simulador Cormix). Como caso de estudo, escolheu-se uma refinaria de petróleo cujo efluente é lançado no rio Atibaia (Paulínia/SP). O comportamento do efluente foi avaliado em um trecho de 1000 m após o ponto de lançamento. Os resultados demonstraram que a medição da condutividade elétrica é uma técnica adequada para a avaliação da dispersão de efluentes líquidos de refinaria de petróleo, pois apresentam alta condutividade elétrica e, com isso, há um forte contraste entre os valores do efluente e do rio. Além disso, outros parâmetros de qualidade da água do rio seguiram comportamento de dispersão semelhante ao da condutividade. A pluma de dispersão gerada pelo simulador computacional apresentou uma elevada concordância com os dados obtidos em campo. Nesse sentido, a simulação computacional pode ser uma ferramenta útil para a avaliação da dispersão do efluente considerando-se cenários hipotéticos, e para projetos de emissários.

Geostatistical Mapping of Outfall Plume Dispersion Data gathered with an autonomous underwater vehicle

The main purpose of this study was to examine the applicability of geostatistical modeling to obtain valuable information for assessing the environmental impact of sewage outfall discharges. The data set used was obtained in a monitoring campaign to S. Jacinto outfall, located off the Portuguese west coast near Aveiro region, using an AUV. The Matheron’s classical estimator was used the compute the experimental semivariogram which was fitted to three theoretical models: spherical, exponential and gaussian. The cross-validation procedure suggested the best semivariogram model and ordinary kriging was used to obtain the predictions of salinity at unknown locations. The generated map shows clearly the plume dispersion in the studied area, indicating that the effluent does not reach the near by beaches. Our study suggests that an optimal design for the AUV sampling trajectory from a geostatistical prediction point of view, can help to compute more precise predictions and hence to quantify more accurately dilution. Moreover, since accurate measurements of plume’s dilution are rare, these studies might be very helpful in the future for validation of dispersion models.

Geostatistical Prediction of Ocean Outfall Plume Characteristics Based on an Autonomous Underwater Vehicle

Geostatistics has been successfully used to analyze and characterize the spatial variability of environmental properties. Besides giving estimated values at unsampled locations, it provides a measure of the accuracy of the estimate, which is a significant advantage over traditional methods used to assess pollution. In this work universal block kriging is novelty used to model and map the spatial distribution of salinity measurements gathered by an Autonomous Underwater Vehicle in a sea outfall monitoring campaign, with the aim of distinguishing the effluent plume from the receiving waters, characterizing its spatial variability in the vicinity of the discharge and estimating dilution. The results demonstrate that geostatistical methodology can provide good estimates of the dispersion of effluents that are very valuable in assessing the environmental impact and managing sea outfalls. Moreover, since accurate measurements of the plume’s dilution are rare, these studies might be very helpful in the future to validate dispersion models.

Dispersion Modelling of Natural Radionuclides 238U, 232Th and 40K Released from Coal-fired Power Plants Operations

The aim of this work was to simulate the radionuclides dispersion in the surrounding area of a coal-fired power plant, operational during the last 25 years. The dispersion of natural radionuclides (236Ra, 232Th and 40K) was simulated by a Gaussian plume dispersion model with three different stability classes estimating the radionuclides concentration at ground level. Measurements of the environmen-tal activity concentrations were carried out by γ-spectrometry and compared with results from the air dispersion and deposition model which showed that the stabil-ity class D causes the dispersion to longer distances up to 20 km from the stacks.

Comparing the performance of geostatistical models with additional information from covariates for sewage plume characterization

In this work, kriging with covariates is used to model and map the spatial distribution of salinity measurements gathered by an autonomous underwater vehicle in a sea outfall monitoring campaign aiming to distinguish the effluent plume from the receiving waters and characterize its spatial variability in the vicinity of the discharge. Four different geostatistical linear models for salinity were assumed, where the distance to diffuser, the west-east positioning, and the south-north positioning were used as covariates. Sample variograms were fitted by the Mat`ern models using weighted least squares and maximum likelihood estimation methods as a way to detect eventual discrepancies. Typically, the maximum likelihood method estimated very low ranges which have limited the kriging process. So, at least for these data sets, weighted least squares showed to be the most appropriate estimation method for variogram fitting. The kriged maps show clearly the spatial variation of salinity, and it is possible to identify the effluent plume in the area studied. The results obtained show some guidelines for sewage monitoring if a geostatistical analysis of the data is in mind. It is important to treat properly the existence of anomalous values and to adopt a sampling strategy that includes transects parallel and perpendicular to the effluent dispersion.

Numerical methods and tangible interfaces for pollutant dispersion simulation

Dissertação apresentada para obtenção do Grau de Doutor em Engenharia do Ambiente, pela Universidade Nova de Lisboa, Faculdade de Ciências e Tecnologia

Application of a non isotropic turbulence model to stable atmospheric flows and dispersion over 3D topography

A non isotropic turbulence model is extended and applied to three dimensional stably stratified flows and dispersion calculations. The model is derived from the algebraic stress model (including wall proximity effects), but it retains the simplicity of the "eddy viscosity" concept of first order models. The "modified k-epsilon" is implemented in a three dimensional numerical code. Once the flow is resolved, the predicted velocity and turbulence fields are interpolated into a second grid and used to solve the concentration equation. To evaluate the model, various steady state numerical solutions are compared with small scale dispersion experiments which were conducted at the wind tunnel of Mitsubishi Heavy Industries, in Japan. Stably stratified flows and plume dispersion over three distinct idealized complex topographies (flat and hilly terrain) are studied. Vertical profiles of velocity and pollutant concentration are shown and discussed. Also, comparisons are made against the results obtained with the standard k-epsilon model.

Wind-tunnel modelling of dispersion from a scalar area source in urban-like roughness

A wind-tunnel study was conducted to investigate ventilation of scalars from urban-like geometries at neighbourhood scale by exploring two different geometries a uniform height roughness and a non-uniform height roughness, both with an equal plan and frontal density of λ p = λ f = 25%. In both configurations a sub-unit of the idealized urban surface was coated with a thin layer of naphthalene to represent area sources. The naphthalene sublimation method was used to measure directly total area-averaged transport of scalars out of the complex geometries. At the same time, naphthalene vapour concentrations controlled by the turbulent fluxes were detected using a fast Flame Ionisation Detection (FID) technique. This paper describes the novel use of a naphthalene coated surface as an area source in dispersion studies. Particular emphasis was also given to testing whether the concentration measurements were independent of Reynolds number. For low wind speeds, transfer from the naphthalene surface is determined by a combination of forced and natural convection. Compared with a propane point source release, a 25% higher free stream velocity was needed for the naphthalene area source to yield Reynolds-number-independent concentration fields. Ventilation transfer coefficients w T /U derived from the naphthalene sublimation method showed that, whilst there was enhanced vertical momentum exchange due to obstacle height variability, advection was reduced and dispersion from the source area was not enhanced. Thus, the height variability of a canopy is an important parameter when generalising urban dispersion. Fine resolution concentration measurements in the canopy showed the effect of height variability on dispersion at street scale. Rapid vertical transport in the wake of individual high-rise obstacles was found to generate elevated point-like sources. A Gaussian plume model was used to analyse differences in the downstream plumes. Intensified lateral and vertical plume spread and plume dilution with height was found for the non-uniform height roughness

Evaluating the structure and magnitude of the ash plume during the initial phase of the 2010 Eyjafjallajökull eruption using lidar observations and NAME simulations

The Eyjafjallajökull volcano in Iceland erupted explosively on 14 April 2010, emitting a plume of ash into the atmosphere. The ash was transported from Iceland toward Europe where mostly cloud-free skies allowed ground-based lidars at Chilbolton in England and Leipzig in Germany to estimate the mass concentration in the ash cloud as it passed overhead. The UK Met Office's Numerical Atmospheric-dispersion Modeling Environment (NAME) has been used to simulate the evolution of the ash cloud from the Eyjafjallajökull volcano during the initial phase of the ash emissions, 14–16 April 2010. NAME captures the timing and sloped structure of the ash layer observed over Leipzig, close to the central axis of the ash cloud. Relatively small errors in the ash cloud position, probably caused by the cumulative effect of errors in the driving meteorology en route, result in a timing error at distances far from the central axis of the ash cloud. Taking the timing error into account, NAME is able to capture the sloped ash layer over the UK. Comparison of the lidar observations and NAME simulations has allowed an estimation of the plume height time series to be made. It is necessary to include in the model input the large variations in plume height in order to accurately predict the ash cloud structure at long range. Quantitative comparison with the mass concentrations at Leipzig and Chilbolton suggest that around 3% of the total emitted mass is transported as far as these sites by small (<100 μm diameter) ash particles.

Dispersion of a point-source release of a passive scalar through an urban-like array for different wind directions

The dispersion of a point-source release of a passive scalar in a regular array of cubical, urban-like, obstacles is investigated by means of direct numerical simulations. The simulations are conducted under conditions of neutral stability and fully rough turbulent flow, at a roughness Reynolds number of Reτ = 500. The Navier–Stokes and scalar equations are integrated assuming a constant rate release from a point source close to the ground within the array. We focus on short-range dispersion, when most of the material is still within the building canopy. Mean and fluctuating concentrations are computed for three different pressure gradient directions (0◦ , 30◦ , 45◦). The results agree well with available experimental data measured in a water channel for a flow angle of 0◦ . Profiles of mean concentration and the three-dimensional structure of the dispersion pattern are compared for the different forcing angles. A number of processes affecting the plume structure are identified and discussed, including: (i) advection or channelling of scalar down ‘streets’, (ii) lateral dispersion by turbulent fluctuations and topological dispersion induced by dividing streamlines around buildings, (iii) skewing of the plume due to flow turning with height, (iv) detrainment by turbulent dispersion or mean recirculation, (v) entrainment and release of scalar in building wakes, giving rise to ‘secondary sources’, (vi) plume meandering due to unsteady turbulent fluctuations. Finally, results on relative concentration fluctuations are presented and compared with the literature for point source dispersion over flat terrain and urban arrays. Keywords Direct numerical simulation · Dispersion modelling · Urban array

Processes controlling atmospheric dispersion through city centres

We develop a process-based model for the dispersion of a passive scalar in the turbulent flow around the buildings of a city centre. The street network model is based on dividing the airspace of the streets and intersections into boxes, within which the turbulence renders the air well mixed. Mean flow advection through the network of street and intersection boxes then mediates further lateral dispersion. At the same time turbulent mixing in the vertical detrains scalar from the streets and intersections into the turbulent boundary layer above the buildings. When the geometry is regular, the street network model has an analytical solution that describes the variation in concentration in a near-field downwind of a single source, where the majority of scalar lies below roof level. The power of the analytical solution is that it demonstrates how the concentration is determined by only three parameters. The plume direction parameter describes the branching of scalar at the street intersections and hence determines the direction of the plume centreline, which may be very different from the above-roof wind direction. The transmission parameter determines the distance travelled before the majority of scalar is detrained into the atmospheric boundary layer above roof level and conventional atmospheric turbulence takes over as the dominant mixing process. Finally, a normalised source strength multiplies this pattern of concentration. This analytical solution converges to a Gaussian plume after a large number of intersections have been traversed, providing theoretical justification for previous studies that have developed empirical fits to Gaussian plume models. The analytical solution is shown to compare well with very high-resolution simulations and with wind tunnel experiments, although re-entrainment of scalar previously detrained into the boundary layer above roofs, which is not accounted for in the analytical solution, is shown to become an important process further downwind from the source.

The plume also rises: trajectories of pheromone plumes issuing from point sources in an orchard canopy at night

The trajectories of pheromone plumes in canopied habitats, such as orchards, have been little studied. We documented the capture of male navel orangeworm moths, Amyelois transitella, in female-baited traps positioned at 5 levels, from ground level to the canopy top, at approximately 6 m above ground, in almond orchards. Males were captured in similar proportions at all levels, suggesting that they do not favor a particular height during ranging flight. A 3-D sonic anemometer was used to establish patterns of wind flow and temperature at 6 heights from 2.08 to 6.65 m in an almond orchard with a 5 m high canopy, every 3 h over 72 h. The horizontal velocity of wind flow was highest above the canopy, where its directionality also was the most consistent. During the time of A. transitella mating (0300–0600), there was a net vertical displacement upward. Vertical buoyancy combined with only minor reductions in the distance that plumes will travel in the lower compared to the upper canopy suggest that the optimal height for release of pheromone from high-release-rate sources, such as aerosol dispensers (“puffers”), that are deployed at low densities (e.g., 3 per ha.) would be at mid or low in the canopy, thereby facilitating dispersion of disruptant throughout the canopy. Optimal placement of aerosol dispensers will vary with the behavioral ecology of the target pest; however, our results suggest that current protocols, which generally propose dispenser placement in the upper third of the canopy, should be reevaluated.