Next Article in Journal
A Rigid Cuckoo Search Algorithm for Solving Short-Term Hydrothermal Scheduling Problem
Next Article in Special Issue
Considering Condensable Particulate Matter Emissions Improves the Accuracy of Air Quality Modeling for Environmental Impact Assessment
Previous Article in Journal
Adsorption/Desorption Capability of Potassium-Type Zeolite Prepared from Coal Fly Ash for Removing of Hg2+
Previous Article in Special Issue
Can Carbon Trading Policies Promote Regional Green Innovation Efficiency? Empirical Data from Pilot Regions in China
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Sustainability
Volume 13
Issue 8
10.3390/su13084276
sustainability-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Use of Simulated and Observed Meteorology for Air Quality Modeling and Source Ranking for an Industrial Region

by
Awkash Kumar
1,2,*,
Anil Kumar Dikshit
1 and
Rashmi S. Patil
1
1
Environmental Science and Engineering Department, Indian Institute of Technology, Bombay, Mumbai 400 076, Maharashtra, India
2
SAGE, Sustainable Approach for Green Environment, Powai, Mumbai 400 076, Maharashtra, India
*
Author to whom correspondence should be addressed.
Sustainability 2021, 13(8), 4276; https://doi.org/10.3390/su13084276
Submission received: 3 February 2021 / Revised: 2 April 2021 / Accepted: 9 April 2021 / Published: 12 April 2021
(This article belongs to the Special Issue Air Quality Characterisation and Modelling)
Browse Figures
Versions Notes

Abstract

:
The Gaussian-based dispersion model American Meteorological Society/Environmental Protection Agency Regulatory Model (AERMOD) is being used to predict concentration for air quality management in several countries. A study was conducted for an industrial area, Chembur of Mumbai city in India, to assess the agreement of observed surface meteorology and weather research and forecasting (WRF) output through AERMOD with ground-level NOx and PM10 concentrations. The model was run with both meteorology and emission inventory. When results were compared, it was observed that the air quality predictions were better with the use of WRF output data for a model run than with the observed meteorological data. This study showed that the onsite meteorological data can be generated by WRF which saves resources and time, and it could be a good option in low-middle income countries (LIMC) where meteorological stations are not available. Also, this study quantifies the source contribution in the ambient air quality for the region. NOx and PM10 emission loads were always observed to be high from the industries but NOx concentration was high from vehicular sources and PM10 concentration was high from industrial sources in ambient concentration. This methodology can help the regulatory authorities to develop control strategies for air quality management in LIMC.

1. Introduction

Urbanization-related issues have become very prominent across the world [1,2,3], especially in developing countries like India where cities have started facing an acute air pollution problem due to urbanization [4]. Many Indian megacities, such as Delhi, Mumbai, Bangalore, and Kolkata, are witnessing increasing health problems due to rapid increase in air pollution [5]. The total health cost due to air pollution for Mumbai was about USD 8 billion for the year 2012 [6]. This problem becomes particularly complex to resolve in urban areas because of diverse emission sources such as vehicles, industries, bakeries, hotels, diesel generating sets and combustion of solid fuels in the domestic sector.
Air quality monitoring networks have been installed at various locations in many cities. Also, installation and operation of a large number of air quality monitoring stations need considerable financial resources from government which may not be supported in low-middle income countries (LIMC). This monitoring data is increasingly used to communicate the existing status of air quality. However, it doesn’t contribute to the understanding of sources and meteorological factors. Although the observed data represent air quality status for a particular location only, the use of dispersion models can provide information about much larger areas. Further, modeling helps in the determination of concentration plots on spatial and temporal scales and contributions from different types of source for air pollutants [7,8,9,10,11]. The dispersion model can also be used to identify pollution sources with the help of emission inventory [12,13]. This is very useful in making rational management strategies [7,14,15,16,17,18,19,20]. A dispersion model can also determine the contribution of various sources in a region whereas a receptor model determines source contribution at a particular location [21]. Data required for dispersion models include emission inventory, geographical data, and meteorological data of the region [22]. Data availability, especially meteorological data, is an important factor for the assessment of air quality in LMIC because running a meteorological monitoring station requires resources. The use of poor-quality meteorological data in air quality models may contain significant adverse effect model output quality [23,24]. Meteorological data is generally taken from a nearby meteorological station and is used for the study region. The results of air quality model may have significant error despite advanced computer technology, and various techniques like numerical modeling techniques, performance evaluations of state-of-the-art [23,24,25,26,27]. A survey has been done for the air quality and meteorological model [28].
The observed meteorological data from a monitoring station may not give good performance in air quality modeling for the urban industrial region where several emission sources are present at multiple heights and variation in topography. An alternative is to generate onsite meteorological data using a meteorological model which could be an effective option in LMIC. A study was conducted on the coupling of American Meteorological Society/ Environmental Protection Agency Regulatory Model (AERMOD) with Weather Research and Forecasting (WRF) model in Pune city (India) for a single pollutant PM10 [29]. However, their predicted concentration obtained by the WRF-generated meteorology and observed values have not been compared and contribution of various sources in the study region has not been estimated. Short term air quality forecasting also has been carried out using WRF forecasted meteorology and AERMOD for five days for Chembur region [30]. In these studies, the requirement of horizontal homogeneous hourly surface and upper meteorological data has been fulfilled from WRF model for AERMOD. The main objective of this study was to generate onsite meteorological data at mesoscale using WRF model and compare the results with observed meteorological data. Then, we proposed to use both the data in air quality modeling and to evaluate the option of making WRF coarse resolution output feasible in LMIC. This study was also continued to rank the contribution of emissions and ambient concentrations from sources for NOx and PM10. This will be useful for air quality management of the urban area for regulation purposes [31] in LMIC.

2. Study Area

