北京松山国家级自然保护区地处延庆县西北部，森林覆盖率接近90%。地处燕山山脉的军都山，坐标为40°29′9″~40°33′35″N，115°43′44″~115°50′22″E^[12]。属于暖温带大陆性季风气候区，年平均气温为8.5 ℃，最高和最低气温分别为39.0 ℃和-27.3 ℃，年降水量与蒸发量分别为493.0 mm和1 772.0 mm，年平均日照时数为2 836.3 h。随海拔由低向高依次分布有山地褐土、山地棕壤和山地草甸土^[12]，土层厚度为20~55 cm。保护区总面积6 264.00 hm²。通过对调查数据的处理，在松山国家级自然保护区一共发现胡桃楸斑块65块，胡桃楸林总计为396.47 hm²。主要分布在大庄科村、兰角沟等地区的沟谷内阴坡一侧，海拔高度为700~1 500 m，以天然纯林为主，萌生情况普遍，少量伴生黑枣Diospyros lotus，榆树Ulmus pumila和大叶白蜡Fraxinus rhynchophylla等。林龄组成以中龄林为主，成熟林和近熟林次之，幼龄林最少。

采用典型性样方法^[13]，于2015年7月至8月，在兰角沟、大庄科村附近、公路两侧主要胡桃楸林分布区选取立地条件相近、林龄(以40年生为林龄期限，用生长锥取得树龄数据)不同的典型林分设置20 m × 20 m的乔木样方24个(图 1)，并在每个乔木样方四角及中心位置各做5 m × 5 m灌木样方1个和1 m × 1 m草本样方1个。样方内记录：① 生态因子，如海拔、坡度、坡向、坡位和土壤情况；② 乔木树种的名称、胸径、高度以及整个群落的郁闭度；③ 灌木和草本植物的名称、数量、高度及盖度^[14]。

图  1  胡桃楸林样方分布图

Figure 1.  Distribution pattern of Juglans mandshurica forests

根据样方内植物个体高度、盖度、密度以及频度，计算灌木和草本物种的重要值。公式为：灌木植物重要值=(相对密度+相对高度+相对盖度)/3；草本植物重要值=(相对频度+相对密度+相对高度+相对盖度)/4^[5]；相对频度=该种的频度/所有种的频度之和，相对密度=该种的个体数/所有种的个体数总和，相对高度=该种的平均高度/所有种的平均高度总和，相对盖度=该种的盖度/全部种的盖度总和^{[5, 15-17]}。总结前人研究成果，根据多样性指数对反映群落物种多样性状况的实际应用效果，本研究选取了以下4种多样性指数^{[5, 15-17]}。Simpson指数(D)：$D = 1 - \sum\limits_{i = 1}^S {P_i^2}$ 。Shannon-Wiener指数(H)：$H = - \sum\limits_{i = 1}^S {{P_i} \times \ln {P_i}}$ 。Pielou均匀度指数(J)：J=H/lnS。Margalef丰富度指数(M)：M=(S-1)/lnN。其中：S为样方内物种数；N为样方内所有物种个体数总和；P_i为种i的相对重要值，等于种i的重要值占样方中所有种的重要值之和的比例。

利用SPSS 22.0软件进行均值分析，并对不同林龄胡桃楸林植物多样性指数进行单因素方差分析；利用SigmaPlot 12.0对结果进行图形化处理；利用GS+9.0软件对群落物种多样性空间分布进行半变异函数与方差模型分析。

调查记录胡桃楸林下植物43科78种，以菊科Asteraceae，蔷薇科Rosaceae，豆科Leguminosae，毛茛科Ranunculaceae和百合科Liliaceae等为主，其中木本植物18科26种，草本植物25科52种。同时还对胡桃楸林内灌木层和草本层植物的重要值进行了统计(表 1与表 2)。

