The status of ambient air quality in terms of PM10, PM2.5, SO2, NO2, NO, CO, BTEX over different years is given in Fig. 2 for Chandrapur and in Fig. 3 for Nagpur.
5.1 Chandrapur
Average PM10 and PM2.5 concentration is almost constant throughout the study years with variability in maximum concentration (refer Fig. 2). Further, PM10 is higher than its threshold limit during all the years, whereas average PM2.5 approaches the threshold concentration. Coarse particulate pollution is significant and fine particulate matter pollution is moderate as average PM2.5/PM10 ratio is observed to be < 0.5. Opencast mining activities cause coarse particulate matter emissions. Combustion due to household heating, coal-fired power plants and transportation emits fine particulate matter. As evident, the dominance of coarse particulate matter suggests the dominance of mining activity followed by a moderate contribution of fine particulate matter by the power plant and other combustion-related emissions. Unlike particulate matter, the monotonous trend is observed in most of the gaseous pollutants. CO concentration is observed to be decreasing, whereas NO2, SO2 and O3 concentrations are increasing. The non-monotonous trend is observed in benzene, toluene, ethylbenzene and xylene. All the gaseous pollutants such as CO, NO2, O3, SO2, benzene, toluene, ethylbenzene and xylene are well below their standard limits.
Hourly diurnal variations in all the pollutant concentrations in Fig. 3 show that average concentration does not vary significantly, however, the maximum concentration and upper quartiles significantly vary diurnally. The diurnal variation in PM10, PM2.5 and CO concentration is bi-modal with a high concentration in the evening and night hours. The concentrations peak during late-evening and night hours. A peak is again observed in the morning hours i.e. 7–10 am with troughs during day-time. During the morning hours and evening (till late evening hours), the movement of vehicular traffic is usually observed. The coincidence of maxima and minima in PM10, PM2.5 and CO with vehicular traffic movement is therefore observed. Calm atmospheric conditions led to high concentration during early night hours. Like PM10, PM2.5 and CO, the pollutants NO, benzene, toluene and xylene concentrations show similar but moderate variations. NO2 concentration does not show any significant variations. Ethylbenzene shows higher concentration during late noon hours lasting till evening. As expected, O3 concentration shows higher levels in day-time as compared to evening and night-time. A strong correlation of O3 concentration with solar radiation is observed. Considerable variations are observed in SO2 concentration minimums and maximums along with upper quartiles. Day-time high SO2 concentration suggests the possible contribution of emissions from the power plant.
Monthly variations in all the parameters (refer Fig. 4) show that particulate matter is lower in monsoon months (mid-June till September) as compared to winter (December-February) and summer months (march-mid June). Rainfall scavenging effect is visible for PM10 and PM2.5 concentration. CO concentration shows similar variations as PM concentration. NO2 concentration is significantly high in monsoon months, even higher than the concentration in winter and summer months. Washout effect in wet months is not evident in NO2 concentration. Average NO concentration does not show significant monthly variation, whereas upper quartile of NO concentration shows variations with a high concentration in winter months (January February and December) and post-monsoon months (October and November). Summer months show relatively lower concentration than monsoon months. Benzene, toluene, ethylbenzene and xylene show similar variations as NO2 and NO concentration. Ethylbenzene, however, shows the lower concentration in post-monsoon months. Like NO concentration, average O3 concentration does not show variations but upper quartile of O3 concentration shows significant variations with a high concentration in February-July and October.
There is a strong positive correlation between BTEX (refer Table 1), which are also observed to be correlated positively with CO, PM10, PM2.5, NO and NO2. A moderate correlation of CO is observed with PM10, PM2.5 and strong correlation with BTEX, NO and NO2, which suggests a similarity in their emission sources. O3 is negatively correlated with relative humidity only. Although O3 data with an hourly resolution is not correlated with other pollutants, the correlation was assessed for other resolutions like monthly and hourly averages. Correlation of monthly O3 time series is however observed to be not significantly correlated with its precursors. Hourly variation time series showed the correlation of O3 with CO, NO, NO2 and BTEX as -0.48,-0.92,-0.11 and − 0.37 which suggests that O3 formation mechanism is driven by both VOCs (here BTEX), CO and NOx. On a finer time scales or at high resolution, O3 formation mechanism can be understood, whereas on the low resolution, the effect of confounding factors such as seasonality distorts the correlation of O3 with its precursor variables. This suggests that the sources of these pollutants are similar in nature. The correlation of particulate matter with ozone precursors indicates the similarity in the sources of these pollutants.
Table 1
Correlation of Various Air Pollutant Concentrations in Chandrapur
Species | Ben | Tol | EBen | X | WD | WS | Temp | CO | PM10 | PM2.5 | RH | SR | O3 | NO | NO2 |
Ben | 1 | | | | | | | | | | | | | | |
Tol | 0.09 | 1 | | | | | | | | | | | | | |
Eth-B | 0.90 | 0.88 | 1 | | | | | | | | | | | | |
X | 0.69 | 0.20 | 1 | 1 | | | | | | | | | | | |
WD | -0.01 | -0.01 | -0.03 | -0.01 | 1 | | | | | | | | | | |
WS | 0.02 | -0.04 | -0.03 | 0.01 | 0.19 | 1 | | | | | | | | | |
Temp | -0.01 | -0.04 | 0.03 | -0.03 | 0.12 | 0.11 | 1 | | | | | | | | |
CO | -0.02 | 0.09 | 0.01 | -0.01 | -0.09 | -0.14 | -0.07 | 1 | | | | | | | |
PM10 | 0.00 | 0.07 | 0.01 | 0.01 | -0.14 | -0.16 | -0.15 | 0.37 | 1 | | | | | | |
PM2.5 | -0.01 | 0.05 | 0.03 | 0.00 | -0.11 | -0.13 | -0.13 | 0.26 | 0.53 | 1 | | | | | |
RH | -0.03 | -0.02 | 0.01 | -0.02 | 0.21 | -0.09 | 0.04 | -0.09 | -0.25 | -0.06 | 1 | | | | |
SR | | | | | | | | | | | | 1 | | | |
O3 | -0.02 | -0.03 | 0.00 | -0.02 | -0.03 | -0.06 | 0.10 | -0.04 | -0.10 | -0.08 | -0.01 | | 1 | | |
NO | -0.01 | 0.07 | 0.02 | 0.01 | 0.00 | -0.17 | -0.09 | 0.33 | 0.29 | 0.26 | 0.01 | | -0.15 | 1 | |
NO2 | 0.02 | 0.05 | -0.01 | -0.02 | -0.03 | -0.13 | 0.09 | 0.18 | 0.12 | 0.03 | 0.05 | | 0.19 | -0.05 | 1 |
Ben: Benzene, Tol: Toluene, Eth-B: Ethyl Benzene, X: Xylene, WD: Wind direction, WS: Wind speed, Temp: Temperature, RH: Relative humidity, SR: Solar radiation |
5.2 Nagpur
Like Chandrapur, Average PM10 and PM2.5 concentration in Nagpur is also almost constant throughout the years with variability in the upper quartiles (refer Fig. 5). PM10 is higher than its threshold limit during all the years, whereas average PM2.5 approaches the threshold concentration. PM2.5/PM10 ratio is observed to be 0.5–0.68 during four years suggesting the presence of the fine particulate matter. Traffic emissions are dominant in Nagpur followed by household combustion and power plant emissions causing the significance of fine particulate matter. The trend is not monotonous in the gaseous pollutants such as NO2, SO2, O3, benzene, toluene, ethylbenzene and xylene, whereas CO concentration is decreasing like Chandrapur. All the gaseous pollutants except toluene and xylene are well below the CPCB threshold.
Like Chandrapur, hourly variations in all the pollutant concentrations show that average concentration do not vary significantly, however, the maximum concentration and upper quartile significantly vary diurnally (refer Fig. 6). A similar bi-modal diurnal behaviour of PM10 and PM2.5 concentration is observed in Nagpur like Chandrapur with a high concentration in the evening and night hours. As compared to Chandrapur, however, the bi-modality is not pronounced. The concentrations peak during late-evening and night hours. During the morning hours and evening (till late evening hours), the movement of vehicular traffic is usually observed. The maxima and minima coincided with the traffic movement extent. Benzene, toluene, ethylbenzene and xylene show similar variations with high concentration during late-evening hours till 11’o clock in the night. This can be attributed to prevailing calm weather conditions in the late evening to early night hours. The similarity in the variation pattern of CO, NO, NO2 and BTEX compounds indicates their similarity in the source and dispersion pattern (Tiwari et al., 2010).The reactive species are generally not abundant during day time due to their involvement in photochemical activity (Tiwari et al., 2010). The diurnal variability in O3 concentration concurs well with the solar radiation variability. Like Chandrapur, day-time high SO2 concentration suggests the possible contribution of emissions from the power plant.
Monthly variations plot (Fig. 7) shows that particulate matter is lower in monsoon months due to the washout effect of rains. CO, NO and O3 concentration does not show significant monthly variations and almost remain constant throughout the year. Higher NO concentration is observed in August. Summer months show relatively lower concentration than monsoon months. Lower O3 concentration is observed in monsoon months. NO2 concentration is initially lower in July and August but then increase in September and again get reduced in October. Benzene, toluene, ethylbenzene and xylene show similar variations as NO2 and NO concentration.
Correlation of pollutants throws some insight into the variability in their emission sources. It can be seen from Table 2 that in Nagpur, O3 is not correlated significantly with any of the precursor variables. PM10 is correlated well with PM2.5. Benzene is correlated significantly with toluene, xylene and ethylbenzene. With other pollutants like CO, NO and NO2 also, a good correlation of benzene is observed. PM10 and PM2.5 are however moderately correlated with benzene. The correlation of BTEX compounds with CO, NO and NO2 suggest their common source. CO is correlated well with NO and NO2 and moderately with PM10. CO is mainly emitted from the vehicular emissions, which is, therefore, the common source of these pollutants in the area. SO2 is not correlated with any of the pollutants. O3 is moderately correlated with temperature and solar radiation which suggests the photochemical reactions are governing the O3 levels in the area.
Table 2
Correlation of Various Air Pollutant Concentrations in Nagpur
Species | B | Eth-B | X | T | WD | WS | Temp | CO | PM10 | PM2.5 | RH | SR | O3 | NO | NO2 | SO2 |
Ben | 1 | | | | | | | | | | | | | | | |
Eth-B | 0.88 | 1 | | | | | | | | | | | | | | |
X | 0.90 | 0.87 | 1 | | | | | | | | | | | | | |
T | 0.92 | 0.83 | 0.96 | 1 | | | | | | | | | | | | |
WD | 0.19 | 0.21 | 0.19 | 0.22 | 1 | | | | | | | | | | | |
WS | -0.04 | -0.04 | -0.04 | -0.02 | 0.00 | 1 | | | | | | | | | | |
Temp | -0.27 | -0.22 | -0.22 | -0.15 | -0.04 | 0.10 | 1 | | | | | | | | | |
CO | 0.86 | 0.76 | 0.83 | 0.86 | 0.24 | -0.01 | -0.11 | 1 | | | | | | | | |
PM10 | 0.65 | 0.51 | 0.57 | 0.57 | 0.09 | -0.09 | -0.20 | 0.53 | 1 | | | | | | | |
PM2.5 | 0.52 | 0.35 | 0.41 | 0.40 | 0.02 | -0.07 | -0.24 | 0.38 | 0.74 | 1 | | | | | | |
RH | 0.07 | 0.09 | 0.08 | 0.06 | 0.10 | 0.04 | -0.57 | 0.00 | -0.12 | -0.07 | 1 | | | | | |
SR | -0.24 | -0.23 | -0.23 | -0.21 | -0.11 | 0.11 | 0.55 | -0.18 | -0.11 | -0.14 | -0.38 | 1 | | | | |
O3 | -0.14 | -0.16 | -0.11 | -0.13 | -0.18 | -0.02 | 0.34 | -0.12 | 0.10 | 0.07 | -0.47 | 0.20 | 1 | | | |
NO | 0.71 | 0.73 | 0.75 | 0.69 | 0.22 | -0.06 | -0.20 | 0.80 | 0.40 | 0.26 | 0.07 | -0.22 | -0.14 | 1 | | |
NO2 | 0.64 | 0.60 | 0.68 | 0.67 | 0.19 | -0.09 | -0.09 | 0.71 | 0.45 | 0.30 | -0.11 | -0.27 | -0.12 | 0.64 | 1 | |
SO2 | 0.03 | 0.00 | 0.02 | 0.03 | -0.19 | 0.00 | 0.22 | 0.06 | 0.12 | 0.09 | -0.35 | 0.17 | 0.26 | 0.01 | 0.09 | 1 |
Ben: Benzene, Eth-B: Ethyl Benzene, X: Xylene, WD: Wind direction, WS: Wind speed, Temp: Temperature, RH: Relative humidity, SR: Solar radiation |
5.3 BTEX ratios
The ratio of BTEX compounds has been studied in the literature to assess the nature of emissions in the area. Xylene to ethylbenzene ratio (X/Eth-B) has been used to assess the contribution of distant or local emission sources. A high ratio suggests the local emissions, whereas a low ratio indicates the contribution from distant sources (Tiwari et al., 2010). The average ratio in Chandrapur is observed to be 4 and in Nagpur, it is 4.65 suggesting the local hydrocarbon source adding to the BTEX concentration in the two cities. Toluene to benzene ratio (T/B) has been used in the literature to assess the relative contribution of BTEX sources (Breton et al., 2020). A low T/B ratio suggests contribution mainly derived from vehicles, whereas high values indicate the contribution other than the vehicles. T/B ratio is observed to be 7.73 and 4.63, for Chandrapur and Nagpur respectively. Lower T/B ratio in Nagpur suggests the contribution of vehicular emissions (Lee et al., 2002; Garg et al., 2019). Higher T/B ratio in Chandrapur suggests the closeness of point sources (Tiwari et al., 2010).
5.4 Application of random forest model
Random forest model is applied with O3 concentration as the response variable and NO2, NO, SO2, CO, BTEX, PM10 and PM2.5 as independent variables. With 70:30 as the training and testing ratio and 10 fold cross-validations, the results are obtained. The model uses R2, Root Mean Square Error (RMSE) and Mean Absolute Percentage Error (MAPE) as the performance measures for its evaluation. The prediction results show R2 = 0.92, RMSE = 8.9 and MAPE = 11.6 for Chandrapur and R2 = 0.96, RMSE = 6.7 and MAPE = 9.4 for Nagpur, which suggests that the predictions are quite reasonable. The feature importance assessment matrix which is the in-built feature of the random forest model is shown in Fig. 8. It can be seen that PM10, NO and CO along with solar radiation are observed to be highly important variable governing the O3 dynamics in Chandrapur. Meteorological variables are relatively less important than precursor variables. BTEX and PM2.5 are the least important variables. In Nagpur, wind direction, relative humidity, temperature and toluene are important. NO2 is more important than CO in governing the O3 dynamics. In general, BTEX is least important in governing the O3 variability in both the cities. The primary pollutants such as CO, NO and NO2 are more important and modulate the O3 dynamics. In Nagpur, meteorology plays an important role in O3 formation, whereas in Chandrapur, only solar radiation is more important amongst meteorological variables.
Further, to attribute the levels of pollutants to emission sources, a qualitative analysis is carried out for the two time periods as; with and without emission activities. In India, due to COVID19, the complete lockdown was initiated during March 22, 2020 to April 20, 2020. During this period, only essential services were retained with the restriction on vehicular movement. Heavy-duty vehicles were allowed only to carry the essential services loads that too were plying outside the city centres. Traffic-related activities, construction, solid waste burning were completely banned. The only allowed activities related to the air pollution emissions were power plant and household emissions from combustion. The contention of analyzing the data during the lockdown and during the normal period is that if the banned activities are contributing to the air pollution in the area, their shutting down will considerably reduce the air pollution concentrations during the lockdown period. No significant reduction in the air pollution levels may also suggest the contribution from running activities such as power plants and domestic heating. A comparison analysis of pollutant concentration is therefore carried for the lockdown period during 2020 and the corresponding pollutant concentration during 2019 for the same period i.e. March 24, 2019 to April 20, 2019.
It can be seen from Table 3 that benzene, toluene, xylene has reduced by 70–88% in Chandrapur and by 77–98% in Nagpur during the lockdown period. This suggests that traffic-related emissions significantly contributed to the BTEX in these two cities. In Chandrapur, PM10 and PM2.5 reduced by about 37% and 18% and in Nagpur, both reduced by about 51% and 54%, respectively. Although the reduction in PM is significant at 5% level of significance in both the cities, it is not pronounced like BTEX compounds. This indicates that apart from the traffic emissions, construction, power plant emissions and domestic combustion also add to the PM levels in the cities. CO is reduced by 30% in Chandrapur and by 60% in Nagpur. NO and NO2 has reduced considerably in Nagpur, but in Chandrapur, NO2 rather has increased and reduction in NO is not significant. O3 is observed to be increased in the lockdown period in Chandrapur, whereas in Nagpur the significant reduction of about 50% is observed. This suggests the influence of traffic emissions in Nagpur and dominance of non-traffic related emissions, mainly power plant and mining activities in Chandrapur.
Table 3
Pollutant Concentrations during COVID-19 Lockdown Period and the Corresponding Period in 2019
Species | Chandrapur | | Nagpur |
2020 | 2019 | %reduction | significance | | 2020 | 2019 | %reduction | significance |
Benzene | 0.1 | 0.34 | 70.6 | * | | 1.32 | 5.76 | 77.1 | * |
Toluene | 0.2 | 1.72 | 88.4 | * | | 0.35 | 22.61 | 98.5 | * |
Eth-Benzene | - | - | - | | | 4.12 | 4.89 | 15.7 | ns |
Xylene | 0.2 | 1.29 | 84.5 | * | | 1.67 | 16.46 | 89.9 | * |
PM10 | 81.8 | 129.63 | 36.9 | * | | 47.47 | 96.60 | 50.9 | * |
PM2.5 | 38.2 | 46.72 | 18.2 | * | | 23.18 | 50.01 | 53.6 | * |
CO | 0.3 | 0.43 | 30.2 | * | | 0.46 | 1.15 | 60.0 | * |
NH3 | 75.8 | 25.20 | - | * | | 32.23 | 41.11 | 21.6 | * |
NO | 3.4 | 3.83 | 11.2 | | | 3.46 | 13.81 | 74.9 | * |
NO2 | 13.3 | 12.47 | - | | | 21.85 | 43.97 | 50.3 | * |
SO2 | 25.1 | 8.52 | - | * | | - | 18.72 | - | |
O3 | 22.2 | 20.59 | | * | | 58.14 | 116.03 | 49.9 | * |
* significant at 0.05 level of significance, ns-not significant |
5.5 Heath risk assessment
Hazard quotient is obtained for different pollutants based on the assumption of chronic exposure of the population. Assuming the exposure time (ET) of 24 hours/day, duration of exposure (DE) as 30 years and exposure frequency (EF) of 312 days/year, exposure duration (ED) is computed. The averaging time (AT) is estimated by multiplying the duration of exposure and 365 days/year. Average inhalation rate (IR) for an adult is 20 m3/day (Van den Berg, 1995). Bodyweight of average Indian male is 70 kg and the average Indian female is 60 kg. Average body weight of 65 kg is therefore considered. Reference concentration for different pollutants is considered based on their respective standard provided by CPCB (2009). The reference concentration for BTEX compounds, NO2, SO2, CO, O3, PM10 and PM2.5 is considered as 5, 40, 50, 2, 100, 60 and 40 µg m− 3, respectively. It can be seen from Table 4 that hazard quotient for PM2.5 and PM10 is relatively higher than other pollutants and ranges between 0.33–0.64 for Chandrapur. All other pollutants have hazard quotient < 0.17 for Chandrapur. Higher hazard quotient is observed for toluene in Nagpur ranging from 0.40–1.35 during 2016–2019. For PM2.5 and PM10, it ranges between 0.29–0.59. Comparing the hazard quotient of Nagpur with Chandrapur, it is higher for ethylbenzene, xylene, benzene and NO2 in Nagpur than in Chandrapur. In a study on the morbidity and mortality due to respiratory and cardiovascular diseases, the excess number of mortality and morbidity in both the cities is observed to be increasing (Maji et al., 2016). Overall, the hazard index is 1.22–1.67 for Chandrapur and 3.07–4.35 for Nagpur. Hazard index > 1 in both cities suggests the health hazard to the residents living in the area.
Table 4
Hazard Quotient of Various Pollutants at Chandrapur and Nagpur
Pollutant | Chandrapur | Nagpur |
2016 | 2017 | 2018 | 2019 | 2016 | 2017 | 2018 | 2019 |
CO | 0.09 | 0.09 | 0.07 | 0.05 | 0.07 | 0.05 | 0.15 | 0.13 |
O3 | 0.03 | 0.03 | 0.05 | 0.05 | 0.15 | 0.09 | 0.11 | 0.15 |
NO | 0.06 | 0.04 | 0.02 | 0.04 | 0.08 | 0.10 | 0.11 | 0.16 |
NO2 | 0.07 | 0.10 | 0.13 | 0.10 | 0.36 | 0.18 | 0.25 | 0.34 |
SO2 | 0.04 | 0.05 | 0.05 | 0.06 | 0.06 | 0.07 | 0.07 | 0.06 |
PM2.5 | 0.48 | 0.41 | 0.33 | 0.30 | 0.29 | 0.49 | 0.32 | 0.31 |
PM10 | 0.61 | 0.64 | 0.54 | 0.49 | 0.38 | 0.52 | 0.42 | 0.38 |
Benzene | 0.01 | 0.01 | 0.05 | 0.01 | 0.14 | 0.12 | 0.36 | 0.31 |
Toluene | 0.04 | 0.10 | 0.04 | 0.06 | 0.64 | 0.40 | 1.35 | 1.03 |
Eth-Benzene | 0.06 | 0.10 | - | - | 0.23 | 0.34 | 0.30 | 0.34 |
Xylene | 0.17 | 0.08 | 0.14 | 0.05 | 0.67 | 1.06 | 0.91 | 0.81 |
Hazard Index (HI) | 1.65 | 1.67 | 1.42 | 1.22 | 3.07 | 3.43 | 4.35 | 4.01 |