The study area, Chembur, represents an industrial site of Mumbai city in India with global coordinates 19.05° N and 72.89° E. This area covers M East and M West wards of Municipal Corporation of Greater Mumbai, which is one of the financial centers of India as shown in Figure 1a. Chembur has a population of 1.2 million. It measures 6.5 km east-to-west and 8.45 km north-to-south, as shown in Figure 1c. This region has marine alluvium type of soil and North-South running basalt hills to its South [32]. The topographical features have been shown in Figure 1d in the Universal Transverse Mercator coordinate system. The elevation is maximum at the central part of the study area and minimum along the boundary of the Eastern study area. The elevation ranges from 1 to 200 m. The elevation just above the location of Rashtriya Chemicals and Fertilizers Limited (RCFL) is 100 to 200 m. Major industries in this area are Bharat Petrochemical Corporation Limited (BPCL), Hindustan Petroleum Refinery Corporation Limited (HPCL), Tata Thermal Power Corporation Limited (TPCL) and RCFL. Containers and heavy-duty vehicles from this area use the Port Trust Road, Mahul Road and Ramakrishna Chemburkar Marg (R C Marg). Road conditions are poor due to the continuous movement of heavy vehicles. The residential areas spread over the north boundary of the study area has a residential zone comprising Chheda Nagar (between point 2 and 3) and Shramjivi Nagar (Left side area of point 2), the south boundary is adjacent to the Tata Thermal Power Plant. The west boundary lies by RCFL and Mahul, and the east boundary is aligned with Shahyadri Nagar and Prayag Nagar.
Around twenty years ago, Chembur was one of the most polluted regions in Mumbai. With the sustained effort and pressure from authorities and industries for implementing a series of control measures, the region has witnessed an improvement in air quality. In the last two decades, the region characteristics have improved due to the closure of many industries, but residential development and vehicular density has increased [33]. Chembur still needs appropriate air quality studies for developing management strategies, as its ambient air quality is poor when compared with the National Ambient Air Quality Standard 2009, Central Pollution Control Board (CPCB) New Delhi (India). CPCB has published a document, giving a Comprehensive Environmental Pollution Index (CEPI) score for various industrial regions in the country. Chembur has a score of 69.19 CEPI in this report [34]. This score shows that this region should be rated as a severely polluted area. Hence, the region requires better understanding of air quality processes so that effectiveness of the action plan can be realized.

3. Methodology and Data

The schematic data flow of the study has been shown in Figure 2. AERMOD requires emission inventory and nine hourly meteorological parameters (wind speed, wind direction, rain fall, temperature, humidity, pressure, ceiling height, global horizontal radiation and cloud cover) as the input data. These meteorological parameters were generated from the WRF model for the year. Also, the meteorological parameters were observed at RCFL for the same time period and both data sets were compared. Prediction of concentrations using an air quality model (AERMOD) was carried out with the observed meteorological data and WRF generated data. Meteorological parameters were prepared in columns and temporal resolution was prepared in rows of a spreadsheet. This spreadsheet was processed in AERMET which is a pre-processor of AERMOD. The terrain data at 90 m resolution of Shuttle Radar Topography Mission (SRTM) was used in AERMAP which is also the pre-processor of AERMOD. Then, the AERMOD model was used to predict concentration of NOx and PM10 as shown in Figure 2. Also, comparisons of both the models, WRF and AERMOD, were done. The metrological model, the setup of parameterization of variables, dispersion model, emission load, and observations have been presented section-wise.

4. Meteorological Model

The mesoscale model, Advanced Research WRF model version 3.2, has been used in this study [35]. This model is designed to assist both atmospheric research and operational forecasting needs [36]. NCEP FNL (Final) Operational Global Analysis data have been used as an input for WRF, which are on 1.0 ×1.0 degree grids prepared operationally every six hours. This product is from the Global Data Assimilation System (GDAS), which continuously collects observational data from the Global Telecommunications System (GTS) and other sources for analyses [37]. It is a limited area, non-hydrostatic primitive equation model with multiple options for various physical parameterization schemes. This version employs Arakawa C-grid staggering for the horizontal grid and a fully compressible system of equations. A terrain following hydrostatic pressure coordinates with vertical grid stretching is implemented vertically. The time split integration uses a third order Runge Kutta scheme with smaller time step for acoustic and gravity wave modes. The WRF model physical options used in this study consist of the WRF model Single Moment 6-class simple ice scheme for microphysics, the Kain–Fritsch scheme for the cumulus convection parameterization, and the Yonsei University planetary boundary layer scheme. The rapid radiative transfer model and the Dudhia scheme are used for longwave and shortwave radiation, respectively, while the Noah land surface model has been selected. All these parameterizations constitute a well-tested suite of schemes over the Indian region [38,39,40]. The model domain extends between 71° E to 81° E zonally and 11° N to 21° N meridionally, consisting of 100 by 100 grid points with 25 km grid spacing as shown in Figure 1b. The model was run from 1st January to 31st December of the year. The model has 28 vertical levels with the top of the model at 10 hPa. Topography as well as snow cover information have been obtained from the United States Geological Survey. The meteorological parameters have been extracted from the WRF model at ground level. The WRF model has been run at 25 km resolution which provides time series meteorological parameters for a specific period at a particular location. In this study, 9 hourly meteorological parameters (cloud cover, temperature, pressure, relative humidity, wind direction, wind speed, ceiling height, rainfall, and global horizontal radiation) have been simulated for the year using WRF. WRF gives output in network common data format, and GRADS v 2.2 is post processing software which reads the network common data format. The output from WRF was fed to Grads 2.2 to generate digital hourly meteorological data to arrange in an Excel spreadsheet. The AERMET required data in excel spreadsheet or other formats. The input in AERMET was given in excel spreadsheet which was prepared using the data obtained by Grads 2.2. Here, hourly data for each meteorological parameter are provided in different columns. The spreadsheet meteorological data were imported in AERMET which is pre-processor of AERMOD.

5. Dispersion Model

