海伦市土壤主要微量元素空间分布特征

李秋燕，张一鹤，魏明辉，贺鹏飞

(1.中国地质调查局 沈阳地质调查中心，辽宁 沈阳 110034；2.自然资源部 黑土地演化与生态效应重点实验室，辽宁 沈阳 110034；3.辽宁省黑土地演化与生态效应重点实验室，辽宁 沈阳 110034)

0 引言

土壤中的硼、铜、锰、钼、锌(B、Cu、Mn、Mo、Zn)等元素是植物生长所必需的微量元素。微量元素作为辅酶或者活化剂参与多种酶的合成，调节植物体内的氧化还原作用，其含量影响植物体的生长[1]。土壤中缺乏微量元素会导致农作物产量和品质下降，进而通过食物链威胁人体健康[2]。近年来，国内外众多学者运用地统计学与GIS相结合的方法，对不同地区土壤中微量元素的空间变异规律进行研究，发现影响土壤中微量元素含量及分布的因素主要有自然因素(成土母质、土壤理化性质、土壤类型、气候等条件)和人为因素(耕作方式、施肥、种植制度等)[3-6]。

海伦是中国黑土硒都，是国家重要的商品粮基地。目前对于海伦土壤质量状况的研究基本都是以土壤养分以及土壤硒元素为主，有关微量元素的空间变异及空间分布规律的研究相对较少。本文以海伦地区表层土壤(0～20 cm)为研究对象，采用GIS与地统计学相结合的方法，以了解当地土壤微量元素丰缺情况，分析土壤中微量元素的空间分布情况，为当地农业产业区划，提高农作物产量和品质提供数字依据。

1 研究区概况

海伦市位于黑龙江省中部、绥化市北部，地处松嫩平原东北端，小兴安岭西麓，地形为丘陵、漫岗，平均海拔239 m。该区属中温带大陆性气候，无霜期近120 d，有效积温2 200～2 400 ℃，年降水量500～600 mm。土壤类型主要为黑土、草甸土、暗棕壤，还有少量的白浆土和水稻土，总面积为4 679 km2。

2 材料与方法

2.1 样品采集与测试

采样区域为海伦市全域。样品采集和分析测试工作均严格按照《多目标地球化学调查规范(1∶250 000)》[7]执行。按照2 km× 2 km规则格网采样法进行土壤样品采集，采样密度为1个点/km2，采样深度0～20 cm，每个土壤样品均由3个子样点组成，每4 km2组合成一个样品分析，共获取表层土壤分析样品1 175件。土壤样品的分析测试由黑龙江省地质矿产实验测试研究中心完成，分析方法配套方案及检出限见表1。以国家一级土壤标准物质(GBW)系列进行准确度、精密度监控,随机抽查异常点监控分析质量，土壤样品准确度、精密度均为100%，重复样合格率100%，数据可靠。

表1 各项指标的分析方法与检出限Table 1 The analysis method and detection limit of target elements

2.2 数据处理

文中采用SPSS25.0统计软件进行描述性统计和正态分布检验，对不服从正态分布的数据，则需要对其剔除异常值，以达到较好的地统计学插值效果。采用GS+10.0软件进行半方差函数的计算、理论模型的拟合。利用ArcGIS10.4软件中的Geostatistical Analyst模块进行趋势分析和克里金插值。

3 结果与讨论

3.1 土壤微量元素含量统计特征

海伦市土壤微量元素含量统计结果如表2所示。由表2可知，研究区土壤B、Cu、Mn、Mo、Zn平均含量分别为30.79×10-6、22.75×10-6、732×10-6、0.58×10-6和67.02×10-6。与松嫩平原背景值相比，B、Mn、Mo含量低于松嫩平原背景值，Cu、Zn与松嫩平原背景值相当。与全国土壤背景值相比，研究区土壤B和Mo严重低于全国土壤背景值，Zn较低于全国背景值，Mn含量接近全国土壤背景值的1.5倍，相对富集，Cu含量与全国土壤背景值相比差别不大。

变异系数(CV)是表征元素在空间分布上变异程度的指标，在一定程度上可以反应出元素受到人为干扰的程度[8]。一般认为，CV<10%为弱变异，10%～100%为中等强度变异，CV≥100%为强变异[9]。从表2可以看出，B、Mn、Mo、Zn变异系数大于10%，均属于中等强度变异，说明这4种微量元素在一定程度上受到了人为因素的影响，如施肥、灌溉、喷洒农药等；Cu变异系数低于10%，属于弱变异，说明海伦市Cu分布受人为影响不大。

表2 土壤微量元素含量统计分析Table 2 Description statistics of trace elements in surface soil

3.2 土壤微量元素含量趋势分析

