A wind tunnel study of dense gas dispersion in a neutral boundary layer over a rough surface

Measurements of the vertical entrainment velocity into two-dimensional dense gas plumes over fully rough surfaces were carried out as part of a co-operative research programme with wind tunnel facilities in the USA. This paper presents results obtained for neutral boundary layer conditions in the EnFlo wind tunnel at the University of Surrey; a companion paper treats the stable boundary layer case. Entrainment velocities, W_E, were deduced from the streamwise development of the concentration field, non-dimensionalised with respect to the friction velocity in the undisturbed flow, u^*, and correlated with the plume Richardson number, Ri^*. Results for Richardson numbers in the range Ri^*<15 were found to be well fitted by the empirical expression: W_E/u^*=0.65/(1+0.2Ri^*). Flow visualisation studies showed layered plume structures with a sharp upper interface at higher Richardson numbers and in this regime turbulent motion below the interface became progressively more intermittent as Ri^* increased. Measured turbulence levels collapsed within such high Richardson number plumes and flow and dispersion were significantly affected by molecular processes. Up-welling above the source was observed when the emission speed exceeded the approach flow friction velocity, though there was no clear evidence that this affected plume behaviour away from the immediate vicinity of the source.     

A wind tunnel study of dense gas dispersion in a stable boundary layer over a rough surface

Measurements of the vertical entrainment velocity into two-dimensional dense gas plumes over fully rough surfaces were carried out as part of a co-operative research programme with wind tunnel facilities in the USA. This paper presents results obtained for stable boundary layer conditions in the EnFlo wind tunnel at the University of Surrey; a companion paper treats the neutral boundary layer case. Mean velocity and temperature, turbulent normal and shear tresses, temperature fluctuations and heat fluxes were measured and used to demonstrate that a moderately stable atmospheric boundary layer had been successfully simulated in the tunnel. Entrainment velocities, W_E, were then deduced from the streamwise development of the concentration field, non-dimensionalised with respect to the friction velocity in the undisturbed flow, u^*, and correlated with the plume Richardson number, Ri^*. Higher non-dimensional entrainment speeds, W_E/u^*, were observed for Ri^*>5 in the stable boundary layer than in the neutral boundary layer, the difference growing with increasing Richardson number. Emission velocity ratios, W₀/u^*, were however larger in the stable experiments, and exceeded one at about Ri^*=18. Entrainment in the stable boundary layer appeared therefore to be more sensitive to emission velocity ratio than in the neutral case. Entrainment behaviour for Ri^*⩽5 followed that found in the neutral boundary layer. In this regime, use of the neutral boundary layer entrainment speed correlation is unlikely to lead to the over-prediction of plume dilution rates in moderately stable boundary layers.     

Statistical and diagnostic evaluation of a new-generation urban dispersion modelling system against an extensive dataset in the Helsinki area

We have developed a modelling system for predicting the traffic volumes, emissions from stationary and vehicular sources, and atmospheric dispersion of pollution in an urban area. This paper describes a comparison of the NO_x and NO₂ concentrations predicted using this modelling system with the results of an urban air quality monitoring network. We performed a statistical analysis to determine the agreement between predicted and measured hourly time series of concentrations at four permanently located and three mobile monitoring stations in the Helsinki Metropolitan Area in 1996–1997 (at a total of ten urban and suburban measurement locations). At the stations considered, the so-called index of agreement values of the predicted and measured time series of the NO₂ concentrations vary between 0.65 and 0.82, while the fractional bias values range from −0.29 to +0.26. In comparison with corresponding results presented in the literature, the agreement between the measured and predicted datasets is good, as indicated by these statistical parameters. The seasonal variations of the NO₂ concentrations were analysed in terms of the relevant meteorological parameters. We also analysed the difference between model predictions and measured data diagnostically, in terms of meteorological parameters, including wind speed and direction (the latter separately for two wind speed classes), atmospheric stability and ambient temperature, at two monitoring stations in central Helsinki. The modelling system tends to overpredict the measured NO₂ concentrations both at the highest (u⩾6 m s⁻¹) and at the lowest wind speeds (u<2 m s⁻¹). For higher wind speeds, the modelling system overpredicts the measured NO₂ concentrations in certain wind direction intervals; specific ranges were found for both monitoring stations considered. The modelling system tends to underpredict the measured concentrations in convective atmospheric conditions, and overpredict in stable conditions. The possible physico-chemical reasons for these differences are discussed.     

