河北平山片麻岩山区微地形表层土壤养分异质性分析

袁 振, 魏松坡, 贾黎明, 张亚雄, 刘龙龙

(北京林业大学 省部共建森林培育与保护教育部重点实验室, 北京100083)

河北平山片麻岩山区微地形表层土壤养分异质性分析

袁 振, 魏松坡, 贾黎明, 张亚雄, 刘龙龙

(北京林业大学 省部共建森林培育与保护教育部重点实验室, 北京100083)

以河北平山片麻岩山区7种微地形(坡顶、塌陷、巨石背阴、缓台、陡坎、谷坡、U形沟)及原状坡的表层土壤养分测定数据为基础，对比分析不同微地形及原状坡之间表层土壤养分的异质性，为片麻岩山区植被恢复提供科学依据。结果表明：(1) 微地形、坡向、风化层厚度均影响了表层土壤养分的变化，而微地形是引起这种变化的主要因子。(2) 微地形的表层土壤养分特征多数优于原状坡。其中，U形沟、塌陷和巨石背阴处有机质含量最大(68.94 g/kg，60.71 g/kg，62.77 g/kg)，显著高于其他微地形和原状坡，坡顶处最小(22.58 g/kg)。(3) 全氮与全磷含量规律基本一致，都在U形沟、塌陷和巨石背阴处最大，缓台和坡顶处较小。全钾含量在坡顶、陡坎、谷坡处最大(11.88 g/kg，12.56 g/kg，11.67 g/kg)，塌陷、巨石背阴处较小(9.39 g/kg，9.69 g/kg)。(4) 碱解氮含量在U形沟、塌陷和巨石背阴处最大(6.87 mg/kg，6.54 mg/kg，5.93 mg/kg)，显著高于其他微地形及原状坡，坡顶处最小(3.75 mg/kg)。U形沟处速效磷含量最大(3.55 mg/kg)，塌陷、巨石背阴、陡坎和谷坡次之，这5种微地形显著高于坡顶和原状坡，坡顶处最小(1.27 mg/kg)。速效钾含量在坡顶、陡坎达到最大(60.09 mg/kg，58.27 mg/kg)，这2种微地形显著高于其他微地形和原状坡。(5) 通过对7种微地形和原状坡的3个坡向的表层土壤养分进行灰色关联分析，得出U形沟、塌陷和巨石背阴的土壤养分条件最好，坡顶和原状坡的土壤养分条件最差。

片麻岩山区; 微地形; 土壤养分; 异质性; 灰色关联分析

土壤是地球上植物生长的基础，其理化性质受成土母质、地形条件、气候条件、人为活动等因素的综合作用[1]，在大尺度上，这些因素促使土壤理化性质具有较强的空间异质性。而在小尺度上，土壤的理化性质则因微地形的变化表现出一定的空间规律性[2-5]。微地形通过地表的起伏变化，导致局部的光照、热量、土壤水分、养分等资源进行再分配[6-8]。因此，微地形的变化势必对土壤的理化性质产生影响，而土壤的理化性质的差异又会导致不同微地形及原状坡的植物群落特征和功能产生异质性[9-11]。因此，研究不同微地形下土壤的空间异质性，对于指导困难立地植被恢复具有重要意义。

已有学者通过对微地形与土壤理化性质之间的关系进行相关研究，得出微地形对土壤的理化性质具有一定的影响[12-13]。高凯等[14]系统分析了浑善达克沙地不同微地形下土壤营养物质含量的变化规律，发现微地形对土壤营养物质含量有显著影响。张宏芝等[15]对黄土区自然恢复状态下的微地形土壤化学性质进行研究发现，不同微地形间的土壤化学性质有显著差异，并通过灰色关联分析得出，塌陷和缓台的土壤化学性质最优。赵秀芳等[16]分析了微地形坡面(坡顶、坡腰、坡底)对滨海吹填土土壤养分的影响，发现3种微地形间土壤养分含量差异显著。刘孝阳等[17]对平朔矿区微地形养分的研究也表明微地形的变化对土壤养分有很大的影响。罗光杰等[18]，张忠华等[19]研究发现，喀斯特山区地面岩石裸露所形成的石面、石沟、石缝等微地形均对土壤养分的空间分布产生重要影响。

