Open access

# Mathematical Modeling of Air Pollutants: An Application to Indian Urban City

Written By

P. Goyal and Anikender Kumar

Submitted: 15 October 2010 Published: 05 July 2011

DOI: 10.5772/16840

From the Edited Volume

## Air Quality - Models and Applications

Edited by Dragana Popovi?

Chapter metrics overview

View Full Metrics

## 1. Introduction

Continuous development and increase of population in the urban areas, a series of problems related to environment such as deforestation, release of toxic materials, solid waste disposals, air pollution and many more, have attracted attention much greater than ever before. The problem of air pollution in cities has become so severe that there is a need for timely information about changes in the pollution level. The air pollution dispersion is a complex problem. It covers the pollutant transport and diffusion in the atmosphere. The pollutant dispersion in the atmosphere depends on pollutant features, meteorological, emission and terrain conditions. Physical and mathematical models are developed to describe the air pollution dispersion. Physical models are small scale representations of the atmospheric flow carried out in wind tunnels. Mathematical models are divided in to statistical and deterministic models. Statistical models are based on analysis of past monitoring air quality data. Deterministic models are based on a mathematical description of physical and chemical processes taking place in the atmosphere. These models are based on mathematical equations, express conservation laws of mass, momentum and energy. Both the models are discussed in this chapter.

Statistical models are also divided into linear and non-linear models. Several studies based on the statistical models have been carried out in different regions to identify local meteorological conditions, most strongly associated with air pollutants concentration to forecast the air quality (McCollister & Willson, 1975; Aron & Aron, 1978; Lin, 1982; Aron, 1984; Katsoulis, 1988; Robeson & Steyn, 1990). Many of the previous studies (Sanchez et. al., 1990; Mantis et al., 1992; Milionis & Davies, 1994) analyzed the meteorological conditions associated with high pollutant concentration. These studies usually produced qualitative or semi quantitative results and shed light on the relation between the meteorological conditions and pollutant concentrations. Shi & Harrison, 1997 developed a linear regression model for the prediction of NOx and NO2 in London. A linear regression model was used by Cogliani, 2001 for air pollution forecast in cities by an air pollution index highly correlated with meteorological variables. Since the relation between air pollutants and meteorological variables is not linear, some non-linear models i.e., Neural Network can also be used to forecast the pollutant concentrations (Bozner et al., 1993; Comrie, 1997).

The deterministic models are divided in to Eulerian, Lagrangian and Gaussian models. In this chapter we discuss only the Eulerian analytical models. The atmospheric diffusion equation (Seinfeld, 1986) has long been used to describe the dispersion of airborne pollutants in a turbulent atmosphere. The use of analytical solutions of this equation was the first and remains the convenient way for modeling the air pollution problems (Demuth, 1978). Air dispersion models based on analytical solutions posses several advantages over numerical models as all the influencing parameters are explicitly expressed in a mathematically closed form. Analytical models are also useful in examining the accuracy and performance of numerical models. In practice, most of the estimates of dispersion are based on the Gaussian plume model, which assumes the constant wind speed and turbulent eddies with height. Hinrichsen (1986) compared a non-Gaussian model, in which wind speed and turbulence, are not constant with height and observed that non Gaussian model agreed better than Gaussian model with the observed data.

Several efforts have also been made for the development of non-Gaussian models of point and line sources, since observational studies show that the wind speed and eddy diffusivity vary with vertical height above the ground (Stull, 1988). Analytical solutions of the advection diffusion equation, with wind speed and vertical eddy diffusivity both as power function of vertical height, bounded by Atmospheric Boundary Layer (ABL) are well known for point and line sources (Seinfeld, 1986; Lin & Hildemann, 1996). Taylor’s (1921) analysis and statistical theory suggest that the eddy diffusivity depends on the downwind distance from the source (Arya, 1995). The advection diffusion equation has also solved analytically with wind speed as function of height and eddy diffusivity as a function of downwind distance from the source (Sharan & Modani, 2006). Thus in general, the eddy diffusivity should be a function of both vertical as well as downwind distance (Mooney & Wilson, 1993). Recently (Sharan & Kumar, 2009) formulate the advection diffusion equation considering the wind speed as a function of vertical height and eddy diffusivity as a function of both vertical height and downwind distance applicable only for point source release in reflecting boundary condition. However, Dirichlet (total absorption), Neumann (total reflection) and mixed boundary conditions are also appropriate for calculating the actual ground-level concentration of air pollutants. In addition to these, the few studies have been made for analytical solution of the advection diffusion equation for area sources. Park & Baik (2008) have solved the advection diffusion equation analytically for finite area source with wind speed and vertical eddy diffusivity as power function of vertical height in unbounded region

The objective of this chapter is to formulate and use the statistical (linear and non-linear) and Eulerian analytical models for prediction/forecast of air pollutants released from point, line and area sources. The analytical models are developed by using four different set of boundary conditions. The model with reflecting boundary condition is used for urban city Delhi, the capital of India and is validated by the observed values of concentration of Respirable Suspended Particulate Matter (RSPM).

## 2. Mathematical models

### 2.1. Statistical models

The main role of statistical models is to analyze past monitored air quality data. They are divided into linear and non-linear models. Linear Models as Multiple Linear Regression (MLR) can be used to make a linear empirical relationship between air pollutants and meteorological variables. The methodology of MLR is explained briefly in next section.

#### 2.1.1. Multiple Linear Regression (MLR) models

A forecast/prediction of air pollutants can be made though regression equation in which unknown variable can be expressed as a function of certain number of known variables. There is one dependent variable to be predicted in relation to the two or more independent variables. The general form of MLR can be expressed as

Y =   b 1 + b 2 X 2 + ........... + b k X k +   e E1

where Y is dependent variable, X2, X3…......., Xk are independent variables, b1, b2…......., bk are linear regression parameters and e is an estimated error term, which is obtained from independent random sampling from the normal distribution with mean zero and constant variance. The purpose of regression modeling is to estimate the b1, b2.........., bk, which can be made using minimum square error technique.

The (Eq. (1)) can also be written as

Y   =   X   b   +   e E2

where Y = [ Y 1 Y 2 : : Y n ] , X = [ 1 X 21 X 31 .. .. X k1 .. .. .. 1 X 2n X 3n .. .. X kn ] , b = [ b 1 b 2 : : b k ] and e = [ e 1 e 2 : : e n ] .

Here Y is an n x 1, X is an n x k, b is a k x 1 and e is an n x 1 matrix.

The solution of above equation can be obtained as b = ( X X ) 1 ( X Y ) using minimum square error technique. Further the F-test has been performed to determine whether a relationship exists between the dependent variables and regressors. The t-test is also performed in order to determine the potential value of each of the regressors’ variables in the regression model. The resulting model can be used to predict future values.

If the relationship between air pollutants and meteorological variables is not linear, non-linear models like Artificial Neural Network can be used for treating the non-linear relationship. The model’s characteristics are described in next section.

#### 2.1.2. Artificial Neural Network

The Artificial Neural Network (ANN) represents an alternative methodology to conventional statistical modeling because of their computational efficiency and generalization ability. ANN models are mathematical models inspired by the biological neurons. The use of ANN as mentioned in the literature is an effective alternative to more traditional statistical techniques for forecasting time series. ANN can be trained to approximate virtually any smooth, measurable and highly nonlinear functions between input and output and requires no prior knowledge to the nature of this relationship (Gardner & Dorling, 1998) and can also be trained to generalize, when presented with new and unseen data. ANNs are made up of interconnected processing elements called neurons or nodes that are arranged in the layers. These layers include an input layer, one or more hidden layers and an output layer which are connected to each neuron of the next layer by the weights. The number of hidden layer is selected based upon the problem complexity. The number of neuron in the input and output layer is problem specific. The information transfer is allowed only to the next consecutive layer. Each node of the hidden layer receives incoming signals from the nodes of the input layer. Each input value is weighted and based on its relative importance before entering the hidden layer. The total input signal NET is calculated as

NET = i W i X i E3

