High contributions of fossil sources to more volatile organic aerosol

Abstract

Sources of particulate organic carbon (OC) with different volatility have rarely been investigated, despite the significant importance for better understanding of the atmospheric processes of organic aerosols. In this study we develop a radiocarbon-based (${}^{\mathrm{14}}\mathrm{C}$) approach for source apportionment of more volatile OC (mvOC) and apply to ambient aerosol samples collected in winter in six Chinese megacities. mvOC is isolated by desorbing organic carbon from the filter samples in helium (He) at 200 ${}^{\circ }$C in a custom-made aerosol combustion system for ${}^{\mathrm{14}}\mathrm{C}$ analysis. Evaluation of this new isolation method shows that the isolated mvOC amount agrees very well with the OC1 fraction (also desorbed at 200 ${}^{\circ }$C in He) measured by a thermal–optical analyzer using the EUSAAR_2 protocol. The mvOC, OC and elemental carbon (EC) of 13 combined PM${}_{\mathrm{2.5}}$ samples in six Chinese cities are analyzed for ${}^{\mathrm{14}}\mathrm{C}$ to investigate their sources and formation mechanisms. The relative contribution of fossil sources to mvOC is $\mathrm{59}\pm \mathrm{11}$ %, consistently larger than the contribution to OC ($\mathrm{48}\pm \mathrm{16}$ %) and smaller than that to EC ($\mathrm{73}\pm \mathrm{9}$ %), despite large differences in fossil contributions in different cities. The average difference in the fossil fractions between mvOC and OC is 13 % (range of 7 %–25 %), similar to that between mvOC and EC (13 %, with a range 4 %–25 %). Secondary OC (SOC) concentrations and sources are modeled based on the${}^{\mathrm{14}}\mathrm{C}$-apportioned OC and EC and compared with concentrations and sources of mvOC. SOC concentrations ($\mathrm{15.4}\pm \mathrm{9.0}$ $\mathrm{\mu }$g m${}^{-\mathrm{3}}$) are consistently higher than those of mvOC ($\mathrm{3.3}\pm \mathrm{2.2}$ $\mathrm{\mu }$g m${}^{-\mathrm{3}}$), indicating that only a fraction of SOC is accounted for by the more volatile carbon fraction desorbed at 200 ${}^{\circ }$C. The fossil fraction in SOC is 43 % (10 %–70 %), lower than that in mvOC (59 %, with a range of 45 %–78 %). Correlation between mvOC and SOC from nonfossil sources (mvOC${}_{\mathrm{nf}}$ vs. SOC${}_{\mathrm{nf}}$) and from fossil sources (mvOC${}_{\mathrm{fossil}}$ vs. SOC${}_{\mathrm{fossil}}$) is examined to further explore sources and formation processes of mvOC and SOC.

1 Introduction

Aerosol particles are of importance for atmospheric chemistry and physics and exert crucial influence on the climate system, air quality and human health (Fuzzi et al., 2015). The carbonaceous fraction of aerosol comprises a large fraction of the fine aerosol mass (20 %–80 %) (Cao et al., 2007; Tao et al., 2017). Carbonaceous aerosols are often operationally subdivided into organic carbon (OC) and elemental carbon (EC). OC consists of a large variety of organic species that cover a great range of volatilities and are not or only weakly light absorbing, while EC is nonvolatile, resistant to chemical transformation and strongly light absorbing (Pöschl, 2005). It should be noted that since no clear distinction between OC and EC exists, OC and EC are operationally defined based on measurement techniques (Petzold et al., 2013). When a thermal–optical method is used to separate and determine OC and EC, EC is described as a thermally refractory carbon continuum and OC is weakly refractory (Pöschl, 2003, 2005; Petzold et al., 2013).

EC is exclusively emitted as primary aerosols from incomplete combustion of biomass (e.g., wood, crop residues and grass) and fossil fuels (e.g., coal, gasoline and diesel). OC can be of both primary and secondary origin. Primary OC (POC) is emitted directly from nonfossil (e.g., biomass burning, biogenic emissions and cooking) and fossil sources. Secondary OC (SOC) is formed in the atmosphere via atmospheric oxidation of volatile organic compounds (VOCs) (Jacobson et al., 2000; Hallquist et al., 2009; Bond et al., 2013). High concentrations of aerosol particles have been observed in many cities, especially in China, where carbonaceous aerosols dominate particle mass concentrations, with SOC responsible for a large fraction of OC (Guo et al., 2014; Huang et al., 2014; Fang et al., 2017). Thus, understanding the sources, formation and transformation of EC, OC and SOC is important for better understanding of air pollution.

For aerosol source apportionment, radiocarbon (${}^{\mathrm{14}}\mathrm{C}$) analysis of OC and EC has been used to quantify their fossil and nonfossil fractions, based on the fact that emissions from fossil sources are ${}^{\mathrm{14}}\mathrm{C}$-free due to long-time decay, whereas nonfossil emissions contain the contemporary ${}^{\mathrm{14}}\mathrm{C}$ content (e.g., Szidat et al., 2004, 2006; Dusek et al., 2013, 2017). Earlier ${}^{\mathrm{14}}\mathrm{C}$ measurements of OC and EC found that fossil sources contribute more to EC than OC (Heal, 2014 and references therein). However, direct ${}^{\mathrm{14}}\mathrm{C}$ measurements of SOC are not possible yet, due to the technical and conceptual difficulties in isolating the SOC fraction from filter samples. In positive-matrix-factorization-based (PMF) source apportionment of an aerosol mass spectrometer (Jayne et al., 2000) dataset, oxidized organic aerosols (OOA, also referred to as secondary organic aerosol, SOA) can be separated into semi-volatile OOA (SV-OOA) and low-volatility OOA (LV-OOA) based on their volatility (e.g., Huffman et al., 2009a; Wang et al., 2017). Volatility is an important physical property as it determines the partitioning between gas and particulate phases of organic species (Donahue et al., 2006, 2009). Direct measurement of OA volatility is challenging because OA is a mixture of thousands of individual organics spanning a wide range of volatilities (Donahue et al., 2011). Different approaches have been used to estimate the OA volatility. For example, using a thermodenuder (TD) coupled with an aerosol mass spectrometer (Burtscher et al., 2001; Wehner et al., 2002), the volatility of different OA components, such as hydrocarbon-like OA (HOA), biomass burning OA (BBOA) and OOA, can be estimated for ambient aerosols. Source and ambient studies indicate that BBOA and HOA are generally more volatile than OOA. Meanwhile, the volatility of BBOA can be quite variable, depending on type of biomass and the combustion conditions, and either higher or lower than that of HOA (e.g., Grieshop et al., 2009c; Huffman et al. 2009a; Paciga et al., 2016; Cao et al., 2018). For example, Huffman et al. (2009b) measured the volatility of primary OA from biomass burning and found that the more volatile BBOA was generally dominated by smoldering combustion, while the less volatile BBOA was more influenced by flaming combustion. Paciga et al. (2016) found that the BBOA was more volatile than HOA in an urban background site in Paris, with 50 % semi-volatile organic components (SVOCs with effective saturation concentrations $C^{*}$ of 1–100 $\mathrm{\mu }$g m${}^{-\mathrm{3}}$ range) and 30 % low-volatility organic compounds (LVOCs with $C^{*}$ of 10${}^{-\mathrm{3}}$–0.1 $\mathrm{\mu }$g m${}^{-\mathrm{3}}$). Gaston et al. (2016) found that levoglucosan (i.e., the biomass burning tracer) was semi-volatile and had the potential to form low-volatility products during the chemical aging of ambient biomass burning aerosol, applying high-resolution time-of-flight chemical ionization mass spectrometer (HRTOF-CIMS) coupled to a Filter Inlet for Gases and AEROsols (FIGAERO) instrument for aerosol in the northwestern US.

In some recent studies, a thermal desorption approach has been used for analysis of filter samples at lower and higher temperatures as an indicator of volatility (Holzinger et al., 2013; Lopez-Hilfiker et al., 2014, 2015; Timkovsky et al., 2015; Masalaite et al., 2017, 2018). Even though the desorption temperature is not a direct measure of the particle volatility, OC desorbed from the filter at lower temperature tends to be more volatile, whereas the less volatile (i.e., more refractory) OC tends to be desorbed at higher temperatures. Vodička et al. (2015) and Masalaite et al. (2017) found that OC in urban samples was desorbed at lower temperatures compared to coastal, forest and rural background samples in Eastern Europe. Both studies suggest that OC is more volatile in urban areas close to the emission sources. Keller and Burtscher (2017) found that after aging of biomass burning emissions, less OC desorbed at lower temperatures. The reduction in volatility after aging agrees with previous TD aerosol mass spectrometer studies (Grieshop et al., 2009a, b). Even though the thermal desorption and TD measurements are different methods, results of both methods are usually in qualitative agreement and show similar trends. Earlier studies (e.g., Grieshop et al., 2009c; Chen et al., 2007) show similar trends qualitatively. Ma et al. (2016) made first efforts to link the OC fractions desorbed at different temperatures from filter samples on a thermal-optical analyzer to the volatility basis set, showing that OC fractions desorbed at lower temperatures (e.g., 140 and 280 ${}^{\circ }$C) are semi-volatile.

Both POC and SOC can contribute to the OC fraction desorbed at lower temperatures. A number of previous studies have found that at least a part of POC from various sources is semi-volatile, including wood burning, gasoline and diesel vehicle exhausts and cooking (Lipsky and Robinson, 2006; Shrivastava et al., 2006; Robinson et al., 2007; Grieshop et al., 2009b, c; May et al., 2013a, b). Moreover, a significant part of freshly formed SOC is semi-volatile and will contribute to the more volatile OC fraction desorbed at lower temperatures (Holzinger et al., 2010; Salo et al., 2011; Meusinger et al., 2017; Gkatzelis et al., 2018). The relative contributions from fossil and nonfossil sources to OC can be dependent on the volatility. For example, recent studies during winter in Lithuania have shown that vehicular sources are associated with OC fractions desorbed at lower temperatures, while biomass burning (as a nonfossil source) contributes more strongly to less volatile OC fractions desorbed at higher temperatures (Masalaite et al., 2017, 2018).

Although an increasing number of studies has shown that OC from different emission sources and/or at different aging status may have different volatility, a quantitative study of the sources of OC with different volatility is scarce. In this study, wintertime fine particulate matter (PM${}_{\mathrm{2.5}}$, particles with aerodynamic diameter $\mathit{<}\mathrm{2.5}$ $\mathrm{\mu }$m) samples from six Chinese megacities are studied. Sources of carbonaceous aerosols, including different carbon fractions such as OC, EC and SOC, are estimated by ${}^{\mathrm{14}}\mathrm{C}$ source apportionment. In addition, we operationally define a more volatile OC fraction (mvOC) as the carbon fraction desorbed from the filter samples at 200 ${}^{\circ }$C in helium (He) and investigate sources of mvOC. The source contribution to mvOC is measured unambiguously by ${}^{\mathrm{14}}\mathrm{C}$.

2 Methods

2.1 Sampling

Simultaneous sampling was made during winter 2013/2014 in two northern (Beijing and Taiyuan) and three southern (Chongqing, Guangzhou and Shanghai) Chinese cities (Supplement Fig. S1). The samples from Xi'an in northern China were collected during winter 2015/2016. These sites were selected to represent urban-scale concentrations and located in the university or research center campus, $\mathit{>}\mathrm{100}$ m from local sources such as main roadways (Table S1 in the Supplement). The 24 h integrated PM${}_{\mathrm{2.5}}$ samples were collected from 10:00 to 10:00 the next day (local standard time, LST). Filter samplers were deployed on rooftops about 6 to 20 m above ground level. In Xi'an, Beijing, Guangzhou and Shanghai, PM${}_{\mathrm{2.5}}$ samples were collected on pre-baked (780 ${}^{\circ }$C for 3 h) quartz fiber filter (QM-A, Whatman Inc., Clifton, NJ, USA, 20.3 cm $\times \mathrm{25.4}$ cm) using a high-volume sampler (TE-6070 MFC, Tisch Inc., Cleveland, OH, USA) at a flow rate of 1.0 m${}^{\mathrm{3}}$ min${}^{-\mathrm{1}}$. In Taiyuan and Chongqing, samples were collected on pre-baked 47 mm Whatman quartz microfiber filters (QM-A) using mini-volume samplers at a flow rate of 5 L min${}^{-\mathrm{1}}$ (Airmetrics, OR, USA). After sampling, all filters were packed in a pre-baked aluminum foil (450 ${}^{\circ }$C for 3 h), sealed in polyethylene bags and stored in a freezer at $-\mathrm{18}$ ${}^{\circ }$C until analysis.

