重要生态功能区土壤化学元素的空间分布模式<br/>——以大兴安岭松岭区为例

重要生态功能区土壤化学元素的空间分布模式
——以大兴安岭松岭区为例

2022-01-06于俊博周传芳梁中恺孙彦峰姜平胡宸

中国地质调查 2021年6期

于俊博，周传芳，梁中恺，孙彦峰，姜平，胡宸

(中国地质调查局哈尔滨自然资源综合调查中心，黑龙江哈尔滨 150086)

0 引言

大兴安岭是我国重要的生态功能区，也是敏感的生态脆弱区，同时大兴安岭山区是我国纬度最高的冻土区。近年来，国外学者以地球化学视角研究高纬度冻土地区的土壤问题，例如土壤中微量元素的浓度[1-2]、土壤中的重金属污染[3-5]、人类活动对土壤的影响[6]、气候变化引起冻土融化从而导致土壤中有机质发生变动[7]等。而国内对于大兴安岭冻土区土壤的研究主要侧重于土壤中的微生物和真菌群落[8-11]、土壤中的温室气体[12-13]以及土壤对植被生长的影响[14-15]等方面，而从地球化学视角分析土壤中化学元素分布的研究较少。本文依托大兴安岭山区生态地质调查项目[16]，在大兴安岭地区开展生态地球化学调查，重点研究了大兴安岭北段松岭区土壤中化学元素的空间分布模式，为区域生态保护与修复提供基础支撑。

1 研究区概况

研究区位于黑龙江省大兴安岭松岭区，处于大兴安岭北段，属于高纬度寒温带湿润区，夏短冬长，结冰期一般在7个月左右，根据加格达奇气象局日气温统计数据，计算得出2000—2019年平均气温为0.23 ℃，最低气温为-35.2 ℃。降雨期主要集中在7—9月，年平均降雨量约460 mm。

区内主要出露晚三叠世—早侏罗世混合花岗岩、中侏罗世二长花岗岩和花岗闪长岩、早白垩世正长花岗岩、二长花岗岩、花岗斑岩、花岗闪长岩等，零星分布早志留世—中泥盆世泥鳅河组大理岩、晚侏罗世玛尼吐组英安岩与英安质凝灰岩以及第四纪松散堆积物[17](图1(a))。

(a) 研究区地质简图 (b) 研究区土壤类型分布图

研究区土壤类型主要以暗棕壤和棕壤为主，伴有少量沼泽土和草甸土[18]。植被以针阔混交林为主，优势树种为兴安落叶松，次优势树种为白桦和蒙古栎，低矮草本植物主要以小叶章为主[19]。

2 样品采集与分析测试

根据大兴安岭生态地质调查的工作要求以及多目标区域地球化学调查规范[20]，结合松岭区的自然地理条件，通过试验确定对森林、沼泽生态系统地质调查样品采集的具体方法。在研究区范围内每1 km2采集1个表层土壤样品，采样深度0～20 cm，样品质量1 kg，共采集284件。每4 km2采集1个深层土壤样品，采样层位为C层(成土母质层)土壤，采样深度45～150 cm，共采集72件(图1(b))。土壤样品为多点组合采集，样品风干后按实际采样点数送样测试。在样品采集过程中，尽量避开农田、河边、道路等人为干扰较为严重的区域。

本研究分析测试了19项土壤样品的组分指标，测试工作由黑龙江省地质矿产实验测试研究中心完成。SiO2、K2O、P、Pb、Cl、Br、CaO、Zn、Na2O由 X射线荧光光谱法(XRF)测定； Ni由电感耦合等离子体发射光谱法(ICP-OES)测定； Cd由电感耦合等离子体质谱法(ICP-MS)测定； As、Hg由原子荧光法(AFS)测定； F由离子选择性电极(ISE)和玻璃电极法测定； I由光度法(COL)测定； N和有机质由容量法(VOL)测定； S和全C由高频燃烧-红外吸收光谱法测定。各元素分析数据的准确度和精密度均满足《DZ/T 0258—2014多目标区域地球化学调查规范(1∶250 000)》[20]的质量要求，分析结果可信。