Dispersion model uses emission inventory, geographical and meteorological data to predict concentration at the receptor’s point in the study region. The format of the input data varies with different models. There are many specific models for vehicular and industrial sources as well as for a variety of sources [19,41,42,43]. Industrial Source Complex (ISC3), developed by the United State Environmental Protection Agency, is a steady-state Gaussian plume model which can be used to assess pollutant concentrations from a wide variety of sources associated with an industrial source complex [44]. ISC3 operates in both long-term (ISCLT3) and short-term (ISCST3) models. ISCST3 model is the regulatory model in India and it has been used in many case studies [45]. Later on, it was updated to AERMOD whose performance was appreciable as they added some advanced algorithms to get more accurate results [46]. The air quality model that we use, AERMOD, has been applied to evaluate dispersion of several pollutants, including PM10, HCN, SO2, SF6, and VOCs and is recommended widely by regulatory authorities [47,48,49,50,51].
The study area (as given Figure 1c) was given in AERMOD and emission locations were digitized according to real earth surface reference and quantities of emissions were put based on estimated emission inventory. Therefore, there is no resolution concept for emission inventory in this study. The meteorological data output from the WRF model was processed in AERMET and its output was fed into AERMOD. The pre-procedure AERMAP of AERMOD calculates representative terrain-influence height, also referred to as the terrain height scale, at a receptor in modeling of air quality. Cartesian uniform gridded receptors were given, apart from discrete receptors and all receptors were at 2 m height. Anemometer height was 10 m and surface roughness length was 1 m in this model run. Building downwash terminology was not considered. AERMOD calculates concentration for each hour using hourly meteorological data for each pollutant separately.

6. Emission and Concentration Data

Emission load has been computed for point sources (specifically 36 stacks of BPCL, 30 stacks of HPCL, 4 stacks of Tata Power and 17 stacks of RCFL), line sources (the 6 roads of Chembur), and area sources (e.g., bakeries, hotels and restaurants, crematorium and domestic sector). These area sources were taken from a previous study (CSIR-NEERI) [52] for M East and M West wards, where area sources emission load has been computed. These sources are for the region of Chembur (M East and West Ward) where domestic sectors are available. Industrial emission data are collected from industries and vehicular emission inventory are prepared based on field survey data for the study period. Vehicular emission rate is estimated using the actual number of vehicles in unit time, emission factors and vehicle kilometer travelled [53]. The percentage of contribution is estimated from type of the sources after making emission inventory of the region. PM10 and NOx emissions are mainly caused by industries (94% and 64%, respectively). Remaining 36% of NOx emission is contributed by vehicles (18%), domestic sources (17%) and others (1%). The emission factors of vehicles are available for particulate matter, and this has also been taken as PM10. Figure 3 depicts emission loads of NOx and PM10 from various types of sources in the study area. The observed concentration data were collected for NOx and PM10 at industrial sites, i.e., HPCL and BPCL at a height of about 3.5–4.0 m. Continuous (hourly) ambient air quality monitoring was done at these sites using Telydene instrument for NOx and PM10.

7. Results and Discussion

The results of WRF were used as inputs in estimation of concentration by AERMOD. NCEP FNL (Final) Operational Global Analysis data was used to process as an input for WRF. It produced 30 meteorological parameters for the required time resolution and study period. AERMET (pre-processor of AERMOD) required nine meteorological parameters and these were extracted, out of 30 meteorological parameters for air quality modeling. After this, the contributions from various sources in ambient concentration were estimated in this study.

8. Validation of Wind and Temperature Time Series

WRF generated the nine meteorological parameters hourly which were used in air quality modeling. The validation of meteorological output from WRF was done. In the validation of WRF, output time series temperature and wind are significant for air pollution. Therefore, the point of validation of temperature was conducted using the hourly temperature of the year of WRF with observed temperature data of RCFL industry, Chembur. A fair estimate of the dispersion of pollutants in the atmosphere is possible based on the frequency distribution of wind direction as well as wind speed. Wind transports pollutants from various sources, causing turbulent mixing and diluting pollution. Boundary layer cumulus clouds vent pollution into the free troposphere, and temperature and humidity levels in the boundary layer affect chemical reactions and the rates at which many dangerous compounds are formed.

8.1. Validation of Temperature

Figure 4 shows a comparison of surface temperature of the study area, between WRF simulation and observations, and it can be seen that they are in good agreement. The average percentage error between simulated and observed temperature is −5.3%. The standard deviation of percentage error is 10.4.

8.2. Validation of Wind

Wind data was derived using WRF model at the height of 10 m from the ground for a whole year, and wind rose of Chembur, which was simulated by (a) WRF model and (b) RCFL observed data for the year, was plotted (Figure 5). Maximum wind persistence corresponds from the west and south-west direction in WRF, while in RCFL observed data, the wind was found to blow in the west and north-west direction for most of the time. Observed wind data had more periods of calmness than the data simulated by WRF. Observed wind rose had 62% calm condition whereas simulated wind rose had 4% calm condition over a period. The root mean square error and mean bias error between predicted and observed wind speed were 4.05 and 3.29, respectively. The root mean square error and mean bias error between predicted and observed temperature were 3.69 and −1.5, respectively. The observed data of wind was collected at RCFL industry at the height of 6 m, and it represents the microscale domain. Chembur has considerable variation in topography as shown Figure 1d. The topography changes after a few meters and this causes the wind to divert. Consequently, this variation of topography may cause the mismatch of wind rose of WRF and observed data.
Modeling was performed with both the meteorological parameters of WRF output and RCFL observed data. WRF performance of wind is poor due to the surface inhomogeneity. However, for the purpose of dispersion modeling, we consider WRF as a good representation of the mesoscale flow. The observed wind at RCFL is not representative of the entire Chembur region for dispersion modeling. This may be because the maximum emission is from industrial sources in this region, which are at an elevated height. Hence, for these sources, the mesoscale meteorology generated by the WRF model may be more appropriate than the observed microscale meteorology.

9. Validation of NOx and PM10

