四川省沐川縣典型土壤剖面重金屬分布特征和生態風險評價

廖碧霞，沈文靈，賀靈

（1. 中國地質科學院地球物理地球化學勘查研究所，河北 廊坊 065000； 2. 桂林理工大學，廣西 桂林 541006）

土壤是人類生存和發展的寶貴資源。在社會經濟快速發展的今天，土壤環境質量受到多種因素的影響，耕地重金屬污染問題在世界范圍內引起了廣泛關注[1-2］。耕地重金屬含量明顯增加，主要是由于農業生產、大氣沉降、廢物處理等人為活動［3］所致。一方面，耕作土壤中重金屬的積累影響了土壤的理化性質，導致土壤養分元素循環失衡，直接影響了耕作土壤的生產力，導致作物產量和品質下降。另一方面，重金屬污染范圍廣、持續時間長、易積累難分解，可通過直接接觸或食物鏈傳遞對人類健康構成威脅［4-7］。例如，汞的積累可導致自身免疫性疾病以及肺和腎衰竭。長期攝入過量Cd 導致前列腺增生性病變、骨折、腎功能障礙、肺癌和肺腺癌等不良影響［8］。中國約2.786×109m2農業土壤被Cd 污染［9］。 Kumar 等［10］收集了1991—2018年印度重金屬污染土壤數據，所有土壤類型Cd 平均含量都超過限量值，Cd 的潛在生態風險值和修正后的潛在生態風險值均大于320，反映出較高的生態風險。研究土壤中重金屬的含量并評估其風險，對于避免重金屬對生態環境和人類健康造成危害具有重要意義。

近年來國內外學者在不同區域，開展自然條件、工礦業及交通發達等不同條件下，土壤重金屬污染及生態風險評價等研究。周亞龍等［11］發現雄安新區受周邊企業生產活動的影響，部分根系土壤樣品As、Cd、Cu、Pb 和Zn 含量超過農用地土壤污染風險篩選值標準(GB 15618—2018)，其超標率分別為23.33%、96.67%、33.33%、33.33% 和10.00%。宋綿等［12］評估了河北阜平縣表層土壤重金屬污染特征發現As 和Cd 超過可接受致癌風險水平(As 為10-5，Cd 為10-6)，具有一定的致癌風險。Huang 等［13］發現九龍江流域水稻土中重金屬具有較高的潛在生態風險，其中鎘和汞的貢獻率最高，分別為59.4%和26.2%。Barrena-González 等［14］研究發現西班牙埃斯特雷馬杜拉地區的土壤養分表現出較高的變異性，80%以上的養分集中在土壤上層和較淺的土層。其他化學性質，如pH 值表現出相反的行為，在深層具有較小的變異性和較高的值。謝龍濤等［15］研究表明土壤重金屬含量與土壤有機質和pH 值密切相關。Ye 等［16］利用聚合增強樹分析表明，土壤pH 和有機碳是控制水稻土Cd 生物有效性的主要因素。

四川省沐川縣位于中國西南部，是典型的土壤重金屬元素地質高背景區。目前，對該區域內土壤剖面重金屬垂直分布及其生態風險的研究認識還不夠全面。本文選取沐川縣3 個剖面的38 個土壤樣品中的7 種重金屬(As、Cd、Cu、Zn、Hg、Ni、Pb)為對象，利用電感耦合等離子體質譜法、氫化物發生原子熒光光譜法、高頻燃燒紅外吸收法等方法測定樣品中7 種重金屬元素含量等指標，采用地質累積指數(Geoaccumulation Index，Igeo)和漢克森潛在生態風險指數法(Potential Ecological Risk Index，RI)，研究重金屬在土壤中的垂向分布特征，探究沐川縣土壤重金屬與土壤養分元素的關系及其污染程度和潛在生態風險，擬為地方政府對土壤資源的安全利用、土壤污染的防治與修復提供依據。

1 研究區概況

研究區位于四川盆地西南邊緣小涼山余脈、烏蒙山區西北部，地處岷江、大渡河、金沙江的腹心地帶，隸屬于沐川縣。地貌主要以平壩、低山和低中山為主［17］。氣候屬亞熱帶濕潤季風氣候。主要作物為玉米、茶葉等。土壤類型多為紫色土和黃壤［18］，區內土地利用方式以旱地、茶園為主，耕地較為分散，主要為山間梯田。

在地質構造上，沐川縣位于揚子淮地臺四川臺坳川中臺拱西南緣的弧形褶帶，西南接涼山凹陷帶、東北連威遠臺凸，縣內斷裂少但褶皺較發育。區內出露地層主要有：侏羅系蓬萊鎮組(J3p)、遂寧組(J3sn)、沙溪廟組(J2s)和自流井組(J1z)；三疊系須家河組(T3x)和垮洪洞組—須家河組并層(T3k-x)；三疊系雷口坡組(T2l)；三疊系飛仙關組—嘉陵江組—雷口坡組并層(T1-2f-l) 和飛仙關組—嘉陵江組并層(T1f-j)；二疊系宣威組(P2x)；二疊系峨眉山玄武巖組(Pe)；白堊系窩頭山組(K1w)(圖1) 。蓬萊鎮組(J3p)、遂寧組(J3sn)、沙溪廟組(J2s) 和自流井組(J1z) ：泥巖、砂巖和粉砂巖；須家河組(T3x) 和垮洪洞組—須家河組并層(T3k-x)：砂巖、粉砂巖、泥巖及煤層組成的沉積旋回；雷口坡組(T2l)：粉砂巖與白云巖、泥質灰巖互層，夾黑色碳質頁巖；飛仙關組—嘉陵江組—雷口坡組并層(T1-2f-l) 和飛仙關組—嘉陵江組并層(T1f-j)：白云巖、灰巖、頁巖；宣威組(P2x)：黃綠色砂、粉砂巖，夾泥巖及煤；峨眉山玄武巖組(Pe)：灰綠色致密、杏仁狀玄武巖，夾苦橄巖、凝灰質砂泥巖、煤線及硅質巖；窩頭山組(K1w)：磚紅色塊狀細粒長石石英砂巖。

圖1 研究區位置和地質簡圖Fig. 1 Location and geological map of the study area.

2 實驗部分

2.1 樣品采集

按照不同的地質背景和土壤類型布設三條土柱剖面(YS、PS、GS)。采用手動取樣鉆連續分段采樣，樣柱長度控制為10cm/件。每件樣品質量不少于500g。取深部土樣時，需防止上部土壤混入。土柱剖面深度視土層發育情況而定，從3 個土壤剖面共采集38 個土壤樣品，采樣點描述見表1。土壤樣品經自然風干后去除石塊、碎屑、塑料等雜質，研磨后過10 目篩，混勻、稱重并送往實驗室。

表1 研究區不同類型土柱剖面取樣點概況Table 1 Sampling points of different types of soil column profiles in the study area.

2.2 樣品測試

以《土地質量地球化學評價規范》(DZ/T 0295—2016)為依據 ，采用電感耦合等離子體質譜法(ICP-MS，ICap Qc 型，美國ThermoFisher 公司)測定Cd、Cu、Ni、Pb、Zn 含量；氫化物發生原子熒光光譜法(HG-AFS，XGY-2020A 型，廊坊開元高技術開發公司)測定As 含量；X 射線熒光光譜法(XRF，PW 4400/40 型，荷蘭帕納科公司) 測定P、K2O 含量；氧化燃燒-氣相色譜法(GC，EA3000 型，意大利Eurovetorspa 公司)測定N 含量；冷蒸氣原子熒光光譜法(CV-AFS，XGY-2020A 型，廊坊開元高技術開發公司)測定Hg 含量；高頻燃燒紅外吸收法(IR，COREY-200 型，四川普瑞晟科技有限公司) 測定有機碳(Corg) 含量；電位法(POT，PHB-4 型，上海儀電科學儀器有限公司)測定土壤pH 值。各指標分析測試的檢出限見表2。

表2 各指標分析測試檢出限Table 2 Detection limit of each index analysis.

分析質量控制嚴格執行《土地質量地球化學評價規范》(DZ/T 0295—2016)和《多目標區域地球化學調查規范》(DZ/T 0258—2014) 。樣品分析過程中，按照規范插入國家一級標準物質(GBW07401、GBW07408、GBW07423～GBW07426) 和重復樣，用于監控分析測試的準確度和精密度。各元素的分析檢出限、準確度和精密度均達到或優于DZ/T 0258—2014 的要求，分析數據質量可靠。

2.3 污染程度和潛在生態風險評價方法

采用地質累積指數(Igeo)［19］評價土壤重金屬污染程度，選取潛在生態風險指數(RI)［20］評價潛在生態風險。Igeo的計算公式如下：

式中：Cn為土壤中目標元素的實測值，Bn為土壤中目標元素的背景值。在本研究中，Bn選取四川省土壤背景值［21］。k是一個常數，用于修正不同地方巖石組成差異引起的背景值的變化(k一般賦值為1.5)。

Igeo與污染程度之間的對應關系為：Igeo≤0 表示無污染； 0＜Igeo≤1 表示輕度至中等污染； 1＜Igeo≤2表示中等污染； 2 ＜Igeo≤3 表示中等至強污染； 3＜Igeo≤4 表示強污染； 4＜Igeo≤5 表示強至極嚴重污染；Igeo＞5 表示極嚴重污染。

式中：RI為所有重金屬元素個體潛在生態風險的總和；Ei為單個元素的潛在生態風險指數；Tm為生物毒性因子。依據 Hakanson［20］和徐爭啟等［22］的研究，給出的毒性響應系數為：Hg=40，Cd=30，As=10，Pb=Cu=Ni=5，Zn=1 和Cr=2。

潛在生態危害系數(Ei) 從低到高分為5 個等級，用于描述某一污染物(元素) 的生態危害程度。不同的Ei值對應的生態危害程度劃分為：Ei＜40 為輕微；40≤Ei≤80 為中等；80≤Ei≤160 為強；160≤Ei≤320 為很強；Ei≥320 為極強。

潛在生態危害指數(RI) 是描述某一點多個污染物(元素) 潛在生態危害系數的綜合值，分為4 個等級。RI＜150 表示輕微；150≤RI≤300 表示中等；300≤RI≤600 表示強；RI≥600 表示很強。

2.4 數據處理與圖件制作

采用Excel 2016 進行土壤重金屬含量、污染指數、生態風險指數的統計和計算，利用ArcGIS 繪制研究區地質簡圖，采用SPSS 26 進行Pearson 相關分析和箱型圖制作，后期用 CorelDRAW2019 進行修飾。

3 結果與討論

3.1 土柱剖面重金屬含量

三個剖面重金屬(As、Cd、Cu、Hg、Ni、Pb、Zn)平均含量對比如圖2。

圖2 三個采樣點土壤剖面重金屬平均含量對比Fig. 2 Comparison of average heavy metal contents in soil profiles at three sampling points. The black dotted line indicates the background value of national soil, and the red dotted line indicates the background value of Sichuan soil.

剖面YS 土壤中，As、Cd、Cu、Hg、Ni、Pb、Zn含量范圍分別為(mg/kg)：18.07 ～22.81、0.22～0.73、22.96 ～31.13、0.088 ～0.148、30.32 ～36.85、37.32～53.52 和84.54 ～104.10；變異系數分別為0.07、0.36、0.09、0.17、0.06、0.11、0.07。

剖面PS 土壤中，As、Cd、Cu、Hg、Ni、Pb、Zn含量范圍分別為(mg/kg)：6.09～8.35、0.21～0.49、24.87～31.71、0.02～0.04、45.02～51.98、29.57～42和 107.70～125.60；變異系數分別為0.09、0.31、0.08、0.23、0.04、0.11、0.05。

剖面YS 和PS 土壤中，Cd 呈高度變異，Hg、Pb呈中等變異，As、Cu、Pb、Zn 呈低變異程度。多種重金屬含量多呈低變異性，表明As、Cu、Pb、Zn 主要來源于母質。

剖面GS 土壤中，As、Cd、Cu、Hg、Ni、Pb、Zn含量范圍分別為(mg/kg)：0.95～10.64、0.06～0.31、14.39～31.01、0.038～0.119、22.46～48.34、22.28～43.68 和59.54～105.80；變異系數分別為0.71、0.47、0.19、0.41、0.30、0.24、0.17。As、Cd、Hg 呈高變異性，表明該研究區這三種金屬含量變異性高，空間分散程度大，可能受到人類活動的影響。重金屬的空間分布可能是多種污染源共同作用的結果，但更多的污染源信息需要通過不同的統計分析來深入挖掘。

As 含量在剖面YS 高于全國和四川背景值，在剖面PS、GS 均低于全國和四川背景值。三個剖面的Cd 含量都高于全國和四川背景值，Cu 含量均低于四川背景值；剖面YS 和GS 的Hg 含量高于全國和四川背景值，剖面PS 的Hg 含量顯著低于全國和四川背景值；Ni 含量在剖面PS 遠高于背景值，在剖面YS 接近四川背景值，高于全國背景值，在剖面GS 高于全國背景值，低于四川背景值；Pb 含量在三個剖面均高于背景值；Zn 含量在剖面YS、PS 都高于背景值，在剖面GS 接近全國土壤背景值，低于四川背景值。三個剖面Zn 含量最高，Cd、Hg 含量最低。這7 種重金屬除了Cu，其余金屬含量都高于國家和四川土壤背景值，表明這些金屬在沐川縣土壤中較富集。成土母質是重金屬富集的重要驅動機制［23］，不同成土母質發育的土壤，其元素含量必定具有一定差異［7］。整體上看，不同土地利用方式下，玉米地(剖面YS、PS)重金屬(Cd、Cu、Pb、Zn)含量高于茶園地(剖面GS)。

3.2 土壤剖面重金屬垂向分布

As、Cd、Cu、Hg、Ni、Pb、Zn 等重金屬元素在三個剖面上的分布如圖3 所示。YS 剖面Cd、Cu、Zn、Pb 隨土層深度的增加而降低，由于作物根系主要分部局在表土中，相對較高的重金屬含量可能會引起人們對糧食安全的擔憂；Hg 含量在50cm 以上時隨著剖面深度的增加而降低，50cm 以下隨土層深度增加而升高；As 和Ni 變化不明顯。在剖面PS 中，7 種重金屬含量均以表土最高，隨土層深度的增加略有所下降，有可能是自然變化造成的［24］。在剖面GS 中，Ni、Cu、Zn 隨深度增加而升高，As、Hg、Pb隨深度增加而降低；Cd 變化不規則，波動較大，可能受成土母質的影響，在剖面GS，0～40cm 為灰色黏質土，80～90cm 為青灰色淤泥，120～130cm 為頁巖、泥巖；一些重金屬(Cd、Cu、Zn、Pb)會隨著工業活動、化肥/農藥的施用和汽車尾氣排放在表層土壤中逐漸積累［24-26］。本研究中Cd、Cu、Zn、Pb 總含量沿土壤剖面的變化也支持了這一點。

圖3 三個采樣點土壤剖面重金屬垂向分布特征Fig. 3 Vertical distribution characteristics of heavy metals in soil profiles at three sampling points.

3.3 土壤剖面有機碳、養分元素及pH 垂向分布

三個剖面土壤中有機碳、N、P、K2O 和pH 值的垂向分布如圖4 所示。

圖4 三個采樣點土壤剖面有機碳(Corg)、氮(N)、磷(P)、氧化鉀(K2O)含量和pH 的垂向分布特征Fig. 4 Vertical distribution characteristics of organic carbon (Corg), nitrogen (N), phosphorus (P), potassium oxide (K2O) contents and pH in soil profiles at three sampling points.

土壤剖面pH 值為4.48～8.46，在剖面GS 呈酸性，在剖面YS、PS 呈堿性，可能與土壤類型有關。土壤有機碳含量在0.2%～3.37%之間，隨土層深度的增加而降低，在剖面GS 變化最明顯。Cd、Pb 含量的最高值出現在表層，這一發現可能是由于不同層次土壤理化性質的影響。N 含量隨土層深度增加而降低，變化范圍在496.94～2323.94mg/kg。P 含量為227.46～1315.8mg/kg，在剖面YS 變化最明顯。K2O 含量范圍為1.6%～3.55%，在剖面YS、PS 差異不明顯，在剖面GS，80cm 以下含量明顯增加。

3.4 重金屬濃度與土壤養分的相關關系

通過Pearson 相關分析，確定了不同土壤類型三個剖面中重金屬濃度與土壤養分的相關關系(表3)。在剖面YS，Cd、Cu 與N、P、Corg呈極顯著(p＜0.01)正相關；Hg、Zn 與K2O 呈顯著正相關(p＜0.05)，與pH 呈極顯著負相關(p＜0.01)；Pb 與N、P、Corg呈極顯著正相關，與K2O 呈極顯著負相關；As 和Ni 與N、P、K2O、Corg、pH 無相關性。在剖面PS，As、Cd、Cu、Hg、Zn 與N、P、Corg呈正相關(As 和P 除外)，與pH 呈負相關；Ni 與K2O 呈正相關。在剖面GS，Cd與N、P、K2O、Corg、pH 之間的關系不密切，Cu 與P、K2O 呈正相關；As、Hg、Pb 與N、Corg呈極顯著正相關，Zn、Ni 與N、Corg呈負相關。一個潛在的原因可能是土壤有機質對重金屬的吸附，土壤有機質對Cd和Pb 等金屬有很大吸附力［27］。但是，腐植酸和胡敏素［28］可能會降低土壤中某些金屬的含量，因此較高的有機質含量可能不會固定更多的金屬，本研究中在茶園剖面(GS)，Zn 和Ni 含量與土壤有機碳(Corg)呈顯著負相關。氮、磷、鉀作為肥料的標志元素，在農業生產過程中極易在表層土壤中累積?；屎娃r藥不僅會將重金屬直接帶入土壤，還會影響pH、有機質等其他土壤因子，從而增加重金屬含量［29-30］。大量研究表明，磷肥中含有一定量重金屬。磷肥中重金屬含量高低與磷礦及其來源有關。美國磷肥中 Cd 含量在 4～100mg/kg ，摩洛哥磷肥中 Cd含量在 10～24mg/kg［31］。本研究中玉米地Cd 與氮磷呈顯著正相關，剖面YS 相關系數分別為0.813 和0.947；剖面PS 相關系數分別為0.845 和0.747。

表3 三個采樣點土壤剖面重金屬與土壤養分指標的Pearson 相關性Table 3 Pearson correlation between heavy metals in soil profiles and soil nutrient indicators at three sampling points.

一般來說，土壤pH 值與土壤中重金屬含量呈負相關關系，因為pH 值不僅影響土壤顆粒對重金屬的靜電吸附，而且會破壞母質的惰性部分，從而影響土壤中重金屬含量［3，32］。較酸性的環境有利于金屬作為土壤中的不穩定組分加入土壤，也限制了金屬離子與土壤礦物質或有機物的結合。

由表3 可知，在YS 剖面有機碳(Corg)與Cd、Cu、Pb 呈極顯著正相關，與Hg 呈極顯著負相關。在PS剖面Corg與As、Cd、Cu、Hg 呈極顯著正相關，與Zn呈正相關，Cd、Hg 與Corg相關程度非常高，相關系數分別達0.934 和0.955(圖5)；Cd、Hg 含量與土壤pH的相關系數(圖5)分別為-0.964 和-0.944。土壤有機質和pH 值是影響重金屬在土壤中遷移的重要因素。pH 和有機碳兩個變量均能較準確地預測土壤重金屬總含量，兩者的分布可能對該地區土壤重金屬的分布有一定的控制作用，Khaledian 等［33］的數據說明了這一點。

圖5 紫色黏質土(PS 剖面) Cd 與Corg (a)、Hg 與Corg (b)、Cd 與pH 值(c)、Hg 與pH 值(d)相關關系Fig. 5 The correlation between Cd and Corg (a), Hg and Corg (b), Cd and pH (c), Hg and pH (d) in purple clay soil (profile PS).

在剖面GS，有機碳(Corg)與As、Hg、Pb 呈極顯著正相關，與Ni 呈極顯著負相關，與Zn 呈顯著負相關。由于Hg 對土壤有機質官能團具有較高的親和力，土壤有機碳通常與Hg 濃度具有很強相關性［34］，與研究結果一致。Ahmad 等［35］認為有機碳增加了土壤溶液中的重金屬，阻止了重金屬與其他離子的絡合，從而增加了土壤中重金屬濃度?？梢?，有機碳對土壤中重金屬的保留起了重要作用。這也解釋了表層土壤含有較高的Cd、Cu、Pb，因為這部分的有機碳含量相對較高。一般而言，細顆粒土壤比表面積較大，有機質含量較高，因此對重金屬(樣體)離子的結合能力較高［36-37］。研究發現［16］，土壤pH 值在決定Cd 在土壤結合位點上的吸附以及Cd 在土壤溶液中的形態、溶解度和流動性方面起著最重要的作用。

使用Pearson 相關分析可以確定重金屬之間的相關關系，提供其來源和傳輸的信息，如果重金屬之間存在顯著的正相關關系，說明其來源相同，途徑相似［38-39］，各采樣點土壤Pearson 相關分析結果見表4。As 與Cd、Hg、Pb 呈極顯著正相關，Cd 與Cu、Zn、Pb 呈極顯著正相關(p＜0.01)，表明相似的地球化學行為或輸入來源可能與農藥和化肥的使用有關。本研究的結果與以前的一些研究一致。例如，在東北德惠的農業土壤中，Zn 和Cu 也存在較高的相關性［25］，這兩種元素往往受到農業生產人為輸入的控制。As、Hg 與Pb 之間呈較強的正相關，提示它們可能有共同來源。但是Ni、Pb 之間沒有發現統計上的相關性，這可能表明它們的來源不同。

揮發性風味物質中的雙乙酰、乙醛和乳酸對酸奶風味貢獻突出，并對酸奶的氣味和口感起到重要作用。酸奶中β-半乳糖苷酶與乳酸脫氫酶對乳酸的生成影響很大，可以通過調節酶活控制酸奶后酸化問題。β-半乳糖苷酶、α-乙酰乳酸脫羧酶、葡萄糖激酶、乳糖通透酶和丙酮酸激酶的酶活對酸奶中雙乙酰、乙醛的產量有很大相關性。探究酸奶中的酶與揮發性風味物質的關系，有助于酸奶風味物質的形成機理研究。

3.5 土壤污染評價結果

圖6 給出0～10cm、30～40cm、70～80cm、100～110cm 各層重金屬的Igeo值，可以代表這些重金屬在整個垂向剖面上的變化。所有土壤剖面中Cu、Zn的Igeo值均小于0，表明研究區土壤未受到這些重金屬的污染。Cd 在四個深度的Igeo值有明顯降低，除了在剖面GS 的Igeo值小于1，在剖面YS、PS 均大于1，說明玉米地(YS、PS)的Cd 污染程度較茶園地(GS)嚴重，這可能是由于耕作條件的差異，且YS 點表土Igeo值在2～3 之間，表現為中等-強污染。綜上所述，由于種植不同作物不可避免地會改變土壤的物理和化學性質，不同農區土壤重金屬污染水平差異很大。旱地和茶園地土壤的基本理化特征不同，從而影響重金屬的化學轉化和生物有效性［40］。金屬的有效性取決于土壤有機質［41］，土壤有機質是土壤系統中最重要的組成部分之一，因為它對陽離子交換能力有重要影響，陽離子交換能力使土壤具有保留陽離子金屬的能力［33］。特別是在污染環境中，土壤中金屬的數量和有效性可以隨著有機質和總有機碳的增加而增加。As、Hg、Ni、Pb 在四個深度的Igeo值均小于1 且接近于0，表明土壤污染輕微，可能是人為輸入造成的，也可能是自然變化造成的。土壤中金屬的天然來源與土壤母質有關，而人類對土壤金屬的輸入與各種人類活動有關。一般來說，由于過度使用化肥和農藥、廢水灌溉和大氣沉降，常規農業做法導致土壤中重金屬富集。Zhao 等［42］研究發現，化肥和糞肥的使用使重金屬(Cd、Cu、Pb 和Zn)含量增加約3%/年。在很多情況下，金屬污染物首先滯留在土壤表層，然后由不定期的淋濾水(特別是大量的雨水、灌溉水等)攜帶向下運移［36］。底土的污染主要是由表層土壤的金屬垂直運移引起的。同時，土壤中重金屬的遷移能力取決于其與土壤組分的相互作用［32］。

圖6 三個采樣點剖面土壤0～10cm 深度(a)、30～40cm 深度(b)、70～80cm 深度(c)和100～110cm 深度(d)的重金屬地質累積指數( Igeo)Fig. 6 Geo-accumulation indexes ( Igeo ) of heavy metals in profile soils at 0-10cm depth (a), 30-40cm depth (b), 70-80cm depth (c)and 100-110cm depth (d) at the three sampling sites.

重金屬的污染程度由高到低依次為：Cd＞Hg＞As＞Pb＞Ni。以白云巖為主的三疊系雷口坡組(T2l)和以泥巖為主的侏羅系蓬萊鎮組(J3p)中Cd的Igeo值較高，以砂巖為主的三疊系須家河組(T3x)中Cd 的Igeo值較低，說明重金屬的富集與成土母質關系較為密切。陳文軒等［43］基于福建、廣東、海南、浙江、湖南、 陜西、甘肅、河南、重慶、山西、天津、內蒙古以及安徽等地農田土壤重金屬實測數據，發現各行政區農田土壤Cd 與Hg 的地質累積指數較高，與本次研究結果相似。Yang 等［44］分析了中國402個工業用地和1041 個農業用地土壤中重金屬的含量，確定Cd、Pb 和 As 為優先控制的重金屬。各重金屬的Igeo值分析表明，Ni 值的Igeo值最低，污染輕微，而Cd 的Igeo值最高，Hg 次之，污染嚴重?？梢源_定Cd 和 Hg 元素為農田土壤優先控制重金屬。由于Cd、Hg 對人體的毒性，可能對人體健康構成更大的威脅，亟待治理。

3.6 潛在生態風險評估

三個剖面10、30、60、90、110、130、140cm 土層As、Cd、Cu、Hg、Ni、Pb、Zn 的值結果列于表5。除了Cd 和Hg，大部分重金屬(Ei＜40) 的Ei值依次降低：Cd＞Hg＞As＞Pb＞Ni＞Cu＞Zn，說明這些重金屬的污染程度較低。根據風險等級的描述，剖面YS 的Cd、Hg 潛在風險指數最高，在0～140cm 深度存在較大的Cd、Hg 生態風險(80≤Ei＜160)，其中表層土壤Cd 存在很強生態風險(160≤Ei＜320)。表明該地區的鎘污染源可能受以往農業活動的影響，包括施用化肥和農藥。Wei 等［45］研究也表明陜西省石泉縣農田土壤中Cd、Hg 富集程度高，毒性系數高，生態風險高。在剖面PS，Cd 在 0～30cm 土壤表現為很強生態危害(160≤Ei＜320)，在60～110cm表現為強生態危害(80≤Ei＜160)。在剖面GS，表層土壤Cd、Hg 存在中等生態風險(40≤Ei＜80)。RI值在剖面YS 上0～10cm 表現為強生態風險(300≤RI＜600)，在30～140cm 表現為中等生態風險(150≤RI＜300)。在剖面PS，在0～30cm 中等生態危害(150≤RI＜300)，在60～110cm 表現為輕微生態危害(RI＜150)。在剖面GS，在0～130cm 表現為生態危害相對較弱。土壤各重金屬Ei值隨深度的增加而降低，與Igeo污染評價結果一致。三個剖面Cd 的Ei值相對較高，表明在控制重金屬污染方面需要特別關注。張小敏等［46］對中國農田土壤重金屬分布特征研究表明，南方地區由于工業活動和礦區開采，導致土壤Cd 累積明顯強于北方地區。綜上所述，可以發現三個剖面表層土壤生態危害指數(RI)較高，且Cd 的貢獻最大。因此，為控制和緩解土壤重金屬污染，今后應進一步加強對沐川縣每個鄉鎮進行土地質量改善和監測，規范和調整工業和農業活動。

表5 三個采樣點土壤剖面重金屬潛在生態風險指數Table 5 Potential ecological risk index of heavy metals in soil profiles of three sampling points.

4 結論

選擇四川省沐川縣采集土壤剖面樣品，采用電感耦合等離子體質譜法、氫化物發生原子熒光光譜法、高頻燃燒紅外吸收法等方法測定了樣品中7 種重金屬元素，土壤養分元素氮、磷、鉀，有機碳含量和pH 值，結合地質累積指數法和潛在生態風險指數法，對研究區的土壤剖面重金屬含量、生態風險特征進行了分析。根據土壤垂直剖面數據結果可以得出，重金屬含量趨于在表層土壤富集，隨深度增加其含量降低。Igeo值與Ei值亦隨地層深度的增加而降低。研究區各指標的相關關系分析結果顯示，土壤重金屬含量與有機碳呈顯著正相關，其原因應與有機質吸附作用有關。土壤重金屬含量與pH 值呈顯著負相關，因為pH 不僅影響土壤顆粒對重金屬的靜電吸附，而且會破壞母質的惰性部分，從而影響土壤中重金屬含量。Igeo分析表明，旱地土壤(種植玉米，剖面YS)受到Cd 的污染，且旱地(剖面YS)的污染程度比旱地(剖面PS)和茶園地(剖面GS)嚴重。因為旱地(剖面YS)位于三疊系雷口坡組(T2l)，巖性為粉砂巖與白云巖、泥質灰巖互層，夾黑色碳質頁巖。表明研究區內土壤中重金屬含量及生態風險與地質背景密切相關。潛在生態風險分析表明，Cd、Hg 對玉米地土壤(剖面YS、PS)均具有較大的生態風險，對茶園地(剖面GS)生態風險相對較弱。

通過本次研究，較為客觀地得出四川省沐川縣土壤剖面重金屬的累積和潛在生態風險狀況，以及與土壤養分元素的關系，能夠為研究區重金屬污染防范提供科學支撐。耕地重金屬負荷量的增加勢必會造成農產品中重金屬含量的增加，嚴重地將導致農產品中的重金屬超標問題，因此要加強關注與當地土壤重金屬直接相關的作物中的重金屬含量。重金屬、pH 和土壤養分元素之間的多元統計分析有待進一步研究。針對研究結果，建議深入研究重金屬在土壤中的累積特征，合理評估其生態危害，保證土地的安全利用。

致謝：特別感謝中國地質科學院地球物理地球化學勘查研究所孫彬彬、候樹軍在樣品采集中付出的辛勤勞動。

Distribution Characteristics and Ecological Risk Assessment of Heavy Metals in Typical Soil Profiles of Muchuan County, Sichuan Province

LIAO Bixia1,2，SHEN Wenling1,2，HE Ling1*
（1. Institute of Geophysical and Geochemical Exploration, Chinese Academy of Geological Sciences, Langfang 065000, China；
2. Guilin University of Technology, Guilin 541006, China）

(1) Most of the heavy metal content in the profile is higher than the national and Sichuan background values,indicating a higher geological background.

(2) The heavy metal content is closely related to soil nutrients and physicochemical properties, positively correlated with organic carbon content, and negatively correlated with pH value.

(3) The accumulation of Cd and Hg is relatively high, and there is a certain potential ecological risk, while the ecological risks of other elements can be ignored.

ABSTRACTBACKGROUND:Soil is a precious resource for human survival and social development. The quality of the soil environment is impacted by a variety of issues due to the social economy’s rapid expansion, and the issue of heavy metal contamination in farmed land has garnered great attention globally. Heavy metals in soil pose a severe risk to the security of agricultural products and public health due to their persistence, latency, and ease of entry into the food chain.In recent years, many scholars have carried out research on soil heavy metal pollution and ecological risk assessment under different conditions such as natural conditions, industrial and mining industries and developed transportation in different regions. Zhou et al.[11]found that Xiong’an New Area was affected by the production activities of surrounding enterprises. The contents of As, Cd, Cu, Pb and Zn in some root soil samples exceeded the screening value standard for soil pollution risk of agricultural land (GB 15618—2018), and the exceeding ratios were 23.33%, 96.67%, 33.33%, 33.33% and 10.00%, respectively. Song et al.[12]evaluated the characteristics of heavy metal pollution in the surface soil of Fuping County, Hebei Province, and found that As and Cd exceeded the acceptable carcinogenic risk level (As is 10-5, Cd is 10-6). Kumar et al.[10]collected data on heavy metalcontaminated soils in India from 1991 to 2018. The average Cd content of all soil types exceeded the limit values,and the potential ecological risk values of Cd were greater than 320, reflecting a higher ecological risk. For the heavily polluted soil, according to the different pollution situation in our country, the remediation measures are taken according to local conditions. However, due to the wide area of contaminated soil and the complex composition of pollution sources, the current soil remediation work still faces huge problems.OBJECTIVES:To study the vertical distribution characteristics of heavy metals in soil, the relationship between soil heavy metals and soil nutrient elements, as well as the degree of pollution and potential ecological risks.METHODS:The contents of Cd, Cu, Ni, Pb, Zn were measured using inductively coupled plasma-mass spectrometry (ICP-MS); As content was determined by hydride generation atomic fluorescence spectrometry (HGAFS); P and K2O contents were determined by X-ray fluorescence spectrometry (XRF); N content was determined by oxidation combustion gas chromatography (GC); Hg content was determined by cold vapor atomic fluorescence spectrometry (CV-AFS); Organic carbon content was determined by high-frequency combustion infrared absorption method (IR); potentiometric method (POT) was used to measure soil pH value. Statistical analysis and calculation of soil heavy metal content, pollution index, and ecological risk index were conducted using Excel 2016. Pearson correlation analysis was conducted using SPSS 26, and the degree of soil heavy metal pollution was evaluated using the geoaccumulation index (Igeo). Potential ecological risk index (RI) values were selected to evaluate potential ecological risks.RESULTS:The average contents of As, Cd, Cu, Hg, Ni, Pb, and Zn in the soil of YS plot were 20.8mg/kg,0.35mg/kg, 26.38mg/kg, 0.121mg/kg, 33.29mg/kg, 42.37mg/kg, and 94.47mg/kg, respectively; The average contents of As, Cd, Cu, Hg, Ni, Pb, and Zn in the soil of PS plot were 7.21mg/kg, 0.32mg/kg, 28.32mg/kg,0.028mg/kg, 47.34mg/kg, 33.29mg/kg, and 116.45mg/kg, respectively; The average contents of As, Cd, Cu, Hg, Ni,Pb, and Zn in the soil of GS plot were 5.42mg/kg, 0.16mg/kg, 22.38mg/kg, 0.08mg/kg, 31.8mg/kg, 30mg/kg, and 75.03mg/kg, respectively. The concentrations of As, Cd, Hg, Ni, Pb, and Zn were higher than the national and Sichuan soil background values, indicating that these metals were relatively enriched in the soil of Muchuan County.The relationship between seven heavy metals at different soil depths was evaluated through Pearson correlation analysis (seen in Table 4). There was a significant positive correlation between heavy metals, indicating their widespread homology. In the PS profile, the correlation between Cd, Hg and organic carbon was very high, with correlation coefficients of 0.934 and 0.955, respectively (Fig.5); As, Cd, Cu, Hg, Zn showed a highly significant negative correlation with pH, and the correlation between Cd, Hg content and soil pH was shown in Fig.5, with correlation coefficients of -0.964 and -0.944, respectively. The content of heavy metals in soil was closely related to organic carbon and pH value, which should be attributed to the adsorption of organic matter and the fact that pH not only affected the electrostatic adsorption of heavy metals by soil particles, but also damaged the inert part of the parent material. Soil organic matter and pH value are important factors affecting the migration of heavy metals in soil. The surface soil had a high content of organic matter, multiple adsorption sites, and a high soil pH value, which reduced the solubility of heavy metals and thus the metal migration rate.Soil pollution assessment results.TheIgeovalues of Cu and Zn in all soil profiles were less than 0, indicating that the soil in the study area was not contaminated by these heavy metals. TheIgeovalue of Cd at four depths was significantly reduced. Except that theIgeovalue at GS point was less than 1, YS and PS were greater than 1,indicating that the Cd pollution degree of corn land (YS, PS) was more serious than that of tea garden land (GS).This may be due to the difference of tillage conditions, and theIgeovalue of surface soil at YS point was between 2 and 3, showing moderate-strong pollution. TheIgeovalues of As, Hg, Ni and Pb at four depths were all less than 1 and close to 0, indicating that the soil pollution was slight, which may be caused by human input or natural changes.In general, conventional agricultural practices lead to the enrichment of heavy metals in soils due to excessive use of fertilizers and pesticides, wastewater irrigation and atmospheric deposition. Zhao et al.[42]found that use of fertilizers and manure increased the content of heavy metals (Cd, Cu, Pb, and Zn) by approximately 3% per year.The order of heavy metal pollution degree from high to low is Cd>Hg>As>Pb>Ni.Potential ecological risk assessment.According to the description of risk level, the YS plot had the highest potential risk index for Cd and Hg, and there was a significant ecological risk of Cd and Hg at depths of 0-140cm(80≤Ei＜160), among which the surface soil Cd had a strong ecological risk (160≤Ei＜320). It indicates that Cd pollution sources in the region may be affected by past agricultural activities, including fertilizers and pesticides.The soil Cd of PS plot exhibited strong ecological hazards (160≤Ei＜320) at the depth of 0-30cm while exhibiting strong ecological hazards (80≤Ei＜160) at 60-110cm. The Cd and Hg in surface soil at the GS plot site had moderate ecological risks (40≤Ei＜80). The value ofRIshowed a strong ecological risk (300≤RI＜600) at 0-10cm of the YS plot, and a moderate ecological risk (150≤RI＜300) at 30-140cm. Moderate ecological hazards (150≤RI＜300) were present in the PS plot, while mild ecological hazards (RI＜150) were present at 60-110cm. The ecological hazards of GS plot at 0-130cm were relatively weak. TheEivalues of heavy metals in soil decreased with the increase of depth,which was consistent with the evaluation results ofIgeopollution. TheEivalues of Cd in the three profiles were relatively high, indicating that special attention should be paid to the control of heavy metal pollution.CONCLUSIONS:According to the results of soil vertical profile data, it can be concluded that heavy metal content tends to accumulate in the surface soil, and its content decreases with increasing depth. TheIgeovalue andEivalue also decrease with the increase of formation depth. The geoaccumulation index and potential ecological risk analysis indicate that Cd poses significant ecological risks to the local soil, and appropriate measures should be taken to strengthen pollution prevention and control in the area to avoid harm to human health. The content of heavy metals is closely related to soil nutrients and physicochemical properties, positively correlated with organic carbon content,and negatively correlated with pH value. According to the research results, it is suggested to carry out further research on the accumulation of heavy metals in soil, rationally assess its ecological harm, and ensure the safe use of land.

KEY WORDS：soil；heavy metals；content distribution；pollution assessment；inductively coupled plasma-mass spectrometry