2.2 Determination of carbon fractions by thermal-optical analysis

For the PM${}_{\mathrm{2.5}}$ samples collected in Xi'an, filter pieces of 1.5 cm${}^{\mathrm{2}}$ were taken for OC and EC analysis using a carbon analyzer (Model 5L, Sunset Laboratory, Inc., Portland, OR, USA) following the thermal-optical transmittance protocol EUSAAR_2 (European Supersites for Atmospheric Aerosol Research; Cavalli et al., 2010). The EUSAAR_2 protocol defines OC1 as the carbon fraction that desorbs in helium (He) at 200 ${}^{\circ }$C for 2 min. OC1 is compared with more volatile OC (mvOC) also extracted in He at 200 ${}^{\circ }$C but for 5 min using our aerosol combustion system (ACS) (Dusek et al., 2014). No charred OC is observed by the transmittance signal in the OC1 stage during the thermal–optical analysis. For PM${}_{\mathrm{2.5}}$ samples collected in Beijing, Taiyuan, Chongqing, Guangzhou and Shanghai, filter pieces of 0.5 cm${}^{\mathrm{2}}$ were used to measure OC and EC using a Desert Research Institute (DRI) Model 2001 Thermal/Optical Carbon Analyzer (Atmoslytic Inc., Calabasas, CA, USA), following the IMPROVE_A thermal–optical reflectance (TOR) protocol (Chow et al., 2007). OC fractions of EUSAAR_2 protocol are desorbed in He at 200 ${}^{\circ }$C (OC1), 300 ${}^{\circ }$C (OC2), 450 ${}^{\circ }$C (OC3) and 650 ${}^{\circ }$C (OC4). Different from EUSAAR_2 protocol, IMPROVE_A (Interagency Monitoring of Protected Visual Environments) protocol defines OC fractions at stepwise temperature of 140, 280, 480 and 580 ${}^{\circ }$C in He. Details of the carbon fraction measurements were described in our previous work (Cao et al., 2013; Zenker et al., 2017).

2.3

Radiocarbon (${}^{\mathrm{14}}\mathrm{C}$) measurements of mvOC, OC and EC

2.3.1

Sample selection for ${}^{\mathrm{14}}\mathrm{C}$ analysis

For ${}^{\mathrm{14}}\mathrm{C}$ analysis, we selected the samples carefully to cover periods of low and high total carbon (TC, the sum of OC and EC) and PM${}_{\mathrm{2.5}}$ concentrations to get samples representative of various pollution conditions that did occur in each city. Two or three composite samples from each city representing high (H), medium (M) and low (L) TC concentrations were selected for ${}^{\mathrm{14}}\mathrm{C}$ analysis of EC, OC and mvOC (Fig. S2). Each composite sample consists of two to four 24 h filter pieces with similar TC loadings and air mass backward trajectories (Fig. S2, Table S1).

2.3.2 Extraction of mvOC, OC and EC

Three separate extractions were performed for mvOC, OC and EC on our aerosol combustion system (ACS) (Dusek et al., 2014). The ACS consists of a combustion tube and a purification line. Aerosol filter pieces placed on the filter boat are heated at different temperatures in pure He or ${\mathrm{O}}_{\mathrm{2}}$ and oxidized through the platinum catalyst in the combustion tube. The resulting ${\mathrm{CO}}_{\mathrm{2}}$ is isolated and separated from other gases (e.g., ${\mathrm{NO}}_x$, water vapor) in the purification line. The purification line is equipped with an oven filled with copper grains and silver wire heated at 650 ${}^{\circ }$C to remove ${\mathrm{NO}}_x$ and liberated halogen, a U-type tube cooled with a dry ice–ethanol mixture for water removal, and a flask containing phosphorous pentoxide (${\mathrm{P}}_{\mathrm{2}}{\mathrm{O}}_{\mathrm{5}}$) for removal of any trace water. The amount of purified ${\mathrm{CO}}_{\mathrm{2}}$ is manometrically quantified and subsequently stored in flame-sealed glass ampoules.

mvOC is desorbed by heating the filter pieces at 200 ${}^{\circ }$C in He for 5 min. After introducing the filter, the ovens are flushed with He for 10 min. Pilot tests show that the flushing time (10, 15 or 60 min) before heating the filter pieces does not affect the desorbed amount of mvOC (Table S2). OC is combusted by heating filter pieces at 375 ${}^{\circ }$C in pure ${\mathrm{O}}_{\mathrm{2}}$ for 10 min. EC is extracted from a separate filter piece after removing OC completely. First water-soluble OC is removed from the filter through water extraction (Dusek et al., 2014) to minimize the charring of organic materials (Yu et al., 2002; Zhang et al., 2012). Then, by heating the filter pieces at 375 ${}^{\circ }$C in ${\mathrm{O}}_{\mathrm{2}}$ for 10 min, most water-insoluble OC can be removed. Subsequently, the oven temperature is increased to 450 ${}^{\circ }$C for 3 min to remove most refractory OC that is left on the filter. However, during this step some less refractory EC might be lost. Finally, the remaining EC is combusted by heating at 650 ${}^{\circ }$C in ${\mathrm{O}}_{\mathrm{2}}$ for 5 min (Dusek et al., 2014, 2017; Zenker et al., 2017).

Contamination introduced by the isolation procedure is determined by following exactly the same isolating procedures with either an empty filter boat or with preheated filters (at 650 ${}^{\circ }$C in ${\mathrm{O}}_{\mathrm{2}}$ for 10 min). The contamination introduced by the combustion process yields $\mathrm{0.72}\pm \mathrm{0.44}$ $\mathrm{\mu }$gC EC, $\mathrm{0.85}\pm \mathrm{0.49}$ $\mathrm{\mu }$gC OC and $\mathrm{0.52}\pm \mathrm{0.31}$ $\mathrm{\mu }$gC mvOC per extraction. Compared with our sample size of 30–391 $\mathrm{\mu }$gC EC, 30–445 $\mathrm{\mu }$gC OC and 15–121 $\mathrm{\mu }$gC mvOC, the blanks are $\le \mathrm{7}$ % of the mvOC sample amount and $\mathit{<}\mathrm{3}$ % of OC and EC sample amount and thus relatively small compared to our sample sizes. The mvOC concentrations reported in this study are corrected for contamination by subtraction. OC and EC concentrations are measured following thermal-optical protocols using carbon analyzers (Sect. 2.2), thus are not affected by the isolation procedure using the ACS. For ${}^{\mathrm{14}}\mathrm{C}$ values (Eq. 1), the contamination is assessed but not used for further ${}^{\mathrm{14}}\mathrm{C}$ data correction for mvOC, OC and EC because the corrections for the small blanks will not introduce large uncertainties to the data, as explained in Supplement Sect. S3.

Two standards with known ${}^{\mathrm{14}}\mathrm{C}{/}^{\mathrm{12}}\mathrm{C}$ ratios are combusted using the ACS as quality control for the combustion process: an oxalic acid standard and anthracite. Small amounts of solid standard materials are directly put on the filter boat of the combustion tube and heated in ${\mathrm{O}}_{\mathrm{2}}$ at 650 ${}^{\circ }$C for 10 min. In the further ${}^{\mathrm{14}}\mathrm{C}$ analysis, the ${\mathrm{CO}}_{\mathrm{2}}$ derived from combustion of the standards is treated exactly like the samples. Therefore, the contamination introduced by the combustion process can be estimated from the deviation of measured values from the nominal values of the standards (Table S3) (Dusek et al., 2014). The contamination inferred in this indirect way is below 1.5 $\mathrm{\mu }$gC per extraction, which is slightly higher than the directly measured contamination of OC and EC separately but in the range of a TC contamination. It is also relatively small compared to the size of OC (30–445 $\mathrm{\mu }$gC) and EC samples (30–391 $\mathrm{\mu }$gC) in this study and thus can be neglected (Sect. S3). The ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$ (Eq. 1) of the contamination can also be indirectly inferred and varied between 0.2 and 0.6 for standards measured on the ACS over the last 2 years. This is broadly within the range of our samples and implies that a correction for ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$ will even be considerably smaller than 3 %.

2.3.3

${}^{\mathrm{14}}\mathrm{C}$ measurements by accelerator mass spectrometer (AMS)

AMS measurements were conducted at the Centre for Isotope Research (CIO) at the University of Groningen. The extracted ${\mathrm{CO}}_{\mathrm{2}}$ is released from the glass ampules and directed to a gas inlet system (Ruff et al., 2007), where the sample is diluted with He to 5 % ${\mathrm{CO}}_{\mathrm{2}}$ (Salazar et al., 2015). The ${\mathrm{CO}}_{\mathrm{2}}/\mathrm{He}$ mixture is directly fed into the Cs sputter ion sources of the Mini Carbon Dating System (MICADAS) AMS at a constant rate (Synal et al., 2007).

The ${}^{\mathrm{14}}\mathrm{C}{/}^{\mathrm{12}}\mathrm{C}$ ratio of an aerosol sample is usually reported relative to an oxalic acid standard (OXII) and expressed as fraction modern (${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$). The ${}^{\mathrm{14}}\mathrm{C}{/}^{\mathrm{12}}\mathrm{C}$ ratio of the standard is related to the unperturbed atmosphere in the reference year of 1950 by multiplying it by a factor of 0.7459 (Mook and Van Der Plicht, 1999; Reimer et al., 2004): 1$${\mathrm{F}}^{\mathrm{14}}\mathrm{C}={\displaystyle \frac{{(}^{\mathrm{14}}\mathrm{C}{/}^{\mathrm{12}}\mathrm{C}{)}_{\mathrm{sample},[-\mathrm{25}]}}{\mathrm{0.7459}\times {(}^{\mathrm{14}}\mathrm{C}{/}^{\mathrm{12}}\mathrm{C}{)}_{\mathrm{OXII},[-\mathrm{25}]}}},$$ where the ${}^{\mathrm{14}}\mathrm{C}{/}^{\mathrm{12}}\mathrm{C}$ ratio of the sample and standard are both corrected for machine background and normalized to ${\mathit{\delta }}^{\mathrm{13}}\mathrm{C}$ $=-\mathrm{25}$ ‰. Aerosol carbon from fossil sources has ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$ $=$ 0 due to the extinction of ${}^{\mathrm{14}}\mathrm{C}$ after long-time decay. Aerosol carbon from contemporary (or nonfossil) sources should have ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$ $\sim \mathrm{1}$ in an atmosphere undisturbed by human influences. However, the ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$ values of contemporary carbon are higher than 1 due to the nuclear bomb tests that nearly doubled the ${}^{\mathrm{14}}{\mathrm{CO}}_{\mathrm{2}}$ in the atmosphere in the 1960s. The ${}^{\mathrm{14}}{\mathrm{CO}}_{\mathrm{2}}$ produced by the nuclear bomb tests has been taken up by oceans and the biosphere and diluted by the ${}^{\mathrm{14}}\mathrm{C}$-free ${\mathrm{CO}}_{\mathrm{2}}$ emitted by the fossil fuel burning. Currently, ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$ of the atmosphere ${\mathrm{CO}}_{\mathrm{2}}$ is $\sim \mathrm{1.04}$ (Levin et al., 2010).

The ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$ values are corrected for cross contamination (also known as memory effect) (Wacker et al., 2010) using alternate measurements of OXII and ${}^{\mathrm{14}}\mathrm{C}$-free material as gaseous standards. Correction for instrument background (Salazar et al., 2015) is done by subtracting the memory-corrected ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$ values of the ${}^{\mathrm{14}}\mathrm{C}$-free standard. Finally, the values are normalized to the average value of the (memory- and background-corrected) OXII standards. All standards used for the corrections are measured on the same day as the samples.

2.4 Estimation of source contributions to different carbon fractions

${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$ of EC, OC and mvOC (i.e., ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\left(\mathrm{EC}\right)}$, ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\left(\mathrm{OC}\right)}$ and ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\left(\mathrm{mvOC}\right)}$, respectively) were directly measured. We define “more refractory organic carbon” (mrOC) as the difference between OC and mvOC. ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$ of mrOC (${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\left(\mathrm{mrOC}\right)}$) was calculated by isotope mass balance:

