土壤中的氮元素不仅是保障植物正常生长的必需元素, 也是衡量土壤肥力的重要指标之一。土壤全氮是指土壤中各种形态氮素含量的总和, 准确掌握土壤全氮含量的空间分布格局及其变异特征是区域合理利用土地资源、进行精准施肥的重要前提^[1-2]。与此同时, 土壤全氮的空间变异受到自然和人为等多种因素的共同影响^[3], 众多学者也纷纷研究指出, 地形^[4]、植被^[5-6]、土地利用方式^[7]以及土壤侵蚀^[8]等诸多因素均显著地影响着土壤全氮的空间分布格局。

随着精准农业的发展和科学利用土地的要求日益提升, 找到控制土壤全氮空间分布的主导因素对揭示土壤全氮分布规律起着极其重要的作用, 在全球气候变化的大背景下, 深层土壤氮库的氮储量十分巨大^[9], 并占据着极为重要的地位, 然而类似的研究多集中于对表层土壤全氮的研究中^[10-12], 土壤全氮的垂直分布特征同样需要更多的关注。Dwivedi等^[13]研究表明受到土壤深度变化的影响, 不同土层感知外界环境变化的敏感程度也存在显著地差异。Fierer等^[14]也指出深层土壤对温度和养分的变化相比表层土壤更加敏感。由此可见, 各环境因素在不同土层上对土壤全氮含量的作用程度也必然呈现加大差异, 土壤全氮的垂直分布及其影响因素仍然是当前土壤氮库的研究的焦点问题。

因此, 本研究以赤峰市敖汉旗为研究对象, 以实测数据为基础, 分析0—100 cm深度范围内5层土壤全氮含量的分布特征, 结合地理信息系统与地统计学, 对敖汉旗土壤全氮含量的空间分布进行预测, 揭示影响各层土壤全氮空间变异的主导因子。旨在为县域尺度的土壤氮库研究提供科学参考和基础数据。

1 材料与方法 1.1 研究区概况

研究区位于内蒙古自治区赤峰市东南部的敖汉旗(119°30′—120°54′E、41°42′—42°02′N), 全旗总面积约为8300 km², 地势起伏多变、总体呈南高北低特征, 海拔为300—1250 m；多年平均降水量在310—460 mm之间, 且由南向北递减；年均蒸发量为2000—2600 mm, 年平均气温为6℃, 冬季寒冷干燥, 夏季温热、降雨集中, 属于温带半干旱大陆性气候。敖汉旗处于欧亚干草原区, 地带性植被以疏林草原为主, 从南到北呈现出由森林和森林草原逐渐向干草原过渡的规律, 全旗主要分布有4个土类, 南部山地主要为棕壤和褐土, 中部黄土丘陵及黄土漫岗主要为栗钙土, 北部沙地主要以风沙土为主。敖汉旗土地利用类型多样属于农牧交错地带, 以农为主, 农牧林结合的经济类型区, 其中, 林地面积为3945 km², 全旗形成带网片、乔灌草相结合的防护林体系；耕地面积为1776 km², 约占全旗面积的21.39%；草地面积为1695 km², 占全旗面积的20.41%。

1.2 样地选择与样品采集

在野外实地调查的基础上, 综合考虑土壤、植被、地形地貌以及土地利用方式等因素, 于2014年8月在敖汉旗选取典型样点, 确保所选样地能够充分反映研究区的基本特征, 共确定182个典型样点(表 1, 图 1)。

确定样点后, 记录样地的坐标及地形地貌、土壤类型、土地利用方式等基本信息。去除样地内土壤表层的植被、枯落物等杂质, 挖掘标准土壤剖面, 按照0—20 cm, 20—40 cm, 40—60 cm, 60—80 cm, 80—100 cm划分5层, 由下至上分层取样, 每层取3个重复, 装入无菌袋带回实验室待测。土壤全氮含量采用半微量开氏蒸馏法测定；土壤机械组成采用Mastersizer3000激光粒度分析仪测定, 参照美国农业部土壤粒径分级标准, 按土壤颗粒大小确定样品粘粒(< 0.002 mm)、粉粒(0.05—0.002 mm)和砂粒(2—0.05 mm)的百分含量^[15]；土壤容重采用环刀法测定。