我国片麻岩山地分布广泛，形成的风化物或土壤保水保肥能力差，是开展植被恢复较为困难的立地之一[20]。平山县是太行山地区典型的片麻岩山区，由于岩石风化、坡面径流，在片麻岩山区内形成许多形形色色的微地形。微地形一般是指小尺度的地形变化，Kikuchi[21]，Nagamatsu等[22]将丘陵地区微地形分为顶坡、上部边坡、谷头凹地、下部边坡、麓坡、泛滥性阶地和谷床等7类。朱清科等[23]将黄土高原地区坡面内地表的起伏形成的微地形分为浅沟、切沟、塌陷、缓台、陡坎等5类。本研究根据片麻岩山区的地形特征，通过广泛踏查，将片麻岩山区微地形分为坡顶、塌陷、巨石背阴、缓台、陡坎、谷坡和U形沟等7类。由于地形的差异，使光照、土壤水分、养分及风化层厚度等小生境在局部范围内产生微小的变化，从而导致地上植物群落特征和类型的不同，对生态系统结构的形成具有重要影响[24]。

目前，已有大量有关微地形的研究，但多集中于微地形对植物群落特征的影响[25-26]，对于微地形土壤理化性质的相关研究较少[15]，针对片麻岩山区微地形与土壤理化性质的研究尚未见报道。因此，本研究以河北平山自然恢复状态下的片麻岩区微地形为研究对象，系统分析片麻岩山区不同微地形条件下表层土壤养分含量的异质性，以期为区域片麻岩山区植被恢复提供科学依据。

1 材料与方法

1.1 研究区自然概况

平山县位于河北省西部太行山中段东麓。地理位置为东经113°31′—114°15′，北纬38°9′—38°47′，是太行山地区典型的片麻岩山区，全县地貌属于丘陵山地类型，海拔变化于111～2 281 m。该地区的母岩为变质深厚的片麻岩和页岩等，片麻岩山区的主要特点是：土壤干旱，土层浅薄、土壤沙性强，侵蚀严重。片麻岩是变质岩，因此具有比其他岩石容易风化的特点，一般易松动，易破碎，地表的松散固体物质丰富，未成土的半风化岩层厚度可达30～60 cm。山顶土层瘠薄，植被的退化往往产生水土流失，导致基岩裸露。年降水量约500 mm，集中于7—9月，一年中大部分时间多干旱少降雨。由于植被稀疏，土层较薄，降水多以地表径流的形式输出。现有的天然植被较稀疏，覆盖度较低且分布不均匀，植被主要以灌木和草本为主，沟谷地带有人工种植的乔木。灌木主要有：酸枣(Zizyphusjujuba)，荆条(Vitexnegundovar.heterophylla)，薄皮木(Leptodermisoblonga)等。草本植物主要有：铁杆蒿(Artemisiasacrorum)，白羊草(Bothriochloaischaemum)，达呼里胡枝子(Lespedezadavurica)，阿尔泰狗娃花(Heteropappusaltaicus)等。

1.2 片麻岩山区微地形类型划分

微地形一般指小尺度的地形变化。片麻岩表层风化程度深，未成土的半风化层厚度达40 cm左右，岩体有裂隙易松动、崩塌，容易在坡面形成大小不等、形状各异的局部小地形[20]。本文参考朱清科等[23]对黄土高原、杨永川等[24]对丘陵地区微地形的划分，通过广泛调查和分析，将河北平山片麻岩微地形划分为：坡顶、塌陷、巨石背阴、缓台、陡坎、谷坡、U形沟等7种(图1)，各自特征见表1。

图1 片麻岩山区微地形

1.3 研究方法

1.3.1 样地设置 于2015年7—9月，选取自然恢复的片麻岩区域，对该区域微地形的分布进行调查，选取阳坡、阴坡、半阴半阳坡，在每个坡向上根据微地形的实际大小进行样方设置，同时取与微地形相邻的原状坡作为对照样方。共计设置样方142个。