The annual PM10 and NOx concentration contour plots for all sources of the study domain are shown in Figure 6a,b respectively. A comparison of simulated concentration using WRF output with the observed concentration of PM10 and NOx are in Table 1, and the model was well-compared for this study area. The root mean square error and mean bias error between predicted and observed concentrations for NOx were 1.76 and 0.063, respectively, while they were 0.41 and 0.83 for PM10, respectively. The standard deviations for NOx at BPCL and HPCL were 33.6 and 30.2, respectively, and the standard deviations for PM10 at BPCL and HPCL were 16.4 and 12.4, respectively. It was seen that the values obtained through air quality modeling were closer to the observed concentrations with the mesoscale meteorology than the surface level meteorology. The model results using observed surface meteorology were high. Modeled values were in good agreement with the observed values at both locations for NOx, but for PM10, simulated concentrations were lesser than the observed concentration at HPCL. This can be due to the vehicular congestion and resuspended particles. The model performed well with mesoscale meteorology after all the sources and the entire region were considered. As Chembur has immense variation in topography, land use and geographical structure, as shown in Figure 1d, microscale meteorology varies with these factors. Further, mesoscale meteorology has been used for other analyses for air quality. The contours were plotted using the model for NOx and PM10 concentrations based on one-hour average values for one year for all sources in the study area. The same analysis was repeated only for industrial sources and vehicular emission sources separately. It was observed that the maximum concentration of PM10 was 71.8 µg/m3.This concentration was observed near the Chembur Gaothan area where the vehicular congestion was more intense. Around BPCL and HPCL area, PM10 concentration was around 50 µg/m3. The minimum concentration of PM10 was less than 42 µg/m3, and this concentration was observed near the southern boundary of the study region, which has been represented by deep violet color. The PM10 modeling was carried out using 30 µg/m3 as the background concentration. This concentration was estimated when modeled concentrations matched with observed concentrations. This background concentration also includes resuspended particulate matter (RSPM), which is induced by vehicular congestion and other factors. In PM10 modeling, this background concentration can be taken as a correction factor. The maximum concentration of NOx was observed to be 53 µg/m3 near the Ghatkopar–Mankhurd link road. This can be due to the vehicular congestion in Deonar Village, BPCL and HPCL area. Annual minimum concentration was less than 10 µg/m3 along the eastern, western and southern boundaries of the study area, which has been shown in deep violet color.
Emission load does not represent the rank-wise contribution of sources to the ambient concentration of the region. Hence, modeling was carried out for industrial sources, vehicle sources, and low duty diesel vehicles (LDDVs) to observe the relative source-wise contribution to the ambient air quality for future scope of implementation of control strategies and environmental management.

9.1. Results of Industrial Sources

For industrial emission sources modeling, four industries (BPCL, HPCL, RCFL and TPCL) were considered in Chembur. In this study, NOx and PM10 emissions were modeled for the year to find out the dominant source in the study domain.

9.2. Contribution of NOx and PM10 Concentration by Industries

NOx and PM10 emission load were 64% and 94%, respectively from industries in Chembur (Figure 3). However, it contributes less to ambient concentration in the study area. The southern part of the study area is dominated by industrial sources, and due to meteorology, the maximum concentration of NOx is 6.2 µg/m3, seen at HPCL. Also, NOx concentration is 4.8 µg/m3 in the southern part of BPCL. The maximum concentration of PM10 is 35 µg/m3 at BPCL and 33 to 34 µg/m3 at HPCL and RCFL. Table 2 shows the comparison of the simulated concentration of NOx and PM10 for industrial sources only and ambient simulated concentration of NOx and PM10 for this study area, respectively.

9.3. Results of Line Sources

In line source modeling, six roads in Chembur area have been considered. In the present study, NOx and PM10 emissions have been modeled for the year to find out the dominant sources in the study domain. Vehicular emission varies with time such as morning peak, evening peak, off peak and the lean peak of the day.

9.4. NOx and PM10 Concentration Contribution by Vehicles

NOx and PM10 emission loads from vehicles in Chembur are 17% and 3%, respectively (Figure 3). Nevertheless, these are contributing considerably to ambient concentration because they are ground emission sources. The northern part of the study area is dominated by vehicular sources and high density of vehicles. The maximum concentration of NOx was 43 µg/m3 at Chedda Nagar and 40.1 µg/m3 at Chembur Naka. The southern part of the study area has an inconsequential effect on vehicular pollution (only 2–5 µg/m3). At Chheda Nagar, the concentration of ambient air quality from all the sources was 54 µg/m3, while the concentration by vehicular sources was 41 µg/m3. Chheda Nagar is in the northern region of the study area, and this part is affected less by industrial emissions. The northern part of the study domain is highly populated compared to the southern part. The southern part has lesser contribution from vehicles. At Chheda Nagar, the maximum concentration of PM10 from vehicles was 37.84 µg/m3 and from the other sources was 70.8 µg/m3.

9.5. NOx and PM10 Concentration Modeling by Diesel Car and LDDV

Northern and western corners in this analysis are dominated by diesel cars and light duty diesel vehicles (LDDVs), while the other areas are almost free from NOx pollutant. The maximum concentration of NOx was 12.5 µg/m3 and it was found at Chheda Nagar. The entire area of Shramjivi Nagar showed NOx concentration of 7-9 µg/m3 sourced from diesel cars and LDDVs. Concentration contribution in ambient air quality from vehicles was 43.1 µg/m3. Thus, it can be concluded that diesel cars and LDDVs are contributing to one fourth of the line source emission. In PM10 concentration modeling, northern and western corners are affected by diesel cars and LDDVs. The maximum PM10 concentration of 32.8 µg/m3 was observed at Everard Nagar (Point 1 area) and 32.5 µg/m3 at Chembur Naka from diesel cars and LDDVs with 30 µg/m3 background concentration. Thus, it can be concluded that diesel cars and LDDVs are contributing to around 17% to line sources and 25% to the total concentration.

10. Summary and Conclusions

