Urban air pollution modelling and measurements of boundary layer height

An urban field trial has been undertaken with the aim of assessing the performance of the boundary layer height (BLH) determination of two models: the Met Office Unified Model (UM) and a Gaussian-type plume model, ADMS. Pulsed Doppler lidar data were used to measure mixing layer height and cloud base heights for a variety of meteorological conditions over a 3 week period in July 2003. In this work, the daily growth and decay of the BLH from the lidar data and model simulations for 5 days are compared. The results show that although the UM can do a good job of reproducing the boundary layer growth, there are occasions where the BLH is overestimated by 30–100%. Within dispersion models it is the BLH that effectively limits the height to which pollution disperses, so these results have very important implications for pollution dispersion modelling. The results show that correct development of the boundary layer in the UM is critically dependant on morning cloud cover. The ADMS model is used routinely by local authorities in the UK for local air-quality forecasting. The ADMS model was run under three settings; an ‘urban’ roughness, a ‘rural’ roughness and a ‘transition’ roughness. In all cases, the ‘urban’ setting over estimated the BLH and is clearly a poor predictor of urban BLH. The ‘transition’ setting, which distinguishes between the meteorological data input site and the dispersion modelling site, gave the best results under the well mixed conditions of the trial.     

A study of the variation of urban mixed layer heights

The AERMET model is used to estimate hourly mixing heights during the Joint URBAN (2003) experiment in Oklahoma City, Oklahoma. AERMET is a simple 2-D model that requires only routine meteorological observations and an early morning atmospheric sounding to estimate convective boundary layer (CBL) growth. Estimated mixing heights are compared with observed mixing heights measured during Joint URBAN 2003. Observed CBL heights are derived from profiler data using a peak signal-to-noise ratio method. The method of deriving mixing heights from profiler data is validated using daily atmospheric sounding data. Estimated mixing heights using AERMET show good agreement with observations on days of varying temperature and cloud cover. AERMET was able to estimate the rapid boundary layer growth observed in the late morning and early afternoon hours during highly convective conditions. CBL heights of over 3000 m are observed in sounding data during the late afternoon. Estimated CBL heights of over 3000 m during the late afternoon agreed well with observations from the sounding and profiler data.     

Some Factors that Affect the Deposition Rates of Sulfur Dioxide and Similar Gases on Vegetation

The deposition of sulfur dioxide on growing vegetation is affected by diverse environmental factors, many of which undergo large diurnal and spatial variations. The aerodynamic resistance to vertical transfer in the surface boundary layer can be formulated in terms of the friction velocity, height of observation, vertical heat flux, and surface roughness. Also important are the resistance in the air layer closest to the surface elements and, in dry vegetation, the average stomatal resistance of the plants. The latter variable is among the most difficult to estimate, but over many agricultural field crops like those in the midwestern U.S., a typical minimum value of average stomatal resistance to SO₂ transfer is about 0.7 s cm^-1, as is indicated by various experimental data. The deposition velocity can be estimated as the inverse of the sum of the resistances of the layers, necessarily down to where the concentrations are zero; in the surface boundary layer, any of the various resistances might be dominant. Above the surface layer, the micrometeorological relationships are known with less certainty, but reasonable approximations indicate that during unstable conditions the resistance to transfer is very small at heights of several tens of meters and during stable conditions the aerodynamic resistance is very large aloft.     

Long-term changes in the tropospheric column ozone from the ozone soundings over Europe

The tropospheric column of ozone is analyzed from the measurements of the vertical profile of ozone by balloon-born ozonesondes. The data base includes ∼16,000 ozone profiles collected above six European stations—three of them have begun the ozonesoundings since 1970. We present a trend analysis (with data up to 2005) focusing on detection of the long-term tropospheric ozone variability over the European network. The ozone time series have been examined separately for each station and season (DJF, MAM, JJA, SON) using a flexible trend model. A trend component of the model is taken as a smooth curve without a priori defined shape. A large increase in the European tropospheric ozone since the beginning of the 1970s (net change of ∼10% in summers and ∼30% in winters) and a kind of stabilization in the early 1990s have been corroborated by the study. This pattern comes from the most extensive data set of ozonesoundings over Hohenpeissenberg and Payern. The decadal differences in the trend pattern between these and other European stations are disclosed. The results of a stepwise regression model using various characteristics of the ozone and temperature profiles as explanatory variables for the tropospheric column ozone (TCO₃) variations show that the ozone changes may be reconstructed using the ozone mixing ratio at 500 hPa, the thermal tropopause (TT) height, and the difference between ozonepause and TT heights. It appears that the last two factors induce 20–30% of the net TCO₃ change over Hohenpeissenberg in the 1970–2004 period.     

