華北平原某地農(nóng)田土壤重金屬污染特征及來源分析

趙永超 龔巍崢 任岳 張孝娟

摘要 選取華北平原某地（全國重要的有色金屬產(chǎn)業(yè)基地）作為研究區(qū)，采集377個表層土壤樣品（0～20 cm），分析測定樣品中Pb、Cd、Cr、As、Hg 5種重金屬的含量。運用單因子指數(shù)法、內梅羅綜合指數(shù)法和潛在生態(tài)風險指數(shù)法評價研究區(qū)農(nóng)田土壤中重金屬的污染及生態(tài)風險狀況，運用因子分析的方法對研究區(qū)農(nóng)田土壤中重金屬元素的來源進行分析。結果表明，研究區(qū)農(nóng)田土壤中Pb、Cd、As含量超過農(nóng)用地土壤污染風險篩選值（GB 15168—2018）的點位占比分別為6.10%、44.30%和1.33%;基于農(nóng)用地土壤污染風險篩選值，研究區(qū)農(nóng)田土壤中Cd、Pb的污染等級為重度污染，As的污染等級為中度污染，Cr、Hg為安全級;基于農(nóng)用地土壤污染風險篩選值，研究區(qū)農(nóng)田土壤中Cd為中等生態(tài)風險，Pb、Cr、As、Hg的生態(tài)風險較低;綜合潛在生態(tài)危害指數(shù)的均值為56.64，總體看研究區(qū)農(nóng)田土壤處于輕微生態(tài)危害水平，其中Cd是研究區(qū)農(nóng)田土壤中最主要的污染因子和生態(tài)風險因子，貢獻率達80%以上;相關分析和主成分分析結果表明，研究區(qū)農(nóng)田土壤中As主要來源于自然源，受成土母質的控制;Cd、Pb主要來源于人類活動中的工農(nóng)業(yè)源;Hg主要來源于人類活動中工業(yè)源;Cr受成土因素和人為活動的共同影響。

關鍵詞 土壤重金屬;含量測定;污染特征;來源解析;生態(tài)風險評價

中圖分類號 X53? 文獻標識碼 A? 文章編號 0517-6611（2022）05-0055-06

doi：10.3969/j.issn.0517-6611.2022.05.015

開放科學（資源服務）標識碼（OSID）：

Analysis of Characteristics and Sources of Heavy Metal Pollution in Farmland Soil in the North China Plain

ZHAO Yong-chao， GONG Wei-zheng， REN Yue et al

（Henan Institute of Geo-Environment Exploration， Zhengzhou，Henan 450000）

Abstract A certain place in the North China Plain （an important nonferrous metal industry base in the country） was selected as the research area， 377 surface soil samples （0-20 cm） were collected， and the contents of 5 heavy metals （Pb， Cd， Cr， As， and Hg） were analyzed and determined.Single factor pollution index， the Nemerov pollution index method and the potential ecological risk index method were used to evaluate the pollution and ecological risk status of heavy metals in the study area. The factor analysis method was used to analyze the sources of heavy metal elements in farmland soils in the study area.The results showed that the percentages of points where the content of Pb， Cd and As in the farmland soil of the study area exceed the soil pollution risk screening value （GB 15168—2018） of the farmland in the study area were 6.10%， 44.30% and 1.33% respectively.Based on the soil pollution risk screening value of farmland， the pollution level of Cd and Pb in the farmland soil of the study area was severe pollution， the pollution level of As was moderate pollution， and the pollution level of Cr and Hg was safety level.Based on the soil pollution risk screening value of farmland， Cd in the farmland soil of the study area was a medium ecological risk， and the ecological risk of Pb， Cr， As， and Hg was low.The average comprehensive potential ecological risk index was 56.64， indicating that the farmland soil in the study area was at a slight ecological hazard level.Among them， Cd was the most important pollution factor and ecological risk factor in the farmland soil in the study area， with a contribution rate of more than 80%.The results of correlation analysis and principal component analysis show that As in the farmland soil of the study area was mainly derived from natural sources and was controlled by soil-forming parent material;Cd and Pb were mainly derived from industrial and agricultural sources in human activities;Hg was mainly derived from industrial sources in human activities;Cr was affected by both soil-forming factors and human activities.

Key words Soil heavy metals;Content determination;Pollution characteristic;Source apportionment;Ecological risk assessment

