基于地统计学的耕层土壤有效态微量元素空间变异特征研究.docxVIP

基于地统计学的耕层土壤有效态微量元素空间变异特征研究

基于地统计学的耕层土壤有效态微量元素空间变异特征研究

引言

土壤有效态微量元素是维持土壤肥力、保障作物生长的重要基础，其含量高低直接影响作物产量与品质。锌（Zn）、铁（Fe）、锰（Mn）、铜（Cu）、硼（B）等微量元素作为作物生长发育的必需营养元素，尽管需求量较低，但在植物体内参与酶促反应、光合作用、蛋白质合成等关键生理过程，缺乏或过量均会导致作物生长障碍甚至减产[1]。耕层土壤（0-20cm）是作物根系主要活动层，其有效态微量元素的空间分布特征对精准施肥、土壤改良及农业可持续发展具有重要意义[2]。

地统计学作为研究空间变异性的核心方法，通过半方差函数、克里金插值等工具，可量化变量的空间自相关性与分布规律，为区域土壤养分管理提供科学依据[3]。本研究以XX省XX市XX县为研究区，基于网格化采样数据，运用地统计学方法分析耕层土壤有效态Zn、Fe、Mn、Cu、B的空间变异特征，揭示其影响因素，为该区域微量元素精准管理提供理论支撑。

研究区域概况

研究区位于XX省XX市XX县（东经116°32′-117°05′，北纬34°18′-34°45′），总面积1286km2，属暖温带半湿润季风气候，年均气温14.6℃，年均降水量656mm，无霜期210天。地形以平原（占62%）、岗地（占25%）和丘陵（占13%）为主，土壤类型主要为褐土（占48%）、潮土（占37%）和砂姜黑土（占15%），耕地面积约7.2万hm2，主要种植小麦、玉米、大豆等作物，是XX省重要的粮食生产基地[4]。研究区农业历史悠久，长期施用化肥、有机肥，土壤微量元素含量受成土母质、耕作方式、施肥习惯及地形等因素综合影响，空间变异特征复杂。

研究方法

1.样品采集与测定

2023年10月（作物收获后），研究区按照“网格法”布设采样点，网格间距500m，共布设312个采样点（图1，注：此处无图）。每个采样点采用“S”形布点法采集0-20cm耕层土壤，去除石砾、植物根系等杂物，四分法保留1kg样品，标记后带回实验室。

样品经风干、研磨、过2mm尼龙筛后，测定有效态微量元素含量：有效态Zn、Fe、Mn、Cu采用DTPA浸提-原子吸收光谱法（AAS）测定[5]；有效态B采用沸水浸提-姜黄素比色法测定[6]。每批样品设置2个空白样和10%的平行样，确保测定精度（相对标准偏差RSD5%）。

2.数据处理与地统计分析

（1）描述性统计

采用SPSS26.0软件计算各微量元素的最小值（Min）、最大值（Max）、平均值（Mean）、标准差（SD）、变异系数（CV）、偏度（Skewness）、峰度（Kurtosis），分析数据分布特征。变异系数评价变异程度：CV≤10%为弱变异，10%CV≤100%为中等变异，CV100%为强变异[7]。

（2）地统计分析

使用GS+9.0软件进行半方差函数拟合，计算块金值（C0）、基台值（C0+C）、块金基台比（C0/(C0+C)）、变程（A）等参数，量化空间自相关性。块金基台比25%为强空间相关，25%-75%为中等空间相关，75%为弱空间相关[8]。理论模型选择依据决定系数（R2）和残差平方和（RSS），R2越大、RSS越小，拟合效果越好。

基于最优半方差函数模型，采用ArcGIS10.8软件进行普通克里金插值，绘制各微量元素空间分布图，识别高值区、低值区分布规律。

结果与分析

1.描述性统计特征

研究区耕层土壤有效态微量元素描述性统计结果见表1。各微量元素含量存在较大差异：有效态B含量最高（平均0.48mg/kg），其次是Fe（平均9.32mg/kg）、Mn（平均5.76mg/kg）、Zn（平均1.35mg/kg）、Cu（平均0.92mg/kg）。变异系数以B最大（CV=68.3%），Mn次之（CV=52.1%），Fe最小（CV=31.2%），均属于中等变异（10%CV≤100%），表明受随机因素（如施肥、耕作）和结构性因素（如母质、地形）共同影响[9]。

偏度系数中，Zn（1.23）、Cu（1.15）为右偏分布，B（0.87）接近正态分布；Fe（-0.62）、Mn（-0.51）为左偏分布，峰度系数均大于3，表明数据存在尖峰厚尾特征，需进行对数转换以满足地统计学正态假设[10]。转换后各微量元素偏度绝对值均1，峰度接近3，符合正态分布要求（P0.05）。

表1耕层土壤有效态微量元素描述性统计结果

|元素|Min(mg/kg)|Max(mg/kg)|Mean(mg/kg