1.3.2 微立地因子测定 GPS测定样方经纬度、海拔，手持罗盘仪测坡向、坡度，同时记录坡位。土壤硬度采用中山式土壤硬度计测定，每个样方测量5次，求平均值，单位为mm。风化层厚度采用钢钎法，垂直样方插入，用钢卷尺测量深度，每个样方随机测量5次，求平均值[27]。每个样方取0—20 cm土层土壤1 kg，3次重复，风干后，过2 mm筛测定直径大于2 mm的石砾含量，计算表土层含石率。

1.3.3 土样采集及分析 采集土样时先去除地表的凋落物，在每个样方用土钻随机钻取5个0—20 cm表层土样，充分混匀后放入采集袋中带回实验室，风干研磨后，测定土壤养分。土壤有机质的测定采用重铬酸钾容量法—外加热法；全氮的测定采用凯氏定氮法；全磷的测定采用NaOH熔融，钼蓝比色法；全钾的测定采用NaOH熔融，火焰光度法；碱解氮的测定采用碱解扩散法；速效磷的测定采用0.5 mol/L的NaHCO3浸提，钼蓝比色法；速效钾的测定采用NH4OAc浸提，火焰光度法。

1.3.4 灰色关联分析法 对于两个系统之间共有的影响因素，其随时间或不同对象而变化的关联性大小的量度，称为关联度。在系统发展过程中，若两个因素变化的趋势具有一致性，即同步变化程度较高，即可谓二者关联程度较高；反之则较低。因此，灰色关联分析方法[28]，是根据因素之间发展趋势的相似或相异程度，作为衡量因素间关联程度的一种方法。在进行数据列关联分析时，必须先确定参考数列，然后比较其他数列同参考数列的接近程度，这样才能对其他数列进行比较，进而做出判断。

设x0={x0(k)k=1，2，…，n}为参考数列(又称母数列)，xi={xi(k)k=1，2，…，n}，(i=1，2，…，m)为比较数列(又称子数列)。比较数列的所有指标对应于参考数列的所有指标的关联系数λi(k)为：

式中：α是分辨系数，取值(0～1)，一般取值0.5；Δi(k)=|x0(k)-xi(k)|，为第k个指标x0与xi的绝对差。其中，min minΔi(k)为两级最小差，max maxΔi(k)为两级最大差。最终我们得到的是各比较数列与参考数列在各点的关联系数值λi(k)。

从关联系数的计算方法来看，得到的是各比较数列与参考数列在各点的关联系数值，结果较多，信息过于分散，不便于比较，为此有必要将每1个比较数列各个指标的关联系数集中体现在1个值上，这个数值就是灰色关联度，平均值法是最常用的方法。灰色关联度常记作Ri，灰色关联度越大，表明二者的几何曲线形状越接近，即发展变化趋势越接近。

1.3.5 数据处理 运用Excel 2007进行数据的整理，运用SPSS 22.0进行Pearson相关分析，检验表层土壤养分与微立地因子之间的相关关系；单因素方差分析(one-way ANOVA)检验不同坡向上，微地形之间表层土壤养分的差异，并采用LSD最小显著差异法进行多重比较(a=0.05)。运用灰色关联分析法对7种微地形及原状坡进行排序。

2 结果与分析

2.1 微立地因子与表层土壤养分指标的相关性分析

由表2可知，微地形与有机质含量、全磷含量呈极显著正相关，与全氮含量、速效磷含量、碱解氮含量显著正相关，与全钾含量呈极显著负相关；风化层厚度与有机质含量呈极显著正相关，与全磷含量、全氮含量、速效磷含量显著正相关，与全钾含量显著负相关；坡向与有机质含量、全磷含量、全氮含量显著正相关。相关分析表明，微地形、风化层厚度、坡向等微立地因子均影响着表层土壤养分含量的异质性，而其中微地形因子是引起这种异质性的主要因子。

2.2 不同坡向微地形表层土壤养分含量

