Blog

Jan 07, 2024

Thermodynamical framework for effective mitigation of high aerosol loading in the Indo

Scientific Reports volume 13, Article number: 13667 (2023) Cite this article 357 Accesses 2 Altmetric Metrics details The Indo-Gangetic Plain (IGP) experiences severe air pollution every winter, with

Scientific Reports volume 13, Article number: 13667 (2023) Cite this article

357 Accesses

2 Altmetric

Metrics details

The Indo-Gangetic Plain (IGP) experiences severe air pollution every winter, with ammonium chloride and ammonium nitrate as the major inorganic fractions of fine aerosols. Many past attempts to tackle air pollution in the IGP were inadequate, as they targeted a subset of the primary pollutants in an environment where the majority of the particulate matter burden is secondary in nature. Here, we provide new mechanistic insight into aerosol mitigation by integrating the ISORROPIA-II thermodynamical model with high-resolution simultaneous measurements of precursor gases and aerosols. A mathematical framework is explored to investigate the complex interaction between hydrochloric acid (HCl), nitrogen oxides (NOx), ammonia (NH3), and aerosol liquid water content (ALWC). Aerosol acidity (pH) and ALWC emerge as governing factors that modulate the gas-to-particle phase partitioning and mass loading of fine aerosols. Six "sensitivity regimes" were defined, where PM1 and PM2.5 fall in the "HCl and HNO3 sensitive regime", emphasizing that HCl and HNO3 reductions would be the most effective pathway for aerosol mitigation in the IGP, which is ammonia-rich during winter. This study provides evidence that precursor abatement for aerosol mitigation should not be based on their descending mass concentrations but instead on their sensitivity to high aerosol loading.

High aerosol loading is a significant cause of millions of premature deaths around the world, and mitigating air pollution is a major concern for researchers over the globe1,2. However, to regulate or mitigate high aerosol loading, it needs to be measured, monitored, and investigated thoroughly. Despite extensive efforts, the understanding of the physical, chemical, and thermodynamical properties of atmospheric constituents has not yet reached a point where aerosol mitigation can be done precisely and optimally3,4,5,6. There is a growing need to better understand the aerosol properties as they severely affect the ecosystem, human health, and the environment.

The Indo-Gangetic Plain (IGP) is one of the most polluted regions in the world7,8,9,10. The mass loading of PM1 and PM2.5 often exceeds 400–600 µg m−3 for the short-term (few hours), particularly during traffic rush hours and nighttime during winter11,12. Studies have shown that the organic fraction is generally responsible for more than 50% of PM1 mass globally, and the inorganic fraction of fine aerosol is composed of sulfate, nitrate, and ammonium (SNA)13,14,15,16. But in contrast, few studies have also shown that more than half of total aerosol loading is inorganic in nature during peak pollution episodes (Table S1), pointing out the necessity of investigating the role of physicochemical and thermodynamical properties of atmospheric constituents in contrasting atmospheric conditions. For instance, Gani et al.17 have shown that during peak pollution periods of January 2018, the inorganic fractions contributed nearly 60% of the total PM1 mass loading in Delhi. These inorganic aerosols are primarily composed of chloride, sulfate, nitrate, and ammonium (CSNA)18,19.

The observed hourly chloride concentrations many times exceed 100 µg m−3, considered to be among the highest reported anywhere in the world20. The gas-phase ammonia (NH3) is also found to be very high, which significantly affects the secondary aerosol formation in winter. Several recent studies have investigated the sensitivity of aerosols to the reduction in precursor gases in China, the USA, and Europe4, but studies conducted in India have not yet been performed. Extensive efforts are needed to improve our scientific understanding of the effective regulation of aerosol loading in the IGP.

In this study, we aim to provide new mechanistic insights into aerosol formation and mitigation by integrating the ISORROPIA-II thermodynamic equilibrium model with the dataset of precursor gases (HCl, HNO3, and NH3) and inorganic constituents (Cl−, NO3−, SO42−, Na+, NH4+, K+, Ca2+, and Mg2+) of PM1 and PM2.5 acquired using the first deployment of the MARGA-2S instrument in the IGP (Supplementary text S1). We explore a mathematical framework to investigate the sensitivity of gas-to-particle partitioning of aerosols to different parameters like gaseous precursor concentrations (HCl, HNO3, and NH3), pH, and ALWC, using sigmoidal curves and "sensitivity regimes" of aerosols. To our knowledge, this is the first attempt to investigate the thermodynamic control of aerosols over the Indian region, and the "thermodynamical roadmap" advanced in this study could provide effective and targeted mitigation strategies in the IGP.

The ISORROPIA-II model was run in the "forward" mode using the hourly measured total (gas + particle) species concentrations instead of only particle-phase concentrations as input, as it is considered to be more accurate21,22. As shown in Figure S1, the predicted PM1 NH4+ (r = 0.97), NO3− (r = 0.93), Cl− (r = 0.98), PM2.5 NH4+ (r = 0.87), NO3− (r = 0.98), Cl− (r = 0.99), and gas-phase NH3 (r = 0.93) showed very good correlation with the measured concentrations, confirming the reliability of the ISORROPIA-II model simulations at a range of temperatures (278–298 K).

The data points are color-coded by ambient temperatures to investigate how temperature variability can alter the ISORROPIA-II predicted gas and particle-phase concentrations. The pH and ALWC of PM1 and PM2.5 were estimated during 19 December 2017 to 10 February 2018. The predicted pH of PM1 varied between 2.2 to 5.6, and the mean PM1 pH (average ± SD) was 4.5 ± 0.5. The PM2.5 pH ranged from 2.5 to 6.5, with a mean value of 4.6 ± 0.5. The predicted pH of PM2.5 was similar to the measured PM2.5 pH of 4.6 ± 0.5 over Delhi, estimated during the Winter Fog EXperiment (WiFEX) campaign period of 2015–16, underscoring the reliability and accuracy of the ISORROPIA-II model simulation20,23.

In spite of the high reliability of the ISORROPIA-II model simulations (Figure S1), it is essential to investigate the uncertainties in observations while using them as inputs to model simulations. There could be many sources of errors in the measurements, like gas-aerosol collection efficiency, internal standard variability, and reproducibility of the chromatogram peak integration in the MARGA-2S. These errors sometimes could propagate to cause significant uncertainties in the model simulations.

