Measurement report: The 4-year variability and influence of the Winter Olympics and other special events on air quality in urban Beijing during wintertime

Comprehensive measurements are vital to obtain big enough datasets for better understanding the complex atmosphere and further improving the air quality. To investigate the 4-year variation of air quality and the influences of special events (Beijing Winter Olympics, COVID lockdown and Chinese New Year) on it during the wintertime in polluted urban air, we conducted comprehensive observations in Beijing, China, during 1 January–20 February, in the years from 2019 to 2022. The mass concentration of PM2.5 and its composition (organics, nitrate, sulfate, ammonium, chloride and black carbon) and the number size distributions of particles (down to ∼ 1 nm) and ions, gaseous pollutants (CO, NOx , SO2, O3) and condensable vapors (sulfuric acid and oxygenated organic molecules), as well as meteorological parameters, were simultaneously measured. The days before 22 January without any special events in each year were selected to investigate the 4-year variability of air quality. We found that the concentrations of CO, NOx , total oxygenated organic molecules (OOMs), total PM2.5, organics, chloride and black carbon and the number concentration of sub-3 nm particles (N1.3−3) showed similar variations, decreasing from 2019 to 2021 and then increasing in 2022. For SO2, however, its concentration decreased year by year due to the significant emission reduction, further leading to the decrease of gaseous sulfuric acid and particulate sulfate from 2019 to 2022. O3 concentration showed an opposite 4-year variation compared with NOx . Meanwhile, both the oxygen and nitrogen contents of oxygenated organic molecules increased year by year, implying that not only the oxidation state of those compounds increased, but also NOx was involved more efficiently in their formation processes. With higher sulfuric acid concentrations and new particle formation (NPF) frequencies in 2021 than in 2022, and with the lowest concentrations of background aerosols and the lowest ambient temperatures in 2021, N1.3−3 was still the lowest in 2021. Unlike N1.3−3, the ion concentrations in both 0.8–2 and 2–4 nm size ranges were higher in 2021 than in the other years. Then, the days after 4 February were chosen to explore the influence of special events. The non-event days within this date Published by Copernicus Publications on behalf of the European Geosciences Union. 6664 Y. Guo et al.: Influence of special events on air quality in urban Beijing during wintertime range in 2019 and 2021 were chosen as the reference period. Due to the favorable meteorological conditions together with reductions in anthropogenic emissions, there were basically no haze events during the Olympics. Therefore, CO, NOx , SO2, total OOMs, accumulation-mode particles (N100−1000), and total PM2.5 and its composition were much lower, while ion concentrations were much higher compared with the reference period. Although there was also emission reduction during COVID, especially for NOx , the enhancement of secondary inorganic aerosol formation, together with unfavorable meteorological conditions, caused severe haze events during this period. Hence, CO, total OOMs and all PM2.5 compositions during COVID increased dramatically compared with the reference period. Influenced by SO2, condensation sink and sunlight, sulfuric acid concentration was found to be comparable between the Olympics and the reference period but was lower during COVID and Chinese New Year. Additionally, N1.3−3 was almost at the same level during different periods, indicating that the special events only had little impact on the NPF processes. These results provide useful information to the development of more targeted pollution control plans.


