Source apportionment of carbonaceous aerosols in Beijing with radiocarbon and organic tracers: insight into the differences between urban and rural sites

. Carbonaceous aerosol is a dominant component of ﬁne particles in Beijing. However, it is challenging to appor-tion its sources. Here, we applied a newly developed method which combined radiocarbon ( 14 C) with organic tracers to apportion the sources of ﬁne carbonaceous particles at an urban (IAP) and a rural (PG) site of Beijing. PM 2 . 5 ﬁlter samples (24 h) were collected at both sites from 10 November to 11 December 2016 and from 22 May to 24 June 2017. 14 C was determined in 25 aerosol samples (13 at IAP and 12 at PG) representing low pollution to haze conditions. Biomass burning tracers (levoglucosan, mannosan, and galactosan) in the samples were also determined using gas chromatography–mass spectrometry (GC-MS). Higher contributions of fossil-derived OC (OC f ) were found at the urban site. The OC f / OC ratio decreased in the summer samples (IAP: 67.8 ± 4.0 % in winter and 54.2 ± 11.7 % in summer; PG: 59.3 ± 5.7 % in winter and 50.0 ± 9.0 % in summer) due to less consumption of coal in the warm season. A novel extended Gelencsér (EG) method incorporating the 14 C and organic tracer data was developed to estimate the fossil and non-fossil sources of primary and secondary OC (POC and SOC). It showed that fossil-derived POC was the largest contributor to OC (35.8 ± 10.5 % and 34.1 ± 8.7 % in wintertime


Introduction
Carbonaceous aerosols, often one of the most abundant components (20 %-80 %) in atmospheric aerosol particles, have a crucial impact on the global climate, air quality, and human health (He et al., 2001;Huang et al., 2014;Jimenez et al., 2009;Song et al., 2007;Zhou et al., 2018).The total content of carbonaceous aerosols (i.e., total carbon, TC) can be divided into organic carbon (OC) and elemental carbon (EC) according to their physical, chemical, and optical properties.The source of EC is from incomplete combustion of fossil fuel or biomass, while OC mainly originates from primary emissions from sources such as coal combustion, traffic emissions, cooking, and biomass burning, as well as from gas-particle conversion (Yang et al., 2016).It is very challenging to quantify the contributions from different sources to OC and EC because of the limited information on the sources, atmospheric loading, and composition of organic aerosols (Huang et al., 2014).Radiocarbon ( 14 C) analysis is a powerful tool for the quantification of fossil and non-fossil contributions to carbonaceous aerosols, as non-fossil sources contain a high contemporary 14 C content, while the fossil fractions are free of 14 C (Zotter et al., 2014;Bernardoni et al., 2013;Liu et al., 2013;Szidat et al., 2004Szidat et al., , 2006Szidat et al., , 2009)).A previous study in north-east China found a dominant fossilfuel contribution to EC (76 ± 11 %) and found that non-fossil sources are major contributors to OC (66 ± 11 %) (Zhang et al., 2016).Non-fossil sources of OC were major contributors to the fine particle pollution in Beijing during the APEC summit (Liu et al., 2016a).Moreover, clear seasonal trends of non-fossil and fossil source contributions to waterinsoluble OC (WINSOC) and water-soluble OC (WSOC) were found.Non-fossil sources were the major contributor (59 %) to WINSOC in summer and autumn, whereas fossil fuel emissions were predominant in winter and spring (Liu et al., 2013).Proportions of non-fossil sources in TC and WSOC associated with biogenic emissions increased during spring and summer, with maxima (85 % and 117 %, respectively) in May (Pavuluri et al., 2013).However, 14 C measurements do not permit direct discrimination of specific sources (e.g.biomass burning or secondary OC, SOC) of modern carbon.A combination with other techniques gives further insight into the characteristics of SOC (Minguillon et al., 2011;Szidat et al., 2006Szidat et al., , 2009;;Yttri et al., 2011).For example, applying Latin hypercube sampling with different OC/EC ratios, relative contributions of primary and secondary organic carbon were estimated (Zhang et al., 2016). 14C analysis combined with AMS-PMF (positive matrix factorization analysis of data from an online Aerodyne aerosol mass spectrometer) data has contributed to the identification of sources from primary emissions and secondary formation (Barrett et al., 2015;Zhang et al., 2018Zhang et al., , 2017;;Vlachou et al., 2018).
The Gelencsér method provides a first-order source apportionment of organic aerosol from fossil fuel combustion, biomass burning, biogenic emissions, and secondary organic aerosol, using measurements of specific organic tracers emitted by fossil and non-fossil sources and their OC/EC ratios derived from literature (Gelencsér et al., 2007).This method was derived in a European context, but for China the inclusion of food cooking and coal combustion is required.Based upon this consideration, an extended Gelencsér (EG) method that includes quantification of fossil and contemporary EC and OC by 14 C analysis has been developed in this study.The diversity of fuel types and combustion conditions makes the selection of OC/EC ratios for biomass burning difficult due to large uncertainties.For non-fossil sources of SOC, quantification by the method of Gelencsér et al. (2007) is totally dependent on the source apportionment of OC from biomass burning, and thus a cautious selection of ratios has been adopted in this study.
Biomass burning is an important source of both EC and OC, which can affect large areas of the world through long-range transport (Andreae and Merlet, 2001).It is also a key component when applying source apportionment by the Gelencsér method.As levoglucosan (1,6-anhydro-β-Dglucopyranose, LG) is an almost specific biomass burning tracer as the main pyrolysis product from cellulose (Puxbaum et al., 2007;Simoneit et al., 1999), concentrations of OC from biomass burning can be obtained by multiplying LG with suitable OC/LG ratios (Gelencsér et al., 2007;Zdrahal et al., 2002).However, the wide range of OC/LG ratios associated with changes in the biofuel types and combustion conditions cause great uncertainty in the estimation (Cheng et al., 2013;Fu et al., 2012;Gelencsér et al., 2007).To mitigate these differences from types of material and the burning conditions, a typical ratio of 12.2-12.5(12.3 on average) was documented in Andreae and Merlet (2001) by considering the biofuels of savanna and crop residues.This ratio is widely accepted in many studies (Andreae and Merlet, 2001;Fu et al., 2012;T. Zhang et al., 2008;Zhang et al., 2007) and has been used to estimate the contributions of biomass burning in Beijing, which ranged from 8 % to 50 % (Cheng et al., 2013;Kang et al., 2018;Liu et al., 2017;T. Zhang et al., 2008;Zhang et al., 2017).However, a single ratio is not representative of all local conditions, and more specific methods are needed. 14C analysis can provide accurate concentrations of EC from biomass burning, assuming that EC from non-fossil sources is exclusively from biomass burning.Introducing EC nf concentrations from 14 C analysis into the Gelencsér method of using OC/LG ratios provides valuable extra information and allows the method to be extended to include other sources.
Beijing, capital of China, has experienced severe PM 2.5 pollution for decades.This has been the subject of extensive research.However, few studies have sought to differentiate the fossil and non-fossil sources of SOC, even though they provide key information on the precursor and formation mechanisms of SOC.In this study, measurements of PM 2.5 , OC, and EC, along with biomass burning tracers, were conducted simultaneously at urban and rural sites of Beijing in the winter of 2016 and summer of 2017. 14C measurements of EC, OC, WINSOC, and WSOC were carried out subsequently on filter samples to enable source apportionment of fossil vs. non-fossil sources.A novel extended Gelencsér method combining 14 C analysis has been developed to explore the source apportionment of OC and EC, with SOC from fossil and non-fossil sources being quantified.The source apportionment results were compared with those by chemical mass balance (CMB).Correlations of WINSOC and WSOC with different sources of OC were also performed to study the formation mechanisms of SOC.

