1,475
Views
17
CrossRef citations to date
0
Altmetric
Research Article

Characteristics of chemical composition and role of meteorological factors during heavy aerosol pollution episodes in northern Beijing area in autumn and winter of 2015

, , , , , , , & show all
Article: 1347484 | Received 21 Jan 2017, Accepted 20 Jun 2017, Published online: 26 Jul 2017

Abstract

Heavy aerosol pollution episodes (HPEs) usually start from late autumn and become more serious in winter in Beijing and its vicinity (BIV). In this study, we examine the reasons for the formation and change of HPEs in the areas of northern BIV. The size-resolved chemical components of PM1 and meteorological conditions were investigated during HPEs in autumn and winter of 2015. Stable regional atmosphere and higher atmospheric condensation processes associated with southerly and lower speed wind led to the formation of HPEs. After the start of these HPEs, the concentration of fine particles increased more than twice in several hours. ~80% of the ‘explosive’ growth in PM mass can be considered as a positive feedback of meteorological factors that come from even more stable atmosphere and larger condensation rate of water vapour, which was derived from the interaction between formed aerosol pollution and the meteorological factors within boundary layer. Nitrate was the largest fraction of PM1 in autumn, and the most significantly increased component during HPEs relative to clean period during both of autumn and winter. The proportion of organic aerosol (OA) was similar to that of nitrate in autumn, but its rise in HPE was much smaller, mainly because of the high concentration of OA existed during clean periods. Compared with the largest increase of nitrate, the similar uplift was found for ammonium production, showing that a large amount of ammonium was mainly formed by the combination of in HPEs, rather than . In addition to the lower southerly wind carrying pollutants from southern part of BIV, westerly wind from central Inner Mongolia and north Shanxi can also bring air pollutants originating from coal combustion, contributing to the heavy pollution in the northern BIV area in winter, and resulting in higher sulphate, nitrate and OA masses.

1. Introduction

In recent decades, Chinese aerosol pollution has become increasingly serious. This is partially due to large emissions of air pollutants from anthropogenic activities, e.g. traffic, industry, and power plants (Zhang et al., Citation2008; Sun et al., Citation2013), but is also influenced by increased bad weather conditions (Choobari et al., Citation2012; Zhang et al., Citation2015a). In particular, heavy aerosol pollution episodes (HPEs) have often occurred and formed over a short period since January of 2013. During typical HPEs in the Beijing and its vicinity (BIV), mass concentrations of PM2.5 can double in several hours, which may be resulted from previously un-described chemical mechanisms, but the interaction between aerosol pollutants and the conditions of boundary layer may also play an important role. Some studies have found that if an atmospheric pollution model does not take into account the impact of atmospheric pollution on the stability of the low layer of atmosphere and their induced feedback effect, the model will seriously underestimate the heavy pollution process (Wang et al., Citation2014b; Zheng et al., Citation2015). Some studies have reported that model simulation can lead to a difference of monthly PM2.5 mean concentration of 10–30% if we do not consider the two-way interaction (Wang et al., Citation2014a). This difference in the hourly maximum value of PM2.5 in the North China Plain can reach to 30–70% (Gao et al., Citation2015). Numerical simulations have shown that the turbulent diffusion led by the interaction between pollutants and the boundary layer decreased by about 52% in a HPE in 2008 and the mixed layer height decreased by about 33%. However, wind speed within the boundary layer increased by about 13%, and the positive feedback from the boundary layer changes correspondingly increased the surface mean concentration of PM2.5 during the HPE by about 13%, but there was no effect on the beginning and end of the HPE (Wang et al., Citation2015a, Citation2015b).

The formation and change in secondary aerosol particles are highly affected by the meteorological conditions during HPEs (Zhang et al., Citation2015a). Investigation of the source and characterization of fine aerosol particle can help us to understand heavy aerosol pollution.

The chemical complexity and labile nature of atmospheric PM strongly favours real-time instrumental analysis techniques that characterize pertinent physical and chemical properties without having to collect, store, and transport samples (Canagaratna et al., Citation2007). High resolution time of flight aerosol mass spectrometry (HR-TOF-AMS) has been widely used in recent years because it has high temporal and spectral resolution, and can quickly and accurately measure the particle mass concentration, chemical composition and particle size distribution. To date, regional research using AMS measurements on fine aerosol have mainly focused on Europe (Schneider et al., Citation2013), North America (Aiken et al., Citation2009) and East Asia (Takegawa et al., Citation2006; Zhang et al., Citation2014). Zhang summarized the AMS observations of 37 sites in the northern hemisphere mid-latitude region, including urban areas, rural areas and background stations in remote areas, in terms of the characteristics of chemical composition, mass concentration and size distribution in these different areas (Zhang et al., Citation2007). In China, AMS-based field campaigns and studies have mainly focused on urban areas such as Beijing (Zhang et al., Citation2013b), Shanghai (Tang et al., Citation2014) and Guangzhou (Huang et al., Citation2011). Zhang et al. conducted an AMS observation experiment in the Beijing urban area in 2008 and examined the characteristics of chemical composition and seasonal changes of PM1 (Zhang et al., Citation2013b). However, knowledge about sources, and changing processes of background aerosol in other areas, especially during HPEs in North China (Quan et al., Citation2014), remains limited. Therefore, the regional HPEs have become the focus of intensive research activities.

In this study, we examined the chemical composition of background NR-PM1 (non-refractory PM1) in the BIV area during HPEs in autumn and winter of 2015. In addition, we analysed the parameterized index linking air pollution with meteorological conditions (PLAM) (Zhang et al., Citation2009) and other meteorological elements to evaluate specific meteorological conditions and their feedback on HPEs. These approaches will provide insights into the reasons behind aggravated heavy pollution processes in the first few hours of HPEs in BIV.

2. Method

2.1. Site and experiment descriptions

