Open access

Urban Air Pollution Modeling

Written By

Anjali Srivastava and B. Padma S. Rao

Submitted: 15 October 2010 Published: 05 July 2011

DOI: 10.5772/16776

From the Edited Volume

Air Quality - Models and Applications

Edited by Dragana Popovi?

Chapter metrics overview

7,072 Chapter Downloads

View Full Metrics

1. Introduction

All life form on this planet depends on clean air. Air quality not only affects human health but also components of environment such as water, soil, and forests, which are the vital resources for human development.

Urbanization is a process of relative growth in a country’s urban population accompanied by an even faster increase in the economic, political, and cultural importance of cities relative to rural areas. Urbanization is the integral part of economic development. It brings in its wake number of challenges like increase in population of urban settlement, high population density, increase in industrial activities (medium and small scale within the urban limits and large scale in the vicinity), high rise buildings and increased vehicular movement. All these activities contribute to air pollution. The shape of a city and the land use distribution determine the location of emission sources and the pattern of urban traffic, affecting urban air quality (World Bank Reports 2002). The dispersion and distribution of air pollutants and thus the major factor affecting urban air quality are geographical setting, climatological and meteorological factors, city planning and design and human activities.

Cities in the developing countries are characterized by old city and new development. The old cities have higher population density, narrow lanes and fortified structures.

In order to ensure clean air in urban settlements urban planning and urban air quality management play an important role. New legislations, public awareness, growth of urban areas, increases in power consumption and traffic pose continuous challenges to urban air quality management. UNEP (2005) has identified niche areas Urban planning need to primarily focus on as:

  • Promotion of efficient provision of urban infrastructure and allocation of land use, thereby contributing to economic growth,

  1. managing spatial extension while minimizing infrastructure costs,

  2. improving and maintaining the quality of the urban environment and

  3. Prerecording the natural environment immediately outside the urban area.

Air quality modelling provides a useful support to decision making processes incorporating environmental policies and management process. They generate information that can be used in the decision making process. The main objectives of models are: to integrate observations, to predict the response of the system to the future changes, to make provision for future development without compromising with quality


2. Urban air quality

The urban air is a complex mixture of toxic gases and particulates, the major source is combustion of fossil fuels.Emissions from fossil fuel combustion are reactive and govern local atmospheric chemistry.Urban air pollution thus in turn affect global troposphere chemistry and climate (e.g. tropospheric O3 and NOX budgets, radiative forcing by O3 and aerosols).

Sources of air pollutants in urban area, their effect and area of concern are summarized in Table 1.

SourcePollutantsEffectsArea of concern
Large number of vehiclesParticulate matters (PM10, PM2.5), Lead (Pb), Sulphur dioxide (SO2), Oxides of nitrogen (NOx), Ozone (O3), Hydro carbons (HCs), Carbon monoxide (CO), Hydrogen fluoride (HF), Heavy metals (e.g. Pb, Hg, Cd etc.)Human Health (acute and chronic)

Local, Regional and Global
Use of diesel powered vehicle in large number
Use of obsolete vehicles in large number
Large number of motor cycles/three wheelers (2 stroke and three stroke)Ecosystem
(acute and chronic)
Local, Regional and Global
Unpaved and/or poorly maintained street
Open burning
Inadequate infrastructureGreenhouse gas emissionGlobal
Low quality of fuel/fuel adulteration
Little emission control & technology in industry
Presence of industries (e.g. ceramic, brick works, agrochemical factory)Acid rainGlobal
Waste incinerationStratospheric ozone depletionGlobal
Limited dry deposition of pollutantsLong-range transportGlobal

Table 1.

Urban sources of air pollutants, their effect and area of concern.

Urban air pollution involves physical and chemical process ranging over a wide scale of time and space. The urban scale modeling systems should consider variations of local scale effects, for example, the influence of buildings and obstacles, downwash phenomena and plume rise, together with chemical transformation and deposition. Atmospheric boundary layer, over 10 to 30 km distances, governs the dispersion of pollutants from near ground level sources. Vehicular emissions are one the major pollution source in urban areas. Ultrafine particles are formed at the tailpipe due to mixing process between exhaust gas and the atmosphere. Processes at urban scale provide momentum sink, heat and pollutant source thereby influencing the larger regional scale (up to 200 km). Typical domain lengths for different scale models is given in table 2.

ModelTypical Domain ScaleTypical resolutionMotion Example
Micro scale200x200x100m5mMolecular diffusion, Molecular viscosity
Mesoscale (urban)100x100x5km2kmEddies, small plumes, car exhaust, cumulus clouds
Regional1000x1000x10km36kmGravity waves, thunderstorms, tornados, cloud clusters, local winds, urban air pollution
Synoptic(continental)3000x3000x20km80 kmHigh and low pressure system, weather fronts, tropical storms, Hurricanes Antarctic ozone hole,
Global65000x65000x20km4° x 5°Global wind speed, rossby (planetary) waves stratospheric ozone reduction Global worming