The total incoming signal is then passed through a non-linear transfer function F to produce the outgoing signal F (NET) of the node.

F (NET) = 1 1 + e NET E4

The output signal of a hidden node is finally passed to the nodes of the next layer (hidden or output), where a similar procedure takes place. There are several transfer functions available such as pure linear, hyperbolic tangent, sigmoid etc. Transfer function plays a key role in training process of neural network because the ANN produce different results sensitive to its transfer function (Wassermann, 1989). The process of optimizing the connection weights is known as training or learning of ANN. This is equivalent to the parameter estimation phase in the conventional statistical models. Iterative techniques are used to get the best values for connection weights by minimizing the performance function i.e. error between model output and the provided target values. The trained network is then used for the forecasting/prediction purpose. ANN has the capability to recognize the patterns in the time series data presented to it and is thus useful in many types of pattern recognition problems.

## 2.2. Deterministic models

Deterministic models are based on a mathematical description of physical and chemical processes taking place in the atmosphere. These models are divided into different categories on the basis of source characteristics as point, line and area sources or on the basis of topography of the region as flat or complex terrain. These models can also be classified on the basis of size of the field they are describing:

1. Short distance (distance from source less than 30-50 km)

2. Mesoscale models (concentration fields of the order of hundreds of kms)

3. Continental or planetary circulation models

Finally, models can also be classified on the basis of the time resolution of the concentration produced:

1. Episodic models (temporal resolution of less than an hour)

2. Short-time models (temporal resolutions grater than or equal to an hour and less than or equal to 24h)

3. Climatologically models (with resolution grater than 24h, generally seasonal or annual) (Tirabassi, 2010).

These models are divided into three categories according to different approaches as Eulerian, Lagrangian and Gaussian. Eulerian approach is based on a fixed spatial-temporal grid. The basic equation used in Eulerian air pollution dispersion models is derived from the equation of the pollutant molecular diffusion:

C t = U . C + D 2 C + S E5

where: C- concentration of pollutant in the atmospheric air,

U- wind speed vector of the components u, v, w,

D- Molecular diffusion coefficient,

S- Represents the sources and sinks of pollutant in the atmosphere i.e. its emission, removal from the atmosphere by dry and wet deposition, chemical reactions,

2- Laplacian.

Eq. (5) has been modified for the turbulent flow in the atmosphere. Modifications include the averaging procedure and closure procedure. The wind speed is expressed as the sum of two components, mean and turbulent: U = U ¯ + U and the same can be made for C as C = C ¯ + C in (Eq. (5)) and hypothesizing a wind with divergence nil:

C ¯ t = U ¯ . C ¯ - . C U ¯ + D 2 C ¯ + S ¯ E6

where C U ¯ is turbulent concentration flux.

The simplest closure method of (Eq. (6)) is a local first order closure in which K-theory is used. Assuming that turbulent concentration flux is proportional to the gradient of the average concentration, the following relation is obtained as:

C U ¯ = K C ¯ E7

where K (3x3) is turbulent diffusion coefficient.

When K tensor is diagonal, molecular diffusion is negligible and C(x, y, z, t) represents the concentration of a non-reactive pollutant ( S ¯ = S) . Thus, the (Eq. (6)) can be written as

C ¯ t = U ¯ . C ¯ + .K C ¯ + S E8

Eq. (8) can be solved for C analytically or numerically if input data for U, K and S are provided with initial and boundary conditions. Exact solution can be obtained by analytical methods, while numerical methods give the only approximate solutions. In this section we will discuss only the Eulerian analytical models.

### 2.2.1. Eulerian analytical solutions of the advection diffusion equation for point source in different atmospheric boundary conditions

Analytical solutions of the above equation are of fundamental importance in understanding and describing physical phenomena. The deterministic models for the dispersion of pollutants in atmosphere, based on the advection diffusion equation and K-theory as (Eq. (8)), can be written as:

C t + u C x + v C y + w C z = x ( K x C x ) + y ( K y C y ) + z ( K z C z ) + S E9

where Kx, Ky and Kz are the eddy diffusivities along x, y and z directions respectively.

The following assumptions are made in solution of Eq. (9):

• Steady-state condition is considered (i.e., C t = 0 ).

• The vertical velocity component (w) is neglected in comparison to horizontal velocity components (u and v).

• x-axis is oriented in the direction of mean wind (i.e., u=U, v=0).

• Downwind diffusion is neglected in comparison to transport due to mean wind

i.e.,

| U ( C x ) | | ( x ) ( K x ( C x ) ) | E10
.

Applications of these assumptions in (Eq. (9)), leads the steady state advection-diffusion equation for dispersion of a non reactive contaminate released from continuous source as (Seinfeld, 1986):

U C x = y ( K y C y ) + z ( K z C z ) E11

where x, y, and z are coordinates in the along-wind, cross wind and vertical directions respectively. C is the mean concentration of pollutants and U is the mean wind speed in downwind direction. Ky and Kz are eddy diffusivities of pollutants in the crosswind and vertical directions respectively.

Eq.(10) is solved with the following boundary conditions, in which, h is the top of the inversion/mixed layer:

Dirchlet Boundary Condition (total absorption)

C(x,y,z) = 0 at z = 0 C(x,y,z) = 0 at z = h E12

Neumann Boundary Condition (total reflection)

K z C(x,y,z) z = 0 at z = 0 K z C(x,y,z) z = 0 at z = h E13

Mixed (type-I) Boundary Condition

K z C(x,y,z) z = 0 at z = 0 C(x,y,z) = 0 at z = h E14

Mixed (type-II) Boundary Condition

C(x,y,z) = 0 at z = 0 K z C(x,y,z) z = 0 at z = h E15

The pollutant decays in cross wind direction:

C ( x ,   y ,   z )   a s   y ± E16

The pollutant is released from an elevated point source of strength Qp located at the point (0, ys, zs),

U   C   ( 0 ,   y ,   z )   =   Q p δ (y y s ) δ (z z s ) E17

where δ is the Dirac delta function.

The transport of contaminant emitted from a source primarily depends on the wind speed U. The formulations of the commonly used dispersion models in air quality studies assume wind speed to be constant. However, it is well known that wind speed increases with height in the lower part of the atmospheric boundary layer (Arya, 1999). The height dependent wind speed can be expressed as

U   ( z )   = az α , a = U(z r )z r α E18

where U (zr) is the wind speed at reference height zr and α depends on atmospheric stability.

In formulation of dispersion models, Kz is parameterized as a function of height z only (Lin & Hildemann, 1996; Park & Baik, 2008). However, based on the Taylor’s analysis and statistical theory, it is revealed that the eddy diffusivity depends on the downwind distance x (Arya, 1995).Thus, the eddy diffusivity can be a function of x and z both (Mooney & Wilson, 1993; Sharan & Kumar, 2009). The modified form of Kz(x, z) is given as:

K z ( x ,   z )   = K z (z) f(x) E19

where K z (z) is the form of eddy diffusivity depending on z and f(x) is assumed to be function of x. K z (z) is parameterized as a power law profile in z:

K z (z) = bz β , b   = K z (z r )z r β E20

where K z (z r ) is the value of K z at height z = zr and β depends on atmospheric stability.

Using Taylor’s hypothesis, the lateral eddy diffusivity can be represented by (Huang, 1979; Brown et al., 1997)

K y (x,z) = 1 2 U(z) y 2 (x) dx E21

where σ y is the standard deviation of concentration distribution in the crosswind direction.

Based on the analysis of Prairie Grass and some other historical tracer experiments of atmospheric dispersion, Irwin et al., 2007 concluded that the ground level crosswind concentration profile of dispersing plume on average is well characterized as having a Gaussian shape, which is well predicted by all atmospheric transport and diffusion models, regardless of their sophistication. Thus, by assuming the Gaussian concentration distribution in crosswind direction (Huang, 1979; Irwin et al., 2007), the steady state concentration of a pollutant released from point source in a three dimensional domain can be described as

C   ( x ,   y ,   z )   =   C   ( x ,   z ) exp( y 2 /2σ y 2 ) σ y E22

where C (x, z) is the crosswind integrated concentration.

