Open access peer-reviewed chapter

Air Pollution Dispersion Using Coupled AERMOD-WRF Modeling System and Generation of Gridded Emission Inventory of NOX over Nagpur

Written By

Rahul Boadh, Yogendra Kumar Rajoria and Vallampati Ramachandra Prasad

Submitted: 12 September 2022 Reviewed: 22 September 2022 Published: 29 October 2022

DOI: 10.5772/intechopen.108230

From the Edited Volume

Boundary Layer Flows - Modelling, Computation, and Applications of Laminar, Turbulent Incompressible and Compressible Flows

Edited by Vallampati Ramachandra Prasad, Valter Silva and João Cardoso

Chapter metrics overview

182 Chapter Downloads

View Full Metrics


This study, a work used to distinguish the dispersion patterns of air pollutants, is designs particularly NOX because of vehicular and modern sources over a quickly creating metropolitan city, Nagpur (21.15° N, 79.09° E), India, during winter and summer representative months, that is, January and April respectively in 2009. Utilizing the emission variables of industry and various vehicles, 1 km x 1 km high-resolution Gridded Emission Inventory (GEI) has produced over Nagpur city. Here, two PBL schemes form WRF model viz., local (MYJ) and non-local (YSU) schemes considered for obtaining the PBL as well as meteorological parameters. Both the schemes MYJ and YSU representing the boundary layer parameters are reasonably well and measured for generating the input parameter for AERMOD model. During the winter and summer month’s fundamentally unique scattering patterns of NOX have been observed. After carefully examining the NOX concentration from the AERMOD model, the reasonable well-by integrated YSU schemes of WRF were found to be superior to the other schemes. The current review advocates the utility of the created GEI of NOX with coupled AERMOD-WRF offline model for air quality appraisal over Nagpur.


  • gridded emission inventory
  • air quality
  • PBL schemes
  • WRF model

1. Introduction

Air pollution is a not unserious risk to the climate influencing human’s health [1, 2, 3], and horticulture and normal biological systems [4]. There are many wellsprings of air pollution, viz., usual sources, transportation, horticultural, and generation of power and industrial industries. According to World Health Organization (WHO) standards, among many Indian urban communities with populace of over a million, air contamination levels surpassed the edge. It is revealed by WHO [5] that nearly 500,000 premature deaths are in India due to only indoor air pollution. Mesoscale atmospheric flows play significant character in dispersion of air pollutants, local air quality problems, and transference [6, 7, 8, 9, 10, 11, 12]. Emission estimates and trends during 1990 to 2000 and its implications over Delhi have been reported [13]. AERMOD [14, 15, 16], a Gaussian plume air dispersion model with explicit representation of PBL processes, is widely being used for the assessment of air pollution dispersion of various pollutants and estimation of ground-level concentrations (GLC) as well. WRF-AERMOD coupling is employed for the methodology to study the dispersion patterns of PM10 over Pune [17]. Gurjar et al. [18] have evaluated emission and air quality in selected mega-cities of different countries. Mohan et al. [19] implemented AERMOD dispersion model in standalone mode for the exposure assessment of particulate matter pollution over Delhi. Goyal et al. [20] have developed grid-based mobile source emission inventory over Delhi using IVE modeling tool. It is appraised the enactment of WRF, for the situation usage as the meteorological pre-processor aimed at investigative modeling of air quality using AERMOD over Cuba [21]. A sensitivity study of different PBL schemes of parameterization schemes of WRF model coupled with AERMOD NOX dispersion over a shoreline city of Visakhapatnam [22] and over Gurugram [23] was conducted.

The good activist resolution of various parameters such as upper air meteorological and surface parameters in standalone mode is required by AERMOD model. In maximum Indian cities, this type of the PBL data is absent and over Nagpur such data are not available. The high resolution (3 km) has been utilized for integrating the meteorological variables from WRF model, which is required in air quality simulation. The result obtained is used with offline coupler [22]. The presented coupler in the present work originates the PBL and surface meteorological parameters over Nagpur and generates required input variables for the integration of AERMOD. It is obviously true that PBL processes assume a significant part in the scattering of air toxins and is crucial for utilizing reasonable PBL plot in WRF-ARW model. It has been reported that YSU followed by MYJ has shown better simulation PBL processes after validating it with available observations over Nagpur region [22, 23].

In this study, AERMOD [24, 25] and a Mesoscale weather prediction model WRF-ARW [26, 27] with YSU and MYJ schemes are employed for calculation of NOX dispersion and estimation of ground-level concentrations (GLC), due to vehicular traffic during two contrast seasons such as winter and summer during the year 2009 over Nagpur city. January has been chosen as a representative month of winter and April as summer month following the criteria of India Meteorological Department, India. In this study, GEI with high resolution (1 km × 1 km) using emission factors of different types of vehicular traffic over the study area has been developed. An additional NOX emission from a thermal power plant (TPP) located in the north direction which is 11 km away from Nagpur city has also been considered. The model-simulated NOX concentrations are validated using the monitored observations by CPCB to assess the presentation of the displaying framework with two distinct PBL plans during summer and winter.