It can be seen that, at almost all times during the study, the concentrations of the measured parameters were much higher than the range of errors, causing them to be unlikely the primary source of bias in the analysis. It is crucial to note that the analytic error in monitoring individual species cannot be evaluated in the absence of specifically designed quality-control experiments. However, to ensure the quality of the data, observational hours were excluded from the study when measurements did not meet the quality assurance and quality check (QA/QC) criteria19. As a result, only 1100 h of correct observational datasets were used for the model simulations, and other datasets are not included in the analysis of this paper.

The major errors inherent with thermodynamic equilibrium modeling are derived from two assumptions, i.e., internal mixing and gas-aerosol equilibria. Among different factors that contribute to the uncertainty in the model simulations, input parameters like gas concentrations, aerosol composition, temperature, and relative humidity are considered to be the most crucial. Uncertainties in these input parameters could cause biases in predicted parameters like aerosol acidity (pH), aerosol liquid water content (ALWC), and gas-aerosol phase partitioning (ε).

To quantitively investigate this, we carried out simulations with different sets of gas and particle-phase inputs. The analysis has also considered the limit of detection of all gas and aerosol species by the MARGA instrument (Table S3). The statistical variability was evaluated by checking the mean bias (MB), normalized mean bias (NMB), and root mean square error (RMSE) following Eqs. (S7–S9). As many previous studies in India considered fixed concentrations of HCl + Cl− (25 µg m−3) as model input due to a lack of gas-phase HCl observations in their studies24,25, we also fixed the input concentrations to 25 µg m−3. To make the analysis more robust, we fixed the concentrations of all three gas + aerosol constituents, i.e., HCl + Cl−, HNO3 + NO3−, and NH3 + NH4+. The parameters in Table S4 illustrate the uncertainty in aerosol acidity (pH) and aerosol liquid water content (ALWC) due to these fixed input concentrations.

As ALWC governs the presence of ions and the water uptake of different hygroscopic inorganic constituents, it is crucial to investigate the uncertainty in ALWC due to different assumptions in model input, especially in a region like Delhi, where ALWC is very high during winter, and the heterogeneous and multiple-phase reactions govern during peak-pollution events. The fixed concentration of HCl + Cl− = 25 µg m−3 resulted in NMB of − 8.34% and − 26.61% for PM1 and PM2.5 ALWC, respectively. The RMSE is 149.69 µg m−3 and 264.98 µg m−3 for fixed total chloride input concentrations. The RMSE in ALWC due to fixed total nitrate (HNO3 + NO3− = 25 µg m−3) and total ammonium (NH3 + NH4+  = 25 µg m−3) concentrations are 46.8 µg m−3, 62.49 µg m−3 and 93.70 µg m−3, 307.81 µg m−3 for PM1 and PM2.5 respectively. These high uncertainties could significantly alter the predicted gas-aerosol phase partitioning and particle-phase loading in the atmosphere.

The aerosol acidity (pH) of PM1 and PM2.5 also shows significant bias due to assumed fixed input concentrations (Table S4). The RMSE of PM1 and PM2.5 pH due to fixed total chloride input concentration is 0.42 and 0.55 respectively, which can significantly alter the phase-partitioning and particle-phase chloride loading.

Furthermore, we have extended our efforts by conducting an uncertainty analysis in the model simulations to account for the inherent measurement uncertainties associated with the MARGA-2S trace gas and aerosol species measurements. This undertaking involves the consideration of three distinct input values, encompassing the average and 95% confidence interval upper and lower limits, for all species under investigation in the "forward" mode model simulations.

Of notable importance are the primary input parameters (95% confidence intervals indicated in parentheses): TCl = 21.78 (20.59–22.98) µg m−3, TNO3 = 19.74 (19.17–20.31) µg m−3, TNH4 = 43.43 (42.28–45.58) µg m−3, SO42− = 11.36 (10.99–11.73) µg m−3 for PM1, and TCl = 39.71 (37.65–41.77) µg m−3, TNO3 = 31.65 (30.80–32.50) µg m−3, TNH4 = 58.02 (56.62–59.42) µg m−3, and SO42− = 20.34 (19.52–21.15) µg m−3 for PM2.5. We have included the uncertainties associated with the model-derived aerosol liquid water content (ALWC) and aerosol acidity (pH) values for both PM1 and PM2.5 in Table S4. These results collectively emphasize the robustness and reliability of the ISORROPIA-II model simulations.

In regions like Europe, where HCl and Cl− concentrations are very low, the assumption of fixed HCl + Cl− may not have a grave impact, but in a region like Delhi, where observed HCl + Cl− concentrations varied between 2.11–112.21 µg m−3 and 2.31–186.31 µg m−3 for PM1 and PM2.5 respectively, the consideration of actual ambient concentrations in ISORROPIA-II model simulations seems to be of utmost importance. The fixed concentration of 25 µg m−3 cannot accurately represent real-world conditions. Previous studies have also emphasized that when gas-phase data are not available, running ISORROPIA-II in the forward mode, with only aerosol concentrations as input, results in the repartitioning of ammonia in the model, causing bias in the predicted pH and ALWC3,21. These sensitivity analyses with fixed input concentrations of total chloride, nitrate, and ammonium illustrate the necessity of accounting for the actual gas and aerosol species concentrations instead of fixed input concentrations, especially in the IGP during winter.

As the high ammonium chloride concentration is a special characteristic of air pollution in India, and particle-phase chloride and nitrates are the major contributors to the inorganic fraction of the particulate matter burden in Delhi, it is crucial to explore their variability in prevailing meteorological conditions17,19,24. But under wintertime ambient meteorological conditions, pure hydrochloric acid, nitric acid, or water particles are difficult to form due to the higher vapor pressure of HCl and HNO3. Instead, ammonium nitrate (NH4NO3) and ammonium chloride (NH4Cl) form with significantly lower vapor pressure and has longer atmospheric residence time compared to the corresponding precursor gases like HNO3, HCl, and NH326. The average chloride partitioning ratio ε(Cl−) of PM1 and PM2.5 was 0.93 ± 0.09 and 0.96 ± 0.07 respectively (Table S2), implying the dominant presence of chloride in the particle phase during winter. The ε(Cl−) was 0.4 at RH ≤ 50%, which sharply increased to 0.95 at RH ≥ 80%, showing the enhanced phase-partitioning to the particulate chloride phase. The increased phase-partitioning of total chloride (HCl + Cl−) in highly humid conditions can further promote chloride formation caused by increased ALWC. The enhanced ALWC increases pH by dilution, further increasing total chloride partitioning and significantly increasing PM1 and PM2.5 chloride in a positive feedback loop22,27.