Introduction
Beijing and the surrounding areas suffer from poor air quality associated with high particulate matter (PM) concentrations, especially during wintertime. To improve the air quality, extensive studies have been carried out to understand the sources, formation and evolution of air pollutants in Beijing Sun et al., 2014;Zheng et al., 2015;Shang et al., 2020;Cheng et al., 2016). Primary emission sources, including traffic, cooking, coal combustion and biomass burning, have been identified (Du et al., 2022b;Cai et al., 2020;. Of primary emissions, those from combustions were found to be the major contributors to haze, especially during the heating season (Cheng et al., 2013; Y. L. Sun et al., 2013Sun et al., , 2014. Compared with primary emissions, secondary aerosols have been shown to play a more important role in haze formation, with contributions higher than 60 % of PM 1 mass concentration in Beijing (Shang et al., 2020;Huang et al., 2014;. Secondary organic aerosols (SOAs), nitrate (NO 3 ) and sulfate (SO 4 ) originate from their gaseous precursors, e.g., volatile organic compounds (VOCs), nitrogen oxides (NO x ) and sulfur dioxide (SO 2 ). VOCs, NO x and SO 2 further produce oxygenated organic molecules (OOMs), nitric acid (HNO 3 ) and sulfuric acid (H 2 SO 4 ), respectively, via homogeneous reactions, which can contribute to particle growth by condensation (Bianchi et al., 2019;Trostl et al., 2016;Yue et al., 2010). Most recently, Nie et al. (2022) reported that the condensation of OOMs contributed to > 30 % of SOA in Beijing. Meanwhile, heterogeneous reactions were found to promote the formation of secondary aerosols. Dissolved SO 2 can be oxidized into SO 4 by atmospheric oxidants, e.g., H 2 O 2 , O 3 and NO 2 . Heterogeneous reactions between HNO 3 and NH 3 in daytime and hydrolysis of N 2 O 5 during nighttime are the main pathways forming NO 3 (Xue et al., 2014;. Furthermore, liquid-phase-related SOAs were found to increase during severe haze when relatively humidity was high Zhao et al., 2019). These results highlight the importance of reducing not only the primary particulate emissions but also the anthropogenic gaseous precursors to suppress secondary aerosol formation on regional scales Du et al., 2021).
Regional transport and meteorological conditions also play roles in the evolution of haze episodes. Air masses coming from the south of Beijing often bring pollutants to this area and hence favor the rapid buildup of particle masses (Ma et al., 2017;Zheng et al., 2015;Sun et al., 2015). It has been found that 40 % of PM 2.5 could originate from regional transport in an annual scale (Ge et al., 2018), and under very unfavorable meteorological conditions, it could contribute even up to 80 % of the total particle mass during a single haze episode (Sun et al., 2016). In addition, haze events occur typically under stagnant conditions . Under such circumstances, high aerosol concentrations tend to delay the onset of precipitation , and thus, the relative humidity (RH) within the boundary layer is often rather high, which further promotes the heterogeneous reactions producing secondary aerosols (Y. Sun et al., 2013). Additionally, results have shown that when RH increased above 60 %, particles were found to be in a liquid state, thereby accelerating secondary formation .
The air quality in Beijing has improved substantially over the past decade, especially after introducing the "Action Plan for Air Pollution Prevention and Control" in 2013 (Wang et al., 2020b;Lu et al., 2020;. However, the particulate pollution in Beijing still exceeds the national air quality standards . In addition, although the concentration of SO 2 decreased significantly, the level of NO x remains still high, which results in NO 3 becoming the main contributor of secondary inorganic aerosols (SIAs) (Xie et al., 2020). Furthermore, concurrent with the decreasing PM concentrations, ozone is rapidly becoming a yearround air pollution problem in China .
In addition to long-term actions, Beijing has imposed strict short-term emission reductions during highly visible international events, such as the Beijing Summer Olympics in 2008 Shou-bin et al., 2009;Okuda et al., 2011;Schleicher et al., 2012), the Asia-Pacific Economic Cooperation (APEC) summit in 2014 Sun et al., 2016) and the Victory Day parade in 2015 . More recently, after the start of the COVID-19 outbreak, the Chinese government carried out strong restrictions (COVID lockdown) in order to prevent the spread of the virus. Reduced anthropogenic activities resulted in decreased emissions, yet severe pollution episodes were still observed (Le et al., 2020). In preparation for the Beijing Winter Olympics, the Chinese government authorized necessary actions, successfully improving air quality during this period. Additionally, annual Chinese New Year celebrations are typically associated with reductions in anthropogenic emissions during the 7 d holiday (Tan et al., 2009;Lin and McElroy, 2011).
All these special events mentioned above have different characteristics in regard to emission reductions. During the COVID lockdown, traffic was reduced, but cooking emissions and industrial activities continued at least partially (Shi and Brasseur, 2020). During the Beijing Winter Olympics, heavy industrial activities were regulated to improve the air quality. During the Chinese New Year, a large proportion of Beijing migrant workers and students from other cities leave the city and return to their hometowns. Therefore, anthropogenic emissions from, e.g., traffic, cooking and industry are reduced substantially. Thus, a comparison between COVID lockdown, the Beijing Winter Olympics and Chinese New Year can provide a unique chance to investigate the response of air pollutants to different emission reduction actions.
In this study, we examined the wintertime (1 January to 20 February) air quality in urban Beijing during a 4-year period from 2019 to 2022. We utilized comprehensive observations of air pollutants in both gas phase (including carbon monoxide, sulfur dioxide, nitrogen oxides, ozone, sulfuric acid, and oxygenated organic molecules) and particle phase (including number size distributions of particles and ions and the mass concentrations of PM 2.5 and its composition), as well as of meteorological conditions (including temperature, RH, UVB radiation, wind speed and direction, and boundary layer height) from the Aerosol and Haze Laboratory operated by the Beijing University of Chemical Technology (AHL-BUCT; Liu et al., 2020). The objectives of this study are (1) to investigate the 4-year variabilities of different pollutants and to understand their connections with emissions and meteorological conditions and (2) to examine the characteristics of the atmospheric pollution cocktail during different short-term special events (Beijing Winter Olympics, COVID lockdown and Chinese New Year periods) associated with substantial emission reductions. This study provides helpful information on pollution characteristics of urban Beijing both on a long-term scale and for short-term special periods with different emission reductions, which can give guidance for making targeted and sustainable emission control plans.

Site description
Beijing is located in the northwestern part of the North China Plain. There are two mountains and a sea nearby. The Taihang Mountain is to the west, and the Yanshan is to the northwest of Beijing, and the nearest sea, the Bohai Sea, is ∼ 150 km to the east. On the south of Beijing, there are several megacities, such as Baoding and Shijiazhuang. Our measurements were conducted in downtown Beijing, the west campus of Beijing University of Chemical Technology (39.95 • N, 116.31 • E). The instruments are located on the fifth floor of a teaching building, which is ∼ 20 m above the ground level. The station can be considered as a representative urban site, and more detailed information can be found elsewhere Yan et al., 2021;Guo et al., 2021;Du et al., 2022a).
The concentrations of neutral sulfuric acid (SA) and oxygenated organic molecules (OOMs) were measured by a chemical ionization atmospheric pressure interface longtime-of-flight (CI-APi-TOF; Aerodyne Research, Inc.) mass spectrometers using nitrate (NO − 3 ) as the reagent ion. Detailed configurations and working parameters of this nitrate chemical ionization mass spectrometer (CIMS) can be found elsewhere Guo et al., 2022b). The calibration of SA was implemented by introducing a known amount of gaseous SA produced by the reaction of SO 2 and OH radical formed by UV photolysis of water vapor, which is similar to the method in previous literature (Kürten et al., 2012). The calibration factor of sulfuric acid was 6.07-7.47 × 10 9 cm −3 / (normalized cps) from 2019 to 2022. For the quantification of OOMs, a mass-dependent transmission method was used, and details of this approach are described elsewhere . The concentration of each molecular OOM species is calculated as follows: of detected signal of that OOM, either as a neutral molecule or as a de-protonated ion (OOM-H) − ; the denominator is the sum of all measured reagent ions; C is the calibration factor of H 2 SO 4 ; and T OOM is the relative transmission coefficient. The mass concentration of PM 2.5 was measured with a tapered element oscillating microbalance dichotomous ambient particulate monitor (TEOM 1405-DF; Thermo Fisher Scientific Inc, USA). Black carbon (BC) in PM 2.5 was measured using a seven-wavelength aethalometer (AE33; Magee Scientific Crop.) (Drinovec et al., 2015). The non-refractory chemical compositions of fine particles (NR-PM 2.5 ), including organics (Org), sulfate (SO 4 ), nitrate (NO 3 ), ammonium (NH 4 ), and chloride (Chl), were measured using an online time-of-flight aerosol chemical speciation monitor (ToF-ACSM; Aerodyne Research Inc. U.S.) equipped with a PM 2.5 aerodynamic lens and a standard vaporizer (Jayne et al., 2000;Drewnick et al., 2005;Cai et al., 2020).

Calculation of averaged oxygen and nitrogen numbers of total OOMs
The fraction-weighted oxygen and nitrogen numbers of OOMs are calculated based on the following equations: where n is the number of OOM molecules, i is one specific OOM molecule, nO i (nN i ) is the number of oxygen (nitrogen) in OOM i , fraction i is the number fraction that OOM i takes and t is one certain moment. Then, nO and nN could reflect the averaged oxygen and nitrogen number of total OOMs at each moment. Here, we use nO to generally reflect the oxidation state of total OOMs, and the higher the nO, the more oxidized the total OOMs.