将敖汉旗2014年的(Landsat 8)遥感影像重新配准校正后, 在ENVI 5.1软件和ArcGIS下完成归一化植被指数(NDVI)的计算。

采用的计算公式为：

(1)

式中, LR表示红光波段的行星反射值；LNIR表示近红外波段的行星反射值。

1.3 数据分析 1.3.1 地统计学

地统计学(Geostatistics)也称为地质统计学, 于50年代初开始形成, 在法国著名统计学家G.Matheron的大量理论研究工作基础上形成一门新的统计学分支。地统计学是以变异函数理论和结构分析为基础, 在有限区域内对区域化变量进行无偏最优估计的一种方法^[16]。半变异函数又称半方差函数, 是地统计分析的特有函数。区域化变量Z(x)在点x和x+h处的值Z(x)与Z(x+h)差的方差的一半称为区域化变量Z(x)的半变异函数, 记为γ(h), 2γ(h)称为变异函数。其数学表达式为：

(2)

式中, h为两样本点空间距离；γ(h)是h的半方差函数值；N(h)是间隔距离等于h的样本点的对数；Z(x_i)为空间位置点x_i处指标的实测值；Z(x_i+h)为空间位置点x_i+h处指标的实测值。

基于地统计学原理, 采用GS+7.0软件完成半方差函数的计算和理论模型的拟合, 并结合ArcGIS10.0软件中Geostatistical Analyst模块, 对采样点进行克里格插值, 生成土壤全氮含量的空间分布图。

1.3.2 典范对应分析(Canonical correspondence analysis, CAA)

典范对应分析(CCA)是由对应分析发展而来的一种多元数据测序方法, 其通过检验两个变量的线性组合之间的相关性, 将因变量矩阵映射到描述变量矩阵中, 并在排序模型中插入回归模型, 使观测因子和环境因子的关系在排序图中直观表达, 同时可用CCA分析结果进行变量影响的贡献分离, 并准确估计自变量的贡献。其数学表达式为：

(3)

式中, x_i代表站点i的排序得分, b₀代表截距, b_j是环境变量j的回归系数, q代表环境变量的数量, 而Z_ij是站点i的环境变量j的值。

CCA排序图中箭头连线的长度表示土壤有机碳变异与环境因子的相关性, 连线越长, 相关性越大；箭头连线与排序轴的夹角表示该环境因子与排序轴相关性的大小, 夹角越小, 相关性越高；各箭头连线之间的夹角表示各环境因子间的相关性, 夹角越小, 相关性越高。研究采用CCA法分析土壤全氮含量对环境变量的响应关系, 并通过蒙特卡罗交互验证发法判断环境因素对土壤全氮影响到显著水平, CCA是使用R 3.0.1软件中的Vegan包执行^[17]。

本研究中的全部基础统计分析均在R3.0.1软件下完成。

2 结果与分析 2.1 土壤全氮含量的垂直分布特征

对研究区内182个样地数据进行统计分析, 结果如表 2所示, 研究区0—100 cm深度土壤全氮的分布范围在0.43—0.68 g/kg之间, 其平均含量随着土壤深度的增加呈现逐级递减的垂直分布规律；方差分析表明, 0—40 cm深度的土壤全氮含量显著高于40 cm深度以下的土壤全氮含量(P < 0.05)。这主要是由于表层土壤能够最直接地获取地表枯落物以及动物残体分解后的养分补给, 因此, 土壤表层全氮含量显著高于其他土层。参照我国土壤养分分级标准, 敖汉旗土壤全氮含量整体处于相对亏缺的水平(Ⅴ：0.5—0.75 g/kg), 其中60 cm深度以下的土壤全氮含量处于严重亏缺的水平(Ⅵ：< 0.5 g/kg)^[18]。就数据分布特征而言, 各层土壤全氮含量均呈现为轻度的正向右偏态分布特征, 峰度值略低。变异系数(CV)分析得知, 各层土壤全氮含量均处在0.45—0.50之间, 并为表现出明显的分布差异, 均属于中等程度的变异^[19]。

