3. 1 Data Analysis
To select the input variables for PMF modeling, the signal-to-noise (S/N) ratios were calculated (Paatero and Hopke, 2003). Variables with very low S/N ratio (≤0.2) were excluded from the PMF analysis, whereas variables with S/N ratio between 0.2 and 2 were downweighted by a factor of 2 or 3. In this study, 21 species (OC, EC, SO4 2-, NO3-, NH4 +, Al, Br, Ca, Cl, Cu, Fe, K+, K, Mg, Mn, Na, Ni, Pb, Si, Ti, and Zn) were selected for PMF modelling, and the weak variable (Cu) was downweighted. In the case of two bad variables (Br and Ni), these variables were exclude from the PMF analysis in principle due to bad effect for the PMF result. However, these variables were included with adjusted error values because of two variables are important marker species for identification of source. Table 1 shows arithmetic mean, standard deviation, geometric mean, minimum and maximum value, number of BDL, and S/N ratio for individual species during the sampling periods at sampling site.
Table 1.
Summary statistics for the PM10 (μg/m3) and species concentrations (ng/m3) at the sampling site.
|
A.M |
S.D |
G.M |
Min. |
Max. |
No. of BDL(%) |
S/N Ratio |
PM10 |
62.75 |
47.39 |
49.42 |
5.60 |
358.60 |
- |
- |
OC |
9582.92 |
5344.04 |
8173.11 |
700.00 |
28050.00 |
0 (0.0) |
- |
EC |
2830.74 |
1612.57 |
2361.55 |
461.77 |
10700.00 |
0 (0.0) |
- |
SO4 2- |
8518.38 |
8043.36 |
5639.76 |
461.17 |
49663.50 |
0 (0.0) |
- |
NO3 - |
7577.40 |
6829.82 |
5154.49 |
509.32 |
48878.06 |
0 (0.0) |
- |
NH4 + |
5830.95 |
4810.58 |
4161.21 |
621.21 |
29129.71 |
0 (0.0) |
- |
Al |
854.55 |
1313.64 |
533.72 |
55.58 |
11542.83 |
10 (6.1) |
160.0 |
Br |
170.38 |
143.15 |
116.83 |
44.37 |
294.76 |
160 (97.6) |
0.1 |
Ca |
699.82 |
742.24 |
515.76 |
98.21 |
6108.98 |
1 (0.6) |
2617.1 |
Cl |
712.20 |
804.70 |
379.56 |
39.63 |
3360.38 |
41 (25.0) |
61.5 |
Cu |
39.29 |
19.52 |
34.65 |
8.71 |
98.29 |
93 (56.7) |
1.2 |
Fe |
1075.24 |
1208.82 |
809.01 |
165.02 |
10935.01 |
0 (0.0) |
- |
K+ |
977.79 |
707.00 |
635.34 |
2.37 |
4231.81 |
0 (0.0) |
- |
K |
751.47 |
689.51 |
555.12 |
68.14 |
5331.03 |
1 (0.6) |
2302.8 |
Mg |
270.12 |
311.19 |
205.97 |
76.70 |
2153.77 |
49 (29.9) |
21.4 |
Mn |
50.09 |
41.29 |
41.19 |
11.80 |
315.66 |
45 (27.4) |
13.0 |
Na |
413.01 |
278.16 |
347.51 |
93.92 |
1697.47 |
0 (0.0) |
- |
Ni |
9.14 |
3.33 |
8.77 |
7.21 |
12.98 |
155 (94.5) |
0.1 |
Pb |
130.60 |
47.88 |
121.99 |
53.96 |
238.07 |
6 (3.7) |
143.0 |
Si |
2275.93 |
3663.28 |
1395.91 |
110.43 |
32283.23 |
0 (0.0) |
- |
Ti |
123.56 |
129.41 |
98.55 |
32.22 |
918.04 |
62 (37.8) |
9.5 |
Zn |
147.85 |
97.73 |
115.70 |
2.23 |
436.96 |
16 (9.8) |
36.1 |
3. 2 Determination of the Number of Factors
In order to estimate the optimal number of factors, the mathematical PMF diagnostics was explored. The PMF diagnostics (model error, Q, rotational ambiguity, rotmat, etc.) were based on Lee et al., 1999. In this study, the scaled residual matrix, Q value, and the rotmat matrix (indicates the rotational uncertainty, p×p matrix of standard deviations of rotational coefficients, where p is the number of sources) were used to determine the number of factors. The optimal number of sources was determined to be nine factors based on examination of the scaled residuals and the Q value. In order to explore the rotational freedom, PMF used a parameter, FPEAK, to control the addition and subtraction of the factors. The FPEAK forces PMF to add one g vector to another and to subtract the corresponding f factors from each other, thereby yielding more physically realistic solutions. The Q values were plotted against the FPEAK value to explore the rotational space where only small changes in the Q values are observed. It is also possible to use pairwise scatter plots of the g vectors to help define the FPEAK value (Hwang and Hopke, 2007). In this study, FPEAK values between -1.0 and 1.0 were examined and the value of FPEAK=0.0 provided the most physically meaningful solution.
Also, the parameter FKEY was used to obtain the reasonable source profiles. If specific species in the source profiles do not seem to be realistic based on comparison with measured source profiles and prior analysis of similar data, it is possible to pull values toward zero to obtain the reasonable source profile using the FKEY matrix. The FKEY matrix of integer values has the same dimension as F matrix, where F matrix is a p×m source profile matrix. The details are reported in previous studies (Hwang and Hopke, 2011; Zhao and Hopke, 2004). In this study, values of all elements in the FKEY matrix were set to zero, except for a value of 1 for OC in field burning source.
3. 3 Source Identification and Apportionment
To estimate source apportionment and source profiles in actual units, scaling coefficients were obtained using a multiple linear regression against the measured PM10 mass. Fig. 2 presents the source profiles (value±standard deviation) obtained for the nine factors from the PMF solution at the sampling site in Seoul, where a thin bar denotes explained variation value. Fig. 3 shows the temporal variations of contributions from each source. Table 2 provides a comparison of seasonal contributions for each source and Fig. 4 shows the average source contributions for the whole sampling period. Also, the average source contributions for weekdays and weekend days are presented in Fig. 5.
Fig. 2
Source profiles of the resolved sources measured at the sampling site.
Fig. 3
Temporal variation of source contributions for the sampling site constructed using the PMF model.
Table 2.
Average seasonal source contributions for the sampling site during sampling period.
|
Winter |
Spring |
Summer |
Fall |
AVG. |
μg/m3 |
% |
μg/m3 |
% |
μg/m3 |
% |
μg/m3 |
% |
μg/m3 |
% |
Secondary Nitrate |
8.93 |
12.3 |
7.78 |
10.7 |
2.38 |
6.3 |
3.44 |
7.1 |
5.35 |
9.3 |
Motor Vehicles |
12.63 |
17.3 |
9.81 |
13.5 |
8.21 |
21.7 |
8.64 |
17.7 |
9.53 |
16.6 |
Road Salt |
9.33 |
12.8 |
3.94 |
5.4 |
0.14 |
0.4 |
1.90 |
3.9 |
3.33 |
5.8 |
Industry |
2.81 |
3.9 |
2.48 |
3.4 |
2.01 |
5.3 |
3.44 |
7.1 |
2.80 |
4.9 |
Airborne Soil |
10.41 |
14.3 |
19.58 |
27.0 |
1.91 |
5.0 |
6.15 |
12.6 |
9.87 |
17.2 |
Aged Sea Salt |
3.77 |
5.2 |
5.40 |
7.4 |
1.24 |
3.3 |
3.29 |
6.8 |
3.58 |
6.2 |
Field Burning |
3.44 |
4.7 |
1.86 |
2.6 |
2.35 |
6.2 |
5.27 |
10.8 |
3.45 |
6.0 |
Secondary Sulfate |
9.13 |
12.5 |
12.28 |
16.9 |
12.46 |
32.9 |
5.32 |
10.9 |
9.27 |
16.2 |
Road Dust |
12.37 |
17.0 |
9.46 |
13.0 |
7.19 |
19.0 |
11.27 |
23.1 |
10.15 |
17.7 |
Sum |
72.82 |
100.0 |
72.59 |
100.0 |
37.89 |
100.0 |
48.71 |
100.0 |
57.33 |
100.0 |
Fig. 4
Comparison of the average seasonal contributions for each source.
Fig. 5
The average source contributions for weekdays and weekend days at the sampling site.
The species contributing to the first source included NO3 -, SO4 2-, and OC. This profile was identified as secondary nitrate and contributed 9.3% (5.35 μg/m3) to the total PM10 mass concentration. In general, secondary nitrate is known to be seasonal with high contributions in winter season because lower temperatures and high humidity help the formation of secondary nitrate particles (Seinfeld and Pandis, 1998). The seasonal average mass contributions of secondary nitrate show a peak in winter (winter 8.93 μg/m3>spring 7.78 μg/m3>fall 3.44 μg/m3>summer 2.38 μg/m3). To test if the seasonal variation of secondary nitrate is statistically significant, ANOVA was used. There were significant seasonal variations (P value<0.0001) on the contribution of secondary nitrate with spring and winter peaks. The secondary nitrate showed no significant differences in the mean contributions between the weekdays (5.35 μg/m3) and the weekends (5.37 μg/m3) (Fig. 5). The NPR plot for the secondary nitrate source is presented in Fig. 6. This plot indicated secondary nitrate sources located to the west, south, southwest, and north-northeast as shown in Fig. 6.
Fig. 6
Nonparametric regression (NPR) results based on source contributions obtained by the PMF analysis of data from the sampling site.
The major marker species contributing to the second source profile included OC, EC, Ca, Fe, K, Zn, and Pb, and this profile was classified as motor vehicle source. The two types of motor vehicle sources (i.e. gasoline and diesel emissions) were not separated in this study, though it may be assumed that the contribution of gasoline vehicles is higher than the contribution of diesel vehicles, as land use around the sampling site was mixed commercial and residential, with gasoline engines being far more common than diesel engines. The peak seasonal mass contribution of motor vehicle source was the summer (21.7%, 8.21 μg/m3), the fall (17.7%, 8.64 μg/m3) and winter (17.3%, 12.63 μg/m3) contributions were higher than the spring contribution (13.5%, 9.81 μg/m3). This may be the result of the lower mixedlayer heights in the winter relative to the rest of the year, and thus, a reduced dispersion of these ground level emissions (Hwang and Hopke, 2006). However, in the case of high summer contributions, it will be needed to obtain an understanding of the reason. As shown in Fig. 5, motor vehicle contributions showed higher weekday values (10.35 μg/m3) than those on the weekends (6.78 μg/m3). This result shows that motor vehicle was mainly from vehicles primarily operating in weekdays. The NPR plot indicated motor vehicle sources located to the northeast, southeast, west, and south as shown in Fig. 6 and this results showed good agreement with Heo et al. (2009).
The third source was classified as a road salt source with high abundances of Cl, OC, EC, NO3 -, Si, Al, and Na. The peak seasonal mass contribution of road salt was in the winter (12.8%, 9.33 μg/m3) because road salt released as a deicer on roads in snowfall season. The weekday contribution of road salt was higher than the weekend contribution, similar to the increased contribution of motor vehicles to the weekday samples. Also, the NPR results of road salt showed similar results to the motor vehicle source.
The fourth source was determined to be industry related source. The major species contributing to this source included NO3 -, NH4 +, Na, Si, Ca, Fe, and Cu. This source contributed 4.9% (2.80 μg/m3) to the total PM10 mass concentration. In the case of this source, the fall contribution (7.1%, 3.44 μg/m3) was higher than other seasonal contributions (winter 3.9%, 2.81 μg/m3; spring 3.4%, 2.48 μg/m3; summer 5.3 %, 2.01vμg/m3). This sampling site is located about 35 km west of a number of industrial complexes (Incheon area) and about 34 km southwest of the Sihwa and Banwol national industrial complexes, which include metal processing, non-ferrous metal smelting, and petroleum chemical process facilities (Gyeonggi-do, 2012). The industry source NPR plot showed good agreement with local source direction, and therefore local industrial activities are likely influence to the industry source of PM10 in sampling area.
The fifth source was identified as airborne soil with high contributions of Si, Al, Fe, K, Ca, Mg, and Ti. The temporal variation of the source contribution plot showed very strong contributions at the end of March, particularly March 27 and 31, 2007 (Fig. 3). It is suggests that the air mass was transported from the Eastern Mongolia and the Gobi Desert and influence of long range transport of Asian dust. Actually, Asian dust episodes have been observed on March 27 and 31, 2007 (KMA, 2007). The peak seasonal mass contribution of airborne soil source was in the spring (27.0%, 19.58 μg/m3).
The sixth source profile was assigned to aged sea salt. The Na and SO4 2-were major species contributing to the aged sea salt source along with minor species such as K, Ca, Mg, and Fe. Although the main species in sea salt are known to be Na, Cl, SO4 2-, K, and Ca (Hopke, 1985), only Na showed a high contribution in association with SO4 2-. The Cl was depleted because NaCl was converted into Na2SO4 and NaNO3as a result of reactions of NaCl with gaseous H2SO4 and gaseous HNO3, respectively (Hwang and Hopke, 2006; Seinfeld and Pandis, 1998). The aged sea salt source peaked in the spring (spring 7.4%, 5.40 μg/m3>fall 6.8%, 3.29 μg/m3>winter 5.2%, 3.77 μg/m3> summer 3.3%, 1.24 μg/m3). This result suggests that the sea salt particles were transported from the Ocean (Yellow Sea) to the west of the site in spring time due to the dominance of wind from the west in the spring time, which is the typical wind direction.
The species associated with the seventh source included OC, EC, SO4 2-, NH4 +, K+, K, Si, Ca, and Fe, and this profile was classified as field burning (included wood combustion and biomass burning) and contributes 6.0% (3.45 μg/m3) to the total PM10 mass concentration. The peak seasonal mass contribution of field burning was the fall (10.8%, 5.27 μg/m3). This source maybe caused by residential wood burning and illegal field burning, such as field burning after harvesting, biofuel burning for heating and cooking, and forest fires that occurred outside of the Seoul area (Heo et al., 2009). There was no significant difference between the concentrations measured on weekdays versus those measured on weekends (Fig. 5).
The eighth source was classified as a secondary sulfate with high abundance of SO4 2-, NH4 +, and NO3 -, and contributed 16.2% (9.27 μg/m3) to the total PM10 mass concentration. The average seasonal mass contributions show its peak contribution to be in summer (32.9%, 12.46 μg/m3) (Table 2) because formation of secondary sulfate was enhanced by the increased photochemical reaction during the summer time with strong solar intensity. The SO2 emitted during various combustion processes was converted into H2SO4 to form SO4 2- as a result of its photochemical reactions (US EPA, 1999). The NH4 + had relatively strong correlations between SO4 2- and NO3 - (0.8 and 0.7), these would presumably exist in the ambient air as secondary aerosols in a sulfate form such as (NH4)2SO4 and (NH4)3H(SO4)2 as a consequence of homogeneous and heterogeneous reactions (Finlayson-Pitts and Pitts, 2000). As shown in Fig. 6, the NPR plot indicated secondary sulfate source located to the southwest and south. In this area, the largest sources of SO2 were ships and a coal-fired power plant. In the case of coal-fired power plant, with a production capacity of 4,000 MW, emitted 4485 ton, 5616 ton, 6040 ton of SO2 during 2005-2007 (KEWP, 2012). Also, the near shore Yellow Sea area where there are active shipping lanes. Many ships burning high sulfur residual oil such that this area can contribute to the secondary sulfate source.
The final source was identified as road dust source. This source profile was similar to the airborne soil source profile, however, these source profiles showed some differences for Zn and Pb. The seasonal average mass contributions of road dust shows a peak in fall (23.1%, 11.27 μg/m3). The seasonal variation of road dust showed similar to the motor vehicle source. This result suggests that road dust particles due to resuspended in the ambient air from operation of motor vehicle.
Fig. 7 presents a comparison of the predicted PM10 contributions from all of the identified sources with measured PM10 concentrations. The PMF resolved sources effectively reproduced the measured values (R2=0.93) and account for most of the variation in the PM10 concentrations (slope=0.79).
Fig. 7
Comparison of the predicted total PM10 mass concentrations from the PMF analysis with measured PM10 mass concentrations for the sampling site.
3. 4 Comparison of the Source Apportionment Studies for PM10
In this study, the average contribution of each source was estimated using the PMF model for PM10 in the Seoul, Korea. The source apportionment studies for PM10 were performed beginning of the 1990 in the Korea, therefore reported source apportionment studies for PM10 have been very limited. Several previous studies for the metropolitan area (i.e. Seoul, Gyeonggido, Incheon-si) have identified a total of 10 sources such as fuel combustion (included oil and coal combustion), soil, industry, Fe-related, vehicle (included gasoline and diesel), cement, incinerator, secondary aerosol (included secondary sulfate and secondary nitrate), field burning, and sea salt source, respectively.
Table 3 shows the comparison of PM10 source contributions for the metropolitan area in Korea using the various receptor models. The comparison of source contributions used only identified same sources between this study and previous studies. In the case of soil source contribution for this study, this contribution included contribution of road dust source. The contribution of the soil source was within the range of 16.8 %-35.9%, these results show differences of contribution presumably due to various sampling site, existence of construction works, and inclusion of road dust. The contribution of the industry related source was within range of 3.9%-8.6%, as quite a low contribution to the total PM10 mass concentration, and field burning source contribution showed within range of 7.2%-18.3 %. The contribution of motor vehicle source estimated within range of 7.9%-18.0%. As mentioned above, the KMOE established the “Special Act on Metropolitan Air Quality Improvement” in 2005 in order to improve the PM10 annual concentration in Seoul metropolitan area. It included reduction strategies and management system for vehicle emissions thus vehicle source contributions will be expected lower. Lastly, the contribution of secondary aerosol source estimated within range of 5.1%-17.0%. In particular, in the case of secondary aerosol, the study of the intensive source apportionment for emitted from metropolitan area in Korea and long range transport from China will be the subject of a future study.
Table 3.
Comparison of source contributions (%) for the metropolitan area in Korea using the various receptor models.
asampling site (Seoul), TTFA model
bsampling site (Suwon), CMB model
csampling site (Suwon), TTFA model
dsampling site (Suwon), PMF model
esampling site (Seoul and Incheon), PMF model