宋清泉，徐夕博，吳泉源，楊 華，王翠秀，顧振飛，徐明雪

(1. 山東省國(guó)土空間數(shù)據(jù)和遙感技術(shù)研究院，中國(guó) 濟(jì)南 250002；2. 山東師范大學(xué)地理與環(huán)境學(xué)院，中國(guó) 濟(jì)南 250358；3. 北京師范大學(xué)地理科學(xué)學(xué)部，中國(guó) 北京 100875；4. 山東揚(yáng)華地理信息有限公司，中國(guó) 濟(jì)南 250102)

重金屬通常指密度大于5 g·cm-3的金屬或類金屬，此類污染物因隱蔽性強(qiáng)、毒性高、不易降解等特點(diǎn)嚴(yán)重危及生態(tài)系統(tǒng)的安全[1]。土壤中的重金屬在累積過(guò)程中會(huì)導(dǎo)致土壤養(yǎng)分流失，影響作物的品質(zhì)和產(chǎn)量，進(jìn)而降低土壤和糧食作物的質(zhì)量。另一方面，人體在經(jīng)皮膚接觸和吸入等途徑進(jìn)入循環(huán)系統(tǒng)后，對(duì)人體健康構(gòu)成威脅。例如，人體內(nèi)過(guò)多的鉛積累會(huì)損害人體神經(jīng)系統(tǒng)，破壞骨骼造血系統(tǒng)而導(dǎo)致貧血[2]。因此，為維持生態(tài)系統(tǒng)穩(wěn)定和保護(hù)人類身體健康，科學(xué)評(píng)估土壤重金屬的污染狀態(tài)、來(lái)源分配及潛在生態(tài)環(huán)境風(fēng)險(xiǎn)是極其必要的。

土壤重金屬是地殼的重要組成部分，直接去除土壤中的重金屬元素非常困難，費(fèi)時(shí)費(fèi)力，而有針對(duì)性地預(yù)防是處理土壤重金屬富集問(wèn)題的最佳策略。李春芳等[3]將表層土壤重金屬含量與深層未受人類干擾背景值進(jìn)行對(duì)比，分析并評(píng)估重金屬富集狀態(tài)及特征。SUN等[4]和徐夕博等[5]采用多元統(tǒng)計(jì)技術(shù)(相關(guān)分析、因子分析和聚類分析)挖掘重金屬相關(guān)變量信息，根據(jù)因子的貢獻(xiàn)度分析土壤重金屬的來(lái)源。SONG等[6]利用地理信息系統(tǒng)(GIS)繪制重金屬空間分布熱點(diǎn)圖，建立空間分布與工業(yè)生產(chǎn)布局之間的聯(lián)系，發(fā)掘潛在的污染源，幫助判別人類活動(dòng)對(duì)土壤重金屬分布的干擾強(qiáng)度。陳志凡等國(guó)內(nèi)外學(xué)者[7-10]采用環(huán)境風(fēng)險(xiǎn)評(píng)估模型探究了研究區(qū)內(nèi)土壤重金屬的潛在生態(tài)環(huán)境風(fēng)險(xiǎn)，量化了土壤重金屬的來(lái)源，并驗(yàn)證了PMF受體模型用于污染源解析的可靠性。上述研究解決了重金屬來(lái)源的識(shí)別及量化問(wèn)題，但由于不同重金屬生態(tài)毒性各異，貢獻(xiàn)度最大的污染源對(duì)應(yīng)環(huán)境風(fēng)險(xiǎn)未必最高，因此，對(duì)不同來(lái)源的土壤重金屬進(jìn)行潛在環(huán)境風(fēng)險(xiǎn)差異性及貢獻(xiàn)比率分析有利于確定污染源治理的優(yōu)先程度。

圣井鎮(zhèn)地處山東省工業(yè)強(qiáng)市濟(jì)南市區(qū)與章丘市交界處，鐵路公路縱橫交錯(cuò)。截至2015年，圣井鎮(zhèn)同時(shí)擁有4處煤礦和多家建材化工廠，工業(yè)發(fā)展與環(huán)境保護(hù)的矛盾日益突出，以圣井鎮(zhèn)為研究區(qū)，科學(xué)精準(zhǔn)地評(píng)估圣井鎮(zhèn)重金屬環(huán)境風(fēng)險(xiǎn)狀態(tài)十分必要。

1 材料與方法

1.1 研究區(qū)概況

