于蕾　魯成秀　劉玉真　劉伏　成杰民

摘要：以山東省基本農(nóng)田土壤為研究對象，采用地積累指數(shù)法和潛在生態(tài)風(fēng)險指數(shù)法分析了土壤中Cu、Zn、Pb、Cd的污染程度。由潛在生態(tài)風(fēng)險指數(shù)可以看出山東省農(nóng)田土壤重金屬潛在生態(tài)風(fēng)險級別褐土為B級，潛在生態(tài)危害程度中等；潮土為A級，潛在生態(tài)危害程度輕微；棕壤為B級，潛在生態(tài)危害程度中等。就單個元素來看，Cd元素的潛在生態(tài)危害程度最高，各元素的潛在生態(tài)危害程度為Cd>Pb>Cu>Zn。由地累積指數(shù)法可以看出Pb的風(fēng)險級別最大；其次是Cu元素；Cd、Zn的污染程度相比較而言較輕，風(fēng)險級別低，各元素的污染程度和風(fēng)險級別為Pb>Cu>Cd>Zn。

關(guān)鍵詞：土壤；重金屬；污染評價；地積累指數(shù)法；潛在生態(tài)風(fēng)險指數(shù)法

中圖分類號：X53 文獻(xiàn)標(biāo)識碼：A 文章編號：0439-8114（2014）13-3010-04

Assessing the Potential Ecological Risks of Heavy Metals in Farmland Soils in Shandong Province

YU Lei，LU Cheng-xiu，LIU Yu-zhen，LIU Fu，CHENG Jie-min

（College of Population Resources and Environment， Shandong Normal University， Jinan 250014， China）

Abstract： Using the basic farmland of Shandong Province served as object， the potential ecological risks of Cu and Zn， Pb， Cd in the soil was evaluated by using the index of potential ecological risk and the index of geoaccumulation.Results based on the index of potential ecological risk showed that the potential ecological risk of moisture soil was at B level， indicating that the ecological damage was moderate. The potential ecological risk of brunisolic soil were at A level， indicating that the ecological damage was not serious； while the potential ecological risk of brown earth was at B level， indicating that the ecological damage was moderate. The potential ecological risk degrees of heavy metals were ranked in order of Cd>Pb>Cu>Zn. Results based on index of geoaccumulation showed that the potential ecological risks of heavy metals were ranked in order of Pb>Cu>Zn>Cd.

Key words： soil； heavy metals； pollution assessment； index of geoaccumulation； index of potential ecologicalrisk

近年來，我國農(nóng)業(yè)生產(chǎn)在快速發(fā)展的同時，農(nóng)業(yè)生態(tài)環(huán)境也遭受著嚴(yán)重的污染和破壞[1]。調(diào)查表明，我國污灌區(qū)被重金屬所污染的土地面積已達(dá)污灌區(qū)面積的64.8%，所以農(nóng)村生態(tài)被稱為“中國環(huán)保的短板”[2]，分析土壤重金屬元素含量對研究人為活動對土壤質(zhì)量的影響以及合理開發(fā)和利用土地資源具有重要意義[3]。根據(jù)農(nóng)業(yè)部對全國污灌區(qū)進(jìn)行的調(diào)查表明，在我國大約140萬hm2的污水灌溉區(qū)中，已經(jīng)遭受重金屬污染的土地面積占到污水灌區(qū)面積的64.8%，具體為輕度污染的占46.7%，中度污染的占9.7%，而嚴(yán)重污染的占8.4%[4]。由農(nóng)田土壤及作物的重金屬污染所引起的潛在健康風(fēng)險引起了國內(nèi)外學(xué)者的廣泛關(guān)注[5-7]。對重金屬進(jìn)行生態(tài)風(fēng)險評價的方法很多，其中常用的有地積累指數(shù)法及潛在生態(tài)風(fēng)險指數(shù)法等。地積累指數(shù)法主要對沉積物或土壤中的重金屬污染程度及其分級情況進(jìn)行定量評價[8，9]。潛在生態(tài)風(fēng)險指數(shù)法可以將生物毒性、生態(tài)危害與污染物濃度有機結(jié)合起來，從而綜合反映重金屬對生態(tài)環(huán)境的影響潛力[10]。本研究以山東省典型農(nóng)田土壤為對象，于2009-2010年對山東省90%以上的棕壤、褐土、潮土等主要土壤類型進(jìn)行調(diào)查，并在此基礎(chǔ)上采用地累積指數(shù)法和潛在生態(tài)風(fēng)險指數(shù)法對山東省典型農(nóng)田土壤重金屬的生態(tài)風(fēng)險進(jìn)行評價，從而為采用何種方法對污染土壤進(jìn)行科學(xué)管理、修復(fù)、治理并防止污染進(jìn)一步發(fā)展提供科學(xué)依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