Shangdianzi (SDZ), one of the regional stations of Global Atmospheric Watch (GAW) programme, located at Shangdianzi village, Gaoling Town of Miyun County (117.07°E, 40.39°N, 287 m a.s.l). It is surrounded by hilly vegetation with limited local anthropogenic pollutant sources. On the territory, it belongs to the northern part of North China Plain, about 55 and 150 km north-east of the urban areas of Miyun and Beijing, respectively. The observations at SDZ are representative for atmospheric conditions of the BIV (Xu, Citation1992; An et al., Citation2012) (Figs. and ).

Fig. 1. Location of Shangdianzi (SDZ) regional GAW station, and back trajectories per hour of air mass in SDZ during autumn.

Fig. 1. Location of Shangdianzi (SDZ) regional GAW station, and back trajectories per hour of air mass in SDZ during autumn.

Fig. 2. Location of Shangdianzi (SDZ) regional GAW station, and back trajectories per hour of air mass in SDZ during winter.

Fig. 2. Location of Shangdianzi (SDZ) regional GAW station, and back trajectories per hour of air mass in SDZ during winter.

This field campaign conducted at two periods, one is from 31 October to 6 November 2015, another is from 16 December to 27 December 2015, which covered autumn and wintertime respectively. A HR-TOF-AMS (Aerodyne Research Inc. Boston, TX, USA) was utilized to measure the size-resolved mass concentration of chemical components (organics, sulphate, nitrate, ammonium and chloride) in non-refractory PM1 (NR-PM1) (Decarlo et al., Citation2006; Kimmel et al., Citation2011). A twin differential mobility particle sizer (TDMPS) were also used to measure aerosol particle number concentration (Shen et al., Citation2014). An aerosol particle mass analyzer (Malloy et al., Citation2009) was used to measure the density of the aerosol. In addition, hourly meteorological data including wind direction, wind speed, temperature, humidity, air pressure, precipitation as well as trace gas data (including NOx, SO2, CO, and O3) were measured at the same time. Except the automatic weather station, all the other instruments were set in the second-floor laboratory. Atmospheric air was introduced to the instruments through a PM10 inlet placed on the roof of the room, with a flow rate of 16.7 l/min. The sample air was dried by an auto-regenerated dryer system (Tuch et al., Citation2009), and the relative humidity inside the inlet was less than 30%.

2.2. Instrument operation and calibration

The HR-ToF-AMS has been described in detail previously (Decarlo et al., Citation2006) and improves upon previous versions of the AMS (Canagaratna et al., Citation2007) by using a custom high-resolution TOFMS (Tofwerk) (Aiken et al., Citation2008). The main advantage of the HR-ToF-AMS over previous versions of AMS is the much improved ability of identification and separation of isobaric ions (especially for m/z < 100) that have the same nominal mass but are slightly different in exact mass due to differences in elemental composition (Decarlo et al., Citation2006). As a result, the high resolution mass spectral data can provide valuable information on the elemental composition (e.g. C, H, O and N) and thus OM/OC ratio of organic aerosol (OA) (Aiken et al., Citation2008). The TOFMS has two operating modes, utilizing two ion optical configurations, V and W modes. Because the W mode has the highest mass resolution, it is used whenever the data are not signal-limited (Aiken et al., Citation2008).

During the campaign, the HR-ToF-AMS operated in a cycle of 3 modes every 4 min, including: 1 min V-mode to obtain the mass concentrations of the non-refractory species; 1 min W-mode to obtain high resolution mass spectral data; 1 min separate PToF (particle time-of-flight) mode to determine size distributions of species under the V-mode. Moreover, a series of calibrations for HR-TOF-AMS, including flow rate calibration, ionization efficiency calibration and particle size calibration were performed before the campaign (Allan et al., Citation2003; Jimenez et al., Citation2003). A default collection efficiency factor of 0.5 is recommended for the sulphate aerosols internally mixed with other inorganic species and oxidized organic species (Canagaratna et al., Citation2007; Middlebrook et al., Citation2012). In the cases of externally mixed aerosols, the particle size, composition, humidity and temperature of inlet might impact the CE values (Allan et al., Citation2004; Weimer et al., Citation2006). The EL (Dva) has been discussed as a factor impact the AMS quantitatively measures mass loadings. For particles falling within the vacuum aerodynamic size range for which the standard lens was designed (60–600 nm), EL(dva) is 100% (Jayne et al., Citation2000). For typical accumulation mode particles (∼50 to 1000 nm diameter), the EL is close to 100%. And the correlation between AMS data and other collocated instrumentation is within ±25%, even over long time periods and over a variety of chemical compositions (Takegawa et al., Citation2005; Canagaratna et al., Citation2007). In this case, it is reliable to determine the CE by comparing the AMS and other collocated instrumentation data. The TDMPS data have been used to verify the CE and correct AMS data by several previous works (Zhang et al., Citation2014, Citation2015b).

The TDMPS system includes two differential mobility analyzers and two condensation particle counters (TSI Model 3776 and 3772) was deployed to detect the size distributions of particles with mobility diameters ranging from 3 to 800 nm. The cutoff diameter for AMS is about 1um in vacuum aerodynamic diameter (dva), while the TDMPS system reports the mobility diameter (dm). Generally, the particle density can be considered as a link between dm and dva for spherical particles. Based on our collocated measurement of aerosol density, the average density in this campaign varied 1.5–1.7 g/cm3. In order to compare the same size range of particles to correct the CE, the number concentrations data of particles less than 600 nm in Dm (approximately equivalent of 1000 nm in Dva) were selected to calculate the mass concentration by the formula as follows:(1)

In the formula, ρ is the average aerosol density, D represents electric mobility size. ND is the number concentration of particle whose electric mobility size equal to D (Shen et al., Citation2014). The calculated PM1 mass concentration was fitting with the mass concentration measured by AMS. As a result, collection efficiency factor of 0.38 was employed to correct the mass concentration of AMS.