Division of different periods
In this study, we focus on the 4-year variability of air quality in winter in Beijing (abbreviated as 4-year variability in the following text), as well as the influence of special events, including the Beijing Winter Olympics (Olympics), COVID lockdown (COVID) and Chinese New Year (CNY) during the 4 years from 2019 to 2022. Thus, the whole winter period from 1 January to 20 February was divided into two separate ones ( Fig. S1 and Table 1): 1. Periods of 4-year variability. When the 4-year variability is investigated, there should be no disturbance from any of the special events, and the chosen days should be from the same periods of different years. Therefore, the period of 4-year variability covered days from 1 to 22 January of each year, lasting for 22 d.

2.
Periods of special events. To ensure comparable durations for special events, days ranging from 4 to 20 February were chosen as the special event periods. The days within this date range in 2020 and 2022 are referred to as COVID and Olympics periods, respectively. Since CNY holidays in 2022 interfered with the Olympics, only CNY holidays in 2019 and 2021 were used and combined together as the CNY period. Finally, for a better understanding of the special event effect, a reference period including the non-event days in 2019 and 2021 was chosen as the "base" period.

Results and discussion
3.1 The 4-year variability of air quality

Meteorological conditions
Meteorological conditions are strongly interlinked with air quality. Therefore, we first focus on the 4-year changes of meteorological parameters. The periods covered the days from 1 to 22 January of each year. As shown in Fig. 1, the local wind distributions in 2019 and 2020 were quite similar, following the typical diurnal wind pattern in Beijing induced by the mountain-valley breeze. From midnight to noon, the wind mainly blew from the northwest, and the wind speed (WS) was quite low (median ∼ 0.5 m s −1 ). From afternoon, however, the wind direction turned to a southerly wind. In 2021, the wind was mostly from the west and northwest without a clear diurnal variation. In 2022, the wind blew mostly from the southeast from midnight to noon but came from the northeast from the afternoon to midnight. Usually, the air masses coming from the south or east possess higher temperatures and bring more water vapor as well as pollutants, whereas air masses coming from the north or west are typically colder with a lower RH and lower pollutant concentrations Zhong et al., 2018). The year of 2021 witnessed the lowest median temperature (−2.0 • C,    Fig. 2b) compared with the other 3 years, which might be partly due to the wind coming mainly from the west and northwest throughout the day. Meanwhile, the overall wind speed in 2021 was also the highest (median 0.75 m s −1 , Figs. 1 and 2d). As the typical diurnal pattern of wind direction in Beijing is induced by the mountain-valley breeze, which usually only occurs during rather stagnant synoptic conditions, the lack of diurnal pattern of wind direction and the higher wind speeds in 2021 suggest that the synoptic meteorological conditions might have been stronger during this period in 2021 than in the other years. Similarly, BLH in 2021 was also the highest throughout the day ( Fig. 2e and f). During the other 3 years, the overall values of the temperature, wind speed and BLH were comparable to each other, whereas the RH was higher in 2020 and 2022, which could partly be attributed to more southern and eastern winds in these 2 years. In addition, differences in cloud conditions could have a minor effect on UVB together with varying pollution conditions. Higher values of RH may facilitate the hygroscopic growth of aerosol particles (Hodas et al., 2014;Cheng et al., 2016) and change the particle phase state, thereby likely promot- ing secondary aerosol formation via heterogeneous reactions (Hung et al., 2016;Chao et al., 2020;Shiraiwa et al., 2017;Sun et al., 2018). Furthermore, precipitation can increase the wet deposition of almost all gaseous and particulate species, influencing their concentrations in the atmosphere. However, since there was no reliable measurement of precipitation at our site, detailed conclusions on this matter cannot be drawn. It should be mentioned that the intensity of snow precipitation in winter in Beijing is generally small, except during special case(s), so the overall effects of precipitation in winter in Beijing are not as large as those in summer in Beijing or in other humid areas.

Trace gases
The 4-year variations of CO and NO 2 were generally opposite to the variation of BLH, which decreased from 2019 to 2021 and rebounded in 2022, resulting in the lowest concentrations in 2021 ( Fig. 3a and c; median values are 308 and 12.8 ppb for CO and NO 2 , respectively). Although the relative changes of them were larger than BLH, the variations of them should be partially attributed to the change in BLH or the change in atmospheric diffusion capacity. Fur-ther, the winds in 2021 came less frequently from polluted sectors, which might also play a role. For NO, it decreased from 2019 (median 10.8 ppb) to 2021 (median 5.5 ppb) and was around the same level in 2021 and 2022 (Fig. 3b). NO is mainly a primary emission, and in urban Beijing, it mainly comes from vehicle emissions. Although the BLH was lower and the wind was more polluted in 2022 than in 2021, NO levels between these 2 years are nearly the same, suggesting that traffic emissions might have reduced in 2022. Despite the influence of different wind direction and speed, SO 2 concentration monotonically decreased from 2019 (median 3.06 ppb) to 2022 (median 0.68 ppb). In addition, this trend is quite different from that of BLH, indicating that the drop of SO 2 is not caused by the variations of meteorology conditions but by the decline in SO 2 emissions. Our result is in line with the long-term variation of SO 2 in China, where SO 2 decreased from 24.8 ppb in 2013 to 8.0 ppb in 2017 for the Beijing-Tianjin-Hebei area after the clean air action in 2013 (Y. Wang et al., 2020b). The 4-year variation of O 3 was generally opposite to CO and NO 2 , and its concentration was the highest in 2021 (15.4 ppb, Fig. 3d). This highest O 3 could be partly associated with the weakened titration effect of NO. Meanwhile, the level of O 3 decreased with the increase of PM 2.5 (Fig. S6), suggesting that O 3 formation might be suppressed during haze episodes within the studied period.