It has been shown that aerosols with higher chloride mass fraction uptake more water than those with less chloride mass fraction under prevailing RH conditions24,28. This is due to the co-condensation of HCl, NH3, and water, as the gas-phase HCl gets dissolved in aerosol water, dissociates, and then equilibrates with ammonia to form ammonium chloride, stabilizing chloride in the particle phase26,29. This particle-phase chloride can absorb even more water from the air, leading to enhanced growth of aerosols into fog droplets during winter and augmenting the particle mass loading24. These results demonstrate the role of ALWC in the phase-partitioning of PM1 and PM2.5 aerosols, which needs to be thoroughly investigated for a deeper understanding of the complex thermodynamical control of high aerosol loading.

The average ε(NO3−) of PM1 and PM2.5 was 0.83 ± 0.11 and 0.89 ± 0.08 respectively, showing the dominance of particle-phase nitrate over gas-phase HNO3 during winter. The high particle-phase nitrate concentrations significantly impact the total aerosol loading as the presence of more ammonium nitrate reduces the deliquescence relative humidity (DRH) for single-component aerosols and mutual deliquescence relative humidity (MDRH) for multicomponent aerosols resulting in the formation of more secondary aerosols in polluted conditions30,31,32. DRH and MDRH are the critical values of relative humidity, above which the water uptake of single component and multicomponent aerosols drastically increases, leading to increased surface area for heterogeneous reactions, enhanced aqueous phase reaction rates, and uptake coefficients of trace gases31,33. This, in turn, governs the atmospheric lifetime of both gas and aerosol phases.

To comprehend the importance of DRH or MDRH, we can see that sea salt particles contain highly hygroscopic salt magnesium chloride (MgCl2) having DRH of ~ 32%, due to which sea salt particles deliquesce at very low RH. This makes the marine atmosphere hazier than continental air at the same RH and ambient temperature (T). So, the presence of salts with less DRH or MDRH will make the atmosphere hazier than those with less hygroscopic salts like ammonium chloride and ammonium nitrate (38). As the efflorescence RH (ERH) of these aerosols are generally low and rarely reached by the ambient RH during winter, these species primarily remain in the aqueous or metastable phase34.

It is to be noted that, the solubility of ammonia is fairly weak, highly temperature-dependent, and greatly influenced by the aerosol pH and effective Henry's law constant. The average ε(NH4+) for PM1 and PM2.5 was 0.42 ± 0.17 and 0.55 ± 0.15 respectively, attributed to the substantial excess ammonia left in the gas-phase. Due to this, several studies have highlighted that pollution control strategies in India should be directed at reducing high NH3 concentrations25,35. Here, in this study, we did not investigate the effect of pH and ALWC on the phase partitioning of sulfate ε(SO42−). This is because the sulfate was considered to be majorly present in the ammonium sulfate form, which has a very low volatility in the wintertime ambient temperature range (38). So, this study focuses on investigating the phase partitioning of semi-volatile species like ammonium chloride and ammonium nitrate during winter.

Due to the significantly different particulate matter (PM) burden, PM chemical composition, particle number size distribution, and meteorological conditions between Delhi and other parts of the globe, these phase-partitioning ratios in Delhi are in significant contrast24,28. For instance, ε(NO3−) was reported to be 0.26 ± 0.15 and 0.39 ± 0.16 in the USA, and ε(NH4+) was 0.2 ± 0.1 in China, but in India, these ratios are significantly different (Table S2)36,37. The extremely high ALWC in ammonia-rich Delhi triggers the water uptake, causing an enhanced gradient of water activity and water affinity between aerosols and their surroundings, which acts as the major driving force for enhanced secondary inorganic aerosol (SIA) formation during winter.

We estimated the sensitivity of the phase partitioning of nitrate [ε(NO3−)], chloride [ε(Cl−)], and ammonium [ε(NH4+)] to pH, ALWC, and T37,38,39 following the methodology given in Guo et al.37 as:

and,

where γ is the activity coefficient of protonated species in the aqueous medium, and ALWCi is the water associated with the inorganic constituents (µg m−3). H* is the equilibrium constant of HNO3, HCl, and NH3 adopted from 40,41 using molality-based units of mol2 kg−2 atm−142,43. R is the universal gas constant (8.314 J K−1 mol−1), and the value 0.987 is for the transformation of 1 atm to 1 bar. The equations describe the HNO3–NO3−, NH3–NH4+, and HCl–Cl− partitioning, and the estimated values of ALWCi and T were used to evaluate the phase-partitioning of ε(NO3−), ε(Cl−), and ε(NH4+) at different pH regimes.

Figure 1a–d show the variability of gas-to-particle partitioning of nitrate and chloride with pH, ALWC, and T following Eqs. (1–3). Three prominent zones are shown, in which ε(Cl−) and ε(NO3−) vary between the complete gas-phase (ε ~ 0%) to the complete particle-phase (ε ~ 100%). In region I, ε(Cl−) and ε(NO3−) asymptotically approach 0, and the total species primarily remain in the gas-phase. In region III, ε asymptotically approaches 1, and the whole species is in the particle phase, whereas in region II, ε varies between 0 and 1, and the species remain as a mixture of gas and particle-phase. A thermodynamical sweet spot, pH50 is defined, where ε(Cl−) and ε(NO3−) are 0.5, and total chloride and nitrate remain 50% in the gas-phase and 50% in the particle phase.

The gas-to-particle partitioning ratio of (a,c) ε(NO3−), (b,d) ε(Cl−) with pH are shown by the sigmoidal (S) curves. Three zones are defined, where ε(NO3−) and ε(Cl−) vary between the complete gas-phase (ε ~ 0; blue-color zone, region I) and complete particle-phase (ε ~ 100%; blue-color zone, region III). A characteristic pH (pH50) has been defined where total nitrate and chloride remained 50% in the gas-phase and 50% in the particle-phase (gray-color zone, region II). The 1100 hourly data points are color coded with ALWC (1a, 2b) and temperature (1c, 2d) to investigate the impact of ALWC and temperature on the phase-partitioning ratio (ε). The ALWC color bar shows how positively ALWC impacts the phase-partitioning, whereas the temperature color bar shows its inverse relationship with the particle-phase loading over a location.