2. Study area

Nagpur (21.15° N, 79.09° E) is the biggest city in focal of India and the second capital of the territory of Maharashtra. Nagpur is the quickly developing city and third most crowded city after Mumbai and Pune in Maharashtra. Additionally, it is the Center for advancement, industrialization, urbanization, and business action (Figure 1). The main sources of air pollution are vehicular, and industrial and domestic sources. Over Nagpur city, according to the National Ambient Air Quality Standards (NAAQS)—2008, suspended particulate matter (SPM), and respirable suspended particulate matter (RSPM) fundamentally high, the degree of nitrogen dioxide (NO2) is moderate [28]. NOX is of specific worry as it is radiated from industrial and vehicle cycles that consume non-renewable energy sources. NOX draws consideration as it responds with volatile organic compounds (VOCs) to shape ground-level ozone and it likewise responds with different constituents in air to frame nitrates and corrosive vapor sprayers. Nitrogen dioxide disturbs the nose and throat, and it seems to build vulnerability to respiratory diseases separated from nearby wellsprings of air contamination. Territorial and metropolitan air quality peculiarities are affected by a few different factors like the local meteorology.

Figure 1.

Study area (30 km × 30 km) over Nagpur [monitoring locations (ML1 = .MIDC industrial area, ML2 = MIDC office, ML3 = Govt. polytechnic college, ML4 = Institute of Engineering, ML5 = Maskasath, and ML6 = NEERI lab)].

According to climate classification by Köppen, the dry climate and tropical wet climate with dry conditions found for throughout year over Nagpur. Nagpur is situated at a height of 312.42 meters above sea level. Over Nagpur region, the winter season lies from late November to January and the temperature can goes below 10°C. Summer lies from late March to June and it is extremely hot and May being the hottest month over Nagpur region. On May 22, 2013, highest temperature (47.9°C) was recorded in Nagpur city.


3. Data and methodology

The data required for running WRF-ARMOD coupled modeling system and methodology of emission inventory and prediction of pollutant concentrations over Nagpur city using integrated WRF-ARMOD coupler are described in the following sections.

3.1 Data

In the present section, the mandatory data for offline coupled modeling system, gridded data of emission inventory from area and point sources over Nagpur, and ambient air quality monitored by CPCB are described in detail.

3.1.1 Data required for WRF-ARW model

In the current review, here accessible 6-h global final analysis (FNL) utilized with 1.0° × 1.0° grid resolution produced by the National Center for Environmental Prediction (NCEP’s) global forecast system (GFS) as initial conditions, during study period in January and April 2009.

3.1.2 Emission data

In this study, two types of emission data have been used. One is the vehicular data along with exhaust emissions considered as an area emission and the other is thermal power plant (TPP) data considered as point source over the Nagpur city. The CPCB has six monitored stations over Nagpur city. They are MIDC industrial area, MIDC office, Government Polytechnic College, Institute of Engineering, Maskasath, and NEERI Lab. The NOX-monitored data are accessible as 24 h average from CPCB, designed for all above six monitoring locations.

The data are not uniform w. r. to time at all monitoring locations; that is, the monitored data at each location are available on different days. Hence, in the current study, the comparison of observed and predicted concentrations has been done at those monitoring locations where data are available for the days, the concentrations are estimated.

3.1.3 Estimating number and types of vehicles pursuing over Nagpur

In this study, a number of vehicles plying in Nagpur city are calculated using RTO (Regional Transport Office) report for Maharashtra state [29]. Based on this report, there is an increase in the number of vehicles registered at the rate of 7.21% per year. Accordingly, we have estimated the number of vehicles at Nagpur city from the base year 2011. The yearly vehicular growth during the last 16 years is shown in Figure 2. According RTO, vehicles registered before 15 years from the current year will not be allowed to play on the road.

Figure 2.

Yearly vehicle growth for the past 16 years over Nagpur city.

We considered this as assumption that the vehicles registered in 1994 were not plying on the roads in 2009. For year 2009 (base year), diverse types of vehicles with pollutant emission factor for NOX obtained from Automotive Research Association of India (ARAI) are shown in Table 1 [30]. The TPP located at Koradi area of Nagpur is 11 km away from Nagpur city zero miles, which is considered as point source.

Vehicle typeNumber of vehiclesFractionEmission factor (gkm−1)
Motor Cycles (Petrol)226,558113,2790.15
Scooters (Petrol)170,70285,3510.14
Mopeds (Petrol)164,82782,4130.02
Personal Cars (Petrol)52,26626,1330.09
Multi Utility Vehicles (Diesel)16,89284460.67
Personal Cars (Diesel)15457730.28
Three Wheelers (Petrol)953447670.16
Three Wheelers (Diesel)357817890.51
Busses (Diesel)17728866.53
Heavy Commercial Vehicles (Diesel)867243369.3
Light Commercial Vehicles (Diesel)823841192.12

