环境科学  2024, Vol. 45 Issue (6): 3502-3511   PDF    
基于APCS-MLR模型的开封市公交站周边灰尘重金属源解析及健康风险评估
段海静1,2,3, 申浩欣1, 彭超月1, 任翀1, 王艳锋1, 刘德新1,2,3, 王玉龙1,2,3, 郭瑞超1,2,3, 马建华1,2,3     
1. 河南大学地理与环境学院, 开封 475004;
2. 河南大学黄河中下游数字地理技术教育部重点实验室, 开封 475004;
3. 河南大学土壤重金属污染控制与修复工程研究中心, 开封 475004
摘要: 为揭示城市交通系统对城市生态环境质量的影响, 选择受交通影响扰动较强的环境指示物——公交站地表灰尘作为研究对象, 利用电感耦合等离子体质谱仪(ICP-MS)和电感耦合等离子体发射光谱仪(ICP-ASE)分别测定灰尘中8种重金属(V、Cr、Co、Ni、Cu、Zn、Cd和Pb)含量, 应用地累积指数和污染负荷指数法分析灰尘重金属的污染程度和空间分布特征, 通过定性(相关性分析、主成分分析)和定量[绝对因子得分-多元线性回归模型(APCS-MLR)]相结合的方法开展公交站附近地表灰尘重金属来源探讨, 应用克里格空间插值法明晰重金属不同来源的空间分布特征, 利用美国环境保护署提出的健康风险评价模型来评价人体健康风险. 结果表明, 开封市公交站地表灰尘重金属ω(V)、ω(Cr)、ω(Co)、ω(Ni)、ω(Cu)、ω(Zn)、ω(Cd)和ω(Pb)平均值依次为:68.36、59.73、5.81、19.34、40.10、208.32、1.01和49.46 mg·kg-1, 灰尘中重金属(Cd、Zn、Pb、Cu、Cr)含量平均值均高于开封市周边灰尘背景值, 分别是背景值的3.37、2.70、2.01、1.95和1.28倍;8种重金属的地累积指数顺序为:Cd > Zn > Pb > Cu > Cr > V > Ni > Co, 其中Cd、Zn、Cu和Pb属于轻度污染水平, 其他元素为无污染;源解析结果显示, Cr、Co和Ni为自然源元素, Cu、Zn、Pb和Cd为交通源元素, V则来源于工业-自然混合源. APCS-MLR结果表明, 4种来源的平均贡献率依次为自然源为34.17%, 交通源为29.84%, 工业-自然混合源为14.64%, 未知源为21.35%, 其中交通源贡献率空间分布规律与交通量及公交线路密集度分布趋势一致. 由健康风险评价可知, 儿童的总致癌风险指数和总非致癌风险指数均高于成人, Cr为主要的非致癌因子, Cd为主要的致癌因子, 自然源和交通源分别对非致癌风险和致癌风险的贡献率最高.
关键词: 灰尘重金属      污染评价      源解析      健康风险      公交车站     
Source Apportionment and Health Risk Assessment of Heavy Metals in Dust Around Bus Stops in Kaifeng City Based on APCS-MLR Model
DUAN Hai-jing1,2,3 , SHEN Hao-xin1 , PENG Chao-yue1 , REN Chong1 , WANG Yan-feng1 , LIU De-xin1,2,3 , WANG Yu-long1,2,3 , GUO Rui-chao1,2,3 , MA Jian-hua1,2,3     
1. College of Geography and Environmental Science, Henan University, Kaifeng 475004, China;
2. Key Laboratory of Geospatial Technology for the Middle and Lower Yellow River Regions, Ministry of Education, Henan University, Kaifeng 475004, China;
3. Engineering Research Center for Control & Remediation of Soil Heavy Metal Pollution, Henan University, Kaifeng 475004, China
Abstract: In order to reveal the influence of urban transportation systems on the quality of urban ecological environment, this study selected surface dust from bus stops, which is strongly disturbed by transportation, as the research object. The contents of eight heavy metals (V, Cr, Co, Ni, Cu, Zn, Cd, and Pb) in the dust were determined through inductively coupled plasma mass spectrometry (ICP-MS) and inductively coupled plasma atomic emission spectroscopy (ICP-ASE). The spatial distribution characteristics and pollution levels of the eight heavy metals in the dust were analyzed using the geo-accumulation index method. A combined qualitative (correlation analysis and principal component analysis) and quantitative (absolute principal component scores-multiple linear regression model (APCS-MLR)) method was used to explore the sources of heavy metals in surface dust near bus stops. The spatial distribution characteristics of heavy metals from different sources were elucidated using the Kriging interpolation method. The health risk assessment model proposed by the United States Environmental Protection Agency was used to evaluate the human health risks. The results showed that the average values of ω(V), ω(Cr), ω(Co), ω(Ni), ω(Cu), ω(Zn), ω(Cd), ω(Pb), and ω(As) in the bus stop surface dust were 68.36, 59.73, 5.81, 19.34, 40.10, 208.32, 1.01, and 49.46 mg·kg-1, respectively. The concentrations of heavy metals (Cd, Zn, Pb, Cu, and Cr) in the dust were all higher than the background values in the surrounding dust, exceeding them by 3.37, 2.70, 2.01, 1.95, and 1.28 times, respectively. The order of the geo-accumulation index for the eight heavy metals was Cd > Zn > Pb > Cu > Cr > V > Ni > Co, with Cd, Zn, Cu, and Pb in the dust indicating mild pollution levels and the others showing no pollution. The source analysis results showed that Cr, Co, and Ni were natural sources, whereas Cu, Zn, Pb, and Cd were traffic sources, and V was derived from a combination of industrial and natural sources. The APCS-MLR results indicated that the average contribution rates of the four sources were as follows: natural source (34.17%), traffic source (29.84%), industrial-natural mixed source (14.64%), and unknown source (21.35%). The spatial distribution map of the contribution rate of the traffic source was consistent with the trends of traffic volume and bus route density distribution. According to the health risk assessment, the cancer risk and non-cancer risk for children were both higher than those for adults. Cr was the main non-cancer factor, and Cd was the main cancer-causing factor. Natural and traffic sources contributed the most to non-cancer risk and cancer risk, respectively.
Key words: dust heavy metal      pollution assessment      source apportionment      health risk      bus stops     

