基于空间统计学的坡耕地土壤空间soc分析.docxVIP

基于空间统计学的坡耕地土壤空间soc分析 土壤中的有机碳（si）是土壤肥料的重要物质，影响着土壤的物理、化学和生物性质。土壤侵蚀一方面造成侵蚀的景观部位SOC含量减少,同时迁移的SOC在低洼的景观部位发生累积。与土壤侵蚀过程相对应,SOC含量变化通过土壤团聚体的破坏和重新形成来影响土壤侵蚀和温室气体CO2的释放。为了弄清楚土壤再分布过程对大气CO2的净效应,需要研究不同景观单元SOC的空间分布特征。SOC作为土壤的一个特性值,是一个不均一和变化的连续体,即使在土壤质地相同的区域内,SOC在不同空间位置上也具有明显差异,即SOC存在明显的空间变异性。经典Fister统计方法在研究土壤属性的变化时通常假定研究对象是纯随机变量,可以无限次重复实验或大量观测,样本之间完全独立且服从正态分布。近年研究表明许多土壤参数在—定范围内存在空间相关性。因此仅运用经典的Fister统计方法已不能满足土壤空间异质性研究的需要。 地统计学是近几十年来发展起来的一门新的空间分析方法,适于区域化变量空间分布特征的描述。并且地统计学提供了一种最优的空间插值方法,它是基于对空间结构特征充分认识的基础上,优于其它插值方法(如几何插值、趋势面分析、普通滑动平均等),广泛应用于土壤和生态因子分布图的绘制。 近年来对SOC空间分异及影响因素的研究多集中于平坦区域不同土地利用方式和土壤管理措施方面,而对坡耕地土壤侵蚀和再分布过程引起的SOC空间分异研究较少。本文根据地统计学理论和方法(半方差函数和Kriging插值),分析东北黑土典型漫岗地形SOC在小尺度上的空间分异,用以说明土壤侵蚀和再分布过程对SOC空间分异的影响。 1 材料和方法 1.1 研究区地貌类型 采样区位于东北黑土区的中部,行政区属于吉林省德惠市松花江镇西南5 km的沫石河小流域(44°43′N,125°51′1″E)。研究区为中温带半湿润气候,年均温4.4℃,年均降水528.3 mm,降水主要集中在6—8月。研究区面积84 855.6 m2。采样区位于小流域的北坡,地形为典型的漫岗台地,包括各种典型的地貌类型:坡顶、坡肩、坡背、坡脚和坡足。利用日本Topcon公司生产的地球同步技术卫星自动聚焦水准仪GTSAF601(Geostationary Technology Satellite Auto focus 601)进行采样点位置和坡度的测量,得出坡顶坡度在0°～4.08°,平均2.6°;坡肩坡度最大,2.44°～5°,平均为3.53°;坡背在0.82°～3.47°,平均为2.23°;坡脚和坡足部位坡度最小,在0°～2.35°之间,平均为1.14°。土壤属于典型黑土亚类中的中层黑土。采样区为农田,玉米连作、顺向打垄和耕翻一直是该区的主要土壤管理方式。 1.2 土柱的土壤划分 在250 m×100 m范围内的典型黑土坡耕地上,以10 m间隔进行规则网格采样。利用直压式取土器(内径为5 cm)取表层30 cm土壤,取出的土柱分为0—10,10—20和20—30 cm层。每层采样数为168个,共504个土样。土壤样品磨碎过100目筛,SOC含量用重铬酸钾—硫酸容量法测定。 1.3 半方差图的绘制 许多学者对半方差函数及克立格插值理论都作了详细的描述。本文应用STATISTICA 6.0 统计软件对数据进行Kolomogorov—Smironov(K—S)正态检验和基本统计分析,以GS+for windows 5.0地统计学软件包进行半方差计算和半方差图的绘制。为了定量地确定研究对象的空间自相关性以及空间内插计算,GS+软件提供了5种数学模型即线性模型、带阈值的线性模型、球面模型、指数模型和高斯模型供半方差图拟合选择。 在此基础上利用Sufer7.1软件进行块段克立格(block Kriging)插值。半方差函数γ(h)一般有3个重要参数:块金值(C0)、基台值(C0+C)和变程(a)。另外分维数(D)也是评价区域化变量变异程度的一个重要参数。 2 结果与讨论 2.1 采用kolmolov—数据的正态分布检验和最大分析尺度的确定 半方差函数通过测定区域化变量分隔等距离的样点间的差异来研究变量的空间相关性和空间结构。进行空间相关分析的变量必须满足正态分布,并且由随机抽样的方式而获得。数据的非正态分布会使方差函数产生比例效应,抬高基台值和块金值,降低估计精度,使某些潜在的特征表现不明显,甚至会掩盖其固有的结构。 为了检验这些实验数据是否符合半方差函数分析的条件,我们用Kolmogorov—Smironov(K—S)正态分布检验概率(PK—S)对它们进行了检验。检验时取显著性水平α=0.05,若PK—S0.05,则认为数据服从正态分布。由表1可以得出,0—10,10—20和20—30 cm层SOC含量在0.05水平