2024年6月14日发(作者:慕容翠柏)
第44卷第7期
2024年4月
ACTAECOLOGICASINICA
生态学报
Vol.44,No.7
Apr.,2024
DOI:10.20103/j.stxb.202305060941
史恭发,徐诺,牛钊倩,孙炜伦,王好,石广振,王玲.内蒙古大兴安岭东部林下植物生物多样性海拔差异.生态学报,2024,44(7):3004⁃3015.
Mountains,InnerMongolia.ActaEcologicaSinica,2024,44(7):3004⁃3015.
ShiGF,XuN,NiuZQ,SunWL,WangH,ShiGZ,WangL.AltitudinaldifferencesofunderstoryplantbiodiversityineasternGreaterXing′an
内蒙古大兴安岭东部林下植物生物多样性海拔差异
史恭发,徐 诺,牛钊倩,孙炜伦,王 好,石广振,王 玲
∗
东北林业大学园林学院,哈尔滨 150040
摘要:海拔影响生物多样性分异,为探究寒温带山岭地区海拔对林下植物多样性分布影响,选择内蒙古大兴安岭东部地区,采用
样地调研法,设置E1(200—350m),E2(350—500m),E3(500—650m),E4(650—800m),E5(800—950m),E6(950—1100m)
6个不同海拔梯度,调研165个林地样点以了解不同海拔森林群落物种组成,研究林下植物生物多样性海拔差异性、地形因素、
乔木层优势物种与林下植物生物多样性的相关性。研究得出:(1)调研区林下植物共277种隶属于53科135属,灌木32种、草
白桦⁃欧亚绣线菊⁃绣线菊群落(E4)、落叶松⁃白桦⁃兴安杜鹃⁃越橘群落(E5)、落叶松⁃白桦⁃越橘⁃偃松群落(E6);(3)林下植物群
落生物多样性随海拔上升呈现明显的单峰格局。生物多样性草本层>灌木层,灌木层随海拔升高呈上升趋势,草本层呈先下降
后升高再下降的波动下降趋势;(4)草本层物种替换速率远高于灌木层,E3—E4与E2—E3梯度内灌木层与草本层物种替换速
率分别达到峰值;灌木层与草本层相邻海拔梯度间相似度指数均呈先下降后升高再下降的波动形式;(5)地理因素及优势乔木
种不同能够显著影响林下植物α多样性,灌木层α多样性与海拔、纬度呈正相关,与经度、坡度呈负相关。草本层与海拔呈负
相关,与纬度、坡向呈正相关,海拔对草本层α多样性影响高于灌木层。郁闭度与灌木层α多样性呈正相关,与草本层α多样
性呈负相关;落叶松增多可提升灌木层生物多样性,黑桦会降低灌木层物种丰富度,而山杨会降低草本层物种分布的均匀度。
研究结果为保护与利用内蒙古大兴安岭林下植物,丰富寒温带山岭地区林下植物多样性海拔分布理论提供科学依据。
关键词:内蒙古大兴安岭;林下植物群落;生物多样性;海拔
本245种;(2)森林乔灌群落随海拔依次为蒙古栎⁃黑桦⁃榛⁃胡枝子群落(E1—E2)、落叶松⁃白桦⁃榛⁃胡枝子群落(E3)、落叶松⁃
AltitudinaldifferencesofunderstoryplantbiodiversityineasternGreaterXing′an
Mountains,InnerMongolia
SHIGongfa,XUNuo,NIUZhaoqian,SUNWeilun,WANGHao,SHIGuangzhen,WANGLing
∗
CollegeofLandscapeArchitecture,NortheastForestryUniversity,Harbin150040,China
Abstract:Inordertoexploretheeffectofaltitudeonthediversityofundergrowthplantsinthecoldtemperatemountain
area,weselectedtheeasternGreaterXing′anMountainsinInnerMongoliabyusingtheplotinvestigationmethod,andset
(800—950m),andE6(950—1100m).Weinvestigated165forestsitestounderstandthecompositionofforest
communitiesatdifferentaltitudes,andstudiedtherelationshipbetweenthealtitudedifferenceinunderstoryplant
sixdifferentaltitudegradients:E1(200—350m),E2(350—500m),E3(500—650m),E4(650—800m),E5
biodiversity,topographicfactors,anddominantspeciesinthearboreallayerandunderstoryplantbiodiversity.Theresults
245herbs.(2)TheforesttreeandshrubcommunitieswithdifferentaltitudeswereintheorderofQuercusMongolia,Betula
heterophylla,Lespedezabicolorcommunities(E3),Larixgmelina,Betulaplatyphylline,Spiraeamedia,Spiraeasilicified
基金项目:国家科技基础资源调查专项(2019FY100500);中央高校基本科研业务费专项资金(2572023CT18)
收稿日期:2023⁃05⁃06; 网络出版日期:2024⁃01⁃12
∗通讯作者Correspondingauthor.E⁃mail:wanglinghlj@126.com
showedthat:(1)therewere277understoryplantspeciesbelongingto53familiesand135genera,including32shrubsand
daturic,Corylusheterophylla,Lespedezabicolorcommunities(E1—E2),Larixgmelina,Betulaplatyphylline,Corylus
http://www.ecologica.cn
7期 史恭发 等:内蒙古大兴安岭东部林下植物生物多样性海拔差异
3005
communities(E4),Larixgmelina,Betulaplatyphylla,Rhododendrondauricum,Vacciniumvitis⁃idaeacommunities(E5),
Larixgmelinii,Betulaplatyphylla,Vacciniumvitis⁃idaea,Pinuspumilacommunities(E6).(3)Thebiodiversityof
understoryplantcommunitiesshowedanobviouslyunimodalpatternwiththeincreasesofaltitude.Theherblayerwashigher
fluctuatingdecreasetrendoffirstdecreasingandthenincreasinganddecreasing.(4)Thereplacementrateofspeciesin
thantheshrublayer.Theshrublayershowedanincreasingtrendwiththeincreasingaltitude,whiletheherblayershoweda
herbaceouslayerwasmuchhigherthanthatinshrublayer,withthereplacementrateinshrubandherbaceouslayerpeaking
thenincreased,andfinallydecreased.(5)Thegeographicalfactorsanddifferentdominanttreespeciescouldsignificantly
andnegativelycorrelatedwithlongitudeandslope.Theherbaceouslayerwasnegativelycorrelatedwithaltitude,but
inE3—E4andE2—E3gradients,respectively.Thesimilarityindexbetweenshrublayerandherbagelayerdecreasedfirst,
affecttheαdiversityofunderstoryplants.Theαdiversityofshrublayerwaspositivelycorrelatedwithaltitudeandlatitude,
positivelycorrelatedwithlatitudeandslopeaspect.Theeffectofaltitudeonαdiversityofherbaceouslayerwashigherthan
thatofshrublayer.Thecanopydensitywaspositivelycorrelatedwiththeαdiversityofshrublayerandnegativelycorrelated
withtheαdiversityofherblayer.Larixgmeliniiincreasedbiodiversityinshrublayer,Betuladahuricadecreasedspecies
richnessinshrublayer,andBetuladahuricadecreasedspeciesdistributionevennessinherbaceouslayer.Ourresults
providedscientificbasisfortheprotectionandutilizationofundergrowthplantsintheGreaterXing′anMountainsandthe
theoryofaltitudinaldistributionofunderstoryplantdiversityinthecoldtemperatemountains.
KeyWords:GreaterXing′anMountainregionofInnerMongolia;understoryplantcommunity;bio⁃diversity;altitude
pattern
降水等多方面环境因子的复合变化,对群落分布与物种多样性产生重要影响
[3]
。山地是具有一定海拔和坡
度变化的地面
[4]
,其具有浓缩的环境梯度、高度异质化的生境
[5]
,通常具有明显的生物多样性与群落的垂直
地带分异
[6]
,始终是生态学研究热点
[7]
,山地生物的垂直变化是水平梯度的近千倍
[8]
,山地垂直带被认为是
水平地带的“微缩模型”
[9]
。植被的垂直带性虽在形态上与水平地带性具有一定的相似性,但绝非等同水平
地带的重复
[4,10]
,山地生态系统植物多样性随海拔上升一般有5种形式,先升后降的单峰格局
[11]
与单调下降
的格局较为普遍
[12]
,此外还存在单调升高、先降低后升高、和不具相关性3种格局
[12]
。根据温度和水分差
异,我国山地垂直自然带划分为大陆性和季风性两类带谱体系
[13]
,不同地区、不同生活型的物种多样性海拔
分布格局具有差异性,研究的区域不同得到的结果不同,因此对不同地区植物多样性海拔格局展开研究十分
必要
[12]
。
不同尺度区域下,植物群落分布的影响因子不同,小尺度区域植物群落主要由地形、土壤等方面影响,大
生物多样性的海拔差异是生态学研究的重要问题
[1—2]
,海拔作为综合性驱动因子,能够引起光照、温度、
尺度区域主要由气候、地势等方面决定
[14]
。许多单一山体按一定方向呈带状延伸形成山岭,相比于单个山
体,山岭地区具有经纬度跨度更大的水平地带性与干湿度地带性,具有海拔变化梯度更大的垂直地带性,因此
具有更丰富的“三维空间”分布现象,此外还具有人类干扰程度较小,对气候环境反应灵敏等特点,彼此之间
相互交叉过渡,更容易形成生物多样性的聚集地。对于山岭地区植被随海拔梯度的分布规律相对于单一山体
更能得出山地植被垂直分布的一般规律性。我国植物多样性海拔梯度的研究集中于单独山体
[15—18]
,并且近
一个世纪山地植被地带性垂直分布变化格局的相关研究集中于我国横断山区、青藏高原、太白山及秦岭区
域
[13]
,十分缺乏东北山岭地区植被地带性垂直分布格局的相关研究。
大兴安岭位于大陆性气候和温带季风气候区分界线与我国地势二三级阶梯分界线,具有较大的海拔差
异。内蒙古大兴安岭林区占整个大兴安岭的46%,森林面积8.37万km
2[19]
并以天然林为主,包括110万hm
2
未开发的原始森林,森林覆盖率78.39%,居国有林区之首,在东北地区具有丰富的植物资源。林下植物是森
林生态系统的重要组成部分,在维持生态系统稳定性、促进群落演替更新、保证人工林生物多样性与生态功能
http://www.ecologica.cn
3006
生 态 学 报 44卷
可持续发展等方面具有重要作用
[20—21]
。林下植物蕴含大量药用、食用、观赏等方面的经济植物类群
[22—24]
,在
推动当地林下产业经济发展方面具有重大潜力。了解林下植物的物种组成和生物多样性梯度分布是进一步
开展保护与开发内蒙古大兴安岭林下植物群落的先期工作。
关系
[26—27]
,林下经济植物资源分布与评价
[22,28—31]
,落叶松林下植物多样性
[32]
与火烧迹地植物多样性
[33]
等方
面。本文选取内蒙古大兴安岭东部林下植物多样性海拔分布差异为研究对象,通过200—1100m连续分布样
地调研方法减少因局部小气候等因素产生的影响,根据重要值、α多样性指数与拟合曲线分析、β多样性指
数,林下植物多样性与地理因素及优势乔木种的相关性等方法研究了内蒙古大兴安岭东部森林植被物种组成
与林下植物多样性沿海拔梯度的分布规律,将灌木层与草本层相比较,分析两种群落层次间差异,考虑乔木层
优势物种比重对林下植物多样性的影响,旨在深入了解内蒙古大兴安岭东部森林群落随海拔分布、林下灌木
与草本物种多样性的海拔分布、地理因子与乔木优势物种组成对林下植物多样性的影响等问题。丰富了我国
寒温带山岭森林植被垂直地带性分布的研究内容,填补了我国东北地区山岭植被垂直地带性分布的研究空
缺,为了解大陆性和季风性两类带谱过渡地区垂直自然带谱分布特征及其规律,保护与利用内蒙古大兴安岭
林下植物,开展不同海拔梯度人工林营造等工作提供科学依据。
1 研究区概况
研究区位于内蒙古大兴安岭林区东部莫力达瓦达斡尔族自治旗(莫旗)与鄂伦春族自治旗(鄂旗),为北
半球中高纬度地区(48°64′—51°13′N,122°05′—125°14′E),位于大兴安岭东南坡。调研区属寒温带大陆性
季风气候,冬寒夏暖,昼夜温差大,1月份平均气温
-
24—
-
31℃,7月份平均气温16—20℃,平均年降水量
300—500mm,河流较多,以嫩江水系为主;地势沿东南至西北方向海拔迅速升高,坡度逐渐增加,海拔164—
1484m之间,具有较大落差,以山地、丘陵、平原等地貌为主。
2 研究方法
2.1 样地设置与调研方法
调研使用典型取样与随机取样相结合的样地调研方法,在内蒙古大兴安岭东部莫旗及鄂旗两地沿海拔
200—1100m范围分布的森林设置6个不同海拔梯度样带:E1:200—350m,E2:350—500m,E3:500—650m,
E4:650—800m,E5:800—950m,E6:950—1100m。根据实际可达性,选取165个森林样点开展实地调查
(图1)。
调查方法以王纯华等人
[26,34—36]
为参考,根据本调研目的做出调整。每个样点在半径500m范围内设置
目前对于内蒙古大兴安岭植物研究多集中于典型森林群落多样性的纬度格局
[25]
,多样性与群落结构间
3个标准地(即3个重复),尽可能涵盖样点区域不同林地位置与方向以获取全面的物种分布信息,每个标准
地由1个20m
×
20m乔木样方及其东北角1个5m
×
5m灌木样方、4个分布于4角的1m
×
1m草本样方组
成,总调研495个标准地,调研面积198000m
2
。记录每个样点的经纬度、海拔、坡度、坡向、郁闭度、林型等信
息(表1)。坡向数据以丛静等人
[35]
方法进行赋值,数字越大越向阳干热。乔木、灌木与草本样方中记录物种
名、株数等;对各调研样点生境、林下植物均拍照记录。现场鉴定参考第二版《东北植物检索表》
[37]
,借鉴《额
尔古纳自然保护区生物多样性》植物名录辅助鉴定
[38]
,最终植物正名参考中国植物志(网址:http://www.
2.2 数据处理
iplant.cn/frps)。调研时间为2022年6月—2022年9月。
运用AdobeLightroomClassic管理样地植物照片,辅助植物鉴定与统计;运用OfficeExcel(2019)汇总数据
并计算α与β多样性;使用IBMSPSSStatistics26进行单因素方差分析、Pearson相关性分析与线性回归分析;
运用Origin(2021)生成α多样性与海拔散点分布图并使用SimpleFit插件生成α多样性与海拔的拟合方程。
运用最小显著差异法检验群落相应条件下多样性指数在置信区间(P<0.05)水平上的差异显著性。
http://www.ecologica.cn
7期 史恭发 等:内蒙古大兴安岭东部林下植物生物多样性海拔差异
3007
图1 调研样点
Fig.1 Surveysamplepoint
表1 研究区6个海拔梯度调研样地特征
Table1 Characteristicsofsixaltitudegradientsurveysamplesitesinthestudyarea
海拔等级
Altitudelevel
E1
E2
E3
E4
E5
E6
海拔梯度
Altitudegradient/m
200—350
350—500
500—650
650—800
800—950
950—1100
样点数量
Numberofsample
points
30
35
34
27
28
11
平均乔木高
Averagetree
height/m
9.24
9.68
平均郁闭度
Averagecanopy
density
0.54
0.62
0.67
0.70
0.66
0.50
平均灌木高
Averageshrub
height/m
0.95
0.99
1.08
0.94
0.80
0.27
12.21
14.99
13.87
9.48
2.2.1 重要值
[39]
V
=
(相对盖度
+
相对频度
+
相对密度)/3
2.2.2 α多样性
(1)物种丰富度指数(Patrick指数)(R)
[20,40]
,指数值越高说明物种丰富度越高。
R
=
S
http://www.ecologica.cn
3008
生 态 学 报 44卷
(2)物种多样性指数Shannon⁃Wiener(H′),指数值越高说明物种多样性越高
[33,40]
。
H′
=-
(3)辛普森指数(优势度指数)Simpson指数值(D),指假设在无限个群落随机抽取样本,样本中两个不同
P
i
lnP
i
∑
i
=
1
S
0—1之间。
物种个体相遇的概率为多样性测度
[20,26]
,指数值越高说明物种多样性越高,而物种优势度越低。数值范围在
D
=
1
-
P
i
2
∑
i
=
1
S
内分布越均匀。数值范围在0—1之间。
(4)均匀度指数Pielou指数(J),表示与物种丰富度无关的均匀度指数
[7,40]
,指数值越高说明物种在样方
S表示研究系统中物种总数,P
i
表示第i个物种的相对密度(P
i
=
N
i
/N,N为所存在群落全部种的个数之
J
=
H′/lnS
2.2.3 β多样性
和,N
i
代表第i个种的个体数)。
(1)Cody指数β
C
[1]
β
C
=
g(H)为沿海拔梯度增加的物种数目;l(H)为沿海拔梯度失去的物种数目,β
C
数值越高表示海拔梯度间物
种替换速率越快。
(2)群落相似度指数C
S
[1]
C
S
=
g
(
H
)
+
l(H)
2
2c
a
+
b
c为两群落共有的物种数,a和b分别为两群落各自物种数,C
S
数值越高表示海拔梯度间物种相似度越
高。C
S
值在0.8—1.0之间为高度相似,0.6—0.8为中度相似,0.4—0.6为低度相似,0.2—0.4为非常低的相
似,0—0.2表示海拔梯度间物种几乎没有相似之处。
3 结果分析
3.1 林下植物物种组成
调研区林下植物共277种隶属于53科135属,灌木32种、草本245种。按照植物类群划分,蕨类植物共
有4科5属6种、裸子植物共有2科3属3种、被子植物共有47科127属268种(其中单子叶植物共有9科
14属27种,双子叶植物共有38科113属241种)。以属与种数量为依据排序前十科为菊科(24属63种)、蔷
薇科(14属35种)、豆科(7属16种)、伞形科(8属15种)、毛茛科(4属14种)、桔梗科(3属10种)、牻牛儿
苗科(1属9种)、石竹科(4属9种)、堇菜科(1属7种)、天门冬科(4属7种)共70属185种,占调研样地属
种的51.85%与66.79%。以群落内植株数量为依据,灌木层群落以蔷薇科(50.61%)、桦木科(25.95%)、豆科
(15.17%)、杜鹃花科(6.64%)等为优势科,草本层以莎草科(32.44%)、蔷薇科(15.57%)、杜鹃花科(8.79%)、
天门冬科(7.09%)、菊科(6.05%)、鸢尾科(5.17%)等为优势科,不同海拔梯度常见科的种类与占比有所不
同。由于实际分类精确问题,参考丁茂等人
[39]
方法,将林下莎草科(Cyperaceae)薹草属(Carexsp.)植物统一
标记为薹草属。
3.2 林地优势物种组成及其重要值沿海拔梯度变化
随着海拔梯度的不断升高,乔灌草层内优势物种重要值发生明显变化(表2)。乔木层中落叶松、白桦呈
上升趋势,蒙古栎、黑桦与山杨呈下降趋势;灌木层中榛、胡枝子集中分布于E1—E3梯度,并随海拔升高呈下
降趋势。刺蔷薇、越橘、西伯利亚刺柏、杜香、偃松等随海拔呈上升趋势,西伯利亚刺柏、偃松等仅分布于E5—
http://www.ecologica.cn
7期 史恭发 等:内蒙古大兴安岭东部林下植物生物多样性海拔差异
3009
E6梯度。土庄绣线菊、珍珠梅、绣线菊、兴安杜鹃、笃斯越橘等呈先上升后下降趋势;草本层红花鹿蹄草、石生
单花鸢尾、草问荆、小玉竹、唐松草、褐苞蒿、柳兰、矮山黧豆等随海拔上升呈下降趋势。
表2 内蒙古大兴安岭林区东部不同海拔梯度林地优势物种组成及其重要值
Table2 ImportantvaluesofunderstoryplantsintheeasternpartoftheGreaterXing′anMountainsforestregion,InnerMongolia
林层
Forestlayer
乔木层Treelayer
物种
Species
落叶松Larixgmelinii
蒙古栎Quercusmongolica
白桦Betulaplatyphylla
黑桦Betuladahurica
山杨Populusdavidiana
灌木层Shrublayer榛Corylusheterophylla
刺蔷薇Rosaacicularis
胡枝子Lespedezabicolor
越橘Vacciniumvitis⁃idaea
土庄绣线菊Spiraeapubescens
珍珠梅Sorbariasorbifolia
绣线菊Spiraeasalicifolia
兴安杜鹃Rhododendrondauricum
西伯利亚刺柏Juniperuscommunis
杜香Rhododendrontomentosum
笃斯越橘Vacciniumuliginosum
偃松Pinuspumila
草本层Fieldlayer薹草属Carexsp.
蚊子草Filipendulapalmata
红花鹿蹄草Pyrolaasarifolia
东方草莓Fragariaorientalis
地榆Sanguisorbaofficinalis
单花鸢尾Irisuniflora
大叶章Deyeuxiapurpurea
铃兰Convallariakeiskei
草问荆Equisetumpratense
莓叶委陵菜Potentillafragarioides
小玉竹Polygonatumhumile
唐松草Thalictrumaquilegiifolium
褐苞蒿Artemisiaphaeolepis
北悬钩子Rubusarcticus
石生悬钩子Rubussaxatilis
柳叶蒿Artemisiaintegrifolia
路边青Geumaleppicum
柳兰Chamerionangustifolium
舞鹤草Maianthemumbifolium
矮山黧豆Lathyrushumilis
重要值Importancevalue/%
9.11
53.11
9.52
18.47
10.42
46.52
3.29
22.45
—
0.96
0.49
0.49
—
—
—
—
—
9.64
4.71
—
4.26
5.34
7.76
4.52
8.80
2.86
8.28
2.46
2.49
0.74
1.04
—
0.87
4.50
1.30
—
1.31
E1
6.18
56.31
12.75
16.63
6.51
45.00
3.36
35.81
1.11
1.42
4.41
0.52
—
—
—
—
—
11.39
8.72
6.17
6.44
6.09
9.03
4.24
7.25
4.99
3.62
3.11
2.59
0.78
2.27
1.47
0.70
1.72
1.35
0.72
2.05
E2
47.27
16.53
20.48
3.47
6.13
22.00
11.24
15.14
4.15
11.01
6.14
4.74
5.35
—
0.15
1.84
—
13.03
9.60
9.35
6.39
4.20
5.55
5.17
7.42
5.96
2.56
5.40
2.75
1.27
3.00
2.41
2.10
1.27
2.53
1.81
2.68
E3
64.49
2.71
28.47
1.03
6.19
8.55
11.07
1.67
6.00
15.63
8.37
10.44
7.40
—
3.40
2.11
—
12.55
18.87
12.99
6.73
5.29
4.94
4.75
5.21
5.18
1.20
3.51
1.49
1.77
3.42
3.25
0.42
1.18
2.73
1.90
1.25
E4
61.84
1.35
29.98
0.27
2.14
6.35
9.75
—
11.88
6.82
5.96
3.58
15.92
2.81
6.05
4.15
3.45
12.39
17.70
15.95
6.30
5.27
2.73
8.94
4.52
3.42
1.11
2.06
1.12
0.86
1.91
4.10
0.19
0.14
1.88
2.78
0.40
E5
76.03
0.83
23.97
—
—
—
13.75
—
29.36
—
—
4.51
9.13
10.82
6.12
1.23
19.15
11.93
5.15
19.19
5.96
3.36
0.83
7.38
3.75
3.26
—
—
0.80
—
0.82
—
0.29
—
1.45
2.91
—
E6
悬钩子、舞鹤草等随海拔升高呈上升趋势,铃兰、莓叶委陵菜、路边青等随海拔呈下降趋势,而薹草属、蚊子草、
E1:200—350m,E2:350—500m,E3:500—650m,E4:650—800m,E5:800—950m,E6:950—1100m;乔木层列举重要值排名前5名物种,灌
木层列举重要值排名前12名物种,草本层列举重要值排名前20名物种
林下植物多样性与林地上层植被类型组成密切相关
[41]
,根据物种重要值差异将同一海拔梯度内乔灌群
http://www.ecologica.cn
3010
生 态 学 报 44卷
落类型进行划分
[36]
,在200—1100m海拔高度内,以150m为1个海拔梯度,海拔由低至高乔灌群落依次为:
蒙古栎⁃黑桦⁃榛⁃胡枝子群落(E1—E2)、落叶松⁃白桦⁃榛⁃胡枝子群落(E3)、落叶松⁃白桦⁃欧亚绣线菊⁃绣线菊
群落(E4)、落叶松⁃白桦⁃兴安杜鹃⁃越橘群落(E5)、落叶松⁃白桦⁃越橘⁃偃松群落(E6)。
3.3 林下植物α多样性沿海拔梯度变化
α多样性是反映群落丰富度与均匀度的综合指标。研究发现(图2),灌木层α多样性Patrick、Shannon⁃
Wiener与Simpson指数呈波动上升趋势,E4梯度Patrick与Simpson指数分别达到最高值4.62与0.79,
Shannon⁃Wiener指数则于E6梯度呈现最高值0.683。Pielou指数波动较大,随海拔梯度上升呈现“W”型趋
势;草本层α多样性随海拔梯度升高变化基本一致,呈现先下降后升高再下降的波动下降趋势,各α多样性
指数在E3海拔梯度均达到最高;同一海拔梯度草本层α多样性指数均高于灌木层,说明研究区林下草本生
物多样性大于灌木层,并且灌木层与草本层α多样性随海拔梯度的变化形式不同,灌木层总体呈现波动上升
趋势,草本层呈现先下降后升高再下降的波动下降趋势,相比草本,灌木对于较高海拔梯度环境适应性更强。
Shannon⁃Wiener指数能综合地描述物种多样性水平
[17]
,以Shannon⁃Wiener指数为评判标准,灌木层与草本层
物种多样性分别在E6与E3海拔梯度达到峰值。
图2 林下植物灌木层与草本层α多样性指数海拔变化特征
Fig.2 Shrubsandherbaceouslayersofunderstoryplantsαaltitudevariationcharacteristicsofdiversityindex
进一步对林下植物群落α多样性指数随海拔升高变化绘制散点图,并使用多项式生成相应的拟合曲线
(图3),拟合曲线均为二项式。调研区林下植物群落Patrick、Shannon⁃Wiener与Simpson指数均呈现较为明显
地先升高后降低的单峰形式曲线。Patrick指数拟合曲线,E1—E2梯度内上升,至E3梯度达到峰值,E4—E6
梯度内迅速下降(P<0.01),Shannon⁃Wiener指数拟合曲线,E1—E2梯度缓慢上升、E3—E6梯度缓慢下降(P<
Pielou指数拟合曲线呈现随海拔升高而降低的趋势(P>0.05)。
0.05)。Simpson指数拟合曲线在E1梯度内呈缓慢上升趋势,在E2—E6梯度呈现缓慢下降趋势(P<0.05),
http://www.ecologica.cn
7期 史恭发 等:内蒙古大兴安岭东部林下植物生物多样性海拔差异
3011
图3 林下植物α多样性指数与海拔的拟合方程
Fig.3 Understoryplantαfittingequationbetweendiversityandaltitude
3.4 林下植物β多样性沿海拔梯度变化
3.4.1 林下植物Cody指数(β
C
)随海拔梯度变化
灌木物种最丰富、E2梯度草本物种最丰富。灌木层与草本层相邻海拔梯度间Cody指数(β
C
)值均先上升后下
草本层在前两个相邻海拔梯度间达到44.5与49.5,随后逐渐下降,说明在E1与E2、E2与E3海拔梯度间,草
本层物种具有较高的替换速率,并在E2与E3梯度间达到峰值。灌木层与草本层均在物种种类最丰富的海
在E1—E6海拔梯度分别统计到16、12、22、25、15、19种灌木与135、166、151、116、65、28种草本,E4梯度
降(表3),灌木层E3—E4海拔梯度间达到相邻海拔梯度间最高值10,此区间段内灌木物种替换速率最高。
拔梯度间出现β
C
极值,此外,草本层β
C
值均远高于灌木层,说明草本层对于海拔升高产生的环境异质性较灌
3.4.2 林下植物相似度指数(C
S
)随海拔梯度变化
木层更加敏感,物种替换速率更高。
E1与E2、E3与E4、E3与E5、E4与E5海拔梯度间均呈现中等相似水平;在E1与E3、E2与E3、E1与E4、E2
E5、E3与E5、E4与E5、E5与E6海拔梯度间呈现低度相似水平,其余海拔梯度间呈现非常低的相似水平。
灌木层与草本层相邻海拔梯度间相似度指数C
S
均呈先下降后升高再下降的波动形式(表3)。灌木层在
与E4、E3与E6、E4与E6海拔梯度间均呈现低度相似水平,其余海拔梯度间呈非常低的相似水平。草本层在
E1与E2、E2与E3、E3与E4海拔梯度间呈现中度相似水平,E1与E3、E1与E4、E2与E4、E1与E5、E2与
3.5 林下植物多样性与地理因素及优势乔木种的关系
地理因素、优势林木组成与林下植物多样性相关性分析表明(表4),灌木层Patrick与Shannon⁃Wiener指
数与海拔呈显著正相关(P<0.05),相关性为0.191与0.225,与经度呈显著负相关(P<0.05),相关性为
-
0.161
与
-
0.163,与坡度呈显著负相关(P<0.05),相关性为
-
0.238与
-
0.155;灌木层Patrick与Simpson指数与纬度
呈显著正相关(P<0.05),相关性为0.368与0.180,Shannon⁃Wiener指数与纬度呈极显著正相关(P<0.01),相
http://www.ecologica.cn
3012
生 态 学 报 44卷
关性为0.261;灌木层α多样性与坡向呈负相关但未达到显著性水平(P>0.05)。草本层Patrick、Shannon⁃
Wiener与Pielou指数与海拔呈极显著负相关(P<0.01),相关性为
-
0.314、
-
0.244与
-
0.428;草本层Patrick与
Pielou指数与纬度呈显著正相关(P<0.05),相关性为0.204、0.253,Shannon⁃Wiener指数与纬度呈极显著正相
关(P<0.01),相关性为0.147;草本层Patrick与Shannon⁃Wiener指数与坡向呈显著正相关(P<0.05),相关性
为0.335与0.245。说明灌木层α多样性随海拔与纬度的增高而升高,随经度与坡度的增高而下降。草本层
α多样性随海拔升高而下降,随纬度增高与坡向朝阳而上升。海拔变化对草本层α多样性影响高于灌木层,
坡度的升高对灌木层α多样性产生明显影响,而坡向则对草本层α多样性产生明显影响。
Table3 Codyindexofunderstoryplantcommunity(β
C
)similarityindex(C
S
)
指数
Index
β
C
海拔
Altitude
E1
E2
E3
E4
E5
C
S
E6
E1
E2
E3
E4
E5
E6
1.000.71
1.00
0.49
0.47
1.00
灌木层Shrublayer
9.00
7.00
0.00
E3
9.50
10.00
0.00
9.50
E4E1E2E5E6E1E2E3
表3 林下植物Cody指数(β
C
)与相似度指数(C
S
)
草本层Herblayer
63.50
49.50
0.00
56.50
57.00
37.50
0.00
E4
57.00
60.50
53.00
37.50
0.00
0.43
0.48
0.51
0.59
1.00
E5
61.50
72.00
62.50
45.00
23.50
0.00
0.25
0.26
0.30
0.38
0.49
1.00
E6
0.004.00
0.00
10.00
11.00
5.00
5.50
0.00
0.38
0.35
0.65
0.74
1.00
11.50
12.50
7.50
9.00
5.50
0.00
0.33
0.38
0.52
0.47
0.59
1.00
0.0044.50
0.00
0.44
0.54
0.62
1.00
1.000.70
1.00
0.56
0.69
1.00
0.55
0.59
0.72
1.00
E1:200—350m,E2:350—500m,E3:500—650m,E4:650—800m,E5:800—950m,E6:950—1100m
表4 林下植物群落物种多样性与地理因素及优势乔木种Pearson相关性
Table4 Correlationbetweenspeciesdiversityofunderstoryplantcommunityandgeographicalfactorsandpearsondominanttreespecies
Patrick指数
Patrickindex
Shannon⁃Wiener指数
Shannon⁃Wienerindex
灌木层
Shrub
layer
草本层
Herb
layer
Simpson指数
Simpsonindex
Pielou指数
Simpsonindexindex
草本层
Herb
layer
项目
Project
地理因素
Geographicfactors
海拔/m
经度/(°)
纬度/(°)
坡度/(°)
坡向
优势林木组成
Dominantforestcomposition
郁闭度Canopydensity
落叶松L.gmelinii/%
蒙古栎Q.mongolica/%
白桦B.platyphylla/%
黑桦B.dahurica/%
山杨P.davidiana/%
∗:P<0.05;∗∗:P<0.01
0.191
∗
-
0.161
∗
0.368
∗
-
0.238
∗
-
0.010
灌木层
Shrub
layer
∗
-
0.314
∗
草本层
Herb
layer
-
0.127
0.108
0.204
∗
0.335
∗
-
0.163
∗
-
0.155
∗
-
0.002
0.183
∗
0.171
∗
0.027
-
0.148
-
0.141
0.225
∗
∗
0.261
∗
∗
-
0.244
∗
灌木层
Shrub
layer
0.096
0.077
0.057
草本层
Herb
layer
-
0.028
-
0.060
0.080
0.035
0.163
灌木层
Shrub
layer
0.049
0.033
-
0.136
0.133
∗
-
0.428
∗
∗
0.147
∗
0.180
∗
-
0.079
-
0.088
0.096
0.024
0.070
0.070
0.180
∗
-
0.130
-
0.074
-
0.036
0.012
-
0.052
0.094
0.253
∗
∗
-
0.413
∗
0.153
∗
-
0.229
-
0.015
∗
-
0.118
0.227
∗
-
0.362
∗
-
0.103
0.130
0.039
-
0.039
-
0.094
0.245
∗
-
0.155
∗
-
0.144
-
0.058
0.072
0.115
-
0.143
-
0.136
-
0.022
-
0.026
0.055
0.034
-
0.049
-
0.035
0.042
0.055
0.065
-
0.146
-
0.158
∗
0.108
0.096
-
0.074
-
0.054
-
0.014
-
0.187
∗
-
0.158
灌木层Patrick指数与郁闭度呈极显著正相关(P<0.01),相关性为0.153,Shannon⁃Wiener指数与郁闭度
呈显著正相关(P<0.05),相关性为0.183;草本层Patrick、Shannon⁃Wiener与Simpson指数则与郁闭度呈显著
负相关(P<0.05),相关性为
-
0.413、
-
0.362与
-
0.155。说明郁闭度的升高对于灌木层与草本层α多样性存在
相反的影响,灌木层α多样性随郁闭度升高略有升高,草本层α多样性随郁闭度升高而下降明显。林木组成
http://www.ecologica.cn
7期 史恭发 等:内蒙古大兴安岭东部林下植物生物多样性海拔差异
3013
同样对林下植物多样性产生影响,林间落叶松数量占比与灌木层Patrick指数呈极显著正相关(P<0.01),相
关性为0.227,与灌木层Shannon⁃Wiener和Simpson指数呈显著正相关(P<0.05),相关性为0.171与0.180;林
间黑桦数量占比与灌木层Patrick指数呈显著负相关(P<0.05),相关性为
-
0.187;林间山杨数量占比与草本层
Pielou指数呈显著负相关(P<0.05),相关性为
-
0.158。说明调研区落叶松增多可提升灌木层生物多样性,黑
桦增多会降低灌木层物种丰富度,而山杨的增多会降低草本层物种分布的均匀度。
4 讨论
调研区为内蒙古大兴安岭的东南坡,地势沿东南至西北方迅速升高,随纬度升高逐渐升高,为寒温带气候
区内典型的山岭地区,拥有广袤的森林地带。林下植物是森林生态系统的重要组成部分,研究寒温带山岭地
区林下植物生物多样性随海拔差异变化有助于了解寒温带山岭地区植物随海拔梯度分布规律,丰富生物多样
性海拔格局理论,掌握内蒙古大兴安岭地区林下植物资源现状。调研区林下植物共277种隶属于53科135
属,灌木32种、草本245种,与张喜亭等人
[42]
对大兴安岭多布库尔国家级自然保护区调研的269种植物相
似。调研区森林乔灌群落随海拔梯度不断变化,乔木层由低海拔梯度以蒙古栎、黑桦为主,逐渐变为以落叶
松、白桦为主,灌木层由以榛、胡枝子为主过渡到兴安杜鹃、越橘等为主,与邻近地区类似海拔梯度如鄂伦春旗
多布库尔保护区(400—600m)
[42]
、呼中森林自然保护区(800—1200m)
[43]
乔灌群落构成具有较高的相似性。
林下植物调查中发现了野大豆(Glycinesoja)、大花杓兰(Cypripediummacranthos)、蒙古黄芪(Astragalus
危、重点保护植物,使得了解与保护内蒙古大兴安岭东部林下植物生物多样性更有必要。
林下植物α多样性随海拔梯度呈现先上升后下降趋势,与常见的生物多样性随海拔升高的单峰格局相
符
[17,36,44]
。林下草本层为林下植物生物多样性的主要贡献层,其α多样性高于灌木层,这与岳永杰
[25]
、张殷
波等人
[45]
研究一致。林下草本层与灌木层α多样性的海拔差异性存在明显不同,随海拔上升,草本层α多样
性均先下降,在E2区间内达到谷值,至E3梯度到达最高峰后持续下降的变化趋势,与传统的单峰格局有所
不同,根据实际调研情况,E2梯度多位于莫旗与鄂旗的东南部,此地区人口相对较多,森林内部常出现大量农
田区域,一方面严重破坏了森林连续度,另一方面农田耕作造成明显的水土流失、农药污染,加之此地区的放
牧、采山等活动,推测其为E2梯度林下草本群落α多样性下降的重要原因之一。灌木层Patrick、Shannon⁃
Wiener与Simpson指数随海拔梯度升高呈逐渐上升趋势,而Pielou指数呈波动下降趋势,灌木层相比于草本
层更能适应高海拔梯度的生境。
灌木层Cody指数在E2与E3、E3与E4相邻两海拔梯度值逐渐升高,在E3与E4海拔区域内灌木物种替
换速率最高,灌木也在此区间内由植株较高,叶片较大的榛、胡枝子、珍珠梅等,转变为了植株较矮,叶片较小
的土庄绣线菊、绣线菊、刺蔷薇等植物,随着海拔进一度升高,越橘、杜香、笃斯越橘、偃松等植叶片更小的灌木
逐渐变为优势物种。草本层则具有更高的Cody指数,在E2与E3区间内物种替换速率最高。
大尺度区域内,植物群落多样性主要受到气候、地势的影响,小尺度区域易受到区域小气候、土壤等方面
的影响
[14]
。相比于单独山体的植物多样性海拔差异性研究
[15—18]
,山岭地区植物多样性不仅会受到小区域山
体坡向、坡度的影响,同时受到经度与纬度与海拔相结合的多重影响。林下灌木层多样性与海拔、纬度呈现正
相关、与经度呈现负相关,根据相关性,这三种因素对于灌木层生物多样性影响排序为纬度>海拔>经度,调研
区随着纬度的上升海拔逐渐升高。草本层与海拔呈负相关,与纬度呈正相关,结合调研区地势随经纬度变化
的情况,这种相关性符合前面的结果,值得注意的是,纬度同时与灌木层、草本层均呈正相关,推测一方面这是
由于调研区东南部低纬度地区人口密集,农耕放牧活动对森林生态造成的破坏导致林下植物生物多样性的下
降,与另一方面调研区西北部纬度较高地区相较于东南部具有更复杂的地形起伏更高的环境异质性,并且良
好的封山育林环境促使林下生物多样性上升导致的。从行政区域表现则为莫旗林下植物多样性低于鄂旗。
坡度方面,平缓地区灌木多样性较高。坡向方面,草本层多样性阴坡小于阳坡,而对灌木层多样性影响不明
http://www.ecologica.cn
membranaceusvar.mongholicus)、紫点杓兰(Cypripediumguttatum)、樱草(Primulasieboldii)等多种国家珍稀濒
3014
生 态 学 报 44卷
显。郁闭度的增加一定程度上能提高灌木层多样性,却显著降低草本层多样性,鉴于草本层对于林下植物多
样性的突出贡献,建议定期对过密林分间伐以获取较好的林下植物多样性,间伐强度与频率有待于进一步研
究
[46]
。此外,不同的乔木组成能够提供不同的林下生存环境,从而产生不同的林下植物群落结构,调研区内
落叶松的加入能够一定程度上增加灌木层多样性,并有降低草本层多样性的趋势,这可能是相对于阔叶林,落
叶松林下凋落物分解速率较慢,养分补充缓慢,表层土壤组成不同导致的。
本研究主要针对海拔这一因素,探讨内蒙古大兴安岭地区林下植物多样性海拔差异性变化。林下植物多
样性受到多种环境异质性影响,调研时还发现向阳山坡通常具有较多种类的林下植物,林缘林下植物多样性
普遍较高,较潮湿或离水源较近的林地物种乔木与灌木丰富度往往较高等规律,不同的土壤、不同的林型
[47]
、
不同坡向与坡度
[48]
,不同季节、不同林龄、不同育林手段对林下植物多样性均可能产生影响。在今后的研究
应重视不同因素综合对林下植物的影响,并增加不同调研季节,考虑气候因素,全面分析,对于掌握林下植物
多样性的多方面规律具有重要意义。
5 结论
内蒙古大兴安岭东部林下植物共277种隶属于53科135属,灌木32种、草本245种;随海拔梯度升高,
森林乔灌群落依次为蒙古栎⁃黑桦⁃榛⁃胡枝子群落(E1—E2)、落叶松⁃白桦⁃榛⁃胡枝子群落(E3)、落叶松⁃白桦⁃
欧亚绣线菊⁃绣线菊群落(E4)、落叶松⁃白桦⁃兴安杜鹃⁃越橘群落(E5)、落叶松⁃白桦⁃越橘⁃偃松群落(E6);林
下植物物种多样性随海拔上升呈现明显的单峰格局。草本层物种多样性普遍高于灌木层,灌木层生物多样性
随海拔升高呈上升趋势,草本层呈现先下降后升高再下降的波动下降趋势;E3—E4与E2—E3梯度内灌木层
与草本层物种替换速率分别达到峰值,草本层物种替换速率远高于灌木层;灌木层与草本层相邻海拔梯度间
相似度指数均呈先下降后升高再下降的波动形式;地理因素及优势乔木种不同能够显著影响林下植物α多样
性,灌木层与海拔、纬度呈正相关,与经度、坡度呈负相关,草本层与海拔呈负相关,与纬度、坡向呈正相关,海拔
对草本层α多样性影响高于灌木层。郁闭度与灌木层α多样性呈正相关,与草本层α多样性呈负相关;落叶松
增多可提升灌木层生物多样性,黑桦增多会降低灌木层物种丰富度,而山杨的增多会降低草本层物种分布的均
匀度。
参考文献(References):
[1] 牟长城,倪志英,李东,孙晓新,陈加利.长白山溪流河岸带森林木本植物多样性沿海拔梯度分布规律.应用生态学报,2007,18(5):
[2] MåsvikenJ,DalerumF,CousinsSAO.ContrastingaltitudinalvariationofalpineplantcommunitiesalongtheSwedishMountains.Ecologyand
[3] OldfatherMF,BrittonMN,PapperPD,KoontzMJ,HalburMM,DodgeC,FlintAL,FlintLE,AckerlyDD.Effectsoftopoclimatic
[4] 刘华训.我国山地植被的垂直分布规律.地理学报,1981,36(3):267⁃279.
7170⁃7180.
complexityonthecompositionofwoodyplantcommunities.AoBPlants,2016,8:plw049.
[5] 孙然好,陈利顶,张百平,傅伯杰.山地景观垂直分异研究进展.应用生态学报,2009,20(7):1617⁃1624.
Evolution,2020,10(11):4838⁃4853.
945⁃952.
[6] 杨朗生,刘兴良,刘世荣,潘红丽,冯秋红,李迈和.卧龙巴郎山川滇高山栎群落植物生活型海拔梯度特征.生态学报,2017,37(21):
[7] 金慧,赵莹,赵伟,尹航,秦立武,崔凯峰,李晓杰,黄祥童.长白山牛皮杜鹃群落物种多样性的海拔梯度变化及相似性.生态学报,
[8] JeníkJ.H.WalterVegetationoftheearthinrelationtoclimateandtheeco⁃physiologicalconditions.FoliaGeobotanicaetPhytotaxonomica,1975,
[9] KörnerC.Whyarethereglobalgradientsinspeciesrichness?mountainsmightholdtheanswer.TrendsinEcology&Evolution,2000,15(12):
[10] UlfMolau,方静.低纬与高纬地区的山地生物多样性格局.AMBIO⁃人类环境杂志,2004,33(S1):25⁃29,57.
1⁃7.
513⁃514.
[11] 姚兰,艾训儒,吕世安,冯广,刘峻城,黄永涛.湖北星斗山天然次生林的群落类型、结构与物种多样性特征.林业科学,2015,51(11):
[12] 朱源,康慕谊,江源,刘全儒.贺兰山木本植物群落物种多样性的海拔格局.植物生态学报,2008,32(3):574⁃581.
[13] 梁红柱,刘丽丽,付同刚,高会,李敏,刘金铜.山地植被垂直分布变化格局研究进展与述评.中国生态农业学报:中英文,2022,30
10(1):100.
2015,35(1):125⁃133.
http://www.ecologica.cn
7期
(7):1077⁃1090.
史恭发 等:内蒙古大兴安岭东部林下植物生物多样性海拔差异
3015
[14] 刘星宏,张青青,徐海量,张鹏,张广鹏,李宏.北疆盐碱地植物群落空间分布及物种多样性.生态学报,2021,41(4):1501⁃1513.
[16] 冯相艳,赵文智,蔺鹏飞,王川.祁连山北坡主要木本植物功能性状及其海拔分异.生态学报,2022,42(23):9726⁃9735.
报,2021,41(3):1148⁃1157.
[15] 陈绪辉,林碧华,罗敏贤,刘宝,刘金福,郑世群.福建灵源山森林群落林下灌草生态位特征.森林与环境学报,2023,43(2):152⁃159.
[17] 李梦佳,何中声,江蓝,谷新光,晋梦然,陈博,刘金福.戴云山物种多样性与系统发育多样性海拔梯度分布格局及驱动因子.生态学
[18] 向琳,陈芳清,官守鹏,王玉兵,吕坤.井冈山鹿角杜鹃群落灌木层功能多样性及其随海拔梯度的变化.生态学报,2019,39(21):
[19] 连承杰.内蒙古大兴安岭林区林下经济发展现状与对策.中国林业经济,2021(1):43⁃45.
76⁃83.
8144⁃8155.
[20] 王群,李东海,杨小波,曾润娟,夏丹,王豪,戚春林,郝杰威.海南霸王岭橡胶林林下植物多样性特征研究.林业资源管理,2022(5):
[21] 谭许脉,张文,肖纳,莫雪青,高冠女,尤业明,明安刚,黄雪蔓.杉木林改造成乡土阔叶林对林下植物物种组成和多样性的影响.生态
[22] 王纯华.大兴安岭林区乔木树种与林下伴生经济植物分布情况.现代农业科技,2020(13):120,122.
学报,2022,42(7):2931⁃2942.
[23] 王纪坤,兰国玉,陈帮乾,谢贵水.海南橡胶林林下药用植物资源调查.南方农业学报,2013,44(3):391⁃396.
[24] 陈树彪,郑野,李秀英.大兴安岭呼中林业局林下主要经济植物资源的调查.内蒙古林业调查设计,2011,34(5):27⁃28,124.
[26] 张建宇.大兴安岭森林植物多样性、群落结构特征及耦合关系分析[D].哈尔滨:东北林业大学,2018.
[29] 冶连华,吕金娃.内蒙古大兴安岭林区野生宿根花卉资源调查.河北林业科技,2013(2):25⁃29.
37,134.
版,2020,41(4):21⁃26.
[25] 岳永杰,苏志成,韩君,格日乐高娃,乌云珠拉.内蒙古大兴安岭林区典型森林群落物种多样性研究.内蒙古农业大学学报:自然科学
[27] 崔佳佳,铁牛.大兴安岭北部森林群落结构及植物多样性特征研究.西北林学院学报,2021,36(2):24⁃30.
[28] 李桂茹.内蒙古大兴安岭林区野生经济植物资源评价.内蒙古林业调查设计,2011,34(1):98⁃99,120.
[30] 王翠敏.内蒙古大兴安岭林区兴安落叶松林下分布主要经济植物的经济价值及开发与利用.内蒙古林业调查设计,2014,37(2):35⁃
[31] 魏麒苏.内蒙古大兴安岭北部原始林区植物资源评价.内蒙古林业调查设计,2020,43(2):32,63⁃64.
41⁃44.
[32] 张重岭,周本志,兰学君,杜德刚.内蒙古大兴安岭北部原始林区兴安落叶松林植物多样性的调查.内蒙古林业调查设计,2003,26(1):
[33] 王丽红,辛颖,邹梦玲,赵雨森,刘双江,任清胜.大兴安岭火烧迹地植被恢复中植物多样性与生物量分配格局.北京林业大学学报,
[34] 王纯华.大兴安岭野生经济植物资源调查方法.现代农业科技,2020(10):122,125.
30⁃35.
2015,37(12):41⁃47.
[35] 丛静,尹华群,卢慧,宿秀江,杨敬元,李迪强,张于光.神农架保护区典型植被的物种多样性和环境解释.林业科学,2013,49(05):
[36] 董雪,李永华,辛智鸣,段瑞兵,姚斌,包岩峰,张正国,刘源.河西走廊西段荒漠戈壁灌木群落物种多样性的海拔格局.林业科学,
[37] 傅沛云.东北植物检索表.第2版.北京:科学出版社,1995:1⁃906.
自然保护区生物多样性研究.北京:科学出版社,2015:1⁃150.
2021,57(2):168⁃178.
[38] 赵敏,郑宝江,孙海,吕连宽,赵宝军,胡建民,周明,龚文峰,许青,崔岱宗,吴玉环,任强,陶雷,邵景文,黄世臣.额尔古纳国家级
[39] 丁茂,郑标,王丽,张金铭,周守标,王智.鹞落坪国家级自然保护区草本植物多样性及空间分布.生态学报,2022,42(19):8015⁃8030.
[40] 黄甫昭,丁涛,李先琨,郭屹立,王斌,向悟生,文淑均,李冬兴,何运林.弄岗喀斯特季节性雨林不同群丛物种多样性随海拔的变化.
[41] 张涵丹,康希睿,邵文豪,杨旭,张建锋,刘学全,陈光才.不同类型杉木人工林林下草本植物多样性特征.生态学报,2021,41(6):
[42] 张喜亭,张建宇,肖路,陈胜仙,仲召亮,高薇,王文杰.大兴安岭多布库尔国家级自然保护区植物多样性和群落结构特征.生态学报,
[43] 孙玉成,张喜亭,季倩如,张建宇,王文杰.大兴安岭泰加林植物多样性特征调查与分析.林业资源管理,2022(6):124⁃130.
化植物群落多样性和地上生物量沿海拔梯度的变化特征.生态学报,2022,42(9):3640⁃3655.
2022,42(1):176⁃185.
[44] 王采娥,黄梅,王文银,李子好,张涛,马林,白彦福,王彦龙,施建军,龙瑞军,刘玉,王晓丽,马玉寿,尚占环.三江源区高寒坡地退
[45] 张殷波,秦浩,孟庆欣,张峰,唐志尧.太行山森林群落物种多样性空间格局及其影响因素.应用与环境生物学报,2022,28(2):
[46] 龚固堂,牛牧,慕长龙,陈俊华,黎燕琼,朱志芳,郑绍伟.间伐强度对柏木人工林生长及林下植物的影响.林业科学,2015,51(4):
[47] 杨光辉,秦树林,金光泽.小兴安岭五种林型早春草本植物多样性及其环境解释.生态学报,2023,43(3):1234⁃1246.
8⁃15.
[48] 崔博亮,安彦明,张起鹏,蒲萌,陈克龙.坡度对亚高寒草甸狼毒群落植物多样性的影响.草业科学,2021,38(5):823⁃834.
331⁃338.
2118⁃2128.
生态学报,2016,36(14):4509⁃4517.
http://www.ecologica.cn
2024年6月14日发(作者:慕容翠柏)
第44卷第7期
2024年4月
ACTAECOLOGICASINICA
生态学报
Vol.44,No.7
Apr.,2024
DOI:10.20103/j.stxb.202305060941
史恭发,徐诺,牛钊倩,孙炜伦,王好,石广振,王玲.内蒙古大兴安岭东部林下植物生物多样性海拔差异.生态学报,2024,44(7):3004⁃3015.
Mountains,InnerMongolia.ActaEcologicaSinica,2024,44(7):3004⁃3015.
ShiGF,XuN,NiuZQ,SunWL,WangH,ShiGZ,WangL.AltitudinaldifferencesofunderstoryplantbiodiversityineasternGreaterXing′an
内蒙古大兴安岭东部林下植物生物多样性海拔差异
史恭发,徐 诺,牛钊倩,孙炜伦,王 好,石广振,王 玲
∗
东北林业大学园林学院,哈尔滨 150040
摘要:海拔影响生物多样性分异,为探究寒温带山岭地区海拔对林下植物多样性分布影响,选择内蒙古大兴安岭东部地区,采用
样地调研法,设置E1(200—350m),E2(350—500m),E3(500—650m),E4(650—800m),E5(800—950m),E6(950—1100m)
6个不同海拔梯度,调研165个林地样点以了解不同海拔森林群落物种组成,研究林下植物生物多样性海拔差异性、地形因素、
乔木层优势物种与林下植物生物多样性的相关性。研究得出:(1)调研区林下植物共277种隶属于53科135属,灌木32种、草
白桦⁃欧亚绣线菊⁃绣线菊群落(E4)、落叶松⁃白桦⁃兴安杜鹃⁃越橘群落(E5)、落叶松⁃白桦⁃越橘⁃偃松群落(E6);(3)林下植物群
落生物多样性随海拔上升呈现明显的单峰格局。生物多样性草本层>灌木层,灌木层随海拔升高呈上升趋势,草本层呈先下降
后升高再下降的波动下降趋势;(4)草本层物种替换速率远高于灌木层,E3—E4与E2—E3梯度内灌木层与草本层物种替换速
率分别达到峰值;灌木层与草本层相邻海拔梯度间相似度指数均呈先下降后升高再下降的波动形式;(5)地理因素及优势乔木
种不同能够显著影响林下植物α多样性,灌木层α多样性与海拔、纬度呈正相关,与经度、坡度呈负相关。草本层与海拔呈负
相关,与纬度、坡向呈正相关,海拔对草本层α多样性影响高于灌木层。郁闭度与灌木层α多样性呈正相关,与草本层α多样
性呈负相关;落叶松增多可提升灌木层生物多样性,黑桦会降低灌木层物种丰富度,而山杨会降低草本层物种分布的均匀度。
研究结果为保护与利用内蒙古大兴安岭林下植物,丰富寒温带山岭地区林下植物多样性海拔分布理论提供科学依据。
关键词:内蒙古大兴安岭;林下植物群落;生物多样性;海拔
本245种;(2)森林乔灌群落随海拔依次为蒙古栎⁃黑桦⁃榛⁃胡枝子群落(E1—E2)、落叶松⁃白桦⁃榛⁃胡枝子群落(E3)、落叶松⁃
AltitudinaldifferencesofunderstoryplantbiodiversityineasternGreaterXing′an
Mountains,InnerMongolia
SHIGongfa,XUNuo,NIUZhaoqian,SUNWeilun,WANGHao,SHIGuangzhen,WANGLing
∗
CollegeofLandscapeArchitecture,NortheastForestryUniversity,Harbin150040,China
Abstract:Inordertoexploretheeffectofaltitudeonthediversityofundergrowthplantsinthecoldtemperatemountain
area,weselectedtheeasternGreaterXing′anMountainsinInnerMongoliabyusingtheplotinvestigationmethod,andset
(800—950m),andE6(950—1100m).Weinvestigated165forestsitestounderstandthecompositionofforest
communitiesatdifferentaltitudes,andstudiedtherelationshipbetweenthealtitudedifferenceinunderstoryplant
sixdifferentaltitudegradients:E1(200—350m),E2(350—500m),E3(500—650m),E4(650—800m),E5
biodiversity,topographicfactors,anddominantspeciesinthearboreallayerandunderstoryplantbiodiversity.Theresults
245herbs.(2)TheforesttreeandshrubcommunitieswithdifferentaltitudeswereintheorderofQuercusMongolia,Betula
heterophylla,Lespedezabicolorcommunities(E3),Larixgmelina,Betulaplatyphylline,Spiraeamedia,Spiraeasilicified
基金项目:国家科技基础资源调查专项(2019FY100500);中央高校基本科研业务费专项资金(2572023CT18)
收稿日期:2023⁃05⁃06; 网络出版日期:2024⁃01⁃12
∗通讯作者Correspondingauthor.E⁃mail:wanglinghlj@126.com
showedthat:(1)therewere277understoryplantspeciesbelongingto53familiesand135genera,including32shrubsand
daturic,Corylusheterophylla,Lespedezabicolorcommunities(E1—E2),Larixgmelina,Betulaplatyphylline,Corylus
http://www.ecologica.cn
7期 史恭发 等:内蒙古大兴安岭东部林下植物生物多样性海拔差异
3005
communities(E4),Larixgmelina,Betulaplatyphylla,Rhododendrondauricum,Vacciniumvitis⁃idaeacommunities(E5),
Larixgmelinii,Betulaplatyphylla,Vacciniumvitis⁃idaea,Pinuspumilacommunities(E6).(3)Thebiodiversityof
understoryplantcommunitiesshowedanobviouslyunimodalpatternwiththeincreasesofaltitude.Theherblayerwashigher
fluctuatingdecreasetrendoffirstdecreasingandthenincreasinganddecreasing.(4)Thereplacementrateofspeciesin
thantheshrublayer.Theshrublayershowedanincreasingtrendwiththeincreasingaltitude,whiletheherblayershoweda
herbaceouslayerwasmuchhigherthanthatinshrublayer,withthereplacementrateinshrubandherbaceouslayerpeaking
thenincreased,andfinallydecreased.(5)Thegeographicalfactorsanddifferentdominanttreespeciescouldsignificantly
andnegativelycorrelatedwithlongitudeandslope.Theherbaceouslayerwasnegativelycorrelatedwithaltitude,but
inE3—E4andE2—E3gradients,respectively.Thesimilarityindexbetweenshrublayerandherbagelayerdecreasedfirst,
affecttheαdiversityofunderstoryplants.Theαdiversityofshrublayerwaspositivelycorrelatedwithaltitudeandlatitude,
positivelycorrelatedwithlatitudeandslopeaspect.Theeffectofaltitudeonαdiversityofherbaceouslayerwashigherthan
thatofshrublayer.Thecanopydensitywaspositivelycorrelatedwiththeαdiversityofshrublayerandnegativelycorrelated
withtheαdiversityofherblayer.Larixgmeliniiincreasedbiodiversityinshrublayer,Betuladahuricadecreasedspecies
richnessinshrublayer,andBetuladahuricadecreasedspeciesdistributionevennessinherbaceouslayer.Ourresults
providedscientificbasisfortheprotectionandutilizationofundergrowthplantsintheGreaterXing′anMountainsandthe
theoryofaltitudinaldistributionofunderstoryplantdiversityinthecoldtemperatemountains.
KeyWords:GreaterXing′anMountainregionofInnerMongolia;understoryplantcommunity;bio⁃diversity;altitude
pattern
降水等多方面环境因子的复合变化,对群落分布与物种多样性产生重要影响
[3]
。山地是具有一定海拔和坡
度变化的地面
[4]
,其具有浓缩的环境梯度、高度异质化的生境
[5]
,通常具有明显的生物多样性与群落的垂直
地带分异
[6]
,始终是生态学研究热点
[7]
,山地生物的垂直变化是水平梯度的近千倍
[8]
,山地垂直带被认为是
水平地带的“微缩模型”
[9]
。植被的垂直带性虽在形态上与水平地带性具有一定的相似性,但绝非等同水平
地带的重复
[4,10]
,山地生态系统植物多样性随海拔上升一般有5种形式,先升后降的单峰格局
[11]
与单调下降
的格局较为普遍
[12]
,此外还存在单调升高、先降低后升高、和不具相关性3种格局
[12]
。根据温度和水分差
异,我国山地垂直自然带划分为大陆性和季风性两类带谱体系
[13]
,不同地区、不同生活型的物种多样性海拔
分布格局具有差异性,研究的区域不同得到的结果不同,因此对不同地区植物多样性海拔格局展开研究十分
必要
[12]
。
不同尺度区域下,植物群落分布的影响因子不同,小尺度区域植物群落主要由地形、土壤等方面影响,大
生物多样性的海拔差异是生态学研究的重要问题
[1—2]
,海拔作为综合性驱动因子,能够引起光照、温度、
尺度区域主要由气候、地势等方面决定
[14]
。许多单一山体按一定方向呈带状延伸形成山岭,相比于单个山
体,山岭地区具有经纬度跨度更大的水平地带性与干湿度地带性,具有海拔变化梯度更大的垂直地带性,因此
具有更丰富的“三维空间”分布现象,此外还具有人类干扰程度较小,对气候环境反应灵敏等特点,彼此之间
相互交叉过渡,更容易形成生物多样性的聚集地。对于山岭地区植被随海拔梯度的分布规律相对于单一山体
更能得出山地植被垂直分布的一般规律性。我国植物多样性海拔梯度的研究集中于单独山体
[15—18]
,并且近
一个世纪山地植被地带性垂直分布变化格局的相关研究集中于我国横断山区、青藏高原、太白山及秦岭区
域
[13]
,十分缺乏东北山岭地区植被地带性垂直分布格局的相关研究。
大兴安岭位于大陆性气候和温带季风气候区分界线与我国地势二三级阶梯分界线,具有较大的海拔差
异。内蒙古大兴安岭林区占整个大兴安岭的46%,森林面积8.37万km
2[19]
并以天然林为主,包括110万hm
2
未开发的原始森林,森林覆盖率78.39%,居国有林区之首,在东北地区具有丰富的植物资源。林下植物是森
林生态系统的重要组成部分,在维持生态系统稳定性、促进群落演替更新、保证人工林生物多样性与生态功能
http://www.ecologica.cn
3006
生 态 学 报 44卷
可持续发展等方面具有重要作用
[20—21]
。林下植物蕴含大量药用、食用、观赏等方面的经济植物类群
[22—24]
,在
推动当地林下产业经济发展方面具有重大潜力。了解林下植物的物种组成和生物多样性梯度分布是进一步
开展保护与开发内蒙古大兴安岭林下植物群落的先期工作。
关系
[26—27]
,林下经济植物资源分布与评价
[22,28—31]
,落叶松林下植物多样性
[32]
与火烧迹地植物多样性
[33]
等方
面。本文选取内蒙古大兴安岭东部林下植物多样性海拔分布差异为研究对象,通过200—1100m连续分布样
地调研方法减少因局部小气候等因素产生的影响,根据重要值、α多样性指数与拟合曲线分析、β多样性指
数,林下植物多样性与地理因素及优势乔木种的相关性等方法研究了内蒙古大兴安岭东部森林植被物种组成
与林下植物多样性沿海拔梯度的分布规律,将灌木层与草本层相比较,分析两种群落层次间差异,考虑乔木层
优势物种比重对林下植物多样性的影响,旨在深入了解内蒙古大兴安岭东部森林群落随海拔分布、林下灌木
与草本物种多样性的海拔分布、地理因子与乔木优势物种组成对林下植物多样性的影响等问题。丰富了我国
寒温带山岭森林植被垂直地带性分布的研究内容,填补了我国东北地区山岭植被垂直地带性分布的研究空
缺,为了解大陆性和季风性两类带谱过渡地区垂直自然带谱分布特征及其规律,保护与利用内蒙古大兴安岭
林下植物,开展不同海拔梯度人工林营造等工作提供科学依据。
1 研究区概况
研究区位于内蒙古大兴安岭林区东部莫力达瓦达斡尔族自治旗(莫旗)与鄂伦春族自治旗(鄂旗),为北
半球中高纬度地区(48°64′—51°13′N,122°05′—125°14′E),位于大兴安岭东南坡。调研区属寒温带大陆性
季风气候,冬寒夏暖,昼夜温差大,1月份平均气温
-
24—
-
31℃,7月份平均气温16—20℃,平均年降水量
300—500mm,河流较多,以嫩江水系为主;地势沿东南至西北方向海拔迅速升高,坡度逐渐增加,海拔164—
1484m之间,具有较大落差,以山地、丘陵、平原等地貌为主。
2 研究方法
2.1 样地设置与调研方法
调研使用典型取样与随机取样相结合的样地调研方法,在内蒙古大兴安岭东部莫旗及鄂旗两地沿海拔
200—1100m范围分布的森林设置6个不同海拔梯度样带:E1:200—350m,E2:350—500m,E3:500—650m,
E4:650—800m,E5:800—950m,E6:950—1100m。根据实际可达性,选取165个森林样点开展实地调查
(图1)。
调查方法以王纯华等人
[26,34—36]
为参考,根据本调研目的做出调整。每个样点在半径500m范围内设置
目前对于内蒙古大兴安岭植物研究多集中于典型森林群落多样性的纬度格局
[25]
,多样性与群落结构间
3个标准地(即3个重复),尽可能涵盖样点区域不同林地位置与方向以获取全面的物种分布信息,每个标准
地由1个20m
×
20m乔木样方及其东北角1个5m
×
5m灌木样方、4个分布于4角的1m
×
1m草本样方组
成,总调研495个标准地,调研面积198000m
2
。记录每个样点的经纬度、海拔、坡度、坡向、郁闭度、林型等信
息(表1)。坡向数据以丛静等人
[35]
方法进行赋值,数字越大越向阳干热。乔木、灌木与草本样方中记录物种
名、株数等;对各调研样点生境、林下植物均拍照记录。现场鉴定参考第二版《东北植物检索表》
[37]
,借鉴《额
尔古纳自然保护区生物多样性》植物名录辅助鉴定
[38]
,最终植物正名参考中国植物志(网址:http://www.
2.2 数据处理
iplant.cn/frps)。调研时间为2022年6月—2022年9月。
运用AdobeLightroomClassic管理样地植物照片,辅助植物鉴定与统计;运用OfficeExcel(2019)汇总数据
并计算α与β多样性;使用IBMSPSSStatistics26进行单因素方差分析、Pearson相关性分析与线性回归分析;
运用Origin(2021)生成α多样性与海拔散点分布图并使用SimpleFit插件生成α多样性与海拔的拟合方程。
运用最小显著差异法检验群落相应条件下多样性指数在置信区间(P<0.05)水平上的差异显著性。
http://www.ecologica.cn
7期 史恭发 等:内蒙古大兴安岭东部林下植物生物多样性海拔差异
3007
图1 调研样点
Fig.1 Surveysamplepoint
表1 研究区6个海拔梯度调研样地特征
Table1 Characteristicsofsixaltitudegradientsurveysamplesitesinthestudyarea
海拔等级
Altitudelevel
E1
E2
E3
E4
E5
E6
海拔梯度
Altitudegradient/m
200—350
350—500
500—650
650—800
800—950
950—1100
样点数量
Numberofsample
points
30
35
34
27
28
11
平均乔木高
Averagetree
height/m
9.24
9.68
平均郁闭度
Averagecanopy
density
0.54
0.62
0.67
0.70
0.66
0.50
平均灌木高
Averageshrub
height/m
0.95
0.99
1.08
0.94
0.80
0.27
12.21
14.99
13.87
9.48
2.2.1 重要值
[39]
V
=
(相对盖度
+
相对频度
+
相对密度)/3
2.2.2 α多样性
(1)物种丰富度指数(Patrick指数)(R)
[20,40]
,指数值越高说明物种丰富度越高。
R
=
S
http://www.ecologica.cn
3008
生 态 学 报 44卷
(2)物种多样性指数Shannon⁃Wiener(H′),指数值越高说明物种多样性越高
[33,40]
。
H′
=-
(3)辛普森指数(优势度指数)Simpson指数值(D),指假设在无限个群落随机抽取样本,样本中两个不同
P
i
lnP
i
∑
i
=
1
S
0—1之间。
物种个体相遇的概率为多样性测度
[20,26]
,指数值越高说明物种多样性越高,而物种优势度越低。数值范围在
D
=
1
-
P
i
2
∑
i
=
1
S
内分布越均匀。数值范围在0—1之间。
(4)均匀度指数Pielou指数(J),表示与物种丰富度无关的均匀度指数
[7,40]
,指数值越高说明物种在样方
S表示研究系统中物种总数,P
i
表示第i个物种的相对密度(P
i
=
N
i
/N,N为所存在群落全部种的个数之
J
=
H′/lnS
2.2.3 β多样性
和,N
i
代表第i个种的个体数)。
(1)Cody指数β
C
[1]
β
C
=
g(H)为沿海拔梯度增加的物种数目;l(H)为沿海拔梯度失去的物种数目,β
C
数值越高表示海拔梯度间物
种替换速率越快。
(2)群落相似度指数C
S
[1]
C
S
=
g
(
H
)
+
l(H)
2
2c
a
+
b
c为两群落共有的物种数,a和b分别为两群落各自物种数,C
S
数值越高表示海拔梯度间物种相似度越
高。C
S
值在0.8—1.0之间为高度相似,0.6—0.8为中度相似,0.4—0.6为低度相似,0.2—0.4为非常低的相
似,0—0.2表示海拔梯度间物种几乎没有相似之处。
3 结果分析
3.1 林下植物物种组成
调研区林下植物共277种隶属于53科135属,灌木32种、草本245种。按照植物类群划分,蕨类植物共
有4科5属6种、裸子植物共有2科3属3种、被子植物共有47科127属268种(其中单子叶植物共有9科
14属27种,双子叶植物共有38科113属241种)。以属与种数量为依据排序前十科为菊科(24属63种)、蔷
薇科(14属35种)、豆科(7属16种)、伞形科(8属15种)、毛茛科(4属14种)、桔梗科(3属10种)、牻牛儿
苗科(1属9种)、石竹科(4属9种)、堇菜科(1属7种)、天门冬科(4属7种)共70属185种,占调研样地属
种的51.85%与66.79%。以群落内植株数量为依据,灌木层群落以蔷薇科(50.61%)、桦木科(25.95%)、豆科
(15.17%)、杜鹃花科(6.64%)等为优势科,草本层以莎草科(32.44%)、蔷薇科(15.57%)、杜鹃花科(8.79%)、
天门冬科(7.09%)、菊科(6.05%)、鸢尾科(5.17%)等为优势科,不同海拔梯度常见科的种类与占比有所不
同。由于实际分类精确问题,参考丁茂等人
[39]
方法,将林下莎草科(Cyperaceae)薹草属(Carexsp.)植物统一
标记为薹草属。
3.2 林地优势物种组成及其重要值沿海拔梯度变化
随着海拔梯度的不断升高,乔灌草层内优势物种重要值发生明显变化(表2)。乔木层中落叶松、白桦呈
上升趋势,蒙古栎、黑桦与山杨呈下降趋势;灌木层中榛、胡枝子集中分布于E1—E3梯度,并随海拔升高呈下
降趋势。刺蔷薇、越橘、西伯利亚刺柏、杜香、偃松等随海拔呈上升趋势,西伯利亚刺柏、偃松等仅分布于E5—
http://www.ecologica.cn
7期 史恭发 等:内蒙古大兴安岭东部林下植物生物多样性海拔差异
3009
E6梯度。土庄绣线菊、珍珠梅、绣线菊、兴安杜鹃、笃斯越橘等呈先上升后下降趋势;草本层红花鹿蹄草、石生
单花鸢尾、草问荆、小玉竹、唐松草、褐苞蒿、柳兰、矮山黧豆等随海拔上升呈下降趋势。
表2 内蒙古大兴安岭林区东部不同海拔梯度林地优势物种组成及其重要值
Table2 ImportantvaluesofunderstoryplantsintheeasternpartoftheGreaterXing′anMountainsforestregion,InnerMongolia
林层
Forestlayer
乔木层Treelayer
物种
Species
落叶松Larixgmelinii
蒙古栎Quercusmongolica
白桦Betulaplatyphylla
黑桦Betuladahurica
山杨Populusdavidiana
灌木层Shrublayer榛Corylusheterophylla
刺蔷薇Rosaacicularis
胡枝子Lespedezabicolor
越橘Vacciniumvitis⁃idaea
土庄绣线菊Spiraeapubescens
珍珠梅Sorbariasorbifolia
绣线菊Spiraeasalicifolia
兴安杜鹃Rhododendrondauricum
西伯利亚刺柏Juniperuscommunis
杜香Rhododendrontomentosum
笃斯越橘Vacciniumuliginosum
偃松Pinuspumila
草本层Fieldlayer薹草属Carexsp.
蚊子草Filipendulapalmata
红花鹿蹄草Pyrolaasarifolia
东方草莓Fragariaorientalis
地榆Sanguisorbaofficinalis
单花鸢尾Irisuniflora
大叶章Deyeuxiapurpurea
铃兰Convallariakeiskei
草问荆Equisetumpratense
莓叶委陵菜Potentillafragarioides
小玉竹Polygonatumhumile
唐松草Thalictrumaquilegiifolium
褐苞蒿Artemisiaphaeolepis
北悬钩子Rubusarcticus
石生悬钩子Rubussaxatilis
柳叶蒿Artemisiaintegrifolia
路边青Geumaleppicum
柳兰Chamerionangustifolium
舞鹤草Maianthemumbifolium
矮山黧豆Lathyrushumilis
重要值Importancevalue/%
9.11
53.11
9.52
18.47
10.42
46.52
3.29
22.45
—
0.96
0.49
0.49
—
—
—
—
—
9.64
4.71
—
4.26
5.34
7.76
4.52
8.80
2.86
8.28
2.46
2.49
0.74
1.04
—
0.87
4.50
1.30
—
1.31
E1
6.18
56.31
12.75
16.63
6.51
45.00
3.36
35.81
1.11
1.42
4.41
0.52
—
—
—
—
—
11.39
8.72
6.17
6.44
6.09
9.03
4.24
7.25
4.99
3.62
3.11
2.59
0.78
2.27
1.47
0.70
1.72
1.35
0.72
2.05
E2
47.27
16.53
20.48
3.47
6.13
22.00
11.24
15.14
4.15
11.01
6.14
4.74
5.35
—
0.15
1.84
—
13.03
9.60
9.35
6.39
4.20
5.55
5.17
7.42
5.96
2.56
5.40
2.75
1.27
3.00
2.41
2.10
1.27
2.53
1.81
2.68
E3
64.49
2.71
28.47
1.03
6.19
8.55
11.07
1.67
6.00
15.63
8.37
10.44
7.40
—
3.40
2.11
—
12.55
18.87
12.99
6.73
5.29
4.94
4.75
5.21
5.18
1.20
3.51
1.49
1.77
3.42
3.25
0.42
1.18
2.73
1.90
1.25
E4
61.84
1.35
29.98
0.27
2.14
6.35
9.75
—
11.88
6.82
5.96
3.58
15.92
2.81
6.05
4.15
3.45
12.39
17.70
15.95
6.30
5.27
2.73
8.94
4.52
3.42
1.11
2.06
1.12
0.86
1.91
4.10
0.19
0.14
1.88
2.78
0.40
E5
76.03
0.83
23.97
—
—
—
13.75
—
29.36
—
—
4.51
9.13
10.82
6.12
1.23
19.15
11.93
5.15
19.19
5.96
3.36
0.83
7.38
3.75
3.26
—
—
0.80
—
0.82
—
0.29
—
1.45
2.91
—
E6
悬钩子、舞鹤草等随海拔升高呈上升趋势,铃兰、莓叶委陵菜、路边青等随海拔呈下降趋势,而薹草属、蚊子草、
E1:200—350m,E2:350—500m,E3:500—650m,E4:650—800m,E5:800—950m,E6:950—1100m;乔木层列举重要值排名前5名物种,灌
木层列举重要值排名前12名物种,草本层列举重要值排名前20名物种
林下植物多样性与林地上层植被类型组成密切相关
[41]
,根据物种重要值差异将同一海拔梯度内乔灌群
http://www.ecologica.cn
3010
生 态 学 报 44卷
落类型进行划分
[36]
,在200—1100m海拔高度内,以150m为1个海拔梯度,海拔由低至高乔灌群落依次为:
蒙古栎⁃黑桦⁃榛⁃胡枝子群落(E1—E2)、落叶松⁃白桦⁃榛⁃胡枝子群落(E3)、落叶松⁃白桦⁃欧亚绣线菊⁃绣线菊
群落(E4)、落叶松⁃白桦⁃兴安杜鹃⁃越橘群落(E5)、落叶松⁃白桦⁃越橘⁃偃松群落(E6)。
3.3 林下植物α多样性沿海拔梯度变化
α多样性是反映群落丰富度与均匀度的综合指标。研究发现(图2),灌木层α多样性Patrick、Shannon⁃
Wiener与Simpson指数呈波动上升趋势,E4梯度Patrick与Simpson指数分别达到最高值4.62与0.79,
Shannon⁃Wiener指数则于E6梯度呈现最高值0.683。Pielou指数波动较大,随海拔梯度上升呈现“W”型趋
势;草本层α多样性随海拔梯度升高变化基本一致,呈现先下降后升高再下降的波动下降趋势,各α多样性
指数在E3海拔梯度均达到最高;同一海拔梯度草本层α多样性指数均高于灌木层,说明研究区林下草本生
物多样性大于灌木层,并且灌木层与草本层α多样性随海拔梯度的变化形式不同,灌木层总体呈现波动上升
趋势,草本层呈现先下降后升高再下降的波动下降趋势,相比草本,灌木对于较高海拔梯度环境适应性更强。
Shannon⁃Wiener指数能综合地描述物种多样性水平
[17]
,以Shannon⁃Wiener指数为评判标准,灌木层与草本层
物种多样性分别在E6与E3海拔梯度达到峰值。
图2 林下植物灌木层与草本层α多样性指数海拔变化特征
Fig.2 Shrubsandherbaceouslayersofunderstoryplantsαaltitudevariationcharacteristicsofdiversityindex
进一步对林下植物群落α多样性指数随海拔升高变化绘制散点图,并使用多项式生成相应的拟合曲线
(图3),拟合曲线均为二项式。调研区林下植物群落Patrick、Shannon⁃Wiener与Simpson指数均呈现较为明显
地先升高后降低的单峰形式曲线。Patrick指数拟合曲线,E1—E2梯度内上升,至E3梯度达到峰值,E4—E6
梯度内迅速下降(P<0.01),Shannon⁃Wiener指数拟合曲线,E1—E2梯度缓慢上升、E3—E6梯度缓慢下降(P<
Pielou指数拟合曲线呈现随海拔升高而降低的趋势(P>0.05)。
0.05)。Simpson指数拟合曲线在E1梯度内呈缓慢上升趋势,在E2—E6梯度呈现缓慢下降趋势(P<0.05),
http://www.ecologica.cn
7期 史恭发 等:内蒙古大兴安岭东部林下植物生物多样性海拔差异
3011
图3 林下植物α多样性指数与海拔的拟合方程
Fig.3 Understoryplantαfittingequationbetweendiversityandaltitude
3.4 林下植物β多样性沿海拔梯度变化
3.4.1 林下植物Cody指数(β
C
)随海拔梯度变化
灌木物种最丰富、E2梯度草本物种最丰富。灌木层与草本层相邻海拔梯度间Cody指数(β
C
)值均先上升后下
草本层在前两个相邻海拔梯度间达到44.5与49.5,随后逐渐下降,说明在E1与E2、E2与E3海拔梯度间,草
本层物种具有较高的替换速率,并在E2与E3梯度间达到峰值。灌木层与草本层均在物种种类最丰富的海
在E1—E6海拔梯度分别统计到16、12、22、25、15、19种灌木与135、166、151、116、65、28种草本,E4梯度
降(表3),灌木层E3—E4海拔梯度间达到相邻海拔梯度间最高值10,此区间段内灌木物种替换速率最高。
拔梯度间出现β
C
极值,此外,草本层β
C
值均远高于灌木层,说明草本层对于海拔升高产生的环境异质性较灌
3.4.2 林下植物相似度指数(C
S
)随海拔梯度变化
木层更加敏感,物种替换速率更高。
E1与E2、E3与E4、E3与E5、E4与E5海拔梯度间均呈现中等相似水平;在E1与E3、E2与E3、E1与E4、E2
E5、E3与E5、E4与E5、E5与E6海拔梯度间呈现低度相似水平,其余海拔梯度间呈现非常低的相似水平。
灌木层与草本层相邻海拔梯度间相似度指数C
S
均呈先下降后升高再下降的波动形式(表3)。灌木层在
与E4、E3与E6、E4与E6海拔梯度间均呈现低度相似水平,其余海拔梯度间呈非常低的相似水平。草本层在
E1与E2、E2与E3、E3与E4海拔梯度间呈现中度相似水平,E1与E3、E1与E4、E2与E4、E1与E5、E2与
3.5 林下植物多样性与地理因素及优势乔木种的关系
地理因素、优势林木组成与林下植物多样性相关性分析表明(表4),灌木层Patrick与Shannon⁃Wiener指
数与海拔呈显著正相关(P<0.05),相关性为0.191与0.225,与经度呈显著负相关(P<0.05),相关性为
-
0.161
与
-
0.163,与坡度呈显著负相关(P<0.05),相关性为
-
0.238与
-
0.155;灌木层Patrick与Simpson指数与纬度
呈显著正相关(P<0.05),相关性为0.368与0.180,Shannon⁃Wiener指数与纬度呈极显著正相关(P<0.01),相
http://www.ecologica.cn
3012
生 态 学 报 44卷
关性为0.261;灌木层α多样性与坡向呈负相关但未达到显著性水平(P>0.05)。草本层Patrick、Shannon⁃
Wiener与Pielou指数与海拔呈极显著负相关(P<0.01),相关性为
-
0.314、
-
0.244与
-
0.428;草本层Patrick与
Pielou指数与纬度呈显著正相关(P<0.05),相关性为0.204、0.253,Shannon⁃Wiener指数与纬度呈极显著正相
关(P<0.01),相关性为0.147;草本层Patrick与Shannon⁃Wiener指数与坡向呈显著正相关(P<0.05),相关性
为0.335与0.245。说明灌木层α多样性随海拔与纬度的增高而升高,随经度与坡度的增高而下降。草本层
α多样性随海拔升高而下降,随纬度增高与坡向朝阳而上升。海拔变化对草本层α多样性影响高于灌木层,
坡度的升高对灌木层α多样性产生明显影响,而坡向则对草本层α多样性产生明显影响。
Table3 Codyindexofunderstoryplantcommunity(β
C
)similarityindex(C
S
)
指数
Index
β
C
海拔
Altitude
E1
E2
E3
E4
E5
C
S
E6
E1
E2
E3
E4
E5
E6
1.000.71
1.00
0.49
0.47
1.00
灌木层Shrublayer
9.00
7.00
0.00
E3
9.50
10.00
0.00
9.50
E4E1E2E5E6E1E2E3
表3 林下植物Cody指数(β
C
)与相似度指数(C
S
)
草本层Herblayer
63.50
49.50
0.00
56.50
57.00
37.50
0.00
E4
57.00
60.50
53.00
37.50
0.00
0.43
0.48
0.51
0.59
1.00
E5
61.50
72.00
62.50
45.00
23.50
0.00
0.25
0.26
0.30
0.38
0.49
1.00
E6
0.004.00
0.00
10.00
11.00
5.00
5.50
0.00
0.38
0.35
0.65
0.74
1.00
11.50
12.50
7.50
9.00
5.50
0.00
0.33
0.38
0.52
0.47
0.59
1.00
0.0044.50
0.00
0.44
0.54
0.62
1.00
1.000.70
1.00
0.56
0.69
1.00
0.55
0.59
0.72
1.00
E1:200—350m,E2:350—500m,E3:500—650m,E4:650—800m,E5:800—950m,E6:950—1100m
表4 林下植物群落物种多样性与地理因素及优势乔木种Pearson相关性
Table4 Correlationbetweenspeciesdiversityofunderstoryplantcommunityandgeographicalfactorsandpearsondominanttreespecies
Patrick指数
Patrickindex
Shannon⁃Wiener指数
Shannon⁃Wienerindex
灌木层
Shrub
layer
草本层
Herb
layer
Simpson指数
Simpsonindex
Pielou指数
Simpsonindexindex
草本层
Herb
layer
项目
Project
地理因素
Geographicfactors
海拔/m
经度/(°)
纬度/(°)
坡度/(°)
坡向
优势林木组成
Dominantforestcomposition
郁闭度Canopydensity
落叶松L.gmelinii/%
蒙古栎Q.mongolica/%
白桦B.platyphylla/%
黑桦B.dahurica/%
山杨P.davidiana/%
∗:P<0.05;∗∗:P<0.01
0.191
∗
-
0.161
∗
0.368
∗
-
0.238
∗
-
0.010
灌木层
Shrub
layer
∗
-
0.314
∗
草本层
Herb
layer
-
0.127
0.108
0.204
∗
0.335
∗
-
0.163
∗
-
0.155
∗
-
0.002
0.183
∗
0.171
∗
0.027
-
0.148
-
0.141
0.225
∗
∗
0.261
∗
∗
-
0.244
∗
灌木层
Shrub
layer
0.096
0.077
0.057
草本层
Herb
layer
-
0.028
-
0.060
0.080
0.035
0.163
灌木层
Shrub
layer
0.049
0.033
-
0.136
0.133
∗
-
0.428
∗
∗
0.147
∗
0.180
∗
-
0.079
-
0.088
0.096
0.024
0.070
0.070
0.180
∗
-
0.130
-
0.074
-
0.036
0.012
-
0.052
0.094
0.253
∗
∗
-
0.413
∗
0.153
∗
-
0.229
-
0.015
∗
-
0.118
0.227
∗
-
0.362
∗
-
0.103
0.130
0.039
-
0.039
-
0.094
0.245
∗
-
0.155
∗
-
0.144
-
0.058
0.072
0.115
-
0.143
-
0.136
-
0.022
-
0.026
0.055
0.034
-
0.049
-
0.035
0.042
0.055
0.065
-
0.146
-
0.158
∗
0.108
0.096
-
0.074
-
0.054
-
0.014
-
0.187
∗
-
0.158
灌木层Patrick指数与郁闭度呈极显著正相关(P<0.01),相关性为0.153,Shannon⁃Wiener指数与郁闭度
呈显著正相关(P<0.05),相关性为0.183;草本层Patrick、Shannon⁃Wiener与Simpson指数则与郁闭度呈显著
负相关(P<0.05),相关性为
-
0.413、
-
0.362与
-
0.155。说明郁闭度的升高对于灌木层与草本层α多样性存在
相反的影响,灌木层α多样性随郁闭度升高略有升高,草本层α多样性随郁闭度升高而下降明显。林木组成
http://www.ecologica.cn
7期 史恭发 等:内蒙古大兴安岭东部林下植物生物多样性海拔差异
3013
同样对林下植物多样性产生影响,林间落叶松数量占比与灌木层Patrick指数呈极显著正相关(P<0.01),相
关性为0.227,与灌木层Shannon⁃Wiener和Simpson指数呈显著正相关(P<0.05),相关性为0.171与0.180;林
间黑桦数量占比与灌木层Patrick指数呈显著负相关(P<0.05),相关性为
-
0.187;林间山杨数量占比与草本层
Pielou指数呈显著负相关(P<0.05),相关性为
-
0.158。说明调研区落叶松增多可提升灌木层生物多样性,黑
桦增多会降低灌木层物种丰富度,而山杨的增多会降低草本层物种分布的均匀度。
4 讨论
调研区为内蒙古大兴安岭的东南坡,地势沿东南至西北方迅速升高,随纬度升高逐渐升高,为寒温带气候
区内典型的山岭地区,拥有广袤的森林地带。林下植物是森林生态系统的重要组成部分,研究寒温带山岭地
区林下植物生物多样性随海拔差异变化有助于了解寒温带山岭地区植物随海拔梯度分布规律,丰富生物多样
性海拔格局理论,掌握内蒙古大兴安岭地区林下植物资源现状。调研区林下植物共277种隶属于53科135
属,灌木32种、草本245种,与张喜亭等人
[42]
对大兴安岭多布库尔国家级自然保护区调研的269种植物相
似。调研区森林乔灌群落随海拔梯度不断变化,乔木层由低海拔梯度以蒙古栎、黑桦为主,逐渐变为以落叶
松、白桦为主,灌木层由以榛、胡枝子为主过渡到兴安杜鹃、越橘等为主,与邻近地区类似海拔梯度如鄂伦春旗
多布库尔保护区(400—600m)
[42]
、呼中森林自然保护区(800—1200m)
[43]
乔灌群落构成具有较高的相似性。
林下植物调查中发现了野大豆(Glycinesoja)、大花杓兰(Cypripediummacranthos)、蒙古黄芪(Astragalus
危、重点保护植物,使得了解与保护内蒙古大兴安岭东部林下植物生物多样性更有必要。
林下植物α多样性随海拔梯度呈现先上升后下降趋势,与常见的生物多样性随海拔升高的单峰格局相
符
[17,36,44]
。林下草本层为林下植物生物多样性的主要贡献层,其α多样性高于灌木层,这与岳永杰
[25]
、张殷
波等人
[45]
研究一致。林下草本层与灌木层α多样性的海拔差异性存在明显不同,随海拔上升,草本层α多样
性均先下降,在E2区间内达到谷值,至E3梯度到达最高峰后持续下降的变化趋势,与传统的单峰格局有所
不同,根据实际调研情况,E2梯度多位于莫旗与鄂旗的东南部,此地区人口相对较多,森林内部常出现大量农
田区域,一方面严重破坏了森林连续度,另一方面农田耕作造成明显的水土流失、农药污染,加之此地区的放
牧、采山等活动,推测其为E2梯度林下草本群落α多样性下降的重要原因之一。灌木层Patrick、Shannon⁃
Wiener与Simpson指数随海拔梯度升高呈逐渐上升趋势,而Pielou指数呈波动下降趋势,灌木层相比于草本
层更能适应高海拔梯度的生境。
灌木层Cody指数在E2与E3、E3与E4相邻两海拔梯度值逐渐升高,在E3与E4海拔区域内灌木物种替
换速率最高,灌木也在此区间内由植株较高,叶片较大的榛、胡枝子、珍珠梅等,转变为了植株较矮,叶片较小
的土庄绣线菊、绣线菊、刺蔷薇等植物,随着海拔进一度升高,越橘、杜香、笃斯越橘、偃松等植叶片更小的灌木
逐渐变为优势物种。草本层则具有更高的Cody指数,在E2与E3区间内物种替换速率最高。
大尺度区域内,植物群落多样性主要受到气候、地势的影响,小尺度区域易受到区域小气候、土壤等方面
的影响
[14]
。相比于单独山体的植物多样性海拔差异性研究
[15—18]
,山岭地区植物多样性不仅会受到小区域山
体坡向、坡度的影响,同时受到经度与纬度与海拔相结合的多重影响。林下灌木层多样性与海拔、纬度呈现正
相关、与经度呈现负相关,根据相关性,这三种因素对于灌木层生物多样性影响排序为纬度>海拔>经度,调研
区随着纬度的上升海拔逐渐升高。草本层与海拔呈负相关,与纬度呈正相关,结合调研区地势随经纬度变化
的情况,这种相关性符合前面的结果,值得注意的是,纬度同时与灌木层、草本层均呈正相关,推测一方面这是
由于调研区东南部低纬度地区人口密集,农耕放牧活动对森林生态造成的破坏导致林下植物生物多样性的下
降,与另一方面调研区西北部纬度较高地区相较于东南部具有更复杂的地形起伏更高的环境异质性,并且良
好的封山育林环境促使林下生物多样性上升导致的。从行政区域表现则为莫旗林下植物多样性低于鄂旗。
坡度方面,平缓地区灌木多样性较高。坡向方面,草本层多样性阴坡小于阳坡,而对灌木层多样性影响不明
http://www.ecologica.cn
membranaceusvar.mongholicus)、紫点杓兰(Cypripediumguttatum)、樱草(Primulasieboldii)等多种国家珍稀濒
3014
生 态 学 报 44卷
显。郁闭度的增加一定程度上能提高灌木层多样性,却显著降低草本层多样性,鉴于草本层对于林下植物多
样性的突出贡献,建议定期对过密林分间伐以获取较好的林下植物多样性,间伐强度与频率有待于进一步研
究
[46]
。此外,不同的乔木组成能够提供不同的林下生存环境,从而产生不同的林下植物群落结构,调研区内
落叶松的加入能够一定程度上增加灌木层多样性,并有降低草本层多样性的趋势,这可能是相对于阔叶林,落
叶松林下凋落物分解速率较慢,养分补充缓慢,表层土壤组成不同导致的。
本研究主要针对海拔这一因素,探讨内蒙古大兴安岭地区林下植物多样性海拔差异性变化。林下植物多
样性受到多种环境异质性影响,调研时还发现向阳山坡通常具有较多种类的林下植物,林缘林下植物多样性
普遍较高,较潮湿或离水源较近的林地物种乔木与灌木丰富度往往较高等规律,不同的土壤、不同的林型
[47]
、
不同坡向与坡度
[48]
,不同季节、不同林龄、不同育林手段对林下植物多样性均可能产生影响。在今后的研究
应重视不同因素综合对林下植物的影响,并增加不同调研季节,考虑气候因素,全面分析,对于掌握林下植物
多样性的多方面规律具有重要意义。
5 结论
内蒙古大兴安岭东部林下植物共277种隶属于53科135属,灌木32种、草本245种;随海拔梯度升高,
森林乔灌群落依次为蒙古栎⁃黑桦⁃榛⁃胡枝子群落(E1—E2)、落叶松⁃白桦⁃榛⁃胡枝子群落(E3)、落叶松⁃白桦⁃
欧亚绣线菊⁃绣线菊群落(E4)、落叶松⁃白桦⁃兴安杜鹃⁃越橘群落(E5)、落叶松⁃白桦⁃越橘⁃偃松群落(E6);林
下植物物种多样性随海拔上升呈现明显的单峰格局。草本层物种多样性普遍高于灌木层,灌木层生物多样性
随海拔升高呈上升趋势,草本层呈现先下降后升高再下降的波动下降趋势;E3—E4与E2—E3梯度内灌木层
与草本层物种替换速率分别达到峰值,草本层物种替换速率远高于灌木层;灌木层与草本层相邻海拔梯度间
相似度指数均呈先下降后升高再下降的波动形式;地理因素及优势乔木种不同能够显著影响林下植物α多样
性,灌木层与海拔、纬度呈正相关,与经度、坡度呈负相关,草本层与海拔呈负相关,与纬度、坡向呈正相关,海拔
对草本层α多样性影响高于灌木层。郁闭度与灌木层α多样性呈正相关,与草本层α多样性呈负相关;落叶松
增多可提升灌木层生物多样性,黑桦增多会降低灌木层物种丰富度,而山杨的增多会降低草本层物种分布的均
匀度。
参考文献(References):
[1] 牟长城,倪志英,李东,孙晓新,陈加利.长白山溪流河岸带森林木本植物多样性沿海拔梯度分布规律.应用生态学报,2007,18(5):
[2] MåsvikenJ,DalerumF,CousinsSAO.ContrastingaltitudinalvariationofalpineplantcommunitiesalongtheSwedishMountains.Ecologyand
[3] OldfatherMF,BrittonMN,PapperPD,KoontzMJ,HalburMM,DodgeC,FlintAL,FlintLE,AckerlyDD.Effectsoftopoclimatic
[4] 刘华训.我国山地植被的垂直分布规律.地理学报,1981,36(3):267⁃279.
7170⁃7180.
complexityonthecompositionofwoodyplantcommunities.AoBPlants,2016,8:plw049.
[5] 孙然好,陈利顶,张百平,傅伯杰.山地景观垂直分异研究进展.应用生态学报,2009,20(7):1617⁃1624.
Evolution,2020,10(11):4838⁃4853.
945⁃952.
[6] 杨朗生,刘兴良,刘世荣,潘红丽,冯秋红,李迈和.卧龙巴郎山川滇高山栎群落植物生活型海拔梯度特征.生态学报,2017,37(21):
[7] 金慧,赵莹,赵伟,尹航,秦立武,崔凯峰,李晓杰,黄祥童.长白山牛皮杜鹃群落物种多样性的海拔梯度变化及相似性.生态学报,
[8] JeníkJ.H.WalterVegetationoftheearthinrelationtoclimateandtheeco⁃physiologicalconditions.FoliaGeobotanicaetPhytotaxonomica,1975,
[9] KörnerC.Whyarethereglobalgradientsinspeciesrichness?mountainsmightholdtheanswer.TrendsinEcology&Evolution,2000,15(12):
[10] UlfMolau,方静.低纬与高纬地区的山地生物多样性格局.AMBIO⁃人类环境杂志,2004,33(S1):25⁃29,57.
1⁃7.
513⁃514.
[11] 姚兰,艾训儒,吕世安,冯广,刘峻城,黄永涛.湖北星斗山天然次生林的群落类型、结构与物种多样性特征.林业科学,2015,51(11):
[12] 朱源,康慕谊,江源,刘全儒.贺兰山木本植物群落物种多样性的海拔格局.植物生态学报,2008,32(3):574⁃581.
[13] 梁红柱,刘丽丽,付同刚,高会,李敏,刘金铜.山地植被垂直分布变化格局研究进展与述评.中国生态农业学报:中英文,2022,30
10(1):100.
2015,35(1):125⁃133.
http://www.ecologica.cn
7期
(7):1077⁃1090.
史恭发 等:内蒙古大兴安岭东部林下植物生物多样性海拔差异
3015
[14] 刘星宏,张青青,徐海量,张鹏,张广鹏,李宏.北疆盐碱地植物群落空间分布及物种多样性.生态学报,2021,41(4):1501⁃1513.
[16] 冯相艳,赵文智,蔺鹏飞,王川.祁连山北坡主要木本植物功能性状及其海拔分异.生态学报,2022,42(23):9726⁃9735.
报,2021,41(3):1148⁃1157.
[15] 陈绪辉,林碧华,罗敏贤,刘宝,刘金福,郑世群.福建灵源山森林群落林下灌草生态位特征.森林与环境学报,2023,43(2):152⁃159.
[17] 李梦佳,何中声,江蓝,谷新光,晋梦然,陈博,刘金福.戴云山物种多样性与系统发育多样性海拔梯度分布格局及驱动因子.生态学
[18] 向琳,陈芳清,官守鹏,王玉兵,吕坤.井冈山鹿角杜鹃群落灌木层功能多样性及其随海拔梯度的变化.生态学报,2019,39(21):
[19] 连承杰.内蒙古大兴安岭林区林下经济发展现状与对策.中国林业经济,2021(1):43⁃45.
76⁃83.
8144⁃8155.
[20] 王群,李东海,杨小波,曾润娟,夏丹,王豪,戚春林,郝杰威.海南霸王岭橡胶林林下植物多样性特征研究.林业资源管理,2022(5):
[21] 谭许脉,张文,肖纳,莫雪青,高冠女,尤业明,明安刚,黄雪蔓.杉木林改造成乡土阔叶林对林下植物物种组成和多样性的影响.生态
[22] 王纯华.大兴安岭林区乔木树种与林下伴生经济植物分布情况.现代农业科技,2020(13):120,122.
学报,2022,42(7):2931⁃2942.
[23] 王纪坤,兰国玉,陈帮乾,谢贵水.海南橡胶林林下药用植物资源调查.南方农业学报,2013,44(3):391⁃396.
[24] 陈树彪,郑野,李秀英.大兴安岭呼中林业局林下主要经济植物资源的调查.内蒙古林业调查设计,2011,34(5):27⁃28,124.
[26] 张建宇.大兴安岭森林植物多样性、群落结构特征及耦合关系分析[D].哈尔滨:东北林业大学,2018.
[29] 冶连华,吕金娃.内蒙古大兴安岭林区野生宿根花卉资源调查.河北林业科技,2013(2):25⁃29.
37,134.
版,2020,41(4):21⁃26.
[25] 岳永杰,苏志成,韩君,格日乐高娃,乌云珠拉.内蒙古大兴安岭林区典型森林群落物种多样性研究.内蒙古农业大学学报:自然科学
[27] 崔佳佳,铁牛.大兴安岭北部森林群落结构及植物多样性特征研究.西北林学院学报,2021,36(2):24⁃30.
[28] 李桂茹.内蒙古大兴安岭林区野生经济植物资源评价.内蒙古林业调查设计,2011,34(1):98⁃99,120.
[30] 王翠敏.内蒙古大兴安岭林区兴安落叶松林下分布主要经济植物的经济价值及开发与利用.内蒙古林业调查设计,2014,37(2):35⁃
[31] 魏麒苏.内蒙古大兴安岭北部原始林区植物资源评价.内蒙古林业调查设计,2020,43(2):32,63⁃64.
41⁃44.
[32] 张重岭,周本志,兰学君,杜德刚.内蒙古大兴安岭北部原始林区兴安落叶松林植物多样性的调查.内蒙古林业调查设计,2003,26(1):
[33] 王丽红,辛颖,邹梦玲,赵雨森,刘双江,任清胜.大兴安岭火烧迹地植被恢复中植物多样性与生物量分配格局.北京林业大学学报,
[34] 王纯华.大兴安岭野生经济植物资源调查方法.现代农业科技,2020(10):122,125.
30⁃35.
2015,37(12):41⁃47.
[35] 丛静,尹华群,卢慧,宿秀江,杨敬元,李迪强,张于光.神农架保护区典型植被的物种多样性和环境解释.林业科学,2013,49(05):
[36] 董雪,李永华,辛智鸣,段瑞兵,姚斌,包岩峰,张正国,刘源.河西走廊西段荒漠戈壁灌木群落物种多样性的海拔格局.林业科学,
[37] 傅沛云.东北植物检索表.第2版.北京:科学出版社,1995:1⁃906.
自然保护区生物多样性研究.北京:科学出版社,2015:1⁃150.
2021,57(2):168⁃178.
[38] 赵敏,郑宝江,孙海,吕连宽,赵宝军,胡建民,周明,龚文峰,许青,崔岱宗,吴玉环,任强,陶雷,邵景文,黄世臣.额尔古纳国家级
[39] 丁茂,郑标,王丽,张金铭,周守标,王智.鹞落坪国家级自然保护区草本植物多样性及空间分布.生态学报,2022,42(19):8015⁃8030.
[40] 黄甫昭,丁涛,李先琨,郭屹立,王斌,向悟生,文淑均,李冬兴,何运林.弄岗喀斯特季节性雨林不同群丛物种多样性随海拔的变化.
[41] 张涵丹,康希睿,邵文豪,杨旭,张建锋,刘学全,陈光才.不同类型杉木人工林林下草本植物多样性特征.生态学报,2021,41(6):
[42] 张喜亭,张建宇,肖路,陈胜仙,仲召亮,高薇,王文杰.大兴安岭多布库尔国家级自然保护区植物多样性和群落结构特征.生态学报,
[43] 孙玉成,张喜亭,季倩如,张建宇,王文杰.大兴安岭泰加林植物多样性特征调查与分析.林业资源管理,2022(6):124⁃130.
化植物群落多样性和地上生物量沿海拔梯度的变化特征.生态学报,2022,42(9):3640⁃3655.
2022,42(1):176⁃185.
[44] 王采娥,黄梅,王文银,李子好,张涛,马林,白彦福,王彦龙,施建军,龙瑞军,刘玉,王晓丽,马玉寿,尚占环.三江源区高寒坡地退
[45] 张殷波,秦浩,孟庆欣,张峰,唐志尧.太行山森林群落物种多样性空间格局及其影响因素.应用与环境生物学报,2022,28(2):
[46] 龚固堂,牛牧,慕长龙,陈俊华,黎燕琼,朱志芳,郑绍伟.间伐强度对柏木人工林生长及林下植物的影响.林业科学,2015,51(4):
[47] 杨光辉,秦树林,金光泽.小兴安岭五种林型早春草本植物多样性及其环境解释.生态学报,2023,43(3):1234⁃1246.
8⁃15.
[48] 崔博亮,安彦明,张起鹏,蒲萌,陈克龙.坡度对亚高寒草甸狼毒群落植物多样性的影响.草业科学,2021,38(5):823⁃834.
331⁃338.
2118⁃2128.
生态学报,2016,36(14):4509⁃4517.
http://www.ecologica.cn