key: cord-0713553-0qgxqlsk
authors: Bai, Kaixu; Wu, Can; Li, Jianjun; Li, Ke; Guo, Jianping; Wang, Gehui
title: Characteristics of Chemical Speciation in PM(1) in Six Representative Regions in China
date: 2021-04-07
journal: Adv Atmos Sci
DOI: 10.1007/s00376-020-0224-2
sha: 71ceef71c85c97efa715bb2869a49ccf2553be9b
doc_id: 713553
cord_uid: 0qgxqlsk

A better knowledge of aerosol properties is of great significance for elucidating the complex mechanisms behind frequently occurring haze pollution events. In this study, we examine the temporal and spatial variations in both PM(1) and its major chemical constituents using three-year field measurements that were collected in six representative regions in China between 2012 and 2014. Our results show that both PM(1) and its chemical compositions varied significantly in space and time, with high PM(1) loadings mainly observed in the winter. By comparing chemical constituents between clean and polluted episodes, we find that the elevated PM(1) mass concentration during pollution events should be largely attributable to significant increases in organic matter (OM) and inorganic aerosols like sulfate, nitrate, and ammonium (SNA), indicative of the critical role of primary emissions and secondary aerosols in elevating PM(1) pollution levels. The ratios of PM(1)/PM(2.5) are found to be generally high in Shanghai and Guangzhou, while relatively low ratios are seen in Xi’an and Chengdu, indicating anthropogenic emissions were more likely to accumulate in forms of finer particles. With respect to the relative importance of chemical components and meteorological factors quantified via statistical modeling practices, we find that primary emissions and secondary aerosols were the two leading factors contributing to PM(1) variations, though meteorological factors also played important roles in regulating the dispersion of atmospheric PM.

Aerosols are a critical component of the Earth's atmosphere and play important roles in deteriorating air quality and changing climate (Xin et al., 2015; Guo et al., 2017; Li et al., 2017 Li et al., , 2019 Lou et al., 2019; Zhao et al., 2020) . Due to the light extinction effect, aerosols in the atmosphere can significantly degrade the visibility (Charlson, 1969; Appel et al., 1985; Shi et al., 2014) . By scattering and absorbing solar and terrestrial radiation, atmospheric aerosols can significantly alter the Earth's radiation budget, which in turn changes the climate through complex cloud-aerosol-precipitation interactions (Guo et al., 2017; Li et al., 2019; Zhao et al., 2020) . Moreover, atmospheric aerosols have negative impacts on human health by causing respiratory and cardiovascular diseases, especially those in the forms of fine particulate matter (Ebenstein et al., 2017; Burnett et al., 2018; Zou et al., 2019; Yue et al., 2020) .

Aerosols in the ambient environment commonly exist in the form of atmospheric particulate matter (PM), which is a compound of multiple minerals and materials originating from nature and anthropogenic emissions. Given different compounds, the climatic and environmental effects of aerosols are thus determined by both physical properties (e.g., particle size) and chemical constituents of PM (Pósfai and Buseck, 2010; Yang et al., 2011; Xin et al., 2015; Zhang et al., 2018) . It is known that mass concentrations of PM and chemical compositions are determined by many factors including primary emissions, chemical reactions, depositions, as well as meteorological conditions that control both local dispersion and long-range transport of PM Wang et al., 2014 Wang et al., , 2019 . Measuring chemical components in PM is thus crucial for understanding the formation of atmospheric aerosols, especially the complex mechanisms behind frequently occurred haze pollution events (e.g., Wang et al., 2016 Wang et al., , 2017 Wang et al., , 2020 . Observational networks have been extensively established to monitor the physical, chemical, and optical properties of atmospheric aerosols worldwide, such as the Aerosol Robotic Network (AERONET) (Holben et al., 1998) and Global Atmosphere Watch: Aerosols (GAW) (WMO, 2001) . Such networks provide ample observational data to help improve understanding of aerosol properties around the world. Similar networks were also established in China. In 2011, the Chinese Academy of Sciences (CAS) launched the Campaign on Atmospheric Aerosol Research network of China (CARE-China) . This is also the first comprehensive research platform for observing atmospheric aerosols in the country, aiming to provide in situ measurements of anthropogenic aerosol emissions to better understand the climatic and environmental effects of aerosols in China (Xin et al., 2015) .

By taking advantage of such invaluable in situ measurements, complex atmospheric haze pollution mechanisms in China have been better elucidated (Guo et al., 2014; Huang et al., 2014; Wang et al., 2020) . For instance, based on chemical data observed in representative mega cities across China, Yang et al. (2011) found that PM 2.5 mass concentrations and chemical compositions varied substantially over space and that the observed severe PM 2.5 pollution events in eastern China were largely attributable to local formation/production and regional transport of the secondary aerosols. Utilizing field measurements of aerosol properties at an urban site in Beijing during January and February 2015, Zamora et al. (2019) found that frequently occurring wintertime haze events in Beijing were attributable to the efficient nucleation and secondary aerosol growth under high gaseous precursor concentrations and stagnant air conditions. Likewise, by taking advantage of comprehensive measurements and modeling simulations, Huang et al. (2020) found that the haze events during the COVID-19 lockdown in eastern China were driven by the formation of substantial secondary aerosols due to increases in atmospheric oxidizing capacity.

Given the adverse effect on public health, PM 2.5 mass concentration and its speciation have been extensively examined in the literature Zamora et al., 2019) . In the present study, emphasis is given to even finer aerosols with an aerodynamic size less than 1.0 μm, which is also known as submicron PM (abbreviated as PM 1 hereafter). Given this smaller size, PM 1 has higher extinction efficiency and larger negative health impacts (Shi et al., 2014) . Therefore, evaluating the variability of PM 1 and its chemical speciation in space and time is of great significance for understanding of its important roles in haze formation and visibility degradation. By using three-year-long aerosol measurements that were sampled in six representative cities in China, here we aim to first examine variations in PM 1 and its major chemical compositions in space and time, and then to explore their potential roles in the formation of severe fine particle pollution events. The following science questions are to be answered: (1) What chemical constituents were dominant factors in elevating PM 1 pollution levels in the selected regions? (2) Was there a common mechanism for the formation of PM 1 pollution events across these regions? and (3) Did meteorological conditions play important roles in modulating PM 1 pollution levels?