The mathematical formulation of C(x,z) is obtained by substituting the wind profile (Eq. (17)), diffusivity profile (Eq. (15)) and boundary conditions (Eqs. (11-16)) in the advection-diffusion equation (Eq. (10)) as:

1 f ( x ) C x = b a z α z ( z β C z ) E23

provided f(x) ≠ 0,x ∈(0, ∞).

By using the separation of variables technique, the solution of the (Eq. (22)) is assumed in the form:

C ( x, z ) = X ( x ) Z ( z ) E24

which transforms the (Eq. (22)) into two following ordinary differential equations by taking 2 as the separation constant:

dX dx + f ( x ) λ 2 X = 0 E25
d dz ( z β dZ dz ) + λ 2 ( a b ) z α Z = 0 E26

Eq. (24) has the solution:

X ( x ) = A exp ( λ 2 0 x f ( s ) ds ) E27

where A is an arbitrary constant.

The solution of Eq. (25) is obtained in different boundary conditions as follows:

• Eq. (25), along with the following boundary condition corresponding to Eq. (11):

Z   = 0                         a t   z   =   0 ,   h E28

represents a Sturm-Liouville eigen value problem.

The solution of (Eq. (25)) with boundary condition (Eq. (27)) is zero for =0.

For a non-zero value of λ, the transformation of the variables

t = z ( α β + 2 2 ) E29
and
Z ( z ) = z ( 1 β 2 ) G ( t ) E30

t 2 d 2 G dt 2 + t dG dt + ( k 2 t 2 μ 2 ) G = 0 E31

where

k = a / b λ ( α β + 2 ) / 2 a n d μ = ( 1 β ) ( α β + 2 ) E32

Eq. (30) is the Bessel’s equation, whose solution can be given as:

G ( t ) = B 1 J μ ( k t ) + B 2 J μ ( k t ) E33

where Jμ and J are the Bessel’s functions of first kind of order μ and -μ respectively. From (Eqs. (28), (29) and (32)), we have

Z ( z ) = z 1 β 2 { B 1 J μ ( k z α β + 2 2 ) + B 2 J μ ( k z α β + 2 2 ) } E34

Application of the boundary condition (Eq. (27)) at z = 0 in (Eq. (33)) yields B2 = 0 and the condition at z = h (Eq. (27)) gives rise:

J μ ( k h α β + 2 2 ) = 0 E35

The corresponding eigen functions are:

Z n ( z ) = z 1 β 2 J μ ( k n z α β + 2 2 ) E36

The general solution of (Eq. (22)) is obtained by using (Eqs. (26), (34) and (35)) as:

C ( x, z ) = z 1 β 2 n = 1 A n J μ ( k n z ( α β + 2 ) / 2 ) exp ( b ( α β + 2 ) 2 k n 2 4 a 0 x f ( s ) ds ) E37

where A1, A2……An are the unknown coefficients. The (Eq. (36)) represents the concentration distribution C through the Fourier-Bessel series (Abramowitz & Stegun, 1972) corresponding to a set of eigen functions Zn.

Estimation of the coefficients An’s for crosswind integrated concentration:

The source at x = 0, Eq. (16) gives:

az α [ z 1 β 2 n = 1 A n J μ ( k n z ( α β + 2 ) / 2 ) ] = Q p δ ( z z s ) E38

The coefficients An’s are estimated using the orthogonal property of eigen functions (Abramowitz & Stegun, 1972).

Multiplying (Eq. (37)) by z ( 1 β ) / 2 J μ ( k m z ( α β + 2 ) / 2 ) for m ≥ 1, and integrating it with respect to z from 0 to h, we get:

A n = Q p α β + 2 a h α β + 2 ( z s ) 1 β 2 J μ ( k n z s ( α β + 2 ) / 2 ) J μ + 1 2 ( k n h ( α β + 2 ) / 2 ) E39

Substituting the expression for An’s for n ≥ 1 in (Eq. (36)), the final solution is obtained as:

C ( x, z ) = Q p [ α β + 2 a h α β + 2 ( z z s ) ( 1 β ) / 2 n = 1 J μ [ γ n ( z / h ) ( α β + 2 ) / 2 ] J μ [ γ n ( z s / h ) ( α β + 2 ) / 2 ] J μ + 1 2 ( γ n ) × exp ( b ( α β + 2 ) 2 γ n 2 4 a h α β + 2 0 x f ( s ) ds ) ] , E40

in which γ n is given as:

J μ ( γ n ) = 0 E41

The solution of Eq. (25) for crosswind integrated concentration by substituting the wind profile (Eq. (17)) and diffusivity profile (Eq. (18)) in (Eq. (10)) for Neumann boundary condition (Eq. (12)) is obtained from earlier study by Sharan & Kumar (2009) as:

C(x, z) = Q p [ α + 1 ah α + 1 + α β + 2 ah α β + 2 (zz s ) (1 β)/2 × n = 1 J μ n (z/h) β + 2)/2 ]J μ n (z s /h) β + 2)/2 ] J μ 2 n ) × exp( b(α β + 2) 2 γ n 2 4ah α β + 2 0 x f(s)ds ) ] , E42

in which γ n is given as:

J μ + 1 n ) =   0 E43

The (Eq. (25)), along with the following boundary conditions corresponding to (Eq. (13)):

b z β d Z d z = 0 a t   z   =   0                     a n d                   Z = 0   a t   z   =   h E44

The solution of (Eq. (25)) with boundary condition (Eq. (43)) is zero for =0.

For a non-zero value of λ, the boundary condition (Eq. (43)) at z = 0 in (Eq. (33)) yields B1 = 0 and the condition at z = h (Eq. (43)) gives rise:

J μ ( k h α β + 2 2 ) = 0 E45

The corresponding eigen functions are:

Z n ( z ) = z 1 β 2 J μ ( k n z α β + 2 2 ) E46

The general solution of (Eq. (22)) is obtained by using (Eqs. (26), (44) and (45)) as:

C ( x, z ) = z 1 β 2 n = 1 A n J μ ( k n z ( α β + 2 ) / 2 ) exp ( b ( α β + 2 ) 2 k n 2 4 a 0 x f ( s ) ds ) E47

where A1, A2…… are the unknown coefficients.

The (Eq. (46)) represents the concentration distribution C as the Fourier-Bessel series (Abramowitz & Stegun, 1972) corresponding to a set of eigen functions Zn.

Estimation of the coefficients An’s for crosswind integrated concentration

The source condition at x = 0 (Eq. (16)), gives:

az α [ z 1 β 2 n = 1 A n J μ ( k n z ( α β + 2 ) / 2 ) ] = Q p δ ( z z s ) E48

The coefficients An’s are estimated using the orthogonal property of eigen functions (Abramowitz & Stegun, 1972).

Multiplying (Eq. (47)) by z ( 1 β ) / 2 J μ ( k m z ( α β + 2 ) / 2 ) , m ≥ 1, and integrating it with respect to z from 0 to h, we get:

A n = Q p α β + 2 a h α β + 2 ( z s ) 1 β 2 J μ ( k n z s ( α β + 2 ) / 2 ) J μ + 1 2 ( k n h ( α β + 2 ) / 2 ) E49

Substituting the expression for An’s, n ≥ 1 into (Eq. (48)), the final solution of an elevated point source at crosswind is obtained as:

C ( x, z ) = Q p [ α β + 2 a h α β + 2 ( z z s ) ( 1 β ) / 2 n = 1 J μ [ γ n ( z / h ) ( α β + 2 ) / 2 ] J [ γ n ( z s / h ) ( α β + 2 ) / 2 ] J μ + 1 2 ( γ n ) × exp ( b ( α β + 2 ) 2 γ n 2 4 a h α β + 2 0 x f ( s ) ds ) ] , E50

in which γ n is given as:

J μ ( γ n ) = 0 E51

The (Eq. (25)), along with the following boundary conditions corresponding to (Eq. (14)):

Z   =   0   a t   z   =   0   a n d b z β d Z d z = 0 a t   z = h E52

