1. 引言
贵州遵义地区,是贵州省农业生产重要基地,针对开展无公害农产品产地现状的评价,需建立土壤环境背景值的调查,这是研究保护指导土地利用的重要指标。在无公害农业产地评价如何掌握、确立土壤重金属环境背景值的计算 [1],要涉及到土壤环境背景值 [2] 的计算。我国碳酸盐岩地区自然过程是造就镉等元素超标的主要原因之一,导致重金属元素的背景值不一致,如表层镉含量范围值(0.015~2.977 mg/kg) [3]。贵州主要农业耕作集中区,需要应用土壤环境背景值的方法对产地进行评价,因而提出应用环境地球化学基线(environmental geochemical baseline),能描述地球表层物质中化学物质(元素)浓度的自然变化 [4],可以用于该地区土壤基线与背景范围,为农业无公害环境评价提供一个重要参考依据 [5]。环境地球化学基线是建立在当前地球表层环境的状态档案,并提供监测环境变化的数据库,研究矿物及化学元素的自然变化,因此用环境地球化学基线,既能反映当前农业生产区背景,又能提供监测变化的数据库和自然空间变异情况 [6]。
本文通过对贵州遵义地区农业表层889个土壤中Hg、Cd、Pb、Cr、As的调查与分布,通过环境地球化学基线分布以及地质累积指数评价与生态风险指数评价,为建立区域尺度土壤评价中环境质量的评价提供数据与信息。
2. 材料和方法
2.1. 研究区域概况
研究区域遵义地区位于贵州北部,地处:东经106˚17'22~107˚25'25,北纬27˚13'15~28˚04'09。国土总面积4092.66 km2,其中耕地面积1212 km2,占总面积的30%。土壤类型以黄壤土类、水稻土类、紫色土类为主。遵义地处大娄山山脉南侧,属黔中丘原和黔北山地过渡地带。属亚热带季风气候,平均有霜期为9.5天,无霜期长,雨量充沛,年平均降雨1035毫米。
遵义地区是全省商品粮、油菜籽、生猪、烤烟生产基地县,“全国无公害农产品示范基地县”,全县城镇化水平高于全省平均水平。
2.2. 采样布点
根据重金属分布的空间变异性和城郊土壤空间差异十分明显的特点,采用非均匀性布点方法,共采集样本889个地理位置图见图1。
![](//html.hanspub.org/file/1-1750821x10_hanspub.png)
Figure 1. Distribution of sample collection map
图1. 样品采集区分布示意图
采样点选在被采土壤类型特征明显,地形相对平坦、稳定、植被良好的地点。人为干扰大,失去土壤的代表性,不设采样点。采样点以剖面发育完整、层次较清楚,不在水土流失严重或表土被破坏处设采样点,不在多种土类、多种母质母岩交错分布、面积较小的边缘地区布设采样点。在100亩耕地内随机确定一个10 m × 10 m正方形作为取样点。该样点的土壤样品从10 m × 10 m正方形4个顶点和中心共5处各采集1 kg表土(0~20 cm),均匀混合后,用四分法选取1 kg作为该点样品,同时用GPS测定正方形中心的地理坐标并编号记录。
2.3. 样品分析
土壤样品采用美国国家环保局相关标准(USE-PA)抽提消煮 [7],用原子荧光测定砷,用冷原子吸收法测定Hg、As,石墨炉原子吸收光谱仪(AAS Vario 6)测定Cd、Pb、Cr。分析过程加入国家标准土壤样品(GSS-2、GSS-5)进行分析质量控制,分析样品重复数10%~15%。
2.4. 数据分析
2.4.1. 统计分析
采用Excel和SPSS数据处理软件对所测的数据进行统计分析。
2.4.2. 环境地球化学基线分析
1) 相对累积总量分析方法
地质表层农业土壤中重金属元素的浓度值呈对数正态分布。在相对累积密度与元素浓度的双对数分布图中,分布曲线拐点处元素的浓度值就是该元素背景与异常的分界线。在小于分界点的元素浓度数据的平均值加2倍标准差的控制线,通常是元素的背景值范围 [8]。
2) 相对累积频率分析方法
该方法采用十进制坐标,而累积频率与元素浓度的分布曲线则有两个拐点。含量较低的拐点(称拐点1)代表元素浓度的上限(基线范围),小于样品元素浓度的平均值或中值即可以作为基线值。含量较高的拐点(称为拐点2)则代表异常的下限(人类活动影响的部分)。这二者之间的部分可能与人类活动有关,也可能无关。若分布曲线近似呈直线,则所测样品的浓度本身就代表了背景范围(基线) [9] [10] [11] [12] [13]。
2.4.3. 地质累积指数
地质累积指数(Geoaccumulation Index, Igeo)通常称为Muller指数 [14] [15]。是20世纪60年代晚期在欧洲发展起来的广泛应用于研究沉积物中重金属元素污染程度的定量指标 [16],其表达式如下:
(1)
修正指数,通常用来表征沉积特征、岩石地质及其他影响。地质累积指数可分为7个级别 [16],不同的级别分别代表不同的重金属元素污染程度(表1)。
![](Images/Table_Tmp.jpg)
Table 1. The contamination degree based on geoaccumulation index
表1. 地质累积指数不同级别代表的重金属元素污染程度
2.4.4. 潜在生态风险评价方法
采用瑞典科学家Hakanson提出的潜在生态风险指数法 [17] (Hakanson, 1980),根据Hakanson法,研究区域内土壤中第i种重金属的单项生态危害系数(The average single ecological risk index) Ei及土壤中多种重金属的综合潜在生态危害指数 [17] (The complex ecological risk index) Ri可分别表示为:
单项生态危害系数:
(2)
复合生态风险指数:
(3)
式中,Pi为单因子污染指数,Ti为毒性响应系数。
因子Ti为重金属i的毒性响应系数,它主要反映重金属的毒性水平和环境对重金属污染的敏感程度。在本研究中,采用Hakanson所制定的标准化重金属毒性响应系数为评价依据。相应重金属的生物毒性响应系数为:Hg (40);Cd (30);As (10);Pb (5);Cr (2)。根据Hakanson法,土壤重金属生态危害程度的划分标准见表2。
![](Images/Table_Tmp.jpg)
Table 2. Relationship between degree of pollution and potential ecological risk assessment
表2. 潜在生态危害指标与污染程度关系
3. 结果与讨论
3.1. 土壤重金属含量分析
土壤重金属含量监测结果统计分析
遵义地区农业土壤中Hg、Cd、Pb、Cr、As土壤重金属浓度分布图见图2,从图中可见,Cd与As呈明显的非均匀分布特点,而各金属元素之间地理位置无明显相关。按Sturges提出的经验公式[K = 1 + (LgN/lg2)]来确定分布组数 [4],本次采样数为889个,Hg、Cd、Pb、Cr、As分别分为11组进行分布频数统计,结果如图3,可以看出,土壤中Hg、Cd、Pb、Cr、As含量均不符合正态分布,经过对数转换后均符合正态分布。因此,五种重金属不能用算术平均值和算术标准差来表示其含量,要用几何平均值和几何标准差来表示其含量。
![](//html.hanspub.org/file/1-1750821x19_hanspub.png)
![](//html.hanspub.org/file/1-1750821x20_hanspub.png)
Figure 3. Distribution diagram of heavy metal content and logarithmic conversion frequency in geological surface agricultural soil
图3. 地质表层农业土壤重金属含量与其对数转换频数分布图
土壤中五种重金属含量均符合对数正态分布,正态分布图见图3,这是开展统计分析的必要前提条件。经异常值检验处理后土壤重金属的描述性统计见表3,结果表明,用平均值加减三倍标准差替换异常值后 [18] [19],五种重金属的含量分布在峰度与偏度上显著下降,平均值和中位数更为接近,标准差和变异系数也明显减小。从表3可以看出,异常值检验后Cr已经符合正态分布,Hg、Cd、Pb、As含量呈偏态分布,经对数转化以后呈现正态分布更好,在研究的五种重金属中,Pb的变异系数最大,为0.72,其次是Hg、Cd分别为0.66、0.61,Cr、As相对要小为0.36、0.32。
各重金属元素在各乡镇的分布情况,依据表3,对遵义地区土壤中各元素最大值及最小值分布情况做描述表4,对该县采样各镇土壤中重金属平均含量排序及平均含量对比情况分别作描述如表4和图4。
从表4可见,遵义地区Hg含量范围为0.02~2.76 mg/kg,平均为0.182 mg/kg;Cd的含量范围为0.099~2.07 mg/kg,平均含量为0.356 mg/kg;Pb含量范围在5.401~224.64 mg/kg之间,平均34.15 mg/kg;Cr含量范围为8.25~140.27 mg/kg,平均值50.58 mg/kg;As含量范围为2.05~70.56 mg/kg,平均值为10.92 mg/kg。
![](Images/Table_Tmp.jpg)
Table 3. Descriptive statistics of heavy metal content in surface soil in Zunyi area after outlier test
表3. 异常值检验后遵义地区土壤重金属含量的描述统计
![](Images/Table_Tmp.jpg)
Table 4. The maximum and minimum of heavy metal content of surface soil in Zunyi area
表4. 遵义地区土壤中各元素最大值及最小值分布表
![](//html.hanspub.org/file/1-1750821x21_hanspub.png)
Figure 4. Double logarithmic distribution curve of heavy metals concentration and relative cumulative density in geological surface agricultural soil
图4. 地质表层农业土壤重金属元素浓度与相对累积密度的双对数分布曲线
3.2. 遵义地区地质表层土壤重金属环境地球化学基线的确定
3.2.1. 表层农业土壤重金属的基线范围
图4为农业土壤重金属元素浓度与相对累积密度的双对数分布曲线图,途中拐点(黑点所示)对应是元素浓度值,也就是该元素背景与异常的分界值,在小于分界点的元素浓度数据的平均值加2倍标准差是该元素的控制线范围,也是该区域表层土壤的基线值范围(见表5)。
![](Images/Table_Tmp.jpg)
Table 5. Baseline range of heavy metals in surface agricultural soil
表5. 表层农业土壤重金属的基线范围(mg/kg)
3.2.2. 确定地质表层农业重金属的基线值
图5为地质表层农业土壤重金属元素含量概率函数曲线图,可见分布曲线存在2个值,结合基线范围(图5)拐点1代表元素浓度的上限(基线范围),样品中小于此值的元素浓度的平均值或中值可作为基线值,图5中拐点2代表异常值下限(人为活动的影响部分),因此选择此值作为地值表层农业土壤中该元素的基线值,计算出背景与异常的分界值(见表6)。
![](//html.hanspub.org/file/1-1750821x22_hanspub.png)
Figure 5. Curve of probability function of heavy metals content in agricultural surface soil
图5. 地质表层农业土壤重金属元素含量概率函数曲线图
![](Images/Table_Tmp.jpg)
Table 6. Baseline of heavy metals in agricultural surface soil
表6. 表层农业土壤重金属的基线值(mg/kg)
由表6可见,Cd的背景与异常分界值(0.299 mg/kg)与《农业环境监测实用手册》 [20] 中提供的贵州表层土壤Cd的数据(0.659 mg/kg)差异较大,难以为无公害评价背景值所应用。由表5可见,表层土壤重金属基线的平均值Cd (0.227 mg/kg),全国各省(市、区)土壤表层镉样品数据均属于近似对数正态分布,因此以几何均值和几何标准差 [21] 计算其结果,全国各省土壤表层Cd背景含量范围值 [3] 为0.2086 mg/kg,结果相符。
3.3. 地质累积污染分析
1) Hg的地质累积污染分析
Hg的基线值(0.090 mg/kg)和Igeo计算公式,获得Hg的地质累积指数为0.431。对表层土壤的样品进行分析,发现24.09%的样品数的地质累积指数Igeo为l级,60.74%的表层土壤样品数为2级,12.43%的样品数为3级,2.02%的样品数为4级,0.502%的样品数为5级,0.22%的样品数为6级。
2) Cd的地质累积污染分析
通过Cd的基线值(0.147 mg/kg)和Igeo计算公式,获得遵义县表层农业土壤Cd地质累积指数为0.691。对采集的样品分析发现,13.84%样品数的Cd地质累积指数为l级,59.17%的样品数为2级,21.37%的样品数为3级,5.17%的样品数为4级,0.45%为5级。
![](Images/Table_Tmp.jpg)
Table 7. Geoaccumulation index of heavy metals content in each countryside in Zunyi area
表7. 遵义县各乡镇土壤重金属地质累积污染指数
3) Pb、Cr、As的地质累积污染分析
Pb、Cr、As的基线值(分别为20.39 mg/kg、42.78 mg/kg和6.54 mg/kg),采用地质累积污染指数Igeo的计算方法,获得Pb、Cr、As的地质累积污染指数分别为0.159、−0.343和0.155 (表7)。结合表7,发现该地区表层农业土壤Cr元素处于无污染状态,Pb、As处于无污染到中度污染状态。
对于Pb,分析该地区所有样品发现,46.42%样品数的地质累积指数为l级,43.97%的样品数为2级,8.59%的样品数为3级,1.02%的样品数为4级。对于Cr,74.58%样品数的地质累积指数为l级,24.41%的样品数为2级,1.01%的样品数为3级。对于As,40.16%样品数的地质累积指数为l级,47.13%的样品数为2级,10.91%的样品数为3级,1.8%的样品数为4级。Pb、Cr、As这3种元素均有90%以上样品数,为1级和2级,与潜在生态危害对遵义县表层农业土壤中Pb、Cr、As的污染判别情况的判识大致相同。
4) 单项生态危害系数分析
单项生态危害系数Ei排序为Cd > Hg > As > Pb > Cr,可以看出土壤中较高水平的Cd、Hg在构成表层土壤环境污染的同时,也构成了相应的生态危害。这揭示了重金属元素含量高的地区伴随着高的生态风险。该地区Hg、Cd、Pb、Cr、As的地质累积污染指数分别为0.431、0.691、0.159、−0.343和0.155,表明该地区整体表层农业土壤Hg、Cd、Pb、As处于轻度污染状态,Cr无污染。
3.4. 潜在生态危害结果分析
根据表8可见,Cd的单项生态危害系数较高,中等达28.23%,该地区Pb、Cr、As的单项生态危害系数均较低,对生态造成的影响较小。
![](Images/Table_Tmp.jpg)
Table 8. Potential ecological risk index evaluation (%)
表8. 潜在生态风险指数评价(%)
根据图6遵义县表层农业土壤重金属元素的复合生态风险指数均小于150,处于轻微生态危害级别。各个采样乡镇的复合生态风险指数均小于150,处于轻微生态危害级别。各个采样乡镇的复合生态风险指数排序:茅栗镇 > 马蹄镇 > 南白镇 > 尚嵇镇 > 新民镇 > 苟江镇 > 三合镇 > 石板镇 > 乐山镇 > 乌江镇 > 枫香镇 > 三岔镇 > 鸭溪镇 > 龙坑镇 > 泮水镇。
![](//html.hanspub.org/file/1-1750821x23_hanspub.png)
Figure 6. Ei and Ri of heavy metals content in agricultural surface soil in each countryside in Zunyi area
图6. 遵义县各乡镇表层农业土壤重金属单项生态危害系数及复合生态风险指数
4. 结论
贵州喀斯特低温成矿区表现Hg、Cd、Pb等元素异常,用相对累积频率分析方法确定该地区表层土壤重金属的基线与背景范围为:Hg 0.062~0.178,Cd 0.117~0.337,Pb 8.78~50.86,运用地球化学基线、CDF曲线可直观在图上看到每一个数据,并且反映异常值较为敏感,也反映出该地区Hg、Cd、Pb等元素的污染趋势。当地土壤自然背景值高与人为活动污染影响有很大关系,用土壤环境地球化学基线与土壤背景的历史值比较,能反映该地区农业土壤中污染趋势,为当地农业无公害环境评价提供一个明确土壤重金属的背景值范围。
用地质累积指数(Igeo)的污染程度计算遵义区域尺度土壤重金属污染程度级别,Hg、Cd、Pb等污染程度与级别分别是:Hg 24.07%为1级,60.74%为2级;Cd 13.84%为1级,80.54%为2级;Pb 46.42%为1级,43.97为2级。采用生态风险指数对遵义地区土壤重金属污染的程度与级别评价,分别是Hg 88.98%为轻微,Cd 65.02%为轻微,Pb 100%为轻微,且均小于150,属于轻微,两者相对独立,对比两种计算结果,获得一致结果,为当地无公害评价提供重要依据。
基金项目
国家自然科学基金项目(21667009);黔合支撑[2018] 2775。
NOTES
*通讯作者。