2$${\mathrm{F}}^{\mathrm{14}}{\mathrm{C}}_{\left(\mathrm{mrOC}\right)}={\displaystyle \frac{{\mathrm{F}}^{\mathrm{14}}{\mathrm{C}}_{\left(\mathrm{OC}\right)}\times \mathrm{OC}-{\mathrm{F}}^{\mathrm{14}}{\mathrm{C}}_{\left(\mathrm{mvOC}\right)}\times \mathrm{mvOC}}{\mathrm{OC}-\mathrm{mvOC}}}.$$ ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\left(\mathrm{EC}\right)}$ was converted to the fraction of biomass burning ($f_{\mathrm{bb}}$(EC)) by dividing with an ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$ value representative of typical biomass burning emissions (${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\mathrm{bb}}$). Analogously, the fraction of nonfossil OC, mvOC and mrOC (i.e., $f_{\mathrm{nf}}$(OC), $f_{\mathrm{nf}}$(mvOC) and$f_{\mathrm{nf}}$(mrOC), respectively) were estimated from their corresponding ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$ values and ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$ of nonfossil sources (${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\mathrm{nf}}$). ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\mathrm{bb}}$ and ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\mathrm{nf}}$ are estimated as $\mathrm{1.10}\pm \mathrm{0.05}$ and $\mathrm{1.09}\pm \mathrm{0.05}$ (Lewis et al., 2004; Mohn et al., 2008; Palstra and Meijer, 2014), respectively, based on tree-growth models and the assumption that wood burning dominates biomass burning. ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\mathrm{bb}}$ for EC is slightly bigger than ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\mathrm{nf}}$ for OC because, besides biomass burning, biogenic emissions also contribute to OC but have a smaller ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$ than that of biomass burning. The estimation of ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\mathrm{bb}}$ and ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\mathrm{nf}}$ has been reported in our previous study (Ni et al., 2018).

Then, $f_{\mathrm{bb}}$(EC) was used to determine EC mass concentrations from nonfossil biomass (EC${}_{\mathrm{bb}}$) and fossil fuel combustion (EC${}_{\mathrm{fossil}}$): $$\begin{array}{}\text{3}& {\displaystyle }{\mathrm{EC}}_{\mathrm{bb}}=\mathrm{EC}\times f_{\mathrm{bb}}\left(\mathrm{EC}\right),\text{4}& {\displaystyle }{\mathrm{EC}}_{\mathrm{fossil}}=\mathrm{EC}-{\mathrm{EC}}_{\mathrm{bb}}.\end{array}$$ Analogously, mass concentrations of OC, mvOC and mrOC from nonfossil sources (OC${}_{\mathrm{nf}}$, mvOC${}_{\mathrm{nf}}$ and mrOC${}_{\mathrm{nf}}$, respectively) and fossil sources (OC${}_{\mathrm{fossil}}$, mvOC${}_{\mathrm{fossil}}$ and mrOC${}_{\mathrm{fossil}}$, respectively) were determined.

Secondary OC (SOC) includes SOC from fossil (SOC${}_{\mathrm{fossil}}$) and nonfossil sources (SOC${}_{\mathrm{nf}}$): 5$$\mathrm{SOC}={\mathrm{SOC}}_{\mathrm{fossil}}+{\mathrm{SOC}}_{\mathrm{nf}}.$$ Fraction fossil in total SOC ($f_{\mathrm{fossil}}$(SOC)) was formulated as follows: 6$$f_{\mathrm{fossil}}\left(\mathrm{SOC}\right)={\displaystyle \frac{{\mathrm{SOC}}_{\mathrm{fossil}}}{{\mathrm{SOC}}_{\mathrm{fossil}}+S{\mathrm{OC}}_{\mathrm{nf}}}}.$$ For the following calculations we assume that SOC${}_{\mathrm{nf}}$ can be approximated by OC from nonfossil sources excluding primary biomass burning (OC${}_{\mathrm{o},\mathrm{nf}}$; OC other nonfossil). In principle, OC${}_{\mathrm{o},\mathrm{nf}}$ includes SOC${}_{\mathrm{nf}}$ and nonfossil primary OC from vegetative detritus, bioaerosols, resuspended soil organic matter, or cooking. But the contributions from plant detritus, bioaerosols and spores to PM${}_{\mathrm{2.5}}$ are likely small (Song et al., 2006; Hu et al., 2010; Guo et al., 2012). If cooking sources are significant then this assumption results in an upper limit of SOC${}_{\mathrm{nf}}$. 7$${\mathrm{SOC}}_{\mathrm{nf}}\cong {\mathrm{OC}}_{\mathrm{o},\mathrm{nf}}={\mathrm{OC}}_{\mathrm{nf}}-{\mathrm{POC}}_{\mathrm{bb}},$$ where OC${}_{\mathrm{o},\mathrm{nf}}$ was calculated by the difference between OC${}_{\mathrm{nf}}$ and primary OC from biomass burning (POC${}_{\mathrm{bb}}$). POC${}_{\mathrm{bb}}$ was calculated by multiplying EC${}_{\mathrm{bb}}$ with the primary $\mathrm{OC}/\mathrm{EC}$ ratio of biomass burning ($r_{\mathrm{bb}}$): 8$${\mathrm{POC}}_{\mathrm{bb}}={\mathrm{EC}}_{\mathrm{bb}}\times r_{\mathrm{bb}}.$$ SOC${}_{\mathrm{fossil}}$ was calculated by subtracting primary fossil OC (POC${}_{\mathrm{fossil}}$) from OC${}_{\mathrm{fossil}}$. POC${}_{\mathrm{fossil}}$ was estimated by multiplying EC${}_{\mathrm{fossil}}$ with primary $\mathrm{OC}/\mathrm{EC}$ ratio of fossil fuel combustion ($r_{\mathrm{fossil}}$): $$\begin{array}{}\text{9}& {\displaystyle }{\mathrm{SOC}}_{\mathrm{fossil}}={\mathrm{OC}}_{\mathrm{fossil}}-{\mathrm{POC}}_{\mathrm{fossil}},\text{10}& {\displaystyle }{\mathrm{POC}}_{\mathrm{fossil}}={\mathrm{EC}}_{\mathrm{fossil}}\times r_{\mathrm{fossil}}.\end{array}$$ Fossil sources in China are almost exclusively from coal combustion and vehicle emissions, thus $r_{\mathrm{fossil}}$ was estimated as follows: 11$$r_{\mathrm{fossil}}=r_{\mathrm{coal}}\times p+r_{\mathrm{vehicle}}\times (\mathrm{1}-p),$$ where $r_{\mathrm{coal}}$ and $r_{\mathrm{vehicle}}$ are the primary $\mathrm{OC}/\mathrm{EC}$ ratio of coal combustion and vehicle emissions, respectively. The $r_{\mathrm{bb}}$, $r_{\mathrm{coal}}$ and $r_{\mathrm{vehicle}}$ varies with the fuel types and properties, combustion conditions, sampling and analysis methods, etc. Best estimate of $r_{\mathrm{bb}}$($\mathrm{4}\pm \mathrm{1}$; average $\pm$ standard deviation), $r_{\mathrm{coal}}$($\mathrm{2.38}\pm \mathrm{0.44}$) and $r_{\mathrm{vehicle}}$($\mathrm{0.85}\pm \mathrm{0.16}$) is done through a literature search and described in our earlier study (Ni et al., 2018). $p$ is the fraction of EC from coal combustion (EC${}_{\mathrm{coal}}$) in EC${}_{\mathrm{fossil}}$.

We used two different methods to estimate $p$. (1) Since both coal combustion and vehicle emissions do not contain ${}^{\mathrm{14}}\mathrm{C}$, they can not be distinguished by ${}^{\mathrm{14}}\mathrm{C}$ measurements alone. Therefore, $p$ is randomly chosen from 0 to 1, that is no constraint on $p$ values. (2) EC from coal combustion is on average more enriched in the stable carbon isotope ${}^{\mathrm{13}}$C compared to vehicle emissions. Therefore, complementing ${}^{\mathrm{14}}\mathrm{C}$ results of EC with measurements of the ${}^{\mathrm{13}}\mathrm{C}{/}^{\mathrm{12}}\mathrm{C}$ ratios of EC (expressed as ${\mathit{\delta }}^{\mathrm{13}}\mathrm{C}$${}_{\mathrm{EC}}$ in Eq. S1; Sect. S1) allows separation of EC${}_{\mathrm{fossil}}$ into EC${}_{\mathrm{coal}}$ and EC from vehicle emissions (EC${}_{\mathrm{vehicle}}$). Samples taken from Beijing, Shanghai and Guangzhou using high-volume samplers had enough material for analysis of both ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\mathrm{EC}}$ and ${\mathit{\delta }}^{\mathrm{13}}\mathrm{C}$${}_{\mathrm{EC}}$, thus EC${}_{\mathrm{coal}}$ and EC${}_{\mathrm{vehicle}}$ are separated as described in detail in Sect. S2. In brief, the fraction of coal combustion and vehicle emissions in EC was calculated from measured ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\mathrm{EC}}$ and ${\mathit{\delta }}^{\mathrm{13}}\mathrm{C}$${}_{\mathrm{EC}}$ for ambient EC combined with the source signatures. Bayesian Markov chain Monte Carlo (MCMC) calculations were used to account for the uncertainties in the source signatures and the measurement uncertainties (Andersson, 2011; Andersson et al., 2015). The results of the Bayesian calculations are the posterior probability density functions (PDFs) for the relative contributions of each source to EC ($f_{\mathrm{coal}}$ and $f_{\mathrm{vehicle}}$; Fig. S3). The $p$ was formulated as follows: 12$$p={\displaystyle \frac{f_{\mathrm{coal}}}{f_{\mathrm{fossil}}}}={\displaystyle \frac{f_{\mathrm{coal}}}{f_{\mathrm{coal}}+f_{\mathrm{vehicle}}}}.$$ The PDF of $p$ was derived from the PDF of $f_{\mathrm{coal}}$ and $f_{\mathrm{vehicle}}$ and shown in Fig. S4.

To propagate uncertainties, a Monte Carlo simulation with 10 000 individual calculations was conducted. For each individual calculation, ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\left(\mathrm{EC}\right)}$, ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\left(\mathrm{OC}\right)}$, ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\left(\mathrm{mvOC}\right)}$, and concentrations of EC, OC and mvOC were randomly chosen from a normal distribution that is symmetric around the measured values with the experimental uncertainties as standard deviation (SD; Table S4). For ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\mathrm{bb}}$, ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\mathrm{nf}}$, $r_{\mathrm{bb}}$, $r_{\mathrm{coal}}$ and $r_{\mathrm{vehicle}}$, random values were chosen from a triangular frequency distribution, which has its maximum at the central value and at 0 at the lower limit and upper limit. For $p$ ranging from 0 to 1 (no ${}^{\mathrm{13}}\mathrm{C}$ constraints), $p$ was randomly chosen from a uniform distribution. For $p$ constrained by ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$ and ${\mathit{\delta }}^{\mathrm{13}}\mathrm{C}$ using MCMC (hereafter, ${}^{\mathrm{13}}$C-constrained $p$), random values from the respective PDF of $p$ were used (Fig. S4). In this way 10 000 different estimation of $f_{\mathrm{bb}}$(EC),$f_{\mathrm{nf}}$(OC), $f_{\mathrm{nf}}$(mvOC), $f_{\mathrm{nf}}$(mrOC), EC${}_{\mathrm{bb}}$, EC${}_{\mathrm{fossil}}$, OC${}_{\mathrm{nf}}$, OC${}_{\mathrm{fossil}}$, mvOC${}_{\mathrm{nf}}$, mvOC${}_{\mathrm{fossil}}$, mrOC${}_{\mathrm{nf}}$, mrOC${}_{\mathrm{fossil}}$, SOC${}_{\mathrm{nf}}$, SOC${}_{\mathrm{fossil}}$, SOC and $f_{\mathrm{fossil}}$(SOC) were calculated (Tables S5, S6, S7, S8). The derived average represents the best estimate and the SD represents the combined uncertainties.