2.2.1 不同微地形表层土壤有机质含量 由表3可知，在阴坡，U形沟的有机质含量最大为76.12 g/kg，塌陷和巨石背阴次之，这3种微地形有机质含量显著高于其他微地形和原状坡，坡顶的有机质含量最小，仅为26.53 g/kg。在半阴半阳坡，U形沟的有机质含量最大为71.39 g/kg，是含量最少的坡顶的2.94倍。塌陷和巨石背阴的次之，介于60.00～66.00 g/kg。在阳坡，U形沟、巨石背阴和塌陷的有机质含量较大且显著高于其他微地形和原状坡，分别是含量最小的坡顶的3.50倍、3.30倍和3.14倍。以上结果表明，无论在阴坡、半阴半阳坡和阳坡，U形沟、塌陷和巨石背阴处的有机质含量最丰富，坡顶处最低。

表2 微立地因子与表层土壤养分之间的相关分析

注：**表示在0.01水平上极显著相关，*表示在0.05水平上显著相关。

2.2.2 不同微地形表层土壤氮素含量 由表3可知，在阴坡，全氮含量在U形沟达到最大，为3.33 g/kg，塌陷和巨石背阴次之，介于2.80～3.20 g/kg，缓台的全氮含量最小，仅为1.35 g/kg，其他微地形及原状坡间的全氮含量差异不显著。U形沟的碱解氮含量最高，达到7.02 mg/kg，塌陷和巨石背阴次之分别为6.65 mg/kg和6.03 mg/kg，这3种微地形碱解氮含量显著高于其他微地形和原状坡，其中U形沟的碱解氮含量是坡顶的1.80倍。在半阴半阳坡，全氮含量在塌陷达到最大，为2.49 g/kg，各微地形及原状坡间差异不显著。U形沟、塌陷和巨石背阴的碱解氮含量显著高于其他微地形和原状坡，其中塌陷碱解氮含量最大，达到6.76 mg/kg，是含量最少的坡顶的1.82倍。在阳坡，全氮含量在U形沟达到最大，为2.89 g/kg，巨石背阴和塌陷次之，在2.50 g/kg左右，缓台的全氮含量最小，仅为1.47 g/kg。碱解氮含量在U形沟达到最大，为6.91 mg/kg，塌陷和巨石背阴次之，介于5.80～6.20 mg/kg，这3种微地形的碱解氮含量显著高于其他微地形和原状坡。其中，U形沟的碱解氮含量是坡顶的1.88倍。由此可以得出，氮素和有机质含量在各个微地形中的变化规律基本一致，土壤氮素含量取决于土壤有机质的积累和分解能力的强弱。全氮含量在U形沟、塌陷处最高，缓台处最低；碱解氮含量在U形沟、塌陷处最高，坡顶处最低。

2.2.3 不同微地形表层土壤磷素含量 由表3可知，在阴坡，全磷含量在U形沟、塌陷和巨石背阴处较大，且这3种微地形的全磷含量显著高于其他微地形和原状坡，坡顶的全磷含量最低仅为1.62 g/kg。U形沟的速效磷含量最高，达到3.56 mg/kg，塌陷、巨石背阴、陡坎和谷坡的速效磷含量次之，介于2.90～3.15 mg/kg，这5种微地形的速效磷含量显著高于坡顶和原状坡，坡顶的速效磷含量最小仅为1.28 mg/kg。在半阴半阳坡，塌陷的全磷含量最大，为3.10 g/kg，U形沟和巨石背阴次之，且3种微地形全磷含量显著高于其他微地形和原状坡，坡顶的全磷含量最小，仅为1.59 g/kg。速效磷含量在U形沟达到最大，为3.60 mg/kg，谷坡、陡坎、塌陷和巨石背阴次之，介于2.80～3.15 mg/kg，坡顶的速效磷含量最小，仅为1.32 mg/kg。在阳坡，全磷含量在U形沟达到最大，为3.13 g/kg，塌陷和巨石背阴次之，在2.90 g/kg左右，这3种微地形的全磷含量显著高于其他微地形和原状坡。速效磷含量在U形沟达到最大，为3.49 mg/kg，坡顶的速效磷含量最小仅为1.22 mg/kg。以上结果说明，U形沟处全磷和速效磷含量都是最高的，坡顶处最低。