Table 1.

Different types of vehicles based on ARAI-2007 and emission factor over Nagpur city in 2009.

For the calculation of NOX emissions from the area sources, the Nagpur city (9 km × 8 km) area of 72 km2 has been considered with 1 km × 1 km grid resolution. The selected area covers all six monitoring stations of CPCB, where urban activities take place due to different sources, and predominantly vehicular exhaust is the main source of air pollution in Nagpur city. Considering this criterion, emissions from 64 grid cells were used as input into the model.

The calculation of emissions from vehicles is based on the emission factors for different types of vehicles according to ARAI 2007, vehicles explicit kind [29], the distance went by a specific vehicle, and their appropriation in light of the sort of the fuel utilized. The methodology to obtain emission from vehicular estimates was based on an earlier work [13, 19, 31];


The emission of the compound (i) is denoted by Ei, per type of vehicles numbers denoted by Vehj and different vehicle type (j) traveled distance defined by Dj in a grid box. This study has assumed based on the previous study according to ARAI, Ei;j; km the emission factor (g km−1) of different vehicles of compound (i) and vehicle type (j) per driven kilometer. The GEI of NOX over Nagpur city is shown in Figure 3. Based on the previous study, this study assumed only 50% vehicles are moving on road in per grid box. According of ARAI emission factors, vehicular emissions of NOX are designed for different types of the vehicles in each 1 km × 1 km grid box.

Figure 3.

(a) Nagpur city map (9 km × 8 km) and (b) emission inventory (gs−1 m−2) over Nagpur city.

3.2 Methodology

3.2.1 Air pollution dispersion model

AERMOD [14] is a consistent state Gaussian plume model valuable for the calculation of contamination scattering relevant to metropolitan and provincial region and several sources (volume, point and area) of emissions [14, 15, 16].

AERMOD expects the likelihood circulation work for the toxin concentration to be Gaussian in the both horizontal and vertical steady boundary layers. Then, at that point, in present instance of the convective boundary layer, while the flat scattering is Gaussian, the upward dissemination is addressed by a bi-Gaussian likelihood conveyance function.

In light of perceptions of meteorological predictions using a similarity connection, the rising profiles of the essential components, such as wind and temperature, are constructed.

AERMOD represents the vertical in-homogeneity of the PBL in its scattering estimation by averaging the boundaries of real PBL into successful boundaries of an identical homogeneous PBL. The AERMOD meteorological preprocessor, recognized as AERMET, calculates the PBL boundaries, viz., friction velocity, temperature scale, Monin-Obukhov length, surface heat flux, convective speed scale, blending tallness by involving nearby surface attributes as unpleasantness of surface, and Bowen proportion in mix with the standard meteorological perceptions [22, 32]. AERMOD is an analytic model and uses limit layer boundaries alongside routine meteorological information as contribution for assessing the contamination focuses. The overproduction of concentrations of pollutants is because of the lower reversal/blending stature under solid stable circumstances anticipated by WRF.

For predicting 24 hourly normal NOX concentration utilizing metropolitan scattering choice, the model has been utilized in this study. In this work, two different domains were chosen: Domain-1 with 30 km × 30 km with 1-km grid resolution considering the entire city including TPP and Domain-2 a finer 12 km × 12 km with 1-km grid resolution without considering TPP emissions. This exercise has been conducted so as to assess whether TPP emission does have any influence in the concentration of NOX over Nagpur. The consequences of the integrated concentration of NOX are approved against CPCB information on 24-h noticed concentration of NOX over Nagpur.

3.2.2 Weather research and forecast model

WRF version 3.6 mesoscale model established by the National Center for Atmospheric Research (NCAR) is utilized in the existing study [26, 27, 33]. This model contains the bother amounts of tension, scalars (cloud water, water vapor mixing ratio etc.), turbulent kinetic energy, surface pressure, geo-potential, three-dimensional wind, and potential temperature and includes the prognostic variables. WRF is also completely squeezable with non-hydrostatic equations. From the literature, one can say that PBL cycles and land surface affect the incorporation of winds, choppiness, and other state factors in the lower atmosphere.

The two PBL schemes YSU and MYJ both have shown better performance than the other schemes of WRF-ARW chosen based on the works over the Nagpur region and three nested domains used in study (Figure 4) [22, 25]. The complete model over view consists of model domain; different parameterizations used are presented in Table 2.

Figure 4.

WRF-ARW model domain used over Nagpur.

PBL scheme1) YSU [34]
2) MYJ [35, 36]
Vertical coordinatesTerrain-following hydrostatic pressure
vertical coordinate with 27 vertical levels
Covered area13.8°-27.4° N and 71.9°-86.2° E
Interval6 hrs
Spatial differencing scheme6th order center differencing
MicrophysicsEta microphysics [37]
Horizontal grid systemArakawa-C grid
Map projectionMercator
Surface layer parameterizationNoah land Surface Scheme
Integration time step90 sec
Time integration scheme3rd-order Runga-Kutta Scheme
Long-wave radiationRRTM scheme radiation [38]
Short-wave radiationDudhia scheme [39]
ResolutionDomain1: 27 km × 27 km
Domain2: 9 km × 9 km
Domain3: 3 km × 3 km
Grid sizeDomain1: (60× 60) × 27
Domain2: (91× 91) × 27
Domain3: (112× 112) × 27
Cumulus parameterizationKain-Fritsch scheme [40]

