全 文 ：北川震后滑坡体表层土壤养分的空间变化∗
麦积山１　 赵廷宁２　 郑江坤１∗∗　 史常青２
（ １四川农业大学林学院， 成都 ６１１１３０； ２北京林业大学水土保持学院， 北京 １０００８６）
摘　 要　 以北川擂鼓镇凤凰山震后滑坡体表层土壤为研究对象，结合网格取样和室内分析，
运用地统计学和经典统计学方法对其养分空间变化进行分析．结果表明： 研究区土壤有效磷
（ＡＰ）属高变异强度，其他养分均属中变异强度．ＡＰ 和土壤有机碳（ＳＯＣ）的块基比分别为
２７．９％和 ２８．８％，表现为中等空间相关性；土壤全氮（ＴＮ）、全磷（ＴＰ）、全钾（ＴＫ）、水解性氮
（ＡＮ）和速效钾（ＡＫ）的块基比分别为 ２０．０％、２４．３％、１１．１％、１１．２％和 ２２．７％，均表现为强烈的
空间相关性．土壤 ＴＰ 的变程最大，为 １２３２．７ ｍ；其次为 ＡＮ（５４１．２７ ｍ）、ＴＮ（４６８．３５ ｍ）、ＴＫ
（１３６．０ ｍ）、ＡＫ（１２８．７ ｍ）、ＡＰ（１１６．６ ｍ）、ＳＯＣ（９３．５ ｍ） ．各土壤养分随海拔的变化规律不明
显，而从滑坡区、过渡区至未破坏区则呈明显增加趋势．其中，滑坡区表层土壤 ＴＰ 含量较未破
坏区减少１０．３％，减幅最小，而 ＡＰ 减幅最大，达 ７９． ７％；过渡区各养分含量占未破坏区的
３１．１％～８７．２％，降幅最小和最大的仍为 ＴＰ 和 ＡＰ．地震滑坡扰动是导致土壤养分空间变化的
主要原因．
关键词　 “５·１２”汶川地震； 滑坡体； 土壤养分； 空间变异； 地统计学
文章编号　 １００１－９３３２（２０１５）１２－３５８８－０７　 中图分类号　 Ｓ１５３．６　 文献标识码　 Ａ
Ｓｐａｔｉａｌ ｖａｒｉａｂｉｌｉｔｙ ｏｆ ｓｕｒｆａｃｅ ｓｏｉｌ ｎｕｔｒｉｅｎｔｓ ｉｎ ｔｈｅ ｌａｎｄｓｌｉｄｅ ａｒｅａ ｏｆ Ｂｅｉｃｈｕａｎ Ｃｏｕｎｔｙ， Ｓｏｕｔｈ⁃
ｗｅｓｔ Ｃｈｉｎａ， ａｆｔｅｒ ５·１２ Ｗｅｎｃｈｕａｎ Ｅａｒｔｈｑｕａｋｅ． ＭＡＩ Ｊｉ⁃ｓｈａｎ１， ＺＨＡＯ Ｔｉｎｇ⁃ｎｉｎｇ２， ＺＨＥＮＧ
Ｊｉａｎｇ⁃ｋｕｎ１， ＳＨＩ Ｃｈａｎｇ⁃ｑｉｎｇ２ （ １Ｃｏｌｌｅｇｅ ｏｆ Ｆｏｒｅｓｔｒｙ， Ｓｉｃｈｕａｎ Ａｇｒｉｃｕｌｔｕｒａｌ Ｕｎｉｖｅｒｓｉｔｙ， Ｃｈｅｎｇｄｕ
６１１１３０， Ｃｈｉｎａ； ２Ｃｏｌｌｅｇｅ ｏｆ Ｓｏｉｌ ａｎｄ Ｗａｔｅｒ Ｃｏｎｓｅｒｖａｔｉｏｎ， Ｂｅｉｊｉｎｇ Ｆｏｒｅｓｔｒｙ Ｕｎｉｖｅｒｓｉｔｙ， Ｂｅｉｊｉｎｇ
１０００８６， Ｃｈｉｎａ） ． ⁃Ｃｈｉｎ． Ｊ． Ａｐｐｌ． Ｅｃｏｌ．， ２０１５， ２６（１２）： ３５８８－３５９４．
Ａｂｓｔｒａｃｔ： Ｂａｓｅｄ ｏｎ ｇｒｉｄ ｓａｍｐｌｉｎｇ ａｎｄ ｌａｂｏｒａｔｏｒｙ ａｎａｌｙｓｉｓ， ｓｐａｔｉａｌ ｖａｒｉａｂｉｌｉｔｙ ｏｆ ｓｕｒｆａｃｅ ｓｏｉｌ ｎｕ⁃
ｔｒｉｅｎｔｓ ｗａｓ ａｎａｌｙｚｅｄ ｗｉｔｈ ＧＳ＋ ａｎｄ ｏｔｈｅｒ ｓｔａｔｉｓｔｉｃｓ ｍｅｔｈｏｄｓ ｏｎ ｔｈｅ ｌａｎｄｓｌｉｄｅ ａｒｅａ ｏｆ Ｆｅｎｇｈｕａｎｇ Ｍｏｕｎ⁃
ｔａｉｎ， Ｌｅｉｇｕ Ｔｏｗｎ， Ｂｅｉｃｈｕａｎ Ｃｏｕｎｔｙ． Ｔｈｅ ｒｅｓｕｌｔｓ ｓｈｏｗｅｄ ｔｈａｔ ｅｘｃｅｐｔ ｆｏｒ ｈｉｇｈ ｖａｒｉａｂｉｌｉｔｙ ｏｆ ａｖａｉｌａｂｌｅ
ｐｈｏｓｐｈｏｒｕｓ， ｏｔｈｅｒ ｓｏｉｌ ｎｕｔｒｉｅｎｔｓ ｅｘｈｉｂｉｔｅｄ ｍｏｄｅｒａｔｅ ｖａｒｉａｂｉｌｉｔｙ． Ｔｈｅ ｒａｔｉｏｓ ｏｆ ｎｕｇｇｅｔ ｔｏ ｓｉｌｌ ｏｆ ｔｈｅ ｓｏｉｌ
ａｖａｉｌａｂｌｅ ｐｈｏｓｐｈｏｒｕｓ ａｎｄ ｓｏｉｌ ｏｒｇａｎｉｃ ｃａｒｂｏｎ ｗｅｒｅ ２７．９％ ａｎｄ ２８．８％， ｒｅｓｐｅｃｔｉｖｅｌｙ， ｓｈｏｗｉｎｇ ｍｏｄｅ⁃
ｒａｔｅ ｓｐａｔｉａｌ ｃｏｒｒｅｌａｔｉｏｎ， ｗｈｉｌｅ ｔｈｅ ｒａｔｉｏｓ ｏｆ ｎｕｇｇｅｔ ｔｏ ｓｉｌｌ ｏｆ ｔｈｅ ｔｏｔａｌ ｎｉｔｒｏｇｅｎ （２０．０％）， ｔｏｔａｌ ｐｈｏｓ⁃
ｐｈｏｒｕｓ （２４．３％）， ｔｏｔａｌ ｐｏｔａｓｓｉｕｍ （１１．１％）， ａｖａｉｌａｂｌｅ ｎｉｔｒｏｇｅｎ （１１．２％）， ａｎｄ ａｖａｉｌａｂｌｅ ｐｏｔａｓｓｉｕｍ
（２２．７％） ｓｕｇｇｅｓｔｅｄ ｓｔｒｏｎｇ ｓｐａｔｉａｌ ｃｏｒｒｅｌａｔｉｏｎ． Ｔｈｅ ｔｏｔａｌ ｐｈｏｓｐｈｏｒｕｓ ｈａｄ ｔｈｅ ｍａｘｉｍｕｍ ｒａｎｇｅ
（１２３２．７ ｍ）， ｆｏｌｌｏｗｅｄ ｂｙ ａｖａｉｌａｂｌｅ ｎｉｔｒｏｇｅｎ （５４１．２７ ｍ）， ｔｏｔａｌ ｎｉｔｒｏｇｅｎ （４６８．３５ ｍ）， ｔｏｔａｌ ｐｏｔａｓｓｉ⁃
ｕｍ （１３６．０ ｍ）， ａｖａｉｌａｂｌｅ ｐｏｔａｓｓｉｕｍ （１２８．７ ｍ）， ａｖａｉｌａｂｌｅ ｐｈｏｓｐｈｏｒｕｓ （１１６．６ ｍ）， ａｎｄ ｓｏｉｌ ｏｒｇａ⁃
ｎｉｃ ｃａｒｂｏｎ （９３．５ ｍ）． Ｓｏｉｌ ｎｕｔｒｉｅｎｔｓ ｈａｄ ｎｏ ｓｉｇｎｉｆｉｃａｎｔ ｖａｒｉａｔｉｏｎ ｗｉｔｈ ｔｈｅ ｉｎｃｒｅａｓｅ ｏｆ ａｌｔｉｔｕｄｅ， ｂｕｔ
ｇｒａｄｕａｌｌｙ ｉｎｃｒｅａｓｅｄ ｆｒｏｍ ｔｈｅ ｌａｎｄｓｌｉｄｅ ａｒｅａ， ｔｈｅ ｔｒａｎｓｉｔｉｏｎ ａｒｅａ， ｔｏ ｔｈｅ ｌｉｔｔｌｅ⁃ｉｍｐａｃｔｅｄ ａｒｅａ． Ｔｈｅ ｔｏｔａｌ
ａｎｄ ａｖａｉｌａｂｌｅ ｐｈｏｓｐｈｏｒｕｓ ｃｏｎｔｅｎｔｓ ｏｆ ｔｈｅ ｌａｎｄｓｌｉｄｅ ａｒｅａ ｄｅｃｒｅａｓｅｄ ｂｙ １０．３％ ａｎｄ ７９．７％ ｃｏｍｐａｒｅｄ ｔｏ
ｔｈａｔ ｏｆ ｔｈｅ ｌｉｔｔｌｅ⁃ｉｍｐａｃｔｅｄ ａｒｅａ， ｒｅｓｐｅｃｔｉｖｅｌｙ． Ｔｈｅ ｓｏｉｌ ｎｕｔｒｉｅｎｔ ｃｏｎｔｅｎｔｓ ｉｎ ｔｈｅ ｔｒａｎｓｉｔｉｏｎ ａｒｅａ ａｃｃｏｕｎ⁃
ｔｅｄ ｆｏｒ ３１．１％－８７．２％ ｏｆ ｔｈａｔ ｏｆ ｔｈｅ ｌｉｔｔｌｅ⁃ｉｍｐａｃｔｅｄ ａｒｅａ， ｗｉｔｈ ｔｈｅ ｌｅａｓｔ ａｎｄ ｌａｒｇｅｓｔ ｒｅｄｕｃｔｉｏｎｓ ｉｎ ｔｏ⁃
ｔａｌ ａｎｄ ａｖａｉｌａｂｌｅ ｐｈｏｓｐｈｏｒｕｓ ｃｏｎｔｅｎｔｓ， ｒｅｓｐｅｃｔｉｖｅｌｙ． Ｔｈｅ ｄｉｓｔｕｒｂａｎｃｅ ｂｙ ｅａｒｔｈｑｕａｋｅ ｌａｎｄｓｌｉｄｅ ｗａｓ
ｔｈｅ ｄｏｍｉｎａｎｔ ｒｅａｓｏｎ ｆｏｒ ｔｈｅ ｓｐａｔｉａｌ ｖａｒｉａｂｉｌｉｔｙ ｏｆ ｓｕｒｆａｃｅ ｓｏｉｌ ｎｕｔｒｉｅｎｔｓ．
Ｋｅｙ ｗｏｒｄｓ： “５·１２” Ｗｅｎｃｈｕａｎ ｅａｒｔｈｑｕａｋｅ； ｌａｎｄｓｌｉｄｅ； ｓｏｉｌ ｎｕｔｒｉｅｎｔ； ｓｐａｔｉａｌ ｖａｒｉａｂｉｌｉｔｙ； ｇｅｏｓｔａ⁃
ｔｉｓｔｉｃｓ．
∗林业公益性行业科研专项（２０１１０４１０９）、中国博士后科学基金面上项目（２０１２Ｍ５１１９３８）和四川省高等学校建设项目资助．
∗∗通讯作者． Ｅ⁃ｍａｉｌ： ｊｉａｎｇｋｕｎｚｈｅｎｇ＠ １２６．ｃｏｍ
２０１５⁃０３⁃０５收稿，２０１５⁃０９⁃１１接受．
应 用 生 态 学 报　 ２０１５年 １２月　 第 ２６卷　 第 １２期　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　
Ｃｈｉｎｅｓｅ Ｊｏｕｒｎａｌ ｏｆ Ａｐｐｌｉｅｄ Ｅｃｏｌｏｇｙ， Ｄｅｃ． ２０１５， ２６（１２）： ３５８８－３５９４
　 　 “５·１２”汶川地震诱发了大量的滑坡泥石流等
衍生地质灾害，据统计，大约有 ２００００处各种类型的
滑坡发生［１］ ．滑坡造成大量的土壤流失和植被摧毁，
改变土壤养分空间格局和土壤通气透水等理化性
质，导致滑坡后的坡面土壤质量严重下降、植被资源
衰退甚至灭绝等一系列生态环境问题．目前，对滑坡
的研究多集中在形成机理、治理措施和滑坡体土壤
理化性质等方面［２－４］，而对其土壤养分空间异质性
及其性质特征的研究较少．Ｃｈｅｎｇ 等［５］提出，地震诱
发的滑坡对土壤的养分、容重、孔隙度等理化性质均
有重大影响．因此，开展震后滑坡体土壤养分空间分
布格局的研究不仅是植被恢复的重要方向和途
径［６］，也可为生态恢复治理提供依据．土壤是时空连
续的变异体，不论是在小尺度还是大尺度上均存在
高度的空间变异性［７］ ．对土壤空间异质性的研究不
仅是分析土壤养分和性质空间分布规律的有效方
法，也是揭示空间格局、生态过程和生态功能三者之
间关系的重要途径［８］ ．国内外对土壤养分空间异质
性进行了大量研究［９－１１］，Ｂｏｅｒｎｅｒ 等［１２］研究了不同
扰动程度的 ３个地块中土壤无机氮、有机碳、有效磷
的空间格局，强调量化和认识空间格局可作为植被
恢复和预测植被变化的先决条件． Ｚｈａｎｇ 等［１３］研究
了东北黑土区土壤养分的空间变异特征，提出结构
因素为空间相关性的主要原因．郭晓敏等［１４］研究发
现，毛竹林土壤养分随机变异小于结构性变异，施肥
时长对土壤养分空间均匀分布有调和作用．本研究
以“５·１２”汶川特大地震衍生的北川县擂鼓镇凤凰
山典型滑坡面为研究对象，设置过渡区和未破坏的
柏木林地作为对照，格网采样结合地统计学，系统分
析了土壤养分空间变异特征，旨在为地震灾后植被
恢复提供一定的理论依据．
１　 研究区域与研究方法
１􀆰 １　 研究区概况
研究区属亚热带湿润季风气候，多年平均气温
１５．６ ℃，平均降雨量 １３９９ ｍｍ，降雨集中在 ６—９ 月，
占全年降雨量的 ７１％ ～ ７６％，且空间分布不均．擂鼓
镇全境皆山，土壤质地以砾石土为主，其次为壤土和
粘土，粗骨性强，地质条件差，水文条件复杂，山体岩
石破碎，沟谷谷坡较大，是滑坡崩塌泥石流易发区．
全镇各类森林面积 １００００ ｈｍ２左右，其中天然林
７３３３ ｈｍ２，人工林 ２６６６ ｈｍ２，森林覆盖率 ７６％．采样
区位于 ３１°４６′２４″ Ｎ、１０４°２５′２４″ Ｅ，滑坡前为柏木
（Ｃｐｕｒｅｓｓｕｓ ｆｕｎｅｂｒｉｓ）林，郁闭度高．地震后部分山体发
表 １　 样地基本情况
Ｔａｂｌｅ １　 Ｂａｓｉｃ ｉｎｆｏｒｍａｔｉｏｎ ｏｆ ｔｈｅ ｐｌｏｔｓ
样地类型
Ｐｌｏｔ ｔｙｐｅ
坡度　 　 　 　
Ｓｌｏｐｅ　 　 　 　
（°） 　 　 　 　
土壤颗粒组成
Ｓｏｉｌ ｐａｒｔｉｃｌｅ ｃｏｍｐｏｓｉｔｉｏｎ （％）
粘粒
Ｃｌａｙ
粉粒
Ｓｉｌｔ
砂粒
Ｓａｎｄ
Ⅰ ２７ １０．３ ５．６ ８４．１
Ⅱ ２５ １２．１ ７．７ ８０．２
Ⅲ ２６ １３．１ １０．０ ７６．９
Ⅰ： 滑坡区 Ｌａｎｄｓｌｉｄｅ ａｒｅａ； Ⅱ： 过渡区 Ｔｒａｎｓｉｔｉｏｎ ｚｏｎｅ； Ⅲ： 未破坏区
Ｆａｉｌｕｒｅ ｚｏｎｅ．下同 Ｔｈｅ ｓａｍｅ ｂｅｌｏｗ．
生剧烈滑坡，植被及表层土壤流失殆尽，地表裸露，
坡面几乎全为砾石；２０１０ 和 ２０１１ 年在当地政府和
一些国际组织的资助下，对滑坡体采取了工程固坡
和灌草措施，有效减少了该区的水土流失（表 １，图
２）．其中，样地Ⅰ为滑坡区，２０１１ 年按株行距 ２ ｍ×
４ ｍ人工栽植刺槐幼苗；样地Ⅱ是滑坡区和未破坏
区之间的过渡区域，存在少量面蚀和沟蚀，主要分布
草灌类植物；样地Ⅲ为未破坏区，未受滑坡影响，植
被保持原状．
１􀆰 ２　 样品采集及分析
２０１３年 ７月，利用森林罗盘仪（ＤＱＬ⁃１）和海拔
仪（Ｍａｇｅｌｌａｎ⁃５００Ｅ），基于网格法布置 ４５ 个样点（图
１），然后在每个样点附近再选取 ２ 个样点，分别用
土刀取表层 ２０ ｃｍ×２０ ｃｍ×１０ ｃｍ 砾土混合后带回
实验室，样品风干后，将土壤中的植物落叶、残根、砾
石、动物残体等挑出，过筛研磨后的土壤待测．同时，
在每个样点用手持 ＧＰＳ 定位，并记录下海拔、经纬
度、坡度、砾石裸露率、人为干扰等信息．测定指标包
括：全氮 （ ＴＮ）、全磷 （ ＴＰ ）、全钾 （ ＴＫ）、碱解氮
（ＡＮ）、有效磷（ＡＰ）、速效钾（ＡＫ）、有机碳（ＳＯＣ），
方法参照森林土壤分析方法［１５］ ．
１􀆰 ３　 数据处理
鉴于采样的随机性和人为活动的干扰，本文剔
除样本均值加减 ３ 倍标准差以外的异常值，保持半
方差函数的稳定性．用 Ｋ⁃Ｓ 检验法（显著性水平 α ＝
０．０５）检验数据是否符合正态分布，保证数据符合
Ｋｒｉｇｉｎｇ空间插值要求．地统计学在土壤养分空间变异
的研究中应用广泛，其方法和原理参见文献［１６－１７］．
图 １　 研究区采样点布设
Ｆｉｇ．１　 Ｓｏｉｌ ｓａｍｐｌｉｎｇ ｌｏｃａｔｉｏｎｓ ｉｎ ｔｈｅ ｒｅｓｅａｒｃｈ ａｒｅａ．
９８５３１２期　 　 　 　 　 　 　 　 　 　 　 　 　 　 麦积山等： 北川震后滑坡体表层土壤养分的空间变化　 　 　 　 　 　
本研究采用 ＧＳ＋Ｖｅｒｓｉｏｎ ９．０ 分析半方差函数、块金
值、块基比等参数；土壤养分 Ｋｒｉｇｉｎｇ 等值线绘制在
ＡｒｃＧＩＳ １０．０ 中完成．方差等统计分析在 ＳＰＳＳ １７．０
完成．
块基比［Ｃ０ ／ （Ｃ０＋Ｃ）］为空间相关度，表示随机
部分引起的空间异质性占系统总变异的比例，如果
该比值低，表明其空间异质性主要由结构性因素造
成，相反则说明随机部分在引起空间异质性程度中
起主要作用［１８］ ．以往研究表明，当 Ｃ０ ／ （Ｃ０＋Ｃ）＜２５％
时，变量具有强烈的空间相关性；Ｃ０ ／ （ Ｃ０ ＋ Ｃ）在
２５％～７５％时，为中等空间相关性；而 Ｃ０ ／ （Ｃ０ ＋Ｃ） ＞
７５％时，变量空间相关性很弱［１９］ ．变程（Ａ）反映各变
量之间只有在变程范围内才具有空间相关性，在变
程之外则被认为是相互独立的［２０］ ．半变异函数的决
定系数（Ｒ２）越大、残差（ＲＳＳ）越小，说明函数曲线
的拟合程度越好［２１］ ．
２　 结果与分析
２􀆰 １　 土壤养分的统计性分析
由表 ２可知，整个采样区土壤 ＴＮ、ＴＰ、ＴＫ、ＡＮ、
ＡＰ、ＡＫ 和 ＳＯＣ 含量范围分别为 ０． １２ ～ ５􀆰 ９０
ｇ·ｋｇ－１、０．０１ ～ １． ４０ ｇ·ｋｇ－１、０． ０７ ～ ０􀆰 ６４ ｇ·ｋｇ－１、
１５􀆰 ６１～４６９．６４ ｍｇ·ｋｇ－１、０．１９～１．６３ ｍｇ·ｋｇ－１、１６．９０～
８７９􀆰 ７０ ｍｇ·ｋｇ－１和 １６．００～１４２．５２ ｇ·ｋｇ－１ ．土壤 ＴＮ、
ＳＯＣ、ＡＫ含量较为丰富，其平均含量分别为 １􀆰 ８６
ｇ·ｋｇ－１、４４．５７ ｇ·ｋｇ－１、１５６．７０ ｍｇ·ｋｇ－１；土壤 ＡＮ
含量为 １１５．５５ ｍｇ·ｋｇ－１，属于中等水平；土壤 ＴＰ、
ＴＫ和 ＡＰ 属极缺乏水平，其平均值分别为 ０􀆰 ３６
ｇ·ｋｇ－１、０．３４ ｇ·ｋｇ－１和 １． ２２ ｍｇ·ｋｇ－１ ．变异系数
（ＣＶ）能够反映数据的离散程度，可划分为：弱变异
性（＜１０％）；中等变异性（１０％ ～ １００％）；强变异性
（＞１００％） ［２２］ ．土壤 ＡＰ 的变异系数超过了 １００％，属
于高变异强度，其余土壤养分含量的变异系数介于
３８􀆰 ２％～９０．２％之间，均属于中等变异强度，各土壤
养分变异程度大小依次为：ＡＰ ＞ＡＫ＞ＡＮ＞ＴＮ＞ＴＰ ＞
ＳＯＣ＞ＴＫ．结合偏度、峰度和 Ｋ⁃Ｓ 检验，土壤 ＴＫ 和
ＡＫ符合正态分布，其他土壤养分数据经对数转换
后均服从正态分布，满足地统计学分析的条件．
２􀆰 ２　 土壤养分的空间结构特征
利用地统计学方法分析各土壤养分的空间变化
特征．由表 ３可知，土壤 ＴＰ 的最优拟合为指数模型，
土壤 ＴＫ、ＡＰ、ＡＫ 和 ＳＯＣ 的最优拟合为球状模型，
ＴＮ和 ＡＮ的最优拟合模型为高斯模型．各土壤养分
指标拟合模型的拟合程度均较高，都能较充分地反
映土壤养分的空间特征，其中，土壤 ＴＰ 的决定系数
最小，表明它在空间变化上相对不稳定．
研究区土壤养分变程排序为：ＳＯＣ＜ＡＰ＜ＡＫ＜ＴＫ
＜ＴＮ＜ＡＮ＜ＴＰ，其中，ＴＰ 的变程最大，达到 １２３２．７ ｍ，
其空间连续性较好，需要更小尺度的研究．土壤 ＴＫ、
ＡＰ、ＡＫ和 ＳＯＣ 的空间连续性尺度较相近，分别为
１３６、１１６．６、１２８􀆰 ７和 ９３．５ ｍ，可视为其生态过程在相
同尺度上起作用，而 ＳＯＣ的变程最小，仅为 ９３􀆰 ５ ｍ，
其空间连续性较差．块金值反映采样尺度和试验误
差引起的空间变异，表 ３ 中全部土壤养分的块金值
均较小且都为正值，即存在较小的采样与试验误
差［２３］ ．土壤 ＡＰ 和 ＳＯＣ 的块基比分别为 ２７． ９％和
２８􀆰 ８％，为中等强度空间自相关．土壤 ＴＮ、ＴＰ、ＴＫ、
ＡＮ和 ＡＫ的块基比均＜２５％，表现出强烈的空间自
相关性，说明其变异主要是因为结构性因素的影响，
如地形、土壤类型、群落结构、土壤母质等［２４］ ．综上
所述，空间自相关部分的异质性在土壤养分的空间
异质性中占主要部分，土壤养分具有重要的空间结
构，变程代表它的尺度大小．
２􀆰 ３　 土壤养分空间分布格局
由图 ２ 可知，整个研究区域左侧的土壤 ＴＮ、
ＴＫ、ＡＮ、ＡＰ、ＡＫ、ＳＯＣ含量均较右侧低，总体呈现出
表 ２　 土壤养分描述性统计特征
Ｔａｂｌｅ ２　 Ｄｅｓｃｒｉｐｔｉｖｅ ｓｔａｔｉｓｔｉｃｓ ｏｆ ｓｏｉｌ ｎｕｔｒｉｅｎｔｓ ｉｎ ｔｈｅ ｓｔｕｄｙ ｓｏｉｌ （ｎ＝４５）
土壤养分
Ｓｏｉｌ ｎｕｔｒｉｅｎｔ
平均值
Ｍｅａｎ
标准差
Ｓｔａｎｄａｒｄ ｄｅｖｉａｔｉｏｎ
变异系数
ＣＶ （％）
范围
Ｒａｎｇｅ
偏度
Ｓｋｅｗ
峰度
Ｋｕｒｔｏｓｉｓ
Ｋ⁃Ｓ检验
Ｋ⁃Ｓ ｔｅｓｔ
ＴＮ （ｇ·ｋｇ－１） １．８６ １．４４ ７７．４ ０．１２～５．９０ １．１７ ０．４９ ０．０１９
ＴＰ （ｇ·ｋｇ－１） ０．３６ ０．２７ ７５．０ ０．０１～１．４０ １．９０ ４．１４ ０．０３２
ＴＫ （ｇ·ｋｇ－１） ０．３４ ０．１３ ３８．２ ０．０７～０．６４ －０．０４ ０．１５ ０．９９５
ＡＮ （ｍｇ·ｋｇ－１） １１５．５５ ９５．３０ ８２．５ １５．６１～４６９．６４ １．７９ ３．５３ ０．０４７
ＡＰ （ｍｇ·ｋｇ－１） １．２２ ２．４０ １９６．７ ０．１９～１．６３ ５．９８ ３８．０４ ０
ＡＫ （ｍｇ·ｋｇ－１） １５６．７０ １４１．３８ ９０．２ １６．９０～８７９．７０ ３．３４ １５．２５ ０．０５９
ＳＯＣ （ｇ·ｋｇ－１） ４４．５７ ２９．３０ ６５．７ １６．００～１４２．５２ １．７５ ２．７６ ０．０１０
ＴＮ： 全氮 Ｔｏｔａｌ ｎｉｔｒｏｇｅｎ； ＴＰ： 全磷 Ｔｏｔａｌ ｐｈｏｓｐｈｏｒｕｓ； ＴＫ： 全钾 Ｔｏｔａｌ ｐｏｔａｓｓｉｕｍ； ＡＮ： 水解性氮 Ｈｙｄｒｏｌｙｚａｂｌｅ ｎｉｔｒｏｇｅｎ； ＡＰ： 有效磷 Ａｖａｉｌａｂｌｅ ｐｈｏｓ⁃
ｐｈｏｒｕｓ； ＡＫ： 速效钾 Ａｖａｉｌａｂｌｅ ｐｏｔａｓｓｉｕｍ； ＳＯＣ： 土壤有机碳 Ｓｏｉｌ ｏｒａｇｎｉｃ ｃａｒｂａｎ．下同 Ｔｈｅ ｓａｍｅ ｂｅｌｏｗ．
０９５３ 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 应　 用　 生　 态　 学　 报　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 ２６卷
表 ３　 土壤养分指标半方差模型及其参数
Ｔａｂｌｅ ３　 Ｓｅｍｉｖａｒｉｏｇｒａｍ ｔｈｅｏｒｅｔｉｃａｌ ｍｏｄｅｌｓ ａｎｄ ｐａｒａｍｅｔｅｒｓ ｏｆ ｓｏｉｌ ｎｕｔｒｉｅｎｔｓ ｉｎｄｉｃｅｓ
土壤养分
Ｓｏｉｌ ｎｕｔｒｉｅｎｔ
模型
Ｍｏｄｅｌ
块金值
Ｃ０
基台值
Ｃ０＋Ｃ
块金值 ／基台值
Ｃ０ ／ （Ｃ０＋Ｃ）
变程
Ａ （ｍ）
决定系数
Ｒ
残差
ＲＳＳ
ＴＮ （ｇ·ｋｇ－１） 高斯模型 Ｇａｕｓｓ ｍｏｄｅｌ ０．４７６ ２．３８３ ０．２００ ４６８．３５ ０．８２３ ０．０４０
ＴＰ （ｇ·ｋｇ－１） 指数模型 Ｅｘｐｏｎｅｎｔｉａｌ ｍｏｄｅｌ ０．０５０ ０．２０６ ０．２４３ １２３２．７０ ０．２８９ ＜０．００３
ＴＫ （ｇ·ｋｇ－１） 球状模型 Ｓｐｈｅｒｉｃａｌ ｍｏｄｅｌ ０．００２ ０．０１８ ０．１１１ １３６．００ ０．８６７ ＜０．００１
ＡＮ （ｍｇ·ｋｇ－１） 高斯模型 Ｇａｕｓｓ ｍｏｄｅｌ ０．２９１ ２．５９２ ０．１１２ ５４１．２７ ０．７８８ ０．０４５
ＡＰ （ｍｇ·ｋｇ－１） 球状模型 Ｓｐｈｅｒｉｃａｌ ｍｏｄｅｌ ０．１６８ ０．６０２ ０．２７９ １１６．６０ ０．４００ ０．１１７
ＡＫ （ｍｇ·ｋｇ－１） 球状模型 Ｓｐｈｅｒｉｃａｌ ｍｏｄｅｌ ０．１１７ ０．５１５ ０．２２７ １２８．７０ ０．８８７ ＜０．００９
ＳＯＣ （ｇ·ｋｇ－１） 球状模型 Ｓｐｈｅｒｉｃａｌ ｍｏｄｅｌ ０．０９０ ０．３１４ ０．２８８ ９３．５０ ０．３９２ ０．０２５
自左向右逐渐增大的趋势．土壤 ＴＮ含量分布呈斑块
状，高值区在研究区的右侧，且表现出随海拔降低而
下降的现象，主要是因为上坡的植被结构为柏木林
图 ２　 研究区土壤养分的空间分布
Ｆｉｇ．２　 Ｓｐａｔｉａｌ ｄｉｓｔｒｉｂｕｔｉｏｎ ｏｆ ｓｏｉｌ ｎｕｔｒｉｅｎｔｓ ｉｎ ｔｈｅ ｒｅｓｅａｒｃｈ ａｒｅａ．
ＴＮ： 全氮 Ｔｏｔａｌ ｎｉｔｒｏｇｅｎ； ＴＰ： 全磷 Ｔｏｔａｌ ｐｈｏｓｐｈｏｒｕｓ； ＴＫ： 全钾 Ｔｏｔａｌ
ｐｏｔａｓｓｉｕｍ； ＡＮ： 水解性氮 Ｈｙｄｒｏｌｙｚａｂｌｅ ｎｉｔｒｏｇｅｎ； ＡＰ： 有效磷 Ａｖａｉｌａ⁃
ｂｌｅ ｐｈｏｓｐｈｏｒｕｓ； ＡＫ： 速效钾 Ａｖａｉｌａｂｌｅ ｐｏｔａｓｓｉｕｍ； ＳＯＣ： 土壤有机碳
Ｓｏｉｌ ｏｒａｇｎｉｃ ｃａｒｂａｎ．
和草、灌植物，而下坡主要是草、灌植物，且植被覆盖
度上坡高于下坡，更有利于养分的积累，所以海拔越
低 ＴＮ含量越低．土壤 ＴＰ 含量的高值集中出现于整
个研究区中的偏左侧，在该区左侧呈由上到下递增
分布，但在水平分布上无明显递变规律．土壤 ＴＫ 含
量在 ０．０１ ～ ０．２８ ｇ·ｋｇ－１的面积占总面积的 ５０％左
右，从左至右呈现高低值交错变化的分布规律，其最
小值与 ＡＫ含量最小值均出现在研究区左上角．土
壤 ＡＮ值呈从右向左递减的带状分布，含量最高处
主要集中在研究区右侧，与 ＳＯＣ 分布规律相似，这
是由于水解性氮与有机质之间通常存在显著的正相
关关系．土壤 ＡＰ 值分布较为均匀，变化缓和，从左
至右呈较为规律的递增趋势，其含量在 ０．０９ ～ ０􀆰 ４２
ｍｇ·ｋｇ－１区间的面积约占总面积的 ２０％．土壤 ＡＫ空
间异质性明显，在研究区右侧表现出随海拔的降低
而递减的规律，而在左侧则无明显变化规律，其最高
值位于研究区的右上角．土壤 ＳＯＣ 含量的分布规律
较为明显，为从右上到左下递减的带状分布．
２􀆰 ４　 滑坡对土壤养分含量的影响
由表 ４ 可知，研究区土壤养分含量整体表现为
从样地Ⅰ到Ⅲ递增趋势，而随海拔升高无明显规律．
样地Ⅰ中，ＴＮ、ＴＰ、ＴＫ、ＡＮ、ＡＰ、ＡＫ和 ＳＯＣ含量依次
为样地Ⅲ的 ３０．９％、８９．７％、５０．０％、２８．４％、２０．３％、
２７．３％和 ４０．９％，样地Ⅱ的相应养分含量依次占样
地Ⅲ的 ６４．６％、８７．２％、７１．７％、７０．５％、３１．１％、４９．８％
和 ６２． ２％．土壤 ＴＮ 含量在海拔梯度 ６１０ ～ ６３０ ｍ
（Ｈ２）和 ６５０～ ６７０ ｍ（Ｈ４）上从样地Ⅰ到Ⅲ递增，分
别增加了 ６．２和 ３．８倍，且在 Ｈ２ 梯度上差异显著．土
壤 ＴＰ 值未表现出明显变化规律．ＴＫ和 ＡＫ含量在 ５
个海拔梯度均呈现出从样地Ⅰ到Ⅲ递增现象．ＡＮ值
在海拔梯度 Ｈ２、６３０ ～ ６５０ ｍ（Ｈ３）和 Ｈ４ 上从样地Ⅰ
到Ⅲ递增．ＡＰ 含量除在海拔梯度 Ｈ４以外，其余海拔
梯度均从样地Ⅰ到Ⅲ递增． ＳＯＣ 含量在海拔梯度
５８０～６１０ ｍ（Ｈ１）、Ｈ２ 、Ｈ４上从样地Ⅰ到Ⅲ呈现出明
１９５３１２期　 　 　 　 　 　 　 　 　 　 　 　 　 　 麦积山等： 北川震后滑坡体表层土壤养分的空间变化　 　 　 　 　 　
表 ４　 不同海拔和样地类型的土壤养分特征
Ｔａｂｌｅ ４　 Ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ｏｆ ｓｏｉｌ ｎｕｔｒｉｅｎｔｓ ｉｎ ｄｉｆｆｅｒｅｎｔ ａｌｔｉｔｕｄｅｓ ａｎｄ ｐｌｏｔｓ
土壤养分
Ｓｏｉｌ ｎｕｔｒｉｅｎｔ
样地类型
Ｓａｍｐｌｅ
ｔｙｐｅ
海 拔 Ａｌｔｉｔｕｄｅ （ｍ）
５８０～６１０ ６１０～６３０ ６３０～６５０ ６５０～６７０ ６７０～７１０ 均值 Ｍｅａｎ
ＴＮ Ⅰ ０．６７Ａａ ０．７６Ｂａ １．１３Ａａ ０．９９Ａｂ ０．８４Ｂａ ０．８８Ａ
（ｇ·ｋｇ－１） Ⅱ １．３５Ａａ １．３４Ｂａ ０．９０Ａａ ２．６８Ａａ ２．９３Ａａ １．８４Ｂ
Ⅲ ０．８６Ａａ ４．７１Ａｂ ２．７７Ａｂ ３．７８Ａｂ ２．１２ＡＢｂ ２．８５Ｃ
均值 Ｍｅａｎ ０．９６ａ ２．２７ａ １．６ａ ２．４８ａ １．９６ａ －
ＴＰ Ⅰ ０．４５Ａａ ０．４７Ａａ ０．１９Ａａ ０．３６Ａａ ０．３０Ａａ ０．３５Ａ
（ｇ·ｋｇ－１） Ⅱ ０．２４Ａａ ０．３１Ａａ ０．２１Ａａ ０．３０Ａａ ０．６６Ａａ ０．３４Ａ
Ⅲ ０．４１Ａａ ０．５９Ａａ ０．３２Ａａ ０．３１Ａａ ０．３３Ａａ ０．３９Ａ
均值 Ｍｅａｎ ０．３６ａ ０．４６ａ ０．２４ａ ０．３２ａ ０．４３ａ －
ＴＫ Ⅰ ０．２７Ａａ ０．２３Ａａ ０．２４Ａａ ０．２５Ａａ ０．１８Ａａ ０．２３Ａ
（ｇ·ｋｇ－１） Ⅱ ０．３５Ａａ ０．３８Ｂａ ０．３３Ａａ ０．２６Ａａ ０．３２ＡＢａ ０．３３Ｂ
Ⅲ ０．５０Ａａ ０．４６Ｂａ ０．４６Ｂａ ０．４６Ａａ ０．４３Ｂａ ０．４６Ｃ
均值 Ｍｅａｎ ０．３７ａ ０．３６ａ ０．３４ａ ０．３２ａ ０．３１ａ －
ＡＮ Ⅰ ４５．０４Ａａ ３４．２４Ａａ ６２．５４Ａａ ５６．６８Ａａ ４９．１２Ａａ ４９．５２Ａ
（ｍｇ·ｋｇ－１） Ⅱ １３７．７８Ａａ １３８．９６Ｂａ ８４．６０Ａａ １１６．５０Ａａ １３６．７５Ａａ １２２．９２Ａ
Ⅲ ８５．２６Ａａ ２３１．９６Ｂａ １６０．８８Ｂａ ２８４．２４Ｂａ １０９．０８Ａａ １７４．２８Ｂ
均值 Ｍｅａｎ ８９．３５ａ １３５．０５ａ １０２．６７ａ １５２．４７ａ ９８．３２ａ －
ＡＰ Ⅰ ０．５５Ａａ ０．４７Ａａ ０．３４Ａａ ０．５ＡＢａ ０．６Ａａ ０．４９Ａ
（ｍｇ·ｋｇ－１） Ⅱ ０．６９Ａａ ０．７９Ａａ ０．７４Ｂａ ０．９０Ａａ ０．６５Ａａ ０．７５Ａ
Ⅲ ２．２Ａａ １．３３Ａａ １．３９Ｂａ ０．８４Ａａ ６．２８Ｂａ ２．４１Ｂ
均值 Ｍｅａｎ １．１５ａ ０．８６ａ ０．８２ａ ０．７５ａ ２．５１ａ －
ＡＫ Ⅰ ７０．１０Ａａ ７９．１Ａａ ８２．５Ａａ ８０．０１Ａａ ５０．８Ａａ ７２．５１Ａ
（ｍｇ·ｋｇ－１） Ⅱ １１６．０２Ａａ １２５．８Ａａ １２１．１０Ａａ １５６．２２ＡＢａ １４１．６ＡＢａ １３２．１５Ａ
Ⅲ １２８．５Ａａ ２７０．４Ｂａ ２２２．９２Ａａ ３０９．３２Ｂａ ３９６Ｂａ ２６５．４３Ｂ
均值 Ｍｅａｎ １０４．８７ａ １５８．４３ａ １４２．１７ａ １８１．８５ａ １９６．１３ａ －
ＳＯＣ Ⅰ ２６．９Ａａ ２２．５２Ａａ ３１．６Ａａ ２７．５Ａａ ２６．１Ａａ ２６．９４Ａ
（ｇ·ｋｇ－１） Ⅱ ３１．１Ａａ ２９．１２Ａａ ２９．３Ａａ ５３．２Ａａ ６１．９Ｂａ ４０．９２Ａ
Ⅲ ３１．７Ａａ ９４．１４Ｂｂ ５６．３Ａａｂ ９７．６Ａｂ ４９．５ＡＢａｂ ６５．８４Ｂ
均值 Ｍｅａｎ ２９．９ａ ４８．５９ａ ３９．０７ａ ５９．４３ａ ４５．８３ａ －
不同大写字母和小写字母分别表示不同样地类型和不同海拔段间差异显著（Ｐ＜０．０５） Ｄｉｆｆｅｒｅｎｔ ｕｐｐｅｒｃａｓｅ ｌｅｔｔｅｒ （ ｓ） ｗｉｔｈｉｎ ｅａｃｈ ｃｏｌｕｍｎ ｉｎｄｉｃａｔｅｄ
ｓｉｇｎｉｆｉｃａｎｔｌｙ ｄｉｆｆｅｒｅｎｔ ａｔ ０．０５ ｌｅｖｅｌ ａｍｏｎｇ ｄｉｆｆｅｒｅｎｔ ｓａｍｐｌｅ ｔｙｐｅｓ， ｄｉｆｆｅｒｅｎｔ ｌｏｗｅｒｃａｓｅ ｌｅｔｔｅｒ （ ｓ） ｗｉｔｈｉｎ ｅａｃｈ ｒｏｗ ｉｎｄｉｃａｔｅｄ ｓｉｇｎｉｆｉｃａｎｔｌｙ ｄｉｆｆｅｒｅｎｔ ａｔ ０．０５
ｌｅｖｅｌ ａｍｏｎｇ ｄｉｆｆｅｒｅｎｔ ａｌｔｉｔｕｄｅ ｇｒａｄｉｅｎｔｓ．
显递增规律．比较各样地养分均值可以得出，土壤
ＴＮ和 ＴＫ含量分别在样地Ⅰ、Ⅱ和Ⅲ之间均差异性
显著．但土壤 ＴＰ 含量则在各样地之间均无显著性差
异．除 ＴＮ、ＴＫ和 ＴＰ 以外，其余各养分含量在样地Ⅰ
和Ⅱ之间均不存在显著性差异，但都与样地Ⅲ之间
差异显著．而全部养分均值在各海拔段之间差异均
不显著．
３　 讨　 　 论
研究区土壤养分整体含量水平较低，变异程度
较大，土壤 ＡＰ 为高强度变异，其余养分都为中等强
度变异．主要由于该区山体岩石破碎，水文条件复
杂，加之受地震和滑坡的影响，改变了地质地貌和植
被分布等．土壤养分有限或缺乏直接导致植物对土
壤有效养分的吸收与利用，影响滑坡体上植被的恢
复，进而不利于水土保持，加大养分流失．这将严重
制约该区生态环境修复．因此，应结合植物措施和布
设工程措施来改善局部环境，为植被生长创造有利
的条件．
土壤各养分指标的变异函数模型拟合均较好，
只有土壤 ＴＰ 的变异函数模型拟合度较低，空间变
异结构较差．土壤养分空间异质性是随机因素和结
构因素综合作用的结果，随机因素促使其空间相关
性减弱，结构性因素加强其空间相关性［２５］ ．土壤 ＡＰ
和 ＳＯＣ表现为中等强度空间自相关性，表明它是受
结构性因子和随机因子共同作用的结果．土壤 ＴＮ、
ＴＰ、ＴＫ、ＡＮ和 ＳＯＣ均为强烈空间自相关，其变异主
要是因为地震和滑坡改变了地形、土壤类型重新分
布，以及植被结构或母质分布等结构性因素．研究区
土壤 ＴＰ 的变程最大，具有较好的空间连续性，表明
环境因素在较大尺度上控制着 ＴＰ 的空间异质
性［２６］，其余各养分变程的大小顺序为：ＳＯＣ＜ＡＰ＜ＡＫ
２９５３ 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 应　 用　 生　 态　 学　 报　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 ２６卷
＜ＴＫ＜ＴＮ＜ＡＮ．土壤 ＴＫ、ＡＰ、ＡＫ、ＳＯＣ的变程较小，分
别为 １３６、１１６．６、１２８．７和 ９３．５ ｍ，说明这几种养分的
空间连续性尺度较小．
从 Ｋｒｉｇｉｎｇ 插值图可以看出，研究区除土壤 ＴＰ
以外，所有土壤养分含量从滑坡区至未破坏区均表
现为递增趋势．这与吴聪等［４］研究结果较一致，即滑
坡将导致各土壤养分含量大幅降低，但不同区域同
一养分含量的降幅存在差异．造成本研究土壤养分
这种分布规律的原因是滑坡对其坡面地表扰动极
大，影响了土壤理化性质，改变了坡面微地貌，大量
表层土壤及其养分损失，土层浅薄，植被只有人工栽
种的刺槐幼林，水土和养分流失严重，生态环境质量
大大降低．过渡区受滑坡影响较小，土壤中石砾含量
较大，但坡面表层土壤含量仍大于滑坡区，且坡面分
布有较多草、灌植物，在一定程度上能防止土壤颗粒
和养分的流失，故而各养分含量高于滑坡区．未破坏
区植被覆盖率高，能有效防止表层土壤物质及其养
分的流失，同时，林下大量的枯落物为有机质积累创
造了良好的条件，因此该样地区域养分含量最高．而
土壤 ＴＰ 特殊的分布规律是由于土壤 ＴＰ 主要来源
于基岩风化，并通过植物的表聚作用在表层土壤积
累［２７］ ．研究区中偏左的位置受滑坡破坏生产大量碎
石，通过风化作用使供给植物吸收的磷来源增多，加
之此处草本植物稀少，减少了磷素消耗，所以 ＴＰ 含
量没有表现出从左至右增加的规律．海拔对土壤养
分影响显著．张伟等［６］研究表明，土壤养分随海拔增
加存在“倒挂现象”，而本研究中土壤养分含量在海
拔范围内无明显变化规律．这一方面因为滑坡对研
究区生态系统造成了极大扰动，致使土壤养分空间
分布格局复杂多样；另一方面，因滑坡发生时间短，
土壤养分分布规律与海拔之间的关系还未明显形
成．该地震滑坡体地形破碎，地表破坏严重，土壤养
分变异强度较大，空间变异明显，各土壤养分都具有
明显的空间自相关格局．因此，探明土壤养分空间分
布特征能有效地指导滑坡区的生态恢复重建．
参考文献
［１］　 Ｘｕ Ｃ （许　 冲）， Ｄａｉ Ｆ⁃Ｃ （戴福初）， Ｘｕ Ｘ⁃Ｗ （徐锡
伟）． Ｗｅｎｃｈｕａｎ Ｅａｒｔｈｑｕａｋｅ⁃ｉｎｄｕｃｅｄ Ｌａｎｄｓｌｉｄｅｓ： Ａｎ
ｏｖｅｒｖｉｅｗ． Ｇｅｏｌｏｇｉｃａｌ Ｒｅｖｉｅｗ （地质评论）， ２０１０， ５６
（６）： ８６０－８７４ （ｉｎ Ｃｈｉｎｅｓｅ）
［２］　 Ｊｉｎ Ｚ⁃Ｆ （金增芳）． Ｙｕｓｈｕ ｓｅｉｓｍｉｃ ｌａｎｄｓｌｉｄｅ ｃｈａｒａｃｔｅｒｉｓ⁃
ｔｉｃｓ ａｎｄ ｆｏｒｍａｔｉｏｎ ｍｅｃｈａｎｉｓｍ ａｎｄ ｉｔｓ ｃｏｎｔｒｏｌ ｍｅａｓｕｒｅｓ．
Ｂｕｉｌｄｉｎｇ Ｓｔｒｕｃｔｕｒｅ （建筑结构）， ２０１３， ４３（ ｓｕｐｐｌ．２）：
６８６－６８９ （ｉｎ Ｃｈｉｎｅｓｅ）
［３］　 Ｚｈａｎｇ Ｄ⁃Ｇ （张德罡）． Ｅｆｆｅｃｔ ｏｆ ｄｅｆｏｒｅｓｔａｔｉｏｎ ａｎｄ ｌａｎｄ⁃
ｓｌｉｄｅ ｏｎ ｓｏｉｌ ｏｆ Ｒｈｏｄｏｄｅｎｄｒｏｎ ｓｈｒｕｂ ｅｃｏｓｙｓｔｅｍ ｉｎ ｔｈｅ
ｅａｓｔｅｒｎ Ｑｉｌｉａｎ Ｍｏｕｎｔａｉｎｓ． Ａｃｔａ Ｐｒａｔａｃｕｌｔｕｒａｅ Ｓｉｎｉｃａ （草
业学报）， ２００２， １１（３）： ７２－７５ （ｉｎ Ｃｈｉｎｅｓｅ）
［４］　 Ｗｕ Ｃ （吴　 聪）， Ｗａｎｇ Ｊ⁃Ｎ （王金牛）， Ｌｕ Ｔ （卢 　
涛）， ｅｔ ａｌ． Ｅｆｆｅｃｔ ｏｆ Ｗｅｎｃｈｕａｎ ｅａｒｔｈｑｕａｋｅ ｏｎ ｓｏｉｌ ｐｈｙ⁃
ｓｉｃａｌ ａｎｄ ｃｈｅｍｉｃａｌ ｐｒｏｐｅｒｔｉｅｓ ｏｆ ｔｈｅ Ｌｏｎｇｍｅｎ Ｍｏｕｎｔａｉｎ，
ｓｏｕｔｈｗｅｓｔｅｒｎ Ｃｈｉｎａ． Ｃｈｉｎｅｓｅ Ｊｏｕｒｎａｌ ｏｆ Ａｐｐｌｉｅｄ ＆ Ｅｎｖｉ⁃
ｒｏｎｍｅｎｔａｌ Ｂｉｏｌｏｇｙ （应用与环境生物学报）， ２０１２， １８
（６）： ９１１－９１６ （ｉｎ Ｃｈｉｎｅｓｅ）
［５］　 Ｃｈｅｎｇ Ｓ， Ｙａｎｇ Ｇ， Ｙｕ Ｈ， ｅｔ ａｌ． Ｉｍｐａｃｔｓ ｏｆ Ｗｅｎｃｈｕａｎ
Ｅａｒｔｈｑｕａｋｅ⁃ｉｎｄｕｃｅｄ ｌａｎｄｓｌｉｄｅｓ ｏｎ ｓｏｉｌ ｐｈｙｓｉｃａｌ ｐｒｏｐｅｒ⁃
ｔｉｅｓ ａｎｄ ｔｒｅｅ ｇｒｏｗｔｈ． Ｅｃｏｌｏｇｉｃａｌ Ｉｎｄｉｃａｔｏｒｓ， ２０１２， １５：
２６３－２７０
［６］　 Ｚｈａｎｇ Ｗ （张　 伟）， Ｃｈｅｎ Ｈ⁃Ｓ （陈洪松）， Ｗａｎｇ Ｋ⁃Ｌ
（王克林）， ｅｔ ａｌ． Ｓｐａｔｉａｌ ｖａｒｉａｂｉｌｉｔｙ ｏｆ ｓｏｉｌ ｎｕｔｒｉｅｎｔｓ ｏｎ
ｈｉｌｌｓｌｏｐｅ ｉｎ ｔｙｐｉｃａｌ ｋａｒｓｔ ｐｅａｋ⁃ｃｌｕｓｔｅｒ ｄｅｐｒｅｓｓｉｏｎ ａｒｅａｓ．
Ｔｒａｎｓａｃｔｉｏｎｓ ｏｆ ｔｈｅ Ｃｈｉｎｅｓｅ Ｓｏｃｉｅｔｙ ｏｆ Ａｇｒｉｃｕｌｔｕｒａｌ Ｅｎｇｉ⁃
ｎｅｅｒｉｎｇ （农业工程学报）， ２００８， ２４（１）： ６８－７３ （ ｉｎ
Ｃｈｉｎｅｓｅ）
［７］　 Ｘｉｏｎｇ Ｈ⁃Ｆ （熊汉锋）， Ｗａｎｇ Ｙ⁃Ｈ （王运华）． Ｓｐａｔｉａｌ
ｖａｒｉａｂｉｌｉｔｙ ｏｆ ｓｏｉｌ ｎｕｔｒｉｅｎｔｓ ｉｎ ｗｅｔｌａｎｄ ｏｆ Ｌｉａｎｇｚｉ Ｌａｋｅ．
Ｐｌａｎｔ Ｎｕｔｒｉｔｉｏｎ ａｎｄ Ｆｅｒｔｉｌｉｚｅｒ Ｓｃｉｅｎｃｅ （植物营养与肥料
学报）， ２００５， １１（５）： ５８４－５８９ （ｉｎ Ｃｈｉｎｅｓｅ）
［８］　 Ｅｔｔｅｍａ ＣＨ， Ｗａｒｄｌｅ ＤＡ． Ｓｐａｔｉａｌ ｓｏｉｌ ｅｃｏｌｏｇｙ． Ｔｒｅｎｄｓ ｉｎ
Ｅｃｏｌｏｇｙ ＆ Ｅｖｏｌｕｔｉｏｎ， ２００２， １７： １７７－１８３
［９］　 Ｚｈａｎｇ Ｙ⁃Ｍ （张玉铭）， Ｍａｏ Ｒ⁃Ｚ （毛任钊）， Ｈｕ Ｃ⁃Ｓ
（胡春胜）， ｅｔ ａｌ． Ｓｐａｔｉａｌ ｖａｒｉａｂｉｌｉｔｙ ｏｆ ｆａｒｍｌａｎｄ ｓｏｉｌ ｎｕ⁃
ｔｒｉｅｎｔｓ ａｔ Ｔａｉｈａｎｇ ｐｉｅｄｍｏｎｔ． Ｃｈｉｎｅｓｅ Ｊｏｕｒｎａｌ ｏｆ Ａｐｐｌｉｅｄ
Ｅｃｏｌｏｇｙ （应用生态学报）， ２００４， １５（１１）： ２０４９－２０５４
（ｉｎ Ｃｈｉｎｅｓｅ）
［１０］　 Ｊｉａｎｇ Ｙ （姜　 勇）， Ｚｈａｎｇ Ｙ⁃Ｇ （张玉革）， Ｌｉａｎｇ Ｗ⁃Ｊ
（梁文举）， ｅｔ ａｌ． Ｓｐａｔｉａｌ ｖａｒｉａｂｉｌｉｔｙ ｏｆ ｓｏｉｌ ｎｕｔｒｉｅｎｔｓ ｉｎ
ｃｕｌｔｉｖａｔｅｄ ｓｕｒｆａｃｅ ｓｏｉｌ Ｓｕｊｉａｔｕｎ Ｄｉｓｔｒｉｃｔ， Ｓｈｅｎｙａｎｇ Ｃｉｔｙ．
Ｃｈｉｎｅｓｅ Ｊｏｕｒｎａｌ ｏｆ Ａｐｐｌｉｅｄ Ｅｃｏｌｏｇｙ （应用生态学报），
２００３， １４（１０）： １６７３－１６７６ （ｉｎ Ｃｈｉｎｅｓｅ）
［１１］　 Ｓｏｎｇ Ｘ （宋　 轩）， Ｌｉ Ｌ⁃Ｄ （李立东）， Ｋｏｕ Ｃ⁃Ｌ （寇
长林）， ｅｔ ａｌ． Ｓｏｉｌ ｎｕｔｒｉｅｎｔ ｄｉｓｔｒｉｂｕｔｉｏｎ ａｎｄ ｉｔｓ ｒｅｌａｔｉｏｎｓ
ｗｉｔｈ ｔｏｐｏｇｒａｐｈｙ ｉｎ Ｈｕａｎｇｓｈｕｉ Ｒｉｖｅｒ ｄｒａｉｎａｇｅ ｂａｓｉｎ． Ｃｈｉ⁃
ｎｅｓｅ Ｊｏｕｒｎａｌ ｏｆ Ａｐｐｌｉｅｄ Ｅｃｏｌｏｇｙ （应用生态学报），
２０１１， ２２（１２）： ３１６３－３１６８ （ｉｎ Ｃｈｉｎｅｓｅ）
［１２］　 Ｂｏｅｒｎｅｒ ＲＥＪ， Ｓｃｈｅｒｚｅｒ ＡＪ， Ｂｒｉｎｋｍａｎ ＪＡ． Ｓｐａｔｉａｌ ｐａｔ⁃
ｔｅｒｎｓ ｏｆ ｉｎｏｒｇａｎｉｃ Ｎ， Ｐ ａｖａｉｌａｂｉｌｉｔｙ， ａｎｄ ｏｒｇａｎｉｃ Ｃ ｉｎ ｒｅ⁃
ｌａｔｉｏｎ ｔｏ ｓｏｉｌ ｄｉｓｔｕｒｂａｎｃｅ： Ａ ｃｈｒｏｎｏｓｅｑｕｅｎｃｅ ａｎａｌｙｓｉｓ．
Ａｐｐｌｉｅｄ Ｓｏｉｌ Ｅｃｏｌｏｇｙ， １９９８， ７： １５９－１７７
［１３］　 Ｚｈａｎｇ ＸＹ， Ｓｕｉ ＹＹ， Ｚｈａｏ ＸＤ， ｅｔ ａｌ． Ｓｐａｔｉａｌ ｖａｒｉａｂｉｌｉｔｙ
ｏｆ ｎｕｔｒｉｅｎｔ ｐｒｏｐｅｒｔｉｅｓ ｉｎ ｂｌａｃｋ ｓｏｉｌ ｏｆ ｎｏｒｔｈｅａｓｔ Ｃｈｉｎａ．
Ｐｅｄｏｓｐｈｅｒｅ， ２００７， １７： １９－２９
［１４］　 Ｇｕｏ Ｘ⁃Ｍ （郭晓敏）， Ｎｉｕ Ｄ⁃Ｋ （牛德奎）， Ｇｕｏ Ｘ （郭
熙）， ｅｔ ａｌ． Ｓｔｕｄｙ ｏｎ ｓｐａｔｉａｌ ｖａｒｉａｂｉｌｉｔｙ ｏｆ ｓｏｉｌ ｎｕｔｒｉｅｎｔｓ
ｆｏｒ Ｍｏｓｏ ｂａｍｂｏｏ ｆｏｒｅｓｔ ｉｎ Ｆｅｎｇｘｉｎ． Ｐｌａｎｔ Ｎｕｔｒｉｔｉｏｎ ａｎｄ
Ｆｅｒｔｉｌｉｚｅｒ Ｓｃｉｅｎｃｅ （植物营养与肥料学报）， ２００６， １２
（３）： ４２０－４２５ （ｉｎ Ｃｈｉｎｅｓｅ）
［１５］ 　 Ｚｈａｎｇ Ｗ⁃Ｒ （张万儒）． Ｆｏｒｅｓｔ Ｓｏｉｌ Ａｎａｌｙｓｉｓ Ｍｅｔｈｏｄ．
Ｂｅｉｊｉｎｇ： Ｃｈｉｎａ Ｓｔａｎｄａｒｄｓ Ｐｒｅｓｓ， ２０００ （ｉｎ Ｃｈｉｎｅｓｅ）
［１６］　 Ｗｅｂｓｔｅｒ Ｒ， Ｏｌｉｖｅｒ ＭＡ． Ｇｅｏｓｔａｔｉｓｔｉｃｓ ｆｏｒ　 Ｅｎｖｉｒｏｎｍｅｎｔａｌ
Ｓｃｉｅｎｔｉｓｔｓ． Ｃｈｉｃｈｅｓｔｅｒ： Ｊｏｈｎ Ｗｉｌｅｙ ａｎｄ Ｓｏｎｓ， ２００７
［１７］　 Ｌｉ Ｌ⁃Ｌ （李亮亮）， Ｙｉ Ｙ⁃Ｌ （依艳丽）， Ｌｉｎｇ Ｇ⁃Ｘ （凌国
３９５３１２期　 　 　 　 　 　 　 　 　 　 　 　 　 　 麦积山等： 北川震后滑坡体表层土壤养分的空间变化　 　 　 　 　 　
鑫）， ｅｔ ａｌ． Ｕｔｉｌｉｚａｔｉｏｎ ｏｆ ｇｅｏｓｔａｔｉｓｔｉｃｓ ｉｎ ｓｏｉｌ ｓｐａｔｉａｌ ｖａｒｉａ⁃
ｂｉｌｉｔｙ． Ｃｈｉｎｅｓｅ Ｊｏｕｒｎａｌ ｏｆ Ｓｏｉｌ Ｓｃｉｅｎｃｅ （土壤通报），
２００５， ３６（２）： ２６５－２６８ （ｉｎ Ｃｈｉｎｅｓｅ）
［１８］　 Ｇａｎ Ｈ⁃Ｈ （甘海华）， Ｐｅｎｇ Ｌ⁃Ｙ （彭凌云）． Ｓｐａｔｉａｌ ｖａ⁃
ｒｉａｂｉｌｉｔｙ ｏｆ ｎｕｔｒｉｅｎｔｓ ｉｎ ｃｕｌｔｉｖａｔｅｄ ｓｏｉｌｓ ｏｆ Ｘｉｎｈｕｉ Ｄｉｓ⁃
ｔｒｉｃｔ， Ｊｉａｎｇｍｅｎ Ｃｉｔｙ． Ｃｈｉｎｅｓｅ Ｊｏｕｒｎａｌ ｏｆ Ａｐｐｌｉｅｄ Ｅｃｏｌｏｇｙ
（应用生态学报）， ２００５， １６（８）： １４３７－１４４２ （ｉｎ Ｃｈｉ⁃
ｎｅｓｅ）
［１９］　 Ｌｉｕ Ｌ （刘　 璐）， Ｚｅｎｇ Ｆ⁃Ｐ （曾馥平）， Ｓｏｎｇ Ｔ⁃Ｑ （宋
同清）， ｅｔ ａｌ． Ｓｐａｔｉａｌ ｈｅｔｅｒｏｇｅｎｅｉｔｙ ｏｆ ｓｏｉｌ ｎｕｔｒｉｅｎｔｓ ｉｎ
Ｋａｒｓｔ ａｒｅａ’ｓ Ｍｕｌｕｎ Ｎａｔｉｏｎａｌ Ｎａｔｕｒｅ Ｒｅｓｅｒｖｅ． Ｃｈｉｎｅｓｅ
Ｊｏｕｒｎａｌ ｏｆ Ａｐｐｌｉｅｄ Ｅｃｏｌｏｇｙ （应用生态学报）， ２０１０， ２１
（７）： １６６７－１６７３ （ｉｎ Ｃｈｉｎｅｓｅ）
［２０］　 Ｗａｎｇ Ｊ （王 　 军）， Ｆｕ Ｂ⁃Ｊ （傅伯杰）， Ｑｉｕ Ｙ （邱 　
扬）， ｅｔ ａｌ． Ｓｐａｔｉａｌ ｈｅｔｅｒｏｇｅｎｅｉｔｙ ｏｆ ｓｏｉｌ ｎｕｔｒｉｅｎｔｓ ｉｎ
ｓｍａｌｌ ｃａｔｃｈｍｅｎｔ ｏｆ ｔｈｅ Ｌｏｅｓｓ Ｐｌａｔｅａｕ． Ａｃｔａ Ｅｃｏｌｏｇｉｃａ
Ｓｉｎｉｃａ （生态学报）， ２００２， ２２（８）： １１７３－ １１７８ （ ｉｎ
Ｃｈｉｎｅｓｅ）
［２１］　 Ｆａｎ Ｆ⁃Ｊ （范夫静）， Ｓｏｎｇ Ｔ⁃Ｑ （宋同清）， Ｈｕａｎｇ Ｇ⁃Ｑ
（黄国勤）， ｅｔ ａｌ． Ｃｈａｒａｃｔｅｒｉｓｔｉｃｓ ｏｆ ｓｐａｔｉａｌ ｖａｒｉａｔｉｏｎ ｏｆ
ｓｏｉｌ ｎｕｔｒｉｅｎｔｓ ｉｎ ｓｌｏｐｉｎｇ ｆｉｅｌｄ ｉｎ ａ ｇｏｒｇｅ ｋａｒｓｔ ｒｅｇｉｏｎ，
ｓｏｕｔｈｗｅｓｔ Ｃｈｉｎａ． Ｃｈｉｎｅｓｅ Ｊｏｕｒｎａｌ ｏｆ Ａｐｐｌｉｅｄ Ｅｃｏｌｏｇｙ （应
用生态学报）， ２０１４， ２５（１）： ９２－９８ （ｉｎ Ｃｈｉｎｅｓｅ）
［２２］　 Ｈｕ Ｋ⁃Ｌ （胡克林）， Ｌｉ Ｂ⁃Ｇ （李保国）， Ｌｉｎ Ｑ⁃Ｍ （林
启美）， ｅｔ ａｌ． Ｓｐａｔｉａｌ ｖａｒｉａｂｉｌｉｔｙ ｏｆ ｓｏｉｌ ｎｕｔｒｉｅｎｔ ｉｎ ｗｈｅａｔ
ｆｉｅｌｄ． Ｔｒａｎｓａｃｔｉｏｎｓ ｏｆ ｔｈｅ Ｃｈｉｎｅｓｅ Ｓｏｃｉｅｔｙ ｏｆ Ａｇｒｉｃｕｌｔｕｒａｌ
Ｅｎｇｉｎｅｅｒｉｎｇ （农业工程学报）， １９９９， １５（３）： ３３－３８
（ｉｎ Ｃｈｉｎｅｓｅ）
［２３］　 Ｓｕ Ｓ⁃Ｊ （苏松锦）， Ｌｉｕ Ｊ⁃Ｆ （刘金福）， Ｈｅ Ｚ⁃Ｓ （何中
声）， ｅｔ ａｌ． Ｔｈｅ ｓｐａｔｉａｌ ｈｅｔｅｒｏｇｅｎｅｉｔｙ ｏｆ ｓｏｉｌ ｎｕｔｒｉｅｎｔｓ ｉｎ
ａ ｍｉｄ⁃ｓｕｂｔｒｏｐｉｃａｌ Ｃａｓｔａｎｏｐｓｉｓ ｋａｗａｋａｍｉｉ ｎａｔｕｒａｌ ｆｏｒｅｓｔ．
Ａｃｔａ Ｅｃｏｌｏｇｉｃａ Ｓｉｎｉｃａ （生态学报）， ２０１２， ３２（１８）：
５６７３－５６８２ （ｉｎ Ｃｈｉｎｅｓｅ）
［２４］　 Ｌｉａｎｇ Ｚ⁃Ｌ （梁中龙）， Ｙｕａｎ Ｚ⁃Ｙ （袁中友）， Ｌｉｎ Ｘ⁃Ｔ
（林兴通）， ｅｔ ａｌ． Ｓｐａｔｉａｌ ｖａｒｉａｂｉｌｉｔｙ ｏｆ ｎｕｔｒｉｅｎｔｓ ｉｎ
Ｇｕａｎｇｚｈｏｕ ｓｕｂｕｒｂ． Ｃｈｉｎｅｓｅ Ｊｏｕｒｎａｌ ｏｆ Ｓｏｉｌ Ｓｃｉｅｎｃｅ （土
壤通报）， ２００６， ３７（３）： ４１７－４２１ （ｉｎ Ｃｈｉｎｅｓｅ）
［２５］　 Ｇｕｏ Ｘ⁃Ｄ （郭旭东）， Ｆｕ Ｂ⁃Ｊ （傅伯杰）， Ｍａ Ｋ⁃Ｍ （马
克明）， ｅｔ ａｌ． Ｓｐａｔｉａｌ ｖａｒｉａｂｉｌｉｔｙ ｏｆ ｓｏｉｌ ｎｕｔｒｉｅｎｔｓ ｂａｓｅｄ
ｏｎ ｇｅｏｓｔａｔｉｓｔｉｃｓ ｃｏｍｂｉｎｅｄ ｗｉｔｈ ＧＩＳ： Ａ ｃａｓｅ ｓｔｕｄｙ ｉｎ
Ｚｕｎｇｈｕａ Ｃｉｔｙ ｏｆ Ｈｅｂｅｉ Ｐｒｏｖｉｎｃｅ． Ｃｈｉｎｅｓｅ Ｊｏｕｒｎａｌ ｏｆ Ａｐ⁃
ｐｌｉｅｄ Ｅｃｏｌｏｇｙ （应用生态学报）， ２０００， １１（４）： ５５７－
１５６３ （ｉｎ Ｃｈｉｎｅｓｅ）
［２６］　 Ｚｈａｎｇ Ｗ （张　 伟）， Ｌｉｕ Ｓ⁃Ｊ （刘淑娟）， Ｙｅ Ｙ⁃Ｙ （叶
莹莹）， ｅｔ ａｌ． Ｓｐａｔｉａｌ ｖａｒｉａｂｉｌｉｔｙ ｏｆ ｓｏｉｌ ｎｕｔｒｉｅｎｔｓ ａｎｄ ｉｎ⁃
ｆｌｕｅｎｃｉｎｇ ｆａｃｔｏｒｓ ｉｎ ｔｙｐｉｃａｌ ｋａｒｓｔ ｖｉｒｇｉｎ ｆｏｒｅｓｔ． Ｔｒａｎｓａｃ⁃
ｔｉｏｎｓ ｏｆ ｔｈｅ Ｃｈｉｎｅｓｅ Ｓｏｃｉｅｔｙ ｏｆ Ａｇｒｉｃｕｌｔｕｒａｌ Ｅｎｇｉｎｅｅｒｉｎｇ
（农业工程学报）， ２０１３， ２９（１）： ９３－ １０１ （ ｉｎ Ｃｈｉ⁃
ｎｅｓｅ）
［２７］　 Ｌｉｕ Ｃ⁃Ｑ （刘丛强）． Ｂｉｏｇｅｏｃｈｅｍｉｃａｌ Ｐｒｏｃｅｓｓｅｓ ａｎｄ Ｓｕｒ⁃
ｆａｃｅ Ｍａｔｅｒｉａｌ Ｃｉｒｃｕｌａｔｉｏｎ： Ｔｈｅ Ｓｏｕｔｈｗｅｓｔ Ｋａｒｓｔ Ｓｏｉｌ⁃
Ｖｅｇｅｔａｔｉｏｎ Ｓｙｓｔｅｍ Ｃｙｃｌｅ． Ｂｅｉｊｉｎｇ： Ｓｃｉｅｎｃｅ Ｐｒｅｓｓ， ２００９
（ｉｎ Ｃｈｉｎｅｓｅ）
作者简介　 麦积山，男，１９９２ 年生，硕士研究生．主要从事水
土保持研究． Ｅ⁃ｍａｉｌ： ｍａｉｊｉｓｈａｎ６６６＠ １６３．ｃｏｍ
责任编辑　 杨　 弘
４９５３ 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 应　 用　 生　 态　 学　 报　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 ２６卷