The solution of (Eq. (25)) with (Eq. (51)) is zero for λ = 0.

For a non-zero value of λ, the boundary condition (Eq. (51)) at z = 0 in (Eq. (33)) yields B2 = 0 and the condition at z = h (Eq. (51)) gives rise:

J μ 1 ( k h α β + 2 2 ) = 0 E53

The corresponding eigen functions are:

Z n ( z ) = z 1 β 2 J μ ( k n z α β + 2 2 ) E54

The general solution of (Eq. (22)) is obtained by using (Eqs. (25), (52) and (53)) as:

C ( x, z ) = z 1 β 2 n = 1 A n J μ ( k n z ( α β + 2 ) / 2 ) exp ( b ( α β + 2 ) 2 k n 2 4 a 0 x f ( s ) ds ) E55

where A1, A2…… are the unknown coefficients. The Eq. (54) represents the concentration distribution C as the Fourier-Bessel series (Abramowitz & Stegun, 1972) corresponding to a set of eigen functions Zn.

Estimation of the coefficients An’s for crosswind integrated concentration

The source condition at x = 0 (Eq. (16)), gives:

az α [ z 1 β 2 n = 1 A n J μ ( k n z ( α β + 2 ) / 2 ) ] = Q p δ ( z z s ) E56

The coefficients An’s are estimated using the orthogonal property of eigen functions (Abramowitz & Stegun, 1972).

Multiplying (Eq. (55)) by z ( 1 β ) / 2 J μ ( k m z ( α β + 2 ) / 2 ) , m ≥ 1, and integrating it with respect to z from 0 to h, we get:

A n = Q p α β + 2 a h α β + 2 ( z s ) 1 β 2 J μ ( k n z s ( α β + 2 ) / 2 ) J μ 2 ( k n h ( α β + 2 ) / 2 ) E57

Substituting the expression for An’s, n ≥ 1 into (Eq. (54)), the final solution of an elevated point source at crosswind is obtained as:

C ( x, z ) = Q p [ α β + 2 a h α β + 2 ( z z s ) ( 1 β ) / 2 n = 1 J μ [ γ n ( z / h ) ( α β + 2 ) / 2 ] J μ [ γ n ( z s / h ) ( α β + 2 ) / 2 ] J μ 2 ( γ n ) × exp ( b ( α β + 2 ) 2 γ n 2 4 a h α β + 2 0 x f ( s ) ds ) ] , E58

in which γ n is given as:

J μ 1 ( γ n ) = 0 E59

### 2.2.2. Point source models

The steady state three dimensional solution of a point source can be obtained from (Eq. (21)) and f(x), which is expressed as a linear function of downwind distance, i.e., f(x) = γ U ¯ x (Sharan & Modani, 2006), in different boundary conditions:

• The solution in Dirichlet Boundary condition is given as

C(x, y, z) = Q p σ y [ α β + 2 ah α β + 2 (zz s ) (1 β)/2 × n = 1 J μ n (z/h) β + 2)/2 ]J μ n (z s /h) β + 2)/2 ] J μ + 1 2 n ) × exp( U ¯ ( x-x s ) 2 β + 2) 2 γ n 2 8ah α β + 2 ) ] × exp [ (y y s ) 2 y 2 ] , E60

where μ = (1 β)/(α β + 2), J μ is the Bessel function of order μ and γ n 's are obtained from the equation:

J μ n ) =   0 E61

The solution in Neumann Boundary condition is given as

C(x, y, z) = Q p σ y [ α + 1 ah α + 1 + α β + 2 ah α β + 2 (zz s ) (1 β)/2 × n = 1 J μ n (z/h) β + 2)/2 ]J μ n (z s /h) β + 2)/2 ] J μ 2 n ) × exp( U ¯ ( x-x s ) 2 β + 2) 2 γ n 2 8ah α β + 2 ) ] × exp [ (y y s ) 2 y 2 ] , E62

where μ = (1 β)/(α β + 2), J μ is the Bessel function of order is μ and γ n 's are obtained from the equation:

J μ + 1 n ) =   0 E63

The solution in Mixed (type-I) Boundary condition is given as

C(x, y, z) = Q p σ y [ α β + 2 ah α β + 2 (zz s ) (1 β)/2 × n = 1 J μ n (z/h) β + 2)/2 ]J μ n (z s /h) β + 2)/2 ] J μ + 1 2 n ) × exp( U ¯ ( x-x s ) 2 β + 2) 2 γ n 2 8ah α β + 2 ) ] × exp [ (y y s ) 2 y 2 ] , E64

where μ = (1 β)/(α β + 2), J μ is the Bessel function of order μ and γ n 's are obtained from the equation:

J μ n ) =   0 E65

The solution in Mixed (type-II) Boundary condition is given as

C(x, y, z) = Q p σ y [ α β + 2 ah α β + 2 (zz s ) (1 β)/2 × n = 1 J μ n (z/h) β + 2)/2 ]J μ n (z s /h) β + 2)/2 ] J μ 2 n ) × exp( U ¯ ( x-x s ) 2 β + 2) 2 γ n 2 8ah α β + 2 ) ] × exp [ (y y s ) 2 y 2 ] , E66

where μ = (1 β)/(α β + 2), J μ isthe Bessel function of order is μ and γ n 's are obtained from the equation:

J μ 1 n ) =   0 E67

### 2.2.3. Correspondence between present solutions to the earlier solutions

In present solution, wind speed is a function of z and vertical eddy diffusivity is considered a function of x and z as Kz(x, z) = bz β f(x). If f(x) is equal to 1.0, the solutions of Eqs. (59), (61), (63) and (65) are become same as the solutions obtained by Lin & Hildemann (1996). When wind speed is power law function of z and vertical eddy diffusivity is function of x only, i.e. Kz(x, z) = f(x) = K(x) (for β = 0 and b =1.0), the expression for the two dimensional concentration i.e., Eq. (41) is same as solution by Sharan & Modani (2006).

### 2.2.4. Line source models

A line source can be considered as a superposition of point sources. The solution for finite line source can be obtained by integrating point source solution from ys = y1 to y2 with unit source strength Q with the same γ n 's as in point source in different boundary conditions:

• The solution in Dirichlet Boundary condition is given as

C(x, y, z) = Q 2 [ α β + 2 ah α β + 2 (zz s ) (1 β)/2 × n = 1 J μ n (z/h) β + 2)/2 ]J μ n (z s /h) β + 2)/2 ] J μ + 1 2 n ) × exp( U ¯ ( x-x s ) 2 β + 2) 2 γ n 2 8ah α β + 2 ) ] × [ erf( y y 1 2 σ y ) erf( y y 2 2 σ y ) ] E68

The solution in Neumann Boundary condition is given as

C(x, y, z) = Q 2 [ α + 1 ah α + 1 + α β + 2 ah α β + 2 (zz s ) (1 β)/2 × n = 1 J μ n (z/h) β + 2)/2 ]J μ n (z s /h) β + 2)/2 ] J μ 2 n ) × exp( U ¯ ( x-x s ) 2 β + 2) 2 γ n 2 8ah α β + 2 ) ] × [ erf( y y 1 2 σ y ) erf( y y 2 2 σ y ) ] E69

The solution in Mixed (type-I) Boundary condition is given as

C(x, y, z) = Q 2 [ α β + 2 ah α β + 2 (zz s ) (1 β)/2 × n = 1 J μ n (z/h) β + 2)/2 ]J μ n (z s /h) β + 2)/2 ] J μ + 1 2 n ) × exp( U ¯ ( x-x s ) 2 β + 2) 2 γ n 2 8ah α β + 2 ) ] × [ erf( y y 1 2 σ y ) erf( y y 2 2 σ y ) ] E70

The solution in Mixed (type-II) Boundary condition is given as

C(x, y, z) = Q 2 [ α β + 2 ah α β + 2 (zz s ) (1 β)/2 × n = 1 J μ n (z/h) β + 2)/2 ]J μ n (z s /h) β + 2)/2 ] J μ 2 n ) × exp( U ¯ ( x-x s ) 2 β + 2) 2 γ n 2 8ah α β + 2 ) ] × [ erf( y y 1 2 σ y ) erf( y y 2 2 σ y ) ] , E71