山東省地處黃河下游，位于東徑114°36′-122°43′，北緯34°22′-38°33′之間，土地總面積15.7萬km2，其中耕地面積為733.5萬hm2。

山東省主要土壤類型有棕壤、褐土、潮土和鹽土等土壤類型。其中褐土占全省土壤總面積的18.16%、潮土占41.10%、棕壤占30.66%，總計約90%。棕壤、褐土、潮土為山東省主要土壤類型[11]，同時也是本研究農(nóng)田土壤的3種類型。

1.2 樣品采集與測定

按照土壤類型和作物種植品種分布及土壤肥力高、中、低分別采樣，采用全球定位系統(tǒng)進(jìn)行全省范圍內(nèi)的精確布設(shè)代表性采樣點60個（其中褐土25個，潮土16個，棕壤19個）（圖1），采集農(nóng)田耕層土壤（0～20 cm），風(fēng)干，磨細(xì)，過篩，備用。土壤中Cu、Zn的測定采用火焰原子吸收分光光度法[12]；土壤中Pb、Cd的測定采用KI-MIBK萃取火焰原子吸收分光光度法[13]。

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

1.3.1 地累積指數(shù)法 地積累指數(shù)（Index of geo-accumulation）又稱Mull指數(shù)，地積累指數(shù)法考慮了元素相對于自然本底值的富集性，主要側(cè)重于從自然角度對土壤進(jìn)行評價[10] 。

計算公式如下：Igeo=log2Cn/（K·Bn）

式中，Cn為實測重金屬元素的含量，mg/kg；Bn為當(dāng)?shù)爻练e物中重金屬元素含量的地球化學(xué)背景值，mg/kg；K為考慮到各地成巖作用不同引起背景值波動所設(shè)定的常數(shù)，K=1.5。地累積指數(shù)法分級標(biāo)準(zhǔn)見表1。

1.3.2 潛在生態(tài)風(fēng)險指數(shù)法 潛在生態(tài)風(fēng)險指數(shù) （The potential ecologicalrisk index） 法則考慮了各重金屬元素的毒性，更側(cè)重于從生物和人的角度對土壤進(jìn)行評價[14]。計算公式如下：

式中， RI為多種重金屬元素的潛在生態(tài)風(fēng)險指數(shù)； Eir為第i種重金屬元素的潛在生態(tài)風(fēng)險指數(shù)； Cif為第i種重金屬元素的污染系數(shù)；Ci為所測樣品中第i種重金屬元素含量的實測值，mg/kg；Cin為第i種重金屬元素含量的背景值，mg/kg；Tir為第i種重金屬元素的毒性響應(yīng)參數(shù)[14]。潛在生態(tài)風(fēng)險指數(shù)法分級標(biāo)準(zhǔn)見表2。

2 結(jié)果與分析

2.1 地累積指數(shù)法評價結(jié)果

山東省農(nóng)田土壤60個采樣點的重金屬污染地累積指數(shù)不同風(fēng)險級別的頻數(shù)及比例如表3所示。

根據(jù)地累積指數(shù)法分級標(biāo)準(zhǔn)可知，山東省典型土壤中Zn、Cu、Pb、Cd等元素多數(shù)樣點在無污染至中等—強污染范圍內(nèi)。其中，Pb的污染最重，其中污染程度達(dá)到強—極嚴(yán)重污染和強污染的采樣點各有1個，風(fēng)險級別分別為5級和4級。另外有21.7%的采樣點達(dá)中等—強污染的程度，值得重視。其次是Cu元素，有5.0%的采樣點達(dá)中等—強污染的污染程度，3級風(fēng)險，28.3%的采樣點達(dá)中等污染程度。Cd、Zn的污染程度相對較輕，分別有40.0%和26.7%的采樣點土壤達(dá)到中等污染程度，其余為無污染或輕度—中等污染程度，風(fēng)險級別較低。