3 Results

3.1 Method evaluation and quality control for mvOC extraction

The separation of OC and EC for ${}^{\mathrm{14}}\mathrm{C}$ analysis using our aerosol combustion system (ACS) were thoroughly evaluated by Dusek et al. (2014, 2017). It is thus necessary to validate the new extraction method for isolating mvOC. The reproducibility of the extracted mvOC amount was tested for two independent test filters with mvOC loadings of 6 and 18 $\mathrm{\mu }$g cm${}^{-\mathrm{2}}$, respectively (Fig. S5). The coefficient of variation was determined as a measure of reproducibility. The reproducibility was found to be $\sim \mathrm{5}$ % ($n=\mathrm{9}$).

Figure 1

(a) OC1 loading ($\mathrm{\mu }$g cm${}^{-\mathrm{2}}$) measured with a thermal–optical method (EUSAAR_2 protocol using a Sunset analyzer) and mvOC loading measured using our custom-made aerosol combustion system. The outliers are shown in red circle and square with their sample name on the top $x$ axis. (b) ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\left(\mathrm{mvOC}\right)}$ of winter-M and winter-H with high and low recoveries with respect to OC1. The red numbers in brackets in panels (a) and (b) give the respective outliers' mvOC recoveries to OC1.

[Figure omitted. See PDF]

Since carbon fractions (e.g., OC1) at different desorption temperatures have mostly been measured using the EUSAAR_2 protocol in many previous studies (e.g., Vodička et al., 2015; Keller and Burtscher, 2017), our goal is to define the mvOC fraction as representative of OC1. Therefore, the mvOC is desorbed at 200 ${}^{\circ }$C, the same temperature as used for OC1 in the EUSAAR_2 protocol. However, the extracted amounts of the ACS system might differ due to different heating rates and length of the temperature step. The winter samples from Xi'an as well as the two test filters described above are used to compare mvOC concentrations from the ACS system to OC1. For most samples, excellent agreement was found between mvOC and OC1 (Fig. 1a) and most data points fall close to the $\mathrm{1}:\mathrm{1}$ line. However, there are three data points deviating largely from the $\mathrm{1}:\mathrm{1}$ line (red circle and square in Fig. 1). The two red squares represent the mvOC extraction for sample winter-H and winter-M using larger filter pieces (i.e., more mvOC in $\mathrm{\mu }$gC per extraction). With larger filter pieces, the mvOC ($\mathrm{\mu }$gC cm${}^{-\mathrm{2}}$ filter area) extracted by the ACS system is significantly lower than that measured by the Sunset analyzer. The recoveries of these two outliers are 0.59 for sample winter-H and 0.74 for winter-M, calculated by dividing mvOC mass by OC1 mass. For winter-H, the low recovery is also repeatable with the same filter area (Fig. S5). The low recoveries for large filter pieces may result from the lower temperature ($\mathit{<}\mathrm{200}$ ${}^{\circ }$C) towards the ends of the filter boat. At the relatively low temperature of 200 ${}^{\circ }$C, only a 3.5 cm length at the center of a 12 cm combustion tube is maintained at 200 ${}^{\circ }$C and outside of that 3.5 cm the temperature is lower than 200 ${}^{\circ }$C, e.g., $\sim \mathrm{170}$ ${}^{\circ }$C at the end of the combustion tube as measured by a thermocouple. When filter pieces are large and placed outside the centered 3.5 cm, the desorption temperature for part of the filter pieces will lower than 200 ${}^{\circ }$C, leading to lower desorbed mvOC amount. Another possibility is saturations of the catalyst (platinum, Pt) in the combustion line of ACS system. Pt works as catalyst by collecting oxygen atoms on the surface as has been demonstrated by direct observation of an ultra-thin oxygen layer on the Pt surface by a microscope (Spronsen et al., 2017). This is used to oxidize CO and hydrocarbons to ${\mathrm{CO}}_{\mathrm{2}}$ in a reducing atmosphere. Thus, for large sample amounts it is possible that the oxygen on the catalyst could not be sufficient to oxidize all desorbed CO and hydrocarbons to ${\mathrm{CO}}_{\mathrm{2}}$. However, we observed recoveries near 100 % for mvOC amounts up to 120 $\mathrm{\mu }$g, which was higher than the total amount of desorbed mvOC for the winter-H sample. Therefore, limited catalyst capacity is not the likely explanation for the low recoveries. For subsequent experiments we consequently placed the filter pieces carefully in the 3.5 cm, 200 ${}^{\circ }$C section of the combustion tube and avoided stacking multiple filter pieces to ensure a desorption temperature of 200 ${}^{\circ }$C and sufficient He supply.

To examine the effect of the low recoveries of mvOC on the ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\left(\mathrm{mvOC}\right)}$, we compare the ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\left(\mathrm{mvOC}\right)}$ of winter-M and winter-H samples for both high and low mvOC recoveries. ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\left(\mathrm{mvOC}\right)}$ for low recoveries is roughly 0.05 (absolute value, Fig. 1b) higher than for high recoveries, which is non-negligible compared to the measurement uncertainty. In addition, to validate the measured ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\left(\mathrm{mvOC}\right)}$ for the combined sample winter-M, we also extracted mvOC separately from those three filters that were combined for the composite winter-M sample. Figure 2 shows the ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\left(\mathrm{mvOC}\right)}$ of combined winter-M and those of the three individual filter samples. From the ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\left(\mathrm{mvOC}\right)}$ and mvOC mass of the individual filters, we can estimate the expected ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\left(\mathrm{mvOC}\right)}$ of the combined winter-M sample using the isotope mass balance equation.$$The ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\left(\mathrm{mvOC}\right)}$ of winter-M calculated from individual filter pieces is 0.524, which is quite similar to the measured ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\left(\mathrm{mvOC}\right)}$ of $\mathrm{0.529}\pm \mathrm{0.007}$ for the combined winter-M with recovery close to 1.

Figure 2

${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\left(\mathrm{mvOC}\right)}$ (shown in circle) and mvOC loading (vertical bar) of the three daily filters (“filter1”, “filter2” and “filter3” shown in dark green) that comprise winter-M and the composite winter-M (shown in blue in the gray area) with mvOC recovery of $\sim \mathrm{1}$. The dashed dark green line indicates the expected mvOC loading and the horizontal solid green line indicates the ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\left(\mathrm{mvOC}\right)}$ for the combined winter-M (weighted mean $\pm$ standard deviation; $\mathrm{0.524}\pm \mathrm{0.028}$). The latter is calculated from the mvOC loading ($\mathrm{\mu }$g cm${}^{-\mathrm{2}}$) and ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\left(\mathrm{mvOC}\right)}$ of the three daily filters using the isotope mass balance equation. ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\left(\mathrm{mvOC}\right)}$ uncertainties are indicated but are too small to be visible.

[Figure omitted. See PDF]

We conclude that ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\left(\mathrm{mvOC}\right)}$ of samples with mvOC recoveries of $\sim \mathrm{1}$ are reliable and used in the following discussion. The mvOC recovery of sample winter-L (red circle in Fig. 1) is also low (0.51), and we could not repeat it due to the limited filter material. For sample winter-L, we take the OC1 concentrations as mvOC and the measured ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\left(\mathrm{mvOC}\right)}$ values but assign a bigger absolute uncertainty of 0.05, due to its low mvOC recoveries. This is based on the difference in ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\left(\mathrm{mvOC}\right)}$ for winter-H and winter-M with low and high mvOC recoveries, which is roughly 0.05 (Fig. 1b).

3.2 mvOC, OC and EC concentrations

Figure 3 shows the concentrations of mvOC, OC and EC and the mvOC contributions to OC (%) for the selected samples in the six studied Chinese megacities. mvOC and OC concentrations averaged $\mathrm{3.3}\pm \mathrm{2.2}$ $\mathrm{\mu }$g m${}^{-\mathrm{3}}$ (range of 0.7–7.4 $\mathrm{\mu }$g m${}^{-\mathrm{3}}$) and $\mathrm{30.0}\pm \mathrm{13.8}$ $\mathrm{\mu }$g m${}^{-\mathrm{3}}$ (8.8–50.4 $\mathrm{\mu }$g m${}^{-\mathrm{3}}$), respectively. EC concentrations ranged from 2.5 to 14.8 $\mathrm{\mu }$g m${}^{-\mathrm{3}}$, with an average of $\mathrm{6.9}\pm \mathrm{3.6}$ $\mathrm{\mu }$g m${}^{-\mathrm{3}}$. High TC concentrations were found in Taiyuan (60 $\mathrm{\mu }$g m${}^{-\mathrm{3}}$ for sample Taiyuan-H), Chongqing (59 $\mathrm{\mu }$g m${}^{-\mathrm{3}}$ for Chongqing-H), Beijing (57 $\mathrm{\mu }$g m${}^{-\mathrm{3}}$ for Beijing-H) and Xi'an (57 $\mathrm{\mu }$g m${}^{-\mathrm{3}}$ for winter-H) in descending order. Of these cities, Chongqing is located in southern China, where there is no official heating season using coal in winter. This study nevertheless indicates severe pollution of carbonaceous aerosols in Chongqing. TC concentrations in the other southern Chinese cities (Shanghai and Guangzhou) were much lower than those in Chongqing (Fig. 3).

Figure 3