where erf is the error function defined by

erf (a) = 2 π 0 a e t 2 dt

The solution of an infinite line source of source strength Q can be obtained by integrating point source along the crosswind direction ys = - to , which is different than the finite line source formulation, since exp( ( y y s ) 2 /2σ y 2 ) σ y dy s =1.

### 2.2.5. Area source models

The concentration of air pollutants due to a finite area source at (x, y, z) is calculated as a superposition of finite line sources extending from x1 to x2 in x direction is obtained as

C(x, y, z) = x 1 x 2 C(x x s , y, z)dx s E72

where C(x-xs, y, z) is equivalent to the finite line source as obtained in section 2.2.4. However, the source strength Q is replaced by Qa in (Eqs. 67, 68, 69 and 70) with same γ n 's as in point source.

For an infinite area source with uniform strength Qa, the solution is obtained as:

C(x, z) = x 1 x 2 C(x x s , z)dx s E73

where C(x-xs, z) is equivalent to the crosswind integrated concentration. However, the source term Qp is replaced by Qa in (Eqs. 39, 41, 49 and 57) with same γ n 's as in point source.

Eqs. 71 and 72 are the solutions for finite and infinite area sources and can be solved numerically.

If wind speed and vertical eddy diffusivity are considered as function of z, i.e. Kz(x, z) = bz β (f(x) =1.0) with unbounded region (h) and ground level source (zs0) in (Eq. (71)), the expression obtained for the concentration is same as in Park & Baik, 2008.

## 3. Case study of Delhi

A case study of Delhi has been attempted through the application of two different categories namely statistical and analytical models to fulfill the objective of the chapter.

### 3.1. Application of statistical models

First of all the daily Air quality index (AQI) as a comprehensive assessment of air quality concentration of criteria pollutants namely Respirable Suspended Particulate Matter (RSPM), Sulphur dioxide (SO2), Nitrogen dioxide (NO2) and Suspended Particulate Matter (SPM) has been calculated at ITO (a busiest traffic intersection) for a period of seven years (2000-2006), monitored continuously by Central Pollution Control Board (CPCB). A method of US Environmental Protection Agency (USEPA) has been used for estimating the AQI, in which the sub-index and breakpoint pollutant concentrations depend on Indian National Ambient Air Quality Standard (NAAQS). There are primarily two steps involved in formulating an AQI: first the formation of sub-indices of each pollutant, second the aggregation (breakpoints) of sub indices. The Breakpoint concentration of each pollutant, used in calculation of AQI, is based on Indian NAAQS and results of epidemiological studies indicating the risk of adverse health effects of specific pollutants. It has been noticed that different breakpoint concentrations and different air quality standards has been reported in literature (Environmental Protection Agency, 1999). In India, to reflect the status of air quality and its effects on human health, the range of index values has been designated as “Good (0-100)”, “Moderate (101-200)”, “Poor (201-300)”, “Very Poor (301-400)” and “Severe (401-500)” (Nagendra et al. (2007)) as shown in Table 1.

The formula (EPA, (1999)) used to calculate AQI for four criteria pollutants RSPM, SO2, NO2 and SPM from 2000-2006 is given below:

I P = [ ( I Hi I Lo ) ( BP Hi BP Lo ) ] ( C P BP Lo ) + I Lo E74

where IP = the AQI for pollutant ‘p’,

CP = actual ambient concentration of the pollutant ‘p’,

BPHi = the breakpoint in Table 1 that is grater than or equal to Cp,

BPLo = the breakpoint in Table 1 that is less than or equal to Cp,

IHi = the subindex value corresponding to BPHi,

ILo = the sub index value corresponding to BPLo.

 SI.No.Index valuesDescriptorSO2 (24-h avg.)NO2 (24-h avg.)RSPM (24-h avg.)SPM (24-h avg.) 10-100Gooda 0-800-800-1000-200 2101-200Moderateb 81-36781-180101-150201-260 3201-300Poorc 368-786181-564151-350261-400 4301-400Very Poord 787-1572565-1272351-420401-800 5401-500Severe e "/1572"/1272"/420"/800 "/1572"/1272"/420"/800 Index values Descriptor SO2 (24-h avg.) NO2 (24-h avg.) RSPM (24-h avg.) SPM (24-h avg.) 1 0-100 Gooda 0-80 0-80 0-100 0-200 2 101-200 Moderateb 81-367 81-180 101-150 201-260 3 201-300 Poorc 368-786 181-564 151-350 261-400 4 301-400 Very Poord 787-1572 565-1272 351-420 401-800 5 401-500 Severe e "/1572"/1272"/420"/800 "/1572 "/1272 "/420 "/800

### Table 1.

ll the values of SO2, NO2, RSPM and SPM are in µg/m3.a. Good: Air quality is acceptable; however, for some pollutants there may be a moderate health concern for a very small number of people.b. Moderate: Members of sensitive groups may experience health effects.c Poor: Members of sensitive groups may experience more serious health effects.d Very poor: Triggers health alter, everyone may experience more serious health effects.e Severe: Triggers health warnings of emergency conditions.Propose sub-index and breakpoint pollutant concentration for Indian-AQI.

The overall AQI is now determined on the basis of the AQI for above pollutant ‘p’ and highest among them is declared as the overall AQI for that day.

The above estimated daily AQI along with meteorological variables like daily maximum temperature (tmax), minimum temperature (tmin), daily temperature range (difference between daily maximum and minimum temperature, trange), average temperature (tavg), wind speed (wsp), wind direction index (wdi), relative humidity (rh), vapor pressure (vp), station level pressure (slp), rainfall (rf), sunshine hours (ssh), cloud cover (cc), visibility (v) and radiation (rd), monitored at Safdarjung airport by Indian Meteorological Department (IMD), Delhi, have been used in statistical models to forecast daily AQI one day in advance.

This study has been carried out for four different seasons namely summer (March, April, May), monsoon (June, July, August), post monsoon (September, October, November) and winter (December, January, February). The location ITO has been chosen in the present study due to various reasons: (i.) a busiest traffic intersection, (ii.) air pollutants concentration of RSPM, SO2, NO2 and SPM is monitored continuously and (iii.) the meteorological station Safdarjung airport is within 10 km radius. These models include previous day’s AQI and meteorological variables as input and yield daily forecasting of AQI. The input and output is normalized between -1 to +1 using the minimum and maximum of the time series before any preprocessing. Forecast of daily AQI in Delhi has been obtained by two models MLR and ANN independently. The performance of both the models has been assessed with respect to the statistical parameters.

#### 3.1.1. AQI by MLR model

The following MLR equations for different seasons are resulted through training of AQI and meteorological data of 2000-2005 using SPSS software for summer, monsoon, post monsoon and winter respectively:

[ A Q I ]   =   0.0478   +   0.504   x   [ A Q I d 1 ]   0.079   x   [ r h ]   + 0.126   x   [ t m a x ]     0.068   x   [ c c ] E75
[ A Q I ]   =   0.181   +   0.599   x   [ A Q I d 1 ]     0.282   x   [ r h ]     0.128   x   [ v ]     0.155   x   [ t m i n ] E76
[ A Q I ]   =   0.324   +   0.537   x   [ A Q I d 1 ]   +   0.573 x   [ s l p ]     0.112 x   [ v p ]   + 0.070   x   [ s s h ]     0.135   x   [ v ]   +   0.066   x   [ t m a x ] E77
[ A Q I ]   =   0.171   +   0.503   x   [ A Q I d 1 ]     0.191 x   [ v ]     0.115   x   [ c c ]   0.169 x   [ w s p ]   0.157   x   [ r h ]   + 0.151   x   [ r f ] E78

The previous day’s AQI is the common variable in all four equations. The above equations have also been used to forecast the daily AQI of 2000-2005 for all four seasons, which have been compared with the observed AQI, used as trained data of the corresponding seasons during 2000-2005 and are shown graphically Figs.1 (a), (b) (c) and (d).