3 地球化学参数特征

对土壤样品数据进行分析处理，本研究采用X±3S为临界值，剔除异点。计算剔除后的土壤中各组分的最大值、最小值、中位数、平均值、标准离差、变异系数等地球化学参数特征见表1。从变异程度看，研究区表层土壤中Cl、F、Ni、Pb、Zn、SiO2、CaO、K2O的变异系数小于25%，呈弱变异性； As、Cd、Hg、I、N、P、S、Na2O、全C的变异系数为25%～50%，呈中度变异性； Br的变异系数为51%，呈高度变异性。研究区深层土壤中Cl、F、Pb、SiO2、CaO、K2O的变异系数小于25%，呈弱变异性； As、Br、Cd、Hg、I、N、Ni、P、S、Zn、Na2O的变异系数在25%～50%，呈中度变异性；全C和有机质的变异系数在50%～75%，呈高度变异性。

表1 土壤样品地球化学参数特征

由于中值比平均值抗干扰性更好，因此采用中值作为研究区各元素在土壤中含量的估计值[22-25]。环境富集系数(q)=表层土壤元素含量/深层土壤元素含量[26]，它可以直观地反映土壤中元素的分散与富集规律。如图2所示，As、F、Pb、Zn、SiO2、CaO、Na2O、K2O的q值在0.85～1.15之间，与其自然土壤中元素含量接近，Br、Cd、Cl、Hg、I、Ni的q值在1.15～1.5之间呈弱富集，Na2O、P、S、全C的q值在1.5～4之间呈强富集。

图2 研究区元素组分环境富集系数

与中国大兴安岭弧盆系表层土壤中元素指标含量值[21]相比(表1)，研究区表层土壤中Cl、Hg、N、P、Pb、S、Zn、Na2O、全C等指标的含量更高，As、Br、Cd、F、I、Ni、SiO2、CaO、K2O等指标的含量更低。

4 元素的空间分布及分布模式

土壤是地球圈层中重要的组成部分，它受成土母质、气候、生物、地形、时间等因素的影响，同时土壤也控制着物质、能量和信息的转化[27-30]。对土壤中元素的空间分布来说，一般情况下，表层土壤更容易受生物、地表水以及人类活动等影响，使其元素分布模式发生改变。而深层土壤中元素的分布模式更多受成土母质的影响[31]。本文使用GeoIPAS软件对土壤数据进行处理分析，并利用克里金插值法对数据进行离散数据网格化，从而绘制元素含量分布图。

4.1 全C、N、P、S的分布模式

C、N、P、S是植物生长所需的必需元素，所有的氨基酸中均含有N，而P、S、C则是糖类、脂类等物质的主要化学成分。研究区内全C、N、P、S在表层土壤中呈显著富集，且高值区的面积明显比深层土壤中的面积大，但深层土壤中各元素的高值区位置在表层土壤中得到了继承。在大兴安岭山区的森林和沼泽生态系统中，地表腐殖质较厚，生物地球化学活动频繁，致使全C、N、P、S等生物必需元素在表层土壤中累积(图3)。

(a) 表层(左)与深层(右)土壤全C元素含量分布(b) 表层(左)与深层(右)土壤N元素含量分布

4.2 K、Ca、Na、Si的分布模式

从K2O、CaO、Na2O、SiO2的分布模式来看，CaO在表层土壤中含量较高，而K2O、Na2O、SiO2含量明显低于在深层土壤中的含量(图4)。CaO和Na2O在表层土壤中高值区的分布和深层土壤中高值区的分布一致，但表层土壤中高值区的强度和面积更小，低值区的强度和面积更大。

