全 文 ：书黄土丘陵区撂荒群落土壤速效养分
空间变异性研究
隋媛媛１，杜峰１，２，张兴昌１，２
（１．西北农林科技大学资源与环境学院，陕西 杨凌７１２１００；２．中国科学院水土保持与生态环境研究中心，陕西 杨凌７１２１００）
摘要：采用经典统计学与地统计学相结合的方法，对黄土丘陵区５块不同撂荒年限样地速效养分在群落尺度上的
空间分布特征进行了分析，结果表明，１）与黄土高原地区整体养分水平相比，撂荒样地土壤养分贫瘠，０～２０和
２０～４０ｃｍ土层速效养分含量在群落演替中期较高，并且０～２０ｃｍ养分含量明显高于２０～４０ｃｍ；２）速效养分的空
间变异主要为结构性变异，撂荒６～７年样地土壤表层的速效钾和硝态氮以及撂荒２５年样地２０～４０ｃｍ硝态氮具
有中等空间相关性，其他撂荒群落土壤０～２０和２０～４０ｃｍ的速效养分均具有高度空间相关性；３）随植被群落撂
荒演替年限的增加，０～２０ｃｍ土层速效钾和硝态氮空间变异性在群落演替前期样地中最大，２０～４０ｃｍ土层速效
钾和硝态氮空间变异性分别在群落演替中期和前、中期样地中最大。０～２０ｃｍ土层速效磷和铵态氮空间变异性在
群落演替中期样地中最大，２０～４０ｃｍ土层速效磷和铵态氮空间变异性分别在群落演替前期和中期样地中最大。
在不同撂荒年限样地中，４种速效养分具有不同空间自相关性和空间变异性，说明控制这几种养分空间分布的生态
学主导过程可能有差异。
关键词：黄土丘陵区；撂荒群落；速效养分；空间变异性
中图分类号：Ｓ１５８．３；Ｘ１７１．４ 文献标识码：Ａ 文章编号：１００４５７５９（２０１１）０２００７６０９
土壤是时空连续的变异体，具有高度的空间异质性［１，２］。一些学者［３，４］利用土壤学家Ｆｉｓｈｅｒ所创立的传统统
计方法对土壤性状的空间分布进行了定量描述，但由于此方法要求样本中的取值相互独立，给土壤空间变异的研
究带来了很大的局限性。２０世纪６０年代 Ｍａｔｈｅｒｏｎ创立了地统计学理论［５］，为土壤空间变异性研究提供了有效
的方法，掀起了国外学者对土壤空间变异性研究的热潮［６８］。目前，土壤空间变异性研究已成为土壤科学的重要
研究内容和热点［９１１］。在土壤属性的空间变异研究中，因土壤养分是土壤肥力的重要指标，能够提供植物生长所
必需的营养元素，与植物群落的类型、分布和动态有着十分密切的关系，因而土壤养分的空间变异研究较多［１２１９］。
土壤养分的空间分布是多种生态学过程综合影响的结果，如降水的淋洗（移）、植物的吸收与表层富集等，因
而土壤养分的空间异质性会受到群落主要植物的分布格局、生态型和养分吸收利用偏好等的影响。有研究表明
随着植被群落的演替，土壤养分含量有所提高［２０］，演替过程中因群落组成和主要植物空间分布格局的变化，土壤
养分的空间分布也会发生相应的变化。初玉等［２１］研究表明，小叶锦鸡儿（犆犪狉犪犵犪狀犪犿犻犮狉狅狆犺狔犾犾犪）灌丛土壤有机
质和全氮空间自相关显著；张璞进等［２２］对藏锦鸡儿（犆犪狉犪犵犪狀犪狋犻犫犲狋犻犮犪）群落土壤资源空间异质性分析研究表明，
藏锦鸡儿在发育过程中形成的“沃岛效应”能够引起有机质和全磷的空间异质性分布。Ｔｓｅｇａｙｅ和 Ｈｉｌ［２３］的研究
发现，植物生长对养分的吸收利用可以影响并加剧土壤养分空间变异；焦峰等［２４］研究了中大尺度上黄土丘陵区
退耕地土壤养分变异特征，结果表明，随植被群落演替恢复，土壤养分先降低再增加并最后趋于稳定，土壤养分的
变异性在前期和中期较大；陈伏生和曾德慧［２５］对科尔沁退化草场土壤养分的空间分析表明，土壤有机质、全氮、
全磷和速效钾都具有明显的空间结构特征。黄土丘陵地区坡耕地比例高［２６］，水土流失严重，土壤养分的贫瘠是
造成土地和植被退化的主要原因，分析它们的空间异质性，理解其空间格局，对植被群落的恢复过程具有重要的
意义。
７６－８４
２０１１年４月
草 业 学 报
ＡＣＴＡＰＲＡＴＡＣＵＬＴＵＲＡＥＳＩＮＩＣＡ
第２０卷 第２期
Ｖｏｌ．２０，Ｎｏ．２
收稿日期：２０１００６０４；改回日期：２０１００９０６
基金项目：中国科学院“西部之光”资助项目（２００７ＺＤ０２）和中国博士后基金资助项目（２００７０４１１１４６）资助。
作者简介：隋媛媛（１９８６），女，吉林东丰人，在读硕士。Ｅｍａｉｌ：ｓｕｉｙｕａｎ＠ｎｗｓｕａｆ．ｅｄｕ．ｃｎ
通讯作者。Ｅｍａｉｌ：ｄｕｆｅｎｇ＠ｍｓ．ｉｓｗｃ．ａｃ．ｃｎ
本研究以黄土丘陵区不同撂荒年限典型群落为研究对象，利用传统统计学与地统计学理论相结合的方法，分
析群落尺度上撂荒样地土壤速效养分变化情况及其空间分布异质性，据此推测影响撂荒地不同形态土壤养分的
生态学过程，为黄土高原地区植被恢复和生态建设的可持续发展提供科学依据。
１ 材料与方法
１．１ 样地自然概况
研究区位于陕北黄土高原丘陵区得安塞县高桥乡，地理位置Ｅ１０９°１２′３２″，Ｎ３６°３８′５２″。年均日照时数
２３００～２５７０ｈ，年均气温７．７～１０．６℃，无霜期为１４２～１７５ｄ，干燥度１．２，≥０℃年积温３７３３．１℃，≥１０℃年积
温３１７０．３℃，年均降水量４９０．５～６６３．３ｍｍ。植被类型为森林草原地带，土壤为黄土母质发育成的黄绵土，轻
壤。研究区退耕地主要优势种有猪毛蒿（犃狉狋犲犿犻狊犻犪狊犮狅狆犪狉犻犪）、铁杆蒿（犃狉狋犲犿犻狊犻犪狊犪犮狉狅狉狌犿）、茭蒿（犃狉狋犲犿犻狊犻犪
犵犻狉犪犾犱犻犻）、长芒草（犛狋犻狆犪犫狌狀犵犲犪狀犪）、达乌里胡枝子（犔犲狊狆犲犱犲狕犪犱犪狌狉犻犮犪）和白羊草（犅狅狋犺狉犻狅犮犺犾狅犪犻狊犮犺犪犲犿狌犿）等。
１．２ 样地调查及采样方法
首先在野外全面踏查的基础上，根据植物群落组成、结构和对当地居民的访问调查结果，选定无人为干扰或
人为干扰较少的５块不同撂荒年限的典型群落作为调查样地，所选样地耕作时间长达５０年以上。２００９年８月，
按照１ｍ×１ｍ网格法进行规则取样，取样层次为０～２０和２０～４０ｃｍ两层，样地面积大小为８ｍ×８ｍ。由于网
格法取样间隔距离单一，规则取样后随机补充７个采样点。每块样地实际采样１４２个，即（６４＋７）×２个，５块样
地共采集７１０个土样。同年９月对各样地进行群落调查，同时使用便携式ＧＰＳ测定各样地经纬度、坡向、坡度和
海拔（表１）。
表１ 调查样地基本情况
犜犪犫犾犲１ 犛狌犿犿犪狉狔狅犳犻狀狏犲狊狋犻犵犪狋犻狅狀狊犻狋犲
样地编号
Ｐｌｏｔ
Ｎｏ．
撂荒年限
Ａｂａｎｄｏｎｅｄ
ａｇｅｓ（年Ｙｅａｒｓ）
经纬度
Ｌｏｎｇｉｔｕｄｅａｎｄ
ｌａｔｉｔｕｄｅ
坡向
Ｓｌｏｐｅｘｐｏｓｕｒｅ
坡度
Ｓｌｏｐｅ
ｄｅｇｒｅｅ（°）
海拔
Ａｌｔｉｔｕｄｅ
（ｍ）
群落类型
Ｃｏｍｍｕｎｉｔｙ
ｔｙｐｅｓ
Ａ ３ １０９°１１′３８″Ｅ ３６°３９′３９″Ｎ 阳坡Ｓｕｎｎｙ ５ １２４０ 猪毛蒿＋狗尾草犃．狊犮狅狆犪狉犻犪＋犛犲狋犪狉犻犪狏犻狉犻犱犻狊
Ｂ ３～４ １０９°１１′５５″Ｅ ３６°３９′２４″Ｎ 阴坡Ｓｈａｄｙ ３ １２９０ 猪毛蒿＋长芒草犃．狊犮狅狆犪狉犻犪＋犛．犫狌狀犵犲犪狀犪
Ｃ ６～７ １０９°１１′５５″Ｅ ３６°３９′２７″Ｎ 半阳坡 Ｈａｌｆｓｕｎｎｙ ２３ １２８０ 长芒草＋达乌里胡枝子犛．犫狌狀犵犲犪狀犪＋犔．犱犪狌狉犻犮犪
Ｄ １４～１５ １０９°１１′５１″Ｅ ３６°３９′２２″Ｎ 阳坡Ｓｕｎｎｙ ２５ １２６０ 铁杆蒿＋达乌里胡枝子犃．狊犪犮狉狅狉狌犿＋犔．犱犪狌狉犻犮犪
Ｅ ２５ １０９°１１′５３″Ｅ ３６°３９′４２″Ｎ 半阳坡 Ｈａｌｆｓｕｎｎｙ １０ １２７０ 达乌里胡枝子＋白羊草犔．犱犪狌狉犻犮犪＋犅．犻狊犮犺犪犲犿狌犿
１．３ 测定项目与方法
土壤样品采集完成后放在干燥、阴凉处自然风干，然后装入密封袋带回实验室，研磨过１ｍｍ筛后立即对
０～２０和２０～４０ｃｍ两土层共７１０个土样进行养分指标分析测定，每个土样做２个重复，从土样采集到速效养分
测定全部过程在１个月内完成。土壤养分指标测定均采用常规土壤化学分析法［２７］，速效钾以１０∶１的水土比用
ＮＨ４ＯＡＣ浸提后，采用原子吸收光谱法测定；速效磷按照１０∶１的水土比用０．５ｍｏｌ／ＬＮａＨＣＯ３ 浸提，采用钼锑
抗比色法测定；铵态氮和硝态氮以５∶１的水土比用 ＫＣｌ浸提，连续流动分析仪测定，两者相加之和为土壤速
效氮。
１．４ 数据处理
本研究采用ＳＰＳＳ１６．０计算各养分变量的描述性统计特征值并对数据进行正态分布检验；利用ＧＳ＋７．０地
统计学软件进行半方差函数及相关参数分析，并使用Ｏｒｉｇｉｎ绘图软件作图。
１．５ 半方差函数
半方差函数是地统计学中定量研究和解释土壤空间异质性的工具函数，首先通过网格取样法得到土壤特性
含量值，根据半方差函数计算公式（１）得到实验半方差函数值［５］，并将其与半方差函数理论模型拟合，依据决定系
７７第２０卷第２期 草业学报２０１１年
数最大原则选取最佳拟合模型，并且得出半方差函数的重要参数，进行土壤养分空间异质性分析。
γ（犺）＝ １２犖（犺）∑
犖（犺）
犻＝１
［犣（狓犻）－犣（狓犻＋犺）］２ （１）
式中，犺为空间分隔距离，犖（犺）为以犺为间距的所取样点的成对数目；犣（狓犻）和犣（狓犻＋犺）分别为空间点狓犻和狓犻＋
犺的养分测定值；γ（犺）为半方差函数。
１．６ 半方差函数模型及参数
图１ 半方差图的结构示意图
犉犻犵．１ 犛狋狉狌犮狋狌狉犲犱犻犪犵狉犪犿狅犳狊犲犿犻狏犪狉犻犪狀犮犲
以变异函数γ（犺）为Ｙ轴，采样间隔ｈ为Ｘ轴，可
绘成变异函数曲线图（图１），根据变异函数曲线的变
化趋势可以拟合成不同的理论模型，如球状模型
（ｓｐｈｅｒｉｃａｌ）、高斯模型（ｇａｕｓｓｉａｎ）、指数模型（ｅｘｐｏｎｅｎ
ｔｉａｌ）。通过半方差函数的模型拟合得到半方差函数极
为重要的参数值，其中块金值Ｎｕｇｇｅｔ（Ｃ０）是间距为０
时的半方差，表示由随机因素引起的变异；Ｃ为结构方
差，表示由系统因素引起的变异；基台值Ｓｉｌ （Ｃ＋Ｃ０）
表示包括随机性变异和空间结构性变异在内的总变
异，基台值越高，说明系统总的空间变异性越高；结构
方差与基台值之比Ｃ／（Ｃ＋Ｃ０）表示结构性变异占总
变异的比值，该比值越大说明结构性空间变异越强，如
果该比值接近０，说明区域化变量在整个尺度上具有
恒定的变异，并且该比值还可以揭示区域化变量的空间相关程度，按照Ｃａｍｂａｒｄｅｌａ等［１２］的划分标准，当Ｃ／（Ｃ＋
Ｃ０）＞７５％，说明变量具有强烈的空间相关性，当Ｃ／（Ｃ＋Ｃ０）＜２５％时，空间相关性较弱，当Ｃ／（Ｃ＋Ｃ０）处于２５％
～７５％时，说明具有中等空间相关性；变程Ｒａｎｇｅ（ａ）为最大相关距离，指在某种观测尺度下，空间相关性的作用
范围，当变量测定点间距小于该值时，说明变量具有空间自相关性，反之则不存在。
２ 结果与分析
２．１ 黄土丘陵区撂荒群落土壤速效养分含量分布及变异情况
统计平均值既能说明土壤特性的大小，又能描述特征数据的集中趋势，结果显示（表２），５块研究样地０～２０
ｃｍ的速效钾含量为６８．５３～９１．８３ｍｇ／ｋｇ，２０～４０ｃｍ的含量为５６．７４～７９．０９ｍｇ／ｋｇ，随着撂荒群落弃耕年限
的增长，０～２０和２０～４０ｃｍ土壤速效钾含量具有相同的变化趋势，均先降低后升高并趋于稳定，在撂荒１４～１５
年样地含量最大；０～２０ｃｍ土壤速效磷含量为０．５４～２．５３ｍｇ／ｋｇ，２０～４０ｃｍ的含量为０．３２～１．９７ｍｇ／ｋｇ，０～
２０和２０～４０ｃｍ土壤速效磷含量随撂荒年限增长，先逐渐升高，并在撂荒１４～１５年样地出现１个小回落，在撂
荒２５年样地含量最大；０～２０ｃｍ硝态氮含量为１．５４～２．７２ｍｇ／ｋｇ，２０～４０ｃｍ含量为１．３５～２．５２ｍｇ／ｋｇ，土壤
０～２０ｃｍ硝态氮含量随撂荒年限变化而上下波动，在撂荒１４～１５年样地含量最大，２０～４０ｃｍ硝态氮在撂荒
３～４年样地含量最低，然后随撂荒年限延长逐渐升高，在撂荒２５年样地含量最高；０～２０ｃｍ氨态氮含量为
２．４０～４．３４ｍｇ／ｋｇ，２０～４０ｃｍ含量为２．２４～４．５１ｍｇ／ｋｇ，０～２０和２０～４０ｃｍ铵态氮随撂荒年限延长均先增
高，在撂荒１４～１５年样地回落到最低值后又升高，但０～２０ｃｍ含量在撂荒２５年样地最大，２０～４０ｃｍ含量在撂
荒６～７年样地最大。在垂直方向上，硝态氮含量在撂荒２５年样地表现为０～２０ｃｍ高于２０～４０ｃｍ，铵态氮在
撂荒６～７年样地也为２０～４０ｃｍ含量高于０～２０ｃｍ含量，其他速效养分含量在不同撂荒年限群落均为土壤０～
２０ｃｍ含量高于２０～４０ｃｍ含量。变异系数是描述土壤特性参数空间变异性程度的指标，依据Ｎｉｅｌｓｅｎ分级标
准［２８］，当变异系数≤１０％时为弱变异性，当１０％＜变异系数＜１００％时为中等变异性，当变异系数≥１００％时为强
变异性。２０～４０ｃｍ速效磷在撂荒３年样地的变异系数为１０５．８％＞１００％（表２），属于强变异，０～２０ｃｍ硝态氮
在撂荒３年样地的变异系数为１３８．４％＞１００％，表现为强变异性，而其他撂荒年限群落土壤速效养分空间分布
均为中等变异。
８７ ＡＣＴＡＰＲＡＴＡＣＵＬＴＵＲＡＥＳＩＮＩＣＡ（２０１１） Ｖｏｌ．２０，Ｎｏ．２
表２ 不同撂荒年限群落土壤速效养分含量统计特征值
犜犪犫犾犲２ 犛狋犪狋犻狊狋犻犮狊狅犳狊狅犻犾犪狏犪犻犾犪犫犾犲狀狌狋狉犻犲狀狋狊犻狀狅犾犱犳犻犲犾犱犮狅犿犿狌狀犻狋犻犲狊狑犻狋犺犱犻犳犳犲狉犲狀狋犪犫犪狀犱狅狀犿犲狀狋犪犵犲狊
分析项目
Ｉｔｅｍｓ
样地编号
ＰｌｏｔＮｏ．
平均值 Ｍｅａｎ
０～２０
ｃｍ
２０～４０
ｃｍ
最大值 Ｍａｘ
０～２０
ｃｍ
２０～４０
ｃｍ
最小值 Ｍｉｎ
０～２０
ｃｍ
２０～４０
ｃｍ
标准差ＳＤ
０～２０
ｃｍ
２０～４０
ｃｍ
变异系数ＣＶ（％）
０～２０
ｃｍ
２０～４０
ｃｍ
分布类型 Ｄｉｓｔｒｉｂｕｔｉｏｎｔｙｐｅ
０～２０
ｃｍ
２０～４０
ｃｍ
速效Ｋ Ａ ８３．６１ ６９．１３ １０７．２０ ８７．００ ５９．１０ ４７．５０ １０．２１ ７．０６ １２．２１ １０．２２ Ｎ Ｎ
ＡｖａｉｌａｂｌｅＫ Ｂ ６８．５３ ５６．７４ １０５．００ ８２．４０ ５０．７０ ４７．９０ ８．９８ ６．０３ １３．１０ １０．６３ Ｎ Ｎ
（ｍｇ／ｋｇ） Ｃ ７４．０７ ７１．０７ １１１．６０１１７．５０ ４０．９０ ３８．６０ ２１．７４ ２２．７２ ２９．３５ ３１．９７ ＮＮ ＮＮ
Ｄ ９１．８３ ７９．０９ １７１．７０１３７．４０ ７０．２０ ６９．７０ １５．８０ ８．９４ １７．２０ １１．３０ ｎ ｎ
Ｅ ９０．９２ ６９．７７ １８８．１０１０６．３０ ６０．４０ ４６．２０ ２３．６３ １３．４３ ２５．９９ １９．２５ Ｎ Ｎ
速效Ｐ Ａ ０．５４ ０．３２ １．４３ １．９８ ０．１０ ０．０２ ０．３７ ０．３４ ６７．９８ １０５．８０ ｎ ｎ
ＡｖａｉｌａｂｌｅＰ Ｂ ０．８０ ０．４８ ２．２４ １．４６ ０．０７ ０．０４ ０．３７ ０．２４ ４６．８８ ５１．１３ ｎ ｎ
（ｍｇ／ｋｇ） Ｃ １．７３ １．３６ ５．８８ ３．３０ ０．６５ ０．５３ ０．７９ ０．６２ ４５．７１ ４６．０１ Ｎ Ｎ
Ｄ １．１５ ０．８８ ２．４８ １．８０ ０．７０ ０．６０ ０．２７ ０．２２ ２３．２４ ２５．０８ Ｎ ｎ
Ｅ ２．５３ １．９７ ５．５８ ６．４０ １．２３ ０．８２ ０．７３ ０．７８ ２８．８４ ３９．３０ Ｎ ｎ
ＮＯ３－Ｎ Ａ ２．３１ １．７６ ２８．０３ １１．７３ ０．３１ ０．６７ ３．２０ １．６８ １３８．４０ ９５．５４ ｎ ｎ
（ｍｇ／ｋｇ） Ｂ １．５４ １．３５ ４．３３ ４．００ ０．９０ ０．６６ ０．４７ ０．５１ ３０．４６ ３７．６０ ｎ ｎ
Ｃ １．７２ １．５５ ２．６０ ２．７７ ０．５１ ０．４１ ０．４５ ０．５８ ２６．２０ ３７．５５ Ｎ Ｎ
Ｄ ２．７２ １．９３ ６．９９ ３．１３ １．１４ １．１４ １．１０ ０．４６ ４０．４６ ２３．６４ Ｎ Ｎ
Ｅ ２．１２ ２．５２ ６．２２ ９．８８ ０．５８ ０．６０ ０．７９ １．６１ ３７．４２ ６３．７７ Ｎ ｎ
ＮＨ４＋Ｎ Ａ ２．４７ ２．３０ ３．６０ ８．９５ １．８１ １．６７ ０．４０ ０．８９ １６．３４ ３８．８１ Ｎ ＮＮ
（ｍｇ／ｋｇ） Ｂ ３．２７ ３．０３ ５．２９ ６．５７ ２．３０ １．４９ ０．６７ ０．７５ ２０．６２ ２４．６５ Ｎ ＮＮ
Ｃ ４．２８ ４．５１ １０．６０ ９．１９ １．６２ ０．８７ １．９１ ２．１４ ４４．５７ ４７．４０ Ｎ Ｎ
Ｄ ２．４０ ２．２４ ３．６７ ３．７３ １．７０ １．７２ ０．４０ ０．３５ １６．５０ １５．５７ Ｎ Ｎ
Ｅ ４．３４ ３．４０ ６．６７ ５．３０ ０．８９ １．２１ ０．９５ ０．７７ ２２．００ ２２．６０ Ｎ Ｎ
ｎ表示服从自然对数正态分布；Ｎ表示服从正态分布；ＮＮ表示服从近似正态分布。
ｎ：Ｌｏｇｎｏｒｍａｌｄｉｓｔｒｉｂｕｔｉｏｎ；Ｎ：Ｎｏｒｍａｌｄｉｓｔｒｉｂｕｔｉｏｎ；ＮＮ：Ｎｅａｒｎｏｒｍａｌｄｉｓｔｒｉｂｕｔｉｏｎ．
２．２ 黄土丘陵区撂荒群落土壤速效养分空间变异的地统计学分析
图２和图３是研究区不同弃耕年限样地土壤速效养分的半方差函数图及理论曲线，根据拟合出的半方差函
数模型得到半方差函数重要参数（表３），由半方差函数参数表可知，不同弃耕年限样地的０～２０和２０～４０ｃｍ土
壤速效养分半方差函数具有不同的变异函数模型，决定系数较大，说明实验半方差与理论模型相关性较好。结构
方差与基台值之比表明，０～２０ｃｍ速效钾在撂荒６～７年样地具有中等空间相关性，其结构方差与基台值之比为
７０％＜７５％，而其他年限撂荒样地在０～２０和２０～４０ｃｍ的比值为７９．８％～９９．９％，＞７５％，具有强烈的空间相
关性，０～２０ｃｍ速效钾在研究区５块样地的空间变异大小为撂荒３年＝撂荒３～４年＞撂荒２５年＞撂荒１４～１５
年＞撂荒６～７年样地，而２０～４０ｃｍ速效钾在所有弃耕年限样地的比值为７９．８％～９３．１％，＞７５％，具有强烈
的空间相关性，各样地的变异大小为撂荒１４～１５年＞撂荒６～７年＞撂荒３～４年＞撂荒２５年＞撂荒３年；０～
２０ｃｍ的速效磷在各撂荒年限样地的比值为７７．４％～９９．９％，表现为强烈的空间相关性，各样地变异大小为撂荒
１４～１５年＞撂荒３～４年＞撂荒２５年＞撂荒６～７年＞撂荒３年，２０～４０ｃｍ的比值为８２．１％～９９．９％，也具有
强烈的空间相关性，变异大小为撂荒３年＞撂荒１４～１５年＞撂荒６～７年＞撂荒２５年＞撂荒３～４年；０～２０ｃｍ
硝态氮在撂荒６～７年样地结构半方差与基台值之比为７４．８％，表现为中等空间相关性，在其余各撂荒样地比值
均接近１，说明硝态氮含量具有强烈的空间相关性，变异大小无显著差异，而２０～４０ｃｍ硝态氮在撂荒２５年样地
的比值为５０％，属于中等相关性，其他样地的比值均在８６．７％～９９．７％，表现为强烈的空间相关性，变异大小顺
９７第２０卷第２期 草业学报２０１１年
序为撂荒３年＞撂荒１４～１５年＞撂荒３～４年＞撂荒６～７年＞撂荒２５年；０～２０ｃｍ铵态氮在所研究各样地的
结构半方差与基台值之比为７９．５％～９９．９％，具有强烈的空间相关性，变异大小顺序为撂荒１４～１５年＞撂荒３
～４年＞撂荒６～７年＞撂荒２５年＞撂荒３年，而２０～４０ｃｍ铵态氮在各样地的比值为８０．８％～９９．９％，表现为
强烈的空间相关性，变异大小顺序为撂荒１４～１５年＞撂荒３～４年＞撂荒３年＞撂荒６～７年＞撂荒２５年。通
过对结构方差和基台值的比值分析可知，在研究区不同撂荒年限群落０～２０和２０～４０ｃｍ土壤速效养分的空间
变异主要属于结构性变异，由土壤母质、地形、气候等结构性因素引起；并且随着植被群落撂荒年限的不同，结构
性变异在总变异中所占的比值不同。
图２ 撂荒群落土壤速效钾、速效磷养分含量半方差函数曲线图
犉犻犵．２ 犛犲犿犻狏犪狉犻狅犵狉犪犿狅犳狊狅犻犾犪狏犪犻犾犪犫犾犲犓犪狀犱犪狏犪犻犾犪犫犾犲犘犻狀犪犫犪狀犱狅狀犲犱狅犾犱犳犻犲犾犱狊
图３ 撂荒群落土壤铵态氮、硝态氮养分含量半方差函数曲线图
犉犻犵．３ 犛犲犿犻狏犪狉犻狅犵狉犪犿狅犳狊狅犻犾狀犻狋狉犪狋犲狀犻狋狉狅犵犲狀犪狀犱犪犿犿狅狀犻狌犿狀犻狋狉狅犵犲狀犻狀犪犫犪狀犱狅狀犲犱狅犾犱犳犻犲犾犱狊
０８ ＡＣＴＡＰＲＡＴＡＣＵＬＴＵＲＡＥＳＩＮＩＣＡ（２０１１） Ｖｏｌ．２０，Ｎｏ．２
表３ 不同撂荒年限群落土壤速效养分含量地统计分析
犜犪犫犾犲３ 犌犲狅狊狋犪狋犻狊狋犻犮狊犪狀犪犾狔狊犻狊狉犲狊狌犾狋狊狅犳狊狅犻犾犪狏犪犻犾犪犫犾犲狀狌狋狉犻犲狀狋狊犻狀狅犾犱犳犻犲犾犱犮狅犿犿狌狀犻狋犻犲狊狑犻狋犺犱犻犳犳犲狉犲狀狋犪犫犪狀犱狅狀犿犲狀狋犪犵犲狊
分析
项目
Ｉｔｅｍｓ
样地
编号
Ｐｌｏｔ
Ｎｏ．
函数模型
Ｍｏｄｅｌｔｙｐｅ
０～２０
ｃｍ
２０～４０
ｃｍ
块金值
Ｎｕｇｇｅｔ（Ｃ０）
０～２０
ｃｍ
２０～４０
ｃｍ
基台值
Ｓｉｌ （Ｃ＋Ｃ０）
０～２０
ｃｍ
２０～４０
ｃｍ
Ｃ／（Ｃ＋Ｃ０）
０～２０
ｃｍ
２０～４０
ｃｍ
变程
Ｒａｎｇｅ（ａ）
０～２０
ｃｍ
２０～４０
ｃｍ
决定系数Ｃｏｅｆｆｅｃｉｅｎｔ
ｏｆｄｅｔｅｒｍｉｎａｔｉｏｎ（ｒ２）
０～２０
ｃｍ
２０～４０
ｃｍ
速效Ｋ Ａ 高斯Ｇａｕｓｓｉａｎ 球状Ｓｐｈｅｒｉｃａｌ ０．１０００ １８．７０００ ９６．２０００ ９２．７０００ ０．９９９００．７９８０１．３０００１４．８５００ ０．７０７０ ０．７６６０
ＡｖａｉｌａｂｌｅＫ Ｂ 高斯Ｇａｕｓｓｉａｎ 球状Ｓｐｈｅｒｉｃａｌ ０．１０００ ４．１０００ ８３．９０００ ３７．２０００ ０．９９９００．８９００２．２２００ ３．４６００ ０．７２００ ０．６３４０
Ｃ 高斯Ｇａｕｓｓｉａｎ 高斯Ｇａｕｓｓｉａｎ３０６．００００ ７４．００００１０２２．９０００１０５８．９０００ ０．７００００．９３００１３．３５００ ９．２１００ ０．７３８０ ０．９６７０
Ｄ 球状Ｓｐｈｅｒｉｃａｌ 球状Ｓｐｈｅｒｉｃａｌ ３８．００００ ９．９０００４８６．９０００ １４３．７０００ ０．９２２００．９３１０１４．５７００１３．７２００ ０．６４００ ０．６０００
Ｅ 指数Ｅｘｐｏｎｅｎｔｉａｌ高斯Ｇａｕｓｓｉａｎ ４０．００００ ４４．００００６３６．４０００ ３９７．４０００ ０．９３７００．８８９０５．５５００１０．６５００ ０．７６８０ ０．９４７０
速效Ｐ Ａ 高斯Ｇａｕｓｓｉａｎ 球状Ｓｐｈｅｒｉｃａｌ ０．０３００ ０．０００１ ０．１３１４ ０．１０９２ ０．７７４００．９９９０１．４３７６ １．５３００ ０．７２４０ ０．５６７０
ＡｖａｉｌａｂｌｅＰ Ｂ 球状Ｓｐｈｅｒｉｃａｌ 指数Ｅｘｐｏｎｅｎｔｉａｌ０．００５０ ０．０１３０ ０．１５６０ ０．０７３０ ０．９６８００．８２１０２．５６００ ８．３７００ ０．４６００ ０．８５９０
Ｃ 球状Ｓｐｈｅｒｉｃａｌ 指数Ｅｘｐｏｎｅｎｔｉａｌ０．２１６０ ０．０３００ １．００５０ ０．９２００ ０．７８５００．９７００１０．７４００２４．２１００ ０．７０３０ ０．９０００
Ｄ 高斯Ｇａｕｓｓｉａｎ 高斯Ｇａｕｓｓｉａｎ ０．０００１ ０．０００１ ０．０７００ ０．０４９０ ０．９９９００．９９８０１．８９００ １．６６２８ ０．６８８０ ０．７４４０
Ｅ 球状Ｓｐｈｅｒｉｃａｌ 球状Ｓｐｈｅｒｉｃａｌ ０．０９４０ ０．１５００ ０．８５００ １．２４３０ ０．８８９００．８８００１１．６１００１４．６６００ ０．８１７０ ０．７４３０
ＮＯ３－Ｎ Ａ 高斯Ｇａｕｓｓｉａｎ 指数Ｅｘｐｏｎｅｎｔｉａｌ０．０１００ ０．０１００ ９．７１００ ３．５５７０ ０．９９９００．９９７０２．００９２ ９．２７００ ０．３３００ ０．５１９０
Ｂ 指数Ｅｘｐｏｎｅｎｔｉａｌ球状Ｓｐｈｅｒｉｃａｌ ０．０００１ ０．０１３０ ０．２２６２ ０．３３００ １．０００００．９６１０５．０４００ ７．９０００ ０．６３４０ ０．８３８０
Ｃ 指数Ｅｘｐｏｎｅｎｔｉａｌ高斯Ｇａｕｓｓｉａｎ ０．０７１０ ０．０７９０ ０．２８２８ ０．５９６０ ０．７４８００．８６７０１３．５９００ ８．６４２０ ０．７０７０ ０．９４５０
Ｄ 指数Ｅｘｐｏｎｅｎｔｉａｌ球状Ｓｐｈｅｒｉｃａｌ ０．００１０ ０．００２７ １．１７３０ ０．２１９４ ０．９９９００．９８８０３．６６００ ３．２０００ ０．４９５０ ０．８５５０
Ｅ 高斯Ｇａｕｓｓｉａｎ 高斯Ｇａｕｓｓｉａｎ ０．００１０ １．３４９０ ０．６２００ ２．６９９０ ０．９９８００．５０００１．４９００ ３．４６４０ ０．４９９０ ０．５４８０
ＮＨ４＋Ｎ Ａ 球状Ｓｐｈｅｒｉｃａｌ 高斯Ｇａｕｓｓｉａｎ ０．０５３０ ０．０６５０ ０．２５７０ ２．１４００ ０．７９５００．９７００１１．９６００１１．２８００ ０．７１６０ ０．６８５０
Ｂ 高斯Ｇａｕｓｓｉａｎ 高斯Ｇａｕｓｓｉａｎ ０．００１０ ０．００１０ ０．４２８０ ０．６３９０ ０．９９８００．９９８０１．６９７０ ０．７１００ ０．７３１０ ０．３６６０
Ｃ 线性Ｌｉｎｅａｒ 线性Ｌｉｎｅａｒ ０．６７００ ０．６７００ ５．５３２０ ６．９０１０ ０．８８０００．９０４０６．７２６６ ６．７２６６ ０．８８７０ ０．９４２０
Ｄ 高斯Ｇａｕｓｓｉａｎ 高斯Ｇａｕｓｓｉａｎ ０．０００１ ０．０００１ ０．１５２０ ０．１２６０ ０．９９９００．９９９０１．４７２０ １．４０３０ ０．７４９０ ０．８９００
Ｅ 指数Ｅｘｐｏｎｅｎｔｉａｌ指数Ｅｘｐｏｎｅｎｔｉａｌ０．１５００ ０．１４００ ０．９３００ ０．７１７０ ０．８３７００．８０８０４．６０００ ８．５２００ ０．５２００ ０．６７６０
变程反应区域化变量空间相关范围的大小，与观测尺度以及在取样尺度上影响土壤养分的各种生态过程相
互作用有关。不同撂荒年限样地速效钾在０～２０和２０～４０ｃｍ土层的变程分别为１．３０～１４．５７和３．４６～１４．８５
ｍ；速效磷在０～２０和２０～４０ｃｍ土层的变程分别为１．４４～１１．６１和１．５３～２４．２１ｍ；硝态氮在０～２０和２０～４０
ｃｍ的变程分别为１．４９～１３．５９和３．２０～９．２７ｍ；铵态氮在０～２０和２０～４０ｃｍ土层的变程为１．４７～１１．９６和
０．７１～１１．２８ｍ（表３）。不同撂荒年限样地土壤速效养分变程有较大差别，说明植被演替恢复对速效养分的空间
分布变异性具有一定影响作用，４种速效养分在相同样地相同深度土层的变程也不同，说明４种速效养分的空间
自相关尺度不同，它们的生态过程在不同尺度上起作用。
３ 结论与讨论
黄土丘陵区沟壑纵横，坡耕地居多，水土流失极为严重。此外，由于研究区样地在撂荒前期耕作时间较长，广
种薄收，施肥不足，撂荒地养分一般较为贫瘠。依据黄土高原土壤养分含量分级［２９］，研究区不同年限撂荒样地
０～２０ｃｍ速效钾含量属于中等水平，２０～４０ｃｍ属于低水平，速效磷和速效氮含量均属于极低水平，说明该研究
区弃耕样地土壤肥力处于低的水平，这与刘梦云等［３０］的研究成果相一致。
在研究区植被群落撂荒演替过程中，０～２０和２０～４０ｃｍ土壤速效养分在演替中后期含量较高。王俊明和
张兴昌［３１］的研究也表明在退耕草地演替前期土壤肥力并不会有明显好转。原因主要是在撂荒演替前期植被覆
盖率低，流失和淋洗较为严重，但随着植被群落演替的进行，群落生物量不断增加，有机质积累越来越多，在矿化
作用和表层富积共同作用下，土壤养分含量又逐渐升高。同时，０～２０ｃｍ土壤速效养分显著高于２０～４０ｃｍ养
１８第２０卷第２期 草业学报２０１１年
分含量，与董莉丽和郑粉莉［３２］的研究结果一致，可能与撂荒前施肥及表层富积作用有关［３３，３４］。
对不同撂荒年限样地土壤速效养分空间分布进行异质性分析，实验半方差函数与理论模型相关较显著。撂
荒６～７年样地０～２０ｃｍ土壤速效钾和硝态氮以及撂荒２５年样地２０～４０ｃｍ硝态氮具有中等空间相关性，其他
撂荒群落０～２０和２０～４０ｃｍ土壤速效养分均具有强烈的空间相关性。
随植被群落撂荒演替年限的增加，研究区样地０～２０和２０～４０ｃｍ土壤速效养分在群落演替不同时期空间
异质性分布具有明显差异。０～２０ｃｍ土壤速效钾空间变异性在演替前期最大，演替中期和后期较小，２０～４０ｃｍ
异质性则在演替中期大于演替前期和后期；０～２０ｃｍ速效磷空间变异在演替中期最大，２０～４０ｃｍ变异性在演替
前期和中期比演替后期大；０～２０ｃｍ硝态氮空间变异性在演替前期最大，２０～４０ｃｍ变异性在演替前期和中期最
大；０～２０和２０～４０ｃｍ铵态氮空间分布异质性均在演替中期最大。由于土壤养分的空间分布异质性既受到人
为因素的干扰（施肥，耕作方式），又受到非人为因素影响（雨水冲刷和植被的土壤养分效应）；０～２０ｃｍ土壤速效
钾和硝态氮分布在群落恢复演替前期异质性大，可能是由于农田在耕种过程中施肥不均匀等人为因素引起的，但
随着植被群落的演替恢复，植被的空间分布趋向于匀质化，造成多数种群斑块消融，从而土壤养分空间分布也趋
向于匀质化，使得演替中后期空间变异性减小，但２０～４０ｃｍ速效钾在演替中期异质性较高，主要是演替中期的
建群种为长芒草，在２０～４０ｃｍ土层，根系对养分的聚集吸收使得土壤养分空间异质性最大，这与杜峰等［３５］的研
究结果一致。速效磷和铵态氮在０～２０和２０～４０ｃｍ土壤空间变异性均在群落演替中期较大，原因可能是随着
弃耕年限的延长，植被群落逐渐得到恢复，植被生物量不断增加并在表层的富积，经过风力和降水的冲刷作用，使
得枯枝落叶聚集在植株根部，矿化后养分含量升高，同时随着植物根系的生长，对周围养分的聚集能力越来越强，
特别是本研究区群落演替中期的主要植被为长芒草（多年生密丛禾草）和铁杆蒿（草本状半灌木），其根系对养分
有较强的吸收和聚集能力，可以形成养分斑块即“沃岛”效应［３６］，造成土壤养分在空间上的分布和利用不均匀。
植被恢复演替对土壤养分空间异质性分布具有影响作用，而土壤养分资源空间分布异质性的变化对黄土丘
陵区植被恢复建设具有推动作用，因而土壤特性空间变异研究可为土壤质量的恢复与改善，生态环境的修复与重
建提供理论依据。在我国水土流失最为严重的黄土丘陵地区，植被恢复是解决此类问题的重要途径之一，土壤空
间变异的研究可以为选择合理的植被类型和结构提供参考。同时，本研究中速效养分空间变异变程在撂荒３年、
６～７年、１４～１５年和２５年样地均有超过最大采样距离的情况，说明在一些撂荒年限样地研究尺度较小，在以后
的研究中应该适当放大采样间距。
参考文献：
［１］ ＨｕｇｇｅｔｔＲＪ．Ｓｏｉｌｃｈｒｏｎｏｓｅｑｕｅｎｃｅｓ，ｓｏｉｌｄｅｖｅｌｏｐｍｅｎｔ，ａｎｄｓｏｉｌｅｖｏｌｕｔｉｏｎ：Ａｃｒｉｔｉｃａｌｒｅｖｉｅｗ［Ｊ］．Ｃａｔｅｎａ，１９９８，３２：１５５１７２．
［２］ ＴｒａｎｇｍａｒＢＢ，ＹｏｓｔＲＳ，ＵｅｈａｒａＧ．Ａｐｐｌｉｃａｔｉｏｎｏｆｇｅｏｓｔａｔｉｓｔｉｃｓｔｏｓｐａｔｉａｌｓｔｕｄｉｅｓｏｆｓｏｉｌｐｒｏｐｅｒｔｉｅｓ［Ｊ］．ＡｄｖａｎｃｅｓｉｎＡｇｒｏｎｏ
ｍｙ，１９８５，３８：４５９４．
［３］ 陈效民，吴华山，孙静红．太湖地区农田土壤中铵态氮和硝态氮的时空变异［Ｊ］．环境科学，２００６，２７（６）：１２１７１２２２．
［４］ 张春霞，郝明德，王旭刚，等．黄土高原沟壑区小流域土壤养分分布特征［Ｊ］．水土保持研究，２００３，１０（１）：７８８０．
［５］ 孙洪全．地统计学及其应用［Ｍ］．北京：中国矿业大学出版社，１９９０．
［６］ ＬｉＪＷ，ＲｉｃｈｔｅｒＤａｎｉｅｌｄｅＢ，ＭｅｎｄｏｚａＡ，犲狋犪犾．Ｅｆｆｅｃｔｓｏｆｌａｎｄｕｓｅｈｉｓｔｏｒｙｏｎｓｏｉｌｓｐａｔｉａｌｈｅｔｅｒｏｇｅｎｅｉｔｙｏｆｍａｃｒｏａｎｄｔｒａｃｅｅｌｅ
ｍｅｎｔｓｉｎｔｈｅＳｏｕｔｈｅｒｎＰｉｅｄｍｏｎｔＵＳＡ［Ｊ］．Ｇｅｏｄｅｒｍａ，２０１０，１５６：６０７３．
［７］ ＯｒｗｉｎＫＨ，ＢｅｒｔｒａｎＪＥ，ＣｌｏｕｇｈＴＪ，犲狋犪犾．Ｓｈｏｒｔｔｅｒｍｃｏｎｓｅｑｕｅｎｃｅｓｏｆｓｐａｔｉａｌｈｅｔｅｒｏｇｅｎｅｉｔｙｉｎｓｏｉｌｎｉｔｒｏｇｅｎｃｏｎｃｅｎｔｒａｔｉｏｎｓ
ｃａｕｓｅｄｂｙｕｒｉｎｅｐａｔｃｈｅｓｏｆｄｉｆｆｅｒｅｎｔｓｉｚｅｓ［Ｊ］．ＡｐｐｌｉｅｄＳｏｉｌＥｃｏｌｏｇｙ，２００９，４２：２７１２７８．
［８］ ＬｏｅｃｋｅＴＤ，ＲｏｂｅｒｔｓｏｎＧＰ．Ｓｏｉｌｒｅｓｏｕｒｃｅｈｅｔｅｒｏｇｅｎｅｉｔｙｉｎｔｅｒｍｓｏｆｌｉｔｔｅｒａｇｇｒｅｇａｔｉｏｎｐｒｏｍｏｔｅｓｎｉｔｒｏｕｓｏｘｉｄｅｆｌｕｘｅｓａｎｄｓｌｏｗｓ
ｄｅｃｏｍｐｏｓｉｔｉｏｎ［Ｊ］．ＳｏｉｌＢｉｏｌｏｇｙ＆Ｂｉｌｃｈｅｍｉｓｔｒｙ，２００９，４１：２２８２３５．
［９］ 庞夙，李廷轩，王永东，等．土壤速效氮、磷、钾含量空间变异特征及其影响因子［Ｊ］．植物营养与肥料学报，２００９，１５（１）：
１１４１２０．
［１０］ 邢永强，李金荣，李金玲，等．土壤中硝态氮的空间变异研究［Ｊ］．灌溉排水学报，２００８，２７（３）：１０６１０８．
［１１］ 徐海燕，聂宜民，赵文武，等．寿光市土壤养分特性及空间变异性分析［Ｊ］．农业现代化研究，２００９，３０（１）：１１８１２１．
２８ ＡＣＴＡＰＲＡＴＡＣＵＬＴＵＲＡＥＳＩＮＩＣＡ（２０１１） Ｖｏｌ．２０，Ｎｏ．２
［１２］ ＣａｍｂａｒｄｅｌａＣＡ，ＹｏｓｔＲＳ，ＷａｄｅＭＫ，犲狋犪犾．ＦｉｅｌｄｓｃａｌｅｖａｒｉａｂｉｌｉｔｙｏｆｓｏｉｌｐｒｏｐｅｒｔｉｅｓｉｎｃｅｎｔｒａｌＩｏｗａｓｏｉｌｓ［Ｊ］．ＳｏｉｌＳｃｉｅｎｃｅ
ＳｏｃｉｅｔｙｏｆＡｍｅｒｉｃａＪｏｕｒｎａｌ，１９９４，５８：１５０１１５１１．
［１３］ ＭｅｇｒａｗＴ．ＦｅｒｔｉｌｉｔｙｖａｒｉａｂｉｌｉｔｙｉｎｔｈｅＭｉｎｎｅｓｏｔａｒｉｖｅｒｖａｌｅｙｗａｔｅｒｓｈｅｄｉｎ１９９３，ａｓｄｅｔｅｒｍｉｎｅｄｆｒｏｍｇｒｉｄｔｅｓｔｉｎｇｒｅｓｕｌｔｏｎ
５２０００ａｃｒｅｓｉｎｃｏｍｍｅｒｃｉａｌｆｉｅｌｄｓ［Ｊ］．ＳｉｔｅＳｐｅｃｉｆｉｃＭａｎａｇｅｍｅｎｔｆｏｒＡｇｒｉｃｕｌｔｕｒａｌＳｙｓｔｅｒｍｓＡＳＡＣＳＳＡＳＳＳＡ，１９９５，（２）：１６７
１７４．
［１４］ 李广，侯扶江．临泽盐化草甸土壤理化性质的空间格局分析［Ｊ］．草业学报，２００３，１２（２）：４６５３．
［１５］ 李香珍，张淑敏，邢学荣．小叶锦鸡儿灌丛引起的植物生物量和土壤化学元素含量的空间变异［Ｊ］．草业学报，２００２，
１１（１）：１３．
［１６］ 张有山，林启美，秦耀东，等．大比例尺区域土壤养分空间变异定量分析［Ｊ］．华北农学报，１９９８，１３（１）：１２２１２８．
［１７］ 刘冬碧，余常兵，熊桂云，等．大比例尺度下土壤的养分特征及其空间变异性研究［Ｊ］．华中农业大学学报，２００４，２３（５）：
５２４５２７．
［１８］ 杨玉玲，田长彦，盛建东，等．灌淤土壤有机质、全量氮磷钾空间变异性初探［Ｊ］．干旱地区农业研究，２００２，２０（３）：２６３０．
［１９］ 朱益玲，刘洪斌，江希流．江津市紫色土中Ｎ、Ｐ养分元素区域空间变异性研究［Ｊ］．环境科学，２００４，２５（１）：１３８１４３．
［２０］ 杜峰，梁宗锁，徐学选，等．陕北黄土丘陵区撂荒草地群落生物量及植被土壤养分效应［Ｊ］．生态学报，２００７，２７（５）：１６７３
１６８２．
［２１］ 初玉，杨慧玲，朱选伟，等．浑善达克沙地小叶锦鸡儿灌丛的空间异质性［Ｊ］．生态学报，２００５，２５（１２）：３２９４３３００．
［２２］ 张璞进，杨稢，宋炳煜，等．藏锦鸡儿群落土壤资源空间异质性［Ｊ］．植物生态学报，２００９，３３（２）：３３８３４６．
［２３］ ＴｓｅｇａｙｅＴ，ＨｉｌＲＬ．Ｉｎｔｅｎｓｉｖｅｔｉｌａｇｅｅｆｆｅｃｔｓｏｎｓｐａｔｉａｌｖａｒｉａｂｉｌｉｔｙｏｆｓｏｉｌｔｅｓｔ，ｐｌａｎｔｇｒｏｗｔｈ，ａｎｄｎｕｔｒｉｅｎｔｕｐｔａｋｅｍｅａｓｕｒｅ
ｍｅｎｔ［Ｊ］．ＳｏｉｌＳｃｉｅｎｃｅ，１９９８，１６３（２）：１５５１６５．
［２４］ 焦峰，温仲明，焦菊英，等．黄土丘陵区退耕地土壤养分变异特征［Ｊ］．植物营养与肥料学报，２００５，１１（６）：７２４７３０．
［２５］ 陈伏生，曾德慧．科尔沁沙地退化草场土壤养分的空间结构分析［Ｊ］．草业学报，２００４，１３（１）：３９４４．
［２６］ 山仑，徐炳成．黄土高原半干旱地区建设稳定人工草地的探讨［Ｊ］．草业学报，２００９，１８（２）：１２．
［２７］ 鲍士旦．土壤农化分析［Ｍ］．北京：农业出版社，１９８６：３３８９．
［２８］ ＮｉｅｌｓｅｎＤＲ，ＢｏｕｍａＪ．ＳｏｉｌＳｐａｔｉａｌＶａｒｉａｂｉｌｉｔｙ［Ｍ］．ＰｕｄｏｃＷａｇｅｎｉｎｇｅｎ，１９８５：２３０．
［２９］ 贾恒义，彭琳，彭祥林，等．黄土高原地区土壤养分资源分区及其评价［Ｊ］．水土保持学报，１９９４，８（３）：２２２８．
［３０］ 刘梦云，寇宝平，常庆瑞，等．安塞小流域土壤养分分布特征研究［Ｊ］．西北农林科技大学学报（自然科学版），２００２，３０（６）：
２１２４．
［３１］ 王俊明，张兴昌．退耕草地演替过程中的碳储量变化［Ｊ］．草业学报，２００９，１８（１）：１８．
［３２］ 董莉丽，郑粉莉．黄土丘陵区不同土地利用类型下土壤酶活性和养分特征［Ｊ］．生态环境，２００８，１７（５）：２０５０２０５８．
［３３］ 董文斌，马玉寿，董全民，等．退耕还（林）草多年生栽培草地土壤养分特征研究［Ｊ］．草业科学，２０１０，２７（４）：４６５０．
［３４］ 李晓东，魏龙，张永超，等．土地利用方式对陇中黄土高原土壤理化性状的影响［Ｊ］．草业学报，２００９，１８（４）：１０３１１０．
［３５］ 杜峰，梁宗锁，徐学选，等．陕北黄土丘陵区撂荒群落土壤养分与地上生物量空间异质性［Ｊ］．生态学报，２００８，２８（１）：１３
２２．
［３６］ 秦伟，朱清科，刘中奇，等．黄土丘陵沟壑区退耕地植被自然演替系列及其植物物种多样性特征［Ｊ］．干旱区研究，２００８，
２５（４）：５０７５１３．
３８第２０卷第２期 草业学报２０１１年
犛狆犪狋犻犪犾犺犲狋犲狉狅犵犲狀犲犻狋狔狅犳犪狏犪犻犾犪犫犾犲狊狅犻犾狀狌狋狉犻犲狀狋狊犻狀犪犫犪狀犱狅狀犲犱狅犾犲犳犻犲犾犱
犮狅犿犿狌狀犻狋犻犲狊犻狀狋犺犲犔狅犲狊狊犎犻犾狔犚犲犵犻狅狀
ＳＵＩＹｕａｎｙｕａｎ１，ＤＵＦｅｎｇ１，２，ＺＨＡＮＧＸｉｎｇｃｈａｎｇ１，２
（１．ＣｏｌｅｇｅｏｆＲｅｓｏｕｒｃｅｓａｎｄＥｎｖｉｒｏｎｍｅｎｔ，ＮｏｒｔｈｗｅｓｔＵｎｉｖｅｒｓｉｔｙｏｆＡｇｒｉｃｕｌｔｕｒｅａｎｄＦｏｒｅｓｔｒｙ，Ｙａｎｇｌｉｎｇ
７１２１００，Ｃｈｉｎａ；２．ＩｎｓｔｉｔｕｔｅｏｆＳｏｉｌａｎｄＷａｔｅｒＣｏｎｓｅｒｖａｔｉｏｎ，Ｙａｎｇｌｉｎｇ７１２１００，Ｃｈｉｎａ）
犃犫狊狋狉犪犮狋：Ｕｓｉｎｇｃｌａｓｓｉｃａｌｓｔａｔｉｓｔｉｃｓａｎｄｇｅｏｓｔａｔｉｓｔｉｃｓ，ｔｈｅｓｐａｔｉａｌｈｅｔｅｒｏｇｅｎｅｉｔｙｏｆｓｏｉｌａｖａｉｌａｂｌｅｎｕｔｒｉｅｎｔｓｉｎｆｉｖｅ
ｔｙｐｉｃａｌｏｌｄｆｉｅｌｄｃｏｍｍｕｎｉｔｉｅｓｗｉｔｈｄｉｆｆｅｒｅｎｔａｂａｎｄｏｎｍｅｎｔａｇｅｓｉｎＬｏｅｓｓＨｉｌｙＲｅｇｉｏｎｗｅｒｅｄｅｍｏｎｓｔｒａｔｅｄａｔｃｏｍ
ｍｕｎｉｔｙｓｃａｌｅ．Ｔｈｅｒｅｓｕｌｔｓｉｎｄｉｃａｔｅｄｔｈａｔ：１）Ｉｎｔｈｅｃａｓｅｏｆｎｕｔｒｉｅｎｔｌｅｖｅｌｓ，ｏｌｄｆｉｅｌｄｓｉｎｔｈｅｓｔｕｄｙａｒｅａａｒｅａｔｔｒｉｂ
ｕｔｅｄｂｙｒｅｌａｔｉｖｅｐｏｏｒｓｏｉｌｓ，０－２０ａｎｄ２０－４０ｃｍａｖａｉｌａｂｌｅｎｕｔｒｉｅｎｔｓｉｎｃｒｅａｓｅｄｉｎｍｉｄｄｌｅａｂａｎｄｏｎｍｅｎｔｓｔａｇｅｏｆ
ｓｕｃｃｅｓｓｉｏｎ，ｎｕｔｒｉｅｎｔｃｏｎｔｅｎｔｉｎｓｕｒｆａｃｅｓｏｉｌｗａｓｈｉｇｈｅｒｔｈａｎｔｈａｔｉｎｓｕｂｌａｙｅｒ．２）Ｔｈｅｖａｒｉａｂｉｌｉｔｙｏｆａｖａｉｌａｂｌｅｎｕ
ｔｒｉｅｎｔｓｗａｓｍａｉｎｌｙｃｏｍｐｏｓｅｄｏｆｓｔｒｕｃｔｕｒａｌｖａｒｉａｔｉｏｎ．Ｓｏｉｌａｖａｉｌａｂｌｅｐｏｔａｓｓｉｕｍａｎｄｎｉｔｒａｔｅｎｉｔｒｏｇｅｎｉｎ０－２０ｃｍ
ｏｆｔｈｅ６－７ｙｅａｒｓｏｌｄｆｉｅｌｄ，ａｎｄｓｕｂｓｕｒｆａｃｅｓｏｉｌｎｉｔｒａｔｅｎｉｔｒｏｇｅｎｉｎ２０－４０ｃｍｏｆｔｈｅ２５ｙｅａｒｓｏｌｄｆｉｅｌｄｈａｄｉｎｔｅｒ
ｍｅｄｉａｔｅｓｐａｔｉａｌｃｏｒｒｅｌａｔｉｏｎ．Ｉｎｏｔｈｅｒｏｌｄｆｉｅｌｄｓ，ａｖａｉｌａｂｌｅｎｕｔｒｉｅｎｔｓｉｎｓｕｒｆａｃｅａｎｄｓｕｂｓｕｒｆａｃｅｓｏｉｌｈａｄｓｔｒｏｎｇｓｐａ
ｔｉａｌｃｏｒｒｅｌａｔｉｏｎ．３）０－２０ｃｍａｖａｉｌａｂｌｅｐｏｔａｓｓｉｕｍａｎｄｎｉｔｒａｔｅｎｉｔｒｏｇｅｎｓｈｏｗｅｄｔｈｅｓｔｒｏｎｇｅｓｔｓｐａｔｉａｌｖａｒｉａｂｉｌｉｔｙｉｎ
ｅａｒｌｉｅｒａｂａｎｄｏｎｍｅｎｔｓｔａｇｅｏｆｓｕｃｃｅｓｓｉｏｎ；２０－４０ｃｍａｖａｉｌａｂｌｅｐｏｔａｓｓｉｕｍａｎｄｎｉｔｒａｔｅｎｉｔｒｏｇｅｎｓｈｏｗｅｄｔｈｅｓｔｒｏｎ
ｇｅｓｔｓｐａｔｉａｌｖａｒｉａｂｉｌｉｔｙｉｎｍｉｄｄｌｅａｎｄｅａｒｌｙｍｉｄｄｌｅｓｔａｇｅｏｆｓｕｃｃｅｓｓｉｏｎｒｅｓｐｅｃｔｉｖｅｌｙ；Ｓｐａｔｉａｌｖａｒｉａｂｉｌｉｔｙｏｆａｖａｉｌａｂｌｅ
ｐｈｏｓｐｈｏｒｕｓａｎｄａｍｍｏｎｉｕｍｎｉｔｒｏｇｅｎｉｎｓｕｒｆａｃｅｌａｙｅｒｗａｓｓｔｒｏｎｇｅｓｔｉｎｍｉｄｄｌｅａｂａｎｄｏｎｍｅｎｔｓｔａｇｅｏｆｓｕｃｃｅｓｓｉｏｎ；
Ｓｐａｔｉａｌｖａｒｉａｂｉｌｉｔｙｏｆａｖａｉｌａｂｌｅｐｈｏｓｐｈｏｒｕｓａｎｄａｍｍｏｎｉｕｍｎｉｔｒｏｇｅｎｉｎｓｕｂｓｕｒｆａｃｅｌａｙｅｒｗａｓｓｔｒｏｎｇｅｓｔｉｎｅａｒｌｙ
ｍｉｄｄｌｅａｎｄｍｉｄｄｌｅａｂａｎｄｏｎｍｅｎｔｓｔａｇｅｏｆｓｕｃｃｅｓｓｉｏｎｒｅｓｐｅｃｔｉｖｅｌｙ．Ｒｅｇａｒｄｉｎｇａｂｏｕｔｔｈｅｉｎｄｅｘｏｆｓｐａｔｉａｌｈｅｔｅｒｏｇｅ
ｎｅｉｔｙ，ｓｉｇｎｉｆｉｃａｎｔｄｉｆｆｅｒｅｎｃｅｓｗｅｒｅｆｏｕｎｄａｍｏｎｇｔｈｅ５ｏｌｄｆｉｅｌｄｃｏｍｍｕｎｉｔｉｅｓ，ｗｈｉｃｈｉｍｐｌｉｅｓｔｈａｔｔｈｅｐｒｉｎｃｉｐａｌｅｃ
ｏｌｏｇｉｃａｌｐｒｏｃｅｓｓｅｓｏｆｔｈｅｆｏｕｒａｖａｉｌａｂｌｅｓｏｉｌｎｕｔｒｉｅｎｔｓａｒｅａｓｓｏｒｔｅｄ．
犓犲狔狑狅狉犱狊：ＬｏｅｓｓＨｉｌｙＲｅｇｉｏｎ；ａｂａｎｄｏｎｅｄｏｌｅｆｉｅｌｄｃｏｍｍｕｎｉｔｉｅｓ；ａｖａｉｌａｂｌｅｓｏｉｌｎｕｔｒｉｅｎｔｓ；ｓｐａｔｉａｌｈｅｔｅｒｏｇｅｎｅｉｔｙ
４８ ＡＣＴＡＰＲＡＴＡＣＵＬＴＵＲＡＥＳＩＮＩＣＡ（２０１１） Ｖｏｌ．２０，Ｎｏ．２