以黄河三角洲地区典型的盐碱地分布地——黄河岛为研究区,结合传统统计学、地统计学及GIS技术,研究盐碱地不同土层土壤有机质和全氮的空间变异特征,并探讨两者的相关性。结果表明:1) 盐碱地不同土层土壤有机质和全氮含量总体偏低,具有中等空间变异性;2) 不同土层土壤有机质含量及0~20 cm土层土壤全氮含量具有一定的空间相关性,变异受结构性因素和随机性因素影响,20~40 cm土层土壤全氮含量空间相关性较弱,变异受随机性因素影响大;3) 不同土层土壤有机质和全氮含量分布状况相似,片状相间分布,东西方向呈先增加后减少趋势,南北方向呈减少趋势;4)不同土层土壤有机质和全氮含量具有强正相关性,0~20 cm和20~40 cm土层土壤有机质含量估算全氮含量的最佳模型分别为S型模型和幂函数模型,回归方程分别为Y = e0.099−7.338/X、Y = 0.064*X0.849。 Taking typical distribution of saline-alkali land in the Yellow River Delta-the Yellow River Island as research area, applying traditional statistics, Geostatistics and GIS technology, the research ana-lyzed the spatial variability characteristics of soil organic matter and total nitrogen in different soil layers in saline-alkali land and explored the correlation between them. The results showed that: 1) The contents of soil organic matter and total nitrogen in different soil layers were generally low and with moderate spatial variability; 2) Soil organic matter in different soil layers and soil total nitro-gen in the soil layers of the depth in 0~20 cm had a certain spatial correlation, and the variation was affected by structure factors and random factors, and the spatial correlation of soil total nitrogen in the 20~40 cm soil layer was relatively weak, and the variation is greatly affected by random factors; 3) The distributions of soil organic matter and total nitrogen in different soil layers are similar. They were distributed in patches, the overall trend was increased and then decreased in the direction of east-west, and decreased in the direction of north-south; 4) There was a strong positive correlation between soil organic matter and total nitrogen content in different soil layers. The best models used for estimating total nitrogen content by soil organic matter content were the S-type model and the power function model respectively in the depth of 0~20 cm and 20~40 cm soil layers, the regression equations were Y = e0.099−7.338/X、Y = 0.064*X0.849.
以黄河三角洲地区典型的盐碱地分布地——黄河岛为研究区,结合传统统计学、地统计学及GIS技术,研究盐碱地不同土层土壤有机质和全氮的空间变异特征,并探讨两者的相关性。结果表明:1) 盐碱地不同土层土壤有机质和全氮含量总体偏低,具有中等空间变异性;2) 不同土层土壤有机质含量及0~20 cm土层土壤全氮含量具有一定的空间相关性,变异受结构性因素和随机性因素影响,20~40 cm土层土壤全氮含量空间相关性较弱,变异受随机性因素影响大;3) 不同土层土壤有机质和全氮含量分布状况相似,片状相间分布,东西方向呈先增加后减少趋势,南北方向呈减少趋势;4)不同土层土壤有机质和全氮含量具有强正相关性,0~20 cm和20~40 cm土层土壤有机质含量估算全氮含量的最佳模型分别为S型模型和幂函数模型,回归方程分别为Y = e0.099−7.338/X、Y = 0.064*X0.849。
盐碱地,有机质,全氮,空间特征,相关性
Huan Yu, Xue Chen, Tianying Zhong, Jian Liu, Ailing Wang*
College of Resources and Environment, Shandong Agricultural University, Tai’an Shandong
Received: Mar. 22nd, 2021; accepted: Apr. 16th, 2021; published: Apr. 23rd, 2021
Taking typical distribution of saline-alkali land in the Yellow River Delta-the Yellow River Island as research area, applying traditional statistics, Geostatistics and GIS technology, the research analyzed the spatial variability characteristics of soil organic matter and total nitrogen in different soil layers in saline-alkali land and explored the correlation between them. The results showed that: 1) The contents of soil organic matter and total nitrogen in different soil layers were generally low and with moderate spatial variability; 2) Soil organic matter in different soil layers and soil total nitrogen in the soil layers of the depth in 0~20 cm had a certain spatial correlation, and the variation was affected by structure factors and random factors, and the spatial correlation of soil total nitrogen in the 20~40 cm soil layer was relatively weak, and the variation is greatly affected by random factors; 3) The distributions of soil organic matter and total nitrogen in different soil layers are similar. They were distributed in patches, the overall trend was increased and then decreased in the direction of east-west, and decreased in the direction of north-south; 4) There was a strong positive correlation between soil organic matter and total nitrogen content in different soil layers. The best models used for estimating total nitrogen content by soil organic matter content were the S-type model and the power function model respectively in the depth of 0~20 cm and 20~40 cm soil layers, the regression equations were Y = e0.099−7.338/X, Y = 0.064*X0.849.
Keywords:Saline-Alkali Land, Organic Matter, Total Nitrogen, Spatial Characteristics, Correlation
Copyright © 2021 by author(s) and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/
土壤有机质是存在于土壤中所含碳的有机物,能改善土壤的物理性质及土壤养分供应能力。土壤全氮是土壤中各种形态氮素含量的总和,反映土壤供氮潜力。两者在土壤养分及植物生长发育过程中占据着重要地位,是土壤质量和土壤肥力的评价指标,因此了解区域土壤有机质和全氮含量的变化,探究其分布状况及空间变异特征对土壤碳氮库的调控和管理以及科学指导农业生产、合理利用土地等具有重要作用。实验室化验分析过程中,全氮的测定步骤繁杂,预处理麻烦,误差较大;而有机质的测定步骤较为简单,耗时短,结果也十分稳定。众多研究已表明土壤有机质与全氮关系密切,因此以简单易测的土壤有机质含量估算全氮含量具有重要意义。多年来,各界学者对不同区域、不同植被类型、不同深度等的土壤有机质和全氮进行了大量空间变异特征及相关性的研究:郑然等对农牧交错带冀北栗钙土的耕层土壤有机质和全氮的空间变异特征进行了研究 [
长期以来盐碱地作为耕地补给的重要来源,在缓和人地矛盾,促进农业高效稳定发展和保障国家粮食安全等方面发挥着重要作用 [
黄河岛(37˚57′N, 118˚4′E)位于山东省滨州市无棣县东北部,距无棣县城、滨州市区均六十公里左右,与滨州港隔河相望,东北濒临渤海,南倚黄河,四面环水,总面积25.9 km2。由于河流淤积及海潮溢流的影响,黄河岛在南北方向形成了独特的条状地形,地势南高北低,造型独特,土壤类型多为滨海盐化潮土,土壤贫瘠,含盐量高,以芦苇、柽柳、苜蓿、碱蓬等耐盐碱植物居多。属温带海洋性气候,气温起伏较小,冬无严寒,夏无酷暑,气候适宜,野生动物种类繁多。近年来,黄河岛积极进行盐碱地改良和产业开发,涉及林草间作、生态林场、现代渔业、休闲旅游等领域,各产业互惠互养,协调发展。
2018年10月外出取样,参考遥感影像,持GPS定位仪,记录采样点的经纬度和海拔高度,考虑研究区地形地貌、土地利用状况、植被分布情况等,分地类用途按照全面性、均匀性、有代表性的原则划定1 m × 1 m、5 m × 5 m、10 m × 10 m的样方后采用五点取样法取土,共取0~20 cm、20~40 cm两深度的土壤样品各51个。将样品带回实验室,经过自然风干、剔除杂质、研钵研磨过100目标准筛后装袋封存,用于土壤有机质和全氮含量的测定。其中有机质含量的测定利用改良版重铬酸钾-硫酸溶液油浴法 [
1) 将野外采样时GPS定位仪记录的采样点经纬度坐标、测定的有机质和全氮含量导入arcgis10.3中,然后将采样点转化为以米为单位的平面坐标系,叠加研究区边界线,形成黄河岛采样点分布图(图1),本文的研究只包含有样点分布的陆地区域,水域不涉及。
2) arcgis10.3用于土壤有机质和全氮含量的趋势分析、半方差函数分析、空间变异特征分析等;spss23用于土壤有机质和全氮含量的统计指标(最大值、最小值、标准差、平均值、变异系数等)的分析、相关性分析、回归分析等。
通过对51个采样点不同土层土壤有机质和全氮含量进行非参数检验(K-S检验),发现其渐进显著性概率P均大于显著性水平0.05,说明不同土层土壤有机质和全氮含量均服从正态分布,具有良好的统计学意义。
图1. 黄河岛土壤采样点分布图
土壤指标 | 土层深度(cm) | 最大值(g/kg) | 最小值(g/kg) | 平均值(g/kg) | 标准差(g/kg) | 变异系数(%) | Pk-s |
---|---|---|---|---|---|---|---|
有机质 | 0~20 | 18.85 | 4.31 | 11.05 | 4.15 | 37.56 | 0.152 |
20~40 | 14.25 | 2.51 | 7.23 | 3.65 | 50.48 | 0.204 | |
全氮 | 0~20 | 0.92 | 0.24 | 0.54 | 0.19 | 35.19 | 0.103 |
20~40 | 0.82 | 0.11 | 0.35 | 0.17 | 48.57 | 0.121 |
表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,离散程度较小。
根据全国第二次土壤普查养分分级标准 [
趋势分析是用三维视图进行数据拟合,以反映数据总体变化趋势及空间分布情况 [
图2. 不同土层土壤有机质和全氮含量的变化趋势图
半方差函数 [
土壤指标 | 土层深度/cm | 函数模型 | 块金值 | 基台值 | 块金系数/% | 变程/m |
---|---|---|---|---|---|---|
有机质 | 0~20 | 指数模型 | 5.661 | 13.676 | 0.414 | 1710.22 |
20~40 | 球面函数 | 4.555 | 10.997 | 0.414 | 2018.21 | |
全氮 | 0~20 | 高斯模型 | 0.029 | 0.034 | 0.853 | 4370.86 |
20~40 | 高斯模型 | 0.011 | 0.024 | 0.458 | 814.983 |
表2. 半方差函数模型参数
图3. 不同土层土壤有机质和全氮含量的半方差函数γ~h图
结合不同土层土壤有机质和全氮含量趋势分析和半方差函数分析的结果选择最佳拟合模型,使用普通克立金法 [
由图4,研究区不同土层土壤有机质呈片状相间分布,含量较高的区域集中于东南部,含量较低的区域集中于北部,通过两土层对比分析,20~40 cm土层的低含量区较0~20 cm土层的低含量区大,集中分布于整个中北部,而0~20 cm土层在中北部区域既有低含量区也包含中含量区,两土层高含量分布区基本相似,原因在于该区域积极探索中度盐碱区的改良措施,引进了耐盐碱作物苜蓿品种,长时间的苜蓿种植使得该区域土体质地和构型明显改善,土壤质量也得到显著提升。研究区不同土层土壤全氮的空间分布情况相似,含量较高的区域集中于南部,含量较低的区域集中于中北部,通过两土层对比分析,0~20 cm土层高含量区分布于东南部,而20~40 cm土层主要分布于南部且面积较0~20 cm土层大,两土层中、低含量区的分布基本相似。土壤全氮的空间分布与有机质空间分布呈现出一定的相关性,即有机质含量高的区域全氮含量也相对较高,原因在于,一方面土壤中氮素的主要源头有两个,一是外在因素,即人为施入氮肥,二是内在因素,即土壤有机质进行分解提高,另一方面是因为土壤有机质中含有有机胶体,其具有强大的吸附能力,可以提高土壤的保肥能力。该分析结果与趋势分析结果基本一致,不同土层土壤有机质和全氮含量自东向西先增加后减少,自南向北不断减少,空间分布规律性较强,该区域可以根据有机质和全氮的含量状况及空间分布情况针对性的进行分区管理,提高土地利用率,充分发挥土地生产潜力。
图4. 不同土层土壤有机质和全氮含量空间分布图
本文建立0~20 cm、20~40 cm两土层因变量(Y)土壤全氮含量关于自变量(X)土壤有机质含量的散点图(图5)。图5显示0~20 cm和20~40 cm土层的相关系数分别为0.793**、0.884**,呈显著正相关且随着土层深度加深,相关性不断提高。
图5. 不同土层的散点图
本文利用不同土层土壤有机质含量选取线性、对数、幂、S型和指数5种模型分别对全氮含量进行估算(表3)。根据“相关性越大,精度越高,模型越优,Sig值越小,显著性越强”的基本原理,确定0~20 cm、20~40 cm土层土壤有机质含量估算全氮含量的最佳模型分别为S型模型和幂函数模型,回归方程分别为为Y = e0.099−7.338/X、Y = 0.064*X0.849。
模型 | 0~20 cm | 20~40 cm | ||
---|---|---|---|---|
回归方程 | R2 | 回归方程 | R2 | |
线性 | Y = 0.037X + 0.137 | 0.629 | Y = 0.041X + 0.042 | 0.782 |
对数 | Y = −0.339 + 0.379*ln(X) | 0.672 | Y = −0.167 + 0.279*ln(X) | 0.751 |
幂 | Y = 0.078*X0.804 | 0.717 | Y = 0.064*X0.849 | 0.792 |
S | Y = e0.099−7.338/X | 0.743 | Y = e−0.325−4.715/X | 0.755 |
指数 | Y = 0.216*e0.077*X | 0.657 | Y = 0.124*e0.125*X | 0.781 |
表3. 土壤有机质含量估算全氮含量的不同模型建模精度
备注:X:有机质,Y:全氮
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)。
于 欢,陈 雪,钟天颖,刘 健,王瑷玲. 盐碱地土壤有机质和全氮的空间变异特征及相关性分析——以黄河岛为例Analysis of Spatial Variability Characteristics and Correlation of Soil Organic Matter and Total Nitrogen in Saline-Alkali Land—The Case Study of the Yellow River Island[J]. 农业科学, 2021, 11(04): 332-340. https://doi.org/10.12677/HJAS.2021.114048