2.2.4 不同微地形表层土壤钾素含量 由表3可知，在阴坡，全钾含量在陡坎达到最大，为13.90 g/kg，坡顶次之，而塌陷、U形沟和巨石背阴处含量较小，在9.50 g/kg左右。速效钾含量在坡顶达到最大，为62.87 mg/kg，陡坎次之为59.67 mg/kg，且这2种微地形的速效钾含量显著高于其他微地形和原状坡。在半阴半阳坡，全钾含量在坡顶达到最大，但各个微地形及原状坡间无显著差异。速效钾含量在陡坎最大，达到57.19 mg/kg，坡顶次之为56.73 mg/kg，这2种微地形的速效钾含量显著高于其他微地形和原状坡。在阳坡，全钾含量在陡坎达到最大，为13.11 g/kg，坡顶次之，塌陷的最小，仅为8.86 g/kg。速效钾含量在坡顶达到最大，为60.67 mg/kg，是含量最少的U形沟的1.40倍，陡坎次之为57.96 mg/kg，这2种微地形的速效钾含量显著高于其他微地形和原状坡。由此可以得出，坡顶、陡坎处全钾含量和速效钾含量最优，U形沟、塌陷处的最差，这与有机质和氮素含量的情形恰恰相反。

表3 不同坡向各类微地形间表层土壤养分含量(平均值±标准差)

注：不同的小写字母表示同一坡向上不同微地形间差异显著(p<0.05)。

2.3 基于灰色关联的微地形表层土壤养分评价分析

在不同坡向上，不同微地形的土壤有机质、全效养分、速效养分含量可以看出，土壤养分的各指标在不同的微地形间变化趋势不是一致的，为了定量评价不同微地形及原状坡对土壤养分改良的效果，选取土壤有机质、全氮、全磷、全钾、碱解氮、速效磷、速效钾等指标中的最大值作为参考数列，以不同微地形的土壤有机质、全效养分、速效养分等指标测定值作为比较数列，进行无量纲化处理，然后采用灰色关联分析法，计算各微地形的关联度值，关联度越大，表示微地形的比较数列与参考数列的变化趋势越接近，即该微地形对提高土壤养分的含量，改良土壤肥力的效果较好。

利用灰色关联分析方法对片麻岩山区微地形的土壤养分情况进行分析排序，结果见表4。可以看出：阴坡，U形沟的关联度值最大为0.892 4，塌陷次之，原状坡最小为0.494 2；半阴半阳坡，基本与阴坡规律一致，U形沟最大为0.852 5，坡顶最小为0.493 0；阳坡，U形沟最大为0.877 2，巨石背阴次之，坡顶最小为0.473 2。在这3个坡向上，U形沟、塌陷和巨石背阴的关联度值最高，坡顶和原状坡的最低。说明U形沟、塌陷和巨石背阴的土壤养分条件最好，而坡顶和原状坡的土壤养分条件较差。

表4 各个坡向上微地形及原状坡的灰色关联度

3 结论与讨论

