Journal cover Journal topic
Atmospheric Chemistry and Physics An interactive open-access journal of the European Geosciences Union
Journal topic
Atmos. Chem. Phys., 18, 14637–14651, 2018
https://doi.org/10.5194/acp-18-14637-2018
Atmos. Chem. Phys., 18, 14637–14651, 2018
https://doi.org/10.5194/acp-18-14637-2018

Research article 12 Oct 2018

Research article | 12 Oct 2018

# Aerosol chemistry and particle growth events at an urban downwind site in North China Plain

Aerosol chemistry and particle growth events at an urban downwind site in North China Plain
Yingjie Zhang1,2, Wei Du1,3, Yuying Wang4, Qingqing Wang1, Haofei Wang5, Haitao Zheng1, Fang Zhang4, Hongrong Shi6, Yuxuan Bian7, Yongxiang Han2, Pingqing Fu1,3, Francesco Canonaco8, André S. H. Prévôt8, Tong Zhu9, Pucai Wang6, Zhanqing Li4, and Yele Sun1,2,3,10 Yingjie Zhang et al.
• 1State Key Laboratory of Atmospheric Boundary Layer Physics and Atmospheric Chemistry, Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100029, China
• 2Collaborative Innovation Center on Forecast and Evaluation of Meteorological Disasters, Nanjing University of Information Science & Technology, Nanjing 210044, China
• 3College of Earth Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
• 4College of Global Change and Earth System Science, Beijing Normal University, Beijing 100875, China
• 5Key Laboratory of 3D Information Acquisition and Application of Ministry of Education, Capital Normal University, Beijing 100048, China
• 6Key Laboratory of Middle Atmosphere and Global Environment Observation, Institute of Atmospheric Physics, Chinese Academy of Sciences, Beijing 100029, China
• 7State Key Laboratory of Severe Weather, Chinese Academy of Meteorological Sciences, Beijing 100081, China
• 8Laboratory of Atmospheric Chemistry, Paul Scherrer Institute, Villigen PSI 5232, Switzerland
• 9College of Environmental Sciences and Engineering, Peking University, Beijing 100871, China
• 10Center for Excellence in Regional Atmospheric Environment, Institute of Urban Environment, Chinese Academy of Sciences, Xiamen 361021, China

Correspondence: Yele Sun (sunyele@mail.iap.ac.cn)

Abstract

The North China Plain (NCP) has experienced frequent severe haze pollution events in recent years. While extensive measurements have been made in megacities, aerosol sources, processes, and particle growth at urban downwind sites remain less understood. Here, an aerosol chemical speciation monitor and a scanning mobility particle sizer, along with a suite of collocated instruments, were deployed at the downwind site of Xingtai, a highly polluted city in the NCP, for real-time measurements of submicron aerosol (PM1) species and particle number size distributions during May and June 2016. The average mass concentration of PM1 was 30.5 (±19.4) µg m−3, which is significantly lower than that during wintertime. Organic aerosols (OAs) constituted the major fraction of PM1 (38 %), followed by sulfate (25 %) and nitrate (14 %). Positive matrix factorization with the multilinear engine version 2 showed that oxygenated OA (OOA) was the dominant species in OA throughout the study, on average accounting for 78 % of OA, while traffic and cooking emissions both accounted for 11 % of OA. Our results highlight that aerosol particles at the urban downwind site were highly aged and mainly from secondary formation. However, the diurnal cycle also illustrated the substantial influence of urban emissions on downwind sites, which are characterized by similar pronounced early morning peaks for most aerosol species. New particle formation and growth events were also frequently observed (58 % of the time) on both clean and polluted days. Particle growth rates varied from 1.2 to 4.9 nm h−1 and our results showed that sulfate and OOA played important roles in particle growth during clean periods, while OOA was more important than sulfate during polluted events. Further analyses showed that particle growth rates have no clear dependence on air mass trajectories.

1 Introduction