Table 2.

Typical domain length for different scale model

Piringer et al., 2007, have demonstrated that atmospheric flow and microclimate are influenced by urban features, and they enhance atmospheric turbulence, and modify turbulent transport, dispersion, and deposition of atmospheric pollutants. Any urban scale modeling systems should consider effects of the various local scales, for example, the influence of buildings and obstacles, downwash phenomena and plume rise, chemical transformation and deposition. The modelling systems also require information on emissions from various sources including urban mobile pollution sources. Simple dispersion air quality pollution transport models and complex numerical simulation model require wind, turbulence profiles, surface heat flux and mixing height as inputs. In urban areas mixing height is mainly influenced by the structure heights and construction materials, in terms of heat flux. Oke (1987, 1988, 1994), Tennekes (1973), Garrat (1978, 1980), Raupach et al (1980) and Rotach (1993, 1995) divided the Atmospheric Boundary Layer within the urban structures into four sub layers (Figure 1).

Figure 1.

Boundary- layer structure over a rough urban built- up area A daytime situation is displayed where Z I denotes the mixed layer height. Modefied after Oke ( 1988) and Rotach (1993).

In urban establishments anthropogenic activities take place between the top of highest building and the ground. People also live in this area. The layer of atmosphere in this volume is termed as Urban Canopy. The thermal exchanges and presence of structures in urban canopy modify the air flows significantly and this makes the atmospheric circulations in urban canopy highly complex. The heterogeneity of urban canopies poses a challenge for air quality modeling in urban areas. The importance of various parameters in different models for urban atmosphere study is given in Table 3. Figure 2 shows the flow and scale lengths within an urban boundary layer, UBL.

ParameterAir QualityUrban ClimatologyUrban Planning
Wind speedVery importantImportantVery Important
Wind DirectionVery importantImportantVery Important
ImportantExtremely importantVery Important
Pollutant ConcentrationExtremely importantVery importantVery important
Turbulent FluxesVery importantVery importantVery important

Table 3.

Ranking of parameters in different applications for urban air environment

Figure 2.

Schematic diagram showing processes, flow and scale lengths within an urban boundary layer, UBL. This is set in the context of the planetary boundary layer, PBL, the urban canopy layer, UCL, and the sky view factor, SVF, a measure of the degree to which the sky is obscured by surrounding buildings at a given point which characterises the geometry of the urban canopy. Ref:. Meteorology applied to urban pollution problems-Final report COST Action 715. Dementra Ltd Publishers

Vehicles are one of the important pollution sources in urban areas. Maximum exposure to local public is from this source and thus they form important receptor group. Pollutant dispersion of vehicular pollution is at street scale and is the smallest scale in urban environment. Hosker (1985) showed that flows in street canyon are like recirculating eddy driven by the wind flow at the top with a shear layer which separates the above canyon flows from those within it. In deep street canyons the primary vortex does not extend to the ground but a weak contra rotating vortex is formed near the ground and is relatively shallow (Figure 3). Pavageau et al (2001) demonstrated that wind directions which are not normal to the street axis cause variations in the flow. The real geometry of the street canyon and the mean flow and turbulence generated by vehicles within the canyon also affect the recirculating flow.

Concentrations of pollutants at a receptor are governed by advection, dispersion and deposition. Air pollutants can be divided into two main categories namely conventional air pollutants and Hazardous Air Pollutants (HAPs). Conventional air pollutants include particulate matters, sulphur dioxide, nitrogen dioxide, carbon monoxide, particles, lead and the secondary pollutant ozone. HAPs include Volatile Organic Compounds, toxic metals

Figure 3.

Air flow pattern in a Street Canyon

and biological agents of many types. All pollutants are not emitted in significant quantities. Secondary pollutants like some VOCs, carbonyls and ozone are formed due to chemical transformation in air. These reactions are often photochemical.

The important components of air quality modelling are thus,

  1. Knowledge of sources and emissions

  2. Transport, diffusion and parametrisation

  3. Chemical transformations

  4. Removal process

  5. Meteorology

Understanding contribution from various sources to air quality is the key for effective management of the air quality. Air quality models offer a useful tool in comprehending these issues. These models evaluate the relationship between air pollutant emissions and their resulting concentration in the ambient air. Commonly used air quality models are: 1) Conceptual Models 2) Emission Models 3) Meteorological Models, 4) Chemical Models, 5) Source Oriented Models and 6) Receptor Models.


3. Air quality model classification

Air quality models cover either separately or together atmospheric phenomena at various temporal and spatial scales. Urban air models generally focus from local (micro- tens of meters to tens of kilometers) to regional (meso) scale. Models can be broadly divided into two types namely physical and mathematical.

Physical models involve reproducing urban area in the wind tunnel. Scale reduction in the replica and producing scaling down actual flows of atmospheric motion result in limited utility of such models. Moreover these are economically undesirable.