作者簡介 趙永超（1989—），男，河南平頂山人，助理工程師，碩士，從事重金屬污染農(nóng)田修復與耕地質量提升研究。*通信作者，高級工程師，碩士，從事污染土壤修復等研究。

收稿日期 2021-06-15

2014年環(huán)境保護部發(fā)布的《全國土壤污染狀況調查公報》顯示，全國土壤環(huán)境狀況總體不容樂觀，部分地區(qū)土壤污染較重，耕地土壤點位超標率為19.4%，主要污染物為Cd、Ni、Cu、As、Hg、Pb、滴滴涕和多環(huán)芳烴。

土壤中的重金屬一方面來源于土壤母質的成土過程，另一方面來源于人類活動。不同成土母質形成的土壤中重金屬的含量差異較大，如成土母質為海河沉積物和花崗巖的土壤，Pb含量高于其他母質發(fā)育的土壤，花崗巖發(fā)育形成的土壤Cd含量低于其他母質發(fā)育的土壤[1];不同人類活動對重金屬的貢獻也存在一定的差異，如Cd主要來源于化肥、農(nóng)藥和電鍍等行業(yè)“三廢”的排放[2]，Pb主要來源于汽車尾氣和煤炭的燃燒等。重金屬元素在土壤中長期的積累一方面可以抑制植物和微生物的生長活動導致土壤肥力的下降，從而造成農(nóng)作物產(chǎn)量大幅降低和品質的下降[3];另一方面重金屬元素在土壤中的遷移和在食物鏈中的富集還會對區(qū)域內地下水質量安全、生態(tài)安全及人類的生存發(fā)展產(chǎn)生直接的影響[4]。因此，揭示土壤中重金屬的來源和生態(tài)風險，可以為區(qū)域土壤環(huán)境質量評價、污染防治和安全防護提供借鑒。

冶煉區(qū)和礦區(qū)被認為是土壤重金屬污染的多發(fā)區(qū)[5]，筆者選取華北平原某地（全國重要的有色金屬產(chǎn)業(yè)基地）作為研究區(qū)，分析測定土壤中Pb、Cd、Cr、As、Hg 5種重金屬的含量，在此基礎上，對區(qū)域土壤重金屬的污染程度、潛在生態(tài)風險進行評價，并運用主成分分析的方法辨識土壤重金屬的來源，以期為研究區(qū)農(nóng)用地的安全利用、土壤環(huán)境評價與管理提供依據(jù)。

1 材料與方法

1.1 研究區(qū)概況

研究區(qū)位于華北平原西部，地勢西高東低，呈階梯狀展布，絕對高程大部分在70～100 m，西部為低山，西南部是丘陵區(qū)，中東部為沖積、洪積平原，耕地面積16 596.77 hm2。屬于暖溫帶半濕潤大陸性季風氣候，夏季暖熱多雨，冬季寒冷干燥，四季分明，年平均氣溫14.9 ℃，全年平均降雨量為606.1 mm，主要土壤類型為褐土。

1.2 樣品采集與分析

于玉米收獲前進行土壤樣品的采集。利用ArcGIS 10軟件按照500 m×500 m的規(guī)則網(wǎng)格進行土壤采樣點的布設，同時綜合考慮種植方式、土壤類型和重點行業(yè)企業(yè)（重點行業(yè)企業(yè)周圍加密布點）分布狀況，共布設377個采樣點;在各采樣點周圍按照多點采樣的方式，采集表層（0～20 cm）土壤，等量混合后，裝入自封袋。

采集到的土壤樣品經(jīng)室溫自然風干后，去除其中的雜質，用木棍碾碎，使其全部通過2 mm的尼龍篩后充分混合待測。Cd、Pb采用石墨爐原子吸收分光光度法（GB/T 17141—1997）測定，Cr采用火焰原子吸收分光光度法（HJ 491—2009）測定，Hg、As采用原子熒光分光光度法（GB/T 22105—2008）測定。

隨機選取采集到的10個土壤樣品，測定土壤pH、有機質、全氮、全磷、全鉀等基本理化性質，測定方法參照魯如坤[6]的《土壤農(nóng)業(yè)化學分析方法》，具體測定結果為pH 7.51～7.80、有機質2.06～9.76 g/kg、全氮0.473～0.908 g/kg、全磷1.78～2.22 g/kg、全鉀 18.4～26.5 g/kg。

1.3 土壤重金屬污染評價方法