Minimum turbulence assumptions and u* and L estimation for dispersion models during low-wind stable conditions

The U.S. Environmental Protection Agency (EPA) short-distance dispersion model, AERMOD, has been shown to overpredict by a factor of as much as 10 when compared with observed concentrations from continuous releases at the Oak Ridge, TN (OR), and Idaho Falls, ID (IF), field experiments during stable periods when wind speeds often dropped below 1 m/sec. Some of this overprediction tendency can be reduced by revising AERMOD's meteorological preprocessor's parameterizations of the friction velocity, u _* , during low-wind stable conditions, thus increasing the calculated σ _v and σ _w and hence the lateral and vertical dispersion rates. Observations show that as the mean wind speed approaches zero at night, there is always significant σ _v and σ _w over time periods of 15 to 60 min, while standard Monin–Obukhov Similarity Theory (MOST) predicts that σ _v and σ _w will approach zero. This paper focuses on the u _* estimation methods and the minimum turbulence (σ _v and σ _w ) assumptions in AERMOD (beta option 4) and two widely used U.S. operational dispersion models, AERMOD (v12345) and SCICHEM. The U.S. EPA has provided results of its tests with the OR and IF data, with its base AERMOD version and its December 2012 modified versions, which assume adjustments to the low-wind u _* and increases in the minimum σ _v parameterization. SCICHEM has relatively small mean bias for both data sets. The revised AERMOD shows much less mean bias, agreeing more with SCICHEM.

Implications:

Suggestions are made for improvements to dispersion models such as AERMOD to correct overpredictions during light-wind stable conditions. Methods for estimating u_*, L, and the minimum turbulence parameters (σ_v and σ_w) are reviewed and compared. SCICHEM and the current operational version and an optional beta version (December 2012) of AERMOD are evaluated with tracer data from low-wind stable field experiments in Idaho Falls and Oak Ridge. It is seen that the operational version of AERMOD overpredicts by a factor of 2 to 10, while the optional beta version of AERMOD and SCICHEM have much less bias.     

Source Apportionment of PM10 and PM25 in Five Chilean Cities Using Factor Analysis

ABSTRACT

Chile is a fast-growing country with important industrial activities near urban areas. In this study, the mass and elemental concentrations of PM₁₀ and PM_2.5 were measured in five major Chilean urban areas. Samples of particles with diameter less than 10 um (PM₁₀) and 2.5 um (PM_2.5) were collected in 1998 in Iquique (northern Chile), Valparaiso, Vina del Mar, Rancagua (central Chile), and Temuco (southern Chile). Both PM₁₀ and PM₂₅ annual mean concentrations (PM₁₀: 56.9-77.6 u,g/m³; PM₂₅: 22.4-42.6 u.g/m³) were significantly higher than the corresponding European Union (EU) and U.S. Environmental Protection Agency (EPA) air quality standards. Moreover, the 24-hr PM₁₀ and PM₂₅ U.S. standards were exceeded infrequently for some of the cities (Rancagua and Valparaiso).     

Intercomparison of nocturnal mixing height estimate methods for urban air pollution modelling

One of the most important meteorological input parameters for three-dimensional photochemical air pollution models is the mixing height h, which has a strong influence on the shape and intensity of the vertical diffusivity K_z and, as a consequence, on ground-level air concentrations of primary and secondary pollutants. A number of indirect algorithms for the estimate of h in nocturnal, stable conditions, when the mixing is dominated by mechanical turbulence, are reviewed and compared with mixing heights derived from wind (SODAR) and temperature (RASS) profiles measured in the Milan urban area during spring and summer 1996. Mixing heights derived from temperature soundings correlate positively with those derived from wind soundings only when a stable layer is superimposed to a quasi-adiabatic layer, while the correlation is very weak in the presence of a ground-based inversion. In general, indirect algorithms perform very poorly if compared with RASS-based estimates, and reasonably well if compared with SODAR-based estimates. Among the others, Benkley and Schulman (1979, Journal of Applied Meteorology 18, 772–780) method, which makes use of wind speed observed at 10 m height, and Nieuwstadt (1984, Boundary-Layer Meteorology 30, 31–55), which makes use of friction velocity and Monin–Obukhov length, give the best results.     