Condensable vapors
SA is a key contributor to the formation and initial growth of aerosol particles (Sipilä et al., 2010;Paasonen et al., 2010;Kirkby et al., 2011;Yao et al., 2018;Yan et al., 2021), and OOMs are the main drivers of the particle growth and formation of secondary organic aerosol (SOA) (Ehn et al., 2014;Nie et al., 2022). Therefore, for further investigation of number size distributions and mass concentrations of particles, the 4-year changes of SA and OOMs were first analyzed. As shown in Fig. 4, the concentrations of sulfuric acid monomer (SA1; 4.5-9.6×10 5 cm −3 , Fig. 4a) and dimer (SA2; 3.1×10 3 -2.8×10 4 cm −3 , Fig. 4b) generally decreased from 2019 to 2022. This decreasing trend is similar to that of SO 2 , and therefore, the decline of SA1 and SA2 was probably mainly caused by the drop of the precursor SO 2 . Additionally, SA2 in 2021 was higher than 2020 and 2022. Figure 5 shows that although the data points in 2021 stand out in the SA2 vs. SA1 plot, they lie together with the data points of the other 3 years in the SA2 vs.
The 4-year variation of the total OOM concentration was generally opposite to that of the BLH. Specifically, the median OOM concentration in 2019 (2.4×10 7 cm −3 ) was ∼ 8.6 times that in 2021 (2.8 × 10 6 cm −3 ) (Fig. 4c), while median and peak BLH values in 2019 (360 and ∼ 940 m for median and peak respectively) were only ∼ 0.7-0.8 times those in 2021 (497 and ∼ 1200 m for median and peak respectively) ( Fig. 2e and f). Thus, the reasons behind this drastic yearly change of the OOM concentration cannot be solely attributed to changes in the BLH. Under polluted conditions, regional transportation of lower-volatility OOMs could result in higher concentration of OOMs (Guo et al., 2022b). In addition, accumulation of precursor VOCs, and possibly an enhancement of heterogeneous reactions (Riedel et al., 2015), may promote the chemical formation of OOMs. However, due to the lack of VOC measurements and inadequate understanding of heterogeneous reactions, it is hard to conclude which one is the determining factor for the observed OOM concentration changes. Additionally, the total OOM concentrations at same PM 2.5 levels generally decreased from 2019 to 2022 (Fig. S7), which implies that the OOM concentrations have reduced in recent years. In terms of the OOM composition, Fig. 4d and e show that  both oxygen and nitrogen contents of OOMs increased year by year. None of the meteorology parameters showed such a behavior, so changes in the chemical pathways are likely the main cause. The increase in the oxygen content indicates an enhancement of the overall oxidation state of OOMs, and the increase in the nitrogen content is related to an enhancement of the NO x involvement in OOM formation. Although there are no direct VOC measurements, studies have shown that the VOC levels are continuously decreasing in the region . Therefore, by assuming that the level of oxidants stays the same, the oxidant-to-VOC ratio increases, thus likely leading to more oxidation of each VOC molecule. In addition, long-term observations have shown that the concentration decrease of NO x is slower than that of VOCs . As a result, the NO x : VOC ratio should have increased with time, which may partially explain the increased OOM nitrogen content. Furthermore, the oxygen content decreases as PM 2.5 increases. This is not surprising, as OOMs with lower oxygen contents possess higher volatilities and thus are more easily to be re-evaporated back to the atmosphere and transported along with PM 2.5 . The nitrogen content first increases with an increasing PM 2.5 level and then decreases when PM 2.5 exceeds 150 µg cm −3 . Higher PM 2.5 levels seem to be associated with higher NO x concentrations (Fig. S6), which promotes the formation of organonitrates, leading to a higher nitrogen content. But the reason why the nitrogen content decreases under serve haze remains a puzzle, deserving further studies. In order to evaluate the OOM contribution to SOA formation through condensation, the parameter [ELVOCs + LVOCs] × CS was used as a simple surrogate for OOM condensation flux (Fig. 4f), where ELVOCs and LVOCs stand for extremely low-volatility organic compounds and low-volatility organic compounds respectively.
[ELVOCs + LVOCs] × CS was the highest in 2019 and the lowest in 2021, with 2020 and 2022 lying between and being comparable. This suggests that the SOA formation potential of OOMs was highest in 2019 and lowest in 2021, which is consistent with the organic aerosol concentration in Fig. 10. Figure 6 shows the 4-year changes of number concentrations of particles in different size ranges. The number concentration of sub-3 nm particles (N 1.3−3 ) decreased from 2019 to 2021 but increased slightly from 2021 to 2022. The influencing factors behind this variability are discussed in the following text. Number concentrations in the nucleation mode (3-25 nm, N 3−25 ) and Aitken mode (25-100 nm, N 25−100 ), generally decreased from 2019 to 2022. In particular, N 3−25 was significantly lower in 2022 than in the other years, which could be attributed to fewer nucleated particles growing into larger sizes when new particle formation (NPF) events occurred in 2022 (as shown in Fig. S2). In contrast, the OOM concentration was higher in 2022 than in 2021 (Fig. 4c), so reasons for the less efficient growth of the nucleation-mode particles in 2022 should be explored further. The accumulation-mode particle number concentration (N 100−1000 ) was significantly lower in 2021 than in other years, suggesting a suppressed growth of particles from the Aitken mode to larger sizes in measured air masses. This is consistent with the fact that the wind was mainly from the northwest in 2021, bringing clean air with low gaseous precursor concentrations. Thus, the overall low N 100−1000 in 2021 could also be associated with the more frequent clean air episodes, with PM 2.5 ≤ 35 µg cm −3 in 2021 (67 %, Fig. S3a).

Particles and ions
The 4-year variation of N 1.3−3 was contributed by several factors associated with both sources and sinks of these particles. In urban Beijing, as reported in previous studies Deng et al., 2021Deng et al., , 2020Du et al., 2022a), the formation and growth of new particles are influenced by the SA concentration, background aerosol population, amine concentration and ambient temperature, among other things. The SA concentration decreased from 2019 to 2022, and the PM 2.5 mass concentration, which is strongly associated with the coagulation caused by background aerosols, was the lowest in 2021 (Fig. 9). Also, as shown in Fig. 2, the ambient temperature was the lowest in 2021, causing favorable conditions for the formation of new particles because small clusters evaporate less efficiently when temperatures are lower. However, with such favorable conditions, the N 1.3−3 was still the lowest in 2021. The intensities of NPF events were lower in 2021 than 2022, while the NPF event frequency was higher in 2021. A possible reason for this observation is amines, the concentrations of which could have been reduced in 2021 due to limited traffic emissions (indicated by NO x ). However, the influence of this factor is uncertain without simultaneous measurements. The opposite 4-year change of BLH can also partly explain that of N 1.3−3 .
Unlike the number concentration of sub-3 nm particles, the ion concentrations in both 0.8-2 (N 0.8−2, ions ) and 2-4 nm (N 2−4, ions ) size ranges were higher in 2021 than in the other years (Fig. 7), possibly due to variations in ion source strengths. The ion concentration was higher in the size range of 0.8-2 nm than in the size range of 2-4 nm, which is consistent with previous studies . Figure S8 shows the levels of N 1.3−3 and N 0.8−2, ions under different pollution levels among those 4 years. We found that N 1.3−3 decreased with the increasing PM 2.5 levels, while sub-2 nm ions showed relatively constant concentrations under different pollution levels. Both sub-3 nm particles and sub-2 nm ions are easily to be scavenged by larger particles, and ions are more vulnerable considering their higher diffusivity. The relatively constant concentrations of sub-2 nm ions under different pollution levels might result from the ion source strengths increasing as the PM 2.5 concentration increased, which needs further exploration. Figure 8 shows that there was a positive and relatively strong association between the concentrations of 2-3 nm particles and the concentrations of SA1 and SA2 in every year, probably because the formation and growth of nucleated particles are strongly dependent on SA. A weak positive association was also observed between the 2-4 nm ion concentration and SA1 and SA2 concentration. Since 2-4 nm ions originate from collisions between small ions and particles, their concentration is expected to be correlated with SA2, especially as the concentrations of 2-3 nm particles were strongly correlated with concentrations of SA2.