Concentrations of mvOC, mrOC and EC ($\mathrm{\mu }$g m${}^{-\mathrm{3}}$) and contribution of mvOC to total OC (%) in three cities in northern China (Xi'an, Beijing and Taiyuan) and three in southern China (Shanghai, Chongqing, Guangzhou). mrOC concentrations are calculated by the difference between OC and mvOC. For details, see Table S4.

[Figure omitted. See PDF]

The fraction of mvOC in total OC ($\mathrm{mvOC}/\mathrm{OC}$ in Fig. 3) gives an indication of OC volatility. The mvOC contributed on average $\mathrm{10.5}\pm \mathrm{3.3}$ % to OC, ranging from 3 % to 15 %. The $\mathrm{mvOC}/\mathrm{OC}$ varies between samples within the same city and between cities, indicating complicated sources and atmospheric processing of OC. The variations might also be partially attributed to the different protocols used for OC quantification. OC in Xi'an is measured with the EUSAAR_2 protocol (up to 650 ${}^{\circ }$C for desorbing OC), whereas the IMPROVE_A protocol (up to 580 ${}^{\circ }$C for OC) was used for the other five cities. However, Han et al. (2016) found that the absolute OC concentrations determined by EUSAAR_2 do not differ much from those determined by IMPROVE_A ($\mathrm{22.6}\pm \mathrm{12.0}$ $\mathrm{\mu }$g m${}^{-\mathrm{3}}$ vs. $\mathrm{19.7}\pm \mathrm{10.7}$ $\mathrm{\mu }$g m${}^{-\mathrm{3}}$) for 1-year PM${}_{\mathrm{2.5}}$ samples in Xi'an during 2012/2013. Because of the small differences of OC between the two protocols, we think the comparison of OC concentrations and $\mathrm{mvOC}/\mathrm{OC}$ amongst the six cities is justified. The rest of OC ($\sim \mathrm{90}$ %) was contributed to the mrOC. mrOC concentrations averaged $\mathrm{26.8}\pm \mathrm{12.0}$ $\mathrm{\mu }$g m${}^{-\mathrm{3}}$, ranging from 7.9 $\mathrm{\mu }$g m${}^{-\mathrm{3}}$ for sample Guangzhou-L to 43.1 $\mathrm{\mu }$g m${}^{-\mathrm{3}}$ for Beijing-H.

Direct comparison of our results with previous works is somewhat difficult because different thermal–optical protocols were used. The averaged mvOC concentration ($\mathrm{3.3}\pm \mathrm{2.2}$ $\mathrm{\mu }$g m${}^{-\mathrm{3}}$) in winter for the six studied sites in China is higher than winter concentrations at an urban background site ($\mathrm{1.6}\pm \mathrm{1.7}$ $\mathrm{\mu }$g m${}^{-\mathrm{3}}$), a rural background site ($\mathrm{0.7}\pm \mathrm{0.6}$) in Prague (Vodička et al., 2015), where 200 ${}^{\circ }$C in He was also applied to desorb this OC fraction using EUSAAR_2 protocol. The $\mathrm{mvOC}/\mathrm{OC}$ ratio in our study is smaller than that of ambient samples from various other locations and much smaller than that of fresh source samples. For example, the urban background site in Prague had a $\mathrm{mvOC}/\mathrm{OC}$ ratio of 28 % and the rural background site of 17 % (Vodička et al., 2015). The contribution of OC1 to total OC was as high as $\sim \mathrm{60}$ % for primary biomass burning measured by EUSAAR_2 thermal–optical protocol (Keller and Burtscher, 2017). For vehicle emissions, the OC fractions desorbed at 140 and 280 ${}^{\circ }$C are the major OC fractions measured by the IMPROVE_A protocol, contributing $\sim \mathrm{30}$ % and $\sim \mathrm{20}$ % to total OC in tunnels in Taiwan, respectively (Zhu et al., 2014, 2010). Using the same protocol for OC analysis, Tian et al. (2017) found that OC fraction desorbed at 140 and 280 ${}^{\circ }$C contributed $\sim \mathrm{13}$ % and $\sim \mathrm{20}$ %, respectively, to primary OC from residential coal combustion in China. The mvOC desorbed at 200 ${}^{\circ }$C falls between 140 and 280 ${}^{\circ }$C, thus the mvOC fraction in total OC should be higher than the faction of OC desorbed at 140 ${}^{\circ }$C in total OC and lower than the fraction of OC desorbed up to 280 ${}^{\circ }$C in total OC.

OC from different emission sources has different volatility and the atmospheric processing can also alter its volatility. Keller and Burtscher (2017) found that aging reduces the volatility of OC from biomass burning; i.e., the contribution of OC1 fraction to total OC decreases from 60 % to 25 % after aging. The photochemical processing of OA can lead to accumulation of carbon in the more refractory organic fraction and also larger organic compounds (Masalaite et al., 2017). The mvOC fraction in OC of our ambient samples is much smaller than that of the primary sources, suggesting that atmospheric aging of OC plays an important role on modifying the volatility of OC.

3.3 Nonfossil and fossil fractions of different carbon fractions

Figure 4a shows the fraction of nonfossil carbon in mvOC, OC and EC ($f_{\mathrm{nf}}$(mvOC), $f_{\mathrm{nf}}$(OC) and $f_{\mathrm{bb}}$(EC), respectively). There are no considerable changes in $f_{\mathrm{nf}}$(mvOC), $f_{\mathrm{nf}}$(OC) and $f_{\mathrm{bb}}$(EC) between polluted days (“H” samples) and clean days (“L” samples) within each site (Fig. 4a), despite the very different concentrations of carbonaceous aerosols (Fig. 3). However, $f_{\mathrm{nf}}$(mvOC), $f_{\mathrm{nf}}$(OC) and $f_{\mathrm{bb}}$(EC) varied significantly among different sites: the lowest values are always found in Taiyuan and the highest are found in Chongqing. This implies that different pollution patterns exist in individual Chinese cities. The smallest $f_{\mathrm{nf}}$(mvOC), $f_{\mathrm{nf}}$(OC) and $f_{\mathrm{bb}}$(EC) being in Taiyuan suggests that fossil sources are the main contributor to mvOC, OC and EC, whereas the largest values being in Chongqing shows that the nonfossil contributions to mvOC, OC and EC are evidently higher in Chongqing than in other sites.

Figure 4

(a) Fraction of nonfossil ($f_{\mathrm{nf}}$) in mvOC, OC and EC. $f_{\mathrm{nf}}$(mrOC) is calculated by the differences between OC and mvOC using isotope mass balance equation. (b) $\mathrm{mvOC}/\mathrm{OC}$ ratios for nonfossil, fossil and total OC.

[Figure omitted. See PDF]

Ranging from 22 % to 55 %, $f_{\mathrm{nf}}$(mvOC) ($\mathrm{41}\pm \mathrm{11}$ %) is consistently smaller than $f_{\mathrm{nf}}$(OC) ($\mathrm{52}\pm \mathrm{16}$ %, with a range of 29 %–77 %) and larger than $f_{\mathrm{bb}}$(EC) ($\mathrm{27}\pm \mathrm{9}$ %, with a range of 10 %–42 %), despite their variations among the cities. The absolute difference in the nonfossil fractions between mvOC and OC is 13 % (7 %–25 %), similar to that between mvOC and EC (13 %, with a range of 4 %–25 %). Consistently smaller $f_{\mathrm{nf}}$(mvOC) than $f_{\mathrm{nf}}$(OC) suggests that mvOC is more fossil (less nonfossil) than the total OC. Liu et al. (2017) also found that ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$ of OC desorbed at lower temperature (up to 200 ${}^{\circ }$C) in He was 0.389, smaller than ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$ of total OC desorbed in the He phase of EUSAAR_2 protocol (${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$ $=\mathrm{0.495}$, up to 650 ${}^{\circ }$C) for a single test sample collected in winter in Xinjiang, China. Agrios et al. (2016) had the qualitatively similar finding that ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$ of OC extracted in ${\mathrm{O}}_{\mathrm{2}}$ at 200 ${}^{\circ }$C was smaller than ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$ of OC extracted in ${\mathrm{O}}_{\mathrm{2}}$ at higher temperatures for samples collected from an urban and rural site in Switzerland and from the Los Angeles Basin, USA. We can unambiguously conclude from ${}^{\mathrm{14}}\mathrm{C}$ measurements of mvOC and OC that mvOC is more fossil than OC in six Chinese cites. The fraction of nonfossil carbon in mrOC ($f_{\mathrm{nf}}$(mrOC)) is calculated by the differences between OC and mvOC using the isotope mass balance. Since mvOC is only a small fraction of OC, $f_{\mathrm{nf}}$(mrOC) is very similar to $f_{\mathrm{nf}}$(OC).

Primary OC from vehicular emissions is generally more volatile (i.e., less refractory) than OC from biomass burning, which is in line with the ambient results presented. For example, Grieshop et al. (2009c) constrained the volatility distribution of primary OA from a diesel engine and wood burning using measurements of TD coupled to aerosol mass spectrometer and found that OA from wood burning was less volatile than from diesel exhaust. Chen et al. (2007) measured $\mathrm{OC}/\mathrm{EC}$ from fresh biomass burning emissions and found that high-temperature OC fractions (desorbed from 450 to 550 ${}^{\circ }$C) are the major fraction of OC, in contrast to gasoline and diesel exhausts, where OC fractions desorbed at temperature lower than 250 ${}^{\circ }$C are more abundant (Watson et al., 1994; Chow et al., 2004). Zhu et al. (2010, 2014) also found that OC desorbed at 140 ${}^{\circ }$C is the major fraction of OC in fresh vehicle emissions in tunnel experiments. A more recent study of Masalaite et al. (2017) found that OC desorbed at 200 ${}^{\circ }$C in He has higher contributions from vehicular emissions than OC desorbed at higher temperatures (250–350 ${}^{\circ }$C) at three sites (urban, coastal and forest) in Lithuania, based on ${\mathit{\delta }}^{\mathrm{13}}\mathrm{C}$ of OC desorbed at a different temperature step. However, the conclusions of this study remained qualitative because ${\mathit{\delta }}^{\mathrm{13}}\mathrm{C}$ of organic carbon can be changed significantly by isotopic fractionation during atmospheric processing, whereas the ${}^{\mathrm{14}}\mathrm{C}$ source apportionment used in our study is independent of atmospheric processing.

By comparing ($\mathrm{mvOC}/\mathrm{OC}$)${}_{\mathrm{fossil}}$ (i.e., the fraction of mvOC${}_{\mathrm{fossil}}$ in OC${}_{\mathrm{fossil}}$) and ($\mathrm{mvOC}/\mathrm{OC}$)${}_{\mathrm{nf}}$, we can get a measure of the volatility of OC from fossil and nonfossil sources, respectively: 13$$\begin{array}{rl}{\displaystyle \frac{(\mathrm{mvOC}/\mathrm{OC}{)}_{\mathrm{fossil}}}{(\mathrm{mvOC}/\mathrm{OC}{)}_{\mathrm{nf}}}}& ={\displaystyle \frac{{\mathrm{mvOC}}_{\mathrm{fossil}}/{\mathrm{mvOC}}_{\mathrm{nf}}}{{\mathrm{OC}}_{\mathrm{fossil}}/{\mathrm{OC}}_{\mathrm{nf}}}}\\ & ={\displaystyle \frac{\frac{\mathrm{1}}{f_{\mathrm{nf}}\left(\mathrm{mvOC}\right)}-\mathrm{1}}{\frac{\mathrm{1}}{f_{\mathrm{nf}}\left(\mathrm{OC}\right)}-\mathrm{1}}}.\end{array}$$ As$f_{\mathrm{nf}}$(mvOC) is smaller than $f_{\mathrm{nf}}$(OC), this is equivalent to ($\mathrm{mvOC}/\mathrm{OC}$)${}_{\mathrm{fossil}}$ >($\mathrm{mvOC}/\mathrm{OC}$)${}_{\mathrm{nf}}$ (Eq. 13). That is, OC from fossil sources contains a larger, more volatile fraction than OC from nonfossil sources. This is true for all studied cities, despite the variation in $f_{\mathrm{nf}}$(mvOC) and $f_{\mathrm{nf}}$(OC) in different cities. ($\mathrm{mvOC}/\mathrm{OC}$)${}_{\mathrm{fossil}}$($\mathrm{14}\pm \mathrm{6.6}$ %, with a range of 3.7 %–28 %) is consistently higher and more variable than ($\mathrm{mvOC}/\mathrm{OC}$)${}_{\mathrm{nf}}$ ($\mathrm{7.5}\pm \mathrm{2.9}$ %, with a range of 2.6 %–11.6 %; Fig. 4b). ($\mathrm{mvOC}/\mathrm{OC}$)${}_{\mathrm{fossil}}$ and ($\mathrm{mvOC}/\mathrm{OC}$)${}_{\mathrm{nf}}$ in general tracks the variation in $\mathrm{mvOC}/\mathrm{OC}$, except that ($\mathrm{mvOC}/\mathrm{OC}$)${}_{\mathrm{fossil}}$ for Chongqing (28 % for Chongqing-H and 21 % for Chongqing-L) is much higher than the other sites (averaged 12 %). When comparing the $\mathrm{mvOC}/\mathrm{OC}$ within each city, it is found that ($\mathrm{mvOC}/\mathrm{OC}$)${}_{\mathrm{fossil}}$ changes more strongly between H and L samples than ($\mathrm{mvOC}/\mathrm{OC}$)${}_{\mathrm{nf}}$, indicated by the bigger absolute differences of ($\mathrm{mvOC}/\mathrm{OC}$)${}_{\mathrm{fossil}}$ between H and L samples than that of ($\mathrm{mvOC}/\mathrm{OC}$)${}_{\mathrm{nf}}$. Correlations between $f_{\mathrm{nf}}$(mvOC) and $f_{\mathrm{nf}}$(OC) and between $f_{\mathrm{nf}}$(mvOC) and $f_{\mathrm{bb}}$(EC) are examined. $f_{\mathrm{nf}}$(mvOC) correlates more closely with $f_{\mathrm{nf}}$(OC) ($r^{\mathrm{2}}=\mathrm{0.89}$) than with $f_{\mathrm{bb}}$(EC) ($r^{\mathrm{2}}=\mathrm{0.62}$) (Fig. 5).

Figure 5

Correlation between $f_{\mathrm{nf}}$(mvOC), $f_{\mathrm{nf}}$(OC), $f_{\mathrm{nf}}$(mvOC) and $f_{\mathrm{bb}}$(EC) for Xi'an, Beijing, Taiyuan, Shanghai, Chongqing and Guangzhou during winter.

[Figure omitted. See PDF]

It is likely that some of mvOC are of secondary origin, since SOC is formed in the atmosphere and subsequently condenses onto the aerosol particles (Hallquist et al., 2009; Jimenez et al., 2009). It has been shown that SOC formed in chambers initially desorbed at relatively low temperatures (e.g., 200 ${}^{\circ }$C) (Holzinger et al., 2010; Salo et al., 2011; Meusinger et al., 2017; Gkatzelis et al., 2018). However, there is also evidence that the formation of highly oxidized low-volatility compounds can occur in SOA formation (Ehn et al., 2014). We thus compared concentrations and sources of SOC to those of mvOC. SOC concentrations and $f_{\mathrm{fossil}}$(SOC) are estimated based on the ${}^{\mathrm{14}}\mathrm{C}$-apportioned OC and EC in combination with $p$ values of 0–1 (Sect. 2.4). Where ${\mathit{\delta }}^{\mathrm{13}}\mathrm{C}$ measurements were available for high-volume filter samples (Beijing, Shanghai and Guangzhou), we compare the SOC concentrations and $f_{\mathrm{fossil}}$(SOC) derived from $p$ values of 0–1 with $p$ values constrained by ${\mathit{\delta }}^{\mathrm{13}}\mathrm{C}$. In these three cities, SOC concentrations and $f_{\mathrm{fossil}}$(SOC) without constraints on $p$ value (0–1) is very similar to that with constraints on $p$ value using ${\mathit{\delta }}^{\mathrm{13}}\mathrm{C}$, only with relatively larger uncertainties (Fig. 6). This indicates that the unconstrained $p$ values do not lead to significant bias in SOC concentrations and $f_{\mathrm{fossil}}$(SOC) and the main gain in using ${\mathit{\delta }}^{\mathrm{13}}\mathrm{C}$ is currently a decrease in uncertainty.

Figure 6

(a) Observed mvOC concentrations (shown as a green bar) and estimated SOC concentrations by ${}^{\mathrm{14}}\mathrm{C}$-apportioned OC and EC in combination with $p$ values without ($p=\mathrm{0}$–1) and with constraints by ${\mathit{\delta }}^{\mathrm{13}}\mathrm{C}$ (${}^{\mathrm{13}}\mathrm{C}$-constrained $p$), shown as solid orange bars and patterned orange bars, respectively. Error bars indicate the uncertainties of mvOC and SOC concentrations calculated by Monte Carlo error propagation. (b) For Chongqing, lower $p$ values of 0–0.5 are used to estimate SOC concentrations and are compared with $p$ values of 0–1. (c) Fraction of fossil ($f_{\mathrm{fossil}}$) in mvOC (blue square) and SOC. $f_{\mathrm{fossil}}$ of SOC is estimated using $p$ values without constraints ($p=\mathrm{0}$–1) are shown as red circles and with constraints by ${\mathit{\delta }}^{\mathrm{13}}\mathrm{C}$ (${}^{\mathrm{13}}\mathrm{C}$-constrained $p$) as purple triangles, respectively. For Chongqing, lower $p$ values of 0–0.5 are used (purple circles). Interquartile ranges (25th–75th percentile) of the median $f_{\mathrm{fossil}}$(SOC) are shown as vertical bars, uncertainties of $f_{\mathrm{fossil}}$(mvOC) are indicated but are too small to be visible.

[Figure omitted. See PDF]

When choosing $p$ randomly from 0 to 1, the calculated median $f_{\mathrm{fossil}}$(SOC) is smaller than 0 for samples in Chongqing, reflecting negative SOC${}_{\mathrm{fossil}}$ values. Since SOC${}_{\mathrm{fossil}}$ is calculated by subtracting POC${}_{\mathrm{fossil}}$ from OC${}_{\mathrm{fossil}}$, this indicates POC${}_{\mathrm{fossil}}$ is overestimated in Chongqing. Because the primary $\mathrm{OC}/\mathrm{EC}$ ratio is higher for coal ($r_{\mathrm{coal}}$) than vehicle ($r_{\mathrm{vehicle}}$) emissions, an overestimate of POC${}_{\mathrm{fossil}}$ is the direct result of overestimated $p$ values, i.e. a too high fraction of coal burning in EC${}_{\mathrm{fossil}}$ (Eq. 11). Chongqing is in southern China, where there is no official heating using coal in winter and smaller contribution of coal combustion is expected than in northern China. Further, earlier studies found that vehicle emissions contributed more to EC than coal combustion in Chongqing (Chen et al., 2017). We thus tried to lower $p$ values to 0–0.5. Lower $p$ values lead to higher $f_{\mathrm{fossil}}$(SOC) (Fig. 6c) and, to a lesser extent, higher SOC concentrations (Fig. 6b). For Chongqing, we took $p$ from 0 to 0.5 as a conservative estimate. For other cities except Chongqing, we kept $p$ values of 0–1 due to the very similar $f_{\mathrm{fossil}}$(SOC) and SOC concentrations derived from $p$ values with and without constraints by ${\mathit{\delta }}^{\mathrm{13}}\mathrm{C}$, especially for the two southern cities in China, e.g., Shanghai and Guangzhou.

$f_{\mathrm{fossil}}$(mvOC) ranged from 45 % to 78 %, with an average of 59 %, whereas $f_{\mathrm{fossil}}$(SOC) ranged from 10 % to 70 % and averaged 43 %. $f_{\mathrm{fossil}}$(mvOC) is always larger than $f_{\mathrm{fossil}}$(SOC) (Fig. 6c), showing that mvOC is more fossil than SOC. However, both SOC and mvOC have a larger contribution of fossil sources than the total OC. The smallest differences between $f_{\mathrm{fossil}}$(SOC) and $f_{\mathrm{fossil}}$(mvOC) are found in Beijing for both samples Beijing-H and Beijing-L, with an absolute difference of 3 % and 2 %, respectively, within the uncertainty of the measurements. Much higher differences are found in other cities, with the absolute differences ranging from 13 % to 35 %, well outside the measurement uncertainties.

Median values for SOC, SOC${}_{\mathrm{fossil}}$ and SOC${}_{\mathrm{nf}}$ are very close to their mean values (Fig. S6, Tables S7, S8). To compare with mvOC, here we use mean ($\pm$SD) as well for SOC, SOC${}_{\mathrm{fossil}}$ and SOC${}_{\mathrm{nf}}$. The lowest SOC concentrations are found to be 2.95 $\mathrm{\mu }$g m${}^{-\mathrm{3}}$ for Guangzhou-L and the highest are found to be 35.9 $\mathrm{\mu }$g m${}^{-\mathrm{3}}$ for Beijing-H. SOC concentrations ($\mathrm{15.4}\pm \mathrm{9.0}$ $\mathrm{\mu }$g m${}^{-\mathrm{3}}$) are consistently higher than those of mvOC ($\mathrm{3.3}\pm \mathrm{2.2}$ $\mathrm{\mu }$g m${}^{-\mathrm{3}}$, Fig. 6a), showing that the contribution of SOC to total OC is considerably larger than that of mvOC. The major fraction of SOC has desorption temperatures above 200 ${}^{\circ }$C and falls into our operationally defined mrOC. SOC${}_{\mathrm{fossil}}$ ranged from 0.90 to 25.3 $\mathrm{\mu }$g m${}^{-\mathrm{3}}$, with an average of $\mathrm{6.8}\pm \mathrm{6.4}$ $\mathrm{\mu }$g m${}^{-\mathrm{3}}$, which is consistently higher than mvOC${}_{\mathrm{fossil}}$ ranging from 0.5 to 5.4 $\mathrm{\mu }$g m${}^{-\mathrm{3}}$ with an average of $\mathrm{2.1}\pm \mathrm{1.4}$ $\mathrm{\mu }$g m${}^{-\mathrm{3}}$ (Table S6). mvOC${}_{\mathrm{nf}}$ is smaller than SOC${}_{\mathrm{nf}}$ as well (Fig. S7), averaged mvOC${}_{\mathrm{nf}}$ concentrations ($\mathrm{1.5}\pm \mathrm{1.1}$ $\mathrm{\mu }$g m${}^{-\mathrm{3}}$) are only one-fifth of that SOC${}_{\mathrm{nf}}$($\mathrm{8.5}\pm \mathrm{5}$ $\mathrm{\mu }$g m${}^{-\mathrm{3}}$).

Figure 7

Correlation between (a) OC${}_{\mathrm{nf}}$ and mvOC${}_{\mathrm{nf}}$ (blue cross), (b) POC${}_{\mathrm{bb}}$ and mvOC${}_{\mathrm{nf}}$ (orange diamond), SOC${}_{\mathrm{nf}}$ and mvOC${}_{\mathrm{nf}}$ (green circle), (c) POC${}_{\mathrm{fossil}}$ and mvOC${}_{\mathrm{fossil}}$(diamond in dark orange), and SOC${}_{\mathrm{fossil}}$ and mvOC${}_{\mathrm{fossil}}$(circle in dark green) for Xi'an, Beijing, Taiyuan, Shanghai, Chongqing and Guangzhou during winter. In panel (c), the data points marked in red are from sample Beijing-H with the considerably higher SOC${}_{\mathrm{fossil}}$ $/$ POC${}_{\mathrm{fossil}}$ ratio (3.1) than the other samples (an average of 0.8) and are not included in the regression.

[Figure omitted. See PDF]

We examined the correlation between mvOC and SOC from nonfossil sources (mvOC${}_{\mathrm{nf}}$ vs. SOC${}_{\mathrm{nf}}$) and from fossil sources (mvOC${}_{\mathrm{fossil}}$ vs. SOC${}_{\mathrm{fossil}}$). Comparing absolute concentrations between different cities and pollution levels is expected to result in reasonably positive correlations for all species because it can be expected that all carbon fractions are higher in polluted cities than in clean cities. However, it is still interesting to compare the correlation coefficients relative to one another. mvOC${}_{\mathrm{nf}}$ includes primary and secondary mvOC from nonfossil sources. Good correlation was found between mvOC${}_{\mathrm{nf}}$ and SOC${}_{\mathrm{nf}}$($r^{\mathrm{2}}=\mathrm{0.87}$, bootstrapped 95 % confidence interval of [0.21, 0.94]; Fig. 7), which is not significantly different from the correlation between mvOC${}_{\mathrm{nf}}$ and POC${}_{\mathrm{bb}}$ ($r^{\mathrm{2}}=\mathrm{0.71}$, bootstrapped 95 % confidence interval of [0.60, 0.99]). SOC${}_{\mathrm{nf}}$ is estimated by subtracting POC${}_{\mathrm{bb}}$ from OC${}_{\mathrm{nf}}$, where POC${}_{\mathrm{bb}}$ is estimated by multiplying EC${}_{\mathrm{bb}}$ with the $\mathrm{OC}/\mathrm{EC}$ ratio of fresh biomass burning plumes ($r_{\mathrm{bb}}$) (Eq. 7). Due to the large uncertainty in $r_{\mathrm{bb}}$, separation between POC${}_{\mathrm{bb}}$ and SOC${}_{\mathrm{nf}}$ is rather uncertain and the $r^{\mathrm{2}}$ for a correlation of mvOC${}_{\mathrm{nf}}$ and SOC${}_{\mathrm{nf}}$ is therefore affected by $r_{\mathrm{bb}}$. For example, with lower $r_{\mathrm{bb}}$ of 3 (mean), the correlation between SOC${}_{\mathrm{nf}}$ and mvOC${}_{\mathrm{nf}}$ is stronger ($r^{\mathrm{2}}=\mathrm{0.92}$) than for our best estimate of $r_{\mathrm{bb}}$(4) ($r^{\mathrm{2}}=\mathrm{0.87}$). For $r_{\mathrm{bb}}$ of 5, it is slightly weaker ($r^{\mathrm{2}}=\mathrm{0.76}$, Fig. S8). OC${}_{\mathrm{nf}}$, the sum of POC${}_{\mathrm{bb}}$ and SOC${}_{\mathrm{nf}}$, is well constrained by ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\mathrm{OC}}$ and not affected by $r_{\mathrm{bb}}$. mvOC${}_{\mathrm{nf}}$ correlates strongly with OC${}_{\mathrm{nf}}$ ($r^{\mathrm{2}}=\mathrm{0.91}$), suggesting strong impacts on the variability of mvOC${}_{\mathrm{nf}}$ from nonfossil emissions, including the secondary formation from nonfossil sources and primary biomass burning emissions.

SOC${}_{\mathrm{fossil}}$ is produced by the oxidation of VOCs from fossil origins, including coal combustion and vehicle emissions, whereas mvOC${}_{\mathrm{fossil}}$ can have both primary and secondary sources. No correlation is found between mvOC${}_{\mathrm{fossil}}$ and SOC${}_{\mathrm{fossil}}$ ($r^{\mathrm{2}}=\mathrm{0.01}$, Fig. 7c), implying that in cities or pollution conditions characterized by high concentrations of fossil SOC, mvOC${}_{\mathrm{fossil}}$ is not concurrently elevated. However, as demonstrated by the good correlation between POC${}_{\mathrm{fossil}}$ and mvOC${}_{\mathrm{fossil}}$ ($r^{\mathrm{2}}=\mathrm{0.69}$; Fig. 7c) we conclude that primary fossil sources may play an important role on the variation in mvOC${}_{\mathrm{fossil}}$. Larger $f_{\mathrm{fossil}}$(mvOC) than $f_{\mathrm{fossil}}$(SOC) also indicates that fossil sources are associated more with mvOC than SOC. In Fig. 7c there are two outlier data points from sample Beijing-H with strongly elevated SOC${}_{\mathrm{fossil}}$ and lower POC${}_{\mathrm{fossil}}$ (shown in red circle and diamond, respectively). The highest SOC${}_{\mathrm{fossil}}$ $/$ POC${}_{\mathrm{fossil}}$ ratio (3.1) is found in Beijing-H, suggesting the largest fossil SOC formation during polluted days in Beijing compared to clean days in Beijing (Beijing-L with SOC${}_{\mathrm{fossil}}$ $/$ POC${}_{\mathrm{fossil}}$ ratio of 1.4) and samples in other cities (average SOC${}_{\mathrm{fossil}}$ $/$ POC${}_{\mathrm{fossil}}$ ratio of 0.7). This is also found in an early study that during winter 2013, SOC${}_{\mathrm{fossil}}$ $/$ POC${}_{\mathrm{fossil}}$ ratios in Beijing (with an average of 4.2) were much higher than those in Xi'an, Shanghai and Guangzhou (average SOC${}_{\mathrm{fossil}}$ $/$ POC${}_{\mathrm{fossil}}$ ratio of 1.3) (Zhang et al., 2015). This pollution event of Beijing-H is thus not included in the regression.

4 Discussion

In this study, samples are collected on single bare quartz (BQ) filters. mvOC may be influenced by adsorption of organic vapors to the quartz fibers (positive sampling artifact). In other words, the measured mvOC in the earlier discussion includes particulate mvOC and gaseous carbon adsorbed to the quartz filters that can also be desorbed from the filters up to 200 ${}^{\circ }$C. This artifact can be estimated by using a backup quartz filter behind the BQ filter (quartz behind quartz, QBQ). The QBQ filter is assumed to be exposed to the same organic vapor concentrations as the front BQ filter and the carbon measured on the QBQ filter is mainly due to adsorbed organic vapor and can serve as an estimate of the positive artifact (Chow et al., 2010; Zhu et al., 2012).

Figure 8

(a) The predicted particulate fraction of OC fractions (i.e., OC${}_{\mathrm{140}{}^{\circ }\mathrm{C}}$, OC${}_{\mathrm{280}{}^{\circ }\mathrm{C}}$, OC${}_{\mathrm{480}{}^{\circ }\mathrm{C}}$, OC${}_{\mathrm{580}{}^{\circ }\mathrm{C}}$, where the subscript refers to the desorption temperature of the IMPROVE_A protocol) as function of OC concentrations (Sect. S4). Hollow shapes are used to indicate samples collected by high-volume samplers (HiVol) at a flow rate of 1.0 m${}^{\mathrm{3}}$ min${}^{-\mathrm{1}}$ and filled shapes by mini-volume samplers (MiniVol) at a flow rate of 5 L min${}^{-\mathrm{1}}$. Particulate fraction of mvOC ($X_{\mathrm{p}\mathrm{\_}\mathrm{mvOC}}$, blue cross) is estimated by interpolating the particulate fraction of OC${}_{\mathrm{140}{}^{\circ }\mathrm{C}}$ and OC${}_{\mathrm{280}{}^{\circ }\mathrm{C}}$. The blue solid line is the fitted line for the $X_{\mathrm{p}\mathrm{\_}\mathrm{mvOC}}$. (b) The $X_{\mathrm{p}\mathrm{\_}\mathrm{mvOC}}$ for each sample corresponding to the fitted line in (a) for each sample; (c, d) fraction of nonfossil in mvOC collected on the single bare quartz (BQ) filter ($f_{\mathrm{nf}}$(mvOC), blue line and solid circle). The blue area indicates the possible range of the nonfossil fraction of particulate carbon after correction for the positive artifact ($f_{\mathrm{nf}}$(p_mvOC)), using an upper limit estimate for the positive artifact. The upper and lower limits of the range correspond to the extreme assumption that the ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$ of the adsorbed organic vapor is 50 % lower and 50 % higher than measured ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\left(\mathrm{mvOC}\right)}$, respectively. $f_{\mathrm{nf}}$(p_mvOC) is compared with $f_{\mathrm{nf}}$(OC) (solid green triangle) in (c) and $f_{\mathrm{nf}}$(SOC) (red circle) in (d). Interquartile range (25th–75th percentile) of the median $f_{\mathrm{fossil}}$ (SOC) is shown as vertical bars.