The use of wind fields in a land use regression model to predict air pollution concentrations for health exposure studies

A methodology is developed to include wind flow effects in land use regression (LUR) models for predicting nitrogen dioxide (NO₂) concentrations for health exposure studies. NO₂ is widely used in health studies as an indicator of traffic-generated air pollution in urban areas. Incorporation of high-resolution interpolated observed wind direction from a network of 38 weather stations in a LUR model improved NO₂ concentration estimates in densely populated, high traffic and industrial/business areas in Toronto-Hamilton urban airshed (THUA) of Ontario, Canada. These small-area variations in air pollution concentrations that are probably more important for health exposure studies may not be detected by sparse continuous air pollution monitoring network or conventional interpolation methods. Observed wind fields were also compared with wind fields generated by Global Environmental Multiscale-High resolution Model Application Project (GEM-HiMAP) to explore the feasibility of using regional weather forecasting model simulated wind fields in LUR models when observed data are either sparse or not available. While GEM-HiMAP predicted wind fields well at large scales, it was unable to resolve wind flow patterns at smaller scales. These results suggest caution and careful evaluation of regional weather forecasting model simulated wind fields before incorporating into human exposure models for health studies. This study has demonstrated that wind fields may be integrated into the land use regression framework. Such integration has a discernable influence on both the overall model prediction and perhaps more importantly for health effects assessment on the relative spatial distribution of traffic pollution throughout the THUA. Methodology developed in this study may be applied in other large urban areas across the world.     

An inexpensive method to estimate CO and NOx emissions from mobile sources

CO and NO_x measurements from mobile sources at two urban locations in Córdoba City, Argentina, were used to develop a very simple method to estimate emission from these sources. This development was possible because primary urban air pollution in Córdoba comes mostly from mobile sources and because a field measurement campaign was conducted by the city government during 1995–1996 that has allowed us to have a complete and valuable data bank. Air concentrations of CO, NO_x as well as physical, and meteorological variables were measured at two urban sites with two monitoring stations. We compared the measured CO and NO_x air concentration data with the predictions of a method that uses regression analysis to estimate the emission factor from the mobile sources. The agreement is good, considering the simplicity of the approach.     

A measurement campaign in a street canyon in Helsinki and comparison of results with predictions of the OSPM model

In 1997, a measuring campaign was conducted in a street canyon (Runeberg St.) in Helsinki. Hourly mean concentrations of CO, NO_x, NO₂ and O₃ were measured at street and roof levels, the latter in order to determine the urban background concentrations. The relevant hourly meteorological parameters were measured at roof level; these included wind speed and direction, temperature and solar radiation. Hourly street level measurements and on-site electronic traffic counts were conducted throughout the whole of 1997; roof level measurements were conducted for approximately two months, from 3 March to 30 April in 1997. CO and NO_x emissions from traffic were computed using measured hourly traffic volumes and evaluated emission factors. The Operational Street Pollution Model (OSPM) was used to calculate the street concentrations and the results were compared with the measurements. The overall agreement between measured and predicted concentrations was good for CO and NO_x (fractional bias were −4.2 and +4.5%, respectively), but the model overpredicted the measured NO₂ concentrations (fractional bias was +22%). The agreement between the measured and predicted values was also analysed in terms of its dependence on wind speed and direction; the latter analysis was performed separately for two categories of wind velocity. The model qualitatively reproduces the observed behaviour very well. The database, which contains all measured and predicted data, is available for further testing of other street canyon dispersion models. The dataset contains a larger proportion of low wind speed cases, compared with other available street canyon measurement datasets.     

An experimental investigation of some methods of estimating turbulence parameters for use in dispersion models