地表灰尘(灰尘)是一种物质组成和来源复杂的固体颗粒物, 其粒径(介于1 ~ 1 000 μm)较小[1], 比表面积较大, 易吸附重金属、多环芳烃等污染物. 由于灰尘受外力作用影响会反复发生“沉降-扬起”等行为, 极易通过皮肤接触、手口摄入和呼吸吸入等暴露途径进入人体, 危害人类健康, 引起了学者们的广泛关注[2 ~ 5]. 有研究显示[6, 7], 汽车尾气排放、车体的自身磨损及路面材料老化等均会释放大量含有重金属的污染物, 尤其是车辆拥堵及行进缓慢的区域, 频发的刹车、启动等行为更易出现重金属污染峰值[8]. 城市公交站作为车辆停靠的重要场所[9, 10], 其周边灰尘重金属污染现状对城市生态环境质量具有重要的指示意义[11]. 随着我国经济的飞速发展, 汽车拥有量逐年攀升, 交通拥堵现象时有发生, 机动车辆运行对城市生态环境扰动强烈[12], 人类长期接触重金属会损害肺、肾脏或其他器官, 易产生脑血管和心血管等疾病[13], 因此开展公交站地表灰尘重金属污染特征研究, 对于城市生态环境质量改善和人群健康具有重要的现实意义.

我国学者对城市公交客运系统周围地表灰尘重金属开展了相关研究[12, 14 ~ 20]:田东凡等[12]和周潇[14]分别对西安市汽车站、客运站地表灰尘重金属进行了污染评估及来源识别;李军等[15]研究了兰州快速公交系统(BRT)沿线站台灰尘重金属含量特征、污染程度、健康风险及可能的污染源. 大部分学者则重点关注公交站附近地表灰尘:如王可等[16]利用生物提取实验法模拟人体胃肠道环境, 研究石家庄市公交站点灰尘重金属的生物可给性及对人体的健康风险;Wu等[17]通过定性源解析识别西峰市公交站地表灰尘重金属来源;李小飞等[18]研究了福州市公交枢纽站地表灰尘重金属含量、来源及健康风险;武永福[19]应用定性源解析法识别庆阳市公交站地表灰尘重金属来源. 仅李春艳等[20]运用正定矩阵分解模型(PMF)定量解析天水市公交站地表灰尘重金属的4类源:自然源、自然-交通混合源、建筑废弃物污染源和燃煤-交通混合源. 综上所述, 有关公交客运系统周围地表灰尘重金属研究主要围绕重金属含量[12, 14 ~ 20]、生物可给性[16]、污染特征[12, 14, 15, 17, 19, 20]、健康风险[14 ~ 16, 18]以及以相关性分析[14, 18]、主成分分析[14, 19]、聚类分析[19]和因子分析[19]为主的定性源解析, 较少开展以PMF[14, 20]、绝对因子得分-多元线性回归模型(APCS-MLR)[21]和多变量受体模型(UNMIX)[22]等定量源解析. 定性源解析仅能识别污染源种类, 无法给出污染源贡献率, 更无法定量污染源对人体健康风险的贡献率. 精准识别污染来源是科学有效控制污染物的重要手段, 将定量源解析与健康风险模型相结合的灰尘重金属风险评估方法, 有利于确定优先控制污染源, 精准制定污染管控措施, 对于降低管理部门资金投入、有效控制灰尘污染物对人体健康影响至关重要. 因此, 开展以污染源为导向的重金属健康风险评估模型已成为学者们定量识别污染风险的重要手段[2, 4].

开封市作为国务院首批公布的历史文化名城之一, 市内旅游景点达17处之多, 景点位置相对集中, 大多位于城墙以内的老城区. 城区内道路狭窄, 车辆拥堵严重, 人流、车流密集, 其城市生态环境质量引起了学者们的广泛关注[23 ~ 28]. 但以往对开封市细颗粒物研究主要集中于公园[23]、驾校[24]和学校[25, 26]等封闭空间的地表灰尘、大气降尘[27]和叶面尘[28]等, 较少开展开放性空间(公交车站)灰尘重金属的系统研究. 因此, 本研究拟以开封市公交站地表灰尘重金属为例, 运用地累积指数和污染负荷指数法评估地表灰尘重金属污染特征, 开展定性定量源解析, 同时将APCS-MLR与健康风险评价模型相结合对开封市公交站地表灰尘重金属进行源导向健康风险评估, 以期为开封城市生态环境质量提升和人群健康管理提供科学支撑.

1 材料与方法 1.1 研究区概况

开封市位于河南省中部偏东, 紧邻省会城市郑州, 其经纬度为34°11′45″~35°01′20″N, 113°52′15″~115°15′42″E. 截至2022年, 市区人口密度由2011年的723人·km-2增长至2022年的767人·km-2, 2021年接待游客人数达7 296.5万人次[29]. 研究区主要位于开封市主城区内:以城墙为界, 城墙北侧为文教区(黄河水利职业技术学院、河南大学和开封文化职业技术学院等高校);城墙南侧滨河路和陇海铁路间为火车站、汽车中心站等客运中心, 也是车流量、人流量密集区;城墙外侧东南方向主要为老工业区(分布着化肥厂、化工厂、流量仪表厂、炼锌厂和药厂等企业);城墙西侧为城市未来发展区, 城墙西南侧为新兴工业区;城墙内侧偏北为开封市文化旅游核心区——宋都皇城旅游度假区(龙亭公园、清明上河园、翰园碑林和万岁山森林公园等景区), 而城墙内侧西门大街和中山路则为核心商业区, 其他区域多为商业和居民混合区. 开封市拥有郑开轻轨、高铁、高速路、郑开大道和郑汴物流通道等完备的交通网络, 极大地推动了当地旅游业的快速发展. 现阶段, 市内公共交通网络的运营线路63条, 公交运营线路总里程693 km, 拥有公交车852辆[30], 极大地提高了公共交通的通行速度. 开封属于温带季风气候, 冬季寒冷干燥, 夏季高温多雨, 年均温度为14℃左右, 年均降雨量约为670 mm. 开封城市表层土壤主要为潮土, 其pH介于4.08 ~ 8.43之间, 均值为7.62, 土壤呈弱碱性[31].