圣井鎮(zhèn)地處濟(jì)南市城區(qū)東側(cè)，是章丘市連接濟(jì)南市區(qū)的重要通道，地理坐標(biāo)介于36°38′32″N～36°43′28″N，117°21′15″E～117°27′55″E(圖1)，占地面積約為52.61 km2。圣井鎮(zhèn)為溫帶大陸性氣候，四季分明，雨熱同期，地形分布主要為平原，最高峰為危山，海拔約為205 m。土壤類型以酸性坡積物母質(zhì)發(fā)育形成的褐土為主，持水保肥能力強(qiáng)，農(nóng)作物主要種植小麥和玉米[11]。圣井鎮(zhèn)交通便利，廠礦企業(yè)多分布在道路兩側(cè)。此外圣井鎮(zhèn)煤炭、石灰石和黏土等礦物儲(chǔ)量豐富，是章丘市重要的產(chǎn)煤基地。

圖1 研究區(qū)及采樣點(diǎn)示意圖

1.2 樣品采集與測(cè)試

在實(shí)驗(yàn)室內(nèi)ArcGIS 10.2數(shù)字電腦底圖的基礎(chǔ)上，結(jié)合土地利用類型圖、地質(zhì)構(gòu)造圖和工業(yè)分布圖，充分考慮道路通達(dá)性后完成采樣點(diǎn)的預(yù)設(shè)。在2018年7月共采集到25處表層土壤(0～20 cm)樣本，在采樣過(guò)程中利用多點(diǎn)混合法收集約1 kg土壤樣品，裝入聚乙烯密封袋中，送往實(shí)驗(yàn)室進(jìn)行化學(xué)元素測(cè)試，并將實(shí)際采樣的位置用手持GPS進(jìn)行記錄，具體如圖1所示。在實(shí)驗(yàn)室中，去除土塊雜草等自然雜物后，土壤樣品在室溫條件下(25 ℃)自然風(fēng)干，過(guò)0.25 mm篩后用瑪瑙研缽研磨細(xì)化，最終通過(guò)0.074 mm孔徑篩。在化學(xué)檢測(cè)階段，首先用HNO3-HCl-HF-HClO4進(jìn)行消煮，為接下來(lái)As，Cd，Co，Cr，Cu，Hg，Mn，Ni，Pb和Zn的測(cè)定做準(zhǔn)備。Mn，Cd和Co采用石墨爐原子吸收分光光度法(GE-AAS)測(cè)定，原子熒光光譜法(HG-AFS)用于As和Hg的測(cè)定，其余元素均采用等離子體原子發(fā)射光譜法(ICP-OES)進(jìn)行測(cè)定。在重金屬元素測(cè)試過(guò)程中，分析測(cè)試試劑選用優(yōu)級(jí)純，測(cè)試質(zhì)量的精確性采用國(guó)家一級(jí)地球化學(xué)標(biāo)準(zhǔn)物質(zhì)(GBW系列)進(jìn)行控制，10種重金屬的回收率保持在±10%，并對(duì)所有測(cè)試樣品進(jìn)行隨機(jī)抽樣檢查和異常點(diǎn)檢查，所得結(jié)果經(jīng)檢查達(dá)到優(yōu)良水平[12]。

1.3 數(shù)據(jù)處理

1.3.1 因子分析法 實(shí)驗(yàn)室中獲得的As等10種重金屬元素含量在SPSS 20.0軟件中對(duì)平均值、范圍、標(biāo)準(zhǔn)差、峰度和偏度描述性統(tǒng)計(jì)特征進(jìn)行計(jì)算，同時(shí)運(yùn)用R語(yǔ)言的psych包中對(duì)上述重金屬元素含量進(jìn)行因子分析。

1.3.2 受體模型 受體模型是一種基于多元統(tǒng)計(jì)技術(shù)的定量化污染源因子的分析技術(shù)，正交矩陣分解模型(Positive Matrix Factorization，PMF)是由美國(guó)環(huán)境保護(hù)署(EPA)推薦使用的一種可靠的受體模型，PMF模型首先將元素濃度矩陣分解為因子貢獻(xiàn)和因子殘差矩陣，用基于重金屬來(lái)源的先驗(yàn)知識(shí)確定不同因子的貢獻(xiàn)率，公式表示如下：

(1)

