李 樋,李隨民,王 軼,劉創脫
1. 河北地質大學, 河北 石家莊 050031;2. 中國地質環境監測院,北京 100081
隨著當前經濟社會的迅速發展,土壤重金屬污染的問題已經得到了社會各界廣泛關注,土壤作為生態系統的重要組成部分,與人類的生命活動密切相關[1-2]。土壤環境背景值作為判斷土壤是否受到重金屬污染的標準目前在很多的研究中得到了廣泛的應用。土壤環境背景值的確定需要找到未受到或很少受到人類行為活動影響的土壤,然而,在地球表面幾乎所有地方都已經受到了人類行為活動的影響,很難找到未受環境影響的土壤,所以國外就有學者提出了地球化學基線的概念。
地球化學基線(Geochemical baseline)最早出現自國際地球化學填圖項目(IGCP259)和全球地球化學基線項目(IGCP360)中,這一概念一直未被很好的定義。在國際地球化學填圖項目中,地球化學基線被定義為地球表層物質中化學物質(元素)濃度的自然變化[3],隨著對地球化學基線研究的深入,其定義也進一步得到明確:地球化學基線將某一地區或數據集合作為參照某一元素在特定物質中(土壤、沉積物、巖石)的自然豐度,并可以表述為區分地球化學背景和異常的單一的基線[4]。對于該問題,雖然國際上尚未形成具有共識的明確定義,但在大尺度上已經對土壤環境地球化學基線的研究方法取得共識并形成了初步的標準。目前,在中、小尺度上確定研究區基線值并準確實現其空間表達仍需深入研究,以便提高陸表環境過程與變化研究的可靠性,同時有利于進行精細格網上的有效調控與管理。
地球化學基線的確定的方法有標準化法、傳統統計法、地球化學對比法,傳統統計法又包含相對累積頻率法、平均值法、盒須法等。研究將選取標準化法[5-6]和累積頻率法[7-8]兩種常用的地球化學基線確定方法。
研究區位于內蒙古自治區通遼市開魯縣東來鎮,西遼河南岸,總面積100 km2。大地構造位置屬華北板塊北部大陸邊緣[9],地層從太古界到新生界在研究區內均有出露。研究區全部被第四系覆蓋(圖1),第四系土壤是研究的主要對象。
研究區交通十分便利,京通鐵路和國道111線穿過全境,鐵路設有貨運站臺。國道111線由東向西經通遼市與遼寧省相接。
研究區屬半干旱大陸季風氣候區,四季分明,春季多風干旱,夏季雨熱同期,秋季溫涼少雨,冬季寒冷漫長。主要農作物為玉米,為一年一作,除玉米外,研究區還分布有零星花生、大豆、向日葵等農作物種植區。

圖1 研究區區域地質圖Fig.1 Regional geological map of the study area
以內蒙古東來地區為研究區域,嚴格按照《多目標區域地球化學調查規范》的要求進行采樣:表層樣采集時除去地表植物凋落物,采集地表至地下20 cm深處的連續土柱,一共采集表層土壤樣品468件(圖2)。樣品的采集方法為刻槽取樣法,以保證土壤樣品上下的均勻性。每處采樣點要在部署采樣點周邊30~40 m范圍內另外采集2處同類型土壤樣品組合成1件樣品,樣品的記錄和GPS定點均以中間取樣點為準。樣品原始總重量不低于1 kg,以保證干樣重量大于800 g,自然風干后,用木棒壓碎,去除雜物,以確保自然干燥過20目篩后,獲取的粒級樣品重量過500 g。
樣品送至天津市地質礦產測試中心實驗室進一步的處理,測試項包括重金屬元素在內共計23項指標,樣品的檢測方法和檢出限如下(表1),Cu、Ni、Zn、Fe測試方法為X熒光光譜法,Cd、Pb測試方法為等離子質譜法。

