Descriptive statistics of Death, Meteorological data and PM2.5 composition
Tables 1 and 2 shows the descriptive statistics of explanatory variables for the four cities. Due to the largest population in SL, daily death counts were much larger in SL while those in the other cities were quite similar to each other. The proportions of CVD and RD death in total death in SL and DJ were comparable to those of the whole country (21.4% and 10.8%) whereas the CVD and RD death in GJ (26.3%, 15.8%) and US (30.8% and 15.4%) accounted for higher proportions of total death.
Table 1
Summary of daily mortality data and meteorological conditions in SL, DJ, GJ and US during 2014-2018.
SL
|
Avg. (S.D.)
|
Min.
|
P25
|
Median
|
P75
|
Max.
|
Unit
|
DJ
|
Avg. (S.D.)
|
Min.
|
P25
|
Median
|
P75
|
Max.
|
NA
|
108 (14)
|
69
|
98
|
107
|
115
|
182
|
counts
|
NA
|
17 (5)
|
5
|
14
|
17
|
20
|
33
|
CVD
|
24 (6)
|
8
|
20
|
23
|
27
|
44
|
CVD
|
4 (2)
|
0
|
2
|
3
|
5
|
11
|
RD
|
11 (4)
|
2
|
8
|
10
|
13
|
29
|
RD
|
2 (2)
|
0
|
1
|
2
|
3
|
8
|
Temp
|
13.3 (10.8)
|
-14.8
|
3.8
|
14.8
|
22.9
|
33.7
|
℃
|
Temp
|
13.7 (10.2)
|
-12.7
|
4.7
|
14.9
|
22.5
|
33.4
|
RH
|
59.4 (14.8)
|
21.8
|
48.4
|
59.5
|
69.5
|
99.8
|
%
|
RH
|
69.8 (13.8)
|
27.5
|
60.0
|
70.9
|
79.3
|
99.3
|
BP
|
1006.2 (7.9)
|
981
|
999.8
|
1006.5
|
1012.6
|
1026.8
|
hPa
|
BP
|
1008.6 (7.9)
|
985.7
|
1002.3
|
1008.9
|
1015.0
|
1029.1
|
GJ
|
|
|
|
|
|
|
|
US
|
|
|
|
|
|
|
NA
|
19 (5)
|
6
|
16
|
18
|
22
|
36
|
counts
|
NA
|
13 (4)
|
2
|
9.25
|
12
|
14
|
28
|
CVD
|
5 (3)
|
0
|
3
|
4
|
6
|
13
|
CVD
|
4 (2)
|
0
|
2
|
3
|
4
|
13
|
RD
|
3 (2)
|
0
|
1
|
2
|
3
|
11
|
RD
|
2 (2)
|
0
|
0
|
1
|
2
|
8
|
Temp
|
14.7 (9.4)
|
-9.5
|
6.3
|
15.75
|
22.7
|
32.0
|
℃
|
Temp
|
14.7 (8.7)
|
-8.2
|
7.2
|
15.7
|
21.8
|
31.8
|
RH
|
69.0 (15.2)
|
23.9
|
58.2
|
69.8
|
79.5
|
99.0
|
%
|
RH
|
65.9 (19.4)
|
16.8
|
50.9
|
68.7
|
81.3
|
99.8
|
BP
|
1008.2 (7.8)
|
985.3
|
1001.9
|
1008.4
|
1014.4
|
1027.6
|
hPa
|
BP
|
1008.7 (8.1)
|
982.5
|
1002.6
|
1008.6
|
1014.6
|
1030
|
Table 2
Summary of PM2.5 mass and chemical constituents concentration of study cities during 2014-2018.
SL
|
Avg. (S.D.)
|
Min.
|
P25
|
Median
|
P75
|
Max.
|
unit
|
DJ
|
Avg. (S.D.)
|
Min.
|
P25
|
Median
|
P75
|
Max.
|
PM2.5
|
27.7 (17.2)
|
1.0
|
16.5
|
23.6
|
33.8
|
153.4
|
㎍/㎥
|
PM2.5
|
29.2 (16)
|
1.71
|
18.1
|
26.5
|
36.7
|
117.2
|
SO42−
|
5.4 (4.7)
|
0.0
|
2.3
|
4.1
|
6.8
|
39.1
|
SO42−
|
5.2 (3.7)
|
0.2
|
2.5
|
4.2
|
6.9
|
28.3
|
NO3−
|
6.1 (5.9)
|
0.0
|
1.7
|
4.3
|
8.4
|
35.6
|
NO3−
|
5.4 (5.2)
|
0.1
|
1.5
|
3.6
|
7.8
|
31.2
|
Cl−
|
0.3 (0.4)
|
0.0
|
0.1
|
0.2
|
0.5
|
2.6
|
Cl−
|
0.4 (0.4)
|
0.0
|
0.1
|
0.2
|
0.5
|
2.8
|
Na+
|
0.1 (0.1)
|
0.0
|
-
|
0.1
|
0.2
|
1.0
|
Na+
|
0.2 (0.3)
|
0.0
|
0.1
|
0.2
|
0.3
|
2.5
|
NH4+
|
3.9 (3.2)
|
0.0
|
1.7
|
3.0
|
5.1
|
24.0
|
NH4+
|
3.6 (2.6)
|
0.1
|
1.8
|
3.1
|
4.8
|
18.0
|
K+
|
0.2 (0.2)
|
0.0
|
-
|
0.1
|
0.2
|
1.0
|
K+
|
0.2 (0.2)
|
0.0
|
0.1
|
0.1
|
0.2
|
1.5
|
OC
|
3.8 (1.9)
|
0.3
|
2.5
|
3.4
|
4.7
|
16.8
|
OC
|
4.5 (2.1)
|
0.6
|
3.0
|
4.2
|
5.7
|
16.5
|
EC
|
1.3 (0.6)
|
0.1
|
0.8
|
1.1
|
1.6
|
4.5
|
EC
|
1.4 (0.7)
|
0.3
|
0.8
|
1.2
|
1.7
|
4.9
|
S
|
3.1 (2.3)
|
0.1
|
1.5
|
2.4
|
3.9
|
23.3
|
S
|
3.9 (3.9)
|
0.1
|
1.5
|
2.6
|
4.8
|
40.1
|
Ca
|
65.8 (69.6)
|
1.5
|
29.4
|
48.5
|
80.9
|
1105.0
|
ng/㎥
|
Ca
|
58 (73.8)
|
1.5
|
16.5
|
36.4
|
75.5
|
1113
|
Ti
|
9.8 (7.8)
|
0.5
|
5.4
|
8.3
|
12.1
|
106.5
|
Ti
|
7.7 (8)
|
0.5
|
3.3
|
5.8
|
9.5
|
78.6
|
V
|
3.7 (4.6)
|
0.3
|
0.3
|
1.8
|
5.0
|
31.6
|
V
|
3.0 (2.6)
|
0.3
|
1.1
|
2.3
|
4.2
|
16.9
|
Cr
|
1.3 (1)
|
0.3
|
0.3
|
1.1
|
1.7
|
7.5
|
Cr
|
1.2 (1.2)
|
0.3
|
0.3
|
0.9
|
1.6
|
10.1
|
Mn
|
12.3 (8.5)
|
0.3
|
6.5
|
10.1
|
15.7
|
56.1
|
Mn
|
10 (6.6)
|
0.3
|
5.2
|
8.6
|
13.4
|
39.9
|
Fe
|
195 (112.2)
|
20.9
|
121
|
172.8
|
242.9
|
1074.5
|
Fe
|
183 (113.9)
|
16.6
|
103.7
|
160
|
234.3
|
1002.7
|
Ni
|
1.6 (1.6)
|
0.2
|
0.4
|
1
|
2.2
|
10.1
|
Ni
|
1.6 (1.4)
|
0.2
|
0.7
|
1.3
|
2.0
|
16.8
|
Cu
|
7.5 (5.2)
|
1.0
|
3.8
|
6.6
|
9.8
|
35.0
|
Cu
|
7 (4.4)
|
1.0
|
3.8
|
6.1
|
9.3
|
33.6
|
Zn
|
65.5 (43.6)
|
2.5
|
35.4
|
54
|
82.7
|
302.3
|
Zn
|
50.7 (32)
|
2.2
|
27.2
|
43.9
|
66
|
254.0
|
As
|
4.0 (3.3)
|
0.3
|
1.7
|
3.4
|
5.6
|
25.7
|
As
|
2.7 (2.6)
|
0.3
|
0.6
|
1.9
|
4.0
|
20.0
|
Se
|
1.2 (1.2)
|
0.3
|
0.3
|
0.9
|
1.6
|
8.6
|
Se
|
1.8 (1.4)
|
0.3
|
0.8
|
1.5
|
2.5
|
8.5
|
Br
|
8.9 (7)
|
0.1
|
4.2
|
6.9
|
11.0
|
63.1
|
Br
|
7.3 (5.8)
|
0.1
|
3.2
|
5.9
|
9.8
|
48.6
|
Pb
|
21.8 (15.7)
|
0.7
|
11.9
|
18.2
|
27.7
|
132.9
|
Pb
|
18.8 (13.9)
|
0.7
|
8.7
|
15.3
|
24.7
|
115.3
|
Table 2
GJ
|
Avg. (S.D.)
|
Min.
|
P25
|
Median
|
P75
|
Max.
|
unit
|
US
|
Avg. (S.D.)
|
Min.
|
P25
|
Median
|
P75
|
Max.
|
PM2.5
|
26.9 (15.8)
|
3.3
|
16.2
|
23.5
|
33.5
|
132.9
|
㎍/㎥
|
PM2.5
|
20.4 (11.8)
|
2
|
11.7
|
17.8
|
26.1
|
88.3
|
SO42−
|
5.4 (3.9)
|
0.3
|
2.6
|
4.3
|
7.0
|
32.2
|
SO42−
|
4.1 (3.1)
|
0.4
|
1.9
|
3.1
|
5.3
|
24.8
|
NO3−
|
4.6 (4.8)
|
0.1
|
1.4
|
3.0
|
6.3
|
41.1
|
NO3−
|
3.0 (3.0)
|
0.1
|
0.9
|
2.0
|
3.9
|
21.8
|
Cl−
|
0.6 (0.5)
|
0.0
|
0.2
|
0.5
|
0.8
|
3.0
|
Cl−
|
0.2 (0.2)
|
0.0
|
0.1
|
0.2
|
0.3
|
1.4
|
Na+
|
0.2 (0.1)
|
0.0
|
0.1
|
0.1
|
0.2
|
1.3
|
Na+
|
0.1 (0.1)
|
0.0
|
0.1
|
0.1
|
0.2
|
0.8
|
NH4+
|
3.7 (2.5)
|
0.2
|
1.9
|
3.1
|
4.7
|
19.3
|
NH4+
|
2.5 (1.8)
|
0.2
|
1.2
|
2.1
|
3.3
|
14.5
|
K+
|
0.2 (0.2)
|
0.0
|
0.1
|
0.1
|
0.2
|
3.1
|
K+
|
0.1 (0.2)
|
0.0
|
0.0
|
0.1
|
0.1
|
2.0
|
OC
|
4.1 (2.2)
|
0.2
|
2.5
|
3.7
|
5.3
|
17.5
|
OC
|
3.1 (1.6)
|
0.1
|
1.9
|
2.8
|
4.1
|
9.5
|
EC
|
1.1 (0.6)
|
0.1
|
0.8
|
1.0
|
1.4
|
3.7
|
EC
|
0.8 (0.4)
|
0.1
|
0.4
|
0.7
|
1.0
|
3.4
|
S
|
2.3 (1.7)
|
0.1
|
1.2
|
1.9
|
3.0
|
13.0
|
S
|
5.2 (3.9)
|
0.5
|
2.5
|
4.0
|
6.9
|
27.6
|
Ca
|
58.1 (81.8)
|
1.5
|
18.9
|
33.0
|
63.2
|
912.6
|
ng/㎥
|
Ca
|
51.9 (51)
|
1.5
|
23.5
|
36.3
|
61.4
|
581.4
|
Ti
|
8.5 (9.7)
|
0.5
|
3.6
|
5.6
|
9.7
|
103.9
|
Ti
|
9.1 (8.8)
|
0.5
|
4.0
|
6.8
|
11.2
|
89.8
|
V
|
3.5 (3.1)
|
0.3
|
1.3
|
2.8
|
4.8
|
23.5
|
V
|
7.7 (11)
|
0.3
|
1.1
|
3.1
|
9.1
|
86.6
|
Cr
|
0.8 (0.7)
|
0.3
|
0.3
|
0.3
|
1.1
|
7.1
|
Cr
|
2.0 (1.3)
|
0.3
|
1.0
|
1.7
|
2.8
|
9.0
|
Mn
|
11.9 (7.9)
|
0.3
|
6.4
|
10.3
|
15.5
|
52.6
|
Mn
|
18.6 (14.7)
|
0.3
|
7.9
|
15.3
|
25.3
|
131.0
|
Fe
|
177 (118.3)
|
12.1
|
103.4
|
149.1
|
216.9
|
1129.8
|
Fe
|
197 (128.4)
|
10.3
|
109
|
169.8
|
255.9
|
1177.6
|
Ni
|
1.3 (1.1)
|
0.2
|
0.5
|
1.0
|
1.7
|
9.8
|
Ni
|
3.0 (3.6)
|
0.2
|
0.8
|
1.6
|
3.6
|
25.8
|
Cu
|
4.0 (3.6)
|
1.0
|
2.1
|
3.3
|
5.1
|
46.3
|
Cu
|
6.7 (4.6)
|
1.0
|
3.8
|
5.7
|
8.5
|
56.0
|
Zn
|
57.1 (35.3)
|
2.4
|
32.2
|
48.7
|
73.2
|
230.4
|
Zn
|
64.3 (53.1)
|
0.6
|
29.9
|
52
|
81.6
|
487.5
|
As
|
3.7 (2.2)
|
0.3
|
2.1
|
3.3
|
4.7
|
16.0
|
As
|
5.9 (7.5)
|
0.3
|
1.5
|
3.2
|
7.1
|
64.9
|
Se
|
1.3 (1.1)
|
0.3
|
0.3
|
1.0
|
1.8
|
6.6
|
Se
|
1.5 (1.6)
|
0.3
|
0.3
|
1.0
|
2.1
|
10.8
|
Br
|
8.8 (6.1)
|
0.3
|
4.4
|
7.3
|
11.9
|
43.5
|
Br
|
8.0 (5.7)
|
0.4
|
4.1
|
6.5
|
10.6
|
55.4
|
Pb
|
19.7 (17.3)
|
0.7
|
7.5
|
15.0
|
26.9
|
115.4
|
Pb
|
17.7 (15.4)
|
0.7
|
7.1
|
13.4
|
23.6
|
121.7
|
During the study period, the daily average temperature (℃) ranged from -14.8 to 33.7 (SL), -12.7 to 33.4 (DJ), -9.5 to 32.0 (GJ), -8.2 to 31.8 (US), which shows substantial temporal variation in all cities while the five-year average temperatures (℃) were 13.3, 13.7, 14.7 and 14.7 in SL, DJ, GJ and US respectively. The daily average relative humidity (%) ranged 21.8-99.8, 27.5-99.3, 23.9-99.0 and 16.8-99.8 in SL, DJ, GJ and US respectively.
The average concentrations of PM2.5 mass at all sites exceeded both World Health Organization guidelines (5 µg/m3 of annual mean concentration) and U.S. Environment Protection Agency standards (15 µg/m3 of annual mean concentration). This suggests the population in four cities have been exposed to potential health risk during 2014-2018. It is notable that Ulsan showed approximately 25-30% lower PM2.5 concentration than the other cities.
PM2.5 chemical compositions showed considerable spatial variability. Secondary inorganic aerosols (SIA; i.e., SO42-, NO3-, NH4+) constructed the total mass at a range of 46.9% (US)-55.6% (SL), which indicates secondary aerosol formation were a important factor in PM2.5. In addition, the SO42-/NO3- ratio was much lower in SL (0.88) than those in DJ (0.97), GJ (1.16) and US (1.38), which implies the dominant gas-phase pollutants were different across the cities. Carbonaceous species, reported to have second largest proportions in PM2.5 mass in South Korea (Son et al. 2012; Bae et al. 2020), accounted for 18.2% (SL), 20.1% (DJ), 19.6% (GJ) and 18.4% (US) respectively. The proportions of trace elements were much higher in US (27.4%) than in SL (12.5%), DJ (14.5%) and GJ (10.0%), considered to reflect the industrial characteristics of US.
Source apportionment PM2.5
A total of 1,438, 1,278, 1,419 and 1,211 PM2.5 speciation datasets were used for identifying PM2.5 sources and their contributions for SL, DJ, GJ and US respectively. We chose the number of factors based on the evaluation of model results including Q-value, residual distribution and the coefficient of determination and we confirmed source profiles in the study cities were physically meaningful and understandable.
Nine-factor solution for DJ and ten-factor solution for SL, GJ and US solution were drawn while source contributions were quite different among cities although the resolved factors were identical. Sources and their contributions to PM2.5 are listed in Table 3 and source profiles for SL, DJ, GJ and US obtained from PMF model are displayed in Figures S1-S4 respectively.
Table 3
Sources with their contributions to PM2.5 in SL, DJ, GJ and US during 2014-2018
Source
|
SL
|
DJ
|
GJ
|
US
|
µg/m3
|
%
|
µg/m3
|
%
|
µg/m3
|
%
|
µg/m3
|
%
|
Secondary nitrate
|
6.4
|
20.3
|
5.7
|
19.5
|
5.2
|
19.5
|
3.6
|
18.9
|
Secondary sulfate
|
6.9
|
21.8
|
5.7
|
19.7
|
5.7
|
21.3
|
4.5
|
23.9
|
Mobile
|
7.1
|
22.3
|
6.3
|
21.7
|
5.5
|
20.5
|
4.7
|
24.9
|
Biomass burning
|
1.4
|
4.3
|
0.8
|
2.7
|
0.9
|
3.2
|
0.7
|
3.6
|
District heating
|
2.8
|
8.7
|
2.6
|
8.9
|
3.3
|
12.2
|
0.7
|
3.7
|
Soil
|
0.8
|
2.6
|
1.2
|
4.2
|
1.1
|
4.2
|
1.0
|
5.3
|
Industry
|
1.4
|
4.4
|
1.9
|
6.6
|
1.2
|
4.3
|
1.7
|
8.8
|
Coal combustion
|
3.4
|
10.7
|
3.1
|
10.7
|
2.6
|
9.7
|
0.7
|
3.5
|
Oil combustion
|
0.9
|
2.9
|
1.8
|
6.1
|
0.6
|
2.2
|
0.6
|
3.4
|
Aged sea salt
|
0.7
|
2.1
|
-
|
-
|
0.8
|
2.9
|
0.7
|
3.9
|
In SL, secondary nitrate, secondary sulfate and mobile source made up the majority (64.4%) of the PM2.5 mass concentration, consistent with the past studies conducted for SL (Heo et al. 2009; Moon et al., 2010). Mobile was the largest contributor followed by secondary sulfate, secondary nitrate and coal combustion, which implies sources related to secondary formation were dominant for PM2.5. Especially sulfur related sources, secondary sulfate and coal combustion, were thought to be affected by several huge coal-fired power plants (e.g., Taean 6100 MW, Dangjin 6040 MW, Yeongheung 5080 MW) located in the southwest part of the city. These results are supported by previous studies which also identified those coal power stations as attributable factors to PM2.5 in SL by analyzing polycyclic aromatic hydrocarbons (Kang et. al. 2020) and executing dispersion modelling (Kim et. al. 2016).
Similar to SL, PM2.5 mass concentration in DJ was mostly contributed by four factors; mobile (21.7%), secondary sulfate (19.7%), secondary nitrate (19.5%) and coal combustion (10.7%) while aged sea salt source was not solely resolved due to the location in deep land. In GJ secondary sulfate (21.3%) was the largest contributor followed by mobile (20.5%) and secondary nitrate (19.5%) source and those major factors in total explained 61.3% of PM2.5 mass concentration.
In US, the two largest factors were mobile (24.9%) and secondary sulfate (23.9%) sources which explained PM2.5 mass bit more than in other cities, assumed to reflect the characteristics of the city where both huge heavy industrial complex and industrial ports exist. The annual emission of SOx in the city during 2014-2017 was 48,433 ton on average more than two times of the national average annual emission (19,516 ton). In the same context the contribution of industry (8.8%) was much bigger than that in the other cities and it was fourth-largest factor in the city.
Association of PM 2.5 chemical constituents and source contributions with mortality
Using PMF-modeled source contributions as well as PM2.5 speciation data, we analyzed the effects of the PM2.5 chemical constituents of PM2.5 and sources on cause-specific mortality. Overall, the IQR increase in the concentration of each constituents influenced all cause-specific mortality whereas only a few constituents were significantly associated with the mortality with different strength depending on city and cause of death, which shows risk heterogeneity among various PM2.5 constituents. In addition, the impacts of source contributions on mortality also varied among the cities. The associations of cause-specific mortality with both PM2.5 constituents and sources are summarized in Tables S1-S2 respectively and Table 4 shows the highest RR of mortality significantly associated with chemical constituents.
Table 4
The highest RR and 95% CI of the significant associations of cause-specific mortality with PM2.5 constituents in SL, DJ, GJ and US.
City
|
Cause
|
Constituents
|
Lags
|
Constituents
|
Lags
|
Constituents
|
Lags
|
Constituents
|
Lags
|
SL
|
NA
|
K+ (1.016)
(1.002, 1.031)
|
7
|
Ti (1.011)
(1.003, 1.019)
|
7
|
EC (1.015)
(1.000, 1.031)
|
2
|
-
|
-
|
CVD
|
Ti (1.016)
(1.001, 1.032)
|
7
|
Fe (1.023)
(1.003, 1.044)
|
7
|
As (1.030)
(1.003, 1.058)
|
1
|
-
|
-
|
RD
|
Pb (1.071)
(1.013, 1.128)
|
0
|
-
|
-
|
-
|
-
|
-
|
-
|
DJ
|
NA
|
EC (1.051)
(1.004, 1.100)
|
2
|
Fe (1.039)
(1.009, 1.069)
|
2
|
-
|
-
|
-
|
-
|
CVD
|
OC (1.146)
(1.057, 1.239)
|
4
|
V (1.080)
(1.007, 1.154)
|
1
|
Mn (1.097)
(1.027, 1.167)
|
1
|
Ni (1.042)
(1.005, 1.080)
|
1
|
Zn (1.076)
(1.005, 1.148)
|
1
|
Pb (1.091)
(1.006, 1.177)
|
1
|
-
|
|
-
|
-
|
RD
|
OC (1.219)
(1.090, 1.354)
|
5
|
As (1.190)
(1.007, 1.384)
|
0
|
-
|
|
-
|
-
|
GJ
|
NA
|
Ti (1.015)
(1.001, 1.029)
|
3
|
-
|
-
|
-
|
-
|
-
|
-
|
RD
|
SO42− (1.085)
(1.009, 1.163)
|
1
|
NH4+ (1.071)
(1.000, 1.143)
|
1
|
Ti (1.048)
(1.009, 1.088)
|
5
|
V (1.109)
(1.040, 1.179)
|
1
|
Fe (1.061)
(1.001, 1.120)
|
5
|
Ni (1.100)
(1.033 1.170)
|
1
|
Zn (1.078)
(1.010, 1.146)
|
1
|
-
|
-
|
US
|
NA
|
SO42− (1.047)
(1.004, 1.091)
|
4
|
OC (1.056)
(1.001, 1.113)
|
4
|
Mn (1.060)
(1.018, 1.103)
|
0
|
-
|
-
|
CVD
|
OC (1.126)
(1.009, 1.248)
|
4
|
-
|
-
|
-
|
-
|
-
|
-
|
In SL, the chemical constituents showing significant associations were as follows: K+ (RR 1.016, 95% CI 1.002-1.032), EC (1.015, 1.000-1.031) and Ti (1.011, 1.003-1.019) with NA mortality, Ti (1.016, 1.001-1.032), Fe (1.023, 1.003-1.044) and As (1.030, 1.003-1.058) with CVD mortality and Pb (1.071, 1.013-1.128) with RD mortality. Ti and Fe, mostly originated from soil dust, are known to generate reactive oxygen species which lead to cell damage (Moreno et al. 2019). As and Pb, markers of coal combustion, are toxic even at low levels and they are known as carcinogen to human by International Agency for Research on Cancer. Significant associations of mortality with sources were mostly found for CVD mortality: district heating (1.032, 1.007-1.057), coal combustion (1.034, 1.003-1.065), soil (1.016, 1.000-1.032). Overall, the results in SL showed that the health effects were closely linked to transition metals from both natural and combustion-derived sources. Transition metals have impacts on cardiovascular via both direct and indirect pathways inducing oxidative stress and inflammation (Mills et al. 2009). In addition, respiratory system is subject to damage by those metals which trigger oxidative stress in lung alveoli (Donaldson et al. 1996).
In DJ we found PM2.5 constituents increased the RR of all types of mortality but most of the significant associations were related to CVD mortality. Constituents showing significant associations were EC (1.051, 1.004-1.100) and Fe (1.039, 1.009-1.069) with NA mortality, OC (1.146, 1.057-1.239), V (1.080, 1.007-1.154), Mn (1.097, 1.027-1.167), Ni (1.042, 1.005-1.080), Zn (1.076, 1.005-1.148) and Pb (1.091, 1.006-1.177) with CVD mortality and OC (1.219, 1.090-1.354) and As (1.190, 1.007-1.384) with RD mortality. Beside transition metals mentioned earlier, OC and EC have been consistently reported to influence on various indicators of cardiovascular disease such as blood pressure, heart rate variability (Huang et al. 2012; Schneider et al. 2012; Wu et al. 2013) which can lead to cardiovascular health outcomes (Bell et al. 2009; Ito et al. 2011). Accordingly, sources with carbonaceous species and transition metals were thought to be major factors associated with mortality in DJ.
Significant associations in GJ were mostly found for RD mortality, which was uniquely observed among the four cities. Chemical species significantly associated with RD mortality were ionic species, SO42- (1.085, 1.009-1.163), NH4+ (1.071, 1.000-1.143), as well as transition metals including Ti (1.048, 1.009-1.088), V (1.109, 1.040-1.179), Fe (1.061, 1.001-1.120), Ni (1.100, 1.033-1.170), Zn (1.078, 1.010-1.146). While several previous studies revealed the significant association SO42- with the risk of mortality, there are not explicit biological mechanisms of SO42- (Ueda et al. 2016). Nevertheless, some plausible theories have been suggested to explain the positive association of SO42- with adverse health effects. First, particle acidity by SO42- may change the pulmonary toxicity of other PM2.5 constituents or physical properties from their own toxicity (Dreher 2000). Second, catalyzation of metals into more bioavailable forms is another possible explanation for the high associations of SO42- with health effects (Lay et al. 2001). Regarding source contributions, sources of the species with significant association (e.g., coal combustion, mobile) were revealed to increase the RD mortality in GJ.
In US, SO42- (1.047, 1.004-1.091), OC (1.056, 1.001-1.113) and Mn (1.060, 1.018-1.103) for NA mortality and OC (1.126, 1.009-1.248) for CVD mortality increased the RR significantly. In addition, secondary sulfate (1.052, 1.007- 1.098), mobile (1.053, 1.010-1.098) and industry (1.059, 1.017-1.102) sources were strongly associated with NA or RD mortality. These results were in consistent with the source profiles from PMF modelling which identified higher contribution of mobile and secondary sulfate and the composition of PM2.5 compositional characteristic with the highest SO42-/NO3- ratio among four cities.