式中，gij表示第i個(gè)土壤樣品在第j個(gè)源的污染貢獻(xiàn)量，fkj代表第k個(gè)污染物在第j個(gè)污染源的貢獻(xiàn)，uij為濃度實(shí)測(cè)值的不確定性，eij表示模型的不確定性。此外，PMF模型構(gòu)建時(shí)，還需考慮樣品元素屬性的不確定性(Unc)，具體計(jì)算如下：

(2)

(3)

式中，cij表示重金屬元素實(shí)測(cè)濃度，θ和MDL分別為標(biāo)準(zhǔn)差和重金屬濃度的檢出限。計(jì)算過(guò)程中，需多次計(jì)算選取最優(yōu)結(jié)果，其中信噪比(Signal to noise,S/N)是評(píng)價(jià)輸入數(shù)據(jù)質(zhì)量的重要參數(shù)。S/N≥2表示數(shù)據(jù)變異程度高，質(zhì)量較好；反之說(shuō)明數(shù)據(jù)變異程度低，無(wú)法得到準(zhǔn)確結(jié)果[13]。

1.3.3 環(huán)境風(fēng)險(xiǎn)評(píng)價(jià) 環(huán)境風(fēng)險(xiǎn)評(píng)價(jià)主要基于潛在生態(tài)風(fēng)險(xiǎn)指數(shù)[14]綜合分析計(jì)算得來(lái)，環(huán)境風(fēng)險(xiǎn)評(píng)價(jià)指數(shù)構(gòu)建過(guò)程中須綜合考慮沉積轉(zhuǎn)換規(guī)律、污染物毒性及環(huán)境交互特性，計(jì)算過(guò)程如下：

fi=ci/Bi，

(4)

Ei=Ti·fi，

(5)

(6)

式中，fi為重金屬元素含量實(shí)測(cè)值與背景值的比值；ci為重金屬含量實(shí)測(cè)值；Bi為元素含量背景值；Ei為單種重金屬元素的環(huán)境風(fēng)險(xiǎn)指數(shù)；Ti表示單種重金屬元素的毒性響應(yīng)參數(shù)，劉碩等[15]和ZHANG等[16]研究表明，Hg，Cd，As，Co，Cu，Ni，Pb，Cr，Mn 和 Zn的Ti值分別為40，30，10，5，5，5，5，2，1和1；RI為多種重金屬綜合環(huán)境風(fēng)險(xiǎn)強(qiáng)度。

2 結(jié)果與討論

2.1 土壤重金屬含量特征分析

表1為研究區(qū)表層土壤重金屬含量的描述性統(tǒng)計(jì)特征?？砂l(fā)現(xiàn)As，Cd，Co，Cr，Cu，Hg，Mn，Ni，Pb和Zn 10種重金屬含量分別超出背景值[17]60.48%，54.55%，17.73%，27.14%，39.39%，136.67%，15.87%，42.09%，9.02%和27.13%，但含量均未超過(guò)農(nóng)用地土壤污染風(fēng)險(xiǎn)篩選值(GB 15618—2018)[18]，處于風(fēng)險(xiǎn)管控值以內(nèi)。其中Cd元素的極大值已接近篩選值，表明在局部區(qū)域土壤受到嚴(yán)重的人為擾動(dòng)或成土本底影響。10種重金屬元素的中值均超過(guò)背景值，表明研究區(qū)內(nèi)大多數(shù)土壤中的重金屬元素存在富集現(xiàn)象，除Mn和Zn以外，其余元素都出現(xiàn)了正偏現(xiàn)象，表明這些元素受到人為干擾偏離了原有的正態(tài)分布。統(tǒng)計(jì)數(shù)據(jù)的波動(dòng)強(qiáng)度通常采用變異系數(shù)進(jìn)行表征，其中Hg的變異系數(shù)(CV)為0.82，已經(jīng)達(dá)到高度變異狀態(tài)(CV>0.36)[19]，此種元素含量分布在外部活動(dòng)影響下處于高度波動(dòng)狀態(tài)，這也與峰度值的特征分析相一致。研究區(qū)土壤的pH值為7.48，表現(xiàn)為弱堿性，一定程度上會(huì)抑制重金屬的活性，有助于降低土壤重金屬的環(huán)境風(fēng)險(xiǎn)。

2.2 主成分分析

