1. 华侨大学化工学院环境科学与工程系, 厦门 361021;
2. 核工业北京地质研究院分析测试研究中心, 北京 100029
收稿日期: 2016-01-08; 修订日期: 2016-03-29
基金项目: 国家自然科学基金项目(21377042,21477042);福建省自然科学基金项目(2014J01159);华侨大学2014级研究生科研创新能力培育项目(435)
Characteristics and Source Apportionment of Metals in the Dustfall of Quanzhou City
1. Department of Environmental Science and Engineering, College of Chemical Engineering, Huaqiao University, Xiamen 361021, China;
2. Center of Analysis, Beijing Research Institute of Uranium Geology, Beijing 100029, China
大气降尘是城市主要污染因子之一[1~3].随着城市化不断发展,城市大气降尘量及降尘中各种污染物不断增加.大气降尘可作为载体,将其携带的重金属等污染物输入到地表,影响陆地及水体生态系统[4].重金属不能被微生物降解,在环境中具有累积性,并可通过食物链威胁人们的健康与安全,对生态环境构成潜在威胁[5~7].
泉州市位于我国福建省东南沿海,是海峡西岸经济繁荣带和厦漳泉经济开发区的重要组成部分.随着泉州市经济的迅速发展,化石燃料的消耗量不断增加;同时,汽车保有量的激增使机动车尾气排放量剧增,泉州市空气质量呈下降趋势.目前对泉州市大气降尘中重金属元素的污染特征、生态风险评价等方面开展了一些研究,得出Cu、Zn、Pb、Cd、Cr等的污染较严重[8, 9].然而以上研究中涉及的金属元素种类不多且未对其具体来源及其贡献率开展系统地深入研究,为弥补这些方面的不足,本文分析了泉州市不同功能区大气降尘中的23种金属元素含量及铅锶同位素组成,利用富集因子法、地累积指数法和潜在生态风险指数法对金属元素富集程度、污染水平及生态风险进行评价,用多元统计分析结合铅锶同位素示踪进行污染源解析,以期为泉州市大气污染防治提供科学依据.
1 材料与方法
1.1 样品采集和预处理
2007年11月~2008年1月,在天气保持晴朗干燥至少3 d后,选择晴朗无大风的天气,在泉州市6大不同功能区布设采样点,共采集33个大气降尘样品,包括:交通繁忙区(J2~J4、J15、J3-1、J16),商业区(J8~J12、J14), 工业区(J17、J18、J22、J231~J233、J40、J41),居住区(J1、J5~J7、J13、J27、J30、J52、J53),风景区(J19、J34)和农业区(J38、J39).对每个采样点,在高度为1.5~2.0 m的一定范围内(半径约10 m的范围)用毛刷扫集玻璃、物件表面上的降尘于洁净信封中.所采集的样品在35℃烘干,过150目尼龙筛,筛下样装入洁净自封袋保存备用.
1.2 样品的分析测定
降尘中的Hg采用王水进行水浴消解,用原子荧光光谱仪(AFS,AF-640型,北京瑞利公司)测定;其它金属元素采用HCl+HNO3+HF微波消解,Fe、Mn用原子吸收光谱仪(AAS,TAS-986型,北京普析通用公司)测定,其它元素(Zn、Ni、Pb、Cu、Cr、V、Co、Cd、Sc、Li、Rb、Sr、Ba、Y、Bi、Th、U、Cs、Ti和Ca)用电感耦合等离子体质谱仪(ICP-MS,ELEMENT XR型,德国Finnigan-MAT公司)测定.实验过程中,每批样品处理时,均随机抽取一个样品作3次平行实验,同时做全程试剂空白,并以土壤标准样(ESS-3,中国环境监测总站)做全程质量监控. Hg测定结果的相对标准偏差(RSD)小于2.5%,Fe、Mn测定结果的RSD均小于5.0%,其它金属元素测定结果的RSD均小于10.0%.
Pb、Sr同位素比值分析中,样品的处理及测定均在核工业北京地质研究院分析测试研究中心同位素超净实验室中完成.同位素分离、纯化及测定方法见文献[10]. Pb、Sr同位素比值的测定精度分别优于0.005%、0.001%,满足Pb、Sr同位素示踪所要求的精度.
2 结果与讨论
2.1 泉州市大气降尘中金属元素含量水平
泉州市不同功能区大气降尘金属元素含量见表 1.降尘中不同金属元素含量有较大差异,Cd、Hg、Co、Ni、Zn、Pb、Ca、V、Cr、Cu、Bi等的变异系数(CV)较大,说明它们的来源较复杂. V、Ni、Zn、Sr、Cd、Bi、Fe在工业区的平均含量高于其他功能区,这主要与工业区内工业废气及粉尘的排放有关.我国目前还没有针对大气降尘中金属含量的标准,本文用福建省海岸带土壤元素背景值[11]作对比. Cd、Hg、Zn、Ca、Pb、Cu、Ni、Sr的平均含量分别是背景值的46.7、29.5、11.7、7.6、6.6、5.5、4.9和4.3倍,说明它们的富集程度相对较高.
表 1
(Table 1)
表 1 泉州市不同功能区大气降尘中金属元素含量/mg·kg-1
Table 1 Concentrations of metals in the dustfall from different functional areas in Quanzhou City/mg·kg-1
位置 |
样号 |
Li |
Sc |
Ti |
V |
Cr |
Mn |
Co |
Ni |
Cu |
Zn |
Rb |
Sr |
Y |
Cd |
Cs |
Ba |
Pb |
Bi |
Th |
U |
Fe |
Hg |
Ca |
交通繁忙区 |
J2 |
24.07 |
8.682 |
3 220 |
58.46 |
286.7 |
1 206 |
23.35 |
74.38 |
174.7 |
666.5 |
129.7 |
223.2 |
28.98 |
1.846 |
3.282 |
723.0 |
338.6 |
1.010 |
15.90 |
5.495 |
26 583 |
1.565 |
44 590 |
J3 |
27.46 |
8.702 |
3 045 |
62.59 |
104.4 |
944.0 |
14.07 |
64.01 |
111.3 |
1 108 |
124.7 |
240.4 |
27.44 |
2.045 |
4.169 |
713.4 |
553.1 |
1.617 |
13.37 |
4.931 |
22 251 |
4.163 |
64 420 |
J4 |
25.14 |
7.312 |
3 137 |
49.18 |
186.7 |
1124 |
12.07 |
102.1 |
224.0 |
722.4 |
113.6 |
249.8 |
22.09 |
1.691 |
3.273 |
1165 |
174.9 |
1.650 |
9.165 |
3.420 |
27 302 |
1.977 |
89 936 |
J15 | 23.60 |
9.216 |
3 375 |
61.70 |
510.1 |
1311 |
14.38 |
93.63 |
169.6 |
543.1 |
144.8 |
202.3 |
31.42 |
1.512 |
3.713 |
745.9 |
189.4 |
1.214 |
15.20 |
5.117 |
31 194 |
4.128 |
19 415 |
J3-1 |
24.00 |
9.577 |
3 167 |
51.30 |
108.0 |
896.7 |
147.0 |
78.04 |
211.9 |
521.8 |
122.9 |
205.6 |
30.82 |
2.417 |
3.571 |
705.3 |
328.2 |
1.721 |
15.96 |
5.230 |
24 205 |
0.602 |
64 420 |
J16 |
23.85 |
9.452 |
3 220 |
51.38 |
75.82 |
904.0 |
137.1 |
68.32 |
214.3 |
528.1 |
125.9 |
194.5 |
30.70 |
3.246 |
3.502 |
699.7 |
272.4 |
1.363 |
15.28 |
5.179 |
22 021 |
4.814 |
22 823 |
平均 |
24.69 |
8.824 |
3 194 |
55.77 |
212.0 |
1064 |
57.99 |
80.09 |
184.3 |
681.7 |
126.9 |
219.3 |
28.57 |
2.126 |
3.585 |
792.0 |
309.4 |
1.429 |
14.15 |
4.895 |
25 593 |
2.875 |
50 934 |
商业区 |
J8 |
32.78 |
8.741 |
3 380 |
63.32 |
195.7 |
1134 |
12.98 |
102.9 |
196.6 |
616.5 |
125.6 |
270.0 |
27.62 |
2.432 |
4.222 |
791.6 |
231.6 |
2.508 |
15.15 |
5.422 |
25 199 |
0.634 |
64 040 |
J9 |
25.30 |
8.562 |
3 248 |
65.08 |
122.5 |
1343 |
18.44 |
90.42 |
293.6 |
1 187 |
123.7 |
273.1 |
30.96 |
1.670 |
3.428 |
1039 |
468.8 |
2.167 |
15.62 |
5.598 |
40 693 |
5.113 |
67 170 |
J10 |
23.52 |
7.336 |
3 178 |
63.52 |
168.8 |
958.0 |
9.721 |
92.28 |
188.6 |
836.2 |
115.2 |
230.4 |
22.99 |
2.061 |
3.282 |
683.7 |
209.3 |
1.718 |
9.581 |
3.422 |
21 694 |
4.557 |
43 697 |
J11 |
33.47 |
8.114 |
2 995 |
64.00 |
162.0 |
977.5 |
16.79 |
84.15 |
186.3 |
647.7 |
137.1 |
239.3 |
25.27 |
1.729 |
4.498 |
848.9 |
161.0 |
2.964 |
13.32 |
4.350 |
29 565 |
0.130 |
32 700 |
J12 |
23.15 |
7.630 |
2 927 |
55.41 |
159.3 |
960.1 |
11.04 |
62.25 |
118.2 |
574.9 |
119.0 |
230.6 |
25.51 |
1.429 |
3.231 |
714.9 |
166.6 |
1.454 |
11.55 |
3.561 |
26 996 |
0.214 |
41 947 |
J14 |
28.76 |
9.598 |
2 947 |
57.22 |
105.1 |
1 124 |
120.5 |
78.46 |
163.3 |
484.2 |
122.8 |
276.1 |
24.77 |
4.522 |
3.387 |
599.8 |
250.8 |
3.022 |
9.404 |
4.233 |
20 854 |
4.985 |
74 696 |
平均 |
27.83 |
8.330 |
3 112 |
61.43 |
152.2 |
1 083 |
31.58 |
85.07 |
191.1 |
724.4 |
123.9 |
253.3 |
26.19 |
2.307 |
3.675 |
779.7 |
248.0 |
2.305 |
12.44 |
4.431 |
27 500 |
2.605 |
54 042 |
风景区 |
J19 |
12.74 |
7.305 |
2 516 |
47.10 |
10.15 |
904.9 |
6.410 |
3.684 |
15.53 |
101.4 |
140.9 |
134.6 |
20.78 |
2.222 |
1.608 |
477.7 |
40.88 |
0.802 |
22.70 |
5.740 |
24 590 |
0.543 |
5 558 |
J34 |
17.84 |
7.747 |
1 997 |
54.50 |
21.31 |
782.9 |
4.496 |
11.51 |
30.13 |
406.0 |
135.1 |
161.9 |
26.82 |
2.157 |
8.800 |
773.5 |
153.1 |
1.360 |
16.21 |
6.892 |
22 545 |
0.976 |
551 |
平均 |
15.29 |
7.526 |
2 257 |
50.80 |
15.73 |
843.9 |
5.453 |
7.596 |
22.83 |
253.7 |
138.0 |
148.2 |
23.80 |
2.190 |
5.204 |
625.6 |
96.99 |
1.081 |
19.46 |
6.316 |
23 568 |
0.760 |
3 054 |
工业区 |
J17 |
26.26 |
9.083 |
3 525 |
39.11 |
93.80 |
750.0 |
9.693 |
51.99 |
157.1 |
3357 |
103.4 |
341.4 |
19.92 |
5.281 |
4.939 |
857.0 |
204.9 |
3.528 |
26.78 |
4.455 |
28 026 |
1.970 |
8 098 |
J18 |
25.46 |
9.413 |
3 800 |
39.63 |
87.86 |
818.2 |
10.39 |
46.03 |
136.1 |
3569 |
100.6 |
333.6 |
21.60 |
5.199 |
4.493 |
868.0 |
213.1 |
3.487 |
26.86 |
6.156 |
42 028 |
1.450 |
63 944 |
J22 |
28.18 |
9.226 |
2 363 |
105.6 |
178.6 |
1 026 |
14.22 |
163.9 |
110.2 |
610.0 |
88.85 |
350.8 |
15.16 |
2.536 |
5.310 |
1 028 |
149.8 |
3.156 |
10.34 |
3.432 |
66 631 |
0.663 |
12 402 |
J231 |
21.24 |
10.48 |
4 462 |
252.0 |
237.5 |
1 300 |
27.49 |
318.6 |
121.6 |
1 725 |
163.8 |
729.4 |
36.10 |
3.494 |
2.301 |
1 688 |
833.0 |
3.102 |
20.46 |
7.486 |
145 584 |
0.918 |
9 056 |
J232 |
10.70 |
11.97 |
2 163 |
72.72 |
35.23 |
627.5 |
17.33 |
81.03 |
121.7 |
4 068 |
122.4 |
178.5 |
17.86 |
6.790 |
3.620 |
793.3 |
169.4 |
1.448 |
11.52 |
3.348 |
32 130 |
3.638 |
11 327 |
J233 |
27.56 |
13.95 |
2 539 |
235.4 |
81.37 |
1 354 |
25.33 |
248.2 |
117.7 |
3766 |
101.7 |
404.2 |
21.10 |
3.466 |
4.672 |
1 839 |
268.8 |
4.470 |
10.31 |
3.544 |
105 228 |
3.872 |
13 268 |
J40 |
28.70 |
10.09 |
1 947 |
58.30 |
25.33 |
1 346 |
8.060 |
44.81 |
32.19 |
369.4 |
29.75 |
258.6 |
17.95 |
1.390 |
2.254 |
569.3 |
89.84 |
2.056 |
4.563 |
1.916 |
48 010 |
0.503 |
20 100 |
J41 |
23.56 |
9.746 |
1 979 |
79.28 |
39.53 |
1 239 |
11.13 |
56.17 |
53.55 |
452.6 |
65.75 |
369.8 |
16.06 |
1.528 |
3.924 |
812.1 |
111.7 |
2.770 |
7.739 |
3.188 |
57 643 |
1.837 |
78 341 |
平均 |
23.96 |
10.49 |
2 847 |
110.2 |
97.40 |
1 058 |
15.45 |
126.3 |
106.3 |
2240 |
97.01 |
370.8 |
20.72 |
3.710 |
3.939 |
1 057 |
255.1 |
3.002 |
14.82 |
4.190 |
65 660 |
1.856 |
27 067 |
居住区 |
J27 |
25.40 |
8.786 |
3 044 |
157.2 |
149.9 |
1 200 |
9.588 |
149.1 |
204.9 |
708.0 |
160.5 |
518.8 |
32.46 |
3.820 |
9.086 |
1 866 |
573.6 |
5.972 |
19.46 |
7.120 |
125 489 |
2.384 |
14 494 |
J30 |
20.46 |
8.526 |
2 025 |
71.64 |
53.95 |
776.1 |
8.020 |
22.29 |
98.15 |
831.8 |
112.2 |
249.0 |
21.16 |
2.966 |
5.616 |
689.9 |
290.2 |
1.820 |
13.07 |
4.140 |
25 942 |
0.503 |
44 836 |
J13 |
28.77 |
6.986 |
2 424 |
49.00 |
307.5 |
700.0 |
15.03 |
232.6 |
155.9 |
816.9 |
115.8 |
211.9 |
19.95 |
3.954 |
4.051 |
911.2 |
291.0 |
1.264 |
8.265 |
3.432 |
7 585 |
4.730 |
12 768 |
J52 |
25.36 |
7.006 |
1 824 |
64.60 |
62.07 |
630.7 |
9.226 |
50.73 |
75.83 |
295.0 |
91.1 |
174.7 |
14.78 |
1.512 |
2.070 |
684.1 |
296.8 |
1.610 |
13.12 |
3.890 |
22 155 |
0.125 |
35 298 |
J53 |
24.24 |
9.026 |
1 859 |
85.32 |
68.41 |
1 405 |
12.08 |
65.15 |
73.59 |
604.6 |
95.7 |
257.2 |
17.44 |
2.726 |
5.198 |
922.7 |
389.4 |
3.456 |
12.84 |
4.720 |
58 850 |
0.186 |
31 850 |
J1 |
23.09 |
7.240 |
4 370 |
45.89 |
206.5 |
1 419 |
10.61 |
101.8 |
62.17 |
582.5 |
117.9 |
175.8 |
21.63 |
1.615 |
3.038 |
873.7 |
327.0 |
0.922 |
9.030 |
4.020 |
28 599 |
0.665 |
30 564 |
J7 |
24.12 |
6.867 |
2 852 |
47.01 |
90.09 |
954.3 |
9.481 |
41.27 |
59.42 |
448.0 |
124.6 |
153.1 |
19.91 |
1.290 |
3.770 |
693.1 |
143.8 |
0.995 |
9.004 |
3.893 |
17 912 |
1.195 |
20 817 |
J5 |
26.80 |
10.66 |
3 376 |
68.09 |
80.20 |
1 233 |
11.60 |
41.37 |
56.12 |
229.0 |
144.0 |
285.0 |
38.29 |
11.38 |
3.848 |
650.8 |
271.6 |
0.967 |
19.24 |
5.779 |
47 714 |
0.481 |
75 777 |
J6 |
25.98 |
9.702 |
3 309 |
60.42 |
71.76 |
1 282 |
13.82 |
48.47 |
69.90 |
265.3 |
153.5 |
227.0 |
34.87 |
1.040 |
3.417 |
719.0 |
118.5 |
0.868 |
18.18 |
5.659 |
44 726 |
0.740 |
28 122 |
平均 |
24.91 |
8.311 |
2 787 |
72.13 |
121.2 |
1 067 |
11.05 |
83.63 |
95.11 |
531.2 |
123.9 |
250.3 |
24.50 |
3.367 |
4.455 |
890.0 |
300.2 |
1.986 |
13.58 |
4.739 |
42 108 |
1.223 |
32 725 |
农业区 |
J38 |
28.24 |
10.55 |
2 698 |
40.30 |
16.33 |
1 209 |
4.608 |
15.39 |
15.53 |
227.2 |
77.09 |
131.2 |
11.49 |
0.468 |
3.402 |
518.1 |
51.92 |
1.574 |
10.75 |
2.690 |
31 253 |
0.601 |
9 053 |
J39 |
28.86 |
9.786 |
2 220 |
106.3 |
37.03 |
1 186 |
9.022 |
14.51 |
30.75 |
306.6 |
143.1 |
116.0 |
14.10 |
1.118 |
4.324 |
1 641 |
127.1 |
3.192 |
11.66 |
3.702 |
59 878 |
0.381 |
60.70 |
平均 |
28.55 |
10.17 |
2 459 |
73.31 |
26.68 |
1 198 |
6.815 |
14.95 |
23.14 |
266.90 |
110.08 |
123.61 |
12.80 |
0.793 |
3.863 |
1 080 |
89.53 |
2.383 |
11.20 |
3.196 |
45 566 |
0.491 |
4 557 |
土壤背景值 |
20.80 |
11.90 |
4 200 |
79.70 |
40.70 |
507.0 |
9.480 |
17.40 |
22.40 |
83.60 |
160.0 |
60.00 |
23.00 |
0.060 |
6.440 |
423.0 |
39.00 |
0.640 |
27.28 |
4.660 |
34 000 |
0.063 |
4 600 |
Min |
10.70 |
6.867 |
1 824 |
39.11 |
10.15 |
627.5 |
4.496 |
3.684 |
15.53 |
101.4 |
29.75 |
116.0 |
11.49 |
0.468 |
1.608 |
477.7 |
40.88 |
0.802 |
4.563 |
1.916 |
7 585 |
0.125 |
60.70 |
Max |
33.47 |
13.95 |
4 462 |
252.0 |
510.1 |
1419 |
147.0 |
318.6 |
293.6 |
4 068 |
163.8 |
729.4 |
38.29 |
11.38 |
9.086 |
1 866 |
833.0 |
5.972 |
26.86 |
7.486 |
145 584 |
5.113 |
89 936 |
Mean |
24.81 |
9.001 |
2 883 |
75.23 |
125.4 |
1061 |
23.79 |
84.83 |
122.7 |
975.0 |
118.0 |
260.6 |
23.88 |
2.805 |
4.039 |
888.0 |
256.4 |
2.158 |
13.99 |
4.580 |
41 245 |
1.856 |
35 035 |
CV/% |
20.03 |
17.21 |
23.39 |
65.22 |
59.10 |
21.56 |
88.66 |
84.52 |
55.70 |
113.3 |
23.17 |
47.37 |
27.88 |
74.58 |
38.41 |
40.13 |
61.42 |
53.94 |
37.50 |
28.94 |
73.93 |
83.70 |
68.56 |
|
表 1 泉州市不同功能区大气降尘中金属元素含量/mg·kg-1
Table 1 Concentrations of metals in the dustfall from different functional areas in Quanzhou City/mg·kg-1
|
2.2 降尘中金属元素污染评价
2.2.1 富集因子法
富集因子(EF)可用以评价元素富集程度,也可用来区分污染来源.其计算公式为:
|
(1) |
式中,Ci为金属元素i的含量(mg·kg-1),Cn为参比元素的含量(mg·kg-1).根据参比元素选择条件[12],可用作参比元素的有Al、Zr、Li、Ti、Sc、Cs、Ca、Fe、Rb等,本研究以Li、Ti、Sc、Cs、Fe为候选参比元素,对候选参比元素与其余元素之间进行相关性分析,发现Ti、Fe与其余元素的相关性相对较好,而Ti的CV相对更小,且其化学稳定性较好,因此本文选Ti作为参比元素.若大气颗粒物中某元素的EF接近于1,可认为土壤尘(自然源)是主导源;若EF>5,则可认为非土壤源(即人类活动)对其有显著贡献[13].
表 2为泉州市各功能区大气降尘中金属元素的平均EF值. Cd和Hg在各功能区的EF值均较高(13.15~93.74),说明这两元素富集较严重. Cd在不同功能区大气降尘中的EF值顺序:工业区(93.74)>居住区(86.75)>风景区(68.71)>商业区(52.30)>交通繁忙区(46.70)>农业区(23.70),最大值出现在居住区的采样点J5(235.95),这可能与该采样点附近大量建筑活动(建筑设备和施工车辆)有关.风景区两个采样点位于清源山,Cd的EF值偏高可能与山中交通干道产生的扬尘和汽车尾气有关. Hg在不同功能区大气降尘中的EF值顺序:交通繁忙区(59.91)>商业区(55.60)>工业区(48.50)>居住区(30.28)>风景区(23.49)>农业区(13.15). Zn在各功能区也均有富集,其平均EF值顺序:工业区(39.68)>商业区(11.62)>交通繁忙区(10.81)>居住区(10.64)>风景区(6.14)>农业区(5.58). Pb在除农业区外的其他功能区均有富集,尤其是居住区和交通繁忙区(EF>10). Cu、Ca、Ni在交通繁忙区、商业区、工业区和居住区均存在一定富集(EF值6.05~15.82).此外,Cr在交通繁忙区(6.74)和商业区(5.04),Sr在商业区(5.70)、工业区(9.26)和居住区(6.70),Bi在工业区(7.25)、居住区(5.27)和农业区(6.63)存在一定富集;Co只在交通繁忙区存在一定富集(EF=8.05).其余元素在各功能区均没有富集.不同功能区的大气降尘中金属元素富集程度有所不同.总体来看,农业区受到的人为污染程度最小.
表 2
(Table 2)
表 2 泉州市各功能区大气降尘金属富集因子
Table 2 Enrichment factors of metals in the dustfall from different functional areas of Quanzhou City
元素 |
交通繁忙区 |
商业区 |
风景区 |
工业区 |
居住区 |
农业区 |
Li |
1.56 |
1.81 |
1.41 |
1.85 |
1.95 |
2.37 |
Sc |
0.97 |
0.95 |
1.20 |
1.43 |
1.13 |
1.47 |
V |
0.92 |
1.04 |
1.21 |
2.10 |
1.49 |
1.66 |
Cr |
6.74 |
5.04 |
0.76 |
3.35 |
4.56 |
1.17 |
Mn |
2.75 |
2.88 |
3.11 |
3.41 |
3.30 |
4.07 |
Co |
8.05 |
4.64 |
1.06 |
2.51 |
1.88 |
1.28 |
Ni |
6.05 |
6.57 |
0.87 |
10.69 |
7.60 |
1.48 |
Cu |
10.81 |
11.43 |
1.99 |
7.04 |
6.96 |
1.84 |
Zn |
10.81 |
11.62 |
6.12 |
39.68 |
10.64 |
5.58 |
Rb |
1.04 |
1.05 |
1.62 |
0.90 |
1.22 |
1.22 |
Sr |
4.82 |
5.70 |
4.71 |
9.26 |
6.70 |
3.53 |
Y |
1.63 |
1.54 |
1.98 |
1.36 |
1.64 |
0.97 |
Cd |
46.70 |
52.30 |
68.71 |
93.74 |
86.75 |
23.70 |
Cs |
0.73 |
0.77 |
1.65 |
0.98 |
1.13 |
1.05 |
Ba |
2.47 |
2.48 |
2.87 |
3.82 |
3.37 |
4.62 |
Pb |
10.53 |
8.52 |
5.00 |
8.77 |
12.75 |
4.12 |
Bi |
2.95 |
4.88 |
3.28 |
7.25 |
5.27 |
6.63 |
Th |
0.68 |
0.61 |
1.32 |
0.76 |
0.80 |
0.71 |
U |
1.38 |
1.28 |
2.58 |
1.30 |
1.62 |
1.20 |
Fe |
0.99 |
1.09 |
1.30 |
2.93 |
1.94 |
2.38 |
Hg |
59.91 |
55.60 |
23.49 |
48.50 |
30.28 |
13.15 |
Ca |
14.74 |
15.82 |
1.13 |
9.90 |
11.55 |
1.54 |
|
表 2 泉州市各功能区大气降尘金属富集因子
Table 2 Enrichment factors of metals in the dustfall from different functional areas of Quanzhou City
|
2.2.2 地累积指数法
Müller[14]于1969年提出地累积指数(Igeo)法,用于研究水体沉积物中金属元素污染程度,现被广泛用于沉积物、土壤及街道灰尘的金属污染评价研究,其计算公式:
|
(2) |
式中,Ci为降尘中重金属元素i的实测含量;Bi为i的环境背景值,本研究采用福建省海岸带土壤元素背景值[11].考虑造岩运动引起的背景值波动设定了常数k(一般取k=1.5).
根据Igeo值,元素污染水平可分为如下6级: Igeo < 0,无污染;0≤Igeo < 1,轻污染;1≤Igeo < 2,偏中度污染;2≤Igeo < 3,中度污染;3≤Igeo < 4,偏重度污染;4≤Igeo < 5,重度污染;Igeo≥5,严重污染.
图 1为泉州市大气降尘中各金属元素的Igeo分布.将各元素平均Igeo值与Igeo分级标准对比可知,Th、Cs、Ti、Rb、Sc、V、U、Y、Fe、Li没有受到污染;Co、Ba、Mn、Cr、Bi污染程度较轻;Ni、Sr、Cu、Ca、Pb受到偏中度污染;Zn受到了中度污染;Hg则污染偏重,Cd达到了重度污染级别,其中部分采样点降尘中Cd、Hg、Zn的Igeo值超过5,达到严重污染程度.
2.2.3 潜在生态危害指数法
潜在生态危害指数法是Hakanson[15]根据毒物性质及环境行为特点,从沉积学角度提出的对土壤或沉积物中毒物污染进行评价的方法.根据这一方法,某一区域大气降尘中第i种重金属的潜在生态危害系数(Eri)及沉积物中多种重金属的潜在危害指数(RI)可分别表示为:
|
(3) |
|
(4) |
式中,Ci为大气降尘中金属元素i的实测浓度;Cni为元素i的背景值,本文选用福建省海岸带土壤元素背景值;Tri为i的毒性系数,反映其毒性水平和生物对其敏感程度.采用该法评价泉州市不同功能区大气降尘中V、Cr、Mn、Co、Ni、Cu、Zn、Cd、Pb、Hg、Ti的潜在生态危害水平,上述各元素的毒性系数分别取2、2、1、5、5、5、1、30、5、40、1[16].
本研究根据元素种类及其毒性系数对各生态风险程度所对应的RI限值,依据文献[8]进行适当调整后的评价标准见表 3.
表 3
(Table 3)
表 3 潜在生态风险评价指标和等级划分
Table 3 Evaluation index and grade division of potential ecological risk
单元素潜在生态 危害系数(Eri) |
多元素潜在生态 危害指数(RI) |
潜在生态 风险程度 |
<40 |
<110 |
生态危害低 |
40~80 |
110~220 |
生态危害中等 |
80~160 |
220~440 |
生态危害较高 |
160~320 |
≥440 |
生态危害高 |
≥320 |
|
生态危害极高 |
|
表 3 潜在生态风险评价指标和等级划分
Table 3 Evaluation index and grade division of potential ecological risk
|
泉州市大气降尘中各元素的平均Eri表现为Cd(1402)>Hg(1178)>Pb(32.87)>Cu(27.40)>Ni(24.38)>Co(12.55)>Zn(11.66)>Cr(6.16)>Mn(2.09)>V(1.89)>Ti(0.69).各元素对大气降尘中重金属潜在生态风险指数贡献见图 2.可见,Cd和Hg的贡献较大(分别为51.93%和43.63%).不同功能区大气降尘重金属的RI值范围为640.1~6073.1,均达到了极高生态风险程度,不同功能区大气降尘中多元素综合潜在生态风险顺序为:工业区(3171.7)>交通繁忙区(3045.6)>商业区(2943.9)>居住区(2566.6)>风景区(1607.0)>农业区(742.1).
2.3 金属元素污染来源分析
2.3.1 聚类分析
聚类分析是对类似的对象进行分类组合的分析过程,可用于重金属来源判别.采用离差平方和法(Ward's method)和欧氏距离法(Euclidean distance)将泉州市大气降尘中金属元素进行R型变量聚类分析,结果如图 3所示.降尘中金属可分为三类:第一类为V、Fe、Ba、Bi、Ni、Sr、Pb、Cs、Sc、Zn、Cd,其中V、Ni多用于工业生产[17](如电镀、合金生产和钢铁制造等),Zn为冶金尘特征元素,同时橡胶轮胎的磨损也会产生Zn[18],Pb和Cd是机动车尾气中的重要元素[19].故第一类代表工业和交通污染源;第二类为Th、U、Rb、Y、Ti,变异系数小且几乎都无富集,代表土壤自然源;第三类为Li、Mn、Cu、Hg、Cr、Co、Ca,其中Cr、Mn、Hg、Co是典型的燃煤排放物[20, 21],故第三类代表燃煤源.综上可知,工业、交通、燃煤和土壤尘是泉州市大气降尘中金属元素的主要来源.
2.3.2 同位素示踪污染来源
按前述讨论,泉州市工业区、交通繁忙区、商业区和居住区大气降尘中重金属的综合潜在生态风险相对更高,Pb、Sr也主要在上述4个功能区相对富集,本研究测定了上述4功能区中共20个降尘样品的Pb、Sr同位素组成(见表 4),并结合潜在污染源的Pb、Sr同位素特征,分析泉州市大气降尘中Pb、Sr主要来源.
表 4
(Table 4)
表 4 泉州大气降尘中Pb、Sr同位素组成
Table 4 Isotope ratios of Pb and Sr in the dustfall of Quanzhou City
功能区 |
样号 |
206Pb/207Pb |
206Pb/208Pb |
87Sr/86Sr |
交通繁忙区 |
J2 |
1.1614 |
0.47417 |
0.710330 |
J3 |
1.1602 |
0.47397 |
0.710362 |
J4 |
1.1687 |
0.47545 |
0.710183 |
商业区 |
J8 |
1.1714 |
0.47665 |
0.710165 |
J10 |
1.1673 |
0.47539 |
0.710244 |
J11 |
1.1700 |
0.47642 |
0.710373 |
J14 |
1.1658 |
0.47472 |
0.709986 |
工业区 |
J17 |
1.1723 |
0.47651 |
0.710852 |
J18 |
1.1752 |
0.47760 |
0.707550 |
J22 |
1.1712 |
0.47588 |
0.709485 |
J232 |
1.1759 |
0.47612 |
0.709316 |
J40 |
1.1653 |
0.47448 |
0.709910 |
J41 |
1.1430 |
0.47097 |
0.709409 |
居住区 |
J27 |
1.1659 |
0.47456 |
0.710336 |
J30 |
1.1701 |
0.47581 |
0.710646 |
J13 |
1.1641 |
0.47467 |
0.710160 |
J53 |
1.1646 |
0.47364 |
0.710316 |
J5 |
1.1642 |
0.47435 |
0.710194 |
J6 |
1.1916 |
0.48515 |
0.710159 |
J7 |
1.1690 |
0.47598 |
0.710652 |
|
表 4 泉州大气降尘中Pb、Sr同位素组成
Table 4 Isotope ratios of Pb and Sr in the dustfall of Quanzhou City
|
(1)铅同位素来源示踪不同铅源具有不同的Pb同位素比值特征,且比值不受环境中物理化学过程的影响,因此Pb同位素比值已被广泛应用于物源识别[22].潜在污染源组分的数据来自文献[23],它们之间的206Pb/207Pb、206Pb/208Pb相差较大,因此可有效示踪泉州市大气降尘铅污染来源.以206Pb/208Pb对206Pb/207Pb作图(见图 4),发现大气降尘与燃煤尘和土壤母质(自然源)接近,而离汽车尾气尘相对较远.工业区样品J41(206Pb/207Pb和206Pb/208Pb值分别为1.1430、0.47097)完全落在燃煤尘范围,因该采样点位于南埔火电厂燃烧设备附近,明显受到燃煤影响.居住区样品J6位于郊区,靠近土壤母质,说明该点的大气降尘主要受控于土壤扬尘.
定量解析是目前污染源解析的一大难题.铅同位素指纹分析技术在定量解析铅污染源时具有独特的优势[24].根据Monna等[25]提出的三元混合模型,计算燃煤、汽车尾气和土壤母质对泉州市不同功能区大气降尘铅来源的贡献率.计算公式如下:
|
(5) |
|
(6) |
|
(7) |
式中,RS和NS分别为降尘样品的206Pb/208Pb和206Pb/207Pb;R1、R2、R3分别为3个主要污染源的206Pb/208Pb;N1、N2、N3分别为3个主要污染源的206Pb/207Pb;f1、f2、f3对应为3个主要污染源的贡献率(%).计算结果见表 5.从中可见,泉州市大气降尘中的Pb主要来源于土壤母质和燃煤尘,平均贡献率分别为45.45% (范围29.41%~64.00%)和45.13% (范围22.53%~60.48%);汽车尾气平均贡献率为9.42% (范围3.13%~13.47%).
表 5
(Table 5)
表 5 不同铅来源相对贡献率的计算结果
Table 5 Calculation result of relative contribution of different lead sources/%
功能区 |
样号 |
汽车尾气 |
燃煤 |
自然源 |
交通繁忙区 |
J2 |
10.29 |
57.71 |
32.00 |
J3 |
10.64 |
59.94 |
29.41 |
J4 |
8.80 |
43.51 |
47.69 |
平均值 |
9.91 |
53.72 |
36.37 |
商业区 |
J8 |
12.59 |
32.20 |
55.21 |
J10 |
10.23 |
44.72 |
45.04 |
J11 |
12.98 |
34.85 |
52.16 |
J14 |
8.10 |
50.96 |
40.94 |
平均值 |
10.98 |
40.68 |
48.34 |
工业区 |
J17 |
10.67 |
32.99 |
56.34 |
J18 |
13.47 |
22.53 |
64.00 |
J22 |
8.23 |
38.69 |
53.08 |
J232 |
3.80 |
34.36 |
61.84 |
J40 |
7.24 |
53.17 |
39.59 |
平均值 |
8.68 |
36.35 |
54.97 |
居住区 |
J27 |
6.94 |
52.25 |
40.81 |
J30 |
9.24 |
39.83 |
50.93 |
J13 |
9.88 |
52.22 |
37.90 |
J53 |
3.13 |
60.48 |
36.39 |
J5 |
7.87 |
54.82 |
37.31 |
J7 |
11.58 |
38.97 |
49.45 |
平均值 |
8.11 |
49.76 |
42.13 |
|
总均值 |
9.42 |
45.13 |
45.45 |
|
表 5 不同铅来源相对贡献率的计算结果
Table 5 Calculation result of relative contribution of different lead sources/%
|
(2)锶同位素来源示踪与铅同位素比值类似,锶同位素比值也可作为环境指示剂[26].利用所研究对象中87Sr/86Sr和1/Sr作图,可以很好地示踪其中Sr的来源[27].泉州市大气降尘87Sr/86Sr和1/Sr分布如图 5所示,所有样品点远离土壤母质,大部分样品点87Sr/86Sr值落在燃煤尘和汽车尾气中87Sr/86Sr值范围内,表明交通和燃煤是泉州市大气降尘中锶的主要来源,而几乎不受土壤母质的影响.
3 结论
(1)泉州市大气降尘中不同金属元素含量空间分布差异较大. Cd、Hg、Zn、Ca、Pb、Cu、Ni、Sr的富集程度及污染水平相对较高.所有功能区的多元素综合潜在生态风险程度极高, 由高到低顺序为工业区>交通繁忙区>商业区>居住区>风景区>农业区;Cd和Hg的潜在生态危害程度极高, 二者对综合潜在生态风险指数的贡献达95.56%.
(2)聚类分析结果表明: V、Fe、Ba、Bi、Ni、Sr、Pb、Cs、Sc、Zn、Cd主要来源于工业污染和交通排放, Th、U、Rb、Y、Ti主要来源于土壤尘, Li、Mn、Cu、Hg、Cr、Co、Ca主要来源于燃煤.铅同位素三元混合模型计算得到降尘中铅主要受土壤母质(范围29.41%~64.00%, 平均45.45%)、燃煤尘(范围22.53%~60.48%, 平均45.13%)和汽车尾气(范围3.13%~13.47%, 平均9.42%)的影响;87Sr/86Sr和1/Sr分布图显示降尘中锶来源主要为燃煤尘和汽车尾气.