In this study, aerosol measurements collected from the CARE-China network at six distinct regions during the period of 2012 to 2014 were used. These six sampling sites were deployed in different geographic regions representative of diverse emission types given different population size and urbanization levels. Table 1 gives a brief description of geographic location and the ambient background (rural or urban site) of these six sampling sites. Among them, the site over Qinghai Lake (abbreviated as QHL hereafter) is referred to as a background site since it was deployed in the Tibetan Plateau where it was less influenced by human activity. In contrast, the other five sites were all deployed in urban regions, thereby representing a pollution regime that was highly influenced by anthropogenic emissions (Ren et al., 2018) .

Size-segmented atmospheric PM was collected using nine-stage Anderson samplers (Anderson Series 20−800, United States) at an airflow rate of 28.3 L min -1 with cutoff sizes of 9.0, 5.8, 4.7, 3.3, 2.1, 1.1, 0.7 and 0.4 μm. Each set of aerosol samples were measured twice monthly, with every data collection task lasting for three days at five urban sites, while six days were used at the QHL site to guarantee an adequate number of valid samples. More detailed description of sampling collection and storage can be found in the literature like Ren et al. (2018) and Xin et al. (2015) .

A DRI 2001 thermal/optical carbon analyzer was used to measure organic carbon (OC) and elemental carbon (EC) compositions in size-resolved particles at each sampling site by following the Interagency Monitoring of Protected