3. Results and discussion

3.1. Large increase of chemical components, fine aerosol particles and the role of meteorological factor during HPEs

In this paper, we are actually interested in the feedback effect of meteorological conditions on the ‘explosive’ growth phase of PM in a few hours to 10 some hours after the formation of a pollution event with PM mass concentration exceeding the level of clean period (CP). In order to explore these, five HPEs were selected during autumn (from 31 October 2015 to 6 November 2015) and winter (from 17 December 2015 to 27 December 2015), respectively, when one can find five ‘explosive’ growth phases of PM during the early stage of pollution formation.

3.1.1. Autumn

The mass concentrations with the change of particle size for organics, nitrate, sulphate, ammonium measured by HR-TOF-AMS during autumn, as well as number concentration of PM1 measured by TDMPS, and various meteorological factors were presented in Fig. . From 31 October to 6 November, the average NR-PM1 mass concentration was 32.3 μg/m3. The maximum concentration of 122.2 μg/m3 occurred on 4 November. The average PM1 mass concentrations were 69.5 μg/m3 during the autumn HPEs. The results show that the main composition of PM1 was OA, nitrate, sulphate, ammonium and chloride, and their average mass concentrations were 21.3, 25.1, 9.4, 12.3, and 1.2 μg/m3. In the autumn HPEs, the average median particle sizes of OA, nitrate, sulphate, ammonium and chloride were 472, 645, 531, 597 and 531 nm, respectively. The median particle size of OA was significantly lower than the other species, which indicated that OA might be fresher, or its growth was limited by hygroscopicity (Meier et al., Citation2009), oxidation, or mixing in the process of ageing so that its ageing degree was lower than those of nitrate, sulphate, and ammonium. The ageing of inorganic aerosol in transmission, touch and accumulation might be a major cause of its larger median particle size. The chloride was too less relative to other components, so this research did not discuss it in detail.

Fig. 3. Time series of mass concentration with the change of size of organics (a), nitrate (b), sulphate (c), ammonium (d) measured by AMS and number concentration of PM1 measured by TDMPS (e) at SDZ regional GAW station, PM2.5 measured at Beijing Wanliu station (f), meteorological characteristics such as wind direction and speed (g), PLAM and relative humidity (h) during autumn.

Fig. 3. Time series of mass concentration with the change of size of organics (a), nitrate (b), sulphate (c), ammonium (d) measured by AMS and number concentration of PM1 measured by TDMPS (e) at SDZ regional GAW station, PM2.5 measured at Beijing Wanliu station (f), meteorological characteristics such as wind direction and speed (g), PLAM and relative humidity (h) during autumn.

Table shows the changes of every component in the autumn HPEs. The average mass concentrations of organics, nitrate, sulphate, ammonium and PM1 were 5.7, 2.1, 0.97, 1.2 and 10.3 μg/m3, respectively, in clean periods, from 31 October to 2 November, when the mass concentration was less than 15 μg/m3, which is usually a lower mode value of mass concentration during the observation period, and is defined as the threshold of mass concentration for clean periods. Several new particle formation (NPF) events were observed in this study during the late autumn of 2015, however, only the NPF event started at 12:00 on 2 November was followed by two persistent HPEs. During HPEs, nitrate was the largest fraction of PM1 (36–38%), and the most significant increased component during HPEs relative to clean period (factors of 9.2 to 15, Table ). The proportion of OA in PM1 was similar to that of nitrate, but the growth rate in HPE was much smaller, mainly because of the high concentration of OAs during cleaning period. Compared with the largest increase of nitrate, the similar increase was found for ammonium product, showing that a large amount of ammonium was formed by the combination of a large amount of in HPEs. We analysed the correlations among the and , . The results show that good correlation between and , and , inferring that NH4NO3 and (NH4)2SO4 are all the major components for ammonium product. But the correlation coefficient between the and was 0.995 and 0.993 during autumn and winter, respectively, which are higher than the correlation between the and , with correlation coefficient of 0.968 and 0.973 during autumn and winter. These results showed that ammonium was more likely in combination with nitrate instead of sulphate in this study. Moreover, a previous study showed that higher RH may enhance NH4NO3 generation via aqueous-phase processes (Li et al., Citation2014). The average RH during autumn HPEs was 85%, much higher than clean period (44%), which may lead to more NH4NO3 in pollution. Although the neutralization of acidic sulphate is favoured over the formation of ammonium nitrate, but due to the significant reduction of SO2 emissions by power plant desulphurization since 2002 in China (Cao et al., Citation2010). In the presence of and lack of , ammonium nitrate was the largest fraction of PM1 in this study. Less ammonium sulphate can lead to more ammonium nitrate being formed due to sufficient NH3 left. Generally, ammonium nitrate is formed in areas characterized by high NH3 and HNO3 concentrations and low temperatures (Pathak et al., Citation2004).

Table 1. Mass concentration and increased percentage of chemical components of PM1 during autumn HPEs.

The key here is using one of the comprehensive meteorological index (PLAM), which mainly indicates the regional atmospheric stability and air condensation ability, to estimate the relative contribution of meteorological factor on the increasing PM2.5 in the first several hours when the pollution is formed. The PLAM was established based on a set of weather sensitive parameters such as air pressure, air temperature, winds, relative humidity, evaporability, stability, and effective parameter associated with the contribution of air pollution, respectively; and the PLAM was proved to be able to quantitatively distinguish the relative contribution of meteorological factor in the total concentration change of PM (Zhang et al., Citation2009).