Mathematical models use either use statistics to analyse the available data or mathematical representation of all the process of concern. The second type of mathematical models is constrained by the ability to represent physical and chemical processes in equations without assumptions.

Statistical model are simple but they do not explicitly describe causal relationships and they cannot be extrapolated beyond limits of data used in their derivation. Thus dependence on past data becomes their major weakness. These cannot be used for planning as they cannot predict effect of changes in emissions.

3.1. Eulerian and lagrangian models

Eulerian approach has been used to predict air pollutant concentrations in urban areas. The space domain (geographical area or air volume), are divided into "small" squares (two-dimensional) or volumes (three-dimensional), i.e. grid cells. Thus Eulerian models are sometimes called "grid models". Equidistant grids are normally used in air pollution modeling. Then the spatial derivatives involved in the system of Partial Differential Equations are discretized on the grid chosen. The transport, diffusion, transformation, and deposition of pollutant emissions in each cell are described by a set of mathematical expressions in a fixed coordinate system. Chemical transformations can also be included. Long range transport, air quality over entire air shed, that is, large scale simulations are mostly done using Eulerian models. Reynolds (1973), Shir and Shieh (1974) applied Eulerian model for ozone and for SO2 concentration simulation in urban areas, and Egan (1976) and Carmichael (1979) for regional scale sulfur. Holmes and Morawska (2006) used Eulerian model to calculate the transport and dispersion over long distances. The modeling studies by Reynolds (1973) on the Los Angeles basin formed the basis of the, the well-known Urban Air shed Model-UAM. Examples of Eulerian models are CALGRID model and ARIA Regional model or the Danish Eulerian Hemispheric Model (DEHM).

Lagrangian Model approach is based on calculation of wind trajectories and on the transportation of air parcels along these trajectories. In the source oriented models the trajectories are calculated forward in time from the release of a pollutant-containing air parcel by a source (forward trajectories from a fixed source) until it reaches a receptor site. And in receptor oriented models the trajectories are calculated backward in time from the arrival of an air parcel at a receptor of interest (backward trajectories from a fixed receptor). Numerical treatment of both backward and forward trajectories is the same. The choice of use of either method depends on specific case. As the air parcel moves it receives the emissions from ground sources, chemical transformations, dry and wet depositions take place. If the models provide average time-varying concentration estimates along the box trajectory then Lagrangian box models have been used for photochemical modeling. The major shortcoming of the approach is the assumption that wind speed and direction are constant throughout the Physical Boundary Layer. As compared to the Eulerian box models the Lagrangian box models can save computational cost as they perform computations of chemical and photochemical reactions on a smaller number of moving cells instead of at each fixed grid cell of Eulerian models. Versions of EMEP (European Monitoring and Evaluation Programme) are examples of Lagrangian models. These models assume pollutants to be evenly distributed within the boundary layer and simplified exchange within the troposphere is considered.

3.2. Box models

Box models are based on the conservation of mass. The receptor is considered as a box into which pollutants are emitted and undergo chemical and physical processes. Input to the model is simple meteorology. Emissions and the movement of pollutants in and out of the box is allowed. The air mass is considered as well mixed and concentrations to be uniform throughout. Advantage of the box model is simple meteorology input and detailed chemical reaction schemes, detailed aerosol dynamics treatment. However, following inputs of the initial conditions a box model simulates the formation of pollutants within the box without providing any information on the local concentrations of the pollutants. Box models are not suitable to model the particle concentrations within a local environment, as it does not provide any information on the local concentrations, where concentrations and particle dynamics are highly influenced by local changes to the wind field and emissions.

3.3. Receptor models

Receptor modeling approach is the apportionment of the contribution of each source, or group of sources, to the measured concentrations without considering the dispersion pattern of the pollutants. The starting point of Receptor models is the observed ambient concentrations at receptors and it aims to apportion the observed concentrations among various source types based on the known source profile (i.e. chemical fractions) of source emissions. Mathematically, the receptor model can be generally expressed in terms of the contribution from ‘n’ independent sources to ‘p’ chemical species in ‘m’ samples as follows:


Where Cik is the measured concentration of the kth species in the ith sample, aik is the concentration from the jth source contributing to the ith sample, and fjk is the kth species fraction from the jth source. Receptor models can be grouped into Chemical mass balance (CMB), Principal Component Analysis (PCA) or Factor analysis, and Multiple Linear Regression Analysis (MLR) and multivariate receptor models.

The Chemical Mass Balance (CMB) Receptor Model used by Friedlander, 1973 uses the chemical and physical characteristics of gases and particulate at source receptor to both identify the presence of and to quantify source contributions of pollutants measured at the receptor. Hopke (1973, 1985) christened this approach as receptor modelling. The CMB model obtains a least square solution to a set of linear equation, expressing each receptor concentration of a chemical species as a linear sum product of source profile species and source contributions. The output to the model consists of the amount contributed by each source type to each chemical species. The model calculates the contribution from each source and uncertainties of those values. CMB model applied to the VOC emissions in the city of Delhi and Mumbai (Figure 4 ) shows that emissions from petrol pumps and vehicles at traffic intersection dominate.