The standard deviations of wind fluctuations in the horizontal and vertical directions, σ_θ and σ_φ, are now used in some pollution dispersion models to estimate the plume spread parameters σ_y and σ_z. Methods exist for estimating σ_θ and σ_φ when direct measurements are unavailable, using routine weather observations or wind measurements and temperature profiles from meteorological towers. In this paper such estimates are compared with direct measurements made at a height of 56 m, for a sampling time of 1 h, for a range of meteorological conditions. The work was carried out at a site in relatively irregular terrain. This was flat to rolling with a mixed surface cover within 1 km of the tower, with hills rising beyond that distance. Profile measurements were made with robust instruments rather than research grade sensors.Estimates of σ_φ made during the daytime agreed well with measurements, with a bias in the estimates of less than 0.4°. The r.m.s. differences between estimates and measurements were 1.1° (profile method) and less than 2° (routine weather observations method). Daytime σ_θ estimates were generally too low (bias 5–6°), although they were positively correlated with the measurements. At night σ_θ, was severely underestimated, and σ_φ was also underestimated.     

Maximum Ground-Level Concentrations with Downwash—Analysis

ABSTRACT

Equations derived previously for critical downwind distance x , wind speed u , and plume rise z , the values that produce maximum ground-level concentrations (MGLC) X_c under downwash conditions, have been solved. Tables of %_c, x_c, u_c, and z_c, and graphs of the relationships among u_c and z_c for a range of stack heights h_s, and building heights h_b, are presented. Results for two types of sources— a turbine and a reciprocating engine—are discussed. Some comparisons are made to the U.S. Environmental Protection Agency's (EPA) SCREEN3 model.     

Unsteady characteristics of the dispersion process in the vicinity of a pig barn. Wind tunnel experiments and comparison with field data

Two complementary methods, field experiments and physical modelling in a wind tunnel, have been used to investigate the dispersion of tracer-gas released from the ventilation system of a pig barn, under near-neutral stability conditions. In both cases, concentration fluctuations were measured and the deduced statistical results were compared. The choice of data processing applied to the time series of concentration was motivated by special issues in the assessment of odour annoyances: “where, how often, how long and how strong does it smell?” These features were described by the mean concentration distribution, the intermittency factor, the persistence and the 90-percentile. The good agreement between field and wind tunnel data confirmed the ability to replicate in wind tunnel the unsteady properties of a dispersion process, if the unsteady turbulent behaviour of the atmospheric boundary layer was properly modelled.A parametrical study of the influence on the dispersion process of the ratio between the exhaust velocity from the stack and the wind speed was then performed in wind tunnel. The fundamental outcome was that the near-field dispersion process under neutral stability conditions, despite the strong influence of the building wake, was for the most part driven by the meandering behaviour of the plume and not so much by the diffusion process.This study was also focused on the influence of the averaging time on the statistical results. The scatter generated by using dimensionless averaging times 200<T_a^*<400 (used during field experiments) instead of T_a^*→∞ (averaging time to ensure reproducible statistic results) was quantified in the wind tunnel. A degree of representativity of the results obtained from short-term samples, compared to fully converged statistical results was therefore assessed.     

Spatial distribution of source locations for particulate nitrate and sulfate in the upper-midwestern United States

Two back-trajectory analysis methods designed to be used with multiple site data, simplified quantitative transport bias analysis (SQTBA) and residence time weighted concentration (RTWC), were applied to nitrate and sulfate concentration data from two rural sites (the Mammoth Cave National Park and the Great Smoky Mountain National Park) and five urban sites (Chicago, Cleveland, Detroit, Indianapolis, and St. Louis) for an intensive investigation on the spatial patterns of origins for these two species in the upper-midwestern area. The study was made by dividing the data into five categories: all sites and all seasons, rural sites in summer, rural sites in winter, urban sites in summer, and urban sites in winter. A general conclusion was that the origins of the nitrate in these seven sites were mainly in the upper-midwestern areas, while the sulfate in these seven sites were mainly from the Ohio and Tennessee River Valley areas. The upper-midwestern areas are regions of high ammonia emissions rather than high NO_x emissions. In the winter, metropolitan areas showed the highest nitrate emission potential suggesting the importance of local NO_x emissions. In the summer, ammonia emissions from fertilizer application in the lower midwestern area made a significant contribution to nitrate in the rural sites of this study. The impact of the wind direction prevalence on the source spatial patterns was observed by comparing the urban and rural patterns of the summer. The differences between the results of two methods are discussed and suggestions for applying these methods are also provided.     