因子分析基于主分量提取技術(shù)，可以將復(fù)雜信息變量轉(zhuǎn)換到平面空間上的線性不相關(guān)變量，分析復(fù)雜變量在各個(gè)分量上的載荷，識(shí)別其在綜合變量上的信息量，有助于判斷土壤重金屬的來(lái)源[20]。表2為旋轉(zhuǎn)后的各主要分量載荷及貢獻(xiàn)率，提取出3個(gè)主分量且總方差貢獻(xiàn)率達(dá)到72.75%，能夠解釋復(fù)雜重金屬變量的大多數(shù)信息，符合預(yù)期要求?？梢钥闯觯谝恢鞣至?PC1)由As，Co，Cr，Cu，Mn和Ni組成，第二主分量(PC2)由Cd，Pb和Zn組成，第三主分量(PC3)為Hg元素。

表2 旋轉(zhuǎn)后分量載荷及貢獻(xiàn)率表

As，Co，Cr，Cu，Mn和Ni在PC1上的載荷均大于0.7，PC1可以反映上述元素的絕大部分信息。圣井鎮(zhèn)的成土母巖類型以石灰?guī)r為主，以此母巖為基礎(chǔ)的風(fēng)化土壤中As，Co，Cu和Mn的含量會(huì)高于其他類型土壤[21]。Mn是地殼組成的重要元素之一，自然土壤中具有較高的豐度，Co和Mn常伴生在同種礦物中，以自然來(lái)源為主。MIC等[22]的研究得出Co和Mn在同一主成分中，歸為成土母質(zhì)來(lái)源。LIANG等[23]在對(duì)漣源市礦區(qū)土壤的研究中發(fā)現(xiàn)Cr和Ni主要來(lái)源于成土母質(zhì)，空間變異性也比較低，故可將As，Co，Cr，Cu，Mn和Ni歸為自然來(lái)源。

Cd和Pb的載荷分別為0.742和0.787，Zn的載荷為0.565，這3種元素的信息主要集中在PC2。Cd，Pb和Zn具有一定的同源特性，是工業(yè)三廢排放的典型元素，Cd和Zn位屬同族，具有類似的化學(xué)特征，在石灰?guī)r風(fēng)化發(fā)育的土壤中含量更高。圣井鎮(zhèn)交通便利，交通污染物散落至土壤中導(dǎo)致Pb和Zn含量較高。已有研究表明[24,25]，交通運(yùn)輸活動(dòng)中的尾氣排放和輪胎磨損也會(huì)加劇土壤Pb和Zn的富集。YANG等[26]的研究表明在密集道路兩旁的塵土中通常含有高濃度的Zn，繼而通過(guò)大氣沉降作用進(jìn)入土壤。此外，圣井鎮(zhèn)分布著大量高標(biāo)準(zhǔn)農(nóng)田，在農(nóng)業(yè)生產(chǎn)活動(dòng)中Cd作為農(nóng)藥的有效成分被廣泛使用，且在我國(guó)利用率比較低，僅有約30%被吸收利用，含Cd農(nóng)藥流失至土壤中易產(chǎn)生元素富集[27,28]。綜上，Cd，Pb和Zn在土壤中受到母質(zhì)和農(nóng)業(yè)生產(chǎn)多種因素共同影響，屬于混合來(lái)源。

Hg在PC3上載荷達(dá)到0.813，可以反映該主成分上的主要信息。本研究中Hg的變異系數(shù)最大，受到的外部干擾強(qiáng)度最高，元素含量變化范圍和空間變異程度最大。張倩等[29]的研究表明Hg主要產(chǎn)生在燃煤化學(xué)轉(zhuǎn)換過(guò)程中，通過(guò)粉煤灰大氣沉降作用降落至土壤產(chǎn)生富集，圣井鎮(zhèn)內(nèi)建有多處煤炭基地，煤炭生產(chǎn)和使用易使Hg在土壤中富集。蔣靖坤等[30]和劉巍等[31]也在研究中證實(shí)工業(yè)、電力和生活消費(fèi)中的燃煤通過(guò)化學(xué)轉(zhuǎn)化、遷移和沉降作用將Hg元素釋放至環(huán)境中，在土壤中積累。綜上，可將Hg判識(shí)為工業(yè)活動(dòng)影響下的人為來(lái)源。

2.3 PMF因子分析