The 1100 hourly data points are color coded with ALWC (1a, 1b) and temperature (1c, 1d) to investigate the impact of ALWC and temperature on the phase-partitioning ratio (ε). The ALWC color bar shows how positively ALWC impacts the phase-partitioning, whereas the temperature color bar shows its inverse relationship with the particle-phase loading over a location. The ALWC color-coded data points show the variability in ε due to variation in ALWC at a constant pH, indicating the significance of ALWC in modulating the particle-phase loading. From Fig. 1a, it can be seen that ε(NO3−) can vary between 70 and 100% at pH ≈ 4, due to variability in ALWC.

The red line is the sigmoidal curve (S) fitting, which depicts the sensitivity of ε to pH and ALWC over a location. The estimated average pH of PM1 and PM2.5 of 4.49 and 4.58, respectively (“Formulation of "sensitivity regimes" of aerosols in the IGP”), falls on the flat side of the S curves in the blue-color zone. Here, chloride and nitrate remain almost exclusively in the particle-phase, and ammonia remains primarily in the gas-phase. The possibility that emerges as the pathway to modulate high aerosol concentrations is reducing precursors like HCl and HNO3, which is feasible as it would not adversely affect agricultural productivity and potentially impact the ecosystem.

We explored a thermodynamically consistent mathematical framework to reduce ammonium chloride and ammonium nitrate concentrations in the IGP. Following Eqs. (5–13), and using pH and ALWC as coordinates, we defined six "sensitivity regimes" in Fig. 2, where aerosols are sensitive to HCl, HNO3, and NH3 perturbations. As chloride dominates the inorganic mass fraction of fine aerosol in the IGP, we specifically investigate the aerosol mass sensitivity to HCl emissions to define an "HCl sensitive regime". This is significantly different from studies conducted in China, USA, and Europe where researchers have investigated aerosol sensitivity to only HNO3 and NH3 and did not investigate the "HCl sensitive regime"4,5,6.

The conceptual thermodynamical framework with the "coordinates" being pH and aerosol liquid water content (ALWC). Aerosols are sensitive to HCl, HNO3, and NH3 concentrations. The black, blue, and red lines define the "characteristic pH" where chemical regimes are separated by a pre-defined threshold value of ε = 10%. Above pH′, pH″, and pH‴ aerosols are sensitive to the variation in precursors like HCl, HNO3, and NH3, and below which aerosols are deemed insensitive to variations in precursors. The six chemical regimes of (a) PM1 and (b) PM2.5 aerosols are shown in different color shades, with the observational data points colored by mass concentrations (µg m−3) of PM1 and PM2.5 monitored during the Winter Fog Experiment (WiFEx) field campaign of 2017–18.

The "sensitivity regimes" are shown in different colors, where 1100 hourly observational data points are plotted to check the instantaneous response of aerosol loading to HCl, HNO3, and NH3 variability. Figure 2a,b show PM1 and PM2.5 aerosols to remain in the gray shaded region, where aerosols respond proportionally to changes in the HCl and HNO3 emissions but tend to be insensitive to NH3 emissions. The off-white region is defined as an "HCl and HNO3 sensitive regime", indicating that HCl and HNO3 reduction would be the most effective pathway in controlling aerosol pollution over IGP. The blue-shaded region is defined as an "HCl, HNO3, and NH3 sensitive regime", where aerosols are sensitive to HCl, HNO3, and NH3. It can be seen that very few observational data points fall in this regime, in contrast to the USA, where most of the aerosols fall in this regime4. During winter, ammonia concentration is much higher than HCl and HNO3 in IGP, but fine aerosols are not sensitive to NH3 variations. Instead, HCl and HNO3 are by far the limiting factors in aerosol loading, which should be controlled by controlling the major HCl and NOx emissions over IGP.

In Fig. 2a,b, all the data points are color-coded with PM1 and PM2.5 concentrations, respectively, demonstrating that higher aerosol loading is often associated with higher ALWC. Interestingly, the ALWC usually ranges between tens to hundreds of micrograms per cubic meter in climatic regions like China and USA4,5,6, but over IGP the ALWC is an order of magnitude higher, sometimes reaching ~ 2000 to 2400 µg m−3 for PM1 and PM2.5. This high ALWC significantly influences the loading of SIA and causes a reduction in visibility as evidenced by Henry's law, which shows that particles with high ALWC would take up more gaseous pollutants, and the equilibrium would lead to an increase in their water content and the formation of more secondary aerosols like ammonium chloride and ammonium nitrate43,44,45,46.

In this work, we explore the physical, chemical, and thermodynamical processes influencing aerosol loading by considering the dominant impact of multiphase and heterogeneous chemical processes on the aerosol growth processes in the wintertime polluted atmosphere in the IGP. From a thermodynamic perspective, it can be seen that acidic gases are first absorbed on the surface layer of aerosols, and heterogeneous reactions rapidly occur at the surface resulting in the rapid increase in secondary aerosol mass concentrations. The newly formed particle mass then gets dispersed through the liquid phase in high ALWC, where multiphase reactions govern. The wintertime ambient atmosphere thermodynamically favors the condensation of species like nitrate and chloride. The vapor pressure of HNO3 and HCl decreases with temperature exponentially, which reacts with lofted NH3 to stabilize the nitrate and chloride in the particle phase47. The higher surface-area-to-volume ratio of PM1 than PM2.5 suggests the heterogeneous chemistry in PM1 to be more crucial than PM2.5.

This study significantly differs from studies conducted in India and other parts of the world. For example, recent studies conducted in Europe, the USA, and China suggested that ammonia reduction is more cost-effective than NOx reduction and would be the most effective pathway to reduce aerosol loading4,38,48. But in this study, we argue that the sensitivity and effectiveness of the adopted reduction mechanisms are more crucial than the cost-effectiveness. If aerosols are not sensitive to the reductions of a specific precursor, then the cost-effectiveness would not assist in developing an effective mitigation policy.