采用單因子指數(shù)法和內梅羅綜合污染指數(shù)法對研究區(qū)農(nóng)田土壤（表層）中重金屬的污染狀況進行評價。

單項污染指數(shù)可以反映不同重金屬元素的污染程度，其計算公式如下：

P i=C i/S i（1）

式中，P i為土壤中重金屬i的單項污染指數(shù);

C i為土壤中重金屬i的實測含量（mg/kg）;

S i為農(nóng)用地土壤重金屬i的風險篩選值（mg/kg），參考GB 15168—2018規(guī)定的農(nóng)用地土壤污染風險篩選值，具體見表1。

單因子污染指數(shù)僅能反映土壤中單個污染物的污染程度，為綜合反映不同污染物的整體污染狀況，該研究采用內梅羅綜合指數(shù)法。內梅羅綜合污染指數(shù)可以全面反映各污染物對土壤的不同作用，并且可以突出高濃度污染物對土壤環(huán)境質量的影響，計算公式如下[7]：

PI=P max2+P ave22（2）

式中，PI為內梅羅綜合指數(shù);

P max為單因子污染指數(shù)的最大值;

P ave為單因子污染指數(shù)的平均值。單因子污染指數(shù)法和內梅羅綜合污染指數(shù)法污染等級劃分標準如表2所示。

1.4 土壤重金屬生態(tài)風險評價方法

采用潛在生態(tài)風險指數(shù)法對研究區(qū)農(nóng)田土壤（表層）中重金屬的潛在生態(tài)風險進行評價。該方法由瑞典科學家Hakanson[8]創(chuàng)建，該方法綜合考慮了環(huán)境中重金屬的生態(tài)效應、環(huán)境效應和毒理學效應，定量計算環(huán)境中重金屬的潛在危害程度，被廣泛應用于環(huán)境風險評價中，其計算公式如下：

Ei r=Ti r×Ci f（3）

RI=Ei r=（Ti r×Ci f）（4）

Ci f=Ci/Ci n（5）

式中，Ei r為土壤或沉積物中元素i的潛在生態(tài)危害指數(shù);

RI為土壤或沉積物中多種重金屬元素潛在生態(tài)危害指數(shù);

Ci f為重金屬元素i的富集系數(shù);Ci為重金屬元素i的實測含量;

Ci n為重金屬元素i的風險篩選值，參考GB 15618—2018規(guī)定的農(nóng)用地土壤污染風險篩選值;

Ti r為元素i的毒性響應系數(shù)，參考已有文獻資料Pb、Cd、Cr、As、Hg分別為5、30、2、10、40[9]。Ei r和RI的分級標準見表3。

1.5 數(shù)據(jù)統(tǒng)計分析

運用SPSS 19.0軟件對研究區(qū)農(nóng)田土壤中5種重金屬的含量進行描述性統(tǒng)計分析，以分析研究區(qū)不同重金屬的含量特征;采用相關分析和主成分分析法對研究區(qū)土壤中重金屬的來源進行解析。

2 結果與分析

2.1 土壤重金屬含量特征

運用SPSS軟件對研究區(qū)農(nóng)田土壤中Pb、Cd、Cr、As、Hg的含量進行統(tǒng)計分析，分析結果如表4所示。從表4可以看出，研究區(qū)農(nóng)田土壤中Pb、Cd、Cr、As、Hg含量的平均值分別為63.52、0.97、61.48、11.40、0.08 mg/kg，Pb、Cd、As、Hg含量的平均值超過河南省土壤背景值[10]，超標倍數(shù)依次為1.85、13.92、0.16、2.20。土壤中Cr含量的平均值低于河南省土壤背景值[10]，說明土壤中Cr沒有明顯的積累。Pb、Cd含量的平均值大于其中位數(shù)，說明研究區(qū)農(nóng)田土壤中Pb、Cd的含量整體偏向最大值方向;Cr、As、Hg含量的平均值與中位數(shù)接近，說明研究區(qū)農(nóng)田土壤中Cr、As、Hg的含量整體偏向最小值方向。偏度和峰度的數(shù)值可大致反映數(shù)據(jù)偏離正態(tài)分布的程度，研究區(qū)農(nóng)田土壤中Pb、Cd、Cr、As、Hg的偏度和峰度絕對值均大于0，說明研究區(qū)農(nóng)田土壤中重金屬濃度的分布為非標準正態(tài)分布。研究區(qū)農(nóng)田土壤中Pb、Cr、As、Hg含量的平均值低于農(nóng)用地土壤污染風險篩選值（GB 15168—2018），Cd含量的平均值高于風險篩選值，說明研究區(qū)農(nóng)田土壤中Cd受人類活動的影響較明顯。土壤中Cd的最大值是農(nóng)用地土壤風險篩選值的20倍，說明Cd在部分點位發(fā)生了明顯的積累，這可能與研究區(qū)內存在重點污染源密切相關。