灌木层中，葎叶蛇葡萄Ampelopsis humulifolia，小叶鼠李Rhamnus parvifolia，红花锦鸡儿Caragana rosea和山楂叶悬钩子Rubus crataegifolius在不同林龄林下都有分布。红花锦鸡儿与小叶鼠李在4种林龄林下都有较大优势，而山楂叶悬钩子和葎叶蛇葡萄则存在差异。荆条Vitex negundo var. heterophylla和山杏Armeniaca sibirica分布在幼龄林及中龄林中，毛叶丁香Syringa pubescens，小花溲疏Deutzia parviflora，暴马丁香Syringa reticulata var. mandshurica和土庄绣线菊Spiraea pubescens等在近熟林与成熟林中分布，细叶小檗Berberis poiretii和大花溲疏Deutzia grandiflora仅在成熟林中可见。整体而言，幼龄林及中龄林的灌木种类较少，随着林龄的增加，灌木层植物种类逐渐增加，存在明显的变化趋势。保护区内胡桃楸林下灌木群落以小花溲疏、毛叶丁香、小叶鼠李、红花锦鸡儿等为主要优势种。

草本层中鸭跖草Commelina communis在不同林龄下均有分布，但重要值存在差异。不同林龄的胡桃楸林下草本群落组成差异显著，随着林龄增加，大油芒Spodiopogon sibiricus，披针叶薹草Carex lancifolia，蝙蝠葛Menispermum dauricum等逐渐被北重楼Paris verticillata，鼠掌老鹳草Geranium sibiricum，蓝萼香茶菜Rabdosia japonica等取代，形成较稳定的顶级群落。优势种上，幼龄林和中龄林主要是披针叶薹草，大油芒，草乌Aconitum carmichaeli和龙芽草Agrimonia pilosa等，近熟林和成熟林则以蓝萼香茶菜和白屈菜Chelidonium majus等为主。整体而言，随着林龄的增加，草本层植物种类减少的趋势明显。

通过单因素方差分析，对不同林龄胡桃楸林下灌木层物种多样性进行差异显著性分析。结果显示：不同演替阶段的灌木层Simpson指数存在显著差异(P＜0.05)，随着林龄增加，灌木层物种变化明显，指数值不断增加(图 2和表 3)；而Shannon-Wiener多样性指数和均匀度指数分别呈极显著(P＜0.001) 和显著(P＜0.05) 的差异，随着林龄的增加而逐渐增加。这表明，随着演替阶段的深入，灌木的物种多样性表现出不断递增的规律。

图  2  不同林龄胡桃楸林下灌木层物种多样性指数

Figure 2.  Species diversity index of shrub layers under forest at different ages of manchurian walnut plantations

单因素方差分析显示，不同林龄下草本层的Simpson多样性指数存在显著差异(P＜0.05)，并随林龄的增加，呈现逐渐降低的规律(图 3与表 3)；Shannon-Wiener多样性指数差异更为显著(P＜0.001)，随着林龄增加，多样性指数下降的趋势更加明显；不同演替阶段的草本层均匀度指数也存在显著差异(P＜0.001)，但变化趋势是上升的。这表明，随着演替阶段的深入，草本植物有逐渐退出生态系统的趋势，最后形成一个稳定的顶级群落，这种变化规律与Bohanek等人提出的观点是一致的^[17]。

图  3  不同林龄胡桃楸林下草本层物种多样性指数

Figure 3.  Species diversity index of herb layers under forest at different ages of manchurian walnut plantations

通过单因素方差分析，灌木层的物种丰富度差异极显著(P＜0.001)，与幼龄林相比，中龄林灌木层物种丰富度增加了20%以上，近熟林比中龄林又增加了近43%(图 4)，最后在成熟林阶段达到了最高值，但与近熟林相比，成熟林的物种丰富度增加比例不到9%。这表明：从胡桃楸群落演替早期到达到稳定状态的期间，物种丰富度先快速增加，到中期增速达到最快，而到了后期阶段，随着群落逐渐稳定，植物物种的替换速率也随之下降。草本层的变化规律与灌木层呈现完全相反的趋势。随着胡桃楸林龄的增加，林下草本层植物的种类逐渐减少(图 4)。不同林龄胡桃楸林下草本层物种丰富度存在极显著差异(P＜0.001)，在演替初期，大量的草本植物出现，形成早期的先锋群落，但随着林龄增加，群落中草本植物不断退出群落，尽管有部分新的植物进入，但在整体上，最后的稳定状态的顶级群落与早期相比，草本层物种丰富度下降非常明显。

