Characterization and Source Identification of Elements and Water ‐ Soluble Ions in Submicrometre Aerosols in Brno and Šlapanice ( Czech Republic )

Submicrometre aerosol particles  (particulate matter, PM1) were collected  in  two Czech cities (Brno and Šlapanice) during week campaigns in winter and summer of 2009 and 2010. The aerosols  were  analysed  for  14  elements  and  12  water‐soluble  ions  using  inductively  coupled plasma–mass  spectrometry  and  ion  chromatography  techniques.  The  average  PM1  mass concentration was 14.4 and 20.4 μg m−3 in Brno and Šlapanice, respectively. Most of the analysed elements  and  ions  exhibit  distinct  seasonal  variability with  higher  concentrations  in winter  in comparison  to summer. The determined elements and  ions  together accounted  for about 29% of total PM1 mass, ranging between 16% and 44%. Ion species were the most abundant components in collected aerosols, accounting for 27.2% of mass of PM1 aerosols, and elements accounted for 1.8% of mass of PM1 aerosols. One‐day backward trajectories were calculated using the Hysplit model to analyse  air masses  transported  towards  the  sampling  sites. The Pearson  correlation  coefficients between  individual  PM1  components  and  PM1 mass  and  air  temperature were  calculated.  To identify the main aerosol sources, factor analysis was applied. Six factors were identified for each locality. The  following sources of PM1 particles were  identified  in Brno: a municipal  incinerator, vehicle  exhausts,  secondary  sulphate,  a  cement  factory,  industry  and  biomass  burning.  The identified sources  in Šlapanice were as  follows: a combustion source, coal combustion, a cement factory, a municipal incinerator, vehicle exhausts and industry.