PCA and MLR are statistical models and both PMF and UNMIX are advanced multivariate receptor models that determine the number of sources and their chemical compositions and contributions without source profiles. The data in PMF are weighted by the inverse of the measurement errors for each observation. Factors in PMF are constrained to be nonnegative. PMF incorporates error estimates of the data to solve matrix factorization as a constrained, weighted least-squares problem (Miller et al., 2002; Paatero, 2004).

Geometrical approach is used in UNMIX to identify contributing sources. If the data consist of ‘m’ observations of ‘p’ species, then the data can be plotted in a p-dimensional data space, where the coordinates of a data point are the observed concentrations of the species during a sampling period. If n sources exist, the data space can be reduced to a (n-1) dimensional space. An assumption that for each source, some data points termed as edge points exist for which the contribution of the source is not present or small compared to the other sources.

Figure 4.

Category wise Contribution to Total VOCs at Mumbai and Delhi based on CMB results(Ref: Anjali Srivastava 2004, 2005)

UNMIX algorithm identifies these points and fits a hyperplane through them; this hyperplane is called an edge. If n sources exist, then the intersection of n-1 of these edges defines a point that has only one contributing source. Thus, this point gives the source composition. In this way, compositions of the n sources are determined which are used to calculate the source contributions (Henry, 2003).

3.4. Computational fluid dynamic models

Resolving the Navier-Stokes equation using finite difference and finite volume methods in three dimensions provides a solution to conservation of mass and momentum. Computational fluid dynamic (CFD) models use this approach to analyse flows in urban areas. In numerous situation of planning and assessment and for the near-sources region, obstacle-resolved modeling approaches are required. Large Eddy Simulations (LES) models explicitly resolve the largest eddies, and parameterize the effect of the sub grid features. Reynolds Averaged Navier Stokes (RANS) models parameterize all the turbulence, and resolve only the mean motions. CFD (large eddy simulation [LES] or Reynolds-averaged Navier-Stokes [RANS]) model can be used to explicitly resolve the urban infrastructure. Galmarini et al., 2008 and Martilli and Santiago,2008, used CFD models to estimate spatial averages required for Urban Canopy Parameters. Using CFD models good agreement in overall wind flow was reported by field Gidhagen et al. (2004).They also reported large differences in velocities and turbulence levels for identical inputs.

3.5. The Gaussian steady-state dispersion model

The Gaussian Plume Model is one of the earliest models still widely used to calculate the maximum ground level impact of plumes and the distance of maximum impact from the source. These models are extensively used to assess the impacts of existing and proposed sources of air pollution on local and urban air quality. An advantage of Gaussian modeling systems is that they can treat a large number of emission sources, dispersion situations, and a receptor grid network, which is sufficiently dense spatially (of the order of tens of meters). Figure 5 shows a buoyant Gaussian air pollutant dispersion plume. The width of the plume is determined by σy and σz, which are defined by stability classes(Pasquill 1961; Gifford Jr. 1976)

Figure 5.

A buoyant Gaussian air pollutant dispersion plume

The assumptions of basic Gaussian diffusion equations are:

  1. that atmospheric stability and all other meteorological parameters are uniform and constant throughout the layer into which the pollutants are discharged, and in particular that wind speed and direction are uniform and constant in the domain;

  2. that turbulent diffusion is a random activity and therefore the dilution of the pollutant can be described in both horizontal and vertical directions by the Gaussian or normal distribution;

  3. that the pollutant is released at a height above the ground that is given by the physical stack height and the rise of the plume due to its momentum and buoyancy (together forming the effective stack height);

  4. that the degree of dilution is inversely proportional to the wind speed;

  5. that pollutant material reaching the ground level is reflected back into the atmosphere;

  6. that the pollutant is conservative, i.e., not undergoing any chemical reactions, transformation or decay.

The spatial dynamics of pollution dispersion is described by the following type of equation in a Gaussian model:

C(x,y,z;He)=Q2πuσyσz                       [exp(y22σy2){exp((zHe)22σz2)+exp((z+He)22σz2)}]E2


C(x, y, z) : pollutant concentration at. point ( x, y, z );

U: wind speed (in the x "downwind" direction, m/s)

Σ: represents the standard deviation of the concentration in the x and y direction, i.e., in the wind direction and cross-wind, in meters;

Q: is the emission strength (g/s)

He: is the effective stack height, see below.

From the above equation, the concentration in any point ( x, y, z ) in the model domain, from a constant emission rate source, in steady state can be calculated.

Plume rise equations have been developed by Briggs (1975). The effective stack height (physical stack height plus plume rise) depends on exit velocity of gas, stack diameter, average ambient velocity, stack gas temperature and stability of atmosphere



