An optimized tracer-based approach for estimating organic carbon emissions from biomass burning in Ulaanbaatar, Mongolia

The impact of biomass burning (BB) on atmospheric particulate matter of < 2.5 μm diameter (PM2.5) at Ulaanbaatar, Mongolia, was investigated using an optimized tracer-based approach during winter and spring 2017. Integrated 24 h PM2.5 samples were collected on quartz-fiber filters using a 30 L min−1 air sampler at an urban site in Ulaanbaatar. The aerosol samples were analyzed for organic carbon (OC) and elemental carbon (EC), anhydrosugars (levoglucosan, mannosan, and galactosan), and water-soluble ions. OC was found to be the predominant species, contributing 64 % and 56 % to the quantified aerosol components in PM2.5 in winter and spring, respectively. BB was identified as a major source of PM2.5, followed by dust and secondary aerosols. Levoglucosan /mannosan and levoglucosan /K ratios indicate that BB in Ulaanbaatar mainly originated from the burning of softwood. Because of the large uncertainty associated with the quantitative estimation of OC emitted from BB (OCBB), a novel approach was developed to optimize the OC / levoglucosan ratio for estimating OCBB. The optimum OC / levoglucosan ratio in Ulaanbaatar was obtained by regression analysis between OCnon-BB (OCtotal– OCBB) and levoglucosan concentrations that gives the lowest coefficient of determination (R2) and slope. The optimum OC / levoglucosan ratio was found to be 27.6 and 18.0 for winter and spring, respectively, and these values were applied in quantifying OCBB. It was found that 68 % and 63 % of the OC were emitted from BB during winter and spring, respectively. This novel approach can also be applied by other researchers to quantify OCBB using their own chemical measurements. In addition to OCBB, sources of OCnon-BB were also investigated through multivariate correlation analysis. It was found that OCnon-BB originated mainly from coal burning, vehicles, and vegetative emissions.