The PLAM index has been developed to diagnose the air quality in Beijing and other cities derived from a correlation between PM2.5 and relevant weather elements based on the data between 2007 and 2008 (Wang et al., Citation2013a). Because significant correlations exist between PM2.5 concentrations and the variation in the sensitive meteorological parameters such as atmospheric condensation function, adaptive weight parameter β in both summer and winter, which indicate the meteorological conditions that affect the concentration of PM2.5. Fig. displays the correlation analysis between the 24-hforecast of PLAM over the 731 days from 1 January 2007 to the 31 December 2008 and the PM2.5 observations (indicated by the solid dotted line) in Beijing (Wang et al., Citation2013a, Citation2013b). There is also a linear correlation, indicated that it is meaningful for analysing, diagnosing and forecasting the impact of meteorological conditions on AQ based on PM2.5. The rise of PM2.5 due to meteorological factors can be expressed by the magnitude of the rise of PLAM (Zhang et al., Citation2015a)

Fig. 4. Correlation between PLAM index 24-h forecast and PM2.5 during 2007–2008 (CitationWang et al., 2013a).

Fig. 4. Correlation between PLAM index 24-h forecast and PM2.5 during 2007–2008 (CitationWang et al., 2013a).

PLAM index started its rise from a very low value at ~12:00 of 2 November to its first higher platform value of 60 in several hours till the early morning of 3 November, which was accompanied by the emergence of small south wind, climbing RH, as well as PM2.5 mass concentrations at Beijing urban area increased from 40 to ~120 μg/m3. At the same time, this NPF event was basically completed (Maso et al., Citation2005). We called this period as primary stage of pollution formation (P-HPE or NPF in Fig. ).

After PM2.5 pollution was formed and maintained at a high level for some time, still associating with southerly wind (Fig. ), one can find that RH was climbing towards its next summit and the atmosphere was even more stable with value of PLAM rose from 50 to its peak value of 100, which had exceeded a threshold value of 80 of PLAM leading to have more secondary aerosol formed to have visibility smaller than 10 km, and normally corresponding to heavy haze episode in BIV (Zhang et al., Citation2013a). We called this stage as HPE1. The HPE1 was from 9 am of 3 November to 10 am of 4 November. Since the PLAM values almost have a linear relationship with the mass concentrations of PM, the increased per cent of PLAM was 100% (from 50 to 100), while the observed total uplift percentage of PM2.5 mass was 122% (from 90 to 200) during these several hours period of HPE1, which implicating that about 80% of this PM mass increase during HPE1 can be considered as a positive feedback effect that come from even more stable atmosphere that derived from the interaction between formed pollution and the boundary layer (BL) structure in the lower atmosphere; Since the rise of PLAM also means that the condensation rate of water vapour will increase, more water vapour will condense on the surface of the aerosol to increase liquid water content, which will increase the formation of the secondary aerosol through the liquid phase reaction under the favourable meteorological condition provided by this feedback effect. The remaining ~20% of the increase can be attributable to the more chemical conversion. What’s the new mechanism of chemical conversion to form more secondary aerosol is worth further study.

Such a positive feedback effect of atmosphere on the PM pollution caused by the interaction of the pollution and the BL occurred again during the subsequent period, called HPE2 in Fig. . It is also a period of about half a day, stared from 10 am of 4 November to 8 am of 5 November. During this episode, PLAM rose from 70 to its highest value of 100, increasing ~40%, and PM2.5 increased from 180 to 250 with similar percentage.

3.1.2. Winter

The mass concentrations with the change of particle size for organics, nitrate, sulphate, ammonium measured by AMS during winter, as well as number concentration of PM1 measured by TDMPS, and various meteorological factors were presented in Fig. . From 17 December to 27 December, the average NR-PM1 mass concentration was 44.2 μg/m3. The maximum concentration of 192.8 μg/m3 occurred on 25 December. The average PM1 mass concentrations were 72.5 μg/m3 during the winter HPEs. The average mass concentration of OA, nitrate, sulphate, ammonium and chloride were 21.3, 25.1, 9.4, 12.3 and 1.2 μg/m3. In the winter HPEs, the average median particle sizes of organic, nitrate, sulphate, ammonium and chloride particles were 569, 629, 731, 598 and 629 nm, respectively. Compared with autumn, the winter average particle size of OA and sulphate was larger, indicating that the secondary organic and sulphate aerosols were more aged than the other particles. The nitrate particle size reduced because the formation of nitrate particles is associated with the condition of low temperature and humidity.

Fig. 5. Time series of mass concentration with the change of size of organics (a), nitrate (b), sulphate (c), ammonium (d) measured by AMS and number concentration of PM1 measured by TDMPS (e) at SDZ regional GAW station, PM2.5 measured at Beijing Wanliu station (f), meteorological characteristics such as wind direction and speed (g), PLAM and relative humidity (h) during winter.

Fig. 5. Time series of mass concentration with the change of size of organics (a), nitrate (b), sulphate (c), ammonium (d) measured by AMS and number concentration of PM1 measured by TDMPS (e) at SDZ regional GAW station, PM2.5 measured at Beijing Wanliu station (f), meteorological characteristics such as wind direction and speed (g), PLAM and relative humidity (h) during winter.

Table shows the changes of every component in the winter HPEs. During the wintertime in this campaign, the average mass concentrations of organics, nitrate, sulphate, ammonium and PM1 were 4.3, 0.90, 1.1, 0.68 and 7.1 μg/m3, respectively, in clean periods, from 16 December and 18 December. After a NPF event ended in the early hours of 19 December 2015, three HPEs (called HPE3, HPE4 and HPE5) were observed thereafter (Fig. ). Similar to autumn, nitrate was also the most significant increase component during HPEs relative to clean period (factors of 19 to 24, Table ), and ammonium had similar rise amplitude also. But the difference is that organics was the largest fraction of PM1 (34–42%).

Table 2. Mass concentration and increased percentage of chemical components of PM1 during winter HPEs.