1.2 样品采集与处理

根据车流量、人群活动情况及居民区分布情况, 灰尘样点主要布设于金明大道、大梁路、宋城路、汴京路和中山路等城市主干道沿线的公交站点附近(图 1). 采样前一周无降雨, 采样时天气晴朗无风. 在每个公交站点附近, 用塑料铲和毛刷采集2个灰尘子样, 然后将每个站所收集的灰尘进行充分混合作为该公交站的代表样, 每份样品重约100 g, 放入聚乙烯自封袋中带回实验室, 共获得灰尘样品91个.

图 1 采样点分布示意 Fig. 1 Geographical distributions of sampling sites

样品自然风干后, 用镊子挑拣出植物残体、石子砖块等杂质, 用研钵研磨后过100目(0.149 mm)尼龙筛, 按编号装入聚乙烯袋中. 样品采用HNO3-HClO4-HF体系消解, 利用电感耦合等离子体质谱仪(ICP-MS, Xseries 2)和电感耦合等离子体原子发射光谱仪(ICAP-ASE)测定优先控制重金属(V、Cr、Co、Ni、Cu、Zn、Cd和Pb)的含量[32]. 在分析过程中加入黄棕壤(GSS-3)进行质量控制[23].

1.3 评价方法 1.3.1 灰尘重金属污染评价方法

运用Muller提出的地累积指数(Igeo)评价地表灰尘重金属的污染状况, 计算公式[33]如下:

(1)

式中, Ci表示灰尘中重金属含量(mg·kg-1);Si表示灰尘背景值(mg·kg-1).

污染负荷指数(PLI)是Tomlinson等[34]提出的能够反映重金属污染水平的评价方法, 计算公式[34]如下:

(2)
(3)

式中, CFi为灰尘重金属i的污染系数, Ci为灰尘重金属i的实测值, Cn为灰尘重金属n的背景值, PLI为某样点多种重金属的污染负荷指数.

1.3.2 APCS-MLR模型

APCS-MLR模型是利用主成分分析(PCA)获取绝对主成分因子得分(APCS), 再将APCS作为自变量, 重金属元素含量为因变量进行多元线性回归分析, 得出不同污染源贡献率, 其公式[21]如下:

(4)

式中, Cii种重金属的实测含量;bi为多元回归的常数项, aim为污染源mi种重金属的回归系数;APCSim为污染源mi种重金属的绝对主因子得分, n为因子个数.

(5)

式中, PCim为污染源m对重金属元素i的贡献率;|bi|为多元回归常数项的绝对值;为污染源mi种重金属绝对主因子得分的平均值.

(6)

式中, PCim为未被识别来源(未知源)的贡献率.

1.3.3 健康风险评价

本研究利用美国环境保护署提出的健康风险评价模型[35]分析开封市公交站地表灰尘重金属的健康风险. 通过手口直接摄入、呼吸吸入和皮肤接触这3种途径进入人体, 各种暴露途径的暴露量计算公式如下[36]

(7)
(8)
(9)

式中, C为重金属的实测含量(mg·kg-1), ADDing、ADDinh、ADDderm分别代表手口直接摄入、呼吸吸入和皮肤接触途径的重金属暴露量, mg·(kg·d)-1, 参数见文献[19, 23, 35, 37 ~ 40].

人体健康风险分为致癌风险和非致癌风险, 其计算公式如下[41]

(10)
(11)

式中, HI为总非致癌风险指数, HQij、ADDij和RfDij分别为第i种污染物通过第j种途径的非致癌风险单项指数(风险商)、暴露量及参考剂量;TCR为总致癌风险指数, CRij和SFij分别为单一致癌风险指数和斜率因子, 参数见文献[17, 23].

2 结果与讨论 2.1 灰尘重金属含量

从开封公交站灰尘8种重金属元素含量可知(表 1), ω(V)、ω(Cr)、ω(Co)、ω(Ni)、ω(Cu)、ω(Zn)、ω(Cd)和ω(Pb)平均值依次为68.36、59.73、5.81、19.34、40.10、208.32、1.01和49.46 mg·kg-1. 与开封市周边灰尘背景值[42]相比, 除V、Co和Ni低于背景值外, 其他重金属(Cd、Zn、Pb、Cu和Cr)平均值分别为背景值的3.37、2.70、2.01、1.95和1.28倍. 与黄河下游潮土区土壤元素背景值[43]相比, 除Ni低于背景值外, 其他重金属(Cd、Zn、Pb、Cu和Cr)平均值分别为背景值的11.22、3.20、2.43、1.87和1.11倍. 5种重金属(V、Co、Ni、Cu和Pb)均未超过《土壤环境质量 农用地土壤污染风险管控标准(试行)》(GB 15618-2018)(pH > 7.5)的风险筛选值[44], Cd为风险筛选值的1.68倍, 但未超过风险管制值, 表明Cd存在一定的污染风险. 变异系数(CV)由高到低依次为:V > Cr > Cd > Cu > Zn > Pb > Ni > Co, 其中, Co属于低等变异水平, Ni和Pb属于中等变异水平, 而V、Cr、Cu、Zn和Cd则达到高度变异水平, 且V、Cr、Cu和Cd的变异系数超过了50%, 表明其空间分异性较大, 受环境影响较强烈[45].

表 1 公交站地表灰尘重金属描述性统计1) Table 1 Descriptive statistics of heavy metals in surface dust of bus stops

2.2 灰尘重金属污染特征

