Abstract
Fine particulate matter (PM2.5, aerodynamic diameter ≤2.5 µm) impacts the climate, reduces visibility and severely influences human health. The Indo-Gangetic Plain (IGP), home to about one-seventh of the world’s total population and a hotspot of aerosol loading, observes strong enhancements in the PM2.5 concentrations towards winter. We performed high-resolution (12 km × 12 km) atmospheric chemical transport modeling (WRF-Chem) for the post-monsoon to winter transition to unravel the underlying dynamics and influences of regional emissions over the region. Model, capturing the observed variations to an extent, reveals that the spatial distribution of PM2.5 having patches of enhanced concentrations (≥100 µgm−3) during post-monsoon, evolves dramatically into a widespread enhancement across the IGP region during winter. A sensitivity simulation, supported by satellite observations of fires, shows that biomass-burning emissions over the northwest IGP play a crucial role during post-monsoon. Whereas, in contrast, towards winter, a large-scale decline in the air temperature, significantly shallower atmospheric boundary layer, and weaker winds lead to stagnant conditions (ventilation coefficient lower by a factor of ~4) thereby confining the anthropogenic influences closer to the surface. Such changes in the controlling processes from post-monsoon to winter transition profoundly affect the composition of the fine aerosols over the IGP region. The study highlights the need to critically consider the distinct meteorological processes of west-to-east IGP and changes in dominant sources from post-monsoon to winter in the formulation of future pollution mitigation policies.
Similar content being viewed by others
Introduction
The Indo-Gangetic Plain (IGP) region is considered to be a global hotspot of elevated aerosol loading, and is among few regions of the world, which have been experiencing enhancements1,2,3,4. Atmospheric processes occurring in this region have been shown to affect atmospheric composition, chemistry and climate over regional as well as global scales5,6,7,8. Fine particulate matter (PM2.5) in this region is estimated to cause severe health implications including premature mortalities9,10,11,12. Numerous efforts based on in-situ and satellite-based observations have revealed strong enhancements in the PM2.5 over this region towards winter, typically every year13,14,15,16. PM2.5 concentrations over the IGP region are observed to exceed the standards of World Health Organization (WHO) as well as the National Ambient Air Quality Standards (NAAQS) of India, attributed to the combined effects from a variety of anthropogenic and biomass-burning emissions and meteorological conditions17,18,19,20,21,22,23,24. The accumulation of aerosols in this region is further supported by the topography of Himalaya that extends from northwest to southeast and forms a corridor along the IGP.
The changes in meteorology and atmospheric dynamics can offset the benefits availed from reduction in anthropogenic emissions. Such meteorological influences, besides the elevated emission source strengths and upwind biomass-burning, highly complicate the detailed understanding of aerosol behavior over the IGP18,22,23,24. The implementation of an odd-even scheme for vehicular traffic aiming to reduce its emissions in Delhi, India’s capital territory, did not lead to a significant improvement in the air quality25,26, which highlighs the crucial need of studies on the spatio-temporal distribution of PM2.5 over the complex environment of IGP region. Surface-based measurements did reveal the impacts of meteorology on aerosol properties over this region, however, such point observations are representative of a limited area and hence need to be corroborated with model results27 for portraying the regional picture. In this direction, a global model has been applied to explore the influences of meteorological conditions on the spatio-temporal distribution of PM2.5 over this region, recently23.
Regional modeling studies, exploiting the potential of high-resolution input emissions together with a well-resolved meteorology28,29,30,31,32, are essential to unravel the effects of dynamics and emissions (biomass-burning versus anthropogenic) on the widespread PM2.5 build-up across the IGP from post-monsoon to winter transition. In this direction, here we performed high-resolution (12 km × 12 km) simulations using the Weather Research and Forecasting model coupled with Chemistry (WRF-Chem)33, configured on the basis of the previous evaluations of meteorology, wind patterns, and chemical fields over the South Asian region30,32,34,35 (Table S1). Figure 1a shows the study region (also the domain of the WRF-Chem model), and the location of ground-based observations. Model simulations are performed for the period of October–December 2016, on the basis of previously reported extremely high fine particulate loading during the post-monsoon and winter over the IGP region36,37,38. Besides the reference simulation, called as “WRF-Chem”, two additional simulations have been performed by alternatively switching off the anthropogenic (anthro_off) and biomass burning (fire_off) emissions in the model.
Results
Aerosol distribution: Model versus observations
Aerosol optical depth (AOD) at 550 nm is simulated to be higher over the IGP region (0.5–1) in the results of the standard simulation, which is in good agreement with the satellite-based observations from moderate resolution imaging spectroradiometer (MODIS) and multiangle imaging spectroradiometer (MISR) (Fig. S1). The domain-wide mean AOD values simulated by the model (e.g. 0.27 ± 0.14 in October) are found to be comparable with that in the MODIS and MISR observations (0.29 ± 0.21 and 0.24 ± 0.14, respectively for the same period). Some differences are, however, seen not only between model and satellite, but also within different satellite-based instruments, as reported for retrievals of other trace species over this region29. Nevertheless, the spatial variability between model and satellite observations is found to be correlated with r values of 0.67–0.74 (with MODIS) and 0.71–0.85 (with MISR) during different months over the period of study. Further, model captures the variations in the AOD (r = ~0.6) at Kanpur and Jaipur stations as obtained from the Aerosol Robotic Network (AERONET) observations in this region, to a reasonable extent (Fig. 1b) with underestimated magnitudes. AOD uncertainties in the model are also found to be associated with the dust transport from the domain boundaries (based on additional sensitivity simulations with varying dust in the boundary conditions, not shown here). Nevertheless, the dust transport from domain boundaries is seen to cause lesser impacts on the surface level PM2.5 in the simulations (also see the section-Methods). Additionally, the model simulated day-to-day variations in the PM2.5 concentrations correlate with the ground-based observations at Delhi, Lucknow, Kanpur, Varanasi, and Patna stations in the IGP region with r values in the range of 0.33–0.73 (Figs. 1c and S2). However, model tends to underestimate the PM2.5 levels with Normalized Mean Biases (NMB) in the range of 21–35%, except at an industrial site Kanpur in the central IGP where analogous to the AOD, PM2.5 levels were underestimated more strongly (NMB = 59%). Model is seen to reasonably capture the temporal variations and magnitudes, and stronger bias at one of the stations is also within the range of model bias (40–60%) reported in an earlier evaluation of this model over this region during winter28. The underestimations of aerosol loadings by model are suggested to be resulting from the uncertainties associated with the input emissions28,39 and fine-scale dynamics30,40, besides some small-scale fires and local sources such as trash and wood burning unaccounted or underestimated in the inputs. Furthermore, models can have a dilution effect over grid boxes as compared to the observations sampling directly the urban and rural air masses at a location/point, which we have tried to minimize by performing simulations at higher spatial resolution (12 km × 12 km). Overall, the WRF-Chem model shows the ability to simulate the distribution of aerosols over northern Indian region, besides some differences in the magnitudes attributed to the aforementioned factors.
The model, capturing the observed variations, reveals strong heterogeneity in the PM2.5 distribution over the northern Indian subcontinent (Fig. 2) with elevated levels (100–200 µgm−3) over the IGP region. During October, high PM2.5 levels (≥100 µgm−3) were most pronounced over the western IGP including Delhi, Haryana, and Punjab. Additionally, in the same period, similarly high PM2.5 concentrations are simulated over the eastern parts of the IGP. A uniform extension of high PM2.5 concentrations (100–200 µgm−3) is evident across the entire IGP region during November–December. The model also simulates a distinct pool of strongly enhanced PM2.5 (100–150 µgm−3) over the eastern IGP (Eastern Uttar Pradesh, and Bihar provinces of India) throughout the study period, despite lower emissions than the northwest and southeast IGP region. Such feature was earlier seen in the satellite-based observations and called as the ‘aerosol pool’ or ‘Bihar pollution pool’2,36. As the model is able to simulate the spatio-temporal characteristics of the PM2.5 distribution over this region, we analyze model results in greater detail to investigate the influences of various underlying processes and factors such as the meteorology, dynamics, and emissions.
Influences of meteorology and dynamics
Model simulated meteorological fields reveal strong spatio-temporal variability over the northern Indian subcontinent during October-December (Fig. 3). Higher temperature (28–33 °C) were simulated over the Punjab and Sindh regions of northwest IGP during October, whereas, temperature over other parts of the IGP were generally lower (25–28 °C), except over some parts of the eastern IGP (West Bengal, India and Bangladesh). Temperature shows a decline in November (20–25 °C) and December (15–20 °C) with nearly similar distribution over the western and central IGP, while the eastern IGP experiences slightly higher temperatures than other regions during November-December. These changes in the temperature are correlated with the changes in the boundary layer dynamics across the region, with higher boundary layer height (BLH) during October (>300 m), deeper BLH being confined to the western IGP over Pakistan (>450 m), while, shallower over the eastern IGP (<450 m). Changes in the meteorological and dynamical conditions lead to the suppression of boundary layer to as low as 200–250 m over the western IGP and 250–350 m over central and eastern IGP in November. The western IGP depicts BLH in the range of 200–250 m during December, while the BLH over the central IGP is 200–300 m and eastern IGP is 250–350 m. Model simulates weak north-westerlies (~1 ms−1) during October with corresponding ventilation coefficients41 (VC) of about 600–1000 m2s−1 over central IGP and a stronger atmospheric dispersion (>1000 m2s−1) surrounding India’s capital region of Delhi. A slight increase in the wind speed (1.5–2.5 ms−1) is revealed from October to November, besides an enhanced wind zone in the northern India (2.5–3.5 ms−1). With decrease in the wind speed from November to December, a significant shallowing of the boundary layer drastically reduces the ventilation to the minimum (VC: < 400 m2s−1) over most of the IGP.
WRF-Chem simulated wind vectors describe the atmospheric dynamic that facilitates the accumulation of fine aerosols over the IGP. Divergence over the northwestern IGP (region of strong biomass-burning emissions) and concurrently weaker winds over the central and eastern IGP allow the accumulation of fine aerosols from northwest, further downwind (central and eastern IGP), favoring its mixing with local anthropogenic influences. However, during December, reduction in BLH along with wind speeds brings down the ventilation coefficient over the IGP enhancing the PM2.5 dramatically. An analysis of vertical distribution of PM2.5 and wind speed (Fig. S3) shows that the vertical extent of higher aerosol concentrations is reduced during December with a concurrent reduction in the wind speed from surface to ~875 hPa. A weakening of winds towards December is seen aloft (Fig. S4) resulting in less efficient dispersion of aerosols accumulated in the planetary boundary layer over the IGP.
Meteorology and regional-scale dynamics is suggested to play the crucial role in the widespread PM2.5 build up across the IGP in December. Such changes in the meteorological conditions reducing drastically the ventilation across the IGP are typically seen every year (e.g. Fig. S5). Aerosol-radiation interactions are suggested to reduce the boundary layer evolution over this region by reducing the incoming solar radiation reaching surface, and such a feedback can further increase the aerosol concentrations near the surface42. The relative importance of the biomass-burning emissions compared to the anthropogenic emissions is analyzed next, based on two sensitivity simulations.
Influences of emissions: Biomass-burning versus anthropogenic
Biomass-burning emissions, in particular from the open burning of the agricultural residues, significantly enhance the aerosol loading over the IGP region during the study period, in addition to strong regional anthropogenic emissions18,22,24. Satellite-based observations (Fig. S6) clearly show strong biomass-burning activities over the northwest IGP during October, whereas, scattered fire activities with lower number of the total events are observed during December. Such distributions of the fire events are not specific for the year 2016 and that fire activities over this region are found to be similar in other years as well (fire activities for years 2015 and 2017 are also shown in the Fig. S6).
Figure 4 shows the percentage reductions in the PM2.5 concentrations due to switching off the biomass-burning and anthropogenic emissions as compared to the reference simulation. A significant reduction in the PM2.5 (by up to ~50–60%) is simulated over the northwest IGP region during October in absence of the biomass-burning emissions, e.g., over Haryana (40–50%), Punjab (50–60%), followed by Delhi National Capital Region (NCR) and the adjoining areas. Model reveals that biomass-burning emissions influenced the PM2.5 over larger areas of Punjab and the entire north-central Indian region during November. Analysis of back air trajectories demonstrates that the fire influences are efficiently transported to Delhi during October as well as November (Fig. S7). Whereas, airmasses may not pick maximum effects from these fires to the central and eastern IGP regions in October; but with the change in synoptic-scale winds in November fire effects are transported across the IGP; and aerosols are also dispersed to the central India and surrounding wider region (Figs. S7 and 4). Interestingly during December, when enhanced PM2.5 is most widespread, only minimal reduction (less than 10%) is simulated over IGP, without biomass-burning emissions. Nevertheless, fires are found to play major influences over the hills located in the north of IGP in India and parts of Nepal during December which is found to be in agreement with satellite-retrieved fire observations (Fig. S6). The effects of observed forest fires in the sub-Himalaya during December are most pronounced within that region itself with lesser effects being transported to the IGP (Figs. S7 and 4). By switching off the anthropogenic emissions, very large reductions (80–100%) in PM2.5 are simulated over nearly the entire IGP region, except the northwest IGP (40–70%) during the post-monsoon, complementing the strong influence of the biomass-burning emissions (noticed in the left panel of Fig. 4). Anthropogenic emissions are found to influence the entire IGP region by 90–100% during December.
Individual influences caused by switching off the biomass-burning and anthropogenic emissions on the PM2.5 levels are further combined to compute their relative effects (normalizing with the sum of their individual effects). Figure 5 shows the day-to-day variations in the relative effects from the regional biomass-burning and anthropogenic emissions on PM2.5 variations at selected locations (Delhi, Kanpur, and Varanasi) in the IGP from west to east. Model results exhibit strong effects of fire emissions on PM2.5 at Delhi (daily values up to 55.4%) during the mid-October to mid-November, attributed to the upwind biomass-burning (Figs. S6 and S7). Such influences from biomass burning are also considerable during this period at other locations in the IGP, nevertheless, effects are found to be lower (up to 36.3%) than that over Delhi. During this period of peak fire influences, the mean relative effects from the biomass burning are estimated to be 30.2%, 19.6%, and 9.4% at Delhi, Kanpur, and Varanasi respectively. However, the impact of biomass-burning decreases gradually from mid-November onwards at all the stations attaining minimal values (~5% and less) during December. Concurrently, the relative effect of anthropogenic emissions at all the stations increased to 95% or more during December.
Further, the influences of anthropogenic and biomass burning emissions on the chemical constituents of PM2.5, e.g., nitrate (NO3−), sulphate (SO4−2), ammonium (NH4+), organic matter (OM) and elemental carbon (EC) are examined. The differences in the concentrations of these species in absence of fire (Fig. S8) and anthropogenic emissions (Fig. S9) are computed as compared to those in the reference simulation. Reductions by up to ~20 µg/kg dry air (written as µg/kg hereafter) in NO3− are simulated in absence of the fire emissions over northwest IGP during October, which spread to larger region during November. Similarly, considerable reductions in the OM (by up to ~50 µg/kg) and EC (by ~10 µg/kg) are also simulated in absence of fires. These results clearly depict that composition of PM2.5 aerosols over the IGP region is significantly affected by the biomass burning emissions. The contributions of biomass burning emissions to NO3−, NH4+, OM, and EC are seen to drop dramatically during December. In contrast, the secondary inorganic species (NO3−, SO4−2, and NH4+) typically show major influences from the anthropogenic emissions (Fig. S9). Observations showing strong role of secondary sources at megacity Delhi during colder conditions43 also highlight the importance of meteorology and regional-scale air chemistry, besides emissions in governing the air quality over the IGP region.
Discussion
Our study shows that both the meteorological processes (winds, temperature, boundary layer mixing) as well as relative influences of emissions (biomass-burning versus anthropogenic) change drastically during the post-monsoon to winter transition, leading to the widespread enhancement in PM2.5 levels across the IGP region. In addition to previous studies, e.g.15,28,38, here the inter-dependence of PM2.5 distribution, meteorology, and atmospheric chemistry is explored further in impacting the regional air quality based on the high-resolution modeling of PM2.5 and constituents. The findings have important implications for the air quality in this densely populated region, which gets severely degraded during post-monsoon and winter seasons. Based on temporal changes in the dominance of biomass-burning and anthropogenic emissions, the study addresses their relative roles in building-up short-term dense haze over the IGP, elucidated in earlier ground-based studies18,20,22. Our current findings are also in line with the wintertime abundance of volatile organic compounds over IGP29 – which can act as aerosol precursors. In addition, the South Asian outflow to the oceanic regions is also strongest during the winter and that an aerosol build up here influences the atmospheric composition as well as climate over a larger region44,45. Therefore, the strategies, also taking into account the complex atmospheric dynamics besides the changes in emissions, are essential to reduce the aerosol build up and subsequently to mitigate their potential impacts on the air quality and climate, including the hydrological cycle. It should be noted that the focus of this study has been on the widespread and persistent enhancement which is clearly observed in the monthly mean distributions, although local sporadic emissions such as fireworks in the festive events have also been shown to increase the severity over Delhi22,46. Reductions in biomass-burning emissions would significantly improve the air quality of the IGP region during the post-monsoon19,21,22,37,47,48 but will have smaller effect in winter. A much more aggressive reduction of anthropogenic emissions over a wider area of IGP is required in the stagnant atmospheric conditions of winter to mitigate the PM2.5 enhancement caused by the meteorology and dynamics of the region. Strong dependence of surface PM2.5 on the meteorology simulated here using a regional model is in agreement with an earlier study based on global modeling over this region23. Among various anthropogenic sectors, residential energy usage is suggested to have larger effects on the annual mean PM2.5 in the region31. Future policies considering these would better complement the recent initiatives of the Indian government towards a reduction of household usage of biomass cook fuels49,50. Moreover, this study extends the possibilities of further evaluating health impacts due to changes in the individual aerosol constituents, originating from various biomass-burning and anthropogenic sources, besides the impact of overall PM2.511,31.
Our study describes the widespread PM2.5 pollution over northern India during post-monsoon to winter and highlights the changing contributions from biomass-burning vs anthropogenic sources during the transition period. We also highlight a west-east gradient in terms of dominance of sources, with western cities dominated by biomass-burning sources and eastern cities dominated by anthropogenic sources during post-monsoon; this gradient weakens with the arrival of winter when the entire region from west to east is dominated by anthropogenic sources. We argue that a generic emission reduction policy will not yield desired results, and advocate for a season-based source-focused mitigation policy involving a reduction of anthropogenic emissions ubiquitously taking unfavorable regional dynamics into account, to improve the air quality and minimize the climate impacts.
Methods
Weather Research and Forecasting model coupled with Chemistry (WRF-Chem)33,51, an online regional model version 3.8.1 has been used to simulate the meteorology and spatio-temporal distribution of PM2.5 over the northern Indian region. Model domain is centered at 76°E, 29°N and has 240 grid points in the east-west direction and 147 grid points in the north-south direction, along with 51 vertical levels. The simulations are conducted at a horizontal resolution of 12 km × 12 km. Anthropogenic emissions are included from the Southeast Asia Composition, Cloud, Climate Coupling Regional Study (SEAC4RS) for gas-phase species, except the NH3, which together with aerosols are included from the Hemispheric Transport of Air Pollution (HTAP) inventory52. The spatial resolution of these emissions is 0.1° × 0.1o. Biomass-burning emissions are based on the NCAR Fire Inventory (FINN)53 and biogenic emissions are calculated online using the Model of Emissions of Gases and Aerosols from Nature (MEGAN)54. Initial and lateral boundary conditions for the meteorological and chemical fields are included from the Era Interim and MOZART-455 respectively. To limit the errors in simulated meteorology and dynamics, model is nudged to the Era Interim at all the vertical levels. Regional Acid Deposition Model - 2nd generation (RADM2)56 coupled with the Modal Aerosol Dynamics Model for Europe/Secondary Organic Aerosol Model (MADE/SORGAM) aerosol module57,58 has been used to simulate the gas-phase chemistry and aerosols in the model. Application of this chemical mechanism showed small bias (~7%) in PM2.5 simulation over Europe59,60. The schemes opted to parameterize different processes in the model (Table S1) are based on earlier studies over this region using the WRF-Chem model32,34,35. A reduction in dust transport to 25% from the domain boundaries has been made resulting in relatively better agreement with the observations of AOD at 0.5 µm, nevertheless, no large changes were seen in the PM2.5 at the surface. Besides the reference simulation, two sensitivity simulations have been conducted (Table 1) by switching off the anthropogenic and biomass-burning emissions, one by one.
The observations of Aerosol optical depth (AOD) have been obtained from the MODIS (Moderate resolution Imaging Spectroradiometer) instrument onboard Aqua satellite, and MISR (Multi-angle Imaging Spectroradiometer) instrument onboard Terra satellite. Ground-based observations of AOD from the AERONET (Aerosol Robotic Network; https://aeronet.gsfc.nasa.gov/)61, at two stations (Kanpur and Jaipur) have been used. Fire count observations are obtained from the MODIS instrument to assess the biomass burning events. Ground-based observations of PM2.5 (µg m−3) at different stations in the IGP region (coordinates of stations provided in the Table S2) are obtained from the Central Pollution Control Board (CPCB), India (https://app.cpcbccr.com/ccr/#/caaqm-dashboard-all/caaqm-landing/data). The hourly datasets needed a screening for removal of fill values and also for very high values at few hours as discussed elsewhere23, performed before computing the daily average values. The HYSPLIT (HYbrid Single Particle Lagrangian Integrated Trajectory) model62 with meteorological inputs from the Global Data Assimilation System has been used to simulate the backward airmass trajectories.
References
Babu, S. S. et al. Trends in aerosol optical depth over Indian region: Potential causes and impact indicators. Journal of Geophysical Research-Atmospheres 118, 11794–11806, https://doi.org/10.1002/2013jd020507 (2013).
Di Girolamo, L. et al. Analysis of Multi-angle Imaging SpectroRadiometer (MISR) aerosol optical depths over greater India during winter 2001-2004. Geophysical Research Letters 31, https://doi.org/10.1029/2004gl021273 (2004).
Jethva, H., Satheesh, S. K. & Srinivasan, J. Seasonal variability of aerosols over the Indo-Gangetic basin. Journal of Geophysical Research-Atmospheres 110, https://doi.org/10.1029/2005jd005938 (2005).
Moorthy, K. K., Babu, S. S., Manoj, M. R. & Satheesh, S. K. Buildup of aerosols over the Indian Region. Geophysical Research Letters 40, 1011–1014, https://doi.org/10.1002/grl.50165 (2013).
Lelieveld, J. et al. The South Asian monsoon-pollution pump and purifier. Science 361, 270, https://doi.org/10.1126/science.aar2501 (2018).
Ojha, N. et al. Variabilities in ozone at a semi-urban site in the Indo-Gangetic Plain region: Association with the meteorology and regional processes. Journal of Geophysical Research-Atmospheres 117, https://doi.org/10.1029/2012jd017716 (2012).
Ramanathan, V. et al. Atmospheric brown clouds: Impacts on South Asian climate and hydrological cycle. Proceedings of the National Academy of Sciences of the United States of America 102, 5326–5333, https://doi.org/10.1073/pnas.0500656102 (2005).
Tripathi, S. N., Dey, S., Tare, V. & Satheesh, S. K. Aerosol black carbon radiative forcing at an industrial city in northern India. Geophysical Research Letters 32, https://doi.org/10.1029/2005gl022515 (2005).
Conibear, L., Butt, E. W., Knote, C., Arnold, S. R. & Spracklen, D. V. Residential energy use emissions dominate health impacts from exposure to ambient particulate matter in India. Nature Communications 9, https://doi.org/10.1038/s41467-018-02986-7 (2018).
David, L. M. et al. Premature Mortality Due to PM2.5 Over India: Effect of Atmospheric Transport and Anthropogenic Emissions. GeoHealth 3, 2–10, https://doi.org/10.1029/2018gh000169 (2019).
Ghude, S. D. et al. Premature mortality in India due to PM2.5 and ozone exposure. Geophysical Research Letters 43, 4650–4658, https://doi.org/10.1002/2016gl068949 (2016).
Lelieveld, J., Evans, J. S., Fnais, M., Giannadaki, D. & Pozzer, A. The contribution of outdoor air pollution sources to premature mortality on a global scale. Nature 525, 367–+, https://doi.org/10.1038/nature15371 (2015).
Murari, V., Kumar, M., Barman, S. C. & Banerjee, T. Temporal variability of MODIS aerosol optical depth and chemical characterization of airborne particulates in Varanasi, India. Environmental Science and Pollution Research 22, 1329–1343, https://doi.org/10.1007/s11356-014-3418-2 (2015).
Nair, V. S. et al. Wintertime aerosol characteristics over the Indo-Gangetic Plain (IGP): Impacts of local boundary layer processes and long-range transport. Journal of Geophysical Research-Atmospheres 112, https://doi.org/10.1029/2006jd008099 (2007).
Ram, K., Sarin, M. M. & Tripathi, S. N. Temporal Trends in Atmospheric PM2.5, PM10, Elemental Carbon, Organic Carbon, Water-Soluble Organic Carbon, and Optical Properties: Impact of Biomass Burning Emissions in The Indo-Gangetic Plain. Environmental Science & Technology 46, 686–695, https://doi.org/10.1021/es202857w (2012).
Tiwari, S. et al. Diurnal and seasonal variations of black carbon and PM2.5 over New Delhi, India: Influence of meteorology. Atmospheric Research 125, 50–62, https://doi.org/10.1016/j.atmosres.2013.01.011 (2013).
Bikkina, S. et al. Air quality in megacity Delhi affected by countryside biomass burning. Nature Sustainability 2, 200–205, https://doi.org/10.1038/s41893-019-0219-0 (2019).
Guttikunda, S. K., Nishadh, K. A. & Jawahar, P. Air pollution knowledge assessments (APnA) for 20 Indian cities. Urban Climate 27, 124–141, https://doi.org/10.1016/j.uclim.2018.11.005 (2019).
Jethva, H. et al. Connecting Crop Productivity, Residue Fires, and Air Quality over Northern India. Scientific Reports 9, https://doi.org/10.1038/s41598-019-52799-x (2019).
Liu, T. J. et al. Seasonal impact of regional outdoor biomass burning on air pollution in three Indian cities: Delhi, Bengaluru, and Pune. Atmospheric Environment 172, 83–92, https://doi.org/10.1016/j.atmosenv.2017.10.024 (2018).
Sarkar, S., Singh, R. P. & Chauhan, A. Crop Residue Burning in Northern India: Increasing Threat to Greater India. Journal of Geophysical Research-Atmospheres 123, 6920–6934, https://doi.org/10.1029/2018jd028428 (2018).
Sawlani, R. et al. The severe Delhi SMOG of 2016: A case of delayed crop residue burning, coincident firecracker emissions, and atypical meteorology. Atmospheric Pollution Research 10, 868–879, https://doi.org/10.1016/j.apr.2018.12.015 (2019).
Schnell, J. L. et al. Exploring the relationship between surface PM2.5 and meteorology in Northern India. Atmospheric Chemistry and Physics 18, 10157–10175, https://doi.org/10.5194/acp-18-10157-2018 (2018).
Sharma, S. K. et al. Source Apportionment of PM2.5 in Delhi, India Using PMF Model. Bulletin of Environmental Contamination and Toxicology 97, 286–293, https://doi.org/10.1007/s00128-016-1836-1 (2016).
Chandra, B. P. et al. Odd-even traffic rule implementation during winter 2016 in Delhi did not reduce traffic emissions of VOCs, carbon dioxide, methane and carbon monoxide. Current Science 114, 1318–1325, https://doi.org/10.18520/cs/v114/i06/1318-1325 (2018).
Sharma, S. K. et al. Study on Ambient Air Quality of Megacity Delhi, India During Odd-Even Strategy. Mapan-Journal of Metrology Society of India 32, 155–165, https://doi.org/10.1007/s12647-016-0201-5 (2017).
Moorthy, K. K., Satheesh, S. K. & Kotamarthi, V. R. Evolution of aerosol research in India and the RAWEX-GVAX: an overview. Current Science 111, 53–75, https://doi.org/10.18520/cs/v111/i1/53-75 (2016).
Bran, S. H. & Srivastava, R. Investigation of PM2.5 mass concentration over India using a regional climate model. Environmental Pollution 224, 484–493, https://doi.org/10.1016/j.envpol.2017.02.030 (2017).
Chutia, L. et al. Distribution of volatile organic compounds over Indian subcontinent during winter: WRF-chem simulation versus observations. Environmental Pollution 252, 256–269, https://doi.org/10.1016/j.envpol.2019.05.097 (2019).
Kumar, R. et al. What controls the seasonal cycle of black carbon aerosols in India? Journal of Geophysical Research-Atmospheres 120, 7788–7812, https://doi.org/10.1002/2015jd023298 (2015).
Reddington, C. L. et al. Exploring the impacts of anthropogenic emission sectors on PM2.5 and human health in South and East Asia. Atmospheric Chemistry and Physics 19, 11887–11910, https://doi.org/10.5194/acp-19-11887-2019 (2019).
Sharma, A. et al. WRF-Chem simulated surface ozone over south Asia during the pre-monsoon: effects of emission inventories and chemical mechanisms. Atmospheric Chemistry and Physics 17, 14393–14413, https://doi.org/10.5194/acp-17-14393-2017 (2017).
Grell, G. A. et al. Fully coupled “online” chemistry within the WRF model. Atmospheric Environment 39, 6957–6975, https://doi.org/10.1016/j.atmosenv.2005.04.027 (2005).
Kumar, R., Naja, M., Pfister, G. G., Barth, M. C. & Brasseur, G. P. Simulations over South Asia using the Weather Research and Forecasting model with Chemistry (WRF-Chem): set-up and meteorological evaluation. Geoscientific Model Development 5, 321–343, https://doi.org/10.5194/gmd-5-321-2012 (2012).
Kumar, R. et al. Simulations over South Asia using the Weather Research and Forecasting model with Chemistry (WRF-Chem): chemistry evaluation and initial results. Geoscientific Model Development 5, 619–648, https://doi.org/10.5194/gmd-5-619-2012 (2012).
Kumar, M. et al. Long-term aerosol climatology over Indo-Gangetic Plain: Trend, prediction and potential source fields. Atmospheric Environment 180, 37–50, https://doi.org/10.1016/j.atmosenv.2018.02.027 (2018).
Rajput, P., Sarin, M., Sharma, D. & Singh, D. Characteristics and emission budget of carbonaceous species from post-harvest agricultural-waste burning in source region of the Indo-Gangetic Plain. Tellus Series B-Chemical and Physical Meteorology 66, https://doi.org/10.3402/tellusb.v66.21026 (2014).
Sen, A. et al. Variations in particulate matter over Indo-Gangetic Plains and Indo-Himalayan Range during four field campaigns in winter monsoon and summer monsoon: Role of pollution pathways. Atmospheric Environment 154, 200–224, https://doi.org/10.1016/j.atmosenv.2016.12.054 (2017).
Ramachandran, S., Kedia, S. & Sheel, V. Spatiotemporal characteristics of aerosols in India: Observations and model simulations. Atmospheric Environment 116, 225–244, https://doi.org/10.1016/j.atmosenv.2015.06.015 (2015).
Singh, N. et al. Boundary layer evolution over the central Himalayas from radio wind profiler and model simulations. Atmospheric Chemistry and Physics 16, 10559–10572, https://doi.org/10.5194/acp-16-10559-2016 (2016).
Zhang, Y., Liu, Z. H., Lv, X. T., Zhang, Y. & Qian, J. Characteristics of the Transport of a Typical Pollution Event in the Chengdu Area Based on Remote Sensing Data and Numerical Simulations. Atmosphere 7, https://doi.org/10.3390/atmos7100127 (2016).
Bharali, C., Nair, V. S., Chutia, L. & Babu, S. S. Modeling of the Effects of Wintertime Aerosols on Boundary Layer Properties Over the Indo Gangetic Plain. Journal of Geophysical Research-Atmospheres 124, 4141–4157, https://doi.org/10.1029/2018jd029758 (2019).
Gani, S. et al. Submicron aerosol composition in the world’s most polluted megacity: the Delhi Aerosol Supersite study. Atmospheric Chemistry and Physics 19, 6843–6859, https://doi.org/10.5194/acp-19-6843-2019 (2019).
Lawrence, M. G. & Lelieveld, J. Atmospheric pollutant outflow from southern Asia: a review. Atmospheric Chemistry and Physics 10, 11017–11096, https://doi.org/10.5194/acp-10-11017-2010 (2010).
Lelieveld, J. et al. The Indian Ocean Experiment: Widespread air pollution from South and Southeast Asia. Science 291, 1031–1036, https://doi.org/10.1126/science.1057103 (2001).
Mukherjee, T. et al. Increasing Potential for Air Pollution over Megacity New Delhi: A Study Based on 2016 Diwali Episode. Aerosol and Air Quality Research 18, 2510–2518, https://doi.org/10.4209/aaqr.2017.11.0440 (2018).
Cusworth, D. H. et al. Quantifying the influence of agricultural fires in northwest India on urban air pollution in Delhi, India. Environmental Research Letters 13, https://doi.org/10.1088/1748-9326/aab303 (2018).
Rastogi, N., Singh, A., Sarin, M. M. & Singh, D. Temporal variability of primary and secondary aerosols over northern India: Impact of biomass burning emissions. Atmospheric Environment 125, 396–403, https://doi.org/10.1016/j.atmosenv.2015.06.010 (2016).
Chowdhury, S. et al. Indian annual ambient air quality standard is achievable by completely mitigating emissions from household sources. Proceedings of the National Academy of Sciences of the United States of America 116, 10711–10716, https://doi.org/10.1073/pnas.1900888116 (2019).
Venkataraman, C. et al. Source influence on emission pathways and ambient PM2:5 pollution over India (2015-2050). Atmospheric Chemistry and Physics 18, 8017–8039, https://doi.org/10.5194/acp-18-8017-2018 (2018).
Fast, J. D. et al. Evolution of ozone, particulates, and aerosol direct radiative forcing in the vicinity of Houston using a fully coupled meteorology-chemistry-aerosol model. Journal of Geophysical Research-Atmospheres 111, https://doi.org/10.1029/2005jd006721 (2006).
Janssens-Maenhout, G. et al. HTAP_v2.2: a mosaic of regional and global emission grid maps for 2008 and 2010 to study hemispheric transport of air pollution. Atmospheric Chemistry and Physics 15, 11411–11432, https://doi.org/10.5194/acp-15-11411-2015 (2015).
Wiedinmyer, C. et al. The Fire INventory from NCAR (FINN): a high resolution global model to estimate the emissions from open burning. Geoscientific Model Development 4, 625–641, https://doi.org/10.5194/gmd-4-625-2011 (2011).
Guenther, A. et al. Estimates of global terrestrial isoprene emissions using MEGAN (Model of Emissions of Gases and Aerosols from Nature). Atmospheric Chemistry and Physics 6, 3181–3210, https://doi.org/10.5194/acp-6-3181-2006 (2006).
Emmons, L. K. et al. Description and evaluation of the Model for Ozone and Related chemical Tracers, version 4 (MOZART-4). Geoscientific Model Development 3, 43–67, https://doi.org/10.5194/gmd-3-43-2010 (2010).
Stockwell, W. R., Middleton, P., Chang, J. S. & Tang, X. Y. The 2nd Generation Regional Acid Deposition Model Chemical Mechanism For Regional Air-Quality Modeling. Journal of Geophysical Research-Atmospheres 95, 16343–16367, https://doi.org/10.1029/JD095iD10p16343 (1990).
Ackermann, I. J. et al. Modal aerosol dynamics model for Europe: Development and first applications. Atmospheric Environment 32, 2981–2999, https://doi.org/10.1016/s1352-2310(98)00006-5 (1998).
Schell, B., Ackermann, I. J., Hass, H., Binkowski, F. S. & Ebel, A. Modeling the formation of secondary organic aerosol within a comprehensive air quality model system. Journal of Geophysical Research-Atmospheres 106, 28275–28293, https://doi.org/10.1029/2001jd000384 (2001).
Georgiou, G. K. et al. Air quality modelling in the summer over the eastern Mediterranean using WRF-Chem: chemistry and aerosol mechanism intercomparison. Atmospheric Chemistry and Physics 18, 1555–1571, https://doi.org/10.5194/acp-18-1555-2018 (2018).
Tuccella, P. et al. Modeling of gas and aerosol with WRF/Chem over Europe: Evaluation and sensitivity study. Journal of Geophysical Research-Atmospheres 117, https://doi.org/10.1029/2011jd016302 (2012).
Holben, B. N. et al. AERONET - A federated instrument network and data archive for aerosol characterization. Remote Sensing of Environment 66, 1–16, https://doi.org/10.1016/s0034-4257(98)00031-5 (1998).
Rolph, G., Stein, A. & Stunder, B. Real-time Environmental Applications and Display sYstem: READY. Environmental Modelling & Software 95, 210–228, https://doi.org/10.1016/j.envsoft.2017.06.025 (2017).
Acknowledgements
We acknowledge use of various WRF-Chem preprocessor tools provided by the Atmospheric Chemistry Observations and Modeling Lab (ACOM) of NCAR. Initial and boundary condition data for meteorological fields were obtained from the ECMWF. We acknowledge use of MOZART4 global model output available at https://www.acom.ucar.edu/wrf-chem/mozart.shtml. HTAP anthropogenic emissions were obtained from http://edgar.jrc.ec.europa.eu/htap_v2/index.php?SECURE=123. Authors are grateful to Yafang Cheng (MPI-C) for providing SEAC4RS emission. Use of observational data from the Central Pollution Control Board (CPCB), India is acknowledged. MODIS fire observations are obtained from ftp://fuoco.geog.umd.edu/modis/C6/mcd14ml/. We acknowledge data centers at NASA for providing the MODIS, MISR and AERONET data. We are also grateful to the Principal Investigators and staff for establishing and maintaining the AERONET sites used in this study. We acknowledge NOAA Air Resources Laboratory (ARL) for the HYSPLIT model (https://www.arl.noaa.gov/hysplit/hysplit/) used in the study. Authors acknowledge the use of Cobra HPC, which is a part of the Max Planck Computing and Data Facility (MPCDF) (https://www.mpcdf.mpg.de/), and Vikram HPC at the PRL, Ahmedabad. N. Ojha acknowledges valuable support and encouragement from D. Pallamraju and Anil Bhardwaj. Constructive comments and suggestions from the anonymous reviewers are gratefully acknowledged. SSG acknowledge the partial financial support from DST - IGCS at IIT Madras (DST; sanction number: DST/CCP/CoE/141/2018C).
Author information
Authors and Affiliations
Contributions
N.O., S.S.G. and A.P. conceived the idea. N.O. and A.S. performed the model simulations. N.O., A.S., and I.G. performed the model, AERONET, and satellite data analysis with inputs from T.U.A., M.K., S.S., A.P., N.S. and S.S.G. N.O. wrote the manuscript with input from all the co-authors
Corresponding authors
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Ojha, N., Sharma, A., Kumar, M. et al. On the widespread enhancement in fine particulate matter across the Indo-Gangetic Plain towards winter. Sci Rep 10, 5862 (2020). https://doi.org/10.1038/s41598-020-62710-8
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-020-62710-8