PM 2.5 and its composition
The PM 2.5 mass concentration generally decreased from 2019 to 2021 while rebounding in 2022 to reach a similar level as in 2020 (Fig. 9). In 2019, although the median value was similar to the values in 2020 and 2022, the average PM 2.5 mass concentration was higher than in 2020 and 2022. This is because there were more severe haze events in 2019 (PM 2.5 > 150 µg cm −3 , 7 %, Fig. S3a), which increased the average PM 2.5 mass concentration considerably. The lowest  PM 2.5 mass concentration coincides with the highest BLH in 2021 ( Fig. 2e and f). In addition, winds were mainly from the northwest in 2021, bringing clean air masses with low PM loadings, and hence the frequency of PM 2.5 ≤ 35 µg cm −3 was the highest in 2021 (67 %, Fig. S3). These results indicate the important influence of meteorological conditions on air quality. PM 2.5 concentrations are also driven by primary emissions and secondary production, discussed in more detail below.
The mass concentrations of Chl, BC and Org showed similar 4-year variations to the PM 2.5 mass concentration, with a decrease from 2019 to 2021 and then an increase in 2022 (Fig. 10). Related mainly to primary sources, the mass concentrations of Chl and BC had similar 4-year variations to CO, a tracer for combustion sources. The mass concentration of Org had a similar 4-year variation to Chl and BC, which is because primary emissions like combustions partly contribute to the Org in the particle phase. Another important contributor to Org is secondary formation processes. Previous studies suggest that OOMs, which also decreased from 2019 to 2021 and then increased in 2022, can contribute a substantial fraction to secondary organic aerosol mass concentrations , Therefore, the 4-year variation of Org can be attributed to both primary emissions and secondary production processes.
However, unlike Chl, BC and Org, sulfate showed the lowest concentration in 2022. The generally decreasing trend of the sulfate mass concentration is consistent with the decreasing trend of the SO 2 concentration from 2019 to 2022, reflecting the significant reduction of SO 2 emissions that China has made during recent years, especially in winter.
Nitrate and ammonium had their lowest concentrations in 2021 with an increase in 2022, but, unlike Chl, BC and Org, they had higher concentrations in 2020 than in 2019. The similar 4-year changes of nitrate and ammonium are mainly because ammonium tends to be associated with nitrate and sulfate in the particle phase, and the contribution of nitrate becomes dominant due to the significant reduction of SO 2 emissions. With lower NO x in 2020 than 2019 ( Fig. 3b and c), the higher NO 3 concentrations in 2020 than those in 2019 suggest that more nitrogen oxides were transformed into particulate nitrate. Figure S12 shows that at almost all PM 2.5 levels, the mass concentration of nitrate in 2020 was higher than that in 2019, while NO 2 and NO concentrations were generally lower in 2020. One possible reason for this is the much higher RH in 2020, leading to more aerosol liquid water content, which promotes the uptake of HNO 3 and the hydrolysis of N 2 O 5 (Y. Wang et al., 2020a).
Despite some 4-year variations of the mass fractions of the different aerosol component, Org always contributed the largest fraction to the mass concentration of PM 2.5 in every year, followed by nitrate (Fig. 11a). The nitrate fraction generally increased from 20 % in 2019 to 28 % in 2022, suggesting a more and more important role of nitrate in PM pollution. The important role of the organics and nitrate in the PM 2.5 composition has also been found in previous studies in Beijing Sun et al., 2020;Hu et al., 2021). Figure 11b shows that when the PM 2.5 mass concentration increased from smaller than 35 µg m −3 to larger than 150 µg m −3 , the mass fraction of secondary inorganic aerosol (SIA), including nitrate, sulfate and ammonium, increased, while the mass fraction of organics decreased in every year. This phenomenon was mainly driven by the increase of the nitrate fraction. Especially in 2021 and 2022 when the PM 2.5 mass concentration was between 35 and 150 µg m −3 , the mass fraction of nitrate was the largest among the different components. Our results indicate the important role of nitrate in haze formation and the urgent need for significant reductions in NO x emissions in order to improve the air quality in Beijing.
A number of studies have investigated the yearly changes of chemical compositions of particulate matter in Beijing in recent years Xu et al., 2019;Zhou et al., 2019;Zhang et al., 2020;Hu et al., 2021). Due to the implementation of air quality regulations, most of the compounds in particulate matter have decreased significantly. For example, chloride of PM 1 decreased by 65 %-89 % from 2011-2012 to 2017-2018 in urban Beijing, measured using an ACSM . In this study, the organics concentration also decreased by 37 %-70 %. However, the nitrate concentration did not decrease but even increased. The nitrate concentration was also found to increase by ∼ 4 % during wintertime from 2007 to 2017 in a previous study . Additionally, many studies found that nitrate gradually becomes the major contributor to the chemical composition of particulate matter, overweighing organics, especially during severe haze periods Hu et al., 2021). These results are consistent with ours, indicating the importance of reducing NO x concentration in Beijing.

Meteorological conditions
There are three notable meteorological characteristics for the special event periods. First, the wind direction distribution during reference, CNY and COVID periods were generally similar and followed the typical diurnal wind pattern in Beijing (Fig. 12a-c). Specifically, from midnight to noon, winds mostly blew from the northwest and sometimes also from the northeast, while from the afternoon to evening, the wind direction turned to the south or southeast. During the Olympics, however, the wind pattern was not as strong as in the other events (Fig. 12d). From midnight to noon, the dominant wind direction was the southwest, instead of northwest as during the other events, while from the afternoon to midnight, there was no prevailing wind direction. This feature could have affected the air quality positively during the Olympics as there were less frequent southerly winds which are typically associated with more polluted air masses. Second, the temperatures during the Olympics (median −0.6 • C) and CNY (median −0.5 • C) were lower than during reference (median 2.5 • C) and COVID (median 2.6 • C) periods (Fig. 13a). The temperature variation may cause changes in the rate of gas-and particle-phase reactions as well as in the volatility of oxygenated organic molecules and thereby may have a slight influence on the air quality. Third, the RH during the COVID period (median 52 %) was higher than during the other three periods (median 24 %-27 %) (Fig. 13b).
And this high RH may change the phase state of particles and increase secondary aerosol formation through heterogeneous reactions. Meanwhile, the wind speeds (Figs. 12 and 13d) and BLH levels ( Fig. 13e and f) during the four periods were comparable to each other, so the vertical mixing of air masses was likely relatively similar among different periods.