Visual Environments (IMPROVE) thermal/optical reflectance (TOR) protocol (Chow et al., 2004 (Chow et al., , 2007 . Meanwhile, soluble inorganic salts including five cations (Na + , NH 4 + , K + , Mg 2+ , and Ca 2+ ) and three anions (Cl -, , and ) were analyzed by ion chromatography (IC), which was equipped with a separation column and a suppressor (Xin et al., 2015) . A detailed descriptions of instrumentation and analytic procedures to measure these compositions can be found in Xin et al. (2015) .

In this study, PM with aerodynamic diameter no more than 1.1 μm were aggregated to represent PM 1 . Mass concentrations of PM 1 and five major chemical compositions of OC, EC, sulfate ( ), nitrate ( ), and ammonium ( ), were extracted and analyzed. To account for unmeasured atoms, organic mass (OM) was derived by multiplying a commonly suggested scale factor of 1.4 to OC concentration, though this scale is subject to many factors and may vary with seasons (Turpin and Lim, 2001; Bae et al., 2006; Yang et al., 2011) .

Meteorological conditions are believed to play important roles in modulating ambient PM concentrations other than primary emissions since they can modulate the dispersion and accumulation of PM (Zheng et al., 2017; An et al., 2019) . Four primary meteorological factors that are highly associated with aerosol variations-air temperature (T), wind speed (WS), relative humidity (RH), and boundary layer height (BLH)-were hereby used. All these data were retrieved from the Meteorological Information Comprehensive Analysis and Process System (MICAPS), except BLH, which was retrieved from the ERA-Interim reanalysis that was provided by the European Centre for Medium-Range Weather Forecasts (ECMWF). 3-h MICAPS data at the sampling sites were averaged for every 24-hour to approximate daily averages.

To better understand the haze formation mechanisms over the selected regions, we attempted to resolve statistical associations between PM 1 and chemical species as well as meteorological factors at each individual site. The ultimate goal is to evaluate the relative importance of each factor contributing to PM 1 loading levels to determine the leading influential factors in each city. The random forest (RF) method was hereby used to statistically link PM 1 mass concentration to its five major chemical species and four meteorological factors. The reason to use the random forest method is twofold. First, RF has a good generalization performance. Second, in contrast to many "black-box" models, RF has a unique variable importance estimation capacity (Ho, 1998; Altmann, et al., 2010) .

The relative importance is calculated by estimating the response of out-of-bag error to the random permutation of out-of-bag data across each variable. Therefore, the larger the increase in error budget, the more important the variable. Such a metric is thus oftentimes used to assess the prediction power of the given predictor, enabling us to determine whether the given variable should be included in the modeling process or not. This capacity helps us better understand the intrinsic associations between the response variable and various explanatory factors (Bai et al., 2019a, b) . For the configuration of RF models, 500 regression trees were used in each random forest model, while each tree was grown on a bootstrap sample with a random subset of predictors selected at each split. In this study, our focus is not on the RF modeling performance given limited samples (as large as of 72 data pairs given two monthly sampling frequency). Rather, we aimed to evaluate the relative importance of each explanatory variable, i.e., the relative contribution of chemical species and meteorological factors in determining PM 1 mass concentration.

3.1. Major chemical species in PM 1 Figure 1 presents the temporal evolutions of PM 1 and five major chemical constituents that were sampled at each individual site during the field campaign between 2012 and 2014. As shown, concentrations of PM 1 and its chemical constituents varied significantly in space and time, with the max- Table 1 . Geolocation and ambient background information of six aerosol sampling sites used in this study. The abbreviation of each city/region was shown in brackets.

Longitude ( imum mean PM 1 concentration being triple that of the minimum (29.5−98 μg m −3 ). Noteworthy is that there were two large outliers (>120 μg m −3 ) in the PM 1 time series of QHL, and this could be due to possible measurement errors. According to the magnitude of mean PM 1 concentration, these six sampling regions could be generally grouped into three primary pollution regimes as clean (QHL: 29.5±25.1 μg m −3 ), polluted (SH and GZ: 54.2±23.3 and 57.3±20.9 μg m −3 ), and severely polluted (XA, CD, and UMQ: 88.1±34.7, 89.7±32.1, and 98±69.5 μg m −3 ), respectively. Temporally, PM 1 varied even more drastically, as high PM 1 loading (>100 μg m −3 ) appeared mainly in the winter, whereas low values were seen in the summer. Meanwhile, high PM 1 concentrations were oftentimes accompanied by a substantial increase in chemical constituents such as organic aerosols, sulfates and nitrates, revealing the significance of secondary aerosols in elevating ambient PM 1 loading. In order to better reveal the spatial and temporal variability of chemical constituents in PM 1 in these six regions, we compared the seasonal mean relative fraction of each chemical species in PM 1 . The results are presented in Fig. 2 . It is seen that total carbon (OC plus EC) and SNA (the sum of sulfate, nitrate, and ammonium) only accounted for a small portion (<50% in most regions) of PM 1 mass concentration, especially at the QHL site (21.4%). Given low primary emis-sions around the QHL site due to limited human activity, the low fraction of five major chemical constituents thus implies higher contribution of other components like mineral elements from manmade or natural primary sources to the observed PM 1 mass concentration. However, note- worthy is that the mean PM 1 concentration at the QHL site is even larger than some urban sites around the world, e.g., Vienna (14.9±7.7 μg m −3 ) (Gomišček et al., 2004) . Due to the data availability issue, the element analysis was not performed here to better examine the constituents of PM 1 at the QHL site; hence, reasons for this effect remain unclear.

An inter-comparison between five chemical constituents revealed that OM had the most significant contribution to PM 1 mass concentration since OM alone accounted for a mass fraction as large as of 10%−42% of PM 1 throughout the year at these six sampling sites. Such a high fraction of OM was also found in other cities like Beijing and Tianjin (Shao et al., 2018; Khan et al., 2021) . Temporally, this high fraction was even more impressive in the autumn and winter, particularly in cities like Xi'an and Guangzhou where OM increased more drastically during the wintertime than in the summer. This could be attributable to extensive primary emissions from biomass burning and coal combustion in the winter (Li et al., 2014; Wu et al., 2020) . In Guangzhou, the increase in OM could be ascribed to both primary and secondary organic aerosols from extensive traffic emissions (e.g., vehicles and ships), since the primary organic pollutants can be easily oxidized to form SOA through gas-phase oxidation chemistry in a wet ambient environment Qin et al., 2017) .

The SNA only accounted for 7%−23% (a mean value of 17.8%±5.9%) of PM 1 concentration. Notably, sulfate constituted much higher fractions (58%±6.6% of the SNA) than the other two counterparts (25.6%±4.3% for nitrate and 16.4%±3% for ammonium). In contrast, the EC constituted a much lower fraction with a percentage ranging from only 0.9%−6.5%. Spatially, Guangzhou had the largest EC fraction (6.5%) while the smallest fraction (<1%) was observed at the QHL site. Since the field measurements at the QHL were considered as the ambient background values, these inter-comparation results collectively imply that the higher PM 1 mass concentration in the other five major cities could be attributable largely to extensive anthropogenic emissions from primary sources like vehicles and heating related biomass and coal combustion. Figure 3 compares the seasonal mean mass ratio of OC

to EC (OC/EC) and to ( ). Noteworthy is that QHL had an extremely large mean OC/EC ratio (22.5) among these six regions, followed by Chengdu (7.4), Ürümqi (6.3), Shanghai (6.0), Xi'an (5.8), and Guangzhou (4.2). Meanwhile, the OC/EC ratio in Guangzhou was found to vary little across seasons, appearing to imply a relatively important role of primary emissions. This could be linked to considerable diesel and low quality bunker fuel powered ships in the Pearl River delta . On the contrary, the extremely large mean OC/EC ratio at the QHL site was mainly observed during summer and fall periods. Our further investigations showed that OC at the QHL had larger mean concentration (2.6 μg m −3 ) during summer and fall than in winter and spring (1.9 μg m −3 ). Given near constant low EC values during the sampling period, we speculate that the observed large OC/EC ratio at the QHL site might be associated with increased OC like biogenic sources (volatile organic compounds) due to higher temperature. Among the five major cities, the largest OC/EC ratio was observed by all in the winter, indicating more significant primary emissions during the wintertime.

The mass ratio of / was also calculated and compared across these six regions. This ratio has been frequently used as an indicator of the relative contribution of mobile sources against stationary sources, with a high mass ratio of / tending to imply a more important role of mobile sources over stationary sources and vice versa (Arimoto et al., 1996; Yao et al., 2002; Yang et al., 2011; Shao et al., 2018) . Among these six regions, all annual mean / ratios were lower than 1.0, with the lowest ratio observed at the QHL (0.36), followed by Guangzhou (0.37), Ürümqi (0.41), Shanghai (0.61), Chengdu (0.62), and Xi'an (0.71). Noteworthy is that the / ratio in Guangzhou was the second lowest, which was almost comparable to that of QHL. Given the consensus of larger amounts of vehicles and ships in Guangzhou and Shanghai, the relatively low / ratios in these two cities are thus insufficient to conclude a less important role for mobile sources, since many other contributing factors and even the degree of representativeness of the sampled data were not well accounted for. For instance, the largest Fig. 3 . Same as in Fig. 2 but for seasonal mean mass ratios of (a) OC to EC and (b) to .

ratios in Xi'an might be more likely associated with extensive coal combustion there (Zhang, 2014) , especially in the winter, and thus it could also yield large amounts of wintertime nitrate as well (Wang et al., 2018) . Overall, more evidence such as emission inventory should be applied in the future to better examine the relative contribution of mobile and stationary sources to PM 1 pollution levels. Nevertheless, a long-term analysis of / ratio could still be performed in the future to evaluate the shift of the energy consumption structure over the study region (Shao et al., 2018) .

The / ratios were also found to exhibit significant seasonal variability, with the largest ratio mainly observed in the winter, except in Ürümqi where the largest ratio emerged in the autumn, while the lowest ratio was observed by all sites in the summer (Fig. 3b ). Among these six regions, the largest ratios tripled that of the lowest in cities such as Xi'an, Shanghai, and Guangzhou, particularly in Xi'an where the ratio even exceeded one in the winter. Inter-comparisons of / ratios between summer and winter imply an increased contribution from stationary NO − 3 SO 2− 4 sources other than mobile sources in the winter, since the latter normally show no significant variation across seasons in urban areas. This could be due to the high vulnerability of volatilization of nitrate salts such as NH 4 NO 3 in the summer as well as the fast formation of nitrates under cold and wet weather conditions in the winter (Hu et al., 2011; Liu et al., 2014) . On the other hand, nitrate loss could be also a possible factor influencing the / ratios. It is known that quartz exhibits much higher nitrate loss than other filter media such as Teflon and Nylon (Chow et al., 2005) , especially for samples from subtropical areas during summer. In other words, the nitrate loss could be uneven between different sampling sites and seasons. Figure 4 compares the annual cycle of EC, OM, and SNA in PM 1 mass in these six representative regions. Noteworthy is that the sum of EC, OM, and SNA at QHL only accounted for <30% of PM 1 mass throughout the year, except in December and January, when there was a marked increase in OM. Nevertheless, the average ratio was still low (24.7%), even after accounting for other inorganic ions such as K + , Na + , Ca 2+ , Mg 2+ , Cl − , and F − . This implies the spe- Fig. 4 . Monthly mean relative abundance of elemental carbon (EC), organic mass (OM), and secondary inorganic species (denoted as SNA, which is represented by the sum of sulfate, nitrate, and ammonium) in PM 1 mass at six cities across China during 2012−14. A commonly applied factor of 1.4 was applied to convert organic carbon to OM to account for unmeasured atoms in organic species. Note that observations at the QHL station were temporally discontinuous during the field campaigns due to the absence of valid measurements in several months. cies considered in this study were insufficient to account for PM 1 mass measured at QHL. Conversely, the sum of EC, OM, and SNA constituted more than 80% of PM 1 mass during wintertime in Guangzhou (January) and Xi'an (October, November, and January), especially the latter two components (i.e., OM and SNA), accounting for the vast majority of wintertime PM 1 mass concentration. Overall, the largest fraction was more likely to be observed in January due to a substantial increase in OM and SNA, particularly OM, which increased from about 20% in December to nearly 50% in January in Guangzhou (Fig. 4f) . These results clearly revealed significant variations in PM 1 and its speciation in space and time, highlighting the importance of continuous field campaigns, particularly a comprehensive analysis of elements, to better understand the variability of chemical compositions in PM 1 across China.

To better examine the associations between chemical compositions and PM 1 mass concentrations, particularly their relative roles in elevating PM 1 loading levels, we compared mass concentration of each major chemical composition in PM 1 at each sampling location between distinct pollution episodes. We define clean and polluted episodes according to the relative levels of PM 1 mass concentration in each specific region as days with the smallest and the largest 20% of PM 1 , respectively. Detailed comparison results were presented in Table 2 .

It is found that PM 1 concentration on polluted days was, on average, 2.6~6.5 times the magnitude of that during clean periods, with the largest increase observed in Ürümqi (6.5), followed by Qinghai (4.6), Shanghai (3.1), Guangzhou (2.8), Chengdu (2.7), and Xi'an (2.7). This also implies PM 1 concentration in Ürümqi experienced the most significant variation over the time covered by the sampling period. On the contrary, PM 1 concentration in Chengdu exhibited a relatively low variability throughout the year. In regard to chemical constituents, the mass concentration of some species during polluted episodes increased by a factor as large as 10 times of that during clean episodes.

During polluted periods, in Qinghai, sulfate almost tripled in concentration compared to that on clean days, while OM, EC, and ammonium doubled their concentrations. In contrast, nitrate maintained a nearly constant concentration. This implies the high PM 1 concentration at QHL was dominated by primary combustion emissions, since emissions from vehicles were negligible around the sampling region. In Ürümqi, there were substantial increases in concentrations of five major chemical compositions as ammonium, nitrate, sulfate, OM, and EC increased by a factor of 11.9, 10.9, 9.5, 4.4, 2.7, respectively. Such a salient increase in these five major compositions inevitably elevated the mass concentration of PM 1 . Given the considerable increase in concentrations of SNA, we may largely ascribe the drastic elevation in PM 1 mass concentration during polluted days in Ürümqi to extensive anthropogenic primary emissions. On the other hand, this also could be attributable to enhanced nitrate formation as well as aerosol hygroscopicity due to increased RH during wintertime in Ürümqi resulting from extensive snowfalls, since a similar effect was also observed in north China (Shao et al., 2018; Sun et al., 2020) .

Xi'an and Chengdu, two inland mega cities, were found to exhibit similar variations in concentrations of both PM 1 and major chemical compositions. Moreover, their data magnitudes were also comparable. Specifically, PM 1 concentration on polluted days in both cities was 2.7 times of that on clean days. In terms of the five major chemical constituents studied, Xi'an had larger concentrations of EC, OM, nitrate, and ammonium than Chengdu during polluted episodes, especially the former two species as their concentrations in Xi' an significantly exceeded those in Chengdu. These results suggest that primary emissions played a more important role in elevating PM 1 pollution levels in Xi'an, whereas PM 1 pollution in Chengdu could be attributable to larger emissions from stationary sources (e.g., coal combustion) given higher Table 2 . Comparisons of mean mass concentration of PM 1 and chemical compositions (units: μg m −3 ) during clean and polluted episodes. Here clean and polluted days were defined as days with the lowest and the largest 20% PM 1 mass concentration during the study time period, respectively. N denotes the number of samples in each episode and uncertainty of each factor was characterized by plus and minus one standard deviation from the mean value. The two mega cites in the coastal regions, Shanghai and Guangzhou, were both found to suffer from similar PM 1 pollution levels given comparable mean PM 1 concentration. It is indicative that PM 1 concentration during polluted days almost tripled that on clean days in both cities. During polluted episodes, concentration of all species except EC was observed to increase significantly in both cities. Noteworthy is that all the remaining four major species were found to increase even more markedly in Shanghai, especially OM (+10.6 μg m −3 ) and nitrate (+8.2 μg m −3 ). This implies more important roles of mobile sources in elevating PM 1 pollution levels in Shanghai during haze days. Similar results have also been revealed in the literature (e.g., Shi et al., 2014; Qiao et al., 2016) .

Compared to the other four species, EC constituted a relatively low fraction in urban areas. Also, the absolute concentration exhibited small variation between two regimes. This indicates that the phenomenal increase in PM 1 concentration during polluted days should be attributable more to other species such as OM, nitrate and sulfate as their concentrations were found to increase drastically. This finding also corroborates the critical role of the fast formation of secondary aerosols via aqueous reaction in exacerbating regional fine particle pollution levels during polluted periods (Wang et al., 2017 (Wang et al., , 2019 Yao et al., 2018; Dang and Liao, 2019) .

In light of the increases in mean concentrations of PM 1 and five major chemical species, we find that mean PM 1 concentration had increased by 45~178 μg m −3 during polluted periods compared to that on clean days (Table 2) . However, the accumulated increases in concentration of the five major chemical constituents only accounted for a small fraction of the total increased PM 1 concentration. This implies that the observed high PM 1 loading also should be attributable to factors other than these five major chemical constituents. In other words, it is challenging to resolve PM 1 pollution mechanisms simply based on these five chemical species since they only accounted for a small fraction of the observed increase in PM 1 concentration.

Since severe haze pollution events are oftentimes observed in the winter, we further compared mass concentrations of PM 1 and five major chemical constituents during polluted days in winter and other seasons. As shown in Table 3 , higher mean PM 1 concentration was observed during wintertime in Xi'an, Shanghai, and Guangzhou. In Ürümqi, there was a noticeable increase in concentration of OM, sulfate, and ammonium in the winter, especially OM (+12.61 μg m −3 ) as the mean wintertime concentration was 1.5 times of that in other seasons. This implies more extensive primary emissions in Ürümqi in the winter and that the elevated OM could be attributable to wetter weather due to extensive wintertime snowfalls that boost the formation of SOA.

In Xi'an, all five major chemical compositions were found to have larger mass concentration during the wintertime except OM, especially EC (+4.01 μg m −3 ), implying the wintertime PM 1 pollution in Xi'an could be ascribed to more extensive primary emissions. In contrast, both PM 1 and the five major chemical constituents were observed to have higher concentration in Chengdu during the wintertime, with the most significant increase found in OM (+12.31 μg m −3 ), followed by nitrate (+5.45 μg m −3 ), ammonium (+1.53 μg m −3 ) and EC (+1.36 μg m −3 ). This could be due to the enhanced formation of secondary aerosols during the wintertime. In contrast, little difference was observed in chemical compositions between seasons in Guangzhou. Noteworthy is that there was a ubiquitous increase in ammonium during wintertime in these five major cities, though the magnitude was relatively small. These diverse results clearly highlight the complexity behind the formation mechanism of PM 1 pollution events as factors like the formation of SOA and aerosol hygroscopicity could all result in a significant increase in PM 1 loading in the atmosphere (Khan et al., 2021) .

In addition to chemical speciation, PM 1 mass concentration was modulated by meteorological conditions, which are believed to play an important role in governing the disper- Table 3 . Same as in Table 2 sion of atmospheric particles and its gaseous precursors. Figure 5 shows the correlation between PM 1 mass concentration and the four meteorological factors of air temperature, wind speed, RH, and BLH, respectively. Negative correlations were observed between PM 1 mass concentration and air temperature (except at the QHL where exhibited a positive correlation), especially in Ürümqi, Xi'an, and Shanghai, where negative correlations were even statistically significant at the 95% confidence interval. High PM 1 concentrations were mainly observed during cold days, whereas low values were prevalent on hot days. This effect was even more pronounced at Ürümqi and Xi' an, in large part possibly due to the wintertime heatingrelated extensive coal combustion in these two cities. In addition, stagnant weather was more likely to occur in the winter, which can also elevate PM 1 loading (Shao et al., 2018) . For regions like Guangzhou, temperature appeared to have no significant association with PM 1 concentration. This is because temperature is normally high in Guangzhou throughout the year and there is almost no extensive wintertime heating related coal combustion there.

Wind speed exhibited a negative correlation with PM 1 mass concentration in the five major cities (Fig. 5b) , which is also in line with our expectation, since wind modulates the diffusion and transport of regional air pollutants. Nevertheless, significant correlations were found only in Ürümqi and Chengdu, indicating wind played critical roles in modulating regional PM 1 concentration levels in these two cities.

Humidity is also a critical factor since it determines the nucleation and hygroscopic growth of particles through gasliquid reaction process, especially the aqueous-phase oxidation of SO 2 (Wang et al., 2016 (Wang et al., , 2020 Li et al., 2017; Zheng et al., 2017; Sun et al., 2020) . Given the complex process of gas-liquid reaction, the relationship between humidity and PM 1 mass concentration is oftentimes nonlinear and nonstationary, making it challenging to characterize the statistical associations between them. A statistically significant positive correlation was found between RH and PM 1 in Ürümqi, with high pollution levels being more likely to be observed under wetter conditions with RH > 70% (Fig. 5c) . A recent study also revealed distinct responses of primary and secondary species in PM 1 and PM 2.5 to RH due to changes in aerosol hygroscopicity and phase states (Sun et al., 2020) .

BLH is another key factor modulating the dispersion and accumulation of air pollutants in the troposphere in addition to wind speed (Guo et al., 2019; Lou et al., 2019) , since a low planetary boundary layer can suppress the vertical and horizontal diffusion of aerosol, in turn exacerbating the haze pollution Miao et al., 2018; Zhu et al., 2018; Bai et al., 2019a) . Negative correlations were found between BLH and PM 1 mass concentration in five major cities, with statistically significant correlations (95% confidence interval) observed in Ürümqi (R = -0.35), Xi'an (R = -0.48), and Shanghai (R = -0.58), respectively. It is understood that severe haze pollution events oftentimes are accompanied by a low BLH (<200 m), causing both fine particles and precursors to be suppressed in the lower troposphere near the ground, and accelerating the accumulation of air pollutants and higher PM 1 loading in turn. Table 4 compares the seasonal mean meteorological conditions during polluted days between winter and other seasons. The results indicate that there were significant yet heterogeneous differences in meteorological conditions between distinct seasons. Generally, the synoptic weather in the winter was oftentimes cold, calm (lower wind speed), along with a low BLH, and even wetter in some regions (e.g., Ürümqi, Chengdu, and Shanghai). Evidently, such a weather condition is favorable to the accumulation of air pollutants (low wind speed and BLH) and the formation of secondary aerosols (cold and wet) (Li et al., 2017; Wang et al., 2017; An et al., 2019) . Noteworthy is that synoptic weather appears to have little effect on PM 1 pollution levels in Chengdu despite the similarity in PM 1 concentration and meteorological conditions (except temperature) across seasons. This could be attributable to the basin topography of Chengdu, where air pollutants cannot be easily diluted under a relatively calm weather regime, thereby yielding a near constant ambient PM 1 loading therein. Figure 6 compares the relative importance of each factor in predicting PM 1 concentration levels at each individual sampling site. As shown, leading factors (i.e., factors with larger relative importance) differed across study regions. Specifically, EC was the leading factor in predicting PM 1 mass concentration at the QHL, implying the variability of PM 1 concentration was highly associated with the loading of EC. Since this region was deemed as the ambient background given limited human activity, it is reasonable to con-clude that PM 1 loading in this region is largely modulated by EC from primary emissions occurring there. In Ürümqi, EC and OC were found to play more important roles, indicative of larger contribution of primary sources to PM 1 variations in this region.

In Xi'an, OC played the predominant role compared to other factors, followed by RH, EC, and nitrate. This is indicative of the significant role of secondary aerosols in modulating PM 1 mass concentration there. In contrast, OC and SNA played comparable roles in Chengdu, followed by wind speed. Among them, sulfate had a larger importance than the others, implying more important roles played by stationary sources. Moreover, wind was also a key factor in modulating the dispersion of PM, since air pollutants cannot be easily diluted in Chengdu due to the basin topography. In Shanghai, OC and wind speed were two leading factors, followed by BLH, indicating the critical role of secondary aerosols and meteorological conditions like temperature inversion in modulating PM 1 concentration in this region.

In spite of the low abundance of ammonium and nitrate as compared to other three species, these two species were found to play the most critical roles in modulating the variations in PM 1 mass concentration in Guangzhou (Fig. 6f) . Conversely, RH was found to play the least important role, because Guangzhou is located in Pear River estuary, which has a humid subtropical climate, thereby the humidity is high year round. In other words, humidity would not be a key factor limiting the formation of secondary aerosols in a humid region. Ammonium was found to play the most critical role in modulating the variations in PM 1 mass concentration, followed by nitrate. This could be related to the neutralization of by since , , and display positive feedback interactions among these constituents, resulting in elevated PM 1 concentration in turn (Shao et al., 2018) . The second leading role of nitrate further corroborates the critical role of mobile sources in contributing to PM 1 pollution levels in Guangzhou, which is also in line with previous findings . These results clearly indicate the relatively important role of secondary aerosols in modulating PM 1 loading variability in Guangzhou. To better understand PM pollution in these six regions, we also compared the mass concentration of PM and five major chemical compositions between PM 1 and PM 2.5 (using PM 2.1 as a proxy). Given the simultaneous sampling, we calculated their mass concentration ratios (PM 1 /PM 2.5 ) and the results are summarized in Table 5 . It is seen that PM 1 accounted for 66%−75% mass concentration of PM 2.5 that were measured in these six regions, which is slightly higher than that in Tianjin (63%) as discovered in Khan et al. (2021) . Such a high PM 1 /PM 2.5 ratio implies PM 2.5 in these six regions was largely dominated by small size particles due to combustion processes and secondary aerosols rather than natural sources and mechanical processes (Khan et al., 2021) . Meanwhile, this also highlights the notorious negative health impacts of PM 2.5 since nearly 70% of its composition is from submicron particles.

Spatially, PM 1 accounted for a larger fraction of PM 2.5 mass concentration in the two coastal megacities of Shang-hai (0.73) and Guangzhou (0.75), where primary emissions from vehicles and ships are normally more extensive than in the two inland megacities of Xi'an (0.66) and Chengdu (0.69). Such a spatial distribution of PM 1 /PM 2.5 ratio is also in line with previous studies in Chen et al. (2018) as high ratios were mainly observed in Southeastern and Central China while low ratios were found in Western China where natural sources are also important contributions. Moreover, this ratio varied with pollution regimes and humidity. For instance, Qiao et al. (2016) found a relatively high PM 1 /PM 2.5 ratio on clean days whereas an elevated ratio of PM 1−2.5 /PM 2.5 on hazy days. Sun et al. (2020) also revealed that the changes in PM 1 /PM 2.5 ratios as a function of RH were largely different for primary and secondary aerosol species due to distinct aerosol hygroscopicity and phase states.

In regard to the five major chemical compositions, we found that most major PM 2.5 constituents in Guangzhou were PM 1 given higher ratios of chemical species, except nitrate and ammonium, which had the largest fraction in Ürümqi (0.72) and Qinghai (0.94), respectively. This Table 5 . Same as in Table 2 implies PM 2.5 pollution in Guangzhou is mainly dominated by submicron fine particle pollution events, followed by Shanghai given the second highest ratio. On the contrary, Xi'an and Chengdu are two cities having relatively low fractions in chemical compositions. This result suggests even more diverse emission sources in these two inland cities, since a larger fraction of chemical compositions in PM 2.5 was contributed by coarser particles. Overall, our results indicate that the main sources of PM 1 and PM 2.5 in the two coastal cities were dominated by combustion processes and secondary aerosols while natural sources and mechanical emissions could contribute significantly in the two inland cities.

Based on three-year-long aerosol samplings from field campaigns, we have examined the variations of chemical components in PM 1 in six representative regions across China. The PM 1 mass concentration and its chemical constituents were compared with respect to pollution levels, seasons, and meteorological conditions. The results revealed that both PM 1 and its chemical constituents varied significantly in space and time, with high PM 1 loading observed mainly in the winter and low values in the summer. Meanwhile, a substantial increase in OM and SNA was oftentimes observed in elevating PM 1 mass concentration, particularly OM during wintertime. Compared to other species, OM had a larger overall abundance in PM 1 , indicative of a more important role of primary emission and the formation of secondary aerosols in exacerbating PM 1 pollution levels during wintertime.

NO − 3 SO 2− 4 Comparisons of OC/EC and / ratios across regions revealed that PM 1 loading in Guangzhou tended to be modulated largely by primary emissions given lower ratios of both factors throughout the year. The largest ratios of both factors in Xi'an were observed mainly in the winter, with the lowest observed in the summer. This effect could be associated with the seasonality of both emission intensity and meteorological conditions. Moreover, elevated concentrations of SNA and OM were oftentimes observed during polluted days, particularly nitrate and sulfate, though their increments differed greatly across the regions. This implies important roles of both primary emission and secondary particle formation in producing haze events. In spite of a small magnitude, a ubiquitous increase in ammonium was observed in the five major cities during wintertime, and this was found to play the most important role in modulating the variability of PM 1 mass concentration in Guangzhou.

Meteorological factors were shown to play important roles in modulating the variability of PM 1 mass loading. Our modeling practices have clearly revealed that ambient PM 1 mass concentration was largely regulated by primary emissions (e.g., Qinghai and Ürümqi) and secondary aerosol formations (e.g., Xi'an, Shanghai, and Chengdu). Meanwhile, the importance of meteorological factors was found to vary across regions under different pollution regimes. Overall, the results in this study clearly indicate that there is no common driving mechanism for haze events due to diverse sources and influential factors. Due to the data issue, an elemental analysis was not performed in the current study, but this should be included in future field campaigns to more fully investigate the driving factors behind haze events.

Permutation importance: A corrected feature importance measure

Severe haze in northern China: A synergy of anthropogenic emissions and atmospheric processes

Visibility as related to atmospheric aerosol constituents

Relationships among aerosol constituents from Asia and the North Pacific during PEM-West A

Seasonal estimation of organic mass to organic carbon in PM 2.5 at rural and urban locations in New York state

Diagnosing atmospheric stability effects on the modeling accuracy of PM 2.5 /AOD relationship in eastern China using radiosonde data

Advancing the prediction accuracy of satellite-based PM 2.5 concentration mapping: A perspective of data mining through in situ PM 2.5 measurements

Global estimates of mortality associated with long-term exposure to outdoor fine particulate matter

Atmospheric visibility related to aerosol mass concentration

Spatiotemporal variation of PM 1 pollution in China

Equivalence of elemental carbon by thermal/optical reflectance and transmittance with different temperature protocols

Loss of PM 2.5 nitrate from filter samples in central

The IMPROVE_ A temperature protocol for thermal/optical carbon analysis: Maintaining consistency with a long-term database

Severe winter haze days in the Beijing−Tianjin−Hebei region from 1985 to 2017 and the roles of anthropogenic emissions and meteorology. Atmospheric Chemistry and Physics

New evidence on the impact of sustained exposure to air pollution on life expectancy from China's Huai River Policy

Declining frequency of summertime local-scale precipitation over eastern China from 1970 to 2010 and its potential link to aerosols

Shift in the temporal trend of boundary layer height in China using long-term (1979− 2016) radiosonde data

Elucidating severe urban haze formation in China

The random subspace method for constructing decision forests

AERONET-A federated instrument network and data archive for aerosol characterization

Hygroscopicity and evaporation of ammonium chloride and ammonium nitrate: Relative humidity and size effects on the growth factor

High secondary aerosol contribution to particulate pollution during haze events in China

Enhanced secondary pollution offset reduction of primary emissions during COVID-19 lockdown in China

2021: Chemical characterization and source apportionment of PM 1 and PM 2.5 in Tianjin, China: Impacts of biomass burning and primary biogenic sources

Comparison of abundances, compositions and sources of elements, inorganic ions and organic compounds in atmospheric aerosols from Xi'an and New Delhi, two megacities in China and India

Aerosol and boundary-layer interactions and impact on air quality

East Asian study of tropospheric aerosols and their impact on regional clouds, precipitation, and climate (EAST-AIR CPC )

Aerosol hygroscopicity derived from size-segregated chemical composition and its parameterization in the North China Plain

Elucidating the relationship between aerosol concentration and summertime boundary layer structure in central China

On the relationship between aerosol and boundary layer height in Summer in China under different thermodynamic conditions

Unraveling the relationships between boundary layer height and PM 2.5 pollution in China based on fouryear radiosonde measurements

Nature and climate effects of individual tropospheric aerosol particles

Simultaneous monitoring and compositions analysis of PM 1 and PM 2.5 in Shanghai: Implications for characterization of haze pollution and source apportionment

Impacts of traffic emissions on atmospheric particulate nitrate and organics at a downwind site on the periphery of Guangzhou

Seasonal variation and size distribution of biogenic secondary organic aerosols at urban and continental background sites of

Characterizing remarkable changes of severe haze events and chemical compositions in multi-size airborne particles

Airborne submicron particulate (PM 1 ) pollution in Shanghai, China: Chemical variability, formation/dissociation of associated semi-volatile components and the impacts on visibility

2020: Chemical differences between PM 1 and PM 2.5 in highly polluted environment and implications in air pollution studies

Species contributions to PM 2.5 mass concentrations: Revisiting common assumptions for estimating organic mass

Evolution of aerosol chemistry in Xi'an, inland China, during the dust storm period of 2013-Part 1: Sources, chemical forms and formation mechanisms of nitrate and sulfate

Persistent sulfate formation from London Fog to Chinese haze

Particle acidity and sulfate production during severe haze events in China cannot be reliably inferred by assuming a mixture of inorganic salts. Atmospheric Chemistry and Physics

Fast sulfate formation from oxidation of SO 2 by NO 2 and HONO observed in Beijing haze

Distinct ultrafine-and accumulation-mode particle properties in clean and polluted urban environments

New particle formation in China: Current knowledge and further directions

Strategy for the implementation of the global atmosphere watch programme

The characteristics of atmospheric brown carbon in Xi'an, Inland China: Sources, size distributions and optical properties

The campaign on atmospheric aerosol research network of China: CARE-China

Characteristics of PM 2.5 speciation in representative megacities and across China. Atmospheric Chemistry and Physics

Atmospheric new particle formation from sulfuric acid and amines in a Chinese megacity

The water-soluble ionic composition of PM 2.5 in Shanghai and Beijing

2020: Stronger policy required to substantially reduce deaths from PM 2.5 pollution in China

Wintertime aerosol properties in Beijing

Characteristics of the chemical components of aerosol particles in the various regions over China

Evidence of major secondary organic aerosol contribution to lensing effect black carbon absorption enhancement

Aerosol characteristics and impacts on weather and climate over the Tibetan Plateau

Analysis of influential factors for the relationship between PM 2.5 and AOD in Beijing. Atmospheric Chemistry and Physics

Mixing layer height on the North China Plain and meteorological evidence of serious air pollution in southern Hebei

Air pollution intervention and life-saving effect in China

Acknowledgements. We are grateful to three anonymous reviewers and editors for their copious and constructive comments and suggestions. This work was financially supported by National Key R&D Plan (NO. 2017YFC0210000), National Natural Science Foundation of China (NO. 41701413), National Key R&D Plan (No. 2017YFC0212703) and Strategic Priority Research Program of the Chinese Academy of Sciences (NO. XDB05020401). Meteorological data were acquired from the Meteorological Information Comprehensive Analysis and Process System (air temperature, relative humidity, and wind speed), and ERA-Interim reanalysis (boundary layer height) that was provided by the European Centre for Medium-Range Weather Forecasts.