微地形作为反映地表起伏变化的指标，会影响径流的产生、流向、汇流和径流量[29]，是困难立地地区植被恢复的重要影响因子，它对温度、养分、水分、光照等生态因子有再分配的作用，会导致小气候、小环境的异质性[30]。相关性分析表明，微立地因子中的微地形、风化层厚度、土壤硬度、坡向都影响有机质含量、全磷含量、全氮含量、全钾含量的变化。其中，微地形的不同是引起这个变化的首要因素。因此，微地形作为主导因子与其他立地因子共同影响着微地形养分含量的空间异质性，这与前人研究的结论一致[31]。本研究通过分析片麻岩山区不同坡向上各个微地形及原状坡间土壤养分含量的差异性得出，U形沟、塌陷和巨石背阴处有机质含量最大，显著高于其他微地形和原状坡，坡顶处最小。全氮与全磷含量规律基本一致，都在U形沟、塌陷和巨石背阴处最大，缓台和坡顶处较小。全钾含量在坡顶、陡坎、谷坡处最大，塌陷、巨石背阴处较小。碱解氮含量在U形沟、塌陷和巨石背阴处最大，显著高于其他微地形及原状坡，坡顶处最小。U形沟处速效磷含量最大，塌陷、巨石背阴、陡坎和谷坡次之，这5种微地形显著高于坡顶和原状坡，坡顶处最小。速效钾含量在坡顶、陡坎达到最大，这2种微地形显著高于其他微地形和原状坡。从整体上看，U形沟、塌陷和巨石背阴处土壤有机质、氮素、磷素的含量较高，而钾素含量最高的微地形是陡坎和坡顶，这与有机质、全氮、全磷的情况相反，说明陡坎和坡顶有利于钾素的富集[15]，这与前人的研究结论基本一致[32-34]。究其原因，有以下两点：首先，由于U形沟、塌陷等微地形有利于坡面径流的汇集，土壤的水分和养分含量比较丰富[35-36]，这在一定程度上改善了土壤的微环境，使土壤中微生物的数量和活动性增强[37]，有助于土壤中的枯落物及植物残体的分解，进而使U形沟、塌陷和巨石背阴处土壤的养分含量高于其他微地形及原状坡。其次，与微地形的植物群落特征有关，U形沟、塌陷和巨石背阴处植物群落的数量特征和多样性优于其他微地形和原状坡[25]，在土壤表层残留有大量的凋落物和死亡残体，这些物质最终分解会导致这些微地形的土壤养分含量增加。同时，植被在一定程度上能有效地减少水土流失，根系也能改善土壤的物理性质，能够使土壤的养分得到保存，这可能也是U形沟、塌陷等微地形土壤养分含量高的一个原因。

利用灰色关联分析方法对片麻岩山区微地形的养分情况进行分析排序得出：在3个坡向上，U形沟、塌陷和巨石背阴的关联度值都是最高的，坡顶和原状坡都较低。说明U形沟、塌陷和巨石背阴的养分条件最好，坡顶和原状坡的养分条件较差，这与张宏芝等[15]在黄土区的研究结论较一致。片麻岩山区是太行山典型的困难立地，土壤养分含量的高低及空间分布格局，直接影响着土壤的生产力和植被恢复的途径和方向，揭示微地形的土壤养分异质性规律，有助于指导片麻岩山区的植被恢复和生态重建。

[1] 张伟,刘淑娟,叶莹莹,等.典型喀斯特林地土壤养分空间变异的影响因素[J].农业工程学报,2013,29(1):93-101.

[2] 范夫静,宋同清,黄国勤,等.西南峡谷型喀斯特坡地土壤养分的空间变异特征[J].应用生态学报,2014,25(1):92-98.

[3] 宋述军,李辉霞,张建国.黄土高原坡地单株植物下的微地形研究[J].山地学报,2003,21(1):106-109.

[4] 赵海霞,李波,刘颖慧,等.皇甫川流域不同尺度景观分异下的土壤性状[J].生态学报,2005,25(8):2010-2018.

[5] 卢纪元.陕北植物群落结构与微地形土壤养分耦合关系研究[D].北京:北京林业大学,2016.

[6] Raulings E J, Morrisk, Roache M C, et al. The importance of water regimes operating at small spatial scales for the diversity and structure of wetland vegetation[J]. Freshwater Biology, 2010,55(3):701-715.

[7] 董云龙,张德罡,陈建纲,等.东祁连山高寒草地不同生境条件下植物群落结构特征及多样性比较:以甘肃省天祝抓喜秀龙乡为例[J].草地学报,2014,22(3):481-487.

[8] Jessia E, Andraw H. Seedling emergence from seed banks of tidal fresh water wetlands:Response to inundation and sedimentation[J]. Aquatic Botany,2004,78(3):243-254.

[9] 卢纪元,朱清科,陈文思,等.陕北黄土区植被特征对坡面微地形的响应[J].中国水土保持科学,2016,14(1):53-60.