Trace gases
Since the restrictions during the COVID, Olympics and CNY periods focus on different sections of social life and production activities, trace gases may have different responses to different periods. During the Olympics, the southerly wind was not that frequent, so this period was the least influenced by polluted air masses. Consequently, CO, NO, NO 2 and SO 2 were the lowest, and O 3 was the highest during the Olympics compared with other three periods. During COVID, however, the southerly wind appeared with the highest frequency, leading to the most severe pollution transportation. Thus, the CO concentration was the highest (849 ppb). It can also be found that, despite the most polluted air masses during COVID, the concentrations of NO, NO 2 and SO 2 were lower than the reference period, suggesting that the restrictions on NO, NO 2 and SO 2 during that period were quite efficient. In addition, compared with NO 2 , the restriction on NO and SO 2 was much more obvious. NO 2 has large sources from secondary formation and is able to transport over long distances (Zhu et al., 2021;Tan et al., 2022); hence, the transportation of NO 2 likely partly compensates for the effect of the NO 2 emission control. Under the synergistic influences of pollution transport and COVID control measures, the O 3 concentration was only slightly higher during COVID than during the reference period. As mentioned before, during CNY, the wind direction had a pattern similar to that during the reference period, and the wind speed as well as BLH was only slightly higher than during the reference period. Therefore, meteorological conditions were similar during these two periods. Consequently, CO was only a bit lower during CNY than during the reference period. For NO and NO 2 , however, their concentrations were much lower during CNY than during the reference period. This is likely due to the reduction of NO x emissions from vehicles, as a large fraction of people went to their hometowns from Beijing. Lower NO x levels also led to higher O 3 levels. The SO 2 concentration was only slightly lower during CNY than during the reference period. During the winter heating period in urban Beijing, SO 2 should mainly come from the combustion of fossil fuel Meng et al., 2016), so the reduced traffic emissions had a minor influence on SO 2 concentrations. In addition, from New Year's Eve to the next following days, SO 2 can also be emitted from fireworks (Foreback et al., 2022). We then investigated these gases under different PM levels. When PM 2.5 exceeded 35 µg m −3 , SO 2 during COVID was much lower than during other three periods (Fig. S10), and this feature was associated with higher RH (Fig. S9) and higher particulate sulfate concentrations (Fig. S13). Therefore, the conversion of gaseous SO 2 to particulate sulfate was likely enhanced during COVID when PM 2.5 was larger than 35 µg m −3 , leading to the low concentration of gaseous SO 2 . Additionally, O 3 and PM 2.5 showed a non-linear relationship; O 3 concentration first decreased and then increased with increasing PM 2.5 , for the period from 4 to 20 February (Fig. S10), which is consistent with previous studies (Y. Wang et al., 2020b;Zhao et al., 2020). As the level of O 3 is affected by many factors, such as meteorological conditions, regional transport (Lin et al., 2019;Zhao et al., 2021;Liu et al., 2019;Ge et al., 2012), and photochemical processes associated with solar radiation, NO x and VOCs (Tan et al., 2018;, further analysis is highly needed for delving into the reasons behind such behavior of O 3 .
In summary, in urban Beijing, concentrations of gaseous pollutants are controlled by their emission sources and the meteorological conditions. Thus, during the Olympics, favored by the clean air mass, CO, NO, NO 2 and SO 2 were dramatically reduced. During COVID, under the influence of extremely polluted air masses, the strict restrictions on people's movement which affected traffic emissions and probably other production activities, CO concentration increased a lot, while NO and SO 2 concentrations reduced, and NO 2 reduced slightly. During CNY, under the influence of unfavorable meteorological conditions, reduced traffic emission and probably some other reduced production activities, CO and SO 2 concentrations did not decrease, while NO x reduced.

Condensable vapors
The responses of SA and OOMs to the three restriction periods were quite different. For SA, the concentrations of SA1 Olympics and CNY periods. The black value inside each box is the median value of corresponding parameter. The gray values are the normalized median values to the lowest median value in each subplot, which is used for quantitative variation comparison between different parameters. Please note that for UVB, only the daytime (08:00-16:00) dataset was used. and SA2 were the lowest during COVID (median value are 4.2×10 5 and 5.7×10 3 cm −3 for SA1 and SA2, respectively) and were comparable in the other three periods (median values are 6.0-8.7 × 10 5 and 1.5-1.9 × 10 4 cm −3 for SA1 and SA2, respectively) ( Fig. 15a and b). Since SA is mainly produced from SO 2 oxidation by OH radicals (Finlayson-Pitts and Pitts, 2000) and is primarily scavenged by condensation onto particles Guo et al., 2021;Yang et al., 2021), the lowest concentration during COVID could be attributed to the highest condensation sink (Fig. S4b) accompanied by the low SO 2 concentration (Fig. 14e). Despite the lowest SO 2 concentrations during the Olympics, SA1 was not low (Fig. 15a), and SA2 was even the highest (Fig. 15b). This was probably caused by the synergistic effects of decreased condensation sink (Fig. S4b) and increased atmospheric oxidation capacity during the Olympics (which could be partially indicated by the highest O 3 , Fig. 14d). Meanwhile, compared with COVID and reference periods, the ratios of SA2 / SA1 during the Olympics and CNY were higher, suggesting that the cluster formation efficiency of sulfuric acid was enhanced during those two periods.
The concentration of the total OOMs (median values are 1.5-2.8 × 10 7 cm −3 , Fig. 15d) was the highest during COVID. Compared with SA, the production and loss processes of OOMs are much more complex: they can be generated from the oxidation of VOCs (volatile organic compounds) and OVOCs (oxygenated volatile organic compounds) or be evaporated from the particle phase (Kohli and Davies, 2021;Yli-Juuti et al., 2017;Wilson et al., 2015). When it comes to the loss processes, OOMs can be consumed by further oxidation processes or through absorption onto particle surfaces via pure condensation or reactive uptake. Among all these factors, precursor VOCs and OVOCs could have great impacts. Since such precursors tend to come along with pollution Niu et al., 2022), OOM concentration also showed a strong association with the pollution level (Fig. S11). This is probably why the OOM concentration was the highest during COVID when severe haze (PM 2.5 > 75 µg m −3 ) was most frequent and the lowest during the Olympics when there was almost no severe haze (Fig. 3b). For OOMs themselves, Figs. 15 and S11 together show that under the same PM 2.5 level, both the averaged oxygen and nitrogen numbers during the Olympics were the highest, which implies that not only the oxidation state of OOMs was enhanced, but also the involvement of NO x was more effective during the Olympics.

