Substantial changes in gaseous pollutants and chemical compositions in fine particles in the North China Plain during the COVID-19 lockdown period: anthropogenic vs. meteorological influences

The rapid response to the COVID-19 pandemic led to unprecedented decreases in economic activities, thereby reducing the pollutant emissions. A random forest (RF) model was applied to determine the respective contributions of meteorology and anthropogenic emissions to the changes in air quality. The result suggested that the strict lockdown measures significantly decreased primary components such as Cr (−67 %) and Fe (−61 %) in PM2.5 (p < 0.01), whereas the higher relative humidity (RH) and NH3 level and the lower air temperature (T ) remarkably enhanced the production of secondary aerosol, including SO2− 4 (29 %), NO − 3 (29 %), and NH+4 (21 %) (p < 0.05). The positive matrix factorization (PMF) result suggested that the contribution ratios of secondary formation (SF), industrial process (IP), biomass burning (BB), coal combustion (CC), and road dust (RD) changed from 36 %, 27 %, 21 %, 12 %, and 4 % before the COVID-19 outbreak to 44 %, 20 %, 20 %, 9 %, and 7 %, respectively. The rapid increase in the contribution ratio derived from SF to PM2.5 implied that the intermittent haze events during the COVID-19 period were characterized by secondary aerosol pollution, which was mainly contributed by the unfavorable meteorological conditions and high NH3 level.