Introduction
Organic aerosol (OA) contributes a significant fraction (10 %-90 %) of atmospheric particulate matter (PM), which can affect human health and air quality (Jimenez et al., 2009;Maenhaut et al., 2011;Fu et al., 2012;Allan et al., 2014;Chen et al., 2018). An understanding of the sources of PM is highly relevant for air-quality remediation. Biomass burning (BB) is a major source of organic carbon (OC) in PM 2.5 (PM with aerodynamic diameter ≤ 2.5 µm), and it may become more significant in the future as air-quality regulations restrict other anthropogenic emissions (Davy et al., 2011;Allan et al., 2014;Sullivan et al., 2019). Coal combustion, thermal power plants, and traffic emissions also make potential contributions to the OC content of PM (Watson et al., 2001a, b;Pei et al., 2016;Deshmukh et al., 2019;, modifying PM characteristics such as hygroscopicity, light-attenuating properties, and health impacts (Jung et al., 2009;Sullivan et al., 2019). Previous studies have observed that the toxicity of PM 2.5 increases with the oxidation potential of BB species because of the water-soluble fraction of OC (Verma et al., 2014).
Ulaanbaatar, with a population of about 1 million, is an atmospheric pollution "hotspot" because of its topography, being situated in the Tuul River valley and surrounded by the Khentii Mountains, with a high elevation (1300-1949 m a.s.l. -above sea level) and large variations in temperature (−28 to +16 • C) and relative humidity (17.7 %-72.7 %; Table 1; Batmunkh et al., 2013;Jung et al., 2014). As the world's coldest capital city during winter, it requires additional fuel for space heating. The topography and low-temperature conditions cause an increase in PM concentrations, which are exacerbated by low wind speeds and atmospheric temperature inversions (Jung et al., 2010).
Half of the residents in Ulaanbaatar live in 160 000 gers (traditional Mongolian dwellings) (Guttikunda and Jawahar, 2014). Biomass is used as fuel for cooking and heating in many low-income gers in Ulaanbaatar. The common tree species in Mongolia are larch, pine, cedar, spruce, and birch; these are mostly softwood (http://www.fao.org/3/w8302e/ w8302e05.htm; http://www.fao.org/3/a-am616e.pdf, last access: 17 December 2019). Each ger burns an average of 3 m 3 of wood per year (Guttikunda, 2008;Zhamsueva et al., 2018). Organic carbon (OC) has severe effects on human health and global climate change (Sun et al., 2019). But there are very few estimates of OC emitted from biomass burning (OC BB ) in Ulaanbaatar. Few studies have investigated the chemical characteristics of aerosol in Ulaanbaatar (Jung et al., 2010;Davy et al., 2011;Batmunkh et al., 2013), with none examining the contribution of OC BB and the type of biomass. Therefore, this study estimated the appropriate In this study, we quantified the BB tracers levoglucosan, mannosan, galactosan, K + , and other chemical species. Potential sources of PM 2.5 were identified by principal component analysis (PCA), with levoglucosan / K + and levoglucosan / mannosan ratios being used to identify major biomass types. OC BB can be quantified from OC / levoglucosan ratios and levoglucosan concentrations in PM. However, uncertainties of OC BB are high because OC / levoglucosan ratios can vary depending on fuel type, burning conditions, and burning location (Duan et al., 2004;Cheng et al., 2013;Jung et al., 2014). Therefore, it is required to determine the most suitable OC / levoglucosan ratio of BB emissions for estimating the appropriate concentration of OC BB . Here, for the first time, optimized OC / levoglucosan ratios were investigated for estimating concentrations of OC BB during winter and spring. OC non-BB sources were also investigated using multivariate correlation analysis with ions and elemental carbon (EC).

Sampling site and aerosol sampling
Aerosol sampling was carried out in Ulaanbaatar during the winter (17 January to 3 February) and spring (17 April to 4 May) of 2017, with 24 h periods commencing daily at 11:00 LT (local time). An aerosol sampler was installed on the rooftop of the National Agency for Meteorology and Environmental Monitoring station in Ulaanbaatar (47 • 92 N, 106 • 90 E, Fig. 1), 10 m a.g.l. (above ground level). The sampling site was located at 8-10 km from two coal-based thermal power plants to the west (Chung and Chon, 2014). PM 2.5 samples were collected on 47 mm diameter quartz-fiber filters (Pall Life Sciences, USA) using an aerosol sampler (Murata Keisokuki Service, Japan) at a flow rate of 30 L min −1 . A field blank filter was collected during winter (n = 1) and spring (n = 1). The quartz-fiber filter was loaded in the sampler for 5 min without operating a pump. The concentration of all chemical analytes has been corrected using the blank filter concentration. Sampled filters were wrapped in aluminum foil and heated at 550 • C for 12 h to remove adsorbed impurities before use and stored at −20 • C before and after sampling.

Filter analysis
One-fourth of each quartz-fiber filter sample was extracted in 10 mL of ultrapure water (resistivity 18.2 M, total OC content < 1 ppb) under ultrasonication for 30 min. The water extract was then filtered using a syringe filter (Millipore, Millex-GV, 0.45 µm) and stored at 4 • C pending analysis. Water-soluble cations (K + , Na + , Ca 2+ , Mg 2+ , and NH + 4 ) were quantified by an ion chromatograph (Dionex ICS 5000, Thermo Fisher Scientific, USA). Water-soluble cations were separated using an IonPac CS-12A column (Thermo Fisher Scientific, USA) with 20 mM of methanesulfonic acid as an eluent at a flow rate of 1.0 mL min −1 . Water-soluble anions (Cl − , NO − 3 , and SO 2− 4 ) were separated using an IonPac AS-15 column (Thermo Fisher Scientific, USA) with 40 mM of KOH as an eluent at a flow rate of 1.2 mL min −1 . The detection limits for major inorganic ions based on 3σ of blanks were 0.01, 0.01, and 0.03 µg m −3 for NO − 3 , SO 2− 4 , and NH + 4 , respectively.
Aerosol samples were analyzed for OC and EC using a thermal optical OC-EC analyzer (Sunset Laboratory Inc. Forest Grove, OR, USA) with laser-transmittance-based correction of pyrolysis. Details of the analyzer and qualitycontrol parameters are reported elsewhere (Jung et al., 2014). In brief, 1.5 cm 2 punch samples of the quartz-fiber filter were placed in a quartz dish inside the thermal desorption oven of the analyzer. OC and EC were quantified using a temperature program developed by the US National Institute for Occupational Safety and Health (NIOSH) in an inert atmosphere (100% He) and in an oxidizing atmosphere (98 % He + 2 % O 2 ), respectively. Detection limits of OC and EC were 0.04 and 0.01 µg C m −3 , and their analytical uncertainties were 1.3 % and 3.7 %, respectively.

Conditional probability function
The conditional probability function (CPF) calculates the probability that a source is located within a particular wind direction sector, : where n is the number of times that the wind passed through direction sector , and m is the number of times that the source contribution peaked while the wind passed through sector (Ashbaugh et al., 1985). To use CPF with the Ulaanbaatar data, the 24 h averaged source contribution data have been applied to all 1 h wind direction averages recorded at the site for each date. The angular interval was set at 10 • . To calculate m , the 75th percentiles of source contribution concentrations were counted. CPF is useful in determining the direction of a source from a receptor site; however, it cannot determine the actual location of the source.

Principal component analysis
In order to identify the source groupings of chemical species in PM 2.5 , principal component analysis (PCA) was applied. PCA is done using a commercially available software package (SPSS, version 10.0). PCA applies projection dimension reduction methods, converting several concentration sets into significant sets of columns (principal components, PCs) without damaging the original data. PCA is a widely used statistical technique to quantitatively identify a small number of independent factors among the species concentrations, which can explain the variance of the data, by using the eigenvector decomposition of a matrix of pair-wise correlations. PCA with varimax rotation and the retention of principal components having eigenvalues > 1.0 was used to identify major species associated with different sources. It was widely used for the identification of pollution sources in the atmosphere (Fang et al., 2003;Nirmalkar et al., 2015).

Chemical characteristics of PM 2.5 and source identification
Mass concentrations of carbonaceous aerosol, BB tracers, and water-soluble ions in PM 2.5 samples collected at Ulaanbaatar during winter and spring of 2017 are summarized in Table 1. OC contributed 64 ± 5.1 % and 56 ± 6.0 % of the quantified aerosol components in PM 2.5 in winter and spring, respectively (Table 1). Average concentrations of OC during winter were 5 times those obtained in spring (Fig. 2). Previously, OC has been observed as a major component in PM 2.5 in Ulaanbaatar during the winter period (Jung et al., 2010;Batmunkh et al., 2013). This may be attributed to additional BB emissions for home heating and temperature inversions with low wind speeds (average wind speed of 1.43 ± 0.73 m s −1 ; Table 1 and Fig. 3a). OC concentrations decreased with increasing wind speed during winter (Fig. 3a) and spring (Fig. 3b) over all air temperature ranges. The inverse relationship between OC and wind speed during winter ( Fig. 3a) and spring (Fig. 3b) suggests a predominance of local sources, with higher wind speeds flushing air pollutants out of the area, whereas low wind speeds allow them to accumulate (Khan et al., 2010;Wang et al., 2018). The average concentration of EC during winter (1.71 ± 0.58 µg m −3 ) was higher than that in spring (1.11 ± 0.42 µg m −3 ) (Table 1), consistent with general urban observations in cities in China  and India (Panda et al., 2016). During both winter and spring, EC concentrations at the study site were lower and had different trends compared to those observed in a suburban site (2.3±1.0 and 3.1± 1.5 µg m −3 , respectively) and an urban site (2.3 ± 1.0 and 3.3 ± 1.2 µg m −3 , respectively) in Shanghai, China (Feng et al., 2009).
The potential source direction of EC during winter and spring was west as shown in Fig. 5; this can be explained by the influence of emissions from thermal power plants. The correlation of EC was strong with Ca 2+ during spring as shown in Fig. 4. CPF analysis suggested that the poten- tial source direction of EC and Ca 2+ was similar (Fig. 5). High abundances of Ca 2+ and EC are observed from stack emissions of coal-fired thermal power plants (Pei et al., 2016;Zhang et al., 2015). Thus, EC and Ca 2+ in Ulaanbaatar might be strongly related to emissions from thermal power plants.
Daily concentrations of levoglucosan, mannosan, and galactosan have similar trends during winter and spring ( Fig. 2), possibly because of the combustion of similar biomass fuels in both seasons. Changes in the concentrations of these BB tracers might be attributed to changes in the relative proportions of cellulose and hemicellulose in different biomass fuels (Zhu et al., 2015;Nirmalkar et al., 2015). Concentrations of anhydrosugars were 4 times higher in winter than in spring (Table 1)    heating requirements in winter. The higher relative humidity (58.5 %-72.7 %) and lower temperature (−10.5 to −27.8 • C; Table 1) in winter can also contribute to longer atmospheric residence times due to increased levoglucosan stability (Lai et al., 2014). Higher concentrations of BB tracers in winter than spring have previously been observed in Beijing, China (Liang et al., 2016), and were attributed to meteorological conditions similar to those in Ulaanbaatar. Further, it was observed that during winter the ambient temperature was consistently low (less than −10 • C; Fig. 3a) in Ulaanbaatar; therefore, residential biomass burning occurred continuously for space heating. Thus, there would be no dependence of emission rates of levoglucosan and mannosan on ambient temperature during the sampling period. The concentration of levoglucosan and mannosan may be influenced by local wind speed, as shown in Fig. 3a, rather than ambient temperature. The average K + concentration (0.08 ± 0.05 µg m −3 ) in this study is significantly lower than the K (0.32 µg m −3 ) observed in Ulaanbaatar during -2008(Davy et al., 2011. Among water-soluble ions, SO 2− 4 (9.7 ± 3.4 µg m −3 ) was the most dominant PM 2.5 species during winter, followed by NH + 4 (6.2 ± 2.4 µg m −3 ) and NO − 3 (4.2 ± 1.7 µg m −3 ), whereas SO 2− 4 (1.9 ± 0.5 µg m −3 ) was the dominant species during spring, followed by Ca 2+ (0.9±0.4 µg m −3 ) and NH + 4 (0.7 ± 0.3 µg m −3 ). The total SO 2− 4 + NH + 4 + NO − 3 content accounted for 27 % and 23 % of the total measured chemical species during winter and spring, respectively ( Fig. 2 and Table 1). SO 2− 4 is the most prevalent water-soluble ion in PM 2.5 in Wuhan, Guangzhou, and Tianjin (China) due to industrial emissions and coal burning (Gu et al., 2011;Tao et al., 2014;Huang et al., 2016;Pei et al., 2016). This suggests that the higher SO 2− 4 concentration in Ulaanbaatar may be attributable to emissions from the three major coal-fired thermal power plants near the study site.
The atmospheric concentrations of OC (11-17 µg m −3 ) and levoglucosan (0.46-0.73 µg m −3 ) were higher for samples collected during 27-30 April 2017 than on almost all remaining days in spring (Fig. 2b). Backward atmospheric trajectories based on the Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model provided by the US National Oceanic and Atmospheric Administration (NOAA) Air Resources Laboratory (ARL) indicate that during those days, air masses originated from a region where a significant number of fires were detected (US Fire Information for Resource Management System -FIRMS; National Aeronautics and Space Administration -NASA; Fig. 6a and b). Thus, the elevated OC and levoglucosan concentrations during 27-30 April might be influenced by the longrange transport of BB from north of Mongolia.

Principal component analysis
Principal component analysis (PCA) is a useful tool for reducing the dimensionality of large aerosol datasets to principal components using varimax rotation for source identification (Cao et al., 2005;Lin et al., 2018;Nirmalkar et al., 2019). Four principal components (PCs) in winter and three in spring were identified with eigenvalues > 1 after varimax rotation, explaining 96 % and 92 %, respectively, of the total variance (Tables 2 and 3). The PCs were categorized on the basis of the loadings of chemical components as follows. In winter, PC1 includes BB characterized by high loadings of levoglucosan, mannosan, and galactosan; PC2 includes dust characterized by Ca 2+ and Mg 2+ content; PC3 includes secondary formation characterized by SO 2− 4 , NO − 3 , and NH + 4 content; and PC4 includes fossil fuel combustion characterized by EC. In spring, PC1 includes BB (levoglucosan, mannosan, and galactosan); PC2 includes dust (Ca 2+ and Mg 2+ ) and fossil fuel combustion (EC); and PC3 includes secondary formation (SO 2− 4 , NO − 3 , and NH + 4 ). The PCA results show that the chemical components of PM 2.5 in Ulaanbaatar were mainly affected by BB during winter and spring. Further, OC was primarily influenced by BB because it correlated well with the total variance of PC1 during winter (0.82; Table 2) and spring (0.77; Table 3).

Relationship among BB tracers
The correlations among the three BB tracers levoglucosan, mannosan, and galactosan are shown in Fig. 7a (winter) and Fig. 7b (spring). The correlations between levoglucosan and mannosan and between levoglucosan and galactosan are strong during winter (R 2 = 0.99 for both pairs) and spring (R 2 = 0.95 and 0.83, respectively; Fig. 7a and b).
Concentrations of levoglucosan and OC are strongly correlated during both winter (R 2 = 0.78) and spring (R 2 = 0.86; Fig. 8a), suggesting that a major fraction of OC might originate from BB in Ulaanbaatar. The similar strong correlation and steep slope observed in OC-levoglucosan plots for PM collected in Chiang Mai Province (Thailand) and Daejeon (Korea) were attributed mainly to BB (Jung et al., 2014;Thepnuan et al., 2019).
Fine-mode K + is considered to be a biomass burning tracer in previous studies (Louie et al., 2005;Deshmukh et al., 2011;Cheng et al., 2013). The moderate correlation between levoglucosan and K + concentrations (R 2 = 0.68) in winter indicates that they are produced from similar sources (Fig. 8b), with BB contributing most of the K + . However, the correlation between levoglucosan and K + was weak in spring (R 2 = 0.49; Fig. 8b). Because K + is typically emitted at a higher mass fraction in flaming-phase combustion compared to smoldering (Lee et al., 2010), smoldering combustion tends to have a higher levoglucosan / K + emission ratio compared to flaming combustion (Schkolnik et al., 2005;Gao et al., 2003). A high levoglucosan / K + ratio was observed during winter (8.92) compared to spring (4.21) at this site. Thus, a weak correlation between levoglucosan and K + concentrations at Ulaanbaatar in spring can be explained by mixed burning conditions such as smoldering and flaming.       Fig. 9a) and spring (R 2 = 0.73; Fig. 9b), suggesting that they might originate from similar sources. Because most of the aerosol particles emitted from BB belong to PM 2.5 , the correlation between OC and K + as well as levoglucosan suggests that BB is one of the potential sources of OC in winter and spring. Because biomass fuel is burned in traditional stoves with no pollution control devices in Ulaanbaatar (Batmunkh et al., 2013), soil and ash particles are entrained in convective processes and uplifted in the atmosphere together with smoke particles (Deshmukh et al., 2011;Nirmalkar et al., 2019).

Tracing the source of BB aerosol
OC is a major contributor of the quantified aerosol components in PM 2.5 in Ulaanbaatar during spring and winter (Table 1). To quantify the OC BB , it is necessary to identify the BB fuel type. Several investigators used levoglucosan / mannosan and levoglucosan / K + ratios to identify BB fuel types (Puxbaum et al., 2007;Cheng et al., 2013;Jung et al., 2014;Chen et al., 2018;Thepnuan et al., 2019).
The levoglucosan / mannosan ratio is source-specific and can be used to identify BB fuel types due to the unique cellulose and hemicellulose compositions of different biomass fuels (Zhang et al., 2007;Cheng et al., 2013). A previous study suggested that the levoglucosan / mannosan ratio is strongly dependent on wood type, rather than on the site where the wood is grown (Cheng et al., 2013). Therefore, the levoglucosan / mannosan ratio was used to trace the type of wood burned during winter and spring for indoor heating and cooking purposes. Previous studies have used levoglucosan / mannosan ratios to investigate the BB fuel types (Cheng et al., 2013;Jung et al., 2014).
Therefore, softwood burning seems to be the major source of BB aerosol in Ulaanbaatar during both winter and spring, consistent with previously reported softwood-burning emissions from fireplaces in northern and southern regions of the USA (Fine et al., 2001, from household combustion in Zhengzhou, China , and from stovewood combustion in the mid-European region (Austria; Schmidl et al., 2008).   was used for estimating OC BB at four background sites in Europe (Puxbaum et al., 2007). Later, a mean value of 11.2 for the OC / levoglucosan ratio derived from ratios ranging 4.5-24.6 was used for estimating OC BB in the UK (Harrison et al., 2012). However, such estimates may not be accurate as the OC / levoglucosan ratio is highly variable in BB emissions. For example, the average OC / levoglucosan ratio from softwood burning (23.8) is much higher than that of hardwood burning (7.35) Schmidl et al., 2008); differences are more than tenfold among studies of softwood-burning OC / levoglucosan ratios Hays et al., 2002;Engling et al., 2006;Iinuma et al., 2007;Gonçalves et al., 2010). Combustion conditions may also significantly influence OC / levoglucosan ratios. For example, the OC / levoglucosan ratio varied by a factor of about 7 between burning the same wood (Loblolly pine) in a fireplace (27.6; and in a stove (3.4; Fine et al., 2004b). Therefore, it is necessary to optimize the OC / levoglucosan ratio for use in estimating OC BB . This study has used an optimized OC / levoglucosan ratio to estimate the precise concentration of OC BB for the Ulaanbaatar study site. We have used a range of different OC / levoglucosan ratios obtained from previous literature (Fig. 11) for a regression analysis with measured levoglucosan concentrations to estimate the optimum OC / levoglucosan ratio ( Fig. 12a and b). First, candidate OC BB (Fig. 11) in this study was estimated from OC / levoglucosan ratios for softwood burning in previous chamber experiments (Cheng et al., 2013;Schauer et al., 2001;Hays et al., 2002;Fine et al., 2001Fine et al., , 2004aEngling et al., 2006;Iinuma et al., 2007;Schmidl et al., 2008;Gonçalves et al., 2010, Fig. 11) and the measured levoglucosan concentration at this site. Second, the OC non-BB concentration was calculated by subtracting OC BB from the corresponding total OC. If calculated OC non-BB does not contain OC BB , both the regression slope and R 2 between OC non-BB versus levoglucosan will be close to zero. As shown in Fig. 12a and b, the lowest R 2 and regression slope were observed with OC / levoglucosan ratios of 27.6 and 18.0 in winter and spring, respectively. Thus, the optimized OC / levoglucosan ratios for our site were determined to be 27.6 and 18.0 in winter and spring, respectively.

Optimization of OC
During winter the higher optimum ratio of OC / levoglucosan might be due to incomplete combustion during smoldering phenomena. As smoldering fires are characterized by lower temperatures and thus have lower combustion efficiency, they release more un-combusted condensable products, resulting in the production of more unbroken organic compounds (Engling et al., 2006). Smoldering combustion generally leads to increased emissions of volatile organic compounds (VOCs) and particulate organic matter (OM) (Obrist et al., 2007). In contrast, the relatively lower optimum ratio of OC / levoglucosan during spring might be due to the higher combustion efficiency during flaming phenomena.  The OC BB concentrations at the Ulaanbaatar study site were calculated from the optimized OC / levoglucosan ratios (winter: 27.6 and spring: 18.0) and measured levoglucosan concentrations. The OC BB concentration was estimated to be 33.1±11.9 µg C m −3 (range 16.0-58.5 µg C m −3 ) and 5.64± 3.29 µg C m −3 (range 0.57-13.1 µg C m −3 ), accounting for 68 % and 63 % of the total OC in winter and spring, respectively (Fig. 13). The average of previously published OC / levoglucosan ratios, 10.1±7.9 (range 1.90-27.6), gives an estimated OC BB concentration of 12.1 ± 4.4 µg C m −3 (range 5.9-21.4 µg C m −3 ) and 3.2 ± 1.8 µg C m −3 (0.32-7.34 µg C m −3 ) in winter and spring, respectively. Their values are 2.7 (winter) and 1.8 (spring) times lower than values estimated using our optimized OC / levoglucosan ratio.
Our estimated contribution of OC BB was higher than that in Daejeon, South Korea (24 %-68 % of total OC, mean 45 % ± 12 %; Jung et al., 2014), and Beijing, China (50 % of total OC; Cheng et al., 2013), where BB aerosols are produced mainly by the burning of crop residues. The contribution of OC BB to total OC is 57 % and 31 % during the heating (average temperature 0.6 • C) and non-heating (average temperature 14 • C) seasons in Krynica-Zdrój, Poland (Klejnowski et al., 2017), significantly lower than that of Ulaanbaatar during both winter (average temperature −21 • C) and spring (average temperature 6 • C). Such high concentrations of OC BB in Ulaanbaatar and Krynica-Zdrój are likely due to intense wood burning for heating during winter.

Tracing sources of OC non-BB
A high concentration of OC non-BB was found during winter compared to spring (Fig. 13). Elevated OC non-BB could be attributed to enhanced emissions from combustion and favorable meteorological conditions (cold temperatures and inversion conditions, etc.) during the winter. There is a strong correlation between OC non-BB and SO 2− 4 , NH + 4 , and K + in winter and OC non-BB and NO − 3 , Na + , K + , Mg 2+ , Ca 2+ , and EC in spring (Table 4). The residential combustion of coal emits significant amounts of OC, EC, and inorganic species (SO 2− 4 and metals) due to incomplete combustion and lack of pollution control devices (Garcia et al., 1992;Li et al., 2016;Watson et al., 2001a, b). Garcia et al. (1992) studied emissions of volatile organic compounds from coal burning and vehicle engines.
In Ulaanbaatar, the use of wood and coal for cooking and heating, as well as emissions from old vehicles are reported as potential sources of OC (Batmunkh et al., 2013;Zhamsueva et al., 2018). The three thermal power plants in Ulaanbaatar are point sources for emissions of carbonaceous aerosol (Batmunkh et al., 2013), burning ∼ 5 million tons of coal per year (Batmunkh et al., 2013). High concentrations of anions (SO 2− 4 and NO − 3 ) and cations (NH + 4 and Na + ) are reported in China (Zhou et al., 2003), the USA (Caiazzo et al., 2013), Brazil (Flues et al., 2002, India (Guttikunda and Jawahar, 2014), Korea (Park and Kim, 2004;Park et al., 2015), and Spain (Alastuey et al., 1999) near coal-fired ther-mal power plants. Emissions of volatile organic compounds from vegetation have also been observed in previous studies (Fehsenfeld et al., 1992;Shao et al., 2001;Acton et al., 2016). The correlations of OC non-BB with ions and EC are thus likely due to volatile organic compounds emitted from coal burning, vehicles, and vegetative emissions.

Conclusions
BB was identified as a major source of the quantified aerosol components in PM 2.5 in Ulaanbaatar, Mongolia, during the winter and spring of 2017 based on PCA. OC was the major component of the quantified aerosol components during the entire sampling period for the winter and spring. For the determination of OC BB , the fuel type must be identified, and levoglucosan / mannosan and levoglucosan / K + ratios obtained from previous studies and our on-site measurements were used for this purpose.
Softwood burning was identified as a major source of OC BB . However, OC / levoglucosan ratios from softwood burning are highly variable, and an optimum ratio was derived by regression analysis between daily concentrations of OC non-BB and levoglucosan, yielding values of 27.6 and 18.0 for winter and spring, respectively. The application of these ratios indicates that 68 % and 63 % of the OC originated from BB during winter and spring, respectively, which is about double that estimated using average values of previous studies. The atmospheric concentration of OC BB was higher in winter than in spring, mainly due to additional BB for heating and cooking. BB aerosols in Ulaanbaatar originate mainly from local softwood burning. The approach developed here may be applied elsewhere to screen regionspecific OC / levoglucosan ratios for estimating appropriate atmospheric concentrations of OC BB , aiding the establishment of BB control measures.
Data availability. The data used in this study are available from the corresponding author upon request (jsjung@kriss.re.kr).
Author contributions. JJ and TB designed the study and carried out the fieldwork. JJ performed chemical analyses and quality-control measures. JN wrote the paper under the guidance of JJ. All authors commented on and discussed the paper.