The Fig. 1 reflects that the trained values of AQI are matching well with observed values (calculated directly from formulation Eq. (73)). The same set of Eqs. (74-77) have been used for forecasting the daily AQI in all four seasons of the year 2006 which have been shown graphically in Figs. 2 (a), (b), (c) and (d) for summer, monsoon, post monsoon and winter seasons respectively. The observed values of AQI of year 2006 for each season have also been plotted in the Fig. 2 in order to validate the forecasted values of AQI.

The quantitative analysis of comparison of forecasted and observed values of AQI has been made through statistical parameters, which are summarized in Table 2. The NMSE and coefficient of determination (R2) are found as (0.0094, 0.5718) in summer season which are followed by (0.0369, 0.5247) in winter; (0.0629, 0.3913) in monsoon and (0.1287, 0.3021) in post monsoon seasons, showing good performance of the model as the ideal values of NMSE and R2 are 0 and 1 respectively. The values of R2 in four different seasons reflect that the model’s forecasted and observed AQI could be correlated explained by the selected input variables as approximately 57% in summer, 52% in winter, 39% in monsoon and 30% in post monsoon seasons. However, the fractional bias shows the under-prediction in all the seasons except Monsoon.

 S.N. Season 2006 RMSE NMSE coefficient ofdetermination Fractional Bias 1 Summer 33.13 0.0094 0.5718 0.0126 2 Monsoon 72.99 0.0629 0.3913 -0.1246 3 Post Monsoon 99.18 0.1287 0.3021 0.1972 4 Winter 64.39 0.0369 0.5247 0.1573

### Table 2.

Comparison of MLR model predicted and observed values in years 2000-2005 and year 2006.

#### 3.1.2. AQI by ANN model

The ANN model for different seasons is developed through training of AQI and meteorological data of 2000-2005 using MATLAB. The weights of a network are iteratively modified to minimize the total mean squared error between the desired target and actual output values. The above model has also been trained using observed daily AQI (same values as in MLR) of 2000-2005 for all four seasons. The comparison of observed and trained AQI has been shown graphically in Figs.3 (a), (b) (c) and (d).

The Fig. 3 reflects that the trained values of AQI are showing the same trend as observed values. The same trained architectures have been used for forecasting the daily AQI in all four seasons of the year 2006, which have been shown graphically in Figs. 4 (a), (b), (c) and (d) for summer, monsoon, post monsoon and winter seasons respectively. The observed values of AQI of year 2006 for each season have also been plotted in the Fig. 4 in order to validate the forecasted values of AQI.

The comparison of forecasted and observed values of AQI has been made through statistical parameters in Table 3, which revels that the NMSE and coefficient of determination (R2) are found as (0.0118, 0.4976) in summer season which are followed by (0.0277, 0.3986) in winter; (0.0540, 0.3619) in monsoon and (0.0718, 0.3816) in post monsoon seasons, showing good performance of the model as the ideal values of NMSE and R2 are 0 and 1 respectively. The values of R2 in four different seasons reflect that the model’s forecasted and observed AQI could be correlated explained by the selected input variables as approximately 50% in summer, 40% in winter, 38% in post monsoon and 36% in monsoon seasons. However, the fractional bias shows the under-prediction in all the seasons except Monsoon.

 S.N. Season 2006 RMSE NMSE Coefficient ofdetermination Fractional Bias 1 Summer 36.88 0.0118 0.4976 0.0269 2 Monsoon 64.27 0.0540 0.3619 -0.0224 3 Post Monsoon 76.44 0.0718 0.3816 0.1355 4 Winter 56.52 0.0277 0.3986 0.1330

### Table 3.

Comparison of ANN model predicted and observed values in years 2000-2005 and year 2006

#### 3.1.3. Comparison of MLR and ANN models

On the basis of above exercise, it can be seen that ANN model is giving better results than MLR model with respect to RMSE and NMSE as shown in Tables 2 and 3, although, the R2 is showing almost same values in both the models throughout all four seasons.

### 3.2. Application of analytical models

The above discussed analytical model with Neumann boundary condition is applied to simulate the hourly concentration of RSPM in the month of January (representative of winter season) 2008 due to point, line and area sources. The most important parameters of the models are emission inventory and meteorological variables, which are pre-process according to the model’s requirement.

A gridded emission inventory of RSPM has been developed over an area of 26 km x 30 km of Delhi. The total area has been divided into 195 square grids of size 2 km x 2 km. Emission of RSPM has been estimated in each grid due to all anthropogenic sources viz., domestic, industries, power plants and vehicles for the year 2008 using the primary and secondary data. The emission of RSPM from domestic sector has been calculated on the basis of fuel consumption data and the emission factor of the corresponding fuel. The emission estimation from industrial sources has been made by using the data obtained from concerned agencies / department. The emission of RSPM from power plants has been calculated on the basis of information given by Delhi Pollution Control Committee (DPCC). However, estimation of vehicular emission in each grid are made by using the following mathematical formulation, based on total number of registered vehicles, emission factors, vehicle kilometer traveled (VKT) of each type of vehicle, different type of roads:

E i =   N i x   e f i x   V K T i E79

where Ni is the number of vehicle of ith category, Ei is the emission of pollutant by vehicle of ith category (g/day), efi is emission factor of ith category vehicle for the pollutant (g/km-vehicle). VKTi is vehicle kilometer travel per day for ith category vehicle (km/day).

Further, the total emission of pollutant has been calculated by summing the emissions from each category of vehicles:

E = i E i E80

The emissions of RSPM due to vehicular sources in Delhi for the year 2008 have been apportioned into each grid w. r. t. to the road lengths.

The gridded emission inventory shows the spatial distribution of emissions of RSPM due to all types of sources (domestic, industries, power plants and vehicles) in Fig. 5. The emission from domestic sources has been distributed uniformly in all the grids and the emission from industries and traffic intersections have been apportioned according to their locations. The emissions of RSPM due to three power plants namely Indraprastha, Badarpur and Rajghat are superimposed in the grids as per their locations. It has been found that vehicles are the major source, contributing 57% of RSPM among all the estimated sources, which is followed by power plants 29% of total RSPM. While the contributions of domestic and other industrial sources are 8% and 6% respectively. The Fig.5 reflects that grids, in which power plants are located in addition to traffic intersections and industries, have higher values of emissions compared to others. This emission inventory is used as input to the model.

From an atmospheric pollution perspective, the most important season in Delhi is the winter lasting from December to February. This period is dominating by cold, dry air and ground based inversion with low wind conditions (<1 ms-1), which increase the concentration of pollutants (Anfossi et al., 1990). For practical reasons, the January month of 2008 is used in this case study. The hourly meteorological data, measured at Safdarjung Airport from India Meteorological Department (IMD) is used as second input file to the models. There are some traditional methods, which are used to determination, the stability classes to produce the input in the model. The some traditional methods based on (i.) the wind velocity, sun radiation intensity and cloud cover; (ii.) measurement of the wind direction fluctuations; (iii.) the vertical temperature gradient; (iv.) the Richardson number can used to determination, the stability classes.

 Day Night Incoming Solar Radiation Cloud Cover Surface wind speed at 10m (m/s) Strong Moderate Slight Mostly MostlyOvercast Clear <2 A A-B B E F 2-3 A-B B C E F 3-5 B B-C C D E 5-6 C C-D D D D "/6 C D D D D

### Table 4.

ource: D.B. Turner. Workbook of Atmospheric Dispersion Estimate. USEPA 999-AP-26. U.S. Environmental Protection Agency, Washington, D.C, 1969.Note: A, strongly unstable; B, unstable; C, weakly unstable; D, neutral; E, weakly stable and F, stable.Table 4. Stability Classification based on wind velocity, sun radiation intensity and cloud cover.

 Stability Class Stability class of Pasquill σ θ (in degree) Vertical Temperature Gradient (0C/m 10-2) Richardson Number at 2 m Very unstable A 25.0 <-1.9 -0.9 Moderately unstable B 20.0 -1.9 to -1.7 -0.5 Slightly unstable C 15.0 -1.7 to -1.5 -0.15 Neutral D 10.0 -1.5 to -0.5 0 Slightly stable E 5.0 -0.5 to 1.5 0.4 Moderately stable F 2.5 1.5 to 4.0 0.8