土壤重金屬元素的變異系數(shù)可以代表重金屬元素分布的均勻程度和受人類活動影響的程度。研究區(qū)農(nóng)田土壤中Pb、Cd和Hg的變異系數(shù)大于1，屬于強變異水平;Cr、As的變異系數(shù)介于0.1～1.0，屬于中等變異水平;上述結果表明研究區(qū)農(nóng)田土壤中Pb、Cd和Hg的分布不均勻，可能受人類活動的影響較大。

2.2 土壤重金屬污染狀況評價

基于《土壤環(huán)境質量 農(nóng)用地土壤污染風險管控標準》（GB 15618—2018）中規(guī)定的各項污染物的風險篩選值，計算各點位5種重金屬的單項污染指數(shù)，依據(jù)表2的劃分標準對計算結果進行統(tǒng)計分析，結果如表5所示。從表5可以看出，研究區(qū)農(nóng)田土壤中Cr、Hg的單項污染指數(shù)的最大值分別為0.470、0.590，均小于1，依據(jù)劃分標準Cr、Hg屬于非污染等級。Pb、Cd、As的單因子污染指數(shù)分別為0.070～4.430、0.100～20.000、0.070～3.240，總體看Pb、Cd、As這3種重金屬的污染程度高于Cr、Hg。統(tǒng)計結果表明研究區(qū)農(nóng)田土壤中Pb、Cd、As達到輕度及以上污染等級的點位數(shù)占比分別為6.10%、44.30%、1.33%，其中Cd達到中度及重度污染等級的點位數(shù)占比較大，達到21.75%，說明研究區(qū)農(nóng)田土壤中Cd的污染程度已經(jīng)遠高于其他重金屬，需要引起高度重視。

基于《土壤環(huán)境質量 農(nóng)用地土壤污染風險管控標準》（GB 15618—2018）中規(guī)定的各項污染物的風險篩選值，以綜合污染指數(shù)進行評價，結果表明（表6），研究區(qū)農(nóng)田土壤中Pb、Cd、Cr、As、Hg的綜合污染指數(shù)（PI）分別為3.14、14.19、0.38、2.31、0.42，結合表2的分級標準，研究區(qū)農(nóng)田土壤中Cd、Pb的污染等級為重度污染，As的污染等級為中度污染，Cr、Hg為安全級，上述結果說明研究區(qū)農(nóng)田土壤已受到一定程度的重金屬污染。其中Cd的綜合污染指數(shù)最大，說明研究區(qū)農(nóng)田土壤中Cd的污染最為突出，這與單項污染指數(shù)推斷的結果一致。

對各點位的綜合污染指數(shù)進行計算，依據(jù)表2的分類標準對各點位的污染程度進行劃分，并對分類結果進行統(tǒng)計分析，結果發(fā)現(xiàn)，研究區(qū)內有51.20%的點位受到污染，其中有14.33%的點位呈中度至重度污染。研究區(qū)所有調查點位PI的平均值為1.24，依據(jù)綜合污染指數(shù)分級標準，研究區(qū)農(nóng)田土壤污染等級為輕度污染，說明研究區(qū)農(nóng)田土壤環(huán)境質量已經(jīng)開始惡化，應引起高度重視。

2.3 土壤重金屬的潛在生態(tài)風險評價

以《土壤環(huán)境質量 農(nóng)用地土壤污染風險管控標準》（GB 15618—2018）中規(guī)定的各項污染物的風險篩選值作為評價參比值，采用潛在生態(tài)風險指數(shù)法對研究區(qū)農(nóng)田土壤中重金屬的潛在生態(tài)風險進行評價，結果如表7所示。從表7可以看出，研究區(qū)農(nóng)田土壤中Pb、Cd、Cr、As、Hg這5種重金屬的潛在生態(tài)危害指數(shù)（Ei r）分別為1.87、48.75、049、4.56、0.97，依據(jù)表3的分級標準，研究區(qū)農(nóng)田土壤中Pb、Cr、As、Hg的生態(tài)風險較低，均處于輕微生態(tài)危害水平;Cd為中等生態(tài)風險。研究區(qū)農(nóng)田土壤中Cd各危害水平的點位占比分別為60.75%（輕微生態(tài)危害水平）、24.93%（中等生態(tài)危害水平）、9.55%（強生態(tài)危害水平）、3.18%（很強生態(tài)危害水平）、1.59%（極強生態(tài)危害水平）。