[Figure omitted. See PDF]

Without the QBQ backup quartz filters in this study, we estimate the particulate fraction of mvOC ($X_{\mathrm{p}\mathrm{\_}\mathrm{mvOC}}$) on the BQ filters using the OC fractions measured by the thermal–optical carbon analyzer in combination with an empirical volatility basis set (VBS) model, following the method of Ma et al. (2016) (Sect. S4). The VBS bins the organic aerosol (OA) compounds according to the effective saturation concentrations ($C^{*}$). For example, if $C^{*}$ is equal to the concentrations of OA ($C_{\mathrm{OA}}$), then 50 % of OA mass will be in the particle phase (Donahue et al., 2006). Ma et al. (2016) provide empirical $C^{*}$ values for each desorption temperature corresponding to the measured particulate fraction on the filter, determined by the comparison of BQ and QBQ filters. They show that these empirical $C^{*}$ values (10${}^{\mathrm{1.6}}$ for OC${}_{\mathrm{140}{}^{\circ }\mathrm{C}}$, 10${}^{\mathrm{1.1}}$ for OC${}_{\mathrm{280}{}^{\circ }\mathrm{C}}$, 10${}^{\mathrm{0.6}}$ for OC${}_{\mathrm{480}{}^{\circ }\mathrm{C}}$, 10${}^{\mathrm{0.1}}$ for OC${}_{\mathrm{580}{}^{\circ }\mathrm{C}}$, where the subscript refers to the desorption temperature of the IMPROVE_A protocol) are useful to predict particulate fractions for various OC loadings and sampling environments. Assuming $\mathrm{OA}/\mathrm{OC}$ mass ratio is 1.8 for urban environment (Xing et al., 2013), the particulate fraction of the collected OC for each of the four IMPROVE_A temperature steps $i(X_{\mathrm{p}i}$) can be estimated as a function of OC concentration, as shown in Fig. 8a. The particulate fraction of OC increases considerably from OC${}_{\mathrm{140}{}^{\circ }\mathrm{C}}$ to OC${}_{\mathrm{580}{}^{\circ }\mathrm{C}}$ and with increasing OC concentration. The desorption temperature for mvOC is 200 ${}^{\circ }$C, falling within the 140 ${}^{\circ }$C for OC${}_{\mathrm{140}{}^{\circ }\mathrm{C}}$ and 280 ${}^{\circ }$C for OC${}_{\mathrm{280}{}^{\circ }\mathrm{C}}$. We thus estimated the $X_{\mathrm{p}\mathrm{\_}\mathrm{mvOC}}$ will fall within the $X_{\mathrm{p}\mathrm{1}}$ for OC${}_{\mathrm{140}{}^{\circ }\mathrm{C}}$ and $X_{\mathrm{p}\mathrm{2}}$ for OC${}_{\mathrm{280}{}^{\circ }\mathrm{C}}$. Particulate fraction of mvOC is estimated by interpolating the particulate fraction of OC${}_{\mathrm{140}{}^{\circ }\mathrm{C}}$ and OC${}_{\mathrm{280}{}^{\circ }\mathrm{C}}$ (Fig. 8a). $X_{\mathrm{p}\mathrm{\_}\mathrm{mvOC}}$ ranged from 0.54 for sample winter-L with the lowest OC concentration of 15 $\mathrm{\mu }$g m${}^{-\mathrm{3}}$ to 0.79 for sample Beijing-H with the highest OC concentration of 50 $\mathrm{\mu }$g m${}^{-\mathrm{3}}$ (Fig. 8b). However, for our study we consider this as a lower estimate of particulate fraction of mvOC (i.e., higher limits of the positive artifact).