圖2 研究區表層土壤采樣點位圖Fig.2 Bitmap of surface soil sampling points in the study area

表1 重金屬元素檢出限及分析方法Table 1 Detection limit and analysis method of heavy metal elements
注:ICP-MS:等離子質譜法;XRF:X熒光光譜法。
將獲得的實驗數據在Spss 20.0里進行描述性統計分析,結果顯示:Cu、Cd、Ni、Pb、Zn含量范圍分別為6.4~38.8 mg/kg、0.025~0.344mg/kg、6.3~45.5 mg/kg、1.3~31.4 mg/kg、11.3~103.6 mg/kg(表2),參照開魯土壤背景值[9-10]和GB15618-2018,Cu、Cd、Ni、Pb、Zn 5種重金屬元素含量的平均值均高于背景值,依次是背景值的1.92倍、2.07倍、1.49倍、1.18倍、2.09倍,5種重金屬元素含量的最大值低于國家標準值。變異系數反映了元素數據間的離散程度和樣本的空間變異性大小,一般認為,變異系數<10%為弱變異,10%~100%中等強度變異,變異系數≥100%為強變異[11-12]。Cu、Cd、Ni、Zn、Pb5種重金屬元素的變異系數分別為33.2%、33.6%、30.6%、29.8%、30.3%,屬于中等強度變異,空間變異相對顯著,可能與人類的行為活動有關。Fe元素的變異系數為29%,小于5種重金屬元素的變異系數,說明Fe元素受人為擾動較小,較穩定。

表2 土壤指標描述統計特征及相關標準Table 2 Descriptive statistics and related standards of soil index
注:元素單位均為mg/kg;*,Ni元素的土壤背景值引自(周妍姿等, 2015);國家標準采用《中華人民共和國土壤環境質量標準GB15618—2018(試行)》。
標準化法的原理是選取地球化學過程中的惰性元素作為標準,用活性元素(污染元素)與惰性元素的相關性來判斷活性元素的富集情況,通過建立二者之間的線型回歸方程[13-14],即基線模型如下:
Cm=aCN+b
(1)
式(1)中Cm和CN分別表示樣品中活性(污染)元素和標準元素的實測濃度;a、b表示回歸系數。
通過對數據的統計分析及數據處理,可以得到回歸方程的系數a、b的值,根據惰性(標準)元素的平均值,可以得到研究區活性(污染)元素平均預測值,即活性元素的地球化學基線值B,模型如下:
(2)
式(2)中Bm表示活性元素m的基線值,CN表示研究區標準元素含量的平均值。
采用標準化方法確定地球化學基線時,首先對標準化元素進行選取及確定,在人為引起的重金屬輸入量低于自然來源時,用Fe元素作為標準化元素也是很好的選擇[15-16],研究區Fe元素較穩定,變異系數為0.29(表2),受人為擾動較小,因此標準化元素選擇Fe元素。

表3 東來地區表層土壤的地球化學基線模型及各元素基線Table 3 Geochemical baseline models and element baselines of surface soils inDonglai area
研究區土壤中Fe元素與Cu、Cd、Ni、Pb、Zn元素具有較好的相關性(圖3)。因此,可以由公式(1)(2)得出Cu、Cd、Ni、Pb、Zn的模型方程及基線值(表3),確定的Cu、Cd、Ni、Zn、Pb五種重金屬元素地球化學基線值依次為21.689 mg/kg、0.143 mg/kg、25.784 mg/kg、17.795 mg/kg、64.425 mg/kg。