[10] 杨更强,严成,陈军纪,等.砾漠区微地形下梭梭群落植被特征[J].东北林业大学学报,2015,43(5):64-67.

[11] 戚德辉,温仲明,王红霞,等.黄土丘陵区不同功能群植物碳氮磷生态化学计量特征及其对微地形的响应[J].生态学报,2016,36(20):1-11.

[12] Lin H, Wheeler D, Bell J, et al. Assessment of soil spatial variability at multiple scales[J]. Ecological Modelling, 2005,182(3):271-290.

[13] 刘旻霞,王刚.高山草甸坡向梯度上植物群落与土壤中的N, P化学计量学特征[J].兰州大学学报:自然科学版,2012,48(3):70-75.

[14] 高凯,张丽娟,于永奇,等.沙地土壤pH值、碳、氮、磷含量对微地形变化的响应[J].水土保持通报,2016,36(1):88-92.

[15] 张宏芝,朱清科,赵磊磊,等.陕北黄土坡面微地形土壤化学性质[J].中国水土保持科学,2011,9(5):20-25.

[16] 赵秀芳,张清,王振宇.微地形对天津滨海吹填土土壤理化性质和植被状况的影响[J].土壤通报,2014,45(2):281-285.

[17] 刘孝阳,周伟,白中科,等.平朔矿区露天煤矿排土场复垦类型及微地形对土壤养分的影响[J].水土保持研究,2016,23(3):6-12.

[18] 罗光杰,李阳兵,王世杰,等.岩溶山区景观多样性变化的生态学意义对比:以贵州4个典型地区为例[J].生态学报,2011,31(14):3882-3889.

[19] 张忠华,胡刚,祝介东,等.喀斯特常绿落叶阔叶混交林物种多度与丰富度空间分布的尺度效应[J].生态学报,2012,32(18):5663-5672.

[20] 李保国,张金柱,郭素平.河北省太行山片麻岩山地景观生态经济资源特征与利用研究[J].河北农业大学学报,2005,28(4):14-17.

[21] Kikuchi T. Vegetation and Land Forms [M]. Tokyo:University of Tokyo Press, 2001.

[22] Nagamatsu D, Miura O. Soil disturbance regime in relation to micro-scale landforms and its effects on vegetation structure in a hilly area in Japan[J]. Plant Ecology, 1997,133(2):191-200.

[23] 张宏芝,朱清科,王晶,等.陕北黄土坡面微地形土壤物理性质研究[J].水土保持通报,2011,31(6):55-58.

[24] 杨永川,达良俊,由文辉.浙江天童国家森林公园微地形与植被结构的关系[J].生态学报,2005,25(11):2830-2840.

[25] 王晶,朱清科,秦伟,等.陕北黄土区封禁流域坡面微地形植被特征分异[J].应用生态学报,2012,23(3):694-700.

[26] 汝海丽,张海东,焦峰,等.黄土丘陵区微地形对草地植物群落结构组成和功能特征的影响[J].应用生态学报,2016,27(1):25-32.

[27] 史薪钰,陈梦华,齐国辉,等.片麻岩山地土壤团聚体分布及养分特征[J].水土保持学报,2015,29(5):172-175.

[28] 唐柄哲,何丙辉,闫建梅.川中丘陵区土地利用方式对土壤理化性质影响的灰色关联分析[J].应用生态学报,2016,27(5):1445-1452.

[29] 塔娜,王健,张慧荟,等.黄土耕作坡面溅蚀过程中微地形响应特征[J].水土保持通报,2016,36(1):110-114.

[30] Forrester D, Bauhus J, Connell M. Competition in thinned Silvertop Ash(EucalyptussieberiL. Johnson)stands from early coppice growth[J]. Forest Ecology and Management, 2003,174(1):459-475.

[31] 卜耀军,朱清科,包耀贤,等.陕北黄土区微地形土壤质量指标变异性及敏感性[J].水土保持学报,2014,28(3):153-157.

[32] 李亨伟,胡玉福,邓良基,等.川中丘陵区紫色土微地形下有机质空间变异特征[J].土壤通报,2009,40(3):552-554.