图  4  胡桃楸林下植物物种丰富度指数

Figure 4.  Species richness index of undergrowth under forest of manchurian walnut plantations

从表 4物种多样性半变异函数模型及方差分析结果可知，松山胡桃楸林4种多样性指数中，均匀度指数和Simpson多样性指数为球面模型，Shannon-Wiener多样性指数和丰富度指数则更符合指数模型，半变异函数模型能够很好地反映物种多样性指数的空间结构特性。在分析模型中块金值表征由于取样误差和小于取样尺度上随机引起的空间异质性，基台值越大，总的空间异质性越高；块金与基台的比值表示由系统变量引起的空间相关性，比值小于0.25表示空间相关性强烈，0.25~0.75表示中等的空间相关性，大于0.75表示空间相关性较弱^[18]；4种多样性指数比值都为0.25~0.75，因此研究区的物种多样性为中等的空间相关性。

同时也说明：在大于研究尺度，小于变程的范围内结构性因子所占比例为0.584~0.681，随机性因子所占比例为0.319~0.416，说明尽管结构性因子仍然占主导地位，但诸如人为砍伐、旅游开发等人为活动引起的空间异质性起着非常重要的作用，人为活动的干扰对胡桃楸植物多样性有着显著影响。

灌木层中红花锦鸡儿、小叶鼠李等在各林分都属于优势度较大的树种，草本层中大油芒、披针叶薹草等属于优势植物。随着林龄的增加，草本植物物种数逐渐下降，可能是由于林下生存空间及养分的竞争加剧，竞争力差的植物逐渐淘汰出局^[19-21]。葎叶蛇葡萄、东亚唐松草、鸭跖草等几乎在胡桃楸各个林龄阶段都有出现，说明这些物种非常适应胡桃楸林下的环境^[22]。由于不同物种对特定的林下生境的适应能力各不相同，群落中的优势种类随演替阶段不断发生变化。整体而言，演替早期胡桃楸草本群落主要以披针叶薹草、大油芒、草乌、龙芽草等为优势种，中后期则以蓝萼香茶菜、白屈菜等为优势种；灌木群落以小花溲疏、毛叶丁香、小叶鼠李、红花锦鸡儿等为优势种。

杨立学等^[10]认为：胡桃楸林龄可以对植物群落的完整性、抗逆性及演替的进程等状况进行真实反映。随着林龄的增加，胡桃楸林下灌木层物种多样性指数呈现连续增加的趋势，而草本层的Simpson指数、Shannon-Wiener指数、物种丰富度指数呈现持续递减的规律，均匀度指数则与之相反，表现为连续递增趋势。这表明，随着胡桃楸的干扰逐渐增强，林下灌木层和草本层的响应并不一致，适应能力较强的灌木在群落中稳定下来，适应能力较弱的草本植物则有退出群落的趋势；而幼龄林林下植被的均匀度指数明显不如成熟林高，可能是因为到了成熟林阶段，林分结构逐渐稳定下来，外部环境相对比较固定，对林下植物生长与分布的限制作用趋于减弱。

胡桃楸不同演替阶段植物多样性之所以存在显著差异，主要源于在以胡桃楸为优势树种的林分中，随着林龄的增加，胡桃楸在群落中的优势地位增强，造成林下植物对有限的空间和资源的竞争更为激烈，只有适应能力更强的个体才得以生存，部分植物逐渐趋于淘汰、消亡，最终完成胡桃楸林下群落演替过程。

自然状态下，胡桃楸的化感作用使其群落保持优势，但林龄的变化会直接导致林地中胡桃醌的含量变化，最终影响胡桃楸林群落结构的维持^{[11, 23]}。松山胡桃楸林多分布于沟谷、道路两旁的区域，这些地区是保护区发展生态旅游的重点区域，景区设施与道路的建设，加上村民对胡桃楸种子的采集，造成很多地方胡桃楸天然实生更新困难。人为活动对胡桃楸林植物多样性产生了显著影响，胡桃楸的自然更新与维持受到了威胁，因此，科学合理的规划景区建设及道路的开发，有效地对当地居民的种子采集等活动进行管理，就显得十分重要。