During the three HPEs in wintertime (Fig. ), one can see the continuous accumulation and becoming more serious aerosol pollution processes, and their interactions with meteorological condition (Fig. ). From 2 pm of 18 December to 6 am of 19 December (marked as P-HPE in Fig. ), a NPF event occurred and eventually laid the foundation of particle for the subsequent constantly aggravated pollution episodes. Along with the PLAM started its rise from a very low value to near its threshold value of 80 in several hours during and after the NPF event, the mass concentration of PM2.5 also increased to 150 μg/m3 and maintain at this level, which led to the primary stage of pollution formation. Within the thereafter HPE3, from 12 am of 19 December to 9 pm of 21 December, the atmosphere was becoming more and more stable alone with the emergence and aggravation of pollution. This can be found that the PLAM value rose from the initial around 60 to its first high value of 170 with ~180% increase, with synchronously increase ~220% of PM2.5 mass concentration at Beijing urban area (PM2.5 mass rose from the initial ~110 μg/m3 and climbing to the first high value of about 350 μg/m3), which implicating that still ~80% of this PM mass increase during HPE3 can be considered as a positive feedback effect that come from interaction between formed pollution and the lower atmosphere. And the remaining ~20% of the increase can be attributed to the role of chemical mechanisms in this stage in BIV. This result is similar to those during late autumn previous in this paper.

In other HPEs (HPE4, HPE5) in wintertime, the estimates of the meteorological contribution on the increasing PM2.5 mass were also around 80%. In order to provide further evidence to support the value of 80% of increased PM2.5 be attributable to the contribution from meteorological factors in accordance with the way of statistical analysis. The correlation between increase of PLAM and PM2.5 during all HPEs is presented in Fig. . The correlation coefficient of 0.82 between the increased PM2.5 and PLAM can be considered as the evidence to support the claim of PLAM explaining 80% of PM2.5 increase above.

Fig. 6. Correlation between the PLAM and PM2.5 during all HPEs in autumn and winter of 2015.

Fig. 6. Correlation between the PLAM and PM2.5 during all HPEs in autumn and winter of 2015.

An interesting phenomenon is that southerly wind dominated in the polluted episodes (Fig. ), exhibiting again the pollutants from southern part of BIV had obvious contribution to aerosol pollution in Beijing area. However, during HPE5 (see in Fig. ), the prevailing wind is westerly wind that is also able to bring pollutants to the BIV area, the source will be discussed in details later.

3.2. Sources of pollutants

To explore the regional sources that influence on PM1 loading and composition during the HPEs, we used the Traject plugin of Meteoinfo (Wang et al., Citation2009) to produce back trajectories of air masses arriving in the SDZ area. We calculated 24-h back trajectories starting at the SDZ ground level (117.07°E, 40.39°N, 287 m a.s.l) and marked the trajectories in colour from shallow to deep according to the time sequence. This allows us to obtain air mass trajectories before, during and after pollution episodes.

The main source of anthropogenic SO2 is coal combustion, which is mainly from the stationary sources. NOx is mainly from the high temperature combustion of motor vehicles using either gasoline or diesel fuel, but also from coal combustion. Based on Chinese emissions inventory developed by Tsinghua University in 2010 (the Multi-resolution Emissions Inventory for China, MEIC) (He and Zhang, Citation2012), the ratios of NOx/SO2 from stationary and mobile sources were 0.98 and 28.7, respectively. NOx and SO2 are the major precursors of the secondary inorganic components and . So the / ratio usually can be used to roughly determine the stationary and mobile sources.

In the autumn clean period that before the start of pollution episodes, main air masses came from north-west region in Inner Mongolia, containing a large proportion of organics with lower / ratio of 2.3 (Fig. ). During the HPEs (HPE1 and HPE2), the air mass was mainly from the south of the station with increased / ratio of 3.6, showing the increased nitrate from southern Beijing and Hebei province. The ratio of / was 2.3 on clean days and 3.6 on the HPEs, both of which were close to the ratio of NOx/SO2 emitted from stationary coal combustion sources. According to the result, the contribution of stationary sources such as coal combustion, biomass burning and industry emissions was around 80–90% in the autumn HPEs. Therefore, our results suggest that coal combustion still was the main source of the increased NOx and nitrate during the HPEs at SDZ. Comparing the nitrate mass in PM10 among the stations within BIV in 2013 (Zhang et al., Citation2015a), the nitrate concentrations in urban Beijing (~20 μg/m3) were not higher than those in other stations (19–22 μg/m3 for Zhengzhou, Gucheng, and XiAn). However, the contribution from motor vehicle emission is the largest in Beijing (NBS-Beijing, Citation2013). This result shows that nitrate over the region does not only originate from motor vehicle emission but also from coal combustion. Even in BIV, about 28% of the emission amount of NOx can still be attributed to coal combustion. In other areas in China, about 50% of nitrate mass can be attributed to coal combustion (Cao et al., Citation2010). Globally, NOx is mainly derived from natural sources (Miyazaki et al., Citation2012). However, in urban areas, NOx mainly comes from anthropogenic activities such as oil and coal combustion. Given the rapid development of China’s economy in recent years, the amount of gasoline and diesel oil consumption has significantly increased over the past few years (Lu et al., Citation2011), even if the Chinese government started promoting the ‘low nitrogen combustion technology’ for denaturalization in power generation since 2012. The installation capacity reached only around 70% until the middle of 2013, except in BIV, which can probably explain the high nitrate concentration in this study. In addition, the highest hygroscopic of nitrate (Che et al., Citation2016) can be explained by the large increase in nitrate in the HPEs at this study. During the ‘after pollution’ stage (Fig. ), air mass came mainly from the north-east of the station with fewer pollutants associated.