[33] 徐明,张健,刘国彬,等.黄土丘陵区不同植被恢复模式对沟谷地土壤碳氮磷元素的影响[J].草地学报,2015,23(1):62-68.

[34] 张伟,陈洪松,王克林,等.典型喀斯特峰丛洼地坡面土壤养分空间变异性研究[J].农业工程学报,2008,24(1):68-73.

[35] 赵荟,朱清科,秦伟,等.黄土高原干旱阳坡微地形土壤水分特征研究[J].水土保持通报,2010,30(3):64-68.

[36] 邝高明,朱清科,刘中奇,等.黄土丘陵沟壑区微地形对土壤水分及生物量的影响[J].水土保持研究,2012,19(3):74-77.

[37] 白爱芹,傅伯杰,曲来叶,等.重度火烧迹地微地形对土壤微生物特性的影响:以坡度和坡向为例[J].生态学报,2013,33(17):5201-5209.

DifferentiationCharacteristicsofTopsoilNutrientsoftheMicro-topographyinGneissMountainousAreaofPingshanCounty,HebeiProvince

YUAN Zhen, WEI Songpo, JIA Liming, ZHANG Yaxiong, LIU Longlong

(Province-MinistryCo-constructedKeyLaboratoryofSilvicultureandConservation,MinistryofEducation,BeijingForestryUniversity,Beijing100083,China)

Based on the determination data of topsoil nutrients in gneiss mountainous area in Pingshan County of Hebei Province, we analyzed the differences of topsoil nutrients on seven micro-topographies including slope crest, collapse, stone shadow, platform, scarp, brae and U-gully, with the undisturbed slope as the control, in order to provide scientific reference for the revegetation of gneiss mountainous area. The results showed that: (1) according to correlation analysis, micro-topography, aspect and the weathered layer thickness influenced the change of topsoil nutrients, among which micro-topography was the main factor; (2) the topsoil nutrient conditions on most of the micro-topographies were better than the undisturbed slope; U-gully, collapse, stone shadow contained significant higher organic matter contents (68.94 g/kg, 60.71 g/kg, 62.77 g/kg) than other micro-topographies and the undisturbed slope, and the lowest on slope crest (22.58 g/kg); (3) changes of total nitrogen content and total phosphorus content are similar, and the thire contents were highest on U-gully, collapse, stone shadow, the lowest on platform, slope crest; total potassium contents on slope crest, scarp, brae were the highest (11.88 g/kg, 12.56 g/kg, 11.67 g/kg), while the least contents (9.39 g/kg, 9.69 g/kg) were observed on collapse, stone shadow; (4) U-gully, collapse, stone shadow contained significantly higher available nitrogen contents (6.87 mg/kg, 6.54 mg/kg, 5.93 mg/kg) than the other micro-topographies and the undisturbed slope, the lowest on slope crest (3.75 mg/kg). U-gully had the highest available phosphorus content (3.55 mg/kg), followed by collapse, stone shadow, scarp and brae; contents of available phosphorus on the five micro-topographies were significantly higher than the undisturbed slope and slope crest(1.27 mg/kg) which had the lowest available phosphorus content; available potassium contents on slope crest (60.09 mg/kg) and scarp (58.27 mg/kg) were significantly higher than other micro-topographies and the undisturbed slope; (5) through the grey correlation analysis on the topsoil nutrients of seven micro-topographies, it can be concluded that the conditions of topsoil nutrients on U-gully, collapse, stone shadow were best, worst on slope crest and the undisturbed slope.

gneiss mountainous area; micro-topography; soil nutrients; spatial variability; grey correlation analysis

S151.9

A

1005-3409(2017)06-0084-07

2016-12-01

2016-12-18

“十二五”国家科技支撑计划项目(2015BAD07B02)困难立地植被恢复技术研究与示范

袁振(1990—),男,河南项城人,硕士研究生,主要从事困难立地植被恢复。E-mail:279938085@qq.com

贾黎明(1968—),男,山西忻州人,教授,博士生导师,主要从事森林培育理论与技术。E-mail:jlm@bjfu.edu.cn