The aim of the study was to generate onsite meteorological data for usage in air quality modeling to see the feasibility of coarse resolution of WRF output in LMIC. It generated onsite and real time meteorological data, which was fed in AERMET, the pre-processor of the dispersion model AERMOD. AERMOD calculated concentrations for NOx and PM10 using compiled emission inventory for all the sources, namely industries and vehicles, of the study area. Air quality modeling results showed that in this particular case, the meteorological data from WRF output at mesoscale performed better than the observed meteorological data. WRF output could be a good option which may represent a better meteorology for the purpose of dispersion modeling. Also, this may be because industrial sources have a significant contribution in the region. The results were used to identify the critical areas and relative contribution of various sources to ambient air quality. The use of WRF model is very economical in resources and time.
The ambient concentration load of the study area is shown in Figure 7 for NOx and PM10. Here, for ambient concentration, the vehicular emissions are dominating in the region because these are ground level sources. Based on the study, following conclusions can be drawn:
⮚
Amongst total emissions, PM10 emission load was 3%, and NOx emission load was 17% from vehicles. Industr sources contributed 64% and 94% of NOx and PM10 load, respectively. The domestic sector contributed significantly to NOx emission as 18% of total emission load.
⮚
NOx emission load of industries was 64% of the total emission load, but it contributed only 25–30% of NOx concentration in the ambient air.
⮚
There was 94% contribution to total emission load of PM10 from industries in the study domain, but only 57% contribution to ambient air quality level.
⮚
NOx emission contribution from vehicles was 17% of total emission, but in ambient air quality it contributed only 26% of the total because it is a ground level source.
⮚
Vehicular PM10 emission contribution was 3% of the total emission load, but in ambient air quality it contributed 25% of the total ambient PM10 concentration.
⮚
At ambient concentration of NOx, diesel cars and LDDVs contributed one fourth of the line sources for this study domain.
Micro-meteorology may vary a lot due to the topography of the region. Topographic features could be one of the limitations of the use of WRF output meteorology for air quality modeling because it may not capture the high-rise buildings in the region. In AERMOD, only one meteorological profile can be used but meteorology may not be uniform for the entire region. It has been found many times that the observed meteorological parameters carry some measurement errors [54,55,56].
This will be very useful in the forecast, implementation, control, and management strategies for improving air quality and also for performing a risk analysis of different types of sources in the region as future scope. Also, this work can be useful to regulatory authorities to develop a framework for air quality management in LMIC. The shortage of research was (a) the meteorological data were available at only one location which can be observed at other locations and comparison can be done for the same and (b) suspended dust can be estimated and incorporated in emission inventory. Various physics option parameters to set up the WRF model and simulation of the microscale meteorology with a comparison of observed meteorology could be possible future researches. Health benefit analysis can also be done by estimating population exposure with air quality [57].

Author Contributions