Principal Plume Dispersion Models: TVA Power Plants

The body of information presented in this paper is directed to those individuals who may be concerned with principal plume dispersion models at coal-burning power plants. About 20 years of comprehensive field surveillance and documentation of dispersion of power plant emissions for a varied range of unit sizes, stack heights, and meteorological conditions have determined the Tennessee Valley Authority’s interpretation of principal plume dispersion models. TVA’s experience indicates that as unit sizes are increased and taller stacks are constructed, the plume dispersion model associated with maximum surface concentrations changes. Maximum surface concentrations for principal plume dispersion models were approximately equal for the early small plants. However, the coning model was considered the critical plume dispersion model because the frequency of recurrence of surface concentrations from this model was appreciably greater than other models.

There were progressive changes because of an increase in unit sizes and stack heights; the magnitude of maximum surface concentrations from the coning model decreased, and the magnitude (relative to the coning model) of concentrations from the inversion breakup model increased. However, with plumes from newer and larger units with higher stacks, the trapping dispersion model became prominent. Finally, by the time unit size had increased to 900 mw and stack height to about 245 meters, as at Bull Run Power Plant, the magnitude of surface concentrations associated with trapping had increased to such a degree that it became the critical dispersion model identified with power plants of this size.     

Recording the Response of Plants to Various Air Pollutants

A discussion of the methods used to determine the most economic design of chimney for a new thermal power station or large industrial plant is presented, with the objective that ground level concentration of pollutants will be kept at a minimum. Attention is paid to the geography and climatology of the site, with special reference to the frequency and height of inversions and the prevailing wind direction and speed.

A method is illustrated in using a large thermal power station as an example. The maximum sulfur dioxide concentrations at ground level are computed for several chimney heights and gas exit velocities. The values of these sulfur dioxide concentrations, the capital cost of the chimney, the pumping costs, and the gas pressures within the chimney are considered in selecting a suitable chimney height and a gas exit velocity which will meet most economically the stated objective.

The paper deals primarily with chimneys for industrial or power boiler plant of maximum continuous rating greater than 450 million Btu/hr (about 450,000 lbs of steam/hr), or to chimneys serving furnaces burning fuel at a maximum rate greater than 50,000 lbs/hr of coal, or 80,000 lbs/hr of oil. For chimneys serving plant with smaller heat inputs, chimney selection by reference to Clean Air Act 1956, Memorandum on Chimney Heights is suggested.     

Validation of ADMS against wind tunnel data of dispersion from chemical warehouse fires

Comparisons are presented of the predictions of the atmospheric dispersion modelling system (ADMS) and wind tunnel data for plume dispersion from chemical warehouse fires. The focus of the comparisons is dispersion from structurally intact buildings with open roofs and dispersion of plumes flush with the ground without obstacles, however, dispersion from building shells and doors is also considered. Both buoyancy driven and momentum driven flows are treated, although emphasis is on buoyancy driven flows as these are generally more likely to occur in warehouse fires. The study shows that the ADMS building module is able to reproduce many of the features of dispersion observed in the wind tunnel. These include a recirculating region behind the building in which material may be trapped, a main wake which brings material down towards the surface, and appropriate sensitivity to the buoyancy and momentum of the emitted material, and the location of sources on the building roof. The comparisons suggest that the ADMS building model can be used to predict dispersion from the stages of fire development studied. The precise level of agreement depends (but not in a systematic way) on the buoyancy flux parameter F_B, the momentum flux parameter F_M and the number of roof lights. There are some significant differences between the wind tunnel boundary layer and the simulated atmospheric boundary layer in ADMS which have to be considered when making wind tunnel model comparisons. These relate mainly to the near surface where the wind tunnel underestimates turbulent velocities, the boundary layer height which in the wind tunnel corresponds to an atmospheric boundary layer depth of 82.5 m (atmospheric boundary layers are frequently an order of magnitude deeper), and the boundary layer top where the ADMS boundary layer is capped by an inversion and has low turbulence levels whereas the wind tunnel boundary layer has higher levels of turbulence and no capping inversion.     