Atmospheric aerosols can reduce visibility (Zhang et al., 2010), have adverse effects on human health (Chen et al., 2013), and also affect radiative forcing directly by absorbing and scattering solar radiation and indirectly by modifying cloud formation and properties (Boucher et al., 2013). According to the latest report on global urban air quality by the World Health Organization (WHO), the top 10 most polluted cities in China are all located in the North China Plain (NCP). The concentration of particulate matter (PM) less than or equal to 2.5 µm in diameter (PM2.5) in Xingtai was ranked the highest in China in 2014 (http://apps.who.int/gho/data/view.main.AMBIENTCITY2016?lang=en, last access: 24 September 2017). Although extensive studies have characterized the formation mechanisms and evolution processes of haze in the NCP (Quan et al., 2011; Zhao et al., 2013; Yang et al., 2015; Li et al., 2017), high anthropogenic emissions and stagnant meteorological conditions are the major factors leading to severe PM pollution (R. Zhang et al., 2015; Fu and Chen, 2017; Guo et al., 2014; Sun et al., 2014). Mitigating air pollution in the NCP remains a challenge. One reason is the complexity of ambient aerosols, which have largely different compositions and come from different sources from different regions and cities.

Xingtai, one of the most polluted cities in China, had an urgent front-burner environmental problem. PM sources are dynamic and include local emissions, e.g., biomass burning, traffic, and cooking emissions, and the transport of pollutants from upwind (east and south) polluted areas (Fu et al., 2014). As a result, both local and regional sources contribute to high concentrations of PM, leading to certain uncertainties in air quality control. Although the annual average concentration of PM2.5 in Xingtai decreased from 160 µg m−3 in 2013 to 87 µg m−3 in 2016 (the data are from four monitoring sites in urban Xingtai that was released by the China National Environmental Monitoring Centre), it still far exceeds the Chinese National Air Quality Standard (35 µg m−3 for an annual average) and that of the WHO (10 µg m−3). In addition, concentrations of NO2 and CO changed little from 2013 to 2016 while that of SO2 decreased substantially (Fig. S1). As a response to the changes in precursors, aerosol particle composition may also change significantly. Therefore, characterization of the composition, sources, and processes of PM in regions near Xingtai is important to do so that effective strategies for future air quality improvements can be provided. Previous studies carried out in Xingtai have investigated the frequency of haze events (Fu et al., 2014), ammonia emissions (Zhou et al., 2015), and the sources of PM2.5 (L. Wang et al., 2015). The results highlighted the importance of both local (especially industrial) and regional sources of air pollution in Xingtai. However, real-time characterizations of aerosol composition and particle number size distributions have not yet been reported. A recent study conducted in a similarly polluted city, Handan, which is approximately ∼50 km south of Xingtai, showed significant contributions of coal and biomass combustion to haze formation in winter (Li et al., 2017). However, aerosol characteristics in Xingtai are not well known and the impacts of urban emissions on downwind sites also remain poorly understood.

Figure 1(a) Location of the sampling site. (b) Average diurnal evolution of wind vector. The pie chart in (a) shows the average aerosol composition for the entire study. The two arrows in (a) show the daytime and nighttime prevailing wind directions.

In this study, an Aerodyne aerosol chemical speciation monitor (ACSM) along with a suite of collocated instruments was deployed at the downwind site of Xingtai from 30 April to 20 June 2016 to characterize aerosol chemistry and particle growth events in spring and summer. The mass concentrations, chemical composition, and temporal and diurnal variations of submicron aerosol (PM1) species are characterized, and the sources of organic aerosols (OA) are investigated with positive matrix factorization (PMF) and a FLEXible PARTicle dispersion model (FLEXPART) analysis. Also, particle growth events and their relationship to aerosol chemistry are discussed.

2 Experimental methods

## 2.1 Sampling site

Xingtai is located in the central-south part of the Beijing–Tianjin–Hebei region with the Taihang Mountains to the west (Fig. 1a). In this work, all measurements were made at the Xingtai National Meteorological Basic Station (XNMBS), a suburban site located approximately 17 km northwest of Xingtai City (37.18 N, 114.37 E; 180 m a.s.l.), from 30 April to 20 June 2016. The sampling site was influenced by mountain–plain winds during the study period. As shown in Fig. 1b, the wind direction showed clear day and night patterns with prevailing south-southeasterly winds in the day and west-northwesterly winds at night. The average wind speed was 4 m s−1 and the average temperature was 22.6 C during the study period.

## 2.2 Measurements

All instruments were placed in a container at the sampling site. The non-refractory PM1 (NR-PM1) chemical components including sulfate (SO4), nitrate (NO3), ammonium (NH4), chloride (Chl), and organics (Org) were measured in situ by an ACSM at a time resolution of 5 min. The ACSM was operated in the same way as in previous studies (Sun et al., 2016b; Zhang et al., 2016). A PM2.5 cyclone (Model: URG-2000-30ED) was supplied in front of the sampling inlet to remove coarse particles larger than 2.5 µm. The ambient air was drawn into the container through a 0.5 in. (outer diameter) stainless steel tube at a flow rate of 3 L min−1 using an external pump, of which ∼0.1 L min−1 was sub-sampled into the ACSM. The sampling height was approximately 2 m, and the particle residence time in the sampling tube was about 5 s. Aerosol particles were then dried by a silica gel diffusion dryer before sampling into the ACSM. Before the campaign, the ACSM was calibrated with pure ammonium nitrate particles following the standard protocols in Ng et al. (2011b). Because the ACSM does not measure refractory black carbon (BC), a seven-wavelength Aethalometer (model AE33, Magee Scientific Corp.; Drinovec et al., 2015) was used to measure BC.

The size-resolved particle number concentrations in the size range from 15 to 685 nm were measured in situ by a condensation particle counter (CPC, model 3775, TSI) equipped with a long differential mobility analyzer (DMA, model 3081A, TSI). The time resolution is 5 min. The total number concentrations (7–2000 nm) were measured by a mixing condensation particle counter (MCPC, model 1720, Brechtel). Other collocated measurements included the light extinction of dry PM2.5 at 630 nm measured by a cavity attenuated phase shift extinction monitor (CAPS PMext; Massoli et al., 2010); the mass concentration of PM2.5 measured by a six-channel particle counter (manufactured by Graywolf); and gaseous species of CO, NO, NOx, SO2, and O3 measured by gas analyzers (manufactured by ECOTECH). Meteorological parameters including ambient temperature (T), relative humidity (RH), wind speed (WS), wind direction (WD), precipitation, and solar radiation were also measured at the same site by the Xingtai Meteorological Administration.

## 2.3 Data analysis

### 2.3.1 ACSM data analysis

ACSM data were analyzed using the ACSM standard software (version 1.5.3.0) within the Igor Pro software environment (Wave Metrics, Inc., Oregon, USA). The default relative ionization efficiencies for all species except NH4 were used in the study. That for NH4 was determined from the ionization efficiency calibration. A collection efficiency (CE) of 0.5 was used to account for the incomplete detection of aerosol species, mainly due to particle bounce at the vaporizer (Matthew et al., 2008). The CE can be composition dependent and especially sensitive to the fraction of ammonium nitrate (NH4NO3) as well as affected by particle acidity and RH (Middlebrook et al., 2012). In this study, SO4 dominated inorganic species and the average contribution of NH4NO3 was 18 % (maximum: 38 %), which would not affect CE substantially. Aerosol particles were slightly acidic, as indicated by the average ratio (∼0.7) of measured NH4 to predicted NH4 that is required to fully neutralize SO4, NO3, and Chl (Zhang et al., 2007), which are also not acidic enough to affect the CE. Default relative ionization efficiencies (RIEs) were used except for ammonium (5.0) and sulfate (0.98), which were determined from pure ammonium nitrate and ammonium sulfate, respectively. In addition, to reduce the influence of RH on CE, a silica gel diffusion dryer was deployed to keep the RH in the sampling line below 40 %. In fact, the differences in mass concentrations were less than 5 % between composition-dependent CE and a constant CE of 0.5 in this study.

Figure 2Mass spectral profiles (a) and time series of the mass concentrations of three OA factors (b), i.e., HOA, COA, and OOA. Time series of the mass concentrations of BC and sulfate (right axis) are also shown.

Figure S2 shows the comparison between the total PM1 mass (equal to NR-PM1+ BC) and particle volume concentrations measured by the SMPS. Particle volume concentrations were highly correlated with total PM1 mass concentrations (r2=0.77, slope = 0.51). We then estimated the particle density using the chemical composition of PM1 (Du et al., 2017). The average density during the study period was 1.5 g cm−3. Assuming spherical particles, the calculated SMPS mass reports 75 % of the total PM1 mass. Such a difference may be caused by measurement uncertainties between different instruments, the effects of particle shape, and the uncertainties in estimating particle density.

### 2.3.2 Positive matrix factorization (PMF) analysis

To determine the sources of OA, ACSM mass spectra were processed using the Multilinear Engine version 2 (ME-2) algorithm implemented with the toolkit called Source Finder (Canonaco et al., 2013). The so-called a value approach allows the user to introduce a priori information in the form of known factor profiles or time series to obtain a unique solution and thus reduce the rotational ambiguity of the PMF algorithm. The mass spectra and error matrices of OA were prepared according to the procedures detailed by Ulbrich et al. (2009) and Zhang et al. (2011). Given the interference of the internal standard of naphthalene at mz 127–129 and the low signal-to-noise ratios of larger ions, we only considered mz values below 120 in this study. A reference hydrocarbon-like OA (HOA) profile, which is the average of multiple ambient data sets (Ng et al., 2011), and a reference cooking OA (COA) profile in Beijing (Sun et al., 2013) were introduced to constrain the model performance with a-values varying from 0 to 1. Following the guidelines presented by Canonaco et al. (2013) and Crippa et al. (2014), an optimal solution involving three factors with an a-value of 0.2 was accepted. Some important criteria for selecting the optimal solution with a values varying from 0 to 1 are shown in Figs. S3–S7. The mass spectra and time series of three OA factors are shown in Fig. 2.

The HOA factor has a similar mass spectrum to that of freshly emitted traffic or other fossil-fuel combustion aerosols with major peaks at mz equal to 41, 43, 55, and 57. HOA was moderately correlated with BC (r2=0.42). The COA resolved in Xingtai had a mz 55∕57 ratio of 2.3, within the range of values for COA ($\sim >\mathrm{1.2}$) (Mohr et al., 2012). We further evaluated the factors of HOA and COA assuming that BC is predominantly from traffic emissions while the contribution from cooking emissions is minor. POA from the two-factor solution was highly correlated (r2>0.66), with BC between 01:00 and 10:00 when cooking emissions were not significant (Fig. S8). The ratios of POA  BC were also the lowest during this period, suggesting the dominant contribution of HOA to POA. We then used the average ratio of POA  BC (0.62) during this period to estimate the concentrations of HOA and COA. The estimated HOA and COA both contributed 11 % to OA on average, consistent with the results of the ME-2 analysis. This suggested that the results from ME-2 analysis are reasonable.

The mass spectrum of oxygenated OA (OOA) is characterized by a prominent peak at mz 44 (23.6 % of the total OOA signal), which has also been reported in previous studies. In addition, OOA was highly correlated with sulfate (r2=0.75) and moderately correlated with nitrate (r2=0.54), suggesting that OOA is a surrogate of secondary OA (SOA; Fig. 2). We also performed a PMF analysis by applying the PMF2 algorithm to the ACSM-measured OA. Although the two-factor solution identifies a primary OA (POA) and an OOA, solutions with three to five factors show a splitting and mixing of factors. Therefore, the ME-2 algorithm was used in this study.

Figure 3Time series of (a) temperature (T, in black) and relative humidity (RH, in magenta), (b) wind direction (WD, in orange) and wind speed (WS, in black), (c) particle extinction coefficient (Ext., in orange) and precipitation (Precip., in purple), (d) particle number size distribution, (e) mass concentrations of Org, NO3, SO4, NH4, Chl, and BC, and (f) mass fractional contribution of chemical species to total PM1. Polluted events (PEs) and clean periods (CPs) are marked as shaded orange and blue areas, respectively.

## 2.4 Source region analysis

The footprints of the selected episodes were determined using backward simulations from FLEXPART, a Lagrangian transport and dispersion model (Stohl et al., 2005). The model calculated the 36 h backward trajectories of 10 000 particles released every hour from the sampling site at a height of 180 m above sea level. The meteorological data driving the model were simulated by version 3.4 of the Weather Research and Forecasting (WRF) model with a 1 h time resolution and a 10 km spatial resolution. The WRF model was driven by initial and boundary conditions from National Centers for Environmental Prediction global reanalysis data.

The 72 h back trajectories at a height of 500 m were calculated every hour using the Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model (Stein et al., 2015) at the XNMBS. To investigate the chemical characteristics of aerosols from different source regions, a cluster analysis was then performed on the trajectories and three clusters were identified according to their similarities in spatial distributions.

In addition, non-parametric wind regression (NWR; Petit et al., 2017) was performed to evaluate the sources of local emissions and regional transport for PM1 aerosol species and OA factors. The NWR plots represent the probability that a specific compound or source is located in a certain wind direction.

3 Results and discussion

## 3.1 Aerosol composition and temporal variations

The temporal variations of PM1 aerosol species and meteorological parameters (RH, T, WS, WD, and precipitation) are shown in Fig. 3. The average mass concentration of PM1 (equal to NR-PM1+ BC) was 30.5 µg m−3 and ranged from 0.2 to 140.1 µg m−3. Several pollution episodes usually lasting ∼2–3 days were observed during the study period, e.g., on 9–11, 17–23, 28–31 May, and 2–4 June. These pollution episodes were quickly cleaned, mainly by wet scavenging. The temporal variations varied differently among different chemical species. Organics showed dramatic variations, ranging from 0.01 to 101.5 µg m−3, and comprised the major fraction of PM1 for most of the time in this study. High concentration peaks of organics were frequently observed, likely due to the influences of local emissions. By contrast, sulfate concentrations obviously increased and remained relatively high during the pollution events, suggesting the important role of regional transport at the downwind site of Xingtai. The average PM2.5 mass concentration was 45.2 µg m−3. Although the average PM2.5 mass concentration was 15 % lower than that (53.3 µg m−3) measured at the urban sites in Xingtai, it exceeded the Chinese National Ambient Air Quality Standards by 29 %. These results suggest that the urban downwind sites also experience similar PM pollution events as the urban sites. For example, the daily PM2.5 exceeded the Chinese National Ambient Air Quality Standards 24 % of the time during this study. PM1 was highly correlated with PM2.5 (r2=0.95) and on average, comprised ∼64 % of PM2.5. The average mass concentration of PM1 is close to that measured in Xinzhou (35 µg m−3), a city in central China (Wang et al., 2016), but lower than that measured in 2013 in Xianghe (73 µg m−3), a rural site near Beijing (Sun et al., 2016a; Fig. S9). One possible explanation is the significant improvement in air quality during the last 4 years (Fig. S1).

On average, OA was the largest component of PM1, accounting for 38 % of the total PM1 mass, followed by sulfate (25 %), nitrate (14 %), ammonium (10 %), and BC (10 %; Fig. 1a). POA (equal to HOA + COA) and SOA (OOA) accounted for 22 % and 78 %, respectively, of the total OA mass. Together, ∼18 % of PM1 was comprised of primary-related materials (POA + BC) and 82 % was from secondary formation (${\mathrm{NO}}_{\mathrm{3}}+{\mathrm{SO}}_{\mathrm{4}}+{\mathrm{NH}}_{\mathrm{4}}+\mathrm{Chl}+\mathrm{SOA}$), indicating that aerosol particles from secondary aerosol formation processes dominated at the downwind site of Xingtai. Compared with aerosol composition in megacities in the NCP, e.g., Beijing (Hu et al., 2016), aerosol composition in this study showed substantially higher contributions of SOA (29 %) and BC (10 %), while the contributions of NO3 (14 %) and COA (4 %) were low due to the lesser amount of local traffic and cooking emissions. The nitrate contribution was similar to that observed at a suburban site in Xinzhou (Wang et al., 2016), but the SO4 contribution was relatively low (25 % vs. 32 %). One reason is that the higher RH (70 % vs. 52 %) in Xinzhou facilitated the formation of SO4. The PM2.5 concentration during this study period was more than twice as low as that in winter (151 µg m−3, Fig. S10), suggesting different sources between spring–summer and winter. For example, the winter season has significantly enhanced coal combustion emissions.

We also investigated the compositional differences between clean periods and polluted events (Fig. 3). Secondary inorganic aerosols (SIA) including SO4 (26 % vs. 21 %), NO3 (15 % vs. 9 %), and NH4 (10 % vs. 7 %) showed enhanced contributions to PM1 during polluted events, while a corresponding decrease was found for organics. For example, the contribution of SIA increased by 14 % from 37 % during clean periods to 51 % during polluted events, while organics and BC decreased by 9 % and 3 %, respectively. Although SOA dominated the OA composition during both polluted events and clean periods, SOA also contributed more to OA during polluted events (78 % vs. 72 %). Comparatively, POA (HOA and COA) was relatively more important than SOA during clean periods. These results suggest that PM at the downwind site of Xingtai was mainly affected by regional transport and secondary formation.

Figure 4Average diurnal cycles of chemical species of PM1 and OA factors during the entire study period, as well as polluted events (PEs) and clean periods (CPs).

Figure 5Bivariate polar plots of PM1 species as a function of wind speed and wind direction: (a) PM1, (b) Org., (c) SO4, (d) NO3, (e) NH4, (f) Chl, (g) BC, (h) HOA, (i) COA, and (j) OOA. The color scales in (a–j) range from 0 to 44, 17, 11, 7.2, 4.5, 1.9, 4.1, 2.0, 2.7, and 12 µg m−3, respectively.

## 3.2 Diurnal patterns

Aerosol species showed distinctly different diurnal patterns in this study (Fig. 4), indicating that the sources and formation processes of PM pollutants were different. The diurnal cycle of PM1 was characterized by peaks at  10:00 and  21:00 local time (LT). The first peak in the late morning was consistent with that of gas pollutants including CO, SO2, and NOx (Fig. S12). Through comparisons with the diurnal cycle of PM2.5 at the urban site in Xingtai (Fig. S13) and daily variations in WD (Fig. 1), the morning peak was mainly associated with the transport of pollutants from urban sites located to the southeast.

The high concentration of PM1 at night was clearly associated with enhanced primary emission concentrations, e.g., HOA and COA, although the shallower boundary layer may have also played a role. It is interesting to note that bext did not show a pronounced nighttime peak as did PM2.5, indicating that the extinction coefficients of primary aerosols were smaller than those of secondary aerosol species (Q. Wang et al., 2015). The diurnal pattern of PM1 only showed a pronounced peak at night during clean periods, consistent with those of HOA and COA. While the diurnal variations support the enhanced roles of primary emissions for PM during clean periods, they also indicate the lesser influences from urban emissions during daytime.

The diurnal pattern of organics overall resembles that of PM1 and was characterized by two pronounced peaks between 09:00–12:00 and 19:00–21:00 LT. However, the three OA factors showed different diurnal cycles. The average diurnal cycle of HOA showed a small morning peak and a pronounced nighttime peak. While the two peaks were comparable during polluted periods, only the pronounced nighttime peak was observed during clean periods. As shown in Fig. S11a, winds were dominantly from the south-southeast during polluted periods and mainly from the west-northwest during clean periods. These results suggest that the high morning HOA peak was mainly caused by the transport of pollutants from urban sites, while local traffic emissions became an important source of HOA at night. This is also consistent with what is seen in the corresponding bivariate polar plot (Fig. 5h), which shows high concentrations of HOA in regions to the south-southeast and to the north. The decrease in HOA during the day was mainly associated with the rising planetary boundary layer height.

The diurnal pattern of COA (Fig. 4h) was also similar to that observed in megacities (Sun et al., 2013; Crippa et al., 2013; Elser et al., 2016; Y. J. Zhang et al., 2015; Hu et al., 2016; Xu et al., 2016) and was characterized by two peaks around the mealtime hours, reflecting the influence of cooking emissions. However, the COA concentration in this study was much lower than that reported in megacities (e.g., peak concentration: 2.5 µg m−3 vs. 14 µg m−3 in Beijing; Hu et al., 2016), suggesting much lower cooking emissions at the suburban site. The corresponding bivariate polar plot further shows that the high concentration of COA mainly originated from restaurants and inhabitants' activities to the southeast of the site (Fig. 5i).

Figure 6Time series of particle number concentrations for (a) N7–15 calculated from the differences between MCPC and SMPS measurements, (b) N15–40 (15–40 nm), (c) N40–100 (40–100 nm), (d) N100–685 (100–685  nm), and (e) all particles, N15–685 (15–685 nm). The gap in (a) is mainly due to the malfunction of MCPC during this period.

Overall, the diurnal cycle of OOA (Fig. 4i) was flat during both clean and polluted periods in this study, reflecting regional characteristics. This is also consistent with the wide distribution of OOA seen in the corresponding bivariate polar plot (Fig. 5j). Secondary inorganic species of nitrate and sulfate had different diurnal profiles. Nitrate had a pronounced diurnal cycle with much higher concentrations at night than during the day. As shown in Fig. 4d, the nitrate concentration decreased from 6.0 to 2.6 µg m−3 between 11:00 and 19:00 LT, which was mainly due to the evaporative loss of particulate ammonium nitrate particles due to high temperatures. Such diurnal cycles have been observed many times during the summer in megacities, e.g., Beijing (Sun et al., 2012), Nanjing (Ge et al., 2017), and Lanzhou (Xu et al., 2016), as well as those on other continents (e.g., Lanz et al., 2007). Nitrate also showed an increase in the early morning after sunrise (Fig. 4d). While the transport from urban sites played a role, this increase was mainly caused by daytime photochemical production when T was not high enough to substantially affect gas partitioning. During daytime clean periods, the diurnal cycle of nitrate was flat, with higher concentrations at night. Sulfate showed a much smoother diurnal evolution compared with nitrate (Fig. 4c), reflecting the regional characteristics of sulfate. Two small peaks were observed during the day. While the first peak between 09:00 and 11:00 LT was most likely from urban transport, consistent with regional sources from the southeast indicated by the corresponding bivariate polar plot (Fig. 5c), the second one was more likely from daytime photochemical production. Chloride accounted for a small fraction of the PM1 mass, yet it showed a strong diurnal cycle with a pronounced peak in the morning (Fig. 4e). This peak was noticeably similar to those of CO and SO2 (Fig. S11a and c), suggesting that the dominant source was from combustion emissions in the southeast.

## 3.3 Particle number size distributions

Figure 6 shows the time series of total number concentration (15–685 nm, N15–685) and three different modes including the small Aitken mode (15–40 nm, N15–40), the large Aitken mode (40–100 nm, N40–100), and the accumulation mode (100–400 nm, N100–685), as well as particles in the range of 7–15 nm (N7–15) that were calculated from the difference between MCPC and SMPS measurements. The corrections for diffusion loss and multiple charges have been applied in SMPS data analysis. The average total number concentration was 11 200 ± 5800 cm−3, which is comparable to that measured in Shangdianzi (12 000 cm−3; Shen et al., 2011) and Yufa (10 200 cm−3; Peng et al., 2014) and slightly higher than that observed in Beijing (10 100 cm−3; Du et al., 2017). N15–685 showed a pronounced diurnal cycle with a clear increase during the day (Fig. 7e). A further analysis highlights that this increase was mainly driven by small and large Aitken mode particles, indicating the impacts of new particle formation and growth on the diurnal variations in particle number (Fig. S14). For example, the small Aitken-mode particles and ultrafine particles showed rapid daytime increases after sunrise by more than a factor of 5 during both polluted events and clean periods. Note that the number concentrations of N7–15 during clean periods was much higher than during polluted events (6600 cm−3 vs. 3300 cm−3, on average), which is comparable to the number concentration of small Aitken-mode particles. These results show that (1) new particle formation was much stronger on clean days than on polluted days, and (2) new particle formation also occurred on polluted days, with newly formed particles growing quickly due to the higher condensation sink (CS; 0.05 s−1). About 49 % of the total particle number was made up of large Aitken-mode particles. However, these particles accounted for a small fraction of the total volume concentration (5 %). We also note that the pronounced nighttime peak in N40–100 coincidently agrees with that of COA during clean periods, suggesting the influence of cooking emissions on large Aitken-mode particles. We calculated the particle number size distribution for two nights that experienced significant cooking emission events, i.e., 2 and 13 May. Narrow single-mode distributions peaking at ∼60 nm were seen (Fig. S15), supporting the influence of cooking emissions on large Aitken-mode particles. The diurnal cycle of N100–685 was relatively flat except for a small morning peak (Fig. 7d). This is also consistent with results from the corresponding bivariate polar plot (Fig. S16d). The number concentration of N100–685 during polluted events was more than a factor of 3–4 times that during clean periods, which contributes toward the major difference in particle number characteristics between polluted events and clean periods. One reason is the higher CS during polluted events, which facilitated the growth of particles.

Figure 7Diurnal cycles of particle number concentration for (a) N7–15 calculated from the differences between MCPC and SMPS measurements, (b) N15–40 (15–40 nm), (c) N40–100 (40–100 nm), (d) N100–685 (100–685  nm), and (e) all particles, N15–685 (15–685 nm). Overall mean cycles are shown as black lines. Mean cycles for polluted events (PEs) and clear periods (CPs) are shown as red and blues lines, respectively.

Figure 8 shows that the average particle number concentration had a bimodal size distribution with the geometric mean diameter (GMD) peaking at 46 and 106 nm. Note that the peak diameter for the entire study was ∼62 nm, which is higher than that observed in urban Beijing (45 nm; Du et al., 2017). The higher CS (0.36 s−1 vs. 0.29 s−1) could be one of the reasons leading to more water vapor condensing on preexisting particles at the rural site. This is also consistent with the shift in the peak diameter from 57 nm during clean periods to 88 nm during polluted events. Although similar bimodal size distribution modes peaking at ∼46 nm (41 nm) and ∼109 nm (106 nm) were observed during clean periods and polluted events, respectively, the relative contributions of the two modes were largely different. While the particle number concentration was dominated by the small mode (56 %) during clean periods, the large mode was more important during polluted events, accounting for 73 % of the particle number concentration (Table 1). These results confirm the different roles of different mode particles between clean periods and polluted events.

Table 1Comparison of aerosol properties and meteorological parameters between polluted events and clean periods.

Figure 8Average particle number size distributions during (a) the entire study, (b) polluted events, and (c) clean periods.

Figure 9Average diurnal evolution of particle number size distributions and aerosol composition for new particle growth events during (a) the entire study, (b) polluted events, and (c) clean periods. The black solid lines in the top three panels show the diurnal cycles of CS. The circles and squares show the GMD from the log-normal fitting from this study and in Beijing (Du et al., 2017), respectively. The average diurnal cycles of aerosol species and CO-normalized aerosol species during polluted events and clean periods are shown in (b) and (c).

## 3.4 Particle growth events

New particle growth events (NPEs) were frequently observed during the study period. As shown in Fig. 9a, particle growth typically started at  08:00 LT and ended at midnight with an increase in GMD from ∼25 to ∼60 nm. This result is consistent with those previously reported for rural sites in the NCP (Wang et al., 2013) and urban sites, e.g., Beijing (Du et al., 2017). Note that the growth sizes of particles were overall larger than those observed in urban Beijing (from ∼22 to ∼55 nm), likely indicating a stronger aging process at the suburban site. The growth of particles tracked the diurnal cycle of CS, which showed a continuous increase from early morning to midnight. Although ACSM only has a good transmission for particles within the size range between 70 and 500 nm in diameter (Jayne et al., 2000), some particles from the small Aitken mode might not be detected by the ACSM. Simultaneous comparison between aerosol chemical composition and particle size distribution during the growth period make some sense. Aerosol composition seemed to significantly change during the growth period. As shown in Fig. 9a, OOA and sulfate were the only two species whose contributions increased, going from 26 % to 33 % and 27 % to 33 %, respectively, during the growth period (10:00–18:00 LT). Although the increases in sulfate and OOA were partly due to the decreases in nitrate and chloride because of the evaporative loss in the afternoon, the CO-normalized sulfate and OOA (SO4∕CO and OOA∕CO) also showed increases. Note that CO here was subtracted by a background value of 0.068 ppm that was calculated as the average of the lowest 5 % data in this study. These results show that both sulfate and OOA played important roles in daytime particle growth.

Figure 10(a) Particle growth rates (GRs) and the corresponding aerosol composition during the growth period, (b) the relationship between GR and CS, color-coded by OOA/PM1, and (c) the relationship between GR and sulfate concentration, color-coded by SO4∕PM1. The numbers over the circles in (b) represent different source regions (Fig. 11) and the triangles represent data without information about composition.

We also examined particle growth events on polluted days and clean days. As shown in Fig. 9b and c, the growth process of particles on polluted days started at  12:00 LT with the GMD increasing from ∼38 to ∼57 nm in 6 h. Particle growth started earlier ( 10:00 LT) on clean days with the GMD increasing from ∼27 to ∼41 nm. Particle growth on polluted days was faster than on clean days, which was likely due to the higher CS on polluted days. This is also consistent with more significant increases in SOA∕CO and sulfate∕CO on polluted days. The increase in CS between 09:00 and 12:00 LT on polluted days was associated with the corresponding increases in most aerosol species. However, the CO-normalized aerosol species did not show such an increase during this period suggesting that aerosol species and gaseous species were from the same air mass, i.e., transported in from urban sites. Such a pattern was not observed on clean days. Particle growth at night was more clearly seen on clean days than on polluted days, consistent with the increase in OOA∕CO. However, such an increase was not observed for sulfate∕CO on clean days, suggesting that OOA played a more important role than sulfate in particle growth at night.

Figure 11Particle GR and sulfate concentrations during the days shown in Fig. 10a with different footprints: (a) north, (b) southwest, and (c) east-southeast.

Figure 12Average composition of aerosol particles and particle number concentration for three different clusters.

We further calculated the particle growth rate (GR) of each growth event that lasted more than 3 h (Fig. 10a and Table S1). The particle growth rates were calculated using Eq. (1).

$\begin{array}{}\text{(1)}& \text{GR}=\mathrm{\Delta }{D}_{\text{m}}/\mathrm{\Delta }t,\end{array}$

where Dm is the geometric mean diameter from the log-normal fitting of each size distribution and ΔDm is the increase in diameter during the growth period of Δt. The particle GR varied from 1.2 to 4.7 nm h−1, which generally falls within the range of values (1–20 nm h−1) based on observations from around the world (Kulmala et al., 2004; Yu et al., 2017). As indicated in Fig. 10b, GR was positively correlated with CS for most of the time, mainly for the fraction of OOA higher than 30 % (r2=0.61, Fig. S17a). As CS increased from ∼0.01 to 0.05 s−1, GR increased from 3 to 5 nm h−1. Higher CS is usually associated with lower particle GR due to faster consumption of condensable vapors. The positive relationship between GR and CS might indicate that the source of condensable vapors contributes to the increase in CS prior to the observation. Another possible explanation is that heterogeneous surface chemistry was more important than CS in growing aerosol particles in highly polluted environment (Kulmala et al., 2017). Figure 10b also showed some low GRs with high CS, which were characterized by low contributions of OOA. These results support the importance of the involvement of OOA in particle growth (Wu et al., 2017). Figure 10c further shows that GR was positively correlated with the concentration of sulfate during periods with low sulfate mass loadings (<3µg m−3, r2=0.42, Fig. S17b), while periods with higher concentrations of sulfate had lower GRs. This is corroborated by Fig. S18, showing that GR decreased as the sulfate contribution increased. By contrast, GR was positively correlated with OOAPM1 for most of the time. These results highlight that OOA played a more important role in particle growth than sulfate, although sulfate was also important during periods with low mass loadings. We further checked the dependence of GR on source region, but no clear relationship was found. For example, GR varied from 1.6 to 3.2 nm h−1 during periods with air masses from the north, 1.8–4.9 nm h−1 from the southwest, and 1.8–4.6 nm h−1 from the east-southeast (Fig. 11). This shows that GR had no clear dependence on air masses from different regions, although sulfate concentrations showed much difference.

## 3.5 Aerosol composition and particle number concentrations from different source regions

Figure 12 shows the average composition and particle number distributions from different source regions. The air masses in cluster 3 (C3, 39 % of the time) were mainly from the southeast while the other two clusters were mainly from the northwest. The average PM1 concentration for C3 was 42 µg m−3, which is ∼30 % higher than that from the other two source regions. The high mass loadings for this cluster suggest that source regions to the southeast were responsible for the high PM pollution at the sampling site. This is also supported by the bivariate polar plots of aerosol species with high concentration regions in the southeast (Fig. 5). Among all clusters, the two dominant species in PM1 were SOA (26 %–28 %) and sulfate (26 %–29 %). Overall, the aerosol bulk composition for the three clusters from different source areas was similar. These results suggest that aerosol particles were relatively well mixed over the region around Xingtai. Particle number concentrations showed more differences among the different clusters. The two clusters from the northwest were both dominated by large Aitken-mode particles, on average accounting for 49 % and 51 % of the total particles, respectively. Although large Aitken-mode particles dominated the total particle number for C3, we also observed large increases in accumulation-mode particles (39 %) compared with the other two clusters (28 %–30 %). These results suggest that air masses from the northwest were relatively clean, which led to a more frequent occurrence of new particle formation and growth events, while those from the southeast with higher PM loadings tended to form more large particles due to the high CS.

4 Conclusions

We presented an analysis of aerosol chemistry and particle growth events at an urban downwind site in the NCP during May and June of 2016 using real-time measurements from an ACSM, an SMPS, and a suite of collocated instruments. Our results showed that the PM1 level in spring and summer (30.5 µg m−3) was much lower than that during wintertime when coal combustion emissions were enhanced. Similar to previous studies with a focus on the NCP, aerosol composition at the downwind site of Xingtai was dominated by organics (38 %), 78 % of which was identified as secondary OA according to the ME-2 analysis. Secondary aerosols (i.e., SNA + SOA) accounted for 78 % of PM1, highlighting the major source of secondary formation and regional transport at the downwind site. Local sources of PM from traffic and cooking emissions accounted for less than 10 % of the total PM1 mass. We also observed daytime transport from urban sites leading to similar pronounced early morning peaks for most aerosol species. New particle growth events were frequently observed (58 % of the time) during the study period. By linking the GR with aerosol composition, we found that OOA and sulfate were two major species affecting the growth of particles. In particular, both OOA and sulfate played important roles in particle growth during clean periods, while OOA was more important than sulfate during polluted events. This is also supported by the decrease in GR as the sulfate contribution increased. A further analysis showed that particle growth rates have no clear dependence on air mass trajectories.

Data availability
Data availability.

The data in this study are available from the authors upon request (sunyele@mail.iap.ac.cn).

Supplement
Supplement.

Author contributions
Author contributions.

YS and ZL designed research. YZ, WD, YW, QW, HW, HZ, FZ, HS, YB, YH, PF, TZ, PW, and YS performed research. FC and AP shared the data analysis software. YZ, WD, and YS analyzed the data and wrote the paper. YZ and WD have contributed equally to this work.

Competing interests
Competing interests.

The authors declare that they have no conflict of interest.

Special issue statement
Special issue statement.

Acknowledgements
Acknowledgements.

This study was supported by the National Key Project of Basic Research (2013CB955801, 2014CB447900) and the National Natural Science Foundation of China (41575120).

Edited by: Hang Su
Reviewed by: two anonymous referees

References

Boucher, O., Randall, D., Artaxo, P., Bretherton, C., Feingold, G., Forster, P., Kerminen, V., Kondo, Y., Liao, H., and Lohmann, U.: Clouds and Aerosols: in: Climate Change 2013: The Physical Science Basis Contribution of Working Group I to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change, edited by: Tignor, M., Allen, S. K., Boschung, J., Nauels, A., Xia, Y., Bex, V., Midgley, P. M., Cambridge University Press, Cambridge, UK, New York, NY, USA, 2013.

Canonaco, F., Crippa, M., Slowik, J. G., Baltensperger, U., and Prévôt, A. S. H.: SoFi, an IGOR-based interface for the efficient use of the generalized multilinear engine (ME-2) for the source apportionment: ME-2 application to aerosol mass spectrometer data, Atmos. Meas. Tech., 6, 3649–3661, https://doi.org/10.5194/amt-6-3649-2013, 2013.

Chen, Y., Ebenstein, A., Greenstone, M., and Li, H.: From the cover: evidence on the impact of sustained exposure to air pollution on life expectancy from China's Huai River policy, P. Natl. Acad. Sci. USA, 110, 12936–12941, 2013.

Crippa, M., DeCarlo, P. F., Slowik, J. G., Mohr, C., Heringa, M. F., Chirico, R., Poulain, L., Freutel, F., Sciare, J., Cozic, J., Di Marco, C. F., Elsasser, M., Nicolas, J. B., Marchand, N., Abidi, E., Wiedensohler, A., Drewnick, F., Schneider, J., Borrmann, S., Nemitz, E., Zimmermann, R., Jaffrezo, J.-L., Prévôt, A. S. H., and Baltensperger, U.: Wintertime aerosol chemical composition and source apportionment of the organic fraction in the metropolitan area of Paris, Atmos. Chem. Phys., 13, 961–981, https://doi.org/10.5194/acp-13-961-2013, 2013.

Crippa, M., Canonaco, F., Lanz, V. A., Äijälä, M., Allan, J. D., Carbone, S., Capes, G., Ceburnis, D., Dall'Osto, M., Day, D. A., DeCarlo, P. F., Ehn, M., Eriksson, A., Freney, E., Hildebrandt Ruiz, L., Hillamo, R., Jimenez, J. L., Junninen, H., Kiendler-Scharr, A., Kortelainen, A.-M., Kulmala, M., Laaksonen, A., Mensah, A. A., Mohr, C., Nemitz, E., O'Dowd, C., Ovadnevaite, J., Pandis, S. N., Petäjä, T., Poulain, L., Saarikoski, S., Sellegri, K., Swietlicki, E., Tiitta, P., Worsnop, D. R., Baltensperger, U., and Prévôt, A. S. H.: Organic aerosol components derived from 25 AMS data sets across Europe using a consistent ME-2 based source apportionment approach, Atmos. Chem. Phys., 14, 6159–6176, https://doi.org/10.5194/acp-14-6159-2014, 2014.

Drinovec, L., Mocnik, G., Zotter, P., Prévôt, A. S. H., Ruckstuhl, C., Coz, E., Rupakheti, M., Sciare, J., Müller, T., Wiedensohler, A., and Hansen, A. D. A.: The “dual-spot” Aethalometer: an improved measurement of aerosol black carbon with real-time loading compensation, Atmos. Meas. Tech., 8, 1965–1979, https://doi.org/10.5194/amt-8-1965-2015, 2015.

Du, W., Zhao, J., Wang, Y., Zhang, Y., Wang, Q., Xu, W., Chen, C., Han, T., Zhang, F., Li, Z., Fu, P., Li, J., Wang, Z., and Sun, Y.: Simultaneous measurements of particle number size distributions at ground level and 260 m on a meteorological tower in urban Beijing, China, Atmos. Chem. Phys., 17, 6797–6811, https://doi.org/10.5194/acp-17-6797-2017, 2017.

Elser, M., Huang, R.-J., Wolf, R., Slowik, J. G., Wang, Q., Canonaco, F., Li, G., Bozzetti, C., Daellenbach, K. R., Huang, Y., Zhang, R., Li, Z., Cao, J., Baltensperger, U., El-Haddad, I., and Prévôt, A. S. H.: New insights into PM2.5 chemical composition and sources in two major cities in China during extreme haze events using aerosol mass spectrometry, Atmos. Chem. Phys., 16, 3207–3225, https://doi.org/10.5194/acp-16-3207-2016, 2016.

Fu, G. Q., Xu, W. Y., Yang, R. F., Li, J. B., and Zhao, C. S.: The distribution and trends of fog and haze in the North China Plain over the past 30 years, Atmos. Chem. Phys., 14, 11949–11958, https://doi.org/10.5194/acp-14-11949-2014, 2014.

Fu, H. and Chen, J.: Formation, features and controlling strategies of severe haze-fog pollutions in China, Sci. Total Environ., 578, 121–138, https://doi.org/10.1016/j.scitotenv.2016.10.201, 2017.

Ge, X., He, Y., Sun, Y., Xu, J., Wang, J., Shen, Y., and Chen, M.: Characteristics and formation mechanisms of fine particulate nitrate in typical urban areas in China, Atmosphere, 8, 62, https://doi.org/10.3390/atmos8030062, 2017.

Guo, S., Hu, M., Zamora, M. L., Peng, J., Shang, D., Zheng, J., Du, Z., Wu, Z., Shao, M., Zeng, L., Molina, M. J., and Zhang, R.: Elucidating severe urban haze formation in China, P. Natl. Acad. Sci. USA, 111, 17373–17378, https://doi.org/10.1073/pnas.1419604111, 2014.

Hu, W., Hu, M., Hu, W., Jimenez, J. L., Yuan, B., Chen, W., Wang, M., Wu, Y., Chen, C., Wang, Z., Peng, J., Zeng, L., and Shao, M.: Chemical composition, sources, and aging process of submicron aerosols in Beijing: contrast between summer and winter, J. Geophys. Res.-Atmos., 121, 1955–1977, https://doi.org/10.1002/2015JD024020, 2016.

Jayne, J. T., Leard, D. C., Zhang, X., Davidovits, P., Smith, K. A., Kolb, C. E., and Worsnop, D. R.: Development of an aerosol mass spectrometer for size and composition analysis of submicron particles, Aerosol Sci. Tech., 33, 49–70, 2000.

Kulmala, M., Vehkamäki, H., Petäjä, T., Dal Maso, M., Lauri, A., Kerminen, V. M., Birmili, W., and McMurry, P. H.: Formation and growth rates of ultrafine atmospheric particles: a review of observations, J. Aerosol Sci., 35, 143–176, https://doi.org/10.1016/j.jaerosci.2003.10.003, 2004.

Kulmala, M., Kerminen, V. M., Petaja, T., Ding, A. J., and Wang, L.: Atmospheric gas-to-particle conversion: why NPF events are observed in megacities?, Faraday Discuss., 200, 271–288, https://doi.org/10.1039/C6FD00257A, 2017.

Lanz, V. A., Alfarra, M. R., Baltensperger, U., Buchmann, B., Hueglin, C., and Prévôt, A. S. H.: Source apportionment of submicron organic aerosols at an urban site by factor analytical modelling of aerosol mass spectra, Atmos. Chem. Phys., 7, 1503–1522, https://doi.org/10.5194/acp-7-1503-2007, 2007.

Li, H., Zhang, Q., Zhang, Q., Chen, C., Wang, L., Wei, Z., Zhou, S., Parworth, C., Zheng, B., Canonaco, F., Prévôt, A. S. H., Chen, P., Zhang, H., Wallington, T. J., and He, K.: Wintertime aerosol chemistry and haze evolution in an extremely polluted city of the North China Plain: significant contribution from coal and biomass combustion, Atmos. Chem. Phys., 17, 4751–4768, https://doi.org/10.5194/acp-17-4751-2017, 2017.

Li, Y. J., Sun, Y., Zhang, Q., Li, X., Li, M., Zhou, Z., and Chan, C. K.: Real-time chemical characterization of atmospheric particulate matter in China: A review, Atmos. Environ., 158, 270–304, 2017.

Massoli, P., Kebabian, P. L., Onasch, T. B., Hills, F. B., and Freedman, A.: Aerosol light extinction measurements by cavity attenuated phase shift (CAPS) spectroscopy: laboratory validation and field deployment of a compact aerosol particle extinction monitor, Aerosol Sci. Tech., 44, 428–435, 2010.

Matthew, B. M., Middlebrook, A. M., and Onasch, T. B.: Collection efficiencies in an aerodyne aerosol mass spectrometer as a function of particle phase for laboratory generated aerosols, Aerosol Sci. Tech., 42, 884–898, 2008.

Middlebrook, A. M., Bahreini, R., Jimenez, J. L., and Canagaratna, M. R.: Evaluation of composition-dependent collection efficiencies for the aerodyne aerosol mass spectrometer using field data, Aerosol Sci. Tech., 46, 258–271, 2012.

Mohr, C., DeCarlo, P. F., Heringa, M. F., Chirico, R., Slowik, J. G., Richter, R., Reche, C., Alastuey, A., Querol, X., Seco, R., Peñuelas, J., Jiménez, J. L., Crippa, M., Zimmermann, R., Baltensperger, U., and Prévôt, A. S. H.: Identification and quantification of organic aerosol from cooking and other sources in Barcelona using aerosol mass spectrometer data, Atmos. Chem. Phys., 12, 1649–1665, https://doi.org/10.5194/acp-12-1649-2012, 2012.

Ng, N. L., Canagaratna, M. R., Jimenez, J. L., Zhang, Q., Ulbrich, I. M., and Worsnop, D. R.: Real-time methods for estimating organic component mass concentrations from aerosol mass spectrometer data, Environ. Sci. Tech., 45, 910–916, https://doi.org/10.1021/es102951k, 2011.

Peng, J. F., Hu, M., Wang, Z. B., Huang, X. F., Kumar, P., Wu, Z. J., Guo, S., Yue, D. L., Shang, D. J., Zheng, Z., and He, L. Y.: Submicron aerosols at thirteen diversified sites in China: size distribution, new particle formation and corresponding contribution to cloud condensation nuclei production, Atmos. Chem. Phys., 14, 10249–10265, https://doi.org/10.5194/acp-14-10249-2014, 2014.

Petit, J. E., Favez, O., Albinet, A., and Canonaco, F.: A user-friendly tool for comprehensive evaluation of the geographical origins of atmospheric pollution: wind and trajectory analyses, Environ. Modell. Softw., 88, 183–187, https://doi.org/10.1016/j.envsoft.2016.11.022, 2017.

Quan, J., Zhang, Q., He, H., Liu, J., Huang, M., and Jin, H.: Analysis of the formation of fog and haze in North China Plain (NCP), Atmos. Chem. Phys., 11, 8205–8214, https://doi.org/10.5194/acp-11-8205-2011, 2011.

Shen, X. J., Sun, J. Y., Zhang, Y. M., Wehner, B., Nowak, A., Tuch, T., Zhang, X. C., Wang, T. T., Zhou, H. G., Zhang, X. L., Dong, F., Birmili, W., and Wiedensohler, A.: First long-term study of particle number size distributions and new particle formation events of regional aerosol in the North China Plain, Atmos. Chem. Phys., 11, 1565–1580, https://doi.org/10.5194/acp-11-1565-2011, 2011.

Stein, A. F., Draxler, R. R., Rolph, G. D., Stunder, B. J. B., Cohen, M. D., and Ngan, F.: NOAA's HYSPLIT atmospheric transport and dispersion modeling system, B. Am. Meteorol. Soc., 96, 2059–2077, https://doi.org/10.1175/bams-d-14-00110.1, 2015.

Stohl, A., Forster, C., Frank, A., Seibert, P., and Wotawa, G.: Technical note: The Lagrangian particle dispersion model FLEXPART version 6.2, Atmos. Chem. Phys., 5, 2461–2474, https://doi.org/10.5194/acp-5-2461-2005, 2005.

Sun, Y., Wang, Z., Dong, H., Yang, T., Li, J., Pan, X., Chen, P., and Jayne, J. T.: Characterization of summer organic and inorganic aerosols in Beijing, China with an aerosol chemical speciation monitor, Atmos. Environ., 51, 250–259, https://doi.org/10.1016/j.atmosenv.2012.01.013, 2012.

Sun, Y., Wang, Z. F., Fu, P. Q., Yang, T., Jiang, Q., Dong, H. B., Li, J., and Jia, J. J.: Aerosol composition, sources and processes during wintertime in Beijing, China, Atmos. Chem. Phys., 13, 4577–4592, https://doi.org/10.5194/acp-13-4577-2013, 2013.

Sun, Y., Jiang, Q., Wang, Z., Fu, P., Li, J., Yang, T., and Yin, Y.: Investigation of the sources and evolution processes of severe haze pollution in Beijing in January 2013, J. Geophys. Res.-Atmos., 119, 4380–4398, https://doi.org/10.1002/2014jd021641, 2014.

Sun, Y., Jiang, Q., Xu, Y., Ma, Y., Zhang, Y., Liu, X., Li, W., Wang, F., Li, J., Wang, P., and Li, Z.: Aerosol characterization over the North China Plain: haze life cycle and biomass burning impacts in summer, J. Geophys. Res.-Atmos., 121, 2508–2521, https://doi.org/10.1002/2015JD024261, 2016a.

Sun, Y., Chen, C., Zhang, Y. J., Xu, W. Q., Zhou, L. B., Cheng, X. L., Zheng, H. T., Ji, D. S., Li, J., Tang, X., Fu, P. Q., and Wang, Z. F.: Rapid formation and evolution of an extreme haze episode in Northern China during winter 2015, Sci. Rep.-UK, 6, 27151, https://doi.org/10.1038/srep27151, 2016b.

Ulbrich, I. M., Canagaratna, M. R., Zhang, Q., Worsnop, D. R., and Jimenez, J. L.: Interpretation of organic components from Positive Matrix Factorization of aerosol mass spectrometric data, Atmos. Chem. Phys., 9, 2891–2918, https://doi.org/10.5194/acp-9-2891-2009, 2009.

Wang, L., Wei, Z., Wei, W., Fu, J. S., Meng, C., and Ma, S.: Source apportionment of PM2.5 in top polluted cities in Hebei, China using the CMAQ model, Atmos. Environ., 122, 723–736, https://doi.org/10.1016/j.atmosenv.2015.10.041, 2015.

Wang, Q., Sun, Y., Jiang, Q., Du, W., Sun, C., Fu, P., and Wang, Z.: Chemical composition of aerosol particles and light extinction apportionment before and during the heating season in Beijing, China, J. Geophys. Res.-Atmos., 120, 12708–12722, https://doi.org/10.1002/2015JD023871, 2015.

Wang, Q., Zhao, J., Du, W., Ana, G., Wang, Z., Sun, L., Wang, Y., Zhang, F., Li, Z., Ye, X., and Sun, Y.: Characterization of submicron aerosols at a suburban site in central China, Atmos. Environ., 131, 115–123, https://doi.org/10.1016/j.atmosenv.2016.01.054, 2016.

Wang, Z. B., Hu, M., Sun, J. Y., Wu, Z. J., Yue, D. L., Shen, X. J., Zhang, Y. M., Pei, X. Y., Cheng, Y. F., and Wiedensohler, A.: Characteristics of regional new particle formation in urban and regional background environments in the North China Plain, Atmos. Chem. Phys., 13, 12495–12506, https://doi.org/10.5194/acp-13-12495-2013, 2013.

Wu, Z.: Gain Insight into Chemical Components Driving New Particle Growth on a Basis of Particle Hygroscopicity and Volatility Measurements: a Short Review, Curr. Pollut. Rep., 3, 175–181, https://doi.org/10.1007/s40726-017-0064-6, 2017.

Xu, J., Shi, J., Zhang, Q., Ge, X., Canonaco, F., Prévôt, A. S. H., Vonwiller, M., Szidat, S., Ge, J., Ma, J., An, Y., Kang, S., and Qin, D.: Wintertime organic and inorganic aerosols in Lanzhou, China: sources, processes, and comparison with the results during summer, Atmos. Chem. Phys., 16, 14937–14957, https://doi.org/10.5194/acp-16-14937-2016, 2016.

Yang, Y. R., Liu, X. G., Qu, Y., An, J. L., Jiang, R., Zhang, Y. H., Sun, Y. L., Wu, Z. J., Zhang, F., Xu, W. Q., and Ma, Q. X.: Characteristics and formation mechanism of continuous hazes in China: a case study during the autumn of 2014 in the North China Plain, Atmos. Chem. Phys., 15, 8165–8178, https://doi.org/10.5194/acp-15-8165-2015, 2015.

Yu, H., Ren, L., and Kanawade, V. P.: New Particle Formation and Growth Mechanisms in Highly Polluted Environments, Curr. Pollut. Rep., 3, 245–253, 10.1007/s40726-017-0067-3, 2017.

Zhang, Q., Jimenez, J. L., Worsnop, D. R., and Canagaratna, M.: A case study of urban particle acidity and its influence on secondary organic aerosol, Environ. Sci. Tech., 41, 3213–3219, https://doi.org/10.1021/es061812j, 2007.

Zhang, Q., Jimenez, J. L., Canagaratna, M. R., Ulbrich, I. M., Ng, N. L., Worsnop, D. R., and Sun, Y.: Understanding atmospheric organic aerosols via factor analysis of aerosol mass spectrometry: a review, Anal. Bioanal. Chem., 401, 3045–3067, https://doi.org/10.1007/s00216-011-5355-y, 2011.

Zhang, Q. H., Zhang, J. P., and Xue, H. W.: The challenge of improving visibility in Beijing, Atmos. Chem. Phys., 10, 7821–7827, https://doi.org/10.5194/acp-10-7821-2010, 2010.

Zhang, R., Wang, G., Guo, S., Zarnora, M. L., Ying, Q., Lin, Y., Wang, W., Hu, M., and Wang, Y.: Formation of urban fine particulate matter, Chem. Rev., 115, 3803–3855, https://doi.org/10.1021/acs.chemrev.5b00067, 2015.

Zhang, Y., Sun, Y., Du, W., Wang, Q., Chen, C., Han, T., Lin, J., Zhao, J., Xu, W., Gao, J., Li, J., Fu, P., Wang, Z., and Han, Y.: Response of aerosol composition to different emission scenarios in Beijing, China, Sci. Total Environ., 571, 902–908, https://doi.org/10.1016/j.scitotenv.2016.07.073, 2016.

Zhang, Y. J., Tang, L. L., Wang, Z., Yu, H. X., Sun, Y. L., Liu, D., Qin, W., Canonaco, F., Prévôt, A. S. H., Zhang, H. L., and Zhou, H. C.: Insights into characteristics, sources, and evolution of submicron aerosols during harvest seasons in the Yangtze River delta region, China, Atmos. Chem. Phys., 15, 1331-1349, https://doi.org/10.5194/acp-15-1331-2015, 2015.

Zhao, X. J., Zhao, P. S., Xu, J., Meng,, W., Pu, W. W., Dong, F., He, D., and Shi, Q. F.: Analysis of a winter regional haze event and its formation mechanism in the North China Plain, Atmos. Chem. Phys., 13, 5685–5696, https://doi.org/10.5194/acp-13-5685-2013, 2013.

Zhou, Y., Cheng, S., Lang, J., Chen, D., Zhao, B., Liu, C., Xu, R., and Li, T.: A comprehensive ammonia emission inventory with high-resolution and its evaluation in the Beijing-Tianjin-Hebei (BTH) region, China, Atmos. Environ., 106, 305–317, https://doi.org/10.1016/j.atmosenv.2015.01.069, 2015.