Aerosol sampling
Daily PM 2.5 (particles with aerodynamic diameter less than 2.5 µm) samples were collected at an urban site (39.98 • N, 116.39 • E; Institute of Atmospheric Physics, IAP) and a rural site (40.17 • N, 117.05 • E, Pinggu, PG) in Beijing during a winter campaign (10 November-11 December 2016) and a summer campaign (22 May-24 June 2017) as part of the Atmospheric Pollution and Human Health in a Chinese megacity (APHH China) programme; further information on the sampling sites is available in Shi et al. (2019).The urban site is a typical urban background site but may be subject to multiple local influences such as cooking emissions from nearby restaurants.The rural site is located in Pinggu District, close to a village surrounded by farmland.It is ∼ 60 km north-east of urban Beijing at the junction of Beijing, Tianjin, and Hebei provinces.A two-lane road is about 200-300 m north of the sampling site, but its traffic volume is relatively low.
High-volume air samplers (Tisch, USA) with a flow rate of 1.1 m 3 min −1 were collected on pre-combusted (450 • C, 6 h) quartz filters (Pallflex, 8 in × 10 in).Field blanks were collected by placing filters onto the filter holder for a few minutes without pumping before and after the campaign.After sampling, each exposed or blank filter was wrapped individually with aluminium foil and stored at −20 • C in the dark prior to analysis.Details of the sample collection are described elsewhere (Shi et al., 2019).

OC, EC, and major inorganic ions
OC and EC mass concentrations were determined with the DRI2015 carbon analyser with the EUSAAR_2 (European Supersites for Atmospheric Aerosol Research) transmittance protocol.Replicate analyses were conducted once every 10 samples.Blank samples (corresponding to 0.40 and 0.01 µg m −3 for OC and EC) were analysed to correct the sample results.The limits of detection of OC and EC were estimated to be 0.03 and 0.05 µg m −3 .Details of the OC/EC measurement method are described elsewhere (Paraskevopoulou et al., 2014).Major ions including SO 2− 4 , NO − 3 , NH + 4 , Na + , K + , and Cl − were determined on water extracts using an ion chromatograph (Dionex, Sunnyvale, CA, USA), with detection limits less than 0.01 µg m −3 .The uncertainties for OC and EC were less than 10 % and less than 5 % for inorganic ions (Xu et al., 2020).

Biomass burning tracers
The methodology to determine biomass burning tracers, including levoglucosan, mannosan, and galactosan, is described elsewhere (Fu et al., 2016).Recoveries for target compounds were better than 80 % as obtained by spiking standards to pre-combusted quartz filters followed by extraction and derivatization.Field blank filters were analysed using the procedure used for the samples above, but no target compounds were detected.Duplicate analyses showed analytical errors of less than 15 %.