Table 2.

Overview of the WRF-ARW modeling system.

3.2.3 Offline coupler of WRF and AERMOD models

Hourly meteorological observations of surface and upper air remain mandatory to AERMOD. The model acknowledges a solitary station information and expects that the condition of the weather is evenly homogenous over the whole study region. In any case, the sort of meteorological information, which is fundamental to process the necessary boundary layer variables that can fill in as contribution to AERMOD, is not sweeping the study region. Subsequently, in the current work, we have utilized the offline coupler to incorporate the weather research forecast and AERMOD models. Presented coupler directly generates meteorological input files mandatory in AERMOD model [22] from the surface and PBL parameters over Nagpur from the output of WRF model.


4. Results and discussion

The current section approval of GLCs and scattering patterns of NOX acquired by utilizing the coordinated coupler of WRF-AERMOD are examined.

4.1 Simulation and dispersion of NOX

In segment, the coupled WRF-AERMOD displaying framework reenacted dispersion patterns and GLCs of NOX utilizing the different information (presented in Section 3) and Emission Inventory (EI); approval practices utilizing the accessible checked perceptions from factual investigation are introduced. The meteorological information obtained from WRF-ARW by utilizing PBL plans YSU (consequently alluded to as Model-1) and MYJ (hereafter alluded to as Model-2) are utilized in the coupled modeling framework for foreseeing NOX concentrations over the review area. The exhibition assessments of integrated NOX concentrations by Model-1 and Model-2 are learned through the accessible perceptions. The investigation and approval of the obtained GLC of NOX subsequent to coordinating the modeling framework more than two chosen domains (Domain-1 and Domain-2) of Nagpur city are presented in Section 3.2.1 with the accessible observations (given in Section 3.1.3).

4.2 Validation and GLCs of NOX

For the model combination, we utilized EI of area sources over Nagpur (without TPP in Domain-2) for integrating the concentrations of NOX pollutant. Concentrations of NOX pollutant in January and April 2009 anticipated by Model-1 and Model-2 with the accessible CPCB checked information are portrayed in Figures 5 and 6, individually.

Figure 5.

Validation of model predicted 24-hour NOX concentration (μgm−3) with CPCB observation in January 2009 over Nagpur (with the monitoring location number given in Figure 1).

Figure 6.

Validation of model predicted 24-hour NOX concentration (μgm-3) with CPCB observation in April 2009 over Nagpur (with the monitoring location number given in Figure 1).

The convergence of NOX anticipated by Model-1 shifts 52 to 110 μgm−3, though by Model-2 fluctuates between 53 and 180 μgm−3, while the noticed focus differs from 28 to 54 μgm−3 during winter. On the correlation with the observed qualities, over-forecast of re-enacted convergence of NOX by both Model-1 and Model-2 is seen (Figure 5). This might be credited to that reality the more grounded stable soundness conditions and low breezes reproduced by the WRF-ARW model prompting more catching of contamination interface of NOX.

During summer month, the anticipated fixation by Model-1 and Model-2 shifts between 20 to 85 μgm−3 and 28 to 80 μgm−3 separately. During same month, CPCB perception of NOX fixations fluctuates somewhere in the range of 18 and 61 μgm−3. On the nearby assessment of the outcomes presented in Figure 6 uncovers, the scope of fixations by the Model-1 and Model-2 is practically identical with the perceptions, with generally better execution by Model-1. The primary explanation of the over expectation of the NOX focus could be because of the exclusion of metropolitan geology in the current review since this sort of information getting is extremely intense in India. The metropolitan regions would adjust the thermo-dynamical design of the barometrical limit layer and changes the strength definition. Because of the avoidance of metropolitan elements, there are potential outcomes of over assessment of stable circumstances by the models and prompts more catching of NOX results in over forecast more. To find out the best execution of the models, factual examination is directed as given in the accompanying section.

4.3 Statistical analysis

The model presentation was assessed utilizing classical statistical [41] parameters such as geometric variance (VG), normalized root mean square error (NMSE), fractional bias (FB), mean absolute error (MAE), correlation coefficient (CC), root mean square error (RMSE), and fraction of two (FAC2) (prediction and observations). The results for the performance of statistical indicators have been shown in Table 3.

Parameter (ideal value)January 2009 (Winter)April 2009 (Summer)
YSU schemeMYJ schemeYSU schemeMYJ scheme
MAE (0)48.6177.3818.5221.85
RMSE (0)54.3487.9024.7328.92
CC (1)0.120.360.550.50
NMSE (0)0.991.940.350.44
FB (0)0.801.030.340.42
FAC2 (1)0.430.320.710.66