(a) 表层(左)与深层(右)土壤CaO含量分布(b) 表层(左)与深层(右)土壤K2O含量分布

研究区内基岩主要以花岗岩类为主，花岗岩风化成土的过程也是碱金属水解淋失的过程[32]。大兴安岭地区7、8月为雨季，表层土壤含水量大，而深层土壤温度低且含水量较少，因此Na2O在表层土壤中更易溶于水，并在水动力下流失。花岗岩中SiO2含量较高，从表层土壤—深层土壤—基岩，SiO2的含量随深度的增加而升高。表层土壤在地表水和风化作用的影响下，SiO2和K2O易于流失，这是表层土壤对深层土壤中SiO2和K2O的分布模式继承相对较弱的主要原因。

4.3 F、Cl、Br、I的分布模式

研究区表层土壤中Cl、Br、I的含量分别比深层土壤高43.8%、46.8%和15.9%，而F的含量比深层土壤低1.29%(表1)。F在表层土壤中继承了深层土壤高值区的位置，但表层土壤中的高值强度和范围有所减少。Cl、Br、I在表层土壤的高值区同样继承了深层土壤高值区的位置，但是表层土壤中Cl、Br、I的高值强度和范围都有明显的增加，表明Cl、Br、I在表层土壤的富集有一部分原因是受基岩影响。通过计算有机质与Br、I含量的相关性，得出有机质与Br的相关系数为0.709，有机质与I的相关系数为0.641，说明Br、I与有机质之间为正相关，Br和I除受深层土壤高值的影响，在表层土壤的富集也受到土壤上部腐殖质层的影响(图5)。

(a) 表层(左)与深层(右)土壤F元素含量分布(b) 表层(左)与深层(右)土壤Cl元素含量分布

4.4 重金属元素的分布模式

密度大于4.50 g/cm3的金属一般被认为是重金属。在天然金属中密度大于4.50 g/cm3的金属元素有54种，其中工业上认定Cu、Pb、Zn、Sn、Ni、Co、Sb、Hg、Cd、Bi等为重金属元素[31-33]。

在研究区内重点选择As、Cd、Hg、Ni、Pb、Zn进行研究。单因子指数法和内梅罗综合污染指数法是评价重金属污染的常用手段，但单纯的数理统计方法不能完全解释地质现象，因此本文除了采用单因子指数法和内梅罗综合污染指数法，还结合重金属元素分布模式图探讨研究区内重金属元素的分布模式[34-35]。

单因子污染指数计算公式为

。

(1)

式中：Pi为表层土壤中i类重金属元素的污染指数；Ci为i类重金属元素的实际测量值，10-6；Si为大兴安岭弧盆系表层土壤i类重金属元素的背景值，10-6。

内梅罗综合污染指数计算公式为

。

(2)

如表2所示，As、Cd、Ni的单因子指数小于1，表明研究区表层土壤未受As、Cd、Ni污染，Hg、Pb、Zn单因子指数大于1，表明Hg、Pb、Zn在研究区表层土壤中轻微超标。研究区表层土壤内梅罗综合指数为1.26，表明研究区内表层土壤重金属元素轻微超标[33-34]。

表2 表层土壤重金属元素单因子指数与内梅罗综合指数

如图6所示，Pb、Ni、Zn在表层土壤中的高值区均较好地继承了深层土壤中高值区的位置和规模，As在表层土壤中的大部分高值区继承了深层土壤中高值区的位置，但表层土壤个别高值区的规模和强度要比深层土壤高值区的规模和强度大。Cd在表层土壤中的高值区继承了深层土壤高值区的位置，但高值区的强度没有深层土壤强，高值区面积较深层土壤大。Hg在表层土壤的高值区强度高且范围大，结合单因子指数和重金属元素分布模式图可以看出，除了深层土壤对表层土壤的影响外，表层土壤中Hg的富集可能受地表腐殖质、森林火灾、人类活动等方面影响[36]。