2.2 土壤全氮含量的空间变异特征

地统计学通过半方差函数模型能够定量地描述土壤全氮空间变异的结构性因素与随机性因素, 准确地分析土壤全氮的空间变异结构。模型拟合结果显示, 不同土壤深度内土壤全氮含量在各个方向上均无显著差异, 无需考虑数据的各向异特征, 对比不同模型的拟合精度, 最终得出0—100 cm深度土壤全氮的半方差函数理论模型及相关参数(表 3)。

半变异函数的拟合结果显示, 0—100 cm各层土壤全氮均适用于Spherical模型, 其拟合决定系数R²在0.71—0.82之间, 拟合残差处于较低水平, 在0.105—0.421之间, 表明拟合具有较好的拟合效果, 能够反映土壤全氮的空间变异特征。块金基台比(C₀/(C₀+ C))是反映土壤有机碳数据间空间自相关性的重要指标。研究区各层土壤全氮的C₀/(C₀+ C)处于50.5%—73.4%之间, 均属于中等强度的空间相关性, 结构性因素与随机性因素共同作用影响着土壤全氮的空间变异结构, 而随机性因素主要表现为人为活动、土地利用方式转变等, 并在研究区全氮的空间变异中占据了主要地位；结构性因素主要表现为地形、土壤母质、降雨分布格局等。随着土壤深度的增加, 随机因素对全氮空间变异的影响作用逐渐减弱, 其空间最大自相关距离(变程)也随着土壤深度的增加而逐渐减小, 这说明在深层土壤中土壤全氮受到外界干扰较弱, 相对表层土壤仍保持着较高的内在自相关关系。

就土壤全氮的空间分布格局而言(图 2), 研究区0—100 cm深度土壤全氮含量的分布趋势基本一致, 均呈现出南高北低的分布特征, 在各层土壤中全氮的低值区均集中呈片状分布在研究区的东北部, 高值区呈岛状分散于研究区的南部区域, 且40 cm以上的土壤全氮南北分布差异相对于深层土壤更加明显。

2.3 环境因素对土壤全氮含量的影响

典范对应分析(CCA)可以有效地分析影响因子对不同深度土壤全氮的影响(表 4, 图 3), CCA分析表明, 第一排序轴主要为NDVI, 海拔, 坡度对土壤全氮含量的综合反映, 对土壤全氮含量的综合解释达到91.24%, 其中, 海拔和NDVI与第一排序轴具有显著相关性(P < 0.05)。第二排序轴主要是土壤容重和土壤粘粒含量的综合反映, 土壤粘粒含量与第二排序轴均有显著的相关性(P < 0.05), 二者对土壤全氮含量的解释达到3.63%。第一排序轴和第二排序轴对土壤全氮的累积解释达到95.17%, 说明第一排序轴和第二排序轴反映了环境因子与土壤全氮的大部分信息, 两排序轴均与土壤全氮含量显著相关(P < 0.01), 可以较好地解释环境因素对土壤全氮的影响。

坡度和土壤容重与土壤全氮含量成负相关关系, NDVI, 土壤粘粒含量, 海拔与土壤全氮成正相关关系。各环境因素对土壤全氮含量的影响程度由高到低表现为海拔NDVI>粘粒含量>土壤容重>坡度。其中土壤粘粒和NDVI与表层土壤的相关性更高, 海拔对60 cm以下土壤全氮的影响更高。

3 讨论 3.1 土壤全氮的空间分布格局

土壤全氮的空间分布格局是在多种因素共同作用下而形成, 在本研究区内土壤全氮的水平分布呈现出较为明显的异质性, 各土壤深度上全氮的低值区域均位于研究区东北部, 高值区域主要位于中南部, 并呈现为由北向南逐步增加的分布趋势。这与研究区植被、土壤分布格局密切相关, 敖汉旗地处科尔沁沙地南缘, 其南北部土壤、植被分布差异较为明显, 而敖汉旗土壤由北到南为风沙土、栗钙土、褐土、棕壤, 受风沙活动的影响, 研究区北部的土壤肥力水平整体偏低；植被由北到南分布依次为沙生植被、旱生草本、森林植被^[20], 土壤和植被综合影响着土壤全氮的水平分布。小尺度上, 受植被类型的影响, 不同林份类型下的土壤全氮的分布呈现较大的差别^[21]；大尺度上, 植被群落表现出明显的地域性特征, 土壤全氮受到植被地带性分布的影响更加明显^[22]。刘庆生等^[23]研究指出内蒙古地区的地带性土壤及植被与土壤全氮的分布紧密相关, 尤其是植被覆盖度、植物种类数量、土壤有机质等分布均与土壤全氮显著相关, 这均验证了本研究的结论。