Radiocarbon ( 14 C) analysis
The 14 C in total carbon (TC), water-insoluble TC (WIN-STC), and EC was determined on 25 (13 from IAP and 12 from PG) time-integrated high-volume PM 2.5 quartz fibre (QF) filter samples.Samples collected during both haze and non-haze days were selected in winter to better understand the pollution sources.PM 2.5 concentrations on 22 November and 1 December at IAP and PG sites were lower than 75 µg m −3 and regarded as non-haze air days, in contrast to other wintertime samples collected during haze pollution days.During the summer, typical samples were selected, with PM 2.5 concentrations of 42.5 ± 26.5 and 42.7 ± 21.2 µg m −3 at IAP and PG, respectively.The concentrations of PM 2.5 , EC, and OC and the corresponding nonfossil fractions of these selected days are shown in Table 2.
The method of 14 C measurement of carbonaceous aerosols has been described elsewhere (Agrios et al., 2015;Levin et al., 2010;Szidat et al., 2014;Vlachou et al., 2018;Zhang et al., 2012Zhang et al., , 2016;;Zotter et al., 2014) (Zhang et al., 2012;Szidat et al., 2014).The EC fraction was separated by an OC/EC analyser (Model 4L, Sunset Laboratory, USA) with the use of the Swiss_4S protocol (Zhang et al., 2012), which was coupled online with the MICADAS (Agrios et al., 2015).Each filter sample was extracted with water before the measurements to minimize the charring effect during the separation of EC from the WINSOC.
14 C results were expressed as fractions of modern (f M ), i.e., the fraction of the 14 C/ 12 C ratio of the sample related to that of the reference year 1950.The data analysis was carried out accounting for the blank correction (one field blank per site was analysed, not relevant for EC), decay of 14 C since the 1950s, nuclear bomb correction, charring of WINSOC (∼ 1 %), and EC yield after OC removal (IAP: 62 ± 6 %; PG: 76 ± 8 %) (Zhang et al., 2012(Zhang et al., , 2016;;Zotter et al., 2014).
Non-fossil fractions (f NF ) were determined from their corresponding f M values and reference values for pure nonfossil sources by Different values for the f NF,ref were applied for the bomb peak correction (Levin, et al. 2010).For EC, the f NF,ref is 1.10 ± 0.05 (Lewis et al., 2004;Palstra and Meijer, 2014), given that biomass burning is assumed to be the only nonfossil source of EC.For OC, it is calculated as where f M,bb and f M,bio are from biomass burning and biogenic sources, respectively, which are 1.10 ± 0.05 and 1.023 ± 0.015 (Lewis et al., 2004;Zotter et al., 2014), while p bio and p bb are the proportions of biogenic source and biomass burning, respectively, which are 0.1 and 0.9 in winter and 0.5 and 0.5 in summer (Levin et al., 2010).Analogously, the non-fossil fractions of OC, WSOC, and WINSOC (f NF,OC , f NF,WSOC , and f NF,WINSOC ) were calculated by following a mass-balance-like approach: where TC and EC are the concentrations of total and elemental carbon, respectively, and f NF,TC , f NF,EC , and f NF,WINSTC are the non-fossil fractions of TC, EC, and WINSTC, respectively.The fraction of fossil-fuel sources was calculated by The uncertainties were determined by error propagation.The mass concentration errors were assumed to be 10 % for EC and 6 % for OC and TC (typical values for EUSAAR2) (Zhang et al., 2016).

Extended Gelencsér method including 14 C data
An extended Gelencsér method including 14 C data was developed to quantify the fossil and non-fossil sources of primary and secondary OC (POC and SOC) along with OC from biomass burning and cooking (OC bb and OC ck ).The equations for the extended Gelencsér method are listed in Table 1.
The detailed selection of the OC/EC ratios will be discussed in Sect.3.2.

Chemical mass balance (CMB) model and AMS/ACSM-PMF analysis
Results on the same sets of samples from a chemical mass balance (US EPA CMB8.2) model and AMS/ACSM-PMF analysis (positive matrix factorization analysis of data from an online Aerodyne aerosol mass spectrometer at IAP and Aerosol Chemical Speciation Monitor at PG) were compared with 14 C-based source apportionment.Details on the experimental details and data analyses can be found in Xu et al. (2021) and Wu et al. (2020).The CMB utilizes a linear least-squares solution considering both uncertainties in source profiles and ambient measurements to ensure reliable fitting results.In order to better represent the source characteristics, the source profiles applied in this model were mostly from local studies in China (Cai et al., 2017;Wang et al., 2009;Y. X. Zhang et al., 2008;Zhang et al., 2007;Zhao et al., 2015), except vegetative detritus (Rogge et al., 1993;Wang et al., 2009).Details of the selection of the organic marker species can be found in Yin et al. (2010Yin et al. ( , 2015)).Experimental details on the AMS and ACSM-PMF method can be found elsewhere (Ng et al., 2011;Sun et al., 2016;Xu et al., 2019).The ACSM data were analysed for the mass concentrations and size distributions of nonrefractory submicron aerosol (NR-PM 1 ) species using the high-resolution data analysis software package PIKA (Sun et al., 2020).Positive matrix factorization was performed on high-resolution mass spectra of V mode and W mode to retrieve potential organic aerosol (OA) factors from different sources (Paatero and Tapper, 1994;Ulbrich et al., 2009).The organic matter / organic carbon (OM / OC) factor used for cooking OA (COA) is 1.38 (Xu et al., 2019).

Results and discussion
3.1 Overall results

Characteristics of PM 2.5 , OC and EC concentrations
Mass concentrations of PM 2.5 , OC, EC, and biomass burning tracers are shown in Fig. 1 and are summarized in Table 2, along with the meteorological conditions during the observation campaign in IAP and PG.The average concentrations of PM 2.5 were 91.2 ± 63.7 and 99.7 ± 77.8 µg m −3 at the IAP and PG sites, respectively, in winter, and 30.2 ± 14.8 and 27.5 ± 12.9 µg m −3 at the IAP and PG sites in summer.The highest 24 h concentration in winter is 239.9 µg m −3 (IAP) and 294.3 µg m −3 (PG), and more than 53 % and 46 % of the days have daily PM 2.5 concentrations higher than the Chinese air quality standard (PM 2.5 concentrations exceeding 75 µg m −3 are defined as haze conditions) in IAP and PG during the observation period.In summer, the air quality was improved, with PM 2.5 concentrations ranging from 12.2 to 78.8 and 11.6 to 70.3 µg m −3 in IAP and PG, respectively.Organic carbon and elemental carbon (OC and EC) are important constituents of PM 2.5 , accounting for 30.9 ± 9.3 % and 43.6 ± 17.9 % of PM 2.5 mass at the IAP and PG sites in winter and 26.8 ± 9.2 % and 37.3 ± 12.6 % in summer.The concentrations of EC showed a strong correlation with OC at both sites during the winter and summer (Table S1 in the Supplement).The average EC concentrations for 14 C analysis varied from 3.8 ± 2.1 µg m −3 (IAP) and 5.4 ± 2.6 µg m −3 (PG) in winter and 1.1 ± 0.3 µg m −3 (IAP) and 2.0 ± 0.7 µg m −3 (PG) in summer for the urban and rural sites, respectively.The mass concentration of OC for 14 C analysis was 4.1-44.9and 12.1-85.0µg m −3 at IAP and PG in winter and 4.7-12.7 and 6.2-17.9µg m −3 in summer.The selected samples are well representative as their concentrations were very close to those from the whole campaign.The average OC concentration for 14 C in winter was 3.2 and 4.3 times higher than in summer at IAP and PG, respectively.The OC/EC ratios at IAP and PG were in the range of 4.1-14.9and 6.2-14.6 in winter and 4.6-14.8and 4.4-28.3 in summer, which were all higher than 2.0 or 1.1 (Chow et al., 1996;Castro et al., 1999), suggesting an important contribution from secondary organic carbon (SOC).

Fossil and non-fossil sources of EC and OC based on radiocarbon ( 14 C) analysis
Figure 2 shows the absolute concentrations of fossil and nonfossil fractions of OC and EC and the relative contributions.EC f refers to EC from coal combustion and liquid fossil fuel (i.e., mainly vehicle emissions) and EC nf from biomass burning (Gray and Cass, 1998).As shown in Fig. 2, high concentrations of EC f and EC nf were found in both urban and rural sites in winter, suggesting elevated emissions from primary fossil and non-fossil sources like coal combustion and biomass burning.Residential coal consumption and biomass burning are still important in winter, especially in rural areas, due to intensive heating activities in the cold season.EC f contributed 7.6 ± 2.1 % and 6.0 ± 1.4 % to TC at the IAP and PG sites in winter and 6.9 ± 1.6 % and 8.9 ± 2.6 % in summer.Most of the fossil fractions of EC are within the range of previous studies in urban Beijing (67 %-96 % of EC) (Liu et al., 2020(Liu et al., , 2017;;Zhang et al., 2016Zhang et al., , 2017)).Higher contributions of EC f were found on the polluted days in wintertime ( OC/EC 6.9 ± 2.4 7.1 ± 0.7 10.8 ± 5.9 7.6 ± 2.2 7.1 ± 1.9 9.0 ± 1.9 9.3 ± 2.5 8.0 ± 1.0 9.0 ± 6.7 6.0 ± 1.3 LG  LG, MN, and GA referred to levoglucosan, mannosan, and galactosan, respectively; WS, wind speed; RH, relative humidity, T , temperature; NA = not available.due to the regional open burning activities during the postharvest period of wheat, which is common in North China (Li et al., 2016;Yan et al., 2006).
OC from fossil sources (OC f ) arises mainly from coal combustion and vehicle emissions, while OC from non-fossil sources (OC nf ) comes mainly from biomass burning, biogenic emissions, and cooking.Although fossil sources are the main contributor to OC (67.8 ± 4.0 % and 59.3 ± 5.7 % at IAP and PG sites in winter, compared to 54.2 ± 11.7 % and 50.0 ± 9.0 % in summer), relative contributions of OC nf were significantly increased in summer.These results are similar to previous studies in urban Beijing, which found a contribution of 66 ± 11 % ranging from 45 % to 82 % in PM 1 collected in 2013/14 (Zhang et al., 2015(Zhang et al., , 2017)).The increased contribution of OC nf in summer is likely to due to the enhancement of biomass burning and biogenic emissions (both primary and through SOC), as well as the decline in emissions from fossil sources such as coal burning, considering that cooking OC emission is unlikely to change with season.
Fossil and non-fossil sources of WSOC and WINSOC were also quantified by the 14 C measurement.Among the four fractions, WINSOC f had the highest contribution to TC in winter (37.9 ± 4.5 % at IAP and 36.2 ± 4.7 % at PG), followed by WSOC f (22.1 ± 5.2 %) at IAP and WINSOC nf (20.8 ± 3.4 %) at PG.In summer, the fraction of WINSOC f fell to 24.4 ± 9.3 % at IAP and 21.8 ± 4.6 % at PG, accompanied by increased fossil and non-fossil contributions from WSOC fractions.The increase of WSOC f fractions in summer implied an enhanced contribution from oxidized volatile organic compounds (VOCs) and aged primary fossil-derived OC, and WSOC nf is probably associated with biomass burning and secondary OC.Moreover, WINSOC was dominated by fossil sources at both sites in winter (IAP: 72.6 ± 3.6 %; PG: 63.4 ± 6.5 %), while the non-fossil fractions significantly increased from winter to summer.Similarly, WSOC was mainly fossil-derived in winter, while it tends to be accumulated in non-fossil fractions in summer (fossil fraction in WSOC, IAP: 60.5 ± 6.6 % in winter, 50.8 ± 12.3 % in summer; PG: 51.6 ± 8.7 % in winter, 47.7 ± 17.4 % in summer).The contribution of OC nf to OC increased with both nonfossil fractions in WINSOC and WSOC increasing.Even though WINSOC nf and WSOC nf cannot be attributed specifically to biomass burning, biogenic emissions, cooking, or secondary formation, it is likely that biogenic-derived POC and SOC make a pronounced contribution.Details of contributions of each fraction to primary and secondary OC will be discussed in Sect.3.2.

Source apportionment by an extended Gelencsér method
Gelencsér et al. ( 2007) reported a method for the source apportionment of carbonaceous aerosol into fractions from biomass burning, road traffic, and secondary organic aerosol, applicable to Europe where these are the dominant sources.
In order to use the same methodological concepts in China, the method required extending to include two further sources: coal combustion and cooking.To do so, the 14 C data were utilized.

Biomass burning
Levoglucosan (LG) is a typical biomass burning tracer, as the main pyrolysis product from cellulose.Much higher concentrations of LG were observed in the winter (311 ± 193 ng m −3 at IAP and 634 ± 483 ng m −3 at PG) than those in the summer (27.9 ± 29.6 ng m −3 at IAP and 74.0 ± 34.2 ng m −3 at PG).In addition, LG concentrations at the rural site were higher than those at the urban site in both winter and summer.This pattern is consistent with previous measurements in Table 3 (Chen et al., 2018;Kang et al., 2018;Li et al., 2018;Liu et al., 2016b;Salma et al., 2017;Sullivan et al., 2019;Yan et al., 2019;Zhu et al., 2017).The Pearson correlations of LG with PM 2.5 , OC, and EC at IAP and PG are shown in Table S1.During winter, LG correlated well with PM 2.5 , OC, and EC at PG, with correlation coefficients of 0.89, 0.89, and 0.81, respectively.These are higher than those at IAP (correlation coefficients of 0.56, 0.60, and 0.74, respectively), suggesting a more significant influence of biomass burning upon PM 2.5 in PG.The correlation coefficients in summer were much lower than those in winter for both sites, evidencing a reduced contribution of biomass burning activities to PM 2.5 .During the wintertime, the increasing use of biofuel for heating exacerbates the biomass pollution, and the stable atmospheric conditions also enhance the accumulation of LG (Shi et al., 2019).Compared with widely used cleaner fossil energy (i.e., natural gas, electricity, and liquefied petroleum gas) and renewable energy (i.e., solar energy) available in urban areas, rural households are still largely using straw and wood for cooking and heating (Hou et al., 2017).Moreover, open burning of crop residues during the post-harvest months (May to July and October to November in North China) in rural areas is still frequently performed in spite of prohibition by the government (Chen et al., 2017;Li et al., 2016).The water-soluble potassium ion (K + ) has been used as a biomass burning tracer previously, due to its good relationship with LG.However, in this study, the Pearson correlation coefficients of K + with LG are 0.51 and 0.86 in winter and 0.85 and 0.51 in summer for IAP and PG (Table S1), indicating other non-biomass sources of K + .Indeed, the sources of K + in the atmosphere are diverse, including sea salt, cooking, dust, coal combustion, and waste incineration, which makes K + less suitable as a biomass burning tracer (Zhang, et al., 2010).
According to Gelencsér et al. (2007), EC from biomass burning (EC bb ) can be derived by multiplying the LG concentrations first by (OC/LG) bb to give OC bb and then by (EC/OC) bb .However, the (OC/LG) bb ratio is highly variable depending on the type of material and the conditions of burning, and the ratios can vary by orders of magnitude.
The ratios of levoglucosan to mannosan (LG/MN) and to galactosan (LG/GA) can help us to infer the (OC/LG) bb and (EC/OC) bb ratios from certain types of biomass fuel (Kawamura et al., 2012).Figure 3 shows the source pro- LG/GA from different types of biomass burning emissions (Cheng et al., 2013;Sun et al., 2019a, b; LG, MN, and GA refer to levoglucosan, mannosan, and galactosan, respectively), including those measured in PM 2.5 samples at IAP and PG during winter and summer.The range of LG/MN is 7.66-13.41and 5.48-22.40for IAP in winter and summer and 2.45-60.20 and 3.27-17.71for PG in winter and summer, while the range of LG/GA is 3.98-6.92,2.78-101.43,1.97-34.66,and 1.51-13.39,respectively.file for LG/MN and LG/GA ratios measured in emissions from controlled biomass burning experiments in previous studies (Cheng et al., 2013;Sun et al., 2019a, b).There is a clear boundary of LG/MN ratios (∼ 10) between softwood and hardwood burning (Kawamura et al., 2012).Using LG/GA ratios can help to identify the burning of straws, woods, needles, and grasses, where most of the LG/GA ratios from burning of hardwood, straws, and grasses are higher than 10, quite different from those of needle plants.LG/MN and LG/GA ratios cannot be used to separate the burning of straws from hardwood, but the burning of hardwood can be neglected considering the local forest types around Beijing (Cheng et al., 2013).
The measured LG/MN and LG/GA ratios in this study (Fig. 3) implied that the use of biofuels was mainly a mixture of softwood and crop straws.Compared with ratios in summer, LG/MN and LG/GA ratios in winter have a narrower range at both IAP and PG sites, and most of them converge at values lower than 10, suggesting wood burning may be dominant in winter, and the contribution of straw burning increases in summer.It has been reported that firewood burning emissions in Beijing represent 47 %-90 % of the total biomass burning and could contribute more than 80 % in winter (Zhou et al., 2017).Although the LG/MN and LG/GA ratios of maize burning and wheat burning are similar and cannot be used to distinguish the two, research into the yield of main crops and farming practices showed that the main straw burning in Beijing and nearby provinces may be associated with maize straw in October to November and wheat straw in May to July (Zhang et al., 2019).
These results suggest that softwood burning and straw burning are the main sources of aerosols from biomass burning in Beijing.Thus, EC from softwood burning and straw burning can be calculated as follows: where f wood represents the fraction of LG from softwood burning, and f straw represents the fraction of LG from straw burning.f straw = 1 − f wood , neglecting other sources of biomass burning.f wood can be expressed as As f wood should be in the range of 0-1, it can be used as the limits of selected EC/OC and OC/LG ratios from softwood and straw.Details of the ratio selection can be found in the Supplement.Once values of f wood are confirmed, OC from softwood (OC wood ) and straw burning (OC straw ) can be obtained by The mass concentrations of OC bb and the contributions of OC wood and OC straw are shown in Fig. 4. The average concentration of OC bb in IAP winter was 2.7 ± 1.3 µg m −3 , with a contribution of 10.6 ± 1.7 % to total OC, which is about half that in PG winter (4.8 ± 2.4 µg m −3 , 10.4 ± 1.5 %).The OC bb concentrations fell in summer (0.6 ± 0.7 µg m −3 at IAP, 2.0 ± 0.8 µg m −3 at PG); the contributions are 6.5 ± 5.2 % and 17.9 ± 3.5 %, respectively.The rural site always has a higher OC bb , which is consistent with our previous discussions.OC wood normally dominated OC bb , while https://doi.org/10.5194/acp-21-8273-2021Atmos.Chem.Phys., 21, 8273-8292, 2021 the increase of the OC straw fraction may be attributed to the local open burning activities.The contributions of OC bb in this study are slightly different from previous studies (Table 3), with lower percentage contributions in the winter sampling period.However, if applying a value of 12.2, a typical OC/LG ratio widely used in previous studies (Andreae and Merlet, 2001;Fu et al., 2012;T. Zhang et al., 2008;Zhang et al., 2007), contributions of OC bb to OC would be 17.7 % and 26.3 % for IAP and PG in winter and 6.1 % and 10.8 % in summer, respectively.Compared with the results in this study, the ratio 12.2 may significantly overestimate the contributions of biomass burning in winter and underestimate the contributions in summer.
In addition, coal combustion may also emit LG, accounting for about 7 % of LG in Beijing (Yan et al., 2018;Y. X. Zhang et al., 2008).Considering that the reported LG/MN and LG/GA ratios of coal combustion (both are in the range of 5-10) are close to those in our measurements, there may be a contribution from coal burning to levoglucosan, but if this amounts to 7 % of levoglucosan, the percentage contribution of biomass burning to OC in winter falls to 10.3 ± 1.6 % and 10.2 ± 1.5 % at IAP and PG, respectively.
Besides the local emissions, regional transport of air masses can elevate the local OC and EC concentrations.To study the influence of regional transport on biomass burning particles, back trajectories and fire spots are plotted in Figs.S1 and S2 in the Supplement.As shown in Fig. S1, air masses at IAP and PG in winter originated from Hebei, Shanxi, and Inner Mongolia, and the fire spots showed low intensity of open burning activities along the air mass transport path during the measurement period.It suggests a less important effect of open burning through regional transport.During the summer, the open burning activities were greatly increased due to the harvesting of wheat (Zhang et al., 2019), which is confirmed by the fire spot distribution in Fig. S2.
Air masses from Hebei, Liaoning, and Shandong provinces, which contain particles from wheat straw burning, may have enhanced the concentrations of OC and EC in Beijing.Considering the open burning of straw is linked with a sudden increase of LG in the ambient atmosphere, the high LG concentrations accompanied with low fire spot intensity suggests a strong local emission, while days with high fire spot intensities may also be affected by regional transport.Combining the analysis of fire spots intensity and the OC wood fractions can help to identify the influence of local emissions and regional transport.Details of regional transport and sources of biomass burning are summarized in Table S2 in the Supplement.

Other sources of OC and consistency with CMB and AMS/ACSM-PMF model results
More detailed source apportionment of OC can be achieved by combining 14 C analysis with primary OC/EC ratios for each source.POC f can be determined from EC f with primary fossil-fuel OC/EC emission ratios (POC/EC) f .However, the (POC/EC) f ratios in previous studies (1.12-2.08 in winter, 0.40-0.77 in summer; Zhang et al., 2017) give much lower POC f values compared to CMB results (Fig. S5 in the Supplement), even though mostly good correlations were found.
In reality, (POC/EC) f ratios vary greatly according to combustion conditions, fuel types and even measurement method for OC and EC (Chow et al., 2001;Han et al., 2016), and it is very hard to accurately predict the (POC/EC) f ratio for a given area.Hence, we used the lowest (OC/EC) f ratios (OC/EC) f,min as the (POC/EC) f to estimate POC f .Due to the limited number of samples for 14 C analysis, the measured lowest (OC/EC) f ratios may be higher than the ratios for the whole sampling period, which will result in an overestimation of POC f .It is necessary to evaluate (OC/EC) f,min ratios for the whole sampling period.The evaluation method is described in the Supplement.In the same way, primary OC from non-fossil sources (POC nf ) can be calculated from EC nf and lowest (OC/EC) nf ratios; therefore concentrations of secondary OC from fossil sources (SOC f ), non-fossil sources (SOC nf ), and OC from cooking (OC ck ) can be obtained by the equations in Table 1.The averaged source apportionment results are presented in Table 4.Primary fossil-derived OC is mainly from coal combustion and traffic emissions in China.However, it cannot be distinguished by 14 C analysis.OC/EC ratios from coal combustion and traffic emissions are dependent on various factors, such as the types of coal, stoves, engines, the vehicle operating modes, and test method.Typical OC/EC ratios of coal combustion and traffic emissions in Beijing are 2.38 ± 0.44 and 0.85 ± 0.16, respectively (Ni et al., 2018).An upper limit of POC from traffic emissions (POC tr ) can be obtained by multiplying EC f by the (OC/EC) tr ratio (0.85 ± 0.16), considering all EC f to come from traffic emissions.A lower limit of POC from coal combustion (POC cc ) is obtained by

Gasoline vehicle
1.3 ± 2.1 0.0 ± 0.0 0.9 ± 1.8 11.4 ± 3.9 1.9 ± 1.5 8.7 ± 5.7 0.1 ± 0.2 0.6 ± 0.3 Industrial CC 5.5 ± 3.9 0.6 subtracting POC tr from POC f .Such calculation shows that POC cc dominated POC at both sites in winter and summer campaigns.The maximum contribution of POC tr to OC was 7.3 % and 5.7 % in winter and 6.8 % and 8.9 % in summer, for IAP and PG, respectively, and POC cc contributed at least 28.5 % and 28.4 % to OC for IAP and PG in winter and 22.2 % and 20.1 % in summer.This is a relatively crude method for source apportionment of primary OC from fossil and non-fossil sources.Ni et al. (2018Ni et al. ( , 2019) ) reported δ 13 C signatures of biomass burning, coal combustion, and traffic emissions, as well as the OC/EC ratios from previous literature.By combining stable carbon isotopic composition analysis of EC with 14 C analysis, the proportions of coal combustion and traffic emission to EC can be derived using Bayesian statistics.The introduction of stable carbon isotopic analysis is suggested as a way to improve our EG method.
Further comparisons with results from application of a CMB method and from application of PMF to ACSM data were conducted to understand the uncertainties in source apportionment from different methods.The source contributions to OC at the IAP and PG sites in winter and summer from the CMB model (Xu et al., 2021;Wu et al., 2020) are presented in Table 4.In brief, seven primary OC sources were apportioned, including emissions from vegetative detritus, biomass burning, cooking, gasoline vehicles, diesel engines, industrial coal combustion, and residential coal combustion, along with other (secondary) OC.Among these sources, coal combustion (the total of residential and industrial coal combustion) accounted for 32.6 % to OC at IAP in winter and 40.0 % to OC at PG in summer, while other OC dominated OC at IAP in summer and at PG in winter, with contributions of 48.2 % and 32.5 %, respectively.
For comparison, OC from gasoline vehicles, diesel engines, industrial coal combustion, and residential coal combustion resolved by the CMB model is summed up as POC f , and OC from vegetative detritus, biomass burning, and cooking is summed up as POC nf .Correlations of different OC sources from the extended Gelencsér method (EG method) and from the CMB model are shown in Fig. 5. Good correlations were found for POC f , POC nf , SOC, and OC bb , despite the combination of sites and seasons (R 2 = 0.96, 0.74, 0.85, and 0.91).The EG method reported lower POC f , POC nf , and OC bb values than those from CMB, with slopes of 0.77, 0.66, and 0.53, respectively.More specifically, OC bb by the EG method is 51 % of that by CMB in winter but 1.33 times higher than CMB in summer (Fig. S6 in the Supplement).The main discrepancy within the apportionment of OC bb is caused by different parameters for the calculations.As the CMB model used source profiles from three major types of cereal straw (wheat, corn, and rice) and two types of wood (pine and mixed wood), it may lead to an overestimation of OC from straw burning.Closer values of POC nf and SOC were found between the two methods in summer, when samples almost all belong to the non-haze period (Figs. S6 and S7 in the Supplement).This indicates that the EG method may perform better when OC concentrations are low.Although poor agreement of OC ck was found between the EG method and the CMB model, the former correlated better with results from the application of PMF to AMS/ACSM data (slope = 0.74, R 2 = 0.61  in the quantification of OC ck .ACSM-PMF may overestimate OC ck by approximately 2 times (Reyes-Villegas, et al., 2018;Yin et al., 2015), whereas CMB may not be sensitive enough to the source profile of cooking aerosols (Abdullahi et al., 2018).Overall, the EG method resolves primary and secondary sources of OC well.Time variations of OC fractions by the EG method and their average contributions are shown in Fig. 6.POC f was the largest contributor to OC at both sites through winter and summer.Comparable contributions of POC f were observed at the urban and rural sites, which reached to 35.8 ± 10.5 % and 34.1 ± 8.7 % in wintertime and fell to 28.9 ± 7.4 % and 29.1 ± 9.4 % in summer, respectively.Pronounced POC f in wintertime implied a significant elevation of coal combustion and traffic emissions.Fossil and non-fossil sources of SOC are distinguished by the EG method in this study for the first time.Average contributions to OC from SOC f are higher in winter.They decreased from 32.0 ± 12.5 % to 25.2 ± 7.6 % at IAP and from 25.2 ± 10.4 % to 21.0 ± 14.4 % at PG from winter to summer.The contributions of SOC nf are slightly greater in summer (18.0 ± 2.9 %, 22.0 ± 17.6 % for IAP in winter and summer and 16.9 ± 10.8 %, 21.7 ± 16.1 % for PG in winter and summer, respectively).Significant contributions of SOC f in the winter sampling period indicated a greater fraction of OC f from ageing and oxidation.The elevated contributions of SOC nf (as a percentage) in the summer sampling period may be assigned to the reduced coal combustion and enhanced biogenic-derived SOC formation.Similar variations of SOC f and OC onf (all OC from non-fossil sources excluding OC bb ) between winter and summer were found in the urban area of Beijing by Zhang et al. (2016).The total SOC accounts for 50.0 ± 12.3 % and 42.0 ± 11.7 % of OC for IAP and PG site in winter, demonstrating the important role of secondary formation processes, especially at the urban site.The average contributions of OC ck were 3.6 ± 2.7 % and 13.4 ± 5.8 % in winter for IAP and PG and 17.4 ± 12.5 % and 10.4 ± 6.7 % in summer, close to those estimated in previous studies (19 ± 4 %; Zhang et al., 2017).The slightly lower value of the OC ck contribution in winter at the IAP site (Fig. 5) was due to OC bb being the overwhelming contributor to POC nf .
In the winter sampling campaign at IAP, POC f was the biggest contributor to OC, followed by SOC f .Both of them were significantly enhanced during haze periods, while the non-fossil fractions, OC bb , OC ck , and SOC nf , did not show much difference between haze and non-haze periods.This implies the haze pollution at IAP in winter was elevated by the accumulation of coal combustion and traffic emissions under favourable weather conditions.The formation of secondary OC associated with coal combustion and traffic emissions was increased during the haze period.In the winter campaign at PG, POC f and SOC f were the top two contributors to OC; however, the contribution of POC f and SOC f did not increase much in the haze period.In contrast, the fractions of SOC nf increased substantially on 3 and 4 December, on which days there were found to be open burning activities in surrounding areas (shown from the fire spots on Fig. S2).This shows that a large proportion of OC bb was transformed to secondary OC during the transport of biomass burning aerosols to the receptor sites.In summer, the sudden increase of the SOC nf fraction on 16 and 17 June at both sites was accompanied by an increase of PM 2.5 and OC concentrations.This is likely due to the open burning activities in surrounding areas.The enhancement of (OC/EC) nf ratios and WSOC nf fractions also suggested secondary OC formation through oxidation of primary non-fossil sources.The standard deviations appear small but obscure the marked differences between seasons.Also, the time series in Fig. 6 shows substantial day-to-day variations in the source contributions within a season but still suggests that meteorological drivers play a major role in determining daily concentrations.

Correlations of WINSOC and WSOC with POC and SOC
In order to better understand the origins and formation mechanism of OC, the correlations between WINSOC and WSOC and between POC and SOC, OC bb and OC ck were investigated (Fig. 7).The regression slopes and correlation coefficients among them are summarized in Table S6.WINSOC f has usually been seen as a proxy for primary fossil-derived OC in many previous studies (Liu et al., 2016c;Miyazaki et al., 2006) S6.
of 1.54, 1.27, 0.99, and 0.99 for IAP and PG in winter and summer (corresponding R 2 is 0.96, 0.89, 0.91, and 0.43).The high WINSOC f / SOC f ratios implied a nonnegligible fraction of WINSOC f in SOC f .Moreover, the ratios of WINSOC f / POC f decreased from winter to summer compared to the WSOC f / POC f ratios increasing, indicating a non-negligible fraction of WSOC f in POC f in summer.WSOC nf and WINSOC nf show good correlations with SOC nf , with larger WINSOC nf / SOC nf ratios in winter.The lower water solubility of SOC f and SOC nf in winter may be due to them originating from the less oxidized semi-volatile POC from wood burning and anthropogenic emission at low temperatures (Favez et al., 2008;Sciare et al., 2011).Weak photochemical activity would also lead to the formation of less oxidized SOC, which is more water-insoluble (Donahue et al., 2006;Robinson et al., 2007).Significant contributions of WINSOC to SOC have been reported in France (Sciare et al., 2011) and Switzerland (Zhang et al., 2016).OC bb correlated well with WSOC nf (R 2 of 0.94, 0.92, 0.93, and 0.93, at IAP and PG in winter and summer) as it is mainly composed of polar and highly oxygenated compounds (Miyazaki et al., 2006).However, more pronounced WINSOC nf in OC bb was found in winter, especially at the rural site.It seems that more water-insoluble fractions were observed in primary OC (POC f , POC nf , OC bb , and OC ck ) at the rural site in both winter and summer.It implies the emitted primary OC at the rural site was probably fresher and hence less aged and oxidized.On the other hand, the source emission profile at rural sites may be different from urban sites, with more heavy-duty diesel trucks with a high content of water-insoluble OC emitted in rural areas.

Conclusions
Measurements of PM 2.5 , OC, EC, and biomass burning tracers were conducted at both urban and rural sites of Beijing from 10 November to 11 December 2016 and from 22 May to 24 June 2017, accompanied by the 14 C analysis of 25 selected samples.On most days, fossil sources dominated EC at IAP and PG in winter and summer, with contributions of 45.9 %-71.7 % at IAP and 48.2 %-76.6 % at PG.The fossil sources of OC contribute 34.7 %-75.0 % and 39.3 %-66.9 % for IAP and PG, with non-fossil fractions of OC elevated in summer.An extended Gelencsér method using the 14 C measurements was applied for the first time to estimate fos-sil and non-fossil sources of primary and secondary OC, as well as OC from biomass burning and cooking (POC f , SOC f , POC nf , SOC nf , OC bb , and OC ck , respectively).Fossilderived POC is a major contributor during winter and summer at both sites.Fossil-derived SOC contributed more in winter, especially at the urban site, with average contributions (to OC) of 32.0 ± 12.5 % and 25.2 ± 7.6 % for IAP and 25.2 ± 10.4 % and 21.0 ± 14.4 % for PG in winter and summer, respectively.The contribution of SOC nf increased in summer, which is probably associated with formation from biogenic emissions.A study of relationships among levoglucosan, mannosan, and galactosan showed that biomass burning was mainly from softwood combustion and straw burning.The extended Gelencsér method using 14 C data provided a more robust calculation of OC bb .The contributions of OC bb to OC were 10.6 ± 1.7 % and 10.4 ± 1.5 % for IAP and PG in winter and 6.5 ± 5.2 % and 17.9 ± 3.5 % for IAP and PG in summer.Correlations among WINSOC and WSOC and POC and SOC showed that WINSOC f and WINSOC nf were the main components of POC f and POC nf , respectively.However, large fractions of WINSOC were found in both SOC f and SOC nf , especially at the rural site, and the contributions of water-insoluble OC decreased from winter to summer, with more WSOC formed under favourable conditions in summer.Although derived from a limited number of samples, our study reflected the different formation mechanisms of SOC between winter and summer and between the urban and rural area.It also confirms the feasibility of a new approach of direct source apportionment of carbonaceous aerosol, which was found to compare generally well with the commonly used chemical mass balance and AMS/ACSM-PMF methods.
Author contributions.ZS and RMH conceived the research.TVV and DL conducted the aerosol sampling and laboratory-based OC/EC analyses.DL and LL carried out the GC-MS analysis.PF supervised GC-MS laboratory work.AV, VM, and GS carried out the 14 C analysis.SS and ASHP supervised the 14 C analysis.XW and JX conducted the CMB modelling at PG and IAP sites, respectively.SH conducted the analysis of the extended Gelencsér (EG) method incorporating 14 C data.YS provided the AMS-PMF data.SH and DL drafted the paper.DS conducted the PMF modelling.
Competing interests.The authors declare that they have no conflict of interest.Special issue statement.This article is part of the special issue "In-depth study of air pollution sources and processes within Beijing and its surrounding region (APHH-Beijing) (ACP/AMT interjournal SI)".It is not associated with a conference.

Figure 1 .
Figure 1.Time series of PM 2.5 and its major components at IAP and PG during winter (left) and summer (right).

Figure 2 .
Figure 2. Time series of concentrations of WSOC f , WSOC nf , WINSOC f , WINSOC nf , EC f , and EC nf (left) and their relative contributions to TC (right) during winter in 2016 and summer in 2017 at IAP and PG.

Figure 3 .
Figure 3. Scatter plot of LG/MN vs.LG/GA from different types of biomass burning emissions(Cheng et al., 2013; Sun et al., 2019a, b;   LG, MN, and GA refer to levoglucosan, mannosan, and galactosan, respectively), including those measured in PM 2.5 samples at IAP and PG during winter and summer.The range of LG/MN is 7.66-13.41and 5.48-22.40for IAP in winter and summer and 2.45-60.20 and 3.27-17.71for PG in winter and summer, while the range of LG/GA is 3.98-6.92,2.78-101.43,1.97-34.66,and 1.51-13.39,respectively.

Figure 4 .
Figure 4. Concentrations of OC from softwood (OC wood ) and OC from straw (OC straw ) at IAP (a) and PG (b) and variations of OC wood fractions.

Figure 5 .
Figure5.Correlations of OC sources from the extended Gelencsér method with those from the CMB model.EG denotes the extended Gelencsér method, (a) primary OC from fossil sources, (b) primary OC from non-fossil sources, (c) secondary OC, (d) OC from biomass burning, (e) OC from cooking, (f) correlations of OC ck from extended Gelencsér method and AMS/ACSM-PMF model (AMS for IAP and ACSM for PG).Good correlations were found within the apportionments of POC f , POC nf , SOC, and OC bb between the CMB and EG method and large discrepancies of OC ck .The PMF model did not show certain source apportionment results from fossil, non-fossil, or secondary OC; i.e.only OC ck concentrations were compared.

Figure 6 .
Figure 6.Time variations of OC source apportionment results by the extended Gelencsér method (upper) and the fractions of each source (i.e., POC f , SOC f , OC bb , OC ck , and SOC nf ) in OC based on the extended Gelencsér method (lower).f: fossil fuel sources, nf: non-fossil sources, bb: biomass burning, ck: cooking.The box denotes the 25th (lower line), 50th (middle line), and 75th (top line) percentiles; the solid squares within the box denote the mean values; the end of the vertical bars represents the 10th (below the box) and 90th (above the box) percentiles; and the solid dots denote maximum and minimum values.

Figure 7 .
Figure 7. Correlations of WINSOC and WSOC with POC and SOC at IAP and PG sites in winter and summer.The slopes and correlation coefficients are summarized in TableS6.

Table 2 .
Statistical summary of concentrations, ratios, and meteorological parameters at IAP and PG sites during winter and summer campaigns.

Table 3 .
Comparison of LG concentrations and OC from biomass burning in this study and related literature.

Table 4 .
Source contribution estimates (µg m −3 ) for OC at IAP and PG in winter and summer.