土壤重金屬在PMF因子上的載荷是重金屬來(lái)源定量分配的重要依據(jù)[13]。模型在運(yùn)行20次后完成因子分配，在對(duì)Q值進(jìn)行綜合測(cè)度之后且S/N>2，元素預(yù)測(cè)值與實(shí)測(cè)值相關(guān)系數(shù)r2均大于0.5，擬合程度較好，PMF模型解析得到的各元素因子貢獻(xiàn)率如圖2所示。在As，Cd，Co，Cr，Cu，Mn，Ni，Pb和Zn中，因子1和因子2所占比例最大，因子2在Hg中所占比例最大，結(jié)合主成分分析結(jié)果，因子2應(yīng)為工業(yè)來(lái)源所占載荷。Cd，Pb和Zn在載荷的比例分配中，因子1所占權(quán)重最大，故因子1可歸為混合來(lái)源，而As，Co，Cr，Cu，Mn和Ni的因子3約等于因子1，結(jié)合主成分分析結(jié)果可將因子3判識(shí)為自然母質(zhì)來(lái)源。圖3為PMF模型得到的不同來(lái)源對(duì)土壤重金屬貢獻(xiàn)率，其中工業(yè)來(lái)源占比22.3%，母質(zhì)和混合來(lái)源分別占38.71%和38.99%。

圖2 土壤重金屬不同因子解析結(jié)果

圖3 不同來(lái)源對(duì)土壤重金屬貢獻(xiàn)度

PMF模型與因子分析法均是進(jìn)行土壤重金屬源解析的重要手段。一方面，因子分析法可以與PMF模型解析的結(jié)果進(jìn)行驗(yàn)證對(duì)比，提高源解析工作的準(zhǔn)確度。另一方面，因子分析法所確定的3個(gè)主成分，可以定性地得出不同土壤重金屬元素的主要來(lái)源和總的污染源貢獻(xiàn)度，但不能準(zhǔn)確解析出各個(gè)污染源的貢獻(xiàn)比率和各個(gè)元素對(duì)污染源的貢獻(xiàn)度，PMF模型可以通過(guò)調(diào)整因子數(shù)和Q誤差值，完成對(duì)污染源貢獻(xiàn)率的定量分析。但PMF模型在不確定度計(jì)算中會(huì)產(chǎn)生一定的誤差，影響因子分析結(jié)果，所以PMF模型與因子分析法的結(jié)合可以快速準(zhǔn)確地完成土壤重金屬的定量源解析。

2.4 環(huán)境風(fēng)險(xiǎn)定量解析評(píng)價(jià)

環(huán)境生態(tài)風(fēng)險(xiǎn)指數(shù)主要用于評(píng)估區(qū)域內(nèi)的綜合潛在生態(tài)環(huán)境風(fēng)險(xiǎn)[29]。Ei和RI分別用于評(píng)價(jià)單種和多種重金屬元素所帶來(lái)的潛在環(huán)境風(fēng)險(xiǎn)，環(huán)境生態(tài)風(fēng)險(xiǎn)等級(jí)的確定應(yīng)根據(jù)污染物的類別和毒性響應(yīng)系數(shù)進(jìn)行劃分，參照HAKANSON[14]和周旭等[32]的研究，結(jié)合本研究區(qū)的實(shí)際情況，風(fēng)險(xiǎn)等級(jí)劃分標(biāo)準(zhǔn)如表3所示。研究區(qū)內(nèi)各重金屬元素總環(huán)境生態(tài)指數(shù)值RI為186.78，整體上處于中度等級(jí)的環(huán)境生態(tài)風(fēng)險(xiǎn)。圖4為各個(gè)重金屬的環(huán)境生態(tài)指數(shù)，As，Cd，Co，Cr，Cu，Hg，Mn，Ni，Pb和Zn環(huán)境生態(tài)風(fēng)險(xiǎn)指數(shù)均值分別為16.05，47.15，5.89，2.54，6.97，94.35，0.01，7.10，5.45和1.27，對(duì)應(yīng)元素由高到低依次為：Hg，Cd，As，Ni，Cu，Co，Pb，Cr，Zn，Mn。研究區(qū)內(nèi)Hg的潛在生態(tài)風(fēng)險(xiǎn)等級(jí)最高，16處(64%)采樣點(diǎn)處于中度生態(tài)風(fēng)險(xiǎn)等級(jí)，6處(24%)采樣點(diǎn)處于較強(qiáng)生態(tài)風(fēng)險(xiǎn)，3處(12%)采樣點(diǎn)表現(xiàn)出了強(qiáng)等級(jí)的潛在生態(tài)風(fēng)險(xiǎn)，總體達(dá)到較強(qiáng)生態(tài)風(fēng)險(xiǎn)狀態(tài)。研究區(qū)內(nèi)布設(shè)的20處(80%)表層土壤樣點(diǎn)中Cd元素處于中度生態(tài)風(fēng)險(xiǎn)等級(jí)，其余樣點(diǎn)中Cd未出現(xiàn)強(qiáng)等級(jí)的生態(tài)風(fēng)險(xiǎn)，總體上Cd為中度生態(tài)風(fēng)險(xiǎn)；其余元素均處于輕微生態(tài)風(fēng)險(xiǎn)等級(jí)。