各點位綜合潛在生態(tài)危害指數(shù)（RI）的統(tǒng)計結果顯示，研究區(qū)有93.10%的點位危害水平較低，有6.90%的點位處于中等及以上生態(tài)危害水平。研究區(qū)各點位RI的平均值為56.64，結合表3的分級標準，研究區(qū)農(nóng)田土壤的生態(tài)風險為輕微生態(tài)危害。

為評價單個重金屬元素的潛在生態(tài)危害指數(shù)對綜合生態(tài)危害指數(shù)的貢獻率，取單項潛在生態(tài)危害指數(shù)與綜合潛在生態(tài)危害指數(shù)的比值，各元素對綜合潛在生態(tài)危害的貢獻率如圖1所示。從圖1可以看出，Cd對綜合潛在生態(tài)危害的貢獻率最大，達86.07%，是研究區(qū)農(nóng)田土壤生態(tài)風險的主要來源[11]。

2.4 土壤重金屬來源分析

相關性分析是指對2個或多個具備相關性的變量元素進行分析，從而衡量2個變量因素的相關密切程度。土壤中各種元素的含量受母質、環(huán)境條件和人為活動的共同影響，來源的相似性會導致土壤中不同元素

間存在相互聯(lián)系[12]，鑒于此相關性分析就被廣泛應用于土壤重金屬元素的來源辨識[13]，若元素間相關性顯著或極顯著，則表明這些元素間一般具有同源關系或者呈現(xiàn)復合污染的現(xiàn)象[14-17]。

2.4.1 相關性分析。

從表8可以看出，研究區(qū)農(nóng)田土壤中Pb與Cd、Pb與Cr、Pb與As、Pb與Hg、Cd與As、Cd與Hg、Cr與As在0.01水平（雙側）上呈顯著正相關，說明這些元素間可能存在著一定的同源關系。已有的研究表明，土壤中Cr主要受成土母質的影響，因此Cr多被認為是自然來源的元素，與其相關性較高的元素被認為受成土母質的影響較大[18]，研究區(qū)農(nóng)田土壤中Cr與Pb、As呈顯著正相關，說明Cr、Pb、As可能具有一定的同源性，而Cr與Pb的相關系數(shù)較小，因此研究區(qū)農(nóng)田土壤中Pb是否存在其他來源還需要進一步驗證。土壤中Hg與Pb、Cd、Cr、As的相關系數(shù)均不高于0.30，說明Hg和Pb、Cd、Cr、As相關性不大，從而說明Hg與Pb、Cd、Cr、As的來源可能不一致[19]。為更加準確地判斷土壤中各重金屬元素的來源，通過主成分分析來進行進一步的判斷。

2.4.2 主成分分析。

從主成分分析的結果（表9～10和圖2）可以看出，通過主成分分析共辨別出3個主成分，累計解釋了總方差的84.589%，基本可以解釋5種元素的大部分信息。

第一主成分（PC1）的方差貢獻率為37.219%，Pb、Cd在PC1上有較大的載荷，分別為0.906、0.894，由前述分析可知，研究區(qū)農(nóng)田土壤中Pb、Cd含量的平均值高于該地區(qū)背景值，2種元素含量的最大值分別為河南省背景值的33.77和184.62倍;且前人的研究發(fā)現(xiàn)，這2種元素常常受到人類活動的影響，其中Pb主要來自化石燃料的燃燒和金屬冶煉行業(yè)“三廢”的排放，Cd主要來自農(nóng)業(yè)生產(chǎn)資料（化肥、農(nóng)藥等）的施用及冶煉、電鍍行業(yè)“三廢”的排放[20]。研究區(qū)作為華北平原主要的農(nóng)業(yè)生產(chǎn)大縣，近年來隨著農(nóng)業(yè)精細化程度的不斷提高，化肥和農(nóng)藥的施用量大增;另外該地區(qū)工業(yè)基礎較好，是全國重要的金屬冶煉區(qū)，化肥、農(nóng)藥和工業(yè)“三廢”中攜帶的Cd、Pb通過直接或間接的方式最終進入到土壤中，造成Pb、Cd在土壤中的累積。因此PC1代表了化肥農(nóng)藥及工業(yè)生產(chǎn)等人為來源，被認為是“工農(nóng)業(yè)源因子”。