While studies in India, like Gunthe et al.24 have used a fixed concentration of total chloride (i.e., HCl + Cl− = 25 µg m−3) for the ISORROPIA-II simulations. This average value does not describe the actual variability of ambient chloride concentrations during the wintertime pollution episodes in Delhi. Instead of this average value, we used measured HCl + Cl− concentrations that varied between 2.11–112.21 µg m−3 and 2.31–186.31 µg m−3 for PM1 and PM2.5, respectively. Thermodynamic model simulations with these high input concentrations helped to capture the different atmospheric concentrations, like clear atmosphere and dense foggy conditions20.

In addition, the total nitrate (HNO3 + NO3−) of PM1 (2.33–59.96 µg m−3) and PM2.5 (2.49–85.58 µg m−3) also showed considerable variability during the study period. The total ammonia (NH3 + NH4+) varied between 2.15–120.98 and 3.22–159.39 µg m−3 for PM1 and PM2.5, respectively. Usage of these measured concentrations as the model input can make the thermodynamic model calculations more reliable, capturing the ambient variability more accurately.

To emphasize how sensitive model simulations are to the gas and aerosol phase inputs, the thermodynamic model can reallocate some small amount of NH4+ into NH3 + H+, to make aerosol pH, aerosol liquid water content (ALWC) and phase partitioning ratio self-consistent26. For example, the transfer of just 1 nmol m−3 of NH4+ into NH3 and H+ can be enough to reduce the pH by several units. To conclude, the results that we present illustrate the importance of using thermodynamically consistent sensitivity analysis to effectively address the particulate matter pollution mitigation problem in the Indian region.

This study presents a thermodynamically consistent "roadmap" for effective aerosol mitigation in the Indo-Gangetic Plain (IGP). The proposed framework categorically considers aerosol acidity (pH) and aerosol liquid water content (ALWC) as the governing parameters that modulate the mass loading of aerosols. The mass loading of PM1 and PM2.5 increases with RH, attributing to the water uptake by hygroscopic constituents and enhanced multiphase reactions in highly humid conditions resulting in visibility reduction.

The identified high ammonium chloride concentration is a special characteristic of air pollution in India, and our results show that aerosols with higher chloride mass fraction absorb more water from the air, triggering the positive feedback between water uptake and enhanced growth of aerosols. This augments the particle mass loading and demonstrates the role of ALWC in the growth of PM1 and PM2.5 aerosols.

It is crucial to note that the major errors inherent with thermodynamic equilibrium modeling are derived from two assumptions, i.e., internal mixing and gas-aerosol equilibria, and the analytic error in monitoring individual species cannot be evaluated in the absence of specifically designed quality-control experiments.

This study is in significant contrast to any previous studies conducted in India. For instance, due to the lack of observations, studies like Gunthe et al.24, and Chen et al.28 have used fixed concentrations (25 µg m−3) of HCl + Cl− concentrations for the ISORROPIA-II simulations, which could make the thermodynamic model calculations less constrained. As thermodynamic equilibrium significantly depends on both gas and particle phase concentrations, it was much needed to use simultaneously measured gas and aerosol concentrations as model input. So, we used simultaneously measured total chloride (HCl + Cl−), total nitrate (HNO3 + NO3−), and total ammonia (NH3 + NH4+) concentrations as model input, which can assist in capturing the ambient variability more accurately. The sigmoidal curves and the sensitivity analysis in Fig. 1a–d show the significance of model-derived parameters like pH and ALWC in modulating the phase-partitioning and particle-phase loading of chloride and nitrate over a location like Delhi.

The sensitivity analyses and statistical parameters like Mean Bias (MB), Normalized Mean Bias (NMB), and Root Mean Square Error (RMSE) presented in this study illustrate the necessity of accounting the actual ambient measurements of gas and aerosol species in the ISORROPIA-II thermodynamic model simulations. For example, the high NMB and RMSE values in pH and ALWC show that the fixed HCl + Cl− = 25 µg m−3 cannot represent the real-world variabilities of 2.11–112.21 µg m−3 and 2.31–186.31 µg m−3 for PM1 and PM2.5, respectively.

In this study, we investigate the sensitivity of ammonium chloride and ammonium nitrate loading to reductions in gaseous precursors like HCl, HNO3, and NH3 as a way to reduce PM1 and PM2.5 mass loading. Six "sensitivity regimes" are defined as (a) "HNO3 sensitive", (b) "HCl and HNO3 sensitive", (c) "HCl, NH3 and HNO3 sensitive", (d) "HNO3 and NH3 sensitive", (e) "NH3 sensitive", and (f) "insensitive" to explore the sensitivity of gaseous precursors in the formation of fine aerosols. The application of the framework to the observational dataset shows PM1 and PM2.5 aerosols to fall in the "HCl and HNO3 sensitive regime", emphasizing that HCl and HNO3 reductions would be the most effective pathway to reduce aerosol loading over IGP.

The thermodynamic control suggests that although NH3 concentrations are much higher than HCl and HNO3, the NH3 reduction should not be explicitly prioritized to mitigate aerosols. This is in disparity with the existing aerosol mitigation strategies over IGP, which are mostly ineffective and thermodynamically unfavorable49,50. Our study differs from similar studies in China, despite the fact that there is a very similar range of phase-partitioning characteristics in winter51.

To the best of our knowledge, we are the first to examine the thermodynamic situation via the ε values that show mitigating HCl and HNO3 is more effective, not ammonia for IGP in wintertime. We also stress that it is important to determine which precursors the aerosol mass concentrations are sensitive to for determining a cost-effective mitigation policy. If aerosols are not sensitive to the reductions of a specific precursor, then the cost-effectiveness would not assist in developing an effective mitigation policy. Following the framework advanced in this study, India may develop its own and more precise solution to its severe air pollution problem during winter.

Factors that strongly influence aerosol properties are aerosol acidity (pH) and aerosol liquid water content (ALWC)47,51. The direct monitoring of aerosol pH is highly challenging since an established analytical method of directly determining aerosol pH does not exist26,52. Previous studies used several proxy methods like H+ ion concentration, ion balance, ammonium-to-sulfate ratio, and cation-to-anion ratio to indirectly infer the fine particle pH53. Pye et al.53 showed that these methods could not estimate aerosol acidity accurately, and a small bias in estimated pH may cause substantial errors in the chemical and thermodynamical properties of aerosols21,54. To accurately estimate pH and ALWC, numerous studies have used thermodynamic models like E-AIM, MARS, and EQUISOL. Past studies have shown that ISORROPIA-II can predict pH and ALWC with great accuracy54,55,56.