Table 3.

Statistical analysis of NOX concentrations during winter and summer months over Nagpur city.

It is found from these results that the RMSE and MAE of Model-1 are less as compared to Model-2 during winter month and summer month, as reported in previous studies [19, 23, 42, 43]. The statistical parameters, for example, FB and NMSE would give the necessary dependability rules in the model execution assessment. Chang and Hanna [43] examined the correlation of execution north of five field tests for three models specifically AERMOD, AMDS, and ISC-3. Hanna et al. [44] contemplated and observed and getting error VG and FAC2 separately 2.9 and 0.46 from AERMOD. This is comparable to model-1 error for both of the month as compared to model-2.

The performance of AERMOD model is examined over Delhi [19]. They observed that RMSE error is higher for each of the seven stations, and the scope of the error lies between 90.05 and 245.66. In the current work, RMSE error is sensibly well in contrast to the past review and this is given affirmation that model-1 is given better performance in contrast to the model-2.

Less values of FB and NMSE indicate that the overall deviation is less by Model-1 compared to Model-2. These boundaries uncover that exhibition of Model-1 is superior to the Model-2 during winter and summer months. FB demonstrates how well the calculation creates the normal qualities around the normal upsides of observed factors. The ideal value is zero of FB, but it can range from −2 to 2. An absolute value of 0.67 corresponds to a prediction inside an influence of two of observation [45]. Interestingly, correlation coefficient of Model-2 is to found better as compared to Model-1 for during winter and summer months. Yet, in light of the in general factual assessment, the presentation of expectation of NOX concentrations utilizing the meteorological factors integrated from Model-1 is superior to Model-2. Large values of MAE, RMSE, NMSE, FB, and VG1 indicate overprediction of concentrations by both models due to low wind conditions.

4.4 Spatial dispersion of GLC NOX concentrations

In the current segment, spatial scattering examples of GLC NOX concentrations are obtained from WRF-AERMOD modeling system utilizing two definition plans; specifically, YSU and MYJ are introduced over Domain-1 and Domain-2 over Nagpur city (as presented in Section 3.2.1). The coupled demonstrating framework has been incorporated for the entire long periods of winter and summer months over the picked spaces at whatever point the CPCB monitoring of GLC perceptions is accessible during winter- and summer-represented months. In this chapter, the aftereffects of 1 day each from these 2 months, January 8 and April 21, 2009 are just introduced.

The dispersion of 24-hourly averaged NOX concentrations for 30 km × 30 km domain on January 8 by utilizing PBL schemes of YSU and MYJ is represented in Figure 7a–e along with wind roses. On January 8, 2009, the most extreme GLCs are anticipated by MYJ plot (Figure 7b) contrasted with the YSU (Figure 7a), in light of the fact that MYJ (Figure 7d) re-enacted low breezes than that of YSU as displayed in Figure 7c. The extent of winds obtained from MYJ is lesser as compared with YSU.

Figure 7.

Spatial distribution of 24-hour predicted NOX GLCs (μgm-3) (a) YSU, (b) MYJ (c) wind roses of YSU (d) wind roses of MYJ and (e) OBS wind roses on January 8, 2009, over Nagpur by using 30 km ×30 km domain.

This clearly shows that the winds stay improved arrested by YSU upon an evaluation by the observations (Figure 7e) with MYJ. Maximum concentrations of NOX by using both schemes (YSU and MYJ) inside the city stay only due to vehicular emission. From the concentration patterns, it is clearly seen that nearby is nope influence of TPP-generated NOX over Nagpur as seen in Figure 7a and b.

The NOX concentrations of April 21, 2009 are depicted in Figure 8ae along with wind roses. The similar patterns of the wind roses are shown by both schemes YSU and MYJ, respectively (as seen in Figure 8c and d) as shown in observations (Figure 8e). The wind simulations of MYJ scheme are showing the low winds in comparison with YSU scheme but are in better arrangement with the observation. The GLCs of NOX predicted by both schemes YSU and MYJ are presented in Figure 8a and b, respectively. MYJ predicted higher concentration compared to YSU, because MYJ simulated low wind compared with YSU. This is on the grounds that the impact of wind speed is contrarily corresponding to the centralizations of the pollutants [46, 47].

Figure 8.

Spatial distribution of 24-hour predicted NOX GLCs (μgm-3) (a) YSU, (b) MYJ (c) wind roses of YSU (d) wind roses of MYJ, and (e) OBS wind roses on April 21, 2009 over Nagpur by using 30 km ×30 km domain.

In comparison to the winter day, there is significantly less NOX emission from TPP over Nagpur during this summer day. Comparative sort of results is seen during every one of the times of the both study months. Henceforth, we can infer that the impact of the power plant emissions of NOX concentrations is negligible over Nagpur city for both winter and summer months. Based on this result, a fine resolution of 12 km × 12 km domain was selected for predicting the GLCs of NOX over the Nagpur city for these two contrasting months and results during January 8, 2009 and April 21, 2009 are presented.