Introduction
In December 2019, a cluster of pneumonia cases with unknown etiology were first reported in Wuhan and quickly spread around the world (F. . The continuous global outbreak of the coronavirus disease , declared as a public health emergency of international concern by the World Health Organization, resulted in unprecedented public health responses in many countries including lockdown, travel restrictions, and quarantines (Griffiths and 8678 R. Li et al.: The impact of COVID-19 lockdown on air quality in Tangshan, China local air quality. S. Feng et al. (2020) confirmed that the COVID-19 lockdown led to more than a 70 % reduction in NO x emissions in many large cities over China. Correspondingly, the concentrations of PM 2.5 and NO 2 decreased by 35 % and 60 %, respectively (Shi and Brasseur, 2020). The natural experiment provided an unprecedented opportunity to explore the potential for emission reduction and the corresponding response of air quality.
A growing body of studies assessed the response of PM 2.5 and gaseous pollutants to the COVID-19 lockdown and found that these stringent restrictions resulted in significant decreases in these pollutant (e.g., PM 2.5 , NO 2 , and CO) concentrations (Miyazaki et al., 2020;Marlier et al., 2020). However, some haze events still occurred during this period, especially in East China.  employed the chemical transport models (CTMs) to infer that these extraordinary findings might be attributable to enhanced secondary pollution. Understanding that the formation mechanism of puzzle haze events depending on CTMs alone might not be very robust, it was highly imperative to perform more field observations to analyze the temporal variations in chemical compositions, especially the secondary ions (e.g., SO 2− 4 , NO − 3 ) in PM 2.5 before and after the COVID-19 outbreak, and then to validate these inferences.
To date, only several field observations analyzed the temporal variations in chemical components in fine particles during the COVID-19 lockdown period. Chang et al. (2020) observed that a remarkably enhanced nitrate formation in the Yangtze River Delta (YRD) counteracted the decreases in primary components in fine particles, which was in good agreement with the modeling result drawn by . In contrast, Xu et al. (2020) found that the marked decreases in fine particle concentrations in Lanzhou during the COVID-19 lockdown period was mainly contributed by the lower production rate for secondary aerosols. Under the condition of similar emission control measures, the polarized conclusion might be associated with the local meteorology. He et al. (2017) demonstrated that meteorology might explain more than 70 % of variances in daily average pollutant levels over China during 2014-2015. Additionally, X.  also revealed that the release of primary pollutants and the generation of reactive semi-volatile products partitioned between gas and aerosol phases were strongly dependent on the temperature and relative humidity (RH). Thus, in order to accurately assess the effects of lockdown measures on air quality and to reveal the key driver of the haze paradox, it was necessary to isolate the contribution of meteorology. Unfortunately, to date, the respective contributions of emission and meteorology to chemical compositions in PM 2.5 during the COVID-19 period have not been quantified yet in most pioneering studies (Chang et al., 2020;Xu et al., 2020). Moreover, the comparison of source contributions to chemical compositions between pre-lockdown and post-lockdown were scarcely per-formed. Such knowledge is critical to design effective PM 2.5 mitigation strategies in the near future.
As a heavily industrialized region, the North China Plain (NCP) possesses many energy-intensive industries including coal-fired power plants, non-ferrous smelting industries, textiles, building materials, chemical engineering, and papermaking industries (Ren et al., 2011). Due to these intensive industrial emissions, the NCP suffered from poor air quality and frequent aerosol pollution in the past decades Luo et al., 2017). Nevertheless, these strict lockdown measures during the COVID-19 period inevitably led to dramatic decreases in industrial emissions, and thus a study about the response of chemical compositions to emission reduction in the heavy-pollution city might be more sensible.
Here, we selected the typical industrial city (Tangshan) in the NCP to determine the concentrations of gaseous pollutants and chemical compositions in PM 2.5 during 1 January-31 March 2020 and then to analyze their temporal variations before and after the COVID-19 outbreak. Additionally, a machine-learning approach was applied to separate the contributions of emission reduction and meteorology to the temporal variabilities in chemical compositions and gaseous pollutants. Finally, the source apportionment was performed based on the meteorology-normalized datasets to compare the source difference for these pollutants before and after the COVID-19 lockdown.

Field observation
Hourly gaseous pollutants and PM 2.5 chemical compositions including water-soluble ions and trace elements were measured using online instruments during 1 January-31 March 2020 at a supersite in Tangshan. The supersite is located in a commercial region without short-distance industrial emissions ( Fig. 1). SO 2 , NO 2 , and CO concentrations were determined by the ultraviolet fluorescence analyzer (TEI model 43i from Thermo Fisher Scientific Inc., USA), chemiluminescence trace gas analyzer (TEI model 42i from Thermo Fisher Scientific Inc., USA), and the correlation infrared absorption analyzer (TAPI model 300E, USA) . The mass concentration of PM 2.5 was determined using an oscillating balance analyzer (TH-2000Z, China) . The NH 3 concentration and watersoluble ions including sulfate (SO 2− 4 ), nitrate (NO − 3 ), ammonium (NH + 4 ), sodium ion (Na + ), and chloridion (Cl − ) were monitored with a gas and aerosol collector combined with ion chromatography (GAC-IC; TH-PKU-303, China) Zheng et al., 2019). Organic carbon (OC) and elemental carbon (EC) were measured using an OC-EC analyzer (model RT-4, Sunset Laboratory Inc., Tigard, Oregon, USA). Nine trace elements including Hg, Pb, K, Ca, Cr, Cu, Fe, Ni, and Zn were determined by an online multielement analyzer (model Xact 625, Cooper Environment Service, USA). The quality assurance of SO 2 , NO 2 , CO, and PM 2.5 was conducted based on HJ 630-2011 specifications. For the quality assurance of NH 3 and water-soluble ions, the concentration gradients of anion and cation standard solutions were set based on the pollution levels of target species, and correlation coefficients of the calibration curve must be higher than 0.99. Additionally, a standard sample was collected each day, and the relative standard deviation for the reproducibility test must be less than 5 %. The online device agreed well with the result determined by filter sampling coupled with inductively coupled plasma mass spectrometry (ICP-MS) and inductively coupled plasma atomic emission spectroscopy (ICP-AES).