The estimated pH and ALWC of PM1 and PM2.5 are used as coordinates to represent aerosols in different "sensitivity regimes" where aerosol mass is sensitive to HNO3, HCl, and/or NH3 variability4. Three parameters Ψ, Ω, and φ are defined as:

Substituting these three parameters in Eqs. (1–3), the partitioning ratio of ε(NO3−), ε(Cl−), and ε(NH4+) can be expressed in terms of Ψ, Ω, φ, and as a function of ALWC, like

and

To check the sensitivity of these partitioning fractions to aerosol pH and ALWC, we assume a characteristic partitioning ratio (ε) value of 0.1 (10%) and define the threshold value of ε(NO3−), ε(Cl−), and ε(NH4+) as α, β, and γ4. Above the threshold value (ε = 10%) the PM1 and PM2.5 aerosols are sensitive to HNO3, HCl, or NH3 emissions.

Based on these "characteristic" threshold values, we define "characteristic acidity" for nitrate, (pH'), chloride (pH") and ammonium (pH"'), respectively, as

and

which vary logarithmically with ALWCi.

Using these characteristic values as coordinates, six "sensitivity regimes" are defined as:

Regime I: pH < pH′, pH″ and pH > pH‴, where aerosol mass concentration is not sensitive to HNO3, HCl or NH3 change and defined as "insensitive",

Regime II: pH > pH′, pH < pH″ and pH > pH‴, where aerosol mass concentration is sensitive to HNO3 and not sensitive to HCl and NH3 perturbations and defined as "HNO3 sensitive",

Regime III: pH > pH′, pH″ and pH > pH‴, where aerosol mass concentration is sensitive to HCl, HNO3 and not sensitive to NH3 perturbations and defined as "HCl and HNO3 sensitive",

Regime IV: pH > pH′, pH″ and pH < pH‴, where aerosol mass concentration is sensitive to HCl, HNO3, and NH3 perturbations and defined as "HCl, NH3 and HNO3 sensitive",

Regime V: pH > pH′ and pH < pH″, pH‴, where aerosol mass concentration is sensitive to HNO3 and NH3 and not sensitive to HCl and defined as "HNO3 and NH3 sensitive", and

Regime VI: pH < pH′, pH″ and pH < pH‴, where aerosol mass concentration is sensitive to NH3 and not sensitive to HCl and HNO3 and defined as "NH3 sensitive".

The ISORROPIA-II thermodynamic equilibrium model code is available at http://isorropia.epfl.ch. The data used to prepare this manuscript can be found at https://doi.org/10.17605/OSF.IO/6HGS7.

Ghude, S. D. et al. Premature mortality in India due to PM25 and ozone exposure. Geophys. Res. Lett. 43(9), 4650–4658 (2016).

Article ADS CAS Google Scholar

World Health Organization. Ambient Air Pollution: A Global Assessment of Exposure and Burden of Disease (World Health Organization, 2016).

Google Scholar

Guo, H. et al. Effectiveness of ammonia reduction on control of fine particle nitrate. Atmos. Chem. Phys. 18(16), 12241–12256 (2018).

Article ADS CAS Google Scholar

Nenes, A., Pandis, S. N., Weber, R. J. & Russell, A. Aerosol pH and liquid water content determine when particulate matter is sensitive to ammonia and nitrate availability. Atmos. Chem. Phys. 20(5), 3249–3258 (2020).

Article ADS CAS Google Scholar

Nenes, A. et al. Aerosol acidity and liquid water content regulate the dry deposition of inorganic reactive nitrogen. Atmos. Chem. Phys. Discuss. 10, 25 (2020).

Google Scholar

Ibikunle, I. et al. Fine particle pH and sensitivity to NH3 and HNO3 over summertime South Korea during KORUS-AQ. Atmos. Chem. Phys. Discuss. 20, 1–42 (2020).

Google Scholar

Ravishankara, A. R., David, L. M., Pierce, J. R. & Venkataraman, C. Outdoor air pollution in India is not only an urban problem. Proc. Natl. Acad. Sci. 117(46), 28640–28644 (2020).

Article ADS CAS PubMed PubMed Central Google Scholar

Ghude, S. D. et al. Reductions in India’s crop yield due to ozone. Geophys. Res. Lett. 41(15), 5685–5691 (2014).

Article ADS Google Scholar

Sinha, B. & Sinha, V. Source apportionment of volatile organic compounds in the northwest Indo-Gangetic Plain using a positive matrix factorization model. Atmos. Chem. Phys. 19(24), 15467–15482 (2019).

Article ADS Google Scholar

Kumar, R. et al. How will air quality change in South Asia by 2050?. J. Geophys. Res. Atmos. 123(3), 1840–1864 (2018).

Article ADS Google Scholar

Ghude, S. D. et al. Winter fog experiment over the Indo-Gangetic plains of India. Curr. Sci. 20, 767–784 (2017).

Article Google Scholar

Jena, C. et al. Performance of high resolution (400 m) PM2.5 forecast over Delhi. Sci. Rep. 11(1), 1–9 (2021).

Article Google Scholar

Cappa, C. D., Lovejoy, E. R. & Ravishankara, A. R. Evidence for liquid-like and nonideal behavior of a mixture of organic aerosol components. Proc. Natl. Acad. Sci. 105(48), 18687–18691 (2008).

Article ADS CAS PubMed PubMed Central Google Scholar

Hodzic, A. et al. Rethinking the global secondary organic aerosol (SOA) budget: Stronger production, faster removal, shorter lifetime. Atmos. Chem. Phys. 16(12), 7917–7941 (2016).

Article ADS CAS Google Scholar

Shrivastava, M. et al. Recent advances in understanding secondary organic aerosol: Implications for global climate forcing. Rev. Geophys. 55(2), 509–559 (2017).

Article ADS Google Scholar

Pye, H. O. et al. Coupling of organic and inorganic aerosol systems and the effect on gas–particle partitioning in the southeastern US. Atmos. Chem. Phys. 18(1), 357–370 (2018).

Article ADS CAS PubMed PubMed Central Google Scholar