Conceptual, methodology development, software operation and writing, A.K., review, data procurement and guidance, A.K.D. and R.S.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Ethical review and approval are not required for the study as the research does not involve humans.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data available on request.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Huff, G.; Angeles, L. Globalization, industrialization and urbanization in Pre-World War II Southeast Asia. Explor. Econ. Hist. 2011, 48, 20–36. [Google Scholar] [CrossRef]
  2. DeFries, R.; Pandey, D. Urbanization, the energy ladder and forest transitions in India’s emerging economy. Land Use Policy 2011, 27, 130–138. [Google Scholar] [CrossRef]
  3. Duh, J.-D.; Shandas, V.; Chang, H.; George, L.A. Rates of urbanisation and the resiliency of air and water quality. Sci. Total. Environ. 2008, 400, 238–256. [Google Scholar] [CrossRef]
  4. Patankar, A.; Trivedi, P. Monetary burden of health impacts of air pollution in Mumbai, India: Implications for public health policy. Public Health 2011, 125, 157–164. [Google Scholar] [CrossRef] [PubMed]
  5. World Health Organization (WHO). World Health Statistics, Department of Public Health, Environmental and Social Determinants of Health; WHO Press: Geneva, Switzerland, 2014. [Google Scholar]
  6. Kumar, A.; Gupta, I.; Brandt, J.; Kumar, R.; Kumar, A.; Patil, R.S. Air quality mapping using GIS and economic evalu-ation of health impact for Mumbai city, India. J. Air Waste Manag. Assoc. 2016, 66, 470–481. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Banerjee, T.; Barman, S.; Srivastava, R. Application of air pollution dispersion modeling for source-contribution assessment and model performance evaluation at integrated industrial estate-Pantnagar. Environ. Pollut. 2011, 159, 865–875. [Google Scholar] [CrossRef]
  8. Kumar, A.; Dikshit, A.K.; Fatima, S.; Patil, R.S. Application of WRF Model for Vehicular Pollution Modelling Using AERMOD. Atmos. Clim. Sci. 2015, 5, 57–62. [Google Scholar] [CrossRef] [Green Version]
  9. Kumar, A.; Patil, R.S.; Dikshit, A.K.; Islam, S.; Kumar, R. Evaluation of control strategies for industrial air pollution sources using American Meteorological Society/Environmental Protection Agency Regulatory Model with simulated meteorology by Weather Research and Forecasting Model. J. Clean. Prod. 2016, 116, 110–117. [Google Scholar] [CrossRef]
  10. Kumar, A.; Patil, R.S.; Dikshit, A.K.; Kumar, R. Comparison of predicted vehicular pollution concentration with air quality standards for different time periods. Clean Technol. Environ. Policy 2016, 18, 2293–2303. [Google Scholar] [CrossRef]
  11. Madala, S.; Satyanarayana, A.; Srinivas, C.; Kumar, M. Mesoscale atmospheric flow-field simulations for air quality modeling over complex terrain region of Ranchi in eastern India using WRF. Atmos. Environ. 2015, 107, 315–328. [Google Scholar] [CrossRef]
  12. Karagulian, F.; Belis, C.A.; Dora, C.F.C.; Prüss-Ustün, A.M.; Bonjour, S.; Adair-Rohani, H.; Amann, M. Contributions to cities’ ambient particulate matter (PM): A systematic review of local source contributions at global level. Atmos. Environ. 2015, 120, 475–483. [Google Scholar] [CrossRef]
  13. Thunis, P.; Degraeuwe, B.; Pisoni, E.; Trombetti, M.; Peduzzi, E.; Belis, C.; Wilson, J.; Clappier, A.; Vignati, E. PM2.5 source allocation in European cities: A SHERPA modelling study. Atmos. Environ. 2018, 187, 93–106. [Google Scholar] [CrossRef]
  14. Abhijith, K.; Gokhale, S. Passive control potentials of trees and on-street parked cars in reduction of air pollution exposure in urban street canyons. Environ. Pollut. 2015, 204, 99–108. [Google Scholar] [CrossRef] [PubMed]
  15. Gulia, S.; Shrivastava, A.; Nema, A.K.; Khare, M. Assessment of Urban Air Quality around a Heritage Site Using AERMOD: A Case Study of Amritsar City, India. Environ. Model. Assess. 2015, 20, 599–608. [Google Scholar] [CrossRef]
  16. Jiménez-Guerrero, P.; Jorba, O.; Baldasano, J.M.; Gassó, S. The use of a modelling system as a tool for air quality management: Annual high-resolution simulations and evaluation. Sci. Total. Environ. 2008, 390, 323–340. [Google Scholar] [CrossRef]
  17. Kumar, A.; Ketzel, M.; Patil, R.S.; Dikshit, A.K.; Hertel, O. Vehicular pollution modeling using the operational street pollution model (OSPM) for Chembur, Mumbai (India). Environ. Monit. Assess. 2016, 188, 349. [Google Scholar] [CrossRef] [PubMed]
  18. Kumar, A.; Patil, R.S.; Dikshit, A.K.; Kumar, R.; Brandt, J.; Hertel, O. Assessment of impact of unaccounted emission on ambient concentration using DEHM and AERMOD in combination with WRF. Atmos. Environ. 2016, 142, 406–413. [Google Scholar] [CrossRef]
  19. Mumovic, D.; Crowther, J.; Stevanovic, Z. Integrated air quality modelling for a designated air quality management area in Glasgow. Build. Environ. 2006, 41, 1703–1712. [Google Scholar] [CrossRef]
  20. Kumar, A.; Patil, R.S.; Dikshit, A.K.; Kumar, R. Assessment of Spatial Ambient Concentration of NH3 and its Health Impact for Mumbai City. Asia J. Atmos. Environ. 2019, 13, 11–19. [Google Scholar] [CrossRef]
  21. Patil, R.S.; Kumar, R.; Menon, R.; Kumar, M.; Sethi, V. Development of particulate matter speciation profiles for major sources in six cities in India. Atmos. Res. 2013, 132–133, 1–11. [Google Scholar] [CrossRef]
  22. Mohan, M.; Bhati, S.; Sreenivas, A.; Marrapu, P. Performance Evaluation of AERMOD and ADMS-Urban for Total Suspended Particulate Matter Concentrations in Megacity Delhi. Aerosol Air Qual. Res. 2011, 11, 883–894. [Google Scholar] [CrossRef] [Green Version]
  23. Seaman, N.L. Meteorological modeling for air-quality assessments. Atmos. Environ. 2000, 34, 2231–2259. [Google Scholar] [CrossRef]
  24. Sistla, G.; Zhou, N.; Hao, W.; Ku, J.-Y.; Rao, S.; Bornstein, R.; Freedman, F.; Thunis, P. Effects of uncertainties in meteorological inputs on urban airshed model predictions and ozone control strategies. Atmos. Environ. 1996, 30, 2011–2025. [Google Scholar] [CrossRef]
  25. Russell, A.; Dennis, R. NARSTO critical review of photo-chemical models and modelling. Atmos. Environ. 2000, 34, 2283–2324. [Google Scholar] [CrossRef]
  26. Gilliam, R.C.; Hogrefe, C.; Rao, S.T. New methods for evaluating meteorological models used in air quality applications. Atmos. Environ. 2006, 40, 5073–5086. [Google Scholar] [CrossRef]
  27. Borge, R.; Alexandrov, V.N.; Del Vas, J.J.; Lumbreras, J.; Rodríguez, E. A comprehensive sensitivity analysis of the WRF model for air quality applications over the Iberian Peninsula. Atmos. Environ. 2008, 42, 8560–8574. [Google Scholar] [CrossRef]
  28. Kumar, A.; Patil, R.S.; Dikshit, A.K.; Kumar, R. Application of WRF Model for Air Quality Modelling and AERMOD—A Survey. Aerosol Air Qual. Res. 2017, 17, 1925–1937. [Google Scholar] [CrossRef] [Green Version]
  29. Kesarkar, A.P.; Dalvi, M.; Kaginalkar, A.; Ojha, A. Coupling of the Weather Research Forecasting Model with AERMOD for pollutant dispersion modelling. A case study for PM10 dispersion over the Pune City. Atmos. Environ. 2006, 41, 1976–1988. [Google Scholar] [CrossRef]
  30. Kumar, A.; Patil, R.S.; Dikshit, A.K.; Kumar, R. Application of AERMOD for short-term air quality prediction with forecasted meteorology using WRF model. Clean Technol. Environ. Policy 2017, 19, 1955–1965. [Google Scholar] [CrossRef]
  31. Myers-Cook, T.; Mallard, J.; Mao, Q. Development of a WRF-AERMOD Tool for Use in Regulatory Application. In Proceedings of the 16th Conference on Air Pollution Meteorology, American Meteorological Society, Atlanta, GA, USA, 16–21 January 2010; pp. 1–10. [Google Scholar]
  32. Ground Water Information, Greater Mumbai District, Maharashtra; 1618/DB/20; Ministry of Water Resources, Central Ground Water Board: New Delhi, India, 2009.
  33. Gupta, I.; Kumar, R. Trends of particulate matter in four cities in India. Atmos. Environ. 2006, 40, 2552–2566. [Google Scholar] [CrossRef]
  34. Comprehensive Environmental Pollution Index (CEPI); Central Pollution Control Board, Ministry of Environment and Forests: New Delhi, India, 2010.
  35. Skamarock, W.C.; Klemp, J.B.; Dudhia, J.; Gill, D.O.; Barker, D.M.; Duda, M.G.; Huang, X.-Y.; Wang, W.; Powers, J.G. A Description of the Advanced Research WRF Version 3; NCAR Technical Note NCAR/TN-45; National Center for Atmospheric Research: Boulder, CO, USA, 2008. [Google Scholar]
  36. Henmi, T.; Flanigan, R.; Padilla, R. Development and Application of an Evaluation Method for the WRF Mesoscale Model; ARL-TR3657; Army Research Laboratory: Adelphi, MD, USA, 2005. [Google Scholar]
  37. National Center for Atmospheric Research (NCAR). University Corporation for Atmospheric Research (UCAR), Research Data Achieve (RDA), USA. Available online: https://rda.ucar.edu/datasets/ds083.2/ (accessed on 15 January 2020).
  38. Kumar, P.; Kishtawal, C.M.; Pal, P.K. Impact of satellite rainfall assimilation on Weather Research and Forecasting model predictions over the Indian region. J. Geophys. Res. Atmos. 2014, 119, 2017–2031. [Google Scholar] [CrossRef]
  39. Kumar, P.; Bhattacharya, B.K.; Nigam, R.; Kishtawal, C.M.; Pal, P.K. Impact of Kalpana-1 derived land surface albedo on short-range weather forecasting over the Indian subcontinent. J. Geophys. Res. Atmos. 2014, 119, 2764–2780. [Google Scholar] [CrossRef] [Green Version]
  40. Kumar, P.; Kishtawal, C.M.; Pal, P.K. Skill of regional and global model forecast over Indian region. Theor. Appl. Clim. 2016, 123, 629–636. [Google Scholar] [CrossRef]
  41. Sax, T.; Isakov, V. A case study for assessing uncertainty in local-scale regulatory air quality modeling applications. Atmos. Environ. 2003, 37, 3481–3489. [Google Scholar] [CrossRef]
  42. Singh, N.P.; Gokhale, S. A method to estimate spatiotemporal air quality in an urban traffic corridor. Sci. Total. Environ. 2015, 538, 458–467. [Google Scholar] [CrossRef] [PubMed]
  43. Thunis, P.; Georgieva, E.; Pederzoli, A. A tool to evaluate air quality model performances in regulatory applications. Environ. Model. Softw. 2012, 38, 220–230. [Google Scholar] [CrossRef]
  44. User’s Guide for Industrial Source Complex (ISC3) Dispersion Models, Volume I—User Instructions; EPA-454/B-95-003a; United State Environmental Protection Agency, Office of Air Quality Planning and Standards, Emissions, Monitoring, and Analysis Division, Research Triangle Park: Piedmont, NC, USA, 1999.
  45. Bandyopadhyay, A. Prediction of ground level concentration of sulfur dioxide using ISCST3 model in Mangalore industrial region of India. Clean Technol. Environ. Policy 2009, 11, 173–188. [Google Scholar] [CrossRef]
  46. Cimorelli, A.J.; Perry, S.G.; Venkatram, A.; Weil, J.C.; Paine, R.J.; Wilson Robert, B.; Lee, R.F.; Peters, W.D.; Brode, R.W.; Paumier, J.O. AERMOD: Description of Model Formulation; EPA-454/R-03-004; USEPA: Washington, DC, USA, 2004. [Google Scholar]
  47. Bhardwaj, K.S. Examination of the Sensitivity of Land Use Parameters and Population on the Performance of the AERMOD Model for and Urban Area. Civil Engineering. Master’s Thesis, University of Toledo, Toledo, OH, USA, 2005. [Google Scholar]
  48. Venkatram, A.; Brode, R.; Cimorelli, A.; Lee, R.; Paine, R.; Perry, S.; Peters, W.; Weil, J.; Wilson, R. A complex terrain dispersion model for regulatory applications. Atmos. Environ. 2001, 35, 4211–4221. [Google Scholar] [CrossRef]
  49. Venkatram, A.; Isakov, V.; Yuan, J.; Pankratz, D. Modeling dispersion at distances of meters from urban sources. Atmos. Environ. 2004, 38, 4633–4641. [Google Scholar] [CrossRef]
  50. Zou, B.; Zhan, F.B.; Wilson, J.G.; Zeng, Y. Performance of AERMOD at different time scales. Simul. Model. Pract. Theory 2010, 18, 612–623. [Google Scholar] [CrossRef]
  51. Zou, B.; Wilson, J.G.; Zhan, F.B.; Zeng, Y. Spatially differentiated and source-specific population exposure to ambient urban air pollution. Atmos. Environ. 2009, 43, 3981–3988. [Google Scholar] [CrossRef]
  52. Air Quality Assessment, Emission Inventory and Source Apportionment Studies: Mumbai; National Environmental Engineering Research Institute, CPCB: New Delhi, India, 2010.
  53. Air Quality Monitoring Project-Indian Clean Air Programme (ICAP). Emission Factor development for Indian Vehicles; Automotive Research Association of India CPCB-MoEF: New Delhi, India, 2008. [Google Scholar]
  54. Fedoseeva, E.V. An Estimate of the Error of Measurements of Radio Brightness Temperature in Radio-Heat Meteorological Parameters with Background Noise Compensation. Meas. Tech. 2015, 57, 1463–1468. [Google Scholar] [CrossRef]
  55. Finkelstein, P.L.; Sims, P.F. Sampling error in eddy correlation flux measurements. J. Geophys. Res. 2001, 106, 3503–3509. [Google Scholar] [CrossRef]
  56. Ulbrich, C.W.; Lee, L.G. Rainfall Measurement Error by WSR-88D Radars due to Variations in Z–R Law Parameters and the Radar Constant. J. Atmos. Ocean. Technol. 1999, 16, 1017–1024. [Google Scholar] [CrossRef]
  57. Sonawane, N.V.; Patil, R.S.; Sethi, V. Health benefit modelling and optimization of vehicular pollution control strategies. Atmos. Environ. 2012, 60, 193–201. [Google Scholar] [CrossRef]