利用式(1)以开封市周边灰尘背景值为参比计算开封市公交站地表灰尘重金属地累积指数(Igeo), 结果见表 2. Igeo排序依次为:Cd(0.96) > Zn(0.74) > Pb(0.34) > Cu(0.20) > Cr(-0.35) > V(-1.04) > Ni(-1.06) > Co(-1.34). Cr、V、Ni和Co的平均Igeo均 < 0, 处于无污染等级, Cd、Zn、Pb和Cu的平均Igeo介于0 ~ 1之间, 处在轻度污染水平. Cd的累积程度最大, 有40.66%的样点为偏中度污染, 有2个采样点处于偏重度污染, 占比为2.20%;对于Zn, 35.16%采样点处于偏中度污染, 1.10%处于中度污染;Pb处在轻度污染的样点数占比最高, 为68.13%, Co所有样点均为无污染水平.

表 2 地表灰尘重金属的Igeo1) Table 2 The Igeo of heavy metals in surface dust

利用式(3)以开封市周边灰尘背景值为参比计算开封市公交站地表灰尘重金属污染负荷指数(PLI), 由研究区灰尘重金属PLI空间分布可知(图 2), 除东南部存在中度污染水平之外, 其余地区均为无污染及轻度污染水平, 所有样点的PLI介于0.80 ~ 2.39之间, 其中属于无污染样点数9个, 占比为9.89%, 轻度污染样点数80个, 占比为87.91%, 中度污染样点数2个, 占比为2.20%. 研究区东南部的PLI较高, 可能是该区域为开封市老工业区与老旧小区叠加的区域[28]. 从整体上看, PLI从西北部向东南部递增, 呈带状分布.

图 2 污染负荷指数空间分布 Fig. 2 Spatial distribution of pollution load index

2.3 灰尘重金属定性源解析 2.3.1 相关性分析

利用皮尔逊相关系数对地表灰尘中重金属进行相关性分析(表 3), 在0.01水平下, V和Cd的相关性系数为0.592, 表明具有一定的同源性;Cr、Co、Ni两两之间显著正相关, Zn、Pb、Cd呈显著相关性, Cu和Zn、Pb呈显著正相关, 显著相关性表明它们之间有相似的地球化学行为或相同的来源.

表 3 各重金属之间的相关系数1) Table 3 Correlation coefficients among heavy metals

2.3.2 主成分分析

利用SPSS 26对数据进行标准化处理, 并且进行KMO检验和巴特利特球度检验, 巴特利特球度检验结果为0.000(小于0.05), KMO度量值检验结果为0.51(大于0.5), 说明数据可以进行主成分分析. 以特征值大于1的主成分为主因子, 本数据共提取3个主因子(表 4), 且累计方差贡献率为75.064%.

表 4 地表灰尘重金属元素含量的旋转成分矩阵 Table 4 Rotated component matrix of heavy metal elements contents in surface dust

主因子1(PC1)的方差贡献率为37.080%, 旋转载荷高的重金属为Cr、Co和Ni, 其权重系数分别为0.876、0.785和0.870. Co和Ni通常与成土母质有关[14, 46], 由表 1可知, Co和Ni平均含量低于背景值, 变异系数分别为14.24%和22.89%, 表明受人为影响因素较小;研究表明[14], Cr也受成土母质影响, Cr的平均值与背景值相近, 且Igeo显示Cr为无污染, 表明受人为影响较小, 所以Cr、Co和Ni元素的来源可能和自然有关, 因此因子1为自然源.

主因子2(PC2)的方差贡献率为22.476%, 其中Cu、Zn和Pb的旋转载荷较高, 其权重系数分别为0.768、0.797和0.796. 重金属Cu来源于汽车零部件的磨损[47];Zn被大量应用于润滑剂和轮胎的生产中[48, 49], 机动车尾气排放或轮胎磨损会带来大量含有Zn的灰尘;李军等[15]对兰州BRT沿线站台灰尘分析认为Pb与交通活动有关, Pb也与轮胎和刹车磨损有关[7], 因此因子2为交通源.

主因子3(PC3)的方差贡献率为15.508%, 其中载荷高的是V和Cd, 其权重系数分别为0.935和0.740. 研究表明[6], 生产汽车零部件和路面交通标志线的制作、煤炭石油燃烧等都会产生Cd, 从而推断Cd来源于交通活动. V的产生和化石燃料燃烧有关[22, 50], 研究区内存在老工业区, 推断因子3为交通-工业混合源.

2.4 APCS-MLR分析

利用式(5)和式(6)计算开封市公交站地表灰尘重金属源贡献率, 运用APCS-MLR模型识别开封市公交站地表灰尘重金属来源(图 3), 并结合克里格空间插值法得到来源贡献率分布(图 4). 由图 3可知, 源1中的Cr、Co和Ni的贡献率最高, 分别为60.16%、50.91%和82.90%. 由图 4(a)可知, 源1贡献率在27.17% ~ 51.77%, 空间分布相对均匀, 东北出现的峰值可能与采样点S35的Cr含量较高有关, 需进一步研究, 且Cr、Co、Ni来源和成土母质有关, 因此推断源1为自然源.

图 3 地表灰尘重金属污染源贡献率和平均贡献率 Fig. 3 Contribution rate and average contribution rate of heavy metal pollution sources in surface dust

图 4 地表灰尘重金属污染源贡献率空间分布 Fig. 4 Spatial distribution of contribution rate of heavy metal pollution sources in surface dust

图 3可知, Cu、Zn、Pb和Cd在源2中有较高的贡献率, 分别为73.79%、57.50%、55.64%和32.23%, 由图 4(b)可知, 源2贡献率在15.96% ~ 35.90%之间, 高值区在中部偏东和南部, 其峰值与主干道以及公交车的运行方向一致, 高值区位于城区入口及火车站和客运站附近, 且Cu、Zn、Pb和Cd受到交通活动影响, 因此推断源2为交通源.

V是源3的主要贡献元素, 贡献率为69.92%, 由图 4(c)可知, 源3贡献率在8.68% ~ 21.44%之间, V的高值区分布在东南部, 该区域与老工业区重叠, 结合图 2(研究区东南部的PLI较高), 可初步识别V主要受工业活动影响. 由图 4(c)还可知, 研究区北部存在V的次高区, 该区域主要位于万岁山森林公园、西北湖湿地公园、龙庭公园等景区附近, 也有研究表明[13], V也受成土母质和自然风化影响;因此初步判定V可能受工业活动和自然因素的双重影响, 推断源3为工业-自然混合源.