就不同的土壤類型來看（表4），褐土中Cu、Zn風(fēng)險級別為1級，Cd、Pb為2級，各元素的風(fēng)險程度依次為Pb>Cd>Cu>Zn；潮土中Cu、Cd、Zn為1級風(fēng)險，Pb為2級，各元素的風(fēng)險程度依次為Pb>Cu>Cd>Zn；棕壤中Cd、Zn為0級風(fēng)險，Pb為1級風(fēng)險，Cu的風(fēng)險級別為2級，各元素的風(fēng)險程度依次為Cu>Pb>Zn>Cd。

2.2 潛在生態(tài)風(fēng)險指數(shù)法評價結(jié)果

研究區(qū)農(nóng)田表層土壤中各元素的單項潛在生態(tài)風(fēng)險指數(shù)和綜合潛在生態(tài)風(fēng)險指數(shù)（表5）顯示，所有采樣點的Cu和Zn元素的潛生態(tài)風(fēng)險指數(shù)均小于40，風(fēng)險級別為A，潛在生態(tài)危害程度輕微；對于Pb元素，占總數(shù)3.30%的采樣點其潛在生態(tài)風(fēng)險指數(shù)大于80但小于160，風(fēng)險級別為C，生態(tài)危害程度強，占總數(shù)10%的采樣點，其潛在生態(tài)風(fēng)險指數(shù)大于40小于或等于80，潛在生態(tài)風(fēng)險級別為B，潛在生態(tài)危害程度中等，其余監(jiān)測點的潛在生態(tài)風(fēng)險指數(shù)均小于或等于40，屬A級風(fēng)險級別，對生態(tài)有輕微危害，全省所有監(jiān)測點平均潛在生態(tài)風(fēng)險級別為A級；對于Cd元素，占總數(shù)3.30%的采樣點，其潛在生態(tài)風(fēng)險指數(shù)大于160，風(fēng)險級別D級，潛在生態(tài)危害程度極強，A、B、C三個級別采樣點所占比例分別為20.00%、36.70%和40.00%，全省平均潛在生態(tài)風(fēng)險指數(shù)為82.78，大于80，屬于C級，對生態(tài)具有強污染。由此可看出，Cd污染較為嚴(yán)重，各元素的潛在生態(tài)危害程度為Cd>Pb>Cu>Zn。

綜合多元素，從綜合潛在生態(tài)風(fēng)險指數(shù)（表6）來看，山東省基本農(nóng)田土壤中褐土和棕壤潛在生態(tài)風(fēng)險級別為B級，潛在生態(tài)危害程度中等，潮土的潛在生態(tài)風(fēng)險級別為A級，潛在生態(tài)危害程度輕微。

3 結(jié)論與討論

從地累積指數(shù)可以看出，Pb的污染最重；其次是Cu元素，有5.0%的采樣點達(dá)中等—強污染的程度，3級風(fēng)險，28.3%的采樣點達(dá)中等污染程度。Cd、Zn的污染程度相對較輕，分別有40%和26.7%的采樣點土壤達(dá)到中等污染程度，其余為無污染或輕度—中等污染，風(fēng)險等級較低。

根據(jù)地累積指數(shù)法，就不同的土壤類型來看，褐土中各元素的風(fēng)險程度為Pb>Cd>Cu>Zn；潮土中各元素的風(fēng)險程度為Pb>Cu>Cd>Zn；棕壤中各元素的風(fēng)險程度為Cu>Pb>Zn>Cd。

研究區(qū)農(nóng)田表層土壤中各元素的單項潛在生態(tài)風(fēng)險指數(shù)和綜合潛在生態(tài)風(fēng)險指數(shù)顯示，所有采樣點的Cu和Zn元素的潛在生態(tài)風(fēng)險指數(shù)均小于40，潛在風(fēng)險級別為A級，潛在生態(tài)危害程度輕微，其中Zn元素的潛在生態(tài)風(fēng)險指數(shù)范圍為0.23～4.70，Cu元素的潛在生態(tài)風(fēng)險指數(shù)范圍為2.31～36.84；對于Pb元素，潛在生態(tài)風(fēng)險指數(shù)范圍為3.03～136.23；對于Cd元素，潛在生態(tài)風(fēng)險指數(shù)范圍為9.88～173.43。由此可看出，Cd元素的潛在生態(tài)危害最大，各元素的潛在生態(tài)危害程度為Cd>Pb>Cu>Zn。