Meteorological conditions at the release site during the two tracer experiments

The state of the boundary layer at the release site during the two tracer experiments is described, based on measurements of atmospheric turbulence carried out by a sonic anemometer, profiles of horizontal and vertical wind by a SODAR as well as frequent radiosonde releases. The boundary layer height is derived from radiosoundings, modelled and discussed. The study is meant as background for discussions of the tracer behaviour near the release site, and possible influence of the local meteorology at the release site on the long-range dispersion. The difference in the initial meso-meteorological conditions between the two experiments is discussed and some of its consequences for the dispersion of the released tracer are clarified. It is concluded that the first experiment is well suited for validation of long-range transport and dispersion models, whereas the second experiment calls for further model development. It is concluded that the micrometeorological measurements are very important as background information for the understanding of plume behaviour on small as well as long scales.     

Estimating micrometeorological inputs for modeling dispersion in urban areas during stable conditions

We examine the performance of three methods to estimate the surface friction velocity and the Monin–Obukhov (MO) length in stable conditions. Estimates from these methods are compared with measurements made at two urban sites: the Wilmington site located in the middle of an urban area, and the VTMX site located on a sloping, smooth area in Salt Lake City. The first method uses the mean wind at a single height (Single U or SU), the second uses the wind speed at a single level and the temperature difference between two levels (U delta T or UDT), and the third method uses two levels of wind speed and temperature (delta U delta T or DUDT). The performance of the SU and UDT methods in estimating u_* are comparable. The SU method yields better estimates of the MO length than the UDT method does. The DUDT method performs poorly in estimating both u_* and L. The major conclusions of this study are that (1) measurements of mean winds and temperatures at one or two levels at an urban location can provide adequate estimates of micrometeorological variables required in modeling dispersion in the stable boundary layer, and (2) methods based on using differences in temperatures and velocities between two levels can provide unreliable estimates of these variables because these differences can be overwhelmed by inevitable uncertainties in the measurement of mean variables.     

A hybrid plume model for local-scale atmospheric dispersion

This paper presents a local-scale dispersion model, based on atmospheric boundary layer scaling theory. In the vicinity of the source, Gaussian equations are used in both the horizontal and vertical directions. After a specified transition distance, gradient transfer theory is applied in the vertical direction, while the horizontal dispersion is still assumed to be Gaussian. The dispersion parameters and eddy diffusivity are modelled in a form, which facilitates the use of a meteorological pre-processor. We present a novel model of the vertical eddy diffusivity (K_z), which is a continuous function of height in various atmospheric scaling regions. The model also includes a treatment of the dry deposition of gases and particulate matter. The accuracy of the numerical model was analysed by comparing the model predictions with two analytical solutions; the numerical deviations from these solutions were less than 2% for the computational regime. The model has been tested against the Kincaid experimental field data. The agreement of the predictions and the data is good on the average, although the internal variation of the predictions versus data scatter plot is substantial.     

Lidar-assisted measurement of PM10 emissions from agricultural tilling in California's San Joaquin Valley – Part I: lidar

Vertical profiling with point samplers is an accepted method for quantifying the fluxes of PM₁₀ from non-point fugitive dust sources, but is limited by uncertainty in estimates of the actual height of the dust plume, especially for plumes that exceed the highest sampling height. Agricultural land preparation operations in the San Joaquin Valley were monitored using upwind–downwind vertical PM₁₀ profiles and data collected during the first successful experiment to include light detection and ranging (lidar), in 1998, were analyzed to provide modeling criteria for the 1996 and 1997 data. A series of six comprehensive PM₁₀ tests with concurrent lidar data was examined to: (a) develop a framework for analyzing upwind–downwind point PM₁₀ concentration profiles of land preparation operations (disking, listing, root cutting, and ripping) and (b) identify conditions under which the field sampling strategies affect the reproducibility of PM₁₀ concentration measurements. Lidar data were used to verify that the plume heights and shapes extrapolated from the point sampler vertical profiles adequately described the plumes. The shortcomings of the vertical profiling technique and lidar methods are discussed in the light of developing efficient robust methods for accurate PM₁₀ emissions quantification from complex non-point sources.     