Similar to autumn, the clean air masses were mainly from the north of the station, especially from the north-west, before and after pollutions in winter (Fig. ), associated with the averaged / ratios of 1.35. During the HPEs, the polluted air masses were mainly from the south of the station, especially from the southwest. The averaged / ratio was 1.82. Compared with autumn, the contribution of stationary sources was higher, more than 90%, especially the coal combustion. In addition, the biggest difference with autumn is that there is a part of polluted air masses coming from the west of the station, associated with relative higher mass of sulphate, nitrate and OA, and the 0.91 of / ratio value. This exhibits that the westerly wind from central Inner Mongolia and north Shanxi can also carry air pollutants mostly derived from coal combustion, contributing to the heavy pollution in the northern BIV area in winter, resulted in higher sulphate, nitrate and OA mass.

4. Conclusions

In this study, we conducted an investigation of size-resolved chemical components of PM1 from a site in northern Beijing area during HPEs in autumn and winter of 2015. The HPEs were always formed under relatively stable weather conditions with a lower speed and southerly wind and started climbing relative humidity; once the formation of pollution, it will further stabilize the atmosphere, increase condensation rate of water vapour on the surface of the aerosol to have more secondary aerosol formation through the liquid phase reaction, and feedback to a substantial increase in the concentration of fine aerosol particles. This feedback effect of meteorological conditions is considered to contribute approximately 80% to the ‘explosive’ increase of aerosol masses within the first several hours to 10 hours after the formation of aerosol pollution.

Nitrate was the most significantly increased component during HPEs relative to clean period. The mass concentration of nitrate was ~9.2 to 15 times higher during HPEs in autumn, and the rise was even higher in wintertime HPEs with factors of 19–24. Nitrate and OAs were the largest fractions of PM1 during HPEs in autumn and winter, respectively. Ammonium showed a similar percentage increase with nitrate during HPEs whether in winter or autumn, exhibiting that a large amount of ammonium was formed by the combination of in HPEs.

At the northern site of BIV, regional transport with air masses coming from the urbanized region of Jing (Beijing)–Jin (Tianjin)–Ji (alias of Hebei province) and its nearby surrounding regions (3JNS) exhibited a major contribution to nitrate, sulphate and OA, especially in autumn and winter. In winter, the westerly wind from central Inner Mongolia and northern Shanxi province can also carry air pollutants that have an influence on the PM1 mass in BIV area. Coal-combustion was still a major source of these pollutants.

Disclosure statement

No potential conflict of interest was reported by the authors.

Funding

This research was supported by National Key Project of MOST [grant number 2016YFC0203306] and Specific Team Fund of the Jiangsu Collaborative Innovation Center for Climate Change.

Acknowledgments

We sincerely thank the support from staffs from Shangdianzi station of CAWNET, CMA.