圖4 各元素潛在生態(tài)環(huán)境風(fēng)險(xiǎn)結(jié)果

表3 潛在生態(tài)風(fēng)險(xiǎn)劃分標(biāo)準(zhǔn)

基于PMF模型的不同重金屬來(lái)源對(duì)各個(gè)采樣點(diǎn)環(huán)境風(fēng)險(xiǎn)貢獻(xiàn)率和對(duì)潛在生態(tài)環(huán)境風(fēng)險(xiǎn)貢獻(xiàn)率分別如圖5和圖6所示。總的來(lái)看，不同土壤樣點(diǎn)間重金屬工業(yè)源和混合源帶來(lái)的環(huán)境風(fēng)險(xiǎn)變異程度最大，采樣點(diǎn)14和20重金屬來(lái)源中混合源帶來(lái)的潛在環(huán)境風(fēng)險(xiǎn)最大，而母質(zhì)源所引發(fā)的環(huán)境風(fēng)險(xiǎn)不足4%，而在采樣點(diǎn)10，11和21中，工業(yè)來(lái)源的重金屬主導(dǎo)了該樣點(diǎn)的環(huán)境風(fēng)險(xiǎn)，這也與環(huán)境生態(tài)風(fēng)險(xiǎn)的綜合評(píng)價(jià)相驗(yàn)證。以土壤的風(fēng)化母質(zhì)為重金屬來(lái)源的樣點(diǎn)，其環(huán)境風(fēng)險(xiǎn)的整體變異程度小于工業(yè)源和混合來(lái)源。工業(yè)源的重金屬對(duì)土壤生態(tài)環(huán)境風(fēng)險(xiǎn)貢獻(xiàn)率最大，約為49.69%，其次為混合來(lái)源29.11%。土壤重金屬的最大來(lái)源成土母質(zhì)部分，所產(chǎn)生的潛在環(huán)境風(fēng)險(xiǎn)卻最小，約為21.20%。土壤重金屬的工業(yè)來(lái)源部分僅占22.3%，卻對(duì)整體生態(tài)環(huán)境風(fēng)險(xiǎn)貢獻(xiàn)占比49.69%，應(yīng)引起相關(guān)部門關(guān)注。

圖5 不同來(lái)源對(duì)各采樣點(diǎn)環(huán)境風(fēng)險(xiǎn)貢獻(xiàn)率

圖6 不同來(lái)源對(duì)環(huán)境風(fēng)險(xiǎn)貢獻(xiàn)率

3 結(jié)論

1)表層土壤中As，Cd，Co，Cr，Cu，Hg，Mn，Ni，Pb和Zn含量分別超出元素背景值60.48%，54.55%，17.73%，27.14%，39.39%，136.67%，15.87%，42.09%，9.02%和27.13%，但未超出農(nóng)用地土壤污染風(fēng)險(xiǎn)篩選值，土壤重金屬出現(xiàn)一定程度富集現(xiàn)象。

2)土壤重金屬的來(lái)源可分為3類：As，Co，Cr，Mn和Ni受成土母質(zhì)影響，屬于自然來(lái)源；Hg主要來(lái)源于工業(yè)排放，屬人為來(lái)源；Cd，Pb和Zn受到交通和母質(zhì)雙重因素控制，屬混合來(lái)源。其中工業(yè)來(lái)源占比22.3%，母質(zhì)和混合來(lái)源分別占38.71%和38.99%。

3)研究區(qū)整體上處于中度生態(tài)風(fēng)險(xiǎn)等級(jí)，工業(yè)來(lái)源的重金屬對(duì)土壤生態(tài)環(huán)境風(fēng)險(xiǎn)貢獻(xiàn)率最大，占比49.69%，其次為混合來(lái)源29.11%，成土母質(zhì)部分產(chǎn)生的潛在環(huán)境風(fēng)險(xiǎn)最小，約為21.20%。土壤重金屬最小的工業(yè)來(lái)源部分卻產(chǎn)生了較大的生態(tài)環(huán)境風(fēng)險(xiǎn)貢獻(xiàn)率，需引起相關(guān)部門關(guān)注。