Our sampling duration (24 h) is much longer than that of Ma et al. (2016), which is 45 min or 1 h. With longer sampling duration, the loading on the quartz filter of both particulate OC and adsorbed organic vapors increase initially; however, the latter will reach a saturation point and level out. Thus, the percent bias by the adsorbed organics is smaller for longer sampling duration (Turpin et al., 2000). For example, Turpin et al. (1994) found that the loading of adsorbed organics reached around 1.7 $\mathrm{\mu }$g cm${}^{-\mathrm{2}}$ and remained constant after 200 min of sampling time, but the aerosol loading increased linearly with the sampling time.

The empirically determined $C^{*}$ is affected by the sampling conditions, e.g., temperature. $C_i^{*}$ is defined at 300 K, higher than the average temperature in winter in China. The lower temperature drives the semi-volatile OC into the particle phase and leads to higher particulate fraction of each OC fraction (Donahue et al., 2006).

The adsorption of organic vapors is also influenced by the face velocity (the volumetric flow rate divided by the exposed filter area). As the face velocity decreases, the adsorption of organic gaseous species increases (McDow and Huntzicker, 1990; Turpin et al., 2000; Viana et al., 2006). In this study, most of samples were collected using high-volume samplers (HiVol) with a flow rate of 1 m${}^{\mathrm{3}}$ min${}^{-\mathrm{1}}$, including samples collected in Xi'an, Beijing, Guangzhou and Shanghai. Only samples in Taiyuan and Chongqing were collected using mini-volume samplers (MiniVol) at a flow rate of 5 L min${}^{-\mathrm{1}}$. The face velocities for the HiVol and MiniVol sampler were 42 and 7 cm s${}^{-\mathrm{1}}$, respectively. Turpin et al. (1994) found that the apparent OC concentrations on the front quartz filter were 22 % greater with a face velocity of 20 cm s${}^{-\mathrm{1}}$ than those with a face velocity of 40 cm s${}^{-\mathrm{1}}$. A more recent study of Kim et al. (2016) found that the positive artifacts for OC using sampler with the face velocity of 42 cm s${}^{-\mathrm{1}}$ is 7.2 %, which is smaller than 12.2 % using sampler of a face velocity of 24 cm s${}^{-\mathrm{1}}$ in spring and fall of 2014 in Korea. We do not take this into consideration and apply the same $C_i^{*}$ for HiVol and MiniVol samples, due to the lake of $C_i^{*}$ for different face velocities.