Deweathered model development
The air pollutants were influenced by the combined effects of meteorological conditions and emissions. In order to quantify the contributions of anthropogenic emissions, the impacts of meteorological conditions should be removed. In our study, a random forest (RF) approach was employed to serve as the site-specific modeling platform . All gaseous pollutants and chemical compositions in PM 2.5 were regarded as the dependent variables. The meteorological parameters including wind speed (WS), wind direction (WD), air temperature (T ), relative humidity (RH), precipitation (Prec), air pressure (P ), and time predictors (year, day of year (DOY), day of week (DOW), hour) served as the independent variables. The original dataset was randomly classified into a training dataset (90 % of input dataset) for developing the RF model, and the remaining one was treated as the test dataset. After the building of the RF model, the deweathered technique was applied to predict the air pollutant level at a specific time point (e.g., 1 January 2020 at 12:00 local time). The differences in original pollutant concentrations and deweathered pollutant concentrations were regarded as the concentrations contributed by meteorology. Some statistical indicators including R 2 value, RMSE, and MAE were regarded as the major criteria to evaluate the modeling performance. In our study, the RF model with the R 2 value lower than 0.50 was treated as the unreliable result and cannot reflect the impacts of emission and meteorology on air pollutants accurately because more than 50 % of the variability in the training model cannot be appropriately explained. After the model evaluation, only the species with cross-validation R 2 values larger than 0.50 were selected to assess the respective contributions of emission and meteorology to their concentrations.

Source apportionment
The positive matrix factorization (PMF 5.0) model version was used to perform the PM 2.5 source apportionment. The deweathered gaseous pollutants and chemical compositions in PM 2.5 were incorporated into the model. The objective of PMF is to resolve the issues of chemical mass balance between measured concentration of each species and its source contributions by decomposing the input matrix into factor contribution and factor profile. The detailed equation is shown in Eqs. (1) and (2). Briefly, the basic principle of PMF is to calculate the least object function Q when the g ik must be a positive definite matrix based on Eq. (2) (Chen et al., 2014;Sharma et al., 2016).
where x ij and e ij represent the concentration and uncertainty in the j th species, respectively; g ik represents the contribution ratio of the kth source to the ith sample; f kj represents the ratio of the j th species in the kth source; and e ij indicates the residual of the j th species in the ith sample. The uncertainties associated with factor profiles were evaluated using three error calculation methods including the bootstraps (BS) method, displacement (DISP) analysis, and the combination method of DISP and BS (BS-DISP). For the BS method, 100 runs were performed, and the result has been believed to be valid since all of the factors showed a mapping of above 90 %. DISP analysis also confirmed that the solution was considered to be stable because the observed drop in the Q value was less than 0.1 %, and no factor swap occurred. For the BS-DISP analysis, the solution has been verified to be useful because the observed drop in the Q value was less than 0.5 %. Furthermore, both of the results from BS and BS-DISP did not suggest any asymmetry or rotational ambiguity for all of the factors (Manousakas et al., 2017;Brown et al., 2015).