Gani, S. et al. Submicron aerosol composition in the world’s most polluted megacity: The Delhi Aerosol Supersite study. Atmos. Chem. Phys. 19(10), 6843–6859 (2019).

Article ADS CAS Google Scholar

Ghude, S.D. et al. WiFEX: Walk into the Warm Fog over Indo-Gangetic Plain Region. Bull. Am. Meteorol. Soc. 104(5), E980–E1005 (2023).

Article Google Scholar

Acharja, P. et al. Characterization of atmospheric trace gases and water soluble inorganic chemical ions of PM1 and PM2.5 at Indira Gandhi International Airport, New Delhi during 2017–18 winter. Sci. Total Environ. 729, 138800 (2020).

Article ADS CAS PubMed Google Scholar

Acharja, P. et al. Enhanced secondary aerosol formation driven by excess ammonia during fog episodes in Delhi, India. Chemosphere 20, 133155 (2021).

Google Scholar

Hennigan, C. J., Izumi, J., Sullivan, A. P., Weber, R. J. & Nenes, A. A critical evaluation of proxy methods used to estimate the acidity of atmospheric particles. Atmos. Chem. Phys. 15(5), 2775–2790 (2015).

Article ADS CAS Google Scholar

Bian, Y. X., Zhao, C. S., Ma, N., Chen, J. & Xu, W. Y. A study of aerosol liquid water content based on hygroscopicity measurements at high relative humidity in the North China Plain. Atmos. Chem. Phys. 14(12), 6417–6426 (2014).

Article ADS CAS Google Scholar

Ali, K. et al. Characterization and source identification of PM2.5 and its chemical and carbonaceous constituents during Winter Fog Experiment 2015–16 at Indira Gandhi International Airport, Delhi. Sci. Total Environ. 662, 687–696 (2019).

Article ADS CAS PubMed Google Scholar

Gunthe, S. S. et al. Enhanced aerosol particle growth sustained by high continental chlorine emission in India. Nat. Geosci. 14(2), 77–84 (2021).

Article ADS CAS Google Scholar

Gopinath, A. K. et al. Complex interplay between organic and secondary inorganic aerosols with ambient relative humidity implicates the aerosol liquid water content over India during wintertime. J. Geophys. Res. Atmos. 127(13), e2021JD036430 (2022).

Article ADS CAS Google Scholar

Murphy, J. G. et al. Observational constraints on particle acidity using measurements and modelling of particles and gases. Faraday Discuss. 200, 379–395 (2017).

Article ADS CAS PubMed Google Scholar

Su, H., Cheng, Y. & Pöschl, U. New multiphase chemical processes influencing atmospheric aerosols, air quality, and climate in the anthropocene. Acc. Chem. Res. 53(10), 2034–2043 (2020).

Article CAS PubMed PubMed Central Google Scholar

Chen, Y. et al. Ammonium chloride associated aerosol liquid water enhances haze in Delhi, India. Environ. Sci. Technol. 56(11), 7163–7173 (2022).

Article ADS CAS PubMed PubMed Central Google Scholar

Liu, Y. et al. Ammonium nitrate promotes sulfate formation through uptake kinetic regime. Atmos. Chem. Phys. Discuss. 20, 1–39 (2021).

CAS Google Scholar

Bertram, T. H. & Thornton, J. A. Toward a general parameterization of N2O5 reactivity on aqueous particles: The competing effects of particle liquid water, nitrate and chloride. Atmos. Chem. Phys. 9(21), 8351–8363 (2009).

Article ADS CAS Google Scholar

Song, S. et al. Fine-particle pH for Beijing winter haze as inferred from different thermodynamic equilibrium models. Atmos. Chem. Phys. 18(10), 7423–7438 (2018).

Article ADS CAS Google Scholar

Davis, R. D., Lance, S., Gordon, J. A., Ushijima, S. B. & Tolbert, M. A. Contact efflorescence as a pathway for crystallization of atmospherically relevant particles. Proc. Natl. Acad. Sci. 112(52), 15815–15820 (2015).

Article ADS CAS PubMed PubMed Central Google Scholar

Shi, G. et al. pH of aerosols in a polluted atmosphere: Source contributions to highly acidic aerosol. Environ. Sci. Technol. 51(8), 4289–4296 (2017).

Article ADS CAS PubMed Google Scholar

Stanier, C. O., Khlystov, A. Y., Chan, W. R., Mandiro, M. & Pandis, S. N. A method for the in situ measurement of fine aerosol water content of ambient aerosols: The Dry-Ambient Aerosol Size Spectrometer (DAASS) Special Issue of Aerosol Science and Technology on Findings from the Fine Particulate Matter Supersites Program. Aerosol Sci. Technol. 38(1), 215–228 (2004).

Article ADS CAS Google Scholar

Ali, U., Faisal, M., Ganguly, D., Kumar, M. & Singh, V. Analysis of aerosol liquid water content and its role in visibility reduction in Delhi. Sci. Total Environ. 20, 161484 (2023).

Article ADS Google Scholar

Nah, T. et al. Characterization of aerosol composition, aerosol acidity, and organic acid partitioning at an agriculturally intensive rural southeastern US site. Atmos. Chem. Phys. 18(15), 11471–11491 (2018).

Article ADS CAS Google Scholar

Guo, H. et al. Fine-particle water and pH in the southeastern United States. Atmos. Chem. Phys. 15(9), 5211–5228 (2015).

Article ADS CAS Google Scholar

Xu, J. et al. Importance of gas-particle partitioning of ammonia in haze formation in the rural agricultural environment. Atmos. Chem. Phys. 20(12), 7259–7269 (2020).

Article ADS CAS Google Scholar

Pitzer, K. S. & Mayorga, G. Thermodynamics of electrolytes. II. Activity and osmotic coefficients for strong electrolytes with one or both ions univalent. J. Phys. Chem. 77(19), 2300–2308 (1973).

Article CAS Google Scholar

Clegg, S. L., Brimblecombe, P. & Wexler, A. S. Thermodynamic model of the system H+−NH4+−SO42-−NO3−−H2O at tropospheric temperatures. J. Phys. Chem. A 102(12), 2137–2154 (1998).

Article CAS Google Scholar

Clegg, S. L. & Brimblecombe, P. Solubility of ammonia in pure aqueous and multicomponent solutions. J. Phys. Chem. 93(20), 7237–7248 (1989).