The forecasted concentration on January 8, 2009 from both (YSU and MYJ) schemes along with wind roses is presented in Figure 9. The determined NOX concentration forecasted by PBL schemes YSU and MYJ (Figure 9a and b respectively) on January 8 is mostly in western part of the city, because the predominant wind directions are east northeasterly and southeasterly (Figure 9c and d). The predicted concentration of NOX by both (YSU scheme and MYJ scheme) the schemes on April 21, 2009 is presented in Figure 10. The predominant wind direction is northwesterly (Figure 10c and d) and the results (the maximum concentration of the NOX) predicted by YSU (Figure 10a) and MYJ (Figure 10b) schemes are in southeastern quadrant of the Nagpur city.

Figure 9.

Spatial distribution of 24-hour predicted NOX GLCs (μgm−3) (a) YSU, (b) MYJ (c) wind roses of YSU (d) wind roses of MYJ and (e) OBS wind roses on January 8, 2009 over Nagpur by using 12 km × 12 km domain.

Figure 10.

Spatial distribution of 24-hour predicted NOX GLCs (μgm−3) (a) YSU, (b) MYJ (c) wind roses of YSU (d) wind roses of MYJ and (e) OBS wind roses on April 21, 2009 over Nagpur by using 12 km ×12 km domain.

Similar kind of signatures of maximum concentrations is noticed during study period of these two contrasting months. These varieties of scattering of concentrations of NOX are critical during these two differentiating months’ outcomes in particular unique air quality over Nagpur. In light of the model inferred most extreme NOX concentrations, not many spots of the Nagpur city are distinguished as area of interest areas meaning decayed air quality and are displayed in Table 4.

Hotspot locations over Nagpur City based on higher NOX concentrationsAverage NOX Concentration (μgm−3) of January 2009 (Winter)Average NOX Concentration (μgm−3) of April 2009 (Summer)
Sita Buldi921054144
Central Bazar Road80903825
North Ambazari Road65873228
Anjuman Square78963024
Liberty Square62842926
Sankar Nagar Square78983428
Low College Square75843723
Katol Road65752822
Police Com. Office73963526
Medical Chowk48523436
Baidwanath Chowk35482934
Great Nag Road54623335
NMC Office Mahal62673829
Sakkardara Chowk38433031
Rajendra Nagar29352129

Table 4.

Identified deteriorated air quality (hotspot) locations over Nagpur City.

The recognized hotspot sites are Sita Buldi, Central bazaar road, North Ambazari road, Anjumna square, Liberty square, Low College square, Katol road, police commissioner office, and Sankar Niger square located in the western part to Nagpur city. These locations have high concentrations in January month and low concentrations in April month. Other locations such as medical chowk, Baidwanath chowk, Great nag road, NMC office Mahal, Sakkardara road, and Rajendra nagar are in north and southeastern part of the Nagpur city. These locations have higher concentrations of NOX in January and low concentrations in April 2009. These changes in concentrations in the different locations of the Nagpur city during these two representative months of summer and winter season have discovered the monthly variability of deteriorated air quality locations in relations to extreme concentrations of NOX.


5. Conclusions

WRF-ARW model is combined with AERMOD for processing GLCs of NOX over Nagpur city. Scattering of NOX by utilizing YSU and MYJ planetary limit layer plans of WRF-ARW model over Nagpur has been assessed.

Accurate assessment of the air pollution dispersion as well as air quality issues requires high-resolution emission inventory and air pollution dispersion models with better representation of the PBL parameters. One of the contributions emanated from this study is the preparation of gridded emission inventory using the vehicular traffic data over Nagpur city with 1-km resolution. The offline coupled modeling system of WRF-AERMOD reasonably simulated GLCs concentrations of NOX during both representative months of winter and summer seasons over Nagpur city.