References

  • Aiken, A. C., DeCarlo, P. F., Kroll, J. H., Worsnop, D. R., Huffman, J. A. and co-authors. 2008. O/C and OM/OC ratios of primary, secondary, and ambient organic aerosols with high-resolution time-of-flight aerosol mass spectrometry. Environ. Sci. Technol. 42, 4478–4485.10.1021/es703009q
  • Aiken, A. C., Salcedo, D., Cubison, M. J. and Huffman, J. A. 2009. Mexico City aerosol analysis during MILAGRO using high resolution aerosol mass spectrometry at the urban supersite – Part 1: Fine particle composition and organic source apportionment. Atmos. Chem. Phys. 9, 6633–6653.10.5194/acp-9-6633-2009
  • Allan, J. D., Bower, K. N., Coe, H., Boudries, H., Jayne, J. T. and co-authors. 2004. Submicron aerosol composition at Trinidad Head, California, during ITCT 2K2: Its relationship with gas phase volatile organic carbon and assessment of instrument performance. J. Geophys. Res. Atmos. 109, 1085–1096.
  • Allan, J. D., Jimenez, J. L., Williams, P. I., Rami, A. M., Bower, K. N. and co-authors. 2003. Quantitative sampling using an aerodyne aerosol mass spectrometer 1. Techniques of data interpretation and error analysis. J. Geophys. Res. Atmos. 0345, 107–118.
  • An, X., Zhou, L., Yao, B., Xu, L. and Ma, L. 2012. Analysis on source features of halogenated gases at Shangdianzi regional atmospheric background station. Atmos. Environ. 57, 91–100.10.1016/j.atmosenv.2012.04.042
  • Canagaratna, M. R., Jayne, J. T., Jimenez, J. L., Allan, J. D., Alfarra, M. R. and co-authors. 2007. Chemical and microphysical characterization of ambient aerosols with the aerodyne aerosol mass spectrometer. Mass Spectrom. Rev. 26, 185–222.10.1002/(ISSN)1098-2787
  • Cao, G. L., An, X. Q., Zhou, C. H., Ren, Y. Q. and Tu, J. 2010. Emission inventory of air pollutants in China. Zhongguo Huanjing Kexue 30, 900–906.
  • Che, H. C., Zhang, X. Y., Wang, Y. Q., Zhang, L., Shen, X. J. and co-authors. 2016. Characterization and parameterization of aerosol cloud condensation nuclei activation under different pollution conditions. Sci. Rep. 6, 9369.10.1038/srep24497
  • Choobari, O. A., Zawar-Reza, P. and Sturman, A. 2012. Feedback between windblown dust and planetary boundary-layer characteristics: Sensitivity to boundary and surface layer parameterizations. Atmos. Environ. 61, 294–304.10.1016/j.atmosenv.2012.07.038
  • DeCarlo, P. F., Kimmel, J. R., Trimborn, A., Northway, M. J., Jayne, J. T. and co-authors. 2006. Field-deployable, high-resolution, time-of-flight aerosol mass spectrometer. Anal. Chem. 78, 8281–8289.10.1021/ac061249n
  • Gao, Y., Zhang, M., Liu, Z., Wang, L., Wang, P. and co-authors. 2015. Modeling the feedback between aerosol and meteorological variables in the atmospheric boundary layer during a severe fog-haze event over the North China Plain. Atmos. Chem. Phys. 15, 1093–1130.
  • He, K. and Zhang, Q. 2012. Multi-resolution Emission Inventory for China (MEIC): Model framework and 1990–2010 anthropogenic emissions. In: Proceedings of the AGU Fall Meeting. Lisbon.
  • Huang, X. F., He, L. Y., Hu, M., Canagaratna, M. R., Kroll, J. H. and co-authors. 2011. Characterization of submicron aerosols at a rural site in Pearl River Delta of China using an Aerodyne High-Resolution Aerosol Mass Spectrometer. Atmos. Chem. Phys. 11, 1865–1877.10.5194/acp-11-1865-2011
  • Jayne, J. T., Worsnop, D. R., Kolb, C. E., Leard, D., Davidovits, P. and co-authors. 2000. Development of an aerosol mass spectrometer for size and composition analysis of submicron particles. Aerosol Sci. Technol. 33, 49–70.10.1080/027868200410840
  • Jimenez, J. L., Jayne, J. T., Shi, Q., Kolb, C. E., Worsnop, D. R. and co-authors. 2003. Ambient aerosol sampling using the aerodyne aerosol mass spectrometer. J. Geophys. Res. 108, 447–457.
  • Kimmel, J. R., Farmer, D. K., Cubison, M. J., Sueper, D., Tanner, C. and co-authors. 2011. Real-time aerosol mass spectrometry with millisecond resolution. Int. J. Mass Spectrom. 303, 15–26.10.1016/j.ijms.2010.12.004
  • Li, Y., Schwandner, F. M., Sewell, H. J., Zivkovich, A., Tigges, M. and co-authors. 2014. Observations of ammonia, nitric acid, and fine particles in a rural gas production region. Atmos. Environ. 83, 80–89.10.1016/j.atmosenv.2013.10.007
  • Lu, Z., Zhang, Q. and Streets, D. G. 2011. Sulfur dioxide and primary carbonaceous aerosol emissions in China and India, 1996–2010. Atmos. Chem. Phys. 11, 9839–9864.10.5194/acp-11-9839-2011
  • Malloy, Q. G. J., Nakao, S., Qi, L., Austin, R., Stothers, C. and co-authors. 2009. Real-time aerosol density determination utilizing a modified scanning mobility particle sizer—aerosol particle mass analyzer system. Aerosol Sci. Technol. 43, 673–678.10.1080/02786820902832960
  • Maso, M. D., Kulmala, M., Riipinen, I. and Wagner, R. 2005. Formation and growth of fresh atmospheric aerosols: Eight years of aerosol size distribution data from SMEAR II, Hyytiälä, Finland. Boreal Environ. Res. 10, 323–336.
  • Meier, J., Wehner, B., Massling, A. and Birmili, W. 2009. Hygroscopic growth of urban aerosol particles in Beijing (China) during wintertime: A comparison of three experimental methods. Atmos. Chem. Phys. 9, 6865–6880.10.5194/acp-9-6865-2009
  • Middlebrook, A. M., Bahreini, R., Jimenez, J. L. and Canagaratna, M. R. 2012. Evaluation of composition-dependent collection efficiencies for the aerodyne aerosol mass spectrometer using field data. Aerosol Sci. Technol. 46, 258–271.10.1080/02786826.2011.620041
  • Miyazaki, K., Eskes, H. J. and Sudo, K. 2012. Global NOx emission estimates derived from an assimilation of OMI tropospheric NO2 columns. Atmos. Chem. Phys. 11, 31523–31583.
  • NBS-Beijing. 2013. Beijing Statistical Yearbook. Beijing.
  • Pathak, R. K., Yao, X. H. and Chan, C. K. 2004. Sampling artifacts of acidityand ionic species in PM2.5. Environ. Sci. Technol. 38, 254–259.10.1021/es0342244
  • Quan, J., Tie, X., Zhang, Q., Liu, Q., Li, X. and co-authors. 2014. Characteristics of heavy aerosol pollution during the 2012–2013 winter in Beijing, China. Atmos. Environ. 88, 83–89.10.1016/j.atmosenv.2014.01.058
  • Schneider, J., Favez, O., George, C. and Baltensperger, U. 2013. Atmospheric chemistry and physics characterization of aerosol chemical composition with aerosol mass spectrometry in Central Europe: An overview. Proc. Natl. Acad. Sci. 72, 794–798.
  • Shen, X., Sun, J., Zhang, Y., Zhang, L., Zhou, H. and co-authors. 2014. The characteristics of particle number size distribution under typical meteorological conditions at Shangdianzi Regional Station in Beijing. Adv. Meteorol. Sci. Technol. 1, 31–37.
  • Sun, Y. L., Wang, Z. F., Fu, P. Q., Yang, T., Jiang, Q. and co-authors. 2013. Aerosol composition, sources and processes during wintertime in Beijing, China. Atmos. Chem. Phys. 13, 4577–4592.10.5194/acp-13-4577-2013
  • Takegawa, N., Miyakawa, T., Kondo, Y., Jimenez, J. L., Zhang, Q. and co-authors. 2006. Seasonal and diurnal variations of submicron organic aerosol in Tokyo observed using the Aerodyne aerosol mass spectrometer. J. Geophys. Res. Atmos. 111, 1937–1952.
  • Takegawa, N., Miyazaki, Y., Kondo, Y., Komazaki, Y., Miyakawa, T. and co-authors. 2005. Characterization of an aerodyne aerosol mass spectrometer (AMS): Intercomparison with Other Aerosol Instruments. Aerosol Sci. Technol. 39, 760–770.10.1080/02786820500243404
  • Tang, L. L., Zhang, Y. J., Sun, Y. L., Hongxia, Y. U., Zhou, H. C. and co-authors. 2014. Components and optical properties of submicron aerosol during the lasting haze period in Nanjing. Chinese J. 59, 1955–1966.
  • Tuch, T. M., Haudek, A., Müller, T. and Nowak, A. 2009. Design and performance of an automatic regenerating adsorption aerosol dryer for continuous operation at monitoring sites. Atmos. Meas. Tech. 2, 417–422.
  • Wang, J. Z., Gong, S. L., Zhang, X. Y., Yang, Y. Q., Hou, Q., Zhou, C. H. and co-authors. 2013a. A parameterized method for air-quality diagnosis and its applications. Adv. Meteorol. 2012.
  • Wang, Z., Li, J., Wang, Z., Yang, W., Tang, X. and co-authors. 2014a. Modeling study of regional severe hazes over mid-eastern China in January 2013 and its implications on pollution prevention and control. Sci. China Earth Sci. 57, 3–13.10.1007/s11430-013-4793-0
  • Wang, H., Shi, G. Y., Zhang, X. Y., Gong, S. L., Tan, S. C. and co-authors. 2015a. Mesoscale modelling study of the interactions between aerosols and PBL meteorology during a haze episode in China Jing–Jin–Ji and its near surrounding region – Part 2: Aerosols’ radiative feedback effects. Atmos. Chem. Phys. 15, 3277–3287.10.5194/acp-15-3277-2015
  • Wang, J., Wang, S., Jiang, J., Ding, A., Zheng, M. and co-authors. 2014b. Impact of aerosol-meteorology interactions on fine particle pollution during China’s severe haze episode in January 2013. Environ. Res. Lett. 9, 94002–94008( 94007).
  • Wang, J., Wang, Y., Liu, H., Yang, Y., Zhang, X. and co-authors. 2013b. Diagnostic identification of the impact of meteorological conditions on PM2.5 concentrations in Beijing. Atmos. Environ. 81, 158–165.10.1016/j.atmosenv.2013.08.033
  • Wang, H., Xue, M., Zhang, X. Y., Liu, H. L., Zhou, C. H. and co-authors. 2015b. Mesoscale modeling study of the interactions between aerosols and PBL meteorology during a haze episode in Jing–Jin–Ji (China) and its nearby surrounding region – Part 1: Aerosol distributions and meteorological features. Atmos. Chem. Phys. 15, 3257–3275.10.5194/acp-15-3257-2015
  • Wang, Y. Q., Zhang, X. Y. and Draxler, R. R. 2009. TrajStat: GIS-based software that uses various trajectory statistical analysis methods to identify potential sources from long-term air pollution measurement data. Environ. Model. Softw. 24, 938–939.10.1016/j.envsoft.2009.01.004
  • Weimer, S., Drewnick, F., Hogrefe, O., Schwab, J. J., Rhoads, K. and co-authors. 2006. Size-selective nonrefractory ambient aerosol measurements during the Particulate Matter Technology Assessment and Characterization Study – New York 2004 Winter Intensive in New York City. J. Geophys. Res. Atmos. 111, 2187–2199.
  • Xu, X. 1992. Characterization of the regional summer background aerosol of Beijing, China. Acta Meteorol. Sinica 50, 310–319.
  • Zhang, Y. H., Hu, M., Zhong, L. J., Wiedensohler, A., Liu, S. C. and co-authors. 2008. Regional integrated experiments on air quality over Pearl River Delta 2004 (PRIDE-PRD2004): Overview. Atmos. Environ. 42, 6157–6173.10.1016/j.atmosenv.2008.03.025
  • Zhang, Q., Jimenez, J. L., Canagaratna, M. R., Allan, J. D., Coe, H. and co-authors. 2007. Ubiquity and dominance of oxygenated species in organic aerosols in anthropogenically-influenced Northern Hemisphere midlatitudes. Geophys. Res. Lett. 34, 173–180.
  • Zhang, X. Y., Sun, J. Y., Wang, Y. Q., Weijun, L. I., Zhang, Q. and co-authors. 2013a. Factors contributing to haze and fog in China. Chinese J. 58, 1178–1187.
  • Zhang, Y. M., Sun, J. Y., Zhang, X. Y., Shen, X. J., Wang, T. T. and co-authors. 2013b. Seasonal characterization of components and size distributions for submicron aerosols in Beijing. Sci. China Earth Sci. 56, 890–900.10.1007/s11430-012-4515-z
  • Zhang, X. Y., Wang, Y. Q., Lin, W. L., Zhang, Y. M., Zhang, X. C. and co-authors. 2009. Changes of atmospheric composition and optical properties over BEIJING-2008 Olympic Monitoring Campaign. Bull. Amer. Meteor. Soc 90, 584–595.
  • Zhang, X. Y., Wang, J. Z., Wang, Y. Q., Liu, H. L., Sun, J. Y. and co-authors. 2015a. Changes in chemical components of aerosol particles in different haze regions in China from 2006 to 2013 and contribution of meteorological factors. Atmos. Chem. Phys. 15, 12935–12952.10.5194/acp-15-12935-2015
  • Zhang, Y. M., Zhang, X. Y., Sun, J. Y., Hu, G. Y., Shen, X. J. and co-authors. 2014. Chemical composition and mass size distribution of PM1.0 at an elevated site in central east China. Atmos. Chem. Phys. 14, 15191–15218.
  • Zhang, Y. W., Zhang, X. Y., Zhang, Y. M., Shen, X. J., Sun, J. Y. and co-authors. 2015b. Significant concentration changes of chemical components of PM1 in the Yangtze River Delta area of China and the implications for the formation mechanism of heavy haze-fog pollution. Sci. Total Environ. 538, 7–15.10.1016/j.scitotenv.2015.06.104
  • Zheng, G. J., Duan, F. K., Ma, Y. L., Cheng, Y., Zheng, B. and co-authors. 2015. Exploring the severe winter haze in Beijing. Atmos. Chem. Phys. Discuss. 14, 17907–17942.