Article CAS Google Scholar

Nenes, A., Pandis, S. N. & Pilinis, C. Continued development and testing of a new thermodynamic aerosol module for urban and regional air quality models. Atmos. Environ. 33(10), 1553–1560 (1999).

Article ADS CAS Google Scholar

Patwardhan, V. S. & Kumar, A. A unified approach for prediction of thermodynamic properties of aqueous mixed-electrolyte solutions. Part I: Vapor pressure and heat of vaporization. AIChE J. 32(9), 1419–1428 (1986).

Article CAS Google Scholar

Ravishankara, A. R. Heterogeneous and multiphase chemistry in the troposphere. Science 276(5315), 1058–1065 (1997).

Article CAS Google Scholar

Hanson, D. R., Ravishankara, A. R. & Solomon, S. Heterogeneous reactions in sulfuric acid aerosols: A framework for model calculations. J. Geophys. Res. Atmos. 99(D2), 3615–3629 (1994).

Article ADS CAS Google Scholar

Rood, M. J., Shaw, M. A., Larson, T. V. & Covert, D. S. Ubiquitous nature of ambient metastable aerosol. Nature 337(6207), 537–539 (1989).

Article ADS CAS Google Scholar

Seinfeld, J., & Pandis, S., Atmospheric Chemistry and Physics. 1997. New York (2008).

Gu, B. et al. Abating ammonia is more cost-effective than nitrogen oxides for mitigating PM2.5 air pollution. Science 374(6568), 758–762 (2021).

Article ADS CAS PubMed Google Scholar

Sharma, V. S., et al. Breathing cleaner air: Ten scalable solutions for Indian Cities. The Energy and Resources Institute. A self-organized task force report for the World Sustainable Development Summit, New Delhi (2016).

Ganguly, T., Selvaraj, K. L. & Guttikunda, S. K. National Clean Air Programme (NCAP) for Indian cities: Review and outlook of clean air action plans. Atmos. Environ. X 8, 100096 (2020).

CAS Google Scholar

Weber, R. J., Guo, H., Russell, A. G. & Nenes, A. High aerosol acidity despite declining atmospheric sulfate concentrations over the past 15 years. Nat. Geosci. 9(4), 282–285 (2016).

Article ADS CAS Google Scholar

Rindelaub, J. D. et al. Direct measurement of pH in individual particles via Raman microspectroscopy and variation in acidity with relative humidity. J. Phys. Chem. A 120(6), 911–917 (2016).

Article CAS PubMed Google Scholar

Pye, H. O. et al. The acidity of atmospheric particles and clouds. Atmos. Chem. Phys. 20(8), 4809–4888 (2020).

Article ADS CAS PubMed PubMed Central Google Scholar

Nault, B. A. et al. Chemical transport models often underestimate inorganic aerosol acidity in remote regions of the atmosphere. Commun. Earth Environ. 2(1), 1–13 (2021).

Article Google Scholar

Fountoukis, C. & Nenes, A. ISORROPIA II: A computationally efficient thermodynamic equilibrium model for K+–Ca2+–Mg2+–NH4+–Na+–SO42−–NO3−–Cl−–H2O aerosols. Atmos. Chem. Phys. 7(17), 4639–4659 (2007).

Article ADS CAS Google Scholar

Zhai, S. et al. Control of particulate nitrate air pollution in China. Nat. Geosci. 14(6), 389–395 (2021).

Article ADS CAS Google Scholar

Download references

The Indian Institute of Tropical Meteorology (IITM) is funded by the Ministry of Earth Sciences (MoES), Government of India. We thank the Director, IITM, for the continuous support and encouragement. The authors are grateful to the Ministry of Earth Sciences (MoES), Government of India, for the support throughout the Winter Fog Experiment (WiFEX) field campaign. The lead author's fellowship was supported by the National Supercomputing Mission (NSM) program grant at C-DAC, and we are grateful to the Executive Director and the Director-General of C-DAC. The National Center for Atmospheric Research (NCAR) is sponsored by the National Science Foundation (NSF).

Indian Institute of Tropical Meteorology, Ministry of Earth Sciences, Pune, India

Prodip Acharja, Sachin D. Ghude, Gaurav Govardhan & Kaushar Ali

Laboratoire des Sciences du Climat et de l’Environnement, LSCE, CNRS, Gif-sur-Yvette, France

Prodip Acharja & Jean-Eudes Petit

Department of Earth and Environmental Sciences, Indian Institute of Science Education and Research Mohali, Sahibzada Ajit Singh Nagar, Punjab, India

Baerbel Sinha & Vinayak Sinha

National Center for Atmospheric Research, Boulder, CO, 80307, USA

Mary Barth & Rajesh Kumar

Savitribai Phule Pune University, Pune, 411007, India

Rachana Kulkarni

Engineering and Applied Science, Ontario Technical University, Oshawa, ON, Canada

Ismail Gultepe

Civil and Environment Eng and Earth Sciences, University of Notre Dame, Notre Dame, IN, USA

Ismail Gultepe

Ministry of Earth Science, Lodhi Road, New Delhi, India

Madhavan Nair Rajeevan

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

P.A. and S.D.G. conceived the idea and conducted the research; P.A. performed the ISORROPIA-II model simulations and analysis, and wrote the initial draft of the manuscript; B.S., M.B., V.S., R.K., and I.G. contributed to the interpretation of the results and provided critical feedback to shape the analysis, research, and manuscript. P.A., S.D.G., and J.-E.P., contributed to the thorough revision of the manuscript and addressed the responses to the reviewer's comments in detail; B.S., M.B., R.K., V.S., R.K., K.A., G.G., I.G., J.-E.P., and M.N.R. formulated the research and edited the manuscript. All authors reviewed the manuscript.

Correspondence to Sachin D. Ghude or Baerbel Sinha.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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 licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence 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 licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Acharja, P., Ghude, S.D., Sinha, B. et al. Thermodynamical framework for effective mitigation of high aerosol loading in the Indo-Gangetic Plain during winter. Sci Rep 13, 13667 (2023). https://doi.org/10.1038/s41598-023-40657-w

Download citation

Received: 21 November 2022

Accepted: 16 August 2023

Published: 22 August 2023

DOI: https://doi.org/10.1038/s41598-023-40657-w

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.