Impact of using prognostic and objective wind fields on the photochemical modeling of Athens, Greece

A three-dimensional Eulerian photochemical model is used to follow the dynamics of ozone, NO_x, and CO over the Athens area, for 25 May 1990, the day considered in the APSIS project. A unique aspect of this work lies in the study of the impacts of the wind field preparation methods on the concentrations predicted by the model. Three sets of wind fields are developed. The first one used is derived from a prognostic meteorological model. The second one is calculated from available wind observations using objective: methods. For these two cases, a previous day is simulated, using the same conditions, to develop preconditioned initial conditions for the following day. For the third simulation, again two days are simulated, this time using the observed winds for each of the two days modeled. The predictions using the prognostically derived and the objective analysis wind fields are significantly different, particularly for the primary pollutants. Comparing predictions to the observations did not favor any particular method of wind field preparation. In this case, when using the prognostically derived field, the simulations are very sensitive to boundary conditions. In contrast, when using the wind fields constructed by objective methods, the simulations became most sensitive to emissions and initial conditions. This comes directly from the different residence times in the domain, which are governed by the wind speed.     

The dispersion of the Buncefield oil fire plume: An extreme accident without air quality consequences

The dispersion of pollutants from the huge Buncefield oil depot fire that occurred on 11 December 2005 is simulated using a regional Eulerian chemistry-transport model. We analyse the transport and mixing of the fire plume. We show that the hot plume never reached the ground. Instead, it pierced the thin wintertime boundary layer and was injected into the free troposphere at higher altitudes. This is in agreement with data from many air quality stations. This high injection was fortunate because the fine aerosol particles (PM₁₀) mass column generated by the fire smoke exceeded that of ordinary pollution by an order of magnitude. Our regional chemistry-transport modelling is able to predict the early development of the plume dispersion, as shown by a qualitative comparison between simulated PM₁₀ columns and a satellite image obtained by the EOS-TERRA-MODIS sensor.If the accident had occurred in summer when boundary layers are much deeper and convective, a severe degradation in air quality due to PM₁₀ could have occurred, as shown by a sensitivity simulation assuming a similar fire during one of the hottest days of August 2003. The modelled impact of the fire on regional and European air quality levels strongly depends on the altitude reached by the buoyant plume, as shown by a set of sensitivity simulations with variable injection heights. However, in all cases we found that the fire only affected surface aerosol concentrations without increasing photochemical pollution.     

The effect of humidity and state of water surfaces on deposition of aerosol particles onto a water surface

Deposition processes of particles with dry diameter larger than about 10 μm are dominated by gravitational settling, while molecular diffusion and Brownian motion predominate the deposition processes of particles smaller than 0.1 μm in dry diameter. Many air pollution derived elements exhibit characteristics common to sub-micron particles. The objective of the present study is to examine the effects of meteorological conditions within the turbulent transfer layer on the deposition velocity of particles with dry diameter between 0.1 and 1 μm. It is for these sub-micron particles that particle growth by condensation in the deposition layer, the broken water surface effect and the enhanced transfer process due to atmospheric turbulence in the turbulent transfer layer play important roles in controlling the particle deposition velocity. Results of the present study show that the `dry air’ assumption of Williams’ model is unrealistic. Effects of ambient air relative humidity and water surface temperature cannot be ignored in determining the deposition velocity over a water surface. Neglecting effects of ambient air relative humidity and water surface temperature will result in defining atmospheric stability incorrectly. It is found that the largest effect of air relative humidity on deposition velocity occurs at an air–water temperature difference corresponding to the point of `displaced neutral stability'. For a given wind speed of U=5 m s⁻¹ the additive effects of water surface temperature, T_w, changes from 5 to 25°C and ambient air relative humidity variations from 85 to 60%, respectively, lead to a maximum difference in v_d of about 20%. For a higher wind speed of 10 m s⁻¹, however, the corresponding change in v_d reduces to less than 5%. This is further confirmation that wind speed is one of the strongest variables that governs the magnitude of v_d. The present study also found that the broken surface transfer coefficient, k_bs, given as a multiple of the smooth surface transfer coefficient, k_ss, is physically more meaningful than assigning it a constant value independent of particle size. The method used in this study requires only a single level of atmospheric data coupled with the surface temperature measurement. The present method is applicable for determining deposition velocity not only at the conventional measurement height of 10 m but also at any other heights that are different from the measurement height.     

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.     

