曾文豪,石 慰,唐一思,郑维艳,曹坤芳,*
1 广西大学林学院,广西森林生态与保育重点实验室,南宁 530004 2 广西大学, 亚热带农业生物资源保护与利用国家重点实验室,南宁 530004
生物多样性的形成和维持机制一直是生态学研究的核心问题之一[1- 3]。很长时间里,生态学家基于环境因素对群落多样性分布格局的形成进行了大量的探索,取得许多重要的发现。生物的代谢和生长离不开温度,物种的分布范围受环境最低温的限制[4];一般情况下,物种多样性随温度的升高而升高[5]。生态学代谢理论认为,较高的环境温度能够促进物种的新陈代谢,进而提高地区的物种丰富度[6]。土壤是植物生长的基础,土壤的质量与物种丰富度的关系密切相关[7- 8]。任学敏等在秦岭太白山巴山冷杉(Abiesfargesii)-糙皮桦(Betulautilis)混交林的调查中发现物种丰富度随土壤pH值升高而增加,随岩石盖度升高而减小[9]。Dornbush等研究发现高草草原土壤厚度的增加有利于提高物种丰富度[10]。海拔作为环境因子的一个集合体,它综合反映着环境温度、降雨、光照和土壤等梯度的变化,这些因子影响群落物种组成以及系统发育结构。卢孟孟等对云南哀牢山亚热带森林的研究中发现群落系统发育结构随着海拔梯度的升高由聚集逐渐变为发散或聚集程度降低[11]。Culmsee在研究马来西亚亚热带森林时发现,从低海拔到高海拔系统发育结构由发散逐渐转向聚集[12]。虽然针对不同区域的研究呈现的结果有所差异,但是这些研究足以说明环境因子对物种的形成及群落的系统发育结构产生影响。目前,群落系统发育学被广泛用于群落的构建机制和进化历史研究,在一个群落的演化过程中,如果群落的构建是以环境选择为主要机制,其系统发育结构表现为聚集;如果是以物种竞争为主要构建机制,群落系统发育结构表现为离散;如果是以随机过程为主要构建机制,则群落的系统发育结构呈随机状态[13- 14]。
我国的亚热带地区受海洋季风影响,光热水资源丰富、土壤肥沃,地带性植被为常绿阔叶林,其生态系统生产力高[15]。广西处于低纬度地区,自北向南占据着亚热带(包括中亚热带和南亚热带)及逐渐向热带过渡的区域,依次分布着中亚热带常绿阔叶林带、南亚热带常绿季雨林带、北热带季节性雨林带三个植被带;由东向西降水减少,广西西部地区为湿润的气候区向半湿润气候区的过渡区[16]。
喀斯特地貌以碳酸盐为主的岩溶构造,多表现为峰林、峰丛、洼地、峡谷等地貌特征;具有土层浅薄偏碱性、富钙镁,岩石裸露,土被不连续,土壤蓄水能力低等特点[17- 18]。复杂独特的喀斯特地质气候条件构成了许多的适合不同植物生存的小生境,决定了其植物群落的特殊性[19- 20]。通过对喀斯特森林和非喀斯特森林树木物种多样性的研究,揭示两种生境森林群落的形成机制的差别,对森林资源的保护和利用提供有价值参考[21]。独特的喀斯特生境几十年来吸引了众多学者对其森林的研究,但大多数都是针对个别保护区开展的研究,比较少的从较大区域尺度上开展对喀斯特和非喀斯特森林群落构建机制、系统发育多样性之间差异的探讨[22- 24]。本研究以广西地区喀斯特与非喀斯特山地森林为研究对象,结合系统发育学的在生态学中的应用探讨下列两个科学问题:(1)海拔≤1000 m喀斯特和非喀斯特森林以及海拔>1000 m的非喀斯特森林之间乔木植物科组成的差异。 (2)这三类生境的森林物种丰富度和系统发育多样性的差异及其与生长地点年平均温度的关系,讨论不同生境山地森林群落的主要构建机制。
广西地处华南,104°26′—112°04′E,20°54′—26°24′N,陆地总面积23.7万km2,其中喀斯特分布面积占广西总面积35%,主要分布于桂西南、桂西北、桂中和桂东北(图1)。广西地区水热同期且丰沛,各地年日照时数1169—2219 h,各地年均气温16.5—23.1℃,累年极端最高温33.7—42.2℃,累年极端最低气温为-8.4—2.9℃,日均温≥10℃的积温在4995—8328℃之间。常年降水量1086.3—2754.5 mm,季节分配不均[16]。
本研究的山地森林群落数据主要来源于《广西植被志要》[25],将收集到的样地分为三类生境,一是海拔≤1000 m非喀斯特样地,二是海拔>1000 m非喀斯特样地,三是海拔≤1000 m喀斯特样地。各个生境样地信息如表1,样地地理分布如图1。从国家气象网收集到的广西地区90个气象站30年气象数据,建立年平均温度与经度(G)、纬度(L)和海拔高度(H)的回归模型(年平均温=72.794-0.736×L-0.31×G-0.005×H,R2=0.938,P<0.001),以此为基础,结合每个样地点的经度、纬度和海拔计算出每个样地点的年平均温度[26]。
表1 样地数据收集概况
图1 样地分布图(部分样地重叠)Fig.1 The distribution map of community sample plots (Partially overlapping)
2.2.1 群落物种多样性和系统发育多样性分析
考虑到部分样方面积的不一致性,采用Gleason指数来表示物种多样性的测度[27],比较和分析广西地区喀斯特和非喀斯特森林树木物种多样性的差异。采用Faith系统发育多样性指数PD(phylogenetic diversity)来分析群落系统发育多样性[28],使用协方差分析的方法检验不同生境群落Gleason物种多样性指数与系统发育多样性线性回归线的差异。Gleason指数公式如下:
式中,A表示所研究的面积(m2);S表示面积A内的种数(个)。
2.2.2 系统发育树的构建和系统发育指数的选择
将收集到的样地树木的拉丁名名录输入到Phylomatic平台[29](http://phylodiversity.net/phylomatic/),利用Phylomatic中数据库构建物种系统发育树,该平台以被子植物分类系统APG Ⅲ 为基础,整合Zanne的进化树信息,输出带枝长的系统进化树[30]。净亲缘关系指数(net relatedness index, NRI)和最近分类群指数(nearest taxon index, NTI)被广泛用于研究群落的系统发育结构[13,31]。当NRI、NTI>0,说明群落物种组成在系统发育结构上聚集,群落由亲缘关系比较近的物种组成,生境过滤是群落构建的主要因素;当NRI、NTI<0,说明群落系统发育结构发散,群落由亲缘关系较远的种构成,竞争排斥是群落构建的主要控制因素;当NRI、NTI=0,说明群落的系统发育结构呈随机状态,环境因素和生物因素共同维持群落生物多样性[13]。使用t-检验分析不同生境间指数是否存在显著性差异。NRI和NTI计算公式如下:
式中,MPD表示物种间的平均成对系统发育距离,MPDsample表示实际观察到的MPD值,MPDrandsample代表随机产生的999个随机组合的群落的MPD值,SD(MPDrandsample)为这些随机值的标准差。MNTD表示平均最近种间系统发育距离,MNTDsample表示群落内亲缘关系最近的物种之间的平均系统发育距离,MNTDrandsample表示随机产生的999个随机组合的群落的MNTD值,SD(MNTDrandsample)为这些随机值的标准差。
将样地群落按照前面所述3种不同生境统计群落的科组成及其物种数情况。发现所研究群落样地共有树木936个种,按照中国植物志分类系统,他们隶属于98科361个属,其中从科的组成上看广西地区海拔≤1000 m的非喀斯特群落中共有87个科,占所有研究群落植物科的88.8%,其中含有4种以上的有樟科、壳斗科、茜草科、桑科、山矾科等44科,共计709个种,占所研究群落所有物种种数的75.7%。海拔>1000 m非喀斯特森林共有61科,占所研究群落所有树木科的62.2%,其中含4种以上的有樟科、壳斗科、山茶科、杜鹃花科等20科,共计360个种,占所研究群落所有物种种数的38.5%。喀斯特森林共有70科,占所研究群落所有植物科的71.4%,其中含有4种以上的有大戟科、樟科、豆科、桑科、茜草科等26科,共计308个种,占所研究群落所有物种种数的32.9%(图2)。
图2 广西海拔≤1000 m非喀斯特、海拔>1000 m非喀斯特和海拔≤1000 m喀斯特天然林科与物种数组成及其所占比例Fig.2 Size of families (%) containing over 4 species of 3 habitats: non-karst with elevation of ≤1000 m, non-karst with elevation of >1000 m and karst with elevation of ≤1000 m科前面数字代表他们所包含的物种数
3种生境比较来看,海拔≤1000 m非喀斯特森林含有4种以上的植物科比海拔>1000 m的非喀斯特中山森林多了24科,其中樟科、壳斗科、山茶科、杜鹃花科等都是两种生境的重要组成科。大戟科为喀斯特森林含有种数最多的科,其中出现在非喀斯特森林中的山茶科和杜鹃花科等没有在喀斯特森林出现。
海拔>1000 m非喀斯特中山生境的Gleason指数的平均值显著高于另外两种生境的森林,海拔≤1000 m非喀斯特和喀斯特两种生境的Gleason指数没有显著差异。其中,喀斯特森林样地分布范围比海拔≤1000 m非喀斯特森林要窄,前者主要分布于年平均温18—22℃的生境之间,后者主要分布于年平均温15—23℃的生境之间,海拔大于1000 m非喀斯特森林样地分布范围的年平均温比低海拔两种生境都要低,且其Gleason指数与年平均温表现为一定的反相关(图3)。
图3 三种生境森林样地Gleason多样性指数与年平均温散点图Fig.3 Lack of correlation between Gleason diversity index and mean annual temperature in the forest communities in all of the three habitats黒点R2=0.2%,P=0.74;红点R2=2.3%,P=0.02; 绿点R2=20.1%,P=0.0003。小图展示三种生境的平均多样指数和标准差,不同字母表示差异显著(P<0.01)
图4 三种生境森林样地系统发育多样性PD指数与Gleason相关散点图Fig.4 The correlation between Gleason diversity index and phylogenetic diversity index in the forest communities in all of the three habitats 黒点R2=90%,P<0.0001;红点R2=95%, P<0.0001;绿点R2=94%, P<0.0001
3种生境的森林群落系统发育多样性PD指数与Gleason指数都表现为显著的线性正相关关系。协方差分析结果表明PD指数与Gleason指数的关系的回归线斜率在3种生境森林之间的没有表现出显著差异,截距表现出显著差异(P<0.05)。在相同的Gleason指数情况下,海拔≤1000 m喀斯特森林样地系统发育多样性比海拔≤1000 m非喀斯特森林样地离散程度小,海拔≤1000 m非喀斯特森林样地系统发育多样性PD指数比>1000 m的要高(图4)。
从净亲缘关系指数NRI来看,除了小部分样地表现为系统发育发散外,3种生境其他大部分样地均表现为系统发育聚集。喀斯特森林和中山地带的非喀斯特森林绝大部分的样地表现为系统发育聚集,而低海拔非喀斯特森林有约25.7%的样地表现为系统发育离散(图5)。
图5 净亲缘关系指数(NRI)、最近分类群指数(NTI)与年平均温度和海拔的变化Fig.5 Variation in NRI and NTI with (A and B) mean annual temperature (MAT) and elevation A图:黒点P=0.47;红点R2=5%, P=0.0007;绿点P=0.99; B图:黒点P=0.08;红点R2=5.7%, P=0.0003;绿点P=0.26
从最近分类群指数NTI来看,3种生境的大部分样地都表现为群落系统发育聚集(图5)。但是与系统发育结构NRI不同的是两种非喀斯特生境的森林群落有部分样地的系统发育结构(NTI)聚集程度比喀斯特森林的要大。
海拔>1000 m非喀斯特森林的Gleason指数总体上都要比海拔≤1000 m非喀斯特和喀斯特两种森林的Gleason指数稍高,同时海拔≤1000 m和>1000 m两种非喀斯特生境的森林含有4种以上的植物科组成上有所不同,并且前者含有4种以上的科比后者要多,造成这差异的原因可能有: 第一,在一定海拔范围内,海拔>1000 m的环境因子并不能对物种的分布造成限制影响,何艳华等在对山西省五鹿山森林物种多样性海拔格局研究中发现在一定范围内物种丰富度随着海拔的升高而逐渐增加,在海拔1500 m左右达到最大[32]。卢孟孟等在对云南哀牢山亚热带森林的调查同样呈现在一定海拔范围内物种丰富度随着海拔的升高而逐渐增加的结果,并在海拔1700 m左右达到最大[11]。朱彪等在对广西猫儿山的森林调查中发现类似的结果,物种丰富度在海拔1350 m以下变化不明显,但在1350 m以上随着海拔的升高物种丰富度呈下降趋势[33]。在我们的研究中,海拔大于1000 m生境的样地中,海拔范围在1000 m至1500 m的样地占了72.1%,海拔范围在1000 m至1700 m的样地占了91.8%。第二, 由于较低海拔范围内,受人为干扰较大,降水较少,土壤元素含量较低,可能会导致物种丰富度较低。随着海拔的升高,降水增加,气温虽然降低但是不太寒冷,土壤中养分含量有所增加,导致适合生存的物种增多[32]。海拔>1000 m非喀斯特森林样地Gleason指数与年均温呈反相关的关系,这与前人对物种丰富度在一定范围内随温度的升高而升高的基本认识相勃[34],造成这一结果的原因可能与中山地区一些特殊的小生境有关。海拔≤1000 m非喀斯特森林和喀斯特森林在科属组成有所不同,前者种数最多的科为樟科,后者为大戟科,其中出现在非喀斯特森林的山茶科和杜鹃花科等都没有在喀斯特森林中出现。造成这种差异的原因可能是由于喀斯特特殊的地质和地形所致,喀斯特基质漏水,土层薄,岩石裸露,容易缺水,并且土壤富积钙,中性或弱碱性,使得嫌钙植物不能生长,而一些喜酸土、耐酸土植物也不能在喀斯特生境分布。
协方差分析表明三种生境森林群落PD指数与Gleason指数的回归线的截距显著差异。总体上看来Gleason指数相同的情况下海拔≤1000 m喀斯特森林样地系统发育多样性PD指数比海拔≤1000 m非喀斯特生境的森林样地离散程度小, 海拔≤1000 m非喀斯特森林样地系统发育多样性PD指数比>1000 m的要高,推测海拔≤1000 m非喀斯特森林生境比海拔>1000 m的非喀斯特生境更适合不同亲缘关系的树种生存。
从净亲缘关系指数(NRI)和最近分类群指数(NTI)来看,3种生境的大部分样地都倾向于系统发育聚集。这说明生境过滤对于这3种生境的群落构建起主要作用,使得发育关系相近物种得以生存,从而表现出系统发育聚集[35]。但是,无论从NRI指数还是NTI指数来看,海拔≤1000 m喀斯特(NRI 3.8%,NTI 15%),海拔>1000 m非喀斯特(NRI 11.5%,NTI 8.2%)两种生境森林群落发生系统发育发散的比例都比海拔≤1000 m非喀斯特(NRI 25.7%,NTI 20.4%)低,揭示生境过滤在这两个生境中表现的更为明显。从净亲缘关系指数(NRI)来看,海拔≤1000 m非喀斯特森林样地呈现系统发育聚集比例(74.3%)比海拔>1000 m非喀斯特(88.5%)生境要小,这与前人的研究结果类似,之前的研究表明,系统发育聚集程度随着海拔的升高而逐渐加强或表现为单峰模式[36];而喀斯特森林的群落系统发育结构聚集程度要比其他两种生境强(96.3%),这可能与喀斯特生境选择压力大相关。从最近分类群指数(NTI)来看,得出的结果与净亲缘关系指数(NRI)有所差别,其中非喀斯特海拔>1000 m的系统发育聚集程度(91.8%)要比另外两种生境要强。在前人的研究中也有类似的结果,随着海拔的升高,温度、水分和土壤等都发生变化,系统发育结构逐渐聚集[37]。
本研究通过分析广西地区喀斯特和非喀斯特山地森林群落的物种多样性和系统发育结构发现,海拔>1000 m非喀斯特生境物种丰富度比另外两种生境要高;海拔≤1000 m的喀斯特和非喀斯特两种生境物种丰富度没有显著差异,说明喀斯特的生境虽比非喀斯特生境严峻,但喀斯特形成许多水热条件差异的小生境使不同类型的物种得以生长,从而保持喀斯特生境也具有较高的物种丰富度。同时,在3种生境中大部分样地的森林群落都表现出系统发育聚集的现象,较少样地表现出发散或随机的现象,说明生境过滤是广西地区天然林群落构建的主要驱动因素。研究结果揭示中山地带的森林树木物种多样性与低海拔森林相似,喀斯特森林具有许多特有植物,表明它们与低海拔非喀斯特生境的森林一样具有高的生物资源保护价值。