Impact of reduced anthropogenic emissions during COVID-19 on air quality in India

To mitigate the impacts of the pandemic of coronavirus disease 2019 (COVID-19), the Indian government implemented lockdown measures on 24 March 2020, which prohibited unnecessary anthropogenic activities, thus leading to a significant reduction in emissions. To investigate the impacts of this lockdown measure on air quality in India, we used the Community Multi-Scale Air Quality (CMAQ) model to estimate the changes of key air pollutants. From pre-lockdown to lockdown periods, improved air quality is observed in India, indicated by the lower key pollutant levels such as PM2.5 (−26 %), maximum daily 8 h average ozone (MDA8 O3) (−11 %), NO2 (−50 %), and SO2 (−14 %). In addition, changes in these pollutants show distinct spatial variations with the more important decrease in northern and western India. During the lockdown, our results illustrate that such emission reductions play a positive role in the improvement of air quality. Significant reductions of PM2.5 concentration and its major components are predicted, especially for secondary inorganic aerosols that are up to 92 %, 57 %, and 79 % for nitrate (NO−3 ), sulfate (SO 2− 4 ), and ammonium (NH+4 ), respectively. On average, the MDA8 O3 also decreases 15 % during the lockdown period although it increases slightly in some VOC-limited urban locations, which is mainly due to the more significant reduction of NOx than VOCs. More aggressive and localized emission control strategies should be implemented in India to mitigate air pollution in the future.


Introduction
India, the second-most populous country in the world, has been suffering from severe air pollution along with rapid urbanization and industrialization in recent decades (Karambelas et al., 2018), and 13 Indian cities were among the world's top 20 most polluted cities according to the World Health Organization (WHO) (WHO, 2018). High-level pollution leads to health risks and ecosystem damage, which caused 1.24 million deaths in India in 2017 (Balakrishnan et al., 2019) and a great loss of crops (Oksanen et al., 2013;Lal et al., 2017). To mitigate air pollution, the Indian government has been promoting effective emission control strategies such as the conversion of fossil fuels to clean fuels in the nationwide Clean India Mission (CIM). However, such long-term or short-term reduction strategies seem to show insufficiency in the restoration of ambient air quality (Beig et al., 2013;Purohit et al., 2019;Banerjee et al., 2017).
Due to the pandemic of coronavirus disease 2019 (COVID-19), nationwide or partial lockdown measures have been implemented in many countries (Chintalapudi et al., 2020;Dantas et al., 2020;Ehrlich et al., 2020). The Indian government declared corresponding bans since the detection of the first confirmed case on 30 January 2020. Then, to counter the fast contagion of COVID-19, a 3-week nationwide lockdown was imposed in India on 24 March, which was extended until 30 June. The lockdown measures miti-gate the impact of COVID-19 on Indian health infrastructure and it also helped in curbing the rate of the spread of this infectious disease among people (Pai et al., 2020;Anderson et al., 2020). Because of the prohibition of industrial activities and mass transportation, anthropogenic emissions showed a tremendous reduction. As well as this, several studies showed that dramatic emission reductions had an enormous impact on the formation of air pollution and positively influenced air quality (Isaifan, 2020;Bao and Zhang, 2020;Gautam, 2020). Thus, the lockdown also provides a valuable opportunity to assess the changes in air pollutants with significantly reduced anthropogenic emissions in a short time.
Conspicuous reductions in concentrations of pollutants were also claimed in different regions (Otmani et al., 2020;Dantas et al., 2020;Nakada and Urban, 2020). Most Indian studies claimed the greatest reduction of particulate matter with an aerodynamic diameter of less than 2.5 µm (PM 2.5 ), up to 50 % Mahato et al., 2020;Sharma et al., 2020). However, an increase in ozone (O 3 ) concentrations was observed (Collivignarelli et al., 2020;Sicard et al., 2020) and severe air pollution events still occurred after large emission reductions due to unfavorable meteorological conditions . Moreover, another analysis showed that the effects of lockdown during the COVID-19 pandemic on PM 2. 5 and O 3 pollution levels were smaller than the expected response to the enacted stay-at-home order (Bujin et al., 2020). Hence, the significance and impacts of lockdown measures are still not well understood.
Therefore, it is significant to understand the mechanisms involved in air pollution formation before and after dramatic emission changes comprehensively, in addition to the comparison of air pollution levels. Mahato et al. (2020) concluded that air quality in India from 24 March to 14 April improved sharply according to the change of the National Air Quality Index, especially for Delhi. Srivastava et al. (2020) reported the concentrations of primary air pollutants were drastically lowered as a result of emission reductions. Kumari and Toshniwal (2020) also stated that the concentration of key pollutants such as PM 2.5 in both Delhi and Mumbai shows a decreasing trend. These studies pointed out that the air quality was improved during the lockdown period compared with the period before lockdown and depends on the duration of the lockdown Mor et al., 2021). As well as this, compared with the same period in previous years, Gautam (2020) claimed that aerosol concentration levels are at their lowest in the last 20 years during lockdown based on satellite data. Selvam et al. (2020) stated that the Air Quality Index (AQI) was improved by 58 % in Gujarat state in western India during lockdown (24 March-20 April 2020) compared to 2019. Kabiraj and Gavli (2020) concluded that the mean concentration of PM 2.5 decreased by 42.25 % from January to May in 2020 compared with 2019. Similarly, Das et al. (2020) also showed that great reductions of PM 2.5 were found across cities in the Indo-Gangetic Plain (IGP) compared with 2018 and 2019. How-ever, the role of meteorological conditions and chemical reactions involving changes in air quality is not clear from these observation-based studies, which only showed the phenomenon of concentration reduction and switch of major primary pollutants mainly in urban cities. Further, the number of monitoring stations in the country is way below the guidelines given by the governing bodies and not uniformly distributed, which results in observation data limitations in India .
In this study, the Community Multi-Scale Air Quality (CMAQ) model was used to investigate changes in air pollutants before the lockdown (from 21 February to 23 March 2020) and during the lockdown (from 24 March to 24 April 2020) periods throughout the Indian region. We explored the synergetic impacts from the meteorological conditions and anthropogenic emissions during the pre-lockdown and lockdown periods. As well as this, we directly quantified the change in air quality during the lockdown due to the reduced anthropogenic emissions by comparing the differences between Case 1 (without emission reductions) and Case 2 (with emission reductions). The model performance was evaluated by comparing the simulation results with the observation data, which is collected by the Central Pollution Control Board (CPCB). This study has important implications for developing control strategies to improve air quality in India.