正态分布检验是选择适宜插值模型的主要依据，正态分布的数据经地统计学插值时得到的效果更好，所以需要对偏分态的数据进行变换使其变成正态分布。采用Kolmogorov-Smirnov正态性检验方法[P(k-s)≥0．05]来检验数据是否服从正态分布，按照算术平均值加(减)3倍算术标准偏差(X±3S)的方法进行迭代剔除离群值后，研究区土壤中5种微量元素符合正态分布。

全局趋势分析可以反应土壤微量元素在研究区域上的总体变化趋势，忽略局部变异情况，从宏观上揭示研究区总体土壤微量元素变化规律[1]。利用ArcGIS中的Geostatistical Analyst模块对研究区5种微量元素进行趋势效应分析，结果如图1所示。其中X轴表示正东方向，Y轴表示正北方向，Z轴表示各样点数值的大小，横向投影面上的曲线表示EW向的全局性趋势效应变化，纵向投影面上的曲线表示SN向的全局性趋势效应变化。

图1 土壤微量元素趋势分析Fig.1 Trend analysis of soil trace elements

从图中可以看出，研究区B、Cu在EW方向和SN方向上均具有二阶趋势效应，并且都呈先增加后减少的“n”型趋势；Mn、Mo在EW和SN方向上也都具有二阶效应，但是与B、Cu不同的是，Mn、Mo由西向东呈现出一直增加的趋势，由南到北呈先减少后增加的“u”型趋势；Zn在EW和SN方向上均为“u”型二阶趋势效应。

3.3 土壤微量元素含量的空间变异结构特征

土壤微量元素的变异函数反映了不同距离的观测值之间的变化[12]。变异函数的计算需要数据符合正态分布，否则可能存在比例效应[13]。因此，对海伦地区土壤中5种微量元素经迭代剔除异常值后的数据进行半方差分析。

块金值(C0)反映的是最小抽样尺度以下变量的变异性及测量误差，属于随机变异；基台值(C0+C)表示变量的最大变异程度；块金值与基台值的比值[C0/(C0+C)]表示空间相关度，又称为块金效应，表明系统变量的空间相关性的程度，反映了空间变异影响因素中自然因素和人为因素的影响程度[14]。当[C0/(C0+C)]<25%时，说明系统具有强烈的空间相关性，受自然因素影响较多，人为因素影响较少；25%≤[C0/(C0+C)]≤75%，表明系统具有中等的空间相关性；若[C0/(C0+C)]>75%说明系统空间相关性很弱，变量受人为因素影响较大[15]。变程(a)指半变异函数的取值由初始的块金值达到基台值时，采样点的距离，表示元素的空间自相关范围。变程越大，说明土壤中该元素的均一性越强；变程越小，则意味着土壤中该元素的均一性越弱[12]。

根据表3可以看出，B和Zn符合指数模型，Cu符合球状模型，Mn、Mo符合高斯模型。从块金效应来看，5种微量元素块金效应都大于25%，最小值为52.5%，说明研究区5种微量元素空间分布受人为因素和自然因素双重影响。Mn、Mo的块金效应大于75%，空间相关性较弱，说明这两种元素受人为因素影响较多，其余3种元素块金效应在25%～75%之间，属于中等空间相关。研究区土壤微量元素含量的变程存在较大差异，B的变程较大，为633 km，表明在较大范围内存在空间自相关性，空间分布的均一性较强，Zn变程较小，表明在较小范围内存在自相关性，农作物对微量元素的吸收偏好也可能会影响土壤微量元素含量的变程大小[4]。

表3 海伦市土壤微量元素半方差理论模型与参数Table 3 Theoretical model and parameter of semi-variance function for soil trace elements in Hailun City

3.4 土壤微量元素空间分布特征

在趋势分析和半方差分析的基础上，选用普通克里金插值法绘制各元素的空间分布图(图2)。由图2可以看出，研究区内各微量元素含量都表现出较为明显的空间分布格局，与趋势分析结果显著一致。B、Cu呈现出由中部向四周逐渐递减的趋势；Mn、Mo含量呈现明显的带状特征，含量最高值出现在研究区的东北部，由东北部向西南部逐渐递减，变程较大；Zn在东西方向上高值区和低值区分布无明显特征，北部含量高于南部。

图2 表层土壤微量元素含量空间分布Fig.2 Spatial distribution of trace elements in surface soils

依据《土地质量地球化学评价规范》(DZ/T 0295—2016)对各微量元素含量进行等级划分[16]，研究区土壤养分单指标分等统计结果见表4。土壤中B明显缺乏，较缺乏(四等)和缺乏(五等)等级占到了研究区面积的98.03%；Mo相对缺乏，丰富(一等)和较丰富(二等)等级仅占到了研究区面积的22.76%；Cu、Zn以中等(三等)为主，分别占研究区面积的71.98%、56.79%；Mn主要以丰富(一等)和较丰富(二等)为主，较丰富(二等)以上占研究区土地总面积的80.87%。

表4 土壤养分单指标评价结果统计Table 4 Statistics of single index evaluation results of soil nutrients