Sustainability 13 04276 g001 550
Figure 1. (a) Mumbai Area, (b) WRF Domain (c) Study Area Chembur (Downloaded from Google Earth) (d) Terrain Map. Note: Figure (d) is given in Universal Transverse Mercator coordinate system. 1-Eastern Express Highway Road, 2-NGA Marg, 3- Ghatkopar-Mankhurd Road, 4-V N Purav Road, 5-R C Marg, 6- B D Patil Marg, A-RCFL, B-BPCL, C-HPCL, D-TPCL.
Figure 1. (a) Mumbai Area, (b) WRF Domain (c) Study Area Chembur (Downloaded from Google Earth) (d) Terrain Map. Note: Figure (d) is given in Universal Transverse Mercator coordinate system. 1-Eastern Express Highway Road, 2-NGA Marg, 3- Ghatkopar-Mankhurd Road, 4-V N Purav Road, 5-R C Marg, 6- B D Patil Marg, A-RCFL, B-BPCL, C-HPCL, D-TPCL.
Sustainability 13 04276 g001
Sustainability 13 04276 g002 550
Figure 2. Schematic Flow for the Study.
Figure 2. Schematic Flow for the Study.
Sustainability 13 04276 g002
Sustainability 13 04276 g003 550
Figure 3. Emission load of NOx and PM10 in Chembur (in kg/day).
Figure 3. Emission load of NOx and PM10 in Chembur (in kg/day).
Sustainability 13 04276 g003
Sustainability 13 04276 g004 550
Figure 4. Comparison of Hourly Temperature of WRF with Observed Data.
Figure 4. Comparison of Hourly Temperature of WRF with Observed Data.
Sustainability 13 04276 g004
Sustainability 13 04276 g005 550
Figure 5. Annual Wind Rose Simulated by (a) WRF and (b) RCFL Observed Data.
Figure 5. Annual Wind Rose Simulated by (a) WRF and (b) RCFL Observed Data.
Sustainability 13 04276 g005
Sustainability 13 04276 g006 550
Figure 6. Annual Concentration of (a) NOx and (b) PM10 (in µg/m3) using observed meteorology for Chembur Area.
Figure 6. Annual Concentration of (a) NOx and (b) PM10 (in µg/m3) using observed meteorology for Chembur Area.
Sustainability 13 04276 g006
Sustainability 13 04276 g007 550
Figure 7. Comparison between emission load and ambient concentration load of NOx and PM10 in Chembur area.
Figure 7. Comparison between emission load and ambient concentration load of NOx and PM10 in Chembur area.
Sustainability 13 04276 g007
Table
Table 1. Monthly and Annual Comparison of Simulated Concentration with Ambient Observed Concentration (µg/m3).
Table 1. Monthly and Annual Comparison of Simulated Concentration with Ambient Observed Concentration (µg/m3).
Month-LocationNOxPM10
SimulatedObserved% ErrorSimulatedObserved% Error
Jan-BPCL292611.55362−14.5
Jan-HPCL362263.65565−15.4
Feb-BPCL36345.954540.0
Feb-HPCL432853.65768−16.2
Mar-BPCL28273.75356−5.4
Mar-HPCL392556.05664−12.5
Apr-BPCL3134−8.860559.1
Apr-HPCL2728−3.6655420.4
May-BPCL2132−34.45158−12.1
May-HPCL282321.75556−1.8
Jun-BPCL1516−6.3513545.7
Jun-HPCL2126−19.25558−5.2
Jul-BPCL820−60.0494314.0
Jul-HPCL1518−16.753506.0
Aug-BPCL720−65.0494314.0
Aug-HPCL18175.9554814.6
Sep-BPCL221915.8514124.4
Sep-HPCL262123.854540.0
Oct-BPCL392556.0524223.8
Oct-HPCL403225.055550.0
Nov-BPCL312429.255533.8
Nov-HPCL2326−11.55059−15.3
Dec-BPCL282512.056560.0
Dec-HPCL1521−28.64962−21.0
Annual-BPCL24.625.2−3.452.849.88.6
Annual-HPCL27.623.914.254.957.8−3.9
Table
Table 2. Comparison of simulated concentration from industries with ambient simulated concentration of NOx and PM10.
Table 2. Comparison of simulated concentration from industries with ambient simulated concentration of NOx and PM10.
PollutantLocationSimulated Concentration for Industries (µg/m3)Ambient Simulated Concentration (µg/m3)Contribution of Industrial Source to Ambient Air Quality
NOxBPCL4.828.217%
HPCL6.121.728%
PM10BPCL35.556.962%
HPCL33.252.164%
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kumar, A.; Dikshit, A.K.; Patil, R.S. Use of Simulated and Observed Meteorology for Air Quality Modeling and Source Ranking for an Industrial Region. Sustainability 2021, 13, 4276. https://doi.org/10.3390/su13084276

AMA Style

Kumar A, Dikshit AK, Patil RS. Use of Simulated and Observed Meteorology for Air Quality Modeling and Source Ranking for an Industrial Region. Sustainability. 2021; 13(8):4276. https://doi.org/10.3390/su13084276

Chicago/Turabian Style

Kumar, Awkash, Anil Kumar Dikshit, and Rashmi S. Patil. 2021. "Use of Simulated and Observed Meteorology for Air Quality Modeling and Source Ranking for an Industrial Region" Sustainability 13, no. 8: 4276. https://doi.org/10.3390/su13084276

APA Style

Kumar, A., Dikshit, A. K., & Patil, R. S. (2021). Use of Simulated and Observed Meteorology for Air Quality Modeling and Source Ranking for an Industrial Region. Sustainability, 13(8), 4276. https://doi.org/10.3390/su13084276

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop