1. 引言
土壤有机质是存在于土壤中所含碳的有机物,能改善土壤的物理性质及土壤养分供应能力。土壤全氮是土壤中各种形态氮素含量的总和,反映土壤供氮潜力。两者在土壤养分及植物生长发育过程中占据着重要地位,是土壤质量和土壤肥力的评价指标,因此了解区域土壤有机质和全氮含量的变化,探究其分布状况及空间变异特征对土壤碳氮库的调控和管理以及科学指导农业生产、合理利用土地等具有重要作用。实验室化验分析过程中,全氮的测定步骤繁杂,预处理麻烦,误差较大;而有机质的测定步骤较为简单,耗时短,结果也十分稳定。众多研究已表明土壤有机质与全氮关系密切,因此以简单易测的土壤有机质含量估算全氮含量具有重要意义。多年来,各界学者对不同区域、不同植被类型、不同深度等的土壤有机质和全氮进行了大量空间变异特征及相关性的研究:郑然等对农牧交错带冀北栗钙土的耕层土壤有机质和全氮的空间变异特征进行了研究 [1];张一扬等以靖西市植烟区土壤为研究对象进行有机质含量时空变异特征的研究 [2];张雅馥等以山西省朔州市平鲁区安家岭3号井工矿塌陷稳沉区为研究对象,对土壤有机质和全氮的垂直变异性和水平变异性进行研究 [3];杨海明等对丰镇市2007~2013年测定的农田土壤有机质和全氮含量进行了相关性分析并求得线性拟合模型 [4];巩万合等发现农业区不同、土壤类型不同,有机质和全氮的相关系数及回归方程各不相同 [5]。
长期以来盐碱地作为耕地补给的重要来源,在缓和人地矛盾,促进农业高效稳定发展和保障国家粮食安全等方面发挥着重要作用 [6],但是针对盐碱地土壤有机质和全氮的分布状况、空间变异特征及相关性的研究较少。因此,本文以黄河三角洲地区典型的盐碱地分布地——黄河岛为研究区,通过野外采样、实验室测定,分析该区域不同土层土壤有机质和全氮含量分布状况及空间变异特征,以期为黄河岛合理利用土地资源和科学管理农业生产提供依据,同时进行不同土层土壤有机质和全氮含量的相关性分析,确定最佳估算模型,以期通过有机质含量能够准确地估算出全氮含量,减少之后相关分析的工作量,并为其他盐碱地区域提供参考。
2. 材料与方法
2.1. 研究区概况
黄河岛(37˚57′N, 118˚4′E)位于山东省滨州市无棣县东北部,距无棣县城、滨州市区均六十公里左右,与滨州港隔河相望,东北濒临渤海,南倚黄河,四面环水,总面积25.9 km2。由于河流淤积及海潮溢流的影响,黄河岛在南北方向形成了独特的条状地形,地势南高北低,造型独特,土壤类型多为滨海盐化潮土,土壤贫瘠,含盐量高,以芦苇、柽柳、苜蓿、碱蓬等耐盐碱植物居多。属温带海洋性气候,气温起伏较小,冬无严寒,夏无酷暑,气候适宜,野生动物种类繁多。近年来,黄河岛积极进行盐碱地改良和产业开发,涉及林草间作、生态林场、现代渔业、休闲旅游等领域,各产业互惠互养,协调发展。
2.2. 土壤样品采集和测定
2018年10月外出取样,参考遥感影像,持GPS定位仪,记录采样点的经纬度和海拔高度,考虑研究区地形地貌、土地利用状况、植被分布情况等,分地类用途按照全面性、均匀性、有代表性的原则划定1 m × 1 m、5 m × 5 m、10 m × 10 m的样方后采用五点取样法取土,共取0~20 cm、20~40 cm两深度的土壤样品各51个。将样品带回实验室,经过自然风干、剔除杂质、研钵研磨过100目标准筛后装袋封存,用于土壤有机质和全氮含量的测定。其中有机质含量的测定利用改良版重铬酸钾-硫酸溶液油浴法 [7],全氮含量的测定利用全自动凯氏定氮法 [8]。
2.3. 数据处理
1) 将野外采样时GPS定位仪记录的采样点经纬度坐标、测定的有机质和全氮含量导入arcgis10.3中,然后将采样点转化为以米为单位的平面坐标系,叠加研究区边界线,形成黄河岛采样点分布图(图1),本文的研究只包含有样点分布的陆地区域,水域不涉及。
2) arcgis10.3用于土壤有机质和全氮含量的趋势分析、半方差函数分析、空间变异特征分析等;spss23用于土壤有机质和全氮含量的统计指标(最大值、最小值、标准差、平均值、变异系数等)的分析、相关性分析、回归分析等。
3. 结果与分析
3.1. 土壤有机质和全氮的描述性统计特征
通过对51个采样点不同土层土壤有机质和全氮含量进行非参数检验(K-S检验),发现其渐进显著性概率P均大于显著性水平0.05,说明不同土层土壤有机质和全氮含量均服从正态分布,具有良好的统计学意义。
![](//html.hanspub.org/file/12-2181583x10_hanspub.png)
Figure 1. Soil sampling point distribution map in the Yellow River Island
图1. 黄河岛土壤采样点分布图
![](Images/Table_Tmp.jpg)
Table 1. Descriptive statistical analysis of soil organic matter and total nitrogen content in different soil layers
表1. 不同土层土壤有机质和全氮含量的描述性统计分析
由表1,研究区0~20 cm土层土壤有机质含量变化范围为4.31~18.85 g/kg,均值为11.05 g/kg,标准差为4.15 g/kg,离散程度一般,全氮含量变化范围为0.24~0.92 g/kg,均值为0.54 g/kg,标准差为0.19 g/kg,离散程度较小;20~40 cm土层土壤有机质含量变化范围为2.51~14.25 g/kg,均值为7.23 g/kg,标准差为3.65 g/kg,离散程度一般,全氮含量变化范围为0.11~0.82 g/kg,均值为0.35 g/kg,标准差为0.17 g/kg,离散程度较小。
根据全国第二次土壤普查养分分级标准 [9],研究区不同土层土壤有机质和全氮含量总体偏低,随土层深度的加深,有机质和全氮含量不断减少。变异系数 [10] 能够反映土壤空间变异的基本情况,研究区0~20 cm土层土壤有机质和全氮含量的变异系数分别为37.56%和35.19%,20~40 cm土层土壤有机质和全氮含量的变异系数分别为50.48%和48.57%,在空间上均表现为中等变异性,且随土层深度加深,有机质和全氮含量的变异系数不断增大。
3.2. 土壤有机质和全氮的空间变异特征
3.2.1. 趋势分析
趋势分析是用三维视图进行数据拟合,以反映数据总体变化趋势及空间分布情况 [11],可为变异函数模型的获取和插值图的绘制提供基础。本文对研究区土壤样点建立X、Y、Z(经度、纬度、有机质(全氮)含量)空间坐标系,将不同土层土壤样点有机质(全氮)含量投影到东西方向和南北方向的平面上,通过投影点拟合出最优趋势线,最终得到不同土层土壤有机质和全氮含量趋势图(图2),其中X轴代表东西方向,Y轴代表南北方向,即绿线与蓝线分别表示土壤样点有机质(全氮)含量东西方向与南北方向的变化趋势。由图2,研究区0~20 cm土层土壤有机质含量和全氮含量在东西方向为凸线,呈现先增加后减少的趋势,在南北方向为直线,呈现持续减少的趋势;20~40 cm土层土壤有机质含量和全氮含量在东西方向为凸线,呈现先增加后减少的趋势,在南北方向为凹线,呈现缓慢减少的趋势。分析发现,研究区不同土层土壤有机质含量和全氮含量的变化趋势相似:东西方向呈先增加后减少趋势,南北方向呈减少趋势,但是变化幅度存在差异。
![](//html.hanspub.org/file/12-2181583x11_hanspub.png)
Figure 2. Change trend graph of soil organic matter and total nitrogen content in different soil layers
图2. 不同土层土壤有机质和全氮含量的变化趋势图
3.2.2. 半方差函数分析
半方差函数 [12] 是研究土壤养分空间变异性常用的函数,它不仅能反映土壤养分的空间变异规律,而且是克里金插值成图的基础前提。本文利用arcgis10.3的平均最邻近分析和欧式距离分析确定半方差函数的最佳步长和步数,在地统计分析模块中选择不同模型分别对土壤有机质和全氮含量进行拟合,以标准均方根误差越接近于1,标准平均误差、均方根误差和平均误差值越小作为选择最佳拟合模型的判定依据 [13],最终确定指数函数和高斯函数分别为0~20 cm土层土壤有机质和全氮含量的最佳拟合模型,球面函数和高斯函数分别为20~40 cm土层土壤有机质和全氮含量的最佳拟合模型,由半方差函数模型相关参数(表2)绘制出不同土层半方差函数图(图3)。据表2和图3分析知,研究区不同土层土壤有机质和全氮含量均存在块金值,这说明在野外采样和实验处理过程中有所误差 [14];0~20 cm和20~40 cm土层土壤有机质含量的块金系数 [15] 均为0.414,介于25%~75%之间,说明不同土层土壤有机质含量具有一定的空间相关性,变异是由结构性因素和随机性因素共同作用的结果,0~20 cm土层土壤全氮含量的块金系数为0.853,大于75%,说明其空间相关性较弱,变异受随机性因素影响较大,20~40 cm土层土壤全氮含量的块金系数为0.458,介于25%~75%之间,说明其具有一定的空间相关性,变异受结构性因素和随机性因素共同影响。
![](Images/Table_Tmp.jpg)
Table 2. Semi-variance function model parameters
表2. 半方差函数模型参数
![](//html.hanspub.org/file/12-2181583x12_hanspub.png)
Figure 3. Semi-variance function γ~h diagram of soil organic matter and total nitrogen content in different soil layers
图3. 不同土层土壤有机质和全氮含量的半方差函数γ~h图
3.2.3. 空间分布特征分析
结合不同土层土壤有机质和全氮含量趋势分析和半方差函数分析的结果选择最佳拟合模型,使用普通克立金法 [16] 未采样点的土壤有机质和全氮含量进行插值,以点推面,采用自然断点法分类,得到研究区不同土层土壤有机质和全氮含量的空间分布图(图4)。
由图4,研究区不同土层土壤有机质呈片状相间分布,含量较高的区域集中于东南部,含量较低的区域集中于北部,通过两土层对比分析,20~40 cm土层的低含量区较0~20 cm土层的低含量区大,集中分布于整个中北部,而0~20 cm土层在中北部区域既有低含量区也包含中含量区,两土层高含量分布区基本相似,原因在于该区域积极探索中度盐碱区的改良措施,引进了耐盐碱作物苜蓿品种,长时间的苜蓿种植使得该区域土体质地和构型明显改善,土壤质量也得到显著提升。研究区不同土层土壤全氮的空间分布情况相似,含量较高的区域集中于南部,含量较低的区域集中于中北部,通过两土层对比分析,0~20 cm土层高含量区分布于东南部,而20~40 cm土层主要分布于南部且面积较0~20 cm土层大,两土层中、低含量区的分布基本相似。土壤全氮的空间分布与有机质空间分布呈现出一定的相关性,即有机质含量高的区域全氮含量也相对较高,原因在于,一方面土壤中氮素的主要源头有两个,一是外在因素,即人为施入氮肥,二是内在因素,即土壤有机质进行分解提高,另一方面是因为土壤有机质中含有有机胶体,其具有强大的吸附能力,可以提高土壤的保肥能力。该分析结果与趋势分析结果基本一致,不同土层土壤有机质和全氮含量自东向西先增加后减少,自南向北不断减少,空间分布规律性较强,该区域可以根据有机质和全氮的含量状况及空间分布情况针对性的进行分区管理,提高土地利用率,充分发挥土地生产潜力。
![](//html.hanspub.org/file/12-2181583x13_hanspub.png)
Figure 4. Spatial distribution map of soil organic matter and total nitrogen content in different soil layers
图4. 不同土层土壤有机质和全氮含量空间分布图
3.3. 土壤有机质和全氮的相关性分析
3.3.1. 相关性分析
本文建立0~20 cm、20~40 cm两土层因变量(Y)土壤全氮含量关于自变量(X)土壤有机质含量的散点图(图5)。图5显示0~20 cm和20~40 cm土层的相关系数分别为0.793**、0.884**,呈显著正相关且随着土层深度加深,相关性不断提高。
备注:**表示在为0.01级别,相关性显著。
Figure 5. Scatter plot of different soil layers
图5. 不同土层的散点图
3.3.2. 最佳估算模型分析
本文利用不同土层土壤有机质含量选取线性、对数、幂、S型和指数5种模型分别对全氮含量进行估算(表3)。根据“相关性越大,精度越高,模型越优,Sig值越小,显著性越强”的基本原理,确定0~20 cm、20~40 cm土层土壤有机质含量估算全氮含量的最佳模型分别为S型模型和幂函数模型,回归方程分别为为Y = e0.099−7.338/X、Y = 0.064*X0.849。
![](Images/Table_Tmp.jpg)
Table 3. Accuracy of different models for estimating total nitrogen content of soil organic matter content
表3. 土壤有机质含量估算全氮含量的不同模型建模精度
备注:X:有机质,Y:全氮
4. 讨论与结论
1) 描述性统计分析表明,黄河岛0~20 cm土层土壤有机质含量和全氮含量的均值分别为11.05 g/kg、0.54 g/kg,变异系数分别为37.56%、35.19%,20~40 cm土层土壤有机质含量和全氮含量的均值分别为7.23 g/kg、0.35 g/kg,变异系数分别为50.48%、48.57%。两土层土壤有机质和全氮含量总体偏低,在空间上均表现为中等变异性。
2) 地统计分析表明,黄河岛不同土层土壤有机质含量和全氮含量的变化趋势相近,东西方向呈先增加后减少趋势,南北方向呈减少趋势,但是变化幅度存在差异。在0~20 cm土层,指数函数和高斯函数分别为土壤有机质和全氮含量的最佳拟合模型,该土层中土壤有机质含量具有一定的空间相关性,变异由结构性因素和随机性因素共同影响,土壤全氮含量空间相关性较弱,变异受随机性因素影响较大;在20~40 cm土层,球面函数和高斯函数分别为土壤有机质和全氮含量的最佳拟合模型,该土层土壤有机质和全氮含量均具有一定的空间相关性,变异由结构性因素和随机性因素共同影响。
3) 普通克里金空间插值结果表明,黄河岛不同土层土壤有机质和全氮的空间分布情况具有相似性,均呈片状相间分布,含量较高的区域集中于东南部和南部,含量较低的区域集中于北部。其中20~40 cm土层有机质低含量区较0~20 cm低含量区大,分布于整个中北部,而0~20 cm土层在中北部区域既有低含量区也包含中含量区;20~40 cm土层全氮高含量区较0~20 cm土层高含量区大,两土层中、低含量区的分布基本相似。该分析结果与趋势分析结果基本一致。
4) 相关性分析表明,黄河岛0~20 cm、20~40 cm土层土壤有机质含量与全氮含量的相关系数分别为0.793、0.884,呈显著正相关,且随土层深度加深,相关性增强,多种模型比较分析,确定0~20 cm、20~40 cm土层土壤有机质含量估算全氮含量的最佳模型分别为S型模型与幂函数模型,回归方程分别为Y = e0.099−7.338/X、Y = 0.064*X0.849。
基金项目
山东省自然科学基金(ZR2019MD014);山东省重点研发项目(2017CXGC0308)。
NOTES
*通讯作者。