兩種方法都得出Zn元素的污染程度最低，但是對于其他3種元素的結(jié)果均不相同，這是各方法的要求不同造成的，具體采用何種方法應(yīng)根據(jù)研究目的而定。根據(jù)潛在生態(tài)風(fēng)險指數(shù)法的評價結(jié)果，Cd元素的潛在生態(tài)危害最大，但是根據(jù)地累積指數(shù)法的評價結(jié)果，只有40%的采樣點土壤Cd達(dá)到中等污染程度，其余屬無污染或輕度—中等污染程度，風(fēng)險等級較低。在成杰民等[15]對Cu、Cd、Pb、Zn的積累速率的計算中發(fā)現(xiàn)，4種重金屬元素中雖然Cd的積累速率非常低，但由于其本身原始含量就較低，其年變化速率卻高于Cu、Zn，僅次于Pb，這從另一方面說明了Cd存在較大潛在風(fēng)險。賈琳等[1]在對山東禹城農(nóng)田土壤的研究中同樣發(fā)現(xiàn)其土壤中Hg和Cd潛在生態(tài)危害指數(shù)較大，存在較大的潛在生態(tài)風(fēng)險。因為禹城為典型施肥區(qū)，其畜禽養(yǎng)殖和污灌以及城市化進(jìn)程是造成土壤中Cd含量超過原有背景值的主要因素。

農(nóng)田土壤的質(zhì)量與人類的生產(chǎn)活動密切相關(guān)，因此對于農(nóng)田土壤重金屬污染的危害應(yīng)多從人類和生物的角度考慮，對毒性的研究要多加注意。潛在生態(tài)風(fēng)險指數(shù)法不僅可以反映在一定環(huán)境中的全部污染物的影響，并且通過潛在生態(tài)危害指數(shù)的計算指出了其中應(yīng)該特別注意的物質(zhì)，所以對于污染的控制非常重要[16]。由此來看，采用潛在生態(tài)風(fēng)險指數(shù)法對農(nóng)田土壤重金屬污染進(jìn)行評價更適合此次的研究目的。

參考文獻(xiàn)：