就土壤全氮的垂直分布而言, 随土壤深度的增加全氮的累积量也逐渐降低, 这也与其他学者的相关研究结果相一致^[24], 受地表植被的影响, 表层土壤最先受到枯落物分解后的补给, 因此表层土壤的氮素累积最丰富。就其空间变异特征而言, 受土壤深度的影响, 表层土壤的空间变异受随机性因素的影响更加明显, 同时深层土壤的空间自相关距离逐渐增加, 刘合满等^[25]同样研究指出表层土壤受到外界侵蚀、温度、湿度以及植被条件的差异, 导致其氮素的矿化速度在空间上的不均一性更加明显, 而深层土壤在深厚土层的封闭下, 处在相对均质稳定的环境中, 因此更易于保持空间自相关关系。这也解释了表层土壤的空间变异性更强的原因。

3.2 影响土壤全氮空间分布的主导因素

本研究指出在不同土层上影响土壤全氮空间分布的主导因素也存在较大差异, 土壤粘粒和NDVI与表层土壤的相关性更高, 海拔对60 cm以下土壤的影响更高。NDVI反映了地表植被的覆盖及生长状态, 不同的植被类型的土壤氮矿化速率差异显著, 特别是在表层土壤的表现更加明显^[26], 苏静等^[27]研究也指出土壤团聚体和植被恢复对土壤氮素的累积都发挥了积极作用, 从土壤侵蚀角度分析, 表层植被覆盖和良好的土壤团粒结构都能够有效控制土壤全氮随降雨径流的流失^[28], 而深层土壤更接近基岩和土壤母质, 同时受压实作用的影响土壤的团粒结构较差, 植被对深层土壤的补给以及微生物的活动相对表层土壤均较差, 这也就使得土壤粘粒和NDVI难以在深层土壤与土壤全氮建立更加紧密的联系。

研究区海拔在300—1250 m之间, 研究表明海拔和土壤全氮含量呈现正相关关系, 这一结论也和和其他学者的研究相符^[29]。而本研究发现深层土壤全氮含量和海拔的相关性更加明显, 这可能与海拔直接影响土壤温度有密切关系^[30], 特别在本研究区海拔相对高差近1000 m, 深层土壤温度的升高积极影响着土壤微生物数量、种类, 台喜生等^[31]研究指出土壤硝化细菌、反硝化细菌和固氮细菌在高海拔处土壤发育与低海拔处深层土壤的早期发育相类似, 氮循环细菌数量的变化受到海拔主导下植被和土壤理化因子的共同作用。这也说明了海拔在影响深层土壤氮累积中具有明显的作用。

4 结论

研究区0—100 cm土壤深度范围内的全氮的平均含量在0.43—0.68 g/kg之间, 随着土壤深度的增加逐级递减, 各层土壤全氮含量均呈现出南高北低的分布特征, 低值区集中呈片状分布在研究区的东北部, 高值区呈岛状分散于研究区的南部区域。

各层土壤全氮均属于中等强度的空间相关性, 结构性因素与随机性因素共同作用影响着土壤全氮的空间变异结构, 随着土壤深度的增加, 随机因素对全氮空间变异的影响作用逐渐减弱, 其空间最大自相关距离也随着土壤深度的增加而逐渐减小。

各环境因素对土壤全氮含量的影响程度由高到低表现为海拔> NDVI >粘粒含量>土壤容重>坡度。其中土壤粘粒和NDVI与表层土壤的相关性更高, 海拔对60 cm以下土壤的影响更高。