Particles and ions
As shown in Fig. 16, N 1.3−3 was almost at the same level during the reference, COVID, Olympics and CNY periods. N 3−25 was lower during CNY period, and N 25−100 was slightly lower during the Olympics than during the other periods. However, N 100−1000 was significantly lower during the Olympics compared with the other periods.
The small variability of N 1.3−3 between the different periods indicates that the special events may have only little effect on the formation of sub-3 nm particles, i.e., new particle formation processes, consistent with a previous study . However, for the accumulation particles, the significantly lower number concentrations during the Olympics indicates a great impact of the reduction of anthropogenic emissions on large particles during this period. Also, the frequent wind direction from the north indicates the transport of clean air mass during the Olympics.
Ion concentrations in the different size ranges were significantly higher during the Olympics compared with the other periods (Fig. 17). During COVID, the ion concentrations were the lowest. Such differences between the different periods basically follow the 4-year changes of ion concentrations. As discussed in Sect. 3.1.5, the ion concentrations were higher in 2022 compared with the other years, which is consistent with the higher ion concentrations during the Olympics in 2022 than during the COVID, CNY and reference periods.

PM 2.5 and its composition
The median mass concentration of PM 2.5 was the lowest during the Olympics and the highest during COVID, as compared with the other periods (Fig. 18). During the CNY period, the median mass concentration of PM 2.5 was the second highest. The differences in the PM 2.5 mass concentration between the different periods were not due to different values of the BLH, since the BLH was similar between the COVID and Olympics periods and the highest during the CNY period (Fig. 13e).
Favorable meteorological conditions for transport of clean air mass during the Olympics are one of the main reasons leading to the lowest PM 2.5 mass concentration. In contrast, unfavorable meteorological conditions during the COVID and CNY periods worsened air quality during these two periods. As shown in Fig. 12, the wind cycle was obviously different during the Olympics than during the other periods. In urban Beijing, surrounded by mountains to its west, north and northeast, there is a typical wind cycle, with winds tending to come from the north or northwest in the morning and from the south or southeast in the evening. Such a cycle was obvious during the CNY and COVID periods, allowing the transport of polluted air masses to urban Beijing from the south or southeast, thereby leading to severe haze events.
Another reason for the much lower mass concentrations of PM 2.5 during the Olympics compared with the reference period was the restriction measures on reducing anthropogenic emissions during this time period. Consistently, there were much lower concentrations of primary particles, i.e., BC and Chl, during the Olympics than during the reference period (Fig. 19). Although strict restrictions were also performed during COVID, BC and Chl concentrations were slightly higher than during the reference period. This should not be interpreted as emission control not being effective during COVID. In fact, BC during COVID was lower than during the reference at the same PM 2.5 levels (Fig. 13f). The higher average BC concentration during COVID was associated with the higher frequency of haze episodes during COVID. During CNY, the relatively high primary aerosol concentrations seem to be associated with fireworks.
Secondary aerosol formation is important in understanding PM 2.5 levels as this process is often the main contributor to the particle mass concentration. The much higher fraction of SIA compared to the other components during COVID (Fig. 20a) indicates the enhancement of atmospheric oxidation capacity, which might explain the more severe haze events during this period compared with the other periods. As shown in Fig. 20b, the mass fraction of SIA was much higher during COVID, and when the PM 2.5 mass concentration was higher than 35 µg m −3 , the mass fraction of SIA reached almost 70 %. The much higher values of RH during COVID could have promoted heterogeneous reactions to produce higher SIA concentrations (Fig. 13b). In addition, nitrate showed higher concentrations and mass fractions when PM 2.5 increased from 35 µg m −3 to larger concentrations during COVID than other periods, although there was a reduction in NO x emissions due to strict restrictions on vehicles this period. This phenomenon indicates that secondary formation of nitrate was enhanced and thus contributed to the formation of haze during COVID. In the case of sulfate, although the SO 2 concentration was low during COVID, the overall mass concentration and mass fraction of sulfate were higher than during the other periods, which might be associated with the much higher values of RH during COVID. With higher mass concentrations and fractions of nitrate and sulfate, the ammonium was correspondingly higher during COVID than during the other periods.