[1] 賈 琳，楊林生，歐陽竹，等.典型農(nóng)業(yè)區(qū)農(nóng)田土壤重金屬潛在生態(tài)風(fēng)險評價[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報，2009，28（11）：2270-2276.

[2] 農(nóng)村生態(tài)：中國環(huán)保的短板[EB/OL].http：//www.chinaenvironment.com/action/Topic/ti_an/ViewNews.aspx？id=1025，2007-03-14.

[3] 羅真富，譚德軍，謝洪斌，等.重慶長壽湖周邊地區(qū)土壤重金屬污染評價[J].湖北農(nóng)業(yè)科學(xué)，2012，51（1）：30-34.

[4] 陳志良，仇榮亮.重金屬污染土壤的修復(fù)技術(shù)[J].環(huán)境保護(hù)，2002，29（6）：21-23.

[5] BERNHARD A Z， CHE F I， MIKE J M， et al. Heavy metals in soils and crops in southeast Asia. l. Peninsular Malaysia[J].Environmental Geochemistry and Health，2004，26：343-357.

[6] MAPANDA F， MANGWAYANA E N， GILLER K E， et al. Uptake of heavy metals by vegetables irrigated using wastewater and the subsequent risks in Harare， Zimbabwe[J].Physics Chemistry of the Earth，2007，32：1399-1405.

[7] 叢 源，鄭 萍，陳岳龍，等.北京市農(nóng)田生態(tài)系統(tǒng)土壤重金屬元素的生態(tài)風(fēng)險評價[J].地質(zhì)通報，2008，27（5）：681-688.

[8] 盧 瑛，龔子同，張甘霖.南京城市土壤Pb的含量及其化學(xué)形態(tài)[J].環(huán)境科學(xué)學(xué)報，2002，22（2）：156-160.

[9] 賈振邦，周 華，趙智杰，等.應(yīng)用地積累指數(shù)法評價太子河沉積物中重金屬污染[J].北京大學(xué)學(xué)報（自然科學(xué)版），2000，36（4）：525-530.

[10] HAKANSON L. An ecological risk index for aquatic pollution control. A sediment logical approach[J]. Water Research，1980，14：975-1001.

[11] 山東省土壤肥料工作站.山東土壤[M].北京：中國農(nóng)業(yè)出版社，1994.

[12] GB/T17138-1997，土壤質(zhì)量銅、鋅的測定 火焰原子吸收分光光度法[S].

[13] GB/T 17140-1997，土壤質(zhì)量鉛、鎘的測定 KI-MIBK萃取火焰原子吸收分光光度法[S].

[14] 于云江，胡林凱，楊 彥，等.典型流域農(nóng)田土壤重金屬污染特征及生態(tài)風(fēng)險評價[J].環(huán)境科學(xué)研究，2010，23（12）：1523-1527.

[15] 成杰民，張麗娜.近20年山東省典型農(nóng)田土壤中重金屬含量變化特征及積累速率估計[J].土壤，2013，45（1）：99-104.

[16] 賈振邦，梁 濤，林健枝，等.香港河流重金屬污染及潛在生態(tài)危害研究[J].北京大學(xué)學(xué)報（自然科學(xué)版），1997，33（4）：485-492.

農(nóng)田土壤的質(zhì)量與人類的生產(chǎn)活動密切相關(guān)，因此對于農(nóng)田土壤重金屬污染的危害應(yīng)多從人類和生物的角度考慮，對毒性的研究要多加注意。潛在生態(tài)風(fēng)險指數(shù)法不僅可以反映在一定環(huán)境中的全部污染物的影響，并且通過潛在生態(tài)危害指數(shù)的計算指出了其中應(yīng)該特別注意的物質(zhì)，所以對于污染的控制非常重要[16]。由此來看，采用潛在生態(tài)風(fēng)險指數(shù)法對農(nóng)田土壤重金屬污染進(jìn)行評價更適合此次的研究目的。

參考文獻(xiàn)：

[1] 賈 琳，楊林生，歐陽竹，等.典型農(nóng)業(yè)區(qū)農(nóng)田土壤重金屬潛在生態(tài)風(fēng)險評價[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報，2009，28（11）：2270-2276.

[2] 農(nóng)村生態(tài)：中國環(huán)保的短板[EB/OL].http：//www.chinaenvironment.com/action/Topic/ti_an/ViewNews.aspx？id=1025，2007-03-14.

[3] 羅真富，譚德軍，謝洪斌，等.重慶長壽湖周邊地區(qū)土壤重金屬污染評價[J].湖北農(nóng)業(yè)科學(xué)，2012，51（1）：30-34.

[4] 陳志良，仇榮亮.重金屬污染土壤的修復(fù)技術(shù)[J].環(huán)境保護(hù)，2002，29（6）：21-23.

[5] BERNHARD A Z， CHE F I， MIKE J M， et al. Heavy metals in soils and crops in southeast Asia. l. Peninsular Malaysia[J].Environmental Geochemistry and Health，2004，26：343-357.

[6] MAPANDA F， MANGWAYANA E N， GILLER K E， et al. Uptake of heavy metals by vegetables irrigated using wastewater and the subsequent risks in Harare， Zimbabwe[J].Physics Chemistry of the Earth，2007，32：1399-1405.

[7] 叢 源，鄭 萍，陳岳龍，等.北京市農(nóng)田生態(tài)系統(tǒng)土壤重金屬元素的生態(tài)風(fēng)險評價[J].地質(zhì)通報，2008，27（5）：681-688.

[8] 盧 瑛，龔子同，張甘霖.南京城市土壤Pb的含量及其化學(xué)形態(tài)[J].環(huán)境科學(xué)學(xué)報，2002，22（2）：156-160.

[9] 賈振邦，周 華，趙智杰，等.應(yīng)用地積累指數(shù)法評價太子河沉積物中重金屬污染[J].北京大學(xué)學(xué)報（自然科學(xué)版），2000，36（4）：525-530.

[10] HAKANSON L. An ecological risk index for aquatic pollution control. A sediment logical approach[J]. Water Research，1980，14：975-1001.

[11] 山東省土壤肥料工作站.山東土壤[M].北京：中國農(nóng)業(yè)出版社，1994.

[12] GB/T17138-1997，土壤質(zhì)量銅、鋅的測定 火焰原子吸收分光光度法[S].

[13] GB/T 17140-1997，土壤質(zhì)量鉛、鎘的測定 KI-MIBK萃取火焰原子吸收分光光度法[S].

[14] 于云江，胡林凱，楊 彥，等.典型流域農(nóng)田土壤重金屬污染特征及生態(tài)風(fēng)險評價[J].環(huán)境科學(xué)研究，2010，23（12）：1523-1527.

[15] 成杰民，張麗娜.近20年山東省典型農(nóng)田土壤中重金屬含量變化特征及積累速率估計[J].土壤，2013，45（1）：99-104.

[16] 賈振邦，梁 濤，林健枝，等.香港河流重金屬污染及潛在生態(tài)危害研究[J].北京大學(xué)學(xué)報（自然科學(xué)版），1997，33（4）：485-492.

農(nóng)田土壤的質(zhì)量與人類的生產(chǎn)活動密切相關(guān)，因此對于農(nóng)田土壤重金屬污染的危害應(yīng)多從人類和生物的角度考慮，對毒性的研究要多加注意。潛在生態(tài)風(fēng)險指數(shù)法不僅可以反映在一定環(huán)境中的全部污染物的影響，并且通過潛在生態(tài)危害指數(shù)的計算指出了其中應(yīng)該特別注意的物質(zhì)，所以對于污染的控制非常重要[16]。由此來看，采用潛在生態(tài)風(fēng)險指數(shù)法對農(nóng)田土壤重金屬污染進(jìn)行評價更適合此次的研究目的。

參考文獻(xiàn)：

[1] 賈 琳，楊林生，歐陽竹，等.典型農(nóng)業(yè)區(qū)農(nóng)田土壤重金屬潛在生態(tài)風(fēng)險評價[J].農(nóng)業(yè)環(huán)境科學(xué)學(xué)報，2009，28（11）：2270-2276.

[2] 農(nóng)村生態(tài)：中國環(huán)保的短板[EB/OL].http：//www.chinaenvironment.com/action/Topic/ti_an/ViewNews.aspx？id=1025，2007-03-14.

[3] 羅真富，譚德軍，謝洪斌，等.重慶長壽湖周邊地區(qū)土壤重金屬污染評價[J].湖北農(nóng)業(yè)科學(xué)，2012，51（1）：30-34.

[4] 陳志良，仇榮亮.重金屬污染土壤的修復(fù)技術(shù)[J].環(huán)境保護(hù)，2002，29（6）：21-23.

[5] BERNHARD A Z， CHE F I， MIKE J M， et al. Heavy metals in soils and crops in southeast Asia. l. Peninsular Malaysia[J].Environmental Geochemistry and Health，2004，26：343-357.

[6] MAPANDA F， MANGWAYANA E N， GILLER K E， et al. Uptake of heavy metals by vegetables irrigated using wastewater and the subsequent risks in Harare， Zimbabwe[J].Physics Chemistry of the Earth，2007，32：1399-1405.

[7] 叢 源，鄭 萍，陳岳龍，等.北京市農(nóng)田生態(tài)系統(tǒng)土壤重金屬元素的生態(tài)風(fēng)險評價[J].地質(zhì)通報，2008，27（5）：681-688.

[8] 盧 瑛，龔子同，張甘霖.南京城市土壤Pb的含量及其化學(xué)形態(tài)[J].環(huán)境科學(xué)學(xué)報，2002，22（2）：156-160.

[9] 賈振邦，周 華，趙智杰，等.應(yīng)用地積累指數(shù)法評價太子河沉積物中重金屬污染[J].北京大學(xué)學(xué)報（自然科學(xué)版），2000，36（4）：525-530.

[10] HAKANSON L. An ecological risk index for aquatic pollution control. A sediment logical approach[J]. Water Research，1980，14：975-1001.

[11] 山東省土壤肥料工作站.山東土壤[M].北京：中國農(nóng)業(yè)出版社，1994.

[12] GB/T17138-1997，土壤質(zhì)量銅、鋅的測定 火焰原子吸收分光光度法[S].

[13] GB/T 17140-1997，土壤質(zhì)量鉛、鎘的測定 KI-MIBK萃取火焰原子吸收分光光度法[S].

[14] 于云江，胡林凱，楊 彥，等.典型流域農(nóng)田土壤重金屬污染特征及生態(tài)風(fēng)險評價[J].環(huán)境科學(xué)研究，2010，23（12）：1523-1527.

[15] 成杰民，張麗娜.近20年山東省典型農(nóng)田土壤中重金屬含量變化特征及積累速率估計[J].土壤，2013，45（1）：99-104.

[16] 賈振邦，梁 濤，林健枝，等.香港河流重金屬污染及潛在生態(tài)危害研究[J].北京大學(xué)學(xué)報（自然科學(xué)版），1997，33（4）：485-492.