### Table 5.

ource: Zannetti, P., Air Pollution Modelling, Computational Mechanics Publications, Southampton, U.K. and Van Nostrand Reinhold, New York, 1990.Note: σ θ is the standard deviation of horizontal wind direction. Classification of Atmospheric Stability based on wind direction fluctuations, vertical temperature gradient and the Richardson number.

In this study atmospheric stability is measured on the basis of surface wind speed, cloud cover and solar insolation (strong, moderate, slight). These stability are classified according to Pasquill's stability classes of A, B, C, D, E and F, which range from extremely unstable to extremely stable as given by Turner (1969). The dispersion parameters have also been calculated based on the stability parameters as discussed in section 2.

In the above solutions discussed in section 2, the exponent α has relationship with Pasquill’s Stability Classes (Hanna et al., 1982) and β = 1 α is based on the Schmidt’s conjugate law.

 Pasquill’s Stability Classes A B C D E F α for Urban 0.15 0.15 0.20 0.25 0.40 0.60

### Table 6.

Approximate Correspondence between Pasquill’s Stability Classes and alpha ( α ) (Hanna, 1982).

In these formulations f(x) is expressed as a linear function of downwind distance as f(x) = γ U ¯ x . In this function U ¯ is average velocity and γ is turbulence parameter. This turbulence parameter γ is parameterized as the square of turbulent intensity using Taylor statistical theory of diffusion γ = ( σ w U ¯ ) 2 (Arya, 1995; 1999). Turbulent intensity can be expressed as ( σ w U ¯ ) = tan ( σ φ ) , where σ φ is the standard deviation of vertical wind direction in radians. For small σ φ , the turbulent intensity which is also depends on atmospheric stability is approximated as (Arya, 1999): ( σ w U ¯ ) = σ φ .

 Pasquill’s Stability Classes A B C D E F σ φ ≥ 11.5 10.0-11.5 7.8-10.0 5.0-7.8 2.4-5.0 < 2.4

### Table 7.

Approximate Correspondence between Pasquill’s Stability Classes and Turbulence parameters σ φ (Arya, 1999).

In the above solutions, standard deviation of concentration distribution in crosswind direction is represented by a power of downwind distance y = Rxr (Seinfeld, 1986). These R and r are constants as depending on the atmospheric stability.

 Source Coefficient Stability Class A B C D E F (Turner, 1969; Martin, 1976) R r 0.443 0.894 0.324 0.894 0.216 0.894 0.141 0.894 0.105 0.894 0.071 0.894

### Table 8.

Approximate Correspondence between Pasquill’s Stability Classes and Gaussian plume dispersion parameters (Seinfeld, 1986).

The 24 hourly averaged ground level concentration of RSPM from point, line and area sources for Jan 2008 has been obtained by using the emission inventory and meteorological data as input parameters to the models. The power plants are considered as the point sources. Most of the analytical methods for predicting the concentration from stack need the plume rise height ( Δh ) of stack (power plants). There are numerous methods for calculating the plume rise height; in this study it has been calculated, given by Carpenter et al., 1970:

Δh = 114 G ( I ) 1/3 u E81

where I = gVsd2(Ts-Ta)/4Ta, m4s-3; G=1.58-41.4 ( Δθ Δz ) , dimensionless; ( Δθ Δz ) = potential temperature gradient, oK/m; the constant 114 in this equation has unit of m2/3; Vs is the stack gas exit velocity, in m/s; d is the stack exit diameter in m; u is the speed at the stack exit in m/s; Ts is the stack gas temperature; Ta is the environmental temperature and g is the gravitational acceleration, m/s2. This method is especially useful because the potential temperature gradient factor is used for adjusting the different stability conditions.

The spatial distribution of RSPM concentration obtained from models (Eqs. (61), (67) and (71)) with background concentration 40 µg/m3 for Delhi (Kansal et al., 2011) is shown in Fig. 2 in the form of isopleths, which indicates the hot spots of RSPM ranging from 300-800 µg/m3 at different locations namely I.T.O., Nizamuddin, Badarpur power station, I.P. powers station, Rajghat power station, AIIMS, Sirifort and Dhaula Khuan, which are major traffic intersections or power plants. This figure also shows the higher concentration near the Shahzada Bagh, which has maximum number of industries.

The concentrations of RSPM, cumulative sum of concentrations due to point, line and area sources, is evaluated against observed concentrations obtained from CPCB and NEERI, at different locations, as shown in Table 9.

 Location Observed (µg/m3) Model predicted(µg/m3) National AmbientAir Quality Standard (µg/m3) Pitampura 353.00 411.30 100.00 Sirifort 374.00 449.45 100.00 Janakpuri 279.00 345.49 100.00 Shahzada Bagh 478.74 603.01 100.00 Sarojini Nagar 363.77 401.28 100.00

### Table 9.

Comparison of 24 hourly averaged predicted and observed concentration of RSPM at different locations in Delhi in Jan 2008.

The table 9 shows that the models predictions are higher than observed as well as NAAQS. However, these values are well within a factor of two of observed values, which satisfies the criteria of Chang and Hanna (1982) for assessing the performance of the model with conclusion that the models are performing satisfactory.

## 4. Conclusion

In the present study, the statistical (linear and non-linear) and analytical dispersion models of air pollutants released from point, line and area sources are discussed. Air quality index has been forecasted using MLR and ANN models. Performance of both the modes has been compared and observed that ANN model is doing better than MLR.

The analytical models are formulated by considering the wind speed as a power law profile of vertical height above the ground and vertical eddy diffusivity as an explicit function of downwind distance and vertical height in different boundary conditions. A case study of Delhi has been made for predict 24 hourly concentration of RSPM through the application models with Neumann boundary condition in the month of Jan 2008. The input parameters namely emission inventory and meteorological variables are pre-processed as per the requirement of the model. The different types of primary and secondary sources of RSPM due to vehicular, domestic, industries and power plants have been used in emission inventory for year 2008. Some traditional methods for determining the meteorological field are also discussed in this chapter, which are used as an input to these analytical models. The analytical models are evaluated with observed concentration at different locations in Delhi obtained from CPCB and NEERI, which show that the concentration levels obtained from the models are always high in comparison to the observed values and NAAQS. However, the models are performing satisfactory. Although the present models have the limitation as the longitudinal diffusion is neglected in comparison to the advection and are not considering the wind directions at different vertical heights. These models can be used for other Indian urban cities.

## 5. Appendix A

The statistical measures, which have been used for statistical evaluation of the performance of models has been given by Chang & Hanna (2004) as follows:

### 5.1. Coefficient of correlation (R)