A study of surface ozone and the relation to complex wind flow in Hong Kong

Ozone measurements made from 5 sites in Hong Kong have been analyzed, including those from one upwind, one downwind, and three urban locales. The data are analyzed in terms of the seasonal and diurnal trends. A subset of data in autumn is further analyzed to study the relationship between the ozone spatial pattern and wind flow as well as other meteorological parameters. The results show that averaged ozone levels at most sites exhibit maxima in autumn, which appears to be a unique feature for eastern Asia. On average the daily maximum 1-h concentrations are found to be higher in the western (normally downwind) site than those on the eastern side and in urban areas. Examination of surface wind patterns and other meteorological parameters suggest that elevated ozone concentrations on the western side occur during the days with intense solar radiation, light winds, and in the presence of a unique wind circulation. The wind reversal in the western parts under the “convergence” flow is believed to be an important cause of the high-ozone events observed there. Such wind flow may re-circulate/transport nearby urban plumes (in this case the Hong Kong–Shenzhen urban complex). Examination of chemical data from the western site has shown that averaged afternoon SO₂ to NO_x ratios on days with wind reversal are larger than those of typical urban Hong Kong and that a significant SO₂ enhancement was clearly indicated on several occasions. The SO₂ enhancement may be interpreted as being the evidence to suggest the contribution of regional sources and/or Hong Kong’s power plants (both containing high SO₂). A case study has shown that when moderately strong northwesterly wind prevails, elevated ozone and SO₂ can be transported to western Hong Kong from the inner Pearl Delta region. This study has also indicated that under the impact of ENE winds the eastern side of Hong Kong is not frequently affected by the re-circulating ozone plumes present in the western side.     

Uranium uptake and depuration in the aquatic invertebrate Chironomus tentans

Evaluation of aqueous uranium (U) uptake and depuration in larvae of the midge Chironomus tentans were investigated in two separated experiments. First, a static-renewal experiment was performed with 10-d old C. tentans larvae exposed to 300 μg U/L. The animals steadily accumulated U (K_u = 20.3) approaching steady-state conditions (BAF = 56) in approximately 9-11 d. However, accumulated U was readily depurated (K_d = 0.36) with U tissue concentration decreasing rapidly within 3 d of the larvae being placed in clean water (t_1/2 = 1.9 d). Also, the growth of C. tentans larvae appeared to decrease after 6-11 d of U exposure, probably due to the reallocation of resources into U detoxification mechanisms. However, growth significantly increased once C. tentans were transferred to clean water. A separate short-term experiment was performed to evaluate the possible mechanism of U uptake in this invertebrate. Results suggested a passive mechanism of U uptake coupled with an active mechanism of U depuration but no details related to the type of mechanisms or pathway was investigated.     

Measurements and analyses of nitrogen oxides and ozone in the yard and on the roof of a street-canyon in Suzhou

The concentrations of air pollutants such as nitrogen oxides and ozone characterised by very fast chemical reactions can significantly vary within urban street-canyon due to the short distances between sources and receptor. With the primary objective to analyse this issue, NO, NO₂, NO_x, O₃, BTX, and wind flow field were continuously measured for 1 week at two heights (a street-level yard and a 25-m-high rooftop) in an urban canyon in Suzhou (China). The yard ozone concentrations were found to be up to six times lower than on the roof. Different frequency distributions (FD), dynamical and chemical processes of the pollutant variations from yard to roof are discussed to explain the findings. The predominant factors for the dissimilar pollutant vertical diffusion at the two measurement locations were associated to dissimilar fluid-dynamic and heterogeneous removal effects that likely induced dissimilar ozone chemical processes relative to NO_x and BTX precursors.     

Prediction of gas diffusion in complicated terrain by a potential flow model