Data collection
We used observed hourly PM 2.5 , O 3 , carbon monoxide (CO), and nitrogen dioxide (NO 2 ) data from 21 February to 24 April 2020 from the CPCB online database (https://app. cpcbccr.com/ccr/#/caaqm-dashboard-all/caaqm-landing, last access: 20 January 2021), which is widely applied in previous studies (Kumar, 2020;Sharma et al., 2020;Srivastava et al., 2020;Shehzad et al., 2020). The CPCB database provides data quality assurance (QA) or quality control (QC) programs by establishing strict procedures for sampling, analysis, and calibration (Gurjar et al., 2016). As well as this, the observed daily averages of PM 2.5 and maximum daily 8 h average ozone (MDA8 O 3 ) have been further calculated to analyze the change in air quality during the pre-lockdown (from 21 February to 23 March 2020) and lockdown (from 24 March to 24 April 2020). The satellite-observed NO 2 and formaldehyde (HCHO) column number density datasets are from the Sentinel-5 Precursor TROPOspheric Monitoring Instrument (S-5P TROPOMI) (https://scihub.copernicus.eu, last access: 7 February 2021). As well as this, we effectively removed the pixels with a QA value less than 0.75 for NO 2 tropospheric column density and 0.5 for HCHO from the datasets to exclude the interferences such as clouds and snow or ice (Apituley et al., 2018).

Model description
This study applied CMAQ (Byun and Schere, 2006) version 5.0.2 with an updated SAPRC-11 photochemical mechanism (Carter, 2011;Hu et al., 2016) and aerosol module (AERO6) (Binkowski and Roselle, 2003) to simulate air pollution across India with a horizontal resolution of 36 km × 36 km (117 × 117 grid cells). Figure 1 shows the simulation domain with positions of the main Indian cities. The simulation was conducted from 21 February to 23 March as a pre-lockdown and 24 March to 24 April as a lockdown period.
The Weather Research & Forecasting model (WRF) version 3.6.1 was utilized to generate meteorology fields driven by the latest FNL (Final) Operational Global Analysis data. Anthropogenic emissions were from the monthly data from the Emissions Database for Global Atmospheric Research (EDGAR) version 4.3 (http://edgar.jrc.ec.europa.eu/ overview.php?v=431, last access: 20 September 2020). The monthly emissions from different source sectors were divided into six major groups, residential, industrial, agriculture, on-road, off-road, and energy, before being adjusted from the base year of 2010 to 2019 based on population and economic growths similar to Guo et al. (2017), and the adjustment factors are shown in Tables S1-S3 in the Supplement. Weekly and diurnal profiles were used to convert monthly emissions to hourly inputs, and the US EPA's SPE-CIATE 4.3 source profiles were used to speciate total particulate matter (PM) and volatile organic compounds (VOCs) to model species (Wang et al., 2014).
The biogenic emissions were derived from The Model of Emissions of Gases and Aerosols from Nature (MEGAN) version 2.1 (Guenther et al., 2012), and the emissions from

Emission reduction during COVID-19
Due to the COVID-19 lockdown, human activities were limited and related anthropogenic emissions were reduced. Different sources were used to obtain changes in anthropogenic emissions from different sectors in comparison to 2019. For the sectors of on-road and off-road, the vehicle emission changes were based on the number of registered vehicles verified from the article (ET Bureau, 2020). The changes in energy demand were obtained from official data released by the Power System Operation Corporation (POSOCO) (Abdi, 2020). Residential and agricultural emissions remain unchanged due to a lack of sufficient information.
For the industrial sector, we classify the Indian industries into three different classes based on the degree of air pollution caused (https://www.indianmirror.com/ indian-industries/environment.html, last access: 7 February 2021) (Table S4) including very polluting (VP), medium polluting (MP), and low polluting (LP) industries. The pollution index (PI) of any industry is a number ranging from 0 to 100 and the increasing value of PI denotes the increasing degree of pollution load from the industry. As well as this, CPCB, State Pollution Control Boards (SPCBs), and the Ministry of Environment, Forest and Climate Change (MoEFCC) have finalized the criteria on the range of PI for the purpose of categorization of the industrial sector (https: //pib.gov.in/newsite/printrelease.aspx?relid=137373, last access: 7 February 2021) (Table 1).
Based on the above definition of the VP, MP, and LP industry, the emissions before lockdown can be expressed as follows: where S VP , S MP , and S LP are 1, 0.6, and 0.4 as the assigned scores, and N VP-pre , N MP-pre , and N LP-pre are the number of each category industry before the lockdown. Similarly, the emissions during the lockdown are as follows: where N VP-lock , N MP-lock , and N LP-lock are the number of functioning industries during the lockdown. Therefore, the percentage reduction of industrial emissions can be calculated as follows: In this study, two sensitivity simulations were conducted during the lockdown periods. Case 1 assumes business as usual with the same emissions as in 2019, while Case 2 adjusts anthropogenic emissions using factors obtained above for different sectors ( Table 2). The differences between Case 2 and Case 1 can be assumed as the effects of COVID-19 lockdowns.
3 Results and discussion

WRF-CMAQ model validation
Meteorology plays an important role in emissions, transport, deposition, and formation of air pollutants . Hence, the performance of WRF is validated to ensure accurate air pollution simulation against available observation from the National Climate Data Center (NCDC). There are more than 1300 stations within the simulation domain with hourly observations. The variables considered include temperature at 2 m above the surface (T2), wind speed (WS), wind direction (WD), and relative humidity (RH). Table S5 shows the statistics of mean observation and mean prediction of meteorological parameters, along with mean bias (MB), gross error (GE), and root mean squared error (RMSE), which are compared to benchmarks suggested by Emery et al. (2001). All the statistical indexes are listed in Table S6.
In general, the WRF model performance is similar to previous studies in India (Kota et al., 2018). For the prelockdown and lockdown period, predicted T2 was underestimated with MB values of −1.5 and −1.2 K, respectively. The GE values for WS were 1.7 % (pre-lockdown) and 1.8 % (lockdown), satisfying the suggested criteria of 2.0 %, and RMSE was slightly over the criteria. The MB values for WD were 3.2 and 2.6 • during the two periods, which are within the criteria of ±10 • . The GE and RMSE for WD were slightly out of the benchmarks. The under-predicted RH was also observed in this study, which was reported in other Asian studies . Those statistic values that did fall in the benchmark were mainly due to the resolution (36 km) applied in this study compared to the finer resolution (4-12 km) suggested in Emery et al. (2001) and Sahu et al. (2020). Table S7 shows the model performance of MDA8 O 3 , PM 2.5 , CO, and NO 2 in five major cities in India including Delhi, Mumbai, Chennai, Hyderabad, and Bengaluru. For PM 2.5 , after excluding some abnormally high values of greater than 300 µg m −3 , the averaged mean fractional bias (MFB) (−0.48) and mean fractional error (MFE) (0.61) values in all the five urban cites met the criteria limits of ±0.6 and 0.75 claimed by the EPA (2007). And the recommended criteria are commonly used for validating air quality model performance in the Indian region Kota et al., 2018). For O 3 , a cut-off value of 40 ppb is applied, which is based on EPA's recommendations (EPA, 2005). As well as this, the model was able to reproduce the variation trends of observed hourly O 3 in all these major cities, although slightly over-estimations have occurred. And averaged MFB (−0.05) and MFE (0.25) values of O 3 also satisfy the benchmarks of ±0.15 and 0.30 set by the EPA (2005) in most of these cities with Chennai and Hyderabad exceeding the limits slightly. The performance of PM 2.5 , NO 2 , O 3 , and CO in these urban areas was also similar to Kota et al. (2018), which could provide robust results for the following air quality study.
To further validate modeled HCHO and NO 2 , we compared our simulated results with satellite-observed data during pre-lockdown and lockdown periods (Fig. S1). The CMAQ predicted vertical column densities (VCDs) of tropospheric NO 2 and HCHO were calculated using Eq. (4) (Eskes et al., 2020).
where n equals 17 as the number of vertical layers in the model (with the highest layer height of ∼ 10 km), C i means species concentration (ppm), H i represents each layer height (m), and α is the coefficient for converting units from ppm to molec. cm −2 . The predicted regional distribution of tropospheric column NO 2 and HCHO is similar to satellite observations. Overall, HCHO and NO 2 are higher in eastern and northern India than in other regions. And their variation trends from CMAQ and TROPOMI are consistent so that NO 2 decreases while HCHO increases during the lockdown. We also acknowledge that the uncertainty of the emission inventory and chemical mechanism in the modeling may affect the simulated results (Dominutti et al., 2020;Kitayama et al., 2019). As well as this, the observed peak values of PM 2.5 in each city also decrease appreciably (up to 57 %) during the lockdown period. On 24 March, the first day of lockdown, a significant drop in PM 2.5 concentration due to the emission reduction of primary pollutants is observed (Fig. S2). However, most of the PM 2.5 concentrations are still above the WHO annual guideline values of 10 µg m −3 (WHO, 2016) during the lockdown period, with peak values over 60 µg m −3 occasionally. Figure 3 shows the temporal variation of MDA8 O 3 in these five cities. The predicted MDA8 O 3 is consistent in trend with observation values on most days, while simulated concentrations are overall higher, particularly in Hyderabad. The observed average MDA8 O 3 during lockdown is higher than that before the lockdown in Delhi (2 %), Hyderabad (12 %), and Bengaluru (2 %). This is likely due to the fact that O 3 formation in these cities is under VOC control , and a nitrogen oxide (NO x ) reduction leads to an O 3 increase by enhanced hydrogen oxide radical (HO x ) concentrations (Zhao et al., 2017). The increase in monthly average T2 from before the lockdown (281.0 K) to lockdown (285.1 K) could also lead to an increase in O 3 (Chen et al., 2019). In contrast, the observed average MDA8 O 3 during lockdown is reduced compared with the pre-lockdown period in both Mumbai (−35 %) and Chennai (−13 %). This could be caused by a much larger reduction in emissions as Mumbai and Chennai with high urbanization and industrialization are the most affected areas. Specifically, more stringent lockdown measures may be implemented in Mumbai than we assumed, which accounted for more than a fifth of infections in India (Mukherjee, 2020). Figure 4 shows the comparison of predicted air pollutants before and during the lockdown throughout India. Generally, decreases in key pollutants including particulate matter with an aerodynamic diameter of less than 10 µm (PM 10 ) (−16 %), PM 2.5 (−26 %), MDA8 O 3 (−11 %), NO 2 (−50 %), and sulfur dioxide (SO 2 ) (−14 %) are calculated across India. Changes in these pollutants present distinct regional variations. In northern and western India, the lower levels of these pollutants are observed during the lockdown, with reductions of PM 2.5 and PM 10 of up to 79 %. In particular, the most significant decreases during the lockdown are found in the populated, industrialized, and polluted IGP region. The average PM 2.5 even drops from approximately 35-70 µg m −3 (pre-lockdown) to 15-40 µg m −3 (lockdown) in these regions because local emissions are generally the largest contributor (38 %-78 %) to PM 2.5 in India (David et al., 2019). However, increases in these key pollutants are found mainly in the northeastern, eastern, and parts of southern India.

Changes in air quality from pre-lockdown to lockdown periods
As well as this, changes in PM 2.5 also show prominent differences in the rural and urban areas. In India, rural areas have different emission sources from urban areas and are less influenced by lockdown measures (Garaga et al., 2020). In megacities such as Delhi, the predicted concentrations of PM 2.5 declined during the lockdown, which is consistent with previous results (Kumari and Toshniwal, 2020;Chauhan and Singh, 2020). For instance, a reduction of over 60 % in PM 2.5 is estimated in Delhi and Ahmedabad. However, increases in PM 2.5 (∼ 20 %) are observed in the remote northeastern part of India. Variations in near-surface mete- orological factors during lockdown also play an important role in PM 2.5 changes. As is shown in Fig. S3, lower PM 2.5 in urban areas during lockdown (Fig. 4) may be attributed to the decrease in RH and increase in planetary boundary layer (PBL) height, while the decrease in precipitation and WS allows PM 2.5 to accumulate in some rural areas (Schnell et al., 2018;Le et al., 2020).
As gaseous precursors of major components to PM 2.5 (Jain et al., 2020), concentrations of NO 2 and SO 2 also decrease significantly in most regions by up to 90 % and 87 %, respectively. However, their levels increase in parts of eastern and southern India, thus leading to higher levels of PM 2.5 and PM 10 in the same regions. MDA8 O 3 is also rising in eastern India by the highest increasing rate of 29 %, while a 30 % reduction is observed in northern and western India. Although significant reductions are found in O 3 precursor emissions throughout India during the lockdown, the MDA8 O 3 has not shown a comparable decrease, which is affected by meteorological conditions such as an increase in temperature and decrease in RH (Fig. S3). Higher temperatures speed up photochemical processes that produce O 3 , while a higher RH re- duces them (Chen et al., 2019;Zhao et al., 2017;Ali et al., 2012).
In summary, the decrease in PM 2.5 , PM 10 , NO 2 , and SO 2 and the increase in MDA8 O 3 during lockdown are consistent with previous results (Srivastava et al., 2020;Mahato et al., 2020). In the case of Delhi, compared with the previous studies, the PM 2.5 reduction (34 %) is comparable with 35 % reported by Chauhan and Singh (2020), with less than 53 % stated by Mahato et al. (2020) and 49 % calculated by Kumari and Toshniwal (2020) during the first phase of lockdown (from 24 March to 15 April 2020). These differences may be caused by the considered duration of the lockdown period. Our study concerns the later lockdown period (after 15 April 2020) when there is an increase in traffic flow and some relaxation of lockdown measures (Kumar, 2020). Moreover, the different characteristics of these air pollutants in rural and urban areas have not been investigated comprehensively in previous studies. Kumari and Toshniwal (2020) also concluded that concentrations of PM 10 , PM 2.5 , and SO 2 tended to rise in Singrauli (rural area, located in central India) during the lockdown, contrary to the results of Delhi and Mumbai. Therefore, our results have important implications for the study of air quality changes and their regional distribution across India and indicate more strident emission reduction policies should be implemented across India, especially in the later phases of lockdown and in rural areas.

Effects of emission reductions on PM 2.5 during the lockdown
There are significant changes in PM 2.5 between the lockdown and pre-lockdown periods. Moreover, we directly quantify the change in PM 2.5 during the lockdown. Figure 5 shows the differences in major PM 2.5 components during the lockdown period with (Case 2) and without (Case 1) control measures. Major components of PM 2.5 including nitrate (NO − 3 ), sulfate (SO 2− 4 ), ammonium (NH + 4 ), elemental carbon (EC), primary organic aerosols (POAs), and secondary organic aerosols (SOAs) decreased significantly in Case 2 compared to Case 1, indicating the positive effects of emission reduction. Primary components of PM 2.5 (EC and POAs) are lowered by an average of 37 % and 14 %, respectively. EC is usually emitted from combustion sources and a drastic decrease of up to 74 % directly reflected the impact of emission reductions from industry and transportation. Secondary inorganic aerosols (SIAs) including NO − 3 , SO 2− 4 , and NH + 4 and SOAs accounted for most of the PM 2.5 bulk mass (39 %) and showed greater decreases than primary components. Moreover, the spatial distribution of SIAs is similar to PM 2.5 in that the reduction is more significant in the north of India where the decrease in NO − 3 , SO 2− 4 , and NH + 4 are up to 92 %, 57 %, and 79 % respectively. The largest reduction of NO − 3 , averaging 62 %, resulted from transportation reduction, and the SO 2− 4 reduction (averaging 31 %) is likely due to the falling industry emissions (Gawhane et al., 2017;Wang et al., 2020). On average, NH + 4 and SOAs are decreased by 41 % and 14 %, respectively. The significant decrease in NH + 4 cannot be attributed to the absence of reduced agricultural emissions in the simulation but may be due to the relatively reduced (NH 4 ) 2 SO 4 and NH 4 NO 3 in the CMAQ chemistrytransport model (Fountoukis and Nenes, 2007). By contrast, compared with VOCs, an important precursor of SOAs, the smaller reduction of SOAs may be related to the weakening of the atmospheric oxidizing capacity (AOC), which plays an important role in the formation of SOAs (Feng et al., 2019). As well as this, the reduction of NO x may lead to an increase in SOAs offsetting some of the influence by the reduction in VOC emissions (Kroll et al., 2020). Figure 6 shows the predicted response of changes in concentration of primary PM 2.5 (PPM) and secondary components to the reduced emissions of related precursors in Delhi, Mumbai, Kolkata, Bengaluru, Hyderabad, Chennai, Ahmedabad, and Lucknow. Generally, all species decreased with the reduced emissions and the great sensitivity of PM 2.5 component concentrations to emissions showed the important role of meteorology and the effectiveness of stringent measures to reduce emissions.
On average, NO − 3 shares the largest reduction of 77 % mainly driven by the decrease in its gaseous precursor NO x (71 %). At least a 27 % decrease in SO 2− 4 is found in each city caused by the largest reduction of SO 2 (averaged 59 %). An average reduction of over 70 % in NO x and NO − 3 may still be related to the reduction of vehicles. And SOAs dropped by an average of 18 % because of the lack of precursors due to the emission reduction of VOCs (29 %). Due to the reduction of emitting precursors, the concentration reduction of PM 2.5 secondary components is less than that of primary components. The ratios of PPM reduction in emission (averaged 39 %) are larger than the reduction in concentration (averaged 43 %) in five selected cities. In particular, a 7 % reduction in emission of PPM caused a 43 % decline in its concentration in Hyderabad. Emissions of EC and organic carbon (OC) have also been reduced by a certain proportion resulting in a similar or greater reduction in concentrations.
The response of concentration to emissions in all cities presented a nonlinear change that has been confirmed previously by Zhao et al. (2017), which is related to various meteorological conditions . For example, in Lucknow, PPM, EC, OC, SO 2 , NO x , and VOCs decreased by 14 %, 25 %, 8 %, 39 %, 55 %, and 11 % respectively, while the concentration of PPM, EC, OC, SO 2− 4 , NO − 3 , and SOAs dropped by 21 %, 32 %, 12 %, 43 %, 78 %, and 18 %. As well as this, the concentration response to emission reduction is likely to be more prominent in highly polluted and industrialized areas. The highest reductions in PPM and these secondary components of PM 2.5 happened in Ahmedabad (an industrial city located in western India) with high vehicular populations. Bengaluru, a major southern Indian city, is considered as one of the cleaner Indian major cities because of its low PM 2.5 concentrations with no heavy industries (Guttikunda et al., 2019). Consequently, the reduction in PM 2.5 and its major components (especially for secondary compo- nents) in Bengaluru is not as significant as in Ahmedabad although a similar reduction in emissions is observed.

Effects of emission reductions on O 3 during the lockdown
We investigated the changes of MDA8 O 3 and its major precursors NO x and HCHO during the lockdown period. HCHO is one of the major contributors to total VOC reactivity (Zhang et al., 2012;Steiner et al., 2008). It also has a strong correlation with VOC (R 2 up to 0.93) (Fig. S4) and performs well when validated by comparing with satelliteobserved data. As a result, HCHO is used as a good proxy in the model for the total VOCs, consistent with previous studies such as Palmer et al. (2003). Figure 7 shows that MDA8 O 3 , NO x , and HCHO decreased all over India. The average reduction rates of MDA8 O 3 , NO x , and HCHO are approximately 15 %, 50 %, and 15 %, respectively. For both Case 1 and Case 2, the higher levels of MDA8 O 3 are in eastern India (over 60 ppb, Case 1) in which the higher NO x is also observed (over 12 ppb, Case 1) during the lockdown. Compared to PM 2.5 , no significant north-south differences are found in the change of O 3 . The NO x concentration has the greatest reduction that is mostly driven by the large cutting of energy emission by 26 %, which is consistent with the decline of India's electricity consumption (9.2 %) (Reuters, 2020). Figure S5 shows the O 3 production sensitivity (O 3 / NO y ) in India during the lockdown, which is considered as an indicator of O 3 sensitivity to NO x and VOCs (Sillman, 1995;Sillman and He, 2002). As well as this, O 3 / NO y < 6 indicates that O 3 formation is VOC-limited, O 3 / NO y > 8 indicates NO x -limited, and intermediate values are transitional.
In India, NO x -limited regimes are found in vast areas from both Case 1 and Case 2, which was also reported in previous studies (Mahajan et al., 2015). As a result, the large reduction of NO x leads to decreased MDA8 O 3 in most Indian regions. Compared to Case 1, the VOC-limited area expands mainly in the northwest and south of India from Case 2 during the lockdown. Simultaneously, the rise of MDA8 O 3 (averaged 5 % and up to 21 %) is found sporadically in these VOC-limited areas in which more significant decreases in NO x (compared with VOCs) reduce the O 3 consumption (NO + O 3 = NO 2 + O 2 ) and enhance HO x concentrations, resulting in an increase in O 3 levels. It may also indicate that the increase in O 3 is amplified regionally by the expansion of the VOC-limited regimes due to the lockdown. The changes in emissions and concentrations of MDA8 O 3 , HCHO, and NO x showed a non-linear response. In Delhi, a 76 % reduction in NO x emissions resulted in a 77 % reduction in its concentration, while a 29 % reduction in HCHO resulted in only an 11 % reduction. In a megacity like Delhi, about 7 million vehicles and many fossil-fuel-based plants lead to high NO x emissions, and local restricted transportation and industrial activities during lockdown could lead to a significant reduction of primary NO x emissions (Sharma et al., 2016). The concentration of NO x is appreciably highly sensitive to a primary NO x emission reduction. However, the VOCs emission reduction resulting from the lockdown is relatively less than NO x in each city. And most of the reduction of HCHO concentration is less than that of emission reduction, which is different from NO x , which indicated that the change of HCHO concentrations is not dominated by primary HCHO emission reduction.

Conclusion
Compared with before the lockdown, observed PM 2.5 during the lockdown in Delhi, Mumbai, Chennai, Hyderabad, and Bengaluru shows an overall decrease. In contrast, MDA8 O 3 increases in three of these cities. The comparison of predicted air pollutants across India before and during the lockdown shows distinct regional characteristics. The most significant reductions of PM 2.5 and PM 10 (up to 79 %) are observed in most of northern and western India including all these megacities. However, increases in MDA8 O 3 (up to 29 %) and other key pollutants are reported in northeastern, east- ern, and parts of southern India covering most of the rural areas. As well as this, it can be concluded that the synergetic impact from the meteorological conditions and anthropogenic emissions plays an important role in those increases from pre-lockdown to lockdown periods.
The drastic decline in PM 2.5 and its major components during the lockdown period in Case 2 compared with Case 1 shows the positive impacts of emission control measures, especially for SIAs. During the lockdown, the decrease in MDA8 O 3 (averaged 15 %) occurs in most regions in India, which is attributed to the lower emissions of NO x (48 %) and VOCs (6 %) that are precursors of O 3 . Our results demonstrate that the strident emission controls due to the lockdown have mitigated air pollution in India. However, more stringent mitigation measures are needed to achieve effective control of air pollution from secondary air pollutants and their components, particularly in rural areas. We also find the scattered increases in MDA8 O 3 (up to 21 %) in some urban lo-cations in the VOC-limited areas due to the emission reduction. This indicates that a more localized control policy with the consideration of the O 3 sensitivity regime should be implemented in India to improve the air quality, especially for secondary pollutants such as O 3 .
Data availability. The datasets used in the study can be accessed from websites listed in the references or by contacting the corresponding authors (peng.ce.wang@polyu.edu.hk and zhanghl@fudan.edu.cn).
Author contributions. MZ conducted the modeling and led the writing of the paper. AK carried out the data collection and initial analysis. SZ, JS, and JM assisted with the data analysis. MX and SHK assisted with the interpretation of the results and the writing of the paper. HZ and PW designed the study, discussed the results, and edited the paper.