Coefficient of correlation (R) is relative measure of the association between the observed and predicted values. It can vary from 0 (which indicates no correlation to +1 (which indicates perfect correlation). A value of R close to 1.0 implies good agreement between the observed and predicted values i.e. good model performance.

R = ( C o C ¯ o ) ( C p C ¯ p ) ¯ σ C p σ C o E82

### 5.2. Coefficient of determination (R2)

Coefficient of determination (R2), which is the square of coefficient of correlation, determines the proportion of variance that can be explained by the model.

### 5.3. Root Mean Square Error (RMSE)

RMSE, is a measure of the differences between values predicted by a model and the observed values and is expressed as follows:

RMSE = ( C o C p ) 2 ¯ E83

### 5.4. Normalized Mean Square Error (NMSE)

NMSE, as a measure of performance, emphasizes the scatter in the entire data set and is defined as follows:

NMSE = ( C o C p ) 2 ¯ C ¯ o . C ¯ p E84

The normalization by C ¯ o . C ¯ p ensures that NMSE will not be biased towards models that over predict or under predict. Ideal value for NMSE is zero. Smaller values of NMSE denote better model performance.

### 5.5. Fractional Bias (FB)

It is a performance measure known as the normalized or fractional bias of the mean concentrations:

FB = ( C ¯ o C ¯ P ) 0.5 ( C ¯ o + C ¯ P ) E85

where:

Cp: model predictions,

Co: observations,

Overbar ( C ¯ ) : average over the dataset, and

σ C : standard deviation over the data set.

## References

1. 1. Abramowitz M. Stegun I. A. 1972 Handbook of Mathematical Functions with Formulas, Graphs and Mathematical Tables. 9th printing. Dover Publications, New York.
2. 2. Anfossi D. Brisasca G. Tinarelli G. 1990 Simulation of atmospheric diffusion in low wind speed meandering conditions by a Monte Carlo dispersion model. IL Nuovo Climento. 13 C, 995 1006 .
3. 3. Aron R. Aron I. M. 1978 Statistical forecasting models: I. Carbon monoxide concentrations in the Los Angeles basin. Journal of Air Pollution Control Association 28 681 684 .
4. 4. Aron R. 1984 Models for estimating current and future sulphur dioxide concentrations in Taipei. Bulletin of Geophysics 25 47 52 .
5. 5. Arya S. P. 1995 Modeling and parameterization of near source diffusion in weak winds. Journal of Applied Meteorology, 34 1112 1122 .
6. 6. Arya S. P. 1999 Air Pollution Meteorology and Dispersion. Oxford University Press, New York.
7. 7. Boznar M. Lesjak M. Mlakar P. 1993 A neural network-based method for short-term predictions of ambient SO2 concentrations in highly polluted industrial areas of complex terrain. Atmospheric Environment B27(2), 221-230.
8. 8. Brown M. J. Arya S. P. Snyder W. H. 1997 Plume descriptors derived from a non- Gaussian concentration model. Atmospheric Environment 31 183 189 .
9. 9. Carpenter S. B. Montgomery T. Leavitt J. M. Colbaugh W. C. Thomas F. W. 1970 Principal plume dispersion models, TVA power plants. 63rd Annual Meeting, Air Pollution Control Association, June 1970.
10. 10. Chang J. C. Hanna S. R. 2004 Air quality model performance evaluation. Meteorology and Atmospheric Physics, 87 167 196 .
11. 11. Cogliani E. 2001 Air pollution forecast in cities by an air pollution index highly correlated with metrological variables. Atmospheric Environment 35 2871 2877 .
12. 12. Comrie A. C. 1997 Comparing neural networks and regression models for ozone forecasting. Journal of air and waste management association 47 653 663 .
13. 13. Demuth C. I. 1978 A contribution to the analytical steady solution of the diffusion equation for line sources. Atmospheric Environment, 12 1255 1258 .
14. 14. EPA 1999 Air quality index Reporting Final Rule 1999 -Federal Register, Part III, CFR Part 58.
15. 15. Gardner M. W. Dorling 1999 Neural network modelling and prediction of hourly NOx and 2 concentrations in urban air in London. Atmospheric Environment 33 709 719 .
16. 16. Hanna S. R. Briggs G. A. Hosker R. P. Jr 1982 Handbook on Atmospheric Diffusion. Atmospheric Turbulence and Diffusion Laboratory National Oceanic and Atmospheric Administration.
17. 17. Hinrichsen K. 1986 Comparison of four analytical dispersion models for near-surface releases above a grass surface. Atmospheric Environment 20 29 40 .
18. 18. Huang C. H. 1979 A theory of dispersion in turbulent shear flow. Atmospheric Environment, 13 453 463 .
19. 19. Irwin J. S. Petersen W. B. Howard S. C. 2007 Probabilistic characterization of atmospheric transport and diffusion. Journal of Applied Meteorology, 46 980 993 .
20. 20. Kansal A. Khare M. Sharma C. S. 2011 Air quality modelling study to analyse the impact of the World Bank emission guidelines for thermal power plants in Delhi. Atmospheric Pollution Research, 2 (2011), 99 EOF 105 EOF .
21. 21. Katsoulis B. D. 1988 Some meteorological aspects of air pollution in Athens, Greece. Meteorology and Atmospheric Physics, 39 203 212 .
22. 22. Lin G. Y. 1982 Oxidant prediction by discriminant analysis in the South coast air basin of California. Atmospheric Environment, 16 135 143 .
23. 23. Lin J. S. Hildemann L. M. 1996 Analytical solutions of the atmospheric diffusion equation with multiple sources and height-dependent wind speed and eddy diffusivity. Atmospheric Environment, 30 239 254 .
24. 24. Mantis H. T. Repapis C. C. Zerefos C. S. Ziomas J. C. 1992 Assessment of the potential for photochemical air pollution in Athens: a comparison of emissions and air pollutant levels in Athens with those in Los Angeles. Journal of Applied Meteorology, 31 1467 1476 .
25. 25. Martin D. O. 1976 Comment on the change of concentration standard deviation with distance. Journal of Air Pollution Control Association, 26 145 146 .
26. 26. Mc Collister G. M. Wilson K. R. 1975 Linear stochastic models for forecasting daily maxima and hourly concentrations of air pollutants. Atmospheric Environment, 9 417 423 .
27. 27. Milionis A. E. Davies T. D. 1994 Regression and stochastic models for air pollution I. review, comments and suggestions. Atmospheric Environment, 28 (17), 2801 EOF -2810.
28. 28. Mooney C. J. Wilson J. D. 1993 Disagreements between gradient-diffusion and Lagrangian stochastic dispersion models, even for surface near the ground. Boundary Layer Meteorology, 64 291 296 .
29. 29. Park Y. S. Baik J. J. 2008 Analytical solution of the advection-diffusion equation for a ground-level finite area source. Atmospheric Environment, 42 9063 9069 .
30. 30. Robeson S. M. Steyn D. G. 1990 Evaluation and comparison of statistical forecast models for daily maximum ozone concentrations. Atmospheric Environment, 24B, 303 312 .
31. 31. Sanchez M. L. Pascual D. Ramos C. Perez I. 1990 Forecasting particulate pollutant concentrations in a city from meteorological variables and regional weather patterns. Atmospheric Environment, 24A (6), 1509 EOF 1519 EOF .
32. 32. Seinfeld J. H. 1986 Atmospheric Chemistry and Physics of Air Pollution. Wiley-Interscience, New York.
33. 33. Sharan M. Kumar P. 2009 An analytical model for crosswind integrated concentration released from a continuous source in a finite atmospheric boundary layer. Atmospheric Environment, 43 2268 2277 .
34. 34. Sharan M. Modani M. 2006 A two-dimensional analytical model for the dispersion of air-pollutants in the atmosphere with a capping inversion. Atmospheric Environment, 40 3479 3489 .
35. 35. Shi J. P. Harrison R. M. 1997 Regression Modelling of Hourly NOx and 2 concentrations in urban air in London. Atmospheric Environment, 31(24), 4081 EOF -4097.
36. 36. Shiva Nagendra. S. M. Venugopal K. Jones S. L. 2007 Assesment of air quality near traffic intersections in Bangalore city using air quality index. Transportation Research Part D, 12 167 176 .
37. 37. Stull R. B. 1988 An Introduction to Boundary Layer Meteorology. Kluwer Academic, 666
38. 38. Taylor G. I. 1921 Diffusion by continuous movements. Proc. London Math. Soc. Ser 2
39. 39. Tirabassi T. 2010 Mathematical air pollution models: Eulerian models, in Air pollution and turbulence modeling and application, D. Moreira and M. Vilhena, 131 155 ), CRC press, 978-1-4398-1144-3, New York.
40. 40. Turner D. B. 1969 Workbook of Atmospheric Diffusion Estimates. USEPA 999 -AP-26. U.S. Environmental Protection Agency, Washington, D.C.
41. 41. Wassermann P. D. 1989 Neural Computing Theory and Practice. New York Van Nostrand Reinhold.XX, 196 212 .
42. 42. Zannetti P. 1990 Air Pollution Modelling, Computational Mechanics Publications, Southampton, U.K. and Van Nostrand Reinhold, New York.

Written By

P. Goyal and Anikender Kumar

Submitted: 15 October 2010 Published: 05 July 2011