第二主成分（PC2）的方差貢獻率為27.268%，Cr、As在PC2上有較大的載荷，分別為0.700、0.907，前述分析結果顯示研究區(qū)Cr、As的平均值與河南省背景值相當。已有的研究發(fā)現(xiàn)Cr、As受地質背景的影響較大[21]，因此把PC2認為是“自然源因子”。

第三主成分（PC3）的方差貢獻率為20.102%，PC3上有較大載荷的僅有Hg。研究區(qū)農(nóng)田土壤中Hg的平均值為008 mg/kg，遠高于河南省背景值，其最大值是河南省背景值的80.4倍。相關分析的結果顯示Hg與Cr、As呈正相關，但是Hg與Cr、As的相關系數(shù)較低，說明研究區(qū)農(nóng)田土壤中Hg與Cr、As的同源性較低。已有的研究表明，土壤是各種化學離子的最終歸宿[22]，Hg在土壤中累積的主要途徑是大氣干濕沉降[23]，其來源主要包括人類活動中的化石燃料燃燒、冶煉、化工、垃圾焚燒等，研究區(qū)的工業(yè)基礎較好，金屬冶煉企業(yè)在生產(chǎn)過程中產(chǎn)生的大量粉塵等最終都匯聚到土壤中造成了研究區(qū)農(nóng)田土壤中Hg的積累。因此PC3代表了人類活動中的工業(yè)源，被認為是“工業(yè)源因子”。

Cr在PC1和PC2上均有相當?shù)暮奢d，說明研究區(qū)農(nóng)田土壤中Cr是混合來源，被認為是受到自然源和人為源（工農(nóng)業(yè)源因子）的共同影響。于元赫等[24]對黃河下游典型區(qū)域土壤重金屬來源的研究也表明Cr受工農(nóng)業(yè)活動和成土母質的綜合影響。

3 結論

（1）研究區(qū)農(nóng)田土壤中Pb、Cr、As、Hg的平均含量高于河南省土壤背景值，Cr含量的平均值低于河南省土壤背景值。研究區(qū)內Pb、Cd、As含量達到輕度及以上污染等級的點位數(shù)占比分別為6.10%、44.30%和133%。

（2）研究區(qū)農(nóng)田土壤中Pb、Cd、Cr、As、Hg的內梅羅綜合污染指數(shù)分別為3.14、14.19、0.38、2.31、0.42，依據(jù)分級標準，Cd、Pb的污染等級為重度污染，As的污染等級為中度污染，Cr、Hg的污染等級為安全級。

（3）單項潛在生態(tài)風險指數(shù)（Ei r）的結果表明，研究區(qū)農(nóng)田土壤中Cd為中等生態(tài)風險，Pb、Cr、As、Hg的生態(tài)風險較低。綜合潛在生態(tài)危害指數(shù)（RI）的結果表明，研究區(qū)9310%的調查點位的生態(tài)風險較低，6.90%的調查點位為中等及以上危害水平;整個研究區(qū)農(nóng)田土壤綜合潛在生態(tài)危害等級為輕微生態(tài)危害水平。Cd是影響研究區(qū)農(nóng)田土壤綜合潛在生態(tài)風險的主要元素。

（4）相關分析和主成分分析結果表明，研究區(qū)農(nóng)田土壤中5種元素的來源可以分為3類：PC1為工農(nóng)業(yè)源因子，PC1包括Pb、Cd;PC2為自然源因子，PC2包括Cr、As;PC3為工業(yè)源因子，僅包括Hg;Cr在PC1和PC2上均有相當?shù)暮奢d，Cr受成土因素及人為活動的共同控制。

參考文獻

[1] LI Y，GOU X，WANG G，et al.Heavy metal contamination air source in air agricultural soil in central Gansu Province，China[J].Journal of environmental sciences，2008，20（5）：607-612.

[2] 章明奎，王浩，張慧敏.浙東海積平原農(nóng)田土壤重金屬來源辨識[J].環(huán)境科學學報，2008，28（10）：1946-1954.