The impact of emission reduction in air quality
Although the observed concentrations of air pollutants can be applied to analyze the impact of the COVID-19 lockdown, the role of emission reduction in air quality might not be clearly revealed because the meteorological factors were also important variables influencing the air pollutant concentrations. In order to accurately reflect the response of air quality to emission reduction during the COVID-19 lockdown period, the meteorological conditions were isolated by the machine-learning model. In our study, we developed a random forest model to remove the impacts of meteorological conditions on air pollutants. Based on the results in Figs. S1-S5 in the Supplement, the RF models for most of the species showed better performance because their R 2 values were higher than 0.50, and the slope of all of the fitting curves was also close to the R 2 values. However, some other species such as Ag, Cd, and Mg 2+ showed worse predictive performances, and thus these data cannot be utilized to distinguish the impacts of meteorology and emission on the concentrations of these species. Based on the cross-validation R 2 value, the species with an R 2 value higher than 0.50 were applied to assess the contributions of meteorology and emission to the concentrations. The deweathered concentrations of gaseous pollutants and chemical compositions in PM 2.5 are depicted in Figs. 2-5. Compared with the period before COVID-19, the deweathered NH 3 , SO 2 , CO, and NO 2 concentrations decreased by 27 %, 31 %, 32 %, and 42 % after the COVID-19 lockdown period outbreak, respectively (p < 0.05), while the deweathered 8 h O 3 concentration increased by 80 % (p < 0.01). Meanwhile, the normalized- meteorology NH + 4 , NO − 3 , SO 2− 4 , Cl − , PM 2.5 , and Na + and concentrations decreased by 14 %, 27 %, 35 %, 35 %, 38 %, and 47 %, respectively. For trace elements, deweathered Cu, K, Ni, Ca, Pb, Fe, Cr, and Zn levels reduced by 15 %, 23 %, 27 %, 54 %, 59 %, 61 %, 67 %, and 69 %, respectively (p < 0.05). Nevertheless, the deweathered Hg concentration still kept a stable increase at a rate of 6 % compared with the period before the COVID-19 outbreak (p > 0.05).
The deweathered concentrations for most of the pollutants showed significant decreases after the COVID-19 outbreak compared with the period before COVID- . It was assumed that many cities proposed the lockdown measures, which significantly minimized industrial, transportation, and commercial activities. Among all of the pollutants, the deweathered Zn, Cr, Fe, Pb, and Ca experienced decrease rates of more than 50 % due to the lockdown measures. It was well known that Zn, Cr, and Fe originated mainly from the metallurgical industry (Sun et al., 2018;Zhu et al., 2018), while Pb might be derived from coal-fired power plants Meng et al., 2020). During the COVID-19 outbreak, most of the industries have been shut down, and energy production by coal-fired power plants was reduced by one-third (Chang et al., 2020). Based on the adjustment factor estimated by Doumbia et al. (2021), the contributions of industrial activity and the power sector have decreased by 40 % after the COVID-19 outbreak, which was close to the decrease ratios of Zn, Cr, Fe, and Pb concentrations. It should be noted that the deweathered Ca concentration also decreased by more than 50 %. It was well documented that the Ca was often associated with the dust resuspension . In fact, the Ca was known as one of the most abundant elements in the upper continental crust, which likely originated from the fugitive dust Shen et al., 2016). A reduction of more than 70 % in vehicle transportation and domestic flights facilitated the rapid decrease in Ca concentration (Chang et al., 2020). Although the observed K concentration did not show a marked decrease after the COVID-19 lockdown, the deweathered K level suffered from rapid decrease (−22 %) (p < 0.05). It was widely acknowledged that K was considered to be a key fingerprint of biomass burning (H. Zheng et al., 2020), and thus the result suggested that the open biomass burning was also restricted during the period. Both of the deweathered concentrations of OC (−22 %) and EC (−45 %) also experienced remarkable decreases. In our study, both OC and EC concentrations showed significant correlation with K level (p < 0.05), indicating that the restriction of biomass burning also led to the decreases in OC and EC. Additionally, PM 2.5 and some water-soluble ions including deweathered SO 2− 4 and NO − 3 concentrations experienced marked decreases after the COVID-19 lockdown, which was in good agreement with their gaseous precursors. It might be attributable to the rapid decreases in precursor emissions. B. Zheng et al. (2020) verified that the SO 2 emission in the industrial sector and NO x emission in the transportation sector in Hebei Province have decreased by 19 % and 13 %, respectively. The deweathered Na + concentration showed the rapid decrease after the COVID-19 lockdown, which suggested that the Na + in the PM 2.5 of Tangshan was probably derived from waste incineration rather than sea-salt aerosol (Deshmukh et al., 2016).
Although most pollutant concentrations suffered from remarkable decreases, the decrease ratios of deweathered NH 3 and NH + 4 concentrations after the COVID-19 outbreak were far lower than those of many other gaseous pollutants and water-soluble ions. It was attributable to the fact that ambient NH 3 was mainly sourced from the fertilizer application and livestock, which did not show significant decrease during the COVID-19 period (Kang et al., 2016;Doumbia et al., 2021). Although the transportation volume suffered from dramatic decrease, the contribution of transportation to NH 3 was generally less than 5 % (Kang et al., 2016). Furthermore, the contribution of urban waste sources slightly increased after the COVID-19 outbreak, offsetting the effect of traffic outage (Y. . Additionally, it should be noted that the normalized-meteorology 8 h O 3 and Hg concentrations still maintained a stable increase. T.  have confirmed that uncoordinated decreases in NO x and volatile organic compound (VOC) emis-sions (decrease ratio: NO x > VOCs) dominated the 8 h O 3 increase in urban areas because most urban areas belonged to a VOC-limited region. Moreover, the excessive decrease in PM 2.5 from primary emission significantly increased the HO 2 radical concentration on the surface of aerosol, thereby promoting the O 3 formation (Shi and Brasseur, 2020). The minor increase in deweathered Hg level was attributable to the fact that the coal combustion for domestic heating was not restricted during the COVID-19 lockdown period (Zhou et al., 2018). Based on the updated global anthropogenic emission adjustment factor during COVID-19, the contribution of the residential sector to air pollutants did not decrease after the COVID-19 lockdown (Doumbia et al., 2021).

The role of meteorology and potential chemical reactions on air quality
Compared with the observed values, the deweathered concentrations of most pollutants were significantly reduced. Meanwhile, the deweathered decrease ratios of pollutants  were significantly higher than those of observed values (Fig. 6). The result suggested that the meteorology conditions during the COVID-19 lockdown period were not favorable to the pollutant dispersion, as evidenced by some recent studies (Chang et al., 2020;. In our study, six meteorological parameters including WS, WD, T , RH, Prec, and P have been integrated into the random forest model to assess the response of each species to different meteorological variables. The variable importance of each meteorological parameter to all of the species is shown in Figs. 7-10. Among all of the gaseous pollutants, the meteorological conditions played significantly positive roles in NH 3 (62 %) and 8 h O 3 concentrations (80 %) (Fig. 6). As shown in Fig. 7, T was the most important factor for the rapid elevation of NH 3 concentration after the COVID-19 lockdown. It was assumed that the higher T enhanced the emissions of NH 3 from soil and urban wastes and promoted the volatilization of NH 3 from aerosol NH + 4 pools (Y. . In our study, the hourly mean air temperature increased from 0 • before the COVID-19 outbreak to 5 • after the COVID-19 lockdown, which strongly supported the inference. For 8 h O 3 concentration (Fig. 7), T was also treated as the most important variable. On the one hand, the higher T generally enhanced biogenic isoprene emissions, which was the most abundant biogenic VOC and showed the highest ozone formation potential . On the other hand, high T often increased chemical reaction rates and accelerated the O 3 formation . Additionally, WS also played an important role in the 8 h O 3 concentration. Shi et al. (2020) have demonstrated that weaker winds often slowed down the advection and convection of NO x and VOCs, which was beneficial to O 3 formation.
Moreover, the contributions of meteorological conditions to some secondary ions (e.g., SO 2− 4 (29 %), NO − 3 (29 %), and NH + 4 (21 %)) were remarkably higher than those to other ions and some trace elements, suggesting that the chemical reactions and formation pathways of these species were more sensitive to meteorological variations. Deshmukh et al. (2016) confirmed that the high RH promoted the aqueous-phase oxidation of SO 2 and the production of sulfate. Tian et al. (2019) also demonstrated that RH-dependent heterogeneous reactions significantly contributed to the sulfate generation, and the high RH enhanced gas-to aqueousphase dissolution of NH 3 and HNO 3 . These pioneering experiments suggested that secondary aerosols were often formed under the condition of high RH. Very recently, Chang et al. (2020) observed that the nitrate concentration in the YRD experienced unusual increase during the COVID-19 period, while Xu et al. (2020) obtained the opposite result in Lanzhou. It was assumed that the persistent increase in T and decrease in RH in Lanzhou during this period was not beneficial to the generation of secondary aerosol, while the high RH in the YRD significantly elevated local nitrate level. Although air temperature in Tangshan suffered from an increase after the COVID-19 lockdown, RH displayed a rapid increase from 47 % to 57 % during this period. Moreover, the increased O 3 could promote the secondary aerosol formation and partially offset the decreased PM 2.5 compositions triggered by the primary emission reduction . Similar to secondary ions, both OC and EC were also sensitive to RH. It was supposed that high RH could increase the secondary organic aerosol (SOA) levels, which accounted for the major fraction of OC (H. Zheng et al., 2020).
In addition, some trace elements such as Fe, Ni, and Cr were also significantly affected by the meteorological conditions. As shown in Fig. 9, these element concentrations were mainly sensitive to WD. It was assumed that the neighboring industrial points including cement plants and coal-fired power plants could influence the concentrations of trace elements via long-or short-range transport, which was strongly dependent on WD. Following WD, RH was also an important factor for the variation in these trace elements. Under the condition of high RH, Fe and Cr could catalyze the heterogeneous generation of sulfate and nitrate on the mineral or soot surface (Hu et al., 2015).
Unlike the trace elements, water-soluble ions and OC were less sensitive to WD. Major water-soluble ions in PM 2.5 including SO 2− 4 , NO − 3 , and NH + 4 were mainly derived from secondary formation rather than the direct emission (J. Feng et al., 2020;, and thus they were not very sensitive to WD.

The enhanced secondary aerosol formation during the COVID-19 lockdown period
The deweathered chemical compositions suggested that the sulfate and nitrate chemistry changed slightly after the COVID-19 outbreak. The oxidation ratio of sulfate (SOR; the ratio of sulfate concentration and the sum of sulfate and SO 2 concentrations) decreased from 0.26 to 0.22, while the oxidation ratio of nitrate (NOR; the ratio of nitrate concentration and the sum of nitrate and NO 2 concentrations) increased from 0.22 to 0.25 (Table 1). The decreased SOR after the COVID-19 outbreak indicated that the decrease rate of sulfate is higher than that of SO 2 . In contrast, the increased NOR during the COVID-19 lockdown period revealed that the decrease rate of nitrate is lower than that of NO 2 . The increased NOR after the COVID-19 outbreak suggested the consecutive nitrate production, though the NO x emission experienced tremendous reduction, which was in good agreement with the result observed by Chang et al. (2020). It was assumed that the persistently higher observed NH 3 concentration during this period promoted the ammonium nitrate formation though the lower NO x emission (Y. , which also partially explained the abnormal increases in observed concentrations of secondary ions after the COVID-19 outbreak. In general, NH 3 firstly tends to react with H 2 SO 4 to form ammonium sulfate, and then the excess NH 3 participated in the reaction with HNO 3 . However, sulfate concentration suffered from a more dramatic decrease compared with SO 2 , which might be associated with the aerosol acidity during the COVID-19 lockdown period. The ratio of NH + 4 and the sum of SO 2− 4 , NO − 3 , and Cl − , named C/A, was regarded as an indicator to reflect the aerosol acidity. In our study, the C/A value decreased from 0.33 to 0.28 after the COVID-19 outbreak, implicating that the aerosol acidity even showed slight increase during the COVID-19 lockdown period. It was well known that the higher aerosol acidity might prohibit the conversion from SO 2 to sulfate Shao et al., 2019), which yielded the lower SOR.

The impact of the COVID-19 lockdown on source apportionment
The emission control measures inevitably triggered the variation in source apportionment Meng et al., 2020). In the present study, positive matrix factorization (PMF 5.0) was employed to identify the major sources of PM 2.5 in Tangshan before and after the COVID-19 outbreak.
. About three-to nine-factor solutions were examined, and a five-factor solution obtained the lowest Q (robust) and Q (true) values. Additionally, the PMF analysis and error diagnostics also suggested that the result was robust (Tables S2-S4 in the Supplement).

SOR NOR C/A
The source apportionment profiles pre-COVID and post-COVID resolved by PMF are depicted in Fig. 11. For pre-COVID, the first factor contributed 36 % to the total species. The factor was characterized by high levels of NH + 4 (41 %), SO 2− Figure 7. Relative importance of the predictors for the prediction of gaseous pollutants. The match in the figure denotes the variable importance in RF models for various species. DOY, WD, P , RH, Hour, T , DOW, WS, Prec, and Year represent day of year, wind direction, air pressure, relative humidity, hour of the day, air temperature, day of week, wind speed, precipitation, and study year.  erally produced by oxidation of SO 2 and NO x , respectively. The NH + 4 was often formed through the heterogeneous reaction of NH 3 and sulfate or HNO 3 . Thus, the factor was regarded as the secondary formation (SF). The second factor was characterized by high loadings of Zn (47 %), Cr (42 %), Fe (42 %), and Pb (31 %). Cr and Fe were mainly originated from fuel combustion and the metallurgical industry such as chrome plating and steel production (J. , while Pb and Zn were derived from the roasting, sintering, and smelting process for the extraction of Pb and Zn ores (Wu et al., 2012). Therefore, the factor 2 was treated as the industrial process (IP) source. The predominant species in factor 3 included Na + (42 %), K (40 %), OC (35 %), and EC (33 %). K was often regarded as the fingerprint of biomass burning (BB) Zheng et al., 2019), whereas the Na + was generally regarded as the tracer of waste incineration (Alam et al., 2019;Durlak et al., 1997). Hence, factor 3 was treated as the BB source. Tangshan suffered from remarkable increasing usage of biomass fuels for domestic heating in winter, which promoted the emissions of K and Na + . The most abundant species in factor 4 were Hg (75 %), Pb (68 %), K (36 %), Cu (35 %), Cl − (33 %), and SO 2− 4 (27 %). Pb, Hg, and Cu were typical marker elements for coal combustion, and around 56 % of Pb and 47 % of Hg were released from coal combustion (Cheng et al., 2015;Zhu et al., 2020). In northern China, coal-based domestic heating was one of the most important sectors of coal consumption (K. . Dai et al. (2019) also verified that the residential coal combustion was a major source of primary sulfate. Thus, factor 4 was regarded as the coal combustion (CC) source. The last factor was distinguished by high loadings of Fe (46 %), Ni (45 %), and Ca (38 %). Fe and Ca were main elements enriched in upper crust, and Ni was enriched in brake wear and tire wear Figure 10. Relative importance of the predictors for the prediction of OC and EC in PM 2.5 . The match in the figure denotes the variable importance in RF models for various species. DOY, WD, P , RH, Hour, T , DOW, WS, Prec, and Year represent day of year, wind direction, air pressure, relative humidity, hour of the day, air temperature, day of week, wind speed, precipitation, and study year. Figure 11. The comparison of source apportionment for PM 2.5 chemical compositions before (a) and after (b) the COVID-19 outbreak. In our study, five major sources were distinguished based on the PMF model. The color bar denotes the contributions of these sources to each species. SF, IP, BB, CC, and RD represent secondary formation, industrial process, biomass burning, coal combustion, and road dust, respectively.
After the COVID-19 outbreak, the chemical compositions in PM 2.5 were also classified into five sources including SF, IP, BB, CC, and RD. However, the contribution ratios of these sources varied greatly after the implementation of serious lockdown measures. The contribution ratio of IP experienced the largest decrease from 27 % to 20 %, whereas the apportionment of SF showed a marked increase from 36 % to 44 %. The contributions of the other three sources only suffered from slight variations. The rapid decrease in IP contribution might be associated with the shutdown of many industries during the COVID-19 period , while the obvious increase in SF contribution was attributable to more heterogeneous or aqueous reactions of precursors (Chang et al., 2020). For nearly all of the species, the contribution ratios of IP suffered from remarkable decreases after the COVID-19 outbreak. Since the COVID-19 lockdown, the contribution ratios of SF to SO 2− 4 , NO − 3 , and NH + 4 increased from 35 %, 33 %, and 41 % to 48 %, 44 %, and 52 %, respectively. However, the contribution ratios of SF for other species remained relatively stable. It was assumed that SO 2− 4 , NO − 3 , and NH + 4 were mainly produced from secondary formation of precursors (Jiang et al., 2019;Yao et al., 2020), while other species, especially the trace elements, were mainly derived from the primary emission (Y. . Although the COVID-19 pandemic led to the shutdown of many coal-fired power plants and industries and decreased the CC emissions from these sectors (Kraemer et al., 2020), the government-enforced home order might increase the electricity consumption (Venter et al., 2020), which offset the decreases in CC contributions to industrial activities. Therefore, the contribution ratios of CC did not experience dramatic variation after the COVID-19 outbreak.

Conclusions and implications
The lockdown measures led to the shutdown of many industries, in turn resulting in significant decreases in primary components in PM 2.5 . We employed an RF model to determine the respective contributions of meteorology and emission reduction to the variations in gaseous pollutants and PM 2.5 chemical compositions during the COVID-19 lockdown period. The deweathered levels of some trace elements (e.g., Pb (−59 %), Zn (−69 %)) derived from industrial emissions experienced decrease rates of more than 50 % due to the stringent lockdown measures. However, the higher relative humidity (RH) and lower air temperature (T ) significantly prohibited the decreases in water-soluble ion concentrations because they were beneficial to the heterogeneous or aqueous reaction of sulfate and nitrate. Trace elements were very sensitive to wind direction (WD) due to the long-range transport of anthropogenic emissions. Additionally, the contributions of secondary formation to PM 2.5 increased from 36 % to 44 % after the COVID-19 outbreak. The finding also explained that the opposite change trends of the secondary aerosols in East and West China found by previous studies were attributable not only to the large difference in meteorological conditions but also the discrepancy of NH 3 concentration.
In future work, it will be necessary to seek multi-pollutant (e.g., VOC, NO x ) emission control measures to reduce the concentrations of primary and secondary components simultaneously since adverse meteorological conditions coupled with slightly higher oxidation capability especially, in winter, still caused haze formation. Our results also highlight that more NH 3 emission control measures are urgently needed because the excess NH 3 could exacerbate the generation of secondary aerosols. Besides, the generation of primary pollutants was very sensitive to RH and WD. Thus, the primary pollutant emissions from the industries in the upwind direction should be strictly restricted.
In addition, the present study still suffered from some uncertainties. At first, only six meteorological factors were incorporated into the RF model to quantify the contributions of emission and meteorology of air pollutants. In particular, the missing solar radiation could affect the accuracy of 8 h O 3 estimation. Besides, solar radiation could change the concentrations of hydroxyl radicals, thereby affecting the NO − 3 formation. In future work, the solar radiation should be integrated into the model. In addition, some temporal indicators such as hour and DOY were applied to reflect the COVID-19 lockdown intensity because an hourly emission inventory during this period was not available, which should be integrated into the RF model after the development of a real-time emission inventory.
Data availability. The meteorological data are available at http: //data.cma.cn/ (China Meteorological Administration, 2021). The website can be browsed in English http://data.cma.cn/en.
Author contributions. MP and CW collected data and provided technical support. HF and JC designed the study. RL wrote the manuscript. YZ analyzed the data.
Competing interests. The authors declare that they have no conflict of interest.