图 4(d)为未知源, 其贡献率在14.53% ~ 32.87%之间, 高值区主要分布于文教区, 空间分布特征受人类活动影响不明显. 由图 3可知, 未知源对Cr、Co、Cd和Pb贡献率分别为37.08%、46.28%、30.74%和35.52%. 由图 3还可知, 污染源平均贡献率从大到小依次是:自然源(34.17%) > 交通源(29.84%) > 未知源(21.35%) > 工业-自然混合源(14.64%).

2.5 重金属健康风险评估 2.5.1 健康风险评估

利用式(7)~(11)计算开封市公交站地表灰尘重金属健康风险评估, 结果见表 5表 6. 儿童8种重金属的非致癌风险单项指数(HQ)依次为:Cr > V > Pb > Cd > Cu > Ni > Zn > Co, 成人与其一致, 且在3种暴露途径下, HQCr均为最大, 说明Cr为主要的非致癌因子, 儿童和成人各重金属的HI均小于1, 表明不存在非致癌风险. 各重金属不同暴露途径的HQ儿童均大于成人, 表明儿童更容易受到影响, 除Co外, 重金属通过3种暴露途径对不同人群所产生的非致癌风险从大到小均为:手口摄入 > 皮肤接触 > 呼吸吸入, 表明手口摄入是非致癌风险的主要途径.

表 5 地表灰尘重金属不同途径非致癌风险平均值 Table 5 Average non-carcinogenic risks of heavy metals in surface dust in different ways

表 6 地表灰尘重金属不同途径致风险平均值1) Table 6 Average carcinogenic risks of heavy metals in surface dust in different ways

鉴于部分元素缺乏SF值, 因此在进行致癌风险评估时, Cr、Co、Ni仅考虑呼吸吸入途径的致癌风险, Pb仅考虑皮肤接触途径的致癌风险. 由表 6可知, 手口摄入是致癌风险的主要途径, 成人和儿童的手口摄入单项致癌风险指数(CR手口摄入)小于1×10-4, 说明致癌风险在可接受范围内. 5种重金属对成人产生的单项致癌风险指数(CR)均值为:Cd > Cr > Co > Ni > Pb, 儿童为:Cd > Cr > Co > Pb > Ni, 说明Cd为主要的致癌因子, 其中Cd和Pb所产生的CR均为儿童高于成人. CRCr、CRCo、CRNi和CRPb远小于1×10-6, 表明无致癌风险, CRCd在1×10-6 ~ 1×10-4之间, 存在可耐受的致癌风险. 成人和儿童的平均TCR分别为1.80×10-6和3.50×10-6, 表明公交站地表灰尘重金属对成人和儿童存在可耐受的致癌风险, 且儿童是成人的1.94倍. 以往的研究也多是儿童的非致癌和致癌风险高于成人, 儿童比成人对环境污染物的敏感性更高[15], 因此更应该关注儿童的卫生习惯, 儿童应勤洗手, 避免重金属通过手口途径摄入体内.

2.5.2 污染源对健康风险率评估

将APCS-MLR模型与健康风险评价模型相结合, 计算得到不同污染源对健康风险的贡献率[46, 51], 结果见图 5. 对于非致癌风险(NCR), 不同来源对成人和儿童贡献率均为:自然源 > 未知源 > 工业-自然混合源 > 交通源, 自然源对非致癌风险贡献率最大, 这与Heidari等[52]的研究结果一致. 自然源贡献最大可能因为Cr、Co和Ni的RfD值较小, 使得非致癌风险较高[53], Cr是主要的非致癌因子, 未知源对Cr也有37.08%的贡献率, 因此除自然源对Cr贡献大外, 其他源亦有贡献.

图 5 不同污染源对成人和儿童非致癌风险(NCR)和致癌风险(TCR)的贡献率 Fig. 5 Contribution rates of different pollution sources to the non-carcinogenic risks (NCR) and carcinogenic risks (TCR) of adults and children

对于致癌风险(TCR), 成人和儿童的4个污染来源贡献率大小均为:交通源 > 未知源 > 工业-自然混合源 > 自然源, 交通源对致癌风险贡献率最大, 其对成人和儿童的致癌风险贡献率分别为31.00%和31.85%. 这与其他学者的研究结果稍有差异, 马杰等[4]对农用地土壤研究结果表明, 农业源对致癌风险贡献率最大;Huang等[54]对武汉市公园灰尘重金属健康风险评价定量分析结果显示, 自然源对致癌风险贡献率最大;Heidari等[52]研究阿巴斯市郊区道路粉尘表明自然源对非致癌风险和致癌风险贡献率最大. 学者们的结果不同与研究区域周围环境有关, 重金属来源不同, 污染源对健康风险的贡献率亦不同, 这与以往的研究结果一致[51].

3 结论

(1)开封市公交车站地表灰尘重金属含量ω(V)、ω(Cr)、ω(Co)、ω(Ni)、ω(Cu)、ω(Zn)、ω(Cd)和ω(Pb)平均值依次为:68.36、59.73、5.81、19.34、40.10、208.32、1.01和49.46 mg·kg-1, Cd、Zn、Pb、Cu和Cr含量均高于开封市周边灰尘背景值, 其中Cd污染最严重为背景值的3.37倍, 是首要的污染因子.

(2)Cd的地累积指数最大, 处于轻度污染水平. 公交站重金属污染整体处在轻度污染水平, 从西北部向东南部递增, 呈带状分布.

(3)Cr、Co和Ni以自然源为主, 方差贡献率为37.080%;Cu、Zn和Pb以交通源为主, 方差贡献率为22.476%;V和Cd以交通-工业混合源为主, 方差贡献率为15.508%.

(4)开封公交站附近地表灰尘重金属元素主要来自于自然源、交通源、工业-自然混合源和未知源, 其平均贡献率分别为34.17%、29.84%、14.64%和21.35%. 自然源高值区分布在东北部, 交通源高值区与人流量及车流量行进轨迹吻合, 工业-自然混合源高值区位于景区及火车站、汽车站等人流量密集区域, 未知源变化规律不明显.