H: height of stackTG : Temperature of exit gasQ: Volume of exit gas dθ/dz : Temperature Gradient ρ: Density of exit gas CP: Specific heat at constant pressure

Some major air pollution dispersion models in current use

  1. ADMS 3: Developed in the United Kingdom (

  2. AERMOD: Developed in the United States, (

  3. AUSPLUME: Developed in Australia, (

  4. CALPUFF: Developed in the United States, (

  5. DISPERSION2:Developed in Sweden,(

  6. ISC3: Developed in the United States, (

  7. LADM: Developed in Australia, (Physick, W.L,et al, 1994 )

  8. NAME: Developed in the United Kingdom,(

  9. MERCURE: Developed in France, (

  10. RIMPUFF: Developed in Denmark, (

Figure 6.

Air Quality Index of an Industrial Area: Orissa, India

8 regional air quality modeling leading to setting up of air quality index for an industrial area in India is given in Fig 2. This study has resulted in estimating the air assimilative capacity of the region and delineating developmental plans accordingly

3.6. Urban pollution and climate integrated modeling

Integrated air quality modelling systems are tools that help in understanding impacts from aerosols and gas-phase compounds emitted from urban sources on the urban, regional, and global climate. Piringer et al., 2007 have demonstrated that urban features essentially influence atmospheric flow and microclimate, strongly enhance atmospheric turbulence, and modify turbulent transport, dispersion, and deposition of atmospheric pollutants. Numerical weather prediction (NWP) models with increased resolution helps to visualize a more realistic reproduction of urban air flows and air pollution processes.

Integrated models thus link urban air pollution, tropospheric chemistry, and climate. Integration time required is ≥ 10 years for tropospheric chemistry studies in order to consider CH4 and O3 simulation and aerosol forcing assessment. Tropospheric chemistry and climate interaction studies extend the integration time to ≥ 100 years.

Urban air quality and population exposure in the context of global to regional to urban transport and climate change is proposed to be assessed by integrating urbanized NWP and Atmospheric Chemistry (ACT) models (Baklanov et al., 2008; Korsholm et al., 2008). A. A. Baklanov and R. B. Nuterman (2009) sugested a multi-scale modelling system which comprised of downscaling from regional to city-scale with the Environment –HIgh Resolution Limited Area Model (Enviro-HIRLAM) and to micro-scale with the obstacle-resolved Microscale Model for Urban Environment (M2UE). Meteorology governs the transport and transformations of anthropogenic and biogenic pollutants, drives urban air quality and emergency preparedness models; meteorological and pollution components have complex and combined effects on human health (e.g., hot spots, heat stresses); and pollutants, especially urban aerosols, influence climate forcing and meteorological events (precipitation, thunderstorms, etc.), thus this approach is closer to real life scenario. Examples of integrated models are Enviro-HIRLAM: Baklanov and Korsholm, 2007, WRF-Chem: Grell et al., 2005; EMS-FUMAPEX: Forecasting Urban Meteorology, Air Pollution and Population Exposure; CFD (large eddy simulation [LES] or Reynolds-averaged Navier-Stokes [RANS]) models: Galmarini et al., 2008 and Martilli and Santiago., 2008; MIT Integrated Global System Model Version 2 (IGSM2): A.P. Sokolov, C.A. Schlosser, S. Dutkiewicz, S. Paltsev, D.W. Kicklighter,H.D. Jacoby, R.G. Prinn, C.E. Forest, J. Reilly, C. Wang, B. Felzer,M.C. Sarofim, J. Scott, P.H. Stone, J.M. Melillo and J. Cohen., 2005; US EPA and NCAR communities for MM5 (Dupont et al., 2004; Bornstein et al., 2006; Taha et al., 2008), WRF models (Chen et al., 2006); THOR - an Integrated Air Pollution Forecasting and Scenario Management System: National Environmental Research Institute (NERI), Denmark.

The outline of overall methodology of FUMAPEX and MIT interactive chemistry model is shown in Figure 6 and 7. Schematic of couplings between atmospheric model and the land model components of the MIT IGSM2 is given in Figure 8.

Need of integrated models

All of these models have uncertainties associated with them. Chemical transport models, such as Gaussian plume models and gridded photochemical models, begin with pollutant emissions estimates and meteorological observations and use chemical and physical principles to predict ambient pollutant concentrations. Since these models require temporally and spatially resolved data and can be computationally intensive, they can only be used for well-characterized regions and over select time periods. Eulerian grid models are not suitable to assess individual source impacts, unless the emissions from the individual source are a significant fraction of the domain total emissions. This limitation

Figure 7.

General scheme of the FUMAPEX urban module for NWP models.

Figure 8.

Overall Scheme MIT Interactive Chemistry-Climate Model

Figure 9.

Schematic of coupling between the atmospheric model (which also includes linkages to the air chemistry and ocean models) and the land model components of the IGSM2, also shown are the linkages between the biogeophysical (CLM) and biogeochemical (TEM) subcomponents. All green shaded boxes indicate fluxes/storage that are explicitly calculated/tracked by this Global Land System (GLS). The blue shaded boxes indicate those quantities that are calculated by the atmospheric model of the IGSM2.

arises from the assumption that emissions are uniformly mixed within the grid cell, and thus do not properly address the initial growth and dispersion of the pollutants.

Lagrangian plume and puff models account for chemical processes by simple linear transformations in time. These models can track individual source impacts, thus enabling user to outline source specific air pollution control strategies. Considerable differences are observed when concentrations are compared in time and space because of uncertainties in the characterization of the direction of transport that are of the order of the actual plume width. The observed and simulated concentrations for fixed receptors, give estimates of maximum concentration values within a factor of two or three of those observed. These differences are an order of magnitude larger than those observed for estimates of secondary pollutants. Both Eulerian and Lagrangian, models are not suitable to handle inert pollutants and secondary pollutants whose concentrations depend on reaction rates and are photochemical in nature.

Receptor models, such as Positive Matrix Factorization and Chemical Mass Balance (CMB), source apportionment addresses the problem by statistical inference of source contributions to total pollution from observations of ambient air chemical composition. Mass balance methods of source apportionment use linear models with chemical composition vectors of sources as covariates. Knowledge of meteorological variables is not required but may be used to refine the analysis. Knowledge of emission sources is useful for the interpretation of results from statistical-based receptor models and is required by receptor models that use a mass balance approach. Less data and computational resource requirement by Receptor models as compared to chemical transport models, make them more convenient tool for evaluation of ambient pollutant concentrations and pollutant emission inventory. However, their utility for reactive air pollutants is uncertain and questionable. The disadvantage of CMB model arises from its assumptions. such as constant compositions of source emissions over the period of ambient and source sampling; linear additive and unreactive chemical species; identification of all sources contributing to the receptor and knowledge of their emission profile, linearly independent emission profiles.

The urban air quality models requires

  • Good net work ambient air concentrations of pollutants of concern: Geography of the urbanarea, constructed clusters, road network, location of bluidings etc play a major role in dispersion of pollutants. Thus to understand the ambient status of pollutants it is necessary to have sufficient number of monitoring locations to cover the urban sprawl of concern.

  • Micro metereology data: The wind patterns, temperature, humidity alter in urban areas according to anthropogenic activity and architecture

  • Bluilding details: To account for the effect of anthropogenic architecture falling in path of plume, its geometry is required to be known.

  • Knowledge of all sources: All sources and their emission profiles are required to be known to plan for further development in urban area and control of pollutant emission

  • Atmospheric Chemistry: All transformations of emitted chemical species, their reaction rates pathways must be known to account for observed concentration of pollutants.

  • Healthy Impacts: Models need to incorporate health effect of pollutants

None of the models available can handle all the requirements of urban air quality management. Each one focuses of one aspect and thus coupling of different models are required.


4. Further issues to be addressed

COST an intergovernmental framework for European Cooperation in Science and Technology, Europe, addressed issues related to urban air quality models in its action programmes. Cost 728 focussed on enhancing mesoscale meteorological Modeling capabilities for air pollution and Dispersion applications under larger programme of urbanization of meteorological and air quality models. The issues identified for improvements to the state of urbanization of models can be summarized as

  1. Systematic evaluation of urban land surface schemes

  2. Increasing the range of variables observed to ensure as complete a range of evaluation as possible

  3. evaluation over a broad spectrum of conditions (meteorological, morphological, geographical setting, etc.

  4. Testbeds and observatories with different objectives and dataset richness.

  5. A deeper understanding of urban PBL dynamics i.e development of long-term urban test beds in a variety of geographic regions (e.g., inland, coastal, complex terrain) and in many climate regimes, with a variety of urban core types (e.g., deep versus shallow, homogeneous versus heterogeneous).

  6. A framework to address conceptual issue of evaluation of model prediction of the flow within the canopy

  7. User friendly and multifaceted urban databases and enabling technology

  8. Developing core capabilities for advancing urban modeling and boundary layer research

  9. An open database to address issues of availability and sources of high-resolution data sets easily to all with mechanism for its maintenance, upgrading, updating, and archiving.



  1. 1. BaklanovA.KorsholmU. 2007 On-line integrated meteorological and chemical transport modelling: advantages and prospective. In: Preprints ITM 2007: 29th NATO/SPS International Technical Meeting on Air Pollution. Modelling and its Application, 24-28.09.2007, University of Aveiro, Portugal, 2134 .
  2. 2. BaklanovA.KorsholmU.MahuraA.PetersenC.GrossA.-HE. N. V. I. R. O.on-lineI. R. L. A. M.coupledmodelling.ofurban.meteorologyairpollution.Adv Sci. Res., 2 41-46, 2008.
  3. 3. BaklanovA. A.NutermanR. B.Multi-scaleatmospheric.environmentmodeling.forurban.areasAdvances.inScience.Research. 53 EOF57 EOF 5357 2009
  4. 4. BornsteinR.BalmoriR.TahaH.ByunD.ChengB.Nielsen-GammonJ.BurianS.StetsonS.EstesM.NowakD.SmithP. 2006 Modeling the effects of land-use land cover modifications on the urban heat island phenomena in Houston, Texas. SJSU Final Report to Houston Advanced Research Center for Project R-04 -0055, 127
  5. 5. BriggsG. A. 1975 Plume Rise Predictions. Lectures on Air Pollution and Environmental Impact Analysis. D. A. Haugen. Boston, MA, American Meteorology Society: 59111 .
  6. 6. CarmichaelG. R.PetersL. K. 1979 Numerical simulation of the regional transport of SO2 and sulfate in the eastern United States, Proc. 4 th Symp. on turbulence, diffusion and air pollution, AMS 337.
  7. 7. ChenF.TewariM.KusakaH.WarnerT. L. 2006 Current status of urban modeling in the community Weather Research and Forecast (WRF) model. Sixth AMS Symposium on the Urban Environment,Atlanta GA, January 2006.
  8. 8. DupontS.OtteT. L.ChingJ. K. S. 2004 Simulation of meteorological fields within and above urban and rural canopies with a mesoscale model (MM5) Boundary-Layer Meteor., 2004, 113:111-158.
  9. 9. EganB. A.RaoK. S.BassA. 1976 A three dimensional advective-diffusive model for long-range sulfate transport and transformation 7 th ITM, 697, Airlie House.
  10. 10. FriedlanderS. K. 1973 Chemical element balances and identification of air pollution sources, Environmental Science and Technology 7 235240
  11. 11. GalmariniS.J.VinuesaF.MartilliA. 2008 Relating small-scale emission and concentration variability in air quality models, Chapter 1.2, URBANIZATION OF METEOROLOGICAL AND AIR QUALITY MODELS, COST Action 728, 15 May 2008,
  12. 12. GarrattJ. R. 1978 Transfer Characteristics for a Heterogeneous Surface of Large Aerodynamic Roughness. Quart. J. R. Meteorol. Soc., 104: 491502 .
  13. 13. Garratt J.R. 1980 Surface influence upon vertical profiles in the nocturnal boundary layer Boundary-Layer Meteorology 26 1 6980 , DOI: 10.1007/BF00164331
  14. 14. Garratt J.R. 1994 Review: the atmospheric boundary layer, Earth-Science Reviews, 37, 89134 .
  15. 15. GidhagenL.JohanssonC.LangnerJ.OlivaresG. 2004 "Simulation of NOx and ultrafine particles in a street canyon in Stockholm, Sweden." Atmospheric Environment 38 14 20292044 .
  16. 16. Gifford JrF. A. 1976 "Consequences of Effluent Releases." Nuclear Safety 17 1 6886 .
  17. 17. GrellG. A.PeckhamS. E.SchmitzR.Mc KeenS. A.FrostG.SkamarockW. C.EderB. 2005 Fully coupled “online“ chemistry within the WRF model, Atmos. Environ., 39(37), 6957-6975.
  18. 18. HenryR. C. 2003 Multivariate receptor modeling by N-dimensional edge detection. Chemometrics and Intelligent Laboratory Systems 65 179-189.
  19. 19. HolmesN. al. 2005 "Spatial distribution of submicrometre particles and CO in an urban microscale environment." Atmospheric Environment 39 22 39773988 .
  20. 20. HolmesN. S.MorawskaL. 2006 A review of dispersion modelling and its application to the dispersion of particles: An overview of different dispersion models available. Atmospheric Environment, 40 5902-5928.
  21. 21. HopkeP. K. 1985 Receptor Modeling in Environmental Chemistry, Wiley, New York.
  22. 22. HopkeP. K.ed 1991 Receptor Modeling for Air Quality Management, Elsevier, Amsterdam.
  23. 23. HoskerG. 1985 1985 Clin. Phys. Physiol. Meas. 6 131 A uni-directional urethral force gauge Pavageau, M., Rafailidis, S. and Schatzmann, M. (2001) ‘A comprehensive experimental databank for the verification of urban car emission dispersion models’, International Journal of Environment and Pollution, 15, 417425
  24. 24. JacobyH. D.PrinnR. G.ForestC. E.ReillyJ.WangC.FelzerB.SarofimM. C.ScottJ.StoneP. H.MelilloJ. M.CohenJ.Report 124 July 2005
  25. 25. KorsholmU. S.BaklanovA.GrossA.SørensenJ. H.Onthe.importanceof.themeteorological.couplinginterval.indispersion.modelingduring. E. T. E. 1 Atmos. Environ., doi:10.1016/j.atmosenv.2008.11.017, 2008.
  26. 26. MartilliA.SantiagoJ. L. 2008 How to use computational fluid dynamics models for urban canopy parameterizations. Chapter 2.1, URBANIZATION OF METEOROLOGICAL AND AIR QUALITY MODELS, COST Action 728, 15 May 2008,
  27. 27. MillerS. L.AndersonM. J.DalyE. P.MilfordJ. B. 2002 Source apportionment of exposures to volatile organic compounds. I. Evaluation of receptor models using simulated exposure data. Atmospheric Environment 36 3629-3641.
  28. 28. MumovicD.CrowtherJ.StevanovicZ. 2006 Integrated Air Quality Modelling for a Designated Air Quality Management Area in Glasgow, Building and Environment, 41 12 17031712 . doi:10.1016/j.buildenv.2005.07.006
  29. 29. MumovicD.CrowtherJ. 2006 Assessing Urban Air Quality using Microscale CFD Modeling, PHOENICS CFD Newsletter Spring 2006 1010
  30. 30. OkeT. R. 1988 "The urban energy balance," Progress in Physical Geography, 12 471508 .
  31. 31. PaateroP. 2004 User’s Guide for Positive Matrix Factorization Programs PMF2 and PMF3, Part 1: Tutorial. University of Helsinki, Finland.
  32. 32. PasquillF. 1961 "The Estimation of the Dispersion of Windborne Material." Meteorology Magazine 90(1063): 3340 .
  33. 33. PavageauM.RafailidisS.SchatzmannM. 2001A comprehensive experimental databank for the verification of urban car emission dispersion models’, International Journal of Environment and Pollution, 15 417425
  34. 34. PhysickW. L.NoonanJ. A.MacGregor. J. L. 1994 LADM: a Lagrangian Atmospheric Dispersion Model [Monograph], CSIRO Division of Atmospheric Research technical paper.
  35. 35. PiringerM.PetzE.GroehnI.SchaubergerG. 2007 A sensitivity study of separation distances calculated with the Austrian Odour Dispersion Model (AODM), Atmospheric Environment, 41 725-1735.
  36. 36. RaupachM. R.ThomA. S.EdwardsI. 1980A Wind Tunnel Study of Turbulent Flow Close to Regularly Arrayed Rough Surfaces”, Boundary-Layer Meteorol. 18 373397 .
  37. 37. ReynoldsS.RothP.SeinfeldJ. 1973 Mathematical modeling of photochemical air pollution Atm.Env 7.
  38. 38. RotachM. W.1993Turbulenceclose.toa.roughurban.surfacePart. I. I.VariancesgradientsBoundary.LayerMeteorol. 66 75-92,
  39. 39. RotachM. 1995 (1995) ‘Profiles of turbulence statistics in and above an urban street canyon’, Atmospheric Environment, 29 14731486
  40. 40. ShirC. C.ShiehL. J. 1974 A generalized urban air pollution model and its application to the study of SO2-distribution in the St. Louis Metropolitan area, J. Appl. Met. 19 185204 .
  41. 41. 2002: South Asia Urban Air Quality Management Briefing. Note 4 , What Do We Know About Air Pollution?-India Case Study,
  42. 42. 2002: South Asia Urban Air Quality Management Briefing Note 5 , Impact of Traffic Management,
  43. 43. Srivastava, Anjali ; Sengupta B., Dutta S.A., ‘Source apportionment of ambient VOCs in Delhi City Science of The Total Environment’, Volume 343 1-3 , 1 May2005,207220
  44. 44. Srivastava, Anjali ; ‘Source apportionment of ambient VOCS in Mumbai city Atmospheric Environment’, Volume 38, Issuey39, December2004, pp. 6829-6843
  45. 45. StannersD.BourdeauP. .Eds. 1995 Europe’s Environment: The Dobris Assessment, European Environment Agency.Office for Publications of the European Communities, Luxemburg.
  46. 46. SokolovA. P.SchlosserC. A.DutkiewiczS.PaltsevS.KicklighterD. W. 2005Th :The MIT Integrated Global System Model (IGSM),Version 2: Model Description and Baseline Evaluation,
  47. 47. TahaH. 2008 Sensitivity of the urbanized MM5 (uMM5) to perturbations in surface properties in Houston Texas. Boundary-Layer Meteorology, 127 193-218
  48. 48. TennekesH. 1973 A model for the dynamics of the inversion above a convective boundary layer. Journal of Atmospheric Science 30, 550567
  49. 49. THOR- an Integrated Air Pollution Forecasting and Scenario Management System. Available at National Environmental Research Institute (NERI), Denmark
  50. 50. Urban Environnent Management, Tool book United Nations Environment Programme, United Nations Human Settlements Programme. Available at

Written By

Anjali Srivastava and B. Padma S. Rao

Submitted: 15 October 2010 Published: 05 July 2011