Introduction
Atmospheric aerosol (particulate matter, PM) is an important airborne component with various environmental and health effects. Aerosols can deteriorate air quality, affect global climate, reduce visibility, and are involved in smog production [1]. Epidemiological studies showed an association between the concentration of aerosol particles in ambient air and adverse health effects [2][3][4]. The environmental and health effects of atmospheric aerosols depend on the particle size, shape and the chemical composition of particles. The atmospheric aerosol consists of a complex mixture of components including carbonaceous species (elemental and organic carbons), inorganic ions and elements in variable amounts, depending on their location and emission sources.
Receptor models based on the statistical evaluation of PM chemical data acquired at receptor sites are used to identify emission sources of aerosols [29]. The chemical mass balance model assumes knowledge of the composition of the emissions for all relevant sources, however, fulfilling this requirement is often problematic [29]. More widespread are two multivariate models, principal component analysis [30][31][32][33] and positive matrix factorisation [5,8,22,24,[34][35][36], which apportion the sources on the basis of the ambient data from the receptor site alone.
The aims of the study were to fill the gap in the missing information about the content of elements and water-soluble ions in aerosols in the Brno agglomeration; to obtain comprehensive information on the composition of PM1 aerosols in Brno and Šlapanice in a combination with other studies in the same area focused on the characterization of organic compounds in PM1 aerosols in Brno and Šlapanice [11,18,19]; to identify their sources. The paper presents the results of determination of elements and water-soluble ions in submicrometre aerosols (PM1, particles with aerodynamic diameter smaller than 1 μm) collected in two cities in the Czech Republic (Brno and Šlapanice) in Central Europe. The seasonal differences (winter vs. summer) were evaluated and the sources of studied aerosol components were analysed.

Sampling Sites
The PM1 aerosols were sampled in Brno and Šlapanice that represent a large city and a small town in the Czech Republic. Brno, the second largest city in the Czech Republic (370,000 inhabitants), is an industrial and administration centre of Moravia, the eastern part of the Czech Republic. There are various local sources of aerosols in Brno, such as traffic (cars and trams on Veveří street), residential heating and large emission sources, such as heating plant, a municipal waste incinerator (on the eastern outskirts of Brno) and industry, including a foundry plant. Regional sources comprise mainly residential heating in the surrounding villages and a cement factory east of Brno. Šlapanice, a small town (6000 inhabitants), is located 3 km southeast from Brno. Overall, the sources of aerosols in Šlapanice are similar to those in Brno, but there is a difference between the composition of local and regional sources in both locations. Local sources of aerosols in Šlapanice include traffic, residential heating, small industrial factories and brickworks. Aerosols can be transported to Šlapanice from various regional sources, such as the nearby motorway between Brno and Šlapanice, Brno airport southwest of Šlapanice, a municipal waste incinerator northwest of Šlapanice, residential heating in the surrounding villages, a cement factory northeast of Šlapanice, etc. In addition, a large power plant in Hodonín that burns coal and biomass is located about 50 km southeast from Brno and Šlapanice. Moreover, a long-range transport of pollutants from distant areas or neighbouring countries to Brno and Šlapanice cannot be ignored.

Aerosol Sampling
Atmospheric aerosols in the size fraction of PM1 were sampled for 24 h every day over one week in winter and one week in summer of 2009 and 2010 in Brno and Šlapanice (Figure 1) to compare the PM1 composition in the large city and a nearby small town. The sampling of aerosols began every day at 9:00 a.m. Aerosol samples in Brno were collected in an urban locality on the balcony on the first floor (at the height of 8.9 m above ground level and at the distance of 15.6 m from the street Veveří) of the Institute of Analytical Chemistry facing northeast toward the street Veveří (49°12'28.27´´N and 16°35´28.00´´E). In Šlapanice the aerosols were collected in small urban locality in the garden of a family house (49°09´55.92´´N and 16°43´26.18´´E). Sampling locations are located relatively in the centre of both Brno and Šlapanice, and therefore, the influence of different PM sources is expected, although the sampling site in Brno is situated near the street with traffic.
Submicrometre aerosols were collected at each site in parallel using a high-volume (HV) and a low-volume (LV) sampler. The HV sampler (DHA-80, Digitel, 30 m 3 h −1 ) equipped with a PM1 size selective impaction inlet (model DPM01/30/00, Digitel) collected PM1 aerosols on cellulose nitrate filters (150 mm diameter, porosity 3 μm, Sartorius). A total number of 52 samples (24 samples from Brno and 28 samples from Šlapanice) were collected with HV sampler during all campaigns. The LV sampler (1 m 3 h −1 ), consisting of a Teflon coated aluminium cyclone inlet (cut point diameter of 1 μm, model URG-2000-30EHB), and a NILU filter unit (type 9633) collected PM1 aerosols on 47mm Teflon filters (Zefluor, porosity 1 μm, PALL). To eliminate interference of gaseous pollutants, such as SO2, HNO3, NH3 and others, an annular diffusion denuder [37] was placed between the cyclone and Teflon filter. A total number of 56 samples (28 samples from Brno and 28 samples from Šlapanice) were collected with the LV sampler during all campaigns. Meteorological parameters (i.e., temperature and relative humidity) measured by means of a commercial sensor (type T3113, Comet Systems) are given in Table 1. We also present predominant wind directions obtained from one-day backward trajectories calculated using the Hysplit v5.0.0 model. All trajectories are shown in the Supplementary Material (Figure S1-S8) separately for each location and campaign.

Processing and Filter Analysis
Mass concentrations of collected aerosols were determined by weighing filters, using a microbalance M5P (±1 μg; Sartorius). Filters were equilibrated before weighing in an air-conditioned room under constant conditions for 48 hrs (temperature 20 ± 1 °C, relative humidity 50 ± 3%). Static electricity was eliminated with an ionizer prior to weighing (PRX-U, Haug). After weighing, exposed filters were cut, using ceramic scissors, into two equal pieces; each of them was weighted again.
One half of the cellulose-nitrate filters was digested in 4 mL of sub-boiling nitric acid in the UniClever microwave device (Plazmatronika). The decomposed samples were transferred quantitatively along with 4 mL of deionized water into polyethylene scintillation vials (Kartel). The extracts were analysed for the content of 14 selected elements (Al, K, Ca, Fe, Mn, Zn, Cu, Cd, Ba, As, Pb, V, Ni, Sb), employing an inductively coupled plasma-mass spectrometry (model 7500 CE, Agilent). Relative uncertainty of element analysis was in the range of 1 to 3%.
Both halves of the Teflon filters were extracted in 8 mL of deionized water under ultrasonic agitation. The extract of the first half was analysed for the content of seven anions (fluoride, chloride, nitrite, nitrate, sulphate, oxalate, phosphate), while the second half of the filter was analysed for five cations (Na + , K + , NH4 + , Ca 2+ , Mg 2+ ) by means of ion chromatography (ICS-2100, Dionex). Relative uncertainty of water-soluble ion analysis was in the range of 4 to 7%.

Calculation of Air Trajectories
The analysis of air mass transported towards the sampling sites was performed using the Hysplit model v5.0.0 [38,39], by calculation of one-day backward trajectories at 300, 750 and 1500 m above ground level.

Factor Analysis Model
Measured aerosol components can be grouped by their correlations. The components within a particular group are highly correlated among themselves but have relatively small correlations with aerosol components in a different group. In the factor analysis presented in this paper, it is assumed the observed correlations in a given group of aerosol components are influenced by a factor, which for given group corresponds to a source of pollution (such a pollution source can be, for example, vehicle exhausts, coal combustion and others). Factors are unmeasurable and the number of factors is much smaller than the number of aerosol components. The correlation relationship among aerosol components can be described in terms of a few underlying factors.
Suppose the observed concentrations of elements and ions in PM particles with means (expectations) . In our study Let , ) be considered common factors corresponding to the sources of pollution. We then consider the linear model describing how the concentrations of elements and ions can be explained using linear combinations of these common factors . The model can be described by equation (1) where , are the coefficients of linear combinations. They represent the so-called factor loadings indicating how strong is the statistical association between the i-th aerosol component and the j-th factor . The higher the absolute value of is, the higher is the statistical association between the aerosol component and the factor The term stands for random errors. Random errors are supposed to be independent with zero mean and variance , and common factors are assumed to have a zero mean and unit variance. The independence of common factors and random errors ( ) is also considered. From the above assumptions follows Var(Xi) = li1 2 + li2 2 +…+ lim 2 + ψi = hi 2 + ψi, where hi 2 = li1 2 + li2 2 +…+ lim 2 is called the i-th communality, the proportion of the variance of the variable Xi contributed by the m common factors. Various methods [40] that allow us to estimate the parameters of model exist. Here, the principal component method [41], which does not impose restrictions on distribution of analysed variables, is preferred. A crucial task of the factor analysis is the selection of appropriate number of common factors . Since the higher values of lead to common factors that are difficult to interpret, the aim is to choose the smallest possible value , such that a sufficient proportion of variability in the data is explained. For this problem, Kaiser´s criterion and scree plot [41] are considered.
The original factor loads obtained by the principal component method can be difficult to interpret. Thus, the factor rotation is performed for better interpretation of factor loads. A common goal of rotation in this paper is to ensure that each aerosol component loads highly on a single factor and has small-to-moderate loads on the remaining single factors. The factor model with such rotated loads is called the model with a simpler structure. Usually the criterion, which is used to obtain a simpler structure, is the varimax criterion [41].

Characterization of Submicrometre Aerosols
Mass concentrations of PM1 aerosols collected at both sites during all campaigns are given in Table 2 Šlapanice (no snow during the campaign), the PM1 concentration remained at a normal level for this time of a year. The concentrations of aerosols in winter were higher than in summer both in Brno and Šlapanice, which can be attributed to increased anthropogenic emissions in heating season due to household heating of local residents [11,18,19] and to lower mixing layer height due to the lower dispersive capacity of the lower atmospheric layers in winter, which favours the accumulation of pollutants and prevents air convention and the dispersion of pollutants [42,43]. The strong relationship of the PM1 concentration with air temperature was confirmed by a significant negative correlation (p < 0.01) both in Šlapanice (a correlation coefficient R = -0.71, Table S1) and in Brno (−0.64, Table S2). Calculated one-day backward trajectories (Figures S1-S8) indicate possible transport of not only regional air pollution from surrounding villages but also long-range transport of polluted air from more distant localities, such as heavy polluted areas in the Silesian Voivodeship in southern Poland that has been recently identified as an important PM source in winter [12]. The differences in the concentration of PM1 aerosols in corresponding seasons in 2009 and 2010 (Table 2) were probably caused mainly by different meteorological and dispersal conditions in those two years.
The mass of determined elements and ions together accounted for about 29% of total PM1 mass, ranging between 16% and 44%. The rest of the mass probably consists of organic compounds, elemental carbon and water. Previous studies at the same locations showed that organic material and elemental carbon formed on average 37.6% and 7.80% of PM1 mass [11,18,19].  The concentration of PM1 aerosols found in Brno and Šlapanice in 2009 and 2010 are similar to those observed in other study from this site [11] as well as in PM1 from other localities in the Czech Republic [6,23], or from other sites in Europe, such as Birmingham [44], Melpitz [45], Bologna [46], Granada [34], Katowice [47] and Racibórz, whereas the concentration of PM1 in Zabrze [48] was higher than in Brno and Šlapanice.
Two recent studies [21,49] from the same location in Brno do not show any decrease in air pollution with PM1 aerosol. The concentration of PM1 aerosols collected in summer 2014 (i.e., the mean value of 10.8 μg/m 3 ) was similar to the concentration found in the presented study, while the concentration found in winter 2015 (i.e., the mean value of 16.5 μg/m 3 ) [21] was in the middle of the concentrations found in the present study in winter 2009 and 2010. The sampling of PM1 aerosols in winter 2017 was accompanied by a few days of smog event, which resulted in an increased PM concentration in Brno and the concentrations of PM1 aerosols measured in this study (i.e., the mean value of 34.2 μg/m 3 ) [49] was thus even higher than the concentrations found in the present study.

Characterization of Elements
We analysed 14 elements in all submicrometre aerosol samples collected in Brno and Šlapanice. The average element concentrations in PM1 aerosols from both sampling sites are summarized in Table 3 (campaigns of 2009) and Table 4 (campaigns of 2010). In winter 2009, the concentrations of elements in Šlapanice were higher than those in Brno, whereas in all other seasons the concentrations of elements at both localities were comparable. The sum of concentrations of the analysed elements accounted for 1.77% of PM1 mass. In winter, the contribution of elements to PM1 mass was 2.17% (1.70-2.76%) and in summer decreased to 1.37% (0.92-1.97%). In winter, lead and potassium were the two most abundant elements accounting for 65-84% of the total element mass in PM1, while in summer potassium prevailed (27-40% of the total element mass in PM1). The daily changes in element concentrations in 2009 and 2010 are shown in Figures 2 and 3, respectively. As, Cd, Pb and Ni are elements known to be toxic to human health [50,51]. The annual average concentration of these elements, calculated as the average of summer and winter concentrations, do The concentrations of elements in the PM1 samples from both the Brno and Šlapanice sites are comparable with those found in Granada [34], Bologna [46] and Frankfurt [53], but lower than those in samples from other European sites, such as Katowice [47] or Tito Scalo [54]. In other Polish cities, Zabrze and Racibórz, the concentrations of several elements (Al, Mn, Fe, Cu, Zn, As) in PM1 were higher than those in Brno and Šlapanice, while the concentrations of other elements were comparable [48].   The enrichment factors (EFs) of elements in aerosols collected in Brno and Šlapanice were calculated to discriminate the anthropogenic and crustal (i.e., natural) origin of the studied elements. The EFs are defined by the equation EF = (X/R)aerosol/(X/R)crust (2) where X represents the considered element and R is the reference element, while the subscripts aerosol and crust indicate concentrations in PM1 and in the Earth's crust, respectively [55]. The EFs were calculated for Fe as the reference element [21]. Generally, EF < 5 indicates the crustal soil as the predominant source of the element, while the anthropogenic origin of the element may be considered at EF > 100 [56]. The enrichment factors calculated from the element concentrations determined in PM1 aerosols collected in Brno and Šlapanice are shown in Table 5. Pb, Cd and Sb have the highest EFs, with Pb being the most enriched element for PM1 particles both in Brno and Šlapanice followed by Cd, Sb, Zn and As with EF > 100, which indicates their anthropogenic sources. The EFs of Cu, Ba, Ni, K, Mn and V are within the range of 5-100, which suggests that these elements are of both anthropogenic origin and soil contribution. The EF of Al approaches to unity, which indicates crustal soil as a predominant source [54,57]. The EFs of most elements in the aerosols sampled in Šlapanice are larger than enrichment factors of elements in aerosols collected in Brno. The obtained value of the EFs from both localities are similar to those of the recent study in Brno [21], with the exception of V, Mn and Cd, the values of which in this study are approximately twice higher, and the EF of Pb is approximately 20 times higher. Most of the elements exhibit distinct seasonal variability when the concentrations of elements in winter were higher than in summer. High seasonal differences were observed especially for K, Pb, Fe, As, V, Cd, Zn, Ni, Sb and Ba, and less for Mn and Cu. Potassium is produced largely during wood combustion [58]. Other metals can originate from several different sources. Arsenic and cadmium originate mainly from coal combustion and partly from industry. Zn, Mn, Cu, Fe or Pb are produced, next to various industrial sources and traffic, mainly by biomass (wood) burning and coal combustion [59][60][61][62]. Coal combustion and biomass burning are still frequently used for energy production and for residential heating in many European countries [36,63]. Hence, it is plausible that increased concentration of these elements (with exception of Al and Ca) in heating seasons when compared to the rest of the year is largely caused by wood and coal combustion in residential heating. Recent studies [11,15,18,19] and the results from the last census in the Czech Republic in 2011 [64] indicate that the combustion of coal and wood is used for heating only in a small part of households both in Brno (0.46%) and in Šlapanice (1.36%), but the proportion of households in small villages near Šlapanice and Brno using wood or coal for heating is much higher (up to 11%). In addition, elements, such as Cd, As, Mn, Ni, Pb, V, Cu, Zn, Sb and so on, can also originate from emission of a municipal solid waste incinerator [65,66] located directly on the connecting line between the sampling site in Brno and Šlapanice, or from the incineration of waste in households (Sb, Cu, Pb, Sn, Ti, and Zn) [67]. Another possible source of elements, in particular Ca, Zn, Fe, Mn, Pb, Cd, As, Cu and Ni, is a cement plant [68][69][70] located north of Šlapanice and east of Brno. To verify this hypothesis, we calculated one-day backward trajectories for all the sampling campaigns in both the Brno and Šlapanice sampling sites ( Figure S1-S8). Detailed analysis revealed the transport of air masses from all directions, although during individual campaigns certain wind directions prevailed (Table 1). Moreover, the backward trajectory analysis shows the possibility of a long-range transport of air masses to Brno and Šlapanice from areas as far as several hundred kilometres.

Characterization of Water-Soluble Ions
We analysed 12 water-soluble ions in total (i.e., 5 inorganic cations, 6 inorganic anions and oxalate) in all submicrometre aerosol samples collected in Brno and Šlapanice. Average water-soluble ion concentrations in PM1 aerosols from both Brno and Šlapanice are summarized in Table 3 (campaigns of 2009) and Table 4 (campaigns of 2010). Ion species were important constituents of submicrometre aerosols both in Brno and Šlapanice, accounting for 27.5% and 26.8% of PM1 mass in Brno and Šlapanice, respectively. Ion species on average accounted for 27.2% of mass of PM1 aerosols. Ammonium, nitrate and sulphate were the three major ion species, contributing 32.3%, 36.2% and 23.4% of the total ion concentrations, respectively. SO4 2-, NO3 -and NH4 + are generally considered as secondary inorganic aerosol components (SIA; [71]). They derive from gas to particle conversion processes when SO2 is transformed to H2SO4 and nitrogen oxides to HNO3, followed by NH3 neutralization to form (NH4)2SO4 and NH4NO3. SO2 and NOX are products of combustion processes, whereas NH3 originates mainly from anthropogenic sources, such as agricultural activity, industry and traffic [36,72]. The sum of SO4 2-, NO3 -and NH4 + accounted together for 91.9% on average (86.1-97.5%) of total ion concentration and for 25. The majority of the analysed ions show significant seasonal variability with winter concentrations, significantly exceeding those in the summer (Tables 3 and 4). The SIA contribution decreased from a share of 32.4% in PM1 mass in winter to 17.8% in summer. The concentration of NO3 -was much higher in winter than in summer, which may be due to the low thermal stability of ammonium nitrate in the summer, favouring the conversion of particulate ammonium nitrate to gaseous nitric acid and ammonia [1,73]. In contrast, the difference between the winter and summer SO4 2-concentrations was much smaller compared to that of NO3 -, which indicates an active photochemical production of SO4 2-in summer [74].  Increased concentrations of F -, Cl -, NO2 -, SO4 2-, PO4 3-, NO3 -, NH4 + and K + in winter are mostly associated with burning of wood and coal in residential heating [11,18,19,36]. Combustion processes are the main sources of fluoride (i.e., coal), chloride (i.e., wood, coal, solid waste) and phosphate (i.e., coal, wood, traffic). K + , considered an inorganic tracer for biomass burning [14,62], formed the majority (i.e., 73.9%) of the total potassium in submicrometre aerosol. Enhanced winter concentrations of Cl -, F -, PO4 3-and K + were observed especially in Šlapanice, which is likely associated with local and regional combustion of wood and coal in residential heating [11,18,19]. Moreover, Cland Fmay also be emitted from the brickworks [75,76] located directly in Šlapanice. Cland PO4 3-serving as markers of plastic waste combustion may also originate from a large municipal waste incinerator located east of Brno or from the combustion of solid waste by households [67].
The concentrations of oxalate in winter were similar to those during summer campaigns, with the exception of Šlapanice in 2010. Oxalates, considered a major water-soluble organic compound in the aerosols, were reported as both a product of primary emissions from combustion processes (traffic, biomass burning, biogenic activity) and as a secondary product of atmospheric chemistry [77,78]. Oxalate correlated both with sulphate and K + (Table S1 and S2), which suggests that oxalate found in PM1 aerosols in Brno and Šlapanice originated from both biomass burning and secondary oxidation.
The concentration of ions in PM1 in Brno and Šlapanice were comparable with those found in Granada [34] and Katowice [47], while in other Polish cities, Zabrze and Racibórz, the concentrations of several ions (Cl -, K + , Na + , Ca 2+ , Mg 2+ ) were higher and others (NH4 + ) were lower [48]. The concentrations of ions at background station in Melpitz were comparable with the exception of a higher concentration of Mg 2+ and a lower concentration of Ca 2+ in comparison with results of this study [45]. A recent study from Prague reported generally higher concentrations of the analysed Cl -, SO4 2-, NO3 -and NH4 + in comparison with those from Brno and Šlapanice [20]. Recently, Cl -, NO2 -, SO4 2-and NO3 -have been analysed in Brno in PM2.5 particles, using the continuous aerosol sampler [16]. The concentrations of Cland SO4 2-found by the aerosol sampler were comparable with filter-based concentrations in this study, but the concentrations of NO2 -and NO3 -from the aerosol sampler were higher. This difference can be explained by the bimodal distribution of nitrate in PM1 and PM2.5 [6] and by sampling artefacts observed during the sampling of nitrite on filters [16].

Correlation Analysis
Pearson correlation coefficients provide another way to assess the sources of the analysed components in PM1 aerosols. The possible sources can be identified from the correlation matrix by analysing the value, which represents the linear coefficient of correlation between the species. The high correlation of two species suggests their identical sources. The results of correlation analysis of the studied elements and water-soluble ions and their relationship to the mass concentration of PM1 and temperature of air are shown in Table S1 (Šlapanice) and Table S2 (Brno), and are discussed in detail below.
Several elements (V, Cd, As, Fe, Zn, K) and most ions (NO3 -, Cl -, F -, PO4 3-, K + and NH4 + ) both in Šlapanice and Brno (p < 0.01) significantly positively correlated with PM1, but negatively correlated with temperature, which indicates a strong association of the mentioned elements and ions with combustion sources, especially in heating seasons.
Many elements and ions have the same overlapping sources. For example, Zn, Cu, Mn, Pb, As, Cd were found in the emissions from several combustion processes, such as coal combustion, biomass burning, an incinerator or vehicle emissions. The effect of biomass burning on the studied elements and ions was investigated through the correlation with K + , which is used as a marker of biomass (especially wood) burning [14,62]. K + correlated significantly (p < 0.01) in Brno with V, Cd, As, Sb, Cu, Ni, Mn, Fe, Zn, K, NO3 -, Cl -, F -, PO4 3-, NO2 -and NH4 + , and in Šlapanice with V, Cd, Sb, Zn, K, NO3 -, F -, NH4 + (p < 0.01) and with As and PO4 3-(p < 0.05), which suggests that wood burning was one of the major sources of the PM1 constituents during the campaigns in Brno and Šlapanice. This finding is in good agreement with the list of elements found in the emissions from the combustion of various biomass fuels [61,62]. Other elements (i.e., Zn, Cu, Cd and As) are important markers of coal combustion [60]. Significant cross-correlations among these elements (p < 0.01) confirm their common origin in combustion of coal. Moreover, some elements and ions originated from the municipal waste incinerator, such as Sb, Zn, Cd, As, Mn, Ni, Cl -, F -, and PO4 3-correlated significantly with other elements and ions in Brno more than in Šlapanice, which suggests a stronger effect of emissions from the municipal incinerator on the composition of PM1 in Brno than in Šlapanice. Ni did not correlate with any other element or ion in Šlapanice, but in Brno, Ni correlated significantly with V, Mn, Fe, K, NH4 + , PO4 3-(p < 0.01), Cd, NO3 -and Cl -(p < 0.05), which indicates the municipal incinerator as a predominant source of Ni in Brno. A similar difference in the correlations between the Brno and Šlapanice was also observed for Mn. The opposite situation was observed for Pb, which highly correlated with other PM1 components more in Šlapanice (i.e., Cd, Cu, Fe, Zn, K, SO4 2-, NO3 -, Cl -, F -, PO4 3-, NH4 + (p < 0.01), Mn, F -(p < 0.05) than in Brno (Ba, NO3 -(p < 0.01), V (p < 0.05)). No apparent correlation with other constituents in the PM1 samples collected in the Brno and Šlapanice sites was observed for Al and Ca, suggesting that Al and Ca have different sources compared to other elements and ions in PM1, such as soil dust resuspension for Al or emissions from a cement factory for Ca and Ca 2+ .

Source Apportionment
The correlation analysis groups the elements and ions on the basis of high pair correlations. The method provides only general information about the possible sources. Moreover, it does not evaluate the relationship between various sources, so overlapping of sources is problematic. To enhance the accuracy of emission source identification and their relative contribution, the method of factor analysis was also applied. The factor analysis belongs to multivariate statistical methods and is commonly used to identify the sources of PM aerosols [29,31,32,79].
In order to identify the aerosol emission sources, the factors were extracted by using the principal component analysis method (PCA) [41] and further rotated by using the varimax criterion, which achieves a simpler structure of the orthogonal factor model and also better interpretable factors. The optimal number of factors was estimated based on the Kaiser criterion and scree plot [41]. All computations were performed by using the software R version 3.6.1.
The principal component analysis was performed on the concentrations of elements and ions in PM1 aerosols collected in Šlapanice and Brno sites. The six extracted factor loadings from PCA analysis in Šlapanice and Brno dataset are given in Tables 6 and 7, and Figures 6 and 7, respectively. Only factor loadings above 0.1 are shown and values greater than 0.5 are in bold. The last lines in the tables (Variance and Cumul. Var., %) show the proportions of the total data variance explained by each individual factor and total explained variance, respectively.
The first factor was identified as a combustion source. It is dominated by K + , K, NH4 + , V and Cd, followed by NO3 -, oxalate and Zn. Other identified tracers are As, Sb, Pb, PO4 3-and SO4 2-. These can be related to biomass (mainly wood) burning and coal combustion in residential heating [11,18,19,35] or in a power plant burning coal and biomass located 50 km southeast of Šlapanice. The factor may also include local sources, such as engine vehicle emissions (Zn, Pb, Mn) [80] or emissions from brickworks built on a small hill at the north side of Šlapanice.
The second factor was identified as coal combustion. The factor was associated with high loadings of Fe, Mn, Pb and Na + and moderate loadings of Cd, As, Cu, SO4 2-, Cl -, NO3 -and NH4 + . Coal is used for household heating to some degree in Šlapanice and, in particular, in villages around Šlapanice. Moreover, a coal-fired power plant is located in a distance of 50 km southeast of Šlapanice.
The third factor with high loadings of Ca and Ca 2+ may be explained mainly by emissions from a cement factory [68][69][70]81]. The cement factory in Mokrá is situated in a short distance, about 6.5 km northeast of Šlapanice. The fourth factor with high loadings of Cl -, F -, PO4 3-and NO3 -and moderate loadings of Cd, Cu, Zn, Pb, K and Na + was identified as emissions from a municipal waste incinerator [65,66]. The large municipal waste incinerator, located at the southeast side of Brno and at a distance of 5 km northwest of Šlapanice, is used for the entire agglomeration of Brno. It is also necessary to take into account a partial contribution from waste incineration in households in Šlapanice and surrounding villages [67].
The fifth factor identified as vehicle exhausts was characterized with high loadings of Sb and Cu and moderate loadings of Cd, Zn, K, Pb and K + . Cu, Sb and Zn are used as markers of vehicle-related sources [57]. They generally originate from abrasions of tires or break wear. However, Cu, Zn and Sb present in PM1 aerosols, originating probably from the exhaust emissions from diesel and petrol engines as Cu, Zn and Sb, along with other elements, are present directly in the fuel or lubricating oils [36,57,80,82].
The sixth factor with moderate loadings of Ni, Ba, Pb, SO4 2-and NO3 -, represented emissions from industry. Several small production workshops dealing with the processing of metals are located on the outskirts of Šlapanice or nearby.
The first factor was identified as the municipal waste incinerator [65,66] situated about 6 km southeast of the sampling site. The factor is associated with high loadings of Mn, Fe, Cl -, Ni, Zn, V, K, NO3 -and NH4 + and moderate loadings of K + , As, Sb, Cu, Cd, Pb and F -.
The second factor was identified as vehicle exhausts. It is dominated by Sb and Cu and by moderate loadings of Zn, K, Cd, F -, K + , Mn, As, Ba, Cland Fe [35,36,57]. The intensity of traffic, including cars and trams on the Veveří street, is quite high, with frequent queues of standing cars in front of the sampling site due to a nearby intersection.
The third factor representing the secondary aerosol production was identified as a secondary sulphate with SO4 2-and NH4 + as the two main markers. They derive from gas to particle conversion processes from SO2 oxidation and NH3 neutralization [30,57]. SO2 originates largely from coal and biomass combustion, while NH3 results from agricultural and traffic emissions [36,72]. This factor represents regional or long-range transport of aerosols.
The fourth factor was dominated by high loadings of Ca 2+ , Ca and Na + , followed by moderate loadings of oxalate and SO4 2-, and represented the emissions from a cement factory [68][69][70] situated about 14 km east of the Brno sampling site.
The fifth factor was identified as an industry source. It was characterized with high loadings of Pb and Ba and moderate loading of NO3 - [57]. The factor is associated with emissions from various industrial plants around Brno.
The sixth factor with PO4 3-, K + , Cl -, F -, oxalate, NH4 + and Zn, as the main indicators, was related to biomass burning [30,35]. The factor relates mainly to the regional transport of emissions from biomass (mainly wood) burning within residential heating in villages near Brno.

Comparison of Results for Šlapanice and Brno
The major sources of PM1 components identified by PCA in both localities are quite similar, but they differ with their contribution. PCA showed that local pollution from the municipal waste incinerator (accounting for 23%) and vehicle exhausts (18%) prevailed in Brno, while Šlapanice was more burdened by the transport of pollution from regional sources, such as the municipal waste incinerator in Brno (12%), the cement factory in Mokrá (13%) or coal and wood combustion (14 and 23%) in the nearby villages. This hypothesis was confirmed by one-day backward trajectories (Figures S2, S4, S6 and S8) that indicate possible regional transport of air masses to Šlapanice from all directions, although the transport from east is less frequent. Villages, predominantly south or north of Šlapanice could thus be an important source of air pollution originated from biomass and coal combustion in the frame of residential heating. Moreover, there is a coal and biomass burning power plant located southeast of Šlapanice within a distance of 1 day of transport. Trajectories also confirmed the transport of air masses to Šlapanice from a northeast or northwest direction, where the cement factory and municipal incinerator are located. One-day backward trajectories calculated for Brno location (Figures S1, S3, S5 and S7) also indicated possible regional transport of air masses to Brno from all directions, which confirms findings of the factor analysis which identified the cement factory and biomass burning (i.e., within residential heating in surrounding villages or in a biomass-burning power plant located southeast of Brno) as regional sources of air pollution in Brno. Finally, we cannot neglect the contribution of long-range transport of polluted air from remote areas at distances of as far as several hundred kilometres to the pollution of both locations.
It should be borne in mind that, due to the limited sampling period at both locations, the findings concerning the source apportionment are valid only for the sampling period, while for the rest of the year, the PM1 sources or their contribution may vary.

Conclusions
The concentrations of elements and water-soluble ions in PM1 aerosols in Brno and Šlapanice, representing a large city and a small town in the Czech Republic, were compared in heating and non-heating seasons in 2009 and 2010. The average PM1 mass concentration was 19.1 and 9.65 μg m −3 in winter and summer in Brno, respectively, and 30.8 and 10.0 μg m −3 in winter and summer in Šlapanice, respectively. In the winter season, ions formed a significant part of PM1 mass of aerosols, accounting for 34.6% (26.8-41.6%), while in the summer season the contribution of ions to PM1 mass decreased to 19.7% (14.8-24.9%). Ammonium, nitrate and sulphate, the three major ion species, accounted for 91.9% of a total ion concentration and 25.1% of PM1 mass. The contribution of elements to PM1 mass was much smaller: 2.17% (1.70-2.76%) in winter and 1.37% (0.92-1.97%) in summer. A more recent seasonal studies in Brno and Šlapanice [11,18,19] ascribed the rest of PM1 mass to other components: organic material and elemental carbon that accounted for 38.8% and 7.00% of PM1 mass in winter, respectively, and for 36.5% and 8.62% of PM1 mass in summer, respectively.
The differences in the concentrations of PM1 aerosols, elements and water-soluble ions in the corresponding seasons (especially winter) in 2009 and 2010 may have mainly been due to different meteorological and dispersal conditions in those two years. The variations in the concentrations of PM1, elements and water-soluble ions in Brno and Šlapanice during the same season result from a change in the actual emission and meteorological and dispersal situation at the relevant site.
The backward trajectory analysis confirmed that the concentrations of elements and ions in aerosols collected in Brno or Šlapanice do not depend only on local emission sources but are affected significantly also by regional transport of polluted air from various sources both nearby (e.g., surrounding villages, cement factory etc.) and by a long-range transport of polluted air from sources at larger distances from both the studied locations, such as a power plant southeast of Brno and Šlapanice, or even from more distant areas, such as a heavily polluted region in Ostrava or southern Poland situated north of Brno and Šlapanice.
The source apportionment of the PM1 samples collected in Brno and Šlapanice was performed using PCA. The six major sources of PM1 components identified by PCA in both localities are quite similar in composition, although differing in their fractional contribution. Coal and biomass (largely wood) combustion, a municipal waste incinerator, vehicle exhausts, a cement factory and industry were identified as major sources at both localities. Both sampling sites were burdened by both local and regional pollution. The municipal waste incinerator (23%) and vehicle exhausts (17%) identified as the two major sources of PM1 in Brno indicate a predominant effect of local sources in Brno, while Šlapanice was more burdened by the transport of pollution from regional sources, such as the municipal waste incinerator in Brno (12%) or the combustion of wood and coal in the nearby villages (37%). The transport of aerosols from sources at larger distances from Šlapanice (for example a coaland biomass-fired power plant near the border with Slovakia) or a long-range transport of PM1 from neighbouring countries should also be taken into account. However, the short sampling period at both locations restricts the validity of conclusions concerning the sources of PM1 aerosols at both locations only to the sampling period, while in rest of year, the PM1 sources or their contribution may vary.  Table S1: Correlation analysis between PM1, temperature, elements and ions in Šlapanice; Table S2: Correlation analysis between PM1, temperature, elements and ions in Brno.