(5)成人和儿童均不存在非致癌风险, 儿童的致癌风险和非致癌风险均高于成人, Cr为主要的非致癌因子, Cd为主要的致癌因子, 不同污染来源对非致癌风险和致癌风险贡献率依次为:自然源 > 未知源 > 工业-自然混合源 > 交通源和交通源 > 未知源 > 工业-自然混合源 > 自然源.

参考文献
[1] Yap C K, Chew W Y, Tan S G. Heavy metal concentrations in ceiling fan and roadside car park dust collected from residential colleges in Universiti Putra Malaysia, Serdang, Selangor[J]. Pertanika Journal of Tropical Agricultural Science, 2012, 35(1): 75-83.
[2] Wang J H, Wu H H, Wei W X, et al. Health risk assessment of heavy metal(loid)s in the farmland of megalopolis in China by using APCS-MLR and PMF receptor models: taking Huairou district of Beijing as an example[J]. Science of the Total Environment, 2022, 835. DOI:10.1016/j.scitotenv.2022.155313
[3] Orosun M M, Oniku S A, Peter A, et al. Magnetic susceptibility measurement and heavy metal pollution at an automobile station in Ilorin, North-Central Nigeria[J]. Environmental Research Communications, 2020, 2(1). DOI:10.1088/2515-7620/ab636a
[4] 马杰, 葛淼, 王胜蓝, 等. 基于源导向的农用地土壤重金属健康风险评估及优先控制因子分析[J]. 环境科学, 2024, 45(1): 396-406.
Ma J, Ge M, Wang S L, et al. Health risk assessment and priority control factors analysis of heavy metals in agricultural soils based on source-oriented[J]. Environmental Science, 2024, 45(1): 396-406.
[5] Liu Y, Jin T, Yu S H, et al. Pollution characteristics and health risks of heavy metals in road dust in Ma'anshan, China[J]. Environmental Science and Pollution Research, 2023, 30(15): 43726-43739. DOI:10.1007/s11356-023-25303-2
[6] 刘玥, 郭文强, 武晔秋. 大同市城区公园表层土壤重金属污染特征和健康风险及来源解析[J]. 土壤通报, 2023, 54(1): 180-191.
Liu Y, Guo W Q, Wu Y Q, et al. Pollution characteristics, health risks and source analysis of surface soil heavy metals in urban parks of Datong city[J]. Chinese Journal of Soil Science, 2023, 54(1): 180-191.
[7] Liu Z, Du Q Q, Guan Q Y, et al. A Monte Carlo simulation-based health risk assessment of heavy metals in soils of an oasis agricultural region in northwest China[J]. Science of the Total Environment, 2023, 857(3). DOI:10.1016/j.scitotenv.2022.159543
[8] Ewen C, Anagnostopoulou M A, Ward N I. Monitoring of heavy metal levels in roadside dusts of Thessaloniki, Greece in relation to motor vehicle traffic density and flow[J]. Environmental Monitoring and Assessment, 2009, 157(1-4): 483-498. DOI:10.1007/s10661-008-0550-9
[9] 马嘉怡. 呼和浩特市公共交通现状与发展对策研究[D]. 呼和浩特: 内蒙古农业大学, 2022.
Ma J Y. Research on the current situation and development countermeasures of public transport in Hohhot[D]. Hohhot: Inner Mongolia Agricultural University, 2022.
[10] 张春勤. 城市公共交通服务治理模式的比较研究[D]. 上海: 上海交通大学, 2015.
Zhang C Q. Comparative research on the governance mode of urban public transport service[D]. Shanghai: Shanghai Jiao Tong University, 2015.
[11] Pal S K, Roy R. Evaluating heavy metals emission' pattern on road influenced by urban road layout[J]. Transportation Research Interdisciplinary Perspectives, 2021, 10. DOI:10.1016/j.trip.2021.100362
[12] 田东凡, 卢新卫, 景一鸣, 等. 西安市汽车站地表灰尘中重金属的污染评价及其来源识别[J]. 水土保持通报, 2017, 37(1): 34-39.
Tian D F, Lu X W, Jing Y M, et al. Assessment of heavy metal pollution and source identification in surface dust at urban coach stations in Xi'an city[J]. Bulletin of Soil and Water Conservation, 2017, 37(1): 34-39.
[13] Huang J L, Wu Y Y, Sun J X, et al. Health risk assessment of heavy metal(loid)s in park soils of the largest megacity in China by using Monte Carlo simulation coupled with positive matrix factorization model[J]. Journal of Hazardous Materials, 2021, 415. DOI:10.1016/j.jhazmat.2021.125629
[14] 周潇. 西安市客运站地表灰尘中重金属与多环芳烃污染特征及风险研究[D]. 西安: 陕西师范大学, 2019.
[15] 李军, 李开明, 位静, 等. 兰州BRT沿线站台灰尘及其两侧绿化带土壤重金属污染及健康风险评价[J]. 地球与环境, 2022, 50(2): 228-240.
Li J, Li K M, Wei J, et al. Contaminations and health risks of heavy metals from the roadside greenbelt soils and dust along the BRT platform in Lanzhou[J]. Earth and Environment, 2022, 50(2): 228-240.
[16] 王可, 左书梅, 李昆, 等. 石家庄市公交站点尘土重金属的生物可给性及人体健康风险评估[J]. 科学技术与工程, 2020, 20(35): 14737-14742.
Wang K, Zuo S M, Li K, et al. Bioaccessibility and health risk assessment of heavy metals in the dusts in bus stops in Shijiazhuang city[J]. Science Technology and Engineering, 2020, 20(35): 14737-14742. DOI:10.3969/j.issn.1671-1815.2020.35.057
[17] Wu Y F, Lu X W, Zhuang S K, et al. Contamination characteristics and assessment of manganese, zinc, chrome, lead, copper and nickel in bus station dusts of Xifeng, northwest China[J]. Preprints, 2016. DOI:10.20944/preprints201608.0137.v1
[18] 李小飞, 陈志彪, 张永贺, 等. 福州市公交枢纽站地表灰尘重金属含量、来源及其健康风险评价[J]. 环境科学研究, 2013, 26(8): 906-912.
Li X F, Chen Z B, Zhang Y H, et al. Concentrations, sources and health risk assessments of heavy metals in ground surface dust from urban bus terminals of Fuzhou city[J]. Research of Environmental Sciences, 2013, 26(8): 906-912.
[19] 武永福. 黄土高原塬区公交站灰尘污染物富集特征及风险评估——以庆阳市为例[D]. 西安: 陕西师范大学, 2019.
[20] 李春艳, 王新民, 王海, 等. 天水市主城区公交站地表灰尘重金属来源解析及污染评价[J]. 环境科学, 2024, 45(1): 417-428.
Li C Y, Wang X M, Wang H, et al. Source appointment and assessment of heavy metals pollution in surface dust in the main district bus stops of Tianshui city[J]. Environmental Science, 2024, 45(1): 417-428.
[21] 段海静, 马嘉玉, 彭超月, 等. 基于APCS-MLR和PMF模型解析黄河下游文化公园土壤重金属污染特征及来源分析[J]. 环境科学, 2023, 44(8): 4406-4415.
Duan H J, Ma J Y, Peng C Y, et al. Pollution characteristics and source analysis of heavy metals in soils in Yellow River Cultural Park based on APCS-MLR and PMF receptor model[J]. Environmental Science, 2023, 44(8): 4406-4415.
[22] 吴静, 方凤满, 马康, 等. 基于Unmix模型的庐江县土壤重金属空间分布和来源解析[J]. 生态与农村环境学报, 2022, 38(9): 1204-1210.
Wu J, Fang F M, Ma K, et al. Spatial distribution and source analysis of soil heavy metals in Lujiang county based on Unmix model[J]. Journal of Ecology and Rural Environment, 2022, 38(9): 1204-1210.
[23] 段海静, 蔡晓强, 阮心玲, 等. 开封市公园地表灰尘重金属污染及健康风险[J]. 环境科学, 2015, 36(8): 2972-2980.
Duan H J, Cai X Q, Ruan X L, et al. Assessment of heavy metal pollution and its health risk of surface dusts from parks of Kaifeng, China[J]. Environmental Science, 2015, 36(8): 2972-2980.
[24] 陈轶楠, 马建华. 河南省某市驾校地表灰尘重金属污染及健康风险评价[J]. 环境科学学报, 2016, 36(8): 3017-3026.
Chen Y N, Ma J H. Assessment of pollution and health risks of heavy metals in surface dusts from driving schools in a city of Henan, China[J]. Acta Scientiae Circumstantiae, 2016, 36(8): 3017-3026.
[25] 王晓云, 马建华, 侯千, 等. 开封市幼儿园地表灰尘重金属积累及健康风险[J]. 环境科学学报, 2011, 31(3): 583-593.
Wang X Y, Ma J H, Hou Q, et al. Accumulation and health risk assessment of heavy metals in kindergarten surface dust in the city of Kaifeng[J]. Acta Scientiae Circumstantiae, 2011, 31(3): 583-593.
[26] 王洪涛, 闫帅腾, 武锟鹏, 等. 开封市某高校校园地表灰尘重金属分布特征及生态风险评估[J]. 土壤通报, 2023, 54(1): 171-179.
Wang H T, Yan S T, Wu K P, et al. Distribution characteristics and ecological risk assessment of heavy metals in surface dust from a university campus in Kaifeng, China[J]. Chinese Journal of Soil Science, 2023, 54(1): 171-179.
[27] 张忠地, 邵天杰, 卫佩茹. 开封市大气降尘中重金属污染及潜在生态危害评价[J]. 干旱区资源与环境, 2019, 33(11): 156-162.
Zhang Z D, Shao T J, Wei P R. Assessment of heavy metal pollution and potential ecological hazards in atmospheric dust in Kaifeng city[J]. Journal of Arid Land Resources and Environment, 2019, 33(11): 156-162.
[28] 段海静, 邢戎光, 王亚琪, 等. 开封市不同功能区叶面尘重金属含量及潜在生态风险[J]. 环境科学学报, 2019, 39(12): 4226-4234.
Duan H J, Xing R G, Wang Y Q, et al. Content and potential ecological risk of heavy metals in foliar dust of different functional of Kaifeng city[J]. Acta Scientiae Circumstantiae, 2019, 39(12): 4226-4234.
[29] 开封统计局. 开封统计年鉴2022[R]. 开封: 开封统计局, 2022.
[30] 刘杨. 构建绿色出行网络 打造绿色交通城市[N]. 开封日报, 2022-09-21.
[31] 李灿. 开封城市土壤性质、重金属污染及变化分析[D]. 开封: 河南大学, 2007.
Li C. The urban soil properties, heavy metal pollution and its changes in Kaifeng city[D]. Kaifeng: Henan University, 2007.
[32] HJ 766-2015, 固体废物 金属元素的测定 电感耦合等离子体质谱法[S].
HJ 766-2015, Solid waste-determination of metals-inductively coupled plasma mass spectrometry[S].
[33] Muller G. Index of geoaccumulation in sediments of the Rhine River[J]. Geojournal, 1969, 2(108): 108-118.
[34] Tomlinson D L, Wilson J G, Harris C R, et al. Problems in the assessment of heavy-metal levels in estuaries and the formation of a pollution index[J]. Helgoländer Meeresuntersuchungen, 1980, 33(1): 566-575.
[35] USEPA (United States Environmental Protection Agency). Risk assessment guidance for superfund (RAGS): part A[R]. Washington DC: USEPA, 1989.
[36] 金梓函, 张军. 宝鸡市不同时期地表灰尘重金属空间分布及其健康风险研究[J]. 环境科学与管理, 2022, 47(11): 172-177.
Jin Z H, Zhang J. Spatial distribution and health risk of heavy metals in surface dust of different periods in Baoji city[J]. Environmental Science and Management, 2022, 47(11): 172-177. DOI:10.3969/j.issn.1673-1212.2022.11.034
[37] Ferreira-Baptista L, De Miguel E. Geochemistry and risk assessment of street dust in Luanda, Angola: a tropical urban environment[J]. Atmospheric Environment, 2005, 39(25): 4501-4512. DOI:10.1016/j.atmosenv.2005.03.026
[38] 耿雅妮, 梁青芳, 杨宁宁, 等. 宝鸡市城区灰尘重金属空间分布、来源及健康风险[J]. 地球与环境, 2019, 47(5): 696-706.
Geng Y N, Liang Q F, Yang N N, et al. Distribution, sources and health risk assessment of heavy metals in dusts of the urban area of the Baoji city[J]. Earth and Environment, 2019, 47(5): 696-706.
[39] 王宗爽, 武婷, 段小丽, 等. 环境健康风险评价中我国居民呼吸速率暴露参数研究[J]. 环境科学研究, 2009, 22(10): 1171-1175.
Wang Z S, Wu T, Duan X L, et al. Research on inhalation rate exposure factors of Chinese residents in environmental health risk assessment[J]. Research of Environmental Sciences, 2009, 22(10): 1171-1175.
[40] 张文超, 吕森林, 刘丁彧, 等. 宣威街道尘中重金属的分布特征及其健康风险评估[J]. 环境科学, 2015, 36(5): 1810-1817.
Zhang W C, Lv S L, Liu D Y, et al. Distribution characteristics of heavy metals in the street dusts in Xuanwei and their health risk assessment[J]. Environmental Science, 2015, 36(5): 1810-1817.
[41] 马杰, 佘泽蕾, 王胜蓝, 等. 基于蒙特卡罗模拟的煤矸山周边农用地土壤重金属健康风险评估[J]. 环境科学, 2023, 44(10): 5666-5678.
Ma J, She Z L, Wang S L, et al. Health risk assessment of heavy metals in agricultural soils around the gangue heap of coal mine based on Monte Carlo simulation[J]. Environmental Science, 2023, 44(10): 5666-5678.
[42] 马建华, 董运武, 陈彦芳. 开封市周边地区地表灰尘重金属背景值研究[J]. 环境科学学报, 2020, 40(5): 1798-1806.
Ma J H, Dong Y W, Chen Y F. Background values of heavy metals in surface dusts in the vicinity of Kaifeng, Henan province[J]. Acta Scientiae Circumstantiae, 2020, 40(5): 1798-1806.
[43] 夏增禄, 李森照, 李廷芳, 等. 土壤元素背景值及其研究方法[M]. 北京: 气象出版社, 1987.
[44] GB 15618-2018, 土壤环境质量农用地土壤污染风险管控标准(试行)[S].
[45] 卢小慧, 余方中, 范一鸣, 等. 三门峡某铅厂遗留场地土壤重金属空间分布特征及来源解析[J]. 环境科学, 2023, 44(3): 1646-1656.
Lu X H, Yu F Z, Fan Y M, et al. Spatial distribution characteristics and source analysis of heavy metals in soil of a lead plant in Sanmenxia[J]. Environmental Science, 2023, 44(3): 1646-1656.
[46] Zhang W H, Yan Y, Yu R L, et al. The sources-specific health risk assessment combined with APCS/MLR model for heavy metals in tea garden soils from south Fujian province, China[J]. CATENA, 2021, 203. DOI:10.1016/j.Catena.2021.105306
[47] 李军, 李旭, 李开明, 等. 黄河流域农田土壤重金属污染特征及其优先控制源分析[J]. 环境科学, 2024, 45(3): 1724-1738.
Li J, Li X, Li K M, et al. Characteristics and identification priority source of heavy metals pollution in farmland soils in the Yellow River Basin[J]. Environmental Science, 2024, 45(3): 1724-1738.
[48] 张扣扣, 贺婧, 钟艳霞, 等. 基于GIS对宁夏某铜银矿区周边土壤重金属来源解析[J]. 环境科学, 2022, 43(11): 5192-5204.
Zhang K K, He J, Zhong Y X, et al. Identification of soil heavy metal sources around a copper-silver mining area in Ningxia based on GIS[J]. Environmental Science, 2022, 43(11): 5192-5204.
[49] 陈明, 王琳玲, 曹柳, 等. 基于PMF模型的某铅锌冶炼城市降尘重金属污染评价及来源解析[J]. 环境科学, 2023, 44(6): 3450-3462.
Chen M, Wang L L, Cao L, et al. Pollution assessment and source analysis of heavy metals in atmospheric deposition in a lead-zinc smelting city based on PMF model[J]. Environmental Science, 2023, 44(6): 3450-3462.
[50] 王显钦, 费学海, 杨员, 等. 贵阳市花溪城区PM2.5中重金属元素的污染特征、来源及健康风险评价[J]. 环境科学学报, 2023, 43(6): 110-118.
Wang X Q, Fei X H, Yang Y, et al. Pollution characteristics, source apportionment and health risk assessment of heavy metal elements in PM2.5 collected in Huaxi urban areas, Guiyang[J]. Acta Scientiae Circumstantiae, 2023, 43(6): 110-118.
[51] Wang F F, Guan Q Y, Tian J, et al. Contamination characteristics, source apportionment, and health risk assessment of heavy metals in agricultural soil in the Hexi Corridor[J]. CATENA, 2020, 191. DOI:10.1016/j.catena.2020.104573
[52] Heidari M, Darijani T, Alipour V. Heavy metal pollution of road dust in a city and its highly polluted suburb; quantitative source apportionment and source-specific ecological and health risk assessment[J]. Chemosphere, 2021, 273. DOI:10.1016/j.chemosphere.2021.129656
[53] Yang S Y, He M J, Zhi Y Y, et al. An integrated analysis on source-exposure risk of heavy metals in agricultural soils near intense electronic waste recycling activities[J]. Environment International, 2019, 133. DOI:10.1016/j.envint.2019.105239
[54] Huang C C, Cai L M, Xu Y H, et al. A comprehensive approach to quantify the source identification and human health risk assessment of toxic elements in park dust[J]. Environmental Geochemistry and Health, 2023, 45(8): 5813-5827.