[3] 邱孟龍，李芳柏，王琦，等.工業(yè)發(fā)達城市區(qū)域耕地土壤重金屬時空變異與來源變化[J].農(nóng)業(yè)工程學報，2015，31（2）：298-305.

[4] WEN H J，ZHANG Y X，CLOQUET C，et al.Tracing sources of pollution in soils from the Jinding Pb-Zn mining district in China using cadmium and lead isotopes[J].Applied geochemistry，2015，52：147-154.

[5] 趙振，史燁弘，毛小晶，等.某冶煉廠周邊土壤重金屬污染評價分析及源解析研究[J].中國無機分析化學，2020，10（4）：22-27.

[6] 魯如坤.土壤農(nóng)業(yè)化學分析方法[M].北京：中國農(nóng)業(yè)科技出版社，1999.

[7] 呂悅風，孫華.浙北某縣域耕地土壤重金屬空間分異特征、污染評價及來源分析[J].農(nóng)業(yè)環(huán)境科學學報，2019，38（1）：95-102.

[8] HAKANSON L.An ecological risk index for aquatic pollution control.A sedimentological approach[J].Water research，1980，14（8）：975-1001.

[9] 徐爭啟，倪師軍，庹先國，等.潛在生態(tài)危害指數(shù)法評價中重金屬毒性系數(shù)計算[J].環(huán)境科學與技術，2008，31（2）：112-115.

[10] 河南省土壤普查辦公室.河南土壤[M].北京：中國農(nóng)業(yè)出版社，2004.

[11] 穆莉，王躍華，徐亞平，等.湖南省某縣稻田土壤重金屬污染特征及來源解析[J].農(nóng)業(yè)環(huán)境科學學報，2019，38（3）：573-582.

[12] 陳碧珊，蘇文華，羅松英，等.湛江特呈島紅樹林濕地土壤重金屬含量特征及污染評價[J].生態(tài)環(huán)境學報，2017，26（1）：159-165.

[13] 呂占祿，張金良，陸少游，等.某區(qū)生活垃圾焚燒發(fā)電廠周邊及廠區(qū)內土壤中重金屬元素的污染特征及評價[J].環(huán)境科學，2019，40（5）：2483-2492.

[14] 鄒曦，鄭志偉，張志永，等.三峽水庫小江流域消落區(qū)土壤重金屬時空分布與來源分析[J].水生態(tài)學雜志，2012，33（4）：33-39.

[15] 李晉昌，張紅，石偉.汾河水庫周邊土壤重金屬含量與空間分布[J].環(huán)境科學，2013，34（1）：116-120.

[16] 柴世偉，溫琰茂，張云霓，等.廣州市郊區(qū)農(nóng)業(yè)土壤重金屬含量特征[J].中國環(huán)境科學，2003，23（6）：592-596.

[17] GUO G H，WU F C，XIE F Z，et al.Spatial distribution and pollution assessment of heavy metals in urban soils from southwest China[J].Journal of environmental sciences，2012，24（3）：410-418.

[18] LIU Y，MA Z W，LV J S，et al.Identifying sources and hazardous risks of heavy metals in topsoils of rapidly urbanizing East China[J].Journal of geographical sciences，2016，26（6）：735-749.

[19] 劉青枰，鄧時強，趙元藝，等.西藏尕爾勤礦床區(qū)域性土壤重金屬元素地球化學特征及來源解析[J].地球學報，2018，39（4）：481-490.

[20] 周旭，呂建樹.山東省廣饒縣土壤重金屬來源、分布及生態(tài)風險[J].地理研究，2019，38（2）：414-426.

[21] SJN R，HALAMIC' J，PEH Z，et al.Assessment of the natural and anthropogenic sources of chemical elements in alluvial soils from the Drava River using multivariate statistical methods[J].Journal of geochemical exploration，2011，110（3）：278-289.

[22] 鄒建美，孫江，戴偉，等.北京近郊耕作土壤重金屬狀況評價分析[J].北京林業(yè)大學學報，2013，35（1）：132-138.

[23] 王梅，黃標，孫維俠，等.強烈人為作用下城鎮(zhèn)周圍汞的空間變異及其積累遷移規(guī)律[J].土壤學報，2011，48（3）：506-515.

[24] 于元赫，呂建樹，王亞夢.黃河下游典型區(qū)域土壤重金屬來源解析及空間分布[J].環(huán)境科學，2018，39（6）：2865-2874.