圖3 土壤中Fe與5種重金屬元素相關關系圖Fig.3 Relative map of Fe and 5 heavy metal elements in Soil
地質統計分析在地球化學異常和背景的分離研究中被廣泛應用,為地球化學家提供了強大的數據分析和處理工具。環境地球化學研究中不僅需要分離背景和異常,更為重要的是從自然異常中判別人類活動引起的異常。如果將人為活動引起的異常地球化學數據中分離出來,那么地球化學基線自然就被建立了[17]。統計分析不僅可以確定環境地球化學基線,還可以對環境污染進行分析和判別污染的來源[18]。
作為統計分析方法之一的相對累積頻率法采用正常的十進制坐標[19],累積頻率-元素濃度的分布曲線可能有兩個拐點,值較低的點可能代表了元素濃度的上限(基線范圍),小于樣品元素濃度的平均值或中值可以作為基線值;值較高的點可能代表了異常的下限;而二者之間的部分可能與人類活動有關,也可能無關[20]。若分布曲線近似呈直線,則所測樣品的濃度可能本身就代表了背景范圍(基線)[21]。
東來地區土壤樣品中的五種重金屬元素Cu、Zn、Ni、Pb、Cd的相對累積頻率散點圖見(圖4),所有元素的相對累積頻率可以分為三部分,小于70%的部分,代表了樣品的基線濃度范圍;大于70%而小于90%的部分,表示該部分既可能遭到污染,也可能沒有汚染;大于90%的部分,該部分代表的是受到人為擾動的元素的濃度[21]。

圖4 東來地區5種重金屬元素相對累積頻率圖Fig.4 Relative accumulative frequency maps of five heavy metal elements in Donglai area
采用相對累積頻率方法確定的東來地區5種重金屬元素的基線值依次為20.31 mg/kg、0.143 mg/kg、25.64 mg/kg、17.95 mg/kg、66.45 mg/kg(表4)。

表4 相對累積頻率法得到的東來地區基線值Table 4 Baseline values ofDonglai area obtained by relative cumulative frequency method
采用標準化方法和相對累積頻率分析方法確定的合肥地區的土壤環境地球化學基線值,可以看出兩種法確定的基線值相近(表5),說明選取的兩種確定研究區地球化學基線方法是可行的。因此,在后續的研究中,將以這兩種方法的均值作為研究區的地球化學基線值。

表5 東來地區的土壤環境地球化學基線值Table 5 Baseline values of Soil Environmental Geochemistry in Donglai area
基線值與背景值有著不同的概念和意義。基線值表示自然狀態下土壤中化學成分的濃度,包含人類的活動對自然影響的程度,背景值表示未受或很少受到人類行為活動污染的環境要素的平均化學成分,是圈定異常大小的一般含量及一般變化幅度[22]。
東來地區Cu、Cd、Ni、Pb、Zn 5種重金屬元素的地球化學基線值均低于國家標準值(表6),說明研究區的土壤整體未受到重金屬元素的污染。但環境地球化學值均高于開魯縣土壤的背景值,分別為土壤背景值的1.86、2.07、1.49、1.18、2.19倍,高出土壤背景值的部分可能來自人類活動的累積,說明人類活動已經對研究區的土壤質量產生了一定影響。

表6 東來地區土壤中重金屬的基線值及評價標準Table 6 Baseline values and evaluation criteria of heavy metals in soils ofDonglai area
通過對內蒙古東來地區土壤中Cu、Cd、Ni、Pb、Zn共5種重金屬元素采用標準化方法和相對累積頻率方法進行分析研究得出:
(1)內蒙古東來地區土壤中Cu、Cd、Ni、Pb、Zn共5種重金屬元素的地球化學基線分別為20.995 mg/kg、0.143 mg/kg、25.712 mg/kg、17.873 mg/kg和65.438 mg/kg。
(2)與研究區土壤背景值相比,Cu、Cd、Ni、Pb、Zn5種重金屬元素含量分別為背景值的1.86、2.07、1.49、1.18、2.19倍。推斷高出土壤背景值的部分主要來自人類活動造成的累積。
(3)確定研究區土壤環境地球化學基線值,對提高研究陸表環境過程與變化研究的可靠性以及對進行精細格網上的有效調控與管理提供了科學依據。