A numerical model was developed to simulate gaseous diffusion in complicated terrain. This model calculates the air flow as a potential flow by the Boundary Element Method, and gaseous diffusion by an analytical Gaussian equation in the potential flow. Plume spreads σ_y and σ_z are modified by multiregression equations derived from wind tunnel experiments, and the terrain height is elongated depending on the atmospheric stability.First, tracer data from Cinder Cone Butte in the U.S. measured by the U.S.-EPA were predicted by the model in order to examine the prediction accuracy under stable conditions. The averaged ratio of the observed concentration to predicted concentration for 12 runs was better than a factor of 10. Next, tracer data from the Geysers area in the U.S. measured by the U.S.-DOE were used to examine the prediction accuracy under neutral conditions. The ratio of the observed concentration to predicted concentration for two runs under neutral conditions was better than factor of two at most locations, but prediction capability is poor in blocked or separated flow conditions.     

An assessment of the impact of climate change on air quality at two UK sites

Possible effects of climate change on air quality are studied for two urban sites in the UK, London and Glasgow. Hourly meteorological data were obtained from climate simulations for two periods representing the current climate and a plausible late 21st century climate. Of the meteorological quantities relevant to air quality, significant changes were found in temperature, specific humidity, wind speed, wind direction, cloud cover, solar radiation, surface sensible heat flux and precipitation. Using these data, dispersion estimates were made for a variety of single sources and some significant changes in environmental impact were found in the future climate. In addition, estimates for future background concentrations of NO_x, NO₂, ozone and PM₁₀ upwind of London and Glasgow were made using the meteorological data in a statistical model. These showed falls in NO_x and increases in ozone for London, while a fall in NO₂ was the largest percentage change for Glasgow. Other changes were small. With these background estimates, annual-average concentrations of NO_x, NO₂, ozone and PM₁₀ were estimated within the two urban areas. For London, results averaged over a number of sites showed a fall in NO_x and a rise in ozone, but only small changes in NO₂ and PM₁₀. For Glasgow, the changes in all four chemical species were small. Large-scale background ozone values from a global chemical transport model are also presented. These show a decrease in background ozone due to climate change. To assess the net impact of both large scale and local processes will require models which treat all relevant scales.     

Assessment of the uncertainty of using an inverse-dispersion technique to measure methane emissions from animals in a barn and in a small pen

Dispersion techniques are useful for assessing the effectiveness of management practices in mitigating methane (CH₄) emissions from animal operations. In this paper, the accuracy of an inverse-dispersion technique was examined for two common situations: 1) emissions from animals in a barn and 2) emissions from animals in a pen. For both situations, the accuracy of emission estimates generally increased with increasing fetch (scaled to barn height, h, or to pen width, X_s) between the source and the concentration measurement. The gas recovery ratio increased from 0.66 at 5 h to 0.93–1.03 at 10 h–25 h for the barn, and decreased from 1.59 at 0 X_s to 0.99 at 5 X_s and from 1.54 at 0 X_s to 1.09 at 5 X_s for the two small pen scenarios, due to the declining sensitivity of the dispersion model to errors in wind complexity and errors in the assumed source configuration. However, the relative uncertainty of the measured concentration increased with fetch due to the decreasing gas concentrations. Hence, improving the accuracy of emission estimates is a compromise between the source configuration sensitivity and the concentration measurement sensitivity. Fetches of about twice the size of the pen and about 10 times the height of the barn are suggested.     

Dilution of exhaust from a rooftop stack on a cubical building in an urban environment

The dispersion of a tracer gas emitted from a short stack on a cubical building was investigated using field and wind tunnel experiments. Air samples were obtained on the roof and the leeward side of the building. Dilution data were compared with estimates of minimum dilution (D_min) obtained with design formulas of Wilson/Chui/Lamb and Halitsky. The Halitsky model produced conservative predictions of D_min. Estimates of D_min obtained with the Wilson/Chui/Lamb model were in reasonable agreement with the field data. Wind tunnel dilution values were usually within a factor of two of the field data. The accuracy of the wind tunnel data improved as distance from the source increased. Wind tunnel dilution was found to be strongly dependent on the ratio of exhaust speed to wind speed, M, particularly for M values in the range of 2–4. This dependence is believed to be associated with the wind tunnel modelling of the stack exhaust, and is thus probably not a feature of the full-scale situation.