The coarse resolution of the (30 km × 30 km) has shown negligible impact of NOX emission released from Thermal Power Plant in both the contrasting. The GLCs obtained from the finer resolution (12 km × 12 km) over of the Nagpur city has revealed significant dispersion patterns during the 2 months. Deteriorated air quality locations referred to as hotspot locations, based on the maximum NOX concentrations, are identified for Nagpur. Significant variations in the concentrations are noticed during these two seasons. This study advocates the better representation of air pollution dispersion has shown by offline coupled modeling system of WRF-AERMOD and given confidence for its utility in the air quality assessment.


  1. 1. Dockery DW, Pope CA, Xu X, Spengler JD, Ware JH, Fay ME, et al. An association between air pollution and mortality in six U.S. cities. The New England Journal of Medicine. 1993;329(24):1753-1759
  2. 2. Solomon S, editor. Intergovernmental Panel on Climate Change, Summary for policymakers, in Climate Change. In: The Physical Science Basis. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change. New York: Cambridge University Press; 2007. pp. 1-18
  3. 3. Ramanathan V, Feng Y. Air pollution greenhouse gases and climate change: Global and regional perspective. Atmospheric Environment. 2009;43:37-50
  4. 4. Grimm NB, Faeth SH, Golubiewski NE, Redman CL, Wu J, Bai X, et al. Global change and the ecology of cities. Science. 2008;319(5864):756-760
  5. 5. World health organization. Quantifying Selected Major Risks to Health, Chapter 4 from the World Health Report, World Health Organization, Geneva. 2002
  6. 6. Srinivas CV, Venkatesan R, Baskaran R, Rajagopal V, Venkatraman B. Regional scale atmospheric dispersion simulation of accidental release of radionuclides from Fukushima Dai-ichi reactor. Atmospheric Environment. 2012;61:66-84
  7. 7. Angevine WM, Eddington L, Durkee K, Fairall C, Bianco L, Brioude J. Meteorological model evaluation for CalNex 2010. Monthly Weather Review. 2012;140:3885-3906
  8. 8. Fast JD, Gustafson WI, Berg LK, Shaw WJ, Pekour M, Shrivastava M, et al. Transport and mixing patterns over Central California during the carbonaceous aerosol and radiative effects study (CARES). Atmospheric Chemistry and Physics. 2012;12:1759-1783
  9. 9. Baker KR, Misenis C, Obland MD, Ferrare RA, Scarino AJ, Kelly JT. Evaluation of surface and upper air fine scale ARW meteorological modeling of the may and June 2010 CalNex period in California. Atmospheric Environment. 2013;80:299-309
  10. 10. Hu XM, Nielsen GJW, Zhang F. Evaluation of three planetary boundary layer schemes in the WRF model. Journal of Applied Meteorology and Climatology. 2010;49:1831-1844
  11. 11. Hu XM, Klein PM, Xue M. Evaluation of the updated YSU planetary boundary layer scheme with in WRF for wind resource and air quality assessments. Journal of Geophysical Research – Atmospheres. 2013;118:10490-10505
  12. 12. Srikanth M, Satyanarayana ANV, Srinivas CV, Kumar M. Mesoscale atmospheric flow-field simulations for air quality modelling over complex terrain region of Ranchi in eastern India using WRF. Atmospheric Environment. 2015;107:315-328
  13. 13. Gurjar BR, Aardenne JA, Lelieveld J, Mohan M. Emission estimates and trends (1990–2000) for mega city Delhi and implications. Atmospheric Environment. 2004;38:5663-5681
  14. 14. Cimorelli AJ, Perry SG, Venkatram A, Weil JC, Paine RJ, Wilson RB, Lee RF, Peters WD, Brode RW, Paumier JO. AERMOD: Description of model formulation. US Environmental Protection Agency, EPA Report No. 454/R-03-002d, 2004; p 85.
  15. 15. Cimorelli AJ, Perry GS, Venkatram A, Weil JC, Paine RJ, Wilson RB, et al. AERMOD: A dispersion model for industrial source applications. Part 1: General model formulation and boundary layer characterization. Journal of Applied Meteorology. 2005;44:682-693
  16. 16. Perry SG, Cimorelli AJ, Paine RJ, Brode RW, Weil JC, Venkatram A, et al. AERMOD: A dispersion model for industrial source applications. Part II: Model performance against 17 field study databases. Journal of Applied Meteorology. 2005;44:694-708
  17. 17. Kesarkar AP, Dalvi M, Kaginalkar A, Ojha A. Coupling of the weather research and forecasting model with AERMOD for pollutant dispersion modeling, a case study for PM10 dispersion over Pune, India. Atmospheric Environment. 2007;41:1976-1988
  18. 18. Gurjar BR, Butler TM, Lawrence MG, Lelieveld J. Evaluation of emissions and air quality in megacities. Atmospheric Environment. 2008;42:1593-1606
  19. 19. Mohan M, Bhati S, Rao A. Application of air dispersion modeling for exposure assessment from particulate matter pollution in mega City Delhi. Asia-Pacific Journal of Chemical Engineering. 2011;6:85-94
  20. 20. Goyal P, Mishra D, Kumar A. Vehicular emission inventory of criteria pollutants in Delhi. Springer Plus. 2013;2:216
  21. 21. Carbonell LT, Mastrapa GC, Rodriguez YF, Escudero LA, Gacita MS, Morlot AB, et al. Assessment of the weather research and forecasting model implementation in Cuba addressed to diagnostic air quality modeling. Atmospheric Pollution Research. 2013;4:64-74
  22. 22. Boadh R, Satyanarayana A NV, Rama Krishna TVBPS, Madala S. Sensitivity of PBL parameterization schemes of Weather Research Forecasting Model and coupling with AERMOD in the dispersion of NOX over Visakhapatnam (India). Asia-Pacific Journal of Chemical Engineering. 2015. DOI: 10.1002/apj.1876
  23. 23. Sharma N, Singh R, Rojoria YK, Rajendra P, Boadh R. Coupled WRF-AERMOD modeling system by using dispersion of air pollutant and generation of gridded emission inventory of NOX over Faridabad and Gurugram. International Journal of Scientific Research in Mathematical and Statistical Sciences. 2021;8(1):01-12
  24. 24. Boadh R, Satyanarayana ANV, Rama Krishna TVBPS, Madala S. Sensitivity of PBL schemes of WRF-ARW model in simulating boundary layer flow parameters for its application to air pollution dispersion modeling over a Tropical Station. Atmosfera. 2016;29(1):61-81
  25. 25. Boadh R, Satyanarayana ANV, Rama Krishna TVBPS, Subba Rao A, Rajendra P. Numerical simulation of boundary layer flow parameters by using wrf-arw model over a tropical region. Journal of Industrial Pollution Control. 2017;33(1):1148-1154
  26. 26. Skamarock WC. Positive-Definite and Monotonic Limiters for Unrestricted-Time-Step Transport Schemes. Monthly Weather Review. 2006;134: 2241-2250
  27. 27. Skamarock WC, Weisman ML. The impact of positive-definite moisture transport on NWP precipitation forecasts. Monthly Weather Review. 2009;137:488-494
  28. 28. CPCB (Central Pollution Control Board), India. National Ambient Air Quality Standards (NAAQS). Gazette notification, New Delhi. 2009.
  29. 29. MVS. Motor transport statistics of Maharashtra. 2011.
  30. 30. ARAI. Draft report on “Emission Factor development for Indian Vehicles” as a part of Ambient Air Quality Monitoring and Emission Source Apportionment Studies. 2007. p. 1-94
  31. 31. Kandlikar M, Ramachandran G. The causes and consequences of particulate air pollution in urban India: A synthesis of the science. Annual Review of Energy and the Environment. 2000;25:629-684
  32. 32. Boadh R, Satyanarayana ANV, Rama Krishna TVBPS. Assessment of dispersion of oxide of nitrogen using AERMOD model over a tropical industrial region. International Journal of Computer Applications. 2014;90(11):43-50
  33. 33. NCEP-DSS1. NCEP Global Tropospheric Analyses, (1 × 1 deg), from Data Support Division (DSS). NCAR. 2005; Dataset # 083.2. / [Accessed June 2005]
  34. 34. Hong SY, Noh Y, Dudhia J. A new vertical diffusion package with explicit treatment of entrainment processes. Monthly Weather Review. 2006;134:2318-2341
  35. 35. Janjic ZA. The step-mountain eta coordinate model, further developments of the convection, viscous sub layer and turbulence closure schemes. Monthly Weather Review. 1994;122:927-945
  36. 36. Janjic ZI. Non-singular Implementation of the Mellor–Yamada Level 2.5 Scheme in the NCEP. Meso Model. 2002;437:61
  37. 37. Ferrier BS, Lin Y, Black T, Rogers E, DiMego G. Implementation of a new grid-scale cloud and precipitation scheme in the NCEP Eta model. American Meteorological Society: Conference Proceeding; 2002. pp. 280-283
  38. 38. Mlawer EJ, Taubman SJ, Brown PD, Iacono MJ, Clough SA. Radiative transfer for inhomogeneous atmosphere RRTM, a validated correlated-k model for the longwave. Journal of Geophysical Research. 1997;102(D14):16663-16682
  39. 39. Dudhia J. Numerical study of convection observed during the winter monsoon experiment using a mesoscale two-dimensional model. Journal of the Atmospheric Sciences. 1989;46:3077-3107
  40. 40. Kain JS. The Kain–Fritsch convective parameterization an update. Journal of Applied Meteorology. 2004;43(1):170-181
  41. 41. Wilks DS. Statistical Methods in the Atmospheric Science. Second ed. U.S.A: Elsevier Academic Press; 2008. p. 627
  42. 42. Kumar A, Luo J, Bennett GF. Statistical evaluation of lower flammability distance (LFD) using four hazardous release models. Process Safety Progress. 1993;12(1):1-11
  43. 43. Chang JC, Hanna SR. Air quality model performance evaluation. Meteorology and Atmospheric Physics. 2004;87:167-196
  44. 44. Hanna S, Yang R. Evaluations of mesoscale models’ simulations of near-surface winds, temperature gradients, and mixing depths. Journal of Applied Meteorology. 2001;40:1095-1104
  45. 45. Arya SPS. Air Pollution Meteorology and Dispersion. New Yrok: Oxford University Press; 1999. p. 310
  46. 46. Harrison RM. An Introduction to Pollution Science. RSC Publication; 2006. p. 29. DOI: 10.1039/9781847555410
  47. 47. Rao MN, Rao HVN. Air pollution. 1st ed, New Delhi: Tata McGraw-Hill; 1990. p. 15

Written By

Rahul Boadh, Yogendra Kumar Rajoria and Vallampati Ramachandra Prasad

Submitted: 12 September 2022 Reviewed: 22 September 2022 Published: 29 October 2022