Operational ozone forecasts for the region of Copenhagen by the Danish Meteorological Institute

The Danish Meteorological Institute (DMI) has developed an operational forecasting system for ozone concentrations in the Atmospheric Boundary Layer; this system is called the Danish Atmospheric Chemistry FOrecasting System (DACFOS). At specific sites where real-time ozone concentration measurements are available, a statistical after-treatment of DACFOS’ results adjusts the next 48 h ozone forecasts. This post-processing of DACFOS’ forecasts is based on an adaptive linear regression model using an optimal state estimator algorithm. The regression analysis uses different linear combinations of meteorological parameters (such as temperature, wind speed, surface heat flux and atmospheric boundary layer height) supplied by the Numerical Weather Prediction model DMI-HIRLAM. Several regressions have been tested for six monitoring stations in Denmark and in England, and four of the linear combinations have been selected to be employed in an automatic forecasting system. A statistical study comparing observations and forecasts shows that this system yields higher correlation coefficients as well as smaller biases and RMSE values than DACFOS; the present post-processing thus improves DACFOS’ forecasts. This system has been operational since June 1998 at the DMI's monitoring station in the north of Copenhagen, for which a new ozone forecast is presented every 6 h on the DMI's internet public homepage.     

Boundary-layer dynamics and its influence on atmospheric chemistry at Summit,Greenland

Sonic anemometer turbulence measurements were made at Summit, Greenland during summer 2004 and spring 2005. These measurements allow for the characterization of the variability of the atmospheric boundary layer at this site by describing seasonal and diurnal changes in sensible heat flux and boundary layer stability as well as providing estimates of mixing layer height. Diurnal sensible heat fluxes at Summit ranged from −18 to −2 W m⁻² in the spring and from −7 to +10 W m⁻² in the summer. Sustained stable surface layer conditions and low wind speeds occured during the spring but not during the summer months. Unstable conditions were not observed at Summit until late April. Diurnal cycles of convective conditions during the daytime (0700–1700 h local time) were observed throughout July and August. Boundary layer heights, which were estimated for neutral to stable conditions, averaged 156 m for the spring 2005 observations. Comparisons of the boundary layer characteristics of Summit with those from South Pole, Antarctica, provide possible explanations for the significant differences in snowpack and surface-layer chemistry between the two sites.     

Lagrangian sampling of 3-D air quality model results for regional transport contributions to sulfate aerosol concentrations at Baltimore,MD, in summer 2004

We use ensemble-mean Lagrangian sampling of a 3-D Eulerian air quality model, CMAQ, together with ground-based ambient monitors data from several air monitoring networks and satellite (MODIS) observations to provide source apportionment and regional transport vs. local contributions to sulfate aerosol and PM_2.5 concentrations at Baltimore, MD, for summer 2004. The Lagrangian method provides estimates of the chemical and physical evolution of air arriving in the daytime boundary layer at Baltimore. Study results indicate a dominant role for regional transport contributions on those days when sulfate air pollution is highest in Baltimore, with a principal transport pathway from the Ohio River Valley (ORV) through southern Pennsylvania and Maryland, consistent with earlier studies. Thus, reductions in sulfur emissions from the ORV under the EPA's Clean Air Interstate Rule may be expected to improve particulate air quality in Baltimore during summer. The Lagrangian sampling of CMAQ offers an inexpensive and complimentary approach to traditional methods of source apportionment based on multivariate observational data analysis, and air quality model emissions separation. This study serves as a prototype for the method applied to Baltimore. EPA is establishing a system to allow air quality planners to readily produce and access equivalent results for locations of their choice.     

The impact of trajectory starting heights on the MURA trajectory source apportionment (TSA) method

Trajectory source apportionment (TSA) methods have been used in many research projects to attempt to identify the sources of pollution. Hybrid Single Particle Lagrangian Integrated Trajectories (HYSPLIT) is a popular model for use in various TSA methods. One of the options in this model is to choose a starting height. Very little research is available to assist a user in making this choice. This paper evaluates starting heights of 10, 50, 100, 250, and 500 m on the accuracy of the Multi-Receptor (MURA) method using artificial sources for three different simulations. It was found that using ensembles of trajectories in the MURA method appear to average out most of the biases found from different trajectory starting heights up to the 500 m tested.