苏打盐碱草地不同研究尺度土壤pH值空间变异特征_王欣慰.pdf

松嫩平原重度退化盐碱草地为典型的苏打盐碱土，对其不同尺度下的土壤异质性特征的了解是理解其退化过程、对其进行合理有效改造的基础。值的变化是土地盐碱化的一个关键性指标，本文对杜尔伯特蒙古族自治县他拉哈镇典型退化盐碱草地 值的空间异质性在 个尺度上进行了研究，这 个尺度分别为大尺度（.，取样间隔 ）、中尺度（.，取样间隔 ）、小尺度（.，取样间隔 ）和微尺度（，取样间隔.）。结果表明：只有大尺度取样时符合正态分布，其他尺度均不符合正态分布。从大尺度到微尺度的变异系数分别为、.、.和.。中尺度和微尺度下随机变异所占比例均比较小，都表现为强空间自相关的特点，大尺度和小尺度为中等空间自相关。在大尺度取样时对 值的预测不够准确。中尺度取样时预测值有相当的参考价值。小尺度和微尺度取样时预测值精度高，四个研究尺度下松嫩草地土壤 值的变程为 、左右、左右和小于 。综合四个尺度的空间特征，研究地的盐碱化草地形成约为.的大尺度斑块，每个斑块内约有 千个约 的 值变化单位，土壤 值和植物之间在小于 的尺度上相互作用，采取保护和改良措施也应以此尺度下土壤的变化为基础。关键词：尺度；值；空间变异中图分类号：.文献标志码：文章编号：（）基金项目：年度河北省社科基金项目（）。第一作者：王欣慰（），男，硕士，主要从事森林植物资源学方面的研究。通讯作者：王洪峰（），男，高级实验师，主要从事森林植物资源学方面的研究。，（，；，；，）：.（，），.，.，.，.，.（），（），.，.：；松嫩平原生态脆弱，土地盐碱化日趋严重，相比于 世纪 年代，目前该地区盐碱化面积已增加近 万公顷。相比于其他区域，松嫩平原土壤空间质异性更为特殊，通常情况下，其 值变化过程更为缓慢，但因松嫩草地土壤多为苏打盐碱土，值的变化在较小尺度上又会有较为剧烈的变化。对于松嫩草原来说，在土壤性质空间变异特 广 州 化 工 年 月性不明的情况下，对盐碱草地改良前后的评价、精确化改良措施的合理尺度等方面的研究和应用都缺乏坚实的基础。所以对土壤性质的空间变异特征的研究是对此区域生态环境保护和改良的理论基础。土壤的化学性质随着时间和空间的变化并非遵循简单的随机函数关系，而是一种自相关函数关系。自 世纪 年代以来，部分外国学者借用矿产分析的经验，发展了主要以半方差函数和克里格插值等方法为基础的对土壤空间异质性进行分析的方法即地统计学，在地统计学中，对有限区域中区域化变量取值进行无偏最优估计的方法为空间局部估计。同时，这种新方法也被利用在土壤变异的研究中，并受到人们的重视。因为空间变异是尺度的函数，所以在这类研究中，研究尺度的选择是一个十分关键的因素，某一尺度下的研究结论会掩盖其他尺度的规律，这在前人的研究中也得到证实。值的变化是一个苏打盐碱土盐碱化的一个关键性指标，可以很好的反映土壤的理化性质，且能够影响土壤的其他性质。本文在不同尺度下对其统计特性的空间变异性进行了研究，力图找到 值在松嫩草地苏打盐碱土环境下的一些变化规律。材料和方法.研究区概况研究区域位于松嫩平原西部杜尔伯特蒙古族自治县他拉哈镇，地处，平均海拔 ，属温带大陆性季风气候，年均降雨量为 左右，集中在 月，年平均气温.，日平均气温稳定在 以上的日数 ，无霜期 ，封冻期在 月上旬，解冻期在 月上旬。.样品采集选择大尺度（.，取样间隔 ）、中尺度（.，取样间隔 ）、小尺度（.，取样间隔）和微尺度（，取样间隔.）个尺度进行取样研究。所有 个尺度样地均设在杜尔伯特蒙古族自治县他拉哈镇境内盐碱化草地。选择一片盐碱草地，设置 样地为大尺度样地，在此样地上划分 网格，每个网格点设 小样方，取小样方四个顶点和中心点 个点 表层土混合。选 的盐碱化草地为中尺度研究取样区域。在此区内以划分 网格，每个网格点处设 小样方，取小样方四个顶点和中心点 个点 表层土混合。小尺度样地为 长方形样地，在此区内以划分 网格，每个网格点处设 小样方，取小样方四个顶点和中心点 个点 表层土混合。微尺度样地为 长方形样地，在此区内以划分 网格，网格点上选取 个点，每个点取 表层土。.值的测定以电位法测定土壤 值。.数据处理使用 软件进行传统统计分析，软件计算半方差函数，.软件结合 软件进行空间任意点值预测。结果与分析.不同研究尺度 值传统统计分析对不同研究尺度土壤 值进行了传统统计分析（表）。只有大尺度取样时符合正态分布，其他尺度均不符合正态分布。说明在大尺度下，值的变化相对缓和，土壤性质大体均一。而相对较小的尺度下，土壤受放牧、土地开垦、小型湖泊其它人为干扰较剧烈，使其性质变化较复杂。除大尺度外，最大值相差不大，最小值和均值有比较明显的差别，值的上限变化不大，下限变化较大。研究尺度对变异系数影响较大，从大尺度到微尺度的变异系数分别为、.、.和.。大尺度上，土壤主要受气候、水文、土壤母质和大范围土地利用方式的影响，这些因素对土壤 值的影响在不同地段有较大的差异。小尺度和微尺度的变异系数明显大于中尺度，土壤在较小尺度上 值受放牧、微地形、植被等随机影响较大，能显著影响土壤表层的因素多且复杂，使较小尺度上土壤表层 值变化剧烈。表 不同研究尺度 值传统统计分析 尺度样点数分布类型偏度峰度均值最小值最大值变异系数.大尺度.中尺度.小尺度.微尺度.不同研究尺度土壤 值空间变异首先考察不同研究尺度下取样点数值的统计分布特征，进行正态检验，不符合正态分布的进行正态转换。应用 地统计学软件，计算半方差函数，公式为：（）（）（）（）（）式中：（）为半方差函数值；（）为间距为向量 的点对总数；（）为区域化变量 在 处的实测值；（）为与 距离为向量 的样点的值，本研究中半方差函数的有效滞后距都设为最大采样间隔的。软件计算半方差函数，为块金值（）；为结构方差；为基台值（）；为滞后距离；为变程（）。块金值表示随机变异的大小，基台值表示系统内的总变异，但由于自身因素和测量单位的影响较大，故块金值和基台值不能用于比较不同变量间的随机变异，而两者比值则能反映块金方差占总空间异质性变异的大小，表示土壤属性的空间依赖性。当该比值时，空间变量为强烈的空间自相关，在 时，为中等空间自相关，时为弱空间自相关。变程与观测尺度及取样尺度上影响土壤水分的各种生态过程及第 卷第 期王欣慰，等：苏打盐碱草地不同研究尺度土壤 值空间变异特征 其相互作用有关，能表示属性因子空间自相关范围的大小。在变程之内，变量具有空间自相关性，反之则不存在。表 不同研究尺度半方差函数 （）大尺度.中尺度.小尺度.微尺度.图 大研究尺度半方差函数图.图 中研究尺度半方差函数图.图 小研究尺度半方差函数图.图 微研究尺度半方差函数图.不同研究尺度下对半变异函数进行拟合，对比选择最优的拟合模型。大研究尺度、中研究尺度、小研究尺度、微研究尺度的最优拟合模型分别为指数模型、高斯模型和球状模型。模型拟合决定系数较高，分别为.，说明理论模型能较好的反映 值的空间变异特征。不同研究尺度的块金值 基台值为.，中尺度和微尺度为强空间自相关，大尺度和小尺度为中等空间自相关。变程不是一成不变的，不同研究尺度下变程也不同。大研究尺度的变程为 ，中小研究尺度的变程差别不大，分别为.和.，微研究尺度变程为.。变程是空间自相关的范围，一般来说，变程范围内存在空间自相关，变程范围外空间自相关消失。从不同研究尺度的结果看，值的变化存在不同的结构层次。松嫩草地土壤 值的变化基本可分为三个层次：小于 范围内的变化、左右范围内的变化，范围内的变化。在这三个范围间不存在空间自相关（表，图 图）。.不同研究尺度土壤 值预测首先绘制不同研究尺度采样图（图 图），在此基础上以克里格法进行插值计算，并绘制 值等值线图（图 图）。在 个尺度试验区内各随机选 个点，或米尺取得坐标值后，结合 预测数值，然后在这些随机选取的坐标点处采样、测量，对预测值进行检验。图 大研究尺度采样图.图 中研究尺度采样图.广 州 化 工 年 月图 小研究尺度采样图.图 微研究尺度采样图.图 大尺度取样等值线图.图 中尺度取样等值线图.图 小尺度取样预测结果.图 微尺度取样预测结果.在回归分析时可以计算 和残差平方和来标志回归方程和实际值的符合程度，在进行插值计算和在 中进行预测时，无法得到回归分析中那样简单明了的公式，但仍然可以用同样的原理计算 和残差来检验预测的准确程度（表）。在大尺度取样时，预测 值时 小残差平方和很大，说明这种方式的大尺度取样和计算虽然在插值时可以得到比较好的结果（插值计算的.），但是其实对 值的预测不够准确。中尺度的预测图过渡平缓，显示出区域内明显的高 值聚集区成条带分布，和大尺度相比，而 有所提高但仍然比较小，中尺度取样对具体点做预测准确率仍不高。不过大尺度和中尺度预测图残差不大，显示出的总体分布特征还是有相当的参考价值。表 和残差平方和 点数残差平方和大尺度.中尺度.小尺度.微尺度.小尺度和微尺度取样时预测值可信，微尺度取样吻合度更高，微尺度预测图比小尺度预测图过渡更为平缓。小尺度取样预测时 虽然较高但是残差平方和反倒是所有尺度里最高的，这应该是在小尺度上随机因素较大造成的，此结论和计算半方差函数时的结论相一致。也和其他研究者对小尺度的重视相一致。这或许和植物竞争主要发生在小尺度上有关。微尺度预测精度很高，这容易被认为是由于微尺度上土壤性质均一性更好，但其实另有原因，这可以从不同尺度取样时变异系数的比较中发现这一点，微尺度取样时其变异系数明显大于中、小尺度，中小尺度由于尺度加大反而降低了变异。说明有产生变异的主要因素有一个最主要的影响尺度，这个尺度明显小于本研究中的中小尺度而和微尺度比较接近。这个尺度其实是产生变化的真正单位。由于微尺度接近产生主要变异的因素作用的尺度而使其预测精度较高，这也是在微尺度以上对具体点预测能力急速降低的真正原因。结 论从变程上看，四个研究尺度下松嫩草地苏打盐碱土壤 值的变程为 、左右和小于 。在这几个范围内存在空间自相关，在这几个范围间不存在空间自相关。也就是说，不同的因素在不同尺度上的效应使 值的变化形成了这几个尺度层次。值的这几个层次的变化和什么因素相关呢？松嫩地区第 卷第 期王欣慰，等：苏打盐碱草地不同研究尺度土壤 值空间变异特征 黑龙江境内的盐碱土呈小斑块状分布，盐土和碱土总面积约为.，每个斑块的平均面积大约为.。这样每个斑块的直径约为 ，和大尺度下的变程十分一致。这些往往被道路、湖泊、河流、村庄和市镇等分割开来的独立的盐碱化斑块自成空间自相关的一个整体。研究地草地退化严重，被碱斑环绕的植物群落呈岛状分布，这些“岛”上的植物形成许多同心圆环状斑块（很多时候这个圆并不不完整），与碱斑距离不同的圆环状斑块上生长着不同的优势植物。假如将生长同一种优势植物的斑块归为一类，则同一类斑块的宽度有相当的差异。在研究地随机设置了一条样线，样线穿过不同种类的斑块环，统计了同一类斑块中的最窄宽度和“岛”的宽度（表）。斑块的平均最窄宽度（）和微尺度的变程（）很接近，这个小于 的范围很可能是一般情况下土壤和植物相互作用并产生可见效应的最小尺度。中尺度和小尺度的变程很接近，约 左右。在研究地地上植被没有发现明显的 左右变化，岛状植物斑块的宽度明显小于这个尺度，有可能是 个岛状植物斑块的性质变化有较紧密的联系而形成直径 左右的 值变化单位。表 植物斑块宽度统计表（括号内为此斑块内优势植物）（）斑块类型斑块数量斑块最小宽度 类型（虎尾草 ）类型（鹅绒委陵菜 寸苔草 ）类型（驴耳风毛菊 ）类型（碱斑 ）类型（萹蓄蓼虎尾草.）类型（蔓委陵菜寸苔草.）类型（车前 ）类型（黄蒿 ）类型（羊草 ）类型（碱地肤 .）类型（碱蓬 ）斑块平均宽度 岛状植物斑块宽度 总体而言，小尺度和微尺度可以方便的得到精确的预测值，大尺度和中尺度取样可以研究 值的总体分布趋势。研究地的盐碱化草地形成约为.的大尺度斑块，每个斑块内约有 个约 的 值变化单位，而苏打盐碱土壤 值和植物之间相互作用主要发生在小于 的尺度上。参考文献 吴东辉，尹文英，卜照义.松嫩草原中度退化草地不同植被恢复方式下土壤线虫的群落特征生态学报，（）：.王志春，李取生，李秀军，等.松嫩平原盐碱化土地治理与农业持续发展对策 中国生态农业学报，（）：.王春裕.试议新型高蛋白植物杂交酸模的开发利用生态学杂志，（）：.徐仁扣，.某些农业措施对土壤酸化的影响农业环境保护，（）：.王政权.地统计学在生态学中的应用北京：科学出版社，.李艳，史舟，徐建明，等.地统计学在土壤科学中的应用及展望水土保持学报，（）：.刘铁军，荣浩，何京丽，等.地质统计学法在土壤风蚀研究中的计算方法探讨 全国水土保持与荒漠化防治及生态修复交流研讨会论文集：.马宁远，王惠卿，张伟，等.基于地统计学的新疆棉田烟粉虱（）危害动态与时空分布生态学报，（）：.王洪峰，王博文，阎秀峰.植物种植对盐碱草地土壤化学性质改良效果的定量化评价中国草地学报，（）：.刘芳，迟耀斌，王智勇，等.列入生态统计指标体系的潜力分析 以北京地区 测算与空间分析为例生态环境学报，（）：.，.，（）：.史舟，李艳，金辉明.基于方差四叉树法的滨海盐土电导率采样布局研究土壤学报，（）：.胡顺军，康绍忠，宋郁东，等.渭干河灌区土壤水盐空间变异性研究水土保持学报，（）：.王志刚，赵永存，廖启林，等.近 年来江苏省土壤 值时空变化及其驱动力生态学报.，（）：.邵学新，黄标，顾志权，等.长三角经济高速发展地区土壤 时空变化及其影响因素矿物岩石地球化学通报，（）：.鲁如坤.土壤农业化学分析方法北京：中国农业科技出版社，.曾艳，张杨珠.地统计学在土壤性状空间变异性研究中的应用湖南农业科学，（）：，.刘璐，宋同清，彭晚霞，等.木论喀斯特自然保护区土壤微生物生物量的空间格局生态学报，（）：.宋同清，彭晚霞，曾馥平，等.喀斯特木论自然保护区旱季土壤水分的空间异质性应用生态学报，（）：.欧阳资文，彭晚霞，宋同清，等.喀斯特峰丛洼地土壤有机质的空间变化及其对干扰的响应应用生态学报，（）：.刘丽，罗承德，雷波，等.蜀南竹海毛竹林土壤生物肥力质量指标与评价土壤，（）：.张利青，彭晚霞，宋同清，等.云贵高原喀斯特坡耕地土壤微生物量、空间分布生态学报，（）：.曹宏杰，王立民，罗春雨，等.三江平原地区土壤有机碳及其组分的空间分布特征生态环境学报，（）：.张忠华，胡刚，祝介东，等.喀斯特森林土壤养分的空间异质性及其对树种分布的影响植物生态学报，（）：.张娜，王希华，郑泽梅，等.浙江天童常绿阔叶林土壤的空间异质性及其与地形的关系应用生态学报，（）：.，.，（）：.刘振国，李镇清.植物群落中物种小尺度空间结构研究植物生态学报，（）：.，.，：中国科学院南京土壤研究所.中国土壤数据库：.，.：.