Due to the lack of information on how sampling duration, temperature and face velocity affect the $C^{*}$ for each carbon fraction, we treat the particulate fraction in Fig. 8b as a rough estimate. But with long sampling duration, low temperature and high face velocity for most of samples as well as high OC concentrations in winter in this study, the particulate fraction of mvOC on the collected filters should be considerably higher than the estimate derived from VBS. For example, Cheng and He (2015) found that the positive artifacts in OC were lowest ($\sim \mathrm{10}$ %) in winter, much lower than in summer ($\sim \mathrm{23}$ %) in Beijing, China. With the same face velocity of $\sim \mathrm{10}$ cm s${}^{-\mathrm{1}}$, the smaller positive artifacts in winter than that in summer can be attributed to higher OC concentrations (35 $\mathrm{\mu }$g m${}^{-\mathrm{3}}$ in winter vs. 12 $\mathrm{\mu }$g m${}^{-\mathrm{3}}$ in summer) and lower temperatures in winter. We consider the positive sampling artifact in this study is small and the collected mvOC on the bare quartz filter is dominated by the particulate mvOC.

Even though the positive artifact introduced by the adsorbed organic vapors is likely small, it is worth investigating its potential effect on the ${}^{\mathrm{14}}\mathrm{C}$ results. The measured ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\left(\mathrm{mvOC}\right)}$ on the BQ filter can be explained by the ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$ of particulate mvOC (${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{(\mathrm{p}}\mathit{\_}\mathrm{mvOC})$) and adsorbed organic vapors (${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{(\mathrm{g}}\mathit{\_}\mathrm{mvOC})$). The ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{(\mathrm{p}}\mathit{\_}\mathrm{mvOC}$) can be calculated as follows: 14$$\begin{array}{rl}& {\mathrm{F}}^{\mathrm{14}}{\mathrm{C}}_{\left(\mathrm{p}\mathrm{\_}\mathrm{mvOC}\right)}\\ & ={\displaystyle \frac{{\mathrm{F}}^{\mathrm{14}}{\mathrm{C}}_{\left(\mathrm{mvOC}\right)}-(\mathrm{1}-X_{\mathrm{p}\mathrm{\_}\mathrm{mvOC}})\times {\mathrm{F}}^{\mathrm{14}}{\mathrm{C}}_{\left(\mathrm{g}{\mathrm{\_}}_{\mathrm{mvOC}}\right)}}{X_{\mathrm{p}\mathrm{\_}\mathrm{mvOC}}}},\end{array}$$ where $X_{\mathrm{p}\mathrm{\_}\mathrm{mvOC}}$ is the estimated particulate fraction of mvOC, as shown in Fig. 8b, which is a much lower limit, as discussed earlier. (1–$X_{\mathrm{p}\mathrm{\_}\mathrm{mvOC}}$) is thus an estimate of positive artifact. How close the measured ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\left(\mathrm{mvOC}\right)}$ is to the ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\left(\mathrm{p}\mathrm{\_}\mathrm{mvOC}\right)}$ depends on the ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\left(\mathrm{g}\mathrm{\_}\mathrm{mvOC}\right)}$.

Therefore, we conducted a sensitivity study based on the following assumption: the positive artifact is dominated by semi-volatile gases that actively partition between gas and particle phase. As the base case, we assume that ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\left(\mathrm{g}\mathrm{\_}\mathrm{mvOC}\right)}=$ ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\left(\mathrm{mvOC}\right)}$ and consider extreme cases of ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\left(\mathrm{g}\mathrm{\_}\mathrm{mvOC}\right)}=\mathrm{1.5}\times$ ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\left(\mathrm{mvOC}\right)}$ and ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\left(\mathrm{g}\mathrm{\_}\mathrm{mvOC}\right)}=\mathrm{0.5}\times$ ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\left(\mathrm{mvOC}\right)}$, i.e. a variation of $\pm \mathrm{50}$ %. This allows ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$ values for the adsorbed gas phase to be as low as for EC and higher than for OC in most cases. When the lower limit of ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\left(\mathrm{g}\mathrm{\_}\mathrm{mvOC}\right)}$ is used (50 % lower than ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\left(\mathrm{mvOC}\right)}$), the $f_{\mathrm{nf}}$(p_mvOC), calculated by dividing ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\left(\mathrm{p}\mathrm{\_}\mathrm{mvOC}\right)}$ with ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\mathrm{nf}}$, is 0.04 to 0.19 (absolute value) higher than $f_{\mathrm{nf}}$(mvOC) (Fig. 8) but still lower than $f_{\mathrm{nf}}$(OC) and $f_{\mathrm{nf}}$(SOC) for most of samples, at least for highly polluted (H) samples with a higher particulate fraction of mvOC. Consequently, our conclusion that mvOC is more fossil than OC and SOC is still valid for most of samples, especially since this sensitivity study assumes (i) an upper limit of the positive artifact as well as (ii) a very large maximum difference between the ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$ of the adsorbed gas phase and the semi-volatile particle phase. The upper limit of $f_{\mathrm{nf}}$(p_mvOC) corresponds to a semi-volatile gas phase with ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$ values in the range of EC, which is probably not realistic.

5 Conclusions

Radiocarbon measurements are conducted for EC and OC of PM${}_{\mathrm{2.5}}$ samples collected in six Chinese megacities during winter. In addition, the sources of mvOC are measured unambiguously using ${}^{\mathrm{14}}\mathrm{C}$. mvOC is isolated by desorbing organics in He at 200 ${}^{\circ }$C using our custom-made aerosol combustion system with a reproducibility of $\sim \mathrm{5}$ %. This new isolation method is evaluated by comparing with the OC1 fraction of the EUSAAR_2 protocol using a thermal–optical analyzer. The mvOC amount agrees very well with the OC1, despite the different heating rate and length of the temperature step, showing the representativeness of the isolated mvOC. The measured ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\left(\mathrm{mvOC}\right)}$ for a combined sample ($\mathrm{0.529}\pm \mathrm{0.007}$), which consists of three individual daily filters, is quite similar to the weighted average of the individual filters. This further validates our isolation method for mvOC that the ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\left(\mathrm{mvOC}\right)}$ is representative of the combined samples. With long sampling duration, low temperature and high face velocities as well as high OC concentrations in winter in this study, the positive artifacts to mvOC should be relatively small and play a limited role in the ${\mathrm{F}}^{\mathrm{14}}\mathrm{C}$${}_{\left(\mathrm{mvOC}\right)}$.

$\mathrm{mvOC}/\mathrm{OC}$ ranges from 3 % to 15 %, with an average of $\mathrm{10.5}\pm \mathrm{3.3}$ %. The ambient samples in this study have a smaller $\mathrm{mvOC}/\mathrm{OC}$ than ambient samples from other places and have a much smaller mvOC fraction than fresh source samples, likely due to the effects from atmospheric aging processes. $\mathrm{mvOC}/\mathrm{OC}$ from fossil sources averages $\mathrm{14}\pm \mathrm{6.6}$ %, consistently higher and more variable than that from nonfossil sources ($\mathrm{7.5}\pm \mathrm{2.9}$ %). That is, OC from fossil sources contains a larger more volatile fraction than that from nonfossil sources.

The nonfossil contribution to mvOC, OC and EC ($f_{\mathrm{nf}}$(mvOC), $f_{\mathrm{nf}}$(OC), $f_{\mathrm{bb}}$(EC)) varies significantly among samples in different sites: the lowest values are always found in Taiyuan and the highest values are found in Chongqing. This implies that different pollution patterns exist in individual Chinese cities. Despite the variations in $f_{\mathrm{nf}}$(mvOC), $f_{\mathrm{nf}}$(OC) and $f_{\mathrm{bb}}$(EC) in different cities, mvOC is more fossil than OC and less fossil than EC in all locations. The fossil contribution to mvOC, OC and EC are, on average, $\mathrm{59}\pm \mathrm{11}$ %, $\mathrm{48}\pm \mathrm{16}$ % and $\mathrm{73}\pm \mathrm{9}$ %, respectively. Higher fossil contribution to mvOC than to OC is consistent with conclusions from ${\mathit{\delta }}^{\mathrm{13}}\mathrm{C}$ of OC desorbed at different temperatures for ambient samples in other locations. $f_{\mathrm{nf}}$(mvOC) correlates better with $f_{\mathrm{nf}}$(OC) ($r^{\mathrm{2}}=\mathrm{0.89}$) than $f_{\mathrm{bb}}$(EC) ($r^{\mathrm{2}}=\mathrm{0.62}$).

SOC concentrations and sources are modeled based on the ${}^{\mathrm{14}}\mathrm{C}$-apportioned OC and EC. SOC concentrations ($\mathrm{15.4}\pm \mathrm{9.0}$ $\mathrm{\mu }$g m${}^{-\mathrm{3}}$) are consistently higher than those of mvOC ($\mathrm{3.3}\pm \mathrm{2.2}$ $\mathrm{\mu }$g m${}^{-\mathrm{3}}$), showing that the contribution of SOC to total OC is considerably larger than that of mvOC. A major fraction of SOC has desorption temperatures above 200 ${}^{\circ }$C and falls into our operationally defined mrOC. mvOC is more fossil than SOC. The relative contribution of fossil sources to mvOC is 59 % on average (range of 45 %–78 %), which is higher than that to SOC (43 %, with a range of 10 %–70 %).

Good correlation is found between mvOC${}_{\mathrm{nf}}$ and OC${}_{\mathrm{nf}}(r^{\mathrm{2}}=$ 0.91). This indicates strong impacts on the variability of mvOC${}_{\mathrm{nf}}$ from nonfossil emissions, including the secondary formation of OC from nonfossil sources and primary biomass burning. The variation in mvOC${}_{\mathrm{fossil}}$ is mainly explained by a good correlation with POC${}_{\mathrm{fossil}}$ ($r^{\mathrm{2}}=\mathrm{0.69}$) and no correlation is found between mvOC${}_{\mathrm{fossil}}$ and SOC${}_{\mathrm{fossil}}$ ($r^{\mathrm{2}}=\mathrm{0.01}$). This indicates that primary mvOC from fossil sources may play an important role in the variation in mvOC${}_{\mathrm{fossil}}$.

Data availability

All data needed to evaluate the conclusions in this study are present in the paper and the Supplement. Additional data related to this study are available upon request to the corresponding author.

The supplement related to this article is available online at: https://doi.org/10.5194/acp-19-10405-2019-supplement.

Author contributions

UD, RJH, HN and JC conducted the study design. HN, HD, WD and JZ carried out the experimental work. HN, UD, AAB and HAJM analyzed data. HN wrote the manuscript. HN, UD and RJH interpreted data and edited the manuscript. HN prepared display items. All authors commented on and discussed the manuscript.

Competing interests

The authors declare that they have no conflict of interest.

Acknowledgements

The authors acknowledge financial support from the Gratama Foundation. Special thanks are given to Dipayan Paul, Marc Bleeker and Henk Been for their help with the AMS measurements at CIO and to Dicky van Zonneveld for her help with ${}^{\mathrm{14}}\mathrm{C}$ data correction at CIO.

Financial support

This research has been supported by the National Key Research and Development Program of China (grant no. 2017YFC0212701), the National Natural Science Foundation of China (grant nos. 91644219 and 41877408), and the KNAW (grant no. 530-5CDP30).

Review statement

This paper was edited by Rupert Holzinger and reviewed by Soenke Szidat and one anonymous referee.