Summary and conclusions
This study investigated the air quality in the megacity of Beijing during wintertime, with an emphasis on the 4year changes and the influences of special events (COVID, Olympics and CNY). Numerous variables, including meteorological parameters, concentrations of trace gases, gaseous sulfuric acid and oxygenated organic molecules, and number concentrations of atmospheric particles and ions, as well as PM 2.5 and its composition, were systematically analyzed. The comprehensive datasets span the dates from 1 January to 20 February in the years from 2019 to 2022.
In the first part, the 4-year changes of atmospheric parameters were explored, using the datasets from 1 to 22 January. Generally, the meteorological conditions for 2019 and 2020 were similar, where temperature, wind speed and BLH are comparable, and the wind mainly blew from the northwest from night to noon and turned to the south from the afternoon. In 2022, the wind mostly blew from the southeast from midnight to noon and turned to the northeast from the afternoon. The wind speed in that year was a little higher than 2019 and 2020, and the temperature and BLH were slightly lower than 2019 and 2020. In 2021, the wind was mostly from the west or northwest. Thus, the air masses in this year were much cleaner than other 3 years due to the least influence of polluted south air masses. This clean condition in 2021 was also associated with the highest wind speed and the maximum BLH, leading to the strongest atmospheric diffusion capacity. Besides, RH showed an increasing pattern from 2019 to 2022, which may facilitate the formation of secondary aerosols through heterogeneous reactions (Sun et al., 2018).  For conventional gaseous and particulate pollutants, the 4-year changes of CO, NO 2 , total PM 2.5 , organic aerosol, chloride and black carbon were generally opposite to those of BLH which decreased from 2019 to 2021 and increased in 2022. The atmospheric conditions of 2021 were also the cleanest. Thus, the atmospheric diffusion capacity and the overall cleanness of the air masses may control the variations of the above pollutants. Unlike them, SO 2 monotonically decreased from 2019 (3.06 ppb) to 2022 (0.68 ppb) regardless of the changes in meteorological conditions, suggesting that SO 2 emissions and sources decreased year by year. Due to the decline of SO 2 , particulate sulfate also decreased from 2019 (2.2 µg m −3 ) to 2022 (1.1 µg m −3 ). NO in 2022 was around the same as 2021 despite more polluted conditions, indicating that the traffic emission in 2022 might have reduced. For O 3 , it showed an opposite variation with NO 2 as it increased from 2019 (6.3 ppb) to 2021 (15.4 ppb) and dropped in 2022 (8.1 ppb). For nitrate, its contribution to PM 2.5 increased from 20 % in 2019 to 28 % in 2022. Moreover, its contribution to PM 2.5 was more pronounced during polluted conditions (PM > 35 µg m −3 ), indicating the important role of nitrate in haze formation and the urgent need for significant reductions in NO x emissions.
For unconventional parameters, owing to the reduction of SO 2 , the concentrations of gaseous sulfuric acid monomer and dimer generally decreased from 2019 (9.6 × 10 5 and 2.8 × 10 4 cm −3 for sulfuric acid monomer and dimer respectively) to 2022 (4.5×10 5 and 3.1×10 3 cm −3 for sulfuric acid monomer and dimer respectively). For total OOMs, however, its concentration decreased from 2019 to 2021 and increased in 2022. This variation was generally opposite to BLH but with larger extent, suggesting that OOM concentration was controlled both by meteorological conditions and local chemical production. Meanwhile, both oxygen and nitrogen contents of gas-phase OOMs increased year by year, implying that not only the oxidation state of those compounds was increased, but also NO x was involved more efficiently in their formation processes. For particles, it was found that with higher SA concentrations and NPF frequencies in 2021 than in 2022, and with the lowest concentrations of background aerosols and the lowest ambient temperatures in 2021, N 1.3−3 was still the lowest in 2021. This is out of the scope of conventional understanding, and further studies are needed to explore the reasons behind this phenomenon.
In the second part, the influences of special events, including COVID lockdown, the Beijing Winter Olympics and Chinese New Year, were investigated, using the datasets from 4  to 20 February. The wind direction distribution during the reference, COVID and CNY periods were similar in that the wind mostly blew from the northwest or northeast from midnight to noon and turned to the south or southeast in the afternoon. During the Olympics, however, the dominant wind was from the southwest from midnight to noon, while there was no prevailing wind from the afternoon. Thus, the air masses during the Olympics were the cleanest due to less frequent southerly wind. In addition, there was not too much difference in wind speed and BLH for four periods. One thing that needs to be mentioned is that the RH during COVID was much higher than the other three periods, which should promote the growth of secondary aerosols.
During the Olympics, favored by the clean conditions, most trace gases (CO, NO, NO 2 and SO 2 ) and particulate matters (total PM 2.5 , organic aerosol, nitrate, ammonium, chloride and black carbon) were dramatically reduced compared with the reference period. During COVID, influenced by the extremely polluted conditions, CO and all particulate matter (total PM 2.5 , organic aerosol, sulfate, nitrate, ammonium, chloride and black carbon) increased a lot. Despite the severe haze, NO and SO 2 reduced, and NO 2 reduced slightly, suggesting that the restrictions on traffic emission and probably on other production activities were effective in reducing NO x and SO 2 . Further, the high RH during COVID likely facilitated the gas-to-particle conversion of sulfur-and nitrogen-containing compounds, resulting in particularly high concentrations of sulfate (8.5 µg m −3 ), nitrate (16.2 µg m −3 ) and ammonium (8.6 µg m −3 ). During CNY, affected by unfavorable meteorological conditions, reduced traffic emission and probably other reduced production activities, CO, SO 2 , total PM 2.5 , nitrate and black carbon were comparable with the reference period, while NO x reduced; organic aerosol, ammonium and chloride increased; and sulfate increased to a large extent. During that period, O 3 had a non-linear relationship with PM 2.5 that its concentration first decreased and then increased with increasing PM 2.5 , which is different from the period of 1-22 January. This suggests that the synergetic control strategy of PM 2.5 and O 3 should consider the effect of seasonality.
SA and OOMs are secondary products, so their variations were complex and were not directly linked to the changes of meteorology or restriction. The level of SA1 was mainly influenced by CS, SO 2 and UVB (Petäjä et al., 2009;Yang et al., 2021). Compared with the reference period, SA1 was comparable during the Olympics, lower during CNY and much lower during COVID. The concentration of SA2 was mainly controlled by SA1 and CS. Compared with the reference period, SA2 during the Olympics and CNY was comparable, while it was much lower during COVID. As the COVID period was the most polluted, the above results suggest that the formation of SA was suppressed under polluted condition. For total OOMs, however, it had the lowest concentration in the Olympics while the highest concentration in COVID, indicating that the level of total OOMs was higher under polluted conditions. Meanwhile, the oxygen content was the highest during the Olympics while the lowest during COVID, which implies the oxidation state of OOMs likely decreased with increasing pollution level. For atmospheric particles, N 1.3−3 was almost the same during reference, COVID, Olympics and CNY periods, indicating that the special events may only have little impact on new particle formation processes. However, N 100−1000 was significantly lower during the Olympics than other periods, suggesting the significant reduction of anthropogenic emissions during this period.
These results provide useful information on how the air quality in urban Beijing has changed in recent for years and how the air quality was affected by meteorological conditions as well as different emission reduction scenarios. All of these will help in planning more targeted and sustainable long-term pollution control plans. Data availability. Datasets for this paper can be accessed at https://doi.org/10.5281/zenodo.7100748 .
Author contributions. YG, CD, TK, WD and MK designed the research; YG, CD, FZ, YiL, JW, CH, JZ, ZW, JJ, YiZ, TL, YuZ, BS, WM, YoL, CY, and MK collected the data; YG, CD and AO analyzed data with the help from TK, WD, TN, VMK, MX, CY and MK; and YG, CD, TK and WD wrote the paper with inputs from all co-authors.
Competing interests. At least one of the (co-)authors is a member of the editorial board of Atmospheric Chemistry and Physics. The peer-review process was guided by an independent editor, and the authors also have no other competing interests to declare.

Disclaimer.
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Special issue statement.
This article is part of the special issue "Pan-Eurasian Experiment (PEEX) -Part II". It is not associated with a conference.