3.5 土壤微量营养元素空间分布的影响因素

3.5.1 有机质的影响

土壤中有机质的含量影响着土壤对微量元素的吸附能力，有机质的含量主要受到地形地貌、气候、土壤质地、成土母质、土壤类型、土地利用方式等因素的影响[17]。土壤有机质通过对微量元素产生吸附、络合或沉淀作用，影响它们在土壤环境中的迁移转化[18]。土壤有机质与土壤微量元素含量的相关性分析见表5。土壤有机质与B、Mn、Mo、Zn呈显著相关关系，其中与B为显著负相关关系，与Mn、Mo、Zn为显著正相关关系，与Mo、Mn的相关性最大，与Cu相关性不明显。有机质对于土壤中B的吸附和解吸机理至今仍未搞清，有研究表明土壤有机质与B呈正相关关系[18]，也有学者认为土壤中有机质的大量存在是导致B缺乏的一个原因[19]。本研究中土壤有机质与B呈显著负相关关系，有可能是因为吸附在有机质上的B的解吸有滞后现象，造成了它们之间的负相关关系[20]。

表5 土壤微量元素和有机质相关性分析Table 5 Correlation analysis of soil trace elements and organic matter

微量元素间的相关分析表明，除了B和Zn不存在相关关系外，其余元素均存在明显的正相关或负相关关系，并且Mn、Mo、Zn之间存在着显著正相关关系，表明这3种元素存在着共生关系。Cu与Zn存在正相关的原因可能是这两个元素均为亲硫元素，具有相似的地球化学性质[21]。相关性分析可以说明，由于不同元素的地球化学性质不同，微量元素之间具有共生关系，也有拮抗关系。

3.5.2 成土母质的影响

成土母质是影响土壤微量元素含量的重要因素，成分不同的成土母质下，土壤中元素含量差异较明显[22]。由表6可以看出，B的平均含量表现为冲湖积>坡洪积>冲积>冰积;Cu的平均含量表现为坡洪积>冲湖积>冲积>冰积;Mn和Mo的平均含量表现为冰积>冲积>冲湖积>坡洪积;Zn的平均含量表现为冰积>冲湖积>坡洪积>冲积。不同母质发育的土壤，微量元素含量分布不同。Mn、Mo、Zn含量在研究区东北部明显较高，可能原因是研究区北部成土母质为冰积，在冰川沉积作用下，岩石碎屑沉积这些元素较多。

3.5.3 土壤类型的影响

研究区土壤类型为暗棕壤、草甸土和黑土。由表6可以看出，Mn、Mo、Zn的平均含量表现为暗棕壤含量最高，草甸土和黑土中含量相近；Cu的平均含量在3种土壤类型中分布较均匀；B的平均含量表现为黑土>草甸土>暗棕壤。土壤类型对有机质含量影响较大，研究区暗棕壤中土壤有机质的含量为45.1×10-6，高于草甸土和黑土中有机质含量，由于土壤中的有机质易与金属元素形成螯合物，从而截留于土壤表层中[23]，因此使得Mn、Mo、Zn在暗棕壤中含量最高。由于B和有机质为显著负相关关系，所以B在暗棕壤中含量最低。

3.5.4 土地利用类型的影响

土地利用是自然因素和人类活动影响的综合过程，不同的轮作制度、管理模式、施肥方式等都可能会导致土壤微量元素分布的差异[24]。由表6可以看出，B平均含量表现为草地>旱地>沼泽>林地>水田，Cu为水田>旱地>草地>沼泽>林地，Mn为林地>沼泽>草地>水田>旱地，Mo为林地>沼泽>草地>旱地>水田，Zn为林地>沼泽>水田>草地>旱地。Mn和Mo含量在不同土地利用类型下具有较为明显的差异，说明人类活动对这两个微量元素含量影响较大。

表6 不同成土母质、土壤类型和土地利用类型的土壤微量元素平均含量Table 6 Average contents of trace elements in soils of different soil parent materials，soil types and land use patterns

4 结论

1)研究区土壤中微量元素B和Mo含量分别为30.79×10-6、0.58×10-6，远小于松辽平原背景值和全国土壤背景值。B、Mn、Mo、Zn属于中等强度变异，Cu属于弱变异。

2)通过半变异函数模型分析可知，Mn和Mo的块金效应大于75%，属于中等空间自相关，空间变异受到自然因素和人为因素的双重影响，其余3种元素空间自相关较强，空间变异主要受到自然因素的影响。

3)研究区土壤5种微量元素空间分布各异，B和Cu含量由中部向四周逐渐递减，Mo和Mn含量呈现明显的带状特征，由东北部向西南部逐渐递减，Zn含量北部高于南部。5种微量元素之间既具有共生作用也具有拮抗作用。同时，土壤有机质含量、成土母质、土壤类型和土地利用类型对土壤微量元素的含量均有不同程度的影响。